Nonlinear Shallow Shell Model

• Nonlinear shallow shell models are rigorous theories derived as asymptotic limits of 3D nonlinear elasticity to capture large deflections and geometric nonlinearities.
• The models use Γ-convergence, asymptotic expansions, and scaling regimes to derive distinct limit energy functionals and equilibrium equations.
• Applications include snap-through buckling analysis, high-fidelity numerical simulations, and the development of homogenized multiscale models for thin shells.

A nonlinear shallow shell model represents a mathematically rigorous theory describing the large deflections, geometric nonlinearities, and possible material inhomogeneities of thin shells whose midsurfaces display small curvature compared with their lateral extent (“shallow shells”). The leading theoretical frameworks derive these models as asymptotic limits of three-dimensional nonlinear elasticity, showing how the scaling of applied loads and shell geometry select distinct equations and functionals at the limit. Central regimes include the time-dependent von Kármán shell equations and the Marguerre–von Kármán system, with recent analysis and numerical methods enabling both existence/uniqueness proofs and high-fidelity computation. Depending on scaling, material oscillations, and mid-surface convexity, the models admit homogenized extensions, Monge–Ampère constraints, and connections to plate theory.

1. Geometric Foundations and Scaling Regimes

A shallow shell consists of a midsurface $S \subset \mathbb{R}^3$ with small Gaussian curvature ($K \ll 1$). Its reference configuration is defined by a $C^2$-immersion $\theta:\omega\rightarrow\mathbb{R}^3$ for parameter domain $\omega\subset\mathbb{R}^2$, with thickness $h$. The shell domain $$S^h = \{ x + t\, n(x) : x \in S,\ t \in (-h/2, h/2) \}$$ incorporates the normal direction $n(x)$.

Key geometric ingredients include the first fundamental form $a_{\alpha\beta}$, the second fundamental form $b_{\alpha\beta}$, and Christoffel symbols $\Gamma_{\alpha\beta}^{\gamma}$ governing surface covariant derivatives. The shallow shell assumption can be formalized by $\|\theta - \theta_0\|_{C^2(\omega)} = \delta \ll 1$, where $\theta_0(y) = (y,0)$ is a reference flat map.

Scaling regimes depend on the order $\alpha$ such that external normal loads are $\mathcal{O}(h^\alpha)$. The principal cases include:

• $\alpha = 3$: Marguerre–von Kármán regime (Velčić, 2011).
• $\alpha > 3$: linearized shallow shell model.
• $2 < \alpha < 3$: partial lower-bound results.

These regimes dictate both the asymptotic expansions and the specific energy functionals derived from 3D elasticity.

2. Derivation from Nonlinear Elasticity and Γ-Convergence

Starting from the 3D nonlinear elastic energy (frame-indifferent, quadratic-coercive near $SO(3)$), the shell energy per unit thickness is

$$J^h(y) = \int_{S^h} W(\nabla y)\,dx - \int_{S^h} f^h(x)\,y_3(x)\,dx,$$

where $y$ is a deformation, and $W$ admits second-order expansions $Q_3$ and $Q_2$: $Q_2(E) = \min \{ Q_3(F) : F_{\mathrm{tan}} = E \}$ for $2 \times 2$ symmetric tensors.

Compactness arguments (geometric rigidity theorems) enable the extraction of rotation fields and decompositions yielding in-plane and out-of-plane displacements ($u, w$). Proper rescaling delivers bounds and convergence in $W^{1,2}$ and $W^{2,2}$ norms. The limit energies take forms such as: $$E^{MvK}[u,w] = \frac{1}{2}\int_\omega Q_2(E[u,w])\,dx' + \frac{1}{24}\int_\omega Q_2(\nabla^2 w)\,dx' - \int_\omega f_3(x')\,w(x')\,dx',$$ where the nonlinear strain $E[u,w]$ couples $u$, $w$, and mid-surface geometry via terms like $\mathrm{sym}\,\nabla u$, $\nabla w \otimes \nabla f$, and $\frac{1}{2}\nabla w \otimes \nabla w$ (Velčić, 2011, Lewicka et al., 2013).

