Overview

The membrane theory of shells models thin shell structures as surfaces where loads are resisted solely by in-plane membrane stresses, such as tension or compression, while neglecting bending moments and transverse shear forces.[1] This approach treats the shell's middle surface as the primary domain for stress analysis, assuming the material is homogeneous and isotropic with small thickness relative to the radii of curvature.[1] The theory simplifies the three-dimensional elasticity problem into a two-dimensional one by considering only strains in the middle surface and disregarding flexural deformations.[3]

The scope of membrane theory is limited to thin-walled shells under loading conditions that induce predominantly membrane action, including distributed loads like internal pressure, axial forces, or self-weight, where the shell thickness is much smaller than its radius of curvature.[1] In this framework, shells achieve equilibrium through tangential force resultants, behaving analogously to inflated membranes without reliance on flexural rigidity.[3] The theory is particularly applicable to shells of revolution, such as cylindrical or spherical components, and assumes small deflections to linearize the governing relations.[1]

A representative example is a spherical dome under self-weight, where meridional and hoop membrane stresses balance the gravitational load, resulting in tensile forces that maintain structural integrity without significant bending.[1] This illustrates the theory's efficiency in predicting stress distributions for pressure vessel heads or storage tanks.[1] The membrane theory originated with foundational work by Eugenio Beltrami in 1882 on equilibrium equations for thin shells, building on 19th-century developments in elasticity theory and providing simplifications for analyzing curved surfaces.[4][2]

Significance

The membrane theory of shells offers significant engineering advantages by simplifying the analysis of thin-walled structures compared to more comprehensive shell theories that include bending effects, thereby reducing computational complexity and facilitating hand calculations for preliminary design stages.[5] This approach assumes predominant in-plane stress states, neglecting transverse shear and bending moments, which streamlines modeling for structures where curvature effectively distributes loads without significant flexural deformation.[6] As a result, it enables engineers to quickly assess feasibility and optimize geometries in resource-constrained environments.

Economically, membrane theory supports the design of lighter and more efficient structures, such as pressure vessels, arched roofs, and storage tanks, by leveraging the inherent stiffness provided by shell curvature under uniform loads, which minimizes material requirements while maintaining load-bearing capacity. Shells governed by membrane action can span large areas with reduced weight and cost compared to equivalent plate structures, as their geometry enhances resistance to transverse loads through natural load paths, leading to substantial savings in material and construction expenses for long-span applications.[5]

In practice, the theory is essential for thin-walled structures in diverse fields, including aerospace (e.g., aircraft fuselages), civil engineering (e.g., domes and hyperbolic paraboloid roofs), and the chemical industry (e.g., cylindrical storage tanks), where bending effects are minimal and membrane stresses dominate under internal pressure or self-weight.[7] Unlike flat plates, which rely heavily on bending for load resistance and thus require thicker sections, curved shells exploit geometric stiffness for superior load distribution, making membrane theory pivotal in realizing efficient, high-performance designs.[8]

Early Foundations

The membrane theory of shells emerged in the 19th century as a foundational simplification of three-dimensional elasticity theory, applied to thin curved structures where bending moments and transverse shear forces are negligible compared to in-plane (membrane) stresses. This approach treated shells as two-dimensional surfaces in equilibrium, deriving equations from principles of linear elasticity and Hooke's law while focusing on middle-surface deformations. Early formulations assumed small deformations and neglected flexural effects, enabling practical analysis of stress states in thin shells under uniform pressure or similar loads.[2]

Eugenio Beltrami contributed significantly in 1882 by deriving the primary equations for the equilibrium of flexible and inextensible surfaces, establishing the overall mathematical framework for thin shell membrane theory from elasticity principles.[2] Building on this, Léon Lecornu in 1880 independently developed complementary equations for the equilibrium of such surfaces, refining the general form of the membrane equations and solidifying their derivation from elastic principles.[2] These works marked the initial formalization of membrane theory, transitioning from general elasticity to specialized shell analysis.

Early conceptual influences included analogies to tensile structures forming minimal surfaces under uniform tension without bending resistance, later exemplified by soap films in form-finding. Augustus Edward Hough Love further advanced the field in his 1892 treatise A Treatise on the Mathematical Theory of Elasticity, where he formalized elastic membrane theory as an integral component of shell analysis, integrating it with Kirchhoff-Love hypotheses for thin elastic shells.[9]

Initial applications of membrane theory centered on the static equilibrium of surfaces of revolution, such as cylinders and spheres, subjected to axisymmetric loads like internal pressure, providing simplified stress distributions for engineering designs like pressure vessels.[10]

Key Advancements

Between the late 19th and early 20th centuries, membrane theory saw significant extensions. In 1912, H. Reissner applied asymptotic integration methods to analyze unsymmetrical loading on shells of revolution. A. Pucher in 1934 incorporated Airy's stress function to handle arbitrary geometries, improving accuracy for non-symmetric cases. By 1938, V. V. Sokolovskii reduced the equations to canonical form, highlighting their properties for efficient membrane stress distributions.[2]

A pivotal advancement in the membrane theory of shells occurred in 1940 with the publication of Theory of Plates and Shells by Stephen Timoshenko and S. Woinowsky-Krieger, which systematically extended classical membrane formulations to address practical engineering challenges, such as the analysis of cylindrical shells under various loading conditions. This work provided detailed derivations and solutions for membrane stresses in shells of revolution, bridging theoretical elasticity with design applications in structures like pressure vessels and silos, thereby establishing a foundational reference for 20th-century shell engineering.[11]

During the mid-20th century, particularly in the 1950s and 1960s, membrane shell theory was integrated with emerging finite element methods (FEM), enabling the discretization of complex geometries into membrane elements for numerical simulation. Pioneering efforts, such as those by Clough and others, incorporated membrane stiffness matrices into FEM frameworks, allowing efficient computation of in-plane stresses in thin shells without bending effects, which revolutionized the analysis of large-scale structures like domes and tanks in civil and aerospace engineering.[12] This integration facilitated the development of commercial software packages that treated shells as assemblies of triangular or quadrilateral membrane elements, markedly improving accuracy for non-developable surfaces.

In the 1990s, Philippe G. Ciarlet and Véronique Lods provided a rigorous mathematical justification for membrane shell models through asymptotic analysis, demonstrating the convergence of three-dimensional elasticity solutions to two-dimensional membrane equations as shell thickness approaches zero. Their seminal work, including the 1996 paper "Asymptotic analysis of linearly elastic shells. I. Justification of membrane shell equations," employed Gamma-convergence techniques to validate the approximations under elliptic regularity assumptions, addressing long-standing concerns about the theory's foundational rigor. This advancement solidified membrane theory's status in applied mathematics, influencing subsequent derivations for elliptic membrane shell problems.[13]

Post-2000 developments have seen the emergence of discrete shell membrane theories tailored for computational geometry and architectural applications, emphasizing integrable structures. For instance, the 2003 discrete shells model by Grinspun et al. discretized bending and stretching energies on triangular meshes, enabling simulations of flexible thin structures while preserving membrane-dominated behaviors in large deformations. Building on this, Schief's 2014 work on integrable structures in discrete shell membrane theory linked equilibrium equations to discrete O-surface geometry, facilitating the design of deployable architectural shells with exact integrability constraints. These contributions have extended membrane theory into computer graphics and parametric architecture, supporting the creation of optimized, lightweight enclosures.[14]

Core Assumptions

The membrane theory of shells represents a simplified framework within thin shell mechanics, where the structural response is dominated by in-plane membrane forces, neglecting effects that would otherwise require consideration of bending or shear. This approach is justified for thin-walled structures under conditions where out-of-plane deformations are minimal, allowing equilibrium to be maintained solely through tensile or compressive forces in the shell's middle surface.[15]

A primary assumption is the neglect of bending moments and transverse shear stresses, positing that all external loads are balanced exclusively by in-plane membrane force resultants NαβN_{\alpha\beta}Nαβ​. This eliminates the need to account for curvature changes that induce bending, treating the shell as incapable of resisting moments, which is valid for sufficiently thin geometries where such effects are orders of magnitude smaller than membrane stresses.[16][15]

Stresses are assumed to be uniformly distributed through the shell's thickness hhh, leading to the definition of membrane force resultants as the integral of these stresses over the thickness, without variation in the transverse direction. Consequently, only these integrated in-plane forces NαβN_{\alpha\beta}Nαβ​ (normal and shear components) characterize the stress state, simplifying the three-dimensional problem to a two-dimensional one on the reference surface. The theory assumes a linear elastic, homogeneous, and isotropic material behavior, with stresses related to strains via Hooke's law.[16][15]

The theory operates under small deformation assumptions, where strains are infinitesimal, and the shell maintains its undeformed middle surface configuration during loading. This linear elastic framework presumes that displacements do not alter the geometry significantly, allowing strain-displacement relations to be linearized and the material to follow Hooke's law without geometric nonlinearities.[16][15]

These assumptions render the theory particularly suitable for load conditions involving pressures or forces normal to the shell surface that do not produce substantial edge bending effects, such as uniform internal pressure on closed shells of revolution. Under such loading, the in-plane forces efficiently counter the applied normals without invoking bending near boundaries or supports.[15]

Geometric and Kinematic Approximations

In the membrane theory of shells, the thin shell approximation is fundamental, positing that the shell thickness hhh is much smaller than the principal radii of curvature RRR of the middle surface (i.e., h≪Rh \ll Rh≪R).[17] This allows neglect of quadratic or higher-order terms in h/Rh/Rh/R or the through-thickness coordinate z/Rz/Rz/R, simplifying the three-dimensional elasticity problem to a two-dimensional surface model while retaining accuracy for dominant membrane effects.[17] The Kirchhoff-Love hypotheses, adapted from plate theory, are employed but modified for membranes by neglecting inextensional bending and transverse shear deformations, focusing solely on in-plane stretching.[17]

Kinematic descriptions in membrane theory center on the middle surface, where deformations are captured by tangential in-plane displacements uαu_\alphauα​ (along curvilinear directions α=1,2\alpha = 1, 2α=1,2) and normal displacement www, with no independent rotations of the surface normal.[17] This assumes that line elements initially normal to the middle surface remain straight and normal after deformation, but without the rotational degrees of freedom present in full shell theories, as bending stiffness is ignored.[17] Consequently, the kinematics reduce to a displacement field where transverse shear strains vanish, and all deformation occurs through extension of the middle surface.[17]

The reference configuration employs orthogonal curvilinear coordinates (α,β)(\alpha, \beta)(α,β) aligned with the principal curvature directions on the middle surface, incorporating metric coefficients aαa_\alphaaα​ and aβa_\betaaβ​ (often denoted as Lamé parameters AAA and BBB) to account for the surface geometry.[17] The line element on the surface is given by ds2=aα2dα2+aβ2dβ2ds^2 = a_\alpha^2 d\alpha^2 + a_\beta^2 d\beta^2ds2=aα2​dα2+aβ2​dβ2, and Gaussian curvature arises naturally from the second fundamental form involving the principal curvatures 1/Rα1/R_\alpha1/Rα​ and 1/Rβ1/R_\beta1/Rβ​.[17] Through-thickness variations are parameterized by ζ=z/(h/2)\zeta = z/(h/2)ζ=z/(h/2), with −1≤ζ≤1-1 \leq \zeta \leq 1−1≤ζ≤1, but membrane approximations integrate these to focus on surface quantities.[17]

Membrane strain measures εαβ\varepsilon_{\alpha\beta}εαβ​ are derived from the middle surface displacements. For the normal components in principal directions,

and the shear component is

neglecting transverse shear and higher-order curvature couplings beyond the thin shell limit.[15] These extensional strains capture both stretching due to tangential displacements and contributions from normal deflection against the shell's curvature, forming the basis for equilibrium in membrane-dominant loading scenarios.[15]

Equilibrium Conditions

In the membrane theory of shells, equilibrium conditions are derived from the principles of statics applied to an infinitesimal element of the shell surface, ensuring force balance in the tangential and normal directions while neglecting bending moments and transverse shear forces. The theory assumes that the shell resists applied loads solely through in-plane membrane stress resultants, which are forces per unit length denoted as NαβN_{\alpha\beta}Nαβ​ in a general curvilinear coordinate system xαx^\alphaxα on the middle surface.[18]

For a general membrane shell, the equilibrium equations in orthogonal curvilinear coordinates are expressed using the covariant form to account for the surface geometry. The tangential equilibrium components take the form

