Enhanced solution control for physically and geometrically non‐linear problems. Part I—the subplane control approach

Geometrically or physically non‐linear problems are often characterized by the presence of critical points with snapping behaviour in the structural response. These structural or material instabilities usually lead to inefficiency of standard numerical solution techniques. Special numerical procedures are therefore required to pass critical points. This paper presents a solution technique which is based on a constraint equation that is defined on a subplane of the degrees‐of‐freedom (dof's) hyperspace or a hyperspace constructed from specific functions of the degrees‐of‐freedom. This unified approach includes many existing methods which have been proposed by various authors. The entire computational process is driven from only one control function which is either a function of a number of degrees‐of‐freedom (local subplane method) or a single automatically weighted function that incorporates all dof's directly or indirectly (weighted subplane method). The control function is generally computed in many points of the structure, which can be related to the finite element discretization. Each point corresponds to one subplane. In the local subplane method, the subplane with the control function that drives the load adaptation is selected automatically during the deformation process. Part I of this two‐part series of papers fully elaborates the proposed solution strategy, including a fully automatic load control, i.e. load estimation, adaptation and correction. Part II presents a comparative analysis in which several choices for the control function in the subplane method are confronted with classical update algorithms. The comparison is carried out by means of a number of geometrically and physically non‐linear examples. General conclusions are drawn with respect to the efficiency and applicability of the subplane solution control method for the numerical analysis of engineering problems. Copyright © 1999 John Wiley & Sons Ltd.     

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

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 Griffith's 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 © 2013 John Wiley & Sons, Ltd.     

Solution strategies for the delamination analysis based on a combination of local‐control arc‐length and line searches

The paper describes a local‐control arc‐length method which can be combined with various forms of line‐search procedure. In particular, a new ‘double‐line‐search’ method is developed, which significantly improves the solution procedure and turns out to be efficient and robust. Although the potential range of applications is wide, the method is here limited to the finite element analysis of delamination in a laminated composite using a cohesive‐zone model combined with interface elements. Three problems have been analysed and comparisons have been made with experimental results. Copyright © 2003 John Wiley & Sons, Ltd.     

A continuum damage mechanics method for fracture simulation of quasi‐brittle materials

This paper addresses a novel continuum damage‐based method for simulating failure process of quasi‐brittle materials starting from local damage initiation to final fracture. In the developed method, the preset characteristic length field is used to evaluate damage instead of element, which is used to reduce the spurious sensitivity. In addition, damage is only updated in the most dangerous location at a time for considering stress redistribution due to damage evolution, which is used to simulate competitive fracture process. As cases study, representative numerical simulations of two benchmark tests are given to verify the performance of the developed continuum damage‐based method together with a used damage model. The simulation results of the crack paths for two concrete specimens obtained from the developed method matched well with the corresponding experimental results. The results show that the developed continuum damage‐based method is effective and can be used to simulate damage and fracture process of brittle or quasi‐brittle materials. And the simulation results based on the developed method depend only the preset characteristic length field and not grid mesh.     

A new path‐following constraint for strain‐softening finite element simulations

The application of strain‐softening constitutive relations to model the failure modes of real‐life structures is faced to numerical difficulties related to instabilities that appear as sharp snap‐backs of the structural response. A path‐following method has to complement the solution algorithm to achieve convergence despite these critical points. Because of the sharpness of the snap‐backs, it is believed essential that the path‐following constraint distinguish between a purely elastic unloading and a dissipative path. For that purpose, a new constraint based on the maximal value of the elastic predictor for the yield function is proposed. As it is highly non linear, a specific solution algorithm is required. The robustness of this constraint is illustrated by three applications: the study of crack propagations by means of a cohesive zone model, the failure of a structure submitted to nonlocal damage and the simulation of a nonlocal strain‐softening plastic specimen. Copyright © 2004 John Wiley & Sons, Ltd.     

A novel singular node‐based smoothed finite element method (NS‐FEM) for upper bound solutions of fracture problems

It is well known that the lower bound to exact solutions in linear fracture problems can be easily obtained by the displacement compatible finite element method (FEM) together with the singular crack tip elements. It is, however, much more difficult to obtain the upper bound solutions for these problems. This paper aims to formulate a novel singular node‐based smoothed finite element method (NS‐FEM) to obtain the upper bound solutions for fracture problems. In the present singular NS‐FEM, the calculation of the system stiffness matrix is performed using the strain smoothing technique over the smoothing domains (SDs) associated with nodes, which leads to the line integrations using only the shape function values along the boundaries of the SDs. A five‐node singular crack tip element is used within the framework of NS‐FEM to construct singular shape functions via direct point interpolation with proper order of fractional basis. The mix‐mode stress intensity factors are evaluated using the domain forms of the interaction integrals. The upper bound solutions of the present singular NS‐FEM are demonstrated via benchmark examples for a wide range of material combinations and boundary conditions. Copyright © 2010 John Wiley & Sons, Ltd.     

A Gibbs‐energy‐barrier‐based computational micro‐sphere model for the simulation of martensitic phase‐transformations

We introduce a material model for the simulation of polycrystalline materials undergoing solid‐to‐solid phase‐transformations. As a basis, we present a scalar‐valued phase‐transformation model where a Helmholtz free energy function depending on volumetric and deviatoric strain measures is assigned to each phase. The analysis of the related overall Gibbs energy density allows for the calculation of energy barriers. With these quantities at hand, we use a statistical‐physics‐based approach to determine the resulting evolution of volume fractions. Though the model facilitates to take into account an arbitrary number of solid phases of the underlying material, we restrict this work to the simulation of phase‐transformations between an austenitic parent phase and a martensitic tension and compression phase. The scalar model is embedded into a computational micro‐sphere formulation in view of the simulation of three‐dimensional boundary value problems. The final modelling approach necessary for macroscopic simulations is accomplished by a finite element formulation, where the local material behaviour at each integration point is governed by the response of the micro‐sphere model.Copyright © 2014 John Wiley & Sons, Ltd.     

A High‐order extended finite element method for extraction of mixed‐mode strain energy release rates in arbitrary crack settings based on Irwin's integral

An analytical formulation based on Irwin's integral and combined with the extended finite element method is proposed to extract mixed‐mode components of strain energy release rates in linear elastic fracture mechanics. The proposed formulation extends our previous work to cracks in arbitrary orientations and is therefore suited for crack propagation problems. In essence, the approach employs high‐order enrichment functions and evaluates Irwin's integral in closed form, once the linear system is solved and the algebraic degrees of freedom are determined. Several benchmark examples are investigated including off‐center cracks, inclined cracks, and two crack growth problems. On all these problems, the method is shown to work well, giving accurate results. Moreover, because of its analytical nature, no special post‐processing is required. Thus, we conclude that this method may provide a good and simple alternative to the popular J‐integral method. In addition, it may circumvent some of the limitations of the J‐integral in 3D modeling and in problems involving branching and coalescence of cracks. Copyright © 2013 John Wiley & Sons, Ltd.     

Study of crack arrest mechanism and dynamic behaviour using arc‐bottom specimen under impacts

In this paper, a large size trapezoidal open crack with arc‐bottom (TOCAB) configuration specimen was proposed. Experiments were carried out using 0°, 60°, 90° and 120° TOCAB specimens under drop weight impacts with loading rate from 150 to 350 GPa/s. The crack length and crack speed were calculated by using fractal method. It shows that the crack speed corrected by fractal method is more than 12% the original speed calculated using the straight length. The AUTODYN code was employed to simulate crack propagation path and crack speed. The dynamic stress intensity factors (DSIFs) were calculated using ABAQUS code. The experimental and numerical results show that all the three TOCAB specimens have arresting function on running cracks; crack propagation velocity increases with loading rates; the critical DSIFs at initiation are higher than those at propagation. The study presents how to arrest a running crack and how to protect significant structures from being completely damaged.     

A Koiter‐Newton approach for nonlinear structural analysis

A new approach termed the Koiter‐Newton is presented for the numerical solution of a class of elastic nonlinear structural response problems. It is a combination of a reduction method inspired by Koiter's post‐buckling analysis and Newton arc‐length method so that it is accurate over the entire equilibrium path and also computationally efficient in the presence of buckling. Finite element implementation based on element independent co‐rotational formulation is used. Various numerical examples of buckling sensitive structures are presented to evaluate the performance of the method. The examples demonstrate that the method is robust and completely automatic and that it outperforms traditional path‐following techniques. This improved efficiency will open the door for the direct use of detailed nonlinear finite element models in the design optimization of next generation flight and launch vehicles. Copyright © 2013 John Wiley & Sons, Ltd.     

A rapid path‐length searching procedure for multi‐axial fatigue cycle counting

In this paper, a series of advanced searching algorithms have been examined and implemented for accelerating multi‐axial fatigue cycle counting efforts when dealing with large time histories. In a computerized calculation of the path‐length dependent cycle counting method, most of the central processor unit's (CPU) time is spent on searching for the maximum range or distance in a stress or strain space. A brute‐force search is the simplest to implement, and will always find a solution if it exists. However, its cost, in many practical problems, tends to grow exponentially as the size of the loading spectrum increases with a search time measured in the order of O(n²), where n is the number of spectrum data points. In contrast, a form of Andrew's monotone chain algorithm, as demonstrated in this paper, can remarkably reduce the solution time to the order of O(n log n). The effectiveness of the new path‐length searching procedure is demonstrated by a series of worked examples with a varying degree of non‐proportionality in multi‐axial loading history.     

Digital energy grade‐based approach for crack path prediction based on 2D X‐ray CT images of geomaterials

Cracking is an important phenomenon in the failure of geomaterials. The prediction of crack paths is difficult and challenging because of the randomness and uncertainty in the cracking behaviors of geomaterials. In this paper, to predict crack paths based on 2D X‐ray computed tomography (CT) images, a digital energy grade‐based approach is proposed, and the corresponding energy principle is established. Excellent consistency of crack paths is found between the predicted crack path and the real crack path of the specimen. The numerical results indicate that the proposed approach provides a useful way to predict cracking paths in geomaterials. Meanwhile, the stress and displacement fields before and after the specimen fails can be obtained by combining the proposed method with finite element (FE) analysis. In total, this proposed method can be applied not only to monitoring the health of but also to the stability analysis of engineering structures during engineering activities.     

A numerical method to simulate ductile failure of tensile plates with interacting through‐wall cracks

This paper proposes a simple numerical method to simulate ductile failure behaviours of tensile plates with interacting through‐wall cracks. The method is based on the finite element damage analysis using the stress‐modified fracture strain damage model. To validate the proposed method, simulated results are compared with a total of 23 published experimental data of flat tensile plates with interacting through‐wall cracks. Despite its simplicity, the proposed method well predicted the experimental maximum loads of tensile plates with interacting cracks, including the loads for crack coalescence. Systematic analyses are also performed to investigate the effect of the element size used in the finite element damage analysis.     

A heuristic method for obtaining quasi ARL‐unbiased p‐Charts

It is known that control charts based on equal tail probability limits are ARL biased when the distribution of the plotted statistic is skewed. This is the case for p‐Charts that serve to monitor processes on the basis of the binomial distribution. For the particular case of the standard three‐sigma Shewhart p‐Chart, which is built on the basis of the binomial to normal distribution approximation, this ARL‐biased condition is particularly severe, and it greatly affects its monitoring capability. Surprisingly, in spite of this, the standard p‐Chart is still widely used and taught. Through a literature search, it was identified that several, simple to use, improved alternative p‐Charts had been proposed over the years; however, at first instance, it was not possible to determine which of them was the best. In order to identify the alternative that excelled, an ARL performance comparison was carried out in terms of their ARL bias severity level (ARL_BSL) and their In‐Control ARL (ARL₀). The results showed that even the best performing alternative charts would often be ARL‐biased or have nonoptimal ARL₀. To improve on the existing alternatives, the “Kmod p‐Chart” was developed; it offers easiness of use, superior ARL performance, and a simple and effective method for verifying its ARL‐bias condition.     

Path‐following analysis of thin‐walled structures and comparison with asymptotic post‐critical solutions

A path‐following non‐linear elastic analysis for structures composed of assemblages of flat slender elastic panels is presented. The proposed path‐following method employs FEM technology and a kinematical model to analyse these structures using a Koiter asymptotic approach. As a result it is possible to verify the accuracy achieved by the asymptotic method. The proposed mixed path‐following formulation is both efficient and robust with regards to the locking extrapolation phenomenon that strongly affects compatible formulations. The use of an HC finite element makes it possible to avoid the problem of the finite rotations in the space, maintaining a high degree of continuity and making the numeric formulation simple and efficient. Copyright © 2002 John Wiley & Sons, Ltd.     

An efficient solution method for finite element ring‐rolling simulation

Finite element ring‐rolling simulation gives rise to poor conditioned non‐linear equations that require repeated solution. The associated computational costs are extreme making analysis impracticable in industry. This paper is concerned with a solution strategy that addresses this problem and involves the combined use of an arbitrary Lagrangian–Eulerian (ALE) formulation and a successive preconditioned conjugate gradient method (SPCGM). This approach, coupled to a finite element flow formulation, is shown to offer considerable computational savings. Through the combined use of the ALE flow formulation and the SPCGM the stability and condition of the non‐linear systems is enhanced. This purely iterative approach takes advantage of the slowly evolving velocity field and the self‐preconditioning offered by the SPCGM. The performance of the solver is compared against well‐known alternatives for varying problem sizes. The approach is shown to be generic but in particular makes ring‐rolling simulation a more practicable proposition. Copyright © 2000 John Wiley & Sons, Ltd.     

A frequency‐domain formulation of MCE method for multi‐axial random loadings

Many multi‐axial fatigue limit criteria are formalized as a linear combination of a shear stress amplitude and a normal stress. To identify the shear stress amplitude, appropriate conventional definitions, as the minimum circumscribed circle (MCC) or ellipse (MCE) proposals, are in use. Despite computational improvements, deterministic algorithms implementing the MCC/MCE methods are exceptionally time‐demanding when applied to “coiled” random loading paths resulting from in‐service multi‐axial loadings and they may also provide insufficiently robust and reliable results. It would be then preferable to characterize multi‐axial random loadings by statistical re‐formulations of the deterministic MCC/MCE methods. Following an early work of Pitoiset et al., this paper presents a statistical re‐formulation for the MCE method. Numerical simulations are used to compare both statistical re‐formulations with their deterministic counterparts. The observed general good trend, with some better performance of the statistical approach, confirms the validity, reliability and robustness of the proposed formulation.     

A magnitude‐robust control chart for monitoring and estimating step changes for normal process means

Statistical process control charts are intended to assist operators of a usually stable system in monitoring whether a change has occurred in the process. When a change does occur, the control chart should detect it quickly. If the operator can also be provided information that aids in the search for the special cause, then critical off‐line time can be saved. We investigate a process‐monitoring tool that not only provides speedy detection regardless of the magnitude of the process shift, but also supplies useful change point statistics. A likelihood ratio approach can be used to develop a control chart for permanent step change shifts of a normal process mean. The average run length performance for this chart is compared to that of several cumulative sum (CUSUM) charts. Our performance comparisons show that this chart performs better than any one CUSUM chart over an entire range of potential shift magnitudes. The likelihood ratio approach also provides point and interval estimates for the time and magnitude of the process shift. These crucial change‐point diagnostics can greatly enhance special cause investigation. Copyright © 2002 John Wiley & Sons, Ltd.     

An iterative method for coupling of deformation and failure mechanisms on different length scales

An iterative method for coupling of numerical simulations on two length scales is presented. The computations on the microscale and on the macroscale are linked via a suitable macroscopic constitutive law. The parameters of this material law depend on the deformation history and are obtained from simulations using microstructurally representative volume elements (RVEs) subjected to strain paths derived from the associated material points in the macroscopic structure. Thus, different constitutive parameter sets are assigned to different regions of the macrostructure. The microscopic and macroscopic simulations are performed iteratively and interact mutually via the strain paths and the constitutive parameters, respectively. As an example, the strip tension test for a porous material is modelled using the finite element (FE) method. The coupling procedure, the material law and its numerical implementation are described. The method is shown to allow for a detailed simulation of the deformation mechanisms both on the micro‐ and the macroscale as well as for an assessment of their interactions while keeping the computational efforts reasonably low. Copyright © 2000 John Wiley & Sons, Ltd.