Calculating Power Parameters of Rolling Mill Based on Model of Deformation Zone with Four-Roll Passes

: Making “digital twins” for rolling processes and mill equipment should begin with the development of mathematical models of the deformation zone. The deformation zone of two-high ﬂat mill rolling have been studied in detail, relevant models are available in many academic papers. However, the same cannot be said about the most complex deformation zones in stands with multi-roll gauge. Therefore, the task of their reliable mathematical description is of immediate interest. The development of mathematical models is necessary for the design of new wire mills and rolling-drawing units. The combination of rolling in stands with multi-roll gauge and drawing is a promising direction in the production of wire from di ﬃ cult-to-form steels and alloys. Digital models for pressure-based metal treatment are also necessary for calculating the rolling-mill power parameters during the development of new assortments at the operating mills. The models of deformation zones present the basis for developing the multivariable control systems of process conditions of continuous mills. This research is devoted to the study of the deformation zone and the development of a procedure for calculating the power parameters of rolling in a stand with four-roll passes. The solution of these challenges is given using the example of an operating ﬁve-stand wire mill. The authors analysed the known analytical dependencies for calculating the rolling mill force and torque. A mathematical model of the deformation zone and a program for calculating the power parameters have been developed. The paper compares the results obtained from calculations based on analytical dependence and on modelling. A comparison with the experimental parameters obtained at the mill is given. The authors assess the feasibility of using the known formulas and analyse the impact of the front and rear tensions on the power parameters of rolling mill. The problem of developing an automatic tension control system for continuous mills with multi-roll groove is substantiated.


Introduction
Rolling mills for the production of steel wire are classified according to the number of mill-stand rolls, which form the geometry of the deformation zone. Accordingly, these can be categorised by the type of roll groove and the roll transmission path. It is customary to refer to a groove formed by three or more rolls with the axes in the same vertical plane (Figure 1), as a multi-roll. Such a composition provides for a high dimensional precision of the manufactured profile. Therefore, modern small-section section and wire mills for rolling steel and non-ferrous metals are designed using finishing mill groups of stands with multi-roll groove.
The most common grooves are: triangular, formed by three smooth rolls; rectangular, formed by four rolls, Figure 1a; hexagonal, formed by six smooth rolls or three grooved rolls, and so on. These are used in different configurations, building complex systems (for example, triangle-triangle), Figure 1b [1,2].
The conditions for metal deformation in a multi-roll groove have advantages over the conditions for deformation in a groove formed by two rolls. Indeed, two-rolls flat rolling usually results in an increase in the length of the roll (draw-down) and its widening (spreading effect). At the same time, rolling in a multi-roll groove is characterized by simultaneous height/width drafting of the roll and a significant draw-down. The principal stress diagrams for the deformation of both types remain the same, the difference is only in values thereof. In multi-roll groove, the height-transverse drafting is a positive factor that provides for rolling low-ductility alloys, which is crucial. Unlike the rolling in two-roll groove, free forging or rotary swaging, drafting of stock material in multi-roll grooves involves three or four sides, rather than two. Moreover, during stretch rolling, for example, a square-shaped stock material is set for a square groove in the "edgeways" position ( Figure 1a) [3]. Multi-sided compression with a varied-width deformation zone leads to the following features:

•
In the deformation zone, a pattern of all-round compression of the workpiece at high hydrostatic pressure is created, which increases the plasticity of the processed material, also providing for better one-time deformations; • The pattern of the deformed state is changing, which makes it possible to eliminate the crossflow of the material, for better one-time deformations, in order to increase one-time draw-downs; • Intensive all-round compression of the workpiece leads to a higher density, better physical and mechanical properties and structure of the product. Figure 1b shows a diagram of longitudinal rolling of metal in a stand with two three-roll grooves, which form the most closely spaced deformation zones in the driven rolls 1-6 [4]. Longitudinal rolling of the round workpiece 7 is carried out to obtain, for example, the hexagonal metal profile 8.
Multi-roll grooves in comparison with two-roll grooves provide for a better rigidity of the stands, and, consequently, a better dimensional precision of rolled products, lower spreading effect and higher drawing capacity of the stand. In addition, the strain pattern in the deformation zone can be easily changed, bringing it closer to all-round compression. This provides favourable conditions for increasing the partial and total drafting [5]. Therefore, the number of rolling stands in the mill line is almost halved in comparison with two-high rolling and drawing in roller dies [6]. Unlike the rolling in two-roll groove, free forging or rotary swaging, drafting of stock material in multi-roll grooves involves three or four sides, rather than two. Moreover, during stretch rolling, for example, a square-shaped stock material is set for a square groove in the "edgeways" position ( Figure 1a) [3]. Multi-sided compression with a varied-width deformation zone leads to the following features:

