CN103729511B - Method for predicating ingredient segregation degrees in casting process of complex-structure casting - Google Patents

Method for predicating ingredient segregation degrees in casting process of complex-structure casting Download PDF

Info

Publication number
CN103729511B
CN103729511B CN201310746719.5A CN201310746719A CN103729511B CN 103729511 B CN103729511 B CN 103729511B CN 201310746719 A CN201310746719 A CN 201310746719A CN 103729511 B CN103729511 B CN 103729511B
Authority
CN
China
Prior art keywords
casting
segregation
solute
foundry goods
solidification
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.)
Active
Application number
CN201310746719.5A
Other languages
Chinese (zh)
Other versions
CN103729511A (en
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.)
Shanghai Zhongchao Hangyu Precision Casting Technology Co ltd
Zhongchao Hangyu Investment Casting S&t Co
Original Assignee
Shanghai Jiaotong 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 Shanghai Jiaotong University filed Critical Shanghai Jiaotong University
Priority to CN201310746719.5A priority Critical patent/CN103729511B/en
Publication of CN103729511A publication Critical patent/CN103729511A/en
Application granted granted Critical
Publication of CN103729511B publication Critical patent/CN103729511B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Investigating Or Analyzing Materials Using Thermal Means (AREA)

Abstract

The invention provides a method for predicating the ingredient segregation degrees in the casting process of a complex-structure casting. The method includes the steps of (1) calculating to obtain thermophysical parameters of alloy through a thermophysical parameter calculating module of commercialized JMatPro software, (2) carrying out analogue simulation on the casting filling and solidifying processes of the complex-structure casting through ProCAST finite element software, and obtaining temperature field data of various parts of the casting, (3) calculating secondary dendrite arm spacing values of the structural parts of the casting through a secondary dendrite arm spacing calculating module built in the ProCAST software, and (4) calculating MSIs of the structural parts of the complex-structure casting, and judging the micro-segregation degrees of the structural parts of the casting according to the MSIs. According to the method, the finite element software and thermodynamic calculation software are used, the casting practice is combined, the micro-segregation degrees of solute elements of the structural parts of the complex-structure casting are effectively predicted, and the application range is wide.

Description

The Forecasting Methodology of component segregation degree in a kind of Complicated structure casting casting process
Technical field
A kind of the present invention relates to steel industry component segregation analysis technical field, in particular it relates to Complicated structure casting casting The Forecasting Methodology of component segregation degree during making.
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:
m s i = pe l kfo s = 1 8 · 1 2 kd l d s vλ 2 3 t f = 1 8 · a · vλ 2 3 t f , a = 1 2 kd l d s
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:
m = 166 · γd l l n ( c e u t / c 0 ) m ( 1 - k ) ( c 0 - c e u t )
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.
Brief description
The detailed description with reference to the following drawings, non-limiting example made by reading, the further feature of the present invention, Objects and advantages will become more apparent upon:
Fig. 1 is one embodiment of the invention Casting Three-dimensional illustraton of model;
Fig. 2 is one embodiment of the invention foundry goods pictorial diagram;
Fig. 3 is one embodiment of the invention msi segregation exponential forecasting result and laves phase content experimental result comparison diagram.
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 δcc=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
pe l = v l 2 d l - - - ( 7 )
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.
fo s = α = d s t f l 2 - - - ( 8 )
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
m s i = pe l kfo s - - - ( 9 )
WhereinSubstitute into (9)
m s i = pe l kfo s = 1 2 kd l d s vl 3 t f = a · vl 3 t f - - - ( 10 )
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)
m s i = pe l kfo s = 1 8 · 1 2 kd l d s vλ 2 3 t f = 1 8 · a · vλ 2 3 t f , a = 1 2 kd l d s - - - ( 11 )
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:
m = 166 · γd l l n ( c e u t / c 0 ) m ( 1 - k ) ( c 0 - c e u t ) - - - ( 13 )
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:
m s i = pe l kfo s = 1 8 · 1 2 kd l d s vλ 2 3 t f = 1 8 · a · vλ 2 3 t f , a = 1 2 kd l d s - - - ( 11 )
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.

