モード減衰の同定と設定

カテゴリ: 構造解析 | 統合版 2026-04-06
CAE visualization for modal damping theory - technical simulation diagram
モード減衰の同定と設定

モード減衰の同定と設定の理論基礎

減衰とは

🧑‍🎓

先生、構造力学で「減衰」はどう扱いますか?


🎓

減衰は振動のエネルギーが散逸する効果だ。減衰がなければ構造は永遠に振動し続ける。実構造には必ず減衰があり、これが共振時の振幅を有限に保つ。


運動方程式

🎓

減衰付きの運動方程式:


$$ [M]\{\ddot{u}\} + [C]\{\dot{u}\} + [K]\{u\} = \{F(t)\} $$

$[C]$ が減衰マトリクス。$[M]$ と $[K]$ は物理的に明確だが、$[C]$ は一般に不確かさが大きい


🧑‍🎓

減衰マトリクスをどう決めるかが問題なんですね。


🎓

その通り。減衰のモデル化はFEMの動的解析で最も不確かなパラメータだ。


減衰のモデル

🎓

主要な減衰モデル:


1. モード減衰(Modal Damping)

🎓

各モードに減衰比 $\zeta_i$ を直接指定。最も一般的な手法。


$$ \ddot{q}_i + 2\zeta_i \omega_i \dot{q}_i + \omega_i^2 q_i = f_i(t) $$

$q_i$ はモード座標、$\omega_i$ は固有角振動数。


2. レイリー減衰(Rayleigh Damping)

🎓

$[C] = \alpha [M] + \beta [K]$。$\alpha$ と $\beta$ は2つの振動数で $\zeta$ を合わせて決める。


$$ \zeta_i = \frac{\alpha}{2\omega_i} + \frac{\beta \omega_i}{2} $$

3. 構造減衰(Structural Damping)

🎓

ヒステリシス減衰とも呼ばれる。周波数に依存しない減衰で、複素剛性 $[K^*] = K$ で表現。$g$ が構造減衰係数。


🧑‍🎓

3つのモデル、どう使い分けますか?


🎓
モデル適用場面利点欠点
モード減衰モード重畳法各モードに異なる $\zeta$モード解析が前提
レイリー減衰直接積分法(時刻歴解析)時間領域で使える2周波数でしか $\zeta$ を合わせられない
構造減衰周波数応答解析周波数依存なし時間領域では使えない

減衰比の典型値

🎓
構造減衰比 $\zeta$
鋼構造(溶接)0.5〜1%
鋼構造(ボルト接合)1〜2%
RC構造3〜5%
免震構造10〜30%
機械構造1〜3%
複合材構造0.5〜2%
🧑‍🎓

鋼構造で0.5〜1%…非常に小さいですね。


🎓

鋼は内部減衰が小さい。だから鋼構造は共振で大きな振幅が出やすく、減衰の設定が結果に大きく影響する。


まとめ

🧑‍🎓

モード減衰の理論を整理します。


🎓

要点:


  • 減衰は動的解析で最も不確かなパラメータ — 感度分析が必須
  • 3つのモデル — モード減衰、レイリー減衰、構造減衰
  • モード減衰が最も一般的 — 各モードに $\zeta_i$ を指定
  • レイリー減衰は時間領域用 — $\alpha, \beta$ で2周波数を合わせる
  • 減衰比の典型値 — 鋼: 1%、RC: 5%。用途で異なる

🧑‍🎓

減衰の設定で結果が何倍も変わることがあるんですね。共振時の振幅は $1/(2\zeta)$ に比例するから、$\zeta = 1\%$ と $\zeta = 2\%$ で振幅が2倍違う。


🎓

だから減衰は「最も影響が大きく、最も不確かなパラメータ」。減衰に対する感度分析なしに動的解析の結果を信じてはいけない。


Coffee Break よもやま話

減衰比2%という「魔法の数字」

構造減衰比ζ=2%は設計上の慣例として広く使われるが、実際の鋼構造物では0.5〜5%と大幅に異なる。この値はLankford(1954年)が建築物の実測データから統計的に提案した中央値だ。1970年代のUBC(米国建築基準法)に採用されてから「標準値」として定着したが、溶接構造は1%未満、ボルト締結は3〜5%とまったく異なることを忘れてはならない。

モード減衰の同定と設定の数値計算手法

減衰の設定方法

🧑‍🎓

FEMで減衰をどう設定しますか?


Nastranでのモード減衰

```

TABDMP1, 1, CRIT

, 0., 0.02, 100., 0.02, ENDT

```

