High-performance large-scale simulation of multi-stable metastructures

In this paper, we have developed a solver based on the message-passing interface (MPI) to enable rapid large-scale simulation of generic metastructures composed of bi- or multi-stable elements. The in-house solver has been thoroughly validated against a commercial numerical solver (Abaqus) and the well-established serial codes from the previous studies. We can achieve up to 4th-order solution accuracy with fully explicit Runge-Kutta (RK) methods, exceeding what many commercial structural analysis tools provide. With our parallel code dedicated to solving specific problem types, the absolute computational speed can be improved by three orders of magnitude, enabling the investigation of a large parameter space. More importantly, the in-house implementation enables an effective distribution of the computational load following the intrinsic structural periodicity, thus achieving efficient parallel scalability. To demonstrate our code's capability to handle massively large problems previously unattainable with existing solvers, we investigate the amplitude-dependent energy transmissibility of bi-stable metabeams and the stability of the transition wave's propagation speed. The achieved numerical and computational performance gains drastically expand the accessible analysis domains of general nonlinear metamaterial and metastructure architectures, thus opening up the potential to uncover new dynamics and enable practical implementations.