3. 循环载荷表征

多轴疲劳寿命分析的起点，不是材料参数，而是载荷本身。单轴问题中，一个应力幅值加一个均值就能描述一个循环；进入多轴后，载荷不再是一个标量，而是随时间变化的应力/应变张量分量组合。若载荷表征不准确，后续的临界面搜索、损伤累积和寿命预测都会失去意义。本节围绕“如何把真实的多轴载荷变成可计算、可计数、可比较的循环块”展开。

一、恒幅与变幅：从单参数到多通道

恒幅加载指各应力分量按同一频率、固定幅值和固定均值循环，例如拉-扭比例加载中，轴向应力 σ(t)=σ_a sin(ωt)+σ_m，剪应力 τ(t)=τ_a sin(ωt)+τ_m。此时每个循环完全相同，寿命分析只需处理一个代表性循环。变幅加载则幅值或均值随时间变化，如汽车悬架承受的路谱、风电主轴承受的风况序列。多轴变幅的难点在于：各通道的幅值变化可能不同步，某一时刻轴向已进入小循环，扭转仍处于大循环，循环的“边界”不再由单一通道决定。

二、比例与非比例加载路径

比例加载指各应力分量随时间成比例变化，即应力向量方向在应力空间中保持不变，路径为一条直线。典型例子是 σ 与 τ 同相位、同频率，比值恒定。非比例加载则路径偏离直线，最常见的是 90° 相位差拉-扭加载：σ(t)=σ_a sin(ωt)，τ(t)=τ_a sin(ωt−90°)=−τ_a cos(ωt)，合成路径为椭圆。非比例加载会引发附加循环硬化，使相同等效应力幅下的疲劳寿命明显低于比例加载，这是多轴疲劳区别于单轴的核心效应之一。

三、相位差与路径形状

相位差 φ 是描述多轴载荷非比例程度的关键参数。设 σ(t)=σ_a sin(ωt)，τ(t)=τ_a sin(ωt−φ)。当 φ=0° 时，两通道同时到达峰值，路径为直线，属于比例加载；当 φ=90° 时，一个通道达峰时另一个为零，路径为椭圆，非比例度最大；当 φ 介于两者之间时，路径为倾斜椭圆。工程中常用非比例度因子 F 来量化：F = (2/π) × (椭圆短轴/长轴)，F=0 对应比例加载，F=1 对应 90° 相位差。该参数直接影响临界面上的剪应变幅值和附加硬化程度。

四、载荷谱与雨流计数在多轴下的扩展

真实工程载荷多为随机变幅，需通过雨流计数提取循环。单轴雨流计数基于“峰谷值序列”和“四点法”判断循环闭合。多轴扩展的核心思路是：先将多通道应力/应变时程转换为一个等效标量通道，再对该等效通道执行雨流计数。常用等效方式包括：

• 等效应力法：按 von Mises 或 Tresca 准则将多轴应力合成为等效应力时程 σ_eq(t)，再计数。该方法计算简单，但会丢失路径形状信息，对非比例加载偏保守或偏危险。
• 临界面法：先搜索损伤最大的材料平面，将各时刻应力投影到该平面，得到平面上的正应力与剪应力时程，再组合成等效损伤参量（如 Fatemi-Socie、Wang-Brown 参量）后计数。该方法能较好保留非比例效应，是多轴疲劳分析的主流。
• 多轴雨流计数直接法：在应力空间中定义循环计数规则，直接对张量路径进行峰谷识别和循环配对，计算量较大，适用于研究性分析。

以 90° 相位差拉-扭为例，若 σ_a=200 MPa、τ_a=120 MPa，按 von Mises 等效应力幅 σ_eq,a=√(σ_a²+3τ_a²)=√(40000+43200)≈288 MPa。若直接对该等效应力时程计数，会得到一个恒幅循环；但实际路径为椭圆，临界面上的剪应变幅值随平面角度变化，损伤最大平面可能并非 σ_eq 最大方向。因此，非比例加载下推荐采用临界面法，先确定临界平面角度 θ_c，再在该平面上执行雨流计数。

五、表征结果如何用于寿命分析

完成上述表征后，每个循环块应包含：等效损伤参量幅值、均值、循环次数、对应的临界面角度以及相位差信息。这些数据直接输入损伤累积模型（如 Miner 线性累积或非线性连续损伤模型）。需要特别注意的是，非比例加载下不同循环块之间的相互作用可能不可忽略，简单线性叠加会带来偏差。若载荷谱中比例与非比例段交替出现，建议分段表征、分段计数，再按顺序累积损伤，而不是将整段时程混合计数。

总结来说，多轴循环载荷表征的完整链条是：多通道时程 → 路径判断（比例/非比例）→ 相位差量化 → 等效转换或临界面搜索 → 雨流计数 → 循环块输出。每一步的精度都会传递到最终寿命结果，其中临界面角度和相位差处理是误差最敏感的两个环节。