📚 课程目录第4章拓扑优化核心方法

4. 拓扑优化核心方法

在拓扑尺寸形状联合优化中，拓扑优化负责回答“材料该放在哪里”，尺寸优化负责回答“截面该多大”，形状优化负责回答“边界该长成什么曲线”。三者要协同，前提是拓扑优化本身能稳定地产生清晰、可制造、可继续被尺寸与形状模块接手的几何。SIMP、水平集和相场是当前最主流的三类拓扑优化方法，它们对“材料有无”的数学表达方式不同，因而在数值实现、边界描述能力和联合优化中的角色也各不相同。

SIMP：以密度为设计变量的惩罚插值法

SIMP（Solid Isotropic Material with Penalization）把每个有限元单元的“相对密度”ρ_e ∈ [0,1] 作为设计变量，并用幂函数惩罚中间密度，使材料刚度按 E_e = ρ_e^pE₀ 插值，通常取惩罚因子 p=3。ρ_e→1 表示实体，ρ_e→0 表示孔洞，中间密度因刚度收益低而被优化器自然推向两端。其典型优化列式为：

min  c(ρ) = U^T K U
s.t. K(ρ) U = F
     V(ρ) = Σ ρ_e v_e ≤ V*
     0 ≤ ρ_e ≤ 1

数值实现上，SIMP 依赖三个关键环节。第一是灵敏度分析，柔度对密度的导数可由单元应变能直接得到：∂c/∂ρ_e = −p ρ_e^p−1 u_e^Tk₀u_e，配合伴随法即可高效求解。第二是密度过滤，用半径 r_min 的加权平均抑制棋盘格和网格依赖，例如 r_min=1.5～3 倍单元尺寸。第三是优化器，常用 OC 准则更新或 MMA 处理多约束。SIMP 的优点是梯度信息完整、实现成熟、易与尺寸变量耦合；缺点是边界呈灰度过渡，需要额外阈值化或后处理才能得到清晰轮廓。

水平集方法：用隐式函数追踪边界

水平集方法把结构边界表示为高维标量函数 φ(x) 的零等值线：φ>0 为材料，φ<0 为孔洞，φ=0 即边界。边界随速度场 V_n 沿法向演化，满足 Hamilton–Jacobi 方程 ∂φ/∂t + V_n|∇φ| = 0。速度场由形状导数驱动，例如以柔度最小为目标时，V_n 与应变能密度和拉格朗日乘子相关。

数值实现通常采用有限差分或有限元离散 φ，用迎风格式和 CFL 条件稳定推进，并每隔若干步重新初始化 φ 为符号距离函数，防止梯度畸变。水平集的突出优势是边界始终清晰、可直接输出光滑轮廓，天然适合与形状优化衔接；其挑战在于拓扑变化能力受初始孔洞分布限制，且对速度场符号和重新初始化频率较敏感。在联合优化中，它常作为“边界精修器”，把 SIMP 得到的灰度结果转化为显式边界后继续做形状微调。

相场方法：用连续场变量描述材料与孔洞

相场方法引入连续相场变量 c(x) ∈ [0,1]，以 c=1 表示材料、c=0 表示孔洞，并在目标函数中加入界面能项，使材料与孔洞之间的过渡带具有可控宽度 ε。其典型形式为：

min  J = 结构柔度 + λ ∫ [ε|∇c|² + (1/ε)c²(1−c)²] dΩ
s.t. 平衡方程与体积约束

相场法的数值实现与 SIMP 类似，仍可用有限元和梯度优化器，但多了界面能正则项，因此边界宽度由 ε 控制，网格需足够细以解析过渡带。它的优点是拓扑变化自然、边界较清晰、数学上接近连续介质描述；代价是计算量较大，且 ε 与网格尺寸需匹配，否则会出现边界模糊或数值不稳定。在联合优化中，相场常被用于需要平滑边界和可控界面宽度的场景，例如增材制造中的最小特征尺寸控制。

三类方法在联合优化中的分工与衔接

从联合优化的流程看，SIMP 通常承担“拓扑主生成”的角色：它能在固定设计域内自由增减材料，快速给出接近最优的材料分布，并把每个单元的密度传递给尺寸优化模块，用于确定杆件截面或板厚。水平集则更适合“边界精修”：当 SIMP 结果经过阈值化得到初始轮廓后，用水平集演化可以消除灰度过渡带，得到光滑、可参数化的边界，再交给形状优化做局部微调。相场方法介于两者之间，它既能自然改变拓扑，又能通过界面能项控制边界宽度，因此在需要同时兼顾拓扑自由度和制造约束（如最小特征尺寸、悬垂角）的联合优化中很有价值。

实际工程中，一种常见的组合策略是“SIMP 粗拓扑 + 水平集精边界 + 尺寸/形状联合微调”。例如某悬臂梁优化示例：设计域 120 mm × 40 mm，左端固定，右端下角施加 500 N 集中力，体积分数上限 0.4，SIMP 取 p=3、r_min=2 mm，迭代 80 步后得到灰度拓扑；阈值化后转入水平集，以柔度为目标演化 30 步，边界灰度带从约 3 个单元宽收敛到 1 个单元以内；最后对主要杆件截面尺寸和边界控制点做联合优化，柔度可再下降约 5%～8%。这些数值属于教学示例的合理假设，用于说明方法衔接的效果趋势，并非特定软件实测结果。

选择方法时可按需求判断：需要自由拓扑且与尺寸变量同框架，优先 SIMP；需要清晰边界并接形状优化，优先水平集；需要同时控制拓扑变化和界面宽度，优先相场。三者并非互斥，联合优化中常按“粗拓扑—精边界—细尺寸”的顺序组合使用。

无论采用哪种方法，联合优化都要保证设计变量、灵敏度信息和几何描述在模块之间一致传递。SIMP 的密度场、水平集的 φ 场、相场的 c 场，最终都要映射为可被尺寸和形状模块识别的几何参数或单元属性。理解这三类方法的数学本质与数值特性，才能在联合优化中选对主方法、设计好衔接接口，避免因边界模糊、拓扑突变或参数失配导致优化震荡或结果不可制造。