/ SIMP Topology Optimization Simulator Tools Index
JP | EN | ZH
Structural Optimization Tool

SIMP Topology Optimization Simulator

Run 2D topology optimization using the SIMP (density) method live in your browser. Adjust loads, boundary conditions, and volume fraction to intuitively explore optimal material layouts for lightweight structural design.

$$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^*$$
Parameter Settings
Boundary condition preset
Volume fraction V* [%] 40 %
Fraction of the domain that can be material. Lower values = lighter and harder to solve.
Mesh resolution
Load direction
Iteration0 / 50
Compliance
Volume fraction
Change Δρ
SIMP Method Overview
Each element is assigned a density ρ∈[0,1].
Stiffness is set to E = Emin + ρp·E₀ (p=3).
The OC method minimizes compliance
(structural flexibility) subject to the volume constraint.
Compliance
Normalized
Volume fraction
%
Mesh
40×20
elements
Solid elements
count
Density Distribution (Topology Optimization Result)
Cantilever beam Iter: 0
Solid (ρ≈1)
Intermediate density
Void (ρ≈0)
Compliance Convergence
Theory — SIMP Topology Optimization

SIMP Material Model (Penalization)

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

The penalty p=3 suppresses intermediate densities, driving the design toward a clear 0/1 solid-void layout.

Optimization Problem Formulation

$$\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^*$$

Minimizing compliance (work done by external forces) finds the stiffest possible structure.

Sensitivity Analysis (Adjoint Method)

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

Sensitivities are computed from element displacements. Higher-density elements have greater sensitivity (contribution).

OC Method (Optimality Criteria) Update Rule

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

The Lagrange multiplier λ is found by bisection to update densities while satisfying the volume constraint.