Finite Element Analysis and Experimental Evaluation of Residual Stress of Zr-4 alloys Processed through Swaging

: Zirconium alloy has been extensively used as a cladding material in nuclear power reactors due to its low neutron absorption cross section, excellent mechanical properties, and corrosion resistance. The inﬂuence of the swaging parameter, feed rate (0.7, 1.25, 2 m / min) on residual stress induced in Zr-4 alloy is investigated in the present work. A three-dimensional ﬁnite element model was implemented in the Deform 3D software to simulate the rotary swaging (RS) process over a circular rod of Zr-4 alloy. The simulation results based on the 3D framework provide a detailed insight of residual stress, true stress versus true strain and force applied over the rod during the multiple pass swaging process; the results are compared with experimental results. The experimental hole drilling method is used to determine the residual stresses on swaged zirconium alloy at di ﬀ erent feed rates (0.7, 1.25, and 2 m / min). A similar trend of residual stress between experimental and numerical results from the surface to the center on the swaged rod samples is observed. The same magnitude of residual stress at the surface of the swaged Zr-4 rod is also observed. It is found to be compressive at the surface and tensile in the center of the samples, as observed in the present work. numerical methods show similar residual stress magnitudes quantitatively at the surface of the swaged rod. For the 0.7 m / min feed rate, deep hole drilling and numerical methods have shown 6 MPa and − 0.9 MP, respectively. In the case of 1.25 m / min feed rate, deep hole drilling technique showed − 52 MPa and − 41.4 MPa from the numerical method. For the 2 m / min feed rate, deep hole drilling showed − 86 MPa and − 74 MPa with the numerical method. Both numerical and experimental methods predicted a similar trend in variation of tensile residual stresses but they di ﬀ er slightly in their magnitude. These results are for swaged samples at di ﬀ erent feed rates and it suggests that increasing the feed rate increases the compressive residual stresses on the surface while tensile stresses are nearly similar on increasing the feed rates. With these


Introduction
Zirconium alloys processed through the swaging technique are used in nuclear power reactors and chemical industries due to their excellent tensile properties, fracture toughness and corrosion resistance. Swaging is a metal forming process by which the diameter of a rod or tube is reduced by forcing it into a die with the help of a reciprocating blow. The forming dies of the swaging machine are arranged concentrically around the work piece ( Figure 1a). The swaging dies perform high frequency radial movements with short strokes. Rotary swaging (RS) can be performed under hot and cold conditions, depending on the swaged material, required surface quality and final properties [1,2]. Infeed rotary swaging, the work piece is fed axially inside the set of dies but in the recess swaging process, the work piece is not fed axially forward or outward (Figure 1a). The sinking zone is much smaller than the forging zone. Therefore, the infeed swaging process is often utilized in the situation where the work piece is necked from one end. The recess swaging process is often used to reduce the diameter at a certain position of the work piece, or to form a concave profile [1][2][3][4][5][6].  [2] and (b) the zones of the swaging process [3], with permission from Elsevier, 2020.  [2] and (b) the zones of the swaging process [3], with permission from Elsevier, 2020.
The work piece can be classified into three zones when it is inside the dies; they are the sinking zone, the forging zone and the sizing zone ( Figure 1b) [3]. Similarly, the dies can be classified into three zones such as the reduction, calibration and exit zones. The outside part at the sinking zone suffers from triaxial compressive stress as it is in direct contact with the die and therefore, the deformation of this part is more severe. The outside part at the sizing zone suffers from axial tensile stress due to the axial elongation. The axial elongation of the work piece as well as the necking predicted in the simulation are caused by the axial tensile stress by virtue of the deformation occurring in the sizing zone [1]. During the rotary swaging process, the forming of the work piece takes place in the swaging head in small steps by a radial oscillating movement of the tools. A rotary swaging process utilizing an energy-controlled method was discussed in [1]. The deformation of the tube mainly occurs at the sinking zone. The radial action of the dies is generated by rotation of the driven swaging spindle due to the cam shaped outline of the shoe. Through the roll-off-process between the shoes and the cylinder rollers, the dies are moved inwards simultaneously at every overrun of the cams. The cross section of the axially fed work piece becomes reduced. Due to the rotation of the swaging spindle relative to the work piece, a constant circumferential shaping can be obtained [7,8].
Various parameters of swaging can affect the microstructure, mechanical properties and residual stresses on the materials being processed. The die angle, die stroke, feed velocity, and the friction coefficient would affect the deformation of the work piece. Incremental forming processes have been studied in recent years, in most cases using finite element modelling (FEM) due to its high precision. Specifically, FEM is used in rotary swaging in order to study the strain and stress distributions in the work piece [4], especially the influence of the axial feed velocity on the strain field [5], the influence of friction and the die angle on residual stresses [6] and temperature distribution during the process [9,10]. Axisymmetric quadrilateral full integration elements were chosen to develop the finite element analysis (FEA) model of rotary swaging for an Mg alloy [11]. FEA was conducted for the tungsten heavy alloy using FORGE software; the highest strain was observed in the surface regions, which were directly affected by the swaging dies [12]. Using Forge NxT commercial software, FEA was conducted for an elastic-viscoplastic material model [12].
Swaging has led to a substantial increase in the ultimate tensile strength and a decrease in ductility [13,14]. The finite element (FE) model has shown the variation of strain from surface to central axis of swaged Al-Cu clad composite [15]. The residual stress component in the axial direction is much larger than that in other directions and as a result, it becomes one of the basic factors causing transverse surface cracks in the subsequent rotary swaging (RS) passes [16]. The zones of swaging are explicitly shown using ABAQUS software (6.16, Dassault Systèmes ® , Vélizy-Villacoublay, France). The axial reaction was compared using lubricant as well as without lubricant [17]. A combined friction law was used to model the friction in the die-workpiece interface for recess swaging [18]. Forge 2011 3D was used for FEA of copper AISI C11000 tube that suffers triaxial stress during the rotary swaging process [3]. Fluctuations in the stroke height, which are in the range of the tolerance of the components of the swaging head, has led to strong fluctuations in residual stress distribution [2]. The 2D axisymmetric and 3D FE analysis is used in modelling residual stresses and forging load in the cold radial forging process. Metal flow in the deformation gap has been performed using the finite element method in the real swaging process [19]. The Lagrange type incompressibility condition and boundary conditions for velocity were considered as reported in the literature [20]. The local strain distribution calculations were verified using experimentally available data on Inconel 718 [21]. FEM analysis based on plane 2D elements was used for determining the state of stress and analysis of the metal flow in swaging [8].
Zirconium alloy is used in nuclear fission reactor applications due to its superior mechanical properties compared to other structural material like stainless steel [22]. Zirconium alloys have a good combination of strength and ductility, high corrosion and oxidation resistance and a low neutron absorption cross section area [23][24][25][26]. Residual stresses are induced into the materials due to plastic deformation, annealing, phase transformation and welding in material components [27]. Residual stress is an important criterion for choosing the cladding material for nuclear applications because of severe plastic deformation processes such as rotary swaging, extrusion, the rolling process, and hot pressure torsion and dissimilar metal welding used for processing this alloy [28][29][30]. Residual stresses are measured by different methods such as, X-Ray diffraction (XRD), hole drilling and deep hole drilling, and the hole contour method and block removal [31]. XRD has the ability to measure only surface residual stresses, while semi destructive techniques like deep hole drilling are able to measure the residual stress up to 750 mm depth in the material [32]. Various studies have been carried out on other alloys such as titanium alloys, tungsten-based alloys and copper alloys. The effect of rotary swaging and annealing on residual stress for WNiCo powder-based pseudo alloy was predicted using numerical and experimental techniques (XRD technique) [33]. It was reported that swaging at room temperature introduced higher compressive stresses which decreased after the post heat treatment process [13,33]. Residual stress analyses through X-ray diffraction and the elastoplastic self-consistent (EPSC) model were made on cold pilgered zirconium alloy [34]. The contribution and magnitude of the first and second order residual stresses were evaluated using inputs from the modeling [34]. It is possible to determine the basic characteristics of strain distribution, the strain rate, the local strain intensities, and the displacement rates of the work piece, as well as the swaging force. To understand residual stresses, Henriksen et.al conducted a series of tests and stated that grain deformation in the surface layer generates residual stresses [35]. Therefore, according to his study, mechanical deformation is the key factor for residual stress generation in the materials. A parametric study of machining parameters in order to understand the impact on the residual stresses in materials is reported using an arbitrary Lagrangian-Eulerian (ALE) finite element approach [36].

Finite Element Analysis
FE analysis is a consequence of the trend followed by industry in the past few decades towards an improvement of the process and the use of total quality programs. There is an increasing demand in recent times to produce more output with better quality and lower costs. Consequently, the use of computer modelling has become increasingly appealing to industry. Its application provides important Metals 2020, 10, 1281 4 of 15 design information on the shape and integrity of the part, the required force, wear on tools and an understanding of the metal working process, where conventional analytical methods are restricted to the steady-state stage of the process. Thus, conventional methods cannot provide an insight into the dynamic changes that occur during the initial stage of the process. In addition, the simulations are advantageous for their ability to easily vary geometry and boundary conditions, consequently facilitating detailed studies of deformation patterns. The process (billet, tools), once modelled, can be used for a large variety of investigations by simply changing the boundary conditions in the model.
The literatures on determination of residual stress in zirconium alloy by numerical and experimental methods is limited. A 3D finite element model based on plastic as well as elastic-plastic formulation are developed using the commercial finite element software DEFORM 3D (12.1, Scientific Forming Technologies Corporation Columbus, Ohio, Unites States) in the present work to analyse the residual stress in zircaloy-4 (Zr-4 alloy). DEFORM 3D uses the finite element method to solve large deformation problems and an automatic mesh generator (AMG) to automatically provide an optimized remeshing capability.
The implicit FE method was used for the analysis, where equilibrium global equations were simultaneously solved implicitly with the Newton-Raphson (N-R) method. The iterative approach used in the implicit FE method may have issues achieving convergence for non-linear material [37]. Explicit FE results are compared with the results of an implicit formulation and computing times were benchmarked for different problems [38]. Implicit FEA in 3D is computationally expensive, but recent works have shown that it can scale well on large high-performance computing machines or clouds [39].

