Published: 30 June 2018

Mathematical solutions for the flexural analysis of Mindlin’s first order shear deformable circular plates

Charles Chinwuba Ike1
1Department of Civil Engineering, Enugu State University of Science and Technology, Enugu State, Nigeria
Views 478
Reads 227
Downloads 3009

Abstract

In this work, the problem of first order shear deformable solid circular plate under transverse load was solved mathematically. The problem considered was assumed axisymmetric. The plate and loading were considered axisymmetric. The problem was defined as a boundary value problem of a system of differential equations of equilibrium in terms of the stress resultants and the stress – resultants – displacement relations. The set of equations were considered simultaneously to express them in variable separable form. The mathematical technique of separation of variables was then used to obtain solutions for the unknown generalised displacements. Specific problems of clamped edge plates and simply supported edge plates under uniformly distributed load and point load at the centre were considered and solved using the same technique of separation of variables. The mathematical expressions obtained showed that in all cases, the deflection was expressible in terms of flexural and shear components. The maximum deflection was found to occur at the plate centre as is expected from the symmetrical nature of the problem. The shear component of the transverse deflection was found to significantly increase with significant increase in the ratio of the plate thickness to the radius (h/r0).

1. Introduction

Circular plate structures have extensive applications in aerospace, civil, structural, geotechnical, naval, and mechanical engineering [1]. They can be subjected to distributed transverse loads, in-plane loads or dynamic loads. They are classified according to the ratio of their thickness, h to the diameter, D, [2, 3] as: thin plates (D/h> 100), moderately thick plates (20 <D/h< 100), and thick plates (D/h< 3). Circular plates are also classified according to their material properties as: laminated plates, anisotropic plates, orthotropic plates, homogeneous plates, heterogeneous plates and isotropic plates [4-8].

The determination of displacement and deflection functions in elastic circular plates subjected to symmetrical distribution of transverse loading is a solid mechanics/theory of elasticity problem frequently encountered in the analysis and design of structural elements and systems [9]. Plate problems are generally three dimensional (3D) problems of the mathematical theory of elasticity, and the three spatial coordinates are used in their formulation. Three-dimensional plate theories have their mathematical and analytical foundations in the mathematical theory of elasticity; and the solution for 3D problems are mathematical rigorous and demanding. Only few 3D problems have so far been solved [10-14].

Theories that have been used to describe the flexural problem of circular plates include: Kirchhoff-Love plate theory also called the Classical small deformation thin plate theory (CPT) [15-18], Mindlin plate theory [19], Reissner plate theory [20-22], Von Karman plate theory [23], Shimpi’s plate theory [24, 25], Levinson’s plate theory [26], Reddy’s plate theory [27, 28] and Modified/refined plate theories [29].

The Kirchhoff-Love plate theory is based on the Kirchhoff-Love’s hypotheses which are assumptions similar to the Bernoulli-Navier’s hypotheses used in the classical theory of thin (or slender) beams. The fundamental assumptions of the classical small deflection thin plate theory for homogeneous, isotropic elastic plates are [18]:

(i) Straight line initially normal to the middle surface to the plate remains straight and normal to the deformed middle surface of the plate and unchanged in length.

(ii) The transverse displacement w is assumed to be very small. This implies that the slope of the deflection surface is small and hence the square of the slope would be considered negligible in comparison with unity.

(iii) The normal stresses in the radial and tangential (circumferential) directions σrr and σθθ and the in-plane shear stress τrθ are assumed to be zero at the middle surface for small deflection cases (i.e. wh).

(iv) The transverse normal stress σzz is very small as compared to the other stress components and can be disregarded in the Hooke’s generalised stress-strain relations without introducing significant errors i.e. σzz(σrr, σθθ, τrθ).

(v) The middle plane remains unstrained after bending deformation and is a neutral surface (plane).

The above assumptions reduce the plate problem from a 3D mathematical theory of elasticity problem to a two dimensional (2D) approximate theory. Hence only two normal stresses σrr, σθθ and one in-plane transverse shear stress τrθ are the stresses used to define the plate problem. These stresses are functions of the two coordinate variables r, and θ.

The differential equations of equilibrium of circular plates are [6]:

1
Mrrr+1rMrθθ+Mrr-Mθθr-Qr=0,
2
Mrθr+1rMθθθ+2Mrθr-Qθ=0,
3
Qrr+1rQθθ+Qrr+q=0,

where Mrr are radial moments, Mθθ are circumferential moments, Mrθ is the twisting moment, Qr, Qθ are shear force in the r and θ directions, Qr is the radial shear force, Qθ is the circumferential shear force and q is the transverse load distribution.

The differential equation of equilibrium is expressed in terms of the deflection of the middle surface by the fourth order biharmonic equation:

4
4w=22w=qD,

where:

5
2=2r2+1rr+2θ2,

and 2 is the Laplacian operator in the polar coordinate system, 4 is the biharmonic operator. D is the modulus of flexural rigidity of the plate.

For cases of axisymmetric or rotationally symmetric loading, the load, internal stress resultants (bending moment and shear force distributions) and the displacements become functions of one coordinate variable only, namely, r. In such cases, each vertical plane (rz) of the plate will be a principal plane, and hence the twisting moment Mrθ will vanish [6]. Also, the shear force Qθ and the tangential displacement v will be zero. The differential equation of equilibrium of axisymmetric circular Kirchhoff-Love plate will simplify to [30, 31]:

6
dMrrdr+Mrr-Mθθr-Qr=0,
7
dQrdr+Qrr+q=0.

The moment-curvature (moment-deflection) equations become:

8
Mrr=-Dd2wdr2+μrdwdr,
9
Mθθ=-Dμd2wdr2+1rdwdr.

The radial shear force Qr is:

10
Qr=-Dddr2w.

D is related to the Young’s modulus, E and the Poisson’s ratio, μ by:

11
D=Eh312(1-μ2),

where h is the plate thickness.

2 is the Laplacian in axisymmetric polar coordinates:

12
2=d2dr2+1rddr.

