![7d867a57481a3e5b742c950a9b60e28d.png](https://i-blog.csdnimg.cn/blog_migrate/596dc8271946f7158e603586a3cc14ac.png)
当有限元分析中需要使用超弹性材料模型时,工程师通常很少有实质性的数据来帮助他们进行超弹有关的非线性分析。有时幸运的工程师会有一些拉伸或压缩应力-应变实验数据,或者还会有一个简单的剪应力-应变实验数据。正确地处理这些数据是分析非线性超弹模型的关键一步,而对这些数据进行曲线拟合而得出对应超弹模型的材料参数就显得尤为重要。今天我们就来了解一下超弹模型相关的实验数据及其曲线拟合,最后也会在MatEditor软件中给出曲线拟合的实例。
![967a11ef9a9f336c738ca96d20f8dd46.png](https://i-blog.csdnimg.cn/blog_migrate/be34633fdeacbc881c0ae8ab70c3e734.jpeg)
超弹材料的力学实验数据
应变能密度函数中的材料常数确定了超弹性模型的力学响应。为了在超弹性分析中获得正确的结果,需要评估被测材料的模型常数。这些常数通常是根据实验应变-应力数据通过曲线拟合得出的。测试数据通常取较大范围内的几种变形模式。使用至少与有限元分析中所要经历的变形状态一样多的测试数据,拟合出的材料常数才能准确体现材料的力学响应。
对于超弹性材料,可以使用简单的变形测试来拟合材料常数。下图表现了6个不同的变形模式。来自多种变形状态下的数据组合要比单一类型的变形数据更好,更可以得到准确的超弹材料参数。
![f4ba704f1a8f9e4c4db9582d79876092.png](https://i-blog.csdnimg.cn/blog_migrate/41520ebe9e3fcfbd6022d41fdaf4ddc3.jpeg)
尽管存在6种不同的变形状态,但是我们发现在施加静水压力后,以下变形模式会变得相同:单轴拉伸和等双轴压缩,单轴压缩和等双轴拉伸,平面张力和平面压缩。通过这些等效的测试模式,我们只需要这三种独立的变形模式的实验数据。
单轴拉伸
一种简单拉伸测试的应力-应变数据。必须提供这类测试数据,才能模拟计算拉伸应变状态。测试数据中不存在侧向约束。单轴压缩数据可以从等双轴数据中得出。
单轴拉伸下的第一和第二应变不变量为:
![ab28b7886e02b7f6892b673201811880.png](https://i-blog.csdnimg.cn/blog_migrate/b28789943d37ebe342ce521fba799bee.png)
所对应的工程应力与主拉伸率的关系为:
![7b84ee0437d53f20c0ddf89ea0ab103a.png](https://i-blog.csdnimg.cn/blog_migrate/ab1c608565e41fe46bfba746a45fceac.png)
等双轴拉伸
等双轴拉伸测试中的应力-应变数据。通过在两个方向上均匀拉伸,体现出相当于纯压缩的应变状态。单轴压缩数据可以从等双轴数据中得出。
等双轴拉伸下的第一和第二应变不变量为:
![8e5331f32d21f7e32d9dc2630278a2c9.png](https://i-blog.csdnimg.cn/blog_migrate/22233c72f5386cea43f2846b44d9cba3.png)
所对应的工程应力与主拉伸率的关系为:
![a6723e56240c3eae3362595fae312a4b.png](https://i-blog.csdnimg.cn/blog_migrate/ed40f2575fb09d595bbdb2bdc99f4141.png)
纯剪切变形
对于超弹性材料,这通常是指大幅拉伸的样品,但由于该材料几乎不可压缩,因此存在纯剪切状态。
第一和第二应变不变量为:
![07269cc2dd8a79bc46cd253bcf147b50.png](https://i-blog.csdnimg.cn/blog_migrate/be872a5f5fc9acbfdae0bce82140093d.png)
所对应的工程应力与主拉伸率的关系为:
![90977cd4cd75d3256d927d54f670cbb0.png](https://i-blog.csdnimg.cn/blog_migrate/afe5533ec21428171b748981841d12d3.png)
体积变形
此数据用于确定体积模量。对于超弹材料,如果材料可稍加压缩或部件整体受约束,则此数据较为重要。体积模量通常比剪切模量大2-3个数值阶次。对于泡沫材料,体积数据对计算材料的可压缩性十分重要。
超弹模型与曲线拟合
由于每一种超弹模型对应的材料常数数量和物理意义都不一样,需要在这里给出每种超弹模型的应变能。目前,MatEditor已经支持下列常见超弹模型:
Arruda-Boyce模型
![5f488fe80d4b1b834e97050579984424.png](https://i-blog.csdnimg.cn/blog_migrate/dd2f19fb28d31f565582671e35556c3d.png)
其中u是初始剪切模量,lambda是有限网格拉伸,D1是不可压缩参数。
Blatz-Ko模型
![d1681418d3620c9b6ad9eea98421524f.png](https://i-blog.csdnimg.cn/blog_migrate/404740fc335f2c1392d99985cd9205a4.png)
其中u是初始剪切模量。
Gent模型
![e706c726c9529d85fe4c5d7fefecdc82.png](https://i-blog.csdnimg.cn/blog_migrate/7ac66ac9dd3e3957e2f935381d3e1743.png)
其中u是初始剪切模量,Jm是I1-3项的限定值,D1是材料不可压缩参数。
Mooney-Rivlin 2参数模型
![be8d56eee11b3e52719e026e8729364f.png](https://i-blog.csdnimg.cn/blog_migrate/9559467b07fb5718aac0dff3caf721bd.png)
其中,C10,C01,和D1是材料参数。
Mooney-Rivlin 3参数模型
![0d4a0e8e0fa2695d3f0e6747d06d83fb.png](https://i-blog.csdnimg.cn/blog_migrate/dffd0c24c0104d1dd44ec438e28e1882.png)
其中,C10,C01,C11,和D1是材料参数。
Mooney-Rivlin 5参数模型
![075a45afcc43d82331acb89cd084f459.png](https://i-blog.csdnimg.cn/blog_migrate/72c5b212bccbb0513a677fdf5d6fdcea.png)
其中,C10,C01,C20,C11,C02,和D1是材料参数。
Mooney-Rivlin 9参数模型
![db45d26d0efbc5421eb71e18565cfb1f.png](https://i-blog.csdnimg.cn/blog_migrate/67e7a36b63c84566d01c9ae8813bd160.jpeg)
其中,C10,C01,C20,C11,C02,C30,C21,C12,C03 和D1是材料参数。
Neo-Hookean模型
![86f5c29233aad0b05ca04d080f7d7027.png](https://i-blog.csdnimg.cn/blog_migrate/81ef10670d00a3fbb0dcd67e5cb3ec99.png)
其中,u是初始剪切模量,D1是材料不可压缩参数。
Ogden模型
![e45ef89d0e4e5fa90e7541532a4c7845.png](https://i-blog.csdnimg.cn/blog_migrate/0f2b369ee2dc29fd9dde61f4d3d64ad5.png)
其中,N确定多项式阶数,u和a是材料参数,D是不可压缩参数。
多项式模型
![997f797afae0d46bf7c54e26a41a7f89.png](https://i-blog.csdnimg.cn/blog_migrate/33f62100cffd6ea37527a8a1b03f2c3d.png)
其中,N确定多项式阶数,c_ij和D是材料参数。
Yeoh模型
![4a99000d3d28a53c99b14f4b0e143d14.png](https://i-blog.csdnimg.cn/blog_migrate/b22ab39f78d1eb30442578e120a5eb5b.png)
其中,N确定多项式阶数,c_i0和D是材料参数。
根据具体条件,使用者可以通过结合测量得到的单轴拉伸、压缩、双轴拉伸和体积测试数据来更好地估算材料参数。这些测试数据随后可以对每个超弹模型进行曲线拟合。
给定测试数据和具体超弹模型后,材料常数可以通过最小二乘法计算残差得到最优的参数结果,常用Levenberg-Marquardt 求解器来实现最小化算法。拟合曲线是目前找出超弹性材料模型中材料参数最主要的方法之一;但还应考虑给定超弹性材料模型的稳定性。通常采用Drucker稳定性准则来确定材料的稳定性。根据Drucker准则,与增量应力相关的应变能应大于零。如果违反了该准则,材料模型将会不稳定。在实际数值计算中,常用的方法是检查材料刚度矩阵的正定性,以此来确定给出的材料常数是否稳定。
使用MatEditor进行超弹材料曲线拟合
工程材料编辑软件MatEditor已经支持了多种超弹模型材料及其参数拟合。通过输入应力-应变试验数据,就可以得到用于有限元分析的超弹材料常数。MatEditor不仅能够使用试验数据拟合出所选本构函数(应变能函数)的参数,而且还能将本构函数曲线与试验数据(名义应力-应变曲线)绘制在同一图表中,便于对比拟合效果。这样,使用者可以利用实验应力应变数据直接让程序自己拟合出超弹材料模型的参数。
目前MatEditor和WelSim已经支持的超弹曲线拟合模型有:Arruda-Boyce, Mooney-Rivlin, Neo-Hookean,Ogden, Polynomial,和Yeoh模型。本例使用经典的Mooney-Rivlin 9 作为示例,曲线拟合基本步骤如下:
1. 选择超弹性材料
从超弹列表中选择Mooney-Rivlin 9材料属性,添加至材料属性编辑窗口。此时的参数为空白需要使用者输入或者参数拟合。
![7bc4a4a712117f1f82966845f0ba0e99.png](https://i-blog.csdnimg.cn/blog_migrate/281f5ad8cf4c574e0f002ab533dfff2b.jpeg)
2. 从测试数据列表中选择添加单轴、双轴、和剪切三种试验数据,并在表格中输入测试数值,如下图所示单轴拉伸试验数据。MatEditor也支持多个温度下的测试数据。测试数据需要覆盖后续仿真的应变范围。如果数据较多,使用者可以通过文件导入的方式输入应力-应变数据。
![dd61334cf5b2233ba421ad5b2465eb53.png](https://i-blog.csdnimg.cn/blog_migrate/7491807c5ef74e0d846d81ddfc3d78c5.jpeg)
3. 确认测试数据输入正确后。右键点击Mooney-Rivlin 9属性,添加曲线拟合属性。
![51e9d92a99d237f168e537e4df64145a.png](https://i-blog.csdnimg.cn/blog_migrate/5d7b4a589c2465bc72a245bcfa0f8365.png)
4. 右键点击刚添加的曲线拟合属性,计算,并拷贝计算值至材料属性中。
![50c130998e51b953f7e3f7cbba3840bc.png](https://i-blog.csdnimg.cn/blog_migrate/a758d4c22588b376e49ce3285894043d.jpeg)
5. 如果计算并赋值成功,可以查看拟合后的材料参数和曲线状况。
![a0110c828ff97420c468c6f2b1524ff9.png](https://i-blog.csdnimg.cn/blog_migrate/af0c5584adf8f38ff8de3ed2f978525d.jpeg)
6. 对比曲线和测试数据的差异。并决定是否采用这些参数进行后续的分析工作。本例中,已经计算出了Mooney-Rivlin的9个参数。而体积数据用来计算最后一个数据d(不可压缩性因子,与泊松比有关),可以通过体积实验曲线拟合得到。如果没有体积试验曲线,则d=0,这种情况下可能需要手工修改d值,如果知道泊松比u,则可用如下公式计算:d=(1-2u)/(c1+c2)。注意此式是建立在几乎不可压缩(u接近于或等于0.5)的前提下的。
现给出同一组数据下,不同模型曲线拟合的结果。可以直观地理解不同超弹模型曲线拟合的特点。
![69df2065a39bc411e4513c7462422656.png](https://i-blog.csdnimg.cn/blog_migrate/7adb1311b9fb1bffa17b56649dcd26ca.jpeg)
![a7f0e7ea007807ccace64a29eceb792e.png](https://i-blog.csdnimg.cn/blog_migrate/a4584727af06adf16e8cf74549b6abe9.jpeg)
![03fb1d7b4ecc243216ee5150ecd2a9ca.png](https://i-blog.csdnimg.cn/blog_migrate/619e48d5b6567ad7e0defd0aa062f281.jpeg)
![813fdbb5dab638baa696d2d8ca9aca50.png](https://i-blog.csdnimg.cn/blog_migrate/e74c85d49b6156e3934ff8f744300e1f.jpeg)
![e43c0521c143e5cf09fd25f3b4ed7fde.png](https://i-blog.csdnimg.cn/blog_migrate/ec704d1fb3e881ea959a49a292a5890c.jpeg)
![4b5555d488bef3166d87f337a0a0cb1f.png](https://i-blog.csdnimg.cn/blog_migrate/703a3ce85281112ca768c38432f34f5f.jpeg)
![943f7a8f3b540b502dbac8bc0dcf80f7.png](https://i-blog.csdnimg.cn/blog_migrate/26cfcb3cfa398ac2492530ebe50846e5.jpeg)
![4519db1ee92657dcd8c1e935ff74c872.png](https://i-blog.csdnimg.cn/blog_migrate/d69762a851ab95f78fbae8f74e65d5e1.jpeg)
![a9eb35f587bc4aa1bbbdb37af85b7f04.png](https://i-blog.csdnimg.cn/blog_migrate/c33979e334f0744fca8af3b25c40ee52.jpeg)
![838fb2c8d02319216f7617efbb63d29a.png](https://i-blog.csdnimg.cn/blog_migrate/525a999d055e50db1182d88ae44c59d6.jpeg)
![17edfa235f929e805b41c41326955a03.png](https://i-blog.csdnimg.cn/blog_migrate/2c889ec0d9c7691eb1b26d9b1650bfbb.jpeg)
![9982f51bf09d6b54c93e5236bb7ac080.png](https://i-blog.csdnimg.cn/blog_migrate/6489f61655d29dfb23f5e1d6ebf51f83.jpeg)
超弹模型曲线拟合中的注意事项
1. 若要得到精确的材料常数,曲线拟合中输入的应力-应变数据必须涵盖分析中将遇到的完整载荷范围。例如,如果部件除了拉伸和压缩以外,还会遇到剪切应力,则必须输入纯剪切应力-应变数据。简单拉伸本身不足以对这些情况下的材料行为进行建模。如果只拟合单轴数据,并使用计算的参数用于实际双轴变形分析,可能会得到错误的有限元分析结果。在材料参数估计中,最好能结合不同的大变形模式执行曲线拟合,而非只使用一种变形模式。
2. 在大多数情况下,曲线拟合例程无法达到完美的精度。曲线拟合的结果很大程度上取决于所需范围内给定的测试数据。
3. 大多数的超弹模型是基于完全不可压缩理论的(除了Blatz-Ko和Ogden Foam等),所以对于这些模型而言,输入参数时,要注意其破松比应该接近0.5,一般大于0.49。这样才能在有限元计算中得到合理的结果。
本文中,我们讨论了各种用于橡胶材料的常用非线性超弹模型及其曲线拟合。为了将材料参数拟合到材料模型中,我们需要正确的测量数据。我们还分析了一些典型的实验测试、材料模型选择。通过在MatEditor软件示例,演示了如何在非线性超弹模型中直接使用测量数据,并拟合到不同超弹性材料模型的参数。
最后,给出相关软件操作视频,供大家参考。
知乎视频www.zhihu.com