宏微观耦合模拟熔池不同区域中枝晶竞争生长过程

针对熔化焊过程建立了宏微观耦合模型,模拟了熔池内不同区域凝固过程中随机取向枝晶的竞争生长过程. 通过宏观三维有限元模型计算熔池中瞬态的传热传质过程,利用双线性插值算法将凝固参数传递给微观组织模型. 采用元胞自动机法模拟随机取向的枝晶在熔池凝固条件下的竞争生长过程. 模拟结果表明,所建立的微观模型能够精确模拟任意生长取向的枝晶. 凝固条件中最大温度梯度方向对枝晶竞争过程有明显选择作用,生长方向与最大温度梯度方向相同或接近的枝晶在竞争中具有更大优势. 焊缝中的晶粒组织由枝晶簇发展形成,晶粒组织的形貌演变取决于相邻枝晶簇之间的竞争过程,具有择优取向的枝晶簇会逐渐排挤非择优取向的枝晶簇并最终将其阻挡在凝固组织内部,宏观晶粒的取向与其内部枝晶簇的生长方向并不一定相同. 熔池中心线附近区域在焊接过程中具有更小的温度梯度、更大的凝固速率以及更大的局部冷却速率,凝固后可以获得更加细小的焊缝枝晶组织. 枝晶间距的模拟结果与相应凝固条件下的试验数据符合较好. 关键词： 焊接熔池 枝晶形貌 竞争生长 元胞自动机     

液氮冷却条件下激光快速熔凝Ni-28 wt%Sn合金组织演变

研究了在液氮冷却条件下激光快速熔凝Ni-28 wt%Sn亚共晶合金的组织演化过程. 结果显示, 熔池从上至下可以分为三个区域: 表层为平行激光扫描方向的α-Ni转向枝晶区; 中部为近乎垂直于熔池底部外延生长的α-Ni柱状晶区; 底部为少量的残留α-Ni初生相和大量的枝晶间(α-Ni+Ni₃Sn) 共晶组织. 激光熔凝区组织受原始基材组织的影响很大, 熔池中的α-Ni枝晶生长方向受到了热流方向和枝晶择优取向的双重影响. 与基材中存在的层片状、棒状和少量离异(α-Ni+Ni₃Sn)共晶的混合组织相比, 熔池内的共晶组织皆为细小的规则(α-Ni+Ni₃Sn)层片状共晶, 皆垂直于熔池底部外延生长, 并且从熔池顶部至底部, 共晶层片间距逐渐增大. 分别应用描述快速枝晶生长的Kurz-Giovanola-Trivedi 模型和描述快速共晶生长的Trivedi-Magnin-Kurz模型对熔池表层凝固界面前沿的过冷度进行估算, 发现熔池表层α-Ni 枝晶和(α-Ni+Ni₃Sn)层片共晶的生长过冷度在50.4-112.5 K 之间, 远大于相应深过冷凝固(α-Ni+Ni₃Sn) 反常共晶生长的临界过冷度20 K, 这说明文献报道的临界过冷度并不是反常共晶出现的充分条件.     

基于MeshTV界面重构算法的二元合金自由枝晶生长元胞自动机模型

本文基于MeshTV界面重构算法, 发展了二元合金凝固自由枝晶生长的元胞自动机 (cellular automaton, CA) 模型. 通过采用MeshTV界面重构算法, 在细化的界面元胞内重构出了固液界面的位置. 在此基础上, 发展了一种同时适合描述纯物质与合金凝固界面生长的动力学模型. 与非界面重构的CA模型相比, 本文所发展的模型可以在较大的网格尺寸下实现模型的收敛, 同时网格各向异性不明显, 且能够反映界面能各向异性参数ε 对自由枝晶生长的影响. 在ε =0.02时, 通过与描述自由枝晶生长的LGK理论模型相比较, 发现计算的枝晶尖端速度与LGK理论模型的预测符合较好, 而计算的枝晶尖端半径比LGK理论预测值大于约20%.     

纯物质枝晶凝固的元胞自动机模型

