剛体要素
理論と物理
剛体要素とは
先生、FEMの「剛体要素」って何ですか？ 変形しない要素？
そう。剛体要素は変形しない。2つ以上の節点を「剛体的に結合」して、一方の変位・回転が他方に伝わるようにする要素だ。
物理的にはどんな構造を表現するんですか？
数学的な定式化
剛体要素はFEMの多点拘束（MPC: Multi-Point Constraint）として実装される。独立節点（マスターノード）の変位 $\{u_m\}$ に対して、従属節点（スレーブノード）の変位が：
$[T]$ は変位変換マトリクスで、剛体運動の式に基づく。
従属節点は独立節点の変位から「自動的に決まる」んですね。
そう。従属節点の自由度は消去される（condensation）。実質的にDOFが減少する。
剛体要素の種類
| 種類 | Nastran | Abaqus | 特徴 |
|---|---|---|---|
| 剛体結合（全DOF） | RBE2 | RIGID BODY / COUPLING | 全自由度を剛体的に結合 |
| 加重平均（荷重分配） | RBE3 | *DISTRIBUTING COUPLING | 力を分配。剛性を追加しない |
| 剛体面 | RBAR | *RIGID BODY | 2節点の剛体リンク |
RBE2とRBE3は全然違うんですか？
本質的に異なる。RBE2は「剛体結合」で構造の剛性を追加する。RBE3は「荷重分配」で剛性を追加しない。この違いを理解しないと重大なモデル化ミスになる。詳しくはRBE2、RBE3のページで解説する。
剛体要素の注意点
剛体要素は便利だが、使い方を間違えると結果を大きく歪める：
1. 剛体要素で局所的な剛性が過大に — 接続部が不自然に硬くなる
2. 応力特異点を生む — 剛体要素と変形要素の境界で応力が不正確
3. 過拘束 — 不要な自由度まで拘束してしまう
剛体要素の境界で応力が信用できないんですか。
剛体要素に接続する変形要素の最初の1〜2層の応力は参考程度。離れた位置の応力を評価すること。
まとめ
剛体要素の理論を整理します。
要点：
- 変形しない要素 — 節点間を剛体的に結合
- MPCとして実装 — 従属節点のDOFが消去される
- RBE2（剛体結合）とRBE3（荷重分配）は全く異なる — 混同厳禁
- 局所的な剛性過大に注意 — 接続部の応力は不正確
- 異なる要素間の接続に便利 — 梁-ソリッド、シェル-ソリッド
RBE2 vs. RBE3 の理解がFEMモデル化の最も重要なスキルの一つですね。
まさにそう。この2つを正しく使い分けられるかどうかで、FEMエンジニアの力量がわかる。
剛体要素の力学的基礎
剛体要素は内部変形を持たない理想剛体として振る舞う有限要素で、1960年代にターナー（M.J. Turner）らがNASAでの宇宙構造解析向けに導入した。数学的には剛体変換行列Tにより独立節点から従属節点の全6自由度を決定論的に計算する。実際の材料は弾性変形するため、剛体要素は真の剛性に比べ少なくとも10〜100倍の弾性率を模擬していると解釈できる。
各項の物理的意味
- 慣性項（質量項）：$\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で統一 |
数値解法と実装
MPC（多点拘束）の実装
剛体要素のMPCはどうやって実装されるんですか？
2つの方法がある：
1. 従属DOFの消去（変換法）
従属節点のDOFを独立節点のDOFで表現し、全体方程式から消去する。連立方程式のサイズが小さくなる。NastranのRBE2はこの方法。
2. ペナルティ法
剛体拘束を「非常に硬いばね」で近似する。完全な剛体ではないが、実用上十分な精度。ペナルティ値が大きすぎると条件数が悪化。AbaqusのKINEMATICオプションは変換法、PENALTYオプションはペナルティ法。
どちらが良いですか？
変換法（KINEMATIC）が基本推奨。ペナルティ法は接触と組み合わせる場合に使うことがある（接触のペナルティ法と整合させるため）。
ソルバー別の実装
剛体要素のMPCはどうやって実装されるんですか？
2つの方法がある：
従属節点のDOFを独立節点のDOFで表現し、全体方程式から消去する。連立方程式のサイズが小さくなる。NastranのRBE2はこの方法。
剛体拘束を「非常に硬いばね」で近似する。完全な剛体ではないが、実用上十分な精度。ペナルティ値が大きすぎると条件数が悪化。AbaqusのKINEMATICオプションは変換法、PENALTYオプションはペナルティ法。
どちらが良いですか？
変換法（KINEMATIC）が基本推奨。ペナルティ法は接触と組み合わせる場合に使うことがある（接触のペナルティ法と整合させるため）。
| 機能 | Nastran | Abaqus | Ansys |
|---|---|---|---|
| 剛体結合 | RBE2 | RIGID BODY / COUPLING, KINEMATIC | CERIG |
| 荷重分配 | RBE3 | *COUPLING, DISTRIBUTING | RBE3(MPC) |
| 剛体リンク（2節点） | RBAR | *RIGID BODY | MPC184 |
| 剛体面 | RBE2（多節点） | *RIGID BODY | TARGE170 |
剛体要素の接続パターン
典型的な接続パターン：
梁要素とソリッド要素の接続
梁要素は回転DOFを持つが、ソリッド要素は持たない。RBE2/RBE3で接続：
- RBE2 — 梁の端点をマスター、ソリッドの面の節点をスレーブ。断面が剛体的に動く（硬すぎる場合がある）
- RBE3 — 梁の端点をスレーブ（参照点）、ソリッドの面の節点をマスター（荷重分配点）。剛性を追加しない（より現実的）
RBE3のほうが「優しい」接続ですね。
その通り。RBE2は「固い溶接」のような接続、RBE3は「バランスの取れた荷重伝達」のような接続。実構造のボルト接合はRBE3に近いことが多い。
まとめ
剛体要素の数値手法、整理します。
要点：
- 変換法（KINEMATIC）が基本推奨 — DOF消去で正確
- ペナルティ法は接触との組み合わせ用 — 近似的な剛体
- 梁-ソリッドの接続 — RBE2（硬い）vs. RBE3（柔らかい）
- 接続部の応力は不正確 — 1〜2要素離れた位置で評価
剛体要素の実装方式比較
剛体要素の実装には主に①ペナルティ法、②ラグランジュ乗数法、③ダイレクト剛性縮合の3方式がある。ペナルティ法は実装が簡単だが罰則係数の選択が解の精度に影響する。ラグランジュ乗数法は数学的に厳密だが連立方程式のサイズが増大する。Nastranはダイレクト縮合法を採用しており、従属DOFを消去することでシステム方程式のサイズを維持するまま計算できる。
線形要素（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要素あたりの計算コストは増えるため、トータルのコスト対効果で判断する。
実践ガイド
剛体要素の実務適用
剛体要素の実務での使い方を教えてください。
最も一般的な用途を整理しよう。
ボルト接合のモデル化
ボルト接合部を簡略的に表現する場合：
RBE2方式 — ボルト穴の周囲の節点をRBE2で結合。「ボルト穴の部分が剛体」として扱う。簡単だが接合部が硬くなりすぎる。
RBE3＋ばね方式 — ボルト位置にRBE3で荷重を分配し、ばね要素でボルトの軸剛性を表現。より現実的だが設定が複雑。
荷重の分配
「1つの点に荷重を与えて、面に分配する」場面でRBE3が活躍。クレーン荷重をフランジに分配、ジャッキ荷重を構造に分配等。
境界条件の表現
「面の平均変位をゼロにする」「面が平面を保つ」等の境界条件をRBE2/RBE3で表現。圧力容器の端面条件（面が平面を保つ＋軸方向移動は自由）はRBE2で設定する。
実務チェックリスト
剛体要素のチェックリストをお願いします。
「RBE2 vs. RBE3の判断」が最初のチェック項目。これが全てを決めるんですね。
RBE2を使うべき場所にRBE3を使うと剛性不足、RBE3を使うべき場所にRBE2を使うと剛性過大。この判断を間違えると結果が根本的に変わる。
剛体要素によるエンジンマウント解析
自動車エンジンマウント解析では、エンジンブロックを一体の剛体要素（RBE2またはRBARの集合）で模擬し計算コストを大幅削減する手法が用いられる。ホンダの公開技報（2015年）によれば、エンジンを詳細FEMから剛体置換することでモデル規模を約60%削減しながらマウント部の反力精度を5%以内に維持できると報告されている。
解析フローのたとえ
解析の流れは、実は料理とそっくりです。まず材料を買い出し（CADモデルの準備）、下ごしらえをして（メッシュ生成）、火にかけて（ソルバー実行）、最後に盛り付ける（後処理で可視化）。ここで大事な問いかけ——料理で一番失敗しやすい工程はどこでしょう？ 実は「下ごしらえ」なんです。メッシュの品質が悪いと、どんなに優秀なソルバーを使っても結果はめちゃくちゃになります。
初心者が陥りやすい落とし穴
あなたはメッシュ収束性を確認していますか？ 「計算が回った＝結果が正しい」と思っていませんか？ これ、実はCAE初心者が最も陥りやすい罠です。ソルバーは与えられたメッシュで「それなりの答え」を必ず返します。でもメッシュが粗すぎれば、その答えは現実から大きくずれている。最低3段階のメッシュ密度で結果が安定することを確認する——これを怠ると「コンピュータが出した答えだから正しいはず」という危険な思い込みに陥ります。
境界条件の考え方
境界条件の設定は、試験の「問題文を書く」のと同じです。問題文が間違っていたら？ どんなに正確に計算しても答えは間違いますよね。「この面は本当に完全固定なのか」「この荷重は本当に一様分布なのか」——現実の拘束条件を正しくモデル化することが、実は解析全体で最も重要なステップだったりします。
ソフトウェア比較
剛体要素のツール比較
各ソルバーの剛体要素を比較してください。
NastranのRBE2/RBE3が最も実績がありますか？
RBE2/RBE3はNastranの発明と言ってよい。1970年代からの歴史があり、航空宇宙の大規模モデルで膨大な使用実績がある。AbaqusやAnsysの同等機能も本質的には同じだが、名前と設定方法が異なる。
選定ガイド
「名前が違うだけ」とはいえ、設定方法が結構違いますよね。
NastranのRBE2/RBE3カードの記法と、Abaqusの*COUPLINGの記法は全然違う。それぞれのマニュアルをよく読んで、同じ物理を異なる記法で設定できるようになることが重要だ。
剛体要素の各ソルバー実装一覧
剛体要素の主要実装はNastranのRBAR/RBE1/RBE2/RBE3/RROD/RTRPLT、AbaqusのRigid Body/Kinematic Coupling、AnsysのCE/CP/RBE3コマンド群に分類される。Siemens NX Nastran 2022ではAUTORBEという自動剛体要素生成機能が追加され、CADソリッドから直接剛体サブ構造を抽出できるようになった。LS-DYNAでは接触ベースの拘束が剛体要素の代替として多用される。
選定で最も重要な3つの問い
- 「何を解くか」：剛体要素に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
剛体要素の先端トピック
剛体要素に先端研究はありますか？
剛体要素自体はシンプルだが、サブストラクチャリングとマルチスケール解析で重要な役割を果たす。
スーパーエレメントとRBE2/RBE3
スーパーエレメント（部分構造法）では、構造の一部を縮約して「1つの大きな要素」にする。縮約された部分のインターフェースはRBE2/RBE3で接続される。航空宇宙の大規模モデル（数千万DOF）で必須の技術。
デジタルモック
自動車のデジタルモック（全車モデル）では数千のRBE2/RBE3が使われる。ボディ、シャーシ、エンジン、サスペンションの接続を全て剛体要素で表現。接続のモデル化が全車振動特性を支配する。
接触との統合
Abaqusの*RIGID BODYは接触の剛体面（ダイ、パンチ、工具）として使える。プレス成形ではダイを剛体面、ブランクを変形体としてモデル化する。
まとめ
剛体要素の先端トピック、まとめます。
剛体要素は「見えない裏方」だが、大規模解析のインフラストラクチャとして不可欠だ。
剛体要素と柔軟体の連成モデル
マルチボディダイナミクス（MBD）とFEMを連成させる際、剛体要素は柔軟体（CMS：コンポーネントモード合成）の慣性特性を表現するリファレンスノードとして機能する。1990年代にSimPackやADAMS/Flexとの連成解析が普及し、柔軟ロボットアームの先端振動解析では剛体要素ベースのCraig-Bampton縮合が業界標準となった。
トラブルシューティング
剛体要素のトラブル
剛体要素でよくあるトラブルを教えてください。
剛体要素のトラブルは意図しない剛性の追加とDOFの過拘束が二大原因だ。
RBE2で局所剛性が過大
RBE2を使ったら接続部の応力が異常に高くなりました。
RBE2は接続部を完全に剛体化する。本来は変形する部分が剛体になるため、隣接する変形要素に応力が集中する。
対策：
- RBE3に切り替え — 荷重分配のみで剛性を追加しない
- 接続部の応力は評価しない — 1〜2要素離れた位置で評価
- ばね要素で置き換え — 有限の剛性で接続
RBE3で構造が不安定
RBE3を使ったら「特異剛性マトリクス」のエラーが出ました。
RBE3は剛性を追加しない。RBE3だけで構造を支持しようとすると、構造が不安定（剛体移動が拘束されない）になる。
対策：
- RBE3の参照点（スレーブ）に別途拘束を追加
- RBE3は「荷重の分配」に使い、「構造の支持」にはRBE2を使う
DOFの過拘束
剛体要素が多すぎて過拘束になります。
RBE2で全6自由度を拘束すると、接続部の変形が完全に抑制される。不要な自由度（例：回転は自由にしたい）がある場合は、拘束するDOFを指定する。
NastranのRBE2ではCM（拘束する成分）フィールドで123456の一部を指定。全DOFではなく必要なDOFだけ拘束する。
まとめ
剛体要素のトラブル対処、整理します。
この根本原則を忘れなければ、ほとんどのトラブルは防げそうですね。
その通り。RBE2とRBE3の本質的な違いを理解することが、FEMモデル化の最も重要なスキルだ。
剛体要素の過拘束チェック方法
複数の剛体要素が同一節点に異なるDOFを拘束する「過拘束」は解析の致命的エラーとなる。MSC NastranではMPC/RBE要素が絡むと Fatal 9050が頻出する。診断手順は①AUTOSPC出力でSPC力を確認、②PARAM,CHECKOUT,YESでMPC拘束マトリクスを書き出し、③重複DOF番号を検索する3ステップが有効だ。Abaqusでは*NODE PRINTのU出力で局所的な変位ジャンプとして現れる。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——剛体要素の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
なった
詳しく
報告