Abstract

Although laminated materials have been used for decades, their employment has increased nowadays in the last years as a result of the gained confidence of the industry on these materials. This has provided the scientific community many reasons to dedicate considerable amount of time and efforts to address a better understanding of their mechanical behavior. With this objective both, experimental and numerical simulation have been working together to give response to a variety of problems related with these materials.

Regarding numerical simulation, a correct modeling of the kinematics of laminated materials is essential to capture the real behavior of the structure. Moreover, once the kinematics of the structure has been accurately predicted other non-linear phenomena such as damage and/or plasticity process could be also studied.

In consequence, in order to contribute to the constant development of simpler and more efficient numerical tools to model laminated materials, a numerical method for modeling mode II/III delamination in advanced composite materials using one-and two-dimensional finite elements is proposed inthis work. In addition, two finite elements base on a zigzag theory for simulating higly heterogeneous multilayered beams and plates structures are developed here.

The document is written based on results of four papers published in indexed journals. Copies of all these papers are included in Appendix. The main body of this thesis is constituted by Chapters 2 to 4. Chapter 2 deals with the numerical treatment of laminated beams and plates. Chapter 3 presents the formulation of the LRZ beam and the QLRZ plate finite elements based on the Refined Zigzag theory. Finally, the main contribution of this thesis, the LRZ/QLRZ delamination model, is developed in Chaptert 4.

Abstract

This work shows possibilities and limitations of the refined zigzag theory (RZT) that has been used in different structural (beam, plate and shell) finite elements. The refined zigzag theory can deal with composite laminates, adding only one nodal degree of freedom per spatial dimension of the laminate, obtaining very good accuracy. It assumes that the in-plane displacements have a piece-wise linear shape across the thickness depending on the shear stiffness of each composite layer. This paper presents the main aspects of a beam/shell of revolution element used for the numerical simulations. The details of the refined zigzag theory are given also in order to discuss some limitations that occur when dealing with the non-linear phenomenon of delamination. Two examples are presented and discussed, including different inhomogeneities that show the limitations of the RZT for the treatment of partially delaminated beams

Abstract

A method based on the Refined Zigzag Theory (RZT) to model delamination in composite laminated beam structures is presented. The novelty of this method is the use of one-dimensional finite elements to discretize the geometry of the beam. The key property of this beam element, named LRZ, is the possibility to capture the relative displacement between consecutive layers which occurs during delamination. The fracture mode that the LRZ element is capable to predict is mode II. In order to capture the relative displacement using the LRZ element it is necessary to adapt the RZT theory as presented in this paper. The mechanical properties of the layers are modeled using a continuum isotropic damage model. The modified Newton–Raphson method is used for solving the non-linear problem.

The RZT theory, the LRZ finite element and the isotropic damage model are described in the paper. Also, the implicit integrations algorithm is presented. The performance of the LRZ element is analyzed by studying the delamination in a beam for two different laminates, using the plane stress solution as a reference.

Abstract

We present a new 2-noded beam element based on the refined zigzag theory and the classical Euler–Bernoulli beam theory for the static analysis of composite laminate and sandwich beams. The proposed element is able to take into account distortion effects due to shear elastic strains and can predict delamination. The element has four degrees of freedom per node. A C1 cubic Hermite interpolation is used for the vertical deflection while a C0 linear interpolation is employed for the other kinematics variables. The stiffness matrix and the load vector are calculated in explicit form using exact integration. The element is free from shear locking as confirmed with numerical tests on a wide range of the slenderness ratios. Numerical results show the ability of the EEBZ2 element to reproduce accurately the vertical deflection along the beam length and complex zigzag distributions of the axial displacement and the stresses across the thickness. Delamination effects are modeled by incorporating of an additional zigzag function corresponding to the kinematics of a zero thickness layer where delamination occurs. An example showing the capability of the new EEBZ2 element for accurately reproducing delamination effects is presented.

Abstract

In this work a kinematics for laminated beams enriched with a refined formulation ZigZag (RZT), originally presented by Tessler et al. in 2007, introduced in a hierarchical one dimensional type “p” finite element is presented. The finite element employs Lagrange polynomials for the approximation of the degrees of freedom of the ends (nodes) and orthogonal Gram-Schmidt polynomials to the internal degrees of freedoms. This finite element allows a very low discretization, is free of shear locking and behaves very well when the analysis of laminated composites with accurate determination of local stresses and strains at laminar level is necessary.

This element has been validated in the analysis of laminated beams with various sequences of symmetric and asymmetric stacking, studying in each case its accuracy and stability.