Free vibration of thick annular sector plate on Pasternak foundation with general boundary conditions

Huimin Liu1 , Fanming Liu2 , Haoran Bai3 , Ranbing Yang4

1, 2College of Automation, Harbin Engineering University, Harbin, P. R. China

1, 3, 4College of Mechanical and Electrical Engineering, Qingdao Agricultural University, Qingdao, P. R. China

2Corresponding author

Journal of Vibroengineering, Vol. 18, Issue 3, 2016, p. 1692-1706. https://doi.org/10.21595/jve.2016.16717
Received 8 December 2015; received in revised form 12 March 2016; accepted 22 April 2016; published 15 May 2016

Copyright © 2016 JVE International Ltd. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Creative Commons License
Abstract.

In this paper, the free vibration of the thick annular sector plate on Pasternak foundation with general boundary conditions is presented by using the three-dimensional elasticity theory and the improved Fourier series method. The displacements of the thick annular sector plate are expanded to be the three-dimensional (3-D) Fourier cosine series supplemented with several auxiliary functions to eliminate all the relevant discontinuities of the displacements and their derivatives at the edges with general boundary conditions. This method can be universally applied to any kinds of classical boundaries, elastic boundaries and the combinations of both boundaries without any special change in the solution procedure. The excellent accuracy and reliability of the current solutions are demonstrated by numerical examples and comparison of the present results with those available in the literature. Furthermore, numerous new results for thick annular sector plates on Pasternak foundation with general elastic restrains are presented. In addition, comprehensive studies on the effects of the elastic restraint parameters, geometric parameters and elastic foundation coefficients are also reported.

Keywords: thick annular sector plate, three-dimensional elasticity theory, general boundary conditions, Pasternak elastic foundation, free vibration.

1. Introduction

Annular sector plates have wide applications in various engineering constructions as one of structure elements, such as aircraft, space vehicles and structures in military industries. However, the annular sector plate structures are often exposed to severe vibration conditions. Therefore, it’s of particular importance to have a good understanding of the vibration characteristics for the thick annular sector plate in the predesign process of the annular sector plate structures. So far, a huge amount of research efforts has been made on the vibration of sector plates by using different methods and plate theories. Ramakris and Kunukkasseril [1] presented a closed-form analytical solution to analyze the free vibration of an annular sector plate with simply supported radial edges. Irie et al. [2] performed the free vibration of ring-shaped polar-orthotropic sector plate with the classical boundary condition by using the Ritz method. Westmann [3] presented an exact and approximate method for free vibration of sector plates with clamped and simply supported boundary conditions. Wang et al. [4] extended the differential quadrature method to analyze the free vibration of thin sector plates with various sector angles and six combinations of classical boundary conditions. Aghdam et al. [5] studied the free vibration analysis of annular sector plates for the clamped support boundary condition by using the extended Kantorovich method. Harik et al. [6] used the finite-strip method to study the free vibration of sector plates with classical boundary conditions. Seok and Tiersten [7] discussed the free vibration of cantilevered annular sector plates by using the variational approximation procedure. Swaminadham et al. [8] compared the natural frequencies of annular sector plates obtained from the finite element method and the experiment. Li et al. [9] used the curve strip Fourier p-element method to study the free vibration of circular and annular sectorial thin plates with classical boundary conditions. The free vibration of annular sector thin plates with classical boundary conditions was studied by Mukhopadhyay [10, 11] and Srinivasan [12] respectively using the semi-analytical method and the integral equation technique. Shi et al. [13] applied the Rayleigh-Ritz method and the improved Fourier series technique to investigate the free vibration of circular, annular and sector thin plate with general boundary conditions. Kim and Dickinson [14] used one-dimensional (1-D) orthogonal polynomials and Liew and Lam [15] adopted two-dimensional (2-D) orthogonal polynomials as the admissible functions to study the free vibration of annular sector plates by the Rayleigh-Ritz method. Ramaiah and Vijayakumar [16] studied the free vibration of annular sector plates with simply supported radial edges by the combination of the Rayleigh-Ritz method and the coordinate transformation. Mizusawa [17] used the spline finite element to analyze the free vibration of annular sector plates with classical boundary condition. Houmat [18] applied the hierarchical finite element method to study the free vibration of annular sector plates. Mirtalaie and Hajabasi [19] performed free vibration analysis of thin functionally graded thin annular sector plates with classical boundary conditions by the differential quadrature method. Liew and Liu [20] studied the free vibration analysis of moderately thick annular sector plates with classical boundary conditions by the differential quadrature method based on the first-order shear deformation theory. Vaidyanathan et al. [21] applied the Rayleigh-Ritz method to study the static analysis of the annular sector mindlin plate. Mizusawa [22] used the finite element method to study the natural frequencies of thick annular sector plates. Leissa et al. [23] considered the effect of stress singularities on the vibration analysis of thick annular sector plates and presented the corner functions to improve the convergence of the numerical solutions. Xiang et al. [24] applied the Ritz method to study the free vibration of thick annular sector plates. Employing the Ritz method and three-dimensional elastic theory, Liew et al. [25] tackled the vibration problem of thick annular sector plates under various classical boundary conditions. Zhou et al. [26] presented the Chebyshev-Ritz method to analyze the three-dimensional free vibration of annular sector plates under various classical boundary conditions.

According to the literature review mentioned above, most of studies on the problem mentioned in the title are usually based on the classical plate theory and the first-order shear deformation plate theory. The results obtained from the three-dimensional elastic theory are very limited. In addition, it is observed that the existing methods or techniques are only suitable for a particular type of classical boundary conditions, i.e., simply-supported supports, clamped boundaries and free edges, which usually leads to constant modifications of the solution procedures and corresponding computation codes to adapt the different boundary cases. Moreover, to the best of our knowledge, there is no work reported on the free vibration behavior of thick annular sector plates on Pasternak foundation with general boundary conditions based on the three-dimensional elastic theory. The only work available in the open literature is that Bhattacharya and Bhowmic [27] used an approximate solution to study the bending analysis of an annular thin sector plate under the free boundary condition and resting on a Winkler-type elastic foundation based on the classical plate theory. Thus, it is desirable to develop a unified and efficient method which is capable of dealing with more complicated boundary conditions involving arbitrary classical boundaries, elastic edge restraints and the combinations of both conditions and structures resting on a Pasternak foundation.

