CN102982218B - Method for acquiring material flow characteristic at throttling valve on furnace top of blast furnace - Google Patents

Abstract

The invention discloses a method for acquiring a material flow characteristic at a throttling valve on the furnace top of a blast furnace, which comprises the following steps: by basis physical property parameters of furnace charges, calculating material flow Q of the furnace charges when the furnace charges pass through the throttling valve; by a curve fitting method, acquiring a corresponding relation between the opening degree R of the throttling valve and the material flow Q; and by the corresponding relation between the opening degree R of the throttling valve and the material flow Q and Q=M/t, obtaining the opening degree R of the throttling valve, wherein M represents a charge weight and t represents distribution time. The method for acquiring the material flow characteristic at the throttling valve on the furnace top of the blast furnace, which is provided by the invention, not only is beneficial for avoiding the condition that due to distribution segregation, the edge or the center is excessively heavy so as to influence the forward movement effect of the blast furnace, but also is beneficial for selecting a proper and high-quality throttling valve plate material and ensuring long service life and high efficiency use of distribution equipment.

Description

Obtain the method for blast furnace roof throttling valve place stream characteristic

Technical field

The present invention relates to technical field of blast furnace ironmaking, particularly a kind of method obtaining blast furnace roof throttling valve place stream characteristic.

Background technology

The throttling valve of blast furnace roof is the visual plant of roof material distributing, and its Main Function is the discharging speed controlling furnace charge in batch can, ensures that furnace charge flows out with rational speed and flow.Especially at many rings without in clock distributing technique, be ensure charge in official hour or the number of turns, the visual plant of the complete a collection of furnace charge of cloth.If throttling valve controls undesirable to stream, cloth segregation can be caused, cause edge or center overweight, thus affect smooth operation of furnace effect.Therefore, each factories and miness, all competitively by a large amount of commerical tests, determine the rational relation between throttle valve opening and stream amount.

But from the angle of production cost, commerical test certainly will consume a large amount of manpowers and material resources, and the test period is also longer.This is unfavorable for, when furnace charge physical property and granularmetric composition change, carrying out adjusting in time, fast to throttle valve opening and stream magnitude relation.And for without clock cloth, the real-time flowing velocity of throttling valve place furnace charge is also grasp the important parameter without clock cloth regular.And be difficult to obtain accurately furnace charge speed in real time by commerical test.In addition, from the angle of throttle plate selection, the stressing conditions of throttle plate in cloth process is still lacked and analyzes and understanding, be unfavorable for selecting throttle plate material that is appropriate, high-quality, ensure longevity, the effective use of apparatus for distributing.

Summary of the invention

Technical matters to be solved by this invention is to provide the relation of a kind of energy accurately, between quick obtaining throttle valve opening and stream amount, thus avoids cloth segregation, improves the method for the acquisition blast furnace roof throttling valve place stream characteristic of operation of blast furnace effect.

For solving the problems of the technologies described above, the invention provides a kind of method obtaining blast furnace roof throttling valve place stream characteristic, comprising:

Calculated by furnace charge basic physical properties parameter and obtain furnace charge by stream amount Q during throttling valve;

The corresponding relation between throttle valve opening R and stream amount Q is obtained by curve fitting method;

By the corresponding relation between throttle valve opening R and stream amount Q and calculate and obtain throttle valve opening R, wherein, M is batch weight, and t is the cloth time.

Further, described by furnace charge basic physical properties parameter calculate obtain furnace charge comprised by stream amount Q during throttling valve:

After measuring acquisition furnace charge basic physical properties parameter, then by calculating the normal force obtaining furnace charge and be subject to with tangential force F _{s, t};

According to described normal force with tangential force F _s,tthe component of both acquisitions on x, y, z tri-directions, then calculated by Newton second law and obtain the linear acceleration of charging movement on x, y, z direction and angular acceleration;

According to the linear acceleration of charging movement on x, y, z direction and the angular acceleration of charging movement, obtain the angular velocity component of t+ Δ t/2 moment furnace charge on x, y, z tri-directions and linear velocity component by center interpolation calculation;

