CN110348156B - Method for simulating movement of furnace charge particles in blast furnace rotating chute - Google Patents
Method for simulating movement of furnace charge particles in blast furnace rotating chute Download PDFInfo
- Publication number
- CN110348156B CN110348156B CN201910648452.3A CN201910648452A CN110348156B CN 110348156 B CN110348156 B CN 110348156B CN 201910648452 A CN201910648452 A CN 201910648452A CN 110348156 B CN110348156 B CN 110348156B
- Authority
- CN
- China
- Prior art keywords
- bottom plate
- chute
- charge particles
- particles
- rotating
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- C—CHEMISTRY; METALLURGY
- C21—METALLURGY OF IRON
- C21B—MANUFACTURE OF IRON OR STEEL
- C21B7/00—Blast furnaces
- C21B7/18—Bell-and-hopper arrangements
- C21B7/20—Bell-and-hopper arrangements with appliances for distributing the burden
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Chemical & Material Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Manufacturing & Machinery (AREA)
- Computational Mathematics (AREA)
- Materials Engineering (AREA)
- Metallurgy (AREA)
- Organic Chemistry (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a method for simulating movement of furnace charge particles in a rotary chute of a blast furnace, which comprises the following steps: firstly, determining parameters used for simulation; step two, using ANSYS to obtain tensors of mass and rotational inertia of a rigid body in the rotating chute; step three, performing pretreatment by using ANSYS to generate a k file; step four, calculating the k file generated in the step three, and generating a calculation result after calculation; and step five, submitting the calculation result of the k file to graphical interface software for observing the calculation result. The simulation method can obtain the motion state of the furnace charge particles in the rotary chute at any time, and the result is more visual and concrete, so that a certain guiding function is provided for the research of the problems of the motion rule of the furnace charge particles in the chute, the distribution rule of the furnace charge particles in the blast furnace, the abrasion condition of the chute lining plate when the chute lining plate is impacted by the furnace charge particles and the like.
Description
Technical Field
The invention relates to the technical field of simulation of movement of furnace charge particles in a rotary chute of a blast furnace, in particular to a combined simulation method of movement of the furnace charge particles in the rotary chute of the blast furnace based on ANSYS/LS-DYNA.
Background
The bell-less top technology is one of the most widely used blast furnace iron smelting technology in modern steel industry. The rotary chute is an important device of a charging material distributing device of a bell-less blast furnace, is arranged at the throat part of the blast furnace and is used for charging material filling and distributing. During charging, a large amount of charge particles (including coke, iron ore, etc.) are discharged from a charging bucket positioned at the top of the blast furnace, fall into the chute under the action of gravity, and then flow into the blast furnace under the guiding action of the rotary chute. The rotary chute has two basic movement forms of rotation and tilting, the axis of the rotation movement is the central line of the blast furnace, and the tilting movement can change the inclination angle of the chute. In the charging process, the circumferential distribution of the furnace burden in the blast furnace can be enabled to be uniform through the rotation motion of the chute, and the position of the furnace burden along the radial direction of the section of the blast furnace can be adjusted through the tilting motion, so that the radial distribution condition of the furnace burden in the blast furnace can be adjusted. The rotation and the tilting are cooperated to realize various cloth modes such as annular cloth, spiral cloth and the like.
The dynamics research of the burden inside the rotary chute belongs to an important research subject of the movement of the burden of the blast furnace. The method has important guiding significance for revealing the movement rule of the furnace burden, revealing the distribution condition of the furnace burden in the blast furnace, revealing the abrasion condition of a lining plate on the inner wall of a chute when the furnace burden is impacted, and the like. Currently, there are related mathematical theoretical models for the analysis of the movement of charge particles within the chute. However, these mathematical theoretical models are simplified, for example, only a single charge particle is considered, or the cross-sectional shape of the chute is ignored, and there is a certain difference from the actual working situation, which is not enough to reveal the complete movement state of the charge particle inside the chute. Therefore, it is necessary to solve the above problems by means of a numerical simulation method. At present, no report is available on a numerical simulation method for movement of charge particles in a rotating chute. The invention establishes a Method based on a Finite Element Method (FEM) simulation model to analyze the motion condition of charge particles in a chute.
Disclosure of Invention
The invention aims to provide a method for simulating the movement of charge particles in a rotating chute of a blast furnace, which aims to solve the problems.
In order to achieve the purpose, the invention adopts the following technical scheme:
a method for simulating movement of furnace charge particles in a blast furnace rotating chute comprises the following steps:
firstly, determining parameters used for simulation;
step two, using ANSYS to obtain tensors of mass and rotational inertia of a rigid body in the rotating chute;
step three, performing pretreatment by using ANSYS to generate a k file;
step four, calculating the k file generated in the step three, and generating a calculation result after calculation;
and fifthly, submitting the calculation result of the k file to graphical interface software for observing the calculation result.
Further, the parameters to be determined in the first step include shape and size of the rotary chute lining plate, particle size of the charge particles, material parameters of the rotary chute lining plate, material parameters of the charge particles, static friction coefficient and dynamic friction coefficient between the charge particles and the rotary chute lining plate, and static friction coefficient and dynamic friction coefficient between the charge particles.
Further, the second step specifically includes the following steps:
step (1), appointing a unit type;
step (2), defining the material of the rotating chute bottom plate;
step (3), establishing a geometric model of a rotating chute bottom plate;
adding materials to the geometric model of the rotating chute bottom plate, and dividing the grid of the rotating chute bottom plate;
and (5) calculating and displaying the mass and the moment of inertia tensor of the rotating chute bottom plate.
Further, the second step specifically comprises the following steps:
step (1), specifying a unit type: adopting a solid164 as a unit type of a rotating chute bottom plate, and setting by using an ET statement;
step (2), defining the material of the rotating chute bottom plate: defining a material and designating the material type as a rigid body; taking the density, the elastic modulus and the Poisson's ratio of the rotating chute bottom plate as the parameters of the material to be input;
step (3), establishing a geometric model of the rotating chute bottom plate: establishing a geometric model of the bottom plate of the rotary chute according to the size of the bottom plate of the rotary chute, selecting a Z axis as a rotating shaft of the rotary chute, and setting an included angle between the axial direction of the rotary chute and the Z axis according to the requirement of a tilting angle;
adding materials for the geometric model of the rotating chute bottom plate and dividing the grid of the rotating chute bottom plate: appropriately designating the size of the grid of the rotating chute bottom plate, designating the unit type established in the step (1) as the unit type of the rotating chute bottom plate, and dividing the grid of the rotating chute bottom plate by taking the material type established in the step (2) as the material type of the rotating chute bottom plate;
and (5) calculating and displaying tensors of mass and moment of inertia of the rotating chute bottom plate: VSUM statements are adopted to display the tensors of the mass and the moment of inertia of the rotating chute bottom plate, and relevant data are recorded.
Further, the third step specifically includes the following steps:
step a, appointing a unit type;
step b, defining material parameters;
step c, establishing a geometric model;
step d, dividing grids;
e, copying a large amount of furnace charge particles;
step f, generating a component;
step g, defining a contact type;
step h, defining boundary conditions;
step i, setting simulation duration and simulation time step length;
and j, specifying the file name of the k file and generating the k file.
Further, the third step specifically comprises the following steps:
step a, specifying a unit type: the solid164 is adopted as a unit type of simulation, and ET statements are used for setting;
step b, defining material parameters: specifying three materials which are respectively a material of furnace charge particles, a material of a rotary chute lining plate and a material of a rotary chute bottom plate, wherein the material setting of the rotary chute bottom plate is consistent with the material setting when the mass and the moment of inertia tensor of the rotary chute bottom plate are calculated in the step two;
step c, establishing a geometric model: establishing a geometric model of the rotary chute according to the geometric dimension of the rotary chute and the inclination angle of the rotary chute, wherein the geometric model component of the rotary chute is divided into a rotary chute lining plate and a rotary chute bottom plate; establishing a geometric model of the first charge particles at the origin;
step d, grid division: respectively assigning unit types, material types, the length of a line segment to be divided into grids or the size of the grids for the first furnace charge particles, the rotary chute lining plate and the rotary chute bottom plate, and then dividing the grids;
e, duplicating furnace charge particles: determining the position of a large amount of charge particles to be generated at a certain position above the rotary chute lining plate and away from the discharge opening; using the random number function of the third-party computer language to designate the position points where the burden particles are generated, wherein the distance between every two generated burden position points is greater than the particle size of the burden, so that the generated burden particles are not geometrically overlapped to be more consistent with the actual working condition; copying the established first charge particles on each charge position point, wherein a geometric model, a grid, a unit and a node of the first charge particles are copied during copying; modifying the copied real constant number of each charging material particle to ensure that all the charging material particles have different real constant numbers; deleting the geometric model, the unit, the grid and the node of the first charging material particle after copying all the charging material particles;
step f, a generating component: the generated components comprise a rotary chute lining plate, a rotary chute bottom plate and different charge particles;
step g, defining the contact type: setting contact types among the charging material particles, defining the friction coefficient between every two charging material particles, and calling nested loops to define the mutual action of a large number of charging material particles by adopting a loop statement function in an APDL (advanced persistent program language);
step h, defining boundary conditions: adding gravity acceleration along the Z-axis direction to all the furnace charge particles; adding rotation along the Z-axis direction to the rotating chute bottom plate, so that the whole chute assembly rotates around the Z-axis under the driving of the chute bottom plate; a rotation center is defined, a rotational inertia tensor is defined, the values of the rotational inertia tensor and the mass obtained in the second step are input, and the values are loaded on a rotating chute bottom plate; setting the displacement of the rotating chute bottom plate along the Z-axis direction to be zero; estimating the initial velocity of the charging material particles along the Z axis according to the distance between the discharge opening and the charging material particles, and adding the initial velocity to the charging material particles;
step i, setting simulation duration and simulation time step length: setting simulation duration and simulation TIME step length by using a TIME statement and an EDRST statement respectively;
step j, appointing the file name of the k file, generating the k file: and outputting a k file of the simulation model by using an EDWRITE statement.
Further, the third-party computer language in step e is Python or Matlab, etc.
Further, in the fourth step, the k file generated in the third step is submitted to a solver carried by ANSYS or an LS-DYNA solver for calculation, and a calculation result is generated after calculation;
and further, submitting the calculation result of the k file in the fifth step to LS-Prepost of a graphical interface software LS-DYNA for observing the calculation result. The invention has the following beneficial effects:
the method for simulating the movement of the furnace charge particles in the blast furnace rotating chute can be used for simulating the movement of the furnace charge particles in the rotating chute when the rotating chute performs annular distribution. During annular material distribution, the chute rotates around the center line of the blast furnace at a constant speed, and the tilting angle of the chute is kept unchanged. The software adopted by simulation is ANSYS and LS-DYNA, and the finite element method is adopted for establishing the simulation model. In the simulation model, the quantity of the charging particles and the shapes of the charging particles can be set automatically according to different research contents. The method can obtain the motion state of the furnace charge particles in the rotary chute at any time, and the result is more intuitive and concrete, so that a certain guiding function is provided for the research of the problems of the motion rule of the furnace charge particles in the rotary chute, the distribution rule of the furnace charge particles in the blast furnace, the abrasion condition of the lining plate when the lining plate of the rotary chute is impacted by the furnace charge particles and the like.
Drawings
FIG. 1 is a flow chart of the present invention;
FIG. 2 is a schematic representation of the operation of the rotating chute of the present invention;
FIG. 3 is a cross-sectional view of a rotating chute of the present invention;
FIG. 4 is a representation of the mass and moment of inertia tensors of a rotating chute floor in a simulation model obtained through ANSYS;
FIG. 5 is a diagram of different components and their corresponding parameters such as component number, material number, unit type number, real constant type number, etc. established in the present simulation model;
FIG. 6 is a simulation model of a charge particle having a spherical shape of 25 particles;
FIG. 7 is a simulation result of the movement of 25 spherical charge material particles within the rotating chute at a simulation time;
FIG. 8 is a simulation model with 100 spherical charge particles;
FIG. 9 shows the simulation results of the movement of 100 spherical charge particles within the rotating chute at a certain simulation time;
FIG. 10 is a simulation model of a charge particle having a 1-pellet shape;
FIG. 11 shows the simulation results of the movement of 1 charge particle in the rotating chute at a certain simulation time;
FIG. 12 is a simulation model with 25 square charge particles.
Fig. 13 shows the simulation result of the movement of 25 square charge material particles in the rotating chute at a certain simulation time.
Shown in the figure: 1-discharge opening, 2-tilting rotation point, 3-charging material particles, 4-rotary chute, 401-rotary chute lining plate, 402-rotary chute bottom plate, 5-mass of the rotary chute bottom plate, 6-rotary inertia tensor of the rotary chute bottom plate, 7-rotary chute lining plate parameters, 8-rotary chute bottom plate parameters and 9-charging material particle parameters.
Detailed Description
Example 1
A method for simulating movement of furnace charge particles in a blast furnace rotating chute comprises the following steps:
firstly, determining parameters used for simulation;
step two, using ANSYS to obtain tensors of mass and rotational inertia of the rigid body in the rotating chute;
step three, performing pretreatment by using ANSYS to generate a k file;
step four, calculating the k file generated in the step three, and generating a calculation result after calculation;
and fifthly, submitting the calculation result of the k file to graphical interface software for observing the calculation result.
Example 2
In order to more clearly explain the contents and features of the present invention, the contents of the present invention are described in detail below with reference to the accompanying drawings, and the specific flow is shown in fig. 1. In this example, the number of charge particles 3 is 25, and the shape of the charge particles 3 is spherical.
1. Determining parameters for simulation
The operation mode and the cross-sectional view of the rotating chute 4 are respectively shown in fig. 2 and fig. 3, and the parameters to be determined include the shape and size of the rotating chute 4, the particle size of the charge particles 3, the material parameters of the rotating chute lining plate 401, the material parameters of the rotating chute bottom plate 402, the material parameters of the charge particles 3, the static friction coefficient (TP _ SF) and the dynamic friction coefficient (TP _ DF) between the charge particles 3 and the rotating chute lining plate 401, and the static friction coefficient (PP _ SF) and the dynamic friction coefficient (PP _ DF) between the charge particles 3.
2. Using ANSYS to obtain tensor of mass and moment of inertia of rigid body in rotating chute
When a rotating body is processed in LS-DYNA, only a rotation boundary condition is added to the rigid body, and the tensors of the mass and the moment of inertia of the rigid body need to be input. The rotating chute 4 component comprises a rotating chute lining plate 401 and a rotating chute bottom plate 402, wherein the rotating chute lining plate 401 is attached to the rotating chute bottom plate 402, the rotating chute bottom plate 402 is a rigid body, and the rotation of the rotating chute bottom plate 402 drives the whole rotating chute 4 component to rotate. To add a rotation boundary condition to the rotating chute base plate 402, the rotating chute base plate 402 is set as a rigid body, and the tensor of mass and moment of inertia of the rotating chute base plate 402 is input. Therefore, before the preprocessing operation is formally performed, the tensor of mass and moment of inertia of the rotating chute base plate 402 is obtained through ANSYS calculation. In order to accelerate the acquisition speed, an APDL language built in ANSYS is adopted for acquisition, and the specific operation is as follows:
2.1 setting Unit type
The 8-node entity unit type Solid164 suitable for ANSYS and LS-DYNA joint simulation is adopted, and ET statement setting is adopted.
2.2 Material defining the floor of the rotating chute
EDMP statements are used to define a rigid body of material, and parameters of the material imparted to the rotating chute floor 402 include density, modulus of elasticity, poisson's ratio.
2.3 building the geometric model of the rotating chute bottom plate
A geometric model of the rotating chute bottom plate 402 is established according to the size of the rotating chute bottom plate 402, a Z axis is selected as a rotating shaft of the rotating chute 4, and an included angle between the axial direction of the rotating chute 4 and the Z axis is set according to the requirement of a tilting angle.
2.4 mesh for dividing rotating chute bottom plate
Since the material cannot be added to the geometric model separately in ANSYS, the material needs to be added by a method of dividing the grid of the rotating chute base plate, and therefore, the rotating chute base plate model needs to be completely subjected to the operation of grid division. The grid size of the rotating pan floor is suitably specified, the 2.1 established unit type is specified as the unit type of the rotating pan floor 402, and the 2.2 established material type is used as the material type of the rotating pan floor 402, dividing the grid of the rotating pan floor 402.
2.5 calculating and displaying tensor of mass and moment of inertia of rotating chute bottom plate
The VSUM statement shows the mass and moment of inertia tensor for the rotating chute floor 402. As shown in fig. 4, since the rotation axis of the rotating chute 4 in the simulation model is the Z axis, it is necessary to record the inertia tensor based on the origin, and the inertia tensor includes six values (Ixx, iyy, izz, ixy, iyz, izx) in total, and record the related data.
3. Pretreatment in ANSYS using APDL language
The entire simulation model was pre-processed using ANSYS software. In order to accelerate the whole modeling process, an APDL language built in ANSYS is adopted to carry out pretreatment operation on a simulation model, and the model comprises charge particles 3, a rotary chute lining plate 401 and a rotary chute bottom plate 402. The whole pretreatment process comprises the steps of specifying unit types, defining material parameters, establishing a geometric model, dividing grids, generating parts (Part), defining contact between different parts, defining boundary conditions, setting simulation duration and simulation time step length, and generating a k file.
3.1 setting Unit type
All unit types in the invention adopt 8-node entity unit types Solid164, and ET statements are used for setting.
3.2 definition of Material parameters
Three materials are specified, namely the material of the charge particles 3, the material of the rotary chute liner 401 and the material of the rotary chute floor 402. The materials of the charge particles 3 and the rotary chute liner plate 401 are set by adopting an MP statement, and parameters such as density, elastic modulus, poisson's ratio and the like of the related materials are respectively specified. The material of the rotating chute base plate 402 is set using EDMP, and the rotating chute base plate 402 is set to a Rigid body (Rigid), and the density, elastic modulus, poisson's ratio of the base plate are set. It is to be noted here that since the boundary conditions of the rotation of the object in LS-DYNA can be added only to a rigid body, the material of the rotating chute floor 402 is to be defined as a rigid body.
3.3 building geometric models
Respectively establishing geometric models of the rotary chute lining plate 401, the rotary chute bottom plate 402 and the charge particles 3 according to the sizes of the geometric models of the rotary chute lining plate 401, the rotary chute bottom plate 402 and the charge particles 3; selecting a Z axis as a rotating shaft of the rotating chute 4, wherein an included angle between a rotating chute lining plate 401 and rotating chute bottom plate 402 combined body and the Z axis is an inclination angle of the rotating chute 4; a geometric model of the first charge of particles 3 is established at the origin, the charge of particles 3 being reduced to a spherical body.
3.4 partitioning the grid
The cell type, the material type, and the length of the line to be gridded or the size of the gridding are specified for the first charge material particles 3, the rotary chute liner 401, and the rotary chute floor 402, respectively, and then the gridding is divided.
3.5 replication of the bulk of the charge particles
Determining the position of a large number of charge particles 3 to be generated in the simulation model area according to the diameter of the discharge opening 1; using a third language such as Matlab or Python to randomly generate 25 discrete points, wherein the distance between two adjacent points is required to be larger than the particle size of the burden particles 3, so that the generated burden particles 3 are not geometrically overlapped to be more consistent with the actual working condition; copying the established geometric model of the first charge of material particles and the corresponding grids, cells and nodes thereof at the generated 25 discrete point positions by using a VGEN copy statement in ANSYS;
it is to be noted here that, since the material parameters of all charge particles 3 are the same, different solid types are assigned to different charge particles 3 in order not to form all charge particles 3 as one and the same component in the subsequent component-forming operation.
Therefore, after copying one charging particle 3, real constant numbers are modified for the copied charging particles 3 by using an EMODIF statement, so that the real constant numbers of the charging particles 3 are respectively 1-25; the purpose of varying the actual number is to make each charge particle 3 a part after the part is produced, hereinafter, of these charge particles 3; after all the charge particles 3 are copied, the geometric model, mesh, cell and node of the first charge particle 3 are deleted.
3.6 generating Components
Using EDPART, the CREATE statement produces different parts, as shown in fig. 5, one for each charge particle 3.
3.7 setting contact
Contact between charge particles 3 and the rotating chute liner 401 is set using the built-in DO END cycle statement in the APDL language, the definition of contact being defined using the EDCGEN statement, which is shown in table 1, for example:
TABLE 1 definition of contact between charge particles and chute liners
In table 1, the contact type between the charge material particles 3 and the rotary chute lining plate 401 is ASTS (Automatic-Surface-to-Surface), and TP _ SF and TP _ DF are static and dynamic friction coefficients between the charge material particles 3 and the rotary chute lining plate 401, respectively.
The contact of a plurality of charge particles 3 with each other is set using a nested DO END loop statement, the definition of contact being defined using the EDCGEN statement, which is shown in table 2, for example:
TABLE 2 definition of the contact of a plurality of charge particles with one another
In table 2, the contact type of the plurality of charge material particles 3 is ASTS (Automatic-Surface-to-Surface), and PP _ SF and PP _ DF are static and dynamic friction coefficients of the charge material particles 3.
3.8 setting boundary conditions
(1) The term DIM, ACCZ is used to define the gravitational acceleration, and the term load is used to add the gravitational acceleration to all the charge granules 3, where the term cm is used to define the whole charge granules 3 as a component (component) and then the gravitational acceleration is added to the whole component.
(2) Adding a rotation speed to the rotating chute floor 402 using EDLOAD statements;
(3) Using DIM, CENTER statement to add the CENTER of rotation of the rotating chute floor 402, in this case the origin;
(4) Defining the moment of inertia tensor using DIM, II statements, and adding the moment of inertia tensor onto the rotating chute floor 402 using EDIPART statements;
(5) The displacement of the rotary chute bottom plate 402 in the Z direction is zero by using the D statement, so that the displacement of the rotary chute liner plate 401 in the Z direction remains unchanged after being impacted by the charge particles 3;
(6) And estimating the initial velocity of the charge particles 3 along the Z-axis direction according to the distance between the discharge opening 1 and the generated charge particles 3, and adding the initial velocity to the charge particles 3. The loop statement can here be used to add the initial velocity separately for each piece of charge material particles 3. Of course, if the position of the charge material particles 3 produced is at the position of the discharge opening 1, the initial velocity may not be added.
3.9 setting simulation duration and simulation time step
And respectively setting the simulation TIME length and the simulation TIME step length by using the TIME statement and the EDRST statement.
3.10 generating an output k file:
the k-file of the simulation model is output using the EDWRITE statement.
4. And submitting the k file generated in the step 3.10 to ANSYS or LS-DYNA operation to generate a calculation result, wherein the simulation model of the furnace charge particles 3 is shown in the attached figure 6.
5. The results of the post-treatment were observed using LS-PrePost software, and the results of the simulation of the movement of charge particles 3 within the rotating chute 4 at a certain moment are shown in fig. 7.
Example 3
The number of charge particles 3 required for simulation study can be arbitrarily set as required, the simulation model and the simulation calculation result in which the charge particles 3 are set to 100 particles are shown in fig. 8 and 9, and the simulation model and the simulation calculation result in which the charge particles 3 are set to 1 particle are shown in fig. 10 and 11. The simulation of other quantities of the burden particles 3 by adopting the method belongs to the protection scope of the invention.
Example 4
Other shapes of charge material particles 3 may be provided as required, and fig. 12 and 13 show a simulation model of a rotary chute 4 with 25 square charge material particles 3 and the simulation results thereof, respectively. The simulation of the furnace charge particles 3 with other shapes by adopting the method also belongs to the protection scope of the invention.
The material of the charge particles 3 in the invention is only one, and a plurality of charge particles 3 with different material types can be arranged, and the simulation of the movement of the charge particles 3 with a plurality of materials in the rotating chute 4 by adopting the method provided by the invention also belongs to the protection scope of the patent of the invention.
Although the present invention has been described in detail with reference to the above embodiments, it should be understood by those skilled in the art that: modifications and equivalents may be made thereto without departing from the spirit and scope of the invention and it is intended to cover in the claims the invention as defined in the appended claims.
Claims (6)
1. A method for simulating movement of furnace charge particles in a rotary chute of a blast furnace is characterized by comprising the following steps:
firstly, determining parameters used for simulation;
step two, using ANSYS to obtain tensors of mass and rotational inertia of a rigid body in the rotating chute;
step three, performing pretreatment by using ANSYS to generate a k file;
step four, calculating the k file generated in the step three, and generating a calculation result after calculation;
step five, submitting the calculation result of the k file to graphical interface software for observing the calculation result;
the rigid body in the rotary chute is a rotary chute bottom plate, and the second step specifically comprises the following steps:
step (1), appointing a unit type;
step (2), defining the material of the rotating chute bottom plate;
step (3), establishing a geometric model of a rotating chute bottom plate;
adding materials to the geometric model of the rotating chute bottom plate, and dividing the grid of the rotating chute bottom plate;
step (5), calculating and displaying the mass and the moment of inertia of the rotating chute bottom plate;
the second step comprises the following specific steps:
step (1), specifying unit types: the solid164 is adopted as a unit type of the rotating chute bottom plate, and ET statements are used for setting;
step (2), defining the material of the rotating chute bottom plate: defining a material and designating the material type as a rigid body; taking the density, the elastic modulus and the Poisson's ratio of the rotating chute bottom plate as the parameters of the material to be input;
step (3), establishing a geometric model of the rotating chute bottom plate: establishing a geometric model of the rotating chute bottom plate according to the size of the rotating chute bottom plate, selecting a Z axis as a rotating shaft of the rotating chute, and setting an included angle between the axial direction of the rotating chute and the Z axis according to the requirement of a tilting angle;
adding materials to the geometric model of the rotating chute bottom plate and dividing the grid of the rotating chute bottom plate: appropriately specifying the size of the grid of the chute bottom plate, specifying the unit type established in the step (1) as the unit type of the chute bottom plate, taking the material type established in the step (2) as the material type of the chute bottom plate, and dividing the grid of the chute bottom plate;
and (5) calculating and displaying tensors of mass and moment of inertia of the rotating chute bottom plate: displaying the mass and the moment of inertia tensor of the rotating chute bottom plate by adopting VSUM statements, and making a record of relevant data;
the third step specifically comprises the following steps:
step a, appointing a unit type;
step b, defining material parameters;
step c, establishing a geometric model;
step d, dividing grids;
e, copying a large amount of furnace charge particles;
step f, generating a component;
step g, defining a contact type;
step h, defining boundary conditions;
step i, setting simulation duration and simulation time step length;
and j, specifying the file name of the k file and generating the k file.
2. The method for simulating the movement of the charge particles in the interior of the blast furnace rotating chute according to claim 1, wherein the parameters to be determined in the first step include the shape and size of the rotating chute lining plate, the particle size of the charge particles, the material parameters of the rotating chute lining plate, the material parameters of the rotating chute bottom plate, the material parameters of the charge particles, the static friction coefficient and the dynamic friction coefficient between the charge particles and the rotating chute lining plate, and the static friction coefficient and the dynamic friction coefficient between the charge particles.
3. The method for simulating the movement of the charge material particles in the blast furnace rotating chute as claimed in claim 1, wherein the specific operation steps of the third step are as follows:
step a, specifying unit types: the solid164 is adopted as a unit type of simulation, and ET statements are used for setting;
step b, defining material parameters: specifying three materials which are respectively a material of furnace charge particles, a material of a rotary chute lining plate and a material of a rotary chute bottom plate, wherein the material setting of the rotary chute bottom plate is consistent with the material setting when the mass and the moment of inertia tensor of the rotary chute bottom plate are calculated in the step two;
step c, establishing a geometric model: establishing a geometric model of the rotary chute according to the geometric dimension of the rotary chute and the inclination angle of the rotary chute, wherein the geometric model component of the rotary chute is divided into a rotary chute lining plate and a rotary chute bottom plate; establishing a geometric model of the first charge particles at the origin;
step d, grid division: respectively assigning unit types, material types, the length of a line segment to be divided into grids or the size of the grids for the first furnace charge particles, the rotary chute lining plate and the rotary chute bottom plate, and then dividing the grids;
e, duplicating furnace charge particles: determining the position of a large amount of charge particles to be generated at a certain position above the rotary chute lining plate and away from the discharge opening; using the random number function of the third-party computer language to designate the position points where the burden particles are generated, wherein the distance between every two generated burden position points is greater than the particle size of the burden, so that the generated burden particles are not geometrically overlapped to be more consistent with the actual working condition; copying the established first charge particles on each charge position point, wherein a geometric model, a grid, a unit and a node of the first charge particles are copied during copying; modifying the copied real constant number of each furnace charge particle to ensure that all the furnace charge particles have different real constant numbers; deleting the geometric model, the unit, the grid and the node of the first charging material particle after copying all the charging material particles;
step f, a generating component: the generated components comprise a rotary chute lining plate, a rotary chute bottom plate and different charge particles;
step g, defining the contact type: setting the contact type between the furnace charge particles, defining the friction coefficient between every two furnace charge particles, and calling nested circulation to define the mutual action of a large number of furnace charge particles by adopting a circulation statement function in an APDL (advanced persistent threat language);
step h, defining boundary conditions: adding gravity acceleration along the Z-axis direction to all the furnace charge particles; adding rotation along the Z-axis direction to the rotating chute bottom plate, so that the whole chute assembly rotates around the Z-axis under the driving of the chute bottom plate; a rotation center is defined, a rotational inertia tensor is defined, the values of the rotational inertia tensor and the mass obtained in the second step are input, and the values are loaded on a rotating chute bottom plate; setting the displacement of the rotating chute bottom plate along the Z-axis direction to be zero; estimating the initial velocity of the furnace charge particles along the Z axis according to the distance between the discharge opening and the furnace charge particles, and adding the initial velocity to the furnace charge particles;
step i, setting simulation duration and simulation time step length: setting simulation duration and simulation TIME step length by using a TIME statement and an EDRST statement respectively;
step j, appointing the file name of the k file, generating the k file: and outputting a k file of the simulation model by using an EDWRITE statement.
4. The method according to claim 3, wherein said third computer language in step e is Python or Matlab.
5. The method for simulating the movement of the furnace burden particles in the blast furnace rotating chute as claimed in claim 1 or 2, wherein in the fourth step, the k file generated in the third step is submitted to a solver carried by ANSYS or an LS-DYNA solver for calculation, and a calculation result is generated after the calculation.
6. The method for simulating the movement of the charge particles inside the rotating chute of the blast furnace as recited in claim 1 or 2, wherein the calculation result of the k file is submitted to the graphic interface software LS-Prepost of LS-DYNA in the fifth step for observing the calculation result.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910648452.3A CN110348156B (en) | 2019-07-18 | 2019-07-18 | Method for simulating movement of furnace charge particles in blast furnace rotating chute |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910648452.3A CN110348156B (en) | 2019-07-18 | 2019-07-18 | Method for simulating movement of furnace charge particles in blast furnace rotating chute |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110348156A CN110348156A (en) | 2019-10-18 |
| CN110348156B true CN110348156B (en) | 2022-10-14 |
Family
ID=68178777
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910648452.3A Active CN110348156B (en) | 2019-07-18 | 2019-07-18 | Method for simulating movement of furnace charge particles in blast furnace rotating chute |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110348156B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111665172B (en) * | 2020-07-03 | 2021-07-27 | 河海大学 | Method for analyzing movement state of microscopic structure of sand particles at structure interface |
| CN112226557B (en) * | 2020-09-27 | 2022-08-19 | 北京科技大学 | Blast furnace burden distribution process evaluation method and system, readable storage medium and application thereof |
| CN112733326A (en) * | 2020-12-21 | 2021-04-30 | 赣江新区澳博颗粒科技研究院有限公司 | Numerical simulation method for material distribution of rotary chute at top of blast furnace |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2014094410A1 (en) * | 2012-12-20 | 2014-06-26 | 中国科学院近代物理研究所 | Particle flow simulation system and method |
| CN106777807A (en) * | 2017-01-13 | 2017-05-31 | 北京航空航天大学 | A kind of random size distribution 3D finite element modelings of powder metallurgy and emulation mode |
| CN108256270A (en) * | 2018-02-27 | 2018-07-06 | 厦门力祺环境工程有限公司 | Defeated coal movement system structure design method for shaping is emulated based on material test and DEM-CFD |
-
2019
- 2019-07-18 CN CN201910648452.3A patent/CN110348156B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2014094410A1 (en) * | 2012-12-20 | 2014-06-26 | 中国科学院近代物理研究所 | Particle flow simulation system and method |
| CN106777807A (en) * | 2017-01-13 | 2017-05-31 | 北京航空航天大学 | A kind of random size distribution 3D finite element modelings of powder metallurgy and emulation mode |
| CN108256270A (en) * | 2018-02-27 | 2018-07-06 | 厦门力祺环境工程有限公司 | Defeated coal movement system structure design method for shaping is emulated based on material test and DEM-CFD |
Non-Patent Citations (1)
| Title |
|---|
| 单颗粒磨料水射流切割石材的仿真研究;关砚聪等;《石材》;20150105(第01期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110348156A (en) | 2019-10-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110348156B (en) | Method for simulating movement of furnace charge particles in blast furnace rotating chute | |
| Liu et al. | Discrete element modeling of realistic particle shapes in stone-based mixtures through MATLAB-based imaging process | |
| Coetzee | Calibration of the discrete element method | |
| AU2012235962B2 (en) | Method and apparatus for discrete element modeling involving a bulk material | |
| He et al. | A GPU-based coupled SPH-DEM method for particle-fluid flow with free surfaces | |
| Angelidakis et al. | CLUMP: a code library to generate universal multi-sphere particles | |
| Martin et al. | Simulation of sintering using a Non Smooth Discrete Element Method. Application to the study of rearrangement | |
| Lv et al. | Simulation of the operation of a fertilizer spreader based on an outer groove wheel using a discrete element method | |
| Xie et al. | Discrete element modeling and simulation of non-spherical particles using polyhedrons and super-ellipsoids | |
| CN110990911B (en) | Rock crystal structure rapid modeling method based on particle cluster substitution | |
| US20200110911A1 (en) | Particle simulation device, particle simulation method, and particle simulation program | |
| Kodicherla et al. | Investigations of the effects of particle morphology on granular material behaviors using a multi-sphere approach | |
| Ding et al. | Numerical investigation of dynamic characteristics of debris bed formation based on CFD-DEM method | |
| Xiao et al. | Parameter calibration and mixing uniformity of irregular gravel materials in a rotating drum | |
| Xiong et al. | Development of an unresolved CFD-DEM method for interaction simulations between large particles and fluids | |
| Gong et al. | A review on the simulation of aggregate morphologies in mixture performances based on discrete element method | |
| Zhang et al. | Fast and efficient particle packing algorithms based on triangular mesh | |
| CN111079266A (en) | Periodic numerical model generation method and device for heterogeneous rock-soil material | |
| CN105950806A (en) | Method and device for determining low-dust area in furnace top of blast furnace | |
| Jonsén et al. | Validation of tumbling mill charge-induced torque as predicted by simulations | |
| CN111080790B (en) | Concave body generation method and device based on convex hull algorithm | |
| Sun et al. | A resolved SPH-DEM coupling method for analysing the interaction of polyhedral granular materials with fluid | |
| Chen et al. | A measuring method of gravel’s coefficient of restitution and discussion | |
| Baumgartner | Prediction of burden distribution and electrical resistance in submerged arc furnaces using DEM modelling | |
| Wu et al. | Concrete fragmentation modeling using coupled finite element-meshfree formulations |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |