Simulation Analysis and Optimization of Rolling Process of Steel Rim

In order to solve the problems of rolling forming accuracy and fillet thinning of alloy steel rim, a three-dimensional model of three pass rolling process was established, and the influence of different process parameters on forming quality was analyzed by using the finite element software, and the optimal process parameter combination was obtained. On this basis, the simulation results of wheel rim stress and strain for each pass rolling are analyzed, and the particle tracking technology is introduced to analyze the variation rule of stress in each incremental step. Finally, the simulation and experimental results show that the simulation thickness is basically consistent with the actual thickness, which improves the accuracy of rim rolling forming, and further verifies the correctness of rolling process simulation.


Introduction
With the acceleration of automobile lightweight process, alloy steel has been widely used in the production of automobile parts [1]. As the wheel hub is an important part of automobile driving, alloy steel for rim production can greatly reduce the mass of the wheel hub and improve vehicle driving performance [2].
However, due to the high strength and stiffness of alloy steel, the complex rolling process of wheel rim, the difficulty in controlling process parameters, high requirements on forming accuracy and thinning amount, and relatively difficult theoretical analysis. In the process of production and processing, enterprises mainly carry out experiments through "trial and error method", resulting in lation model, used the side guide wheel connecting spring to compress the blank, and optimized the simulation results to improve the engineering applicability. Wang H.L., Li Y.H., et al. [4] [5] established the finite element model of rim rolling forming, and focused on the influence of boundary conditions on the simulation results. At the same time, the metal flow and stress-strain distribution in the contact deformation zone are analyzed. Li D.Y., et al. [6] carried out simulation analysis on the roll bending process of U-beam, and optimized the roll diameter through orthogonal experiment. Tao H.W., Liu Z.W., et al. [7] [8] used numerical simulation technology to study the influence of material characteristics and roller speed on forming quality, and verified the accuracy of the analysis results through experiments. Based on the actual rolling process of the rim, the new modeling method and contact principle are determined. The influence of rolling speed and lateral force on the forming quality is considered in the parameter optimization, and the particle tracking technology is introduced to study the stress change state of the rim in different processing periods. Finally, compared with the thickness of the workpiece in the experimental production, the feasibility of the simulation is verified, which provides theoretical support for the enterprise rim manufacturing.

Deformation Mechanism of Rim Rolling Material
During the multi-step rolling process of automobile rim, the blank undergoes elastic and plastic deformation, and is formed into the final workpiece with the loading movement of the roller. The elasto-plastic deformation process is accompanied by a high degree of nonlinearity, and structural nonlinearity has always been a difficult point in the finite element analysis process. The incremental method is used to deal with the nonlinear problem, and the evolution history of stress and strain during the loading process can be obtained, and the solution accuracy is high, Good stability [9].

Nonlinear Relationship of Rim Elastic-Plastic Finite Element Method
There are many reasons for the non-linear problem of the blank, which can be divided into three main types: 1) The rim is geometrically nonlinear. The rim material undergoes a large deformation under the action of external load, which causes a large change in the force of the structure system. As the load continues to increase, the unit coordinates and structural rigidity in the workpiece change, causing the unit body to produce translation, torsion, elongation and compres- Therefore, it is necessary to numerically integrate the basic equation of infinitesimal increment form within the loading step time, and linearize the basic equation, so that the rim elastic-plastic nonlinear problem is simplified into a series of linear problems.

Elastic-Plastic Finite Element Constitutive Equation of Wheel Rim
The essence of rim elastoplastic constitutive relationship is the stress-strain relationship of the blank, the elastoplastic constitutive equation is calculated by the Von-Mises yield criterion. The yield stress is expressed by and the Von-Mises yield function is as follows The total strain increment is the sum of elastic strain increment and plastic strain increment. The calculation formula is as follows: where, ij dε is the total strain increment, The relationship between yield stress, plastic strain and work hardening coefficient is as follows: The relationship between plastic strain increment and stress state under the relevant flow criterion is as follows: The subsequent yield state can be obtained from Equation (1) as follows: Among them From Equation (5), the relationship between plastic strain increment and stress increment can be obtained as follows: Substituting formulas (4), (7) and (8) into (6), the relationship is as follows: Substituting formulas (5) and (7) into (3), the relationship is as follows: Substituting formula (10) into formula (9), the relationship is as follows: Substituting the above equation into Equation (10), the relationship between the rim elastic-plastic stress increment and strain increment can be obtained as follows:

