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.     

Non‐linear analysis of thin‐walled members of open cross‐section

The paper presents a means of determining the non‐linear stiffness matrices from expressions for the first and second variation of the Total Potential of a thin‐walled open section finite element that lead to non‐linear stiffness equations. These non‐linear equations can be solved for moderate to large displacements. The variations of the Total Potential have been developed elsewhere by the authors, and their contribution to the various non‐linear matrices is stated herein. It is shown that the method of solution of the non‐linear stiffness matrices is problem dependent. The finite element procedure is used to study non‐linear torsion that illustrates torsional hardening, and the Newton–Raphson method is deployed for this study. However, it is shown that this solution strategy is unsuitable for the second example, namely that of the post‐buckling response of a cantilever, and a direct iteration method is described. The good agreement for both of these problems with the work of independent researchers validates the non‐linear finite element method of analysis. Copyright © 1999 John Wiley & Sons, Ltd.     

Crack–surface displacements for cracks emanating from a circular hole under various loading conditions

The purpose of this paper is to calculate and develop equations for crack–surface displacements for two‐symmetric cracks emanating from a circular hole in an infinite plate for use in strip‐yield crack‐closure models. In particular, the displacements were determined under two loading conditions: (1) remote applied stress and (2) uniform stress applied to a segment of the crack surface (partially loaded crack). The displacements were calculated by an integral‐equation method based on accurate stress–intensity factor equations for concentrated forces applied to the crack surfaces and those for remote applied stress or for a partially loaded crack surface. A boundary‐element code was also used to calculate crack–surface displacements for some selected cases. Comparisons made with crack–surface displacement equations previously developed for the same crack configuration and loading showed significant differences near the location where the crack intersected the hole surface. However, the previous equations were fairly accurate near the crack‐tip location. Herein an improved crack–surface displacement equation was developed for the case of remote applied stress. For the partially loaded crack case, only numerical comparisons were made between the previous equations and numerical integration. A rapid algorithm, based on the integral‐equation method, was developed to calculate these displacements. Because cracks emanating from a hole are quite common in the aerospace industry, accurate displacement solutions are crucial for improving life‐prediction methods based on the strip‐yield crack‐closure models.     

A simplified co‐rotational method for quadrilateral shell elements in geometrically nonlinear analysis

This paper presents a simplified co‐rotational formulation for quadrilateral shell elements inheriting the merit of element‐independence from the traditional co‐rotational approach in literature. With the objective of application to nonlinear analysis of civil engineering structures, the authors further simplify the formulation of the geometrical stiffness using the small strain assumption, which is valid in the co‐rotational approach, with the warping effects considered as eccentricities. Compared with the traditional element‐independent co‐rotational method, the projector is neglected both in the tangent stiffness matrix and in the internal force vector for simplicity in formulation. Meanwhile, a quadrilateral flat shell element allowing for drilling rotations is adopted and incorporated into this simplified co‐rotational algorithm for geometrically nonlinear analysis involved with large displacements and large rotations. Several benchmark problems are presented to confirm the efficiency and accuracy of the proposed method for practical applications.     

A robust preconditioning technique for the extended finite element method

The extended finite element method enhances the approximation properties of the finite element space by using additional enrichment functions. But the resulting stiffness matrices can become ill‐conditioned. In that case iterative solvers need a large number of iterations to obtain an acceptable solution. In this paper a procedure is described to obtain stiffness matrices whose condition number is close to the one of the finite element matrices without any enrichments. A domain decomposition is employed and the algorithm is very well suited for parallel computations. The method was tested in numerical experiments to show its effectiveness. The experiments have been conducted for structures containing cracks and material interfaces. We show that the corresponding enrichments can result in arbitrarily ill‐conditioned matrices. The method proposed here, however, provides well‐conditioned matrices and can be applied to any sort of enrichment. The complexity of this approach and its relation to the domain decomposition is discussed. Computation times have been measured for a structure containing multiple cracks. For this structure the computation times could be decreased by a factor of 2. Copyright © 2010 John Wiley & Sons, Ltd.     

Structural damage identification using transfer matrix with lumped crack properties

A novel damage detection scheme is developed for detecting multiple cracks in beams, based on a transfer matrix (TM) approach. Lumped Crack TM of a beam element with multiple cracks is derived based on lumped crack properties. A cracked beam element is assumed as two intact beam elements connected with a hinge or torsional spring. The crack is modelled as an element of zero length and mass, but with elastic properties. Lumped crack approach is simpler for multiple cracks than the possible alternative methods. The state vector at a node includes displacements, forces and moments at that node; when it is multiplied with TM the state vector at the adjacent node can be obtained. The crack identification strategy used here, involves measuring the initial state vector at a node, in the zone of interest. The displacements at the adjacent nodes are measured and these are predicted using TM. Using an optimization algorithm the unknown crack parameters in the TM are solved by minimizing the deviation between measured and predicted displacements. The method is verified using several numerical models as well as experiments with cracked beams. The TM is shown to be suitable for local identification and also fast and accurate compare to other methods.     

