随着汽车工业的发展，轮胎作为其中的一个重要组成部分，对其性能要求越来越高。在轮胎的各项性能中，其疲劳破坏特性越来越受到人们的关注。整胎疲劳破坏是由橡胶基体的破坏为基础的，因此橡胶基体的破坏对整个轮胎的寿命起着决定性的因素。近年来随着有限元技术的不断成熟，用有限单元法来分析橡胶材料的疲劳破坏被各国学者广泛采用。
本论文以 abaqus有限元软件为平台，建立了轮胎橡胶的二维模型，分析了在单向拉伸时裂纹的扩展方向。并将计算结果与试验结果对比，两者吻合的很好。
1橡胶裂纹扩展的理论基础
1.1橡胶材料的本构关系
橡胶是一种非线性的、不可压缩的或近似不可压缩的超弹性材料，其本构关系通常是从应变不变量或基本伸长率表示的应变能密度得到的。rivlin[1提出通用应变能函数
其中 i1、i2、i3 为主伸长比的第一、第二、第三不变量。
本论文采用 yeoh 橡胶材料本构模型，c10。通过单轴拉伸的实验数据获得。
1.2 橡胶断裂理论
由 rilvin 和 thomas 开始，继而有 lake，gent 等人将断裂力学中的 griffith 断裂理论应用于橡胶材料。rilvin 和 thomas 通过大量的实验证明了撕裂能可以有效地解决橡胶大变形试件的断裂问题。断裂力学中的应变能释放率 g，也就是撕裂能，是目前被一致公认的解决橡胶疲劳破坏现象最有效的参数。gjlake 系统地研究了橡胶材料的疲劳断裂问题,提出了能量法,定义裂纹每扩展单位面积所吸收的能量为应变能释放率。
其中 u为储存在试件中的总弹性应变能,a 为裂纹的一个断裂面积。他指出裂纹尖端的应变能释放率 g 的具体表达式为.
其中 k为材料常数，wt为裂纹尖端的应变能密度，为裂纹尖端的有效直径。在线弹性断裂力学中，预测裂纹扩展规律时要考虑裂纹沿什么方向扩展，即确定开裂角。目前被国内外广泛采用的几种判据有最大周向应力理论，能量释放率理论，应变能密度理论。对于橡胶超弹性材料来说裂纹沿着撕裂能释放率最大的方向扩展，由于应变能释放率是应变能密度的函数，对于给定的裂纹长度下橡胶材料的应变能释放率取决于裂纹尖端未开裂部分各个方向的应变能密度，因此裂纹尖端应变能密度的大小决定着裂纹开裂方向。
2 轮胎橡胶有限元模型的建立
橡胶模型试样采用二维板状结构，板的长度是宽度的2倍。橡胶内部裂纹采用经典的 griffith 裂纹来模拟，即在大的弹性体中间预置一个小的穿透裂纹，为了便于计算取裂纹的长度为板宽的一半。
图1 中心裂纹橡胶模型单向拉伸时的网格划分(a)一般网格(b)细化后的网格
由于是对称结构，可取四分之一为研究对象，用 abaqus有限元软件建立单向拉伸有限元模型。考虑到橡胶材料的高度非线性，通过施加位移来代替均布应力荷载。采用粗细不同的两种网格，单元类型采用平面应变八节点杂交单元，裂纹通过约束对称边界上相应节点位移来模拟。在断裂力学中，对于线弹性小应变构件，裂纹尖端要用断裂奇异单元来模拟。由于橡胶是超弹材料，加载后裂尖会发生很大的变形，考虑到有限元计算的收敛性，裂纹尖端可以不用断裂奇异单元。图1为单向拉伸时的有限元模型。
3 有限元计算结果及分析
第一次加载使裂纹尖端单元应力小于橡胶材料单轴拉伸试验时对应的破坏应力;第二次加载使裂尖单元应力超过橡胶单向拉伸试验破坏时对应的应力大小。此时两种不同载荷下对应不同网格的裂纹尖端单元的应变能密度曲线如图2示。
图2单向拉伸时裂尖单元的应变能密度曲线图(a)一般网格(b)细化网格
当所加应力荷载较小时，单元应变能密度最大值出现在与初始裂纹垂直的单元处，裂纹尖端的撕裂能不足以使其开裂。对应第二次加载与第一次不同，裂尖单元的应变能密度最大值出现在与初始裂纹夹角近似165的单元处，应变能释放率沿该方向最大，裂纹在此方向处失稳断裂。对于橡胶材料来说只有当撕裂能超过一定值时，裂纹才会失稳开裂导致试件最终破坏。有限元计算的中心裂纹橡胶模型的裂纹初始扩展力为145n。
图3 裂纹扩展方向示意图
4 橡胶试样的试验分析
4.1试验模型及设备橡胶拉伸试样采用哑铃状其各项尺寸按最新国标 gb/t528-1998 规定的尺寸标定裁减如图4所示。图5是带裂纹的试样模型。试验仪器采用css-44000系列型电子万能拉伸试验机。拉伸试验中用xy函数记录仪记录载荷f和位移 d。考虑到边界效应，一般取预置裂纹长度a和试样宽度w之比为:0.15
