GND直接关联各向同性硬化和运动硬化
参考文献：《Indentation size effects in spherical nanoindentation an alyzed by experiment and non-local crystal plasticity》 在微尺度金属塑性变形中，一个长期存在却又难以在传统模型中准确描述的现象是尺寸效应：随着构件尺度、晶粒尺寸或应变梯度的增大，材料表现出明显高于宏观预测的强度。这一现象的物理根源，普遍认为与几何必要位错（Geometrically Necessary Dislocations, GNDs）密切相关。 该文章围绕这一问题，构建了一套基于应变梯度的晶体塑性理论框架，通过在连续介质尺度中显式引入 GND 的产生、演化及其力学效应，实现了对微尺度强化行为的统一描述。 作者所提出与以往经典的GND模型的差异在于，文章系统地考虑了 GND 对塑性行为的两种关键影响： （1）局部强化效应（短程作用）：GND 与统计存储位错（SSD）共同构成材料内部的位错网络。模型采用 Taylor 型关系，将 GND 密度引入滑移阻力中，从而自然解释了应变梯度增大时材料强度提升的现象。（2）长程内应力效应（非局部作用）：由于 GND 本质上是具有空间分布的位错结构，其不均匀性会产生长程弹性应力场。模型通过引入高阶梯度项，等效构造位错线段，并利用解析位错力学解计算其诱导的内应力。这部分应力作为附加分解剪应力或背应力，直接参与滑移驱动力的计算。作者将上述两类应力直接引入经典的流动方程：由于同时计算了塑性变形梯度的一阶和二阶梯度，因此数值实现难度较大，作者建议使用基于谱方法的中心差分计算对应的梯度。将有限元高斯点上的塑性变形梯度映射至规则网格，在频域中高效计算其空间导数与平滑梯度，从而稳定、精确地获得 GND 相关量。这种方法在保证数值效率的同时，有效抑制了高阶梯度带来的数值噪声问题。作者使用该模型结合AFM和TEM研究了球形压痕的尺寸效应问题.作者的研究表明：（1）在小压痕载荷/小压痕深度下，材料的塑性区尺度受几何强烈约束，GND 密度较高；大压痕载荷下，塑性区持续扩展，SSD 累积占主导。（2）对于球形压头，接触曲率在加载过程中基本保持不变；这意味着由几何约束产生的应变梯度尺度是近似恒定的；因此，GND 密度及其引起的各向同性强化项在不同压痕深度下变化不大。相反，随着压痕加深：塑性应变不断累积；SSD 引起的传统加工硬化持续增强。（3）仅考虑 GND 引起的各向同性强化项，即可清晰再现 ISE；仅考虑 GND 的运动学背应力项，或完全忽略 GND，则尺寸效应显著减弱甚至消失。ISE 的主导机制并非位错长程内应力本身，而是 GND 对滑移阻力的统计强化作用。（4）非局部晶体塑性模型可以同时拟合：载荷–位移曲线，残余压痕形貌，硬度随载荷变化规律；计算得到的 GND/SSD 分布趋势，与 TEM 观察到的位错结构在定性上高度一致。这表明该模型不仅是“可用的”，而且在物理解释层面具有可信度。作者提出的这类方案便于数值实现，对于异质性结构诱导额外强化的分析也有一定的参考意义，使用作者提到的这种谱方法计算GND测试显示计算结果和FEM的MLS分布接近，数值大小会高一些，同时也会更平滑一些。使用FFT方案计算的案例如下：初始包含200个晶粒的多晶模型（拉伸变形2.0%，单元类型CPE4R）等效应力分布：变形结束后的GND分布： GND一阶效应的各向同性阻力：GND二阶效应的背应力 来源：我的博士日记