Kirchhoff’s theory, through satisfactory for thin plates, has been found to be unsatisfactory as the plate thickness increases. Its major limitation is the disregard for shear deformation which renders it incapable of accurately describing the behaviour of moderately thick plates; where shear deformation will be significant. Von Karman’s theory for the large deflection, axisymmetric bending of a circular plate which is either simply supported or clamped at the boundary can be represented by the system of coupled non-linear partial differential equations [23] expressed in terms of two unknowns wr and ϕr:

13
Eh312(1-μ2)2w-1rddrdϕdrdwdr=prh,
14
1E2ϕ+1rdwdrd2wdr2=0,

where:

15
2=1rddrrddr,

and E and μ are the Young’s modulus and Poisson’s ratio, respectively. h is the plate thickness, r is the radial coordinate variable. ϕ is the Airy stress function, pr is the applied radially symmetric load distribution.

The study is focused on first order shear deformable circular plate theory of Mindlin. The theory is capable of describing the behaviours of moderately thick plates where shear deformations will contribute to the flexural behaviour.

1.1. Review of solution methods for the boundary value problem of circular plates

The boundary value problem (BVP) of circular plates is represented generally by differential equations and boundary conditions. They are solved using techniques for solving ordinary and partial differential equations. The problems can also be formulated in weighted residual form or using variational calculus. Thus, the methods used in the literature for the solution of the plate problem include: Separation of variables, Integral transforms, Eigen function expansions, Fourier series, etc. The numerical methods used in the literature include: Galerkin method, Kantorovich method, Ritz method, Finite Element method and boundary element methods [32-35].

1.2. Justification for the Mindlin First order shear deformable theory

The classical two dimensional (2D) Kirchhoff – theory of flexure of elastic circular (thin) plates disregards the axial shear strains and normal stress. Hence the Kirchhoff circular plate theory cannot be expected to give good estimates of the deflection and internal forces in circular plates with a thickness to radius ratio exceeding 0.05. Mindlin [19] considered the effect of rotatory inertia and shear deformation on the flexural behaviour of isotropic circular plates and derived from first principles the governing equations for a more comprehensive 2D theory called the Mindlin first order shear deformation plate theory FSDT for circular plates.

Mindlin’s theory is analogous to the Timoshenko’s theory for beams, which takes into consideration shear deformation in the flexural analysis of beams. Mindlin’s FSDT for circular plates satisfies the stress-strain equations for transverse shear stresses and shear strains by the use of shear modification (correction) factors ks. The shear modification factors depend upon the plate material, geometry, cross-sectional shape, loading, and support/boundary conditions.

However, a fundamental inadequacy of the Mindlin’s FSDT is that all the equations of three dimensional theory of elasticity are not satisfied by the theory, and only deflections of the middle surfaces (planes)of the plates are considered [36]. Despite the obvious limitation, the Mindlin 2D plate theory plays a significant role in the validation of results for deflections and internal stress resultants of thick plate theories [36].

Reddy’s higher order shear deformation plate theory (HSDPT) assumes a quadratic (parabolic) variation of transverse shear stress through the plate thickness such that shear stress free boundary conditions required for the top and bottom surfaces of the plate are automatically satisfied. Hence the Reddy theory does not require shear correction factors.

1.3. Advantages of mathematical solutions over numerical solutions in plate bending analysis

The focus of the present study is to derive mathematical solutions, also called closed form or analytical solutions to the flexural problem of Mindlin circular plates; due to the obvious merits of the mathematical solutions over the numerical or approximate solutions. The advantages of the mathematical solutions are:

(i) The mathematical solutions aim to derive using rigorous mathematical techniques, the “exact” solutions to the governing partial differential equations of the Mindlin circular plate theory, which are valid at all points on the circular plate domain as well as on the boundaries; for given loading and restraint conditions. The solutions derived in a mathematical approach are called ‘exact’ within the limitations of the foundational theory used in formulating the plate problem, in this case, the Mindlin first order shear deformation theory. However, the numerical solutions aim to obtain approximate solutions to the plate flexure problem which may not be valid solutions at all points on the solution domain; nor at the boundaries.

(ii) The derivation of analytical solutions is important since they serve as benchmark solutions for comparing the accuracy of numerical solutions.

(iii) The advantages of mathematical solutions include its better accuracy as compared with numerical solutions and broad generality in solving plate flexure problems with given loads and restraint conditions.

(iv) The mathematical derivation of the solutions provide a systematic framework for solving the flexural problems of Mindlin circular plates with different edge restraint conditions.

Despite the advantages of mathematical solutions, exact analytical solutions for plate bending problems are very limited, even for the classical Kirchhoff-thin plate theory. They are usually mathematically rigorous, and mathematical solutions do not exist for many plate bending problems involving complicated restraint conditions and complex load distributions. For plates of arbitrary shape, the solution of the governing plate equations, even in the classical Kirchhoff thin plate theory is a non-trivial matter and in general, an exact analytical/mathematical solution does not exist [37].

Even for the classical Kirchhoff plate theory, the mathematical rigours of obtaining closed form solutions to the BVP are so complex that mathematical solutions leading to formula for stresses or deflections are only possible/available for a few simple cases of geometry, boundary conditions and loading [38].

1.4. Research aim and objectives

The general aim of this paper is to solve and obtain mathematical solutions for the boundary value problem of first order shear deformable circular plates. The general objectives are:

(i) to solve the governing boundary value problem of flexure of first order shear deformable circular plates under arbitrary distributed load pr and obtain general solutions for the generalised displacements in terms of constants of integration.

(ii) to find mathematical solutions to the flexural problem of first order shear deformable circular plates with clamped edge (r=r0) for the case of (a) uniformly distributed transverse load over the entire plate domain, and (b) point load P0 applied at the centre (r= 0) of the plate.

(iii) to find analytical solutions for the flexural problem of first order shear deformable circular plates with simply supported edges for the case of uniformly distributed transverse load over the entire plate domain.

2. Theoretical framework

Axisymmetric problems of the theory of elasticity are governed by the requirements of the strain displacement relations, the differential equations of equilibrium and the stress-strain law. The differential equations of equilibrium are given by:

