CN105243190A - Method for modeling material surface output shape in material layout process in blast furnace - Google Patents

Abstract

The present application discloses a method for modeling a material surface output shape in a material layout process in blast furnace. The method comprises: setting a chute inclination vector and a material layout laps vector of a material layout matrix; setting a function of the material surface output shape in the material layout process; taking a horizontal plane of the material bottom output shape as a reference to construct an integral constraint of the function of the material surface output shape; setting a base function of a unique lap material surface output of the material layout to construct a material layout output function with different ring bits; and determining a final material surface output shape according to the material layout output function with a last ring bit. According to the method of the present application, an entire outline of the material surface shape output by the material layout matrix can be conveniently and accurately obtained, and image visualization of the material surface output shape in the blast furnace material layout matrix can be realized, so as to reduce a fuel consumption ratio, save energy consumption, ensure stable and smooth in blast furnace conditions, stable in blast furnace yield, and prolonged in blast furnace life, avoid failures such as stuffiness, difficult operation, collapse, material collapse, and the like, and facilitate realization of dynamic distribution control of the blast furnace material output shape, and automation of the whole smelting process of the blast furnace.

Description

A kind of blast furnace material distribution process charge level exports the modeling method of shape

Technical field

The invention belongs to blast furnace ironmaking process control field, particularly relate to the modeling method that a kind of blast furnace material distribution process charge level exports shape.

Background technology

Blast furnace material distribution is a kind of operational means of blast furnace ironmaking procedure regulation furnace charge in the distribution of furnace throat, and furnace charge is burden distribution shape whether reasonable in blase furnace cast iron output of furnace throat, and fuel consumption, the stable smooth operation of prolonging campaign and the working of a furnace has great impact.Burden distribution matrix is the dependent variable that charge level exports shape, refers to the rule that distributing device is followed in cloth process, is made up of according to a certain order chute angle and cloth number of turns sequence.The charge level output shape of blast furnace material distribution process is a distribution function be made up of internal state parameter in airtight smelting environment, but not the data message of single-point.In view of blast furnace material distribution process charge level exports shape to the key effect of blast furnace stable smooth operation, and existing charge level exports the disappearance of shape modeling technological means, build the model of blast furnace material distribution the output of process charge level function from the angle of process control, the control exporting shape for blast furnace material distribution process charge level provides input/output model to be still the difficult problem that in blast furnace ironmaking process, is not captured.

Publication No. is that patent document discloses of CN104133945A a kind ofly establishes the relationship model of blast furnace material distribution controling parameters to charge level, for description blast furnace material distribution model, blanking process model has certain positive role, but significantly shortcoming needs to there are more material surface stock rod height detection data, and the accuracy of charge level altitude information is had higher requirements, and in reality, charge level real time scan detection technique is also a difficult problem.Meanwhile, it have ignored cloth number of rings in burden distribution matrix and charge level is exported to the impact of shape, and under not considering permanent batch weight, dynamic charge level exports the integral constraint of shape simultaneously.Therefore, existing detection technique can not penetrate the scanning distributed intelligence that furnace wall obtains charge level, the output shape of blast furnace material distribution matrix detects by stock rod all the time, can not obtain the overall picture of whole burden distribution shape, and be unfavorable for the control of cloth the output of process shape of charge level.

Summary of the invention

The object of the present invention is to provide a kind of blast furnace material distribution process charge level to export the modeling method of shape, reduce modeling process to the requirement detecting data, reduce fuel ratio, easily and accurately can obtain the overall picture of the whole shape of charge level that burden distribution matrix exports.

In order to solve the problems of the technologies described above, the invention discloses the modeling method that a kind of blast furnace material distribution exports shape of charge level, it is characterized in that, comprise the steps:

Step 1: the chute dip vector α and the cloth number of turns vector κ that build burden distribution matrix

α＝[α ₁,α ₂,…,α _m]∈R ^ma≤α _i≤b，

κ＝[κ ₁,κ ₂,…,κ _m]∈N ^m,

A and b represents the minimum and maximal value of chute tilt adjustable scope respectively, and m represents the total number of rings of cloth, k _ibe i-th be periphery material time in inclination alpha _ithe number of turns that lower chute rotates, k _afor total cloth number of turns;

Step 2: according to described chute dip vector α and cloth number of turns vector κ, the output shape of charge level algorithm arranging cloth process is:

γ(y)＝f(y,α,κ)，0≤y≤r，

0 and r distribution represent blast furnace center and furnace wall, border;

Step 3: using the bottom shape of charge level output for surface level is as with reference to benchmark, build the integral constraint exporting shape of charge level function

$V_a={\∫}_0^{r}2\mathrm{\π}y\mathrm{\γ}\left(y\right)dy={\∫}_0^{r}2\mathrm{\π}yf(y,\mathrm{\α},\mathrm{\κ})dy;$

Step 4: the charge volume building individual pen cloth is:

V _u＝V _a/κ _a；

Step 5: the basis function arranging individual pen cloth output charge level is:

$u(y,\mathrm{\α},\mathrm{\σ})=\mathrm{\ξ}\mathrm{\σ}\exp (-\frac{{(y-g(\mathrm{\α}\left)\right)}^2}{2{\mathrm{\σ}}^2});$

G (α) is drop point function corresponding to cloth inclination angle, and ξ is corrected parameter, and variable σ is form parameter;

Step 6: the integral constraint building unit basis function