Recently, a modified Fourier series technique proposed by Li [28, 29] is widely used in the vibration analysis of plates and shells with general boundary conditions by the Ritz method, e.g. [30-39]. Based on this, the authors apply the modified Fourier series technique together with the Rayleigh-Ritz method to provide a new three-dimensional exact solution for the free vibration of thick annular sector plates on the Pasternak foundation with general boundary conditions in this paper. The three-dimensional elastic theory is adopted to formulate the theoretical model. The general boundary conditions of the annular sector plate are achieved by the technique of applying the artificial springs on the boundaries. Regardless of the boundary conditions, each displacement of the annular sector plate is expressed as an improved Fourier series which is constructed by a three-dimensional (3-D) Fourier cosine series and several auxiliary functions employed to ensure and accelerate the convergence of the series expansions. All the unknown expanded coefficients are treated equally and independently as the generalized coordinates and obtained directly by using the Rayleigh-Ritz method. The good convergence of the present method is checked and the excellent accuracy and reliability are validated compared with the results presented by other contributors. Some new results of thick annular sector plates with elastic boundary conditions and resting on the elastic foundation are presented, which may serve as benchmark solution for future computational methods. In addition, the effects of the boundary conditions, geometric parameters and elastic foundation coefficients on vibration behaviors of the thick annular sector plate are also presented.

2. Theoretical formulations

2.1. Preliminaries

Consider a thick annular sector plate with uniform thickness h, inner radius R0, outer radius R1, width R (R=R1-R0) of plate in the radial direction and sector angle ϕ, as shown in Fig. 1. A coordinate system (r,θ,z) is also shown in the Fig. 1, which will be used in the analysis. u, v and w denote the displacement components in r, θ and z directions.

Fig. 1. Schematic diagram of the 3-D annular sector plates: a) the geometry and coordinates; b) the boundary restraining springs

 Schematic diagram of the 3-D annular sector plates:  a) the geometry and coordinates; b) the boundary restraining springs

According to the small deformation and linear strain-displacement relations, the strain components of thick annular sector plates can be expressed as:

(1)
ε r r = u r ,         ε θ θ = u r + 1 r v θ ,         ε z z = w z ,
γ r θ = v r + 1 r u θ - v r ,         γ r z = u z + w r ,         γ θ z = v z + 1 r w θ .

Based on the Hooke’s law, the corresponding stress-strain relations of the thick annular sector plates can be written as:

(2)
σ r r = Θ + 2 Ξ ε r r + Θ ε θ θ + Θ ε z z ,
σ θ θ = Θ ε r r + Θ + 2 Ξ ε θ θ + Θ ε z z ,
σ z z = Θ ε r r + Θ ε θ θ + Θ + 2 Ξ ε z z ,
τ r θ = Ξ γ r θ ,         τ r z = Ξ γ r z ,         τ θ z = Ξ γ θ z ,

where:

(3)
Θ = E μ ( 1 + μ ) ( 1 - 2 μ ) ,         Ξ = E 2 ( 1 + μ ) ,

in which E and μ are the Young’s modulus and Poisson’s ratio of the thick annular sector plates.

The boundary conditions for a generally supported thick sector plate can be expressed as the following forms based on the force equilibrium relationship on the four edges:

(4)
r = R 0 : σ r r + k r 0 u u = 0 ,         τ r θ + k r 0 v v = 0 ,         τ r z + k r 0 w w = 0 ,
r = R 1 : σ r r - k r 1 u u = 0 ,         τ r θ - k r 1 v v = 0 ,         τ r z - k r 1 w w = 0 ,
θ = 0 : τ r θ + k θ 0 u u = 0 ,         σ θ θ + k θ 0 v v = 0 ,         τ θ z + k θ 0 w w = 0 ,
θ = ϕ : τ r θ - k θ 1 u u = 0 ,         σ θ θ - k θ 1 v v = 0 ,         τ θ z - k θ 1 w w = 0 ,

in which kij means the restraining stiffness of the three groups of liner springs along the boundary edges, of which the superscripts and subscripts indicate the type of springs and the location of the corresponding springs, respectively. For example, kr0u is the stiffness for the liner boundary springs in r direction along the edge r=R0.

2.2. Solution procedure

The strain energy (U) of the thick annular sector plate during vibration can be define as:

(5)
U = 1 2 V ε r r σ r r + ε θ θ σ θ θ + ε z z σ z z + γ r θ τ r θ + γ r z τ r z + γ θ z τ θ z r d r d θ d z .

Substituting Eqs. (1)-(2) into Eq. (5), the strain energy expression of the thick annular sector plate can be written as follows:

(6)
U = 1 2 R 0 R 1 0 ϕ 0 h Θ + 2 Ξ u r 2 + u r + 1 r v θ 2 + w z 2 2 Θ u r + w z u r + 1 r v θ + 2 Θ u r w z + Ξ v r + 1 r u θ - v r 2 + u z + w r 2 + v z + 1 r w θ 2 r d r d θ d z .

The corresponding kinetic energy (T) function of the thick annular sector plate can be given as:

(7)
T = 1 2 R 0 R 1 0 ϕ 0 h ρ u t 2 + v t 2 + w t 2 r d r d θ d z ,

where ρ is the density of the plate and t donates time.

The potential energy stored in the boundary springs is given as:

(8)
U s p = 1 2 0 ϕ 0 h R 0 k r 0 u u   2 + k r 0 v v   2 + k r 0 w w 2 r = R 0 + R 1 k r 1 u u   2 + k r 1 v v   2 + k r 1 w w 2 r = R 1 d θ d z
            + 1 2 R 0 R 1 0 h k θ 0 u u   2 + k θ 0 v v   2 + k θ 0 w w 2 θ = 0 + k θ 1 u u   2 + k θ 1 v v   2 + k θ 1 w w 2 θ = ϕ d r d z .

As mentioned previously, this paper is focused on vibration analysis of the annular sector plate resting on the elastic foundation. As illustrated in Fig. 1, the elastic foundation is achieved by applying the two-parameter elastic foundation (Pasternak) mode with Winkler layer (stiffness Kw) and Shear layer (stiffness Ks). Therefore, the total energy stored by the foundation springs can be given by:

