Open access peer-reviewed chapter

Vibration Analysis of Laminated Composite Variable Thickness Plate Using Finite Strip Transition Matrix Technique and MATLAB Verifications

By Wael A. Al-Tabey

Submitted: September 19th 2013Reviewed: November 8th 2013Published: September 8th 2014

DOI: 10.5772/57384

Downloaded: 2531

1. Introduction

In the past, the model of thin plate on the elastic foundation was mainly used in structural applications. Currently, thin films of metal, ceramic or synthetic materials deposited on the surface of the structural parts of the electronic devices are used to improve their mechanical, thermal, electrical and tribological properties. These thin films of material are considered as thin plates and in these applications, the substrate of thin film can be simulated as an elastic foundation [1-2].

The laminated composite rectangular plate is very common in many engineering fields such as aerospace industries, civil engineering and marine engineering. The ability to conduct an accurate free vibration analysis of plates with variable thickness is absolutely essential if the designer is concerned with possible resonance between the plate and driving force [3].

Ungbhakorn and Singhatanadgid [4] investigated the buckling problem of rectangular laminated composite plates with various edge supports by using an extended Kantorovich method is employed.

Setoodeh, Karami [5] investigated A three-dimensional elasticity approach to develop a general free vibration and buckling analysis of composite plates with elastic restrained edges.

Luura and Gutierrez [6] studied the vibration of rectangular plates by a non-homogenous elastic foundation using the Rayleigh-Ritz method.

Ashour [7] investigated the vibration analysis of variable thickness plates in one direction with edges elastically restrained against both rotation and translation using the finite strip transition matrix technique.

Grossi, Nallim [8] investigated the free vibration of anisotropic plates of different geometrical shapes and generally restrained boundaries. An analytical formulation, based on the Ritz method and polynomial expressions as approximate functions for analyzing the free vibrations of laminated plates with smooth and non-smooth boundary with non classical edge supports is presented.

LU, et al [9] presented the exact analysis for free vibration of long-span continuous rectangular plates based on the classical Kirchhoff plate theory, using state space approach associated with joint coupling matrices.

Chopra [10] studied the free vibration of stepped plates by analytical method. Using the solutions to the differential equations for each region of the plate with uniform thickness, he formulated the overall Eigen value problem by introducing the boundary conditions and continuity conditions at the location of abrupt change of thickness. However this method suffers from the drawback of excessive continuity, as in theory the second and third derivatives of the deflection function at the locations of abrupt change of thickness should not be continuous.

Cortinez and Laura [11] computed the natural frequencies of stepped rectangular plates by means of the Kantorovich extended method, whereby the accuracy was improved by inclusion of an exponential optimization parameter in the formulation.

Bambill et al. [12] subsequently obtained the fundamental frequencies of simply supported stepped rectangular plates by the Rayleigh–Ritz method using a truncated double Fourier expansion.

Laura and Gutierrez [13] studied the free vibration problem of uniform rectangular plates supported on a non-homogeneous elastic foundation based on the Rayleigh–Ritz method using polynomial coordinate functions which identically satisfy the governing boundary conditions.

Harik and Andrade [14] used the “analytical strip method” to the stability analysis of uni-directionally stepped plates. In essence, the stepped plate is divided into rectangular regions of uniform thickness. The differential equations of stability for each region are solved and the continuity conditions at the junction lines as well as the boundary conditions are then imposed.

1.1. The chapter aims

This chapter presents the finite strip transition matrix technique (FSTM) and a semi-analytical method to obtain the natural frequencies and mode shapes of symmetric angle-ply laminated composite rectangular plate with classical boundary conditions (S-S-F-F). The plate has a uniform thickness in x direction and varying thickness h(y) in y direction, as shown in Figure 1. The boundary conditions in the variable thickness direction are simply supported and they are satisfied identically and the boundary conditions in the other direction are free and are approximated. Numerical results for simple-free (S-S-F-F) boundary conditions at the plate edges are presented. The illustrated results are in excellent agreement compared with solutions available in the literature, which validates the accuracy and reliability of the proposed technique.

Figure 1.

A rectangular laminated plate with variable thickness

2. Formulation

The equation of motion governing the vibration of rectangular plate under the assumption of the classical deformation theory in terms of the plate deflection W(x, y, t) is given by:

2MXx222MXYxy+2MYy2=ρh(y)2wot2

