拓扑优化多向软气动执行器 TO-mdiSPAs：仿真实现三臂夹持器多方向弯曲

如何用拓扑优化系统地设计出能在三维空间任意方向弯曲的软气动执行器？

现有软气动执行器大多依赖设计师经验和直觉，缺乏系统化设计方法；大多数执行器只能实现平面弯曲，难以满足三维复杂任务；压力载荷随设计变化，传统拓扑优化难以处理这类设计相关载荷。

此前软气动执行器设计主要基于启发式原则；少数拓扑优化研究要么只针对单向或双向弯曲，要么未同时考虑设计相关压力载荷和制造误差鲁棒性。

采用密度法拓扑优化，引入三场（blueprint、dilated、eroded）鲁棒公式以容忍制造误差；使用带排水项的达西定律建模设计相关气压载荷；目标为最小化blueprint和eroded设计的最大输出位移（即最大化最差输出），并受体积分数和应变能约束；用MMA更新设计变量。设计域为单个气腔对称半截面，离散为80×320个双线性四边形单元，体积分数上限0.20，应变能约束80%，滤波半径7.6。优化得到非传统气腔几何，将其拉伸为三维单元并沿四个方向（前、后、左、右）各排列10个气腔构成臂段，再用三个臂组成夹持器；仿真在Abaqus中进行，采用Ogden超弹性材料模型（μ₁=1.57, α=2, D₁=5×10⁻⁵）。

结果为纯数值仿真，在Abaqus中对三臂夹持器模型进行，未与基线设计定量对比。等压工况（表I）压力为0.150 MPa时，臂1（前+左通道）位移 u_x=1.948 mm, u_y=-15.667 mm, u_z=15.658 mm；臂2（后通道）位移 u_x=0.994 mm, u_y=15.820 mm, u_z=0.004 mm；臂3（前+右通道）位移 u_x=1.948 mm, u_y=-15.667 mm, u_z=-15.658 mm。臂1尾端与正x、y、z轴夹角分别为84.97°、134.80°、45.24°；臂2尾端与正y、x、z轴夹角分别为3.58°、86.40°、89.99°。不等压工况（表II）中，臂1施加前0.050 MPa、左0.075 MPa，位移 u_x=0.009 mm, u_y=-3.114 mm, u_z=5.212 mm；臂2施加前0.100 MPa、右0.150 MPa，位移 u_x=1.131 mm, u_y=-7.809 mm, u_z=-15.751 mm；臂3施加后0.080 MPa、右0.011 MPa，位移 u_x=0.031 mm, u_y=5.676 mm, u_z=-0.630 mm。弯曲方向与压力合力方向接近，证明可通过调节通道和压力大小控制空间弯曲方向。论文只在单一仿真条件下展示，未报告重复试验或统计误差。论文正文中提及等压工况压力为0.015 MPa，但表I为0.150 MPa，存在矛盾，本简报以表I为准。

作者承认尚未制造实物，仅数值仿真，缺少实验验证；未与传统执行器设计进行定量对比；仅使用单一材料模型和单一仿真条件；论文未报告代码开源或制造细节；压力值在正文与表格中不一致（0.015 vs 0.150 MPa），可能影响可重复性。