CN103399779B - A kind of burning clinker of cement process semi-physical emulation platform and emulation mode thereof - Google Patents

Abstract

The present invention relates to a kind of burning clinker of cement process semi-physical emulation platform software emulation device, chromacoder, control system, supervising device be linked in sequence; Its emulation mode comprises the following steps: input process parameter in the burning clinker of cement process dynamic model that supervising device is set up also is sent to control system; Control system exports procedure parameter to chromacoder, carries out data conversion to export in the internal memory of software emulation device again and to store through chromacoder; Software emulation device obtains virtual duty parameter; The virtual duty parameter that software emulation device exports is fed back to control system through chromacoder; The virtual duty parameter of supervising device receiving control system also carries out graphic software platform by human-computer interaction interface.Modeling of the present invention is set out according to quality, energy conservation and reaction macrokinetics ultimate principle, can describe the dynamic rule of the alternate energy of gas-solid in kiln and mass exchange from the angle of mechanism.

Description

A kind of burning clinker of cement process semi-physical emulation platform and emulation mode thereof

Technical field

The present invention relates to burning clinker of cement technical field, especially relate to dynamic modeling and the semi-physical emulation platform of cement burning process.

Background technology

Burning clinker of cement process have nonlinearity, strong coupling and time stickiness feature.Angle based on continuity of operation and security is considered, new staffs training, design of new processes and risk assessment and the evaluation of various alternatives feasibility etc. are all difficult to realize site test.Even if having ready conditions research experiment, also a large amount of wastings of resources can be caused.In addition, clinker burning process is carried out in a rotating rotary kiln, and internal air temperature is up to thousands of degree, and with very large dust, a lot of key parameter cannot directly be measured, clinker quality is caused to control to rely on operating experience, human factor impact is comparatively large, is theoretically unsound.Therefore, exploitation one can substitute actual clinker burning process, is convenient to experimental study and the method for more complete procedure parameter can be provided to be very necessary.

And system simulation experiments can simulate the various operating conditions of actual calcination process, the danger and extreme operating condition that cannot realize can be tested by model engineering, comparatively complete data can also be provided for control algolithm.Simulation can produce the production result under different technology conditions, dispatcher can utilize these data analysis produced, preferred plan is determined on the basis of flow process, also the basic data of these data as theoretical research can be used, and without the need to destroying the continuity of producing.Simulation Test has good controllability, without destructive, can repeatedly repeat, economic, safety, not by the restriction of meteorological condition and space enrironment.

In emulation technology, the degree of confidence of hardware-in-the-loop simulation is the highest, has the plurality of advantages such as validity, repeatability, economy, security, is subject to the great attention of military and civilian all departments.From the viewpoint of system, hardware-in-the-loop simulation allows to access part material object in systems in which, mean part material object to be put and investigate in systems in which, thus parts can be checked in the environment meeting entire system performance index, while reduction modeling complicacy, improve reliability and the development quality of system.Therefore, be necessary that the semi-physical emulation platform developing a set of cement burning process substitutes true clinker burning process.

Summary of the invention

Commerical test is carried out for being difficult to of existing in clinker burning process, operation lacks the problems such as theoretical direction, the present invention proposes a kind of method based on semi-true object emulation technology, material object is connected with clinker burning model, open so as to describing, the noncontinuity process such as parking, quality and energy implementation, discrete logic operates, and new departure design and control program compare.

The technical scheme that the present invention is adopted for achieving the above object is: a kind of burning clinker of cement process semi-physical emulation platform, and software emulation device, chromacoder, control system, supervising device are linked in sequence;

Software emulation device: for setting up burning clinker of cement process dynamic model and carrying out burning clinker of cement emulation, obtain virtual duty parameter and feed back to chromacoder, and the controling parameters of receiving signal conversion equipment;

Chromacoder: the data for realizing software emulation device and control system are changed and real-time Communication for Power;

Control system: the virtual duty parameter transmitted through chromacoder is fed back to supervising device, and exports the procedure parameter of the supervising device of reception to chromacoder, realize the closed-loop control of burning clinker of cement simulation process;

