剛体要素
剛体要素の理論基礎
剛体要素とは
先生、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倍の弾性率を模擬していると解釈できる。
剛体要素の数値計算手法
MPC（多点拘束）の実装
剛体要素のMPCはどうやって実装されるんですか？
2つの方法がある：
1. 従属DOFの消去（変換法）
従属節点のDOFを独立節点のDOFで表現し、全体方程式から消去する。連立方程式のサイズが小さくなる。NastranのRBE2はこの方法。
2. ペナルティ法
剛体拘束を「非常に硬いばね」で近似する。完全な剛体ではないが、実用上十分な精度。ペナルティ値が大きすぎると条件数が悪化。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を消去することでシステム方程式のサイズを維持するまま計算できる。
剛体要素の実務適用
剛体要素の実務適用
剛体要素の実務での使い方を教えてください。
最も一般的な用途を整理しよう。
ボルト接合のモデル化
ボルト接合部を簡略的に表現する場合：
RBE2方式 — ボルト穴の周囲の節点をRBE2で結合。「ボルト穴の部分が剛体」として扱う。簡単だが接合部が硬くなりすぎる。
RBE3＋ばね方式 — ボルト位置にRBE3で荷重を分配し、ばね要素でボルトの軸剛性を表現。より現実的だが設定が複雑。
荷重の分配
「1つの点に荷重を与えて、面に分配する」場面でRBE3が活躍。クレーン荷重をフランジに分配、ジャッキ荷重を構造に分配等。
境界条件の表現
「面の平均変位をゼロにする」「面が平面を保つ」等の境界条件をRBE2/RBE3で表現。圧力容器の端面条件（面が平面を保つ＋軸方向移動は自由）はRBE2で設定する。
実務チェックリスト
剛体要素のチェックリストをお願いします。
「RBE2 vs. RBE3の判断」が最初のチェック項目。これが全てを決めるんですね。
RBE2を使うべき場所にRBE3を使うと剛性不足、RBE3を使うべき場所にRBE2を使うと剛性過大。この判断を間違えると結果が根本的に変わる。
剛体要素によるエンジンマウント解析
自動車エンジンマウント解析では、エンジンブロックを一体の剛体要素（RBE2またはRBARの集合）で模擬し計算コストを大幅削減する手法が用いられる。ホンダの公開技報（2015年）によれば、エンジンを詳細FEMから剛体置換することでモデル規模を約60%削減しながらマウント部の反力精度を5%以内に維持できると報告されている。
剛体要素のソフトウェア比較
剛体要素のツール比較
各ソルバーの剛体要素を比較してください。
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では接触ベースの拘束が剛体要素の代替として多用される。
剛体要素の先端研究
剛体要素の先端トピック
剛体要素に先端研究はありますか？
剛体要素自体はシンプルだが、サブストラクチャリングとマルチスケール解析で重要な役割を果たす。
スーパーエレメントと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出力で局所的な変位ジャンプとして現れる。
なった
詳しく
報告