⚡ electrical engineering

A Non-Intrusive Full-3D Intrinsic Modal Formulation for Geometrically Nonlinear Structural Dynamics

본 논문은 요소 수준의 코드나 비선형 유한 요소 해에 대한 접근 없이도 3D, 쉘 및 빔 서브도메인의 원활한 조립을 가능하게 하며, 회전되지 않은 연속체 방정식으로부터 차수 축소 모델을 도출하는, 기하학적 비선형 구조 역학을 위한 비침습적 전방향 3D 고유 모드 정식화 방식을 제시한다.

원저자: Leonardo Buttà, Francesco Saltari, Michele Pasquali, Franco Mastroddi

게시일 2026-09-07

📖 4 분 읽기☕ 가벼운 읽기

CC BY 4.0

원저자: Leonardo Buttà, Francesco Saltari, Michele Pasquali, Franco Mastroddi

원본 논문은 CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/) 라이선스로 제공됩니다. ✨ 이것은 아래 논문에 대한 AI 생성 설명입니다. 저자가 작성하거나 승인한 것이 아닙니다. 기술적 정확성을 위해서는 원본 논문을 참조하세요. 전체 면책 조항 읽기

현대 항공기들은 점점 더 가늘고 유연해지고 있으며, 비행 중 발생하는 압력에 따라 날개가 휘고 뒤틀리며 하늘을 효율적으로 활공하도록 설계되고 있습니다. 전통적인 공학 모델은 강체 구조에는 잘 작동하지만, 날개가 너무 많이 휘어서 그 형상 자체가 힘이 작용하는 방식을 변화시킬 정도의 상황에서는 한계에 부딪힙니다. 이러한 극한 조건에서 구조 역학의 수학은 믿기 힘들 정도로 복잡해집니다. 엔지니어들은 종종 유연한 날개의 거동을 예측하기 위해 전체 컴퓨터 시뮬레이션에 의존하지만, 이러한 시뮬레이션은 계산량이 너무 많아 새로운 항공기를 제작하는 초기 단계에서 필요한 신속한 테스트와 설계 조정을 수행하기에는 부적합합니다. 따라서 정밀한 전체 시뮬레이션의 정확도를 희생하지 않으면서도, 이러한 비선형적 거동을 예측할 수 있는 더 단순하고 빠른 방법이 절실히 필요합니다.

이 과제를 이해하려면 먼저 재료가 어떻게 변형되는지를 파악해야 합니다. 고체가 휘거나 늘어날 때, 그 내부의 모든 미세한 점들은 움직입니다. 고전 물리학에서 재료의 한 점이 회전하는 것은 대개 그 점이 얼마나 늘어나는지와 직접적으로 연결되어 있어, 물체가 격렬하게 회전할 때 수학적으로 이를 분리해내는 것이 매우 어렵습니다. '고유 역학(intrinsic dynamics)'이라 불리는 다른 접근 방식은 회전과 신축을 별개의 독립적인 요소로 취급합니다. 이 방식은 엔지니어가 공간상의 모든 개별 점의 정확한 위치를 추적할 필요 없이, 오직 속도와 작용하는 힘만을 사용하여 구조의 운동을 기술할 수 있게 해줍니다. 이 방법은 단순한 보(beam) 형태의 구조물에는 매우 성공적으로 적용되어 왔으나, 보나 평판이라는 좁은 정의에 부합하지 않는 복잡한 3차원 형상에 적용될 때는 역사적으로 실패해 왔습니다.

