CN101871903A - Method for determining interface heat exchange coefficient of large-sized steel ingot - Google Patents

Method for determining interface heat exchange coefficient of large-sized steel ingot Download PDF

Info

Publication number
CN101871903A
CN101871903A CN201010195121A CN201010195121A CN101871903A CN 101871903 A CN101871903 A CN 101871903A CN 201010195121 A CN201010195121 A CN 201010195121A CN 201010195121 A CN201010195121 A CN 201010195121A CN 101871903 A CN101871903 A CN 101871903A
Authority
CN
China
Prior art keywords
ingot
heat exchange
exchange coefficient
steel ingot
interface
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
CN201010195121A
Other languages
Chinese (zh)
Other versions
CN101871903B (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.)
Tsinghua University
Original Assignee
Tsinghua 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 Tsinghua University filed Critical Tsinghua University
Priority to CN2010101951218A priority Critical patent/CN101871903B/en
Publication of CN101871903A publication Critical patent/CN101871903A/en
Application granted granted Critical
Publication of CN101871903B publication Critical patent/CN101871903B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Continuous Casting (AREA)
  • Investigating And Analyzing Materials By Characteristic Methods (AREA)

Abstract

The invention relates to a method for determining an interface heat exchange coefficient of a large-sized steel ingot. The method comprises the following steps of: 1) simplifying a contact interface between the steel ingot and an ingot mould into a plurality of test subsections, determining a plurality of test points in each subsection, setting a sensor on each test point, connecting the sensors to a data acquisition system together through respective leads, and connecting an interface heat exchange coefficient inverse calculation system at the output end of the data acquisition system; 2) when molten steel is solidified to form the steel ingot, acquiring temperature or heat flux information during molten steel solidification and transmitting the acquired signals to the data acquisition system through the sensors; 3) after the molten steel is fully solidified into the steel ingot, reading the data of each test point through the data acquisition system, acquiring the actually measured temperature or heat flux information of each test point along with time change, and importing the data and the information into the interface heat exchange coefficient inverse calculation system between the steel ingot and the ingot mould; and 4) solving the interface heat exchange coefficient between the steel ingot and the ingot mould by inverse calculation of the interface heat exchange coefficient inverse calculation system.

Description

