轮胎有限元模型.docxVIP

轮 胎 有 限 元 模 型

I,＝矿＋戊＋戊 (3-4)

4＝人力＋勾戊＋戊矿 (3-5)

八＝界芬芬 (3-6)

由千认为橡胶材料在变形过程中是各向同性而且不可压缩的，所以

沉＝人人人＝1。由公式(3-1)和(3-3)-(3-6)可以得到

, aw 1aw

T=2（勾 一 －一 一 ）－p, i=l,2,3

祝 矿扒

(3-7)

常用的描述橡胶材料力学性能的模型土要有两类：一类是认为橡胶为连续介原的唯象学方法．另一类是基千统计热力学理论。

唯象学的描述方法假设在未变形状态下橡胶为各向同性材料，即长分子链

方向在橡胶中是随机分布的。这种各向同性的假设是用单位体积（弹性）应变能密度来描述橡胶特性的，应变能密度函数可以表示成土伸长率或变形张鱼的三个不变登的函数。

基千统计热力学方法认为弹性体是很多任意取向的长的柔性分子链通过分子间稀疏的交联点组成的分子网格，由千分子间相互作用力很弱，观察到橡胶中的弹性恢复力土要来自橡胶中摘的减少，墒的减少是由千橡胶的伸长使得橡股结构由高度的无序变得有序巴l。由对橡胶中分子链的长度，方向以及结构的统计得到橡胶的本构关系。

常用的本构模型土要有以下儿种：

(I) 高斯统计模型

假设体系热力学能与分子各链构象无关，两个交联点之间的网链为高斯链，其末端距符合高斯分布，交联点固儿在各自的平均位置上，当橡胶变形时，交联点以相同千宏观变形的比例变形。最后，从分子构象摘的改变得到弹性应变能函数为：

I

W=-flk,7（术＋戊＋戊－3)

2

(3-8)

式中，／1为平均单位体积内的网链数，k为Boltzmann常数，T为绝对温度，人为土伸长比。

高斯统计模型是基千末端距远小千分子链的全部伸长度假设，不能够用来描述分子链的大变形伸长过程，只能近似的预测小变形时的情况。

(2) 多项式形式本构模型

在各向同性假设下，橡胶材料应变能密度W可以分解为应变偏昼能和体

积应变能两部分，形式如下：

W=/(f,-3,7,-3)+g(.J－I)

式中，1 72分别 为第一、第一应变不变昼偏里，J为体积压缩比。

令g=

令g=I:I

一(J-l)？，并进行泰勒展开，可得下式

,=,D,.

(3-9)

W=I; 已．7（,－3)＇（72-3／)

＋工一(J-1);2

(3-10)

九“ ,=,D,．

这种形 式为多项式表示的橡胶本构模型，参数N为多项式阶数。对千多项式模型，材料的初始剪切模昼凡、初始体积模昼k。，都取决千多项式的一阶系

数：

μ,,=2(c0,+c0,), k。=21D, (3-11)

对千完全不可压缩材料有D;=O。

假设橡胶材料在未变形状态下其弹性是各向同性的，且完全不可压缩，则完全多项式形式可简化为：

｀W= 2 C,．1（,-3/(/2-3/ (3-12)

｀

r-0．产°

这是比较常用的橡胶类固体的应力应变关系的处理方法，Rivlin(1951)考虑以上假设从纯数学的角度提出了此模型。

(3) Mooney-Rivlin 材料模型45必l

当 不 可压 缩材 料的 多 项式 形式的 本构 模 型中 只取 i=I J=O和i=O

J=I两项时，就可以得到应变不变呈表示的Mooney-Rivlin模型：

w=c0,(1,-3)+c0,(l,-3) (3-13)

这个模型用起来方 便，但是对千轮胎的变形来说是复杂的，这个材料模型不能预测多轴向数据，只能用千相同变形下的材料变形模拟。

(4) Neo-Hookean材料模型

当在Rivlin形式的模型中只取i=l,J=O或在Mooney模型取前一项时，就可得到最简单的不可压缩橡胶本构模型一 Neo-Hookean材料模型旧l:

W=C0,(/,-3) (3-14)

若令c,.=fl灯／2，则与高斯网络统计理论一致。它是最常用的分子统计

学模型，它只有一个常数，尤其是这种材料模型的无条件稳定性，大可以把一种变形条件下的应力应变曲线拟合的材料常数来预测其它变形方式的应力应变曲