镍基超合金因其优异的高温强度、抗蠕变性能以及微观结构稳定性，被广泛用于航空航天系统及先进超临界发电厂的关键部件制造[[1], [2], [3]]。这类合金通常含有一定比例的孪晶界，其变形行为与孪晶界的存在及演化密切相关。一方面，孪晶界可作为位错吸收或传递的通道，从而促进局部塑性变形[[4], [5], [6]]；另一方面，孪晶界也可能阻碍位错运动，导致位错在界面处堆积，进而提高材料的强度[7,8]。此外，位错堆积引发的应力集中还可能进一步诱发晶界裂纹的产生[9,10]。因此，深入理解孪晶界的响应对于提升这类合金的机械性能具有重要意义。
随着微观结构表征技术、数值方法以及基于微观结构的建模方法的不断发展，应变梯度晶体塑性有限元模型已被广泛用于研究多晶材料中的异质变形现象[[11], [12], [13], [14], [15], [16]]。在镍基合金研究中，该模型被广泛应用于分析晶格旋转、滑移系激活以及晶间变形的协调性。Chen等人[17]将原位高温拉伸电子背散射衍射技术与应变梯度晶体塑性有限元模型相结合，研究了因康镍718中常规晶界与孪晶界处的滑移传递机制。他们发现，晶界处的滑移传递不仅取决于晶界错位程度，还与相邻滑移系的几何匹配性以及局部应力状态有关。Hestroffer等人[18]通过三维EBSD、高分辨率数字图像相关技术以及三维晶体塑性建模，确定了镍基合金中晶界三重结点及孪晶界附近区域是滑移局部化的关键位置。Texier等人[19]进一步表明，因康镍718中的显著应变局部化现象往往发生在孪晶界或其邻近区域。其他研究[[20], [21], [22], [23]]也指出，孪晶界附近的局部响应与相邻晶粒的取向、滑移系激活情况以及各滑移系之间的几何匹配性密切相关。这些研究都表明，应变梯度晶体塑性有限元模型能够建立局部晶界尺度响应与宏观力学性能之间的关联，是研究镍基合金异质变形的重要工具。
近年来，将应变梯度晶体塑性有限元模型的模拟结果转换为与实验EBSD兼容的数据格式，为模拟结果与实验数据的直接对比提供了一种新方法。Song等人[24]在应变梯度晶体塑性有限元模型中加入了真实的晶粒形态、晶体学取向以及晶粒尺寸效应，并将模拟得到的晶体学取向数据转换为虚拟EBSD数据，从而预测了铝合金焊接接头在拉伸变形过程中的微观结构演化规律。与传统方法相比，虚拟EBSD能够在与实验EBSD相同的数据框架下，直接分析晶粒旋转、错位以及晶粒尺度的异质变形现象。不过，这种方法尚未被用于研究镍铬钴基合金中孪晶界附近的异质变形与滑移传递现象。
虽然传统的应变梯度晶体塑性有限元模型能够描述由晶粒取向差异引起的滑移系激活与晶格旋转现象，但其本构模型通常缺乏明确的内在材料长度尺度，因此难以准确描述塑性应变梯度、几何上必然产生的位错积累以及晶界附近的额外强化效应[25,26]。离散位错动力学模型虽然能够直接描述单个位错与晶界之间的相互作用，但由于计算成本过高，难以同时兼顾真实的多晶微观结构与计算效率[27,28]。近年来发展的应变梯度晶体塑性有限元模型则引入了滑移梯度以及基于GND的强化机制，因此能够有效描述取向梯度、位错积累以及晶界附近的异质强化现象[[29], [30], [31]]。尽管如此，将应变梯度晶体塑性有限元模型、真实微观结构以及虚拟EBSD重建技术相结合的研究仍然较少，针对异质变形以及孪晶界处的滑移传递与阻塞现象的多尺度分析也还不够充分。
为此，本研究以通过因康镍740的成分优化所得到的镍铬钴基合金为研究对象，结合EBSD表征与单轴拉伸试验，基于真实的初始微观结构建立了准三维应变梯度晶体塑性有限元模型。该模型引入了基于GND的强化机制，用于描述塑性应变梯度、位错积累以及晶界与孪晶界附近的局部强化效应。在此基础上，利用MATLAB开发的重建脚本将模拟得到的晶体学取向与晶格旋转信息转换为虚拟EBSD数据，使得模拟结果能够在与EBSD兼容的框架内进行直观展示与分析。与传统应变梯度晶体塑性有限元模型相比，该方法能够更有效地捕捉由应变梯度以及晶界附近的GND积累所引发的异质强化效应；与离散位错动力学模型相比，则能够高效地模拟具有真实晶粒形态的多晶微观结构。通过该方法，系统分析了拉伸变形过程中滑移系激活、滑移传递以及孪晶界处的滑移阻塞现象，进一步阐明了孪晶界在调控合金强度与延性方面的协同作用。这些研究结果为镍基合金的晶界设计以及微观结构优化提供了理论支撑。