Supervising device: with human-computer interaction interface, for the controling parameters of the virtual duty parameter and control system of monitoring burning clinker of cement process dynamic model.

Described chromacoder adopts data collecting card.

Described control system adopts PLC.

A kind of burning clinker of cement process Hardware In The Loop Simulation Method, comprises the following steps:

Input process parameter in the burning clinker of cement process dynamic model that supervising device is set up also is sent to control system; Control system exports procedure parameter to chromacoder, carries out data conversion to export in the internal memory of software emulation device again and to store through chromacoder;

Software emulation device, by setting up the procedure parameter in interface module invoke memory, obtains virtual duty parameter by running burning clinker of cement process dynamic model; Again realized the real-time Communication for Power of software emulation device and data collecting card by interface module, the virtual duty parameter that software emulation device exports is fed back to control system through chromacoder; The virtual duty parameter of supervising device receiving control system also carries out graphic software platform by human-computer interaction interface.

Described virtual duty parameter is burning zone temperature.

Described procedure parameter comprises ature of coal, kiln type specification and kiln speed.

Described interface module of setting up refers in simulator, to utilize Visual C++ to compile dynamic link library, and dynamic link library comprises gPROMS digital independent output function and Labcard driver routine call function; Input-output function is used for the read/write operation of gPROMS data, realizes the input and output function of data; Labcard driver routine call function, for calling the driver of board, realizes the function of gPROMS data write board and board data write gPROMS.

Described burning clinker of cement process dynamic model of setting up comprises the following steps:

First according to initial conditions such as ature of coal, kiln type specification and kiln speed, the functional relation T of gas solid matter and kiln temperature is obtained according to chemical kinetic reaction formula and material balance formula _l=f (n _i/j, l), T _lfor the kiln temperature of infinitesimal l, n _i/jfor the molar weight of gas solid matter in kiln, l is infinitesimal in kiln;

Wherein, chemical kinetic reaction formula is $\stackrel{\·}{{\mathrm{\Δn}}_{i/j}}=V_{R,l}{\mathrm{\ΠC}}_{i/j}^{\mathrm{\α}}k{\left(\frac{T}{T_0}\right)}^{\mathrm{\β}}\exp (\[-\left(\frac{E}{R}\right)(\frac{1}{T}-\frac{1}{T_0})\])---\left(1\right)$

represent the molar weight of reaction generation or the solid phase components i consumed or gaseous component j in time per unit infinitesimal; V _r,lfor the reaction volume in infinitesimal l; C ^α _i/jfor the effecting reaction concentration of solid phase components i or gaseous component j; α is the order of reaction; K is pre-exponential factor; T is temperature of reaction; T ₀for reference temperature; β is temperature coefficient; E is reaction activity; R is gas law constant;

Material balance formula is $\frac{\∂n_i}{\∂t}=-\frac{\∂{\mathrm{\ν}}_s{n}_i}{\∂l}+\mathrm{\Δ}{\stackrel{\·}{n}}_i---\left(2\right)$

N _irepresent the molar weight of solid phase components i in infinitesimal; T represents the time; v _srepresent solids movement speed, obtain according to kiln type specification and kiln speed; L represents the axial infinitesimal of kiln; represent the molar weight of the solid phase components i that reaction generates or consumes in time per unit infinitesimal;

Then brought into by this functional relation in solid phase energy equilibrium formula (3), gas solid matter molar weight phase in cancellation solid phase energy equilibrium formula (3), formula (3) only comprises a unknown quantity T _l, equation can be separated, and moving model obtains burning zone temperature side line;

Wherein, solid phase energy equilibrium such as formula (3) is

$\frac{\∂q_s}{\∂t}=-\frac{\∂{\mathrm{\ν}}_s{q}_s}{\∂l}+\mathrm{\Δ}{\stackrel{\·}{q}}_{s,r}+\mathrm{\Δ}{\stackrel{\·}{q}}_{s,t}---\left(3\right)$

Q _srepresent solid phase mole sensible heat in infinitesimal, determined by material thermal capacitance and temperature; represent solid phase reaction heat flux; represent and the heat flux that solid phase receives obtain according to ature of coal.

The present invention has following beneficial effect and advantage:

1. the present invention is based on gPROMS and solve burning clinker of cement process, improve the speed solved.

2. modeling is set out according to quality, energy conservation and reaction macrokinetics ultimate principle, can describe the dynamic rule of the alternate energy of gas-solid in kiln and mass exchange from the angle of mechanism.

3. layering of the present invention and modular design method improve readability, the portability of program effectively.

4. the present invention can be applicable to the overall process of the demonstration of product development, design analysis, test assessment, personnel training.

5. the present invention has good controllability, without destructive, repeatability is strong, economic, safety, not by environmental restraint.

6. the present invention facilitates data acquisition and procession, is convenient to the establishment etc. that system and program were studied, carried out to the Theories and methods carried out in data acquisition and procession.

7. the present invention is by setting up the external interface routine call data acquisition card driver of gPROMS, achieves communicating of software emulation device and data collecting card.

Accompanying drawing explanation

Fig. 1 is rotary kiln distributed model structural drawing of the present invention;

Fig. 2 is semi-matter simulating system hardware structure diagram of the present invention;

Fig. 3 is semi-matter simulating system emulation mode process flow diagram of the present invention.

Embodiment

Below in conjunction with drawings and Examples, the present invention is described in further detail.

A kind of burning clinker of cement process semi-physical emulation platform, comprises the software emulation device, chromacoder, control system and the supervising device that are linked in sequence.

Described control system is real programmable logic controller (PLC) (PLC), is connected with data collecting card; Supervising device comprises industrial computer and connected communication processing card.Software emulation device adopts industrial control computer, is connected with chromacoder.Chromacoder gathers I/O board (data collecting card) for message exchange medium with True Data, and complete process signal type is changed; PLC hardware platform adopts actual PLC device, completes closed-loop control and the logic control of control loop, and digital quantity input/output module wherein also can take into account the Function Extension of platform for the signal gathering the software emulation device transmitted through data collecting card; Communication processing card in supervising device is connected with the industrial computer of PLC, supervising device, adopts Profibus communication protocol to carry out exchanges data.

The software that the present invention relates to mainly comprises the human-computer interaction interface in gPROMS, communication interface program, PLC program and the supervising device in software emulation device.Burning clinker of cement process model in software emulation device builds on gPROMS, based on reaction macrokinetics ultimate principle, the conservation of matter and energy conservation relation formula, set up the dynamic model of burning clinker of cement process, equipment and instrument and meter etc., the software emulation device with burning clinker of cement process model replaces actual sinter leaching system; Communication interface program adopts VC++ to compile dynamic link library (DLL), calls Labcard driver program, realizes collection and conversion that software emulation device exports data; PLC plug-in is write and is write by SIMATIC Manager, completes the collection of data, computing and process control; Supervising device adopts WinCC exploit person machine interactive interface (HMI), realizes the supervisory and control of software emulation device and PLC associated control parameters.

Link through I/O plate (data collecting card) between PLC and software emulation device and connect, and connected with the industrial computer of supervising device by communication processing card, according to the status signal that the industrial computer of supervising device sends, state modulator and logic control are carried out to software emulation device.Supervising device input process parameter transmission outputs to software emulation device to PLC after routine processes, and the output quantity of software emulation device shows with the form of numerical value and animation after PLC routine processes in supervising device simultaneously.The software emulation device of described supervising device, PLC, replacement actual production process forms closed loop test system together.

The place that in Cement clinker buring process, clinkering zone is that cement clinker mineral is formed, fire resistive material stands that temperature is the highest, staggered reaction appears in liquid-solid two-phase in kiln, therefore burning zone temperature is the key factor affecting cement clinker quality output, and how controlling, monitor, comparing burning zone temperature is the problem that manufacture of cement must be considered.The present embodiment is that monitored variable illustrates burning clinker of cement process semi-physical emulation platform build process with burning zone temperature.