Modelling
Considering that rotary swaging is a cold forming process, the following assumptions were considered to analyze the processing of Zr-4 alloys: (a) the material is uniform and isotropic, (b) the interface of the dies and the work piece are under isothermal conditions, and (c) the gravity and inertial force are ignored. The friction condition in the die-work piece interface was modelled by the combined (hybrid) friction law, which is the combination of two friction models. The friction model is assumed to be sliding velocity dependent, with a friction coefficient of 0.4, which is commonly used for cold forging conditions. The swaging dies moved radially toward the rod and deformed it. The initial length and diameter of the rod were 500 and 21 mm, respectively. Table 1 shows some process parameters used for the FE modeling of Zircaloy-4. The Lagrangian incremental model is generally used for all the conventional applications like forming, heat transfer or heat-treatment. It is also useful for modelling the transient phase of processes like rolling, machining, extrusion, drawing, cogging etc. The main advantage is due to easier remeshing, which is useful for large scale deformation as well. However, in case of large deformations, the mesh quality quickly degrades, leading to inaccuracy or even instability, and remeshing is inevitable. Mesh generation for such remeshing can not only be time and resource consuming, but there would also be an error associated with the projection of the solution onto the new mesh; it is the major drawback of the arbitrary Lagrangian-Eulerian (ALE) method.
The convergence of an elasto-plastic solution is dependent on the initial guess of the stress-strain state. There are three initial guess solutions available: (1) the elastic solution, which uses purely elastic deformation data, (2) the plastic solution, which uses plastic deformation data to generate the initial guess, and (3) the previous step solution, which uses the elasto-plastic solution from the previous step to generate the initial guess. In most of the cases, the latter, previous step solution seems to give the best convergence.
The Newton-Raphson (N-R) method is used to solve the problem. It is recommended because generally it converges in fewer iteration than the other available methods [37]. However, solutions are more likely to fail to converge with this method than with other methods. Whereas, the direct iteration method is more likely to converge than the N-R method, but to do so, will generally require more iterations. In the Deform software, matrix equations take care of solving especially large-scale deformation problems.

Material
In order for a simulation to achieve a high level of accuracy, it is important to have an understanding of the material properties required to specify a material used in DEFORM 3D. The Zr-4 alloy, with an equivalent stress-strain relationship at room temperature, was used for the analysis. The Create a New Material option is used to define the Zr-4 alloy. The rod was assumed to be elasto-plastic. For elastic data, elastic properties like Young's modulus, Poisson's ratio and the thermal expansion coefficient were used for the deformation analysis. Assuming that there was no internal heat generation, a Young's modulus of 99 GPa and a Poisson's ratio of 0.37 were considered for the simulation. For plastic material, the tabular data format was used. This format is the most versatile format as it can represent any material where flow stress can be given as a function of strain, strain rate, and temperature. This flow stress curve is used to determine the residual stress analysis from the simulation. The Von Mises stress criterion is imposed during simulation. The die and manipulator are considered as rigid bodies for this analysis.

Mesh
The geometry of the elasto-plastic work specimen is meshed with finite element meshes of tetrahedron elements (shown in Figure 2a); the mesh of the work piece (initial diameter of 21 mm and length of 500 mm) consisted of 27,806 elements and 7142 nodes. For plastic material, brick mesh elements with 6500 elements and 7878 nodes are generated. Deform 3D provides the facilities for automatic remesh wherever necessary. Remeshing criteria were selected as it is the most convenient way to handle the remeshing of objects undergoing large plastic deformation. Remeshing is the process of replacing a distorted mesh with a new undistorted mesh and interpolating the field variables (strain, velocity, damage, and temperature etc.) from the old mesh to the new mesh.

Boundary Condition
Contact boundary conditions display inter-object boundary contact conditions on a given object. Contact boundary conditions are applied to nodes of a slave object (work piece), and specified contact between those nodes and the surface of a master object (die and manipulator). In master-slave combination, in case of a single deforming object, the deforming object should always be the slave object. In case of multiple deforming bodies, the object with the finer mesh at the interface of the two objects should be the slave object. Contact must be specified for every pair of deformable objects that may contact each other during the simulation. At the interface between two objects, the constant

Boundary Condition
Contact boundary conditions display inter-object boundary contact conditions on a given object. Contact boundary conditions are applied to nodes of a slave object (work piece), and specified contact between those nodes and the surface of a master object (die and manipulator). In master-slave combination, in case of a single deforming object, the deforming object should always be the slave object. In case of multiple deforming bodies, the object with the finer mesh at the interface of the two objects should be the slave object. Contact must be specified for every pair of deformable objects that may contact each other during the simulation. At the interface between two objects, the constant friction coefficient is specified. A separation criteria was imposed for inter object relationships. The dies are assigned a sliding velocity of 30 mm/s over the work piece.