The simple boundary element method for multiple cracks in functionally graded media governed by potential theory: a three‐dimensional Galerkin approach

The simple boundary element method consists of recycling existing codes for homogeneous media to solve problems in non‐homogeneous media while maintaining a purely boundary‐only formulation. Within this scope, this paper presents a ‘simple’ Galerkin boundary element method for multiple cracks in problems governed by potential theory in functionally graded media. Steady‐state heat conduction is investigated for thermal conductivity varying either parabolically, exponentially, or trigonometrically in one or more co‐ordinates. A three‐dimensional implementation which merges the dual boundary integral equation technique with the Galerkin approach is presented. Special emphasis is given to the treatment of crack surfaces and boundary conditions. The test examples simulated with the present method are verified with finite element results using graded finite elements. The numerical examples demonstrate the accuracy and efficiency of the present method especially when multiple interacting cracks are involved. Copyright © 2005 John Wiley & Sons, Ltd.     

Coupled mechanics and finite element for non‐linear laminated piezoelectric shallow shells undergoing large displacements and rotations

A theoretical framework is presented for analysing the coupled non‐linear response of shallow doubly curved adaptive laminated piezoelectric shells undergoing large displacements and rotations. The formulated mechanics incorporate coupling between in‐plane and flexural stiffness terms due to geometric curvature, coupling between mechanical and electric fields, and encompass geometric non‐linearity effects due to large displacements and rotations. The governing equations are formulated explicitly in orthogonal curvilinear co‐ordinates and are combined with the kinematic assumptions of a mixed‐field shear‐layerwise shell laminate theory. Based on the above formulation, a finite element methodology together with an incremental‐iterative technique, based on Newton–Raphson method is formulated. An eight‐node coupled non‐linear shell element is also developed. Various evaluation cases on laminated curved beams and cylindrical panels illustrate the capability of the shell finite element to predict the complex non‐linear behaviour of active shell structures including buckling, which is not captured by linear shell models. The numerical results also show the inherent capability of piezoelectric shell structures to actively induce large displacements through piezoelectric actuators, by jumping between multiple equilibrium states. Copyright © 2005 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.     

A threshold fatigue crack closure model: Part I – model development

A fatigue crack closure model is developed that includes the effects of, and interactions between, the three closure mechanisms most likely to occur at threshold; plasticity, roughness, and oxide. This model, herein referred to as the CROP model (for Closure, Roughness, Oxide, and Plasticity), also includes the effects of out‐of‐plane cracking and multi‐axial loading. These features make the CROP closure model uniquely suited for, but not limited to, threshold applications. Rough cracks are idealized here as two‐dimensional sawtooths, whose geometry induces mixed‐mode crack‐tip stresses. Continuum mechanics and crack‐tip dislocation concepts are combined to relate crack face displacements to crack‐tip loads. Geometric criteria are used to determine closure loads from crack‐face displacements. Finite element results, used to verify model predictions, provide critical information about the locations where crack closure occurs. The CROP model is verified with experimental data in part II of this paper.     

Exact static element stiffness matrices of non‐symmetric thin‐walled curved beams

An evaluation procedure of exact static stiffness matrices for curved beams with non‐symmetric thin‐walled cross section are rigorously presented for the static analysis. Higher‐order differential equations for a uniform curved beam element are first transformed into a set of the first‐order simultaneous ordinary differential equations by introducing 14 displacement parameters where displacement modes corresponding to zero eigenvalues are suitably taken into account. This numerical technique is then accomplished via a generalized linear eigenvalue problem with non‐symmetric matrices. Next, the displacement functions of displacement parameters are exactly calculated by determining general solutions of simultaneous non‐homogeneous differential equations. Finally an exact stiffness matrix is evaluated using force–deformation relationships. In order to demonstrate the validity and effectiveness of this method, displacements and normal stresses of cantilever thin‐walled curved beams subjected to tip loads are evaluated and compared with those by thin‐walled curved beam elements as well as shell elements. Copyright © 2004 John Wiley & Sons, Ltd.     

A bilinear shell element based on a refined shallow shell model