where α,β=1,2\alpha, \beta = 1, 2α,β=1,2, N;βαβN^{\alpha\beta}{;\beta}N;βαβ​ is the covariant derivative ∂βNαβ+ΓλβαNλβ+ΓμββNαμ\partial\beta N^{\alpha\beta} + \Gamma^\alpha_{\lambda\beta} N^{\lambda\beta} + \Gamma^\beta_{\mu\beta} N^{\alpha\mu}∂β​Nαβ+Γλβα​Nλβ+Γμββ​Nαμ (summed over repeated indices), Γβγα\Gamma^\alpha_{\beta\gamma}Γβγα​ are the Christoffel symbols of the second kind capturing the coordinate curvature, and pαp^\alphapα are the tangential components of the surface load per unit area.[19] The normal equilibrium equation balances the membrane forces against the normal load pnp^npn via the surface curvature:

where bαβb_{\alpha\beta}bαβ​ are the components of the second fundamental form representing the principal curvatures.[19] These equations ensure zero net force on the shell element, with NαβN^{\alpha\beta}Nαβ symmetric (Nαβ=NβαN^{\alpha\beta} = N^{\beta\alpha}Nαβ=Nβα) providing three unknowns resolved by the three equilibrium conditions.[18]

In the specific case of shells of revolution under axisymmetric loading (independent of the azimuthal angle θ\thetaθ), the coordinates simplify to the meridional direction ϕ\phiϕ and the geometry features a radius r(ϕ)r(\phi)r(ϕ) from the axis of revolution to the surface point. The membrane stress resultants reduce to the meridional force NϕN_\phiNϕ​ (acting along the meridian, balancing tangential loads in that direction) and the hoop force NθN_\thetaNθ​ (acting circumferentially, resisting radial expansion).[1] The tangential (meridional) equilibrium equation becomes

where r1r_1r1​ is the principal radius of curvature in the meridional direction, pϕp_\phipϕ​ is the meridional component of the surface load, and the term involving cos⁡ϕ\cos\phicosϕ arises from the geometry of the parallel circle.[1] The normal equilibrium equation is

where r2=r/sin⁡ϕr_2 = r / \sin\phir2​=r/sinϕ is the principal radius of curvature in the hoop direction and pnp_npn​ is the normal surface load (positive outward).[1] These NϕN_\phiNϕ​ and NθN_\thetaNθ​ directly balance the tangential and normal components of the applied loads, such as internal pressure or self-weight, through the shell's curvature.[15]

Boundary conditions at the shell edges (e.g., open ends or junctions) are essential to solve the equilibrium equations, typically specifying the membrane forces NϕN_\phiNϕ​ or essential displacements that relate to strains via kinematic relations.[1] For instance, at an edge with meridional angle ϕ0\phi_0ϕ0​, a prescribed meridional force Nϕ(ϕ0)N_\phi(\phi_0)Nϕ​(ϕ0​) or horizontal/vertical components may be applied, ensuring continuity with adjacent structures or supports; poles require finite NϕN_\phiNϕ​ and NθN_\thetaNθ​ through symmetry.[1] These conditions determine integration constants in the solutions, preventing singularities and maintaining global equilibrium.[15]

Kinematic Relations

To complete the governing equations, the strains εαβ\varepsilon_{\alpha\beta}εαβ​ must satisfy kinematic compatibility conditions derived from the displacement field uαu^\alphauα, ensuring the deformed geometry is consistent. The strain tensor is εαβ=12(uα;β+uβ;α)\varepsilon_{\alpha\beta} = \frac{1}{2} (u_{\alpha;\beta} + u_{\beta;\alpha})εαβ​=21​(uα;β​+uβ;α​), where ; denotes covariant derivative. For solvability, especially in closed shells, the compatibility equation ε;βααβ−R γαβαεγβ=0\varepsilon^{\alpha\beta}{;\beta\alpha} - R^\alpha{\ \gamma\alpha\beta} \varepsilon^{\gamma\beta} = 0ε;βααβ​−R γαβα​εγβ=0 must hold (with RRR the Riemann tensor related to Gaussian curvature). In axisymmetric cases, this simplifies to an ordinary differential equation linking εϕ\varepsilon_\phiεϕ​ and εθ\varepsilon_\thetaεθ​.[1]

Constitutive Relations

In the membrane theory of shells, constitutive relations link the membrane stress resultants to the strains through material laws derived from three-dimensional elasticity under simplifying assumptions. These relations assume a thin shell where transverse normal stress is negligible, leading to a plane stress state (σ33=0\sigma_{33} = 0σ33​=0) and integration of stresses through the thickness to obtain two-dimensional resultants.[15]

For isotropic, linearly elastic materials, the constitutive relations follow from Hooke's law in plane stress, relating in-plane stresses to strains and then integrating to yield membrane forces. The stress-strain relations are

(sum on γ=1,2\gamma = 1,2γ=1,2), where EEE is Young's modulus, ν\nuν is Poisson's ratio, σαβ\sigma_{\alpha\beta}σαβ​ are in-plane stresses, and εαβ\varepsilon_{\alpha\beta}εαβ​ are in-plane strains (Greek indices denote surface coordinates 1,2). Integrating through the shell thickness hhh gives the membrane stress resultants Nαβ=∫−h/2h/2σαβ dz=Eh1−ν2[(1−ν)εαβ+νδαβεγγ]N_{\alpha\beta} = \int_{-h/2}^{h/2} \sigma_{\alpha\beta} , dz = \frac{E h}{1 - \nu^2} \left[ (1 - \nu) \varepsilon_{\alpha\beta} + \nu \delta_{\alpha\beta} \varepsilon_{\gamma\gamma} \right]Nαβ​=∫−h/2h/2​σαβ​dz=1−ν2Eh​[(1−ν)εαβ​+νδαβ​εγγ​], assuming constant thickness and mid-surface strains.[15][20]

In a more general framework, the constitutive relations take the tensor form Nαβ=CαβγδεγδN_{\alpha\beta} = C_{\alpha\beta\gamma\delta} \varepsilon_{\gamma\delta}Nαβ​=Cαβγδ​εγδ​, where CαβγδC_{\alpha\beta\gamma\delta}Cαβγδ​ is the 2D tangential elasticity tensor that incorporates material anisotropy and accounts for shell curvature effects through the geometry. This form arises from the symmetry of the stress and strain tensors and reduces to the isotropic case when CαβγδC_{\alpha\beta\gamma\delta}Cαβγδ​ reflects isotropy.[15]

Extensions to orthotropic materials, common in composite shells, replace isotropic constants with direction-dependent properties such as principal moduli Eϕ,EθE_\phi, E_\thetaEϕ​,Eθ​, Poisson's ratios νϕθ,νθϕ\nu_{\phi\theta}, \nu_{\theta\phi}νϕθ​,νθϕ​, and shear modulus GϕθG_{\phi\theta}Gϕθ​. The full orthotropic relations use a stiffness matrix, but in principal directions for shells of revolution (meridional ϕ\phiϕ and circumferential θ\thetaθ), an uncoupled approximation simplifies to Nϕ=EϕhεϕN_\phi = E_\phi h \varepsilon_\phiNϕ​=Eϕ​hεϕ​ and Nθ=EθhεθN_\theta = E_\theta h \varepsilon_\thetaNθ​=Eθ​hεθ​ when Poisson effects are neglected, focusing on dominant directional stiffness.[21]

Strain-Displacement Relations

In the membrane theory of shells, the strain-displacement relations provide the kinematic foundation for linking the displacements of the shell midsurface to the resulting membrane strains, under the assumption of negligible transverse shear and bending effects. These relations are derived within the framework of differential geometry on curved surfaces, where the shell is modeled as a two-dimensional manifold embedded in three-dimensional space. The tangential displacements are denoted by the vector u=(uα)\mathbf{u} = (u_\alpha)u=(uα​) in curvilinear coordinates (α,β)(\alpha, \beta)(α,β) on the midsurface, and www represents the normal displacement. The membrane strain tensor ϵαβ\epsilon_{\alpha\beta}ϵαβ​ captures the in-plane extensions and shears, incorporating the effects of both tangential stretching and the rotation induced by normal deflection due to surface curvature.

The explicit form of the membrane strain tensor in covariant components is given by

where the semicolon denotes the covariant derivative with respect to the surface metric, bαβb_{\alpha\beta}bαβ​ is the second fundamental form (curvature tensor) describing the principal curvatures, and www is the scalar normal displacement. This expression linearizes the Green-Lagrange strain tensor for small deformations, separating the symmetric part of the tangential displacement gradient from the curvature-induced term that arises as the normal displacement alters the surface geometry. The covariant derivative ensures invariance under coordinate transformations, making the formulation applicable to arbitrary shell geometries. This standard relation originates from asymptotic expansions of three-dimensional elasticity for thin structures and is fundamental to linear shell theories.[22]

For shells of revolution, which possess axial symmetry and are commonly analyzed using meridional angle ϕ\phiϕ and azimuthal angle θ\thetaθ as coordinates, the strain-displacement relations simplify into principal directions aligned with the generators and parallels. Here, uϕu_\phiuϕ​ denotes the displacement tangent to the meridian, uθu_\thetauθ​ the circumferential displacement, RϕR_\phiRϕ​ the principal radius of curvature along the meridian, RθR_\thetaRθ​ along the parallel, and rrr the distance from the axis of revolution (parallel radius, r=Rϕsin⁡ϕr = R_\phi \sin \phir=Rϕ​sinϕ for a sphere). For the axisymmetric case (uθ=0u_\theta = 0uθ​=0, no θ\thetaθ-dependence), the meridional strain is

while the hoop strain is

These expressions account for the extension along the curving meridian and the change in circumference of the parallel circles, respectively, with the cot⁡ϕ\cot \phicotϕ term reflecting the geometry of revolution (e.g., for a sphere of radius RRR, r=Rsin⁡ϕr = R \sin \phir=Rsinϕ, Rϕ=Rθ=RR_\phi = R_\theta = RRϕ​=Rθ​=R, yielding uniform ϵϕ=ϵθ=w/R\epsilon_\phi = \epsilon_\theta = w / Rϵϕ​=ϵθ​=w/R for pure radial displacement uϕ=0u_\phi = 0uϕ​=0, w=w =w= constant). They reduce to ordinary differential equations solvable along the meridian.[23]

Under pure membrane action, where bending stiffness is neglected, an inextensibility condition may apply in certain approximations, such as for developable surfaces or pre-stressed membranes. In this context, the change in Gaussian curvature KKK of the deformed midsurface is intrinsically related to the membrane strains through the Gauss-Theorema egregium, which states that KKK is determined solely by the first fundamental form (metric tensor) of the surface. The linearized compatibility equations for the strains ϵαβ\epsilon_{\alpha\beta}ϵαβ​ ensure that the induced metric change preserves surface integrability, linking ∂2ϵαβ/∂xα∂xβ\partial^2 \epsilon_{\alpha\beta} / \partial x^\alpha \partial x^\beta∂2ϵαβ​/∂xα∂xβ (in appropriate components) to ΔK\Delta KΔK, the variation in Gaussian curvature. This geometric constraint is crucial for solvability in membrane problems without bending moments, as it enforces that the strains correspond to a realizable deformation field.

The aforementioned relations rely on the small strain approximation, valid when the gradient of displacements satisfies ∣∇u∣≪1|\nabla \mathbf{u}| \ll 1∣∇u∣≪1 and higher-order terms like (∇u)2(\nabla \mathbf{u})^2(∇u)2 or products involving w/Rw / Rw/R (where RRR is a characteristic radius of curvature) are negligible compared to unity. This linearization simplifies the exact nonlinear kinematic relations from three-dimensional continuum mechanics to a first-order theory, enabling closed-form solutions for many shell geometries while maintaining accuracy for moderately thin structures under small deflections. Constitutive relations then link these strains to membrane stress resultants, but such connections are addressed separately.[24]

Stress Resultants

In the membrane theory of shells, stress resultants represent the integrated in-plane forces per unit length acting on the shell's middle surface, obtained by integrating the three-dimensional stress components through the shell thickness while neglecting bending and shear effects. The components of the membrane stress resultant tensor, NαβN_{\alpha\beta}Nαβ​, are defined as

where σαβ\sigma_{\alpha\beta}σαβ​ denotes the in-plane stress tensor components in the curvilinear coordinates (α,β)(\alpha, \beta)(α,β) tangent to the shell surface, hhh is the shell thickness, and zzz is the coordinate normal to the middle surface. This formulation assumes that σαβ\sigma_{\alpha\beta}σαβ​ varies negligibly through the thickness, justified by the thin-shell approximation (h≪Rh \ll Rh≪R, with RRR the principal radius of curvature) and the exclusion of transverse normal stresses and moments, leading to a pure state of tangential membrane stress.[25]

In orthogonal curvilinear coordinates aligned with the principal curvature directions, the stress resultants simplify to principal components: the meridional resultant N1N_1N1​ (or NϕN_\phiNϕ​) acting along the meridian curve, and the circumferential (hoop) resultant N2N_2N2​ (or NθN_\thetaNθ​) perpendicular to it, with the shear component N12N_{12}N12​ (or NϕθN_{\phi\theta}Nϕθ​) coupling them in non-orthogonal cases. For instance, in a thin-walled cylindrical shell under uniform internal pressure ppp and radius rrr, the hoop stress resultant is Nθ=prN_\theta = p rNθ​=pr, while the longitudinal resultant is Nϕ=pr/2N_\phi = p r / 2Nϕ​=pr/2, illustrating how geometry dictates the distribution of these forces to maintain equilibrium. These principal resultants capture the load-carrying capacity in the tangential plane, with their magnitudes depending on the shell's Gaussian curvature and applied loads.[25][26]

The paths of stress resultants in shells can be visualized as "force flow" lines tracing the tangential transmission of loads across the surface, analogous to streamlines in fluid mechanics, which highlight regions of tension, compression, or shear without invoking three-dimensional stress fields. This representation underscores the membrane theory's focus on in-plane force balance, where resultants follow geodesic paths on the developable surface for minimal energy. Seminal formulations of these concepts trace to Love's integration of elasticity principles for thin shells.[25]

For shells with variable thickness h(α,β)h(\alpha, \beta)h(α,β), the integration for NαβN_{\alpha\beta}Nαβ​ is performed locally over the varying hhh, often approximated by weighting with an average thickness or assuming piecewise constancy to simplify computations; however, classical analyses typically presume uniform hhh to maintain analytical tractability. This assumption aligns with the theory's origins in Love's work and subsequent refinements by Flügge and Timoshenko for practical engineering applications.[25][26]

Shells of Revolution

Shells of revolution in membrane theory are axisymmetric structures formed by rotating a smooth curve, known as the generator or meridian, around a fixed axis in its plane, resulting in a surface of revolution. The geometry is typically parameterized using the arc length sss along the meridian or the angle ϕ\phiϕ measured from the axis of rotation, where ϕ\phiϕ increases along the generator. The principal radii of curvature are the meridional radius R1(s)R_1(s)R1​(s) (or R1(ϕ)R_1(\phi)R1​(ϕ)), which is the radius of curvature of the meridian itself, and the parallel or hoop radius R2(s)R_2(s)R2​(s) (or R2(ϕ)R_2(\phi)R2​(ϕ)), which is the distance from the axis of rotation to the point on the surface, measured perpendicular to the meridian. These radii determine the local curvature and are essential for formulating the equilibrium equations, as the shell's resistance to loads depends on its double curvature. For instance, in a spherical shell, R1=R2=RR_1 = R_2 = RR1​=R2​=R (constant), while in more general cases like ellipsoids or paraboloids, they vary along the meridian.[1]

Under axisymmetric loading conditions, which preserve rotational symmetry about the axis, membrane theory applies to predict the in-plane stress resultants without considering bending effects. Typical load types include normal pressure pn(s)p_n(s)pn​(s) distributed symmetrically over the surface, self-weight acting vertically with components resolved into normal and meridional directions, or axial forces applied along the generator; torsional or shear loads in the circumferential direction are excluded to maintain axisymmetry. The normal load pnp_npn​ is perpendicular to the shell surface, while any meridional tangent load pϕp_\phipϕ​ lies in the plane of the meridian. For self-weight of uniform density γ\gammaγ, the components are pn=−γcos⁡ϕp_n = -\gamma \cos \phipn​=−γcosϕ and pϕ=γsin⁡ϕp_\phi = \gamma \sin \phipϕ​=γsinϕ, ensuring the load vector aligns with gravity. Equilibrium is achieved solely through membrane forces, assuming the shell thickness is small compared to the radii, and no transverse shear or moments are present.[1][27]

The stress distribution in such shells is derived by integrating the equilibrium equations along the generator. The meridional stress resultant Nϕ(s)N_\phi(s)Nϕ​(s) is obtained by resolving forces in the meridional direction and integrating from the apex or a known boundary condition, often incorporating the total vertical load above a given section to balance the resultant. Specifically, for axisymmetric cases, NϕN_\phiNϕ​ satisfies

where r=R2sin⁡ϕr = R_2 \sin \phir=R2​sinϕ is the radial distance from the axis. The hoop stress resultant NθN_\thetaNθ​ is then found from the normal equilibrium equation:

yielding Nθ=R2(pn−NϕR1)N_\theta = R_2 \left( p_n - \frac{N_\phi}{R_1} \right)Nθ​=R2​(pn​−R1​Nϕ​​), which includes terms from the integrated meridional forces and local loading. This relation highlights how the hoop stress depends on the ratio of radii R1/R2R_1 / R_2R1​/R2​ and accumulated effects along the meridian. In general, NϕN_\phiNϕ​ carries the primary vertical resolution of loads, while NθN_\thetaNθ​ provides lateral stability through curvature.[1][27]

A representative example is the parabolic dome under uniform vertical loading, such as self-weight or snow approximated as uniform, where the meridian follows a parabolic profile r=kz2r = k z^2r=kz2. In this case, the meridional stress resultant NϕN_\phiNϕ​ is compressive, starting near zero at the apex and increasing in magnitude toward the base, where support reactions dominate. Such distributions are critical for design, as they indicate regions prone to buckling or requiring reinforcement, and are solved using stress functions tailored to the parabolic geometry.[25]

Cylindrical and Spherical Shells

In membrane theory, cylindrical shells under uniform internal pressure are analyzed assuming an infinite length to neglect end effects, focusing solely on in-plane stress resultants. The hoop stress resultant NθN_\thetaNθ​ is given by Nθ=prN_\theta = p rNθ​=pr, where ppp is the internal pressure, rrr is the radius, and the thickness hhh is implicit in the resultant definition as force per unit length.[28] The longitudinal stress resultant NzN_zNz​ is Nz=pr2N_z = \frac{p r}{2}Nz​=2pr​, reflecting equilibrium along the axis.[28] For finite-length cylinders, pure membrane theory neglects edge effects, providing a uniform stress distribution away from boundaries.[28]

Spherical shells under internal pressure exhibit isotropic membrane stresses due to geometric symmetry. The meridional and hoop stress resultants are equal, Nϕ=Nθ=pr2N_\phi = N_\theta = \frac{p r}{2}Nϕ​=Nθ​=2pr​, for thin-walled vessels where the radius rrr greatly exceeds thickness hhh.[29] This uniform biaxial tension arises from balancing the pressure force on a hemispherical section against the resultant forces in the wall, resulting in stresses half those of the hoop stress in a comparable cylinder.[29]

In contrast, spherical shells under uniform external pressure develop compressive membrane stresses Nϕ=Nθ=−pR2N_\phi = N_\theta = -\frac{p R}{2}Nϕ​=Nθ​=−2pR​, but membrane theory does not account for stability, and buckling instability occurs at much lower pressures. For shallow spherical caps (small angular extent α≪1\alpha \ll 1α≪1), the Donnell–Mushtari–Vlasov (DMV) theory, a simplified nonlinear shallow shell theory, is used to analyze buckling, post-buckling, thermal instability, vibrations, and static bistability (two stable states: natural and everted). The DMV theory neglects certain nonlinear terms (e.g., square of normal rotation in in-plane strains) and assumes short-wavelength deformations relative to the radius of curvature, suitable for small to moderate deflections and shallow geometries. For perfect shells, it yields the classical critical buckling pressure pcr=2E3(1−ν2)(hR)2p_{cr} = \frac{2E}{\sqrt{3(1-\nu^2)}} \left( \frac{h}{R} \right)^2pcr​=3(1−ν2)​2E​(Rh​)2, where EEE is Young's modulus, ν\nuν is Poisson's ratio, hhh is thickness, and RRR is radius.[30][31]

Practical applications of membrane theory dominate in boiler design, where cylindrical shells withstand internal steam pressure primarily through hoop and longitudinal stresses, guiding thickness selection for safety.[1] For spherical geometries, such as storage tanks or balloons, the theory predicts uniform membrane stresses under pressure, enabling efficient material use in scenarios like liquefied gas containment where bending is minimal.[1]

Analytical Solutions

In the membrane theory of shells, analytical solutions are obtained through direct integration or series expansions for cases with simple geometries and symmetric or periodic loads, providing closed-form expressions for stress resultants without bending effects.

For shells of revolution under axisymmetric loading, integration techniques exploit the meridional equilibrium equation to determine the meridional stress resultant NϕN_\phiNϕ​. Considering the free body diagram of the shell portion from the apex to arc length sss along the meridian, the resultant Nϕ(s)N_\phi(s)Nϕ​(s) balances the integrated component of the applied load in the direction parallel to the axis of revolution: Nϕ(s)=−12πrsin⁡α∫0s2πr′pz ds′N_\phi(s) = -\frac{1}{2\pi r \sin\alpha} \int_0^s 2\pi r' p_z , ds'Nϕ​(s)=−2πrsinα1​∫0s​2πr′pz​ds′, where rrr is the radius of the parallel circle at sss, α\alphaα is the angle between the meridian tangent and the axis, and pzp_zpz​ is the axial load component per unit area. The hoop stress resultant NθN_\thetaNθ​ is then found from the normal equilibrium equation: Nϕr1+Nθr2=pn\frac{N_\phi}{r_1} + \frac{N_\theta}{r_2} = p_nr1​Nϕ​​+r2​Nθ​​=pn​, with r1r_1r1​ and r2r_2r2​ the principal radii of curvature and pnp_npn​ the normal load component. This "marching" approach along the meridian yields exact solutions for geometries defined parametrically, such as spheres or cones, under uniform pressure or self-weight.[](Timoshenko and Woinowsky-Krieger 1959)[](Fl%C3%BCgge 1967)

Exact solutions are particularly notable for hyperbolic shells under self-weight, where the geometry approximates an inverted catenary to achieve pure compression. For such a shell with line load www due to self-weight, the meridional force simplifies to Nϕ=−wr/tan⁡αN_\phi = -w r / \tan\alphaNϕ​=−wr/tanα, with rrr the parallel radius and α\alphaα the angle of the tangent to the meridian with the horizontal; the hoop force NθN_\thetaNθ​ follows from normal balance, ensuring no bending moments throughout. This form arises from the equilibrium of vertical forces, where the shell shape satisfies the differential equation for uniform thickness under gravity, as applied in cooling tower design.[](Zingoni 1997)

For cylindrical shells subjected to non-axisymmetric loads, series methods decompose the loading into a Fourier series in the circumferential coordinate θ\thetaθ: p( x, θ)=∑n=0∞(ancos⁡nθ+bnsin⁡nθ)eikxp(\ x,\ \theta) = \sum_{n=0}^\infty (a_n \cos n\theta + b_n \sin n\theta) e^{ikx}p( x, θ)=∑n=0∞​(an​cosnθ+bn​sinnθ)eikx or similar modal forms along the axial coordinate xxx. Substituting into the governing partial differential equation reduces the problem to a set of ordinary differential equations for each mode nnn, solvable analytically using exponential or hyperbolic functions for infinite cylinders, or boundary matching for finite lengths. These solutions capture circumferential variations efficiently for wind or thermal loads.[](Timoshenko and Woinowsky-Krieger 1959)

Boundary value problems in membrane theory, such as specified NϕN_\phiNϕ​ or displacements at shell edges, are addressed via direct integration for axisymmetric cases or the method of characteristics for the resulting hyperbolic partial differential equations in general geometries. For instance, in shells of revolution, edge conditions like zero radial displacement lead to characteristic curves along which Riemann invariants propagate, allowing step-by-step solution from boundaries inward; compatibility with constitutive relations ensures unique stress fields.[](Niordson 1984)

Numerical Approaches

