Homogenization of rod-like metamaterials as a special Cosserat rod

본 논문은 단순한 격자 구조부터 복잡한 음의 푸아송비 구조 및 인공 근육 구조에 이르기까지 다양한 수치 예시를 통해 검증된, 주기적으로 조립된 막대형 메타물질의 비선형 구성 응답과 강성을 유도하기 위해 기하학적으로 정확한 특수 코시세 막대 이론에 기반한 변분 동질화 기법을 제시한다.

핵심

상상해 보세요. 단일 고체 물질로 만들어진 긴 유연한 튜브가 아니라, 미세한 막대들이 연결된 복잡하고 반복적인 패턴으로 이루어진 튜브를 말입니다. 마치 원통으로 감긴 미세한 사다리나 사슬 링크 울타리와 같은 구조죠. 저자들은 이를"막대형 메타물질"이라고 부릅니다.

저자들이 해결하려는 문제는 다음과 같습니다. 이 긴 튜브 전체가 어떻게 구부러지고, 늘어나고, 비틀리는지 알고 싶다면, 단순히 하나의 작은 막대만 보면 안 됩니다. 수천 개의 막대로 이루어진 전체 네트워크가 어떻게 상호작용하는지 살펴봐야 합니다. 긴 튜브의 모든 단일 막대를 시뮬레이션하는 것은 해변의 모든 모래 알갱이를 세어 바람에 따라 해변이 어떻게 움직이는지 이해하려는 것과 같습니다. 이는 너무 많은 컴퓨터 성능을 필요로 합니다.

저자들은 전체 튜브의 거동을 예측하기 위해 아주 작은 대표 조각만 연구하는 똑똑한 방법, 즉"레시피"를 제안합니다. 간단한 비유로 설명해 보겠습니다.

1."마법 같은 줌인"(동질화)

메타물질을 거대한 반복 벽지 패턴이라고 생각해 보세요. 벽 전체를 분석하는 대신 벽지 한 칸만 보면 됩니다 (이를 RVE, 즉 대표 부피 요소라고 합니다).

저자들의 비법은 긴 튜브 전체를 당기거나 비틀면, 그 작은 사각형 조각도 매우 특정한 나선형 방식으로 당겨지거나 비틀린다고 가정하는 것입니다. 이를나선형 변형이라고 부릅니다. 스프링을 잡아당기는 것을 상상해 보세요. 코일들이 단순히 길어지는 것뿐만 아니라 약간 회전합니다. 저자들은 이 작은 사각형 조각이 그 정확한 나선 운동을 모방하도록 강제함으로써, 전체를 시뮬레이션하지 않고도 긴 튜브 전체가 어떻게 반응할지 알아낼 수 있음을 깨달았습니다.

2."완전히 유연한"막대들

대부분의 컴퓨터 모델은 막대를 강철 자처럼 뻣뻣하고 변하지 않는 것으로 취급합니다. 하지만 실제 생활, 특히 이러한 미세 메타물질에서는 막대가 동시에 구부러지고, 늘어나며, 전단 (옆으로 미끄러짐) 될 수 있습니다. 변형이 거대할 때조차도 그렇습니다.

저자들은특수 코시레트 막대라는 특수한 수학적 모델을 사용합니다.

• 비유: 삶은 스파게티 조각을 상상해 보세요. 구부러질 수도 있고, 약간 늘어날 수도 있으며, 비틀릴 수도 있습니다. 이제 그 스파게티가 원형으로 구부러지거나 길이의 두 배로 늘어날 때조차도 이러한 모든 일을완벽하게그리고정확하게해낼 수 있는 재료로 만들어졌다고 상상해 보세요. 이것이 바로 그들의 모델이 하는 일입니다. 단순히 근사하는 것이 아니라, 굽힘과 비틀림의 정확한 기하학을 포착합니다.

3."댄스 플로어"규칙 (경계 조건)

작은 사각형 조각이 거대한 반복 튜브의 일부처럼 행동하도록 하기 위해, 저자들은 그 사각형의 가장자리들이 서로 어떻게 소통하는지에 대한 일련의 규칙을 고안해야 했습니다.

