4. 材料循环本构模型

在高低周联合分析中，材料循环本构模型决定了应力—应变滞回环的形状、循环软化/硬化趋势以及平均应力松弛行为。若本构模型选择不当，即使载荷谱和疲劳损伤模型正确，寿命预测也可能出现数倍偏差。本节从四种常用模型出发，说明其数学结构、物理含义、适用场景，以及在联合分析中如何用同一组试验数据完成参数标定。

一、四种模型的定位与数学结构

各向同性硬化描述屈服面尺寸随累积塑性应变扩大或缩小，不改变屈服面中心位置。其背应力为零，屈服函数为 f = |σ| − R − σ_y0，其中 R 为各向同性硬化变量，常用饱和形式 R = Q(1 − e^(−b·p))，p 为累积塑性应变。该模型能表现循环硬化或软化，但无法描述包辛格效应引起的反向屈服提前。

随动硬化则保持屈服面尺寸不变，让屈服面中心随塑性流动移动，背应力 α 记录加载历史方向。线性随动硬化取 dα = (2/3)C·dε^p，反向加载时屈服应力为 σ_y0 + α，从而正确给出滞回环的“尖角”与包辛格效应。但线性随动硬化给出的滞回环过窄，循环稳定后平均应力不松弛，与很多金属的实际行为不符。

Chaboche模型是随动硬化的多背应力叠加形式：α = Σ α_i，每个分量满足 dα_i = (2/3)C_i·dε^p − γ_i·α_i·dp。多个分量分别对应短程快速饱和、中程演化和长程线性段，能同时拟合单调曲线、滞回环形状和循环松弛。工程中常用 2～3 个背应力分量，配合一个各向同性分量 R 描述循环软化/硬化。

Ramberg-Osgood关系是另一条路线，直接以总应变描述应力：ε = σ/E + (σ/K)^(1/n)。它不追踪加载历史，而是给出循环稳定滞回环的骨架曲线，常用于高周疲劳的应变—寿命曲线配套、缺口局部应力应变近似（Neuber 修正）以及低周疲劳的稳定环描述。其优点是参数少、标定快；缺点是无法直接描述循环瞬态、平均应力松弛和加载路径相关效应。

二、参数标定的实际操作

标定应遵循“先弹性、再屈服、后硬化”的顺序，避免参数间相互补偿。以某低合金钢为例（以下数值为教学示例，非实测数据）：E = 210 GPa，ν = 0.3，σ_y0 = 350 MPa。单调拉伸曲线在塑性应变 1% 处应力约 480 MPa，在 5% 处约 560 MPa，据此可先估计各向同性饱和值 Q ≈ 120 MPa，b ≈ 8。

Chaboche 背应力分量的标定通常分三步。第一分量取较大的 γ₁（如 800～1500）和适中的 C₁，用于拟合反向加载后 0.1%～0.3% 塑性应变内的快速硬化；第二分量取中等 γ₂（如 100～300），拟合滞回环中段的曲率；第三分量取 γ₃ 接近 0（如 5～20），提供大应变下的线性硬化斜率。一个可用的初值是 C₁ = 60000 MPa，γ₁ = 1000；C₂ = 30000 MPa，γ₂ = 200；C₃ = 8000 MPa，γ₃ = 10。标定后需检查：在 Δε/2 = 0.5% 的稳定环上，模型峰值应力与试验峰值应力偏差应控制在 5% 以内，且卸载段与再加载段的应力—应变曲线形状一致。

Ramberg-Osgood 参数的获取更直接。取循环稳定后的滞回环顶点，以塑性应变 ε_p = ε − σ/E 为横坐标、应力 σ 为纵坐标，在双对数坐标中拟合直线，斜率即 n，截距给出 K。对上述示例，若稳定环在 ε_p = 0.2% 时 σ = 420 MPa，在 ε_p = 1% 时 σ = 520 MPa，则 n = ln(520/420)/ln(1/0.2) ≈ 0.134，K = 420/(0.002^0.134) ≈ 950 MPa。该组参数可直接用于高周疲劳的局部应力应变修正，也可作为低周分析中稳定环的快速校核。

三、在高低周联合分析中的使用策略

高低周联合载荷谱通常包含两类循环：小幅值高周循环（Δε/2 < 0.2%，寿命 > 10⁴）和大幅值低周循环（Δε/2 > 0.5%，寿命 < 10⁴）。建议采用“双模型分工”策略：低周段用 Chaboche 模型进行逐循环瞬态分析，捕捉循环软化、平均应力松弛和棘轮应变；高周段用 Ramberg-Osgood 稳定环配合 Basquin/Coffin-Manson 方程进行损伤累积。两者在过渡幅值处需做一致性检查——即用 Chaboche 模型算出的稳定环峰值应力与 R-O 关系预测值之差应小于 3%。

四、标定结果检查表