基于ABAQUS的橡胶减振器非线性有限元分析

0 引言

    橡胶减振器因其阻尼性能好、结构简单、价格低廉等优点，在航天、船舶、汽车等工程领域得到广泛应用。但由于橡胶的复杂的非线性力学行为，缺乏准确的材料模型来描述其特性，使得橡胶减振器的选型和设计主要依靠工程经验和试验方法，耗费大量人力物力。随着有限元技术和计算机技术的发展，作为先进的非线性有限元软件一ABAQUS已经具备超弹性、粘弹性橡胶模型和体积不可压缩材料大变形的有限元计算功能，这为橡胶减振器的有限元分析提供了极大的便利。

1 橡胶减振器非线性有限元分析

1.1 橡胶减振器非线性特性

    橡胶减振器非线性主要来源于几何非线性、材料非线性、边界非线性。其中几何非线性指橡胶为体积不可压缩材料，其在工况中会出现大弹性形变，已不能用小形变理论分析；材料非线性指橡胶材料为粘弹性材料，其应力-应变曲线为非线性的，且动态特性与频率、振幅、温度相关；边界非线性指橡胶材料与金属的接触过程中，其边界条件在分析过程中发生变化。在ABAQUS中，几何非线性即大变形可在分析步(Step)中选择开启，开启后在分析过程中会考虑到几何非线性;材料库中有丰富的材料可供选择，对于橡胶材料可选择超弹性、粘弹性；ABAQUS中在相互作用（Interaction)内定义各部分之间的相互作用，对于边界非线性可定义橡胶与金属支架之间的接触，通过接触算法实现边界条件非线性和力的传递。

1.2 橡胶材料参数的确定

    在有限元方法中，常把橡胶近似视为不可压缩材料，使用超弹性本构模型描述其弹性能力，其中的Mooney-Rivlin本构模型如下

    其中，I1、I2、I3是应变不变量，C10、C01、D1是由材料决定的常数。当D1=0时，橡胶完全不可压缩。Mooney-Rivlin模型可以仅使用C10、C01和D1三个常数，较精确地刻画在橡胶元件工作于中小应变时的弹性能力。

    材料常数C10、C01和D1的确定，需要进行一系列的复杂试验。由超弹性材料常数C10/C01与橡胶的邵氏硬度HS有如图2的经验关系。而对于可自由变形的橡胶元件，其剪切模量G(MPa)与材料常数C10、C01有关系：

    则对于选定硬度的橡胶材料，通过图确定C10/C01的值，在代入式(2)中，可以求解出材料常Cl0、C01。

    图1 材料常数与邵氏硬度的经验关系

    本文所用橡胶减振器橡胶为某型号硅橡胶，其弹性剪切模量为2MPa，邵氏硬度为50。通过经验公式得到材料参数如表1。