• 문제: 나선형 계단의 한 조각을 잘라내면, 위쪽 가장자리가 아래쪽 가장자리와 완벽하게 정렬되지 않습니다.
• 해결책: 그들은"나선형 경계 조건"을 만들었습니다. 작은 사각형 조각의 왼쪽 가장자리가 오른쪽 가장자리와 손을 잡고 있다고 상상해 보세요. 하지만 오른쪽은 나선형 계단의 단계처럼 약간 회전하고 이동되어 있습니다.
• 혁신: 이전 방법들은 작고 부드러운 움직임만 다룰 수 있었습니다. 저자들의 새로운 규칙은 튜브가 프렌치 도넛처럼 비틀리거나 실처럼 얇아질 때까지 늘어날 때에도 작동합니다. 이는"기하학적으로 정확"한 것으로, 모양이 얼마나 기괴해지더라도 정확도를 잃지 않습니다.

4."조인트"와"접착제"

그 작은 사각형 조각 내부에서 막대들은 조인트에서 연결됩니다.

• 강성 조인트: 일부 조인트는 초강력 접착제와 같습니다. 연결점에서 막대들은 서로 상대적으로 움직일 수 없습니다.
• 수학: 저자들은 컴퓨터가 최소한의 에너지만 사용하면서 그 작은 사각형 내의 모든 막대의 최상의 위치를 계산하도록 시스템을 구성했습니다. 이때 조인트가 연결된 상태를 유지하고"나선형 계단"규칙이 준수되도록 합니다.

5. 그들이 발견한 것 (결과)

작은 조각에 대한 수학을 푼 후, 그들은 전체 튜브가 어떻게 행동할지 예측할 수 있었습니다. 그들은 다양한 모양으로 이를 테스트했습니다.

• 십자형과 사각형 모양: 그들은 수학이 작동하는지 증명하기 위해 간단한 모양 (막대로 만든 플러스 기호나 사각형) 으로 시작했습니다. 그들은 짧은 두꺼운 막대라면 늘어남이나 전단 여부가 매우 중요하다는 것을 발견했습니다. 반면 매우 얇고 긴 막대라면 기존의 더 간단한 수학으로도 충분했습니다.
• 나선형 (스프링) 막대: 그들은 이미 스프링처럼 휘어진 막대 (나선) 로 만든 사각형을 살펴보았습니다.
  • "J 자형"늘어남: 이 재료를 당겼을 때, 처음에는 스프링이 풀리는 것처럼 부드럽다가 곧게 펴지면서 매우 뻣뻣해졌습니다. 이로 인해"J 자형"곡선이 생성됩니다. 이는 근육과 같은 생물학적 조직이 행동하는 방식과 정확히 일치하므로, 저자들은 이것이인공 근육에 사용될 수 있다고 언급했습니다.
  • 연화되는 굽힘: 구부렸을 때, 재료가 구부러질수록 더 부드러워졌습니다. 이는 연결된 스프링 막대가 평면에서 비틀려 힌지처럼 작용했기 때문입니다.
• 음의 푸아송 비 튜브: 그들은 당기면 넓어지는 (벌집 구조처럼) 중공 튜브를 모델링했습니다.
  • 그들은 막대의 각도를 변경함으로써 튜브를 측면으로는 매우 유연하게 (구부리기에 좋음) 만들면서 압축에는 매우 뻣뻣하게 (혈관을 열어두기에 좋음) 만들 수 있음을 보였습니다.
  • 그들은 이러한 구조들을"단축"(확장될 때 짧아지는 현상) 을 피하도록 조정할 수 있다고 지적했습니다. 이는 동맥을 지탱하는 메쉬 튜브인심혈관 스텐트에서 흔한 문제입니다.

요약

저자들은 메타물질을 위한"보편적 번역기"를 구축했습니다. 그들은 복잡한 3 차원 막대 네트워크를 단일 막대의 단순하고 매끄러운 수학적 설명으로 변환하는 방법을 개발했습니다. 이를 통해 엔지니어들은 미세한 내부 패턴을 조정하여로봇 팔, 인공 근육, 의료용 스텐트와 같은 복잡하고 유연한 재료를 설계할 수 있게 되었습니다. 최종 제품이 어떻게 구부러지고 늘어날지 정확히 알 수 있으므로, 모든 설계 변경 사항마다 슈퍼컴퓨터 시뮬레이션을 실행할 필요가 없습니다.

기술 요약

기술 요약: 특수 코시 rod 로서의 막대형 메타물질 균질화

