一种基于自由能二阶导数判据筛选Cu-Ni-Mn调幅分解合金成分的方法

本发明属于铜镍合金成分设计与组织调控，具体涉及一种基于自由能二阶导数判据筛选cu-ni-mn调幅分解合金成分的方法，尤其涉及一种通过自由能数据库计算cu基fcc固溶体gibbs自由能曲线，并根据自由能二阶导数判据筛选可发生调幅分解的cu-ni-mn合金成分区间的方法。

背景技术：

1、铜镍合金具有良好的低温韧性、耐腐蚀性能、加工性能和组织稳定性，在低温装备、发动机密封件、航空航天、海洋工程及复杂服役环境中具有重要应用价值。尤其是在液氧超低温环境下，特种铜镍合金不仅需要具有较高强度，还需要保持良好的塑性、冲击韧性和组织稳定性。因此，开发适用于液氧超低温环境的高强韧特种铜镍合金，是铜基合金材料设计的重要方向。

2、传统cu-ni系合金主要依赖固溶强化提高强度，但单纯固溶强化的强化能力有限。若通过大量添加强化元素或引入粗大第二相提高强度，又可能导致合金低温韧性下降，并增加裂纹萌生风险。对于液氧超低温服役环境而言，材料内部粗大析出相、连续晶界析出相或明显成分偏聚均可能成为低温损伤源。因此，如何在不引入大量粗大脆性相的前提下提高铜镍合金强度，是特种铜镍合金设计中的关键问题。

3、调幅分解是一种由过饱和固溶体在热力学不稳定区内发生的连续成分起伏过程。与传统形核-长大型析出不同，调幅分解通常不需要明显形核势垒，可形成尺度细小、分布连续、界面弥散的成分调制组织。该类组织可通过成分起伏、弹性应变场和细尺度界面调制效应提高合金强度，同时避免粗大第二相对低温韧性产生显著不利影响。

4、从热力学角度看，调幅分解能否发生取决于固溶体自由能曲线对微小成分扰动的响应。当固溶体处于稳定区时，微小成分起伏会导致体系自由能升高，成分扰动不能自发放大；当固溶体处于调幅分解热力学不稳定区时，微小成分起伏能够降低体系自由能，成分扰动可自发放大，进而形成连续成分调制组织。

5、对于设定成分扫描路径下的固溶体自由能曲线，可通过gibbs自由能对成分变量的二阶导数判断其稳定性。当自由能曲线二阶导数大于0时，体系处于成分扰动稳定区；当自由能曲线二阶导数小于0时，体系处于调幅分解热力学不稳定区，对应成分区间具有发生调幅分解的可能性。

6、然而，现有cu-ni-mn调幅分解合金的设计方法存在以下技术缺陷：

7、第一，现有技术多依赖经验试错和大量热处理实验，通常需要在宽成分范围内制备数十种合金，再通过显微组织观察和性能测试筛选可发生调幅分解的成分及热处理制度。该方法实验周期长、样品消耗大，导致合金开发周期长达数月，研发成本居高不下。

8、第二，现有技术缺乏明确的物理判据指导成分设计，多以具体合金成分或热处理工艺为主，难以精确定量给出调幅分解不稳定区的起止成分点，也无法明确区分稳定区与不稳定区，导致成分设计存在较大模糊性和不确定性。

9、第三，现有技术通常仅针对特定温度或特定成分路径进行设计，缺乏对不同热处理温度或服役温度的灵活适配能力，也难以适应不同成分扫描路径的设计需求。同时，现有技术多依赖特定商业热力学软件，技术门槛高且扩展性受限，难以推广至其他铜镍基合金体系。

10、第四，现有技术在建立热力学数据库时数据来源单一，计算精度受限，且往往需要复杂的专业软件支持，普通研发人员难以快速掌握和实施，技术普及难度大。

11、第五，现有技术缺乏理论预测与实验验证的闭环设计流程，通常先进行大量实验再总结经验，难以在实验制备前快速判断合金是否具备发生调幅分解的热力学条件，预测结果的可靠性缺乏有效验证。

12、第六，现有技术虽然公开了部分cu-ni-mn合金的成分和制备方法，但未能给出经过系统验证的、可直接指导生产的工艺参数窗口，如时效温度范围对调幅分解组织形成能力的影响缺乏系统研究，难以满足工程化应用对工艺可操作性的要求。

