2. 连续介质力学基础回顾

在材料损伤断裂有限元分析中，所有数值结果最终都落在“应力—应变—位移”这三类场量上。若对连续介质力学的基本方程没有清晰认识，后续看到损伤变量、断裂能、单元删除或裂纹扩展结果时，很容易只记住操作而无法判断结果是否合理。本节把有限元离散之前必须掌握的力学基础串成一条逻辑链：如何描述变形、如何描述内力、二者通过什么关系联系、整体必须满足什么平衡，以及如何把连续问题转化为可求解的代数问题。

一、应力与应变的度量

描述一点附近的变形，最基本的是位移场 u(x)。小变形假设下，应变张量由位移梯度对称化得到：

ε_ij = 1/2 ( ∂u_i/∂x_j + ∂u_j/∂x_i )

其中 ε_xx、ε_yy、ε_zz 为线应变，γ_xy = 2ε_xy 为工程剪应变。小变形理论适用于多数金属结构在损伤萌生前的阶段；当裂纹尖端附近出现大转动或大应变时，应改用有限变形度量，如格林—拉格朗日应变 E = 1/2(FᵀF − I)，其中 F 为变形梯度。

应力度量同样需要区分。柯西应力 σ 定义在当前构形上，是真实应力；第一类皮奥拉—基尔霍夫应力 P 把力映射到参考构形，不对称；第二类皮奥拉—基尔霍夫应力 S 与格林应变功共轭，对称且便于本构积分。在损伤力学中，若损伤变量定义在参考构形上，通常采用 S 与 E 的功共轭对；若在商业软件中查看结果，则多数输出的是柯西应力。参照约定必须全文统一：本文默认小变形下使用 σ 与 ε，大变形下明确标注所用度量。

二、平衡方程与几何方程

平衡方程来自动量守恒。静力条件下，域内每一点满足：

∂σ_ij/∂x_j + f_i = 0        （域内）
σ_ij n_j = t_i               （边界上）

其中 f_i 为体积力，n_j 为外法向，t_i 为面力。平衡方程是强形式，要求应力场处处可微；有限元并不直接求解它，而是转化为等效积分弱形式。

几何方程给出应变与位移的关系，即上文 ε_ij 表达式。它保证变形协调，即由应变积分回位移时不出现矛盾。在有限元中，位移是节点自由度，应变由位移插值函数求导得到，因此几何方程天然满足协调性；但若使用减缩积分，可能出现沙漏模式，需要额外控制。

三、本构关系

本构关系连接应力与应变，是材料行为的核心。线弹性下：

σ_ij = C_ijkl ε_kl

各向同性材料只需弹性模量 E 与泊松比 ν。进入损伤阶段后，常用连续损伤力学框架，引入损伤变量 D（0 表示完好，1 表示完全失效），有效应力 σ̃ = σ/(1−D)，本构变为：

σ = (1 − D) C : ε

损伤演化通常由等效塑性应变、最大主应变或断裂能释放率驱动。例如延性损伤可写为 D = (ε_p − ε_p^0)/(ε_p^f − ε_p^0)，其中 ε_p^0 为损伤起始等效塑性应变，ε_p^f 为失效等效塑性应变。这些参数必须与单元特征长度关联，否则网格敏感。

四、虚功原理与有限元离散

虚功原理把平衡方程与力边界条件转化为积分形式：对任意虚位移 δu，有

∫_V σ : δε dV = ∫_V f · δu dV + ∫_S t · δu dS

左端为内虚功，右端为外虚功。该式不要求应力场二阶可微，降低了对解的光滑性要求，是有限元法的出发点。

离散时，把位移场用节点位移 d 与形函数 N 插值：u = N d，应变 ε = B d，其中 B 为应变—位移矩阵。代入虚功原理并消去任意虚位移，得到单元平衡方程：

K_e d_e = f_e^ext − f_e^int
K_e = ∫_V Bᵀ C B dV
f_e^int = ∫_V Bᵀ σ dV

组装后得到整体方程 K d = F。线性问题一次求解；损伤断裂问题中 K 随 D 演化，需增量迭代，每步更新损伤、判断单元是否失效，并重新组装刚度。

五、损伤断裂分析中的关键衔接

把上述基础用于损伤断裂，需要特别注意三点。第一，本构中的 D 使切线刚度不断退化，Newton 迭代可能不收敛，应控制增量步长并采用线搜索或粘性正则化。第二，单元失效后若直接删除，会造成刚度矩阵奇异或能量不守恒，工程中常用断裂能 G_f 与特征长度 l_e 定义等效失效应变 ε_f = 2G_f/(σ_y l_e)，使能量耗散与网格尺寸弱相关。第三，裂纹尖端应力奇异，位移型有限元需很细网格；若关注裂纹路径，可结合内聚力单元或相场法，把不连续问题正则化。

本节给出的方程与离散关系，是后续所有损伤模型、断裂准则和求解设置的共同基础。理解 B 矩阵、K 组装和 f_int 的含义，才能在软件中正确设置材料参数、控制网格并解释收敛失败的原因。