이 연구에서 로마 사피엔차 대학교의 연구팀은 이 강력하고 단순화된 접근 방식을 어떤 복잡한 형상을 가진 3차원 물체에도 적용할 수 있도록 확장하는 새로운 수학적 프레임워크를 개발했습니다. 그들은 회전을 벗겨내어 순수한 변형을 드러내는 과정처럼, 회전을 분리하여 일반적인 고체의 운동을 기술하는 방법을 만들어냈습니다. 이렇게 함으로써, 그들은 변위와 거대한 회전이라는 번거로운 변수들을 완전히 배제하고, 오직 속도와 내부 힘만을 사용하여 구조의 거동을 설명하는 일련의 방정식들을 도출했습니다. 결정적으로, 이 새로운 공식은 '연역적(deductive)'입니다. 즉, 처음부터 해당 물체가 보인지 판인지 가정하지 않습니다. 대신, 고체 블록에 대한 일반적인 규칙에서 시작하여, 만약 보나 판에 대한 규칙을 적용한다면 방정식이 자연스럽게 알려진 형태의 정확한 식으로 단순화됨을 보여줍니다. 이는 동일한 수학적 엔진이 하나의 구조물 안에서 고체 블록, 얇은 껍데기(shell), 또는 긴 보, 혹은 이 세 가지가 혼합된 형태를 모두 처리할 수 있음을 의미합니다.

연구진은 이 복잡한 방정식들을 훨씬 작고 다루기 쉬운 모델로 축소했습니다. 그들은 구조물의 단순한 선형 분석으로부터 계산된 '모드(modes)'라고 불리는 표준 진동 패턴 위로 운동을 투영함으로써 이 작업을 수행했습니다. 이 단계가 바로 이 방법이 '비침입적(non-intrusive)'인 이유입니다. 즉, 엔지니어들이 원래의 복잡한 모델을 구축하는 데 사용된 소프트웨어의 깊숙한 독점 코드를 파헤칠 필요가 없습니다. 대신, 이 방법은 모든 구조 분석 소프트웨어가 생성하는 표준 출력값인 진동 모드의 형상과 그에 수반되는 힘만을 필요로 합니다. 이 제한된 데이터를 바탕으로, 연구팀은 전체 구조의 핵심적인 비선형 거동을 포착하는 새롭고 압축된 모델을 구축했습니다. 이 모델은 구동 속도가 매우 빨라 수천 번을 실행할 수 있으므로, 예비 항공 공학에서 사용되는 신속한 설계 루프에 적합합니다.

그들의 방법이 효과적임을 증명하기 위해, 연구팀은 13만 개 이상의 세부 지점을 가진 상세한 운송용 항공기 날개 컴퓨터 모델을 대상으로 테스트를 진행했습니다. 연구진은 이 가상의 날개에 점진적으로 하중을 가하여 비행 중 발생하는 극심한 굽힘 현상을 시뮬레이션했습니다. 그리고 그들은 이 새로운 빠른 모델의 결과를 정확도의 황금 표준으로 간주되는 고성능의 무거운 시뮬레이션 결과와 비교했습니다. 새로운 모델은 날개가 전체 길이의 거의 40%에 달할 정도로 휘어졌을 때조차도, 끝단 변위를 약 3~5%의 오차 내로 매우 정밀하게 추적하며 움직임을 예측했습니다. 또한 연구진은 날개 내부의 힘을 점검하였으며, 단순화된 모델이 복잡한 3차원 응력 패턴을 완벽에 가깝게 재현하여 무거운 시뮬레이션의 각 항목과 일치함을 확인했습니다.

본 연구는 모델의 복잡성을 줄이는 과정에서 원래 구조의 디테일을 얼마나 유지하느냐에 따라 모델의 정확도가 달라짐을 입증했습니다. 기본적인 방식으로 모델의 복잡성을 줄였을 때는 결과가 양호했으나 완벽하지는 않았습니다. 그러나 원래 구조의 동적 거동을 더 많이 보존하는 더 정교한 축소 기법을 사용했을 때, 새로운 모델은 무거운 시뮬레이션과 거의 구별할 수 없을 정도로 정밀해졌습니다. 이는 이 방법이 단순히 대략적인 근사치가 아니라, 기존 기술의 엄격한 일반화임을 확인시켜 줍니다. 이 방법은 단순한 보 모델의 속도와 전체 3차원 시뮬레이션의 정확성 사이의 간극을 성공적으로 메움으로써, 차세대 유연하고 고성능인 항공기를 설계하기 위한 강력한 새로운 도구를 제공합니다.

