Nonlinear Mechanics and Predictable Bifurcation of Multi-Cell Kresling Origami Chains

本文通过系统分析随层数增加的平衡分支，建立了一个用于预测多胞 Kresling 折纸链非线性力学与分叉行为的预测框架，最终通过对临界点的几何控制，实现了可编程机械超材料的反向设计。

要点

想象一个由折叠纸张制成的长长的、类似手风琴般的管状物，类似于一种名为 Kresling 的日本折纸图案。当你向下按压这个管子的顶部时，它不仅仅是变短，还会发生扭转。这篇论文探讨了将许多个这样的“纸质单元”堆叠在一起组成一条长链时，它们会发生什么，以及在被挤压时如何表现。

以下是这篇论文的叙述，通过简单的概念进行了拆解：

1. 基础构建块：一个扭转的纸质单元

把单个 Kresling 单元想象成一个由三角形组成的微型中空圆柱体。它有一个特殊的属性：如果你向下按压它，它就会想要扭转。

• 形状至关重要： 论文表明，单个单元的行为高度取决于它的形状。具体来说，它取决于初始形状的“扭曲程度”（折痕的角度）以及它的高度与宽度的比例。
• 四种性格类型： 根据这些形状，研究人员发现单个单元具有四种不同的“性格”（或区域）：
  1. 一心一意型（One-Track Mind）： 它只有一个稳定的形状。如果你按压它，它只会平滑地缩短。
  2. 分裂人格（非对称型）： 它可以处于两种不同的稳定形状，但这两个形状并不是镜像对称的。
  3. 分裂人格（对称型）： 它可以处于两个镜像对称的稳定形状，其中包含一个中间的“悬浮”状态，在该状态下它感觉不到应力。
  4. 富有弹性型（Stretchy One）： 它主要倾向于保持高耸，但也可以瞬间切换到一种拉伸后的形状（不过论文主要关注挤压而非拉伸）。

2. 连锁反应：堆叠单元

研究人员接着问道：“如果我们把两个、三个甚至 n 个这样的单元堆叠在一起，会发生什么？”

想象一叠这样的纸杯。当你向下按压顶部时：

• 双单元堆叠： 如果有两个完全相同的单元，它们可能会表现得截然不同。其中一个可能完全塌陷，而另一个保持挺立；或者它们两者同时塌陷。论文精确地描绘了它们何时会步调一致行动，以及何时会“各行其道”表现出差异。
• 三单元堆叠： 当有三个单元时，情况变得更加复杂。它们可以分成不同的组（例如：两个塌陷，一个挺立；或者三个都表现得各不相同）。研究人员发现，随着增加更多的单元，可能出现的“跳变”时刻也会增加，从而创造出一种稳定性与不稳定性交织的复杂舞蹈。

3. “跳变”与“切换”

这篇论文对 分叉（Bifurcation） 非常感兴趣。用通俗的话说，这就像是在路口遇到了分叉路。

• 突跳（The Snap-Through）： 当你向下按压时，链条会达到一个点，必须选择一条路径。
• 突跳现象： 有时链条是稳定的，但随后只需再加一点点压力，它就会突然“跳变”进入一种新形状。这就像按下易拉罐拉环，它会先承受住压力，然后突然翻转过去。
• 研究人员发现，在链条中，这些跳变并不是同时发生的。它们是按顺序发生的：一个单元跳变，接着下一个，再下一个。这创造了一个能量吸收的“阶梯”过程，这对于需要吸收冲击力的应用（如汽车的溃缩区，尽管论文并未明确提出此应用，但描述了其力学机制）非常有用。

4. 魔术技巧：预测未来

研究这些链条最难的部分在于，随着单元数量的增加，数学计算变得极其混乱。这就像试图预测单片叶子在风暴中的路径，然后又试图预测整片森林的叶子如何随风起舞。

研究人员开发了一种 广义策略（数学上的魔术技巧）：

• 他们意识到，即使是在一个拥有 100 个单元的长链中，这些单元在任何给定时刻也只能以有限数量的“状态”（形状）存在。
• 与其追踪每一个单独的单元，不如对它们进行分组。例如，他们可能会说：“好的，有 4 个单元处于状态 A，有 1 个单元处于状态 B。”
• 通过这种方式，他们仅通过观察单个单元的行为，就能预测整个庞大链条的行为。他们发现，“跳变”点会以非常规则的间隔发生，就像梯子的台阶一样。

5. 大局观：利用不稳定性进行设计

通常情况下，工程师试图制造那些不会摇晃或突跳的物体。这篇论文颠覆了这个想法。它表明我们可以 设计不稳定性。

通过仔细选择折痕的角度和尺寸（几何结构），我们可以精确地告诉链条何时跳变、跳变多少次以及最终会变成什么形状。

• 逆向设计（Inverse Design）： 与其先造出一个链条再观察它的表现，不如你可以直接说：“我想要一个能在特定压力下跳变三次的链条”，然后数学会告诉你该如何建造它。

总结

这篇论文是一份关于复杂的、扭转的、会跳变的折纸链条的地图。它告诉我们：