16
σrrr+1rτrθθ+1r(σrr-σθθ)=0,
17
τrθr+1rτθθθ+2τrθr=0,

where σrr and σθθ are the radial and circumferential stresses and τrθ is the shear stress.

The strain-displacement relations for axisymmetric small displacement problems are given by the three equations:

18
εrr=urr,
19
εθθ=1ruθθ+urr,
20
2εrθ=γrθ=1rurθ+uθr-uθr,

where εrr is the radial strain (normal strain in the radial direction), εθθ is the circumferential normal strain, εrθ is the shear strain in the rθ direction, ur is the radial component of displacement, uθ is the circumferential component (tangential component) of the displacement.

The strain-stress relations in the axisymmetric (polar) coordinates are for plane stress conditions:

21
εrr=1Eσrr-μσθθ,
22
εθθ=1Eσθθ-μσrr,
23
εzz=-μEσrr+σθθ,
24
2εrθ=2(1+μ)Eτrθ=τrθG=γrθ,

and for plane strain conditions:

25
εrr=1+μE1-μσrr-μσθθ,
26
εθθ=1+μEμσrr+1-μσθθ,
27
γrθ=τrθG=2(1+μ)τrθE,

where E is the Young’s modulus of elasticity, μ is the Poisson’s ratio, and G is the shear modulus. G is related to the Young’s modulus E by:

28
G=E2(1+μ).

3. Governing equations

The differential equations of equilibrium of first order shear deformable circular plates are given by the system of three partial differential equations in terms of the stress resultants, Mrr, Mθθ, Mrθ, Qr and Qθ as follows:

29
-1rrrMrr+Mrθθ-Mθθ+Qr=0,
30
-1rrrMrθ+Mθθθ+Mrθ+Qθ=0,
31
-1rrrQr+θθθ-p=0.

The stress resultants for first order shear deformable circular plates are expressed in terms of the displacements w0, φr and φθ as:

32
Mrr=Dφrr+μrφr+1rφθθ,
33
Mθθ=Dμφrr+1rφr+1rφθθ,
34
Mrθ=D(1-μ)rφrθ+rφθr-φθ,
35
Qr=ks-h/2h/2τrz dz=ksGhφr+w0r,
36
Qθ=ks-h/2h/2τθz dz=ksGhφθ+1rw0θ,

where ks is the shear correction factor.

For axisymmetric problems:

37
φθ=0,
38
Mrθ=0,
39
Qθ=0.

Consequently, the differential equations of equilibrium of axially symmetric problems simplify to the system of two partial differential equations:

40
-1rrrMrr-Mθθ+Qr=0,
41
-1rddr(rQr)-p(r)=0.

For axially symmetric problems of first order shear deformable circular plates, the internal stress resultants (bending moments Mrr, Mθθ and shear force distributions Qr) are related to the generalized displacement functions φr, w0 as follows:

42
Mrr=Ddφrdr+μφrr,
43
Mθθ=Dμdφrdr+φrr,
44
Qr=ksGhφr+dw0dr.

4. General solution of the boundary value problem of first order shear deformable circular plates under transverse flexure

The differential equation of equilibrium Eq. (41) is observed to be variable separable. We thus seek a solution using the method of separation of variables. Eq. (41) is rewritten as:

45
-1rddrr Qr=pr.

Separation of variables, and integration yields:

46
d(rQr)=rQr=-r prdr+c0,

where c0 is a constant of integration.

From Eq. (40), we have:

47
rr Mrr-Mθθ=rdMrrdr+Mrr-Mθθ=rQr.

Substitution of Eqs. (42) and (43) into Eq. (47) yields after simplification:

48
Drd2φrdr2+dφrdr-φrr=rQr=Drddrdφrdr+φrr=rQr=Dr1rddr(rφr).

Using Eq. (46) we have:

49
Drddrdφrdr+φrr=-rprdr+c0.

Solution by successive integration yields:

50
Dφr=-1rr1rr prdrdrdr+c0r42logr-1+rc12+c2r,

where c1 and c2 are constants of integration.

Thus, we obtain:

51
Dφr=-I'r+r42logr-1c0+r2c1+c2r,

where:

52
I'r=1rr1rr prdrdrdr,
53
Irr=1rr1rr prdrdrdrdr.

The solution for w0 is obtained using Eqs. (44) and (46).

Thus:

54
r Qr=rksGhφr+dw0dr=-rprdr+c0,
55
ksGhdw0dr=-1rrprdr+c0r-ksGhφr,
56
ksGhdw0dr=-1rrprdr+c0r-ksGhD-I'r+c0r42logr-1+rc12+c2r.

Integrating both sides of Eq. (56) with respect to r, we obtain:

57
ksGhw0(r)=-1rrp(r)drdr
+c0rdr-ksGhD-I'r+c0r42logr-1+c1r2+c2rdr.

Thus:

58
ksGhw0(r)=-1rrp(r)drdr+c0logr
- ksGhD-Ir+ c0r24logr - 1+ r2c14+c2logr+c3,

where c3 is a fourth constant of integration.

The equations obtained for φ0r and w0r given respectively by Eqs. (51) and (58) are the general solutions for the boundary value problem of first order shear deformable circular plates subject to the transverse distributed load pr. The general solution is obtained in terms of four unknown constants of integration c0, c1, c2 and c3. These four constants of integration can be determined from the essential and natural boundary conditions for particular problems of first order shear deformable circular plates.

This paper considers three particular problems/cases, namely: (i) solid first order shear deformable circular plate with clamped edges and under uniformly distributed transverse load (ii) solid first order shear deformable circular plate with clamped edges and under point load P0 applied at the centre (iii) simply supported solid shear deformable circular plate under uniformly distributed load p0.

4.1. Case 1: First order shear deformable solid circular plate with clamped edge and under uniform transverse load distribution

For bounded solutions for the rotation φr, φr is required to be non-singular at r= 0, hence:

59
c2=0.

