CN103729511A - 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
CN103729511A
CN103729511A CN201310746719.5A CN201310746719A CN103729511A CN 103729511 A CN103729511 A CN 103729511A CN 201310746719 A CN201310746719 A CN 201310746719A CN 103729511 A CN103729511 A CN 103729511A
Authority
CN
China
Prior art keywords
casting
foundry goods
solute
arm spacing
dendrite arm
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
Application number
CN201310746719.5A
Other languages
Chinese (zh)
Other versions
CN103729511B (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

Images

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
The present invention relates to steel industry component segregation analysis technical field, particularly, relate to the Forecasting Methodology of component segregation degree in a kind of Complicated structure casting casting process.
Background technology
Component segregation is aluminium alloy intrinsic characteristic while solidifying.When the solubleness of solute element in solid phase is less than its solubleness in liquid phase, solid-liquid interface will be got rid of solute in liquid phase; And when the solubleness of solute element in solid phase is greater than the solubleness in liquid phase, solute, from liquid phase to solid-state diffusion, will cause solid-liquid interface forward position to occur the poor district of solute.This reallocation due to solute in alloy graining process causes the uneven components after alloy graining, is component segregation.Conventionally, the solubleness of solute element in liquid phase is than high in solid phase, so solute has the solidifying phase of disengaging and be pulled to the tendency in the liquid of forward position, and end is solidified the grown dendritic arm of liquid that part is rich in solute and cut apart, forming component segregation region.Component segregation is one of major defect of foundry goods, reduces alloy solid solution strengthening effect, and meanwhile, the low melting point precipitated phase that interdendritic forms easily becomes formation of crack in the welding process of foundry goods, reduces mechanical property and the welding performance of foundry goods.
Known according to classical solidification theory, under general alloy graining condition, the solute redistribution behavior of micro-scale and the microsegregation of generation not only depend on that the physical parameter of alloy is (as solid-state diffusion coefficient D s, Liquid Diffusion Coefficient D l, solute distribution coefficient k) and external process condition as freezing rate V(or local solidification time t f), thermograde G, but also closely related with factors such as the growth pattern of dendrite, dendrite morphologies.
The simplest analytic model is Lever analytic model and the Scheil analytic model under partial balancing's curing condition, and other is all correction type analytic model.Consideration based on back diffusion, lever model and Scheil model are two kinds of limiting cases.
Lever model is a kind of equilibrium freezing model, supposes that solute is in all fully diffusions of solid-liquid two-phase,
C L=C 0[1-(1-k)f s] -1
Formula 1. in, C lfor liquid phase solute concentration, C 0for initial solute concentration, f sfor fraction solid, k is solute distribution coefficient.
Scheil model is also referred to as non-equilibrium lever model, wherein suppose solute in liquid phase fully diffusion and solid phase without diffusion, in model, do not consider that forming core is excessively cold, the formation of dendritic arm alligatoring, macroscopic material transmission and pore,
C L=C 0(1-f s) (k-1)
Obviously, as fraction solid f s=1 o'clock C ltrend towards infinity, so Scheil model can not accurately calculate to final solute concentration.
H.D.Brody & M.C, Flemings considers the antidiffusion of solute element in process of setting, has
C L=C 0[1-(1-2αk)f s] (k-1)/(1-2αk)
Formula 3. in, α is nondimensional solute invasin, or claims nondimensional diffusion time (Fourier number), it has represented the diffusion boundary layer thickness δ in solid phase swith the proportionate relationship of system dimension L, can be expressed as:
α = D s t f L 2
When α=0 is nothing diffusion in solid phase, 3. formula is reduced to Scheil model;
When α=0.5,3. formula is reduced to lever model.
The people such as W.Kurz revise α value, propose with α ' replacement α,
α ′ = 2 α [ 1 - exp ( - 1 α ) ] - exp ( - 1 2 α )
Although scholars have studied most advanced and sophisticated cold, the impact that the solute assigning processs such as cold-peace coagulated volume variation are crossed at eutectic interface of coagulating property, the solid-back diffusion in process of setting, the limited diffusion of liquid, dendritic arm alligatoring, dendrite that the alloying component of multicomponent alloy changes energetically with excessivelying in recent years, set up in some process of setting solute concentration with microsegregation analytic model and the numerical model of fraction solid increase.But the microsegregation model that existing solute distributes is all based upon the variation that in research process of setting, solute concentration increases with fraction solid, can only represent interfacial concentration
Figure BDA0000450257850000023
with fraction solid f sfuntcional relationship, can not determine accordingly the component distributing of final solid phase, and must consider this dynamic parameter of fraction solid, and can not judge for a large and complex structure part different structure position, by the segregation degree under many factors acting in conjunction, cannot instruct the Foundry Production of large and complex structure foundry goods.
Summary of the invention
For defect of the prior art, the Forecasting Methodology that the object of this invention is to provide component segregation degree in a kind of Complicated structure casting casting process, this Forecasting Methodology can predict that parts with complex structures different structure position, in the microsegregation degree of different curing conditions, has directive significance to the Foundry Production of large and complex structure foundry goods.
For realizing above object, the invention provides the Forecasting Methodology of component segregation degree in a kind of Complicated structure casting casting process, the method comprises:
(1) obtain calculating thermal physical property parameter
Utilize the thermophysical parameter computing module of commercialization JMatPro software to calculate the thermal physical property parameter that obtains alloy, parameter comprises solute balance partition factor k, Liquid Diffusion Coefficient D l, solid-state diffusion coefficient D s, alloy initial solidification temperature T l, eutectic reaction temperature T s;
(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 the temperature field data at each position of foundry goods, and data comprise setting rate V, the local setting time t of molten metal f; Or local setting time t f, the setting rate V method that also can combine with sunykatuib analysis by precision-investment casting obtains.
(3) obtain the secondary dendrite arm spacing λ of each structure position place of foundry goods 2
Utilize secondary dendrite arm spacing computing module that ProCAST software carries to calculate the secondary dendrite arm spacing λ of each structure position of foundry goods 2; Or, secondary dendrite arm spacing λ 2the method that also can combine with sunykatuib analysis by precision-investment casting is obtained.
(4) calculate each structure position place segregation index M of foundry goods SI value
The Solidification Parameters obtaining from ProCAST: the setting rate V of molten metal, local setting time t f, secondary dendrite arm spacing λ 2, and the thermophysical parameter of JMatPro calculating: solute balance partition factor k, Liquid Diffusion Coefficient D l, solid-state diffusion coefficient D sin substitution microsegregation deciding degree index M SI formula, and by the ViewCAST module of ProCAST, calculate the segregation index M SI of each structure position of Complicated structure casting, according to the size of the segregation index M SI value of each structure position of Complicated structure casting, judge the microsegregation degree of each structure position of foundry goods: SI is larger for segregation index M, microsegregation is more serious, solidifies the interdendritic precipitated phase content forming latter stage just higher; Wherein:
MSI = Pe L kF o S = 1 8 · 1 2 k D L D S V λ 2 3 t f = 1 8 · A · V λ 2 3 t f , a = 1 2 k D L D S
In formula: Pe lfor dimensionless solute P é let number, characterize in the liquid phase of solid-liquid interface forward position by the determined driving force of solute degree of supersaturation; Fo sfor dimensionless Fourier number, also become the solid-back diffusion factor; K is solute balance partition factor; D lfor Liquid Diffusion Coefficient; D sfor solid-state diffusion coefficient; V is setting rate; L is feature diffusion length; t ffor local solidification time;
Figure BDA0000450257850000032
for comprising the physical parameter of alloy material, depend on the selection of material, for the different structure position on Same Part, A can think constant.
Preferably, in step (2), secondary dendrite arm spacing λ 2computing formula be
λ 2=(Mt f) 1/3
Wherein:
M = 166 · ΓD l ln ( C eut / C 0 ) m ( 1 - k ) ( C 0 - C eut )
In formula, Γ is Gibbs-Thomson coefficient, C 0for initial solute concentration, C eutliquid phase solute concentration while there is eutectic reaction latter stage for solidifying, m is alloy liquid phase line slope, M is for following the relevant parameter of material.
More preferably, M value is calculated by 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. acquisition.
Preferably, at ProCAST, calculate in panel, according to different alloys, input corresponding initial solidification temperature T l, eutectic reaction temperature T sand the M value calculating, can calculate the secondary dendrite arm spacing value λ of each structure position of foundry goods 2.
Compared with prior art, the present invention has following beneficial effect:
The present invention, by finite element software ProCAST and Thermodynamic Calculation Software JMatPro, can effectively predict the microsegregation degree of the solute element of each structure position of Complicated structure casting; By the application of the inventive method, not only can dissect foundry goods and can predict the segregation degree of a certain solute element of each structure position of foundry goods, for optimizing and revising casting technique and Technology for Heating Processing provides foundation, and the method can be to other solute element segregation degree of multicomponent alloy prediction, applied range.
Accompanying drawing explanation
By reading the detailed description of non-limiting example being done with reference to the following drawings, it is more obvious that other features, objects and advantages of the present invention will become:
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.
Embodiment
Below in conjunction with specific embodiment, the present invention is described in detail.Following examples will contribute to those skilled in the art further to understand the present invention, but not limit in any form the present invention.It should be pointed out that to those skilled in the art, without departing from the inventive concept of the premise, can also make some distortion and improvement.These all belong to protection scope of the present invention.
In the present embodiment, microsegregation deciding degree index M SI's (Micro-segregation Index) determines based on following consideration.
Theoretical by means of widely used Buckingham π non-dimension analysis in mathematical physics, solute in process of setting is distributed and carries out non-dimension analysis.The research of diffusion flux based on solid-liquid interface both sides liquid phase and solid phase solute atoms, the behavior that the antidiffusion alloy that the present invention considers the limited diffusion of solid-liquid interface forward position liquid phase solute and solid phase solute solidifies solute redistribution is studied.The employing P é let number of solid-liquid interface forward position solute diffusion represents, P é let number is dendrite end radius of curvature R and liquid phase solute Equivalent Boundary layer thickness δ cratio, liquid phase solute Equivalent Boundary layer thickness δ cc=2D l/ V) represent that in diffusion length, contained whole solutes equal the contained solute total amount of infinite boundary layer.Solute concentration gradient due to dendrite end, can see approx it is the solute concentration gradient in the spheroid forward position of growing up, according to the way that solves in diffusion field around the globular crystal of growing, if only consider radial diffusion, the globular crystal of growing boundary layer thickness around equals the radius of ball, and the solute diffusion length L of dendrite end is approximately equal to the radius R of ball.Therefore, P é let number can be expressed as
Pe L = VL 2 D L - - - ( 7 )
P é let number has characterized in the liquid phase of solid-liquid interface forward position by the determined driving force of solute degree of supersaturation.
Owing to usually there is eutectic in actual process of setting, liquid concentration can not continue to increase, and in order to describe truly the component distributing of solidification end, must consider the antidiffusion of solute atoms in solid phase.According to the principle of mass conservation, in solid phase, the speed degree of solute atoms antidiffusion depends on dimensionless group α, also referred to as Fo s, it has characterized the ability of solute atoms antidiffusion in solid phase.
Fo S = α = D S t f L 2 - - - ( 8 )
According to the solute mass conservation, the segregation index M SI (Micro-segregation Index) that characterizes the solute segregation degree in process of setting can be expressed as
MSI = Pe L k Fo S - - - ( 9 )
Wherein Pe L = VL 2 D L , Fo s = D S t f L 2 , Substitution (9)
MSI = Pe L k Fo S = 1 2 K D L D S VL 3 t f = A · VL 3 t f - - - ( 10 )
Wherein
Figure BDA0000450257850000056
be the physical parameter that comprises alloy material, depend on the selection of material.For the different structure position on Same Part, A can think constant.As seen from formula (10), the segregation index M SI value of the judgement microsegregation degree that the present invention proposes size, not only follows the thermal physical property parameter solute balance partition factor k of alloy, Liquid Diffusion Coefficient D l, solid-state diffusion coefficient D sbe correlated with, and also have local solidification condition setting rate V, feature diffusion length L and the local solidification time t of alloy frelevant.
For general casting process, dendrite is to wait axle dendrite mode to grow.Now, secondary dendrite arm spacing λ 2it is the characteristic length of solute diffusion.Therefore, half of getting secondary dendrite arm spacing is diffusion length,
Figure BDA0000450257850000061
Substitution formula (10)
MSI = Pe L kF o S = 1 8 · 1 2 k D L D S V λ 2 3 t f = 1 8 · A · V λ 2 3 t f , a = 1 2 k D L D S - - - ( 11 )
Be below the present invention's one application the present embodiment:
The present embodiment provides the Forecasting Methodology of component segregation degree in a kind of Complicated structure casting casting process, and the K4169 nickel base superalloy large-scale complex thin-wall foundry goods of casting is carried out to the prediction of microsegregation degree.The large-area thin plate that this casting design contains different-thickness, variable cross section, the structures such as cylinder.As shown in Figure 1, the foundry goods pictorial diagram that cast obtains as shown in Figure 2 for Casting Three-dimensional illustraton of model.The alloying element adding due to K4169 nickel base superalloy is many, and alloying level is high, and solidification temperature range is wide, and the solidifying segregation of alloying element is serious.Microsegregation is one of major defect of K4169 nickel base superalloy foundry goods, and the mechanical property of foundry goods, welding performance and processing characteristics are had to important impact.In K4169 nickel base superalloy casting solidification process, the solubleness of Nb solute atoms in matrix is little, in process of setting, solid-liquid interface will be got rid of a large amount of Nb solute atomss to liquid phase, while being cooled to 1180 ℃, there is γ+Laves eutectic reaction, in interdendritic, form Laves phase.If due to alloy graining process component segregation, solidify the Laves phase amount forming latter stage more, no matter their form and how distributing, because Laves has fixed the intensified elements such as a large amount of Nb in mutually, thereby reduce solid solution strengthening effect.Meanwhile, also reduced hardening constituent γ ' or γ " quantity, weakened precipitation strengthening effect, cause alloy mechanical property to worsen.Therefore, eliminate solute element Nb segregation and Laves is a gordian technique difficult problem for K4169 nickel base superalloy castings production mutually always.
(1) obtain calculating thermal physical property parameter
By JMatPro software, K4169 alloy is carried out to PHASE DIAGRAM CALCULATION, and in conjunction with 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 and calculates each thermophysical parameter used, and each parameter is as follows:
Γ=3.65×10 -7(km),k Nb=0.48,D l,Nb=3×10 -9(m 2s -1),D s,Nb=2.82×10 -13(m 2s -1),C 0,Nb=4.316(wt.pct),C eut,Nb=19.1(wt.pct),m l=-10.9,T l=1360℃,T s=1180℃;
(2) obtain the temperature field data at each position of foundry goods
Application ProCAST finite element software carries out analogue simulation to casting filling, the process of setting of foundry goods as shown in Figure 1, obtains the temperature field data at each position of foundry goods, as the local setting time t of molten metal f, setting rate V;
Or local setting time t f, the setting rate V method that also can combine with sunykatuib analysis by precision-investment casting obtains.For fixing K4169 alloy, the cast model based on fixing and full form casting process condition, carry out actual investment cast and obtain foundry goods.In casting process, at foundry goods feature structure position, place in advance high temperature resistance tungsten-rhenium thermal couple, measure in real time the cooling curve of casting process molten metal.On the one hand, in conjunction with JMatPro software, K4169 alloy is carried out the beginning temperature of solidification T of PHASE DIAGRAM CALCULATION acquisition land solidify end temp T s, can obtain local setting time t fwith cooling velocity R, local setting time t fbe molten metal and start to solidify the time of experiencing to solidifying end, and cooling velocity R is by formula R=-(T s-T l)/t fobtain.On the other hand, in conjunction with the equidistant position of the branched thermopair of layout, can obtain the thermograde G of appointed part, then according to formula R=-GV, can obtain the setting rate V of appointed part molten metal.Finally the actual measured value of thermopair and the analogue value are compared, if both are identical, think the local setting time t that this simulation obtains f, setting rate V is the local setting time t of foundry goods f, setting rate V, if instead two not identical, adjustment model boundary condition arranges, then simulates, until both are identical.
(3) obtain the secondary dendrite arm spacing λ of each structure position place of foundry goods 2
Utilize secondary dendrite arm spacing computing module that ProCAST software carries to calculate the λ of each structure position of foundry goods 2value.The M value that obtains as calculated K4169 alloy is 19.43(um 3s -1); Secondary dendrite arm spacing λ 2computing formula be
λ 2=(Mt f) 1/3 (12)
Wherein:
M = 166 · ΓD l ln ( C eut / C 0 ) m ( 1 - k ) ( C 0 - C eut ) - - - ( 13 )
In formula (13), Γ is Gibbs-Thomson coefficient, C 0for initial solute concentration, C eutliquid phase solute concentration while there is eutectic reaction latter stage for solidifying, m is alloy liquid phase line slope, M is for following the relevant parameter of material.At ProCAST, calculate in panel, according to different alloys, input corresponding initial solidification temperature T l, eutectic reaction temperature T s, and the M value calculating according to formula (13), can calculate the secondary dendrite arm spacing λ of each structure position of foundry goods 2value.
Or, secondary dendrite arm spacing λ 2the method that also can combine with sunykatuib analysis by precision-investment casting is obtained.For fixing K4169 alloy, the cast model based on fixing and full form casting process condition, carry out actual investment cast and obtain foundry goods, and foundry goods is carried out to subdivision, by each structure position intercepting metallographic specimen, adopts 10gCuCl 2the mix reagent of+100 alcohol+100ml hydrochloric acid slightly corrodes, and then carries out metallographic observation statistics secondary dendrite arm spacing λ 2mean value, the actual secondary dendrite arm spacing λ finally subdivision foundry goods being obtained 2value and the analogue value contrast, if both identical or errors are in 5%, think that this analog result is the secondary dendrite arm spacing λ of foundry goods 2value, if instead two not identical, adjustment model boundary condition arranges, then simulates, until both are identical or in error range.
(4) calculate each structure position place segregation index M of foundry goods SI value
The temperature field data of obtaining from ProCAST: the setting rate V of molten metal, local setting time t f, secondary dendrite arm spacing λ 2, and the thermophysical parameter of the calculating of JMatPro: solute balance partition factor k, Liquid Diffusion Coefficient D l, solid-state diffusion coefficient D s, in substitution formula (11), and by the ViewCAST module of ProCAST, calculate the MSI segregation index of each structure position of Complicated structure casting; According to the size of each structure position MSI segregation index, judge the microsegregation degree of each structure position of foundry goods; Wherein:
MSI = Pe L kF o S = 1 8 · 1 2 k D L D S V λ 2 3 t f = 1 8 · A · V λ 2 3 t f , a = 1 2 k D L D S - - - ( 11 )
In formula (11): Pe lfor dimensionless solute P é let number, characterize in the liquid phase of solid-liquid interface forward position by the determined driving force of solute degree of supersaturation; Fo sfor dimensionless Fourier number, also become the solid-back diffusion factor; K is solute balance partition factor; D lfor Liquid Diffusion Coefficient; D sfor solid-state diffusion coefficient; V is setting rate; L is feature diffusion length; t ffor local solidification time.
Predict the outcome and obtain the subdivision result contrast (comparing result is as shown in Figure 3) of foundry goods with investment cast, visible, predict the outcome and coincide with subdivision result, segregation index is larger, and microsegregation is more serious, and the interdendritic Laves phase content that solidifies formation in latter stage is just higher.
The present invention, by finite element software ProCAST and Thermodynamic Calculation Software JMatPro, in conjunction with foundry practice, can effectively predict the microsegregation degree of the solute element of each structure position of Complicated structure casting; By the application of the inventive method, not only can dissect foundry goods and can predict the segregation degree of a certain solute element of each structure position of foundry goods, for optimizing and revising casting technique and Technology for Heating Processing provides foundation, and the method can be to other solute element segregation degree of multicomponent alloy prediction, applied range.
Above specific embodiments of the invention are described.It will be appreciated that, the present invention is not limited to above-mentioned specific implementations, and those skilled in the art can make various distortion or modification within the scope of the claims, and this does not affect flesh and blood of the present invention.

Claims (6)

1. a Forecasting Methodology for component segregation degree in Complicated structure casting casting process, is characterized in that, the method comprises the steps:
(1) obtain calculating thermal physical property parameter
The thermal physical property parameter that utilizes the thermophysical parameter computing module acquisition alloy of commercialization JMatPro software, parameter comprises solute balance partition factor k, Liquid Diffusion Coefficient D l, solid-state diffusion coefficient D s, alloy initial solidification temperature T l, eutectic reaction temperature T s;
(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 the temperature field data at each position of foundry goods, and data comprise setting rate V, the local setting time t of molten metal f; Or local setting time t f, the method that combines with sunykatuib analysis by precision-investment casting of setting rate V obtains;
(3) obtain the secondary dendrite arm spacing λ of each structure position place of foundry goods 2
Utilize secondary dendrite arm spacing computing module that ProCAST finite element software carries to calculate the secondary dendrite arm spacing λ of each structure position of foundry goods 2; Or secondary dendrite arm spacing λ 2the method combining with sunykatuib analysis by precision-investment casting is obtained;
(4) calculate each structure position place segregation index M of foundry goods SI value
The Solidification Parameters obtaining from ProCAST finite element software: the setting rate V of molten metal, local setting time t f, secondary dendrite arm spacing λ 2, and the thermophysical parameter that obtains of JMatPro software: solute balance partition factor k, Liquid Diffusion Coefficient D l, solid-state diffusion coefficient D sin substitution microsegregation deciding degree index M SI formula, and by the ViewCAST module of ProCAST finite element software, calculate the segregation index M SI of each structure position of Complicated structure casting, according to the size of the segregation index M SI value of each structure position of Complicated structure casting, judge the microsegregation degree of each structure position of foundry goods: SI is larger for segregation index M, microsegregation is more serious, solidifies the interdendritic precipitated phase content forming latter stage just higher; Wherein:
MSI = Pe L kF o S = 1 8 · 1 2 k D L D S V λ 2 3 t f = 1 8 · A · V λ 2 3 t f , a = 1 2 k D L D S
In formula: Pe lfor dimensionless solute P é let number, characterize in the liquid phase of solid-liquid interface forward position by the determined driving force of solute degree of supersaturation; Fo sfor dimensionless Fourier number, also become the solid-back diffusion factor; K is solute balance partition factor; D lfor Liquid Diffusion Coefficient; D sfor solid-state diffusion coefficient; V is setting rate; L is feature diffusion length; t ffor local solidification time;
Figure FDA0000450257840000021
for comprising the physical parameter of alloy material, depend on the selection of material, for the different structure position on Same Part, A thinks constant.
2. the Forecasting Methodology of component segregation degree in a kind of Complicated structure casting casting process according to claim 1, is characterized in that described local setting time t f, the method that combines with sunykatuib analysis by precision-investment casting of setting rate V obtains, be specially: for fixing K4169 alloy, cast model based on fixing and casting technique condition, carry out actual casting and obtain foundry goods, in casting process, at foundry goods feature structure position, place in advance high temperature resisting thermocouple, measure in real time the cooling curve of casting process molten metal, on the one hand, in conjunction with JMatPro software, K4169 alloy is carried out to phasor and obtain beginning temperature of solidification T land solidify end temp T s, can obtain local setting time t fwith cooling velocity R, wherein: local setting time t fbe molten metal and start to solidify the time of experiencing to solidifying end, and cooling velocity R is by formula R=-(T s-T l)/t fobtain; On the other hand, in conjunction with the position of branched thermopair, obtain the thermograde G of appointed part, then according to formula R=-GV, obtain the setting rate V of appointed part molten metal;
Finally the actual measured value of thermopair and the analogue value are compared, if both are identical, think the local setting time t that this simulation obtains f, setting rate V is the local setting time t of foundry goods f, setting rate V, if instead two not identical, adjustment model boundary condition arranges, then simulates, until both are identical.
3. the Forecasting Methodology of component segregation degree in a kind of Complicated structure casting casting process according to claim 1, is characterized in that described secondary dendrite arm spacing λ 2the method combining with sunykatuib analysis by precision-investment casting is obtained, be specially: for fixing K4169 alloy, cast model based on fixing and full form casting process condition, carry out actual investment cast and obtain foundry goods, foundry goods is carried out to subdivision, by each structure position intercepting metallographic specimen, adopt 10gCuCl 2the mix reagent of+100 alcohol+100ml hydrochloric acid slightly corrodes, and then carries out metallographic observation statistics secondary dendrite arm spacing λ 2mean value, the actual secondary dendrite arm spacing λ finally subdivision foundry goods being obtained 2value and the analogue value contrast, if both identical or errors are in 5%, think that this analog result is the secondary dendrite arm spacing λ of foundry goods 2value, if instead two not identical, adjustment model boundary condition arranges, then simulates, until both are identical or in error range.
4. the Forecasting Methodology of component segregation degree in a kind of Complicated structure casting casting process according to claim 3, is characterized in that, in step (2), and secondary dendrite arm spacing λ 2computing formula be
λ 2=(Mt f) 1/3
Wherein:
M = 166 · ΓD l ln ( C eut / C 0 ) m ( 1 - k ) ( C 0 - C eut )
In formula, Γ is Gibbs-Thomson coefficient, C 0for initial solute concentration, C eutliquid phase solute concentration while there is eutectic reaction latter stage for solidifying, m is alloy liquid phase line slope, M is for following the relevant parameter of material.
5. the Forecasting Methodology of component segregation degree in a kind of Complicated structure casting casting process according to claim 4, it is characterized in that, M value is calculated by 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. acquisition.
6. according to the Forecasting Methodology of component segregation degree in a kind of Complicated structure casting casting process described in claim 1-5 any one, it is characterized in that, at ProCAST, calculate in panel, according to different alloys, input corresponding initial solidification temperature T l, eutectic reaction temperature T sand the M value calculating, can calculate the secondary dendrite arm spacing value λ of each structure position of foundry goods 2.
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 true CN103729511A (en) 2014-04-16
CN103729511B 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)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104014768A (en) * 2014-06-24 2014-09-03 哈尔滨理工大学 Numerical modeling method for magnesium alloy dendritic structure
CN104881588A (en) * 2015-06-19 2015-09-02 哈尔滨理工大学 Ingot casting macrosegregation numerical simulation method
CN105320804A (en) * 2014-08-01 2016-02-10 通用汽车环球科技运作有限责任公司 Material property predictor for cast aluminum alloys
CN105354372A (en) * 2015-10-22 2016-02-24 上海交通大学 Prediction simulation method for segregation in steel ingot
CN108108529A (en) * 2017-12-01 2018-06-01 东方电气集团东方汽轮机有限公司 A kind of reverse calculation algorithms of the easy measurement cast interface coefficient of heat transfer
CN110111861A (en) * 2019-05-24 2019-08-09 上海交通大学 The prediction technique of fire check in a kind of magnesium, aluminium alloy castings process of setting
CN110263418A (en) * 2019-06-17 2019-09-20 哈尔滨理工大学 A kind of body centred cubic alloy microsegregation Numerical Predicting Method
CN111539094A (en) * 2020-04-10 2020-08-14 苏州大学 Numerical simulation method for segregation of MIG (Metal-inert gas) welding area of dissimilar welding wire
CN111872324A (en) * 2020-06-23 2020-11-03 上海交通大学 Parameter acquisition method for casting solidification simulation and gridding design method of casting system
CN113263154A (en) * 2021-05-14 2021-08-17 北京理工大学珠海学院 Method, device, equipment and storage medium for predicting metal dendrite spacing
CN114637954A (en) * 2022-03-25 2022-06-17 宁夏中欣晶圆半导体科技有限公司 Method for calculating axial distribution of carbon content of crystal bar

Citations (2)

* 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
CN102289542A (en) * 2011-07-12 2011-12-21 哈尔滨理工大学 Method for macrosegregation numerical simulation of casting

Patent Citations (2)

* 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
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》, vol. 67, no. 1, 15 January 2012 (2012-01-15) *
YONGZHI ZHANG等: ""Precipitation evolution of Al–Zr–Yb alloys during isochronal aging"", 《SCRIPTA MATERIALIA》, vol. 69, no. 6, 30 September 2013 (2013-09-30) *
王恺等: ""合金化元素对K4169高温合金组织稳定性影响的研究现状"", 《材料导报》, vol. 27, no. 5, 10 May 2013 (2013-05-10) *
胡红军等: ""ProCAST软件在铸造凝固模拟中的应用"", 《材料科学与工艺》, vol. 14, no. 3, 30 June 2006 (2006-06-30) *

Cited By (20)

* 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
CN104014768A (en) * 2014-06-24 2014-09-03 哈尔滨理工大学 Numerical modeling method for magnesium alloy dendritic structure
CN105320804A (en) * 2014-08-01 2016-02-10 通用汽车环球科技运作有限责任公司 Material property predictor for cast aluminum alloys
CN104881588A (en) * 2015-06-19 2015-09-02 哈尔滨理工大学 Ingot casting macrosegregation numerical simulation method
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
CN105354372A (en) * 2015-10-22 2016-02-24 上海交通大学 Prediction simulation method for segregation in steel ingot
CN108108529A (en) * 2017-12-01 2018-06-01 东方电气集团东方汽轮机有限公司 A kind of reverse calculation algorithms of the easy measurement cast interface coefficient of heat transfer
CN108108529B (en) * 2017-12-01 2021-07-06 东方电气集团东方汽轮机有限公司 Inverse calculation method for simply and conveniently measuring heat exchange coefficient of casting interface
CN110111861A (en) * 2019-05-24 2019-08-09 上海交通大学 The prediction technique of fire check in a kind of magnesium, aluminium alloy castings process of setting
CN110111861B (en) * 2019-05-24 2020-12-22 上海交通大学 Prediction method for thermal cracks in solidification process of magnesium and aluminum alloy castings
CN110263418A (en) * 2019-06-17 2019-09-20 哈尔滨理工大学 A kind of body centred cubic alloy microsegregation Numerical Predicting Method
CN110263418B (en) * 2019-06-17 2022-10-21 哈尔滨理工大学 Body-centered cubic alloy microsegregation numerical prediction method
CN111539094A (en) * 2020-04-10 2020-08-14 苏州大学 Numerical simulation method for segregation of MIG (Metal-inert gas) welding area of dissimilar welding wire
CN111539094B (en) * 2020-04-10 2023-04-18 苏州大学 Numerical simulation method for segregation of MIG (Metal-inert gas) welding area of dissimilar welding wire
CN111872324A (en) * 2020-06-23 2020-11-03 上海交通大学 Parameter acquisition method for casting solidification simulation and gridding design method of casting system
CN113263154A (en) * 2021-05-14 2021-08-17 北京理工大学珠海学院 Method, device, equipment and storage medium for predicting metal dendrite spacing
CN113263154B (en) * 2021-05-14 2022-11-25 北京理工大学珠海学院 Method, device, equipment and storage medium for predicting metal dendrite spacing
CN114637954A (en) * 2022-03-25 2022-06-17 宁夏中欣晶圆半导体科技有限公司 Method for calculating axial distribution of carbon content of crystal bar
CN114637954B (en) * 2022-03-25 2023-02-07 宁夏中欣晶圆半导体科技有限公司 Method for calculating axial distribution of carbon content of crystal bar

Also Published As

Publication number Publication date
CN103729511B (en) 2017-01-25

Similar Documents

Publication Publication Date Title
CN103729511A (en) Method for predicating ingredient segregation degrees in casting process of complex-structure casting
CN105057642B (en) Casting crystalline grain organizes the formation of the analogy method of correlation values
Takaki et al. Phase-field lattice Boltzmann simulations of multiple dendrite growth with motion, collision, and coalescence and subsequent grain growth
CN106407623B (en) Welding process Microstructure evolution analogy method under transient state based on Phase Field
Chen et al. Thermal and solute diffusion in α-Mg dendrite growth of Mg-5wt.% Zn alloy: A phase-field study
CN107309543A (en) A kind of laser weld pools dendritic growth analogy method
Ahmadein et al. Analysis of macrosegregation formation and columnar-to-equiaxed transition during solidification of Al-4 wt.% Cu ingot using a 5-phase model
CN104014768B (en) A kind of method of magnesium alloy arborescent structure numerical simulation
Böttger et al. Relationship between solidification microstructure and hot cracking susceptibility for continuous casting of low-carbon and high-strength low-alloyed steels: a phase-field study
CN105665684B (en) A kind of method of casting crystalline grain tissue values prediction
CN107092754A (en) A kind of alloy grain tissue values Forecasting Methodology
Liu et al. Advances in multi-scale modeling of solidification and casting processes
CN110245449A (en) A kind of magnesium alloy cast part uneven components Numerical Predicting Method
CN103521731B (en) Determine the method for tundish exit concentration under the continuous pouring of steel ingot many bags Variable Composition
CN107243611B (en) A kind of bloom continuous casting solidification end list roller depressing position determines method
WU et al. Modeling of equiaxed and columnar dendritic growth of magnesium alloy
Cho et al. Modeling of microstructural evolution in squeeze casting of an Al-4.5 mass% Cu alloy
CN114139466A (en) Method for calculating permeability of alloy solidification mushy zone
Zhu et al. Crack formation of camshaft castings: hot tearing susceptibility and root causes
Wang et al. Phase-field simulations of dendrite morphologies and selected evolution of primary α-Mg phases during the solidification of Mg-rich Mg–Al-based alloys
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
Peng et al. Analysis on fluid permeability of dendritic mushy zone during peritectic solidification in a temperature gradient
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

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

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.

TR01 Transfer of patent right