• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.2
• /
• pp.354-362
• /
• 1995

The finite element modeling is used to study the buckling and postbuckling behavior of composite laminates with an embedded delamination. Degenerated shell element and rigid beam element are utilized for the finite element modeling. In the nonlinear finite element formulation, the updated Lagrangian description method based on the second Piola-Kirchhoff stress tensor and the Green strain tensor is used. The buckling and postbuckling behavior of composite laminates with a delamination are investigated for various delamination sizes, stacking sequences, and boundary conditions. It is shown that the buckling load and postbuckling behavior of composite laminates depend on the buckling model which is determined by the delamination size, stacking sequence and boundary condition. Also, results show that introduction of couplings can reduce greatly the buckling load.

• Advances in Computational Design
• /
• v.4 no.1
• /
• pp.33-41
• /
• 2019

In this study, 2-mm Al/Cu bimetallic laminates were produced using asymmetric roll bonding (RB) process. The asymmetric RB process was carried out with thickness reduction ratios of 10%, 20% and 30% and mismatch rolling speeds 1:1, 1:1.1 and 1:1.2, separately. For various experimental conditions, finite element simulation was used to model the deformation of bimetallic Al/Cu laminates. Specific attention was focused on the bonding strength and bonding quality of the interface between Al and Cu layers in the simulation and experiment. The optimization of mismatch rolling speed ratios was obtained for the improvement of the bond strength of bimetallic laminates during the asymmetric RB process. During the finite element simulation, the plastic strain of samples was found to reach the maximum value with a high quality bond for the samples produced with mismatch rolling speed 1:1.2. Moreover, the peeling surfaces of samples around the interface of laminates after the peeling test were studied to investigate the bonding quality by scanning electron microscopy.

• Structural Engineering and Mechanics
• /
• v.14 no.3
• /
• pp.245-262
• /
• 2002

The free vibration analysis of stiffened laminated composite plates has been performed using the layered (zigzag) finite element method based on the first order shear deformation theory. The layers of the laminated plate is modeled using nine-node isoparametric degenerated flat shell element. The stiffeners are modeled as three-node isoparametric beam elements based on Timoshenko beam theory. Bilinear in-plane displacement constraints are used to maintain the inter-layer continuity. A special lumping technique is used in deriving the lumped mass matrices. The natural frequencies are extracted using the subspace iteration method. Numerical results are presented for unstiffened laminated plates, stiffened isotropic plates, stiffened symmetric angle-ply laminates, stiffened skew-symmetric angle-ply laminates and stiffened skew-symmetric cross-ply laminates. The effects of fiber orientations (ply angles), number of layers, stiffener depths and degrees of orthotropy are examined.

• Steel and Composite Structures
• /
• v.20 no.3
• /
• pp.693-718
• /
• 2016

This paper presents a finite element procedure based on Bridging multi-scale method (BMM) in order to incorporate the effect of local/cross-sectional deformations (e.g., flange local buckling and web crippling) on the global behaviour of thin-walled members made of fibre-reinforced polymer composite laminates. This method allows the application of local shell elements in critical regions of an existing beam-type model. Therefore, it obviates the need for using computationally expensive shell elements in the whole domain of the structure, which is otherwise necessary to capture the effect of the localized behaviour. Consequently, highly accurate analysis results can be achieved with this method by using significantly smaller finite element model, compared to the existing methods. The proposed method can be used for composite polymer laminates with arbitrary fibre orientation directions in different layers of the material, and under various loading conditions. Comparison with full shell-type finite element analysis results are made in order to illustrate the efficiency and accuracy of the proposed technique.

• Steel and Composite Structures
• /
• v.28 no.1
• /
• pp.111-121
• /
• 2018

It has been argued that fracture energy of composite laminates depends on their thickness and number of layers. In this paper a modified direct energy balance approach (DEBA) has been developed to evaluate the mode-II shear fracture energy for E-glass/Epoxy laminates from finite element model at an arbitrary thickness. This approach considers friction and damage/plasticity deformations using cohesive zone modeling (CZM) and nonlinear finite element modeling. The presence of compressive stress and resulting friction was argued to be a possible cause for the thickness dependency of fracture energy. In the finite element modeling, CZM formulation has been developed with bilinear cohesive constitutive law combined with friction consideration. Also ply element have been developed with shear plastic damage model. Modified direct energy balance approach has been proposed for estimation of mode-II shear fracture energy. Experiments were performed on laminates of glass epoxy specimens for characterization of material parameters and determination of mode-II fracture energies for different thicknesses. Effect of laminate thickness on fracture energy of transverse crack tension (TCT) and end notched flexure (ENF) specimens has been numerically studied and comparison with experimental results has been made. It is shown that the developed numerical approach is capable of estimating increase in fracture energy due to size effect.