For solid circular plates that are not subject to point load at the centre, r= 0, the shear force at the centre is required to vanish. This yields from Eq. (46):

60
c0=0.

For clamped solid circular first order shear deformable plates of radius r0, the boundary conditions of the clamped edge r=r0 are:

61
w0r=r0=0,
62
φr(r=r0)=0.

For uniformly distributed transverse load of intensity p0, p(r)=p0:

63
Ir=1rr1rrp0drdrdrdr=p01rr1rrdrdrdrdr=p0r464,
64
I'r=p0r316.

The solution for w0r becomes:

65
ksGhw0r=-1rrp0drdr-ksGhDp0r464+r2c14+c3,
66
ks Ghw0r=-p0r24-ksGhDp0r464+r2c14+c3.

Similarly, the solution for φr becomes:

67
Dφr=-p0r316+rc12.

Using the boundary condition Eq. (62) in Eq. (67) we obtain:

68
c1=p0r028.

Application of the boundary condition Eq. (61) in Eq. (66) yields:

69
c3=p0r024+ksGhDp0r0464.

Substitution of c1 and c3 into the equation for w0r and φr give the solutions for this case as:

70
w0(r)=p0r0464D1-rr022+p0r024ksGh1-rr02
=p0r0464D1-2rr02+rr04+p0r024ksGh1-rr02,
71
φrr=p0r02r16D1-rr02,
72
φrr=p0r0316Drr01-rr02=p0r0316Drr0-rr03.

The transverse deflection w0r can be decomposed into two components, a flexural (bending) component w0f(r) and a shear component w0s(r) as follows:

73
wr=w0fr+w0sr,

where:

74
w0fr=p0r0464D1-rr022,
75
w0sr=p0r024ksGh1-rr02.

4.1.1. Maximum deflection

The maximum deflection wmax is found to occur at the centre of the plate where r= 0, and is given by:

76
wmax=p0r0464D+p0r024ksGh,
77
wmax=p0r0464D1+83ks(1-μ)hr02.

For ks= 5/6, and μ= 0.30:

78
wmax=p0r0464D1+81.75hr02.

Table 1 shows the maximum deflection of first order shear deformable solid circular plates with clamped edge for the case of uniformly distributed load computed for varying values of h/r0.

Table 2 shows the difference between wmax the maximum deflection obtained by the first order shear deformable circular plate theory and the classical Kirchhoff-Love circular plate theory for various values of h/r0 for the case of clamped edge and uniformly distributed load.

Table 1Maximum deflection of first order shear deformable circular plates with clamped edges and under uniformly distributed load p0

hr0
wmax×p0r04D
0.005
1.562679×10-2
0.01
1.563214×10-2
0.05
1.580357×10-2
0.10
1.633929×10-2
0.2
1.848214×10-2
0.3
2.220536×10-2

Table 2Difference between the maximum deflections of first order shear deformable circular plates and the Kirchhoff-Love circular plate with clamped edges (case of uniformly distributed load)

h/r0
(wmaxFSDP-wmaxCPT)
Difference %
0.005
1.142857×10-4
0.0114
0.01
4.57142×10-4
0.0457
0.05
1.142857×10-2
1.1429
0.1
4.57142×10-2
4.5714
0.2
18.2857×10-2
18.2857
0.3
41.142857×10-2
41.1429

The bending moment and shear force distributions are found by substitution of Eq. (72) into Eqs. (42) and (43) as:

79
Mrrr=p0r02161+μ-3+μrr02,
80
Mθθr=p0r02161+μ-1+3μrr02.

The bending stresses are found as:

81
σrrr,z=3p0r02z4h31+μ-3+μrr02.
82
σθθr,z=3p0r02z4h31+μ-1+3μrr02.

The shear stress distribution τrz is given by:

83
τrz(r)=Gγrz=Gφr+dw0(r)dr
=Gp0r0464D-4rr02+4r3r04+p0r024ksGh-2rr02+p0r0316Drr0 -rr03
=Gp0r0316D-rr0+rr03+rr0-rr03-p0r2ksGh=-p0r2ksh.

The shear force distribution Qrr is from Eq. (75):

84
Qrr=ksGhdw0sdr=ksGhp0r024ksGh-2rr02=-p0r2=ksh τrz=k-h/2h/2τrz dz.

The deflection and rotation functions satisfy the boundary conditions at the clamped edge since upon substitution of r=r0 into the equations, Eqs. (61) and (62) are satisfied. They also satisfy the conditions of symmetry since φr(r=0)=0.

4.2. Case 2: First order shear deformable solid circular plate with clamped edge and point load P0 applied at the centre, r= 0

For first order shear deformable solid circular plate (FSDP) subject to a point load at the centre, enforcement of boundary conditions yield:

85
c0=-P02π.

The requirement of boundedness of w0r also require that:

86
c2=0.

Then, the displacements become:

87
Dφr=r42logr-1-P02π+rc12,
88
ksGhw0(r)=-1rr p(r)drdr-P02πlogr
-ksGhD-Ir-P02πr24logr-1+r24P04π2logr0-1+c3.

Enforcement of the boundary conditions Eqs. (61) and (62) in Eqs. (87) and (88) yield the following constants of integration:

89
c1=P02π2logr0-1,
90
c3=P0logr02π+ksGhDP0r0216π.

Hence the solutions for transverse deflections and rotation functions become:

91
w0r=P0r0216πD1-rr02+2rr02logerr0-P02πksGhlogerr0=w0fr+w0sr,
92
φr=-P0r4πDlogerr0,
93
w0(r)=P0r0216πD1-rr02+2rr02logerr0-P0h32πksh6D(1-μ)logerr0
=P0r0216πD1-rr02+2rr02logerr0-P0h212πDks1-μlogerr0.

4.3. Case 3: First order shear deformable solid circular plate with simply supported edge and under uniformly distributed transverse load

The essential and natural boundary conditions of the simple supports at the edge r=r0 are:

94
w(r=r0)=0,
95
Mrr(r=r0)=0.

At the centre, r=0:

96
(r Qr)r=0=0.

From Eq. (93) it is found that:

97
c0=0.