Where Wis the transverse deflection, ρ = the density per unit area of the plate and h(y) is the plate thickness at any point. The bending and the twisting moments in terms of displacements are given by:

MX=D112wox2D122woy22D162woxyMY=D122wox2D222woy22D262woxyMXY=D162wox2D262woy22D662woxy}

The flexural rigidities Dijof the plate are given by:

Dij=13h3(y)ho3k=1n[(Q¯ij)]k(hok3hok13),i,j=1,2,3,.........

Where hokis the distance from the middle-plane of the plate according to hoto the bottom of the hothlayer as shown in Figure 1. And Qijk¯are the plane stress transformed reduced stiffness coefficients of the lamina in the laminate Cartesian coordinate system. They are related to reduced stiffness coefficients of the lamina in the material axes of lamina Qijkby proper coordinate relationships they can be expressed in terms of the engineering notations as:

Qij=[Q11Q12Q13Q12Q22Q23Q13Q23Q66]=[E11(1υ12υ21)υ21E11(1υ12υ21)0υ21E11(1υ12υ21)E22(1υ21υ12)000G12]

Where E11, E22are the longitudinal and transverse young's moduli parallel and perpendicular to the fiber orientation, respectively and G12is the plane shear modulus of elasticity, υ12and υ21are the poisson's ratios. Thus, the governing partial differential equation of laminated composite rectangular plate with variable thickness as shown in Figure 1 is reduced to:

D114wox4+4D164wox3y+2(D12+2D66)4wox2y2+4D264woxy3+D224woy4=ρh(y)2wot2

Or in contraction form:

D11Wxxxx+4D16Wxxxy+2(D12+2D66)Wxxyy+4D26Wxyyy+D22Wyyyy=ρh(y)Wtt

The substitution of equation (3) into equation (6) given the governing Partial differential equation:

D11{2x2(h3(y)ho3W,xx)}+2(D12+2D66){2xy(h3(y)ho3W,xy)}+D16{2x2(h3(y)ho3W,xy)}+4D26{2y2(h3(y)ho3W,xy)}+D22{2y2(h3(y)ho3W,yy)}=moh(y)hoWtt

Equation (7) may be written as:

D11h3(y)ho3Wxxxx+(2(D12+2D66)ho3)h3(y)yWxxy+(2(D12+2D66)ho3)h3(y)Wxxyy+D16h3(y)ho3Wxxxy+(4D26ho32h3(y)y2)Wxy+4D26ho3h3(y)Wxyyy+8D26ho3h3(y)yWxyy+(D22ho32h3(y)y2)Wyy+D22ho3h3(y)Wyyyy+2D22ho3h3(y)yWyyy=moh(y)hoWtt

The equation of motion (8) can be normalized using the non-Dimensional variables ξ and η as follows :

ψ11a4Wξξξξ+2ψ2h3(η)1a2bh3(η)ηWξξη+2ψ21a2b2Wξξηη+ψ31a3bWξξξη+4ψ41ab3Wξηηη+1ab4ψ4h3(η)2h3(η)η2Wξη+8ψ4h3(η)1ab2h3(η)ηWξηη+1b21h3(η)2h3(η)η2Wηη+1b4Wηηηη+2h3(η)1b3h3(η)ηWηηη=moD22ho2h2(η)Wtt

Where β=abis the aspect ratio, ξ=xa, η=yb, ψ1=D11D22, ψ2=(D12+2D66)D22, ψ3=D16D22and ψ4=D26D22.

3. Method of solution

The displacement W(ξ,η,t)=W(ξ,η)eiωtcan be expressed in terms of the shape functionXi(ξ), chosen a prior; and the unknown function Yi(η)as:

W(ξ,η,t)=i=0NXi(ξ)Yi(η)eiωt

The most commonly used is the Eigen function obtained from the solution of beam free vibration under the prescribed boundary conditions at ξ=0 and ξ=1.

The free vibration of a beam of length a can be described by the non-Dimensional differential equation:

(2FT1μi(FR1FR2μi2)(μi4+FR1FT1)sinμi+2FT1μi(FR1+FR2)(μi4+FR1FT1)μicosμi)A3

Where EIis the flexural rigidity of the beam. The boundary conditions for free edges beam as shown in Fig. 2 are:

at ξ=0 and ξ=1

2Xi(ξ)ξ2=03Xi(ξ)ξ3=0}

