# Thermal Buckling of Piezoelectric Composite Beam

Authors: S. Yazdani; Y. Kiani; M. Jabbari; M. R. Eslami
Journal: ISRN Mechanical Engineering
 (2011)

---

## Abstract

Buckling analysis of laminated composite beams with piezoelectric
layers subjected to thermal loading and constant voltage is studied. The material properties are assumed to be homogeneous in any layer through the beam thickness. The first-order beam theory and nonlinear
strain-displacement relation are used to obtain the governing equations of the composite beam. The beam is assumed under uniform type of thermal loading and various types of boundary conditions. For each case of boundary conditions, closed-form solutions are obtained. The effects of the applied actuator voltage, beam geometry, and boundary conditions on the buckling temperature are investigated.


---

## Body

## 1. Introduction

Static and dynamic analysis for multilayer composite structures have been well established for various engineering applications during the last decades. Brush and Almroth [1] have a general treatment on the subject of structural stability, including beams, plates, and shells. Wang et al. [2] presented the closed-form solutions for buckling of beams, plates, and shells based on the classical, first-order, and higher-order displacement theories under compressive loads. Eslami and Shariyat [3, 4] used the improved equations to obtain the elastic, plastic, and creep buckling of thin cylindrical shells under different mechanical loading conditions. Analytical solutions of refined beam theories are developed to study the buckling behavior of cross-ply rectangular beams with arbitrary boundary conditions [5]. Kolakowski et al. [6] presented a modal interactive buckling of thin-walled composite beam columns regarding distortional deformations. Buckling analysis of cross-ply laminated beams with general boundary conditions by Ritz method is studied by Aydogdu [7].If the membrane stresses due to a temperature distribution in a composite laminate are compressive and sufficiently large, equilibrium may become unstable, and thermal buckling may occur. In recent years, many studies have focused on the analysis of the thermal buckling and postbuckling responses of composite laminates. Eslami et al. [8] obtained the thermoelastic buckling of thin cylindrical shells under a number of practical thermal loadings. Shear deformation effects on thermal buckling of cross-ply composite laminates have been studied by Mannini [9]. In this paper, thermal buckling of symmetric and antisymmetric cross-ply composite laminates is investigated. The first-order shear deformation theory in conjunction with the Rayleigh-Ritz method is used for the evaluation of the thermal buckling parameters of structures.Jordan canonical form solution for thermally induced deformation of cross-ply laminated composite beams has been presented by Khdeir and Reddy [10]. Also, Khdeir [11] studied the thermal buckling of thick, moderately thick, and thin cross-ply laminated beams subjected to uniform temperature distribution. He presented the exact analytical solutions of refined beam theories to obtain the critical buckling temperature of cross-ply beams with various boundary conditions. Li and Song [12] studied the large thermal deflections of Timoshenko beams under transversely nonuniform temperature rise. Thermal buckling analysis of cross-ply laminated composite beams with general boundary conditions is presented by Aydogdu [13]. The study is concerned with the thermal buckling analysis of cross-ply laminated beams subjected to different sets of boundary conditions. The analysis is based on a three-degrees-of-freedom shear deformable beam theory. The governing equations are obtained by means of the minimum energy principle. Thermal buckling load optimization of angle-ply symmetrically laminated composite beams is studied by Topal [14]. The objective of the optimization problem is to maximize the critical thermal buckling load of the laminated beams, and the fibre orientation is considered as the design variable.Advanced structures with integrated self-monitoring and control capabilities are increasingly becoming important due to the rapid development of smart structure and mechanical systems. Bailey and Hubbard [15] reported vibration control of a piezoelectric beam with a simplified beam model. Recently, discrete layer theories are utilized for the analysis of composite structures with piezoelectrics in order to fully consider the effects of the transverse shear and variable in-plane displacements [16]. Tzou et al. [17, 18] proposed the mathematical modelling of nonlinear thermopiezoelastic laminates and investigated the static and dynamic control of beams and plates. Abramovich [19] presented the closed-form solutions for deflection control of laminated composite beams with piezoceramic layers. In his study, the three coupled equations of motion of a general nonsymmetric piezolaminated composite beam subjected to axial and lateral traction, and its corresponding boundary conditions are derived using a variational approach. The static shape control is performed using either continuous piezoceramic layers or patches embedded or bonded to the surface of the beam structure. Closed-form solutions for the bending angle and the axial lateral displacements along the beam are presented for various configurations of layup, boundary conditions, and mechanical loading. Waisman and Abramovich [20] studied the active stiffening of laminated composite beams using piezoelectric actuators. The present study deals with the stiffening effects of a smart piezolaminated composite beam. The structure consists of piezoceramic layers or patches bonded on the surface of the beam. The analysis considers the linear piezoelectric constitutive relations and the first-order shear deformation theory. Aldraihem and Khdeir [21, 22] presented the exact deflection solutions of beams with shear piezoelectric patches and actuators. Jerome and Ganesan [23] developed a generalized plane strain finite element formulation to predict the critical buckling voltage and temperature of a piezo composite beam. Akhras and Li [24] proposed the three-dimensional thermal buckling analysis of piezoelectric antisymmetric angle-ply laminates using finite layer method.In this paper, the thermal buckling of piezoelectric laminated composite beams is studied. The first-order shear deformation beam theory is employed, and the closed-form solutions are presented for different types of boundary conditions.


