寸法最適化 — CAE用語解説
寸法最適化
理論と物理
寸法最適化の基本概念
寸法最適化って、具体的に何を変えているんですか？形状そのものは変えないんですよね？
その通りです。形状のトポロジーや外形は固定し、設計変数として定義した「寸法」を変化させます。例えば、トラスの各棒材の断面積、板厚、梁の高さや幅などです。目的は、質量最小化や応力均一化、固有振動数制御などです。制約条件として、応力が材料の降伏強度345MPaを超えない、たわみがスパンの1/500以下、などが典型的です。
「寸法」を変数として扱うということは、設計変数の数は最初から決まっているんですか？
はい、最適化問題の規模は初期設計で決まります。10個の板厚を変数にすれば10次元の最適化問題です。ここがトポロジー最適化との大きな違いで、後者は要素の有無が変数なので、設計領域内で最適な形状が「創生」されます。寸法最適化は、既存の設計案のパラメータチューニングに近いですね。
最適化の数学的な定式化はどうなりますか？
一般的な定式化は以下の通りです。
数値解法と実装
感度解析と最適化アルゴリズム
最適化アルゴリズムはどう選ぶんですか？勾配法と直接探索法の違いは？
変数が多く（数十〜数百）、目的・制約関数が滑らかな場合、勾配ベースの方法が効率的です。例えば、MMA (Method of Moving Asymptotes) やSQP (Sequential Quadratic Programming) がよく使われます。これらの方法は、目的関数や制約関数の設計変数に対する感度
感度解析って、FEMの結果からどうやって計算するんですか？
構造最適化で応力制約の感度を求める場合、随伴変数法が一般的です。全体剛性方程式
設計変数が離散値、例えばJIS規格の板厚（1.6, 2.0, 2.3mm）しか選べない場合は？
それは離散変数最適化問題になり、難易度が跳ね上がります。実務では二段階アプローチを取ることが多いです。まずは連続変数として最適化し、得られた最適寸法に最も近い規格サイズを選択します。その後、その離散値で強度や振動を再評価し、問題があれば一つ上のサイズに切り上げます。専用の離散最適化アルゴリズム（遺伝的アルゴリズム等）もありますが、計算コストが非常に高くなります。
実践ガイド
ワークフローと設計変数の定義
実際に寸法最適化を始めるとき、最初に何を決めるべきですか？
まず「設計変数」を明確に定義します。例えば、自動車ドアの内板を5つのパネルに分け、それぞれの板厚を変数とする。次に「設計変数リンク」を設定します。左右対称なら、右側の板厚を左側の変数にリンクさせ、変数を減らします。上下限は、製造可能な最小板厚0.7mm（自動車ボディ）と、剛性確保のための最大板厚2.3mmなど、実用的な範囲で設定します。
目的関数と制約条件の設定で、陥りやすい失敗はありますか？
よくあるのは「制約条件が多すぎる」または「矛盾する」場合です。例えば、「質量を最小化せよ」「かつ、第1次固有振動数を50Hz以上にせよ」「かつ、最大応力を200MPa以下にせよ」「かつ、最大たわみを2mm以下にせよ」と全てを厳しい制約にすると、実行可能解が存在せず最適化が失敗します。まずは最重要制約（例えば応力）だけを課し、最適解を見てから、他の条件を緩和しながら追加していくのが現実的です。
最適化が収束した結果をどう評価すれば信頼できますか？
必ず「最適解」のモデルで通常のFEM解析を実行し、結果を独立に検証します。最適化プロセス中の近似（例えば、応力制約を特定点のみで評価している等）があるためです。また、感度がほぼゼロでない「アクティブな制約」を確認します。例えば、板厚が下限値0.7mmに張り付き、かつその部位の応力が制約値の200MPaぴったりなら、それは有効な制約です。板厚が下限なのに応力が100MPaしかなければ、さらに薄くできる可能性があります。
ソフトウェア比較
各ソルバーの特徴と適用事例
Ansys、Abaqus、Altair OptiStructでは、寸法最適化のアプローチに違いはありますか？
大きな違いがあります。OptiStructは元来、トポロジー/形状/寸法最適化に特化したソルバーで、特に寸法最適化では複数の荷重ケースや動的応答制約を効率的に扱えます。Abaqusは基本解析が主で、最適化にはIsightなどのプロセス自動化ツールと連携するか、Abaqus/CAEの「パラメトリック研究」を利用します。Ansys Workbenchでは「Goal Driven Optimization」モジュールや「DesignXplorer」が使われ、応答曲面法を用いて最適化するのが特徴です。
応答曲面法とは何ですか？直接FEMで評価するのとどう違うんですか？
応答曲面法(RSM)は、設計空間内の限られたサンプル点（例えば、各変数の最小、中央、最大値の組合せ）でFEM解析を実行し、その結果から目的関数や制約関数を近似する多項式（応答曲面）を構築します。その後、最適化はこの軽量な近似モデル上で実行されるため、非常に高速です。Ansys DesignXplorerはこれを多用します。一方、OptiStructに代表される勾配法は、各反復で実際のFEM解析と感度解析を行うため、より正確ですが、反復回数分の計算コストがかかります。
無償・低価格のツールで寸法最適化できるものはありますか？
はい、あります。例えば、オープンソースのCalculixには最適化モジュールはありませんが、Pythonスクリプトでパラメータを変化させてCalculixを繰り返し実行し、SciPyの最適化ライブラリ（例えば`scipy.optimize.minimize`）で制約付き最適化を行うことは可能です。また、SALOMEプラットフォームとCode_Asterを組み合わせ、`EFICAS`というモジュールを使って最適化設定を行う方法もあります。ただし、商用ソフトのようなGUIや高度な感度解析、収束管理はほぼ手作業になるため、技術的ハードルは高いです。
トラブルシューティング
よくあるエラーと対策
最適化を実行したら「設計可能解が見つかりません」というエラーが出ます。まず何を疑うべきですか？
まず、初期設計が制約条件を満たしているか確認してください。初期モデルですでに応力オーバーやたわみ過大があれば、そもそも実行可能領域の外からスタートしていることになります。次に、制約条件が矛盾していないかチェックします。例えば、「板厚を1.0mm以上2.0mm以下」としつつ、「質量を5kg以下」かつ「一階固有振動数を100Hz以上」という制約は、物理的に達成不可能かもしれません。制約を一つずつ緩和（100Hz→80Hz）して試す必要があります。
最適化が振動して、なかなか収束しません。板厚の値が反復ごとに大きく変動します。
これは「制約条件が急峻に変化する」場合に起こります。例えば、ある板厚をわずかに下回ると、その部分の座屈モードが突然現れ、座屈荷重が劇的に低下する場合などです。対策は二つあります。1) 移動漸近線の初期値や減衰係数（MMAアルゴリズムの場合）を調整してステップ幅を小さくする。2) 問題の性質を変える。座屈制約があるなら、板厚の下限を実質的に上げるか、座屈制約そのものを外して、最適化後に別途座屈検証を行う。
得られた最適解が「局所最適解」に陥っている気がします。大域的最適解を探す方法は？
勾配法は基本的に局所解を探すので、初期値を変えて複数回実行する「マルチスタート」が現実的な対策です。OptiStructなら`OPTIMIZE`カードで`DESMAX`を増やし、初期設計をランダムに変えて再実行できます。より確実なのは、大域的最適化アルゴリズム（遺伝的アルゴリズム、粒子群最適化など）を使うことですが、計算コストは数百倍になります。実務では、設計者の経験から「あり得そうな」複数の初期案を用意し、それぞれから最適化をスタートさせるのが効率的です。
最適化後の形状が、製造や組み立てを考慮していない、変なものになることは？
寸法最適化では形状は変わらないので「変な形状」にはなりませんが、製造不可能な結果にはなります。例えば、隣接するパネルの板厚が最適化で0.8mmと2.2mmと大きく異なり、溶接や接着に支障が出る場合です。これを防ぐには「製造制約」を追加します。OptiStructでは`DEQCON`カードで複数の板厚変数を等しくする、または`DLINK`カードで比例関係を持たせることができます。Ansysでは「パラメータ相関」を設定します。要は、最適化問題に「現実の縛り」を事前に組み込むことが肝心です。
関連トピック
なった
詳しく
報告