全モードに $\zeta = 2\%$ を設定。振動数範囲ごとに異なる $\zeta$ も可能。

Abaqusでのモード減衰

```

*MODAL DAMPING

1, 50, 0.02

```

モード1〜50に $\zeta = 2\%$ を一括設定。

レイリー減衰の設定

🎓

$\alpha$ と $\beta$ の決定方法:


2つの振動数 $f_1, f_2$ で $\zeta_1 = \zeta_2 = \zeta$ とする場合:


$$ \alpha = \frac{2\omega_1 \omega_2 \zeta}{\omega_1 + \omega_2}, \quad \beta = \frac{2\zeta}{\omega_1 + \omega_2} $$

🧑‍🎓

$f_1$ と $f_2$ はどう選びますか?


🎓

着目する振動数範囲の下限と上限。例えば地震応答で1〜10 Hzを対象とするなら $f_1 = 1$ Hz, $f_2 = 10$ Hz。この範囲外では減衰比がずれるので注意。


Abaqusでのレイリー減衰

```

*DAMPING, ALPHA=0.5, BETA=0.001

```

減衰の同定

🧑‍🎓

実構造の減衰比はどう測定しますか?


🎓

1. ハンマー加振法 — インパルスハンマーで加振し、加速度を計測

2. 加振器法 — 正弦波/ランダム加振で周波数応答関数(FRF)を取得

3. 減衰の同定 — FRFの半値幅法、または曲線フィッティング


🎓

半値幅法:共振ピークの振幅が $1/\sqrt{2}$ になる2つの振動数 $f_1, f_2$ から:


$$ \zeta = \frac{f_2 - f_1}{2f_n} $$

🧑‍🎓

簡単に測定できるんですね。


🎓

原理はシンプルだが、実験の精度(加振点、計測点、ノイズ処理)が結果に影響する。複数の方法で交差検証するのが安全だ。


まとめ

🧑‍🎓

減衰設定の数値手法、整理します。


🎓