(9)
U F = 1 2 0 ϕ R 0 R 1 K w w 2 + K s w r 2 + K s 1 r w θ 2 z = 0 r d r d θ .

To determine the solutions, the Rayleigh-Ritz method is employed. The energy function defined by Lagrangian function is as:

(10)
L = U + U s p + U F - T .

In the Rayleigh-Ritz method, the solutions can be obtained by minimizing the energy function with respect to the coefficients of the admissible functions. So how to choose the proper admissible functions is the core of present method, for it acts as the most important factor in determining the accuracy and stability of the Rayleigh-Ritz method. It is expected that the displacement functions can be expressed in the form of Fourier series expansions because Fourier functions is of a complete set and exhibits good numerical stability. However, the conventional Fourier series expression generally have a convergence problem along the boundary edges except for a few simple boundary conditions. Thus, if not uniformly convergent, the derivatives of a Fourier series cannot be obtained simply through term-by-term differentiations. To overcome these problems, the displacement function will be here expressed as a more robust form of Fourier series expansion:

(11)
u s , θ , z , t = m = 0 n = 0 q = 0 A m n q c o s λ R m s c o s λ ϕ n θ c o s λ h q z + m = 0 n = 0 a m n 1 ζ z 1 z + a m n 2 ζ z 2 z c o s λ R m s c o s λ ϕ n θ + m = 0 q = 0 a m q 3 ζ θ 1 θ + a m q 4 ζ θ 2 θ c o s λ R m s c o s λ h q z + n = 0 q = 0 a n q 5 ζ s 1 s + a n q 6 ζ s 2 s c o s λ ϕ n θ c o s λ h q z e j ω t ,
(12)
v s , θ , z , t = m = 0 n = 0 q = 0 B m n q c o s λ R m s c o s λ ϕ n θ c o s λ h q z + m = 0 n = 0 b m n 1 ζ z 1 z + b m n 2 ζ z 2 z c o s λ R m s c o s λ ϕ n θ + m = 0 q = 0 b m q 3 ζ θ 1 θ + b m q 4 ζ θ 2 θ c o s λ R m s c o s λ h q z + n = 0 q = 0 b n q 5 ζ s 1 s + b n q 6 ζ s 2 s c o s λ ϕ n θ c o s λ h q z e j ω t ,
(13)
w s , θ , z , t = m = 0 n = 0 q = 0 C m n q c o s λ R m s c o s λ ϕ n θ c o s λ h q z + m = 0 n = 0 c m n 1 ζ z 1 z + c m n 2 ζ z 2 z c o s λ R m s c o s λ ϕ n θ + m = 0 q = 0 c m q 3 ζ θ 1 θ + c m q 4 ζ θ 2 θ c o s λ R m s c o s λ h q z + n = 0 q = 0 c n q 5 ζ s 1 r + c n q 6 ζ s 2 r c o s λ ϕ n θ c o s λ h q z e j ω t ,

where λRm=mπ/R,λϕn=nπ/ϕ,λhq=qπ/h and Amnq,Bmnq,Cmnq are the Fourier coefficients of three-dimensional Fourier series expansions for the displacements functions, respectively. amn1, amn2, amq3, amq4 , anq5, anq6 , bmn1, bmn2, bmq3, bmq4 , bnq5, bnq6 , cmn1, cmn2, cmq3, cmq4 , cnq5, cnq6 are the supplemented coefficients of the auxiliary functions. As mentioned before, introducing these supplementary terms into the Fourier series is to remove any potential discontinuities of the original displacements and their derivatives throughout the entire plate structure including the boundaries and then to effectively enhance the convergence of the results. According to the three-dimensional elasticity theory, it is required that at least two-order derivatives of the three displacement functions exist and are continuous at any point on the annular sector plates. It can be proven mathematically that the series expansion given in Eqs. (11)-(13) is able to expand and uniformly converge to any function Ξ(r,θ,z)C1 for V:0,R×0,ϕ×0,h. Also, the series can be easily differentiated, through term-by-term, to obtain uniformly convergent series expansions for up to the second-order derivatives. Since the series expression has to be truncated numerically, the proposed solution should be treated as a solution with arbitrary precision. In actual calculations, we truncate the infinite series expression to M, N, Q to obtain the results with acceptable accuracy.

Substituting Eqs. (6)-(9) and (11)-(13) into Eq. (10) and performing the Rayleigh-Ritz procedure with respect to each unknown coefficient, the equation of motion for plates can be yielded and is given in the matrix form:

(14)
K u u K u v K u w K u v T K v v K v w K u w T K v w T K w w - ω 2 M u u 0 0 0 M v v 0 0 0 M w w G u G v G w = 0 ,

where G is the vector of the unknown coefficients. K and M are the generalized stiffness and mass matrices of the thick sector plate, respectively. Obviously, the frequencies (or eigenvalues) of thick annular sector plate can be readily obtained by solving the Eq. (14) and by substituting the corresponding eigenvectors into series representations of displacement components the mode shapes are acquired then.

3. Numerical results and discussion

3.1. Convergence and boundary spring study

Table 1 shows the first six non-dimensional frequency parameters Ω=ωR12/ρh/D12/π2 (D=Eh3/121-μ2) of the simply supported thick annular sector plate with different truncation numbers M×N×Q. The geometry parameters of the annular sector plate are given as: R0/R1= 0.4, h/R= 1/6, ϕ= 90°. From Table 1, it is obvious that the current method has an excellent convergence, and is sufficiently accurate even when only a small number of terms are included in the series expression. When the truncated numbers are changed from 12×12×7 to 18×18×8, the maximum difference of the frequency parameters is within 0.1 % for the worst case, which can be accepted. Furthermore, by comparing the results with those published by Zhou et al. [26], we can find a very good agreement between the two results, although different displacement admissible functions and solution procedures were used in the literature. Moreover, more accurate results may be obtained by more truncated numbers, but the computational cost will increase meanwhile. Therefore, for the sake of both accuracy and computational cost, the truncated number of the displacement expressions will be uniformly selected as M×N×Q= 12×12×7 in the following numerical examples.

