Open Access Library Journal
Vol.04 No.01(2017), Article ID:73734,17 pages
10.4236/oalib.1103247

Novel Analytic Model for the Projected Contact Zone Based on the Flow Line Element Method in Alloyed Bar Rolling

Yonggang Dong, Heng Zhu, Jianfeng Song*

College of Mechanical Engineering, Yanshan University, Qinhuangdao, China

Copyright © 2017 by authors and Open Access Library Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

Received: December 21, 2016; Accepted: January 20, 2017; Published: January 23, 2017

ABSTRACT

For obtaining the profile and the projected area of contact zone exactly during alloyed bar rolling by Round-Oval-Round pass sequence, the analytic model for the length contact and contact boundary curve were built firstly by considering the influence of the spread of the outgoing workpiece on the effective height of outgoing workpiece and roll mean radius, and then the contact surface was discretized by finite flow line elements. Moreover, the radius equation and bite angle equation of different flow lines were derived and they were all expressed as the function of the position angle, then any flow line on the contact surface can be determined since the position angle has been given. Finally, since the analytic equation for the projected area of contact surface was hard to be integrated directly, the analytic model was proposed by summing up the area of discretized parts on the contact surface. Based on the analytical model of contact boundary and flow line element, 3-dimension contact surface was rebuilt by mathematical software, and the validity of analytic model was examined by the bar rolling experiments and the numerical simulation of alloyed bar rolling by rigid-plastic FEM software. Compared with the existing models, the precision of the projecting area of contact zone was improved obviously. So, it can be applied in alloyed bar rolling to predict the projected area of contact zone and velocity of outgoing workpiece exactly.

Subject Areas:

Mechanical Engineering

Keywords:

Alloyed Bar, Contact Zone, Projected Area, Round-Oval-Round, Flow Line Element

1. Introduction

Compared with products rolled with the oval-square-oval pass sequence, the bar or rod rolled with round-oval-round pass sequence has a better surface quality and mechanical performance, and the round-oval-round pass sequence is the most common roll pass in bar or rod continuous rolling recently. Since the characteristics of non-uniform distribution of stain, stress and flow velocity on the deformation zone, it is difficult to analyze the process of alloyed bar rolling in oval-square-oval pass sequence accurately [1] [2] [3] .

In the past few years, the research on the alloyed bar rolling process was carried out by the simulation methods and experiments. References [4] [5] [6] [7] [8] studied the strain, stress and geometry of deformation zone of bar rolling by the numerical simulation methods based on the FEM software, and references [9] [10] [11] analyzed the effect of roll gap, roll profile and rolling speed on the wear and exit section area by the single-pass rolling and multi-pass rolling experiments and FE analysis. References [12] [13] [14] studied the rolling force and rolling torque of bar rolling and proposed analytic models for estimating the force energy parameters. Because the contact surface between the rolling workpiece and shaped roll is not a cylindrical surface like flat roll rolling but a three- dimension curve with complex boundary, the length and boundary of contact surface has to be simplified without taking the spread of rolling workpiece into account for obtaining projected area easily. Based on these simplifications, the geometry and projected area of contact zone were calculated by the conventional equation for billet rolling and the roll force and roll torque was derived on these results, and the results of roll force and roll torque were just an approximate value for bar rolling in round-oval-round pass sequence. Furthermore, the profile of contact boundary and the projected area of contact zone influence the results of velocity field distribution, stress field distribution and strain field distribution directly. Therefore, it is indispensable for analyzing the bar or rod continuous rolling process to determine the profile of contact boundary and the projected area of contact zone exactly.

It is not easy to obtain the precise projecting area of contact zone and the accuracy profile of contact boundary since the contour of contact zone is not regular and difficult to be defined. So, at present the geometry and projected area of contact zone is mostly calculated by the simplified methods, such as the empirical equation [2] , the graphical solving method [3] .

Shinokura and Takai [15] [16] [17] gave the spread formula and projected area of the deformed workpiece in round-oval-round pass sequence. However, this formula was derived on the base of rolling experiments of plain carbon steel. Y. Lee proposed a mathematical model for predicting the surface profile of deformed workpiece and mean roll radius in Round-Oval-Round pass sequence [18] [19] [20] .

The profile of contact boundary and the projecting area of contact zone are correlated to the solution of the rolling force, the velocity field distribution, stress field distribution and strain field distribution directly. So, it is very important to determine the profile of contact boundary and the projected area of contact zone exactly.