•
In the deformation zone, a pattern of all-round compression of the workpiece at high hydrostatic pressure is created, which increases the plasticity of the processed material, also providing for better one-time deformations; • The pattern of the deformed state is changing, which makes it possible to eliminate the crossflow of the material, for better one-time deformations, in order to increase one-time draw-downs; • Intensive all-round compression of the workpiece leads to a higher density, better physical and mechanical properties and structure of the product. Figure 1b shows a diagram of longitudinal rolling of metal in a stand with two three-roll grooves, which form the most closely spaced deformation zones in the driven rolls 1-6 [4]. Longitudinal rolling of the round workpiece 7 is carried out to obtain, for example, the hexagonal metal profile 8.
Multi-roll grooves in comparison with two-roll grooves provide for a better rigidity of the stands, and, consequently, a better dimensional precision of rolled products, lower spreading effect and higher drawing capacity of the stand. In addition, the strain pattern in the deformation zone can be easily changed, bringing it closer to all-round compression. This provides favourable conditions for increasing the partial and total drafting [5]. Therefore, the number of rolling stands in the mill line is almost halved in comparison with two-high rolling and drawing in roller dies [6].
The three-roll milling process in two-groove stands ( Figure 1b) is a vital technology for the production of round bars and hexagonal sections used in construction. It provides better mechanical properties of profiles due to large structural yielding [7]. Publications [4,8,9] are devoted to the study of such modes. They investigate the processes of symmetric and asymmetric rolling. The analysis of metal deformation was performed using the DEFORM-3D software. The key results are the substantiation of the drafting patterns and recommendations for the practical use of stands in processing lines of combined rolling-drawing units.
A significant drawback of stands with three-roll grooves is the low load capacity in terms of rolling force and torque supplied to the rolls. This reduces the deformation value of the workpiece in each stand and limits the options for rolling materials with high resistance to deformation [2]. Therefore, rolling stands with four-roll pass have found the greatest practical application.
The first information about a universal mill where the four-roll pass included two horizontal driven and two vertical non-driven rolls, dates back to 1854 [10]. The scientists from the Chelyabinsk Polytechnic Institute (currently South-Ural State University, Russia) were the first to develop stand designs for universal section mills containing replaceable cassettes with three-and four-roll grooves [11,12]. A special feature of these stands was that the work rolls are rotated from contact with two driven backup rolls. At the end of the 1970s, such universal mills began to operate at the Izhstal enterprise (Izhevsk city). Stands of a similar design have been developed by the Research Institute of Metallurgical Engineering (VNIImetmash).
Over the past decades, there have been significant changes in the design of stands and technological schemes for rolling, as well as in drafting modes. Significant experience has been accumulated in operating mills that were manufactured at various metallurgical plants. Constantly growing market requirements for metal goods and products dictate the need to expand the range of rolled assortment.
One of the successfully operated units is a five-stand wire mill with four-roll grooves built at the Beloretsk Metallurgical Plant [13]. Currently, its owners explore the possibility of producing the wire rod with over 8 mm thickness, which is outside the design parameters of the rolling mill. This predetermines the need to improve the methodology for calculating the rolling force, the shaft torque and the power of the actuating unit. Such studies would make it possible, firstly, to assess the possibilities of expanding the product range, and secondly, to determine the limits imposed by the electric drive on the technological process.
In addition to determining the power parameters, the improvement of their calculation methods will provide an opportunity to simulate the rolling process on promising multi-stand mills, taking into account the interaction of electromechanical systems through the metal [14]. This challenge cannot be solved by using well-known software packages based on finite-element models (DEFORM, ABAQUS, ANSYS, Q-FORM, etc.), since it requires runtime assessment of multivariable systems.
The calculation of electric and force parameters is relevant either for the construction of new mills or at the extension of the rolled product range of the existing ones. The analysis of stand force interaction is important for improving automated electric drives and the systems of automatic technological parameters control.
Modern stand design procedures are based on the accumulated experience of using various groove systems. They solve the problem of obtaining the required precision, the stable continuous rolling process, and the required drawing capacity [15]. This determines the relevance of returning to the formation of mathematical models of deformation zones in stands with multi-roll grooves.

Study Object Properties
The object of research is the aforementioned five-stand wire mill with four-roll passes. The mechanical units of this mill were manufactured by SKET (Magdeburg, Germany). It is designed for rolling process wire of 2.5-3.5 mm diameter using 6.5-8.0 mm diameter initial wire rod. Electric drives and control systems have been developed by experts from the Magnitogorsk Mining and Metallurgical Institute (currently, the Magnitogorsk State Technical University).
The layout of the equipment is shown in Figure 2 [13]. The stands are divided into roughing (with tension roll passes) and finishing (with finishing passes). The main work on total deformation of the metal is carried out in the tension roll passes, while the finishing passes are used to form the geometry and size of the finished profile. Table 1 shows options for four-roll passes rolled on the mill.
The rolled stock is fed from a reel or a coil. There is a welding device for welding the ends of the wire rod during continuous rolling. The Wire Accumulator 5 eliminates shock loads on the decoiler if the turns get overlapped. If necessary, the High-Frequency Inductor 6 can heat the stock material to temperatures of warm rolling. The layout of the equipment is shown in Figure 2 [13]. The stands are divided into roughing (with tension roll passes) and finishing (with finishing passes). The main work on total deformation of the metal is carried out in the tension roll passes, while the finishing passes are used to form the geometry and size of the finished profile. Table 1 shows options for four-roll passes rolled on the mill.
The rolled stock is fed from a reel or a coil. There is a welding device for welding the ends of the wire rod during continuous rolling. The Wire Accumulator 5 eliminates shock loads on the decoiler if the turns get overlapped. If necessary, the High-Frequency Inductor 6 can heat the stock material to temperatures of warm rolling.

Description Diagram Draw-Down Ratio
Circle-square µ = 1.57 Circle-octagon µ = 1.11 Octagon-square µ = 1.66 Octagon-octagon µ = 1.17 Octagon-circle µ = 1.05 Constant change of angles during the transition from groove to groove, and extra heating due to intense deformations, create a favourable temperature distribution over the metal section, providing high quality of wire surface and microstructure [5,16]. It is also an advantage that an operator can feed stock material with different cross-sections into the same grooves. This ensures that a finished profile has the target size at different values of the total draw-down. This property demonstrates a great flexibility of four-roll pass compared to two-and three-roll grooves.  The layout of the equipment is shown in Figure 2 [13]. The stands are divided into roughing (with tension roll passes) and finishing (with finishing passes). The main work on total deformation of the metal is carried out in the tension roll passes, while the finishing passes are used to form the geometry and size of the finished profile. Table 1 shows options for four-roll passes rolled on the mill.
The rolled stock is fed from a reel or a coil. There is a welding device for welding the ends of the wire rod during continuous rolling. The Wire Accumulator 5 eliminates shock loads on the decoiler if the turns get overlapped. If necessary, the High-Frequency Inductor 6 can heat the stock material to temperatures of warm rolling.

