Formulation and Application of the Initially Rigid Cohesive Finite Element Method

The presented procedure for cohesive crack propagation is based on an adaptive finite element (FE) implementation, which enables the introduction of cohesive surfaces in dependence on the current crack state. In contrast to already existing formulations, the focus of the present model lies on failure processes that can be described at quasi-static conditions within an implicit framework. Furthermore, an extension for mesh independent crack propagation in terms of an additional mesh adaptive formulation is presented. By the evaluation of the failure criterion considering the preferred crack direction, a new crack tip coordinate is computed and the discretization is accordingly adjusted. The remaining mesh is modified for the new boundary representation. The application of the proposed method is shown by the numerical investigation of a concrete fracture specimen from an experimental research project. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Adaptive extended isogeometric analysis based on PHT-splines for thin cracked plates and shells with Kirchhoff–Love theory

In this paper, a posteriori error estimation and mesh adaptation approach for thin plate and shell structures of through-the-thickness crack is presented. This method uses the extended isogeometric analysis (XIGA) based on PHT-splines (Polynomial splines over Hierarchical T-meshes), which is abbreviated as XIGA-PHT. In XIGA-PHT, the isogeometric displacement approximation is locally enriched with enrichment functions, which efficiently capture the displacement discontinuity across the crack face as well as the stress singularity in the vicinity of the crack tip. On the one hand, the rotational degrees of freedom (RDOFs) are not required in Kirchhoff–Love theory, which drastically reduces the complexity of enrichment mode and computational scale for crack analysis. On the other hand, the PHT-splines basis functions can automatically satisfy the requirement of C¹-continuity for the Kirchhoff–Love theory. Moreover, the PHT-splines facilitate the local refinement, which is the deficiency of NURBS-based isogeometric formulations. The local refinement is highly suitable for adaptive analysis. The stress recovery-based posteriori error estimator combined with the superconvergent patch recovery (SPR) technique is used to evaluate the approximate local discretization error. A new strategy for selecting enriched recovered functions in the enriched areas was proposed. Special functions extracted from the asymptotic stress solutions are applied to obtain the recovered stress field in the enriched area. The results of stress intensity factors or J-integral values obtained by the adaptive XIGA-PHT are compared with reference solutions. Several thin plate and shell illustrative examples demonstrate the effectiveness and accuracy of the proposed adaptive XIGA-PHT.     

Configurational-Force-Based Adaptive FE Solver for a Phase Field Model of Fracture

The computational modeling of failure mechanisms in solids due to fracture based on sharp crack discontinuities suffers in situations with complex crack topologies. This can be overcome by diffusive crack modeling, based on the introduction of a crack phase field as outlined in [1, 2]. Following these formulations, we outline a thermodynamically consistent framework for phase field models of crack propagation in elastic solids, develop incremental variational principles and, as an extension to [1, 2], consider their numerical implementations by an efficient h-adaptive finite element method. A key problem of the phase field formulation is the mesh density, which is required for the resolution of the diffusive crack patterns. To this end, we embed the computational framework into an adaptive mesh refinement strategy that resolves the fracture process zones. We construct a configurational-force-based framework for h-adaptive finite element discretizations of the gradient-type diffusive fracture model. We develop a staggered computational scheme for the solution of the coupled balances in physical and material space. The balance in the material space is then used to set up indicators for the quality of the finite element mesh and accounts for a subsequent h-type mesh refinement. The capability of the proposed method is demonstrated by means of a numerical example. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

On Configurational-Force-Driven Crack Propagation and Adaptive Mesh Refinement in Finite Inelasticity

We introduce a consistent variational framework for inelasticity at finite strains, yielding dual balances in physical and material space as the Euler equations. The formulation is employed for the simultaneous usage of configurational forces as both driving forces for crack propagation as well as h-adaptive mesh refinement. The theoretical basis builds upon a global balance of internal and external power, where the mechanical response is exclusively governed by two scalar functions, the free energy function and a dissipation potential. The resulting variational structure is exploited in the context of fracture mechanics and yields evolution equations for internal variables. In the discrete setting, we present a geometry model fully separated from the finite element mesh structure that represents structural changes of the material configuration due to crack propagation. Advanced meshing algorithms provide an optimal discretization at the crack tip. Local and global criteria are obtained via error estimators based on configurational forces being interpreted as indicators of an energetic misfit due to an insufficient discretization. The numerical handling is decomposed into a staggered algorithm scheme for the dual set of equilibrium equations in material and physical space and efficient mesh generation tools. Exemplary numerical examples are considered to illustrate the method and to underline the effects of inelastic material behaviour in the presented context. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Adaptive finite element analysis of structures under transient dynamic loading using modal superposition

