剛体動力学
剛体動力学の理論基礎
剛体動力学
先生、剛体動力学とFEMの動的解析はどう違いますか？
FEMの動的解析は構造の弾性変形を扱う。剛体動力学は変形しない物体の運動（並進＋回転）を追跡。機構（メカニズム）の運動シミュレーションが目的。
運動方程式
ニュートン-オイラー方程式：
$m$: 質量、$[I]$: 慣性テンソル、$\mathbf{F}$: 力、$\mathbf{M}$: モーメント。
拘束条件（ジョイント）
ジョイントで剛体間の相対運動を拘束：
| ジョイント | 自由度 | 例 |
|---|---|---|
| 固定（Weld） | 0 | 溶接接合 |
| 回転（Revolute） | 1（回転） | ヒンジ、軸受 |
| 並進（Prismatic） | 1（並進） | スライダー |
| 円筒（Cylindrical） | 2 | ピストン |
| 球（Spherical） | 3（回転3） | ボールジョイント |
| 自由 | 6 | 拘束なし |
まとめ
オイラーの剛体方程式は1758年に提案された
剛体の回転運動を記述する「オイラー方程式」はレオンハルト・オイラーが1758年に発表した『Novi Commentarii academiae scientiarum Petropolitanae』に掲載された。慣性テンソルと角速度ベクトルの関係を主慣性軸系で簡潔に表現したこの式は、現代のマルチボディ剛体ダイナミクスの数値積分の核心だ。コマやジャイロスコープの歳差運動もオイラー方程式から直接導出でき、宇宙機の姿勢安定解析の基礎として今も教科書第1章に登場する。
剛体動力学の数値計算手法
MBDのソルバー
まとめ
Verlet積分はゲームエンジンと分子動力学を変えた
剛体MBD数値積分のうち「Störmer-Verlet法」（1907年のStormerの論文、1967年にVerletが分子動力学に再適用）はエネルギー保存性が高いシンプレクティック積分の代表格だ。計算コストが低く長時間積分で誤差が蓄積しにくいため、VelocityVerlet変形がUnity・Unreal Engineの物理エンジンに採用されている。CAEの剛体MBDでは接触なしの保守系問題では今も有効だが、非保守力（減衰）が大きい問題では高次Runge-Kuttaの方が精度面で優れる。
剛体動力学の実務適用
剛体動力学の実務
自動車のサスペンション運動学、ロボットの関節運動、エンジンのクランク機構、展開構造の運動。
実務チェックリスト
輸送包装の落下試験解析は剛体MBDが最速
電子機器・家電の輸送包装設計における落下試験解析（1.2m高さからの自由落下）は、有限要素法では計算コストが高いため剛体MBDが実務の主流だ。Apple iPhone（iPhone 12以降）のパッケージ設計では、剛体MBD＋ポリウレタンフォームの非線形バネモデルで衝撃加速度を300Gから120G以下に収める設計を行い、Apple Park内部テスト施設での落下試験との誤差を±15%以内に抑えていることをAppleエンジニアの学会発表が示唆している。
剛体動力学のソフトウェア比較
ツール
Modelicaベースの剛体MBDは教育・研究機関で急拡大
オープン標準言語Modelicaに基づくオープンソースMBD環境OpenModelica（リンショーピング大学、1997年〜）は研究・教育機関で急速に普及している。商用では Dymola（Dassault Systèmes）がModelicaの最大手実装であり、ボルボ社のハイブリッドパワートレイン制御とMBD剛体ダイナミクスを単一モデルで連成するシステムシミュレーションに使っていることが2020年MathWorks技術交流会のVolvo発表で紹介された。
剛体動力学の先端研究
剛体動力学の先端
ジャイロスコープ効果を利用した船舶安定化装置
大型クルーズ船の横揺れ（ローリング）低減に剛体回転ダイナミクスが活用されている。Seakeepingジャイロ安定装置（Gyrostabilizer）は重量数トンのフライホイールを高速回転させ、コリオリ力によるカウンタートルクで横揺れを最大70%削減する。Mitsubishi Shipbuilding（三菱造船）は2019年に全長150mのフェリーへジャイロスタビライザーを搭載し、5°の波高環境での横揺れをロール角5°以内に抑えたMBD解析結果を国際船舶工学会議で発表した。
剛体動力学のトラブル対応
トラブル
慣性テンソルの符号ミスで剛体が爆発的回転する
剛体MBDで質量・慣性テンソルを手動入力した場合、慣性積（Ixy, Ixz, Iyz）の符号を誤ると非物理的な回転加速度が生じてシミュレーションが発散する。特にCADからのMass Properties読み込み時に座標系変換（右手系⇔左手系）が自動変換されず符号反転が起きるケースが多い。Adamsでは「Mass Properties Verify」コマンドで全剛体の慣性テンソルの正定値性を確認できるが、デフォルトでは実行されないため初期チェックの癖をつけることが重要だ。
なった
詳しく
報告