一种计算轮胎胎面压缩对滚动阻力贡献的方法、设备和计算机程序技术

本发明专利技术涉及轮胎计算机智能设计技术领域，尤其涉及一种计算轮胎胎面压缩对滚动阻力贡献的方法、设备和计算机程序。本发明专利技术通过计算过程即可获得计算轮胎胎面部位由于压缩变形产生的滚动阻力以及由于弯曲变形产生的滚动阻力。实现对胎面部位在不同变形模式下的滚动阻力分析，为降低轮胎滚动阻力提供分析手段。提供分析手段。提供分析手段。

【技术实现步骤摘要】
一种计算轮胎胎面压缩对滚动阻力贡献的方法、设备和计算机程序


[0001]本专利技术涉及轮胎计算机智能设计
，尤其涉及一种计算轮胎胎面压缩对滚动阻力贡献的方法、设备和计算机程序。

技术介绍

[0002]轮胎作为汽车的必备组成部件，其在服役过程中由于本身橡胶材料的粘弹性质会引起滚动能量损失，即滚动阻力，滚动阻力将消耗一定数量的电能，降低电动汽车的行驶里程，低滚动阻力轮胎早已成为轮胎行业的共识性追求。影响轮胎滚动阻力的因素主要有材料、结构、路面条件和使用工况等，如何调整结构设计以使轮胎的滚动阻力最小是各个轮胎企业的关键核心技术。轮胎滚动时胎面的能量损失占整个轮胎滚动阻力的40％～80％，是轮胎所有部件中占比最多的部件，通过调控胎面材料内部的能量分配能显著降低轮胎的滚动阻力。胎面材料的变形主要分为压缩变形和弯曲变形，压缩变形是由于轮胎与路面接触产生的压力导致，弯曲变形是由于轮胎周向圆形环与路面接触后的曲率变化导致。对于轮胎设计工程师而言，在设计和分析轮胎时需要获得两种变形各引起轮胎滚动阻力的数值大小，以确定和调整设计方向。然而，目前测试方法只能获得整个轮胎的滚动阻力值，无法获得具体的胎面压缩和弯曲所占数值(或比例)，而仿真方法则较多关注轮胎整体和各部件的滚动阻力，而没有区分胎面不同变形模式对滚动阻力的影响，导致轮胎设计工程师在改进轮胎结构时缺乏指导依据。

技术实现思路

