Buckling of shear deformable laminated composite plates

A shear deformable finite element is developed for the buckling analysis of laminated composite plates. The finite element formulation is based on Mindlin's theory in which shear correction factors are derived from the exact expressions for orthotropic materials. A variety of problems on uniaxial and shear bucklings of laminated composite plates are solved. The effects of material properties, plate aspect ratio, length-to-thickness ratio, number of layers and lamination angle on the buckling loads of symmetrically and antisymmetrically laminated composite plates are investigated. Optimal lamination arrangements of layers for maximizing the buckling loads of the plates are determined.     

A two-step optimization scheme for maximum stiffness design of laminated plates based on lamination parameters

The purpose of this paper is to propose an effective solution scheme of simultaneous optimization design of layup configuration and fiber distribution for maximum stiffness design of laminated plates. Firstly, a numerical analysis of the lamination parameters feasible region for a laminated plate consisting of various given number of ply groups (each ply group may have different thickness and all the fibers in one ply group are orientated in an identical direction) is carried out, and it is found that the feasible region based on only a few ply groups is very close to the overall one determined by infinite plies. Therefore, it is suggested that the feasible region of lamination parameters of a laminated plate could be approximately determined by the layup configuration of least ply groups. Secondly, a two-step simultaneous optimization scheme of layup configuration and fiber distribution for maximum stiffness design of laminated plates is proposed. Accordingly, by using ply thickness, fiber orientation angle and fiber volume fraction in a laminated plate of least ply groups as design variables, the optimal lamination parameters for maximum stiffness is obtained. Then, taking the optimal lamination parameters as the design objective, a detailed layup design optimization is implemented by considering some limitations on manufacturing, such as preset ply thickness, and specific fiber orientation angle and a limited maximum number of consecutive plies in the same fiber orientation. Numerical examples are also presented to validate the proposed two-step optimization scheme.     

Three-dimensional analysis for rectangular 1–3 piezoelectric fiber-reinforced composite laminates with the interdigitated: electrodes under electromechanical loadings

In this paper, an exact analysis for the rectangular composite laminated plate consisting of 1–3 piezoelectric fiber-reinforced composite layer and orthotropic elastic composite layer subjected to the electric field with the interdigitated electrode (IDE) and arbitrary loads is presented without any simplification. The solution of the derived governing differential equations is obtained through the power series expansion and Fourier expansion methods. An illustrative analysis is carried out to investigate the influence of the stiffness anisotropy, free-strain anisotropy and electric field on the shear curvature of the laminated plate with 1–3 active fiber-reinforced composite layer. The numerical results show that, the magnitude of shear curvature gradually increases as the stiffness anisotropy and the free-strain anisotropy increase, and increasing the electrical field also leads to an increase in shear curvature. Results presented here can be used to assess various approximate theories and enhance the understanding of the static and dynamic response behavior of the 1–3 piezoelectric composite structures.     

The vibration damping of laminated plates

A method previously developed for determination of elastic and damping parameters of orthotropic plates (McIntyre, M. E. and Woodhouse, J., Acta Metall., 1988, 36, 1397–1416) was applied to laminated composite plates. The necessary theory is summarised, and the predictions of laminate theory compared with experimental results for three CFRP laminated plates with different constructions. It is also shown that laminate theory can be inverted, to obtain the ply properties from measurements on the laminated plate. This can sometimes afford a good way to obtain the necessary calibration data on the material properties of the plies.     

Off-axis transformation of the composite laminate stiffness properties

Most composite structures are orthotropic with respect to the major structural loading direction, i.e. the 0 ° fibres are along the principal bending axes of an aircraft wing (the spar line). The laminate stiffness properties are given with the laminate orthotropic axes aligned with the structural axes. If the orthotropic axes are not aligned to the structural axes, then the laminate stiffness properties are generated from classical laminated plate theory with individual ply angles rotated through the appropriae angle of transformation. In this discussion, the transformation of a laminate's stiffness from the on-axis position to an off-axis position is accomplished in one step. Two typical examples are shown to illustrate where and why such a transformation is used.     

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.     

