Size-dependent elastostatic behavior of non-centrosymmetric solids via a mixed variational finite element method

기술 요약: 혼합 변분 유한 요소법을 이용한 비중심대칭 고체의 크기 의존적 탄성 정역학적 거동

문제 정의
고전 탄성학 및 일관적 커플 응력 이론(C-CST)은 중심대칭 재료에 광범위하게 적용되어 왔으나, 비중심대칭 고체의 크기 의존적 탄성 정역학적 응답은 C-CST 프레임워크 내에서 거의 탐구되지 않았다. 압전 현상과 같은 현상으로 알려진 비중심대칭 결정 구조는 힘-응력(force-stress)과 변형률(strain) 사이, 그리고 커플-응력(couple-stress)과 평균 곡률(mean curvature) 사이를 결합하는 3차 구성 텐서를 갖는다. 본 논문은 마이크로 및 나노 스케일에서 크기 의존적 효과가 중요해지는 비중심대칭 재료를 분석하기 위한 계산 도구와 이론적 정식화의 부재를 다룬다. 주요 과제는 변위장의 고차 연속성을 요구하지 않으면서도 힘-응력, 커플-응력, 변형률, 그리고 공학적 평균 곡률 사이의 열역학적 및 연속체 역학적 상호작용을 완전히 설명할 수 있는 정식화를 개발하는 것이다.

방법론
저자들은 C-CST에 기반한 새로운 혼합 변분 유한 요소법(FEM)을 개발하였다. 방법론은 다음 단계로 진행된다:

1. 이론적 정식화: 구동 방정식은 준정적(quasistatic) 조건과 미소 변형을 가정하여 C-CST로부터 유도된다. 구성 관계는 압전 현상과의 유사성을 이용하여 일반화된 헤크만(Heckmann) 도표 접근법을 통해 설정된다. 이는 힘-응력 및 커플-응력을 각각 변형률 및 평균 곡률과 결합하는 3차 구성 텐서($Q$ 또는 $P$)를 도입한다. 이 정식화는 강성 기반 텐서와 컴플라이언스(compliance) 기반 텐서를 구분하며, 다양한 실험적 제약 조건(예: 일정 변형률 vs 일정 응력) 하에서의 관계를 정의한다.
2. 혼합 변분 원리: 변위의 2차 미분을 포함하는 표준 변위 기반 정식화에서 요구되는 $C^1$ 연속성의 계산적 어려움을 피하기 위해, 저자들은 혼합 변분 원리를 유도한다. 이 원리는 변위($u_i$)와 커플-응력($\mu_m$)을 주요 변수로 취급한다. 운동학적 제약을 만족시키기 위해 라그랑주 승수(Lagrange multipliers)를 도입함으로써, 저자들은 "정적 혼합 퍼텐셜 에너지 원리(Principle of Stationary Mixed Potential Energy)"를 유도한다. 이는 변위 및 커플-응력장에 대한 연속성 요구를 $C^0$로 감소시킨다.
3. 유한 요소 구현: 2차원 평면 변형(plane strain) 이산화가 저차 삼각형 요소를 사용하여 구현된다.
   • 변위: 표준 3노드 선형 형상 함수가 사용된다.
   • 커플-응력: 커플-응력 벡터를 표현하기 위해 저자들은 에지(edge) 기반 형상 함수를 채택한다. 중심대칭 사례에서는 커플-응력 벡터가 발산이 0(divergence-free)이지만, 비중심대칭 재료는 일반적으로 요소 에지를 따라 접선 및 법선 성분을 모두 필요로 한다. 이는 에지당 두 개의 자유도를 발생시킨다.
   • 수치 해법: 약형(weak form)은 대칭 및 반대칭 응력 성분을 포함하는 선형 대수 방정식 시스템으로 이산화된다. 시스템은 MATLAB의 희소 행렬(sparse matrix) 기법을 사용하여 해결된다.

주요 기여

