A yield function in the stress resultant space of geometrically exact beams based on the elastoplastic cross-sectional warping problem has been proposed by Herrnböck et al. (Comput Mech, 67(3):723–742, 2021). This plasticity framework has been extended with a hardening tensor to model the kinematic hardening effects in Herrnböck et al. (Comput Mech, 71(1):1–24, 2022). While this framework provides scaling for the yield surface in ideal plasticity, scaling in hardening plasticity has not yet been explored. This paper focuses on the numeric modelling of hardening beams and beam assemblies at different geometric scales. Discretization effects from the introduction of plasticity into the geometrically exact beam model are demonstrated. Furthermore, the effects of scaling are explored, and a method to mitigate undesirable effects in order to achieve a size-agnostic formulation is proposed. Consistent geometric scaling is demonstrated for two alternative scaling approaches of the yield function.

Lightweight materials used for impact mitigation must be able to resist impact and absorb the maximum amount of energy from the impactor. Auxetic materials have the potential to achieve high resistance by drawing material into the impact zone and providing higher indentation and shear resistance. However, these materials must be artificially designed, and the large deformation dynamic effects of the created structures must be taken into consideration when deciding on a protection concept. Despite their promise, little attention has been given to understanding the working mechanisms of high-rate and finite deformation effects of architected auxetic lattice structures. This study compares the static and dynamic elastic properties of different auxetic structures with a honeycomb structure, a typical non-auxetic lattice, at equivalent mass and stiffness levels. In this study, we limit the investigation to elastic material behavior and do not consider contact between the beams of the lattices. It is demonstrated that the equivalent static and dynamic properties of individual lattices at an undeformed state are insufficient to explain the variations observed in impact situations. In particular, the initial Poisson's ratio does not determine the ability of a structure to resist impact. To gain a thorough comprehension of the overall behavior of these structures during localized, high rate compression, the evolution of the elastic tangent properties under compression and shear deformation was monitored, leading to a more profound understanding. Observations made in one configuration of stiffness and mass are replicated and analyzed in related configurations.

A sequential nonlinear multiscale method for the simulation of elastic metamaterials subject to large deformations and instabilities is proposed. For the finite strain homogenization of cubic beam lattice unit cells, a stochastic perturbation approach is applied to induce buckling. Then, three variants of anisotropic effective constitutive models built upon artificial neural networks are trained on the homogenization data and investigated: one is hyperelastic and fulfills the material symmetry conditions by construction, while the other two are hyperelastic and elastic, respectively, and approximate the material symmetry through data augmentation based on strain energy densities and stresses. Finally, macroscopic nonlinear finite element simulations are conducted and compared to fully resolved simulations of a lattice structure. The good agreement between both approaches in tension and compression scenarios shows that the sequential multiscale approach based on anisotropic constitutive models can accurately reproduce the highly nonlinear behavior of buckling-driven 3D metamaterials at lesser computational effort.