Model-Informed Joint Material-Structural Optimization of Hard-Magnetic Soft Materials

기술 요약: 경자성 연성 재료의 모델 기반 통합 재료-구조 최적화

문제 정의
경자성 연성 재료(hard-magnetic soft materials, hMSMs)는 유연한 엘라스토머 매트릭스 내에 경자성 입자가 매립된 복합체이다. 이 재료들은 소프트 로보틱스 및 적응형 구조물에 적합한 접촉 없는 자기장 구동 변형을 가능하게 한다. 그러나 hMSM의 예측 가능한 설계를 위해서는 두 가지 상호 연결된 과제에 직면한다:

1. 구성 모델링(Constitutive Modeling): 정확한 예측을 위해서는 자기 입자의 부피 분율이 유효 탄성 계수를 어떻게 변화시키는지 설명하는 효과적인 구조-물성 관계가 필요하다. 기존 문헌에는 동일한 입자 농도에 대해 상당히 다른 강성 값을 예측할 수 있는 다양한 모델들(강성 입자 포함, 자기 일관 모델, 구속 운동학 모델 등)이 존재한다. 또한, 변형 에너지 밀도 함수의 선택은 수치적 견고성이나 예측된 변형에 영향을 미칠 수 있다.
2. 통합 최적화(Joint Optimization): 최적 설계에는 구조적 밀도, 자기 입자 분포, 그리고 잔류 자화 방향의 동시 제어가 필요하다. hMSM에 대한 기존의 위상 최적화 접근 방식은 자기 재료의 존재를 강성 및 구동 강도를 조절하는 연속적인 부피 분율 필드로 취급하기보다는, 자기화 여부(자화됨 vs 자화되지 않음)를 이산적인 선택으로 취급하는 경우가 많다. 또한, 자기 변수들은 구조적 재료가 존재하는 곳에서만 의미를 갖기 때문에 이러한 필드들을 물리적으로 결합하는 것은 어렵다.

방법론
저자들은 구성 모델 선택과 통합 최적화를 모두 다루는 통일된 프레임워크를 개발하였다.

• 통합 구성 프레임ik: 저자들은 고전적인 강성 입자 관계(Guth, Mooney, Kerner), Hill의 자기 일관 이론(강성 입자 및 근사 비압축성 매트릭스에 특화됨), 그리고 구속 운동학 모델(Lopez-Pamies, Alshammari)을 공통의 유효 전단 탄성률 형태인 $G(\phi)$ (여기서 $\phi$는 입자 부피 분율)로 재구성하였다. 이를 통해 7가지 후보 전단 탄성률 관계식을 도출하였다. 이들은 세 가지 변형 에너지 밀도 함수(두 가지 압축성 neo-Hookean 형태인 NH1, NH2와 St. Venant–Kirchhoff 모델인 SVK)와 결합되어 21개의 후보 구성 모델을 생성한다.
• 모델 민감도 분석: 캔틸레버 빔, 바퀴, 그리퍼, 3D 빔 등의 수치 벤치마크를 사용하여, 예측된 변형이 변형 에너지 함수의 선택과 $G(\phi)$의 선택에 따라 어떻게 민감하게 반응하는지 평가한다.
• 구성 모델 선택: Garcia-Gonzalez 등(2021)의 실험적 단축 응력-변형 데이터로부터 가장 대표적인 $G(\phi)$ 관계를 선택한다. 두 가지 지표가 사용된다: (1) 추출된 전단 탄성률 값에 대한 상대 오차, (2) 전체 실제 응력-변형 곡선에 대한 상대 평균 제곱근 오차(rRMSE).
• 통합 재료-구조 최적화 정식화: 세 가지 설계 필드가 매개변수화된 위상 최적화 프레임워크를 개발하였다: 구조적 밀도($\rho$), 자기 입자 부피 분율($\phi$), 그리고 잔류 자화 방향($\theta_{Br}$).
  • 결합 메커니면: 물리적 일관성을 보장하기 위해, 자기 에너지와 유효 계수는 자기 기여도가 공극 영역($\rho \approx 0$)에서 사라지도록 보간된다. 자기 입자 부피 분율 $\phi$는 선택된 $G(\phi)$ 관계를 통해 자기 상호작용 에너지와 유효 전단 탄성률을 모두 스케일링하는 연속적인 필드로 작용한다.
  • 구현: 본 프레임워크는 밀도에 대한 SIMP형 페널티 적용, 모든 필드에 대한 헬름홀츠 필터링, 그리고 Method of Moving Asymptotes (MMA) 최적화기를 위한 어드조인트 기반 그래디언트 계산을 사용하는 FEniTop 라이브러리와 FEnicsx 기반의 유한 요소 구현을 활용한다.