• 비중심대칭 고체에 대한 첫 적용: 본 연구는 C-CST를 사용하여 비중심대칭 고체의 크기 의존적 탄성 정역학적 거동을 계산적으로 탐구한 첫 사례이다.
• 새로운 혼합 변분 프레임워크: 이 논문은 (중심대칭 재료에 대해 이전에 사용되었던) 혼합 변분법을 쌍방향 이방성(비중심대칭) 사례로 확장하여, 힘/커플-응력과 변형/곡률 사이의 완전한 결합을 성공적으로 처리하였다.
• 구성 모델링: 저자들은 비중심대칭 재료의 열역학적 관계와 구성 텐서를 명시적으로 정의하며, 테스트 중에 일정하게 유지되는 변수에 따라 재료 매개변수가 어떻게 달라지는지 명확히 한다.
• 요소 개발: 이 정식화는 일반적인 비중심대칭 대칭성에 필수적인 커플-응력의 접선 및 법선 성분을 모두 포착할 수 있는 특정 삼각형 에지 요소를 도입한다.

결과
방법론은 세 가지 비중심대칭 재료(입방정계 탄화규소(SiC), 육방정계 SiC, 삼방정계 석영)를 포함하는 네 가지 계산 예제를 통해 검증되었다. 문제에는 굽힘 하의 프리즘형 캔틸레버, 축하중 하의 프리즘형 바, 전단 하의 깊은 보, 그리고 복합 하중 하의 곡선형 보가 포함되었다.

• 수렴성: 혼합 FEM 해는 우수한 수렴성을 보였다. 상대적 에너지 노름은 굽힘 문제에서는 1보다 약간 낮은 수렴 속도를, 축 하중 문제에서는 거의 1에 가까운 수렴 속도를 보였으며, 이는 선형 형상 함수의 사용과 일치한다.
• 크기 의존적 강성: 중심대칭 재료에서 관찰된 바와 같이, 무차원 고유 길이 척도($l/h$)가 증가함에 따라 구조물은 더 단단해졌으며, 결국 포화 영역에 도달했다. 굽힘 지배 문제(예: 캔틸레버)는 상당한 강성 증폭(최대 ~43배)을 보인 반면, 전단 지배 문제(깊은 보)는 최소한의 증폭(<2배)을 보였다.
• 비중심대칭 효과:
  • 결합 매개변수 ($l_p$): 비중심대칭 결합 텐서와 관련된 두 번째 길이 척도 매개변수($l_p$)는, 특히 $l_p$가 $2l$에 접근함에 따라 강성을 더욱 증가시키는 것으로 나타났다.
  • 이방성: 응답은 결정 방향과 하중 평면에 따라 크게 의존했다. 예를 들어, 육방정계 SiC의 경우 $x-y$ 평면에서는 축 하중과 횡방향 곡률 사이의 결합이 유의미했으나, $x-z$ 평면에서는 무시할 만했다.
  • 선형 횡방향 변위: 축 압축 하중 하에서 비중심대칭 재료는 임계 좌굴 하중 없이 선형 횡방향 변위(굽힘)를 나타낸다는 것이 중요한 발견이다. 이는 결합 매개변수의 부호에 따라 변위의 부호가 결정되는 역압전 효과와 유사한 순수 선형 역학 효과이다.
  • 응력 비대칭성: 힘-응력 텐서의 반대칭 부분이 0이 아니므로, 불균일한 전단 응력 성분($\sigma_{xy} \neq \sigma_{yx}$)과 고전 탄성학과는 구별되는 독특한 응력 분포를 초래한다.

의의 및 주장
본 논문은 비중심대칭 탄성 고체의 크기 의존적 거동을 분석하기 위한 최초의 계산 프레임워크를 제공한다고 주장한다. 저자들은 축 하중과 횡방향 곡률 사이의 관찰된 선형 결합이, 규명되지 않은 커플-응력 구성 매개변수($l_p$ 및 $P$ 텐서 성분)를 측정하기 위한 물리적 실험의 잠재적 메커니즘을 제공한다고 강조한다. 연구는 얕은 캔틸레버 굽힘 테스트가 고유 길이 척도 $l$을 추정하는 데 효과적인 반면, 축 하중 테스트는 비중심대칭 결합 매개변수를 분리해낼 수 있음을 시사한다. 저자들은 현재의 작업이 선형 탄성 정역학에 집중하고 있지만, 개발된 프레임워크가 진동, 과도 동역학, 3차원 문제, 그리고 중심대칭 결정에서의 잠재적인 대변형 유도 대칭성 깨짐을 포함한 대변형 연구를 위한 토대를 마련한다고 결 결론짓는다.