Background technology
Component segregation is intrinsic characteristic during aluminium alloy solidification.When dissolubility in solid phase for the solute element is less than it in liquid
During dissolubility in phase, solid liquid interface will exclude solute in liquid phase;And work as dissolubility in solid phase for the solute element and be more than
During dissolubility in liquid phase, solute will lead to solid liquid interface forward position the poor area of solute from liquid phase to solid-state diffusion.This by
In alloy graining process, the reallocation of solute leads to the uneven components after alloy graining, as component segregation.Generally, solute
Element dissolubility in the liquid phase is higher than in solid phase, and therefore solute has disengaging solidifying phase to be pulled to inclining in the liquid of forward position
To the dendritic arm that the liquid rich in solute for the end solidification part is grown is split, forming component segregation zones.Component segregation is foundry goods
One of major defect, reduce alloy solid solution strengthening effect, meanwhile, the low melting point precipitated phase that interdendritic is formed, in the weldering of foundry goods
Easily become formation of crack in termination process, reduce mechanical property and the welding performance of foundry goods.
According to classical solidification theory, under the conditions of general alloy graining, the solute redistribution behavior of micro-scale
And the microsegregation producing depends not only on the physical parameter of alloy (as solid phase diffusion welding ds, Liquid Diffusion Coefficient dl, molten
Matter partition coefficient k) and external process condition such as freezing rate v (or local solidification time tf), thermograde g, but also with
The factors such as the growth pattern of dendrite, dendrite morphology are closely related.
Simplest analytical model is lever analytical model and scheil analytical model under local equilibrium's curing condition,
Other is all correction type analytical model.Based on the consideration of reverse diffusion, lever model and scheil model are two kinds of limiting cases.
Lever model is a kind of equilibrium freezing model it is assumed that solute all fully spreads solid-liquid is biphase, that is,
cl=c0[1-(1-k)fs]-1①
Formula 1. in, clFor liquid phase solute concentration, c0For initial solute concentration, fsFor fraction solid, k is solute distribution coefficient.
Scheil model is also referred to as non-equilibrium lever model, wherein assumes that solute fully spreads in the liquid phase and solid phase no expands
Dissipate, in model, do not account for that forming core is supercool, the formation of dendritic arm roughening, macroscopic material transmission and pore, that is,
cl=c0(1-fs)(k-1)②
Obviously, as fraction solid fsC when=1lTrend towards infinity, therefore scheil model can not be dense to final solute
Degree is accurately calculated.
H.d.brody&m.c, flemings consider the back-diffusion of solute element in process of setting, that is, have
cl=c0[1-(1-2αk)fs](k-1)/(1-2αk)③
Formula 3. in, α be the nondimensional solutes accumulation factor, or claim nondimensional diffusion time (fourier number), it represent
Diffusion boundary layer thickness δ in solid phasesWith the proportionate relationship of system dimension l, it is represented by:
When α=0 is no to spread in solid phase, 3. formula is reduced to scheil model;
When α=0.5,3. formula is reduced to lever model.
W.kurz et al. is corrected to α value, proposes to replace α with α ', that is,
Although scholars have studied the coagulating property of alloying component change of multicomponent alloy, in process of setting energetically in recent years
Solid-back diffusion, the limited diffusion of liquid, dendritic arm roughening, dentrite tip is supercool, eutectic interface undercooling and coagulated volume change
Affecting it is established that solute concentration parses with the microsegregation that fraction solid increases in some process of setting Deng solute assigning process
Model and numerical model.But the microsegregation model of existing solute distribution be built upon during research solidification solute concentration with
The change that fraction solid increases, can only represent interfacial concentrationWith fraction solid fsFunctional relationship, can not determine therefrom that
The component distributing of whole solid phase, and must take into this dynamic parameter of fraction solid, and not can determine that large complicated for one
Constitutional detail different structure position, by the degree of segregation under many factors collective effect it is impossible to instruct large and complex structure foundry goods
Foundry Production.
Content of the invention
For defect of the prior art, it is an object of the invention to provide composition in a kind of Complicated structure casting casting process
The Forecasting Methodology of degree of segregation, this Forecasting Methodology can predict parts with complex structures different structure position in different curing conditions
Microsegregation degree, has directive significance to the Foundry Production of large and complex structure foundry goods.
For realizing object above, the present invention provides a kind of prediction of component segregation degree in Complicated structure casting casting process
Method, the method includes:
(1) obtain calculating thermal physical property parameter
Calculate the thermal physical property parameter obtaining alloy, ginseng using the thermophysical parameter computing module of commercialization jmatpro software
Number includes solute balance partition coefficient k, Liquid Diffusion Coefficient dl, solid phase diffusion welding ds, alloy initial solidification temperature tl, eutectic
Reaction temperature ts;
(2) obtain the temperature field data at each position of foundry goods
Application procast finite element software carries out analogue simulation to the casting filling of Complicated structure casting, process of setting, obtains
Take the temperature field data at each position of foundry goods, data includes setting rate v of molten metal, local setting time tf;Or local
Setting time tf, the method that also can be combined with sunykatuib analyses by precision-investment casting of setting rate v obtained.
(3) obtain secondary dendrite arm spacing λ at each structure position of foundry goods2
Calculate the secondary dendrite of each structure position of foundry goods using the secondary dendrite arm spacing computing module that procast software carries
Spacing λ2;Or, secondary dendrite arm spacing λ2Also can be obtained by the method that precision-investment casting is combined with sunykatuib analyses.
(4) calculate segregation exponent m si value at each structure position of foundry goods
The Solidification Parameters obtaining from procast: setting rate v of molten metal, local setting time tf, secondary dendrite arm spacing
λ2, and the thermophysical parameter that jmatpro calculates: solute balance partition coefficient k, Liquid Diffusion Coefficient dl, solid phase diffusion welding
Ds, substitutes in microsegregation deciding degree exponent m si formula, and the viewcast module by procast, calculates labyrinth
Segregation exponent m si of each structure position of foundry goods, the size of the segregation exponent m si value according to each structure position of Complicated structure casting is sentenced
Determine the microsegregation degree of each structure position of foundry goods: segregation exponent m si is bigger, and microsegregation is more serious, solidification latter stage is formed
It is higher that interdendritic separates out phase content;Wherein:
In formula: pelFor dimensionless solute p é let number, characterize in the liquid phase of solid liquid interface forward position by solute degree of supersaturation institute really
Fixed driving force;fosFor dimensionless Fourier number, also referred to as the solid-back diffusion factor;K is solute balance partition coefficient;dlFor
Liquid Diffusion Coefficient;dsFor solid phase diffusion welding;V is setting rate;L is characterized diffusion length;tfFor local solidification time;For comprising the physical parameter of alloy material, depending on the selection of material, for the different structure on Same Part
Position, a may be considered constant.
Preferably, in step (2), secondary dendrite arm spacing λ2Computing formula be
λ2=(mtf)1/3
Wherein:
In formula, γ is gibbs thomson coefficient, c0For initial solute concentration, ceutThere is eutectic reaction for solidification latter stage
When liquid phase solute concentration, m is alloy liquid phase line slope, and m is with the related parameter of material.
It is highly preferred that m value calculates with each parameter from Thermodynamic Calculation Software jmatpro or list of references
l.nastac.numerical modeling of solidification morphologies and segregation
Patterns in cast dendritic alloys [j] .acta mater.vol.17.pp:4253-4262.1999. obtains
?.
Preferably, calculate in panel in procast, corresponding initial solidification temperature t is inputted according to different alloysl, eutectic
Reaction temperature tsAnd the m value calculating, secondary dendrite arm spacing value λ of each structure position of foundry goods can be calculated2.
Compared with prior art, the present invention has a following beneficial effect:
The present invention, by finite element software procast and Thermodynamic Calculation Software jmatpro, can effectively predict labyrinth
The microsegregation degree of the solute element of each structure position of foundry goods;By the application of the inventive method, not only can be without dissection
Foundry goods can predict the degree of segregation of a certain solute element of each structure position of foundry goods, for optimizing and revising at casting technique and heat
Science and engineering skill provides foundation, and the method can be to the prediction of other solute element degree of segregation of multicomponent alloy, range of application
Extensively.
Specific embodiment
With reference to specific embodiment, the present invention is described in detail.Following examples will be helpful to the technology of this area
Personnel further understand the present invention, but the invention is not limited in any way.It should be pointed out that the ordinary skill to this area
For personnel, without departing from the inventive concept of the premise, some deformation can also be made and improve.These broadly fall into the present invention
Protection domain.
In the present embodiment, the determination of microsegregation deciding degree exponent m si (micro-segregation index) is based on
Considered below.
Theoretical by means of buckingham π non-dimension analysis widely used in mathematical physics, to solute in process of setting
Distribution carries out non-dimension analysis.Diffusion flux research based on solid liquid interface both sides liquid phase and solid phase solute atoms, the present invention is comprehensive
Close the back-diffusion considering the limited diffusion of solid liquid interface forward position liquid phase solute and solid phase solute to alloy graining solute redistribution
Behavior is studied.The employing p é let number of solid liquid interface forward position solutes accumulation represents, p é let number is dendrite end radius of curvature r
With liquid phase solute Equivalent boundary layer thickness δcRatio, liquid phase solute Equivalent boundary layer thickness δc(δc=2dl/ v) represent diffusion length
In degree, contained whole solutes are equal to the solute total amount contained by infinite boundary layer.Due to the solute concentration gradient of dendrite end,
The solute concentration gradient in the spheroid forward position being to grow up can approx be look at, according to the asking of diffusion field around the globular crystal growing
Solution method, if only considering radial diffusion, the boundary layer thickness around globular crystal growing is equal to the radius of ball, i.e. dendrite end
Solutes accumulation length l in portion is approximately equal to the radius r of ball.Therefore, p é let number is represented by
P é let number characterizes in the liquid phase of solid liquid interface forward position driving force determined by solute degree of supersaturation.
Due to eutectic usually occurring in actual process of setting and liquid concentration can not possibly continue to increase, in order to truly
The component distributing of description solidification end is it is necessary to consider the back-diffusion of solute atoms in solid phase.According to the principle of mass conservation, in solid phase
The speed degree of solute atoms back-diffusion depends on dimensionless group α, also referred to as fos, it is anti-that it characterizes solute atoms in solid phase
The ability of diffusion.
According to Solute mass conservation, characterize the segregation exponent m si (micro- of the solute segregation degree in process of setting
Segregation index) it is represented by
WhereinSubstitute into (9)
WhereinIt is the physical parameter comprising alloy material, depending on the selection of material.For Same Part
On different structure position, a may be considered constant.Knowable to formula (10), judgement microsegregation degree proposed by the present invention
Segregation exponent m si value size, not only with the thermal physical property parameter solute balance partition coefficient k of alloy, Liquid Diffusion Coefficient dl, solid phase
Diffusion coefficient dsCorrelation, and also local solidification condition setting rate v of alloy, feature diffusion length l and local solidification when
Between tfRelated.
For general casting process, dendrite is grown with isometry branch crystal type.Now, secondary dendrite arm spacing λ2It is solute
The characteristic length of diffusion.Therefore, the half taking secondary dendrite arm spacing is diffusion length, that is,
Substitution formula (10)
It is below the present invention one application the present embodiment:
The present embodiment provides a kind of Forecasting Methodology of component segregation degree in Complicated structure casting casting process, to being cast as
The k4169 nickel base superalloy large-scale complex thin-wall foundry goods of type carries out the prediction of microsegregation degree.The design of this foundry goods is not containing
The large-area thin plate of stack pile, variable cross-section, the structure such as cylinder.Casting Three-dimensional illustraton of model is as shown in figure 1, the foundry goods of cast acquisition is real
Thing figure is as shown in Figure 2.Because the alloying element that k4169 nickel base superalloy adds is many, alloying level is high, solidification temperature range
Width, the solidifying segregation of alloying element is serious.Microsegregation is one of major defect of k4169 nickel base superalloy foundry goods, to casting
The mechanical property of part, welding performance and processing characteristics have important impact.K4169 nickel base superalloy casting solidification process
Middle nb solute atoms dissolubility in the base is little, and in process of setting, solid liquid interface will exclude a large amount of nb solute atoms to liquid phase,
When being cooled to 1180 DEG C, there is γ+laves eutectic reaction, form laves phase in interdendritic.If due to alloy graining process
Component segregation, the laves phase amount that solidification latter stage is formed is more, regardless of their form and distribution, because solid in laves phase
The intensified elements such as a large amount of nb are determined, thus reducing solid solution strengthening effect.Meanwhile, decrease hardening constituent γ ' or γ " quantity,
Reduce precipitation strengthening effect, cause alloy mechanical property to deteriorate.Therefore, solute element nb segregation and laves phase are eliminated always
It is the key technology difficulty of k4169 nickel base superalloy castings production.
(1) obtain calculating thermal physical property parameter
By jmatpro software, PHASE DIAGRAM CALCULATION is carried out to k4169 alloy, and combine list of references
l.nastac.numerical modeling of solidification morphologies and segregation
Patterns in cast dendritic alloys [j] .acta mater.vol.17.pp:4253-4262.1999. obtains
Each thermophysical parameter used by calculating, each parameter is as follows:
γ=3.65 × 10-7(km), knb=0.48, dl,nb=3 × 10-9(m2s-1), ds,nb=2.82 × 10-13(m2s-1),
c0,nb=4.316 (wt.pct), ceut,nb=19.1 (wt.pct), ml=-10.9, tl=1360 DEG C, ts=1180 DEG C;
(2) obtain the temperature field data at each position of foundry goods
Application procast finite element software carries out emulating mould to the casting filling of foundry goods as shown in Figure 1, process of setting
Intend, obtain the temperature field data at each position of foundry goods, such as local setting time t of molten metalf, setting rate v;
Or local setting time tf, the side that also can be combined with sunykatuib analyses by precision-investment casting of setting rate v
Method obtains.For fixing k4169 alloy, based on fixing cast model and full form casting process condition, carry out actual fusible pattern
Casting obtains foundry goods.In casting process, it is pre-placed high temperature resistance tungsten-rhenium thermal couple in casting features structure position, in real time measurement casting
Make the cooling curve of process molten metal.On the one hand, carry out PHASE DIAGRAM CALCULATION acquisition in conjunction with jmatpro software to k4169 alloy to open
Beginning solidification temperature tlAnd solidification end temp ts, you can obtain local setting time tfWith rate of cooling r, local setting time tf
It is molten metal to start to solidify to solidification end the experienced time, and rate of cooling r passes through formula r=- (ts-tl)/tfObtain
?.On the other hand, in conjunction with the position of the branched thermocouple of equidistant layout, you can obtain thermograde g of appointed part, then root
According to formula r=-gv, you can obtain setting rate v of appointed part molten metal.Finally by the actual measured value of thermocouple and simulation
Value is compared, if both identical then it is assumed that this simulation obtain local setting time tf, setting rate v be foundry goods office
Domain setting time tf, setting rate v, if instead then two; would differ, adjust model boundary condition setting, then be simulated, directly
Identical to both.
(3) obtain secondary dendrite arm spacing λ at each structure position of foundry goods2
Calculate the λ of each structure position of foundry goods using the secondary dendrite arm spacing computing module that procast software carries2Value.Warp
The m value calculating acquisition k4169 alloy is 19.43 (um3s‐1);Secondary dendrite arm spacing λ2Computing formula be
λ2=(mtf)1/3(12)
Wherein:
In formula (13), γ is gibbs thomson coefficient, c0For initial solute concentration, ceutThere is eutectic for solidification latter stage
Liquid phase solute concentration during reaction, m is alloy liquid phase line slope, and m is with the related parameter of material.Calculate in panel in procast,
Corresponding initial solidification temperature t is inputted according to different alloysl, eutectic reaction temperature ts, and the m value being calculated according to formula (13),
The secondary dendrite arm spacing λ of each structure position of foundry goods can be calculated2Value.
Or, secondary dendrite arm spacing λ2Also can be obtained by the method that precision-investment casting is combined with sunykatuib analyses.Pin
To fixing k4169 alloy, based on fixing cast model and full form casting process condition, carry out actual model casting and cast
Part, carries out subdivision to foundry goods, each structure position is intercepted metallographic specimen, using 10gcucl2+ 100 ethanol+100ml hydrochloric acid
Mix reagent carries out slight erosion, then carries out metallographic observation statistics secondary dendrite arm spacing λ2Meansigma methodss, finally subdivision is cast
The actual secondary dendrite arm spacing λ that part obtains2Value is contrasted with the analogue value, if within both identical or errors are 5%, recognizing
It is the secondary dendrite arm spacing λ of foundry goods for this analog result2Value, if instead then two, differ, adjustment model boundary condition sets
Put, then be simulated, until both are identical or in range of error.
(4) calculate segregation exponent m si value at each structure position of foundry goods
Temperature field data obtaining from procast: setting rate v of molten metal, local setting time tf, secondary dendrite
Spacing λ2, and the thermophysical parameter of the calculating of jmatpro: solute balance partition coefficient k, Liquid Diffusion Coefficient dl, solid-state diffusion
Coefficient ds, substitute in formula (11), and the viewcast module by procast, calculate each structure position of Complicated structure casting
Msi segregation index;According to the size of each structure position msi segregation index, judge the microsegregation of each structure position of foundry goods
Degree;Wherein:
In formula (11): pelFor dimensionless solute p é let number, characterize in the liquid phase of solid liquid interface forward position by solute degree of supersaturation
Determined by driving force;fosFor dimensionless Fourier number, also referred to as the solid-back diffusion factor;K distributes system for solute balance
Number;dlFor Liquid Diffusion Coefficient;dsFor solid phase diffusion welding;V is setting rate;L is characterized diffusion length;tfFor local solidification
Time.
Predict the outcome and model casting obtains the subdivision Comparative result (comparing result is as shown in Figure 3) of foundry goods it is seen then that predicting
Result and subdivision result are coincide, and segregation index is bigger, and microsegregation is more serious, and the interdendritic laves phase that solidification latter stage is formed contains
Amount is higher.
The present invention, by finite element software procast and Thermodynamic Calculation Software jmatpro, in conjunction with foundry practice, can have
The microsegregation degree of the solute element of the effect prediction each structure position of Complicated structure casting;By the application of the inventive method, no
Only can be without dissecting the degree of segregation that foundry goods can predict a certain solute element of each structure position of foundry goods, for optimizing and revising casting
Make technique and Technology for Heating Processing provides foundation, and the method can be to other solute element degree of segregation of multicomponent alloy
Prediction, applied range.
Above the specific embodiment of the present invention is described.It is to be appreciated that the invention is not limited in above-mentioned
Particular implementation, those skilled in the art can make various modifications or modification within the scope of the claims, this not shadow
Ring the flesh and blood of the present invention.