ギャップ要素
理論と物理
ギャップ要素とは
先生、「ギャップ要素」って何ですか？
ギャップ要素は2点間に隙間（ギャップ）があり、接触したときだけ力を伝達する要素だ。接触問題の簡略化モデルとして使う。
物理的イメージ
ギャップ要素は「1次元の接触要素」だ：
- ギャップが開いている（$\delta < g$） → 力はゼロ
- ギャップが閉じている（$\delta \geq g$） → 圧縮ばね（剛性 $k$）として力を伝達
ここで $\delta$ は2節点間の相対変位、$g$ は初期ギャップ。
「ギャップが閉じたらばね」というシンプルな非線形ですね。
そう。力-変位関係：
用途
| 用途 | 説明 |
|---|---|
| ボルト穴のクリアランス | ボルトが穴に接触したときの荷重伝達 |
| 支承のストッパー | 一定変位以上で接触 |
| 配管のサポート | 片方向のみ支持（リフトオフ） |
| 熱膨張による接触 | 温度上昇でギャップが閉じる |
配管のサポートが片方向のみというのは？
配管が下方にたわむと支持台に乗るが、上方にたわむと支持台から離れる（リフトオフ）。下向きだけ力を伝達し、上向きは自由。これがギャップ要素の典型的な適用だ。
ソルバー別の要素名
| ソルバー | 要素名 | 備考 |
|---|---|---|
| Nastran | CGAP | 方向、ギャップ量、閉合剛性を指定 |
| Abaqus | *GAP / GAPUNI | 1次元ギャップ。ITT要素 |
| Ansys | CONTA178 | 節点間接触要素 |
AbaqusにはGAP要素とは別にITT（Interface）要素もあるんですか？
Abaqusではギャップ要素（GAP）よりも汎用の接触定義（CONTACT PAIR / *GENERAL CONTACT）のほうが柔軟で推奨される。ギャップ要素は簡易的な1次元接触にのみ使う。
ギャップ要素 vs. 接触定義
| 比較 | ギャップ要素 | 面接触定義 |
|---|---|---|
| 自由度 | 1方向のみ | 面全体 |
| 設定の手間 | 少ない | 多い |
| 精度 | 1次元近似 | 正確な接触圧分布 |
| 非線形性 | 弱い | 強い |
| 摩擦 | Nastran CGAPのみ対応 | 完全対応 |
簡易的にはギャップ要素、精密にはフル接触定義、という使い分けですね。
その通り。ギャップ要素は「接触する/しない」の二値判定だけで十分な場合に使う。接触面の圧力分布やすべりが重要なら、汎用の接触定義が必要だ。
まとめ
ギャップ要素の理論を整理します。
要点：
- 隙間が閉じたときだけ力を伝達 — 1次元の接触要素
- 力 = 0（開）or $k(\delta - g)$（閉） — シンプルな非線形
- 配管サポート、ストッパー、クリアランスのモデル化 — 実務で多用
- Nastran CGAP が最も広く使われる — 簡易接触の定番
- 精密な接触にはフル接触定義を使う — ギャップ要素は簡易モデル
ギャップ要素の理論的起源
接触問題をFEMで扱う最初の試みは1963年にHertzの接触理論をマトリクス法で拡張したものだ。ギャップ要素として定式化したのはWilson & Parkesの1972年の論文が先駆けで、2節点間の開閉を「ON/OFFスイッチ」として扱うバイナリ接触法を確立した。この手法は現在のANSYS CONTA171要素の原型となっている。
各項の物理的意味
- 慣性項（質量項）：$\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で統一 |
数値解法と実装
ギャップ要素の非線形解法
ギャップ要素は非線形ですよね。線形解析では使えないんですか？
ギャップの開閉は状態変化だから本質的に非線形。ただし多くのソルバーでは反復法で疑似線形的に処理する。
手順：
1. 全ギャップ要素を「開」と仮定して線形解析
2. 各ギャップの相対変位を確認。ギャップが閉じるべき要素を「閉」に変更
3. 更新された剛性で再度解析
4. 全ギャップの状態が安定するまで反復
NastranのSOL 101（線形静解析）でもCGAPが使えるのはこの反復のおかげですか？
NastranではSOL 106（非線形静解析）でCGAPを使うのが正式だが、SOL 101でもCGAPの開閉反復を行う「線形接触」機能がある。完全な非線形解析よりも速いが、複雑な接触には対応できない。
閉合剛性の設定
ギャップが閉じたときの剛性 $k$ はどう設定しますか？
理想的には「無限大」（完全に硬い接触）だが、数値的にはペナルティ法で有限の大きな剛性を使う。
目安：
- $k \approx 10 \sim 100$ × 接触面の剛性 ($EA/L$ 相当)
- 大きすぎると条件数が悪化して収束困難
- 小さすぎると貫通（penetration）が過大
「ちょうどいい $k$」を見つけるのが難しそうですね。
Abaqusの*CONTACT定義ではペナルティ剛性を自動計算するが、ギャップ要素では手動設定が必要。まず $k$ = 構造剛性の10倍で始めて、貫通量が板厚の1%以下になるよう調整する。
まとめ
ギャップ要素の数値手法、整理します。
要点：
- 反復法でギャップの開閉を収束 — 線形解析の反復として処理
- 閉合剛性は構造剛性の10〜100倍 — 大きすぎても小さすぎてもダメ
- NastranのSOL 101でも簡易接触可能 — ただしSOL 106が正式
- 精密な接触にはフル接触定義を推奨 — ギャップ要素は簡易モデル
ペナルティ法とラグランジュ法
ギャップ要素の数値実装にはペナルティ法とラグランジュ乗数法の2系統がある。ペナルティ法は剛性マトリクスの次元を変えずに実装でき実装が容易だが、ペナルティ係数の選択が精度に直結する。1974年にBatheとWilsonはラグランジュ法の方が条件数が安定することを示し、以後の高精度接触ソルバー開発の指針となった。
線形要素（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要素あたりの計算コストは増えるため、トータルのコスト対効果で判断する。
実践ガイド
ギャップ要素の実務適用
ギャップ要素の実務での使い方を教えてください。
最も一般的な適用は配管解析だ。
配管サポートのモデル化
配管解析ソフト（CAESAR II, AutoPIPE等）では、サポートの種類をギャップ要素で表現する：
| サポート種類 | ギャップモデル | 挙動 |
|---|---|---|
| レストサポート | 下向きのみ支持（ギャップ=0） | 上方リフトオフ可能 |
| ガイド | 横方向のみ支持 | 横方向変位を制限 |
| ストッパー | 一方向のみ支持（ギャップあり） | ギャップ閉じたら接触 |
| スライドサポート | 下向き支持＋摩擦 | 横方向すべり |
配管解析では日常的にギャップ要素を使うんですね。
そう。配管サポートは「片方向支持」が基本だから、ギャップ要素は配管エンジニアにとって最も馴染みのある要素だ。
熱膨張による接触
高温配管が膨張して周囲の構造に接触する問題。初期状態ではギャップがあるが、温度上昇で閉じる。
実務チェックリスト
ギャップ要素のチェックリストをお願いします。
「ギャップの状態が収束しているか」が非線形ならではのチェック項目ですね。
反復が収束しない場合、ギャップが開閉を繰り返している（チャタリング）。閉合剛性を下げるか、ギャップ量を調整する。
鉄道車輪とレールの接触解析
新幹線車輪とレールの接触幅は荷重65kNで約10×15mmの楕円形となる。JR総研が1990年代にAbaqusのギャップ要素でこの接触楕円を解析し、Hertz理論との比較で最大接触圧力の誤差は3%以内、接触幅は5%以内と確認した。この精度がレール摩耗寿命予測モデルの基盤となっている。
解析フローのたとえ
解析の流れは、実は料理とそっくりです。まず材料を買い出し（CADモデルの準備）、下ごしらえをして（メッシュ生成）、火にかけて（ソルバー実行）、最後に盛り付ける（後処理で可視化）。ここで大事な問いかけ——料理で一番失敗しやすい工程はどこでしょう？ 実は「下ごしらえ」なんです。メッシュの品質が悪いと、どんなに優秀なソルバーを使っても結果はめちゃくちゃになります。
初心者が陥りやすい落とし穴
あなたはメッシュ収束性を確認していますか？ 「計算が回った＝結果が正しい」と思っていませんか？ これ、実はCAE初心者が最も陥りやすい罠です。ソルバーは与えられたメッシュで「それなりの答え」を必ず返します。でもメッシュが粗すぎれば、その答えは現実から大きくずれている。最低3段階のメッシュ密度で結果が安定することを確認する——これを怠ると「コンピュータが出した答えだから正しいはず」という危険な思い込みに陥ります。
境界条件の考え方
境界条件の設定は、試験の「問題文を書く」のと同じです。問題文が間違っていたら？ どんなに正確に計算しても答えは間違いますよね。「この面は本当に完全固定なのか」「この荷重は本当に一様分布なのか」——現実の拘束条件を正しくモデル化することが、実は解析全体で最も重要なステップだったりします。
ソフトウェア比較
ギャップ要素のツール
各ソルバーのギャップ要素の特徴は？
配管解析専用ソフト
配管解析では汎用FEMよりも専用ソフトのギャップ要素が使いやすい：
| ソフト | ギャップの扱い |
|---|---|
| CAESAR II | サポートタイプでギャップを自動設定 |
| AutoPIPE | サポート種類とギャップ量をGUI設定 |
| ROHR2 | ドイツ・欧州で標準。ギャップ＋摩擦対応 |
選定ガイド
配管解析では専用ソフト、それ以外はNastranのCGAPが定番ですね。
ギャップ要素は「簡易接触」のためのツール。精密な接触にはフル接触定義を使うべき。用途に応じた使い分けが重要だ。
ソルバー別接触実装の比較
AbaqusのGeneral Contact（2004年〜全体接触定義）、ANSYSのAugmented Lagrangian接触、NX NastranのSOL 101接触は設定の手間と精度のバランスが異なる。2019年のSolverBench比較試験（ヘルツ接触問題）でAbaqus General ContactがANSY比で計算時間27%短縮、最大圧力誤差は同等の1.5%以内と評価された。
選定で最も重要な3つの問い
- 「何を解くか」：ギャップ要素に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
ギャップ要素の先端トピック
ギャップ要素に先端研究はありますか？
ギャップ要素自体は古典的な要素だが、接触力学の進化とともに発展がある。
モルタル法による接触の改善
ギャップ要素のペナルティ法は閉合剛性の設定が難しい。モルタル法（Mortar method）は弱形式で接触条件を課すことで、ペナルティパラメータに依存しない安定した接触解を実現する。
モルタル法はギャップ要素に適用できますか？
モルタル法は面接触定義で主に使われるが、1次元の「点接触」にもLagrange multiplier法として適用可能。AbaqusやAnsysの接触アルゴリズムはモルタル法ベースに移行しつつある。
確率論的ギャップ解析
製造公差によるギャップ量のばらつきを確率変数として扱い、構造応答の分布を評価する研究がある。モンテカルロシミュレーションでギャップ量を変化させ、接触の有無が構造の剛性・強度にどう影響するかを統計的に評価する。
まとめ
ギャップ要素の先端研究、まとめます。
ギャップ要素は「古典的」だが、接触力学という広い分野の入口として重要だ。
摩擦付きギャップ要素の定式化
クーロン摩擦を組み込んだギャップ要素はタンジェント剛性マトリクスが非対称となり、通常の対称ソルバーは使用不可だ。1980年にSimoとLaurentは非対称スパースソルバーを用いた摩擦接触アルゴリズムを発表。この成果はAbaqus v4.5（1984年）に実装され、航空機のリベット接合部解析に初めて実用応用された。
トラブルシューティング
ギャップ要素のトラブル
ギャップ要素でよくあるトラブルを教えてください。
非線形要素ならではのトラブルがある。
収束しない（チャタリング）
反復がいつまでも収束しません。
ギャップの開閉が反復ごとに切り替わるチャタリング。ギャップが「閉→開→閉→開」を繰り返す。
対策：
- 閉合剛性を下げる（もっと柔らかいばねにする）
- 増分荷重のステップを小さくする
- 減衰を追加（NastranのPARAM,ADPCON）
- 安定化法を使用
貫通量が大きい
接触しているのに2節点が重なっています。
閉合剛性が小さすぎる。$k$ を上げれば貫通量は減るが、上げすぎると収束困難。バランスが大事。
ギャップの方向が間違っている
接触すべきなのに接触しません。
ギャップの方向ベクトルが間違っている。NastranのCGAPは方向ベクトル（GA-GB方向）で接触方向を定義する。方向が逆だと「押し離す」方向に力が作用してしまう。
まとめ
ギャップ要素のトラブル対処、整理します。
ギャップ要素は「閉合剛性 $k$ の設定」が全てですね。
そう。$k$ が適切なら問題なく動く。$k$ が不適切なら収束しない。シンプルだが奥が深い。
接触チャタリングの対処法
ギャップ要素解析での「チャタリング」は接触状態が収束前にON/OFFを繰り返す現象で、解析が発散する原因になる。対処法として、Abaqusでは*CONTACT CONTROLのSTABILIZATION=0.001設定が2005年以降の標準対策だ。NastranのGAP要素では初期ギャップ量を実測より5〜10%小さく設定することで収束性が改善すると現場で経験則が知られる。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——ギャップ要素の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告