RBE3加重平均要素
理論と物理
RBE3とは
先生、RBE3はRBE2とどう違うんですか？
RBE3は剛性を追加しない荷重分配要素だ。RBE2が「剛体結合」なのに対し、RBE3は「加重平均」。この違いがFEMモデル化で最も重要な区別の一つだ。
動作原理
RBE3の動作：
- 参照点（reference point）の変位 = 周囲の独立節点の変位の加重平均
- 参照点に作用する力 = 周囲の独立節点に重み付きで分配
数学的に：
マスターとスレーブの関係がRBE2と逆ですか？
RBE3では独立節点が周囲のノード（荷重を受ける側）で、従属節点が参照点（荷重を与える点）。RBE2とは独立/従属が逆だ。
重要な違い：
| 特性 | RBE2 | RBE3 |
|---|---|---|
| 剛性の追加 | あり（無限大） | なし |
| 独立節点 | マスター（1点） | 周囲ノード（多点） |
| 従属節点 | スレーブ（多点） | 参照点（1点） |
| 物理的イメージ | 溶接接合 | 吊り荷重の分配 |
「吊り荷重の分配」がイメージしやすいです。1点で吊った荷物の重さがロープを通じて複数の支持点に分配される。
完璧なイメージだ。RBE3は「柔らかいロープで吊っている」ようなもの。力は分配されるが、支持構造の剛性は変わらない。
なぜRBE3が重要か
なぜRBE3がRBE2より推奨されることが多いんですか？
実構造の接合部は完全な剛体ではないからだ。ボルト接合やピン接合は有限の剛性を持つ。RBE2で結合すると接合部が無限に硬くなり、結果が非現実的になる。RBE3は剛性を変えずに力だけ伝達するから、実構造に近い。
例：クレーン荷重をフランジに伝達する場合
- RBE2 → フランジが剛体化。フランジの変形が出ない。周囲に応力集中
- RBE3 → フランジの変形はそのまま。荷重が分配されるだけ。現実的
重み係数
RBE3の「重み」はどう設定しますか？
重み $w_i$ は力の分配比率を決める。全て同じ重み（$w_i = 1$）なら均等分配。節点によって重みを変えれば不均等分配も可能。
実務的には全て $w_i = 1$（均等分配）が最も一般的。荷重分布が不均一な場合は、節点の支配面積に比例した重みを設定する。
まとめ
RBE3の理論を整理します。
要点：
- 荷重分配要素 — 剛性を追加しない。RBE2との根本的な違い
- 参照点の変位 = 周囲の加重平均 — 力は重み付きで分配
- 独立/従属がRBE2と逆 — 周囲ノードが独立、参照点が従属
- 荷重の分配にはRBE3を使うべき — RBE2は剛性過大の原因
- 重み $w_i = 1$（均等分配）が標準 — 必要に応じて不均等も可
「荷重分配にはRBE3」。これがFEMモデル化の鉄則ですね。
そう。RBE2 vs. RBE3の選択はFEMで最も間違えやすく、最も影響が大きい設定だ。この違いを理解していないエンジニアは、FEMの結果を信用すべきでない。
RBE3の加重平均定式
RBE3（Rigid Body Element 3）はNastranに1970年代に追加された荷重分配要素で、独立節点群の変位加重平均を参照節点の変位に等しくする拘束を設ける。RBE2と異なり剛性を追加しないため「ゼロ剛性要素」とも呼ばれる。重み係数Wiは面積・長さ・一定値などで指定でき、非均一荷重分布の模擬に威力を発揮する。
各項の物理的意味
- 慣性項（質量項）：$\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で統一 |
数値解法と実装
RBE3の各ソルバーでの実装
RBE3は各ソルバーでどう設定しますか？
Nastran
```
RBE3, 200, , 5000, 123456, 1.0, 123, 2001, 2002,+
+, 2003, 2004
```
- 200: 要素ID
- 5000: 参照点（従属節点）
- 123456: 参照点の拘束DOF
- 1.0: 重み
- 123: 独立節点のDOF
- 2001〜2004: 独立節点
Abaqus
```
*COUPLING, CONSTRAINT NAME=rbe3_1, REF NODE=5000
*DISTRIBUTING
slave_surface, 1, 6
```
Ansys
```
RBE3, 5000, , 2001, 1.0, UX, UY, UZ
```
Abaqusでは*COUPLINGのDISTRIBUTINGオプションがRBE3に相当するんですね。
そう。Abaqusの記法は「KINEMATIC = RBE2、DISTRIBUTING = RBE3」と覚えるとよい。
RBE3の注意点
RBE3を使うときの注意点は？
RBE3の参照点（従属節点）は別途拘束しない限り自由に動く。RBE3だけでは構造を支持できない。
例えば、RBE3の参照点に荷重を与え、周囲の独立節点は全て自由（SPCなし）だと、構造全体が剛体移動する。RBE3は力を分配するだけで、構造の安定性は保証しない。
支持はRBE2（または直接SPC）で、荷重分配はRBE3で、という使い分けですね。
まさにそう。支持 = RBE2/SPC、荷重 = RBE3。この組み合わせが実務の基本パターンだ。
RBE3のDOF設定
独立節点のDOFは何を指定すべきですか？
通常は123（並進3方向）を指定。回転DOF（456）は独立節点がシェルや梁で回転DOFを持つ場合に指定する。ソリッド要素の節点は回転DOFを持たないから、123のみ。
ソリッド要素に456を指定したらどうなりますか？
ソルバーによっては警告なしに無視されるか、エラーになる。ソリッド要素のノードに対しては123のみを指定すること。
まとめ
RBE3の実装詳細、整理します。
要点：
- Nastran: RBE3カード、Abaqus: *COUPLING DISTRIBUTING — 記法は異なるが同じ機能
- 参照点は別途拘束が必要 — RBE3だけでは構造が不安定
- 独立節点のDOF — ソリッドは123のみ、シェル/梁は123456
- 重み $w = 1$（均等）が標準 — 不均等分配も可能
RBE3の荷重分配アルゴリズム
RBE3の参照節点に集中荷重Fを加えると、従属節点iへの分配荷重はFi = (Wi × Ai / ΣWj×Aj) × Fで計算される。Aiは各節点の寄与面積で、均一重みWi=1の場合は単純な節点数分割になる。ヘキサコア方程式ソルバーを用いるとRBE3処理時間が従来LAPACKベース比で約1/3に短縮されることがSiemens NX Nastran 2021リリースノートに記載されている。
線形要素（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要素あたりの計算コストは増えるため、トータルのコスト対効果で判断する。
実践ガイド
RBE3の実務適用
RBE3の典型的な使い方を教えてください。
1. 集中荷重の面分配
1つの点に作用する力（クレーン荷重、装置荷重等）をフランジやブラケットの面に分配する。参照点に力を与え、RBE3で面の節点に分配。
2. ベアリング反力の分配
ベアリングの反力を軸受け穴の周囲の節点に分配。1点集中ではなく穴の周囲に面分配される。
3. 梁要素の結果をシェルに伝達
全体モデル（梁要素）の結果をサブモデル（シェル/ソリッド）に伝達する際、梁の端点の力をRBE3でサブモデルの面に分配する。
4. 振動解析の質量分配
集中質量（CONM2）をRBE3経由で面に接続。質量の力が面に分配されるが、構造の剛性は変わらない。
RBE2 vs. RBE3 の判断フロー
RBE2とRBE3のどちらを使うか迷ったときの判断フローを教えてください。
1. 「この接続部は物理的に剛体か？」 → Yes: RBE2 / No: RBE3
2. 「荷重の分配が目的か？」 → Yes: RBE3
3. 「面が平面を保つ必要があるか？」 → Yes: RBE2
4. 「接続部の変形を見たいか？」 → Yes: RBE3（剛性を追加しない）
5. 迷ったら → RBE3（安全側。剛性を追加しないほうが保守的）
「迷ったらRBE3」がデフォルト。RBE2は「本当に剛体が必要な場合のみ」。
この判断フローを頭に入れておくだけで、FEMモデル化の品質が劇的に向上する。
実務チェックリスト
RBE3のチェックリストをお願いします。
RBE3は「荷重分配に使い、支持には使わない」。シンプルだけど大事な原則ですね。
RBE2とRBE3の使い分けを完璧に理解することは、FEMエンジニアの最も重要な基本スキルだ。これを間違えると全ての解析結果が信用できなくなる。
RBE3による分布荷重のモデル化
風洞試験で得た翼面圧力分布をFEMに入力する際、RBE3を使って試験点群から有限要素節点への荷重分配を行う手法が航空宇宙業界で確立されている。ボーイング787の翼構造解析でもこの手法が採用されており、CFD圧力データ（約50万点）をFEM節点（約3万点）に変換する際にRBE3の面積比重みが精度の鍵となっている。
解析フローのたとえ
解析の流れは、実は料理とそっくりです。まず材料を買い出し（CADモデルの準備）、下ごしらえをして（メッシュ生成）、火にかけて（ソルバー実行）、最後に盛り付ける（後処理で可視化）。ここで大事な問いかけ——料理で一番失敗しやすい工程はどこでしょう？ 実は「下ごしらえ」なんです。メッシュの品質が悪いと、どんなに優秀なソルバーを使っても結果はめちゃくちゃになります。
初心者が陥りやすい落とし穴
あなたはメッシュ収束性を確認していますか？ 「計算が回った＝結果が正しい」と思っていませんか？ これ、実はCAE初心者が最も陥りやすい罠です。ソルバーは与えられたメッシュで「それなりの答え」を必ず返します。でもメッシュが粗すぎれば、その答えは現実から大きくずれている。最低3段階のメッシュ密度で結果が安定することを確認する——これを怠ると「コンピュータが出した答えだから正しいはず」という危険な思い込みに陥ります。
境界条件の考え方
境界条件の設定は、試験の「問題文を書く」のと同じです。問題文が間違っていたら？ どんなに正確に計算しても答えは間違いますよね。「この面は本当に完全固定なのか」「この荷重は本当に一様分布なのか」——現実の拘束条件を正しくモデル化することが、実は解析全体で最も重要なステップだったりします。
ソフトウェア比較
RBE3のソルバー比較
RBE3の各ソルバーでの実装差は？
Abaqusの*DISTRIBUTING COUPLINGはNastranのRBE3と完全に同じですか？
概念は同じだが、Abaqusの実装は面ベースで均等分配が基本。NastranのRBE3のほうが節点ごとの重みの細かい制御が得意。ただし実用上の差は小さい。
選定ガイド
RBE3はどのソルバーでも基本機能。ツールの差よりもエンジニアの判断力が大事。
その通り。RBE3 の「使い方の判断」は全ソルバー共通の知識であり、FEMの基本教養だ。
各ソルバーのRBE3相当要素
NastranのRBE3に相当する機能として、AbaqusのDistributing Coupling（*COUPLING, DISTRIBUTING）、AnsysのForce Distributed Constraint（RBE3コマンド）、Altair OptiStructのRBE3がある。SIMULIA 2023バージョンからはDistributing CouplingにPROJECTIONオプションが追加され、曲面上の節点群への荷重分配がより正確になった。実務上RBE3はRBE2の「柔らかい代替品」として応力集中を回避したい箇所に多用される。
選定で最も重要な3つの問い
- 「何を解くか」：RBE3加重平均要素に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
RBE3の先端トピック
RBE3に先端研究はありますか？
RBE3自体はシンプルだが、モデル間の荷重伝達とマルチスケール解析で重要な役割を果たす。
サブモデリングでのRBE3
グローバルモデルからサブモデルへの荷重伝達にRBE3を使う。グローバルモデルの反力をRBE3でサブモデルの面に分配する。変位境界条件よりも力の分配が自然な場合がある。
RBE3の動的特性
動的解析でRBE3を使う場合、参照点に集中質量（CONM2）を配置することが多い。RBE3は剛性を追加しないので、質量だけが構造に付加される。固有振動モードに影響する質量を正しく分配できる。
RBE2で質量を分配したら構造の剛性が変わってしまいますもんね。
まさにそう。「質量を追加したいが剛性は変えたくない」→ CONM2 + RBE3 の組み合わせが最適解。
加重平均の高度化
通常のRBE3は変位の加重平均だが、応力の加重平均や温度の加重平均に拡張する研究もある。マルチフィジックス連成での荷重伝達に応用される。
まとめ
RBE3の先端トピック、まとめます。
RBE3は「荷重の伝達」という基本的な役割で、FEMのモデル化を支えている。
RBE3の特異解と対策
RBE3は剛性を付与しないため、従属DOFが参照節点のDOFに対して一次独立でない場合に特異解が生じる。1990年代にNASAのTurbullらが「回転DOFのみ拘束したRBE3に並進荷重を加える」と特異になる問題を報告し、独立節点群がコプレーナー（同一平面）配置の場合の最低6節点配置ルールを提案した。Nastranでは診断メッセージ番号5293として警告が出る。
トラブルシューティング
RBE3のトラブル
RBE3でよくあるトラブルを教えてください。
RBE3はRBE2より安全な要素だが、特有のトラブルがある。
構造が不安定（特異剛性マトリクス）
RBE3を使ったら「特異剛性マトリクス」のエラーが出ました。
RBE3は剛性を追加しない。RBE3の参照点だけでは構造を支持できない。
対策：
- 構造の支持はSPC（拘束）またはRBE2で行う
- RBE3は荷重分配のみに使う
- RBE3の参照点に荷重を与える場合、構造の他の部分で支持されていることを確認
荷重分配が不均一
RBE3で均等に分配したつもりが、力の分配が偏っています。
RBE3の分配は重みと距離に依存する。全ての独立節点の重みが同じでも、参照点からの距離が異なれば力の分配は不均一になる。
対策：
- 重みを面積比例に設定（各節点の支配面積に比例）
- 独立節点の配置を参照点から概ね等距離にする
- 分配された力をソルバーの出力で確認
RBE2との混同
RBE3を使うべきところでRBE2を使ってしまいました。
FEMで最もよくあるミス。以下の症状で判別：
| 症状 | RBE2で起きる | RBE3で起きない |
|---|---|---|
| 接続部が硬すぎる | ○ | — |
| たわみが過小 | ○ | — |
| 接続部に応力集中 | ○ | — |
| 全体剛性が変わる | ○ | — |
「たわみが過小」「接続部の応力集中」があったらRBE2→RBE3への切り替えを検討すべきですね。
その通り。全ての解析者が覚えるべきルール：荷重分配にはRBE3、剛体結合にはRBE2。
まとめ
RBE3のトラブル対処、整理します。
RBE2 vs. RBE3の正しい使い分け…これがFEMの最も重要なスキルだと実感しました。
FEMの技術は要素理論、メッシュ、ソルバーと多岐にわたるが、実務で最もインパクトが大きいのは「接続のモデル化」だ。RBE2/RBE3の使い分けを制するものがFEMを制する。
RBE3の回転DOF欠落エラー
RBE3でUM（従属DOF）に回転自由度を含めない場合、参照節点のモーメントが正しく分配されず、解析結果が非物理的な変位分布を示す。この問題はNastranのWarning 5291として出力される。対策は従属節点群のUMリストに123456（全DOF）を明示的に記載するか、RBE3を2段重ねにしてモーメントアームを確保することだ。Abaqus Distributingカップリングでは自動的に処理される。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——RBE3加重平均要素の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告