SIMP法とは何ですか？

設計領域の各要素に0〜1の密度を持たせ、材料の性質を密度のp乗で表して最適な配置を探す方法です。pはペナルティ係数で、3程度がよく使われます。中間の密度は材料を使うわりに性能が低くなるので、最適化の結果が0か1に近づきます。

ペナルティ係数を1にするとどうなりますか？

中間の密度が損をしなくなり、灰色のぼやけた結果になります。40×40の放熱経路の例では、熱の通りにくさは一様に材料を置いた場合の50%まで下がりましたが、66%の要素が中間の密度で、0と1に丸めると一様配置の約19倍まで悪化しました。そのままでは作れない設計です。

マルチフィジックスで目的を組み合わせるときの注意は？

構造の剛性と放熱のように単位も大きさも違う目的を重み付きで足すと、値の大きい目的だけが効いてしまいます。それぞれを初期値で割って正規化してから重みを付け、重みを変えた複数の結果（パレート解）を比べて選ぶのが基本です。

マルチフィジックストポロジー最適化 — SIMP法で構造・熱・流体を同時最適化

カテゴリ: 連成解析 | 2026-09-30 改稿

Multiphysics topology optimization showing coupled structural-thermal-fluid density distribution

マルチフィジックストポロジー最適化 ― 構造・熱・流体を同時に考慮した材料分布の最適解

理論：材料を置くか置かないかを連続値で探す

概要

🙋

発熱する部品を支えるブラケットを、軽くて放熱もよい形にしたいです。トポロジー最適化で構造と熱を同時に最適化できますか？

🎓

できるよ。トポロジー最適化は、設計する領域を細かい要素に分け、それぞれに材料を置くか置かないかを決めて、最も良い形を探す方法だ。0か1かの選び方は組み合わせが膨大なので、SIMP法では密度を0〜1の連続値にして、材料の性質を密度のべき乗で表す。構造の剛性と熱の通りやすさの両方を目的に入れれば、両方を考えた形が得られる。ただし、連続値にしたことで生じる『中間の密度』の扱いには注意がいるんだ。

SIMP法の材料モデル

$$ k(\rho_e) = k_{min} + \rho_e^{\,p}\,(k_0 - k_{min}), \qquad 0 \le \rho_e \le 1 $$

$\rho_e$ は要素の密度、$k_0$ は材料の熱伝導率（構造ならヤング率）、$k_{min}$ は計算を安定させる小さな値、$p$ はペナルティ係数。$p = 3$ なら密度0.5の要素は材料を半分使うのに性能は1/8しかない。この『損』があるので、最適化は中間の密度を避けて0か1に寄っていく。

目的と制約

$$ \min_{\boldsymbol\rho}\ C = \mathbf{f}^{\mathsf T}\mathbf{T} \quad \text{s.t.}\quad \mathbf{K}(\boldsymbol\rho)\,\mathbf{T} = \mathbf{f},\ \ \frac{1}{N}\sum_e \rho_e \le V^* $$

熱の問題では、発熱 $\mathbf{f}$ と温度 $\mathbf{T}$ の積（熱コンプライアンス）を最小にすると、平均的に温度の低い、熱を通しやすい配置になる。$V^*$ は使える材料の割合。構造の場合は荷重と変位の積（コンプライアンス）で、剛性の高い配置を探す。

Coffee Break よもやま話

1988年の論文から自動車の部品へ

トポロジー最適化は1988年、ベンソーと菊池による均質化法の論文で本格的に始まった。材料に細かな穴があいた状態を考えて、その穴の大きさを最適化する方法だった。その後、より簡単なSIMP法が広まり、2000年代には市販のソフトに組み込まれて、自動車や航空機の軽量化に使われるようになった。金属の積層造形が普及すると、最適化で得られた複雑な形をそのまま作れるようになり、熱と構造、流体を組み合わせた最適化も実用の段階に入っている。