$F(\mathrm{\α},\mathrm{\σ})={\∫}_0^{r}2\mathrm{\π}y\mathrm{\σ}\exp (-\frac{{(y-g(\mathrm{\α}\left)\right)}^2}{2{\mathrm{\σ}}^2})dy=V_u;$

Step 7: according to the iterative solution of the form parameter σ of the integral constraint solution procedure 5 of step 6;

Step 8: build first ring cloth output function by burden distribution matrix:

${\mathrm{\γ}}_{l_1}(y,{\mathrm{\α}}_1,{\mathrm{\σ}}_1)={\mathrm{\σ}}_1\exp (-\frac{{(y-g({\mathrm{\α}}_1\left)\right)}^2}{2{\mathrm{\σ}}_1^2});$

Step 9: by burden distribution matrix build the second ring to m be periphery material output function be:

$\begin{array}{c}{\mathrm{\γ}}_{l_2}\left(y\right)=f(y,\{{\mathrm{\α}}_1,{\mathrm{\κ}}_1\},\{{\mathrm{\α}}_2,{\mathrm{\κ}}_2\left\}\right)\\ \begin{array}{c}.\\ .\\ .\end{array}\\ {\mathrm{\γ}}_{l_m}\left(y\right)=f(y,\{{\mathrm{\α}}_1,{\mathrm{\κ}}_1\},\{{\mathrm{\α}}_2,{\mathrm{\κ}}_2\},...,\{{\mathrm{\α}}_m,{\mathrm{\κ}}_m\left\}\right)\end{array};$

Step 10: determine that final charge level exports shape and is according to the cloth output function of last ring position:

$\mathrm{\γ}\left(y\right)=f(y,\mathrm{\α},\mathrm{\κ})={\mathrm{\γ}}_{l_m}\left(y\right)=f(y,\{{\mathrm{\α}}_1,{\mathrm{\κ}}_1\},\{{\mathrm{\α}}_2,{\mathrm{\κ}}_2\},...,\{{\mathrm{\α}}_m,{\mathrm{\κ}}_m\left\}\right).$

Compared with prior art, the present invention has following outstanding substantive distinguishing features:

Method is simple, reduce modeling process to the requirement detecting data number, the control exporting shape for charge level provides computation model, comprise chute dip vector and the cloth number of turns vector of setting burden distribution matrix, the output shape of charge level function of cloth process is set, using the bottom shape of charge level output for surface level is as reference benchmark, build the integral constraint exporting shape of charge level function, the basis function of charge level is exported by arranging individual pen cloth, build the cloth output function of different rings position, and determine that final charge level exports shape according to the cloth output function of last ring position, easily and accurately can obtain the overall picture of the whole shape of charge level that burden distribution matrix exports, achieve the graph visualization of the output shape of charge level of blast furnace material distribution matrix, reduce fuel ratio, save energy consumption, ensure working of a furnace stable smooth operation, blast furnace stable yields, blast furnace is lengthened the life, blast furnace is avoided to suppress wind, be difficult to walk, cave in and collapse the appearance of the faults such as material, be conducive to realizing blast furnace burden and export the DYNAMIC DISTRIBUTION control of shape and even the robotization of whole blast furnace ironmaking process.

Accompanying drawing explanation

Accompanying drawing described herein is used to provide a further understanding of the present invention, forms a part of the present invention, and schematic description and description of the present invention, for explaining the present invention, does not form inappropriate limitation of the present invention.In the accompanying drawings:

Fig. 1 is the schematic diagram of the output shape of charge level in blast furnace material distribution of the present invention;

Fig. 2 is discrete α angle of the present invention and the relativity of corresponding drop point x result of calculation;

The unit of Fig. 3 corresponding to dynamic deformation parameter of the present invention exports the relativity figure that basis function and definite value unit export basis function;

Fig. 4 is the comparison diagram that blast furnace material distribution of the present invention exports charge level function monocycle cloth;

Fig. 5 is the comparison diagram that blast furnace material distribution of the present invention exports charge level function monocycle cloth and polycyclic distributing;

Fig. 6 is the comparison diagram of two polycyclic distributings of blast furnace material distribution Output matrix charge level function of the present invention.

Embodiment

Blast furnace material distribution process charge level of the present invention exports the modeling method of shape, comprises the steps:

Refer to Fig. 1 to Fig. 6, blast furnace material distribution process parameter data is as following table:

The volume V of charge batch weight a	30m 3
		Stream initial velocity v 0	1.2m/s
Cloth period of time T a	75s
		The length H of Y-piece or passage y	1.8m
Collision loss coefficient η	0.8
		Chute length l o	4.0m
Coefficientoffrictionμ	0.5
		Gravity acceleration g	9.8m/s 2
Stockline degree of depth h	1.0m
		Furnace charge decline resistance F m	0.2mg
Chute anglec of rotation speed ω	0.84rad/s
		Chute inclination maximum b	45 degree
Chute inclination minimum a	10 degree

According to blast furnace material distribution process parameter data, by chute inclination angle in the equivalent discretize of the interval precision with 0.01 of a and b, calculated the drop point x of α respectively by blast furnace material distribution equation of locus, be the distance at furnace charge drop point distance blast furnace center, α is the cloth angle of chute:

$v_1=\sqrt{2{\mathrm{gH}}_y+v_0^2};$

v ₂＝ηv ₁cosα；

