1. 联合优化问题定义
结构优化在工程实践中通常被拆成三个层次：拓扑优化决定材料“在哪里”，尺寸优化决定构件“多粗多厚”，形状优化决定边界“长成什么样”。三者单独使用都能取得效果，但真实工程需求往往要求它们同时满足。本节先建立三者协同的数学表述，再说明设计变量之间的耦合关系，最后给出典型工程需求的判断依据。
一、三个层次的数学表述
设设计域为 Ω，材料分布用密度场 ρ(x) 表示，ρ=1 为实体、ρ=0 为空。以最小柔顺度（最大刚度）为例，单一拓扑优化的标准形式为：
min C(ρ) = fᵀu
s.t. K(ρ)u = f
V(ρ) = ∫ρ dΩ ≤ V*
0 ≤ ρ ≤ 1其中 K 为刚度矩阵，u 为位移，f 为载荷，V* 为体积上限。尺寸优化把 ρ 换成显式截面参数，例如杆件截面积 Aₑ 或板厚 tₑ，约束仍是体积或质量：
min C(A) = fᵀu
s.t. K(A)u = f
Σ AₑLₑ ≤ V*
A_min ≤ Aₑ ≤ A_max形状优化则把边界节点坐标 x_b 作为变量，网格随边界移动而变形：
min C(x_b) = fᵀu
s.t. K(x_b)u = f
g(x_b) ≤ 0 (几何/应力约束)三者协同的联合优化，就是把上述变量合并为一个设计向量：
min C(ρ, A, x_b) = fᵀu
s.t. K(ρ, A, x_b)u = f
V(ρ, A, x_b) ≤ V*
ρ ∈ {0,1} 或 [0,1]
A_min ≤ A ≤ A_max
x_b ∈ X_feasible关键点在于：目标函数和约束同时依赖三类变量，刚度矩阵 K 是它们的复合函数，因此不能简单地把三个子问题串行求解后叠加。
二、设计变量的耦合关系
三类变量并非独立。拓扑变量 ρ 决定哪些区域存在材料，尺寸变量 A 只在 ρ 非零的区域有意义；形状变量 x_b 改变边界位置，又会反过来影响 ρ 的分布和 A 的取值。这种耦合可以用“变量作用域嵌套”来理解：拓扑给出可行域，尺寸在可行域内分配截面，形状在可行域边界上做微调。
耦合还体现在灵敏度上。以柔顺度对某变量的导数为例，拓扑灵敏度为 ∂C/∂ρ = -uᵀ(∂K/∂ρ)u，尺寸灵敏度为 ∂C/∂A = -uᵀ(∂K/∂A)u，形状灵敏度为 ∂C/∂x_b = -uᵀ(∂K/∂x_b)u。三者形式一致，但 ∂K/∂ρ、∂K/∂A、∂K/∂x_b 的物理含义不同：前者是材料有无，中者是截面刚度贡献，后者是网格几何映射。若忽略这种差异，直接套用同一套过滤和更新策略，容易出现棋盘格、网格畸变或尺寸振荡。
三、典型工程需求与判断依据
什么时候必须做联合优化，而不是只做其中一层？可以从三个信号判断。
信号一：载荷路径不明确。如果初始设计域内没有明显的传力路径，仅靠尺寸优化无法决定材料该放在哪里，必须先做拓扑优化。例如卫星支架在多个工况下受力方向不同，材料布局本身是未知量。
信号二：截面与布局相互制约。当拓扑给出的杆件布局合理，但各杆截面差异很大，且截面变化会改变整体刚度分布时，尺寸变量不能事后单独调整。例如桁架桥的杆件截面与节点位置共同决定弯矩分布，需要与形状变量同步更新。
信号三：边界形状影响局部应力。如果结构存在应力集中或疲劳敏感区，边界圆角、过渡曲线等形状变量会显著改变峰值应力，而拓扑和尺寸无法替代这种几何光顺。例如涡轮盘榫槽边缘，形状微调可降低应力集中系数 10%～20%（示例量级，实际依工况而定）。
一个可操作的判断流程是：先看设计域是否已知。若未知，必须引入拓扑变量；再看截面是否可调。若可调且影响全局刚度，引入尺寸变量；最后看边界是否允许移动。若允许且应力/疲劳是控制指标，引入形状变量。三者同时满足，才构成联合优化问题。
四、一个简化的教学示例
考虑一个 200 mm × 100 mm 的悬臂梁设计域，左端固定，右端下角受集中力 F=1000 N。材料弹性模量 E=210 GPa，泊松比 0.3，体积上限为设计域的 40%。
若只做拓扑优化，得到的是黑白密度分布，杆件边界呈锯齿状，无法直接加工。若只做尺寸优化，需要预先给定杆件布局，布局不合理时优化空间有限。若只做形状优化，材料分布不变，只能微调边界。
联合优化的做法是：拓扑变量控制材料分布，尺寸变量为每根杆件分配截面积，形状变量调整杆件连接处的过渡曲线。迭代过程中，体积约束始终保持在 40%，目标柔顺度从初始设计的 2.8 J 逐步下降到约 1.1 J（示例数值，用于说明下降趋势，非实测结果）。判断收敛的依据是：连续 5 次迭代柔顺度变化小于 0.5%，且体积约束满足、无棋盘格、边界光顺。
这个示例说明：联合优化的价值不在于某一层做得更好，而在于三层变量在同一个约束下互相让步，最终得到一个既省材料、又满足刚度、还能加工的设计。