织构分析:性能应用推导

从 ODF 到弹性与磁性:张量旋转、均匀化和低阶系数

推导为什么晶体性质需要先旋转再平均,解释 Voigt、Reuss 与球谐系数加权的意义,并用立方晶体展示有织构也未必产生某种各向异性。

阅读设置

先记住这三点

  1. 标量、二阶张量、四阶张量的旋转与平均规则不同。
  2. 均匀应变平均刚度,均匀应力平均柔度,两者不能互换。
  3. 性能只读取某些取向统计,并不意味着 ODF 的其余细节不存在。
01

从一支导热方向开始

设晶体坐标中热流 qᶜ=−Kᶜ∇Tᶜ,R 将晶体系转到试样系。因为 qˢ=Rqᶜ、∇Tᶜ=Rᵀ∇Tˢ,代入即得 qˢ=−RKᶜRᵀ∇Tˢ,所以 Kˢ=RKᶜRᵀ。两个指标各带一个 R,这是二阶张量旋转的来源,不是记忆规则。

若单晶为立方对称,普通线性导热张量必须是 Kᶜ=kI,因此 RKᶜRᵀ=kI。无论织构多强,单靠这一二阶单晶性质都不会产生方向差。这不排除孔隙、层状组织、第二相或晶界造成宏观导热各向异性;它提醒我们不能笼统说“有织构就使所有性能各向异性”。

02

把同样过程推到弹性

晶体系线弹性写 σᶜᵢⱼ=Cᶜᵢⱼₖₗεᶜₖₗ。应力与应变都按二阶张量旋转,代入并逐指标整理得到 Cˢᵢⱼₖₗ=RᵢₐRⱼᵦRₖ𝚌Rₗ𝚍Cᶜₐᵦ𝚌𝚍。四个指标各乘一个 R。若采用 6×6 Voigt 记号,还须区分工程剪切和张量剪切;直接把一个 3×3 R 左右乘 6×6 刚度矩阵在维数上就不成立。

立方单晶弹性一般有 C11、C12、C44 三个独立常数,而各向同性仅需要两个;只有满足 C44=(C11−C12)/2 时才达到弹性各向同性。由此可理解立方晶体的导热可能各向同性,弹性却仍有明显方向性,两者并不矛盾。

03

Voigt 与 Reuss 各自假设什么

若假设所有晶粒承受同一应变 ε,晶粒应力为 Cˢ(R):ε,平均后 Σ=〈Cˢ〉:ε,所以 CV=∫Cˢ(R)f(R)dμ。若假设各晶粒同一应力 Σ,则 ε(R)=Sˢ(R):Σ,平均柔度 SR=∫Sˢ(R)f(R)dμ,得到 CR=SR⁻¹。这里 S=C⁻¹ 指四阶张量在对称应变空间中的逆。

一般有 〈C〉⁻¹≠〈C⁻¹〉,所以先平均再取逆与先取逆再平均不是同一操作。Voigt 和 Reuss 在相应线弹性均匀化条件下给出能量界限;Hill 平均常取 CH=(CV+CR)/2 作为估计,不是由此证明真实组织必然在算术中点。真实空间组织和相互作用可能需要自洽或全场方法。

04

用两个晶粒的标量类比算清平均顺序

设两个等体积分量的一维模量为 100、200 GPa。统一应变下,EV=(100+200)/2=150 GPa;统一应力下,SR=(1/100+1/200)/2=0.0075 GPa⁻¹,ER=133.33 GPa。差异来自不同的场分配,不是某个平均算错。

这一例子只是说明代数顺序,不能把各向异性晶粒沿某轴的 E 直接平均,冒充完整三维弹性均匀化。完整计算先平均四阶刚度或柔度,再在沿 a 的单轴应力、横向自由条件下得到 1/E(a)=aᵢaⱼSijkl aₖaₗ。若横向应变受约束,测得的是另一种约束模量,不能继续套用这一自由拉伸公式。

05

为何低阶 ODF 系数能控制一些性能

在正交归一基 {Tₙ} 下,写 f=ΣcₙTₙ、M=ΣmₙTₙ,则〈M〉=∫Mf dμ=Σmₙcₙ,实基时不需复共轭;复基应按所用内积处理共轭。非单位归一的书中广义球函数会留下相应范数因子,例如 1/(2l+1)。所以不能把不同教材的系数直接混用。

如果 M 对某些高阶基函数的投影 mₙ 为零,那些 cₙ 对该线性平均没有贡献。四阶弹性张量旋转和截断的磁晶各向异性能是典型的有限阶信息入口。非线性塑性、局部失效和晶粒连通则不自动满足相同结论。低阶参数可用于某种性能,不等于足以重建任意细节的织构。

06

磁晶能的三个方向手算

立方晶体截断磁晶能写 E=K1(α1²α2²+α2²α3²+α3²α1²)+K2α1²α2²α3²。方向余弦满足 Σαi²=1。沿 [100] 得 E100=0;沿 [110],α1²=α2²=1/2,得 E110=K1/4;沿 [111] 得 E111=K1/3+K2/27。给定常数后就能比较这些方向,而不必仅背诵易轴名称。

对相应温度成分下 K1 为正、K2 不改变排序的铁基材料,[100] 是有利方向。Goss {110}〈001〉使一个这种方向平行 RD。若晶体在 (001) 面内偏离易轴小角 θ,则 E≈K1cos²θsin²θ≈K1θ²,说明小角散布通过均方偏离影响这一能量近似,角度要用弧度。

07

织构到磁性能之间还缺什么

上述能量描述单晶磁化方向代价,不能直接等同于铁损。交流磁场中的总损耗还包含磁滞、涡流及其他动态贡献,依赖厚度、电阻率、磁畴、应力和频率等。EBSD 可测取向分布,为磁晶能或更完整磁模型提供输入,却不直接测得铁损。

毛书印刷页 197–200 的硅钢与 Fe–Ni 实例适合用来联系易轴和工艺;铝箔实例则把取向联系到选择性腐蚀。两者提醒读者:先找单晶层面的方向性机制,再写出如何平均和有哪些其他变量,最后才讨论材料性能,不能从一张 ODF 跳过这些环节。

08

附表该如何使用,而不是逐个数字背下来

毛书附录将反演和广义球函数中固定的 α′、K、Q、B 系数列到 lmax=22,用途是减少重复计算。使用前先锁定晶体和试样对称性、基函数相位、归一化与索引约定。附表中 l=12 等阶出现多列独立对称函数,不能把同阶误当成只有一个系数。

可靠的学习方式是先用低阶多项式检查,例如 P2(x)=(3x²−1)/2,令 x=cosα 后得到 P2(cosα)=1/4+3cos2α/4;这展示角度函数如何变成余弦系数。然后再学习本书的对称归一组合。教材附表不是实验数据,也不是无需约定就能直接塞进任意软件的通用常数表。

来源与延伸阅读

以下链接用于核对概念、标准状态或原始文献。本页为中文教学整理,不复制受版权保护的标准正文。