相同的膨胀，首先要说明相同的是哪一层
把一张受刺激会改变尺寸的薄片放在桌上，它可能整体弯曲，也可能出现许多细皱纹，或者仍然平坦却带着应力。仅说“材料膨胀了多少”还没有固定问题：面内各方向偏好的距离、沿厚度两侧的差异、边界能否跟着伸长，都可能不同。
2021年Siéfert等用预拉伸乳胶双层和热激活打印PLA制作自然度量平坦、自然曲率方向却随位置变化的薄带，观察到光滑环形、应力局部集中和分面等形态。2026年Maor与Mora从三维弹性出发证明了一组薄带极限结果，解释哪些几何类有形态转变，哪些没有。两者提供实验与严格理论的不同证据，不能把2026年论文当作一组新实验。12
本篇先固定自然几何研究厚度和宽度，再单独改变边界与支撑研究起皱。两组最小模型不拼成一张通用于所有薄片的相图。
自然度量和自然曲率，分别在要求什么
自然度量$\bar a$规定材料面内偏好的距离；自然曲率张量$\bar b$规定它偏好怎样弯。中面实际形状$f$产生实际度量$a_f$和曲率$b_f$。各向同性薄片的能量可写成
这里$Y$是三维杨氏模量，$t$是厚度；式子省略泊松比与张量收缩的系数，只说明尺度。弯曲刚度为$B=Yt^3/[12(1-\nu^2)]$，面内模量量级为$Yt$，所以在相同特征长度下，越薄越倾向用弯曲减少面内应变。但曲率差与面内应变有不同单位，不能只比较$t$与$t^3$的数值。
实际度量与实际曲率必须来自同一张空间曲面，不能独立任选。对平坦自然度量$\bar a=I$，无应变构形须满足Gauss条件
以及Codazzi条件
自然目标如果违反前者，称Gauss不相容；若满足前者却违反后两式，称Codazzi不相容。比如$\bar b=\kappa I$要求两个方向都弯，乘积$\kappa^2$非零，与平坦度量冲突。另一个目标$\bar b_{11}=k_0+\beta y$、其余分量为零，行列式处处为零，但$\partial_y\bar b_{11}=\beta\ne0$，仍不能在整片上同时实现。1
均匀、沿厚度也均匀的各向同性膨胀是一个必要反例：若边界自由、自然曲率为零，整片可以平坦地增大尺寸，既不必弯曲也不必起皱。若夹住它不许增大，留下的压缩则来自边界限制。
固定$\bar a,\bar b$也不等于固定原始双层膨胀百分比。双层的自然曲率通常随厚度与层间应变差共同变化；扫描$t$时要维持同一$\bar b$，就须相应调整制造或激活条件。本文的固定自然几何比较明确采用前一种控制。
一个能走完的例子：弯成双曲率，还是牺牲一个方向
取长$L$、宽$w$的浅薄带，$\bar a=I$、$\bar b=\kappa I$固定。用两个近似常量曲率$a,b$描述高度
这是本库自建的小斜率试探族，尚不包含任意褶皱、扭转、接触或完整有限转角。先放松平面位移、忽略长带端部层，得到单位自然体积、再除以$Y$的能量
第一项是为实现非零曲率乘积而付出的拉伸，第二项是偏离自然曲率的弯曲代价。这里每一项都无量纲，乘回$YtLw$才是物理总能量。
系数$1/1440$可以独立核。长带内部的纵向膜力为
横向膜力和剪切为零；纵向合力也为零。将$N_{xx}^2/(2Yt)$沿宽度平均，再除以$Yt$，即得到$w^4a^2b^2/1440$。$N_{xx}$的单位是力/长度，不是三维应力。本库还从实际多项式平面位移求应变、用高斯积分核同一结果。4
令
约化能量为
它可以完整分类驻点。对称解$x=y=s$满足
非对称解则满足
在$\nu=0.3$时，本模型的分支临界值为$\chi_c=10.4518$。小$\chi$时，两曲率接近自然值；大$\chi$时，一项接近$(1+\nu)\kappa$，另一项接近零，趋向单向弯曲。这不是完全无能量的放松：自然曲率的另一方向仍未满足。
本文对全部解析驻点比较能量，并以和$s=x+y$、积$p=xy$将能量配方检验全局选择；这里的“全局”限于这两个曲率变量的试探族。真实薄片还有族外变形模式，曲率轴的取向也可能由端部与边界选定，不能据此宣称完整弹性问题的唯一最低能形状。
一个有数值的控制是$\kappa=0.2$、$w=0.08$、$L=1$，长度以同一参考单位计。$t=0.001$时，最低曲率约为$(0.19638,0.19638)$；$t=0.00005$时，约为$(0.24867,0.01133)$。自然度量、自然曲率、宽度与杨氏模量均不变。或者固定$t$增宽，也会增大$\chi$。这解释的是变形代价如何竞争，不是所有膨胀薄片都会出现相同临界数。
端部与初态，为什么仍要另核
上式的膜力在纵向端部并不逐点为零，因而它是长带内部近似，不能当有限自由端薄片的完整解。
本库保持二次高度场不变，让两项平面位移在最高3、6、10、14、20、26次的多项式空间内自由放松，边界平面位移不固定。空间扩大时拉伸能单调下降。最高26次结果中，长宽比20的系数是长带值的0.96315倍，长宽比2.5时为0.70513倍。端部松弛在较宽短片上明显改变系数。4
这是一组有限试探空间的上界；即使20次与26次结果接近，也不证明已经求到完整连续体最低值。它足以说明，忽略端部得到的10.4518不应直接贴到有限样品上。
初态检查还发现更直接的陷阱。在$\chi=15$时，严格对称初态能让BFGS返回“成功”，停在$x=y=0.553903$；梯度接近零，但Hessian最低特征值为−0.693966，是鞍点。微小破坏对称后，求解得到$(1.116231,0.183769)$或互换解，能量由0.838630降到0.766410。
48次多初态运行中，39次通过驻点与非负Hessian检查并达到约化最低能，7次为驻点鞍点，2次未达到驻点容差；全部保留。求解器“失败”也可能只是精度判据，而“成功”不足以证明稳定。本文没有把这些数值优化轨迹冒充真实材料的时间动力学或能垒测量。
2026严格结果，支持到哪一步
Maor与Mora使用三维预应变度量、满足转动不变性和正定弹性控制的能量密度，从三维变形推到薄带极限；能量按体积平均。长度已无量纲化，基本条件是$t\ll w\ll1$，不是把有单位的毫米厚度直接与毫米宽度的平方比较。1
对自然度量平坦、沿中线有Gauss不相容的薄带，最低体积平均能量的尺度为
窄带极限要求$w^2\ll t\ll w$；宽带结论使用先$t\to0$、再$w\to0$的迭代极限。两种极限中的中线曲率目标不同，给出形态转变。上面的简化模型具有相同的尺度竞争，但其系数、分支形状和阈值是本库模型结果，不能说由严格定理逐数值推出。
Codazzi不相容的情况不同。若沿中线的纵向自然曲率$\bar b_{11}$处处非零，原定理给出两种极限均为$t^2w^2$尺度，没有这一领先阶的转变。对于$\bar b_{11}=k_0+\beta y$、$k_0\ne0$的例子，柱面可以精确满足平坦度量和中线曲率，仅在离开中线后偏离目标。这个柱面试探解的弯曲代价为
这是明确可积分的上界，与原尺度一致；本文不称它是该问题的精确最小值。不能从“仍然不相容”推出“一定有和Gauss情形同样的转变”。
若$\bar b_{11}$在中线上消失，转变才有可能，但这只是必要条件，不是充分条件。原文§5.4对不同Codazzi几何保留不同程度的未决：一例已证明能量尺度发生变化，但最优幂次仍未确定；另例提出试探形状和预期尺度，未证明该尺度最优。文中的$w^{8/3}$、$w^{7/2}$等讨论含条件推断，不能一律宣传为已经完成的严格临界律。
这些定理控制近似最低能序列的渐近行为，不直接给出有限厚度误差棒、局部能垒、摩擦、接触、动态到达时间或全局无自穿透约束下的最终外形。数学极限和实际可达需要在相同材料、边界及观测协议下接起来。
同样的多余长度，为什么有时是一道拱，有时是细皱
起皱另有必要条件：某个方向出现不能自由释放的压缩，且存在允许的面外模式。自然几何不相容可以提供压缩，边界也可以提供；单纯“更薄”不足以说明是哪一种。
Cerda与Mahadevan的原始研究展示了纵向拉伸的聚乙烯薄片也会出现横向皱纹：夹持抑制端部横向收缩，造成非均匀面内应力。其理论用弯曲和有效支撑的竞争选波长，支撑既可是真基底，也可来自另一方向的张力。所谓“张力下起皱”，不等于局部所有方向都在拉伸。3
本库用一个更简单、可完整核的对照：长度$L$的薄条，限制横向曲率为零，端点高度为零、转角自由；固定投影跨度，使自然长度比允许跨度多出比例$\delta$。其局部线性支撑满足压力$p=-Kh$。令$Y_2=Yt$、$B=Yt^3/[12(1-\nu^2)]$，取$h=A\sin(qx)$，$q=m\pi/L$，每单位投影面积能量为
其中$A^2q^2/4$是面外起伏收纳的平均多余长度。最优振幅满足
在此模型中，失稳门槛是
没有支撑时$K=0$，最低门槛对应$m=1$，即整体拱起。$K>0$时，长波需要较大高度，会增加支撑代价；短波曲率大，增加弯曲代价。连续波数近似给出
有限端点再从整数$m$中选择。在这套小斜率模型内，加入多个正弦模式也不会降低最低门槛：由正交性，多余长度的收纳量相加，各模式的弯曲与支撑代价也是逐项相加，所以最省的做法是将收纳量分配给$\delta_c$最小的模式，简并时可有混合。
若$\delta$低于所有门槛，平坦态仍是本模型最低能，但带着压缩应变$-\delta$；这里“仍然被拉伸”更准确地说是仍保留面内弹性应变，其符号须实际记录。
取$Y=1$MPa、$\nu=0.3$、$L=20$mm、$\delta=0.002$，支撑$K=10^5$N/m³。本库实际检查160个整数模式：
| 厚度 | 最低门槛$\delta_c$ | 实际半波数 | 振幅 |
|---|---|---|---|
| 20μm | 0.0008561 | 22 | 19.57μm |
| 50μm | 0.0013534 | 11 | 29.43μm |
| 100μm | 0.0019343 | 7 | 14.75μm |
| 200μm | 0.0027126 | 0，平坦 | 0 |
若保持20μm厚度却撤掉支撑，最低模式改为1，拱高约569.28μm。若进一步允许投影跨度随自然长度增长，则$\delta=0$，平坦无应变状态的能量恰为零。
数值模量、支撑和尺寸是本库示例输入，没有对特定实验拟合。定形的局部常量支撑、铰支端、零横向曲率与小斜率都是条件；真实半空间基底的支撑会依赖波长，不能直接套用这里的四分之一次方及其系数。该模型也没有给自由二维薄片的完整褶皱分类。
实验观察怎样接回最低能和可达性
2021年原薄带实验的自然曲率来自沿方向收缩的双层，且带本来具有面内测地曲率。它是Codazzi几何的实验，不是本文$\bar b=\kappa I$的字面复现。主文方法记录乳胶贴合、PLA打印层厚和加热收缩；乳胶径向曲率的变化只在$R\gg w$时近似忽略。2
作者的有限元采用超过3500个三角形、梯度下降和多个初态；对于分面数还改变了初始周期扰动。原文明确承认实验与数值只代表局部能量极小，数值未惩罚表面相交。本文没有取得其逐点数据或重跑原有限元，因此只用它核实际制造、形态观察与求解限制，不将图形相似当作全局最低能证明。
这给出有区分力的推进方式：固定自然几何和边界，分别测面内应变、曲率与分项能量，再改变厚度或宽度；把解除夹持、改变支撑与改变自然几何分别操作。若只换初态就出现不同持久形态，应核局部稳定和能垒；若多初态求解一致，也只能说明已检查的初态，没有自动排除未探索分支。
较稳的机制是厚度改变弯曲与面内代价比例，宽度改变相容条件在整带上施加的限制，边界与支撑改变可以怎样释放压缩。本文得到的是这些机制在明确模型和原始研究中的条件分界。一般自然度量、有限厚度误差、接触、不可逆松弛和动态路径仍须针对具体材料处理。
来源与复核
<a id="source-s1"></a> [1] Maor & Mora (2026). Rigorous Analysis of Shape Transitions in Frustrated Elastic Ribbons, Archive for Rational Mechanics and Analysis 250, 96，2026-10-03。已取得60页发表PDF与固定作者2503.11779v3，核定义、材料假设、主定理与§5.3—5.4相关证明/未决；不称逐行重证全部60页。
<a id="source-s2"></a> [2] Siéfert, Levin & Sharon (2021). Euclidean Frustrated Ribbons, Physical Review X 11, 011062。原始实验/理论/模拟，发表全文已取。核自然几何、形态和§IV制造/有限元方法及局部极小、无交叉惩罚的原文限制；原数据与代码未复现。
<a id="source-s3"></a> [3] Cerda & Mahadevan (2003). Geometry and Physics of Wrinkling, Physical Review Letters 90, 074302。作者站PDF已取；前四页是原论文，末附一页解读不作为一手证据。核Eq1—6、图1与边界/近似注释；未复算其特定二维实验。
<a id="source-s4"></a> [4] 本库实际计算 (2026)。 全部输入与结果。research/calculations中的frustrated_ribbon_shape_selection.py、_audit.py、_figures.py实际执行；48次求解含失败/鞍点、有限平面位移试探空间、整数模式和直接积分均留记录。自建模型、独立核算与原作者研究分开；原件和读取范围保存在本篇evidence及sources。