Description Diagram Draw-Down Ratio
Circle-square µ = 1.57 Circle-octagon µ = 1.11 Octagon-square µ = 1.66 Octagon-octagon µ = 1.17 Octagon-circle µ = 1.05 Constant change of angles during the transition from groove to groove, and extra heating due to intense deformations, create a favourable temperature distribution over the metal section, providing high quality of wire surface and microstructure [5,16]. It is also an advantage that an operator can feed stock material with different cross-sections into the same grooves. This ensures that a finished profile has the target size at different values of the total draw-down. This property demonstrates a great flexibility of four-roll pass compared to two-and three-roll grooves. The layout of the equipment is shown in Figure 2 [13]. The stands are divided into roughing (with tension roll passes) and finishing (with finishing passes). The main work on total deformation of the metal is carried out in the tension roll passes, while the finishing passes are used to form the geometry and size of the finished profile. Table 1 shows options for four-roll passes rolled on the mill.
The rolled stock is fed from a reel or a coil. There is a welding device for welding the ends of the wire rod during continuous rolling. The Wire Accumulator 5 eliminates shock loads on the decoiler if the turns get overlapped. If necessary, the High-Frequency Inductor 6 can heat the stock material to temperatures of warm rolling.

Description Diagram Draw-Down Ratio
Circle-square µ = 1.57 Circle-octagon µ = 1.11 Octagon-square µ = 1.66 Octagon-octagon µ = 1.17 Octagon-circle µ = 1.05 Constant change of angles during the transition from groove to groove, and extra heating due to intense deformations, create a favourable temperature distribution over the metal section, providing high quality of wire surface and microstructure [5,16]. It is also an advantage that an operator can feed stock material with different cross-sections into the same grooves. This ensures that a finished profile has the target size at different values of the total draw-down. This property demonstrates a great flexibility of four-roll pass compared to two-and three-roll grooves. The layout of the equipment is shown in Figure 2 [13]. The stands are divided into roughing (with tension roll passes) and finishing (with finishing passes). The main work on total deformation of the metal is carried out in the tension roll passes, while the finishing passes are used to form the geometry and size of the finished profile. Table 1 shows options for four-roll passes rolled on the mill.
The rolled stock is fed from a reel or a coil. There is a welding device for welding the ends of the wire rod during continuous rolling. The Wire Accumulator 5 eliminates shock loads on the decoiler if the turns get overlapped. If necessary, the High-Frequency Inductor 6 can heat the stock material to temperatures of warm rolling.

Description Diagram Draw-Down Ratio
Circle-square µ = 1.57 Circle-octagon µ = 1.11 Octagon-square µ = 1.66 Octagon-octagon µ = 1.17 Octagon-circle µ = 1.05 Constant change of angles during the transition from groove to groove, and extra heating due to intense deformations, create a favourable temperature distribution over the metal section, providing high quality of wire surface and microstructure [5,16]. It is also an advantage that an operator can feed stock material with different cross-sections into the same grooves. This ensures that a finished profile has the target size at different values of the total draw-down. This property demonstrates a great flexibility of four-roll pass compared to two-and three-roll grooves. The layout of the equipment is shown in Figure 2 [13]. The stands are divided into roughing (with tension roll passes) and finishing (with finishing passes). The main work on total deformation of the metal is carried out in the tension roll passes, while the finishing passes are used to form the geometry and size of the finished profile. Table 1 shows options for four-roll passes rolled on the mill.
The rolled stock is fed from a reel or a coil. There is a welding device for welding the ends of the wire rod during continuous rolling. The Wire Accumulator 5 eliminates shock loads on the decoiler if the turns get overlapped. If necessary, the High-Frequency Inductor 6 can heat the stock material to temperatures of warm rolling.

Description Diagram Draw-Down Ratio
Circle-square µ = 1.57 Circle-octagon µ = 1.11 Octagon-square µ = 1.66 Octagon-octagon µ = 1.17 Octagon-circle µ = 1.05 Constant change of angles during the transition from groove to groove, and extra heating due to intense deformations, create a favourable temperature distribution over the metal section, providing high quality of wire surface and microstructure [5,16]. It is also an advantage that an operator can feed stock material with different cross-sections into the same grooves. This ensures that a finished profile has the target size at different values of the total draw-down. This property demonstrates a great flexibility of four-roll pass compared to two-and three-roll grooves. The layout of the equipment is shown in Figure 2 [13]. The stands are divided into roughing (with tension roll passes) and finishing (with finishing passes). The main work on total deformation of the metal is carried out in the tension roll passes, while the finishing passes are used to form the geometry and size of the finished profile. Table 1 shows options for four-roll passes rolled on the mill.
The rolled stock is fed from a reel or a coil. There is a welding device for welding the ends of the wire rod during continuous rolling. The Wire Accumulator 5 eliminates shock loads on the decoiler if the turns get overlapped. If necessary, the High-Frequency Inductor 6 can heat the stock material to temperatures of warm rolling.

Description Diagram Draw-Down Ratio
Circle-square µ = 1.57 Circle-octagon µ = 1.11 Octagon-square µ = 1.66 Octagon-octagon µ = 1.17 Octagon-circle µ = 1.05 Constant change of angles during the transition from groove to groove, and extra heating due to intense deformations, create a favourable temperature distribution over the metal section, providing high quality of wire surface and microstructure [5,16]. It is also an advantage that an operator can feed stock material with different cross-sections into the same grooves. This ensures that a finished profile has the target size at different values of the total draw-down. This property demonstrates a great flexibility of four-roll pass compared to two-and three-roll grooves. Constant change of angles during the transition from groove to groove, and extra heating due to intense deformations, create a favourable temperature distribution over the metal section, providing high quality of wire surface and microstructure [5,16]. It is also an advantage that an operator can feed stock material with different cross-sections into the same grooves. This ensures that a finished profile has the target size at different values of the total draw-down. This property demonstrates a great flexibility of four-roll pass compared to two-and three-roll grooves.
This study is devoted to the analysis of the calculated and experimental rolling forces (rolling pressure of the metal). The paper analyses the results of calculating the mill torque according to the earlier developed analytical dependencies. The results of comparison with similar values obtained in the course of modelling are presented. The conclusions of these studies can be applied to any units where production lines include stands with multi-roll grooves.

