1. 先说思路轮胎仿真为什么要分三步走用ABAQUS做轮胎仿真很多新手上来就想直接建好轮胎模型加个速度让它滚起来。结果呢要么接触半天不收敛要么应力分布乱成一团根本不知道哪一步出了问题。我一开始也踩过这个坑后来跟着老工程师捋清楚流程才发现轮胎仿真的核心套路就三个字过盈、充气、滚动传荷。这三步不是随便分的每一步都有它的物理意义和数值意义缺一步后面的结果就可能失真。轮胎能不能正常运转首先靠的是胎圈与轮辋之间的过盈配合。没有过盈量轮胎充气后容易滑脱轮辋也带不动胎体。过盈装配在有限元里相当于给接触面预加一个法向压缩量让胎圈“咬”在轮辋上为后续充气和滚动提供初始的接触压力和摩擦基础。你可以把它想象成用力把自行车外胎套进钢圈最后那一下“勒”住的感觉这就是过盈配合。接着是充气。充气压力是轮胎工作状态的来源胎压建立起来后胎体内部产生预张应力整体刚度也随之变大。如果跳过充气直接滚动胎体和帘线层就没有正常的工作应力状态算出来的接地印痕、侧偏特性和耐久性能全是错的。所以在工程仿真中充气是一个独立的静力分析步必须先把胎压加载到设计值保证轮胎“鼓起来”再进行下面的工况分析。最后才是滚动传荷。所谓滚动传荷本质是轮胎在滚动过程中载荷从轮轴经过胎体、帘线层传递到胎面再通过接触压力分布到地面的过程。这一阶段要看的是接地面的压力分布、下沉量、胎肩应力以及不同工况自由滚动、制动、驱动下的力学响应。滚动分析需要在已充气的状态下施加轴荷和运动边界条件让它跑起来。这三步互相依赖过盈为充气提供装配应力充气为滚动提供初始状态。如果你直接滚动相当于让一个瘪胎在地上蹭物理上就不成立。数值上也容易触发接触穿透和网格畸变。所以老老实实按三步走是轮胎仿真最稳妥、也最好排查问题的路径。在求解器选型上我建议过盈和充气阶段使用ABAQUS/Standard的隐式静力学。原因很简单这两个阶段变形相对可控接触状态稳定隐式求解可以不考虑惯性效应收敛性好结果精度也高。滚动阶段则有两种选择一是用Standard中的稳态传输分析Steady-State TransportSST适合做自由滚动和稳态工况效率高二是用Explicit显式分析适合动态工况、瞬态响应以及大滑移或者复杂路面。SST的优势是避免了滚动过程中的惯性振荡计算量小Explicit更能还原真实的瞬态过程但耗时和稳定性就要看建模功力了。工程上耐久、磨损、参数扫描多用SSTNVH和瞬态冲击才上Explicit。2. 关键设置逐个说模型、材料、接触、边界2.1 几何模型2D还是3D怎么选轮胎仿真建模的第一步是几何。工程上最常用的思路是先做2D轴对称模型轮胎断面做过盈和充气求解稳定后把变形后的2D模型旋转拉伸成3D再进行滚动分析。这样做的原因很实际过盈和充气阶段是轴对称工况用2D模型能极大压缩计算量调试参数也快。等旋转成3D后再在周向上加密网格模拟接地和滚动。如果你直接建3D模型做全流程也不是不行但前期调试会非常痛苦。网格稍密一点过盈接触一对一次完整计算可能就要跑几小时甚至十几小时。对称模型跑一个充气只要几分钟等参数都调好了再转3D效率翻倍。这也是为什么绝大多数轮胎仿真教材和工程案例都采用这种方式。几何来源可以是你自己画的轮胎断面图或者逆向扫描的CAD模型。最关键的是保证断面线条封闭、无重叠、无微小的裂缝否则网格生成时会出奇形怪状的单元算到一半就崩。导入ABAQUS之前建议先用CAD软件清理一遍几何把圆角、倒角这些非关键特征简化掉尤其是胎圈区域接触面要尽量光滑连续。另一个细节是单位制。轮胎仿真中我习惯用mm、N、MPa、t吨、s这套单位制。比如密度要写成7.8e-9吨每立方毫米压强用MPa力用N。很多人建模时单位忘了换算结果力气大了十的六次方应力云图直接爆掉。检查单位是我每次建模前必做的一步。2.2 材料定义橡胶和帘线一个都不能糊弄轮胎不是均质体它由橡胶基体和帘线增强层构成。橡胶是不可压缩超弹性材料在ABAQUS里常用*HYPERELASTIC定义本构模型选Yeoh或Ogden。Yeoh模型参数少、拟合稳定性好适合轮胎类的工程计算Ogden模型精度更高但参数标定要求也更高需要完整的单轴、双轴、平面拉伸实验数据。如果你手里只有简单的材料实验报告建议从Yeoh开始。Yeoh模型的参数通常是C10、C20、C30再加一个D参数控制可压缩性。我给一组典型值供参考C100.5C20-0.02C300.01D10.001单位MPa。这些值只是普通胎面橡胶的近似实际项目必须根据你的材料实验数据重新标定。注意超弹性材料在定义时一定要打开*HYPERELASTIC的“全积分”选项否则大变形下容易体积锁定。帘线层用的是嵌入单元方法Embedded Element或者Rebar。后者更专业一点可以在壳单元里定义多层帘线指定帘线间距、角度和材料。在2D轴对称模型里帘线层可以用REBARLAYER定义在橡胶单元内部也可以单独建立膜单元再通过EMBEDDED ELEMENT嵌入到橡胶基体中。用嵌入的方式在3D转换和接触处理上更省事但计算开销略高。帘线材料通常是高强度钢丝或者聚酯定义成线弹性材料即可杨氏模量在MPa量级钢丝可以取200000 MPa聚酯取5000 MPa左右泊松比0.30.4。注意帘线方向是有角度的胎冠区一般为0度圆周方向胎侧和胎圈区域可能有多角度层建模时必须填对角度否则滚动后的受力状态会严重失真。2.3 接触设置过盈量和摩擦系数怎么给过盈配合是轮胎仿真最容易出问题的一环。ABAQUS里模拟过盈有两种常见办法第一种是直接在装配时让胎圈内径比轮辋外径小一点接触对建立后程序会自动识别干涉再通过接触算法逐步消除干涉第二种是使用接触对的“Interference”选项直接指定需要消除的过盈量在分析步里逐渐施加。实践上我更喜欢用第二种因为可控性好。具体操作是建立接触对后在接触属性里选Interference fit输入过盈量数值然后在Step里选择“Ramp over step”或“Instantaneous”。我一般选Ramp让接触压力平滑上升避免一步到位导致单元翻折。轮胎和轮辋之间的过盈量通常取0.51.5 mm具体看实际装配要求别贪大过盈量太大收敛极难。接触属性里除了法向行为切向需要定义摩擦系数。胎圈和轮辋之间是干摩擦夹持摩擦系数一般取0.50.8。这里注意ABAQUS的摩擦默认是“Lagrange”还是“penalty”手写建议使用罚函数法并设置一个合理的弹性滑移容差比如0.005 mm。容差太小接触容易震荡太大摩擦行为就失真了。接触算法里我还建议打开接触面上的“Adjust only master”选项避免从面节点初始穿透的时候被强制拉回导致应力振荡。2.4 分析步和边界条件从过盈一路滚到传荷建议至少建立三个分析步Step-1过盈装配。只施加接触过盈轮辋固定轮胎不施加其他载荷。NLGEOM必须打开因为这个阶段位移量虽然不大但接触和大变形效应明显。Step-2充气。在胎腔表面施加均布压力比如乘用车胎0.25 MPa商用车可能0.81.0 MPa。这个分析步要稳定增量步可以自动。Step-3滚动传荷。施加垂向轴载荷和滚动速度。如果是稳态传输用*STEADY STATE TRANSPORT同时定义转动角速度和前进速度或者让地面移动带动轮胎旋转如果是Explicit需要先建立一个准静态的短时间加载步再进入滚动状态。边界条件上轮辋内表面通常在过盈步开始前就固定所有自由度。充气时胎腔压力方向要朝外用“Surface traction”或“Pressure”均可压力值始终垂直于胎腔表面。滚动时的垂向载荷一般施加在轮轴参考点上并通过*COUPLING将轮轴与胎圈节点耦合起来。滚动速度不是乱给的自由滚动工况需要同时指定角速度和平动速度使接触区的剪切应力满足自由滚动条件驱动和制动工况则通过施加额外的驱动力矩或制动力矩来实现。这里特别提一下“稳态传输”中的滚动方向。ABAQUS SST分析要求你把转动轴定义清楚可以用局部坐标系自动计算角速度。接触地面设为刚体平面固定不动轮胎刚体运动由参考点控制。算出来的结果能直接看接地压力沿周向的分布非常方便。这种情况下你不必给轮胎施加转动位移软件内部会组装出稳态滚动状态省事不少。3. 实操流程走一遍从建模到后处理3.1 快速搭建2D断面模型我以最常见的子午线轮胎为例走一遍完整流程。打开ABAQUS/CAE先在Part模块创建二维轴对称部件类型选“Axisymmetric”变形体。用草图绘制轮胎断面主要几何包括胎冠、胎侧、胎圈、胎腔、轮辋座。如果你的几何来自外部文件导入后务必用“Partition”把断面分割成规则区域尤其是胎圈包裹区域和帘线层区域否则网格质量完全不受控。断面画好后赋予材料。橡胶用超弹性帘线层用Rebar或者嵌入膜单元。在Section里创建剖面指派给对应区域。这里有个经验轮胎断面网格不要用二次三角形单元推荐用一次四边形减缩积分杂交单元CAX4RH。为什么是杂交因为橡胶近似不可压缩普通单元在变形时会产生体积锁定杂交单元专门解决这个痛点。减缩积分则能防止剪切锁死。网格尺寸在胎侧等薄壁区域至少保证两层单元胎冠可以稍微粗一点但胎肩过渡区要加密。3.2 过盈和充气分析步的设置细节装配后定义接触。主面选轮辋外表面从面选胎圈内表面。主面网格要粗于从面这样接触计算更稳定。接触属性里切向给摩擦系数0.6法向“Hard”接触允许分离但不允许穿透。然后在Step-1的接触定义里设置过盈量我通常取1.0 mm加载方式选Ramp。充气压力加载之前先要确定胎压作用面。胎腔表面是轮胎内侧在Load模块创建Pressure选择胎腔表面数值按设计胎压填写。这里有个关键操作压力面的法向必须指向外侧。如果不放心可以在加载后先跑一个充气步在可视化里看胎体是不是均匀向外膨胀一旦发现某个方向凹进去了多半是压力方向反了赶紧改。求解时打开大变形。增量步初始0.01最小1e-8最大0.1允许自动增量。如果遇到过盈不收敛一个技巧是先把过盈量分再小一点比如0.2 mm逐步增加到目标值虽然是麻烦一点但比死磕一个分析步省时间得多。充气完成后最好检查一下接触力。用History Output输出总接触力CFN看看过盈配合后接触力是否随充气压力平滑变化。正常情况下充气后胎圈与轮辋之间的接触力会略有降低但不应变为零。如果变成一个很小的值甚至正负跳变说明过盈量不够或者轮辋形状匹配有问题后续滚动的稳定性会很差。3.3 从2D扩展到3D滚动模型得到充气后的变形结果在ABAQUS里用“Generate Part”将2D变形网格旋转成3D。具体操作是在Job模块提交一个只包含过盈和充气的分析完成后进入Visualization在主变量列表里选位移U当前帧显示的是最后一步的变形结果。然后在Interaction模块或Mesh模块用“Generate 3D mesh”功能旋转角可以选90度或360度周向网格数量根据你要分析的工况来定。自由滚动、接地压力分布建议至少3040个周向单元做瞬态冲击能加密到60以上。转3D之后原先的轴对称单元会变成对应的实体单元CAX4RH转成C3D8RH。注意转出来的3D模型已经带着充气后的应力状态这是好事但你要确保导入3D模型时把初始应力也导进去了。ABAQUS CAE会自动生成一个带有“初始应力状态”的关键词文件一般没问题。如果不放心可以在3D模型里重新做一次充气静力分析让应力稳定下来代价是多花一点时间但心理踏实。滚动模型的装配包括轮胎、地面和轮辋。地面用解析刚体轮辋用离散刚体主从面还是老规矩刚体做主面轮胎做从面。轮胎与地面的摩擦系数在干路面取0.81.0雨天工况另说。轮胎与轮辋的接触沿用2D模型里的设置但过盈已经有了不需要再给一次否则等于双重过盈应力就爆了。3.4 滚动分析的关键参数选择和后处理技巧滚动分析步的设置我建议先做一个短的静力步只施加垂向载荷让轮胎和地面接触建立稳定然后再切换到稳态传输或显式滚动步。垂向载荷施加在轮轴参考点上可以加载到目标轴荷比如5 kN用Ramp方式缓慢加载。加载太快会导致接触振荡尤其是显式分析还可能在轮胎上激起虚假的振动模态。在稳态传输分析步里需要定义地面不动轮胎转动和移动的速度。自由滚动工况下旋转角速度和前进速度不是独立给定的而是满足一个关系接触区的纵向力合力和力矩为零。ABAQUS允许你通过“自由滚动”边界条件自动求解这个状态操作上在SST分析步选择“free rolling”选项给定一个前进速度软件就会迭代求解对应的角速度非常实用。不过要注意自由滚动只能用于对称或者没有侧偏的工况。如果有侧偏角就得指定角速度和前进速度并让软件算侧向力。用Explicit做滚动的话轮胎前进速度不宜给太高否则稳定时间步非常小计算量吓人。一般可以降低速度比如把实际车速30 m/s降为10 m/s同时保持低速下惯性效应不大。如果你需要看到高频特性再考虑真实速度。显式分析还要考虑质量缩放但不要一上来就放大十倍那样惯性力会掩盖真实接触力。推荐把质量缩放目标稳定增量设在1e-6左右观察整体动能占比保持在总能量的5%以内即可。后处理时我最常看的是主变量S应力和CPRESS接触压力。在Visualization模块的Field Output里主变量列表选S然后勾选Mises分量当前帧一般显示最后一步。若想看滚动一周过程中某个截面的应力变化可以输出多条路径的历史变量也可以查看不同分析帧的云图。注意主变量“S”和“Current frame”的区别前者是场输出变量后者是时间帧索引。很多新手问“选中的主变量在当前帧为什么不变化”多半是只输出了最后一帧而没开整个历史帧输出。我在Step模块的Output里习惯把Field Output的频率设为“Last frame”之外的“Every 1/20 frame”这样重启动后能看到整个滚动过程的应力演变排查问题方便得多。4. 典型问题排查与优化建议4.1 过盈和充气步老是不收敛怎么办遇到不收敛先别急着加密网格。我用一套固定排查流程通常能快速定位问题。第一步看Msg文件搜索“negative eigenvalue”或“too many attempts”确认卡在哪个增量步。第二步可视化该增量步的应力/位移云图看是不是某个局部区域出现了畸变或接触穿透。第三步检查接触状态是不是接触从面初始穿透过大。如果初始穿透超过网格尺寸的10%就应该减小过盈量或者调整网格。过盈不收敛还有个常见原因是主从面方向反了。ABAQUS接触对中主面法线方向要指向从面所在侧否则软件会认为接触面法向相反根本无法建立接触。检查一下接触对的法线箭头如果方向反了直接翻转再试。材料参数也容易导致不收敛。如果超弹性参数让体积模量接近完全不可压在隐式分析中会产生极高的体积应变导致单元畸变。所以一定要把D参数设置为一个小但非零的值给材料一点点可压缩性数值稳定性会好很多。我遇到算崩的时候把D1从0.0001改到0.001增量和收敛性立刻改善。4.2 滚动时接触压力像“麻点”一样跳动如果滚动后接触压力云图不是平滑的梯形分布而是密集的不规则斑点大概率是轮胎与地面的网格不匹配导致的接触噪声。解决办法有三一是把地面网格尺寸调的比胎面网格粗糙一点主面不要太细二是在接触属性里增加一个“Contact surface smoothing”选项让从面接触压力更光滑三是检查胎面是否是橡胶层如果胎面单元太薄在弯曲时会产生伪应力此时需要增加胎面区域的网格层数。另外接触压力抖动可能来自减缩积分单元的沙漏模式。如果使用C3D8R单元检查一下沙漏能是否占总内能的比例过高。如果超过5%就需要加密网格或者改用C3D8RH带杂交单元同时开启增强沙漏控制。我在胎面区域一般不用减缩积分换成全积分杂交单元接触压力会干净很多。4.3 滚动后轮胎越跑越扁应力场不对称有时候你发现充气形状挺正常但滚动几帧后胎体局部出现异常压缩左右不对称。这通常不是材料问题而是你的负载边界没有正确耦合到轮轴上。用*COUPLING时耦合节点不要选胎圈最外层节点要选胎圈与轮辋座接触的中心层不然力矩传递路径是歪的一加载就产生侧向分力。还有一点载荷施加方向必须指向全局Y轴负方向假设轮胎装在X轴上如果参考点约束错了自由度轮胎会“别”着转应力自然不对。另一个可能导致不对称的原因是稳态传输中的转动轴方向。ABAQUS SST要求你指定转动轴方向的坐标系通常用圆柱坐标系。如果转动轴定义为全局Z轴而轮胎实际绕X轴转结果就是滚动方向错误。检查一下参考点的局部坐标系方向以及刚体运动的约束这个错误我犯过不止一次现在每次建模都会在预定义场里先画一个转动轴向量确认方向。4.4 显式滚动时间太长如何优化求解时长显式滚动仿真时间成本高这是公认的。优化手段我常试过几个效果最好的是“减小仿真总时长的物理速度”。比如实际滚动0.1秒速度30 m/s那么轮胎走了3米网格必须足够长才能覆盖3米路面网格数量巨大。但如果我降低速度到3 m/s总时长也可以缩短到0.05秒轮胎只走0.15米只要确保该段距离能覆盖接地印迹长度的几倍即可。这样时间增量步基本不变但总增量步数大幅减少计算时间缩短。配合质量缩放把目标稳定时间增量固定到5e-7左右注意不要让总质量增加超过2%。如何判断质量缩放是否过度了在输出里增加一个变量“ALLKE”动能如果动能峰值占内能的比例超过5%说明惯性效应太大了结果不可信。我一般控制在全过程中动能/内能2%稳妥。4.5 后处理里看什么、怎么看滚动分析完成后我习惯先看三个结果一是接地印痕的CPRESS云图判断接触形状是否符合轮胎特征二是轮胎断面的S.Mises应力分布寻找胎肩、胎圈易损区三是轴心的垂向位移下沉量验证载荷—变形曲线是否符合经验值。如果CPRESS云图呈U形或者双峰分布说明胎压偏低或胎冠刚度不足如果呈锥形单峰偏前说明轮胎在驱动或制动工况下纵向滑移分布不均。操作后处理时先把主变量设为S当前帧选择最后一步云图显示Mises应力并勾选“Deformed Variable”为U变形比例设1:1这样能更直观看到滚动状态下的形状。再创建一条路径Path沿胎冠周向输出该路径上的接触压力曲线用来对比不同工况。这些结果写报告时直接用就行比云图更能说明问题。最后关于轮胎仿真我个人踩过的坑再啰嗦几句轮胎仿真入门并不难难的是把每一步的物理背景吃透。你如果只是照着教程点鼠标也许能跑通一两个例子但换一个胎型号、换一个工况照样抓瞎。我个人的建议是先从2D断面开始把过盈和充气算明白记录不同压力下的下沉量和接触力再去碰3D滚动。别一上来就幻想直接做个复杂动态工况那是把前面的所有误差堆在一起最后只能看到一团废数据。做轮胎仿真要多和试验对标。仿真结果是辅助设计不是代替试验。你算出来的接地印痕长度、胎侧变形、胎肩应力都该有对应的实测值。哪怕只是拿尺子量一下真实轮胎的断面尺寸对比一下仿真模型的轮廓偏差也能及时发现建模时的几何错误。很多人仿真结果离了大谱往往不是求解器不行而是几何或材料参数一开始就错了。滚动传荷这一步对网格质量、时间步长、接触稳定性都是三重考验。正常情况下调试一个3D滚动模型从最初建模到得到可信结果消耗掉一两周时间是家常便饭。所以别气馁卡住了就回头检查思路该分步算就分步算该简化就简化。只要过盈、充气、滚动传荷这个大框架不乱每个坑都有解。希望这篇文章能帮你少走几个弯路。后面有机会我再专门写写不同胎压和速度下的滚动响应对比或者轮胎磨损仿真的那些破事。咱们车接着飙。