PaperlayerView exact revision on arXiv

Source-linked AI summary

Inverse origami inspired by insect morphogenesis: geometric design of wrinkle structures for curved surface formation

Kentaro Morikawa, Yuuki Shimizu, Nen Saito, Yasuhiro Inoue

arXiv:2610.34903v1cond-mat.softcs.GRphysics.app-phphysics.bio-ph

TL;DR

Deriving a wrinkle structure from a prescribed target surface is the central design problem in forming curved surfaces through unfolding. This paper models wrinkle unfolding as a unidirectional stretching field and derives wrinkle geometry from target geometry; its designs are nearly isometric to targets, with area-weighted median strain-like error of 0.2–1.7% after amplitude optimization.

  • Problem

    The central design problem is deriving a wrinkle structure that encodes a prescribed target surface.

    Source: Discussion · p. 7

  • Method

    The method represents wrinkle unfolding as a unidirectional stretching field and derives wrinkle geometry from the target surface.

    Source: The relationship between wrinkle… · p. 2, Discussion · p. 7

  • Results

    0.2–1.7% area-weighted median strain-like error after amplitude optimization, reducing median error to 32–67% of pre-optimization values across five target shapes.

    Source: Demonstrations · p. 7

  • Takeaways & Limitations

    The principle provides a common framework for designing deployable structures and programmable materials and generating testable hypotheses about insect morphogenesis.

    Source: Discussion · p. 10

    The framework is restricted to wrinkles represented as height functions over a plane, limiting its scope on general curved substrates.

    Source: Discussion · p. 7

Abstract

from arXiv · show

Insect exoskeletons abruptly acquire complex three-dimensional forms during metamorphosis from larva to pupa. This transformation occurs through the unfolding of preformed wrinkle structures inside the larval body under internal pressure. Inspired by this process, we investigate an inverse-origami strategy for fabricating curved surfaces. Whereas conventional origami produces three-dimensional objects by folding a flat sheet, this strategy produces curved surfaces by unfolding preformed wrinkle structures. The central design problem is to derive a wrinkle structure from a prescribed target surface. Our key idea is to regard a wrinkle structure as a unidirectional stretching field, which specifies at each point the stretch produced when the structure unfolds. We first derive such a field from the target surface geometry and then translate it into spatially varying wrinkle directions and spacing, so that the wrinkle structure unfolds into the target. The resulting numerical designs approximate the intrinsic metric of the target surface. We demonstrate the method on multiple target surfaces through numerical construction, mechanical unfolding simulations, and 3D-printed prototypes. This work provides a biologically inspired geometric design principle for encoding prescribed surface shapes in wrinkle structures, as well as a framework for designing and fabricating deployable curved surfaces.

Inverse origami principle

The inverse-origami principle treats wrinkle unfolding as a unidirectional stretch field and encodes a target surface by deriving that field from a suitable parameterization. It translates the field into a wrinkle pattern whose induced metric matches the target under stated compatibility conditions.

The relationship between wrinkle structures and surfaces.

  • On the plane, a wrinkle structure is a unidirectional stretch field: stretch is perpendicular to ridges, and its magnitude tracks wrinkle fineness under constant amplitude.In wave terminology, stretch direction aligns with the wave-number vector, whose magnitude corresponds to wrinkle fineness.Source: The relationship between wrinkle… · pp. 2–3
  • A piecewise-linear periodic function enables sharp creases, allowing isometric folding without the curvature change associated with smooth bends; smoothly varying phase yields smooth wrinkle regions between creases.For the hemisphere example, a smooth bend creates a region of negative Gaussian curvature, whereas sharp creases preserve curvature away from crease curves.Source: The relationship between wrinkle… · p. 3

Encoding target surfaces into wrinkle structures.

  • The target is parameterized so the Jacobian singular values satisfy σ1 ≥1 and σ2 =1; the principal stretch yields s, then k, phase ϕ, and wrinkle height w.The stretch direction is the principal direction associated with σ1, and its stretch ratio is σ1.Source: Encoding target surfaces into wrinkle… · p. 4

Verification of isometry via the induced metric.

  • When a suitable parameterization and globally compatible phase exist, the constant-amplitude wrinkle surface has the target’s induced metric wherever f is differentiable, with piecewise-smooth isometry across creases.The parameterization must satisfy σ1 ≥1 and σ2 =1; phase compatibility includes the necessary conditions around singularities and non-simply-connected regions.Source: Verification of isometry via the… · pp. 4–5

Numerical construction of wrinkle structures for general surfaces

The method converts a target surface into a unidirectional stretching field, then constructs a wrinkle structure from that field. A postprocessing amplitude adjustment alleviates height discontinuities at singular points and reduces deviation from isometry across most of every tested surface.

Numerical construction of wrinkle structures for general surfaces

  • A local/global parameterization constrains each triangle to a unidirectional stretch, yielding the stretching field used to encode the target surface.The constraint sets the larger singular value to at least 1 and the smaller to 1.Source: Numerical construction of wrinkle… · p. 5
  • The construction converts the stretching field into a wave-number field and uses a complex phase function to handle non-integrable regions.It minimizes the phase-field residual while allowing zeros that introduce singularities, then derives wrinkle heights from the phase and remeshes to place crease lines on mesh edges.Source: Numerical construction of wrinkle… · p. 5

Amplitude optimization to resolve discontinuities at singular points.

  • Singular points in the wave-number field can make the phase and wrinkle height discontinuous, motivating position-dependent amplitude adjustment.Such points occur where wrinkle ridges terminate; the adjustment relaxes the constant-amplitude assumption used to derive the wave-number field.Source: Amplitude optimization to resolve… · p. 5
  • Amplitude optimization reduces deviation from isometry over most of the surface area for every target shape.The adjustment also alleviates height discontinuities at phase singularities by updating vertex heights and allowing position-dependent wrinkle amplitude.Source: Amplitude optimization to resolve… · pp. 5–6

Evaluation methods for the designed wrinkle structures

The wrinkle designs are evaluated for isometry, mechanical unfoldability, agreement with the target surface, and physical implementation through 3D-printed prototypes.

Evaluation of isometry.

  • Isometry is evaluated per triangle using the Frobenius norm between its Jacobian and the closest rigid rotation; area-weighted median and 95th-percentile strain-like errors summarize deviations.The strain-like error is the largest deviation of principal stretch ratios from unity, reported as percentage strain; area weighting reduces the influence of near-degenerate triangles at phase singularities.Source: Evaluation of isometry · p. 6

Unfolding simulation of the wrinkle structure.

  • Unfolding is simulated by minimizing a discrete shell energy under internal pressure, with area and angle terms maintaining isometry and bending driving wrinkle deployment.Pressure acts only during the first n1 steps, then is switched off; it inflates the structure toward a target lumen volume.Source: Unfolding simulation of the wrinkle… · p. 6

Evaluation of agreement between the unfolded shape and the target.

  • Agreement with the target is measured after rigid ICP registration without rescaling, using bidirectional unsigned point-to-surface distances normalized by the target bounding-box diagonal.Source: Evaluation of agreement between the… · p. 6
  • All five target shapes are physically implemented by multimaterial 3D printing, with creases thinned as flexible hinges and regions requiring inversion printed in soft material.The remaining regions are printed in stiff material, and the structures are unfolded manually.Source: Evaluation of agreement between the… · p. 6

Demonstrations

The demonstrations show that inverse-origami wrinkle structures can encode varied target surfaces and unfold into shapes that closely resemble them, both in simulation and physical prototypes.

Demonstrations

  • The demonstrations covered cap, hat, overhang, horn-like, and frill surfaces; the workflow derived a stretching field, computed wrinkle phase and height, then optimized amplitude and simulated unfolding.The targets span dome-like forms, an overhanging region, a biologically motivated horn-like shape, and a wavy negatively curved rim.Source: Demonstrations · p. 6
  • After amplitude optimization, area-weighted median strain-like error ε was 0.2–1.7% (32–67% of its initial value), while unfolded-shape mean distance was 0.06–2.1% of the target bounding-box diagonal.Residual strain-like error concentrated near phase singularities; simulations used coarsened meshes, so distances include remeshing error. Frill crest placement is not determined by the intrinsic metric.Source: Demonstrations · p. 7
  • 3D-printed structures designed at A = 0.12 reproduced the characteristic features of all five targets after manual unfolding, but the frill reverted toward its folded state upon release.The designs' metric error at A = 0.12 was comparable to that at A = 0.02.Source: Demonstrations · p. 7

Discussion

The discussion presents inverse origami as a geometric method for deriving preformed wrinkle structures from target surfaces, while qualifying its exactness, biological interpretation, and physical realization. It also identifies extensions needed for broader geometric and engineering applications.

Discussion

  • The stretching field specifies geometric information, not developmental mechanisms, and should not be identified directly with cell division, tissue growth, or stress.At the beetle-horn tip, division orientation is thought to play a minor role, while spatial variation in division frequency and areal change may influence shape (5, 25). The horn-like design only qualitatively resembles the primordium (4, 28); testing that hypothesis requires measured shapes and quantitative comparisons.Source: Discussion · p. 7
  • Exactness requires a target parameterization with singular values σ1 ≥ 1 and σ2 = 1, plus a globally compatible phase field; otherwise, the construction is approximate.For constant amplitude, the constructed surface has the target’s induced metric where differentiable, with piecewise-smooth isometry along creases; amplitude optimization modifies the construction, especially near singularities.Source: Discussion · p. 7
  • The framework represents wrinkles as height functions over a plane, limiting its applicability to nearly planar substrates; wrinkle direction also becomes ill-conditioned as stretch approaches one.Extending the method to curved substrates while consistently treating singularities, variable amplitude, and field regularity remains a challenge.Source: Discussion · p. 7
  • For the frill, the metric does not uniquely determine crest placement or number, and unfolding simulations assess mechanical feasibility and visualize deployment rather than establish isometry.Bending stiffness, boundary conditions, and discretization affect the crests and troughs; metric error is the more appropriate measure of isometry.Source: Discussion · pp. 7–8
  • Unlike conventional origami, the method derives wrinkle geometry from the target; unlike pop-up printing, it derives the crease pattern rather than requiring the user to prescribe it.The discussion positions the principle as a shared framework for deployable structures, programmable materials, and testable hypotheses about insect morphogenesis.Source: Discussion · pp. 9–10
  • The method derives wrinkle geometry, including crease patterns, from a prescribed target surface by converting a unidirectional stretching field into a preformed wrinkle surface.Under compatibility conditions, the wrinkle surface is isometric to the target, and unfolding realizes the effective stretch geometrically without material extension.Source: Discussion · pp. 7, 9
  • The 3D-printed structures demonstrate fabrication and unfolding, but do not establish physical isometry or general stability; the frill reverted toward its folded state when released.Stiff material was used in regions not intended to deform to resist distortion from finite-thickness crease inversion. Robust self-deployment requires joint design of wrinkle geometry, stiffness, thickness, boundary constraints, and loading.Source: Discussion · pp. 9–10
  • The numerical pipeline was implemented in custom C++ code, and unfolding simulations used simulated annealing on meshes coarsened to 104 triangles.The passage directs readers to the SI Appendix for stage-specific energy functions, update rules, mesh resolutions, iteration counts, and simulation parameters.Source: Discussion · p. 10

SI Appendix

The SI Appendix provides supporting derivations, implementation and evaluation details, and 3D-printing parameters summarized in the main text.

Supporting derivations and methods

  • The appendix derives the isometry condition and details the numerical construction pipeline and unfolding-simulation energy functional.Source: SI Appendix · p. 11

Evaluation and fabrication

  • It also specifies the quantitative evaluation procedures and 3D-printing parameters.Source: SI Appendix · p. 11

Numbering conventions

  • Equation, figure, and table numbers prefixed with “S” refer to the appendix, while unprefixed equation numbers refer to the main text.Source: SI Appendix · p. 11

Derivation of the isometry condition

The isometry condition is derived by matching the target surface’s pullback metric to the metric induced by the wrinkle structure. Equating them yields Eq. 9 of the main text.

Derivation of the isometry condition

  • The target metric is pulled back through its parameterization and expressed using coefficients E, F, and G, with the parameterization constrained by σ2 = 1.The resulting metric depends on the stretching σ1 and wrinkle-direction angle α.Source: Derivation of the isometry condition · p. 11
  • For the wrinkle structure φW(u, v) = (u, v, w(u, v)), the induced metric is du2 + dv2 + dw2, with w(u, v) = A f(ϕ(u, v)).Source: Derivation of the isometry condition · p. 11
  • Equating the target and wrinkle-structure metrics gives the isometry condition, which is Eq. 9 of the main text.The matching double signs in Eqs. (S5) and (S6) are required.Source: Derivation of the isometry condition · p. 11

Numerical construction of the wrinkle structure · Parameterization for the unidirectional stretching field.

The wrinkle parameterization is computed by an alternating local/global optimization constrained to allowed triangle transformations. Its resulting triangle transformations yield the unidirectional stretching field.

Parameterization for the unidirectional stretching field.

  • The method minimizes an area-weighted objective over vertex positions and per-triangle transformations constrained to the allowed set M.The parameterization φ maps the domain D to the surface S, with each triangle’s Jacobian constrained by Eq. 10.Source: Parameterization for the… · p. 11
  • The local/global algorithm alternates between fixing vertex positions to find triangle transformations and fixing transformations to find vertex positions.It extends as-rigid-as-possible surface modeling.Source: Parameterization for the… · p. 11
  • The initial vertex positions are obtained by conformal parameterization.Source: Parameterization for the… · p. 11
  • In the local phase, the Jacobian’s singular value decomposition is used to determine the triangle transformation update.Source: Parameterization for the… · p. 11
  • In the global phase, the method solves for updated vertex positions with the triangle transformations fixed.The supplied passage also gives a vertex-indexed condition equal to zero, but does not expose further equation details.Source: Parameterization for the… · p. 12
  • For each triangle, the resulting transformation yields the unidirectional stretch s_i = σ1,i v1,i.The parameterization is obtained by repeating the local and global phases.Source: Parameterization for the… · p. 12

