Łukasz Kaczmarczyk

Three-dimensional brittle fracture: configurational-force-driven crack propagation

SUMMARY This paper presents a computational framework for quasi-static brittle fracture in three-dimensional solids. The paper sets out the theoretical basis for determining the initiation and direction of propagating cracks based on the concept of configurational mechanics, consistent with Griffiths theory. Resolution of the propagating crack by the FEM is achieved by restricting cracks to element faces and adapting the mesh to align it with the predicted crack direction. A local mesh improvement procedure is developed to maximise mesh quality in order to improve both accuracy and solution robustness and to remove the influence of the initial mesh on the direction of propagating cracks. An arc-length control technique is derived to enable the dissipative load path to be traced. A hierarchical hp-refinement strategy is implemented in order to improve both the approximation of displacements and crack geometry. The performance of this modelling approach is demonstrated on two numerical examples that qualitatively illustrate its ability to predict complex crack paths. All problems are three-dimensional, including a torsion problem that results in the accurate prediction of a doubly-curved crack. Copyright

A micromechanics-enhanced finite element formulation for modelling heterogeneous materials

Abstract In the analysis of composite materials with heterogeneous microstructures, full resolution of the heterogeneities using classical numerical approaches can be computationally prohibitive. This paper presents a micromechanics-enhanced finite element formulation that accurately captures the mechanical behaviour of heterogeneous materials in a computationally efficient manner. The strategy exploits analytical solutions derived by Eshelby for ellipsoidal inclusions in order to determine the mechanical perturbation fields as a result of the underlying heterogeneities. Approximation functions for these perturbation fields are then incorporated into a finite element formulation to augment those of the macroscopic fields. A significant feature of this approach is that the finite element mesh does not explicitly resolve the heterogeneities and that no additional degrees of freedom are introduced. In this paper, Hybrid-Trefftz stress finite elements are utilised and performance of the proposed formulation is demonstrated with numerical examples. The method is restricted here to elastic particulate composites with ellipsoidal inclusions but it has been designed to be extensible to a wider class of materials comprising arbitrary shaped inclusions.

Efficient numerical analysis of bone remodelling.

This paper presents a formulation for the three-dimensional numerical simulation of mechanically regulated bone adaptation. Attention is focussed on a phenomenologically-based approach to bone remodelling that can be used as a computationally efficient tool to provide insight into the overall response of bone to mechanical loading. A discretisation approach is developed based on a hybrid finite element formulation where displacement, stress and density fields are approximated independently. The paper also discusses a solution algorithm tailored for shared memory multi-core computers. The performance of the model is demonstrated by two numerical examples.

Energy consistent framework for continuously evolving 3D crack propagation

This paper presents an enhanced theoretical formulation and associated computational framework for brittle fracture in elastic solids within the context of configurational mechanics, building on the authors’ previous paper, Kaczmarczyk et al. (2014). The local form of the first law of thermodynamics provides an equilibrium condition for the crack front, expressed in terms of the configurational forces. Applying the principle of maximal energy dissipation, it is shown that the direction of the crack propagation is given by the direction of the configurational forces. In combination with a fracture criterion, these are utilised to determine the position of the continuously evolving crack front. This exploitation of the crack front equilibrium condition leads to a completely new, implicit, crack propagation formulation. A monolithic solution strategy is adopted, solving simultaneously for both the material displacements (i.e. crack extension) and the spatial displacements. The resulting crack path is resolved as a discrete displacement discontinuity, where the material displacements of the nodes on the crack front change continuously, without the need for mesh splitting or the use of enrichment techniques. In order to trace the dissipative loading path, an arc-length procedure is adopted that controls the incremental crack area growth. In order to maintain mesh quality, smoothing of the mesh is undertaken as a continuous process, together with face flipping, node merging and edge splitting where necessary. Hierarchical basis functions of arbitrary polynomial order are adopted to increase the order of approximation without the need to change the finite element mesh. Performance of the formulation is demonstrated by means of three representative numerical simulations, demonstrating both accuracy and robustness.