Solid elements with rotational degrees of freedom: Part 1—hexahedron elements

This is the first of a two part paper on three-dimensional finite elements with rotational degrees of freedom (DOF). Part I introduces an 8-node solid hexahedron element having three translational and three rotational DOF per node. The corner rotations are introduced by transformation of the midside translational DOF of a 20-node hexahedron element. The new element produces a much smaller effective band width of the global system equations than does the 20-node hexahedron element having midside nodes. A small penalty stiffness is introduced to augment the usual element stiffness so that no spurious zero energy modes are present. The new element passes the patch test and demonstrates greatly improved performance over elements of identical shape but having only translational DOF at the corner nodes.     

Solid elements with rotational degrees of freedom: Part II—tetrahedron elements

This is the second part of a two part paper on three-dimensional finite elements with rotational degrees of freedom (DOF). Part II introduces a solid tetrahedron element having 3 translational and 3 rotational DOF per node. The corner rotations are introduced by transformation of the midside translational DOF of a 10-node tetrahedron element. To further enhance the element performance a least squares strain extraction technique is also implemented to develop the stiffness matrix with a desired field. The strain smoothing improves performance without causing a loss in generality. As with the hexahedron in Part I, the element stiffness is augmented with a small penalty stiffness to eliminate any possible spurious zero energy modes. The new tetrahedron element passes the patch test and demonstrates much improved performance over the 4-node translational DOF only (constant strain) tetrahedron element.     

Modified and Trefftz unsymmetric finite element models

The unsymmetric finite element method employs compatible test functions but incompatible trial functions. The pertinent 8-node quadrilateral and 20-node hexahedron unsymmetric elements possess exceptional immunity to mesh distortion. It was noted later that they are not invariant and the proposed remedy is to formulate the element stiffness matrix in a local frame and then transform the matrix back to the global frame. In this paper, a more efficient approach will be proposed to secure the invariance. To our best knowledge, unsymmetric 4-node quadrilateral and 8-node hexahedron do not exist. They will be devised by using the Trefftz functions as the trial function. Numerical examples show that the two elements also possess exceptional immunity to mesh distortion with respect to other advanced elements of the same nodal configurations.     

Alternative reduced integration avoiding spurious modes for 8-node quadrilateral and 20-node hexahedron finite elements

Integrating the isoparametric 8-node quadrilateral and the 20-node hexahedron elements with Gauss integration based on the 3 point rule produces stiff elements. The excessive stiffness is mainly due to locking phenomenon. One remedy to partly remove locking consists in using reduced integration. Mostly, 2 × 2 or 2 × 2 × 2 integration, respectively, is employed. The lower order integration introduces spurious element modes, however. These modes may deteriorate solutions for finite element models. To overcome this drawback alternative reduced integration procedures are presented. A 5-point rule for the quadrilateral is described. 9-point and 21-point procedures are introduced for the hexahedron. The performance of these procedures is studied by some test problems.     

带转角自由度的四边形六面体单元

本文提出了三个带转角自由度单元，其中一个平面四边形单元，两个空间六面体单元。对平面单元每个结点有两个线位移自由度、一个转角自由度；对空间单元，每个结点有三个线位移自由度、三个转角自由度。这些单元列式简单，其中两个无多余零能模式，数值计算表明，它们的计算精度高。     

Decohesion finite element with enriched basis functions for delamination

In order to address the stringent mesh size requirement of the cohesive element, currently 0.5 mm or less, the mixed-mode interface finite element is enriched with the analytical solution of an idealized beam on elastic foundation. Both the interface and the solid continuum elements are enriched; the partition of unity (PU) method is utilized to obtain the enhanced interpolation, implemented with two user elements (UEL) in the commercial package ABAQUS/Standard. The new formulation yields predictions in excellent agreement with theoretical and experimental results for a typical mode I delamination benchmark, with elements 5 mm in size.     

Modelling delamination onset and growth in pin loaded composite laminates

A three-dimensional (3D) finite element (FE) model is created with cohesive zone elements (CZE) to simulate a mechanically fastened [0°/90°]_s pin-loaded joint in a composite laminate. The model incorporates fully integrated solid elements in the pin-loaded area to accurately capture the high stress gradients. Contact based cohesive elements with a bilinear traction–separation law are inserted between the layers to capture the onset and growth of delamination. The stress distribution around the pin-loaded hole was verified with the widely used cosine stress distribution model. Results from the FE model show that delamination damage initiated at the point of maximum average shear stress at the 0°/90° interface. The delaminated area develops an elliptical shape which grows in a non-self similar manner with increasing pin displacement. It is concluded that a progressive damage model should be included to provide a full understanding of the failure sequence, work that the authors are currently engaged with.     

On improved hybrid finite elements with rotational degrees of freedom

