CN102982218A - 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, 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 material properties of flow

Technical field

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

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 of furnace charge in the control batch can, guarantees that furnace charge flows out with rational speed and flow.Especially in encircling without the clock distributing technique, be to guarantee charge in official hour or the number of turns, the visual plant of the complete a collection of furnace charge of cloth more.If throttling valve is undesirable to the material current control, can cause the cloth segregation, cause edge or center overweight, thereby affect the smooth operation of furnace effect.Therefore, each factories and miness are determined the rational relation between throttle valve opening and the material flow all competitively by a large amount of commerical tests.

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 size composition change, and throttle valve opening and material discharge relation are carried out in time, adjust fast.And for without the clock cloth, the real-time flowing velocity of throttling valve place furnace charge also is the important parameter of grasping without the clock cloth regular.And be difficult to obtain accurately in real time furnace charge speed by commerical test.In addition, from the angle of throttle plate selection, still lack for the stressing conditions of throttle plate in the cloth process and to analyze and understanding, be unfavorable for selecting throttle plate material appropriate, high-quality, guarantee longevity, the effective use of apparatus for distributing.

Summary of the invention

Technical matters to be solved by this invention provide a kind of energy accurately, the relation between quick obtaining throttle valve opening and material flow, thereby avoid the cloth segregation, improve the method that properties of flow is expected at blast furnace roof throttling valve place of obtaining of operation of blast furnace effect.

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

Material flow Q when obtaining furnace charge by throttling valve by furnace charge basic physical properties calculation of parameter;

By the corresponding relation between curve fitting method acquisition throttle valve opening R and the material flow Q;

Reach by the corresponding relation between throttle valve opening R and the material flow Q

Calculate and obtain throttle valve opening R, wherein, M is batch weight, and t is the cloth time.

Further, described material flow Q when obtaining furnace charge by throttling valve by furnace charge basic physical properties calculation of parameter comprises:

After measuring acquisition furnace charge basic physical properties parameter, obtain the normal force that furnace charge is subject to by calculating again With tangential force F _{S, t}

According to described normal force

With tangential force F _{S, t}Obtain both at x, y, the component on three directions of z calculates by Newton second law and obtains charging movement at x, y, the linear acceleration on the z direction and angular acceleration;

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

Furnace charge is at x constantly according to t+ Δ t/2, and y, the angular velocity component on three directions of z and linear velocity component calculate and obtain furnace charge in t+ Δ t/2 displacement constantly, thereby calculate the different constantly displacements of furnace charge;

Obtain furnace charge by the material flow of throttling valve according to the displacement situation of calculating different constantly furnace charges $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, described displacement situation according to calculating different constantly furnace charges obtains furnace charge by the material flow of throttling valve

Comprise:

According to the constantly displacement of furnace charge of difference, remaining furnace charge amount P in the batch can when obtaining each moment t _t, as known t ₁, t ₂Furnace charge amount P in the moment batch can _T1, P _T2The time, the material flow 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 ₁) situation under, calculate different material flow Q constantly ₁, Q ₂, Q ₃Q _t, then be in the average material flow by throttling valve of furnace charge in the situation of a certain joint aperture by calculate obtaining throttling valve and be

Further, the corresponding relation that obtains between throttle valve opening R and the material flow Q by curve fitting method comprises:

Be in the situation of a certain aperture furnace charge according to described throttling valve by the material flow of throttling valve

Be in the situation of different opening furnace charge by calculate obtaining throttling valve and draw out curve map by the material flow of throttling valve;

The pass of adopting curve fitting method to get between throttle valve opening R and the material flow Q according to described curve map is Q=x-yR+zR ², wherein, x, y, z are the constant that curve obtains.

Further, described material flow Q when obtaining furnace charge by throttling valve by furnace charge basic physical properties calculation of parameter also comprises:

According to normal force

Calculate the stressing conditions on the throttle plate unit area, when furnace charge and throttle plate contacted, furnace charge can be subject to the extruding force from valve plate, and same furnace charge also can produce an opposite direction to valve plate, and size is

Reacting force, by the power that calculate to obtain to be subject on the valve plate unit area be $P=\frac{F_{n,t}}{\mathrm{\π}{\left(\frac{d_{\mathrm{pi}}}{2}\right)}^2}.$

