[논문 리뷰] Volumetric expressions of the shape gradient of the compliance in structural shape optimization

Matteo Giacomini, Olivier Pantz외 1명|arXiv (Cornell University)|2017. 01. 20.

Topology Optimization in Engineering참고 문헌 31인용 수 8

한 줄 요약

이 논문은 선형 탄성에서 강성의 형상 기울기의 부피 표현을 순수 이동 및 이중 혼합 유한요소 형식을 사용하여 유도한다. 수치 시뮬레이션을 통해 표면 기반 방법에 비해 부피 기반 방법이 메쉬 변형 상황에서도 더 정확하고 안정적인 내림방향을 제공함을 보이며, 메쉬 왜곡과 기울기 근사 오차로 인한 형상 최적화의 과제를 부각시킨다.

ABSTRACT

In this article, we consider the problem of optimal design of a compliant structure under a volume constraint, within the framework of linear elasticity. We introduce the pure displacement and the dual mixed formulations of the linear elasticity problem and we compute the volumetric expressions of the shape gradient of the compliance by means of the velocity method. A preliminary qualitative comparison of the two expressions of the shape gradient is performed through some numerical simulations using the Boundary Variation Algorithm.

연구 동기 및 목표

• 형상 최적화에서 경계에만 의존하는 표현을 피하기 위해 선형 탄성에서 강성 최소화를 위한 부피 기반 형상 기울기 표현을 개발하는 것.
• 선형 탄성 문제의 순수 이동 및 이중 혼합 형식에서 유도된 부피 기반 형상 기울기의 수치 성능를 비교하는 것.
• 이산 부피 기반 형상 기울기를 사용할 때 경계 변화 알고리즘(Boundary Variation Algorithm)을 통해 계산된 내림방향의 안정성과 정확도를 조사하는 것.
• 메쉬 왜곡과 기울기 근사 오차로 인한 수치적 불안정성, 예를 들어 음수 방향 도함수에도 불구하고 내림방향이 아닌 행동을 보이는 현상을 규명하는 것.
• 구형 내림방향 알고리즘과 레벨셋 방법을 사용하여 구조 형상 최적화의 강건성을 향상시키기 위한 동기를 제공하는 것.

제안 방법

• 변분 프레임워크 내에서 속도 방법을 사용하여 강성 기능의 형상 기울기를 유도하고, 도메인 매핑을 통해 경계 적분을 부피 적분으로 변환하는 것.
• 이동이 주요 미지수인 순수 이동 형식의 선형 탄성 문제를 적용하고, 이에 해당하는 부피 기반 형상 기울기 표현을 유도하는 것.
• 헬링거-라이스너 형식과 약하게 강제된 응력 대칭성을 가진 변형을 포함한 두 가지 이중 혼합 형식을 도입하며, 이들 모두가 부피 기반 형상 기울기 표현을 제공하는 것.
• 계산된 형상 기울기를 기반으로 도메인을 반복적으로 변형하는 데 사용되는 경계 변화 알고리즘(Boundary Variation Algorithm, BVA)을 적용하며, 메쉬 업데이트는 내림방향에 의해 유도된다.
• 사후 오차 추정기법과 인증된 내림방향 알고리즘(Certified Descent Algorithm, CDA)을 사용하여 이산 형상 기울기에서 비내림방향 행동을 탐지하고 수정하는 것.
• 다양한 초기 기하 구조를 가진 2차원 테스트 케이스에서 수치 시뮬레이션을 수행하여 순수 이동 형식과 이중 혼합 형식의 수렴성과 안정성을 비교하는 것.

실험 결과

연구 질문

• RQ1 선형 탄성에서 강성의 형상 기울기를 표면 적분을 피하는 순수 부피 형태로 어떻게 표현할 수 있는가?
• RQ2 순수 이동 형식과 이중 혼합 형식에서 유도된 부피 기반 형상 기울기를 사용할 때의 수치적 장점과 단점은 무엇인가?
• RQ3 경계 변화 알고리즘이 음수 방향 도함수를 가질 때에도 진정한 내림방향을 생성하지 못하는 이유는 무엇인가?
• RQ4 메쉬 왜곡과 속도 장의 낮은 정(regularity)이 형상 최적화 알고리즘의 수렴성에 어느 정도 영향을 미치는가?
• RQ5 이산 부피 기반 형상 기울기를 사용할 때 인증된 내림방향 알고리즘이 신뢰성 있는 형상 최적화를 향상시킬 수 있는가?

주요 결과

• 순수 이동 및 이중 혼합 형식에서 유도된 부피 기반 형상 기울기 표현은 연속적인 설정에서 고전적인 표면 기반 형상 기울기와 수학적으로 동일하다.
• 수치 시뮬레이션 결과, 이중 혼합 형식이 순수 이동 형식보다 더 안정적이고 정확한 내림방향을 제공하는 것으로 나타났으며, 특히 초기 반복 단계에서 두드러졌다.
• ⟨dₕL(Ω), θₕ⟩ < 0와 같은 음수 방향 도함수에도 불구하고, 페널티 강성 기능 L(Ω)이 몇 차례 반복에서 증가하는 경우가 있었으며(예: 표 1c에서 10.64 → 10.66), 이는 θₕ가 진정한 내림방향이 아니었음을 시사한다.
• 순수 이동 형식을 사용한 테스트 케이스 (c)에서, 강성은 30회 반복 동안 10.64에서 4.66으로 감소했지만, 알고리즘이 진동 행동을 보이며 내림방향의 불안정성을 시사했다.
• 테스트 케이스 (d)에서 이중 혼합 형식은 더 나은 수렴 행동을 보였으며, 강성은 30회 반복 동안 11.01에서 3.37로 감소했지만 여전히 진동이 관찰되었다.
• 이 연구는 표준 BVA가 기울기 근사 부정확성이나 메쉬 왜곡으로 인해 실패할 수 있으므로, 이산 형상 기울기의 내림성질을 검증하기 위해 인증된 내림방향 알고리즘의 필요성을 확인한다.