Kirchhoff plate theory

Q: Professor, is Kirchhoff plate theory the two-dimensional version of Euler-Bernoulli beam theory?

Precisely. Just as the Euler-Bernoulli beam assumes 'cross-sections remain orthogonal to the neutral axis,' Kirchhoff plate theory assumes 'straight lines normal to the mid-surface remain straight and normal to the mid-surface after deformation.'

Q: $\nabla^4$ is a fourth-order differential operator. It's the 2D version of the beam equation $EI w'''' = q$, correct?

That's correct. Since $\nabla^4 = \nabla^2(\nabla^2)$: $$ \nabla^4 w = \frac{\partial^4 w}{\partial x^4} + 2\frac{\partial^4 w}{\partial x^2 \partial y^2} + \frac{\partial^4 w}{\partial y^4} $$

Category: Structural Analysis | Integrated 2026-04-06

CAE visualization for plate kirchhoff theory - technical simulation diagram

Kirchhoff Plate Theory

Kirchhoff plate theory: Theoretical Foundations

What is Kirchhoff Plate Theory?

🙋

Professor, is Kirchhoff plate theory the 2D version of Euler-Bernoulli beam theory?

🎓

Exactly. Just as the Euler-Bernoulli beam assumes "cross-sections remain perpendicular to the neutral axis," Kirchhoff plate theory assumes "straight lines through the thickness remain straight and perpendicular to the mid-surface after deformation."

Fundamental Assumptions

🎓

Kirchhoff's Assumptions:

1. Straight Normal Assumption — Straight lines perpendicular to the mid-surface before deformation remain perpendicular to the mid-surface after deformation.

2. Non-Extending Normal Assumption — Strain in the thickness direction $\varepsilon_{zz} = 0$.

3. Zero Transverse Shear Strain — $\gamma_{xz} = \gamma_{yz} = 0$.

🙋

Assumption 3 is the same as in Euler-Bernoulli beam theory. Shear deformation is neglected.

🎓

Correct. This assumption means the rotation angles are determined by the derivatives of the out-of-plane deflection:

$$ \theta_x = -\frac{\partial w}{\partial y}, \quad \theta_y = \frac{\partial w}{\partial x} $$

The only degree of freedom is the deflection $w(x,y)$; the rotation angles are not independent variables.

Governing Equation

🎓

The biharmonic equation for plate deflection $w(x,y)$:

$$ D\nabla^4 w = q(x,y) $$

Here, $D = Et^3/(12(1-\nu^2))$ is the plate bending rigidity, and $q$ is the out-of-plane distributed load.

🙋

$\nabla^4$ is a fourth-order differential operator. It's the 2D version of the beam equation $EI w'''' = q$, right?

🎓

Exactly. Since $\nabla^4 = \nabla^2(\nabla^2)$:

$$ \nabla^4 w = \frac{\partial^4 w}{\partial x^4} + 2\frac{\partial^4 w}{\partial x^2 \partial y^2} + \frac{\partial^4 w}{\partial y^4} $$

Bending Moments and Shear Forces

🎓

Internal force components:

$$ M_x = -D\left(\frac{\partial^2 w}{\partial x^2} + \nu\frac{\partial^2 w}{\partial y^2}\right) $$

$$ M_y = -D\left(\frac{\partial^2 w}{\partial y^2} + \nu\frac{\partial^2 w}{\partial x^2}\right) $$

$$ M_{xy} = -D(1-\nu)\frac{\partial^2 w}{\partial x \partial y} $$

🙋

Moments are determined by the second derivatives of $w$. Bending moment is obtained by differentiating deflection twice. It has the same structure as beams.

🎓

Yes. Bending stresses in the plate are:

$$ \sigma_x = \frac{12 M_x z}{t^3}, \quad \sigma_y = \frac{12 M_y z}{t^3} $$

Maximum stress occurs at the plate surfaces ($z = \pm t/2$). It has the same structure as the beam formula $\sigma = My/I$.

Applicability Range

🙋

Up to what plate thickness can Kirchhoff plate theory be used?

🎓

Thin plates are a prerequisite. A rule of thumb is $b/t > 20$ ($b$: shorter side of the plate, $t$: thickness).

Similar to Timoshenko beam theory:

$b/t > 20$: Kirchhoff plate theory is sufficient.
$10 < b/t < 20$: Consider Mindlin plate theory.
$b/t < 10$: Use Mindlin plate theory or solid elements.

🙋

Is Mindlin plate theory the 2D version of Timoshenko beam theory?

🎓

Exactly. Kirchhoff plate = 2D version of EB beam, Mindlin plate = 2D version of Timoshenko beam. The difference is whether shear deformation is considered.

