基于非线性有限元的橡胶隔振器压缩刚度的研究

2014年11月机械设计与制造工程NOV2014第43卷第11期MACHINEDESIGNANDMANUFACTURINGENGINEERINGVOL43NO11DOI103969/JISSN2095509X201411016基于非线性有限元的橡胶隔振器压缩刚度的研究秦洪艳三江学院机械工程学院，江苏南京210012摘要借助ANSYS软件完成了压缩状态下矩形橡胶隔振器的非线性有限元分析，总结了矩形橡胶隔振器压缩刚度计算公式中形状系数的确定方法，并与传统经验公式进行了比较。关键词矩形橡胶隔振器形状系数非线性有限元中图分类号TB535文献标识码A文章编号2095509X201411007403由于橡胶材料具有良好的弹性性能和缓冲功能，工程机械设计中通常采用橡胶隔振器达到支承、限位、隔振及减噪的目的。为了方便设计计算，通常将复杂形状的橡胶隔振器看成是几个简单形体的叠加。国内外对典型形状的橡胶隔振器刚度计算方法进行了研究，一些机构也总结出计算典型橡胶隔振器刚度的经验公式。但橡胶属于典型的非线性弹性材料，其使用性能随橡胶配方、使用环境、隔振器造型及受力变形量而变化，因此传统的经验公式具有一定的局限性。面积AL与自由面积AF的比值N有关，如图1所示。形状系数的表达式为MFN，NAL/AF3本文借助有限元分析软件ANSYS对压缩状态下的矩形橡胶隔振器进行非线性有限元分析，从而得到矩形橡胶隔振器压缩刚度计算公式中形状系数的表达式，此表达式同样适用于其他形状的橡胶隔振器。1矩形橡胶隔振器压缩刚度的计算橡胶作为一种非线性弹性材料，其基本的物理性能不会由于外表形状的改变而改变，但制成隔振器后它的弹性模量不再是橡胶材料的弹性模量E，而是表观弹性模量EA。表观弹性模量与弹性模量图1矩形橡胶块对于矩形橡胶隔振器来说NLB2LBH式中L为矩形隔振器的长B为矩形隔振器的宽。只要确定了形状系数MFN的表达式，橡胶隔振器的压缩刚度可以通过式2计算。而橡胶材料的形状、硬度、变形量都会对形状系数M产生影响。上海橡胶研究所总结的形状系数经验公式是的关系为M115N2N225N34EAME11适用于N02时日本机械学会总结的形状系数经验公式是因此，橡胶隔振器压缩刚度KX的计算公式为ALMM1219N25KXHE2适用于正方体，N较小时式中AL为橡胶块的约束面积H为矩形橡胶块的M1039N2273N36高度M为形状系数，其大小与橡胶隔振器的约束适用于正方体，N较大时收稿日期20141020作者简介秦洪艳1981，女，山东潍坊人，三江学院讲师，硕士，主要从事车辆振动方面的研究。2014年第11期秦洪艳基于非线性有限元的橡胶隔振器压缩刚度的研究M1331093N2适用于长度远大于宽度的矩形体72利用有限元法计算矩形橡胶隔振器的形状系数对标准橡胶试样进行单轴拉伸和压缩试验，可以得到该橡胶材料的应力应变曲线，图2是邵氏硬度HA65时的橡胶试样的平均应力应变曲线，该曲线用于计算隔振器的弹性模量E，以及确定橡胶材料的MOONEYIVLIN常数。图2HA65的橡胶试样的平均应力应变曲线借助ANSYS软件对压缩状态下的矩形橡胶隔振器进行非线性有限元分析，可得到某一压缩量下的载荷位移关系，具体步骤如下A在ANSYS软件中建立矩形橡胶隔振器的有限元模型，包括定义单元类型、创建节点和单元，同时定义材料性质等。B输入对橡胶材料标准试样的实测实验数据，通过橡胶材料的应力应变曲线确定MOONEYIVLIN常数的个数，如果橡胶材料的应力应变曲线没有拐点反弯点，则MOONEYIVLIN常数为2个有1个拐点，MOONEYIVLIN常数为5个有2个拐点，MOONEYIVLIN常数为9个。图2所示的应力应变曲线只有1个拐点，因此MOONEYIVLIN常数为5个，借助ANSYS软件计算MOONEYIVLIN常数2。C对矩形橡胶隔振器进行加载，设置载荷步和各项参数。D求解并进行结果的后处理。图3为ANSYS后处理器中显示的矩形橡胶隔振器压缩量为10时的变形情况。通过时间历程后处理器定义了隔振器顶面位移和顶面承受载荷2个变量，得到图4所示的载荷位移关系曲线。图3压缩量在10的矩形橡胶隔振器的变形情况图4载荷位移关系曲线通过载荷位移曲线可以求橡胶隔振器的表观弹性模量EA，该压缩量下，橡胶材料的弹性模量E可以通过拉伸压缩试验得到的应力应变曲线进行求解，并利用公式1求解压缩刚度的形状系数M。改变隔振器的外形尺寸，可以得到若干组MFN的对应关系。橡胶隔振器用于减振设计时，通常压缩变形量控制在15以下。分别计算压缩量为5，10，12和15时橡胶隔振器的刚度特性，得到了HA65时的矩形橡胶隔振器在不同压缩量下的表观弹性模量EA，见表1。表1HA65时橡胶材料在不同压缩量下的表观弹性模量EA压缩量5101215EA/MPA621643650657遵循前述方法与步骤得到多组MN的对应关系。借助OIGIN或MATLAB软件，利用最小二乘法求解形状系数表达式MABNCN2，得到不同压缩量下的MN关系，见表2。752014年第43卷机械设计与制造工程表2HA65时矩形橡胶隔振器在不同压缩量时的MN对应关系压缩量拟合公式5M1062N203N210M1045N246N212M1041N261N215M102N286N2将非线性有限元法获得的矩形橡胶隔振器在不同压缩量下形状系数的表达式与经验公式4和经验公式5的计算结果进行比较，结果如图5所示，横坐标为N，纵坐标为形状系数M。很明显，由日本机械协会的经验公式5计算的形状系数M的值误差较大，不适合此隔振器当02N035时，由上海橡胶研究所的经验公式4计算的形状系数M的值误差较小，但随着N值的增大，误差也越来越大。另外，压缩变形量在15以下时，图5非线性有限元法得到的形状系数同传统经验公式计算结果的比较该隔振器压缩刚度形状系数变化不大。为了简化计算，在误差允许范围内，取压缩量为10的形状系数的表达式M1045N246N2作为压缩量在15以下的矩形橡胶隔振器压缩刚度形状系数的通用表达式。3结束语橡胶隔振器的使用性能受橡胶的配方、硬度、温度以及隔振器的外形尺寸、施加载荷、变形量等许多因素影响3。传统橡胶隔振器压缩刚度的形状系数是通过对类似模型进行试验得到，然后再进行产品试验，纠正设计计算的误差，这种方法得到的经验公式不仅误差大，而且浪费大量的人力物力。目前ANSYS、ADINA、ABAQUS等有限元软件已具备橡胶超弹性特性材料和体积近似不可压缩材料的大变形有限元计算功能，借助这些专用软件可以实现橡胶隔振器的有限元仿真设计。应用有限元软件对橡胶隔振器进行静态和动态特性分析，不仅可以得到较为精确的压缩刚度，而且可以实现产品的优化设计，既提高了隔振器的设计质量，又减少了试验次数，节约了试验耗材。参考文献1严济宽机械振动隔离技术M上海上海科学技术文献出版社，19852772902王利荣，吕振华橡胶隔振器有限元建模技术及静态弹性特性分析J汽车工程，200264804853秦洪艳，张红圆柱形橡胶垫刚度经验公式的研究J青岛农业大学学报，200735961ESEARCHONCOMPRESSIONSTIFFNESSOFRUBBERISOLATORBASEDONNONLINEARFINITEELEMENTQINHONGYANSANJIANGUNIVERSITY，JIANGSUNANJING，210012，CHINAABSTRACTITINTRODUCESTHENONLINEARFINITEELEMENTANALYSISUSEDTOTHECOMPRESSEDRECTANGULARRUBBERISOLATORWITHTHEHELPOFANSYSSOFTWARE，DESCRIBESTHECOMPUTI