Claims (3)

1. in a kind of Complicated structure casting casting process the Forecasting Methodology of component segregation degree it is characterised in that the method includes Following steps:
(1) obtain calculating thermal physical property parameter
Thermophysical parameter computing module using commercialization jmatpro software obtains the thermal physical property parameter of alloy, and parameter includes molten Matter equilibrium distribution 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 casting The temperature field data at each position of part, data includes setting rate v of molten metal, local setting time tf;Or local solidifies Time tf, the method that 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 Branch of each structure position of foundry goods using the secondary dendrite arm spacing computing module that procast finite element software carries Intergranular is away from λ2;Or secondary dendrite arm spacing λ2Obtained 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
Solidification Parameters obtaining from procast finite element software: setting rate v of molten metal, local setting time tf, secondary Dendrite interval λ2, and the thermophysical parameter that jmatpro software obtains: solute balance partition coefficient k, Liquid Diffusion Coefficient dl, solid Phase diffusion coefficient ds, substitute in microsegregation deciding degree exponent m si formula, and by procast finite element software Viewcast module, calculates segregation exponent m si of each structure position of Complicated structure casting, according to each structure of Complicated structure casting The size of the segregation exponent m si value at position judges the microsegregation degree of each structure position of foundry goods: segregation exponent m si is bigger, micro- Sight segregation is more serious, and it is higher that the interdendritic that solidification latter stage is formed separates out phase content;Wherein:
m s i = pe l kfo s = 1 8 · 1 2 kd l d s vλ 2 3 t f = 1 8 · a · vλ 2 3 t f , a = 1 2 kd l d s
In formula: pelFor dimensionless solute p é let number, characterize in the liquid phase of solid liquid interface forward position and expand determined by solute degree of supersaturation Scattered driving force;fosFor dimensionless Fourier number, also referred to as the solid-back diffusion factor;K is solute balance partition coefficient;dlExpand for liquid phase Scattered coefficient;dsFor solid phase diffusion welding;V is setting rate;tfFor local setting time;For comprising alloy material Physical parameter, depending on the selection of material, for the different structure position on Same Part, a is considered constant;
Described local setting time tf, the method that combined with sunykatuib analyses by precision-investment casting of setting rate v obtain, tool Body is: for fixing k4169 alloy, based on fixing cast model and casting technique condition, carries out actual casting and is cast Part, in casting process, is pre-placed high temperature resisting thermocouple in casting features structure position, measures casting process molten metal in real time Cooling curve, on the one hand, with reference to jmatpro software, k4169 alloy is carried out with phasor and obtain beginning solidification temperature tlAnd solidification knot Shu Wendu ts, you can obtain local setting time tfWith rate of cooling r, wherein: local setting time tfIt is molten metal to start to coagulate Gu terminating the experienced time to solidification, and rate of cooling r passes through formula r=- (ts-tl)/tfObtain;On the other hand, in conjunction with many Prop up the position of thermocouple, obtain thermograde g of appointed part, obtain the solidifying of appointed part molten metal further according to formula r=-gv Gu speed v;
Finally the actual measured value of thermocouple and the analogue value are compared, if both identical then it is assumed that this simulation obtain office Domain setting time tf, setting rate v be foundry goods local setting time tf, setting rate v, if both differ, adjust mould Type boundary condition is arranged, then is simulated, until both are identical;
Described secondary dendrite arm spacing λ2Obtained by the method that precision-investment casting is combined with sunykatuib analyses, particularly as follows: for solid Fixed k4169 alloy, based on fixing cast model and full form casting process condition, carries out actual model casting and obtains foundry goods, Subdivision is carried out to foundry goods, each structure position is intercepted metallographic specimen, using 10gcucl2+ 100ml 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 both differ, adjusts model boundary condition setting, then It is simulated, until both are identical or in range of error.
2. in a kind of Complicated structure casting casting process according to claim 1 component segregation degree Forecasting Methodology, its It is characterised by, in step (3), secondary dendrite arm spacing λ2Computing formula be
λ2=(mtf)1/3
Wherein:
m = 166 · γd l l n ( c e u t / c 0 ) m ( 1 - k ) ( c 0 - c e u t )
In formula, γ is gibbs-thomson coefficient, c0For initial solute concentration, ceutThere is liquid during eutectic reaction for solidification latter stage Phase solute concentration, m is alloy liquid phase line slope, and m is with the related parameter of material.
3. in a kind of Complicated structure casting casting process according to claim 2 component segregation degree Forecasting Methodology, its It is characterised by, 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.
CN201310746719.5A 2013-12-30 2013-12-30 Method for predicating ingredient segregation degrees in casting process of complex-structure casting Active CN103729511B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310746719.5A CN103729511B (en) 2013-12-30 2013-12-30 Method for predicating ingredient segregation degrees in casting process of complex-structure casting

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310746719.5A CN103729511B (en) 2013-12-30 2013-12-30 Method for predicating ingredient segregation degrees in casting process of complex-structure casting