Residual Stress Analysis using FEM
Prediction and control of residual stresses and stress concentration in the materials is essential for understanding the fracture toughness and fatigue life of the components. Residual stresses in the work piece are periodically introduced during in various thermo-mechanical forming processes. To predict residual stress values from the simulations of the swaged Zr-4 alloy, 20 measuring points were chosen along the diameter, using the point-tracking technique. The dies are taken away from the tube. These steps provide significant time to recover the elastic part of the formulation. The stress remaining in the tube after the relaxation period was assumed to be the residual stress.

Experimental-Residual Stress Measurement
The through thickness residual stresses are measured by using the deep hole drilling technique. The main advantage of using this technique is due to the fact that it can measure the residual stress fields through the thickness of the swaged rod samples. The tests are carried out in four steps: 1. A through reference hole is drilled at the measured location using gun drilling. 2. The measurement of the reference hole diameter is made at pre-decided angles for each increment in the depth direction. 3. A core of material containing the reference hole is trepanned free of the rest of the material using an electric discharge machine (EDM) (Makino, Mason, Ohio, United States). 4. The reference hole diameter is again measured to calculate the magnitude of residual stress fields from the amount of deformation observed in the inner core column of trepanning process.

Residual Stress Analysis Using FEM
Prediction and control of residual stresses and stress concentration in the materials is essential for understanding the fracture toughness and fatigue life of the components. Residual stresses in the work piece are periodically introduced during in various thermo-mechanical forming processes. To predict residual stress values from the simulations of the swaged Zr-4 alloy, 20 measuring points were chosen along the diameter, using the point-tracking technique. The dies are taken away from the tube. These steps provide significant time to recover the elastic part of the formulation. The stress remaining in the tube after the relaxation period was assumed to be the residual stress.

Experimental-Residual Stress Measurement
The through thickness residual stresses are measured by using the deep hole drilling technique. The main advantage of using this technique is due to the fact that it can measure the residual stress fields through the thickness of the swaged rod samples. The tests are carried out in four steps:

1.
A through reference hole is drilled at the measured location using gun drilling.

2.
The measurement of the reference hole diameter is made at pre-decided angles for each increment in the depth direction.

3.
A core of material containing the reference hole is trepanned free of the rest of the material using an electric discharge machine (EDM) (Makino, Mason, Ohio, United States). 4.
The reference hole diameter is again measured to calculate the magnitude of residual stress fields from the amount of deformation observed in the inner core column of trepanning process.
The following Equations (1) and (2) were used to determine the residual stresses (σ y , σ z ): where, ε y and ε z are measured strains in the y and z directions, respectively; ∆d y and ∆z are the change in diameters in the y and z directions, respectively; d y and d z are the initial diameters, in the y and z directions, respectively; E is Young's modulus; and σ y and σ z are the stresses in the y and z directions, respectively [31,40]. The through thickness residual stress measurement is explained in detail in our earlier work [40]. The Zr-4 alloys are produced by the induction melting process and swaged up to a 16.2 mm rod from 21 mm rod samples in 5 passes (NFC Hyderabad, India). The composition of Zr-4 alloy is shown in Table 2. The through thickness residual stresses measured by the deep hole drilling technique of 25% (18.2 mm rod) swaged samples (0.7 m/min, 1.25 m/min and 2 m/min) and 40% (16.2 mm rod) swaged samples (0.7 m/min). The residual stresses along both the directions (axial and radial directions) peak near the middle point in the thickness direction and accordingly the stress patterns are represented for the swaged bars, as shown in Figure 2b. Samples are polished by emery papers because surface defects may affect the strain relaxation readings. Young's modulus and Poisson's ratio of Zr-4 alloy are 99 GPa and 0.37, respectively [23,24]. The tensile properties (yield strength and tensile strength) of Zr-4 alloys are taken from our earlier work [23]. The scanning electron microscopy-based electron backscattered diffraction (EBSD) (FEI, Hillsboro, OH, USA) was carried out on 25% of the samples at different feed rates for kernel average misorientation (KAM) mapping. The samples were polished with up to 2000 µm emery papers, followed by alumina and colloidal cloth polishing. Electropolishing was done using 80% methanol and 20% perchloric acid for 17 s at −30 • C.

FEM Results
A rod with initial diameter of 21 mm and length of 500 mm was considered for plastic deformation, with single pass to reduce the diameter of the rod by 1 mm. The final length is calculated as 513.21 mm. The true σ versus the true ε curve is plotted using origin software for the swaged piece and is shown in Figure 3a. and z directions, respectively [31,40]. The through thickness residual stress measurement is explained in detail in our earlier work [40].
The Zr-4 alloys are produced by the induction melting process and swaged up to a 16.2 mm rod from 21 mm rod samples in 5 passes (NFC Hyderabad, India). The composition of Zr-4 alloy is shown in Table 2. The through thickness residual stresses measured by the deep hole drilling technique of 25% (18.2 mm rod) swaged samples (0.7 m/min, 1.25 m/min and 2 m/min) and 40% (16.2 mm rod) swaged samples (0.7 m/min). The residual stresses along both the directions (axial and radial directions) peak near the middle point in the thickness direction and accordingly the stress patterns are represented for the swaged bars, as shown in Figure 2b. Samples are polished by emery papers because surface defects may affect the strain relaxation readings. Young's modulus and Poisson's ratio of Zr-4 alloy are 99 GPa and 0.37, respectively [23,24]. The tensile properties (yield strength and tensile strength) of Zr-4 alloys are taken from our earlier work [23]. The scanning electron microscopy-based electron backscattered diffraction (EBSD) (FEI, Hillsboro, OH, USA) was carried out on 25% of the samples at different feed rates for kernel average misorientation (KAM) mapping. The samples were polished with up to 2000 µm emery papers, followed by alumina and colloidal cloth polishing. Electropolishing was done using 80% methanol and 20% perchloric acid for 17 s at −30 °C.