13、因此，亟需建立一种开发周期短、研发成本低、判据物理意义明确、成分边界可量化、方法灵活通用、工艺适配性强、数据来源多元、技术门槛低、预测与实验闭环验证、工艺参数具体可行的cu-ni-mn调幅分解合金成分筛选方法，用于指导cu-ni-mn及其他cu-ni-x调幅分解强化型铜镍合金的开发。

技术实现思路

1、一种基于自由能二阶导数判据筛选cu-ni-mn调幅分解合金成分的方法，包括以下步骤：

2、步骤1：建立或调用cu-ni-mn体系自由能数据库

3、建立或调用用于计算cu基fcc固溶体在不同cu、ni和mn含量条件下gibbs自由能的数据库。所述自由能数据库包括以下参数中的一种或多种：cu、ni、mn纯组元的gibbs自由能参数；cu-ni二元相互作用参数；cu-mn二元相互作用参数；ni-mn二元相互作用参数；以及cu-ni-mn三元相互作用参数。

4、所述自由能数据库由以下数据中的一种或多种构建或修正得到：实验热力学数据、相平衡数据、热分析数据、第一性原理计算数据、团簇展开计算数据、声子振动自由能计算数据。

5、基于所述自由能数据库，cu基fcc固溶体的gibbs自由能按式（1）计算：

6、gfcc = xcu*gcufcc + xni*gnifcc + xmn*gmnfcc + r*t*[xcu*ln(xcu) + xni*ln(xni) +xmn*ln(xmn)] + gexcess                      （1）

7、其中，xcu、xni和xmn分别为cu、ni和mn的摩尔分数，gcufcc、gnifcc和gmnfcc分别为cu、ni和mn纯组元在fcc结构下的gibbs自由能，r为气体常数，t为计算温度，gexcess为过量gibbs自由能。

8、过量gibbs自由能gexcess按式（2）计算：

9、gexcess = xcu*xni*lcu,ni + xcu*xmn*lcu,mn + xni*xmn*lni,mn + xcu* xni* xmn*lcu,ni,mn                         (2)

10、式（2）中，lcu,ni、lcu,mn、lni,mn和lcu,ni,mn分别为cu-ni、cu-mn、ni-mn二元相互作用参数和cu-ni-mn三元相互作用参数；各相互作用参数可为温度相关参数，也可由多项式形式表示。

11、步骤2：设定计算温度并沿预设成分扫描路径计算gibbs自由能曲线

12、根据目标热处理温度或目标服役温度设定计算温度t（t 为热力学温度，单位：k）。所述设定温度根据目标热处理温度或目标服役温度确定，示例性为450℃，但不局限于450℃。

13、在设定的计算温度t下，沿预设cu-ni-mn成分扫描路径，计算cu基fcc固溶体的gibbs自由能，获得gibbs自由能随成分变量变化的自由能曲线g=g(x,t)，其中x为预设成分扫描路径上的成分变量。

14、所述预设成分扫描路径为以下任一种：固定mn含量扫描ni含量；固定ni含量扫描mn含量；固定ni/mn比例扫描总合金元素含量；或者在cu-ni-mn三元成分空间中选取的任一连续成分路径。

15、步骤3：计算自由能曲线的二阶导数

16、计算步骤2所得自由能曲线g=g(x,t)对成分变量x的二阶导数d2g/dx2。二阶导数可通过解析求导方式计算，或通过有限差分方式近似计算。

17、采用有限差分计算时，按式（3）进行计算：

18、d2g/dx2 ≈ [g(x+δx,t) - 2g(x,t) + g(x-δx,t)] / (δx)2     (3)

19、式（3）中，δx为预设成分扫描路径上的成分步长。

20、步骤4：根据二阶导数判据判定调幅分解热力学不稳定区

21、当步骤3计算得到的二阶导数满足d2g/dx2＜0时，判定该成分点或相邻成分区间处于调幅分解热力学不稳定区，具有发生调幅分解的可能性；当d2g/dx2＞0时，判定该成分区间处于成分扰动稳定区。

22、步骤5：筛选得到cu-ni-mn调幅分解合金成分区间

23、根据步骤4判定的调幅分解热力学不稳定区，筛选得到具有调幅分解可能性的cu-ni-mn合金成分区间。

24、优选地，所述cu-ni-mn调幅分解合金成分区间按质量百分比计包括：ni 18～22%、mn 18～22%，余量为cu及不可避免杂质。

25、步骤6：选取候选合金成分并通过实验验证调幅分解组织形成能力

