Tensional wrinkling of thin elastic sheets with two circular holes

技术摘要：含两个圆孔的薄弹性片在张力作用下的起皱

问题陈述
本研究探讨了含两个圆孔的薄线性弹性片中由张力诱导的起皱现象。虽然经典的拉梅问题（无限大薄片中单个圆孔在径向张力作用下）已为人熟知，但多个孔的存在引入了几何对称性破缺和新的长度尺度。作者解决了“双极拉梅问题”，即一片薄板在两个半径为 $R_{in}$ 的圆孔边界处承受均匀法向张力 $T_{in}$，而在无穷远处承受各向同性张力 $T_{out}$。主要目标是确定起皱开始的临界张力比 $\gamma = T_{in}/T_{out}$，分析由此产生的起皱的空间范围和取向，并理解两个孔之间的相互作用如何相较于单孔情况改变不稳定性。

方法论
本研究结合了分析建模、渐近分析和实验验证：

1. 分析建模：作者采用双极坐标 $(\tau, \xi)$ 对两个孔的几何形状进行映射。他们利用艾雷应力函数 $\psi$ 推导了屈曲前（线性弹性）状态的精确解。通过满足双调和方程和边界条件，他们获得了应力分量的级数解。
2. 稳定性分析：起皱的 onset 通过识别屈曲前状态中首次出现的主压应力（$\sigma_2 < 0$）来确定。临界张力比 $\gamma_{cr}$ 被确定为归一化孔间距 $\zeta = d/R_{in}$ 的函数。
3. 渐近分析：
   • 大间距 ($\zeta \gg 1$)：作者进行 $1/\zeta$ 的幂级数展开，推导 $\gamma_{cr}$ 的渐近公式。他们进一步利用帕德近似（Padé approximant）将该近似的适用范围扩展到较小的 $\zeta$。
   • 小间距 ($\zeta \ll 1$)：将 $\zeta^{1/2}$ 的幂级数展开拟合到数值解，以描述孔接近接触时的行为。
   • 起皱范围：通过识别压应力的零等值线来分析近阈值起皱。利用张力场理论和摄动分析研究远阈值行为，以确定起皱区域的空间范围。
4. 实验验证：实验使用旋涂聚苯乙烯（PS）薄膜（厚度 $\approx 166$ nm）漂浮在水 - 气界面上进行。使用朗缪尔槽通过改变表面活性剂浓度来控制外表面张力（$T_{out}$），而内张力（$T_{in}$）保持在纯水的张力不变。这使得能够精确调节 $\gamma$。通过视频成像记录了不同孔间距（$\zeta = 0.5, 1, 1.5, 3$）下的起皱阈值。

主要贡献与结果

• 临界阈值 ($\gamma_{cr}$)：本研究提供了屈曲前应力场的精确解析表达式，并推导了作为孔间距 $\zeta$ 函数的临界张力比 $\gamma_{cr}$。
  • 临界比是 $\zeta$ 的单调递增函数。
  • 当 $\zeta \to \infty$（孔相距很远）时，$\gamma_{cr} \to 2$，恢复了经典的拉梅结果。
  • 当 $\zeta \to 0$（孔接触）时，$\gamma_{cr} \to 1$。
  • 与单孔情况相比，第二个孔的存在始终降低了不稳定性阈值，促进了更早的起皱。
• 渐近近似：作者提供了显式的有理近似（通过帕德近似）和幂级数拟合，能够准确预测整个 $\zeta$ 范围内的 $\gamma_{cr}$， bridging 了单孔极限和接触孔极限之间的差距。
• 起皱成核与形态：
  • 成核：起皱首先在孔边界上最靠近另一个孔的点成核（在双极坐标中具体对应 $\xi = \pi$，即孔的“内侧”）。
  • 演化：随着 $\gamma$ 增加，起皱区域扩大。对于中等间距，起皱最初出现在孔的相对侧，并最终在中心区域合并。对于大间距，局部行为类似于单孔情况，但全局对称性被打破。
  • 远阈值：在远阈值机制下，起皱区域延伸得比近阈值分析预测的要远得多。最大起皱范围出现在连接两个孔中心的轴线上，确定了这是应力传递的主方向。
• 实验一致性：不同 $\zeta$ 值下的起皱阈值实验测量结果与理论预测显示出极好的一致性（误差范围在 $\pm 5\%$ 以内）。实验证实，薄膜厚度、液体密度和重力对阈值没有显著影响，验证了模型对张力和几何形状的关注。

意义与主张
该论文主张，几何对称性破缺，特别是引入第二个孔，从根本上改变了张力诱导起皱的力学机制。主要发现包括：

1. 应力传递：与单孔拉梅问题类似，起皱的存在允许应力信号在更长的距离上传递（随距离倒数衰减，而非距离平方衰减）。然而，双极几何结构引入了沿连接孔轴线方向的优先传递方向。
2. 力学传感类比：作者提出，这些发现为细胞力学传感提供了力学类比。正如单个细胞可以产生在远处可检测的应力场一样，两个“细胞”（由孔代表）之间的相互作用创造了一个复杂的应力景观，其中力学不稳定性（起皱）放大并引导力的传递。这为理解多个细胞如何通过 compliant 基底进行力学相互作用提供了一个初步模型。
3. 几何控制：本研究证明了起皱的模式和 onset 可以通过几何参数（孔的大小和间距）进行精确控制，为设计具有特定应力重分布特性的薄膜系统提供了一种机制。

该工作被呈现为理解薄片中多孔相互作用的基础步骤，超越了理想化的单缺陷模型，转向具有对称性破缺的更真实系统。