拓扑优化（SIMP法）
拓扑优化（SIMP法）的理论基础
什么是拓扑优化？
老师，什么是拓扑优化？
拓扑优化是优化设计区域内材料的有无（0/1）。自动决定在哪里开孔、在哪里保留材料。由Bendsøe & Kikuchi于1988年提出。
SIMP法
SIMP（Solid Isotropic Material with Penalization）是最广泛使用的拓扑优化方法。为每个单元分配设计变量 $\rho_e$（0~1的密度）：
通过 $p$（惩罚指数，通常$p = 3$）抑制中间密度，使其趋近于0/1。
优化问题
典型的数学表述：
是"在将材料控制在$V^*$以下的同时，找到最刚硬的结构"对吧。
正是如此。通过FEM计算每次迭代的位移→计算灵敏度（各单元密度变化时目标函数的变化）→更新密度→迭代直至收敛。
总结
SIMP法的"SIMP"由Bendsoe（1989）命名
拓扑优化的代表方法SIMP（Solid Isotropic Material with Penalization），是Bendsoe（1989年）在简化Bendsoe & Kikuchi（1988年）的均匀化法后提出的方法。将每个单元的密度ρ设为0~1的连续变量，通过刚度Eρ^p表示，惩罚参数p抑制中间密度，从而得到接近0或1的清晰材料布局。名称的由来是Bendsoe在1989年的论文标题中写了"Solid Isotropic..."，后来被缩写成首字母缩写。
数值解法与实现
SIMP法的算法
1. 将所有单元的初始密度设为 $\rho = V^*/V_{total}$
2. 用FEM计算位移和应力
3. 计算灵敏度 $\partial C / \partial \rho_e$（伴随法）
4. 更新密度（OC法或MMA法）
5. 迭代直至收敛（通常50~200次迭代）
求解器
总结
无密度过滤的SIMP会产生棋盘格图案
众所周知，如果不使用密度过滤执行SIMP拓扑优化，会出现相邻单元交替为ρ=0/1的"棋盘格图案（checkerboard pattern）"这一数值病态。Bourdin（2001年）提出的Helmholtz PDE（偏微分方程）滤波器可以自然地控制最小构件尺寸（rmin），现已作为标准实现被集成到OptiStruct、Tosim和ABAQUS中。通过将rmin的设置与制造约束（最小壁厚、拔模斜度）对应起来，可以同时管理设计性和可制造性。
拓扑优化（SIMP法）拓扑优化（SIMP法）实践指南
拓扑优化的实务
汽车轻量化（支架、悬架臂）、航空航天（结构件）、3D打印（自由形状）。
实务检查清单
空客A380的机翼安装支架是SIMP优化的杰作
空客A380的客舱天花板面板安装支架（2006年首飞）作为使用OptiStruct进行SIMP拓扑优化设计的零件在业界非常有名。相比传统手工设计品，在满足疲劳寿命约束的同时重量减轻了30%，并获得了Altair的EngineeringImpact奖。目前OptiStruct作为标准的拓扑优化工具在空客所有机型中使用，每年有超过1000个零件优化通过此工具实施。