CN104461691B - 界面层互扩散优化时效成形/扩散复合组织的相场模拟方法 - Google Patents

界面层互扩散优化时效成形/扩散复合组织的相场模拟方法 Download PDF

Info

Publication number
CN104461691B
CN104461691B CN201410744317.6A CN201410744317A CN104461691B CN 104461691 B CN104461691 B CN 104461691B CN 201410744317 A CN201410744317 A CN 201410744317A CN 104461691 B CN104461691 B CN 104461691B
Authority
CN
China
Prior art keywords
phase
interface
diffusion
atomic
precipitation
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201410744317.6A
Other languages
English (en)
Other versions
CN104461691A (zh
Inventor
张静
陈铮
王永欣
卢艳丽
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Northwestern Polytechnical University
Original Assignee
Northwestern Polytechnical University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN201410744317.6A priority Critical patent/CN104461691B/zh
Publication of CN104461691A publication Critical patent/CN104461691A/zh
Application granted granted Critical
Publication of CN104461691B publication Critical patent/CN104461691B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

一种界面层互扩散优化时效成形/扩散复合组织的相场模拟方法,该方法针对时效成形/扩散复合界面层组织取向粗化、沿晶无沉淀带等界面不连续性现象,分析取向性决定的沉淀相取向粗化形貌、取向演化机制、沿晶无沉淀带形成规律,探明界面层应力对原子穿越结合面扩散的促进作用,得出界面层上应力作用引发的沉淀行为变化,提出消除界面层不连续沉淀相的原理,优化组织。本发明用于消除时效成形/扩散复合过程应力诱发界面层动态不连续沉淀相,以获得弥散分布、交叉排列的动态沉淀组织,提升材料性能。

Description

界面层互扩散优化时效成形/扩散复合组织的相场模拟方法
技术领域
本发明涉及塑性成形领域的计算机模拟和组织优化,具体涉及界面层互扩散优化时效成形/扩散复合组织的相场模拟方法。
背景技术
时效成形/扩散复合是在时效成形基础上,受层状复合材料启发发展而来的一种复合材料制备和构件成形一体化技术。时效成形,又称蠕变时效成形(creep ageforming)、热压罐成形,与应力时效,蠕变时效相通,本质是时效蠕变变形机制和时效硬化提升强度同时发生的过程,能同时获得动态析出组织和实现精确成形。但是两个过程结合,二者相互影响,蠕变变形应力影响抑制沉淀相形核,但加速长大和粗化、诱发沉淀相定向排列,为组织各向同性的控制带来了新的问题,其中以沉淀相粗大、沿晶无沉淀带形成,致使界面不连续性对性能损害最为严重。
沉淀组织的稳定性、取向、分布、颗粒度是决定材料性能的关键要素。时效成形的高温下应力时效作用,原子取向扩散致使沉淀相取向排列、应力诱导沉淀相粗化严重、出现沿晶无沉淀相带(PFZ,phase free zone),致使组织不均匀,沉淀相不连续分布。机械性能显示应力时效和无应力时效相比拉伸性能稍低,屈服强度、延伸率明显降低。实验观察表明,沉淀相在应力时效下定向排列,其取向与拉应力方向平行,与压应力方向垂直;高应力区域,颗粒粗化严重。时效成形的高温时效使溶质原子扩散加剧,沉淀相在位错处异相形核和晶界处不连续分布,加剧了沉淀相析出速度和过时效。温度越高、时效时间越长,晶粒越粗大,PFZ越明显,沉淀相不连续分布越明显。PFZ的出现原因与高温时效沉淀相的长大和粗化消耗大量溶质原子导致沉淀相相周围基体中溶质原子贫化有关。软质的PFZ也是性能降低的一个重要原因。
发明内容
本发明的目的在于提供了界面层互扩散优化时效成形/扩散复合组织的相场模拟方法,它通过追踪应力时效全过程组织从无序固溶态到沉淀相析出、长大、粗化过程的微观形貌演化,解析沉淀相析出机制、取向粗化机制、晶间无沉淀带形成规律,探明两相弹性模量差、沉淀相结构、相界结构等微观相界失配应力场和时效成形宏观应力场诱发原子穿越结合面扩散的促进作用,得出结合层上沉淀行为变化规律,达到消除结合面不连续沉淀相的目的。
本发明的技术解决方案是:
一种界面层互扩散优化时效成形/扩散复合组织的相场模拟方法,其特殊之处在于,该方法是:
1)用原子尺度的相场扩散方程表征时效成形/扩散复合全过程界面层沉淀相的形核、长大、粗化过程的形貌演化;该形貌演化通过把原子占位的数据信息转化为微结构演化的图形信息来实现;
2)分析上述图形信息,观察界面层沉淀组织形貌演化,用原子占位定量表征形貌演化过程,明晰不同纳观相界失配应力场和宏观弹性应力场作用下原子取向扩散规律,得出原子取向扩散导致的沉淀相取向粗化和沿晶无沉淀带等界面层不连续性,明晰应力场对界面层不连续性的影响规律;
3)分析上述形貌演化过程中同相沉淀相、异相沉淀相间界面结构及其演化,用原子占位定量表征不同界面结构边界处、界面两侧的演化,分析原子簇聚和贫化机理,得到沉淀相相界稳定性和迁移规律,明晰沉淀相长大、取向粗化以及沿晶无沉淀带形成机制;
4)结合2)形貌演化和3)界面演化,加以原子取向扩散定量分析,得到影响界面层不连续性的因素,得出消除结合面不连续沉淀相的规律。
上述界面层互扩散优化时效成形/扩散复合组织的相场模拟方法,其特征在于,该方法具体包括:
分为两个步骤进行,第一,数据信息转化为图形信息;第二,分析图形信息得出有益于消除界面层不连续性的结果;
第一,数据信息转化为图形信息:
(1)依据微观相场理论编制方程,并求解;
(2)设置初始的变量值,环境变量:温度、成分、外加应力;固有参量:晶格常数,弹性常数,原子间相互作用势,热起伏;计算参量:格点数,迭代步数,迭代步长;
(3)方程求解过程是在倒易空间下进行的,并间隔一定步数转为正空间做一次判定,判定条件为:占位几率在0~1之间,程序继续执行;占位几率在0~1区间之外,计算终止,返回修改参数;
(4)计算结束得到一组原子占位几率值数值,把原子占位几率转化为图形信息;
第二,分析图形信息得出有益于消除界面层不连续性的结果:
(1)绘制不同温度、成分、应力条件下的形貌演化图,得到环境变量对沉淀相组织形貌影响规律;
(2)分析形貌演化图,得到沉淀相形核孕育期、形核、长大、粗化规律;沉淀相稳定性、弥散分布或取向排列的条件;同相、异相间界面关系、稳定性规律;
(3)绘制沉淀颗粒的中心至相界的原子占位变化曲线,分析得到沉淀相沉淀机制;
(4)针对具体沉淀相,可绘制组分原子在亚晶格格点位置的时间演化曲线,得到该沉淀相亚晶格位置组分变化规律,反位缺陷演化规律,原子扩散通量和路径;
(5)解析沉淀相异相、同相间界面关系及其稳定性,分析界面迁移机制,得到沉淀相长大时原子迁移规律;绘制界面两侧溶质、溶剂的原子分布和取向扩散,得到界面迁移的微观机制。
上述求解方程的方法使用的数值方法是欧拉迭代法。
上述沉淀相形貌与温度、应力、组分均有关系;其中宏/微耦合应力作用下,沉淀相取向排列,倾向于形成针状、棒状形貌。
上述宏观耦合应力指时效成形的外加应力;所述微观耦合应力指异相沉淀相间的晶格错配应力,其中晶格错配应力与沉淀相错配度、相取向、弹性常数及分布均有关系。
上述原子占位定量表征形貌演化和晶界迁移是通过追踪在某个格点位置所有组元的占位变化,得到该位置原子占位的动态变化规律;进而把选定区域放大,可定量得出沉淀相、相界面、界面两侧包括空间和时间的原子扩散规律,通过原子的扩散路径,分析组元簇聚、贫化、富化特点,解析沉淀相取向排列形成规律、界面结构稳定性规律。
该发明的有益效果在于:
该方法通过追踪应力时效全过程组织从无序固溶态至沉淀相析出、长大、粗化过程的微观形貌演化,解析沉淀相析出机制、取向生长和粗化规律、晶间无沉淀带形成,探明两相弹性模量差、沉淀相结构、相界结构等微观相界失配应力场和蠕变成形宏观应力场诱发原子穿越结合面扩散的促进作用,得出结合层上沉淀行为变化规律,达到消除结合面连续沉淀相的目的。
附图说明
图1是本发明的技术路线图;
图2是无应力时效(左)/应力时效(右)的形貌对照图;
图3是应力时效温度不同时的形貌图;
图4是沉淀相中的原子占位时间演化曲线;
图5是有序沉淀相长程序参数随时间的演化;
图6是应变能对组分原子占位的影响;
图7是界面形态和稳定性。
具体实施方式
下面结合实施例对本发明做进一步的描述。
图1是本发明的技术路线图。按照该技术路线图设定点的思路和步骤,可达到消除界面层不连续性的目的。首先,微观相场动力学方程的无需预先设定沉淀相,从秒级至数百小时的时间尺度,可得出GP区、过渡相、稳定相的形核、长大、粗化的全过程,具有物理原理优势。从软件条件角度,开发的宏观弹性应力场、微观相界失配应力场表示为自由能方程式的场变量函数,建立的宏观弹性应力场、纳观相界失配应力场函数与微扩散公式耦合的动力学方程,适合于含纳米尺度颗粒的各向异性弹塑性非均质体系,建立了模拟动态沉淀的条件。计算结果与实验对照,具有较好的一致性。
图2是无应力时效(左)/应力时效(右)的形貌对照图。图中有两种面心立方衍生结构的沉淀相A和B。无应力时效时,沉淀相弥散排列,形状为不规则圆形或椭圆形。考虑时效成形宏观应力和微观弹性错配应力后,沉淀相呈明显的取向排列。该现象说明,应力诱发原子定向扩散,致使沉淀相取向粗化,呈取向排列形貌。
时效成形/扩散复合的温度,对沉淀相形貌影响很大,变换温度的情况下,其组织形貌图如图3所示。图3选了六个温度,(a)-(f)的温度越来越高,从(b)图开始,取向粗化排列的沉淀相之间,包括同相之间、异相之间,有无沉淀相带出现,其宽度在一定的温度范围变化不大,在较高温度下,逐渐增大,如(f)图所示,沉淀颗粒粗化严重,带状取向分布明显,无沉淀带较宽。
图4是Ni3Al有序沉淀相中原子占位几率演化曲线。曲线可分为三个阶段,开始时期的平直对应的孕育期,原子上升或下降阶段对应的形核和长大,以及随后的又趋于平直对应粗化。在该相中Ni原子占据面心位置,Al原子占位顶角位置。孕育期原子占位几率几乎没有变化;随时效时间延长,在面心位置,Ni的占位开始增加而该位置的Al则下降,在顶角位置,Al原子的占位开始增加而该位置的Ni占位则呈下降趋势;继续保温,原子占位趋于平衡,但是Al在面心、Ni在顶角都有一定的占据几率。
图5是该相的长程序参数演化曲线,用于表征沉淀相有序度,解析沉淀相沉淀机制。观察该图,随时效时间延长,有序度增加,当有序度增加到趋于1时,逐渐平衡。序参数呈大范围小起伏变化,可以判定该相的沉淀机制为失稳分解。
应力时效的原子占位减去无应力时效的原子占位,得到的一组数据绘制图6,用于表征应力对原子占位的影响。(a)图表示面心立方面心位置,(b)图表示面心立方顶角位置。可以看出在两个位置,三种合金元素占位均受应变能影响,尤其是在顶角位置。
图7是一个平衡的64×64矩阵格点微结构形貌图。该图列举了三种类型的相与相的界面类型:第一,B结构(100)面和A结构的异相界面(如 a,b,c),该例中,B和A的相界就是B和A结构的共用面——(100)面,原子有序排列;第二,B结构(001)面和A结构的异相界面(如 d,e),该例中,B和A的相界不是某个相完整的晶面,而是2~4个原子厚度的弥散过渡界面,界面处原子无序排列;第三,B结构同相界面,即B结构的(100)面和另一B结构的(001)面的界面(如 g, h), g,h是沿Y轴排列变体和沿Z轴排列变体的B结构同相相界,和异相变体一样,这种界面有2~4个原子厚度的弥散过渡界面,界面处原子无序排列。其中f的B结构(100)面和(001)面均与A结构相邻。
以上所述是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也视为本发明的保护范围。

