======================================================================
1) 表 3 误差复算 & 冲击速度
E= 8 J  v0=3.47 m/s  位移误差=+1.8%  吸能误差=-23.6%  吸能比(exp)=0.424
E=10 J  v0=3.89 m/s  位移误差=-3.9%  吸能误差=-13.1%  吸能比(exp)=0.421
E=12 J  v0=4.26 m/s  位移误差=+1.1%  吸能误差=-5.7%  吸能比(exp)=0.406
正文写 10 J 实验吸能 4.12 J: 3.66 J 对应误差 -11.2%, 3.27 J 对应 -20.6%  (与 -13.1%/-22.3% 不符 → 4.12 为笔误, 应为 4.21)
======================================================================
2) Puck IFF 在固定断裂角 θ=53° 下的单轴强度（表 1 公式原样）
S23A = 75.3 MPa, mu_nt = 0.287
(a) 纯横向拉伸: 预测失效 σ22 = 111.7 MPa  (输入 S22t = 57.7)
    → 若 θ 恒取 53°，拉伸 IFF 起始被推迟约 1.94 倍；Puck 原理论中纯拉伸断裂面为 θ=0°
(b) 纯横向压缩 [表 1 原式  S23A + mu*σnn]: 预测失效 |σ22| = 128.8 MPa (输入 S22c = 199.8)
(b) 纯横向压缩 [摩擦增强式 S23A - mu*σnn]: 预测失效 |σ22| = 199.8 MPa (输入 S22c = 199.8)
======================================================================
3) 裂纹带模型单元尺寸上限 l_max = 2 G E / S²  (N-mm-MPa)
纤维拉伸: l_max = 6.44 mm
纤维压缩: l_max = 6.91 mm
基体拉伸: l_max = 3.57 mm
基体压缩: l_max = 1.04 mm
1 mm×1 mm×0.2 mm 单元: 体积立方根 = 0.585 mm；面内边长 1 mm
======================================================================
4) 内聚区长度与界面罚刚度
IMS/977-2 (表2): l_cz,I = 4.80 mm, l_cz,II = 6.27 mm, 1 mm 网格下内聚区内单元数 ≈ 4.8 / 6.3
T300/YH69 (表4): l_cz,I = 0.60 mm, l_cz,II = 1.41 mm, 1 mm 网格下内聚区内单元数 ≈ 0.6 / 1.4
Turon 建议 K ≥ 50·E3/t  (t=0.2 mm): 2.48e+06 N/mm³  vs 论文 5e+03 N/mm³
Turon 建议 K ≥ 50·E3/t  (t=2.4 mm): 2.06e+05 N/mm³  vs 论文 5e+03 N/mm³
厚向柔度: 铺层 4.85e-04 mm/MPa, 23 个界面 4.60e-03 mm/MPa → 等效 E3 = 944 MPa (铺层 E3=9900)
横向剪切柔度: 铺层 9.70e-04, 界面 4.60e-03 mm/MPa → 等效 G13 = 862 MPa (铺层 4950)
======================================================================
5) 显式积分稳定时间步估计 Δt ≈ L_e / c_d
面内 E1, L=1 mm: c = 9860 m/s, Δt ≈ 1.01e-07 s
厚向 E3, L=0.2 mm: c = 2472 m/s, Δt ≈ 8.09e-08 s
2.5 ms 加载步 ≈ 30901 增量步 (未计质量缩放/内聚单元)
======================================================================
6) 网格研究（Fig.11 读图峰值，kN）与细化比
细化比 r21 = 1.333, r32 = 1.500 (非恒定)
ε21 = +0.10, ε32 = -0.50, ε21/ε32 = -0.20 → 振荡收敛, 标准 GCI 不适用
粗网格相对细网格峰值偏差 ≈ +5.2%  (读图)
======================================================================
7) 式(1) 剪切非线性系数核对：对 Fig.2(c) 读图点做三次拟合
γ 以 % 计: τ ≈ 0.801γ³ -10.34γ² +51.39γ  (MPa)  → 初始 G12 ≈ 5.14 GPa
γ 为无量纲应变: c1 ≈ 8.010e+05, c2 ≈ -1.034e+05, c3 ≈ 5.139e+03 MPa
论文表 2: c1 = 10.2, c2 = -0.5, c3 = 0.1 MPa —— 在上述两种单位下均无法复现曲线
  用表 2 系数(γ以%计) 算 τ(1.0%) = 9.8 MPa
  用表 2 系数(γ以%计) 算 τ(5.0%) = 1263.0 MPa