Calculate according to the angular velocity component of t+ Δ t/2 moment furnace charge on x, y, z tri-directions and linear velocity component and obtain the displacement of furnace charge in the t+ Δ t/2 moment, thus calculate the displacement of not furnace charge in the same time;

The stream amount of furnace charge by throttling valve is obtained according to the misalignment calculating not furnace charge in the same time $Q=\frac{\underset{`1}{\overset{t}{\mathrm{\Σ}}}{Q}_t}{t}.$

Further, described furnace charge basic physical properties parameter comprises:

The density of furnace charge, particle diameter, the coefficient of sliding friction, coefficient of rolling friction, elastic modulus, Poisson ratio.

Further, the described misalignment according to calculating not furnace charge in the same time obtains the stream amount of furnace charge by throttling valve comprise:

According to the displacement of not furnace charge in the same time, remaining furnace charge amount P in batch can when obtaining each moment t _t, as known t ₁, t ₂furnace charge amount P in moment batch can _t1, P _t2time, the stream amount that just can obtain in the unit interval is $Q_1=\frac{P_{t2}-P_{t1}}{t_2-t_1};$

At Fixed Time Interval (t ₂-t ₁) when, calculate not stream amount Q in the same time ₁, Q ₂, Q ₃q _t, then obtain furnace charge when throttling valve is in a certain joint aperture by calculating by the average stream amount of throttling valve be

Further, the corresponding relation obtained between throttle valve opening R and stream amount Q by curve fitting method is comprised:

When being in a certain aperture according to described throttling valve, furnace charge is by the stream amount of throttling valve obtain furnace charge when throttling valve is in different opening and, by the stream amount of throttling valve, draw out curve map by calculating;

The curve fitting method pass obtained between throttle valve opening R and stream amount Q is adopted to be Q=x-yR+zR according to described curve map ², wherein, x, y, z is the constant that curve obtains.

Further, described by furnace charge basic physical properties parameter calculate obtain furnace charge also comprised by stream amount Q during throttling valve:

According to normal force calculate the stressing conditions in throttle plate unit area, when furnace charge and throttle plate contact, furnace charge can be subject to the extruding force from valve plate, and same furnace charge also can produce a direction on the contrary to valve plate, and size is reacting force, obtain the power that valve plate unit area is subject to by calculating be $P=\frac{F_{n,t}}{\mathrm{\π}{\left(\frac{d_{\mathrm{pi}}}{2}\right)}^2}.$

A kind of method obtaining blast furnace roof throttling valve place stream characteristic provided by the invention, first measures the basic physical properties parameter such as density, particle diameter, elastic modulus, Poisson ratio, friction factor of furnace charge by experiment.Again according to stressed to furnace charge of Newton interpolation algorithm-central difference method, movement velocity carries out iterative, obtain not stream amount in the same time, the relation between throttle valve opening and stream amount is obtained again according to curve fitting method analysis, finally, according to the concrete distribution of force situation in acting force-opposition relation determination throttling valve valve plate unit area, the relation between the throttle valve opening of arbitrary shape and stream amount accurately can be judged by aforesaid operations, the concrete stressing conditions of throttling valve valve plate and real-time stream speed, not only be conducive to avoiding cloth segregation, cause edge or center overweight and affect smooth operation of furnace effect, also help and select appropriately, the throttle plate material of high-quality, ensure the longevity of apparatus for distributing, effective use.

Accompanying drawing explanation

The concrete operation step process flow diagram of the method for the acquisition blast furnace roof throttling valve place stream characteristic that Fig. 1 provides for the embodiment of the present invention;

Fig. 2 is the curve synoptic diagram in the embodiment of the present invention between throttle valve opening and stream amount.

Embodiment

See Fig. 1, a kind of method obtaining blast furnace roof throttling valve place stream characteristic that the embodiment of the present invention provides, comprises following step:

Step S1: calculated by furnace charge basic physical properties parameter and obtain furnace charge by stream amount Q during throttling valve.

Step S2: calculate the corresponding relation obtained between throttle valve opening R and stream amount Q by curve fitting method and charge M and cloth time t.

Step S3: by the corresponding relation between throttle valve opening R and stream amount Q and calculate and obtain throttle valve opening R, wherein, M is batch weight, and t is the cloth time.

Wherein, step S1 is comprised by stream amount Q during throttling valve by furnace charge basic physical properties parameter calculating acquisition furnace charge:

Step S11: measure the density p obtaining furnace charge by the method first by Physical Experiment _i, particle diameter R _i, friction coefficient mu, coefficient of rolling friction μ _r, elastic modulus E, after the basic physical properties parameters such as Poisson ratio ν, then obtain the normal force that is subject to of furnace charge by calculating with tangential force F _s,t;

To the physical parameter that obtains be measured, and substitute into following formula and can calculate the power that furnace charge is subject to:

$F_{n,t}=\underset{j}{\mathrm{\Σ}}(2{(\frac{1-v_i^2}{E_i}+\frac{1-v_j^2}{E_j})}^{-1}\mathrm{\Δu}+2\sqrt{m2{(\frac{1-v_i^2}{E_i}+\frac{1-v_j^2}{E_j})}^{-1}}\frac{\mathrm{\Δu}}{\mathrm{\Δt}})$ （1-1）

$F_{s,t}=\underset{j}{\mathrm{\Σ}}(8{(\frac{2(2-v_i)(1+v_i)}{E_i}+\frac{2(2-v_j)(1+v_j)}{E_j})}^{-1}\mathrm{\Δu}+2\sqrt{m8{(\frac{2(2-v_i)(1+v_i)}{E_i}+\frac{2(2-v_j)(1+v_j)}{E_j})}^{-1}}\frac{\mathrm{\Δu}}{\mathrm{\Δt}})$

（1-2）

$\mathrm{\Δu}=|R_i+R_j-\sqrt{{(x_i-x_j)}^2+{(y_i-y_j)}^2+{(z_i-z_j)}^2}|$ （1-3）

$m=\frac{4}{3}\mathrm{\π}{R}_i^3$ （1-4）

F _s,t>μ|F _n，t| （1-5）

In formula

f _{s, t}---the normal direction that t furnace charge i is subject to and tangential force, N;

Δ u---overlapping amount between normal direction particle, m;

Δ t---time step, s;

V _i, v _j---the Poisson ratio of particle i, j;

E _i, E _j---the Young modulus of particle i, j, GPa;

R _i, R _j---the radius of particle i, j, m;

X _i, x _j, y _i, y _j, z _i, z _j---the locus of particle i, j;

The quality of m---particle, kg.

The coefficient of sliding friction of μ---charging movement.

Step S12: first according to described normal force with tangential force F _s,tthe component of both acquisitions on x, y, z tri-directions, then calculated by Newton second law and obtain the linear acceleration of charging movement on x, y, z direction and angular acceleration;

By its component on x, y, z tri-directions can be obtained

$\left[\begin{array}{c}{a}_{\mathrm{ix}}^{\left(t\right)}\\ a_{\mathrm{iy}}^{\left(t\right)}\\ a_{\mathrm{iz}}^{\left(t\right)}\end{array}\right]=\stackrel{\‾}{m_i}\left[\begin{array}{c}{F}_{n,x}^{\left(t\right)}\\ F_{n,y}^{\left(t\right)}\\ F_{n,z}^{\left(t\right)}\end{array}\right]$ （1-6）

$a_{\mathrm{\ω}}=\frac{\mathrm{\Σ}(-R_i\underset{j}{\mathrm{\Σ}}{F}_{s,t}-\frac{3}{8}{\mathrm{\μ}}_r|F_{n,t}\left|\right)}{\frac{8}{15}\mathrm{\π}{\mathrm{\ρ}}_i{R}_i^5}$ （1-7）

---the acceleration on x, y, z direction of charging movement, m/s ^-2;

A _ω---the angular acceleration of charging movement, Rad/s ²;

---particle i is in the normal force of t on x, y, z tri-directions;

M _i---the quality of particle, kg.

μ _r---the coefficient of rolling friction of particle;

R _i---the radius of particle, m;

π---circular constant, 3.1415926;

ρ _i---the density of particle, kg/m ³.

Step S13: according to the linear acceleration of charging movement on x, y, z direction and the angular acceleration of charging movement, obtains the angular velocity component of t+ Δ t/2 moment furnace charge on x, y, z tri-directions and linear velocity component by center interpolation calculation.

According to center method of interpolation, calculate the speed of t+ Δ t/2 moment particle i.

$\left[\begin{array}{c}{\stackrel{\·}{x}}_i^{(t+\mathrm{\Δt}/2)}\\ {\stackrel{\·}{y}}_i^{(t+\mathrm{\Δt}/2)}\\ {\stackrel{\·}{z}}_i^{(t+\mathrm{\Δt}/2)}\\ \end{array}\right]=\left[\begin{array}{c}{\stackrel{\·}{x}}_i^{(t-\mathrm{\Δt}/2)}\\ {\stackrel{\·}{y}}_i^{(t-\mathrm{\Δt}/2)}\\ {\stackrel{\·}{z}}_i^{(t-\mathrm{\Δt}/2)}\\ \end{array}\right]+\left[\begin{array}{c}{a}_{\mathrm{ix}}^{\left(t\right)}\\ a_{\mathrm{iy}}^{\left(t\right)}\\ a_{\mathrm{iz}}^{\left(t\right)}\end{array}\right]\mathrm{\Δt}$ （1-8）

（1-9）

$\sqrt{{\left(a_{\mathrm{\ωx}}^{\left(t\right)}\right)}^2+{\left(a_{\mathrm{\ωy}}^{\left(t\right)}\right)}^2+{\left(a_{\mathrm{\ωz}}^{\left(t\right)}\right)}^2}=a_{\mathrm{\ω}}$ （1-10）

In formula:

---particle i is engraved in x when t+ Δ t/2, the linear velocity component on y, z tri-directions;

---particle i is engraved in x when t-Δ t/2, the linear velocity component on y, z tri-directions;

---particle i is engraved in x when t+ Δ t/2, the angular velocity component on y, z tri-directions;

---particle i is engraved in x when t-Δ t/2, the angular velocity component on y, z tri-directions;

---the linear acceleration on x, y, z direction of charging movement, m/s ^-2;

A _{ω x} ^(t), a _{ω y} ^(t), a _{ω z} ^(t)---the angular acceleration of charging movement on x, y, z direction, Rad/s ^-2;

Δ t---time step, s;

M _i---the quality of particle i, kg;

Step S14: calculate according to the angular velocity component of t+ Δ t/2 moment furnace charge on x, y, z tri-directions and linear velocity component and obtain the displacement of furnace charge in the t+ Δ t/2 moment, thus calculate the displacement of not furnace charge in the same time.

The displacement of t+ Δ t unit can be obtained by the speed in t+ Δ t/2 moment.

$\left[\begin{array}{c}{x}_i^{(t+\mathrm{\Δt})}\\ y_i^{(t+\mathrm{\Δt})}\\ z_i^{(t+\mathrm{\Δt})}\end{array}\right]=\left[\begin{array}{c}{x}_i^{\left(t\right)}\\ y_i^{\left(t\right)}\\ z_i^{\left(t\right)}\end{array}\right]+\mathrm{\Δt}\left[\begin{array}{c}{\stackrel{\·}{x}}_i^{(t+\mathrm{\Δt}/2)}\\ {\stackrel{\·}{y}}_i^{(t+\mathrm{\Δt}/2)}\\ {\stackrel{\·}{z}}_i^{(t+\mathrm{\Δt}/2)}\\ \end{array}\right]$ （1-11）

（1-12）

---particle i is engraved in x when t+ Δ t/2, the speed component on y, z tri-directions;

---particle i is at the translational component of t+ Δ t on x, y, z tri-directions;

---particle i is at the translational component of t on x, y, z tri-directions.

---particle i is at the angular displacement component of t+ Δ t on x, y, z tri-directions;

---particle i is at the angular displacement component of t on x, y, z tri-directions;

Step S15: obtain the stream amount of furnace charge by throttling valve according to the misalignment calculating not furnace charge in the same time $Q=\frac{\underset{`1}{\overset{t}{\mathrm{\Σ}}}{Q}_t}{t}.$