Figure 2.

The two free edges beam strip in ξ-direction

In this paper, the beam shape function in ξ-direction is considered as a strip element of the plate and the flexural rigidity EIof the beam can replaced by (1υ2)D22and for υ=0.3, it can be just approximated by ED22. The solution of the beam equation is given as:

Xi(ξ)=A1sin(μiξ)+A2cos(μiξ)+A3sinh(μiξ)+A4cosh(μiξ)

One can obtain the following system of homogenous linear equations by satisfying the boundary conditions (12) at ξ=0 and ξ=1.

Φi=sinhμi+sinμicoshμicosμiXi(ξ)=sin(μiξ)sinh(μiξ)+1Φi(cos(μiξ)cosh(μiξ))}

The different value of μiare the roots of equation:

2cos(μi)cosh(μi)+cos2(μi)+sin2(μi)sinh2(μi)+cosh2(μi)sinh(μi)cosh(μi)+sin(μi)cosh(μi)sinh(μi)cos(μi)sin(μi)cos(μi)=0

The roots of equation (15) are represented in the recurrence form:

μi=(i+0.5)π, i= 0, 1, 2, 3, ..........

The substitution of equation (10) into equation (9), multiplying both sides by Xj(x)and after some manipulation, we can find:

i=0Nj=0Mβ4f3(η)Yi,ηηηη+2β3af1(η)f3(η)Yi,ηηη+(2ψ2β2f3(η)cijaij+8ψ4β2af1(η)f3(η)bijaij+β2a2f2(η)f3(η))Yi,ηη+(2ψ2βaf1(η)f3(η)cijaij+ψ3βf3(η)dijaij+4ψ4βa2f2(η)f3(η)bijaij+4ψ4β3f3(η)bijaij)Yi,η+(ψ1f3(η)eijaijλ2)Yi=0

Where λ2=moω2a4D22, f1(η)=1h3(η)h3(η)η, f2(η)=1h3(η)2h3(η)η2, f3(η)=ho2h2(η),

aij=01XiXjdξ,

bij=01XjXi,ξdξ,

cij=01XjXi,ξξdξ,

dij=01XjXi,ξξξdξ

and eij=01XjXi,ξξξξdξ.

From the orthogonality of the beam Eigen function, aij=eij=0for ij, this is true for all boundary conditions except for plates having free edges in the ξ-direction.

The system of fourth order partial differential equations in equation (17) can be reduced to a system of first order homogeneous ordinary differential equations:

ddη{Yk}ij=[Ai]k{Yk}ij

And after some manipulation, the governing differential equation (17) will become:

i=0Nj=0MEijYi////+(O1)ij(O0)ijYi///+(O2)ij(O0)ijYi//+(O3)ij(O0)ijYi/+(O4)ijλ2(O0)ijYi=0

Where the frame denotes differentiation with respect to η.

Where: (O0)ij=β4t1(η)Eij, (O1)ij=2β3at2(η)Eij, (O2)ij=(2ψ2β2t1(η)cijaij+8ψ4β2at2(η)bijaij+β2a2t3(η))(O3)ij=(2ψ2βat2(η)cijaij+ψ3βt1(η)dijaij+4ψ4βa2t3(η)bijaij+4ψ4β3t1(η)bijaij), (O4)ij=ψ1t1(η)eijaij,

[Eij]=i×jUnit matrix,

i= 0, 1, 2, 3, ……….,N, j= 0, 1, 2, 3, ……….,M

where the coefficients of the matrix [Ai]kin equation (18), in general, are functions of η and the Eigen value parameter λ. The vector Ykis given by:

Yk=[Y¯1Y¯2Y¯iY¯N]

Where:

Y¯i=[YiYi/Yi//Yi///]

Solving the above system of first order ordinary differential equations using the transition matrix technique yields, at any strip element (i) with boundaries (i-1) and (i) to,

{Yi}j=[Bi]j{Yi1}j

Where [Bi]jis called the transition matrix of the strip element (i), which can be obtained using the method of system linear differential equations of the strip element (i) in equation (18) (the exact solution of (ODE)).

Following the same procedure, the above boundary conditions (equations (12)) can be written. The simple boundary conditions at η=0 and η=1 as shown in Figure 3 are:

Figure 3.

The two edges clamped variable thickness beam strip in η-direction

The boundary conditions at η=0 and η=1 can be expressed as:

wo=0D12a22woξ2D22b22woη22D26ab2woξη=0}

