CN114921598A - Blast furnace top furnace burden movement track modeling method and system - Google Patents
Blast furnace top furnace burden movement track modeling method and system Download PDFInfo
- Publication number
- CN114921598A CN114921598A CN202210454570.2A CN202210454570A CN114921598A CN 114921598 A CN114921598 A CN 114921598A CN 202210454570 A CN202210454570 A CN 202210454570A CN 114921598 A CN114921598 A CN 114921598A
- Authority
- CN
- China
- Prior art keywords
- chute
- furnace
- burden
- movement
- charge
- 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
Images
Classifications
-
- C—CHEMISTRY; METALLURGY
- C21—METALLURGY OF IRON
- C21B—MANUFACTURE OF IRON OR STEEL
- C21B5/00—Making pig-iron in the blast furnace
-
- C—CHEMISTRY; METALLURGY
- C21—METALLURGY OF IRON
- C21B—MANUFACTURE OF IRON OR STEEL
- C21B7/00—Blast furnaces
- C21B7/18—Bell-and-hopper arrangements
- C21B7/20—Bell-and-hopper arrangements with appliances for distributing the burden
-
- C—CHEMISTRY; METALLURGY
- C21—METALLURGY OF IRON
- C21B—MANUFACTURE OF IRON OR STEEL
- C21B2300/00—Process aspects
- C21B2300/04—Modeling of the process, e.g. for control purposes; CII
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P10/00—Technologies related to metal processing
- Y02P10/25—Process efficiency
Landscapes
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Manufacturing & Machinery (AREA)
- Materials Engineering (AREA)
- Metallurgy (AREA)
- Organic Chemistry (AREA)
- Blast Furnaces (AREA)
Abstract
The invention discloses a modeling method and a system for movement tracks of furnace burden on the top of a blast furnace, which establish a static coordinate system which is static relative to the blast furnace and a moving coordinate system which rotates together with a chute, establish a mathematical model for discharge of a burden throttling valve according to the discharge mode of the burden at a blanking tank, establish a mathematical model for free fall of the burden at a central throat, establish an initial movement mathematical model of the burden on the chute after the burden collides with the chute according to the collision recovery coefficient of the burden and the chute, establish a mathematical model for chute movement of the burden in the chute and obtain a mathematical model for movement tracks of the burden from the throttling valve to the tail end of the chute of single-particle burden according to the mathematical model for chute movement, thereby obtaining the mathematical model for movement tracks of the burden on the top of the whole blast furnace, solving the technical problem of low calculation precision of the movement tracks of the burden on the top of the blast furnace in the existing blast furnace, realizing the calculation of the movement tracks of the burden on the top of the blast furnace in different initial movement states, the calculation precision of the furnace charge motion trail is improved.
Description
Technical Field
The invention mainly relates to the technical field of blast furnace smelting, in particular to a method and a system for modeling a movement track of furnace burden at the top of a blast furnace.
Background
Blast furnace iron making is a key process of a steel production flow and is a production process of continuous blast, batch distribution, periodic iron tapping, accompanied complex physical and chemical reactions, severe substance and phase change conversion, high-strength energy transfer and transfer. The three-dimensional shape of the blast furnace top charge level is the main basis for optimizing and adjusting the gas flow distribution of the blast furnace, and is also the important basis for assisting a blast furnace smelting expert to find abnormal furnace conditions in time, appropriately adjust the charge distribution system and avoid further deterioration of the furnace conditions. And the movement track of furnace top furnace burden of the blast furnace directly determines the three-dimensional shape of the burden surface, and the successful prediction of the three-dimensional movement track of the furnace top furnace burden has great significance for promoting the stable and smooth operation of the blast furnace, improving the utilization rate of coal gas and ensuring the quality of molten iron. Because of the severe environments of closed blast furnace top, high temperature and high pressure, high dust, weak light and the like, the conventional mechanical trial rod, radar trial rod, laser trial rod and industrial endoscope are difficult to stably detect the charge level shape for a long time in the environments. Therefore, the movement situation of the furnace burden from the blanking tank to the burden surface needs to be analyzed based on the movement mechanism of the furnace burden, a three-dimensional mathematical model of the movement track of the furnace burden is established, and the movement track and the drop point position of the furnace burden under different blast furnace burden distribution operation parameters are calculated.
The furnace burden is conveyed from the storage bin to the interior of the blast furnace and then passes through the charging bucket, the discharging bucket, the central throat pipe and the rotary chute on the top of the blast furnace in sequence, and finally falls to the charge level to form a new charge level shape. On the top of the serial tank type bell-less blast furnace, the movement process of the charging material from the charging bucket to the discharging bucket almost has no influence on the charging material distribution operation. Therefore, the model of the motion track of the charge material on the top of the blast furnace can be divided into five parts. Firstly, a throttle valve discharging model away from the tail end of a discharging tank; then a free fall motion model at the central throat; then, a model of collision with the chute is carried out; next, a three-dimensional spiral motion model on the rotating chute is obtained; and finally, an oblique throwing motion model of the furnace throat dead zone. At present, a lot of furnace burden motion trajectory models are available, but the furnace burden is only considered as mass points to be subjected to stress analysis, and the obtained furnace burden motion trajectory is a motion trajectory of the mass points and cannot represent the position distribution and speed distribution conditions of the furnace burden when the furnace burden leaves a chute. Therefore, the modeling method for the furnace top furnace burden motion trail of the blast furnace is provided, the accurate calculation of the furnace top furnace burden motion trail of the blast furnace is realized, and the modeling method has important significance for the accurate material distribution operation of the blast furnace.
The invention discloses a blast furnace bell-less multi-ring distribution mathematical model, which comprehensively considers the influence of factors such as a stockline position, the change of the movement distance of furnace materials on a chute along with the change of the inclination angle of the chute and the like, quantitatively analyzes the initial distribution of the furnace materials in the furnace, provides a model of the collision speed of the furnace materials with the chute through a central throat pipe, a model of the movement speed of the furnace materials leaving the chute, a model of the movement stress of the furnace materials in a throat dead zone and the like, and realizes the bell-less multi-ring distribution of the blast furnace.
However, the material flow is regarded as mass points for modeling, the obtained motion trail of the material flow is a single-point motion trail, and the position and the speed distribution condition of the whole material flow at the tail end of the chute are difficult to obtain; meanwhile, the falling point of the furnace burden on the charge level is also a single point, so that the distribution condition of the furnace burden on the charge level is difficult to obtain.
The invention discloses a modeling method and a modeling system for a furnace charge movement locus on a U-shaped chute of a blast furnace, and the modeling method and the modeling system are mainly used for the furnace charge movement locus on the U-shaped chute, a static coordinate system and a moving coordinate system are established, the stress condition of the furnace charge is analyzed in the U-shaped chute, and a mathematical model of the movement locus of the furnace charge relative to the U-shaped rotating chute is obtained according to Newton's second law. After the initial movement state of the particles on the U-shaped chute is set, the movement track of the furnace burden on the U-shaped chute can be obtained.
However, in the invention, the initial motion state of the furnace charge particles in the chute is set artificially, the influence of the initial position of the furnace charge at the throttle valve on the motion trajectory of the furnace charge is not considered, and the calculation of the motion trajectory of the furnace charge at the top of the blast furnace is difficult to realize.
Disclosure of Invention
The modeling method and the system for the movement track of the furnace top furnace burden of the blast furnace solve the technical problem that the calculation precision of the movement track of the furnace top furnace burden of the blast furnace is low.
In order to solve the technical problem, the modeling method for the movement track of the furnace top burden of the blast furnace provided by the invention comprises the following steps:
establishing a static coordinate system which is static relative to the blast furnace and a movable coordinate system which rotates together with the chute;
establishing a furnace burden throttling valve discharge mathematical model according to the discharge mode of the furnace burden at the blanking tank;
establishing a free falling mathematical model of the furnace charge in the central throat, wherein the free falling mathematical model is used for acquiring the speed and the position of the furnace charge in a static coordinate system before the furnace charge collides with the chute;
establishing an initial motion mathematical model of the furnace burden on the chute after the furnace burden collides with the chute according to a collision recovery coefficient of the furnace burden colliding with the chute, wherein the collision recovery coefficient comprises a normal recovery coefficient and a radial recovery coefficient;
according to an initial movement mathematical model of the furnace burden on the chute after the furnace burden collides with the chute, establishing a chute movement mathematical model of the furnace burden in the chute;
and acquiring a furnace charge movement locus mathematical model of the single-particle furnace charge from the throttle valve to the tail end of the chute according to the chute movement mathematical model, and further acquiring a movement locus mathematical model of the whole furnace top furnace charge.
Further, establishing a static coordinate system which is static relative to the blast furnace and a moving coordinate system which rotates together with the chute comprises:
establishing a static coordinate system which is static relative to the blast furnace according to the right-hand rule, wherein the origin of the static coordinate system is the intersection point between the horizontal connecting line of the two connecting points of the chute and the central throat pipe and the symmetrical axis of the blast furnace;
according to the right-hand rule, sequentially rotating the static coordinate system around two coordinate axes in the static coordinate system to obtain a dynamic coordinate system rotating together with the chute;
and solving a coordinate transformation matrix between the static coordinate system and the moving coordinate system.
Further, establishing a mathematical model of the free fall of the charge in the central throat comprises:
obtaining the coordinate of the furnace charge in a static coordinate system when the furnace charge leaves the throttling valve according to the furnace charge throttling valve discharge mathematical model;
according to the coordinates of the furnace charge in a static coordinate system when leaving the throttling valve, a free fall mathematical model of the furnace charge in the central throat is established, and the specific calculation formula of the free fall mathematical model of the furnace charge in the central throat is as follows:
wherein r is 1 And v 1 Respectively the speed and the position of the furnace charge in a static coordinate system before the collision with the chute, x 0 、y 0 、h a Respectively the coordinates of the X-axis, the Y-axis and the Z-axis in the static coordinate system when the burden leaves the throttling valve, h w Is the distance between the tail end of the central throat pipe and the collision point of the furnace charge and the chute, v 0 The Z-axis speed of the furnace charge leaving the throttle valve, g is the gravity acceleration, e is the chute tilting distance, beta 0 And gamma 0 Respectively the horizontal rotation angle and the inclination angle theta of the chute relative to the initial position when the furnace burden collides with the chute 0 Is an included angle between the furnace burden and the symmetrical axis of the chute, and
and R is the chute radius.
Further, according to a collision recovery coefficient of the furnace burden colliding with the chute, establishing an initial motion mathematical model of the furnace burden on the chute after the furnace burden collides with the chute comprises the following steps:
the speed and the position of the furnace burden in a static coordinate system at the moment before the furnace burden collides with the chute are converted into the moving speed and the moving position in a moving coordinate system;
obtaining normal phase vectors of collision points on the chute according to the movement position of the furnace charge in the moving coordinate system before the furnace charge collides with the chute;
obtaining the incident speed of the furnace burden before the furnace burden collides with the chute according to the movement speed of the furnace burden before the furnace burden collides with the chute in the moving coordinate system;
calculating an included angle between a normal phase vector of a collision point on the chute and an incident speed of the furnace burden before the furnace burden collides with the chute to obtain a first included angle;
obtaining an initial speed constant of the furnace burden moving on the chute after the furnace burden collides with the chute according to a collision recovery coefficient and a first included angle when the furnace burden collides with the chute;
and calculating the emergent speed of the furnace burden after the furnace burden collides with the chute according to the initial speed constant of the furnace burden moving on the chute after the furnace burden collides with the chute.
Further, a specific formula for calculating the exit speed of the furnace burden after the furnace burden collides with the chute is as follows:
wherein, v' 1 Is incident velocity v 'of furnace burden before collision with chute' 2 The emergence speed of the furnace charge after colliding with the chute, | | · | | represents the modular length operation, e n Is the normal coefficient of restitution, theta, of the charge when it collides with the chute int Is a first angle θ out Is an included angle between a normal phase vector of a collision point on the chute and the emergent speed of the furnace charge after the collision with the chute, and
and n is a normal phase vector of a collision point on the chute.
Further, according to the initial velocity constant of the movement on the chute after the furnace burden collides with the chute, the calculation formula for calculating the emergence velocity of the furnace burden after the collision with the chute is as follows:
v 2 ′＝av 1 ′+bn，
wherein, v' 1 Is incident velocity v 'of furnace burden before collision with chute' 2 The speed is the emergent speed of the furnace burden after the furnace burden collides with the chute, a and b are initial speed constants of the furnace burden moving on the chute after the furnace burden collides with the chute, and n is a normal phase vector of a collision point on the chute.
Further, according to the initial movement mathematical model of the furnace burden on the chute after the furnace burden collides with the chute, the establishment of the chute movement mathematical model of the furnace burden in the chute comprises the following steps:
obtaining an initial movement position and an initial movement speed of the furnace burden on the chute according to an initial movement mathematical model of the furnace burden on the chute after the furnace burden collides with the chute;
respectively analyzing the relative motion and the traction motion of the furnace burden relative to the chute according to the initial motion position and the initial motion speed of the furnace burden on the chute and the point-based composite motion, further analyzing the absolute motion of the furnace burden in the blast furnace, and establishing a chute motion mathematical model of the furnace burden in the chute according to Newton's second law;
and obtaining the movement speed and the movement position of the burden from the throttling valve to the tail end of the chute according to the chute movement mathematical model.
Further, according to the chute motion mathematical model, after obtaining the tail end motion speed and the tail end motion position of the burden from the throttling valve to the chute tail end, the method also comprises the following steps:
establishing a furnace throat empty area movement mathematical model, wherein the calculation formula of the furnace throat empty area movement mathematical model is as follows:
wherein S is x 、S y And S z Represents the moving distance of the furnace charge on the X axis, the Y axis and the Z axis of the throat area, t is the moving time of the furnace charge from the tail end of the chute to the charge level, v cx ,v cy ,v cz Respectively representing the moving speed of the burden material leaving the tail end of the chute;
obtaining the position of a falling point of the furnace burden according to the furnace throat empty area movement mathematical model, wherein the calculation formula of the position of the falling point of the furnace burden is as follows:
wherein r is xc 、r yc And r zc Respectively, representing the position of movement of the charge material away from the end of the chute.
The invention provides a modeling system for furnace top furnace burden movement locus of a blast furnace, comprising:
the device comprises a memory, a processor and a computer program which is stored on the memory and can run on the processor, and the steps of the modeling method for the movement locus of the furnace top burden of the blast furnace provided by the invention are realized when the processor executes the computer program.
Compared with the prior art, the invention has the advantages that:
the invention provides a modeling method and a system for movement tracks of furnace burden on the top of a blast furnace, which establish a static coordinate system which is static relative to the blast furnace and a moving coordinate system which rotates together with a chute, establish a discharge mathematical model of a burden throttling valve according to the discharge mode of the burden at a blanking pot, establish a free falling mathematical model of the burden at a central throat, establish an initial movement mathematical model of the burden on the chute after the burden collides with the chute according to the collision recovery coefficient of the burden and the chute, establish a chute movement mathematical model of the burden in the chute according to the initial movement mathematical model of the burden on the chute after the burden collides with the chute, and obtain a movement track mathematical model of the burden from the throttling valve to the tail end of the chute of single-particle burden according to the chute movement mathematical model, thereby obtaining the movement track mathematical model of the burden on the top of the whole blast furnace, and solving the technical problem of low calculation precision of the movement tracks of the burden on the top of the blast furnace, the movement track calculation of the furnace material particles at the top of the blast furnace in different initial movement states is realized by calculating the initial movement state of the single-particle furnace material and the movement state of the rotating chute, so that the calculation precision of the movement track of the furnace material is greatly improved, the function of a traditional material flow movement track model can be realized, the position and the speed distribution condition of the whole material flow at the tail end of the chute can be solved according to the movement tracks of the particles in different initial movement states, and the falling point position of the whole material flow on the material surface of the blast furnace can be solved.
Specifically, the invention aims to provide a method and a system for modeling a movement track of furnace burden at the top of a blast furnace, which are used for accurately calculating the movement track of the furnace burden at the top of the blast furnace. The movement of the furnace burden on the top of the blast furnace is in the form of material flow, the current furnace burden movement trajectory mainly treats the furnace burden as mass points to be subjected to stress analysis, and the obtained furnace burden movement trajectory is only related to the initial speed of the furnace burden, the chute inclination angle, the chute rotation speed, the collision recovery coefficient and the model coefficient of the furnace burden and the chute in contact. Meanwhile, the traditional material flow motion trajectory model considers that the initial motion position of the burden on the chute is only related to the chute tilting distance and the chute tilting angle, and the motion speed is only along the chute longitudinal direction. In fact, the initial movement state of the charge on the chute after its impact with the chute is not only related to the initial movement state of the charge particles when they leave the throttle valve, but also to the movement state of the chute during its free fall in the central throat. Therefore, the accuracy of the burden motion trail predicted by the traditional burden motion trail model is low. Aiming at the defects of a traditional furnace charge movement track model, the invention provides a mathematical modeling method for a furnace top furnace charge movement track of a blast furnace based on coordinate transformation. Meanwhile, the furnace burden motion model is subjected to vectorization analysis, so that the complexity of stress analysis in the furnace burden motion process is reduced, and the calculation accuracy of the furnace burden motion track is improved.
The invention aims to comprehensively consider the initial motion state of single-particle furnace burden and the motion state of a rotary chute and realize the motion track calculation of the furnace burden particles in different initial motion states on the top of a blast furnace. The invention can realize the function of the traditional material flow motion track model, and can calculate the position and the speed distribution condition of the whole material flow at the tail end of the chute according to the motion tracks of particles in different initial motion states, thereby calculating the position of the falling point of the whole material flow on the material surface of the blast furnace.
The invention provides a method for modeling movement tracks of furnace top burden of a serial tank type bell-less blast furnace based on coordinate transformation, wherein the model built by the method comprises a throttle valve discharge model, a central throat pipe descending model, a burden and chute collision model, a rotary chute movement model and a furnace throat empty area inclined throwing movement model; the method comprises the steps of establishing a static coordinate system and a moving coordinate system, calculating the descending height of the furnace burden in a central throat pipe based on a coordinate transformation matrix, and further calculating the three-dimensional movement speed of the furnace burden after the furnace burden collides with a chute; the method considers the influence of the initial motion state of the charging particles and the motion state of the chute on the movement track of the charging, can calculate the position distribution and the speed distribution of the material flow at the tail end of the chute, meets the distribution requirement of the serial tank type bell-less blast furnace, and has great application value.
The key points of the invention comprise:
(1) the mathematical modeling method for the furnace top furnace charge movement locus of the serial-tank type bell-less blast furnace based on coordinate transformation is provided, so that the movement locus of the furnace charge on the furnace top of the blast furnace is accurately calculated;
(2) a static coordinate system which is static relative to the blast furnace and a moving coordinate system which moves relative to the blast furnace are constructed, and the three-dimensional movement position and speed of a collision point in the chute when the furnace burden collides with the chute are calculated based on a coordinate transformation matrix;
(3) vectorization modeling is carried out on the movement process of the furnace burden, and the relation among the relative movement, the involved movement and the absolute movement of the furnace burden in the chute is analyzed on the basis of the compound movement of points, so that the complexity of stress analysis of the furnace burden is reduced, and the accuracy of a mathematical model of the movement locus of the furnace burden is improved;
(4) the calculation of the charge motion tracks of the charge particles in different initial motion states at the top of the blast furnace is realized, and the model can obtain the position distribution and the speed distribution of the charge at the tail end of the chute.
Drawings
FIG. 1 is a schematic view of the top of a tandem-type bell-less blast furnace according to a second embodiment of the present invention;
FIG. 2 is a diagram illustrating a relationship between a static coordinate system and a moving coordinate system according to a second embodiment of the present invention;
FIG. 3 is a flowchart of a modeling method of a movement trajectory of furnace top burden of a blast furnace according to a second embodiment of the present invention;
FIG. 4 is a schematic view showing the position of a burden relative to a chute according to a second embodiment of the present invention;
FIG. 5 is a schematic view of the O ' X ' Z ' plane of the charge particles in the chute according to the second embodiment of the present invention;
FIG. 6 is a schematic view of the O ' Y ' Z ' plane of the charge particles in the chute according to the second embodiment of the present invention;
FIG. 7 is a flow chart of a blast furnace top burden movement locus modeling method according to a third embodiment of the invention;
fig. 8 is a structural block diagram of a blast furnace top charge movement trajectory modeling system according to an embodiment of the present invention.
Reference numerals:
1. weighing a material tank; 2. a throttle valve; 3. a central throat; 4. rotating the chute; 5. a throat; 6. material level; 7. charge particles; 8. a stationary coordinate system; 9. a stream; 10. a memory; 20. a processor.
Detailed Description
In order to facilitate an understanding of the invention, the invention will be described more fully and in detail below with reference to the accompanying drawings and preferred embodiments, but the scope of the invention is not limited to the specific embodiments below.
The embodiments of the invention will be described in detail below with reference to the drawings, but the invention can be implemented in many different ways as defined and covered by the claims.
Example one
The modeling method for the movement track of the furnace top burden of the blast furnace provided by the embodiment of the invention comprises the following steps:
step S101, establishing a static coordinate system which is static relative to the blast furnace and a moving coordinate system which rotates together with the chute;
s102, establishing a furnace burden throttling valve discharge mathematical model according to a discharge mode of furnace burden at a blanking tank;
step S103, establishing a free fall mathematical model of the furnace burden in the central throat, wherein the free fall mathematical model is used for acquiring the speed and the position of the furnace burden in a static coordinate system before the furnace burden collides with a chute;
step S104, establishing an initial motion mathematical model of the furnace burden on the chute after the furnace burden collides with the chute according to the collision recovery coefficient of the furnace burden colliding with the chute;
step S105, establishing a chute motion mathematical model of the furnace burden in the chute according to an initial motion mathematical model of the furnace burden on the chute after the furnace burden collides with the chute;
and S106, acquiring a charge movement track mathematical model of the single-particle charge from the throttling valve to the tail end of the chute according to the chute movement mathematical model, and further acquiring a movement track mathematical model of the whole blast furnace top charge.
The modeling method of the furnace top burden movement track of the blast furnace provided by the embodiment of the invention comprises the steps of establishing a static coordinate system which is static relative to the blast furnace and a moving coordinate system which rotates together with a chute, establishing a burden throttling valve discharge mathematical model according to the discharge mode of burden at a blanking pot, establishing a free falling mathematical model of the burden at a central throat, establishing an initial movement mathematical model of the burden on the chute after the burden collides with the chute according to the collision recovery coefficient of the burden and the chute, establishing a chute movement mathematical model of the burden in the chute according to the initial movement mathematical model of the burden on the chute after the burden collides with the chute, and obtaining the burden movement track mathematical model of single-particle burden from the throttling valve to the tail end of the chute according to the chute movement mathematical model, thereby obtaining the movement track mathematical model of the whole blast furnace top burden and solving the technical problem of low calculation precision of the burden movement track of the blast furnace top burden of the existing blast furnace, the movement track calculation of the furnace material particles at the top of the blast furnace in different initial movement states is realized by calculating the initial movement state of the single-particle furnace material and the movement state of the rotating chute, so that the calculation precision of the movement track of the furnace material is greatly improved, the function of a traditional material flow movement track model can be realized, the position and the speed distribution condition of the whole material flow at the tail end of the chute can be solved according to the movement tracks of the particles in different initial movement states, and the falling point position of the whole material flow on the material surface of the blast furnace can be solved.
Specifically, the movement of the furnace burden on the top of the blast furnace is in the form of material flow, the current furnace burden movement trajectory mainly treats the furnace burden as a mass point for stress analysis, and the obtained furnace burden movement trajectory is only related to the initial speed of the furnace burden, the chute inclination angle, the chute rotation speed, the collision recovery coefficient and the model coefficient of the furnace burden and the chute in contact. Meanwhile, the traditional material flow motion trajectory model considers that the initial motion position of the burden on the chute is only related to the chute tilting distance and the chute tilting angle, and the motion speed is only along the chute longitudinal direction. In fact, the initial movement of the charge on the chute after its impact with the chute is not only related to the initial movement of the charge particles as they leave the throttle valve, but also to the movement of the chute during its free fall in the central throat. Therefore, the accuracy of the burden motion trail predicted by the traditional burden motion trail model is low. Aiming at the defects of a traditional furnace charge movement locus model, the embodiment of the invention provides a coordinate transformation-based mathematical modeling method for the furnace top furnace charge movement locus of a blast furnace. Meanwhile, the furnace burden motion model is subjected to vectorization analysis, so that the complexity of stress analysis in the furnace burden motion process is reduced, and the calculation accuracy of the furnace burden motion track is improved.
Example two
The embodiment of the invention discloses a mathematical model of a furnace top charge movement locus of a serial tank type bell-less blast furnace based on coordinate transformation, and figure 1 is a schematic diagram of the furnace top of the serial tank type bell-less blast furnace of the embodiment, and comprises a weighing charging bucket 1, a throttle valve 2, a central throat pipe 3, a rotary chute 4, a furnace throat 5, a charge level 6, charge particles 7, a static coordinate system 8 and a material flow 9. Fig. 2 is a relation diagram of a static coordinate system and a moving coordinate system. FIG. 3 is a mathematical model solving implementation step of the movement locus of the single-particle furnace burden at the top of the blast furnace, which comprises the following steps:
(1) establishing a static coordinate system OXYZ which is static relative to the blast furnace and a movable coordinate system O 'X' Y 'Z' which rotates together with the chute according to the right-hand rule, and solving a coordinate transformation matrix between the static coordinate system and the movable coordinate system;
(2) establishing a furnace charge throttling valve discharge mathematical model according to the discharge mode of the furnace charge at the blanking tank;
(3) calculating the total height of the furnace burden which leaves the throttle valve in a certain initial motion state and falls off the central throat pipe when the furnace burden collides with the chute according to the coordinate transformation matrix, and establishing a free fall mathematical model of the furnace burden on the central throat pipe;
(4) converting the motion information of the furnace burden in the static coordinate system into a dynamic coordinate system, and establishing an initial motion mathematical model of the furnace burden on the chute after the furnace burden collides with the chute according to the speed loss coefficient of the furnace burden colliding with the chute;
(5) on the basis of the known initial movement position and movement speed of the furnace burden on the chute, respectively analyzing the relative movement of the furnace burden relative to the chute and the traction movement of the chute based on the point composite movement, further analyzing the absolute movement of the furnace burden in the blast furnace, and establishing a mathematical model of the furnace burden in the chute according to Newton's second law;
(6) according to the mathematical model of the movement track of the single-particle furnace burden from the throttling valve to the chute tail end, the movement track of the furnace burden on the blast furnace top when the furnace burden particles with different initial movement positions and movement speeds leave the throttling valve is analyzed, and the mathematical model of the movement track of the whole furnace burden is formed.
(7) And analyzing the stress condition of the furnace burden at the furnace throat of the blast furnace by taking the movement position and the movement speed of the furnace burden at the tail end of the chute as initial conditions, and establishing a mathematical model of the movement of the furnace burden in the empty area of the furnace throat.
The specific implementation scheme is as follows:
(1) establishing a static coordinate system, a dynamic coordinate system, a static coordinate system, a coordinate transformation matrix and a computer readable medium
In order to describe the movement track of furnace charge on the top of the bell-less blast furnace, a static coordinate system OXXYZ which is static relative to the blast furnace and a moving coordinate system O 'X' Y 'Z' which moves relative to the blast furnace are established according to the right-hand coordinate system rule.
(a) Establishing a static coordinate system OXYZ:
taking the symmetrical center line of the blast furnace as the Z axis of a static coordinate system, and taking the vertical downward direction as the positive direction of the Z axis of the static coordinate system; the intersection point between the horizontal connecting line of the two connecting points of the U-shaped rotating chute and the central throat pipe and the symmetrical axis of the blast furnace is used as the origin O of the static coordinate system; a horizontal line which is vertical to the OZ and intersects with the point O is an X axis, and the horizontal right direction is taken as the positive direction of the X axis of the static coordinate system; the straight line perpendicular to the plane OXZ and intersecting point O is the Y-axis, and the straight line perpendicular to the paper surface is taken as the positive direction of the Y-axis of the static coordinate system, specifically as shown by the coordinate system OXYZ in FIG. 1.
(b) Establishing a moving coordinate system OX ' Y ' Z ':
according to the actual operation of the blast furnace chute, any position of the blast furnace chute in the operation process can be reached through two times of rotation in the original static state. According to the right-hand rule, the rotation angle beta is set to rad with the positive direction of OZ axis as the positive direction of rotation around the axis, the OX axis as the initial position, and the rotation coordinate system OX is reached once r Y r Z r (ii) a Then on the basis of the first rotation, with OY r Positive axial direction as positive direction of the second rotation around the axis, in OZ r The axis is the initial position of the second rotation, the rotation angle γ, unit is rad, and reaches the moving coordinate system O 'X' Y 'Z', and the specific coordinate transformation diagram is shown in fig. 2.
(c) Calculating a coordinate transformation matrix:
when transforming the coordinate system, the coordinate system is rotated around a certain axis, and only the plane perpendicular to the rotation axis is actually located. Therefore, the relation of the furnace charge in the static coordinate system and the dynamic coordinate system is
Wherein (x, y, z) is the motion state information of the furnace charge in the static coordinate system, and (x ', y ', z ') is the motion state information of the furnace charge in the dynamic coordinate system.
Which respectively represent coordinate transformation matrices when the coordinate system is rotated about the Y-axis and the Z-axis.
(2) Throttle valve discharge mathematical model
The discharge mode of the charging material at the charging bucket is a funnel type, and can be calculated by a hydraulics continuity equation, which is described as follows in a static coordinate system:
wherein Q is the mass flow rate of the burden leaving the throttling valve, the unit is kg/S, rho is the bulk density of the burden, and S is the projection plane of the throttling valveProduct, L s The length of the periphery of the throttle valve is shown, and d is the diameter of the furnace charge. The position of the charge in the static coordinate system when leaving the throttle is expressed as:
r 0 ＝(x 0 ,y 0 ,h a ) T (3)
wherein x is 0 、y 0 、h a Respectively, the X-, Y-, and Z-axis coordinates in the static coordinate system as the particle exits the throttle valve.
(3) Central throat free-falling mathematical model
When the throttle valve is opened, the charge particles are driven by gravity to move in a V mode 0 Moving vertically downwards, falling on the rotating chute after passing through the central throat, the position where the burden collides with the chute can be expressed as
r 1 ＝(x 0 ,y 0 ,h w ) T (4)
Wherein x is 0 、y 0 X-axis and Y-axis coordinates, h, respectively, of the particles as they leave the throttle valve w Is the distance from the tail end of the central throat pipe to the collision point of the furnace burden and the chute. If the burden collides with the chute, the chute is horizontally rotated by beta relative to the initial position 0 Is inclined by gamma 0 Then h is w Can be described as:
wherein e is the chute tilting distance, R is the chute radius,
is the angle between the charge particles and the chute axis of symmetry when theta 0 Positive values are given for the negative half axis of the Y' axis. The charge is only subjected to the action of gravity during the movement of the central throat, so that the speed of movement of the charge immediately before impact with the chute is expressed as:
wherein v is 0 Is the Z-axis velocity of the charge as it leaves the choke, and g is the acceleration of gravity.
(4) A furnace charge and chute collision mathematical model;
when the furnace burden collides with the chute, the chute rotates beta around the Z axis and gamma around the Y axis relative to the static coordinate system. Thus, the velocity v of the charge at the moment before it collides with the chute 1 And position r 1 Expressed as follows in the moving coordinate system:
wherein v is 1 ' represents the moving speed of the furnace charge in a moving coordinate system when the furnace charge is about to collide with the chute, r 1 ' denotes the position of the point of impact of the burden with the chute in the moving coordinate system, i.e. the initial position of the movement of the burden on the chute. Because the horizontal direction velocity of the furnace charge is ignored when the furnace charge moves in the central throat, the incident velocity direction before the furnace charge collides with the chute can be V 1 The direction of' indicates. The curved surface of the chute can be expressed as f (x ', y ', z ') in the moving coordinate system, and then the collision point r on the chute 1 The normal phase vector of' is expressed as:
an included angle theta between the incident speed of the furnace burden before the furnace burden collides with the chute and the normal vector of the chute int Expressed as:
the emergent speed of the furnace charge after the collision with the chute is set as v 2 ′＝(v x,2 ′,v y,2 ′,v z,2 ') the magnitude of the exit velocity of the charge is expressed as
Wherein v' 2,x 、v′ 2,y And v' 2,z The magnitudes of the velocity components of the exit velocity in the X, Y, and Z axes are shown, respectively. The included angle between the emergent speed of the furnace charge after collision with the chute and the normal vector of the collision point is expressed as
The normal speed and the tangential speed of the furnace burden after the furnace burden collides with the chute have loss, and the speed after the collision is related to the collision recovery coefficient. The collision recovery coefficient of the charge with the chute is expressed as:
wherein e n Is a normal coefficient of restitution when the furnace charge collides with the chute, is related to the material of the collision body and is a fixed value, e t For the tangential restitution coefficient when furnace charge and chute collide, relevant with the material and the collision angle of collision body, can express as:
e t ＝1-μ(1+e n )cotθ int (14)
wherein mu is the sliding friction coefficient between the furnace burden and the chute. The angle between the exit velocity and the normal vector of the collision point can be found according to the equations (13) and (14), and is expressed as:
therefore, the exit velocity of the furnace burden after the furnace burden collides with the chute can be obtained according to the formulas (12) and (15) and is expressed as follows:
simultaneously, the exit velocity after furnace charge and chute collision and the incident velocity before the collision and the normal vector of collision point are in the coplanar, satisfy promptly:
v 2 ′＝av 1 ′+bn (17)
converting equation (17) into a scalar expression:
according to the known parameters and the formulas (10), (11) and (18), the initial speed constants a, b and v of the furnace burden moving on the chute after the furnace burden collides with the chute can be obtained 2 ′＝(v x，2 ′,v y，2 ′,v z，2 ′)。
(5) Rotary chute motion mathematical model
When the chute rotates, the burden is subjected to gravity, supporting force, friction force, Coriolis force and the like on the chute, and besides the direction and the magnitude of the gravity are unchanged, the other four forces change along with the movement position of the burden on the chute. For this purpose, a point-based compound motion method is proposed for analyzing a three-dimensional motion mathematical model of the furnace charge in the rotating chute.
(a) Relative movement:
the position of charge P relative to the chute is shown in fig. 4. The position of charge P relative to the chute may be represented by P (X ', θ), where X' represents the projection of charge P on the axis O 'X' in the moving coordinate system O 'X' Y 'Z', as shown in FIG. 5, and θ is the angle between charge P and the symmetry axis of the chute, as shown in FIG. 6. The radial dimension of the charge material P relative to the moving coordinate system O 'X' Y 'Z' is expressed as:
wherein R isThe radius of the chute is such that,
the segregation angle of the furnace burden on the chute is defined as theta which is segregated in a positive angle in the negative direction of the Y 'axis, segregated in a negative angle in the positive direction of the Y' axis, and e is the tilting distance of the chute. In the moving coordinate system O 'X' Y 'Z', the relative velocity of the charge P is expressed as:
the relative acceleration of charge P is expressed as:
(b) and (3) carrying out traction movement:
in the moving coordinate system, the pipeline points of the moving coordinate system superposed with the moving points are deleted to form the connecting points, and the vector diameter of the connecting points relative to the moving coordinate system is the same as the vector diameter of the furnace charge P relative to the moving coordinate system and is r'. The sagittal of the tie point relative to the static coordinate is expressed as:
r e ＝r o′ +r′ (22)
wherein r is o′ Represents the position of the origin O' of the moving coordinate system within the stationary coordinate system, here zero. The formula (22) is subjected to time derivation, and the involvement velocity v of the moving point can be obtained e Expressed as:
v e ＝ω×r′ (23)
where ω is the chute rotational angular velocity. The formula (23) is derived from the time to obtain the involved acceleration a of the moving point e Expressed as:
a e ＝a×r′+ω×ω×r′ (24)
wherein a is the chute rotation angular acceleration.
(c) Absolute motion:
when the chute rotates, the absolute motion of the moving point is equal to the vector sum of the relative acceleration, the involved acceleration and the Coriolis acceleration, and is expressed as:
a a ＝a r +a e +a c (25)
wherein a is c Is the Coriolis acceleration of charge P, expressed as:
a c ＝2ω×v r (26)
by substituting (25) into equations (21), (24), and (26), the absolute acceleration can be obtained:
wherein
(d) And (3) stress analysis:
in the moving coordinate system, the furnace burden is subjected to gravity, supporting force and friction force. The weight force to which the charge is subjected can be expressed as:
wherein m is the mass of the furnace charge P in kg, g is the acceleration of gravity, m/s 2 . The supporting force to which the charge is subjected is expressed as:
wherein F N The supporting force applied to the charge P. The sliding friction experienced by the charge is expressed as:
wherein mu is a dynamic friction factor between the furnace burden P and the U-shaped chute. Therefore, the furnace burden is in a moving coordinate systemInternal applied resultant force F Σ Expressed as:
F Σ ＝G′+F N ′+F f ′ (31)
according to Newton's second law, the relation between the absolute acceleration and the resultant external force of the charging material P can be obtained, and the relation is expressed as follows:
F Σ ＝ma a (32)
substituting equations (27) and (31) into equation (32) and reducing:
wherein
The movement position, the speed and the acceleration of the furnace burden on the chute at different moments can be solved by adopting a four-order Runge-Kutta algorithm. Furthermore, the position r of the burden on the chute at the moment i can be obtained according to the formulas (19) and (20) i ' sum velocity v ri . Then the position and the speed of the furnace charge in the static coordinate system can be obtained according to the formula (7):
(6) and (4) a motion mathematical model of the furnace throat dead zone.
When no gas flow influence exists, the furnace burden is only under the action of vertical downward gravity in the furnace throat dead zone, so that the furnace burden does uniform motion on an X axis and a Y axis and does uniform accelerated motion on a Z axis, and the three-dimensional motion distance of the furnace burden in the dead zone is expressed as follows:
where t is the time of movement of the charge material from the end of the chute to the charge level, v cx ,v cy ,v cz Respectively representing the flow of charge material out of the chute endsThe speed of movement of the end. The location of the drop point of the charge may be expressed as:
wherein r is xc ,r yc ,r zc Respectively, showing the position of movement of the charge material away from the end of the chute.
(7) Mathematical model for furnace charge movement locus
The charging materials at different positions of the throttling valve form different charging material movement tracks on the top of the blast furnace, and the set of the movement tracks of the material flow formed by all particles on the top of the blast furnace is called a material flow movement track mathematical model.
The modeling method of the furnace top burden movement track of the blast furnace provided by the embodiment of the invention comprises the steps of establishing a static coordinate system which is static relative to the blast furnace and a moving coordinate system which rotates together with a chute, establishing a burden throttling valve discharge mathematical model according to the discharge mode of burden at a blanking pot, establishing a free falling mathematical model of the burden at a central throat, establishing an initial movement mathematical model of the burden on the chute after the burden collides with the chute according to the collision recovery coefficient of the burden and the chute, establishing a chute movement mathematical model of the burden in the chute according to the initial movement mathematical model of the burden on the chute after the burden collides with the chute, and obtaining the burden movement track mathematical model of single-particle burden from the throttling valve to the tail end of the chute according to the chute movement mathematical model, thereby obtaining the movement track mathematical model of the whole blast furnace top burden and solving the technical problem of low calculation precision of the burden movement track of the blast furnace top burden of the existing blast furnace, the movement track of the furnace charge particles in different initial movement states at the top of the blast furnace is calculated by calculating the initial movement state of the single-particle furnace charge and the movement state of the rotary chute, so that the calculation precision of the movement track of the furnace charge is greatly improved, the function of a traditional material flow movement track model can be realized, meanwhile, the position and the speed distribution condition of the whole material flow at the tail end of the chute can be solved according to the movement tracks of the particles in different initial movement states, and the falling point position of the whole material flow at the charge level of the blast furnace can be further solved.
Further, the modeling method for the movement track of the furnace top of the serial tank type bell-less blast furnace based on coordinate transformation comprises a throttle valve discharge model, a central throat pipe descending model, a furnace material and chute collision model, a rotary chute movement model and a furnace throat area inclined throwing movement model, calculates the descending height of the furnace material at the central throat pipe based on a coordinate transformation matrix by establishing a static coordinate system and a dynamic coordinate system so as to further calculate the three-dimensional movement speed of the furnace material after collision with the chute, considers the influence of the initial movement state of the furnace material particles and the movement state of the chute on the movement track of the furnace material, can calculate the position distribution and the speed distribution condition of the material flow at the tail end of the chute, meets the distribution requirement of the serial tank type bell-less blast furnace, and has great application value.
EXAMPLE III
2650m of the third embodiment of the invention 3 The large-scale bell-less blast furnace is an experimental platform, a charging motion trajectory mathematical model is applied to the top of the serial tank type bell-less blast furnace to estimate and calculate the charging motion, and a charging motion trajectory calculation model shown in fig. 3 is constructed. Referring to fig. 7, specific steps of the embodiment for specifically completing the calculation of the movement trajectory model of the furnace top burden of the serial tank type bell-less blast furnace are as follows:
step 1: 2650m according to domestic 3 Specific physical parameters of the large-scale bell-less blast furnace, initializing a calculation model, wherein the calculation model comprises the following physical parameters: the diameter D of the central throat pipe, the height H of the central throat pipe, the tilting distance e of the chute, the radius R of the chute, the length L of the chute and the normal recovery coefficient e of the furnace burden when colliding with the chute n And tangential coefficient of restitution e t The sliding friction coefficient mu between the furnace burden and the chute and the like; initial movement state of charge particles: initial motion position and velocity; initial position of chute: the initial position of the furnace burden under the moving coordinate system and the static coordinate system; initial movement state of the chute: initial movement speeds under a moving coordinate system and a static coordinate system; initial position of charge: initial positions under a moving coordinate system and a static coordinate system; initial movement state of the charge: initial speeds under a moving coordinate system and a static coordinate system; the single-step operation time h of the furnace burden is set as 0 at the moment;
Step2：initial movement state of the input particles, i.e. initial movement position r of the particles at the throttle valve 0 Inputting the motion state parameters of the U-shaped chute, including the initial motion position (beta) of the chute when the particles leave the throttle valve 0 ,γ 0 ) Horizontal rotation angular velocity omega of chute 1 And angular acceleration a 1 Angular velocity ω of chute inclination rotation 2 And angular acceleration a 2 ；
Step 3: according to the opening size of the throttle valve, calculating different movement positions of the burden at the throttle valve { r } 1 ,r 2 ...r n N represents the number of charge particles at the choke. Initiating the position r of the movement of the particles leaving the throttle valve i Let i equal to 1;
step 3: calculating the initial movement speed v of the charge i when leaving the throttle valve according to the material flow discharge model 0 ；
Step 4: according to the initial movement state of charge i when leaving the throttle valve, including the initial position
And initial velocity
And the movement state of the chute, including the initial position (beta) 0 ,γ 0 ) And the rotational speed (omega) 1 ,ω 2 ) Calculating the time t consumed by the process that the charge particles leave the throttle valve and are placed in the chute collision process 0 Further, the total height h of the particle falling in the central throat is determined a +h b The speed of movement v of the particles at the point of impact 1 Position of collision point r 1 Position of chute (. beta.) 1 ,γ 1 )。
Step 5: according to the position of the charge when it collides with the chute
And velocity
And the position (beta) of the chute at the time of collision 1 ,γ 1 ) Calculate outPosition of furnace charge in moving coordinate system
And velocity
And according to the collision mathematical model of the furnace charge and the chute, calculating the position of the furnace charge after the collision of the furnace charge and the chute
And velocity
Step 6: by the position of the charge
And velocity
Solving the position of the furnace charge leaving the chute for the initial motion state of the furnace charge on the chute according to a motion mathematical model of the furnace charge on the chute
And velocity
And calculating the position of the furnace charge in the static coordinate system when the furnace charge leaves the chute according to the coordinate transformation matrix
And velocity
Step 7: by location of charge
And velocity
Solving the position of a falling point of the furnace burden on a charge level according to an inclined throwing motion model of the furnace burden in the dead zone for the initial motion state of the furnace burden in the throat dead zone.
Step 8: judging whether i is larger than n, if so, turning to Step 9; otherwise, go to Step 3;
step 9: and outputting the furnace burden at the position of the falling point of the blast furnace burden surface, and finishing.
Referring to fig. 8, the blast furnace top burden movement trajectory modeling system provided by the embodiment of the present invention includes:
the modeling method comprises a memory 10, a processor 20 and a computer program stored on the memory 10 and capable of running on the processor 20, wherein the processor 20 realizes the steps of the modeling method for the movement locus of the blast furnace top burden material provided by the embodiment when executing the computer program.
The specific working process and working principle of the modeling system for the furnace top furnace burden motion trail of the blast furnace in the embodiment can refer to the working process and working principle of the modeling method for the furnace top furnace burden motion trail of the blast furnace in the embodiment.
The present invention has been described in terms of the preferred embodiment, and it is not intended to be limited to the embodiment. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Claims (9)
1. A modeling method for a furnace top charging material movement track of a blast furnace is characterized by comprising the following steps:
establishing a static coordinate system which is static relative to the blast furnace and a moving coordinate system which rotates together with the chute;
establishing a furnace burden throttling valve discharge mathematical model according to the discharge mode of the furnace burden at the blanking tank;
establishing a free falling mathematical model of the furnace charge in a central throat, wherein the free falling mathematical model is used for acquiring the speed and the position of the furnace charge in a static coordinate system before the furnace charge collides with a chute;
according to a collision recovery coefficient of the furnace charge colliding with the chute, establishing an initial motion mathematical model of the furnace charge on the chute after the furnace charge collides with the chute, wherein the collision recovery coefficient comprises a normal recovery coefficient and a radial recovery coefficient;
according to an initial movement mathematical model of the furnace burden on the chute after the furnace burden collides with the chute, establishing a chute movement mathematical model of the furnace burden in the chute;
and acquiring a furnace charge movement track mathematical model of the single-particle furnace charge from the throttle valve to the tail end of the chute according to the chute movement mathematical model, and further acquiring the movement track mathematical model of the whole furnace top furnace charge of the blast furnace.
2. The method of claim 1, wherein establishing a static coordinate system that is stationary relative to the furnace and a moving coordinate system that rotates with the chute comprises:
establishing a static coordinate system which is static relative to the blast furnace according to the right-hand rule, wherein the origin of the static coordinate system is the intersection point between the horizontal connecting line of the two connecting points of the chute and the central throat pipe and the symmetrical axis of the blast furnace;
according to the right-hand rule, sequentially rotating the static coordinate system around two coordinate axes in the static coordinate system to obtain a dynamic coordinate system rotating together with the chute;
and solving a coordinate transformation matrix between the static coordinate system and the moving coordinate system.
3. The modeling method of the trajectory of the blast furnace top charge according to claim 1, wherein establishing a mathematical model of the free fall of the charge in the central throat comprises:
obtaining coordinates of the furnace charge in a static coordinate system when the furnace charge leaves the throttling valve according to the furnace charge throttling valve discharge mathematical model;
according to the coordinates of the furnace charge in a static coordinate system when leaving the throttling valve, a free fall mathematical model of the furnace charge in the central throat is established, and the specific calculation formula of the free fall mathematical model of the furnace charge in the central throat is as follows:
wherein r is 1 And v 1 Respectively the speed and the position of the furnace charge in a static coordinate system before the collision with the chute, x 0 、y 0 、h a Respectively the coordinates of the X-axis, the Y-axis and the Z-axis in the static coordinate system when the burden leaves the throttling valve, h w Is the distance between the tail end of the central throat pipe and the collision point of the furnace charge and the chute, v 0 The Z-axis speed of the furnace charge when leaving the throttle valve, g is the gravity acceleration, e is the chute tilting distance, beta 0 And gamma 0 Respectively the horizontal rotation angle and the inclination angle theta of the chute relative to the initial position when the furnace burden collides with the chute 0 Is the included angle between the furnace charge and the symmetrical axis of the chute, and
r is the chute radius.
4. The modeling method for the movement locus of the furnace top burden of the blast furnace as claimed in claim 1 or 3, wherein the step of establishing the initial movement mathematical model of the burden on the chute after the burden collides with the chute according to the collision recovery coefficient of the burden and the chute comprises the following steps:
the speed and the position of the furnace burden in a static coordinate system at the moment before the furnace burden collides with the chute are converted into the moving speed and the moving position in a moving coordinate system;
obtaining normal phase vectors of collision points on the chute according to the movement position of the furnace charge in the moving coordinate system before the furnace charge collides with the chute;
obtaining the incident speed of the furnace burden before the furnace burden collides with the chute according to the movement speed of the furnace burden before the furnace burden collides with the chute in the moving coordinate system;
calculating an included angle between a normal phase vector of a collision point on the chute and an incident speed of the furnace burden before the furnace burden collides with the chute to obtain a first included angle;
obtaining an initial speed constant of the furnace burden moving on the chute after the furnace burden collides with the chute according to the collision recovery coefficient of the furnace burden colliding with the chute and the first included angle;
and calculating the emergent speed of the furnace burden after the furnace burden collides with the chute according to the initial speed constant of the furnace burden moving on the chute after the furnace burden collides with the chute.
5. The modeling method for the motion trajectory of furnace top furnace burden of the blast furnace as defined in claim 4, wherein the specific formula for calculating the exit velocity of the furnace burden after the furnace burden collides with the chute is as follows:
wherein, v' 1 Is incident velocity v 'before furnace burden collides with chute' 2 The emergence speed of the furnace burden after the furnace burden collides with the chute, | | · | | represents the operation of mold length taking, e n Is the normal coefficient of restitution, theta, of the charge when it collides with the chute int Is a first angle theta out Is an included angle between a normal phase vector of a collision point on the chute and an emergent speed of the furnace burden after the collision with the chute, and
and n is a normal phase vector of a collision point on the chute.
6. The modeling method for the movement locus of the furnace top furnace burden of the blast furnace as claimed in claim 5, wherein the calculation formula for calculating the exit velocity of the furnace burden after the collision with the chute according to the initial velocity constant of the furnace burden moving on the chute after the collision with the chute is as follows:
v 2 ′＝av 1 ′+bn，
wherein, v' 1 Is incident velocity v 'before furnace burden collides with chute' 2 The speed is the emergent speed of the furnace burden after the furnace burden collides with the chute, a and b are initial speed constants of the furnace burden moving on the chute after the furnace burden collides with the chute, and n is a normal phase vector of a collision point on the chute.
7. The modeling method of the movement trajectory of the furnace top burden of the blast furnace as defined in claim 6, wherein the establishing of the mathematical model of the chute movement of the burden within the chute based on the mathematical model of the initial movement of the burden on the chute after the burden collides with the chute comprises:
obtaining an initial movement position and an initial movement speed of the furnace burden on the chute according to an initial movement mathematical model of the furnace burden on the chute after the furnace burden collides with the chute;
respectively analyzing the relative motion and the involved motion of the furnace burden relative to the chute based on the point composite motion according to the initial motion position and the initial motion speed of the furnace burden on the chute, further analyzing the absolute motion of the furnace burden in the blast furnace, and establishing a chute motion mathematical model of the furnace burden in the chute according to Newton's second law;
and obtaining the movement speed and the movement position of the furnace charge from the throttling valve to the tail end of the chute according to the chute movement mathematical model.
8. The modeling method of blast furnace top charge motion profile according to claim 7, further comprising after obtaining the tip end motion velocity and the tip end motion position of the charge from the choke valve to the chute tip end based on said mathematical model of chute motion:
establishing a furnace throat empty area movement mathematical model, wherein the calculation formula of the furnace throat empty area movement mathematical model is as follows:
wherein S is x 、S y And S z Represents the moving distance of the furnace charge on the X axis, the Y axis and the Z axis of the throat area, t is the moving time of the furnace charge from the tail end of the chute to the charge level, v cx ,v cy ,v cz Respectively representing the moving speed of the furnace burden leaving the tail end of the chute;
obtaining the position of a falling point of the furnace burden according to the furnace throat empty area movement mathematical model, wherein the calculation formula of the position of the falling point of the furnace burden is as follows:
wherein r is xc 、r yc And r zc Respectively, representing the position of movement of the charge material away from the end of the chute.
9. A blast furnace top charge motion profile modeling system, the system comprising:
memory (10), processor (20) and computer program stored on the memory (10) and executable on the processor (20), characterized in that the steps of the method according to any of the preceding claims 1 to 8 are implemented when the computer program is executed by the processor (20).
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210454570.2A CN114921598B (en) | 2022-04-27 | 2022-04-27 | Blast furnace top burden motion trail modeling method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210454570.2A CN114921598B (en) | 2022-04-27 | 2022-04-27 | Blast furnace top burden motion trail modeling method and system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114921598A true CN114921598A (en) | 2022-08-19 |
| CN114921598B CN114921598B (en) | 2023-05-09 |
Family
ID=82806472
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210454570.2A Active CN114921598B (en) | 2022-04-27 | 2022-04-27 | Blast furnace top burden motion trail modeling method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114921598B (en) |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105243190A (en) * | 2015-09-14 | 2016-01-13 | 内蒙古科技大学 | Method for modeling material surface output shape in material layout process in blast furnace |
| JP2016017197A (en) * | 2014-07-08 | 2016-02-01 | 新日鐵住金株式会社 | Discharge behavior estimation method and discharge behavior estimation system of raw material from furnace top bunker of blast furnace |
| CN106011350A (en) * | 2016-07-05 | 2016-10-12 | 东北大学 | Method and system for estimating charge surface shapes in blast furnace charge distribution process |
| CN110551861A (en) * | 2019-09-23 | 2019-12-10 | 中冶赛迪重庆信息技术有限公司 | method, system, equipment and storage medium for characterizing blast furnace burden surface shape |
| CN111666683A (en) * | 2020-06-05 | 2020-09-15 | 中冶赛迪重庆信息技术有限公司 | Method and system for predicting distribution state of furnace burden in blast furnace |
| CN112176136A (en) * | 2020-09-24 | 2021-01-05 | 中南大学 | Method and system for modeling movement locus of furnace charge on U-shaped chute of blast furnace |
| CN113656966A (en) * | 2021-08-17 | 2021-11-16 | 北京神网创新科技有限公司 | Blast furnace bell-less top on-line material distribution model simulation method |
| CN114139438A (en) * | 2021-10-26 | 2022-03-04 | 中冶南方工程技术有限公司 | Blast furnace burden track model construction method |
-
2022
- 2022-04-27 CN CN202210454570.2A patent/CN114921598B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2016017197A (en) * | 2014-07-08 | 2016-02-01 | 新日鐵住金株式会社 | Discharge behavior estimation method and discharge behavior estimation system of raw material from furnace top bunker of blast furnace |
| CN105243190A (en) * | 2015-09-14 | 2016-01-13 | 内蒙古科技大学 | Method for modeling material surface output shape in material layout process in blast furnace |
| CN106011350A (en) * | 2016-07-05 | 2016-10-12 | 东北大学 | Method and system for estimating charge surface shapes in blast furnace charge distribution process |
| CN110551861A (en) * | 2019-09-23 | 2019-12-10 | 中冶赛迪重庆信息技术有限公司 | method, system, equipment and storage medium for characterizing blast furnace burden surface shape |
| CN111666683A (en) * | 2020-06-05 | 2020-09-15 | 中冶赛迪重庆信息技术有限公司 | Method and system for predicting distribution state of furnace burden in blast furnace |
| CN112176136A (en) * | 2020-09-24 | 2021-01-05 | 中南大学 | Method and system for modeling movement locus of furnace charge on U-shaped chute of blast furnace |
| CN113656966A (en) * | 2021-08-17 | 2021-11-16 | 北京神网创新科技有限公司 | Blast furnace bell-less top on-line material distribution model simulation method |
| CN114139438A (en) * | 2021-10-26 | 2022-03-04 | 中冶南方工程技术有限公司 | Blast furnace burden track model construction method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114921598B (en) | 2023-05-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112176136B (en) | Method and system for modeling movement locus of furnace charge on U-shaped chute of blast furnace | |
| CN113139275B (en) | Blast furnace throat temperature estimation method based on multilayer ore-coke ratio distribution model | |
| Mio et al. | Analysis of traveling behavior of nut coke particles in bell-type charging process of blast furnace by using discrete element method | |
| Zhao et al. | Comprehensive mathematical model for particle flow and circumferential burden distribution in charging process of bell-less top blast furnace with parallel hoppers | |
| Hattori et al. | Development of burden distribution simulation model for bell-less top in a large blast furnace and its application | |
| Nag et al. | Development of material trajectory simulation model for blast furnace compact bell-less top | |
| Mitra | Modeling of burden distribution in the blast furnace | |
| CN114921598A (en) | Blast furnace top furnace burden movement track modeling method and system | |
| CN112226557B (en) | Blast furnace burden distribution process evaluation method and system, readable storage medium and application thereof | |
| JPH0586410A (en) | Method for controlling charge of raw material in furnace top of blast furnace | |
| CN113656966A (en) | Blast furnace bell-less top on-line material distribution model simulation method | |
| Ichida et al. | Radial distribution of burden descent velocity near burden surface in blast furnace | |
| Fu et al. | Burden descending and redistribution in a blast furnace | |
| JP6747232B2 (en) | Blast furnace raw material mixing ratio estimation method | |
| Zhou et al. | Motion trajectory mathematical model of burden flow at the top of bell-less blast furnace based on coordinate transformation | |
| JP2727563B2 (en) | Blast furnace operation method | |
| Chen et al. | DEM simulation of burden circumferential distribution of blast furnace with parallel hoppers | |
| EP3989013A1 (en) | Method for controlling process, operation guidance method, method for operating blast furnace, method for manufacturing molten iron, and device for controlling process | |
| Terui et al. | Estimation of particle segregation behavior in ore-coke mixed layer using screening layer model | |
| CN114065597A (en) | Mathematical method for representing iron ore reflow characteristics in blast furnace reflow zone | |
| JPH07179916A (en) | Bell-less type furnace top charging device for vertical furnace | |
| Bayat | Optimization of Lime Injection in EAF Systems by Using Pneumatic Sidewall Injectors | |
| Qiu et al. | DEM simulation of particle flow in a parallel-hopper bell-less blast furnace charging model | |
| CN118374755A (en) | Continuous tin smelting equipment and method based on immersed top-blown smelting technology | |
| JPH0627284B2 (en) | Blast furnace operation method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |