In continuum mechanics , plate theories are mathematical descriptions of the mechanics of flat plates that draw on the theory of beams . Plates are defined as plane structural elements with a small thickness compared to the planar dimensions.[ 1] The typical thickness to width ratio of a plate structure is less than 0.1.[citation needed ] A plate theory takes advantage of this disparity in length scale to reduce the full three-dimensional solid mechanics problem to a two-dimensional problem. The aim of plate theory is to calculate the deformation and stresses in a plate subjected to loads .
Vibration mode of a clamped square plate
Of the numerous plate theories that have been developed since the late 19th century, two are widely accepted and used in engineering. These are
the Kirchhoff –Love theory of plates (classical plate theory)
The Uflyand-Mindlin theory of plates (first-order shear plate theory)
Kirchhoff–Love theory for thin plates
edit
Deformation of a thin plate highlighting the displacement, the mid-surface (red) and the normal to the mid-surface (blue)
The Kirchhoff –Love theory is an extension of Euler–Bernoulli beam theory to thin plates. The theory was developed in 1888 by Love[ 2] using assumptions proposed by Kirchhoff. It is assumed that a mid-surface plane can be used to represent the three-dimensional plate in two-dimensional form.
The following kinematic assumptions are made in this theory:[ 3]
straight lines normal to the mid-surface remain straight after deformation
straight lines normal to the mid-surface remain normal to the mid-surface after deformation
the thickness of the plate does not change during a deformation.
The Kirchhoff hypothesis implies that the displacement field has the form
u
α
(
x
)
=
u
α
0
(
x
1
,
x
2
)
−
x
3
∂
w
0
∂
x
α
=
u
α
0
−
x
3
w
,
α
0
;
α
=
1
,
2
u
3
(
x
)
=
w
0
(
x
1
,
x
2
)
{\displaystyle {\begin{aligned}u_{\alpha }(\mathbf {x} )&=u_{\alpha }^{0}(x_{1},x_{2})-x_{3}~{\frac {\partial w^{0}}{\partial x_{\alpha }}}=u_{\alpha }^{0}-x_{3}~w_{,\alpha }^{0}~;~~\alpha =1,2\\u_{3}(\mathbf {x} )&=w^{0}(x_{1},x_{2})\end{aligned}}}
where
x
1
{\displaystyle x_{1}}
and
x
2
{\displaystyle x_{2}}
are the Cartesian coordinates on the mid-surface of the undeformed plate,
x
3
{\displaystyle x_{3}}
is the coordinate for the thickness direction,
u
1
0
,
u
2
0
{\displaystyle u_{1}^{0},u_{2}^{0}}
are the in-plane displacements of the mid-surface, and
w
0
{\displaystyle w^{0}}
is the displacement of the mid-surface in the
x
3
{\displaystyle x_{3}}
direction.
If
φ
α
{\displaystyle \varphi _{\alpha }}
are the angles of rotation of the normal to the mid-surface, then in the Kirchhoff–Love theory
φ
α
=
w
,
α
0
.
{\displaystyle \varphi _{\alpha }=w_{,\alpha }^{0}\,.}
Displacement of the mid-surface (left) and of a normal (right)
Strain-displacement relations
edit
For the situation where the strains in the plate are infinitesimal and the rotations of the mid-surface normals are less than 10° the strains-displacement relations are
ε
α
β
=
1
2
(
u
α
,
β
0
+
u
β
,
α
0
)
−
x
3
w
,
α
β
0
ε
α
3
=
−
w
,
α
0
+
w
,
α
0
=
0
ε
33
=
0
{\displaystyle {\begin{aligned}\varepsilon _{\alpha \beta }&={\tfrac {1}{2}}(u_{\alpha ,\beta }^{0}+u_{\beta ,\alpha }^{0})-x_{3}~w_{,\alpha \beta }^{0}\\\varepsilon _{\alpha 3}&=-w_{,\alpha }^{0}+w_{,\alpha }^{0}=0\\\varepsilon _{33}&=0\end{aligned}}}
Therefore, the only non-zero strains are in the in-plane directions.
If the rotations of the normals to the mid-surface are in the range of 10° to 15°, the strain-displacement relations can be approximated using the von Kármán strains. Then the kinematic assumptions of Kirchhoff-Love theory lead to the following strain-displacement relations
ε
α
β
=
1
2
(
u
α
,
β
0
+
u
β
,
α
0
+
w
,
α
0
w
,
β
0
)
−
x
3
w
,
α
β
0
ε
α
3
=
−
w
,
α
0
+
w
,
α
0
=
0
ε
33
=
0
{\displaystyle {\begin{aligned}\varepsilon _{\alpha \beta }&={\frac {1}{2}}(u_{\alpha ,\beta }^{0}+u_{\beta ,\alpha }^{0}+w_{,\alpha }^{0}~w_{,\beta }^{0})-x_{3}~w_{,\alpha \beta }^{0}\\\varepsilon _{\alpha 3}&=-w_{,\alpha }^{0}+w_{,\alpha }^{0}=0\\\varepsilon _{33}&=0\end{aligned}}}
This theory is nonlinear because of the quadratic terms in the strain-displacement relations.
Equilibrium equations
edit
The equilibrium equations for the plate can be derived from the principle of virtual work . For the situation where the strains and rotations of the plate are small, the equilibrium equations for an unloaded plate are given by
N
α
β
,
α
=
0
M
α
β
,
α
β
=
0
{\displaystyle {\begin{aligned}N_{\alpha \beta ,\alpha }&=0\\M_{\alpha \beta ,\alpha \beta }&=0\end{aligned}}}
where the stress resultants and stress moment resultants are defined as
N
α
β
:=
∫
−
h
h
σ
α
β
d
x
3
;
M
α
β
:=
∫
−
h
h
x
3
σ
α
β
d
x
3
{\displaystyle N_{\alpha \beta }:=\int _{-h}^{h}\sigma _{\alpha \beta }~dx_{3}~;~~M_{\alpha \beta }:=\int _{-h}^{h}x_{3}~\sigma _{\alpha \beta }~dx_{3}}
and the thickness of the plate is
2
h
{\displaystyle 2h}
. The quantities
σ
α
β
{\displaystyle \sigma _{\alpha \beta }}
are the stresses.
If the plate is loaded by an external distributed load
q
(
x
)
{\displaystyle q(x)}
that is normal to the mid-surface and directed in the positive
x
3
{\displaystyle x_{3}}
direction, the principle of virtual work then leads to the equilibrium equations
N
α
β
,
α
=
0
M
α
β
,
α
β
−
q
=
0
{\displaystyle {\begin{aligned}N_{\alpha \beta ,\alpha }&=0\\M_{\alpha \beta ,\alpha \beta }-q&=0\end{aligned}}}
For moderate rotations, the strain-displacement relations take the von Karman form and the equilibrium equations can be expressed as
N
α
β
,
α
=
0
M
α
β
,
α
β
+
[
N
α
β
w
,
β
0
]
,
α
−
q
=
0
{\displaystyle {\begin{aligned}N_{\alpha \beta ,\alpha }&=0\\M_{\alpha \beta ,\alpha \beta }+[N_{\alpha \beta }~w_{,\beta }^{0}]_{,\alpha }-q&=0\end{aligned}}}
The boundary conditions that are needed to solve the equilibrium equations of plate theory can be obtained from the boundary terms in the principle of virtual work.
For small strains and small rotations, the boundary conditions are
n
α
N
α
β
o
r
u
β
0
n
α
M
α
β
,
β
o
r
w
0
n
β
M
α
β
o
r
w
,
α
0
{\displaystyle {\begin{aligned}n_{\alpha }~N_{\alpha \beta }&\quad \mathrm {or} \quad u_{\beta }^{0}\\n_{\alpha }~M_{\alpha \beta ,\beta }&\quad \mathrm {or} \quad w^{0}\\n_{\beta }~M_{\alpha \beta }&\quad \mathrm {or} \quad w_{,\alpha }^{0}\end{aligned}}}
Note that the quantity
n
α
M
α
β
,
β
{\displaystyle n_{\alpha }~M_{\alpha \beta ,\beta }}
is an effective shear force.
Stress–strain relations
edit
The stress–strain relations for a linear elastic Kirchhoff plate are given by
[
σ
11
σ
22
σ
12
]
=
[
C
11
C
12
C
13
C
12
C
22
C
23
C
13
C
23
C
33
]
[
ε
11
ε
22
ε
12
]
{\displaystyle {\begin{bmatrix}\sigma _{11}\\\sigma _{22}\\\sigma _{12}\end{bmatrix}}={\begin{bmatrix}C_{11}&C_{12}&C_{13}\\C_{12}&C_{22}&C_{23}\\C_{13}&C_{23}&C_{33}\end{bmatrix}}{\begin{bmatrix}\varepsilon _{11}\\\varepsilon _{22}\\\varepsilon _{12}\end{bmatrix}}}
Since
σ
α
3
{\displaystyle \sigma _{\alpha 3}}
and
σ
33
{\displaystyle \sigma _{33}}
do not appear in the equilibrium equations it is implicitly assumed that these quantities do not have any effect on the momentum balance and are neglected.
It is more convenient to work with the stress and moment resultants that enter the equilibrium equations. These are related to the displacements by
[
N
11
N
22
N
12
]
=
{
∫
−
h
h
[
C
11
C
12
C
13
C
12
C
22
C
23
C
13
C
23
C
33
]
d
x
3
}
[
u
1
,
1
0
u
2
,
2
0
1
2
(
u
1
,
2
0
+
u
2
,
1
0
)
]
{\displaystyle {\begin{bmatrix}N_{11}\\N_{22}\\N_{12}\end{bmatrix}}=\left\{\int _{-h}^{h}{\begin{bmatrix}C_{11}&C_{12}&C_{13}\\C_{12}&C_{22}&C_{23}\\C_{13}&C_{23}&C_{33}\end{bmatrix}}~dx_{3}\right\}{\begin{bmatrix}u_{1,1}^{0}\\u_{2,2}^{0}\\{\frac {1}{2}}~(u_{1,2}^{0}+u_{2,1}^{0})\end{bmatrix}}}
and
[
M
11
M
22
M
12
]
=
−
{
∫
−
h
h
x
3
2
[
C
11
C
12
C
13
C
12
C
22
C
23
C
13
C
23
C
33
]
d
x
3
}
[
w
,
11
0
w
,
22
0
w
,
12
0
]
.
{\displaystyle {\begin{bmatrix}M_{11}\\M_{22}\\M_{12}\end{bmatrix}}=-\left\{\int _{-h}^{h}x_{3}^{2}~{\begin{bmatrix}C_{11}&C_{12}&C_{13}\\C_{12}&C_{22}&C_{23}\\C_{13}&C_{23}&C_{33}\end{bmatrix}}~dx_{3}\right\}{\begin{bmatrix}w_{,11}^{0}\\w_{,22}^{0}\\w_{,12}^{0}\end{bmatrix}}\,.}
The extensional stiffnesses are the quantities
A
α
β
:=
∫
−
h
h
C
α
β
d
x
3
{\displaystyle A_{\alpha \beta }:=\int _{-h}^{h}C_{\alpha \beta }~dx_{3}}
The bending stiffnesses (also called flexural rigidity ) are the quantities
D
α
β
:=
∫
−
h
h
x
3
2
C
α
β
d
x
3
{\displaystyle D_{\alpha \beta }:=\int _{-h}^{h}x_{3}^{2}~C_{\alpha \beta }~dx_{3}}
Isotropic and homogeneous Kirchhoff plate
edit
For an isotropic and homogeneous plate, the stress–strain relations are
[
σ
11
σ
22
σ
12
]
=
E
1
−
ν
2
[
1
ν
0
ν
1
0
0
0
1
−
ν
]
[
ε
11
ε
22
ε
12
]
.
{\displaystyle {\begin{bmatrix}\sigma _{11}\\\sigma _{22}\\\sigma _{12}\end{bmatrix}}={\cfrac {E}{1-\nu ^{2}}}{\begin{bmatrix}1&\nu &0\\\nu &1&0\\0&0&1-\nu \end{bmatrix}}{\begin{bmatrix}\varepsilon _{11}\\\varepsilon _{22}\\\varepsilon _{12}\end{bmatrix}}\,.}
The moments corresponding to these stresses are
[
M
11
M
22
M
12
]
=
−
2
h
3
E
3
(
1
−
ν
2
)
[
1
ν
0
ν
1
0
0
0
1
−
ν
]
[
w
,
11
0
w
,
22
0
w
,
12
0
]
{\displaystyle {\begin{bmatrix}M_{11}\\M_{22}\\M_{12}\end{bmatrix}}=-{\cfrac {2h^{3}E}{3(1-\nu ^{2})}}~{\begin{bmatrix}1&\nu &0\\\nu &1&0\\0&0&1-\nu \end{bmatrix}}{\begin{bmatrix}w_{,11}^{0}\\w_{,22}^{0}\\w_{,12}^{0}\end{bmatrix}}}
The displacements
u
1
0
{\displaystyle u_{1}^{0}}
and
u
2
0
{\displaystyle u_{2}^{0}}
are zero under pure bending conditions. For an isotropic, homogeneous plate under pure bending the governing equation is
∂
4
w
∂
x
1
4
+
2
∂
4
w
∂
x
1
2
∂
x
2
2
+
∂
4
w
∂
x
2
4
=
0
where
w
:=
w
0
.
{\displaystyle {\frac {\partial ^{4}w}{\partial x_{1}^{4}}}+2{\frac {\partial ^{4}w}{\partial x_{1}^{2}\partial x_{2}^{2}}}+{\frac {\partial ^{4}w}{\partial x_{2}^{4}}}=0\quad {\text{where}}\quad w:=w^{0}\,.}
In index notation,
w
,
1111
0
+
2
w
,
1212
0
+
w
,
2222
0
=
0
.
{\displaystyle w_{,1111}^{0}+2~w_{,1212}^{0}+w_{,2222}^{0}=0\,.}
In direct tensor notation, the governing equation is
∇
2
∇
2
w
=
0
.
{\displaystyle \nabla ^{2}\nabla ^{2}w=0\,.}
For a transversely loaded plate without axial deformations, the governing equation has the form
∂
4
w
∂
x
1
4
+
2
∂
4
w
∂
x
1
2
∂
x
2
2
+
∂
4
w
∂
x
2
4
=
−
q
D
{\displaystyle {\frac {\partial ^{4}w}{\partial x_{1}^{4}}}+2{\frac {\partial ^{4}w}{\partial x_{1}^{2}\partial x_{2}^{2}}}+{\frac {\partial ^{4}w}{\partial x_{2}^{4}}}=-{\frac {q}{D}}}
where
D
:=
2
h
3
E
3
(
1
−
ν
2
)
.
{\displaystyle D:={\cfrac {2h^{3}E}{3(1-\nu ^{2})}}\,.}
for a plate with thickness
2
h
{\displaystyle 2h}
.
In index notation,
w
,
1111
0
+
2
w
,
1212
0
+
w
,
2222
0
=
−
q
D
{\displaystyle w_{,1111}^{0}+2\,w_{,1212}^{0}+w_{,2222}^{0}=-{\frac {q}{D}}}
and in direct notation
∇
2
∇
2
w
=
−
q
D
.
{\displaystyle \nabla ^{2}\nabla ^{2}w=-{\frac {q}{D}}\,.}
In cylindrical coordinates
(
r
,
θ
,
z
)
{\displaystyle (r,\theta ,z)}
, the governing equation is
1
r
d
d
r
[
r
d
d
r
{
1
r
d
d
r
(
r
d
w
d
r
)
}
]
=
−
q
D
.
{\displaystyle {\frac {1}{r}}{\cfrac {d}{dr}}\left[r{\cfrac {d}{dr}}\left\{{\frac {1}{r}}{\cfrac {d}{dr}}\left(r{\cfrac {dw}{dr}}\right)\right\}\right]=-{\frac {q}{D}}\,.}
Orthotropic and homogeneous Kirchhoff plate
edit
For an orthotropic plate
[
C
11
C
12
C
13
C
12
C
22
C
23
C
13
C
23
C
33
]
=
1
1
−
ν
12
ν
21
[
E
1
ν
12
E
2
0
ν
21
E
1
E
2
0
0
0
2
G
12
(
1
−
ν
12
ν
21
)
]
.
{\displaystyle {\begin{bmatrix}C_{11}&C_{12}&C_{13}\\C_{12}&C_{22}&C_{23}\\C_{13}&C_{23}&C_{33}\end{bmatrix}}={\cfrac {1}{1-\nu _{12}\nu _{21}}}{\begin{bmatrix}E_{1}&\nu _{12}E_{2}&0\\\nu _{21}E_{1}&E_{2}&0\\0&0&2G_{12}(1-\nu _{12}\nu _{21})\end{bmatrix}}\,.}
Therefore,
[
A
11
A
12
A
13
A
21
A
22
A
23
A
31
A
32
A
33
]
=
2
h
1
−
ν
12
ν
21
[
E
1
ν
12
E
2
0
ν
21
E
1
E
2
0
0
0
2
G
12
(
1
−
ν
12
ν
21
)
]
{\displaystyle {\begin{bmatrix}A_{11}&A_{12}&A_{13}\\A_{21}&A_{22}&A_{23}\\A_{31}&A_{32}&A_{33}\end{bmatrix}}={\cfrac {2h}{1-\nu _{12}\nu _{21}}}{\begin{bmatrix}E_{1}&\nu _{12}E_{2}&0\\\nu _{21}E_{1}&E_{2}&0\\0&0&2G_{12}(1-\nu _{12}\nu _{21})\end{bmatrix}}}
and
[
D
11
D
12
D
13
D
21
D
22
D
23
D
31
D
32
D
33
]
=
2
h
3
3
(
1
−
ν
12
ν
21
)
[
E
1
ν
12
E
2
0
ν
21
E
1
E
2
0
0
0
2
G
12
(
1
−
ν
12
ν
21
)
]
.
{\displaystyle {\begin{bmatrix}D_{11}&D_{12}&D_{13}\\D_{21}&D_{22}&D_{23}\\D_{31}&D_{32}&D_{33}\end{bmatrix}}={\cfrac {2h^{3}}{3(1-\nu _{12}\nu _{21})}}{\begin{bmatrix}E_{1}&\nu _{12}E_{2}&0\\\nu _{21}E_{1}&E_{2}&0\\0&0&2G_{12}(1-\nu _{12}\nu _{21})\end{bmatrix}}\,.}
The governing equation of an orthotropic Kirchhoff plate loaded transversely by a distributed load
q
{\displaystyle q}
per unit area is
D
x
w
,
1111
0
+
2
D
x
y
w
,
1122
0
+
D
y
w
,
2222
0
=
−
q
{\displaystyle D_{x}w_{,1111}^{0}+2D_{xy}w_{,1122}^{0}+D_{y}w_{,2222}^{0}=-q}
where
D
x
=
D
11
=
2
h
3
E
1
3
(
1
−
ν
12
ν
21
)
D
y
=
D
22
=
2
h
3
E
2
3
(
1
−
ν
12
ν
21
)
D
x
y
=
D
33
+
1
2
(
ν
21
D
11
+
ν
12
D
22
)
=
D
33
+
ν
21
D
11
=
4
h
3
G
12
3
+
2
h
3
ν
21
E
1
3
(
1
−
ν
12
ν
21
)
.
{\displaystyle {\begin{aligned}D_{x}&=D_{11}={\frac {2h^{3}E_{1}}{3(1-\nu _{12}\nu _{21})}}\\D_{y}&=D_{22}={\frac {2h^{3}E_{2}}{3(1-\nu _{12}\nu _{21})}}\\D_{xy}&=D_{33}+{\tfrac {1}{2}}(\nu _{21}D_{11}+\nu _{12}D_{22})=D_{33}+\nu _{21}D_{11}={\frac {4h^{3}G_{12}}{3}}+{\frac {2h^{3}\nu _{21}E_{1}}{3(1-\nu _{12}\nu _{21})}}\,.\end{aligned}}}
Dynamics of thin Kirchhoff plates
edit
The dynamic theory of plates determines the propagation of waves in the plates, and the study of standing waves and vibration modes.
The governing equations for the dynamics of a Kirchhoff–Love plate are
N
α
β
,
β
=
J
1
u
¨
α
0
M
α
β
,
α
β
−
q
(
x
,
t
)
=
J
1
w
¨
0
−
J
3
w
¨
,
α
α
0
{\displaystyle {\begin{aligned}N_{\alpha \beta ,\beta }&=J_{1}~{\ddot {u}}_{\alpha }^{0}\\M_{\alpha \beta ,\alpha \beta }-q(x,t)&=J_{1}~{\ddot {w}}^{0}-J_{3}~{\ddot {w}}_{,\alpha \alpha }^{0}\end{aligned}}}
where, for a plate with density
ρ
=
ρ
(
x
)
{\displaystyle \rho =\rho (x)}
,
J
1
:=
∫
−
h
h
ρ
d
x
3
=
2
ρ
h
;
J
3
:=
∫
−
h
h
x
3
2
ρ
d
x
3
=
2
3
ρ
h
3
{\displaystyle J_{1}:=\int _{-h}^{h}\rho ~dx_{3}=2~\rho ~h~;~~J_{3}:=\int _{-h}^{h}x_{3}^{2}~\rho ~dx_{3}={\frac {2}{3}}~\rho ~h^{3}}
and
u
˙
i
=
∂
u
i
∂
t
;
u
¨
i
=
∂
2
u
i
∂
t
2
;
u
i
,
α
=
∂
u
i
∂
x
α
;
u
i
,
α
β
=
∂
2
u
i
∂
x
α
∂
x
β
{\displaystyle {\dot {u}}_{i}={\frac {\partial u_{i}}{\partial t}}~;~~{\ddot {u}}_{i}={\frac {\partial ^{2}u_{i}}{\partial t^{2}}}~;~~u_{i,\alpha }={\frac {\partial u_{i}}{\partial x_{\alpha }}}~;~~u_{i,\alpha \beta }={\frac {\partial ^{2}u_{i}}{\partial x_{\alpha }\partial x_{\beta }}}}
The figures below show some vibrational modes of a circular plate.
mode k = 0, p = 1
mode k = 1, p = 2
The governing equations simplify considerably for isotropic and homogeneous plates for which the in-plane deformations can be neglected and have the form
D
(
∂
4
w
0
∂
x
1
4
+
2
∂
4
w
0
∂
x
1
2
∂
x
2
2
+
∂
4
w
0
∂
x
2
4
)
=
−
q
(
x
1
,
x
2
,
t
)
−
2
ρ
h
∂
2
w
0
∂
t
2
.
{\displaystyle D\,\left({\frac {\partial ^{4}w^{0}}{\partial x_{1}^{4}}}+2{\frac {\partial ^{4}w^{0}}{\partial x_{1}^{2}\partial x_{2}^{2}}}+{\frac {\partial ^{4}w^{0}}{\partial x_{2}^{4}}}\right)=-q(x_{1},x_{2},t)-2\rho h\,{\frac {\partial ^{2}w^{0}}{\partial t^{2}}}\,.}
where
D
{\displaystyle D}
is the bending stiffness of the plate. For a uniform plate of thickness
2
h
{\displaystyle 2h}
,
D
:=
2
h
3
E
3
(
1
−
ν
2
)
.
{\displaystyle D:={\cfrac {2h^{3}E}{3(1-\nu ^{2})}}\,.}
In direct notation
D
∇
2
∇
2
w
0
=
−
q
(
x
,
y
,
t
)
−
2
ρ
h
w
¨
0
.
{\displaystyle D\,\nabla ^{2}\nabla ^{2}w^{0}=-q(x,y,t)-2\rho h\,{\ddot {w}}^{0}\,.}
Uflyand-Mindlin theory for thick plates
edit
In the theory of thick plates, or theory of Yakov S. Uflyand[ 4] (see, for details, Elishakoff 's handbook[ 5] ), Raymond Mindlin [ 6] and Eric Reissner , the normal to the mid-surface remains straight but not necessarily perpendicular to the mid-surface. If
φ
1
{\displaystyle \varphi _{1}}
and
φ
2
{\displaystyle \varphi _{2}}
designate the angles which the mid-surface makes with the
x
3
{\displaystyle x_{3}}
axis then
φ
1
≠
w
,
1
;
φ
2
≠
w
,
2
{\displaystyle \varphi _{1}\neq w_{,1}~;~~\varphi _{2}\neq w_{,2}}
Then the Mindlin–Reissner hypothesis implies that
u
α
(
x
)
=
u
α
0
(
x
1
,
x
2
)
−
x
3
φ
α
;
α
=
1
,
2
u
3
(
x
)
=
w
0
(
x
1
,
x
2
)
{\displaystyle {\begin{aligned}u_{\alpha }(\mathbf {x} )&=u_{\alpha }^{0}(x_{1},x_{2})-x_{3}~\varphi _{\alpha }~;~~\alpha =1,2\\u_{3}(\mathbf {x} )&=w^{0}(x_{1},x_{2})\end{aligned}}}
Strain-displacement relations
edit
Depending on the amount of rotation of the plate normals two different approximations for the strains can be derived from the basic kinematic assumptions.
For small strains and small rotations the strain-displacement relations for Mindlin–Reissner plates are
ε
α
β
=
1
2
(
u
α
,
β
0
+
u
β
,
α
0
)
−
x
3
2
(
φ
α
,
β
+
φ
β
,
α
)
ε
α
3
=
1
2
(
w
,
α
0
−
φ
α
)
ε
33
=
0
{\displaystyle {\begin{aligned}\varepsilon _{\alpha \beta }&={\frac {1}{2}}(u_{\alpha ,\beta }^{0}+u_{\beta ,\alpha }^{0})-{\frac {x_{3}}{2}}~(\varphi _{\alpha ,\beta }+\varphi _{\beta ,\alpha })\\\varepsilon _{\alpha 3}&={\cfrac {1}{2}}\left(w_{,\alpha }^{0}-\varphi _{\alpha }\right)\\\varepsilon _{33}&=0\end{aligned}}}
The shear strain, and hence the shear stress , across the thickness of the plate is not neglected in this theory. However, the shear strain is constant across the thickness of the plate. This cannot be accurate since the shear stress is known to be parabolic even for simple plate geometries. To account for the inaccuracy in the shear strain, a shear correction factor (
κ
{\displaystyle \kappa }
) is applied so that the correct amount of internal energy is predicted by the theory. Then
ε
α
3
=
1
2
κ
(
w
,
α
0
−
φ
α
)
{\displaystyle \varepsilon _{\alpha 3}={\cfrac {1}{2}}~\kappa ~\left(w_{,\alpha }^{0}-\varphi _{\alpha }\right)}
Equilibrium equations
edit
The equilibrium equations have slightly different forms depending on the amount of bending expected in the plate. For the situation where the strains and rotations of the plate are small the equilibrium equations for a Mindlin–Reissner plate are
N
α
β
,
α
=
0
M
α
β
,
β
−
Q
α
=
0
Q
α
,
α
+
q
=
0
.
{\displaystyle {\begin{aligned}&N_{\alpha \beta ,\alpha }=0\\&M_{\alpha \beta ,\beta }-Q_{\alpha }=0\\&Q_{\alpha ,\alpha }+q=0\,.\end{aligned}}}
The resultant shear forces in the above equations are defined as
Q
α
:=
κ
∫
−
h
h
σ
α
3
d
x
3
.
{\displaystyle Q_{\alpha }:=\kappa ~\int _{-h}^{h}\sigma _{\alpha 3}~dx_{3}\,.}
The boundary conditions are indicated by the boundary terms in the principle of virtual work.
If the only external force is a vertical force on the top surface of the plate, the boundary conditions are
n
α
N
α
β
o
r
u
β
0
n
α
M
α
β
o
r
φ
α
n
α
Q
α
o
r
w
0
{\displaystyle {\begin{aligned}n_{\alpha }~N_{\alpha \beta }&\quad \mathrm {or} \quad u_{\beta }^{0}\\n_{\alpha }~M_{\alpha \beta }&\quad \mathrm {or} \quad \varphi _{\alpha }\\n_{\alpha }~Q_{\alpha }&\quad \mathrm {or} \quad w^{0}\end{aligned}}}
Constitutive relations
edit
The stress–strain relations for a linear elastic Mindlin–Reissner plate are given by
σ
α
β
=
C
α
β
γ
θ
ε
γ
θ
σ
α
3
=
C
α
3
γ
θ
ε
γ
θ
σ
33
=
C
33
γ
θ
ε
γ
θ
{\displaystyle {\begin{aligned}\sigma _{\alpha \beta }&=C_{\alpha \beta \gamma \theta }~\varepsilon _{\gamma \theta }\\\sigma _{\alpha 3}&=C_{\alpha 3\gamma \theta }~\varepsilon _{\gamma \theta }\\\sigma _{33}&=C_{33\gamma \theta }~\varepsilon _{\gamma \theta }\end{aligned}}}
Since
σ
33
{\displaystyle \sigma _{33}}
does not appear in the equilibrium equations it is implicitly assumed that it do not have any effect on the momentum balance and is neglected. This assumption is also called the plane stress assumption. The remaining stress–strain relations for an orthotropic material , in matrix form, can be written as
[
σ
11
σ
22
σ
23
σ
31
σ
12
]
=
[
C
11
C
12
0
0
0
C
12
C
22
0
0
0
0
0
C
44
0
0
0
0
0
C
55
0
0
0
0
0
C
66
]
[
ε
11
ε
22
ε
23
ε
31
ε
12
]
{\displaystyle {\begin{bmatrix}\sigma _{11}\\\sigma _{22}\\\sigma _{23}\\\sigma _{31}\\\sigma _{12}\end{bmatrix}}={\begin{bmatrix}C_{11}&C_{12}&0&0&0\\C_{12}&C_{22}&0&0&0\\0&0&C_{44}&0&0\\0&0&0&C_{55}&0\\0&0&0&0&C_{66}\end{bmatrix}}{\begin{bmatrix}\varepsilon _{11}\\\varepsilon _{22}\\\varepsilon _{23}\\\varepsilon _{31}\\\varepsilon _{12}\end{bmatrix}}}
Then,
[
N
11
N
22
N
12
]
=
{
∫
−
h
h
[
C
11
C
12
0
C
12
C
22
0
0
0
C
66
]
d
x
3
}
[
u
1
,
1
0
u
2
,
2
0
1
2
(
u
1
,
2
0
+
u
2
,
1
0
)
]
{\displaystyle {\begin{bmatrix}N_{11}\\N_{22}\\N_{12}\end{bmatrix}}=\left\{\int _{-h}^{h}{\begin{bmatrix}C_{11}&C_{12}&0\\C_{12}&C_{22}&0\\0&0&C_{66}\end{bmatrix}}~dx_{3}\right\}{\begin{bmatrix}u_{1,1}^{0}\\u_{2,2}^{0}\\{\frac {1}{2}}~(u_{1,2}^{0}+u_{2,1}^{0})\end{bmatrix}}}
and
[
M
11
M
22
M
12
]
=
−
{
∫
−
h
h
x
3
2
[
C
11
C
12
0
C
12
C
22
0
0
0
C
66
]
d
x
3
}
[
φ
1
,
1
φ
2
,
2
1
2
(
φ
1
,
2
+
φ
2
,
1
)
]
{\displaystyle {\begin{bmatrix}M_{11}\\M_{22}\\M_{12}\end{bmatrix}}=-\left\{\int _{-h}^{h}x_{3}^{2}~{\begin{bmatrix}C_{11}&C_{12}&0\\C_{12}&C_{22}&0\\0&0&C_{66}\end{bmatrix}}~dx_{3}\right\}{\begin{bmatrix}\varphi _{1,1}\\\varphi _{2,2}\\{\frac {1}{2}}~(\varphi _{1,2}+\varphi _{2,1})\end{bmatrix}}}
For the shear terms
[
Q
1
Q
2
]
=
κ
2
{
∫
−
h
h
[
C
55
0
0
C
44
]
d
x
3
}
[
w
,
1
0
−
φ
1
w
,
2
0
−
φ
2
]
{\displaystyle {\begin{bmatrix}Q_{1}\\Q_{2}\end{bmatrix}}={\cfrac {\kappa }{2}}\left\{\int _{-h}^{h}{\begin{bmatrix}C_{55}&0\\0&C_{44}\end{bmatrix}}~dx_{3}\right\}{\begin{bmatrix}w_{,1}^{0}-\varphi _{1}\\w_{,2}^{0}-\varphi _{2}\end{bmatrix}}}
The extensional stiffnesses are the quantities
A
α
β
:=
∫
−
h
h
C
α
β
d
x
3
{\displaystyle A_{\alpha \beta }:=\int _{-h}^{h}C_{\alpha \beta }~dx_{3}}
The bending stiffnesses are the quantities
D
α
β
:=
∫
−
h
h
x
3
2
C
α
β
d
x
3
{\displaystyle D_{\alpha \beta }:=\int _{-h}^{h}x_{3}^{2}~C_{\alpha \beta }~dx_{3}}
Isotropic and homogeneous Uflyand-Mindlin plates
edit
For uniformly thick, homogeneous, and isotropic plates, the stress–strain relations in the plane of the plate are
[
σ
11
σ
22
σ
12
]
=
E
1
−
ν
2
[
1
ν
0
ν
1
0
0
0
1
−
ν
]
[
ε
11
ε
22
ε
12
]
.
{\displaystyle {\begin{bmatrix}\sigma _{11}\\\sigma _{22}\\\sigma _{12}\end{bmatrix}}={\cfrac {E}{1-\nu ^{2}}}{\begin{bmatrix}1&\nu &0\\\nu &1&0\\0&0&1-\nu \end{bmatrix}}{\begin{bmatrix}\varepsilon _{11}\\\varepsilon _{22}\\\varepsilon _{12}\end{bmatrix}}\,.}
where
E
{\displaystyle E}
is the Young's modulus,
ν
{\displaystyle \nu }
is the Poisson's ratio, and
ε
α
β
{\displaystyle \varepsilon _{\alpha \beta }}
are the in-plane strains. The through-the-thickness shear stresses and strains are related by
σ
31
=
2
G
ε
31
and
σ
32
=
2
G
ε
32
{\displaystyle \sigma _{31}=2G\varepsilon _{31}\quad {\text{and}}\quad \sigma _{32}=2G\varepsilon _{32}}
where
G
=
E
/
(
2
(
1
+
ν
)
)
{\displaystyle G=E/(2(1+\nu ))}
is the shear modulus .
Constitutive relations
edit
The relations between the stress resultants and the generalized displacements for an isotropic Mindlin–Reissner plate are:
[
N
11
N
22
N
12
]
=
2
E
h
1
−
ν
2
[
1
ν
0
ν
1
0
0
0
1
−
ν
]
[
u
1
,
1
0
u
2
,
2
0
1
2
(
u
1
,
2
0
+
u
2
,
1
0
)
]
,
{\displaystyle {\begin{bmatrix}N_{11}\\N_{22}\\N_{12}\end{bmatrix}}={\cfrac {2Eh}{1-\nu ^{2}}}{\begin{bmatrix}1&\nu &0\\\nu &1&0\\0&0&1-\nu \end{bmatrix}}{\begin{bmatrix}u_{1,1}^{0}\\u_{2,2}^{0}\\{\frac {1}{2}}~(u_{1,2}^{0}+u_{2,1}^{0})\end{bmatrix}}\,,}
[
M
11
M
22
M
12
]
=
−
2
E
h
3
3
(
1
−
ν
2
)
[
1
ν
0
ν
1
0
0
0
1
−
ν
]
[
φ
1
,
1
φ
2
,
2
1
2
(
φ
1
,
2
+
φ
2
,
1
)
]
,
{\displaystyle {\begin{bmatrix}M_{11}\\M_{22}\\M_{12}\end{bmatrix}}=-{\cfrac {2Eh^{3}}{3(1-\nu ^{2})}}{\begin{bmatrix}1&\nu &0\\\nu &1&0\\0&0&1-\nu \end{bmatrix}}{\begin{bmatrix}\varphi _{1,1}\\\varphi _{2,2}\\{\frac {1}{2}}(\varphi _{1,2}+\varphi _{2,1})\end{bmatrix}}\,,}
and
[
Q
1
Q
2
]
=
κ
G
h
[
w
,
1
0
−
φ
1
w
,
2
0
−
φ
2
]
.
{\displaystyle {\begin{bmatrix}Q_{1}\\Q_{2}\end{bmatrix}}=\kappa Gh{\begin{bmatrix}w_{,1}^{0}-\varphi _{1}\\w_{,2}^{0}-\varphi _{2}\end{bmatrix}}\,.}
The bending rigidity is defined as the quantity
D
=
2
E
h
3
3
(
1
−
ν
2
)
.
{\displaystyle D={\cfrac {2Eh^{3}}{3(1-\nu ^{2})}}\,.}
For a plate of thickness
H
{\displaystyle H}
, the bending rigidity has the form
D
=
E
H
3
12
(
1
−
ν
2
)
.
{\displaystyle D={\cfrac {EH^{3}}{12(1-\nu ^{2})}}\,.}
where
h
=
H
2
{\displaystyle h={\frac {H}{2}}}
If we ignore the in-plane extension of the plate, the governing equations are
M
α
β
,
β
−
Q
α
=
0
Q
α
,
α
+
q
=
0
.
{\displaystyle {\begin{aligned}M_{\alpha \beta ,\beta }-Q_{\alpha }&=0\\Q_{\alpha ,\alpha }+q&=0\,.\end{aligned}}}
In terms of the generalized deformations
w
0
,
φ
1
,
φ
2
{\displaystyle w^{0},\varphi _{1},\varphi _{2}}
, the three governing equations are
∇
2
(
∂
φ
1
∂
x
1
+
∂
φ
2
∂
x
2
)
=
−
q
D
∇
2
w
0
−
∂
φ
1
∂
x
1
−
∂
φ
2
∂
x
2
=
−
q
κ
G
h
∇
2
(
∂
φ
1
∂
x
2
−
∂
φ
2
∂
x
1
)
=
−
2
κ
G
h
D
(
1
−
ν
)
(
∂
φ
1
∂
x
2
−
∂
φ
2
∂
x
1
)
.
{\displaystyle {\begin{aligned}&\nabla ^{2}\left({\frac {\partial \varphi _{1}}{\partial x_{1}}}+{\frac {\partial \varphi _{2}}{\partial x_{2}}}\right)=-{\frac {q}{D}}\\&\nabla ^{2}w^{0}-{\frac {\partial \varphi _{1}}{\partial x_{1}}}-{\frac {\partial \varphi _{2}}{\partial x_{2}}}=-{\frac {q}{\kappa Gh}}\\&\nabla ^{2}\left({\frac {\partial \varphi _{1}}{\partial x_{2}}}-{\frac {\partial \varphi _{2}}{\partial x_{1}}}\right)=-{\frac {2\kappa Gh}{D(1-\nu )}}\left({\frac {\partial \varphi _{1}}{\partial x_{2}}}-{\frac {\partial \varphi _{2}}{\partial x_{1}}}\right)\,.\end{aligned}}}
The boundary conditions along the edges of a rectangular plate are
simply supported
w
0
=
0
,
M
11
=
0
(
oder
M
22
=
0
)
,
φ
1
=
0
(
oder
φ
2
=
0
)
clamped
w
0
=
0
,
φ
1
=
0
,
φ
2
=
0
.
{\displaystyle {\begin{aligned}{\text{simply supported}}\quad &\quad w^{0}=0,M_{11}=0~({\text{or}}~M_{22}=0),\varphi _{1}=0~({\text{or}}~\varphi _{2}=0)\\{\text{clamped}}\quad &\quad w^{0}=0,\varphi _{1}=0,\varphi _{2}=0\,.\end{aligned}}}
Reissner–Stein static theory for isotropic cantilever plates
edit
In general, exact solutions for cantilever plates using plate theory are quite involved and few exact solutions can be found in the literature. Reissner and Stein[ 7] provide a simplified theory for cantilever plates that is an improvement over older theories such as Saint-Venant plate theory.
The Reissner-Stein theory assumes a transverse displacement field of the form
w
(
x
,
y
)
=
w
x
(
x
)
+
y
θ
x
(
x
)
.
{\displaystyle w(x,y)=w_{x}(x)+y\,\theta _{x}(x)\,.}
The governing equations for the plate then reduce to two coupled ordinary differential equations:
b
D
d
4
w
x
d
x
4
=
q
1
(
x
)
−
n
1
(
x
)
d
2
w
x
d
x
2
−
d
n
1
d
x
d
w
x
d
x
−
1
2
d
n
2
d
x
d
θ
x
d
x
−
n
2
(
x
)
2
d
2
θ
x
d
x
2
b
3
D
12
d
4
θ
x
d
x
4
−
2
b
D
(
1
−
ν
)
d
2
θ
x
d
x
2
=
q
2
(
x
)
−
n
3
(
x
)
d
2
θ
x
d
x
2
−
d
n
3
d
x
d
θ
x
d
x
−
n
2
(
x
)
2
d
2
w
x
d
x
2
−
1
2
d
n
2
d
x
d
w
x
d
x
{\displaystyle {\begin{aligned}&bD{\frac {\mathrm {d} ^{4}w_{x}}{\mathrm {d} x^{4}}}=q_{1}(x)-n_{1}(x){\cfrac {d^{2}w_{x}}{dx^{2}}}-{\cfrac {dn_{1}}{dx}}\,{\cfrac {dw_{x}}{dx}}-{\frac {1}{2}}{\cfrac {dn_{2}}{dx}}\,{\cfrac {d\theta _{x}}{dx}}-{\frac {n_{2}(x)}{2}}{\cfrac {d^{2}\theta _{x}}{dx^{2}}}\\&{\frac {b^{3}D}{12}}\,{\frac {\mathrm {d} ^{4}\theta _{x}}{\mathrm {d} x^{4}}}-2bD(1-\nu ){\cfrac {d^{2}\theta _{x}}{dx^{2}}}=q_{2}(x)-n_{3}(x){\cfrac {d^{2}\theta _{x}}{dx^{2}}}-{\cfrac {dn_{3}}{dx}}\,{\cfrac {d\theta _{x}}{dx}}-{\frac {n_{2}(x)}{2}}\,{\cfrac {d^{2}w_{x}}{dx^{2}}}-{\frac {1}{2}}{\cfrac {dn_{2}}{dx}}\,{\cfrac {dw_{x}}{dx}}\end{aligned}}}
where
q
1
(
x
)
=
∫
−
b
/
2
b
/
2
q
(
x
,
y
)
d
y
,
q
2
(
x
)
=
∫
−
b
/
2
b
/
2
y
q
(
x
,
y
)
d
y
,
n
1
(
x
)
=
∫
−
b
/
2
b
/
2
n
x
(
x
,
y
)
d
y
n
2
(
x
)
=
∫
−
b
/
2
b
/
2
y
n
x
(
x
,
y
)
d
y
,
n
3
(
x
)
=
∫
−
b
/
2
b
/
2
y
2
n
x
(
x
,
y
)
d
y
.
{\displaystyle {\begin{aligned}q_{1}(x)&=\int _{-b/2}^{b/2}q(x,y)\,{\text{d}}y~,~~q_{2}(x)=\int _{-b/2}^{b/2}y\,q(x,y)\,{\text{d}}y~,~~n_{1}(x)=\int _{-b/2}^{b/2}n_{x}(x,y)\,{\text{d}}y\\n_{2}(x)&=\int _{-b/2}^{b/2}y\,n_{x}(x,y)\,{\text{d}}y~,~~n_{3}(x)=\int _{-b/2}^{b/2}y^{2}\,n_{x}(x,y)\,{\text{d}}y\,.\end{aligned}}}
At
x
=
0
{\displaystyle x=0}
, since the beam is clamped, the boundary conditions are
w
(
0
,
y
)
=
d
w
d
x
|
x
=
0
=
0
⟹
w
x
(
0
)
=
d
w
x
d
x
|
x
=
0
=
θ
x
(
0
)
=
d
θ
x
d
x
|
x
=
0
=
0
.
{\displaystyle w(0,y)={\cfrac {dw}{dx}}{\Bigr |}_{x=0}=0\qquad \implies \qquad w_{x}(0)={\cfrac {dw_{x}}{dx}}{\Bigr |}_{x=0}=\theta _{x}(0)={\cfrac {d\theta _{x}}{dx}}{\Bigr |}_{x=0}=0\,.}
The boundary conditions at
x
=
a
{\displaystyle x=a}
are
b
D
d
3
w
x
d
x
3
+
n
1
(
x
)
d
w
x
d
x
+
n
2
(
x
)
d
θ
x
d
x
+
q
x
1
=
0
b
3
D
12
d
3
θ
x
d
x
3
+
[
n
3
(
x
)
−
2
b
D
(
1
−
ν
)
]
d
θ
x
d
x
+
n
2
(
x
)
d
w
x
d
x
+
t
=
0
b
D
d
2
w
x
d
x
2
+
m
1
=
0
,
b
3
D
12
d
2
θ
x
d
x
2
+
m
2
=
0
{\displaystyle {\begin{aligned}&bD{\cfrac {d^{3}w_{x}}{dx^{3}}}+n_{1}(x){\cfrac {dw_{x}}{dx}}+n_{2}(x){\cfrac {d\theta _{x}}{dx}}+q_{x1}=0\\&{\frac {b^{3}D}{12}}{\cfrac {d^{3}\theta _{x}}{dx^{3}}}+\left[n_{3}(x)-2bD(1-\nu )\right]{\cfrac {d\theta _{x}}{dx}}+n_{2}(x){\cfrac {dw_{x}}{dx}}+t=0\\&bD{\cfrac {d^{2}w_{x}}{dx^{2}}}+m_{1}=0\quad ,\quad {\frac {b^{3}D}{12}}{\cfrac {d^{2}\theta _{x}}{dx^{2}}}+m_{2}=0\end{aligned}}}
where
m
1
=
∫
−
b
/
2
b
/
2
m
x
(
y
)
d
y
,
m
2
=
∫
−
b
/
2
b
/
2
y
m
x
(
y
)
d
y
,
q
x
1
=
∫
−
b
/
2
b
/
2
q
x
(
y
)
d
y
t
=
q
x
2
+
m
3
=
∫
−
b
/
2
b
/
2
y
q
x
(
y
)
d
y
+
∫
−
b
/
2
b
/
2
m
x
y
(
y
)
d
y
.
{\displaystyle {\begin{aligned}m_{1}&=\int _{-b/2}^{b/2}m_{x}(y)\,{\text{d}}y~,~~m_{2}=\int _{-b/2}^{b/2}y\,m_{x}(y)\,{\text{d}}y~,~~q_{x1}=\int _{-b/2}^{b/2}q_{x}(y)\,{\text{d}}y\\t&=q_{x2}+m_{3}=\int _{-b/2}^{b/2}y\,q_{x}(y)\,{\text{d}}y+\int _{-b/2}^{b/2}m_{xy}(y)\,{\text{d}}y\,.\end{aligned}}}
Derivation of Reissner–Stein cantilever plate equations
The strain energy of bending of a thin rectangular plate of uniform thickness
h
{\displaystyle h}
is given by
U
=
1
2
∫
0
a
∫
−
b
/
2
b
/
2
D
{
(
∂
2
w
∂
x
2
+
∂
2
w
∂
y
2
)
2
+
2
(
1
−
ν
)
[
(
∂
2
w
∂
x
∂
y
)
2
−
∂
2
w
∂
x
2
∂
2
w
∂
y
2
]
}
d
x
d
y
{\displaystyle U={\frac {1}{2}}\int _{0}^{a}\int _{-b/2}^{b/2}D\left\{\left({\frac {\partial ^{2}w}{\partial x^{2}}}+{\frac {\partial ^{2}w}{\partial y^{2}}}\right)^{2}+2(1-\nu )\left[\left({\frac {\partial ^{2}w}{\partial x\partial y}}\right)^{2}-{\frac {\partial ^{2}w}{\partial x^{2}}}{\frac {\partial ^{2}w}{\partial y^{2}}}\right]\right\}{\text{d}}x{\text{d}}y}
where
w
{\displaystyle w}
is the transverse displacement,
a
{\displaystyle a}
is the length,
b
{\displaystyle b}
is the width,
ν
{\displaystyle \nu }
is the Poisson's
ratio,
E
{\displaystyle E}
is the Young's modulus, and
D
=
E
h
3
12
(
1
−
ν
)
.
{\displaystyle D={\frac {Eh^{3}}{12(1-\nu )}}.}
The potential energy of transverse loads
q
(
x
,
y
)
{\displaystyle q(x,y)}
(per unit length) is
P
q
=
∫
0
a
∫
−
b
/
2
b
/
2
q
(
x
,
y
)
w
(
x
,
y
)
d
x
d
y
.
{\displaystyle P_{q}=\int _{0}^{a}\int _{-b/2}^{b/2}q(x,y)\,w(x,y)\,{\text{d}}x{\text{d}}y\,.}
The potential energy of in-plane loads
n
x
(
x
,
y
)
{\displaystyle n_{x}(x,y)}
(per unit width) is
P
n
=
1
2
∫
0
a
∫
−
b
/
2
b
/
2
n
x
(
x
,
y
)
(
∂
w
∂
x
)
2
d
x
d
y
.
{\displaystyle P_{n}={\frac {1}{2}}\int _{0}^{a}\int _{-b/2}^{b/2}n_{x}(x,y)\,\left({\frac {\partial w}{\partial x}}\right)^{2}\,{\text{d}}x{\text{d}}y\,.}
The potential energy of tip forces
q
x
(
y
)
{\displaystyle q_{x}(y)}
(per unit width), and bending moments
m
x
(
y
)
{\displaystyle m_{x}(y)}
and
m
x
y
(
y
)
{\displaystyle m_{xy}(y)}
(per unit width) is
P
t
=
∫
−
b
/
2
b
/
2
(
q
x
(
y
)
w
(
x
,
y
)
−
m
x
(
y
)
∂
w
∂
x
+
m
x
y
(
y
)
∂
w
∂
y
)
d
x
d
y
.
{\displaystyle P_{t}=\int _{-b/2}^{b/2}\left(q_{x}(y)\,w(x,y)-m_{x}(y)\,{\frac {\partial w}{\partial x}}+m_{xy}(y)\,{\frac {\partial w}{\partial y}}\right){\text{d}}x{\text{d}}y\,.}
A balance of energy requires that the total energy is
W
=
U
−
(
P
q
+
P
n
+
P
t
)
.
{\displaystyle W=U-(P_{q}+P_{n}+P_{t})\,.}
With the Reissener–Stein assumption for the displacement, we have
U
=
∫
0
a
b
D
24
[
12
(
d
2
w
x
d
x
2
)
2
+
b
2
(
d
2
θ
x
d
x
2
)
2
+
24
(
1
−
ν
)
(
d
θ
x
d
x
)
2
]
d
x
,
{\displaystyle U=\int _{0}^{a}{\frac {bD}{24}}\left[12\left({\cfrac {d^{2}w_{x}}{dx^{2}}}\right)^{2}+b^{2}\left({\cfrac {d^{2}\theta _{x}}{dx^{2}}}\right)^{2}+24(1-\nu )\left({\cfrac {d\theta _{x}}{dx}}\right)^{2}\right]\,{\text{d}}x\,,}
P
q
=
∫
0
a
[
(
∫
−
b
/
2
b
/
2
q
(
x
,
y
)
d
y
)
w
x
+
(
∫
−
b
/
2
b
/
2
y
q
(
x
,
y
)
d
y
)
θ
x
]
d
x
,
{\displaystyle P_{q}=\int _{0}^{a}\left[\left(\int _{-b/2}^{b/2}q(x,y)\,{\text{d}}y\right)w_{x}+\left(\int _{-b/2}^{b/2}yq(x,y)\,{\text{d}}y\right)\theta _{x}\right]\,dx\,,}
P
n
=
1
2
∫
0
a
[
(
∫
−
b
/
2
b
/
2
n
x
(
x
,
y
)
d
y
)
(
d
w
x
d
x
)
2
+
(
∫
−
b
/
2
b
/
2
y
n
x
(
x
,
y
)
d
y
)
d
w
x
d
x
d
θ
x
d
x
+
(
∫
−
b
/
2
b
/
2
y
2
n
x
(
x
,
y
)
d
y
)
(
d
θ
x
d
x
)
2
]
d
x
,
{\displaystyle {\begin{aligned}P_{n}&={\frac {1}{2}}\int _{0}^{a}\left[\left(\int _{-b/2}^{b/2}n_{x}(x,y)\,{\text{d}}y\right)\left({\cfrac {dw_{x}}{dx}}\right)^{2}+\left(\int _{-b/2}^{b/2}yn_{x}(x,y)\,{\text{d}}y\right){\cfrac {dw_{x}}{dx}}\,{\cfrac {d\theta _{x}}{dx}}\right.\\&\left.\qquad \qquad +\left(\int _{-b/2}^{b/2}y^{2}n_{x}(x,y)\,{\text{d}}y\right)\left({\cfrac {d\theta _{x}}{dx}}\right)^{2}\right]{\text{d}}x\,,\end{aligned}}}
and
P
t
=
(
∫
−
b
/
2
b
/
2
q
x
(
y
)
d
y
)
w
x
−
(
∫
−
b
/
2
b
/
2
m
x
(
y
)
d
y
)
d
w
x
d
x
+
[
∫
−
b
/
2
b
/
2
(
y
q
x
(
y
)
+
m
x
y
(
y
)
)
d
y
]
θ
x
−
(
∫
−
b
/
2
b
/
2
y
m
x
(
y
)
d
y
)
d
θ
x
d
x
.
{\displaystyle {\begin{aligned}P_{t}&=\left(\int _{-b/2}^{b/2}q_{x}(y)\,{\text{d}}y\right)w_{x}-\left(\int _{-b/2}^{b/2}m_{x}(y)\,{\text{d}}y\right){\cfrac {dw_{x}}{dx}}+\left[\int _{-b/2}^{b/2}\left(yq_{x}(y)+m_{xy}(y)\right)\,{\text{d}}y\right]\theta _{x}\\&\qquad \qquad -\left(\int _{-b/2}^{b/2}ym_{x}(y)\,{\text{d}}y\right){\cfrac {d\theta _{x}}{dx}}\,.\end{aligned}}}
Taking the first variation of
W
{\displaystyle W}
with respect to
(
w
x
,
θ
x
,
x
)
{\displaystyle (w_{x},\theta _{x},x)}
and
setting it to zero gives us the Euler equations
(1)
b
D
d
4
w
x
d
x
4
=
q
1
(
x
)
−
n
1
(
x
)
d
2
w
x
d
x
2
−
d
n
1
d
x
d
w
x
d
x
−
1
2
d
n
2
d
x
d
θ
x
d
x
−
n
2
(
x
)
2
d
2
θ
x
d
x
2
{\displaystyle {\text{(1)}}\qquad bD{\frac {\mathrm {d} ^{4}w_{x}}{\mathrm {d} x^{4}}}=q_{1}(x)-n_{1}(x){\cfrac {d^{2}w_{x}}{dx^{2}}}-{\cfrac {dn_{1}}{dx}}\,{\cfrac {dw_{x}}{dx}}-{\frac {1}{2}}{\cfrac {dn_{2}}{dx}}\,{\cfrac {d\theta _{x}}{dx}}-{\frac {n_{2}(x)}{2}}{\cfrac {d^{2}\theta _{x}}{dx^{2}}}}
and
(2)
b
3
D
12
d
4
θ
x
d
x
4
−
2
b
D
(
1
−
ν
)
d
2
θ
x
d
x
2
=
q
2
(
x
)
−
n
3
(
x
)
d
2
θ
x
d
x
2
−
d
n
3
d
x
d
θ
x
d
x
−
n
2
(
x
)
2
d
2
w
x
d
x
2
−
1
2
d
n
2
d
x
d
w
x
d
x
{\displaystyle {\text{(2)}}\qquad {\frac {b^{3}D}{12}}\,{\frac {\mathrm {d} ^{4}\theta _{x}}{\mathrm {d} x^{4}}}-2bD(1-\nu ){\cfrac {d^{2}\theta _{x}}{dx^{2}}}=q_{2}(x)-n_{3}(x){\cfrac {d^{2}\theta _{x}}{dx^{2}}}-{\cfrac {dn_{3}}{dx}}\,{\cfrac {d\theta _{x}}{dx}}-{\frac {n_{2}(x)}{2}}\,{\cfrac {d^{2}w_{x}}{dx^{2}}}-{\frac {1}{2}}{\cfrac {dn_{2}}{dx}}\,{\cfrac {dw_{x}}{dx}}}
where
q
1
(
x
)
=
∫
−
b
/
2
b
/
2
q
(
x
,
y
)
d
y
,
q
2
(
x
)
=
∫
−
b
/
2
b
/
2
y
q
(
x
,
y
)
d
y
,
n
1
(
x
)
=
∫
−
b
/
2
b
/
2
n
x
(
x
,
y
)
d
y
n
2
(
x
)
=
∫
−
b
/
2
b
/
2
y
n
x
(
x
,
y
)
d
y
,
n
3
(
x
)
=
∫
−
b
/
2
b
/
2
y
2
n
x
(
x
,
y
)
d
y
.
{\displaystyle {\begin{aligned}q_{1}(x)&=\int _{-b/2}^{b/2}q(x,y)\,{\text{d}}y~,~~q_{2}(x)=\int _{-b/2}^{b/2}y\,q(x,y)\,{\text{d}}y~,~~n_{1}(x)=\int _{-b/2}^{b/2}n_{x}(x,y)\,{\text{d}}y\\n_{2}(x)&=\int _{-b/2}^{b/2}y\,n_{x}(x,y)\,{\text{d}}y~,~~n_{3}(x)=\int _{-b/2}^{b/2}y^{2}\,n_{x}(x,y)\,{\text{d}}y.\end{aligned}}}
Since the beam is clamped at
x
=
0
{\displaystyle x=0}
, we have
w
(
0
,
y
)
=
d
w
d
x
|
x
=
0
=
0
⟹
w
x
(
0
)
=
d
w
x
d
x
|
x
=
0
=
θ
x
(
0
)
=
d
θ
x
d
x
|
x
=
0
=
0
.
{\displaystyle w(0,y)={\cfrac {dw}{dx}}{\Bigr |}_{x=0}=0\qquad \implies \qquad w_{x}(0)={\cfrac {dw_{x}}{dx}}{\Bigr |}_{x=0}=\theta _{x}(0)={\cfrac {d\theta _{x}}{dx}}{\Bigr |}_{x=0}=0\,.}
The boundary conditions at
x
=
a
{\displaystyle x=a}
can be found by integration by parts:
b
D
d
3
w
x
d
x
3
+
n
1
(
x
)
d
w
x
d
x
+
n
2
(
x
)
d
θ
x
d
x
+
q
x
1
=
0
b
3
D
12
d
3
θ
x
d
x
3
+
[
n
3
(
x
)
−
2
b
D
(
1
−
ν
)
]
d
θ
x
d
x
+
n
2
(
x
)
d
w
x
d
x
+
t
=
0
b
D
d
2
w
x
d
x
2
+
m
1
=
0
,
b
3
D
12
d
2
θ
x
d
x
2
+
m
2
=
0
{\displaystyle {\begin{aligned}&bD{\cfrac {d^{3}w_{x}}{dx^{3}}}+n_{1}(x){\cfrac {dw_{x}}{dx}}+n_{2}(x){\cfrac {d\theta _{x}}{dx}}+q_{x1}=0\\&{\frac {b^{3}D}{12}}{\cfrac {d^{3}\theta _{x}}{dx^{3}}}+\left[n_{3}(x)-2bD(1-\nu )\right]{\cfrac {d\theta _{x}}{dx}}+n_{2}(x){\cfrac {dw_{x}}{dx}}+t=0\\&bD{\cfrac {d^{2}w_{x}}{dx^{2}}}+m_{1}=0\quad ,\quad {\frac {b^{3}D}{12}}{\cfrac {d^{2}\theta _{x}}{dx^{2}}}+m_{2}=0\end{aligned}}}
where
m
1
=
∫
−
b
/
2
b
/
2
m
x
(
y
)
d
y
,
m
2
=
∫
−
b
/
2
b
/
2
y
m
x
(
y
)
d
y
,
q
x
1
=
∫
−
b
/
2
b
/
2
q
x
(
y
)
d
y
t
=
q
x
2
+
m
3
=
∫
−
b
/
2
b
/
2
y
q
x
(
y
)
d
y
+
∫
−
b
/
2
b
/
2
m
x
y
(
y
)
d
y
.
{\displaystyle {\begin{aligned}m_{1}&=\int _{-b/2}^{b/2}m_{x}(y)\,{\text{d}}y~,~~m_{2}=\int _{-b/2}^{b/2}y\,m_{x}(y)\,{\text{d}}y~,~~q_{x1}=\int _{-b/2}^{b/2}q_{x}(y)\,{\text{d}}y\\t&=q_{x2}+m_{3}=\int _{-b/2}^{b/2}y\,q_{x}(y)\,{\text{d}}y+\int _{-b/2}^{b/2}m_{xy}(y)\,{\text{d}}y.\end{aligned}}}
^ Timoshenko, S. and Woinowsky-Krieger, S. "Theory of plates and shells". McGraw–Hill New York, 1959.
^ A. E. H. Love, On the small free vibrations and deformations of elastic shells , Philosophical trans. of the Royal Society (London), 1888, Vol. série A, N° 17 p. 491–549.
^ Reddy, J. N., 2007, Theory and analysis of elastic plates and shells , CRC Press, Taylor and Francis.
^ Uflyand, Ya. S.,1948, Wave Propagation by Transverse Vibrations of Beams and Plates, PMM: Journal of Applied Mathematics and Mechanics, Vol. 12, 287-300 (in Russian)
^ Elishakoff ,I.,2020, Handbook on Timoshenko-Ehrenfest Beam and Uflyand-Mindlin Plate Theories , World Scientific, Singapore, ISBN 978-981-3236-51-6
^ R. D. Mindlin, Influence of rotatory inertia and shear on flexural motions of isotropic, elastic plates , Journal of Applied Mechanics, 1951, Vol. 18 p. 31–38.
^ E. Reissner and M. Stein. Torsion and transverse bending of cantilever plates. Technical Note 2369, National Advisory Committee for Aeronautics,Washington, 1951.