采用元胞自动机模型模拟了纯物质凝固微观组织演化行为. 与以往采用KGT模型简化考虑不同,采用界面能量守恒方程来自洽地确定液固界面的生长动力学,考察了模型的网格尺寸相关性和时间相关性,并对模型的稳定性进行了验证. 采用这个模型模拟了纯丁二腈在过冷熔体中的自由枝晶生长,并与Glicksman的等温枝晶生长实验(isothermal dendritic growth experiment,IDGE)和LM-K临界稳定性原理进行了对比. 模拟的形貌再现了实验形貌的特征,并且特征参数与理论结果及实验结果符合得很好. 关键词： 元胞自动机 界面 凝固 临界稳定性原理     

界面能各向异性对定向凝固枝晶生长的影响

采用元胞自动机 (cellular automaton, CA) 模型研究了界面能各向异性对二维定向凝固枝晶生长的影响. 模拟结果显示当晶体的择优生长方向与热流方向一致时, 随着界面能各向异性强度的增大, 凝固组织形态由弱界面能各向异性时的海藻晶转变为强界面能各向异性时的树枝晶. 同时, 界面能各向异性强度会影响稳态枝晶尖端状态的选择, 界面能各向异性越强, 定向凝固稳态枝晶尖端半径越小, 尖端界面前沿的液相浓度和过冷度越小. 稳态枝晶生长的尖端状态选择参数与界面能各向异性强度也存在标度律的指数关系, 而枝晶一次间距则受界面能各向异性强度影响较弱. 当晶体的择优生长方向与热流方向呈-40°夹角时随着界面能各向异性强度的增大, 凝固组织形态由海藻晶逐渐转变为退化枝晶, 后又逐渐转变为倾斜枝晶. 关键词： 元胞自动机 枝晶 界面能各向异性     

横向限制下凝固微观组织演化的相场法模拟

凝固过程中横向限制挡板的存在对晶体微观结构的演化存在重要的影响, 不同性质的横向挡板将产生不同的限制效应, 对最终凝固微观组织形成起着决定性作用. 本文利用非等温相场模型, 定性地模拟了纯金属Ni凝固过程中横向限制的存在对其枝晶微观形貌演化的影响, 研究了不同尺寸及性质的横向挡板对枝晶微观结构形成的影响, 讨论了横向限制对不同初始枝晶间距枝晶形貌发展的作用. 计算结果表明, 横向限制挡板的存在将直接影响凝固过程中微观组织的形貌演化过程并最终改变微观结构. 随着横向挡板间距的减小, 微观组织变化更加明显; 挡板初始温度越低, 枝晶形貌改变越明显; 初始枝晶间距越大, 形貌变化越明显; 不同挡板高度对微观结构具有基本相同的影响. 关键词： 相场模拟 微观组织演化 横向限制     

基于元胞自动机方法的定向凝固枝晶竞争生长数值模拟

通过耦合温度场模型、溶质扩散方程以及枝晶生长动力学方程等重要因素,建立了一种改进的元胞自动机模型.该模型通过采用偏心算法消除网格各向异性,实现了二维尺度上任意角度枝晶生长的模拟,同时适用于模拟三维尺度上枝晶的生长过程.利用建立的模型开展了定向凝固枝晶竞争生长过程的数值模拟.为了体现本模型的有效性,模拟了透明合金的竞争生长过程,并与实验结果符合良好.镍基高温合金汇聚竞争和发散竞争的模拟结果清楚地展现了不同抽拉速度和枝晶优先生长角度下枝晶的竞争生长过程,并且模拟结果与理论模型相符合.三维枝晶生长的模拟结果表明本模型可以用来模拟三维枝晶一次臂间距的调整过程.     

晶体取向对定向凝固枝晶生长的影响

采用类金属透明模型合金丁二腈-1.0 wt%乙醇(SCN-1.0 wt% Eth)合金, 考察了晶体取向对定向凝固过程中晶粒的平界面失稳孕育时间、枝晶形态演化以及枝晶一次间距的影响. 结果表明, 随着枝晶择优生长方向与温度梯度方向夹角的增大, 晶粒的平界面失稳孕育时间增加, 界面的稳定性增强; 对于不同晶体取向的枝晶形态演化, 枝晶择优生长方向与温度梯度方向夹角越大, 枝晶二次臂不对称生长越严重, 同时, 具有生长优势的枝晶二次臂对相邻枝晶的生长的抑制越强烈; 至于不同晶体取向的枝晶一次间距, 随着枝晶择优生长方向与温度梯度方向夹角的增大, 枝晶一次间距增大. 关键词： 定向凝固 平界面失稳 枝晶间距 晶体取向     

