CN104408300A - 时效成形/扩散复合工艺的微观相场建模与分析方法 - Google Patents
时效成形/扩散复合工艺的微观相场建模与分析方法 Download PDFInfo
- Publication number
- CN104408300A CN104408300A CN201410658159.2A CN201410658159A CN104408300A CN 104408300 A CN104408300 A CN 104408300A CN 201410658159 A CN201410658159 A CN 201410658159A CN 104408300 A CN104408300 A CN 104408300A
- Authority
- CN
- China
- Prior art keywords
- field
- diffusion
- atom
- phase
- microcosmic
- 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.)
- Pending
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种时效成形/扩散复合工艺的微观相场建模与分析方法,用于解决现有时效成形/扩散复合工艺导致材料机械性能差的技术问题。技术方案是构建时效成形/扩散复合工艺的微观相场模型,以耦合温度场、浓度场和宏/微观弹性应力场的微观原子扩散方程求解原子占位信息。再将原子占位信息转化成组织形貌演化等二次信息,用于分析温度场、浓度场和弹性应力场对时效成形/扩散复合过程的组织形貌演变和界面结构迁移的影响规律。组织形貌演化表达了时效成形/扩散复合过程的组织变化,原子占位能定量表征组织变化的微观原子扩散机理,温度场、浓度场和弹性应力场等的组织响应度能反推得到消除取向粗化的工艺参数。提升了材料的机械性能。
Description
技术领域
本发明涉及一种微观相场建模与分析方法,特别涉及一种时效成形/扩散复合工艺的微观相场建模与分析方法。
背景技术
时效成形是蠕变成形与动态析出相结合的精确成形技术,在人工时效温度以蠕变变形实现铝合金构件的精确成形,克服了冷成形极限问题,获得动态析出的组织和性能。时效成形大型复杂整体构件取代焊接、铆接、胶接的装配复杂壁板,减少开孔、紧固件等机械连接,既下降应力集中,又明显减轻重量,显著提高服役寿命与可靠性,呈现出高的结构效率,是现代飞机大型复杂整体构件制造工艺技术的一个重要标志。
基于时效成形原理以及铝合金层状复合材料研究,发展了由构件时效成形工艺扩展到复合材料制备和构件成形的一体化技术——扩散复合/时效成形方法,即两种合金分别在相应的温度固溶处理之后,在时效成形设备进行时效成形/扩散复合,蠕变作用下接触表面原子由物理作用距离进入化学作用距离,物理接触面转化为化学结合面;扩散复合阶段,界面层两侧元素互扩散,削弱形成结合面连续沉淀相的倾向;兼实现精确成形,制备出一种合金作为表面层和一种合金作为中心层的复合材料制件,同时实现形变时效的组织性能。
时效成形过程蠕变变形和沉淀强化同时进行[International Journal of MachineTools&Manufacture 51(2011)1–17]。基于统一理论和时效动力学的蠕变本构方程常用于刻画应力时效硬化、蠕变诱发沉淀相演化,位错强化,固溶强化和时效强化、多级时效沉淀相演化[Ho K.C,Lin J,Dean T.A.J.Mater.Process.Tech.2004153-154:122-127]。本构方程在刻画沉淀相演化方面为时效成形提供了有益理论支持,然而,时效成形/扩散复合过程溶质原子滞留在结合层,沉淀相取向粗化严重,沿晶无沉淀带分布其中,组织不均匀性严重损害材料性能,本构方程在空间尺度上无法观察原子级别组织不均匀性,在时间尺度上无法捕捉沉淀相形核、长大和取向粗化的演化过程。
发明内容
为了克服现有时效成形/扩散复合工艺导致材料机械性能差的不足,本发明提供一种时效成形/扩散复合工艺的微观相场建模与分析方法。该方法通过构建时效成形/扩散复合工艺的微观相场模型,以耦合温度场、浓度场和宏/微观弹性应力场的微观原子扩散方程求解原子占位信息。再将原子占位信息转化成组织形貌演化等二次信息,用于分析温度场、浓度场和弹性应力场对时效成形/扩散复合过程的组织形貌演变和界面结构迁移的影响规律。组织形貌演化表达了时效成形/扩散复合过程的组织变化,原子占位能定量表征组织变化的微观原子扩散机理,温度场、浓度场和弹性应力场等的组织响应度能反推得到消除取向粗化的工艺参数。该方法能减少实验试错性,优化时效成形/扩散复合组织,提升了材料的机械性能。
本发明解决其技术问题所采用的技术方案:一种时效成形/扩散复合工艺的微观相场建模与分析方法,其特点是包括以下步骤:
步骤一、构建二维微观相场模型的Langevin形式微观扩散方程。
采用Onsager型离散格点形式的动力学方程描述合金沉淀,采用平均场理论计算自由能函数。
(a)据Onsager扩散方程,几率的变化率与热力学驱动力成正比,即:
表示A原子在t时刻占据晶格位置r的几率,C0为基体平均浓度,L(r-r′)为与单位时间内由格点r′跃迁至的几率有关的常数;T为绝对温度;KB为玻尔兹曼常数;C0为基体平均浓度;F为***的总自由能,平均场近似对离散格点模型中***总自由能的表示式为:
V(r-r′)为原子间有效作用能,由下式给出,
V(r-r′)=VAA(r-r′)+VBB(r-r′)-2VAB(r-r′) (3)
VAA是A原子之间的相互作用势,VAB是A原子和B原子之间的相互作用势。A、B为合金组元。
(b)为简化计算,将面心立方的三维空间在[001]取向上投影,方程(1)在倒易空间的平面投影中,动力学方程为:
ξ(k,t)为随机噪声项,满足涨落-耗散理论的随机项。
其中,
V1,V2分别为第一近邻、次近邻原子间有效交互作用能。
步骤二、时效成形/扩散复合的宏/微观耦合弹性应力。
微观相场模型中描述沉淀相的晶格错配能的微观弹性表达式为,
其中ε0=(ap-a0)/a0为原子尺寸差异引起的原子尺寸失配程度的晶格点阵错配度,cij为立方晶格的弹性常数。
外应力引起的弹性应变能在傅立叶空间下有,
A为常数,a0为外力为0时原子的平均间距,a为原子间距,且a=a0+Δa,Δa为外应力引起的形变量,E<100>、E<001>是[100]和[001]上的弹性模量,且△a=σa0/E。σx为外力在x轴的分量,σy为外力在y轴的分量。
将错配及外应力引起的弹性应变能(8)式引入自由能表达式(2)式中,计算时效成形和扩散复合的宏观应力以及结构演化引入的微观弹性应力的耦合应力作用。
步骤三、方程求解和数据处理。
用欧拉迭代法求解微观扩散方程,计算得到原子占位信息。以原子间作用势,热起伏,投影后的晶格信息,弹性常数等作为输入参量,在倒易空间下求解相场方程得到原子占位几率值,然后变换到实空间用于绘制原子演化图。在计算机求解扩散方程过程,如果占位几率大于1或者小于0,程序终止;如果占位几率是0和1之间的值,程序继续。计算完毕,处理数据过程把实空间的原子占位几率以四维矩阵形式保存,并用于绘制形貌演化图和后续的原子占位分析。通过分析形貌演化图,得到固溶态有序化过程,沉淀相的形核、长大及粗化过程以及时间、成分、耦合应力对微观沉淀组织形貌的影响;进一步分析形貌演化图,得到同相、异相沉淀相的界面结构信息;用原子占位定量表征沉淀相演化的原子扩散规律,界面处原子组成、分布及取向扩散规律。
本发明的有益效果是:该方法通过构建时效成形/扩散复合工艺的微观相场模型,以耦合温度场、浓度场和宏/微观弹性应力场的微观原子扩散方程求解原子占位信息。再将原子占位信息转化成组织形貌演化等二次信息,用于分析温度场、浓度场和弹性应力场对时效成形/扩散复合过程的组织形貌演变和界面结构迁移的影响规律。组织形貌演化表达了时效成形/扩散复合过程的组织变化,原子占位能定量表征组织变化的微观原子扩散机理,温度场、浓度场和弹性应力场等的组织响应度能反推得到消除取向粗化的工艺参数。该方法能减少实验试错性,优化时效成形/扩散复合组织,提升了材料的机械性能。
下面结合附图和具体实施方式对本发明作详细说明。
附图说明
图1是本发明方法的流程图。
图2是本发明实施例中时效成形/扩散复合过程组织演化图。
图3是本发明实施例中时效成形/扩散复合温度对微观组织形貌的影响图。
图4是本发明实施例中时效成形/扩散复合的宏/微观耦合应力作用下不同成分组织的取向生长图。
图5是本发明实施例中面心立方结构的两种不同沉淀相界面结构图。
图6是本发明实施例中合金元素在异相畴界及其两侧的分布及其演化图。
具体实施方式
以下参照图1-6详细说明本发明。
步骤一、二维微观相场模型的Langevin形式微观扩散方程构建。
采用Onsager型离散格点形式的动力学方程描述合金沉淀,采用平均场理论计算自由能函数。可定量描述原子尺度微结构演化。
(a)据Onsager扩散方程可知,几率的变化率与热力学驱动力成正比,即:
表示A原子在t时刻占据晶格位置r的几率,C0为基体平均浓度,L(r-r′)为与单位时间内由格点r′跃迁至的几率有关的常数;T为绝对温度;KB为玻尔兹曼常数;C0为基体平均浓度;F为***的总自由能,平均场近似对离散格点模型中***总自由能的表示式为:
V(r-r′)为原子间有效作用能,由下式给出,
V(r-r′)=VAA(r-r′)+VBB(r-r′)-2VAB(r-r′) (3)
VAA是A原子之间的相互作用势,VAB是A原子和B原子之间的相互作用势。
A、B为合金组元。
(b)为简化计算,将面心立方的三维空间在[001]取向上投影,方程(1)在倒易空间的平面投影中,动力学方程为:
ξ(k,t)为随机噪声项,满足涨落-耗散理论的随机项。
其中,
V1,V2分别为第一近邻、次近邻原子间有效交互作用能。
步骤二、时效成形/扩散复合的宏/微观耦合弹性应力。
微观相场模型中描述沉淀相的晶格错配能的微观弹性表达式为,
其中ε0=(ap-a0)/a0为原子尺寸差异引起的原子尺寸失配程度的晶格点阵错配度,cij为立方晶格的弹性常数。
外应力引起的弹性应变能在傅立叶空间下有,
A为常数,a0为外力为0时原子的平均间距,a为原子间距,且a=a0+Δa,Δa为外应力引起的形变量,E<100>、E<001>是[100]和[001]上的弹性模量,且△a=σa0/E。σx为外力在x轴的分量,σy为外力在y轴的分量。
将错配及外应力引起的弹性应变能(8)式引入自由能表达式(2)式中,计算时效成形和扩散复合的宏观应力以及结构演化引入的微观弹性应力的耦合应力作用。
步骤三、方程求解和数据处理。
用欧拉迭代法求解微观扩散方程,计算得到原子占位信息。以原子间作用势,热起伏,投影后的晶格信息,弹性常数等作为输入参量,在倒易空间下求解相场方程得到原子占位几率值,然后变换到实空间用于绘制原子演化图。在计算机求解扩散方程过程,如果占位几率大于1或者小于0,程序终止;如果占位几率是0和1之间的值,程序继续。计算完毕,处理数据过程把实空间的原子占位几率以四维矩阵形式保存,并用于绘制形貌演化图和后续的原子占位分析。通过分析形貌演化图,可得到固溶态有序化过程,沉淀相的形核、长大及粗化过程以及时间、成分、耦合应力对微观沉淀组织形貌的影响;进一步分析形貌演化图,可得到同相、异相沉淀相的界面结构信息;用原子占位定量表征沉淀相演化的原子扩散规律,界面处原子组成、分布及取向扩散规律。
(1)微观相场方程构建:建立耦合了宏观变形应力和微观弹性错配应力的微观相场方程,在倒易空间下求解该方程,方程的解就是倒易空间下原子占位信息,对该倒易空间下的原子占位信息进行傅里叶变换,得到实空间下的原子占位几率值。该占位几率用于绘制形貌演化图,原子占位演化图曲线,原子扩散,沉淀相组分分析的原始数据。
(2)溢出判定:把倒易空间下的原子占位信息每间隔一定的时间步数(间隔可以没有,也可以选择任意整数计算时间步,一般选择1000步),转化为实空间下的占位几率值进行一次溢出判定,判定依据是原子在某个格点的占位不可能大于1或者小于0,因此当计算结果在0~1之间时,计算符合现实情况,可计算可继续进行,而在0~1区间之外时,程序终止,需修改参数,再进行新一轮的计算。
(3)转化到实空间的占位几率数值存储在一个(n,n,t,x)的四维矩阵中,n,t,x分表表计算的空间尺度、时间尺度,以及组分数目。其中t与时间正相关,但无具体对应关系,t表征从计算开始到终了的状态数目,可根据需求选择存储数目的多少。例如:计算总步数为10 0000步,如果选择每1个计算步存储数据一次数据,则t=10 0000,优点是可观察每个计算步的微小变化,缺点是存储量过大对计算机存储量是个巨大挑战;如果选择每1000步存储一次数据,则t=10 0000/1000,即t=100,表示本次计算每1000步存储一次数据,共存储100组数据。取(n,n,t,x)四维矩阵中的某个组元的占位几率,绘制形貌演化图,即得到一个n×n大小,共t副的一组形貌演化图。分析组织形貌演化信息得到取向排列、定向粗化、界面结构及组分扩散等信息。
具体结果如下:
a.组织形貌。
参照图2、图3、图4的时间、温度、成分相关的形貌演化图,该演化图的空间尺度是128×128格点的二维点阵(也可选择256×256)。图2是时效成形/扩散复合过程组织随时间的演化图,图2(a)是时效成形/扩散复合初期的形貌图,显示了初期大面积的浓度起伏和结构起伏,虽然结构还不能清晰可辨。时效成形过程,热作用和耦合应力作用下,原子运动剧烈,扩散速度快,观察到图2(b)中有沉淀相析出,用A表示,A沉淀相的面心立方结构构型清晰可辩。随时效成形/扩散复合过程时间的延长,A沉淀相长大粗化明显。图2(d)中开始有红色沉淀相在A沉淀相相的边界处析出,用B表示,在边界处吞噬A沉淀相而逐渐长大、粗化,体系中有A和B两种面心立方沉淀相。
图3是温度变化时,时效成形/扩散复合温度对微观组织形貌的影响,图3(a)-(f)为温度逐渐增加时的形貌图。低温下,A和B两种结构沉淀相充满整个体系,且两沉淀相颗粒细密,相互之间呈独立、弥散分布。随温度升高,沉淀相数目越来越少,颗粒半径越来越大,沉淀相逐渐由低温下的单个颗粒向高温下的互联状过渡。温度越低,过冷度越大,形核数量越多,形成致密的两相区;温度升高,晶核的过冷度减小,但长程扩散占据主导地位,形核所需克服的势垒随温度升高而增大,这时的晶核难形成。温度更高时有少量无序相存在,温度越高,无序相占的比重越大。
图4是宏/微观耦合应力作用下不同成分组织的取向生长形貌图。随着合金成分调整,A和B沉淀相的比例发生变化,随着沉淀相比例的变化,沉淀相取向排列也随着变化。当B沉淀相占主导时,A沉淀相单向取向分布,且取向明显,B沉淀相呈筏状在取向分布的A沉淀相间排列,并且粗化严重。随A沉淀相的比例增加,取向分布依然存在,但是取向的方向性不明显,取向成长方向由单向生长转向[001]和[100]两个方向生长。
b.界面和占位。
图5是计算得到的A和B沉淀相异相界面结构,图5中包含五种异相间界面,即(100)B//(200)A(箭头A、B所示)、(100)B//(200)A·1/2[001]A,(002)B//(002)A·1/2[100]A和(002)B//(002)A(箭头C、D和E所示),(002)B//(001)A(箭头F所示)。时效成形过程,A和B沉淀相各自粗化,两相比例发生变化,这一过程的异相界面迁移规律是除(002)B//(001)A之外,其余四种有序畴界均可发生迁移。(100)B//(200)A和(100)B//(200)A·1/2[001]A在迁移前后有序畴界结构保持不变,而(002)B//(002)A和(100)B//(200)A·1/2[001]A则在迁移过程中交替出现。
图6是合金元素在有序畴界及其两侧的分布及其演化图。(a)、(b)、(c)、(d)分别为(100)B//(200)A,(100)B//(200)A·1/2[001]A,(002)B//(002)A·1/2[100]A和(002)B//(002)A,(002)B//(001)A处(对应图5中箭头A、C、E、G所指有序畴界)及其两侧的合金元素分布。图中平行于纵坐标的实线表示有序畴界的初始位置,虚线表示有序畴界迁移后的终止位置。随着时间的进行,虽然有序畴界发生迁移,但合金元素在有序畴界处的偏聚和贫化倾向不发生改变。
Claims (1)
1.一种时效成形/扩散复合工艺的微观相场建模与分析方法,其特征在于包括以下步骤:
步骤一、构建二维微观相场模型的Langevin形式微观扩散方程;
采用Onsager型离散格点形式的动力学方程描述合金沉淀,采用平均场理论计算自由能函数;
(a)据Onsager扩散方程,几率的变化率与热力学驱动力成正比,即:
表示A原子在t时刻占据晶格位置r的几率,C0为基体平均浓度,L(r-r′)为与单位时间内由格点r′跃迁至的几率有关的常数;T为绝对温度;KB为玻尔兹曼常数;C0为基体平均浓度;F为***的总自由能,平均场近似对离散格点模型中***总自由能的表示式为:
V(r-r′)为原子间有效作用能,由下式给出,
V(r-r′)=VAA(r-r′)+VBB(r-r′)-2VAB(r-r′) (3)
VAA是A原子之间的相互作用势,VAB是A原子和B原子之间的相互作用势;A、B为合金组元;
(b)为简化计算,将面心立方的三维空间在[001]取向上投影,方程(1)在倒易空间的平面投影中,动力学方程为:
ζ(k,t)为随机噪声项,满足涨落-耗散理论的随机项;
其中,
V1,V2分别为第一近邻、次近邻原子间有效交互作用能;
步骤二、时效成形/扩散复合的宏/微观耦合弹性应力;
微观相场模型中描述沉淀相的晶格错配能的微观弹性表达式为,
其中ε0=(ap-a0)/a0为原子尺寸差异引起的原子尺寸失配程度的晶格点阵错配度,cij为立方晶格的弹性常数;
外应力引起的弹性应变能在傅立叶空间下有,
A为常数,a0为外力为0时原子的平均间距,a为原子间距,且a=a0+Δa,Δa为外应力引起的形变量,E<100>、E<001>是[100]和[001]上的弹性模量,且Δa=σa0/E;σx为外力在x轴的分量,σy为外力在y轴的分量;
将错配及外应力引起的弹性应变能(8)式引入自由能表达式(2)式中,计算时效成形和扩散复合的宏观应力以及结构演化引入的微观弹性应力的耦合应力作用;
步骤三、方程求解和数据处理;
用欧拉迭代法求解微观扩散方程,计算得到原子占位信息;以原子间作用势,热起伏,投影后的晶格信息,弹性常数等作为输入参量,在倒易空间下求解相场方程得到原子占位几率值,然后变换到实空间用于绘制原子演化图;在计算机求解扩散方程过程,如果占位几率大于1或者小于0,程序终止;如果占位几率是0和1之间的值,程序继续;计算完毕,处理数据过程把实空间的原子占位几率以四维矩阵形式保存,并用于绘制形貌演化图和后续的原子占位分析;通过分析形貌演化图,得到固溶态有序化过程,沉淀相的形核、长大及粗化过程以及时间、成分、耦合应力对微观沉淀组织形貌的影响;进一步分析形貌演化图,得到同相、异相沉淀相的界面结构信息;用原子占位定量表征沉淀相演化的原子扩散规律,界面处原子组成、分布及取向扩散规律。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410658159.2A CN104408300A (zh) | 2014-11-18 | 2014-11-18 | 时效成形/扩散复合工艺的微观相场建模与分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410658159.2A CN104408300A (zh) | 2014-11-18 | 2014-11-18 | 时效成形/扩散复合工艺的微观相场建模与分析方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104408300A true CN104408300A (zh) | 2015-03-11 |
Family
ID=52645931
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410658159.2A Pending CN104408300A (zh) | 2014-11-18 | 2014-11-18 | 时效成形/扩散复合工艺的微观相场建模与分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104408300A (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105956249A (zh) * | 2016-04-25 | 2016-09-21 | 杭州电子科技大学 | 高温下基于相场模型的硅基微结构形变机理研究方法 |
CN106780407A (zh) * | 2017-03-01 | 2017-05-31 | 成都优途科技有限公司 | 一种针对超声图像斑点噪声的去噪***及去噪方法 |
CN109740205A (zh) * | 2018-12-18 | 2019-05-10 | 江苏理工学院 | 基于耦合晶界扩散相场法的多元合金固态相变的模拟方法 |
CN109991387A (zh) * | 2019-03-06 | 2019-07-09 | 华南理工大学 | 一种模拟γ-TiAl合金在非等温条件下的相场分析法 |
CN110008531A (zh) * | 2019-03-15 | 2019-07-12 | 北京应用物理与计算数学研究所 | 一种含间隙原子的无序固溶材料原子结构建模方法 |
CN110033832A (zh) * | 2019-04-10 | 2019-07-19 | 南京理工大学 | 一种应力调控纳米尺度颗粒微观组织的方法 |
CN110717271A (zh) * | 2019-10-09 | 2020-01-21 | 中国原子能科学研究院 | 一种基于指数时间差分格式求解的物质演化模拟方法 |
CN112071371A (zh) * | 2020-08-28 | 2020-12-11 | 武汉大学 | 一种用于物质和材料的计算机仿真模拟方法及装置 |
CN112084622A (zh) * | 2020-07-28 | 2020-12-15 | 浙江师范大学 | 复合材料微观形貌的模拟方法、装置和电子装置 |
CN112464532A (zh) * | 2020-11-27 | 2021-03-09 | 湖南工业大学 | 一种纳米级硬质合金颗粒生长规律研究方法 |
WO2022042522A1 (zh) * | 2020-08-28 | 2022-03-03 | 武汉大学 | 用于物质和材料的计算机仿真模拟方法及装置 |
WO2022247870A1 (zh) * | 2021-05-27 | 2022-12-01 | 复旦大学 | 选通器开关的仿真方法及*** |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102925826A (zh) * | 2012-11-16 | 2013-02-13 | 中国航空工业集团公司北京航空制造工程研究所 | 一种基于热压罐的整体壁板蠕变时效成形工装 |
CN103464536A (zh) * | 2013-08-12 | 2013-12-25 | 西北工业大学 | 一种弹性变形条件下的振动时效成形方法及装置 |
KR20140041121A (ko) * | 2012-09-27 | 2014-04-04 | (주)대우건설 | 대규모 구조물의 성능 기반 구조 해석 방법 |
-
2014
- 2014-11-18 CN CN201410658159.2A patent/CN104408300A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20140041121A (ko) * | 2012-09-27 | 2014-04-04 | (주)대우건설 | 대규모 구조물의 성능 기반 구조 해석 방법 |
CN102925826A (zh) * | 2012-11-16 | 2013-02-13 | 中国航空工业集团公司北京航空制造工程研究所 | 一种基于热压罐的整体壁板蠕变时效成形工装 |
CN103464536A (zh) * | 2013-08-12 | 2013-12-25 | 西北工业大学 | 一种弹性变形条件下的振动时效成形方法及装置 |
Non-Patent Citations (1)
Title |
---|
CHANG XIULI等: "Microscopic Phase-Field Simulation of the Structure Transformation from DO22 to L10(M=1) Phase in Ni64Al21V15 Alloy", 《RARE METAL MATERIALS AND ENGINEERING》 * |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105956249A (zh) * | 2016-04-25 | 2016-09-21 | 杭州电子科技大学 | 高温下基于相场模型的硅基微结构形变机理研究方法 |
CN105956249B (zh) * | 2016-04-25 | 2019-04-02 | 杭州电子科技大学 | 高温下基于相场模型的硅基微结构形变机理研究方法 |
CN106780407A (zh) * | 2017-03-01 | 2017-05-31 | 成都优途科技有限公司 | 一种针对超声图像斑点噪声的去噪***及去噪方法 |
CN106780407B (zh) * | 2017-03-01 | 2024-03-26 | 清远先导科臻医疗科技有限公司 | 一种针对超声图像斑点噪声的去噪***及去噪方法 |
CN109740205A (zh) * | 2018-12-18 | 2019-05-10 | 江苏理工学院 | 基于耦合晶界扩散相场法的多元合金固态相变的模拟方法 |
CN109740205B (zh) * | 2018-12-18 | 2023-07-14 | 江苏理工学院 | 基于耦合晶界扩散相场法的多元合金固态相变的模拟方法 |
CN109991387A (zh) * | 2019-03-06 | 2019-07-09 | 华南理工大学 | 一种模拟γ-TiAl合金在非等温条件下的相场分析法 |
CN110008531A (zh) * | 2019-03-15 | 2019-07-12 | 北京应用物理与计算数学研究所 | 一种含间隙原子的无序固溶材料原子结构建模方法 |
CN110033832A (zh) * | 2019-04-10 | 2019-07-19 | 南京理工大学 | 一种应力调控纳米尺度颗粒微观组织的方法 |
CN110717271A (zh) * | 2019-10-09 | 2020-01-21 | 中国原子能科学研究院 | 一种基于指数时间差分格式求解的物质演化模拟方法 |
CN110717271B (zh) * | 2019-10-09 | 2021-07-30 | 中国原子能科学研究院 | 一种基于指数时间差分格式求解的物质演化模拟方法 |
CN112084622A (zh) * | 2020-07-28 | 2020-12-15 | 浙江师范大学 | 复合材料微观形貌的模拟方法、装置和电子装置 |
CN112084622B (zh) * | 2020-07-28 | 2023-06-23 | 浙江师范大学 | 复合材料微观形貌的模拟方法、装置和电子装置 |
CN112071371A (zh) * | 2020-08-28 | 2020-12-11 | 武汉大学 | 一种用于物质和材料的计算机仿真模拟方法及装置 |
WO2022042522A1 (zh) * | 2020-08-28 | 2022-03-03 | 武汉大学 | 用于物质和材料的计算机仿真模拟方法及装置 |
CN112464532A (zh) * | 2020-11-27 | 2021-03-09 | 湖南工业大学 | 一种纳米级硬质合金颗粒生长规律研究方法 |
WO2022247870A1 (zh) * | 2021-05-27 | 2022-12-01 | 复旦大学 | 选通器开关的仿真方法及*** |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104408300A (zh) | 时效成形/扩散复合工艺的微观相场建模与分析方法 | |
Sagui et al. | Kinetics of phase separation in two-dimensional systems with competing interactions | |
Arsenlis et al. | Modeling the evolution of crystallographic dislocation density in crystal plasticity | |
Koyama et al. | Computer simulation of phase decomposition in Fe–Cu–Mn–Ni quaternary alloy based on the phase-field method | |
Mellbin et al. | Recrystallization and texture evolution during hot rolling of copper, studied by a multiscale model combining crystal plasticity and vertex models | |
Kardan-Halvaei et al. | Crystal plasticity finite element simulation and experimental investigation of the micro-upsetting process of OFHC copper | |
Kords et al. | On the role of dislocation transport in the constitutive description of crystal plasticity | |
Sun et al. | Multi-component phase-field simulation of microstructural evolution and elemental distribution in Fe–Cu–Mn–Ni–Al alloy | |
Shibuta et al. | Micrometer-scale molecular dynamics simulation of microstructure formation linked with multi-phase-field simulation in same space scale | |
CN114267417A (zh) | 基于相场法的耦合弹性能的铁铜锰合金固态相变模拟方法 | |
Weinkamer et al. | Using kinetic Monte Carlo simulations to study phase separation in alloys | |
Simeone et al. | Radiation-induced micro-structures as ground states of a Swift-Hohenberg energy functional | |
Aksimentiev et al. | Morphology of Surfaces in Mesoscopic Polymers, Surfactants, Electrons, or Reaction–Diffusion Systems: Methods, Simulations, and Measurements | |
Berinskii et al. | Application of the pair torque interaction potential to simulate the elastic behavior of SLMoS2 | |
Cai et al. | Stress relaxation aging behavior and constitutive modelling of AA7150-T7751 under different temperatures, initial stress levels and pre-strains | |
Parthasarathy et al. | Finite-temperature stress calculations in atomic models using moments of position | |
Leung et al. | Size dependence of yield strength simulated by a dislocation-density function dynamics approach | |
Schenk et al. | Dislocation densities and velocities within the γ channels of an SX superalloy during in situ high-temperature creep tests | |
Popova et al. | Anomalous diffusion of a tethered membrane: A Monte Carlo investigation | |
Yu | Multiresolution Clustering Analysis of Hierarchical Material Systems | |
He et al. | Effect of elastic strain energy on grain growth and texture in AZ31 magnesium alloy by phase-field simulation | |
Cao et al. | Analyzing effects of temperature gradient and scan rate on metal additive manufacturing microstructure by using phase field-finite element method | |
Lin et al. | Crystal plasticity-inspired statistical analysis of dislocation substructures generated by continuum dislocation dynamics | |
Haghighat et al. | Investigation of limiting grain size and microstructure homogeneity in the presence of second phase particles using the Monte Carlo method | |
Fedotov et al. | Modeling of Deformation and Fracture Process of Layered Nanocomposites |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20150311 |
|
WD01 | Invention patent application deemed withdrawn after publication |