A kind of method of definite interface heat exchange coefficient of large-sized steel ingot
Technical field
The present invention relates to a kind of method of definite coefficient of heat transfer, particularly about a kind of method of definite interface heat exchange coefficient of large-sized steel ingot.
Background technology
At present, the maximization along with great equipments of industrial circle such as electric power, petrochemical industry, metallurgy, boats and ships has also proposed more and more higher requirement to large-scale steel ingot.For example, the manufacturing of million kilowatt nuclear power unit conventional island low pressure integral rotor needs 600 tonnes of steel ingots.In the trial-manufacturing process of large-scale steel ingot, numerical simulation technology is being brought into play more and more important effect.Interfacial heat exchange behavior between steel ingot and ingot mould has determined the solidification mode and the setting time of steel ingot, and has influenced shrinkage cavity shrinkage porosity, segregation and tissue and the performance of steel ingot to a great extent.And the interface coefficient of heat transfer is one of most important parameter of restriction analog result accuracy between steel ingot and ingot mould.
Go into initial stage of ingot mould at pouring molten steel, because molten steel flow and static pressure effect, molten steel is closely to contact with the ingot mould inwall, and the interface coefficient of heat transfer is bigger between molten steel and ingot mould at this moment.And along with the propelling of process of setting, because ingot mould is to solidifying the chilling action of molten steel, molten steel forms along the ingot mould inside surface and solidifies shell, and in cooling procedure subsequently continuous thickening and produce volumetric contraction; Meanwhile, the expanded by heating of ingot mould also can move by generation type wall shift.Therefore, between steel ingot and ingot mould, form air gap usually.The formation time of air gap is with spatial variations, and size of gaps also changes with the space in time.After air gap formed, the heat transmission at interface then became heat conduction, convection current and three sums of radiation by simple heat conduction.More than these complicated factors, cause that the interface coefficient of heat transfer is difficult to determine between steel ingot and ingot mould.And in actual computation, normally adopt simplification measures such as a conjecture value or experimental formula that the interface heat exchange coefficient between steel ingot and ingot mould is determined, but the error calculated of these conjecture values or experimental formula is very big.
Summary of the invention
At the problems referred to above, the purpose of this invention is to provide a kind of simple to operately, can determine the method for interface heat exchange coefficient of large-sized steel ingot fast and accurately.
For achieving the above object, the present invention takes following technical scheme: a kind of method of definite interface heat exchange coefficient of large-sized steel ingot, it may further comprise the steps: 1) contact interface between steel ingot and the ingot mould is reduced to some test segmentations, in each segmentation, determine several test points, and a sensor is set on each test point, the lead of each sensor by separately is connected to a data acquisition system (DAS) jointly, and the output terminal of data acquisition system (DAS) connects an interface heat exchange coefficient inverse system; 2) become in the process of steel ingot at molten steel solidification, temperature or hot-fluid information in each sensor acquisition molten steel solidification process, and the information data that collects is sent in the data acquisition system (DAS); 3) after molten steel is frozen into steel ingot fully, read the data of each test point by data acquisition system (DAS), obtain the time dependent temperature or the hot-fluid information of each test point actual measurement, and import between steel ingot and ingot mould in the coefficient of heat transfer inverse system of interface; 4) interface heat exchange coefficient inverse system estimates the interface heat exchange coefficient of certain time period according to the data that step 3) imports; Carry out simulation of Temperature in the process of setting according to the interface heat exchange coefficient of estimating, to obtain the accounting temperature of test point; Relatively whether accounting temperature and observed temperature coincide, if accounting temperature and observed temperature are misfitted, then find the solution by inverse, proofread and correct interface heat exchange coefficient, and then return the simulation of solidification processing temperature field, up to the result of calculation that obtains to match with observed temperature, promptly accounting temperature and observed temperature are identical, at this moment obtain the interface heat exchange coefficient of this time period, the record interface heat exchange coefficient; Judge that whether setting time t is greater than complete setting time t MaxIf, t＞t Max, then finish to calculate; Otherwise, the calculating that enters next time period.
The inverse method for solving of described step 4) median surface coefficient of heat transfer inverse system is: the starting condition that 1. certain time period is set; Corresponding first time period, then be: the starting condition of input temp field, promptly t=0 Temperature Distribution constantly is made as the pouring molten steel temperature usually, and the initial value of interface heat exchange coefficient, i.e. estimated value; Definition F (h 1..., h I) be
F ( h 1 , . . . , h I ) = Σ i = 1 I Σ j = 1 J Σ k = 1 K ( Y k , j , i - T k , j , i ) 2 - - - ( 1 )
In the formula, h 1..., h IBe the interface heat exchange coefficient of each segmentation, Y K, j, iAnd T K, j, iRepresent temperature measured value and the calculated value of segmentation i at time period j, position k respectively, I is the segmentation sum, and J is the time period sum, and K is the test point sum of a certain segmentation i correspondence; 2. utilize interface heat exchange coefficient
Figure BSA00000147402200022
, find the solution the temperature field of this time period, i.e. steel ingot three-dimensional temperature field simulation, subscript l represents iterations, l=0,1,2 ..., wherein, l=0 is illustrated in the initial of each time period, from estimated value
Figure BSA00000147402200023
The beginning iteration; 3. at interface heat exchange coefficient The basis on increase an amplitude of variation ζ, generally be taken as 10 -3, utilize this new interface heat exchange coefficient
Figure BSA00000147402200025
, find the solution the temperature field of this time period; 4. according to step 2. and the temperature field of 3. finding the solution, calculate interface heat exchange coefficient
Figure BSA00000147402200026
Increment
Figure BSA00000147402200027
:
δ h i l + 1 = Σ i = 1 I Σ j = 1 J Σ k = 1 K ( Y k , j , i - T k , j , i l ) φ k , j , i l Σ i = 1 I Σ j = 1 J Σ k = 1 K ( φ k , j , i l ) 2 - - - ( 2 )
In the formula,
Figure BSA00000147402200029
The expression sensitivity coefficient, it is defined as:
φ k , j , i l = ∂ T k , j , i ∂ h i | h = h i l - - - ( 3 )
And the interface heat exchange coefficient that obtains next iteration step is:
h i l + 1 = h i l + δ h i l + 1 - - - ( 4 )
If 5.
Figure BSA000001474022000212
Satisfy:
&delta; h i l + 1 h i l + 1 < &epsiv; - - - ( 5 )
Then explanation
Figure BSA00000147402200032
Tally with the actual situation, it also is accurately that correspondence is found the solution the temperature field that obtains, and then, returns the calculating that 2. step enters the next time period; Wherein, ε is a decimal, generally is taken as 10 -4If
Figure BSA00000147402200033
Do not satisfy formula (5), then to interface heat exchange coefficient Carry out assignment again, promptly
Figure BSA00000147402200035
, return step and 2. carry out iterative computation, up to satisfying formula (5); If 6. t＞t Max, then finish to calculate; Otherwise, return the calculating that 2. step enters next time period.
In the described step 1), the ingot mould in the contact interface outside between steel ingot and the ingot mould is reduced to five test segmentations, is respectively: 1. with the horizontal surface of contact of indent on steel ingot ingot tail and ingot mould chassis as segmentation I; 2. with the horizontal surface of contact of the side direction surface of contact on steel ingot ingot tail and ingot mould chassis and steel ingot ingot body and ingot mould mould body as segmentation II; 3. with the vertical surface of contact of steel ingot ingot body and ingot mould mould body as segmentation III; 4. with the vertical surface of contact of steel ingot rising head and ingot mould riser buss as segmentation IV; 5. with the horizontal surface of contact of the exothermic mixture bottom surface of pouring into a mould on steel ingot rising head and the steel ingot as segmentation V.
Each test test point quantity in segmentation is 1～5, is 20～50mm apart from the nearest test point of contact interface between steel ingot and the ingot mould and the distance at interface.
Described sensor is a thermopair.
The present invention is owing to take above technical scheme, and it has the following advantages: what 1, the present invention tested is the ingot mould in the contact interface outside between steel ingot and the ingot mould, tests simple relatively; And by the measurement to the temperature variation of limited unique point on the ingot mould, in conjunction with the temperature field simulation of actual steel ingot and ingot mould, but inverse obtains interface heat exchange coefficient, and is easy to operate, and test accurately.2, the present invention is reduced to five sections with the ingot mould in the contact interface outside between steel ingot and the ingot mould, and the interface heat exchange coefficient of diverse location differs less in every section, can simplify and think equal; And the interface heat exchange coefficient between section and the section differs bigger, and then can carry out areal survey.3, the present invention is provided with some test points in each segmentation, arranges a sensor on the position of each test point, and each sensor is connected to data acquisition system (DAS) jointly; Data importing interface heat exchange coefficient inverse system with data acquisition system (DAS) output, by interface heat exchange coefficient inverse system, need not to measure the formation of interface air gap, can determine interface heat exchange coefficient, and the interface heat exchange coefficient of determining tallies with the actual situation and along with the time changes, and then the formation that has solved air gap is difficult to measure, and because the influence of air gap, interface heat exchange coefficient is along with the time changes, and can't accurately measure the problem of interface heat exchange coefficient.The present invention is skillfully constructed, and the interface heat exchange coefficient inverse method for solving that is proposed can be handled the situation at any multistage interface, can be widely used in Castingother processes such as sand casting, die casting and any interface and exist in the heat transfer system of air gap.
Description of drawings
Fig. 1 is steel ingot of the present invention and ingot mould interface heat exchange coefficient segmentation simplified structure synoptic diagram
Fig. 2 is apparatus of the present invention connection diagrams
Fig. 3 is the schematic flow sheet of finding the solution in the interface heat exchange coefficient inverse of the present invention system
Fig. 4 is the idiographic flow synoptic diagram of inverse solving heat exchange coefficient in the interface heat exchange coefficient inverse of the present invention system
Embodiment
Below in conjunction with drawings and Examples the present invention is described in detail.
The present invention is the coefficient of heat transfer that will determine contact interface between steel ingot 1 and the ingot mould 2.As shown in Figure 1, steel ingot 1 of the present invention comprises ingot tail 11, ingot body 12 and rising head 13.Ingot mould 2 of the present invention comprises chassis 21, mould body 22, riser buss 23 and exothermic mixture 24, and exothermic mixture 24 is cast in rising head 13 end faces.Because the inboard of contact interface is molten steel or solidifies shell that difficulty of test is big, the testing cost height between steel ingot 1 and the ingot mould 2; And the outside of contact interface is an ingot mould 2 between steel ingot 1 and the ingot mould 2, tests simply relatively, and therefore, the present invention determines the coefficient of heat transfer of contact interface between steel ingot 1 and the ingot mould 2 by the temperature variation of test contact interface outside ingot mould 2.
The present invention is based on following thought:
The present invention is divided into some test segmentations with the ingot mould 2 in the contact interface outside between steel ingot 1 and the ingot mould 2, each segmentation is provided with some test points, by the temperature variation and articulated system (steel ingot 1 and the ingot mould 2) temperature field simulation of testing each test point, get final product the coefficient of heat transfer that inverse obtains contact interface between steel ingot 1 and the ingot mould 2.
The present invention includes following steps:
1) ingot mould 2 with the contact interface outside between steel ingot 1 and the ingot mould 2 is reduced to following five test segmentations, promptly has the five segment limit face coefficients of heat transfer, and every segment limit face coefficient of heat transfer difference, and all change along with the time;
1. with the horizontal surface of contact of indent on ingot tail 11 and chassis 21 (as a point among the figure to shown in the b point) as segmentation I;
2. with the horizontal surface of contact of the side direction surface of contact on ingot tail 11 and chassis 21 and ingot body 12 and mould body 22 (as b point among the figure to shown in the c point) as segmentation II;
3. with the vertical surface of contact of ingot body 12 and mould body 22 (to shown in the d point, the angle of vertical surface of contact and vertical direction is a taper of ingot as c point among the figure) as segmentation III;
4. with the vertical surface of contact of rising head 13 and riser buss 23 (as d point among the figure to shown in the e point) as segmentation IV:
5. with the horizontal surface of contact of rising head 13 and exothermic mixture 24 bottom surfaces (as e point among the figure to shown in the f point) as segmentation V.
The foundation that above-mentioned segmentation is simplified is: show that by practical experience and numerical simulation the interface heat exchange coefficient of diverse location differs less in each of division section, can simplify and think equal, and section with section between interface heat exchange coefficient differ bigger.
2) as shown in Figure 2, before pouring molten steel, in the ingot mould 2 of segmentation I～IV, determine several test points for every section, and on the position of each test point, arrange a sensor 3; After pouring molten steel, on molten steel, cover one deck exothermic mixture 24, in the exothermic mixture 24 of segmentation V, determine several test points, and on each test point, arrange a sensor 3.The lead of all the sensors 3 by separately is connected to data acquisition system (DAS) 4, the output terminal linkage interface coefficient of heat transfer inverse system 5 of data acquisition system (DAS) 4.
This step needs to determine the susceptibility of test point position and the susceptibility of test point quantity when determining test point.
Determine the susceptibility of test point position: test point is near more apart from the contact interface between steel ingot 1 and the ingot mould 2, and the speed of finding the solution of interface heat exchange coefficient inverse system 5 is fast more, promptly helps determining of interface heat exchange coefficient more.Test point is relevant with tested systems with the distance at interface, for the large-scale steel ingot more than the 100t, considering the needs that measuring accuracy and temperature field simulation and inverse are found the solution, is 20～50mm apart from the nearest test point of contact interface between steel ingot and the ingot mould and the distance at interface.
Determine the susceptibility of test point quantity:,, arrange that a characteristic test point can satisfy the requirement that this segment limit face coefficient of heat transfer inverse is found the solution corresponding to each section for the situation of interface heat exchange coefficient staging treating.And arrange two characteristic test points in the normal direction at this interface or approximate normal direction (being direction of heat flow), the experimental data of Huo Deing helps apace inverse and finds the solution the interface heat exchange coefficient that obtains this section like this.The present invention determines that only unique point quantity is 1～5.
3) become in the process of steel ingot at molten steel solidification, temperature or hot-fluid information that each sensor 3 is gathered in molten steel solidification processes, and the information data that collects is sent in the data acquisition system (DAS) 4.
4) after molten steel is frozen into steel ingot fully, read the data of each test point by data acquisition system (DAS) 4, obtain the temperature or the hot-fluid information of each test point.
5) with the data importing of each test point in steel ingot 1 and 2 interface coefficient of heat transfer inverse systems 5 of ingot mould, utilize interface heat exchange coefficient inverse system 5 to find the solution to obtain the interface heat exchange coefficient of 2 of steel ingot 1 and ingot moulds.
As shown in Figure 3, the flow process of finding the solution in 2 interface coefficient of heat transfer inverse systems 5 of steel ingot 1 and ingot mould is:
Estimate the interface heat exchange coefficient of certain time period according to the time dependent temperature data of each test point actual measurement of data acquisition system (DAS) 4 input; Carry out simulation of Temperature in the process of setting according to the interface heat exchange coefficient of estimating, to obtain the accounting temperature of test point; Relatively whether accounting temperature and observed temperature coincide, if accounting temperature and observed temperature are misfitted, then find the solution by inverse, proofread and correct interface heat exchange coefficient, and then return the simulation of solidification processing temperature field, up to the result of calculation that obtains to match with observed temperature, promptly accounting temperature and observed temperature are identical, at this time also obtained should the time period interface heat exchange coefficient, the record interface heat exchange coefficient; Whether the setting time t that judges this moment is greater than complete setting time t MaxIf, t＞t Max, then finish to calculate; Otherwise, the calculating that enters next time period.
As shown in Figure 4, the specific algorithm of the interface heat exchange coefficient inverse system 5 median surface coefficients of heat transfer is as follows.
1. begin.
2. the starting condition of certain time period is set.Corresponding first time period, then be: the starting condition of input temp field (be t=0 Temperature Distribution constantly, be made as the pouring molten steel temperature usually), and the initial value of interface heat exchange coefficient, i.e. estimated value.
Definition F (h 1..., h I) be
F ( h 1 , . . . , h I ) = &Sigma; i = 1 I &Sigma; j = 1 J &Sigma; k = 1 K ( Y k , j , i - T k , j , i ) 2 - - - ( 1 )
In the formula, h 1..., h IBe the interface heat exchange coefficient of each segmentation, Y K, j, iAnd T K, j, iRepresent temperature measured value and the calculated value of segmentation i at time period j, position k respectively, I is the segmentation sum, and J is the time period sum, and K is the test point sum of a certain segmentation i correspondence.The length of time period j can rule of thumb be chosen, and finds the solution for implicit expression, and the length of time period j can be got greatly relatively; Find the solution for explicit, the length of time period j is then got little relatively.The length of time period j can be definite value, also can dynamically change.
3. utilize interface heat exchange coefficient
Figure BSA00000147402200062
, find the solution the temperature field of this time period, i.e. steel ingot three-dimensional temperature field simulation, subscript l represents iterations, l=0,1,2 ..., wherein, l=0 is illustrated in the initial of each time period, from estimated value
Figure BSA00000147402200063
The beginning iteration.
4. at interface heat exchange coefficient
Figure BSA00000147402200064
The basis on increase an amplitude of variation ζ, generally be taken as 10 -3, utilize this new interface heat exchange coefficient
Figure BSA00000147402200065
, find the solution the temperature field of this time period.
5. according to step 3. and the temperature field of 4. finding the solution, calculate interface heat exchange coefficient
Figure BSA00000147402200066
Increment
Figure BSA00000147402200067
:
&delta; h i l + 1 = &Sigma; i = 1 I &Sigma; j = 1 J &Sigma; k = 1 K ( Y k , j , i - T k , j , i l ) &phi; k , j , i l &Sigma; i = 1 I &Sigma; j = 1 J &Sigma; k = 1 K ( &phi; k , j , i l ) 2 - - - ( 2 )
In the formula,
Figure BSA00000147402200069
The expression sensitivity coefficient, it is defined as:
&phi; k , j , i l = &PartialD; T k , j , i &PartialD; h i | h = h i l - - - ( 3 )
And the interface heat exchange coefficient that obtains next iteration step is:
h i l + 1 = h i l + &delta; h i l + 1 - - - ( 4 )
If 6. Satisfy:
&delta; h i l + 1 h i l + 1 < &epsiv; - - - ( 5 )
Then explanation
Figure BSA000001474022000614
Tally with the actual situation, it also is accurately that correspondence is found the solution the temperature field that obtains, and then, returns the calculating that 3. step enters the next time period.Wherein, ε is a decimal, generally is taken as 10 -4
If
Figure BSA00000147402200071
Do not satisfy formula (5), then to interface heat exchange coefficient Carry out assignment again, promptly
Figure BSA00000147402200073
Return step and 3. carry out iterative computation, up to satisfying formula (5).
If 7. t＞t Max, (t wherein MaxBe complete setting time), then finish to calculate; Otherwise, return the calculating that 3. step enters next time period.
In the various embodiments described above, the sensor 3 that is provided with on each test point can be thermopair, also can be the sensor of other type.
In the various embodiments described above, the solidification processing temperature field stimulation can be found the solution the transient state heat conduction equation of considering liquid-solid-phase changeable based on finite difference method, Finite Element Method or finite volume method.Because actual steel ingot is a 3D shape, therefore adopt the three-dimensional temperature field simulation.
The various embodiments described above only are used to illustrate the present invention, and wherein the structure of each parts, connected mode etc. all can change to some extent, and every equivalents of carrying out on the basis of technical solution of the present invention and improvement all should not got rid of outside protection scope of the present invention.