Table 1. Convergence of frequency parameters Ω for annular sector plate with different truncated number M×N×L

M × N × Q
Mode number
1
2
3
4
5
6
9×9×5
5.9442
7.9134
10.947
13.168
14.536
14.571
9×9×6
5.9434
7.9124
10.945
13.165
14.533
14.571
9×9×7
5.9424
7.9110
10.944
13.163
14.532
14.571
9×9×8
5.9418
7.9103
10.945
13.162
14.529
14.570
10×10×7
5.9392
7.9067
10.939
13.157
14.522
14.568
11×11×7
5.9368
7.9038
10.934
13.153
14.514
14.567
12×12×7
5.9324
7.9009
10.931
13.149
14.509
14.566
14×14×7
5.9316
7.8969
10.925
13.144
14.498
14.564
16×16×8
5.9308
7.8961
10.925
13.144
14.498
14.564
18×18×8
5.9305
7.8955
10.924
13.142
14.496
14.563
Ref. [26]
5.9306
7.8931
10.918
13.141
14.490
14.565

Fig. 2 shows the variation of the lowest three frequency parameters ΔΩ versus the elastic boundary restraint parameters Γλ for thick annular sector plates. The elastic boundary restraint parameters Γλ=Log10kλ/D (λ=u, v and w) are defined as ratios of corresponding spring stiffness to reference bending stiffness D. In addition, a frequency parameter ΔΩ which is defined as the difference of the non-dimensional frequency parameter Ω to those of the elastic boundary restraint parameters Γλλ=u,v,wto 10-8, i.e., ΔΩ=ΩΓλ-ΩΓλ=10-8 is used in the calculation. The geometric parameters of the annular sector plate used are: R0/R1= 0.5, h/R= 0.5, ϕ= 120o. The thick annular sector plate is free at edges θ= constant and under clamped restraint at r=R0, while at r=R1, the annular sector plate is elastically supported by only one group of spring component with stiffness varying from 10-8D to 1010D. It is obvious that in certain range the increase of the elastic restraint parameter leads to the increase of the non-dimensional frequency parameter. Thus, the classical boundary condition and general elastic edges can be obtained by setting the stiffness of springs to proper values. Taking edge r=R0 for example, the corresponding spring stiffness parameters for four types of classical boundary conditions and three types of elastic boundary conditions which are commonly encountered in engineering practices are given as follows: Clamped boundary condition (C): kr0u= 108D, kr0v= 108D, kr0w= 108D; Free boundary condition (F): kr0u= 0, kr0v= 0, kr0w= 0; Hard Simply-support boundary condition (S): kr0u= 108D, kr0v= 0, kr0w= 108D; Soft simply-support boundary condition (S1): kr0u= 0, kr0v= 0, kr0w= 108D; E1 type elastic boundary condition: kr0u= 102D, kr0v= 102D, kr0w= 0; E2 type elastic boundary condition: kr0u= 0, kr0v= 0, kr0w= 102D; E3 type elastic boundary condition: kr0u= 102D, kr0v= 102D, kr0w= 102D. The rationality of defining the boundary conditions in terms of boundary spring parameters will be proved by several examples in the next subsections. For concision, the boundary condition of thick sector plate is expressed by several simple letter strings, such as FCSE1 identifies the thick annular sector plate with F, C, S and E1 boundary conditions at boundaries r=R0, θ= 0, r=R1 and θ=ϕ.

Fig. 2. Variation of the frequency parameters ΔΩ versus the elastic restraint parameters Γλλ=u,v,w for thick sector plate

 Variation of the frequency parameters ΔΩ versus  the elastic restraint parameters Γλλ=u,v,w for thick sector plate
 Variation of the frequency parameters ΔΩ versus  the elastic restraint parameters Γλλ=u,v,w for thick sector plate
 Variation of the frequency parameters ΔΩ versus  the elastic restraint parameters Γλλ=u,v,w for thick sector plate

3.2. Thick annular sector plate with various boundary conditions

In the convergence study, the verification of the present method has been presented in the Table 1. But, it is confined to only one case of the boundary conditions and geometrical parameters. Thus, in this sub-section, the present formulations are applied to study the free vibration of thick sectors with various boundary conditions and geometrical parameters. In Table 2, the comparison of the frequency parameters Ω for annular sector plates with different boundary conditions is presented. The geometry parameters of the annular sector plates are given as: R0/R1= 0.4, ϕ= 90°, h/R= 1/3. Exact 3-D elasticity results provided by Liew et al. [25] and Zhou et al. [26] are involved in the comparison. In addition, the numerical results produced from the ABAQUS based on the FEA method (mesh size:0.02×0.02×0.02) are also shown in Table 2. From the table, we can see that there is a good agreement between the present results and the referential data. The discrepancy is very small and does not exceed 1.79 % for the worst case. The small discrepancy may be attributed to the different solution approaches used. As mentioned before, the mode shapes of the thick annular sector plates can be easily obtained by substituting the corresponding eigenvectors into series representations of displacement components. The comparison of the mode shapes of the thick annular sector plates using the present method and ABAQUS based on FEM is shown in Fig. 3. From the figure, it can be seen that the agreement of the results between these two methods is very good.

Table 2. Comparison of frequency parameters Ω for annular sector plates with various boundary conditions (R0/R1= 0.4, θ= 90°, h/R= 1/3)

Boundary
conditions
Method
Mode sequence
1
2
3
4
5
6
SSSS
Present
2.9972
3.0422
4.5866
4.9811
5.5974
6.7042
Ref. [25]
2.9973
3.0424
4.5870
4.9818
5.5973
6.7040
Ref. [26]
2.9973
3.0424
4.5870
4.9818
5.5973
6.7040
SSFF
Present
1.5619
2.3559
2.7834
2.9136
4.0291
4.6820
Ref. [25]
1.5608
2.3558
2.7829
2.9136
4.0289
4.6820
Ref. [26]
1.5587
2.3556
2.7827
2.9121
4.0288
4.6812
CCFF
Present
3.8822
4.0323
4.7494
5.2950
6.2421
6.7243
Ref. [25]
3.8839
4.0336
4.7411
5.2979
6.2429
6.7281
Ref. [26]
3.8821
4.0318
4.7494
5.2949
6.2415
6.7243
FE M
3.8801
4.0283
4.7408
5.2939
6.2262
6.7218
CCCC
Present
4.3740
5.6925
7.2910
7.5622
8.5620
9.4663
Ref. [25]
4.3755
5.6945
7.2911
7.5643
8.5621
9.4678
Ref. [26]
4.3738
5.5926
7.2904
7.5617
8.5600
9.4663
FEM
4.3701
5.5884
7.2858
7.5534
8.5488
9.4610

