4. 轻量化设计中的拓扑优化原理
拓扑优化是OptiStruct轻量化流程中最具“结构创成”意味的一步：它不要求工程师先给出具体截面尺寸，而是把设计区域当作可增可减的材料场，在给定载荷、边界条件和体积上限下，让算法决定材料应该分布在哪里。理解其原理，关键要抓住四件事：材料如何被参数化（SIMP）、设计变量与响应如何定义、体积约束与目标函数如何共同构成优化问题、以及拔模/挤压/对称等制造工艺约束如何把“数学最优”拉回“可制造”。
一、变密度法SIMP：把“有/无材料”变成连续可导问题
最朴素的拓扑优化是0/1离散问题：每个单元要么保留材料，要么删除。但离散变量不可导，无法用梯度类算法高效求解。SIMP（Solid Isotropic Material with Penalization，固体各向同性材料惩罚模型）的做法是引入一个连续伪密度变量 ρe∈[ρmin, 1]，并让单元刚度按幂次关系缩放：
E_e(ρ_e) = ρ_e^p · E_0
其中：E_0 为实体材料弹性模量
p 为惩罚因子，OptiStruct 常用 p = 3
ρ_min 为防止刚度矩阵奇异的下限，常取 0.01惩罚因子 p 的作用是“惩罚中间密度”：当 p>1 时，ρ=0.5 的单元刚度只有 0.5³=0.125 倍实体刚度，即“半密度”材料性价比极低，算法会自然把密度推向0或1，最终得到接近黑白分明的拓扑。p 太小则灰度单元多、边界模糊；p 太大则容易陷入局部最优、收敛震荡。工程上常采用连续化策略，例如从 p=1 逐步升到 p=3，兼顾稳定与清晰。
需要明确：SIMP 得到的 ρ 是伪密度，不是真实材料密度，它只用于刚度插值。优化结束后必须设定密度阈值（如 ρ≥0.5 视为保留）提取等值面，再重新建模、赋予真实材料属性做尺寸优化与性能验证。
二、设计变量与响应定义
在OptiStruct中，拓扑优化的设计变量通常由 DESVAR 与 TOPOLOGY 卡片定义，核心是 PTYPE 与 MATINIT：前者指定惩罚方式（SIMP或RAMP），后者给出初始密度。设计区域（design domain）与非设计区域（non-design domain）必须提前划分——安装点、密封面、螺栓孔周围通常设为非设计区，保证连接功能不被优化掉。
响应（response）是优化问题的“度量”，常见类型包括：
- 质量/体积响应：
VOLFRAC或MASS，用于约束材料用量； - 柔度响应：
COMPLIANCE，即应变能，越小代表结构越刚； - 位移响应：
DISPLACEMENT，控制关键点变形； - 频率响应：
FREQ，约束一阶模态避免共振； - 应力响应：
STRESS，控制局部强度。
响应定义要与工况绑定。例如白车身扭转刚度工况，可定义扭转角位移响应；悬架控制臂则可定义多个加载方向的柔度响应。响应数量不宜过多，否则约束冲突会导致无可行解。
三、体积约束与目标函数
拓扑优化的标准数学形式可写为：
min C(ρ) = U^T K U （柔度最小化）
s.t. V(ρ)/V_0 ≤ f （体积分数约束）
0 < ρ_min ≤ ρ_e ≤ 1
K U = F （平衡方程）其中 f 是体积分数上限，例如 0.3 表示最多保留原始设计区域的30%材料。目标函数与约束可以互换：若目标是减重，则把质量设为目标、把位移或频率设为约束；若目标是刚度最大化，则把柔度设为目标、体积设为约束。OptiStruct 的 DRESP1 定义响应，DCONSTR 施加约束，DESOBJ 指定目标，三者共同构成优化模型。
体积分数 f 的选择直接影响结果形态。f 过小，结构会退化为细杆或桁架，制造困难；f 过大，优化空间不足，减重收益有限。工程经验是先从 0.3 左右试算，观察材料分布趋势，再结合制造能力调整。
四、制造工艺约束：拔模、挤压与对称
无约束拓扑优化常产生“树枝状”或带倒扣的结构，无法脱模或加工。OptiStruct 提供三类关键工艺约束：
- 拔模约束（DRAW）：指定拔模方向与最小拔模角，保证沿该方向每个单元上方无遮挡，适用于铸造、注塑件；
- 挤压约束（EXTRUDE）：要求材料沿指定方向等截面延伸，适用于铝型材、冲压件；
- 对称约束（SYMMETRY）：强制关于某平面或轴镜像，适用于左右对称的悬架臂、副车架，同时可减少设计变量、提高求解稳定性。
这些约束通过 MANUFACTURING 卡片下的 DRAW、EXTRUDE、SYMMETRY 子卡片设置。需要注意：工艺约束越强，可行设计空间越小，目标值可能变差，但结果可制造性显著提升。实际项目中通常先做无约束拓扑看趋势，再逐步叠加工艺约束做对比。
五、实操要点与结果判读
在OptiStruct中提交拓扑优化，典型流程是：定义设计区域与非设计区域 → 建立工况（约束与载荷）→ 定义响应（柔度、体积分数）→ 设置目标与约束 → 添加制造工艺约束 → 提交求解。求解器默认采用MMA或OC准则迭代，通常30～60步收敛。判断收敛可看目标函数变化率是否小于0.1%，以及密度场是否趋于稳定。
结果判读时，密度云图中红色代表高密度（保留），蓝色代表低密度（可去除）。但不要直接按云图建模：应设定阈值（如0.5）提取光滑等值面，再根据工艺可行性做人工修整——补圆角、加筋、保证最小壁厚。修整后的模型必须重新做有限元分析，验证刚度、模态、应力是否满足原约束。若性能下降超过可接受范围，需回到优化阶段调整体积分数或工艺约束重新计算。