[Paper Review] Continuum Model for Pressure Actuated Cellular Structures

Markus Pagitz, Remco I. Leine|arXiv (Cornell University)|Mar 17, 2014

Advanced Materials and Mechanics12 references3 citations

TL;DR

This paper introduces a coupled continuum and numerical model for optimizing pressure-actuated cellular structures, enabling simultaneous shape and stress optimization by linking continuum geometry to spring stiffnesses and hinge eccentricities via stress constraints. The method achieves full coupling between models, yielding manufacturable, fully-stressed designs that meet target shapes with minimal prototyping iterations.

ABSTRACT

Previous work introduced a lower-dimensional numerical model for the geometric nonlinear simulation and optimization of compliant pressure actuated cellular structures. This model takes into account hinge eccentricities as well as rotational and axial cell side springs. The aim of this article is twofold. First, previous work is extended by introducing an associated continuum model. This model is an exact geometric representation of a cellular structure and the basis for the spring stiffnesses and eccentricities of the numerical model. Second, the state variables of the continuum and numerical model are linked via discontinuous stress constraints on the one hand and spring stiffness, hinge eccentricities on the other hand. An efficient optimization algorithm that fully couples both sets of variables is presented. The performance of the proposed approach is demonstrated with the help of an examples.

Motivation & Objective

• To develop a continuum model that exactly represents the geometric and mechanical behavior of compliant pressure-actuated cellular structures.
• To establish a fully coupled optimization framework between the continuum model and a lower-dimensional numerical model.
• To enable direct manufacturability of optimized designs by ensuring stress and geometric constraints are met simultaneously.
• To reduce prototyping time and cost by generating optimized, fabrication-ready designs directly from the model.

Proposed method

• The continuum model is decomposed into rigid cell corners and elastic cell sides, with hinges modeled as circular cutouts with finite bending and infinite axial stiffness.
• Spring stiffnesses and hinge eccentricities in the numerical model are derived from the continuum model's geometry and stress distribution.
• Maximum hinge and cell side stresses are computed using von Mises yield criteria with a stress reduction factor ρ = 0.9.
• Rotational stiffness is computed via detailed finite element analysis, while axial stiffness is derived from a lower-dimensional model.
• A bilevel optimization approach is used to determine optimal cell corner geometries under equilibrium and stress constraints.
• An efficient optimization algorithm couples state variables of both models, with discontinuous stress constraints handled via iterative refinement.

Experimental results

Research questions

• RQ1 How can a continuum model be constructed to exactly represent the geometric and mechanical behavior of pressure-actuated cellular structures?
• RQ2 What is the optimal way to link the continuum model’s stress and geometry to the spring stiffnesses and hinge eccentricities in the numerical model?
• RQ3 How can both models be fully coupled to satisfy target shapes and stress constraints simultaneously?
• RQ4 Can the resulting optimized design be directly manufactured via rapid prototyping without iterative physical testing?

Key findings

• The optimized structure achieved equilibrium shapes closely matching the target half-circle and full-circle configurations under two distinct pressure sets.
• The optimization converged to the desired target shapes in 25 iterations, but required 41 iterations to simultaneously satisfy shape and maximum stress constraints.
• All hinges and central cell sides were fully stressed (σ_max = 85 MPa) for at least one of the two pressure sets, indicating optimal load utilization.
• The maximum step length in the optimization was limited to Δv_max = 5 mm and Δt_max = 0.6 mm, ensuring stable convergence.
• Convergence exhibited zigzag behavior after 15 iterations due to discontinuous stress constraints, which shift as the pressure set causing maximum stress changes during optimization.
• The final optimized continuum model was geometrically distinct from the initial design, with varying thicknesses and lengths along the structure, confirming the necessity of full coupling.