형상 최적화를 통한 포노닉 밴드갭 구조 설계

초록

본 논문은 연성 포함재와 강성 매트릭스로 구성된 강하게 이질적인 탄성 복합재의 동질화 모델을 이용해 저주파 영역의 음향 밴드갭을 최적화한다. 포함재의 형상을 매개변수화하고, 유효 질량 텐서의 고유값 민감도를 계산해 그래디언트 기반 최적화 문제를 정의한다. 2차원 사례를 통해 강성 제약 하에서 밴드갭 폭을 확대하는 최적 형상이 도출되었다.

상세 분석

이 연구는 두 단계의 핵심 아이디어에 기반한다. 첫 번째는 강한 이질성을 갖는 탄성 복합재를 비스케일(ε) 전개법으로 동질화하여, 유효 탄성 4차 텐서와 유효 질량 2차 텐서를 도출하는 것이다. 특히, 포함재가 매우 연성(ε² 스케일)으로 가정되므로, 유효 질량 텐서는 주파수 의존적인 비대칭 행렬이 되며, 특정 주파수 구간에서 부정정(negative‑definite)으로 전이될 때 밴드갭이 발생한다. 두 번째 아이디어는 이러한 밴드갭을 형상 변수에 대한 민감도 분석을 통해 직접 최적화한다는 점이다. 포함재 경계는 원형 B‑스플라인으로 매개화하여 충분한 C¹ 연속성을 보장하고, 형상 변화에 따른 셀 문제와 고유값 문제의 변분 형태를 이용해 파라미터에 대한 도함수를 구한다. 여기서 핵심은 유효 질량 텐서의 고유값 λᵢ(α) (α는 형상 파라미터)의 도함수 dλᵢ/dα 를 정확히 계산하는 것이며, 이를 위해 셀 문제의 해와 고유모드의 정규화 조건을 동시에 고려한 연쇄법칙을 적용한다.

민감도 식은 두 부분으로 나뉜다. (1) 셀 문제의 계수(탄성 텐서, 밀도)의 변형에 따른 직접 항, (2) 경계 변형에 따른 영역 적분 경계항이다. 후자는 변형률 텐서와 표면 법선의 변화를 포함하며, 레벨셋 방식이 아닌 B‑스플라인 파라미터화 덕분에 해석적 형태로 정리할 수 있다. 이렇게 얻은 전체 그래디언트는 L‑BFGS와 같은 제한 메모리 quasi‑Newton 방법에 바로 투입 가능하며, 최적화 과정에서 강성 제약(유효 탄성 텐서의 최소 고유값 ≥ 지정값)을 라그랑주 승수 혹은 페널티 함수로 구현한다.

수치 실험에서는 2차원 평면 응력 조건 하에 사각형, L‑형, 타원형 등 다양한 초기 형상을 설정하고, 첫 번째와 두 번째 밴드갭을 각각 최대화하는 두 개의 목표 함수를 정의하였다. 결과는 포함재의 부피 비가 증가할수록 밴드갭 하한이 낮아지지만, 형상 최적화를 통해 동일 부피에서도 밴드갭 폭을 크게 늘릴 수 있음을 보여준다. 특히, 강성 제약을 만족하면서도 연성 포함의 “스파이크” 형태가 저주파 밴드갭을 크게 확장시키는 최적 형상으로 수렴하였다. 이는 유효 질량 텐서의 최소 고유값이 더 넓은 주파수 구간에서 음의 값을 갖게 함으로써, 파동 전파를 억제하는 메커니즘과 일치한다.

이 논문은 전통적인 Bloch‑Floquet 기반 밴드갭 분석이 요구하는 전체 브릴루앙 영역 탐색을 회피하고, 동질화 모델만으로도 저주파 영역의 밴드갭을 정확히 예측·최적화할 수 있음을 증명한다. 또한, 형상 파라미터에 대한 민감도 도출 과정을 체계화함으로써, 고차원 설계 변수와 복합 제약을 동시에 다루는 대규모 최적화 문제에도 적용 가능함을 시사한다.