The Rim Finite Element Solution Equation
In finite element analysis, isoparametric elements are used to discretize the solution domain, and the displacement of the element coordinates is expressed by node interpolation: where, t k i X is the coordinate component of node k in the i direction at time t, t k i u is the displacement component of node k in the i direction at time t, k N is the shape function associated with node k, and N is the number of nodes for unit interpolation.
The Piola Kirchhoff stress increment is related to the Green strain increment by the constitutive tensor, and the nonlinear equation described by T. L. can be obtained, Combining Equations (14) and (15) to calculate the displacement derivatives in the integral formulas of the full Lagrangian equation, the finite element matrix form of T. L. scheme is derived as follows [10]: are the conversion matrices of linear strain, nonlinear strain, and displacement, respectively, 0 D are the incremental material characteristic matrix, 0 t S and 0t S are the matrix formed by Piola-Kirchhoff stress.

Rolling Process
Roll forming is an advanced cold forging process, The forming process is complex, including rounding, welding and slagging, flaring, one-pass rolling, two-pass rolling, three-pass rolling, expansion finishing and Inflatable door holes and other processes [11]. However, the deformation of the rim mainly occurred in the process of three-pass rolling. In order to control the forming quality of rim, the rim three-pass rolling process of rim was simulated.

Finite Element Model of Wheel Rim
According to the rim rolling forming principle, the rim three pass rolling model  Figure 1.
The function of the first roll forming is the pre-forming of the rim, which mainly rolls out the transition shape of the deep groove and rounded corners of the rim, and the flared parts on both sides are not deformed, the second pass of roll forming is mainly to shape the groove bottom and fillet part, and pre-shape the flange part at the same time, the third pass of roll forming accurately forms the rim so that its shape and size meet the standard requirements [12]. The schematic diagrams of each pass of roll forming are shown in Figure 2.

Model Parameter Setting
In order to ensure the accuracy of simulation and obtain the accurate mechanical properties parameters of the rim material, the tensile test machine was used    Figure 3 shows the stress-strain relationship obtained during the test, and Table 1 shows the material performance parameters obtained after the tensile test.
In the rolling process, almost no deformation occurs in the mold, which is defined as analytic rigid body [13]. As the amount of deformation of billet is large, solid unit type should be adopted, hexahedral mesh should be divided, unit size is 5 mm, and unit number is 68,496.
In the actual processing, the regulating wheel is connected to the cylinder, and the blank is pressed by air pressure [14]. In the simulation, the effect of air pressure is replaced by the die spring, so that the regulating wheel and the flange are inflexible contacts, and the boundary conditions similar to the actual effect are achieved, as shown in Figure 4.
The rolling time includes feed forming time and progressive forming time, which are related to feed quantity, feed speed and blank rotation speed. The time step is calculated by where, t is the time step, t 1 is the feed forming time, t 2 is the progressive forming time, v is the feed velocity, d is the feed quantity, n is the number of turns of the progressive forming workpiece, n is 1.
The billet rotation speed (w) is calculated based on the equal linear velocity of the rolling contact point.
where, v 1 is the linear velocity of the contact point, r 1 is the roll-down radius, w 1 is the roll-down speed, r 2 is the blank radius, and w 2 is the blank speed.

Analysis of the Influence of Process Parameters on Forming Quality
The main process parameters that affect the quality of rim roll forming include friction coefficient, feed speed and wheel speed. In actual production, the process parameters are mainly selected by the researcher's experience and cannot ensure the optimal combination of parameters. In order to study the relationship between process parameters and forming quality, a single parameter experiment method was used to simulate the rim.

1) The influence of friction on the forming quality
In rolling processing, the workpiece and the mold are rotated by friction, but due to the influence of lubrication conditions and the wear of the mold surface, it is difficult to determine the relationship between the friction and the work- 2) The influence of feed speed on forming quality In actual rolling, if the feed speed is too high, it will cause the blank to run off and the forming will be unstable, if the feed speed is too small, the rolling time will be prolonged and the processing efficiency will be reduced.   in actual processing, the feed speed of the mold is too fast, which will cause the workpiece to produce lateral force and the workpiece will shift. As shown in Figure 6(b) the lateral force curve at different feed speeds, as the feed speed becomes larger, the lateral force also becomes larger, which is more likely to cause the workpiece to shift. Comprehensively considering its influence on the forming thickness and lateral force, multiple simulations and adjustments, a feed speed of 21.89 mm/s can obtain an ideal workpiece quality.
3) The influence of roller speed on forming quality The speed ratio of the roller is matched inversely according to the diameter of the roller [15], which ensures that the linear velocity of the roller and the workpiece are equal on the contact surface, and reduces the forming defects of the workpiece. With the friction coefficient of 0.3 and the feed speed of 21.89mm /s unchanged, the influence of the rotating speed of the upper roller die on the forming quality of the workpiece at 150 r/min, 200 r/min and 250 r/min was studied. Figure 7(a) shows the thickness change curve of measuring points at different rotating speeds. It is found in the figure that with the increase of rotating speed, the thickness of groove, flange and rounded corners in a roll decreases to different degrees, and the thinning of rounded corners is the most serious. In order to ensure the safety of the wheels, the speed of the rollers should be appropriately reduced to control the thinning of the rim. However, if the speed is too small, it will cause uneven forming of the rim, as shown in Figure 7(b) for the rim strain curve under the same section. Comparative analysis found that the deformation shapes of the two sections of the rim were obviously different when the rotating speed was 150 r/min. In particular, the upper section had small strain and deformation, which could not meet the forming requirements. However, the workpiece with uniform thickness and better quality could be obtained when the rotating speed was about 200 r/min.

Analysis of Simulation Results under Parameter Optimization Conditions
Simufact Forming software is used to perform finite element simulation on the rim. In the simulation, the feed distance is the translational distance of the roller, and the analysis results of the speed, feed speed, time step and feed distance are shown in Table 2.
The analysis cloud diagram of the rolling stress and strain of each pass is shown in Figure 8.    ing that the plastic deformation occurs in the groove area is larger.
The second pass roll forming shape, mainly preform the flange area, and further shape the groove and rounded corners. Figure 9 shows the equivalent stress cloud diagram and equivalent strain cloud diagram of the second rolling. The equivalent stress at the end of the blank rolling reaches 714.84 MPa, and the maximum equivalent stress appears in the flange part, and the groove and rounded also produce large stress; The maximum equivalent strain is 0.47, which also occurs at the flange, indicating that the material deformation at this place is also large.
The third pass rolls mainly shape the wheel flange accurately, and makes the groove fillet size meet the requirements of the workpiece. Figure 10 shows the equivalent stress cloud diagram and equivalent strain cloud diagram of the third pass rolls. The stress value at the end of rolling reaches 726.69 MPa, the strain value reaches 0.64, and the large stress and large strain are mainly concentrated at the flange and groove fillet. Since the basic shape of the rim has been formed after the first and second rolls, the third pass rolls are mainly refined, with small deformation and more even stress distribution.
In order to analyze in detail the rolling simulation process of each pass of the rim, the particle tracking technology is introduced into the rolling results to

Roll Test Verification
In order to verify the correctness of the simulated rim forming, the actual production workpiece and the simulated formed workpiece were compared and analyzed in thickness. The rim rolling process is shown in Figure 14. According to the actual rolling equipment 14(a), the rim rolling forming device 14(b) is simplified, which is mainly composed of up roll, down roll, blank, guide wheel, air cylinder and main shaft [16].   Through the rim rolling process, each pass roll forming workpiece is obtained. Figure 15 shows the cross-sectional shape of the actual roll forming rim and the simulated cross-sectional shape. The actual cross-sectional shape of the workpiece is obtained by computer vision measurement. The thickness of the measuring point is mainly selected from the areas that are more sensitive to thickness changes during the rolling process (rounded corners, grooves and flange). Figure 16 shows the comparison curve between actual measured thickness and simulated thickness.
It can be found from Figure 15 that the simulated cross-sectional shape is basically the same as the actual processed cross-sectional shape. The simulation data curve of each pass in Figure 16. is basically consistent with the actual measurement curve, which proves the feasibility of simulation. Among them, the largest thinning area of the roll thickness occurred at the rounded corner, which was subjected to large tensile and compressive stresses, and the deformation was  large, and the thinning rate reached 10%, the flange is subjected to little compressive stress, and under the action of tensile stress, the material flows to the fillet to compensate for excessive thinning of fillet; The largest thickness thinning area of the second roll occurs at the wheel rim, mainly because of the preforming of the wheel rim, the stress is large, the deformation is large, and there is no material supplement in this area, the thinning rate reaches 10.8%, and the thinning rate of other positions is small; After three-roll forming, the thickness distribution is relatively uniform, and the largest thinning area occurs at the rounded corner, and the thinning rate is about 9.7%.

Conclusions
1) The increase of the friction coefficient is beneficial to the forming of the workpiece, but after the friction coefficient is increased to 0.3, the effect on the forming quality is not obvious. The increase of the feed speed will increase the thickness of the measuring point, which is conducive to forming, but if the feed speed is too high, the lateral force will increase, which will easily cause the workpiece to shift. The reduction of the speed is beneficial to control the thinning of the rim, but if the speed is too small, the formation of the rim will be uneven.
2) The simulation results show that first rolling rounded corner and groove produce greater stress and strain. Second rolling pre-formed flange, where the stress and strain are both large, and the maximum equivalent stress appears at the flange, and the groove and rounded corner also produce large stress. Third rolling rim is accurately formed, the stress distribution is more uniform, and the fillet deformation is relatively small.
3) Comparative analysis of experimental and simulation results. The simulated forming thickness is basically consistent with the experimentally measured thickness, which proves the feasibility of simulation. The design optimization of new products can be carried out on this basis, which has certain reference value for the production and processing of enterprises.

Authors Contributions
Conceptualization