一张距离蓝图，还不是唯一的三维成品
设想把一张平片加热后变成指定球冠，或者一张人脸。可以先计算每一小块应当伸多少、缩多少，让相邻材料点之间的距离接近目标表面。但相同距离可能对应不同空间嵌入；材料还可能弯向另一面，或停在局部稳定形态。制造出来的局部伸缩也不一定等于设计值。
因此，本篇把“可靠”限定为：在预先指定的目标、材料驱动、边界朝向、初态及误差范围内，激活后的实际形态持续接近目标。形态误差须说明是否允许刚体配准、翻面和材料点重排；不能只看一张最好照片，也不能只报度量误差。
两个原始研究提供不同入口。Yang等用二元交联区域的半色调编码调节热响应水凝胶，展示球冠、鞍面及局部纹理。1 Aharoni等把液晶弹性体的定向伸缩用于复杂目标，结合近似逆设计与实验；作者也明确说明一般精确解的存在性未证。3 这些是有条件的成功实例，不能升级为任意材料、任意目标和误差下的可靠保证。
先问材料允许编程什么
把平片的自然度量记为$\bar a$，它规定材料点在放松时偏好的距离；自然曲率$\bar b$还规定偏好的弯曲。实际形态须同时满足空间几何约束，并支付偏离两者的拉伸与弯曲能量。
理想局部各向同性伸缩只有一个标量$\lambda(u,v)$：
面积比为$\lambda^2$，不是$\lambda$。水凝胶分区收缩可以据此作有效连续近似，但真实二元域也改变模量、方向响应和溶剂输运；原论文还报告过固化、扩散和图案方向的影响。因此，本库下面的标量模型是明确的均匀化假设，不是原全部制造过程。
液晶弹性体不同。原研究在目标温度使用固定的平行、垂直伸长比$\lambda,\lambda^{-\nu}$，通过面内取向$\theta$及平片参数化设计
在固定材料坐标中，旋转取向只能旋转主轴，不能任意改两条主伸长；原算法同时改变二维顶点坐标，寻找目标面各三角形映射的奇异值接近$\lambda,\lambda^{-\nu}$。不能拿固定坐标下一点的特征值不合，直接宣布所有重新参数化都无解；也不能把它当任意两个可独立指定的伸缩场。3
原local/global方法交替做局部SVD匹配与全局坐标求解，并通过网格细化、重启寻找近似解。作者对所试曲面报告度量RMS误差低于1%，同时明确保留一般精确存在性问题。本文核发表主文及固定作者v2附录的方法，没有执行作者完整Matlab算法，也没有取得有效发表版SI，故不把人脸图样当本库独立重建。3
从原始表开始，区分设计函数和实测校准
本库取得水凝胶发表主文、SI和源工作簿。Fig. 6两块各九行表给出归一化半径、面收缩比及G1–G9灰度级。令表中横坐标为无量纲$x$，两块分别符合
其中$c_+=0.493961173$、$c_-=0.315198651$。实际逐行代入，最大绝对差分别$4.12\times10^{-10}$和$4.60\times10^{-10}$。24 这核的是表值与设计函数一致，不是九次独立形态测量，也不证明空间形态精确相同。
SI第14页公式印成$r^2/R$，而正文和表变量使用$r/R$。本篇采用无量纲$x$并由表核式，不把这个有量纲记号的打印歧义直接当精确物理输入。负曲率设计的鞍形也不因径向编码就自动具有旋转对称性。
源工作簿Fig. S12还有另一组FM收缩校准：G1为0.295802、G9为0.500996；Fig. 6对应值为0.315199、0.493961。G1–G9同级表值的最大相对差约6.15%。两表作用和采样口径不能仅凭同级标签视为相同重复样本，也没有足够信息给这差异作唯一实验归因。下面限定使用Fig. 6九级表的范围；换成G0、其他校准、工艺或温度，资源边界就会改变。
走完一个球冠的逆解和正向几何核验
取初始单位圆片、$0\le r\le1$，设计
自然线元为$ds^2=\lambda^2(dr^2+r^2d\theta^2)$。它的高斯曲率
为正常数。一个显式球冠嵌入是
其球半径$R_s=\sqrt c/(2\sqrt\beta)$，直接计算可得$\rho'^2+z'^2=\lambda^2$、$\rho^2=\lambda^2r^2$。这是实际走完的等距几何见证；尚不是有限厚度的能量极小证明。
限定Fig. 6的$\eta_{\min}=0.315198651$到$\eta_{\max}=c$，这族设计要求
相应球冠高/底半径$z(1)/\rho(1)=\sqrt\beta$最多约0.50185。这是所选圆片、径向各向同性设计族及九级资源表的边界，不是所有曲面的不可能定理。
本库取$\beta=0.2$为可达连续目标，另取$\beta=0.8$为超范围对照。前者理想高度约0.261927；后者需要边缘面收缩比约0.15246，低于所限0.31520，不能靠在同一区间里优化灰度补救。
对可达目标，再把面收缩量选成最近的九级表值，并以径向宽度0.02平滑；对超范围目标，先截到区间再同样平滑。以正向几何关系
重建。所有负根号值先记录失败，绝不截成零来假造见证。平滑采用保单调插值，避免普通三次样条的平台过冲被误认为物理障碍。4
定义度量误差为原圆片面积权重$2r\,dr$下的相对RMS，形态误差为同一材料点的三维位置RMS除以目标高度；固定圆心、边界法向及正面，不允许翻面或重标材料点。九级平滑编码的度量误差约2.00%，形态误差约11.94%，边缘高度约0.18142；超范围截幅对照分别约48.76%、35.39%。这些是本库明确量化/平滑规则的结果，不是原打印工艺的误差预测。
初次计算把S12标题行误纳为数值，且普通样条在平台产生微小正导数，均在正式评分前留快照纠正，重算并独立核验。未将这两种数值问题解释为材料失败。
正向几何通过，还要过有限厚度与分支选择
为独立检验下一环节，保持上述理想$\rho(r)$，只允许高度变为$qz(r)$。$q=1$是目标球冠，$q=-1$是反向球冠，$q=0$是平的。这是一个受限Rayleigh–Ritz形态族；固定径向投影，未开放所有面内位移、非轴对称模式和接触。以下的极小与稳定性都只在这个族内成立。
在$\bar b=0$、常模量$Y=1$、泊松比0.3的薄片能量中，只有径向膜应变非零：
用自然面积积分，膜能为$hA(q^2-1)^2$，弯曲能为$h^3B(q)$；$B$由实际第二基本形式在自然正交基下的分量积分，程序保留非线性曲率，没有把它随意写成常数乘$q^2$。因此几何目标$q=\pm1$可恰好无拉伸，却仍付弯曲代价，有限厚度的极小会略变平。
实际求所有驻点、Hessian和两侧标量极小：$h=0.02$时$q\approx\pm0.999129$，$h=0.1$时$q\approx\pm0.977491$。两分支度量及能量相同，固定正面和边界法向时却是两个不同设计输出。数据另存$h=0.2,0.4$的形式延拓，但这些点不当薄片理论对真实厚体的相变预测。4
再加明确的外部选择能$-\eta hAq$，取$\eta=0.02$、$h=0.02$，正分支变为更低能，但负分支仍是局部极小。这是外加操作资源的示例，不冒充原水凝胶已具备这种压力/偏置，也不因最低能变了就说任意初态都能自行到达。
LCE原实验对复杂人脸也遇到薄极限不足及不同区域弯向不同方向的问题。作者引入上下表面的轻微取向差来诱导自然曲率，但该资源受限：严格面内锚定下，所诱导的曲率张量有一正一负的主值，方向由已定取向限制，不能任意指定完整$\bar b$。固定作者v2附录还说明，额外放大取向差的试样常更不规则，展示的是较好结果；不应把一张人脸及重复热循环当所有制品成功率。3
本库另核一个吸收系数后的曲率投影：允许张量为$kR_\theta\begin{pmatrix}0&1\\1&0\end{pmatrix}R_\theta^T$，只优化标量$k$。单位球面曲率张量$I$的最佳投影为零、残差$\sqrt2$；鞍形张量$\mathrm{diag}(1,-1)$在合适取向可以匹配，在另一取向却不能。它说明可控曲率子空间的限制，不复算原具体扭转角系数或证明完整形态稳定。
为什么小度量误差可以有大曲率误差
度量决定距离，曲率还依赖空间导数。给$\beta=0.2$的伸缩场加
用均匀乘性误差$0.003/\sqrt2$作伸缩RMS近似匹配控制。各项度量误差都约0.424%，但导数随$m$放大。直接计算内禀曲率、并用独立的圆周长度/弧长导数式核合，得到：
| 扰动 | 度量RMS误差 | 高斯曲率相对RMS误差 | 轴对称等距见证 |
|---|---|---|---|
| 均匀控制 | 0.4247% | 0.423% | 通过 |
| $m=1$ | 0.4243% | 13.1% | 通过 |
| $m=4$ | 0.4243% | 205% | 失败 |
| $m=20$ | 0.4243% | 5110% | 失败 |
失败条件是某处$|\lambda+r\lambda'|>\lambda$，令上述$z'$不能取实值。这不是所有三维构形都不存在的证明；真实薄片可以拉伸、起皱或采取非轴对称形态。高频时厚度与特征尺度的比也会变大，理想局部驱动及薄极限更需重新检验。本篇没有在这些失败点伪报实际机械平衡形态或实验失效率。
这也解释为何只约束逐点百分误差不够。相同RMS误差可以是平缓标定偏差，也可以是细尺度图案误差；二者对曲率和弯曲代价不同。真实稳健设计至少需要误差幅度、相关长度、厚度和局部刚度共同进入正向计算。
怎样把设计推进到可靠成形
| 要回答的问题 | 应检查的量 | 本篇支持到哪一步 |
|---|---|---|
| 目标在资源范围内吗 | 允许伸长/自然曲率、温度、厚度、边界和材料域 | 九级径向球冠族给出明确边界；不作任意目标存在性证明 |
| 逆解是否匹配目标几何 | 度量误差、奇异值、网格与坐标变换 | 原近似算法条件核；本库解析球冠和表函数核合 |
| 正向会形成什么 | 所有允许变形的拉伸/弯曲与接触 | 几何见证与受限Ritz族实算，未做完整壳FEM |
| 目标是否稳定且可达 | 全部允许模式、分支、初态、驱动路径 | 本库一参数族稳定性及弱偏置保留亚稳态；完整稳定性未证 |
| 制造误差下是否仍成功 | 同口径形态误差、幅度/相关长度、重复样本 | 明确离散及频谱控制，未报告真实生产成功率 |
实际推进时，先固定目标的尺寸、朝向、容许误差和驱动条件，再求材料允许的逆解；随后必须回到正向力学和加载路径，最后在有依据的误差分布与重复试样上验收。若失败，需要知道是资源越界、几何近似不足、有限厚度、分支选择还是制造偏差，而不是只把优化迭代次数加大。
这与前七篇形成直接联系：T125/T129要求先核整体相容和错配组织；T126要求开放合适形态并计厚度；T127/T128说明形状与缺陷自由度会改能量竞争；T130要求路径与加载条件；T131要求自然状态和历史载体匹配。这里的综合答案有边界：可靠逆设计是一条需要逐环核验的链，度量匹配只是其中一环。
本篇不重做T117的训练任务容量、T121的主动方向反馈或T124的后续适应。父库晶体结构预测及有效理论篇作方法边界入口；材料离线形态逆设计有自己的驱动、本构和承重证据。未决包括任意LCE目标的精确存在性、作者最终SI与完整代码执行、全部壳模式稳定性，以及真实水凝胶制造的误差统计、输运与模量耦合。
来源与复算
<a id="source-s1"></a>[1] 原始发表研究。 Yang, H.; Li, H.; Zhang, J.; Liu, T.; Qi, H. J.; Sun, H. Halftone-encoded 4D printing of stimulus-reconfigurable binary domains for cephalopod-inspired synthetic smart skins. Nature Communications (2025)。DOI/全文。发表PDF/HTML取得，核Fig6增长、二元域/模量/制造与温度方法。
<a id="source-s2"></a>[2] 原始SI和源数据。 同论文SI、出版方Source Data。SI约21.5 MB、源工作簿约23.8 MB取得；核S8/S12–S15相关段及SI第14页原图，实际读取Fig6、S12相关表。其他表和视频未作本篇全面再分析。
<a id="source-s3"></a>[3] 原始发表研究与作者旧版补充。 Aharoni, H.; Xia, Y.; Zhang, X.; Kamien, R. D.; Yang, S. Universal inverse design of surfaces with thin nematic elastomer sheets. PNAS 115, 7206–7211 (2018)。DOI，发表PMC，固定1710.08485v2。发表主文HTML和含SM的作者v2 PDF取得，核local/global、近似存在性、曲率投影、制作与试样限制。发表SI入口仅得1816字节验证页，不冒充已取有效最终附件。
<a id="source-s4"></a>[4] 本库执行。 输入表、曲线、结果。research/calculations/reliable_surface_inverse_design.py及_audit.py、_figures.py。原表18点、球冠/离散/截幅正向见证、六误差频率、六厚度/弱偏置驻点及曲率投影执行；88项独立积分、导数、标量最小和边界核通过，三图目读。原件、快照和日志在本篇evidence；不称原LCE全算法、原壳FEM或新实验复现。