計算例：放熱経路の最適化

条件

正方形の板を40×40の要素に分け、全体で一様に発熱させ、左辺の中央の小さな部分だけを0 ℃に保つ（そこから熱を逃がす）。熱を通す材料を面積の40%だけ置けるとして、熱コンプライアンスを最小にした。比べる基準は、同じ量の材料を一様に広げた配置（熱伝導率0.4）とした。

結果

条件	最適化の途中の値（基準比）	中間の密度の要素	0/1に丸めた設計（基準比）
p = 1、フィルタ半径1.5、150回	50%	66%	1,905%（約19倍）
p = 3、フィルタ半径1.5、150回	75%	51%	121%
p = 3、フィルタ半径1.2、300回	71%	41%	96%

🙋

p = 1のほうが途中の値はずっと良いのに、丸めると使い物にならないんですね。

🎓

そう。p = 1では中間の密度が損をしないので、材料を薄く広く配った灰色の結果になる。熱は通るけど、『熱伝導率が半分の材料』は現実には作れない。0と1に丸めると細い経路が途切れて、一様配置の19倍も悪くなった。p = 3では中間の密度が減り、丸めたあとも一様配置より少し良い設計が残った。この小さな例でも灰色がまだ4割残っているように、実際には細かいメッシュ、密度をくっきりさせる処理（投影）、ペナルティを少しずつ上げる工夫が要る。最適化の結果は、必ず0と1に丸めた形でもう一度解析して確かめよう。

構造と熱を組み合わせる

$$ J = w_1\,\frac{C_{struct}}{C_{struct,0}} + w_2\,\frac{C_{therm}}{C_{therm,0}} $$

単位も大きさも違う目的をそのまま足すと、値の大きい方だけが最適化される。初期値で割って正規化し、重み $w_1, w_2$ を変えて複数の設計を作る。重みの違う結果を並べると、剛性を少し犠牲にして放熱を大きく改善できる点などが見えてくる（パレート解）。熱で膨張して生じる応力のように、2つの物理が直接つながる場合は、連成した解析を最適化の中で解く必要がある。

実務での進め方

手順

設計できる領域と、固定する領域（取り付け面、発熱部）を分ける。
荷重・発熱・境界条件を、実際の使われ方の最も厳しい組み合わせで与える。
ペナルティ・フィルタ・投影の設定で、灰色の少ない結果にする。
結果を0/1に丸め、滑らかな形状に作り直して、元の条件で解析し直す。
製造の制約（積層造形の張り出し、鋳造の抜き方向など）を入れて再度最適化する。

Coffee Break よもやま話

「最適化したら細い枝だらけになった」

放熱の最適化を細かいメッシュで行うと、木の枝のように細かく分かれた形が出てくることがある。熱を集めるには理にかなった形だが、メッシュを細かくするほど枝も細かくなり、作れない太さになってしまう。フィルタ半径を実際の最小の肉厚に合わせて決めると、メッシュによらない太さの枝になる。フィルタは数値の安定のためだけでなく、作れる形にするための設計条件でもあるんだ。

よくある間違い

間違いと対策

間違い	起きること	対策
灰色の結果をそのまま評価	性能を過大評価	0/1に丸めて再解析
ペナルティが小さい	作れない中間材料	p = 3程度・継続法
フィルタ半径を適当に決める	細すぎる枝・メッシュ依存	最小肉厚から決める
目的を正規化せず足す	片方だけ最適化	初期値で正規化
荷重条件が1つだけ	別の条件で弱い形	複数の条件を入れる

🙋

トポロジー最適化を自分でも試したいです。

🎓

トポロジー最適化シミュレーターとSIMP法ツールで最適化の進み方、ヒートシンク設計ツールで放熱の形の効果を試せるよ。関連する連成の話題は分離型と一体型、マルチフィジックスの不確かさ、縮約モデルの記事にまとめてある。