基于能量限制器(energy limiter)概念的动力脆性断裂非局部损伤区自适应更新技术

《COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING》：A novel nonlocal damage zone adaptation technique for dynamic brittle fracture based on the energy limiter concept

　　

摘要：本文提出一种新颖且通用的数值技术，用于在基于能量限制器(energy limiter)的梯度损伤(gradient damage)框架内更新求解域及相关的非局部扩散裂纹，旨在提高裂纹扩展分析中计算效率。鉴于脆性材料中的裂纹通常在相对于整体结构较窄的区域内萌生和演化，所提出的算法仅在活跃损伤区(active damage zone)内求解非局部损伤演化方程，该区域外的未损伤区被排除在损伤场更新之外，从而避免不必要的计算。具体而言，这种加速处理方法通过在弥散裂纹(smeared crack)推进过程中调整参与梯度损伤演化律的自由度(degree of freedom, DOF)数目，自适应地重新定义活跃求解域，从而在不损害正则化效益或非局部损伤建模框架精度的前提下显著降低计算成本。为进一步提升性能，该算法与显式动力有限元(explicit dynamic finite element)格式相结合，通过直接时间积分更新位移、速度和加速度，无需迭代求解器。二维和三维数值算例表明该方法稳健，能准确预测裂纹萌生、扩展、分叉及汇合，且与常规全域模拟相比大幅缩短求解时间。本研究为降低梯度损伤分析的计算成本、使其更适用于大规模断裂模拟开辟了新方向。

本文解读对象为 Hung Thanh Tran 发表于《COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING》的论文"A novel nonlocal damage zone adaptation technique for dynamic brittle fracture based on the energy limiter concept"。

一、研究背景与问题提出

动态脆性断裂(dynamic brittle fracture)在时变载荷下的准确高效数值模拟是计算固体力学中的经典难题。连续介质损伤力学(Continuum Damage Mechanics, CDM)通过弥散裂纹(smeared crack)表征开裂，避免了显式追踪不连续面，其中局部(local)损伤模型存在严重的网格敏感性(mesh sensitivity)和病态局部化(pathological localization)问题。经典的梯度增强损伤(gradient-enhanced damage)模型及与之理论相近的相场断裂(phase-field, PF)模型通过引入特征长度尺度使损伤在有限宽度内扩散，恢复了边值问题的适定性(well-posedness)并保证网格客观(mesh-objective)结果。Tran 等人此前提出的基于能量限制器(energy limiter)概念的隐式梯度增强损伤(implicit gradient-enhanced damage)公式，将非局部损伤场作为由非局部演化方程控制的主未知量，并在本构中嵌入能量损伤阈值(energy damage threshold)以抑制伪损伤(spurious damage)并实现窄区内准确损伤局部化，其物理可确定的特征长度尺度由总能量限制器和断裂能(fracture energy, G_f)给出，已在准静态和动力脆性断裂中验证有效性。然而该公式在每个加载步均需在全域(domain Ω)求解梯度增强损伤方程，对于大规模问题计算量极为昂贵。现有加速策略如自适应网格加密(Adaptive Mesh Refinement, AMR)、全局—局部(global-local)或多尺度耦合、区域分解(domain decomposition)或保留全域损伤方程求解、或修改离散网格、或引入复杂多级耦合与空间剖分，实现难度随三维一般裂纹增长而增大。因此，研究人员开展了以下研究：在固定网格上，利用脆性断裂裂纹萌生与扩展局限于窄带区的特点，提出仅在自适应更新的活跃域(active domain / active damage zone)内求解非局部损伤演化方程，未损伤区跳过损伤场更新，不改变网格与原梯度损伤列式，保留正则化属性和精度，并与集中质量矩阵(lumped mass matrix)的显式动力时间积分结合，形成高效计算框架，通过二维及三维典型动力断裂算例验证。研究表明该方法可在百万自由度(DOF)规模动力断裂模拟中实现高达约70%的运行时间缩减且保持精度，为大规模梯度损伤与相场断裂模拟的加速提供了简便易实现的途径。

二、主要关键技术方法概述