문제 제기
막대형 메타물질은 미세구조 단위 (일반적으로 막대들의 네트워크) 를 단일 방향으로 주기적으로 조립하여 형성된 구조로, 한 차원의 거시적 기하학이 나머지 두 차원에 비해 현저히 길어지는 특징을 가집니다. 이러한 구조는 심혈관 스텐트와 인공 근육부터 로봇 구동기에 이르기까지 다양한 응용 분야에서 활용됩니다. 그러나 이러한 이질적 구조의 직접적인 전체 규모 수치 시뮬레이션은 계산 비용이 너무 많이 들어 실행이 불가능한 경우가 많으며, 표준 균질화 기법들은 종종 한계에 직면합니다. 구체적으로, 1 차 균질화 체계는 막대형 구조에서 전형적으로 발생하는 큰 변형률 기울기 (예: 굽힘과 비틀림) 를 수반하는 큰 변형을 포착하기에 부족합니다. 또한, 횡방향에서 길이 척도의 분리가 성립하지 않을 때, 표준 3 차원 균질화 접근법은 크기 효과를 고려하지 못합니다. 따라서 미세 및 거시 규모 모두에서 유한 전단, 신장, 굽힘, 비틀림을 고려하여 이러한 메타물질의 비선형 구성 응답을 임의의 큰 변형 하에서 정확하게 모델링할 수 있는 프레임워크가 필요합니다.

방법론
저자들은 막대형 메타물질을 거시 규모에서 "특수 코시 rod"로 모델링하는 계산 균질화 체계를 제안하며, 미세구조 단위 (대표 체적 요소 또는 RVE) 는 기하학적으로 정확한 코시 막대들의 네트워크로 모델링합니다.

1. 기하학적으로 정확한 운동학: 거시적 유효 막대와 미시적 구성 막대 모두 기하학적으로 정확한 특수 코시 rod 이론을 사용하여 모델링됩니다. 이 공식은 고전적 빔 이론의 소변형 가정들을 피하면서 유한 굽힘, 비틀림, 전단, 그리고 축방향 신장을 고려합니다.
2. 나선형 코시 - 본 (HCB) 규칙 확장: 이 방법은 나선형 코시 - 본 규칙을 확장합니다. 이 방법은 메타물질이 거시 규모에서 아크 길이 따라 균일한 변형률 장을 받는다고 가정합니다. 이 가정은 전체 구조 문제를 단일 RVE 문제로 축소할 수 있게 합니다.
3. 나선형 주기적 경계 조건: 주요 혁신은 RVE 에 기하학적으로 정확한 비선형 매핑을 유도하여 나선형 주기적 경계 조건을 적용하는 것입니다. 이 매핑은 거시적 변형률 매개변수 ($v_0, k_0$) 를 기반으로 RVE 의 오른쪽 경계의 병진 및 회전 자유도를 왼쪽 경계와 연관시킵니다. 이전의 선형 전이 매핑과 달리, 이 공식은 임의의 큰 거시적 변형률에서도 유효합니다.
4. 제약된 에너지 최소화: RVE 문제는 탄성 에너지의 제약 최소화 문제로 공식화됩니다. 제약 조건에는 다음이 포함됩니다:
   • 내부 조인트 제약: 마스터 - 슬레이브 접근법을 사용하여 RVE 내 막대들 사이의 강체 연결을 강제합니다.
   • 나선형 경계 제약: HCB 규칙에서 유도된 주기성을 강제합니다.
   • 전역 제약: RVE 의 강체 병진 및 회전을 제한하여 고유한 해를 보장하고 최소화 동안 고정된 거시적 변형률 매개변수를 유지합니다.
5. 이산화 및 해법: RVE 에너지는 "나선형 요소"를 사용하여 이산화되며, 여기서 각 세그먼트는 균일한 변형률 장을 받는다고 가정하여 세그먼트 자체가 나선형이 됩니다. 결과적으로 생성된 비선형 약형식은 deal.II 플랫폼에서 자체 개발 코드를 사용하여 해결됩니다.
6. 거시적 특성 유도: 균질화된 막대의 내부 접촉력, 모멘트, 그리고 강성 텐서에 대한 폐쇄형 표현식은 최소화된 RVE 에너지 밀도를 거시적 변형률 매개변수에 대해 미분하여 유도됩니다.