Shown in accompanying drawing 1 and accompanying drawing 2, a kind of for burning clinker of cement process semi-physical emulation platform, hardware comprises EVOC industrial control computer EIC-2403, PLC adopts Siemens S7-400 system, I/O board adopts and grinds magnificent PCL-711B data collecting card, and communication processing card adopts Siemens CP5611 network interface card.Software platform is gPROMS, Visual studio, SIMATIC Manager and WinCC.

1) software emulation device

Burning zone temperature is the temperature in rotary kiln in one section of region, but not certain any temperature.Burning zone temperature in nsp kiln can be accurately calculated by setting up macro-kinetics, and under utilizing the dynamic mathematical models set up to analyze different operating mode further, as ature of coal, the impacts of change on burning zone temperature such as kiln type specification and kiln speed, thus lay the first stone for understanding cement rotary kiln burning zone temperature variation characteristic in depth.

Described dynamic mathematical models comprise burning clinker of cement process model and rotary kiln device model, adopt gPROMS to describe the Wuli-Shili-Renli system approach of its correspondence and process and equipment material consumption and energy consumption.Mainly comprise material balance, heat balance and need electric power model.

Because cement rotary kiln is a rotating cylindrical tube, gas-particle two-phase is radially evenly distributed, and can think indifference, and solid phase and gas phase all exist larger difference vertically, so the present invention sets up rotary kiln distributed model vertically, shown in Fig. 1.Can sketch in kiln as raw material are through laser heating, the decomposition (about about 5%) of part thereof, solid phase reaction, liquid-phase sintering reaction and fuel combustion reaction etc., as shown in table 1.The speed responded is such as formula shown in (1).

Main chemical reactions in table 1 cement kiln system

Wherein, k ₀represent reaction pre-exponential factor, E represents energy of activation, C ₂s is dicalcium silicate, C ₃s is tricalcium silicate, C ₃a is tricalcium aluminate, C ₄aF is tetra calcium aluminoferrite.

