Finite Element Analysis of Composite Laminated

In the research for sustainable construction, cross-laminated timber (CLT) has gained popularity and become a widely used engineered timber product. However, there are few numerical studies of the structural behaviour of CLT. Among other issues, the orthotropic properties of CLT complicate finite element analysis (FEA). This paper presents a finite element model (FEM) to predict the structural behaviour of CLT beams subjected to sustained flexural loading. This numerical model includes a material model based on the orthotropic material properties of different timber species. Furthermore, the orientation and the properties of each layer are considered. Most of the previous studies simulate CLT beams as a homogeneous material. However, in this work the CLT beam is modelled as a composite material made up of five layers with different orientations and properties. Bonded contacts are used to define the interaction between layers. In addition, nonlinearities, such as large displacement, are used to simulate the behaviour of CLT beams. The model provides the loaddisplacement relationship and stress concentration. Tsai-Wu failure criteria is used in the simulation to predict the failure modes of the CLT beams studied.


Introduction
Cross-laminated timber (CLT) beams have become a popular structural element. In the early 1990s, CLT was developed in Austria, Germany and Switzerland. The advantages of CLT include shorter building time and a high strength ratio [1]. CLT beams are composed of wood layers, arranged and glued crosswise as illustrated in Figure 1. CLT beams are commonly made of an odd number of lamellas which form a solid panel. The alternating grain orientation, laminar structure, and capacity to bear in-and out-of-plane loads [2] make it suitable as a wall and floor element.
Due to its orthogonal grain direction and to the anisotropy of wood, the mechanical behaviour of CLT beams is complex. The elastic mechanical properties are different along principal axes, which are the axial direction (grain), the radial direction, and the circumferential direction (i.e., strength and stiffness in the axial direction are larger than those in radial and circumferential direction) [3]. Therefore, to design CLT elements (walls and slabs) the mechanical properties must be understood. Various analytical methods are used to study the mechanical properties and the mechanical behaviour of CLT for both in-and out-of-plane loads [4,5]. Some authors consider CLT as a laminated composite material. Blass [6] developed a methodology for designing CLT panels based on the composite theory. Most of the existing studies are based on 2D plate theories [7]. At present, research is aimed to develop 3D analyses to better understand CLT behaviour.
This study aims to develop a finite element method (FEM) to predict the deformation and the failure mechanism of a CLT panel under perpendicular loads to plane. Experiments conducted in the University of Edinburgh at ambient temperature [8] were used to compare the finite element (FE) results.

Finite Element Analysis (FEA)
In order to compare and verify the accuracy of the FE model developed, the boundary conditions at ambient temperature applied in the model were the same as in work [8]. A 3D FE model of one of the CLT beams studied by [8] was developed in ANSYS 18.1.
This FE model analyzed the structural response of a CLT beam under a sustained load. The CLT beam used to develop the model has a 5-layer configuration. The mechanical properties of timber were taken from previous research [8].
Considering the CLT beam as a composite element, it was modelled using the composite module of ANSYS [9]. The FE model developed used solid elements since each of the plies was considered as a thick element. Elastic material properties were assigned in each orthogonal direction of the layer. Other assumptions are: layers are perfectly bonded together and the material properties of each layer are constant through the thickness.
The configuration of the model in different layers was the same as the real configuration. The 5layer configuration used was: 20(l) + 20(p) + 20(l) + 20(p) + 20(l) mm. Where l and p means longitudinal and perpendicular grain direction, as shown in Figure 2.
The CLT model was 1000 mm long, 300 mm wide, and 100 mm thick, as shown in Figure 3. In order to reduce the computational cost, a simplified model using conditions of symmetry is used, (Figure 4). The symmetry conditions were applied in the longitudinal direction. These conditions reduce the size of the numerical model.  Each ply was modelled using 8-node solid elements (SOLID 186). The main features of the mesh for the model are shown in Table 1. Boundary conditions applied in this simulation are the same as in experimental tests. These conditions are the following:
Loading applied at the top surface with a roller.  The Tsai Wu failure criterion was used in ACP-Post to evaluate the margin to failure. This criterion is useful to evaluate mechanical behaviour in timber elements due to its different strengths in compression and tension [10].

Results
The numerical and the experimental results are compared to validate the numerical model. The displacement applied to the beam generates a reaction load ( Figure 5). This load is equal to the load applied in the experimental tests. The ACP-Post module plots the stresses and failure modes of the numerical beam model. As this model was developed with solid elements, the stress is plotted at the top and bottom of each layer.
The stress distribution is shown in Figure 6, which has the same form as the stress diagrams for CLT beams in other works [2,4].

Conclusions
This work represents the first step to analyze the structural behaviour of a CLT beam by using finite elements. This procedure uses ANSYS Composite PrePost and provides the tools to study CLT beam as a composite element. Furthermore, the flexibility of the model in the design of the fibre orientation of each layer makes it possible to assess the ultimate limit state of the CLT beams.