基于Abaqus的非晶颗粒铜晶体塑性仿真

关键词：Abaqus；晶体塑性；铜基复合材料；增材制造；疲劳

一、文章简要介绍

金属强化的老路子是往里撒颗粒挡位错，但晶体颗粒硬碰硬，界面处应力集中，裂纹从这冒头，强度上去塑性就塌。这篇Nature Communications换了个思路（DOI：10.1038/s41467-025-58893-1）：激光粉末床熔融打印时把B4C微颗粒原位熔碎成平均47纳米的非晶颗粒，散进铜基体。非晶没有晶格，位错打进去直接被吸收。

L-PBF原位纳米碎裂原理与微观结构表征

图1 激光粉末床熔融原位纳米碎裂原理与结构表征

铜晶粒也细，平均0.39微米，颗粒含量越高晶粒越细。这套复合材料强度破1 GPa、延伸率约10%，比晶体颗粒版本翻倍；10的7次方周次下疲劳极限还有抗拉强度的七成，压过纳米晶和梯度铜。

不同颗粒含量下的EBSD晶粒尺寸演变

图2 颗粒含量越高晶粒越细：EBSD晶粒尺寸演变

二、仿真流程与建模细节详解

颗粒怎么碎成非晶的，靠熔池仿真还原：流体仿真追踪B4C在熔池里的熔化与碎裂，网格0.025微米，碎到比网格还小的液滴交给离散颗粒方法接管；光线追迹算出混合粉末床的激光吸收率。能量密度从86.9加到180.6 J/mm³，颗粒从5微米碎到纳米级，快速凝固把非晶结构冻住。

能量密度驱动的颗粒碎裂机制与熔池仿真

图3 体积能量密度驱动颗粒碎裂与熔池仿真

力学这边建了双相晶体塑性模型：铜基体用率相关有限应变晶体塑性，12个{111}<110>滑移系，滑移率按幂律走，潜硬化矩阵非对角元取1.4；B4C颗粒用率无关J2塑性加各向同性硬化。整套本构写成VUMAT子程序在Abaqus/Explicit里求解，参数由原位同步辐射XRD的晶格应变标定，仿真曲线与实验对得上。分子动力学补一刀：位错被非晶颗粒吸收，遇晶体颗粒则先钉扎再绕过。

CPFE应力分布、晶格应变对比与MD位错交互

图4 CPFE应力与晶格应变对比及MD位错交互

三、关键结果与发现

应力配分是关键：颗粒真应力扛到1 GPa以上，替基体分担大部分载荷，两相变形协调。疲劳这关也过了：应力比0.1的拉拉疲劳，10的7次方周次疲劳极限超过抗拉强度的70%，高周段甩开纳米晶铜和梯度铜。

拉伸曲线、性能对比与疲劳寿命

图5 拉伸性能、横向对比与疲劳寿命曲线

机理看TEM：晶体颗粒周边位错塞积，非晶颗粒周围没有位错堆积，变形更均匀。断口也印证：非晶版本韧窝深，晶体版本颗粒界面有微裂纹。非晶还得扛热：退火到673 K硬度不掉，一旦晶化就明显软化。

位错与晶体/非晶颗粒交互的TEM观察

图6 位错在晶体颗粒处钉扎、在非晶颗粒处无塞积

断口对比与退火晶化对性能的影响

图7 断口形貌对比与退火晶化影响

四、我们提供的服务

晶体塑性与疲劳仿真我们按项目交付：Abaqus UMAT/VUMAT本构开发、晶体塑性RVE建模与参数标定、原位XRD与EBSD对标、疲劳寿命预测，也接增材制造工艺仿真。

如有需要，欢迎通过公众号"320科技工作室"联系我们。