Reynolds–Finan 异质核 FEM 分析 —— 审查与修改报告

§0先说结论

但有三个层级的问题，其中第一层是会被直接拒稿的（见页首三卡）。

TIER A数值验证问题（必须修，改完反而加强结论）

1.1　主结果选错了离散化 —— 最该改的一条

把 Table S2 的五个离散化各自算了组间变化：

两族的行为：

• P1 族还在漂：−37.9 → −40.9 → −42.3，每次加密还在往下走 3.0 pp、1.4 pp，外推极限约 −43.6 ~ −44.5%。
• P2 族已经收敛：粗网格 −44.64%、中等网格 −44.51%，只差 0.13 pp，而且是在网格量差 3.4 倍的情况下。

也就是说：coarse_h1p25_P2（14 646 dof）比 main_h0p45_P1（22 593 dof）又便宜又准。现在把最不收敛的 P1 结果当头条，是把结论往小里报了 2 个百分点，还给审稿人送了一把刀。

1.2　高渗组的 h 收敛是非单调的，"passes h-refinement" 不成立

Control：0.060352 → 0.065003 → 0.066064，单调，增量比 4.38，收敛阶 p≈2.89，GCI≈0.59%。这条没问题。

Hyperosmotic：0.037475 → 0.038412 → 0.038092。增量符号翻转（+0.000937 → −0.000320）。这是振荡收敛，Richardson/GCI 无法计算。现在报的"高渗组 h 加密差 0.84%"是振荡的巧合，不是收敛度量。

1.3　P1 与 P2 外推到了不同极限 —— 真信号，不是噪声

Control 的 P1 三点 Richardson 极限是 0.066378，而 P2 在中等网格就给出 0.067416，差 1.54%。同一个连续问题的两族离散化不可能收敛到不同极限。只有两种解释：

1. P2 还没收敛（但 P2 粗/中两网格只差 0.13 pp，不像）；
2. P1 和 P2 求解的根本不是同一个连续问题——因为染色质材料场是按积分点/单元采样 50³ 体数据赋值的，换阶次、换网格 = 换材料分布。

第 2 条几乎肯定是主因，理由见下条。

1.4　网格根本分辨不了模型的两个关键尺度

从 dof 反推（P1 向量场，3 dof/节点，椭球体积 604.6 µm³）：

• 壳层：0.12 µm 的 lamina，最细网格的单元也比它大 3.6 倍。壳不可能是被解析的三维实体层。它到底是什么单元（膜单元？面上的 surface energy？厚度被放大的实体层？）全文没有一处说明。"current-thickness-equivalent resultant proxy" 的定义也没有公式。审稿人无法评估头条数字。
• 染色质：输入是 50³ 体素（约 0.30 × 0.22 × 0.14 µm/体素），而主网格间距 0.43 µm。材料场被欠采样一个数量级，10 个灰度组在网格上严重混叠。这直接导致 §1.3 的极限不一致，也导致"h 加密收敛"在数学上不成立——加密的同时材料场也在变，你不是在收敛到一个固定问题。

1.5　残差判据的三处问题

1. 文本说 <1e-6，Table S1 写 1e-07——两个门槛不一致，选一个。
2. 用的是绝对残差范数。绝对残差随网格、载荷、单位制变化，没有参照量就没有意义。必须报相对残差（除以外载荷范数或首次 Newton 残差）。这是数值审稿人的标准问题。
3. coarse_h1p25_P1 / Hyperosmotic 的残差是 2.67e-8，比其余所有算例（~1e-10）大 100 倍，而全部 UQ 都跑在这个粗网格上。得说明这是不是普遍情况——建议把 141 对 UQ 的残差分布（直方图/最大值）放进 Figure S5。

1.6　P1 线性四面体 + 近不可压核 = 体积锁定

核的 K/G = 30 → ν = 0.4835（近不可压）。线性四面体是近不可压有限变形弹性里最差的单元，体积锁定是教科书级别的问题。这正好解释了：为什么 P1/P2 同网格差 3.58%/2.69%（p 误差 > h 误差），为什么 P1 族收敛慢而 P2 族一到粗网格就稳。

1.7　完全没有稳定性检查