For non singular solutions for the displacements, the expression for φr gives the condition for bounded solutions as:

98
c2=0.

Hence, for simply supported FSDP under uniformly distributed transverse load:

99
ksGhw0r=-1rrprdrdr-ksGhD-Ir+r2c14+c3,
100
Dφr=-I'r+r2logr-r4c0+r2c1+c2r,
101
Dφr=-I'r+r2c1,
102
Dφ'r=-I''r+c12,
103
φ'r=1D-I''(r)+c12,
104
Mrr=DμrD-I'(r)+r2c1+1Dc12-I''(r)=μrrc12-I'(r)+c12-I''r,
105
Mrrr=μrrc12-pr316+c12-3p0r216.

The force boundary condition yields:

106
Mrr(r=r0)=μr0r0c12-p0r0316+c12-3p0r0216=0.

Solving:

107
c1=p0r0283+μ1+μ.

The deflection boundary condition at the simply supported edge gives:

108
w(r=r0)=-p0r024+p0r04ksGh64D-ksGhr024Dp0r0283+μ1+μ+c3=0.

Solving:

109
c3=ksGhDp0r04645+μ1+μ+p0r024.

Thus, the solutions for the displacements become:

110
w0r=p0r0464D5+μ1+μ-rr04-2rr023+μ1+μ+p0r024ksGh1-rr02
=w0fr+w0sr,
111
φr=p0 r0316Drr03+μ1+μ-rr02.

4.3.1. Maximum deflection

The maximum deflection wmax occurs at the plate centre, r= 0, and is found as:

112
wmax=wr=0=p0 r0464D5+μ1+μ+p0 r024ks Gh,
113
wmax=p0 r0464D5+μ1+μ+p0 r024ks hG,
114
wmax=p0 r04D5+μ1+μ164+124(1-μ)kshr02.

For μ= 0.3, ks= 5/6:

115
wmax=p0r04D5.31.3164+124(0.7)56hr02.

4.3.2. Maximum rotation

The maximum rotation (φmax) occurs at the edge r=r0 is given by:

116
φmax=φrr=r0=p0 r038D1+μ.

4.3.3. Bending moments

The bending moment distribution Mrrr is given by:

117
Mrr=μc12-p0r216+c12-3p0r216=μc12-μp0r216+c12-3p0r216,
118
Mrrr=c121+μ-p0r2163+μ=p0 r02163+μ-p0r2163+μ,
119
Mrrr=p03+μ16r02-r2.

The bending moment at the centre of the plate is given by:

120
Mrrr=0=p0 r023+μ16.

5. Validation of the solution for simply supported circular Mindlin plate using the Ritz finite element method

The total potential energy functional Π for a Mindlin’s first order shear deformable circular plate is given by:

121
Π=12A-h/2h/2(σrrεrr+σθθεθθ+τrθγrθ+τrzγrz+τzθγzθ)dzrdrdθ-Apwrdrdθ,

where A is the cross-sectional area (domain) of the circular plate given by:

0rr0, 0θ2π.

From the principle of virtual displacements:

122
δΠ=0,
123
AMrrδφrr+Mrθ1rδφrθ+δφθr-φθr+Mθθδφrr+1rδφθθ
+ Qrδφr+δw0r+Qθδφθ+1rδw0θ-pδw0rdrdθ=0,

where:

124
Qr=ks-h/2h/2τzθdz,
125
Qr=ks-h/2h/2τzrdz.

The weak form of the virtual work statement is:

126
02π0r0Dddrφr(r)+μφr(r)rdδφr(r)dr+Dμdφr(r)dr+φr(r)rδφr(r)r
+ksGhφr(r)+dw0(r)drδφr(r)+dδw0(r)dr-p(r)δw0(r)rdrdθ=0,
127
2π0r0Ddφr(r)dr+μφr(r)rdδφrdr+Dμdφr(r)dr+φr(r)rδφr(r)r
+ksGhφr(r)+dw0(r)drδφr(r)+dδw0(r)dr-p(r)δw0(r)rdr=0,
128
0r0Ddφrdr+μφrrδdφrdr+Dμdφrdr+φrrδφrr+ksGhφr+dw0dr
×δφr+δdw0(r)dr-p(r)δw0(r)rdr=0.

Let:

129
w0r=i=1Maifir,
130
φrr=j=1Ncjgjr.

Then:

131
dφr(r)dr=j=1Ncjdgj(r)dr,
132
dw0(r)dr=i=1Maidfi(r)dr,
133
0r0D cjdgjdr+μ cj gjδ ckdgkdr
+Dμ cjdgjdr+1r cjgj(r)1rδ ckgk
+ksGh cjgj+ aidfidrδ ckgk+δ akdfkdr-p(r)δ aifirdr=0.

Equating the coefficients of δai and δci to zero, and this corresponds to the requirement that Π=0 with respect to ai and Π=0 with respect to ci we have the following Ritz variational statements:

134
j=1MksGh0r0dfi(r)drdfj(r)drrdraj+j=1MksGh0r0dfidrgj(r)rdrcj
-0r0prfirr dr=0,
135
ksGhj=1M0r0gk(r)dfj(r)drr draj+j=1NksGh0r0gk(r)gj(r)rdr
+D0r0dgk(r)drdgj(r)dr+μrgk(r)dgj(r)dr+dgkdrgj(r)+1r2gk(r)gj(r)rdrcj=0,

where k=1,2, ..., N.

In matrix form, the Ritz variational equations are:

136
XijYijYijTZijaicj=Pi0,

in which:

137
Xij=ksGh0r0dfi(r)drdfj(r)drrdr,
138
Xij=ksGh0r0fi'rfj'rrdr,
139
Yij=ksGh0r0dfirdrgjrrdr,
140
Yij=ksGh0r0fi'rgjrrdr,
141
Pi=0r0p(r)fi(r)rdr,
142
Zkj=D0r0dgk(r)drdgj(r)dr+μrgk(r)dgj(r)dr+dgk(r)drgj(r)+1r2gk(r)gj(r)rdr
+ ksGh0r0gk(r)gj(r)rdr.

