基于rebar单元的载重子午线轮胎模型建立及验证

1、基于rebar单元的载重子午线轮胎模型建立及验证作者：应世洲陈方王国林来源：PT现代塑料随着计算机技术的飞速发展，以有限元法为核心的计算机辅助工程软件(CAE软件)在轮胎设计和分析中得到迅速应用世界各大轮胎公司普遍利用CAE软件在计算机上建立数字轮胎模型，模拟轮胎静态、动态过程及轮胎的实际使用伏况，从而在设计阶段就获得轮胎的各种静力学及动力学特性。轮胎是由橡胶和骨架材料组成的一种十分复杂的柔性结构轮胎模型建立得是否准确对进一步的分析起着至关重要的影响，因此建模的方法很重要本文采用有限元商业软件ABAQUS基于rebar单元建立载重子午线轮胎12.00R20的有限元模型通过对比分析试验测得的载荷

2、一下沉量曲线和仿真模拟得到的曲线，进一步确认该建模方法的可行性。1有限元模型的建立1.1模型简化为了提高计算效率，对模型作了如下简化:(1)忽略防擦线、标志线等。(2)忽略横向花纹的影响;(3)忽略0度带束层与第3带束层间的间隙(见图1);图1带束层位置示意(4)由于气密层很薄，如果也对其进行网格划分，则网格的质量较差，因此将该层与内衬层合并为一层划分网格;(5)考虑到胎肩垫胶和胎冠基部胶用的是同一种胶料，因此将这两处合并进行网格划分。1.2帘线一橡胶复合材料模型建立为了分析帘线一橡胶复合材料的力学性能，本研究采用在橡胶实体单元中定义rebar单元的方法，其中实体单元模拟橡胶的力学性能,reb

3、ar单元模拟钢丝帘线的力学性能，这样大大方便了模型的建立。在ABAQUS中有两种方法定义rebar单元，一种是直接将rebar单元定义在橡胶实休单元内，如图2所示，图中橡胶基体单元的节点编号依次为1,2,3和4(图中实线表示轴对称单元的边，虚线是为了表示帘线铺设方向而画的辅助线)，由节点编号顺序所决定的单元边的编号为edgel,edge2,edge3和edge4,帘线铺层和橡胶基体单元的交线位置是可以定义的，该交线还规定了正方向以便定义rebar的方位，因为交线总是要与基体单元的两条边相交，根据交点所在的单元边的编号定义交线的正方向，即边的编号由小到大的方向为正方向。teMt方粧甬蛊考正为甸图

4、3局部坐标系的建立另一种方法是先将rebar单元定义在面单元上，再将面单元嵌人到相应的橡胶实体单元内，如图3所示。图中1方向是整体坐标轴x轴向该面单元投影得到的，图中所示的局部坐标满足右手规则。在模型中还可以定义钢丝帘线的帘线角、帘线间距、单根钢丝的截面积及钢丝截面中心距面单元中面的距离，如图4所示，帘线角度为1方向与实际钢丝帘线铺设位置的夹角，即所要定义的rebar角。周向世bif甬默诉或用户宦丈的齒的揑影力离而单无中葩nwh聽中春拒蚯单元申閒的跑离图4rebar单元的有关参数定义经过比较分析发现，运用第2种方法最大的优点在于可以减少橡胶实体单元网格划分的困难，因为在橡胶实体单元网格划分中帘

5、布层的位置由面单元来准确模拟。其次，可以分别定义橡胶基体单元和钢丝帘线的材料参数，因此本研究均采用第2种方法模拟钢丝帘线，材料参数由试验得到。1.3橡胶本构模型选择目前，橡胶的本构模型有很多种，大致分为两类，一类是唯像理论模型，另一类是分子网络理论模型。这些模型已经被一些大型商业软件(如ABAQUS和ANSYS等采用，但是在选用橡胶材料本构模型时存在下述间题。有些模型在拟合小变形范围的试验数据时效果较好，而另一些模型在拟合大变形时效果较好，从而说明不同的本构模型对橡胶材料力学特性的表征效果是不同的。利用同一材料模型、选取不同变形范围的试验数据进行拟合得到的材料参数也有明显的区别，表明本构模型的

6、表征效果还与变形范围有关。轮胎中橡胶的应变一般不大于50%，但是不少研究者对橡胶大变形范围内的力学特性也很重视，因此，本研究首先对材料做单轴拉伸试验，运用ABAQUS软件对试验得到的材料参数进行评枯，在综合考虑不同变形范围内的曲线拟合效果后决定采用Yeoh模型。1.4单元选择用二维模型分析轮胎装配和充气工况，几何条件和载荷条件是轴对称的，而且帘线的铺设形式也是轴对称的，因此可以用平面轴对称有限元模型来模拟。由于带束层和胎体层端部帘线的影响，导致轮胎绕旋转轴扭转变形，因此针对四边形单元和三角形单元分别采用CGAX4H和CGAX3H轴对称单元来模拟。1.5接触模拟本研究采用直接约束法处理接触间题。

7、直接约束法是追踪物体的运动轨迹，一旦探测出发生接触，便将接触所需的运动约束和节点作为边界条件直接施加在产生接触的节点上。这种方法对接触的描述精度高，具有普遍适应性.模拟时将轮辆和路面定义为解析刚体。采用解析刚体能更好地模拟刚体形状，从而使模拟出来的刚体表面更光滑，降低接触噪声。2模型分析及验证2.1轮胎有限元分析模型本模型以12,OOR20载重子午线轮胎为原型，使用8.51型标准平底轮辆，额定充气压力为840KPa，单胎额定载荷为36554N,所建立的二维和三维模型如图5和6所示，二维模型单元数为501个，节点数为560个，三维模型沿圆周方向划分100等份，单元总数为50100个。图5二维轮胎

8、有限元模型图6三维轮胎有限元模型2.2模型分析模型验证的方式有很多种，目前主要有测试轮胎与路面的接触压力和面积及轮胎充气断面宽等方法，本研究采用测试轮胎载荷一下沉量曲线的方法，这也是目前比较常用的验证方法之一。首先把轮胎装上轮辆并按额定气压充气，然后装到载荷下沉量试验机上。试验前先测试轮胎充气外半径，然后测试不同载荷下轮惘中心到接地中心的距离，经过换算得到在不同载荷下的下沉量。有限元分析和试验结果对比如图7所示，误差最大为5.9500，最小为3.99符合工程研究的要求，可以认为所建立的模型是准确的。爲L.020-有鶴元祐i*一试爺俏亠蹄304050TMS/niin2.3轮胎变形用二维模型分析轮胎装配和充气工况，轮胎在充气压力的作用下产生轴对称变形，因此可以用轮胎断面来表示轮胎整体的变形情况，如图8所示.由图8可以看出，胎圈处的变形较大，主要是由装配引起的，并在胎趾处离轮辆翘起一定距离;在胎冠处沿着1方向产生一定的径向位移，该位移主要是由充气的。图8二维模型轮胎变形轮胎在接地区所受的载荷是非轴对称的，载荷影响区域也主要集中在接地区，如图9所示。由图9可以看出，在额定静负荷下，轮胎与路面接触区域为离接地面中心两侧各150范围由于接地影响而产生明显非轴对称变形的范围大约为离接地面中心两侧各30度。图9三维模型轮胎变形3结论经过对比分析由试验和仿真模拟得到的载荷一下沉量曲线发现，