According to the displacement of not furnace charge in the same time, remaining furnace charge amount P in batch can when obtaining each moment t _t, as known t ₁, t ₂furnace charge amount P in moment batch can _t1, P _t2time, the stream amount that just can obtain in the unit interval is Q ₁.

$Q_1=\frac{P_{t2}-P_{t1}}{t_2-t_1}$ （1-13）

At Fixed Time Interval (t ₂-t ₁) when, calculate a series of stream amount Q ₁, Q ₂, Q ₃q _t, then by calculating, to obtain when throttling valve is in a certain joint aperture furnace charge by the average stream amount of throttling valve be Q.

$Q=\frac{\underset{`1}{\overset{t}{\mathrm{\Σ}}}{Q}_t}{t}$ （1-14）

Therefore, can by the material flow rate calculation mode of furnace charge during a certain aperture by throttling valve, calculate and to obtain when throttling valve is in different opening furnace charge by the stream amount of throttling valve, according to stream amount Q corresponding in different opening (R) situation, draw the graph of a relation of discharging flow Q and aperture (R).

Step S16: according to the numerical density situation calculated in step S11, calculate the stressing conditions in throttle plate unit area.

When furnace charge and throttle plate contact, furnace charge can be subject to the extruding force from valve plate, and same furnace charge also can produce a direction on the contrary to valve plate, and size is reacting force, its computing formula is see step S11.By calculating the power obtaining and valve plate unit area is subject to:

$P=\frac{F_{n,t}}{\mathrm{\π}{\left(\frac{d_{\mathrm{pi}}}{2}\right)}^2}$ （1-15）

By calculating the stressing conditions in valve plate unit area, and then also just there is known the stressed of valve plate, being conducive to selecting throttle plate material that is appropriate, high-quality, ensureing longevity, the effective use of apparatus for distributing.

Step S2 calculates by curve fitting method and charge M and cloth time t the corresponding relation obtained between throttle valve opening R and stream amount Q and comprises:

Step S21: when being in same aperture according to described throttling valve, furnace charge is by the stream amount of throttling valve obtain furnace charge when throttling valve is in different opening and, by the stream amount of throttling valve, draw out curve map by calculating;

Step S22: adopt the curve fitting method pass obtained between throttle valve opening R and stream amount Q to be Q=x-yR+zR according to described curve map ², wherein, x, y, z is the constant that curve obtains.

Step S3 by corresponding relation between throttle valve opening R and stream amount Q and obtain throttle valve opening R, wherein M is batch weight, and t comprises the cloth time:

Step S31: in a practical situation, charge M and cloth time t known, just can calculate the material flow of furnace charge

Step S32: according to Q=x-yR+zR ²and obtain throttle valve opening R.

A kind of method obtaining blast furnace roof throttling valve place stream characteristic provided by the invention, first measures the basic physical properties parameter such as density, particle diameter, elastic modulus, Poisson ratio, friction factor of furnace charge by experiment.Again according to stressed to furnace charge of Newton interpolation algorithm-central difference method, movement velocity carries out iterative, obtain not stream amount in the same time, the relation between throttle valve opening and stream amount is obtained again according to curve fitting method analysis, finally, according to the concrete distribution of force situation in acting force-opposition relation determination throttling valve valve plate unit area, the relation between the throttle valve opening of arbitrary shape and stream amount accurately can be judged by aforesaid operations, the concrete stressing conditions of throttling valve valve plate and real-time stream speed, not only be conducive to avoiding cloth segregation, cause edge or center overweight and affect smooth operation of furnace effect, also help and select appropriately, the throttle plate material of high-quality, ensure the longevity of apparatus for distributing, effective use.

Below by specific embodiment, the present invention is further described.

Embodiment one:

Steps A 11: by measuring the size obtaining furnace charge basic physical properties parameter and batch can, then stove calculates by material basic physical properties parameter the stressing conditions F obtaining particle.

For coke by Physical Experiment can record coke basic physical properties as shown in Table 1, density (ρ), coefficient of rolling friction (μ _r), the coefficient of sliding friction (μ), Poisson ratio (ν), elastic modulus (E) and batch can size.

Table 1

The center first setting batch can exit is the initial point (0,0,0) of three-dimensional coordinate, and an existing position is in coke granule, the normal direction that it is subject to and the tangential force of (1.69,1.73,2.87) f _s,tbe respectively:

$F_{n,t}=286968,$ F _s，t＝526400

Steps A 12: the acceleration calculating particle according to Newton second law.

By the normal direction that obtains in steps A 11 and tangential force f _s,tdecompose x, acceleration of motion y, z tri-directions solving coke obtains:

$a_{\mathrm{ix}}^{\left(t\right)}=1006.1,$ $a_{\mathrm{iy}}^{\left(t\right)}=2487.7,$

a _ωx ^(t)＝286.2，a _ωy ^(t)＝1415，a _ωz ^(t)＝1543

Steps A 13: by calculating the speed obtaining particle.

According to the acceleration that steps A 12 obtains, solve line of motion speed and the angular velocity of coke according to central difference method, obtain:

${\stackrel{\·}{x}}_i^{(t+\mathrm{\Δt}/2)}=11.7,$ ${\stackrel{\·}{y}}_i^{(t+\mathrm{\Δt}/2)}=13.5,$ ${\stackrel{\·}{z}}_i^{(t+\mathrm{\Δt}/2)}=58.8$

Steps A 14: by calculating the displacement obtaining particle.

$x_i^{(t+\mathrm{\Δt})}=0.00017,$ $y_i^{(t+\mathrm{\Δt})}=0.0002,$ $z_i^{(t+\mathrm{\Δt})}=0.00088$

Steps A 15: according to calculating the furnace charge misalignment of gained, remaining furnace charge amount P in batch can when just can judge each moment t _t.So, as known t ₁, t ₂furnace charge amount P in moment batch can _t1, P _t2time, just can obtain the stream amount Q in the unit interval ₁.

$Q_1=\frac{P_{t2}-P_{t1}}{t_2-t_1}$

At Fixed Time Interval (t ₂-t ₁) when, just can obtain a series of stream amount Q ₁, Q ₂, Q ₃q _t, then the Q that averages.

$Q=\frac{\underset{`1}{\overset{t}{\mathrm{\Σ}}}{Q}_t}{t}$