## 2. Formulation of Problem

Consider a laminated composite beam with lengthL, width c, and total thickness h. The rectangular Cartesian coordinates is used such that the x axis is along the length of the beam on its middle surface and z is measured from the middle surface and is positive upward, as shown in Figure 1. The analysis is based on the first-order beam theory. The displacement field for the beam is u¯ and w¯, which is based on Timoshenko beam theory, can be written as (1)u¯(x,z)=u(x)+zϕ(x),w¯(x,z)=w(x),
where u and w are the axial and lateral displacements of a point on the midplane and ϕ is the bending rotation of the normal to the mid plane. The normal strain ɛx and the transverse shear strain γxz at any point in the laminate are (2)ɛx=∂u¯∂x+12(∂w¯∂x)2=u′+12w′2+zϕ′,γxz=∂u¯∂z+∂w¯∂x=ϕ+w′,
where a  ′  stands for a derivation respect to x. When piezo composite beam is subjected to thermal load, the force and moment equations are written as [11, 19](3)[NxMxQxz]=[A11B110B11D11000A55][u′+12w′2ϕ′ϕ+w′]-[NTMT0]-[NEME0],
where in this equation (4)Nx=∫-h/2h/2cσxdz,Mx=∫-h/2h/2cσxzdz,Qxz=∫-h/2h/2cτxzdz,σx and τxz being the normal and shear stresses, respectively. Thermal force and thermal moment are (5)NT=c∑n=1N∫zn-1znQ11nαxnΔTdz,MT=c∑n=1N∫zn-1znQ11nαxnΔTzdz.
Here, N is the number of layers and αx is the axial coefficient of thermal expansion. Terms NEand ME are the piezoelectric force and moment and are (6)NE=c∑n=1Na(Q11)anVnd31n,ME=c2∑n=1Na(Q11)anVnd31n(2zan+han).
A subscript a stands for quantities associated with piezoelectric layers. Here, Vn is the applied actuator voltage to the surface of nth piezoelectric layer nsd d31 is the piezoelectric constant. Also, A11, B11, D11, and A55 are the usual extensional, bending-extension, bending, and transverse shear stiffness coefficients defined as (7)A11=c∫-h/2h/2Q¯11dz=c∑n=1N(Q¯11)n(zn+1-zn),B11=c∫-h/2h/2Q¯11zdz=c2∑n=1N(Q¯11)n(zn+12-zn2),D11=c∫-h/2h/2Q¯11z2dz=c3∑n=1N(Q¯11)n(zn+13-zn3),A55=ck∫-h/2h/2Q¯55zdz=ck∑n=1N(Q¯55)n(zn+1-zn),
where N is the number of layers, k is a shear correction factor, and Q¯11 and Q¯55 are the transformed material constants given by (8)Q¯11=Q11cos⁡4⁡θ+Q22sin4⁡θ+2(Q12+2Q66)sin2⁡θcos⁡2⁡θ,Q¯55=G13cos⁡2⁡θ+G23sin2θ.
The angle θ is the angle between the fibre direction and longitudinal axis (x axis) of the beam, and the constants Q11, Q12, Q22, and Q66 are (9)Q11=E111-ν12ν21,Q22=E221-ν12ν21,Q12=E11ν121-ν12ν21,Q66=G12.
Using the principle of minimum total potential energy, the governing equations for the displacement field of (1) are derived in [1, 11, 19] as(10)dNxdx=0,dMxdx-Qxz=0,dQxzdx+Nxd2wdx2=0.Figure 1
Geometry of cross-ply composite beam.In this paper, it is assumed that the thermal load is uniform. The equilibrium equations in terms of the displacement components are obtained by substituting (3) into  (10)(11)A11(u′′+w′w′′)+B11ϕ′′=0,B11(u′′+w′w′′)+D11ϕ′′-A55(ϕ+w′)=0,A55(ϕ′+w′′)+Nxw′′=0.


