An enriched element‐failure method (REFM) for delamination analysis of composite structures

This paper develops an enriched element‐failure method for delamination analysis of composite structures. This method combines discontinuous enrichments in the extended finite element method and element‐failure concepts in the element‐failure method within the finite element framework. An improved discontinuous enrichment function is presented to effectively model the kinked discontinuities; and, based on fracture mechanics, a general near‐tip enrichment function is also derived from the asymptotic displacement fields to represent the discontinuity and local stress intensification around the crack‐tip. The delamination is treated as a crack problem that is represented by the discontinuous enrichment functions and then the enrichments are transformed to external nodal forces applied to nodes around the crack. The crack and its propagation are modeled by the ‘failed elements’ that are applied to the external nodal forces. Delamination and crack kinking problems can be solved simultaneously without remeshing the model or re‐assembling the stiffness matrix with this method. Examples are used to demonstrate the application of the proposed method to delamination analysis. The validity of the proposed method is verified and the simulation results show that both interlaminar delamination and crack kinking (intralaminar crack) occur in the cross‐ply laminated plate, which is observed in the experiment. Copyright © 2009 John Wiley & Sons, Ltd.     

3D‐based hp‐adaptive first‐order shell finite element for modelling and analysis of complex structures—Part 2: Application to structural analysis

The paper presents a 3D‐based adaptive first‐order shell finite element to be applied to hierarchical modelling and adaptive analysis of complex structures. The main feature of the element is that it is equipped with 3D degrees of freedom, while its mechanical model corresponds to classical first‐order shell theory. Other useful features of the element are its modelling and adaptive capabilities. The element is assigned to hierarchical modelling and hpq‐adaptive analysis of shell parts of complex structures consisting of solid, thick‐ and thin‐shell parts, as well as of transition zones, where h, p and q denote the mesh density parameter and the longitudinal and transverse orders of approximation, respectively. The proposed hp‐adaptive first‐order shell element can be joined with 3D‐based hpq‐adaptive hierarchical shell elements or 3D hpp‐adaptive solid elements by means of the family of 3D‐based hpq/hp‐ or hpp/hp‐adaptive transition elements. The main objective of the first part of our research, presented in the first part of the paper, was to provide non‐standard information on the original parts of the element algorithm. Here we describe the second part of the research, devoted to the methodology and results of the application of the element to various plate and shell problems. The main objective of this part is to verify algorithms of the element and to show its usefulness in modelling and adaptive analysis of shell and plate parts of complex structures. In order to do that, there is a presentation of the results of a comparative analysis of model plate and shell problems using the classical and our elements, and equidistributed and integrated Legendre shape functions. For the plate problem a comparison of the results obtained from the adaptive and non‐adaptive analysis is also included. Additionally, some advantages of the application of our element are shown through a comparative analysis of p‐convergence of the thin plate problem and an adaptive analysis of the exemplary complex structure. Copyright © 2006 John Wiley & Sons, Ltd.     

hp‐Generalized FEM and crack surface representation for non‐planar 3‐D cracks

A high‐order generalized finite element method (GFEM) for non‐planar three‐dimensional crack surfaces is presented. Discontinuous p‐hierarchical enrichment functions are applied to strongly graded tetrahedral meshes automatically created around crack fronts. The GFEM is able to model a crack arbitrarily located within a finite element (FE) mesh and thus the proposed method allows fully automated fracture analysis using an existing FE discretization without cracks. We also propose a crack surface representation that is independent of the underlying GFEM discretization and controlled only by the physics of the problem. The representation preserves continuity of the crack surface while being able to represent non‐planar, non‐smooth, crack surfaces inside of elements of any size. The proposed representation also provides support for the implementation of accurate, robust, and computationally efficient numerical integration of the weak form over elements cut by the crack surface. Numerical simulations using the proposed GFEM show high convergence rates of extracted stress intensity factors along non‐planar curved crack fronts and the robustness of the method. Copyright © 2008 John Wiley & Sons, Ltd.     

