地震応答解析
構造物が地震動を受けたときの動きを力学的に予測・評価する手法を地震応答解析という。設計用入力地震動に対し、変位・速度・加速度、部材力、損傷度などを時間領域または周波数領域で求め、耐震安全性や機能維持性を検討する。解析は一次設計の保守的評価を補完し、免震・制振の効果検証、既存建築の補強方針立案、重要設備の機能喪失リスク管理などに活用される。数理モデルの仮定、入力地震動の選定、減衰や非線形の表現が結果を支配するため、目的に応じたモデル化と検証が不可欠である。
定義と目的
対象構造物を質点系または連続体としてモデル化し、地盤からの強制入力に対する応答を数値的に解くことが中核である。目的は、最大応答や累積損傷の把握、限界状態（層間変形角、塑性率、残留変形）への到達判定、機器・配管・天井など非構造部材の健全性確認にある。性能設計では、地震動レベルごとに許容応答を設定し、解析結果が要求性能を満たすかを判断する。
モデル化の基本
簡略化の基本は1自由度（SDOF）と多自由度（MDOF）である。建築骨組なら各層を質点化し、せん断型・曲げせん断型・連続梁モデルなどを用いる。機械設備や橋梁では、シェル・ソリッド要素を併用する有限要素モデルが一般的である。ばね・ダンパ・質量で表す集中定数系は計算が速く、形状依存の局所応力を扱う場合は分布定数系が有利である。
入力地震動の扱い
観測記録、模擬波（位相特性を考慮した合成）、設計用基準波が利用される。場の不確実性を考慮して複数波を採用し、規模・継続時間・卓越周期のばらつきを含める。応答スペクトル適合は、目標スペクトルに対しフーリエ成分を修正して加速度時刻歴を作る実務的手順である。入力は原則として自由地表面で与え、基礎ばねやインピーダンスで地盤–構造相互作用を表現する。
時間領域と周波数領域
時間領域では加速度時刻歴を直接与え、運動方程式を逐次積分して最大・累積応答を得る。非線形・履歴特性を素直に扱える。一方、周波数領域では線形系を前提に伝達関数やパワースペクトル密度を用い、期待値や分散など統計量を評価する。長時間応答の平均的傾向の把握や感度解析に向く。
数値積分法
代表的手法はNewmark-β法とWilson-θ法である。Newmark-β法はβ、γの選択で無条件安定を実現し、構造解析で広く用いられる。Wilson-θ法は擬似加速度法の一種で、高周波成分の数値減衰を調整しやすい。高精度が必要な高速機器ではRunge-Kutta法が選ばれることもある。時間刻みは固有周期の十分小さい分割（例：最短固有周期の1/20〜1/50）とし、数値安定と応答忠実度の両立を図る。
減衰モデルの設定
比例減衰の代表はレイリー減衰で、低次と高次の2点で目標減衰比に一致させる。粘性ダンパを明示要素として入れる場合は、系全体の等価減衰が過大にならないよう整合性を確認する。材料・接合の非線形に伴うエネルギー散逸は履歴減衰として現れ、粘性仮定のみでは表しきれないため、必要に応じて骨格曲線とルールを与える。
非線形挙動の表現
降伏後の塑性化、付着すべり、開閉口や接触・摩擦、座屈後挙動などをモデル化する。骨組では塑性ヒンジ法またはファイバーモデル、RCではコンクリートの圧壊と鉄筋の降伏・座屈、鋼ではバイリニアやキンメン硬化などを採用する。繰返し載荷での劣化やピンチングも、累積損傷評価に重要である。
性能評価指標
- 層間変形角：機能維持と非構造部材の損傷抑制に直結する。
- 塑性率・累積塑性変形：靭性需要の把握に用いる。
- 最大応答加速度：機器・配管・天井の耐震設計に必要。
- 残留変形：補修・継続使用判断の重要指標。
- エネルギー指標：入力・散逸・ひずみエネルギーの収支で損傷度を見る。
免震・制振デバイスの解析
免震では積層ゴムやすべり支承の水平柔化と復元力を正しく表す。速度依存を持つ粘性ダンパ、履歴ダンパ、TMDの設計では、周波数帯域と目標減衰比に合わせてパラメータ同定を行う。装置の温度・長期劣化・ばらつきは等価モデルの係数に反映し、感度解析でロバスト性を確かめる。
検証とモデル同定
常時微動や加振試験から固有振動数・減衰比・モード形を推定し、数値モデルを調整する。ARX/ARMAやサブスペース法による同定は、観測ノイズ下でも安定した推定を与える。観測記録に対する再現解析は、設計時の不確実性低減と継続的なモデル改善に有用である。
数値安定と収束管理
非線形反復では収束基準、最大反復回数、ステップ制御を適切に定める。剛性の強非線形や接触を含む系では、接線剛性の更新則と接触判定の滑らかさが計算の安定を左右する。メッシュは応力集中や境界近傍を局所的に細分化し、要素歪みやアスペクト比を監視する。
周辺技術とデータ活用
FFTによる周波数成分解析、短時間フーリエ変換やWaveletでの時変特性評価、PSDからの確率論的設計が有効である。大量シミュレーションではパラメトリックスタディやモンテカルロ法で不確実性を評価し、クラウド計算資源で計算負荷を捌く。監視センサの時系列と解析モデルを統合するデジタルツインは、地震直後の健全度推定と早期復旧計画に資する。
関連する入力資料の品質
加速度計の設置方向、ベースラインドリフト、飽和や欠測は応答に直結する。前処理としてトレンド除去、フィルタリング、ベースライン補正を行い、卓越周期の改変や長周期成分の過小評価に注意する。設計用地震動はサイト特性と目標スペクトルの一致を確認する。
よくある誤り
- レイリー減衰の当てはめ周波数が不適切で高次モードが過減衰になる。
- 免震・制振要素の速度・温度依存を無視し、等価粘性を固定してしまう。
- 時刻歴の時間刻みが粗く、ピーク応答と衝撃力を過小評価する。
- 地盤–構造相互作用を無視し、基礎固定の仮定で過剛評価となる。
- 観測記録1波のみで設計判定を行い、ばらつきの影響を見落とす。
以上の留意点を踏まえ、目的整合のモデル化、適切な入力選定、妥当な減衰・非線形表現、検証データによる同定を組み合わせることで、地震応答解析は性能設計・維持管理・リスク評価の実務に有効な意思決定根拠を与える。
コメント（β版）