주요 기여
본 논문은 세 가지 주요한 신규 기여를 식별합니다:

1. 이중 규모 기하학적 정확성: 미시 규모와 거시 규모의 막대 모두 기하학적으로 정확한 특수 코시 rod 로 모델링되어, 두 규모 모두에서 굽힘과 비틀림뿐만 아니라 유한 전단과 신장을 포착합니다.
2. 기하학적으로 정확한 거시 - 미시 매핑: 병진 및 회전 자유도에 나선형 주기적 경계 조건을 적용하기 위한 완전한 비선형 매핑이 제안됩니다. 이 매핑은 큰 거시적 변형률에서도 유효하여 이전의 선형 전이 매핑의 한계를 해결합니다.
3. 임의의 미시 규모 구성 거동: 이 프레임워크는 RVE 내 구성 막대들의 임의의 구성 거동을 허용하여 선형 탄성 가정을 넘어섭니다.

결과 및 수치 예시
저자들은 다양한 RVE 복잡성을 가진 여러 수치 예시를 사용하여 체계를 검증했습니다:

• 십자형 및 정사각형 RVE: 기존 분석 공식 및 문헌과 결과를 비교 검증하기 위해 간단한 평면 RVE(십자형 및 정사각형 형태) 가 사용되었습니다. 이 연구는 미시 규모 막대에서 전단과 신장을 포함하는 것이 거시적 강도에 특히 큰 세장비를 가진 막대의 경우 상당한 영향을 미친다는 것을 보여주었습니다. 결과는 키르히호프 (비신장/비전단) 막대에 대한 분석적 한계와 일치했습니다.
• 나선형 막대 RVE: 나선형 구성 막대로 구성된 더 복잡한 정사각형 RVE 가 분석되었습니다. 결과는 굽힘 지배적 거동에서 신장 지배적 거동으로 전환되는 J 자형 거시적 응력 - 신장 응답을 보여주었습니다. 이 연구는 나선형 회전 수와 피치와 같은 매개변수가 이러한 응답을 어떻게 조절할 수 있는지 강조했습니다. 굽힘 분석은 교차 연결 나선형 막대의 면외 굽힘에 의해 유발되는 연화 거동을 드러냈습니다.
• 음의 푸아송비 튜브형 구조: 이 체계는 스텐트와 관련된 음의 푸아송비 튜브형 메타물질에 적용되었습니다. 이 모델은 설계 각도 $\theta_a$에 따른 거시적 굽힘 강성, 푸아송 비 (단축 거동), 그리고 방사형 강성의 의존성을 성공적으로 포착했습니다. 또한 큰 굽힘 하에서 단면의 타원화 현상인 바지에 (Bazier) 효과를 포착하는 것을 보여주었습니다.
• FE2 검증: 제안된 균질화 체계 (FE2 접근법) 와 십자형 RVE 메타물질의 전체 규모 시뮬레이션 간의 비교는 선형 구성 모델이 실패하는 큰 비틀림 모멘트 하에서도 방법의 정확성을 검증하는 좋은 일치를 보여주었습니다.

의의 및 주장
본 논문은 1 차 및 선형 균질화 체계의 한계를 극복하는 막대형 메타물질 균질화를 위한 강력한 프레임워크를 제공한다고 주장합니다. 두 규모에서 기하학적으로 정확한 특수 코시 rod 이론을 활용하고 비선형 거시 - 미시 전이 매핑을 유도함으로써, 이 방법은 큰 변형 하에서 이러한 구조들의 비선형 구성 응답을 정확하게 예측합니다.

저자들은 그들의 접근법이 일반적이며 다른 RVE 유형 (예: 다공성 고체, 오리가미 구조) 으로 확장될 수 있으며, 탄소성, 점탄성, 그리고 다물리 결합 (열, 화학, 또는 전기 탄성) 과 같은 복잡한 미시 규모 현상을 통합할 수 있다고 강조합니다. 이 연구는 미세구조 매개변수를 조절하여 인공 근육이나 전개형 스텐트와 같은 특정 거시적 특성을 가진 맞춤형 메타물질을 설계하기 위한 기초 역할을 합니다. 저자들은 겸손하게도 현재 작업은 RVE 가 막대 네트워크라고 가정하지만, 유도된 경계 조건은 다른 미세구조 표현에 맞게 적응될 수 있다고 언급합니다.