小块各自舒服,不保证拼起来也舒服
一根天然长度为1的弹簧,可以在长度1时没有拉力。另一根也可以如此。但如果它们连接的位置、角度和闭合方式要求几项距离同时满足,局部的“最舒服”未必能拼成一个允许的整体构形。
这就是几何不相容要问的问题:是否存在一组实际位置,使所有给定的自然几何同时实现? 它依赖允许在哪个空间变形、哪些角度能转、边界能否移动,以及是否允许切开、重叠、缺陷或塑性变化。不能只看到应力就反推必有内部错配,也不能只数到环路就断言不能放松。
2025年的不相容图研究把这种冲突写成可分解的闭环量;2026年的生长薄片研究则显示,即使局部几何条件满足,向更大范围延续某种空间构形仍会遇到障碍。两类对象需要不同的相容条件,本篇先走完一个可完整求解的弹簧例子,再接回原研究。13
三根弹簧:力平衡之后,为什么还有储能
设三个节点按顺序位于一条导轨上,位置为$0,x,L$。节点之间各有一根理想轴向弹簧:两段短连接的天然长度均为1,跨两端的连接天然长度为1.8;刚度分别为$k_1=1,k_2=4,k_3=2$。导轨只允许共线运动,端点在轴向没有外加负荷。固定第一个节点只是去掉整体平移。
弹性能为
若三条边都处于天然长度,前两条要求$L=2$,第三条却要求$L=1.8$。这组共线约束不能同时满足。
本库实际求解$\partial_xU=\partial_LU=0$,得到
三条边的有符号拉力依次为$-0.114286,-0.114286,+0.114286$,负号表示受压。每个节点的合力为零,但最低弹性能仍为0.0114286。它已经找到了导轨允许的最低能位置,并非“还没有给它足够时间放松”。这里长度、刚度分别以$\ell_0,k_0$为单位,能量单位为$k_0\ell_0^2$,不是某种材料的实测量。4
这里保存的是静态弹性能,理想保持不需要持续供能。原不相容图研究用电池与电流说明数学分解,但电路的耗散功率由电源维持,不能直接当作静态材料的持续能耗。
对这个单环结构,令天然闭合差$\Delta=1+1-1.8=0.2$,可以独立推出
不相容给出不能消失的闭合差;刚度则决定怎样分摊这份差,以及储能多大。若把其中一条边切开,剩余两条边可各自达到天然长度。但切边已经改变了材料和连接,释放的能量不能假装没有破坏或操作代价。
同一组天然长度,离开导轨后可以没有应力
上面的应力结论有一个承重条件:节点必须共线。
允许中间节点横向移动后,天然长度$1,1,1.8$满足三角不等式。两端放在$(0,0)$、$(1.8,0)$,中间放在$(0.9,\sqrt{0.19})$,三条边便都精确达到天然长度,弹性能为零。刚度没有改变,改变的是允许的构形集合。