Newton 收敛 + J>0 不等于平衡是稳定的。你有一个薄壳（0.12 µm），UQ 里还出现了整体平均 resultant 穿越零（→ 净受压）的算例。真实 lamina 在净压缩下会皱曲/屈曲，这个模型（光滑 Newton、对称网格、无缺陷、无接触）既不会捕捉屈曲也不会告诉你它错过了分叉点。

1.8　缺代码本身的验证

自写的 scikit-fem 工作流，带乘法生长分解 F_e = F·F_g⁻¹ 和初应力残差化——这两块最容易写错，但全文只有"projected MINRES 与 augmented direct 在粗 P2 上机器精度比对"（那只验证了线性解法器，没验证残差/切线）。

TIER B力学建模与结论逻辑问题（最致命的一层）

2.1　"141 个点全部卸载"不是稳健性证据，是构造使然

这是最担心的一条，也是审稿人一定会问的。

模型里壳的张力来自两项：初始预张力 0.05 mN/m 和 375 Pa 探针。现在施加各向同性自然收缩 F_g：

• coupling = 0 → 壳自然构型不变，被核往里拽 → 壳面积减小 → 张力降；
• coupling = 0.7 → 壳自然构型也缩 → 张力更降。

在参数盒的任何一点，几何上都不存在张力升高的路径。所以 64/64 和 77/77 不翻号是必然的，跟"稳健"没关系。把它当作 Main Figure 3 的核心内容，等于用一整张主图证明了一件由模型定义决定的事。

更要命的是 §5 里的表述读起来像发现，其实是定义："All 64 paired Latin-hypercube points and all 77 parameter-domain boundary/anchor points predict a lower Hyperosmotic global mean"。

2.2　−42% 里约 76% 是那个纯假设的预张力在松弛

换算一下（t = 0.12 µm）：

预张力占 Control 端点的 76%。而且 UQ 的 pretension 范围是 0.025–0.1，从来没有包含 0，所以这个分解从头到尾没做过。

2.3　缺"均质核对照" —— 异质染色质可能对结论完全无关

论文的输入卖点是 Reynolds 的 10 级异质染色质图。但端点是外表面的面积加权全局平均。核内异质性往外表面传递时极可能被平均掉。

如果把核换成体积平均（或 Voigt/Reuss）的均质材料，端点还是 −44.5%，那么整套数字化图像、灰度标定、10 级映射的工作对头条结论毫无贡献。审稿人一定会问。

2.4　"shear probe" 与实际施加的载荷自相矛盾

正文写 "self-equilibrated 375 Pa nominal shear probe"、Table S4 写 probe_shear_stress_pa，但 §9 第 2 节写的是 "逐点牵引为 t₀ = P̄·N"。

t₀ = P̄N 是沿法向的牵引，即压力，不是剪切。各向同性压力加在 K/G=30 的近不可压体上基本什么也不会发生（这也许正是端点被预张力主导的原因）。

量级参照：200 mOsm 的渗透压差在 310 K 下约 515 kPa，比 375 Pa 探针大约 1 400 倍。你没有把渗透压当外压加（这是对的），但正文最好主动点明这个量级对比，否则读者会困惑。

2.5　预张力的参考构型不自平衡

带预张力的闭合壳需要核内有平衡压力（类似压力容器 T = pR/2）。只往壳里塞 T₀ = 0.05 mN/m 而没有配平，参考构型不处于平衡，第一次 Newton 就会把它松掉，实际实现的预张力不等于 0.05。

正文说 "residualized reference-configuration initial stress"——如果 "residualized" 指减掉 u=0 处的内力残差（常见技巧），那这个初应力就成了只进切线不进平衡的虚拟量，而它偏偏占端点的 76%。

2.6　把 lamina 当 ν=0.35 的可压缩三维实体是错的物理

从 Table S2 反推可确认最小 J_e 出现在壳里：min J_e(shell) × J_g(shell) = 0.69278 × 0.725517^0.35 = 0.61918——与表里 minimum_total_jacobian 精确吻合，说明 coupling 按 J_g,shell = J_nat^0.35 实现（这条内部一致，没问题）。

