圧力容器の線形解析
理論と物理
圧力容器の力学
先生、圧力容器の応力解析って特別な理論が必要ですか?
圧力容器は回転シェルの力学が基礎になる。薄肉の円筒や球殻に内圧がかかる問題で、膜理論(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)」を正式に認めている。
各項の物理的意味
- 慣性項(質量項):$\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$)に対して:
全応力から膜と曲げを引いた残りがピーク応力。応力集中やノッチ効果がここに入る。
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年版)にベンチマークとして掲載されている。
線形要素(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の出力を自動的に行うことができない。エンジニアの判断が不可欠だ。特に「この応力は一次か二次か」の判断は経験が必要で、圧力容器エンジニアの最も重要なスキルだ。
水素タンクの高圧解析
燃料電池車(FCV)のトヨタMIRAI(2014年発売)に搭載された70MPa高圧水素タンクは、CFRPフィラメントワインディング+アルミライナー構造で内圧解析に軸対称FEMと平面応力シェルの複合モデルを使用した。ASME VIII Div.3基準で定める設計破裂圧力(175MPa)に対するFEM予測は176.2MPaで実験破裂圧の99.3%に一致した。
解析フローのたとえ
解析の流れは、実は料理とそっくりです。まず材料を買い出し(CADモデルの準備)、下ごしらえをして(メッシュ生成)、火にかけて(ソルバー実行)、最後に盛り付ける(後処理で可視化)。ここで大事な問いかけ——料理で一番失敗しやすい工程はどこでしょう? 実は「下ごしらえ」なんです。メッシュの品質が悪いと、どんなに優秀なソルバーを使っても結果はめちゃくちゃになります。
初心者が陥りやすい落とし穴
あなたはメッシュ収束性を確認していますか? 「計算が回った=結果が正しい」と思っていませんか? これ、実はCAE初心者が最も陥りやすい罠です。ソルバーは与えられたメッシュで「それなりの答え」を必ず返します。でもメッシュが粗すぎれば、その答えは現実から大きくずれている。最低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で圧力容器を解析する場合のソルバー選定は?
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つの問い
- 「何を解くか」:圧力容器の線形解析に必要な物理モデル・要素タイプが対応しているか。例えば、流体では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)
既存の圧力容器の健全性評価はどうするんですか?
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%低く予測した実例がある。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——圧力容器の線形解析の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告