2. The Existing Models for Determining the Projecting Area of Contact Zone

2.1. Graphical Solving Method [3]

As can be seen in Figure 1, the boundary curve of contact zone was simplified as a intersection line between the revolution surface of roll pass and the profile of incoming workpiece, and the intersection line was obtained by linking the intersect in point one by one. So the projection of contact zone was shown as the hatching zone.

2.2. Analytic Methods Based on the Two Hypothetic Curve Functions for the Contact Boundary

Shinokura and Takai [15] proposed a formula to estimate the maximum contact length L and projected area Ap of contact zone in alloyed bar rolling by ignoring the spread of outgoing workpiece.

The size of oval groove, round groove and corresponding incoming workpiece were shown in Figure 2(a) and Figure 2(b) respectively. The equivalent height of incoming workpiece H 0 ¯ and equivalent height of outgoing workpiece H m 0 ¯ were simplified as

H 0 ¯ = A 0 A s 0 2 C Z 0 (1)

H m 0 ¯ = A 0 A s 0 A h 2 C Z 0 = A e 0 2 C Z 0 (2)

where A 0 is the section area of incoming workpiece, A s 0 , A h , A e 0 are the non-effective reduction area, effective reduction area and effective exit section

Figure 1. The projected area by a graphical solving method.

Figure 2. Definition of pass profile and the incoming workpiece. (a) Round-oval pass; (b) Oval-round pass.

area respectively, ( C y 0 , C z 0 ) is the intersection point between the profile of incoming workpiece and the roll pass.

In round-oval pass rolling, A h and A s 0 may be obtained by

A h / 2 = arctan ( C Y 0 / C Z 0 ) R a 2 [ - C Y 0 C Y 0 ( R 1 2 y 2 D z ) d y C Y 0 C Z 0 ] (3)

A s 0 / 2 = arctan ( C Z 0 / C Y 0 ) R a C Y 0 C Z 0 . (4)

In oval-round pass rolling, A h and A s 0 may be given by

A h / 2 = π 8 W max H p A s 0 2 arccos ( C Z 0 / R Y 0 ) R g 2 C Y 0 C Z 0 (5)

A s 0 / 2 = arcsin ( C Z 0 / R 1 ) R 1 2 ( D Z + C Y 0 ) C Z 0 (6)

where R 1 is the radius of the oval groove, R s is the radius of the curvature of the incoming cross-section, D Z is the distance along the Z-axis direction between the origin coordinate

Then the maximum contact length L max was obtained by

L max = R m 0 ( H 0 ¯ H m 0 ¯ ) = R m 0 A h 2 C Z 0 (7)

where R m 0 is the mean roll radius at the entrance-section, it can be shown as

R m 0 = R min + H p 2 H m 0 ¯ 2 . (8)

For obtaining the projected area of contact zone in alloy bar rolling, Shinokura and Takai expressed the contact boundary curve by the function L max 1 x 2 C y 2 according to the empirical data, and the projected area was shown as

A p = 0 C y L max 1 x 2 C y 2 d x = π 2 L max C y . (9)

As shown in Figure 2, C y is the coordinate of critical point on the contact boundary along the spread direction. ( C y , C z ) is the coordinates of the critical point on the contact boundary at the exit-section, and it can be determined by the reference [16] [17] .

On the base of Equation (9) Y.lee given another contact boundary curve function L max ( 1 x C y ) 1 / m to modify this equation and the projected area was shown as

A p = 0 C y L max ( 1 x C y ) 1 / m d x = 3 2 L max C y ( m = 1 3 ) . (10)

Moreover, as shown in Equation (9) and Equation (10), two hypothetical functions for indicating the contact boundary curve were given directly without any reasoning and any derivation process. Although the error between the results of two equations are not obvious when the size of rolling workpiece is small enough and then the contact length L and Cy is small enough, the absolute error of these two equations will be considerable and it should not be ignored in large diameter bar rolling. So it is not precise enough for these two semi-analytic models to calculate the projected area of multi-pass alloyed bar rolling, then an analytic model should be built to predict the projected area accurately.

In summing up these models for projected area, the Equation (2), which is based on the graphical solving method, does not take the influence of the spread and the contact boundary status of deformed workpiece into account. Therefore, the Graphical solving method can be just used as an estimating value when the spread of deformed workpiece is small enough (Figure 2).