但结论是：壳被体积压缩了 31%（J_e = 0.693）。真实 lamina 在面内近乎面积不可压、通过减薄响应，不会损失 31% 体积。用 ν = 0.35 的可压缩 Neo-Hookean 实体做 0.12 µm 的 lamina，会产生大量伪体积弹性应变，直接污染 resultant proxy。

2.7　材料参数表内部矛盾

2.8　物理上被略过、但必须至少写进 limitation 的四件事

1. 染色质压实后不变刚：F_g 把核压了 27%，但 G 场保持不变。Finan 那篇的核心正是非线性渗透行为和 321 mOsm 的"变紧阈值"。压实必然增刚。至少跑一个 G ∝ ρⁿ（或 G(J_e)）的敏感性算例。
2. 用线性 Ponder 关系外推到 580 mOsm：论文标题就叫 Nonlinear Osmotic Properties，你却用单一 φ_s 的线性关系外推。必须论证 580 mOsm 仍在适用范围内，或改用 Finan 的非线性形式，或把 φ_s 的不确定度真正传播进 UQ（而不是手挑 0.65–0.85）。
3. 无孔弹性/无时间：渗透响应是水的输运问题（秒–分钟尺度），核是孔弹性/黏弹性的。超弹性平衡态 + 强加自然体积删掉了所有速率效应。直接关系到实验：Flipper 在加压后第几分钟测？没到平衡则模型假设与实验条件不匹配。
4. 核是悬空的：无细胞质、无细胞骨架、无 LINC 耦合。真实高渗休克下整个细胞在缩，骨架在给核加载，这可能主导实际方向。

2.9　各向同性 F_g 与实测形态不符

强加的是各向同性自然收缩，但贴壁细胞核在高渗下主要是变扁（高度方向优先塌陷）。各向异性 F_g 会显著改变壳 resultant 分布。而 §11 恰好提出用"核体积/轴比"做湿实验读出——那就更应该把轴比作为模型的可检验预测做出来。

2.10　端点量的选择有缺陷

"面积加权的有符号 resultant 平均"，在有大面积受压区时会正负抵消（你们自己就撞上了一个穿零算例）。一个壳可以有巨大的等量拉压而平均为零。这也是为什么不得不引入"绝对差"作为对抗性目标。

TIER CUQ 设计与可证伪性

3.1　这个预测无法被证伪 —— 最深的科学问题

如果参数盒里不存在任何一点给出"加载"，那么 Flipper 的方向性检验就没有区分力：它无法把你的模型和一个平凡假设（"东西缩了就松了"）区分开。一个不可能错的预测不是检验。

3.2　UQ 全部跑在 710 顶点网格上，幅度带已知偏差

141 对（282 次求解）全在 coarse_h1p25，而已知该网格把名义点的卸载低估 6.7 pp（−37.9% vs 收敛的 −44.5%）。所以 "P5/中位/P95 = −65.5% / −27.4% / −9.1%" 整体带系统偏差。

好消息是偏差方向保守：粗网格系统性低估卸载幅度（名义点 −37.9 → −40.9 → −42.3；最不利顶点 corner_101011 从 −4.74% 到中等网格 −5.69%，幅度也变大）。

3.3　UQ 参数盒漏掉了影响最直接的参数

resultant N = σ·t，壳厚 t 与端点成正比，但 0.12 µm 被固定死、零不确定度，provenance 还标着 assumption。它不会在相对变化里抵消，因为两组的当前厚度不同。

同样缺席的：几何（半轴、长宽比）、核-壳界面耦合（真实是 LAD 部分耦合，比 "shell shrinkage coupling" 更重要）、K/G 比、灰度→模量映射的形状（下一条）。

3.4　灰度标定只约束了 3 个分位点

用低/中/高相比例 57%/31%/12% 做单调标定，这个想法很好，应该在正文更突出。但它只钉住 3 个分位点，箱内映射形状（S 形 vs 分段线性 vs 阶梯）完全自由，而模量跨度是 23.5 倍（170–4000 Pa）。

另外：出版图的灰度不是定量 DNA 密度（显示 gamma、亮区截断、有损压缩）；8 张 z 面板线性插值到 50 层是很强的 z 分辨率外推；遮盖 panel A 标记后那块区域填了什么值也没写。

