Rapid Kirigami Simulation using the Bar & Hinge Approach
Raj Pradip Khawale, Elaheh Mehdizadeh, John Brigham, Evgueni T. Filipov
cond-mat.other, cs.CE
2026-08-14
把剪纸离散成杆、弯曲铰和三点扭转弹簧。对照Abaqus开口面积差3.8%–6.8%，对照Mylar实验泊松比和面外位移差不到4%；10×10胞元从81645秒降到166秒。
剪纸（kirigami）在薄板上切开缝，平面板就能做大变形、可展开结构和可调刚度超材料。工程上要扫切缝图案，几乎只能上壳单元有限元。壳单元能给出应力场，但薄板大变形、屈曲、剪切闭锁叠在一起，网格一细就慢，还经常不收敛。解析模型快，却绑在特定胞元上。近年也有 GAN、深度学习当代理，训练数据从哪来、切缝一换还能不能用，都是硬约束。
缺的是一个能接任意切缝、只追整体变形和反力、算得动参数扫描的降阶模型。
密歇根大学与匹兹堡大学把折纸社区用了十多年的 bar and hinge（杆-铰）框架接到剪纸网格上。薄板不再当连续壳，拆成三类弹簧：
节点只有平动：3D 全局刚度 \(3n \times 3n\)，2D 是 \(2n \times 2n\)。材料按线弹性，几何非线性用位移控制的 Newton–Raphson，步长自适应。薄板屈曲对初始缺陷敏感，预处理会在 z 向加幅值约 \(10^{-2}\) 倍特征尺寸的模态扰动。网格走 MATLAB PDE Toolbox，按最大边长 \(L{\max}\) 控制，切口附近自动加密。
对照跑在 Apple M2、64 GB 内存上。四组周期性胞元横向拉伸 20%（切缝最多的 10×10 只拉到应变 0.17）：
| 胞元 | 杆铰节点 / 秒 | Abaqus 节点 / 秒 | 开口面积差 |
| 2×5 | 7542 / 89.4 | 67325 / 1144 | 4% |
| 5×10 | 11581 / 102.1 | 69453 / 1835 | 4.5% |
| 5×15 | 14068 / 185.9 | 71929 / 4352 | 3.8% |
| 10×10 | 13704 / 165.9 | 75899 / 81645 | 6.8% |
简单几何大约快 13–23 倍；10×10 那组 Abaqus 跑了 22.7 小时，杆铰 166 秒，大约 492 倍。论文对外口径写「至少十倍」，把收敛顺利的算例也算进去了。杆铰整体略刚，硬结构上更明显。
实验端：单切口圆板对照 Chaudhary 等人的拉伸屈曲；彩虹纹和蜂窝纹 Mylar 激光切割，Mark-10 拉到应变 0.5。泊松比、最大面外位移、开口面积分别差不到 4%、4%、5%。力-应变在小应变贴得紧，大应变再次偏刚。非周期、相交切口在 12% 拉伸下也能出整体构型，没有另给误差表。
四个用得上的场景：75 组拉丁超立方切缝扫渗透率，M2 Max 上约 5 小时，渗透率大致随开口面积三次方走；三侧面爬行器上三角、圆、梯形切缝的各向异性摩擦，三角切最能拉开方向差；方-菱非互易超材料（不对称角 \(\theta=\pi/16\)，Mylar \(E=3.5\,\mathrm{GPa}\)，边长 6.75 mm），\(F<10\,\mathrm{N}\) 就看得出 \(u{21}-u{12}\)；可极化超材料扫扭转角，单次约 20 秒。
做剪纸可展开结构、软体爬行器、建筑表皮开口、面内超材料的人，如果目标是整体变形和反力而不是切口尖端应力，这套离散弹簧能把「扫几百个几何」从有限元农场变成单机任务。MATLAB 代码开源。
它不靠神经网络代理，换切缝不用重训；也不是新的连续介质理论，就是把折纸杆铰接到剪纸三角网上。和 GAN 全场预测比，它慢一些，但能吐刚度矩阵和内力，还能做特征值和屈曲。\(Fs\)、\(FB\) 是全局常数，换材料、厚度和网格密度时要自己核一遍。
论文自己列了：局部应力集中和边缘效应抓不住，目前不能估断裂或屈服；剪切刚度系统性偏高（杆网用轴向变形近似剪切，一个全局 \(Fs\) 抹不平）；有网格收敛问题；只有 MATLAB 脚本，没有封装好的界面；材料非线性、塑性、率相关都没做。
对照并不完全公平。Abaqus 节点数是杆铰的 5–9 倍，加速里有一部分来自更粗的离散。10×10 拉到近 500 倍，很大一块是壳单元大变形收敛差，不是渐近复杂度。杆铰几乎处处略刚，高应变更明显，拿它估峰值载荷会偏危险侧。爬行器摩擦用的是经验式 \(\muf=\mu0+Kf\cdot h\cdot\rho\cdot Sf\)，不是接触仿真。\(Fs\)、\(FB\) 只在标定集上调过，泛化范围没有单独测。