3. A Novel Analytic Model for the Projected Area of Contact Surface

3.1. Modification of Contact Length Model

As shown in Figure 3, if the spread of outgoing workpiece was not negelected, the effective section area A e and the equivalent width 2 C y , height H m ¯ and mean roll radius R m are totally different from the corresponding parameters A e 0 , 2 C y 0 , H m 0 ¯ and R m 0 in Equation (2). ( C y , C z ) can be calculated by references [19] [20] . So, the contact length modified as

L m = R m ( H 0 ¯ H m ¯ ) = ( R min + H p / 2 G / 2 H m ¯ / 2 ) ( H 0 ¯ H m ¯ ) . (11)

In round-oval pass rolling, the equivalent height of outgoing workpiece was shown as

H m ¯ 2 = A e 4 C y = C y C y ( R 1 2 y 2 D z ) d y C y = R 1 2 ( 2 C y R 1 cos ( arcsin C y R 1 ) + 2 arcsin C y R 1 ) 4 C y D z . (12)

In oval-round pass rolling, the equivalent height of outgoing workpiece was

Figure 3. Parameters in modified model of contact length. (a) Round-oval pass; (b) Oval-round pass.

expressed as

H m ¯ 2 = A e 4 C y = arctan ( C y / C z ) R + g 2 C y C z C y . (13)

3.2. Deformation Zone Geometry and the Boundary of Contact Zone

As can be seen in Figure 4 and Figure 5, coordinate axes x, y, z are chosen to be the directions of bar length, width and height, respectively, with the origin of the x-axis at the midpoint of the entry plane. The x-y and x-z are plane of symmetry. The curve equation of contact boundary S can be obtained approximately by interpolating between the point ( 0,0, H 0 ) and the point ( L , C y , C z ) . The coordinates of any point on the contact boundary S are shown as ( x s , C y s , C z s ) .

3.3. The Definition of the Flow Line and Flow Line Element on the Contact Surface

The flow plane was defined as a set of eccentric continuous cylindrical surfaces having almost straight generators parallel to the roll axis. Assuming that at any cross section along the roll bite the bar height deformation is uniform, the deformation zone was constituted by a set of flow plane, and the contact surface was constituted by a set of flow line element f α . The flow line element f α is a set of concentric circular arc with different radius R and different bite angle θ α s , which center is attached on the roll axis. Since the roll radius in the roll pass is different and the height of incoming workpiece is different along the y-axis direction, the bite angle θ α s , along the whole contact boundary, is not a constant but a variable which changes with a different position angle α or a different roll radius R .

Figure 4. Geometry and flow line element on the contact surface in oval pass. (a) Geometry of contact surface; (b) Flow line on the contact surface.

Figure 5. Geometry and flow line element on the contact surface in round pass. (a) Geometry of contact surface; (b) Flow line on the contact surface.

3.3.1. The Radius of Flow Line R and the Bite Angle θαs in the Round-Oval Pass Sequence

The radius of flow line f α was expressed as

R = R min + R 1 ( 1 cos α ) ( 0 α arcsin C y R 1 or 0 α arccos C z + D z R 1 ) . (14)

The three-dimension coordinates of a random point on the contact surface was shown as

