A new method of crack propagation X-FEM reliable, robust and powerful
by D. Colombo, University of Manchester
A new method for crack growth simulation, based on the X-FEM which avoids to explicitly mesh the crack, is available in PROPA_FISS. This method, called "geometrical", is based on an implicit representation of the crack by means of two level set functions, as in the case of the "upwind" and "simplex" methods already available. However, unlike these ones, the update of the level sets consequent to the propagation is based on a geometric calculation and not on the resolution of Hamilton-Jacobi partial differential equations. The new method differs also from the "maillage" method available in PROPA_FISS, which calculates the two level set functions locating the crack by means of a mesh representing the crack surface and independent from the structural mesh.
![]() |
With this new method, crack propagation simulation gets more stabile, more robuste and much faster. In the 3D mixed mode I + II + III propagation test shown in figures 1 to 3, the geometrical method is 2 to 7 times faster than the "maillage" method. Unlike this one, the calculation time for the geometrical method is independent of crack extension and is constant throughout the whole propagation simulation. It is also 3 to 5 times faster than the "upwind" method with calculated crack fronts, and crack surfaces that are much smoother. In this case, the "simplex" method, second available method based on the integration of Hamilton-Jacobi equations to obtain the level sets evolution, fails to converge.
![]() |
The localisation of the computational domain where the the crack level sets are updated, which was introduced for the “simplex” and “upwind” methods, is not recommended with this new method because it can be time consuming, except for very large meshes. Finally, the regular auxiliary grid used to solve Hamilton-Jacobi update equations is no longer necessary, which greatly simplifies the implementation of the propagation process.
Finally cases of crack front splitting and merging during the propagation are automatically taken into account by this new method, without any intervention from the user, as for the other methods of crack propagation.
![]() |
![]() |
This work is the subject of a publication: D. Colombo. An implicit geometrical approach to level sets update for 3D non planar X-FEM crack propagation. Computer Methods in Applied Mechanics and Engineering, published online.