滚动橡胶轮胎模态仿真到底怎么做？

详细信息

作者优秀: 优秀教师/意见领袖/博士学历/特邀专家/独家讲师
平台推荐: 内容稀缺

1年前浏览5112

导读：随着企业对研发技术的愈来愈重视，对提高设计仿真研发效率、实现知识工程化的需求愈加旺盛。虚拟仿真是轮胎设计研发的重要内容，是增强企业产品竞争力的关键之一。达索系统SIMULIA提供针对轮胎智能化仿真的完整解决方案，包括基于Abaqus的非线性分析，基于二次开发的快速建模和自动网格，以及结合Abaqus和SIMULIA各个工具进行动力学、噪声等分析，为轮胎智能化仿真提供了强大技术支撑。3月9日14时，由达索系统主办《探索轮胎智能化仿真的关键技术》将在仿真秀官网和APP同步直播，笔者强烈推荐工程师朋友参加，详情见后文。

一、写在文前

滚动轮胎模态仿真实际上是在轮荷加载的基础之上的重启动分析。轮胎在稳态滚动过程中,会受到预加载荷、惯性力以及轮胎和地面的摩擦力的影响,这些力会对整个系统的刚度矩阵和阻尼矩阵产生影响,导致非对称性。故不能采用常规方法对动力学方程进行解耦,必须用复模态来解耦，所以滚动轮胎的模态仿真其实是复模态的的提取。

在abaqus的稳态滚动中,轮胎实际上是不滚动的,只是内部材料的流动（欧拉-拉格朗日法）。在提取复模态之前,必须保证轮胎滚动的转速和线速度相匹配，故需先进行roll tire仿真调试：

roll tire计算中，先给定轮胎线速度、转动角速度然后提取轮胎轮心的力矩M,当轮胎稳态滚动的时候, 轮胎轮心的力矩M应该为0。在实际操作中,需要不断的调节定义的ω值,使最终绕Y向的力矩M在[-10,10]之内。

进行roll tire计算时，首先进行step1二维轮胎充气仿真计算，然后进行step2rev旋转3D轮胎生成及轮荷加载计算，在此基础上进行Free roll计算Inp文件的编写，进行计算，查看轮胎轮心的力矩M判断轮胎是否处于稳态滚动状态。下图为step1.inp以及step2rev.inp运行结果图：

二、roll tire计算

step3roll_tire.inp文件编写如下:

*HEADING（重启动分析） Step3: roll_tire from step2rev.inp*RESTART,READ*STEP,INC=500,NLGEOM=YES,UNSYMM=YES（几何非线性打开；非对称打开）4: roll_tire at 60km/h*STEADY STATE TRANSPORT,LONG TERM,INERTIA=YES（动态分析，不考虑粘塑性影响，惯性打开） 0.5, 1.0, 1E-6, 1.0*TRANSPORT VELOCITY（初始角速度=V/r） TIRE-RIM,57.0776*MOTION,TYPE=VELOCITY,TRANSLATION（线速度V） NODE_TIRE,1, （方向为1）,16666.67（60公里每小时化为毫米/S）*CHANGE FRICTION,INTERACTION=SRIGID（改变摩擦系数）*FRICTION,SLIP=0.02 1.0*NODE PRINT,FREQ=1*EL PRINT,FREQ=1*OUTPUT,HISTORY,FREQ=1*NODE PRINT,FREQUENCY=999,NSET=TIRE-RIM RF, （输出轮心的力矩）*END STEP

将上述所编写的inp文件命名为step3roll_tire.inp并保存至step2rev.inp运算结果所在文件夹，在文件夹空白处单击右键选择“在此处打开命令窗口”，输入:

  abaqus job=step3roll_tire  oldjob=step2rev int cpus=4

点击回车运行

计算完成后可出现以下几个文件：

其中：

提交计算文件的语法错误信息在.dat文件中查看；

提交计算文件的计算步长在.sta文件中查看；

所得计算结果的.odb文件可在Hyper view或Abaqus中查看；

因我们需要的是轮胎轮心的力矩M：

故用UE打开step3roll_tire.dat文件查看NODE OUTPUT：

因我们需要的是绕Y向的力矩M（RM2）在[-10,10]之内，这时才算稳态滚动，而输出值RM2为-9.5877E 05故需调整ω值,使最终绕Y向的力矩M在[-10,10]之内。在此大家自己调试即可。

三、滚动轮胎模态提取

在ABAQUS/STANDARD中,复模态的提取采用子空间投影法,而且在复模态提取之前,一般先进行具有对称刚度矩阵的无阻尼系统的特征值和特征模态计算。

滚动轮胎模态提取inp文件编写如下：

  *HEADING Step4:Caculation complex mode from step3roll_tire.inp*RESTART,READ*STEP,PERTURBATION5:static modes analysis*FREQUENCY,EIGENSOLVER=LANCZOS,NORMALIZATION=DISPLACEMENT（实模态提取）,0.1,300.,,,*ACOUSTIC FLOW VELOCITY,ROTATION（定义流体声腔的角速度） NODE_AIR,57.0776*BOUNDARY,OP=NEW ROAD,1,6（路面全约束）*END STEP***STEP,NLGEOM=YES 6:Complex modes analysis,spindle free*COMPLEX FREQUENCY,FRICTION DAMPING=YES,PROPERTY EVALUATION,NORMALIZATION=DISPLACEMENT（复模态提取） ,*ACOUSTIC FLOW VELOCITY,ROTATION NODE_AIR,57.0776*END STEP

将上述所编写的inp文件命名为step4rolltire_mode.inp并保存至step3roll_tire.inp运算结果所在文件夹，在文件夹空白处单击右键选择“在此处打开命令窗口”，输入：

  abaqus job=step4rolltire_mode oldjob=step3roll_tire int cpus=4

点击回车运行

计算完成后可出现以下几个文件：

其中：

提交计算文件的语法错误信息在.dat文件中查看；

提交计算文件的计算步长在.sta文件中查看；

所得计算结果的.odb文件可在Hyper view或Abaqus中查看；

现用UE打开step4rolltire_mode.dat文件可得滚动轮胎的实模态频率如下：

复模态频率如下：

（注：step4rolltire_mode.inp是基于上文step3roll_tire.inp进行运算的，没有进行ω值的调试，故上复模态频率不作为参考依据，下只呈现实模态振型云图，大家可后续进行ω值的,使绕Y向的力矩M在[-10,10]之内再进行复模态频率的提取）

在Abaqus中打开step4rolltire_mode的ODB文件可得：

以上就是本次文章的所有内容，其中Hyper mesh输出的2D.inp文件以及以前课程编写的step1.inp文件、step2rev.inp文件以及本节课编写的step3roll_tire.inp文件、step4rolltire_mode.inp文件我都会以附件的形式给大家，有什么问题点击文尾阅读原文进入笔者专栏文章下方留言，我会酌情回复。

文章末尾点赞和在看，且分享本文到朋友圈前30名用户，将获得本文模型和资料，送完即止。

四、探索轮胎智能化仿真的关键技术

3月9日14时，由达索系统主办“探索轮胎智能化仿真的关键技术”报告会，将针对“轮胎”这一特定研究对象，将介绍将Python语言与Abaqus软件结合开发的轮胎网格自动化划分工具，首次在该行业引入机器视觉和深度学习算法进行快速建模和分析预测，逐步帮助企业实现轮胎仿真工作从传统手工建模、划分网格，向设计仿真一体化和智能化设计方向发展。

以下是直播安排