$\stackrel{\·}{{\mathrm{\Δn}}_{i/j}}=V_{R,l}{\mathrm{\ΠC}}_{i/j}^{\mathrm{\α}}k{\left(\frac{T}{T_0}\right)}^{\mathrm{\β}}\exp (\[-\left(\frac{E}{R}\right)(\frac{1}{T}-\frac{1}{T_0}))---\left(1\right)$

Wherein, represent the molar weight of reaction generation or the solid phase components i consumed or gaseous component j in time per unit infinitesimal; V _r,lfor the reaction volume in infinitesimal l; C ^α _i/jfor the effecting reaction concentration of solid phase components i or gaseous component j, relevant with temperature, partial pressure and solid phase particle diameter; α is the order of reaction; K is pre-exponential factor; T is temperature of reaction; T ₀for reference temperature; β is temperature coefficient; E is reaction activity; R is gas law constant.

Along the Changing Pattern of kiln axis, partial differential mol balance equation is adopted to describe quality transmission and energy transferring in kiln to describe burning zone temperature in kiln.

For solid phase, the material balance of component i is such as formula shown in (2).

$\frac{\∂n_i}{\∂t}=-\frac{\∂{\mathrm{\ν}}_s{n}_i}{\∂l}+\mathrm{\Δ}{\stackrel{\·}{n}}_i---\left(2\right)$

Wherein, n _irepresent the molar weight of solid phase components i in infinitesimal; T represents the time; v _srepresent solids movement speed, relevant with kiln type and kiln speed etc., by angle of repose, kiln rotating speed, kiln body diameter, kiln inclination angle etc. is determined; L represents the axial infinitesimal of kiln; represent the molar weight of the solid phase components i that reaction generates or consumes in time per unit infinitesimal.

Solid phase energy equilibrium is such as formula shown in (3).

$\frac{\∂q_s}{\∂t}=-\frac{\∂{\mathrm{\ν}}_s{q}_s}{\∂l}+\mathrm{\Δ}{\stackrel{\·}{q}}_{s,r}+\mathrm{\Δ}{\stackrel{\·}{q}}_{s,t}---\left(3\right)$

Wherein, q _srepresent solid phase mole sensible heat in infinitesimal, determined by material thermal capacitance and temperature; represent solid phase reaction heat flux; represent the heat flux that solid phase receives, comprise gas phase and conduct heat to solid phase to solid phase heat transfer, kiln wall, provide primarily of coal combustion, and the thermal discharge of coal is by effect on coal quality.

For gas phase, suppose that in kiln, gas is flow at high speed.With solid-phase ratio, the accumulation of gas phase can be ignored, and therefore the material balance of gaseous component j is such as formula shown in (4).

$0=-\frac{\∂{\mathrm{\ν}}_g{n}_j}{\∂l}+\mathrm{\Δ}{\stackrel{\·}{n}}_j---\left(4\right)$

Wherein, v _grepresent gas phase velocity, determined by ideal gas continuity equation; n _jrepresent the molar weight of gaseous component j in infinitesimal; represent the molar weight of the gaseous component j that reaction generates or consumes in time per unit infinitesimal.

Gas phase energy equilibrium is such as formula shown in (5).

$0=-\frac{\∂{\mathrm{\ν}}_g{q}_g}{\∂l}+\mathrm{\Δ}{\stackrel{\·}{q}}_{g,t}+\mathrm{\Δ}{\stackrel{\·}{q}}_{g,r}---\left(5\right)$

Wherein, q _ggas phase mole sensible heat in expression infinitesimal, is determined by material thermal capacitance and temperature or temperature dependence is determined; represent gas-phase reaction heat flux; represent all heat fluxs being delivered to gas phase, comprise coal dust firing heat, gas phase to solid phase heat transfer, kiln wall to gaseous heat transfer.

Cement rotary kiln is a rotating metallic cylinder, maintains the power demand of kiln rotation such as formula shown in (6).

$N=\frac{1}{\mathrm{\η}}(N_1+N_2)=\frac{1}{\mathrm{\η}}(0.55{\mathrm{nLR}}^3{\sin }^3\mathrm{\α}+1.19frGn\frac{D}{d})---\left(6\right)$

Wherein, N represents needs power; N ₁represent useful power; N ₂represent friction horsepower; When η represents mechanical efficiency airtight speed reducer drive, η=0.90; N represents kiln rotating speed; L represents kiln cylinder body length; R represents the effective inside radius of kiln; α represents the half at material center angle in kiln; F represents to rub and examines coefficient, and rolling bearing gets 0.001, and thin oil lubricant sliding bearing gets 0.018, and grease lubrication sliding bearing gets 0.06; γ represents support roller axle journal radius; G represents the quality sum of kiln body rotating part and support roller and axle; D represents wheel belt overall diameter; D support roller diameter.

To sum up, improve boundary condition according to process characteristic and (comprise kiln hood place solid flow and temperature graded rule, gas-phase feed flow, composition and temperature vertically, kiln tail place's gas phase flow rate and temperature be graded rule, solid phase feed rate, composition and temperature etc. vertically) and starting condition (comprising gas phase and solid formation doses, composition and temperature etc. in t=0 moment kiln), Closure equation, sets up the gPROMS dynamic model of cement clinker burning system.

In addition, set up external interface module at software emulation device, function be realize gPROMS data input and output and call I/O Labcard driver program, data directly can output to board, and board data also can be delivered directly to gPROMS.External interface module (i.e. external interface program) adopts dynamic link library form to write, and comprises the data input-output function of gPROMS, Labcard driver routine call function, and the Visual C++ statement of measure for fulfill function; Wherein input-output function is used for the read/write operation of gPROMS data, realizes the input and output function of data; Labcard driver routine call function, for calling the driver of board, realizes the function of gPROMS data write board and board data write gPROMS.

Procedure parameter in software emulation device invoke memory or when carrying out real-time Communication for Power with chromacoder, first operating software simulator, start external interface program, modulating output result exported in the Installed System Memory that external interface program specifies through data input-output function, then Labcard driver routine call function shares parameter in this internal memory, is transferred to Labcard driver program and writes board; Simultaneously, the standard industrial process signals that PLC exports to be converted to after numerical signal in Labcard driver program writing system internal memory through board, then this internal memory of external interface routine call, by in the gPROMS program of parameter read-in software emulation device in internal memory, realize the real-time Communication for Power of software emulation device and chromacoder.

2) chromacoder

The present embodiment adopts and grinds conversion and the output task that magnificent PCL-711B data collecting card realizes software emulation device burning zone temperature signal.Utilize Visual C++ to compile dynamic link library (DLL), set up the external interface program FPI of gPROMS, calling data capturing card driver program realizes ature of coal, the input of semaphore and the outputs of burning zone temperature signal such as kiln type specification and kiln speed.Interface routine comprises two functions: burning zone temperature digital quantity signal and the input of output software simulator are converted into the process variable such as the ature of coal of 0-5V standard industrial process signals, kiln type specification and kiln speed.

3) PLC hardware platform