3.5　"顶点覆盖不能证明连续内部" —— 这个自我限制其实可以解决

3.6　统计表述

• 对确定性设计报 "5/50/95 百分位" 会被读成概率分布，即使声明了也一样。改报 min/median/max 或全分布图。
• −164.35% 这种百分比没有意义（分子穿零）。全文把 "traditional percentage change" 降为次要，主指标统一改绝对差（mN/m）——你们在对抗性设计里已经用了绝对差，这个逻辑应该贯穿全文。
• 那 6 个落在 ±10% 带里的算例是信息量最大的，应列出它们的参数组合（那才是结论的边界），而不是只给一个计数。

§4表述、精度与格式

4.1　有效数字虚高（几乎每个数字都超报）

离散化误差 ≥ 3.6%（同网格 P1/P2），跨离散化的组间变化散布 6.7 pp，但正文在报：

规则：报告精度不能超过验证精度。这是审稿人最容易抓、也最伤信任的一类问题。

4.2　0.735031 与 0.725517 不是同类量

前者是外体积比（含收缩较少的壳），后者是核的自然体积比。拿它们说 "close" 是在比不同的东西。应该比同类：核的平衡体积比 vs 0.726。

4.3　语言与格式

• 全文中英混排（§1、§6、§9、§10 大段中文）。投稿版必须统一英文。
• um → µm；phi_s → φ_s；1e-6 → 10⁻⁶；V/V0=phi_s+(1-phi_s)C0/C 应排成公式。
• 参考文献 1–5 在同一段里连排，需要拆成条目。
• prospective_validation_pending、corner_101011 这类程序内部标识符不应出现在正文散文里。
• Nature 主刊正文不用 "1. What was modelled" 这种编号章节标题。
• "traditional percentage change" 应为 "relative change" / "percent change"。
• 载荷命名 shear/normal 自相矛盾见 §2.4。
• 附带说明：正则粗提文本时看到 Table S4 有 XML 片段外溢，复核后确认是解析器把 9 个自闭合 <w:p .../> 标签吃掉造成的，docx 文件本身格式完好，不用管。

4.4　可复现性

• "fixed-seed" 但没写种子数值。
• 需要完整环境锁定（不只 scikit-fem 12.0.2，还有 numpy/scipy/BLAS 版本与线程数——MINRES 结果对这些敏感）。
• DOI 占位符的处理方式（§10 结尾那句禁令）是对的，保留。

§5Flipper 验证方案的问题

5.1　观测量不匹配，且方向关系未建立

Flipper-TR 的寿命主要反映脂膜张力/脂质有序度，而你们明确声明 proxy 不是脂膜张力。那么 "lamina resultant 下降 → Flipper 寿命下降" 这个方向关系本身就是一个未经证明的假设。现在的逻辑是：用假设 X 去检验预测 Y，而 X 和 Y 的联系没人验证过。

5.2　高渗休克对 Flipper 有多重混杂

• 膜曲率变化（核变小 → 曲率升高）本身改变探针构象；
• 高渗直接改变脂质堆积/有序度，与张力无关；
• 胞内拥挤度、折射率、黏度改变影响 FLIM 拟合；
• 温度、成像时长。

一个方向性变化不足以判定力学假说。

5.3　更便宜的第一道验证被漏掉了

模型预测的三轴轴比变化是共聚焦直接可测的，成本远低于 FLIM，而且是模型的非平凡输出（体积比是输入，轴比不是）。现在只在 §11 一句带过，既不在 claims 矩阵里，也不在任何图里。

图审计第二轮：9 张图逐张核对（2026-08-11）

Main Fig 1–3 与 S1–S6 全部与正文/Table S2 数字一致（−42.3%、0.0661/0.0381、V_nat=0.726/V_eq=0.735、corner_101011 的 −0.002078/−0.002486、Spearman 各值均对上），没有图文矛盾。但图里有六个此前只能推测、现在有了直接证据的问题：

F1　压缩面积 28.6% → 62.6%（Fig S4）—— "全局卸载"的实相是大半个壳转入受压

