从基础理论到大变形有限元分析,全面掌握真实应力应变的物理意义、数学推导及工程应用价值。
在材料力学和金属塑性加工领域,真实应力应变(True Stress and True Strain)是描述材料在大变形阶段力学行为最准确的物理量。与传统的工程应力应变不同,它考虑了材料在受力过程中几何尺寸(截面积和长度)的实时变化。
真实应力(σ_T)定义为瞬时载荷除以试样的瞬时横截面积。
真实应变(ε_T)定义为在变形过程中每一瞬间的微小应变增量之和(即对数应变)。
当金属材料受到拉力作用时,随着变形的增加,其横截面会逐渐缩小,长度会显著增加。如果在计算应力时仍使用初始截面积,将严重低估材料内部的实际承载能力。特别是在进入塑性阶段后期,这种误差可能高达数倍。因此,理解并正确使用真实应力应变公式对于航空航天、汽车碰撞安全及精密制造至关重要。
为了精确描述材料的塑性流动规律,我们需要建立严格的数学模型。以下推导基于体积不变假设(Volume Constancy),这是金属塑性变形的一个基本前提。
设 F 为瞬时载荷,A 为瞬时截面积,则真实应力 σ_T 为:
F / A
其中,若已知初始截面积 A₀ 和初始长度 L₀,根据体积不变原理 A·L = A₀·L₀,可得 A = A₀·L₀/L。代入上式可推导出与工程应力的关系:
σ_T = σ_E (1 + ε_E)
此处 σ_E 为工程应力,ε_E 为工程应变。该公式仅在均匀变形阶段(颈缩发生前)严格成立。
工程应变 ε_E = (L - L₀) / L₀,它假设整个变形过程发生在同一尺度上,这在微小变形下适用,但在大变形下会产生累积误差。
真实应变 ε_T 采用对数定义,累加每一微段的变形:
ε_T = ∫(dL / L) = ln(L / L₀)
同样利用体积不变原理,可以将其转换为截面积的函数:
ε_T = ln(A₀ / A)
这一特性使得真实应变具有可加性,非常适合描述多道次轧制或复杂成型工艺中的总变形量。
在塑性变形阶段,大多数金属材料的真实应力应变关系遵循幂律硬化规律,这是有限元模拟中最常用的本构模型之一:
σ_T = K (ε_T)^n
其中:
K 为强度系数(Strength Coefficient),代表材料在单位真实应变下的应力水平。
n 为应变硬化指数(Strain Hardening Exponent),表征材料抵抗进一步塑性变形的能力。n 值越大,材料均匀变形能力越强。
许多初学者容易混淆这两个概念。下表详细列出了它们在定义、适用范围及物理意义上的本质区别。
| 比较维度 | 工程应力应变 (Engineering) | 真实应力应变 (True) |
|---|---|---|
| 定义基准 | 基于初始尺寸 (A₀, L₀) | 基于瞬时尺寸 (A, L) |
| 计算公式 | σ = F/A₀ ε = ΔL/L₀ |
σ_T = F/A ε_T = ln(L/L₀) |
| 颈缩后表现 | 应力下降(因名义面积不变,载荷减小) | 应力持续上升(因实际面积急剧减小) |
| 物理意义 | 设计校核用(偏保守) | 材料本构研究、失效分析用 |
| 可加性 | 不可加 | 可加(对数性质) |
在拉伸试验的应力 - 应变曲线上,工程曲线会在达到抗拉强度后开始下降,最终断裂。然而,真实应力应变曲线在颈缩发生后,由于截面积减小的速率远大于载荷下降的速率,其应力值会继续单调上升,直到材料断裂。这反映了材料内部位错运动导致的加工硬化效应并未停止。
真实应力应变公式不仅仅是教科书上的理论,它在现代工业生产中有着广泛的应用。我们将通过选项卡形式展示三个典型场景。
在汽车正面碰撞中,车身结构件(如纵梁、吸能盒)会发生巨大的塑性变形以吸收动能。此时,材料已进入大变形阶段,工程应力数据完全失效。
在深冲(Deep Drawing)或拉延工艺中,板料经历复杂的拉压组合变形。模具间隙的设定直接依赖于材料的硬化行为。
主流 CAE 软件(ABAQUS, LS-DYNA, ANSYS)在进行显式动力学或隐式静态分析时,要求用户输入真实的塑性数据。
获取准确的真实应力应变数据并非易事,需要严谨的实验流程和数据处理步骤。以下是标准化的操作流程:
选取符合 ASTM E8 或 GB/T 228 标准的圆柱或板状试样。精确测量原始直径 d₀ 和标距 L₀,建议使用千分尺在不同位置测量取平均值,以减少误差。
在标距范围内安装高精度引伸计。注意:在接近屈服点时,需小心切换或断开引伸计以防损坏,或使用非接触式光学引伸计(Digital Image Correlation, DIC)全程跟踪。
启动万能试验机,以恒定速度加载。采集系统需同步记录载荷 F 和位移 ΔL。采样频率建议在屈服阶段高于 100Hz,均匀变形阶段不低于 10Hz。
利用专用软件或 Excel 宏,将 F-ΔL 数据转换为工程应力 - 应变曲线。随后,针对塑性段应用真实应力应变公式进行转换。对于颈缩后的数据,需结合最小截面法或 DIC 技术修正截面积。
使用幂律方程 σ = Kεⁿ 对转换后的数据进行回归拟合,提取 K 值和 n 值,用于后续的材料建模。
A: 在弹性变形阶段,材料的变形量极小(通常小于 0.2%),截面积的变化微乎其微。此时 (1+ε) ≈ 1,工程应力与真实应力的数值差异在工程允许误差范围内(<1%),因此为了简化计算,两者通常视为相等。
A: 如果缺乏引伸计,只能依靠十字头位移,但这包含了机器柔度和夹具滑移误差,数据不可靠。一种折衷方案是使用光栅引伸计或非接触式视频引伸计。若仅有载荷 - 位移数据,且必须估算,可利用体积不变假设结合断口最小直径反推,但这仅适用于断裂后的估算。
A: 理论上,随着颈缩加剧,真实应力趋向于无穷大。但实际上,材料内部会出现微孔洞聚集(Micro-void coalescence),导致有效承载面积急剧下降,应力无法无限增加。断裂时的真实应力通常对应于材料的极限抗拉强度之后的某个高点,具体取决于材料的韧性。