Application of the quadrilateral area coordinate method:a new element for laminated composite plate bending problems

Recently, some new quadrilateral finite elements were successfully developed by the Quadrilateral Area Coordinate （QAC） method. Compared with those traditional models using isoparametric coordinates, these new models are less sensitive to mesh distortion. In this paper, a new displacement-based, 4-node 20-DOF （5-DOF per node） quadrilateral bending element based on the first-order shear deformation theory for analysis of arbitrary laminated composite plates is presented. Its bending part is based on the element AC-MQ4, a recent-developed high-performance Mindlin-Reissner plate element formulated by QAC method and the generalized conforming condition method; and its in-plane displacement fields are interpolated by bilinear shape functions in isoparametric coordinates. Furthermore, the hybrid post-rocessing procedure, which was firstly proposed by the authors, is employed again to improve the stress solutions, especially for the transverse shear stresses. The resulting element, denoted as AC-MQ4-LC, exhibits excellent performance in all linear static and dynamic numerical examples. It demonstrates again that the QAC method, the generalized conforming condition method, and the hybrid post-processing procedure are efficient tools for developing simple, effective and reliable finite element models.

Optimization of buckling load for laminated composite plates using adaptive Kriging-improved PSO:A novel hybrid intelligent method

An effective hybrid optimization method is proposed by integrating an adaptive Kriging(A-Kriging)into an improved partial swarm optimization algorithm(IPSO)to give a so-called A-Kriging-IPSO for maximizing the buckling load of laminated composite plates(LCPs)under uniaxial and biaxial compressions.In this method,a novel iterative adaptive Kriging model,which is structured using two training sample sets as active and adaptive points,is utilized to directly predict the buckling load of the LCPs and to improve the efficiency of the optimization process.The active points are selected from the initial data set while the adaptive points are generated using the radial random-based convex samples.The cell-based smoothed discrete shear gap method(CS-DSG3)is employed to analyze the buckling behavior of the LCPs to provide the response of adaptive and input data sets.The buckling load of the LCPs is maximized by utilizing the IPSO algorithm.To demonstrate the efficiency and accuracy of the proposed methodology,the LCPs with different layers(2,3,4,and 10 layers),boundary conditions,aspect ratios and load patterns(biaxial and uniaxial loads)are investigated.The results obtained by proposed method are in good agreement with the literature results,but with less computational burden.By applying adaptive radial Kriging model,the accurate optimal resultsebased predictions of the buckling load are obtained for the studied LCPs.

Free Vibration Analysis of Symmetrically Laminated Composite Plates on Elastic Foundation and Coupled with Stationary Fluid

Free vibration analysis of symmetrically laminated composite plates resting on Pasternak elastic support and coupled with an ideal, incompressible and inviscid fluid is the objective of the present work. The fluid domain is considered to be infinite in the length direction but bounded in the depth and width directions. In order to derive the eigenvalue equation, Rayleigh-Ritz method is applied for the fluid-plate-foundation system. The efficiency of the method is proved by comparison studies with those reported in the open literature. At the end, parametric studies are carried out to examine the impact of different parameters on the natural frequencies.

Time-varying nonlinear dynamics of a deploying piezoelectric laminated composite plate under aerodynamic force

Using Reddy’s high-order shear theory for laminated plates and Hamilton’s principle, a nonlinear partial differential equation for the dynamics of a deploying cantilevered piezoelectric laminated composite plate, under the combined action of aerodynamic load and piezoelectric excitation, is introduced. Two-degree of freedom(DOF)nonlinear dynamic models for the time-varying coefficients describing the transverse vibration of the deploying laminate under the combined actions of a first-order aerodynamic force and piezoelectric excitation were obtained by selecting a suitable time-dependent modal function satisfying the displacement boundary conditions and applying second-order discretization using the Galerkin method. Using a numerical method, the time history curves of the deploying laminate were obtained, and its nonlinear dynamic characteristics,including extension speed and different piezoelectric excitations, were studied. The results suggest that the piezoelectric excitation has a clear effect on the change of the nonlinear dynamic characteristics of such piezoelectric laminated composite plates. The nonlinear vibration of the deploying cantilevered laminate can be effectively suppressed by choosing a suitable voltage and polarity.

The Hamilton System and Hamilton Type Generalized Variational Principle for the Laminated Composite Plates and Shells

By introducing the Hamilton theory and algorithms into the problems of laminated composite plates andshells, the Hamiltion type generalized variational principle is established, and the Hamilton canonical equations andthe boundary conditions for the static and elastoplastic analysis of composite plates are presented. With thetransformation of phase variables, the Hamilton canonical equations and their boundary conditions for thecylindrical shells and doubly curved shells in the curvilinear coordinate are given.

