《Mechanics of Materials》：Quantitative prediction of size effects using an advanced micromorphic approach
编辑推荐：
本研究针对高阶应变梯度塑性理论在预测小尺度材料尺寸效应时存在的能耗处理不当、尺寸依赖性高估等问题，采用一种基于多重级数分解的塑性变形微态模型，对受限压缩薄层、L梁单晶弯曲-扭转以及循环剪切加载等多种复杂条件下的尺寸效应进行了系统的有限元模拟与定量评估。结果表明，该模型能精确预测取向相关的屈服响应，并与离散位错动力学模拟及实验观测数据吻合良好，显著提升了小尺度晶体材料尺寸效应的定量建模能力。
随着微电子、生物医疗器件、涂层技术和航空航天等领域的迅猛发展，微型化器件的应用日益广泛。然而，一个令人头疼的现象也随之而来：当材料的尺寸缩小到微米甚至纳米尺度时，其力学行为会与宏观“大块头”截然不同，通常表现出显著的“尺寸效应”，即材料强度会随着特征尺寸的减小而反常地增高。比如，一根微米级的金属丝，其屈服强度可能远高于同种材料的宏观块体。这种现象通常归因于非均匀塑性变形产生的几何必需位错（GNDs），而传统的塑性理论由于缺乏内禀长度尺度，对此束手无策。
为了捕捉这种尺寸效应，应变梯度塑性（SGP）理论应运而生。其中，高阶SGP理论通过引入额外的自由度和平衡方程，成为目前最常用的方法。但现有理论在预测实践中仍面临两大核心挑战：一是对“耗散”处理的争议，高阶过程带来的能量损失如何准确建模以避免非物理的“弹性间隙”现象；二是对尺寸依赖性强度的预测普遍偏强。多数SGP模型预测的屈服应力与特征尺寸的尺度因子（s）在1到2之间，而实验和离散位错动力学（DDD）模拟却一致显示尺度因子通常小于或等于1，著名的霍尔-佩奇（Hall-Petch）关系给出的指数是0.5。这种预测与观测的鸿沟，促使研究者们不断寻求新的理论框架。
微态方法为此提供了一条有前景的路径。它通过引入与物理塑性场弱耦合的辅助（微态）变量，并考虑其梯度而非物理塑性场的梯度，为灵活调控本构方程中的梯度贡献及最终的尺寸依赖性提供了可能。最近，一种基于塑性变形、融入了多重高阶运动学项级数分解的新型微态模型被开发出来，它能够通过罚参数约束微态变量，灵活控制尺度效应，并同时包含了高阶能量的和耗散的作用，从而避免了弹性间隙。
为了系统评估这一先进模型的定量预测能力，Anjan Mukherjee, Mohamed Jebahi, Mohamed Abatour, Samuel Forest等研究人员在《Mechanics of Materials》期刊上发表了一项研究。他们通过理论分析、有限元模拟，并将结果与实验观测及DDD预测进行系统对比，全面考察了该模型在多种载荷条件和试样几何下的表现。研究首先通过一个理想化的受约束剪切层问题进行解析研究，阐明了与高阶能量和耗散效应相关的尺度效应。随后，将模型应用于预测受限压缩下薄金属层的取向依赖性响应。有限元模拟显示，对于垂直于加载方向的层，模型预测了显著的一阶和高阶尺寸效应，而对于45°倾角的层，尺寸效应则可忽略，这与文献中的实验观察一致。接着，模型被用于捕捉承受弯曲-扭转载荷的L梁单晶试样中的尺寸效应，通过同时纳入能量和耗散高阶贡献，模型预测的尺度指数（0.36）与实验观测值（0.38）吻合良好。最后，模型被用来复现DDD在循环剪切载荷下获得的复杂尺寸依赖响应，预测的尺度指数为0.26，与近期DDD研究报告的范围（0.2, 0.3）完全一致。这些结果充分证明了所采用模型在多种复杂场景下强大的定量预测能力。
本研究主要采用了三种关键技术方法：一是理论建模与解析推导，基于塑性变形微态模型，对受约束剪切层问题进行解析，探究高阶能量和耗散项的尺度效应；二是有限元模拟，利用商业软件Abaqus通过用户子程序UEL和UMAT实现了该微态模型，模拟了受限压缩、L梁弯曲-扭转等复杂工况；三是系统性的对比验证，将模型预测结果与来自已发表文献的实验数据（如Mu等人关于压缩的实验，Zhang等人关于L梁弯扭的实验）以及离散位错动力学（DDD）模拟结果（如Amouzou-adoun等人的循环剪切研究）进行定量对比，以评估模型的准确性。
3. 约束剪切层分析
通过一个一维的理想化受约束剪切问题，研究人员对模型的尺度效应进行了解析研究。他们发现，屈服应力与层厚（h）的尺度关系（σy∝ h-s）中，尺度因子s强烈依赖于高阶能量和耗散参数的比值。当高阶能量主导时，模型会预测较强的尺寸效应（s≈1）；而当耗散贡献显著时，尺寸效应会减弱（s<1）。这一解析结果为理解和调控模型预测的尺寸依赖性强度提供了清晰的物理图像。
4. 实验和DDD结果的数值预测
- •
4.1 受限压缩中的取向依赖尺寸效应：研究模拟了不同晶体取向的薄金属层在刚性模具间的压缩。结果显示，对于晶体取向使滑移面垂直于加载方向的层（即[111]取向），模型成功预测了显著的尺寸强化效应，这来自于一阶和高阶贡献的共同作用。而对于滑移面与加载方向成45°的层（即[112]取向），模型预测的尺寸效应可忽略。这一预测与Mu等人的实验结果定性一致，而传统的高阶SGP理论在不进行特殊边界条件修正的情况下难以复现这种取向依赖性。
- •
4.2 L梁弯曲-扭转中的尺寸效应：研究人员应用模型模拟了带有钝化（passivated）和未钝化（unpassivated）表面的L形单晶梁在弯曲-扭转载荷下的响应。通过引入单一的高阶运动学分解并纳入耗散贡献，模型成功地复现了Zhang等人在实验中观察到的屈服应力尺寸依赖性。对于未钝化的梁，模型预测的尺度因子（s=0.36）与实验值（s=0.38）非常接近。模拟还表明，未钝化梁的强化主要源于边界处的GNDs积累，而钝化梁的额外强化则归因于表面约束对位错运动的阻碍。
- •
4.3 循环剪切下的尺寸效应：最后，研究团队将模型应用于Amouzou-adoun等人通过DDD模拟获得的循环剪切载荷数据。在合理标定材料参数后，该微态模型成功地复现了DDD预测的循环应力-应变响应及随试样尺寸变化的屈服应力。模型预测的尺度因子为0.26，完全落在近期DDD研究所报告的范围（0.2, 0.3）之内。这证明了该模型即使在涉及位错重排和湮灭的复杂循环载荷条件下，仍能定量捕捉尺寸效应。
本研究通过全面的解析、数值和对比分析，证实了基于塑性变形和多重级数分解的先进微态模型在预测小尺度材料尺寸效应方面具有强大的定量能力。该模型不仅能够复现传统SGP理论难以描述的取向相关尺寸效应，还能在弯曲-扭转、循环剪切等多种复杂加载条件下，与实验和离散位错动力学模拟结果取得良好的一致。其成功的关键在于：1）通过微态变量和罚参数实现了对梯度效应的灵活控制，缓解了传统SGP对尺寸效应的高估问题；2）采用高阶运动学分解（而非应力分解）来纳入耗散，避免了非物理的“弹性间隙”；3）基于广义幂律缺陷能势计算高阶应力，能够模拟复杂的非线性运动学响应。这项工作凸显了该松弛微态实现方案在定量建模小尺度晶体材料尺寸效应方面的有效性和广阔应用前景，为微纳米器件的力学设计和可靠性评估提供了更精准的理论工具。