5.1. Simply supported solid first order shear deformable circular plate

Shape functions for simply supported solid first order shear deformable circular plates are:

143
f1r=1-rr0,
144
f2r=1-r2r02,
145
g1r=rr0,

where f1, f2 and g1 satisfy the essential boundary conditions:

146
w0(r=r0)=0,
147
φr(r=0)=0,
148
X11=ksGh0r0f1'(r)f1'rrdr,
149
f1'r=-1r0,
150
X11=ksGh0r0-1r0-1r0 rdr=ksGhr020r0rdr=
X11=ksGhr02r220r0=ksGhr02r022=ksGh2,
151
X12=ksGh0r0f1'(r)f2'rr dr,
152
X12=ksGh0r0-1r0-2rr02rdr=2ksGhr030r0r2dr,
153
X12=2ksGhr03r330r0=2ksGhr03r033=2ksGh3,
154
Y11=ksGh0r0-1r0rr0rdr,
155
Y11=-ksGhr020r0r2dr=-ksGhr02r330r0=-ksGhr02r033=-ksGhr03,
156
Y21=ksGh f2'g1rdr,
157
Y21=ksGh0r0-2rr02rr0rdr,
158
Y21=-2ksGhr030r0r3dr=-2ksGhr03r440r0=-2ksGhr03r044=-ksGhr02,
159
X21=ksGh0r0f2'(r)f1'rdr,
160
X21=ksGh0r0-2rr02-1r0rdr=2ksGhr030r0r2dr=2ksGhr03r330r0,
161
X21=2ksGhr03r033=23ksGh,
162
X22=ksGh0r0f2'(r)f2'(r)rdr,
163
X22=ksGh0r0-2rr02-2rr02rdr,
164
X22=4ksGhr040r0r3dr=4ksGhr04r440r0=4ksGhr04r044=ksGh,
165
Z=D0r01r01r0+μrrr01r0+1r0rr0+1r2rr0rr0rdr+ksGh0r0rr0rr0rdr,
166
Z=D0r01r02+μr2rr02+1r02rdr+ksGhr020r0r3dr,
167
Z=Dr020r0(2+2μ)r dr+ksGhr020r0r3dr,
168
Z=(2+2μ)Dr020r0rdr+ksGhr020r0r3dr,
169
Z=(2+2μ)Dr02r220r0+ksGhr02r440r0,
170
Z=1+μD+ksGhr02r044=1+μD+ksGh4r02,
171
Z=ksGh123r02+12(1+μ)DksGh,
172
P1=0r0p01-rr0rdr=p00r01-rr0rdr=p00r0r-r2r0dr,
173
P1=p0r22-r33r00r0=p0r022-r033r0,
174
P1=p0r022-r023=p0r026,
175
P2=0r0p01-r2r02rdr=p00r0r-r3r02dr,
176
P2=p0r22-r44r020r0=p0r022-r024=p0r024.

The Ritz variational equations become in matrix form:

177
ksGh1268-4r0812-6r0-4r0-6r03r02+12(1+μ)DksGha1a2c1=p0r0212230,
178
68-4r0812-6r0-4r0-6r03r02+12(1+μ)DksGha1a2c1=p0r02ksGh230.

The undetermined displacement parameters are:

179
a1=0,
180
a2=p0r02244ksGh+r02D(1+μ),
181
c1=2p0r033D(1+μ).

The Ritz variational solution is:

182
w2r=p0r0424D1+μ1-r2r02+p0r026ksGh1-r2r02,
183
φ1r=2p0r033D1+μrr0.

The Ritz variational solutions are only approximate.

5.2. Validation of solution for clamped circular Mindlin plate

The solutions obtained in this study for Mindlin’s first order shear deformable solid circular plate with clamped edge and under uniform transverse load distribution are exactly the same as solutions obtained by Ike et al. [31] who used the method of direct integration of the biharmonic equation for circular Mindlin plates, and the enforcement of appropriate boundary conditions.

6. Discussion

In this work, the governing partial differential equations of equilibrium for axially symmetric flexural problems of first order shear deformable circular plates were solved analytically to obtain closed form mathematical solutions for the unknown generalised displacements w0r and φrr. The governing equations, Eqs. (40) and (41), for axially symmetric first order shear deformable (FSD) circular plate problems were first expressed in terms of the unknown generalised displacements with the aid of the differential equations relating the stress resultants and the generalised displacements Eqs. (42-44). The resulting differential expressions, now in variable separable form were integrated to obtain the general solutions for the unknown displacements as Eqs. (50) and (58). The general solution for the rotation function φrr contained three unknown constants of integration; while the general solutions for the deflection w0r contained four unknown constants of integration. Mathematical solutions for w0r and φrrwere then sought for three specific cases.

For first order shear deformable solid circular plate with clamped edge and under uniformly distributed transverse load over the entire plate domain, boundedness conditions for the rotation, and the need for the shear force at the centre to vanish were used to find two constants of integration as Eqs. (59) and (60). The boundary conditions at the clamped edge – Eqs. (61) and (62) – were used to obtain the other two constants of integration as Eqs. (68) and (69). The deflection function was thus obtained as Eq. (70) and the rotation function as Eq. (72). It was observed that the deflection function could be decomposed into a flexural component, given by Eq. (74) and a shear component given by Eq. (75). The maximum deflection was found to occur at the centre of the plate, in line with the demands of the symmetry of the plate and the loading. The maximum deflection was found as Eqs. (76) or (77). For the particular case where ks=5/6 and μ=0.30, the maximum deflection was obtained as Eq. (78). The maximum deflection of first order shear deformable circular plate with clamped edge for the case of uniformly distributed transverse load are shown in Table 1 for various values of the ratio of the plate thickness to the radius (h/r0). The difference between the maximum deflection values obtained using the first order shear deformable circular plate with clamped edge for the case of uniformly distributed load were calculated for varying ratios of h/r0 and shown in Table 2. An examination of Table 2 shows that the contribution of shear deformation to the overall bending deformation of the plate increases as the plate thickness increases. For thin plates, the difference between the maximum transverse deflection obtained using the first order shear deformable circular plate and the classical Kirchhoff-Love circular thin plate theory is insignificant. However, as the plate thickness increases, the difference becomes significant.