주요 결과

• 구성 민감도: 변형 에너지 밀도 함수(NH1, NH2, SVK)는 고려된 구동 문제에서 예측된 변형에 비교적 작은 영향을 미친다. 반면, 유효 전단 탄성률 관계 $G(\phi)$의 선택은 특히 자기 재료가 높은 변형 영역과 겹치는 경우 변형에 상당한 영향을 미친다. 자기 재료가 주요 굽힘 영역(예: 팁 자화 빔)에서 떨어져 국부화되어 있을 때는 모델 의존성이 미미하다.
• 모델 선택: 실험 데이터를 바탕으로, $G(\phi)$에 대한 Mooney 관계식이 추출된 전단 탄성률 값과 전체 응력-변형 응답 모두에서 가장 우수한 종합적 일치를 보였다. NH2 변형 에너지 함수는 수치적 견고성 측면에서 선택되었다.
• 최적화 성능: 통합 최적화 프레임워크는 세 가지 뚜렷한 사례에 대해 비직관적인 설계를 성공적으로 도출하였다:
  1. 회전 액추에이터: 고정된 바퀴 기하 구조에서 $\phi$와 $\theta_{Br}$를 최적화함으로써, 직관적인 균일 설계 대비 접선 변위를 59.7% 증가시켰다.
  2. 직선 운동 액추에이터: 고정된 가위 모양 기하 구조에서 $\phi$와 $\theta_{Br}$를 최적화하여, 수직 방향의 기생 변위를 억제하면서 목표한 수평 운동을 달성하였다.
  3. 복원 빔: $\rho$, $\phi$, $\theta_{Br}$를 동시에 최적화하여, 반대되는 기계적 하중 하에서도 자기 구동을 통해 원래 형상으로 복원할 수 있는 빔을 제작하였다. 이 설계는 자기 재료를 구조적으로 점유된 영역 내에 자연스럽게 국부화시켰다.

의의 및 주장
본 논문은 hMSM 설계를 위한 체계적이고 모델 기반의 경로를 제공한다고 주장한다. 주요 기여는 다음과 같다:

1. 다양한 구조-물성 관계를 직접 비교할 수 있는 통합 유효 전단 탄성률 프레임워크를 제공한다.
2. 자기 및 기계적 필드가 겹칠 때 예측 정확도를 위해 유효 전단 탄성률 관계가 결정적인 반면, 변형 에너지 형태는 테스트된 사례들에서 덜 민감하다는 증거를 제시한다.
3. 자기 입자 부피 분율을 구조적 밀도 및 자화 방향과 결합된 연속적인 설계 변수로 취급하고, 에너지 보간을 통해 저밀도 영역에서의 물리적 결합 문제를 해결하는 통합 재료-구조 최적화 정식화를 제안한다.
4. 재료 조성과 구조적 위상을 동시에 최적화하는 것이 고정된 기하 구조의 자기 설계보다 우수한 성능을 낸다는 것을 입증한다.

저자들은 본 프레임키가 hMSM에 대해 시연되었으나, 관련 에너지 항을 교체함으로써 다른 분야 반응성 엘라스토머(예: 유전체, 자기 유변 또는 형상 기억 복합체)에도 적용 가능한 계산 구조를 가지고 있음을 언급하였다. 재현성과 확장을 지원하기 위해 구현체는 오픈 소스 CEADpx/top_optim 저장소에 공개되어 있다.