Fig. 3. Comparison of the mode shapes for a thick annular sector plate (R0/R1= 0.5, h/R1= 1/3) by using the present method and ABAQUS

 Comparison of the mode shapes for a thick annular sector plate  (R0/R1= 0.5, h/R1= 1/3) by using the present method and ABAQUS

a) Annular sector plate with CCFF boundary condition by using present method

 Comparison of the mode shapes for a thick annular sector plate  (R0/R1= 0.5, h/R1= 1/3) by using the present method and ABAQUS

b) Annular sector plate with CCFF boundary condition by using ABAQUS

 Comparison of the mode shapes for a thick annular sector plate  (R0/R1= 0.5, h/R1= 1/3) by using the present method and ABAQUS

c) Annular sector plate with CCCC boundary condition by using present method

 Comparison of the mode shapes for a thick annular sector plate  (R0/R1= 0.5, h/R1= 1/3) by using the present method and ABAQUS

d) Annular sector plate with CCCC boundary condition by using ABAQUS

Table 3. Fundamental frequency parameters Ω for annular sector plates with various boundary conditions, thickness ratios and sectors angle (R0/R1= 0.5)

h R 1
ϕ
Boundary conditions
CCCC
CFCS
SFSC
CS1CS1
CE1CE1
SE2SE2
CE2CE2
E3E3E3E3
E2E3E2E3
1/3
30°
5.6742
3.5112
3.1694
3.5035
3.3599
2.8146
3.4943
1.6672
1.6657
60°
3.9846
3.4131
2.7038
3.7407
3.3659
2.6999
3.4812
1.3533
1.3494
120°
3.4976
3.3832
2.5976
3.4512
3.3707
2.6402
3.4170
1.1585
1.1498
210°
3.4090
3.3773
2.5803
3.4002
3.3738
2.6144
3.3937
1.0988
1.0905
270°
3.3953
3.3767
2.5789
3.3914
3.3747
2.6068
3.3880
1.0866
1.0789
1/2
30°
3.9274
2.4495
2.3414
2.3376
2.3625
2.3330
2.3336
1.8747
1.8691
60°
2.8068
2.3956
2.0360
2.4948
2.3642
2.1581
2.4922
1.4852
1.4724
120°
2.4642
2.3740
1.9731
2.4285
2.3661
2.0537
2.4147
1.2778
1.2523
210°
2.3942
2.3691
1.9711
2.3867
2.3671
2.0132
2.3846
1.2165
1.2192
270°
2.3832
2.3685
1.9651
2.3796
2.3674
2.0029
2.3782
1.2041
1.1775
1
30°
2.0323
1.2443
1.3632
1.1701
1.2411
1.7257
1.1697
1.6958
1.6742
60°
1.4723
1.2531
1.1850
1.2481
1.2396
1.3435
1.2476
1.2706
1.2401
120°
1.2962
1.2427
1.1601
1.2801
1.2392
1.2068
1.2782
1.1108
1.0501
210°
1.2556
1.2395
1.1553
1.2512
1.2392
1.1732
1.2508
1.0696
1.0535
270°
1.2482
1.2393
1.1457
1.2461
1.2391
1.1662
1.2461
1.0617
1.0386

As mentioned previously, most of the existing contributions were restricted to the annular sector plate subjected to a limited set of classical supports. Thus, numerous new results of thick annular sector plates with various boundary conditions and geometry parameters are presented in Table 3. From the Table 3, it is obvious that the frequency parameters of annular sector plates generally decrease with the sector angle ϕ increasing under the classical boundary conditions. However, for the case under the elastic boundary condition, the influence of sector angles on the frequency parameters of annular sector plates becomes more complex. Fig. 4 shows the variation of the first four frequency parameters Ω of the annular sector plate against the sector angle ϕ with different boundary conditions. It can be seen from Fig. 4 that the result under CE2CE2 is different from that with other boundary conditions. Under the CE2CE2 condition, the frequency parameters increase rapidly as the sector angle ϕ varies from 0° to 30° and then decrease monotonously as the sector angle ϕ keeps increasing. But for the rest of the boundary conditions, the frequency parameters monotonously decrease as the sector angle ϕ increases. Thus, according to the analysis above, the fundamental frequencies of thick annular sector plate are strongly influenced by boundary conditions and geometry parameters.

Fig. 4. The variation of frequency parameters Ω of thick annular sector plates (R0/R1= 0.5, h/R1= 0.5) with different sector angles: a) CCCC; b) SSSS; c) CE2CE2; d) E1E2E1E2

 The variation of frequency parameters Ω of thick annular sector plates (R0/R1= 0.5, h/R1= 0.5) with different sector angles: a) CCCC; b) SSSS; c) CE2CE2; d) E1E2E1E2
 The variation of frequency parameters Ω of thick annular sector plates (R0/R1= 0.5, h/R1= 0.5) with different sector angles: a) CCCC; b) SSSS; c) CE2CE2; d) E1E2E1E2
 The variation of frequency parameters Ω of thick annular sector plates (R0/R1= 0.5, h/R1= 0.5) with different sector angles: a) CCCC; b) SSSS; c) CE2CE2; d) E1E2E1E2
 The variation of frequency parameters Ω of thick annular sector plates (R0/R1= 0.5, h/R1= 0.5) with different sector angles: a) CCCC; b) SSSS; c) CE2CE2; d) E1E2E1E2

3.3. Thick annular sector plate resting on Pasternak foundations

