圧力容器の線形解析
圧力容器の線形の理論基礎
圧力容器の力学
先生、圧力容器の応力解析って特別な理論が必要ですか?
圧力容器は回転シェルの力学が基礎になる。薄肉の円筒や球殻に内圧がかかる問題で、膜理論(membrane theory)と曲げ理論(bending theory)の2段階で考える。
薄肉理論(膜理論)
薄肉円筒の応力式を教えてください。
内圧 $p$、内径 $D$、板厚 $t$ の薄肉円筒:
フープ応力は軸方向応力の2倍! だから圧力容器は周方向に裂けるんですね。
まさにそう。ソーセージを茹でると縦に裂けるのと同じ理屈だ。圧力容器の設計ではフープ応力が支配的であり、板厚は主にフープ応力で決まる。
薄肉球殻では:
球殻のほうが応力が低い(円筒のフープ応力の半分)。だから球形タンクは円筒より効率的だが、製造が難しい。
厚肉理論(ラメの問題)
薄肉公式はどの程度の板厚まで使えますか?
$D/t > 20$ 程度が目安。$D/t < 20$ では厚肉理論(ラメの式)が必要:
厚肉では $\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で応力を出すだけでなく、それを正しく分類する能力が圧力容器エンジニアに求められる。
圧力容器設計規格の歴史
圧力容器の安全設計基準は1914年にAMSE(米国機械学会)が「Boiler and Pressure Vessel Code(BPVC)」初版を発行したことで体系化された。その契機は1905年のマサチューセッツ州ブロックトンのBoiler爆発事故(死者58名)だ。現在のASME BPVC Section VIII Division 2(2017年版)は有限要素法に基づく「設計・解析(DBA)」を正式に認めている。
圧力容器の線形の数値計算手法
FEMによる圧力容器解析
圧力容器をFEMで解析する場合、どの要素を使いますか?
モデル化の選択肢は3つある。
| モデル | 要素タイプ | 適用場面 |
|---|---|---|
| 軸対称 | CAX8R等 | 円筒胴、鏡板、軸方向ノズル |
| シェル | S4R, S8R等 | ノズル接続部、全体モデル |
| ソリッド | C3D20R等 | ノズル接続部の詳細、溶接部 |
軸対称が使えるなら最も効率的ですね。
そう。圧力容器の大部分は軸対称だから、まず軸対称で解析し、非軸対称な部分(ノズル、支持脚)だけ3次元で解析するのが標準的なアプローチだ。
応力分類線(Stress Classification Line, SCL)
FEMの結果からASMEの応力分類はどうやりますか?
応力分類線(SCL)を設定して、線上の応力を膜成分と曲げ成分に分離する。SCLは板厚方向に直線を引き、その上の応力分布を積分する。
板厚方向の応力分布 $\sigma(y)$($y$: 板厚位置、$-t/2$ 〜 $t/2$)に対して:
全応力から膜と曲げを引いた残りがピーク応力。応力集中やノッチ効果がここに入る。
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が便利
- 不連続部のメッシュを十分細かく — ナックル、ノズル、溶接部
ラメの応力式とFEM検証
厚肉円筒の内圧応力はラメ(Gabriel Lamé, 1852年)の厳密解が存在し、FEM検証問題として古典的に使われる。内径50mm・外径100mm・内圧10MPaの鋼管では、最大ラジアル方向応力がラメ解で13.3MPaとなる。この問題をAbaqus CAX4要素で解くと、要素数20で誤差0.3%以内に収束し、有名なFEMテキスト(Zienkiewicz 2000年版)にベンチマークとして掲載されている。
圧力容器の線形の実務適用
圧力容器設計の実務フロー
圧力容器の設計は実務でどのように行われますか?
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の出力を自動的に行うことができない。エンジニアの判断が不可欠だ。特に「この応力は一次か二次か」の判断は経験が必要で、圧力容器エンジニアの最も重要なスキルだ。
水素タンクの高圧解析
燃料電池車(FCV)のトヨタMIRAI(2014年発売)に搭載された70MPa高圧水素タンクは、CFRPフィラメントワインディング+アルミライナー構造で内圧解析に軸対称FEMと平面応力シェルの複合モデルを使用した。ASME VIII Div.3基準で定める設計破裂圧力(175MPa)に対するFEM予測は176.2MPaで実験破裂圧の99.3%に一致した。
圧力容器の線形のソフトウェア比較
圧力容器設計のツール
圧力容器の設計にはどんなツールが使われていますか?
規格計算ソフトとFEMの2段構成が一般的だ。
規格計算ソフト(Design by Rule)
| ソフト | 対応規格 | 特徴 |
|---|---|---|
| PVElite | ASME VIII Div.1/2, EN 13445 | 最も広く使われている。コード計算+簡易FEM |
| Compress | ASME VIII Div.1/2 | ASME専門。コードチェック自動化 |
| PASSA/Vessel | EN 13445 | 欧州規格に強い |
| AutoPIPE Vessel | ASME, EN | 配管解析との統合 |
PVEliteが業界標準ですか?
グローバルではそうだ。胴、鏡板、ノズル、フランジ、支持脚の規格計算を一貫して行える。内蔵の簡易FEM(Nozzle Pro相当)でノズル不連続応力も評価可能。
FEM(Design by Analysis)
FEMで圧力容器を解析する場合のソルバー選定は?
Ansysの応力線形化機能が圧力容器には有利ですね。
Ansys WorkbenchのLinearized Stress結果は圧力容器エンジニアに非常に人気がある。SCLを視覚的に設定し、膜・曲げ・ピークの各成分をASME形式で自動出力する。Abaqusでも可能だが、より手動の操作が必要だ。
特殊ツール
選定ガイド
まとめると?
規格計算ソフトで大部分をカバーし、FEMは検証や複雑部位に使うんですね。
そう。圧力容器設計ではFEMは「補助ツール」であり、規格の理解と応力分類のスキルが本質だ。ツールが正しい答えを出しても、エンジニアが応力分類を間違えたら設計は成り立たない。
圧力容器解析のソルバー別認証
原子力圧力容器(原子炉圧力容器など)の解析には、NRC(米国原子力規制委員会)が認証した検証済みFEMコードが要求される。Abaqus(Dassault Systèmes)、ANSYS Mechanical、ADINA(Adina R&D)はすべてNRC NUREG/CR-7103(2012年)認証を取得しており、それぞれASME VIII検証問題で99%以上の精度を文書化している。NASTRANは構造規範には含まれるが原子力特化認証は非取得だ。
圧力容器の線形の先端研究
圧力容器解析の先端トピック
圧力容器の解析技術はどう進化していますか?
3つの方向が活発だ。
弾塑性崩壊解析(DBA-Plastic)
弾性応力分類ではなく、直接塑性解析で設計する手法があると聞きました。
ASME Div. 2 Part 5.2.3 の弾塑性解析法だ。弾完全塑性モデル(降伏応力 = 1.5S)で2倍の設計荷重をかけて、解が収束すれば合格。応力分類が不要になる画期的な手法だ。
応力分類が不要! それは大きなメリットですね。
応力分類の曖昧さ(SCLの位置、一次/二次の判断)を排除できる。ただしFEMの非線形解析が必要で、計算コストは弾性解析より大きい。近年は計算機の性能向上で実用的になりつつある。
Fitness-for-Service(FFS)
既存の圧力容器の健全性評価はどうするんですか?
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 bar | FCV |
| Type IV | 樹脂ライナー+全面FRP巻き | 〜700 bar | FCV(最軽量) |
Type IVはFRP(炭素繊維)でほぼ全ての荷重を担うんですね。
まとめ
圧力容器の先端トピック、まとめます。
圧力容器は「枯れた技術」に見えるが、水素社会や高経年プラントの課題で新たな解析ニーズが生まれている。
応力線形化と疲労評価
圧力容器の疲労寿命は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解析の最終チェックリストをお願いします。
「薄肉公式との整合確認」が出発点。これが合わなければ先に進めないですね。
その通り。単純な理論解と合わせるのはFEM解析の鉄則だ。圧力容器では薄肉公式がその最も基本的な検証手段になる。
ノズル開口部の応力集中ミス
圧力容器のノズル取り合い部(貫通孔)は応力集中係数Kt=2.5〜4に達するが、粗いシェルメッシュでは集中応力を大幅に過小評価する。ASME PTC 60のガイドラインでは、開口部半径rに対しメッシュサイズをr/5以下とすること、かつ要素は二次要素を用いることを推奨している。この基準を守らず3倍粗いメッシュを使った場合、ピーク応力を約40%低く予測した実例がある。
関連トピック
なった
詳しく
報告