概念与定义
von Mises 准则(亦称 J2 等效应力屈服判据、von Mises 屈服准则)是一类广泛使用的屈服判定理论,用于描述金属材料在复杂三维应力条件下从弹性或近似弹性行为进入塑性变形的倾向。其基本出发点是:屈服主要与“偏应力”有关,而与纯体积应力(球应力)无关或关系较弱。因而可以将三维应力状态通过一个标量量——von Mises 等效应力(与 J2 相关)——进行归并,再与材料屈服强度进行比较。
在工程计算中,该准则常用于有限元分析中的屈服判定、塑性区识别与安全系数评估。它特别适用于各向同性金属的工程估算场景,尤其在塑性屈服区对比试验数据时表现较稳定。同时,它也是与其他常见准则(如 Tresca 准则)进行对照研究与实践选择的重要基准。
1.1 屈服准则的基本思想
屈服准则的任务,是把“材料何时开始出现不可逆变形”的条件,从复杂的应力状态简化为可计算的判据。von Mises 准则采取的思路是:将应力分解为球应力与偏应力,并认为屈服由偏应力主导。于是,不同空间方向的复杂应力,只要它们对应的偏应力强度相同,就会预测相似的屈服倾向。
这种归并带来的直接好处是,三维应力张量(包含多个独立分量)可以被映射到一个等效标量上,便于在数值方法中逐点计算并做阈值比较。
1.2 von Mises 准则的核心表述
von Mises 准则通常表述为:当 von Mises 等效应力达到材料的屈服强度时,材料进入屈服状态。若采用等效应力记号 \(\sigma_\text{vM}\),屈服判据可写作
\[ \sigma_\text{vM} \ge \sigma_y \]
其中 \(\sigma_y\) 为材料屈服强度。对有限元实现而言,这一判断会在每个单元积分点进行,从而得到屈服区域分布。
1.3 J2(第二不变量)与等效应力
von Mises 等效应力与应力张量的第二不变量 \(J_2\) 密切相关。\(J_2\) 衡量偏应力分量的“能量型强度”,其定义基于偏应力张量与自身的二次型组合。通过 \(J_2\) 的开方并乘以比例系数,可以得到与实验屈服行为相对应的等效应力。
因此,在工程实现中常见流程是:先从应力张量求出偏应力,再计算 \(J_2\),最后换算为 \(\sigma_\text{vM}\),再与 \(\sigma_y\) 比较。
1.4 与偏应力(deviatoric stress)的关系
偏应力(deviatoric stress)是指把球应力去除后的那部分应力,其张量的迹为零。von Mises 准则强调:屈服与偏应力的强度有关,而球应力对屈服的影响不在该准则的主导假设之内。直观理解是,单纯改变材料的“整体体积受压/受拉”并不等同于触发屈服;真正促使塑性起始的更像是剪切形变与畸变程度的增强。
这种“偏应力控制屈服”的假设,使得 von Mises 准则在各向同性金属中具有较强的实用性。
数学表述:从应力张量到等效应力
von Mises 等效应力的计算依赖于应力张量的分解与不变量。由于判据最终基于不变量,它具有坐标变换下的客观性:无论你用哪套坐标系表示同一物理应力状态,只要偏应力相同,等效应力结果应一致。
2.1 应力不变量与坐标无关性
应力不变量是由应力张量构造、对坐标旋转不发生变化的标量量。它们使得屈服判据可以避免依赖特定坐标轴取向,从而更适合工程中“任意受力姿态”的分析。
2.1.1 第一不变量与球应力
第一不变量通常与应力张量的迹有关。设应力张量为 \(\boldsymbol{\sigma}\),则
\[ I_1 = \text{tr}(\boldsymbol{\sigma}) = \sigma_{xx}+\sigma_{yy}+\sigma_{zz} \]
球应力(平均应力)与第一不变量直接对应。球应力反映“体积平均”的作用水平,von Mises 准则在计算中对其不设主要作用。
2.1.2 第二不变量 J2 的定义
第二不变量 \(J_2\) 在 von Mises 理论中用于表征偏应力强度。若记偏应力张量为 \(\boldsymbol{s}\),则
\[ J_2 = \frac{1}{2}\,\boldsymbol{s}:\boldsymbol{s} \]
其中 “\(:\)” 表示双点积,即将张量对应分量相乘并求和。该构造体现了“与偏应力二次强度相关”的思想。
2.2 偏应力张量的构造
偏应力张量由应力张量减去球应力得到。球应力为 \(\frac{1}{3}I_1\boldsymbol{I}\),因此
\[ \boldsymbol{s} = \boldsymbol{\sigma} - \frac{1}{3}I_1\boldsymbol{I} \]
其中 \(\boldsymbol{I}\) 为单位张量。偏应力张量满足 \(\text{tr}(\boldsymbol{s})=0\),因此它专注于畸变部分。
2.3 von Mises 等效应力公式
von Mises 等效应力从 \(J_2\) 或偏应力二次型推得。经典形式为
\[ \sigma_\text{vM} = \sqrt{3J_2} \]
该表达式将“偏应力二次强度”映射到与单轴屈服相匹配的标量尺度。
2.3.1 三维主应力形式
当以主应力表示时,设三主应力为 \(\sigma_1,\sigma_2,\sigma_3\),von Mises 等效应力可写为
\[ \sigma_\text{vM} = \sqrt{\frac{(\sigma_1-\sigma_2)^2+(\sigma_2-\sigma_3)^2+(\sigma_3-\sigma_1)^2}{2}} \]
该形式强调:等效应力由主应力之间的差异决定,反映剪切/畸变相关的特征。
2.3.2 应力分量(工程计算)形式
在以应力分量进行工程计算时,可通过偏应力分量或直接使用应力张量分量的等价表达获得结果。以常见的工程实现思路为例,会先计算偏应力分量,再利用 \(J_2=\frac{1}{2}\boldsymbol{s}:\boldsymbol{s}\) 得到等效应力。
在有限元软件的输出或二次后处理里,通常会直接给出 von Mises 应力或提供可基于应力分量重建的计算口径。无论采用哪种写法,关键在于保持与所定义的 \(\sigma_\text{vM}=\sqrt{3J_2}\) 口径一致。
2.4 屈服判据的比较形式
除“等效应力与屈服强度比较”的常用形式外,也可采用“屈服函数”的表达。令
\[ f = \sigma_\text{vM} - \sigma_y \]
则屈服判定等价于 \(f \ge 0\)。在数值算法中,屈服函数常用于迭代求解或塑性流动/硬化模型的判据触发;而在工程评估中,则常以一次计算的最大值或分布图形式呈现。
几何与物理直观
von Mises 准则不仅是计算公式,也可以用“屈服面”与“主应力空间”来理解。其几何图像有助于把复杂应力状态直观化,并解释为什么它更偏向“偏应力控制”的机制。
3.1 屈服面的“距离”理解
在应力空间中,给定等效应力的集合可形成屈服面的几何形状。由于 von Mises 等效应力由 \(\sqrt{3J_2}\) 给出,屈服面常被理解为与偏应力强度相关的等值面。在这种视角下,材料是否屈服取决于当前“偏应力强度”是否越过由 \(\sigma_y\) 定义的边界。
“距离”的比喻并不等同于空间欧氏几何的字面距离,但可以用来表达:等效应力越大,越接近或超过屈服边界。
3.2 主应力空间中的表达
在主应力空间里,不同应力状态对应不同主应力组合。由于 von Mises 等效应力与主应力差的平方和相关,屈服趋势主要由主应力之间的差异所驱动。若三主应力相等,则对应纯球应力状态,主应力差为零,von Mises 等效应力也随之为零,体现了“纯体积平均不触发屈服”的直观假设。
3.3 典型应力状态的对应关系
以下对应关系用于帮助理解不同典型加载下等效应力与屈服判据的含义。实际工程中材料、几何与边界条件的组合可能更复杂,但这些例子体现了 von Mises 准则的核心偏应力逻辑。
3.3.1 纯拉伸
在单轴拉伸条件下,主应力中只有一个分量占主导。由于等效应力与主应力差有关,该状态将给出与材料单轴屈服强度相匹配的等效应力水平。换言之,von Mises 准则在此类加载下能直接复现以 \(\sigma_y\) 作为屈服阈值的工程标定思想。
3.3.2 纯剪切
在纯剪切中,主应力一般表现为成对的正负对称关系。由于主应力差显著,von Mises 等效应力会反映出剪切引起的畸变强度。其结果常用于解释为什么金属在一定剪切水平下会产生塑性流动,并把剪切效应归并到等效拉伸强度的框架中。
3.3.3 双轴拉伸/压缩的判读
双轴拉伸或双轴压缩会改变主应力之间的差异结构。若两方向的拉压关系导致偏应力强度增加,则 von Mises 等效应力随之增大,从而可能达到屈服阈值。需要注意的是:双轴状态不仅涉及“幅值”,也涉及“符号组合”与主应力差的变化,因此不能仅凭直观“谁更大”直接判断,仍应通过等效应力计算或基于屈服面理解来判读。
工程应用与计算流程
在工程软件或自编程序中,von Mises 准则的使用通常围绕“应力→偏应力→J2→等效应力→屈服判断”这一链条完成。其可靠性依赖于材料参数口径与应力结果的质量。
4.1 单元应力到等效应力的映射(如有限元)
在有限元分析中,求得每个单元积分点处的应力张量 \(\boldsymbol{\sigma}\) 后,后处理即可计算 von Mises 等效应力。典型映射步骤为:
- 由 \(\boldsymbol{\sigma}\) 计算应力迹 \(I_1\);
- 得到偏应力张量 \(\boldsymbol{s}=\boldsymbol{\sigma}-\frac{1}{3}I_1\boldsymbol{I}\);
- 计算 \(J_2=\frac{1}{2}\boldsymbol{s}:\boldsymbol{s}\);
- 得到 \(\sigma_\text{vM}=\sqrt{3J_2}\);
- 与材料屈服强度 \(\sigma_y\) 对比并形成判定结果或云图。
由于屈服过程往往发生在局部,应力的离散方式(如积分点数量、插值阶次)会影响等效应力峰值,因此需要结合网格与结果稳定性进行解释。
4.2 屈服判定与安全系数
屈服判定通常使用“等效应力是否达到阈值”的准则。若希望给出安全系数,可采用
\[ n = \frac{\sigma_y}{\sigma_\text{vM}} \]
当 \(n\ge 1\) 时认为满足不屈服条件;当 \(n<1\) 则表示可能发生屈服。工程实践中还会结合最危险点的取值(如最大 \(\sigma_\text{vM}\))或考虑一定的统计/平均策略,来形成更稳健的评估结论。
4.3 材料参数:屈服强度的选取
von Mises 准则本身需要一个屈服阈值 \(\sigma_y\)。该参数应来自相应工况下的材料数据,如单轴拉伸屈服强度或规范规定的等效屈服指标。由于不同材料体系与数据来源可能存在口径差异,参数选取的关键在于“一致性”:计算中采用的等效应力定义与屈服参数应匹配。
4.3.1 单轴拉伸屈服条件
最常见的数据基础是单轴拉伸试验的屈服强度。将该值作为 \(\sigma_y\) 后,von Mises 准则可以在多轴状态下对塑性起始进行工程估算。其原因在于该准则在逻辑上把偏应力强度标定到与单轴屈服等效的尺度。
4.3.2 压缩/拉伸数据差异的处理思路
一些材料可能表现出拉伸与压缩屈服强度不一致。von Mises 准则属于各向同性且“偏应力控制”的模型,若仅使用单一屈服强度参数,可能忽略该差异。工程上常见的处理思路是:根据材料试验获得对应的屈服强度口径,选择更保守或更符合工况的等效屈服指标;或在更精细的本构/屈服模型中引入拉压不对称的描述。无论采取何种做法,核心仍是保证参数与模型假设的一致性。
4.4 常见工程注意事项
von Mises 准则易于实现,但其工程结果的可信度受前处理、后处理与边界条件影响。注意事项主要集中在量纲、数值误差与适用范围。
4.4.1 单位与量纲一致性
计算中等效应力与屈服强度必须使用一致的单位体系。例如应力以 MPa 表示时,\(\sigma_y\) 也应以 MPa 给出。若单位混用(如应力以 Pa、屈服以 MPa),不仅会改变屈服判断,还可能导致安全系数失真。
4.4.2 网格与应力离散误差对结果的影响
屈服通常发生在应力集中的位置,例如孔洞边缘、几何突变与接触区域。有限元网格不足或局部加密不合理时,\(\sigma_\text{vM}\) 的峰值可能被低估或高估,从而影响屈服位置判断。工程实践中常通过网格收敛性检查、局部加密以及对峰值结果的敏感性分析来降低不确定性。
4.4.3 线弹性范围内的合理性边界
von Mises 准则用于判断屈服起始点时依赖“弹性应力场的可得性”。在已经明显进入大范围塑性后的场景,仅用屈服判定而不进行塑性本构更新,通常难以得到真实的响应演化。因此,在使用该准则进行早期评估或相对对比时,应把它视为“屈服触发与趋势判断”的工具;若进入非线性阶段,仍需配合弹塑性分析与合适的屈服/硬化模型。