Claims (6)

1.一种界面层互扩散优化时效成形/扩散复合组织的相场模拟方法,其特征在于,该方法是:
1)用原子尺度的相场扩散方程表征时效成形/扩散复合全过程界面层沉淀相的形核、长大、粗化过程的形貌演化;该形貌演化通过把原子占位的数据信息转化为微结构演化的图形信息来实现;
2)分析上述图形信息,观察界面层沉淀组织形貌演化,用原子占位定量表征形貌演化过程,明晰不同纳观相界失配应力场和宏观弹性应力场作用下原子取向扩散规律,得出原子取向扩散导致的沉淀相取向粗化和沿晶无沉淀带的界面层不连续性,明晰应力场对界面层不连续性的影响规律;
3)分析上述形貌演化过程中同相沉淀相、异相沉淀相间界面结构及其演化,用原子占位定量表征不同界面结构边界处、界面两侧的演化,分析原子簇聚和贫化机理,得到沉淀相相界稳定性和迁移规律,明晰沉淀相长大、取向粗化以及沿晶无沉淀带形成机制;
4)结合2)形貌演化和3)界面演化,加以原子取向扩散定量分析,得到影响界面层不连续性的因素,得出消除结合面不连续沉淀相的规律。
2.根据权利要求1所述界面层互扩散优化时效成形/扩散复合组织的相场模拟方法,其特征在于,该方法具体包括:
分为两个步骤进行,第一,数据信息转化为图形信息;第二,分析图形信息得出有益于消除界面层不连续性的结果;
第一 数据信息转化为图形信息:
(1)依据微观相场理论编制方程,并求解;
(2)设置初始的变量值,环境变量:温度、成分、外加应力;固有参量:晶格常数,弹性常数,原子间相互作用势,热起伏;计算参量:格点数,迭代步数,迭代步长;
(3)方程求解过程是在倒易空间下进行的,并间隔一定步数转为实空间做溢出判定,判定条件为:占位几率在0~1之间,无溢出,程序继续执行;占位几率在0~1区间之外,计算终止,返回修改参数;
(4)计算结束得到一组原子占位几率数值,并把该组数值转化为图形信息;
第二 分析图形信息得出有益于消除界面层不连续性的结果:
(1)绘制不同温度、成分、应力条件下的形貌演化图,得到环境变量对沉淀相组织形貌影响规律;
(2)分析形貌演化图,得到沉淀相形核孕育期、形核、长大、粗化规律;沉淀相稳定性、弥散分布或取向排列的条件;同相、异相间界面关系、稳定性规律;
(3)绘制沉淀颗粒的中心至相界的原子占位变化曲线,分析得到沉淀相沉淀机制;
(4)针对具体沉淀相,绘制组分原子在亚晶格格点位置的时间演化曲线,得到该沉淀相亚晶格位置组分变化规律,反位缺陷演化规律,原子扩散通量和路径;
(5)解析沉淀相异相、同相间界面关系及其稳定性,分析界面迁移机制,得到沉淀相长大时原子迁移规律;绘制界面两侧溶质、溶剂的原子分布和取向扩散,得到界面迁移的微观机制。
3.根据权利要求2所述界面层互扩散优化时效成形/扩散复合组织的相场模拟方法,其特征在于:所述求解方程的方法使用的数值方法是欧拉迭代法。
4.根据权利要求1~3任一所述界面层互扩散优化时效成形/扩散复合组织的相场模拟方法,其特征在于:所述沉淀相形貌与温度、应力、组分均有关系;其中宏观/微观耦合应力作用下,沉淀相呈取向排列,倾向于形成针状、棒状形貌。
5.根据权利要求4所述界面层互扩散优化时效成形/扩散复合组织的相场模拟方法,其特征在于:所述宏观耦合应力指时效成形的外加应力;所述微观耦合应力指异相沉淀相间的晶格错配应力,其中晶格错配应力与沉淀相错配度、相取向、弹性常数及分布均有关系。
6.根据权利要求5所述界面层互扩散优化时效成形/扩散复合组织的相场模拟方法,其特征在于:所述原子占位定量表征形貌演化和晶界迁移是通过追踪在某个格点位置所有组元的占位变化,得到该位置原子占位的动态变化规律;进而把选定区域放大,定量得出沉淀相、相界面、界面两侧包括空间和时间的原子扩散规律,分析组元簇聚、贫化、富化特点,解析沉淀相取向排列形成规律、界面结构稳定性规律。
CN201410744317.6A 2014-12-09 2014-12-09 界面层互扩散优化时效成形/扩散复合组织的相场模拟方法 Expired - Fee Related CN104461691B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410744317.6A CN104461691B (zh) 2014-12-09 2014-12-09 界面层互扩散优化时效成形/扩散复合组织的相场模拟方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410744317.6A CN104461691B (zh) 2014-12-09 2014-12-09 界面层互扩散优化时效成形/扩散复合组织的相场模拟方法

