본 논문에서는 역전파 방법 기반 자동미분법을 이용하여 설계민감도를 구하고 이를 응력제한조건을 고려한 위상최적설계에 적용 하였다. 응력제한조건이 있는 위상최적화문제는 특이점(singularity)과 응력의 국부성(local nature of stress constraint)문제, 그리고 설 계 변수에 대한 비선형성의 문제를 포함하고 최적해를 얻기가 매우 힘들다. 특이점 문제를 해결하기 위해서 응력 완화(stress relaxation) 기법을 사용하였고, 응력의 국부성을 해결하기 위해 p-norm을 이용한 전역 응력치를 제한조건에 사용하였다. 설계 변수에 대한 비선 형성을 극복하기 위해 해석적인 방법으로 정확한 설계민감도를 구하는 것이 중요하다. 위상최적설계에서 기존에는 보조변수방법 (adjoint variable method)을 사용하여 빠르고 정확한 설계민감도를 구했지만, 설계민감도를 해석적으로 구해야 하고, 보조평형방정식 을 추가로 풀어야 하는 어려움이 있다. 이를 해결하기 위해서 인공신경망에서 최적 가중치(weights)와 편차(biases)를 구할 때 쓰이는 역전파 기법을 이용하여 설계민감도를 구하고 이를 응력제한조건을 고려한 위상최적설계에 적용하였다. 역전파 기법은 자동미분에 쓰이는 기법으로 목적함수나 제한조건에 대한 설계민감도를 별도의 수식유도 없이 간단하게 구할 수 있는 장점이 있다. 또한, 미분값 을 구하는 역전파의 과정이 보조평형방정식을 푸는 것보다 계산시간이 빠르고 해석적 방법으로 구한 설계민감도와 같은 정확도를 보 여준다

레벨셋 기법과 위상민감도를 이용하여 선형 탄성 구조물에 대하여, 초기 설계형상에 의존성이 없는 위상 및 형상 최적설계 기법을 개발하였다. 레벨셋 기법에서는 복잡한 위상 형상변화를 쉽게 다루기 위해 초기 영역은 고정한 채 레벨셋 함수로 표현되는 암시적 이동경계로 경계를 표현한다. 해밀턴-자코비(H-J) 방정식과 수치적으로 강건한 기법인 ‘up-wind scheme’은 컴플라이언스 목적함수를 최소화시키고 허용체적 제약조건을 만족시키면서, 초기 암시적 경계를 법선 속도장에 따라 최적의 형상으로 이끌어 낸다. 점근적인 정규화 개념에 근거하여, 구멍의 반지름을 0으로 접근시켜 형상 미분의 극한을 취한 위상민감도를 고려하였다. 최적조건으로부터 유도된 라그란지안의 감소 방향을 이용하여 H-J 방정식을 갱신하기 위한 속도장을 결정하였다. 개발한 방법에서는 위상민감도로부터 얻어지는 지표를 이용하여 구멍을 언제든지 어디에서나 생성가능하기 때문에 초기 구멍이 최적 형상을 얻기 위해 요구되지 않는다는 사실을 확인하였다. 또한 효율적인 최적화 과정을 위해서는 구멍 생성을 위한 조정변수의 적절한 선택이 중요함을 확인하였다.

A level set based topological shape optimization method for nonlinear structure considering hyper-elastic problems is developed. To relieve significant convergence difficulty in topology optimization of nonlinear structure due to inaccurate tangent stiffness which comes from material penalization of whole domain, explicit boundary for exact tangent stiffness is used by taking advantage of level set function for arbitrary boundary shape. For given arbitrary boundary which is represented by level set function, a Delaunay triangulation scheme is used for current structure discretization instead of using implicit fixed grid. The required velocity field in the actual domain to update the level set equation is determined from the descent direction of Lagrangian derived from optimality conditions. The velocity field outside the actual domain is determined through a velocity extension scheme based on the method suggested by Adalsteinsson and Sethian(1999). The topological derivatives are incorporated into the level set based framework to enable to create holes whenever and wherever necessary during the optimization.