圧力容器の線形解析

カテゴリ: 構造解析 | 統合版 2026-04-06
CAE visualization for pressure vessel theory - technical simulation diagram
圧力容器の線形解析

理論と物理

圧力容器の力学

🧑‍🎓

先生、圧力容器の応力解析って特別な理論が必要ですか?


🎓

圧力容器は回転シェルの力学が基礎になる。薄肉の円筒や球殻に内圧がかかる問題で、膜理論(membrane theory)と曲げ理論(bending theory)の2段階で考える。


薄肉理論(膜理論)

🧑‍🎓

薄肉円筒の応力式を教えてください。


🎓

内圧 $p$、内径 $D$、板厚 $t$ の薄肉円筒:


$$ \sigma_\theta = \frac{pD}{2t} \quad \text{(フープ応力)} $$
$$ \sigma_z = \frac{pD}{4t} \quad \text{(軸方向応力)} $$

🧑‍🎓

フープ応力は軸方向応力の2倍! だから圧力容器は周方向に裂けるんですね。


🎓

まさにそう。ソーセージを茹でると縦に裂けるのと同じ理屈だ。圧力容器の設計ではフープ応力が支配的であり、板厚は主にフープ応力で決まる。


🎓

薄肉球殻では:


$$ \sigma = \frac{pD}{4t} \quad \text{(全方向一様)} $$

球殻のほうが応力が低い(円筒のフープ応力の半分)。だから球形タンクは円筒より効率的だが、製造が難しい。


厚肉理論(ラメの問題)

🧑‍🎓

薄肉公式はどの程度の板厚まで使えますか?


🎓

$D/t > 20$ 程度が目安。$D/t < 20$ では厚肉理論(ラメの式)が必要:


$$ \sigma_r = \frac{p_i a^2}{b^2 - a^2} \left(1 - \frac{b^2}{r^2}\right) $$
$$ \sigma_\theta = \frac{p_i a^2}{b^2 - a^2} \left(1 + \frac{b^2}{r^2}\right) $$

🧑‍🎓

厚肉では $\sigma_r$ も無視できないんですね。内面で $\sigma_r = -p$(圧縮)。


🎓

そう。薄肉では $\sigma_r \approx 0$ と仮定するが、厚肉では板厚方向の応力勾配が重要になる。内面のフープ応力は薄肉公式より高くなる。


不連続応力

🧑‍🎓

圧力容器の「不連続応力」とは何ですか?


🎓

これが圧力容器解析の核心だ。一様な円筒や球殻だけなら膜理論で十分だが、実際の圧力容器には:


  • 胴と鏡板の接続部
  • ノズル接続部
  • 板厚の変化部
  • 支持部

これらの幾何学的不連続で、膜理論だけでは満たせない適合条件を満足するために局所的な曲げ応力が発生する。これが不連続応力だ。


🧑‍🎓

膜応力の上に、局所的な曲げ応力が重畳するんですね。


🎓

そう。不連続応力は不連続部から離れると急速に減衰する。減衰距離は概ね $\sqrt{Rt}$ のオーダー($R$: 半径、$t$: 板厚)。この「影響範囲」の中では膜理論だけでは不正確で、FEMによる曲げ解析が必要だ。


ASME規格の応力分類

🧑‍🎓

ASME規格では応力をどう分類するんですか?


🎓

ASME BPVC Section VIII Div. 2 Part 5 の応力分類:


応力カテゴリ記号意味許容値
一般膜応力$P_m$容器全体に作用する膜応力$S$(許容応力)
局所膜応力$P_L$不連続部の局所膜応力$1.5S$
一次曲げ応力$P_b$荷重平衡に必要な曲げ応力$1.5S$($P_L + P_b$で)
二次応力$Q$自己制限性のある変位起因応力$3S$($P_L + P_b + Q$で)
ピーク応力$F$局所的な応力集中疲労評価に使用
🧑‍🎓

膜応力が一番厳しくて、二次応力は3倍まで許容されるんですか。