Using the assumed solution, equation (10) the boundary conditions can be given by the following equations:

At η=0 and η=1

Yi=02Yiη2=i=0N2ψ4bijβaijYiηψ5cijβ2aijYi}

Or in contraction form:

Yi=0Yi//=i=0NCF1bijaijYi/CF2cijaijYi}

Where CF1=2ψ4β, CF2=ψ5β2, ψ5=D12D22

The solution is found using 2N initial vectors Y0at η=0. Equation (22) is applied across the stripped plate until the final end at η=1 is reached. Thus, 2N solutions Si, where i= 0, 1, 2, 3, ……….,N, can be obtained. The true solutions [S]can be written as a linear combination of these solutions [7]:

[S]=i=12NCiSi

Where Ciare arbitrary constants. These constants can be determined by satisfying 2N boundary conditions at η=1 [7]. The matrix [S]forms a standard Eigen value problem.

4. Numerical results and discussion

In this section, some numerical results are presented for symmetrically laminated, angle-ply variable thickness rectangular plate with simple support in the variable thickness direction and free in the other direction. The designation (S-S-F-F) means that the edges x=0, x=a, y=0, y=b are free, free, simple supported and simple supported respectively. The plates are made up of five laminates with the fiber orientations [θ, - θ, θ, - θ, θ] and the composite material is Graphite/Epoxy, of which mechanical properties are given in Table 1. The Eigen frequencies obtained are expressed in terms of non-dimensional frequency parameter λ=(ρhoω2a4/D22)1/2. To illustrate the solution, a plate with linear variable thickness, h(y)is used (see Appendix A).

h(η)=1+Δη

Where Δis the tapered ratio of plate given by Δ=(hbh0)/h0, (h0)is the thickness of the plate at η=0 and (hb)is the thickness of the plate at η=1. A convergence investigation is carried out for a uniform plate and for plate of variable thickness (Δ=0.5)with aspect ratio β=(0.5,1.0). By varying the harmonic numbers of the series solution in equation (10). The results are shown in Table 2. It is found that excellent agreement and stable and fast convergence can be achieved with only a few terms of series solution (N= 3 to 5).

MaterialE1, (GPa)E2, (GPa)G12, (GPa)υ12E2/ E1G12/ E1
Graphite/Epoxy1388.967.10.3250.8

Table 1.

Material properties of unidirectional composite

In order to validate the proposed technique, a comparison of the results with some results available for other numerical methods [15] for uniform laminated plates with simple support in the y-direction and free in the other direction. The first six natural frequencies of such uniform laminated plates are depicted in Table 2.

Δ = 0.0
Nλ1λ2λ3λ4λ5λ6
170.421270.7012140.4421173.5211180.6231235.6753
270.421270.7012140.4421173.5211180.6231235.6753
370.288270.5827140.2496173.2098180.2833235.3197
470.288270.5827140.2496173.2098180.2833235.3197
570.288270.5827140.2496173.2098180.2833235.3197
Ref*70.30270.604140.255173.218180.287235.322

Table 2.

Comparison of the first six natural frequencies of symmetric angle-ply uniform laminated square plates (θ=45), (β=1.0)

*Y.K. Cheung and D. Zhou [15].


Table 3 and Table 4 shows a convergence analysis of the first six frequencies parameters of symmetrically angle-ply five laminates [45/-45/45/-45/45] variable thickness plate with tapered ratio (Δ=0.5)and with aspect ratio β=(0.5,1.0)with simple support in the y-direction and free in the other direction (S-S-F-F).

Figure 4 and Figure 5 show the mode shapes of the first six fundamental frequencies of the above plate. Figure 4 and Figure 5 both are divided into two graphics. The first one shows the mode shapes of the plate in surface form and the other shows the mode shapes of the plate in surface contour form. All simulation results and graphics were obtained using MATLAB software.

Δ = 0.5
β = 0.5
Nλ1λ2λ3λ4λ5λ6
180.217782.5621155.9665188.6633194.6253251.7333
280.217782.5621155.9665188.6633194.6253251.7333
379.862582.0025155.3232188.1111194.1002251.2035
479.862582.0025155.3232188.1111194.1002251.2035
579.862582.0025155.3232188.1111194.1002251.2035