Firstly, the accessibility of applying the present method to study the vibration analysis of thick annular sector plates resting on elastic foundations need checking. The frequency parameters Ω for the thick annular sector plates with Winkler foundation stiffness are listed in Table 4. The results obtained from ABAQUS based on the FEA method (mesh size:0.02×0.02×0.02) are tabulated in the table for comparison. The geometry parameters of the annular sector plates are: R0/R1= 0.5, ϕ= 90°, h/R1= 0.5. For general purpose, the non-dimensional foundation parameters are used for numerical results as follows: KW=Kw×R14/D and KS=Ks×R12/D. In Table 4, the Winkler foundation stiffness KW is set as 10, 100 and 1000. From the table a consistent agreement of present results and referential data is clear. The discrepancy is very small and doesn’t exceed 0.62 % for the worst case. In addition, the table shows that the increasing of Winkler foundation stiffness contributes to the increase of frequencies of the sector plate.

As mentioned before, the results of the thick annular sector plate resting on Winkler/Pasternak foundation are limited in the published literature. Thus, the authors will present some new results of thick sector plates resting on the Winkler/Pasternak foundation with various boundary conditions and foundation coefficients using the present method in the next example. These results can be served as the benchmark solutions for future computational methods in this field. The fundamental frequency parameters Ω for thick annular sector plates with different boundary conditions and elastic foundation coefficients are shown in Table 5.

Table 4. Comparison of frequency parameters Ω for annular sector plates with various boundary conditions and thickness ratios (R0/R1= 0.5, θ= 90°, h/R1= 0.5)

K W
Mode No.
CCCC
SSSS
CFCF
Present
FEM
Error (%)
Present
FEM
Error (%)
Present
FEM
Error (%)
10
1
2.5834
2.5754
0.32
1.2600
1.2601
0.00
2.3956
2.3929
0.12
2
3.0541
3.0424
0.39
2.0197
2.0191
0.03
2.4442
2.4400
0.19
3
3.1901
3.1852
0.16
2.2451
2.2421
0.13
2.5511
2.5462
0.20
4
3.7672
3.7479
0.51
2.6542
2.6532
0.03
2.7350
2.7335
0.07
100
1
2.7241
2.7148
0.34
1.2615
1.2617
0.01
2.5485
2.5462
0.14
2
3.1412
3.1357
0.18
2.0201
2.0191
0.05
2.5493
2.5479
0.08
3
3.2157
3.2022
0.43
2.3627
2.3556
0.31
2.5811
2.5770
0.14
4
3.8435
3.8355
0.21
2.6541
2.6532
0.03
2.8276
2.8138
0.49
1000
1
3.1545
3.1479
0.21
1.2649
1.2650
0.01
2.5516
2.5478
0.15
2
3.1774
3.1722
0.18
2.0212
2.0207
0.02
2.8683
2.8592
0.32
3
3.5032
3.5079
0.14
2.6422
2.6378
0.16
2.9306
2.9249
0.20
4
4.0165
4.0041
0.31
2.6541
2.6532
0.03
3.0724
3.0643
0.26

Table 5. Fundamental frequency parameters Ω for annular sector plates with various boundary conditions and elastic foundations coefficients

h R
K W
K S
Boundary conditions
CCCC
CFCS
SFSC
CS1CS1
CE1CE1
SE2SE2
E1E2E1E2
E1E3E1E3
E2E3E2E3
0.2
10
0
7.5804
7.1688
3.9167
7.4145
7.1302
3.7846
0.3665
0.3691
0.8004
10
7.8922
7.4653
4.4038
7.7283
7.4141
4.2609
0.5063
0.5064
1.3482
100
10.179
9.6151
7.3462
10.029
9.4916
7.1434
0.5233
0.5232
2.9215
100
0
7.6385
7.2306
4.0272
7.4736
7.1925
3.8994
1.0235
1.0243
1.2501
10
7.9478
7.5245
4.5023
7.7849
7.4737
4.3623
1.0862
1.0862
1.6532
100
10.222
9.6601
7.4041
10.072
9.5374
7.2032
1.0943
1.0943
3.0736
1000
0
8.1947
7.8205
4.2929
8.0411
7.7863
4.8962
3.1714
3.1722
3.2791
10
8.4825
8.0916
5.3853
8.3301
8.0457
5.2705
3.2224
3.2224
3.4445
100
10.635
10.099
7.9599
10.490
9.9838
7.7746
3.2256
3.2256
4.3054
0.5
10
0
4.5435
4.2766
3.0554
4.4225
4.3022
3.0418
0.4503
0.4708
1.1416
10
4.9677
4.6832
3.5515
4.8467
4.6866
3.5116
0.6457
0.6476
1.1824
100
6.6658
6.0255
4.7135
6.6375
6.3919
5.1375
0.7088
0.7088
2.4056
100
0
4.6366
4.3774
3.1784
4.5164
4.4025
3.1653
1.0305
1.0426
1.4802
10
5.0457
4.7672
3.6476
4.9243
4.7718
3.6095
1.1492
1.1502
1.5154
100
6.6985
6.0338
4.7145
6.6703
6.4277
5.1706
1.1883
1.1633
3.0536
1000
0
5.4060
5.1586
3.3396
5.2773
5.2112
3.9422
2.8205
2.8493
3.1352
10
5.7036
5.1698
4.3963
5.5505
5.4585
3.9443
3.0376
3.0425
3.1797
100
6.9935
6.0735
4.7191
6.9663
6.7705
5.4056
3.0715
3.0754
3.9706

From the tables, it is obvious that the variation of the coefficients has significant effect on frequency parameters of the sector plate. Besides, it also can be found that for the same elastic foundations coefficients, the variation scope of the frequency parameter Ω of the annular sector plate h/R is different with different thickness ratios. In order to further study the influence of foundation coefficients on vibration characteristics of the thick annular sector plate, the variation of the lowest three frequency parameters ΔΩ versus the Winkler foundation stiffness and shearing layer stiffness for thick sector plate is presented in Fig. 5. It can be easily observed that regardless of the boundary condition and geometric parameters of the sector plate, there exists a range of the elastic foundation coefficients in which the frequency parameter Ω increases while the elastic foundation parameter increases and out of which frequency parameter Ω nearly does not change while the elastic foundation parameter keeps increasing. In addition, the influence scope of the frequency parameters of the sector plate expands obviously as the thickness ratio h/R of the annular sector plate decreases. Based on the above analysis, the influence of the elastic foundation coefficients on vibration characteristic is mainly determined by the boundary conditions and geometric parameters.

