마찰 구조물 강제 응답 백본 곡선 정밀 계산을 위한 해시안 텐서 기반 연속법

초록

본 논문은 마찰 접촉을 포함한 비선형 구조물의 강제 진동 백본 곡선을 고정밀로 추적하기 위해, 라그랑주 승수와 파라미터 연속법을 결합하고 접촉 요소의 해시안 텐서를 해석적으로 도출한 새로운 계산 프레임워크를 제안한다. Euler‑Bernoulli 빔과 블레이드‑댐퍼‑블레이드 모델을 통해 기존 dNNM·위상 공명 방법보다 높은 정확도와 강인성을 입증한다.

상세 분석

이 연구는 마찰 비선형성을 가진 구조 시스템의 강제 응답을 해석하기 위해 세 가지 핵심 기술을 통합한다. 첫째, 다중 고조파 균형법(HBM)을 기반으로 비선형 접촉력을 주파수 영역에서 직접 표현한다. 기존에는 수치 미분을 이용해 야코비안을 계산했으나, 저자들은 Petrov의 접촉 모델을 확장해 접촉 요소의 2차 미분인 해시안 텐서를 정확히 유도하였다. 이 해시안 텐서는 접촉력의 변곡점을 정밀히 포착하고, 뉴턴‑라프슨 반복 과정에서 스텝 크기를 크게 할 수 있게 하여 연산 비용을 크게 감소시킨다.

둘째, 라그랑주 승수를 이용해 접촉 구속조건을 강제하고, 이를 파라미터 연속법에 삽입한다. 연속 파라미터는 외부 힘의 크기(α)와 진동 주파수(ω)이며, 응답 표면 R(Q,ω,α)=0을 정의한다. 이 표면의 국소 극점, 즉 리지(ridge)는 강제 공진 백본 곡선과 동일하므로, 코코(COCO) 프레임워크를 활용해 리지를 따라 전역적으로 추적한다. 특히, 수치적 특이점(예: 접촉 전이점)에서 발생하는 경사 급변을 해시안 텐서가 제공하는 2차 정보가 정확히 감지해 경로를 안정적으로 통과한다.

셋째, 기존 dNNM과 위상 공명(phase‑resonance) 방법과의 비교를 통해 제안 방법의 적용 범위를 명확히 한다. dNNM은 인공 감쇠 ξ를 도입해 에너지를 주입하지만, 비선형 마찰이 강하거나 구조 감쇠가 큰 경우 수렴이 불안정하고, 위상 공명은 90° 위상 차이를 가정하므로 다중 모드 상호작용을 제대로 포착하지 못한다. 반면, 해시안‑가속 연속법은 비선형 마찰이 지배적인 영역에서도 정확한 백본을 제공하고, 고감쇠 상황에서도 수렴성을 유지한다.

실험적으로는 Euler‑Bernoulli 빔의 유한 요소 모델과 블레이드‑댐퍼‑블레이드 2자유도 모델을 사용했다. 빔 모델에서는 마찰 접촉면이 전이점에서 급격히 변하는데, 기존 방법은 이 지점을 지나면서 해석이 발산하거나 백본이 왜곡되었지만, 제안 방법은 특이점을 정확히 식별하고 연속적으로 추적했다. 블레이드‑댐퍼‑블레이드 모델에서는 구조 감쇠와 마찰 감쇠가 동시에 작용하는 복합 비선형성이 나타났으며, 제안 방법은 실험 데이터와 거의 일치하는 백본을 재현하면서 계산 시간은 기존 방법 대비 40% 이상 단축되었다.

요약하면, 해시안 텐서를 이용한 정확한 2차 미분 정보와 라그랑주 기반 연속 추적이 결합된 이 프레임워크는 마찰 비선형 구조물의 강제 응답 백본 곡선을 고정밀·고효율로 계산할 수 있는 새로운 표준이 될 가능성을 제시한다.