{ x = L R sin θ ( 0 θ θ α s ) y = R 1 sin α z = R ( 1 cos θ ) + R 1 cos α D z = R cos θ + R c . (15)

On the symmetry plane of deformed workpiece, the position angle α is 0 and the bite angle θ α s reaches the maximum value θ max

θ max = arccos ( 1 H 0 H p 2 R 0 ) . (16)

At the exit section of deformed workpiece, the bite angle θ α s is 0 and the position angle α reaches the maximum value α max

α max = arccos C z + D z R 1 (17)

According to the Equation (15), the coordinates of points on the contact boundary was expressed as

{ x s = L R sin θ α s y s = R 1 sin α z s = R ( 1 cos θ α s ) + R 1 cos α D z = R cos θ α s + R c . (18)

In the contact zone, the height of the profile at y = 0 along the x-direction is expressed as

H ( x ) = H 0 2 L R 0 x + 1 R 0 x 2 . (19)

The curve equation of contact boundary S can be obtained approximately by interpolating between the point ( 0 , 0 , H 0 ) and the point ( L , C y , C z ) . The coordinate of any point on the contact boundary S is shown as ( x , C y s , C z s ) , and the coordinate C z s can be shown as

2 C z s = H 0 2 ( H 0 2 C z ) x s L + ( H 0 2 C z ) x s 2 L 2 . (20)

Substituting the xs of Equation (18) into Equation (20) yields

2 C z s = ( H 0 2 C z ) R 2 sin 2 θ α s L 2 + 2 C z . (21)

According to the equation C z s = z s and sin 2 θ α s = 1 cos 2 θ α s yields

( H 0 2 C z ) R 2 cos 2 θ α s L 2 + 2 R cos θ α s + ( H 0 2 C z ) R 2 L 2 2 ( R c C z ) = 0 . (22)

Solving the Equation (22) yields

θ α b = arccos L 2 + L L 2 + ( H 0 2 C z ) [ ( H 0 2 C z ) R 2 L 2 2 ( R c C z ) ] R ( H 0 2 C z ) . (23)

Substituting Equation (14) into Equation (23) yields

θ α s = arccos L 2 + L L 2 + ( H 0 2 C z ) [ ( H 0 2 C z ) ( R min + R p ( 1 cos α ) ) 2 L 2 2 ( R c C z ) ] ( R min + R 1 ( 1 cos α ) ) ( H 0 2 C z ) . (24)

3.3.2. The Radius of Flow Line R and the Bite Angle θαs in the Oval-Round Pass Sequence

For the oval-round pass sequence, the radius of flow line f α was expressed as

R = R min + R g ( 1 cos α ) ( 0 α arcsin C y R g or 0 α arccos C z R g ) . (25)

The three-dimensional coordinates of a random point on the contact surface was shown as

{ x = L R sin θ ( 0 θ θ α s ) y = R g sin α z = R ( 1 cos θ ) + R g cos α = R cos θ + R c . (26)

On the symmetry plane of deformed workpiece, the position angle α is 0 and the bite angle θ α s reaches the maximum value θ max

θ max = arccos ( 1 H 0 H p 2 R 0 ) . (27)

At the exit section of deformed workpiece, the bite angle θ α s is 0 and the position angle α reaches the maximum value α max

α max = arccos C z R g . (28)

By the same methods as round-oval pass sequence, the bite angle of oval- round pass sequence was

θ α s = arccos L 2 + L L 2 + ( H 0 2 C z ) [ ( H 0 2 C z ) ( R min + R g ( 1 cos α ) ) 2 L 2 2 ( R c C z ) ] ( R min + R g ( 1 cos α ) ) ( H 0 2 C z ) . (29)

3.3.3. The Projected Area of Contact Zone

According to definition of projected area of contact zone, it should be calculated by the integral equation as

A p = 2 0 α max R sin θ α s d y = 2 0 α max R sin θ α s R 1 cos α d α . (30)

However, it is too complex for Equation (30) to integrate it directly and get a function for projected area. So, the flow line field method based on the discrete law is an effective way to calculate projected area.

Once the position angle of flow line α is given, the roll radius R and the bite angle θ α s against different position angle of flow line can be determined. Then the space position and the length of any flow line on the contact surface can be obtained. If the whole contact zone was discretized into n flow lines with different position angle α i and different bite angle θ α s i , the arc length of any flow line was expressed as

L f i = R i θ α s i = f ( α i ) ( 0 i n ) . (31)

The position angle of ith flow line can be shown as

α i = i n α max ( 0 i n ) . (32)

The projected length of this flow line on plane xoy can be obtained by

L p i = R i sin θ α s i ( 0 i n ) . (33)

Substituting Equation (32) into Equation (14), Equations ((24), (25) and (29)) to replace the position angle α can yield R i and θ α s i in oval pass rolling and round pass rolling respectively, then the projected area of contact zone on the plane xoy can be shown as

A p = 2 0 α max R sin θ α s d y = i = 0 n ( R i sin θ α s i ) C y n . (34)

4. Results and Discussions

The alloyed bar rolling experiments had been accomplished in BEIMAN SPICIAL STEEL CO. LTD, the round workpiece of diameter 171 mm were rolled in one oval pass and one round pass of 22-stand Pomini Rolling Mill. The deformation zone of rolling workpices was obtained by stopping the rolling process when the workpiece was rolled in the oval pass and round pass simultaneously. The rolling schedule is shown as Table 1. The material is structural alloyed steel 41Cr4. The dimension of the groove schedule is respectively shown as Figure 6(a) and Figure 6(b).

As can be seen in Table 2, the parameters for calculating the projecting area of

Table 1. Rolling schedule of the Pomini Rolling Mills.

Table 2. The parameters for calculating the projecting area.

Figure 6. The dimension of the groove schedule. (a) Oval pass (b) Round pass.

Figure 7. Profile of the curve on the contact boundary. (a) Oval-pass; (b) Round-pass.

contact zone were listed one by one.

As shown in Figure 7(a) and Figure 7(b), the profile of curve on the contact boundary in oval pass and round pass were measured and distinguished by a white chalk.

As shown in Figure 8, the contact boundary and the contact surface were obtained by the rigid-plastic FEM software DEFORM5.03. Moreover, the surface profile of incoming workpiece and outgoing workpiece in oval pass and round pass was shown as Figure 9.

Figure 8. Contact boundary and contact zone in oval pass rolling.

Figure 9. Surface profile of incoming workpiece and outgoing workpiece.

As can be seen in Figure 10(a) and Figure 10(b), the 3-dimension contact surface and contact boundary were rebuilt by Matlab 7.0 according to the novel analytic model of contact boundary and contact surface. In Figure 11, the contact boundary from the novel analytic model was compared with that of experimental data and simulation result.

Results of contact length and projected area from the novel model, the calculating results from the existing models, the experimental data and the simulation results were all listed in Table 3. As shown in Table 3, the results of contact length and projected area from different models were listed and compared with corresponding experimental data and simulation results. The results of contact length from the modified analytic model are less than that of empirical formula and the graphical solution, and it is greater than that of Shinokura formula and Y. Lee formula. Moreover, the error of contact length from the modified model is less than that of existing models.

Since the outgoing workpiece of oval pass rolling will be rolled in next round pass as an incoming workpiece and the section profile at the exit of oval pass influence the contact surface of round pass rolling greatly, the prediction error of round pass is obviously greater than the prediction error of oval pass. Moreover,

Figure 10. Rebuilt contact surface based on analytic model by Matlab7.0. (a) Oval-pass; (b) Round-pass.

Figure 11. Profile of the curve on the contact boundary. (a) Oval-pass; (b) Round-pass.

Table 3. Results of different models.

results of projected area from the novel analytic model approaches the experimental data and simulation results very well, and its error is less than any existing models.

5. Conclusions

1) The contact boundary is a complex 3-dimension curve, and its profile is not only concerned with the parameters of pass profile R1, Rg, Rmin, G, Dz and the shape parameters of incoming workpiece H0, Ra, but also influenced by the coordinates of critical point ( C Y , C Z ) ;

2) The modified contact length model is rational because the influence the effective section area of the outgoing workpiece Ae, the critical point ( C Y , C Z ) on the contact boundary, the effective height of outgoing workpiece H m ¯ and the mean roll radius R m , has been taken into account in this model;

3) Based on the different position angle α and bite angle θ α s , the flow line element discretizes the complicated 3-dimension contact surface conveniently and makes it easier to rebuild the contact surface, and it is a good way to analyze the non-uniform stress and strain distribution accurately;