S4 面板 d 标注：对照组压缩面积 28.6%，高渗组 62.6%；P5 最小主合力从 −0.04 恶化到 −0.075 mN/m。也就是说高渗态下壳的大部分面积处于最小主合力受压。这直接加重两条原有批评：真实 lamina 在这种状态下会皱曲，无屈曲模型报告的压缩应力是伪像（§1.7）；有符号均值的下降有相当一部分来自符号混合而非均匀松弛，端点量必须换（§2.10）。

F2　作者自己的 "locking audit" 面板实锤了锁定：同网格 P1/P2 峰值差 72%（Fig S2）

S2 左下面板（自题 "Near-incompressibility / locking audit"）：mean shell 3.5% 过 10% 门槛，但 resolved shear max 差 ~21%、principal stress max 差 ~72%。锁定不是我的猜测——正文只说"未通过门槛"，图上显示的是差了七倍于门槛。§1.6 的换单元建议从"应该做"升级为"必须做"。

F3　高渗组在载荷因子 ≈0.62 处残差瞬时飙升两个数量级（Fig S1）—— 疑似接近分叉

S1 的 "Intermediate P2 nonlinear convergence" 面板：高渗组残差在载荷因子 0.3–0.62 段持续爬升到 ~5×10⁻⁸（比基线高两个数量级），随后回落。这是求解器路过近失稳区的典型指纹，位置正好落在 §1.7 要求做特征值检查的地方。建议对载荷因子 0.5–0.7 区间输出投影切线最小特征值。

F4　变形量级是 O(1) 的，"探针"绝不是小扰动（Fig S1）

位移色标最大 ~2.1 µm（长半轴 7.5 µm 的 ~28%）；对照组（无收缩、仅探针+预张力）min J 也降到 0.78。量纲上：壳膜刚度 E·t ≈ 2000 Pa × 0.12 µm = 0.24 mN/m，而探针合力 ~P·R/2 ≈ 0.94 mN/m——载荷超过膜刚度，模型工作在大应变区，P1 锁定在此最严重。另外 Fig 2d 标注对照组 V_eq=1.002 > 1：探针在轻微吹胀核——外法向牵引，再次印证 §2.4 的 "shear" 命名与实现不符。

F5　Spearman 面板预告了均质对照的结果：chromatin scale ρ = −0.08（Main Fig 3c）

染色质模量整体 0.5–2 倍缩放与端点变化的相关只有 −0.08（probe 也只有 +0.06）——模量整体翻四倍几乎不动全局端点。这强烈暗示 §2.3 均质对照的结果将是"端点对异质性不敏感"，该对照更加必须做，而且论文叙事大概率要按"异质性只影响局部场"那条路线改写。

F6　残差门槛不一致是实质问题，不是笔误（Fig S5）

S5 面板 b：128 个 UQ 平衡解的最大残差散布到 ~1×10⁻⁷ 附近（若干点 6–9×10⁻⁸），图上门槛线画在 1e-6。按正文的 1e-6 全部通过；按 Table S1 写的 1e-7，若干 UQ 算例贴线甚至越线——两个门槛的矛盾（§1.5）直接改变部分算例是否通过验收，必须统一并重新声明。

可视化小问题

• S1/S2 的收敛图把 coarse P1 → coarse P2 → intermed P1 → intermed P2 → main P1 混排在一根轴上，掩盖了 P1/P2 两族各自的行为（P1 还在漂、P2 已平）；S2 柱状图 y 轴从 0 起，12% 的散布视觉上被压没。应按 h 分组、P1/P2 分两条线、y 轴放大。
• Fig 2c 直方图右缘的蓝色尖峰是 1–99% 截断的堆积 bin，应在图注说明。
• S3 的 max principal stress 点云呈盐椒噪声（相邻点在 −0.5 与 +2.5 kPa 间跳变）——这正是 P1 近不可压逐点应力噪声的直观形态，exploratory 定位正确，但应注明噪声来源是单元阶次而非物理。

代码审计第三轮：源码逐项核对（2026-08-11）

审计对象：submission_ready_nature_optimized_20260810(1).zip 中 05_Source_Data 的 8 个 Python 文件（共 9 650 行，重点 nonlinear_solver.py 2 298 行与主运行脚本 1 586 行）。此前存疑的四件事全部定案，其中一条是我此前批错了，撤回。

