Equivalent localization element for crack band approach to mesh‐sensitivity in microplane model

The paper presents in detail a novel method for finite element analysis of materials undergoing strain‐softening damage based on the crack band concept. The method allows applying complex material models, such as the microplane model for concrete or rock, in finite element calculations with variable finite element sizes not smaller than the localized crack band width (corresponding to the material characteristic length). The method uses special localization elements in which a zone of characteristic size, undergoing strain softening, is coupled with layers (called ‘springs’) which undergo elastic unloading and are normal to the principal stress directions. Because of the coupling of strain‐softening zone with elastic layers, the computations of the microplane model need to be iterated in each finite element in each load step, which increases the computer time. Insensitivity of the proposed method to mesh size is demonstrated by numerical examples. Simulation of various experimental results is shown to give good agreement. Copyright © 2004 John Wiley & Sons, Ltd.     

A fixed‐mesh finite element thermally coupled flow formulation for the numerical analysis of melting processes

In the framework of the finite element method, a temperature‐based thermally coupled flow formulation including phase‐change effects is proposed to study melting processes. The governing equations of the problem, written in terms of its primitive variables, are solved using a generalized streamline operator technique that enables the use of equal interpolation functions for the unknowns: velocity, pressure and temperature. Moreover, a unique fixed finite element mesh is used to avoid the difficulties related to moving meshes. This methodology is applied and assessed in the numerical analysis of a benchmark problem known as the melting process of gallium in a differentially heated recipient using distinct geometric aspect ratios. Copyright © 2001 John Wiley & Sons, Ltd.     

A super‐element for crack analysis in the time domain

A super‐element for the dynamic analysis of two‐dimensional crack problems is developed based on the scaled boundary finite‐element method. The boundary of the super‐element containing a crack tip is discretized with line elements. The governing partial differential equations formulated in the scaled boundary co‐ordinates are transformed to ordinary differential equations in the frequency domain by applying the Galerkin's weighted residual technique. The displacements in the radial direction from the crack tip to a point on the boundary are solved analytically without any a priori assumption. The scaled boundary finite‐element formulation leads to symmetric static stiffness and mass matrices. The super‐element can be coupled seamlessly with standard finite elements. The transient response is evaluated directly in the time domain using a standard time‐integration scheme. The stress field, including the singularity around the crack tip, is expressed semi‐analytically. The stress intensity factors are evaluated without directly addressing singular functions, as the limit in their definitions is performed analytically. Copyright © 2004 John Wiley & Sons, Ltd.     

Two‐scale diffusion–deformation coupling model for material deterioration involving micro‐crack propagation

In this paper, we introduce a two‐scale diffusion–deformation coupled model that represents the aging material deterioration of two‐phase materials involving micro‐crack propagations. The mathematical homogenization method is applied to relate the micro‐ and macroscopic field variables, and a weak coupling solution method is employed to solve the two‐way coupling phenomena between the diffusion of scalar fields and the deformation of quasi‐brittle solids. The macroscopic mechanical behavior represented by the derived two‐scale two‐way coupled model reveals material nonlinearity due to micro‐scale cracking induced by the scalar‐field‐induced deformation, which can be simulated by the finite cover method. After verifying the fundamental validity of the proposed model and the analysis method, we perform a simple numerical example to demonstrate their ability to predict aging material deterioration. Copyright © 2010 John Wiley & Sons, Ltd.     

Coupled meshfree and fractal finite element method for mixed mode two‐dimensional crack problems

This paper presents a coupling technique for integrating the element‐free Galerkin method (EFGM) with the fractal finite element method (FFEM) for analyzing homogeneous, isotropic, and two‐dimensional linear‐elastic cracked structures subjected to mixed‐mode (modes I and II) loading conditions. FFEM is adopted for discretization of the domain close to the crack tip and EFGM is adopted in the rest of the domain. In the transition region interface elements are employed. The shape functions within interface elements which comprise both the EFG and the finite element (FE) shape functions, satisfies the consistency condition thus ensuring convergence of the proposed coupled EFGM–FFEM. The proposed method combines the best features of EFGM and FFEM, in the sense that no special enriched basis functions or no structured mesh with special FEs are necessary and no post‐processing (employing any path independent integrals) is needed to determine fracture parameters, such as stress‐intensity factors (SIFs) and T‐stress. The numerical results show that SIFs and T‐stress obtained using the proposed method are in excellent agreement with the reference solutions for the structural and crack geometries considered in the present study. Also, a parametric study is carried out to examine the effects of the integration order, the similarity ratio, the number of transformation terms, and the crack length to width ratio on the quality of the numerical solutions. A numerical example on mixed‐mode condition is presented to simulate crack propagation. As in the proposed coupled EFGM–FFEM at each increment during the crack propagation, the FFEM mesh (around the crack tip) is shifted as it is to the new updated position of the crack tip (such that FFEM mesh center coincides with the crack tip) and few meshless nodes are sprinkled in the location where the FFEM mesh was lying previously, crack‐propagation analysis can be dramatically simplified. Copyright © 2010 John Wiley & Sons, Ltd.     