1. 形状决定行为： 折痕角度的微小变化会导致链条运动方式的巨大变化。
2. 堆叠产生复杂性： 将它们组合在一起会产生新的跳变和状态切换方式。
3. 一切皆可预测： 即使对于非常长的链条，我们也可以使用一种简化的数学技巧来预测跳变发生的精确位置，从而设计出具有特定、可编程行为的结构。

作者们基本上将一个混乱的、扭转的纸质玩具，变成了一台可预测的、可编程的机器。

技术摘要

技术摘要：多单元 Kresling 折纸链的非线性力学与可预测分叉行为

问题陈述
具有轴向-扭转耦合特性的超材料（如源自 Kresling 折纸图案的结构）在可编程功能（包括多稳态、能量吸收和可调刚度）方面具有显著潜力。然而，这类结构工程设计的核心挑战在于理解其复杂的非线性力学行为，特别是其平衡分支和分叉图。随着 $n$ 层链中组成单元数量的增加，系统会在连续的不稳定性（包括分叉点分叉和极限点不稳定性）下，展现出延伸至后临界状态的复杂平衡路径。微小的几何变化可能导致显著的后屈曲行为差异，使得对可控后临界响应进行逆向设计变得困难。现有模型往往难以追踪大规模系统中相关的平衡路径和稳定性切换，特别是在从单单元扩展到多层组装时，对称性破缺和层间相互作用引入了额外的自由度（DOFs）。

方法论
作者开发了一个通用的框架来模拟 Kresling 折纸链的非线性力学，实现了从单单元到 $n$ 层组装体的演进。

• 建模方法： 研究采用基于桁架的描述方式，将折痕线理想化为可变形的轴向载荷承载元件（垂直和对角桁架），而将三角形面板视为刚体。系统由总势能控制，该能量包含来自桁架伸长/收缩的弹性应变能、总链高的约束项以及强制正高度层的惩罚项。
• 控制方程： 通过将总链高作为连续参数，利用变分原理处理代数系统，从而推导出平衡路径。研究使用连续化与分叉分析软件（AUTO 07P）求解该系统。通过投影后的势能 Hessian 矩阵的正定性来评估稳定性。
• 泛化策略： 为了克服大规模 $n$ 层系统中由于分叉点附近雅可比矩阵奇异化带来的数值挑战，作者提出了一种简化策略。该方法将共享相同平衡状态的单元格进行分组，将 $n$ 层问题简化为等效的具有最多五个不同状态的系统（基于单个双稳态单元的本构行为）。这使得高效构建任意层数链的分叉图成为可能。

主要贡献与结果

1. 单单元本构映射： 研究基于初始旋转角 ($\delta$) 和高度-半径比 ($h^{(0)}/R^{(0)}$) 建立了单 Kresling 单元的相图。该图识别了四个不同的区域（I–IV），其特征在于不同的稳定性行为：

   • 区域 I： 单稳态（未变形状态）。
   • 区域 II： 非对称双稳态（未变形态与近乎塌陷态）。
   • 区域 III： 具有中间态的双稳态（未变形态与中间态）。
   • 区域 IV： 单稳态（未变形态，可能存在一个伸展态）。
     分析表明，随着多边形边数 ($N$) 的增加，相边界趋于稳定，对 $N$ 的敏感度降低。

2. 两层与三层链的分叉分析：

   • 两层链： 两个单元的组装引入了对称性破缺分叉。由区域 I 和 IV 单元构建的链在链级层面保持单稳态，但表现出超临界或亚临界分叉。由区域 II 和 III 单元构建的链获得了额外的稳定状态，从而产生三个稳定的平衡状态。分叉图揭示了取决于单元格所属区域的超临界与亚临界分支、合并点以及跳跃（snap-through）不稳定性等截然不同的模式。
   • 三层链： 扩展到三层增加了复杂性，引入了二次和三次平衡分支。系统展现出更丰富的分叉点集（B1–B11），对应于单元格划分到不同变形状态（例如，2 个单元处于一种状态，1 个单元处于另一种状态）的各种情况。

3. $n$ 层链的广义解法：

   • 作者证明，对于由区域 III 单元组装而成的链，其分叉结构由两个截然不同的阶段（阶段 1 和阶段 2）组成。
   • 可预测轨迹： 一个关键发现是，二次分叉点和合并点沿分叉图中五条不同的线性轨迹排列。这些点的分布呈现出由单单元本构响应的局部极值所决定的均匀水平间距。
   • 解析表达式： 研究推导了这些轨迹以及分叉点间距 ($\Delta \tilde{u}$) 的解析表达式。这些表达式依赖于单单元的临界应变值（$u_{Fmax}, u_{Fmin}$ 等）以及层数 ($n$)，从而实现了无需求解整个链的完整控制方程即可预测分叉点。

意义
本文声称为具有可编程响应的构筑化机械超材料的逆向设计与优化提供了一个基础框架。通过建立几何设计变量（多边形数量、初始扭转角、长径比）与所得平衡路径之间的直接联系，该工作实现了对多层超结构的预测性构建。所提出的泛化策略能够系统地追踪复杂的平衡路径，并识别控制后临界行为的关键点。这种能力被认为是定制机械响应（如多稳态、能量吸收和部署特性）的核心，使 Kresling 类结构能够超越简单的展开与塌陷，实现复杂的、可编程的功能。