4. 悬臂梁拓扑优化实战
悬臂梁是拓扑优化最经典的验证模型：左端固定、右端受集中力，设计域为一个矩形板。它的价值不在于结果本身，而在于它把“几何清理—网格划分—优化设置—结果解读—CAD重构”这条完整链路压缩到最小规模，任何一步出错都会在最终构型上留下明显痕迹。本节以一个 200 mm × 50 mm × 5 mm 的钢板悬臂梁为例，右端下缘施加 500 N 向下载荷，材料取钢（E = 210 GPa，ν = 0.3，ρ = 7.85×10⁻⁹ t/mm³），目标是在体积分数不超过 40% 的前提下最大化刚度。以下参数均为教学示例，实际项目需按真实工况替换。
几何清理：先决定“哪些区域不能动”
拓扑优化的设计域必须是干净的实体，不能带小圆角、倒角、装配孔或导入残留的碎面。清理的核心动作有三类：一是删除对刚度贡献极小但会污染网格的细节，例如 R2 以下的圆角；二是把固定端和加载端预留为非设计域，避免优化器把约束点或加载点“挖掉”导致结果不可用；三是统一单位与坐标系，本例约定 X 向右、Y 向上、Z 为板厚方向，载荷沿 −Y 方向，全文与图示保持一致。
在 HyperMesh 中，通常先做几何拓扑检查，合并自由边与重复面，再用 Geom > Quick Edit 修补缺失面。设计域与非设计域建议用不同 component 分开存放，例如 design_domain 与 non_design，后续在 OptiStruct 中通过 DESVAR 与 DVGRID 或直接以属性区分。若几何来自 CAD，务必先做一次 edges > preview equiv，把容差内的自由边合并，否则网格会出现不连续节点，优化迭代中刚度矩阵会异常。
网格划分：单元尺寸决定构型分辨率
拓扑优化的网格不是越细越好，而是要与最小成员尺寸匹配。经验上，若希望优化结果中杆件最小宽度约为 6 mm，则单元边长取 2 mm 左右，即最小成员尺寸约为 3 个单元宽。本例 200 mm × 50 mm 的板，取 2 mm 单元，得到约 100 × 25 = 2500 个壳单元，计算量适中。单元类型用 CQUAD4，避免大量三角形单元造成刚度偏硬和棋盘格。
网格质量要控制：雅可比大于 0.6，翘曲小于 5°，长宽比小于 5。板厚方向若用壳单元，需在属性中给出 PSHELL 的厚度 T = 5 mm；若用实体单元，则至少两层，否则弯曲刚度会被高估。加载点与约束点附近应保证节点落在几何点上，必要时局部加密，但不要突然过渡，否则会在结果中形成虚假的“硬点”。
优化设置：目标、约束与制造工艺
在 OptiStruct 中，拓扑优化的标准卡片组合是：DESVAR 定义设计变量，TOPDES 或 TOPOLOGY 指定拓扑优化，RESPONSE 定义柔度与体积响应，OBJECTIVE 设为最小化柔度，CONSTRAINT 设为体积分数上限 0.4。若使用 HyperMesh 的 Optimization 面板，对应操作是：Topology > design region 选设计域，Responses > compliance 与 volume，Objectives > minimize compliance，Constraints > volume fraction < 0.4。
制造工艺约束是工程可用性的关键。本例加入最小成员尺寸 6 mm 和 X 轴对称约束。最小成员尺寸过小会得到细密枝杈，无法加工；过大则构型粗笨、刚度下降。对称约束能避免优化器因数值噪声产生上下不对称的“伪最优”。若后续要铸造，还应加入拔模方向约束；若为增材制造，则更关注最小成员与悬垂角。
结果解读：密度云图怎么看
求解完成后，OptiStruct 输出 .h3d 结果，在 HyperView 中查看 Element Densities。密度接近 1 的区域保留材料，接近 0 的区域可去除，中间灰度区是优化器尚未完全收敛的过渡带。判读时重点看三件事：一是载荷路径是否连续，从加载点能否沿材料“走”到固定端；二是是否存在孤岛，即与主传力路径不相连的高密度小块，这类区域在重构时应删除；三是上下是否对称，若明显不对称，需检查对称约束是否生效或网格是否偏斜。
本例的典型结果是：上、下缘形成两条近似平行的弦杆，中间大部分被挖空，靠近固定端材料逐渐增厚，靠近加载端出现斜向过渡。这与悬臂梁的弯矩分布一致——固定端弯矩最大，故材料最多；自由端弯矩趋近于零，故材料最少。若结果中出现大面积灰色（密度 0.4～0.6），说明迭代未收敛，可把 MINDENS 从默认 0.01 适当调整，或增加迭代步数到 60 步再观察。
CAD 重构：从密度场到可制造几何
密度云图不能直接用于加工，必须重构为光滑的实体或壳。常用做法是在 HyperMesh 中用 OSSmooth 导出 STL 或 IGES 等值面，等值面阈值一般取 0.3～0.5，本例取 0.4。导出后不要直接交付，而要回到 CAD 中做三步处理：一是用样条重新拟合边界，去掉锯齿；二是把最小圆角补到 R3 以上，满足加工要求；三是把非设计域与优化区域合并为一个连续实体，检查是否有薄壁或尖角。
重构完成后，建议做一次验证：把重构几何重新划分网格，施加同样的边界条件与载荷，比较其柔度与优化结果的柔度。若两者相差在 10% 以内，说明重构基本保留了优化构型；若明显变差，通常是等值面阈值过高导致杆件被削细，或圆角处理过度削弱了传力路径。这一步是很多初学者容易忽略的闭环，也是从“看起来像”走向“真的能用”的关键。