相场方法模拟铝合金三维枝晶生长

以相场模型为基础，采用宏微观耦合方法和界面捕获液态方法对铝合金枝晶生长进行模拟计算.为解决试样全场微观计算的困难，采取宏微观耦合的计算方法，试样整体计算温度场，而微观组织计算只在一个确定的宏观单元内进行，宏观微观计算交替耦合进行.在不改变相场模型的条件下，提出界面捕获液态计算方法.赋值计算单元界面标志，只对界面处的单元求解相场变量，当枝晶生长时，捕获液态单元为界面从而推进界面，并对捕获到的单元校正相场变量.通过界面捕获液态方法加速相场模型的计算，实现了铝合金试件局部三维单晶粒和多晶粒的模拟.对模拟结果与实 关键词： 相场方法 液态捕获 微观组织 枝晶     

晶体取向对定向凝固平界面失稳行为的影响

采用丁二腈-丙酮透明模型合金研究了不同晶体取向的晶粒在定向凝固条件下的平界面失稳过程.实验选择了三个界面失稳后具有不同生长形态的典型晶粒作为研究对象, 分别为择优生长枝晶、倾斜枝晶和海藻晶.结果表明可发展为择优生长枝晶的晶粒的平界面失稳孕育时间和初始扰动波长最小,海藻晶次之, 倾斜枝晶最大,这与以往的解析结果和相场模拟结果一致. 同时,实验观察发现可发展为择优生长枝晶和倾斜枝晶的晶粒的界面非稳态演化过程与海藻晶显著不同,这表明平界面失稳的非稳态演化过程与晶体取向相关.     

Microstructure characteristics of laser-MIG hybrid welded mild steel

To deepen the understanding of laser-arc hybrid welding, the weld shape and microstructure characteristics of laser-metal inert gas hybrid welded mild steel were analyzed. The results showed typical hybrid weld could be classified as two parts: the wide upper zone and the narrow nether zone, which were defined as arc zone and laser zone, respectively. In the hybrid weld, the microstructure, alloy element distribution and microhardness all have evident difference between laser zone and arc zone. The microstructure of arc zone consists of coarse columnar dendrite and fine acicular dendrite between the columnar dendrites, but that of laser zone is composed of fine equiaxed dendrite in weld center and columnar dendrite around the equiaxed dendrite. Compared to arc zone, laser zone has finer grain size, higher microhardness, smaller alloy element content in the fusion zone and narrower heat affected zone. The discussions demonstrated that the observed difference was caused by the difference of temperature gradient, crystallizing and the effects of arc pressure on the molten pool between laser zone and arc zone.     

Numerical simulation of temperature distribution and TiC growth kinetics for high power laser clad TiC/NiCrBSiC composite coatings

A three dimensional model was proposed to simulate high power laser clad TiC/NiCrBSiC composite coatings on Ti6Al4V alloys. The temperature distribution, temperature curves on different nodes, three dimensional shape and size of TiC melting region, molten pool and heat affected zone (HAZ) of the substrate were obtained. To have a clear physical insight into the phase transformation and microstructure evolution in the coatings during laser cladding process, a theoretical kinetic analysis was performed to elucidate the nucleation, growth velocity, and size of TiC particles on the basis of simulated temperature curves of the molten pool. A good quality TiC/NiCrBSiC composite coating with low dilution rate and excellent metallurgical bond was fabricated under optimal processing parameters using powder mixture of TiC and NiCrBSiC as clad material and cuboid of Ti6Al4V alloys as substrate. To validate the reliability of the proposed model, the theoretical results were compared with the microstructure of the coatings. It shows that these theoretical results are in excellent agreement with the experiment cases.     

Heat transfer and diffusion-controlled kinetics of liquid–solid phase in titanium matrix composite during selective laser melting

This study describes a self-consistent theoretical model of simulating diffusion-controlled kinetics on the liquid–solid phase boundary during high-speed solidification in the melt pool after the selective laser melting (SLM) process for titanium matrix composite based on Ti–TiC system. The model includes the heat transfer equation to estimate the temperature distribution in the melt pool and during crystallization process for some deposited layers. The temperature field is used in a micro region next to solid–liquid boundary, where solute micro segregation and dendrite growth are calculated by special approach based on transient liquid phase bonding. The effect of the SLM process parameters (laser power, scanning velocity, layer thickness and substrate size) on the microstructure solidification is being discussed.     