FEM Results
A rod with initial diameter of 21 mm and length of 500 mm was considered for plastic deformation, with single pass to reduce the diameter of the rod by 1 mm. The final length is calculated as 513.21 mm. The true σ versus the true ε curve is plotted using origin software for the swaged piece and is shown in Figure 3a. During the swaging process, due to the tolerance applied to the dies (3 or 4), the swaged work piece is not uniform, as some flash was produced on the outer surface of the rod as shown in Figure   Figure 3. (a) The true σ versus the true ε of the work piece, obtained from numerical data, and (b) the flash produced over the outer surface of the work piece During the swaging process, due to the tolerance applied to the dies (3 or 4), the swaged work piece is not uniform, as some flash was produced on the outer surface of the rod as shown in Figure 3b. To avoid this problem the work piece is allowed to rotate after each bite/blow. With this rotation, it is possible to have a uniform surface of the swaged rod. The rotation angle per bite of 15 degrees was assumed, which does not lead to a uniform surface. Later, the rotation of 90 degrees was used to ensure a better surface after the operation. Now, considering an elasto-plastic case, the rod of 500 mm length was deformed for three consecutive passes; in the first pass the diameter of the work piece was reduced from 21 mm to 20 mm and in the next two passes, it was reduced to a final diameter of 18.2 mm. The final length obtained after two subsequent passes is 540.671 mm, as shown in Figure 4a. The stress distribution contour over a particular cross section of the rod is displayed after slicing the rod as shown in Figure 4b. Now, considering an elasto-plastic case, the rod of 500 mm length was deformed for three consecutive passes; in the first pass the diameter of the work piece was reduced from 21 mm to 20 mm and in the next two passes, it was reduced to a final diameter of 18.2 mm. The final length obtained after two subsequent passes is 540.671 mm, as shown in Figure 4a. The stress distribution contour over a particular cross section of the rod is displayed after slicing the rod as shown in Figure  4b. A path operation was performed over the swaged work piece for a random cross section, as shown in Figure 5, using the point tracking option of the software to see the σ and ϵ variation across the cross section. Stress/strain generated is the maximum at the surface and gradually decreases towards the center of the work piece. Hence, the surface is the highly stressed region, as can be seen from Figure 4b. From this analysis, the load applied by the die on the work piece has been generated. Figure 6a indicates the Z-load and Figure 6b indicates the Y-load. The load at the top and bottom are almost A path operation was performed over the swaged work piece for a random cross section, as shown in Figure 5, using the point tracking option of the software to see the σ and variation across the cross section. Stress/strain generated is the maximum at the surface and gradually decreases towards the center of the work piece. Hence, the surface is the highly stressed region, as can be seen from Figure 4b. mm and in the next two passes, it was reduced to a final diameter of 18.2 mm. The final length obtained after two subsequent passes is 540.671 mm, as shown in Figure 4a. The stress distribution contour over a particular cross section of the rod is displayed after slicing the rod as shown in Figure  4b. A path operation was performed over the swaged work piece for a random cross section, as shown in Figure 5, using the point tracking option of the software to see the σ and ϵ variation across the cross section. Stress/strain generated is the maximum at the surface and gradually decreases towards the center of the work piece. Hence, the surface is the highly stressed region, as can be seen from Figure 4b. From this analysis, the load applied by the die on the work piece has been generated. Figure 6a indicates the Z-load and Figure 6b indicates the Y-load. The load at the top and bottom are almost From this analysis, the load applied by the die on the work piece has been generated. Figure 6a indicates the Z-load and Figure 6b indicates the Y-load. The load at the top and bottom are almost the same. The maximum value is located where the load is exerted by the right and the minimum value is exerted by the left; as measured along the Z direction.
In Figure 6b, the load exerted by the right die and the left die, which are in the horizontal direction, are the maximum and the top/bottom die exerts minimum load on the work piece as it measures the load in the Y direction. At the same time, all four dies exert the same amount of load over the work piece, to have uniform deformation. The effective stress with respect to time is plotted from elasto-plastic data undergoing multi-pass swaging. From Figure 7, it is clear that the value of stress increases with the number of passes. Actually, the results presented here refer to the working stress during the swaging process and the results give the value of the residual stress after successful completion of the swaging process. From the load curve, it can be observed that at the same time all four dies exert equal amounts of load over the work piece; this leads to a uniform homogeneous deformation in the axial direction of the swaged Zr-4 rod.
Metals 2020, 10, x FOR PEER REVIEW 9 of 16 the same. The maximum value is located where the load is exerted by the right and the minimum value is exerted by the left; as measured along the Z direction. In Figure 6b, the load exerted by the right die and the left die, which are in the horizontal direction, are the maximum and the top/bottom die exerts minimum load on the work piece as it measures the load in the Y direction. At the same time, all four dies exert the same amount of load over the work piece, to have uniform deformation. The effective stress with respect to time is plotted from elasto-plastic data undergoing multi-pass swaging. From Figure 7, it is clear that the value of stress increases with the number of passes. Actually, the results presented here refer to the working stress during the swaging process and the results give the value of the residual stress after successful completion of the swaging process. From the load curve, it can be observed that at the same time all four dies exert equal amounts of load over the work piece; this leads to a uniform homogeneous deformation in the axial direction of the swaged Zr-4 rod. Moreover, for the above cases, the die has been designed for the optimum value of the required diameter of the rod. However, as our analysis is limited to a very small value of diameter reduction, the appropriate consideration of the die geometry would play a significant role. Here, the dies have  In Figure 6b, the load exerted by the right die and the left die, which are in the horizontal direction, are the maximum and the top/bottom die exerts minimum load on the work piece as it measures the load in the Y direction. At the same time, all four dies exert the same amount of load over the work piece, to have uniform deformation. The effective stress with respect to time is plotted from elasto-plastic data undergoing multi-pass swaging. From Figure 7, it is clear that the value of stress increases with the number of passes. Actually, the results presented here refer to the working stress during the swaging process and the results give the value of the residual stress after successful completion of the swaging process. From the load curve, it can be observed that at the same time all four dies exert equal amounts of load over the work piece; this leads to a uniform homogeneous deformation in the axial direction of the swaged Zr-4 rod. Moreover, for the above cases, the die has been designed for the optimum value of the required diameter of the rod. However, as our analysis is limited to a very small value of diameter reduction, the appropriate consideration of the die geometry would play a significant role. Here, the dies have Moreover, for the above cases, the die has been designed for the optimum value of the required diameter of the rod. However, as our analysis is limited to a very small value of diameter reduction, the appropriate consideration of the die geometry would play a significant role. Here, the dies have been designed according to the required diameter (e.g., 21-20, 20-19, etc.), i.e., different die geometry is used for each 1 mm reduction of diameter.
As a post-processing result of this analysis, a 20 mm diameter rod is obtained with an incremental length. Now, instead of a multi pass operation, the obtained 20 mm diameter rod is imported to the next set of simulations using the respective die. Assuming a similar condition of all other parameters involved, the swaging process is completed for the next 1 mm diameter reduction. Repeating this procedure for subsequent steps (e.g., 19 to 18.2 mm and 18.2 to 17.4 mm), a final swaged rod of diameter nearly equal to 16.2 mm is obtained (Figure 8a). True σ versus the curve is obtained from experimental as well as numerical analysis and they are plotted together and shown in Figure 8b.
other parameters involved, the swaging process is completed for the next 1 mm diameter reduction. Repeating this procedure for subsequent steps (e.g., 19 to 18.2 mm and 18.2 to 17.4 mm), a final swaged rod of diameter nearly equal to 16.2 mm is obtained (Figure 8a). True σ versus the ϵ curve is obtained from experimental as well as numerical analysis and they are plotted together and shown in Figure 8b. Our analysis is based on a varying axial feed rate. The variation of effective stress was estimated for different feed rates (0.7 m/min, 1.25 m/min and 2 m/min). Using Deform 3D, the time step is specified as 0.5 s for a bite and accordingly, the axial feed per bite has been calculated and is given as input in the pre-processor. According to the parameters, die velocity and stroke in the radial direction is fixed. The obtained results are shown below.
Residual stresses in the work material are periodically distributed during the forming process. To ensure a stress value from the simulations that describe a homogeneous residual stress state, 20 measuring points were chosen along the diameter; the positions are seen in Figure 5. The residual stresses were extracted after the die was released from the work piece, with subsequent cooling to room temperature. It is observed that with an increasing number of passes, the value of residual stress also increases. This can be observed from the FEM analysis for the different number of passes. The residual stresses were measured along the radial direction of the rod at different cross sections of the swaged work piece.
Numerical simulation was performed to predict the residual stress over the surface as well as over the cross-section of the 25% swaged samples, with different feed rates using DEFORM 3D software. This result (Figure 9a-c) shows the compressive residual stress at the surface and the tensile residual stress at the center in the radial direction of the rod. The obtained results qualitatively match with the experiment data. It may be because of varying parameters involved in the complete swaging process, such as die geometry, die movement velocity, and swaging stroke. In the next step, die geometry was varied according to the required diameter, i.e., 5 nos. of die sets (one die set contains 4 segments) for swaging the bar from ϕ 21 to ϕ 16.  Our analysis is based on a varying axial feed rate. The variation of effective stress was estimated for different feed rates (0.7 m/min, 1.25 m/min and 2 m/min). Using Deform 3D, the time step is specified as 0.5 s for a bite and accordingly, the axial feed per bite has been calculated and is given as input in the pre-processor. According to the parameters, die velocity and stroke in the radial direction is fixed. The obtained results are shown below.
Residual stresses in the work material are periodically distributed during the forming process. To ensure a stress value from the simulations that describe a homogeneous residual stress state, 20 measuring points were chosen along the diameter; the positions are seen in Figure 5. The residual stresses were extracted after the die was released from the work piece, with subsequent cooling to room temperature. It is observed that with an increasing number of passes, the value of residual stress also increases. This can be observed from the FEM analysis for the different number of passes. The residual stresses were measured along the radial direction of the rod at different cross sections of the swaged work piece.
Numerical simulation was performed to predict the residual stress over the surface as well as over the cross-section of the 25% swaged samples, with different feed rates using DEFORM 3D software. This result (Figure 9a-c) shows the compressive residual stress at the surface and the tensile residual stress at the center in the radial direction of the rod. The obtained results qualitatively match with the experiment data. It may be because of varying parameters involved in the complete swaging process, such as die geometry, die movement velocity, and swaging stroke. In the next step, die geometry was varied according to the required diameter, i.e., 5 nos. of die sets (one die set contains 4 segments) for swaging the bar from φ 21 to φ 16. Parameters play a significant role in the variation of residual stresses. A better value of the obtained residual stresses has been noticed in the below figures ( Figure 10) when appropriate die geometry is used. Also, in next step, simulated results have shown a better trend with values matching with the experimental results, obtained by employing 30 mm/sec die movement velocity. Parameters play a significant role in the variation of residual stresses. A better value of the obtained residual stresses has been noticed in the below figures ( Figure 10) when appropriate die geometry is used. Also, in next step, simulated results have shown a better trend with values matching with the experimental results, obtained by employing 30 mm/sec die movement velocity. There is evidence of inhomogeneous deformation of swaged Zr-4 alloy as the feed rate is increased from 0.7 m/min to 2 m/min in the case of the 25% swaged rod, measured in the radial direction. As seen from Figure 10a  Parameters play a significant role in the variation of residual stresses. A better value of the obtained residual stresses has been noticed in the below figures ( Figure 10) when appropriate die geometry is used. Also, in next step, simulated results have shown a better trend with values matching with the experimental results, obtained by employing 30 mm/sec die movement velocity. There is evidence of inhomogeneous deformation of swaged Zr-4 alloy as the feed rate is increased from 0.7 m/min to 2 m/min in the case of the 25% swaged rod, measured in the radial direction. As seen from Figure 10a

