航空学报 > 2022, Vol. 43 Issue (3): 225159-225159   doi: 10.7527/S1000-6893.2022.25159

基于子结构法与损伤识别的周期性结构脆性断裂相场模拟

王博臣1, 侯玉亮2, 夏凉1, 史铁林1   

  1. 1. 华中科技大学 数字制造装备与技术国家重点实验室, 武汉 430074;
    2. 郑州大学 机械与动力工程学院, 郑州 450001
  • 收稿日期:2020-12-25 修回日期:2021-04-08 出版日期:2022-03-15 发布日期:2022-04-02
  • 通讯作者: 夏凉 E-mail:xialiang@hust.edu.cn
  • 基金资助:
    国家自然科学基金(11972166);湖北省杰出青年基金(2020CFA080);中央高校基本科研业务费专项资金(2019kfyXKJC044)

Phase field modeling of brittle fracture of periodic structures based on ubstructuring and damage identification

WANG Bochen1, HOU Yuliang2, XIA Liang1, SHI Tielin1   

  1. 1. State Key Laboratory of Digital Manufacturing Equipment and Technology, Huazhong University of Science and Technology, Wuhan 430074, China;
    2. School of Mechanical and Power Engineering, Zhengzhou University, Zhengzhou 450001, China
  • Received:2020-12-25 Revised:2021-04-08 Online:2022-03-15 Published:2022-04-02
  • Supported by:
    National Natural Science Foundation of China (11972166);Natural Science Foundation for Distinguished Young Scholars of Hubei Province of China (2020CFA080);Fundamental Research Funds for the Central Universities (2019kfyXKJC044)

摘要: 相场断裂模型将裂纹近似描述为全局弥散标量场,通过求解偏微分方程模拟裂纹的成核与扩展,避免了不连续位移场中对裂纹几何形状描述和裂纹尖端追踪的难题;在处理裂纹分叉、多裂纹、非均质材料和三维等复杂断裂模拟问题时具有显著的优势。然而,相场断裂模型需要在裂纹附近进行精细化的有限元网格剖分,计算规模庞大,因此,对计算能力提出了极高的挑战。基于子结构法和损伤识别,发展了一种适用于周期性结构脆性断裂的高效相场模拟方法。首先,根据周期性特征对结构进行分区并对子域内部节点进行子结构静态凝聚,降低有限元模型位移响应的维数;然后,根据离线(Off-line)子域典型测试工况标定弹性应变能阈值,在线(On-line)断裂模拟过程中仅对阈值以上子域的内部位移和断裂相场进行解耦分析。一方面,利用结构的周期性特征,通过子结构静态凝聚降低计算规模;另一方面,根据能量阈值损伤识别,避免非裂纹扩展子域的重复计算消耗。

关键词: 相场, 断裂模型, 计算缩减, 子结构建模, 损伤识别

Abstract: Phase field models for fracturing approximate the crack as a globally diffused scalar field.The crack nucleation and propagation are simulated by solving partial differential equations, which avoids the numerical challenges in describing and tracking cracks in the discontinuous displacement field, and exhibits several advantages in modeling conventionally challenging problems, such as crack bifurcation, multiple cracks, fracture problems within heterogeneous materials and three-dimensional cases.However, fracture modeling based on the phase field method requires a refined finite element mesh near the crack, resulting in enormous computational costs.With substructuring and damage identification, this work proposes a highly efficient phase field method for brittle fracturing simulation of periodic structures.Firstly, the structure is partitioned into identical substructures according to its periodicity, and their internal nodes are statically condensed by means of substructuring.Secondly, an elastic strain energy threshold is identified by representative off-line tests of the substructure.Moreover, only the internal displacement and crack phase field of the substructures whose strain energies are above the threshold are further analyzed during the on-line fracturing simulation.Therefore, not only the structural periodicity contributes to the reduction of computational costs through sbustructuring, but the non-crack propagation substructures are prevented from repetitive FE analysis by virtue of damage identification.

Key words: phase field, fracture modeling, computation reduction, substructuring, damage identification

中图分类号: