TY - JOUR
T1 - PolyDyna
T2 - a Matlab implementation for topology optimization of structures subjected to dynamic loads
AU - Giraldo-Londoño, Oliver
AU - Paulino, Glaucio H.
N1 - Funding Information:
We acknowledge Sandia National Laboratories, a multimission laboratory managed and operated by National Technology and Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International, Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525. We also thank the support provided by the National Science Foundation (NSF) under grant number 1663244. We are grateful to the insightful comments by Emily D. Sanders and Americo Cunha, which contributed to substantial improvements to the paper. The interpretation of the results of this work is solely that by the authors, and it does not necessarily reflect the views of the sponsors or sponsoring agencies.
Funding Information:
We acknowledge Sandia National Laboratories, a multimission laboratory managed and operated by National Technology and Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International, Inc., for the U.S. Department of Energy?s National Nuclear Security Administration under contract DE-NA0003525. We also thank the support provided by the National Science Foundation (NSF) under grant number 1663244. We are grateful to the insightful comments by Emily D. Sanders and Americo Cunha, which contributed to substantial improvements to the paper. The interpretation of the results of this work is solely that by the authors, and it does not necessarily reflect the views of the sponsors or sponsoring agencies.
Publisher Copyright:
© 2021, The Author(s), under exclusive licence to Springer-Verlag GmbH, DE part of Springer Nature.
PY - 2021/8
Y1 - 2021/8
N2 - We present a Matlab implementation for topology optimization of structures subjected to dynamic loads. The code, which we name PolyDyna, is built on top of PolyTop—a Matlab code for static compliance minimization based on polygonal finite elements. To solve the structural dynamics problem, we use the HHT-α method, which is a generalization of the classical Newmark-β method. In order to handle multiple regional volume constraints efficiently, PolyDyna uses a variation of the ZPR design variable update scheme enhanced by a sensitivity separation technique, which enables it to solve non-self-adjoint topology optimization problems. We conduct the sensitivity analysis using the adjoint method with the “discretize-then-differentiate” approach, such that the sensitivity analysis is consistently evaluated on the discretized system (both in space and time). We present several numerical examples, which are explained in detail and summarized in a library of benchmark problems. PolyDyna is intended for educational purposes and the complete Matlab code is provided as electronic supplementary material.
AB - We present a Matlab implementation for topology optimization of structures subjected to dynamic loads. The code, which we name PolyDyna, is built on top of PolyTop—a Matlab code for static compliance minimization based on polygonal finite elements. To solve the structural dynamics problem, we use the HHT-α method, which is a generalization of the classical Newmark-β method. In order to handle multiple regional volume constraints efficiently, PolyDyna uses a variation of the ZPR design variable update scheme enhanced by a sensitivity separation technique, which enables it to solve non-self-adjoint topology optimization problems. We conduct the sensitivity analysis using the adjoint method with the “discretize-then-differentiate” approach, such that the sensitivity analysis is consistently evaluated on the discretized system (both in space and time). We present several numerical examples, which are explained in detail and summarized in a library of benchmark problems. PolyDyna is intended for educational purposes and the complete Matlab code is provided as electronic supplementary material.
KW - Compliance minimization
KW - Elastodynamics
KW - HHT-α method
KW - Newmark-β method
KW - Sensitivity separation
KW - Topology optimization
KW - ZPR update scheme
UR - http://www.scopus.com/inward/record.url?scp=85109188425&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85109188425&partnerID=8YFLogxK
U2 - 10.1007/s00158-021-02859-6
DO - 10.1007/s00158-021-02859-6
M3 - Article
AN - SCOPUS:85109188425
SN - 1615-147X
VL - 64
SP - 957
EP - 990
JO - Structural Optimization
JF - Structural Optimization
IS - 2
ER -