A comparison of stress intensity factors through crack closure integral and other approaches using eXtended element-free Galerkin method

In this paper, a new approach for extracting stress intensity factors (SIFs) by the extended element-free Galerkin method, through a crack closure integral (CCI) scheme, is proposed. The CCI calculation is used in conjunction with a local smoothing technique to improve the accuracy of the computed SIFs in a number of case studies of linear elastic fracture mechanics. The cases involve problems of mixed-mode, curved crack and thermo-mechanical loading. The SIFs by CCI, displacement and stress methods are compared with those based on the M-integral technique reported in the literature. The proposed CCI method involves very simple relations, and still gives good accuracy. The convergence of the results is also examined.

Element-Free Galerkin Modelling of Crack Migration in Composite Laminates

The development of a virtual testing environment, as a cost-effective industrial design tool in the design and analysis of composite structures, requires the need to create models efficiently, as well as accelerate the analysis by reducing the number of degrees of freedom, while still satisfying the need for accurately tracking the evolution of a debond, delamination or crack front. The eventual aim is to simulate both damage initiation and propagation in components with realistic geometrical features, where crack propagation paths are not trivial. Meshless approaches, and the Element-Free Galerkin (EFG) method, are particularly suitable for problems involving changes in topology and have been successfully applied to simulate damage in homogeneous materials and concrete. In this work, the method is utilized to model initiation and mixed-mode propagation of cracks in composite laminates, and to simulate experimentally-observed crack migration which is difficult to model using standard finite element analysis. N

The Behavior of compressively loaded stiffener runout specimens - Part II: Finite element analysis

The termination of stiffeners in composite aircraft structures give rise to regions of high interlaminar shear and peel stresses as the load in the stiffener is diffused into the skin. This is of particular concern in co-cured composite stiffened structures where there is a relatively low resistance to through-thickness stress components at the skin-stiffener interface. In Part I, experimental results of tested specimens highlighted the influence of local design parameters on their structural response. Indeed some of the observed behavior was unexpected. There is a need to be able to analyse a range of changes in geometry rapidly to allow the analysis to form an integral part of the structural design process.

This work presents the development of a finite element methodology for modelling the failure process of these critical regions. An efficient thick shell element formulation is presented and this element is used in conjuction with the Virtual Crack Closure Technique (VCCT) to predict the crack growth characteristics of the modelled specimens. Three specimens were modelled and the qualitative aspects of crack growth were captured successfully. The shortcomings in the quantitative correlation between the predicted and observed failure loads are discussed. There was evidence to suggest that high through-thickness compressive stresses enhanced the fracture toughness in these critical regions.

Modified crack closure integral technique for extraction of SIFs in meshfee methods

A new approach for extracting stress intensity factors (SIFs) by the element-free Galerkin (EFG) class of methods through a modified crack closure integral (MCCI) scheme is proposed. Its primary feature is that it allows accurate calculation of mode I and mode II SIFs with a relatively simple and straightforward analysis even when a coarser nodal density is employed. The details of the adoption of the MCCI technique in the EFG method are described. Its performance is demonstrated through a number of case studies including mixed-mode and thermal problems in linear elastic fracture mechanics (LEFM). The results are compared with published theoretical solutions and those based on the displacement method, stress method, crack closure integral in conjunction with local smoothing (CCI–LS) technique, as well as the M-integral method. Its advantages are discussed.

Finite element modelling of the formation of pre-load damage in cement mantles of orthopaedic joint replacements

Finite element modeling of the formation of pre-loaded damage in cement mantles of orthopaedic joint replacements was presented. The existence of cracking suggested a high level of residual stress. The direction of maximum principal stress vectors corresponded well with the observed crack orientation. Results suggested that cracking depends upon a combination of residual stress, porosity and temperature rise during polymerization.

Prediction of failure and fracture mechanisms of polymeric composites using finite element analysis. Part 1:Particulate filled composites

Micro-mechanical analysis of polymeric composites provides a powerful means for the quantitative assessment of their bulk behavior. In this paper we describe a robust finite element model (FEM) for the micro-structural modeling of the behavior of particulate filled polymer composites under external loads. The developed model is applied to simulate stress distribution in polymer composites containing particulate fillers. Quantitative information about the magnitude and location of maximum stress concentrations obtained from these simulations is used to predict the dominant failure and crack growth mechanisms in these composites. The model predictions are compared with the available experimental data and also with the values found using other methods reported in the literature. These comparisons show the range of the validity of the developed model and its predictive potential.

Prediction of failure and fracture mechanisms of polymeric composites using finite element analysis. Part 2:Fiber reinforced composites

A robust finite element scheme for the micro-mechanical modeling of the behavior of fiber reinforced polymeric composites under external loads is developed. The developed model is used to simulate stress distribution throughout the composite domain and to identify the locations where maximum stress concentrations occur. This information is used as a guide to predict dominant failure and crack growth mechanisms in fiber reinforced composites. The differences between continuous fibers, which are susceptible to unidirectional transverse fracture, and short fibers have been demonstrated. To assess the validity and range of applicability of the developed scheme, numerical results obtained by the model are compared with the available experimental data and also with the values found using other methods reported in the literature. These comparisons show that the present finite element scheme can generate meaningful results in the analysis of fiber reinforced composites.

Crack propagation in non-homogenous materials: Evaluation of mixed-mode SIFs, T-stress and kinking angle using a variant EFG method

A new variant of the Element-Free Galerkin (EFG) method, that combines the diffraction method, to characterize the crack tip solution, and the Heaviside enrichment function for representing discontinuity due to a crack, has been used to model crack propagation through non-homogenous materials. In the case of interface crack propagation, the kink angle is predicted by applying the maximum tangential principal stress (MTPS) criterion in conjunction with consideration of the energy release rate (ERR). The MTPS criterion is applied to the crack tip stress field described by both the stress intensity factor (SIF) and the T-stress, which are extracted using the interaction integral method. The proposed EFG method has been developed and applied for 2D case studies involving a crack in an orthotropic material, crack along an interface and a crack terminating at a bi-material interface, under mechanical or thermal loading; this is done to demonstrate the advantages and efficiency of the proposed methodology. The computed SIFs, T-stress and the predicted interface crack kink angles are compared with existing results in the literature and are found to be in good agreement. An example of crack growth through a particle-reinforced composite materials, which may involve crack meandering around the particle, is reported.