3、拓扑优化入门：基本概念、变密度法原理与优化三要素

各位好，欢迎来到拓扑优化入门这一讲。

说实话，在接触反射镜支撑结构之前，我对拓扑优化的理解还停留在“把材料挖掉一些”的层面。直到有一次，我接手一个口径300mm的反射镜项目，重量指标卡得死死的，传统经验设计怎么都减不下来。那时候我才真正意识到——拓扑优化不是简单的“挖洞”，而是一门在给定空间里寻找最优材料分布的科学。

这一讲，我们就从最基础的概念讲起，把变密度法的原理和优化三要素说清楚。

3.1 什么是拓扑优化？

拓扑优化，说白了就是：在给定的设计区域内，告诉你哪里该有材料，哪里该是空腔。

你想想看，传统的尺寸优化是“把一根梁从10mm改成8mm”，形状优化是“把圆角从R5改成R8”。但拓扑优化不一样——它直接改变结构的连通性和布局，自由度最高，减重潜力也最大。

我在做某型号空间相机的主镜支撑时，初始设计域是一个实心圆柱。拓扑优化跑完后，材料自动汇聚成三条120°均布的筋板，中间完全镂空。那个结果，说实话，比我凭经验画出来的方案轻了40%。

3.2 变密度法原理

变密度法是当前工程中最主流的拓扑优化方法。它的核心假设听起来有点“玄”，但非常实用——把每个有限单元的材料密度当作一个连续变量，取值范围从0到1。

• ρ = 1：该单元有材料，保留
• ρ = 0：该单元无材料，去除
• 0 < ρ < 1：中间密度，灰度单元

但问题来了：密度为0.5的单元，物理上代表什么？半截材料？这显然不合理。所以变密度法引入了一个关键技巧——材料插值模型，最常用的是SIMP模型（Solid Isotropic Material with Penalization）。

E(ρ) = ρ^p × E₀

我曾经犯过一个错误：在早期项目中把惩罚因子设成了1.5，结果跑出来的结果全是灰度单元，根本没法加工。后来才明白，p值至少要取3，才能有效抑制中间密度。

3.3 优化三要素：设计域、约束、目标

任何一个拓扑优化问题，都必须明确三个要素。我习惯把它们称为“拓扑优化的三根柱子”。

3.3.1 设计域

设计域就是你的“画布”。在反射镜支撑中，设计域通常包括：

• 镜体背部的加强筋区域
• 支撑结构（如柔性铰链、背板）的包络空间
• 不允许修改的区域（如镜面光学区、安装接口）要提前设为非设计域

我记得有一次，一个同事把整个镜体都设成了设计域，结果优化结果把镜面也挖掉了——这显然不行。所以，非设计域的设定是第一步，也是最容易出错的一步。

3.3.2 约束

约束是保证结构“能用”的底线。常见的约束包括：

1. 体积约束：保留材料的百分比，比如“只保留30%的材料”
2. 应力约束：最大von Mises应力不超过材料屈服强度
3. 位移约束：镜面变形不超过λ/10（比如632.8nm的1/10）
4. 频率约束：一阶固有频率高于某值（避免共振）
5. 制造约束：最小成员尺寸、拔模方向、对称性等

3.3.3 目标

目标函数是优化器“努力的方向”。在反射镜支撑中，最常见的目标是：

• 最小柔度（即最大刚度）：给定材料用量，让结构最“硬”
• 最小质量：给定刚度要求，让结构最“轻”
• 最小柔度 + 频率约束：兼顾静刚度和动刚度

我个人习惯用最小柔度作为目标，因为它在数学上最稳定，收敛性最好。然后通过调整体积约束来控制质量。

3.4 拓扑优化的知识体系

下面这张图，是我自己总结的拓扑优化入门知识框架。你可以把它当作一张“地图”，后续几讲都会围绕它展开。

3.5 一个简单的入门案例

为了让你更直观地理解，我拿一个简化的反射镜支撑来举例。

问题描述：一个圆形镜体，直径200mm，背面需要设计支撑结构。设计域是镜体背面到安装法兰之间的圆柱体空间。目标是在保留30%材料的前提下，使结构刚度最大。

优化设置：

• 设计域：整个圆柱体（高度20mm）
• 非设计域：镜面区域、安装法兰面
• 约束：体积保留30%
• 目标：最小柔度
• 惩罚因子：p = 3

优化结果通常是一个类似“三脚架”或“六边形蜂窝”的筋板布局。材料集中在从镜面中心到支撑点的传力路径上，其余区域被挖空。

3.6 常见误区与避坑

最后，我总结几个新手容易踩的坑：

1. 设计域给得太小：优化器没有足够的空间去“发挥”，结果跟传统设计差不多。我建议设计域至少比预期结构大20%。
2. 忽略制造约束：拓扑优化结果往往是不规则的曲面，直接加工很困难。一定要加上最小成员尺寸、对称性等制造约束。
3. 目标函数选错：如果目标是“最小质量”，但没加刚度约束，优化器会把所有材料都去掉——结果就是零质量、零刚度，毫无意义。
4. 惩罚因子太小：p值小于2时，灰度单元泛滥，结果无法解释。p=3是工程默认值。