CN102236726A - Method for predicting shrinkage porosity in molten metal solidification process and continuous shrinkage porosity prediction method - Google Patents
Method for predicting shrinkage porosity in molten metal solidification process and continuous shrinkage porosity prediction method Download PDFInfo
- Publication number
- CN102236726A CN102236726A CN2010101603050A CN201010160305A CN102236726A CN 102236726 A CN102236726 A CN 102236726A CN 2010101603050 A CN2010101603050 A CN 2010101603050A CN 201010160305 A CN201010160305 A CN 201010160305A CN 102236726 A CN102236726 A CN 102236726A
- Authority
- CN
- China
- Prior art keywords
- active cell
- unit
- temperature
- molten metal
- shrinkage porosite
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
The invention discloses a method for predicting shrinkage porosity in a molten metal solidification process and a continuous shrinkage porosity prediction method. The method for predicting the shrinkage porosity in the molten metal solidification process comprises the following steps of: (1) acquiring temperature field data of molten metal; (2) selecting a predicted unit as a current unit, and acquiring pressure P, viscosity gradient Gv, equivalent cooling rate R1, feeding distance d and critical value C of the current unit in combination with the temperature field data; and (3) judging according to the acquired P, Gv, R1, d and C values based on a criterion, and if the criterion is met, judging that the current unit is a shrinkage porosity unit, wherein e is a natural napierian base. According to the method, the prediction of the shrinkage porosity in the molten metal solidification process has high accuracy. The method disclosed by the invention can be used for predicting the solidification processes of various crystals or amorphous alloy melts and has a wide application range.
Description
Technical field
The present invention relates to method and the shrinkage porosite continuous prediction method of prediction shrinkage porosite in a kind of molten metal process of setting.
Background technology
The process of setting of foundry goods also is accompanied by the motion of heat transfer, mass transfer except the solid-liquid interface migration, many defectives of foundry goods are all closely related with process of setting.The computer simulation of casting solidification process can this change procedure of pictute, for analyzing casting flaw, prediction casting quality, optimizing casting technique positive effect is arranged.
The application of computing machine in the simulation casting process at present mainly concentrates on aspect following 4: fill type solidification simulation, Shrinkage Prediction, process of setting stress simulation and the simulation of process of setting microstructure.Wherein, Shrinkage Prediction is even more important.
For generation and the degree (shrinkage porosite shape, size) thereof of predicting the foundry goods shrinkage porosite, it is necessary that the mechanism that the consideration shrinkage porosite generates is solidified parsing, but tight parsing is difficulty very, in fact often adopts following several simple and easy methods to predict their generation: 1, equitemperature, etc. the solid rate curve method.(1) solidus temperature method judges that the condition that shrinkage porosite produces is: occur in g
sIn=1 the closed loop; (2) critical solid rate method judges that the condition that shrinkage porosite produces is to occur in g
s=g
ScThe closed loop in.But above method is difficult to when not existing in the closed loop judge.2, temperature gradient method judges that shrinkage porosite generation condition is g
s=1 or g
s=g
ScThe time maximum temperature G≤critical value.This method is easier, but since the critical temperature gradient with shape, cooldown rate and different, so precision of prediction is very low.3, solid rate gradient method.Judge that shrinkage porosite generation condition is g
s=1 or g
s=g
ScThe time maximum solid rate≤critical value.This method relatively is fit to the shrinkage porosite prediction of eutectic alloy process of setting; When the pass of solid rate and temperature is 1 pair 1 and thermograde be equal to.
The criterion of main prediction shrinkage porosite of using is at present
Wherein G is thermograde (K/cm), and R is cooldown rate (K/min), and C is a critical value.But this method is lower to alloy graining process shrinkage porosite accuracy of predicting.
Summary of the invention
In order to overcome in the prior art to the lower problem of alloy graining process shrinkage porosite accuracy of predicting, the invention provides method and the shrinkage porosite continuous prediction method of prediction shrinkage porosite in a kind of molten metal process of setting, this method is to alloy liquid process of setting shrinkage porosite accuracy of predicting height, and applied range.
The method of prediction shrinkage porosite in the molten metal process of setting disclosed by the invention comprises:
(1), obtains the temperature field data of molten metal;
(2), to choose the unit of predicting be active cell, in conjunction with described temperature field data, obtains pressure P, the viscosity gradient G of active cell
v, cooldown rate R of equal value
1And feeding distance d and critical value C;
(3), according to the P, the G that obtain
v, R
1, d and C value, according to criterion:
Judge,, can judge that then active cell is the shrinkage porosite unit if satisfy above-mentioned criterion; Wherein, at the bottom of e is natural logarithm.
The invention also discloses a kind of shrinkage porosite continuous prediction method in addition, comprising:
(a), calculate the temperature field of molten metal;
(b), determine according to the temperature field whether molten metal is about to solidify, if determine it is not to be about to solidify, then return step (a) and continue the accounting temperature field; If determine to be about to solidify, then judge according to preceding method.
By said method, to shrinkage porosite accuracy of predicting height in the molten metal process of setting.And this method can be predicted applied range to the process of setting of various crystal or non-crystaline amorphous metal liquation.
Embodiment
In order to make technical matters solved by the invention, technical scheme and beneficial effect clearer,, the present invention is further elaborated below in conjunction with drawings and Examples.Should be appreciated that specific embodiment described herein only in order to explanation the present invention, and be not used in qualification the present invention.
The method of prediction shrinkage porosite in the molten metal process of setting disclosed by the invention comprises:
(1), obtains the temperature field data of molten metal;
(2), to choose the unit of predicting be active cell, in conjunction with described temperature field data, obtains pressure P, the viscosity gradient G of active cell
v, cooldown rate R of equal value
1And feeding distance d and critical value C;
(3), according to the P, the G that obtain
v, R
1, d and C value, according to criterion:
Judge,, can judge that then active cell is the shrinkage porosite unit if satisfy above-mentioned criterion; Wherein, at the bottom of e is natural logarithm.Be e in the present invention value be 2.71828.
According to said method disclosed by the invention, compare with existing shrinkage porosity criterion, the inventor finds that the pressure P of having introduced another molten metal can improve the accuracy of judging as the criterion influence factor.Simultaneously, the present inventor has proposed with viscosity gradient G according to the analysis-by-synthesis of the various factors in the molten metal process of setting
vAs a big influence factor of criterion, and creationary in conjunction with feeding distance d, increase new parameter e
dThereby, formed above-mentioned criterion
The inventor finds that in the molten metal process of setting, the feeding distance of metal also has considerable influence to its process of setting.As the particular community of metal, the feeding distance of metal is big more, obtains feeding easily more when producing shrinkage defect in the molten metal process of setting.For existing criterion, do not consider feeding distance, can be equivalent to that feeding distance is 0 in the criterion disclosed by the invention, the feeding distance that is equivalent to give tacit consent to metal this moment is 0, is not inconsistent with actual conditions.In the present invention, by increasing new parameter e
dConsider feeding distance, the feeding distance of metal is big more, and the predicted value that calculates by above-mentioned criterion is also big more, and it is big more to meet feeding distance, the easy more practical situation that obtains feeding.Improved precision of prediction greatly.Introducing pressure P and G
vAnd in conjunction with e
dAfter criterion revised, the precision of prediction that shrinkage porosite in the molten metal process of setting is produced had improved greatly.And, than traditional solid rate gradient criterion, the present invention is by adopting in this area not used viscosity gradient criterion, can realize the shrinkage porosite prediction to the process of setting of the metal that comprises various forms such as crystal, amorphous, and cooldown rate of equal value is not subjected to the influence of metallic crystal form, the scope of the application of expansion; Accuracy of predicting has also obtained further raising.Adopt method disclosed by the invention to after the shrinkage porosite situation is judged in the molten metal process of setting, help the technician smelting technology of metal is adjusted.
In above-mentioned criterion, the unit of pressure P is Pa, viscosity gradient G
vUnit be kg/ (m
2.s), cooldown rate R of equal value
1Unit be K/min, the unit of feeding distance d is mm.
In the present invention, the temperature field data can obtain by existing the whole bag of tricks among the step a, as calculating by the thermal conduction study fundamental equation.
Casting solidification process is actually the unstable heat exchanging process between " casting and mold-environment ".All phenomenons in conducting heat in casting solidification process, have been consisted essentially of: heat conduction, thermal convection, heat radiation.After liquid metal injects casting mold, the heat of liquid metal inside arrives the surface of foundry goods by heat conduction, pass to casting mold by heat radiation and heat conduction then, arrive the outside surface of casting mold again by heat conduction, be dispersed in the environment by heat radiation and convection current at last.Meanwhile the temperature of liquid metal constantly descends in the casting mold die cavity, and non-uniform temperature.This unevenness has caused the thermal convection of the liquid metal in the casting mold die cavity.
Respectively heat conduction, thermal convection and heat radiation are described below.
1, heat conduction heat exchange
Heat conduction is called for short heat conduction, and it belongs to transmission of heat by contact, be continuous medium under the situation that does not have each several part relative displacement between the material, the heat transferred that relies on the thermal vibration of microscopic particles such as molecule, atom and free electron to carry out.In closely knit opaque solid, can only rely on heat-conducting mode to transmit heat.
For determining of multidimensional temperature field, must analyze the micro unit in the heat conductor based on energy conservation and Fourier's law, draw the three-dimensional heat conducting differential equation of expression heat conduction phenomenon basic law:
In the formula:
Q
3-heat conduction thermal change value, unit is: J;
The density of ρ-object, unit is: Kg/m
3
λ-coefficient of heat conductivity, unit is: W/mK;
C-specific heat, unit is: J/m
3
The T-temperature, unit is: K;
The t-time, unit is: s;
The coordinate figure of x, y, z-arbitrfary point.
Choose active cell, promptly obtain x, y, the z value of active cell,, can inquire its density p, coefficient of heat conductivity λ and specific heat c then according to concrete metal species.According to the temperature field data in a last moment, obtained the temperature data T of a last moment active cell again, can calculate heat conduction thermal change value by above-mentioned equation.
Wherein, because the temperature of initial time is prior setting, promptly sets molten metal and reach beginning process of setting simulation after certain temperature, this temperature is the temperature of initial time, so be known for the temperature field data of initial time.According to these known initial time temperature field data, can calculate from initial time to next heat conduction thermal change value constantly.
2, thermal convection heat exchange
Convection current is meant between the fluid each several part relative displacement takes place, the caused heat transferred mode of the mutual blending of cold flow object.Convection current only can occur in the fluid, and must be accompanied by heat conduction.
Convection heat transfer is formal what calculate, adopts the formula that newton proposed, that is:
Q
1＝αA(T
f-T
w)
In the formula:
Q
1-thermal convection thermal change value, unit is: J;
Heat under α-unit temperature difference on unit area, unit is: W/m
2K;
The A-heat interchanging area, unit is: m
2
T
fEnvironment temperature in the-thermal convection heat transfer process, unit is: K;
T
wThe temperature of-wall, unit is: K.
In the above-mentioned equation, α can look into by the related tool book and get for certain specific metal material.
T
fActual temperature according to environment in the thermal convection process is set.
T
wThe temperature field data of passing through are directly obtained.
A is for the particular model difference, as all different with cube for right cylinder, but long-pending the obtaining of mold surface of the active cell that all can select for use by calculating.
Can calculate thermal convection thermal change value Q by above-mentioned equation
1
3, heat radiation heat exchange
Directly can not transmit radiation energy between the two articles of contact, at this moment need not to exist any medium between the object yet.Usually the process that object is sent radiant rays and propagation thereof calls radiation, and object calls heat radiation to the interior process that can be converted into external emitted radiation line and propagation thereof of itself, and promptly the object phenomenon of sending radiation energy because of the reason of heat is called heat radiation.Be called the heat radiation heat exchange by sending and accept the heat exchange that radiation energy carries out between the object.
Casting mold is to the heat of surrounding environment heat radiation heat exchange:
In the formula:
Q
2-heat radiation thermal change value, unit is: J;
The blackness of ε-grey body, 0-1;
σ
b-Si Difen-Boltzmann constant is 5.67 * 10
-8(W/m
2K
4);
A
1The area of dissipation of-object, unit is: m
2
T
1The temperature of object in the-heat radiation process, unit is: K;
T
2The temperature of environment in the-heat radiation process, unit is: K.
In the above-mentioned equation, ε can look into by the related tool book and get for certain specific metal material.
T
1Be the temperature of active cell.
T
2The temperature of environment is set and is obtained in the heat radiation process.
A
1For the particular model difference,, obtain but the mold surface of the active cell that all can select for use by calculating is long-pending as all different with cube for right cylinder.
Can calculate heat radiation thermal change value Q by above-mentioned equation
2
Method by thermal conduction study Equation for Calculating temperature field data comprises:
(11), choose arbitrary unit, based on last one constantly temperature field data, draw on this unit constantly a temperature data; Calculate this unit by heat conduction, thermal convection, heat radiation then and be carved into the total amount of heat changing value of current time,, adopt this unit of temperature variation Equation for Calculating to be carved into the temperature change value of current time from last a period of time by this total amount of heat changing value from last a period of time; Again based on last one constantly temperature data and be carved into the temperature change value of current time from last a period of time, obtain the temperature data of this unit current time;
Described temperature variation equation is:
Wherein: Q is the total amount of heat changing value, is aforementioned Q
1, Q
2, Q
3Sum, unit is J, and C is the specific heat of metal, and unit is kJ/kg ℃, and ρ is the density of metal, unit is kg/m
3, V is the volume of active cell, unit is m
3
(12), repeating step (11), to calculating all unit of molten metal, obtain the temperature field data of molten metal by the temperature variation equation.
Wherein, heat conduction thermal change value, thermal convection thermal change value, heat radiation thermal change value can be calculated by preceding method.Heat conduction thermal change value, thermal convection thermal change value, heat radiation thermal change value sum are the total amount of heat changing value.
Those skilled in the art can obtain the said temperature field data by the calculating of above equation.
The pressure P of active cell can obtain by existing the whole bag of tricks, as the molten metal static pressure that can directly adopt active cell as pressure P.After calculation of pressure finishes, remain to follow-up shrinkage porosity criterion and use.
Wherein, described molten metal static pressure can be calculated in the following way:
P＝h×ρ×g
In the formula:
The molten metal static pressure that the P-active cell is subjected to, unit: Pa;
The h-active cell is apart from the distance at the top, molten bath at place, unit: m;
The density of ρ-molten metal, unit: kg/m
3
G-acceleration of gravity, unit: m/s
2
In the said method, the molten bath may be defined as in the molten metal process of setting, the zone that molten metal can be interconnected.Promptly a molten bath can comprise a plurality of unit.The division in molten bath can adopt existing many molten baths partitioning technology to carry out.
In method disclosed by the invention, viscosity gradient G
vAcquisition methods be:
Choose 26 adjacent cells adjacent with active cell, described 26 adjacent cells are 6 unit adjacent with the active cell face, 12 unit and with active cell summit adjacent 8 unit adjacent with the active cell rib, the distance of active cell and each adjacent cells is the distance between active cell central point and the adjacent cells central point;
Obtain the viscosity number of active cell and 26 adjacent cells;
The difference of viscosity number of getting active cell and one of them adjacent cells obtains a quotient divided by the distance of active cell and this adjacent cells, repeats this step, and getting maximum quotient then is viscosity gradient.
In the said method, the method of obtaining the viscosity number of active cell and 26 adjacent cells is: in conjunction with the temperature field data, obtain the temperature data of active cell and 26 adjacent cells respectively, according to the temperature data of active cell and 26 adjacent cells, obtain the viscosity number of active cell and 26 adjacent cells respectively.
For every kind of fixing molten metal of forming, when other conditions were constant, its viscosity under fixing temperature was fixed, and promptly viscosity has fixing curve with respect to temperature.So, obtain the temperature data of active cell and 26 adjacent cells after, can directly read the viscosity of active cell and 26 adjacent cells respectively.As can in conjunction with the temperature-viscograph of this metal material, directly reading the viscosity number of active cell and 26 adjacent cells according to the temperature data of active cell and 26 adjacent cells.Temperature-the viscograph of above-mentioned metal material is for recording by prior art.
In the present invention, cooldown rate R of equal value
1Acquisition methods be: the latent heat that active cell is discharged in a time step is converted into the equivalence reduction value of active cell temperature in this time step, get the actual reduction value sum of this equivalence reduction value and this unit actual temperature, divided by described time step, obtain the cooldown rate R of equal value of active cell then
1Above-mentioned time step is the mistiming of current time and previous moment.
Solidify the key criterion that whether produces shrinkage porosite as weighing molten metal, critical value C is most important.For different alloys, critical value C is also different.Usually, this critical value C is known.Simultaneously, this critical value C also can obtain voluntarily.The acquisition methods of critical value C can be the whole bag of tricks of the prior art, mainly obtains with the method that simulation combines by actual die casting as criterion critical value C.At certain alloy material,, carry out actual die casting based on conditions such as given model and extrusion processes.Observe position of inner shrinkage porosite etc. by foundry goods being carried out subdivision.Then, constantly this critical value is adjusted, set the predicted critical value, adopt basic criterion
Carry out sunykatuib analysis and obtain to comprise and the analog result of shrinkage porosite position the position of analog result with foundry goods inside shrinkage porosite contrasted,, then adjust the predicted critical value, simulate again if the position of shrinkage porosite is inequality in analog result and the actual foundry goods; If analog result is identical with the position of shrinkage porosite in the actual foundry goods, then this predicted critical value is critical value C.For verifying the correctness of this critical value, can analyze according to different models and extrusion process.In the prior art, critical value C can directly obtain by the correlation technique dictionary enquiry.In the present invention, the acquisition methods of critical value C is: combine with simulation by actual die casting, and revise and obtain.
According to the present invention, after obtaining critical value C, the acquisition methods of described feeding distance d is: at fixing metal material, based on fixing model and extrusion process condition, carry out actual die casting, obtain foundry goods; Foundry goods is carried out the position that subdivision is determined the inner shrinkage porosite of foundry goods;
Then,, set the prediction feeding distance, adopt criterion according to the critical value C that obtains
Carry out sunykatuib analysis, obtain to comprise the analog result of shrinkage porosite position, the position of analog result with the inner shrinkage porosite of foundry goods contrasted, if the position of shrinkage porosite is inequality in analog result and the actual foundry goods, then feeding distance is predicted in adjustment, simulates again; If analog result is identical with the position of shrinkage porosite in the actual foundry goods, then this prediction feeding distance is feeding distance d.
Obtain above-mentioned P, G
v, R
1, after d and the C value, according to criterion:
Judge,, can judge that then active cell is the shrinkage porosite unit if satisfy above-mentioned criterion.
For the ease of using in practice, in the present invention, a kind of shrinkage porosite continuous prediction method is also disclosed, comprising:
(a), calculate the temperature field of molten metal;
(b), determine according to the temperature field whether molten metal is about to solidify, if determine it is not to be about to solidify, then return step (a) and continue the accounting temperature field; If determine to be about to solidify, then judge according to the method for prediction shrinkage porosite in the above-mentioned molten metal process of setting.
In said method, the temperature field of having calculated molten metal in step (a) then in step (b), can directly adopt this temperature field as required temperature field data.
In the said method, determine according to the temperature field whether molten metal is about to be solidified as method well known in the art.Molten metal is solidified as phase transition process, and those skilled in the art can judge whether to be about to solidify according to the temperature field data of the molten metal that obtains.
Equally, the temperature field described in the step (a) can obtain by the thermal conduction study Equation for Calculating.Among the present invention, the acquisition of various data and calculating all can be finished by computer installation of the prior art.
According to the present invention, be example with the non-crystaline amorphous metal, its process of setting is carried out shrinkage porosite predict continuously.Concrete grammar is: after Metal Melting finishes, constantly, determine according to the temperature field whether molten metal is about to solidify by the temperature field of thermal conduction study Equation for Calculating molten metal, if not, then return step (a) and continue the accounting temperature field; If determine to be about to solidify, choosing the unit of predicting is active cell, selects the pressure P of the molten metal static pressure of active cell as active cell for use.
Choose 26 adjacent cells adjacent with active cell, described 26 adjacent cells are 6 unit adjacent with the active cell face, 12 unit and with active cell summit adjacent 8 unit adjacent with the active cell rib, the distance of active cell and each adjacent cells is the distance between active cell central point and the adjacent cells central point;
Obtain the viscosity number of active cell and 26 adjacent cells;
The difference of viscosity number of getting active cell and certain adjacent cells obtains a quotient divided by the distance of active cell and this adjacent cells, repeats this step, and getting maximum quotient then is viscosity gradient G
v
Then, the latent heat that active cell is discharged in a time step is converted into the equivalence reduction value of active cell temperature in this time step, get the actual reduction value sum of this equivalence reduction value and this unit actual temperature, divided by described time step, obtain the cooldown rate R1 of equal value of active cell then.
At last, according to the P, the G that obtain
v, R
1, d, and in conjunction with known critical value C, according to criterion:
Judge,, can judge that then active cell is the shrinkage porosite unit if satisfy above-mentioned criterion.
Adopt the method for prediction shrinkage porosite in the molten metal process of setting disclosed by the invention, can realize the accurately predicting of shrinkage porosite situation in the molten metal process of setting.And the applied range of this method can be predicted shrinkage porosite situation in the molten metal process of setting of each metalloid material such as eutectic or amorphous
The above only is preferred embodiment of the present invention, not in order to restriction the present invention, all any modifications of being done within the spirit and principles in the present invention, is equal to and replaces and improvement etc., all should be included within protection scope of the present invention.
Claims (10)
1. the method for prediction shrinkage porosite in the molten metal process of setting comprises:
(1), obtains the temperature field data of molten metal;
(2), to choose the unit of predicting be active cell, in conjunction with described temperature field data, obtains pressure P, the viscosity gradient G of active cell
v, cooldown rate R of equal value
1And feeding distance d and critical value C;
(3), according to the P, the G that obtain
v, R
1, d and C value, according to criterion:
Judge,, can judge that then active cell is the shrinkage porosite unit if satisfy above-mentioned criterion; Wherein, at the bottom of e is natural logarithm.
2. method according to claim 1, wherein, the temperature field data obtain by the thermal conduction study Equation for Calculating in the described step (1).
3. method according to claim 2, wherein, described method by thermal conduction study Equation for Calculating temperature field data comprises:
(11), choose arbitrary unit, based on last one constantly temperature field data, draw on this unit constantly a temperature data; Calculate this unit by heat conduction, thermal convection, heat radiation then and be carved into the total amount of heat changing value of current time,, adopt this unit of temperature variation Equation for Calculating to be carved into the temperature change value of current time from last a period of time by this total amount of heat changing value from last a period of time; Again based on last one constantly temperature data and be carved into the temperature change value of current time from last a period of time, obtain the temperature data of this unit current time;
Described temperature variation equation is:
Wherein: Q is the total amount of heat changing value, and unit is J, and C is the specific heat of metal, and unit is kJ/kg ℃, and ρ is the density of metal, and unit is kg/m
3, V is the volume of active cell, unit is m
3
(12), repeating step (11), to calculating all unit of molten metal, obtain the temperature field data of molten metal by the temperature variation equation.
4. method according to claim 1, wherein, described pressure is the molten metal static pressure of active cell.
5. method according to claim 1, wherein, described viscosity gradient G
vAcquisition methods be:
Choose 26 adjacent cells adjacent with active cell, described 26 adjacent cells are 6 unit adjacent with the active cell face, 12 unit and with active cell summit adjacent 8 unit adjacent with the active cell rib, the distance of active cell and each adjacent cells is the distance between active cell central point and the adjacent cells central point;
Obtain the viscosity number of active cell and 26 adjacent cells;
The difference of viscosity number of getting active cell and one of them adjacent cells obtains a quotient divided by the distance of active cell and this adjacent cells, repeats this step, and getting maximum quotient then is viscosity gradient G
v
6. method according to claim 6, wherein, the method of obtaining the viscosity number of active cell and 26 adjacent cells is: in conjunction with the temperature field data, obtain the temperature data of active cell and 26 adjacent cells respectively, according to the temperature data of active cell and 26 adjacent cells, obtain the viscosity number of active cell and 26 adjacent cells respectively.
7. method according to claim 1, wherein, the cooldown rate R of equal value of described active cell
1Acquisition methods be: the latent heat that active cell is discharged in a time step is converted into the equivalence reduction value of active cell temperature in this time step, get the actual reduction value sum of this equivalence reduction value and this unit actual temperature, divided by described time step, obtain the cooldown rate R of equal value of active cell then
1
8. method according to claim 1, wherein, the acquisition methods of described critical value C is:
At fixing metal material, based on fixing model and extrusion process condition, carry out actual die casting, obtain foundry goods; Foundry goods is carried out the position that subdivision is determined the inner shrinkage porosite of foundry goods;
Then, set the predicted critical value, adopt basic criterion
Carry out sunykatuib analysis, obtain to comprise the analog result of shrinkage porosite position, the position of analog result with foundry goods inside shrinkage porosite contrasted,, then adjust the predicted critical value, simulate again if the position of shrinkage porosite is inequality in analog result and the actual foundry goods; If analog result is identical with the position of shrinkage porosite in the actual foundry goods, then this predicted critical value is critical value C.
9. method according to claim 8, wherein, the acquisition methods of described feeding distance d is: at fixing metal material, based on fixing model and extrusion process condition, carry out actual die casting, obtain foundry goods; Foundry goods is carried out the position that subdivision is determined the inner shrinkage porosite of foundry goods;
Then,, set the prediction feeding distance, adopt criterion according to the critical value C that obtains
Carry out sunykatuib analysis, obtain to comprise the analog result of shrinkage porosite position, the position of analog result with the inner shrinkage porosite of foundry goods contrasted, if the position of shrinkage porosite is inequality in analog result and the actual foundry goods, then feeding distance is predicted in adjustment, simulates again; If analog result is identical with the position of shrinkage porosite in the actual foundry goods, then this prediction feeding distance is feeding distance d.
10. shrinkage porosite continuous prediction method comprises:
(a), calculate the temperature field of molten metal;
(b), determine according to the temperature field whether molten metal is about to solidify, if determine it is not to be about to solidify, then return step (a) and continue the accounting temperature field; If determine to be about to solidify, then judge according to any described method among the claim 1-9.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN 201010160305 CN102236726B (en) | 2010-04-23 | 2010-04-23 | Method for predicting shrinkage porosity in molten metal solidification process and continuous shrinkage porosity prediction method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN 201010160305 CN102236726B (en) | 2010-04-23 | 2010-04-23 | Method for predicting shrinkage porosity in molten metal solidification process and continuous shrinkage porosity prediction method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102236726A true CN102236726A (en) | 2011-11-09 |
| CN102236726B CN102236726B (en) | 2013-01-30 |
Family
ID=44887372
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN 201010160305 Expired - Fee Related CN102236726B (en) | 2010-04-23 | 2010-04-23 | Method for predicting shrinkage porosity in molten metal solidification process and continuous shrinkage porosity prediction method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102236726B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107844852A (en) * | 2017-10-20 | 2018-03-27 | 太原理工大学 | A kind of shrinkage defect Forecasting Methodology for simulating steel-casting sand casting process |
| CN112115634A (en) * | 2020-09-21 | 2020-12-22 | 哈尔滨理工大学 | Three-dimensional numerical prediction method for grain structure in unidirectional solidification process of molten metal |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1996319A (en) * | 2006-12-22 | 2007-07-11 | 山东大学 | Method for prediction of shrinkage porosity in casting solidification simulation |
| CN101239385A (en) * | 2008-03-18 | 2008-08-13 | 上海大学 | Casting method capable of effectively improving quality of aluminum alloy casting parts |
| CN101567019A (en) * | 2009-05-08 | 2009-10-28 | 江苏大学 | Computer simulation method for casting property of particle-reinforced aluminum matrix composite |
| CN101602102A (en) * | 2009-06-17 | 2009-12-16 | 北京航空航天大学 | The solidification process control method of using small external temperature gradient to eliminate shrinkage cavities and porosity in casting |
-
2010
- 2010-04-23 CN CN 201010160305 patent/CN102236726B/en not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1996319A (en) * | 2006-12-22 | 2007-07-11 | 山东大学 | Method for prediction of shrinkage porosity in casting solidification simulation |
| CN101239385A (en) * | 2008-03-18 | 2008-08-13 | 上海大学 | Casting method capable of effectively improving quality of aluminum alloy casting parts |
| CN101567019A (en) * | 2009-05-08 | 2009-10-28 | 江苏大学 | Computer simulation method for casting property of particle-reinforced aluminum matrix composite |
| CN101602102A (en) * | 2009-06-17 | 2009-12-16 | 北京航空航天大学 | The solidification process control method of using small external temperature gradient to eliminate shrinkage cavities and porosity in casting |
Non-Patent Citations (1)
| Title |
|---|
| 郝良品 等: "定向凝固过程的数值模拟进展", 《铸造》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107844852A (en) * | 2017-10-20 | 2018-03-27 | 太原理工大学 | A kind of shrinkage defect Forecasting Methodology for simulating steel-casting sand casting process |
| CN112115634A (en) * | 2020-09-21 | 2020-12-22 | 哈尔滨理工大学 | Three-dimensional numerical prediction method for grain structure in unidirectional solidification process of molten metal |
| CN112115634B (en) * | 2020-09-21 | 2021-04-06 | 哈尔滨理工大学 | Three-dimensional numerical prediction method for grain structure in unidirectional solidification process of molten metal |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102236726B (en) | 2013-01-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Beckermann et al. | Equiaxed dendritic solidification with convection: Part III. Comparisons with NH 4 Cl-H 2 O experiments | |
| Shepel et al. | Numerical simulation of filling and solidification of permanent mold castings | |
| Lewis et al. | Finite element simulation of metal casting | |
| Liu et al. | Simulation of EPS foam decomposition in the lost foam casting process | |
| Zappulla et al. | Multiphysics modeling of continuous casting of stainless steel | |
| Bounds et al. | A computational model for defect prediction in shape castings based on the interaction of free surface flow, heat transfer, and solidification phenomena | |
| Lin et al. | Cooling of a Pr< 1 fluid in a rectangular container | |
| Duan et al. | Experimental measurements for numerical simulation of macrosegregation in a 36-ton steel ingot | |
| Chen et al. | Progress in numerical simulation of casting process | |
| CN102294466B (en) | Method for predicting shrinkage in metal liquid solidification process | |
| Mirbagheri et al. | Simulation of melt flow in coated mould cavity in the casting process | |
| Anglada et al. | Prediction and validation of shape distortions in the simulation of high pressure die casting | |
| Ravindran et al. | Finite element modelling of solidification effects in mould filling | |
| Nzebuka et al. | Numerical modeling of the effect of nozzle diameter and heat flux on the polymer flow in fused filament fabrication | |
| CN102236726B (en) | Method for predicting shrinkage porosity in molten metal solidification process and continuous shrinkage porosity prediction method | |
| US8712750B2 (en) | Molten alloy solidification analyzing method and solidification analyzing program for performing the same | |
| Miłkowska-Piszczek et al. | A comparison of models describing heat transfer in the primary cooling zone of a continuous casting machine | |
| CN106944607A (en) | A kind of inoculant alloy grain structure Numerical Predicting Method | |
| Egole et al. | Micro-macro model for the transient heat and fluid transport in solidification structure evolution during static casting processes | |
| Liu et al. | Stochastic modeling of columnar-to-equiaxed transition in Ti–(45–48 at%) Al alloy ingots | |
| Mahmoudi et al. | An experimental and numerical study on the modelling of fluid flow, heat transfer and solidification in a copper continuous strip casting process | |
| Monde et al. | Shrinkage induced flow during directional solidification of pure substance in a bottom cooled cavity: A study on flow reversal phenomena | |
| CN102072914B (en) | Method for predicting shrinkage in solidifying process of molten metal | |
| Mirbagheri et al. | Modelling of foam degradation in lost foam casting process | |
| Katzarov et al. | Porosity formation in axi-symmetric castings produced by counter-pressure casting method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130130 Termination date: 20170423 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |