共振破壊と免震ビルの数理（固有値と外力周期）

Q: 共振破壊と免震ビルの数理（固有値と外力周期）とは何ですか？

1自由度の強制減衰振動系における運動方程式は次式で与えられる。 m \ddot{x} + c \dot{x} + k x = F_0 \cos(\omega t) 両辺を質量 m で割り、固有角振動数 \omega_0 = \sqrt{k/m}、減衰比 \zeta = \frac{c}{2\sqrt{mk}} と置くと標準形が得られる。 \ddot{x} + 2\zeta\omega_0 \dot{x} + \omega_0^2 x = \frac{F_0}{m} \cos(\omega t) 定常状態での振幅増幅倍率 M(r)（周波数比 r = \omega / \omega_0）は次のようになる。 M(r) = \frac{1}{\sqrt{(1 - r^2)^2 + (2\zeta r)^2}}

Q: 共振破壊と免震ビルの数理（固有値と外力周期）は社会やプログラミングでどのように使われていますか？

建築耐震設計における固有値モード解析と免震シミュレーション: 多層建物の剛性行列と質量行列から固有振動数を求め、予想される地震波スペクトルとの共振危険度を評価する。実際の構造計算CADやCAE解析では、ルンゲ・クッタ法等の数値積分を用いて時刻歴応答解析を行い、各階の最大層間変形角が許容値以下に収まるようオイルダンパーの最適配置を決定する。

地震や風の周波数と建造物の固有振動数が一致したとき、振幅が無限に跳ね上がる物理現象

この単元の発見をシェア

3行要約（Key Takeaways / AIOサマリー）

1: すべての建造物や橋には固有振動数（揺れやすい周波数）が存在し、外力周波数がこれに一致すると共振（Resonance）を起こして破壊に至る
2: タコマナローズ橋の崩壊や超高層ビルを襲う長周期地震動は、線形代数の固有値解析と2階線形微分方程式でモデル化される
3: 現代の免震構造・オイルダンパーは、減衰比 $\zeta$ を高めて共振倍率のピークを数学的に抑え込む仕組み

なぜこれを学ぶのか？ — 現実社会とテクノロジーでの価値

Purpose & Motivation

なぜ生まれたのか？

橋が風で崩落したり地震で超高層ビルだけが激しく揺れる「共振現象」のメカニズムを解明し、人命を守る免震・制震構造を設計するために誕生した。

身につく武器・思考力

建造物や機械の固有振動数を逆算し、ダンパーによる減衰比 ζ をチューニングして破滅的破壊を数理的に防ぐ「構造力学・振動制御の基礎」が身につく。

実社会・ITでの活躍現場

•超高層ビル（スカイツリーやあべのハルカス）のTMD（同調質量ダンパー）制震設計
•橋梁（吊り橋・斜張橋）の耐風・耐震モード解析
•電気自動車（EV）や航空機エンジンの異音・有害振動シミュレーション
•スピーカーや楽器の共鳴箱（アコースティックチェンバー）の音響設計

Layer 1 — 直感・シミュレーター

地震の周波数 $\omega$ をスライダーで変えてビルの固有振動数 $\omega_0$ （1.00Hz）に近づけてみよう。周波数比 $r \approx 1.0$ に達した瞬間、揺れが急増する「共振」を体感できる。次に「免震ダンパー」スライダーを右に引き上げると、共振周波数であってもビルの揺れが安全に抑え込まれる様子を確認しよう。

外力・地震の周波数 $\omega$0.50 Hz

0.2 Hz (長周期)2.0 Hz (高周波)

ビルの固有振動数 $\omega_0$1.00 Hz

0.5 Hz (超高層)1.5 Hz (中低層)

免震ダンパー減衰比 $\zeta$非免震 (0.05)

0.01 (共振危険)0.25 (高減衰免震)

外力と固有振動数の差により振幅は抑制されています

周波数比 $r = \omega / \omega_0 = 0.50$ ｜理論増幅倍率: 1.3倍

地盤の揺れ（地震動 $F\cos\omega t$）

地面変位: 15.0 px

免震装置 (Damper)

直感チェック！ミニクイズ

免震ビルにおいて、減衰ダンパー（ダンピング比 $\zeta$ ）を増やすと共振時の揺れの振幅はどう変化するか？

Layer 2 — 定義・受験・証明

定義

1自由度の強制減衰振動系における運動方程式は次式で与えられる。

m \ddot{x} + c \dot{x} + k x = F_0 \cos(\omega t)

両辺を質量 $m$ で割り、固有角振動数 $\omega_0 = \sqrt{k/m}$ 、減衰比 $\zeta = \frac{c}{2\sqrt{mk}}$ と置くと標準形が得られる。

\ddot{x} + 2\zeta\omega_0 \dot{x} + \omega_0^2 x = \frac{F_0}{m} \cos(\omega t)

定常状態での振幅増幅倍率 $M(r)$ （周波数比 $r = \omega / \omega_0$ ）は次のようになる。

M(r) = \frac{1}{\sqrt{(1 - r^2)^2 + (2\zeta r)^2}}

公式・数理モデル

M(r) = \frac{1}{\sqrt{(1 - r^2)^2 + (2\zeta r)^2}}, \quad \omega_0 = \sqrt{\frac{k}{m}}

証明・導出

共振点における振幅発散の導出：