Constructing the phase field by the stripe pattern method.

The method represents phase as the argument of a complex function because the wave-number field is generally not curl-free. It finds the phase through a constrained eigenvalue problem, then derives vertex heights and remeshes crease lines onto mesh edges.

Constructing the phase field by the stripe pattern method.

  • Because the wave-number field is generally not curl-free, the stripe pattern method represents phase as the argument of a complex function ψ and minimizes an edge-based energy.The edge phase changes use the average wave-number vectors at the edge endpoints.Source: Constructing the phase field by the… · p. 12
  • The energy is a Hermitian quadratic form; imposing ψ∗Bψ = 1 makes its minimization equivalent to a generalized eigenvalue problem.The constraint avoids the trivial solution ψ = 0.Source: Constructing the phase field by the… · p. 12
  • The smallest-eigenvalue eigenvector gives vertex phases ϕi = arg ψi, which determine heights; the mesh is remeshed so crease lines lie on its edges.The eigenvector is computed by the inverse power method.Source: Constructing the phase field by the… · p. 12

Amplitude optimization at singular points.

The section describes an as-rigid-as-possible adjustment of wrinkle heights that uses rigid rotations to preserve isometry to the target surface and alleviate height discontinuities at singular points.

Amplitude optimization at singular points

  • The amplitude adjustment optimizes vertex heights while restricting each triangle’s deformation to rigid rotations, preserving isometry to the target surface.It uses the same as-rigid-as-possible framework as the parameterization.Source: Amplitude optimization at singular… · p. 12
  • The height update alternates local and global phases, with the local phase using singular value decompositions of triangle Jacobians.The global phase updates the vertex heights.Source: Amplitude optimization at singular… · p. 12
  • The iteration yields position-dependent wrinkle amplitude A(u, v) = w(u, v)/f(ϕ(u, v)), alleviating height-function discontinuities at singular points.Source: Amplitude optimization at singular… · p. 12

Energy functional for the unfolding simulation … Physical implementation by 3D printing

The paper simulates pressure-driven wrinkle unfolding by minimizing an energy that balances isometry, bending, and internal pressure, then evaluates the resulting structures and fabricates all five target shapes by 3D printing. Isometry errors are reported as comparable across the two tested amplitudes, before and after amplitude optimization.

Energy functional for the unfolding simulation

  • The unfolding simulation minimizes an energy by simulated annealing, using area and angle deviations to maintain isometry and dihedral deviation as bending energy to unfold wrinkles.Internal-pressure energy drives the lumen volume toward a target volume, but acts only during the first n1 steps to exclude pressure-induced stretching.Source: Energy functional for the unfolding… · pp. 12–13

Numerical parameters for the demonstrations

  • Demonstration meshes were reduced to 10,000 triangles; simulations used fixed-boundary, pressurized steps followed by free-boundary steps without pressure, with a different schedule for overhang and frill.The frill used a free boundary and zero pressure throughout; the overhang used 6.0 × 10^5 initial steps and 4.0 × 10^5 subsequent steps.Source: Numerical parameters for the… · p. 13

Isometry error.

  • Isometry errors were computed from reference-to-wrinkle triangle-map Jacobians, weighted by reference area, and summarized by median and 95th percentile to reduce sensitivity to tiny distorted triangles.Reference faces with area below 10^-9 were excluded.Source: Isometry error · p. 13
  • Isometry error levels are comparable at A = 0.02 and A = 0.12, both before and after amplitude optimization.Table S2 reports area-weighted isometry error statistics for the two amplitudes and optimization conditions; the supplied text gives no individual values.Source: Isometry error · p. 13, Table S2 · p. 13

Agreement between the unfolded shape and the target.

  • Unfolded shapes were registered to targets without scaling, then compared using bidirectional unsigned point-to-surface distances normalized by target bounding-box diagonal.The reported distances include both remeshing and unfolding error because simulations used coarsened meshes.Source: Agreement between the unfolded shape… · p. 14

Physical implementation by 3D printing

  • All five target shapes were physically printed at A = 0.12, using thin flexible TPU crease bands and inversion regions, thicker PETG non-inverting regions, and PLA supports.Each design was scaled to a 100 mm longer bounding-rectangle side; crease bands were 0.4 mm thick and other regions 1.4 mm thick.Source: Physical implementation by 3D printing · p. 14
Loading 2610.34903v1…