$v_3=\sqrt{\begin{array}{c}2{\mathrm{gl}}_0\left(\cos \mathrm{\α}-\mathrm{\μ}\sin \mathrm{\α}\right)\\ +4{\mathrm{\π}}^2{\mathrm{\ω}}^2{l}_0^{2}sin\mathrm{\α}{\left(\sin \mathrm{\α}+\mathrm{\μ}\cos \mathrm{\α}\right)}^2\\ +v_2^2\end{array}};$

$x_y=\frac{{\mathrm{mv}}_3^2{\sin }^2\mathrm{\α}}{mg-F_m}\left(\sqrt{\frac{1}{{\tan }^2\mathrm{\α}}+\frac{2mg-F_m}{{\mathrm{mv}}_3^2{\sin }^2\mathrm{\α}}\left(l_0\left(1-\cos \mathrm{\α}\right)+h\right)}-\frac{1}{\tan \mathrm{\α}}\right);$

x＝g(α)＝l ₀sinα+x _y。

By chute inclination angle after the equivalent discretize of the interval precision with 0.01 of a and b, drop point x discrete data result corresponding to the α angle of discretize as shown in Figure 2.

1 monocycle cloth example:

During monocycle cloth, m=1, the chute dip vector α and the cloth number of turns vector κ that form burden distribution matrix are scalar form, the arbitrary value within the scope of chute tilt adjustable can be got in chute inclination angle, and for the ease of simulation comparison, α gets 15,20,25,30,7 correlative values such as 35,40,45 grades, cloth number of rings was determined by the cloth time cycle, due to cloth time cycle 75s, often enclose 7.5 seconds cloth time, therefore total cloth number of turns k _a=10.

Provide the output shape of charge level function of monocycle cloth process

γ(y)＝f(y,α,κ)，0≤y≤r，

0 and r distribution represent blast furnace center and furnace wall, border.

Assuming that γ _ey () ≡ 0, builds the integral constraint exporting shape of charge level function

$V_a={\∫}_0^{r}2\mathrm{\π}y\mathrm{\γ}\left(y\right)dy={\∫}_0^{r}2\mathrm{\π}yf(y,\mathrm{\α},\mathrm{\κ})dy,$

Determine the charge volume of individual pen cloth: V _u=V _a/ κ _a=3m ³

Provide the basis function that individual pen cloth exports charge level

$u(y,\mathrm{\α},\mathrm{\σ})=\mathrm{\ξ}\mathrm{\σ}\exp (-\frac{{(y-g(\mathrm{\α}\left)\right)}^2}{2{\mathrm{\σ}}^2}),$

G (α) is calculated by cloth equation of locus, corrected parameter ξ=1;

Build the integral constraint of unit basis function

$F(\mathrm{\α},\mathrm{\σ})={\∫}_0^{r}2\mathrm{\π}y\mathrm{\σ}\exp (-\frac{{(y-g(\mathrm{\α}\left)\right)}^2}{2{\mathrm{\σ}}^2})dy=V_u=3,$

Solved the iterative solution of corresponding form parameter σ according to integral constraint by Gauss's Saden that alternative manner; α gets 15, and 20,25,30,35, during 40,45 7 correlative values such as grade, the contrast of the calculation by computer of basis function and definite value form parameter σ=0.2 is as shown in Figure 3.

Because monocycle cloth first ring is also last ring, therefore charge level output function:

$\mathrm{\γ}(y,\mathrm{\α},\mathrm{\σ})=\mathrm{\ξ}\mathrm{\σ}\exp (-\frac{{(y-g(\mathrm{\α}\left)\right)}^2}{2{\mathrm{\σ}}^2});$

Build the constraint condition of charge level output function:

$F(\mathrm{\α},\mathrm{\σ})={\∫}_0^{r}2\mathrm{\π}y\mathrm{\γ}\left(y,\mathrm{\α},\mathrm{\σ}\right)dy={\mathrm{\κ}}_a{V}_u=V_a=30m^3,$

α is solved 15,20,25,30,35,40, the iterative solution of the form parameter σ under 45 grade 7 correlative values according to integral constraint;

During monocycle cloth, charge level exports shape function,

γ(y)＝f(y,α,κ)＝γ(y,α,σ)，

The comparing result of its computer simulation data as shown in Figure 4.

The comparison example of 2 polycyclic distributing Computer Simulations and monocycle:

Monocycle charge level output function relation calculates, and this link emphasis calculates the charge level output function of polycyclic distributing;

By given burden distribution matrix α=[44.5,39.2,37.0,33.4,29.4], κ=[3,4,1,1,1], known m=5, cloth number of total coils k _a=3+4+1+1+1=10;

Provide many rings to be periphery the output shape of charge level function of material process

γ(y)＝f(y,α,κ)，0≤y≤r，

Assuming that γ _ey () ≡ 0, builds the integral constraint exporting shape of charge level function

$V_a={\∫}_0^{r}2\mathrm{\π}y\mathrm{\γ}\left(y\right)dy={\∫}_0^{r}2\mathrm{\π}yf(y,\mathrm{\α},\mathrm{\κ})dy,$

Determine the charge volume of individual pen cloth: V _u=V _a/ κ _a=3m ³

Provide the basis function that individual pen cloth exports charge level

$u(y,\mathrm{\α},\mathrm{\σ})=\mathrm{\ξ}\mathrm{\σ}\exp (-\frac{{(y-g(\mathrm{\α}\left)\right)}^2}{2{\mathrm{\σ}}^2}),$

