Consensus Control Design for 360 MN Extrusion Machine Producing Process

This paper mainly addresses the issue of 360 MN extrusion machine and focuses on the stabilization control of main table attitude. We will first introduce the problem and then model the extrusion machine. As the machine is a multi-input multioutput (MIMO) and strong coupling system, it is challenging to apply existing control theory to design a controller to stabilize the main table attitude. Motivated by recent research in the field of multiagent systems, we design a consensus control protocol for our system and derive our convergence conditions based directly on Routh stability criterion. The advantages of the design are also demonstrated by numerical simulation.


Introduction
Thick-walled seamless tubes are extremely important to the development of some heavy industries.With higher performance of heat-resisting and pressure-resisting, the thickwalled seamless tubes are widely used to transfer liquid in the fields of ships, chemicals, power plants, petroleum, and other industries.36,000-ton (360 MN) extrusion machine is the world's biggest ferrous metal vertical extrusion press machine, as shown in Figure 1, which has been put into production in Baotou, China, in 2009.In Figure 1, the left side figure shows process of squeezing tube billet by driving activity crossbeam upward, the upper right figure is high temperature tube billet formed by extrusion, and the low right figure is finished product.The 360 MN extrusion machine can make large-diameter thick-walled seamless tubes for over 600 MW ultrasupercritical and nuclear power plants.Extrusion is a popular mode for steel production.Through the mold fixed on pressure machine, extrusion table pressures the metal blank into the mold cavity to deform the plastic and to obtain the required size, appearance shape, and certain mechanical properties.During the extrusion, the movement and modeling of the metal blank are controlled by the mold.
The 360 MN extrusion machine uses one-step forming technology to make thick-walled tubes.The materials include high strength alloy steel, stainless steel, refractory alloys, and special metals like uranium, titanium, and niobium.The 360 MN extrusion machine improves the production efficiency greatly.A single extrusion process takes 40∼60 seconds.The 360 MN extrusion machine is typical hydromechatronics integration equipment.The main table is pushed upward by six main cylinders underneath, which extrude out the tube.To generate qualified products, we have to design a control strategy to generate the main table's attitude in horizontal direction.
Researches on the relevant fields of extrusion attract more and more attention in recent years.Tibbetts and Wen put forward that the complex material extrusion process can be viewed as three stages including Material Response, Thermomechanical Response, and Press Response in [1].Moreover, they developed a model of the extrusion process specifically for real-time control and figured out the optimization of extrusion control.Researchers have analyzed means to measure the internal temperature in the cross section of the polymer extrusion which is critical to stabilization process [2].However, due to high-power, long-loop, complex extrusion process characteristics, some new problems come up to 360 MN extrusion machine compared with the ordinary ones.
In this paper, we try to solve the control design problem through a simplified model in which the six main cylinders are regarded as six strong coupling agents.We then introduce consensus control method to keep the balance of main table.Development of consensus control can be traced back to 1960s [3].As time goes on, researchers began to use graph Laplacians and their relevant properties to analyze converge of consensus and agents' position or velocity.Based on the theory of Graph Laplacians, Olfati-Saber et al. provide a theoretical framework to analyze consensus algorithms for multiagent networked systems [4].Recently, researchers pay close attention to branches of consensus control, such as behavior of flocks and swarms [5], nonlinear mathematical models [6], optimization of consensus algorithm [7], the effect of the fuzzy sliding-mode [8], and the communication time-delay effect [9].Moreover, the situation where every agent is a multidimensional dynamic system has been considered by Wen and Chen [10].On the other hand, researchers keep a closer eye on experimental studies.Ren and Chao and Bourgeous proposed a mobile actuator and sensor network platform, which can validate the correctness of Consensus Algorithms for Multivehicle Cooperative Control under different conditions [11].We also utilized a mechanismhydraulic cosimulation platform constituted by Adams and Easy5 software [12], which establishes a simulation platform for the design and evaluation of the control system of a 360 MN extrusion machine.The simulation results of the cosimulation platform capture the device characteristics of mechanical and hydraulic systems.
This paper is organized into seven sections.Section 2 describes the system of the huge extrusion machine, builds a mathematic model, and states the problem when applying the existing control strategy.Section 3 designs a new control strategy based on the consensus control.In Section 4, we establish the convergence condition of the new consensus control strategy.Then, in Section 5, several simulations have been conducted on the mathematic platform to examine the result of the theoretical derivation.In Section 6, more simulations are provided on a proved cosimulation platform which is believed to be closer to the reality.Finally, the conclusion is summarized in Section 7. The hydraulic control system of 360 MN extrusion machine is shown in Figure 3. Four SSI digital position transducers placed at the four corners collect the main table position data and transfer the signal to the central controller.Meanwhile, six pressure transducers in each main cylinder collect the extrusion force signals and transfer it to the controller as well.

Problem Statement
The control system output is the set point value for the servo valve PID controller AN430, which, together with the servo valve, piston, and variable pump, comprises a typical electrohydraulic position control system.
The major disturbance in this control system is the uncertain resistance force during extrusion process, due to temperature asymmetry of the ingot, metal turbulent, inconformity among cylinders, and many others.These disturbances will result in the unbalance of the main table, which will lead to the decline of the products.The control purpose is to ensure the horizontal balance of the main table which is actuated upward by a set of main hydraulic cylinders symmetrically installed around the circular plane.In practice, open loop control is used.The disadvantage of open loop control is obvious: performance will become poor when the system faces unexpected nonlinear resistance.

Mathematical Model of the Machine.
In order to study the control of the 360 MN machine, we introduce a simplified mathematical model of the system.The actual dynamics of the model is much more complicated.We built a more realistic simulation model in Adams and Easy5 to test the control design.The details of the simulation model can be found [12].
Abstract the main table of the extrusion machine into a circular plane.Assume that  points are uniformly distributed around the circular plane.At each point , there is an upward force driven upward.The center of the circle is fixed on the horizontal plane.Define the axis and axis on the circular plane as below.The number 0 point locates on the negative direction of axis .Because the  points are uniformly distributed, assume that the angular interval is  = 2/, so the angular between the vector, composed of the center of circle to the  point, and the negative direction of axis  is  (Figure 4).Suppose the radius of the circle is , the angular deviation of axis  is   , and the angular speed of axis  is   .The angular deviation of axis  is   , and the angular speed of axis  is   .Let the total gravity of the main table be  =  and assume it is equally shared by the  cylinders (Figure 5).Then the torque of the force   and the shared gravity / by cylinder  projected at the axis  is So we can get the torque of all force projects at the axis  as From the angular momentum, we can get where   = (1/2) 2 ; set  = 1, so Similarly, the torque of the force   project at the axis  is The total torque for axis  is From the angular momentum, we can get where   = (1/2) 2 ; set  = 1, so We also have So, the dynamics of the plane under input force at the cylinders can be summarized as which is actually the state equation of the system.In matrix form, we have where When there are random frictions caused by the load, we have where In this extrusion machine case, denote   as the position of the number  cylinder's force point, and denote V  as the velocity of the number  cylinder's force point.According to the conclusion in Appendix A, the observation matrix of the position is where 5 Average

Control Strategy Design
The control design problem is to design controller with low complexity to stabilize the main table described by ( 13) with available measurements given in (14).

