1. Introduction
Zirconium alloy products, due to their physical and mechanical properties, are used in industries in which high requirements are placed on mechanical properties and corrosion resistance in aggressive environments. One of the main areas of use for zirconium alloy products is as construction elements for the production of fuel assemblies (nuclear reactor cores). The use of alloys on baize zirconia for the manufacture of the construction elements of poly rods is related to the small active cross-section for the absorption of thermal neutrons. They are characterized by the following additive properties: impermeability of fission products (compared to steel, they have a 15-fold smaller active cross-section for neutron capture) and resistance to radiation damage; low activity under neutron irradiation; sufficient mechanical strength;a low creep rate at high temperatures; sufficient thermal conductivity;andvery slow reactions with fuel and coolant (water), i.e., significant resistance to corrosion [
1,
2,
3]. In addition to meeting the requirements for ensuring high mechanical and anti-corrosion properties, the products must have a uniform structure without metallurgical discontinuities [
2]. The heterogeneity of the structure in the finished products may lead to stress concentrations, which may cause failure of fuel assemblies, which in turn may eventually lead to accidents [
4].
There are many zirconium-based industrial alloys used in the nuclear industry. In general, the alloys are differentiated by the presence and content of elements such as Nb, Sn, Fe, Cr, Ni, and O. The choice of an alloy is related to the purpose of a given constructional element for a specific type of reactor. PWR zirconium with 1% Nb(M5) is used for making the plugs for cladding tubes of reactor fuel assemblies [
5].
The technology for manufacturing tubes and rods from Zr–1%Nb alloy includes the following operations (
Figure 1) [
2,
6]. During the hot forging process, the cast structure is fragmented, and the cross-section of the forging is reduced. The next stage is hot extrusion, during which there is a significant change in the cross-section. The next stage is cold rolling with intermediate and final annealing, during which a specific structure and the properties of the finished products are formed.
Authors [
6,
7] have stated that, in the process of ingot melting, the formation of a piping cavity in the axis of the ingot and of central porosity is inevitable. The formation of metallurgical discontinuities is related to the gradient of the crystallization rate in relation to the volume of the ingot. After melting, the ingots are subjected to ultrasonic testing to find metallurgical discontinuities. However, it is known from the authors’ experience that there may be central porosities with diameters of up to Ø10 mm in the axial part, as well as metallurgical voids with diameters of Ø1.0 mm at a certain distance from the axis of the ingot.
An analysis of the literature shows that methods for closing foundry voids during hot elongation of steel bars and aluminum and magnesium alloys are known [
8,
9,
10,
11,
12,
13,
14,
15], while no workswere found on closing foundry voids in the elongation operation of zirconium bars with the exception of the authors’ work [
16].
Banaszek et al. [
15] carried out numerical modeling of the closure of foundry voids during the forging operation of Mg alloys bars and physical verification of the obtained results.
In the work [
15], the conditions were determined for closing foundry voids in the forging process of M5 alloy rods with the use of rhombic and flat anvils. To do so, the investigators used a physical simulator of metallurgical processes, the Gleeble 3800, and FORGE software, with which computer simulations of this process were carried out. The use of a combination of computer and physical modeling to determine the parameters of the hot elongation process is effective in terms of materials and time costs [
14,
17,
18].
In the work, the authors showed the results of studies of the influence of the geometry of flat trapezoidal and rhombic trapezoidal anvils on the closure of foundry voids in the process of hot forging of M5 alloy rods. Based on the testing results, recommendations were made for the technology of hot forging processes, ensuring the closing of foundry voids occurring in ingots made of M5 alloy.
2. Purpose and Scope of Work
The objective of the study was to formulate recommendations for the technology of forging elongation that would ensure the closing of foundry voids occurring in zirconium alloy, observed after the vacuum process of double remelting in an electric arc furnace. The authors proposed that the operations of hot forging elongation should be performed with the use of flat trapezoidal and rhombic–trapezoidal tools.
To achieve the goals of the study, tests investigatingthe forging of ingots in two anvil assemblies were carried out. Modeling of the elongation investigation was performed using the Forge program based on FEM, with the distributions of the hydrostatic pressure and effective strain being determined on a cross-section of the deformed zirconium alloy after each elongation process.
To test the forging process, the datafrom plastometric tests of the zirconium alloy were obtained using the GLEBLLE 3800, based on which diagrams showing the relation between stress and real deformation were developed, and the coefficients of the plasticizing stress function were determined.
Based on the results of numerical studies, guidelines for the technology of forging shaped anvils for closing foundry voids were developed.
The test results should help to improve the structural and mechanical properties of zirconium alloy during the appropriate selection of the tool shape and technological parameters of the forging process.
3. Methodology of Experimental Research
The zirconium alloy applied for the investigation was an M5 alloy with the following composition and chemical configuration: zirconium with 1.1% Nb, 0.05% Fe, and 0.6% O.
The impact of thermomechanical parameters of the hot elongation of the M5 alloy rod on the plasticizing stress
σp was determined using a Gleeble 3800 plastometer, 323 NY-355, Poestenkill, NY 12140, USA (
Figure 2). Specimens of the M5 alloy in a crystallized state (
Figure 3), each of a working portion with a diameter of 10 mm and a height of 12 mm, were used for the tests.
The tests were carried out in the temperature range of 770–950 °C [
2,
16]. On the basis of the Zr–Nb phase equilibrium scheme at the 870–950 °C temperature range, the structure of the metal is characterized by grains in the shape of the A2 unit cell. At forging temperatures of 750–870 °C, the structure of Zr is in the zone of α + β [
2]. During the deformation of the zirconium with 1% Nb alloy sample in the specified temperatures, the limit value of the deformation decreases by half [
19,
20,
21,
22]. In this article, the authors show that the forging process of these alloys at temperatures less than 750 °C results in the occurrence of cracks.
Tests were carried out for the deformation speed range from 0.5 s
−1 to 5.0 s
−1, even though the average deformation speed for the hot elongation of the Zr alloys rods on the pressure is 0.5 s
−1 [
16]. However, in local places of the deformed rods in the shaped tools, the deformation speed is slightly higher, which is why it was decided to expand the research area.
The functionality of this process is the ability to obtain relatively large strain values (up to ε = 1.2 [
19]); moreover, it is the most favorable strain state for describing the properties of the alloys during plastic deformation [
17].
Origin software was used to process the final test results. This software allowed for a comprehensive analysis of the experimental data obtained, reducing the number of samples tested to two for a single point with an error limit less than 2–3%.
To use of the real results from plastometric tests of the M5 Zr alloys, an approximation of flow curves
was carried out with the use of the generalizing connection–functions of Henzel A. and SpittelT. [
19]:
where
is stress, Tis temperature,
is strain,
is strain rate, A, and m1–m9 are function coefficients.
The approximation of the testing effects was carried out according to the method [
14] in FORGE@2023 software, Transvalor S.A., E-Golf Park, 950 avenue Roumanille, CS 40237 Biot, 06904 Sophia Antipolis cedex, France.
4. Testing Results
Figure 4,
Figure 5 and
Figure 6 illustrate the
curves of the zirconium alloys deformed at T = 770–950 °C for the deformation speed
from 0.5 s
−1 to 5.0 s
−1. These figures show increases in the temperature of test T from 770 °C to 950 °C, during which the value of the plasticizing stress σ
p decreases about 2.5 times. An increase in the deformation speed from 0.5 s
−1 to 5.0 s
−1 led to an increase in the plasticizing stress σ
p, and at T = 770 °C, this increase amounted to ~29%. At T = 850 °C, this increase was ~36%, and at T = 950 °C, it was 43%. The conducted research showed that the increase in the deformation speed affected the plasticizing stress σ
p, which increased, and in addition, the magnitude of this increase was correlated with the temperature.
The plasticity curves for the tested temperature range varied. For the deformation speed = 0.5 s−1 in the tested temperature range, the curves reached the max. stress value for a programmed deformation value, and as the temperature increased, the max. plasticizing stress value σp moved toward actual smaller values.
The data shown in
Figure 5 demonstrate that the curve obtained at deformation speeds of
= 0.5 s
−1 and
= 5.0 s
−1 are of a different character. The max. value of yield stress during deformation speed
= 0.5 s
−1 was observed for the value of real deformation ε = 0.12, and with a continued increase in deformation, a decrease in the value of yield stress was obtained. For the deformation speed
= 5 s
−1, the value of yield stress increases to the value of real deformation ε = 0.3, where the curve flattening effect takes place.
Figure 5.
The flow curves of the M5 alloy, obtained for the temperature T = 850 °C in the deformation speed range from 0.5 to 5.0 s−1 using the Gleeble 3800 metallurgical simulator.
Figure 5.
The flow curves of the M5 alloy, obtained for the temperature T = 850 °C in the deformation speed range from 0.5 to 5.0 s−1 using the Gleeble 3800 metallurgical simulator.
This difference may be due to the dynamic softening process being delayed with the increase in strain rate. Dynamic recrystallization is inhibited, and the softening of the material proceeds in accordance with the mechanism of dynamic polygonization. At the temperature T = 862 °C, a polymorphic transformation takes place in zircon, and the crystal lattice changes from the hexagonal A3 lattice to the spatially centered cubic lattice A2 [
2]. The non-monotonicity of the curve may be caused by the beginning of the polymorphic transformation [
11].
Figure 6.
The flow curves of the M5 alloy, obtained for T = 950 °C in the deformation speed range at 0.5 to 5.0 s−1 using the Gleeble 3800.
Figure 6.
The flow curves of the M5 alloy, obtained for T = 950 °C in the deformation speed range at 0.5 to 5.0 s−1 using the Gleeble 3800.
As an effect of approximation of the investigational zirconium alloy curves, Hensel–Spittel coefficients were determined for the temperature and speed of the following plastic processing conditions: T = 770–950 °C, and
= 0.5–5.0 s
−1 (
Figure 7).
The multifactorial Equation (1), using the coefficients shown in
Figure 8 in the process of computer simulation, describes the changes in the value of the plasticizing stress of the M5 alloy depending on the thermomechanical parameters of the tested technological processes. The approximation results for the M5 alloy in the temperature range of T = 770–950 °C and deformation rates from 0.5 s
−1 to 5.0 s
−1 are shown in
Figure 8 in the form of three-dimensional graphs.
It can be assumed that the coefficients of the approximating Function (1) are wellchosen if the average error does not exceed 15% [
18]. For the range of parameters examined in the article, the average error of approximation is ~10%.
5. Methodology of Conducted Numerical Research
This article provides a theoretical analysis of the process of elongating Zr alloy bars for two different anvils to determine the possibility of closing foundry voids in Zr alloys. Tools were characterized by various geometrical surfaces of the deformation valley, which resulted in obtaining various directions and orientations of the pressure and friction forces in the deformed Zr alloy in the deformation place. The various directions and orientations of forces influenced the mechanism of welding of the foundry voids inside the deformed bars. Constant parameters of the forging operation were assumed, such as a forging initiation temperature of 950 °C, a value of relative compressionof 35%, and a value of the feed angleof 8 mm/s [
23,
24].
The parameters for the lengthening operation were adopted based on the authors’ experience in this field [
15,
16].
The computer program FORGE was applied to model the forging process in shaped anvils [
25]. Forges enable the runningofthermomechanical simulations of, among other processes, plastic forming processes [
26]. A description of the temperature, force, stress, strain, and thermomechanical and frictional laws used in the investigation appeared in [
21].
For thermal calculations, the Galerkin equation was used, while the strengthening curves were approximated by the Henzel–Spitel equation (
Figure 8).In this paper, to simulate the elongation operation, a thermo–visco–plastic model of the deformed body, which is based on the theory of large plastic deformations, was used. To generate the grid of finite elements, tetrahedral elements with the bases of triangles were used. In the generated model input, a number of nodes equal to 8958 was used for the simulation, while the number of tetrahedral elements adopted for the simulation was 40,760. The coefficient of friction between the surface of the tools and the bars, adoptedaccording to Coulomb’s law, was μ = 0.3. The value of the friction coefficient used for numerical modeling wasnot determined experimentally, but results from forging practice were used for many years.For the free forging of steels and alloys, friction coefficient values between 0.24 and 0.46 are assumed.
The following boundary conditions were assumed during the numerical tests: Indirect heat transfer between material and tool –“Steel Hot Medium.tef”; Friction coefficient between material and tools–“High.tff”;and Heat transfer coefficient between material and environment—“Air. tef”. The heat transfer coefficient connecting the tools and the zirconium alloy was assumed to be α = 10,000 W/m2K and the heat transfer coefficient between the zirconium alloy and the environment to be equal to λ = 10 W/m2K. The environment temperature was 25 °C, and the tool temperature was 250 °C. The starting temperature was 950 °C. In all forging steps, a relative 35% crumple was assumed. The speed of the upper tool was v = 8 mm/s, while the lower tool was assumed to be stationary.
Table 1 shows the values of the boundary and initial conditions used in the numerical model.
In the FORGE computer program, the diffusion model does not exist. The inference of foundry voids closing is based on the values of hydrostatic pressure and the temperature of the rods being elongated. Therefore, in the forging process, the aim is to achieve the maximum possible values of hydrostatic pressure within the foundry voids, as well as to maintain a high temperature close to that at the beginning of forging. Through the process of closing the foundry voids, there are high values of hydrostatic pressure on their right and left sides, while there is no hydrostatic pressure on their top and bottom sides (positive values of average stress) [
13].
Zr alloy bars with a diameter of 100 mm and a length of l = 50 mm were deformed in two compositions of tools in four forging transitions.
The contour and size of the tools used in the investigation are shown in
Figure 9 and
Figure 10, respectively.
This work describes a graph of the forging operation for flat trapezoidal tools (
Figure 9). During in first step, flat trapezoidal tools were used, in which a bar heated to a temperature of 950 °C was deformed with a relative reduction of 35%, and then the bar was rotated 90°. In the next step, with the bar using the same anvils, relative reduction of 35% was again applied, and it was rotated 90° clockwise. Before the third step, the bar was heated to a temperature of 950 °C because, after the second step, the temperature of the bar decreased to 750 °C, and according to the practice of forging Zr alloys, it is impossible to carry out further steps of the forging process. Subsequently, the next steps were already carried out in flat tools with a relative reduction of 35%. After the third pass, the bar was rotated 90° clockwise.
For the rhombic trapezoidal tools (
Figure 10), the procedure was analogous to the forging scheme described above.
The geometry and distribution of the modeled foundry voids on the end face of the model Zr alloy bar are shown in
Figure 11.
The initial model for the forging process was a cylinder with a diameter of 100 mm and a length of 50 mm, in which nine foundry voids were modeled with a length of charge (50 mm); the first axial one had a diameter of 10 mm, and the other eight, with a diameter of 1 mm, were distributed around the circumference of the model charge 25 mm from its axis (
Figure 11). The charge modeled in this way corresponded to the shape of an ingot obtained after double remelting in a vacuum in an arc furnace. Artificially modeled holes simulated internal casting discontinuities, such as central porosity and casting voids. Descriptions of these defects can be found in the literature [
23,
24,
25,
26].
The model inputs and tools were drawn with AutoCad software. The foundry voids were drawn as cylinders inside the model input. Using the various tools in the program, the foundry voids were interpreted by the AC program as holes. The model diagram with artificially simulated foundry voids was exported to the Forge program as a file with the “stl” extension. The same process was performed for all of the tools’ compositions. In Forge, using the stl@meshing and volume@meshing tools, a 2D triangular grid was useful to the input and artificially modeled foundry voids, followed by a 3D tetrahedral grid. In total, for the correct simulation of the closure of foundry voids and their good formation, the built-in tools FoldsDetection andSelfContactwere used. Due to the use of the shown tools during the FEM simulation in the places of welded foundry voids, there were no calculation errors or mesh deterioration. During welding, the foundry voids of the nodal points of the tetrahedral elements in close contact connected; thus, the modeled foundry voids were closed. During modeling, there were no errors in the form of interpenetration of nodes of tetrahedral elements around closed foundry voids.
On the basis of the simulations of the closing process of foundry voids during the forging of bars made of zirconium alloy in flat and rhombic tools, the volume values of the remaining unwelded foundry voids were determined. The procedure was as follows.After the simulation of the forging process was completed, the stlfile with a deformed bar was exported from the Forge program to the RinoCeros program, where the tetrahedral element grid was separated in this program, and the bar contour elements were removed. After this step, only tetrahedral elements that constituted the contours of the unwelded discontinuities remained. The next step was to merge the previously broken grid of tetrahedral elements. In the last step, using the physical parameters tool in the RinoCeros program, the whole volume of foundry voids for individual forging process was obtained.
After the bar forging process was carried out in two anvil compositions, the values of the hydrostatic pressure at the nodal points of the tetrahedra, constituting the contour of the unwelded foundry voids, were determined as the arithmetic means of the occurring stresses.The same procedure was performed for the value of the effective strain. All the results shown in the article are subject to some errors resulting from the need to analyze the numerical and data and the results obtained.
8. The Influence of the Anvil Shape on Changes in the Volume of Discontinuities, Hydrostatic Pressure, and the Effective Strain in the Zr Alloy in the First Two Steps
Figure 44,
Figure 45 and
Figure 46 show diagrams of the total volume of unclosed foundry voids, mean values of hydrostatic pressure, and arithmetic mean values of the effective strain occurring around unclosed foundry voids for the anvils analyzed in the article, divided into the first and second steps.
To determine the influence of the deformation basin shape on the closing of internal foundry voids, the authors also presented the research results in [
18].
By analyzing the data presented in
Figure 44,
Figure 45 and
Figure 46, the authors try to answer the question of how the shape of the hitting surfaces of the tools affects the closing of internal foundry voids.
Figure 44.
The total volume of unclosed foundry voids in the volume of a forged Zr alloy bar for two forging steps in relation to individual tools’ compositions.
Figure 44.
The total volume of unclosed foundry voids in the volume of a forged Zr alloy bar for two forging steps in relation to individual tools’ compositions.
By analyzing the data presented in
Figure 44, it can be seen that the best results in terms of the closing of foundry voids during the forging process of Zr alloy bars after the first and second steps were obtained for rhombic–flat and flat tools [
18]. After the second step, the values of the total volumes of unclosed foundry voids ranged from 66.40 mm
3 to 82.80 mm
3, which constituted only 1.5% to 2% of the volume of the initially modeled discontinuities. It is worth noting that the total starting value of the volume of all modeled foundry voids was 4239 mm
3. The conducted research shows that the closing of foundry voids can also be performed with the use of rhombic trapezoidal and flat trapezoidal anvils. For these tools, after the first step, an observed reduction in the total volume of foundry voids amounted to 105.53 mm
3 and 1181.38 mm
3, respectively. After the second step for these tools, a decrease in the total value of the volume of unclosed foundry voids was noted. The volumes of unclosed foundry voids were 99.64 mm
3 and 157.62 mm
3, respectively, which accounted for 2.3% and 3.7% of the volume of the initially modeled foundry voids.
Figure 45.
The arithmetic mean value of the hydrostatic pressure around the unclosed foundry voids in relation to the individual tools’ compositions.
Figure 45.
The arithmetic mean value of the hydrostatic pressure around the unclosed foundry voids in relation to the individual tools’ compositions.
Figure 45 also presents data on the value of hydrostatic pressure obtained after forging processeswith flat and shaped tools presented in [
18] to compare the hydrostatic pressure values obtained for different anvils. The data presented in
Figure 45 show that, for the assembly of flat–trapezoidaland rhombic–trapezoidal anvils, the value of the hydrostatic pressure after the first step was similar and amounted to 26 MPa and 17 MPa, respectively. After the second step, a decrease in this value was observed for both the first and the second anvils’ composition, while for the rhombic–trapezoidal anvils, the decrease was about 16 MPa.
Based on the analysis of the data shown in
Figure 45, it can be seen that the highest value of hydrostatic pressure after the first and second steps was obtained for the rhombic–flat and flat tools. The high value of the hydrostatic pressure around the foundry voids is not unambiguous in relation to the total value of the closed foundry voids (
Figure 44).
Figure 46.
The arithmetic mean value of the deformation intensity occurring around unclosed foundry voids in relation to individual tool compositions.
Figure 46.
The arithmetic mean value of the deformation intensity occurring around unclosed foundry voids in relation to individual tool compositions.
Figure 46 compares the data on the value of the effective strain obtained after the forging process in flat and shaped tools presented in [
18]. The study of the information in
Figure 46 shows that the most favorable arithmetic mean value of the effective strain after two steps was obtained for flat tools and was equal to 1.70. The arithmetic mean values of the effective strain calculated for the bar forging in the remaining tool compositions, except for the rhombic–trapezoidal tools, had much lower values. In the implementation of the first two steps with the rhombic–trapezoidal anvils, the arithmetic mean value of the deformation intensity after the first transition was 1.44, and after the second transition, it was 1.31. These values were, however, slightly lower than the values obtained for the deformation with flat anvils. This slightly smaller difference resulted in significant inhibition of the welding process of metallurgical discontinuities after the second forging transition. For trapezoidal anvils, the worst results were obtained in terms of closing of the foundry voids (
Figure 44).
Considering the operations of the forging processes of Zr alloy with five different tool compositions, it can be stated that it is favorable, in terms of closed foundry voids, to maintain the effective strain in the area of the closed foundry voids at the level of 1.70 and higher and the hydrostatic pressure at the level of 45 MPa and higher.