Therefore, according to stream amount Q corresponding in different opening (R) situation, draw the graph of a relation of discharging flow Q and aperture (R).

Such as: work as t ₁when=0, P _t1=5609kg; Work as t ₂when=0.15, P _t2=5558kg;

Q ₁=340kg/s just can obtain a series of stream value in this way, then asks its average stream amount Q=347kg/s

Steps A 16: the normal force F being subject to valve plate effect according to the particle calculated in step S11 _{n, t}, according to the relation between acting force and reacting force, calculate the stressing conditions in throttle plate unit area, and then also just there is known the stressed of valve plate.

When steps A 2 is in same aperture according to described throttling valve, furnace charge is by the stream amount of throttling valve obtain furnace charge when throttling valve is in different opening and, by the stream amount of throttling valve, draw out curve map by calculating, as shown in Figure 2, adopt the method for curve to obtain the relational expression of mass velocity (Q) and throttle valve opening (R).

Q=2406.44-142.41*R+2.17257*R ²（1-17）

R=0.9646

As charge M and the cloth time, t was known when, just can calculate the mass rate (Q) of furnace charge:

$Q=\frac{M}{t}$ （1-18）

Simultaneous formula (1-17), (1-18) just can obtain the size of corresponding throttle valve opening (R).

It should be noted last that, above embodiment is only in order to illustrate technical scheme of the present invention and unrestricted, although with reference to example to invention has been detailed description, those of ordinary skill in the art is to be understood that, can modify to technical scheme of the present invention or equivalent replacement, and not departing from the spirit and scope of technical solution of the present invention, it all should be encompassed in the middle of right of the present invention.

