主题031：损伤力学与失效分析

本章介绍了损伤力学与失效分析的基本概念、理论基础、数值方法以及应用实例。损伤力学的基本概念和发展历程。损伤变量的定义和演化。损伤本构关系，如线弹性损伤本构关系、弹塑性损伤本构关系、粘性损伤本构关系等。损伤演化方程，如Kachanov-Rabotnov模型、Lemaitre模型、Chaboche模型等。失效分析方法，如强度准则、损伤失效准则、断裂力学方法等。损伤力学的数值实现，如有限元方法、数值积分

kkchenkx

322人浏览 · 2026-03-06 23:08:21

kkchenkx · 2026-03-06 23:08:21 发布

主题031：损伤力学与失效分析

31.1 损伤力学概述

31.1.1 损伤的基本概念

损伤是指材料或结构在外部载荷、环境因素等作用下，内部产生的微观或宏观缺陷，如微裂纹、微孔洞、位错等，这些缺陷会导致材料的刚度、强度等性能逐渐降低，最终可能导致失效。

损伤的主要特征包括：

渐进性：损伤是一个逐渐累积的过程，从微观缺陷的形成到宏观失效。
不可逆性：损伤一旦产生，就无法完全恢复。
局部性：损伤通常首先在材料的薄弱区域产生，然后逐渐扩展。
耦合性：损伤与应力、应变、温度等因素相互耦合。

31.1.2 损伤力学的发展历程

损伤力学的发展可以分为以下几个阶段：

萌芽阶段（1950s-1960s）：提出了损伤的基本概念和初步模型，如Kachanov的连续性损伤模型。
发展阶段（1970s-1980s）：建立了系统的损伤力学理论框架，提出了多种损伤演化方程。
成熟阶段（1990s至今）：损伤力学与有限元方法相结合，广泛应用于工程实践，同时发展了多尺度损伤模型。

31.1.3 损伤力学的应用领域

损伤力学在工程领域有着广泛的应用，包括：

航空航天：飞机结构的损伤容限设计和寿命预测。
核工业：核反应堆压力容器的损伤分析和安全评估。
土木工程：桥梁、建筑结构的损伤演化分析。
机械工程：机械零件的疲劳损伤分析。
材料科学：新材料的损伤特性研究。

31.2 损伤变量与损伤张量

31.2.1 标量损伤变量

标量损伤变量是最简单的损伤描述方法，它用一个标量 D 来表示材料的损伤程度，其中：

D = 0 表示材料完全无损。
D = 1 表示材料完全失效。

标量损伤变量通常定义为：

$\frac{A}{A_0}$

其中，

A 是材料的有效承载面积。
A_0 是材料的初始面积。

31.2.2 张量损伤变量

张量损伤变量考虑了损伤的各向异性，它用一个二阶张量 D_{ij} 来表示材料在不同方向上的损伤程度。张量损伤变量通常定义为：

$D_{ij} = 1 - \frac{E_{ij}}{E_{ij}^0}$

其中，

$E_{ij}$ 是损伤材料的弹性张量。
$E_{ij}^0$ 是无损材料的弹性张量。

31.2.3 损伤变量的演化

损伤变量的演化是损伤力学的核心内容，它描述了损伤如何随时间、载荷等因素的变化而发展。损伤演化方程通常形式为：

$\dot{D} = f(\sigma, \varepsilon, D, T, \dots)$

其中，

$D˙\dot{D}$ 是损伤演化速率。
$σ\sigma$ 是应力张量。
$ε\varepsilon$ 是应变张量。
T 是温度。

31.3 损伤本构关系

31.3.1 线弹性损伤本构关系

线弹性损伤本构关系假设材料在损伤过程中仍然保持线弹性，其形式为：

$\sigma = (1 - D)E\varepsilon$

其中，

$σ\sigma$ 是应力。
$ε\varepsilon$ 是应变。
$E$ 是无损材料的弹性模量。
$D$ 是标量损伤变量。

对于各向异性损伤，本构关系为：

$\sigma = (I - D) : E : \varepsilon$

其中，

$I$ 是单位张量。
$D$ 是损伤张量。
$:$ 表示张量的双点积。

31.3.2 弹塑性损伤本构关系

弹塑性损伤本构关系考虑了材料的塑性变形和损伤的耦合作用，其形式为：

$\sigma = (1 - D)E(\varepsilon - \varepsilon^p)$

其中，

$εp\varepsilon^p$ 是塑性应变。

弹塑性损伤本构关系通常包括以下几个部分：

弹性关系：描述应力与弹性应变的关系。
塑性流动规则：描述塑性应变的发展。
硬化规律：描述屈服面的演化。
损伤演化方程：描述损伤的发展。

31.3.3 粘性损伤本构关系

粘性损伤本构关系考虑了材料的粘性效应，适用于高温、高应变速率等情况，其形式为：

$\sigma = (1 - D)E(\varepsilon - \varepsilon^p - \varepsilon^v)$

其中，

$εv\varepsilon^v$ 是粘性应变。

31.4 损伤演化方程

31.4.1 Kachanov-Rabotnov模型