Continuum Shape Sensitivity analysis of a mode-I fracture in functionally graded materials

This paper presents a new method for conducting a continuum shape sensitivity analysis of a crack in an isotropic, linear-elastic, functionally graded material. This method involves the material derivative concept from continuum mechanics, domain integral representation of the J-integral and direct differentiation. Unlike virtual crack extension techniques, no mesh perturbation is needed to calculate the sensitivity of stress-intensity factors. Since the governing variational equation is differentiated prior to the process of discretization, the resulting sensitivity equations are independent of approximate numerical techniques, such as the meshless method, finite element method, boundary element method, or others. In addition, since the J-integral is represented by domain integration, only the first-order sensitivity of the displacement field is needed. Several numerical examples are presented to calculate the first-order derivative of the J-integral, using the proposed method. Numerical results obtained using the proposed method are compared with the reference solutions obtained from finite-difference methods for the structural and crack geometries considered in this study.     

A new method for modelling cohesive cracks using finite elements

A model which allows the introduction of displacements jumps to conventional finite elements is developed. The path of the discontinuity is completely independent of the mesh structure. Unlike so‐called ‘embedded discontinuity’ models, which are based on incompatible strain modes, there is no restriction on the type of underlying solid finite element that can be used and displacement jumps are continuous across element boundaries. Using finite element shape functions as partitions of unity, the displacement jump across a crack is represented by extra degrees of freedom at existing nodes. To model fracture in quasi‐brittle heterogeneous materials, a cohesive crack model is used. Numerical simulations illustrate the ability of the method to objectively simulate fracture with unstructured meshes. Copyright © 2001 John Wiley & Sons, Ltd.     

A perturbation method for stochastic meshless analysis in elastostatics

A stochastic meshless method is presented for solving boundary‐value problems in linear elasticity that involves random material properties. The material property was modelled as a homogeneous random field. A meshless formulation was developed to predict stochastic structural response. Unlike the finite element method, the meshless method requires no structured mesh, since only a scattered set of nodal points is required in the domain of interest. There is no need for fixed connectivities between nodes. In conjunction with the meshless equations, classical perturbation expansions were derived to predict second‐moment characteristics of response. Numerical examples based on one‐ and two‐dimensional problems are presented to examine the accuracy and convergence of the stochastic meshless method. A good agreement is obtained between the results of the proposed method and Monte Carlo simulation. Since mesh generation of complex structures can be a far more time‐consuming and costly effort than the solution of a discrete set of equations, the meshless method provides an attractive alternative to finite element method for solving stochastic mechanics problems. Copyright © 2001 John Wiley & Sons, Ltd.     

Finite element simulation of dynamic crack propagation process using an arbitrary Lagrangian Eulerian formulation

In this paper, the arbitrary Lagrangian Eulerian formulation is employed for finite element modelling of dynamic crack propagation problem. The application phase simulation of computational dynamic fracture is applied to model by which the crack propagation history and variation of crack velocity are predicted using the material dynamic fracture toughness. The dynamic solution of problem is accomplished using the implicit time integration method. The convective terms due to mesh‐material motion are taken into account via the convection equation. A robust and efficient mesh motion technique, that its equations need not to be solved at every time step, is employed in Eulerian phase. The mesh connectivity is preserved during the analysis. So, the successive remeshing of model is eliminated. When the dynamic fracture criterion is satisfied for crack growth, the presented algorithm allows the crack to advance by splitting the material particle at the crack tip. The dynamic energy release rate is calculated at each time step to determine dynamic stress intensity factor. The predicted results are compared with those obtained through the experimental study and remeshing technique cited in the literature. The proposed computational algorithm leads to an accurate and efficient simulation of dynamic crack propagation process.     

Some recent developments in numerical modelling of fracture toughness in brittle matrix composites

