structural-analysis-workflow

ポンプ部品の構造解析ワークフロー

このスキルでは、流体荷重を受けるポンプ部品、特にケーシングとインペラの有限要素解析（FEA）を実施するための体系的な手順を説明します。

FEA ワークフロー

1. ジオメトリと材料の定義

ジオメトリの準備：

• ポンプのケーシングまたはインペラの CAD ジオメトリをインポートするか、作成します
• 構造に関係しない形状（小さなフィレット、ボルト穴、外観上の形状）を削除して、ジオメトリを簡略化します
• 重要領域を特定します：高応力領域、圧力境界、取り付け部
• ポンプ軸に沿った座標系を定義します（通常、Z 軸をシャフトに沿って配置します）

材料の定義：

• すべての部品に材料特性を割り当てます
• 必要に応じて、等方性または異方性の特性を定義します
• 熱の影響がある場合は、温度依存の特性を入力します
• 鋳造品または鍛造品の材料配向を考慮します

重要なパラメーター：

• ヤング率（E）
• ポアソン比（ν）
• 降伏強さ（σ_y）
• 引張強さ（σ_u）
• 密度（ρ）
• 熱膨張係数（α）

2. 荷重の適用

圧力荷重：

• ケーシングの流体接触面に内圧を適用します
• 設計点の静圧または最大運転圧力を使用します
• インペラには、CFD 解析から得た圧力分布を適用します
• 疲労解析では圧力脈動を考慮します

遠心力（インペラ）：

• インペラ本体に回転速度を適用します
• 遠心力：F_c = mω²r
• 方向と大きさが正しいことを確認します
• 過速度状態を考慮します（通常、定格速度の 120%）

水力学的荷重：

• インペラに作用するラジアルスラスト
• 圧力差によるアキシャルスラスト
• 流体の運動量変化による接線力
• 設計点外の条件における不釣り合い力

熱荷重：

• 熱解析から得た温度分布
• 部品間の差動膨張
• 起動・停止時の温度勾配

3. 境界条件（拘束）

ケーシングの拘束：

• 取り付け脚に固定支持を設定します
• 該当する場合は対称面を設定します
• 配管接続部（ノズル）では、通常、熱膨張を許容します
• アンカーポイントは、すべての方向に完全拘束します

インペラの拘束：

• ハブをシャフトに拘束します（円筒拘束）
• インペラとシャフトの間に接触を設定します（結合接触または摩擦接触を使用する場合があります）
• セクターのみを解析する場合は、対称条件を設定します

接続の種類：

• 結合：相対運動なし（溶接接合）
• 摩擦：クーロン摩擦モデル（焼きばめ）
• 分離なし：滑りは可能ですが、分離はできない接触
• 摩擦なし：自由に滑る接触

4. メッシュ生成

要素の種類：

• 四面体（Tet10）：汎用、自動メッシュ生成
• 六面体（Hex20）：高精度、構造化領域
• くさび（プリズム）：四面体メッシュと六面体メッシュの間の遷移
• シェル要素：ケーシングの薄肉部

メッシュ品質の基準：

• 要素アスペクト比 < 5:1（理想は < 3:1）
• スキュー < 0.85
• ヤコビアン比 > 0.6
• 厚さ方向に最低 3 要素
• 応力集中部（フィレット、穴、溶接部）ではメッシュを細かくします

メッシュの細分化：

• 重要な形状の周囲に影響範囲を設定します
• 以下の箇所に局所メッシュ制御を設定します：
  • フィレット半径
  • ノズルとシェルの接合部
  • インペラアイ（入口）
  • インペラ羽根とハブの接合部
  • ウェアリングの隙間
• メッシュ収束性の検討を行います（メッシュ密度を上げたときに結果が収束することを確認します）

推奨要素サイズ：

• インペラ羽根の厚さ：最低 4 要素
• ケーシング壁の厚さ：最低 3 要素
• フィレット半径：円弧に沿って最低 6 要素
• 全体の要素サイズ：代表寸法の 5～10%

5. ソルバーの選択

線形静解析：

• 使用する場合：変形が小さく、材料が弾性で、時間依存の影響がない場合
• 適している用途：初期設計の検証、安全率の計算
• 計算が速く、パラメトリックスタディに適しています

非線形解析：

• 使用する場合：大変形、材料の塑性、接触の非線形性がある場合
• 必要な用途：塑性崩壊解析、接触圧力分布の評価
• 計算に時間がかかり、段階的な荷重が必要です

モーダル解析：

• 固有振動数と振動モードを特定します
• ポンプの運転速度や高調波との共振を避けます
• 確認項目：1次危険振動数 > 運転振動数 × 1.5

