Buckling of Laminated and Functionally Graded Plates using Radial Basis Functions

Transcription

1 Paper 223 Buckling of Laminated and Functionally Graded Plates using Radial Basis Functions A.M.A. Neves, A.J.M. Ferreira, E. Carrera 3, M. Cinefra 3 C.M.C. Roque 2, C.M. Mota Soares 4 and R.M.N. Jorge Departamento de Engenaria Mecânica, 2 INEGI Faculdade de Engenaria, Universidade do Porto, Portugal 3 Departament of Aeronautics and Aerospace Engineering Politecnico di Torino, Italy 4 Instituto Superior Técnico, Lisbon, Portugal Civil-Comp Press, 22 Proceedings of te Elevent International Conference on Computational Structures Tecnology, B.H.V. Topping, (Editor, Civil-Comp Press, Stirlingsire, Scotland Abstract A yperbolic sine sear deformation teory is used for te buckling analysis of functionally graded plates, accounting for troug-te-tickness deformations. Te linearized buckling equations and boundary conditions are derived using Carrera s Unified Formulation and furter interpolated by collocation wit radial basis functions. Numerical results demonstrate te accuracy of te present approac. Keywords: plates and sells, collocation, mesless metods, functionally graded plates. Introduction A conventional functionally graded plate (FGM considers a continuous variation of material properties over te tickness direction by mixing two different materials []. Te material properties of te FGM plate are assumed to cange continuously trougout te tickness of te plate, according to te volume fraction of te constituent materials. Analysis of vibrations of FGM plates can be found in Batra and Jin [2], Ferreira et al. [3], Vel and Batra [4], Zenkour [5], and Ceng and Batra [6]. Te analysis of mecanical buckling of FGM structures is less common in te literature. It can be found in Najafizade and Eslami [7], Zenkour [5], Ceng and Batra [6], Birman [8], Javaeri and Eslami [9]. Typically, te analysis of FGM plates is performed using te clasical plate teory (CPT [, ], te first-order sear deformation teory (FSDT [2, 2, 3, 3] or iger-order sear deformation teories (HSDT [4, 5, 3, 6, 3]. Te FSDT gives acceptable results but depends on te sear correction factor wic is ard to find as it depends on many parameters. Tere is no need of a sear correction factor wen using a HSDT but linearized buckling equations are more complicated tan

2 tose of te FSDT. Te Unified Formulation proposed by Carrera [7, 8] made te implementation of suc teories easier. Te use of alternative-to-finite element metods for te analysis of plates, suc as te mesless metods based on collocation wit radial basis functions is atractive due to te absence of a mes and te ease of collocation metods. In recent years, radial basis functions (RBFs sowed excellent accuracy in te interpolation of data and functions. Kansa [9] introduced te concept of solving partial differential equations by an unsymmetric RBF collocation metod based upon te multiquadric interpolation functions. Te autors ave applied successfully te RBF collocation tecnique to te static and dynamic analysis of composite structures [2, 2, 22, 23, 24, 25, 26, 27, 28, 29]. Te present paper considers te tickness stretcing issue on te buckling analysis of FGM plates, by a mesless tecnique based on collocation wit radial basis functions. 2 Problem formulation Te problem consists of a rectangular sandwic plate of plan-form dimensions a and b and uniform tickness. Te co-ordinate system is taken suc tat te x-y plane coincides wit te midplane of te plate (z [, /2]. Te sandwic core is a ceramic material and skins are composed of a functionally graded material across te tickness direction. Te bottom skin varies from a metal-ric surface (z = = to a ceramic-ric surface wile te top skin face varies from a ceramic-ric surface to a metal-ric surface (z = 3 = /2 as illustrated in figure. Te volume fraction of te ceramic pase is obtained from a simple rule of mixtures as: ( p z V c =, z [, ] V c =, z [, 2 ] ( ( p z 3 V c =, z [ 2, 3 ] 2 3 were p is a scalar parameter tat allows te user to define te gradation of material properties across te tickness direction of te skins. Wit tis formulation te interfaces between core and skins disappear. Note tat te core of te present sandwic and any isotropic material can be obtained as a particular case of te power-law function by setting p =. Te volume fraction for te metal pase is given as V m = V c. Te sandwic may be symmetric or non-symmetric about te mid-plane as we may vary te tickness of eac face. Figure 2 sows a non-symmetric sandwic wit volume fraction defined by te power-law ( for various exponents p, in wic top skin tickness is te same as te core tickness and te bottom skin tickness is twice te core tickness. Suc tickness relation is denoted as 2--. A bottom-core-top notation is used. -- means tat skins and core ave te same tickness. Te sandwic plate is subjected to compressive in-plane forces acting on te midplane of te plate. Nxx and N denote te in-plane loads perpendicular to te edges 2

3 3 = /2 z METAL 2 CERAMIC CERAMIC CERAMIC = -/2 METAL Figure : Sandwic wit isotropic core and FGM skins. / z coordinate V c Figure 2: Illustration of a 2-- sandwic wit FGM skins for several volume fractions. 3

4 y N N xx x N N yx Figure 3: Rectangular plate subjected to in-plane forces. x = and y = respectively, and N denotes te distributed sear force parallel to te edges x = and y = respectively (see fig A quasi-3d yperbolic sine plate sear deformation teory 3. Displacement field Te present teory is based on te following displacement field: u(x,y,z,t = u (x,y,t + zu (x,y,t + sin v(x,y,z,t = v (x,y,t + zv (x,y,t + sin u Z (x,y,t (2 v Z (x,y,t (3 w(x,y,z,t = w (x,y,t + zw (x,y,t + z 2 w 2 (x,y,t (4 were u, v, and w are te displacements in te x, y, and z directions, respectively. u, u, u Z, v, v, v Z, w, w, and w 2 are te unknowns to be determined. A constant term is assumed for te transverse displacement component instead of (4 (w = w to investigate te effect of te tickness stretcing on te buckling load. 4

5 3.2 Strains Te strains can be related to te displacement field as: ǫ xx ǫ γ = u + ( w 2 x 2( x 2 v + w 2 u + v + w x x w, γ γ ǫ = u + w z x v + w z w z (5 By substitution of te displacement field in (5, te strains are obtained: ǫ xx ǫ γ = ǫ ( xx ǫ ( γ ( + ǫ (nl xx ǫ (nl γ (nl + z ǫ ( xx ǫ ( γ ( + sin ǫ (Z xx ǫ (Z γ (Z (6 γ γ ǫ = γ ( γ ( ǫ ( +z γ ( γ ( ǫ ( +z2 γ (2 γ (2 ǫ (2 +π cos γ (Z γ (Z ǫ (Z (7 being te strain components obtained as ǫ ( xx ǫ ( γ ( ǫ ( xx ǫ ( γ ( = = γ ( γ ( ǫ ( = γ (2 γ (2 ǫ (2 u x v u + v x u x v u + v x = ; ; u + w x v + w w w 2 x w 2 ; ; ǫ (nl xx ǫ (nl γ (nl ǫ (Z xx ǫ (Z γ (Z = = γ ( γ ( ǫ ( γ (Z γ (Z ǫ (Z = = ( w 2 2( x 2 w 2 w w x u Z x v Z u Z + v Z x u Z v Z w x w 2w 2 (8 (9 ( ( are te non-linear terms tat will originate te linearized buckling equa- were ǫ (nl αβ tions. 5

6 3.3 Elastic stress-strain relations In te case of isotropic functionally graded materials, te 3D constitutive equations can be written as: σ xx σ τ τ τ σ = = C C 2 C 2 C C 44 C 2 C 2 ǫ xx ǫ γ ǫ xx ǫ γ + + C 2 C 2 C 44 C 44 C 33 γ γ ǫ γ γ ǫ (2 Te computation of te elastic constants C ij depends on wic assumption of ǫ we consider. If ǫ =, ten C ij are te plane-stress reduced elastic constants: C = E ν 2; C E 2 = ν ν 2; C 44 = G; C 33 = (3 were E is te modulus of elasticity, ν is te Poisson s ratio, and G is te sear modulus G = E. 2(+ν If ǫ (tickness stretcing, ten te elastic coefficients C ij are tose of te tree-dimensional stress state, given by C = E( ν2 3ν 2 2ν 3, C 2 = E(ν + ν2 3ν 2 2ν 3 (4 C 44 = G, C 33 = E( ν2 3ν 2 2ν 3 (5 3.4 Governing equations and boundary conditions Te governing equations of present teory are derived from te dynamic version of te Principle of Virtual Displacements. Te internal virtual work is δu = Ω { /2 [σ xx ( δǫ ( xx + zδǫ ( xx + sin + σ ( δγ ( + zδγ ( + sin + σ ( δγ ( + zδγ ( + z 2 δγ (2 ( δǫ xx (Z + σ δǫ ( + zδǫ ( + sin ( δγ (Z + σ + π cos δγ ( + zδγ ( + z 2 δγ (2 ( + σ δγ (Z δǫ ( + zδǫ ( δǫ (Z + π cos ] } dz dx dy (6 δγ (Z 6

7 δu = ( Nxx δǫ ( xx + M xx δǫ ( xx + Rxxδǫ Z xx (Z + N δǫ ( + M δǫ ( + Rδǫ Z (Z Ω + N δγ ( + M δγ ( + Rδγ Z (Z + Q δγ ( + M δγ ( + Rδγ 2 (2 + Q δγ ( + M δγ ( + Rδγ 2 (2 + Rδγ Z (Z + Q δǫ ( + M δǫ ( were Ω is te integration domain on plane (x,y and N xx N N = /2 σ xx σ σ dz, Q Q Q = /2 σ σ σ + Rδγ Z (Z dx dy (7 dz (8 M xx M M = /2 z σ xx σ σ dz, M M M = /2 z σ σ σ dz (9 R Z xx R Z R Z = Ω /2 sin { R 2 R 2 σ xx σ σ } = dz, /2 { R Z R Z z 2 { σ σ } = /2 π { cos σ σ (2 } dz. (2 Te external virtual work due to external loads applied to te plate is given as: δv = (p x δu + p y δv + p z δwdx dy Ω ( ( = (p x δu + zδu + sin δu Z + p y δv + zδv + sin δv Z + p z ( δw + zδw + z 2 δw 2 dx dy Te external virtual work due to in-plane forces and sear forces applied to te plate is given as: [ δv = Ω N xx w x δ w x + N w δ w x + N w δ w yx x + N w (22 } dz ] δ w dx dy (23 7

8 being N xx and N te in-plane loads perpendicular to te edges x = and y = respectively, and N and N yx te distributed sear forces parallel to te edges x = and y = respectively. Te virtual kinetic energy is given as: δk = Ω = /2 ρ ( uδ u + vδ v + ẇδẇdz dx dy [ I ( u δ u + v δ v + ẇ δẇ Ω + I ( u δ u + u δ u + v δ v + v δ v + ẇ δẇ + ẇ δẇ + I 2 ( u δ u + v δ v + ẇ δẇ 2 + ẇ δẇ + ẇ 2 δẇ (24 + I 3 (ẇ δẇ 2 + ẇ 2 δẇ + I 4 ẇ 2 δẇ 2 + I 5 ( u δ u Z + u Z δ u + v δ v Z + v Z δ v + I 6 ( u Z δ u Z + v Z δ v Z + I 7 ( u δ u Z + u Z δ u + v Z δ v + v δ v Z ] dx dy were te dots denote te derivative wit respect to time t and te inertia terms are defined as I i = /2 ρz i dz i =,, 2, 3, 4 (25 I 5 = /2 ρ sin dz; I 6 = /2 ρ sin 2 dz; I 7 = /2 ρz sin dz (26 Substituting δu, δv, and δk in te virtual work statement, integrating by parts wit respect to x, y, and t and collecting te coefficients of δu, δu, δu Z, δv, δv, δv Z, 8

9 δw, δw, δw 2, te following governing equations are obtained: δu : N xx x N = δv : N x N = /2 /2 { ( } ρ ü + zü + sin ü Z + p x dz { ( } ρ v + z v + sin v Z + p y dz δw : Q x Q + N xx 2 w x 2 + N 2 w x + N yx 2 w x + N 2 w 2 = /2 δu : M xx x M δv : M x M δw : M x M { ρ (ẅ + zẅ + z 2 ẅ 2 + pz } dz + Q = + Q = + Q = δu Z : RZ xx x RZ + RZ = δv Z : RZ x RZ + RZ = δw 2 : R2 x R2 + 2M = /2 /2 /2 /2 /2 /2 { ( } ρz ü + zü + sin ü Z + zp x dz { ( } ρz v + z v + sin v Z + zp y dz { ρz (ẅ + zẅ + z 2 ẅ 2 + zpz } dz { ρ sin { ρ sin ( ü + zü + sin ( v + z v + sin { ρz 2 ( ẅ + zẅ + z 2 ẅ 2 + z 2 p z } dz ü Z + sin v Z + sin (27 } p x dz } p y dz Te mecanical boundary conditions are: 9

10 grid P Table : Convergence study for te uni-axial buckling load of a simply supported 2-2- sandwic square plate wit FGM skins and p = 5 case using te iger-order teory. δu : n x N xx + n y N = n x Nxx + n y N δv : n x N + n y N = n x N + n y N δw : n x Q + n y Q = n x Q + n y Q δu : n x M xx + n y M = n x Mxx + n y M δv : n x M + n y M = n x M + n y M δw : n x M + n y M = n x M + n y M (28 δu Z : n x Rxx Z + n y R Z = n x RZ xx + n y RZ δv Z : n x R Z + n y R Z = n x RZ + n y RZ δw 2 : n x R 2 + n y R 2 = n x R2 + n y R2 were (n x,n y denotes te unit normal-to-boundary vector. 4 Numerical examples Te iger-order plate teory and collocation wit RBFs are considered for te uniaxial and bi-axial buckling analysis of simply supported functionally graded sandwic square plates (a = b of type C wit side-to-tickness ratio a/ =. Te material properties considered are E m = 7E (aluminum for te metal and E c = 38E (alumina for te ceramic being E = GPa. Te non-dimensional parameter used is P = Pa2 3 E. An initial convergence study wit te iger-order teory was conducted for eac buckling load type using 3 2, 7 2, and 2 2 grids. In table te uniaxial case is sown for te 2-2- sandwic wit p = 5 and in table 2 te bi-axial case is presented for te -2- sandwic wit p =. In te following, results are obtained by considering a grid of 7 2 points, wic seems acceptable by te convergence study. In tables 3 and 4 te critical buckling loads obtained from te present approac wit ǫ and ǫ = are tabulated for various power-law exponents p and tickness ratios. Bot tables include results obtained from classical plate teory (CLPT, firstorder sear deformation plate teory (FSDPT, K = 5/6 as sear correction factor, Reddy s iger-order sear deformation plate teory (TSDPT [6], and Zenkour s

11 grid P Table 2: Convergence study for te bi-axial buckling load of a simply supported -2- sandwic square plate wit FGM skins and p = case using te iger-order teory. sinusoidal sear deformation plate teory (SSDPT [5]. Table 3 refers to te uni-axial buckling load and table 4 refers to te bi-axial buckling load. A good agreement between te present solution and references considered, specially [6] and [5] is obtained. Tis allow us to conclude tat te present iger-order plate teory is good for te modeling of simply supported sandwic FGM plates and tat collocation wit RBFs is a good formulation. Present results wit ǫ = approximates better references [6] and [5] tan ǫ as te autors use te ǫ = approac. Tis study also lead us to conclude tat te tickness stretcing effect as a strong influence on te buckling analysis of sandwic FGM plates as ǫ = gives iger fundamental buckling loads tan ǫ. Te isotropic fully ceramic plate (first line on tables 3 and 4 as te iger fundamental buckling loads. As te core tickness to te total tickness of te plate ratio (( 2 / increases te buckling loads increase as well. We may conclude tat te critical buckling loads decrease as te power-law exponent p increases. By comparing tables 3 and 4 we also conclude tat te bi-axial buckling load of simply supported sandwic square plate wit FGM skins is alf te uni-axial one for te same plate. In figure 4 te first four buckling modes of a simply supported 2--2 sandwic square plate wit FGM skins, p =.5, subjected to a uni-axial in-plane compressive load, using te iger-order plate teory and 7 2 grid is presented. Figure 5 presents te first four buckling modes of a simply supported 2-- sandwic square plate wit FGM skins, p =, subjected to a bi-axial in-plane compressive load. 5 Conclusions An application of a Unified formulation by a mesless discretization is proposed, based on a tickness-stretcing yperbolic sine sear deformation teory tat was implemented for te buckling analysis of functionally graded sandwic plates. Te present formulation was compared wit analytical, mesless or finite element metods and sowed very accurate results. Te effect of ǫ sowed to be significant in suc sandwic problems. Acknowledgments Ana M. A. Neves acknowledges support from FCT grant SFRH/BD/45554/28.

12 P p Teory CLPT FSDPT TSDPT [6] SSDPT [5] present ǫ present ǫ = CLPT FSDPT TSDPT [6] SSDPT [5] present ǫ present ǫ = CLPT FSDPT TSDPT [6] SSDPT [5] present ǫ present ǫ = CLPT FSDPT TSDPT [6] SSDPT [5] present ǫ present ǫ = CLPT FSDPT TSDPT [6] SSDPT [5] present ǫ present ǫ = Table 3: Uni-axial buckling load of simply supported plate of C-type using te igerorder teory and a grid wit 7 2 points. 2

13 P p Teory CLPT FSDPT TSDPT [6] SSDPT [5] present ǫ present ǫ = CLPT FSDPT TSDPT [6] SSDPT [5] present ǫ present ǫ = CLPT FSDPT TSDPT [6] SSDPT [5] present ǫ present ǫ = CLPT FSDPT TSDPT [6] SSDPT [5] present ǫ present ǫ = CLPT FSDPT TSDPT [6] SSDPT [5] present ǫ present ǫ = Table 4: Bi-axial buckling load of simply supported plate of C-type using te igerorder teory and a grid wit 7 2 points. 3

14 .5 eig = eig = eig = eig = Figure 4: First four buckling modes. Uni-axial buckling load of a simply supported 2--2 plate C-type, p =.5, a 7 2 points grid, and using te iger-order teory..5 eig = eig = eig = eig = Figure 5: First four buckling modes. Bi-axial buckling load of a simply supported 2-- plate C-type, p =, a 7 2 points grid, and using te iger-order teory. 4

15 References [] Y. Miyamoto, W.A. Kaysser, B.H. Rabin, A. Kawasaki, and R.G. Ford. Functionally Graded Materials: Design, Processing and Applications. Kluwer Academic Publisers, 999. [2] R.C. Batra and J. Jin. Natural frequencies of a functionally graded anisotropic rectangular plate. Journal of Sound and Vibration, 282(-2:59 56, 25. [3] A. J. M. Ferreira, R. C. Batra, C. M. C. Roque, L. F. Qian, and R. M. N. Jorge. Natural frequencies of functionally graded plates by a mesless metod. Composite Structures, 75(-4:593 6, September 26. [4] S. S. Vel and R. C. Batra. Tree-dimensional exact solution for te vibration of functionally graded rectangular plates. Journal of Sound and Vibration, 272:73 73, 24. [5] A.M. Zenkour. A compreensive analysis of functionally graded sandwic plates: Part 2 buckling and free vibration. International Journal of Solids and Structures, 42(8-9: , 25. [6] Z. Q. Ceng and R. C. Batra. Exact correspondence between eigenvalues of membranes and functionally graded simply supported polygonal plates. Journal of Sound and Vibration, 229: , 2. [7] M. M. Najafizade and M. R. Eslami. Buckling analysis of circular plates of functionally graded materials under uniform radial compression. International Journal of Mecanical Sciences, 44(2: , 22. [8] V. Birman. Buckling of functionally graded ybrid composite plates. Proceedings of te t Conference on Engineering Mecanics, pages 99 22, 995. [9] R. Javaeri and M.R. Eslami. Buckling of functionally graded plates under inplane compressive loading. ZAMM - Journal of Applied Matematics and Mecanics / Zeitscrift fr Angewandte Matematik und Mecanik, 82(4: , 22. [] Syang-Ho Ci and Yen-Ling Cung. Mecanical beavior of functionally graded material plates under transverse load part i: Analysis. International Journal of Solids and Structures, 43(3: , 26. [] Syang-Ho Ci and Yen-Ling Cung. Mecanical beavior of functionally graded material plates under transverse load part ii: Numerical results. International Journal of Solids and Structures, 43(3: , 26. [2] A. M. Zenkour. Generalized sear deformation teory for bending analysis of functionally graded plates. Appl Mat Modell, 3:67 84, 26. [3] Z. Q. Ceng and R. C. Batra. Deflection relationsips between te omogeneous kircoff plate teory and different functionally graded plate teories. Arcive of Mecanics, 52:43 58, 2. [4] L. F. Qian, R. C. Batra, and L. M. Cen. Static and dynamic deformations of tick functionally graded elastic plate by using iger-order sear and normal deformable plate teory and mesless local petrov-galerkin metod. Composites: Part B, 35: , 24. [5] A. J. M. Ferreira, R. C. Batra, C. M. C. Roque, L. F. Qian, and P. A. L. S. Martins. 5

16 Static analysis of functionally graded plates using tird-order sear deformation teory and a mesless metod. Composite Structures, 69(4: , 25. [6] J. N. Reddy. Analysis of functionally graded plates. International Journal for Numerical Metods in Engineering, 47: , 2. [7] E. Carrera. C reissner-mindlin multilayered plate elements including zig-zag and interlaminar stress continuity. International Journal of Numerical Metods in Engineering, 39:797 82, 996. [8] E. Carrera. Developments, ideas, and evaluations based upon reissner s mixed variational teorem in te modelling of multilayered plates and sells. Applied Mecanics Reviews, 54:3 329, 2. [9] E. J. Kansa. Multiquadrics- a scattered data approximation sceme wit applications to computational fluid dynamics. i: Surface approximations and partial derivative estimates. Computers and Matematics wit Applications, 9(8/9:27 45, 99. [2] A. J. M. Ferreira. A formulation of te multiquadric radial basis function metod for te analysis of laminated composite plates. Composite Structures, 59: , 23. [2] A. J. M. Ferreira. Tick composite beam analysis using a global mesless approximation based on radial basis functions. Mecanics of Advanced Materials and Structures, :27 284, 23. [22] A. J. M. Ferreira, C. M. C. Roque, and P. A. L. S. Martins. Analysis of composite plates using iger-order sear deformation teory and a finite point formulation based on te multiquadric radial basis function metod. Composites: Part B, 34: , 23. [23] A.J.M. Ferreira, C.M.C. Roque, R.M.N. Jorge, and E.J. Kansa. Static deformations and vibration analysis of composite and sandwic plates using a layerwise teory and multiquadrics discretizations. Engineering Analysis wit Boundary Elements, 29(2:4 4, 25. [24] A.J.M. Ferreira, C.M.C. Roque, and R.M.N. Jorge. Analysis of composite plates by trigonometric sear deformation teory and multiquadrics. Computers & Structures, 83(27: , 25. [25] A.J.M. Ferreira, R.C. Batra, C.M.C. Roque, L.F. Qian, and R.M.N. Jorge. Natural frequencies of functionally graded plates by a mesless metod. Composite Structures, 75(-4:593 6, 26. Tirteent International Conference on Composite Structures - ICCS/3. [26] A.J.M. Ferreira, C.M.C. Roque, and R.M.N. Jorge. Free vibration analysis of symmetric laminated composite plates by fsdt and radial basis functions. Computer Metods in Applied Mecanics and Engineering, 94(39-4: , 25. [27] A. J. M. Ferreira, C. M. C. Roque, and P. A. L. S. Martins. Radial basis functions and iger-order sear deformation teories in te analysis of laminated composite beams and plates. Composite Structures, 66(-4: , 24. Twelft International Conference on Composite Structures. [28] A.M.A. Neves, A.J.M. Ferreira, E. Carrera, C.M.C. Roque, M. Cinefra, R. M.N. 6

17 Jorge, and C.M.M. Soares. Bending of fgm plates by a sinusoidal plate formulation and collocation wit radial basis functions. Mecanics Researc Communications, In Press, Accepted Manuscript:, 2. [29] A.M.A. Neves, A.J.M. Ferreira, E. Carrera, C.M.C. Roque, M. Cinefra, R. M.N. Jorge, and C.M.M. Soares. A quasi-3d sinusoidal sear deformation teory for te static and free vibration analysis of functionally graded plates. Composites Part B, In Press, Accepted Manuscript:, 2. 7