In this paper, the error estimation and adaptive strategy developed for the linear elastodynamic problem under transient dynamic loading based on the Z–Z criterion is utilized for 2D and plate bending problems. An automatic mesh generator based on “growth meshing” is utilized effectively for adaptive mesh refinement. Optimal meshes are obtained iteratively corresponding to the prescribed domain discretization error limit and for a chosen number of basis modes satisfying modal truncation errors. Numerous examples show the effectiveness of the integrated approach in achieving the target accuracy in finite element transient dynamic analysis.     

A strong discontinuity based adaptive refinement approach for the modeling of crack branching

In the current work, the physical phenomena of dynamic fracture of brittle materials involving crack growth, acceleration and consequent branching is simulated. The numerical modeling is based on the approach where the failure in the form of cracks or shear bands is modeled by a jump in the displacement field, the so called ‘strong discontinuity’. The finite element method is employed with this strong discontinuity approach where each finite element is capable of developing a strong discontinuity locally embedded into it. The focus in this work is on branching phenomena which is modeled by an adaptive refinement method by solving a new sub-boundary value problem represented by a finite element at the growing crack tip. The sub-boundary value problem is subjected to a certain kinematic constraint on the boundary in the form of a linear deformation constraint. An accurate resolution of the state of material at the branching crack tip is achieved which results in realistic dynamic fracture simulations. A comparison of resulting numerical simulations is provided with the experiment of dynamic fracture from the literature. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Adaptive mesh refinement for the control of cost and quality in finite element analysis

Adaptive strategies are a necessary tool to make finite element analysis applicable to engineering practice. In this paper, attention is restricted to mesh adaptivity. Traditionally, the most common mesh adaptive strategies for linear problems are used to reach a prescribed accuracy. This goal is best met with an h-adaptive scheme in combination with an error estimator. In an industrial context, the aim of the mechanical simulations in engineering design is not only to obtain greatest quality but more often a compromise between the desired quality and the computation cost (CPU time, storage, software, competence, human cost, computer used). In this paper, we propose the use of alternative mesh refinement criteria with an h-adaptive procedure for 3D elastic problems. The alternative mesh refinement criteria (MR) are based on: prescribed number of elements with maximum accuracy, prescribed CPU time with maximum accuracy and prescribed memory size with maximum accuracy. These adaptive strategies are based on a technique of error in constitutive relation (the process could be used with other error estimators) and an efficient adaptive technique which automatically takes into account the steep gradient areas. This work proposes a 3D method of adaptivity with the latest version of the INRIA automatic mesh generator GAMHIC3D.     

Simulation of brittle crack growth in functionally graded materials

We consider a thermodynamic consistent framework for crack propagation by applying a dissipation inequality to a time dependent migrating control volume. The direction of crack growth is obtained in terms of material forces as a result of the principle of maximum dissipation. In the numerical implementation a staggered algorithm – deformation update for fixed geometry followed by geometry update for fixed deformation – is employed within each time increment. The corresponding mesh is generated by combining Delaunay triangulation with local mesh refinement. A numerical example with inhomogeneous material properties illustrates the capability of the resulting algorithm. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Simulation of Crack Propagation under Small-Scale Yielding by means of a Non-local GTN-Model

Today, the local approach to fracture is widely applied to simulate the failure of specimens. For ductile damage processes the Gurson-Tvergaard-Needleman model is the quasi-standard. In the last time non-local extensions allowed a mesh-size independent simulation of crack growth. However, most publications dealing with this subject focus upon the convergence regarding global quantities such as the load-displacement relation. Minor attention is paid to the fields directly at the crack tip. Correspondingly, the interrelationship between the intrinsic length of the model and relevant microscopic damage processes at the crack tip is only partly established until now. In the present study the crack propagation is simulated for an implicitly gradient enriched GTN-model within a boundary layer in order to overcome influences of the specimen geometry. The different stages of damage evolution are resolved by a fine mesh. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Four-triangles adaptive algorithms for RTIN terrain meshes

