2. 疲劳基础与HBK分析框架
疲劳分析的起点不是软件操作，而是先判断“该用哪条寿命曲线”。HBK（HBM nCode）把疲劳问题分成两条主线：应力-寿命（S-N）方法适用于高周疲劳，应力水平低于材料屈服强度、循环次数通常在 10⁴ 以上，弹性应变主导；应变-寿命（E-N，ε-N）方法适用于低周疲劳，局部应力已进入塑性，循环次数约在 10²～10⁴ 量级，必须用应变幅描述。选错方法，后面的载荷谱和参数再精细也没有意义。
一、S-N曲线：从材料数据到设计曲线
S-N曲线描述应力幅 Sa 与失效循环数 Nf 的关系，常用幂函数形式表达：
S_a = S_f' · (2N_f)^b
σ_a = σ_f' · (2N_f)^b （应力-寿命，Basquin 形式）
ε_a = σ_f'/E · (2N_f)^b + ε_f' · (2N_f)^c （应变-寿命，Coffin-Manson-Basquin）其中 σf' 为疲劳强度系数，b 为疲劳强度指数（典型 −0.05～−0.12），εf' 为疲劳延性系数，c 为疲劳延性指数（典型 −0.5～−0.7）。在双对数坐标下，S-N 曲线近似为直线；钢材在 10⁶～10⁷ 次附近常出现疲劳极限，曲线转为水平，低于该应力幅理论上不再失效。铝合金等无疲劳极限材料则持续下降，工程上常以 5×10⁸ 次对应的条件疲劳强度作为判据。
实际分析用的是设计曲线，而非材料手册中的均值曲线。HBK 中通过 Material Manager 定义材料时，需要输入上述四个参数，并设置存活率（如 50%、97.7%、99%）与修正方式。常用修正包括：
- Goodman / Gerber 平均应力修正：把非零平均应力折算为等效对称应力幅，Goodman 偏保守，Gerber 偏乐观；
- Haigh 图：以平均应力为横轴、应力幅为纵轴，直观显示修正后的许用区域；
- 表面加工、尺寸、缺口系数：通过
Kf、Kt与表面质量因子修正材料曲线，得到零件曲线。
二、Miner线性累积损伤理论
真实载荷是变幅的，需要把载荷谱分解为若干级恒幅块，再逐级累加损伤。Miner 法则假设：每一级应力幅 Si 作用 ni 次，对应寿命 Ni，则该级损伤为 ni/Ni，总损伤为：
D = Σ (n_i / N_i)
D ≥ 1 → 判定失效（理论临界值）
工程常用临界值：D_cr = 0.3 ~ 0.5（考虑分散性与保守设计）Miner 法则的局限必须清楚：它忽略载荷顺序效应（高-低与低-高序列损伤不同）、忽略低于疲劳极限的循环（实际上小幅循环仍可能造成损伤，尤其在高应力历史之后）、且不区分裂纹萌生与扩展。因此 HBK 中通常配合雨流计数把随机时域信号转换为应力幅-均值-循环数的矩阵，再逐级代入 Miner 累加。
三、HBK 中的模块对应与参数设置
在 HBK nCode DesignLife 中，疲劳分析以“流程链”方式搭建，核心模块与上述理论一一对应：
| 理论环节 | HBK 模块 | 关键参数 |
|---|---|---|
| 载荷输入 | FE Input / Time Series Input | 单位制、通道映射、采样率 |
| 雨流计数 | Rainflow Counting | 计数方式、均值处理、循环矩阵输出 |
| S-N 分析 | Stress Life (SN) | 材料曲线、存活率、Goodman/Gerber 修正、Kf |
| E-N 分析 | Strain Life (EN) | 循环应力-应变曲线、Neuber 修正、SWT 参数 |
| 损伤累加 | Damage Accumulation | Miner 临界值 D_cr、损伤位置输出 |
| 结果查看 | Results / Hot Spot | 损伤云图、寿命云图、危险点定位 |
以 Stress Life 模块为例，典型设置流程是：在 Material 页选择或导入材料，确认 σf'、b 与存活率；在 Loading 页指定载荷通道与缩放因子；在 Analysis 页选择平均应力修正方法（Goodman 或 Gerber），并设定 D_cr（常用 0.3～0.5）；在 Output 页勾选损伤、寿命与安全系数输出。若为低周疲劳，则切换到 Strain Life 模块，额外需要循环应力-应变曲线参数 K'、n'，并选择 Neuber 或 Glinka 缺口修正。
理解这套框架后，后续的载荷谱编辑、雨流计数参数调整和损伤云图判读，才有明确的物理依据。HBK 的价值不在于自动算出数字，而在于把“材料曲线—载荷计数—损伤累加—结果定位”这条链路透明地暴露给分析者，让每一个参数都能追溯到工程假设。