High-performance large-scale simulation of multi-stable metastructures
High-performance large-scale simulation of multi-stable metastructures
复制标题
DOI：
10.1016/j.cpc.2022.108365
复制
发表时间：
2022-04
期刊：
影响因子：
--
通讯作者：
Myungwon Hwang;C. Scalo;A. F. Arrieta
中科院分区：
文献类型：
--
作者：
Myungwon Hwang;C. Scalo;A. F. Arrieta
在本文中，我们开发了一个基于消息传递接口（MPI）的求解器，以实现双稳定或多稳定元素组成的通用元结构的快速大规模模拟。内部求解器已经对商业数值求解器（Abaqus）和先前研究中建立的串行代码进行了彻底验证。我们可以通过完全显式的龙格-库塔（RK）方法实现高达四阶解的精度，超过许多商业结构分析工具提供的精度。我们的并行代码专门用于解决特定类型的问题，绝对计算速度可以提高三个数量级，从而可以研究大的参数空间。更重要的是，内部实现可以根据固有的结构周期性有效地分配计算负载，从而实现高效的并行可扩展性。为了证明我们的代码能够处理以前用现有求解器无法解决的大规模问题，我们研究了双稳定元束的振幅相关能量传输率和过渡波传播速度的稳定性。所取得的数值和计算性能的提高极大地扩展了一般非线性超材料和元结构体系结构的可访问分析领域，从而开辟了发现新动力学和实现实际应用的潜力。程序摘要程序标题：NM³（非线性超材料MPI） solverCPC库链接到程序文件：https://doi.org/10.17632/8f4n99jccf.1Developer's存储库链接：https://github.com/wonnie87/NMCubeLicensing条款：mit编程语言：fortran问题的性质：NM³能够大规模并行模拟强非线性超材料和元结构，包括一维多稳定晶格与耦合摆（离散sine-Gordon模型），具有四次位势的一维晶格（离散的ϕ-4模型）和具有双稳定微观结构的亚梁。求解方法：在NM³中实现了高达四阶的显式龙格-库塔（RK）方法。如果需要无条件的数值稳定性，具有恒定平均加速度的Newmark-β（隐式）方法也可用。附加注释，包括限制和不寻常的功能：运行NM / 3需要安装Python（带NumPy库），MPI和HDF5。使用Python脚本生成输入文件。代码使用MPI系统调用来允许计算进程之间的大规模并行化。该代码使用HDF5文件格式进行数据存储。
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.Program summaryProgram Title:NM̂3 (Nonlinear MetaMaterials MPI) solverCPC Library link to program files:https://doi.org/10.17632/8f4n99jccf.1Developer's repository link:https://github.com/wonnie87/NMCubeLicensing provisions:MITProgramming language:FortranNature of problem:NM̂3 enables massively parallel simulations of strongly nonlinear metamaterials and metastructures, including 1D multi-stable lattice with coupled pendula (discrete sine-Gordon model), 1D lattice with quartic on-site potentials (discreteϕ-4 model), and metabeam with a bi-stable microstructure.Solution method:Up to the 4th-order explicit Runge-Kutta (RK) methods are implemented in NM̂3. The Newmark-β(implicit) method with constant average acceleration is also available if unconditional numerical stability is desired.Additional comments including restrictions and unusual features:Running NM̂3 requires installation of Python (with NumPy library), MPI, and HDF5. A Python script is used to generate input files. The code use MPI system calls to allow a massive parallelization among the compute processes. The code uses HDF5 file format for data storage.