A Finite Element Method for Simulation of Coupled Dynamics of Dislocations and Fracture
Boyang Gu, Adrian Diaz, Yang Li, Youping Chen
arXiv:2511.22676·cond-mat.mtrl-sci·Published 2025-11-27
This work presents a finite element method for simulating dynamic processes that involve the coupled evolution of dislocation motion and crack propagation. The method numerically solves the Concurrent Atomistic-Continuum (CAC) formulation of the conservation of linear momentum. A crystalline material is discretized at the unit-cell level using 6-node prism elements whose geometry allows dislocations and cracks to nucleate and propagate along element facets. Nanoscale simulations of single-crystal Cu, Fe, and Si demonstrate the initiation and propagation of dislocations and cracks, and these results are reproduced by the finite element method in excellent agreement with fully atomistic molecular dynamics simulations. Mesoscale simulations of single-crystal Cu further demonstrate the ability of the method to capture size-dependent brittle and ductile behavior. Under plane-strain conditions the Cu model fractures in a brittle manner, while a fully three-dimensional model exhibits curved and intersecting dislocations that blunt the crack tip and prevent crack propagation, resulting in ductile behavior. The accuracy, efficiency, and applicability of the method are discussed.
TopicsQuantum Chemistry & Force Fields
Tagsmolecular-dynamics
arXiv categoriescond-mat.mtrl-sci, physics.comp-ph
arXiv abstract pagePDF