🎓

二次応力は「自己制限性」(self-limiting)がある。局所的な降伏が起きても、変形が拘束されるため塑性崩壊には至らない。だから膜応力より大きな許容値が設定されている。不連続応力の大部分は二次応力に分類される。


まとめ

🧑‍🎓

圧力容器の理論を整理します。


🎓

要点:


  • 薄肉公式 — $\sigma_\theta = pD/(2t)$。フープ応力が軸方向の2倍
  • 厚肉はラメの式 — $D/t < 20$ で必要。板厚方向の応力勾配
  • 不連続応力 — 幾何学的不連続部での局所曲げ。FEMが必須
  • ASME応力分類 — $P_m, P_L, P_b, Q, F$ の5カテゴリ
  • 二次応力は3Sまで許容 — 自己制限性のある応力

🧑‍🎓

圧力容器の設計は「どこが膜応力でどこが二次応力か」を正しく分類することが鍵なんですね。


🎓

まさにそう。応力分類を間違えると、安全側にも危険側にもなる。FEMで応力を出すだけでなく、それを正しく分類する能力が圧力容器エンジニアに求められる。


Coffee Break よもやま話

圧力容器設計規格の歴史

圧力容器の安全設計基準は1914年にAMSE(米国機械学会)が「Boiler and Pressure Vessel Code(BPVC)」初版を発行したことで体系化された。その契機は1905年のマサチューセッツ州ブロックトンのBoiler爆発事故(死者58名)だ。現在のASME BPVC Section VIII Division 2(2017年版)は有限要素法に基づく「設計・解析(DBA)」を正式に認めている。

各項の物理的意味
  • 慣性項(質量項):$\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で統一

数値解法と実装

FEMによる圧力容器解析

🧑‍🎓

圧力容器をFEMで解析する場合、どの要素を使いますか?


🎓

モデル化の選択肢は3つある。


モデル要素タイプ適用場面
軸対称CAX8R等円筒胴、鏡板、軸方向ノズル
シェルS4R, S8R等ノズル接続部、全体モデル
ソリッドC3D20R等ノズル接続部の詳細、溶接部
🧑‍🎓

軸対称が使えるなら最も効率的ですね。


🎓

そう。圧力容器の大部分は軸対称だから、まず軸対称で解析し、非軸対称な部分(ノズル、支持脚)だけ3次元で解析するのが標準的なアプローチだ。


応力分類線(Stress Classification Line, SCL)

🧑‍🎓

FEMの結果からASMEの応力分類はどうやりますか?


🎓

応力分類線(SCL)を設定して、線上の応力を膜成分と曲げ成分に分離する。SCLは板厚方向に直線を引き、その上の応力分布を積分する。


🎓

板厚方向の応力分布 $\sigma(y)$($y$: 板厚位置、$-t/2$ 〜 $t/2$)に対して:


$$ \sigma_m = \frac{1}{t} \int_{-t/2}^{t/2} \sigma(y) \, dy \quad \text{(膜応力)} $$
$$ \sigma_b = \frac{6}{t^2} \int_{-t/2}^{t/2} \sigma(y) \cdot y \, dy \quad \text{(曲げ応力)} $$
$$ \sigma_{peak} = \sigma_{max} - \sigma_m - \sigma_b \quad \text{(ピーク応力)} $$

🧑‍🎓

全応力から膜と曲げを引いた残りがピーク応力。応力集中やノッチ効果がここに入る。


🎓

SCLの設定位置は解析者が決める。不連続部の近傍に板厚方向に垂直な線を引く。ASMEのガイドライン(Div. 2, Part 5, Annex 5-A)でSCLの設定方法が規定されている。


メッシュ要件

🧑‍🎓

圧力容器の解析ではメッシュにどんな注意が必要ですか?


🎓

