We investigate the influence of curvature and topology on crystalline wrinkling patterns in generic elastic bilayers. Our numerical analysis predicts that the total number of defects created by adiabatic compression exhibits universal quadratic scali
ng for spherical, ellipsoidal and toroidal surfaces over a wide range of system sizes. However, both the localization of individual defects and the orientation of defect chains depend strongly on the local Gaussian curvature and its gradients across a surface. Our results imply that curvature and topology can be utilized to pattern defects in elastic materials, thus promising improved control over hierarchical bending, buckling or folding processes. Generally, this study suggests that bilayer systems provide an inexpensive yet valuable experimental test-bed for exploring the effects of geometrically induced forces on assemblies of topological charges.
A thin shell finite element approach based on Loops 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 arbitr
ary 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.
A finite element program is presented to simulate the process of packing and coiling elastic wires in two- and three-dimensional confining cavities. The wire is represented by third order beam elements and embedded into a corotational formulation to
capture the geometric nonlinearity resulting from large rotations and deformations. The hyperbolic equations of motion are integrated in time using two different integration methods from the Newmark family: an implicit iterative Newton-Raphson line search solver, and an explicit predictor-corrector scheme, both with adaptive time stepping. These two approaches reveal fundamentally different suitability for the problem of strongly self-interacting bodies found in densely packed cavities. Generalizing the spherical confinement symmetry investigated in recent studies, the packing of a wire in hard ellipsoidal cavities is simulated in the frictionless elastic limit. Evidence is given that packings in oblate spheroids and scalene ellipsoids are energetically preferred to spheres.
We investigate the morphology of thin discs and rings growing in the circumferential direction. Recent analytical results suggest that this growth produces symmetric excess cones (e-cones). We study the stability of such solutions considering self-co
ntact and bending stress. We show that, contrary to what was assumed in previous analytical solutions, beyond a critical growth factor, no symmetric textit{e}-cone solution is energetically minimal any more. Instead, we obtain skewed e-cone solutions having lower energy, characterized by a skewness angle and repetitive spiral winding with increasing growth. These results are generalized to discs with varying thickness and rings with holes of different radii.
