Three dimensional virtual crack closure-integral method (VCCM) with skewed and non-symmetric mesh arrangement at the crack front

In this paper, an extension of virtual crack closure-integral method (VCCM) for three dimensional linear fracture mechanics analysis using hexahedron finite elements is presented. In conventional three dimensional VCCM, there are some inherent requirements on the finite element model. They are (i) the faces of finite elements across the crack front have the same areas and (ii) they must be arranged symmetrically across the crack front. In present study, we developed a three dimensional VCCM without such requirements by considering work required to open one element face area whose shape is arbitrary. Though we assume the use of an ordinary 20 node serendipity element, present approach can be applied to other types of hexahedron elements.     

Finite element simulations of crack propagation in functionally graded materials under flexural loading

Cracks in stepped and continuously graded material specimens under flexural loading were investigated via finite element analysis. Calculation of mechanical energy release rates and propagation angles with crack-opening displacement correlation and the local symmetry (K_II = 0) criterion, respectively, provided results most efficiently and accurately, as compared with compliance and J-integral approaches and other deflection criteria. A routine was developed for automatic crack extension and remeshing, enabling simulation of incremental crack propagation. Effects of gradient profile and crack geometry on crack-tip stresses and crack propagation path are examined, and implications of these for optimal design of graded components against failure by fast fracture are discussed.     

On solving large strain hyperelastic problems with the natural element method

In this paper, an extension of the natural element method (NEM) is presented to solve finite deformation problems. Since NEM is a meshless method, its implementation does not require an explicit connectivity definition. Consequently, it is quite adequate to simulate large strain problems with important mesh distortions, reducing the need for remeshing and projection of results (extremely important in three‐dimensional problems). NEM has important advantages over other meshless methods, such as the interpolant character of its shape functions and the ability of exactly reproducing essential boundary conditions along convex boundaries. The α‐NEM extension generalizes this behaviour to non‐convex boundaries. A total Lagrangian formulation has been employed to solve different problems with large strains, considering hyperelastic behaviour. Several examples are presented in two and three dimensions, comparing the results with the ones of the finite element method. NEM performs better showing its important capabilities in this kind of applications. Copyright © 2004 John Wiley & Sons, Ltd.     

Non‐linear analysis of shells with arbitrary evolving cracks using XFEM

A new formulation and numerical procedures are developed for the analysis of arbitrary crack propagation in shells using the extended finite element method. The method is valid for completely non‐linear problems. Through‐the‐thickness cracks in sandwich shells are considered. An exact shell kinematics is presented, and a new enrichment of the rotation field is proposed which satisfies the director inextensibility condition. To avoid locking, an enhanced strain formulation is proposed for the 4‐node cracked shell element. A finite strain plane stress constitutive model based on the logarithmic corotational rate is employed. A cohesive zone model is introduced which embodies the special characteristics of the shell kinematics. Stress intensity factors are calculated for selected problems and crack propagation problems are solved. Copyright © 2004 John Wiley & Sons, Ltd.     

Dynamic soil–structure interaction analysis of layered unbounded media via a coupled finite element/boundary element/scaled boundary finite element model

In this paper, a coupled model based on finite element method (FEM), boundary element method (BEM) and scaled boundary FEM (SBFEM) (also referred to as the consistent infinitesimal finite element cell method) for dynamic response of 2D structures resting on layered soil media is presented. The SBFEM proposed by Wolf and Song (Finite‐element Modelling of Unbounded Media. Wiley: England, 1996) and BEM are used for modelling the dynamic response of the unbounded media (far‐field). The standard FEM is used for modelling the finite region (near‐field) and the structure. In SBFEM, which is a semi‐analytical technique, the radiation condition at infinity is satisfied exactly without requiring the fundamental solution. This method, also eliminates the need for the discretization of interfaces between different layers. In both SBFEM and BEM, the spatial dimension is decreased by one. The objective of the development of this coupled model is to combine advantages of above‐mentioned three numerical models to solve various soil–structure interaction (SSI) problems efficiently and effectively. These three methods are coupled (FE–BE–SBFEM) via substructuring method, and a computer programme is developed for the harmonic analyses of SSI systems. The coupled model is established in such a way that, depending upon the problem and far‐field properties, one can choose BEM and/or SBFEM in modelling related far‐field region(s). Thus, BEM and/or SBFEM can be used efficiently in modelling the far‐field. The proposed model is applied to investigate dynamic response of rigid and elastic structures resting on layered soil media. To assess the proposed SSI model, several problems existing in the literature are chosen and analysed. The results of the proposed model agree with the results presented in the literature for the chosen problems. The advantages of the model are demonstrated through these comparisons. Copyright © 2004 John Wiley & Sons, Ltd.     