Publications (2)

Publication Number Publication Date
CN103729511A CN103729511A (en) 2014-04-16
CN103729511B true CN103729511B (en) 2017-01-25

Family

ID=50453585

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310746719.5A Active CN103729511B (en) 2013-12-30 2013-12-30 Method for predicating ingredient segregation degrees in casting process of complex-structure casting

Country Status (1)

Country Link
CN (1) CN103729511B (en)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104014768B (en) * 2014-06-24 2015-11-18 哈尔滨理工大学 A kind of method of magnesium alloy arborescent structure numerical simulation
US20160034614A1 (en) * 2014-08-01 2016-02-04 GM Global Technology Operations LLC Materials property predictor for cast aluminum alloys
CN104881588B (en) * 2015-06-19 2017-11-10 哈尔滨理工大学 Macroscopic segregation of cast ingot method for numerical simulation
CN105354372B (en) * 2015-10-22 2019-01-22 上海交通大学 A kind of prognosis modelling method being segregated in steel ingot
CN108108529B (en) * 2017-12-01 2021-07-06 东方电气集团东方汽轮机有限公司 Inverse calculation method for simply and conveniently measuring heat exchange coefficient of casting interface
CN110111861B (en) * 2019-05-24 2020-12-22 上海交通大学 Prediction method for thermal cracks in solidification process of magnesium and aluminum alloy castings
CN110263418B (en) * 2019-06-17 2022-10-21 哈尔滨理工大学 Body-centered cubic alloy microsegregation numerical prediction method
CN111539094B (en) * 2020-04-10 2023-04-18 苏州大学 Numerical simulation method for segregation of MIG (Metal-inert gas) welding area of dissimilar welding wire
CN111872324B (en) * 2020-06-23 2021-05-18 上海交通大学 Parameter acquisition method for casting solidification simulation and gridding design method of casting system
CN113263154B (en) * 2021-05-14 2022-11-25 北京理工大学珠海学院 Method, device, equipment and storage medium for predicting metal dendrite spacing
CN114637954B (en) * 2022-03-25 2023-02-07 宁夏中欣晶圆半导体科技有限公司 Method for calculating axial distribution of carbon content of crystal bar

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102289542A (en) * 2011-07-12 2011-12-21 哈尔滨理工大学 Method for macrosegregation numerical simulation of casting

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003033864A (en) * 2001-07-17 2003-02-04 Mitsubishi Heavy Ind Ltd Simulator for casting process and judging method therefor

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102289542A (en) * 2011-07-12 2011-12-21 哈尔滨理工大学 Method for macrosegregation numerical simulation of casting

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Xu Li等."Influence of a weak static magnetic field on the primary dendrite arm spacing of a directionally solidified Ni-based superalloy".《Materials Letters》.2012,第67卷(第1期), *
Yongzhi Zhang等."Precipitation evolution of Al–Zr–Yb alloys during isochronal aging".《Scripta Materialia》.2013,第69卷(第6期), *
王恺等."合金化元素对K4169高温合金组织稳定性影响的研究现状".《材料导报》.2013,第27卷(第5期), *
胡红军等."ProCAST软件在铸造凝固模拟中的应用".《材料科学与工艺》.2006,第14卷(第3期), *

