综述：局部方法在结构材料裂纹扩展行为中的应用

《Materials》：Application of Local Approaches to Crack Development in Structural Materials

　　

科研人员及工程界广泛关注单个构件及复杂结构服役行为的模拟与分析。本综述重点阐述了基于有限元法(FEM)的结构材料建模领域的最新进展。除确定性描述外，潜在应力集中区前缘的状态还借助现代统计学手段加以定义，并通过长度尺度(length scales)概念解释其力学响应。构件中裂纹萌生与扩展的模拟是本研究的重点应用方向，旨在揭示微观结构对宏观性能的局域化影响机制。多尺度分析方法主要分为层级式(hierarchical)和并发式(concurrent)两大类，混合类型虽存在但不在本文讨论范围内。本文主要展示双尺度层级式模拟案例。建模对象涵盖三类材料体系：(i)夹杂物周围、晶界碳化物处以及天然或人工应力集中源附近的裂纹萌生；(ii)含短纤维或长纤维且界面特征显著的复合材料；(iii)晶粒结构表现为短纤维行为的(i)与(ii)组合体系。本文作者团队基于多年数值方法研究积累，特别聚焦于改进型有限元法的开发与应用。

1. 引言

结构构件的安全性要求必须证明其在服役过程中不会发生裂纹萌生乃至扩展，且需纳入环境因素的综合影响。机械载荷、热载荷及外部应力的多种耦合作用可能导致脆性断裂。从断裂力学角度，需重点解决两个关键科学问题：一是小试样测得的断裂韧性数据向大尺寸实际构件的转移适用性(transferability)；二是以可接受概率预测脆性断裂的发生。物理系统的建模通常依赖偏微分方程配合初边值条件，但仅在少数特定几何构形下存在解析解，一般情形下需借助数值方法逼近。数值模拟结果普遍存在误差，因此需要通过与实验数据的对比不断校准以提升精度。有限元法(FEM)的发展历史可追溯至20世纪初Ritz和Galerkin基于变分原理求解偏微分方程的早期工作，经过半个多世纪的发展已形成完备的数学理论体系。过去百余年间，弹性、塑性、断裂等相关理论层出不穷，但在缺乏合理物理驱动的线性化假设（如比例性原则和叠加原理）的工程应用中，这些理论的数学和计算统一面临巨大挑战。固体物理学、数值数学、计算机架构及并行系统的进步使得跨尺度理论构建成为可能，例如认识到位错运动与累积这一微观过程是塑性变形的本质，推动了结构化变形(structured deformation)理论的发展，该理论通过有限应变张量的乘法分解推广了可变形体弹塑性运动学描述，能够同时捕捉宏观尺度的光滑几何变化和亚宏观尺度的非光滑几何变化（即编码连续体中裂纹和缺陷存在的"紊乱"，disarrangements）。针对断裂区的分析，现有方法包括对经典Mazars模型的修正，引入损伤因子、各向异性推广、平滑梯度损伤技术、随机方法及利用分形维数量化裂纹分布等手段，这些方法均需能量正则化(energy regularization)支持。非局部公式可用于评估宏观裂纹系统的潜在萌生，其行为可基于内聚力模型(cohesive model)进行研究，并通过扩展有限元法(XFEM)或单位分解法(partition of unity methods)实现高效数值求解。传统FEM采用奇异单元表征裂尖应力奇异性，而XFEM通过引入渐近裂尖富集函数处理此类奇异性，且在结合内聚力模型时可省略渐近富集函数。相较于传统FEM在裂纹扩展过程中需要反复重新划分网格导致计算成本剧增，XFEM无需重划或仅需少量局部重划，且因网格独立于裂纹位置，避免了FEM中追踪点时间历程的复杂性。近几十年来XFEM被广泛应用于复合材料、压电材料及功能梯度材料等复杂应力响应先进材料的断裂研究，随着材料科学向智能材料和纳米技术发展，XFEM精确模拟微纳尺度断裂行为的能力将在下一代生物相容材料和柔性电子器件等新兴技术领域发挥关键作用。本文在上述背景下，聚焦近年来未被系统性综述覆盖的XFEM及其改进版本在特定分析中的应用进展。全文结构安排如下：第2节阐述计算方法所需的物理基础；第3节总结数学模型构建，详述基于XFEM的方法；第4节介绍其他数值实现途径；第5节展示作者研究团队的相关典型案例；第6节提出受XFEM启发的新技术方法及其与前人工作的关联；第7节总结成果并指出未来研究方向。

2. 物理基础

脆性断裂中裂纹萌生与扩展呈渐进式发展。多数钢材的裂纹成核发生于晶界脆性颗粒（如碳化物）处，由位错在这些颗粒处的堆积引起的应力集中所驱动；复合材料中则发生于界面处。这一现象解释了为何局部塑性变形常先于断裂发生。断裂局部方法的核心是最弱链(weakest-link)概念，即含有大量统计独立体积的材料体失效由其中一个参考体积V_ref的失效引发，该参考体积定义为包含统计显著数量微裂纹的材料体积。在有限元分析中，参考体积需保持恒定，且裂尖前缘区域单元尺寸必须小于V_ref，以确保局部应力计算的网格无关性。裂纹核(crack nuclei, CNs)若在扩展过程中未表现出失稳特征，最终会发生钝化并不再参与断裂起始过程。单位体积内CNs的生成速率由塑性应变水平决定，由于裂纹或缺口前缘局部塑性应变可能存在超过一个数量级的变化，该参数具有显著敏感性，同时还受到温度T的显著影响。从数学角度看，传统模型中认定的参考体积V_ref并非恒定值，其大小受温度和局部塑性应变的共同影响，结构钢中V_ref在裂尖前缘及晶界碳化物处的变尺寸规律已有实验观测佐证。复合材料特别是纤维复合材料的失效分析需采用不同方法。当裂纹垂直于增强纤维方向扩展时，通常涉及四种基本机制：基体开裂、纤维/基体脱粘、纤维断裂及纤维拔出。纤维与基体间的界面是决定复合材料力学及断裂行为的关键因素，其界面断裂韧性直接影响应力在增强相与基体间的传递效率。微观结构的形态特征——包括体积分数、增强相尺寸与空间分布、界面强度及缺陷尺寸——共同决定了损伤演化路径。玻璃纤维复合材料制造过程中施加的浸润剂(sizing)控制着界面性能，改变浸润剂配方会显著影响断裂韧性和强度等力学性能。复合材料设计中，裂纹萌生由强度准则控制，而裂纹扩展与损伤演化则由断裂韧性主导。I型开裂模式下裂纹沿纤维路径扩展时会发生桥接现象，这种纤维桥接区需作为独立机制建模，对复合材料分层及缺口/孔洞周围的裂纹分叉行为至关重要。陶瓷材料作为一类相对新兴的材料体系，可采用类似纤维复合材料的建模策略，其中纤维的影响被晶粒的行为所取代。氮化硅(Si₃N₄)等陶瓷在现代轴承应用中表现优异，但其服役环境往往伴随高温、腐蚀介质及循环接触应力等严苛条件。通过调控微观结构（特别是晶界相改性和晶粒桥接效应）可显著提升陶瓷的断裂韧性与强度。上述几何构型更新的物理表述需建立在经典热力学守恒原理之上，即能量守恒、质量守恒及线动量和角动量守恒（热力学第一定律）。这些原理数学上表述为偏微分方程的初边值问题，并需配合符合热力学第二定律的本构关系以实现合理的工程应用，例如热传导对应的傅里叶定律、多孔介质中污染物扩散对应的菲克定律，以及弹性变形对应的胡克定律（可进一步推广至开尔文或麦克斯韦模型描述粘弹性行为）。本文重点关注时间相关的应变-应力关系，至少需包含弹性、粘性、塑性和损伤分量，且考虑到现代计算策略（包括前述XFEM技术）普遍采用变分、弱形式或积分形式替代经典强形式微分方程。在裂尖框架下，可在平衡变量中引入内部状态变量以描述非牛顿流体及具有内部结构的材料行为。理性热力学以公理化方式处理连续介质热力学，利用客观性(objectivity)和材料对称性约束本构函数，相比标准方法对 constitutive函数的限制更为宽松。此外，材料特性设计及相关实验标定、数学与计算识别工具构成另一类重要问题。广义弹性模型往往参数过多且难以标定，导致反问题不适定(ill-posed)，即鲁棒性和唯一性缺失，微小输入扰动即可引起识别结果的不稳定响应。因此亟需探索能够显著减少参数数量并建立微-宏尺度解析与数值关联的可选表述形式。