• Structural Engineering and Mechanics
• /
• v.25 no.6
• /
• pp.733-751
• /
• 2007

The metis element method (Hung 1978) has been applied to analyse free edge interlaminar stresses and delamination in composite laminates, which are subjected to extension and bending. The paper recalls Lekhnitskii's solution for generalized plane strain state of composite laminate and Wang's singular solution for determination of stress singularity order and of eigen coefficients $C_m$ for delamination problem. Then the formulae of metis displacement finite element in two-dimensional problem are established. Computation of the stress intensity factors and the energy release rates are presented in details. The energy release rate, G, is computed by Irwin's virtual crack technique using metis elements. Finally, results of interlaminar stresses, the three stress intensity factors and the energy release rates for delamination crack in composite laminates under extension and bending are illustrated and compared with the literature to demonstrate the efficiency of the present method.

• Structural Engineering and Mechanics
• /
• v.17 no.5
• /
• pp.727-734
• /
• 2004

An uncoupled computational model for analyzing the hygrothermal dynamic response of composite laminates has been developed. The constitutive equations, expressed in an integral form, and involving relaxation moduli are adopted, to describe the non-aging hygrothermorheologically simple materials. A Prony series represents the relaxation moduli is exploited in order to derive a recursive relationship, and thereby eliminate the storage problem that arises when dealing with material possessing memory. The problem is formulated in a descritized variational form. Mindlin and higher order finite elements are employed for spatial descretization, while the Newmark average acceleration scheme is exploited for temporal descritization. The adopted recursive formula uses only the details of the previous event to compute the details of the current one. Numerical results of the displacement fields of both thin and thick viscoelastic laminates problems are discussed to show up the effectiveness of Mindlin and higher-order shear theories.

• Transactions of the Korean Society of Machine Tool Engineers
• /
• v.11 no.2
• /
• pp.59-64
• /
• 2002

This paper presets experimental and analytical studies of ultimate strength of [$[30_2/-30_2/90]_S$ carbon/epoxy laminates with free-edge delamination under uniaxial tension. We performed tensile teat far laminates with Telflon inserted on interfaces to simulate initial free-edge delamination, The experiment reveals that extensional stiffness of the laminate decreases by the initiation of the delamination, and that strength of the laminate without delamination is smaller than that of the laminates with delamination. Generalized quasi-three delamination finite element analysis, which employs energy release rate and maximum stress criteria, predicts the ultimate strength of the laminates with sufficient accuracy.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.39 no.10
• /
• pp.963-970
• /
• 2015

Laminate structures are used in the automobile, aerospace, and display industries. The advantages of fiber metal laminates are well known. Fiber metal laminates are useful for reducing the weight and improving impact resistance . However, currently, the mechanical properties of fiber metal laminates are not derived. In this paper, we use thickness as a factor for comparing the properties of laminates of various thickness combinations. The properties fiber metal laminates are analyzed using design of experiments. In addition, the finite element method is used to analyze elastic and plastic strains of fiber metal laminates and aluminum plates. The final goal of this paper is to find a suitable finite element model of fiber metal laminates under bending.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.12
• /
• pp.2982-2994
• /
• 1993

Low-velocity impact responses of composite laminates are investigated using the finite element method based on various theories. In two-dimensional nonlinear analysis, a displacement field considering higher order shear deformation and large deflection of the laminate is assumed and a finite element formulation is developed using a C$^{o}$-continuous 9-node plate element. Also, three-dimensional linear analysis based on the infinitesimal strain-displacement assumptions is performed using 8-node brick elements with incompatible modes. A modified Hertzian contact law is incorporated into the finite element program to evaluate the impact force. In the time integration, the Newmark constant acceleration algorithm is used in conjuction with successive iterations within each time step. Numerical results from static analysis as well as the impact response analysis are presented including impact force histories, deflections, strains in the laminate. Impact responses according to two typical low-velocity impact conditions are compared each other.