蠕变和粘塑性是一回事吗？

在固体力学和材料科学中，粘塑性（Viscoplasticity）和蠕变（Creep）都是描述材料率相关（Time-dependent / Rate-dependent）非弹性变形的核心概念。它们经常被放在一起讨论，甚至在某些语境下会有所重叠，但它们的物理关注点、宏观表现和建模侧重点有着明显的区别。

一、 粘塑性与蠕变的核心区别

我们可以从定义、驱动力、应力阈值以及应用场景四个维度来精确区分这两种力学行为：

1. 概念定义不同

• 蠕变 (Creep)：

   指材料在恒定应力（或恒定载荷）下，非弹性应变随时间持续增长的现象。它通常是一个相对缓慢的演化过程，且往往在高温环境下（通常 T > 0.3 ~ 0.5 T_m，其中 T_m 为绝对熔点）更为显著。  

• 粘塑性 (Viscoplasticity)：

   是指材料的塑性变形不仅取决于应力的大小，还强烈取决于应变率（变形速度）。它是“粘性（液体特征）”与“塑性（固体特征）”的结合，强调的是应力-应变关系随加载速率的变化。  
   
  
   

2. 驱动条件与现象侧重

• 蠕变的视角：

   σ = Constant → ε(t) 随时间增加。典型的蠕变曲线通常包含初始蠕变、稳态蠕变和加速蠕变三个阶段。  
   
  
   

• 粘塑性的视角：

   dε/dt = Constant → σ(ε) 发生变化。加载得越快，材料表现出的屈服强度和流动应力通常越高（即率敏感性）。  
   
  
   

3. 屈服极限（Yield Surface）的概念

• 经典蠕变：

   传统的蠕变理论（如金属的高温蠕变）往往不引入明确的屈服面概念。即使应力很低（低于宏观屈服强度），只要时间够长、温度够高，蠕变就会发生。  

• 粘塑性：

   通常保留了屈服面的概念，并引入了“过应力（Overstress）”的想法。当应力低于屈服极限时，材料表现为弹性或粘弹性；一旦应力超过屈服极限，材料就会发生粘塑性流动，且超越屈服面的程度（过应力）决定了塑性应变率的大小。  

统一的视角： 在现代本构理论中，力学界普遍认为蠕变是粘塑性的一种特例（即在恒定应力边界条件下的粘塑性响应）。因此，许多先进的本构模型直接统称为统一粘塑性模型（Unified Viscoplastic Models），用一套统一的方程同时描述蠕变、应力松弛以及速率相关塑性。

二、 常见的蠕变材料模型

蠕变模型主要关注在恒定应力下，蠕变应变率 dε_c/dt 与应力 σ、时间 t 或温度 T 的数学关系。

1. 经验与半经验模型（常用于工程有限元分析）

• 幂律蠕变模型 (Power-Law / Norton-Bailey Creep):

  
   

  这是工业界最经典的稳态蠕变模型。    

  dε_c/dt = A σⁿ

  （其中 A 和 n 是材料常数。如果考虑时间硬化，可扩展为 ε_c = A σⁿ t^m）  

• 双曲正弦模型 (Garofalo / Hyperbolic Sine Law):

  
   

  适用于更广泛的应力范围，能同时模拟低应力下的线性蠕变和高应力下的指数级蠕变。    

  dε_c/dt = A [ sinh(B σ) ]ⁿ

• Andrade 蠕变模型:

  
   

  专门用来描述早期非稳态的瞬态蠕变（Transient Creep），应变与时间的 1/3 次方成正比。  

2. 基于损伤机制的模型

• Kachanov-Rabotnov 蠕变损伤模型:

  
   

  引入了一个连续损伤变量 ω (0 ≤ ω ≤ 1)，用以模拟材料由于内部微空洞演化导致的第三阶段（加速）蠕变直至断裂。  

三、 常见的粘塑性材料模型

粘塑性模型需要定量描述流变应力、应变、应变率以及硬化之间的强耦合关系。

1. 经典过应力模型 (Overstress Models)

• Perzyna 粘塑性模型:

  
   

  最著名的早期模型。它认为粘塑性应变率 dε_vp/dt 与“当前应力超过静屈服面的程度（过应力）”成正比。    

  dε_vp/dt = γ ⟨ ( f(σ) - σ_y ) / σ_y ⟩ⁿ

  （其中 ⟨•⟩ 为 MacCauley 括号，即负数取 0，正数取自身；σ_y 为静屈服应力）  

• Duvaut-Lions 粘塑性模型:

  
   

  另一种形式的过应力模型，它基于当前非率相关（弹性/纯塑性）预测解与实际粘塑性解之间的差距来计算率相关流动，数学边界性质较好，常用于现代大变形有限元算法中。  

2. 经验/半经验现象学模型（常用于冲击、切削与高速变形）

• Johnson-Cook 模型 (J-C Model):

  
   

  工业界（如 ABAQUS、LS-DYNA）中最常用的模型，将应变硬化、应变率强化和温度软化进行了解耦相乘：    

  σ = ( A + B ε_pⁿ ) ( 1 + C ln(dε^*/dt) ) ( 1 - T^*m )

• Khan-Huang-Liang 模型 (KHL Model):

  
   

  修正了 J-C 模型在高应变率下不够准确的缺点，能更好地描述强耦合的率敏感性和工作硬化行为。  

3. 统一粘塑性本构模型（Unified Models - 高级工程应用）

这些模型不再人为区分“纯塑性”和“蠕变”，而是用统一的内部变量（如随体硬化背应力 X 和各向同性硬化变量 R）来描述所有非弹性变形，非常适合模拟航空发动机、核电结构在循环热机械载荷下的蠕变-疲劳交互作用。

• Chaboche 粘塑性模型:

  
   

  基于连续介质热力学框架，引入了非线性随体硬化（多个背应力叠加）和率敏感流动法则，是目前模拟金属材料复杂循环加载、应力松弛和蠕变最强大的模型之一。  

• Bodner-Partom 模型:

  
   

  一种无屈服面（Elastic-Viscoplastic without yield surface）的统一模型。它认为非弹性变形在任何应力下都存在，只是在低应力时微乎其微，无需跟踪复杂的屈服面演化。  

四、 总结对比