4. 疲劳分析前的静力学准备
疲劳寿命评估建立在静力学应力场之上,因此静力学准备阶段的质量直接决定疲劳结果的可靠性。本节围绕四个关键环节展开:网格质量、边界条件与载荷、应力奇异与热点识别、子模型技术。它们并非彼此独立,而是构成一条从“几何离散”到“局部精细化”的递进链条。
一、网格质量:疲劳结果的第一道门槛
疲劳损伤由应力幅值驱动,而应力幅值对网格尺寸极为敏感。在应力梯度大的区域(如圆角、孔边、焊缝趾部),粗网格会显著低估峰值应力,导致寿命预测偏长;反之,过度细化又可能引入数值噪声。工程上通常关注三类指标:
- 单元质量:Jacobian 比值、翘曲度、长宽比。实体单元长宽比建议控制在 5 以内,壳单元翘曲角小于 10°。
- 过渡质量:粗细网格交界处应平滑过渡,避免相邻单元尺寸突变超过 2 倍。
- 厚度方向层数:弯曲主导区域建议至少 3 层二阶单元,以捕捉线性应力梯度。
判断网格是否足够,最实用的方法是网格收敛性检查:将关注点(如孔边最大主应力)随网格加密的变化绘制曲线,当连续两次加密的应力变化小于 5% 时,可认为该区域网格收敛。注意,收敛的是“关注点应力”,而非整体模型。
二、边界条件与载荷设置
边界条件决定应力分布形态,载荷决定应力幅值。疲劳分析中常见的错误是边界条件过于刚硬,人为制造出实际不存在的应力集中。例如,将螺栓连接简化为全约束面,会在孔边产生虚假峰值应力。
合理的做法是:
- 约束位置远离关注区:约束应设置在远离疲劳热点的截面,或采用远端位移、弹性支撑等柔化方式。
- 载荷按实际传力路径施加:集中力应通过刚性元或分布耦合施加到实际承载面,避免单点加载。
- 区分恒幅与谱载荷:恒幅疲劳只需一个静力工况;谱载荷疲劳需要多个工况组合,并明确各工况的循环次数与比例关系。
以一个支架焊缝疲劳为例:若实际工况为 10 kN 拉压循环,约束应设置在支架安装底面,载荷通过销轴孔以分布力施加。此时焊缝趾部的应力幅才是可用于疲劳计算的真实值。
三、应力奇异与热点识别
应力奇异是指应力随网格加密而无限增大的现象,通常出现在尖角、点载荷、刚性约束边缘等位置。奇异点的应力值没有物理意义,不能直接用于疲劳计算。识别方法很简单:加密网格后应力持续显著上升且不收敛,即可判定为奇异。
与奇异相对的是热点:热点是真实存在的应力集中区域,其应力随网格加密会收敛到有限值。疲劳分析真正关心的是热点,而非奇异点。工程上常用以下方式区分与处理:
- 对几何尖角进行倒圆,消除人为奇异;
- 对焊缝采用热点应力法,提取距焊趾一定距离(如 0.4t 或 1.0t)处的应力并线性外推;
- 对无法消除的奇异点,改用名义应力法或结构应力法,避开奇异峰值。
四、子模型技术:从整体到局部的桥梁
当整体模型网格较粗、无法准确捕捉热点应力时,子模型技术是最经济的解决方案。其基本流程是:
- 在整体模型中完成静力分析,获得全局位移场;
- 在关注区域切割出子模型,切割边界应远离热点,避免边界效应影响;
- 将整体模型在切割边界上的位移插值到子模型边界,作为驱动条件;
- 对子模型局部加密,重新求解,提取热点应力。
子模型的关键在于切割边界的位置。边界应选在应力梯度平缓处,通常要求边界到热点的距离大于 3 倍壁厚。若边界太靠近热点,插值误差会直接污染热点应力。验证子模型是否正确的一个实用方法:在切割边界附近取若干点,比较整体模型与子模型的应力值,若偏差小于 5%,说明边界传递可靠。
五、静力学准备检查清单
在将静力结果提交疲劳计算前,建议逐项确认:
| 检查项 | 合格标准 | 常见问题 |
|---|---|---|
| 热点区网格 | 连续加密应力变化 < 5% | 粗网格低估峰值 |
| 约束位置 | 远离热点,无虚假应力集中 | 全约束面制造假峰值 |
| 载荷施加 | 分布力/耦合,符合传力路径 | 单点加载引入奇异 |
| 奇异点处理 | 倒圆或改用名义应力法 | 直接取奇异值算寿命 |
| 子模型边界 | 距热点 > 3t,应力偏差 < 5% | 边界太近,插值失真 |
完成以上准备后,静力分析输出的应力幅、平均应力、应力比等参数才具备进入疲劳求解器的资格。下一阶段将进入疲劳材料参数与损伤模型的建立。