01 引子
橡胶材料是典型的超弹性材料,要获取超弹性材料本构模型(常见有Mooney-Rivlin、Ogden、Yeoh等),一般需要做一系列标准橡胶试验并进行数据拟合。
本例演示了ANSYS对超弹性材料的曲线拟合能力,并通过有限元分析与拉扭试验的对比,验证所建立的本构模型的有效性。
常见的橡胶标准拉伸试验
02 案例介绍
现需要一个本构模型来匹配硫化天然橡胶材料在各种变形模式下的100%工程应变的行为。
本例中,已通过试验(单轴、双轴和平面拉伸试验)获取了橡胶的实验数据。使用这些数据,通过超弹性拟合能力确定本构模型的参数,可以拟合3参、5参和9参的Mooney-Rivlin超弹性模型。
试验数据
同时对橡胶进行了拉扭实验(将条形试件的两端夹入测试仪器中,然后将试样拉伸到原尺寸长度的50%,并将试样的一端扭四圈)。试样与ASTM D1043中规定的试样相似,如下图所示:
拉扭试验条形试件
使用拟合得出的Mooney-Rivlin超弹性模型(5参为例)对拉扭试验就行有限元分析,并与试验结果相对比,据此判断前面拟合得出的本构模型能否反映橡胶材料的真实行为。
模型采用SOLID186单元,两端夹钳区域采用MPC算法绑定到定位点。
有限元模型示意图
按照拉扭试验的加载顺序:
step1:对两端夹持区域施加试件厚度25%的压缩位移,模拟夹具对试件的夹持作用。
step2:通过移动一侧的夹持区域(刚性接触面),同时固定另一侧夹持区域,模拟拉伸到50%的拉伸状况。
step3:通过将一侧夹持区域保持在原位并围绕纵轴旋转另一侧夹持区域绑定的关键点来模拟试件的扭转。
加载顺序
03 关键命令流
通过*create命令创建data文件存储试验数据:
!!!!!生成data文件并存储单轴拉伸试验数据
*create,uniaxial.data
0.342582 2.94685
0.525656 3.90358
0.712876 4.67940
0.902195 5.35192
1.091940 5.98056
*end
使用tbft命令拟合数据:
!定义1号材料为超弹性mooney模型(5参)
tbft,fadd,1,hyper,mooney,5
! 读入(前面定义好的)单轴拉伸试验数据
tbft,eadd,1,unia,uniaxial.data
!读入双轴拉伸试验数据
tbft,eadd,1,biax,equibiax.data
!读入平面试验数据
tbft,eadd,1,shea,planar.data
!进行拟合q
tbft,solve,1,hyper,mooney,5
!应用拟合结果
tbft,fset,1,hyper,mooney,5
tblist,all,all !弹窗列出拟合结果
将夹持面用MPC算法绑定到关键点上:
R,11
real,11
et,11,174
keyopt,11,12,5 ! 绑定接触
keyopt,11,4,2
keyopt,11,2,2 ! MPC算法
et,12,170
keyopt,12,2,1
cmsel,s,fr !选择夹持面
nsla,,1
real,11
type,11
esurf !夹持面覆盖接触单元
type,12
kmesh,10001 !关键点10001定为参考点
04 计算结果
计算模拟出试件夹持、拉伸、扭转行为:
提取扭转过程的扭矩-角度(M-θ)曲线,与试验结果较为吻合,证明前面通过拟合得出的5参Mooney模型可准确预测该橡胶行为。
免责声明:本文系网络转载或改编,未找到原创作者,版权归原作者所有。如涉及版权,请联系删