作为车辆与路面接触的唯一部件,轮胎的许多重要性能都体现在动态的使用条件下,如操控性能,磨耗性能,耐冲击性能,排水性能,等等。目前应用 Abaqus/Standard 进行轮胎的静态仿真以及稳态输运仿真,已经成为轮胎业内常用的手段,如进行轮胎的充气轮廓,静态刚度,侧偏刚度计算等[1-2]。但轮胎产品的开发和性能改进,还需要轮胎动态仿真技术的支持和指导,我们应用显式有限元求解器 Abaqus/Explicit,克服了显式计算方法计算时间较长,结果存在振荡等问题,开发了一整套轮胎动态 FEA 仿真技术,包括轮胎动态加载,瞬态滚动,侧偏仿真,驱/制动仿真,摩擦能仿真等。
轮胎动态 FEA 仿真中的挑战
目前轮胎行业内常用的动态仿真方法为,使用 Abaqus/Standard 求解器中的*Steady State Transport 功能,实现稳态的动态仿真。但是使用 Abaqus/Standard 进行轮胎动态 FEA 仿真存在如下挑战:
较强的非线性效应导致隐式方法收敛困难
Abaqus/Standard 为隐式有限元求解器,求解非线性问题使用迭代的方法,而轮胎的动态仿真中一般存在较强的非线性效应,包括:
材料非线性,轮胎的胶料不仅具有超弹性,而且具有粘弹性,当轮胎承受静态的载荷时,其粘弹性效应可以忽略不计,而当轮胎处于动态的加载条件时,胶料会表现出粘弹性的行为,材料非线性可能会导致计算收敛的困难。
几何非线性,轮胎的分析全部为几何大变形分析,而在一些极端的动态工况下,如大角度的侧偏工况,大滑移率的驱/制动工况,冲击碰撞工况等,轮胎有限元模型可能会发生网格的较大变形,乃至畸变,这可能会导致计算收敛的困难。
边界条件非线性,当考虑轮胎的胎面花纹时,轮胎滚动过程中胎面花纹与路面(或介质,如水)发生接触,且胎面花纹中通常存在较细的沟槽,这些沟槽的内壁也有可能发生自接触,当胎面花纹较为复杂时,模型中接触面的数量也会增加,随着加载状态的改变,接触状态也会发生复杂的改变,这也可能导致计算收敛的困难。
计算收敛的困难表现为两个方面,一是计算因不收敛而中止,无法得到所需工况下的计算结果;二是计算中迭代次数过多,导致计算时间过长,而且总迭代次数可能随模型的不同而发生改变,导致总计算时间不可预计。
较大的模型规模导致隐式迭代方法的效率受到限制
当前使用的轮胎有限元模型规模较大,尤其是当轮胎模型中包含复杂胎面花纹时,可达数十万至上百万的单元规模,每次对模型的迭代求解都会花费较长的时间。
综上所述,使用隐式有限元求解器 Abaqus/Standard,进行轮胎动态有限元仿真时受到较大的限制。Abaqus/Explicit 是一款功能强大完备的显式有限元求解器,最终我们选择Abaqus/Explicit 进行轮胎动态有限元仿真。
显式有限元求解器 Abaqus/Explicit 进行轮胎动态 FEA 仿真存在的挑战
显式有限元求解器 Abaqus/Explicit 进行轮胎动态 FEA 仿真也存在一些挑战[3],由于显式方法存在稳定时间增量的限制,一般其稳定时间增量非常小,可能需要几十万甚至上百万个增量,才能模拟所需要的工况历程,因此计算时间较长;另外使用显式方法时,可能发生结果中存在振荡现象的问题。经过我们努力克服这两个问题,可以使用 Abaqus/Explicit,在可接受的计算时间内,得到稳定准确的计算结果。
应用 Abaqus/Explicit 的动态 FEA 仿真技术介绍
我们应用 Abaqus/Explicit,开发了一整套轮胎显式动态 FEA 仿真技术,包括轮胎动态加载,瞬态滚动,侧偏仿真,驱/制动仿真,摩擦能仿真等。
轮胎建模流程
在轮胎显式动态 FEA 仿真技术中,对轮胎进行了专门的模型简化和建模。
在我们的建模流程中,将轮胎的有限元网格模型分为两个部分:轮胎主体部分和胎面花纹部分。对于轮胎主体部分,首先根据轮胎的材料分布图,使用内部建模工具生成轮胎主体部分 2D 模型,再利用轮胎的轴对称性质,得到轮胎主体部分 3D 模型;对于轮胎的花纹部分,首先根据花纹设计图,使用内部建模工具生成单节距胎面花纹模型,再利用胎面花纹的周期对称性质,得到胎面花纹的 3D 模型。最后将轮胎的主体部分和花纹部分绑定组合,得到轮胎的整体有限元网格模型,如图 1 所示。
图 1. 显式轮胎 FEA 模型建模流程
轮胎加载与瞬态滚动仿真
使用 Abaqus/Explicit 进行轮胎准静态加载,既可以解决隐式方法收敛困难的问题,而且具有较高的计算效率,当处理包含胎面花纹的轮胎模型时其优势更加明显。下表为我们使用Abaqus/Explicit 进行轮胎准静态加载计算的四组轮胎案例的接地印痕结果。
表 1. 四组轮胎案例的试验和计算接地印痕结果
轮胎瞬态侧偏仿真
轮胎在侧偏工况下产生所产生的侧向力大小是其重要的操控性能指标。当轮胎的侧偏角较小时,所产生的侧向力基本随着侧偏角的增大而线性增加,而随着侧偏角的继续增大,侧向力增加的幅度不断减小,直到在某个侧偏角下,侧向力达到峰值,继续增大侧偏角,侧向力反而会随之下降,如图 2 所示。
图 2. 侧向力随侧偏角的关系示意图
图 3 为自由滚动下的轮胎接地印痕,图 4 为 6°侧偏角下的接地印痕,从图中可看出,在侧偏角较大时,轮胎的接地印痕形状相对于自由滚动时将发生明显的改变,且胎面将产生明显的变形。
图 3. 自由滚动下的轮胎接地印痕
图 4. 侧偏角为 6°时的接地印痕
分别对轮胎施加不同的侧偏角,可得到各侧偏角下轮胎的稳态侧向力结果,如图 5 所示。将所有侧偏角下的侧向力结果进行汇总,可得到轮胎的侧向力与侧偏角关系曲线,如图 6所示。
图 5. 不同侧偏角下轮胎的稳态侧向力结果示意图
图 6. 轮胎的侧向力与侧偏角关系曲线示意图
轮胎瞬态驱/制动力仿真
轮胎在驱动制动工况下所产生的驱/制动力是其重要的操控性能指标。一般用滑移率SR来定义义轮胎驱/制动的工况,
(为轮胎当前的平动速度所对应的自由滚动角速其中,Ω为轮胎当前的转动角速度,度。当SR>0时轮胎处于驱动的工况,而当SR<0时轮胎处于制动的工况。当轮胎的滑移率较小时,所产生的驱/制动力基本随着滑移率的增大而线性增加,而随着滑移率的继续增大,驱/制动力增加的幅度不断减小,直到在某个滑移率下,驱/制动力达到峰值,继续增大滑移率,驱/制动力反而会随之下降,如图7所示。
图 7. 制动力随滑移率的关系示意图
图8为自由滚动下的轮胎接地印痕,图9为滑移率为-15%的制动工况下的轮胎接地印痕,
图 8. 自由滚动下的轮胎接地印痕示意图
图 9. 滑移率为-15%的制动工况下的轮胎接地印痕示意图
分别对轮胎施加不同的滑移率,可得到各滑移率下轮胎的稳态驱/制动力结果,如图 10所示。将所有滑移率下轮胎的稳态驱/制动力结果进行汇总,可得到轮胎的驱/制动力与滑移率关系曲线,如图 11 所示。
图 10. 不同滑移率下轮胎的稳态制动力结果示意图
图 11. 轮胎的驱/制动力与滑移率关系曲线示意图
轮胎摩擦能仿真
轮胎的磨耗性能是其重要使用性能之一,使用摩擦能仿真技术可以计算在轮胎的某种使用工况下,轮胎胎面在与地面接触的过程中,所产生的摩擦能密度(单位面积上的摩擦能量)的分布情况。使用摩擦能密度结果可有助于预测轮胎的偏磨耗现象。
轮胎的偏磨耗现象可能由多种因素导致,除轮胎本身的结构设计和花纹设计原因之外,还有可能由于不当的车轮定位参数导致,图 12 和图 13 为自由滚动状态和处于某定位参数的轮胎(左前轮位)示意图,由于定位参数的差异,导致轮胎的静态接地印痕产生一定的差异,如图 14 所示。而通过进行摩擦能密度的仿真计算,可以发现当轮胎处于不受驱/制动力的滚动状态时,由于车轮的定位参数,轮胎处于轻微的侧偏和侧倾工况,就可能导致轮胎的胎面摩擦能密度结果产生很大的差异。图 15 为摩擦能计算结果,从结果可以看出,不合理的定位参数将导致轮胎出现偏磨耗,其内侧磨耗速率远高于外侧,将明显影响轮胎的使用寿命。
图 12. 处于不同定位参数的轮胎示意图(后视图),上图为自由滚动轮胎,下图为处于某定位参数下的轮胎
图 13. 处于不同定位参数的轮胎示意图(俯视图),
上图为自由滚动轮胎,下图为处于某定位参数下的轮胎
图 14. 处于不同定位参数的轮胎接地印痕计算结果示意图
图 15. 处于不同定位参数的轮胎摩擦能密度计算结果示意图
其他仿真
除此之外,应用 Abaqus/Explicit 进行的轮胎动态 FEA 仿真技术,还可应用于轮胎方冲击计算,水滑计算,雪地牵引计算等。
总结
应用显式有限元求解器 Abaqus/Explicit,开发了一整套轮胎动态 FEA 仿真技术,包括轮胎动态加载,瞬态滚动,侧偏仿真,驱/制动仿真,摩擦能仿真等。将显式动态计算结果与试验结果及隐式计算结果进行对比验证,通过与试验吻合较好的计算结果,加深对实际机理的理解,并应用于产品实际开发过程,解决若干问题。
资料来源:达索官方
[ABAQUS] Abaqus自带应变应力本构方程吗?Abaqus购买价格咨
2025-06-20
[ABAQUS] ABAQUS壳单元与实体单元怎么设置接触?abaqus广东
2025-06-20
[ABAQUS] Abaqus节点过度约束怎么解决?Abaqus经销商思茂信
2025-06-20
[ABAQUS] Abaqus接触压力不均匀出现的原因?广州Abaqus代理
2025-06-19
[ABAQUS] Abaqus的从面总是穿透怎么处理?Abaqus购买指南
2025-06-19
[ABAQUS] Abaqus/CAE为什么施加应力前后的图重叠了?Abaq
2025-06-19
[ABAQUS] 有限元分析软件Abaqus采购指南-达索代理商思茂信息
2025-06-18
[CST] CST年度更新需额外付费吗?达索 SIMULIA 代理商思
2025-06-18
[有限元知识] 网格细化程度与有限元分析结果准确性有必然的关联性吗?
2025-06-18
[ABAQUS] Abaqus多物理场仿真如何解决复杂工程耦合问题
2025-06-17
2023-08-29
2023-08-24
[ABAQUS] ABAQUS如何建模?ABAQUS有限元分析教程
2023-07-07
[ABAQUS] 有限元分析软件abaqus单位在哪设置?【操作教程】
2023-09-05
[ABAQUS] ABAQUS软件中Interaction功能模块中的绑定接
2023-07-19
[ABAQUS] Abaqus怎么撤回上一步操作?Abauqs教程
2024-05-01
[ABAQUS] abaqus里面s11、s12和u1、u2是什么意思?s和
2023-08-30
[ABAQUS] Abaqus单位对应关系及参数介绍-Abaqus软件
2023-11-20
[ABAQUS] ABAQUS软件教程|场变量输出历史变量输出
2023-07-18
[ABAQUS] Abaqus中的S、U、V、E、CF分别是什么意思?
2024-05-11
[有限元知识] 网格细化程度与有限元分析结果准确性有必然的关联性吗?
2025-06-18
2025-06-16
2025-06-16
[行业资讯] Celanese与SIMULIA合作应对复杂的客户请求
2025-06-16
[行业资讯] 正版Abaqus软件价格及有限元分析软件Abaqus模块解
2025-06-12
2025-06-12
[有限元知识] 如何通过达索代理商购买Abaqus有限元分析软件
2025-06-12
[有限元知识] 收藏!一文搞懂有限元分析的底层逻辑、操作要点和学习路径
2025-06-11
[有限元知识] 有限元分析时,网格尺寸通常如何确定?一般取多大?
2025-06-11
[行业资讯] SIMULIA电磁仿真线上技术汇|专题二:面向汽车电缆设计
2025-06-10
地址: 广州市天河区天河北路663号广东省机械研究所8栋9层 电话:020-38921052 传真:020-38921345 邮箱:thinks@think-s.com
Copyright © 2010-2023 广州思茂信息科技有限公司 All Rights Reserved. 粤ICP备11003060号-2