Adaptive vibration control of a cylindrical shell with laminated PVDF actuator

A generalized higher-order theory describing the mechanical behavior of multi-layered composite plates with arbitrary lamination scheme is proposed. Ritz’s method is employed to determine the kinematic unknowns expressed in a complete polynomial power series of the thickness-wise coordinate whereas the dependence on the in-plane coordinates is such that the functions satisfy all boundary conditions. The correct constitutive laws of a three-dimensional orthotropic elastic continuum are employed for each individual layer. The convergence and accuracy of the computational scheme are investigated by comparing elastic static and buckling results with analytical or finite element solutions for complex cross- and angle-ply laminates. For further validation of the theory, laminated plates under a transverse pressure are investigated for technically relevant lamination schemes and the associated deformation and stress results are compared with those obtained through FE calculations.     

Elastic stability of skew laminated composite plates subjected to biaxial follower forces

The present study investigates the elastic stability of skew laminated composite plates subjected to biaxial inplane follower forces by the finite element method. The plate is assumed to follow first-order shear deformation plate theory (FSDPT). The kinetic and strain energies of skew laminated composite plate and the work done by the biaxial inplane follower forces are derived by using tensor theory. Then, by Hamilton's principle, the dynamic mathematical model to describe the free vibration of this problem is formed. The finite element method and the isoparametric element are utilized to discretize the continuous system and to obtain the characteristic equations of the present problem. Finally, natural vibration frequencies, buckling loads (also the instability types) and their corresponding mode shapes are found by solving the characteristic equations. Numerical results are presented to demonstrate the effects of those parameters, such as various inplane force combinations, skew angle and lamination scheme, on the elastic stability of skew laminated composite plates subjected to biaxial inplane follower forces.     

Optimal design of laminated composite plates using a global optimization technique

A multi-start global optimization technique is used to investigate the lamination arrangements of laminated composite plates designed for maximum stiffness. The multi-start global optimization technique which originated from the concept of minimizing the potential energy of a moving particle in a conservative force field is extended to the optimal design of laminated composite plates in which the strain energies of the plates are minimized. The optimization algorithm has been proved to be efficient and effective in producing the global optima. Numerical examples of the selection of optimal lamination arrangements of symmetrically laminated composite plates with different aspect ratios subject to different loading conditions are given. The results show that aspect ratio, loading condition and material property can affect the optimal lamination arrangement.     

A generalized higher-order theory for buckling of thick multi-layered composite plates with normal and transverse shear strains

The onset of buckling in square laminated multi-layered composite plates, subject to unidirectional in-plane loads, is investigated within the framework of a generalized higher-order shear deformation theory suitable to capture significant transverse shear and thickness-wise deformation effects. The displacement field is expanded in a Taylor series of the thickness coordinate with arbitrary polynomial degree; in turn, the series coefficients, expressed as a superposition of admissible functions, are determined according to the Rayleigh–Ritz method. Truly higher-order polynomial terms, along with a sufficient number of in-plane admissible functions, are shown to be necessary for convergence towards the fundamental buckling load multiplier. As a by-product, reduced-order models are identified for various plate geometries and lamination schemes. The sensitivity of the lowest buckling load with respect to the nondimensional parameters (the thickness ratio, the ratio between the elastic moduli, the ply angle) is investigated. In particular, the attention is focused on the cross-over phenomenon between the lowest two buckling eigenvalues in multi-layered composite square plates with different lamination schemes. The presented results shed light onto the buckling behavior of thick shear-deformable multi-layered plates.     

Reliability analysis of nonlinear laminated composite plate structures

