3차원 축대칭 물체의 지진 응답: 2D 모델의 한계와 새로운 3D 해석

초록

본 연구는 축대칭 강체를 고정점 주변에서 3차원 회전시키는 동역학 모델을 구축하고, 이를 수치적으로 해석함으로써 기존 2차원(2D) 가정의 한계를 검증한다. 3개의 2차 상미분방정식을 유도하고, 반사인파와 실제 엘센트로 지진 기록을 적용해 2D와 3D의 전복 가속도와 응답 차이를 비교하였다. 결과는 단순한 반사·전파 파형에서는 두 모델이 유사하지만, 실제 지진 데이터에서는 수평·수직 성분이 복합적으로 작용해 2D 모델이 과소평가하는 현상을 보였다. 따라서 3D 동역학 해석이 필요함을 강조한다.

상세 분석

이 논문은 기존 연구가 3차원 실제 물체를 축대칭이라는 특수한 대칭성에 의해 2차원 평면 문제로 환원한다는 전제를 비판한다. 저자들은 강체의 회전축이 고정점에 고정된 상황을 가정하고, Euler 각(θ, φ, ψ)으로 물체의 자세를 기술한다. 축대칭 물체의 관성 모멘트는 두 개(I₁=I₂, I₃)만 필요하므로, 라그랑지안에 중력·지진 가속도(수평·수직 성분)와 구속조건을 포함시켜 변분 원리를 적용하였다. 그 결과, 회전 각에 대한 세 개의 2차 상미분방정식이 도출되며, 각 방정식은 상호 비선형 결합항(예: sinθ·φ̇·ψ̇ 등)을 포함한다. 이러한 비선형성은 특히 지진 가속도의 고주파 성분이 강할 때 회전·구르기 모드가 급격히 전이되는 현상을 야기한다.

수치해석에서는 4차 Runge‑Kutta와 적응형 시간 스텝 제어를 이용해 방정식을 통합했으며, 초기 조건은 물체가 완전 평형 상태에 있을 때의 작은 각도 편차로 설정하였다. 먼저 반사인파(반사·전파) 형태의 가속도 펄스를 적용해 최소 전복 가속도(critical overturning acceleration)를 구했는데, 2D 모델과 3D 모델이 거의 일치하였으나 특정 주파수 구간(특히 물체 고유 진동수와 근접한 0.8–1.2 Hz)에서는 3D 모델이 전복을 더 일찍 예측하였다. 이는 회전축이 수평면을 벗어나면서 발생하는 복합 구르기 모드가 2D 모델에 포함되지 않기 때문이다.

실제 지진 데이터(엘센트로 1940년 기록)를 적용했을 때는 상황이 크게 달라졌다. NS(북‑남) 성분과 UD(상‑하) 성분이 동시에 작용하면서 물체는 회전·구르기·전복이 복합적으로 진행되었다. 3D 해석에서는 수직 가속도가 전복에 미치는 영향이 무시할 수 없을 정도로 커졌으며, 2D 모델이 예측한 전복 가속도보다 약 15 % 낮은 값으로 전복이 발생했다. 이는 2D 모델이 수직 성분을 단순히 무시하거나 보조적인 효과로만 고려했기 때문에 발생한 오차이다.

결론적으로, 축대칭 강체의 지진 응답을 정확히 예측하려면 3차원 비선형 동역학을 고려해야 하며, 특히 실제 복합 지진파형에서는 수평·수직 성분의 상호작용이 전복 위험을 크게 변화시킨다. 이러한 결과는 지진 위험 평가와 aseismatic 설계에 중요한 시사점을 제공한다.