4. 常用拓扑优化方法对比

拓扑优化的本质，是在给定设计域、载荷与约束下，寻找材料“该放在哪里、该去掉哪里”的最优分布。不同方法的分歧，集中在两个问题上：用什么变量描述材料的有无，以及如何让边界随迭代演化。理解这两点，就能把 SIMP、水平集、相场、BESO 与 MMC 放在同一张坐标系里比较，而不是死记各自的名字。

一、五种方法的核心思想

SIMP（固体各向同性材料惩罚）采用单元相对密度 ρ_e∈[0,1] 作为设计变量，通过幂次惩罚把中间密度“推向”0 或 1，常用插值模型为 E_e=E_min+ρ_e^p(E₀−E_min)，p 通常取 3。它把离散的“有/无”问题松弛为连续优化，灵敏度可解析求导，因此能直接接入梯度类算法（如 MMA、OC）。缺点是灰度单元和网格依赖性，需要过滤或投影来抑制。

水平集方法用高维标量函数 φ(x) 的零等值线隐式表达结构边界：φ>0 为材料，φ<0 为空，φ=0 即边界。边界演化由 Hamilton-Jacobi 方程驱动，速度场来自形状导数。它的最大优势是边界始终清晰、无灰度，且能自然处理边界形状与曲率相关的问题；代价是每次迭代需重新初始化符号距离函数，且对初始孔洞布局较敏感，容易陷入局部解。

相场方法引入连续相场变量 c(x)∈[0,1] 描述材料与空洞，用扩散界面代替尖锐边界，通过 Cahn-Hilliard 型演化方程在“界面能”与“结构柔度”之间求平衡。它天然允许拓扑自动成核与合并，数值上比水平集更稳健，但界面宽度参数需要标定，且计算量偏大。

BESO（双向渐进结构优化）直接以离散单元“存在/删除”为变量，按灵敏度排序，逐步删除低效单元、并在需要处添加材料。它概念直观、结果黑白分明，非常适合工程快速迭代；但灵敏度排序的离散更新容易产生棋盘格与网格依赖，通常需配合灵敏度过滤与历史平均。

MMC（移动可变形组件）把结构显式表示为若干可变形组件的并集，每个组件由中心、长度、角度、宽度等几何参数控制。设计变量是几何参数而非单元密度，因此优化结果天然是显式几何、边界光滑、可直接参数化建模，便于与 CAD 对接；难点在于组件数量与初始布局的设定，以及组件重叠区域的灵敏度处理。

二、光机结构选型：先看约束，再选方法

光机结构（如反射镜背部支撑、镜筒桁架、柔性铰链、光学平台蜂窝板）的拓扑优化，与一般承力结构相比有三个特殊约束：一是刚度与热稳定性并重，常需同时控制柔度、热变形与一阶频率；二是可制造性要求高，减材加工、增材制造或铸造工艺对最小壁厚、拔模、悬垂有硬性限制；三是几何需可参数化交付，最终要落到 CAD 与图纸。选型建议如下。

• 以刚度/频率为目标、设计域规则、追求成熟稳定：优先 SIMP。它灵敏度解析、算法生态最完善，配合密度过滤（半径约 1.5～2 倍单元尺寸）与 Heaviside 投影，可得到较清晰的桁架式传力路径。适合镜筒、支撑背板等大规模连续体问题。
• 边界形状敏感、要求无灰度且边界光滑：选水平集。例如反射镜轻量化孔的形状优化、柔性铰链轮廓设计，边界清晰可直接提取曲线。但需给定合理的初始孔洞，否则易停在局部解。
• 需要自动成核、拓扑变化剧烈、追求数值稳健：选相场。适合初始布局难以预判、孔洞需自发产生与合并的场景，但计算成本高，需先标定界面宽度。
• 快速工程迭代、结果要黑白分明、便于直接制造：选 BESO。对桁架、支腿、加强筋布局尤其直观，配合灵敏度过滤与历史平均可抑制棋盘格。
• 最终要交付参数化 CAD、组件式结构（如支撑杆、肋板、柔性单元）：选 MMC。几何参数即设计变量，优化结果可直接驱动 CAD 重建，适合光机中由杆、板、铰链组成的装配式结构。

三、一个可复现的选型判断示例

设某口径 300 mm 反射镜背部支撑结构，设计域 300 mm×300 mm×40 mm，材料为碳化硅（E≈410 GPa，ρ≈3.1 g/cm³），约束为背部三点固定，载荷为自重加 5 ℃ 均匀温升，目标是最小化镜面 RMS 位移并约束一阶频率不低于 800 Hz、最小壁厚不小于 3 mm。可按以下顺序判断：

1. 设计域规则、单元规模大、需多目标（柔度+热+频率），首选 SIMP，用加权柔度或归一化多目标，密度过滤半径取 2 倍单元尺寸（若单元 2 mm，则过滤半径 4 mm）。
2. 若优化后出现大量灰度、边界难以提取，切换到水平集，以 SIMP 结果为初始 φ 场，做边界精修。
3. 若结构最终由若干支撑杆/肋组成、需参数化交付，改用 MMC，把每根支撑杆设为一个组件，优化其中心、长度、角度与宽度。
4. 若迭代周期要求短、只需快速给出加强筋布局，用 BESO，体积分数从 0.4 逐步降到 0.2，每步进化率取 1%～2%。

判断结果是否可用，重点看三条：目标是否收敛且无剧烈振荡；结构是否满足最小壁厚与工艺约束；灰度或中间密度单元占比是否可接受（SIMP 一般希望低于 5%～10%）。若灰度偏高，先检查过滤半径与惩罚因子 p，再考虑换方法，而不是盲目增加迭代步数。

把五种方法放在“设计变量—边界表达—结果特征”这条链上看，选型就变成了一道有依据的判断题：连续密度类（SIMP、相场）擅长探索，隐式边界类（水平集）擅长精修，离散类（BESO）擅长快速迭代，显式几何类（MMC）擅长交付。光机结构往往需要其中两到三种接力使用，而不是只押一种。