Finite strain analysis of limestone / basaltic magma interaction and fracture: Low order mixed tetrahedron and remeshing

In this investigation, we use a recent constitutive framework and remeshing technique for tetrahedra to analyze the pressure-driven crack propagation of limestone intruded by basaltic magma. Limestone is represented by an elasto-plastic capped Drucker-Prager model with an hypoelastic term in order to account for inelastic effects from plastic signature. Kinematic hardening is considered for limestone, whereas magma is modeled by means of a compressible Bingham fluid. Classical limit surfaces of the capped model are used to initiate the crack events in the limestone. Propagation is performed by a local remeshing technique with mesh smoothing for uncracked elements. Arnold's MINI element is used to avoid locking in the quasi-incompressible case. Verification (for limestone) examples and a crack propagation example in 3D are performed. A mesh convergence study is performed, exhibiting very promising results.