A four‐node shell finite element of arbitrary quadrilateral shape is developed and applied to the solution of static and vibration problems. The element incorporates five generalized degrees of freedom per node, namely the three displacements of the curved middle surface and the two rotations of its normal vector. The stiffness properties of the element are defined using isoparametric principles in a local co‐ordinate system with axes approximately parallel to the edges of the element. The formulation is based on a modern, refined variant of the shallow shell models found from the classical books on shell theory. In addition, the bending behavior of the element is improved with numerical modifications, which include mixed interpolation of the membrane and transverse shear strains. The numerical experiments show that the element is able to compete in accuracy with the highly reputable bilinear elements of the commercial codes ABAQUS and ADINA. The new formulation even outperforms its commercial rivals in problems with strong layers such as vibration problems or problems with concentrated loads. Copyright © 2009 John Wiley & Sons, Ltd.     

A stabilized co‐rotational curved quadrilateral composite shell element

A new curved quadrilateral composite shell element using vectorial rotational variables is presented. An advanced co‐rotational framework defined by the two vectors generated by the four corner nodes is employed to extract pure element deformation from large displacement/rotation problems, and thus an element‐independent formulation is obtained. The present line of formulation differs from other co‐rotational formulations in that (i) all nodal variables are additive in an incremental solution procedure, (ii) the resulting element tangent stiffness is symmetric, and (iii) is updated using the total values of the nodal variables, making solving dynamic problems highly efficient. To overcome locking problems, uniformly reduced integration is used to compute the internal force vector and the element tangent stiffness matrix. A stabilized assumed strain procedure is employed to avoid spurious zero‐energy modes. Several examples involving composite plates and shells with large displacements and large rotations are presented to testify to the reliability, computational efficiency, and accuracy of the present formulation. Copyright © 2011 John Wiley & Sons, Ltd.     

Measurement of discontinuous displacement/strain using mesh-based digital image correlation

In this study, a method for evaluating discontinuous displacement and strain distributions using digital image correlation (DIC) is proposed. A finite element mesh-based DIC method is used for measuring displacements while taking into account displacement and strain discontinuities. Smoothed displacements are thus obtained, and strains are computed from the measured displacements using the finite element mesh again. Discontinuous strains can be obtained by the proposed method using a split finite element mesh. The effectiveness of this method is validated by applying it to measure the displacement and strain in a triaxially woven fabric composite containing numerous free boundaries, to measure displacements around a crack and the displacement and strain around the interface between dissimilar materials. Results show that the discontinuous displacement and strain distributions can be measured by the proposed method. The proposed method is expected to be applicable for the experimental evaluations of various structures and members, including displacement and strain discontinuities such as free boundaries, cracks, and interfaces.     

A robust non‐linear mixed hybrid quadrilateral shell element

In the paper a non‐linear quadrilateral shell element for the analysis of thin structures is presented. The variational formulation is based on a Hu–Washizu functional with independent displacement, stress and strain fields. The interpolation matrices for the mid‐surface displacements and rotations as well as for the stress resultants and strains are specified. Restrictions on the interpolation functions concerning fulfillment of the patch test and stability are derived. The developed mixed hybrid shell element possesses the correct rank and fulfills the in‐plane and bending patch test. Using Newton's method the finite element approximation of the stationary condition is iteratively solved. Our formulation can accommodate arbitrary non‐linear material models for finite deformations. In the examples we present results for isotropic plasticity at finite rotations and small strains as well as bifurcation problems and post‐buckling response. The essential feature of the new element is the robustness in the equilibrium iterations. It allows very large load steps in comparison to other element formulations. Copyright © 2005 John Wiley & Sons, Ltd.     

A new approach to the dynamic analysis of structures using fixed frequency dynamic stiffness matrices

The dynamic analysis of structures by the standard finite element method introduces additional inaccuracies into the solution which are not present when the method is used for static analyses. These inaccuracies can arise from two sources: (i) the element formulation and (ii) the reduction of the size of the matrices by a static condensation (i.e. using the Guyan method^1,2). The errors in both cases are caused by neglecting frequency-dependent terms in the functions relating the displacements at any point in the structure to the displacements at certain fixed points (i.e. nodes in the element formulation and ‘masters’ in the condensation). A new method of solution is proposed in this paper in which the frequency-dependent terms are retained implicitly by using dynamic stiffness matrices defined at a number of fixed frequencies. The dynamic stiffness matrices may be condensed efficiently to a relatively small number of master degrees-of-freedom using a front solution algorithm. The final stage in the solution uses these matrices to synthesize a high-order eigenvalue problem. A method of solving such an eigenvalue problem, of arbitrary order, is described in a separate paper.³ Numerical examples are given to show the accuracy and efficiency of the proposed method compared with conventional methods of solution.     