기술 요약: 기하학적 비선형 구조 동역학을 위한 비침입적 전방향 3D 고유 모드 정식화(Full-3D Intrinsic Modal Formulation)

문제 정의
현대 항공 및 우주 구조물은 점점 더 가늘고 유연해지고 있으며, 기하학적 비선형성이 유의미한 변형을 겪고 있다. 기존의 선형 모델은 이러한 영역에서 실패하는 반면, 완전 비선형 유한요소(FE) 시뮬레이션은 최적화, 안정성 분석, 트림 스윕(trim sweeps)과 같이 다수의 쿼리가 필요한 설계 작업에서 계산 비용이 너무 높다. 기존의 차수 축소 모델(ROM)은 이분법적 한계에 직면해 있다: 고유 빔 이론(예: Hodges, Palacios, Cea)은 정확한 대변형 회전 구조를 보존하고 일정한 비선형 계수를 산출하지만, 국부적인 3차원 유연성을 무시한 채 1차원 운동학에 국한된다. 반대로 STEP/ICE와 같은 비침입적 방법은 임의의 격자를 처리할 수 있지만, 비선형 정적 FE 솔루션을 통한 오프라인 시스템 식별에 의존하며, 이는 계산 비용이 많이 들고 로드 케이스에 따라 정확도가 달라지는 문제가 있다. 따라서 고유 빔 이론의 정확한 대회전 구조를 유지하면서도, 요소 수준의 코드나 비선형 FE 솔버에 접근할 필요 없이 일반적인 3차원 연속체에 적용 가능한 압축된 비선형 ROM이 필요하다.

방법론
저자들은 다음과 같은 단계를 통해 연속체 역학과 스펙트럼 축소를 연결하는 연역적 프레임워크를 제안한다:

고유 연속체 정식화(Intrinsic Continuum Formulation): 일반적인 3D 연속체의 역학을 준 극분해(quasi-polar decomposition)인 $F = \Lambda H$ 를 사용하여 회전되지 않은(unrotated) 형상에서 정식화한다. 여기서 $\Lambda$ 는 유한 회전을 일반화된 신장(stretch) $H$ 와 분리하는 독립적인 자세 필드(proper-orthogonal)이다. 이러한 디커플링을 통해 강체 회전이 없는 형상에서 응력-변형률 공액(stress-strain conjugates)을 정의할 수 있다.
혼합 변분 기술(Mixed Variational Statement): 운동 방정식은 Hu–Washizu 유형의 혼합 변분 원리로부터 유도된다. 속도/스핀, 일반화된 신장, 응력이 독립적인 필드로 취급된다. 결정적으로, 각운동량의 평형은 응력 대칭성에 대한 사전 가정으로 부과되는 것이 아니라, 독립적인 스핀 필드에 대한 정지 조건으로서 변분적으로 회복된다. 이는 변위 및 회전 변수 없이 속도와 변형률에 관한 1차 방정식 세트를 생성한다.
차원 축소(Dimensional Reduction): 이 정식화는 동일한 일반적 약형(weak form) 내에서 운동학적 제약 조건(예: 빔의 강체 단면 또는 판의 강체 두께 방향 섬유)을 부과함으로써 차원 축소에 대해 연역적이다. 이를 통해 클래식한 고유 빔 및 판 이론이 새로운 평형 법칙을 상정하지 않고도 특수한 사례로서 회복된다.
스펙트럼 투영(Spectral Projection): 혼합 약형은 회전되지 않은 형상 위에 정의된 선형 모드 기저로 투영된다. 이를 통해 스펙트럼 좌표( $\nu_m$ 속도형 및 $q_m$ 변형률형)에 대한 $2N$ 개의 비선형 상미분 방정식 시스템이 도출된다. 기하학적 및 자이로스코픽 비선형성은 상수 형태의 3차 결합 텐서( $\Gamma_1, \Gamma_2, \Gamma_3$ )로 수렴한다.
비침입적 조립(Non-Intrusive Assembly): 축소된 계수들은 표준 선형 고유 모드 해석의 출력값만을 사용하여 완전히 계산된다. 스펙트럼 형상(속도, 스핀, 응력, 변형률, 곡률)은 모드 형상과 요소 수준의 응력/힘 회복으로부터 유도된다. 이 방법은 요소 정식화나 비선형 FE 솔버에 대한 접근을 요구하지 않는다. 또한 이 방법은 서로 다른 하위 도메인(3D 고체, 쉘, 빔)이 독립적으로 이상화되어 가산적으로 조립되는 공간적 분할 모델을 지원한다.

