双轴拉伸下近α钛合金中高能{13‾41}孪晶界及层级孪晶行为研究
双轴拉伸下近α钛合金中高能{13‾41}孪晶界及层级孪晶行为研究
1. 研究背景
近α钛合金（如TA31）因其高比强度、优良断裂韧性、耐腐蚀性与可焊性，广泛应用于航空航天结构件及深海装备。在实际服役过程中，这些部件往往承受多轴复杂应力状态，例如深海潜航器耐压壳便处于典型的双轴加载状态。然而，长期以来材料的结构设计与性能评估仍主要依赖单轴拉伸试验数据，对双轴乃至多轴应力状态下的变形行为与微观机制缺乏系统认识。
六方结构（HCP）金属由于对称性低、独立滑移系有限，孪生在c轴应变协调中扮演关键角色。已有研究表明，双轴加载可能改变滑移与孪生之间的竞争关系，甚至激活非传统孪晶模式。然而，关于双轴加载下钛合金中孪晶行为的报道存在明显分歧，且对高指数、高能量孪晶界的形成及其在应变协调中的作用尚缺乏深入理解。
本工作选取工业级TA31近α钛合金，系统研究其在双轴拉伸下的力学响应与微观变形机制，尤其聚焦于首次发现的高能{13‾41}孪晶及其引发的层级孪晶行为，旨在揭示双轴应力状态对孪生活化与应变分配的新调控机制。
2. 材料及试验
研究选用热轧+退火态TA31钛合金板材，厚度约42 mm，具有明显的双相组织（α相占比~96%，平均晶粒尺寸~20 μm）及近基面织构。设计并优化了十字形双轴拉伸试样，采用中心减薄与应力释放孔结构以确保断裂起始于中心区域。双轴试验在立式双轴试验机上进行，沿轧向（RD）与横向（TD）按1:1载荷比加载，同步采用数字图像相关（DIC）获取中心区域全场应变演化。
由于十字形试样不存在均匀标距段，无法直接计算真应力-应变曲线。为此，作者引入“逆分析方法”，即以实测载荷-时间历程作为有限元模型边界条件，将DIC获取的中心应变与模拟结果迭代对比，直至收敛，从而准确反推双轴加载下的真应力-应变响应。这一方法有效克服了双轴试验中应力应变直接测量的难题，具有重要推广价值。
微观表征方面，综合运用EBSD、TKD与TEM，对变形后试样中心区域进行多尺度晶体学分析，并辅以第一性原理DFT计算，定量评估孪晶界能。界面位错与台阶结构则通过HRTEM结合Burgers回路映射进行解析。
3. 结果与讨论
3.1 力学响应：双轴加载下各向异性显著增强
图2（a）真应力-应变曲线对比（单轴vs双轴）
通过逆分析获得的真应力-应变曲线（图2a）显示，双轴加载下RD与TD方向的屈服强度分别提升至约950 MPa与1070 MPa，显著高于单轴拉伸。同时，TD方向的塑性应变仅约2.1%，远低于RD方向的11.7%，表明双轴加载显著放大了材料固有的各向异性。有限元分析与DIC结果一致表明，中心区域应力集中沿TD方向呈带状分布，并伴随应力重新分配，最终诱发沿TD的早期断裂。这种各向异性加剧源于双轴拉伸对横向收缩的抑制，以及更高静水应力分量（三轴度约0.66）对塑性流动的约束。
3.2 孪晶类型演变：从单轴到双轴
图3 单轴与双轴变形后EBSD孪晶分布与统计（{10‾12}绿，{11‾22}红）
单轴拉伸下，{10‾12}拉伸孪晶与{11‾22}压缩孪晶共存，界面密度比约3:1。而双轴加载下，{10‾12}孪晶密度显著增加，但{11‾22}压缩孪晶几乎消失。这一差异并非源于文献中常见的应力路径差异，而是由初始织构决定的施密特因子变化所致。初始基面织构中c轴大多偏离ND方向，在双轴加载下中心区域ND方向以压应力为主，不利于压缩孪晶的活化。
3.3 首次发现高能{13‾41}孪晶及层级孪晶结构
图4 TKD图像及极图（{13‾41}主孪晶 + {10‾12}二、三级孪晶）
在双轴试样断口附近的高应变区域，TKD表征首次识别出{13‾41}孪晶（记为T_H）。该孪晶与基体间存在63°[55‾103]取向关系，界面平直且呈锯齿状形貌。更为重要的是，在该高能孪晶内部观察到清晰的多级孪晶结构：次级{10‾12}孪晶（T′）沿主孪晶界优先形核，随后在次级孪晶内再次形成三级{10‾12}孪晶（T″）。极图与取向分析确认，次级孪晶与基体之间还形成了接近{10‾11}孪晶关系的特殊界面。
TEM选区电子衍射（SAED）进一步验证了{13‾41}界面的镜像对称性，以及次级/三级{10‾12}孪晶间的85°[11‾20]取向关系。HRTEM显示，次级孪晶与基体界面两侧(0001)面夹角约59°，非常接近理论{10‾11}孪晶的57.2°。此外，界面附近存在高密度位错与台阶结构，表明该界面正经历应力驱动的晶体学弛豫。
3.4 DFT与界面位错分析：高能界面的本质与演化驱动力
表1 DFT计算的不同孪晶界能对比
DFT计算表明，{13‾41}孪晶界能为0.52 J/m²，显著高于{10‾12}（0.29）、{11‾22}（0.37）等常见孪晶，证实其“高能”本质。但双轴加载下约两倍的静水拉应力显著降低了形核能垒，从而允许该高能孪晶活化。
进一步通过HRTEM结合Pond回路映射，解析了{13‾41}界面上的界面位错，确定其伯氏矢量空间方向为<56‾12>，面内分量为<18‾9‾11>。采用广义施密特因子（GSF）评估不同应力状态对界面迁移的驱动效果，发现双轴加载不仅提升了最优变体的GSF上限（从0.43增至0.49），而且实际应力水平远高于单轴，从而为{13‾41}孪晶的形核与扩展提供了更强的力学驱动力。
3.5 界面演化机制：从形核到近{10‾11}稳定界面
图8–9 次级孪晶形核与生长界面的TEM对比及位错分析
统计显示，次级孪晶在形核初期与基体的界面角通常大于60°，而生长后逐渐降至约60°以下，并趋向{10‾11}孪晶的对称几何。HRTEM进一步表明，生长界面已具备接近{10‾11}孪晶的原子排布，且存在b₃/₂型界面位错，可分解为可动的b₂/₂孪晶位错与
图10 晶体学关系示意图及层级孪晶形成过程图解
综合晶体学分析，作者提出以下演化路径：{13‾41}主孪晶在高应力下形核→其高能界面成为{10‾12}次级孪晶的优先形核位点→在十二个{13‾41}变体及六个{10‾12}变体中，仅特定变体组合满足几何匹配与能量最低原则，使得次级孪晶与基体界面逐渐趋近于稳定的{10‾11}孪晶关系。这一“高能界面→选择性形核→应力驱动弛豫→近稳定界面”的路径，是一种全新的应变协调机制，在室温变形中前所未见。
4. 总结
本研究首次系统揭示了近α钛合金在双轴拉伸下的力学响应与微观孪生行为，主要结论如下：
力学行为：通过逆分析方法准确获取双轴真应力-应变曲线，发现双轴加载显著放大各向异性，RD与TD方向屈服强度分别达950 MPa与1070 MPa，且TD方向塑性极低（ε_f ≈ 2.1%）。
孪晶模式转变：单轴拉伸下{10‾12}+{11‾22}孪晶共存；双轴加载下{11‾22}压缩孪晶被抑制，{10‾12}拉伸孪晶密度显著增加，反映织构与应力状态对孪晶选择的强耦合作用。
高能{13‾41}孪晶首次发现：双轴加载下，高静水应力降低形核能垒，使原本不利的{13‾41}孪晶（γ=0.52 J/m²）得以活化，参与局部应变协调。
层级孪晶与界面弛豫新机制：{13‾41}孪晶界优先诱导{10‾12}次级孪晶形核，并通过应力驱动界面演化，使次级孪晶与基体界面逐渐趋近于{10‾11}稳定孪晶界。这一多级孪晶路径为HCP金属在复杂应力状态下的应变协调提供了全新认知。
本工作不仅为双轴加载下钛合金力学行为提供了高质量实验数据与物理机制解释，更启示我们：在复杂应力状态下，材料可能激活传统上被视为“不可能”的变形模式，这对先进结构材料的本构建模、寿命评估与服役可靠性设计具有重要理论意义与工程价值。