Deep Hole Drilling
For the 25% swaged samples (18.2 mm swaged rod), residual stress results are shown in Figure  11a-c. It is observed that as the diameter decreases through swaging operation, and the magnitude of the residual stresses increases as is evident from Figure 11d for the 16.2 mm rod (0.7 m/min feed rate). This is attributed to high cold work, which would lead to higher flow stress as well as higher residual stress.

Deep Hole Drilling
For the 25% swaged samples (18.2 mm swaged rod), residual stress results are shown in Figure 11a-c. It is observed that as the diameter decreases through swaging operation, and the magnitude of the residual stresses increases as is evident from Figure 11d for the 16.2 mm rod (0.7 m/min feed rate). This is attributed to high cold work, which would lead to higher flow stress as well as higher residual stress. For the 25% swaged alloy, deep hole drilling results indicate the compressive stresses at the surface and tensile stresses at center for 0.7 m/min, 1.25 m/min and 2 m/min feed rates. The magnitude of compressive residual stress increases with increasing feed rates. This is due to the higher amount of deformation induced during higher feed rates. Deep hole drilling measurements show the compressive stress at the surface and the tensile stress at the center for both the axial and radial directions. Figure 12a-c shows the kernel average misorientation of 25% swaged Zr-4 alloys at different For the 25% swaged alloy, deep hole drilling results indicate the compressive stresses at the surface and tensile stresses at center for 0.7 m/min, 1.25 m/min and 2 m/min feed rates. The magnitude of compressive residual stress increases with increasing feed rates. This is due to the higher amount of deformation induced during higher feed rates. Deep hole drilling measurements show the compressive stress at the surface and the tensile stress at the center for both the axial and radial directions. Figure 12a-c shows the kernel average misorientation of 25% swaged Zr-4 alloys at different feed rates. This result shows that the amount of strain increases with increasing feed rates (from 0.7 m/min to 2 m/min). That also supports the trend of surface residual stress measured through experimental and numerical methods. For the 25% swaged alloy, deep hole drilling results indicate the compressive stresses at the surface and tensile stresses at center for 0.7 m/min, 1.25 m/min and 2 m/min feed rates. The magnitude of compressive residual stress increases with increasing feed rates. This is due to the higher amount of deformation induced during higher feed rates. Deep hole drilling measurements show the compressive stress at the surface and the tensile stress at the center for both the axial and radial directions. Figure 12a-c shows the kernel average misorientation of 25% swaged Zr-4 alloys at different feed rates. This result shows that the amount of strain increases with increasing feed rates (from 0.7 m/min to 2 m/min). That also supports the trend of surface residual stress measured through experimental and numerical methods.