3. 数学模型

上述材料及复合材料在机械、热及其他载荷下的应变-应力演化及损伤风险评估，需要在微观、介观、宏观等多尺度分别构建模型并实现合理的尺度衔接。多尺度方法涵盖从头算(ab initio)计算到宏观结构建模的广阔范围。基于XFEM和桥接域方法(Bridging Domain Method, BDM)的双层裂纹扩展建模简化方案展示了尺度衔接的实现路径。当前工程研究中物理、数学及计算资源尚未形成统一框架，尺度衔接收敛性的数学理论发展有其自身规律，仅部分受到工程问题的驱动。以相场模型为例，其核心问题在于当分布函数中引入的长度尺度趋于零时，描述分布式脆性断裂的泛函是否收敛于离散断裂对应的泛函。材料特性均匀化(homogenization)的数学理论同样远未统一，经典结果与新进展足以构成一篇专门综述。工程中成功应用的代表性均匀化方法包括周期双尺度收敛、G-收敛和H-收敛等，但这些严格的数学结果在结合最弱链模型和微观结构随机性的工程研究中往往被简化处理。代表性体积单元(Representative Volume Element, RVE)配合尺度衔接优化技术构成了一种合理折衷方案，可与图像处理或机器学习支持的随机技术相结合。大多数工程研究倾向于通过实验验证而非严格的数学存在性、收敛性证明来评估模型可靠性。针对不同特殊材料类别，XFEM（或类似方法）常与随机建模结合使用，包括考虑裂纹几何不确定性的形状灵敏度分析（基于比例边界有限元法）、曲率-裂纹间距的概率预测建模等。近期研究尝试将XFEM类数值分析与机器学习、人工神经网络及蒙特卡洛模拟等软计算技术融合。为便于比较不同数学与计算途径，本文构建一个简洁的确定性模型问题作为基准：考虑非极性可变形体，采用小应变线性化假设，位移分量均相对于初始零状态定义。设可变形体占据三维欧几里得空间Ω内的区域，时间区间为[0,T]。Ω的Lipschitz边界Γ分为互不相交的两部分Γ_D（位移边界）和Γ_N（力边界）。采用笛卡尔坐标系x_i，t表示时间。引入虚位移δu作为检验函数，基于勒贝格空间、索伯列夫空间和Bochner-索伯列夫空间的抽象函数空间框架为形式化存在性与收敛性证明提供支撑。运动方程基于线动量守恒，引入对称应力张量σ_ij，其与应变ε_ij的关系由本构方程确定。本文依次考察五种材料模型：(1)纯弹性（线性问题）；(2)粘弹性（仍属线性问题）；(3)弹塑性（第一类非线性推广）；(4)微观损伤（第二类非线性推广）；(5)宏观开裂（第三类非线性推广）。纯弹性情形对应胡克定律σ_ij=C_ijklε_kl，四阶刚度张量C_ijkl在完全各向同性情况下仅需拉梅常数或杨氏模量与泊松比即可完全确定，并需满足正定性条件。然而纯弹性模型构成封闭物理系统，无法描述真实工程应用中的能量耗散机制，需引入粘弹性模型。标准线性固体(Standard Linear Solid, SLS)模型作为合理折衷方案，通过引入标量正系数η和E_v构建串联弹簧-粘壶-并联弹簧的力学示意方案，其控制方程为积分-微分形式的本构关系。弹塑性推广通过在总应变中分离出塑性应变分量ε^p_ij，基于关联流动法则引入塑性势F和塑性流动参数λ，并配合屈服准则F(σ, ε^p)≤0及Karush-Kuhn-Tucker互补条件刻画加载与卸载过程。针对材料拉压性能差异（如水泥基基体），可采用Rankine型拉伸准则与Hill型压缩准则的组合。微观损伤的另一种非线性来源是通过损伤因子d∈[0,1]降低刚度，实际应力张量表示为σ=(1-d)C:ε，其中d由包含多个材料参数的复杂算法确定，且需要非局部正则化步骤以保证数学问题的良定性。需注意应力不变量而非具体应力分量被用于d的计算，以保证客观性（与坐标系选择无关）。宏观开裂则需考虑多体接触及位移跳跃，在弱形式中引入沿裂纹面的面力贡献项，对应内聚力模型中的牵引-分离律(traction-separation law)。内聚力方法仅需最少参数即可解决各类材料和载荷下的实际问题，内聚应力σ_c和内聚能G_c是核心参数，可通过实验测定，其中G_c在一阶近似下由J积分的临界值J_c确定。内聚区内材料分离律的形态可为平滑型（如指数型）或包含纤维拔出与桥接效应的复杂形态。

