A Continuum Macro-Model for Bistable Periodic Auxetic Surfaces

A macro-constitutive model for the deformation response of periodic rotating bistable auxetic surfaces is developed. Focus is placed on isotropic surfaces made of bistable hexagonal cells composed of six triangular units with two stable equilibrium states. Adopting a variational formulation, the effective stress-strain response is derived from a free energy function expressed in terms of the invariants of the logarithmic strain. A regularization of the governing equations via a gradient-enhanced first invariant of the logarithmic strain is introduced since the double-well nature of the free energy may result in mathematical ill-posedness and related numerical artifacts, such as mesh sensitivity. Despite this regularization, the numerical scheme may still suffer from divergence issues due to the highly non-linear material behavior. To enhance numerical stability, an artificial material rate-dependency is additionally introduced. Although it does not guarantee solution uniqueness or eliminate mesh sensitivity, it is conjectured to assist the numerical scheme in overcoming snap-backs caused by local non-proportional loading induced by transition fronts. The model is implemented using membrane/shell structural elements and plane stress continuum ones within the ABAQUS finite element suite. Numerical simulations demonstrate the efficacy of the proposed formulation and its implementation.