原共线状态在平面中虽然仍然合力为零,却不稳定:本库计算去掉刚体运动后的Hessian,最低特征值约为−0.24668;独立横向二阶差分得到相同负值。稍微向侧面移动中间节点,就能降低能量。这说明平衡、有储能、稳定和全局最低能是不同判断。
本库从正负横向方向各取初态,共12次非线性求解,都找到零能三角形。零能的依据首先是上面的精确位置见证;多初态求解是数值检查,不是用12次成功代替存在性证明。
如果将第三条边的天然长度改为2.2,结论又不同。任何欧氏空间中的三个点都满足$d_{02}\le d_{01}+d_{12}$,所以$1,1,2.2$不可能同时实现。对允许的三边距离最小化同一二次能量,最优值位于$d_{02}=d_{01}+d_{12}$,为0.0114286;它由$d_{01}=1.114286,d_{12}=1.028571,d_{02}=2.142857$达到。12次平面优化也达到该值。
这里存在的是明确的全局长度障碍,增加一个横向方向也消不掉。它与前一个可以通过转向放松的例子使用相同连接、刚度和节点数量。因此,环路数或自由度数独自不足以判断冲突。
刚度不均、外部夹持和天然错配怎样分开
保持三条刚度为$1,4,2$,把天然长度换成相容的$1,1,2$。导轨内自由平衡时,位置恰为$0,1,2$,各弹簧均无应力。刚度不均本身没有制造应变。
若再把总跨度夹到2.1,内部节点停在1.08,两段短边拉力均为0.08,跨端边拉力为0.20;边界需要提供大小0.28的反力,储能为0.014。撤掉轴向夹持后,这份应力消失。相比之下,天然长度$1,1,1.8$的闭合差在导轨内解除轴向夹持后仍在。
| 相同刚度下的条件 | 导轨内平衡储能 | 有区分力的改变 |
|---|---|---|
| 天然长度相容,端点可移动 | 0 | 刚度不同仍可无应力 |
| 天然长度相容,总跨度夹到2.1 | 0.014 | 解除夹持可恢复零能 |
| 天然长度错配,端点可移动 | 0.0114286 | 解除轴向夹持不足以消除闭合差 |
| 上一组切掉跨端连接 | 剩余弹簧为0 | 连接已改变;不再是同一材料条件 |

实际材料还可能有未平衡、塑性、热膨胀差或对天然状态的错误估计。撤去载荷后仍有应力,排除了当前载荷的一部分解释,却没有单独证明来源。切开后的翘曲也要与边界、厚度、材料损伤和测量参考状态一起解释。
图上的闭合差怎样成为一项不能放松的能量
对标量位置差模型,定义关联矩阵$C$:每条边一行,乘节点位置向量$u$便给出该边的有符号实际差。天然差记作$w$,正刚度对角阵为$K$,则
后一式就是节点力平衡。如果$w$能写成$Cu$,所有偏好差便可以实现;等价地,沿每一条闭环的有符号和都为零。若存在非零闭合差,调节点位置只能去掉其中可由位置差产生的部分。
令$S$的列张成$\ker C^T$,也就是无需外部节点力就能平衡的边力方向。消去节点位置后,有
这是原作者公开旧稿附录中的加权结构,本库按自己的拉力正号约定重新推导并执行;2025接受稿正文以等电阻的标量图说明正交分解。12
在三边例子里,$S$只有一个独立方向,便退化成前面的闭合差/总柔度公式。一般图中不能把任意选出的环路基直接当成正交投影:本库一个双环例中,未经归一化的基给出平方范数0.5275,正确投影为0.0559375,独立力平衡能量为其一半0.02796875。基之间的重叠需要计入。
本库还核了100个随机标量图和100个带真实边方向的线性框架。加权节点最小化与自应力空间计算的拉力最大差约$5.9\times10^{-15}$。这一检查支持公式实现,并没有把标量闭环条件升级成任意弹簧的有限构形定理。
在二维、三维弹簧框架中,长度变化还取决于边的方向,应使用几何相关的相容矩阵$R$。$Ru=w$描述的是给定构形附近的一阶变化;它的自应力空间未必等于普通图的环路空间。原附录明确保留了这项线性近似。前面的三角形例子正说明,转向等有限变化会改变这个判断。
薄片更难:局部度量规则,旋转曲面仍会走到尽头
薄片的自然几何,不再只有若干边长。自然度量$\bar a$规定小距离,天然曲率$\bar b$规定偏好的弯曲;实际构形有自己的$a,b$。局部Gauss与Mainardi–Codazzi–Peterson条件检查两者能否来自局部空间曲面,整体还要满足周期闭合、边界、光滑性与嵌入条件。
Zhang、Moshe与Sharon的2026研究,考虑周向周期的环带:$u\in[0,2\pi)$,$v\in[-v_0,v_0]$,自然度量为
这里自然高斯曲率恒为$K_0>0$,$A$表示赤道周长相对半径$1/\sqrt{K_0}$球面的倍数。若寻找通常的旋转曲面
无拉伸条件要求$z'^2=1-\Phi'^2$。当$A>1$时,这个分支到
便达到$z'=0$;再往外,右侧变成负数,不能继续得到同一类实旋转曲面。这是作者称为geometric horizon的边界。3
本库复算$K_0=1,A=1.4$时,$v_h=0.795603$。到$v_0=0.9$,$z'^2=-0.202658$,但自然度量行列式仍为0.757342,局部距离规则并没有在这里退化。趋近视界时,该分支的一个主曲率趋零,另一个发散;独立积分得到总自然曲率
在两边对称达到视界时恰为$4\pi$。

