论文标题
一种用于脆性断裂相模型的平滑粒子流体动力学方法
A smoothed particle hydrodynamics approach for phase field modeling of brittle fracture
论文作者
论文摘要
裂缝是工程和物理中各种应用的非常具有挑战性且复杂的问题。尽管它在基于网格的数值技术(例如有限元方法(FEM))的背景下进行了广泛的研究,但平滑粒子流体动力学(SPH)社区中的研究活动仍然很少。 SPH是一种基于粒子的数值方法,用于离散连续介质的方程。它的无网状性质使模拟涉及极端变形的断裂场景是理想的选择。但是,为了建模骨折,SPH研究人员主要依赖于临时经验的局部损害模型,凝聚区方法或伪弹力模型,这些模型带有一组缺点和局限性。另一方面,脆性断裂的相位田地模型最近在学术界越来越受欢迎,与以前的方法相比提供了重大改进。这些改进包括从基本断裂理论,非本地性的引入以及在没有可用裂纹路径的先验知识的情况下,对多重裂纹启动,传播,分支和聚结的能力进行了模拟的能力。然而,在SPH中尚未研究骨折的相位场。在本概念验证论文中,我们在SPH的背景下开发并实施了脆性断裂的相位场模型。提供了全面的数学和实施细节,并计算了几个具有挑战性的数值示例,并说明了所提出的方法准确有效地模拟复杂的断裂方案的能力。
Fracture is a very challenging and complicated problem with various applications in engineering and physics. Although it has been extensively studied within the context of mesh-based numerical techniques, such as the finite element method (FEM), the research activity within the Smoothed Particle Hydrodynamics (SPH) community remains scarce. SPH is a particle-based numerical method used to discretize equations of continuum media. Its meshfree nature makes it ideal to simulate fracture scenarios that involve extreme deformations. However, to model fracture, SPH researchers have mostly relied on ad-hoc empirical local damage models, cohesive zone approaches, or pseudo-spring models, which come with a set of drawbacks and limitations. On the other hand, phase field models of brittle fracture have recently gained popularity in academic circles and provide significant improvements compared to previous approaches. These improvements include the derivation from fundamental fracture theories, the introduction of non-locality, and the ability to model multiple crack initiation, propagation, branching, and coalescence, in situations where no prior knowledge of the crack paths is available. Nevertheless, phase field for fracture has not been studied within SPH. In this proof-of-concept paper we develop and implement a phase field model of brittle fracture within the context of SPH. Comprehensive mathematical and implementation details are provided, and several challenging numerical examples are computed and illustrate the proposed method's ability to accurately and efficiently simulate complex fracture scenarios.