Γ-convergence provides lower and upper bounds for the energy, ensuring that minimizers of the 3D energy converge (weakly/strongly) to minimizers of the limiting shell model.

3. Equilibrium Equations and Strain Measures

The equilibrium (Euler–Lagrange) equations derived from the asymptotic energy yield coupled PDEs for the membrane and bending behavior. For a general $S$ parameterized in local coordinates $(\xi^1,\xi^2)$, with in-plane displacement $v^\alpha$ and transverse displacement $w$:

• Membrane strain:

$$\gamma_{\alpha\beta} = \frac{1}{2}(v_{\alpha;\beta} + v_{\beta;\alpha}) + w\,b_{\alpha\beta} + \frac{1}{2} w_{,\alpha} w_{,\beta}$$

• Bending strain:

$$\rho_{\alpha\beta} = -w_{;\alpha\beta} + b_{\alpha}^{\ \gamma} v_{\gamma;\beta} + b_{\beta}^{\ \gamma} v_{\gamma;\alpha}$$

Constitutive laws relate stresses to strains: $$N^{\alpha\beta} = A^{\alpha\beta\gamma\delta}\,\gamma_{\gamma\delta}, \qquad M^{\alpha\beta} = D^{\alpha\beta\gamma\delta}\,\rho_{\gamma\delta}$$ where $A, D$ are shell stiffness tensors.

The strong-form PDE system for time-dependent von Kármán shallow shells reads: $$\begin{aligned} & \rho h v^\alpha_{tt} - (N^{\alpha\beta})_{;\beta} + b^{\alpha}_{\ \beta} (M^{\beta\gamma})_{;\gamma} = f^\alpha \ & \rho h w_{tt} + (M^{\alpha\beta})_{;\alpha\beta} + N^{\alpha\beta} b_{\alpha\beta} = f^3 \end{aligned}$$ with appropriate covariant derivatives; static forms set $v_{tt}, w_{tt} = 0$ (Qin et al., 2018).

4. Boundary Conditions, Functional Analysis, and Existence Theory

Natural boundary conditions for shallow shell models include:

• Clamped: $v = 0$, $w = 0$, $w_{,n} = 0$ on $\partial S$ (the shell boundary).
• Free-edge: $N_{\alpha\beta} n_\beta = 0$, $M_{\alpha\beta} n_\beta = 0$ (no imposed traction or bending).

The admissible function space typically combines $W^{1,2}$ for tangential displacements and $W^{2,2}$ for normal components, e.g.,

$$X_0 = \{ n : n^\alpha \in W^{1,2}_0(\Omega),\ n^3 \in W^{2,2}_0(\Omega) \}$$

(Giang, 2024).

Key existence and uniqueness results (“Direct method of the calculus of variations,” coercivity by Korn-type inequalities on surfaces, and lower semicontinuity) guarantee minimizers under suitably small loads (or for special classes of applied force densities). Uniqueness follows under smallness of force, strict convexity of the energy, and boundedness of the admissible space.

Recent advances extend existence to shells of arbitrary geometry without requiring “small curvature” or “special” Christoffel symbol conditions (Giang, 2024).

5. Asymptotic Limit: Plate Models and Convergence

The asymptotic behavior of shallow shell models as the mid-surface curvature approaches zero (becoming a plate) is rigorously treated via Γ-convergence. For a sequence of immersions $\theta_n \to \theta_0$ (flat), shallow shell minimizers $u_{\theta_n}$ converge in $H^1 \times H^2$ to solutions of the limiting plate model: $$J_0(\eta) = \frac{1}{2} \int_\omega \left( \frac{\epsilon^3}{3} a_0^{\alpha\beta\sigma\tau} \partial_{\sigma\tau} \eta_3 \partial_{\alpha\beta} \eta_3 + \epsilon a_0^{\alpha\beta\sigma\tau} E^0_{\sigma\tau}(\eta) E^0_{\alpha\beta}(\eta) \right) dy - \int_\omega p^i \eta_i dy,$$ where $E^0$ is the flat plate strain (Giang et al., 27 Dec 2025). The main theorem guarantees convergence for arbitrarily large applied loads, provided the sequence of mid-surfaces converges in $C^2$.

