[MATLAB: 6-3.1절] SHAKE 자유장 해석 (mshake.m)

다음 예제(Kramer, 1996)는 동일한 암반 노두 운동(rock outcrop motion)을 사용할 때에도, 기반암(bed rock)을 각각 강체(rigid)와 경암(hard rock)으로 달리 가정하면, 자유장 해석을 통해 산정한 지표면 운동은 차이가 날 수 있음을 보여준다. 따라서 강체 기반암 가정은 지반조건에 따라 실제 계측된 지진과는 큰 차이를 나타낼 수 있음을 알 수 있다.
가장 널리 사용되는 자유장 해석 프로그램인 SHAKE 시리즈는 앞서 설명한 1차원 파동전파이론을 사용하였다. 이 이론을 MATLAB 코드를 이용하여 구현할 수 있다. 이 예제에서는 참고문헌(김재민, 2009)에서 수록된 MATLAB 코드인 mshake를 사용하였으며, 1989년 캘리포니아 로마 프리에타(Loma Prieta) 지진(Giroy No. 2)을 암반 노두 운동으로 사용하였다(자료: Center for Engineering Strong Motion Data: http://strongmotioncenter.org/).

[첨부자료 참조(구조동역학 3판에 보강될 예정)]