Kachanov-Rabotnov模型是最早的连续损伤模型之一，它假设损伤演化与应力有关，其损伤演化方程为：

$\dot{D} = A(\sigma)(1 - D)^{-n}$

其中，

$A(σ)A(\sigma)$ 是与应力有关的函数。
$n$ 是材料常数。

31.4.2 Lemaitre模型

Lemaitre模型是一种基于应变能密度的损伤模型，它假设损伤演化与应变能密度有关，其损伤演化方程为：

$\dot{D} = \frac{C}{S_v}(\sigma_{eq} - \sigma_0)^m$

其中，

$C$ , $m$ 是材料常数。
$S_v$ 是材料的体积弹性模量。
$σeq\sigma_{eq}$ 是等效应力。
$σ0\sigma_0$ 是材料的初始屈服应力。

31.4.3 Chaboche模型

Chaboche模型是一种考虑循环加载的损伤模型，它假设损伤演化与塑性应变幅有关，其损伤演化方程为：

$\dot{D} = C(\Delta\varepsilon^p)^m(1 - D)^{-n}$

其中，

$C$ , $m$ , $n$ 是材料常数。
$Δεp\Delta\varepsilon^p$ 是塑性应变幅。

31.4.4 各向异性损伤模型

各向异性损伤模型考虑了损伤的方向依赖性，其损伤演化方程为：

$\dot{D}_{ij} = f(\sigma_{ij}, D_{ij})$

其中，

$D_{ij}$ 是损伤张量。
$σij\sigma_{ij}$ 是应力张量。

31.5 失效分析方法

31.5.1 强度准则

强度准则是判断材料是否失效的最简单方法，它基于材料的极限强度，当应力达到材料的极限强度时，材料发生失效。常用的强度准则包括：

最大拉应力准则：当最大拉应力达到材料的抗拉强度时，材料发生失效。
最大剪应力准则（Tresca准则）：当最大剪应力达到材料的剪切强度时，材料发生失效。
最大畸变能准则（Mises准则）：当畸变能密度达到材料的临界值时，材料发生失效。

31.5.2 损伤失效准则

损伤失效准则基于损伤变量，当损伤变量达到临界值时，材料发生失效。常用的损伤失效准则包括：

标量损伤失效准则：当标量损伤变量 D 达到临界值 $D_c$ 时，材料发生失效。
张量损伤失效准则：当损伤张量的某个不变量达到临界值时，材料发生失效。

31.5.3 断裂力学方法

断裂力学方法基于裂纹扩展理论，当裂纹扩展到临界尺寸时，材料发生失效。常用的断裂力学方法包括：

线弹性断裂力学：基于应力强度因子的断裂准则。
弹塑性断裂力学：基于J积分、COD等参量的断裂准则。

31.5.4 损伤与断裂的耦合分析

损伤与断裂的耦合分析考虑了损伤的累积和裂纹的扩展之间的相互作用，它将损伤力学和断裂力学相结合，更全面地描述材料的失效过程。

31.6 损伤力学的数值实现

31.6.1 有限元方法

有限元方法是损伤力学数值实现的主要工具，它将结构离散为有限个单元，在每个单元中计算损伤的演化，然后更新单元的刚度矩阵。

损伤力学有限元分析的主要步骤包括：

网格划分：将结构离散为有限个单元。
材料模型定义：选择合适的损伤模型，输入模型参数。
边界条件和载荷施加：定义结构的边界条件和外部载荷。
求解：在每个时间步，求解有限元方程，更新损伤变量和单元刚度矩阵。
后处理：分析结构的损伤演化和失效情况。

31.6.2 数值积分方法

损伤演化方程的求解通常需要使用数值积分方法，如欧拉法、龙格-库塔法等。对于复杂的损伤模型，可能需要使用隐式积分方法来保证数值稳定性。

31.6.3 损伤局部化问题

损伤局部化是指损伤集中在一个狭窄的区域内，形成宏观裂纹的现象。损伤局部化会导致有限元分析中的网格敏感性问题，即计算结果依赖于网格的密度和形状。

为了解决损伤局部化问题，发展了多种方法，如：

非局部损伤模型：考虑损伤的空间相关性。
梯度增强损伤模型：引入损伤梯度项。
粘性正则化：引入粘性效应。
损伤塑性模型：将损伤与塑性耦合。

31.7 Python 实现：损伤力学模型

31.7.1 标量损伤模型

import numpy as np
import matplotlib.pyplot as plt

# 设置中文字体
plt.rcParams['font.sans-serif'] = ['SimHei', 'DejaVu Sans']
plt.rcParams['axes.unicode_minus'] = False