THE MIXED STATE HAMILTONIAN DYNAMIC ELEMENT AND A SEMI-ANALYTICAL SOLUTION FOR THE ANALYSIS OF THICK LAMINATED COMPOSITE PLATES

Through introducing the Laplace transformation in the time direction, the mixed state Hamilton canonical equation and a semi-analytical solution are presented for analyzing the dynamic response of laminated composite plates. This method accounts for the separation of variables, the finite element discretization can be employed in the plane of laminar, and the exact solution in the thickness direction is derived by the state space control method. To apply the transfer matrix method, the relational expression at the top and bottom surface is established. So the general solution in transformation space is deduced by the spot method. By the application of inversion of Laplace transformation, the transient displacements and stresses can be derived.

Thermal buckling behavior of laminated composite plates： a finite-element study

In this paper, the thermal buckling behavior of composite laminated plates under a uniform temperature distribution is studied. A finite element of four nodes and 32 degrees of freedom （DOF）, the bending and mechanical previously developed for buckling of laminated composite plates, is extended to investigate the thermal buckling behavior of laminated composite plates. Based upon the classical plate theory, the present finite element is a combination of a linear isoparametric membrane element and a high precision rectangular Hermitian element. The numerical implementation of the present finite element allowed the comparison of the numerical obtained results with results obtained from the literature： 1） with element of the same order, 2） the first order shear defo^ation theory, 3） the high order shear deformation theory and 4） the three- dimensional solution. It was found that the obtained results were very close to the reference results and the proposed element offers a good convergence speed. Furthermore, a parametrical study was also conducted to investigate the effect of the anisotropy of composite materials on the critical buckling temperature of laminated plates. The study showed that： 1） the critical buckling temperature generally decreases with the increasing of the modulus ratio EL/ET and thermal expansion ratio aT/aL, and 2） the boundary conditions and the orientation angles signifi- cantly affect the critical buckling temperature of laminated plates.

VIBRATION SUPPRESSION OF CANTILEVER LAMINATED COMPOSITE PLATE WITH NONLINEAR GIANT MAGNETOSTRICTIVE MATERIAL LAYERS

In this paper, a nonlinear and coupled constitutive model for giant magnetostrictive materials （GMM） is employed to predict the active vibration suppression process of cantilever laminated composite plate with GMM layers. The nonlinear and coupled constitutive model has great advantages in demonstrating the inherent and complicated nonlinearities of GMM in re- sponse to applied magnetic field under variable bias conditions （pre-stress and bias magnetic field）. The Hamilton principle is used to derive the nonlinear and coupled governing differential equation for a cantilever laminated composite plate with GMM layers. The derived equation is handled by the finite element method （FEM） in space domain, and solved with Newmark method and an iteration process in time domain. The numerical simulation results indicate that the proposed active control system by embedding GMM layers in cantilever laminated composite plate can efficiently suppress vibrations under variable bias conditions. The effects of embedded placement of GMM layers and control gain on vibration suppression are discussed respectively in detail.

Random Vibration Control of Laminated Composite Plates with Piezoelectric Fiber Reinforced Composites

This paper presents an analysis of the active control of random vibration for lami- nated composite plates using piezoelectric fiber reinforced composites （PFRC）. With Hamilton＇s principle and the Rayleigh-Ritz method, the equation of motion for the resulting electromechani- cal coupling system is derived. A velocity feedback control rule is employed to obtain an effective active damping in the suppression of random vibration. The power spectral density and mean- square displacements of the random vibration for laminated plates under different control gains are simulated and the validity of the present control strategy is confirmed. The effect of piezoelec- tric fiber orientation in the PFRC layer on the random vibration suppression is also investigated. The analytical methodology can be expanded to other kinds of random vibration.

Double Hopf bifurcation of composite laminated piezoelectric plate subjected to external and internal excitations

The double Hopf bifurcation of a composite laminated piezoelectric plate with combined external and internal excitations is studied. Using a multiple scale method, the average equations are obtained in two coordinates. The bifurcation response equations of the composite laminated piezoelectric plate with the primary parameter resonance, i.e., 1：3 internal resonance, are achieved. Then, the bifurcation feature of bifurcation equations is considered using the singularity theory. A bifurcation diagram is obtained on the parameter plane. Different steady state solutions of the average equations are analyzed. By numerical simulation, periodic vibration and quasi-periodic vibration responses of the Composite laminated piezoelectric plate are obtained.

Vibration characteristics and flutter analysis of a composite laminated plate with a store

