织构分析：性能应用推导
r 值从哪里来：从宽厚收缩推到 Taylor 变分和在线估计
用功共轭关系解释为什么对横向收缩比求最小功，再给出数值手算、误差传播及在线衍射估计的完整逻辑。
先记住这三点
- r 是塑性应变比，测量方向和应变区间必须一致。
- 最小功条件对应自由侧面的应力条件，并非纯粹的经验求极值。
- 从织构估计 r 需要本构模型，不能把预测称为直接测量。
先从长度、宽度和厚度说起
在板面沿 e₁ 拉伸，e₂ 为板面内宽向，e₃ 为厚向。塑性伸长、收缩应变为 ε₁ᵖ、ε₂ᵖ、ε₃ᵖ，定义 r=ε₂ᵖ/ε₃ᵖ。通常后二者都为负，因此 r 为正。采用对数塑性应变且塑性体积近似守恒时，ε₁ᵖ+ε₂ᵖ+ε₃ᵖ=0，所以 r=−ε₂ᵖ/(ε₁ᵖ+ε₂ᵖ)。
例如课堂算例取 ε₁ᵖ=0.12、ε₂ᵖ=−0.08，则 ε₃ᵖ=−0.04，r=2。不能把“长度增加 12%”直接当作对数应变 0.12，也不能在弹性比例明显时把总应变代替塑性应变。测量通常还应限定均匀变形阶段和所用应变区间，颈缩区的局部几何不符合简单均匀假设。
用一个未知量表示收缩分配
把轴向塑性应变率记为 ε̇>0，令宽向收缩占比为 q，写 D(q)=ε̇diag(1,−q,−(1−q))。这样 trD=0 自动满足。这里先假定加载方向具有所需对称性，宏观塑性剪切可取零，或明确采用零剪切的近似约束。宽厚比 r=q/(1−q)，反过来 q=r/(1+r)。r=1 对应宽厚各承担一半收缩，r=2 对应 q=2/3，即厚向只承担三分之一。
q 常在正常宽厚同时收缩条件下取 0 到 1，但这一范围是所讨论变形方式的条件，不是任意各向异性应变状态的数学定律。若得到 q≥1 或 q≤0，应检查模型是否预测了某方向扩张、求解约束是否合理，不能不加解释地把异常结果截断。
为什么要对 q 最小化
对每个取向 R，在给定 D(q) 下做 Taylor 最小功，得到塑性耗散 ΨR(q)。以体积权重平均得到 Ψ̄(q)=∫ΨR(q)f(R)dμ。若全部活动系共用 τc，可写 Ψ̄=τc ε̇M̄(q)。未知的不是外加轴向伸长，而是材料自由选择的横向收缩分配。
对 q 作虚变化 δq，有 δD=ε̇diag(0,−δq,δq)。功共轭应力给出 δΨ̄=Σ:δD=ε̇(Σ33−Σ22)δq。单轴拉伸的两个自由侧面要求 Σ22=Σ33=0，因此内部最优点须满足 dΨ̄/dq=0，即 dM̄/dq=0。压力不确定性由侧面牵引条件固定；只求导为零尚不足以自行证明两个横向应力都为零。
最小值不能只靠解导数等于零
率无关滑移的最小功是分段线性的凸函数，活动滑移组合变化时可能有折点。此时应要求 0 落入次梯度集合；数值上比较相邻区间与边界的耗散，不能用光滑函数的求导公式硬套。若使用平滑拟合或速率敏感模型，则可求导，但还须验证极值是最小值且位于允许区间。
教学例设 M̄(q)=2.8+1.8(q−0.65)²，则导数 3.6(q−0.65)=0，q*=0.65；二阶导数 3.6>0，得到 r=0.65/0.35≈1.857。此抛物线只用于把求解过程演示清楚，不是任何真实钢板测得或模拟得到的曲线。真实 M̄(q) 必须由取向与滑移约束计算。
加载方向变化怎样进入计算
拉伸方向与 RD 夹角为 α 时，先将加载坐标旋到该方向，再旋转滑移系。在无剪切耦合的对称方向，可继续求 M̄(q,α) 的最小值，再得 r(α)。一般离轴拉伸却可能出现 D12 等剪切分量；应把允许的剪切率也列为未知量，与 q 一起最小化平均功，使各共轭宏观剪应力为零。若仍预先把全部剪切率锁为零，得到的是附加运动学约束下的近似，未必满足自由单轴拉伸。
同样，先对各晶粒独立求最优 q 再平均，一般不等于先平均功再求统一的 q；后者才对应统一应变假设。常用 r̄=(r0+2r45+r90)/4 和 Δr=(r0−2r45+r90)/2。若三值为 1.8、2.2、1.6，便得 r̄=1.95、Δr=−0.50。平均 r 较大与平面方向差较小是不同目标，Δr=0 也不意味着所有中间角度完全相同，更不能单凭它排除六耳、八耳等高阶方向变化。
误差为什么在大 r 时特别敏感
由 r=q/(1−q) 求导，dr/dq=1/(1−q)²。当 q=0.8 时 r=4，q 的微小误差 0.01 会带来一阶近似 dr≈0.25。厚向收缩越小，分母越接近零，求比值越敏感。报告高 r 时尤其需要说明应变测量误差、塑性分量提取和模型假设。
织构误差还会经过 M̄(q) 的最小位置传到 q。若 M̄ 的谷底很平，多个 q 都几乎同功，微小 ODF 或 CRSS 变化即可显著移动最优点。因此不能只输出一个无小数不确定性的 r 值；要说明结果是否对带宽、权重和滑移族参数敏感。
在线估计不是直接看见 r
毛书印刷页 213–216 介绍固定几何下的能量色散衍射：E=hc/λ 与 λ=2d sinθ 联立为 d=hc/(2E sinθ)。固定 θ 后，不同能量峰对应不同面间距；少量方向上同时记录多族反射，可约束低阶织构系数，再经塑性模型估计 r。
这条链条依次是光子能量与强度、晶面/方向分布、受限 ODF 参数、模型性能。少量观测未必能恢复尖锐的完整 ODF，之所以仍可能支持性能估计，是目标对某些低阶信息更敏感。是否足够必须结合目标的灵敏度和验证，不能把“低阶通常有用”说成高阶一定无关。
真实工程案例应该怎样迁移
书中深冲钢案例把 {111} 面接近板面的取向分布与较高 r 联系起来，并明确提醒单晶单滑移表格不能定量代替多晶钢板。学习时可先用其单晶例理解收缩方向，再用本节多晶最小功建立计算；材料等级、间隙原子、硬化和可能的 BCC 非 Schmid 行为都要补上。
与铝板制耳研究对照时，区分单轴 r(α) 与拉深过程中变化的多轴应力路径。r 是一条有用的性能线索，但完整杯形预测还需要成形模型和工艺边界。这样读者才能把一个公式用在它真正支持的问题上。
来源与延伸阅读
以下链接用于核对概念、标准状态或原始文献。本页为中文教学整理，不复制受版权保护的标准正文。