几何弹性 · 生长与形状

非欧薄板

一张按弯曲几何生长的圆盘，为什么厚时保持平坦，薄时却变成穹顶？ 生长可以要求薄片具有平面做不到的弯曲几何，比如让外圈长得比平面所需的短。厚圆盘保持平坦，用拉伸承受这种不匹配；薄于某个临界厚度后，弯曲更省能量，圆盘屈曲成几乎没有拉伸的穹顶。越薄，穹顶越接近球冠。

生长可以给薄片一种平面做不到的目标几何，这里是球面的几何。让圆盘生长：厚圆盘保持平坦、靠拉伸承受；薄圆盘屈曲成几乎没有拉伸的穹顶。屈曲厚度、各项能量和穹顶高度都与无拉伸球冠和文献给出的临界值对照。

低于此厚度，平盘就会屈曲 $h_{\mathrm{B}}=0.299$

长成后能量中拉伸所占的比例 $7.24\,\%$

穹顶高度与无拉伸球冠之比 $0.999$

这张圆盘离开平面时的生长量

K

$K_{\mathrm{B}}=0.0186$

模型假设

圆盘是 Efrati、Sharon 与 Kupferman（J. Mech. Phys. Solids，2009）的轴对称非欧薄板：测地半径从 0.1 到 1.1 的圆环，目标度规具有恒定高斯曲率 K，采用他们式 4.2 的约化能量，并像他们的图 4.1 一样把前置系数取为 1。径向分成 80 个单元，生长的每一帧都是用牛顿法求出的平衡态。只允许旋转对称的形状，因此褶皱、凹坑和负曲率的马鞍形都不在模型范围内；生长是给定的，而不是算出来的。探针是作用在内缘高度上的硬弹簧；松开后，圆盘按过阻尼梯度流弛豫，这是一条通往最近平衡态的路径，而不是带惯性的真实薄片的运动。

橙色是拉伸、蓝色是弯曲、绿色是两者之和，都在圆盘模型上算出；米色虚线是平盘和无拉伸球冠。

更多读数与保存的记录

平盘的拉伸能

0.0143

无拉伸球冠的弯曲能

0.541

两者相等时的交叉厚度

0.163

屈曲厚度除以 K 的平方根

0.299

长成后圆盘的弹性能

2.69\times10^{-5}

长成后穹顶的高度

0.541

无拉伸球冠的高度

0.541

径向单元数

80

JSON 记录包含参数、生长过程和厚度扫描中的每个平衡态、屈曲临界值，并单独列出平盘与无拉伸极限以及文献给出的临界值。

物理教程

生长怎样让平的圆盘弯起来？

背景叶片、花瓣和不均匀溶胀的凝胶片，被生长赋予的是一种几何，而不是一个现成的形状。如果这种几何是弯曲的，平的薄片就不可能具有它，薄片要么拉伸，要么弯离平面。

为什么值得理解无论厚薄，圆盘的每个小单元被要求的长度都一样。可是厚圆盘保持平坦，薄圆盘变成穹顶：厚度决定的只是拉伸和弯曲这两种代价哪一种更便宜。

先抓重点

核心问题: 一张按弯曲几何生长的圆盘，为什么厚时保持平坦，薄时却变成穹顶？
一句话直觉: 拉伸的能量与厚度成正比，弯曲的能量与厚度的三次方成正比。厚圆盘在平面内拉伸来承受弯曲几何；一旦薄于临界厚度，弯成几乎等距的穹顶更省能量，而且越薄，剩下的拉伸能比弯曲能下降得快得多。

核心数学模型

目标几何

d\bar{s}^{2}=dr^{2}+\Phi(r)^{2}\,d\theta^{2},\qquad \Phi(r)=\frac{\sin(\sqrt{K}\,r)}{\sqrt{K}},\qquad K=-\frac{\Phi''}{\Phi}

沿半径的长度不变，但测地半径 $r$ 处的圆环周长应为 $2\pi\Phi(r)$ ，短于 $2\pi r$ 。这正是半径为 $1/\sqrt{K}$ 的球面的度规（Efrati 等，第 4 节）。

拉伸加弯曲

E=\int_{r_{0}}^{r_{1}}\bigl(h\,w_{S}+h^{3}\,w_{B}\bigr)\,\Phi\,dr,\qquad w_{S}=\epsilon_{1}^{2}+\epsilon_{2}^{2}+2\nu\,\epsilon_{1}\epsilon_{2}

Efrati 等的式 4.2，适用于半径为 $\phi(r)$ 、斜率为 $\zeta$ 的旋转面。应变 $\epsilon_{1}=\phi'^{2}(1+\zeta^{2})-1$ 和 $\epsilon_{2}=\phi^{2}/\Phi^{2}-1$ 衡量与目标长度的偏差； $w_{B}$ 是弯曲能密度。这里的圆盘是圆环 $0.1\le r\le 1.1$ 。