## 3. Prebuckling Deformation

The flat prebuckling configurations are assumed. For this purpose, the prebuckling deformation of laminated composite beam should be studied to assure that the beam remains flat under uniform thermal loading. The deformation of a beam prior to buckling may be obtained by solving the equilibrium equations (11) with the nonlinear terms set equal to zero [11](12)A11u′′+B11ϕ′′=0,B11u′′+D11ϕ′′-A55(ϕ+w′)=0,A55(ϕ′+w′′)=0.
Solving these equations, we obtain (13)u=B11A11A55((B112/A11)-D11)b1x22+b5x+b6,w=A55((B112/A11)-D11)b1x36+b2x22+b3x+b4,ϕ=-A55((B112/A11)-D11)b1x22-b2x+b3+b1,Nx=A11b5-B11b2-NT-NE,Mx=A55b1x+B11b5-D11b2-MT-ME,Qxz=A55b1,
where b1 to b6 are constants which have to be determined using the associated boundary conditions. The prebuckling boundary conditions are listed in Table 1.  For each case of boundary conditions, constants b1 to b6 have been evaluated and listed in Table 2,(14)I=A55(B112/A11)-D11,G=A55L-B112A11IL2+D111+I(L2/6)L/2.
From this table, one may obtain that except the Clamped-Clamped and Clamped-Roller laminated composite beams, the other types of boundary conditions under thermal loading initially start to deflect rather than buckling. But the C-C and C-R boundary condition follow the bifurcation type buckling for uniform temperature rise loading.Table 1
Prebuckling boundary conditions for various edge supports. (C indicates clamped, S shows simply supported and R is used for roller edge).
B.C.B.Cs atx=0B.Cs atx=LS-Su=w=Mx=0u=w=Mx=0C-Cu=w=ϕ=0u=w=ϕ=0C-Su=w=ϕ=0u=w=Mx=0C-Ru=w=ϕ=0u=ϕ=dMxdx=0S-Ru=w=Mx=0u=ϕ=dMxdx=0Table 2
Constantsbi for various edge conditions.
B.C.b1b2b3b4b5b6C-C000000C-R000000C-SMTG2L(1+IL26)MT+MEGMT+MEG0-B11A11IMT+MEGL20S-R0-MT+MED11-MT+MED11L000S-S0-MT+MED11MT+MED11L2000


## 4. Stability Equations

