约束与权衡13 分钟

每一小块都能放松,为什么拼成整体后仍然有应力?

每根弹簧都有自己的天然长度,连成整体后却未必能同时满足。本文用三根弹簧的精确解区分边界加载、天然错配与允许构形:一个共线有应力状态可以转成无应力三角形,违反三角不等式的目标则不能靠转向消除。再核原始不相容图研究与2026生长薄片实验、方法和作者明确保留的证明缺口,说明局部相容、整体闭合、空间嵌入与平衡稳定需要分别判断。

小块各自舒服,不保证拼起来也舒服

一根天然长度为1的弹簧,可以在长度1时没有拉力。另一根也可以如此。但如果它们连接的位置、角度和闭合方式要求几项距离同时满足,局部的“最舒服”未必能拼成一个允许的整体构形。

这就是几何不相容要问的问题:是否存在一组实际位置,使所有给定的自然几何同时实现? 它依赖允许在哪个空间变形、哪些角度能转、边界能否移动,以及是否允许切开、重叠、缺陷或塑性变化。不能只看到应力就反推必有内部错配,也不能只数到环路就断言不能放松。

2025年的不相容图研究把这种冲突写成可分解的闭环量;2026年的生长薄片研究则显示,即使局部几何条件满足,向更大范围延续某种空间构形仍会遇到障碍。两类对象需要不同的相容条件,本篇先走完一个可完整求解的弹簧例子,再接回原研究。13

三根弹簧:力平衡之后,为什么还有储能

设三个节点按顺序位于一条导轨上,位置为$0,x,L$。节点之间各有一根理想轴向弹簧:两段短连接的天然长度均为1,跨两端的连接天然长度为1.8;刚度分别为$k_1=1,k_2=4,k_3=2$。导轨只允许共线运动,端点在轴向没有外加负荷。固定第一个节点只是去掉整体平移。

弹性能为

$$ U(x,L)=\frac12k_1(x-1)^2+\frac12k_2(L-x-1)^2+\frac12k_3(L-1.8)^2. $$

若三条边都处于天然长度,前两条要求$L=2$,第三条却要求$L=1.8$。这组共线约束不能同时满足。

本库实际求解$\partial_xU=\partial_LU=0$,得到

$$ x=0.885714,\qquad L=1.857143. $$

三条边的有符号拉力依次为$-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$,可以独立推出

$$ U_{\min}=\frac{\Delta^2}{2(1/k_1+1/k_2+1/k_3)}. $$

不相容给出不能消失的闭合差;刚度则决定怎样分摊这份差,以及储能多大。若把其中一条边切开,剩余两条边可各自达到天然长度。但切边已经改变了材料和连接,释放的能量不能假装没有破坏或操作代价。

同一组天然长度,离开导轨后可以没有应力

上面的应力结论有一个承重条件:节点必须共线。

允许中间节点横向移动后,天然长度$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$,则

$$ U=\frac12(Cu-w)^TK(Cu-w),\qquad C^TK(Cu-w)=0. $$

后一式就是节点力平衡。如果$w$能写成$Cu$,所有偏好差便可以实现;等价地,沿每一条闭环的有符号和都为零。若存在非零闭合差,调节点位置只能去掉其中可由位置差产生的部分。

令$S$的列张成$\ker C^T$,也就是无需外部节点力就能平衡的边力方向。消去节点位置后,有

$$ U_{\min}=\frac12(S^Tw)^T(S^TK^{-1}S)^{-1}(S^Tw). $$

这是原作者公开旧稿附录中的加权结构,本库按自己的拉力正号约定重新推导并执行;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]$,自然度量为

$$ d\bar s^2=\Phi(v)^2du^2+dv^2,\qquad \Phi(v)=\frac{A}{\sqrt{K_0}}\cos(\sqrt{K_0}v). $$

这里自然高斯曲率恒为$K_0>0$,$A$表示赤道周长相对半径$1/\sqrt{K_0}$球面的倍数。若寻找通常的旋转曲面

$$ X(u,v)=(\Phi(v)\cos u,\Phi(v)\sin u,z(v)), $$

无拉伸条件要求$z'^2=1-\Phi'^2$。当$A>1$时,这个分支到

$$ v_h=\frac{\arcsin(1/A)}{\sqrt{K_0}} $$

便达到$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,局部距离规则并没有在这里退化。趋近视界时,该分支的一个主曲率趋零,另一个发散;独立积分得到总自然曲率

$$ \bar K_T=4\pi A\sin(\sqrt{K_0}v_0), $$

在两边对称达到视界时恰为$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$的相同厚度条件。

沿子午线再次切开,作者观察到表面可以重叠、恢复接近无拉伸的球面片。原补充材料给出参数曲面

$$ X_p(u,v)=\frac1{\sqrt{K_0}}\big(\cos(\sqrt{K_0}v)\cos(Au),\cos(\sqrt{K_0}v)\sin(Au),\sin(\sqrt{K_0}v)\big). $$

本库实际微分核对,它产生相同自然度量。但$A=1.4$时,$u=0$和$2\pi$两端并不闭合,且展开超过一整圈会重叠。它因此是切开后允许重叠的无拉伸参数曲面,不是原闭合、无自穿透材料的同条件解。接触厚度和摩擦还会影响真实放松程度。

这个干预支持整体闭合方式参与了应力形成。它没有把“允许切口和重叠”偷换成“不改变任何条件也能放松”。

