SIMP位相最適化シミュレーター シミュレーター集へ
構造最適化ツール

SIMP位相最適化シミュレーター

SIMP法(密度法)による2次元トポロジー最適化をブラウザ上でリアルタイム実行。荷重・境界条件・体積制約を変えて、最適な材料配置を直感的に体験できます。

$$E(\rho_e) = E_{\min} + \rho_e^p(E_0 - E_{\min}),\quad \min_{\boldsymbol\rho}\; \mathbf{f}^T\mathbf{u} \quad \text{s.t.} \sum\rho_e v_e \le V^*$$
パラメータ設定
境界条件预设
体積制約 V* [%] 40 %
材料を使える割合。小さいほど軽量・難解に
メッシュ解像度
荷重方向
イテレーション0 / 50
コンプライアンス
体積分率
変化量 Δρ
SIMP法の概要
各要素に密度 ρ∈[0,1] を与え、
剛性を E = Emin + ρp·E₀ とする(p=3)。
OC法で体積制約を満たしながら
コンプライアンス(柔さ)を最小化。
コンプライアンス
正規化値
体積分率
%
メッシュ
40×20
要素
固体要素数
密度分布(位相最適化結果)
片持ち梁 Iter: 0
固体 (ρ≈1)
中間密度
空洞 (ρ≈0)
コンプライアンス収束
理論 — SIMP位相最適化

SIMP材料モデル(ペナリゼーション)

$$E_e(\rho_e) = E_{\min} + \rho_e^p (E_0 - E_{\min}),\quad p=3$$

ペナルティ $p=3$ で中間密度を抑制し、0/1 の明確な構造に誘導する。

最適化問題の定式化

$$\min_{\boldsymbol\rho}\; c = \mathbf{f}^T\mathbf{u} \quad \text{s.t.}\; \mathbf{K}\mathbf{u}=\mathbf{f},\; \sum_e\rho_e v_e \le V^*$$

コンプライアンス(外力仕事)を最小化 = 最も剛な構造を探索。

感度解析(adjoint法)

$$\frac{\partial c}{\partial \rho_e} = -p\,\rho_e^{p-1}\,\mathbf{u}_e^T \mathbf{K}_e^0 \mathbf{u}_e$$

各要素の変位から感度を計算。密度が高い要素ほど感度(貢献度)が大きい。

OC法(最適性基準法)更新則

$$\rho_e^{new} = \rho_e \cdot B_e^{\eta},\quad B_e = \frac{-\partial c/\partial\rho_e}{\lambda\,\partial V/\partial\rho_e}$$

ラグランジュ乗数 λ を二分法で求め、体積制約を満たすよう密度を更新。

相关工具

カテナリー曲線・ケーブル張力計算機
自重を受けるケーブルのカテナリー形状・張力分布・最大たわみを計算
圧力容器応力(Lamé方程式)
Lamé方程式で圧力容器の応力分布を計算
降伏条件比較ツール — フォンミーゼス・トレスカ
フォンミーゼス・トレスカ・モールクーロンの降伏条件を主応力空間で比較
静力学・力とモーメントのつり合い計算機
2D物体に作用する複数の力のつり合い条件を計算