PLC adopts Siemens S7-400 system, and it mainly completes closed-loop control and the logic control of control loop.According to cement clinker burning system actual condition control flow, utilize SIMATIC Manager software in PLC, write control program and cement clinker burning system is controlled, when software emulation device receives PLC steering order, carry out simulation calculation.

4) supervising device

Supervising device adopts Siemens WinCC exploitation watchdog routine, realize the message exchange between control system, and utilize man-machine interaction (HMI) to realize the supervisory and control of software emulation device input/output argument (burning zone temperature) and PLC associated control parameters (comprising ature of coal, the parameters such as kiln type and kiln speed).Supervising device is connected in series to PLC by CP5611 network interface card, adopts Profibus communication protocol to carry out exchanges data.

Burning zone temperature is the temperature in rotary kiln in one section of region, but not certain any temperature.Burning zone temperature in nsp kiln can be accurately calculated by setting up macro-kinetics, and under utilizing the dynamic mathematical models set up to analyze different operating mode further, as ature of coal, the impacts of change on burning zone temperature such as kiln type specification and kiln speed, thus lay the first stone for understanding cement rotary kiln burning zone temperature variation characteristic in depth.

Semi-matter simulating system general frame as shown in Figure 2.Emulated execution and rate-determining steps are as shown in Figure 3, specific as follows:

Step 1. inputs the process variable such as ature of coal, kiln type specification and kiln speed at man-machine interaction monitoring device and sends, and is sent to PLC through communication processing card;

The process variable such as Step 2. ature of coal, kiln type specification and kiln speed through the control program process of setting, are that standard industrial process signals passes to I/O board by PLC analog quantity/digital output module converts in PLC;

Step 3. standard industrial process signals is that numerical signal is stored into simulation computer internal memory through I/O interface conversion;

Step 4. runs FPI program, and in invoke memory, the numerical signal of process variable such as ature of coal, kiln type specification and kiln speed, runs gPROMS software emulation device, export burning zone temperature side line;

Step 5. calls FPI program again, is that bridge connects software emulation device and I/O board with DLL, is industry standard signal by the numerical signal of the burning zone temperature of modulating output through I/O interface conversion;

The industry standard signal of Step 6. burning zone temperature passes to PLC by analog/digital signal load module, row relax of going forward side by side;

The burning zone temperature signal transmission of Step 7. after PLC process to the human-computer interaction interface of supervising device, visual display;

If Step 8. shows result and expected result is not inconsistent, forward Step 1 to.

Claims (1)

1. a burning clinker of cement process Hardware In The Loop Simulation Method, is characterized in that comprising the following steps:

Input process parameter in the burning clinker of cement process dynamic model that supervising device is set up also is sent to control system; Control system exports procedure parameter to chromacoder, carries out data conversion to export in the internal memory of software emulation device again and to store through chromacoder;

Software emulation device, by setting up the procedure parameter in interface module invoke memory, obtains virtual duty parameter by running burning clinker of cement process dynamic model; Again realized the real-time Communication for Power of software emulation device and data collecting card by interface module, the virtual duty parameter that software emulation device exports is fed back to control system through chromacoder; The virtual duty parameter of supervising device receiving control system also carries out graphic software platform by human-computer interaction interface;