Arbitrary discontinuities in space–time finite elements by level sets and X‐FEM

An enriched finite element method with arbitrary discontinuities in space–time is presented. The discontinuities are treated by the extended finite element method (X‐FEM), which uses a local partition of unity enrichment to introduce discontinuities along a moving hyper‐surface which is described by level sets. A space–time weak form for conservation laws is developed where the Rankine–Hugoniot jump conditions are natural conditions of the weak form. The method is illustrated in the solution of first order hyperbolic equations and applied to linear first order wave and non‐linear Burgers' equations. By capturing the discontinuity in time as well as space, results are improved over capturing the discontinuity in space alone and the method is remarkably accurate. Implications to standard semi‐discretization X‐FEM formulations are also discussed. Copyright © 2004 John Wiley & Sons, Ltd.     

A computational model for analysing interactive buckling and delamination growth in composite structures

In this paper, a unified method is presented: (i) to model delaminated stiffened laminated composite shells; (ii) for synthesising accurate multiple post-buckling solution paths under compressive loading; and (iii) for predicting delamination growth. A multi-domain modelling technique is used for modelling the delaminated stiffened shell structures. Error-free geometrically nonlinear element formulations — a 2-noded curved stiffener element (BEAM2) and a 3-noded shell element (SHELL3) — are used for the finite element analysis. An accurate and simple automated solution strategy based on Newton type iterations is used for predicting the general geometrically nonlinear and postbuckling behaviour of structures. A simple method derived from the 3-dimensionalJ-integral is used for computing the pointwise energy release rate at the delamination front in the plate/shell models. Finally, the influence of post-buckling structural behaviour and the delamination growth on each other has been demonstrated.     

Delamination modeling in doubly curved laminated shells for free vibration analysis using zigzag theory-based facet shell element and hybrid continuity method

We present a finite element (FE) formulation for the free vibration analysis of doubly curved laminated composite and sandwich shells having multiple delaminations, employing a facet shell element based on the efficient third-order zigzag theory and the region approach of modeling delaminations. The methodology, hitherto not attempted, is general for delaminations occurring at multiple interfacial and spatial locations. A recently developed hybrid method is used for satisfying the continuity of the nonlinear layerwise displacement field at the delamination fronts. The formulation is shown to yield very accurate results with reference to full-field three-dimensional FE solutions, for the natural frequencies and mode shapes of delaminated shallow and deep, composite and highly inhomogeneous soft-core sandwich shells of different geometries and boundary conditions, with a significant computational advantage. The accuracy is sensitive to the continuity method used at the delamination fronts, the usual point continuity method yielding rather poor accuracy, and the proposed hybrid method giving the best accuracy. Such efficient modeling of laminated shells with delamination damage will be of immense use for model-based techniques for structural health monitoring of laminated shell-type structures.     

A conformal decomposition finite element method for arbitrary discontinuities on moving interfaces

Interface capturing methods using enriched finite element formulations are well suited for solving multimaterial transport problems that contain weak or strong discontinuities. The conformal decomposition FEM decomposes multimaterial elements of a non‐conforming background mesh into sub‐elements that conform to material interfaces captured using a level set method. As the interface evolves, interfacial nodes move, and background nodes may change material. The present work describes approaches for handling moving interfaces in the context of the conformal decomposition FEM for both weakly and strongly discontinuous fields. Dynamic discretization methods using extrapolation and moving mesh approaches are considered and developed with first‐order and second‐order time integration methods. The moving mesh approach is demonstrated to be a stable method that preserves both weak and strong discontinuities on a variety of one‐dimensional and two‐dimensional test problems, while achieving the expected second‐order error convergence rate in space and time. Copyright © 2014 John Wiley & Sons, Ltd.     

Adaptive modelling of delamination initiation and propagation using an equivalent single‐layer shell approach

A major challenge for crash failure analysis of laminated composites is to find a modelling approach, which is both sufficiently accurate, for example, able to capture delaminations, and computationally efficient to allow full‐scale vehicle crash simulations. Addressing this challenge, we propose a methodology based on an equivalent single‐layer shell formulation which is adaptively through‐the‐thickness refined to capture initiating and propagating delaminations. To be specific, single shell elements through the laminate thickness are locally and adaptively enriched using the extended finite element method such that delaminations can be explicitly modelled without having to be represented by separate elements. Furthermore, the shell formulation is combined with a stress recovery technique which increases the accuracy of predicting delamination initiations. The paper focuses on the parameters associated with identifying, introducing and extending the enrichment areas; especially on the impact of these parameters on the resulting structural deformation behaviour. We conclude that the delamination enrichment must be large enough to allow the fracture process to be accurately resolved, and we propose a suitable approach to achieve this. The proposed methodology for adaptive delamination modelling shows potential for being computationally efficient, and thereby, it has the potential to enable efficient and accurate full vehicle crash simulations of laminated composites. Copyright © 2017 John Wiley & Sons, Ltd.     

含内埋矩形脱层板屈曲分析的传递函数方法

研究了含任意内埋矩形脱层复合材料层合板的屈曲问题。采用一种基于Mindlin一阶剪切理论的条形传递函数方法,将含内埋矩形脱层的复合材料层合板分成含脱层和不含脱层的两种矩形超级单元,然后由各超级单元之间连接结点处的位移连续和力平衡条件得到脱层板屈曲的特征方程,进而得到脱层板的屈曲载荷和屈曲模态。进行参数分析发现,脱层大小、深度、位置以及脱层板的边界条件和复合材料铺层方向对脱层板屈曲载荷的影响较显著。     

Differential quadrature element method for static analysis of shear deformable cross‐ply laminates

In this paper, the two‐dimensional differential quadrature element method (DQEM) is developed for the static analysis of symmetric cross‐ply laminates using the first‐order shear deformation plate theory. In this study, the laminated plate, which may contain different discontinuities in loading, geometry, material, and boundary conditions, is first divided into several simple plate elements and then the differential quadrature method (DQM) is applied to each simple element. Compatibility conditions are derived to connect the plate elements so that the overall matrix equation system for the whole plate is obtained and solved. The reliability of the DQEM for solving the titled problems is examined carefully through convergence and accuracy studies and finally some numerical test examples are given to demonstrate the applicability and flexibility of this method for practical use. The methodology presented here has overcome some critical drawbacks of the global DQM but is different from the Quadrature Element Method (QEM) since only one grid point is employed to represent the interface point. Copyright © 1999 John Wiley & Sons, Ltd.     

Finite element modeling of delamination by layerwise shell element allowing for interlaminar displacements

This paper presents a novel approach of modeling the delamination phenomenon experienced by laminated composite plate and shell structures by using a previously developed layerwise shell finite element in conjunction with some transformations. This layerwise element is formulated by stacking some single-layered shell elements through a transformation of displacements of the mid-surface of a layer to those on the mid-surface of the laminated composite shell structure. It can accurately model the overall displacements and interlaminar stresses of a laminated composite shell structures whose layers are perfectly (rigidly) bonded. The novelty of the present approach, however, lies in the fact that two different transformations are used so that interlaminar displacements as well as interlaminar stresses can be represented in the finite element model. The transformations allow for displacement mismatches across the normal direction of the layer interfaces (the normal mode of delamination) and between layers (the shear mode). As a result, the proposed methodology can be used to model the open and shear modes of delamination. A two-layered simply supported composite beam and a two-layered simply supported cross-ply square plate are then chosen for numerical studies. These examples demonstrate how the present approach can be applied to accurately model delamination phenomena such as shear slip and normal separation. The paper concludes with suggestions for future work.     

3D‐based hp‐adaptive first‐order shell finite element for modelling and analysis of complex structures—Part 1: The model and the approximation

The paper presents a 3D‐based adaptive first‐order shell finite element to be applied to hierarchical modelling and adaptive analysis of complex structures. The main feature of the element is that it is equipped with 3D degrees of freedom, while its mechanical model corresponds to classical first‐order shell theory. Other useful features of the element are its modelling and adaptive capabilities. The element is assigned to hierarchical modelling and hpq‐adaptive analysis of shell parts of complex structures consisting of solid, thick‐ and thin‐shell parts, as well as of transition zones, where h, p and q denote the mesh density parameter and the longitudinal and transverse orders of approximation, respectively. The proposed hp‐adaptive first‐order shell element can be joined with 3D‐based hpq‐adaptive hierarchical shell elements or 3D hpp‐adaptive solid elements by means of the family of 3D‐based hpq/hp‐ or hpp/hp‐adaptive transition elements. The main objective of the first part of our research, presented in this paper, is to provide non‐standard information on the original parts of the element algorithm. In order to do that, we present the definition of shape functions necessary for p‐adaptivity, as well as the procedure for imposing constraints corresponding to the lack of elongation of the straight lines perpendicular to the shell mid‐surface, which is the procedure necessary for q‐adaptivity. The 3D version of constrained approximation presented next is the basis for h‐adaptivity of the element. The second part of our research, devoted to methodology and results of the numerical research on application of the element to various plate and shell problems, are described in the second part of this paper. Copyright © 2006 John Wiley & Sons, Ltd.     

Adaptive methods for thermomechanical coupled contact problems

In this paper an adaptive method for the analysis of thermomechanical coupled multi‐body contact problems is presented. The method is applied to non‐linear elastic solids undergoing finite (thermal) deformations. The contact model considers non‐linear pressure‐dependent heat flux as well as frictional heating in the interface. A time–space‐finite element discretization of the governing equations is formulated including unilateral constraints due to contact. A staggered solution algorithm has been constructed that allows an independent spatial discretization of the coupled subproblems. A posteriori projection‐based error estimators, which enforce implicitly the special boundary conditions due to thermal contact, are used to control the spatial discretization as well as the adaptive time stepping. Numerical examples are presented to corroborate the applicability of the adaptive algorithm to the considered problem type. Copyright © 2004 John Wiley & Sons, Ltd.     

Finite volume analysis of laminated composite plates

A numerical procedure for analysis of general laminated plates under transverse load is developed utilizing the Mindlin plate theory, the finite volume discretization, and a segregated solution algorithm. The force and moment balance equations with the laminate constitutive relations are written in the form of a generic transport equation. In order to obtain discrete counterparts of the governing equations, the plate is subdivided into N control volumes by a Cartesian numerical mesh. As a result, five sets of N linear equations with N unknowns are obtained and solved using the conjugate gradient method with preconditioning. For the method validation, a number of test cases are designed to cover thick and thin laminated plates with aspect ratio (width to thickness) from 4 to 100. Simply supported orthotropic, symmetric cross‐ply, and angle‐ply laminated plates under uniform and sinusoidal pressure loads are solved, and results are compared with available analytical solutions. The shear correction factor of 5/6 is utilized throughout the procedure, which is consistent with test cases used in the reviewed literature. Comparisons of the finite volume method results for maximum deflections at the center of the plate and the Navier solutions obtained for aspect ratios 10, 20, and 100 shows a very good agreement. Copyright © 2016 John Wiley & Sons, Ltd.     

Buckling analysis of rectangular composite plates with rectangular cutout subjected to linearly varying in-plane loading using fem

A numerical study is carried out using finite element method, to examine the effects of square and rectangular cutout on the buckling behavior of a sixteen ply quasi-isotropic graphite/epoxy symmetrically laminated rectangular composite plate [0°/+45°/-45°/90°]_2s, subjected to various linearly varying in-plane compressive loads. Further, this paper addresses the effects of size of square/rectangular cutout, orientation of square/rectangular cutout, plate aspect ratio(a/b), plate length/thickness ratio(a/t), boundary conditions on the buckling bahaviour of symmetrically laminated rectangular composite plates subjected to various linearly varying in-plane compressive loading. It is observed that the various linearly varying in-plane loads and boundary conditions have a substantial influence on buckling strength of rectangular composite plate with square/rectangular cutout.     

Three-dimensional stress analysis of rotating composite beams due to material discontinuities

Material discontinuity could cause in-plane stress gradients that it arises out-of-plane stresses in regions of sudden transition of material properties. A layerwise laminated plate theory is adapted to laminated beams to analyze analytically the three-dimensional stress field at material discontinuities in rotating composite beams. Equations of motion are obtained by using Hamilton’s principle. The beam is divided into two regions with different layups which are joined together to model the region of material discontinuity. The predicted stress distributions at the ply interfaces are shown to be in good agreement with comparative three-dimensional finite element analysis.     

A design study into the delamination behaviour of tapered composites

A fracture mechanics based analysis has been used to predict the tensile delamination load of tapered laminated plates. Simple laminate examples are used to show the effect of dropped ply thickness, number of delaminating surfaces, and dropped ply axial stiffness on the delamination load. Using these trends and acknowledged guidelines, a design is presented for a complex tapered plate with a view to maximising the onset of delamination.     

Adaptive Kronrod–Patterson integration of non‐linear finite element matrices

Efficient simulation of unsaturated moisture flow in porous media is of great importance in many engineering fields. The highly non‐linear character of unsaturated flow typically gives sharp moving moisture fronts during wetting and drying of materials with strong local moisture permeability and capacity variations as result. It is shown that these strong variations conflict with the common preference for low‐order numerical integration in finite element simulations of unsaturated moisture flow: inaccurate numerical integration leads to errors that are often far more important than errors from inappropriate discretization. In response, this article develops adaptive integration, based on nested Kronrod–Patterson–Gauss integration schemes: basically, the integration order is adapted to the locally observed grade of non‐linearity. Adaptive integration is developed based on a standard infiltration problem, and it is demonstrated that serious reductions in the numbers of required integration points and discretization nodes can be obtained, thus significantly increasing computational efficiency. The multi‐dimensional applicability is exemplified with two‐dimensional wetting and drying applications. While developed for finite element unsaturated moisture transfer simulation, adaptive integration is similarly applicable for other non‐linear problems and other discretization methods, and whereas perhaps outperformed by mesh‐adaptive techniques, adaptive integration requires much less implementation and computation. Both techniques can moreover be easily combined. Copyright © 2009 John Wiley & Sons, Ltd.     

Sensitivity analysis with plate and shell finite elements

This paper addresses the problem of calculating sensitivity data by direct methods for isoparametric plate or shell elements. Sensitivity parameters of interest include intrinsic properties such as material modulus and plate thickness, as well as geometry variables which influence the size and shape of a structure. The sensitivity calculation therefore must consider the parametric mapping within an element, as well as the influence of geometric variables on the orientation of an element in space. The methods presented specialize directly to continuum elements, in which the co-ordinate transformation is omitted, or to simple structural members situated arbitrarily in space. Numerical examples are presented which illustrate the accuracy of the proposed techniques, and the effect of discretization error on computed sensitivities.     

Failure analysis of multiple delaminated composite plates due to bending and impact

The present work aims at the first ply failure analysis of laminated composite plates with arbitrarily located multiple delaminations subjected to transverse static load as well as impact. The theoretical formulation is based on a simple multiple delamination model. Conventional first order shear deformation is assumed using eight-noded isoparametric quadratic elements to develop the finite element analysis procedure. Composite plates are assumed to contain both single and multiple delaminations. For the case of impact, Newmark time integration algorithm is employed for solving the time dependent multiple equations of the plate and the impactor. Tsai-Wu failure criterion is used to check for failure of the laminate for both the cases. To investigate the first ply failure, parametric studies are made for different cases by varying the size and number of delaminations as well as the stacking sequences and boundary conditions.