ボルト締結体の線形解析
理論と物理
ボルト締結の力学
先生、ボルト締結体のFEM解析はどう行うんですか？
ボルト締結は構造解析で最も一般的な接合形式だが、FEMのモデル化は意外と難しい。ボルトのプリテンション（初期張力）、被締結体の圧縮、外力による軸力変動、接触面の滑り…全て非線形要素を含む。
ボルト締結の基本力学
ボルト締結の力学の基本はボルト-被締結体の直列ばねモデルだ。
ボルトの軸剛性：$k_b = E_b A_b / L_b$
被締結体の圧縮剛性：$k_c$
プリテンション（初期張力）$F_i$ をかけた後に外力 $F_{ext}$ が作用すると：
ここで $\Phi = k_b / (k_b + k_c)$ は内力係数（荷重導入率）。
外力の一部（$\Phi$倍）だけがボルトに追加荷重として入る。残りは被締結体の圧縮が緩和されることで吸収される。
通常のボルト締結では $\Phi = 0.1 \sim 0.3$。つまり外力の10〜30%しかボルトに入らない。これがプリテンションの効果で、ボルトの疲労寿命を大幅に改善する。
プリテンションが重要なのは疲労のためなんですね。
そう。プリテンションなしでは外力がそのままボルトに入り、疲労で破断する。適切なプリテンションにより、ボルトの応力変動幅を大幅に小さくできる。
VDI 2230 ガイドライン
ボルト締結の設計基準はありますか？
VDI 2230の手計算はボルト1本の軸方向問題に限定される。偏心荷重や複数ボルトの相互作用はFEMが必要になる。
FEMでのモデル化レベル
ボルト締結のFEMモデル化は3つのレベルがある：
| レベル | モデル化 | 精度 | コスト |
|---|---|---|---|
| Level 1 | ばね要素＋RBE2/RBE3 | 低（概算） | 低 |
| Level 2 | 梁要素＋接触 | 中 | 中 |
| Level 3 | ソリッド要素＋接触＋ねじ形状 | 高 | 高 |
Level 3でねじの形状まで入れるんですか？
ねじ山の応力集中や座面圧を精密に評価したいなら必要。ただし通常はLevel 1〜2で十分。Level 3は特殊なケース（疲労評価、高温ボルト等）のみ。
まとめ
ボルト締結の理論を整理します。
要点：
- 直列ばねモデル — ボルト剛性 $k_b$ と被締結体剛性 $k_c$ の比率がポイント
- 内力係数 $\Phi = k_b/(k_b+k_c)$ — 外力の10〜30%がボルトに入る
- プリテンションが疲労寿命を改善 — 応力変動幅の低減
- VDI 2230が設計基準 — 手計算は1本の軸方向問題
- 3つのモデル化レベル — 用途に応じた詳細度
$\Phi$ の値がボルト設計の全てを決めるんですね。
$\Phi$ と初期張力 $F_i$ がわかれば、ボルトの最大軸力、座面圧、疲労寿命が計算できる。FEMでもこの2つを正しく評価することが目標だ。
ボルト締め付け理論の歴史
ボルト軸力（プリテンション）の概念を体系化したのは1938年にドイツのVDI（ドイツエンジニア協会）が発行したVDI 2230規格が初出とされる。この規格は現在も改訂を重ね、VDI 2230:2015が世界標準として自動車・航空産業で用いられている。初版から核心にある「締め付け係数」の概念は80年後も変わらない。
各項の物理的意味
- 慣性項（質量項）：$\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でどう設定しますか？
各ソルバーにプリテンション専用の機能がある。
Nastran
```
$ ボルト梁要素にプリテンション
TEMP(INIT) = 100
TEMP(LOAD) = 200
```
温度荷重でプリテンションを模擬する方法が古典的。最近はSOL 400の*BOLT PRELOAD機能。
Abaqus
```
*STEP
*STATIC
*BOLT LOAD
bolt_section, bolt_mid_surface, 50000.
*END STEP
```
*BOLT LOADはボルトの断面に初期張力を直接与える。最も直感的。
Ansys
```
! Workbench: Bolt Pretension荷重
! APDL: PRETS179要素
SLOAD, bolt_section, PRETEN, 50000.
```
Abaqusの*BOLT LOADが一番わかりやすいですね。ボルトの断面を切って初期張力を与える。
Abaqusの実装では、ボルトの中間断面に「プリテンション節点」を追加し、その節点に軸方向の力を与える。力からロック変位に切り替えることで、外力作用後もプリテンションを保持する。
2ステップ解析
プリテンションと外力は同時に与えるんですか？
2ステップで解析するのが基本：
- ボルトに初期張力を与える
- 被締結体が圧縮される
- ボルト軸力 = $F_i$（設計値）
Step 2: 外力
- プリテンションを「ロック」（長さを固定）
- 外力（内圧、自重、風荷重等）を追加
- ボルト軸力 = $F_i + \Phi \cdot F_{ext}$
「ロック」するのはなぜですか？
Step 2でプリテンションの力をそのまま維持すると、外力で変形してもプリテンション力が一定のまま（物理的に不正確）。ロック（変位固定）に切り替えることで、外力による変形でボルト軸力が自然に変動する。これがVDI 2230のモデルと整合する。
接触の必要性
ボルト締結の解析に接触は必要ですか？
本質的には非線形接触問題だ。被締結面のクランプ力が外力で減少し、最終的に離開（separation）する可能性がある。
ただし：
- 離開しないことが確実な場合 → 線形解析（接触なし）で十分
- 離開の可能性がある場合 → 接触定義が必要（非線形解析）
- 被締結面の滑りが重要な場合 → 摩擦付き接触
プリテンションが十分高ければ離開しない？
そう。$F_i > F_{ext}$ なら離開しない（単純な場合）。この条件が設計で確認済みなら線形解析で十分。VDI 2230の手計算でこの確認を行ってからFEMに進むのが効率的だ。
まとめ
ボルト締結の数値手法、整理します。
要点：
- *BOLT LOAD（Abaqus）が最も直感的 — プリテンション断面に力を直接与える
- 2ステップ解析 — Step 1: プリテンション → Step 2: 外力（ロック）
- 接触は離開の可能性がある場合に必要 — 線形解析は離開なしの前提
- VDI 2230の手計算で事前検証 — 離開条件、$\Phi$の確認
FEMでのプリテンション実装
FEMでボルトプリテンションを再現する手法は主に3種類ある。①初期ひずみ法（最も単純）、②温度荷重による熱収縮法、③Abaqusの*PRETENSIONコマンドによる切断面法だ。切断面法は1990年代にHBBSが開発し、締め付けシーケンスを1ステップずつシミュレートできる点が特徴で、航空機のパネル組立解析に採用されている。
線形要素（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要素あたりの計算コストは増えるため、トータルのコスト対効果で判断する。
実践ガイド
ボルト締結解析の実務フロー
ボルト締結のFEM解析を実務で行う手順を教えてください。
1. VDI 2230の手計算 — プリテンション、内力係数、座面圧の概算
2. FEMモデル化レベルの選択 — Level 1（ばね）〜 Level 3（ソリッド）
3. プリテンションステップの設定 — 初期張力を与えて被締結体を圧縮
4. 外力ステップ — 使用荷重を追加
5. 結果の評価 — ボルト軸力、座面圧、被締結面の離開
6. VDI 2230との整合確認 — FEMと手計算の結果を比較
プリテンションの設定値
プリテンションの値はどう決めますか？
VDI 2230に基づく：
$R_p$ はボルトの0.2%耐力（降伏強度）、$A_s$ は有効断面積。10.9のM16ボルトなら $F_i \approx 110$ kN。
降伏強度の90%まで締めるんですか？
高強度ボルトでは一般的。VDI 2230の設計ではボルトを降伏点付近まで締めることでクランプ力を最大化する。ただし繰り返し締め付けではプリテンションの損失（リラクゼーション）も考慮する。
被締結体の圧縮剛性 $k_c$
$k_c$ はFEMから直接求められますか？
求められる。FEMの結果から：
$\delta_c$ はプリテンションステップでの被締結面の圧縮変位。FEMでは自動的にこの変位が出るから、$k_c$ を逆算できる。
VDI 2230の手計算のk_cと比較できますね。
VDI 2230は被締結体を「圧縮円錐」で近似して $k_c$ を求める。FEMはより正確な応力分布を反映するから、手計算と10〜20%程度の差が出ることがある。FEMが正確な場合が多い。
実務チェックリスト
ボルト締結解析のチェックリストをお願いします。
「VDI 2230との整合確認」が最終チェック。手計算とFEMの相互検証ですね。
ボルト締結はFEMだけでなく手計算との組み合わせで設計する分野だ。どちらか一方だけでは不十分。
シリンダーヘッドの締め付け順
自動車エンジンのシリンダーヘッドボルトは通常10本前後を対角順に3段階で締め付ける。ホンダK型エンジン（2001年〜）の開発では、Abaqusのプリテンション解析でガスケット面圧分布を最適化し、従来比で冷機時のヘッドガスケット漏れリスクを約20%低減したとSAE 2003-01-0114に報告されている。
解析フローのたとえ
解析の流れは、実は料理とそっくりです。まず材料を買い出し（CADモデルの準備）、下ごしらえをして（メッシュ生成）、火にかけて（ソルバー実行）、最後に盛り付ける（後処理で可視化）。ここで大事な問いかけ——料理で一番失敗しやすい工程はどこでしょう？ 実は「下ごしらえ」なんです。メッシュの品質が悪いと、どんなに優秀なソルバーを使っても結果はめちゃくちゃになります。
初心者が陥りやすい落とし穴
あなたはメッシュ収束性を確認していますか？ 「計算が回った＝結果が正しい」と思っていませんか？ これ、実はCAE初心者が最も陥りやすい罠です。ソルバーは与えられたメッシュで「それなりの答え」を必ず返します。でもメッシュが粗すぎれば、その答えは現実から大きくずれている。最低3段階のメッシュ密度で結果が安定することを確認する——これを怠ると「コンピュータが出した答えだから正しいはず」という危険な思い込みに陥ります。
境界条件の考え方
境界条件の設定は、試験の「問題文を書く」のと同じです。問題文が間違っていたら？ どんなに正確に計算しても答えは間違いますよね。「この面は本当に完全固定なのか」「この荷重は本当に一様分布なのか」——現実の拘束条件を正しくモデル化することが、実は解析全体で最も重要なステップだったりします。
ソフトウェア比較
ボルト締結のツール
ボルト締結解析に使えるツールを教えてください。
VDI 2230の計算ツール
| ツール | 特徴 |
|---|---|
| BOLT（VDI 2230計算ソフト） | VDI 2230に準拠した手計算の自動化 |
| Bossard Expert | ボルト選定＋VDI 2230計算 |
| Excel自作シート | VDI 2230の式を実装 |
FEMソルバーの比較
| 機能 | Abaqus | Nastran | Ansys |
|---|---|---|---|
| プリテンション | *BOLT LOAD | SOL 400 BOLT | PRETS179 |
| 使いやすさ | ◎（直感的） | △（DMAP/SOL 400必要） | ○（Workbench） |
| 接触との統合 | ◎ | ○ | ○ |
| 疲労後処理 | *OUTPUT + nCode | — + nCode | Fatigue Tool |
Abaqusの*BOLT LOADが最も使いやすいんですね。
Abaqusはボルトプリテンションの設定がワークフローに自然に組み込まれている。Nastranでは歴史的に温度荷重でプリテンションを模擬する手法が使われてきたが、SOL 400の*BOLTオプションで改善されている。
選定ガイド
手計算（VDI 2230）→ FEM → 疲労評価の3段階がボルト設計のフルフローですね。
そう。ボルト設計は「手計算で大まかに → FEMで精密に → 疲労で寿命評価」という階層構造だ。
主要ソルバーのボルト機能比較
Abaqusの*PRETENSIONとANSYSのPSMESH/PSMNUM（Bolt Pretension）コマンド、NX NastranのMPC拘束法はそれぞれ精度と手間が異なる。2018年にNAFEMS Benchmark Challenge #5として公開された標準問題では、Abaqus切断面法がMPC法より締め付けシーケンス依存の面圧分布を5%高精度に再現したと報告された。
選定で最も重要な3つの問い
- 「何を解くか」：ボルト締結体の線形解析に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
ボルト締結の先端研究
ボルト締結解析の最前線を教えてください。
3つの方向がある。
マルチボルト接合の最適化
複数ボルトのフランジ接合では、各ボルトの軸力が一様にならない。プリテンションのばらつき、フランジの変形、荷重の偏心…全てがボルト間の荷重分布に影響する。
FEMでフランジ全体をモデル化し、各ボルトのプリテンション順序まで含めた解析が行われている。ボルトの締め付け順序（スター型 vs. 順次型）で座面圧分布が変わることがFEMで確認されている。
自己緩み（self-loosening）のシミュレーション
振動環境でのボルトの自己緩みは実務上の大問題。JunkersのDIN 65151試験をFEMで再現する研究がある。ボルト-ナット-被締結体のフル接触モデル（ねじ山まで含む）で、横方向振動によるプリテンション減少を追跡する。
ねじ山までモデル化するのはLevel 3ですね。計算コストが大きそう。
1本のボルトで数十万DOF。しかし自己緩みのメカニズム（ナット面の微小すべり → ねじ山面の相対回転 → 軸力低下）を正確に再現するにはこの詳細度が必要だ。
水素脆化によるボルト破損
高強度ボルト（10.9, 12.9クラス）の遅れ破壊（水素脆化）はFEMで残留応力分布と水素拡散の連成解析で予測する研究がある。ボルトの首下フィレット部の三軸応力が水素のトラップサイトになり、亀裂が発生する。
まとめ
ボルト締結の先端研究、まとめます。
ボルト締結は「枯れた技術」ではなく、高信頼性が求められる場面で最先端の解析が必要な分野だ。
熱負荷時のプリテンション変化
高温環境ではボルトと被締結体の熱膨張差によりプリテンション力が変化する。インコネル718製ボルトとアルミ合金フランジの組み合わせでは、200℃の昇温で初期軸力の最大25%が失われることがある。1999年のGE航空エンジン設計チームはこの「軸力弛緩」のFEM解析を標準化し、設計裕度の根拠として社内規格に組み込んだ。
トラブルシューティング
ボルト締結解析のトラブル
ボルト締結のFEM解析でよくあるトラブルを教えてください。
ボルト締結解析は接触＋プリテンションの複合問題だから、トラブルも複合的だ。
プリテンション後のボルト軸力が設定値と異なる
50 kNのプリテンションを設定したのに、ボルト軸力が45 kNしか出ません。
被締結体の変形によりプリテンションの一部が失われている。薄いフランジや軟らかいガスケットでは、ボルトの軸力が設定値より低くなる。
確認方法：
- ボルト要素の軸力を直接出力
- 被締結面の接触圧の合計を確認
- プリテンション後の被締結体の変形量を確認
これは物理的に正しい結果ですか？
正しい場合もある。ガスケット付きフランジでは、ガスケットの圧縮変形でプリテンションが低下するのは現実に起きる。「設定値と軸力が一致しない = エラー」とは限らない。ただし差が大きすぎる（20%以上）場合はモデル化を確認すべき。
接触が収束しない
プリテンションステップで接触が収束しません。
対策：
1. プリテンションを段階的に与える — 一度に全力ではなく、10%, 50%, 100% と段階的に
2. 初期接触のギャップを除去 — 被締結面の初期ギャップがあると接触確立が困難
3. 接触安定化を使用 — Abaqusの *CONTACT STABILIZATION
4. ペナルティ剛性を下げる — 初期は低いペナルティで接触を確立
被締結面の離開判定
外力を加えた後に被締結面が離開しているかどうかは？
接触圧（CPRESS）がゼロの領域 = 離開。接触圧コンターを見れば一目瞭然。
離開が発生している場合：
- プリテンションが不足 → $F_i$ を上げる
- 偏心荷重が大きい → ボルト配置を見直す
- フランジが薄すぎる → フランジ厚さを増やす
ボルトの応力が降伏を超える
FEMでボルトの応力が降伏点を超えています。大丈夫ですか？
線形解析では弾性を仮定しているから、「応力が降伏点を超えた」という結果は「線形解析の仮定が破れた」ことを意味する。弾塑性解析（非線形）に切り替えるか、プリテンションを下げる必要がある。
VDI 2230ではボルトのプリテンションを降伏点の90%に設定するが、これはボルトの軸方向応力が降伏点に近いことを意味する。この状態に外力が追加されると容易に降伏に達する。
まとめ
ボルト締結のトラブル対処、整理します。
手計算（VDI 2230）とFEMの相互検証が全てのデバッグの基本ですね。
ボルト締結は「手計算で予想→FEMで検証→差異の原因を調査」というループが最も効果的なデバッグ方法だ。
プリテンション解析の収束不良
ボルトプリテンション解析でよくある収束失敗は、初期クリアランスを持つ接触面との相互作用によるものだ。Abaqusでは*CONTACT PAIR定義前に*PRETENSION SECTIONを設定すると解析が発散するケースが報告されており、ステップ順序（まず接触確立→次に締め付け荷重）を守ることが2010年以降のマニュアルで明記されている。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——ボルト締結体の線形解析の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告