A kind of method of obtaining blast furnace roof throttling valve place material properties of flow provided by the invention is measured first the basic physical properties parameters such as density, particle diameter, elastic modulus, Poisson ratio, friction factor of furnace charge by experiment.Stressed to furnace charge according to newton's second law of motion-center differential technique again, movement velocity is carried out iterative, obtain different material flows constantly, analysis obtains throttle valve opening and the relation of material between flow according to curve fitting method again, at last, determine the concrete distribution of force situation on the throttling valve valve plate unit area according to acting force-opposition relation, can accurately judge the throttle valve opening of arbitrary shape and the relation between the material flow by aforesaid operations, the concrete stressing conditions of throttling valve valve plate and real-time material Flow Velocity, not only be conducive to avoid the cloth segregation, cause edge or center overweight and affect the smooth operation of furnace effect, also help and select appropriately, the throttle plate material of high-quality, the longevity of assurance apparatus for distributing, effective use.

Description of drawings

The concrete operation step process flow diagram of the method for obtaining blast furnace roof throttling valve place material properties of flow that Fig. 1 provides for the embodiment of the invention;

Fig. 2 is throttle valve opening and the curve synoptic diagram of material between flow in the embodiment of the invention.

Embodiment

Referring to Fig. 1, a kind of method of obtaining blast furnace roof throttling valve place material properties of flow that the embodiment of the invention provides comprises following step:

Step S1: the material flow Q when obtaining furnace charge by throttling valve by furnace charge basic physical properties calculation of parameter.

Step S2: by the corresponding relation between curve fitting method and charge M and cloth time t calculating acquisition throttle valve opening R and the material flow Q.

Step S3: reach by the corresponding relation between throttle valve opening R and the material flow Q

Calculate and obtain throttle valve opening R, wherein, M is batch weight, and t is the cloth time.

Material flow Q when wherein, step S1 obtains furnace charge by throttling valve by furnace charge basic physical properties calculation of parameter comprises:

Step S11: by measure first the density p of acquisition furnace charge by the method for 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 ν, the normal force that is subject to by calculate obtaining furnace charge again

With tangential force F _{S, t}

With the physical parameter that measures, the following formula of substitution 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 the formula

F _{S, t}---normal direction and tangential force that t moment furnace charge i is subject to, N;

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

Δ t---time step, s;

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

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

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

x _i, x _j, y _i, y _j, z _i, z _j---particle i, the locus of 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, t}Obtain both at x, y, the component on three directions of z calculates by Newton second law and obtains charging movement at x, y, the linear acceleration on the z direction and angular acceleration;

By

Can obtain it at x, y, the component on three directions of z

$\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）

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

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

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

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 charging movement at x, y, the linear acceleration on the z direction and the angular acceleration of charging movement, by the center interpolation calculation obtain t+ Δ t/2 constantly furnace charge at x, y, angular velocity component and linear velocity component on three directions of z.

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

$\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 the formula:

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

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

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

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

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

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

Δ t---time step, s;

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

Step S14: furnace charge is at x constantly according to t+ Δ t/2, and y, the angular velocity component on three directions of z and linear velocity component calculate and obtain furnace charge in t+ Δ t/2 displacement constantly, thereby calculate the different constantly displacements of furnace charge.

Can get the constantly displacement of unit of t+ Δ t by t+ Δ t/2 speed constantly.

$\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, y, the speed component on three directions of z;

---particle i is engraved in x when t+ Δ t, y, the translational component on three directions of z;

---particle i is engraved in x when t, y, the translational component on three directions of z.

---particle i is engraved in x when t+ Δ t, y, the angular displacement component on three directions of z;

---particle i is engraved in x when t, y, the angular displacement component on three directions of z;