A procedure for the reliability analysis of laminated composite plate structures subjected to large deflections under random static loads is presented. The nonlinear analysis of laminated composite plate structures is achieved via a corotational total Lagrangian finite element formulation which is based on the von Karman assumption and first order shear deformation theory. This formulation is applicable for the nonlinear analysis of plate structures with large rotations but moderate deformation and thus accurate enough to predict the behavior of the structures at the point of failure. The reliability assessment of laminated composite plate structures with random strength subjected to random loads is approached by the determination of limit state surfaces in load space. The limit space surfaces are obtained by performing a series of first ply failure analyses following different load paths in load space using the proposed nonlinear structural analysis technique and an appropriate failure criterion. A numerical technique is then proposed to evaluate the reliability of the plate structures. Examples of the reliability analyses of laminated plates with different layer orientations subject to random loads are given for illustration.     

Free vibration analysis of composite sandwich plates based on Reddy's higher-order theory

Two new C⁰ assumed strain finite element formulations of Reddy's higher-order theory are used to determine the natural frequencies of isotropic, orthotropic, and layered anisotropic composite and sandwich plates. The material properties typical of glass fibre polyester resins for the skin and HEREX C70 PVC (polyvinyl chloride) foam materials for the core are used to show the parametric effects of plate aspect ratio, length-to-thickness ratio, degree of orthotropy, number of layers and lamination scheme on the natural frequencies. A consistent mass matrix is adopted in the present formulation. The results presented in this investigation could be useful for a better understanding of the behaviour of sandwich laminates under free vibration conditions and potentially beneficial for designers of sandwich structures.     

Transient response of composite sandwich plates

The transient response of orthotropic, layered composite sandwich plates is investigated by using two new C⁰ four and nine node finite element formulations of a refined form of Reddy's higher-order theory. This refined third order theory accounts for parabolic variation of the transverse shear stresses, and requires no shear correction factors. The assumed strain approach is employed to model both thin and thick plates without any major defects like shear locking and parasitic spurious zero energy modes. A consistent mass matrix formulation is adopted. The Newmark direct integration scheme is used to solve the governing equilibrium equations. The parametric effects of plate aspect ratio, length to thickness ratio, boundary conditions and lamination scheme on the transient response are investigated. The present results are in very close agreement with earlier published results in the literature and can serve as a benchmark for future investigators.     

The effects of geometric parameters on the failure strength for pin-loaded multi-directional fiber-glass reinforced epoxy laminate

A numerical and experimental study was carried out to determine the failure of mechanically fastened fiber-reinforced laminated composite joints. E/glass–epoxy composites were manufactured to fabricate the specimens. Mechanical properties and strengths of the composite were obtained experimentally. Tests have been carried out on single pinned joints in [0/90/0]_s and [90/0/90]_s laminated composites. A parametric study considering geometries was performed to identify the failure characteristics of the pin-loaded laminated composite. Data obtained from pin-loaded laminate tests were compared with the ones calculated from a finite element model (PDNLPIN computer code). Damage accumulations in the laminates were evaluated by using Hashin's failure criteria combined with the proposed property degradation model. Based on the results, ply orientation and geometries of composites could be crucial for pinned laminated composite joints.     

An engineering formula for the stress concentration factor of orthotropic composite plates

An engineering formula for the theoretical stress concentration factor of orthotropic notched plates under tension is provided, as a function of the material elastic constants and the K_t of the corresponding isotropic case. The accuracy and limits of applicability of the new solution are discussed by comparison to data from the literature and results from FE analyses on notched geometries of practical interests. The proposed solution represents a very useful tool to estimate the stress concentration factor of notched orthotropic plates, composite orthotropic laminae, orthotropic unidirectional laminates and homogenised orthotropic composite laminates.     

Optimal lay-up design of variable stiffness laminated composite plates by a layer-wise optimization technique

The optimal lay-up design for the maximum fundamental frequency of variable stiffness laminated composite plates is investigated using a layer-wise optimization technique. The design variables are two fibre orientation angles per ply. Thin plate theory is used in conjunction with a p-element to calculate the fundamental frequencies of symmetrically and antisymmetrically laminated composite plates. Comparisons with existing optimal solutions for constant stiffness symmetrically laminated composite plates show excellent agreement. It is observed that the maximum fundamental frequency can be increased considerably using variable stiffness design as compared to constant stiffness design. In addition, optimal lay-ups for the maximum fundamental frequency of variable stiffness symmetrically and antisymmetrically laminated composite plates with different aspect ratios and various combinations of free, simply supported and clamped edge conditions are presented. These should prove a useful benchmark for optimal lay-ups of variable stiffness laminated composite plates.     