Publications (2)

Publication Number Publication Date
CN104461691A CN104461691A (zh) 2015-03-25
CN104461691B true CN104461691B (zh) 2017-11-21

Family

ID=52907795

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410744317.6A Expired - Fee Related CN104461691B (zh) 2014-12-09 2014-12-09 界面层互扩散优化时效成形/扩散复合组织的相场模拟方法

Country Status (1)

Country Link
CN (1) CN104461691B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109522675B (zh) * 2018-12-10 2022-08-16 桂林电子科技大学 锡基二元共晶合金微观组织的模拟及有限元求解分析方法
CN112084622B (zh) * 2020-07-28 2023-06-23 浙江师范大学 复合材料微观形貌的模拟方法、装置和电子装置

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102930115A (zh) * 2012-11-16 2013-02-13 中国航空工业集团公司北京航空制造工程研究所 基于有限元模具型面回弹补偿的壁板蠕变时效成形方法
CN103678890A (zh) * 2013-11-28 2014-03-26 西北工业大学 一种采用晶体相场方法模拟加热工艺对晶界预熔及熔化影响的方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005038547A (ja) * 2003-07-17 2005-02-10 National Institute For Materials Science 磁気記録材料の組織形成設計方法および組織形成設計装置

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102930115A (zh) * 2012-11-16 2013-02-13 中国航空工业集团公司北京航空制造工程研究所 基于有限元模具型面回弹补偿的壁板蠕变时效成形方法
CN103678890A (zh) * 2013-11-28 2014-03-26 西北工业大学 一种采用晶体相场方法模拟加热工艺对晶界预熔及熔化影响的方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Ni64Al21V15沉淀过程DO22向L10(M=1)结构转变的微观相场模拟;常秀丽等;《稀有金属材料与工程》;20130630;第42卷(第6期);第1113-1117页 *
Ni75Cr25-xAlx合金中L12相和D022相形核孕育期的微观相场模拟;卢艳丽等;《物理学报》;20090630;第58卷;第S319-S326页 *
微观相场法模拟Ni75Al5.3V19.7中L12和D022结构反位缺陷的演化;张静等;《物理学报》;20090131;第58卷(第1期);第631-637页 *