板厚方向のメッシュが最も重要。応力分類のためには板厚方向の応力勾配を正確に捕捉する必要がある。


  • 軸対称モデル — 板厚方向に最低4要素(二次要素なら2要素)
  • ソリッドモデル — 同上
  • シェルモデル — 板厚方向は自動的に考慮される(積分点で応力出力)

🎓

不連続部のメッシュ密度

  • ナックル部(鏡板と胴の接続): 板厚の1/2以下の要素サイズ
  • ノズル接続部: ノズル板厚の1/2以下
  • 溶接部: 溶接脚長の1/3以下

Abaqus

🎓

```

*STEP

*STATIC

*DLOAD

cylinder, P, 10.0 $ 内圧 10 MPa

*END STEP

```


内圧は *DLOAD の P(圧力)オプションで面に直接与える。方向は要素の法線方向。


Nastran

🎓

```

PLOAD4, 100, elem_id, 10.0, , , , , 1

```


PLOAD4カードで面圧を定義。最後のフラグで法線方向を指定。


Ansys

🎓

```

SF, area_num, PRES, 10.0

```


またはWorkbenchのPressure荷重で面を選択して圧力値を入力。


応力線形化ツール

🧑‍🎓

応力分類は手動でやるんですか?


🎓

ポストプロセッサに応力線形化機能があるものが便利だ:


  • Abaqus/CAE — Path機能でSCLを定義し、膜/曲げを自動分離
  • Ansys Workbench — Linearized Stress機能でSCL上の応力を直接出力
  • Nastran — 直接の機能はないが、f06の応力から手動またはスクリプトで計算

🧑‍🎓

Ansys Workbenchの機能が一番使いやすそうですね。


🎓

Ansys WorkbenchのLinearized Stressは圧力容器エンジニアに人気がある。SCLを視覚的に設定して、膜・曲げ・ピークの各成分と等価応力(von Mises/Tresca)を自動出力してくれる。


まとめ

🧑‍🎓

圧力容器の数値解法、整理します。


🎓

要点:


  • 軸対称+3Dサブモデルが効率的 — 全体を3次元にする必要はない
  • 板厚方向に最低4要素 — 応力分類のための勾配捕捉
  • SCLで膜・曲げ・ピークに分離 — ASME Div. 2の応力分類に必須
  • ポストプロセッサの応力線形化機能を活用 — Ansys Workbenchが便利
  • 不連続部のメッシュを十分細かく — ナックル、ノズル、溶接部

Coffee Break よもやま話

ラメの応力式とFEM検証

厚肉円筒の内圧応力はラメ(Gabriel Lamé, 1852年)の厳密解が存在し、FEM検証問題として古典的に使われる。内径50mm・外径100mm・内圧10MPaの鋼管では、最大ラジアル方向応力がラメ解で13.3MPaとなる。この問題をAbaqus CAX4要素で解くと、要素数20で誤差0.3%以内に収束し、有名なFEMテキスト(Zienkiewicz 2000年版)にベンチマークとして掲載されている。

線形要素(1次要素)

節点間を線形補間。計算コストは低いが、応力の精度が低い。せん断ロッキングに注意(低減積分やB-bar法で緩和)。

2次要素(中間節点付き)

曲線的な変形を表現可能。応力精度が大幅に向上するが、自由度は約2〜3倍に増加。推奨:応力評価が重要な場合。

完全積分 vs 低減積分

完全積分:過剰拘束(ロッキング)のリスク。低減積分:アワーグラスモード(零エネルギーモード)のリスク。適材適所で選択。

アダプティブメッシュ

誤差指標(ZZ推定量等)に基づく自動細分化。応力集中部の精度を効率的に向上。h法(要素分割)とp法(次数増加)がある。

ニュートン・ラフソン法

非線形解析の標準的手法。接線剛性マトリクスを毎反復更新。収束半径内で2次収束するが、計算コストが高い。

修正ニュートン・ラフソン法

接線剛性マトリクスを初期値または数反復毎に更新。各反復のコストは低いが、収束速度は線形的。

収束判定基準

