1 屈服准则的基本概念

屈服准则(yield criterion)是材料力学与本构建模中的判据,用于在给定应力状态下判断材料是否仍能保持以弹性为主的可逆响应。一旦满足屈服条件,材料的主要变形机制将转向塑性流动,外载继续增加时,变形将不可完全恢复,并伴随塑性应变累积与相关力学量的重分配。

1.1 屈服与弹塑性分界的意义

在弹性范围内,应力与应变之间通常满足线性或预定的可逆关系;当外载使材料达到临界状态,内部微观滑移、位错运动或微结构重排等机制开始主导宏观响应。屈服准则的价值在于:它把连续的真实材料行为用可计算的“边界条件”加以抽象,使工程设计有限元分析能够明确“从哪里开始算塑性”。

需要强调的是,“屈服点”在实验上往往表现为一个区间而非绝对瞬点;模型中屈服准则给出的是理想边界,实际材料的渐进过渡可通过硬化规律与塑性流动框架来体现。

1.2 判据的输入与输出:应力状态到屈服条件

屈服准则通常以应力张量(或主应力)、以及可能的内部变量(如等效塑性应变、硬化变量、温度等)作为输入。输出不是“材料是否会断裂”的结论,而是一个用于判断屈服与否的标量条件。常见做法是构造屈服函数 \(f\),并以

  • \(f<0\):弹性区(或“未达到屈服”)
  • \(f=0\):屈服发生的边界
  • \(f&gt;0\):理论上不允许的状态(数值算法会通过回算使其满足边界条件)

这样,屈服判定就可以嵌入到增量式本构更新中。

1.3 三维应力下的表征:主应力与等效应力

在三维问题中,应力状态用主应力(\(\sigma_1,\sigma_2,\sigma_3\))更直观。许多经典准则可写成主应力的组合形式;当准则强调“与应力方向无关但与应力偏量/能量等效相关”时,也会引入等效应力(如 von Mises 等效应力)。等效应力将多维应力状态映射为单一尺度,从而把复杂判定简化为对某个标量与临界值的比较。

1.4 工程尺度的“屈服点”与模型参数的关系

工程上常用的屈服应力来自单轴试验或特定加载路径,但在多轴状态下,材料是否屈服取决于所选准则的形式及其参数。不同准则对主应力差、剪应力贡献、体积应力或应力偏量的权重不同,因此同一材料的“屈服点”在不同准则下对应的参数含义会不同。通常需要通过多种实验(单轴、双轴/多轴、剪切或压缩试验)对参数进行标定,使模型能在预期应力路径上给出合理预测。


2 常用屈服准则(经典模型)

经典屈服准则多围绕“剪切引发屈服”或“应变能等效”等思想展开。它们在工程软件与有限元求解中应用广泛,但适用范围与参数可迁移性取决于材料类别(如金属、岩土、摩擦材料)与加载条件。

2.1 Tresca 最大剪应力准则

2.1.1 主应力差与临界剪应力

Tresca 准则把屈服与最大剪应力联系起来。对主应力 \(\sigma_1\ge \sigma_2\ge \sigma_3\),最大剪应力与主应力差有关。Tresca 准则常表述为:当某一对主应力的差达到与材料临界剪应力等价的阈值时发生屈服。该思想直观:屈服由最“强剪切”驱动,而非单纯由法向拉压大小决定。

从几何角度看,Tresca 在主应力空间形成分段平面边界,屈服面呈现棱边结构,因而对主应力差的敏感性较强。

2.2 von Mises 等效应力准则

2.2.1 六面体应力空间与圆柱屈服面的直观

von Mises 准则也常被称为“J2 屈服准则”。它将屈服与偏应力(与体积变形无关的那部分)联系起来。几何直观上,在偏应力空间中,屈服面常表现为圆柱形结构(沿静水压力方向延伸),意味着纯静水压力本身不直接触发屈服。该性质在许多延性金属中体现较好,因此该准则在工程金属塑性建模中非常常见。

2.2.2 等效应力公式与工程应力换算

von Mises 准则通常以等效应力 \(\sigma_\text{eq}\) 来表达。其关键是将 \(\sigma_1,\sigma_2,\sigma_3\) 组合成一个标量,使屈服条件可写为 \[ \sigma_\text{eq} = \sigma_\text{y} \] 其中 \(\sigma_\text{y}\) 为屈服应力(或由硬化关系确定的随塑性演化临界值)。在工程参数标定中,单轴屈服应力与 von Mises 临界值之间常通过等效应力定义建立对应关系,从而实现不同加载路径下的预测一致性

