去年底,一个衬套的耐久性开发项目走到收尾阶段。
台架按标准块谱跑的,五个载荷级。仿真结果出来,裂纹起裂位置落在橡胶与金属骨架的粘接边缘,和台架拆解看到的位置一致;开裂走向和模式也对得上;寿命预测值落在实测值的 1.5 倍分散线以内。
客户随后补发了一段从试验场采集的路谱,希望用同一套模型再算一遍。
同一个件,同一套本构参数,同一张网格,边界条件一个字没动。只把载荷输入从块谱换成路谱。结果出来,寿命预测值往危险方向偏了将近 40%。
先怀疑参数。本构参数核对了两遍,数值没动过。网格加密一级,结果不变。边界条件逐条复对,和上一轮一致。路谱本身也查了,采样率正常,峰值没超材料的强度上限。流程上一处错都找不出来。
示例:原始路谱
示例:块普载荷
问题出在哪了?
我们最后怀疑到一个方向:载荷加载的先后顺序变了。
台架块谱是固定次序跑的,仿真参数隐式地拟合了那个特定次序。换成随机路谱,大载荷和小载荷的排列全部打乱。而橡胶跟金属不一样,它的机械状态会记住之前发生了什么。
部分工程师基于金属疲劳仿真经验进行橡胶寿命预测,Miner线性累积损伤法则假设每个循环的损伤贡献相互独立,与加载顺序无关。但橡胶不服这个假设。
我们顺着这个方向往下挖,发现了一条连资深工程师都不太谈论的逻辑链:Mullins效应 → 能量释放率漂移 → 裂纹扩展加速。
下面展开说。
你省掉的时序,材料都记着呢
金属疲劳用Miner线性累积是个好东西。ABCD四个载荷块,不管怎么排列,算出来的总损伤一样。工程师只管把块谱打出来,顺序不用操心。
橡胶能不能也这么干?
有人做过对照验证,用同一批哑铃形橡胶试样,两组路径的载荷幅值范围一模一样,只换顺序:一组从大应变降到小应变,另一组从小应变升到大应变。后者的实测寿命比前者高了不止一倍。
就这么简单的一个交换。载荷总量不变,排列顺序一变,寿命差了1倍多。
这说明橡胶疲劳不是简单的损伤叠加。某种物理机制在上一刻留下的痕迹,改变了这一刻的损伤效率。
Mullins效应:
大载荷走过的路,后面替你算账
做橡胶仿真的都知道Mullins效应。填充橡胶第一次拉到某个应变,卸载后再拉,刚度掉一截,俗称应力软化。
台架跑几十个循环之后应力幅值基本稳下来,通常比首圈低10%到30%。常规做法是取稳定后的响应做分析输入,Mullins效应到此为止,不再细想了。
循环软化的应力-应变响应图,以及其在五个不同应变水平下峰值应力的演变情况
等幅工况下这样处理是合理的,所有循环都在同一个软化平面上跑。
变幅工况就不一样了。
路上的减振件不是等幅加载。过减速带那一瞬间,载荷从巡航水平猛跳几倍,材料当场进一步软化。然后车回到平路,正常水平的小载荷循环,面对的是一个刚被软化了的结构。局部刚度更低,变形范围更大。
对裂纹尖端来说,这直接意味着什么?
能量释放率变了。
裂纹扩展速率dc/dN跟能量释放率T之间是幂函数关系。T变大一点,裂纹扩展速率就快一截。问题在于这个T不是常数,它被前面的载荷改了。
串起来的逻辑是这样:大载荷(载荷控制)不只自己制造了损伤,还通过Mullins软化,把后续每一个小载荷造成的裂纹扩展都给提速了。这就是加载次序效应的核心机制。
不是Miner算错了,是Miner的底层假设不适用于橡胶。它假设的是每个循环打在一个状态不变的材料上。实际上每个循环都在一个被历史修改过的材料上运行。
Mars和Fatemi早就提过,
只是没被重视
加载次序效应不是新发现。
Mars和Fatemi在2004年的橡胶疲劳综述里就讨论过。变幅载荷下Mullins效应造成的材料状态变化,会显著影响能量释放率,进而造成疲劳寿命预测的偏差。结合工程验证数据,在某些工况下偏差可能超过两倍。
2004年,这篇文章到现在整整22年了。
为什么学界知道,工程界不太提?
一个原因是常规做法确实能cover大部分工况。等幅台架、固定次序的块谱,这些场景下Mullins效应稳定在一个平面上,次序偏差不会暴露。问题只有在客户掏出路谱、打乱加载历史的时候才出现。
另一个原因是处理起来麻烦。要让仿真记住加载历史,意味着不能再用传统两步法走到底。需要换一套底层逻辑。这在传统仿真工具里不太容易实现。
但现在可以做了。Endurica增量疲劳求解器就是解决这个问题来的。后文讲做法,先说怎么判断自己有没有踩坑。
Endurica DT 多工况疲劳累积损伤求解流程
三种最容易中招的产品
不是所有橡胶件都受加载次序影响。如果你的产品载荷基本等幅、峰值不高,传统方法够用。
但下面三种类型跑不掉。
底盘衬套
轴向、径向、扭转三个方向不是同时来的。过大坑时径向载荷冲上去,跑平路时轴向扭转载荷占主导。先过大坑再跑平路,跟先跑平路再过大坑,橡胶的软化程度完全不同。有文献报道,同批件换一下加载次序,裂纹萌生位置不变,疲劳寿命差异明显。
轮胎
胎面每转一圈一个循环。紧急制动时一个大载荷插 进稳态巡航序列里。同一规格轮胎跑稳态和跑混合工况,胎面裂纹扩展速率差1.5到2倍。天然橡胶含量越高,应变结晶越明显,次序敏感性可能更突出。
动力总成悬置
点火激励和路面随机振动同时作用。发动机先高转再怠速,跟先怠速再高转,悬置橡胶的历史软化状态不同,当前循环的刚度特性和局部应变幅值就不同。
三个场景的共同特征:载荷不是稳态的、不是单一方向的、不是等幅的。如果你的仿真用的参数是从等幅循环标定的,疲劳计算用的是传统Miner,那偏差不是意外,是必然。
让仿真记住历史:
三步就够了
传统做法是两步走:FEA算应力应变,提取危险点载荷谱,疲劳后处理。FEA那步用的是初始材料参数,它不知道材料经历几千次循环后会软化成什么样。
要让仿真记住加载历史,逻辑上需要三样东西。
第一,保留时序。不要把路谱打成等幅块谱。大载荷和小载荷谁先来谁后到,这个顺序信息不能丢。丢了就等于默认次序不影响,而我们已经知道影响很大。
第二,材料模型里把Mullins软化带上。(MR模型的Mullins参数r需要根据实验数据标定,通常取值大于1以保证合理的软化趋势。)常规做法是在超弹性模型中加损伤变量,描述材料应变历史下的软化状态。这个不是新理论,Abaqus里就能做,很多工程师只用了超弹性部分,Mullins参数空着没填。
第三,疲劳计算用裂纹扩展法。裂纹扩展法直接算dc/dN,它对能量释放率T(Tearing Energy)的变化天然敏感。加载历史改了T,裂纹扩展速率自动跟着变,不需要加校正系数。
这三条合在一起,就是增量疲劳求解器的底层逻辑。每完成一个载荷循环,更新材料状态(软化程度、裂纹长度、当前刚度),下一个循环在更新后的状态上继续。次序一变,能量释放率的变化就被自动捕获。
从真实路谱到寿命预测的完整仿真工作流
关键认知:不是在S-N方法上打个补丁。是换一套疲劳计算的底层框架。不再假设损伤跟次序无关,而是让每个循环的起点状态由历史决定。
如果你不确定自己的产品受不受加载次序影响,可以先做一个零成本的验证。
拿出你的随机路谱,跑一版原版仿真。然后把同一个路谱拆成大载荷段和小载荷段,交换前后顺序再跑一版。对比两版的裂纹扩展轨迹和最终寿命。
如果寿命差异超过一倍,说明加载次序在你的产品上不可忽略。如果差异很小,你可以放心继续用传统方法。
不增加新的试验成本,只多跑一轮仿真。
总结
从我们的经验看,很多时候不是参数调得不够准,而是计算框架少了一个维度。
对于变幅载荷下的橡胶件,你省掉的时序信息,材料都替你记着呢。损失的是仿真精度,买单的是产品可靠性。
把加载历史补进去,可能参数都不用动,寿命预测就能往前走一大步。
关注易瑞博科技,持续获取橡胶疲劳仿真与材料测试的干货;菜单可按主题回看历史文章。