Rationale for Model Development
The main task of rolling on industrial mills is to obtain a required product profile using minimum operations. At the same time, high productivity provision is also vital. Little attention is paid to equipment loads, metal temperature distribution, friction conditions, and other process parameters. Typically, in complex stands, the desired shape change is achieved by trial and error.
On the basis of the energy theory of the interaction of rolled stock and rolls, scientists of Chelyabinsk Polytechnic Institute under the guidance of Professor Vydrin V.N. developed methods for calculating the kinematic and power parameters for pressure-based metal processing in the studied stands. Technological processes have been developed for the grooving of round sections, the production of high-precision shaped profiles, and the hot rolling of low-ductility alloys in stands with four-roll passes [17,18].
The study of the stressed state of metal in complex deformation zones is covered in the work of Professor Kharitonov V.A. (Magnitogorsk State Technical University) [15,19,20]. In the research, the Professor used the Finite Element method implemented in the DEFORM-3D software package. In addition, this software was used to study the processes of wire rolling and drawing in the research [21,22]. It has been established that the combined method of rolling and drawing provides for a comprehensive processing of the wire cross section. This is because the monolithic die usually processes the surface layers, while the roller die processes the inner layers of the roll material.
Most experts agree that further improvement in wire rolling will be achieved through process control, as happened in sheet rolling production. This development would require the elaboration of budget-friendly high-precision models that can operate autonomously. The speed required for analysing the dynamics of the process indicates that relatively simple models for calculating the technological rolling mill parameters are needed as much as complex finite-element models. In essence, they should be the digital twins developed with due consideration to the results of empirical research. This task is of immediate interest due to the expanding assortment of mills with multi-roll grooves.

Materials and Methods
As it has been mentioned in Section 1, the research presented in the article, should provide solutions for the following problems: 1.
Improvement of the method applied to calculate the metal rolling pressure, moment and power of the drive motor.