Also Published As

Publication number Publication date
CN104461691A (zh) 2015-03-25

Similar Documents

Publication Publication Date Title
Aslan et al. Micromorphic approach to single crystal plasticity and damage
Guo et al. Micromechanical simulation of fracture behavior of bimodal nanostructured metals
Liu et al. The effect of void defect on the evolution mechanisms of dislocations and mechanical properties in nickel-based superalloys by molecular dynamics simulation of real γ/γ′ structures
CN103770935B (zh) 乘波体外形设计方法
Liang et al. Material point method with enriched shape function for crack problems
CN104461691B (zh) 界面层互扩散优化时效成形/扩散复合组织的相场模拟方法
Zhang et al. A molecular dynamics study of intercrystalline crack propagation in nano-nickel bicrystal films with (0 1 0) twist boundary
CN105760631A (zh) 一种实现多样性设计的结构拓扑优化方法
Jiang et al. Molecular dynamics study of microscopic deformation mechanism and tensile properties in AlxCoCrFeNi amorphous high-entropy alloys
Feng et al. Twin and dislocation mechanisms in tensile W single crystal with temperature change: a molecular dynamics study
CN109117604A (zh) 一种镍基单晶合金基于滑移系的Walker模型材料常数的获取方法
Tursun et al. The influence of transition phases on the damage behaviour of an Al/10vol.% SiC composite
Qiang et al. Synergistic strengthening of dislocation and heterodeformation induced in gradient nanograined copper film: A molecular dynamics study
CN108153928A (zh) 一种含夹杂物粉末高温合金裂纹萌生寿命预测方法
Wang et al. Three-dimensional finite element analysis using crystal plasticity for a parameter study of microstructurally small fatigue crack growth in a AA7075 aluminum alloy
CN105447232B (zh) 用于橡胶金属复合止挡大变形橡胶分析的网格生成方法
Weber et al. Thermo-mechanical deformation evolution in polycrystalline Ni-based superalloys by a hierarchical crystal plasticity model
CN104537253A (zh) 一种时效成形预时效过程的微观相场分析方法
Lin et al. Investigation of the effects of point defects on the tensile strength of BCC-Fe using molecular dynamics
Han et al. Research on the effect of micro-voids on the deformation behavior and crack initiation lifetime of titanium alloy under cyclic loading by crystal plasticity finite element method
Hue et al. Atomic-level analysis of unusual mechanical and failure behaviors in compositionally graded nanowires: A molecular dynamics study
CN107219136A (zh) 一种考虑不同减薄率的超塑性成形件疲劳寿命预测方法
Morgan et al. Pressure-driven phase transitions in crystalline nanoparticles: Surface effects on hysteresis
Wang et al. Molecular dynamics investigation of the grain boundary migration hysteresis of nanocrystalline Ni under cyclic shear loading
Liu et al. The fatigue crack growth in hierarchically nano-twinned materials

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20171121

Termination date: 20181209

CF01 Termination of patent right due to non-payment of annual fee