For the case of first order shear deformable solid circular plate with clamped edge and point load P0 applied at the centre, the mathematical solutions for w0r and φrr were obtained as Eqs. (91) and (92), respectively. The deflection w0r in the first order shear deformable circular plate was found to be decomposed into two components, flexural component, and shear component, where the flexural component corresponded to the solution for the Kirchhoff-Love circular thin plate.

It is also observed that:

184
φr=-dw0frdr.

It is observed that the shear component of the deflection is singular and undefined at r= 0, the point of application of the point load P0. It is thus mathematically impossible to determine the maximum deflection since the deflection is singular at the centre of the plate where r= 0.

For the case of first order shear deformable solid circular plate with simply supported edge and under uniformly distributed transverse load over the entire plate region, the general solutions for w0r and φrr were obtained as Eqs. (110) and (111) respectively. The deflection w0r was observed to be decomposed into two components; a flexural component and shear component. The maximum deflection was found to occur at the plate centre, in line with the symmetrical nature of the problem. The maximum deflection was obtained as Eqs. (113) or (114). The maximum deflection was also found in terms of the ratio h/r0 as Eq. (115). The maximum rotation was found to occur at the plate edge and was obtained as Equation (116). The bending moment distribution Mrrr was obtained as Eq. (119). Maximum bending moment Mrrmax was found to occur at the centre of the plate and was obtained as Eq. (120).

The significant contributions of the present study are:

(i) Closed form mathematical solutions were obtained for the general case of Mindlin’s first order shear deformable circular plate subject to an arbitrary distribution of transverse load as Eqs. (50) and (58). The solutions were derived from first principles, and were derived to satisfy the governing equations of equilibrium at all points on the solution domain, which is the two dimensional domain of the circular plate 0rr0,0θ2π.

(ii) Mathematical solutions were derived from the general solution for the particular cases of clamped edge and uniformly distributed load as Eqs. (70) and (72). The solutions were found to satisfy the boundary conditions and were derived such that all points on the circular plate domain satisfied the governing equations. For point load at the centre, the mathematical solutions were found from the general solution as Eqs. (92) and (93).

(iii) Mathematical solutions were similarly derived from the general solution for the particular case of simply supported edge and uniformly distributed load as Eqs. (110) and (111).

7. Conclusions

The following conclusions are made from the study:

(i) The partial differential equations of equilibrium of axisymmetrically loaded first order shear deformable circular plates have been solved mathematically using the method of separation of variables.

(ii) The general solutions obtained for the transverse deflection w0r and the rotation functions for any arbitrary axisymmetric load distribution contained four constants of integration, which are determined for specific cases of edge supports and transverse loads by the application of the appropriate boundary conditions.

(iii) The general solution obtained for the deflection of first order shear deformable solid circular plates fixed at the edge r=r0 was decomposable into the flexural component and a shear component. The flexural component of the deflection in this case was exactly the same as the deflection expression for a Kirchhoff-Love circular plate of the same radius. The contribution of the shear strength component to the bending deformation increased significantly with increase in the ratio of plate thickness to the radius (h/r0).

(iv) The shear force distribution for first order shear deformable solid circular plates vanished when:

185
φr=dw0(r)dr.

(v) The deflections obtained for first order shear deformable solid circular plates with clamped edge (r=r0) and subject to a point load P0 applied at the centre was decomposable into a flexural component and a shear component. The deflections in this case also increased significantly with significant increase in the ratio of the plate thickness to the radius (h/r0).

(vi) For first order shear deformable solid circular plates with simply supported edge (r=r0) and subject to uniformly distributed transverse load, the solution obtained for deflection was decomposed into a flexural component and a shear component. The deflection increased significantly with increase in the ratio of the plate thickness to the radius (h/r0). The flexural component of the bending deflection was found to be exactly identical with the deflections obtained using the Kirchhoff-Love theory for the circular plate with simply supported edge (r=r0).

(vii) For all the cases considered the maximum deflection was found to occur at the plate centre; in line with the requirements of symmetry of the problems. The expressions for the maximum deflection in all the three cases were decomposable into flexural component and shear component.

(viii) For first order shear deformable solid circular plates with simply supported edge (r=r0) the maximum rotation was found to occur at the simply supported edge (r=r0).

(ix) Mathematical (closed form) expressions were obtained for the generalised displacements w0r and φrr for any given load distribution pr and any support (restraint) condition.

(x) Mathematical (closed form) expressions were obtained for the unknown generalised displacements in the boundary value problem for specific cases of fixed edge and simply supported edge and for uniformly distributed transverse load on the plate domain and for point load applied at the centre.