Described interface module of setting up refers in simulator, to utilize Visual C++ to compile dynamic link library, and dynamic link library comprises gPROMS digital independent output function and Labcard driver routine call function; Input-output function is used for the read/write operation of gPROMS data, realizes the input and output function of data; Labcard driver routine call function, for calling the driver of board, realizes the function of gPROMS data write board and board data write gPROMS;

Described burning clinker of cement process dynamic model of setting up comprises the following steps:

First according to initial conditions such as ature of coal, kiln type specification and kiln speed, the functional relation T of gas solid matter and kiln temperature is obtained according to chemical kinetic reaction formula and material balance formula _l=f (n _i/j, l), T _lfor the kiln temperature of infinitesimal l, n _i/jfor the molar weight of gas solid matter in kiln, l is infinitesimal in kiln;

Wherein, chemical kinetic reaction formula is ${\stackrel{\·}{\mathrm{\Δ}n}}_{i/j}=V_{R,l}{\mathrm{\ΠC}}_{i/j}^{\mathrm{\α}}k{\left(\frac{T}{T_0}\right)}^{\mathrm{\β}}\exp (\[-\left(\frac{E}{R}\right)(\frac{1}{T}-\frac{1}{T_0})\])---\left(1\right)$

represent the molar weight of reaction generation or the solid phase components i consumed or gaseous component j in time per unit infinitesimal; V _r,lfor the reaction volume in infinitesimal l; C ^α _i/jfor the effecting reaction concentration of solid phase components i or gaseous component j; α is the order of reaction; K is pre-exponential factor; T is temperature of reaction; T ₀for reference temperature; β is temperature coefficient; E is reaction activity; R is gas law constant;

Material balance formula is $\frac{\∂n_i}{\∂t}=-\frac{\∂{\mathrm{\ν}}_s{n}_i}{\∂l}+\mathrm{\Δ}{\stackrel{\·}{n}}_i---\left(2\right)$

N _irepresent the molar weight of solid phase components i in infinitesimal; T represents the time; v _srepresent solids movement speed, obtain according to kiln type specification and kiln speed; L represents the axial infinitesimal of kiln; represent the molar weight of the solid phase components i that reaction generates or consumes in time per unit infinitesimal;

Then brought into by this functional relation in solid phase energy equilibrium formula (3), gas solid matter molar weight phase in cancellation solid phase energy equilibrium formula (3), formula (3) only comprises a unknown quantity T _l, equation can be separated, and moving model obtains burning zone temperature side line;

Wherein, solid phase energy equilibrium such as formula (3) is

$\frac{\∂q_s}{\∂t}=-\frac{\∂{\mathrm{\ν}}_s{q}_s}{\∂l}+\mathrm{\Δ}{\stackrel{\·}{q}}_{s,r}+\mathrm{\Δ}{\stackrel{\·}{q}}_{s,t}---\left(3\right)$

Q _srepresent solid phase mole sensible heat in infinitesimal, determined by material thermal capacitance and temperature; represent solid phase reaction heat flux; represent and the heat flux that solid phase receives obtain according to ature of coal;

Maintain the power demand of kiln rotation such as formula shown in (6):

$N=\frac{1}{\mathrm{\η}}(N_1+N_2)=\frac{1}{\mathrm{\η}}(0.55{\mathrm{nLR}}^3{\sin }^3\mathrm{\α}+1.19frGn\frac{D}{d})---\left(6\right)$

Wherein, N represents needs power; N ₁represent useful power; N ₂represent friction horsepower; When η represents mechanical efficiency airtight speed reducer drive, η=0.90; N represents kiln rotating speed; L represents kiln cylinder body length; R represents the effective inside radius of kiln; α represents the half at material center angle in kiln; F represents to rub and examines coefficient, and rolling bearing gets 0.001, and thin oil lubricant sliding bearing gets 0.018, and grease lubrication sliding bearing gets 0.06; γ represents support roller axle journal radius; G represents the quality sum of kiln body rotating part and support roller and axle; D represents wheel belt overall diameter; D support roller diameter;

Described virtual duty parameter is burning zone temperature;

Described procedure parameter comprises ature of coal, kiln type specification and kiln speed.