Consensus Control Strategy Design.
With the idea borrowed from the multiagent and consensus control which was briefly reviewed in Appendix B, we now try to design simple yet effective control strategies for our extrusion machine.
The forty-two pumps are divided into six groups, each one of which is assigned to one cylinder, so the cylinders are independent from each other.Then, we can regard the six cylinders as six agents, and they constitute a multiagent system, the topology of which is the physical constraint of the main table.Our purpose of control is maintaining the same velocity and position of the six cylinders during the extrusion process.
In fact, there are six hydraulic cylinders offering the force, and the topology of the agent network is shown in Figure 6.Our control protocol is where where V 0 is the target speed of the cylinders and  is the average position of the six cylinders.Substituting into system equation ( 13) we can get Remark 1.Although our control protocol looks similar to the existing ones, it should be pointed out that the dynamics of the overall system is more challenging than the systems in the consensus control literature.There, the state equations of the agents are independent.In contrast, in our system, the state equations of the agents are coupled.They have to satisfy the geometry constraint enforced by the main table.(See the geometric derivation in the appendix).
Remark 2. When replacing the  with the  0 , which means that the position is set in advance, the topology connection between the agents has been cut down, and the consensus control turns into independent loop PID (actually PI) control strategy.When we analyze the attitude of the plane model, because the form of the two control strategy is the same, the convergence condition and the numerical simulation results are the same too.But in the actual producing process, because the  is more flexible and more close to reality than  0 , the consensus control's effect is also better than the independent loop PID control's effect.(We will prove it in Section 6 in the means of cosimulation platform).

The Convergence Condition Derivation
In this section, we try to find the convergence condition of our control protocol in Section 3. Depending on the number of cylinders, we will discuss three cases,  = 4,  = 4+2, and  = 2 + 1, for any integer .We will derive our convergence conditions based directly on Routh stability criterion.
In general, due to the symmetry of the distribution of forces, the / part of the   is counteracting to the gravity each hydraulic cylinder subjected to.So we can assume a gravity-free environment, which means   =   , to deduce the convergence condition derivation.
Also, the random frictions are not taken into consideration during our derivation.
Under these assumptions, we can get our convergence condition of the consensus control strategy as follows: Next, we are going to prove this conclusion.Assume point 2 locates at the axis , the relative position between point 2 and the center point is   , and the velocity of the point is V  , and assume point  locates at the axis , the relative position between the point and the center point is   , and the velocity of the point is V  .Without loss of generality, we can assume that V  ,   and V  are all greater than 0 (Figure 7).
From the analysis in Appendix A, we know how to calculate the perpendicular distance between any edge point of the circular plane and the horizontal plane.
Assume the point  is the highest point of all the edge points of the circular point from the horizontal plane, with ∠ = , and the initial position is where The velocity is the first derivative of the position, so we can get where According to (18), we have Assume So the coefficient matrix  can be rewritten as follows: So With the Routh stability criterion, which provides the necessary and sufficient conditions for the convergence of the polynomial, we can get Now, we are trying to find a set of solutions to the array of the inequalities above.
According to the inequalities array (26), we can get the result These can meet the conclusion (19).When  > 1, each force   can be decomposed into two mutually perpendicular directions: axis  and axis , and the coefficients are cos  or sin .And the total effect is just like the case when  = 1; also, the calculation steps are the same with  = 1.
So the conclusion (19) is the solution to the array of the inequalities (26) (Figure 8).

𝑁 = 4𝑛 + 2:
Uniform Distribution.In this situation, the distribution of the agents is also Centrosymmetric and Axisymmetric.We cannot obtain   ,   , V  , V  (or say   ,   ,   ,   ) directly.Assume the  point's relative deviation from the center point is  1 , the velocity of the point is V 1 , the  + 1 point's relative deviation from the center point is  2 , and the velocity of the point is V 2 .
We can get 4 + 2 (or say 0), which locate at the axis , and the point's relative deviation from the center point is   = ( 1 −  2 )/2 cos(/2), the velocity of the point is ), and we can get the virtual point, which locates at the axis , point's relative deviation from the center point is   = ( 1 +  2 )/2 sin(/2), and the velocity of the point is Then, the next steps are the same as the case  = 4, and we can get the same equations ( 22), (25), and (26) too.
Assume  = 0; then  = 4 + 2 = 2.In this case According to the inequalities array, we can get the result These can meet the conclusion (19).When  > 0, each force   can be decomposed into two mutually perpendicular directions: axis  and axis , and the coefficients are cos  or sin .And the total effect is just like the case when  = 4; also, the calculation steps are the same with  = 4.
So the conclusion ( 19) is the convergence condition of this case.