力の残差ノルム: $||R|| / ||F_{ext}|| < \epsilon$(一般に $\epsilon = 10^{-3}$〜$10^{-6}$)。変位増分ノルム: $||\Delta u|| / ||u|| < \epsilon$。エネルギーノルム: $\Delta u \cdot R < \epsilon$

荷重増分法

全荷重を一度に負荷せず、小刻みに増加させる。弧長法(Riks法)は荷重-変位関係の極値点を越えて追跡可能。

直接法 vs 反復法のたとえ

直接法は「連立方程式を筆算で正確に解く」方法——確実だが大規模問題では時間がかかりすぎる。反復法は「当て推量を繰り返して正解に近づく」方法——最初は大雑把な答えだが、反復するたびに精度が上がる。辞書で言葉を探すとき、最初のページから順番に探す(直接法)より、見当をつけて開き、前後に調整する(反復法)方が効率的なのと同じ原理。

メッシュの次数と精度の関係

1次要素は「定規で曲線を近似する」——直線の折れ線で表現するため精度に限界がある。2次要素は「フレキシブルカーブ」——曲線的な変化を表現でき、同じメッシュ密度でも格段に精度が向上する。ただし、1要素あたりの計算コストは増えるため、トータルのコスト対効果で判断する。

実践ガイド

圧力容器設計の実務フロー

🧑‍🎓

圧力容器の設計は実務でどのように行われますか?


🎓

ASME BPVC Section VIII では2つの設計方法がある:


Division 1(Design by Rule)

🎓

公式・チャートに基づく設計。FEMは不要。主要な部品(胴、鏡板、ノズル、フランジ)それぞれに設計式がある。シンプルな形状ならこれで十分。


Division 2(Design by Analysis)

🎓

FEMによる応力解析に基づく設計。Part 4(Design by Rule)とPart 5(Design by Analysis)の2段階。Part 5では:


  • 弾性応力解析法 — FEMの線形弾性解析結果を応力分類して許容値と比較
  • 弾塑性解析法 — FEMの非線形弾塑性解析で塑性崩壊荷重を直接評価

🧑‍🎓

Division 1とDivision 2の使い分けは?


🎓

Division 1はシンプルで安全側だが、板厚が厚くなりがち。Division 2はFEMで精密に評価するため板厚を薄くできるが、解析コストがかかる。高圧・大型の容器や軽量化が必要な場合にDivision 2を使う。


ノズル接続部の解析

🧑‍🎓

ノズル接続部は圧力容器の最も重要な部位ですよね。


🎓

ノズル接続部は応力集中が最も大きい部位であり、破壊・疲労の起点になりやすい。


🎓

解析のポイント:

  • WRC 297/537 — ノズル接続部の応力を簡易的に計算する方法。局所応力係数を提供
  • FEMでの詳細解析 — ソリッド要素で溶接形状まで含めてモデル化
  • 応力分類 — ノズル接続部ではSCLの設定位置が特に重要。ASME Annex 5-Aに従う

🧑‍🎓

WRC 297/537とFEMの結果は一致しますか?


🎓

概ね一致するが、10〜20%の差が出ることがある。WRC は無限円筒殻に取り付けたノズルの理論解に基づいているため、端部効果やノズル同士の干渉がある場合は不正確だ。重要な容器ではFEMで検証すべき。


鏡板の設計

🧑‍🎓

鏡板(ヘッド)にはどんな種類がありますか?


🎓
鏡板タイプ形状応力特性コスト
半球形完全な球殻最小の応力最も高い
半楕円形(2:1)楕円回転体ナックル部に不連続応力中程度
皿形(トリスフェリカル)球殻+円錐接線ナックル部に高応力比較的安い
平板平面最大の応力(曲げ支配)最も安い
🧑‍🎓

半楕円形の「ナックル部に不連続応力」って何ですか?


🎓

ナックル(knuckle)は鏡板と胴の接続部で、曲率が急変する部位だ。膜理論ではこの曲率変化を扱えないため、局所的な曲げ応力(不連続応力)が発生する。FEMの軸対称モデルで正確に評価できる。