Claims (2)

1. obtain a method for blast furnace roof throttling valve place stream characteristic, it is characterized in that, comprising:

Calculated by furnace charge basic physical properties parameter and obtain furnace charge by stream amount Q during throttling valve;

The corresponding relation between throttle valve opening R and stream amount Q is obtained by curve fitting method;

By the corresponding relation between throttle valve opening R and stream amount Q and calculate and obtain throttle valve opening R, wherein, M is batch weight, and t is the cloth time;

Wherein, described calculating by furnace charge basic physical properties parameter is obtained furnace charge and is comprised by stream amount Q during throttling valve: after measuring and obtaining furnace charge basic physical properties parameter, then by calculating the normal force F that acquisition furnace charge is subject to _n,twith tangential force F _s,t, " n, t " representation to, " s, t " represents tangential;

According to described normal force F _n,twith tangential force F _s,tthe component of both acquisitions on x, y, z tri-directions, then calculated by Newton second law and obtain the linear acceleration of charging movement on x, y, z direction and angular acceleration;

According to the linear acceleration of charging movement on x, y, z direction and the angular acceleration of charging movement, obtain the angular velocity component of t+ Δ t/2 moment furnace charge on x, y, z tri-directions and linear velocity component by center interpolation calculation;

Calculate according to the angular velocity component of t+ Δ t/2 moment furnace charge on x, y, z tri-directions and linear velocity component and obtain the displacement of furnace charge in the t+ Δ t/2 moment, thus calculate the displacement of not furnace charge in the same time, Δ t represents the time interval (t ₂-t ₁);

