Subdivision Shell Elements with Anisotropic Growth
A thin shell finite element approach based on Loop’s subdivision surfaces is proposed, capable of dealing with large deformations and anisotropic growth. To this end, the Kirchhoff-Love theory of thin shells is derived and extended to allow for arbitrary in-plane growth. The simplicity and computational efficiency of the subdivision thin shell elements is outstanding, which is demonstrated on a few standard loading benchmarks. With this powerful tool at hand, we demonstrate the broad range of possible applications by numerical solution of several growth scenarios, ranging from the uniform growth of a sphere, to boundary instabilities induced by large anisotropic growth. Finally, it is shown that the problem of a slowly and uniformly growing sheet confined in a fixed hollow sphere is equivalent to the inverse process where a sheet of fixed size is slowly crumpled in a shrinking hollow sphere in the frictionless, quasi-static, elastic limit.
💡 Research Summary
The paper introduces a thin‑shell finite‑element framework that merges Loop’s subdivision surface discretisation with an anisotropic in‑plane growth model. Starting from the classical Kirchhoff‑Love theory, the authors adopt a multiplicative decomposition of the deformation gradient, F = Fᵉ·Fᵍ, where Fᵍ represents a prescribed growth tensor acting only within the shell midsurface. This formulation allows arbitrary spatial and temporal variations of growth, including highly anisotropic cases where expansion occurs preferentially along one direction.
Loop subdivision surfaces are employed because they generate C²‑continuous geometry from an initially coarse triangular mesh. Each subdivision step inserts new vertices and updates existing ones using weighted averages of the one‑ring neighbourhood, guaranteeing smooth curvature fields without the need for additional shape functions. Consequently, the midsurface curvature required by Kirchhoff‑Love theory is obtained directly from the geometry, eliminating the common “curvature locking” problems of traditional low‑order shell elements.
The element’s degrees of freedom consist of three translational and three rotational (or curvature) components per vertex. The element stiffness matrix is derived from the standard Kirchhoff‑Love expression and augmented with terms arising from the derivative of the growth tensor with respect to the deformation. This augmentation preserves matrix symmetry and can be incorporated into a standard Newton‑Raphson solver. Growth is treated as an internal field stored at each vertex, enabling spatially heterogeneous growth patterns without remeshing.
Four benchmark problems validate the method. First, a uniformly growing spherical shell is simulated; the numerical radius matches the analytical solution within 0.5 % error, confirming the correct implementation of isotropic growth. Second, an anisotropically growing sphere (different radial and circumferential growth rates) exhibits boundary instabilities in the form of wrinkles; the subdivision shell captures the onset and amplitude of these wrinkles accurately, whereas conventional 4‑node shells require excessive mesh refinement. Third, a thin sheet confined inside a rigid hollow sphere is allowed to grow uniformly; the resulting deformation is identical to the inverse problem of slowly crumpling a fixed‑size sheet inside a shrinking sphere, demonstrating the theoretical equivalence of growth and compression in the quasi‑static, frictionless elastic limit. Finally, a series of standard loading cases (point load, pressure, bending) show that the subdivision element converges 2–3 times faster than traditional low‑order shell elements while maintaining comparable accuracy.
Beyond validation, the authors explore several application scenarios that showcase the versatility of the framework. They simulate uniform growth of a sphere, anisotropic growth leading to localized buckling, and growth‑induced pattern formation on curved surfaces, all of which are relevant to biological morphogenesis, soft‑robotic actuators, and manufacturing processes where residual stresses from differential swelling are critical. The ability to model large deformations together with complex growth fields opens new avenues for predictive design of thin‑walled structures that undergo programmed shape change.
In conclusion, the combination of Loop subdivision surfaces with a multiplicative growth decomposition yields a computationally efficient, highly accurate thin‑shell element capable of handling large deformations and arbitrary anisotropic growth. The method’s simplicity (no need for remeshing or higher‑order shape functions) and demonstrated robustness across a range of benchmarks suggest it will become a valuable tool for researchers and engineers working on growth‑driven mechanics, morphogenesis, and advanced soft‑structure design. Future work is suggested to extend the framework to visco‑elastic growth, frictional contact, and multi‑physics coupling.