Claims (7)

1. the method for a definite interface heat exchange coefficient of large-sized steel ingot, it may further comprise the steps:
1) contact interface between steel ingot and the ingot mould is reduced to some test segmentations, in each segmentation, determine several test points, and a sensor is set on each test point, the lead of each sensor by separately is connected to a data acquisition system (DAS) jointly, and the output terminal of data acquisition system (DAS) connects an interface heat exchange coefficient inverse system;
2) become in the process of steel ingot at molten steel solidification, temperature or hot-fluid information in each sensor acquisition molten steel solidification process, and the information data that collects is sent in the data acquisition system (DAS);
3) after molten steel is frozen into steel ingot fully, read the data of each test point by data acquisition system (DAS), obtain the time dependent temperature or the hot-fluid information of each test point actual measurement, and import between steel ingot and ingot mould in the coefficient of heat transfer inverse system of interface;
4) interface heat exchange coefficient inverse system estimates the interface heat exchange coefficient of certain time period according to the data that step 3) imports; Carry out simulation of Temperature in the process of setting according to the interface heat exchange coefficient of estimating, to obtain the accounting temperature of test point; Relatively whether accounting temperature and observed temperature coincide, if accounting temperature and observed temperature are misfitted, then find the solution by inverse, proofread and correct interface heat exchange coefficient, and then return the simulation of solidification processing temperature field, up to the result of calculation that obtains to match with observed temperature, promptly accounting temperature and observed temperature are identical, at this moment obtain the interface heat exchange coefficient of this time period, the record interface heat exchange coefficient; Judge that whether setting time t is greater than complete setting time t MaxIf, t＞t Max, then finish to calculate; Otherwise, the calculating that enters next time period.
2. the method for a kind of definite interface heat exchange coefficient of large-sized steel ingot as claimed in claim 1, it is characterized in that: the inverse method for solving of described step 4) median surface coefficient of heat transfer inverse system is:
1. the starting condition of certain time period is set; Corresponding first time period, then be: the starting condition of input temp field, promptly t=0 Temperature Distribution constantly is made as the pouring molten steel temperature usually, and the initial value of interface heat exchange coefficient, i.e. estimated value;
Definition F (h 1..., h I) be
F ( h 1 , . . . , h I ) = &Sigma; i = 1 I &Sigma; j = 1 J &Sigma; k = 1 K ( Y k , j , i - T k , j , i ) 2 - - - ( 1 )
In the formula, h 1..., h IBe the interface heat exchange coefficient of each segmentation, Y K, j, iAnd T K, j, iRepresent temperature measured value and the calculated value of segmentation i at time period j, position k respectively, I is the segmentation sum, and J is the time period sum, and K is the test point sum of a certain segmentation i correspondence;
2. utilize interface heat exchange coefficient
Figure FSA00000147402100012
, find the solution the temperature field of this time period, i.e. steel ingot three-dimensional temperature field simulation, subscript l represents iterations, l=0,1,2 ..., wherein, l=0 is illustrated in the initial of each time period, from estimated value
Figure FSA00000147402100013
The beginning iteration;
3. at interface heat exchange coefficient
Figure FSA00000147402100014
The basis on increase an amplitude of variation ζ, generally be taken as 10 -3, utilize this new interface heat exchange coefficient
Figure FSA00000147402100021
, find the solution the temperature field of this time period;
4. according to step 2. and the temperature field of 3. finding the solution, calculate interface heat exchange coefficient
Figure FSA00000147402100022
Increment
Figure FSA00000147402100023
:
&delta; h i l + 1 = &Sigma; i = 1 I &Sigma; j = 1 J &Sigma; k = 1 K ( Y k , j , i - T k , j , i l ) &phi; k , j , i l &Sigma; i = 1 I &Sigma; j = 1 J &Sigma; k = 1 K ( &phi; k , j , i l ) 2 - - - ( 2 )
In the formula,
Figure FSA00000147402100025
The expression sensitivity coefficient, it is defined as:
&phi; k , j , i l = &PartialD; T k , j , i &PartialD; h i | h = h i l - - - ( 3 )
And the interface heat exchange coefficient that obtains next iteration step is:
h i l + 1 = h i l + &delta; h i l + 1 - - - ( 4 )
If 5.
Figure FSA00000147402100028
Satisfy:
&delta; h i l + 1 h i l + 1 < &epsiv; - - - ( 5 )
Then explanation
Figure FSA000001474021000210
Tally with the actual situation, it also is accurately that correspondence is found the solution the temperature field that obtains, and then, returns the calculating that 2. step enters the next time period; Wherein, ε is a decimal, generally is taken as 10 -4
If
Figure FSA000001474021000211
Do not satisfy formula (5), then to interface heat exchange coefficient
Figure FSA000001474021000212
Carry out assignment again, promptly
Figure FSA000001474021000213
Return step and 2. carry out iterative computation, up to satisfying formula (5);
If 6. t＞t Max, then finish to calculate; Otherwise, return the calculating that 2. step enters next time period.
3. the method for a kind of definite interface heat exchange coefficient of large-sized steel ingot as claimed in claim 1 is characterized in that: in the described step 1), the ingot mould in the contact interface outside between steel ingot and the ingot mould is reduced to five test segmentations, is respectively:
1. with the horizontal surface of contact of indent on steel ingot ingot tail and ingot mould chassis as segmentation I;
2. with the horizontal surface of contact of the side direction surface of contact on steel ingot ingot tail and ingot mould chassis and steel ingot ingot body and ingot mould mould body as segmentation II;
3. with the vertical surface of contact of steel ingot ingot body and ingot mould mould body as segmentation III;
4. with the vertical surface of contact of steel ingot rising head and ingot mould riser buss as segmentation IV;
5. with the horizontal surface of contact of the exothermic mixture bottom surface of pouring into a mould on steel ingot rising head and the steel ingot as segmentation V.
4. the method for a kind of definite interface heat exchange coefficient of large-sized steel ingot as claimed in claim 2 is characterized in that: in the described step 1), the ingot mould in the contact interface outside between steel ingot and the ingot mould is reduced to five test segmentations, is respectively:
1. with the horizontal surface of contact of indent on steel ingot ingot tail and ingot mould chassis as segmentation I;
2. with the horizontal surface of contact of the side direction surface of contact on steel ingot ingot tail and ingot mould chassis and steel ingot ingot body and ingot mould mould body as segmentation II;
3. with the vertical surface of contact of steel ingot ingot body and ingot mould mould body as segmentation III;
4. with the vertical surface of contact of steel ingot rising head and ingot mould riser buss as segmentation IV;
5. with the horizontal surface of contact of the exothermic mixture bottom surface of pouring into a mould on steel ingot rising head and the steel ingot as segmentation V.
5. as the method for claim 1 or 2 or 3 or 4 described a kind of definite interface heat exchange coefficient of large-sized steel ingot, it is characterized in that: each test test point quantity in segmentation is 1～5, is 20～50mm apart from the nearest test point of contact interface between steel ingot and the ingot mould and the distance at interface.
6. as the method for claim 1 or 2 or 3 or 4 described a kind of definite interface heat exchange coefficient of large-sized steel ingot, it is characterized in that: described sensor is a thermopair.
7. the method for a kind of definite interface heat exchange coefficient of large-sized steel ingot as claimed in claim 5, it is characterized in that: described sensor is a thermopair.
CN2010101951218A 2010-05-31 2010-05-31 Method for determining interface heat exchange coefficient of large-sized steel ingot Expired - Fee Related CN101871903B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2010101951218A CN101871903B (en) 2010-05-31 2010-05-31 Method for determining interface heat exchange coefficient of large-sized steel ingot

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2010101951218A CN101871903B (en) 2010-05-31 2010-05-31 Method for determining interface heat exchange coefficient of large-sized steel ingot