According to the displacement of not furnace charge in the same time, remaining furnace charge amount P in batch can when obtaining each moment t _t, as known t ₁, t ₂furnace charge amount P in moment batch can _t1, P _t2time, the stream amount that just can obtain in the unit interval is $Q_1=\frac{P_{t2}-P_{t1}}{t_2-t_1};$

At Fixed Time Interval (t ₂-t ₁) when, calculate not stream amount Q in the same time ₁, Q ₂, Q ₃q _t, then obtain furnace charge when throttling valve is in a certain joint aperture by calculating by the average stream amount of throttling valve be

Described furnace charge basic physical properties parameter comprises: the density of furnace charge, particle diameter, the coefficient of sliding friction, coefficient of rolling friction, elastic modulus, Poisson ratio;

The described corresponding relation obtained between throttle valve opening R and stream amount Q by curve fitting method is comprised: when being in a certain aperture according to described throttling valve, furnace charge is by the stream amount of throttling valve $Q=\frac{\underset{`1}{\overset{t}{\mathrm{\Σ}}}{Q}_t}{t},$ Obtain furnace charge when throttling valve is in different opening and, by the stream amount of throttling valve, draw out curve map by calculating;

The curve fitting method pass obtained between throttle valve opening R and stream amount Q is adopted to be Q=x-yR+zR according to described curve map ², wherein, x, y, z is the constant that curve obtains.

2. the as claimed in claim 1 method obtaining blast furnace roof throttling valve place stream characteristic, is characterized in that, described calculating by furnace charge basic physical properties parameter is obtained furnace charge and also comprised by stream amount Q during throttling valve:

According to normal force F _n,t, calculate the stressing conditions in throttle plate unit area, when furnace charge and throttle plate contact, furnace charge can be subject to the extruding force from valve plate, and same furnace charge also can produce a direction on the contrary to valve plate, and size is F _n,treacting force, obtain the power that valve plate unit area is subject to by calculating be d _pifor the particle diameter of furnace charge.