To derive the stability equations, the adjacent-equilibrium criterion is used. Assume that the equilibrium state of a laminated composite beam is defined in terms of the displacement componentsu0, w0, and ϕ0 and the displacement components of a neighboring stable state differ by u1, w1, and ϕ1 with respect to the equilibrium position. Thus, the total displacements of a neighboring state are [1] (15)u=u0+u1,w=w0+w1,ϕ=ϕ0+ϕ1.
Similar to the displacements, the force and moment of a neighboring state may be related to the state of equilibrium as (16)Nx=Nx0+Nx1,Mx=Mx0+Mx1,Qxz=Qxz0+Qxz1.
Here, Nx1, Mx1, and Qxz1 represent the linear parts of the force and moment increments corresponding to u1, w1, and ϕ1. The stability equations may be obtained by substituting (15) and (16) in (3). Upon substitution, the terms in the resulting equations with subscript 0 satisfy the equilibrium conditions and, therefore, drop out of the equations. The remaining terms form the stability equations as(17)dNx1dx=0,dMx1dx-Qxz1=0,dQxz1dx+Nx0d2w1dx2=0.
Using (3) and (15), the force and moment with subscript 1 may be defined by(18)Nx1=A11u1′+B11ϕ1′,Mx1=B11u1′+D11ϕ1′,Qxz1=A55(ϕ1+w1′).
For C-C and C-R composite beam subjected to uniform temperature rise, one may obtain(19)Nx0=-NT-NE,Mx0=-MT-ME.
Combining (17) and (18) by eliminating u1 and ϕ1 provides an ordinary differential equation in terms of w1, which is the stability equation of composite beam under thermal loading(20)d4w1dx4+μ2d2w1dx2=0,
with(21)μ2=-Nx0(D11-(B112/A11))(1+(Nx0/A55)).
When the temperature distribution in composite beam is uniform, the parameter μ is constant, and then the exact solution of (20) is(22)w1(x)=C1sin⁡(μx)+C2cos⁡⁡(μx)+C3x+C4.
Using (17), (18), and (22), the expressions for u1, ϕ1, and Nx1, Mx1, Qxz1 become(23)u1(x)=B11A11(1-μ2(D11-(B112/A11))A55+μ2(D11-(B112/A11)))×μ(C1cos⁡⁡(μx)-C2sin⁡(μx))+C5x+C6,(24)ϕ1(x)=(1-μ2(D11-(B112/A11))A55+μ2(D11-(B112/A11)))×μ(-C1cos⁡⁡(μx)+C2sin⁡(μx))-C3,(25)Nx1(x)=A11C5,Mx1(x)=(D11-B112A11)(1-μ2(D11-(B112/A11))A55+μ2(D11-(B112/A11)))×μ2(C1cos⁡⁡(μx)+C2sin⁡(μx))+B11C5,Qxz1(x)=(μ3(D11-(B112/A11))1+(μ2/A55)(D11-(B112/A11)))×(C1cos⁡⁡(μx)-C2sin⁡(μx)).Constants of these equations (C1 to C6) are obtained using the boundary conditions of the composite beam. To find the minimum value of Nx0 associated with the thermal buckling load, the parameter μ must be minimized. Five types of boundary conditions are assumed for the composite beam. Consider a beam with both edges clamped. The edge conditions of the clamped-clamped composite beam are(26)u1(0)=w1(0)=ϕ1(0)=u1(L)=w1(L)=ϕ1(L)=0.Using (22)–(24) and (26), the constants C1 to C6 must satisfy the system of equations(27)[010100sin⁡(μL)cos⁡⁡(μL)L100B11A11Pμ00001B11A11Pμcos⁡⁡(μL)-B11A11Pμsin⁡(μL)00L1-Pμ0-1000-Pμcos⁡⁡(μL)Pμsin⁡(μL)-1000][C1C2C3C4C5C6]=[000000],
where(28)P=(1-μ2(D11-(B112/A11))A55+μ2(D11-(B112/A11))).
To have a nontrivial solution, the determinant of coefficient matrix must be zero, which yields(29)PμL(2-2cos⁡⁡(μL)+PμLsin⁡(μL))=0.
The smallest positive value of μ which satisfies (29) is μmin=2π/L. Table 3 shows different types of boundary conditions and the minimum values of μ associated with the thermal buckling loads. Now, the critical force for buckling from (21) (except for C-S beam, where the approximate solution from [2] is considered) is(30)Nx0=-μ2(D11-(B112/A11))1+(μ2/A55)(D11-(B112/A11)).
Then with this equation and (19), the buckling force of the beam for all cases of boundary conditions can be written in the form(31)NT+NE=μ2(D11-(B112/A11))1+(μ2/A55)(D11-(B112/A11)).Table 3
Boundary conditions and minimum value ofμ for various edge supports. (C indicates clamped, S shows simply supported, and R is used for roller edge).
B.C.B.Cs atx=0B.Cs atx=Lμmin×LS-Su1=w1=Mx1=0u1=w1=Mx1=0πC-Cu1=w1=ϕ1=0u1=w1=ϕ1=02πC-Su1=w1=ϕ1=0u1=w1=Mx1=04.49341C-Ru1=w1=ϕ1=0u1=dw1dx=dMx1dx+Nx0dw1dx=0πS-Ru1=w1=Mx1=0u1=dw1dx=dMx1dx+Nx0dw1dx=0π2