調和／過渡解析：

• 翼通過周波数における圧力脈動
• 起動時および停止時の過渡現象
• ウォーターハンマー現象

疲労解析：

• 高サイクル疲労（> 10^4 サイクル）：インペラー、シャフト
• 低サイクル疲労：熱サイクル、起動／停止
• S-N 曲線またはひずみ寿命法を使用します

ソルバー設定：

• 収束基準：通常、力／モーメントの収束率 0.1％
• 大変形：薄い羽根を持つインペラーでは ON
• 接触設定：初回の解析ではプログラム制御
• 安定化：非常に薄い形状では必要になる場合があります

6. 結果の解釈

主な結果：

フォン・ミーゼス応力（σ_vm）：

• 延性材料の相当応力
• 材料の降伏強度と比較します
• σ_vm = √[(σ₁-σ₂)² + (σ₂-σ₃)² + (σ₃-σ₁)²] / √2
• 最も一般的に使用される破損基準です

主応力（σ₁、σ₂、σ₃）：

• 最大および最小の垂直応力
• 脆性材料で重要です
• 最大主応力を極限強度と比較します

変形：

• 全変形：変位の全体的な大きさ
• 方向別変形：すきま（インペラーとケーシング間）を確認します
• 重要なすきま：ウェアリングのすきま、スラスト軸受の軸方向すきま

安全率：

• 安全率（FOS）= σ_許容 / σ_実際
• 色分け：FOS < 1（赤／破損）、FOS > 1（安全）

副次的な結果：

せん断応力：

• 溶接部、キー溝、シャフト接続部で重要です
• 最大せん断応力説（トレスカ基準）

ひずみ：

• 弾性ひずみと塑性ひずみ
• ひずみエネルギー密度
• 延性破壊に対する相当塑性ひずみ

接触圧力：

• 焼きばめ接合面
• 軸受面
• シール面

重要なレビュー項目：

1. 応力集中は幾何形状の特徴と一致していますか？
2. 境界条件は現実的ですか（剛体運動がないこと）？
3. 変形パターンは物理的に妥当ですか？
4. 最大応力は想定される位置に発生していますか？
5. 解は力の釣り合いを満たしていますか（反力の合計 = 適用荷重）？

7. 安全率の評価

設計規格と標準：

• ASME Section VIII Division 1：圧力容器（ケーシング）
• ASME Section VIII Division 2：高圧向け、FEA に基づく設計
• API 610：石油産業向け遠心ポンプ
• ISO 13709：石油、石油化学、天然ガス産業向け遠心ポンプ

必要な安全率：

静荷重：

• 降伏に対する安全率：FOS ≥ 1.5（一般機械）
• 降伏に対する安全率：FOS ≥ 2.0（ASME に基づく圧力容器）
• 極限強度に対する安全率：FOS ≥ 2.5（静荷重）
• 極限強度に対する安全率：FOS ≥ 4.0（動荷重／衝撃荷重）

疲労荷重：

• 無限寿命設計：応力 < 耐久限度
• 有限寿命：累積損傷理論（マイナー則）を使用します
• 疲労寿命に対する安全率：通常 2.0 以上

安全率の計算：

延性材料の場合：

FOS_yield = σ_y / σ_vm
FOS_ultimate = σ_u / σ_vm
Design FOS = min(FOS_yield, FOS_ultimate)

脆性材料の場合：

FOS = σ_u / σ₁ (based on maximum principal stress)

材料ごとの調整：

• 鋳鉄：許容応力を 25％低減します
• 溶接継手：継手効率係数（0.7-1.0）を適用します
• 高温：クリープ低減係数を適用します
• 腐食環境：腐食代を追加します

ポンプへの適用

ポンプケーシングの圧力保持

解析の目的：

• ケーシングが最大許容使用圧力（MAWP）に耐えられることを検証します
• ノズルとシェルの接合部の応力を確認します
• フランジのシール性能を評価します
• ケーシングの変形がウェアリングのすきまに及ぼす影響を評価します

重要領域：

• ボリュートのカットウォーター（圧力勾配が最大となる点）
• ノズル補強パッド
• ケーシングの分割面（水平または垂直分割）
• 取付脚の接合部

荷重ケース：

• 最大運転圧力
• 水圧試験圧力（通常、設計圧力の1.5倍）
• 圧力＋熱膨張
• 圧力＋配管荷重（ノズルにかかる力とモーメント）

合格基準：

• ミーゼス応力 ＜ 許容応力
• フランジ面の変形 ＜ ガスケットの許容範囲
• ボルト穴で降伏が発生しないこと
• 最大変形 ＜ ケーシング直径の0.5％