Experimental vs. Numerical
Experimental and numerical results obtained from the FE analysis of residual stresses are plotted together and compared in Figure 13a-c.

Experimental vs. Numerical
Experimental and numerical results obtained from the FE analysis of residual stresses are plotted together and compared in Figure 13a The residual stress predicted by numerical modelling is in tandem with the deep hole drilling technique for the 25% (18.2 mm rod) swaged samples. For the 25% swaged alloy, deep hole drilling results reveal compressive stresses at the surface and tensile stresses at the center. For the 25% swaged alloy, numerical modelling showed (in the radial direction, Z) compressive stresses at the surface and tensile stresses at the center. Both experimental and numerical methods show similar residual stress magnitudes quantitatively at the surface of the swaged rod. For the 0.7 m/min feed rate, deep hole drilling and numerical methods have shown 6 MPa and −0.9 MP, respectively. In the case of 1.25 m/min feed rate, deep hole drilling technique showed −52 MPa and −41.4 MPa from the numerical method. For the 2 m/min feed rate, deep hole drilling showed −86 MPa and −74 MPa with the numerical method. Both numerical and experimental methods predicted a similar trend in variation of tensile residual stresses but they differ slightly in their magnitude. These results are for swaged samples at different feed rates and it suggests that increasing the feed rate increases the compressive residual stresses on the surface while tensile stresses are nearly similar on increasing the feed rates. With these results, we suggest that 1.25 m/min feed rate is good for producing swaged Z-4 alloy, because of the good combination of mechanical properties of tensile strength (535.50 MPa), ductility (19.84%) and fracture toughness (44.37 MPa·m 0.5 ) reported in our earlier work [24]. Further annealing treatment will decrease the tensile residual stress at the center of Zr-4 alloy rod. Numerical results The residual stress predicted by numerical modelling is in tandem with the deep hole drilling technique for the 25% (18.2 mm rod) swaged samples. For the 25% swaged alloy, deep hole drilling results reveal compressive stresses at the surface and tensile stresses at the center. For the 25% swaged alloy, numerical modelling showed (in the radial direction, Z) compressive stresses at the surface and tensile stresses at the center. Both experimental and numerical methods show similar residual stress magnitudes quantitatively at the surface of the swaged rod. For the 0.7 m/min feed rate, deep hole drilling and numerical methods have shown 6 MPa and −0.9 MP, respectively. In the case of 1.25 m/min feed rate, deep hole drilling technique showed −52 MPa and −41.4 MPa from the numerical method. For the 2 m/min feed rate, deep hole drilling showed −86 MPa and −74 MPa with the numerical method. Both numerical and experimental methods predicted a similar trend in variation of tensile residual stresses but they differ slightly in their magnitude. These results are for swaged samples at different feed rates and it suggests that increasing the feed rate increases the compressive residual stresses on the surface while tensile stresses are nearly similar on increasing the feed rates. With these results, we suggest that 1.25 m/min feed rate is good for producing swaged Z-4 alloy, because of the good combination of mechanical properties of tensile strength (535.50 MPa), ductility (19.84%) and fracture toughness (44.37 MPa·m 0.5 ) reported in our earlier work [24]. Further annealing treatment will decrease the tensile residual stress at the center of Zr-4 alloy rod. Numerical results could be improved further by using more refined meshes and minor time steps in order to match closely with experimental results. To improve the die design, die movement velocity with or without lubricant during swaging could also contribute to obviate the detrimental effect of residual stress in the alloy.

