水平集法拓扑优化
水平集法拓扑优化的理论基础
什么是水平集方法
老师，水平集方法和SIMP法有什么区别？
SIMP法通过每个单元的密度（0~1）来表示材料的有无。水平集方法则是通过隐函数 $\phi(\mathbf{x})$ 直接追踪边界（形状）。$\phi > 0$: 有材料，$\phi < 0$: 空洞，$\phi = 0$: 边界。
比SIMP法的边界更清晰吗？
是的。SIMP法存在灰色单元（中间密度）问题，而水平集方法边界始终清晰。不过，它不擅长孔洞成核（生成新孔洞）。
总结
水平集方法是Osher和Sethian（1988年）的发明
水平集方法是斯坦利·Osher和James Sethian（加州大学伯克利分校）于1988年在JCP期刊上发表的"界面追踪法"。最初是作为火焰传播和水面波数值模拟方法开发的，但通过Sethian在2000年代的著作，其在计算机视觉、医学图像处理和拓扑优化中的应用迅速扩展。水平集在结构优化中的应用由Wang等人（2003）和Allaire等人（2004）独立发表，其特点是相比SIMP法能自然地获得平滑边界。
数值解法与实现
水平集方法的FEM
Abaqus TOSCA:
*TOPOLOGY OPTIMIZATION, LEVELSET 设置。用Hamilton-Jacobi方程更新边界。
总结
Hamilton-Jacobi方程驱动水平集界面
在水平集方法拓扑优化中，将材料·空洞界面作为零水平集面进行追踪，通过Hamilton-Jacobi方程（∂φ/∂t + v|∇φ|=0）使界面随时间演化。将形状灵敏度（shape gradient）代入速度场v，界面会自动向改善目标函数的方向移动。NaN传播或速度场符号错误导致的"界面内外反转"错误是实现时最常遇到的问题，带符号距离函数的周期性重新初始化（reinitialization）是稳定化的关键。
水平集法拓扑优化水平集法拓扑优化实践指南
水平集的实务
边界清晰，易于CAD转换，因此适用于3D打印或精密机械零件的优化。
实务检查清单
时尚家具品牌Kartell的水平集优化椅子
意大利高端塑料家具品牌Kartell在其2019年发布的"Masters"系列新型椅子中采用了水平集拓扑优化。针对人体工程学载荷工况（100kg体重就座+横向冲击），使用Altair Inspire（原solidThinking Inspire）进行水平集优化，优化聚碳酸酯椅子的壁厚分布，在同等强度下将零件重量从290g减至210g。"有机流动的形状"也被采纳为商品设计，在米兰设计周成为话题。
水平集法拓扑优化软件与求解器比较
工具
各公司水平集优化实现比较
水平集方法拓扑优化的商业实现先后有Altair OptiStruct（2012年~）、COMSOL Multiphysics 5.4（2018年~）、Simulia Tosca（2020年~）。OptiStruct因其AMOS（自适应形态优化策略）能自动应用制造约束（最小板厚·拔模斜度）而在汽车行业获得好评，被用于丰田制造的转向节臂轻量化设计。
尖端技术
水平集的尖端
水平集方法的起源：1988年的Osher-Sethia论文
拓扑优化的水平集方法基于斯坦利·奥舍和詹姆斯·塞西亚于1988年在Journal of Computational Physics上发表的界面追踪算法。与传统的SIMP法不同，其边界清晰，易于评估制造性，因此被Altair solidThinking用于空客A350钛合金支架（重量从1.2kg减至0.7kg，轻量化42%）的设计。
水平集法拓扑优化常见问题与调试
水平集的故障
水平集方法的数值不稳定性对策
在水平集方法拓扑优化中，"哈密顿-雅可比方程的数值扩散"导致的边界消失是典型的收敛失败模式。COMSOL 6.0实现了每5次优化迭代自动设置重新初始化周期的自适应重新初始化，报告称在与Windchill联动的设计案例中，将风扇叶片拓扑优化的收敛失败率从35%降低至8%。