No abstract

In this article, we propose a thermodynamically consistent phase-field model for thermo-mechanical fracture and provide an open-source implementation of the proposed model using a recently developed finite element toolbox, Gridap in Julia. Here, we have derived the balance equations for the thermo-mechanical fracture by invoking the virtual power principle and determined the constitutive relations for the thermodynamic fluxes based on the satisfaction of the thermodynamic laws. Our proposed formulation provides an equation of temperature evolution that can easily accommodate dissipative effects such as viscous damping. We provide very compact and user-friendly open-source codes for implementing the proposed model using Gridap in Julia that requires very low memory usage and gives a high degree of flexibility to the users in defining weak forms of the governing partial differential equations (PDEs). We have validated the proposed model and its implementation against such standard results available in the literature as crack propagation in the cruciform shape material, single edge notched plate, bi-material beam, and a quenching test.

This article proposes an open-source implementation of a phase-field model for brittle fracture using a recently developed finite-element toolbox, Gridap in Julia. This work exploits the advantages of both the phase-field model and Gridap toolbox for simulating fracture in brittle materials. On one hand, the use of the phase-field model, which is a continuum approach and uses a diffuse representation of sharp cracks, enables the proposed implementation to overcome such well-known drawbacks of the discrete approach for predicting complex crack paths as the need for re-meshing, enrichment of finite-element shape functions, and an explicit tracking of the crack surfaces. On the other hand, the use of Gridap makes the proposed implementation very compact and user-friendly that requires low memory usage, and provides a high degree of flexibility to the users in defining weak forms of partial differential equations. Tests on a single-edge notched plate under tension, an L-shaped panel, a notched plate with a hole, a notched beam under symmetric three-point bending and a notched beam with three holes under asymmetric three-point bending are considered to demonstrate how the proposed Gridap-based phase-field Julia code can be used to simulate fracture in brittle materials.

This article proposes an open-source implementation of a phase-field model for brittle fracture using a recently developed finite element toolbox, Gridap in Julia. The present work exploits the advantages of both the phase-field model and Gridap toolbox for simulating fracture in brittle materials. On one hand, the use of the phase-field model, which is a continuum approach and uses a diffuse representation of sharp cracks, enables the proposed implementation to overcome such well-known drawbacks of the discrete approach for predicting complex crack paths as the need for re-meshing, enrichment of finite element shape functions and an explicit tracking of the crack surfaces. On the other hand, the use of Gridap makes the proposed implementation very compact and user-friendly that requires low memory usage and provides a high degree of flexibility to the users in defining weak forms of partial differential equations. A test on a notched beam under symmetric three-point bending and a set of tests on a notched beam with three holes under asymmetric three-point bending is considered to demonstrate how the proposed Gridap based phase-field Julia code can be used to simulate fracture in brittle materials.

In this paper, we propose a phase-field model for electro-mechanical fracture in brittle materials and provide an open-source implementation of the proposed model using a newly developed finite element (FE) toolbox, Gridap in Julia. Here, we have considered electric potential as an additional kinematic descriptor along with the displacement and phase-field variable. Using the virtual power principle, we have derived force balances for the electro-mechanical forces and a force balance associated with the damage. To incorporate the strong electro-mechanical coupling effect into the model, we have considered the Helmholtz free energy as a function of the electric field vector, strain tensor, phase-field variable, and the gradient of the phase-field variable. The proposed model ensures that the constitutive relations of the thermodynamic fluxes are on the satisfaction of the thermodynamic laws and can readily accommodate dissipative energy effects whenever needed. The proposed model provides complex crack paths as a solution to the governing partial differential equations (PDEs), and thus bypasses the need for re-meshing, enrichment of FE shape functions, and an explicit tracking of the crack surfaces. The use of Gridap makes the FE implementation of the proposed model exceedingly compact and user-friendly requiring very low memory usage and providing a high degree of flexibility to the users in defining weak forms of the governing PDEs. We have validated the proposed model and its implementation against a compact tension test on a single edge notched plate and a set of three-point bending tests on a notched beam made of lead zirconium titanate (PZT)-4 piezoelectric ceramics.

No abstract

No abstract

No abstract