2.
Making up a methodological base to simulate the rolling process at promising multi-stand mills, taking into account the electromechanical system interaction via metal.
To solve these problems, one does not have to conduct a detailed study of metal flow in axial directions and, correspondingly, to solve a 3D problem. That is why it is reasonable to apply one-dimensional model, based on the known dependences of rolling pressure and mill torque calculation. Currently, several dependences are known for determining the mill torque in stands with multi-roll grooves [6,[23][24][25][26]. Some of them are obtained for specific cases, others are more general. In this connection, there was a task of comparative analysis of dependencies in order to determine the most suitable parameters for continuous wire mills. They were developed under the following assumptions [6]: 1.
The hypothesis of plane sections, experimentally confirmed by many studies, is accepted as a working hypothesis. 2.
The relationship between normal and shear stresses is described either by the Amonton-Coulomb dry friction law or by Siebel's law.

3.
The hardening curve is approximated as a straight line. 4.
The friction coefficient along the contact arc is considered constant. 5.
The distribution of specific pressure over the groove width is assumed to be uniform. 6.
The contact arc is replaced by a chord.
The metal hardening curve along the deformation zone is complex, and it is difficult to consider in practical calculations [27,28]. At the calculation of rolling energy and force parameters, the dependence can be linearized. That is why at the derivation of the Equation (1), the hardening along the deformation zone is assumed to be linear [6]. Essentially, this dependence is a linear hardening kinematic model. Such dependences are obtained by the authors at the rolling of various profiles in laboratory conditions and at the existing rolling mill.
These assumptions helped to obtain an equation for determining the average specific rolling pressure of the metal, based on the condition of the equilibrium of forces acting on the elementary volume in the deformation zone [6]: where n-number of rolls in the groove; f -friction coefficient; α = arccos(1 − ∆h/2R)-entering angle; R-roll radius; ∆h-absolute drafting; σ y0 , σ y f -initial and final values of the yield point; the value σ y f depends on the drafting and can be found using the experimental curves given in [16]; µ = Q 0 /Q 1 -draw-down extraction; Q 0 , Q 1 -section of the stock material before and after drafting; σ 0 = F 0 /Q 0 , σ 1 = F 1 /Q 1 -specific rear and front tensions; F 0 , F 1 -rear and front tension; µ−1 -hardening module; Q ca -rolls' contact area, for a circle-square diagram: h 1 -strip height at the stand exit; µ γ -drawing-down in neutral section. Equation (1) was obtained by the monograph [6] authors for rolling at four-roll pass with all driven rolls at uniform drafting from all sides. As a working hypothesis, a plane section hypothesis is taken. Its truth is proved by the results of the research on the regularities of shape changing at the round billet (20 mm+ in diameter) deformation. Even at the increased diameter (up to 50 mm), the deformation of cross-sections is small. The friction coefficient in the arc of contact is taken as constant while the specific pressure distribution along the pass width is taken as uniform.
The cross-section area of the original workpiece Q 0 depends on its shape; the cross-section Q 1 at the roll exit depends on the shape of the groove. To find µ γ in [6], the recommended equation is where ε-relative deformation of the metal; the dependence for its calculation is provided below in the explanation of the Equation (16).
It should be noted that with (m − σ y0 ) < 0 (when rolling in the second, third and subsequent stands), this Equation is unsolvable. Therefore, it is advisable to use the Equation The authors of [6] recommend obtaining the rolling moment for the case of all driven rolls according to the dependence where l = √ R∆h-contact arc length. The authors of the publication [6] also obtained the dependence similar to the Equation (1), however, applicable at the hardening curve approximation by the second order dependence (parabola). The following Equations are recommended to use for its derivation: -for drawing-down in neutral section module; r-co-sign coefficient of q; δ-coefficient given by Equation (2).
-for calculating the average specific pressure of the metal on the rolls where C 0 , C 1 -integration constants: C 0 = σ y0 − σ 0 + m 1 (1 + δ)/2 − 2m 1 δµ; When deriving Equations (6)-(8), Siebel's law of friction was used. Metal forward slip in the deformation zone: The research in [24][25][26] studied the deformation zone using the energy conservation law and the least energy law. For the mill torque in multi-roll groove with all driven rolls, the following equations are obtained where σ ay = (σ y f + σ y0 )/2-average value of the yield point.
Obviously, the results obtained by applying the given formulas may differ significantly from each other and from the true values. However, the necessary analysis of the calculations accuracy has not yet been carried out. There are no recommendations in the literature on the preference of this or that dependency. Therefore, such an analysis is of immediate interest. To solve this challenge, one needs to know the reference parameters-using these for comparison would provide an opportunity to assess the accuracy of calculations. Such parameters are the results of the experiment. However, conducting a series of experiments for various grooves, profiles and rolling conditions is a difficult task. Therefore, the rational way is mathematical modelling. The use of methods of numerical integration provides more accurate application of the dynamic parameters of the deformation zone in comparison with calculations using analytical expressions. These will provide for consideration of the relationships and satisfactory precision of the results, while being less "science intensive" in comparison with the finite-element models in the DEFORM-3D package and others.

Digital Model for Deformation Zone of Stand with Four-Roll Pass
For the elementary volume of the deformation zone (Figure 3), the equation of contact stresses is written in the form [13] where p r , τ-normal (radial) and shear stress; Q x -cross-sectional area of metal with coordinate x; Z x -half the width of the metal contact with roll in the cross-section with the coordinate x; α x -angle, defining the position of the cross-section with the coordinate x; σ x -axial stress. The "+" sign refers to the forward slip zone, the "−" sign to the backward slip zone. To obtain the differential equation of contact stresses, it is advisable to use the hypothesis of plane sections and the plasticity equation in the form where σ dr -deformation resistance; σ 1 , σ 2 , σ 3 -principal stresses.
where r p , τ-normal (radial) and shear stress; x Q -cross-sectional area of metal with coordinate х; x Z -half the width of the metal contact with roll in the cross-section with the coordinate х; x αangle, defining the position of the cross-section with the coordinate x; x σ -axial stress.
The "+" sign refers to the forward slip zone, the "−" sign to the backward slip zone. To obtain the differential equation of contact stresses, it is advisable to use the hypothesis of plane sections and the plasticity equation in the form Further, for a symmetrical four-roll pass with all-driven rolls, in case of equal drafting by vertical and horizontal rolls and with small entering angle, the following equations are valid Taking into account Equations (12)- (14), Equation (11) takes the form The resulting equation serves as the basis for mathematical description of the deformation zone. Notably, its derivation did not take into account the tightening and spreading effect in the direction of the groove diagonal axis, or the inertial force arising from the accelerated motion of the elementary volume along the deformation zone, due to insignificance thereof [29].
The initial conditions are determined by the plasticity equations where σ dr0 , σ dr1 -the values of deformation resistance at the beginning and at the end of the contact arc ( Figure 4); σ A , σ B -stresses from tension forces at points A and B of the contact arc. The relationship between tangential and normal stresses is described by the dry friction law It is expedient to solve the differential Equation (15) according to the following algorithm. Normal stresses in the backward p r,back and forward p r, f orw slip zones are calculated with the step ∆x.
In the backward slip zone, the coordinate x 0 is determined by the equation where l x -contact arc projection on x axis; x 1 -x coordinate in the forward slip zone. volume along the deformation zone, due to insignificance thereof [29].
The initial conditions are determined by the plasticity equations σ , dr1 σ -the values of deformation resistance at the beginning and at the end of the contact arc ( Figure 4); A σ , B σ -stresses from tension forces at points A and B of the contact arc. The relationship between tangential and normal stresses is described by the dry friction law Comparison is performed at each step, p r, f orw and p r,back . If p r,back < p r, f orw , the calculation in the next step is carried out for the backward slip zone, if p r,back > p r, f orw -for the forward slip zone. As a result, a successive approximation to the neutral section is performed (in the neutral section the values of the normal stresses are maximum).
The pressure on the rolls and the mill torque are found by the Equations The average specific pressure of the metal on the rolls is determined by the equation x Q ca -horizontal projection of the roll/metal contact area.
The value of the deformation resistance, taking into account the work-hardening of the metal [16]: where α, β-empirical coefficients; σ s -deformation resistance; ε = Q csa −Q x Q csa ·100%; Q csa -original cross-sectional area of wire rod.
According to Equation (16), the derivative dσ dr /dx: dx . The equation for the forward slip in the deformation zone can be obtained using the equality of per-second volume passing through any section of the deformation zone where Q γ , γ-neutral area and angle. The definition of the contact arc projection l x , as well as the values Z x , tgα x , dQ x /dx, Q x , included in Equation (15), depends on the shape of the grooves formed by the rolls. Below, the authors consider the procedures for determining the calculation formulas using the example of the "circle-incomplete square" groove.

Groove System "Circle-Incomplete Square" (First Stand)
The original wire rod diameter is indicated by D; the square output size is indicated by a (Figures 4  and 5). Contact arc projection l x and tgα x are defined by equations They are derived from the geometric relationships that are explained by the above figures.
The equation for the forward slip in the deformation zone can be obtained using the equality of per-second volume passing through any section of the deformation zone where Q γ , γ -neutral area and angle.
The definition of the contact arc projection x l , as well as the values x Z , included in Equation (15), depends on the shape of the grooves formed by the rolls. Below, the authors consider the procedures for determining the calculation formulas using the example of the "circle-incomplete square" groove.

Groove System "Circle-Incomplete Square" (First Stand)
The original wire rod diameter is indicated by D ; the square output size is indicated by a (Figures 4 and 5). Contact arc projection x l and x tgα are defined by equations They are derived from the geometric relationships that are explained by the above figures. Equations for Q x and Z x are derived according to Figure 5. The cross-sectional area at the deformation zone entry (x = l x ) is equal to the cross-sectional area of the wire rod Area of output cross-section at x = 0: The areas of the triangles AOB and MOC are equal and determined by the equation Area of the COA sector: where Taking into account Equations (18) and (19), Equation (17) takes the form: Cross-section Q x with coordinate x (0 < x < l x ) is shown in Figure 5 by dotted lines. Its area is determined by the dependency where y-coordinate of the point M ( Figure 5) on the plane of metal contact with the rolls (Figure 4): Half the width of the contact plane Z x ( Figure 5): After differentiation of Equation (20) with respect to x and transformations, the following equation is obtained The above Equations are used to find the desired values for calculating the specific pressure according to dependency (15).
Similar results were obtained for the groove systems "incomplete square-irregular octagon" (second stand), "irregular octagon-irregular octagon" (third stand). Figure 6 provides a generalized structural diagram of the developed five-stand mill model shown in Figure 1 [14]. The considered above mathematical description makes it possible to establish the interconnection between metal rolling pressure P, mill torque M and metal forward slip in the deformation zone s in the form of complex non-linear functions of input and output sizes of the strip D and α, friction coefficient f, the yield strength of the initial billet σ S , forward and back tension values

Structure of a Digital Model for a Five-Stand Mill
For the purpose of linearization of these equations in the operational range of coordinate change, the authors transit from the absolute values P i , M i , s i to their increments. In this case, the system (21) takes the form where c 1,i . . . c 17,i -technological coefficients (partial derivatives) , The mathematical description of interstand space includes the interconnection of tension in the i-th space F i (after the i-th stand) and strip elastic elongation. The Equation (22) is written in compliance with Hooke's law [30] F The interconnection between the tension and metal velocity in the interstand space is described by the dependence where E-strip elasticity modulus; Q i -strip cross-section in the i-th space; E i -relative elastic elongation in the i-th spacee; l i -distance between the axes of rolls in adjacent stands. The mathematical descriptions of the system "Frequency converter-asynchronous motor (FC-AM)" as well gear boxes in the implemented model are similar to the known ones [30][31][32].
The developed dynamic model of a continuous rolling mill with four-roll passes as control object taking into account electric drive interconnection by means of rolled metal has been implemented in the form of the SIMULINK visual programming software modules of the software package MATLAB 6.0 [33]. Here the authors do not provide a detailed description of the software. Due to the abundance of information, it can be covered in the scope of a separate publication.
The results of the energy and force rolling parameters simulation are considered below in the Section 6.1. Figure 6. Structural diagram of a mathematical model of a continuous five-stand mill with four-roll passes. Figure 6. Structural diagram of a mathematical model of a continuous five-stand mill with four-roll passes.

Initial Rolling Product Range
The experimental wire rolling in the stands with multi-roll passes was conducted from the steel grades and alloys given in the Table 2. Their rolling is well-reasoned as drawing is either allied with process complexities or characterized by low economic indicators.
To develop the technology and define energy and force parameters for rolling high-strength reinforcement wire, the authors used medium carbon and high carbon steel 60, 70 and 80. During this process the authors estimated a possibility to obtain reinforcement wire from the steel with a lower carbon content and increased cold hardening by means of an overall reduction increase to improve structural behavior.
The experimental research was conducted in laboratory conditions at a double-stand mill by means of rolling short samples as well as at continuous bundle rolling. The research was conducted in the conditions matching the industrial ones as closely as possible. The wire of the 2.5-4 mm diameter was obtained by three passes of 6.5-12 mm rod coil. To comply with the temperature mode occuring at continuous rolling, the intermediate billet was heated before every next pass up to the temperature the metal had after the last stand during the previous pass. Alongside energy and force parameters, the authors measured the following process parameters: rolled stock dimensions after the pass, drawing down for a pass, metal slip, temperature.
By the measurement results, the authors calculated the degree of pass filling, and measured the cross-section dimensions on the pass connectors and the metal heating per passage. For the treatment of experimental data, mathematical statistics methods were used.
The change in contact friction was reached by the application of oil or chalk on the roll operating surface. Additionally, the authors conducted rolling with the emulsion on dry rolls washed with acetone. According to the data [34], the value of the friction coefficient on the rolls lubricated with industrial oil was f = 0.088; at the rolling on dry rolls, f = 0.22. At the change of the friction coefficient in this range, the pressure increment equaled 10% ( Figure 7a). As, due to the change of the emulsion content, the friction coefficient will fluctuate significantly less during one bundle rolling, it can be reasonably concluded that the value of pressure change virtually does not depend on this factor.

Measuring Energy and Force Rolling Parameters
Depending on the stand structure to measure a full rolling pressure, the authors applied various methods [35]:

1.
Defining the deformation of measuring pins installed in the chock body under roll bearing.
Operating and compensating sensors were put on the pin polished edges; the sensors were connected in a half-bridge circuit, four sensors in each leg.

2.
Stand housing deformation measurement. To do this, the authors measured deformations occurring in the roll base and determined the reactions occurring in the bearing assembly. For this purpose, tension meters were put on the bases. They were intended to measure the resulting force of rolling pressure and its horizontal component. The calibration of pressure meters was conducted in the stands with special hydraulic devices allowing for the simultaneous loading of all rolls with equal forces simulating rolling pressure. 3.
The rotating moment was obtained by the motor armature current (at the experimental plant DC current motors are applied). To measure the tension in the interstand space, two methods were used at continuous rolling: measurement of stand movement from the axial force acting on the strip side and tension measurement with the help of a roller type sensor.

4.
Sample temperature at cool and hot rolling was measured by a double-electrode thermocouple with information record. The temperature of rolled stock surface at continuous rolling measured by a photoelectric pyrometer (non-contact method).

Experimental Research Results
The drawing down at rolling in the system "circle-tetragon" is within the range µ = 1.2-1.6. The dependences P = f (µ) at such drawing-downs are provided in Figure 7c, and they are virtually linear. Note that at the increased rolling rate, due to the reduced emulsion amount on the wire, the rolled stock temperature grows, thus causing less rolling pressure of the metal (Figure 7b). This occurs due to heat release at plastic deformation. However, the initial heating has a more significant impact ( Figure 7d); thus, at the temperature of 500 °С, one observes the pressure decreasing by 20-30%.

Mathematic Simulation Result
The developed mathematical model provides details of the distribution of contact stresses along the contact arc and the groove width, as well as the values of the average specific pressure, rolling pressure, mill torque, slip, profile dimensions. For its practical implementation, the authors worked out an algorithm and a program for calculating the above parameters. The calculations of power parameters for rolling profiles of various assortments on a five-stand mill have been performed. The correspondence of the calculated and experimental results has been checked for several profiles.
As an example, below, the authors consider the results of studies for rolling a wire with 11 mm diameter, from 10KP steel, with the "circle-incomplete square" groove system. The mill was not designed for rolling such a profile, but the task was mastered in process of its exploitation. During the experiments, a total rolling pressure of the metal, the mill torque, the drawing-down and the metal flow in the deformation zone were measured. To determine the draw-down, the initial samples were marked every 100 mm and then assessed after the rolling with an instrumental micrometer to determine the post-rolling elongation. The metal flow in the deformation zone was determined by the method of core imprints [36,37].
The friction coefficient for analytical calculations and the hardening curve were determined empirically. For rough rolls, the friction coefficient is 1.2. The hardening curve is described by the equation . Note that at the increased rolling rate, due to the reduced emulsion amount on the wire, the rolled stock temperature grows, thus causing less rolling pressure of the metal (Figure 7b). This occurs due to heat release at plastic deformation. However, the initial heating has a more significant impact ( Figure 7d); thus, at the temperature of 500 • C, one observes the pressure decreasing by 20-30%.

Mathematic Simulation Result
The developed mathematical model provides details of the distribution of contact stresses along the contact arc and the groove width, as well as the values of the average specific pressure, rolling pressure, mill torque, slip, profile dimensions. For its practical implementation, the authors worked out an algorithm and a program for calculating the above parameters. The calculations of power parameters for rolling profiles of various assortments on a five-stand mill have been performed. The correspondence of the calculated and experimental results has been checked for several profiles.
As an example, below, the authors consider the results of studies for rolling a wire with 11 mm diameter, from 10KP steel, with the "circle-incomplete square" groove system. The mill was not designed for rolling such a profile, but the task was mastered in process of its exploitation. During the experiments, a total rolling pressure of the metal, the mill torque, the drawing-down and the metal flow in the deformation zone were measured. To determine the draw-down, the initial samples were marked every 100 mm and then assessed after the rolling with an instrumental micrometer to determine the post-rolling elongation. The metal flow in the deformation zone was determined by the method of core imprints [36,37].
The friction coefficient for analytical calculations and the hardening curve were determined empirically. For rough rolls, the friction coefficient is 1.2. The hardening curve is described by the equation σ x = 41.8 + 3.92ε 0.506 (23) As mentioned above, the hardening curves were obtained experimentally at the rolling of various profiles. For each profile, the experimental non-linear dependence is approximated by the equation similar to (23). Such an approximating equation is provided in the guides [16,29] and other references. Therefore, the authors apply more precise kinematic hardening models as compared with linear approximation dependences.
The calculated and experimental dependencies during tension-free rolling are shown in Figure 8. Legend: Curve 1 (dashed) is plotted according to experimental data; Curve 2 is obtained using a mathematical model. As mentioned above, the hardening curves were obtained experimentally at the rolling of various profiles. For each profile, the experimental non-linear dependence is approximated by the equation similar to (23). Such an approximating equation is provided in the guides [16,29] and other references. Therefore, the authors apply more precise kinematic hardening models as compared with linear approximation dependences.
The calculated and experimental dependencies during tension-free rolling are shown in Figure  8. Legend: Curve 1 (dashed) is plotted according to experimental data; Curve 2 is obtained using a mathematical model.   (1) and (3); Curve 4 was calculated using Equations (3) and (7). Figure 8b shows the dependencies of metal flow in the deformation zone: Curves 3 and 4 are calculated by Equations (4) and (6), respectively. Figure 8c shows similar dependencies for the mill torque: Curves 3 and 4 are calculated according to Equation (5); Curves 5 and 6-according to Equations (9) and (10); Curve 7-according to Equation (8).
Calculations and experimental studies were carried out for same-profile wire-rod milling when the front tension changed from zero to 1500 kgf (15 kN The hardening curve was calculated using Equation (21). The results are shown in Figure 9; the numbering of curves is the same as in Figure 8.   (1) and (3); Curve 4 was calculated using Equations (3) and (7). Figure 8b shows the dependencies of metal flow in the deformation zone: Curves 3 and 4 are calculated by Equations (4) and (6), respectively. Figure 8c shows similar dependencies for the mill torque: Curves 3 and 4 are calculated according to Equation (5); Curves 5 and 6-according to Equations (9) and (10); Curve 7-according to Equation (8).
Calculations and experimental studies were carried out for same-profile wire-rod milling when the front tension changed from zero to 1500 kgf (15 kN). The rear tension was assumed to be zero. Rolling was carried out at values µ = 1.316 and f = 0.12. The hardening curve was calculated using Equation (21). The results are shown in Figure 9; the numbering of curves is the same as in Figure 8. Figure 10 shows similar dependences obtained when changing the rear tension in the same range at zero front tension.
Several observations can be made from analysing the results: 1.
The considered analytical expressions and the developed digital model reliably reflect the patterns of the metal deformation process in a multi-roll groove. By the nature of the change, they correspond to the results of experimental studies. However, the degree of impact of the process parameters on the rolling pressure, the mill torque and the metal flow are different.

2.
The closest coincidence of the calculated and experimental data is provided by the mathematical model. In the assessed cases, the discrepancy ceiling is 10-15%.

3.
From the analytical dependencies considered above, Equation (10) can be recommended for determining the mill torque. It provides acceptable results (the discrepancy at µ = 1.3-1.5 does not exceed 20%) and requires minimal calculations.

4.
It is advisable to calculate the metal flow in the deformation zone using the mathematical model only. The use of analytical formulas gives a result that is inflated by 2-5 times.
Machines 2020, 8, x FOR PEER REVIEW 18 of 22 Figure 10 shows similar dependences obtained when changing the rear tension in the same range at zero front tension.
Several observations can be made from analysing the results: 1. The considered analytical expressions and the developed digital model reliably reflect the patterns of the metal deformation process in a multi-roll groove. By the nature of the change, they correspond to the results of experimental studies. However, the degree of impact of the process parameters on the rolling pressure, the mill torque and the metal flow are different. 2. The closest coincidence of the calculated and experimental data is provided by the mathematical model. In the assessed cases, the discrepancy ceiling is 10-15%. 3. From the analytical dependencies considered above, Equation (10) can be recommended for determining the mill torque. It provides acceptable results (the discrepancy at µ = 1.3-1.5 does not exceed 20%) and requires minimal calculations. 4. It is advisable to calculate the metal flow in the deformation zone using the mathematical model only. The use of analytical formulas gives a result that is inflated by 2-5 times.

Discussion of the Results
The comparative analysis of the processes shown in Figures 9 and 10 shows that an increase in the wire tension (either front or rear) decreases the rolling force. At the same time, the rear tension has a greater impact on the pressure than the front tension. They also have the opposite effect on the mill torque. With an increase in the front tension from 0 to 1500 kgf (Figure 9b), the torque decreases from ~150 kg•m (1.5 kN•m) to zero. This leads to a decrease in the rolling force (Figure 9a), which is quite understandable from the viewpoint of the process "physics".
An increase in the rear tension leads to an increase in the mill torque (Figure 10b), although the Figure 10. Dependencies similar to Figure 9, (a), torque (b) and flow (c) from rear tension.

Discussion of the Results
The comparative analysis of the processes shown in Figures 9 and 10 shows that an increase in the wire tension (either front or rear) decreases the rolling force. At the same time, the rear tension has a greater impact on the pressure than the front tension. They also have the opposite effect on the mill torque. With an increase in the front tension from 0 to 1500 kgf (Figure 9b), the torque decreases from 150 kg·m (1.5 kN·m) to zero. This leads to a decrease in the rolling force (Figure 9a), which is quite understandable from the viewpoint of the process "physics".
An increase in the rear tension leads to an increase in the mill torque (Figure 10b), although the rolling force (Figure 10a) decreases in this case. This is because the rear tension creates a torque that prevents deformation. The increase in the torque of the stand motor provides its compensation.
The results of torque comparison in Figures 6-8 confirm the relevance of the development of a tension control system in the inter-stand spaces of a continuous wire mill with multi-roll grooves. Some results obtained in the development of such a system for cold rolling mills are given in publications [38][39][40].
Similar systems have been developed for a continuous subgroup of the roughing group of stands of a broad-strip hot-rolling mill [41,42]. The results of their studies are presented in [43][44][45][46]. However, for mills with multi-roll groove, the solution to this challenge is more complicated than the solution for mills with two-high rolling. Here, the distinctive feature lies in the fact that direct measurement of tension on mills like the one considered is practically impossible-at the same time, its calculation by the parameters of the electric drive is complicated due to the constant change of groove and roll profiles [47][48][49].
That said, the development of such a system is underway, and the research results will be published. The results of experimental research were applied to check the correctness of theoretical conclusions, calculate and design process, electrical equipment and control systems. This material is abundant and cannot be covered in the scope of one article. It may be considered in the further publications made by the authors.

1.
The authors carried out a complex of theoretical and experimental studies of power parameters for rolling mills with four-roll pass stands. The obtained results can be used as the initial data for development of new assortments for rolled products or in the design of new mills. The paper analyses the accuracy of determining the power parameters by the known analytical dependences, giving recommendations for their application in calculations.

2.
The developed digital model of the deformation zone provides more accurate calculations in comparison with analytical expressions. The discrepancy between the results from the model and from the experimental data does not exceed 10-15%.

3.
When calculating the power of actuating units, the mill torque should be defined using the digital model. For approximate calculations, it is advisable to use Equation (10) as it provides results with acceptable precision while requiring the least amount of calculations.

4.
The results of experimental research are recommended for the application to optimize drafting, and define energy and power and kinematic parameters and the temperature mode impact at the rolling of the extended range of billets at the existing multi-roll mills. 5.
The developed model and program for calculating power parameters are recommended for use in studies of rolling on "conventional" cold rolling mills.

Funding:
The work was carried out with the financial support of the Ministry of Science and Higher Education of the Russian Federation within the framework of a subsidy for financial support for the fulfillment of a state task (fundamental scientific research), contract No. FENU-2020-0020 (2020071GZ).