轮胎转动惯量的有限元计算方法

2016-07-30陈振艺

橡胶科技 2016年11期

陈振艺

（双钱集团上海轮胎研究所有限公司，上海 200245）

轮胎是汽车的重要组成部分，轮胎性能与汽车整体性能密切相关，因此越来越多的汽车企业要求轮胎企业尽可能多地提供轮胎性能指标，其中包括与汽车启动、加速和制动性能密切相关的轮胎转动惯量。

目前我国还没有方便、有效测试轮胎转动惯量的设备和方法，需要通过计算获得轮胎的转动惯量。由于轮胎各个部件形状不规则，估算精度不高，而目前广泛使用的有限元分析软件（如Abaqus等）也没有便于计算轮胎转动惯量的功能，所以需要开发专用软件计算轮胎转动惯量。

1 现有计算方法

目前可以检索到的轮胎转动惯量计算方法有如下3种[1-3]。

采用有限元方法，也就是利用轮胎性能有限元分析得到的单元信息，通过编程计算各个单元的转动惯量来计算整个轮胎的转动惯量[1]。但该方法只涉及到了单元体积的计算方法，对如何计算单元质心以及如何处理可变帘线单元密度问题并未提及。

利用UG三维软件强大的实体分析功能来获取轮胎各个部件的体积和转动惯量数据，并利用对UG三维软件的二次开发来自动获取相关数据，得到整个轮胎的转动惯量[2]。该方法只需输入各个部件的质量，不用计算体积、质心和密度等。该方法简单，编程工作量小，缺点是需要对整个轮胎进行三维建模，而不只是对花纹三维建模，这增加了很多工作量，并无法处理可变帘线单元密度问题。

采用将一个不规则形状物体积分解成多个规则形状物体积之和的方法来计算各个部件的转动惯量[3]。该方法的优点是避免了不规则形状部件质心的计算，计算精度高；缺点是计算方法复杂，编程复杂，也无法处理可变帘线单元密度问题，而且只适用于无花纹的光面轮胎。

2 计算方法确定

对比以上3种轮胎转动惯量的计算方法，有限元方法最方便，且编程难度适中，但主要需要解决3个问题：单元转动惯量的计算、帘线单元密度的处理和花纹沟部分的处理。本研究确定采用有限元方法计算轮胎转动惯量。

2.1 单元转动惯量的计算

转动惯量的计算公式很简单，一个质点的转动惯量（J0）计算公式见式（1）。

式中，m0为质点质量，r0为质点到转动轴的距离。如果一个物体的质量是连续分布的，则这个物体的转动惯量（J）计算公式见式（2）。

式中，ρ为物体密度，r为物体到转动轴的距离，V为物体体积。如果一个单元内的材料密度是均匀的，且单元的形状基本上是规则的，那么对一个单元来说，式（2）计算并不复杂。整个轮胎的转动惯量则等于各个单元转动惯量的总和。

对轮胎而言，单元可以分为3种类型，胎面单元、除胎面以外的其他纯橡胶单元和含帘线单元。胎面单元涉及花纹沟，因此无法在空间上形成连续圆环，需要单个计算，计算量较大；而其他纯橡胶单元和含帘线单元可以形成圆环，计算量较小。含帘线单元涉及到密度处理问题，计算步骤比纯橡胶单元多。

对一个矩形的平面单元来说，如果可以形成一个圆环，那么这个圆环的转动惯量（Jc）计算非常简单，根据式（2）可以得到式（3）。

式中，mc为整个圆环的质量，r1和r2分别为圆环的内、外半径，各参数如图1所示。

图1 圆环转动惯量计算参数示意

Jc的近似计算公式如下：

式中，rc为平面单元质心到轮胎旋转轴的距离，即r1和r2的平均值，Vc为圆环体积，ρc为平面单元密度，Sc为平面单元面积。