## 5. Thermal Loading

Consider a beam under uniform temperature rise. That is, consider a beam at reference temperatureT0. The uniform temperature may be raised to T0+ΔT such that the beam buckles. Substituting (5) and (6) into (31) gives(32)cΔT∑n=1N∫zn-1znQ11nαxndz+c∑n=1Na(Q11)anVnd31n=μ2(D11-(B112/A11))1+(μ2/A55)(D11-(B112/A11)).


## 6. Numerical Result and Discussions

In this section, various combinations of composite beams comprising piezoelectric layers are assumed. General boundary conditions are considered on both sides to determine the critical buckling temperatures.


### 6.1. Aluminium Beam

Consider an aluminium beam with surface-bonded piezoelectric layers. we consider PZT-5A for piezoelectric layers. The beam thickness and length areh=0.01 m and L=0.25 m, and the actuator layer thickness is ha=0.001 m. The shear correction factor is k=5/6. Young's modules, coefficient of thermal expansion, Poisson's ratio, and the shear modules for aluminum are E=72.4 GPa,   α=22.5×10-6/∘C,  ν=0.3, and G=27.8 GPa, respectively [25]. The PZT-5A properties are Ea=63 GPa, αa=0.9×10-6/∘C,  ν12a=0.3, Ga=24.2 GPa, and d31=2.54×10-10 m/V [26]. Five electric loading cases are considered V0=0,  ±200 V,  ±500 V. Here, V0=0V denotes a grounding condition. Figure 2 and Table 4 depict the critical buckling temperature for various types of boundary conditions, and various voltages subjected to the uniform temperature rise. Also, the critical buckling temperature for the S-S and C-R types of boundary conditions are equal and larger than the value related to the S-R beams but lower than C-C and C-S beams.Table 4
Effect of applied voltage onΔTcr, (piezo aluminium beam).
Voltage (V)0200−200500−500S-S94.2829193.9278794.6379793.3952995.17054C-C370.19256369.83751370.54760369.30493371.08018C-S191.38436191.02931191.73940190.49673192.27198C-R94.2829293.92786794.6379693.3952995.17054S-R23.6817023.3266624.0367522.7940824.56932Critical buckling temperatures for piezoelectric aluminium beams with various boundary conditions and various voltages.
(a)(b)(c)(d)(e)


### 6.2. Glass-Epoxy Symmetric Beam

Consider a three-layered cross-ply composite beam (0/90/0), with surface-bonded piezoelectric layers. Also, similar to the previous example, consider PZT-5A for piezoelectric layers. The beam thickness and length areh=0.0045 m, and L=0.25 m and the actuator layer thickness is ha=0.001 m. The shear correction factor is k=5/6. It is assumed that the thickness and the material for all laminae are the same, (glass-epoxy) with the following characteristics [25]:(33)E11=50GPa,E22=15.2GPa,G12=G13=4.7GPa,G23=3.28GPa,α1=6×10-6/∘C,α2=α3=23.3×10-6/∘C,ν12=ν13=0.254,ν23=0.428.Figure3 and Table 5 depict the critical buckling temperature for various types of boundary conditions and various voltages subjected to the uniform temperature rise. The critical buckling temperature for the S-S and C-R types of boundary conditions are equal and larger than the value related to the S-R beams, but lower than C-C and C-S beams.Table 5
Effect of applied voltage onΔTcr (three-layered cross-ply composite beam).
Voltage (V)0200−200500−500S-S153.38994149.04308157.73679142.52281164.25707C-C606.25793601.91107610.60478595.39080617.12506C-S312.23341307.88656316.58027301.36629323.10054C-R153.38994149.04308157.73679142.52281164.25707S-R38.46406234.11721142.8109127.5969349.33119Critical buckling temperature for three layered cross-ply beams with various boundary conditions and various voltages.
(a)(b)(c)(d)(e)


### 6.3. Glass-Epoxy Antisymmetric Beam

Consider an antisymmetric four-layered composite beam (0/90/0/90), with surface-bonded piezoelectric layers. Similar to the previous examples, consider PZT-5A for piezoelectric layers. The beam thickness and length areh=0.004 m and L=0.25 m, and the actuator layer thickness is ha=0.001 m. The shear correction factor is k=5/6. The thickness and the material for all laminae are the same, (glass-epoxy), with material properties given in the previous example. In this example, we first consider one piezoelectric layer on the top surface of the beam, and then with two piezoelectric layers on the top and bottom surfaces of the beam.Figures4 and 5 and Tables 6 and 7 depict the critical buckling temperature for various types of boundary conditions and various voltages subjected to the uniform temperature rise. The critical buckling temperature for the S-S and C-R types of boundary conditions are equal and larger than the values related to the S-R beams but lower than the C-C and C-S beams.Table 6
Effect of applied voltage onΔTcr. (four-layered antisymmetric composite beam with one piezoelectric layer on the top surface of the beam).
Voltage (V)0200−200500−500C-C179.84154177.55059182.13249174.11417185.56892C-R45.2149242.9239747.5058739.4875450.94230Table 7
Effect of applied voltage onΔTcr. (four-layered antisymmetric composite beam with two piezoelectric layers on top and bottom surfaces of the beam).
Voltage (V)0200−200500−500C-C448.83477444.44855453.22102437.86920459.80038C-R113.24407108.85784117.63030102.27848124.20966Critical buckling temperature for antisymmetric four layered beam (0/90/0/90) with one piezoelectric layer on the top surface of the beam with various boundary conditions and various voltages.
(a)(b)Critical buckling temperature for antisymmetric four layered beam (0/90/0/90) with two piezoelectric layers on the top and bottom surfaces of the beam with various boundary conditions.
(a)(b)Figure6 depicts the difference between the buckling temperature for the four-layered antisymmetric beam (0/90/0/90) with one and two piezoelectric layers with various boundary conditions.Figure 6
Difference between the buckling temperature for four layered antisymmetric beam (0/90/0/90) with one and two piezoelectric layers with various boundary conditions.The results show that for this type of piezoelectric layer, the buckling temperature decreases with the increase of the applied voltage and increases with the increase of applied voltage in opposite phase. The changes are, however, small. It should be mentioned that increasing or decreasing the buckling temperature by applying voltage in comparison with the grounding condition depends upon both the sign of applied voltage and the sign of the piezoelectric constant.


