Information Multivariable Pid Decoupling Control Method of Electroslag Remelting Process Based on Improved Particle Swarm Optimization (pso) Algorithm

A mathematical model of electroslag remelting (ESR) process is established based on its technical features and dynamic characteristics. A new multivariable self-tuning proportional-integral-derivative (PID) controller tuned optimally by an improved particle swarm optimization (IPSO) algorithm is proposed to control the two-input/two-output (TITO) ESR process. An adaptive chaotic migration mutation operator is used to tackle the particles trapped in the clustering field in order to enhance the diversity of the particles in the population, prevent premature convergence and improve the search efficiency of PSO algorithm. The simulation results show the feasibility and effectiveness of the proposed control method. The new method can overcome dynamic working conditions and coupling features of the system in a wide range, and it has strong robustness and adaptability.


Introduction
The electroslag remelting (ESR) process is an advanced smelting method to make purified steels based on rudiment steel in order to reduce impurity and get the high-quality steel which is defined by: uniformity, density and crystal in vertical [1].The main purpose of the ESR process is to purify metal OPEN ACCESS and get the ingot with uniform density crystallization.The steel undergoing the ESR process has many advantages, such as high-purity, lower sulfur and inclusion of nonmetal, smooth surface of the ingot, the uniform density crystallization and the uniform metal structure and chemical composition [2,3].
From the control point of view, the ESR process is a typical complex controlled object, which has multi-variable, distributed parameters, nonlinear and strong coupling features.During the early stage of research, the normal control methods include voltage swings, constant current, constant voltage and descending power.The voltage swing control method is adopted to adjust the voltage swing amplitude in order to control the electrode movements and maintain the stability of slag resistance [4].The reference [5] mainly discusses the key technology of melting speed control.The mathematical model of the time-variant system does not analyse the fundamental reasons causing voltage fluctuation at the expense of the current control accuracy in exchange for the smoothly control of the voltage swing.With the development of the control theory and modern optimization algorithm, the intelligent control strategies, which do not strictly depend on a mathematical model, are used in the ESR process.The reference [6] analyzes the ESR model, adopts the cooperative control method for the melting rate and position of remelting electrodes, whose parameters are optimized by the improved genetic algorithm.The reference [7] proposes an intelligent control method for the ESR process, which combines a variable frequency drive control strategy with the artificial neural network theory.The experimental results show that the proposed control strategy improves the accuracy and intelligence level of the ESR furnace, stables the remelting current.
The proportional-integral-derivative (PID) controller is widely used in many industrial processes because it is simple in structure, and has good robustness and high reliability.The self-tuning and optimization of the PID controller parameters is a very important research field [8,9].With the rapid development of swarm intelligent theory, many PID controller parameters self-tuning methods based on the swarm intelligent theory appeared, such as genetic algorithm (GA), particle swarm algorithm (PSO), etc. [10][11][12][13][14]. Genetic algorithm adopts the duplication, crossover and mutation operators, whose performance has a larger dependence on parameters.On the other hand, PSO will easily lead to the local optima and produce early-maturation phenomena.According to dynamic characteristics of the ESR process, this paper establishes the mathematical model of the ESR process, and puts forward a kind of multivariable self-setting PID control strategy based on improved particle swarm optimization algorithm.The simulation results indicate that the proposed control method is feasible and effective.
The paper is organized as follows.In Section 2, the technique flowchart and mathematics model of the ESR process are introduced.In Section 3, the self-tuning of multivariable PID controller parameters based on IPSO is summarized.In Section 4, experiment and simulation results are introduced in details.Finally, the conclusion constitutes the last part.

Technique Flowchart
Electroslag furnace is a complex controlled object, whose production technique is more complex than other steelmaking methods [15,16].There is still no good mathematical model to describe the ESR process.In order to ensure the stability of the ESR process and quality of the remelting metal, the reasonable technique parameters must be chosen.The technique flowchart of electroslag furnace is described in the Figure 1.The operation flow of the ESR process is illustrated as follows: Add slag→Start arc→Melting slag→Casting→Feeding→Remelting end→Thermal insulation→Finished product.The technique of electroslag furnace may be divided into four stages: melting slag stage, exchange electrodes stage, casting stage and feeding stage.
(1) Melting slag stage.It includes an opening stage to make slag and slag melting with low power by using the lower voltage and current.Because the slag pool has not been formed, there will be large current fluctuations.When the slag pool is formed, the current becomes stable.The slag will be melted until all slag materials have been added.From now on, slag melting will be put under a high power to melt all slag and keep a certain temperature.(2) Exchange electrodes.They will be finished by the movement of the robot arm, which should be finished in no more than five minutes to avoid the solidification of the melted electroslag.
After the completion of the electrode exchange, because there will large heat absorption after the cold metal electrodes are inserted into the electroslag pool, the normal casting should not begin.To achieve the melting speed required by the steel casting technique, the current set-point must be increased.(3) Casting stage.The electrode melting velocity in the normal casting stage is the key factor of deciding the solidification of steel ingot.So the biggest allowed electrode melting velocity of the electroslag furnace may be concluded in accord with the related formulas.(4) Feeding stage.The main function of the feeding stage is to remove the pit on the top of the ingot caused by the solidification.