References

  • Civalek O., Oztuk B. Discrete singular convolution algorithm for non-linear transient response of circular plates resting on Winkler-Pasternak foundations with different types of dynamic loading. Indian Journal of Engineering and Material Sciences, Vol. 16, 2009, p. 259-268.
  • Steele C. R., Balach C. D. Introduction to the theory of plates. Stanford University, https//web.stanford.edu/~chasst/Course%20Notes/Introduction%20to%2.
  • Timoshenko S., Woinowsky Krieger S. Theory of Plates and Shells. 2nd Edition, Tata McGraw Hill, 2010.
  • Ventsel E., Krauthammer T. Thin Plates and Shells: Theory, Analysis and Applications. Marcel Dekker Inc., USA, 2004.
  • Szilard R. Theory and Analysis of Plates: Classical and Numerical Methods. Prentice Hall, Upper Saddle River, New Jersey, 2004.
  • Chandrashekhara K. Theory of Plates. Universities Press, India, 2011.
  • Megson T. H. G. An Introduction to Aircraft Structural Analysis. Elsevier Armsterdam Butterworth – Heinemann USA, 2010.
  • Birman Victor Plate Structures: Solid Mechanics and Its Applications. Springer, New York, 2011.
  • Kirstein A. F., Woolley R. M. Symmetrical bending of thin circular elastic plates on equally spaced point support. Journal of Research of the National Bureau of Standards – C Engineering and Instrumentation, Vol. 71, Issue 1, 1967, p. 1-10.
  • Timoshenko S. P., Goodier J. M. Theory of Elasticity. Third Edition, McGraw Hill, New York, 1970.
  • Sokolnikoff I. S. Mathematical Theory of Elasticity. Second Edition, McGraw Hill, New York, 1956.
  • Hazel A. MATH 350211: Elasticity. University of Manchester, http://www.maths.manchester.ac.uk/~ahazel/MATH35021/MATH35021.html.
  • Sadd M. H. Elasticity Theory, Application and Numerics. Elsevier Academic Press, Armsterdam, 2014.
  • Barber J. R. Elasticity. Third Revised Edition, Springer Science and Business Media Derdrecht, The Netherlands, 2010.
  • Kirchhoff G. R. Uber das Gleichgewicht und die Bewegung einer Elastischen Schelbe. Journal Für Die Reine Und Angewandte Mathematik, Vol. 1850, Issue 40, 1850, p. 51-88.
  • Kapadiya H. M., Patel A. D. Review of bending solutions of thin plates. International Journal of Scientific Research and Development, Vol. 3, Issue 3, 2015, p. 1709-1712.
  • Balasubramanian Ashwin Plate Analysis with Different Geometrics and Arbitrary Boundary Conditions. M.Sc. Thesis, Mechanical Engineering Faculty of Graduate School, The University of Texas at Arlington, 2011.
  • Gujar P. S., Ladhane K. B. Bending analysis of simply supported and clamped circular plate. SSRG International Journal of Civil Engineering, Vol. 2, Issue 5, 2015, p. 69-75.
  • Mindlin R. D. Influence of rotary inertia and shear on flexural motions of isotropic elastic plates. Journal of Applied Mechanics, Vol. 18, Issue 1, 1951, p. 31-38.
  • Reissner E. On the theory of bending of elastic plates. Journal of Mathematics and Physics, Vol. 23, 1944, p. 184-191.
  • Reissner E. The effect of transverse shear deformation on the bending of elastic plates. Journal of Applied Mechanics, Vol. 12, 1945, p. 69-77.
  • Reissner E. A twelfth order theory of transverse bending of transverse isotropic plates. Journal of Applied Mathematics and Mechanics, Vol. 63, 1983, p. 285-289.
  • Dickey R. W. Nonlinear bending of circular plates. Journal of Applied Mathematics., Vol. 30, Issue 1, 1976, p. 1-9.
  • Shimpi R. P. Refined plate theory and its variants. AIAA Journal, Vol. 40, Issue 1, 2002, p. 137-146.
  • Shimpi R. P., Patel H. G. A two variable refined plate theory for orthotropic plate analysis. International Journal of Solids and Structures, Vol. 43, Issues 22-23, 2006, p. 6783-6799.
  • Levinson M. An accurate simple theory of statics and dynamics of elastic plates. Mechanics Research Communications, Vol. 7, 1980, p. 343-350.
  • Reddy J. N. A simple higher order theory for laminated composite plates. Journal of Applied Mechanics, Vol. 51, 1984, p. 745-752.
  • Reddy J. N. A refined non-linear theory of plates with transverse shear deformation. International Journal of Solids and Structures, Vol. 20, 1984, p. 881-896.
  • Suetake Y. Plate bending analysis using a modified plate theory. CMES, Vol. 11, Issue 3, 2006, p. 103-110.
  • Ike C. C. Equilibrium approach in the derivation of differential equations for homogeneous isotropic Mindlin plates. Nigerian Journal of Technology, Vol. 36, Issue 2, 2017, p. 346-350.
  • Ike C. C., Nwoji C. U., Ofondu I. O. Variational formulation of Mindlin plate equation and solution for deflections of clamped Mindlin plate. International Journal for Research in Applied Sciences and Engineering Technology, Vol. 5, Issue 1, 2017, p. 340-353.
  • Wang Yun, Xu Rongqiao, Ding Haojiang Three dimensional solution of axisymmetric bending of functionally graded circular plates. Composite Structures, Vol. 92, 2010, p. 1683-1693.
  • Mekalke G. C., Kavade M. V., Deshpande S. S. Analysis of a plate with a circular hole by FEM. Journal of Mechanical and Civil Engineering, 2012, p. 25-30.
  • Gujar P. S., Pendhari S. S. Static bending analysis of an isotropic circular plate using finite element method. International Journal of Modern Trends in Engineering and Research, Vol. 3, Issue 4, 2016, p. 491-497.
  • Aliabadi M. H., Wen P. H. Boundary Element Methods in Engineering and Sciences. Computational and Experimental Methods in Structures, Vol. 4, Imperial College Press, London, 2011.
  • Tseng W. D., Tarn J. Q., Tung L. Y., Liang W. Y. Free vibration of circular plates with various edge boundary conditions by 3D elasticity theory. Journal of Robotics and Mechanical Engineering Research, Vol. 1, Issue 2, 2015, p. 28-35.
  • Sader J. E., White L. Theoretical analysis of the static deflection of plates for atomic force microscope applications. Journal of Applied Physics, Vol. 74, Issues 1, 1993, https://doi.org/10.1063/1.354137.
  • Karamanlidis D., The Le H., Atlure S. N. Mixed finite element models for plate bending analysis. a new element and its applications. Computers and Structures, Vol. 19, Issue 4, 1984, p. 565-581.

Cited by

Mathematical modeling of the solution to the problem of plate bending given various fastenings along the contour and various loads
Serikbay Toibayev | Assemgul Uderbayeva | Zhanar Bidakhmet | Ayaulym Rakhmatulina | Algazy Zhauyt
(2020)

About this article

Received
15 March 2018
Accepted
07 May 2018
Published
30 June 2018
Keywords
first order shear deformable circular plate
shear deformation
flexural deformation
axisymmetrical problem
differential equations of equilibrium