式（3）和（4）的计算结果差异取决于平面单元高度（h）与rc的比值，比值越大，差异也越大。当rc为200时，如果h为4，则差异为0.01%；如果h为6，差异为0.1%；如果h为10，差异为0.25%。因此当h不太大时，式（4）的计算精度是令人满意的。

对于矩形单元式（4）的优势并不明显，但对于任意四边形或者三角形的平面单元，式（4）就显示出明显的优势。只要求出平面单元的面积和质心位置，就能求出圆环的体积，再根据密度就能计算出转动惯量。如果要获得更高的计算精度，可以把一个单元划分成若干个三角形或者四边形，再分别计算后求和即可。

假定三角形的3个顶点坐标分别为（x1，y1），（x2，y2）和（x3，y3），那么三角形面积（SA）计算公式见式（5）。

对密度均匀的三角形来说，质心位置（x0，y0）也就是重心位置，计算公式见式（6）和（7）。

假定四边形的4个顶点坐标分别为（x1，y1），（x2，y2），（x3，y3）和（x4，y4），用两条对角线将四边形划分为4个三角形，分别求出4个三角形的面积S1，S2，S3和S4，以及4个三角形的质心位置（u1，v1），（u2，v2），（u3，v3）和（u4，v4），那么四边形的面积（SB）和质心位置（x00，y00）计算公式见式（8）—（10）。

因此，对于除胎面以外的其他纯橡胶轮胎部件来说，只要处理一个轮胎断面上的平面单元数据，就能计算出该部件的转动惯量，可以不涉及该部件沿轮胎圆周方向的单元划分信息。

2.2 帘线单元密度的处理

纯橡胶单元的密度是固定的，通常从数据文件中读取即可，而帘线单元的密度则要另行处理。除钢丝圈单元外，其他帘线单元有两种处理方法。第1种方法是将帘线单元看作密度均匀的单元，根据帘线的直径、排列密度和密度，以及帘线夹胶的密度计算出平均密度，再计算转动惯量。第2种方法是沿着帘线的排列方向将帘线单元划分成3层，上下两层为纯橡胶，中间层为帘线和橡胶，厚度与帘线直径相同，然后分别计算转动惯量。无论是哪种处理方法，都需要计算密度。

第1种方法中帘线单元的平均密度是动态的，即使帘线规格和排列密度保持不变，轮胎规格不同，因帘布厚度不同其平均密度也会变化。因此在计算时要额外计算单元厚度，即使同一个轮胎部件的单元厚度也会有所不同，其平均密度也不同。第2种方法中间层的密度是固定的，只要根据帘线的直径、排列密度和密度，以及帘线夹胶的密度计算1次就可以了，但是将1个单元划分成3层，前处理工作量会增加，只是计算精度从理论上说比第1种方法要高。

钢丝圈单元的处理与第1种方法类似，需要根据钢丝直径、钢丝排列方式和夹胶厚度计算出平均密度。

同样，对于含帘线部件来说，只要处理一个轮胎断面上的平面单元数据，就能计算出该部件的转动惯量，可以不涉及该部件沿轮胎圆周方向的单元划分信息。

2.3 花纹沟部分的处理

胎面单元的处理也分为两种情况：一是单元划分时已经包含花纹沟信息，二是不包含花纹沟信息。如果单元划分时已经包含了花纹沟信息，那就涉及到沿轮胎圆周方向的胎面单元划分信息，也就是必须逐个计算胎面单元的转动惯量（Jt），近似计算公式见式（11）。

式中，mt为胎面单元质量，rt为胎面单元到轮胎旋转轴的距离，ρt为胎面单元密度，Vt为胎面单元体积。胎面单元大部分为六面体，还有少量五面体，按轮胎各个单元的通用划分方法，这些六面体和五面体基本上都是规则的。也就是说总是存在两个面，这两个面延长的相交线就是轮胎的旋转轴，并且绝大部分情况下这两个面的形状和尺寸基本相同，接近平行，如图2所示的面1和2。

图2 胎面六面体单元体积计算参数示意

六面体单元的Vt近似计算公式见式（12）。

式中，St1和St2分别为面1和2的面积，H是面1质心（x01，y01）到面2质心（x02，y02）的距离，而六面体单元的rt就等于面1和2的质心到轮胎旋转轴距离的平均值。

五面体单元的体积计算方法与六面体类似。如果面1和2的形状或者尺寸相差较大，可以将单元沿轮胎圆周方向划分成若干个小单元分别计算，以提高Vt和rt的计算精度。

如果单元划分时不包含花纹沟信息，也就是当作光面轮胎来划分单元网格，那么胎面单元的转动惯量计算方法与其他纯橡胶单元的计算方法一样，只需要最后减掉花纹沟部分的转动惯量。由于单元划分时不包含花纹沟信息，因此花纹沟的转动惯量计算是很困难的，通常需要用户自行计算，或者由计算软件根据用户输入的数据进行估算。计算花纹沟转动惯量的方法有多种，计算精度各不相同。例如可以只输入花纹的饱和度、花纹沟的平均宽度和平均深度，就能计算出花纹沟的转动惯量。但由于横向花纹沟和纵向花纹沟的形状、宽度、深度、倾角、长度和数量等有差异，因此这种方法很难精确得到花纹沟平均宽度和平均深度数值，由此得到的花纹沟转动惯量的计算值精度不高。也可以根据花纹沟形状和尺寸不同将花纹沟分为若干类，由用户分别计算这些花纹沟的等效长度，然后将各个类型花纹沟的宽度、深度、等效长度和数量等输入计算软件中，即可得到所有花纹沟的转动惯量，且计算值精度有所提高。

3 计算软件的开发

转动惯量的计算过程并不复杂，编程工作量也不大。开发的有限元方法计算轮胎转动惯量软件应包括界面显示、用户数据输入、数据文件处理、转动惯量计算和计算结果输出等功能，软件的调试和修改为开发的主要内容。

界面设计主要考虑便于参数的输入和输出，以及方便软件各种功能的使用。对于参数输入，应尽量减少输入参数的数量，用数据文件替代参数输入。轮胎有限元分析前后处理会产生许多数据文件，包含大量计算转动惯量所需的数据和信息，充分利用这些数据文件可以减少用户输入参数量。确实需要输入参数时，也应尽量提供参数默认值和选项，也就是说，尽量用鼠标点击替代键盘输入。

在数据文件处理设计时，要注意对数据格式的检查，以及对读取到的数据进行检查，判断其数值是否处在合理范围内，这样有助于判断数据文件是否有误。如果要处理多个数据文件，要注意检查这些文件的匹配性，例如有些参数可能在不同的数据文件中都出现，通过对比这些参数，就可以判断这些数据文件是否匹配，防止使用错误的数据文件。同时在处理数据文件时要考虑效率和技巧，尽量读取一遍就能分离出所有需要的数据，避免反复读取数据文件，节省读取数据的时间。

总的来看，转动惯量的计算方法比较明确。相对而言，将一个立体单元划分成几个小立体单元，并计算出这些小立体单元的顶点坐标的方法比较繁琐；同样当帘线单元的厚度不一致时，计算帘线单元的厚度以及密度的方法也比较繁琐。在计算转动惯量时，还可以顺便计算一些其他参数，例如轮胎的外直径、断面宽、体积和质量等，当单元划分包含花纹沟信息时，可以顺便计算出所有花纹沟的体积，这可能对用户有参考价值。

除了输出整个轮胎的转动惯量外，还可以输出一些辅助信息，例如整个轮胎的体积和质量，各部件的转动惯量、体积和质量等，供用户参考。另外还可以列出从数据文件中获得、或者用户自行输入的一些主要数据，例如一个断面的单元和节点数量、整个轮胎的单元数量、各个部件的密度（特别是通过计算得到的帘线单元平均密度）等，这些参数除了可以提供给用户参考以外，还可以帮助用户判断计算结果是否合理、数据文件或者用户输入的参数值是否正确等。