Step S15: obtain furnace charge by the material flow of throttling valve according to the displacement situation of calculating different constantly furnace charges $Q=\frac{\underset{`1}{\overset{t}{\mathrm{\Σ}}}{Q}_t}{t}.$

According to the constantly displacement of furnace charge of difference, remaining furnace charge amount P in the batch can when obtaining each moment t _t, as known t ₁, t ₂Furnace charge amount P in the moment batch can _T1, P _T2The time, the material flow 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 ₁) situation under, calculate a series of material flow Q ₁, Q ₂, Q ₃Q _t, then being in the average material flow by throttling valve of furnace charge in the situation of a certain joint aperture by calculate obtaining throttling valve is Q.

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

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

Step S16: the particle stressing conditions according to calculating among the step S11 calculates the stressing conditions on the throttle plate unit area.

When furnace charge and throttle plate contacted, furnace charge can be subject to the extruding force from valve plate, and same furnace charge also can produce an opposite direction to valve plate, and size is

Reacting force, its computing formula is referring to step S11.By calculating the power that is subject on the acquisition valve plate unit area:

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

By calculating the stressing conditions on the valve plate unit area, so also stressed with regard to known valve plate, be conducive to select throttle plate material appropriate, high-quality, guarantee longevity, the effective use of apparatus for distributing.

Step S2 comprises by the corresponding relation between curve fitting method and charge M and cloth time t calculating acquisition throttle valve opening R and the material flow Q:

Step S21: be in the situation of same aperture furnace charge according to described throttling valve by the material flow of throttling valve

Be in the situation of different opening furnace charge by calculate obtaining throttling valve and draw out curve map by the material flow of throttling valve;

Step S22: the pass of adopting curve fitting method to get between throttle valve opening R and the material flow Q according to described curve map is Q=x-yR+zR ², wherein, x, y, z are the constant that curve obtains.

Step S3 reaches by the corresponding relation between throttle valve opening R and the material flow Q

Obtain throttle valve opening R, wherein M is batch weight, and t comprises the cloth time:

Step S31: in actual conditions, charge M and cloth time t are 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 of obtaining blast furnace roof throttling valve place material properties of flow provided by the invention is measured first the basic physical properties parameters such as density, particle diameter, elastic modulus, Poisson ratio, friction factor of furnace charge by experiment.Stressed to furnace charge according to newton's second law of motion-center differential technique again, movement velocity is carried out iterative, obtain different material flows constantly, analysis obtains throttle valve opening and the relation of material between flow according to curve fitting method again, at last, determine the concrete distribution of force situation on the throttling valve valve plate unit area according to acting force-opposition relation, can accurately judge the throttle valve opening of arbitrary shape and the relation between the material flow by aforesaid operations, the concrete stressing conditions of throttling valve valve plate and real-time material Flow Velocity, not only be conducive to avoid the cloth segregation, cause edge or center overweight and affect the smooth operation of furnace effect, also help and select appropriately, the throttle plate material of high-quality, the longevity of assurance apparatus for distributing, effective use.

Below by specific embodiment the present invention is further specified.

Embodiment one:

Steps A 11: obtain the size of furnace charge basic physical properties parameter and batch can by measurement, then stove obtains the stressing conditions F of particle by material basic physical properties calculation of parameter.

As an example of coke example 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 of setting first the batch can exit is the initial point of three-dimensional coordinate (0,0,0), and an existing position is in the coke granule of (1.69,1.73,2.87), the normal direction that it is subject to and tangential force

F _{S, t}Be respectively:

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

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

With normal direction and the tangential force that obtains in the steps A 11 F _{S, t}Decompose x, y, the acceleration of motion of finding the solution coke on three directions of z 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 that obtains particle.

The acceleration that obtains according to steps A 12 is found the solution line of motion speed and the angular velocity of coke according to the center differential technique, obtains:

${\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 that obtains 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 the furnace charge displacement situation of calculating gained, remaining furnace charge amount P in the batch can in the time of just can judging each moment t _tSo, as known t ₁, t ₂Furnace charge amount P in the moment batch can _T1, P _T2The time, just can obtain the material flow Q in the unit interval ₁

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

At Fixed Time Interval (t ₂-t ₁) situation under, just can obtain a series of material flow Q ₁, Q ₂, Q ₃Q _t, Q again averages.

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

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

For example: work as t ₁=0 o'clock, P _T1=5609kg; Work as t ₂=0.15 o'clock, P _T2=5558kg;

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

Steps A 16: the normal force F that is subject to the valve plate effect according to the particle that calculates among the step S11 _{N, t},, according to the relation between acting force and reacting force, calculate the stressing conditions on the throttle plate unit area, and then also stressed with regard to known valve plate.

Steps A 2 is in the situation of same aperture furnace charge by the material flow of throttling valve according to described throttling valve

Be in the situation of different opening furnace charge by calculate obtaining throttling valve and draw out curve map by the material flow of throttling valve, 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

In the known situation of charge M and cloth time t, 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 unrestricted in order to technical scheme of the present invention to be described, although with reference to example the present invention is had been described in detail, those of ordinary skill in the art is to be understood that, can make amendment or be equal to replacement technical scheme of the present invention, and not breaking away from the spirit and scope of technical solution of the present invention, it all should be encompassed in the middle of the claim scope of the present invention.

Claims (6)

1. a method of obtaining blast furnace roof throttling valve place material properties of flow is characterized in that, comprising:

Material flow Q when obtaining furnace charge by throttling valve by furnace charge basic physical properties calculation of parameter;

By the corresponding relation between curve fitting method acquisition throttle valve opening R and the material flow Q;

Reach by the corresponding relation between throttle valve opening R and the material flow Q

Calculate and obtain throttle valve opening R, wherein, M is batch weight, and t is the cloth time.

2. the method for obtaining blast furnace roof throttling valve place material properties of flow as claimed in claim 1 is characterized in that, described material flow Q when obtaining furnace charge by throttling valve by furnace charge basic physical properties calculation of parameter comprises:

After measuring acquisition furnace charge basic physical properties parameter, obtain the normal force that furnace charge is subject to by calculating again

With tangential force F _{S, t}

According to described normal force

With tangential force F _{S, t}Obtain both at x, y, the component on three directions of z calculates by Newton second law and obtains charging movement at x, y, the linear acceleration on the z direction and angular acceleration;

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

Furnace charge is at x constantly according to t+ Δ t/2, and y, the angular velocity component on three directions of z and linear velocity component calculate and obtain furnace charge in t+ Δ t/2 displacement constantly, thereby calculate the different constantly displacements of furnace charge;

Obtain furnace charge by the material flow of throttling valve according to the displacement situation of calculating different constantly furnace charges $Q=\frac{\underset{`1}{\overset{t}{\mathrm{\Σ}}}{Q}_t}{t}.$

