[논문 리뷰] Overhang penalization in additive manufacturing via phase field structural topology optimization with anisotropic energies

Harald Garcke, Kei Fong Lam외 2명|arXiv (Cornell University)|2021. 11. 28.

Topology Optimization in Engineering인용 수 3

한 줄 요약

이 논문은 적층 제조에서 오버행을 억제하기 위해 비등방성 에너지 함수를 도입한 단계장 최적화 방법을 제안한다. 이는 외부 지지 구조 없이도 자가 지지 구조를 보장한다. 볼록 및 비볼록 비등방성 특성을 통합함으로써, 부분미분 기하학과 점점 가까워지는 분석을 통해 오버행 각도 제약 조건을 강제하고, 특히 강한 비등방성 조건에서 더 기울어진 자가 지지 구조를 생성한다. 이는 캔틸레버 빔과 다리 설계에 대한 수치 시뮬레이션을 통해 입증되었다.

ABSTRACT

A phase field approach for structural topology optimization with application to additive manufacturing is analyzed. The main novelty is the penalization of overhangs (regions of the design that require underlying support structures during construction) with anisotropic energy functionals. Convex and non-convex examples are provided, with the latter showcasing oscillatory behavior along the object boundary termed the dripping effect in the literature. We provide a rigorous mathematical analysis for the structural topology optimization problem with convex and non-continuously-differentiable anisotropies, deriving the first order necessary optimality condition using subdifferential calculus. Via formally matched asymptotic expansions we connect our approach with previous works in the literature based on a sharp interface shape optimization description. Finally, we present several numerical results to demonstrate the advantages of our proposed approach in penalizing overhang developments.

연구 동기 및 목표

• 외부 지지 구조가 필요로 하는 적층 제조에서의 오버행 문제를 해결함으로써 비용과 후처리 복잡성을 감소시키는 것.
• 최소 오버행 각도 제약 조건을 강제함으로써 최적화 과정에서 오버행을 억제할 수 있는 수학적 프레임워크를 개발하는 것.
• 제조 방향에 따라 오버행 영역을 억제함으로써 자가 지지 형상을 선호하는 비등방성 에너지 함수를 도입하는 것.
• 비연속적으로 미분 가능한 비등방성 특성에 대해 부분미분 기하학을 사용한 엄밀한 수학적 분석을 수행하는 것.
• 형상 최적화의 날카운면 형태 모델과 단계장 모델 간의 연결 고리를 형식적로 일치하는 점점 가까워지는 전개를 통해 확립하는 것.

제안 방법

• 설계를 단계장 변수 ϕ ∈ [0,1]로 표현하는 단계장 접근법을 제안하며, 여기서 ϕ = 1은 재료, ϕ = 0은 빈공간을 나타낸다.
• 오버행 각도가 임계 기준 ψ 이하인 영역을 억제하기 위해 비등방성 에너지 함수 γ(ν)를 도입하며, ν는 표면의 단위 법선 벡터이다.
• 강도와 방향성을 제어할 수 있는 볼록 정규화 비등방성 (식 6.2)를 사용하며, 매개변수 α는 비등방성의 강도를 조절한다. α = 1일 경우 등방성 행동을 나타낸다.
• 비볼록성 및 비연속적으로 미분 가능한 비등방성 특성을 다루기 위해 부분미분 기하학을 사용하여 1차 최적성 조건을 유도한다.
• 형상 최적화의 날카운면 형태 모델과의 연결을 위해 형식적로 일치하는 점점 가까워지는 전개를 수행하여, 단계장 모델이 고전적 수식과 일관성을 갖는다는 것을 검증한다.
• 유한 요소 방법을 사용하여 시간 간격을 적응적으로 조정하고 질량 보존 조건을 포함한 PDE 제약 최적화 문제를 수치적으로 해결한다.

실험 결과

연구 질문

• RQ1 적층 제조에서 오버행을 억제하면서도 구조적 안정성을 확보하기 위해 비등방성 에너지 함수를 어떻게 수학적으로 설정할 수 있는가?
• RQ2 볼록 비등방성과 비볼록 비등방성의 영향이 최적화 과정에서 진동하는 경계층(드리핑 효과) 형성에 미치는 영향은 무엇인가?
• RQ3 제안된 비등방성 에너지를 갖는 단계장 모델은 고전적 날카운면 형태 최적화 모델과 어떤 관계가 있는가?
• RQ4 비등방성 매개변수 α의 선택이 캔틸레버 및 다리 설계에서 자가 지지 구조의 기하학적 형태에 미치는 영향은 어느 정도인가?
• RQ5 이 방법은 외부 지지 구조 없이도 자가 지지 구조를 생성할 수 있으며, 이는 등방성 또는 기존 최적화 결과와 비교해 어떻게 다른가?

주요 결과

• 제안된 방법은 특히 강한 비등방성 조건(α = 0.2)에서 내부 브릿지가 거의 수직이 되도록 하여 캔틸레버 빔에서 오버행을 효과적으로 억제한다.
• α = 0.2일 경우, 내부 연결 부품이 없는 구성으로 수렴함으로써 자가 지지 형상에 대한 강한 선호도를 보인다.
• 수치 결과에 따르면, α 값이 작을수록 특정 법선 방향에서 γ(ν) 값이 높아서 상세한 표면 전이 영역(단계장 전이 영역)이 두꺼워지며, 이는 점점 가까워지는 분석의 예측과 일치한다.
• 다리 설계에서는 비등방성 사례(α = 0.7, 0.5)에서 등방성 사례(α = 1)보다 아래쪽면이 더 날카롭고 명확한 모서리를 갖는다. 이는 오버행 제약 조건을 더 잘 이행하고 있음을 시사한다.
• ε를 1/(32π)에서 1/(64π)로 줄임으로써 상세한 표면 전이 영역이 얇아지며, 이는 단계장 모델이 날카운면 근사로 수렴하는 것을 확인한다.
• 비볼록 비등방성 조건에서 드리핑 효과가 효과적으로 억제되기는 하나, 경계 근처에서는 진동하는 행동이 관찰되며, 이는 이전 문헌과 일치한다.