2.3 Mohr-Coulomb 与 Drucker-Prager(岩土/摩擦材料语境

岩土、混凝土某些阶段性响应、土壤或摩擦主导材料常呈现“抗压能力与抗拉/抗剪行为不对称”的特征。相较于金属的等向延性模型,这类材料的屈服与强度包络更贴近剪切破坏与摩擦-黏聚机制,因此常采用 Mohr-Coulomb 或 Drucker-Prager 形式。

2.3.1 适用条件与偏差来源

Mohr-Coulomb 准则基于摩尔圆与剪切破坏思想,通常采用黏聚力与内摩擦角参数来描述强度包络;它在特定应力状态下更贴近实验曲线,但在三维应力空间的严格映射与数值稳定性方面,需要额外处理。Drucker-Prager 准则作为对 Mohr-Coulomb 的三维光滑近似(思想上更便于数值计算),其偏差来源常与具体材料破坏包络的形状不完全一致有关,也与是否存在显著的非关联流动、应变硬化/软化等影响有关。

2.3.2 参数标定与边界形式

这类模型的参数通常来自三轴压缩、直剪、双轴试验或间接强度指标。标定时需明确屈服面(或强度包络)的具体形式:例如在某些实现中,参数可能通过“与某些特定应力路径的强度匹配”来确定。边界形式的选择会影响预测的峰值强度与屈服后演化趋势,因此需要与材料的实际加载工况一致。


3 屈服准则的数学表述与几何解释

屈服准则的通用数学框架是用屈服函数 \(f\) 构造边界,并与弹性区/塑性区的判别规则相联系。几何解释帮助理解屈服面在应力空间中的形状及其对不同应力路径的敏感性。

3.1 屈服函数形式:f(σ, κ)= 0

常见写法为 \[ f(\sigma,\kappa)=0 \] 其中 \(\sigma\) 表示应力张量(或等效主应力组合),\(\kappa\) 表示内部变量集合(如硬化变量、温度相关参数或等效塑性应变)。\(\kappa\) 用于描述屈服面随塑性发展如何演化,使材料从初始屈服到后续硬化阶段的边界迁移成为可计算过程。

3.2 屈服面、弹性区与塑性区

在应力空间中,满足 \(f(\sigma,\kappa)=0\) 的集合构成屈服面。对于给定 \(\kappa\),屈服面将空间分成若干区域。弹性区通常被定义为 \(f&lt;0\),塑性区不允许“超越”屈服边界(在理想塑性框架中,若达到边界则需要塑性流动来把响应约束回 \(f=0\))。在数值实现中,这种回算过程决定了塑性积分算法的正确性。

3.3 主应力空间中的屈服边界

不同准则对应不同的屈服面形状。Tresca 的屈服边界由主应力差触发,形成棱角或平面片段;von Mises 则更平滑,体现出对主应力组合的均衡考虑。对于各向异性材料,屈服面会偏离理想对称形状,可能在不同方向呈现不等强度,从而导致屈服边界随材料取向变化。

3.4 屈服条件的等价变换与不变性要求

为了使模型具有物理一致性,屈服准则通常要求满足某些不变性条件,例如在坐标旋转下的形式保持(客观应力描述)以及与材料对称性相匹配的特征。例如各向同性材料的 von Mises 类准则在应力旋转下应保持等效性;而若涉及压力敏感或各向异性参数,则不变性要求会相应调整。数值实现中还需注意等价变换对计算精度的影响,避免因公式不稳定导致误判。


4 塑性理论中的联动机

屈服准则只是“何时进入塑性”的判据;塑性理论进一步规定“进入塑性后如何演化”。联动机制主要包括流动法则(塑性应变增量与应力关系)、硬化/演化规律(屈服面如何移动或改变)、以及与破坏判据的层次区分。

4.1 流动法则(关联/非关联)

4.1.1 法向量与塑性应变增量方向

连续介质塑性中,塑性应变增量方向常由屈服面法向确定。关联流动意味着塑性势函数与屈服函数相同或成严格对应关系,因此塑性流动与屈服面法向一致。非关联流动则采用不同的塑性势,使得体积应变变化与剪切塑性发展之间的关系更符合某些材料实验现象,例如岩土材料常出现体积变化与剪切方向不完全一致的情况。

4.1.2 选取差异对预测的影响

关联/非关联的选择直接影响体积应变演化、屈服后应力路径以及数值收敛性。在某些情况下,非关联流动能更好描述剪切膨胀或压密趋势,但也会导致模型需要更多参数或更复杂标定。对于金属延性材料,多数情况下采用关联流动或接近关联的形式往往更合适。

4.2 硬化模型与屈服准则的“移动”

硬化决定屈服边界随塑性发展如何变化。若不引入硬化,屈服面在首次达到后将保持不变,无法描述材料随塑性增加的强度提升或软化。

4.2.1 等向硬化与运动硬化

等向硬化常见于各向同性延性描述:屈服面整体“膨胀”或“收缩”,中心位置不变;运动硬化则考虑屈服面在应力空间中平移,使得循环加载下的应力反向更合理。二者都可与相应屈服准则耦合,用于建模循环稳定性和卸载再加载行为。

4.2.2 含应变硬化变量的耦合形式

含应变硬化变量时,屈服函数中的内部变量 \(\kappa\) 会随塑性累积更新。耦合形式通常体现为:等效塑性应变控制临界值的演化;或以更复杂的张量形式描述硬化方向性。耦合越复杂,需要的实验数据与标定工作也越多,且更依赖选取合理的加载路径。

4.3 破坏与屈服的区分(材料响应层次)

屈服描述的是弹塑性边界;破坏或断裂通常属于更高层次的失效判据。两者并非必然同步:某些材料可能在屈服后经历一段稳定塑性流动才断裂;另一些材料可能出现屈服后迅速软化并导致局部化失稳。

4.3.1 屈服后失稳/断裂的不同判据

屈服判据与断裂准则往往独立:屈服用于确定材料进入塑性积分阶段;断裂准则则用于预测承载丧失、裂纹形成或强度衰减。断裂准则可基于损伤累积、应变能或局部应变阈值等思想。将两者混用会导致预测偏差,尤其在软化材料或大变形分析中更为明显。


5 参数确定与实验标定

屈服准则的参数并非凭空给定。由于不同准则对不同应力分量的敏感性不同,参数必须与所选材料类别和加载条件匹配。实验标定通常目标是让模型在多种加载路径下给出一致或可接受的预测。

5.1 单轴拉伸试验与屈服参数

单轴拉伸试验可直接提供弹性范围与屈服起点信息,如屈服应力、弹性模量以及一定阶段的应变硬化曲线。对于与等效应力相关的准则(尤其是 von Mises 类),单轴屈服常可用于确定初始临界值或与硬化变量的起点对应关系。但单轴数据不足以充分确定模型在复杂三轴状态下的预测能力,仍需补充多轴或剪切试验。

5.2 多轴加载试验:张压组合与剪切

多轴试验通过改变主应力比例、应力偏量与体积效应来提供更丰富的信息。例如拉-压组合能反映压力敏感性与应力路径依赖;剪切试验则对最大剪切驱动的准则更关键。通过多组实验,可以对准则参数进行约束,使屈服面形状与硬化演化更贴近真实材料。

5.3 反演标定思路与误差来源

反演标定通常把模型计算结果与实验响应进行对比,通过优化或拟合确定参数。误差来源包括:实验测量噪声、应变测量方法(如是否能准确反映局部应变)、样品尺寸效应、以及模型假设与真实材料微观机制之间的差异。数值上还可能出现过拟合或参数耦合导致的不唯一性,因此需要选择合适的目标函数与约束条件,并进行交叉验证。

5.4 温度与应变率对参数的影响(概念性)

温度与应变率会改变材料流动应力与硬化行为,导致“同一屈服准则形式”下的参数本身随工况变化。概念上,工程实践常通过引入温度/应变率修正项或采用率相关本构模型来处理。标定时需要确保实验覆盖目标温区和速率范围,否则模型外推到极端条件可能失效。


6 工程与计算应用:有限元中的屈服检测

有限元中,屈服准则用于在每个时间增量的积分点处进行“弹性-塑性切换”判断,并驱动随后的塑性本构更新。它既是物理判据,也是数值算法的一部分。

6.1 屈服判定在应力更新中的位置

典型流程是:已知上一增量的应力与状态变量,依据应变增量计算试探应力(弹性预测)。随后计算试探应力下的屈服函数 \(f\)。若试探态落在弹性区,则保持弹性更新;若越过屈服面,则需要引入塑性校正(返回到 \(f=0\) 的一致性条件),并更新塑性应变与内部变量。

6.2 塑性积分与回算流程概览

塑性积分通常以“返算”方式实现:在达到屈服边界后,通过确定塑性增量(或等效塑性乘子)使得屈服函数满足一致性条件。回算过程中还需要更新弹塑性刚度贡献及状态变量(如硬化变量)。不同算法(显式/隐式、径向回归/通用牛顿迭代等)会影响精度与稳定性,但共同目标是确保屈服判定与状态更新一致。

6.3 收敛性与网格/步长对结果的影响

屈服检测与塑性回算往往是非线性环节。时间步长过大或网格过粗会放大局部化误差,使得屈服边界频繁穿越或迭代不稳定,从而影响收敛性与预测精度。工程上常通过自适应步长、合理网格加密、以及改进的数值稳定策略来缓解问题。对于软化或强局部化材料,除了屈服判定还需考虑网格客观性(例如采用非局部或正则化方法)。

6.4 常见实现细节:一致切线与算法稳定性

隐式有限元通常要求提供一致切线刚度(consistent tangent),以改善牛顿迭代收敛。若切线不一致,可能出现收敛变慢、震荡甚至失败。实现中还需处理屈服边界附近的数值容差、主应力排序(以避免不连续)、以及屈服函数的梯度计算。对一些复杂模型,稳健的正则化与边界处理同样重要。


7 适用性、局限与改进方向

屈服准则在许多应用中有效,但其理想化假设不可避免地带来局限。改进方向通常围绕各向异性、材料微观机制、以及更贴近真实渐进行为的描述展开。

7.1 各向同性与各向异性差异

各向同性模型假设材料属性与方向无关,屈服面具有对称性;各向异性模型允许屈服面随材料取向变化,常用于轧制、拉拔等导致的方向性材料行为。各向异性模型通常需要更多参数,并对实验设计提出更高要求,否则会出现可辨识性不足或预测不稳。

7.2 材料微观机制与模型偏离

屈服准则多是宏观等效描述。真实材料屈服可能来自多种微观机制的竞争与协同,例如位错滑移、孪生、损伤萌生等。模型偏离的常见表现包括:在某些应力路径下预测过于乐观或过于保守,或难以复现实验中的非线性拐点。通过引入更丰富的内部变量、硬化规则或流动法则,通常能改善拟合,但会增加参数标定成本。

7.3 大塑性、循环载荷与鼠标点击式“梗”:反复判定导致的工程误读(轻度梗)

在大塑性或强循环载荷条件下,积分点可能频繁在屈服边界附近来回判断,尤其当步长或容差设置不理想时,会出现“屈服判定反复发生”的数值现象。有人形容这种过程像在软件里反复点击鼠标确认弹窗:看似每次都“通过”,但整体工程意义却可能被误读。更严谨的做法是通过改进步长控制、一致切线、以及合理的正则化策略,减少边界抖动,并同时关注整体响应而非单点瞬时判据。

7.4 面向更复杂材料体系的扩展思路

针对更复杂的材料体系,扩展方向包括:引入非局部或梯度增强以处理局部化;结合损伤-塑性耦合实现从屈服到软化的统一描述;使用更通用的各向异性硬化与温湿度耦合模型。选择扩展路线时需在“模型复杂度—可标定性—计算成本”之间权衡。


8 相关概念与术语

8.1 等效应力、等效应变与硬化变量

等效应力把三维应力状态映射为单一尺度;等效应变用于度量塑性累积程度。硬化变量通常代表材料内部结构演化的抽象量,决定屈服面随塑性发展如何变化。二者共同构成屈服判定与塑性积分之间的桥梁。

8.2 端点响应:屈服点、屈服应力与平台现象

屈服点是材料在实验中从弹性主导转向塑性主导的起始特征;屈服应力对应这一转变的定量描述。某些材料在屈服后可能出现平台或较弱硬化段,反映内部结构重排与宏观流动应力变化规律。模型中可通过硬化曲线形状或硬化变量演化来复现这类端点响应。

8.3 屈服准则 vs 屈服准则参数(区分要点)

屈服准则指“形式”(例如 von Mises 或 Tresca 的数学结构与几何含义);屈服准则参数指“具体数值”(如临界应力、硬化系数、摩擦角等)。不同准则形式选择会改变应力状态到屈服判定的映射关系,而参数则调节该映射与实验结果的匹配程度。二者应分别理解与标定。

8.4 与本构模型、损伤/断裂准则的关系

屈服准则通常嵌入弹塑性本构模型中,决定塑性阶段启动与状态更新;损伤与断裂准则则用于描述从屈服之后到失效之前的演化乃至断裂发生。二者在层次上可理解为:屈服是“进入塑性”的门槛,损伤/断裂是“承载能力衰减或失效”的判据。合理分工能提升预测的可解释性与可靠性。