Also Published As

Publication number Publication date
CN103729511A (en) 2014-04-16

Similar Documents

Publication Publication Date Title
CN103729511B (en) Method for predicating ingredient segregation degrees in casting process of complex-structure casting
CN102750425B (en) Simulation method of texture evolvement of heat affected zone in welding process
CN105057642B (en) Casting crystalline grain organizes the formation of the analogy method of correlation values
CN107309543A (en) A kind of laser weld pools dendritic growth analogy method
Pan et al. Modeling of grain selection during directional solidification of single crystal superalloy turbine blade castings
CN105665684B (en) A kind of method of casting crystalline grain tissue values prediction
Zhu et al. Cellular automaton modeling of microporosity formation during solidification of aluminum alloys
CN102426622B (en) Adaptive variable-speed drawing simulation method for production of single-crystal blade
CN107092754B (en) Alloy grain structure numerical value prediction method
Liu et al. Advances in multi-scale modeling of solidification and casting processes
CN103521731A (en) Method for determining tundish outlet concentration through steel ingot multi-tundish concentration-changing continuous pouring
CN101329290A (en) Method for forecasting different constituent elements segregation distribution in the multi-constituent element alloy graining process
Zhu et al. Computational modeling of microstructure evolution in solidification of aluminum alloys
Zhao et al. Effect of the withdrawal rate on the microstructure and thermal durability of a third-generation single crystal superalloy
Cho et al. Modeling of microstructural evolution in squeeze casting of an Al-4.5 mass% Cu alloy
Geng et al. Correlation of crystal orientation, crystal morphology and mechanical properties of directionally solidified Mg-xGd alloys
Liu et al. A phase-field study of the pattern selection between dendrite and seaweed during directional solidification
Danilov et al. Dendritic to globular morphology transition in ternary alloy solidification
CN110993038A (en) Prediction method for dendritic crystal growth in standing molten steel based on parallel computation
Lopez-Galilea et al. Microsegregation and secondary phase formation during directional solidification of the single-crystal Ni-based superalloy LEK94
ZHU et al. Numerical simulation of recalescence of 3-dimensional isothermal solidification for binary alloy using phase-field approach
Tveito et al. The effect of finite microscopic liquid solute diffusion on macrosegregation formation
CN102156184B (en) Method for predicting space between aluminium-silicon alloy eutectic structure lamellas
Ma et al. Evolution of the temperature, microstructures and microsegregation in equiaxed solidification of Al-5 wt.% Cu alloy
An et al. Control Center Segregation in Continuously Cast GCr15 Bloom by Optimization of Solidification Structure

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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20191211

Address after: 201203 room 105g, No.199, GuoShouJing Road, Shanghai pilot Free Trade Zone, Pudong New Area, Shanghai (biological and pharmaceutical innovation building)

Co-patentee after: Wang Jun

Patentee after: Shanghai Jiaotong University Intellectual Property Management Co.,Ltd.

Address before: 200240 Dongchuan Road, Shanghai, No. 800, No.

Patentee before: SHANGHAI JIAO TONG University

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20200720

Address after: 214241 Jiangsu city of Wuxi province Xu Yixing City Industrial Zone

Patentee after: ZHONGCHAO HANGYU INVESTMENT CASTING S&T Co.

Address before: Room 111, NO.201 Qinzhou Road, Xuhui District, Shanghai, 200235

Patentee before: SHANGHAI ZHONGCHAO HANGYU PRECISION CASTING TECHNOLOGY Co.,Ltd.

Effective date of registration: 20200720

Address after: Room 111, NO.201 Qinzhou Road, Xuhui District, Shanghai, 200235

Patentee after: SHANGHAI ZHONGCHAO HANGYU PRECISION CASTING TECHNOLOGY Co.,Ltd.

Address before: 201203 room 105g, No.199, GuoShouJing Road, Shanghai pilot Free Trade Zone, Pudong New Area, Shanghai (biological and pharmaceutical innovation building)

Co-patentee before: Wang Jun

Patentee before: Shanghai Jiaotong University Intellectual Property Management Co.,Ltd.