A set of three new hybrid elements with rotational degrees-of-freedom (d.o.f.'s) is introduced. The solid, 8-node, hexahedron element is developed for solving three-dimensional elasticity problems. This element has three translational and three rotational d.o.f.'s at each node and is based on a 42-parameter. three-dimensional stress field in the natural convected co-ordinate system. For two-dimensional, plane elasticity problems, an improved triangular hybrid element and a quadrilateral hybrid element are presented. These elements use two translational and one rotational d.o.f. at each node. Three different sets of five-parameter stress fields defined in a natural convected co-ordinate system for the entire element are used for the mixed triangular element. The mixed quadrilateral element is based on a nine-parameter complete linear stress field in natural space. The midside translational d.o.f.'s are expressed in terms of the corner nodal translations and rotations using appropriate transformations. The stiffness matrix is derived based on the Hellinger–Reissner variational principle. The elements pass the patch test and demonstrate an improved performance over the existing elements for prescribed test examples.     

On the robustness of finite element procedures based on Virtual Crack Closure Technique and fail release approach for delamination growth phenomena. Definition and assessment of a novel methodology

Numerical procedures based on the combined use of the Virtual Crack Closure Technique and of a fail release approach have been widely used to simulate delamination growth phenomena of composite material structures. This paper starts explaining why this kind of methodologies might not be robust due to mesh and load step size dependency and introduces a novel approach able to cope with the problems identified. Finally the effectiveness and robustness of the proposed procedure, implemented into a commercial finite element software by means of user subroutines, are assessed by comparing the obtained numerical results for a delamination growth phenomenon against literature experimental data on a stiffened panel with a circular embedded delamination under compressive load.     

Element by element a posteriori error estimation of the finite element analysis for three-dimensional elastic problems

An a posteriori error estimation method for finite element solutions for three-dimensional elastic problems is presented based on the theory developed by the authors for two-dimensional problems.¹ The error is estimated for the finite element solutions obtained using three-dimensional 8-node elements with a linear interpolation function in an arbitrary hexahedron. The method is successfully applied to three-dimensional elastic problems. In order to decrease computing time and memory use, the error is estimated element by element. The major difficulty in the element-wise error estimation technique is satisfying the self-equilibrium condition of applied forces, especially in three-dimensional problems. These forces are mainly due to traction discontinuity on the element boundaries. The difficulty is circumvented by employing an element-wise optimal procedure. It is also shown that a very accurate stress solution can be obtained by adding estimated error to the original finite element solutions.     

Formulation and assessment of an enhanced finite element procedure for the analysis of delamination growth phenomena in composite structures

An existing procedure based on the combined use of the Virtual Crack Closure Technique and of a fail release approach for the analysis of delamination growth phenomena in composite structures has been enhanced with a front-tracing algorithm and suitable expressions for the evaluation of the Strain Energy Release Rate when dealing with non-smoothed delamination fronts. The enhanced procedure has been implemented into a commercial finite element software by means of user subroutines and applied to the analysis of a composite stiffened panel with an embedded delamination under compressive load. The effectiveness and robustness of the enhanced procedure have been assessed by comparing literature experimental data and numerical results obtained by using different mesh densities in the damaged area (global/local approach).     

XFEM simulation of delamination in composite laminates

In this paper, the extended finite element method (XFEM) is extended to simulate delamination problems in composite laminates. A crack-leading model is proposed and implemented in the ABAQUS® to discriminate different delamination morphologies, i.e., the 0°/0° interface in unidirectional laminates and the 0°/90° interface in multidirectional laminates, which accounts for both interlaminar and intralaminar crack propagation. Three typical delamination problems were simulated and verified. The results of single delamination in unidirectional laminates under pure mode I, mode II, and mixed mode I/II correspond well with the analytical solutions. The results of multiple delaminations in unidirectional laminates are in good agreement with experimental data. Finally, using a recently proposed test that characterizes the interaction of delamination and matrix cracks in cross-ply laminates, the present numerical results of the delamination migration caused by the coupled failure mechanisms are consistent with experimental observations.     

Experimental investigation and finite element modeling of mixed-mode delamination in a moisture-exposed carbon/epoxy composite

An investigation of the effects of moisture on mixed-mode I/II delamination growth in a carbon/epoxy composite is presented. Experimental quasi-static and fatigue delamination tests were carried out on composite specimens. The quasi-static fracture test results showed that exposure to moisture led to a decrease in mode II and mixed-mode delamination toughness while mode I toughness was enhanced. The fatigue tests revealed an adverse effect of moisture on delamination growth under mixed-mode loadings. Existing delamination criteria and growth rate models were evaluated to determine which ones best predict delamination toughness and growth, respectively, at any given mixed-mode ratio. Quasi-static and fatigue simulations with a cohesive zone-based finite element model that incorporated the selected mixed-mode delamination models were performed and good agreement between experimental and numerical data was shown for dry and moisture-exposed specimens.     

Predicting self-healing strength recovery using a multi-objective genetic algorithm

This study aims to identify the optimum location and distribution of a healing agent within the delamination of a fibre reinforced plastic to ensure effective self-healing by utilising a multi-objective Genetic Algorithm (GA). Two optimisation problems were formulated and addressed with a different set of objectives. A simple finite element (FE) model is used to evaluate the mechanical performance of the healing component. The FE model consists of an idealised delamination region, which allows the direct discretisation of the problem used for the optimisation algorithm. Effective healing locations are found for a specific load case with a healing efficiency of up to 95% for the best performing solution.     

Panel stiffener debonding analysis using a shell/3D modeling technique

A shear loaded, stringer reinforced composite panel is analyzed to evaluate the fidelity of computational fracture mechanics analyses of complex structures. Shear loading causes the panel to buckle. The resulting out-of-plane deformations initiate skin/stringer separation at the location of an embedded defect. The panel and surrounding load fixture were modeled with shell elements. A small section of the stringer foot, web and noodle as well as the panel skin near the delamination front were modeled with a local 3D solid model. Across the width of the stringer foot, the mixed-mode strain energy release rates were calculated using the virtual crack closure technique. A failure index was calculated by correlating the results with a mixed-mode failure criterion of the graphite/epoxy material. The objective was to study the effect of the fidelity of the local 3D finite element model on the computed mixed-mode strain energy release rates and the failure index.     

Analysis of Mode I and Mode II Crack Growth Arrest Mechanism with Z-Fibre Pins in Composite Laminated Joint

This paper presents the numerical study of the mode I and mode II interlaminar crack growth arrest in hybrid laminated curved composite stiffened joint with Z-fibre reinforcement. A FE model of hybrid laminated skin-stiffener joint reinforced with Z-pins is developed to investigate the effect of Z- fibre pins on mode I and mode II crack growth where the delamination is embedded inbetween the skin and stiffener interface. A finite element model was developed using S4R element of a 4-node doubly curved thick shell elements to model the composite laminates and non linear interface elements to simulate the reinforcements. The numerical analyses revealed that Z-fibre pinning were effective in suppressing the delamination growth when propagated due to applied loads. Therefore, the Z-fibre technique effectively improves the crack growth resistance and hence arrests or delays crack growth extension.     

实体退化板单元及其在板的振动分析中的应用

经典板壳单元是由板壳理论构造出来的,而经典的板壳理论是在空间弹性理论的基础上考虑板壳的基本假定得来的。在空间等参数单元的基础上,直接引入板壳的基本假定,修改空间等参数单元的弹性矩阵,从而构造出适合于厚薄板壳分析的20结点实体退化板单元,并将其应用于开口圆柱薄壳的静力分析和厚薄板的固有振动分析。数值算例表明,该单元收敛快,稳定性好,具有较高的精度。此外,该单元还可以用于曲边变厚度板、壳体及层合板的振动分析。     

Environmental effects on mode II fatigue delamination growth in an aerospace grade carbon/epoxy composite

An investigation of the effects of water, hydraulic fluid and deicing fluid exposure on mode II delamination propagation in an aerospace grade composite is presented. All exposed specimens suffered a loss in delamination toughness and an increase in fatigue delamination growth rate, which was particularly significant for deicing fluid exposure. The number of cycles for delamination onset was also reduced by these exposures, although scanning electron micrographs showed no significant differences between the fracture pattern of dry and exposed specimens. It was also shown that environmental effects can be simply accounted for in a cohesive zone based finite element model.     

Mode III fatigue delamination growth of glass fiber reinforced polymer woven laminates at cryogenic temperatures

This paper investigates the cryogenic fatigue delamination behavior of glass fiber reinforced polymer woven laminates under Mode III loading. Fatigue delamination tests were conducted using split cantilever beam specimens at room temperature, liquid nitrogen temperature (77 K) and liquid helium temperature (4 K). A finite element analysis was also employed to calculate the energy release rate. The temperature dependence of the fatigue delamination growth rate vs. energy release rate range is discussed. Fracture surfaces were examined by scanning electron microscopy to identify the delamination mechanisms under fatigue loading. The important conclusion we reach is that the Mode III fatigue delamination growth rates of woven laminates at cryogenic temperatures are lower than that at room temperature.     

A cohesive zone model for predicting delamination suppression in z-pinned laminates

This paper presents a cohesive zone model based finite element analysis of delamination resistance of z-pin reinforced double cantilever beam (DCB). The main difference between this and existing cohesive zone models is that each z-pin bridging force is governed by a traction-separation law derived from a meso-mechanical model of the pin pullout process, which is independent of the fracture toughness of unreinforced laminate. Therefore, two different traction-separation laws are used: one representing the toughness of unreinforced laminate and the other the enhanced delamination toughness owing to the pin bridging action. This approach can account for the large scale bridging effect and avoid using concentrated pin forces, thus removing the mesh dependency and permitting more accurate analysis solution. Computations were performed using a simplified unit strip model. Predicted delamination growth and load vs. displacement relation are in excellent agreement with the prediction by a complete model, and both models are in good agreement with test measured load vs. displacement relation. For a pinned DCB specimen, the unit strip model can reduce the computing time by 85%.