The finite element method (FEM) serves as a primary numerical approach for analyzing membrane shells, particularly for complex geometries where analytical solutions are infeasible. Membrane elements in FEM are typically formulated as thin surfaces capable of carrying in-plane stresses without bending stiffness, often using 3-node triangular or 4- to 6-node quadrilateral elements to discretize the shell midsurface. The element stiffness matrix is derived from the integral K=∫ABTCB dA\mathbf{K} = \int_A \mathbf{B}^T \mathbf{C} \mathbf{B} , dAK=∫A​BTCBdA, where B\mathbf{B}B is the strain-displacement matrix incorporating membrane strains, C\mathbf{C}C is the constitutive matrix relating stresses to strains, and the integration is performed over the element area AAA.[32][33]

Discrete formulations provide an alternative numerical framework for membrane shell equilibrium, leveraging concepts from discrete differential geometry to model integrable discrete shell membranes. These approaches often employ bar networks or quadrilateral "plated" elements inscribed in circles, ensuring equilibrium through discrete analogues of continuous integrability conditions derived from O-surface theory. Such methods are particularly useful for simulating large deformations while preserving geometric constraints inherent to membrane theory.[34]

In practical implementations, commercial software like ANSYS and ABAQUS facilitates membrane shell analysis using dedicated membrane-only elements, which support arbitrary shapes through surface meshing and handle nonlinear effects such as geometric nonlinearity. For instance, ANSYS shell elements can be configured for pure membrane behavior by suppressing bending and shear terms, while ABAQUS membrane elements explicitly model in-plane forces without transverse stiffness. These tools enable efficient computation for engineering applications like pressure vessels or tensile structures.[35][36]

Asymptotic expansions offer a numerical validation strategy for membrane approximations, particularly in semi-membrane theories that blend pure membrane behavior with minor bending effects through series expansions in the small parameter h/Rh/Rh/R, where hhh is shell thickness and RRR is the principal radius of curvature. These methods systematically improve accuracy for moderately thick shells by integrating three-dimensional elasticity equations asymptotically.[37][38]

Validity and Limitations

The membrane theory of shells is valid under conditions where bending deformations are negligible compared to in-plane stretching, primarily for thin shells with a small thickness-to-radius ratio, such as h/R<1/20h/R < 1/20h/R<1/20 (or R/h>20R/h > 20R/h>20), allowing the neglect of moment resultants while retaining extensional stresses.[39] This approximation holds particularly well for uniform or smoothly distributed loads where the characteristic wavelength of loading significantly exceeds the shell radius RRR, as well as for shallow shells under axisymmetric conditions, ensuring equilibrium through hoop and meridional forces without significant curvature changes.[23] For instance, in cylindrical or spherical shells subjected to internal pressure, the theory accurately predicts stress resultants, aligning closely with experimental measurements for thin-walled pressure vessels where radial expansions dominate.[39]

However, the theory fails near edges, supports, and discontinuities, where incompatible boundary conditions induce localized bending moments in boundary layers, typically spanning 10-20% of the shell height and decaying exponentially away from the edge.[23] It also breaks down under compressive loads, as thin shells exhibit buckling instability before reaching significant membrane compression, limiting its use in stability analyses without supplementary bending considerations.[40] In regions dominated by bending, such as near concentrated loads or abrupt geometric changes, the theory overpredicts structural stiffness by ignoring transverse shear and rotary inertia effects.[39]

Further limitations arise for thicker shells where h/R>1/10h/R > 1/10h/R>1/10, as bending contributions become prominent, invalidating the pure membrane state assumption; similarly, dynamic or non-uniform loads, like those from wind on roofs, introduce local bending that causes deviations from predictions.[39] While extensions incorporating bending terms can address these shortcomings, membrane theory remains a foundational approximation best suited to extensional-dominated responses in thin, smoothly loaded structures.[23]

Extensions to Bending Theories

The bending theory of shells extends the membrane theory by incorporating resistance to transverse shear deformation and flexural rigidity, which are crucial for capturing localized effects such as edge bending and boundary disturbances. In this framework, moment resultants MαβM_{\alpha\beta}Mαβ​ (bending and twisting moments) and Kirchhoff shear forces QαQ_\alphaQα​ (effective transverse shears derived from moment equilibrium) are introduced alongside the membrane stress resultants NαβN_{\alpha\beta}Nαβ​. These additions lead to a system of fifth-order partial differential equations (PDEs) governing the shell's deformation, in contrast to the second-order PDEs of pure membrane theory, enabling the analysis of bending-dominated behaviors in thin shells under non-uniform loading or at boundaries.[41]

Semi-membrane approaches represent hybrid models that blend membrane assumptions with selective inclusion of bending stiffness, particularly suited for structures like long cylindrical shells where global behavior is membrane-like but local ring bending occurs at the ends or under concentrated loads. A prominent example is the Donnell-Mushtari-Vlasov (DMV) theory, a simplified nonlinear shallow shell theory that approximates the shell equations by neglecting certain nonlinear terms (e.g., the square of the normal rotation in the in-plane strains) and assuming short-wavelength deformations relative to the radius of curvature. This makes it suitable for small to moderate deflections and shallow geometries (angular extent α ≪ 1), including thin spherical caps (shallow portions of spheres). The theory derives coupled fourth-order PDEs for the normal displacement www and Airy stress function FFF (often denoted φ), typically in cylindrical coordinates. These equations are widely used for buckling, post-buckling, thermal instability, vibrations, and static bistability (two stable states: natural and everted) of spherical caps, as well as buckling and vibration analyses of cylindrical and nearly flat shells, providing a computationally efficient approximation valid when the shell's rise is small compared to its planform dimensions. For perfect spherical shells, it yields the classical critical buckling pressure under uniform external pressure: pcr=2E3(1−ν2)(tR)2p_{\mathrm{cr}} = 2E \sqrt{3(1-\nu^{2})} \left( \frac{t}{R} \right)^{2}pcr​=2E3(1−ν2)​(Rt​)2, where EEE is Young's modulus, ν\nuν is Poisson's ratio, ttt is thickness, and RRR is radius.[42][30]

More rigorous linearized extensions, such as the Koiter and Sanders theories, build on the Kirchhoff-Love hypotheses to formulate bending-inclusive models for shallow and moderately curved shells, emphasizing variational principles for stability and small deformations. Koiter's theory derives a consistent first-order approximation through asymptotic expansion of the three-dimensional elasticity equations, yielding energy functionals that combine quadratic membrane strain energy with bending energy terms proportional to the square of curvature changes, ∫(ϵαβϵαβ+καβκαβ) dA\int ( \epsilon_{\alpha\beta} \epsilon^{\alpha\beta} + \kappa_{\alpha\beta} \kappa^{\alpha\beta} ) , dA∫(ϵαβ​ϵαβ+καβ​καβ)dA, where ϵ\epsilonϵ and κ\kappaκ denote membrane and bending strains, respectively; this approach is particularly effective for post-buckling analysis and error estimation in thin shells. Sanders' nonlinear extension refines this by incorporating moderate rotations and improved strain-displacement relations, avoiding inconsistencies in earlier approximations and facilitating applications to imperfection-sensitive buckling in cylindrical shells. Both theories are valid for shallow geometries where the Gaussian curvature is small, offering enhanced accuracy over DMV for nonlinear regimes without the full complexity of sixth-order shell equations.[43][44]

In transitioning from membrane to bending theories, the membrane solution serves as a zeroth-order approximation for the interior of the shell, with bending perturbations introduced asymptotically to resolve boundary layers—narrow regions near edges where rapid variations in deflection and stresses occur due to compatibility requirements. The characteristic width of these boundary layers scales as Rh\sqrt{R h}Rh​, where RRR is the principal radius of curvature and hhh is the shell thickness, allowing membrane theory to approximate the bulk while bending theory corrects localized effects; this perturbation method, rooted in matched asymptotic expansions, underpins many semi-membrane and shallow shell analyses for efficient yet accurate predictions.[41][45]

Foundational Texts

The foundational texts of membrane theory for shells emerged in the late 19th and early 20th centuries, building on principles of elasticity to describe stress states in thin curved structures where bending moments are negligible. Early contributions focused on the equilibrium of flexible, inextensible surfaces, establishing the core equations for membrane forces. Eugenio Beltrami's 1882 work, Sull’equilibrio delle superficiali flessibili ed inestendibili, provided one of the first mathematical treatments of equilibrium in curved elastic surfaces, deriving equations that relate membrane stresses to surface geometry and applied loads without considering flexural rigidity.[2] This paper laid groundwork for analyzing thin shells under uniform tension, emphasizing in-plane forces that maintain structural integrity in momentless states.[2]

Complementing Beltrami, Léon Lecornu's 1885 publication, Sur l'équilibre des surfaces flexibles et inextendibles in the Journal de l'École Polytechnique, extended these ideas to thin shells by formulating equilibrium conditions for flexible surfaces under distributed loads.[46] Lecornu derived differential equations governing membrane stresses, assuming negligible thickness and no resistance to bending, which proved essential for shells like domes and cylinders where membrane action dominates.[2] Together, these papers represented pioneering efforts in membrane theory, prioritizing simplified stress analysis for engineering applications over full three-dimensional elasticity.[47]

Augustus Edward Hough Love's A Treatise on the Mathematical Theory of Elasticity (1892) marked the first systematic integration of membrane theory into a broader elastic framework for shells. Love built on prior elasticity works by Cauchy and Kirchhoff, introducing approximations for thin shells where normal stresses are ignored and midsurface strains govern behavior.[47] His derivations yielded general equations for membrane forces in curved surfaces, applicable to arbitrary geometries, and established the Love-Kirchhoff hypotheses that remain central to the field.[47] This text shifted focus toward practical computation of deformations in elastic membranes, influencing subsequent shell analyses.

Stephen Timoshenko and S. Woinowsky-Krieger's Theory of Plates and Shells (1940) offered practical formulations tailored for engineers, emphasizing membrane solutions for common geometries like cylinders and spheres.[41] The book derived explicit equations for axisymmetric loading, such as hoop and meridional stresses in spherical shells under internal pressure, assuming pure membrane action without bending.[41] Timoshenko's approach simplified Love's general theory into solvable forms, incorporating boundary conditions for real-world structures like pressure vessels, and highlighted the theory's limitations near edges where bending effects emerge.[47]

Vasily Z. Vlasov's General Theory of Shells and Its Application in Engineering (1958, translated into English in 1961 by NASA) provided a comprehensive set of membrane equations for general surfaces, advancing beyond earlier works by addressing non-developable geometries. Vlasov formulated coupled differential equations relating principal curvatures to stress resultants, enabling analysis of complex shells like hyperboloids, while maintaining the membrane assumption of zero transverse shear.[48] This text synthesized prior contributions into a unified framework, emphasizing applications in aerospace and civil engineering, and introduced semi-momentless approximations for efficient computation.

Modern Reviews

In the late 20th and early 21st centuries, mathematical analyses have provided rigorous justifications for membrane shell models, particularly through asymptotic methods and variational convergence techniques. Philippe G. Ciarlet's "Mathematical Elasticity" series, culminating in Volume III: Theory of Shells (2000), establishes the foundations of nonlinear elastic membrane shells using Gamma-convergence to derive two-dimensional equations from three-dimensional elasticity, emphasizing the limiting behavior as shell thickness approaches zero.[49] This work highlights how membrane approximations capture in-plane stresses effectively for thin structures while neglecting bending, offering a mathematically sound basis for applications in engineering design.

Frithiof I. Niordson's "Shell Theory" (1985) provides a comprehensive overview distinguishing membrane theory from bending-inclusive models, incorporating numerical examples to illustrate stress distributions in cylindrical and spherical shells under various loads.[50] The text underscores the practical utility of membrane assumptions for large-radius shells, where bending effects are minimal, and includes computational strategies for solving equilibrium equations.

Recent advancements extend membrane theory into discrete and computational realms. For instance, W. K. Schief's 2014 paper on integrable structure in discrete shell membrane theory introduces discrete analogues of classical integrable membrane classes, enabling equilibrium-preserving discretizations suitable for computational design of curved structures like architectural shells. Complementing this, Ciarlet and Véronique Lods' 1996 asymptotic analysis justifies membrane shell equations for linearly elastic cases, deriving Koiter-type models via expansion in thickness and proving their validity under ellipticity conditions.

Engineering handbooks continue to adapt membrane theory for industrial applications. Wilhelm Flügge's "Stresses in Shells" (1967, with later reprints) remains influential for practical calculations of membrane stresses in pressure vessels and roofs, providing approximate formulas and charts for axisymmetric cases while noting limitations near boundaries. These modern reviews collectively synthesize theoretical rigor with applied extensions, bridging classical foundations to contemporary computational and design challenges in shell mechanics.