주요 기여

전방향 3D 정식화: 본 논문은 고유 빔 이론의 정확한 기하학적 회전 구조를 유지하는 일반적인 3차원 연속체를 위한 최초의 고유(intrinsic, 속도-변형률) 정식화를 제시한다.
연역적 차원 축소: 고유 빔 및 판 이론이 별개의 가정이 아니라 단일 부모 연속체 진술의 운동학적 제한임을 입증한다.
비침입적 3D ROM: 선형 모드 데이터만을 사용하여 임의의 3D 이산화 모델(복잡한 윙박스 포함)에 대한 비선형 ROM을 구축할 수 있게 함으로써, 식별 기반 방법에서 요구되는 비선형 정적 솔브의 계산 비용을 피한다.
통합 조립: 혼합 차원 이상화(3D, 2D, 1D)로 구성된 구조물을 단일 스펙트럼 프레임워크 내에서 비침입적으로 조립할 수 있게 한다.

결과
본 방법론은 131,982개의 자유도를 가진 쉘 요소로 모델링된 대표적인 유연 수송기 날개(uCRM)를 통해 검증되었다.

계수 검증: 전체 3D 쉘 이산화를 통해 조립된 기하학적 결합 텐서( $\Gamma_2$ )는 공통 주 모드 영역에서 기존의 고유 빔 정식화와 0.97~0.98의 피어슨 상관 계수를 보였으며, 이는 3D 경로가 강체 단면 가정이 성립하는 지점에서 빔 결과를 재현함을 확인시켜 준다.
정적 비선형 검증: 설계 한계치의 최대 5배에 달하는 추종 공력(follower aerodynamic loads) 하에서, ROM은 기하학적 비선형 FE 솔루션(NASTRAN SOL400)과 매우 유사한 평형 경로(팁 변위 최대 38% 반폭)를 재구성하였다. 팁 변위 오차는 3%에서 5% 사이였으며, 응력장 잔차는 1.7% 미만이었다.
동적 충실도: 변형된 평형 상태에 대한 접선 모드 해석 결과, '완전한'(비축소) 및 Craig–Bampton 축소 기저 모두 참조 모델과 비교하여 고유 진동수 및 모드 형상(교차 직교성)에서 높은 충실도를 유지했다. 반면, Guyan–Irons 축소는 변형이 증가함에 따라 성능이 크게 저하되었으며, 이는 대변형 시 잔류 동적 콘텐츠를 유지하는 것의 중요성을 강조한다.

의의 및 주장
본 논문은 비선형 고유 모드 동역학의 일관된 일반화를 제공한다고 주장한다. 일반적인 3D 연속체로부터 모델을 유도함으로써, 본 연구는 고유 스펙트럼 방법의 적용 범위를 빔 형태의 구조를 넘어 국부적 유연성(워핑, 패널 굽힘)이 중요한 복잡한 3차원 기하학적 구조로 확장한다. 이 접근 방식은 정확한 기하학적 비선형 구조를 보존하면서도 전체 비선형 FE 해석에 대한 계산 효율적인 비침입적 대안을 제공한다는 점에서 중요하다. 저자들은 현재 구현이 공간적 형상 회복을 위해 빔 형태의 골격(skeleton)에 의존하고 있으나, 프레임워크 자체는 축소 없이 완전한 3D 스펙트럼 모델을 수용할 수 있어 전통적인 빔 기반 ROM의 범위를 벗어난 국부 좌굴 및 비선형 분기(bifurcation) 현상을 연구할 수 있는 길을 열어준다고 언급한다.

유사한 논문