无拉伸球冠

\phi=\Phi,\qquad H_{0}=\frac{\cos\sqrt{K}r_{0}-\cos\sqrt{K}r_{1}}{\sqrt{K}},\qquad E_{B0}=\tfrac{2}{3}(1+\nu)\bigl(\cos\sqrt{K}r_{0}-\cos\sqrt{K}r_{1}\bigr)

球面的一个环带恰好具有目标长度，所以拉伸为零。它的高度是 $H_{0}$ ；由于其上处处弯曲能密度为 $\tfrac{2}{3}(1+\nu)K$ ，弯曲能为 $E_{B0}$ 。当 $K=1$ 、 $\nu=0.5$ 时， $H_{0}=0.541$ ， $E_{B0}=0.541$ 。

两个厚度

h_{C}=\sqrt{E_{S}^{\infty}/E_{B0}}\approx 0.163,\qquad h_{B}\approx 0.299

令平盘的拉伸能 $h\,E_{S}^{\infty}$ 等于球冠的弯曲能 $h^{3}E_{B0}$ ，得到交叉厚度 $h_{C}$ 。平盘的线性稳定性给出屈曲厚度 $h_{B}$ ；Efrati 等给出约 0.3，并指出 $h_{C}$ 更小。

常见难点

这里薄并不意味着软

容易误解薄圆盘屈曲是因为它太软，撑不住形状。

正确理解每种厚度被要求的弯曲长度都一样。厚度只决定弯曲相对拉伸的价格，这个比值随厚度平方增长，所以薄于临界值时，圆盘选择弯曲。

能量并不平分

容易误解平衡时拉伸和弯曲大致平分能量。

正确理解在能量图中，厚度从 0.037 减到 0.01 时，单位厚度的拉伸能下降 27 倍，从 $4.1\times10^{-5}$ 降到 $1.5\times10^{-6}$ ；弯曲部分只下降 12 倍，从 $5.5\times10^{-4}$ 降到 $4.7\times10^{-5}$ 。Efrati 等认为这种不平分出人意料；正因如此，很薄的薄片表现得好像不能拉伸。

模型省略了什么

容易误解真实的叶片和凝胶圆盘就长这样。

正确理解圆盘是 Efrati、Sharon 与 Kupferman（J. Mech. Phys. Solids，2009）的轴对称非欧薄板：测地半径从 0.1 到 1.1 的圆环，目标度规具有恒定高斯曲率 K，采用他们式 4.2 的约化能量，并像他们的图 4.1 一样把前置系数取为 1。径向分成 80 个单元，生长的每一帧都是用牛顿法求出的平衡态。只允许旋转对称的形状，因此褶皱、凹坑和负曲率的马鞍形都不在模型范围内；生长是给定的，而不是算出来的。探针是作用在内缘高度上的硬弹簧；松开后，圆盘按过阻尼梯度流弛豫，这是一条通往最近平衡态的路径，而不是带惯性的真实薄片的运动。

动手实验

01
让薄圆盘生长
从“薄圆盘长成球冠”开始，看一次生长过程。
你会看到： 圆盘几乎立刻离开平面，目标曲率只有 0.0186 时就开始，最后的穹顶高度是无拉伸球冠的 0.999。能量中只有 7.2% 是拉伸，颜色一直很暗。
02
把它变厚
选“厚圆盘保持平坦”，厚度为 0.45。
你会看到： 圆盘一直长到 K = 1 都没有离开平面，拉伸让它通体发亮。它的厚度高于 0.299 的屈曲厚度，平坦状态始终稳定。
03
停在刚越过临界的位置
选“刚刚越过屈曲临界”，厚度 0.26。
你会看到： 圆盘直到 K = 0.760 才离开平面，然后缓缓升起。到 K = 1 时，穹顶只有球冠的 0.577，拉伸仍占能量的 59.7%。
04
把它翻过来
选“球冠一下翻面”。
你会看到： 探针把内缘往下推。能量升到球冠的约六倍，在穹顶高度约 −0.21 处外圈一翻、能量减半。松开后，圆盘弛豫成能量相同的镜像球冠， $8.48\times10^{-5}$ 。
05
亲手推内缘
把内缘的黄色圆环向下拖过圆盘，然后松手。
你会看到： 松手太早，穹顶会弹回去；推过突跳点，它就翻面停住。你的路径以蓝色画在能量图上。

非欧薄板

生长怎样让平的圆盘弯起来？

核心数学模型

目标几何

拉伸加弯曲

无拉伸球冠

两个厚度

常见难点

这里薄并不意味着软

能量并不平分

模型省略了什么

动手实验

让薄圆盘生长

把它变厚

停在刚越过临界的位置

把它翻过来

亲手推内缘

相关实验与文章

几何弹性：薄结构怎样弯曲、屈曲和突跳