Hybrid Dynamics of Rocking Blocks Beyond Overturning: Saltation Analysis, Bifurcations, and Stability Characterization

기술 요약: 전도 너머의 록킹 블록(Rocking Blocks)에 관한 하이브리드 역학

문제 제기
록킹 블록은 조적 벽체, 조각상, 고가 수조와 같은 강체 구조물의 지진 응답을 이해하는 데 있어 근본적인 모델이다. 전통적으로 전도를 방지하는 것이 주요 공학적 목표였으나, 붕괴 전의 비선형 진동은 다주기 및 카오스적 영역을 포함한 복잡한 거동을 보인다. 이러한 역학을 결정짓는 핵심 요소는 충격 시의 에너지 소산이며, 이는 반발 계수(COR)를 통해 모델링된다.

수십 년 동안 각운동량 보존으로부터 유도된 하우스너(Housner) 모델이 표준으로 사용되어 왔으나, 최근 연구들은 이 모델의 예측과 실험적 충격 거동 사이의 불일치를 지적해 왔다. 마오(Mao) 등이 제안한 운동 에너지 재분배 모델과 같은 대안적 정식화는 특히 세장비가 낮은 블록에 대해 실험 데이터와 더 잘 일치하는 것으로 제안되었다. 그러나 이러한 대안적 반발 정식화를 채택하는 것이 록킹 시스템의 전역적(global) 역학, 즉 어트랙터(attractor)의 안정성, 분기 구조, 그리고 유계 운동(bounded motions)의 접근 가능성에 미치는 영향은 여전히 미개척 영역으로 남아 있다.

방법론
저자들은 조화 기저 진동(harmonic base excitation)을 받는 하이브리드 비매끄러운 동역학계(hybrid non-smooth dynamical system)로서 록킹 블록을 정식화한다. 이 시스템은 질량 $M$, 폭 $2b$, 높이 $2h$를 가진 직사각형 강체로 정의되며, 강체 기초 위에서 미끄러짐 없이 록킹한다. 역학은 활성화된 피벗(pivot)에 따라 전환되는 구간별 자율 미분 방정식(piecewise-autonomous differential equations)에 의해 지배된다.

시스템 분석을 위해 저자들은 두 가지 구별된 반발 모델을 비교한다:

1. 하우스너 모델 ($e_H$): 각운동량 보존으로부터 유도되었으며, 블록의 기하학적 형상(임계각 $\theta_c$)에만 의존한다.
2. 마오 등의 모델 ($e_M$): 속도장의 벡터 기반 분석을 통해 유도되었으며, 접촉 계면에서의 에너지 소산과 탄성 회복을 고려한다.

연구는 다음을 사용하는 포괄적인 수치적 특성화를 수행한다:

• 분기 다이어그램(Bifurcation Diagrams): 가해지는 진폭($\alpha$)이 변화함에 따라 점근적 응답(주기적, 다주기적, 카오스적, 전도적 거동)의 질적 변화를 매핑한다.
• 최대 리아푸노프 지수 ($\lambda_{max}$): 불연속점에서의 충격을 처리하기 위해 살테이션 행렬(saltation matrices)과 QR 재직교화(QR re-orthonormalization)를 사용하는 베네틴-울프(Benettin–Wolf) 알고리즘을 사용하여 계산하며, 이를 통해 카오스적 영역과 안정적 주기적 영역을 식별한다.
• 어트랙터의 유역(Basins of Attraction): 초기 조건 공간을 이산화하여 서로 다른 역학적 결과(감쇠 운동, 지속적 진동, 전도)에 대한 전역적 접근 가능성을 결정한다.

시뮬레이션은 동일한 기하학적 및 진동 파라미터 하에서 두 가지 세장비($h/b = 1$ 및 $h/b = 2$)에 대해 수행되었으며, 이를 통해 반발 법칙이 시스템의 진화에 미치는 영향을 분리하여 조사하였다.

주요 기여 및 결과
본 연구는 반발 모델의 선택이 특히 세장비가 낮은 블록에서 예측된 역학적 지형을 유의미하게 변화시킨다는 것을 입증한다.

• 저세장비 블록($h/b = 1$)에서의 발산:

  • 마오 등의 정식화는 하우스너 모델에 비해 더 높은 반발 값을 예측(더 낮은 에너지 소산)하며, 이는 약 79.85%의 대칭적 백분율 차이를 발생시킨다.
  • 이러한 감소된 소산은 복잡한 진동 거동의 조기 발생과 카오스 역학을 위한 더 넓은 파라미터 범위를 초래한다.
  • 유역 분석: 마오 모델은 안정적인 대칭 진동의 유역을 실질적으로 축소시키며(초기 조건의 약 55%에서 약 20%로 감소), 전도 결과의 빈도를 크게 증가시킨다(약 44%에서 약 80%로 증가). 진동 유역의 경계는 파편화되고 불규칙해지며, 이는 초기 조건에 대한 민감도가 높아졌음을 나타낸다.
  • 두 모델의 역학적 결과에 대한 전반적인 일치도는 $h/b=1$에서 **64.09%**에 불과하며, 이는 초기 조건의 3분의 1 이상이 오직 충격 모델의 사용 여부에 따라 질적으로 다른 운명을 맞이하게 됨을 의미한다.

• 고세장비 블록($h/b = 2$)에서의 수렴:

  • 세장비가 증가함에 따라 두 반발 계수 사이의 차이가 줄어든다(대칭 차이가 약 17%로 감소).
  • 결과적으로, 두 모델의 분기 구조와 리아푸노프 지수 분포는 점진적으로 수렴한다.
  • 마오 모델이 유계 영역 내에 다주기 어트랙터를 도입함에도 불구하고, 모델 간의 전역적 일치도는 **86.93%**로 상승하며, 이는 더 높은 블록이 반발 계수의 특정 정식화에 덜 민감함을 시사한다.

의의 및 주장
본 논문은 반발 계수를 단일 충격에서의 에너지 손실을 지배하는 단순한 국부적 파라미터로 볼 것이 아니라, 위상 공간의 전역적 구조를 근본적으로 변화시킬 수 있는 모델링 가설으로 간주해야 한다고 주장한다.

저자들은 실험적 충격 데이터와 일치시키기 위해 흔히 채택되는 대안적 반발 정식화가 클래식한 하우스너 모델과 비교했을 때 현저히 다른 장기 예측을 생성할 수 있다고 주장한다. 구체적으로, 모델의 선택은 다음 사항에 영향을 미친다:

1. 안정성 임계값 및 역학적 전이(분기)의 위치.
2. 안전한 유계 록킹을 위한 어트랙터의 유역 크기와 견고성.
3. 동일한 진동 조건 하에서의 전도 취약성.

연구는 저세장비 블록의 경우 충격 모델이 시스템의 안전성과 거동을 결정하는 결정적 요인인 반면, 더 높은 블록의 경우 역학적 특징이 반발 모델링의 변화에 대해 더 견고하다는 결론을 내린다. 저자들은 본 분석이 이상적인 강체 가정을 기반으로 하며 재료 변형, 미끄러짐 또는 점탄성 소산을 고려하지 않았음을 언급하며 연구의 범위를 겸허히 제한하였다.