实验与切口提供了什么证据
作者用硅橡胶球壳切出环带,再沿子午线切开、插入球面楔片并重新粘合,增加周长因子$A$。原实验观察到越过相应范围后出现周期凹坑;模拟使用自由边界,并比较了两种自然曲率设定,包括局部Gauss/MCP相容的设定。本文已核作者修订版的制作、有限元和几何方法,未取得原有限元场或独立重做实验。
作者补充图3给出的名义输入为球半径25 mm、极区切口半径10 mm、厚度0.22 mm。代入对称旋转曲面条件,得到名义临界$A\approx1.09109$,对应插入楔角约32.79°。这是本库根据图注输入计算的理论阈值,未从实验照片重拟合转变角。实验$t/R=0.0088$,也不是模拟所述$t\sqrt{K_0}\le0.002$的相同厚度条件。
沿子午线再次切开,作者观察到表面可以重叠、恢复接近无拉伸的球面片。原补充材料给出参数曲面
本库实际微分核对,它产生相同自然度量。但$A=1.4$时,$u=0$和$2\pi$两端并不闭合,且展开超过一整圈会重叠。它因此是切开后允许重叠的无拉伸参数曲面,不是原闭合、无自穿透材料的同条件解。接触厚度和摩擦还会影响真实放松程度。
这个干预支持整体闭合方式参与了应力形成。它没有把“允许切口和重叠”偷换成“不改变任何条件也能放松”。
一条分支失败,还不能排除全部构形
这里需要保留原文自己写明的限制。作者v2正文指出,关于旋转曲面视界的论证针对该分支的扰动,尚未排除全部非对称等距构形,不是完整证明;End Matter也注明高阶非线性等距变形的一步缺少严格证明。
因此,本篇能独立确认的是:显式旋转曲面分支的视界、曲率行为、参考度量在该处仍规则,以及切开后的重叠参数曲面。作者的实验和有限元支持指定条件下的有应力形态;它们不独自证明所有允许构形均无解。若讨论有限厚度材料,还需比较弯曲与拉伸能、初态、局部稳定、接触和其他分支。
$4\pi$也需要条件。对一般有边界表面,Gauss–Bonnet关系包含边界测地曲率项,不能把它直接当所有开放薄片的普适总曲率上限。原补充材料已列出单侧生长、在总曲率$2\pi$附近先到达一个边界视界的例子。两个边界是否同时增长,与形态选择有实际关系。
怎样判断一个材料为什么没有完全放松
最有用的起点是实际列出它偏好的几何、允许的构形和当前约束,再看应力如何响应干预。
解除外载、释放夹持,可以检验边界贡献;在保持天然长度的条件下允许转向,可以检验此前是否把几何自由度锁住;切开或允许缺陷会检验闭合约束,但同时改变材料资源;改变等待时间、温度或加载历史,则有助于识别未平衡与不可逆过程。不同机制可以共存,观察到一种释放现象并不自动排除其余来源。
本文三根弹簧的有限构形给出精确见证或精确能量下限;图与线性框架公式给出指定近似内的可实现条件;薄片研究则把局部相容、整体延续和实际形态之间的缺口暴露出来。它们共同说明,材料的内应力需要在明确的自然状态与允许条件下解释,而非由“有环路”“有缺陷”或“局部都相容”中的任何一句独自决定。