4. 计算实现

针对前述线性模型问题，理论及数值分析存在两种基本途径：Rothe方法和Faedo-Galerkin方法。Rothe方法首先对时间域进行离散，将弱形式下的抛物型方程组转化为一系列椭圆型方程，仍在无限维函数空间中求解，随后通过对空间域的离散（基于XFEM等技术）获得大型稀疏线性代数方程组，非线性情况需通过迭代程序处理。Faedo-Galerkin方法首先离散空间域，将原系统转化为大型稀疏常微分方程组，理论上可通过广义特征向量和特征值求解，但实践中仍需额外的时间离散化。鉴于线性模型下两种方法等价而非线性推广时存在差异，本文倾向于采用Rothe方法进行算法设计说明。将时间区间划分为N个离散时刻t_n，采用隐式（后向）欧拉格式以避免条件稳定性限制，得到时间半离散格式后再对空间域进行离散。标准线性样条及自适应FEM策略需针对XFEM的复杂结构进行修订。三种非线性贡献分别处理：弹塑性贡献在方程右端引入塑性应变相关项，需迭代更新塑性内变量；粘弹性-损伤贡献修改左端刚度矩阵和右端荷载向量，需迭代更新损伤因子；宏观开裂贡献引入位移跳跃相关项，需迭代更新裂纹构型并动态适应多裂纹系统的演化，对XFEM近似提出更高要求。多尺度材料表征——涵盖从电子密度到连续介质层级的层级结构——是当前断裂力学的前沿挑战。脆性及准脆性断裂建模通常采用水平集(level-set)方法，利用跳跃函数和渐近分支函数分别捕捉裂纹面位移间断和裂尖渐近场，富集近似形式为u = ΣN_Iu_I + ΣN_Jψ_JH + ΣN_KΦ_K，其中N_I为标准FEM节点形函数，N_J为裂纹贯穿单元节点，N_K为裂尖富集单元节点，ψ_J和Φ_K为相应富集函数。桥接域方法作为域分解耦合方法，将全域Ω分解为重叠的子域Ω_FEM和Ω_MD，前者采用FEM/XFEM，后者采用细观尺度模型（如分子动力学）。在重叠区引入权重函数α(x)实现能量的连续过渡，总能量泛函表示为E_total = αE_FEM + (1-α)E_MD + λ^T(u_FEM-u_MD)，其中λ为拉格朗日乘子强制两尺度位移协调。并发式和层级式是两类基本的跨尺度计算框架：层级式适用于尺度清晰分离的问题，信息从小尺度逐级传递至大尺度；并发式动态耦合不同尺度，对局部化现象（尺度界限模糊）的分析至关重要。物理一致性的保障依赖于Cauchy-Born规则和Hill-Mandel条件及其推广形式。不确定性量化可通过赋予外载荷、材料参数等随机变量以概率分布来实现，结合蒙特卡洛方法处理大规模稀疏非线性方程组的求解。