class ScalarDamageModel:
    """标量损伤模型"""
    def __init__(self, E, sigma_0, D_c, A, n):
        """
        参数:
        E: 弹性模量 (MPa)
        sigma_0: 初始屈服应力 (MPa)
        D_c: 临界损伤值
        A: 损伤演化参数
        n: 损伤演化参数
        """
        self.E = E
        self.sigma_0 = sigma_0
        self.D_c = D_c
        self.A = A
        self.n = n
        self.D = 0  # 初始损伤
    
    def stress_strain_relation(self, strain):
        """应力-应变关系
        
        参数:
        strain: 应变
        
        返回:
        stress: 应力 (MPa)
        """
        if self.D >= self.D_c:
            return 0
        elastic_stress = self.E * strain
        damage_factor = 1 - self.D
        return damage_factor * elastic_stress
    
    def update_damage(self, stress, strain):
        """更新损伤
        
        参数:
        stress: 应力 (MPa)
        strain: 应变
        
        返回:
        D: 更新后的损伤
        """
        if self.D >= self.D_c:
            return self.D
        
        # 计算损伤演化速率
        if stress > self.sigma_0:
            dot_D = self.A * ((stress - self.sigma_0) / self.sigma_0) ** self.n
        else:
            dot_D = 0
        
        # 更新损伤
        self.D += dot_D
        
        # 限制损伤值
        if self.D > self.D_c:
            self.D = self.D_c
        
        return self.D

31.7.2 Lemaitre损伤模型

class LemaitreDamageModel:
    """Lemaitre损伤模型"""
    def __init__(self, E, nu, sigma_0, D_c, C, m):
        """
        参数:
        E: 弹性模量 (MPa)
        nu: 泊松比
        sigma_0: 初始屈服应力 (MPa)
        D_c: 临界损伤值
        C: 损伤演化参数
        m: 损伤演化参数
        """
        self.E = E
        self.nu = nu
        self.sigma_0 = sigma_0
        self.D_c = D_c
        self.C = C
        self.m = m
        self.D = 0  # 初始损伤
        self.epsilon_p = 0  # 初始塑性应变
    
    def stress_strain_relation(self, strain):
        """应力-应变关系
        
        参数:
        strain: 应变
        
        返回:
        stress: 应力 (MPa)
        """
        if self.D >= self.D_c:
            return 0
        
        # 计算弹性应变
        epsilon_e = strain - self.epsilon_p
        
        # 计算应力
        damage_factor = 1 - self.D
        stress = damage_factor * self.E * epsilon_e
        
        return stress
    
    def update(self, strain):
        """更新状态变量
        
        参数:
        strain: 应变
        
        返回:
        stress: 应力 (MPa)
        D: 损伤
        epsilon_p: 塑性应变
        """
        if self.D >= self.D_c:
            return 0, self.D, self.epsilon_p
        
        # 计算弹性应变
        epsilon_e = strain - self.epsilon_p
        
        # 计算应力
        damage_factor = 1 - self.D
        stress = damage_factor * self.E * epsilon_e
        
        # 塑性修正
        if stress > self.sigma_0:
            # 计算塑性应变增量
            d_epsilon_p = (stress - self.sigma_0) / self.E
            self.epsilon_p += d_epsilon_p
            
            # 重新计算应力
            epsilon_e = strain - self.epsilon_p
            stress = damage_factor * self.E * epsilon_e
            
            # 更新损伤
            sigma_eq = stress  # 假设为单向应力状态
            dot_D = self.C * ((sigma_eq - self.sigma_0) / self.sigma_0) ** self.m
            self.D += dot_D
            
            # 限制损伤值
            if self.D > self.D_c:
                self.D = self.D_c
        
        return stress, self.D, self.epsilon_p

31.7.3 非局部损伤模型

class NonlocalDamageModel:
    """非局部损伤模型"""
    def __init__(self, E, sigma_0, D_c, A, n, l):
        """
        参数:
        E: 弹性模量 (MPa)
        sigma_0: 初始屈服应力 (MPa)
        D_c: 临界损伤值
        A: 损伤演化参数
        n: 损伤演化参数
        l: 非局部特征长度
        """
        self.E = E
        self.sigma_0 = sigma_0
        self.D_c = D_c
        self.A = A
        self.n = n
        self.l = l
        self.D = 0  # 初始损伤
        self.D_nonlocal = 0  # 初始非局部损伤
    
    def stress_strain_relation(self, strain):
        """应力-应变关系
        
        参数:
        strain: 应变
        
        返回:
        stress: 应力 (MPa)
        """
        if self.D_nonlocal >= self.D_c:
            return 0
        elastic_stress = self.E * strain
        damage_factor = 1 - self.D_nonlocal
        return damage_factor * elastic_stress
    
    def update_damage(self, stress, strain, D_local):
        """更新损伤
        
        参数:
        stress: 应力 (MPa)
        strain: 应变
        D_local: 局部损伤
        
        返回:
        D_nonlocal: 更新后的非局部损伤
        """
        if self.D_nonlocal >= self.D_c:
            return self.D_nonlocal
        
        # 计算局部损伤演化速率
        if stress > self.sigma_0:
            dot_D_local = self.A * ((stress - self.sigma_0) / self.sigma_0) ** self.n
        else:
            dot_D_local = 0
        
        # 更新局部损伤
        self.D += dot_D_local
        
        # 计算非局部损伤（简单平均）
        self.D_nonlocal = 0.5 * (self.D + D_local)
        
        # 限制损伤值
        if self.D_nonlocal > self.D_c:
            self.D_nonlocal = self.D_c
        
        return self.D_nonlocal