Table 3.

The first six frequencies parameter of S-S-F-F symmetrically angle-ply laminated [45/-45/45/-45/45] variable thickness plate (Δ=0.5), (β=0.5).

Δ = 0.5
β = 1.0
Nλ1λ2λ3λ4λ5λ6
172.757573.8666143.3334175.4963183.7825240.7621
272.757573.8666143.3334175.4963183.7825240.7621
372.119973.4444142.9019175.0024183.1121240.0159
472.119973.4444142.9019175.0024183.1121240.0159
572.119973.4444142.9019175.0024183.1121240.0159

Table 4.

The first six frequencies parameter of S-S-F-F symmetrically angle-ply laminated [45/-45/45/-45/45] variable thickness plate (Δ=0.5), (β=1.0)

Figure 4.

The mode shapes of the first six fundamental frequencies of the angle-ply symmetrically [45/-45/45/-45/45] laminated variable thickness rectangular plate with S-S-F-F edges, aspect ratioβ=a/b=0.5, tapered ratio Δ=0.5

Figure 5.

The mode shapes of the first six fundamental frequencies of the angle-ply symmetrically [45/-45/45/-45/45] laminated variable thickness rectangular plate with S-S-F-F edges, aspect ratio β=a/b=1.0, tapered ratio Δ=0.5.

5. Concluding remarks

A semi-analytical solution of the free vibration of angle-ply symmetrically laminated variable thickness rectangular plate with classical boundary condition (S-S-F-F) is investigated using the finite strip transition matrix technique (FSTM). The numerical results for uniform angle-ply symmetrically square plate with classical boundary condition (S-S-F-F) is presented and compared with some available results. The results agree very closely with other results available in the literature. It can be observed from Tables 2 and 3 that rapid convergence is achieved with small numbers of N in the series solution. Comparing to other techniques, the finite strip transition matrix (FSTM) proves to be valid enough in this kind of application. In all cases the FSTM method is easily implemented in a computer program a yields a fast convergence and reliable results. Also, the effect of the tapered ratio (Δ)and aspect ratio (β)on the fundamental natural frequencies and the mode shapes for five layers angle-ply symmetrically laminated variable thickness plates has been investigated for two cases of tapered ratio (uniform and variable thickness) and two cases of aspect ratio (square and rectangular). In fact the varying of the thickness and the increase the length (b) about a length (a) tend to increase the natural frequencies and the mode shapes of the laminated plate. The results from this investigation have been illustrated in the three dimensional surface contours for two different aspect ratios.

Appendix (A)

Plate thickness function

In this appendix the derivation of the relation of the plate thickness h(y) in y-direction as shown in the Figure 6 is given.

Figure 6.

The relation of the plate thickness h(y) in y-direction

By similarity between the triangles (ABG) and (ACF):

h(y)=ho(1+yc)

By similarity between the triangles (ABG) and (ADE):

hoc=hbc+b

From equations (28) and (29) the plate thickness relation is:

h(y)=ho+(hbho)by

Where h(y)=hoat y=0,

h(y)=hbat y=b,

h(y)=ho+(hbho)byat y=y,

and h(y)=hat ho=hb

Using the assumed solution, equation (10) The relation between the thickness of the plate h(y) can be given by the following equation:

h(η)=ho+(hbho)η

Appendix (B)

MATLAB code

Composite coefficients (function programs)

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Wael A. Al-Tabey (September 8th 2014). Vibration Analysis of Laminated Composite Variable Thickness Plate Using Finite Strip Transition Matrix Technique and MATLAB Verifications, MATLAB, Kelly Bennett, IntechOpen, DOI: 10.5772/57384. Available from:

Embed this chapter on your site Copy to clipboard

<iframe src="http://www.intechopen.com/embed/matlab-applications-for-the-practical-engineer/vibration-analysis-of-laminated-composite-variable-thickness-plate-using-finite-strip-transition-mat" />

Embed this code snippet in the HTML of your website to show this chapter

chapter statistics

2531total chapter downloads

1Crossref citations

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Image Processing with MATLAB and GPU

By Antonios Georgantzoglou, Joakim da Silva and Rajesh Jena

Related Book

Frontiers in Guided Wave Optics and Optoelectronics

Edited by Bishnu Pal

First chapter

Frontiers in Guided Wave Optics and Optoelectronics

By Bishnu Pal

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More about us