[论文解读] An efficient 146-line 3D sensitivity analysis code of stress-based topology optimization written in MATLAB

Hao Deng, Praveen S. Vulimiri+1 more|arXiv (Cornell University)|Apr 2, 2021

Topology Optimization in Engineering参考文献 68被引用 6

一句话总结

本文提出了一种146行的MATLAB代码，用于基于p-范数全局应力度量和伴随敏感性分析的高效3D应力基拓扑优化。通过与有限差分法对比验证的解析导数，该方法实现了精确的敏感性计算，结合移动渐近线法（Method of Moving Asymptotes），成功实现了复杂3D结构中的体积约束应力最小化。

ABSTRACT

This paper presents an efficient and compact MATLAB code for three-dimensional stress-based sensitivity analysis. The 146 lines code includes the finite element analysis and p-norm stress sensitivity analysis based on the adjoint method. The 3D sensitivity analysis for p-norm global stress measure is derived and explained in detail accompanied by corresponding MATLAB code. The correctness of the analytical sensitivity is verified by comparison with finite difference approximation. The nonlinear optimization solver is chosen as the Method of moving asymptotes (MMA). Three typical volume-constrained stress minimization problems are presented to verify the effectiveness of sensitivity analysis code. The MATLAB code presented in this paper can be extended to resolve different stress related 3D topology optimization problems. The complete program for sensitivity analysis is given in the Appendix and is intended for educational purposes only.

研究动机与目标

• 开发一种简洁、高效且可重用的MATLAB实现，用于3D应力基拓扑优化。
• 推导并实现基于伴随法的p-范数全局应力度量的解析敏感性分析。
• 通过与有限差分近似比较，验证解析敏感性的准确性。
• 展示敏感性代码在求解3D体积约束应力最小化问题中的有效性。
• 为计算力学和拓扑优化领域的研究人员与学生提供一种教育工具。

提出的方法

• 该代码使用三维线性六面体单元进行有限元分析。
• 通过计算单元应力的p-范数，形成用于优化的全局应力度量。
• 采用伴随法进行敏感性分析，以计算应力度量对设计变量的导数。
• 解析导数被推导并直接实现于146行代码中。
• 采用移动渐近线法（Method of Moving Asymptotes, MMA）作为非线性优化器，求解应力最小化问题。
• 该实现结构清晰，便于扩展至其他与应力相关的3D拓扑优化问题。

实验结果

研究问题

• RQ1 如何在MATLAB中实现一种紧凑且高效的3D应力基拓扑优化敏感性分析代码？
• RQ2 与有限差分近似相比，p-范数全局应力度量的解析敏感性精度如何？
• RQ3 所提出的敏感性代码能否有效求解体积约束的3D应力最小化问题？
• RQ4 伴随法如何实现3D拓扑优化中敏感性的高效计算？
• RQ5 该146行代码在多大程度上可作为可重用的教育与研究工具，用于应力基拓扑优化？

主要发现

• 通过伴随法计算的解析敏感性与有限差分近似结果高度一致，验证了其正确性。
• 146行MATLAB代码成功在一个紧凑实现中完成3D有限元分析与p-范数应力敏感性计算。
• 三个基准的体积约束应力最小化问题被有效求解，展示了代码的鲁棒性与高效性。
• 该代码适用于教育用途，并可扩展至其他与应力相关的3D拓扑优化问题。
• 该实现保持了高性能与清晰性，使计算力学领域的研究人员与学生易于使用。
• 采用p-范数全局应力度量可使3D优化问题更加平滑且条件良好。