Overview

The membrane theory of shells models thin shell structures as surfaces where loads are resisted solely by in-plane membrane stresses, such as tension or compression, while neglecting bending moments and transverse shear forces.[1] This approach treats the shell's middle surface as the primary domain for stress analysis, assuming the material is homogeneous and isotropic with small thickness relative to the radii of curvature.[1] The theory simplifies the three-dimensional elasticity problem into a two-dimensional one by considering only strains in the middle surface and disregarding flexural deformations.[3]

The scope of membrane theory is limited to thin-walled shells under loading conditions that induce predominantly membrane action, including distributed loads like internal pressure, axial forces, or self-weight, where the shell thickness is much smaller than its radius of curvature.[1] In this framework, shells achieve equilibrium through tangential force resultants, behaving analogously to inflated membranes without reliance on flexural rigidity.[3] The theory is particularly applicable to shells of revolution, such as cylindrical or spherical components, and assumes small deflections to linearize the governing relations.[1]

A representative example is a spherical dome under self-weight, where meridional and hoop membrane stresses balance the gravitational load, resulting in tensile forces that maintain structural integrity without significant bending.[1] This illustrates the theory's efficiency in predicting stress distributions for pressure vessel heads or storage tanks.[1] The membrane theory originated with foundational work by Eugenio Beltrami in 1882 on equilibrium equations for thin shells, building on 19th-century developments in elasticity theory and providing simplifications for analyzing curved surfaces.[4][2]

Significance

The membrane theory of shells offers significant engineering advantages by simplifying the analysis of thin-walled structures compared to more comprehensive shell theories that include bending effects, thereby reducing computational complexity and facilitating hand calculations for preliminary design stages.[5] This approach assumes predominant in-plane stress states, neglecting transverse shear and bending moments, which streamlines modeling for structures where curvature effectively distributes loads without significant flexural deformation.[6] As a result, it enables engineers to quickly assess feasibility and optimize geometries in resource-constrained environments.

Economically, membrane theory supports the design of lighter and more efficient structures, such as pressure vessels, arched roofs, and storage tanks, by leveraging the inherent stiffness provided by shell curvature under uniform loads, which minimizes material requirements while maintaining load-bearing capacity. Shells governed by membrane action can span large areas with reduced weight and cost compared to equivalent plate structures, as their geometry enhances resistance to transverse loads through natural load paths, leading to substantial savings in material and construction expenses for long-span applications.[5]

In practice, the theory is essential for thin-walled structures in diverse fields, including aerospace (e.g., aircraft fuselages), civil engineering (e.g., domes and hyperbolic paraboloid roofs), and the chemical industry (e.g., cylindrical storage tanks), where bending effects are minimal and membrane stresses dominate under internal pressure or self-weight.[7] Unlike flat plates, which rely heavily on bending for load resistance and thus require thicker sections, curved shells exploit geometric stiffness for superior load distribution, making membrane theory pivotal in realizing efficient, high-performance designs.[8]

Early Foundations

The membrane theory of shells emerged in the 19th century as a foundational simplification of three-dimensional elasticity theory, applied to thin curved structures where bending moments and transverse shear forces are negligible compared to in-plane (membrane) stresses. This approach treated shells as two-dimensional surfaces in equilibrium, deriving equations from principles of linear elasticity and Hooke's law while focusing on middle-surface deformations. Early formulations assumed small deformations and neglected flexural effects, enabling practical analysis of stress states in thin shells under uniform pressure or similar loads.[2]

Eugenio Beltrami contributed significantly in 1882 by deriving the primary equations for the equilibrium of flexible and inextensible surfaces, establishing the overall mathematical framework for thin shell membrane theory from elasticity principles.[2] Building on this, Léon Lecornu in 1880 independently developed complementary equations for the equilibrium of such surfaces, refining the general form of the membrane equations and solidifying their derivation from elastic principles.[2] These works marked the initial formalization of membrane theory, transitioning from general elasticity to specialized shell analysis.

Early conceptual influences included analogies to tensile structures forming minimal surfaces under uniform tension without bending resistance, later exemplified by soap films in form-finding. Augustus Edward Hough Love further advanced the field in his 1892 treatise A Treatise on the Mathematical Theory of Elasticity, where he formalized elastic membrane theory as an integral component of shell analysis, integrating it with Kirchhoff-Love hypotheses for thin elastic shells.[9]

Initial applications of membrane theory centered on the static equilibrium of surfaces of revolution, such as cylinders and spheres, subjected to axisymmetric loads like internal pressure, providing simplified stress distributions for engineering designs like pressure vessels.[10]

Key Advancements

Between the late 19th and early 20th centuries, membrane theory saw significant extensions. In 1912, H. Reissner applied asymptotic integration methods to analyze unsymmetrical loading on shells of revolution. A. Pucher in 1934 incorporated Airy's stress function to handle arbitrary geometries, improving accuracy for non-symmetric cases. By 1938, V. V. Sokolovskii reduced the equations to canonical form, highlighting their properties for efficient membrane stress distributions.[2]

A pivotal advancement in the membrane theory of shells occurred in 1940 with the publication of Theory of Plates and Shells by Stephen Timoshenko and S. Woinowsky-Krieger, which systematically extended classical membrane formulations to address practical engineering challenges, such as the analysis of cylindrical shells under various loading conditions. This work provided detailed derivations and solutions for membrane stresses in shells of revolution, bridging theoretical elasticity with design applications in structures like pressure vessels and silos, thereby establishing a foundational reference for 20th-century shell engineering.[11]

During the mid-20th century, particularly in the 1950s and 1960s, membrane shell theory was integrated with emerging finite element methods (FEM), enabling the discretization of complex geometries into membrane elements for numerical simulation. Pioneering efforts, such as those by Clough and others, incorporated membrane stiffness matrices into FEM frameworks, allowing efficient computation of in-plane stresses in thin shells without bending effects, which revolutionized the analysis of large-scale structures like domes and tanks in civil and aerospace engineering.[12] This integration facilitated the development of commercial software packages that treated shells as assemblies of triangular or quadrilateral membrane elements, markedly improving accuracy for non-developable surfaces.

In the 1990s, Philippe G. Ciarlet and Véronique Lods provided a rigorous mathematical justification for membrane shell models through asymptotic analysis, demonstrating the convergence of three-dimensional elasticity solutions to two-dimensional membrane equations as shell thickness approaches zero. Their seminal work, including the 1996 paper "Asymptotic analysis of linearly elastic shells. I. Justification of membrane shell equations," employed Gamma-convergence techniques to validate the approximations under elliptic regularity assumptions, addressing long-standing concerns about the theory's foundational rigor. This advancement solidified membrane theory's status in applied mathematics, influencing subsequent derivations for elliptic membrane shell problems.[13]

Post-2000 developments have seen the emergence of discrete shell membrane theories tailored for computational geometry and architectural applications, emphasizing integrable structures. For instance, the 2003 discrete shells model by Grinspun et al. discretized bending and stretching energies on triangular meshes, enabling simulations of flexible thin structures while preserving membrane-dominated behaviors in large deformations. Building on this, Schief's 2014 work on integrable structures in discrete shell membrane theory linked equilibrium equations to discrete O-surface geometry, facilitating the design of deployable architectural shells with exact integrability constraints. These contributions have extended membrane theory into computer graphics and parametric architecture, supporting the creation of optimized, lightweight enclosures.[14]

Core Assumptions

The membrane theory of shells represents a simplified framework within thin shell mechanics, where the structural response is dominated by in-plane membrane forces, neglecting effects that would otherwise require consideration of bending or shear. This approach is justified for thin-walled structures under conditions where out-of-plane deformations are minimal, allowing equilibrium to be maintained solely through tensile or compressive forces in the shell's middle surface.[15]

A primary assumption is the neglect of bending moments and transverse shear stresses, positing that all external loads are balanced exclusively by in-plane membrane force resultants NαβN_{\alpha\beta}Nαβ​. This eliminates the need to account for curvature changes that induce bending, treating the shell as incapable of resisting moments, which is valid for sufficiently thin geometries where such effects are orders of magnitude smaller than membrane stresses.[16][15]

Stresses are assumed to be uniformly distributed through the shell's thickness hhh, leading to the definition of membrane force resultants as the integral of these stresses over the thickness, without variation in the transverse direction. Consequently, only these integrated in-plane forces NαβN_{\alpha\beta}Nαβ​ (normal and shear components) characterize the stress state, simplifying the three-dimensional problem to a two-dimensional one on the reference surface. The theory assumes a linear elastic, homogeneous, and isotropic material behavior, with stresses related to strains via Hooke's law.[16][15]

The theory operates under small deformation assumptions, where strains are infinitesimal, and the shell maintains its undeformed middle surface configuration during loading. This linear elastic framework presumes that displacements do not alter the geometry significantly, allowing strain-displacement relations to be linearized and the material to follow Hooke's law without geometric nonlinearities.[16][15]

These assumptions render the theory particularly suitable for load conditions involving pressures or forces normal to the shell surface that do not produce substantial edge bending effects, such as uniform internal pressure on closed shells of revolution. Under such loading, the in-plane forces efficiently counter the applied normals without invoking bending near boundaries or supports.[15]

Geometric and Kinematic Approximations

In the membrane theory of shells, the thin shell approximation is fundamental, positing that the shell thickness hhh is much smaller than the principal radii of curvature RRR of the middle surface (i.e., h≪Rh \ll Rh≪R).[17] This allows neglect of quadratic or higher-order terms in h/Rh/Rh/R or the through-thickness coordinate z/Rz/Rz/R, simplifying the three-dimensional elasticity problem to a two-dimensional surface model while retaining accuracy for dominant membrane effects.[17] The Kirchhoff-Love hypotheses, adapted from plate theory, are employed but modified for membranes by neglecting inextensional bending and transverse shear deformations, focusing solely on in-plane stretching.[17]

Kinematic descriptions in membrane theory center on the middle surface, where deformations are captured by tangential in-plane displacements uαu_\alphauα​ (along curvilinear directions α=1,2\alpha = 1, 2α=1,2) and normal displacement www, with no independent rotations of the surface normal.[17] This assumes that line elements initially normal to the middle surface remain straight and normal after deformation, but without the rotational degrees of freedom present in full shell theories, as bending stiffness is ignored.[17] Consequently, the kinematics reduce to a displacement field where transverse shear strains vanish, and all deformation occurs through extension of the middle surface.[17]

The reference configuration employs orthogonal curvilinear coordinates (α,β)(\alpha, \beta)(α,β) aligned with the principal curvature directions on the middle surface, incorporating metric coefficients aαa_\alphaaα​ and aβa_\betaaβ​ (often denoted as Lamé parameters AAA and BBB) to account for the surface geometry.[17] The line element on the surface is given by ds2=aα2dα2+aβ2dβ2ds^2 = a_\alpha^2 d\alpha^2 + a_\beta^2 d\beta^2ds2=aα2​dα2+aβ2​dβ2, and Gaussian curvature arises naturally from the second fundamental form involving the principal curvatures 1/Rα1/R_\alpha1/Rα​ and 1/Rβ1/R_\beta1/Rβ​.[17] Through-thickness variations are parameterized by ζ=z/(h/2)\zeta = z/(h/2)ζ=z/(h/2), with −1≤ζ≤1-1 \leq \zeta \leq 1−1≤ζ≤1, but membrane approximations integrate these to focus on surface quantities.[17]

Membrane strain measures εαβ\varepsilon_{\alpha\beta}εαβ​ are derived from the middle surface displacements. For the normal components in principal directions,

and the shear component is

neglecting transverse shear and higher-order curvature couplings beyond the thin shell limit.[15] These extensional strains capture both stretching due to tangential displacements and contributions from normal deflection against the shell's curvature, forming the basis for equilibrium in membrane-dominant loading scenarios.[15]

Equilibrium Conditions

In the membrane theory of shells, equilibrium conditions are derived from the principles of statics applied to an infinitesimal element of the shell surface, ensuring force balance in the tangential and normal directions while neglecting bending moments and transverse shear forces. The theory assumes that the shell resists applied loads solely through in-plane membrane stress resultants, which are forces per unit length denoted as NαβN_{\alpha\beta}Nαβ​ in a general curvilinear coordinate system xαx^\alphaxα on the middle surface.[18]

For a general membrane shell, the equilibrium equations in orthogonal curvilinear coordinates are expressed using the covariant form to account for the surface geometry. The tangential equilibrium components take the form

where α,β=1,2\alpha, \beta = 1, 2α,β=1,2, N;βαβN^{\alpha\beta}{;\beta}N;βαβ​ is the covariant derivative ∂βNαβ+ΓλβαNλβ+ΓμββNαμ\partial\beta N^{\alpha\beta} + \Gamma^\alpha_{\lambda\beta} N^{\lambda\beta} + \Gamma^\beta_{\mu\beta} N^{\alpha\mu}∂β​Nαβ+Γλβα​Nλβ+Γμββ​Nαμ (summed over repeated indices), Γβγα\Gamma^\alpha_{\beta\gamma}Γβγα​ are the Christoffel symbols of the second kind capturing the coordinate curvature, and pαp^\alphapα are the tangential components of the surface load per unit area.[19] The normal equilibrium equation balances the membrane forces against the normal load pnp^npn via the surface curvature:

where bαβb_{\alpha\beta}bαβ​ are the components of the second fundamental form representing the principal curvatures.[19] These equations ensure zero net force on the shell element, with NαβN^{\alpha\beta}Nαβ symmetric (Nαβ=NβαN^{\alpha\beta} = N^{\beta\alpha}Nαβ=Nβα) providing three unknowns resolved by the three equilibrium conditions.[18]

In the specific case of shells of revolution under axisymmetric loading (independent of the azimuthal angle θ\thetaθ), the coordinates simplify to the meridional direction ϕ\phiϕ and the geometry features a radius r(ϕ)r(\phi)r(ϕ) from the axis of revolution to the surface point. The membrane stress resultants reduce to the meridional force NϕN_\phiNϕ​ (acting along the meridian, balancing tangential loads in that direction) and the hoop force NθN_\thetaNθ​ (acting circumferentially, resisting radial expansion).[1] The tangential (meridional) equilibrium equation becomes

where r1r_1r1​ is the principal radius of curvature in the meridional direction, pϕp_\phipϕ​ is the meridional component of the surface load, and the term involving cos⁡ϕ\cos\phicosϕ arises from the geometry of the parallel circle.[1] The normal equilibrium equation is

where r2=r/sin⁡ϕr_2 = r / \sin\phir2​=r/sinϕ is the principal radius of curvature in the hoop direction and pnp_npn​ is the normal surface load (positive outward).[1] These NϕN_\phiNϕ​ and NθN_\thetaNθ​ directly balance the tangential and normal components of the applied loads, such as internal pressure or self-weight, through the shell's curvature.[15]

Boundary conditions at the shell edges (e.g., open ends or junctions) are essential to solve the equilibrium equations, typically specifying the membrane forces NϕN_\phiNϕ​ or essential displacements that relate to strains via kinematic relations.[1] For instance, at an edge with meridional angle ϕ0\phi_0ϕ0​, a prescribed meridional force Nϕ(ϕ0)N_\phi(\phi_0)Nϕ​(ϕ0​) or horizontal/vertical components may be applied, ensuring continuity with adjacent structures or supports; poles require finite NϕN_\phiNϕ​ and NθN_\thetaNθ​ through symmetry.[1] These conditions determine integration constants in the solutions, preventing singularities and maintaining global equilibrium.[15]

Kinematic Relations

To complete the governing equations, the strains εαβ\varepsilon_{\alpha\beta}εαβ​ must satisfy kinematic compatibility conditions derived from the displacement field uαu^\alphauα, ensuring the deformed geometry is consistent. The strain tensor is εαβ=12(uα;β+uβ;α)\varepsilon_{\alpha\beta} = \frac{1}{2} (u_{\alpha;\beta} + u_{\beta;\alpha})εαβ​=21​(uα;β​+uβ;α​), where ; denotes covariant derivative. For solvability, especially in closed shells, the compatibility equation ε;βααβ−R γαβαεγβ=0\varepsilon^{\alpha\beta}{;\beta\alpha} - R^\alpha{\ \gamma\alpha\beta} \varepsilon^{\gamma\beta} = 0ε;βααβ​−R γαβα​εγβ=0 must hold (with RRR the Riemann tensor related to Gaussian curvature). In axisymmetric cases, this simplifies to an ordinary differential equation linking εϕ\varepsilon_\phiεϕ​ and εθ\varepsilon_\thetaεθ​.[1]

Constitutive Relations

In the membrane theory of shells, constitutive relations link the membrane stress resultants to the strains through material laws derived from three-dimensional elasticity under simplifying assumptions. These relations assume a thin shell where transverse normal stress is negligible, leading to a plane stress state (σ33=0\sigma_{33} = 0σ33​=0) and integration of stresses through the thickness to obtain two-dimensional resultants.[15]

For isotropic, linearly elastic materials, the constitutive relations follow from Hooke's law in plane stress, relating in-plane stresses to strains and then integrating to yield membrane forces. The stress-strain relations are

(sum on γ=1,2\gamma = 1,2γ=1,2), where EEE is Young's modulus, ν\nuν is Poisson's ratio, σαβ\sigma_{\alpha\beta}σαβ​ are in-plane stresses, and εαβ\varepsilon_{\alpha\beta}εαβ​ are in-plane strains (Greek indices denote surface coordinates 1,2). Integrating through the shell thickness hhh gives the membrane stress resultants Nαβ=∫−h/2h/2σαβ dz=Eh1−ν2[(1−ν)εαβ+νδαβεγγ]N_{\alpha\beta} = \int_{-h/2}^{h/2} \sigma_{\alpha\beta} , dz = \frac{E h}{1 - \nu^2} \left[ (1 - \nu) \varepsilon_{\alpha\beta} + \nu \delta_{\alpha\beta} \varepsilon_{\gamma\gamma} \right]Nαβ​=∫−h/2h/2​σαβ​dz=1−ν2Eh​[(1−ν)εαβ​+νδαβ​εγγ​], assuming constant thickness and mid-surface strains.[15][20]

In a more general framework, the constitutive relations take the tensor form Nαβ=CαβγδεγδN_{\alpha\beta} = C_{\alpha\beta\gamma\delta} \varepsilon_{\gamma\delta}Nαβ​=Cαβγδ​εγδ​, where CαβγδC_{\alpha\beta\gamma\delta}Cαβγδ​ is the 2D tangential elasticity tensor that incorporates material anisotropy and accounts for shell curvature effects through the geometry. This form arises from the symmetry of the stress and strain tensors and reduces to the isotropic case when CαβγδC_{\alpha\beta\gamma\delta}Cαβγδ​ reflects isotropy.[15]

Extensions to orthotropic materials, common in composite shells, replace isotropic constants with direction-dependent properties such as principal moduli Eϕ,EθE_\phi, E_\thetaEϕ​,Eθ​, Poisson's ratios νϕθ,νθϕ\nu_{\phi\theta}, \nu_{\theta\phi}νϕθ​,νθϕ​, and shear modulus GϕθG_{\phi\theta}Gϕθ​. The full orthotropic relations use a stiffness matrix, but in principal directions for shells of revolution (meridional ϕ\phiϕ and circumferential θ\thetaθ), an uncoupled approximation simplifies to Nϕ=EϕhεϕN_\phi = E_\phi h \varepsilon_\phiNϕ​=Eϕ​hεϕ​ and Nθ=EθhεθN_\theta = E_\theta h \varepsilon_\thetaNθ​=Eθ​hεθ​ when Poisson effects are neglected, focusing on dominant directional stiffness.[21]

Strain-Displacement Relations

In the membrane theory of shells, the strain-displacement relations provide the kinematic foundation for linking the displacements of the shell midsurface to the resulting membrane strains, under the assumption of negligible transverse shear and bending effects. These relations are derived within the framework of differential geometry on curved surfaces, where the shell is modeled as a two-dimensional manifold embedded in three-dimensional space. The tangential displacements are denoted by the vector u=(uα)\mathbf{u} = (u_\alpha)u=(uα​) in curvilinear coordinates (α,β)(\alpha, \beta)(α,β) on the midsurface, and www represents the normal displacement. The membrane strain tensor ϵαβ\epsilon_{\alpha\beta}ϵαβ​ captures the in-plane extensions and shears, incorporating the effects of both tangential stretching and the rotation induced by normal deflection due to surface curvature.

The explicit form of the membrane strain tensor in covariant components is given by

where the semicolon denotes the covariant derivative with respect to the surface metric, bαβb_{\alpha\beta}bαβ​ is the second fundamental form (curvature tensor) describing the principal curvatures, and www is the scalar normal displacement. This expression linearizes the Green-Lagrange strain tensor for small deformations, separating the symmetric part of the tangential displacement gradient from the curvature-induced term that arises as the normal displacement alters the surface geometry. The covariant derivative ensures invariance under coordinate transformations, making the formulation applicable to arbitrary shell geometries. This standard relation originates from asymptotic expansions of three-dimensional elasticity for thin structures and is fundamental to linear shell theories.[22]

For shells of revolution, which possess axial symmetry and are commonly analyzed using meridional angle ϕ\phiϕ and azimuthal angle θ\thetaθ as coordinates, the strain-displacement relations simplify into principal directions aligned with the generators and parallels. Here, uϕu_\phiuϕ​ denotes the displacement tangent to the meridian, uθu_\thetauθ​ the circumferential displacement, RϕR_\phiRϕ​ the principal radius of curvature along the meridian, RθR_\thetaRθ​ along the parallel, and rrr the distance from the axis of revolution (parallel radius, r=Rϕsin⁡ϕr = R_\phi \sin \phir=Rϕ​sinϕ for a sphere). For the axisymmetric case (uθ=0u_\theta = 0uθ​=0, no θ\thetaθ-dependence), the meridional strain is

while the hoop strain is

These expressions account for the extension along the curving meridian and the change in circumference of the parallel circles, respectively, with the cot⁡ϕ\cot \phicotϕ term reflecting the geometry of revolution (e.g., for a sphere of radius RRR, r=Rsin⁡ϕr = R \sin \phir=Rsinϕ, Rϕ=Rθ=RR_\phi = R_\theta = RRϕ​=Rθ​=R, yielding uniform ϵϕ=ϵθ=w/R\epsilon_\phi = \epsilon_\theta = w / Rϵϕ​=ϵθ​=w/R for pure radial displacement uϕ=0u_\phi = 0uϕ​=0, w=w =w= constant). They reduce to ordinary differential equations solvable along the meridian.[23]

Under pure membrane action, where bending stiffness is neglected, an inextensibility condition may apply in certain approximations, such as for developable surfaces or pre-stressed membranes. In this context, the change in Gaussian curvature KKK of the deformed midsurface is intrinsically related to the membrane strains through the Gauss-Theorema egregium, which states that KKK is determined solely by the first fundamental form (metric tensor) of the surface. The linearized compatibility equations for the strains ϵαβ\epsilon_{\alpha\beta}ϵαβ​ ensure that the induced metric change preserves surface integrability, linking ∂2ϵαβ/∂xα∂xβ\partial^2 \epsilon_{\alpha\beta} / \partial x^\alpha \partial x^\beta∂2ϵαβ​/∂xα∂xβ (in appropriate components) to ΔK\Delta KΔK, the variation in Gaussian curvature. This geometric constraint is crucial for solvability in membrane problems without bending moments, as it enforces that the strains correspond to a realizable deformation field.

The aforementioned relations rely on the small strain approximation, valid when the gradient of displacements satisfies ∣∇u∣≪1|\nabla \mathbf{u}| \ll 1∣∇u∣≪1 and higher-order terms like (∇u)2(\nabla \mathbf{u})^2(∇u)2 or products involving w/Rw / Rw/R (where RRR is a characteristic radius of curvature) are negligible compared to unity. This linearization simplifies the exact nonlinear kinematic relations from three-dimensional continuum mechanics to a first-order theory, enabling closed-form solutions for many shell geometries while maintaining accuracy for moderately thin structures under small deflections. Constitutive relations then link these strains to membrane stress resultants, but such connections are addressed separately.[24]

Stress Resultants

In the membrane theory of shells, stress resultants represent the integrated in-plane forces per unit length acting on the shell's middle surface, obtained by integrating the three-dimensional stress components through the shell thickness while neglecting bending and shear effects. The components of the membrane stress resultant tensor, NαβN_{\alpha\beta}Nαβ​, are defined as

where σαβ\sigma_{\alpha\beta}σαβ​ denotes the in-plane stress tensor components in the curvilinear coordinates (α,β)(\alpha, \beta)(α,β) tangent to the shell surface, hhh is the shell thickness, and zzz is the coordinate normal to the middle surface. This formulation assumes that σαβ\sigma_{\alpha\beta}σαβ​ varies negligibly through the thickness, justified by the thin-shell approximation (h≪Rh \ll Rh≪R, with RRR the principal radius of curvature) and the exclusion of transverse normal stresses and moments, leading to a pure state of tangential membrane stress.[25]

