什么是应力-应变曲线与材料非线性模型
🙋
老师，这个模拟器里画的应力-应变曲线到底是什么？就是材料被拉断的过程吗？
🎓
简单来说，它就像材料的“体检报告”。横坐标是应变，代表材料被拉长了多少；纵坐标是应力，代表材料内部抵抗变形的力。这条曲线能告诉我们材料是软是硬，是脆是韧。比如在汽车碰撞试验中，工程师就是通过分析不同钢材的应力-应变曲线，来选择既安全又轻量化的材料的。你试着在模拟器里选择“结构钢”和“铸铁”对比一下，马上就能看到它们的巨大差异！
🙋
诶，真的吗？我选了，结构钢在屈服后还有很长一段上升，而铸铁几乎直接就断了！那旁边这些“双线性”、“Ramberg-Osgood”模型又是干嘛的？
🎓
问得好！这些模型就是我们用来“数学描述”那条曲线的工具。双线性模型就像用两根直线去近似，计算简单，工程现场常用在初步设计阶段。而Ramberg-Osgood模型则用一条光滑的曲线来更精确地描述从弹性到塑性的过渡区域，这在分析飞机机翼的疲劳寿命时特别重要。你试着把“材料模型”从线性弹性切换到Ramberg-Osgood，然后拖动“硬化指数n”的滑块，看看曲线的形状是如何连续变化的。
🙋
哦！我看到了，n越大，屈服点附近的“膝盖”弯得越明显。那下面自动算出来的“韧性”又是什么？和曲线有什么关系？
🎓
韧性就是材料在断裂前吸收能量的能力，数值上就等于整条应力-应变曲线下方的面积！面积越大，材料越“韧”，被打断前能做的功越多。比如防弹衣里的材料就需要极高的韧性。你可以试试把“断裂应变”调大，或者把“抗拉强度”调高，观察下方“韧性”数值和曲线下面积的同步增长。这个数值在CAE软件里对预测结构的失效至关重要。
物理模型与关键公式
Ramberg-Osgood 模型：这是一个经典的连续非线性模型，特别擅长描述金属材料从弹性到塑性的平滑过渡。公式由弹性应变和塑性应变两部分相加而成。
$$\varepsilon = \frac{\sigma}{E}+ 0.002\left(\frac{\sigma}{\sigma_y}\right)^n$$
其中，$\varepsilon$ 是总应变，$\sigma$ 是应力，$E$ 是弹性模量（材料刚度），$\sigma_y$ 是屈服强度，$n$ 是硬化指数（控制曲线过渡的陡峭程度）。当 $\sigma = \sigma_y$ 时，塑性应变项恰好为0.002（即0.2%的残余应变，这是工程上定义屈服点的常用标准）。
幂次硬化法则（Hollomon）：这个模型主要用于描述材料在达到屈服后、颈缩前的塑性硬化阶段，形式简洁，在金属成形分析中应用广泛。
$$\sigma = K \varepsilon_p^n$$
其中，$\sigma$ 是真实应力，$\varepsilon_p$ 是塑性应变，$K$ 是强度系数，$n$ 是应变硬化指数。$n$ 值越大，材料在变形过程中强化的能力越强。注意，这里的应变通常指塑性应变部分。
应力-应变曲线的各区域与力学性能
拉伸试验得到的 $\sigma$–$\varepsilon$ 曲线是材料的"指纹"。延性金属呈现以下区域：
- 弹性区：胡克定律 $\sigma = E\varepsilon$ 成立，卸载后恢复，斜率为杨氏模量 $E$。
- 屈服点：开始塑性变形的应力，即屈服强度 $\sigma_y$（软钢有上、下屈服点）。
- 塑性区（加工硬化）：应力随变形上升，峰值为抗拉强度 $\sigma_u$（UTS）。
- 颈缩与断裂：超过 UTS 后局部颈缩，在断裂应变 $\varepsilon_f$（伸长率）处破坏。
曲线下面积为韧性（断裂前吸收的能量），弹性区面积为回弹性。脆性材料（铸铁、混凝土）几乎无屈服，在弹性阶段突然断裂。
0.2% 条件屈服与真应力/真应变
0.2% 条件屈服：对无明显屈服点的材料（如铝合金），取 0.2% 残余应变（$\varepsilon=0.002$）对应的应力为屈服强度——将弹性线平移 $0.002$，与曲线交点即为该值。
名义 vs 真应力/真应变：常规曲线用原始面积与标距计算（名义/工程值）。按瞬时截面计算的真应力/真应变转换如下（颈缩前有效）：
$\sigma_{true} = \sigma_{nom}(1+\varepsilon_{nom}), \qquad \varepsilon_{true} = \ln(1+\varepsilon_{nom})$
真应力大于名义应力，且在 UTS 之后继续上升。塑性分析与有限元采用真应力/真应变。
材料非线性模型
分析中用数学模型逼近真实曲线。本模拟器支持的代表模型：
| 模型 | 公式 | 特点 |
|---|
| 线弹性 | $\sigma = E\varepsilon$ | 仅屈服前，最简单 |
| 双线性（各向同性硬化） | $\sigma=\sigma_y+E_t(\varepsilon-\varepsilon_y)$（屈服后） | 两直线近似，切线模量 $E_t$ |
| Ramberg-Osgood | $\varepsilon=\dfrac{\sigma}{E}+0.002\left(\dfrac{\sigma}{\sigma_y}\right)^{n}$ | 弹塑性过渡平滑，航空材料常用 |
| 幂硬化（Hollomon） | $\sigma = K\,\varepsilon_p^{\,n}$ | 表征塑性区，$n$＝加工硬化指数 |
加工硬化指数 $n$ 越大，均匀伸长越大（成形性越好）；$\sigma_u/\sigma_y$ 比越大越延性。切换预设可比较钢、铝、钛与橡胶状聚合物的曲线形状差异。
现实世界中的应用
汽车安全与轻量化设计：通过对比高强度钢、铝合金和碳纤维复合材料的应力-应变曲线，工程师可以优化车身不同区域的材料。例如，A柱需要高强度和一定韧性（曲线高且面积大），而某些覆盖件则可选用延展性好（曲线长）的铝合金以吸收碰撞能量。
航空航天结构疲劳分析：飞机机翼承受数万次的循环载荷。Ramberg-Osgood模型能精确模拟材料在循环加载下的应力-应变响应，是预测部件疲劳寿命、防止灾难性失效的关键输入，符合ASME等国际规范。
金属板材冲压成形：在制造汽车车门等部件时，需要预测板材在冲压过程中是否会开裂或起皱。幂次硬化模型（Hollomon）的参数（K和n）是CAE冲压仿真软件的核心输入，用于模拟材料的流动和硬化行为。
土木工程与抗震设计：钢筋混凝土结构中，钢筋需要良好的延性（长的塑性平台）以保证地震时结构“裂而不倒”；而混凝土本身的脆性曲线则提醒工程师需要箍筋对其进行约束。模拟器中的“混凝土”预设直观展示了这一特性。
常见误解与注意事项
首先，切勿将“屈服应力等同于材料的最大强度”。实际上，屈服应力是材料开始发生永久变形的参考点，而最大强度是后续加工硬化后达到的“抗拉强度（UTS）”。例如，结构钢的屈服应力可能为350MPa，但其UTS可能超过500MPa。在CAE中若仅输入屈服应力并据此判断“超过此值即失效”，可能会低估构件的承载潜力，导致设计过于保守而重量增加。
其次，避免模型选择的随意泛化。双线性模型虽计算快捷，但盲目用于所有非线性分析存在风险。例如，橡胶及某些高分子材料呈现无明显屈服点的“超弹性”行为。在本工具中选择“橡胶”材料时，可观察到曲线初始阶段即为弯曲形态。对此类材料强行套用双线性模型，将严重误判其变形行为。第一步应是观察材料的内在力学响应。
最后，注意参数“n”的误读。幂律硬化模型中的“n”与Ramberg-Osgood模型中的“n”虽名称相同，物理意义却不同：前者主要表征加工硬化的“难易程度”，后者则主要描述屈服过程的“平缓度”。若在两个模型中均设置“n=0.2”并对比，可明显观察到曲线形状的显著差异。务必查阅CAE软件手册，确认特定材料模型中各参数的物理定义。