Conclusions
Numerical and experimental investigations are carried out to investigate the influence of the feed rate in the radial direction on residual stresses of a swaged zircaloy-4 rod. The following conclusions are made based on the present work: 1.
Numerical and experimental results both show a similar trend of residual stresses from the surface to center on the 25% swaged rod at different feed rates (0.7 m/min, 1.25 m/min and 2 m/min). Almost the same magnitude of residual stress on the surface as well as center of the 25% swaged rod are observed.

2.
Deep hole drilling results show an increase in magnitude of residual stress on increasing the percentage of cold work. This is due to increase in the amount of deformation, which leads to an increase in yield strength at higher cold work on zirconium alloys. 3.
The inhomogeneous deformation observed from the surface to the center (in the radial direction) on the 25% swaged rod (0.7 m/min, 1.25 m/min and 2 m/min) as observed in residual stress magnitude for both experimental and numerical methods.

4.
It can be concluded from the load curve that at the same time all four dies exert an equal amount of load over the work piece, it results in a uniform homogeneous deformation in the axial direction of the swaged Zr-4 rod.

5.
The axial (X) and radial (Z) residual stresses on the surface are compressive and thus, can help in preventing crack propagation and can enhance product life. On the other hand, the stresses towards the center of the rod are tensile and thus, may assist the propagation of possible cracks. Hence, there is a need to undergo heat treatment before use.
Author Contributions: G.S. and B.K., conceptualized, procured and analyzed data, and wrote the manuscript; K.I.V.N. and U.K.A., provided the methodology (the swaging process) and resources for samples; M.M.M., provided the facilities for hole drilling and carrying out the experiments; R.J., directed the research, wrote the manuscript and provided comments in the manuscript. All authors have read and agreed to the published version of the manuscript.