一般的な故障モード：

• カットウォーター部でのケーシング破裂
• 荷重の不均一によるフランジ漏れ
• 溶接止端部の亀裂
• 圧力サイクルによるボルトの破損

インペラーの遠心応力

解析の目的：

• ブレード、シュラウド、ハブの最大応力を計算します
• ブレードとシュラウドの接合部の応力を評価します
• すきまを維持できるよう、たわみを確認します
• ブレード通過周波数での疲労寿命を評価します

重要領域：

• ブレード後縁（曲げモーメントが最大）
• ブレードとハブのフィレット（応力集中部）
• ブレードとシュラウドの接合部（クローズドインペラーの場合）
• ハブ穴（キー溝による応力集中部）

荷重ケース：

• 定格回転速度での遠心力
• 過速度（定格の120％）での遠心力
• 遠心力と油圧の複合荷重
• 流体温度による熱荷重

応力成分：

遠心応力（円周応力）：

σ_hoop = ρω²r² (for thin disk)
where:
  ρ = material density (kg/m³)
  ω = angular velocity (rad/s)
  r = radius (m)

ブレードの曲げ応力：

• ブレード両側の圧力差によって発生します
• ブレード根元（ハブ接合部）で最大になります
• ブレード長と圧力差が大きいほど増加します

合格基準：

• 定格回転速度でのミーゼス応力 ＜ 0.5 × σ_y
• 過速度でのミーゼス応力 ＜ 0.67 × σ_y
• ブレード先端のたわみ ＜ すきまの50％
• 運転速度範囲の±20％以内に固有振動数がないこと

設計の最適化：

• ハブ直径を大きくして応力を低減します
• ハブからシュラウドに向かってブレード厚さを薄くします
• ブレードとハブの接合部にフィレットを追加します（通常、r ≥ 3mm）
• 強度と密度の比が高い材料を使用します

シャフトのたわみ

解析の目的：

• ラジアル荷重とアキシャル荷重によるシャフトのたわみを計算します
• ベアリングのアライメントが維持されることを確認します
• カップリングのアライメント公差を満たすことを確認します
• シャフトの危険速度と運転速度を確認します

重要領域：

• ベアリング間の中央部（たわみが最大）
• インペラー位置（ウェアリングのすきまに影響）
• カップリング位置（ミスアライメントは振動の原因になります）
• ベアリングジャーナル（段差部の応力集中箇所）

荷重ケース：

• ラジアル方向の油圧推力（締切運転時に最大）
• アキシャル方向の油圧推力
• インペラー重量＋流体重量
• 熱膨張

ラジアル方向の油圧推力：

F_radial = K_r × ρ × g × H × D₂ × b₂
where:
  K_r = radial thrust coefficient (0.2-0.4 for single volute)
  H = head (m)
  D₂ = impeller diameter (m)
  b₂ = impeller outlet width (m)

合格基準：

• シャフトの最大たわみ ＜ シャフトスパンの0.0005倍
• インペラー位置でのたわみ ＜ ウェアリングのすきまの25％
• シャフト応力 ＜ 0.3 × σ_y（耐疲労性を確保するため）
• 第1危険速度 ＞ 最大運転速度の1.5倍

たわみの計算：

• 予備推定には梁理論を使用します
• 複雑な形状や荷重分布にはFEAを使用します
• 高速ポンプではジャイロ効果を考慮します

疲労解析

解析の目的：

• 疲労寿命をサイクル数または年数で予測します
• き裂が発生しやすい箇所を特定します
• 応力集中係数を評価します
• 無限寿命または安全な有限寿命を実現する設計を行います

疲労が問題となる箇所：

• インペラ羽根の後縁（圧力脈動）
• 羽根とハブのフィレット部（応力集中）
• シャフトのキー溝（応力集中係数 ~3）
• ベアリング位置にあるシャフトの段差部
• ケーシングの溶接止端部

荷重サイクル：

• 起動／停止サイクル：低サイクル疲労
• 羽根通過周波数：高サイクル疲労
• 圧力脈動：高サイクル疲労
• ローターの不釣り合い：高サイクル疲労

疲労解析手法：

S-N 曲線法（高サイクル疲労）：

• N > 10⁴ サイクルに適用します
• 材料の S-N 曲線（応力と破損までのサイクル数）を使用します
• 平均応力補正（Goodman または Soderberg）を適用します
• Miner の法則を用いて累積損傷を計算します

ひずみ寿命法（低サイクル疲労）：

• N < 10⁴ サイクルに適用します
• 塑性ひずみ範囲を使用します
• Coffin-Manson の式を使用します
• 起動／停止の解析に必要です