In orthogonal curvilinear coordinates aligned with the principal curvature directions, the stress resultants simplify to principal components: the meridional resultant N1N_1N1​ (or NϕN_\phiNϕ​) acting along the meridian curve, and the circumferential (hoop) resultant N2N_2N2​ (or NθN_\thetaNθ​) perpendicular to it, with the shear component N12N_{12}N12​ (or NϕθN_{\phi\theta}Nϕθ​) coupling them in non-orthogonal cases. For instance, in a thin-walled cylindrical shell under uniform internal pressure ppp and radius rrr, the hoop stress resultant is Nθ=prN_\theta = p rNθ​=pr, while the longitudinal resultant is Nϕ=pr/2N_\phi = p r / 2Nϕ​=pr/2, illustrating how geometry dictates the distribution of these forces to maintain equilibrium. These principal resultants capture the load-carrying capacity in the tangential plane, with their magnitudes depending on the shell's Gaussian curvature and applied loads.[25][26]

The paths of stress resultants in shells can be visualized as "force flow" lines tracing the tangential transmission of loads across the surface, analogous to streamlines in fluid mechanics, which highlight regions of tension, compression, or shear without invoking three-dimensional stress fields. This representation underscores the membrane theory's focus on in-plane force balance, where resultants follow geodesic paths on the developable surface for minimal energy. Seminal formulations of these concepts trace to Love's integration of elasticity principles for thin shells.[25]

For shells with variable thickness h(α,β)h(\alpha, \beta)h(α,β), the integration for NαβN_{\alpha\beta}Nαβ​ is performed locally over the varying hhh, often approximated by weighting with an average thickness or assuming piecewise constancy to simplify computations; however, classical analyses typically presume uniform hhh to maintain analytical tractability. This assumption aligns with the theory's origins in Love's work and subsequent refinements by Flügge and Timoshenko for practical engineering applications.[25][26]

Shells of Revolution

Shells of revolution in membrane theory are axisymmetric structures formed by rotating a smooth curve, known as the generator or meridian, around a fixed axis in its plane, resulting in a surface of revolution. The geometry is typically parameterized using the arc length sss along the meridian or the angle ϕ\phiϕ measured from the axis of rotation, where ϕ\phiϕ increases along the generator. The principal radii of curvature are the meridional radius R1(s)R_1(s)R1​(s) (or R1(ϕ)R_1(\phi)R1​(ϕ)), which is the radius of curvature of the meridian itself, and the parallel or hoop radius R2(s)R_2(s)R2​(s) (or R2(ϕ)R_2(\phi)R2​(ϕ)), which is the distance from the axis of rotation to the point on the surface, measured perpendicular to the meridian. These radii determine the local curvature and are essential for formulating the equilibrium equations, as the shell's resistance to loads depends on its double curvature. For instance, in a spherical shell, R1=R2=RR_1 = R_2 = RR1​=R2​=R (constant), while in more general cases like ellipsoids or paraboloids, they vary along the meridian.[1]

Under axisymmetric loading conditions, which preserve rotational symmetry about the axis, membrane theory applies to predict the in-plane stress resultants without considering bending effects. Typical load types include normal pressure pn(s)p_n(s)pn​(s) distributed symmetrically over the surface, self-weight acting vertically with components resolved into normal and meridional directions, or axial forces applied along the generator; torsional or shear loads in the circumferential direction are excluded to maintain axisymmetry. The normal load pnp_npn​ is perpendicular to the shell surface, while any meridional tangent load pϕp_\phipϕ​ lies in the plane of the meridian. For self-weight of uniform density γ\gammaγ, the components are pn=−γcos⁡ϕp_n = -\gamma \cos \phipn​=−γcosϕ and pϕ=γsin⁡ϕp_\phi = \gamma \sin \phipϕ​=γsinϕ, ensuring the load vector aligns with gravity. Equilibrium is achieved solely through membrane forces, assuming the shell thickness is small compared to the radii, and no transverse shear or moments are present.[1][27]

The stress distribution in such shells is derived by integrating the equilibrium equations along the generator. The meridional stress resultant Nϕ(s)N_\phi(s)Nϕ​(s) is obtained by resolving forces in the meridional direction and integrating from the apex or a known boundary condition, often incorporating the total vertical load above a given section to balance the resultant. Specifically, for axisymmetric cases, NϕN_\phiNϕ​ satisfies

where r=R2sin⁡ϕr = R_2 \sin \phir=R2​sinϕ is the radial distance from the axis. The hoop stress resultant NθN_\thetaNθ​ is then found from the normal equilibrium equation:

yielding Nθ=R2(pn−NϕR1)N_\theta = R_2 \left( p_n - \frac{N_\phi}{R_1} \right)Nθ​=R2​(pn​−R1​Nϕ​​), which includes terms from the integrated meridional forces and local loading. This relation highlights how the hoop stress depends on the ratio of radii R1/R2R_1 / R_2R1​/R2​ and accumulated effects along the meridian. In general, NϕN_\phiNϕ​ carries the primary vertical resolution of loads, while NθN_\thetaNθ​ provides lateral stability through curvature.[1][27]

A representative example is the parabolic dome under uniform vertical loading, such as self-weight or snow approximated as uniform, where the meridian follows a parabolic profile r=kz2r = k z^2r=kz2. In this case, the meridional stress resultant NϕN_\phiNϕ​ is compressive, starting near zero at the apex and increasing in magnitude toward the base, where support reactions dominate. Such distributions are critical for design, as they indicate regions prone to buckling or requiring reinforcement, and are solved using stress functions tailored to the parabolic geometry.[25]

Cylindrical and Spherical Shells

In membrane theory, cylindrical shells under uniform internal pressure are analyzed assuming an infinite length to neglect end effects, focusing solely on in-plane stress resultants. The hoop stress resultant NθN_\thetaNθ​ is given by Nθ=prN_\theta = p rNθ​=pr, where ppp is the internal pressure, rrr is the radius, and the thickness hhh is implicit in the resultant definition as force per unit length.[28] The longitudinal stress resultant NzN_zNz​ is Nz=pr2N_z = \frac{p r}{2}Nz​=2pr​, reflecting equilibrium along the axis.[28] For finite-length cylinders, pure membrane theory neglects edge effects, providing a uniform stress distribution away from boundaries.[28]

Spherical shells under internal pressure exhibit isotropic membrane stresses due to geometric symmetry. The meridional and hoop stress resultants are equal, Nϕ=Nθ=pr2N_\phi = N_\theta = \frac{p r}{2}Nϕ​=Nθ​=2pr​, for thin-walled vessels where the radius rrr greatly exceeds thickness hhh.[29] This uniform biaxial tension arises from balancing the pressure force on a hemispherical section against the resultant forces in the wall, resulting in stresses half those of the hoop stress in a comparable cylinder.[29]

In contrast, spherical shells under uniform external pressure develop compressive membrane stresses Nϕ=Nθ=−pR2N_\phi = N_\theta = -\frac{p R}{2}Nϕ​=Nθ​=−2pR​, but membrane theory does not account for stability, and buckling instability occurs at much lower pressures. For shallow spherical caps (small angular extent α≪1\alpha \ll 1α≪1), the Donnell–Mushtari–Vlasov (DMV) theory, a simplified nonlinear shallow shell theory, is used to analyze buckling, post-buckling, thermal instability, vibrations, and static bistability (two stable states: natural and everted). The DMV theory neglects certain nonlinear terms (e.g., square of normal rotation in in-plane strains) and assumes short-wavelength deformations relative to the radius of curvature, suitable for small to moderate deflections and shallow geometries. For perfect shells, it yields the classical critical buckling pressure pcr=2E3(1−ν2)(hR)2p_{cr} = \frac{2E}{\sqrt{3(1-\nu^2)}} \left( \frac{h}{R} \right)^2pcr​=3(1−ν2)​2E​(Rh​)2, where EEE is Young's modulus, ν\nuν is Poisson's ratio, hhh is thickness, and RRR is radius.[30][31]

Practical applications of membrane theory dominate in boiler design, where cylindrical shells withstand internal steam pressure primarily through hoop and longitudinal stresses, guiding thickness selection for safety.[1] For spherical geometries, such as storage tanks or balloons, the theory predicts uniform membrane stresses under pressure, enabling efficient material use in scenarios like liquefied gas containment where bending is minimal.[1]

Analytical Solutions

In the membrane theory of shells, analytical solutions are obtained through direct integration or series expansions for cases with simple geometries and symmetric or periodic loads, providing closed-form expressions for stress resultants without bending effects.

For shells of revolution under axisymmetric loading, integration techniques exploit the meridional equilibrium equation to determine the meridional stress resultant NϕN_\phiNϕ​. Considering the free body diagram of the shell portion from the apex to arc length sss along the meridian, the resultant Nϕ(s)N_\phi(s)Nϕ​(s) balances the integrated component of the applied load in the direction parallel to the axis of revolution: Nϕ(s)=−12πrsin⁡α∫0s2πr′pz ds′N_\phi(s) = -\frac{1}{2\pi r \sin\alpha} \int_0^s 2\pi r' p_z , ds'Nϕ​(s)=−2πrsinα1​∫0s​2πr′pz​ds′, where rrr is the radius of the parallel circle at sss, α\alphaα is the angle between the meridian tangent and the axis, and pzp_zpz​ is the axial load component per unit area. The hoop stress resultant NθN_\thetaNθ​ is then found from the normal equilibrium equation: Nϕr1+Nθr2=pn\frac{N_\phi}{r_1} + \frac{N_\theta}{r_2} = p_nr1​Nϕ​​+r2​Nθ​​=pn​, with r1r_1r1​ and r2r_2r2​ the principal radii of curvature and pnp_npn​ the normal load component. This "marching" approach along the meridian yields exact solutions for geometries defined parametrically, such as spheres or cones, under uniform pressure or self-weight.[](Timoshenko and Woinowsky-Krieger 1959)[](Fl%C3%BCgge 1967)

Exact solutions are particularly notable for hyperbolic shells under self-weight, where the geometry approximates an inverted catenary to achieve pure compression. For such a shell with line load www due to self-weight, the meridional force simplifies to Nϕ=−wr/tan⁡αN_\phi = -w r / \tan\alphaNϕ​=−wr/tanα, with rrr the parallel radius and α\alphaα the angle of the tangent to the meridian with the horizontal; the hoop force NθN_\thetaNθ​ follows from normal balance, ensuring no bending moments throughout. This form arises from the equilibrium of vertical forces, where the shell shape satisfies the differential equation for uniform thickness under gravity, as applied in cooling tower design.[](Zingoni 1997)

For cylindrical shells subjected to non-axisymmetric loads, series methods decompose the loading into a Fourier series in the circumferential coordinate θ\thetaθ: p( x, θ)=∑n=0∞(ancos⁡nθ+bnsin⁡nθ)eikxp(\ x,\ \theta) = \sum_{n=0}^\infty (a_n \cos n\theta + b_n \sin n\theta) e^{ikx}p( x, θ)=∑n=0∞​(an​cosnθ+bn​sinnθ)eikx or similar modal forms along the axial coordinate xxx. Substituting into the governing partial differential equation reduces the problem to a set of ordinary differential equations for each mode nnn, solvable analytically using exponential or hyperbolic functions for infinite cylinders, or boundary matching for finite lengths. These solutions capture circumferential variations efficiently for wind or thermal loads.[](Timoshenko and Woinowsky-Krieger 1959)

Boundary value problems in membrane theory, such as specified NϕN_\phiNϕ​ or displacements at shell edges, are addressed via direct integration for axisymmetric cases or the method of characteristics for the resulting hyperbolic partial differential equations in general geometries. For instance, in shells of revolution, edge conditions like zero radial displacement lead to characteristic curves along which Riemann invariants propagate, allowing step-by-step solution from boundaries inward; compatibility with constitutive relations ensures unique stress fields.[](Niordson 1984)

Numerical Approaches

The finite element method (FEM) serves as a primary numerical approach for analyzing membrane shells, particularly for complex geometries where analytical solutions are infeasible. Membrane elements in FEM are typically formulated as thin surfaces capable of carrying in-plane stresses without bending stiffness, often using 3-node triangular or 4- to 6-node quadrilateral elements to discretize the shell midsurface. The element stiffness matrix is derived from the integral K=∫ABTCB dA\mathbf{K} = \int_A \mathbf{B}^T \mathbf{C} \mathbf{B} , dAK=∫A​BTCBdA, where B\mathbf{B}B is the strain-displacement matrix incorporating membrane strains, C\mathbf{C}C is the constitutive matrix relating stresses to strains, and the integration is performed over the element area AAA.[32][33]

Discrete formulations provide an alternative numerical framework for membrane shell equilibrium, leveraging concepts from discrete differential geometry to model integrable discrete shell membranes. These approaches often employ bar networks or quadrilateral "plated" elements inscribed in circles, ensuring equilibrium through discrete analogues of continuous integrability conditions derived from O-surface theory. Such methods are particularly useful for simulating large deformations while preserving geometric constraints inherent to membrane theory.[34]

In practical implementations, commercial software like ANSYS and ABAQUS facilitates membrane shell analysis using dedicated membrane-only elements, which support arbitrary shapes through surface meshing and handle nonlinear effects such as geometric nonlinearity. For instance, ANSYS shell elements can be configured for pure membrane behavior by suppressing bending and shear terms, while ABAQUS membrane elements explicitly model in-plane forces without transverse stiffness. These tools enable efficient computation for engineering applications like pressure vessels or tensile structures.[35][36]

Asymptotic expansions offer a numerical validation strategy for membrane approximations, particularly in semi-membrane theories that blend pure membrane behavior with minor bending effects through series expansions in the small parameter h/Rh/Rh/R, where hhh is shell thickness and RRR is the principal radius of curvature. These methods systematically improve accuracy for moderately thick shells by integrating three-dimensional elasticity equations asymptotically.[37][38]

Validity and Limitations

The membrane theory of shells is valid under conditions where bending deformations are negligible compared to in-plane stretching, primarily for thin shells with a small thickness-to-radius ratio, such as h/R<1/20h/R < 1/20h/R<1/20 (or R/h>20R/h > 20R/h>20), allowing the neglect of moment resultants while retaining extensional stresses.[39] This approximation holds particularly well for uniform or smoothly distributed loads where the characteristic wavelength of loading significantly exceeds the shell radius RRR, as well as for shallow shells under axisymmetric conditions, ensuring equilibrium through hoop and meridional forces without significant curvature changes.[23] For instance, in cylindrical or spherical shells subjected to internal pressure, the theory accurately predicts stress resultants, aligning closely with experimental measurements for thin-walled pressure vessels where radial expansions dominate.[39]

However, the theory fails near edges, supports, and discontinuities, where incompatible boundary conditions induce localized bending moments in boundary layers, typically spanning 10-20% of the shell height and decaying exponentially away from the edge.[23] It also breaks down under compressive loads, as thin shells exhibit buckling instability before reaching significant membrane compression, limiting its use in stability analyses without supplementary bending considerations.[40] In regions dominated by bending, such as near concentrated loads or abrupt geometric changes, the theory overpredicts structural stiffness by ignoring transverse shear and rotary inertia effects.[39]

Further limitations arise for thicker shells where h/R>1/10h/R > 1/10h/R>1/10, as bending contributions become prominent, invalidating the pure membrane state assumption; similarly, dynamic or non-uniform loads, like those from wind on roofs, introduce local bending that causes deviations from predictions.[39] While extensions incorporating bending terms can address these shortcomings, membrane theory remains a foundational approximation best suited to extensional-dominated responses in thin, smoothly loaded structures.[23]

Extensions to Bending Theories

The bending theory of shells extends the membrane theory by incorporating resistance to transverse shear deformation and flexural rigidity, which are crucial for capturing localized effects such as edge bending and boundary disturbances. In this framework, moment resultants MαβM_{\alpha\beta}Mαβ​ (bending and twisting moments) and Kirchhoff shear forces QαQ_\alphaQα​ (effective transverse shears derived from moment equilibrium) are introduced alongside the membrane stress resultants NαβN_{\alpha\beta}Nαβ​. These additions lead to a system of fifth-order partial differential equations (PDEs) governing the shell's deformation, in contrast to the second-order PDEs of pure membrane theory, enabling the analysis of bending-dominated behaviors in thin shells under non-uniform loading or at boundaries.[41]

Semi-membrane approaches represent hybrid models that blend membrane assumptions with selective inclusion of bending stiffness, particularly suited for structures like long cylindrical shells where global behavior is membrane-like but local ring bending occurs at the ends or under concentrated loads. A prominent example is the Donnell-Mushtari-Vlasov (DMV) theory, a simplified nonlinear shallow shell theory that approximates the shell equations by neglecting certain nonlinear terms (e.g., the square of the normal rotation in the in-plane strains) and assuming short-wavelength deformations relative to the radius of curvature. This makes it suitable for small to moderate deflections and shallow geometries (angular extent α ≪ 1), including thin spherical caps (shallow portions of spheres). The theory derives coupled fourth-order PDEs for the normal displacement www and Airy stress function FFF (often denoted φ), typically in cylindrical coordinates. These equations are widely used for buckling, post-buckling, thermal instability, vibrations, and static bistability (two stable states: natural and everted) of spherical caps, as well as buckling and vibration analyses of cylindrical and nearly flat shells, providing a computationally efficient approximation valid when the shell's rise is small compared to its planform dimensions. For perfect spherical shells, it yields the classical critical buckling pressure under uniform external pressure: pcr=2E3(1−ν2)(tR)2p_{\mathrm{cr}} = 2E \sqrt{3(1-\nu^{2})} \left( \frac{t}{R} \right)^{2}pcr​=2E3(1−ν2)​(Rt​)2, where EEE is Young's modulus, ν\nuν is Poisson's ratio, ttt is thickness, and RRR is radius.[42][30]

More rigorous linearized extensions, such as the Koiter and Sanders theories, build on the Kirchhoff-Love hypotheses to formulate bending-inclusive models for shallow and moderately curved shells, emphasizing variational principles for stability and small deformations. Koiter's theory derives a consistent first-order approximation through asymptotic expansion of the three-dimensional elasticity equations, yielding energy functionals that combine quadratic membrane strain energy with bending energy terms proportional to the square of curvature changes, ∫(ϵαβϵαβ+καβκαβ) dA\int ( \epsilon_{\alpha\beta} \epsilon^{\alpha\beta} + \kappa_{\alpha\beta} \kappa^{\alpha\beta} ) , dA∫(ϵαβ​ϵαβ+καβ​καβ)dA, where ϵ\epsilonϵ and κ\kappaκ denote membrane and bending strains, respectively; this approach is particularly effective for post-buckling analysis and error estimation in thin shells. Sanders' nonlinear extension refines this by incorporating moderate rotations and improved strain-displacement relations, avoiding inconsistencies in earlier approximations and facilitating applications to imperfection-sensitive buckling in cylindrical shells. Both theories are valid for shallow geometries where the Gaussian curvature is small, offering enhanced accuracy over DMV for nonlinear regimes without the full complexity of sixth-order shell equations.[43][44]

In transitioning from membrane to bending theories, the membrane solution serves as a zeroth-order approximation for the interior of the shell, with bending perturbations introduced asymptotically to resolve boundary layers—narrow regions near edges where rapid variations in deflection and stresses occur due to compatibility requirements. The characteristic width of these boundary layers scales as Rh\sqrt{R h}Rh​, where RRR is the principal radius of curvature and hhh is the shell thickness, allowing membrane theory to approximate the bulk while bending theory corrects localized effects; this perturbation method, rooted in matched asymptotic expansions, underpins many semi-membrane and shallow shell analyses for efficient yet accurate predictions.[41][45]

Foundational Texts

The foundational texts of membrane theory for shells emerged in the late 19th and early 20th centuries, building on principles of elasticity to describe stress states in thin curved structures where bending moments are negligible. Early contributions focused on the equilibrium of flexible, inextensible surfaces, establishing the core equations for membrane forces. Eugenio Beltrami's 1882 work, Sull’equilibrio delle superficiali flessibili ed inestendibili, provided one of the first mathematical treatments of equilibrium in curved elastic surfaces, deriving equations that relate membrane stresses to surface geometry and applied loads without considering flexural rigidity.[2] This paper laid groundwork for analyzing thin shells under uniform tension, emphasizing in-plane forces that maintain structural integrity in momentless states.[2]

Complementing Beltrami, Léon Lecornu's 1885 publication, Sur l'équilibre des surfaces flexibles et inextendibles in the Journal de l'École Polytechnique, extended these ideas to thin shells by formulating equilibrium conditions for flexible surfaces under distributed loads.[46] Lecornu derived differential equations governing membrane stresses, assuming negligible thickness and no resistance to bending, which proved essential for shells like domes and cylinders where membrane action dominates.[2] Together, these papers represented pioneering efforts in membrane theory, prioritizing simplified stress analysis for engineering applications over full three-dimensional elasticity.[47]

Augustus Edward Hough Love's A Treatise on the Mathematical Theory of Elasticity (1892) marked the first systematic integration of membrane theory into a broader elastic framework for shells. Love built on prior elasticity works by Cauchy and Kirchhoff, introducing approximations for thin shells where normal stresses are ignored and midsurface strains govern behavior.[47] His derivations yielded general equations for membrane forces in curved surfaces, applicable to arbitrary geometries, and established the Love-Kirchhoff hypotheses that remain central to the field.[47] This text shifted focus toward practical computation of deformations in elastic membranes, influencing subsequent shell analyses.

Stephen Timoshenko and S. Woinowsky-Krieger's Theory of Plates and Shells (1940) offered practical formulations tailored for engineers, emphasizing membrane solutions for common geometries like cylinders and spheres.[41] The book derived explicit equations for axisymmetric loading, such as hoop and meridional stresses in spherical shells under internal pressure, assuming pure membrane action without bending.[41] Timoshenko's approach simplified Love's general theory into solvable forms, incorporating boundary conditions for real-world structures like pressure vessels, and highlighted the theory's limitations near edges where bending effects emerge.[47]

Vasily Z. Vlasov's General Theory of Shells and Its Application in Engineering (1958, translated into English in 1961 by NASA) provided a comprehensive set of membrane equations for general surfaces, advancing beyond earlier works by addressing non-developable geometries. Vlasov formulated coupled differential equations relating principal curvatures to stress resultants, enabling analysis of complex shells like hyperboloids, while maintaining the membrane assumption of zero transverse shear.[48] This text synthesized prior contributions into a unified framework, emphasizing applications in aerospace and civil engineering, and introduced semi-momentless approximations for efficient computation.

Modern Reviews

In the late 20th and early 21st centuries, mathematical analyses have provided rigorous justifications for membrane shell models, particularly through asymptotic methods and variational convergence techniques. Philippe G. Ciarlet's "Mathematical Elasticity" series, culminating in Volume III: Theory of Shells (2000), establishes the foundations of nonlinear elastic membrane shells using Gamma-convergence to derive two-dimensional equations from three-dimensional elasticity, emphasizing the limiting behavior as shell thickness approaches zero.[49] This work highlights how membrane approximations capture in-plane stresses effectively for thin structures while neglecting bending, offering a mathematically sound basis for applications in engineering design.

Frithiof I. Niordson's "Shell Theory" (1985) provides a comprehensive overview distinguishing membrane theory from bending-inclusive models, incorporating numerical examples to illustrate stress distributions in cylindrical and spherical shells under various loads.[50] The text underscores the practical utility of membrane assumptions for large-radius shells, where bending effects are minimal, and includes computational strategies for solving equilibrium equations.

Recent advancements extend membrane theory into discrete and computational realms. For instance, W. K. Schief's 2014 paper on integrable structure in discrete shell membrane theory introduces discrete analogues of classical integrable membrane classes, enabling equilibrium-preserving discretizations suitable for computational design of curved structures like architectural shells. Complementing this, Ciarlet and Véronique Lods' 1996 asymptotic analysis justifies membrane shell equations for linearly elastic cases, deriving Koiter-type models via expansion in thickness and proving their validity under ellipticity conditions.

Engineering handbooks continue to adapt membrane theory for industrial applications. Wilhelm Flügge's "Stresses in Shells" (1967, with later reprints) remains influential for practical calculations of membrane stresses in pressure vessels and roofs, providing approximate formulas and charts for axisymmetric cases while noting limitations near boundaries. These modern reviews collectively synthesize theoretical rigor with applied extensions, bridging classical foundations to contemporary computational and design challenges in shell mechanics.