We present a refinement and coarsening algorithm for the adaptive representation of Right-Triangulated Irregular Network (RTIN) meshes. The refinement algorithm is very simple and proceeds uniformly or locally in triangle meshes. The coarsening algorithm decreases mesh complexity by reducing unnecessary data points in the mesh after a given error criterion is applied. We describe the most important features of the algorithms and give a brief numerical study on the propagation associated with the adaptive scheme used for the refinement algorithm. We also present a comparison with a commercial tool for mesh simplification, Rational Reducer, showing that our coarsening algorithm offers better results in terms of accuracy of the generated meshes.     

The effect of microdamage on fatigue-crack growth

The theory of fatigue-crack growth, based on a synthesis of fracture mechanics and continuum mechanics of microdamage accumulation, is applied to the problem of crack growth under cyclic loading, taking into account the plastic deformations in the tip zones. The model of a thin plastic zone, which is a region of considerable plastic deformations at the crack tip, is supplemented by taking into account the effect of microdamage on the value of the specific fracture work and the limit stresses in the tip zone. Governing equations which describe fatigue-crack growth taking these factors into account are derived. The effect of the material characteristics and the load parameters on the growth rate and the distribution of microdamage in the tip zone and on its extensions is investigated by a computational experiment. Particular attention is given to the initial stage when crack growth may occur abruptly and the growth rate depends substantially on the initial conditions     

Configurational-force-based structural updates and adaptive remeshing strategies

The application of configurational forces in h -adaptive strategies for fracture mechanics and inelasticity is investigated. Starting from a global Clausius-Planck inequality, dual equilibrium conditions are derived by means of a Coleman-type exploitation method. The remaining reduced dissipation inequality is used for the derivation of evolution equations for the internal variables. In fracture mechanics, crack loading conditions as well as a normality rule for the crack propagation are obtained. In the discrete setting, the crack propagation is governed by a configurational-force-driven update of the geometry model. The material balance equation is used to set up a h -adaptive refinement indicator. A relative global criterion is defined used for the decision on mesh refinement. In addition, a criterion on the element level is evaluated controlling the local refinement procedure. The capability of the proposed procedures is demonstrated by means of numerical examples. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

On using meshes of mixed element types in adaptive finite element analysis

Four different automatic mesh generators capable of generating either triangular meshes or hybrid meshes of mixed element types have been used in the mesh generation process. The performance of these mesh generators were tested by applying them to the adaptive finite element refinement procedure. It is found that by carefully controlling the quality and grading of the quadrilateral elements, an increase in efficiency over pure triangular meshes can be achieved. Furthermore, if linear elements are employed, an optimal hybrid mesh can be obtained most economically by a combined use of the mesh coring technique suggested by Lo and Lau and a selective removal of diagonals from the triangular element mesh. On the other hand, if quadratic elements are used, it is preferable to generate a pure triangular mesh first, and then obtain a hybrid mesh by merging of triangles.     

A generalized cohesive element technique for arbitrary crack motion

A computational method for arbitrary crack motion through a finite element mesh, termed as the generalized cohesive element technique, is presented. In this method, an element with an internal discontinuity is replaced by two superimposed elements with a combination of original and imaginary nodes. Conventional cohesive zone modeling, limited to crack propagation along the edges of the elements, is extended to incorporate the intra-element mixed-mode crack propagation. Proposed numerical technique has been shown to be quite accurate, robust and mesh insensitive provided the cohesive zone ahead of the crack tip is resolved adequately. A series of numerical examples is presented to demonstrate the validity and applicability of the proposed method.     

Adaptive Runge-Kutta Discontinuous Galerkin Methods with Modified Ghost Fluid Method for Simulating the Compressible Two-Medium Flow

