CN102750425A - 一种焊接过程热影响区组织演变的模拟方法 - Google Patents
一种焊接过程热影响区组织演变的模拟方法 Download PDFInfo
- Publication number
- CN102750425A CN102750425A CN201210247239XA CN201210247239A CN102750425A CN 102750425 A CN102750425 A CN 102750425A CN 201210247239X A CN201210247239X A CN 201210247239XA CN 201210247239 A CN201210247239 A CN 201210247239A CN 102750425 A CN102750425 A CN 102750425A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- phase
- mfrac
- beta
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 76
- 238000003466 welding Methods 0.000 title claims abstract description 68
- 230000008569 process Effects 0.000 title claims abstract description 48
- 238000004088 simulation Methods 0.000 title claims abstract description 45
- 230000008859 change Effects 0.000 claims abstract description 69
- 238000001816 cooling Methods 0.000 claims abstract description 61
- 238000004364 calculation method Methods 0.000 claims abstract description 49
- 238000009826 distribution Methods 0.000 claims abstract description 29
- 239000012071 phase Substances 0.000 claims description 383
- 239000013078 crystal Substances 0.000 claims description 120
- 230000009466 transformation Effects 0.000 claims description 46
- 230000006911 nucleation Effects 0.000 claims description 37
- 238000010899 nucleation Methods 0.000 claims description 37
- 230000005501 phase interface Effects 0.000 claims description 31
- 229910001069 Ti alloy Inorganic materials 0.000 claims description 23
- 238000009792 diffusion process Methods 0.000 claims description 22
- 230000007704 transition Effects 0.000 claims description 22
- 230000001413 cellular effect Effects 0.000 claims description 18
- 238000004781 supercooling Methods 0.000 claims description 18
- 238000010438 heat treatment Methods 0.000 claims description 17
- 239000000463 material Substances 0.000 claims description 15
- 230000005012 migration Effects 0.000 claims description 15
- 238000013508 migration Methods 0.000 claims description 15
- 238000001994 activation Methods 0.000 claims description 12
- 230000004913 activation Effects 0.000 claims description 12
- 239000007790 solid phase Substances 0.000 claims description 11
- 235000015220 hamburgers Nutrition 0.000 claims description 8
- 239000011159 matrix material Substances 0.000 claims description 8
- 230000007246 mechanism Effects 0.000 claims description 8
- 238000012360 testing method Methods 0.000 claims description 8
- 229910045601 alloy Inorganic materials 0.000 claims description 6
- 239000000956 alloy Substances 0.000 claims description 6
- 230000004888 barrier function Effects 0.000 claims description 6
- 238000007725 thermal activation Methods 0.000 claims description 6
- AZDRQVAHHNSJOQ-UHFFFAOYSA-N alumane Chemical group [AlH3] AZDRQVAHHNSJOQ-UHFFFAOYSA-N 0.000 claims description 5
- 229910052782 aluminium Inorganic materials 0.000 claims description 4
- 230000000694 effects Effects 0.000 claims description 4
- 230000004907 flux Effects 0.000 claims description 4
- 238000010587 phase diagram Methods 0.000 claims description 4
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 claims description 3
- 238000000137 annealing Methods 0.000 claims description 3
- 238000005452 bending Methods 0.000 claims description 3
- 210000002257 embryonic structure Anatomy 0.000 claims description 3
- 230000002349 favourable effect Effects 0.000 claims description 3
- 239000002904 solvent Substances 0.000 claims description 3
- 235000002779 Morchella esculenta Nutrition 0.000 claims 1
- 240000002769 Morchella esculenta Species 0.000 claims 1
- 238000011002 quantification Methods 0.000 abstract description 3
- 239000007787 solid Substances 0.000 abstract 1
- 210000001519 tissue Anatomy 0.000 description 20
- 238000011160 research Methods 0.000 description 6
- 210000005067 joint tissue Anatomy 0.000 description 3
- XKRFYHLGVUSROY-UHFFFAOYSA-N Argon Chemical compound [Ar] XKRFYHLGVUSROY-UHFFFAOYSA-N 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000007789 gas Substances 0.000 description 2
- 238000003913 materials processing Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 230000035882 stress Effects 0.000 description 2
- FGUUSXIOTUKUDN-IBGZPJMESA-N C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 Chemical compound C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 FGUUSXIOTUKUDN-IBGZPJMESA-N 0.000 description 1
- 230000032683 aging Effects 0.000 description 1
- 229910052786 argon Inorganic materials 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000012084 conversion product Substances 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000002050 diffraction method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- CKSRCDNUMJATGA-UHFFFAOYSA-N germanium platinum Chemical compound [Ge].[Pt] CKSRCDNUMJATGA-UHFFFAOYSA-N 0.000 description 1
- 230000002650 habitual effect Effects 0.000 description 1
- 238000011534 incubation Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000013021 overheating Methods 0.000 description 1
- 238000001556 precipitation Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000005204 segregation Methods 0.000 description 1
- 238000007711 solidification Methods 0.000 description 1
- 230000008023 solidification Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000005496 tempering Methods 0.000 description 1
- 238000005382 thermal cycling Methods 0.000 description 1
Images
Landscapes
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
一种焊接过程热影响区组织演变的模拟方法,它涉及一种焊接接头微观组织模拟方法,以解决目前焊缝微观组织演变的定量化,主要基于经验或者半经验的确定性模型或者解析计算,只进行组织含量的计算,而不能动态地反映组织形态、尺寸和分布的问题。本发明的方法是通过以下步骤实现的:步骤一:计算热影响区温度场;步骤二:根据温度场分布,计算不同位置晶粒在β相区以上的晶粒长大过程;步骤三:根据温度场计算获得的冷却速度和β相晶粒分布计算结果,计算连续冷却固态相变;本发明用于焊接接头微观组织模拟。
Description
技术领域
本发明涉及一种焊接接头微观组织模拟方法,具体涉及一种焊接过程热影响区组织演变的模拟方法。
背景技术
焊接热影响区的晶粒长大、淬硬脆化、时效软化或回火软化、脆性相析出等问题,一直都是导致焊接裂纹或接头性能下降的薄弱环节。因此,焊接热影响区组织演化的研究具有非常现实和重要的意义,其中,数值模拟是最有效的研究手段之一。
焊接接头组织模拟是一项涉及空间尺度、时间尺度、多物理场的巨大研究工程,难度很大,以往的焊接微观组织模拟工作侧重在凝固组织和晶粒长大模拟方面,研究成果相对成熟,对于固态相变尤其是焊接热影响区组织动态演变的模拟较少,且有一些问题始终未解决,如:焊缝微观组织演变的定量化研究,目前主要基于经验或者半经验的确定性模型或者解析计算,只进行组织含量的计算,不能动态地反映组织形态、尺寸和分布;数学模型尚不完善,特别是形核模型,一方面模型中的一些参数还取决于经验公式,另一方面没有考虑焊接过程中不同冷却速度引起的形核方式的改变;由于焊接过程具有高温过热和快速冷却的特点,热影响区的相变过程偏离平衡态,然而大多研究只局限于平衡相变或者固态相变一个部分,即组织动态演变和组织历史相关性的模拟罕有报道;目前关于相变的模拟没有考虑相变晶体学特征即惯习长大现象;对于特定的合金如钛合金,尚无工程用三元或多元合金平衡相图、等温转变曲线(TTT曲线)、连续冷却转变曲线(CCT曲线),热力学数据的困乏限制了对组织演变过程的精确描述;目前,还没有工程可用的焊接接头组织模拟可操作软件。
发明内容
本发明的目的是为了解决目前焊缝微观组织演变的定量化,主要基于经验或者半经验的确定性模型或者解析计算,只进行组织含量的计算,而不能动态地反映组织形态、尺寸和分布的问题,提供一种焊接过程热影响区组织演变的模拟方法。
本发明的焊接过程热影响区组织演变的模拟方法是通过以下步骤实现的:
步骤一:计算热影响区温度场:
通过求解导热微分方程,并且附上初始条件和换热边界条件,得到焊接过程的温度分布:
公式一:
将公式一计算的热影响区温度场分布,通过插值方法赋给晶粒长大的计算域,提取热影响区不同位置的热循环曲线,获得热影响区不同位置在加热阶段从低温α到高温β相变的实际温度、在冷却阶段从高温β到低温α相变的实际温度、峰值温度、晶粒长大温度区间的加热速度及冷却速度,其中,加热速度根据加热阶段的热循环曲线在晶粒长大区间的切线斜率获得,冷却速度根据冷却阶段的热循环曲线在晶粒长大区间的切线斜率获得,峰值温度即焊接粗晶区达到的最高温度;晶粒长大完毕后,开始发生高温β到低温α相变,此时冷却速度根据冷却阶段的热循环曲线在相变区间的切线斜率获得;
步骤二:根据温度场分布,计算不同位置晶粒在β相区以上的晶粒长大过程:
(1)、定义初始退火态β相组织,以及初始参数包括焊接工艺参数,热源参数,晶粒长大热力学参数;
(2)、计算晶粒长大速度,包括晶界曲率和界面能:晶界迁移速度V表示为:
公式二:V=MF
其中,M为晶界迁移率,是衡量晶界迁移快慢的参数,F为晶粒长大驱动力,晶粒长大驱动力F由界面能和晶界的弯曲程度求解:
公式三:
其中,γ为晶界的界面能,r为晶界的曲率半径,晶界长大驱动力亦即晶粒长大前后的自由能之差,
小角度晶界的界面能由相邻两晶粒间的取向差决定:
公式四:
其中,θ为任意晶界的取向差,γm为界面能的最大值,θm为界面能对应的晶界取向差,取θm=10°,
基于晶粒长大的原子跃迁机制,求得晶界迁移的速度:
公式五:
其中,晶界迁移率A2为晶粒2接纳原子的概率,N1为单位面积上处于能跳动的有利位置的平均原子数目,υ1为晶粒1中原子跳跃的频率,Vβ为β相的摩尔体积,Na为阿伏伽德罗常数,R为气体常数,T为热力学温度,为原子能量高于激活能的概率,ΔGA为原子从晶粒1分离出来靠热激活获得的激活能势垒;
(3)、建立晶粒长大动力学条件:假设原子内能的分布符合麦克斯威尔-波尔兹曼统计规律,原子内能的计算公式:
公式六:GT=-RTln(RAND)
其中,RAND()表示大于0并且小于等于1的随机数,GT为原子的内能,ln()为自然对数函数,
晶界元胞跨过晶界即长大转变的动力学能量条件为:
公式七:GT,i≥ΔGA-∑ΔGi
其中,GT,i为晶界原子i的内能,∑ΔGi为晶界原子i向临近晶粒转变前后的自由能之差,
在上述方法中,只要存在晶界迁移驱动力,晶粒就可以长大,即晶粒长大的条件为ΔG>0,或者是晶界原子转变前后的自由能与内能变化之和大于0,即热力学能量条件ΔG+ΔGT≥0,通过公式一至公式七建立了晶粒长大的动力学能量条件,即当晶界原子转变前后的内能与自由能变化,大于该原子从当前所属晶粒分离出来的跳跃激活能势垒时,晶粒才会长大,体现了原子跳跃是一个热激活过程;
(4)、计算晶粒长大概率:为了体现晶粒长大过程的不可逆性,只有达到最小能量条件的晶粒长大过程才是最稳定和最该发生的;为体现随机性,以最大晶粒长大速度为参照计算晶界原子的转变概率:
公式八:
其中,Mi为原子i所在晶粒的晶界迁移率,Fi为原子i的晶粒长大驱动力,Mmax为最大晶界迁移率,Fmax为最大晶粒长大驱动力,当晶界原子可以向一个以上临近晶粒转变时,实际上只会向转变概率最大的晶粒迁移;
(5)、显示温度场分布以及焊接接头晶粒形态和大小拓扑形貌;
步骤三:根据温度场计算获得的冷却速度和β相晶粒分布计算结果,计算连续冷却固态相变:
(1)、定义初始参数,包括溶质扩散系数、界面迁移率、冷却速度、过冷度、溶质原子在β和α相中的平衡溶解度,设置焊接热影响区晶粒长大计算结果为初始组织:
利用Thermα-Calc相图计算软件进行热力学计算,获得相变特征温度,即β和α相共存温度区间,以及溶质原子在该温度区间的溶解度,由于焊接过程快速加热和冷却,偏离热力学平衡状态,模拟β相向α相转变需考虑过冷度与特定温度、压力下的稳态平衡相比,溶质在两相中的实际溶解度按照相变平衡温度和相变过冷度计算,采用热模拟试验测量合金的CCT曲线,通过CCT曲线,获得不同冷却速度下连续冷却非平衡相变的起始温度和终止温度,将不同冷却速度的过冷度和相变温度区间、时间作为相变计算的输入条件,设置焊接热影响区晶粒长大计算结果为初始组织;
(2)、计算相变的形核率,设置形核:
晶界形核稳态形核率I稳为:
其中,Nv为单位体积核胚的数量,k为波尔兹曼常数,h为普朗克常量,Q为溶质原子在β相中的热扩散激活能,ΔG*为形核驱动力,
当过冷度较大或者冷却速度较大时,认为形核在相变初期瞬间完成且形核位置从晶界向晶内均匀形核转移;
(3)、计算生长态元胞的α相长大相变驱动力、速度,计算α相前沿生长态元胞的α相转变分数、α相溶质浓度、β相溶质浓度及整个计算域的溶质场分布:
相界面移动速度VP:
公式十:VP=MPFP/Vβ
其中,MP是相界面迁移率,FP为相变驱动力,
当以理想稀溶液处理时,相变驱动力FP为:
公式十一:
其中,Xi是界面处β相中溶质原子铝的摩尔百分数,XAl β为β/α相平衡时β相的溶质浓度;
对于界面或混合控制型相变,由相长大速度获得界面处母相中铝原子的摩尔百分数:
公式十二:
β相溶质扩散方程为:
公式十三:
α相溶质扩散方程为:
公式十四:
根据界面处溶质守恒,获得:
公式十五:
其中,XAl α为β/α相平衡态α相的溶质浓度,VP(XAl α-Xi)表示单位面积单位时间相界面处从β相脱溶析出的溶质原子量,表示界面处β相中溶质梯度引起的扩散流量,为界面处β相的溶质浓度梯度,n为β/α相界面的法向矢量;
因此,对于混合型相变,相界面移动速度VP在x轴方向的差分表达式为:
公式十六:
其中VP,x(i,j)为计算域中坐标为(i,j)的节点在x方向的α相长大速度,Δx为空间步长,VP,x为上一时间步α相尖端在x方向的长大速度,Cβ(i-1,j)为坐标(i-1,j)的节点的β相溶质浓度,fβ(i-1,j)为坐标(i-1,j)的节点的β相体积分数,Cβ(i+1,j)为坐标(i+1,j)的节点的β相溶质浓度,fβ(i+1,j)为坐标(i+1,j)的节点的β相体积分数,Cα(i-1,j)为坐标(i-1,j)的节点的α相溶质浓度,fα(i-1,j)为坐标(i-1,j)的节点的α相体积分数,Cα(i+1,j)为坐标(i+1,j)的节点的α相溶质浓度,fα(i+1,j)为坐标(i+1,j)的节点的α相体积分数;
(4)、计算α相尖端曲率,修正热力学平衡β相溶质浓度:根据Thermo-Calc软件计算可知,当热力学温度一定时,α相的平衡溶解度XAl α确定,此时,当α相具有无限大曲率半径时,溶质原子在β相中的热力学平衡溶解度为XAl β;当α相尖端曲率为时,还需考虑β/α相界面曲率效应,即溶质原子在β相中的热力学平衡溶解度应为即将XAl β修正为对于稀溶液应满足由此:
公式十七:
公式十八:
其中,fα生长态元胞的α相分数,N为元胞自动机的邻居元胞数目,对于冯·诺依曼型邻居关系N=4,对于摩尔型邻居关系型邻居关系N=8;
Gibbs-Thomson系数:
公式十九:
其中,σ为β/α的相界界面能;
生长态元胞α相分数:
公式二十: 并且fα n=fα 0+δfα
其中,δt为计算时间步长,Vx为生长态元胞在x方向的相界面移动速度,Vy为生长态元胞在y方向的相界面移动速度,fα 0为上一时间步的生长态元胞的α相分数,fα n分别为当前时间步的生长态元胞的α相分数,δfα为生长态元胞的α相分数增量;
(5)、计算各向异性相界面迁移率:钛合金从β相区连续冷却时,α相与β相基体之间存在严格的伯格斯(Burgers)晶体学取向关系,以界面迁移率的各向异性表征,在择优生长方向和与其垂直方向的界面迁移率分别记为M1和M2,根据相界迁移的微观物理机制,获得界面迁移率的表达式:
公式二十一:M=M0exp(-Qm/RT)
公式二十二:
其中,M0为指数前因子,Vα为α相的摩尔体积、vD≈5×1013/s为德拜频率;由于M1/M2的大小直接决定析出新相的长细比,取M1/M2=10、100或1000,分别进行析出α相形态的模拟,并且通过扫描电镜观测相同焊接工艺条件下热影响区的组织形态,将析出α相形态的模拟结果与实验结果对比,两者最吻合时采用的M1/M2最佳;
(6)、计算晶内不同择优取向的α相长大形貌:新相长大方向与母相晶粒的晶体位向有关,当母相晶体取向为优先取向时,按照与β相晶粒的晶体取向一致方向长大,当β相晶粒的晶体取向不是优先取向时,则与β相晶粒的晶体取向呈一特定角度长大,由于对于钛合金而言,β-α相变存在多个惯习面和惯习方向,所以每个β相晶粒内α相长大的位相取向并不唯一,因此,α相生长方向与β相的晶界呈一定角度,最终经历连续冷却非平衡相变后,在焊接热影响区形成网篮组织;
(7)、更新元胞状态,显示组织形态及溶质场分布、不同位置不同时间的α相变分数,并以文件形式保存。
本发明具有以下有益效果:
一、本发明的方法是利用计算机在虚拟条件下进行的。采用本发明的焊接热影响区晶粒长大和固态相变组织演变元胞自动机模型,实现了焊接热影响区不同位置尤其是粗晶区的组织演变模拟,能够动态再现焊接全过程的组织形态,定量预测焊接接头最终晶粒尺寸、相含量、组成相的特征尺寸,如片状α相的宽度,α相集束的宽度等,从而建立了工艺-组织与力学性能之间的映射关系,为制定钛合金的焊接工艺,改善焊接接头的显微组织和力学性能提供理论基础,因此,本发明解决了目前焊缝微观组织演变的定量化,主要基于经验或者半经验的确定性模型或者解析计算,只进行组织含量的计算,而不能动态地反映组织形态、尺寸和分布的问题。
二、由于本发明采用元胞自动机法能较精确地模拟材料的介观和微观组织,从而实现固态相变元胞自动机法模拟与介观尺度上晶粒长大元胞自动机法模拟的“无缝连接”,并以模拟晶粒长大为桥梁建立与宏观温度场之间的联系,实现宏观-介观-微观多尺度耦合模拟,揭示焊接过程组织遗传规律,即焊接热影响区晶粒过热粗大对后续接头组织演化规律的影响。
三、本发明通过数值模拟、相变理论和材料热力学多学科结合,在不同尺度和不同位置上再现焊接接头组织演化的形貌,定量描述晶体学、热力学和动力学机制,并且可预测TTT曲线,填充了航空新材料钛合金TTT/CCT图的空白。
四、组织场、宏观应力应变场和温度场三者相互作用,共同影响焊接接头性能,通过组织模拟能定量计算溶质偏析,并能建立微观局部应力场和组织场的耦合模型,进而结合断裂力学和损伤力学,分析产生焊接缺陷如冷裂纹产生的冶金条件和力学条件。
附图说明
图1为具体实施方式一中焊接热影响区组织演变的计算方法的流程图;图2为焊接热影响区温度场分布的模拟结果图;图3为焊接热影响区距离焊缝不同位置的热循环曲线图;图4为具体实施方式一中步骤二的流程图;图5为具体实施方式一中步骤三的流程图;图6为TA15钛合金在焊接热影响区晶粒长大和固态相变组织形貌模拟结果图;图7为TA15钛合金在焊接热影响区不同位置随时间变化的最终相变分数计算结果图;图8为TA15钛合金在焊接热影响区不同大小的晶粒内部不同角度的α相长大形貌图。
具体实施方式
具体实施方式一:结合图1~图5说明本实施方式,本发明的焊接过程热影响区组织演变的模拟方法是通过以下步骤实现的:
步骤一:计算热影响区温度场:
通过求解导热微分方程,并且附上初始条件和换热边界条件,得到焊接过程的温度分布:
公式一:
将公式一计算的热影响区温度场分布,如图2所示,通过插值方法赋给晶粒长大的计算域,提取热影响区不同位置的热循环曲线,如图3所示,获得热影响区不同位置在加热阶段从低温α到高温β相变的实际温度、在冷却阶段从高温β到低温α相变的实际温度、峰值温度、晶粒长大温度区间的加热速度及冷却速度,其中,加热速度根据加热阶段的热循环曲线在晶粒长大区间的切线斜率获得,冷却速度根据冷却阶段的热循环曲线在晶粒长大区间的切线斜率获得,峰值温度即焊接粗晶区达到的最高温度;晶粒长大完毕后,开始发生高温β到低温α相变,此时冷却速度根据冷却阶段的热循环曲线在相变区间的切线斜率获得;
步骤二:根据温度场分布,计算不同位置晶粒在β相区以上的晶粒长大过程:如图4所示,
(1)、定义初始退火态β相组织,以及初始参数包括焊接工艺参数,热源参数,晶粒长大热力学参数;
(2)、计算晶粒长大速度,包括晶界曲率和界面能:晶界迁移速度V表示为:
公式二:V=MF
其中,M为晶界迁移率,是衡量晶界迁移快慢的参数,F为晶粒长大驱动力,
晶粒长大驱动力F由界面能和晶界的弯曲程度求解:
公式三:
其中,γ为晶界的界面能,r为晶界的曲率半径,晶界长大驱动力亦即晶粒长大前后的自由能之差,
小角度晶界的界面能由相邻两晶粒间的取向差决定:
公式四:
其中,θ为任意晶界的取向差,γm为界面能的最大值,θm为界面能对应的晶界取向差,取θm=10°,
基于晶粒长大的原子跃迁机制,求得晶界迁移的速度:
公式五:
其中,晶界迁移率A2为晶粒2接纳原子的概率,N1为单位面积上处于能跳动的有利位置的平均原子数目,υ1为晶粒1中原子跳跃的频率,Vβ为β相的摩尔体积,Nα为阿伏伽德罗常数,R为气体常数,T为热力学温度,为原子能量高于激活能的概率,ΔGA为原子从晶粒1分离出来靠热激活获得的激活能势垒;
(3)、建立晶粒长大条件:假设原子内能的分布符合麦克斯威尔-波尔兹曼统计规律,原子内能的计算公式:
公式六:GT=-RTln(RAND)
其中,RAND()表示大于0并且小于等于1的随机数,GT为原子的内能,ln()为自然对数函数,
晶界元胞跨过晶界即长大转变的动力学能量条件为:
公式七:GT,i≥ΔGA-∑ΔGi
其中,GT,i为晶界原子i的内能,∑ΔGi为晶界原子i向临近晶粒转变前后的自由能之差,
在上述方法中,只要存在晶界迁移驱动力,晶粒就可以长大,即晶粒长大的条件为ΔG>0,或者是晶界原子转变前后的自由能与内能变化之和大于0,即热力学能量条件ΔG+ΔGT≥0,通过公式一至公式七建立了晶粒长大的动力学能量条件,即当晶界原子转变前后的内能与自由能变化,大于该原子从当前所属晶粒分离出来的跳跃激活能势垒时,晶粒才会长大,体现了原子跳跃是一个热激活过程;
(4)、计算晶粒长大概率:为了体现晶粒长大过程的不可逆性,只有达到最小能量条件的晶粒长大过程才是最稳定和最该发生的;为体现随机性,以最大晶粒长大速度为参照计算晶界原子的转变概率:
公式八:
其中,Mi为原子i所在晶粒的晶界迁移率,Fi为原子i的晶粒长大驱动力,Mmax为最大晶界迁移率,Fmax为最大晶粒长大驱动力,当晶界原子可以向一个以上临近晶粒转变时,实际上只会向转变概率最大的晶粒迁移;
(5)、显示温度场分布以及焊接接头晶粒形态和大小拓扑形貌;
步骤三:根据温度场计算获得的冷却速度和β相晶粒分布计算结果,计算连续冷却固态相变:如图5所示,
(1)、定义初始参数,包括溶质扩散系数、界面迁移率、冷却速度、过冷度、溶质原子在β和α相中的平衡溶解度,设置焊接热影响区晶粒长大计算结果为初始组织:
利用Thermo-Calc相图计算软件进行热力学计算,获得相变特征温度,即β和α相共存温度区间,以及溶质原子在该温度区间的溶解度,由于焊接过程快速加热和冷却,偏离热力学平衡状态,模拟β相向α相转变需考虑过冷度与特定温度、压力下的稳态平衡相比,溶质在两相中的实际溶解度按照相变平衡温度和相变过冷度计算,采用热模拟试验测量合金的CCT曲线,通过CCT曲线,获得不同冷却速度下连续冷却非平衡相变的起始温度和终止温度,将不同冷却速度的过冷度和相变温度区间、时间作为相变计算的输入条件,设置焊接热影响区晶粒长大计算结果为初始组织;
(2)、计算相变的形核率,设置形核:
晶界形核稳态形核率I稳为:
其中,Nv为单位体积核胚的数量,k为波尔兹曼常数,h为普朗克常量,Q为溶质原子在β相中的热扩散激活能,ΔG*为形核驱动力,
当过冷度较大或者冷却速度较大时,认为形核在相变初期瞬间完成且形核位置从晶界向晶内均匀形核转移;
(3)、计算生长态元胞的α相长大相变驱动力、速度,计算α相前沿生长态元胞的α相转变分数、α相溶质浓度、β相溶质浓度及整个计算域的溶质场分布:
相界面移动速度VP:
公式十:VP=MPFP/Vβ
其中,MP是相界面迁移率,FP为相变驱动力,
当以理想稀溶液处理时,相变驱动力FP为:
公式十一:
其中,Xi是界面处β相中溶质原子铝的摩尔百分数,XAl β为β/α相平衡时β相的溶质浓度;
对于界面或混合控制型相变,由相长大速度获得界面处母相中铝原子的摩尔百分数:
公式十二:
β相溶质扩散方程为:
公式十三:
α相溶质扩散方程为:
公式十四:
根据界面处溶质守恒,获得:
公式十五:
其中,XAl α为β/α相平衡态α相的溶质浓度,VP(XAl α-Xi)表示单位面积单位时间相界面处从β相脱溶析出的溶质原子量,表示界面处β相中溶质梯度引起的扩散流量,为界面处β相的溶质浓度梯度,n为β/α相界面的法向矢量;
因此,对于混合型相变,相界面移动速度VP在x轴方向的差分表达式为:
公式十六:
其中VP,x(i,j)为计算域中坐标为(i,j)的节点在x方向的α相长大速度,Δx为空间步长,VP,x为上一时间步α相尖端在x方向的长大速度,Cβ(i-1,j)为坐标(i-1,j)的节点的β相溶质浓度,fβ(i-1,j)为坐标(i-1,j)的节点的β相体积分数,Cβ(i+1,j)为坐标(i+1,j)的节点的β相溶质浓度,fβ(i+1,j)为坐标(i+1,j)的节点的β相体积分数,Cα(i-1,j)为坐标(i-1,j)的节点的α相溶质浓度,fα(i-1,j)为坐标(i-1,j)的节点的α相体积分数,Cα(i+1,j)为坐标(i+1,j)的节点的α相溶质浓度,fα(i+1,j)为坐标(i+1,j)的节点的α相体积分数;
(4)、计算α相尖端曲率,修正热力学平衡β相溶质浓度:根据Thermo-Calc软件计算可知,当热力学温度一定时,α相的平衡溶解度XAl α确定,此时,当α相具有无限大曲率半径时,溶质原子在β相中的热力学平衡溶解度为XAl β;当α相尖端曲率为时,还需考虑β/α相界面曲率效应,即溶质原子在β相中的热力学平衡溶解度应为即将XAl β修正为对于稀溶液应满足由此:
公式十七:
其中,为α相尖端尖端曲率,Γ为Gibbs-Thomson系数,α相尖端尖端曲率:
公式十八:
其中,fα生长态元胞的α相分数,N为元胞自动机的邻居元胞数目,对于冯·诺依曼型邻居关系N=4,对于摩尔型邻居关系型邻居关系N=8;
Gibbs-Thomson系数:
公式十九:
其中,σ为β/α的相界界面能;
生长态元胞α相分数:
公式二十: 并且fα n=fα 0+δfα
其中,δt为计算时间步长,Vx为生长态元胞在x方向的相界面移动速度,Vy为生长态元胞在y方向的相界面移动速度,fα 0为上一时间步的生长态元胞的α相分数,fα n分别为当前时间步的生长态元胞的α相分数,δfα为生长态元胞的α相分数增量;
(5)、计算各向异性相界面迁移率:钛合金从β相区连续冷却时,α相与β相基体之间存在严格的伯格斯(Burgers)晶体学取向关系,以界面迁移率的各向异性表征,在择优生长方向和与其垂直方向的界面迁移率分别记为M1和M2,根据相界迁移的微观物理机制,获得界面迁移率的表达式:
公式二十一:M=M0exp(-Qm/RT)
公式二十二:
其中,M0为指数前因子,Vα为α相的摩尔体积、vD≈5×1013/s为德拜频率;由于M1/M2的大小直接决定析出新相的长细比,取M1/M2=10、100或1000,分别进行析出α相形态的模拟,并且通过扫描电镜观测相同焊接工艺条件下热影响区的组织形态,将析出α相形态的模拟结果与实验结果对比,两者最吻合时采用的M1/M2最佳;
(6)、计算晶内不同择优取向的α相长大形貌:新相长大方向与母相晶粒的晶体位向有关,当母相晶体取向为优先取向时,按照与β相晶粒的晶体取向一致方向长大,当β相晶粒的晶体取向不是优先取向时,则与β相晶粒的晶体取向呈一特定角度长大,由于对于钛合金而言,β-α相变存在多个惯习面和惯习方向,所以每个β相晶粒内α相长大的位相取向并不唯一,因此,α相生长方向与β相的晶界呈一定角度,最终经历连续冷却非平衡相变后,在焊接热影响区形成网篮组织;
(7)、更新元胞状态,显示组织形态及溶质场分布、不同位置不同时间的α相变分数,并以文件形式保存。
不同角度α相的长大形貌的计算按照占小红提出的对角线模拟法则进行,(参见《Cellular automaton simulation of grain growth with different orientation angles duringsolidification process》:X.H.Zhan,Y.H.Wei,Z.B.Dong.Journal of Materials ProcessingTechnology.Vol.12(2007):130)。
本发明建立了在焊接工艺下钛合金的非平衡连续冷却相变模型,预测热影响区晶粒长大和粗大β相晶粒在冷却时经历相变过程的最终相变分数以及显微组织状态,并且定量表征冷却速度、界面能各向异性对组织形貌、新相长大动力学和组织转变分数的影响规律,为制定合理的钛合金焊接工艺提供理论指导。
对于建立的晶粒长大和固态相变计算模型,现以目前应用较广泛的TA15钛合金为例,进行CCT曲线测定实验和相关热力学计算和组织演变模拟,该TA15钛合金的材料成分原子百分含量如下表一所示。
表一TA15的化学成分(mole.%)
其热影响区固态相变数值模拟的方法通过一下步骤实现:
步骤a、针对TA15钛合金,采用最新版本的Thermo-Calc软件开展热力学计算,为简便,以Al当量和Mo当量为参考,将TA15钛合金近似为Ti-Al-Mo三元系合金,采用正规稀溶体近似描述β相和α相,结合材料热力学理论,获得Ti-Al-Mo三元系Al和Mo在母相和新相溶解度随着温度的变化,两相共存温度区间为980℃-774℃,由于焊接过程快速加热和冷却,偏离热力学平衡状态,模拟β相向α相转变需考虑过冷度,即焊接条件下系统处于亚稳平衡态,与特定温度、压力下的稳态平衡相比,平衡成分范围或温度都发生了明显的变化;
步骤b、采用Gleeble1500热模拟试验测量TA15钛合金的CCT曲线,通过CCT曲线,定量的计算连续冷却时相变孕育时间,相变开始时间和温度、结束时间和温度,以及最终转变产物的各相体积分数比,并将热循环试验的热循环曲线作为模拟条件,从而提高模拟结果的可信度。试验棒材的规格为Φ6mm×130mm(可接受长度在100mm~150mm之间,粗糙度为Ra3.2,径向正偏差);热循环工艺如下:加热速度为90℃/s、峰值温度为1360℃,所以β相晶粒粗大,冷却速度依次为18.3℃/s、16.7℃/s、9.9℃/s、5.0℃/s、2.5℃/s、1.5℃/s。整个试验过程在氩气保护或者真空状态下进行,以防止高温下试样被氧化,并由铂金-铂锗热电偶插入预先激光加工的Φ1mm×2mm深孔内测量加热区温度变化,测量得到不同冷却速度时试验棒料径向收缩量随温度的变化,通过切线法获得相变起始和终止温度。步骤c、计算β-α相变形核率。晶界形核稳态形核率为:
其中,Nv为单位体积核胚(原子团)的数量。
计算表明,当相变温度在750℃~900℃之间时,形核率在0(nuclei/m/s)~0.01(nuclei/m/s)之间,与温度呈抛物线关系变化,成分固定时,在850℃左右形核率达最大值,温度过高或过低都会造成形核率的下降,形核方式即是均匀形核还是晶界形核,以及是瞬间形核还是连续形核,依据过冷度和冷却速度大小而定;当发生连续冷却相变时,与等温相变形核率计算不同,近似认为形核率与冷却速度呈线性关系变化,由于相变过冷度大,认为形核在相变初期瞬间完成。
步骤d、计算相变速度和溶质场分布:
相界面移动速度VP:VP=MPFP/Vβ
其中,MP是相界面迁移率,FP为相变驱动力。
其中,Xi是界面处母相中铝原子的摩尔百分数,XAl β为母/新相平衡态母相的溶质浓度;
对于界面或混合控制型相变,由相长大速度获得界面处β相中铝原子的摩尔百分数:
β相溶质扩散方程:
α相溶质扩散方程:
其中,Dβ、Dα分别为β相和α相溶质扩散系数,Cβ和Cα为β相和α相的溶质浓度。
根据界面处溶质守恒,获得:
其中,XAl β为β/α相平衡时β相的溶质浓度,等式左侧表示单位面积单位时间相界面处从β相脱溶析出的溶质原子量,为界面处β相中溶质梯度引起的扩散流量,为β相的溶质浓度梯度,n为β/α相界面的法向矢量。
因此,对于混合型相变,速度VP在x轴方向的差分表达式为:
α相尖端尖端曲率:
Gibbs-Thomson系数:
生长态元胞α相分数: 并且fα n=fα 0+δfα
步骤f、计算各向异性相界面迁移率:钛合金从β相区连续冷却时,α相与β相基体之间存在严格的伯格斯(Burgers)晶体学取向关系,因此,以界面迁移率的各向异性表征,在择优生长方向和与其垂直方向的界面迁移率分别记为M1和M2,根据相界迁移的微观物理机制,获得相界面迁移率的表达式:MP=M0exp(-Q/RT)
其中,为指数前因子,Vα为α相的摩尔体积,vD≈5×1013/s为德拜频率;由于M1/M2的大小直接决定析出新相的长细比,取M1/M2=10、100或1000,分别进行析出α相形态的模拟,并且通过扫描电镜观测相同焊接工艺条件下热影响区的组织形态,将析出α相形态的模拟结果与实验结果对比,两者最吻合时采用的M1/M2最佳,为M1/M2=100。随着冷却速度增大,在以上关于界面迁移率的表达式基础上,近似认为相界面迁移率与冷却速度呈线性正比关系。
步骤g、根据TIG焊接工艺下(如表二所示)的温度场和热循环模拟结果,开展焊接热影响区组织演变计算,计算域为200×150μm,空间网格数为400×300和空间步长为0.5μm。焊接热影响区晶粒长大和固态相变组织形貌模拟组织形貌模拟结果如图6所示,焊接热影响区不同位置随时间变化的最终相变分数如图7所示;
表二TIG焊接工艺参数
步骤h、计算晶内不同择优取向的α相长大形貌:固态相变需遵循晶体学位向关系,α相从β相析出长大时,按照与惯习面平行的方向优先选择生长,所以,α相长大方向与β相晶粒的晶体位向有关,当β相晶体取向为优先取向时,按照与β相晶体取向一致方向长大;当β相晶体取向不是优先取向时则与β相晶体取向为一特定角度长大;在模拟时,当设置晶界方向固定为x轴方向,在任意晶粒内的α相长大,即在面心立方体其它晶向的β相内长大,长大方向与x轴呈一定夹角;由于存在多个惯习面和方向,所以每个晶粒内α相长大的取向并不唯一,形成网篮组织,如图8所示。
Claims (1)
1.一种焊接过程热影响区组织演变的模拟方法,其特征在于:所述方法是通过以下步骤实现的:
步骤一:计算热影响区温度场:
通过求解导热微分方程,并且附上初始条件和换热边界条件,得到焊接过程的温度分布:
公式一:
将公式一计算的热影响区温度场分布,通过插值方法赋给晶粒长大的计算域,提取热影响区不同位置的热循环曲线,获得热影响区不同位置在加热阶段从低温α到高温β相变的实际温度、在冷却阶段从高温β到低温α相变的实际温度、峰值温度、晶粒长大温度区间的加热速度及冷却速度,其中,加热速度根据加热阶段的热循环曲线在晶粒长大区间的切线斜率获得,冷却速度根据冷却阶段的热循环曲线在晶粒长大区间的切线斜率获得,峰值温度即焊接粗晶区达到的最高温度;晶粒长大完毕后,开始发生高温β到低温α相变,此时冷却速度根据冷却阶段的热循环曲线在相变区间的切线斜率获得;
步骤二:根据温度场分布,计算不同位置晶粒在β相区以上的晶粒长大过程:
(1)、定义初始退火态β相组织,以及初始参数包括焊接工艺参数,热源参数,晶粒长大热力学参数;
(2)、计算晶粒长大速度,包括晶界曲率和界面能:晶界迁移速度V表示为:
公式二:V=MF
其中,M为晶界迁移率,是衡量晶界迁移快慢的参数,F为晶粒长大驱动力,晶粒长大驱动力F由界面能和晶界的弯曲程度求解:
公式三:
其中,γ为晶界的界面能,r为晶界的曲率半径,晶界长大驱动力亦即晶粒长大前后的自由能之差,
小角度晶界的界面能由相邻两晶粒间的取向差决定:
公式四:
其中,θ为任意晶界的取向差,γm为界面能的最大值,θm为界面能对应的晶界取向差,取θm=10°,
基于晶粒长大的原子跃迁机制,求得晶界迁移的速度:
公式五:
其中,晶界迁移率A2为晶粒2接纳原子的概率,N1为单位面积上处于能跳动的有利位置的平均原子数目,υ1为晶粒1中原子跳跃的频率,Vβ为β相的摩尔体积,Na为阿伏伽德罗常数,R为气体常数,T为热力学温度,为原子能量高于激活能的概率,ΔGA为原子从晶粒1分离出来靠热激活获得的激活能势垒;
(3)、建立晶粒长大动力学条件:假设原子内能的分布符合麦克斯威尔-波尔兹曼统计规律,原子内能的计算公式:
公式六:GT=-RTln(RAND)
其中,RAND()表示大于0并且小于等于1的随机数,GT为原子的内能,ln()为自然对数函数,
晶界元胞跨过晶界即长大转变的动力学能量条件为:
公式七:GT,i≥ΔGA-∑ΔGi
其中,GT,i为晶界原子i的内能,∑ΔGi为晶界原子i向临近晶粒转变前后的自由能之差,
在上述步骤中,只要存在晶界迁移驱动力,晶粒就可以长大,即晶粒长大的条件为ΔG>0,或者是晶界原子转变前后的自由能与内能变化之和大于0,即热力学能量条件ΔG+ΔGT≥0,通过公式一至公式七建立了晶粒长大的动力学能量条件,即当晶界原子转变前后的内能与自由能变化,大于该原子从当前所属晶粒分离出来的跳跃激活能势垒时,晶粒才会长大,体现了原子跳跃是一个热激活过程;
(4)、计算晶粒长大概率:为了体现晶粒长大过程的不可逆性,只有达到最小能量条件的晶粒长大过程才是最稳定和最该发生的;为体现随机性,以最大晶粒长大速度为参照计算晶界原子的转变概率:
公式八:
其中,Mi为原子i所在晶粒的晶界迁移率,Fi为原子i的晶粒长大驱动力,Mmax为最大晶界迁移率,Fmax为最大晶粒长大驱动力,当晶界原子可以向一个以上临近晶粒转变时,实际上只会向转变概率最大的晶粒迁移;
(5)、显示温度场分布以及焊接接头晶粒形态和大小拓扑形貌;
步骤三:根据温度场计算获得的冷却速度和β相晶粒分布计算结果,计算连续冷却固态相变:
(1)、定义初始参数,包括溶质扩散系数、界面迁移率、冷却速度、过冷度、溶质原子在β和α相中的平衡溶解度,设置焊接热影响区晶粒长大计算结果为初始组织:
利用Thermo-Calc相图计算软件进行热力学计算,获得相变特征温度,即β和α相共存温度区间,以及溶质原子在该温度区间的溶解度,由于焊接过程快速加热和冷却,偏离热力学平衡状态,模拟β相向α相转变需考虑过冷度与特定温度、压力下的稳态平衡相比,溶质在两相中的实际溶解度按照相变平衡温度和相变过冷度计算,采用热模拟试验测量合金的CCT曲线,通过CCT曲线,获得不同冷却速度下连续冷却非平衡相变的起始温度和终止温度,将不同冷却速度的过冷度和相变温度区间、时间作为相变计算的输入条件,设置焊接热影响区晶粒长大计算结果为初始组织;
(2)、计算相变的形核率,设置形核:
晶界形核稳态形核率I稳为:
其中,Nv为单位体积核胚的数量,k为波尔兹曼常数,h为普朗克常量,Q为溶质原子在β相中的热扩散激活能,ΔG*为形核驱动力,
当过冷度较大或者冷却速度较大时,认为形核在相变初期瞬间完成且形核位置从晶界向晶内均匀形核转移;
(3)、计算生长态元胞的α相长大相变驱动力、速度,计算α相前沿生长态元胞的α相转变分数、α相溶质浓度、β相溶质浓度及整个计算域的溶质场分布:
相界面移动速度VP:
公式十:VP=MPFP/Vβ
其中,MP是相界面迁移率,FP为相变驱动力,
当以理想稀溶液处理时,相变驱动力FP为:
公式十一:
其中,Xi是界面处β相中溶质原子铝的摩尔百分数,XAl β为β/α相平衡时β相的溶质浓度;
对于界面或混合控制型相变,由相长大速度获得界面处母相中铝原子的摩尔百分数:
公式十二:
β相溶质扩散方程为:
公式十三:
α相溶质扩散方程为:
公式十四:
根据界面处溶质守恒,获得:
公式十五:
其中,XAl α为β/α相平衡态α相的溶质浓度,VP(XAl α-Xi)表示单位面积单位时间相界面处从β相脱溶析出的溶质原子量,表示界面处β相中溶质梯度引起的扩散流量,为界面处β相的溶质浓度梯度,n为β/α相界面的法向矢量;
因此,对于混合型相变,相界面移动速度VP在x轴方向的差分表达式为:
公式十六:
其中VP,x(i,j)为计算域中坐标为(i,j)的节点在x方向的α相长大速度,Δx为空间步长,VP,x为上一时间步α相尖端在x方向的长大速度,Cβ(i-1,j)为坐标(i-1,j)的节点的β相溶质浓度,fβ(i-1,j)为坐标(i-1,j)的节点的β相体积分数,Cβ(i+1,j)为坐标(i+1,j)的节点的β相溶质浓度,fβ(i+1,j)为坐标(i+1,j)的节点的β相体积分数,Cα(i-1,j)为坐标(i-1,j)的节点的α相溶质浓度,fα(i-1,j)为坐标(i-1,j)的节点的α相体积分数,Cα(i+1,j)为坐标(i+1,j)的节点的α相溶质浓度,fα(i+1,j)为坐标(i+1,j)的节点的α相体积分数;
(4)、计算α相尖端曲率,修正热力学平衡β相溶质浓度:根据Thermo-Calc软件计算可知,当热力学温度一定时,α相的平衡溶解度XAl α确定,此时,当α相具有无限大曲率半径时,溶质原子在β相中的热力学平衡溶解度为XAl β;当α相尖端曲率为时,还需考虑β/α相界面曲率效应,即溶质原子在β相中的热力学平衡溶解度应为即将XAl β修正为对于稀溶液应满足由此:
公式十七:
α相尖端尖端曲率:
公式十八:
其中,fα生长态元胞的α相分数,N为元胞自动机的邻居元胞数目,对于冯·诺依曼型邻居关系N=4,对于摩尔型邻居关系型邻居关系N=8;
Gibbs-Thomson系数:
公式十九:
其中,σ为β/α的相界界面能;
生长态元胞α相分数:
公式二十: 并且fα n=fα 0+δfα
其中,δt为计算时间步长,Vx为生长态元胞在x方向的相界面移动速度,Vy为生长态元胞在y方向的相界面移动速度,fα 0为上一时间步的生长态元胞的α相分数,fα n分别为当前时间步的生长态元胞的α相分数,δfα为生长态元胞的α相分数增量;
(5)、计算各向异性相界面迁移率:钛合金从β相区连续冷却时,α相与β相基体之间存在严格的伯格斯(Burgers)晶体学取向关系,以界面迁移率的各向异性表征,在择优生长方向和与其垂直方向的界面迁移率分别记为M1和M2,根据相界迁移的微观物理机制,获得界面迁移率的表达式:
公式二十一:M=M0exp(-Qm/RT)
公式二十二:
其中,M0为指数前因子,Vα为α相的摩尔体积、vD≈5×1013/s为德拜频率;由于M1/M2的大小直接决定析出新相的长细比,取M1/M2=10、100或1000,分别进行析出α相形态的模拟,并且通过扫描电镜观测相同焊接工艺条件下热影响区的组织形态,将析出α相形态的模拟结果与实验结果对比,两者最吻合时采用的M1/M2最佳;
(6)、计算晶内不同择优取向的α相长大形貌:新相长大方向与母相晶粒的晶体位向有关,当母相晶体取向为优先取向时,按照与β相晶粒的晶体取向一致方向长大,当β相晶粒的晶体取向不是优先取向时,则与β相晶粒的晶体取向呈一特定角度长大,由于对于钛合金而言,β-α相变存在多个惯习面和惯习方向,所以每个β相晶粒内α相长大的位相取向并不唯一,因此,α相生长方向与β相的晶界呈一定角度,最终经历连续冷却非平衡相变后,在焊接热影响区形成网篮组织;
(7)、更新元胞状态,显示组织形态及溶质场分布、不同位置不同时间的α相变分数,并以文件形式保存。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210247239.XA CN102750425B (zh) | 2012-07-17 | 2012-07-17 | 一种焊接过程热影响区组织演变的模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210247239.XA CN102750425B (zh) | 2012-07-17 | 2012-07-17 | 一种焊接过程热影响区组织演变的模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102750425A true CN102750425A (zh) | 2012-10-24 |
CN102750425B CN102750425B (zh) | 2014-09-03 |
Family
ID=47030607
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210247239.XA Expired - Fee Related CN102750425B (zh) | 2012-07-17 | 2012-07-17 | 一种焊接过程热影响区组织演变的模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102750425B (zh) |
Cited By (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103246774A (zh) * | 2013-05-13 | 2013-08-14 | 天津大学 | 数值模拟p92钢管焊接热影响区宽度的方法 |
CN104439127A (zh) * | 2014-12-17 | 2015-03-25 | 重庆大学 | 计算品种钢连铸坯柱状晶区凝固组织真实生长方向的方法 |
CN105665684A (zh) * | 2016-04-13 | 2016-06-15 | 哈尔滨理工大学 | 一种铸件晶粒组织数值预测的方法 |
CN105868468A (zh) * | 2016-03-28 | 2016-08-17 | 南京航空航天大学 | 一种基于元胞自动机的新型邻居捕获方法 |
CN106407623A (zh) * | 2016-11-15 | 2017-02-15 | 南京航空航天大学 | 基于相场法的瞬态下焊接过程微观组织演变模拟方法 |
CN108469313A (zh) * | 2018-04-16 | 2018-08-31 | 大连理工大学 | 基于元胞自动机的结晶器铜板温度异常区域检测方法 |
CN105468875B (zh) * | 2015-12-28 | 2018-10-09 | 西北工业大学 | 确定tc17钛合金坯料在热处理炉内加热时间的方法 |
CN108694303A (zh) * | 2018-08-24 | 2018-10-23 | 上汽大众汽车有限公司 | 具有十六边形点焊热影响区结构的碰撞有限元模型建模方法 |
CN108959832A (zh) * | 2018-09-25 | 2018-12-07 | 武汉理工大学 | 用于优化M50NiL钢渗碳淬火工艺的晶粒长大预测模型 |
CN109446728A (zh) * | 2018-12-04 | 2019-03-08 | 燕山大学 | 近α钛合金低倍粗晶组织分布的预测方法 |
CN109522675A (zh) * | 2018-12-10 | 2019-03-26 | 桂林电子科技大学 | 锡基二元共晶合金微观组织的模拟及有限元求解分析方法 |
CN109885895A (zh) * | 2019-01-25 | 2019-06-14 | 南京航空航天大学 | 一种基于分子动力学的材料表面结冰形核过程的监测方法 |
CN110008631A (zh) * | 2019-04-18 | 2019-07-12 | 燕山大学 | 一种铜管轧制和冷却过程中参数模拟方法及其验证方法 |
CN110117710A (zh) * | 2019-05-27 | 2019-08-13 | 燕山大学 | 一种适合于连退机组rcs段的冷却速度预报方法 |
CN110276149A (zh) * | 2019-06-27 | 2019-09-24 | 大连理工大学 | 激光烧蚀金属靶材的径向热影响区宽度求解方法 |
CN110472322A (zh) * | 2019-07-31 | 2019-11-19 | 北京科技大学 | 一种基于热力学和动力学预测珠光体钢显微组织的方法 |
CN110489818A (zh) * | 2019-07-29 | 2019-11-22 | 西安理工大学 | 一种三元合金焊接熔池柱状枝晶生长数值模拟方法 |
CN111037143A (zh) * | 2019-11-28 | 2020-04-21 | 新兴际华集团有限公司 | 基于焊接热循环调控焊缝组织性能的方法 |
CN111988938A (zh) * | 2020-08-20 | 2020-11-24 | 广州创利宝科技有限公司 | 一种通讯设备存放保护设备 |
CN112182806A (zh) * | 2020-10-20 | 2021-01-05 | 同济大学 | 一种力流引导的介观结构设计方法 |
CN112214863A (zh) * | 2020-07-31 | 2021-01-12 | 中国科学院金属研究所 | 一种奥氏体不锈钢焊接热影响区粗晶区特征热循环曲线的获取方法 |
CN112992280A (zh) * | 2021-03-16 | 2021-06-18 | 苏州健雄职业技术学院 | 一种Ti-Al-V三元合金分子动力学α+β双相几何模型构建及势函数拟合方法 |
CN113130019A (zh) * | 2021-03-31 | 2021-07-16 | 西安理工大学 | 一种tc4钛合金焊接熔池固态相变的模拟方法 |
CN114283900A (zh) * | 2021-12-14 | 2022-04-05 | 燕山大学 | 针对近β钛合金低倍粗晶组织分布的预测与调控方法 |
CN115372135A (zh) * | 2022-08-10 | 2022-11-22 | 国标(北京)检验认证有限公司 | 一种测定高温合金晶粒转动角度的方法 |
CN117422205A (zh) * | 2023-12-18 | 2024-01-19 | 天津电力工程监理有限公司 | 基于数字孪生的装配式钢结构变电站建造管理系统和方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101329290A (zh) * | 2008-07-25 | 2008-12-24 | 沈阳工业大学 | 预测多组元合金凝固过程不同组元偏析分布的方法 |
-
2012
- 2012-07-17 CN CN201210247239.XA patent/CN102750425B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101329290A (zh) * | 2008-07-25 | 2008-12-24 | 沈阳工业大学 | 预测多组元合金凝固过程不同组元偏析分布的方法 |
Non-Patent Citations (4)
Title |
---|
X.H.ZHAN ET AL: "Simulation of grain morphologies and competitive growth in weld pool of Ni–Cr alloy", 《JOURNAL OF CRYSTAL GROWTH》 * |
YANHONG WEI ET AL: "Three-dimensional Monte Carlo simulation of discontinuous grain growth in HAZ of stainless steel during GTAW process", 《JOURNAL OF MATERIALS PROCESSING TECHNOLOGY》 * |
占小红: "Ni-Cr二元合金焊接熔池枝晶生长模拟", 《中国博士学位论文全文数据库(电子期刊)-工程科技I辑》 * |
马瑞: "镍基合金焊接熔池凝固过程微观组织模拟", 《中国博士学位论文全文数据库(电子期刊)-工程科技I辑》 * |
Cited By (43)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103246774A (zh) * | 2013-05-13 | 2013-08-14 | 天津大学 | 数值模拟p92钢管焊接热影响区宽度的方法 |
CN103246774B (zh) * | 2013-05-13 | 2016-03-16 | 天津大学 | 数值模拟p92钢管焊接热影响区宽度的方法 |
CN104439127A (zh) * | 2014-12-17 | 2015-03-25 | 重庆大学 | 计算品种钢连铸坯柱状晶区凝固组织真实生长方向的方法 |
CN104439127B (zh) * | 2014-12-17 | 2016-02-03 | 重庆大学 | 计算品种钢连铸坯柱状晶区凝固组织真实生长方向的方法 |
CN105468875B (zh) * | 2015-12-28 | 2018-10-09 | 西北工业大学 | 确定tc17钛合金坯料在热处理炉内加热时间的方法 |
CN105868468A (zh) * | 2016-03-28 | 2016-08-17 | 南京航空航天大学 | 一种基于元胞自动机的新型邻居捕获方法 |
CN105665684A (zh) * | 2016-04-13 | 2016-06-15 | 哈尔滨理工大学 | 一种铸件晶粒组织数值预测的方法 |
CN105665684B (zh) * | 2016-04-13 | 2017-11-10 | 哈尔滨理工大学 | 一种铸件晶粒组织数值预测的方法 |
CN106407623A (zh) * | 2016-11-15 | 2017-02-15 | 南京航空航天大学 | 基于相场法的瞬态下焊接过程微观组织演变模拟方法 |
CN106407623B (zh) * | 2016-11-15 | 2019-10-18 | 南京航空航天大学 | 基于相场法的瞬态下焊接过程微观组织演变模拟方法 |
CN108469313A (zh) * | 2018-04-16 | 2018-08-31 | 大连理工大学 | 基于元胞自动机的结晶器铜板温度异常区域检测方法 |
CN108469313B (zh) * | 2018-04-16 | 2019-06-14 | 大连理工大学 | 基于元胞自动机的结晶器铜板温度异常区域检测方法 |
CN108694303B (zh) * | 2018-08-24 | 2022-04-08 | 上汽大众汽车有限公司 | 有十六边形点焊热影响区结构的碰撞有限元模型建模方法 |
CN108694303A (zh) * | 2018-08-24 | 2018-10-23 | 上汽大众汽车有限公司 | 具有十六边形点焊热影响区结构的碰撞有限元模型建模方法 |
CN108959832A (zh) * | 2018-09-25 | 2018-12-07 | 武汉理工大学 | 用于优化M50NiL钢渗碳淬火工艺的晶粒长大预测模型 |
CN108959832B (zh) * | 2018-09-25 | 2022-12-02 | 武汉理工大学 | 用于优化M50NiL钢渗碳淬火工艺的晶粒长大预测方法 |
CN109446728A (zh) * | 2018-12-04 | 2019-03-08 | 燕山大学 | 近α钛合金低倍粗晶组织分布的预测方法 |
CN109522675A (zh) * | 2018-12-10 | 2019-03-26 | 桂林电子科技大学 | 锡基二元共晶合金微观组织的模拟及有限元求解分析方法 |
CN109885895A (zh) * | 2019-01-25 | 2019-06-14 | 南京航空航天大学 | 一种基于分子动力学的材料表面结冰形核过程的监测方法 |
CN109885895B (zh) * | 2019-01-25 | 2022-02-01 | 南京航空航天大学 | 一种基于分子动力学的材料表面结冰形核过程的监测方法 |
CN110008631A (zh) * | 2019-04-18 | 2019-07-12 | 燕山大学 | 一种铜管轧制和冷却过程中参数模拟方法及其验证方法 |
CN110008631B (zh) * | 2019-04-18 | 2020-11-24 | 燕山大学 | 一种铜管轧制和冷却过程中参数模拟方法及其验证方法 |
CN110117710A (zh) * | 2019-05-27 | 2019-08-13 | 燕山大学 | 一种适合于连退机组rcs段的冷却速度预报方法 |
CN110117710B (zh) * | 2019-05-27 | 2020-09-18 | 燕山大学 | 一种适合于连退机组rcs段的冷却速度预报方法 |
CN110276149A (zh) * | 2019-06-27 | 2019-09-24 | 大连理工大学 | 激光烧蚀金属靶材的径向热影响区宽度求解方法 |
CN110489818A (zh) * | 2019-07-29 | 2019-11-22 | 西安理工大学 | 一种三元合金焊接熔池柱状枝晶生长数值模拟方法 |
CN110489818B (zh) * | 2019-07-29 | 2022-10-21 | 西安理工大学 | 一种三元合金焊接熔池柱状枝晶生长数值模拟方法 |
CN110472322A (zh) * | 2019-07-31 | 2019-11-19 | 北京科技大学 | 一种基于热力学和动力学预测珠光体钢显微组织的方法 |
CN110472322B (zh) * | 2019-07-31 | 2021-07-23 | 北京科技大学 | 一种基于热力学和动力学预测珠光体钢显微组织的方法 |
CN111037143A (zh) * | 2019-11-28 | 2020-04-21 | 新兴际华集团有限公司 | 基于焊接热循环调控焊缝组织性能的方法 |
CN111037143B (zh) * | 2019-11-28 | 2022-09-06 | 新兴际华集团有限公司 | 基于焊接热循环调控焊缝组织性能的方法 |
CN112214863A (zh) * | 2020-07-31 | 2021-01-12 | 中国科学院金属研究所 | 一种奥氏体不锈钢焊接热影响区粗晶区特征热循环曲线的获取方法 |
CN112214863B (zh) * | 2020-07-31 | 2024-01-02 | 中国科学院金属研究所 | 一种奥氏体不锈钢焊接热影响区粗晶区特征热循环曲线的获取方法 |
CN111988938A (zh) * | 2020-08-20 | 2020-11-24 | 广州创利宝科技有限公司 | 一种通讯设备存放保护设备 |
CN112182806A (zh) * | 2020-10-20 | 2021-01-05 | 同济大学 | 一种力流引导的介观结构设计方法 |
CN112992280A (zh) * | 2021-03-16 | 2021-06-18 | 苏州健雄职业技术学院 | 一种Ti-Al-V三元合金分子动力学α+β双相几何模型构建及势函数拟合方法 |
CN113130019B (zh) * | 2021-03-31 | 2024-02-06 | 西安理工大学 | 一种tc4钛合金焊接熔池固态相变的模拟方法 |
CN113130019A (zh) * | 2021-03-31 | 2021-07-16 | 西安理工大学 | 一种tc4钛合金焊接熔池固态相变的模拟方法 |
CN114283900A (zh) * | 2021-12-14 | 2022-04-05 | 燕山大学 | 针对近β钛合金低倍粗晶组织分布的预测与调控方法 |
CN114283900B (zh) * | 2021-12-14 | 2024-09-13 | 燕山大学 | 针对近β钛合金低倍粗晶组织分布的预测与调控方法 |
CN115372135A (zh) * | 2022-08-10 | 2022-11-22 | 国标(北京)检验认证有限公司 | 一种测定高温合金晶粒转动角度的方法 |
CN117422205B (zh) * | 2023-12-18 | 2024-04-19 | 天津电力工程监理有限公司 | 基于数字孪生的装配式钢结构变电站建造管理系统和方法 |
CN117422205A (zh) * | 2023-12-18 | 2024-01-19 | 天津电力工程监理有限公司 | 基于数字孪生的装配式钢结构变电站建造管理系统和方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102750425B (zh) | 2014-09-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102750425B (zh) | 一种焊接过程热影响区组织演变的模拟方法 | |
Bozzolo et al. | Viewpoint on the formation and evolution of annealing twins during thermomechanical processing of FCC metals and alloys | |
Liu et al. | 2D cellular automaton simulation of hot deformation behavior in a Ni-based superalloy under varying thermal-mechanical conditions | |
Jaafar et al. | A review of dendritic growth during solidification: Mathematical modeling and numerical simulations | |
Li et al. | Microstructure simulation on recrystallization of an as-cast nickel based single crystal superalloy | |
Yang et al. | Three dimensional Monte Carlo simulation of grain growth during GTA welding of titanium | |
Boettinger et al. | Solidification microstructures: recent developments, future directions | |
Hoyt et al. | Crystal–melt interfaces and solidification morphologies in metals and alloys | |
Zhang et al. | Hot compression deformation and constitutive modeling of GH4698 alloy | |
Song et al. | Numerical simulation of β to α phase transformation in heat affected zone during welding of TA15 alloy | |
Tian et al. | Experimental and simulation on the grain growth of P/M nickel-base superalloy during the heat treatment process | |
Fang et al. | Modeling the competition between solid solution and precipitate strengthening of alloys in a 3D space | |
Dong et al. | Current development in quantitative phase-field modeling of solidification | |
Wedberg et al. | Modelling flow stress of AISI 316L at high strain rates | |
CN112185488B (zh) | 基于数据驱动多组元镍基高温合金γ`相演化的预测方法 | |
Chandra et al. | Deformation behavior of Nickel-based superalloy Su-263: Experimental characterization and crystal plasticity finite element modeling | |
Han et al. | Temperature-dependent fatigue response of a Fe44Mn36Co10Cr10 high entropy alloy: A coupled in-situ electron microscopy study and crystal plasticity simulation | |
Rafique | Modeling and Simulation of Heat Transfer Phenomena | |
Li et al. | Finite element simulation and process optimization of microstructure evolution in the formation of Inconel 718 alloy bolts | |
Qian et al. | Fatigue behaviour analysis of thermal cyclic loading for through-silicon via structures based on backstress stored energy density | |
Sun et al. | Adapted solute drag model for impurity-controlled grain boundary motion | |
Zhang et al. | Numerical simulation of microstructure evolution during directional solidification process in directional solidified (DS) turbine blades | |
Hemmer et al. | A process model for the heat-affected zone microstructure evolution in duplex stainless steel weldments: Part II. Application to electron beam welding | |
Miller | Heat Extraction and Dendritic Growth during Directional Solidification of Single-Crystal Nickel-base Superalloys. | |
Yan et al. | Numerical simulation of dendrite growth in nickel-based superalloy and validated by in-situ observation using high temperature confocal laser scanning microscopy |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140903 |