運動方程式の特殊解（定常解）を $x_p(t) = X \cos(\omega t - \delta)$ と仮定して微分方程式に代入する。三角関数の加法定理または複素平面上の指数関数表現 $e^{i\omega t}$ を用いて整理すると、振幅 $X$ は次のように求まる。

X = \frac{F_0 / k}{\sqrt{(1 - r^2)^2 + (2\zeta r)^2}}

ここで減衰のない理想系（ $\zeta = 0$ ）を考えると、

X = \frac{F_0 / k}{|1 - r^2|}

となり、外力周波数が系の固有振動数に一致する $r = \omega / \omega_0 \to 1$ の極限で、分母がゼロに近づき振幅 $X \to \infty$ （無限大に発散）する。

多層ビルの場合、質量行列 $M$ と剛性行列 $K$ による多自由度連立方程式 $M \ddot{\mathbf{x}} + K \mathbf{x} = \mathbf{0}$ となり、一般化固有値問題 $(K - \omega^2 M)\mathbf{u} = \mathbf{0}$ の固有値 $\lambda_i = \omega_i^2$ と固有ベクトル $\mathbf{u}_i$ が各階の振動モード形を表す。建築構造設計では、この固有値を事前に計算して地盤の卓越周期と一致しないよう構造剛性を調整している。

押さえておきたいポイント

外力周波数が建造物の固有振動数と一致したとき、エネルギーの蓄積によって振幅が爆発的に増大する現象を「共振（Resonance）」と呼ぶ
1940年のタコマナローズ橋崩壊は、横風によって生じたカルマン渦の周波数が橋のねじれ固有振動数と一致した空気力学的自励振動が原因とされている
高層ビルでは長周期地震動（遠方の大地震でゆっくり長く揺れる波）との共振を防ぐため、TMD（同調質量ダンパー）や粘性オイルダンパーによる $\zeta$ の強化が不可欠である

参照基準・学習指導要領との対応

本トピックの内容・数理モデルは、文部科学省学習指導要領（算数・数学・情報）および国際標準カリキュラムに準拠し、エンジニアリング現場の実装事例と照らし合わせてMath World Lab編集部により検証・作成されています。

Layer 3 — 工学応用・現場実装

建築耐震設計における固有値モード解析と免震シミュレーション

多層建物の剛性行列と質量行列から固有振動数を求め、予想される地震波スペクトルとの共振危険度を評価する。実際の構造計算CADやCAE解析では、ルンゲ・クッタ法等の数値積分を用いて時刻歴応答解析を行い、各階の最大層間変形角が許容値以下に収まるようオイルダンパーの最適配置を決定する。

1// 1自由度建造物の強制減衰振動（ルンゲ・クッタ第4次法）
2function simulateBuildingMotion(
3  naturalFreqHz: number, // 固有振動数 f0 (Hz)
4  dampingRatio: number,  // 減衰比 zeta (0.02: 一般建築, 0.15: 免震ダンパー)
5  forceFreqHz: number,   // 地震周波数 f (Hz)
6  durationSec: number = 5,
7  dt: number = 0.005
8) {
9  const w0 = 2 * Math.PI * naturalFreqHz;
10  const w = 2 * Math.PI * forceFreqHz;
11  const f0 = 10.0; // 地震加振力振幅
12
13  let x = 0; // 変位
14  let v = 0; // 速度
15  let t = 0;
16
17  // 加速度微分方程式: a = -2*zeta*w0*v - w0^2*x + f0*cos(w*t)
18  const deriv = (currX: number, currV: number, currT: number) => {
19    return -2 * dampingRatio * w0 * currV - w0 * w0 * currX + f0 * Math.cos(w * currT);
20  };
21
22  let maxDisplacement = 0;
23
24  while (t < durationSec) {
25    // RK4 4次ルンゲクッタ
26    const k1_v = deriv(x, v, t);
27    const k1_x = v;
28
29    const k2_v = deriv(x + 0.5 * dt * k1_x, v + 0.5 * dt * k1_v, t + 0.5 * dt);
30    const k2_x = v + 0.5 * dt * k1_v;
31
32    const k3_v = deriv(x + 0.5 * dt * k2_x, v + 0.5 * dt * k2_v, t + 0.5 * dt);
33    const k3_x = v + 0.5 * dt * k2_v;
34
35    const k4_v = deriv(x + dt * k3_x, v + dt * k3_v, t + dt);
36    const k4_x = v + dt * k3_v;
37
38    x += (dt / 6) * (k1_x + 2 * k2_x + 2 * k3_x + k4_x);
39    v += (dt / 6) * (k1_v + 2 * k2_v + 2 * k3_v + k4_v);
40    t += dt;
41
42    if (Math.abs(x) > maxDisplacement) {
43      maxDisplacement = Math.abs(x);
44    }
45  }
46
47  return { maxDisplacement, isSafe: maxDisplacement < 25 };
48}
49
50// 検証: 非免震ビルの共振 vs 免震ダンパー追加
51const unshielded = simulateBuildingMotion(1.0, 0.02, 1.0); // 共振・非免震
52const shielded = simulateBuildingMotion(1.0, 0.20, 1.0);   // 共振・免震ダンパー搭載
53
54console.log(`非免震ビル最大変位: ${unshielded.maxDisplacement.toFixed(2)} cm (安全判定: ${unshielded.isSafe})`);
55console.log(`免震ビル最大変位: ${shielded.maxDisplacement.toFixed(2)} cm (安全判定: ${shielded.isSafe})`);