G (α) is calculated by cloth equation of locus, corrected parameter ξ=1;

Build the integral constraint of unit basis function

$F(\mathrm{\α},\mathrm{\σ})={\∫}_0^{r}2\mathrm{\π}y\mathrm{\σ}\exp (-\frac{{(y-g(\mathrm{\α}\left)\right)}^2}{2{\mathrm{\σ}}^2})dy=V_u=3,$

Solved the iterative solution of corresponding form parameter σ according to described integral constraint by Gauss's Saden that iterative learning method; α gets 15, and 20,25,30,35, during 40,45 7 correlative values such as grade, the contrast of the calculation by computer of basis function and definite value form parameter σ=0.2 is as shown in Figure 3.

First ring cloth output function is built by burden distribution matrix:

${\mathrm{\γ}}_{l_1}(y,{\mathrm{\α}}_1,{\mathrm{\σ}}_1)={\mathrm{\σ}}_1\exp (-\frac{{(y-g({\mathrm{\α}}_1\left)\right)}^2}{2{\mathrm{\σ}}_1^2}),$

The constraint condition of first ring cloth is built by burden distribution matrix:

$F({\mathrm{\α}}_1,{\mathrm{\σ}}_1)={\∫}_0^{r}2{\mathrm{\πy\γ}}_{l_1}\left(y,{\mathrm{\α}}_1,{\mathrm{\σ}}_1\right)dy={\mathrm{\κ}}_1{V}_u=9m^3,$

Form parameter σ is solved according to integral constraint ₁iterative solution;

Build the second ring by burden distribution matrix to be periphery material output function to the 5th:

$\begin{array}{c}{\mathrm{\γ}}_{l_2}\left(y\right)=f(y,\{{\mathrm{\α}}_1,{\mathrm{\κ}}_1\},\{{\mathrm{\α}}_2,{\mathrm{\κ}}_2\left\}\right)\\ \begin{array}{c}.\\ .\\ .\end{array}\\ {\mathrm{\γ}}_{l_5}\left(y\right)=f(y,\{{\mathrm{\α}}_1,{\mathrm{\κ}}_1\},\{{\mathrm{\α}}_2,{\mathrm{\κ}}_2\},...,\{{\mathrm{\α}}_5,{\mathrm{\κ}}_{m5}\left\}\right)\end{array};$

By burden distribution matrix build the second ring to the 5th be periphery material constraint condition:

${\∫}_0^{r}2\mathrm{\π}y\left({\mathrm{\γ}}_{l_2}\right(y)-{\mathrm{\γ}}_{l_1}(y\left)\right)dy={\mathrm{\κ}}_2{V}_u=12m^3,$

${\∫}_0^{r}2\mathrm{\π}y\left({\mathrm{\γ}}_{l_3}\right(y)-{\mathrm{\γ}}_{l_2}(y\left)\right)dy={\mathrm{\κ}}_3{V}_u=3m^3,$

${\∫}_0^{r}2\mathrm{\π}y\left({\mathrm{\γ}}_{l_4}\right(y)-{\mathrm{\γ}}_{l_3}(y\left)\right)dy={\mathrm{\κ}}_4{V}_u=3m^3,$

${\∫}_0^{r}2\mathrm{\π}y\left({\mathrm{\γ}}_{l_5}\right(y)-{\mathrm{\γ}}_{l_4}(y\left)\right)dy={\mathrm{\κ}}_5{V}_u=3m^3,$

Row write the piecewise function of i ring:

