CAEとは? — 定義・メリット・活用事例をわかりやすく解説
理論と物理
CAEの定義と基本概念
CAEって、CADで作ったモデルを解析するソフトのことですよね? 具体的に何を「計算」しているんですか?
本質的には、物理現象を記述する「支配方程式」をコンピュータで解いています。例えば構造解析なら、材料の変形と応力を求めるための釣り合い方程式です。具体的には、物体内部の微小な立方体(要素)に働く力を考え、そのつり合い式
方程式をそのまま解けないのはなぜですか? 高校で習う連立方程式みたいにできないんですか?
形状が複雑だと、境界条件が複雑で解析解が得られないからです。例えば、自動車のエンジンブロックに複数のボルトで締結力を加えた時の応力分布は、手計算では絶対に無理です。そこで、領域を小さな三角形や四角形(有限要素)に分割し、各要素で近似的に方程式を立て、全体で数万〜数百万の連立方程式として解く「離散化」を行います。これがFEM(有限要素法)の核心です。
「近似的に」解くということは、答えは厳密ではないんですか? 誤差はどれくらいあるんですか?
その通りで、誤差は必ず生じます。主に「離散化誤差」と「数値誤差」です。メッシュを細かくする(h-refinement)か、要素の次数を上げる(p-refinement)ことで、解は真の解に収束します。実務では、メッシュサイズを半分にして結果が1〜2%しか変わらなくなった時を「メッシュ収束」と判断します。例えば最大応力が350MPaで、メッシュを細かくしても353MPaなら、実用上は収束したとみなします。
数値解法と実装
FEMの実際の計算プロセス
ソフトが「メッシュを切る」と言いますが、その後、内部でどんな行列計算が行われているんですか?
各要素の剛性マトリックス$$ [k_e] $$を組み立てて、全体剛性マトリックス$$ [K] $$を作成します。これに境界条件を適用し、連立方程式
100万×100万の方程式をどうやって現実的な時間で解くことができるんですか? 全部を一気に計算するんですか?
いいえ、直接法(例えばLDLT分解)と反復法(例えば共役勾配法)を使い分けます。直接法は安定ですがメモリ消費が大きく、反復法は大規模問題に適しますが収束性が課題です。Ansys Mechanicalでは、デフォルトで疎行列直接ソルバー「Sparse」が使われ、メモリが足りない場合には反復ソルバー「PCG」が推奨されます。GPUを活用したソルバーも増えており、NVIDIAのCUDAライブラリを使うと、特定の線形ソルバーが10倍以上高速化される事例があります。
非線形解析、例えば接触や塑性変形がある場合は、このプロセスはどう変わるんですか?
荷重や変位を少しずつ増加させる「増分解析」を行い、各増分ステップで釣り合いを反復計算(Newton-Raphson法など)で達成します。式で書くと、$$ [K_T]_i \{\Delta u\} = \{F\}_{ext} - \{F\}_{int}_i $$
実践ガイド
信頼性の高い解析を行うために
解析を始める時、最初に何を確認すべきですか? CADデータを読み込んだらすぐにメッシュを切っていいんですか?
絶対にダメです。まず「幾何チェック」が必要です。インポート時に面の欠損や微小なエッジ、重複面がないか確認します。次に「単位系の一貫性」です。CADがmm、材料データがMPa、密度がtonne/mm^3になっているか。よくある失敗は、長さ単位をmmでモデリングしながら、材料のヤング率を鋼の210GPa (=210,000 MPa)の代わりに210と入力してしまうことです。これでは結果が1000倍違ってきます。
境界条件の設定で、特に気をつけるポイントは? 「固定」と「ピン」はどう使い分けるんですか?
実物をよく観察することです。ボルトで締結されたフランジは「固定」ではなく、「ボルト穴の周りだけ変位を拘束し、面圧による摩擦を考慮する」かもしれません。また、構造物が地面に置かれている場合、垂直方向の変位だけを拘束する「ローラー支承」が適切な場合があります。安易な完全固定は、実際より剛性を高く見積もり、応力を過小評価する危険があります。ISO 10791-7などの工作機械の試験規格には、実機に近い境界条件の設定例が記載されています。
結果が出たら、どうやってその「正しさ」を判断するんですか? 応力が材料の降伏強度より低ければOK?
それだけでは不十分です。少なくとも以下のチェックリストを通します:1) 反力のつり合い(印加荷重の総和と反力の総和が一致するか)、2) メッシュ依存性の確認(メッシュを細かくして結果が大きく変わらないか)、3) 変形モードの物理的妥当性(変形が常識的に考えておかしくないか)。さらに、応力集中部の要素形状(アスペクト比、スキュー角)を確認します。JIS B 9951(構造解析の検証と妥当性確認の指針)が参考になります。
ソフトウェア比較
主要ソフトウェアの特徴と選択基準
Ansys、Abaqus、COMSOLはどれも有名ですが、根本的に何が違うんですか? 値段ですか?
値段も違いますが、哲学と強みが違います。Ansysは「個別最適化された専用ソルバーのスイート」です。構造はMechanical、流体はFluent/CFX、電磁界はMaxwellと、それぞれ独立した最強のソルバーを持ち、System Couplingで連成します。一方、Abaqus/Standard(Implicit)と/Explicit(Explicit)に分かれ、特に非線形接触問題(例えばタイヤの接地)や破壊解析に強い。COMSOL Multiphysicsは「1つの環境で任意の物理場を組み立てられる」のが売りで、独自のPDE(偏微分方程式)も直接記述できます。
自動車業界ではAbaqusがよく使われると聞きますが、なぜAnsysじゃないんですか?
歴史的な経緯と、自動車の主要な非線形問題への適性です。Abaqusはもともと非線形構造解析に特化して開発され、ゴム材料の超弾性モデルや、ボディ部品の大変形・接触解析に強みがあります。また、Explicitソルバーによる衝突解析も業界標準です。トヨタやホンダなどは長年Abaqusを使い、社内の材料ライブラリや解析手順書(ノウハウ)が蓄積されています。ただし、電動化に伴い、Ansys Maxwell(モーター解析)やFluent(電池熱流体)の採用も増えています。
無料や安価なCAEソフト(例えばFreeCADのFEMワークベンチやCalculiX)は実務で使えるレベルですか?
用途によります。CalculiXはAbaqusと類似の入力形式を持ち、線形静解析や単純な非線形解析なら実用可能です。しかし、商用ソフトのような高度な接触アルゴリズム(例えばサーフェス・トゥ・サーフェス接触の安定性)や、豊富な材料モデル(複合材、超弾性、クリープ)、そして何より技術サポートがありません。試作前の最終検証にはリスクが高すぎますが、設計初期のコンセプト検討や学生の学習用には非常に優れています。また、オープンソース故に計算過程を全て追えるのは教育上大きな利点です。
トラブルシューティング
よくあるエラーとその対処法
解析を実行したら、「負のヤコビアン」や「過度の変形」でエラーになりました。どうすればいいですか?
これは非線形解析で最も頻出するエラーです。原因と対策は:1) メッシュ品質:特に変形が大きい領域で、要素が極端に歪んでいないか確認。アスペクト比20以上は危険。2) 材料モデル:大変形を起こすのに線形弾性材料を使っていないか。ゴムなら超弾性モデルが必要。3) 荷重・境界条件:現実にはあり得ない大きな荷重や拘束をかけていないか。まずは荷重を1/10にして実行し、収束するか試します。Abaqusでは「ステビリゼーション」や「自動増分」のパラメータを調整します。
接触解析で、部品がお互いを「貫通」してしまいます。接触設定で「 bonded」や「frictional」は選んでいるのに。
主な原因は3つです。第一に「初期貫通」。解析開始時点でメッシュが既に重なっていると、ソルバーが修正できずに暴走します。Ansysでは「Adjust to Touch」機能で微調整します。第二に「接触面の定義不足」。従属面(slave)のメッシュが主面(master)より粗すぎると、隙間からすり抜けます。第三に「剛体運動」。接触しているだけの部品が、他の方向に拘束されていないと、つり合いが取れずに動き出します。全ての部品に適切な拘束があるか確認してください。
固有値解析(モード解析)で、変形モードが明らかに物理的にあり得ない、グチャグチャな形になります。なぜですか?
それは「剛体モード」が混入している可能性が高いです。構造物が十分に拘束されていない(浮いている)状態で解析すると、変形しない剛体移動(並進・回転)が周波数0Hzのモードとして現れます。ソルバーが数値誤差でこれを低い非ゼロ周波数で計算し、変形モードと一緒に出力してしまうのです。対策は、すべての剛体運動を抑える境界条件を適用すること。例えば、3点支持(1点でXYZ固定、別の1点でYZ固定、さらに別の1点でZ固定)が一般的です。結果を確認する時は、まず最低次のモードの振動数が0に近くないか、変形が剛体移動に見えないかをチェックします。
なった
詳しく
報告