The effects of an external store on the flutter characteristics of a composite laminated plate in a supersonic flow are investigated. The Dirac function is used to formulate the interaction between the plate and the store. The first-order piston theory is used to describe the aerodynamic load. The governing equation of the composite laminated plate with an external store is established based on the Hamilton principle. The mode shapes are constructed by the admissible functions which axe a set of characteristic orthogonal polynomials generated directly by the Gram-Schmidt process, and the boundary constraint is modeled as the artificial springs. The frequency and mode shapes of the plate under different boundaries are determined by the Rayleigh-Ritz method. The validity of the proposed approach is confirmed by comparing the results with those obtained from the finite element method （FEM）. The effects of the mounting position, the center of gravity position and the mounting points spacing of the external store on the flutter boundary are discussed for both the simply supported and cantilever plates, respectively, which correspond to the two installation sites of the external store, i.e., the belly and wings of the aircraft.

NUMERICAL ANALYSIS OF DELAMINATION GROWTH FOR STIFFENED COMPOSITE LAMINATED PLATES

A study of postbuckling and delamination propagation behavior in delaminated stiffened composite plates was presented. A methodology was proposed for simulating the multi-failure responses, such as initial and postbuckling, delamination onset and propagation, etc. A finite element analysis was conducted on the basis of the Mindlin first order shear effect theory and the von-Krmn nonlinear deformation assumption. The total energy release rate used as the criteria of delamination growth was estimated with virtual crack closure technique (VCCT). A self-adaptive grid moving technology was adopted to model the delamination growth process. Moreover, the contact effect along delamination front was also considered during the numerical simulation process. By some numerical examples, the influence of distribution and location of stiffener, configuration and size of the delamination, boundary condition and contact effect upon the delamination growth behavior of the stiffened composite plates were investigated. The method and numerical conclusion provided should be of great value to engineers dealing with composite structures.

ANALYSIS OF INTER-LAMINAR STRESSES FOR COMPOSITE LAMINATED PLATE WITH INTERFACIAL DAMAGE

A constitutive model for composite laminated plates with the damage effect of the intra-layers and inter-laminar interface is presented. The model is based on the general six-degrees-of-freedom plate theory, the discontinuity of displacement on the interfaces are depicted by three shape functions, which are formulated according to solutions satisfying three equilibrium equations, By using the variation principle, the three-dimensional non-linear equilibrium differential equations of the laminated plates with two different damage models are derived. Then, considering a simply supported laminated plate with damage, an analytical solution is presented using finite difference method to obtain the inter-laminar stresses.

ITERATIVE METHOD FOR LARGE DEFLECTION NONLINEAR PROBLEM OF LAMINATED COMPOSITE SHALLOW SHELLS AND PLATES

Based on the results by Wang,in this paper, the iterative method is presented for the study of large deflection nonlinear problem of laminated composite shallow shells and plates. The rectangular laminated composite shallow shells have been analyzed. The results have been compared with the small deflection linear analytical solution and finite element nonlinear solution. The results proved that the solution coincide with small deflection linear analytical solution in the condition of the low loads and finite element nonlinear solution in the condition of the high loads.

Power transmission through double-walled laminated composite panels considering porous layer-air gap insulation

The acoustic behavior of double-walled laminated composite panels consisting of two porous and air gap middle layers is studied within the classical laminated plate theory （CLPT）. Thus, viscous and inertia coupling in a dynamic equation, as well as stress transfer, thermal and elastic coupling of porous material ave based on the Biot theory. In addition, the wave equations are extracted according to the vibration equation of composite layers. The transmission loss （TL） of the structure is then calculated by solving these equations simultaneously. Statistical energy analysis （SEA） is developed to divide the structure into specific subsystems, and power transmission is extracted with balancing power flow equations of the subsystems. Comparison between the present work and the results reported elsewhere shows excellent agreement. The results also indicate that, although favorable enhancement is seen in noise control particularly at high frequencies, the corresponding parameters associated with fluid phase and solid phase of the porous layer are important on TL according to the boundary condition interfaces. Finally, the influence of composite material and stacking sequence on power transmission is discussed.

KRMN-TYPE EQUATIONS FOR A HIGHER-ORDER SHEARDEFORMATION PLATE THEORY AND ITS USE IN THE THERMAL POSTBUCKLING ANALYSIS

arman-type nonlinear large deflection equations are derived occnrding to theReddy's higher-order shear deformation plate theory and used in the thermal postbuckling analysis The effects of initial geometric imperfections of the plate areincluded in the present study which also includes th thermal effects.Simply supported,symmetric cross-ply laminated plates subjected to uniform or nomuniform parabolictemperature distribution are considered. The analysis uses a mixed GalerkinGolerkinperlurbation technique to determine thermal buckling louds and postbucklingequilibrium paths.The effects played by transverse shear deformation plate aspeclraio, total number of plies thermal load ratio and initial geometric imperfections arealso studied.