C1　撤回 §2.4 的指控：探针确实是剪切，名实相符

代码构造均匀对称名义应力张量 σ = τ(e_x⊗e_y + e_y⊗e_x)（τ=375 Pa），施加牵引 t₀ = σ·N 于参考外表面。因 σ 为常张量且对称，∫σN dA = 0、合矩为零——"self-equilibrated shear probe" 的实现是对的，文档 §9 的 "t₀ = P̄·N" 里 P̄ 是应力张量而非压力，是我误读。遗留问题缩小为表述：正文应写出该张量定义，并明确这是参考构型 dead load（代码用 reference normals，两组因此受严格相同的外载——这正是 "same probe" 的实现，值得明说）。

C2　证实 §2.5：预张力是 u=0 时严格零贡献的残差化初应力

first_piola += (F−I)·S₀、切线 += I⊗S₀：参考态残差为零，预张力只随变形激活。因此对照组实际实现的壳张力是求解输出，不等于标称 0.05 mN/m。且 resultant proxy 里把推前初应力 F·S₀·Fᵀ/J 加进了 Cauchy 应力——0.066 这个头条数字确实包含全部预张力贡献，机制分解（建议 8）的必要性升级。同时代码里还有一个正文完全没提的 initial_prestress_heterogeneity 选项（默认 0，未启用，无碍）。

C3　重大利好：均质核对照的代码路径已经存在

core_heterogeneity_model="homogeneous" 是现成分支，走 cfg.core（E=800 Pa、ν=0.45）。建议 9（第二优先级）就是改一个参数重跑两次的事。同时 Table S1 的参数矛盾（§2.7）解密：800 Pa/0.45 是均质回退值，异质运行时被十组 G∈[170,4000]、K/G=30 覆盖——表里必须标注清楚这两套参数各自何时生效。

C4　proxy 公式找到了（正文必须补写）

外表面积分点上：取切向两个主 Cauchy 应力（含推前初应力），乘以当前厚度 = 0.12 µm × 法向厚度伸缩，其中法向伸缩以 J/面积伸缩运动学估计；mean resultant = (σ₁+σ₂)/2 × t_cur，面积加权用当前面积。这印证它是代理量而非穿厚积分（代码内部注释也如此声明）。公式本身自洽，但手稿正文一处都没写——建议 5 维持，照此写入方法学即可。

C5　证实 §1.3/1.4/2.6：逐单元材料赋组 + gmsh 薄层实体壳，且网格生成脚本缺失

染色质是 CellwiseCoreHeterogeneity——每个四面体按 50³ 场赋一个组号，材料场随网格变化（P1/P2 极限不一致的根源证实）。壳是 gmsh 网格 ellipsoid_core_shell.msh 里内缩 0.12 µm 的实体四面体薄层，ν=0.35 体积压缩伪像（§2.6）成立，且薄层四面体在 h=0.45–1.25 µm 下必为极端长宽比单元，单元质量需报告。可复现性缺口：生成 .msh 的脚本不在压缩包里，与 §10 "custom code included" 的声明冲突，投稿前必须补入。

解析分解第四轮：闭式推导（2026-08-11）

这一轮不需要跑任何算例，只用代码里的公式做闭式推导，再用���发表的两个端点数字对账。结论改写了正文该怎么写，也撤回了我前面的一条批评。可运行脚本：17_analytic_membrane_reference.py（纯 numpy，无依赖，直接出下面全部数字）。

A1　平衡用 (F−I)·S₀，报告用 F·S₀·Fᵀ/J —— 两者不是同一个本构声明

各自都站得住，但合在一起意味着：参考构型之所以平衡，是因为从残差里减掉了一个参考态常数项——模型隐含假设存在一个把预应力参考态撑住的虚拟平衡载荷，并把它当 dead load 一路带下去。手稿从未提及。这必须写进方法学：一个受预张力的闭合壳不可能自己平衡，必须有核内压力配平。

A2　预张力对报告值的贡献是一个与变形无关的常数 —— 我的原批评撤回