要点:


  • TABDMP1(Nastran), *MODAL DAMPING(Abaqus — モード減衰の設定
  • $\alpha, \beta$ の決定 — 2振動数での合わせ込み
  • 実験モード解析で減衰比を測定 — 半値幅法が基本
  • 減衰比の範囲に注意 — レイリー減衰は指定範囲外でずれる
  • 実験データがなければ文献値を使い、感度分析する

Coffee Break よもやま話

半値幅法と対数減衰率の実測手順

実験から減衰比を求める半値幅法では、FRF(周波数応答関数)の共振ピーク周波数から±3dB点の幅Δfとピークfrからζ=Δf/(2fr)で計算する。精度は周波数分解能に依存し、Δfの5分の1以下の分解能が必要だ。インパルスハンマ試験では自由減衰波形の包絡線の対数から対数減衰率を直接求める方が精度が高い場合もある。

モード減衰の同定と設定の実務適用

減衰の実務適用

🧑‍🎓

減衰の設定は実務でどう行いますか?


設計コードの減衰値

🎓

設計コードで減衰比が規定されている場合がある:


コード構造減衰比 $\zeta$
建築基準法(日本)RC造5%
建築基準法(日本)S造2%
ユーロコード8RC造5%
ASCE 7全構造5%(デフォルト)
API 617回転機械実測または1%
🧑‍🎓

地震応答ではRC造5%がデフォルトなんですね。


🎓

RC造の5%は「ひび割れを含むRC構造の等価減衰」を意味する。無損傷のRC構造の減衰は1〜2%だが、地震で微小なひび割れが入ることでエネルギー散逸が増加し、等価的に5%になる。


実験データがない場合

🎓
  • 文献値を参考に保守的な値を選択 — 低い減衰比は保守的(振幅が大きくなる)
  • 上下限の2ケースで感度分析 — $\zeta_{low}$ と $\zeta_{high}$ の2ケース
  • 接合部の影響を考慮 — ボルト接合は溶接より減衰が大きい

  • 実務チェックリスト

    🧑‍🎓

    減衰設定のチェックリストをお願いします。


    🎓
    • [ ] 減衰比の値が設計コードまたは実験データに基づいているか
    • [ ] 減衰モデル(モード/レイリー/構造)が解析手法と整合しているか
    • [ ] レイリー減衰の場合、着目振動数範囲で $\zeta$ が適切か(範囲外のずれ確認)
    • [ ] 感度分析($\zeta$ を上下に変えた2ケース)を実施したか
    • [ ] 共振時の振幅が $1/(2\zeta)$ のオーダーと整合するか

    • 🧑‍🎓

      減衰比を $\zeta/2$ にすると共振振幅が2倍…感度分析は必須ですね。


      🎓

      減衰はFEMの動的解析で最大の不確定要素。感度分析を怠ると、結果に意味がない。


      Coffee Break よもやま話

      免震建物のダンパー最適化

      超高層ビルの制振設計では粘性ダンパーや摩擦ダンパーをモーダル減衰比として換算し、目標減衰比(通常3〜5%)を達成するダンパー配置を最適化する。六本木ヒルズ森タワー(238m)では粘性制振ダンパー設置後に低次モードの減衰比が1%から4%に上昇し、地震応答加速度を40%低減させた。

      モード減衰の同定と設定のソフトウェア比較

      減衰設定のツール

      🧑‍🎓

      減衰設定に関するソルバー間の差は?


      🎓
      機能NastranAbaqusAnsys
      モード減衰TABDMP1*MODAL DAMPINGMDAMP
      レイリー減衰PARAM,ALPHA1/2*DAMPING, ALPHA/BETAALPHAD, BETAD
      構造減衰PARAM,G*DAMPING, STRUCTURALDMPSTR
      材料減衰GE on MAT1*DAMPING, STRUCTURAL (mat)MP,DAMP
      粘弾性減衰CVISC*VISCOELASTICPRONY系列
      🧑‍🎓

      材料ごとに異なる減衰を設定できますか?


      🎓

      NastranのMAT1カードのGEフィールドで材料ごとの構造減衰係数を設定できる。鋼とゴムで異なる減衰を持つ構造で有用。AbaqusやAnsysでも同様の機能がある。


      実験モード解析ツール

      ツール特徴
      LMS Test.LabFRFからの減衰同定。FEMとの相関
      ME'scopeモード形状と減衰の可視化
      ModalVIEWExcelベース。手軽なモード解析

      選定ガイド

      🎓
      • モード重畳法 → モード減衰(全ソルバー対応)
      • 直接時刻歴解析 → レイリー減衰($\alpha, \beta$)
      • 周波数応答解析構造減衰($g$)
      • 粘弾性材料 → Prony系列(Abaqusが最も柔軟)
      • 実験データの活用 → LMS Test.Lab + FEMモデル更新

      • 🧑‍🎓

        解析手法によって使える減衰モデルが異なるんですね。


        🎓

        モード減衰はモード重畳法、レイリー減衰は直接積分法、構造減衰は周波数応答。これを間違えると物理的に意味のない結果が出る。


        Coffee Break よもやま話

        MSC Nastranモーダル減衰入力形式

        Nastran SOL 108(周波数応答)ではTABDMP1カードでモーダル減衰比の周波数依存性を定義できる。各モードに個別の減衰比を設定するCRITDAMP=MODEオプションも使え、実測値をそのまま入力可能だ。LMS(現Siemens Simcenter)社はこの機能を活用したテスト-解析相関フローを構築し、試験から1時間以内にFEMの減衰モデルを更新するシステムを自動車向けに提供している。

        モード減衰の同定と設定の先端研究

        減衰の先端研究

        🧑‍🎓

        減衰モデルの最前線を教えてください。


        非粘性減衰(Non-viscous Damping)

        🎓

        レイリー減衰は粘性減衰($[C]\{\dot{u}\}$、速度に比例)を仮定するが、実構造の減衰は非粘性であることが多い。摩擦減衰(接合部の微小すべり)、ヒステリシス減衰(材料の非線形履歴)は粘性モデルでは正確に表現できない。


        🎓

        一般化減衰モデル

        $$ [M]\{\ddot{u}\} + \int_0^t [G(t-\tau)]\{\dot{u}(\tau)\}d\tau + [K]\{u\} = \{F(t)\} $$

        $[G(t)]$ は減衰核関数。粘性減衰は $[G(t)] = [C]\delta(t)$ の特殊ケース。


        接合部の減衰

        🎓

        ボルト接合部の微小すべり(フレッティング)が構造全体の減衰の大部分を占めることがある。接合部減衰モデル(Iwan要素、Jenkins要素)が研究されている。


        🧑‍🎓

        接合部の微小すべりが減衰の主因?


        🎓

        溶接の鋼構造では材料自体の内部減衰が非常に小さい($\zeta < 0.1\%$)。実測の1〜2%の減衰は接合部のフレッティングと非構造材(塗装、ケーブル等)の摩擦が主因。


        能動減衰(Active Damping)

        🎓

        圧電素子やMRダンパーで能動的に減衰を制御する。FEMではセンサー-アクチュエータのフィードバック制御と構造解析を連成する。


        まとめ

        🧑‍🎓

        減衰の先端研究、まとめます。


        🎓
        • 非粘性減衰 — 減衰核関数による一般化。摩擦ヒステリシス
        • 接合部減衰 — Iwan/Jenkins要素。フレッティングの定量化
        • 能動減衰圧電素子・MRダンパーのフィードバック制御

        • 減衰は構造力学の「最もソフトな話題」だが、振動制御の実効性を決める最も重要なパラメータだ。


          Coffee Break よもやま話

          材料損失係数とモーダル減衰の関係

          構造物全体のモーダル減衰比ζnは各部材の材料損失係数ηと部材がモードで蓄えるひずみエネルギーの比率で決まる(モーダルひずみエネルギー法)。高減衰ゴム(η=0.3〜1.0)を少量挟むだけで構造全体の減衰比を大幅に高められるのはこの原理による。制振鋼板(η=0.1〜0.3)を自動車フロアに使うと室内騒音を3〜6dB低減できる。

          モード減衰の同定と設定のトラブル対応

          減衰設定のトラブル

          🧑‍🎓

          減衰設定でよくあるトラブルは?


          共振振幅が非現実的

          🧑‍🎓

          共振での応答が非常に大きくなります。


          🎓

          減衰が小さすぎる。減衰比を確認。$\zeta = 0$(減衰なし)だと共振で振幅が無限大。


          確認:

          • 減衰が設定されているか(設定忘れが最多)
          • 減衰比の値が適切か(文献値・コード値と比較)
          • $\zeta = 0.01$ で振幅が $1/(2 \times 0.01) = 50$ 倍。これが妥当か

          レイリー減衰で低周波/高周波の応答がおかしい

          🎓

          レイリー減衰は $\alpha$ 項が低周波で支配的(減衰大)、$\beta$ 項が高周波で支配的(減衰大)。指定した2周波数の間では $\zeta$ が適切だが、範囲外では減衰が過大になる。


          🧑‍🎓

          どう対処しますか?


          🎓
          • 2周波数を着目範囲の上下限に設定
          • 着目範囲外のモードを結果から除外
          • Extended Rayleigh(3点以上で合わせる拡張版)を検討

          • 構造減衰を時間領域で使ってしまう

            🎓

            構造減衰($g$: ヒステリシス)は周波数領域専用。時間領域(直接積分法)で構造減衰を使うと物理的に意味のない結果(因果律の破れ)が出る。


            🧑‍🎓

            時間領域ではモード減衰かレイリー減衰を使うべきですね。


            🎓

            そう。周波数領域 = 構造減衰OK。時間領域 = レイリー or モード減衰。これを間違えると因果律が破れて非物理的な応答が出る。


            まとめ

            🧑‍🎓

            減衰設定のトラブル対処、整理します。


            🎓
            • 共振振幅過大 → 減衰が設定されているか確認。$\zeta = 0$ は無限大振幅
            • レイリー減衰の範囲外 → 着目範囲の上下限で $\alpha, \beta$ を合わせる
            • 構造減衰を時間領域で使わない → 周波数領域専用
            • 減衰の感度分析は必須 → $\zeta$ を $\times 2$ と $\times 0.5$ で2ケース
            • 減衰設定忘れが最多のミス — 動的解析では必ず確認

            • Coffee Break よもやま話

              モーダル減衰比が実測より低く解析される場合

              FEM解析でのモーダル減衰比は実測より低く見積もられることが多い。これはFEMが接合部(ボルト・リベット・溶接)での摩擦・微小すべりによる散逸エネルギーを表現できないためだ。接合部の減衰を集中減衰要素(C要素)として付加するか、実測減衰比をモーダル解析にそのまま入力するのが実務的な対処法だ。

              関連シミュレーター

              この分野のインタラクティブシミュレーターで理論を体感しよう

              シミュレーター一覧

              関連する分野

              熱解析製造プロセス解析V&V・品質保証
              この記事の評価
              ご回答ありがとうございます!
              参考に
              なった
              もっと
              詳しく
              誤りを
              報告
              参考になった
              0
              もっと詳しく
              0
              誤りを報告
              0
              Written by NovaSolver Contributors
              Anonymous Engineers & AI — サイトマップ
              プロフィールを見る