Finite element analysis was used to study the fracture toughening of a ceramic by a stress induced dilatant transformation of second phase particles. The finite element method was based on a continuum theory which modelled the composite as subcritical material. Transient crack growth was simulated in the finite element mesh by a nodal release technique. The crack's remote tensile opening load was adjusted to maintain the near-tip energy release rate at the level necessary for crack advance. The transformation zone surrounding the crack developed as the crack propagated through the composite. Resistance curves were computed from the analysis; and the results show that during crack advance maximum toughness is achieved before a steady state is reached. The toughening effect of a crack-bridging ductile phase in a brittle material may be predicted if ligament deformation is characterized. A plastically deforming ligament constrained by surrounding elastic matrix material is modelled using finite elements and the relevant toughness enhancement information extracted. Comparison is made to model experiments as well as to toughness measured for technologically important materials. The results suggest that debonding along the interface between the ligament and the matrix may enhance the toughening effect of a ductile phase.     

Dynamic crack propagation analysis of orthotropic media by the extended finite element method

Dynamic crack propagation of composites is investigated in this paper based on the recent advances and development of orthotropic enrichment functions within the framework of partition of unity and the extended finite element method (XFEM). The method allows for analysis of the whole crack propagation pattern on an unaltered finite element mesh, defined independent of the existence of any predefined crack or its propagation path. A relatively simple, though efficient formulation is implemented, which consists of using a dynamic crack initiation toughness, a crack orientation along the maximum circumferential stress, and a simple equation to presume the crack speed. Dynamic stress intensity factors (DSIFs) are evaluated by means of the domain separation integral method. The governing elastodynamics equation is first transformed into a standard weak formulation and is then discretized into an XFEM system of time dependent equations, to be solved by the unconditionally stable Newmark time integration scheme. A number of benchmark and test problems are simulated and the results are compared with available reference results.     

Phase‐field modeling of ductile fracture at finite strains: A robust variational‐based numerical implementation of a gradient‐extended theory by micromorphic regularization

This work provides a robust variational‐based numerical implementation of a phase field model of ductile fracture in elastic–plastic solids undergoing large strains. This covers a computationally efficient micromorphic regularization of the coupled gradient plasticity‐damage formulation. The phase field approach regularizes sharp crack surfaces within a pure continuum setting by a specific gradient damage modeling with geometric features rooted in fracture mechanics. It has proven immensely successful with regard to the analysis of complex crack topologies without the need for fracture‐specific computational structures such as finite element design of crack discontinuities or intricate crack‐tracking algorithms. The proposed gradient‐extended plasticity‐damage formulation includes two independent length scales that regularize both the plastic response as well as the crack discontinuities. This ensures that the damage zones of ductile fracture are inside of plastic zones or vice versa and guarantees on the computational side a mesh objectivity in post‐critical ranges. The proposed setting is rooted in a canonical variational principle. The coupling of gradient plasticity to gradient damage is realized by a constitutive work density function that includes the stored elastic energy and the dissipated work due to plasticity and fracture. The latter represents a coupled resistance to plasticity and damage, depending on the gradient‐extended internal variables that enter plastic yield functions and fracture threshold functions. With this viewpoint on the generalized internal variables at hand, the thermodynamic formulation is outlined for gradient‐extended dissipative solids with generalized internal variables that are passive in nature. It is specified for a conceptual model of von Mises‐type elasto‐plasticity at finite strains coupled with fracture. The canonical theory proposed is shown to be governed by a rate‐type minimization principle, which fully determines the coupled multi‐field evolution problem. This is exploited on the numerical side by a fully symmetric monolithic finite element implementation. An important aspect of this work is the regularization towards a micromorphic gradient plasticity‐damage setting by taking into account additional internal variable fields linked to the original ones by penalty terms. This enhances the robustness of the finite element implementation, in particular, on the side of gradient plasticity. The performance of the formulation is demonstrated by means of some representative examples. Copyright © 2016 John Wiley & Sons, Ltd.     

A self-adaptive finite element approach for simulation of mixed-mode delamination using cohesive zone models

Oscillations observed in the load–displacement response of brittle interfaces modeled by cohesive zone elements in a quasi-static finite element framework are artifacts of the discretization. The typical limit points in this oscillatory path can be traced by application of path-following techniques, or avoided altogether by adequately refining the mesh until the standard iterative Newton–Raphson method becomes applicable. Both strategies however lead to an unacceptably high computational cost and a low efficiency, justifying the development of a process driven hierarchical extension of the discretization used in the process zone of a cohesive crack. A self-adaptive enrichment scheme within individual cohesive zone elements driven by the physics governing the problem, is an efficient solution that does not require further mesh refinements. A two-dimensional mixed-mode example in a general framework with an irreversible cohesive zone law shows that an enriched formulation restores the smoothness of the solution in structures that are discretized in a relatively coarse manner.     

