フレキシブルボディ — CAE用語解説
フレキシブルボディ
先生、多体動力学解析の「フレキシブルボディ」って、変形する部品をMBSに取り込む方法ですか？
理論と物理
フレキシブルボディの定義と物理的背景
「フレキシブルボディ」という用語をよく見かけますが、CAEの文脈では、剛体と具体的にどう違うんですか？
本質的な違いは、変形を考慮するかどうかです。剛体は形状が一切変わらないと仮定しますが、フレキシブルボディは外力や慣性力によって内部に応力とひずみが発生し、形状が変化します。例えば、自動車のドアミラーの振動解析では、ミラー支持部の金属板をフレキシブルボディとしてモデル化し、
変形を計算するということは、メッシュの各節点の位置が時間とともに変化するということですか？ それは計算コストが剛体に比べてかなり高くなるのでは？
その通りです。フレキシブルボディの運動方程式は、剛体の6自由度（並進3、回転3）とは異なり、メッシュの節点数×自由度分の大規模な連立方程式を解く必要があります。支配方程式は、
「フレキシブル」というとゴムのような柔らかいものだけを想像していましたが、金属板でもフレキシブルボディと呼ぶんですね。変形量が微小な場合と大きな場合で、扱い方は変わるんですか？
非常に重要なポイントです。微小変形理論では、ひずみと変位の関係が線形（
数値解法と実装
フレキシブルボディの離散化とモーダル縮約
フレキシブルボディをFEMで離散化する際、特に動解析ではどのようなソルバーを使うのが一般的なんですか？
直接積分法とモーダル法の2つが主流です。直接積分法（例：Newmark-β法、HHT-α法）は、運動方程式を時間ステップごとに直接数値積分します。Ansys Mechanicalの「Transient Structural」解析ではこれが使われます。一方、モーダル法は事前に固有値解析を行い、変形形状を固有モードの重ね合わせで表現します。例えば、
モーダル法で「自由度を削減できる」と言いますが、具体的にどのくらい計算が速くなるんですか？また、すべての変形を正確に表現するには何次モードまで必要なんですか？
例えば、元のFEMモデルの自由度が10,000（節点約3,300個）あった場合、直接積分では10,000個の連立方程式を解きます。一方、モーダル法では、重要な低次モード（例えば1次〜50次）のみを保持し、未知数を50個にまで減らせます。計算時間は数十分の一から数百分の一に短縮可能です。必要なモード数は、加振力の周波数成分によります。自動車の車体振動解析では、ISO 2631-1に基づく人間の感覚振動領域（1〜80 Hz）をカバーするために、80 Hz以下のモードを全て含めるのが一般的です。
モーダル法は線形の範囲でしか使えないと聞きました。現実の製品には接触や大変形があるはずですが、その場合はどうするんですか？
その通り、古典的なモーダル法は線形の仮定が前提です。非線形性が重要な場合、例えばAnsysの「Nonlinear Adaptive」のような手法や、完全な非線形動解析（直接積分法）に頼らざるを得ません。また、マルチボディダイナミクス（MBD）ソフトウェアでは、「フレキシブルボディ」を線形のモーダル表現でインポートし、MBDソルバー内で剛体と組み合わせてシミュレーションしますが、接触や大変形が生じる部分は「剛体」として扱い、その境界条件を慎重に設定するなどの工夫が必要です。
実践ガイド
フレキシブルボディ作成のワークフロー
実際にMBD解析用のフレキシブルボディを作成する手順を教えてください。CADデータから始める場合、何から手を付ければいいですか？
典型的なワークフローは以下の5ステップです。
インターフェースノードを「固定」して固有値解析するのはなぜですか？ 実際の製品は固定されていないのに。
良い質問です。これは「拘束モード」または「Craig-Bampton法」の考え方に基づいています。フレキシブルボディ単体の自由-自由モード（剛体モードを含む）をMBDソフトに渡すと、後で剛体と組み合わせた際に重複した剛体運動が生じ、数値的不安定の原因になります。そこで、結合点を固定した状態での「拘束モード」と、結合点を単位変位させた「静的な変形モード（附着モード）」をセットで出力します。これにより、MBDソフト側で剛体運動と弾性変形を明確に分離して扱えるのです。
メッシュの粗さは結果にどのくらい影響しますか？ また、どこを細かくメッシュするべきかの判断基準は？
メッシュサイズは、対象とする最高次のモード形状を捉えられるかどうかで決まります。経験則として、1つの波長（モードの山と谷）に少なくとも6〜8個の要素が必要です。例えば、板厚1mmの鋼板の曲げ振動で1000Hzのモードを捉えたい場合、その波長を推定し、それに応じたメッシュサイズ（例えば5mm以下）を設定します。特に、インターフェースノード周辺、断面が急激に変化する部分、応力集中が予想される部分はメッシュを細かくする必要があります。最初は中程度のメッシュで解析し、メッシュを半分にした時に固有周波数が1%以内に収束するかを確認する「メッシュ収束性解析」を行うのが確実です。
ソフトウェア比較
主要CAEソフトにおける実装の違い
Ansys、Abaqus、そして専用のMBDソフトであるAdamsでは、フレキシブルボディの扱い方にどんな違いがありますか？
役割が根本的に異なります。
COMSOL Multiphysicsは「マルチフィジックス」を売りにしていますが、フレキシブルボディの扱いはどうなんでしょう？ 振動と音響の連成解析などは得意そうですが。
COMSOLは、単一の統合環境内でフレキシブルボディの変形と他の物理場（音響、熱、流体）との連成解析を直接行える点が強みです。例えば、圧電アクチュエータで駆動するマイクロミラーの解析では、構造（フレキシブルボディ）、圧電、静電場の方程式を一度に連立して解けます。ただし、大規模な機構運動のシミュレーションという点では、AdamsのようなMBD専用ソフトに比べて計算効率が劣る場合があります。COMSOLでモーダル縮約を行い、外部ツールと連携する機能も備えています。
.mnfファイルの中身は何が入っているんですか？ また、ソフト間で互換性はあるんですか？
.mnf (Modal Neutral File) は、MSC Software（Adamsの開発元）が定義したフォーマットで、主に以下のデータを含むテキストファイルです：
トラブルシューティング
よくあるエラーとその対策
MBD解析でフレキシブルボディをインポートしたら、シミュレーションが発散したり、異常に大きな変形が発生したりしました。原因として考えられることは？
最も多い原因は3つです。
固有値解析を実行したら、想定していたはずの「1次曲げモード」の前に、非常に低い周波数（例：0.01 Hz）の奇妙なモードが出てきました。これは何ですか？
それはほぼ間違いなく「剛体モード」です。拘束が不十分で、部品がごくわずかな剛性（数値誤差レベルの剛性）でふらつくモードとして現れます。インターフェースノードに付与した固定拘束が、CADの位置とメッシュの節点位置でずれていないか確認してください。また、複数の部品をボルト結合する場合、結合面全体を「剛体結合（RBE2）」などで適切に拘束し、そのマスターノードを固定する必要があります。
フレキシブルボディを含むMBD解析は計算が重いです。高速化するための実践的なテクニックはありますか？
いくつかの効果的な方法があります。
関連トピック
なった
詳しく
報告