6.6.1 Plate Equations in Polar Coordinates
6.6.1 Plate Equations in Polar Coordinates
h / 2 h / 2 h / 2
Mr z
h / 2
rr dz , M z dz,
h / 2
M r z dz
h / 2
r (6.6.1)
and
h / 2 h / 2
Vr zr dz , V z dz (6.6.2)
h / 2 h / 2
The strain-curvature relations, Eqns. 6.2.27, can be transformed to polar coordinates using
the transformations from Cartesian to polar coordinates detailed in §4.2 (in particular,
§4.2.6). One finds that {▲Problem 1}
2w
rr z
r 2
1 w 1 2 w
z 2
2
(6.6.3)
r r r
1 w 1 2w
r z 2
r r r
2w 1 w 1 2 w
M r D 2 2
2
r r r r
1 w 1 2 w 2w
M D 2
2
(6.6.4)
r r r r 2
1 w 1 2 w
M r D1 2
r r r
2
2 1 1 2 q
2 2 w (6.6.5)
r r r r 2 D
The shear forces in terms of deflection, Eqn 6.4.12, now read {▲Problem 3}
2 w 1 w 1 2 w 1 2 w 1 w 1 2 w
Vr D 2 2 , V D
2 2 2 (6.6.6)
2
r r r r r r r r r r
12 z 12 z 12 z
rr 3
M r , 3 M , r 3 M r (6.6.7)
h h h
and
3Vr z 2 3V z 2
zr 1 , z 1 (6.6.8)
2h h / 2 2h h / 2
The differential equation 6.6.5 can be solved using a method similar to the Airy stress
function method for problems in polar coordinates (the Mitchell solution), that is, a
solution is sought in the form of a Fourier series. Here, however, only axisymmetric
problems will be considered in detail.
d 2w 1 dw
M r D 2
dr r dr
1 dw d 2w
M D (6.6.9)
r dr dr 2
M r 0
2
d2 1 d 1 d d 1 d dw q(r )
2 w r r (6.6.10)
dr r dr r dr dr r dr dr D
d d 2 w 1 dw
Vr D 2 , V 0 (6.6.11)
dr dr r dr
12 z 12 z
rr 3
M r , 3 M , r 0 (6.6.12)
h h
and
3Vr z 2
zr 1 , z 0 (6.6.13)
2h h / 2
Note that there is no twisting moment, so the problem of dealing with non-zero twisting
moments on free boundaries seen with rectangular plate does not arise here.
qr 4 1 1
w A r 2 ln r 1 B r 2 C ln r D (6.6.14)
64 D 4 4
with
qr 3 1 1 1
A r 2 ln r 1 B r C
dw
dr 16 D 4 2 r
2 2
3qr 1 1 1
A 2 ln r 1 B C 2
d w
2
(6.6.15)
dr 16 D 4 2 r
3
d w 3qr 1 1 1
3
A 2C 3
dr 8D 2 r r
and
d 3 w 1 d 2 w 1 dw qr A
Vr D 3 2
2 D (6.6.16)
dr r dr r dr 2 r
There are two classes of problem to consider, plates with a central hole and plates with no
hole. For a plate with no hole in it, the condition that the stresses remain finite at the plate
centre requires that d 2 w / dr 2 remains finite, so A C 0 . Thus immediately one has
Vr qr / 2 . The boundary conditions at the outer edge r a give B and D .
d 2 1 d 1
M 0 D 2 B D1 (6.6.17)
dr r dr 2
M0
w r2 (6.6.18)
2 D1
Consider next the case of clamped plate under uniform loading. The boundary conditions
are that w dw / dr 0 at r a , leading to
qa 2 qa 4
B , D (6.6.19)
8D 64 D
and hence
w
q
64 D
r 2 a2 2
(6.6.20)
The reaction force at the outer rim is Vr (a ) qa / 2 . This is a force per unit length; the
force acting on an element of the outer rim is qa a / 2 and the total reaction force
around the outer rim is qa 2 , which balances the same applied force.
For a simply supported plate, w 0 and M r 0 at r a . Using 6.6.9a, one then has
{▲Problem 5}
3 qa 2 5 qa 4
B , D (6.6.21)
1 8D 1 64 D
and hence
q 5 2
w
64 D 1
a r2 a2 r2 (6.6.22)
The deflection for the clamped and simply supported cases are plotted in Fig. 6.6.1 (for
0.3 ).
r/a
x
0.2 0.4 0.6 0.8 1
0
-1
clamped
64 D
q -2
simply
-3
supported
-4
Consider now the case of a plate subjected to a single concentrated force F at r 0 . The
resultant shear force acting on any cylindrical portion of the plate with radius r about the
plate-centre is 2rVr (r ) . As r 0 , one must have an infinite Vr so that this resultant is
finite and equal to the applied force F. An infinite shear force implies infinite stresses. It
is possible for the stresses at the centre of the plate to be infinite. However, although the
stresses and strain might be infinite, the displacements, which are obtained from the
strains through integration, can remain, and should remain, finite. Although the solution
will be “unreal” at the plate-centre, one can again use Saint-Venant’s principle to argue
that the solution obtained will be valid everywhere except in a small region near where
the force is applied.
Thus, seek a solution which has finite displacement in which case, by symmetry, the slope
at r 0 will be zero. From the general axisymmetric solution 6.6.15a,
dw 1
C (6.6.23)
dr r 0 r r 0
so C 0 .
From 6.6.16
2rVr r 0
2A D F (6.6.24)
Thus A F / 2D and the moments and shear force become infinite at the plate-centre.
The other two constants can be obtained from the boundary conditions. For a clamped
plate, w dw / dr 0 , and one finds that {▲Problem 7}
w
F
16D
a 2 r 2 2r 2 ln r / a (6.6.25)
This solution results in ln(r / a ) terms in the expressions for moments, giving
logarithmically infinite in-plane stresses at the plate-centre.
For a plate with a hole in it, there will be four boundary conditions to determine the four
constants in Eqn. 6.6.14. For example, for a plate which is simply supported around the
outer edge r b and free on the inner surface r a , one has
M r (a ) 0, Fr (a ) 0
(6.6.26)
w(b) 0, M r (b) 0
6.6.4 Problems
1. Use the expressions 4.2.11-12, which relate second partial derivatives in the Cartesian
and polar coordinate systems, together with the strain transformation relations 4.2.17,
to derive the strain-curvature relations in polar coordinates, Eqn. 6.6.3.
2. Use the definitions of the moments, 6.6.1, and again relations 4.2.11-12, together with
the stress transformation relations 4.2.18, to derive the moment-curvature relations in
polar coordinates, Eqn. 6.6.4.
4. Use 6.2.33, 6.4.15-16 to derive the stresses in terms of moments and shear forces,
Eqns. 6.6.7-8.
5. Solve the simply supported solid plate problem and hence derive the constants 6.6.21.
6. Show that the solution for a simply supported plate (with no hole), Eqn. 6.6.22, can be
considered a superposition of the clamped solution, Eqn. 6.6.20, and a pure bending,
by taking an appropriate deflection at the plate-centre in the pure bending case.
7. Solve for the deflection in the case of a clamped solid circular plate loaded by a single
concentrated force, Eqn. 6.6.25.