応力-ひずみ曲線・材料非線形モデルとは
🙋
「材料モデル」って何ですか？線形弾性やバイリニア、Ramberg-Osgoodって、どう使い分けるんですか？
🎓
大まかに言うと、材料の「壊れるまでの強さと柔らかさ」をどう数式で表すか、その方法だよ。例えば自動車の衝突シミュレーションでは、計算が速くて収束しやすいバイリニアモデルが初期設計でよく使われる。上の「材料モデル」ドロップダウンを「バイリニア」に変えてみて。降伏した後、曲線が折れ線になるのがわかるね。
🙋
え、そうなんですか！じゃあ、Ramberg-Osgoodはどんな時に使うんですか？曲線が滑らかですね。
🎓
その通り！降伏点付近のなだらかな変化を精密に再現したい時に使うんだ。例えば、飛行機の翼の継ぎ目など、繰り返し荷重（疲労）を受ける部分の解析だ。ツールで「硬化指数 n」のスライダーを動かしてみて。nが大きいと降伏がよりシャープになる、材料の「硬さ」の感じが変わるのがわかるよ。
🙋
なるほど！「比較材料をオーバーレイ」で他の材料と比べると、コンクリートはすぐに折れてしまうんですね。この違いをCAEソフトではどうやって入力するんですか？
🎓
良いところに気づいたね。実務では、このツールで確認したパラメータをそのままCAEソフトに入力するんだ。例えばANSYSなら「MISO」、LS-DYNAなら「MAT_024」という材料カードに、この曲線を表すデータポイントを入力する。ツールで「構造用鋼」と「アルミ」をオーバーレイ表示して、強さと延びのトレードオフを確認してみよう。それが軽量化設計の第一歩だ。
応力-ひずみ曲線の各領域と機械的性質
引張試験で得られる応力-ひずみ（$\sigma$-$\varepsilon$）曲線は、材料の力学的特性を凝縮した「材料の指紋」です。延性金属では次の領域が現れます。
- 弾性域：フックの法則 $\sigma = E\varepsilon$ が成り立ち、除荷すると元に戻ります。傾きがヤング率 $E$です。
- 降伏点：塑性変形が始まる応力が降伏応力 $\sigma_y$。軟鋼では上降伏点・下降伏点が現れます。
- 塑性域（加工硬化）：変形とともに応力が上昇。最大応力が引張強さ $\sigma_u$（UTS）です。
- ネッキング・破断：UTS以降は局所的にくびれ、破断ひずみ $\varepsilon_f$（伸び）で破壊します。
曲線下の面積は靭性（破壊までに吸収するエネルギー）、弾性域の面積はレジリエンスを表します。脆性材料（鋳鉄・コンクリート）は降伏域がほとんど無く、弾性的に急に破断します。
0.2%耐力と真応力・真ひずみ
0.2%耐力（offset yield）：アルミ合金など明確な降伏点を持たない材料では、永久ひずみ0.2%（$\varepsilon=0.002$）に対応する応力を降伏強度とみなします。弾性線を $\varepsilon=0.002$ だけ平行移動し、曲線との交点を読み取ります。
公称 vs 真応力・真ひずみ：通常の曲線は元の断面積・標点距離で計算した公称（engineering）値です。実際の瞬間断面で評価する真応力・真ひずみは次式で変換します（ネッキング前まで有効）。
$\sigma_{true} = \sigma_{nom}(1+\varepsilon_{nom}), \qquad \varepsilon_{true} = \ln(1+\varepsilon_{nom})$
真応力は公称応力より大きく、真応力-真ひずみ曲線はUTS以降も上昇を続けます。塑性解析やFEMでは真応力-真ひずみを用います。
材料非線形モデルの式
解析では実曲線を数式モデルで近似します。本シミュレーターが対応する代表モデルは次の通りです。
| モデル | 式 | 特徴 |
|---|
| 線形弾性 | $\sigma = E\varepsilon$ | 降伏前のみ。最も単純 |
| バイリニア（等方硬化） | $\sigma=\sigma_y+E_t(\varepsilon-\varepsilon_y)$（降伏後） | 2直線近似。接線剛性 $E_t$ |
| Ramberg-Osgood | $\varepsilon=\dfrac{\sigma}{E}+0.002\left(\dfrac{\sigma}{\sigma_y}\right)^{n}$ | 滑らかな弾塑性遷移。航空材料で多用 |
| べき乗硬化（Hollomon） | $\sigma = K\,\varepsilon_p^{\,n}$ | 塑性域を表現。$n$＝加工硬化指数 |
加工硬化指数 $n$ が大きいほど均一に伸び（成形性が良く）、UTSと降伏応力の比 $\sigma_u/\sigma_y$ が大きいほど延性的です。プリセットを切り替え、鋼・アルミ・チタン・ゴム状ポリマーの曲線形状の違いを比較できます。
よくある質問
Ramberg-Osgoodは降伏点付近の滑らかな遷移を表現できるため、アルミやステンレス鋼など明確な降伏点がない材料に適しています。バイリニアは降伏後の剛性が一定で単純なため、構造用鋼のように明瞭な降伏点を持つ材料や、大変形を伴わない弾塑性解析に向いています。
ツール画面の「材料選択」ドロップダウンで構造用鋼・アルミなどを選ぶと、右側のパラメータ入力欄にデフォルト値が自動反映されます。これらの値は一般的な材料特性に基づいており、必要に応じて手動で変更してシミュレーションできます。
主な原因はパラメータの物理的整合性の欠如です。例えば、降伏応力がヤング率に対して極端に大きい、または硬化指数nが負の値になっていないか確認してください。また、ひずみ範囲が大きすぎると塑性域で計算が発散する場合があります。
参考データとしては有用ですが、そのまま使用する際は注意が必要です。本ツールは単軸・静的負荷を前提とした簡易モデルです。実際のFEM解析では、メッシュサイズや負荷速度、多軸応力状態の影響を考慮し、別途材料試験データでキャリブレーションすることを推奨します。
実世界での応用
自動車車体の衝突安全性解析： 計算コストと収束性に優れるバイリニアモデルが多用されます。エンジンフードやサイドメンバーが衝突時にどのように変形し、エネルギーを吸収するかをシミュレーションし、乗員保護設計に役立てます。
航空機構造部品の疲労寿命評価： 降伏点近傍の応答を精密に再現できるRamberg-Osgoodモデルが適しています。主翼取付部やランディングギアなど、繰り返し荷重を受ける部位のき裂発生寿命を予測する際の材料入力データとして使用されます。
金属板のプレス成形シミュレーション： 大きな塑性変形を伴うため、べき乗硬化則がよく用いられます。自動車ドアパネルなどの複雑形状部品を、割れやシワなく成形できる金型形状と加工条件をCAE上で事前に検討します。
建築・土木構造物の耐震解析： 鉄筋コンクリート柱の応答解析では、鉄筋（バイリニアまたはRamberg-Osgoodモデル）とコンクリート（圧縮側のみ耐力を持つ非対称モデル）を組み合わせた複合材料モデルとして定義され、大地震時の残留変形や崩壊メカニズムの解明に貢献します。
よくある誤解と注意点
まず、「降伏応力は材料の最大強度」と思い込むことです。実は、降伏応力は材料が永久変形し始める目安で、最大強度はその後の加工硬化で上がる「引張強さ（UTS）」です。例えば、構造用鋼は降伏応力が350MPaでも、UTSは500MPaを超えることがあります。CAEで降伏応力だけを入力して「これ以上は壊れる」と判断すると、部材の余力を過小評価して、必要以上に重い設計になる恐れがあります。
次に、モデル選択の安易な一般化。バイリニアモデルは計算が速いからと、すべての非線形解析に使うのは危険です。例えば、ゴムや一部の高分子材料は、降伏点すら明確でない「超弾性」挙動を示します。このツールで「ゴム」を選んでみると、最初からカーブしているのがわかりますよね。こういう材料にバイリニアを当てはめると、変形挙動を大きく見誤ります。材料の本質的な挙動をまず観察することが第一歩です。
最後に、パラメータ「n」の解釈ミス。べき乗硬化則の「n」とRamberg-Osgoodの「n」は名前が同じでも意味が異なります。前者は加工硬化の「しやすさ」、後者は降伏の「なだらかさ」を主に表します。両方のモデルで「n=0.2」に設定して比べてみると、曲線の形状が全く違うことが確認できます。CAEソフトのマニュアルで、その材料モデルにおける各パラメータの物理的定義を必ず確認しましょう。