プレストレスモーダル解析

カテゴリ: 構造解析 | 統合版 2026-04-06
CAE visualization for prestressed modal theory - technical simulation diagram
プレストレスモーダル解析

プレストレスモーダルの理論基礎

プレストレスモーダル解析とは

🧑‍🎓

先生、「プレストレスモーダル解析」って何ですか?


🎓

初期応力(プレストレス)が存在する状態での固有振動数解析だ。構造に引張力や圧縮力がかかっている状態では、固有振動数が変化する。弦を張ると音が高くなるのと同じ原理だ。


🧑‍🎓

ギターの弦ですね。張力を上げると振動数が上がる。


🎓

まさにそう。逆に圧縮力がかかると振動数は下がる。圧縮力が臨界座屈荷重に近づくと振動数がゼロに近づく。振動数がゼロ = 座屈だ。


支配方程式

🎓

プレストレスを含む固有値問題:


$$ ([K_0] + [K_\sigma] - \omega^2 [M])\{\phi\} = \{0\} $$

$[K_\sigma]$ は幾何剛性マトリクス(応力剛性)で、座屈解析と同じ。


🧑‍🎓

座屈は $([K_0] + \lambda [K_\sigma])\{\phi\} = \{0\}$ で、振動は $([K_0] + [K_\sigma] - \omega^2 [M])\{\phi\} = \{0\}$。$[K_\sigma]$ が共通!


🎓

完璧な観察だ。座屈と振動は同じ幾何剛性マトリクスを共有する。引張プレストレスは $[K_\sigma] > 0$ で全体剛性を増し、振動数を上げる。圧縮プレストレスは $[K_\sigma] < 0$ で全体剛性を下げ、振動数を下げる。


適用例

🎓
構造プレストレスの種類振動数への影響
弦・ケーブル引張引張で振動数上昇
回転ディスク遠心力(引張)回転速度で振動数上昇
タービンブレード遠心力回転で振動数が変化
圧縮柱軸圧縮圧縮で振動数低下
プレストレスコンクリート梁引張(PC鋼材)引張で振動数やや上昇
膜構造(テント)面内引張引張で振動数上昇
🧑‍🎓

回転体の振動では遠心力によるプレストレスが重要なんですね。


🎓

タービンブレードや回転ディスクでは遠心力で引張プレストレスが発生し、振動数が上がる。これをスピンソフトニング/ハードニングと呼ぶ。回転速度ごとに振動数が変わるから、各速度での固有振動数を評価する必要がある。


FEMでの手順

🎓

1. 静解析(プリロード) — 初期応力(圧縮、引張、遠心力等)を求める

2. 幾何剛性マトリクスの構成 — 静解析の応力から $[K_\sigma]$ を計算

3. 固有値解析 — $[K_0] + [K_\sigma]$ を剛性として固有振動数を求める


🧑‍🎓

座屈解析とほぼ同じ手順ですね。


🎓

Step 1〜2は全く同じ。Step 3で座屈なら $\lambda$ を求め、振動なら $\omega$ を求める。多くのソルバーでは1回のプリロードステップから座屈と振動の両方を実行できる。


まとめ

🧑‍🎓

プレストレスモーダル解析を整理します。


🎓

要点:


  • 初期応力が振動数を変える — 引張で上昇、圧縮で低下
  • $[K_0] + [K_\sigma] - \omega^2 [M] = 0$ — 幾何剛性が加わった固有値問題
  • 座屈と振動は同じ $[K_\sigma]$ — 座屈点で振動数がゼロ
  • 回転体では遠心力プレストレスが重要 — スピンソフトニング/ハードニング
  • 静解析→固有値解析の2段階 — 座屈解析と同じ手順

🧑‍🎓

「座屈点で振動数がゼロ」という関係は深いですね。振動と座屈が一つの理論で繋がっている。


🎓

構造力学の最も美しい関係の一つだ。VCT(Vibration Correlation Technique)はこの関係を利用して、振動数の変化から座屈荷重を非破壊で予測する。


Coffee Break よもやま話

張った弦の高い音と圧縮ばねの振動数