Finite element analysis of laser irradiated metal heating and melting processes

Laser technology has shown fast growth due to its demands in material processing and manufacturing. Laser material processing includes various applications like cutting, welding, drilling, cladding and surface treatment. In laser surface treatment, the material properties at the surface are altered through surface alloying and transformation hardening. In this study, an enthalpy-based computational model is developed for analyzing laser heating and melting of metals. The solution to the problem is obtained by using a finite element method and validated by comparing the results with that given by an analytical solution to a limiting case problem. A solution algorithm and a computational code are developed to estimate the temperature distribution, solid-liquid interface location and shape and size of the molten pool. The computational model is validated by comparing results with a limiting case analytical model. The study is conducted to analyze the heating rate, the heat affected zone, and the shape and size of the molten pool using a Gaussian laser beam.     

Elimination of turned dendrite in laser multilayer deposition of Rene88DT superalloy on DD3 single crystal

Laser multilayer deposition of Rene88DT superalloy on DD3 single-crystal substrate is conducted.The influences of the crystal orientation of the substrate and the profile of the solid/liquid interface of the molten pool on the deposited microstructure are investigated.A unique strategy is proposed by adjusting the angle between the substrate surface and the substrate crystal orientation.This approach prevents the formation of the turned dendrite,thus obtaining a fully directional microstructure in the deposits.     

Absolute stability of the solidification interface in a laser resolidified Zn--2wt.%Cu hypoperitectic alloy

This paper reports on laser surface remelting experiments performed on a Zn--2wt.%}Cu hypoperitectic alloy by employing a 5kW CW CO₂ laser at scanning velocities between 6 and 1207mm/s. The growth velocities of the microstructures in the laser molten pool were accurately measured. The planar interface structure caused by the high velocity absolute stability was achieved at a growth velocity of 210~mm/s. An implicit expression of the critical solidification velocity for the cellular--planar transition was carried out by nonlinear stability analyses of the planar interface. The results showed a better agreement with the measured critical velocity than that predicted by M--S theory. Cell-free structures were observed throughout the whole molten pool at a scanning velocity of 652~mm/s and the calculated minimum temperature gradient in this molten pool was very close to the critical temperature gradient for high gradient absolute stability (HGAS) of the \eta phase. This indicates that HGAS was successfully achieved in the present experiments.     

镍基合金激光熔覆数值模拟及凝固参数预测

本文采用有限体积法对激光熔覆过程的温度场分布和熔体流动进行了数值模拟.基于CALPHAD相图法计算了基体和粉末的热物理性质,采用三维热源精确预测了凝固过程和温度分布,研究了Marangoni对流对熔池尺寸的影响.在熔池凝固过程中模拟所得出的温度梯度和凝固速度,预测了熔覆层凝固组织的演变趋势,相应的显微组织与实验结果吻合...     

Numerical simulation of dendrite growth and microsegregation formation of binary alloys during solidification process

The dendrite growth and solute microsegregation of Fe-C binary alloy are simulated during solidification process by using cellular automaton method. In the model the solid fraction is deduced from the relationship among the temperature, solute concentration and curvature of the solid/liquid interface unit, which can be expressed as a quadric equation, instead of assuming the interface position and calculating the solid fraction from the interface velocity. Then by using this model a dendrite with 0 and 45 degree of preferential growth direction are simulated respectively. Furthermore, a solidification microstructure and solute microsegregation are simulated by this method. Finally, different Gibbs-Thomson coefficient and liquid solute diffusing coefficient are adopted to investigate their influences on the morphology of dendrite.     

Solute transport and composition profile during direct metal deposition with coaxial powder injection

Direct metal deposition (DMD) with coaxial powder injection allows fabrication of three-dimensional geometry with rapidly solidified microstructure. During DMD, addition of powder leads to the interaction between laser and powder, and also the redistribution of solute. The concentration distribution of the alloying element is very important for mechanical properties of the deposited clad material. The evolution of concentration distribution of carbon and chromium in the molten pool is simulated using a self-consistent three-dimensional model, based on the solution of the equations of mass, momentum, energy conservation and solute transport in the molten pool. The experimental and calculated molten pool geometry is compared for model validation purposes.