### 6.4. Influence of Geometry on Critical Buckling Temperature

Consider three cross-ply composite beams with three layers (0/90/0) that are bonded with two piezoelectric layers on the top and bottom surfaces of the beams. The thickness of the beams areh=0.006 m, h=0.0045 m, and h=0.003 m. The lengths of the beams are equal and is L=0.25 m. The thickness and the material properties for all laminae are the same, (glass-epoxy), and the actuator layer is PZT-5A with thickness ha=0.001 m. The influence of beam geometry on the buckling temperature ΔTcr for various types of boundary conditions under applied voltages is shown in Figure 7. As shown, when the thickness increases, the critical buckling temperature increases for various types of boundary conditions, as expected.Figure 7
Influence of thickness on the buckling temperature.


## 6.1. Aluminium Beam

Consider an aluminium beam with surface-bonded piezoelectric layers. we consider PZT-5A for piezoelectric layers. The beam thickness and length areh=0.01 m and L=0.25 m, and the actuator layer thickness is ha=0.001 m. The shear correction factor is k=5/6. Young's modules, coefficient of thermal expansion, Poisson's ratio, and the shear modules for aluminum are E=72.4 GPa,   α=22.5×10-6/∘C,  ν=0.3, and G=27.8 GPa, respectively [25]. The PZT-5A properties are Ea=63 GPa, αa=0.9×10-6/∘C,  ν12a=0.3, Ga=24.2 GPa, and d31=2.54×10-10 m/V [26]. Five electric loading cases are considered V0=0,  ±200 V,  ±500 V. Here, V0=0V denotes a grounding condition. Figure 2 and Table 4 depict the critical buckling temperature for various types of boundary conditions, and various voltages subjected to the uniform temperature rise. Also, the critical buckling temperature for the S-S and C-R types of boundary conditions are equal and larger than the value related to the S-R beams but lower than C-C and C-S beams.Table 4
Effect of applied voltage onΔTcr, (piezo aluminium beam).
Voltage (V)0200−200500−500S-S94.2829193.9278794.6379793.3952995.17054C-C370.19256369.83751370.54760369.30493371.08018C-S191.38436191.02931191.73940190.49673192.27198C-R94.2829293.92786794.6379693.3952995.17054S-R23.6817023.3266624.0367522.7940824.56932Critical buckling temperatures for piezoelectric aluminium beams with various boundary conditions and various voltages.
(a)(b)(c)(d)(e)


## 6.2. Glass-Epoxy Symmetric Beam

Consider a three-layered cross-ply composite beam (0/90/0), with surface-bonded piezoelectric layers. Also, similar to the previous example, consider PZT-5A for piezoelectric layers. The beam thickness and length areh=0.0045 m, and L=0.25 m and the actuator layer thickness is ha=0.001 m. The shear correction factor is k=5/6. It is assumed that the thickness and the material for all laminae are the same, (glass-epoxy) with the following characteristics [25]:(33)E11=50GPa,E22=15.2GPa,G12=G13=4.7GPa,G23=3.28GPa,α1=6×10-6/∘C,α2=α3=23.3×10-6/∘C,ν12=ν13=0.254,ν23=0.428.Figure3 and Table 5 depict the critical buckling temperature for various types of boundary conditions and various voltages subjected to the uniform temperature rise. The critical buckling temperature for the S-S and C-R types of boundary conditions are equal and larger than the value related to the S-R beams, but lower than C-C and C-S beams.Table 5
Effect of applied voltage onΔTcr (three-layered cross-ply composite beam).
Voltage (V)0200−200500−500S-S153.38994149.04308157.73679142.52281164.25707C-C606.25793601.91107610.60478595.39080617.12506C-S312.23341307.88656316.58027301.36629323.10054C-R153.38994149.04308157.73679142.52281164.25707S-R38.46406234.11721142.8109127.5969349.33119Critical buckling temperature for three layered cross-ply beams with various boundary conditions and various voltages.
(a)(b)(c)(d)(e)