${\mathrm{\γ}}_{l_i}\left(y\right)=\left\{\begin{array}{cc}{\mathrm{\σ}}_i\exp \left(-\frac{{\left(y-g\left({\mathrm{\α}}_i\right)\right)}^2}{2{\mathrm{\σ}}_i^2}\right),& {\mathrm{\γ}}_{l_i}\left(y\right)\≥{\mathrm{\γ}}_{l_{i-1}}\left(y\right)\\ {\mathrm{\γ}}_{l_{i-1}}\left(y\right),& {\mathrm{\γ}}_{l_{i-1}}\left(y\right)\≥{\mathrm{\γ}}_{l_i}\left(y\right)\end{array}\right.,$

Form parameter σ is solved respectively according to integral constraint _iiterative solution;

It is determined by the output shape of charge level function of last ring position that row write final charge level output shape function:

$\mathrm{\γ}\left(y\right)=f(y,\mathrm{\α},\mathrm{\κ})={\mathrm{\γ}}_{l_5}\left(y\right)=f(y,\{{\mathrm{\α}}_1,{\mathrm{\κ}}_1\},\{{\mathrm{\α}}_2,{\mathrm{\κ}}_2\},...,\{{\mathrm{\α}}_5,{\mathrm{\κ}}_5\left\}\right),$

Five rings cloth and the charge level of monocycle cloth when chute inclination alpha=35 ° export the comparing result of shape computer emulated data as shown in Figure 5.

The case of 3 two polycyclic distributing contrasts:

The burden distribution matrix of two contrasts is respectively,

Burden distribution matrix 1: α ₁=[37,34,32,28,21], κ ₁=[3,2,2,2,1];

Burden distribution matrix 2: α ₂=[42.5,40.0,37.5,34.5,31.5], κ ₂=[3,3,2,1,1];

3.1 charge levels calculating burden distribution matrix 1 export shape;

By given two known m=5 of burden distribution matrix, cloth number of total coils k _a=3+2+2+2+1=10;

Provide many rings to be periphery the output shape of charge level function of material process

γ(y)＝f(y,α,κ)，0≤y≤r，

Assuming that γ _ey () ≡ 0, builds the integral constraint exporting shape of charge level function

$V_a={\∫}_0^{r}2\mathrm{\π}y\mathrm{\γ}\left(y\right)dy={\∫}_0^{r}2\mathrm{\π}yf(y,\mathrm{\α},\mathrm{\κ})dy,$

Determine the charge volume of individual pen cloth: V _u=V _a/ κ _a=3m ³

Provide the basis function that individual pen cloth exports charge level

$u(y,\mathrm{\α},\mathrm{\σ})=\mathrm{\ξ}\mathrm{\σ}\exp (-\frac{{(y-g(\mathrm{\α}\left)\right)}^2}{2{\mathrm{\σ}}^2}),$

G (α) is calculated by cloth equation of locus, corrected parameter ξ=1;

Build the integral constraint of unit basis function

$F(\mathrm{\α},\mathrm{\σ})={\∫}_0^{r}2\mathrm{\π}y\mathrm{\σ}\exp (-\frac{{(y-g(\mathrm{\α}\left)\right)}^2}{2{\mathrm{\σ}}^2})dy=V_u=3,$

Solved the iterative solution of corresponding form parameter σ according to described integral constraint by Gauss's Saden that alternative manner; α gets 15, and 20,25,30,35, during 40,45 7 correlative values such as grade, the contrast of the calculation by computer of basis function and definite value form parameter σ=0.2 is as shown in Figure 3.

First ring cloth output function is built by burden distribution matrix:

${\mathrm{\γ}}_{l_1}(y,{\mathrm{\α}}_1,{\mathrm{\σ}}_1)={\mathrm{\σ}}_1\exp (-\frac{{(y-g({\mathrm{\α}}_1\left)\right)}^2}{2{\mathrm{\σ}}_1^2}),$

The constraint condition of first ring cloth is built by burden distribution matrix:

$F({\mathrm{\α}}_1,{\mathrm{\σ}}_1)={\∫}_0^{r}2{\mathrm{\πy\γ}}_{l_1}(y,{\mathrm{\α}}_1,{\mathrm{\σ}}_1)dy={\mathrm{\κ}}_1{V}_u=9m^3,$

Form parameter σ is solved according to integral constraint ₁iterative solution;

Build the second ring by burden distribution matrix to be periphery material output function to the 5th:

$\begin{array}{c}{\mathrm{\γ}}_{l_2}\left(y\right)=f(y,\{{\mathrm{\α}}_1,{\mathrm{\κ}}_1\},\{{\mathrm{\α}}_2,{\mathrm{\κ}}_2\left\}\right)\\ \begin{array}{c}.\\ .\\ .\end{array}\\ {\mathrm{\γ}}_{l_5}\left(y\right)=f(y,\{{\mathrm{\α}}_1,{\mathrm{\κ}}_1\},\{{\mathrm{\α}}_2,{\mathrm{\κ}}_2\},...,\{{\mathrm{\α}}_5,{\mathrm{\κ}}_{m5}\left\}\right)\end{array};$

By burden distribution matrix build the second ring to the 5th be periphery material constraint condition:

${\∫}_0^{r}2\mathrm{\π}y\left({\mathrm{\γ}}_{l_2}\right(y)-{\mathrm{\γ}}_{l_1}(y\left)\right)dy={\mathrm{\κ}}_2{V}_u=6m^3,$

${\∫}_0^{r}2\mathrm{\π}y\left({\mathrm{\γ}}_{l_3}\right(y)-{\mathrm{\γ}}_{l_2}(y\left)\right)dy={\mathrm{\κ}}_3{V}_u=6m^3,$

${\∫}_0^{r}2\mathrm{\π}y\left({\mathrm{\γ}}_{l_4}\right(y)-{\mathrm{\γ}}_{l_3}(y\left)\right)dy={\mathrm{\κ}}_4{V}_u=6m^3,$

${\∫}_0^{r}2\mathrm{\π}y\left({\mathrm{\γ}}_{l_5}\right(y)-{\mathrm{\γ}}_{l_4}(y\left)\right)dy={\mathrm{\κ}}_5{V}_u=3m^3,$

Row write the piecewise function of i ring:

${\mathrm{\γ}}_{l_i}\left(y\right)=\left\{\begin{array}{cc}{\mathrm{\σ}}_i\exp \left(-\frac{{\left(y-g\left({\mathrm{\α}}_i\right)\right)}^2}{2{\mathrm{\σ}}_i^2}\right),& {\mathrm{\γ}}_{l_i}\left(y\right)\≥{\mathrm{\γ}}_{l_{i-1}}\left(y\right)\\ {\mathrm{\γ}}_{l_{i-1}}\left(y\right),& {\mathrm{\γ}}_{l_{i-1}}\left(y\right)\≥{\mathrm{\γ}}_{l_i}\left(y\right)\end{array}\right.,$

Form parameter σ is solved respectively according to integral constraint _iiterative solution;

It is determined by the output shape of charge level function of last ring position that row write final charge level output shape function:

$\mathrm{\γ}\left(y\right)=f(y,\mathrm{\α},\mathrm{\κ})={\mathrm{\γ}}_{l_5}\left(y\right)=f(y,\{{\mathrm{\α}}_1,{\mathrm{\κ}}_1\},\{{\mathrm{\α}}_2,{\mathrm{\κ}}_2\},...,\{{\mathrm{\α}}_5,{\mathrm{\κ}}_5\left\}\right),$

Do by the burden distribution matrix of monocycle cloth the charge level inputted and export the comparing result of shape computer emulated data as shown in Figure 5.

3.2 calculating burden distribution matrixes are α ₂=[42.5,40.0,37.5,34.5,31.5], κ ₂during=[3,3,2,1,1], the charge level of polycyclic distributing exports shape;

By given two known m=5 of burden distribution matrix, cloth number of total coils k _a=3+3+2+1+1=10;

Provide many rings to be periphery the output shape of charge level function of material process

γ(y)＝f(y,α,κ)，0≤y≤r，

Assuming that γ _ey () ≡ 0, builds the integral constraint exporting shape of charge level function

$V_a={\∫}_0^{r}2\mathrm{\π}y\mathrm{\γ}\left(y\right)dy={\∫}_0^{r}2\mathrm{\π}yf(y,\mathrm{\α},\mathrm{\κ})dy,$

Determine the charge volume of individual pen cloth: V _u=V _a/ κ _a=3m ³

Provide the basis function that individual pen cloth exports charge level

$u(y,\mathrm{\α},\mathrm{\σ})=\mathrm{\ξ}\mathrm{\σ}\exp (-\frac{{(y-g(\mathrm{\α}\left)\right)}^2}{2{\mathrm{\σ}}^2}),$

G (α) is calculated by cloth equation of locus, corrected parameter ξ=1;

Build the integral constraint of unit basis function

$F(\mathrm{\α},\mathrm{\σ})={\∫}_0^{r}2\mathrm{\π}y\mathrm{\σ}\exp (-\frac{{(y-g(\mathrm{\α}\left)\right)}^2}{2{\mathrm{\σ}}^2})dy=V_u=3,$

Solved the iterative solution of corresponding form parameter σ according to described integral constraint by Gauss's Saden that alternative manner; α gets 15, and 20,25,30,35, during 40,45 7 correlative values such as grade, the contrast of the calculation by computer of basis function and definite value form parameter σ=0.2 is as shown in Figure 3.

First ring cloth output function is built by burden distribution matrix:

${\mathrm{\γ}}_{l_1}(y,{\mathrm{\α}}_1,{\mathrm{\σ}}_1)={\mathrm{\σ}}_1\exp (-\frac{{(y-g({\mathrm{\α}}_1\left)\right)}^2}{2{\mathrm{\σ}}_1^2}),$

The constraint condition of first ring cloth is built by burden distribution matrix:

$F({\mathrm{\α}}_1,{\mathrm{\σ}}_1)={\∫}_0^{r}2{\mathrm{\πy\γ}}_{l_1}(y,{\mathrm{\α}}_1,{\mathrm{\σ}}_1)dy={\mathrm{\κ}}_1{V}_u=9m^3,$

Form parameter σ is solved according to integral constraint ₁iterative solution;

Build the second ring by burden distribution matrix to be periphery material output function to the 5th:

$\begin{array}{c}{\mathrm{\γ}}_{l_2}\left(y\right)=f(y,\{{\mathrm{\α}}_1,{\mathrm{\κ}}_1\},\{{\mathrm{\α}}_2,{\mathrm{\κ}}_2\left\}\right)\\ \begin{array}{c}.\\ .\\ .\end{array}\\ {\mathrm{\γ}}_{l_5}\left(y\right)=f(y,\{{\mathrm{\α}}_1,{\mathrm{\κ}}_1\},\{{\mathrm{\α}}_2,{\mathrm{\κ}}_2\},...,\{{\mathrm{\α}}_5,{\mathrm{\κ}}_{m5}\left\}\right)\end{array};$

By burden distribution matrix build the second ring to the 5th be periphery material constraint condition:

${\∫}_0^{r}2\mathrm{\π}y\left({\mathrm{\γ}}_{l_2}\right(y)-{\mathrm{\γ}}_{l_1}(y\left)\right)dy={\mathrm{\κ}}_2{V}_u=9m^3,$

${\∫}_0^{r}2\mathrm{\π}y\left({\mathrm{\γ}}_{l_3}\right(y)-{\mathrm{\γ}}_{l_2}(y\left)\right)dy={\mathrm{\κ}}_3{V}_u=6m^3,$

${\∫}_0^{r}2\mathrm{\π}y\left({\mathrm{\γ}}_{l_4}\right(y)-{\mathrm{\γ}}_{l_3}(y\left)\right)dy={\mathrm{\κ}}_4{V}_u=3m^3,$

${\∫}_0^{r}2\mathrm{\π}y\left({\mathrm{\γ}}_{l_5}\right(y)-{\mathrm{\γ}}_{l_4}(y\left)\right)dy={\mathrm{\κ}}_5{V}_u=3m^3,$

Row write the piecewise function of i ring:

${\mathrm{\γ}}_{l_i}\left(y\right)=\left\{\begin{array}{cc}{\mathrm{\σ}}_i\exp \left(-\frac{{\left(y-g\left({\mathrm{\α}}_i\right)\right)}^2}{2{\mathrm{\σ}}_i^2}\right),& {\mathrm{\γ}}_{l_i}\left(y\right)\≥{\mathrm{\γ}}_{l_{i-1}}\left(y\right)\\ {\mathrm{\γ}}_{l_{i-1}}\left(y\right),& {\mathrm{\γ}}_{l_{i-1}}\left(y\right)\≥{\mathrm{\γ}}_{l_i}\left(y\right)\end{array}\right.,$

Form parameter σ is solved respectively according to integral constraint _iiterative solution;

It is determined by the output shape of charge level function of last ring position that row write final charge level output shape function:

$\mathrm{\γ}\left(y\right)=f(y,\mathrm{\α},\mathrm{\κ})={\mathrm{\γ}}_{l_5}\left(y\right)=f(y,\{{\mathrm{\α}}_1,{\mathrm{\κ}}_1\},\{{\mathrm{\α}}_2,{\mathrm{\κ}}_2\},...,\{{\mathrm{\α}}_5,{\mathrm{\κ}}_5\left\}\right),$

Do by the burden distribution matrix of polycyclic distributing the charge level inputted and export the comparing result of shape computer emulated data as shown in Figure 6.

Test findings: the modeling method that blast furnace material distribution process charge level of the present invention exports shape reduces modeling process to the requirement detecting data, and the output shape achieving burden distribution matrix is visual, is conducive to formulation and the adjustment of burden distribution system.Fuel ratio is a principal economic indicators of iron and steel process, represent the fuel sum often smelted one ton of iron and consume, it is an important indicator of a measurement blast furnace energy consumption height, certain ironworks fuel ratio in March 529.1kg, fuel ratio 530kg, revised significantly and regulated burden distribution matrix from May under technological guidance of the present invention April, and May, fuel ratio was 518, energy consumption reduces by 2.2%, significantly achieves energy-saving and cost-reducing effect.This patent obtains the overall picture of the output shape of charge level of burden distribution matrix easily and accurately, achieve the graph visualization of the output shape of charge level of blast furnace material distribution matrix, guarantee working of a furnace stable smooth operation, blast furnace stable yields, blast furnace are lengthened the life, avoid blast furnace to suppress wind, be difficult to walk, cave in and collapse the appearance of the faults such as material, be conducive to realizing blast furnace burden and export the DYNAMIC DISTRIBUTION control of shape and even the robotization of whole blast furnace ironmaking process.

Above-mentioned explanation illustrate and describes some calculated examples of the application, but as previously mentioned, be to be understood that the application is not limited to the form disclosed by this paper, should not regard the eliminating to other embodiments as, and can be used for other combinations various, amendment and environment, and can in application contemplated scope described herein, changed by the technology of above-mentioned instruction or association area or knowledge.And the change that those skilled in the art carry out and change do not depart from the spirit and scope of the application, then all should in the protection domain of the application's claims.

Claims (9)

1. blast furnace material distribution process charge level exports a modeling method for shape, it is characterized in that, comprises the steps:

Step 1: the chute dip vector α and the cloth number of turns vector κ that build burden distribution matrix

α＝[α ₁,α ₂,…,α _m]∈R ^m,a≤α _i≤b，

κ＝[κ ₁,κ ₂,…,κ _m]∈N ^m,

A and b represents the minimum and maximal value of chute tilt adjustable scope respectively, and m represents the total number of rings of cloth, k _ibe i-th be periphery material time in inclination alpha _ithe number of turns that lower chute rotates, k _afor total cloth number of turns;

Step 2: according to described chute dip vector α and cloth number of turns vector κ, the output shape of charge level algorithm arranging cloth process is:

γ(y)＝f(y,α,κ)，0≤y≤r；

0 and r distribution represent blast furnace center and furnace wall, border;

Step 3: using the bottom shape of charge level output for surface level is as with reference to benchmark, build the integral constraint exporting shape of charge level function

$V_a={\∫}_0^{r}2\mathrm{\π}y\mathrm{\γ}\left(y\right)dy={\∫}_0^{r}2\mathrm{\π}yf(y,\mathrm{\α},\mathrm{\κ})dy;$

Step 4: the charge volume building individual pen cloth is:

V _u＝V _a/κ _a；

Step 5: the basis function arranging individual pen cloth output charge level is:

$u\left(y,\mathrm{\α},\mathrm{\σ}\right)=\mathrm{\ξ}\mathrm{\σ}\exp \left(-\frac{{\left(y-g\left(\mathrm{\α}\right)\right)}^2}{2{\mathrm{\σ}}^2}\right);$

G (α) is drop point function corresponding to cloth inclination angle, and ξ is corrected parameter, and variable σ is form parameter;

Step 6: the integral constraint building unit basis function

$F(\mathrm{\α},\mathrm{\σ})={\∫}_0^{r}2\mathrm{\π}y\mathrm{\σ}\exp \left(-\frac{{(y-g(\mathrm{\α}\left)\right)}^2}{2{\mathrm{\σ}}^2}\right)dy=V_u;$

Step 7: according to the iterative solution of the form parameter σ of the integral constraint solution procedure 5 of step 6;

Step 8: build first ring cloth output function by burden distribution matrix:

${\mathrm{\λ}}_{l_1}(y,{\mathrm{\α}}_1,{\mathrm{\σ}}_1)={\mathrm{\σ}}_1\exp \left(-\frac{{(y-g({\mathrm{\α}}_1\left)\right)}^2}{2{\mathrm{\σ}}_1^2}\right);$

Step 9: by burden distribution matrix build the second ring to m be periphery material output function be:

$\begin{array}{c}{\mathrm{\γ}}_{l_2}\left(y\right)=f\left(y,\left\{{\mathrm{\α}}_1,{\mathrm{\κ}}_1\right\},\left\{{\mathrm{\α}}_2,{\mathrm{\κ}}_2\right\}\right)\\ .\\ .\\ .\\ {\mathrm{\γ}}_{l_m}\left(y\right)=f\left(y,\left\{{\mathrm{\α}}_1,{\mathrm{\κ}}_1\right\},\left\{{\mathrm{\α}}_2,{\mathrm{\κ}}_2\right\},...\left\{{\mathrm{\α}}_m,{\mathrm{\κ}}_m\right\}\right)\end{array};$

Step 10: determine that final charge level exports shape and is according to the cloth output function of last ring position:

$\mathrm{\γ}\left(y\right)=f(y,\mathrm{\α},\mathrm{\κ})={\mathrm{\γ}}_{l_m}\left(y\right)=f(y,\{{\mathrm{\α}}_1,{\mathrm{\κ}}_1\},\{{\mathrm{\α}}_2,{\mathrm{\κ}}_2\},...,\{{\mathrm{\α}}_m,{\mathrm{\κ}}_m\left\}\right).$

2. blast furnace material distribution process charge level exports the modeling method of shape according to claim 1, it is characterized in that, described burden distribution matrix is the input of cloth process dynamics, and it is that system exports that charge level exports shape.

3. blast furnace material distribution process charge level exports the modeling method of shape according to claim 1, it is characterized in that, the charge volume in the volume integral that the shaped upper part of described output charge level and bottom shape are formed and feed bin keeps constant.

4. blast furnace material distribution process charge level exports the modeling method of shape according to claim 3, it is characterized in that, the iterative solution of described form parameter σ adopts your alternative manner of Gauss's Saden to solve.

5. blast furnace material distribution process charge level exports the modeling method of shape according to claim 3, and it is characterized in that, described step 8 comprises:

Step 8.1: the constraint condition being built first ring cloth by burden distribution matrix is:

$F({\mathrm{\α}}_1,{\mathrm{\σ}}_1)={\∫}_0^{r}2{\mathrm{\πy\γ}}_{l_1}(y,{\mathrm{\α}}_1,{\mathrm{\σ}}_1)dy={\mathrm{\κ}}_1{V}_u$

Step 8.2: the integral constraint according to step 8.1 solves form parameter σ ₁iterative solution.

6. blast furnace material distribution process charge level exports the modeling method of shape according to claim 5, and it is characterized in that, described step 9 comprises:

Step 9.1: the constraint condition of the material that is periphery by burden distribution matrix structure i-th:

${\∫}_0^{r}2\mathrm{\π}y\left({\mathrm{\γ}}_{l_i}\left(y\right)-{\mathrm{\γ}}_{l_{i-1}}\left(y\right)\right)dy={\mathrm{\κ}}_i{V}_u;$

Step 9.2: row write the piecewise function of i ring:

$y_{l_i}\left(y\right)=\left\{\begin{array}{cc}{\mathrm{\ξ\σ}}_i\exp (-\frac{{(y-g({\mathrm{\α}}_i\left)\right)}^2}{2{\mathrm{\σ}}_i^2}),& {\mathrm{\γ}}_{l_i}\left(y\right)\≥{\mathrm{\γ}}_{{l_i}_{-1}}\left(y\right)\\ {\mathrm{\γ}}_{{l_i}_{-1}}\left(y\right),& {\mathrm{\γ}}_{{l_i}_{-1}}\left(y\right)\≥{\mathrm{\γ}}_{l_i}\left(y\right)\end{array}\right.;$

Step 9.3: the integral constraint according to step 9.1 solves form parameter σ _iiterative solution.

7. blast furnace material distribution process charge level exports the modeling method of shape according to claim 3, and it is characterized in that, the described cloth output function of different rings position is formed with corresponding superposing of integral constraint function by the section unit basis function of drop point site process.

8. blast furnace material distribution process charge level exports the modeling method of shape according to claim 7, it is characterized in that, blast furnace material distribution process charge level exports basis function that shape function exports charge level by the individual pen cloth set and superposes with number of turns weights and form.

9. blast furnace material distribution process charge level exports the modeling method of shape according to claim 7, it is characterized in that, described drop point site process specifically comprises the distance calculating furnace charge drop point distance blast furnace center, and α is the cloth angle of chute:

$v_1=\sqrt{2{\mathrm{gH}}_y+v_0^2;}$

v ₂＝ηv ₁cosα；

$v_3=\sqrt{\begin{array}{c}2{\mathrm{gl}}_0\left(\cos \mathrm{\α}-\mathrm{\μ}\sin \mathrm{\α}\right)\\ +4{\mathrm{\π}}^2{\mathrm{\ω}}^2{l}_0^2\sin \mathrm{\α}{\left(\sin \mathrm{\α}+\mathrm{\μ}\cos \mathrm{\α}\right)}^2\\ +v_2^2\end{array}};$

$x_y=\frac{{\mathrm{mv}}_s^2{\sin }^2\mathrm{\α}}{mg-F_m}\left(\sqrt{\frac{1}{{\tan }^s\mathrm{\α}}+\frac{2mg-F_m}{{\mathrm{mv}}_3^2{\sin }^2\mathrm{\α}}\left(l_0\left(1-\cos \mathrm{\α}\right)+h\right)}-\frac{1}{\tan \mathrm{\α}}\right);$

x＝g(α)＝l ₀sinα+x _y。