基于拓扑优化的多层点阵超材料超刚度设计仿真

于 2026-10-25 10:39:50 发布 4 0 0

关键词：拓扑优化；点阵超材料；TPMS；多层策略；能量吸收

一、文章简要介绍

点阵超材料按基本构型分三类：杆系、板式和壳式。杆系容易在节点应力集中，板式刚度高但可调维度少，壳式（如TPMS）光滑连续、抗屈曲，拓扑变化却有限。这篇清华大学与南洋理工大学团队2024年发表在Nature Communications的工作，把三条路线捏在一起：先用多层策略放大设计空间，再交给拓扑优化去搜，优化结果自动收敛成杆-板-壳复合构型，刚度直接摸到理论上限（DOI：10.1038/s41467-024-47089-8）。

图1 多层策略与拓扑优化技术路线

多层策略放大设计空间，拓扑优化自动收敛出杆-板-壳复合构型。

二、拓扑优化建模

单元基于三周期极小曲面生成，Schwarz P面就是cos X加cos Y加cos Z等于零的等值面。多层策略用缩放、变换、混合三种手段做出P-2、P-4、P-5嵌套构型。拓扑优化以最小化整体柔度为目标：设计变量是单元密度，取值0到1，约束包括平衡方程、面积分数上限和厚度上限，用常微分方程驱动的水平集密度法迭代。这套设计域、目标、约束、迭代的流程，和OptiStruct拓扑优化模块走的是同一条路。

基于拓扑优化的多层点阵超材料超刚度设计仿真的图2

图2 各构型的刚度与强度对比

优化后构型的刚度与屈服强度显著抬升，逼近理论上限。

单元尺寸15毫米，基体材料按铁算，弹性模量210吉帕、屈服强度400兆帕，相对密度控制在20%以内。结果分两层看：单层优化里Opt-Neovius的刚度和屈服强度提升约50%，逼近Voigt上限；多层加优化后Opt-P-4同样接近上限。幂律拟合指数从2.27降到1.75，说明变形模式从弯曲主导转成了拉伸主导——材料利用率完全不同。

三、实验验证与关键发现

仿真靠实验兜底。团队用TPU、PA12尼龙、316不锈钢、AlSi10Mg铝合金四种材料打印试件，INSTRON试验机单轴压缩，刚度排序与仿真一致。放大到4×4×4阵列，Opt-P-4的实测刚度达到原始Schwarz P的3倍以上。

$基于拓扑优化的多层点阵超材料超刚度设计仿真的图3$

图3 四种材料打印试件与应力-应变曲线

四种材料打印试件，应力-应变排序与仿真一致。

基于拓扑优化的多层点阵超材料超刚度设计仿真的图4

图4 4×4×4阵列实验与仿真对比

4×4×4阵列实验：Opt-P-4实测刚度约为原始构型3倍。

吸能场景更实用：以每毫秒2毫米的速度动态压缩20毫秒，压缩量三分之二。原始Schwarz P还有大片区域停在弹性段，材料没用满；Opt-P-4几乎整体屈服，归一化吸能到236%，比原始构型高136%。

图5 动态冲击下的变形与吸能对比

动态压缩下Opt-P-4几乎整体屈服，吸能提升136%。

可调性是这套方法的隐藏价值：P-5的Zener比接近1，几乎各向同性；面积分数与刚度呈线性关系，按梯度组装微单元就能得到功能梯度点阵，论文演示了高孔隙人工骨的填充设计。同一套优化框架还能延伸到声学、静电屏蔽和热-流耦合场景。