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

기술 요약: 전도(Overturning) 너머의 록킹 블록의 하이브리드 역학

문제 제기
록킹 블록(Rocking blocks)은 조적벽, 조각상, 고가 탱크와 같은 강체 구조물의 지진 응답을 이해하기 위한 기초적인 모델이다. 고전적인 하우스너(Housner) 모델이 오랫동안 이러한 시스템을 분석하는 초석 역할을 해왔으나, 최근 연구들은 이 모델의 예측과 실험적으로 관찰된 충격 거동 사이의 불일치를 강조해 왔다. 특히 각운동량 보존에서 유도된 고전적인 반발 계수(COR)는, 특히 세장비(slenderness ratio)가 낮은 블록의 경우 에너지 소산 메커니즘을 정확하게 포착하지 못하는 경우가 많다. 실험 데이터와 더 잘 일치하도록 제안된 Mao 등의 운동 에너지 재분배 모델과 같은 대안적 정식화가 도입되었음에도 불구하고, 이러한 대안적 반발 법칙의 채택이 록킹 시스템의 전역적(global) 비선형 역학에 미치는 영향은 여전히 미개척 영역으로 남아 있다. 본 연구는 반발 모델의 선택이 시스템의 분기 구조(bifurcation structure), 안정성 및 인력권(basin of attraction) 조직에 어떻게 영향을 미치는지 조사하며, 전통적인 전도 임계값 분석을 넘어 복잡한 진동 레짐(regime)을 규명한다.

방법론
본 연구는 조화 기저 진동(harmonic base excitation)을 받는 하이브리드 비매끄러운 동역학계(hybrid non-smooth dynamical system)로서 록킹 블록을 정식화한다. 시스템은 질량 $M$, 폭 $2b$, 높이 $2h$를 가진 직사각형 강체로 정의되며, 미끄러짐 없이 록킹한다. 역학은 각 변위 $\theta$에 따라 전환되는 구간별 미분 방정식에 의해 지배되며, 충격은 $\theta = 0$에서 발생한다.

두 가지 서로 다른 반발 모델이 동일한 기하학적 및 진동 조건 하에서 비교된다:

1. 하우스너 모델 ($e_H$): 반발 계수가 블록의 기하학적 구조(임계각 $\theta_c$)에만 의존하는 고전적 정식화.
2. Mao 등의 모델 ($e_M$): 탄성 회복과 접촉 계면에서의 에너지 소산을 고려하여 속도장의 벡터 기반 분석으로부터 유도된 대안적 정식화.

본 연구는 서로 다른 세장비($h/b = 1$ 및 $h/b = 2$)에 걸친 역학적 특성을 규명하기 위해 세 가지 주요 수치적 도구를 사용한다:

• 분기 다이어그램(Bifurcation Diagrams): 푸앵카레 단면(Poincaré sections)을 사용하여 진폭($\alpha$)이 변화함에 따라 주기적, 다주기적, 카오스적 레짐 사이의 전이를 식별한다.
• 리야푸노프 지수(Lyapunov Exponents): 불연속점에서의 충격을 처리하기 위해 살테이션 행렬(saltation matrices)을 사용하는 Benettin–Wolf 알고리즘을 통해 최대 리야푸노프 지수($\lambda_{max}$)를 계산하여, 안정적인 주기 운동($\lambda_{max} < 0$)과 카오스 역학($\lambda_{max} > 0$)을 구분한다.
• 인력권(Basins of Attraction): 상태 공간을 이산화하여 초기 조건을 점근적 결과(감쇠 운동, 대칭/비대칭 진동, 다주기 응답 또는 전도)에 매핑한다. 이를 통해 유계 운동(bounded motion)의 견고성과 전도 민감도를 정량적으로 평가한다.

주요 결과
비교 결과, 반발 모델의 선택은 예측된 역학적 지형을 유의미하게 변화시키며, 이러한 효과는 세장비에 따라 매우 높은 의존성을 보인다:

• 저세장비 블록에서의 발산 ($h/b = 1$): Mao 등의 정식화는 하우스너 모델에 비해 더 높은 반발 계수(더 낮은 에너지 소산)를 예측한다. 이는 복잡한 진동 거동의 조기 발생과 더 넓은 카오스 역학 관련 파라미터 범위를 초래한다. 결정적으로, 인력권 분석은 극명한 차이를 보여준다: Mao 모델은 안정적인 대칭 진동의 영역을 급격히 감소시키며(55.55%에서 19.64%로), 전도의 빈도를 크게 증가시킨다(44.45%에서 80.36%로). 두 모델 사이의 일치도는 이 기하학적 구조에서 64.09%에 불과하며, 이는 초기 조건의 3분의 1 이상이 충격 법칙에 따라 질적으로 다른 결과를 낳음을 의미한다.
• 고세장비 블록에서의 수렴 ($h/b = 2$): 세장비가 증가함에 따라 두 반발 계수 간의 차이가 줄어든다. 결과적으로, 두 모델에 의해 생성된 역학적 특징들은 점진적으로 수렴한다. Mao 모델이 여전히 유계 영역 내에 다주기 끌개(multiperiodic attractors)를 도입하기는 하지만, 위상 공간의 전반적인 조직과 전도 사건의 비율은 두 정식화 사이에서 더욱 유사해진다.
• 전이의 성격: Mao 모델은 규칙적인 주기 운동의 급격한 상실을 촉진하고 전이를 더 낮은 진동 진폭 쪽으로 이동시킨다. 반면, 하우스너 모델은 전도 전까지 안정성과 카오스의 창(window)이 교차하는 보다 완만한 전이를 보인다.

의의 및 주장
본 논문은 반발 계수가 단순히 국부적인 충격 파라미터가 아니라, 록킹 시스템의 전역적 역학 구조를 근본적으로 변화시킬 수 있는 모델링 가설으로 간주되어야 한다고 주장한다. 본 연구는 대안적인 반발 정식화가 국부적인 충격 응답을 예측하는 데 더 정확할 수 있음에도 불구하고, 고전적인 하우스너 모델과 비교했을 때 현저히 다른 장기적 거동을 생성할 수 있음을 입증한다.

구체적으로, 저세장비 블록의 경우, Mao 등의 정식화를 채택하는 것은 안전한 록킹 응답의 견고성이 감소함을 시사하며, 초기 조건의 작은 섭동이 전도로 이어질 가능성이 더 높음을 보여준다. 반대로, 더 높은 블록의 경우 특정 반발 모델에 대한 민감도가 감소한다. 이러한 결과는 왜 고전적 예측과 실험적 관찰 사이에 불일치가 존재하는지에 대한 역학적 관점을 제공하며, 충격 소산의 처리가 하이브리드 록킹 시스템의 안정성 경계와 끌개의 접근성을 결정하는 핵심 요소임을 강조한다. 본 연구는 Mao 모델이 실험적 충격 데이터에 더 가까운 근사치를 제공하지만, 두 모델 모두에 내재된 이상적인 강체 가정이 실제 구조 시스템에서 발견되는 국부 변형이나 마찰 미끄러짐과 같은 복잡한 접촉 역학을 포착하는 데 한계가 있음을 결론짓는다.