Shape sensitivity analysis in mixed-mode fracture mechanics

 This paper presents a new method for continuum-based shape sensitivity analysis for a crack in a homogeneous, isotropic, and linear-elastic body subject to mixed-mode (modes I and II) loading conditions. The method is based on the material derivative concept of continuum mechanics, domain integral representation of an interaction integral, and direct differentiation. Unlike virtual crack extension techniques, no mesh perturbation is needed in the proposed method to calculate the sensitivity of stress-intensity factors. Since the governing variational equation is differentiated prior to the process of discretization, the resulting sensitivity equations are independent of approximate numerical techniques, such as the finite element method, boundary element method, meshless methods, or others. In addition, since the interaction integral is represented by domain integration, only the first-order sensitivity of the displacement field is needed. Two numerical examples are presented to illustrate the proposed method. The results show that the maximum difference in the sensitivity of stress-intensity factors calculated using the proposed method and reference solutions obtained by analytical or finite-difference methods is less than four percent. Received 19 September 2000     

3D modelling of strong discontinuities in elastoplastic solids: fixed and rotating localization formulations

This paper is concerned with the incorporation of displacement discontinuities into a continuum theory of elastoplasticity for the modelling of localization processes such as cracking in brittle materials. Based on the strong discontinuity approach (SDA) (Computational Mechanics 1993; 12: 277–296) mesh objective 2D and 3D finite element formulations are developed using linear and quadratic 2D elements as well as 8‐noded 3D elements. In the formulation of the finite‐element model proposed in the paper, the analogy with standard formulations is emphasized. The parameter defining the amplitude of the displacement jump within the finite element is condensed out at the material level without employing the standard static condensation technique. This approach results in linearized constitutive equations formally identical to continuum models. Therefore, the standard return mapping algorithm is used to solve the non‐linear equations. In analogy to concepts used in continuum smeared crack models, a rotating formulation of the SDA is proposed in addition to the standard concept of fixed discontinuities. It is shown that the rotating localization approach reduces locking effects observed in analyses based on fixed localization directions. The applicability of the proposed SDA finite‐element model as well as its numerical performance is investigated by means of a three‐dimensional ultimate load analysis of a steel anchor embedded in a concrete block subjected to a shear force. Copyright © 2003 John Wiley & Sons, Ltd.     

On advanced solution strategies to overcome locking effects in strong discontinuity approaches

This paper is concerned with the analysis of locking effects resulting from different orientations of micro‐defects and those of the corresponding macro‐defects. Based on a mixed‐mode material model embedded within the framework of the strong discontinuity approach (SDA), the described locking effect is illustrated by means of a crack analysis of a notched concrete beam. To overcome the deficiency of the proposed finite element model, the original SDA is modified and extended. For that purpose, two different advanced numerical formulations are developed: a rotating crack approach and a multiple crack approach. Restricting the governing equations to the material point level, a standard return‐mapping procedure is applied to the algorithmic formulations of both models. The applicability and the performance of the proposed numerical implementations are investigated by means of a re‐analysis of the two‐dimensional notched concrete beam. Copyright © 2005 John Wiley & Sons, Ltd.     

A mesh objective method for modeling crack propagation using the smeared crack approach

This paper presents a mesh objective method for modeling crack propagation in brittle materials using a conventional finite element formulation. The primary shortcoming of the smeared crack approach is its pathological sensitivity to the mesh orientation, which is manifested by shear locking and stress field misalignment around the crack tip. Such undesirable characteristics preclude the ability to model arbitrary crack propagation at an angle through the mesh. Several techniques are developed to address these shortcomings. First, to preclude shear locking, a modified failure constitutive model is developed, which projects out the spurious stress increments as the crack opens. If a crack exists in an element, a crack tracking algorithm is used to identify the neighboring elements most likely to show crack continuation. This algorithm also identifies a crossover element when a crack passes through adjacent sides of an element. Then, the characteristic element length used in the constitutive equation is changed with the objective of providing the correct failure energy per unit crack length, a procedure called crossover scaling. The examples provided demonstrate that the developed methods work collectively to provide a simple and efficient method for modeling failure in brittle materials without mesh bias.