[논문 리뷰] Shape optimization to decrease failure probability

Matthias Bolten, Hanno Gottschalk외 2명|arXiv (Cornell University)|2017. 05. 16.

Topology Optimization in Engineering참고 문헌 15인용 수 3

한 줄 요약

이 논문은 세라믹 부품의 실패 확률을 최소화하기 위해 Weibull 분포를 확률적 강도 모델로 사용하는 새로운 기울기 기반 형상 최적화 프레임워크를 제안한다. 실패 확률을 PDE 제약 조건이 있는 목적 함수로 설정하고, 유한 요소 메esh에서 애드조인트 방법을 통해 형상 기울기를 계산함으로써 안정적인 수렴을 이룩하였다. 이는 직선 보와 S자형 조인트에 대한 수치 실험을 통해 검증되었으며, 인위적 스무딩 없이도 실패 확률을 크게 감소시켰다.

ABSTRACT

Ceramic is a material frequently used in industry because of its favorable properties. Common approaches in shape optimization for ceramic structures aim to minimize the tensile stress acting on the component, as it is the main driver for failure. In contrast to this, we follow a more natural approach by minimizing the component's probability of failure under a given tensile load. Since the fundamental work of Weibull, the probabilistic description of the strength of ceramics is standard and has been widely applied. Here, for the first time, the resulting failure probabilities are used as objective functions in PDE constrained shape optimization. To minimize the probability of failure, we choose a gradient based method combined with a first discretize then optimize approach. For discretization finite elements are used. Using the Lagrangian formalism, the shape gradient via the adjoint equation is calculated at low computational cost. The implementation is verified by comparison of it with a finite difference method applied to a minimal 2d example. Furthermore, we construct shape flows towards an optimal / improved shape in the case of a simple beam and a bended joint.

연구 동기 및 목표

• 피크 응력 기준에 의존하지 않고, 세라믹 부품의 실패 확률을 직접 최소화하는 형상 최적화 방법을 개발하는 것.
• 특히 Weibull 모델을 활용한 확률적 균열역학을 형상 최적화의 목적 기능으로 통합하는 것.
• 유한 요소를 사용한 '먼저 이산화한 후 최적화' 접근 방식을 통해 실패 확률에 대한 기울기 기반 최적화를 가능하게 하는 것.
• 유한 차분 비교를 통해 계산된 형상 기울기의 정확성과 안정성을 검증하는 것.
• 체적 제약 조건 하에서 직선화된 보와 최적화된 S자형 조인트와 같은 향상된 형상으로 수렴하는지 검증하는 것.

제안 방법

• 실패 확률은 Weibull 분포를 사용하여 모델링되며, 척도 매개변수는 부품 형상과 하중에 따라 달라지며, 형태 매개변수(Weibull 모듈러스)는 재질에 따라 고정된다.
• 목적 기능은 주어진 인장 하중 하에서의 실패 확률이며, 탄성 균열역학과 초기 균열의 점 과정 모델에서 유도된다.
• 선형 탄성 문제의 PDE를 이산화하고 다각형 메쉬에서 변위와 응력을 계산하기 위해 유한 요소 방법을 사용한다.
• 형상 기울기는 애드조인트 방법을 통해 계산되며, 경계 변화에 따른 실패 확률의 도함수를 효율적으로 계산할 수 있다.
• 체적 제약 조건은 체적 기울기의 영공간에 형상 기울기를 투영시켜 일정한 체적을 유지하도록 구현된다.
• 경계 노드를 형상 기울기를 사용하여 갱신하고, 메쉬 품질을 유지하기 위해 반복적으로 메쉬 변형을 적용함으로써 기하학적 흐름을 생성한다.

실험 결과

연구 질문

• RQ1 기울기 기반 형상 최적화 프레임워크를 통해 세라믹 부품의 실패 확률을 효과적으로 최소화할 수 있는가?
• RQ2 이산적인 유한 요소 메쉬에서 애드조인트 방법을 사용해 실패 확률에 대한 형상 기울기를 얼마나 정확하게 계산할 수 있는가?
• RQ3 특히 직관적으로 명백하지 않은 최적 형상이 존재하는 경우, 인위적 스무딩 없이 최적화 과정이 안정적으로 수렴하는가?
• RQ4 최적화 과정 중 실패 확률은 어떻게 변화하며, 상당한 감소를 달성할 수 있는가?
• RQ5 비연속적이거나 내측 굴곡이 있는 형상의 경우 최적화의 수렴성과 안정성을 확보하는 데 있어 수학적·수치적 과제는 무엇인가?

주요 결과

• 체적 제약 조건 하에서 굽은 보의 실패 확률 최적화 알고리즘이 분석적으로 예상되는 최적 형상인 직선 막대로 수렴하는 데 성공하였다.
• S자형 조인트의 경우, 명백한 최적 구조가 없음에도 불구하고 안정적인 기하학적 흐름을 생성하고 실패 확률을 상당히 감소시켰다.
• 애드조인트 방법을 통해 계산된 형상 기울기는 최소한의 2차원 예제에서 유한 차분과의 비교를 통해 높은 수치 정밀도를 확보하였다.
• 인위적 스무딩 없이도 큰 기하학적 변형 조건에서도 최적화 과정이 안정적이고 수렴성을 유지하였다.
• 반복 과정 동안 시각화된 실패 확률 분포는 실패 위험을 지속적이고 단조롭게 감소시키는 일관된 개선을 보였다.
• 피크 응력 기준의 미분 불가능성 문제를 극복하고, 확률적 기계적 무결성을 형상 최적화의 직접적인 目적 기능으로 사용할 수 있음을 입증하였다.