一条分支失败,还不能排除全部构形

这里需要保留原文自己写明的限制。作者v2正文指出,关于旋转曲面视界的论证针对该分支的扰动,尚未排除全部非对称等距构形,不是完整证明;End Matter也注明高阶非线性等距变形的一步缺少严格证明。

因此,本篇能独立确认的是:显式旋转曲面分支的视界、曲率行为、参考度量在该处仍规则,以及切开后的重叠参数曲面。作者的实验和有限元支持指定条件下的有应力形态;它们不独自证明所有允许构形均无解。若讨论有限厚度材料,还需比较弯曲与拉伸能、初态、局部稳定、接触和其他分支。

$4\pi$也需要条件。对一般有边界表面,Gauss–Bonnet关系包含边界测地曲率项,不能把它直接当所有开放薄片的普适总曲率上限。原补充材料已列出单侧生长、在总曲率$2\pi$附近先到达一个边界视界的例子。两个边界是否同时增长,与形态选择有实际关系。

怎样判断一个材料为什么没有完全放松

最有用的起点是实际列出它偏好的几何、允许的构形和当前约束,再看应力如何响应干预。

解除外载、释放夹持,可以检验边界贡献;在保持天然长度的条件下允许转向,可以检验此前是否把几何自由度锁住;切开或允许缺陷会检验闭合约束,但同时改变材料资源;改变等待时间、温度或加载历史,则有助于识别未平衡与不可逆过程。不同机制可以共存,观察到一种释放现象并不自动排除其余来源。

本文三根弹簧的有限构形给出精确见证或精确能量下限;图与线性框架公式给出指定近似内的可实现条件;薄片研究则把局部相容、整体延续和实际形态之间的缺口暴露出来。它们共同说明,材料的内应力需要在明确的自然状态与允许条件下解释,而非由“有环路”“有缺陷”或“局部都相容”中的任何一句独自决定。

来源与依据

  1. S1 原始理论/模拟研究,发表接受稿

    Statistical Mechanics of Frustrated Assemblies and Incompatible Graphs

    José M. Ortiz-Tavárez; Zhen Yang; Nicholas Kotov; Xiaoming Mao · 2025

    2025-04-07 PRL134,147401接受稿;SHA256 0373e31f7083b61660388857a729f2a85029f7d22db570833b2dcb92d120e788 · 相关全文已核

    7页含CHORUS封页;主文Frustration as incompatibility、Incompatible flow graphs、Relationship to continuous elasticity,Eq1–5

    实际核原主机制、标量关联/正交环路、等电阻例与矢量弹性分界。主文大规模MC结果未独立重做;电路耗散功率与静态弹性储能不同资源账。

  2. S2 原始作者旧稿与理论附录

    Statistical mechanics of frustrated assemblies and incompatible graphs — author v1 with appendix

    José M. Ortiz-Tavárez; Zhen Yang; Nicholas Kotov; Xiaoming Mao · 2024

    固定v1 2024-07-25,11页;SHA256 3b7d8bd6b1484dc7a42a24ff78c1984c6234949f18712efda42ebf03fa84513a · 相关全文已核

    Appendix A1–A18,PDF第5–7页:正交基、标量电池图、几何相容矩阵及加权K公式

    实际读相关附录并移植加权结构独立执行;二维以上弹簧只在固定方向一阶近似中对应线性流;未把旧稿附录当2025发表附件逐字同版,MC附录仅定位,不称复现。

  3. S3 PRL发表论文的作者修订版与完整补充材料

    Isometric Incompatibility in Growing Elastic Sheets

    Yafei Zhang; Michael Moshe; Eran Sharon · 2026

    v2 2026-07-28,31页含SM;发表PRL137,058201 DOI10.1103/bczk-bfy5;SHA256 842dbd3acaf7e4cdfc35cd76cb72d2ea4c556ec86160dac7f778dc8acf8a800b · 相关全文已核

    主文显式旋转曲面/非对称论证限制/实验/切口;End Matter;SI I–IV、V.1、VI、VIII、IX,图3名义尺寸,Eq19/20/21/30–32/41–50/59

    实际核相关主文、方法和SM,PDF/HTML交叉。主文承认不排除非对称等距构形,End Matter承认高阶步骤缺证明;4π按对称分支及边界条件限域,SI单侧2π例保留。制作为插球面楔片再粘合,切开后允许重叠;模拟自由边界、2e4–1.3e5网格收敛范围和局部相容b对照已核,未取原FE场或重做实验。图3R25/r10/t.22与generic制作R20分开,理论临界角为本库代入而非实验重新拟合。

  4. S4 本库计算证据

    本库T125加权弹簧、有限构形与薄片显式方程复算

    Codex · 2026

    2026-10-11;原显式模型评价与自建力学对照分开 · 相关全文已核

    research/calculations/geometric_incompatibility.py、_audit.py、_figures.py;assets/results.json及audit-inputs.npz

    三边k1/4/2:相容/夹持/天然错配/切边,共线到平面精确见证,违反三角不等式的精确全局下限;24非线性起点与完整力平衡/预应力Hessian,200随机标量/几何框架加权对偶。原曲面分支视界/曲率/积分和切开映射实际执行,独立差分/求根/积分核。自建模型单位,切边/导轨/接触约束明确,不冒充原MC、FE或真实应力测量。

最后更新 2026-10-11

来源