26、从步骤5筛选得到的cu-ni-mn调幅分解合金成分区间中选取候选合金成分，并通过熔炼、固溶处理和时效处理验证其调幅分解组织形成能力。

27、进一步优选地，从所述cu-ni-mn调幅分解合金成分区间中选取的候选合金成分为cu-20ni-20mn，其中ni含量为20wt.%，mn含量为20wt.%，余量为cu及不可避免杂质。

28、所述熔炼方式为以下任一种：真空感应熔炼、真空电弧熔炼、保护气氛熔炼或感应熔炼。

29、所述固溶处理温度为550～700℃，保温时间为1～5h，固溶处理后进行水淬。

30、所述时效处理温度为400～500℃，时效处理后进行空冷。

31、所述调幅分解组织为富cu区与富ni-mn区交替分布的条幅状、波纹状或连续成分调制组织。

32、验证调幅分解组织形成能力采用以下一种或多种方法：扫描电子显微镜、透射电子显微镜、能谱分析、电子探针、原子探针。

33、相对于现有技术，本专利申请具有如下有益效果：

34、1. 开发周期大幅缩短，研发成本显著降低

35、本发明摒弃了传统试错法需要在宽成分范围内制备数十种合金、进行大量热处理实验的繁琐流程，通过自由能二阶导数判据可在实验前快速预测调幅分解成分区间，将成分筛选周期从数月缩短至数天。同时，本发明大幅减少了高纯金属原料的消耗及熔炼、热处理等实验成本。仅通过一次针对cu-20ni-20mn合金的验证实验即成功证实了预测区间的准确性，避免了大量无效实验的重复开展。

36、2. 预测判据物理意义明确，成分边界可量化

37、本发明基于cahn-hilliard调幅分解理论，采用自由能对成分的二阶导数小于零作为核心判据，物理内涵清晰明确。该判据能够精确定量给出调幅分解不稳定区的起止成分点，450℃下ni含量为15至25个重量百分比的区间，并明确区分稳定区与不稳定区，消除了传统经验设计中的模糊性，为合金成分设计提供了精确的热力学边界。

38、3. 方法灵活通用，工艺适配性强

39、本发明的计算温度可根据实际热处理或服役温度灵活设定。实施例2验证了本发明在400℃、450℃、500℃及550℃下的预测能力，结果表明随着温度升高，调幅分解不稳定区间向低镍含量方向移动，当温度升高至550℃时已无调幅分解发生，该方法可有效指导不同温度下的工艺窗口设计。同时，本发明支持固定某元素含量、固定元素比例等多种成分扫描路径。无论采用固定锰含量、固定镍含量还是固定镍锰比例路径，预测结果均高度自洽，均指向镍和锰各为15至25个重量百分比的成分区间，说明该方法对不同设计路径具有稳健性。此外，本发明不依赖于特定商业热力学软件，仅需建立或调用自由能数据库即可实施，可方便地推广至其他铜镍基合金体系的调幅分解强化型合金设计。

40、4. 数据来源多元，技术门槛低

41、本发明建立自由能数据库的途径多样，可采用实验热力学数据、相平衡数据、第一性原理计算、团簇展开计算等多种来源，不依赖单一数据源，有效保证了计算的准确性和可靠性。在有限差分计算中，差分步长可根据精度需求灵活调整，如取零点五个重量百分比可在计算精度与效率之间取得良好平衡。配合常规编程手段如python或matlab即可实现全部计算，显著降低了技术门槛和实施成本。

42、5. 预测与实验闭环验证，可靠性高

43、本发明不仅提供理论预测，还通过具体的候选合金的熔炼、固溶、时效及组织表征进行闭环验证。对cu-20ni-20mn合金的扫描电子显微镜观察显示形成了条宽约50至150纳米的条幅状组织，能谱分析证实了富铜区与富镍锰区的交替分布，实验结果与自由能二阶导数判据的预测完全吻合。这种理论预测与实验验证相结合的闭环模式，确保了本发明方法的科学性与实用性。

44、6. 工艺参数具体可行，指导价值高

45、本发明给出了经过实验验证的具体工艺参数，包括固溶处理温度为600摄氏度、保温3小时后水淬，时效处理温度为450摄氏度、保温4小时后空冷等。这些参数可直接用于指导生产实践。在400至500摄氏度的时效温度窗口内均可形成调幅分解组织，这为通过调控时效温度来调节组织尺度和优化力学性能提供了明确的操作空间。