Publications (2)

Publication Number Publication Date
CN101871903A true CN101871903A (en) 2010-10-27
CN101871903B CN101871903B (en) 2012-07-04

Family

ID=42996896

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2010101951218A Expired - Fee Related CN101871903B (en) 2010-05-31 2010-05-31 Method for determining interface heat exchange coefficient of large-sized steel ingot

Country Status (1)

Country Link
CN (1) CN101871903B (en)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102507636A (en) * 2011-09-30 2012-06-20 中国科学院金属研究所 Method for measuring interfacial heat transfer coefficient of rapid cooling process of steel
CN103033277A (en) * 2012-12-07 2013-04-10 山东科技大学 Device and method for evaluating relation of interface temperature and interface heat exchange coefficient
CN103995017A (en) * 2014-04-04 2014-08-20 内蒙古科技大学 Experiment method for determining periodic heat transfer coefficient
CN104698030A (en) * 2015-03-27 2015-06-10 中南林业科技大学 Determination method for interface heat transfer coefficient in casting process
CN105445322A (en) * 2015-12-02 2016-03-30 东北大学 Method for measuring interface heat transfer coefficient under pressurization condition
CN108108529A (en) * 2017-12-01 2018-06-01 东方电气集团东方汽轮机有限公司 A kind of reverse calculation algorithms of the easy measurement cast interface coefficient of heat transfer
CN109856183A (en) * 2019-03-25 2019-06-07 上海工程技术大学 A kind of measuring method and device of the metal mold counter-pressure casting solid liquid interface coefficient of heat transfer
CN111855739A (en) * 2020-09-10 2020-10-30 东北大学 Method and system for determining heat exchange coefficient of interface between ingot and casting mold in pressurized solidification process
CN112285151A (en) * 2020-09-21 2021-01-29 北京空间飞行器总体设计部 Method for determining heat transfer coefficient of complex heat transfer component interface based on actual product
CN112464398A (en) * 2020-11-12 2021-03-09 北京适创科技有限公司 Four-dimensional interface heat exchange coefficient model establishing method for casting process

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU464626A1 (en) * 1973-12-03 1975-03-25 Московский Ордена Трудового Красного Знамени Институт Стали И Сплавов The method of controlling the heating of metal products
SU864024A1 (en) * 1979-11-22 1981-09-15 Украинский Ордена Трудового Красного Знамени Научно-Исследовательский Институт Металлов Device for measuring surface temperature, mainly of steel ingots

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU464626A1 (en) * 1973-12-03 1975-03-25 Московский Ордена Трудового Красного Знамени Институт Стали И Сплавов The method of controlling the heating of metal products
SU864024A1 (en) * 1979-11-22 1981-09-15 Украинский Ордена Трудового Красного Знамени Научно-Исследовательский Институт Металлов Device for measuring surface temperature, mainly of steel ingots

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
《大型铸锻件》 20100525 李文胜等 大型钢锭凝固过程三维数值模拟 1-4 1-7 , 第3期 2 *
《航空材料学报》 20100228 张家锋等 FGH96粉末高温合金涡轮盘淬火过程界面换热系数的研究 26-29 1-7 第30卷, 第1期 2 *

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102507636A (en) * 2011-09-30 2012-06-20 中国科学院金属研究所 Method for measuring interfacial heat transfer coefficient of rapid cooling process of steel
CN102507636B (en) * 2011-09-30 2013-11-06 中国科学院金属研究所 Method for measuring interfacial heat transfer coefficient of rapid cooling process of steel
CN103033277A (en) * 2012-12-07 2013-04-10 山东科技大学 Device and method for evaluating relation of interface temperature and interface heat exchange coefficient
CN103033277B (en) * 2012-12-07 2014-12-10 山东科技大学 Device and method for evaluating relation of interface temperature and interface heat exchange coefficient
CN103995017A (en) * 2014-04-04 2014-08-20 内蒙古科技大学 Experiment method for determining periodic heat transfer coefficient
CN103995017B (en) * 2014-04-04 2016-08-17 内蒙古科技大学 A kind of experimental technique measuring cyclical heat transmission coefficient
CN104698030A (en) * 2015-03-27 2015-06-10 中南林业科技大学 Determination method for interface heat transfer coefficient in casting process
CN105445322A (en) * 2015-12-02 2016-03-30 东北大学 Method for measuring interface heat transfer coefficient under pressurization condition
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
CN109856183A (en) * 2019-03-25 2019-06-07 上海工程技术大学 A kind of measuring method and device of the metal mold counter-pressure casting solid liquid interface coefficient of heat transfer
CN109856183B (en) * 2019-03-25 2021-12-17 上海工程技术大学 Method and device for measuring heat exchange coefficient of metal mold counter-pressure casting solid-liquid interface
CN111855739A (en) * 2020-09-10 2020-10-30 东北大学 Method and system for determining heat exchange coefficient of interface between ingot and casting mold in pressurized solidification process
CN112285151A (en) * 2020-09-21 2021-01-29 北京空间飞行器总体设计部 Method for determining heat transfer coefficient of complex heat transfer component interface based on actual product
CN112285151B (en) * 2020-09-21 2024-05-10 北京空间飞行器总体设计部 Complex heat transfer member interface heat exchange coefficient determination method based on actual product
CN112464398A (en) * 2020-11-12 2021-03-09 北京适创科技有限公司 Four-dimensional interface heat exchange coefficient model establishing method for casting process

