軸対称解析
理論と物理
軸対称問題とは
先生、「軸対称」解析って何ですか？ 2次元で3次元の問題を解けるんですか？
そう。形状、材料、荷重、境界条件が全て回転軸に対して対称な場合、3次元問題を2次元断面の解析に帰着できる。円筒座標 $(r, \theta, z)$ で $\theta$ 方向に変化がないという仮定だ。
どんな構造が軸対称ですか？
実はものすごく多い：
- 圧力容器 — 円筒胴、鏡板、ノズル（周方向に一様な場合）
- フランジ・ボルト — 締結体の軸力解析
- ベアリング — 内輪・外輪・転動体の接触
- ピストン・シリンダー — 内燃機関の圧力解析
- パイプの内圧 — 厚肉円筒のラメの問題
- 回転体 — ディスク、フライホイール、タービンロータの遠心力
圧力容器の設計では必須ですね。
圧力容器はASME規格でも軸対称解析が認められている。3次元で圧力容器全体をモデル化するより、軸対称の断面解析のほうが100分の1以下のDOFで同等以上の精度が出る。
支配方程式
軸対称問題の応力成分は？
円筒座標系で4つの応力成分が存在する：
平面応力や平面ひずみは3成分でしたが、軸対称は4成分ですか？
そう。軸対称の特徴はフープ応力 $\sigma_\theta$ が存在すること。$r$ 方向に変位 $u_r$ があると、周方向のひずみは：
$u_r / r$！ 変位を半径で割るだけでひずみが出るんですね。これは幾何学的な効果ですか？
完璧な理解だ。半径 $r$ の円周は $2\pi r$ だから、半径が $u_r$ だけ増えると円周は $2\pi(r + u_r)$ になる。ひずみは $\Delta L / L = u_r / r$。この$1/r$ の項が軸対称解析の核心で、平面応力/ひずみにはない特徴だ。
構成則
軸対称のフックの法則：
$$ \begin{Bmatrix} \sigma_r \\ \sigma_\theta \\ \sigma_z \\ \tau_{rz} \end{Bmatrix} = \frac{E}{(1+\nu)(1-2\nu)} \begin{bmatrix} 1-\nu & \nu & \nu & 0 \\ \nu & 1-\nu & \nu & 0 \\ \nu & \nu & 1-\nu & 0 \\ 0 & 0 & 0 & \frac{1-2\nu}{2} \end{bmatrix} \begin{Bmatrix} \varepsilon_r \\ \varepsilon_\theta \\ \varepsilon_z \\ \gamma_{rz} \end{Bmatrix} $$
3次元のフックの法則と同じ形ですね。4成分だけど実質3次元。
そう。軸対称は「見た目は2次元、中身は3次元」だ。応力状態は完全に3次元で、$\sigma_\theta$ がゼロになることはない（$r = 0$ を除いて）。
厚肉円筒（ラメの問題）
内圧を受ける厚肉円筒の理論解はありますか？
ラメの問題（Lamé problem）は軸対称の最も基本的な理論解だ。内半径 $a$, 外半径 $b$, 内圧 $p$ のとき：
$$ \sigma_r = \frac{p a^2}{b^2 - a^2} \left(1 - \frac{b^2}{r^2}\right) $$
$$ \sigma_\theta = \frac{p a^2}{b^2 - a^2} \left(1 + \frac{b^2}{r^2}\right) $$
内面（$r = a$）でフープ応力が最大で、外面に向かって減少する。内面は引張、外面に近づくにつれ引張が小さくなる。
内面の最大フープ応力：
$$ \sigma_{\theta,max} = p \frac{a^2 + b^2}{b^2 - a^2} $$
軸対称のフックの法則：
3次元のフックの法則と同じ形ですね。4成分だけど実質3次元。
そう。軸対称は「見た目は2次元、中身は3次元」だ。応力状態は完全に3次元で、$\sigma_\theta$ がゼロになることはない（$r = 0$ を除いて）。
内圧を受ける厚肉円筒の理論解はありますか？
ラメの問題（Lamé problem）は軸対称の最も基本的な理論解だ。内半径 $a$, 外半径 $b$, 内圧 $p$ のとき：
内面（$r = a$）でフープ応力が最大で、外面に向かって減少する。内面は引張、外面に近づくにつれ引張が小さくなる。
内面の最大フープ応力：
薄肉（$b \approx a$）の場合は $\sigma_\theta \approx pD/(2t)$（薄肉円筒の公式）に近づく。これはFEMの軸対称解析を検証する最良のベンチマーク問題だ。
まとめ
軸対称解析の理論を整理します。
要点：
- 形状・荷重・材料が回転対称なら2次元で解ける — DOFが100分の1以下に
- 4つの応力成分 — $\sigma_r, \sigma_\theta, \sigma_z, \tau_{rz}$。フープ応力 $\sigma_\theta$ が固有
- $\varepsilon_\theta = u_r / r$ — 幾何学的な周方向ひずみ
- 実質3次元の応力状態 — 見た目は2次元だが中身は3次元
- ラメの問題で検証 — 厚肉円筒の内圧は最良のベンチマーク
圧力容器の設計では軸対称を使わない手はないですね。3次元の100分の1のコストで同等精度。
その通り。ただし非軸対称な荷重（ノズル荷重、地震荷重）や非軸対称な形状（穴、不連続部）がある場合は3次元解析が必要だ。軸対称でスクリーニングし、必要な部分だけ3次元で詳細解析するのが効率的なアプローチだ。
軸対称の歴史的起源
軸対称解析の理論的基盤は1850年代にG.B.エアリーが応力関数を導入したことに遡る。1960年代にClough & Rashidがリング要素を定式化し、フロップ数を3次元解析の約1/10に削減。圧力容器や原子炉格納容器の設計で今も現役の手法だ。
各項の物理的意味
- 慣性項（質量項）：$\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要素はどんな特徴がありますか？
軸対称要素は $(r, z)$ 平面の2次元要素だが、内部の定式化に$2\pi r$ の重みが入る。これは要素を回転軸回りに一周させた体積に対応する。
要素の定式化
軸対称要素の剛性マトリクスの積分：
平面問題の $t$（板厚）の代わりに $2\pi r$（周長）が入る。
$r = 0$（回転軸上）で $2\pi r = 0$ になりますが、問題ないですか？
非常にいい質問だ。$r = 0$ での積分は特異性を持つ。実装上は以下の対応がされている：
- ガウス積分点は $r = 0$ 上には配置されない（積分点は要素内部にある）
- $r = 0$ 上の節点では $\varepsilon_\theta = u_r / r$ が $0/0$ になるが、L'Hôpitalの法則で $\varepsilon_\theta = \varepsilon_r$ となる
- 実用上は問題ないが、$r = 0$ 付近で非常に細かいメッシュを使うと数値的に不安定になることがある
ソルバー別の要素名
| 要素 | Nastran | Abaqus | Ansys |
|---|---|---|---|
| 4節点四角形 | CQUADX4 / CTRIAX | CAX4, CAX4R, CAX4H | PLANE182 (KEYOPT3=1) |
| 8節点四角形 | CQUADX8 | CAX8, CAX8R | PLANE183 (KEYOPT3=1) |
| 3節点三角形 | CTRIAX6 | CAX3 | PLANE182（退化） |
AbaqusのCAXは「Continuum AXisymmetric」の略ですね。CPSが平面応力、CPEが平面ひずみ、CAXが軸対称。
その通り。Ansysでは同じPLANE182/183要素のKEYOPTで切り替える。NastranではCQUADX系の別要素になる。
境界条件の設定
軸対称解析で注意すべき境界条件は？
1. 回転軸（$r = 0$）の拘束
- $r = 0$ 上の節点は $u_r = 0$（半径方向変位ゼロ）を強制
- $u_z$ は自由（軸方向には動ける）
- これを忘れると軸が横に動いてしまい、非物理的な結果になる
2. 対称面の扱い
- 軸方向の対称面（例：円筒の中央断面）では $u_z = 0$
- これにより半モデルにできる
3. 荷重の入力
- 集中力は「$2\pi r$ で割った値」ではなく、断面上の線荷重として入力
- 内圧は面圧として直接入力（$2\pi r$ は自動で考慮される）
荷重入力で間違えやすいのは集中力ですか？
そう。軸対称解析で「100 kNの集中荷重」を $r = 50$ mm の節点に入力すると、実際には周方向に分布した $100 / (2\pi \times 0.05) = 318$ kN/m の線荷重になる。入力値がそのまま「リング状荷重」になることを理解していないと、荷重のオーダーが合わない。
非軸対称荷重の扱い（フーリエ展開）
荷重が軸対称でない場合はどうするんですか？
フーリエ展開（Fourier series expansion）で非軸対称成分を分離できる。荷重を周方向の高調波（$\cos n\theta, \sin n\theta$）に分解し、各高調波に対して2次元解析を行い、最後に重ね合わせる。
NastranではCQUADX + SOL 101 でフーリエ展開解析ができる（ハーモニック要素）。Abaqusでは CAXA 系要素（asymmetric-axisymmetric）で非軸対称荷重を扱える。
ノズル荷重のような局所的な非軸対称荷重もフーリエ展開で扱えますか？
原理的には可能だが、局所荷重をフーリエ級数で表現するには多数の高調波が必要で、計算コストが3次元に近づく。ノズル荷重のような局所問題は3次元サブモデリングのほうが効率的だ。
まとめ
軸対称の数値手法、整理します。
要点：
- $2\pi r$ の重みが剛性マトリクスに入る — 平面要素との本質的な違い
- $r = 0$ 上の節点は $u_r = 0$ — 回転軸の拘束を忘れない
- 荷重はリング状に作用する — 集中荷重の意味が平面問題と異なる
- フーリエ展開で非軸対称荷重に対応 — ただし局所荷重は3Dが効率的
- ハイブリッド要素 — 非圧縮材料（ゴム）の軸対称解析にも必要
軸対称要素の数値積分
軸対称要素の剛性行列はr方向の半径重みを含むため、通常の平面要素と異なるガウス積分則が必要だ。NASTRANのCTRAXA要素は2×2ガウス点を採用し、原点（r=0）近傍での数値不安定を回避する特別処理を1972年版から実装している。
線形要素（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 Div. 2では、Design by Analysisの手法としてFEM軸対称解析が広く使われている。
圧力容器の軸対称解析
圧力容器のどの部分を軸対称で解析しますか？
| 部位 | 軸対称で解析可能か | 注意点 |
|---|---|---|
| 円筒胴 | ○ | 基本的なケース |
| 半球鏡板 | ○ | 胴との接続部の不連続応力 |
| 半楕円鏡板 | ○ | ナックル部の高応力 |
| フランジ | ○ | ボルト荷重を等価リング荷重で |
| ノズル（軸方向） | △ | 軸対称なら○、偏心ノズルは× |
| 支持脚・サドル | × | 非軸対称。3Dが必要 |
鏡板と胴の接続部が重要なんですね。
そう。半楕円鏡板や皿形鏡板のナックル部（曲率が変化する部分）に不連続応力が発生する。この応力は膜理論では計算できず、FEMの軸対称解析で初めて正確に評価できる。
実務的なポイント：
- ナックル部のメッシュは板厚の1/4以下の要素サイズ
- 板厚方向に最低4要素（二次要素なら2要素）
- 溶接部のミスアライメント（目違い）は初期不整として導入
Oリング溝の接触解析
軸対称で接触問題も解けますか？
Oリングの圧縮変形は典型的な軸対称接触問題だ。ゴム（非圧縮材料）のOリングが溝に押し込まれて変形する。
ポイント：
- ゴムの超弾性構成則 — Mooney-Rivlin, Ogden等。線形弾性では不正確
- ハイブリッド要素（CAX4H） — 非圧縮材のロッキング対策
- 接触面の定義 — Oリングと溝面の接触。摩擦あり
- 大変形解析 — Oリングの圧縮率が10〜30%と大きい
2次元で接触を解けるなら、3次元より圧倒的に速いですね。
Oリングの設計では100通り以上のパラメータ（溝深さ、溝幅、Oリング径、材料硬度）を評価することがある。3次元では不可能だが、軸対称なら各ケース数秒で解ける。
回転体の遠心力解析
タービンディスクの遠心力解析も軸対称ですか？
ディスク自体は軸対称だが、ブレードは非軸対称。実務では：
1. ディスクを軸対称で解析 — 遠心力 $\omega^2 r \rho$ を体積力として入力
2. ブレードの遠心力を等価リング荷重として入力 — ブレード本数×遠心力 / $2\pi r$
3. クリティカルな部位（ボルト穴等）は3次元サブモデルで詳細解析
ブレードを等価リング荷重にするのは近似ですが、ディスクの全体応力には十分な精度ですか？
ディスクのボア部（中心穴の周辺）やウェブ部の応力は軸対称で十分正確に出る。ブレード取付部のダブテール溝はローカルに非軸対称だから、3次元のサブモデルが必要。
実務チェックリスト
軸対称解析のチェックリストをお願いします。
「結果が全周分か断面分か」はどういう意味ですか？
軸対称解析の反力は断面上の値だが、全周に換算するには $2\pi r$ を掛ける。例えば $r = 100$ mm の節点の反力が $F_r = 50$ N なら、全周の力は $50 \times 2\pi \times 0.1 = 31.4$ N…ではなく、軸対称要素では通常全周分の値がすでに含まれている。ソルバーによって異なるから、マニュアルで確認すること。
ターボポンプ解析の実務
H-IIAロケットのLE-7Aエンジンのターボポンプ解析では、軸対称モデルが初期設計段階で多用された。3次元モデル構築前に軸対称で板厚・圧力・回転荷重を試計算することで、JAXA内部では設計工数を約40%削減できたと報告されている。
解析フローのたとえ
解析の流れは、実は料理とそっくりです。まず材料を買い出し（CADモデルの準備）、下ごしらえをして（メッシュ生成）、火にかけて（ソルバー実行）、最後に盛り付ける（後処理で可視化）。ここで大事な問いかけ——料理で一番失敗しやすい工程はどこでしょう？ 実は「下ごしらえ」なんです。メッシュの品質が悪いと、どんなに優秀なソルバーを使っても結果はめちゃくちゃになります。
初心者が陥りやすい落とし穴
あなたはメッシュ収束性を確認していますか？ 「計算が回った＝結果が正しい」と思っていませんか？ これ、実はCAE初心者が最も陥りやすい罠です。ソルバーは与えられたメッシュで「それなりの答え」を必ず返します。でもメッシュが粗すぎれば、その答えは現実から大きくずれている。最低3段階のメッシュ密度で結果が安定することを確認する——これを怠ると「コンピュータが出した答えだから正しいはず」という危険な思い込みに陥ります。
境界条件の考え方
境界条件の設定は、試験の「問題文を書く」のと同じです。問題文が間違っていたら？ どんなに正確に計算しても答えは間違いますよね。「この面は本当に完全固定なのか」「この荷重は本当に一様分布なのか」——現実の拘束条件を正しくモデル化することが、実は解析全体で最も重要なステップだったりします。
ソフトウェア比較
軸対称解析のツール
軸対称解析に特に強いツールはありますか？
汎用FEMは全て軸対称に対応しているが、圧力容器設計では専用ツールが便利だ。
圧力容器専用ツール
| ソフト | 特徴 |
|---|---|
| PVElite | 圧力容器のコード計算（ASME, EN 13445等）。FEMではなくルールベース |
| Compress | ASME Section VIIIに準拠した設計計算 |
| Nozzle Pro | ノズル接続部の応力解析。軸対称FEM内蔵 |
| FEPipe | 配管の柔軟性解析と局所応力評価 |
Nozzle Proは軸対称FEMが内蔵されているんですね。
Nozzle ProはWRC 297/537の計算に加えて、軸対称FEMで不連続応力を自動計算する。ASME Div. 2の応力分類（膜応力、曲げ応力）も自動で行う。圧力容器エンジニアにとっては汎用FEMより効率的だ。
汎用FEM
| 観点 | Nastran | Abaqus | Ansys |
|---|---|---|---|
| 軸対称要素 | CQUADX4/8 | CAX4/8, CAXA | PLANE182/183 |
| フーリエ展開 | ハーモニック要素 | CAXA（制限あり） | PLANE25 |
| 非圧縮材の軸対称 | 限定的 | CAX4H, CAX8RH | Hybrid対応 |
| 接触（軸対称） | 対応 | 対応（最も柔軟） | 対応 |
Abaqusのフーリエ展開は「制限あり」ですか？
AbaqusのCAXA要素はフーリエ展開の高調波数に制限がある。NastranやAnsysのハーモニック要素のほうがフーリエ展開には柔軟。ただし多くの場合、非軸対称荷重は3次元サブモデルで扱うほうが実務的だ。
選定ガイド
まとめると？
圧力容器は専用ツールのほうがASMEコードとの整合性が取りやすいんですね。
ASME Div. 2の応力分類（一般膜、局所膜、一次曲げ、二次応力）は汎用FEMの結果から手動で分類する必要がある。専用ツールはこれを自動化してくれるから、設計者の負担が大幅に減る。
主要ソルバーの軸対称対応
Abaqusはリビジョン3.0（1984年）から軸対称要素CAX4RをFull integrationとReduced integrationの両方で提供。MSC NastranはCTRIAX6、Abaqus CAX6MはMidside節点付き6角形を採用し、応力集中部の精度で異なる使い分けが現場で定着している。
選定で最も重要な3つの問い
- 「何を解くか」：軸対称解析に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
軸対称解析の先端トピック
軸対称解析の先端研究にはどんなものがありますか？
軸対称自体は成熟した技術だが、応用面で発展がある。
2.5次元解析（フーリエ展開の高度化）
「2.5次元」解析とは？
例えばベアリングの接触問題。内輪と外輪は軸対称だが、荷重（ラジアル荷重）は非軸対称。フーリエ展開で非軸対称の接触力を分解し、各高調波で軸対称の接触解析を解いて重ね合わせる。
3次元の接触解析を2.5次元で代替できれば、計算コストが大幅に下がりますね。
1〜2桁の高速化が期待できる。特にベアリングやギアのような回転体の接触は、この手法の恩恵が大きい。
軸対称から3次元へのサブモデリング
最近の実務では軸対称グローバルモデル → 3次元ローカルモデルのサブモデリングが標準化しつつある。
手順：
1. 軸対称で圧力容器全体を解析
2. ノズル接続部など非軸対称部位の変位を抽出
3. 3次元のサブモデルの境界条件として適用
4. サブモデルで局所応力を詳細評価
AbaqusやAnsysのサブモデリング機能で自動化できますか？
Abaqusの *SUBMODEL キーワードは軸対称→3次元のサブモデリングに対応している。軸対称モデルの $(r, z)$ 座標を3次元の $(x, y, z)$ に自動マッピングする。Ansysでも Workbench のサブモデリング機能で同様のことが可能。
水素脆化・拡散連成
圧力容器関連の先端トピックは？
水素貯蔵容器（高圧水素タンク）の設計が注目分野だ。水素が材料中に拡散して脆化を引き起こす現象（水素脆化）を応力解析と連成する。
軸対称の応力場→水素の応力駆動拡散→脆化による材料劣化→応力の再配分…という連鎖をFEMで追跡する。圧力容器の大部分は軸対称だから、この連成解析を2次元で効率的に実施できる。
水素社会の実現に直結する研究ですね。
FCV（燃料電池車）の水素タンク（700気圧）やパイプラインの水素輸送で、軸対称のFEM連成解析が設計ツールとして使われている。
まとめ
軸対称の先端トピック、まとめます。
軸対称は「古い手法」ではなく、最先端の連成解析を効率的に実施するための基盤だ。
フーリエ級数軸対称拡張
非軸対称荷重を持つ構造に対し、荷重をフーリエ級数展開して各調和成分を独立した軸対称問題として解く手法が1970年代にWilsonらにより確立された。ADINA-Aソルバーはこの手法で円筒タンクの地震応答解析を行い、3次元解析比で計算時間を1/20に短縮した。
トラブルシューティング
軸対称解析のトラブル
軸対称解析でよくあるトラブルを教えてください。
軸対称特有の盲点がいくつかある。
回転軸上の拘束忘れ
$r = 0$ 上の節点で $u_r = 0$ を忘れるとどうなりますか？
回転軸が横に動くため、応力分布が非物理的になる。特異剛性マトリクスのエラー（Nastranの FATAL 2012 等）が出ることもある。
多くのソルバーは軸対称要素の $r = 0$ 上の節点に自動的に $u_r = 0$ を設定するが、手動で設定が必要なソルバーもある。必ず確認すること。
フープ応力の見落とし
結果で $\sigma_r$ と $\sigma_z$ だけ見て $\sigma_\theta$ を見忘れることがあります。
軸対称では $\sigma_\theta$（フープ応力）が最大応力になることが多い。内圧を受ける円筒では $\sigma_\theta > \sigma_z > \sigma_r$ だ。$\sigma_\theta$ を見落とすと最大応力を過小評価する重大なミスになる。
von Mises応力をチェックする場合も、$\sigma_\theta$ は自動的に含まれるが、応力分類（ASME Div. 2の膜/曲げ分離）では個別に確認する必要がある。
荷重の入力単位ミス
軸対称の荷重入力で注意すべきことは？
ソルバーによって荷重の解釈が異なる：
- Abaqus — 集中力は全周分（$2\pi r$ 分）。$F = 100$ N を入力すると、全周で100 N
- Nastran — FORCE カードは全周分。$F = 100$ N = 全周100 N
- Ansys — ノードの集中力は全周分
全ソルバーで全周分なんですね。
「リング荷重」という意味だ。内圧100 MPaを面として与える場合はどのソルバーでも正しく処理されるが、集中力を与えるときは全周の合計力として入力することを忘れないこと。
$r = 0$ 付近の応力集中
回転軸付近で応力が異常に大きくなることがあります。
$r \to 0$ で $\varepsilon_\theta = u_r / r$ の分母がゼロに近づくため、数値的な精度低下が起きることがある。特にメッシュが粗いと $r = 0$ 付近の要素で応力がばらつく。
対策：
- $r = 0$ 付近のメッシュを十分細かくする
- $r = 0$ 上の応力は「参考値」として扱い、設計判断には使わない
- 必要なら3次元モデルで $r = 0$ 付近を検証
ラメの理論解との比較
FEMの結果がラメの理論解と合わないときは？
確認項目：
1. 開口端か密閉端か — ラメの問題は平面ひずみ（無限長円筒）を仮定。密閉端では端部効果で $\sigma_z$ が変わる
2. 端部条件 — FEMモデルの端部（$z$ 方向の境界条件）が理論と一致するか
3. メッシュ密度 — 肉厚方向に最低4要素（二次要素なら2要素）
4. 薄肉近似との混同 — 薄肉公式 $\sigma_\theta = pD/2t$ はラメの近似。$D/t < 20$ では差が出る
薄肉公式とラメの式の差が出るのは $D/t$ がいくつ以下ですか？
$D/t < 20$ で差が5%以上になる。$D/t < 10$（厚肉）ではラメの式が必須。FEMは当然どちらも正確に解けるが、検証でどちらの理論解と比較するかを間違えないこと。
まとめ
軸対称のトラブル対処、整理します。
フープ応力の見落としが一番怖いですね。$\sigma_\theta$ が最大応力なのに気づかないと、設計が危険側になる。
圧力容器の事故の多くはフープ応力による破裂だ。$\sigma_\theta$ は軸対称解析で最も注目すべき成分であることを常に意識してほしい。
軸上特異点のトラブル
軸対称解析でr=0の対称軸上に節点を配置すると、フープ応力σθがσrと一致しない「疑似応力」が生じる場合がある。ANSYS PLANE25要素では、Ux=0の境界条件を対称軸全節点に強制しないと1980年代初期バージョンで最大15%の誤差が報告されていた。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——軸対称解析の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告