Fig. 5. Variation of the frequency parameters ΔΩ versus the shearing layer stiffness and Winkler foundation coefficients for thick annular sector plate with different boundary conditions: a) CCCC; b) CFCF; c) E1E2E1E2; d) E2E3E2E3

 Variation of the frequency parameters ΔΩ versus the shearing layer stiffness and Winkler foundation coefficients for thick annular sector plate with different boundary conditions:  a) CCCC; b) CFCF; c) E1E2E1E2; d) E2E3E2E3
 Variation of the frequency parameters ΔΩ versus the shearing layer stiffness and Winkler foundation coefficients for thick annular sector plate with different boundary conditions:  a) CCCC; b) CFCF; c) E1E2E1E2; d) E2E3E2E3
 Variation of the frequency parameters ΔΩ versus the shearing layer stiffness and Winkler foundation coefficients for thick annular sector plate with different boundary conditions:  a) CCCC; b) CFCF; c) E1E2E1E2; d) E2E3E2E3
 Variation of the frequency parameters ΔΩ versus the shearing layer stiffness and Winkler foundation coefficients for thick annular sector plate with different boundary conditions:  a) CCCC; b) CFCF; c) E1E2E1E2; d) E2E3E2E3

4. Conclusions

In this paper, based on the three-dimensional elasticity theory, the free vibration characteristics of thick annular sector plate with general boundary conditions resting on elastic foundation are presented. Regardless of the boundary conditions, each of the annular sector plate displacements is constructed in the form of a standard Fourier cosine series supplemented with several auxiliary functions introduced to eliminate all the possible discontinuities of the original displacement function and its derivatives throughout the entire plate space including the boundaries at the edges and to accelerate the convergence of series representation. The excellent accuracy and reliability of the current solutions are demonstrated by numerical examples and the comparison with the results available in the literature and obtained by using the FEA method. Besides, numerous new results for thick annular sector plates on the Pasternak foundation with elastic boundary conditions are presented. The effects of elastic boundary restraint parameters, boundary conditions, foundation coefficients and sector angles on the frequency values are examined and discussed in detail. The following conclusions can be drawn:

1) For the thick annular sector plate with classical boundary conditions, there is a good agreement between the frequency parameters obtained by the present method and the published results. The results show that the present method is able to converge rapidly and owns high reliability and accuracy.

2) The vibration characteristics of the thick annular sector plates are related to not only the boundary conditions but also their self-properties. In addition, the influence of the elastic foundation coefficients on vibration characteristics is mainly determined by the boundary conditions and geometric parameters.

3) Compared with most of the existing computational methods, the present method can be universally applied to a variety of boundary conditions including all the classical cases, elastic restraints and the combinations of both conditions without the need of change in the solution procedure. Numerous new free vibration results for thick annular sector plates with different geometric parameters and elastic restraints are presented.

Acknowledgements

The authors would like to thank the anonymous reviewers for their very valuable comments. The authors also gratefully acknowledge the financial support from the National Natural Science Foundation of China (No. 60834005), National Nature Science Foundation of Shandong Province of China and QingDao Science and Technology Plan Projects under Funds of ZR2013CL021 and 14-2-3-54-nsh respectively.