🎓

2:1半楕円鏡板のナックル部応力は、胴のフープ応力の1.5〜2.0倍になることがある。これがDivision 2でFEMが求められる典型的なケースだ。


実務チェックリスト

🧑‍🎓

圧力容器解析のチェックリストをお願いします。


🎓
  • [ ] 設計圧力・設計温度・腐れ代が正しく設定されているか
  • [ ] 膜応力が薄肉公式と整合しているか(FEMの検証)
  • [ ] 不連続部(ナックル、ノズル)のメッシュは十分か
  • [ ] SCLの設定位置はASME Annex 5-Aに従っているか
  • [ ] 応力分類($P_m, P_L, P_b, Q, F$)は正しいか
  • [ ] 各カテゴリの応力が許容値以内か
  • [ ] 疲労評価が必要な場合、サイクル数と応力範囲を確認したか
  • [ ] Division 1のルールチェックとFEMの結果が矛盾していないか

  • 🧑‍🎓

    応力分類が正しいかどうかが最も判断が難しそうですね。


    🎓

    その通り。応力分類はFEMの出力を自動的に行うことができない。エンジニアの判断が不可欠だ。特に「この応力は一次か二次か」の判断は経験が必要で、圧力容器エンジニアの最も重要なスキルだ。


    Coffee Break よもやま話

    水素タンクの高圧解析

    燃料電池車(FCV)のトヨタMIRAI(2014年発売)に搭載された70MPa高圧水素タンクは、CFRPフィラメントワインディング+アルミライナー構造で内圧解析に軸対称FEMと平面応力シェルの複合モデルを使用した。ASME VIII Div.3基準で定める設計破裂圧力(175MPa)に対するFEM予測は176.2MPaで実験破裂圧の99.3%に一致した。

    解析フローのたとえ

    解析の流れは、実は料理とそっくりです。まず材料を買い出し(CADモデルの準備)、下ごしらえをして(メッシュ生成)、火にかけて(ソルバー実行)、最後に盛り付ける(後処理で可視化)。ここで大事な問いかけ——料理で一番失敗しやすい工程はどこでしょう? 実は「下ごしらえ」なんです。メッシュの品質が悪いと、どんなに優秀なソルバーを使っても結果はめちゃくちゃになります。

    初心者が陥りやすい落とし穴

    あなたはメッシュ収束性を確認していますか? 「計算が回った=結果が正しい」と思っていませんか? これ、実はCAE初心者が最も陥りやすい罠です。ソルバーは与えられたメッシュで「それなりの答え」を必ず返します。でもメッシュが粗すぎれば、その答えは現実から大きくずれている。最低3段階のメッシュ密度で結果が安定することを確認する——これを怠ると「コンピュータが出した答えだから正しいはず」という危険な思い込みに陥ります。

    境界条件の考え方

    境界条件の設定は、試験の「問題文を書く」のと同じです。問題文が間違っていたら? どんなに正確に計算しても答えは間違いますよね。「この面は本当に完全固定なのか」「この荷重は本当に一様分布なのか」——現実の拘束条件を正しくモデル化することが、実は解析全体で最も重要なステップだったりします。

    ソフトウェア比較

    圧力容器設計のツール

    🧑‍🎓

    圧力容器の設計にはどんなツールが使われていますか?


    🎓

    規格計算ソフトとFEMの2段構成が一般的だ。


    規格計算ソフト(Design by Rule)

    ソフト対応規格特徴
    PVEliteASME VIII Div.1/2, EN 13445最も広く使われている。コード計算+簡易FEM
    CompressASME VIII Div.1/2ASME専門。コードチェック自動化
    PASSA/VesselEN 13445欧州規格に強い
    AutoPIPE VesselASME, EN配管解析との統合
    🧑‍🎓

    PVEliteが業界標準ですか?


    🎓

    グローバルではそうだ。胴、鏡板、ノズル、フランジ、支持脚の規格計算を一貫して行える。内蔵の簡易FEM(Nozzle Pro相当)でノズル不連続応力も評価可能。


    FEM(Design by Analysis)

    🧑‍🎓

    FEMで圧力容器を解析する場合のソルバー選定は?


    🎓
    観点AbaqusAnsysNastran
    軸対称解析CAX8R(標準的)PLANE183CQUADX8
    応力線形化Path機能(手動)Linearized Stress(自動)スクリプト必要
    ASME Div.2 Part 5弾性・弾塑性とも対応弾性・弾塑性とも対応弾性のみ実用的
    非線形(弾塑性崩壊)*STATIC, RIKSArc-LengthSOL 106/400
    🧑‍🎓

    Ansysの応力線形化機能が圧力容器には有利ですね。


    🎓

    Ansys WorkbenchのLinearized Stress結果は圧力容器エンジニアに非常に人気がある。SCLを視覚的に設定し、膜・曲げ・ピークの各成分をASME形式で自動出力する。Abaqusでも可能だが、より手動の操作が必要だ。


    特殊ツール

    🎓
    • Nozzle Pro(Paulin Research Group) — ノズル接続部のFEM解析に特化。WRC 297/537との比較機能内蔵
    • FE/Pipe(Paulin Research Group) — 配管の局所応力評価
    • midas FEA NX — 圧力容器・プラント構造に特化したFEM

    • 選定ガイド

      🧑‍🎓

      まとめると?


      🎓
      • ASME Div.1の規格計算 → PVElite / Compress
      • ASME Div.2の弾性応力分類Ansys Workbench(Linearized Stress)
      • ASME Div.2の弾塑性解析Abaqus(Riks法の実績)
      • ノズル接続部の詳細 → Nozzle Pro(専用ツール)
      • 全体の効率 → PVElite(規格計算)+ Ansys/Abaqus(FEM検証)の2段構成

      • 🧑‍🎓

        規格計算ソフトで大部分をカバーし、FEMは検証や複雑部位に使うんですね。


        🎓

        そう。圧力容器設計ではFEMは「補助ツール」であり、規格の理解と応力分類のスキルが本質だ。ツールが正しい答えを出しても、エンジニアが応力分類を間違えたら設計は成り立たない。


        Coffee Break よもやま話

        圧力容器解析のソルバー別認証

        原子力圧力容器(原子炉圧力容器など)の解析には、NRC(米国原子力規制委員会)が認証した検証済みFEMコードが要求される。Abaqus(Dassault Systèmes)、ANSYS Mechanical、ADINA(Adina R&D)はすべてNRC NUREG/CR-7103(2012年)認証を取得しており、それぞれASME VIII検証問題で99%以上の精度を文書化している。NASTRANは構造規範には含まれるが原子力特化認証は非取得だ。

        選定で最も重要な3つの問い

        • 「何を解くか」:圧力容器の線形解析に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
        • 「誰が使うか」:初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車(GUI)とMT車(スクリプト)の違いに似ている。
        • 「どこまで拡張するか」:将来の解析規模拡大(HPC対応)、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。

        先端技術

        圧力容器解析の先端トピック

        🧑‍🎓

        圧力容器の解析技術はどう進化していますか?


        🎓

        3つの方向が活発だ。


        弾塑性崩壊解析(DBA-Plastic)

        🧑‍🎓

        弾性応力分類ではなく、直接塑性解析で設計する手法があると聞きました。


        🎓

        ASME Div. 2 Part 5.2.3 の弾塑性解析法だ。弾完全塑性モデル(降伏応力 = 1.5S)で2倍の設計荷重をかけて、解が収束すれば合格。応力分類が不要になる画期的な手法だ。


        🧑‍🎓

        応力分類が不要! それは大きなメリットですね。


        🎓

        応力分類の曖昧さ(SCLの位置、一次/二次の判断)を排除できる。ただしFEMの非線形解析が必要で、計算コストは弾性解析より大きい。近年は計算機の性能向上で実用的になりつつある。


        Fitness-for-Service(FFS)

        🧑‍🎓

        既存の圧力容器の健全性評価はどうするんですか?


        🎓

        API 579-1/ASME FFS-1(Fitness-for-Service)規格が適用される。腐食、クリープ、亀裂などの劣化がある既存容器を、そのまま使い続けられるかどうかを評価する。


        🎓

        FEMの適用:

        • Level 3評価 — 詳細なFEM解析による健全性評価
        • 亀裂のある容器 — J積分やSIFをFEMで計算し、破壊力学的に評価
        • 腐食減肉 — 実測の板厚分布をFEMモデルに反映して応力評価

        🧑‍🎓

        腐食で板厚が不均一になった容器をFEMで直接解析するんですね。


        🎓

        実測データ(UT測定やレーザースキャン)をFEMメッシュにマッピングする。板厚が場所によって異なるモデルで応力解析し、最薄部の応力が許容値以内かを確認する。


        水素環境での設計

        🎓

        水素貯蔵容器(Type I〜IV)は圧力容器の最もホットな分野だ。


        タイプ構造最大圧力用途
        Type I全金属〜200 bar工業用
        Type II金属ライナー+部分FRP巻き〜300 barバス
        Type III金属ライナー+全面FRP巻き〜700 barFCV
        Type IV樹脂ライナー+全面FRP巻き〜700 barFCV(最軽量)
        🧑‍🎓

        Type IVはFRP(炭素繊維)でほぼ全ての荷重を担うんですね。


        🎓

        Type IVの設計はFRP積層の解析が核心だ。内圧によるフープ応力を繊維が負担し、軸方向はドーム部の形状で荷重を伝達する。軸対称FEMの複合材積層解析(CLT + FEM)が設計ツールとして使われている。


        まとめ

        🧑‍🎓

        圧力容器の先端トピック、まとめます。


        🎓
        • 弾塑性崩壊解析 — 応力分類不要の直接的手法。実用化が進む
        • Fitness-for-Service — 既存容器の継続使用可否をFEMで評価
        • 水素貯蔵容器 — Type IV CFRPタンクの積層最適化

        • 圧力容器は「枯れた技術」に見えるが、水素社会や高経年プラントの課題で新たな解析ニーズが生まれている。


          Coffee Break よもやま話

          応力線形化と疲労評価

          圧力容器の疲労寿命はASME VIIIの応力線形化(Stress Linearization)手順で評価する。断面上で膜応力・曲げ応力・ピーク応力に分解し、それぞれに許容値を設ける。1990年代のABB原子力のPWR圧力容器ノズル解析では、Abaqusで応力線形化パスを200本自動処理するスクリプトを開発し、手動評価比で解析工数を90%削減した事例が核工学ジャーナルに掲載された。

          トラブルシューティング

          圧力容器解析のトラブル

          🧑‍🎓

          圧力容器のFEM解析でよくあるトラブルを教えてください。


          🎓

          圧力容器特有のトラブルは応力分類に関連するものが多い。


          フープ応力が薄肉公式と合わない

          🧑‍🎓

          FEMのフープ応力が $pD/(2t)$ と一致しません。


          🎓

          確認項目:


          1. 端部条件 — 開放端か密閉端かで軸方向応力が変わり、ポアソン効果でフープ応力にも影響

          2. 不連続部の影響範囲内 — 鏡板接続部から $\sqrt{Rt}$ 以内では不連続応力が重畳。十分離れた位置で比較

          3. 厚肉効果 — $D/t < 20$ ではラメの式との比較が必要。内面応力は薄肉公式より高い

          4. 圧力の作用面 — 内面圧力の向きが法線方向(外向き)か確認


          🧑‍🎓

          端部から離れた位置で比較するのがポイントですね。


          🎓

          胴の中央部(端部から十分離れた一様な領域)で薄肉公式と比較する。ここで一致しなければモデル化に問題がある。


          応力分類で迷う

          🧑‍🎓

          「この応力は一次か二次か」がわかりません。


          🎓

          最も難しい判断だ。基本原則:


          🎓

          一次応力($P$) — 荷重平衡に必要な応力。除去すると構造が平衡を保てない。

          • 例: 内圧による膜応力、自重による応力

          二次応力($Q$) — 変位の適合条件から生じる応力。自己制限性がある。

          • 例: 不連続応力(胴と鏡板の接続部)、熱応力

          🎓

          判断のテスト: 「この応力を取り除いても構造は崩壊しないか?」→ Yes なら二次。No なら一次。


          🧑‍🎓

          不連続応力は取り除いても構造は崩壊しない(形状が変わって応力が自動的に再配分される)から二次。なるほど。


          SCLの設定位置が結果に影響する

          🧑‍🎓

          SCLの位置を少し変えると応力分類の結果が変わります。


          🎓

          これはSCL法の本質的な弱点だ。対策:


          • ASME Annex 5-Aのガイドラインに従う — SCLの設定位置が規定されている
          • SCLは板厚方向に直交させる — 曲面上では局所的な法線方向
          • 複数のSCL位置で結果を確認 — 最も厳しい位置を採用
          • 弾塑性解析法の検討 — SCLが不要な手法で検証

          🧑‍🎓

          弾塑性解析法はSCLの問題を根本的に解決するんですね。


          🎓

          そう。応力分類の曖昧さに悩むなら、弾塑性解析法(ASME Div. 2 Part 5.2.3)への移行を検討すべきだ。計算コストは上がるが、判断の曖昧さが排除される。


          圧力の方向ミス

          🧑‍🎓

          内圧を与えたのに応力が圧縮になりました。


          🎓

          圧力の方向が逆になっている。FEMでは圧力は面の法線方向に作用する。法線が内向きの場合、正の圧力は外圧として作用する。


          🎓

          確認:

          • 要素の法線方向を表示して確認
          • 内圧なら変位が外向き(半径増大)になるはず
          • 反力が圧力×面積と整合するか確認

          実務チェックリスト追加

          🧑‍🎓

          圧力容器FEM解析の最終チェックリストをお願いします。


          🎓
          • [ ] 胴の中央部でフープ応力が薄肉公式と整合するか(基本検証)
          • [ ] 圧力の方向(法線方向)は正しいか(変位方向で確認)
          • [ ] 不連続部のメッシュは板厚の1/2以下の要素サイズか
          • [ ] SCLの位置はASME Annex 5-Aに従っているか
          • [ ] 応力分類(一次/二次/ピーク)は物理的に妥当か
          • [ ] 全応力カテゴリが許容値以内か
          • [ ] 疲労評価が必要な場合、ピーク応力を使っているか
          • [ ] 弾塑性解析法での検証も検討したか

          • 🧑‍🎓

            「薄肉公式との整合確認」が出発点。これが合わなければ先に進めないですね。


            🎓

            その通り。単純な理論解と合わせるのはFEM解析の鉄則だ。圧力容器では薄肉公式がその最も基本的な検証手段になる。


            Coffee Break よもやま話

            ノズル開口部の応力集中ミス

            圧力容器のノズル取り合い部(貫通孔)は応力集中係数Kt=2.5〜4に達するが、粗いシェルメッシュでは集中応力を大幅に過小評価する。ASME PTC 60のガイドラインでは、開口部半径rに対しメッシュサイズをr/5以下とすること、かつ要素は二次要素を用いることを推奨している。この基準を守らず3倍粗いメッシュを使った場合、ピーク応力を約40%低く予測した実例がある。

            「解析が合わない」と思ったら

            1. まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
            2. 最小再現ケースを作る——圧力容器の線形解析の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
            3. 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
            4. 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
            関連シミュレーター

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

            シミュレーター一覧

            関連する分野

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