4) The discretizing and summing up method is an efficient way to solve the projected area, and results from this method approach the experimental data and simulating results very well.

Cite this paper

Dong, Y.G., Zhu, H. and Song, J.F. (2017) Novel Analytic Model for the Projected Contact Zone Based on the Flow Line Element Method in Alloyed Bar Rolling. Open Access Library Journal, 4: e3247. http://dx.doi.org/10.4236/oalib.1103247

References

  1. 1. Xu, C. (2009) The Pass Design of Structural Sections. Chemistry Industry Press, Beijing, 79. (In Chinese)

  2. 2. Hu, B. (2010) The Pass Design of Structural Sections. Metallurgical Industry Press, Beijing, 31-39. (In Chinese)

  3. 3. Yuan, Z.-X. and Wang, S.-P. (2008) The Plastic Deformation and Rolling Fundamentals. Metallurgical Industry Press, Beijing, 211-213. (In Chinese)

  4. 4. Inoue, T., Torizuka, S. and Nagai, K. (2005) Novel Rod Rolling Process for Designing Efficiently Ultrafine-Grained Steel Bars Based on Numerical Analysis. Journal of the Japan Institute of Metals, 69, 934-942. https://doi.org/10.2320/jinstmet.69.934

  5. 5. Milenin, A.A., Dyja, H. and Mroz, S. (2004) Simulation of Metal Forming during Multi-Pass Rolling of Shape Bars. Journal of Materials Processing Technology, 153, 108-114.

  6. 6. Aksenov, S.A., Chumachenko, E.N., Logashina, I.V., et al. (2015) Numerical Simulation in Roll Pass Design for Bar Rolling. Metalurgija, 54, 75-78.

  7. 7. Graf, M. and Kawalla, R. (2010) Simulation System for Fast Analysis of Multistage Hot Rolling Processes Strip and Rod/Wire. Steel Research International, 81, 122-125.

  8. 8. Nekrasov, I.I., Karamyshev, A.P., Parshin, V.S., et al. (2011) Studies of the Rolling of Ribbed Bars Using the Deform Software Package. Metallurgist, 55, 167-170. https://doi.org/10.1007/s11015-011-9408-4

  9. 9. Park, H.S. and Lee, Y. (2008) Experimental Study for Roll Gap Adjustment due to Roll Wear in Single-Stand Rolling and Multi-Stand Rolling Test. Journal of Mechanical Science and Technology, 22, 937-945.

  10. 10. Byon, S.M., Na, D.H. and Lee, Y. (2009) Effect of Roll Gap Adjustment on Exit Cross Sec-tional Shape in Groove Rolling—Experimental and FE Analysis. Journal of Materials Processing Technology, 209, 4465-4470. https://doi.org/10.1016/j.jmatprotec.2008.10.045

  11. 11. Lee, Y. (2007) Effect of Rolling Speed on the Exit Cross Sectional Shape in Rod Rolling Process. International Journal of Precision Engineering and Manufacturing, 8, 27-31.

  12. 12. Ragab, A.R. and Samy, S.N. (2006) Evaluation of Estimates of Roll Separating Force in Bar Rolling. Journal of Manufacturing Science and Engineering-Transactions of the ASME, 128, 34-45. https://doi.org/10.1115/1.2120779

  13. 13. Deng, W., Zhao, D.-W., Qin, X.-M., et al. (2010) Linear Integral Analysis of Bar Rough Rolling by Strain Rate Vector. Journal of Iron and Steel Research International, 17, 28-33. https://doi.org/10.1016/S1006-706X(10)60068-4

  14. 14. Lee, Y. and Kim, Y.H. (2001) Approximate Analysis of Roll Force in a Round-Oval-Round Pass Rolling Sequence. Journal of Material Processing Technology, 113, 124- 130. https://doi.org/10.1016/S0924-0136(01)00712-9

  15. 15. Shinokura, T. and Takai, K. (1986) Mathematical Model of Roll Force and Torque in Steel Bar Rolling. Tetsu-To-Hagane, 72, 58-64. (In Japanese)

  16. 16. Lee, Y. and Choi, S. (2000) New Approach for the Prediction of Stress Free Surface Profile of a Workpiece in Rod Rolling. ISIJ International, 40, 624. https://doi.org/10.2355/isijinternational.40.624

  17. 17. Lee, Y., Choi, S. and Kim, Y.H. (2000) Mathematical Model and Experimental Validation of Surface Profile of a Workpiece in Round-Oval-Round Pass Sequence. Journal of Materials Processing Technology, 108, 87. https://doi.org/10.1016/S0924-0136(00)00734-2

  18. 18. Lee, Y. (2001) An Analytical Study of Mean Roll Radius in Rod Rolling. ISIJ International, 41, 1414. https://doi.org/10.2355/isijinternational.41.1414

  19. 19. Dong, Y.G., Zhang, W.Z. and Song, J.F. (2008) An Analytical Model for the Prediction of Cross-Section Profile and Mean Roll Radius in Alloy Bar Rolling. Journal of University of Science and Technology Beijing, 15, 344-351. https://doi.org/10.1016/S1005-8850(08)60065-1

  20. 20. Dong, Y.G., Zhang, W.Z. and Song, J.F. (2006) A New Analytical Model for the Calculation of Mean Roll Radius in Round-Oval-Round Alloy Bar Rolling. ISIJ International, 46, 1458-1466. https://doi.org/10.2355/isijinternational.46.1458