(3) 78 795 Application of Laplace Iteration method to Study of Nonlinear Vibration of laminated composite plates Abstract In this paper, free vibration characteristics of laminated composite plates are investigated. A model is developed for a composite layer based on the consideration of non-linear terms in von-karman s non-linear deformation theory. The governing partial equation of motion is reduced to an ordinary non-linear equation and then solved using LIM method. The variation of frequency ratio of the Isotropic and composite plates is brought out considering parameters such as aspect ratio, fiber arrangements (orientation), number of layers, and modal ratios. Keywords Nonlinear analytical analysis; Nonlinear Vibration; Composite Plate, Classical Plate Theory; Laplace Iteration Method H. Rafieipour *, A. Lotfavar, A. Masroori, E. Mahmoodi Department of Mechanical and Aerospace Engineering, Shiraz University of Technology, 7555-33 Shiraz, Iran Received Nov In revised form 4 Dec * Author email: H.Rafieipour@sutech.ac.ir Lotfavar@sutech.ac.ir INTRODUCTION Laminated composite plates, due to their high specific strength and stiffness have been increasingly used in a wide range of aerospace, mechanical, and civil applications. By tailoring the sequence of the stacks and the thickness of the layers, composite laminates characteristics can be matched to the structural requirements with no difficulty. To use composite laminates efficiently, an accurate knowledge of vibration characteristics is essential. Vibration not only creates excessive noise and wastes energy but also may result in catastrophic failures. These phenomena, when the system operates around its natural frequencies, would be even more disastrous. Many publications have dealt with the linear vibrations of laminated composites. In these cases, the equation of motion is obtained easily and then by a reduction into a generalized eigenvalue problem, frequencies and mode shapes are determined. However, in many working conditions, plates are subjected to large amplitudes, so a nonlinear frequency analysis is required as a result. Obtaining an exact solution for the nonlinear free vibration of composite plates is very difficult due to the complexity of the equation of motion. The first approximate solution developed by Chu and Herman[] in 956, was the start point for so many other numerical methods introduced in the following years such as the finite element method (FEM)[], the discrete singular convolution method (DSC), the strip element method, and the Ritz methods[3].
78 H. Rafieipour et al. / Application of Laplace Iteration method to Study of Nonlinear Vibration of laminated composite plates Singh et al.[4] used direct numerical integration to study non-linear vibration of rectangular laminated composite plates in 99. Using Kirchhoff hypothesizes and von-karman strain-stress relations, they derived governing equations. They also employed harmonic oscillating assumptions and investigated large amplitude vibrations for various arrangements. Singha and Ganapathi[5], examined thin isotropic composite plates exploiting finite element method. They also considered the effects of shear deformation, rotational inertia and in plane inertia, in their formulation. In their study, non-linear matrices derived from Lagrange equations were solved using direct iteration technique and also they examined variations of non-linear frequency ratio versus various amplitudes considering influencing parameters, such as boundary conditions and fiber orientations. Malekzadeh and Karami[6], investigated non-linear vibrations of composite plates considering first order shear deformation theory employing Quadrature method. Effects of various parameters, such as foundations, thickness to length ratio and aspect ratio on frequency ratio were also examined in their study. Although numerical approaches are applicable to a wide range of practical cases, approximate analytical methods provide highly accurate solutions and a deep physical insight. One of the main approximate analytical approaches on nonlinear vibration analysis is Perturbation Method. This method is effective just in solving weakly nonlinear differential equations. Because of the limited application of the perturbation methods, newer approaches have been developed during recent years which are more powerful. For example, Pirbodaghi et al.[7] used the homotopy analysis method (HAM) to obtain an approximate analytical solution for nonlinear vibrations of thin laminated composite plates and showed effectiveness of this method to obtain an accurate solution with relatively less computational effort. Moreover, the energy balance method (EBM) [8], the max-min approach (MMA) [9], the homotopy perturbation method (HPM) [] and so many others have been proved to be strong approximate analytical approaches, obtaining accurate results for nonlinear differential equations []. The LIM method, which was introduced by Rafieipour et al.[] is a very powerful method in solving non-linear differential equations and its effectiveness is proved in studying non-linear vibration of composite beams. It was shown that this is one of the best analytical methods due to the rate of convergence and its accuracy. In the present paper, LIM method is used to obtain approximate analytical solutions for nonlinear vibrations of a thin laminated composite plate. A variational method is employed to derive the governing equations. The nonlinear term added is associated with von-karman strain displacement relation. LIM method is used to achieve nonlinear natural frequency and its excellent accuracy for the wide range of amplitude values is satisfied. A comparison with results in other articles has been done to validate the answers. The impact of layer arrangements and aspect ratio effects on the free vibration are investigated in detail. FORMULATION Consider a thin laminated composite plate of length a, width b, and thickness h as shown in figure. A positive set of coordinate system demonstrated in a way that reference surface is taken as the Latin American Journal of Solids and Structures (3) 78 795
H. Rafieipour et al. / Application of Laplace Iteration method to Study of Nonlinear Vibration of laminated composite plates 783 Latin American Journal of Solids and Structures (3) 78 795 mid plane and its origin is considered at the corner of the plate. The plate is supposed to be simply supported along all its edges. In addition, no slip condition between the layers is assumed. Figure Schematic of the laminated composite plate Assume u, v and w are the displacements of an arbitrary point of the plane in x, y and z axis, u and v to be the corresponding displacements of that point in the mid-plane and ε and κ to be the mid-plain strain and mid-plane curvature respectively; mechanical shear relation can be demonstrated as[3]: ε x ε y ε xy ε x ε y ε xy + z κ x κ y κ xy () Following von-karman s strain-displacement assumptions the in plane strain, shear strain and plane curvatures can be expressed as: ε x ε y ε xy u x + w x v y + w y u y + v x + w x w y () κ x κ y κ xy w x w y w x y (3)
784 H. Rafieipour et al. / Application of Laplace Iteration method to Study of Nonlinear Vibration of laminated composite plates According to Classical Theory of Elasticity, the strain-stress relations for each layer can be derived as: σ σ σ 6 k Q Q Q Q Q 66 k ε ε ε 6 k (4) Where the numbers,, 6 referred to principal axis of each layer. ( Q ij ) k (i,, 6) are the coefficients of the reduced stiffness matrix at the kth layer and are defined as: Q Q E ν ν, Q ν E ν ν ν E ν ν E ν ν, Q 66 G (5) Below stress-strain relations are obtained from the axis transformation of the each layer stressstrain equations referred to global axes: σ x σ y σ xy k Q Q Q 6 Q Q Q 6 Q 6 Q 6 Q 66 k ε x ε y ε xy + z κ x κ y κ xy k (6) And ( Q ij ) k (i,, 6) are plane stress-reduced stiffness coefficients. Constitutive equations which relate force and moment resultants to the strains through an appropriate integration along the thickness can be developed as: N x N y N xy A A A 6 A A A 6 A 6 A 6 A 66 ε x ε y ε xy + B B B 6 B B B 6 B 6 B 6 B 66 κ x κ y κ xy (7) Latin American Journal of Solids and Structures (3) 78 795
H. Rafieipour et al. / Application of Laplace Iteration method to Study of Nonlinear Vibration of laminated composite plates 785 M x M y M xy B B B 6 B B B 6 B 6 B 6 B 66 ε x ε y ε xy + D D D 6 D D D 6 D 6 D 6 D 66 κ x κ y κ xy (8) Where, A ij B ij D ij 3 n ( Q ij ) h k ( k h k ) k n ( Q ij ) h k k h k k n ( Q ij ) h 3 3 k k h k k (9) A ij, B ij, and D ij are called extensional stiffnesses, bending-extensional coupling stiffnesses and bending stiffness, respectively. The equations of motion are derived from Hamilton s principle as follows[4]: t δ U T d t () t by: Where U, the potential energy of the plate and T, the kinematic energy of the plate are given U σ ( xx ε xx + σ yy ε yy + σ xy ε xy )dz dx dy () V a b w dx T ρ i h i dy () Substituting Eqs.() and () into the Hamilton s equation of motion and with an appropriate integration along the thickness and then by equating δu, δ v, and δ w coefficients to zero, one obtains the below relations: Latin American Journal of Solids and Structures (3) 78 795
786 H. Rafieipour et al. / Application of Laplace Iteration method to Study of Nonlinear Vibration of laminated composite plates N x x + N xy y N y y + N xy x M x x + M xy x y + y M y + x N x w x + x + y N y N xy w y + y w y ρh w t N xy w x (3) By assuming simply supported boundary condition, the relations below are considered for displacement equations to satisfy the boundary conditions: u U (t)sin mπ x a v V (t)sin mπ x a w W (t)sin mπ x a Sin nπ y b nπ y Sin b nπ y Sin b (4) And U(t), V(t), and W(t) are the maximum displacements of plate center point along principal axes x, y, and z respectively. U(t) and V(t) can be expressed in terms of W(t) using first two equations of Eq. (3) and then by employing Galerkin method and substituting U(t) and V(t) in terms of W(t) into Eq. (3), the governing equation can be written as: d W(t) dt + α W(t)+ α V (t)+ α 3 W 3 (t) (5) Supposing W max as the maximum vibration amplitude of the plate center, the initial conditions of the center of the plate can be expressed as: W () W max, dw () dt (6) 3 DESCRIPTION OF THE PROPOSED METHOD Using the Laplace Transformation method, an analytical approximated technique is proposed to present an accurate solution for nonlinear differential equations. To clarify the basic ideas of proposed method consider the following second order differential equation, Latin American Journal of Solids and Structures (3) 78 795
H. Rafieipour et al. / Application of Laplace Iteration method to Study of Nonlinear Vibration of laminated composite plates 787 u(t) + N { u(t) } (7) with artificial zero initial conditions and N is a nonlinear operator. Adding and subtracting the term ω ut (), the Eq. (7) can be written in the form u(t) + ω u(t) L { u(t) } f ( u(t) ) (8) where L is a linear operator and ( ()) ω () { ()} f u t u t N u t (9) Taking Laplace transform of both sides of the Eq. (8) in the usual way and using the homogenous initial conditions gives ( s + ω )U (s) I f u(t) { } () where s and I are the Laplace variable and operator, correspondingly. Therefore it is obvious that U (s) I{ f ( u(t) )} G(s) () where Gs s + ω () Now, implementing the Laplace inverse transform of Eq. () and using the Convolution theorem offer where t (3) ut () f u() τ gt ( τ) dτ g(t) I { G(s) } sin ωt ω (4) Substituting Eq. (9) and (4) into (3) gives Latin American Journal of Solids and Structures (3) 78 795
788 H. Rafieipour et al. / Application of Laplace Iteration method to Study of Nonlinear Vibration of laminated composite plates t ut ( ω u( τ) N{ u( τ) }) sin ( ω( t τ) ) dτ ω (5) Now, the actual initial conditions must be imposed. Finally the following iteration formulation can be used [5] t u u + ω u τ + N u τ ω t τ dτ ( { }) sin ( ) (6) n n n ω Knowing the initial approximation u, the next approximations u, n > can be determined from previous iterations. Consequently, the exact solution may be obtained by using: n u lim u (7) n n The method proposed here can be applied in various non-linear problems. However, there is no need for any linearization and neither any small parameter. Also, the obtained approximate solutions converge quickly to the exact one. 4 IMPLEMENTATION OF THE PROPOSED METHOD Rewriting Eq(5) in the standard form of Eq. (8) results in the following equation: where d η() t dt ( η ) + ωη() t f () t (8) ( η ) ω η { η } { } + + f () t () t N (), t N η() t αη() t αη () t αη (), t 3 3 (9) Applying the proposed method, the following iterative formula is formed as: t η ( t) η ( t) f η ( τ) sin ω( t τ) dτ (3) n+ + n ω Eq. (8) will be homogeneous, when f ( η(t) ) has a value of zero. So, its homogeneous. Solution η (t) Acos(ωt) (3) Latin American Journal of Solids and Structures (3) 78 795
H. Rafieipour et al. / Application of Laplace Iteration method to Study of Nonlinear Vibration of laminated composite plates 789 is considered as the zero approximation for using in iterative Eq.(3). Expanding f ( η (τ )) f η (τ ) we have: α A + ω A 3 4 α 3 A3 cos(ωt) 4 α 3 A3 cos(3ωt) α A + cos(ωt) (3) Considering the relation: t ( ω cos(mωτ ) )sin( ω(t τ )) dτ cos(ωt) cos(mωt) ω (m ) t sin(ωt) ω m m (33) The coefficient of the term cos(ωt) in f ( η (τ )) should be vanished in order to avoid secular terms in subsequent iterations. As a result, the first approximation of the frequency can be expressed as: ω α + 3 4 α 3 A (34) Substituting Eq. (3) into (3) and eliminating the secular terms which is the coefficient of f η, results in: cos(ωt) in forcing function 3 η() t 3α A + 96ω A 3α3A cos( ωt) + 6αA cos( ωt) 96ω { 3 + 3α3A cos(3 ωt) 48αA } (35) This is the first approximation of η () t. Replacement of Eq. (35) in Eq. (3) and applying the procedure again results in the second approximation of η() t as ( ωt) ( ωt) ( ωt) ( ωt) ( ωt) ( ωt) Ι +Ι cos +Ι cos +Ι3cos 3 η() t 3 +Ι 4cos4ωt +Ι 5cos5ωt +Ι6cos6ωt 988864 A +Ι 7cos 7 + Ι 8cos 8 + Ι9cos 9 (36) Latin American Journal of Solids and Structures (3) 78 795
79 H. Rafieipour et al. / Application of Laplace Iteration method to Study of Nonlinear Vibration of laminated composite plates Where Ι i are given in Appendix. The coefficient of cos( ω t ) in forcing function assumed to be zero in order to evade the secular terms: ω + β ω + β ω + β ω + β (37) 8 6 4 6 4 where 5 β6 α 5 3A + αa α 3 3 4 3 3 5 β4 α 6 3A α α3a + αα 5 3+ α A αα A 3 3 β α α α α α α 3 3 3 4 8 3 3 7 3 5 6 5 β α 6 3A α α3a + α 9 α3a 3 3 3 3 6 5 5 3 3 79 3A + 5 3A + A 5 3A (38) Solving Eq. (37) estimates the value of ω for the actual natural frequency of the system. 5 RESULTS AND DISCUSSIONS In order to verify the precision of the suggested method, current results were compared with other articles. Table illustrates the frequency ratio of square and rectangular plates. It is observable that our results are in excellent agreement with the results provided by other references. It can be observed that the second displacement coefficient is zero for square plates. Table Frequency ratio of isotropic rectangular and square plates using various methods ( υ.3) W / h..4.6.8..4.6.8 LIM.5.967.56.343.499.673.85.4.364.436 ( a / b ) HAM.5.967.56.34.4988.6698.85.44.353.4347 Ref [4].8.89.743.937.437 ---.753 --- ---.88 LIM.85.9.37.389.554.74.9384.435.3543.5694 ( a / b ) HAM.85.9.36.388.5537.744.9376.45.353.5679 Ref [4].54.98.97.355.54 ---.877 --- ---.4798 Latin American Journal of Solids and Structures (3) 78 795
H. Rafieipour et al. / Application of Laplace Iteration method to Study of Nonlinear Vibration of laminated composite plates 79 In both cases, the difference between the non-linear frequency and linear frequency increased by increasing the amplitude of vibration. This difference increased even more when aspect ratio increased (figure ). Figure The effect of aspect ratio on the nonlinear to the linear frequency ratio In table, the Variation of frequency ratios versus non-dimensional amplitude ratio for symmetrical and non-symmetrical square plate arrangements are shown. The symmetrical arrangement plates have the same frequency ratios values. Table Frequency ratio of composite square plate with different arrangements using various methods E / E 4, G / E.5, ν.5) ( W / h.5.5.75.5.5.75 [ /9 / /9 ].575.575.634...388.438.435.483.6886.688.7679.973.9694.67.659.4.5739.574.8876.8857 3.79 [ /9 /9 / ].59.59.535.89.89.38.3874.387.47.63.67.669.887.889.574.563.355.44.448.734.735.8439 [9 / / /9 ].59.59.535.89.89.38.3874.387.47.63.67.669.887.889.574.563.355.44.448.734.735.8439 The ratio of linear to non-linear frequency for rectangular plates with aspect ratios of and 4 are presented in tables 3 and 4 for different arrangements. The [9///9] arrangement shows the least frequency ratio. For all arrangements excluding the above-mentioned, the frequency ratio increases when aspect ratio increases. Latin American Journal of Solids and Structures (3) 78 795
79 H. Rafieipour et al. / Application of Laplace Iteration method to Study of Nonlinear Vibration of laminated composite plates Table 3 Frequency ratio of composite rectangular plate with different arrangements using various methods E / E 4, G / E.5, ν.5, a / b ) ( W / h.5.5.75.5.5.75 [ /9 / /9 ].535.5.645.6.56.47.4797.46.495.774.7463.7787.78.53.3955.373.478.784.6935 3.45 3.5 3.976 [ /9 /9 / ].97.97.37.446.4458.4674.8538.853.8946.35.34.365.7795.7777 3.666 3.644 3.3634 3.769 3.759 4.6 4.59 4.3949 [9 / / /9 ].3.3.34.6.6.33.577.576.798.44.439.4593.64.69.855.848.8777.9.83.55.493.3396 Table 4 Frequency ratio of composite rectangular plate with different arrangements using various methods E / E 4, G / E.5, ν.5, a / b 4) ( W / h.5.5.75.5.5.75 [ /9 / /9 ].56.499.653.67.9.454.4865.467.4956.783.75.7863.96.64.478.386.433.7444.736 3.745 3.444 3.48 [ /9 /9 / ].657.657.77.5547.5544.5838.4.43.95.579.5694.6489 3.3 3. 3.6857 3.683 3.899 4.56 4.53 4.836 4.879 5. [9 / / /9 ].39.39.36.78.78.75.48.48.688.49.489.443.597.593.789.7885.8479.997.996.7.5.97 Tables 5 to 7 summarize comparison results for different modulus ratios. Table 5: Frequency ratio of composite rectangular plate with different arrangements using various methods E / E, G / E.5, ν.5) ( W / h.5.5.75.5.5.75 [ /9 / /9 ].58.58.556.95.94.3.3936.3933.435.636.63.696.895.8944.73.7.75.46.4597.756.7545.894 [ /9 /9 / ].478.478.498.784.783.97.3666.3664.39.598.593.634.843.8396.4.3.7.378.3769.6595.6578.7553 [9 / / /9 ].478.478.498.784.783.97.3666.3664.39.598.593.634.843.8396.4.3.7.378.3769.6595.6578.7553 Obviously, there is a greater frequency ratio for higher modulus ratios; however symmetric [9///9] shows an opposite behavior. Latin American Journal of Solids and Structures (3) 78 795
H. Rafieipour et al. / Application of Laplace Iteration method to Study of Nonlinear Vibration of laminated composite plates 793 Table 6 Frequency ratio of composite rectangular plate with different arrangements using various methods E / E, G / E.5, ν.5, a / b ) ( W / h [ /9 / /9 ] [ /9 /9 / ] [9 / / /9 ].5.5.75.5.5.75 LIM.5..447.794.4.39.698.93 HAM.53.59.4345.734.993.94.68.94 Ref [4].588..456.74.36.96 LIM.99.35.684.6.46.8769 3.34 3.798 HAM.99.3498.6835.6.466.875 3.98 3.77 Ref [4].986.3579.73.877.9 3.798 LIM.334.7.665.438.633.8396.586.85 HAM.334.7.664.4377.639.839.576.839 Ref [4].35.375.885.479.6795.373 Table 7 Frequency ratio of composite rectangular plate with different arrangements using various methods E / E, G / E.5, ν.5, a / b 4) ( W / h.5.5.75.5.5.75 [ /9 / /9 ].53.5.63.9.99.3.4645.447.4683.747.756.7455.478.48.3579.3348.364.674.65.994.973 3.4 [ /9 /9 / ].36.36.59.47.47.4459.89.8.857.58.57.37.784.768 3.98 3.897 3.788 3.6734 3.679 4.64 4.573 4.796 [9 / / /9 ].36.36.34.4.4.337.65.64.8.487.484.463.685.68.834.87.88.388.378.68.66.3444 It can be concluded from data given in tables above that frequency ratio increases with an increase in aspect ratio and modulus ratio and decreases when number of layers increases. It should be noted that added number of layers will result in a raise of the frequency ratio. 6 CONCLUSION In this study, a composite laminated plate model is established to verify the vibrational behaviors of composite plate. The von-karman s assumption and efficient approximate method (LIM) is employed to derive the nonlinear governing equation of motion. Analytical expressions are presented for nonlinear natural vibration analysis of composite laminated plate. Comparing with other Results, it is shown that the approximate analytical solutions are in very good agreement with the corresponding solutions. The effects of different parameters such as layer arrangements, skew angles and aspect ratio on the natural frequencies of the composite laminated plates are also studied. Results show that for all arrangements but [9///9], the frequency ratio increases when aspect ratio increases. The frequency ratio increases with an increase in aspect ratio and modulus ratio and decreases when number of layers increases. Furthermore, the symmetrical arrangement plates have the same frequency ratios values. Latin American Journal of Solids and Structures (3) 78 795
794 H. Rafieipour et al. / Application of Laplace Iteration method to Study of Nonlinear Vibration of laminated composite plates References [] Chu HN, Herrmann G. Plate Influence of large amplitudes on free flexural vibrations of rectangular plates. ASME J Appl Mech 956;3:53-54. [] Singha MK, Daripa R. Nonlinear vibration of symmetrically laminated composite skew plates by finite element method. International Journal of Non-Linear Mechanics 7;4(9)44-5. [3] Raju IS, Rao GV, Raju KK. Effect of longitudinal or inplane deformation and inertia on the large amplitude flexural vibrations of slender beams and thin plates. Journal of Sound and Vibration 976;49(3):45-4. [4] Singh S, Raju KK, Rao GV. Non-linear vibrations of simply supported rectangular cross-ply plates. Journal of Sound and Vibration 99;4():3-6. [5] Singha MK, Ganapathi M. Large amplitude free flexural vibrations of laminated composite skew plates. International Journal of Non-Linear Mechanics.4;39():79-7. [6] Malekzadeha P, Karami G. Differential quadrature nonlinear analysis of skew composite plates based on FSDT. Engineering Structures 6;8(9)p37-38. [7] Pirbodaghi T, Ahmadian MT, Fesanghary M. On the homotopy analysis method for nonlinear vibration of beams. Mechanics Research Communications 9;36():43-48. [8] Younesian D, Askari H, Saadatnia Z, KalamiYazdi M. Frequency analysis of strongly nonlinear generalized Duffing oscillators using He s frequency amplitude formulation and He s energy balance method. Comput Math Appl ;59:3-38. [9] Yue YS, Lu FM. The max_min approach to a relativistic equation. Computers and Mathematics with Applications 9;58: 3_33. [] He JH. Homotopy perturbation method for solving boundary value problems. Phys Lett 6;A 35:87-88. [] He JH. Some asymptotic methods for strongly nonlinear equations. International Journal of Modern Physics B 6;():4-99. [] Rafieipour H,Lotfavar A, Mansoori MH. New analytical approach to nonlinear behavior study of asymmetrically LCBs on nonlinear elastic foundation under Steady axial and thermal loading. Latin American Journal of Solids and Structures ;9(5):53-545. [3] Reddy J. Mechanics of Laminated Composite Plates Theory and Analysis: CRC, Boca Raton, 997. [4] Rao SS. Vibration of Continuous Systems: John Wiely & Sons, New Jersy, 7. [5] Hesameddini E, Latifizadeh H. Reconstruction of Variational Iteration Algorithms using the Laplace Transform. Int J Nonlinear Sci Numer Simul 9;():365-37. Latin American Journal of Solids and Structures (3) 78 795
H. Rafieipour et al. / Application of Laplace Iteration method to Study of Nonlinear Vibration of laminated composite plates 795 Appendix Ιi coefficients in the second approximation of deflection. 49 4 3 3 3 3 Ι 99943A Π+ Π3-Π4A + - 3+ - 5+ 4 3-4 A+ 5+ 4 3-4 6 Π Π Π Π Π Π Π Π Π 3 5 48 3 496 7 36 Ι 988864A + -5498Π -583584Π -48Π -77444Π + 766768Π A 4 3 5 3 4 ( 3 4 5 3 5 4) 3 3 4 ((-45793Π -37464Π ) Π + 69885Π + 57456Π ) A-4364Π 4 5 3 4 5 5 + 658Π + -595574784Π + 736Π Π -77444Π + 39944Π A + 3 ( 3 4) 5 7 3 Ι 48A Π -3 Π+Π A + - Π3+ -Π4- Π5Π3+Π4 A+ Π4+ Π5 Π3- Π4 64 4 6 5 3 ( ) Ι Π Π Π 854784 Π+ 3877Π -45584Π Π 3 3 5 4 3 3 ( 77444 +5498 5+6935 3) A + A +5868 Π5+48768Π 4 + -844Π -66976Π Π -5443 Π +37656Π A+5Π 3 3 4 5 4 3 5 4 5 49 3 4 3 Ι 4 64384AA Π3- Π3 Π5-6 Π4+ Π3A+ - Π4+ Π5Π3+ Π4 6 45 3 45 3 8 4 3 Π 5 Π 5 Ι 5-48384A Π3+ - Π5- Π4Π3-4Π5 A+ - Π4+ Π3+ 3 9 7 4 8 5 3 4 Ι 6 47456A Π3Π3+ Π5A- Π5+ Π4 6 8 9 Ι 688ΠA +8Π 3 Π+ Π A-945Π 4 7 3 5 5 3 5 Ι 384Π AΠ 8 5 3 Ι 89Π 4 9 5 where α A α α A α α A α α A α A Π Π Π Π Π ω ω ω ω ω 3 4 3 3 3 3,, 4 4 3, 4 4, 5 Latin American Journal of Solids and Structures (3) 78 795