6. Special Nonlinear and Constrained Shell Theories

For models incorporating Monge–Ampère constraints (as in (Lewicka et al., 2013)), shallow shell energies are minimized over fields $v$ with prescribed determinant of the Hessian: $$\min \left\{ \frac{1}{24} \int_\Omega Q_2(\nabla^2 v - \nabla^2 v_0) dx - \int_\Omega f v dx\ :\ \det\nabla^2 v = \det\nabla^2 v_0 \right\}$$ Such constraints arise directly from asymptotic expansions, enforcing a second-order isometry in the limit. Existence of smooth approximating maps for positive right-hand sides follows from advanced regularity and density results.

Other variants include the nonlinear Budiansky–Sanders model (Giang, 2024) and the Donnell–Vlasov–Mushtari–Galimov–Koiter model (Giang et al., 27 Dec 2025), which expand the admissible shell geometries and loading conditions.

7. Homogenized and Multiscale Nonlinear Shallow Shell Models

Homogenized models incorporate periodic material oscillations at scale $\epsilon$, requiring two-scale convergence and relaxation-cell formulae in the von Kármán scaling (Hornung et al., 2012): $$E^{h,\epsilon}(u^h) = \frac{1}{h^4} \int_{S^h} W(x, r(x)/\epsilon, \nabla u^h(x)) dx$$ The limit functional depends on both the shell thickness and material period ratio ($h/\epsilon\to\rho$), with effective quadratic forms $Q_2^\rho$ computed via cell minimization over admissible strains. In strictly convex shells, the membrane term vanishes and the model reduces to pure bending.

The Euler–Lagrange equations inherit this structure, yielding coupled membrane and bending equilibrium in the strong form: $$\text{(i)}\quad \operatorname{Div}_S N = 0, \qquad \text{(ii)}\quad \operatorname{Div}_S\operatorname{Div}_S M + N : S = 0,$$ with $N$ and $M$ as effective stress and moment tensors.

8. Numerical Analysis and Computational Methods

Efficient numerical methods are developed for the von Kármán-type shallow shell equations, forming fourth-order coupled nonlinear biharmonic systems for the transverse deflection and Airy’s stress function. Second-order finite difference discretizations, iterative Picard, Newton, and trust-region dogleg solvers are validated for systems with various boundary conditions, including mixed types yielding boundary singularities. Local corrections and transition functions mitigate these singularities (Ji et al., 2017).

Applications include calculation of critical thermal loads, snap-through bifurcation curves, and continuation methods for computing buckling energy barrier landscapes, showing consistency with analytic scaling laws established for spherical shells under compressive pressure (Baumgarten et al., 2018).

9. Applications, Limitations, and Extensions

Nonlinear shallow shell models are used to study:

• Snap-through buckling and bifurcation phenomena in pressurized spherical shells.
• Dynamics of thin-walled constructions subject to large transverse loads.
• Multiscale effects in shells with periodic material microstructure.
• Limiting cases approaching plate theory.

Limitations include restrictions to normal “dead” loads of order $h^3$, moderate rotation theory, and lack of full large-deformation coupling. Intermediate scaling regimes ($2<\alpha<3$) remain partially treated.

Recent analytic advancements provide existence and uniqueness without small-curvature restrictions and extend results to general shell geometries and loading conditions (Giang, 2024, Giang et al., 27 Dec 2025).

---

Cited arXiv works: (Velčić, 2011, Lewicka et al., 2013, Hornung et al., 2012, Qin et al., 2018, Ji et al., 2017, Baumgarten et al., 2018, Giang, 2024, Giang et al., 27 Dec 2025).