弦楽器の弦に張力をかけると固有振動数が上がる。逆に圧縮荷重をかけると固有振動数は下がり、座屈荷重(Pcr)で0になる。この関係はf²=f₀²(1-P/Pcr)で表され、測定固有振動数から圧縮荷重P(実測が難しいプレストレス量)を逆算できる。この原理は橋梁ケーブルの張力管理に実際に使われている。

プレストレスモーダルの数値計算手法

ソルバー設定

🧑‍🎓

プレストレスモーダル解析の設定を教えてください。


Nastran

```

SOL 103

CEND

SUBCASE 1

LOAD = 100 $ プリロード

METHOD = 10 $ 固有値解析

SPC = 1

```

SOL 103に荷重を設定すると、自動的にプレストレスモーダル解析になる。

Abaqus

```

*STEP

*STATIC

*CLOAD

...

*END STEP

*STEP

*FREQUENCY

20, ,

*END STEP

```

Staticステップの後にFrequencyステップを配置。自動的に前ステップの応力が幾何剛性に反映される。

Ansys

```

/SOLU

ANTYPE, STATIC

PSTRES, ON ! 応力剛性を有効化

SOLVE

FINISH

/SOLU

ANTYPE, MODAL

MODOPT, LANB, 20

SOLVE

```

🧑‍🎓

AnsysではPSTRES, ONが必須でしたよね。座屈のときと同じ。


🎓

PSTRES, ONを忘れると応力剛性が計算されず、プレストレスの影響がゼロの「通常の」固有振動数が出る。座屈と同じ落とし穴。


プリロードの種類