A size-dependent composite laminated skew plate model based on a new modified couple stress theory

In this study, a size-dependent composite laminated skew Mindlin plate model is proposed based on a new modified couple stress theory. This plate model can be viewed as a simplified couple stress theory in engineering mechanics. Governing equations and related boundary conditions are derived based on the principle of minimum potential energy. The Rayleigh–Ritz method is employed to obtain the numerical solutions of the center deflections of simply supported plates with different ply orientations. Numerical results show that the normalized center deflections obtained by the proposed model are always smaller than those obtained by the classical one, i.e. the present model can capture the scale effects of microstructures. Moreover, a phenomenon reveals that the ply orientation would make a significant influence on the magnitude of scale effects of composite laminated plates at micro scale. Additionally, the present model of thick skew plate can be degenerated to the model of Kirchhoff plate based on the modified couple stress theory by adopting the assumptions in Bernoulli–Euler beam and material isotropy.

MULTI-PULSE ORBITS AND CHAOTIC DYNAMICS OF A COMPOSITE LAMINATED RECTANGULAR PLATE

This paper presents an analysis on the nonlinear dynamics and multi-pulse chaotic motions of a simply-supported symmetric cross-ply composite laminated rectangular thin plate with the parametric and forcing excitations. Firstly, based on the Reddy＇s third-order shear deformation plate theory and the model of the yon Karman type geometric nonlinearity, the nonlinear governing partial differential equations of motion for the composite laminated rectangular thin plate are derived by using the Hamilton＇s principle. Then, using the second-order Galerkin dis- cretization, the partial differential governing equations of motion are transformed to nonlinear ordinary differential equations. The case of the primary parametric resonance and 1：1 internal resonance is considered. Four-dimensional averaged equation is obtained by using the method of multiple scales. From the averaged equation obtained here, the theory of normal form is used to give the explicit expressions of normal form. Based on normal form, the energy phase method is utilized to analyze the global bifurcations and multi-pulse chaotic dynamics of the composite laminated rectangular thin plate. The theoretic results obtained above illustrate the existence of the chaos for the Smale horseshoe sense in a parametrical and forcing excited composite laminated thin plate. The chaotic motions of the composite laminated rectangular thin plate are also found by using numerical simulation, which also indicate that there exist different shapes of the multi-pulse chaotic motions for the composite laminated rectangular thin plate.

Nonlinear Vibrations of Two-Span Composite Laminated Plates with Equal and Unequal Subspan Lengths

The nonlinear transverse vibrations of ordered and disordered twodimensional(2D)two-span composite laminated plates are studied.Based on the von Karman’s large deformation theory,the equations of motion of each-span composite laminated plate are formulated using Hamilton’s principle,and the partial differential equations are discretized into nonlinear ordinary ones through the Galerkin’s method.The primary resonance and 1/3 sub-harmonic resonance are investigated by using the method of multiple scales.The amplitude-frequency relations of the steady-state responses and their stability analyses in each kind of resonance are carried out.The effects of the disorder ratio and ply angle on the two different resonances are analyzed.From the numerical results,it can be concluded that disorder in the length of the two-span 2D composite laminated plate will cause the nonlinear vibration localization phenomenon,and with the increase of the disorder ratio,the vibration localization phenomenon will become more obvious.Moreover,the amplitude-frequency curves for both primary resonance and 1/3 sub-harmonic resonance obtained by the present analytical method are compared with those by the numerical integration,and satisfactory precision can be obtained for engineering applications and the results certify the correctness of the present approximately analytical solutions.

Finite element analysis of stress concentrations and failure criteria in composite plates with circular holes

In this study, the （SCF） in cross-and-angle-ply stress concentration factors laminated composite plates as well as in isotropic plates with single circular holes subjected to uniaxial loading is studied. A quadrilateral finite element of four-node with 32 degrees of freedom at each node, previously developed for the bending and mechanical buckling of laminated composite plates, is used to evaluate the stress distribution in laminated composite plates with central circular holes. Based up on the classical plate theory, the present finite element is a combination of a linear isoparametric membrane element and a high precision rectangular Hermitian element. The numerical results obtained by the present element compare favorably with those obtained by the analytic approaches published in literature. It is observed that the obtained results are very close to the reference results, which demonstrates the accuracy of the present element. Additionally, to determine the first ply failure （FPF） of laminated plate, several failure criterions are employed. Finally, to show the effect of E1/E2 ratio on the failure of plates, a number of figures are given for different fiber orientation angles.