3. the method for obtaining blast furnace roof throttling valve place material properties of flow as claimed in claim 2 is characterized in that 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.

4. the method for obtaining blast furnace roof throttling valve place material properties of flow as claimed in claim 3 is characterized in that, described displacement situation according to calculating different constantly furnace charges obtains furnace charge by the material flow of throttling valve $Q=\frac{\underset{`1}{\overset{t}{\mathrm{\Σ}}}{Q}_t}{t}$ Comprise:

According to the constantly displacement of furnace charge of difference, remaining furnace charge amount P in the batch can when obtaining each moment t _t, as known t ₁, t ₂Furnace charge amount P in the moment batch can _T1, P _T2The time, the material flow 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 ₁) situation under, calculate different material flow Q constantly ₁, Q ₂, Q ₃Q _t, then be in the average material flow by throttling valve of furnace charge in the situation of a certain joint aperture by calculate obtaining throttling valve and be

5. the method for obtaining blast furnace roof throttling valve place material properties of flow as claimed in claim 4 is characterized in that, the described corresponding relation that obtains between throttle valve opening R and the material flow Q by curve fitting method comprises:

Be in the situation of a certain aperture furnace charge according to described throttling valve by the material flow of throttling valve

Be in the situation of different opening furnace charge by calculate obtaining throttling valve and draw out curve map by the material flow of throttling valve;

The pass of adopting curve fitting method to get between throttle valve opening R and the material flow Q according to described curve map is Q=x-yR+zR ², wherein, x, y, z are the constant that curve obtains.

6. the method for obtaining blast furnace roof throttling valve place material properties of flow as claimed in claim 3 is characterized in that, described material flow Q when obtaining furnace charge by throttling valve by furnace charge basic physical properties calculation of parameter also comprises:

According to normal force

Calculate the stressing conditions on the throttle plate unit area, when furnace charge and throttle plate contacted, furnace charge can be subject to the extruding force from valve plate, and same furnace charge also can produce an opposite direction to valve plate, and size is

Reacting force, by the power that calculate to obtain to be subject on the valve plate unit area be $P=\frac{F_{n,t}}{\mathrm{\π}{\left(\frac{d_{\mathrm{pi}}}{2}\right)}^2}.$

Images