🎓
プリロード設定方法
軸力集中荷重 or 分布荷重引張ケーブル、圧縮柱
遠心力RFORCE(Nastran), *DLOAD CENTRIFUGAL(Abaqus回転体
温度TEMP荷重熱応力による振動数変化
内圧PLOAD4 / *DLOAD P加圧容器の振動
🧑‍🎓

温度によるプレストレスもあるんですか。


🎓

構造が拘束された状態で温度変化があると熱応力が発生し、これが振動数に影響する。火災時の鉄骨梁では温度上昇で軸圧縮力が発生し、振動数が低下→座屈という過程をたどる。


まとめ

🧑‍🎓

プレストレスモーダルの数値手法、整理します。


🎓

要点:


  • 静解析→固有値解析の2段階 — 座屈と同じ手順
  • Nastran SOL 103に荷重を追加するだけ — 最もシンプル
  • Abaqus: Static→Frequency — 自動的にプレストレス反映
  • Ansys: PSTRES, ON を忘れない — 必須設定
  • 遠心力、温度、内圧もプリロードに含まれる

Coffee Break よもやま話

プレストレスモーダル解析の2ステップ手順

プレストレスを考慮したモーダル解析はStep1(静解析でプレストレス状態取得)→Step2(応力剛性行列を加えた固有値解析)の2ステップで行う。ANSYSではPSTRES,ON設定でStep1の解を自動的にStep2に引き継ぐ。AbaqusではInitial Stress Stiffnessオプションが同等機能を持つ。応力剛性の影響が10%以上の場合にプレストレス考慮が必須となる。

プレストレスモーダルの実務適用

プレストレスモーダルの実務適用

🧑‍🎓

プレストレスモーダル解析は実務でどう使われていますか?


キャンベルダイアグラム(回転機械)

🎓

回転体(タービン、コンプレッサー、ファン)では回転速度ごとに固有振動数が変わる。キャンベルダイアグラムは回転速度 vs. 固有振動数のグラフで、回転速度の整数倍の励振線($f = n \times N_{rpm}/60$)との交点が共振条件。


🧑‍🎓

回転速度を変えながら何回も固有値解析を行うんですか?


🎓

そう。10〜20点の回転速度で固有振動数を求め、プロットする。各点で遠心力プレストレスを計算→固有値解析


ケーブル構造の振動

🎓

吊り橋のケーブルやPC梁のPC鋼材は引張プレストレスで振動数が決まる。ケーブルの振動数測定から逆にプレストレス(張力)を推定する張力推定法がある。


$$ T = 4mL^2 f_n^2 / n^2 $$

🧑‍🎓

振動数を測れば張力がわかる! 便利ですね。


🎓

吊り橋のケーブル張力管理に実用化されている。加速度センサーでケーブルの固有振動数を測り、上式で張力を算出する。


実務チェックリスト

🧑‍🎓

プレストレスモーダルのチェックリストをお願いします。


🎓
  • [ ] プリロード(静解析)が正しく収束しているか
  • [ ] 幾何剛性($[K_\sigma]$)が正しく計算されているか(Ansys: PSTRES, ON)
  • [ ] プレストレスなしの固有振動数と比較して、変化が物理的に妥当か
  • [ ] 引張→振動数上昇、圧縮→振動数低下の方向が正しいか
  • [ ] 回転体ではキャンベルダイアグラムを作成したか
  • [ ] 座屈荷重に近いプレストレスで振動数がゼロに近づくか確認したか

  • 🧑‍🎓

    「引張で上昇、圧縮で低下」の方向チェックが重要ですね。逆だったら何かがおかしい。


    🎓

    この方向チェックは最もシンプルで効果的な検証だ。弦の物理を知っていれば直感的に判断できる。


    Coffee Break よもやま話

    PC橋梁の固有振動数変化の実測管理

    プレストレストコンクリート橋の健全度管理に固有振動数の経時変化を使う技術が普及している。PC鋼材の緊張力が設計値の10%低下すると橋の1次固有振動数が約3〜5%低下することがFEM解析で確認されており、定期計測データから緊張力低下を非破壊で検知できる。国土交通省の橋梁点検マニュアル(2019年改訂)でこの手法の参考値として言及されている。

    プレストレスモーダルのソフトウェア比較

    プレストレスモーダルのツール

    🧑‍🎓

    プレストレスモーダルのソルバー比較は?


    🎓
    機能NastranAbaqusAnsys
    プレストレスモーダルSOL 103 + LOADStatic→FrequencyStatic(PSTRES)→Modal
    回転体(遠心力RFORCE*DLOAD, CENTRIFUGALOMEGA
    キャンベルダイアグラムSOL 107/110PythonスクリプトCampbell Diagram Tool
    温度プレストレスTEMP荷重*TEMPERATUREBF, TEMP
    🧑‍🎓

    Ansysにキャンベルダイアグラム専用のツールがあるんですか。


    🎓

    Ansys Workbenchの回転動力学解析モジュールにCampbell Diagram Toolが含まれている。回転速度を自動的にパラメトリックに変化させて、キャンベルダイアグラムを自動生成する。


    選定ガイド

    🎓
    • 回転機械のキャンベルダイアグラム → Ansys Campbell Diagram Tool
    • 航空宇宙のタービンブレードNastran SOL 103
    • 一般のプレストレスモーダル → 手持ちソルバー(どれでも対応)
    • ケーブルの張力推定 → 手計算 + FEM検証

    • 🧑‍🎓

      プレストレスモーダルはどのソルバーでも基本機能ですね。


      🎓

      そう。キャンベルダイアグラムの自動生成はAnsysが便利だが、手動で回転速度を変えて計算すれば他のソルバーでも可能。


      Coffee Break よもやま話

      NASTRAN SOL103 STATSUB機能

      MSC NastranのSOL 103にSTATSUB(SID)オプションを追加すると、先行する静解析(SIDはサブケース番号)の応力場を自動引き継ぎしてモーダル解析を実行できる。Boeing社はこの機能を使い、翼の空力荷重下でのフラッター解析用固有値をプレストレス考慮付きで算出し、787の型式証明申請に使用している。

      プレストレスモーダルの先端研究

      プレストレスモーダルの先端研究

      🧑‍🎓

      プレストレスモーダルの最前線を教えてください。


      VCT(Vibration Correlation Technique)

      🎓

      「振動数がゼロ = 座屈」の関係を利用して、振動数の低下から座屈荷重を非破壊予測する手法。荷重を段階的に増加しながら振動数を測定し、$f^2$ vs. $P$ のプロットから $f = 0$ になる荷重を外挿する。


      🧑‍🎓

      構造を壊さずに座屈荷重がわかる! シェル座屈のページでも出てきましたね。


      🎓

      VCTはプレストレスモーダル解析の最も重要な応用だ。FEMで $f^2$ vs. $P$ の関係をシミュレーションし、実験と比較して検証する。


      非線形プレストレスモーダル

      🎓

      大変形プレストレス(例:膨らんだ膜構造)の固有振動数。線形のプレストレスモーダルでは不十分で、非線形静解析の後に線形化して固有値を求める。AbaqusのNLGEOM=YES → *FREQUENCYで対応。


      デジタルツインとVCT

      🎓

      構造のデジタルツインで振動数をリアルタイム監視し、座屈余裕($P/P_{cr}$)を推定する。橋梁やタンクの構造健全性モニタリングに応用。


      まとめ

      🧑‍🎓

      プレストレスモーダルの先端研究、まとめます。


      🎓
      • VCT — 振動数変化から座屈荷重を非破壊予測
      • 非線形プレストレスモーダル — 大変形状態の振動
      • デジタルツイン — リアルタイム振動数監視で座屈余裕を推定

      • Coffee Break よもやま話

        スピン試験での遠心プレストレス固有値

        高速回転する翼(タービンブレード・プロペラ)では遠心力によるプレストレスが大きく、静止時と回転時で固有振動数が大幅に異なる。典型的なタービンブレードでは回転数を運転条件(3000〜15000rpm)にすると1次固有振動数が20〜50%上昇する。IHI・Kawasaki Aerospaceはスピン試験と連成解析を組み合わせた固有振動数計測で翼設計を検証している。

        プレストレスモーダルのトラブル対応

        プレストレスモーダルのトラブル

        🧑‍🎓

        プレストレスモーダル解析でよくあるトラブルは?


        プレストレスの効果が出ない

        🧑‍🎓

        プリロードをかけたのに振動数が変わりません。


        🎓

        確認項目:

        1. Ansys: PSTRES, ON を設定したか — 忘れると幾何剛性がゼロ

        2. プリロードの静解析が正しく解けているか — 応力がゼロなら $[K_\sigma] = 0$

        3. 荷重が正しい方向か — 引張のつもりが圧縮になっていないか


        振動数が負($\omega^2 < 0$)になる

        🎓

        $[K_0] + [K_\sigma]$ の全体剛性がゼロ以下になると $\omega^2 < 0$。これは座屈を超えた状態(構造が不安定)を意味する。プリロードが座屈荷重を超えている。


        🧑‍🎓

        対策は?


        🎓

        プリロードを減らすか、構造を補強する。$\omega^2 < 0$ は「構造がこの荷重で持たない」というFEMからの警告だ。


        回転体の結果が回転速度に依存しない

        🎓

        遠心力の設定が間違っている可能性。NastranのRFORCEカードで回転速度と回転軸を正しく指定しているか確認。


        まとめ

        🧑‍🎓

        プレストレスモーダルのトラブル対処、整理します。


        🎓
        • 効果が出ない → PSTRES, ON(Ansys)。プリロードの応力確認
        • $\omega^2 < 0$ → 座屈を超えている。プリロードを減らす
        • 回転体で変化なし → 遠心力の設定(回転速度、軸方向)を確認
        • 方向チェック → 引張で振動数上昇、圧縮で低下

        • Coffee Break よもやま話

          プレストレス解析で固有振動数が高い場合

          プレストレスを考慮した解析でも実測より固有振動数が高い場合、接触部のすべりや隙間の影響が考えられる。プレストレスが想定より低い(リラクゼーション・クリープ)場合もある。まず実測のFRFと解析のFRFのピーク形状を比較し、モード形状が一致しているか確認する。形状が一致していれば境界条件・プレストレス量の見直し、形状が異なればモデル化の根本的な再検討が必要だ。

          関連シミュレーター

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

          シミュレーター一覧

          関連する分野

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