[0003]本专利技术的目的是为了解决上述现有技术存在的问题，即对于轮胎设计工程师而言，在设计和分析轮胎时需要获得两种变形各引起轮胎滚动阻力的数值大小，以确定和调整设计方向。然而，目前测试方法只能获得整个轮胎的滚动阻力值，无法获得具体的胎面压缩和弯曲所占数值(或比例)，而仿真方法则较多关注轮胎整体和各部件的滚动阻力，而没有区分胎面不同变形模式对滚动阻力的影响，导致轮胎设计工程师在改进轮胎结构时缺乏指导依据。
[0004]进而提供一种计算子午线轮胎胎面压缩对滚动阻力贡献的方法，实现对胎面部位在不同变形模式下的滚动阻力分析，为降低轮胎滚动阻力提供分析手段。
[0005]为了实现上述的目的，本专利技术采用了以下的技术方案：
[0006]一种计算轮胎胎面压缩对滚动阻力贡献的方法，该方法包括以下的步骤：
[0007]1)建立轮胎有限元模型，进行二维轴对称充气分析：对轮胎的设计图纸进行网格划分，赋予材料属性，在轮胎内部边界层施加额定的充气压力，进行充气分析；
[0008]2)在充气分析的基础上，使二维轴对称模型旋转一周，生成轮胎三维模型，在距离轮胎表面1mm的位置放置刚性平直路面，固定轮胎轮辋，对刚性路面施加额定载荷，进行轮胎负载分析，输出胎面部位的应力值和应变值；
[0009]3)使轮胎的带束层、胎体及冠带层含有骨架材料的部件保持步骤2)中变形后的位置不变，卸载路面受力，即使刚性路面施加值为0的载荷，输出计算后胎面部位的应力应变值；
[0010]4)利用程序分别提取步骤2)和步骤3)中的轮胎一周处于相同位置不同角度的单元真实应力和应变值分别记作σ值和ε值，以轮胎周向角度为x值，x以弧度为单位，利用400阶三角级数对σ，ε值进行拟合：
[0011][0012][0013]5)利用拟合后的参数数据分别第二步和第三步中胎面橡胶材料单元的在不同方向应力下的能量损耗密度e，tanδ值为橡胶材料的损耗正切，
[0014][0015]将6个应力下的能量损耗密度相加，再乘以单元的体积V
el
，即为单元m的能量损耗W
m
：
[0016][0017]将胎面所有橡胶材料单元的能量损耗相加即为整个胎面部件的能量损耗，再将胎面部件的能量损耗除以轮胎的周长即为胎面部件的滚动阻力R
T
：
[0018][0019]6)根据步骤5)的计算方法，计算获得步骤2)中胎面的滚动阻力为R
T2
，步骤3)中胎面的滚动阻力为R
T3
，R
T3
即为弯曲变形模式下胎面的滚动阻力，R
T2
‑
R
T3
即压缩变形模式下胎面的滚动阻力。
[0020]作为优选，步骤1)中对轮胎材料分布图进行网格划分，划分为三角形或四边形单元，骨架材料划分为2节点一维单元，并对各部件材料赋予材料属性建立轮胎有限元模型。
[0021]作为优选，步骤1)中给轮胎内衬层内表面施加0.25MPa气压，利用Abaqus软件进行充气分析。
[0022]作为优选，步骤2)中在步骤1)二维充气分析的基础上，将轮胎断面旋转360度，并在周向进行断面划分，将一周分为110个断面，将轮胎轮辋固定，给路面施加额定载荷，使路面向轮辋方向移动，进行负载分析，计算完成后，输出所有橡胶材料单元的应力应变值。
[0023]作为优选，步骤3)中位置不变为限制包括胎体、带束层、冠带层在内的骨架材料节点的6个方向自由度，使其不能变形。
[0024]作为优选，步骤4)中利用python程序分别提取步骤2)和步骤3)中中的轮胎一周处
于相同位置不同角度的单元真实应力S11，S22，S33，S12，S13，S23和应变值LE11、LE 22、LE 33、LE 12、LE 13、LE 23，分别记作σ值和ε。
[0025]进一步，本专利技术还公开了所述方法用于降低轮胎滚动阻力中的应用。
[0026]进一步，本专利技术还公开了一种计算机设备，包括存储器、处理器及存储在存储器上的计算机程序，所述处理器执行所述计算机程序以实现上述任意一个技术方案所述方法的步骤。
[0027]进一步，本专利技术还公开了一种计算机可读存储介质，其上存储有计算机程序或指令，该计算机程序或指令被处理器执行时实现上述任意一个技术方案所述方法的步骤。
[0028]进一步，本专利技术还公开了一种计算机程序产品，包括计算机程序或指令，该计算机程序或指令被处理器执行时实现上述任意一个技术方案所述方法的步骤。
[0029]本专利技术由于采用了上述的技术方案，通过上述计算过程即可获得计算轮胎胎面部位由于压缩变形产生的滚动阻力R
T2
‑
R
T3
以及由于弯曲变形产生的滚动阻力R
T3
。实现对胎面部位在不同变形模式下的滚动阻力分析，为降低轮胎滚动阻力提供分析手段。
附图说明
[0030]图1为21550R15轮胎材料分布图；
[0031]图2为21550R15轮胎断面网格及材料部件；
[0032]图3为21550R15轮胎充气变形结果；
[0033]图4为21550R15轮胎三维周向网格划分结果；
[0034]图5为21550R15轮胎三维负载结果及方向标识；
[0035]图6为输出第二步中的轮胎部件单元的部分应力值；
[0036]图7为输出第二步中的轮胎部件单元的部分应变值；
[0037]图8为第三步计算完成载荷释放后胎面变本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种计算轮胎胎面压缩对滚动阻力贡献的方法，其特征在于，该方法包括以下的步骤：1)建立轮胎有限元模型，进行二维轴对称充气分析：对轮胎的设计图纸进行网格划分，赋予材料属性，在轮胎内部边界层施加额定的充气压力，进行充气分析；2)在充气分析的基础上，使二维轴对称模型旋转一周，生成轮胎三维模型，在距离轮胎表面1mm的位置放置刚性平直路面，固定轮胎轮辋，对刚性路面施加额定载荷，进行轮胎负载分析，输出胎面部位的应力值和应变值；3)使轮胎的带束层、胎体及冠带层含有骨架材料的部件保持步骤2)中变形后的位置不变，卸载路面受力，即使刚性路面施加值为0的载荷，输出计算后胎面部位的应力应变值；4)利用程序分别提取步骤2)和步骤3)中的轮胎一周处于相同位置不同角度的单元真实应力和应变值分别记作ζ值和ε值，以轮胎周向角度为x值，x以弧度为单位，利用400阶三角级数对ζ，ε值进行拟合：角级数对ζ，ε值进行拟合：5)利用拟合后的参数数据分别第二步和第三步中胎面橡胶材料单元的在不同方向应力下的能量损耗密度e，tanδ值为橡胶材料的损耗正切，将6个应力下的能量损耗密度相加，再乘以单元的体积V
el
，即为单元m的能量损耗W
m
：将胎面所有橡胶材料单元的能量损耗相加即为整个胎面部件的能量损耗，再将胎面部件的能量损耗除以轮胎的周长即为胎面部件的滚动阻力R
T
：6)根据步骤5)的计算方法，计算获得步骤2)中胎面的滚动阻力为R
T2
，步骤3)中胎面的滚动阻力为R
T3
，R
T3
即为弯曲变形模式下胎面的滚动阻力，R
T2
‑
R
T3
即压缩变形模式下胎面的滚动阻力。2.根据权利要求1所述的一种计算轮胎胎面压缩对滚动阻力贡献的方法，其特征在于，步骤1)中对轮胎材料分布图进行网格划分，划分为三角形或四边形单元，骨架...

【专利技术属性】
技术研发人员：王剑波，刘芳，胡德斌，崔志博，郭磊磊，李进，王毅，
申请(专利权)人：中策橡胶集团股份有限公司，
类型：发明
国别省市：