Non‐linear random response of laminated composite shallow shells using finite element modal method

This paper investigates the large‐amplitude multi‐mode random response of thin shallow shells with rectangular planform at elevated temperatures using a finite element non‐linear modal formulation. A thin laminated composite shallow shell element and the system equations of motion are developed. The system equations in structural node degrees‐of‐freedom (DOF) are transformed into modal co‐ordinates, and the non‐linear stiffness matrices are transformed into non‐linear modal stiffness matrices. The number of modal equations is much smaller than the number of equations in structural node DOF. A numerical integration is employed to determine the random response. Thermal buckling deflections are obtained to explain the intermittent snap‐through phenomenon. The natural frequencies of the infinitesimal vibration about the thermally buckled equilibrium positions (BEPs) are studied, and it is found that there is great difference between the frequencies about the primary (positive) and the secondary (negative) BEPs. All three types of motion: (i) linear random vibration about the primary BEP, (ii) intermittent snap‐through between the two BEPs, and (iii) non‐linear large‐amplitude random vibration over the two BEPs, can be predicted. Copyright © 2006 John Wiley & Sons, Ltd.     

Free vibration analysis of arches using curved beam elements

The natural frequencies and mode shapes for the radial (in‐plane) bending vibrations of the uniform circular arches were investigated by means of the finite arch (curved beam) elements. Instead of the complicated explicit shape functions of the arch element given by the existing literature, the simple implicit shape functions associated with the tangential, radial (or normal) and rotational displacements of the arch element were derived and presented in matrix form. Based on the relationship between the nodal forces and the nodal displacements of a two‐node six‐degree‐of‐freedom arch element, the elemental stiffness matrix was derived, and based on the equation of kinetic energy and the implicit shape functions of an arch element the elemental consistent mass matrix with rotary inertia effect considered was obtained. Assembly of the foregoing elemental property matrices yields the overall stiffness and mass matrices of the complete curved beam. The standard techniques were used to determine the natural frequencies and mode shapes for the curved beam with various boundary conditions and subtended angles. In addition to the typical circular arches with constant curvatures, a hybrid beam constructed by using an arch segment connected with a straight beam segment at each of its two ends was also studied. For simplicity, a lumped mass model for the arch element was also presented. All numerical results were compared with the existing literature or those obtained from the finite element method based on the conventional straight beam element and good agreements were achieved. Copyright © 2003 John Wiley & Sons, Ltd.     

Isogeometric analysis of fluid‐saturated porous media including flow in the cracks

An isogeometric model is developed for the analysis of fluid transport in pre‐existing faults or cracks that are embedded in a fluid‐saturated deformable porous medium. Flow of the interstitial fluid in the porous medium and fluid transport in the discontinuities are accounted for and are coupled. The modelling of a fluid‐saturated porous medium in general requires the interpolation of the displacements of the solid to be one order higher than that of the pressure of the interstitial fluid. Using order elevation and Bézier projection, a consistent procedure has been developed to accomplish this in an isogeometric framework. Particular attention has also been given to the spatial integration along the isogeometric interface element in order to suppress traction oscillations that can arise for certain integration rules when a relatively high dummy stiffness is used in a poromechanical model. © 2016 The Authors. International Journal for Numerical Methods in Engineering Published by John Wiley & Sons Ltd.     

Coupled FEM–DBEM approach on multiple crack growth in cryogenic magnet system of nuclear fusion experiment ‘Wendelstein 7‐X’

At the Max Planck Institute for plasma physics in Greifswald, Germany, the world's largest nuclear fusion experiment of modular stellarator type Wendelstein 7‐X has started plasma operation. The hot hydrogen plasma is confined in a plasma vessel by an electromagnetic field generated by 50 non‐planar and 20 planar superconducting coils. The superconducting coils are encased in cast stainless steel housings. The coils are bolted onto a central support ring and welded together by so called lateral support elements (LSEs). In this paper, a procedure, based on a global–local finite element method (FEM)–dual boundary element method (DBEM) approach, is developed to simulate the propagation of multiple cracks detected in LSEs and undergoing a fatigue load spectrum. The global stress analysis on the superconducting coils is performed by FEM whereas the sub‐modelling approach is adopted to solve the crack propagation in the DBEM environment. The boundary conditions applied on the DBEM submodel are the displacements calculated by the FEM global analysis, in correspondence of the cut surfaces (there are no body forces nor external loads applied on the submodel volume). Two cracks are simultaneously introduced, and a linear elastic fracture mechanics analysis is performed. Results in terms of cracks growth rates and evolving crack shapes are provided, and the residual life of the component is forecast.