## 6.3. Glass-Epoxy Antisymmetric Beam

Consider an antisymmetric four-layered composite beam (0/90/0/90), with surface-bonded piezoelectric layers. Similar to the previous examples, consider PZT-5A for piezoelectric layers. The beam thickness and length areh=0.004 m and L=0.25 m, and the actuator layer thickness is ha=0.001 m. The shear correction factor is k=5/6. The thickness and the material for all laminae are the same, (glass-epoxy), with material properties given in the previous example. In this example, we first consider one piezoelectric layer on the top surface of the beam, and then with two piezoelectric layers on the top and bottom surfaces of the beam.Figures4 and 5 and Tables 6 and 7 depict the critical buckling temperature for various types of boundary conditions and various voltages subjected to the uniform temperature rise. The critical buckling temperature for the S-S and C-R types of boundary conditions are equal and larger than the values related to the S-R beams but lower than the C-C and C-S beams.Table 6
Effect of applied voltage onΔTcr. (four-layered antisymmetric composite beam with one piezoelectric layer on the top surface of the beam).
Voltage (V)0200−200500−500C-C179.84154177.55059182.13249174.11417185.56892C-R45.2149242.9239747.5058739.4875450.94230Table 7
Effect of applied voltage onΔTcr. (four-layered antisymmetric composite beam with two piezoelectric layers on top and bottom surfaces of the beam).
Voltage (V)0200−200500−500C-C448.83477444.44855453.22102437.86920459.80038C-R113.24407108.85784117.63030102.27848124.20966Critical buckling temperature for antisymmetric four layered beam (0/90/0/90) with one piezoelectric layer on the top surface of the beam with various boundary conditions and various voltages.
(a)(b)Critical buckling temperature for antisymmetric four layered beam (0/90/0/90) with two piezoelectric layers on the top and bottom surfaces of the beam with various boundary conditions.
(a)(b)Figure6 depicts the difference between the buckling temperature for the four-layered antisymmetric beam (0/90/0/90) with one and two piezoelectric layers with various boundary conditions.Figure 6
Difference between the buckling temperature for four layered antisymmetric beam (0/90/0/90) with one and two piezoelectric layers with various boundary conditions.The results show that for this type of piezoelectric layer, the buckling temperature decreases with the increase of the applied voltage and increases with the increase of applied voltage in opposite phase. The changes are, however, small. It should be mentioned that increasing or decreasing the buckling temperature by applying voltage in comparison with the grounding condition depends upon both the sign of applied voltage and the sign of the piezoelectric constant.


## 6.4. Influence of Geometry on Critical Buckling Temperature

Consider three cross-ply composite beams with three layers (0/90/0) that are bonded with two piezoelectric layers on the top and bottom surfaces of the beams. The thickness of the beams areh=0.006 m, h=0.0045 m, and h=0.003 m. The lengths of the beams are equal and is L=0.25 m. The thickness and the material properties for all laminae are the same, (glass-epoxy), and the actuator layer is PZT-5A with thickness ha=0.001 m. The influence of beam geometry on the buckling temperature ΔTcr for various types of boundary conditions under applied voltages is shown in Figure 7. As shown, when the thickness increases, the critical buckling temperature increases for various types of boundary conditions, as expected.Figure 7
Influence of thickness on the buckling temperature.


## 7. Conclusion

In this paper, the buckling analysis of composite beams with piezoelectric layers under various types of boundary conditions is investigated. Exact analytical solutions for the critical buckling temperature differences of beams are presented. The following are concluded.(1)
The buckling temperature difference for homogeneous, symmetric composite, and antisymmetric composite beams can be controlled by applying suitable voltage on the actuator layers, but the effect of this control voltage is small.(2)
For composite beams under uniform temperature rise, by increasing the beam thickness, the critical buckling temperature increases for any type of boundary conditions.



---

*Source: 362030-2011-04-10.xml*