応力集中係数：

• 鋭角部：K_t = 2-3
• キー溝：K_t = 2-3
• フィレット：K_t = 1.5-2.5（半径に依存）
• ねじ部：K_t = 2-4

合否判定基準：

• 無限寿命：交番応力 < 疲労限度
• 有限寿命：計算寿命 > 2 × 必要寿命
• 応力振幅に対する安全率 ≥ 1.5
• 寿命に対する安全率 ≥ 2.0

平均応力の影響：

Goodman correction:
σ_a / S_e + σ_m / σ_u = 1 / FOS

where:
  σ_a = alternating stress amplitude
  σ_m = mean stress
  S_e = endurance limit
  σ_u = ultimate tensile strength

材料選定基準

ポンプケーシングの材料

鋳鉄（ASTM A48、A278）：

• 用途：低圧の水用途、非腐食性流体
• 特性：σ_u = 150-400 MPa、低コスト、鋳造性に優れています
• 制約：脆性があり、疲労抵抗性が低く、<250°C に制限されます
• 使用例：自治体の水道、HVAC、圧力 < 20 bar

ダクタイル鋳鉄（ASTM A536）：

• 用途：一般産業用途、中程度の圧力
• 特性：σ_y = 275-550 MPa、σ_u = 400-800 MPa、機械加工性に優れています
• グレード：65-45-12、80-55-06（σ_u-σ_y-伸び）
• 使用例：プロセス水、油の移送、圧力 < 40 bar

炭素鋼（ASTM A216 WCB）：

• 用途：高圧・高温用途
• 特性：σ_y = 250 MPa、σ_u = 485 MPa、溶接性に優れています
• 温度範囲：-29°C～400°C
• 使用例：ボイラー給水、発電所、圧力 < 100 bar

ステンレス鋼 316（ASTM A743 CF-8M）：

• 用途：腐食性流体、海水、化学薬品
• 特性：σ_y = 205 MPa、σ_u = 485 MPa、耐食性に非常に優れています
• 温度範囲：-196°C～400°C
• 使用例：化学処理、海洋、食品・飲料

青銅（ASTM B584）：

• 用途：海水用途、小型ポンプ
• 特性：耐食性に優れ、中程度の強度
• 使用例：海洋、海水淡水化、ポンプ内部部品

インペラの材料

ダクタイル鋳鉄（ASTM A536）：

• 強度とコストのバランスに優れています
• 密度：7,100 kg/m³
• σ_y = 275 MPa（最小値）
• 使用例：一般用途、水

316 ステンレス鋼：

• 耐食性に優れています
• 密度：8,000 kg/m³
• σ_y = 275 MPa（鋳造）、290 MPa（インベストメント鋳造）
• 使用例：化学、食品、医薬品

二相ステンレス鋼（CD4MCu）：

• 強度と耐食性に優れています
• σ_y = 450 MPa
• 高価ですが、長寿命です
• 使用例：海水、強腐食性の化学薬品

ニッケルアルミニウム青銅（NAB）：

• キャビテーション抵抗性に非常に優れています
• σ_y = 240 MPa
• 海水用途に最適です
• 使用例：海洋、海水淡水化

チタン（Ti-6Al-4V）：

• 強度重量比が最も高い材料です
• 密度：4,430 kg/m³
• σ_y = 880 MPa
• 非常に高価です
• 使用例：高速用途、航空宇宙

シャフト材料

炭素鋼（AISI 1045、4140）：

• 汎用
• σ_y = 400-650 MPa
• 被削性に優れています
• 腐食防止が必要です

ステンレス鋼 416、17-4PH：

• 腐食環境向け
• σ_y = 520-1,170 MPa（17-4PH 熱処理済み）
• 疲労抵抗に優れています
• 高荷重には 17-4PH を使用します

合金鋼（AISI 4340）：

• 高出力用途向け
• σ_y = 860-1,380 MPa（熱処理済み）
• 疲労特性に非常に優れています
• 保護コーティングが必要です

材料選定プロセス

ステップ 1：使用条件

• 流体の種類（腐食性、研磨性、清浄）
• 温度範囲
• 圧力レベル
• 環境への露出

ステップ 2：性能要件

• 必要な強度（応力レベル）
• 疲労抵抗（繰り返し荷重）
• 耐摩耗性（滑り接触）
• 熱膨張の適合性

ステップ 3：製造上の考慮事項

• 鋳造と棒材からの機械加工の比較
• 溶接性の要件
• 熱処理の実現可能性
• 検査要件（NDT）

ステップ 4：費用対効果の分析

