|
摘要本文分析了过盈配合的有限元计算时用到超弹性本构时可压缩性对计算结果的影响情况，得到在过盈配合中必须考虑这种可压缩性的结论并分析考虑可压缩性的原因。
1、引言
过盈配合是橡胶工业中的一种常见的配合方式。橡胶为超弹性材料，有限元计算中通常假定为不可压或者几乎不可压。本文首先给出一种不可压橡胶模型过盈配合的理论解，并与ABAQUS计算解进行比较。进一步本文探讨过盈配合中假定橡胶不可压时遇到的问题，提出处理过盈配合中橡胶计算的方法。
2、可压模型理论解与ABAQUS数值解的比较
2.1、理论解
理论解模型如图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-2 | 4.6% | 径向应力S11 | 59.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.30 | 0.57% | 80.0（钢） | 0 | 28.15 | － | 环向应力S22 | 59.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.36 | 5541.00 | -1.52% | 80.0（钢） | 4956.85 | 4957.00 | 0.00% |
3、可压缩模型
橡胶的应变能采用多项式模型时，在静水压力荷载下p与J的关系如下：（图片） 用ABAQUS对这1-4组系数进行评估：（图片） （图片）
图3 不同系数对应的橡胶静水压力下的应力应变关系 将这六种橡胶本构代入第二部分中的算例中进行计算结果如下：（图片）
图4 第6组系数对应的位移图图5 第1组系数对应的位移图 由图4和图5容易看到这两组系数对应的位移差异非常大。表3. 不同系数对过盈配合的影响
- | 橡胶外侧位移U1 | 外部钢环的内侧位移U1（mm） | 过盈面上的应力S11（MPa） | 橡胶内最大S22（MPa） | 橡胶内最大S33（MPa） | 1 | -1.866 | 3.36E-02 | -12.175 | -11.94 | -11.9 | 2 | -1.60 | 2.95E-02 | -106.2 | -106.3 | -106.0 | 3 | -0.7 | 1.21 | -422.5 | -422.5 | -422.5 | 4 | -0.17 | 1.735 | -602 | -602 | -602 | 5 | -0.11 | 1.79 | -621.0 | -621.0 | -621.0 | 6 | -9.6E-2 | 1.825 | -631.7 | -632.0 | -631 |
4、结论
由以上分析可以看到橡胶的可压缩性在过盈配合中起到的作用是非常大的。
例如，使用初始Poisson's Ratio为0.499时，橡胶相对于钢是非常软的，因为这组系数确定的橡胶的本构模型杨氏模量大约为3MPa，体积模量（ABAQUS对系数的评估给出的结果），而钢的杨氏模量为210000MPa，取时。
因此，如果用实验方法来确定橡胶应变能的形式，应对橡胶进行三向应力实验，以确定橡胶应变能中与体积应变有关的各个系数。
参考文献[1] ABAQUS/Thoery用户手册V6.1 HKS.INC 2000
3/14/2005
|