研究人员采用基于能量限制器(energy limiter)概念的隐式梯度增强损伤(implicit gradient-enhanced damage)本构框架，应力按σ=(1?d?)²C:ε计算(d?为非局部损伤变量，C为弹性刚度张量，ε为应变)，引入能量损伤阈值Φ_c抑制伪损伤。运动方程与损伤演化方程通过有限元离散，位移场采用一致质量或集中质量矩阵(lumped mass matrix)，时间积分采用中心差分显式格式(explicit time integration scheme)逐步更新位移、速度、加速度，免除牛顿—拉夫森(Newton-Raphson)类迭代求解。核心创新为提出非局部损伤区自适应更新技术：每时间步先按局部历史变量H(等效应变驱动的历史场)判断是否超过阈值Φ_c，若超过则标记当前损伤单元及邻域(考虑梯度扩散影响范围，通常取特征长度尺度l_c或数倍单元层)加入活跃集(active set)；仅在活跃节点子集上组集并求解非局部损伤演化方程(通常为Helmholtz型偏微分方程以获非局部损伤d?)，其余区域令d保持不变(或d=0若从未激活)。活跃域随裂纹推进动态膨胀与收缩(或维持)，可设置缓冲层(buffer layer)保证梯度算子稳定性。数值验证选用二维/三维标准动力断裂基准：混合模式裂纹扩展(mixed-mode crack propagation)、裂纹分叉(crack branching)、多裂纹扩展(multiple crack growth)及三维冲击诱发断裂(3D impact-induced fracture)，对比全求解域(full-domain)模拟的裂纹形貌、能量耗散、裂尖速度(crack-tip velocity)及CPU时间。

三、论文主体各小节结果与结论

Governing equations（控制方程）

研究人员回顾了能量限制器基梯度损伤模型的控制方程组：动量方程 ρ?v/?t = ?·σ + b (in Ω)，本构 σ = (1?d?)²C : ε，位移边界 u = ū on Γ_u，力边界 σ·n = t? on Γ_t，初始条件 u(t=0)=u₀, v(t=0)=v₀。非局部损伤变量 d? 满足含梯度的 Helmholtz 型演化/平衡方程 ?·(l_c²?d?) + d? = H(x,t)，其中 H 为由等效应变导出的历史场，当 H > Φ_c时损伤方可演化。此节确认了原有列式及显式动力时间积分流程。

Novel gradient damage zone adaptation technique（新型梯度损伤区自适应技术）

研究人员指出传统做法为一旦 max(H_n+1) > Φ_c即在全域 Ω 求解非局部损伤演化方程。由于脆性裂纹仅在窄带内局部化，研究人员设计如下流程：(1) 每步显式更新位移与局部历史场 H；(2) 扫描所有积分点，若 H > α·Φ_c(α 为略小于1的激活系数)或其相邻单元已有 d? > d_min，将该单元及 l_c影响半径内邻接单元标记为活跃；(3) 仅在活跃节点上组集损伤方程刚度阵与载荷向量，求解获得更新后的 d?；(4) 非活跃区 d? 保持原值。此算法不改变网格、不修改原梯度损伤弱形式，实现 DOF 数量随裂纹推进动态削减。研究人员说明该策略可无缝嵌入现有梯度损伤或相场代码，三维通用性强。

Numerical examples（数值算例）

研究人员设置了四个基准问题。(1) 混合模式裂纹扩展：与实验/全求解域结果对比，裂纹轨迹、裂尖速度与能量释放吻合良好，活跃域法计算时间显著减少。(2) 裂纹分叉：在高速动态加载下捕捉到分叉行为，分支角度与文献相符，损伤带宽度受 l_c控制，全域与自适应结果一致。(3) 多裂纹扩展：预置多缺陷引发多裂纹相互作用与汇合，自适应域随多裂纹前缘同步更新，预测与全求解域重合。(4) 三维冲击诱发断裂：Kalthoff 型三点弯或冲击缺口梁三维模型，显示复杂曲面裂纹面，自适应法在数百万 DOF 下获近70% CPU 时间缩减且裂纹形貌、动能—弹性能—耗散能演化匹配参考解。结果表明该自适应技术在 2D/3D 动力脆性断裂中具鲁棒性，精度无损且效率大幅提升。

Concluding remarks（结论）

研究人员总结：提出的基于能量限制器框架的动态断裂非局部损伤区自适应更新算法，通过动态演化的非局部损伤区限定损伤演化方程求解范围，在保留原模型精度与正则化特性的同时显著降低全域求解非局部损伤场的计算开销；结合显式动力求解器，在 2D 与 3D 算例中验证了其稳健性与准确性，最高可实现约70%运行时间缩减，为大尺度梯度损伤及相场断裂模拟提供实用加速手段。

四、讨论与结论翻译

结论翻译：本研究在基于能量限制器(energy limiter)的非局部损伤框架内，提出了一种用于动力断裂模拟的高效通用非局部损伤区更新算法。该方法采用动态演化的非局部损伤区，显著减少了在全域上求解非局部损伤场的计算代价，同时保留了原始模型的精度和正则化特性。结合显式动力求解器，数值结果表明该方法对二维及三维动力脆性断裂问题稳健可靠，能在保持结果准确性的前提下大幅缩短计算时间，为大规模梯度损伤与相场断裂分析提供了可行的加速途径。