• 材料費
• 加工費
• 想定耐用年数
• 保守費用

ステップ 5：規格への適合

• ASME の許容応力
• 温度制限
• 衝撃試験要件（低温）
• 腐食代

応力集中

理論応力集中係数（K_t）

応力集中係数は、局所的な最大応力と公称応力の関係を示します。

σ_max = K_t × σ_nominal

ポンプにおける一般的な応力集中

フィレット：

K_t = 1 + 2√(d/r)
where:
  d = smaller diameter
  r = fillet radius

設計ルール：r ≥ 0.1 × d で K_t ≤ 2.0

段差（直径の変化）：

• r/d = 0.05：K_t ≈ 2.5
• r/d = 0.10：K_t ≈ 2.0
• r/d = 0.20：K_t ≈ 1.6

板の円形穴：

• 単一の穴：K_t = 3.0
• 複数の穴：相互作用により K_t が増加します

シャフトのキー溝：

• プロファイルキー溝：K_t = 2.0-2.5
• スレッドキー溝：K_t = 2.5-3.0
• エンドミル加工キー溝：K_t = 3.0-3.5

ねじ：

• 並目ねじ：K_t = 2.2-3.0
• 細目ねじ：K_t = 2.8-4.0
• 応力除去溝により K_t が低下します

溶接部：

• 溶接止端：K_t = 2.0-3.0
• 溶接ルート：K_t = 2.5-4.0
• 平滑に仕上げると K_t が 30-50% 低下します

応力集中を緩和する方法

フィレット半径を大きくする：

• 半径を 2 倍にする → K_t が 20-30% 低下します
• FEA のパラメトリックスタディで最適化します
• 一定半径より滑らかにつないだ半径のほうが優れています

逃げ溝を追加する：

• 肩部の変化箇所にアンダーカットを設けます
• 致命的断面の応力を低減します
• 軸受のはめあい部によく用いられます

鋭角な角をなくす：

• すべての内角に半径を設けます
• 最小半径 = 2 × 肉厚
• 製造可能性を確認します

形状の変化箇所を移動する：

• 公称応力の高い領域から応力集中箇所を離します
• 曲げモーメントの大きい位置から直径の変化箇所を離します

局所的に材料を追加する：

• 応力集中箇所の局所的な厚さを増やします
• ボスまたは補強パッドを追加します
• 滑らかにつなぎます

表面処理：

• ショットピーニングにより圧縮残留応力が生じます
• 疲労寿命が 2-4 倍向上します
• シャフトやブレードに効果的です

疲労切欠係数（K_f）

疲労解析では、有効応力集中は低減されます。

K_f = 1 + q(K_t - 1)
where:
  q = notch sensitivity (0 to 1)
  K_t = theoretical stress concentration factor

切欠感度は以下に依存します。

• 材料の延性（延性材料ほど低くなります）
• 切欠半径（半径が大きいほど低くなります）
• 材料の結晶粒径

代表値：

• 鋼、r = 0.5 mm：q ≈ 0.9
• 鋼、r = 2.0 mm：q ≈ 0.95
• アルミニウム：q ≈ 0.7-0.8

FEA による検証

局所メッシュの細分化：

• 切欠部のメッシュサイズ < r/4（フィレット半径）
• フィレット円弧に沿って少なくとも 6-8 個の要素を使用します
• メッシュを細分化して応力の収束を確認します

線形化：

• 断面厚さ方向の応力を抽出します
• 膜応力成分と曲げ応力成分を分離します
• 圧力容器規格への適合確認に使用します

サブモデリング：

• 粗いメッシュの全体モデルを作成します
• 応力集中箇所の周囲に細かいメッシュを用いた局所モデルを作成します
• 全体モデルから局所モデルへ変位境界条件を適用します

---

使用ガイドライン

この構造解析ワークフローは、反復して適用する必要があります。

1. 予備設計： 手計算と簡略化したFEAを使用します
2. 設計の改良： 正確な形状と荷重を用いて詳細なFEAを実施します
3. 設計検証： すべての荷重ケースと安全率を考慮した最終解析を行います
4. 疲労評価： 繰り返し荷重が大きい場合に実施します
5. 文書化： 応力レベル、安全率、規格への適合状況を報告します

FEAの結果は、常に以下と照合して検証してください。

• 可能な場合は理論計算
• 現場での実績がある過去の設計
• 業界標準およびガイドライン
• 実施可能な場合は物理試験

ご留意ください： FEAは技術者の判断を支援するツールであり、その代わりになるものではありません。結果の妥当性を常に確認し、境界条件が実際の設置条件および運転条件を反映していることを確かめてください。