1. 课程导入

拓扑优化解决的是“材料放在哪里最合理”的问题:在给定载荷、边界条件和体积上限下,把有限材料分配到最需要的位置,使结构刚度最大、应力更均匀或频率满足要求。但经典拓扑优化往往只关心数学上的最优解,产出的结构常带有大量悬垂面、细长杆、内部空腔和突变截面。这些形态在减材制造中可以通过多轴加工勉强实现,在铸造中可以通过分型与芯子解决,而在增材制造中却会直接转化为支撑材料、残余应力和后处理成本。

增材制造约束拓扑优化(AM-constrained topology optimization)就是在优化模型中显式引入增材制造工艺约束,使优化结果在获得良好力学性能的同时,具备可打印性。它研究的不是“先优化、后修补”的串行流程,而是把悬垂角、最小特征尺寸、自支撑、连通性、热积累等工艺条件嵌入灵敏度分析或过滤步骤,让优化器在迭代过程中主动避开不可制造的结构。

核心区别可以概括为:传统拓扑优化输出“力学最优形状”,增材制造约束拓扑优化输出“力学最优且可稳定打印的形状”。前者是设计问题,后者是设计—工艺协同问题。

为什么需要工艺约束

以激光粉末床熔融(LPBF)为例,当悬垂面与基板夹角小于约 40°~45° 时,熔池下方缺乏足够支撑,容易出现翘曲、挂渣甚至塌陷。若优化结果中存在大量此类区域,就必须添加支撑。支撑带来三个后果:一是材料与打印时间增加;二是支撑与本体连接处难以去除,影响表面质量;三是支撑区域的热历史与本体不同,可能引入新的残余应力。因此,把悬垂角约束前移到优化阶段,比事后加支撑更经济。

另一个典型约束是最小特征尺寸。LPBF 的可分辨最小杆径通常在 0.2~0.5 mm 量级,但真正能稳定成形的杆径往往更大。如果优化结果中出现大量接近分辨率极限的细杆,打印后可能断裂或尺寸偏差过大。通过在密度场中施加最小长度尺度过滤,可以避免这类“理论可行、实际不可靠”的结构。

增材制造约束拓扑优化:从设计输入到可打印结果 设计输入 • 设计域与体积分数 • 载荷与边界条件 • 材料属性(E、ν、σy) • 打印方向与基板位置 • 工艺窗口(激光功率等) • 目标:刚度/频率/应力 优化循环(含工艺约束) 有限元分析 K·u = f 灵敏度分析 ∂c/∂ρₑ 密度更新 OC / MMA 工艺约束投影 / 过滤 悬垂角 · 最小尺寸 · 自支撑 · 连通性 · 热积累 迭代 输出与验证 • 可打印密度场 • STL / 实体模型 • 悬垂角检查 • 最小杆径检查 • 仿真打印验证 • 力学性能复算 典型应用场景与关注点 航空航天支架 减重与刚度并重 悬垂角 ≤ 45° 最小杆径 ≥ 0.8 mm 避免封闭内腔 热交换器/流道 多物理场耦合 流道连通性约束 壁厚均匀性 粉末排出通道 医疗植入体 多孔结构设计 弹性模量匹配 最小孔隙 ≥ 0.3 mm 生物相容性 轻量化连接件 多工况载荷 疲劳寿命约束 装配接口保留 打印方向优化

典型应用与价值

在航空航天领域,支架类零件常采用“减重优先、刚度约束”的优化目标。以某钛合金支架为例,设计域尺寸约 120 mm × 80 mm × 60 mm,材料为 Ti-6Al-4V,弹性模量约 110 GPa,屈服强度约 880 MPa。体积分数设为 0.3,即最终结构只保留设计域 30% 的材料。若不加悬垂角约束,优化结果中会出现大量与基板夹角约 20° 的斜杆,打印时必须加支撑;加入 45° 悬垂角约束后,斜杆角度被“推”到 45° 以上,支撑量可显著减少,同时刚度仅下降约 5%~8%。这个折中正是工艺约束优化的价值所在。

在热交换器和流道设计中,约束重点从悬垂角转向连通性与粉末排出。LPBF 打印内部流道时,未熔化的粉末会残留在封闭腔体内,必须设计排出通道。拓扑优化若只考虑流体压降和换热效率,可能生成封闭或半封闭的流道网络,导致粉末无法清除。因此,优化模型需要加入连通性约束,确保每个空腔都与外部连通,同时控制最小流道直径,避免堵塞。

医疗植入体则更关注多孔结构的孔径与弹性模量匹配。以 Ti-6Al-4V 髋臼杯为例,骨小梁多孔结构的孔径通常设计在 0.3~0.8 mm,孔隙率 60%~80%,弹性模量目标在 1~5 GPa 之间,以接近人体骨组织。拓扑优化在这里的作用不是追求最大刚度,而是在满足孔径和连通性约束的前提下,让应力传递更均匀,减少应力屏蔽效应。

学习路线

本课程按“基础—约束—实现—验证”的逻辑展开,建议按以下顺序推进:

  1. 拓扑优化基础:理解 SIMP 插值模型、密度场、灵敏度分析和优化准则法(OC),能独立写出二维悬臂梁的最小柔度优化代码。
  2. 增材制造工艺约束建模:掌握悬垂角约束、最小长度尺度过滤、自支撑约束和连通性约束的数学表达,理解它们如何嵌入灵敏度分析。
  3. 数值实现:在 MATLAB 或 Python 中实现带过滤的优化循环,处理灰度单元、棋盘格和网格依赖性。
  4. 打印方向与各向异性:理解打印方向对力学性能的影响,学习如何把方向变量引入优化模型。
  5. 结果验证:用悬垂角检查、最小杆径测量和仿真打印验证优化结果,形成“优化—检查—修正”的闭环。
学习本课程不需要预先掌握增材制造设备操作,但需要具备材料力学、有限元方法和基础编程能力。每章会给出可运行的示例参数和判断依据,建议边学边动手实现。

一个最小示例:悬臂梁的工艺约束对比

考虑一个 80 mm × 40 mm 的二维悬臂梁,左端固定,右端下角施加集中力 F = 100 N,体积分数 0.4,材料弹性模量 200 GPa,泊松比 0.3。采用 80 × 40 网格,单元尺寸 1 mm。传统 SIMP 优化(无工艺约束)会得到一棵“树状”结构,主杆与水平方向夹角约 30°,存在明显悬垂。加入 45° 悬垂角约束后,主杆角度被抬升到 45° 以上,结构略微变重,但支撑面积减少约 40%。若同时加入最小杆径 2 mm 约束,细枝被合并,结构更简洁,刚度下降约 3%。这些数值属于教学示例中的合理假设,用于说明约束之间的折中关系,实际结果会随网格、过滤半径和优化参数变化。

判断优化结果是否可用,可以按以下顺序检查:第一,看密度场是否清晰,是否存在大量灰度单元;第二,检查悬垂角分布,统计小于 45° 的面积占比;第三,测量最小杆径,确认不低于工艺下限;第四,检查连通性,确认无封闭空腔;第五,复算刚度或频率,确认力学性能满足设计要求。只有这五项都通过,优化结果才具备进入打印验证的资格。

接下来的章节将从拓扑优化的数学基础开始,逐步引入增材制造约束,最终形成一套可复现、可验证的工艺约束优化流程。