基于本构模型的橡胶隔振器静态力-位移特性分析

基于本构模型的橡胶隔振器静态力-位移特性分析

橡胶材料作为一种高性能的超弹性材料，广泛应用于汽车耳盖的振动装置中，如汽车动力的总悬挂、车身的悬挂、气管的悬挂、底盘套筒等。橡胶材料是一种大变形的非线性超弹性材料,对其进行力-位移分析时,其本构模型的选择和本构模型常数的获取是一项重要的工作。目前可用于橡胶隔振器静态力-位移关系计算的本构模型的较多,如Mooney-Rivlin、Ogden等。不同本构模型的常数不一样,如何合理的选择本构模型和开展哪些实验以获取本构模型的常数是橡胶隔振器静态力-位移分析的重要工作。Seibert等对Arruda-Boyce模型、VanDerWaals模型、Yeoh模型三种本构模型进行对比研究。黄建龙等利用Mooney-Rivlin模型和Yeoh模型对橡胶材料进行有限元分析,介绍两种本构模型材料常数的获取方法。王丽荣等利用Mooney-Rivlin模型和Ogden模型对橡胶隔振器静态力-位移特性进行有限元分析,介绍了单元特性和网格质量对橡胶隔振器静态特性的影响。Gracia等研究了一种应用于工业橡胶材料弹塑性特性分析的overlay模型,并介绍其常数拟合过程。StefanHartmann等研究了一种获取橡胶本构模型材料常数的新方法-光学测量法,确定采用光学测量法获取材料常数的准确性。本文选取了目前广泛应用的几种本构模型来研究橡胶隔振器力-位移关系计算结果的影响因素。为获得各本构模型的材料常数,开展了橡胶材料在不同应力-应变状态(单轴拉伸、等双轴拉伸和平面剪切)和不同最大应变下的应力-应变关系的测试。利用所测得的应力-应变关系,和不同本构模型得到的理论应力-应变关系,根据最小二乘法获得不同本构模型中的材料常数。利用获得的材料常数,对两种橡胶隔振器(橡胶悬置1,2)的静态力-位移特性进行了计算并和实测值进行了对比分析。分析结果表明,利用同一本构模型计算隔振器的静刚度时,对悬置1,在拉压变形的方向,采用三种应力-应变状态拟合得到的本构模型常数计算得到静刚度值,比单轴拉伸应力-应变状态下的计算精度高;对只有剪切变形的方向,用平面剪切应力-应变状态拟合得到的本构模型常数计算得到的静刚度,比三种应力-应变状态下的计算精度高。对悬置2,在同一本构模型下,利用单轴拉伸、等双轴拉伸和平面剪切三种应力-应变状态拟合得到的本构模型常数计算出的静刚度值,比在两种应力-应变状态下计算出的静刚度值相对误差更小。以橡胶悬置2为研究对象,分析了不同本构模型对橡胶隔振器静态力-位移特性计算结果的影响。结果表明,Mooney-Rivlin模型的计算精度最高,其相对误差均小于10%,其次是Arruda-Boyce模型,计算相对误差最大的是Marlow模型。1本构模型及其常数对橡胶隔振器进行静态特性计算时,假定橡胶材料为各向同性的超弹性材料。通过对其应变能密度函数对应变不变量求导,可以得到材料的工程应力与工程应变之间的本构关系。工程应力和工程应变是忽略因施加载荷的增加或者减少而引起的横截面面积的变化所得到的应力、应变。超弹性材料的应变能密度函数有多种形式,如Mooney-Rivlin模型、VanDerWaals模型、Marlow模型、Ogden模型、Yeoh模型、Arruda-Boyce模型、Neo-Hookean模型、Ploynomial模型等。本文只讨论目前广泛应用的前六种本构模型,及其常数的获取。应变能密度函数的一般表达式为:其中,I1、I2、I3分别为一阶、二阶、三阶应变不变量,它们为三个主拉伸比的函数;Ci(i=1,2,…,m)为m个表示超弹性材料剪切特性的常数,dj(j=1,2,…,n)为n个表示超弹性材料压缩特性的常数。I1、I2、I3与超弹性材料的三个主拉伸比λ1、λ2、λ3的关系为:由文献、文献[9-10]中对应的各种应变能密度函数,对三个主拉伸比λ1、λ2、λ3分别求导数,就可以计算出材料在不同变形状态时的工程应力:其中λU、λB和λP分别为测试得到的单轴拉伸、等双轴拉伸和平面剪切收缩率。三种不同的变形方式的应变不变量(I1,I2)与主拉伸比之间的关系,可由式(2)求得。单等轴双拉轴伸拉:伸I:1I1==λ2U2λ+2B2+λU-λ1B-4,,I2I2==λU2-2λB-+22+λUλ4B}(4)平面剪切:I1=λP2+λP-1+1,I2=I1采用最小二乘法拟合实测的应力-应变和由本构模型计算得到的应力-应变,可以求得各本构模型的材料常数:其中,σki为实验测试的工程应力,j为实验类型:单轴拉伸时,j=1;平面剪切时,j=2;等双轴拉伸时j=3。S为计算应力与实验测试应力的相对误差。2橡胶材料应力-应变本构模型根据上节所述的不同本构模型的材料常数获取方法,在对橡胶材料进行应力-应变测试时,对橡胶试件进行单轴拉伸、等双轴拉伸以及平面剪切实验,以获得橡胶材料在不同应力-应变状态和不同最大应变下的应力-应变关系。橡胶材料的应力-应变实验是委托美国AxelProducts,Inc测试的。测试用的橡胶材料的邵氏硬度为50。图2为在不同应力-应变状态和不同最大应变下,实测的工程应力-工程应变关系。在同一应力-应变状态,不同的最大工程应变时,工程应力-工程应变曲线是不重合的,在加载和卸载过程中的应力-应变曲线也是不重合的。这是由于橡胶材料存在Mullins效应的作用:即在对橡胶试件进行实验过程中,在加载和卸载中应力会出现变小的现象。在进行橡胶材料的应力-应变测试时,为保证计算精度以及减小Mullins效应的影响,取循环加载(循环次数一般为5次)的最后一次应力-应变数据为实测的工程应力与工程应变的关系。3模型模型材料的常数3.1等双轴拉伸应变的最大限度由于等双轴拉伸等效单轴压缩,且本文计算的橡胶隔振器的总体变形的应变水平小于0.5,故等双轴拉伸应变的最大值仅取1。为了研究应变状态和最大应变对橡胶隔振器的力-位移计算结果的影响,选取其中不同应力-应变状态和最大应变进行组合,以得到同一本构模型在不同应力-应变状态和最大应变下的本构模型常数。不同应力-应变状态和最大应变组合如表1所示。3.2拉伸应力-应变状态下的拟合效果由图3可见:(1)在单轴拉伸应力-应变状态下,1号、3号组合时,Mooney-Rivlin模型拟合效果最好,Marlow模型拟合测试曲线趋势较好,VanDerWaals模型拟合效果相对较差,其他三种本构模型拟合材料常数效果较好。(2)在等双轴拉伸应力-应变状态下MooneyRivlin模型拟合效果最好,其次是VanDerWaals模型,Ogden模型和Yeoh模型相对于Mooney-Rivlin模型来说,拟合效果相对较差,拟合效果最差的是Marlow模型,特别在3号组合时,Marlow模型拟合曲线基本偏离实测应力-应变曲线。(3)在平面剪切应力-应变状态下,Mooney-Rivlin模型拟合误差最小,其次是Arruda-Boyce模型,拟合效果相对较差的是Marlow模型。综合以上分析结可知,Mooney-Rivlin模型在单轴拉伸、等双轴拉伸和平面剪切三种应力-应变状态下的拟合效果最好,Arruda-Boyce模型和VanDerWaals模型次之,拟合效果最差的是Marlow模型。表2列出了在表1中的九种应力-应变状态和最大应变组合下,各本构模型的材料常数。4该模型常数对橡胶横截面的机械强度评价有影响本文选取两种不同类型的悬置做为研究对象。悬置1主要承受拉压或者剪切变形,而悬置2则同时承受拉压和剪切变形。4.1悬置本构模型和材料常数的确定图4所示的橡胶悬置在X向只承受拉压、Y向只承受剪切。在对此橡胶悬置进行有限元分析时,在X、Y、Z三向均加载5mm的位移,计算所需的外力。用Mooney-Rivlin本构模型,利用表1中的1号、2号、3号、4号、5号和7号、8号、9号应力-应变状态组合和不同的最大应变状态得到的本构模型材料常数,对悬置1进行三向静刚度的有限元计算。计算得到的悬置1在X、Y、Z三个方向的刚度和实测值的对比见表3~表5。由表3可见,在悬置只承受拉压变形的X向,利用1号、2号、3号组合拟合得到的本构模型常数,计算得到的静刚度,比只在单轴拉伸应力-应变状态(7号)下计算的静刚度相对误差小。由表4可见,在悬置只受剪切变形的Y向,利用三种应力-应变状态拟合得到的材料常数,计算得到的静刚度,比只在平面剪切应力-应变状态(9号)下计算得到的静刚度相对误差大。由表5可见,在同时承受压缩和剪切应力的Z向,在单轴、等双轴拉伸两种应力-应变状态下计算的静刚度相对误差最小,而在实测等双轴拉伸应力-应变状态(8号组合)下计算得到的静刚度相对误差最大。4.2本构模型的检验和计算图5所示的悬置2橡胶主簧同时承受拉压变形和剪切变形。在进行力-位移关系计算时,X向加载力的范围从-3500~3500N;Y向加载力的范围从-500~500N;Z向加载力的范围从-1000~1000N。(1)利用Mooney-Riv-lin本构模型,和不同实测应力-应变组合下拟合得到的本构模型材料常数,计算其三向力-位移曲线,计算结果见图6。在X向和Z向,Mooney-Rivlin模型在1号应力-应变组合下计算得到静刚度值与实测值的相对误差较小;在Y向,在3号应力-应变组合下计算的静刚度与实测值的误差较大。(2)利用Arruda-Boyce模型、VanderWaals模型、Ogden模型和Yeoh模型,选取4组应力-应变状态组合,对悬置2进行三向静刚度计算。计算结果分别见表6、表7、表8和表9。由表6~表9可见,由三种应力-应变状态下计算的静刚度值,比用两种应力-应变状态下计算得到的静刚度值更准确。5几种本构模型的计算结果比较在同种应力-应变状态组合下,用不同本构模型对悬置2进行静刚度计算并与实测值对比。选取1号、3号应力-应变组合下的计算结果进行分析,如图6所示。由图6可见,用Mooney-Rivlin模型计算的悬置在X、Y和Z三向静刚度计算值与测试值最接近,其相对误差小于10%。用VanderWaals模型和Marlow模型计算得到的静刚度与实测值相对误差较大;利用Arruda-Boyce计算得到的静刚度与实测值吻合的较好,计算精度较高,大部分的静刚度计算值与实测值的相对误差小于10%。在3号应力-应变组合下,Ogden模型和Yeoh模型计算得到的静刚度与实测值相对误差小于10%。由此可见,对于研究的几种本构模型,利用Mooney-Rivlin模型计算得到的悬置的静刚度计算值与实测值的相对误差最小。其次是Arruda-Boyce模型,误差最大的是Marlow模型和VanderWaals模型。6本构模型的选取根据以上计算结果,可得到如下结论:(1)对于只受拉压变形的橡胶隔振器,用三种应力-应变状态比只用单轴拉伸应力-应变状态下拟合得到的本构模型常数计算得到的静刚度计算值精度高。对于只受剪切变形的橡胶隔振器,用平面剪切应力-应变状态比用三种应力-应变状态拟合得到的本构模型材料常数计算得到的静刚度精度高。(2)对橡胶隔振器进行力-位移特性分析时,需选取合适的超弹性橡胶本构模型。本文的计算结果表明Mooney-Rivlin模型的计算