クリープ座屈 — 理論と支配方程式
クリープ座屈とは
先生、「クリープ座屈」って普通の座屈と何が違うんですか？
通常の座屈は瞬間的に起きる — 荷重が臨界値を超えた瞬間に座屈変形が始まる。一方クリープ座屈は時間の経過とともにゆっくり進行する。荷重が弾性座屈荷重より低くても、長時間かけるとクリープ変形が蓄積して最終的に座屈に至る。
弾性座屈荷重以下でも座屈する！？ それは怖いですね。
クリープは高温環境で材料が時間とともに変形する現象だ。一定の応力下でもひずみが増加し続ける。このクリープひずみの蓄積が構造の形状を徐々に変え、不安定化させるのがクリープ座屈だ。
クリープの基礎
クリープ現象の基本を教えてください。
一定応力 $\sigma$ 一定温度 $T$ でのクリープひずみは3段階で進行する：
1. 第1期クリープ（遷移クリープ） — ひずみ速度が時間とともに減少
2. 第2期クリープ（定常クリープ） — ひずみ速度が一定。最も長い段階
3. 第3期クリープ（加速クリープ） — ひずみ速度が増大し、最終的に破断
定常クリープのひずみ速度はNorton（べき乗）則で表されることが多い：
ここで $A, n$ は材料定数、$Q$ は活性化エネルギー、$R$ はガス定数、$T$ は絶対温度。
$\sigma^n$ で $n$ が3〜8程度の鋼だと、応力が2倍になるとクリープ速度は8〜256倍！ 応力への感度がものすごく高いですね。
その通り。だからクリープ座屈では応力の再配分が重要になる。初期の弾性応力分布が時間とともにクリープ緩和で均一化されていく。この過程で構造の挙動が変わる。
クリープ座屈のメカニズム
クリープ座屈はどうやって起きるんですか？
2つのメカニズムがある。
1. 分岐型クリープ座屈 — 弾性座屈と同様の分岐が、時間遅れで発生する。圧縮応力下でクリープにより曲げ変形が徐々に増大し、ある時点で急激に座屈する。
2. 擬似座屈（creep buckling by deflection amplification） — 初期不整による曲げ変形がクリープで時間とともに増幅される。明確な分岐点はなく、変形が許容値を超えた時点を「座屈」と定義する。
擬似座屈は「変形が大きくなりすぎる」ことが座屈の定義なんですね。
そう。クリープ座屈の「臨界時間」は、変位が初期値の何倍になったかで定義されることが多い。例えば「変位が初期値の5倍になる時間」を臨界時間とする。
臨界時間の概念
「臨界時間」とは具体的に何ですか？
荷重レベル $P/P_{cr}$（弾性座屈荷重に対する比率）に対応する「座屈までの時間」だ。
Hoffの古典的結果（1958年）では、初期不整を持つ柱のクリープ座屈時間：
荷重が弾性座屈荷重に近いほど $t_{cr}$ は短く、荷重が低いほど $t_{cr}$ は長い。
$P/P_{cr} = 0.5$ でも十分な時間が経てば座屈する可能性があるんですか。
理論的にはそうだ。ただし $P/P_{cr}$ が低い場合、$t_{cr}$ が構造の寿命（数十年）を超えることもある。その場合は実用上クリープ座屈は問題にならない。
クリープ座屈が問題になる分野
どんな構造でクリープ座屈が問題になりますか？
コンクリートも常温でクリープするんですか。
コンクリートは常温でもクリープする（乾燥クリープ）。長期的に大きな持続荷重がかかる柱や壁では、クリープによる附加偏心が座屈耐力を低下させる。設計基準（ユーロコード2等）では長期荷重の影響をクリープ係数で考慮している。
まとめ
クリープ座屈の理論、整理します。
要点：
- クリープ座屈は時間依存の座屈 — 弾性座屈荷重以下でも長時間で座屈し得る
- Norton則 $\dot{\varepsilon}_{cr} = A\sigma^n$ — 応力の $n$ 乗に比例するクリープ速度
- 臨界時間 — 荷重レベルに対応する座屈までの時間
- 高温環境の構造で重要 — ボイラー、原子炉、タービン、化学プラント
- コンクリートでも長期荷重でクリープ座屈が問題になる
時間という次元が加わることで、座屈問題が一気に複雑になるんですね。
そう。弾性座屈は「荷重が臨界値を超えるかどうか」の二択だが、クリープ座屈は「いつ座屈するか」という連続的な問題だ。設計寿命との関係で判断する必要がある。
タイタニック号と安全率の教訓
「不沈」と謳われたタイタニック号は、低温でのリベット材の脆性破壊が沈没の一因とされています。現代の破壊力学CAEでは、温度依存の材料特性と応力拡大係数を計算して「その温度で本当に大丈夫か？」を事前に検証できます。技術の進歩は、過去の悲劇から学んだ結果です。
各項の物理的意味
- 慣性項（質量項）：$\rho \ddot{u}$、つまり「質量×加速度」。急ブレーキで体が前に投げ出された経験はありませんか？ あの「持っていかれる感じ」がまさに慣性力です。重い物体ほど動き出しにくく、動き出したら止まりにくい。地震で建物が揺れるのも、地面が急に動いたのに建物の質量が「置いていかれる」から。静解析ではこの項をゼロにしますが、それは「ゆっくり力をかけるから加速度は無視できる」という仮定です。衝撃荷重や振動問題では絶対に省略できません。
- 剛性項（弾性復元力）：$Ku$ や $\nabla \cdot \sigma$。ばねを引っ張ると「戻ろうとする力」を感じますよね？ あれがフックの法則 $F=kx$ であり、剛性項の本質です。では質問——鉄の棒とゴム紐、同じ力で引っ張るとどちらが伸びるでしょうか？ 当然ゴムです。この「伸びにくさ」がヤング率 $E$ であり、剛性を決めます。よくある勘違い：「剛性が高い＝強い」ではありません。剛性は「変形しにくさ」、強度は「壊れにくさ」で、別の概念です。
- 外力項（荷重項）：体積力 $f_b$（重力など）と表面力 $f_s$（圧力、接触力など）。こう考えてみてください——橋の上のトラックの重さは「中身全体にかかる力」（体積力）、タイヤが路面を押す力は「表面だけにかかる力」（表面力）。風圧、水圧、ボルトの締付力…すべて外力です。ここでありがちな失敗：荷重の方向を間違える。「引張」のつもりが「圧縮」になっていた——笑い話に聞こえますが、3D空間で座標系が回転していると実際に起こります。
- 減衰項：レイリー減衰 $C\dot{u} = (\alpha M + \beta K)\dot{u}$。ギターの弦を弾いてみてください。音は鳴り続けますか？ いいえ、徐々に小さくなりますよね。振動エネルギーが空気抵抗や弦の内部摩擦で熱に変わるからです。車のショックアブソーバーも同じ原理——わざと振動エネルギーを吸収して乗り心地を良くしています。もし減衰がゼロだったら？ 建物は地震の後いつまでも揺れ続けることになります。実際にはそうならないので、適切な減衰の設定が重要です。
仮定条件と適用限界
次元解析と単位系
| 変数 | SI単位 | 注意点・換算メモ |
|---|---|---|
| 変位 $u$ | m（メートル） | mm入力時は荷重・弾性率もMPa/N系に統一すること |
| 応力 $\sigma$ | Pa（パスカル）= N/m² | MPa = 10⁶ Pa。降伏応力との比較時に単位系の不一致に注意 |
| 歪み $\varepsilon$ | 無次元（m/m） | 工学歪みと対数歪みの区別に注意（大変形時） |
| 弾性率 $E$ | Pa | 鋼: 約210 GPa、アルミ: 約70 GPa。温度依存性に注意 |
| 密度 $\rho$ | kg/m³ | mm系ではtonne/mm³（= 10⁻⁹ tonne/mm³ for 鋼） |
| 力 $F$ | N（ニュートン） | mm系ではN、m系ではNで統一 |
数値例：片持ち梁の先端荷重（L=1m, 断面50×100mm, 鋼材E=210GPa, P=1kN）
最大たわみ δ = PL³/(3EI) = 1000×1000³/(3×210000×4,166,667) ≈ 0.381 mm 最大応力 σ = PL×(h/2)/I ≈ 12.0 MPa（降伏応力235MPaに対して安全率19.6）
メッシュ密度を変えた収束性の確認：
ポイント：要素数を4倍にしても結果は0.5%しか変わらない→8,000要素で十分収束。これが「メッシュ収束性」の確認です。
構造解析の収束問題や計算コストに課題を感じていませんか？ — Project NovaSolverは、実務者が日々直面するこうした課題の解決を目指す研究開発プロジェクトです。
クリープ座屈の実務で感じる課題を教えてください
Project NovaSolverは、CAEエンジニアが日々直面する課題——セットアップの煩雑さ、計算コスト、結果の解釈——の解決を目指しています。あなたの実務経験が、より良いツール開発の原動力になります。
実務課題アンケートに回答する →