proxy 取切向主 Cauchy 应力 × 当前厚度，而当前厚度按 Nanson 取 λ_n = J/J_s。对预应力项这两个因子恰好抵消：

N_pre = ½·tr_tang(F S₀ Fᵀ/J)·t₀·(J/J_s) = S₀·t₀·(λ₁/λ₂ + λ₂/λ₁)/2

等双轴时恒等于 S₀·t₀ = 0.05 mN/m，且由均值不等式永远不小于它（多出来的部分纯粹是面内各向异性，与体积变化无关）。脚本用 200 个随机变形数值验证：最小值恰好落在 0.050000。

A3　但扣掉这个常数后，弹性合力变号了

三个后果，都比我原来那条批评更重要：

1. 头条百分比被一个假设常数稀释了。−42.34% 的分母有 76% 是 provenance 为 assumption 的加性常数。把预张力在它自己的 UQ 区间（0.025–0.1）里挪一挪，报告的百分比会大幅改变而力学几乎不变——这正好解释了 UQ 里 pretension 的 Spearman ρ=+0.29 从何而来。主指标必须改用绝对差（与 §3.6 建议一致）。
2. 高渗态下壳处于净弹性受压。这与 Fig S4d 的压缩面积 62.6% 独立吻合。真实 lamina 受压会皱曲，而模型不含屈曲——稳定性检查（脚本 05）从"应该做"变成"不做就不能报"。
3. 那个穿零的 UQ 算例不是异常点，是同一现象在更低预张力下的表现。

A4　闭式膜合力：方向性是本构的推论，不是发现

对实现的应变能施加平面应力 σ₃₃=0，解得：

N_elastic = (μ·t₀/g²)·(1 − (λ_t/λ)²)

脚本数值验证：N 在自然态 λ = g 处严格为零（求得零点 0.963631 vs g = 0.963257），且在 λ∈[0.80, 1.10] 上严格单调递增。各向同性自然收缩必然降低 λ，因此必然降低 N。

脚本补算脚本下载

全部脚本只驱动现有代码，不修改任何上游源文件——两个没有命令行开关的选项（均质核、稳定性模态）用运行时 patch _options 实现，上游代码树保持干净可 diff。每个脚本都支持 --dry-run，先看计划再跑。

⬇ 下载完整报告（Word，五轮整合版） ⬇ 下载全部脚本（zip）

06_collect_results.py 已用已发表的 Table S2 数据做回归验证，输出与手工复算逐位一致：Control/P1 观测阶 2.89、GCI 0.59%；Hyperosmotic/P1 因增量变号被拒绝给 GCI；跨离散化散布 −37.91% ~ −44.64%。收敛判定是内置的、不迁就——非单调就不出数并说明理由。

这批脚本改不了的四件事

以下只能改代码或改写正文：（1）材料场随网格变（逐四面体赋组），需网格无关的 L² 投影；（2）壳是 ν=0.35 的可压缩实体薄层，最小 J_e≈0.69 即 31% 体积压缩，应改膜/近不可压公式；（3）proxy 公式与探针应力张量定义正文缺失（代码里都有，照抄即可）；（4）A1 那个虚拟平衡载荷假设必须写明。

第二梯队脚本（已完成）

脚本 14 已在已发表的 141 个设计点上实跑，产出了 §E3/E4 的全部结论。脚本 15 顺带给出一个可证伪的一阶预测：壳厚在 0.06–0.25 µm 内变动会把头条百分比从约 −21% 拉到 −67%（对比 shell_young_pa 的 −17.0% ~ −74.6%），即一个从未被扫过的输入具有顶级杠杆——脚本把这个估计与求解值并排打印，便于证伪。

第五轮脚本落地后的新证据（2026-08-11）

第二梯队 7 个脚本写完并在已发表数据上实跑，又查���五件事。第一件是目前为止最严重的数值问题。

E1　粗网格的染色质相分数超出项目自己的验收门槛 3 倍，而全部 UQ 都跑在这张网格上

投稿包 heterogeneous_two_group_summary.csv 里就带着这个数据，我逐行核验：