Also Published As

Publication number Publication date
CN101871903B (en) 2012-07-04

Similar Documents

Publication Publication Date Title
CN101871903B (en) Method for determining interface heat exchange coefficient of large-sized steel ingot
CN103433448A (en) Method for determining heat flux density of continuous casting crystallizer based on flux film and air gap dynamic distribution
CN104439142B (en) A kind of for detecting Mold liquid level and the method for covering slag liquid slag layer thickness
CN102773443A (en) Method for determining heat transfer coefficient of secondary cooling zones in steel continuous casting process
CN104331629A (en) Uneven distributed calculating method of liquid, solid powder film and air gap thickness of continuous casting crystallizer casting powder
CN103542953B (en) The measuring method in sand mold temperature field in dull and stereotyped magnesium alloy sand mold casting
CN106525466B (en) A kind of Braking System for Multiple Units critical component robust filtering method and system
CN108108529B (en) Inverse calculation method for simply and conveniently measuring heat exchange coefficient of casting interface
CN108759652B (en) A kind of curvature measurement method based on favour stone full-bridge principle
CN105447302A (en) Method for estimating aerodynamic force of motor train unit
CN102205403B (en) Method for detecting local heat flow of continuous casting crystallizer copper plate
CN104001876B (en) The width of air gap of on-line measurement steel ingot and ingot mould and the method for interface heat exchange coefficient
CN113573826B (en) Device and method for estimating thickness of solidified shell in casting mold
CN102407295A (en) Method for indirectly measuring surface temperature change in casting process of large castings
JP4244675B2 (en) Detection method of width direction distribution of molten steel flow velocity in continuous casting mold
CN205263007U (en) Measure device of sand mould feeder modulus
CN102228972A (en) Calculation method for solidification heat transfer process of continuous casting crystallizer
CN112881463B (en) Method for visually processing temperature change of liquid in container
CN111723427B (en) Bridge structure damage positioning method based on recursive feature decomposition
CN111855739A (en) Method and system for determining heat exchange coefficient of interface between ingot and casting mold in pressurized solidification process
CN105259209A (en) Measuring system and computing method for LNG evaporation rate
CN106394539B (en) A kind of state monitoring method and device of high-speed train braking system critical component
JP4508556B2 (en) Monitoring system for liquefied gas carrier
JPH11338898A (en) Method for analyzing flow in molding
CN105301036A (en) Device used for measuring sand mould riser modulus

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: 20120704

Termination date: 20190531

CF01 Termination of patent right due to non-payment of annual fee