In this paper, we investigate using the adaptive Runge-Kutta discontinuous Galerkin (RKDG) methods with the modified ghost fluid method (MGFM) in conjunction with the adaptive RKDG methods for solving the level set function to simulate the compressible two-medium flow in one and two dimensions. A shock detection technique (KXRCF method) is adopted as an indicator to identify the troubled cell, which serves for further numerical limiting procedure which uses a modified TVB limiter to reconstruct different degrees of freedom and an adaptive mesh refinement procedure. If the computational mesh should be refined or coarsened, and the detail of the implementation algorithm is presented on how to modulate the hanging nodes and redefine the numerical solutions of the two-medium flow and the level set function on such adaptive mesh. Extensive numerical tests are provided to illustrate the proposed adaptive methods may possess the capability of enhancing the resolutions nearby the discontinuities inside of the single medium flow region and material interfacial vicinities of the two-medium flow region.     

High-order adaptive finite element-singly implicit Runge-Kutta methods for parabolic differential equations

We describe an adaptive mesh refinement finite element method-of-lines procedure for solving one-dimensional parabolic partial differential equations. Solutions are calculated using Galerkin's method with a piecewise hierarchical polynomial basis in space and singly implicit Runge-Kutta (SIRK) methods in time. A modified SIRK formulation eliminates a linear systems solution that is required by the traditional SIRK formulation and leads to a new reduced-order interpolation formula. Stability and temporal error estimation techniques allow acceptance of approximate solutions at intermediate stages, yielding increased efficiency when solving partial differential equations. A priori energy estimates of the local discretization error are obtained for a nonlinear scalar problem. A posteriori estimates of local spatial discretization errors, obtained by order variation, are used with the a priori error estimates to control the adaptive mesh refinement strategy. Computational results suggest convergence of the a posteriori error estimate to the exact discretization error and verify the utility of the adaptive technique.This research was partially supported by the U.S. Air Force Office of Scientific Research, Air Force Systems Command, USAF, under Grant Number AFOSR-90-0194; the U.S. Army Research Office under Contract Number DAAL 03-91-G-0215; by the National Science Foundation under Grant Number CDA-8805910; and by a grant from the Committee on Research, Tulane University.     

Towards a more efficient dynamic mesh adaptation methodology for continuum discretization in complex engineering problems

A novel and efficient method of adaptive mesh generation, for dynamically adaptive unstructured grids, is proposed. A locally refined triangulation is constructed on a coarse background mesh, subdividing each triangle in the refinement region R into four congruent sub-triangles iteratively, by connecting edge midpoints, until triangles of a prescribed lengthscale are obtained. The unavoidable propagation outside the refinement region R is restricted to a single triangle in the coarse background mesh. The triangles, in the immediate vicinity of region R, are broken down using the concept of iterated function systems, widely used in fractal modeling, by recursive generation of sub-triangles with a gradation towards the region R triangles. A quantitative assessment of the present algorithm proves its superiority over other comparable models reported in the literature. The time cost of the algorithm is linear, and the method can be easily extended to three dimensions.     

Adaptive refinement procedure for the finite difference method

An adaptive refinement procedure consisting of a localized error estimator and a physically based approach to mesh refinement is developed for the finite difference method. The error estimator is a variation of a successful finite element error estimator. The errors are estimated by computing an error energy norm in terms of discontinuous and continuous stress fields formed from the finite difference results for plane stress problems. The error measure identifies regions of high error which are subsequently refined to improve the result. The local refinement procedure utilizes a recently developed approach for developing finite difference templates to produce a graduated mesh. The adaptive refinement procedure is demonstrated with a problem that contains a well-defined singularity. The results are compared to finite element and uniformly refined finite difference results.     

ADAPTIVITY IN SPACE AND TIME FOR MAGNETOQUASISTATICS

This paper addresses fully space-time adaptive magnetic field computations. We describe an adaptive Whitney finite element method for solving the magnetoquasistatic formulation of Maxwell＇s equations on unstructured 3D tetrahedral grids. Spatial mesh re- finement and coarsening are based on hierarchical error estimators especially designed for combining tetrahedral H（curl）-conforming edge elements in space with linearly implicit Rosenbrock methods in time. An embedding technique is applied to get efficiency in time through variable time steps. Finally, we present numerical results for the magnetic recording write head benchmark problem proposed by the Storage Research Consortium in Japan.