超越理想弹性

Learning objectives / 学习目标

完成本课后，你应当能够：

1. 正确阅读工程应力—应变曲线，区分比例极限、弹性极限、屈服、极限抗拉强度、颈缩与断裂。
2. 在均匀变形阶段转换工程应力/应变与真应力/对数应变，说明颈缩后为何需局部截面测量。
3. 用 Maxwell、Kelvin--Voigt 与标准线性固体的时间尺度解释蠕变、应力松弛与滞回。
4. 从 Euler--Bernoulli 假设得到 \(\sigma=-My/I\)、\(EI v''''=q\)，并计算典型梁挠度。
5. 使用带端部有效长度系数的 Euler 屈曲载荷，区分材料强度失效与几何稳定性失效。
6. 把本构关系视为带状态、历史、温度和加载速率的模型，并明确校准域与外推风险。

Prerequisites / 先修知识

• 应力、应变与线弹性：掌握 \(E,G,K,\nu\)、弹性能密度及其小应变边界。
• 能量景观与稳定性：会用平衡点与二阶变化区分稳定和失稳。
• 需要常系数一阶微分方程、截面二次矩和函数曲线面积。

0. Notation and modeling contract / 记号与建模契约

除非另行声明，拉应力为正、梁挠度与载荷符号服从相应小节的坐标约定；同一公式中的工程量、真量和局部裂尖量不得混用。

0.1 Four ways the ideal spring picture can fail / 理想弹簧图景的四种失效

Figure 0.1 — 相同的“变形很大”可能来自四种不同机制，必须先诊断再选模型。

Knowledge 1 — Elastic limit, yield, plasticity, and fracture / 弹性极限、屈服、塑性与断裂

Knowledge ID: P01-U06-L03-K01

1.1 一条曲线上的不同事件

工程应力与工程应变以原始几何量定义：

\[\sigma_e=\frac{F}{A_0}, \qquad \varepsilon_e=\frac{L-L_0}{L_0}.\]

典型韧性金属曲线依次出现：近似线性区；比例极限；卸载仍可恢复的非线性弹性区（有些材料很窄）；屈服与塑性流动；加工硬化；工程应力最大值 UTS；颈缩；断裂。比例极限、弹性极限和屈服强度不是同一个概念，也并非所有材料都有清晰屈服平台。

无明显屈服点时常报告 \(0.2\%\) 偏移屈服强度：作一条斜率为 \(E\)、从 \(\varepsilon=0.002\) 出发的直线；该直线与实验曲线交点所定义的应力称为证明应力（proof stress，或 \(0.2\%\) 偏移屈服强度）。它是约定的可重复指标，不是微观机制突然切换的普适点。

Figure 1.1 — 典型韧性金属工程应力—应变曲线；各标志点代表不同定义。

1.2 工程量与真量

真应力使用瞬时截面积，真应变用连续增量叠加：

\[\sigma_t=\frac{F}{A}, \qquad \varepsilon_t=\int_{L_0}^{L}\frac{d\ell}{\ell} =\ln\frac{L}{L_0}=\ln(1+\varepsilon_e).\]

若塑性流动近似不可压缩且变形仍均匀，\(AL\approx A_0L_0\)，因此

\[\sigma_t\approx\sigma_e(1+\varepsilon_e).\]

颈缩后截面和应变不再均匀，上式不能用试样平均伸长恢复颈部局部真应力；必须测量局部面积，并可能修正三轴应力状态。

Figure 1.2 — 颈缩后工程应力下降不表示颈部材料的局部真应力必然下降。

1.3 卸载、永久应变与韧性

理想弹塑性或带硬化模型中，卸载近似以弹性斜率 \(E\) 返回，应变分解为

\[\varepsilon=\varepsilon^e+\varepsilon^p.\]

卸载到 \(\sigma=0\) 后保留 \(\varepsilon^p\)。屈服强度衡量开始显著塑性流动的门槛；刚度由 \(E\) 衡量；强度、刚度与韧性不可混为一谈。单位体积吸收到断裂的能量近似为整条曲线面积 \(\int_0^{\varepsilon_f}\sigma\,d\varepsilon\)，常称韧性（toughness）。

Figure 1.3 — 塑性加载后的卸载留下永久应变；恢复部分是弹性应变。

1.4 断裂由缺陷与能量共同控制

裂纹尖端会把名义应力放大。线弹性断裂力学的 Mode-I 应力强度因子常写成

\[K_I=Y\sigma\sqrt{\pi a},\]

其中 \(a\) 是从裂纹口到裂尖的裂纹尺度，\(Y\) 是无量纲几何因子。给定试样厚度与几何下发生失稳的临界值写作 \(K_c\)，它可以随厚度和约束状态改变；只有 Mode-I 试样满足平面应变的厚度、裂纹与剩余韧带尺寸有效性要求，并处于小规模屈服时，测得的临界值才可报告为材料的 \(K_{IC}\)。因此，\(K_I\ge K_c\) 是针对给定构型的失稳判据，不能把任意薄试样的 \(K_c\) 都称作 \(K_{IC}\)。

Figure 1.4 — 裂纹尺寸进入失效判据；相同名义应力不保证相同安全裕度。

Worked Example 1 — Engineering versus true measures / 工程量与真量

拉伸试样 \(A_0=200\ \mathrm{mm^2}\)，均匀阶段 \(L/L_0=1.20\)，载荷 \(F=40.0\ \mathrm{kN}\)。假设塑性体积近似不变：

\[\sigma_e=\frac{40.0\times10^3}{200\times10^{-6}} =200\ \mathrm{MPa}, \qquad \varepsilon_e=0.200.\]

\[\varepsilon_t=\ln(1.20)=0.1823, \qquad A=A_0/1.20=166.7\ \mathrm{mm^2}.\]

\[\sigma_t=\frac{40.0\ \mathrm{kN}}{166.7\ \mathrm{mm^2}} =240\ \mathrm{MPa} =\sigma_e(1+\varepsilon_e).\]

真应变小于同一伸长对应的工程应变，因为每一增量用不断增长的当前长度归一化；真应力更大，因为使用缩小的当前面积。

Worked Example 2 — The 0.2% offset construction / 偏移屈服构造

某简化双线性材料：\(E=200\ \mathrm{GPa}\)，在 \(\varepsilon_0=0.00125\)、\(\sigma_0=250\ \mathrm{MPa}\) 后以切线模量 \(H=1.00\ \mathrm{GPa}\) 硬化。塑性段

\[\sigma=250\ \mathrm{MPa}+H(\varepsilon-0.00125).\]

\(0.2\%\) 偏移线为 \(\sigma=E(\varepsilon-0.002)\)。联立得

\[\varepsilon_y=0.003260, \qquad \sigma_{0.2}=252\ \mathrm{MPa}.\]

交点略高于初始拐点，符合正硬化斜率；结果是给定偏移规则下的证明应力，而不是“材料在原子尺度突然屈服”的唯一真值。

Knowledge 2 — Viscoelastic response and hysteresis / 黏弹性响应与滞后

Knowledge ID: P01-U06-L03-K02

2.1 弹簧与黏壶是最小元件

弹簧给 \(\sigma=E\varepsilon\)，黏壶（dashpot）给 \(\sigma=\eta\dot\varepsilon\)。把它们串联得到 Maxwell 模型：应力相同、总应变相加，故

\[\dot\varepsilon=\frac{\dot\sigma}{E}+\frac{\sigma}{\eta}.\]

施加并保持应变 \(\varepsilon_0\) 时，发生指数应力松弛：

\[\sigma(t)=E\varepsilon_0e^{-t/\tau_M}, \qquad \tau_M=\frac{\eta}{E}.\]

保持应力 \(\sigma_0\) 时，Maxwell 模型蠕变为

\[\varepsilon(t)=\frac{\sigma_0}{E}+\frac{\sigma_0}{\eta}t.\]

Figure 2.1 — Maxwell 模型能表示持续流动与应力松弛。

Figure 2.2 — Maxwell 应力松弛：一个时间常数后剩下 \(e^{-1}\)。

2.2 Kelvin--Voigt 与延迟弹性

弹簧与黏壶并联得到 Kelvin--Voigt 模型：应变相同、应力相加，

\[\sigma=E\varepsilon+\eta\dot\varepsilon.\]

对阶跃恒应力 \(\sigma_0\)：

\[\varepsilon(t)=\frac{\sigma_0}{E} \left(1-e^{-t/\tau_K}\right), \qquad \tau_K=\frac{\eta}{E}.\]

它能表示延迟蠕变，却不能表示瞬时弹性跳变或恒应变下有限初始应力松弛。标准线性固体（standard linear solid）增加一个弹簧，可同时给有限瞬时与长期模量，是更真实但仍线性的最小模型之一。

Figure 2.3 — Kelvin--Voigt 蠕变渐近接近 \(\sigma_0/E\)，没有瞬时应变跳跃。

2.3 滞回与时间尺度

周期加载中，黏性耗散使卸载路径不同于加载路径。每周期单位体积耗散能为

\[W_{\mathrm{diss}}=\oint\sigma\,d\varepsilon,\]

即滞回环的有向面积。用材料松弛时间 \(\tau\) 与观察时间 \(T\) 定义 Deborah number

\[\mathrm{De}=\frac{\tau}{T}.\]

\(\mathrm{De}\ll1\) 时材料在观察期间有时间松弛，表现更流体化；\(\mathrm{De}\gg1\) 时内部结构来不及松弛，表现更固体化。这是尺度比较，不是材料的永久标签。

Figure 2.4 — 滞回环面积记录一个循环内不可恢复的机械能。

Worked Example 3 — Stress relaxation / 应力松弛

Maxwell 材料 \(E=2.00\ \mathrm{GPa}\)、\(\eta=200\ \mathrm{GPa\,s}\)，突然施加并保持 \(\varepsilon_0=0.0100\)。

\[\tau_M=\eta/E=100\ \mathrm s, \qquad \sigma(0^+)=E\varepsilon_0=20.0\ \mathrm{MPa}.\]

在 \(t=200\ \mathrm s=2\tau_M\)：

\[\sigma(200)=20.0e^{-2}=2.71\ \mathrm{MPa}.\]

量纲检查：\(\eta/E\) 的单位是秒；指数无量纲。若实验显示长期应力趋于非零平台，单一 Maxwell 模型不足，应考虑标准线性固体等模型。

Knowledge 3 — Bending, buckling, and stability / 弯曲、屈曲与稳定性

Knowledge ID: P01-U06-L03-K03

3.1 Euler--Bernoulli 梁：截面保持平面

取 \(x\) 沿梁向右，截面坐标 \(y\)、挠度 \(v\) 与分布载荷 \(q\) 均以向上为正；定义正弯矩 \(M>0\) 使 \(y>0\) 的纤维受压，并采用与此相容的广义剪力符号，使局部平衡写成 \(M''=q\)。对细长梁、小斜率、线弹性、截面在弯曲后仍保持平面并垂直于中性轴的模型，距中性轴 \(y\) 的轴向应变近似

\[\varepsilon_{xx}=-y\kappa, \qquad \kappa\approx\frac{d^2v}{dx^2}.\]

由 \(\sigma_{xx}=E\varepsilon_{xx}\) 和 \(M=\int_A(-y\sigma_{xx})\,dA\) 得

\[\boxed{\sigma_{xx}=-\frac{My}{I}}, \qquad \boxed{M=EI\kappa}.\]

若采用一致符号约定，横向载荷 \(q(x)\) 与挠度满足

\[EI\frac{d^4v}{dx^4}=q(x).\]

Figure 3.1 — 纯弯曲中轴向应变随离中性轴距离线性变化。

3.2 典型挠度与模型条件

等截面悬臂梁端部集中力 \(P\) 的自由端挠度大小

\[\delta_{\mathrm{tip}}=\frac{PL^3}{3EI}.\]

长度三次方使细长梁对尺度非常敏感；截面二次矩 \(I=\int y^2dA\) 强调把材料布置得远离中性轴。若梁很短厚、剪切变形显著，应改用 Timoshenko 梁；若挠度/转角大，应使用几何非线性理论。

Figure 3.2 — 悬臂梁小挠度曲线；虚线是未变形中性轴。

3.3 Euler 屈曲：强度尚足也可能失稳

理想直、细长、线弹性柱受通过形心的轴向压力时，微小侧向扰动可在临界载荷达到非平凡平衡模式：

\[\boxed{P_{\mathrm{cr}}=\frac{\pi^2EI}{(KL)^2}}.\]

\(K\) 是端部有效长度系数：铰—铰 \(K=1\)，固—固 \(K=0.5\)，固—自由 \(K=2\)，固—铰约 \(K=0.699\)。应选择最弱弯曲轴的 \(I\)。实际初弯曲、偏心、残余应力、材料非线性与端部柔度会使响应在理想临界点前就显著侧弯。

Figure 3.3 — 端部转动与平移约束改变有效屈曲长度和临界载荷。

Worked Example 4 — Cantilever stress and deflection / 悬臂梁应力与挠度

铝悬臂梁 \(L=1.00\ \mathrm m\)，矩形截面宽 \(b=40.0\ \mathrm{mm}\)、高 \(h=80.0\ \mathrm{mm}\)，\(E=70.0\ \mathrm{GPa}\)，端载 \(P=500\ \mathrm N\)。

\[I=\frac{bh^3}{12} =\frac{0.0400(0.0800)^3}{12} =1.707\times10^{-6}\ \mathrm{m^4}.\]

固定端 \(M_{\max}=PL=500\ \mathrm{N\,m}\)，\(c=h/2=0.0400\ \mathrm m\)：

\[\sigma_{\max}=\frac{M_{\max}c}{I} =11.7\ \mathrm{MPa}.\]

\[\delta_{\mathrm{tip}}=\frac{PL^3}{3EI} =1.40\times10^{-3}\ \mathrm m=1.40\ \mathrm{mm}.\]

\(\delta/L=1.40\times10^{-3}\ll1\)，小挠度假设自洽；这并未检查局部应力集中或剪切变形。

Worked Example 5 — End conditions control buckling / 端部条件控制屈曲

钢柱 \(E=200\ \mathrm{GPa}\)、\(I=8.00\times10^{-6}\ \mathrm{m^4}\)、\(L=3.00\ \mathrm m\)。铰—铰时

\[P_{\mathrm{cr}}=\frac{\pi^2(200\times10^9)(8.00\times10^{-6})}{3.00^2} =1.75\ \mathrm{MN}.\]

固—自由 \(K=2\) 时为 \(1.75/4=0.439\ \mathrm{MN}\)；固—固 \(K=0.5\) 时为 \(4(1.75)=7.02\ \mathrm{MN}\)。仅改变理想端部条件就产生 \(16\) 倍跨度，说明不能看到“柱长 \(L\)”就默认 \(K=1\)。

Knowledge 4 — Constitutive relations and material models / 本构关系与材料模型

Knowledge ID: P01-U06-L03-K04

4.1 本构不是守恒定律

动量、角动量和能量平衡具有普适结构；本构关系描述特定材料在特定状态下如何响应。一般可写成

\[\boldsymbol\sigma(t) =\mathcal F\!\left[ \boldsymbol\varepsilon(s),\dot{\boldsymbol\varepsilon}(s),T(s), \alpha(s);\ 0\le s\le t \right],\]

其中 \(\alpha\) 是硬化、损伤等内部变量。最简单各向同性小应变线弹性为

\[\boldsymbol\sigma =2G\boldsymbol\varepsilon +\lambda\operatorname{tr}(\boldsymbol\varepsilon)\mathbf I, \qquad \lambda=\frac{E\nu}{(1+\nu)(1-2\nu)}.\]

塑性模型需屈服函数、流动法则与硬化法则；黏弹模型需时间核或微分方程；超弹模型从应变能函数导出有限变形应力；损伤与断裂模型还需缺陷演化。方程更多并不自动更真实，参数必须可由独立实验识别。

Figure 4.1 — 模型选择由可观测机制驱动，而不是由公式熟悉程度驱动。

4.2 Calibration domain and falsification / 校准域与可证伪性

一个可用本构模型必须说明：

1. 输入与状态： 应变、应变率、温度、湿度、加载方向、内部变量。
2. 参数来源： 哪些试验、哪个应变/速率/温度范围。
3. 客观预测： 不参与拟合的数据是否能被预测。
4. 热力学约束： 弹性能非负，耗散不应违反第二定律。
5. 失效标准： 何时停止相信连续、均匀、无损伤假设。

“曲线拟合得很好”不足以证明机制正确；用高阶多项式拟合单调加载数据，可能在卸载、循环或不同温度下给出荒谬预测。

5. Lesson synthesis / 全课整合

5.1 Model decision table

5.2 Distinguish three kinds of “failure”

• Material yielding: local stress reaches a constitutive threshold; geometry may still be stable.
• Fracture: a crack consumes energy and creates new surface; nominal stress may be below yield.
• Buckling: equilibrium path loses stability; material may remain nearly elastic at onset.

These mechanisms can interact. A buckled column can yield in bending; cyclic plasticity can nucleate cracks; a viscoelastic polymer column can buckle with a time-dependent effective stiffness.

5.3 Model boundaries / 模型边界

• Engineering-to-true conversion assumes uniform deformation and an appropriate volume-change model.
• \(0.2\%\) offset is a reporting convention, not a universal microscopic boundary.
• \(K_c\) 是给定试样厚度、几何和约束下的临界值；只有 Mode-I 平面应变尺寸有效性与小规模屈服条件都满足时，才可把它作为材料常数 \(K_{IC}\)。
• Maxwell and Kelvin--Voigt elements are linear, one-time-scale idealizations; real polymers need spectra of relaxation times and temperature dependence.
• Euler--Bernoulli neglects transverse shear and rotary effects; \(EI\) may vary with \(x\), time, temperature or damage.
• Euler buckling assumes an ideal straight centered column; design strength is normally lower.
• Constitutive parameters are state-dependent data, not timeless constants detached from test conditions.

6. Misconception clinic / 易错概念诊所

7. Practice and retrieval / 练习与检索

A. Conceptual questions

1. 区分比例极限、弹性极限与屈服强度。
2. 为什么颈缩后 \(\sigma_t=\sigma_e(1+\varepsilon_e)\) 的平均转换失效？
3. 材料刚度、强度和韧性分别由什么量描述？
4. Maxwell 与 Kelvin--Voigt 模型分别擅长描述哪类阶跃试验？
5. 为什么滞回环面积具有 \(\mathrm{J/m^3}\) 的单位？
6. 一根柱在材料仍线弹性时为何可能突然侧弯？

B. Quantitative and derivation questions

7. 均匀拉伸到 \(L/L_0=1.10\)，求工程应变与真应变。
8. 工程应力 \(300\ \mathrm{MPa}\)、工程应变 \(0.10\)，假设均匀不可压缩，求真应力。
9. Maxwell 材料 \(E=1.0\ \mathrm{GPa}\)、\(\eta=50\ \mathrm{GPa\,s}\)，求松弛时间及 \(t=\tau\) 时相对初始应力。
10. Kelvin--Voigt 材料在 \(t=2\tau\) 的蠕变应变达到长期值的多少？
11. 悬臂梁端载加倍、长度加倍而 \(EI\) 不变，端挠度变为原来的多少倍？
12. 同一柱由铰—铰改为固—固，理想 Euler 载荷变为多少倍？
13. \(E=70\ \mathrm{GPa}\)、\(I=2.0\times10^{-6}\ \mathrm{m^4}\)、\(L=2.0\ \mathrm m\) 的铰—铰柱，求 \(P_{cr}\)。
14. 裂纹几何不变时，裂纹长从 \(a\) 增至 \(4a\)，达到同一 \(K_{IC}\) 所需名义应力如何变化？

Exit ticket

观察到一根聚合物细柱在恒定压载下缓慢侧弯直至失稳。分别指出材料模型与结构模型中至少一个必须超越“\(\sigma=E\varepsilon\) + 静态 Euler 柱”的因素。

8. References and source notes / 参考资料与使用说明

8.1 Source-to-content map

来源最后核对：2026-08-20。正文与例题均为综合重写；全部示意图和函数图依据本课方程原创绘制。