Transient analysis of FGM and laminated composite structures using a refined 8-node ANS shell element

In this paper, we investigate the vibration analysis of functionally graded material (FGM) and laminated composite structures, using a refined 8-node shell element that allows for the effects of transverse shear deformation and rotary inertia. The properties of FGM vary continuously through the thickness direction according to the volume fraction of constituents defined by sigmoid function, but in this method, their Poisson’s ratios of the FGM plates and shells are assumed to be constant. The finite element, based on a first-order shear deformation theory, is further improved by the combined use of assumed natural strains and different sets of collocation points for interpolation the different strain components. We analyze the influence of the shell element with the various location and number of enhanced membrane and shear interpolation. Using the assumed natural strain method with proper interpolation functions the present shell element generates neither membrane nor shear locking behavior even when full integration is used in the formulation. The natural frequencies of plates and shells are presented, and the forced vibration analysis of FGM and laminated composite plates and shells subjected to arbitrary loading is carried out. In order to overcome membrane and shear locking phenomena, the assumed natural strain method is used. To validate and compare the finite element numerical solutions, the reference solutions of plates based on the Navier’s method, the series solutions of sigmoid FGM (S-FGM) plates are obtained. Results of the present theory show good agreement with the reference solutions. In addition the effect of damping is investigated on the forced vibration analysis of FGM plates and shells.     

Nonlinear flexural analysis of laminated composite plates using RBF based meshless method

The paper presents the nonlinear flexural response of laminated composite plates. The mathematical formulation of the actual physical problem of the laminated composite plate subjected to mechanical loading is presented utilizing higher order shear deformation theory and von-Karman nonlinear kinematics. These nonlinear governing differential equations of equilibrium are linearized using quadratic extrapolation technique. A meshfree technique based on multiquadric RBFs is used for analysis of the problems. Isotropic, orthotropic and laminated composite plates with immovable simply supported and clamped edges are analyzed.     

Thermal buckling and postbuckling behavior of functionally graded carbon nanotube-reinforced composite plates

Thermal buckling and postbuckling behavior is presented for functionally graded nanocomposite plates reinforced by single-walled carbon nanotubes (SWCNTs) subjected to in-plane temperature variation. The material properties of SWCNTs are assumed to be temperature-dependent and are obtained from molecular dynamics simulations. The material properties of functionally graded carbon nanotube-reinforced composites (FG-CNTRCs) are assumed to be graded in the thickness direction, and are estimated through a micromechanical model. Based on the multi-scale approach, numerical illustrations are carried out for perfect and imperfect, geometrically mid-plane symmetric FG-CNTRC plates and uniformly distributed CNTRC plates under different values of the nanotube volume fractions. The results show that the buckling temperature as well as thermal postbuckling strength of the plate can be increased as a result of a functionally graded reinforcement. It is found that in some cases the CNTRC plate with intermediate nanotube volume fraction does not have intermediate buckling temperature and initial thermal postbuckling strength.     

Thermal postbuckling analysis of laminated composite plates by the finite element method

The thermal postbuckling behavior of composite laminated plates subjected to a nonuniform temperature field is investigated by the finite element method. Based on the principle of minimum potential energy, the nonlinear stiffness matrix and geometry matrix are derived. The assumed displacement state over the middle surface of the plate element is expressed as a product of one-dimensional, first-order Hermitian polynomials. An iterative method is employed to determine the thermal postbuckling load. The results of the computations reveal that the thermal postbuckling behavior of composite laminated plates is influenced by lamination angle, plate aspect ratio, modulus ratio and the number of layers.