5. 示例案例

以金属间化合物合金Ti-46Al-0.7Cr-0.1Si-7Nb-0.2Ni（带尖锐V型缺口）为模型材料进行验证。该类金属间化合物因其优异的高温力学性能和低密度被视为最具前景的"复合"材料之一，主要应用于航空航天领域涡轮部件以替代镍基合金。采用XFEM预测裂纹扩展行为，材料本构采用指数型内聚力模型。三点及四点弯曲试验采用截面约3.2×4.1 mm、长度约45 mm的梁式试样，断裂韧性通过金刚石砂轮切割的V型缺口试样测定。系统持续计算裂尖前缘J积分值，裂纹长度由模型中单元尺寸决定（本案例单元边长为0.08 mm）。材料参数为：杨氏模量E=175,000 MPa，泊松比ν=0.3，屈服应力σ_y=300 MPa，内聚应力σ_c=550 MPa，内聚能G_c=120 N/m。通过优化算法标定所需参数后，预测的裂纹扩展阻力曲线(R曲线)与实验数据吻合良好，验证了方法的可靠性。结果表明，塑性区尺寸对裂纹扩展区域的影响与牵引-分离律精度的敏感性具有相似规律，较高塑性水平在一定程度上削弱了微观不均匀性的影响，这为统计方法的引入提供了依据。

6. 展望、未来工作与发展趋势

除经典热机械理论与计算技术外，近年涌现的多种新方法值得关注。尽管大量计算结果已通过实验验证，但针对强解、弱解、甚弱解等存在性及收敛性的形式化数学验证仍远未完善，未来进展可能需要非线性均匀化理论结合松弛微态连续介质模型或广义非平衡可逆-不可逆耦合(GENERIC)热力学框架。机器学习与人工智能技术在纤维增强混凝土破坏数值模拟等领域的应用正在兴起，但当前多数工作缺乏足够规模的经实验验证的科学数据集来构建基于物理机制的预测模型，物理守恒律（质量、动量、能量）仍以偏微分方程的一般数学形式为主导。未来技术发展可归纳为若干方向：一是拓展基于泛函分析能量表述的多层级建模，特别是在低层级采用更接近物理图景的层级式描述，配套研究算法的稳定性、收敛性及物理合理的过渡区确定；二是对新型材料模型开展改进型FEM的数学分析，涵盖静力、准静态、粘性或塑性变形及动力分析，可延伸至分子统计模型及分子动力学层面；三是推广节点数恒定的本征(intrinsic) XFEM方法，配合形函数修正和最小二乘法，解决新表面生成及裂纹扩展的技术难题；四是深化基于牵引-分离律的单层级方法在裂尖前缘局部微小区域过程机理的研究；五是面向复合材料（尤其是纤维复合材料）的两极应用——从水泥基复合材料到纤维增韧陶瓷，发展适配局部区域尺寸与整体几何差异巨大的建模策略；六是融合标准FEM/XFEM与随机过程的单/多层模型，多层建模虽算法复杂但更具潜力，桥接域方法是该类别的代表性实现；七是利用蒙特卡洛方法或其变体消除统计不确定性，通过微孔洞生长分布的统计模型作为裂纹萌生与