𝑁 = 2𝑛 + 1:
Uniform Distribution.In this situation, the distribution of the agents is Axisymmetric.We cannot obtain   ,   , V  , V  (or say   ,   ,   ,   ) directly either.
Assume the  point's relative deviation from the center point is  1 and the velocity of the point is V 1 , and assume the  + 1 point's relative deviation from the center point is  2 and the velocity of the point is V 2 .
And we can get 2 + 1 (or say 0), which locate at the axis , point's relative deviation from the center point is   = ( 1 +  2 )/2 sin(/2), and the velocity of the point is  the center point is   = ( 1 −  2 )/2 cos(/2), and the velocity of the point is 9).Then, the next steps are the same as the case of  = 4, and we can get the same equations ( 22), (25), and (26) too.Assuming  = 0, then  = 2 + 1 = 1.In this case According to the inequalities array, we can get the result These can meet the conclusion (19).When  > 0, each force   can be decomposed into two mutually perpendicular directions: axis  and axis , and the coefficients are cos  or sin .And the total effect is just like the case when  = 4; also, the calculation steps are the same with  = 4.
So the conclusion ( 19) is the convergence condition of this case.
Above all, we have proved that the conclusion ( 19) is the convergence condition of our consensus control strategy.

Simple Mathematical Simulation
In this section, we do some simulation to illustrate the control design and convergence analysis.We will focus on the specific examples; we will do the simulation on the case in which  = 4 and  = 6.
5.1.Square Case.In this section, there are four agents that need to be controlled, and each of them has two variables: the position and the velocity.The initial value of each variable is set as follows: At first, all the feedback coefficients are set equal to 1, which means The simulation results of the position and velocity of each agent are shown in Figure 10.We can learn the convergence time is about 10 s.
Then, let us carry out simulations for parameters violating the convergence conditions established in Section 6.We set the coefficients as follows: and check the position and velocity of each agent.We can see from Figure 11 that, as expected, the system could not keep convergence any more.Moreover, it is the axis  1 +  3 > 0,  1 +  3 > 0 that could not converge, while the other axis still meets the convergence requirement.

The Hexagon Case.
We design the simulation as the case of square.The initial value of each variable is set as At first, we set the velocity feedback coefficients  1 =  2 =  3 =  4 =  5 =  6 = 1 and the position feedback coefficients The simulation results of the position and velocity of each agent are shown in Figure 12.We can learn the convergence time is about 20 s.Then, we check the situations where the convergence conditions are not satisfied.To be more specific, we run simulation for the case where We check the position and velocity of each agent.As we can see from Figure 13, the system could not keep convergence any more.

Simulation on the Cosimulation Platform
One of the biggest differences between the mathematic simulation and the reality is that the simulation situation cannot take the unforeseen circumstances into consideration.Factors like the temperature variation, measurement error, and feedback time delay are ignored while quite usual in the practice.
To make the simulation more realistic, we built a cosimulation platform consisting of Easy5 and Adams [12].In this platform, Adams and Easy5 can exchange data to each other, while solving their own problems.The cosimulation principle is shown in Figure 14.Easy5 uses mechanical dynamic parameters as its input, simulates the responses of hydraulic system, and passes the results to Adams; Adams adopts hydraulic power parameters as its input, computes the responses of mechanical system, and delivers the consequences to Easy5.The communication of these two software programs uses Virtual Server-Client way to transfer data at settable step.The characteristics including validity and authenticity of the cosimulation platform have been proved in [12], and we will use this cosimulation platform to compare our control strategy with open loop control used currently in practice.In the following simulation, the period time of data exchanging between Easy5 and Adams is 0.01 s.Table 1 lists the maximal position error and velocity error of each cylinder, and the bold ones of each column represent maximal position errors or velocity errors of all six cylinders.Form Figures 16 and 17 and Table 1, we can easily find that the consensus control strategy has much better performance than the open loop control; in fact the maximum position error is reduced by one order of magnitude.

