[논문 리뷰] An efficient 146-line 3D sensitivity analysis code of stress-based topology optimization written in MATLAB

Hao Deng, Praveen S. Vulimiri외 1명|arXiv (Cornell University)|2021. 04. 02.

Topology Optimization in Engineering참고 문헌 68인용 수 6

한 줄 요약

이 논문은 p-노름 전역 응력 측정과 양면 감도 분석을 사용한 효율적인 3D 응력 기반 토폴로지 최적화를 위한 146줄의 MATLAB 코드를 제시한다. 이 방법은 유한 차분과 대조되는 분석 도함수를 통해 정확한 감도 계산을 실현하며, 이는 방법의 이동 점근법(MMA)을 사용하여 복잡한 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.

연구 동기 및 목표

• 3D 응력 기반 토폴로지 최적화를 위한 간결하고 효율적이며 재사용 가능한 MATLAB 구현을 개발하는 것.
• 양면 방법을 사용하여 p-노름 전역 응력 측정의 분석 감도 분석을 유도하고 구현하는 것.
• 유한 차분 근사와의 비교를 통해 분석 감도의 정확성을 검증하는 것.
• 감도 코드가 3D에서 체적 제약 조건을 갖는 응력 최소화 문제를 효과적으로 해결할 수 있음을 보여주는 것.
• 계산 역학 및 토폴로지 최적화 분야의 연구자들과 학생들에게 교육적 도구로 활용될 수 있도록 제공하는 것.

제안 방법

• 코드는 3D에서 선형 육면체 요소를 사용하여 유한 요소 분석을 수행한다.
• 요소 응력의 p-노름을 계산하여 최적화를 위한 전역 응력 측정을 형성한다.
• 설계 변수에 대한 응력 측정의 도함수를 계산하기 위해 양면 기반 감도 분석을 사용한다.
• 감도 도함수는 분석적으로 유도되어 146줄의 코드에 직접 구현된다.
• 비선형 최적화기로 방법의 이동 점근법(MMA)을 사용하여 응력 최소화 문제를 해결한다.
• 해당 구현은 명확성과 다른 응력 관련 3D 토폴로지 최적화 문제로의 확장성 확보를 위해 구성되어 있다.

실험 결과

연구 질문

• RQ1 어떻게 간결하고 효율적인 3D 감도 분석 코드를 MATLAB을 사용하여 응력 기반 토폴로지 최적화에 구현할 수 있는가?
• RQ2 p-노름 전역 응력 측정에 대한 분석 감도의 정확도는 유한 차분 근사와 비교하여 어떻게 되는가?
• RQ3 제안된 감도 코드는 체적 제약 조건을 갖는 3D 응력 최소화 문제를 효과적으로 해결할 수 있는가?
• RQ4 양면 방법은 3D 토폴로지 최적화에서 감도 계산을 어떻게 효율적으로 가능하게 하는가?
• RQ5 이 146줄의 코드는 얼마나 널리 응력 기반 토폴로지 최적화의 재사용 가능한 교육적 및 연구 도구로 활용될 수 있는가?

주요 결과

• 양면 방법을 통해 계산된 분석 감도는 유한 차분 근사와 뛰어난 일치를 보이며 정확성을 확인한다.
• 146줄의 MATLAB 코드는 하나의 컴act한 구현 내에서 3D 유한 요소 분석과 p-노름 응력 감도 계산을 성공적으로 수행한다.
• 세 가지 벤치마크 체적 제약 조건을 갖는 응력 최소화 문제를 효과적으로 해결하여 코드의 강력성과 효율성을 입증하였다.
• 코드는 교육적 용도로 적합하며 다른 응력 관련 3D 토폴로지 최적화 문제로의 확장이 가능하다.
• 해당 구현은 높은 성능과 명확성을 유지하여 계산 역학 분야의 연구자들과 학생들에게 접근하기 용이하다.
• p-노름 전역 응력 측정의 사용은 3D에서 매끄럽고 잘 조절된 최적화 문제를 가능하게 한다.