随着汽车工业的持续进步，轮胎作为车辆不可或缺的一部分，其性能标准也愈发严苛。在众多性能指标中，轮胎的疲劳破坏特性尤为引人瞩目，因为它直接关系到轮胎的可靠性和使用寿命。而整胎的疲劳破坏往往以橡胶基体的损坏为起点，这意味着橡胶基体的性能稳定性对轮胎的整体寿命具有决定性的影响。近年来，随着有限元技术的飞速发展，利用有限单元法来研究橡胶材料的疲劳破坏特性已成为学术界和工业界的共同选择，这种方法为轮胎的设计和性能优化提供了有力的技术支撑。
本论文依托于abaqus有限元软件平台，精心构建了轮胎橡胶的二维模拟模型。通过该模型，深入分析了在单向拉伸条件下，轮胎橡胶中裂纹的扩展方向。为了验证模拟结果的准确性，将所得计算结果与实际试验数据进行对比，发现两者高度一致，充分证明了本论文所建立的模型和分析方法的可靠性。
其中 i1、i2、i3 为主伸长比的第一、第二、第三不变量。
本论文采用yeoh橡胶材料本构模型，co10、c20、c30通过单轴拉伸的实验数据获得。
1.2 橡胶断裂理论
在线弹性断裂力学中，预测裂纹扩展规律时要考虑裂纹沿什么方向扩展，即确定开裂角。目前被国内外广泛采用的几种判据有最大周向应力理论，能量释放率理论，应变能密度理论。
对于橡胶超弹性材料来说裂纹沿着撕裂能释放率最大的方向扩展，由于应变能释放率是应变能密度的函数，对于给定的裂纹长度下橡胶材料的应变能释放率取决于裂纹尖端未开裂部分各个方向的应变能密度，因此裂纹尖端应变能密度的大小决定着裂纹开裂方向。
2 轮胎橡胶有限元模型的建立
橡胶模型试样采用二维板状结构，板的长度是宽度的2倍。橡胶内部裂纹采用经典的 griffith 裂纹来模拟，即在大的弹性体中间预置一个小的穿透裂纹，为了便于计算取裂纹的长度为板宽的一半。
图1 中心裂纹橡胶模型单向拉伸时的网格划分(a)一般网格(b)细化后的网格
由于是对称结构，可取四分之一为研究对象，用abaqus有限元软件建立单向拉伸有限元模型。考虑到橡胶材料的高度非线性，通过施加位移来代替均布应力荷载。采用粗细不同的两种网格，单元类型采用平面应变八节点杂交单元，裂纹通过约束对称边界上相应节点位移来模拟。在断裂力学中，对于线弹性小应变构件，裂纹尖端要用断裂奇异单元来模拟。由于橡胶是超弹材料，加载后裂尖会发生很大的变形，考虑到有限元计算的收敛性，裂纹尖端可以不用断裂奇异单元。图1为单向拉伸时的有限元模型。
由于是对称结构，可取四分之一为研究对象，用 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