With Loading.
In this case, we compare the open loop control and the consensus control strategy on the cosimulation platform with load.The load curve is shown in Figure 18.All the other settings, including period time, resistance, and control protocol, are the same as in the previous case.The simulation results are shown in Figure 17.Table 2 lists the maximal position error and velocity error of each cylinder, and the bold ones of each column represent maximal position errors or velocity errors of all six cylinders.Form Figures 19 and 20 and Table 2, we can see that the consensus control strategy has much better performance than the open loop control when the system is loaded; the maximal position error is reduced more than 10 times.

Conclusions
In previous practices, because of the strong coupling problem, the 360 MN extrusion machine usually works under open loop control strategy.To improve its security and performance, we build a simplified control model for   a cosimulation platform.The results suggest that the consensus control method is a more effective way for the industrial process control problem.
Future work will emphasize more the application of the consensus like control strategy in reality; we will try to apply the model to the control system of 360 MN vertical extrusion machine in Baotou.We will also deepen theoretical analysis of the problem to include more realistic situations (e.g., effects of limited communication rate similar to [13]).

A. Attitude Analysis of Circular Plane
As we can see from Figure 21, the plane  is a horizontal plane and line  and line  are the two mutually perpendicular axes of the circular plane .The distances between the horizontal plane and point  and point  are  and ℎ.Draw the extension line of  and meet the horizontal plane at point .Connect the point  and point  then extended to point .Draw the vertical line of point  and point  to line  at the point  and point .
As it is easy to know that the Δ and Δ are congruent and Δ is similar to Δ, Δ and Δ are also similar.

Figure 1 :
Figure 1: The 360 MN vertical extrusion press and the largediameter thick-walled seamless tubes.

2. 1 .
Problem Description.Figure2shows the hydraulic system layout.The main pumps are placed on the B1 floor (−7 meters); the tanks and auxiliary pumps are placed on the B2 floor (−14 meters).Forty-two WEPUKO RX500 Radial Piston Variable pumps drive the extrusion process with maximum working pressure 420 bar.A pump provides 1000 L/min flow at 350 bar or 750 L/min at 420 bar.The movement of the pump's embedded piston causes the oscillation of the cam, making the desired change of radial piston position and, as a result, the desired change of the flow.The piston movement is controlled by a MOOG servo valve.The forty-two pumps are divided into six groups, each one of which is assigned to one 1350 mm diameter main cylinder.

Figure 4 :Figure 5 :
Figure 4: The space structure of the circular plane mathematic model.

Figure 6 :
Figure 6: The graph used for consensus control.

4. 1 .
= 4: Uniform Distribution.The distribution of the agents is both Centrosymmetric and Axisymmetric.

Figure 10 :Figure 11 :
Figure 10: (a) is the position of each agent; (b) is the velocity of each agent.

Figure 12 :Figure 13 :
Figure 12: (a) is the position of each agent; (b) is the velocity of each agent.

Figure 16 :
Figure 16: The position error of each cylinder.

6. 1 .
Without Loading.In this case, we compare the open loop control and the consensus control strategy on the cosimulation platform without load.The resistance forces applied on the four corners are shown in Figure 15.The patterns of the disturbances include constant step functions and Sine functions.Consensus protocol used in this simulation is (17), with target speed of cylinders V 0 = 5 cm/s; all the feedback coefficients are equal in groups  1 =  2 =  3 =  4 =  5 =  6 = 150 and all the feedback coefficients are equal in groups  1 =  2 =  3 =  4 =  5 =  6 = 400.The simulation results are shown in Figure 16.

Figure 17 :
Figure 17: The velocity error of each cylinder.

Figure 18 :
Figure18: The load in the simulation.

Figure 19 :
Figure 19: The position error of each cylinder.

Figure 20 :Figure 21 :
Figure 20: The velocity error of each cylinder.
sin(/2), and we can get the virtual point, which locates at the axis , point's relative deviation from

Table 1 :
Maximal position error and velocity error of each cylinder.

Table 2 :
Maximal position error and velocity error of each cylinder.360 MN extrusion machine and design a consensus control strategy to keep the main table work in a balanced attitude.Necessary and sufficient conditions are established for the convergence of the proposed control strategy based on the Routh stability criterion.The control strategy is also tested and compared with the open loop control strategy on the