References

  1. Ramakrishnan R., Kunukkasseril V. X. Free vibration of annular sector plates. Journal of Sound and Vibration, Vol. 30, Issue 1, 1973, p. 127-129. [Publisher]
  2. Irie T., Yamada G., Ito F. Free vibration of polar-orthotropic sector plates. Journal of Sound and Vibration, Vol. 67, Issue 1, 1979, p. 89-100. [Publisher]
  3. Westmann R. A Note on Free Vibrations of Triangular and Sector Plates. 2012. [Search CrossRef]
  4. Wang X., Wang Y. Free vibration analyses of thin sector plates by the new version of differential quadrature method. Computer Methods in Applied Mechanics and Engineering, Vol. 193, Issue 36, 2004, p. 3957-3971. [Publisher]
  5. Aghdam M., Mohammadi M., Erfanian V. Bending analysis of thin annular sector plates using extended Kantorovich method. Thin-Walled Structures, Vol. 45, Issue 12, 2007, p. 983-990. [Publisher]
  6. Harik I. E., Molaghasemi H. R. Analytical solution to free vibration of sector plates. Journal of Engineering Mechanics, Vol. 115, Issue 12, 1989, p. 2709-2722. [Publisher]
  7. Seok J., Tiersten H. Free vibrations of annular sector cantilever plates. Part 1: out-of-plane motion. Journal of Sound and Vibration, Vol. 271, Issue 3, 2004, p. 757-772. [Publisher]
  8. Swaminadham M., Danielski J., Mahrenholtz O. Free vibration analysis of annular sector plates by holographic experiments. Journal of Sound and Vibration, Vol. 95, Issue 3, 1984, p. 333-340. [Publisher]
  9. Yongqiang L., Jian L. Free vibration analysis of circular and annular sectorial thin plates using curve strip Fourier p-element. Journal of Sound and Vibration, Vol. 305, Issue 3, 2007, p. 457-466. [Publisher]
  10. Mukhopadhyay M. A semi-analytic solution for free vibration of annular sector plates. Journal of Sound and Vibration, Vol. 63, Issue 1, 1979, p. 87-95. [Publisher]
  11. Mukhopadhyay M. Free vibration of annular sector plates with edges possessing different degrees of rotational restraint. Journal of Sound and Vibration, Vol. 80, Issue 2, 1982, p. 275-279. [Publisher]
  12. Srinivasan R., Thiruvenkatachari V. Free vibration of annular sector plates by an integral equation technique. Journal of Sound and Vibration, Vol. 89, Issue 3, 1983, p. 425-432. [Publisher]
  13. Shi X., Shi D., Li W. L., Wang Q. A unified method for free vibration analysis of circular, annular and sector plates with arbitrary boundary conditions. Journal of Vibration and Control, Vol. 22, Issue 2, 2016, p. 442-456. [Publisher]
  14. Kim C., Dickinson S. On the free, transverse vibration of annular and circular, thin, sectorial plates subject to certain complicating effects. Journal of Sound and Vibration, Vol. 134, Issue 3, 1989, p. 407-421. [Publisher]
  15. Liew K., Lam K. On the use of 2-D orthogonal polynomials in the Rayleigh-Ritz method for flexural vibration of annular sector plates of arbitrary shape. International Journal of Mechanical Sciences, Vol. 35, Issue 2, 1993, p. 129-139. [Publisher]
  16. Ramaiah G., Vijayakumar K. Natural frequencies of circumferentially truncated sector plates with simply supported straight edges. Journal of Sound and Vibration, Vol. 34, Issue 1, 1974, p. 53-61. [Publisher]
  17. Mizusawa T. Application of the spline element method to analyze vibration of annular sector plates. Journal of Sound and Vibration, Vol. 149, Issue 3, 1991, p. 461-470. [Publisher]
  18. Houmat A. A sector Fourier p-element applied to free vibration analysis of sectorial plates. Journal of Sound and Vibration, Vol. 243, Issue 2, 2001, p. 269-282. [Publisher]
  19. Mirtalaie S., Hajabasi M. Free vibration analysis of functionally graded thin annular sector plates using the differential quadrature method. Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science, Vol. 225, Issue 3, 2011, p. 568-583. [Publisher]
  20. Liew K. M., Liu F. L. Differential quadrature method for vibration analysis of shear deformable annular sector plates. Journal of Sound and Vibration, Vol. 230, Issue 2, 2000, p. 335-356. [Publisher]
  21. Vaidyanathan S., Busby H., Houser D. A numerical approach to the static analysis of an annular sector Mindlin plate with applications to bevel gear design. Computers and Structures, Vol. 51, Issue 3, 1994, p. 255-266. [Publisher]
  22. Mizusawa T., Kito H., Kajita T. Vibration of annular sector Mindlin plates by the spline strip method. Computers and Structures, Vol. 53, Issue 5, 1994, p. 1205-1215. [Publisher]
  23. Leissa A., McGee O., Huang C. Vibrations of sectorial plates having corner stress singularities. Journal of Applied Mechanics, Vol. 60, Issue 1, 1993, p. 134-140. [Publisher]
  24. Xiang Y., Liew K., Kitipornchai S. Transverse vibration of thick annular sector plates. Journal of Engineering Mechanics, Vol. 119, Issue 8, 1993, p. 1579-1599. [Publisher]
  25. Liew K., Ng T., Wang B. Vibration of annular sector plates from three-dimensional analysis. The Journal of the Acoustical Society of America, Vol. 110, Issue 1, 2001, p. 233-242. [Publisher]
  26. Zhou D., Lo S., Cheung Y. 3-D vibration analysis of annular sector plates using the Chebyshev-Ritz method. Journal of Sound and Vibration, Vol. 320, Issue 1, 2009, p. 421-437. [Publisher]
  27. Bhattacharya A., Bhowmic K. Note on the bending of an annular sector plate resting on an elastic foundation. Journal of Structural Mechanics, Vol. 4, Issue 3, 1976, p. 321-325. [Publisher]
  28. Li W. L. Free vibrations of beams with general boundary conditions. Journal of Sound and Vibration, Vol. 237, Issue 4, 2000, p. 709-725. [Publisher]
  29. Li W. L. Comparison of Fourier sine and cosine series expansions for beams with arbitrary boundary conditions. Journal of Sound and Vibration, Vol. 255, Issue 1, 2002, p. 185-194. [Publisher]
  30. Jin G., Ma X., Shi S., Ye T., Liu Z. A modified Fourier series solution for vibration analysis of truncated conical shells with general boundary conditions. Applied Acoustics, Vol. 85, 2014, p. 82-96. [Publisher]
  31. Jin G., Ye T., Chen Y., Su Z., Yan Y. An exact solution for the free vibration analysis of laminated composite cylindrical shells with general elastic boundary conditions. Composite Structures, Vol. 106, 2013, p. 114-127. [Publisher]
  32. Jin G., Ye T., Jia X., Gao S. A general Fourier solution for the vibration analysis of composite laminated structure elements of revolution with general elastic restraints. Composite Structures, Vol. 109, 2014, p. 150-168. [Publisher]
  33. Jin G., Ye T., Ma X., Chen Y., Su Z., Xie X. A unified approach for the vibration analysis of moderately thick composite laminated cylindrical shells with arbitrary boundary conditions. International Journal of Mechanical Sciences, Vol. 75, 2013, p. 357-376. [Publisher]
  34. Shi D., Wang Q., Shi X., Pang F. A series solution for the in-plane vibration analysis of orthotropic rectangular plates with non-uniform elastic boundary constraints and internal line supports. Archive of Applied Mechanics, Vol. 85, Issue 1, 2015, p. 51-73. [Publisher]
  35. Jin G., Su Z., Shi S., Ye T., Gao S. Three-dimensional exact solution for the free vibration of arbitrarily thick functionally graded rectangular plates with general boundary conditions. Composite Structures, Vol. 108, 2014, p. 565-577. [Publisher]
  36. Jin G., Su Z., Ye T., Jia X. Three-dimensional vibration analysis of isotropic and orthotropic conical shells with elastic boundary restraints. International Journal of Mechanical Sciences, Vol. 89, 2014, p. 207-221. [Publisher]
  37. Jin G., Xie X., Liu Z. The Haar wavelet method for free vibration analysis of functionally graded cylindrical shells based on the shear deformation theory. Composite Structures, Vol. 108, 2014, p. 435-448. [Publisher]
  38. Wang Q., Shi D., Liang Q. Free vibration analysis of axially loaded laminated composite beams with general boundary conditions by using a modified Fourier-Ritz approach. Journal of Composite Materials, 2015. [Search CrossRef]
  39. Wang Q., Shi D., Liang Q., Shi X. A unified solution for vibration analysis of functionally graded circular, annular and sector plates with general boundary conditions. Composites Part B: Engineering, Vol. 88, 2016, p. 264-294. [Publisher]

Cited By

Journal of Composites Science
Ahmad Saood, Zain A. Khan, Mohd T. Parvez, Arshad H. Khan
2021