目标是 Reynolds 报告的 57%/31%/12%，项目自定门槛 ≤0.02。粗网格偏差 0.0629，超标 3.1 倍。而 numerical_acceptance.json 里的门槛键名全是 primary_phase_fraction_*——只对 primary mesh（intermediate）检查，粗网格的失败因此从未浮出水面。后果有两层：

1. 全部 141 对 / 282 个平衡解（UQ + 边界应力设计）都跑在这张网格上（reynolds_finan_epistemic_uq.json 记录 mesh: coarse_h1p25, element_degree: 1）。所以 UQ 不只是幅度有偏，它求解的是一个数字化未通过验收的染色质场。
2. 粗网格还是 P1 三点 Richardson 外推的基点。外推建立在一个材料场超差 3 倍的点上。

E2　精化 §1.3：混叠与锁定是两个可分离的问题，我之前混为一谈

关键观察：同一网格上 P1 与 P2 的相分数逐位相同（元素集合不随阶次改变）。所以：

• 同网格 P1/P2 的 3.58%/2.69% 差异 → 材料场完全一致 → 这是纯锁定，不是混叠。§1.6 的诊断成立且更确定。
• 跨网格才会重采样材料场 → h 加密被污染，尤以粗网格为甚（E1）。

我原文说"P1/P2 外推到不同极限，几乎肯定是材料重采样为主因"——这句要修正为：P1 的 h 序列因基点材料场超差而不可靠，P2 与 P1 的同网格差则是锁定。两个病因分别对应两条不同的修法（换单元 vs 网格无关投影），不能只做一个。

E3　Sobol 分析独立证实了预张力的"假灵敏度"

用 141 个已解算点拟合代理模型（留一交叉验证 R²=0.9585），对绝对差做方差分解：

预张力按已发表 ρ 排第 3，按总效应 Sobol 排第 6，指数为零。这与第四轮的解析结论完全吻合、且方法完全独立：预张力是加性常数，只抬高百分比的分母，对绝对差几乎无作用。发表的 ρ=+0.29 是归一化的假象，不是力学灵敏度。

另外两点：一阶指数之和 0.842，交互作用占方差 15.8%——秩相关完全看不见这部分；shell_coupling 按 Sobol 排第 2，比按 ρ 排的第 4 高一位。

E4　顶点覆盖确实不界定连续内部——而且代理找到了比所有已解顶点更接近翻号的点

单调性检验：shell_young_pa 上有 0.4% 的随机直线出现反转，最大反转 1.07× 交叉验证 RMSE。虽弱，但足以说明极值可以落在盒内部而非顶点——手稿那条"顶点覆盖不证明连续内部"的自我限制是必要的，不是过度谨慎。

连续域搜索给出的最不利点预测值 −0.001884 mN/m，优于最佳已解顶点 corner_101011 的 −0.002078，且加 1σ 后上界触及零。也就是说"全域卸载"这个claim，按代理自身的不确定度衡量并不安全。脚本已生成对应的真实求解命令——必须解出来才能下结论，代理预测不是证据。

E5　两个更小但需要交代的问题

• 残差化预应力在壳自然态留下伪应力。在 F = g·I（壳自身生长态，本应无应力）处，(F−I)·S₀ 不为零，每个切向分量残留 −15.31 Pa = 预张力的 3.67%，折合 0.001837 mN/m ≈ 报告差值 0.027972 的 6.6%。这是残差化公式的固有代价，必须在方法学里声明。
• 过渡带只有 1/6 被复核过，且已发表的 4 个复核全是 intermediate P1——即只加密了 h、没有加密阶次。而 E2 已证同网格 P1/P2 差 3.58% 才是主误差源。所以"已在细网格复核"的说服力比看上去弱。
• 界面耦合确实没有旋钮（核与壳是同一 basis 的单元子集，界面节点就是同一批未知量，无界面积分项）。§3.3 建议里的"界面耦合加入参数盒"需要先改代码改网格，代价远高于壳厚——壳厚是纯 config 字段，先做壳厚。

§6优先级与工作量

第一梯队 —— 不做就会被拒（约 2–3 天，多数是算力）

第二梯队 —— 决定能投什么档次的刊物（1–2 周）

第三梯队 —— 提升科学分量