A kernel free particle‐finite element method for hypervelocity impact simulation

An improved hybrid particle‐finite element method has been developed for the simulation of hypervelocity impact problems. Unlike alternative methods, the revised formulation computes the density without reference to any kernel or interpolation functions, for either the density or the rate of dilatation. This simplifies the state space model and leads to a significant reduction in computational cost. The improved method introduces internal energy variables as generalized co‐ordinates in a new formulation of the thermomechanical Lagrange equations. Example problems show good agreement with exact solutions in one dimension and good agreement with experimental data in a three‐dimensional simulation. Copyright © 2005 John Wiley & Sons, Ltd.     

Boundary element‐free method (BEFM) for two‐dimensional elastodynamic analysis using Laplace transform

In this paper, we present a direct meshless method of boundary integral equation (BIE), known as the boundary element‐free method (BEFM), for two‐dimensional (2D) elastodynamic problems that combines the BIE method for 2D elastodynamics in the Laplace‐transformed domain and the improved moving least‐squares (IMLS) approximation. The formulae for the BEFM for 2D elastodynamic problems are given, and the numerical procedures are also shown. The BEFM is a direct numerical method, in which the basic unknown quantities are the real solutions of the nodal variables, and the boundary conditions can be implemented directly and easily that leads to a greater computational precision. For the purpose of demonstration, some selected numerical examples are solved using the BEFM. Copyright © 2005 John Wiley & Sons, Ltd.     

Mixed finite element method for coupled thermo‐hydro‐mechanical process in poro‐elasto‐plastic media at large strains

A mixed finite element for coupled thermo‐hydro‐mechanical (THM) analysis in unsaturated porous media is proposed. Displacements, strains, the net stresses for the solid phase; pressures, pressure gradients, Darcy velocities for pore water and pore air phases; temperature, temperature gradients, the total heat flux are interpolated as independent variables. The weak form of the governing equations of coupled THM problems in porous media within the element is given on the basis of the Hu–Washizu three‐filed variational principle. The proposed mixed finite element formulation is derived. The non‐linear version of the element formulation is further derived with particular consideration of the THM constitutive model for unsaturated porous media based on the CAP model. The return mapping algorithm for the integration of the rate constitutive equation, the consistent elasto‐plastic tangent modulus matrix and the element tangent stiffness matrix are developed. For geometrical non‐linearity, the co‐rotational formulation approach is utilized. Numerical results demonstrate the capability and the performance of the proposed element in modelling progressive failure characterized by strain localization and the softening behaviours caused by thermal and chemical effects. Copyright © 2005 John Wiley & Sons, Ltd.     

Calibration of the Weibull stress scale parameter, σu, using the Master Curve

This work proposes that the Weibull stress scale parameter, σ_u, increases with temperature to reflect the increasing microscale toughness of ferritic steels caused by local events that include plastic shielding of microcracks, microcrack blunting, and microcrack arrest. The Weibull modulus, m, then characterizes the temperature invariant, random distribution of microcrack sizes in the material. Direct calibration of σ_u values at temperatures over the DBT region requires extensive sets of fracture toughness values. A more practical approach developed here utilizes the so-called Master Curve standardized in ASTM Test Method E1921-02 to provide the needed temperature vs. toughness dependence for a material using a minimum number of fracture tests conducted at one temperature. The calibration procedure then selects σ_u values that force the Weibull stress model to predict the Master Curve temperature dependence of K_Jc values for the material. At temperatures in mid-to-upper transition, the process becomes more complex as fracture test specimens undergo gradual constraint loss and the idealized conditions of high-constraint, small-scale yielding assumed in E1921-02 gradually degenerate. The paper develops the σ_u calibration process to incorporate these effects in addition to consideration of threshold toughness effects and the testing of fracture specimens with varying crack-front lengths. Initial illustrations of the calibration process for simpler conditions, i.e. 1T crack-front lengths, use the temperature dependent flow properties and a range of toughness levels for an A533B pressure vessel steel. Then using the extensive fracture toughness data sets for an A508 pressure vessel steel generated recently by Faleskog et al. [Engng. Fract. Mech., in press], the paper concludes with calibrations of both m and σ_u over the DBT region and assessments of the Master Curve calibration approach developed here.