Kirchhoff Plate Elements in FEM

🙋

Is implementing Kirchhoff plates in FEM difficult?

🎓

Actually, it's very difficult. Because Kirchhoff plate theory involves fourth-order derivatives of $w$, FEM implementation requires $C^1$ continuity (both displacement and rotation angles must be continuous across elements). Standard FEM ($C^0$ continuity) cannot satisfy this.

🙋

Is $C^1$ continuity that difficult?

🎓

Creating polynomial elements that achieve $C^1$ continuity in 2D is challenging. Historically, the Argyris triangle (21 DOF) and Bell triangle (18 DOF) were developed, but they have many degrees of freedom and are not practical. This is why Mindlin plate theory (which only requires $C^0$) became mainstream in FEM.

Summary

🙋

Let me summarize Kirchhoff plate theory.

🎓

Key points:

Bending theory for thin plates neglecting shear deformation — 2D version of EB beam.
$D\nabla^4 w = q$ — Biharmonic equation.
Applicable for $b/t > 20$ — Limited to thin plates.
Requires $C^1$ continuity in FEM — Difficult to implement.
Mindlin plates (only $C^0$ needed) are mainstream in practice — Kirchhoff plate elements are rare.

🙋

So, the theory is elegant, but it yielded the spotlight to Mindlin plates because FEM implementation is difficult.

🎓

Yes. However, Kirchhoff plate theory is the foundation for analytical solutions and is indispensable as a reference solution for verifying FEM results. The Navier solution (double Fourier series solution for rectangular plates) is a classical solution based on Kirchhoff theory.

Coffee Break Trivia

Origins of Kirchhoff Plate Theory

Gustav Kirchhoff established plate bending theory in his 1850 paper "Über das Gleichgewicht und die Bewegung einer elastischen Scheibe." His assumptions (thin plate, preservation of straight normals, neglect of mid-surface strain) remain valid today as "Classical Plate Theory (CPT)." Numerical experiments have confirmed that the error of this theory is less than 1% for forming analysis of rolled aluminum thin sheets (t<2mm).

Computational Methods for Kirchhoff plate theory

Implementation of Kirchhoff Plates in FEM

🙋

How was the $C^1$ continuity problem solved?

🎓

Historically, there are three approaches.

1. High-Order Conforming Elements

🎓

Argyris triangle (21 DOF) and HCT triangle (Hsieh-Clough-Tocher, 12 DOF). They fully satisfy $C^1$ continuity but have many DOF. Academically elegant but low practicality.

2. DKT/DKQ Elements (Discrete Kirchhoff)

🎓

DKT (Discrete Kirchhoff Triangle) uses the Mindlin plate theory framework for discretization and "discretely" enforces Kirchhoff constraints (shear strain = 0) at Gauss integration points.

🙋

Applying Kirchhoff conditions after using Mindlin plate elements... clever.

🎓

DKT was proposed by Batoz, Bathe, and Ho in 1980. It has only 3 nodes and 9 DOF ($w, \theta_x, \theta_y$ per node) and is highly accurate. Nastran's CTRIA3 (bending) internally uses a DKT-type formulation.

🎓

DKQ (Discrete Kirchhoff Quadrilateral) is the quadrilateral version with 4 nodes. Similarly practical.

3. Using Mindlin Plate Elements as Thin Plates

🎓

The most practical approach. Use Mindlin plate elements (including shear deformation); if the plate is thin, shear deformation automatically becomes small. For thin plates, results identical to Kirchhoff plate theory are obtained.

🙋

So, in modern FEM, using "dedicated elements" for Kirchhoff plates is not mainstream; instead, Mindlin plate elements are used as substitutes?

🎓

Exactly. General-purpose shell elements in Abaqus, Ansys, and Nastran are all based on Mindlin (Reissner-Mindlin) theory. For thin plates, they converge to Kirchhoff's theoretical solutions.

Utilizing Analytical Solutions

🙋

How are Kirchhoff plate analytical solutions used for FEM verification?

🎓

The most famous analytical solution is the Navier solution. Deflection for a simply supported rectangular plate under uniform load $q$:

$$ w = \frac{16q}{\pi^6 D} \sum_{m=1,3,5}^{\infty} \sum_{n=1,3,5}^{\infty} \frac{1}{mn\left(\frac{m^2}{a^2}+\frac{n^2}{b^2}\right)^2} \sin\frac{m\pi x}{a} \sin\frac{n\pi y}{b} $$

🙋

A double Fourier series...

Rate this article

Thank you for your feedback!

Helpful

More details

Report error