橡胶工业中有限元计算问题过盈配合.docx

作者：清华大学工程力学系范成业摘要本文分析了过盈配合的有限元计算时用到超弹性本构时可压缩性对计算结果的影响情 况，得到在过盈配合中必须考虑这种可压缩性的结论并分析考虑可压缩性的原因1、引言过盈配合是橡胶工业中的一种常见的配合方式橡胶为超弹性材料，有限元计算中通常假 定为不可压或者几乎不可压本文首先给出一种不可压橡胶模型过盈配合的理论解，并与 ABAQUS 计算解进行比较进一步本文探讨过盈配合中假定橡胶不可压时遇到的问题，提 出处理过盈配合中橡胶计算的方法2、可压模型理论解与ABAQUS数值解的比较2.1、理论解理论解模型如图1，内层为钢，中间不可压橡胶，最外层为钢给出橡胶和橡胶之间的过盈 量求整个结构的应力应变状态假设平面应变状态本构方程：对于钢：对于橡胶：2.1 材料性质：钢：E=210000v=0.3橡胶：C10=0.461312, C20=0.01752, C30=8.8e-05,其余为0,（三次多项式模型，材料不可压缩）2.2.2几何特性如图2所示，R59.50为内层钢的半径和中间层橡胶的内径，R73.00为中间层橡胶的外径,R71.10为外层钢的内径，R80.00为外层钢的外径。

图 2 不可压模型算例几何特征理论解与计算解的比较（理论解由Maple计算得出）表 1 理论解与 ABAQUS 解的比较半径（mm）理论解ABAQUS计算解误差位移59.5-9.2984E-02-9.73152E-24.6%径向应力S1159.5 （钢）-660.51-631.60-4.38%59.5 （橡胶）-660.51-631.60-4.38%73.0 （橡胶）-660.51-631.60-4.38%71.1 （钢）-660.51-664.300.57%80.0 （钢）028.15一环向应力S2259.5 （钢）-660.51-631.20-4.44%59.5 （橡胶）-660.51-631.20-4.38%73.0 （橡胶）-660.51-631.40-4.38%71.1 （钢）5626.365541.00-1.52%80.0 （钢）4956.854957.000.00%3、可压缩模型橡胶的应变能采用多项式模型时，在静水压力荷载下p与J的关系如下:用 ABAQUS 对这 1-4 组系数进行评估：⑷路」跖WCEON图 3 不同系数对应的橡胶静水压力下的应力应变关系将这六种橡胶本构代入第二部分中的算例中进行计算结果如下：图 4 第 6 组系数对应的位移图图 5 第 1 组系数对应的位移图由图4 和图5容易看到这两组系数对应的位移差异非常大。

表 3. 不同系数对过盈配合的影响-橡胶外侧位移U1外部钢环的内侧位移U1 (mm)过盈面上的应力S11 (MPa)橡胶内最大S22(MPa)橡胶内最大S33(MPa)1-1.8663.36E-02-12.175-11.94-11.92-1.602.95E-02-106.2-106.3-106.03-0.71.21-422.5-422.5-422.54-0.171.735-602-602-6025-0.111.79-621.0-621.0-621.06-9.6E-21.825-631.7-632.0-6314、结论由以上分析可以看到橡胶的可压缩性在过盈配合中起到的作用是非常大的例如，使用初始Poisson's Ratio为0.499时，橡胶相对于钢是非常软的，因为这组系数确 定的橡胶的本构模型杨氏模量大约为3MPa,体积模量(ABAQUS对系数的评估给出的结 果)，而钢的杨氏模量为210000MPa，取时因此，如果用实验方法来确定橡胶应变能的形式，应对橡胶进行三向应力实验，以确定橡 胶应变能中与体积应变有关的各个系数参考文献[1] ABAQUS/Thoery 用户手册 V6.1 HKS.INC 2000(end)。