Mathematic Model
The main circuit of the electroslag furnace consists of a power supply whose voltage and current can be adjusted, a consumable electrode, electroslag, crystallizer, ingot casting and bilge chest.The major loop and its equivalent circuit are shown in the Figure 2. The relationship between the voltage and current of electroslag furnace's major loop is described as follows [1]: where V 2 (t) is secondary voltage, I 2 (t) is smelting current, R s is slag resistance, R V + jX l is short impedance, r is filling ratio, v(t) is melting velocity and K 1 is the resistance coefficient of steel grades.The melting velocity of the electroslag furnace is mainly decided by the voltage and current.Its relationship is described as follows: where K 2 is the coefficient of melting velocity.Without taking into account the influences of the environment temperature and the cooling water, if the melting velocity is kept constant under the constant slag resistance, the control given model for melting current are defined as: where K 3 and K 4 are the initial current coefficient and current equilibrium coefficient.In order to keep outlet voltage stable, the secondary voltage should be calculated by the following formula: where K 5 is initial voltage coefficient, K 6 is current equilibrium coefficient and K 7 is the changeable ratio coefficient of the consumable electrode and ingot resistance.
On the premise of the controlled object under steady working conditions, the system transfer functional is decided by the step response.By imposing the step disturbance to the secondary voltage V 2 and the melting current I 2 respectively, the step response curve of the excitation voltage ΔV and the melting velocity v are calculated.On the basis of considering the variable dimension, the transfer function matrix of the electroslag furnace under the given working conditions obtained by the curve fitting method is illustrated as follows: (5)

Multivariable PID Control Strategy
The PID controller is a regulator in accordance with the linear combination of the proportion, differential and integral of the error, which can be described as: where , and e(t) is the feedback error.
According to the dynamic characteristics of ESR process, the paper puts forward a self-tuning multivariable PID control strategy optimized by the IPSO, whose structure is shown in Figure 3.The system is formed by the PID controller, the diagonal matrix decoupling module and the controlled object.The structure of ESR process multi-variable PID controller which this paper proposed is shown in Figure 3.It uses diagonal matrix decoupling control to eliminate the coupling of the original coupling system between each channel, and improve particle swarm algorithm to optimize control parameters of PID controller.

Diagonal Matrix Decoupling
According the above discussion, the system model can be represented as: In order to realize the decoupling control of ESR process, the diagonal matrix decoupling method is adopted, which is described as follows: Thus, Equations ( 10) and ( 11) are obtained based on the transfer functions of the controlled object and the diagonal matrix decoupling strategy.
Finally, the transfer functions after decoupling disposal of ESR process are calculated as follows: 11 2 2.7 5.4 ( ) 3 15 18 8 7 ( ) 15 18 So the system transfer function matrix after decoupling can be expressed as: After decoupling, the original coupling system is decomposed into two independent, insusceptible single-input/single-output (SISO) control channels (System 1 and System 2) shown in Equations ( 15) and ( 16).
The hybrid PSO algorithm is utilized to optimize the parameters of two PID controllers in order to obtain the better control performances.The design of the PID controller is actually a multidimensional function optimization problem in essence, so the improved PSO algorithm adopts real value coding pattern.The multi-variable PID controller of ESR process can be directly coded as follows.
, , , , , The optimization object of the PID Controller parameters is to make the system overall control deviation tend to zero, and have a fast response speed and smaller overshoot.So the four error integral criterions described as follows are used to evaluate the control effect. (

Particle Swarm Optimization Algorithm
Particle swarm optimization algorithm [17][18][19] is a kind of swarm intelligence evolution computing technology originating from the flock foraging behavior.Usually, a particle swarm is composed of M particles and the position of these particles X i (i = 1,2,•••,M) stands for the potential solutions in the D dimension.The possibility of the potential better solution is larger in the neighborhood space of good particles.These particles fly with a speed V i according to the flight experiences of the individuals and groups, which is adjusted dynamically in D search space.The best position which each particle experiences indicates pbest [i], and the best particle position in the whole swarm indicates gbest [g].The velocity and position of the particles are adjusted according to the following formula.
where c 1 and c 2 are given constants called learning factors, and rand() is a random number within range [0,1].
In the course of the evolution of the flight speed V i , the inertia weight w is introduced to increase the particle's motion inertia, expand the search space and improve the convergence of the discussed PSO algorithm.The PSO algorithm with inertia weights is described as follows:

Adaptive Chaotic Migration Mutation Operator
The genetic distance in GA is introduced to PSO algorithm, in which the particles diversity in the particle swarm is determined by the particle genetic distance.On the basis of this, a special mutation operation (adaptive chaotic migration mutation operator) is utilized to process the particles sinking into gathered area [20,21].
When the position of the particle i is X id (t) in the t generation and the global optimal value of particle swarm is P gd (t), the particle genetic distance relative to the P gd (t) is defined as follows: where ( ) i t  is the genetic distance of particle i relative to the overall swarm, ( ) t  is the population mean genetic distance and N pop is population size.All particles, whose genetic distance is less than the average genetic distance, consist of a new particle group S opt , in which the average distance of particles in new swarm is defined as ( ) where |S opt | is the particle number in particle swarm S opt .The smaller ( ) opt t


, the higher the degree of particle collection.Set the valve value of collection degree   , when ( )  , the particle swarm |S opt | is carried out the mutation operation to enhance the diversity of the particle swarm.Before the mutation operation, the particle swarm S opt need to be carried out the chaos optimization.According to the Logistic mapping principle, a chaos optimization system is realized based on the following equation: where z k is the number sequence within the scope of (0,1) and μ is the number sequence within the scope of [3.571448,4].Suppose the disturbance quantity of a chaotic system is Δx = (Δx 1 , Δx 2 , … , Δx n ) and the disturbance range is [−β,β].Each component of z k is carried to the disturbance range: The former |S opt | with minimum fitness replaces the original particle swarm S opt .Then the migration mutation operation is carried out based on the following rules. (3) (1) ( 1) where (1) ( ) x t are different random individuals in the particle swarm S opt ,  is the migration factor, F is the scale factor, t is the current iteration and T is the maximum iteration number.

Algorithmic Procedure
The flowchart that the proposed hybrid PSO algorithm optimizes the PID controller parameters is shown in Figure 4.The learning steps of the proposed algorithm are states as follows [22].
(1) Initialization population.Initialize the population size, random positions, particle velocities, inertia weights, learning factor, collection valve value, chaos perturbation range [−β,β], scale factor F, maximum iteration number, the particle's position range [−x max ,x max ] and the particle flight speed range [−v max ,v max ]. (2) Calculate the fitness of each particle, and determine the optimal value P id (t) of each particle and the optimal value P gd (t) of the whole particle swarm.
(3) According to Equations ( 19) and (20), the particle's position and velocity are updated.Then update the individual optimal value P id (t) of the particle and the optimal values P gd (t) of particle swarm.


   , judge whether the optimal termination condition is satisfied or not.If it is met, the optimization process is finished and output the optimal results.Otherwise, return to step 2. (5) According to Equations ( 24) and (25), the chaotic optimization is processed.Extract the former |S opt | particle with minimum fitness to take the place of the individuals in the original particle swarm S opt .According to Equations ( 26) and ( 27), the chaotic migration mutation operation is carried out.After the mutation operation, return to step 2.

Input-Output Characteristics of ESR Process
Firstly, the simulation experiments are carried out to observe the control results before decoupling, where two input signals X 1 and X 2 are chosen.The coupling response curves with the input signal X 1 = 1 and X 2 = 0 are shown in Figure 5.The coupling response curves with the input signal X 2 = 1 and X 1 = 0 are shown in Figure 6.

PID Controller Optimized by IPSO Algorithm
This paper utilizes the PID diagonal matrix decoupling controller to control the ESR process.The original coupling control system is divided into two independent control loops.The parameters of the PID controller are optimized by the proposed hybrid particle swarm optimization algorithm on two independent control channels.The paper contrasts the simulation control performances of Z-N method and four different kinds of fitness functions (ISE, IAE, ITAE and ITSE).The simulation results are shown in Figure 7 by adopting the proposed decoupling control strategy on the ESR process.

Conclusions
Based on the technological characteristics, a PID decoupling control strategy based on the improved particle swarm optimization algorithm is proposed.The simulation results indicated that the proposed control strategy has many characteristics of good dynamic and steady performance, strong robustness and the adaptability of the various working conditions and anti-interference capacity.

Figure 1 .
Figure 1.Technique flowchart of the ESR furnace.

Figure 2 .
Figure 2. Main circuit and equivalent circuit of ESR process.

Figure 3 .
Figure 3. Configuration of self-tuning PID decoupling control system.

Figure 7 .
Figure 7. Output responses corresponding to different PID tuning methods.(a) Response curves of output 1; (b) Response curves of output 2.

Table 1 .
The statistical datum of the system performance (overshoot, rising time and adjusting time) are shown in Table2.By carrying out the Amplitude comparison of simulation results and performance index with the above discussed the equivalent decoupling controlled models in various fitness functions, it can be seen that the control effect under ITAE fitness function is optimal among other fitness functions, which is suitable to realize the decoupling control of the ESR process.

Table 1 .
Tuned parameters of PID controllers.

Table 2 .
Performance indices of PID controllers.