Next Article in Journal
Comprehensive Assessment from Optimum Biodiesel Yield to Combustion Characteristics of Light Duty Diesel Engine Fueled with Palm Kernel Oil Biodiesel and Fuel Additives
Next Article in Special Issue
Hot Deformation Behavior of Non-Alloyed Carbon Steels
Previous Article in Journal
A Comparison Study on the Strengthening and Toughening Mechanism between Cu-Bearing Age-Hardening Steel and NiCrMoV Steel
Previous Article in Special Issue
Investigation of the Influence of Open-Die Forging Parameters on the Flow Kinetics of AZ91 Magnesium Alloy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Developments towards a Multiscale Meshless Rolling Simulation System

1
Institute of Metals and Technology, Lepi Pot 11, 1000 Ljubljana, Slovenia
2
Faculty of Mechanical Engineering, University of Ljubljana, Aškerčeva 6, 1000 Ljubljana, Slovenia
*
Author to whom correspondence should be addressed.
Materials 2021, 14(15), 4277; https://doi.org/10.3390/ma14154277
Submission received: 24 June 2021 / Revised: 17 July 2021 / Accepted: 26 July 2021 / Published: 30 July 2021
(This article belongs to the Special Issue Hot Deformation and Microstructure Evolution of Metallic Materials)

Abstract

:
The purpose of the present paper is to predict the grain size of steel during the hot-rolling process. The basis represents a macroscopic simulation system that can cope with temperatures, stresses and strains of steel in a complete continuous rolling mill, including reversible pre-rolling and finishing rolling with several tenths of rolling passes. The grain size models, newly introduced in the present paper, are one-way coupled to the macro-scale calculations performed with the slice model assumption. Macroscale solution is based on a novel radial basis function collocation method. This numerical method is truly meshless by involving the space discretization in arbitrarily distributed nodes without meshing. A new efficient node generation algorithm is implemented in the present paper and demonstrated for irregular domains of the slice as they appear in different rolling passes. Multiple grain size prediction models are considered. Grain size prediction models are based on empirical relations. Austenite grain size at each rolling pass as well as the ferrite grain size at the end of rolling are predicted in this simulation. It is also shown that based on the rolling schedule, it is highly likely that recrystallization takes place at each pass throughout a continuous rolling mill. The simulation system is coded as a user-friendly computer application for industrial use based on programming language C# and an open source developer platform NET and runs on regular personal computers the computational time for a typical rolling simulation is usually less than one hour and can thus be straightforwardly used to optimize the rolling mill design in a reasonable time.

1. Introduction

In the present paper, the macroscopic simulation system, coping with a complete continuous rolling mill, including reversible pre-rolling and finishing rolling with several tenths of rolling passes [1,2,3,4,5], is upgraded to deal with grain size prediction. The main aim is to calculate, in addition to the temperature and deformation fields, the austenite grain size at each rolling pass and ferrite grain size at the end. The microscopic grain-size prediction models are one-way coupled to the macroscopic thermo-mechanical simulation. A previously developed 2D numerical model [6], based on the slices that allow only plane strain deformation and are aligned perpendicular to the rolling direction, is used. The slices on which 2D stresses and temperature fields are computed are illustrated in Figure 1.
Shaping any metal, such as in the hot rolling process, requires consideration of large plastic deformation. An overview of related multiscale physics is given in the following paragraphs. The start of this plastic deformation process is triggered when the dislocations can move inside a grain. This threshold is the yield stress, and its value depends on the dislocation density. By applying deformation and heat treatment, dislocation density and the grain sizes face a significant variation.
If the strain rate and temperature are high enough, a different phenomenon, dynamic recrystallization, occurs. The dynamic recrystallization depends on the critical strain. As a result, new grain sizes appear, either partially or fully. However, not all recrystallization takes place immediately after exceeding the critical strain. After a particular time, meta-dynamic crystallization may take place.
When the deformation is small, or there is no deformation, a softening process may occur due to sufficient time at high temperatures. The softening can again lead to partial or complete static recrystallization. Static recrystallization usually leads to larger grains with higher temperatures [7]. The stored energy in a material is proportional to the defect density. Therefore, an increase in grain size decreases the total grain boundary area and as a result, the stored energy is reduced.
The numerical prediction of grain size constitutes an essential aspect of microstructure studies [8]. Smaller grains may be desired since they lead to higher strength because it is harder for dislocations to pass over the grain boundaries during plastic deformation. Or the other way around for the ductility. The steel industry needs to predict the mechanical properties of their final product, which is only possible with the microstructure analysis [9].
In this paper, multiple grain size prediction models [10,11,12,13,14,15,16,17] were coupled with the macroscopic rolling simulation system. The basis of the micro-scale simulation used in the present paper is based on the model predicted by Hodgson and Gibbs [10] with minor adjustments. Multiple grain size prediction models for different types of steels were gathered. These models consist of empirical correlations to determine the critical strain and grain size after meta-dynamic, dynamic or static recrystallization. All the necessary input data for these models were obtained from the macroscopic 2D slice model results. Checking for possible recrystallization type and predicting the corresponding austenite grain size was performed immediately after obtaining the macroscopic results for each slice. After the hot rolling process was completed, the rolled steel was still hot and considered mainly in the austenite phase. Based on the models found in the literature, ferrite grain size prediction could be made as a function of the cooling process. Ferrite transformation temperature is calculated based on the material composition, and in the simulations, the steel is cooled down just below that transformation temperature. The discussion of how to implement the microscopic models can be found in textbooks [1,9]. Since it was not possible to obtain all the correlations for each steel grade separately for all the micro-scale simulation sections, a combination of those relations for various steel types were gathered and also applied. This was just to show that in this way, both macro and micro simulation results of hot rolling could be obtained easily by using a single meshless simulator. A specific model was used for each sub-section of micro simulation, but the necessary parameters were material dependent, and they could also be user-defined. In the results section, the first macro scale simulation results by using 16MnCrS steel were given. Later the same results were used as inputs for multiple micro-scale simulations. Comparison of different grain size prediction models were shown to distinguish the behaviour of each model. Each collocation point in a computational domain represents a small area with a uniform grain structure.
Over the years, a comprehensive rolling simulation system has been developed for industrial use. The supporting steel factory provides rolling schedules in use, and macro-scale simulations have been up to now run by a successful macroscopic hot rolling simulation system [18]. The meshless solution of the macro-scale deformation was chosen here since it is very straightforward to distribute collation nodes over a computational domain. Therefore, it was easy to implement numerically and also did not require any background integration. Meshless methods were multiple times proven to solve thermal [19] and mechanical [20] problems efficiently and accurately. Recently, micro-scale simulation capabilities based on empirical relations found in the literature were added, and a one-way coupling between macro and micro-scale models was formed. The motivation for the present research was understanding and simulating the material behaviour in between the passes, which is usually ignored during a multi-pass rolling. In the majority of the simulations, the initial material definition was used throughout the simulation regardless of the number of deformation steps, time and temperature. This work clearly shows the possibility of recrystallization at each pass separately during hot rolling and how much the material properties may change, based on the change of the grain size.
In the present paper, the possibility of complete or partial recrystallization is modelled. This makes more realistic the assumption of complete recrystallization at each roll pass, as assumed in some publications [17]. Furthermore, the effects of the chemical composition of the steel billet on its macro response and grain structure are also investigated.

2. Materials and Methods

It is pointed out by Sellars [8] that the experimental tests or industrial trials are not enough to understand all the details of the hot rolling process since it is impossible to control or monitor all the parameters of rolling in industrial conditions and even in laboratory conditions. Therefore, the best solution suggested by Sellars [14] is combining the numerical simulations with empirical models based on experiments. A macro-scale numerical prediction of the deformation of steel with heat transfer during hot rolling has been previously developed [6,18]. To achieve the microstructure evolution just as described in Figure 2, we need to know the initial grain configuration and empirical models that depend on the involved steel grade.

2.1. Macro Scale Solution

Thermo-mechanical simulation of rolling is based on 2D slices aligned towards the rolling direction. The slices were simulated one after another, and both thermal and mechanical models were solved in a coupled way
L T σ + b = 0 ,
ρ c p T t = ( k T ) + Q ˙ ,
where L , σ , b are derivative operator matrix, stress vector and body force vector, respectively. ρ , c p , T , k , Q ˙ stand for the density, specific heat, temperature, thermal conductivity and internal heat generation, respectively. The meshless numerical procedure is based on the interpolation of an unknown field θ ( p ) at position p = p x i x + p y i y , interpolated with radial basis shape functions ψ ( p ) over seven neighbouring nodes ( N ) as shown in Figure 3,
θ ( p ) = i = 1 N ψ i ( p ) α i ,
where α i are the collocation coefficients that need to be determined. The partial derivatives of the unknown field can be calculated as
θ ( p ) x j = i = 1 N ψ i ( p ) x j α i .
The temperature in the thermal model and the mechanical model’s displacement field are interpolated with the multi-quadric (MQ) radial basis functions (RBF). The details of the solution of these two models have been previously published in [21].
The boundary condition for the thermal model is of the Robin type,
k T n Γ = h ( T T r e f ) .
where k is thermal conductivity, n Γ is unit normal vector on the surface, h is heat transfer coefficient, and T r e f is the reference temperature. The majority of complications during simulation usually appear at the contact boundaries between the strand and the groove. The heat flux there becomes very high in comparison with the non-contact boundaries. At the contact boundaries, T r e f is considered as the roll’s surface temperature. The boundary conditions for the mechanical model are described both with the prescribed traction ( τ ¯ ) and the prescribed displacement ( u ¯ )
τ ¯ = μ p ,   u ¯ = u ,
where μ is the coefficient of friction, p is stress acting on the boundary, and u is the displacement. In some rare cases, when the groove surface has a unique geometry, prescribed traction boundary condition may not be appropriately satisfied, and might lead to erroneous results. To overcome this issue, during the calculation step at first contact, an artificially sticking boundary condition was applied to some of the boundary points. This resulted in no relative motion between the roll surface and the boundary in one calculation step. Later all the boundary nodes received proper boundary conditions defined in Equation (6) as schematically seen in Figure 4.
Interpolating an unknown field with regular node distribution over a rectangular domain works perfectly with meshless Local Radial basis Function Collocation Method (LRBFCM) and is also easy to implement numerically. However, severe deformation simulations, such as rolling or forging, require redistribution of the computational nodes. In other words, a solution with the regular node distribution becomes impossible in large deformations. An essential advantage of the LRBFCM is that it gives stable results also with scattered node distribution. A comparison of the results is shown later in this paper where a rectangular domain is interpolated with regular and scattered node distribution.
LRBFCM interpolates the temperature field in the thermal model, and the governing equation is solved through explicit time stepping. The time step is chosen as 1 ms. The nodes of the mechanical and thermal models coincide. The solution is obtained through local interpolation for each influence domain without creating a global matrix. In the mechanical model, interpolation of displacement filed over each local influence domain can be assembled into a global solution matrix ( A ) by replacing the local collocation coefficients with the inverse of the local interpolation matrix and the local displacement vector. Hence, the system of equations can be written in the following form
A U = B ,
and solved for global column matrix of displacements U = [ u x 1 , u y 1 , , u x N , u y N ] T with the size of two times the number of collocation nodes ( 2 N ). The adjacent vector B includes the boundary values. As a result of the nonlinear material model of steel, the solution matrix A depends on the current state of displacements. The solution can only be achieved iteratively. In this paper, a direct iteration method was chosen due to its simplicity and to avoid assembling the Jacobian matrix. The local solution of the system of equations leads to a sparse matrix A , and with the help of the Intel math kernel library [22], it can be solved efficiently by direct iteration with iteration index j .
A ( U j ) U j + 1 = B ,

2.2. Micro Scale Solution

In the present work, the critical strain, the time for dynamic, meta-dynamic or static recrystallization and the grain sizes after different types of recrystallization implemented through a sequence of different empirical models. It is assumed that each collocation point in a computational domain represents a small area with a uniform grain structure. Grain size prediction models were obtained for the austenite phase. Solution steps were adopted from Hodgson and Gibbs [10], as shown in Figure 5.

2.2.1. Critical Strain

The micro-scale modelling starts by checking first if the strain field in some of the points, some of the sub-areas or globally exceeds the critical strain value ( ε c )
ε c = A c ( ε ˙ exp Q d R T ) Z p D q ,
where Z denotes the Zener–Hollomon parameter ( s 1 ), used in a multitude of empirical relations, ε ˙ is strain rate, Q d is activation energy ( J / mol   K ), R is ideal gas constant 8.314 J / mol   K , T is temperature in Kelvin and D is the initial austenite grain size diameter. The coefficients in Equation (9) for different steel grades are listed in Table 1. The dynamic or meta-dynamic recrystallization is triggered at the points where the critical strain value is exceeded. This can occur locally or on the global scale of the domain.

2.2.2. Dynamic Recrystallization

The dynamic recrystallization (DRX) occurs when the critical strain is reached, but time is not long enough for meta-dynamic recrystallization. The grain size is predicted from the following equation
D D R X = A d Z r ,
The coefficients of Equation (10) are defined in Table 2 below.

2.2.3. Meta-Dynamic Recrystallization

The meta-dynamic crystallization is also called post dynamic recrystallization. Its effect on the grain size is similar to the dynamic recrystallization; however, there is a time-based criterion that includes strain rate. In this paper, the following definition is used [15],
t M D R X = 1.166 ε 0.8 ,
as developed by Manohar et al. If the current time t at any deformation step exceeds the criterion ( t > t M D R X ), it is considered that the meta-dynamic recrystallization (MDRX) occurs. If this is the case, the grain sizes must be recalculated from the following equation,
D M D R X = A m Z r .
The grain size prediction model for the meta-dynamic and dynamic recrystallization is the same, but the coefficients are different. The constants of Equation (12) are defined in Table 3 below.

2.2.4. Static Recrystallization

There is a time criterion for static recrystallization (SRX) to occur, and in this paper, a model defined by Zhang et al. [16] for high carbon steel is used. It has the following form,
t S R X = 1.944 × 10 4 ε 1.0 D 0 0.6 ( 3.6 ε ˙ ) 0.28 exp ( 6900 R T ) .
when the time during the simulation is larger than ( t > t S R X ) then, static recrystallization may be expected. The grain size after the recrystallization is obtained from
D S R X = A s ε m D 0 r
The coefficients for Equation (14) are defined in Table 4.

2.2.5. Grain Growth

After complete recrystallization, the grains will grow until the next partial or complete recrystallization. In the present paper, the following grain growth model to calculate the new grain size is considered for three types of steel.
D G G 0 = D 0 n + k s t exp ( Q G G R T ) .
The coefficients for Equation (15) are defined in Table 5.

2.2.6. No-Recrystallization Temperature

It is known from the experiments that recrystallization may not occur every time, even though the conditions are met. This limitation comes from the temperature value. If it is below a certain value ( T n r ), recrystallization will not take place. This temperature value is defined from the following equation given in [23]
T n r = A n ε m ε ˙ n t k ,
where A n = 905 , m = 0.045 , n = 0.006 , k = 0.024 and t is the time in seconds. These values are for high strength low alloy steel, experimentally obtained in [24]. The time in Equation (16) is accounted for only during the deformation process.

2.2.7. Ferrite Grain Size Prediction

After the steel leaves the last rolling stand, it is cooled in a controlled way at the cooling bed. During the cooling process the austenite phase transforms into ferrite and perlite, especially for C-Mn steels. In this paper, only the ferrite transformation is considered, and the temperature at which the transformation occurs is called Ferrite Transformation Temperature (FTT). (FTT) for any type of steel in °C is calculated based on the following equation given in [25]:
T F T T = 910 310 [ % C ] 80 [ % M n ] 20 [ % C u ] 15 [ % C r ] 80 [ % M o ] [ % N i ] .
The steel has to cool down until the maximum temperature goes below the FTT, and with that specific cooling rate at that temperature, the ferrite grain size can be predicted based on the following equation.
D F 0 = ( α 0 + α 1 C e q ) + ( α 3 + α 4 C e q ) T ˙ 0.5 + α 4 ( 1 exp ( α 5 D A ) ) .
The necessary coefficients are defined in Table 6. It is important here not to ignore the residual strain ( ε r ) which has an impact on reducing the ferrite grain size. If complete recrystallization occurred at each roll pass, the residual strains were obtained only from the last rolling stand. The following model by Sellars and Beynon [26] is considered in the simulations to calculate the final ferrite grain size.
D F = D F 0 ( 1 0.45 ε r ) .

2.3. Coupling of Micro and Macro Models

Macro-scale simulation is followed by micro-scale simulation for each slice at a time except for the final ferrite grain size, which is calculated at the end. The major steps of the macro and micro-simulation steps and the coupling between them can be seen in Figure 5. The results were visualized for each slice, in particular at the exit from each rolling stand. The temperature, strain and strain rate fields are the inputs of the micro-scale model and were obtained from the macro-scale simulation results. This solution procedure represents a one way coupled system. When the simulation is completed, switching between or rearranging the micro-scale models and recalculating them is possible. Each micro-scale simulation took a few seconds. It was unnecessary to rerun the macro-scale simulation, which is time-consuming, to make a sensitivity study with different micro-scale model parameters. It was also possible to add entirely new micro-scale correlations into the database. There are two different ways of visualization of the micro-scale results. First is the type of recrystallization to see if it is partial or complete and when it happens. The second outcome is the grain size at each slice position, scaled from the minimum to the maximum.

3. Node Positioning

After the continuous casting in a typical steel production line, the strand proceeds to a reheating furnace to remove the cast dendritic structures and homogenize most alloying elements [2]. The strand afterwards deforms in the reversing rolling mill to obtain the desired pre-shape for the continuous rolling mill. The initial rectangular shape for the finishing rolling was not a perfect rectangle, but it had curved corners. We obtained the proper initial shape for use in the finishing rolling by using the reversing rolling mill simulation results. The regular node distribution was straightforward to generate and effective for the perfect rectangular shape. However, the regularly distributed nodes did not correctly provide the appropriate results near the curved edges.
Most of the rolling schedules used in the simulations may be categorized into two groups when rolling the billets. In the first group, oval grooves were used to achieve more oval or round-shaped forms. In the second group, primarily flat rolls were used to get more slab-like shapes. Regular node distribution was adequate for most of the slab simulations; however, regular nodes faced stability issues for a schedule that turns a rectangular billet into a round bar. Previously we used an Elliptic Node Generation (ENG) algorithm [27,28] for positioning the nodes in the domains with curved boundaries. This algorithm successfully maps a rectangular shape to an arbitrary shape with four sides and ensures high orthogonality. The downside of this method represents a mapped corner. This corner is potentially most tricky during any interpolation, just like in a rectangular domain. However, if we use an irregular node arrangement, this corner complexity will be lost completely. This provides a much more extensive deformation range with stable results.

3.1. Regular vs. Irregular Node Positioning

A fast, irregular node generation algorithm [29] is numerically implemented in the rolling simulation system. Additionally, the first row of internal points, just next to the boundaries, is aligned in the opposite direction of the unit normal of the nearest boundary point.

3.2. Node Repelling Algorithm

Due to the different texture of internal points near the boundaries and the rest of the internal points, an iterative node adjustment algorithm was used. In this way, the neighbouring nodes have a more uniform distance to the central nodes. The adjustment is based on a node repelling algorithm. The vector component of displacement Δ u i , due to repelling, is defined as
Δ u i = C j = 1 N max s p i p i j ( p i p i j ) n .
C is a smoothing factor usually taken between 0 and 1. The higher the node density, the smaller it should be. In this paper it is taken as ( 1 / 6 ) . s is also a similar factor taken as 1, and reduced to 0.8 when at least one neighbouring point is located at a boundary. N M A X is a maximum number of the nearest nodes to consider without the boundary nodes. It was taken as 45 during the first iteration, and it was slowly reduced in each iteration until it reached seven nodes inside an influence domain.

4. Numerical Results

In this section, multiple numerical results from the represented multiscale rolling simulation system (Institute of metals and technology, version 7.6.4, Ljubljana, Slovenia) are illustrated. Macro-scale results are compared with regular and irregular collocation nodes and also compared with Finite Element Method (FEM) results. Later, an actual rolling schedule which consists of eight rolling stands, is applied. Micro-scale results are demonstrated in terms of grain size and recrystallization throughout the rolling schedule. The necessary values and definitions are given in Table 7.

4.1. Regular vs. Irregular Node Positioning—Slab

An initial slab with 217 × 45 mm2 is flat-rolled into 34.5 mm height. The material is elastic with Young’s modulus 106 Pa and Poisson’s ratio 0.35. A regular node distribution with 546 nodes and irregular with 571 nodes are used, as shown in Figure 6. The deformation is calculated in 11 steps in both cases, and the coefficient of friction is taken as 0.1.
The results in terms of effective strain can be seen in Figure 7. The results are almost identical; therefore, for a slab that undergoes the flat rolling, irregular node generation has no distinct benefit.

4.2. Regular vs. Irregular Node Distribution—Real Shape

More distinct differences occur when starting from an initial shape with curved corners, as obtained from the reversing rolling mill. The node arrangement is kept regular with ENG [28] and irregular with the algorithm [29], as shown in Figure 8. The regular and irregular node distributions involve 396 and 400 nodes, respectively.
The instability issue at the edge of the first contact is still carried out until the end of the simulation in the regular node arrangement. In Figure 9, it is shown that irregular node arrangement made a better interpolation when the boundary shape diverges from regular (rectangular) geometry. Smoothed corners created problems with regular node arrangement. Comparing with FEM results in Figure 9c, the overshoot near the corner in Figure 9b is very clear. However, this can be overcome by using an irregular node arrangement where the corner information is completely lost as shown in Figure 9a.

4.3. Irregular Meshless Node Distribution vs. FEM

The same flat rolling example as in Section 4.1 is chosen here with an ideal plastic material model. Effective stress–effective strain relation is σ ¯ = 589 ε ¯ 0.2 . A comparison is made with FEM code [30], and a good agreement is found. An example of a shear strain comparison is given in Figure 10 below. As expected, a good agreement has been found between the meshless results with irregular node arrangement and FEM.

4.4. Rolling Schedule with Eight Rolling Stands

In this section, a rolling schedule from the industry with specific grooves are shown in Figure 11 The same was used in all subsequent simulations. The initial size of the billet was 60.48 by 60.48 mm, where the corners are rounded with a radius of 11 mm. An irregular node distribution, the same as demonstrated in Section 4.2, is used. This rolling schedule aimed to obtain a 45 × 30 mm2 end shape within a half millimetre error range.

4.4.1. Effective Strain Results

Effective strain and temperature fields are the most critical inputs of the micro-scale models. The corresponding position and time of each slice have to be known throughout the simulation. Figure 12 represents the effective strain result of slices positioned at the exit of each roll pass. These results have shown the total effective strain through each of the rolling passes. Therefore these fields contain the highest possible values for each pass. These results are also crucial since they have a significant impact on the micro-scale simulation.

4.4.2. Recrystallization Results at Each Pass

In the simulations, the deformation occurred so rapidly that the time for meta-dynamic or static recrystallization was never reached. However, at each pass, soon after the first contact, a complete dynamic recrystallization was always achieved except for the model developed by Yada [12]. In that model, the dynamic recrystallization was partial at each pass. To contrast, the results obtained by Yada’s model are shown in Figure 13. In the three other models, there was always complete dynamic recrystallization.

4.4.3. Recrystallization Results at Each Pass with Yada Model

The results of partial dynamic recrystallization obtained by the rolling simulation system are shown in Figure 13 for slices at each roll pass’s exit.

4.4.4. Grain Size Comparisons of Different Models

A rolling schedule is used as described in Section 4.4 and Figure 11. Grain size simulation results based on the model by Sellars and Whiteman [14] is shown in Figure 14.
The model by Yada [12] is shown in Figure 15, the model by Manohar, et al. [15] is shown in Figure 16 and the model by Hodgson and Gibbs [10] is shown in Figure 17. The grain growth model is obtained in [10]. In micro-scale simulations, each collocation point in a computational domain represents a small area with a uniform grain structure.

4.4.5. Ferrite Grain Size for 16MnCr5

For the calculation of the final grain size in the ferrite phase, 16MnCr5 steel is considered. This steel grade has 0.16% carbon, 1.15% manganese and 0.95% chromium. The ferrite transformation temperature is based on the model given by Hodgson and Gibbs, calculated as 754.15 °C, and the user interface is displayed in Figure 18. The real-time cooling to below FTT took 14 min; however, the computational time took twice as long. The final ferrite grain size and comparison between the two models is shown in Figure 19. A strong dependence on residual strain on the ferrite grain size is observed.

4.4.6. No Recrystallization Temperature

Based on the model expressed by Yang, et al. [24], no-recrystallization temperatures are calculated and compared with minimum and average temperature values for each slice position through the rolling simulation. The result is shown in Figure 20 below, and it is clear that for hot rolling at 1100 °C, the temperature of steel was always above the no-recrystallization temperature.
The positions of roll passes which are used in Figure 14, Figure 15, Figure 16, Figure 17 and Figure 20 are given in Table 8 below. The origin of this coordinate ( z = 0 ) is considered to be the initial slice’s position.

5. Discussion

Micro-scale calculations are based on a multitude of empirical correlations obtained from the experiments in the last four decades. The simulation system is capable of checking for meta-dynamic, dynamic or static recrystallization. When recrystallization occurs, a new grain size is calculated. If there is a full recrystallization, grain growth will start. At the end of rolling, cooling is simulated with the thermal model until it reaches the ferrite transformation temperature. The simulation ends with ferrite grain size prediction. A comparison of different models based on maximum and minimum grain sizes is shown. The accuracy of those results are limited to the accuracy of the input data which is simulated by the macro-scale model and the accuracy of the empirical micro-scale models. It is known from the industrial user experience that the accuracy of the macro-scale results is highly satisfactory. Therefore, if the empirical models are based on enough experiments, the expected grain-related results should be between the maximum and minimum simulated values.
Another outcome of this research is to be able to come up with new grain size prediction models in the future. This is why all the definitions are implemented as user-defined functions. Numerous additional experimental results are needed to obtain this goal.
The simulation system is currently limited to austenite and ferrite grain size predictions and assumes only one phase at a time. This will be improved in the future with multiple initial phase definitions and phase transformation criteria. The micro-scale models are also limited to the adequacy of those coefficients used. All the necessary coefficients for each steel type may not be found. It is a long and complicated experimental process to define all those parameters by analyzing numerous grain size data obtained by either using electron microscope [31] or ultrasound [32].

6. Conclusions

In the present paper, a rolling simulation system capable of predicting the grain structure is described. First, a macroscopic thermo-mechanical simulation is performed, based on a 2D solution over each cross-sectional slice. Displacement, strain, stress and temperature fields are calculated considering the defined process parameters (rolling schedule), heat transfer and material properties. The solution is achieved through a truly meshless radial basis function collocation method. A new irregular node generation algorithm is used, which is shown, by FEM comparison, to perform better on irregular domains. This method requires only the distribution of collocation nodes over a domain, without classical meshing. It is possible to represent the complex domain shapes more accurately with the irregular node distribution as with the regular node distribution. Due to the high possibility of recrystallization at each pass, strain filed during the last rolling stand leads to relatively smaller ferrite grain size. It was shown in the simulations that dynamic recrystallization always took place at each roll pass. Meta-dynamic and static recrystallizations did not occur due to time restrictions. It was also expressed that no-recrystallization temperatures most probably does not stop recrystallization during hot rolling of steel above 1100 °C. A list of all the correlations used in this paper may be found in Appendix A.

Author Contributions

Methodology, supervision and manuscript editing B.Š.; software, validation and writing U.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Slovenian Grant Agency by grant numbers P2-0162 and L2-2609 and Štore-Steel Company (www.store-steel.si).

Data Availability Statement

Data sharing not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

C-MnLow CHigh CC-Mn-VC-Mn-TiC-Mn-NbHSLA
Critical strain A 6.82 × 10−44.76 × 10−4 4.9 × 10−44.9 × 10−4
Q d 312,00066,500 312,000300,000
p 0.130 0.150.17
q 00 10.3
Dynamic recrystallization A 16,00022,600180016,000
r −0.23−0.27−0.15−0.23
Q d 312,000267,100312,000300,000
Meta-dynamic recrystallization A 26,000 160026,00025,000
r −0.23 −0.11−0.23−0.23
Q d 312,000 312,000300,000312,000
Static recrystallization A 1 0.5
m −0.5 −0.67
r 0.4 0.67
Grain growth m 7104.5
k s 1.45 × 10272.6 × 10284.1 × 1028
Q G G −400−437−435
No-recrystallization temperature A 905
m −0.045
n −0.006
k 0.024
Ferrite grain size α 0 −0.422.6 1.4
α 1 6.37−57 0
α 2 24.23 5
α 3 −590 0
α 4 2222 22

References

  1. Lenard, J.G.; Pietrzyk, M.; Cser, L. Mathematical and Physical Simulation of the Properties of Hot Rolled Products, 1st ed.; Elsevier: Oxford, UK, 1999. [Google Scholar]
  2. Lenard, J.G. Primer on Flat Rolling, 1st ed.; Elsevier: Oxford, UK, 2007. [Google Scholar]
  3. Roberts, W.L. Hot Rolling of Steel, 1st ed.; Marcel Dekker Inc.: New York, NY, USA, 1983. [Google Scholar]
  4. Lee, Y.; Choi, S.; Kim, Y.H. Mathematical model and experimental validation of surface profile of a workpiece in round-oval-round pass sequence. J. Mater. Process. Technol. 2000, 108, 87–96. [Google Scholar] [CrossRef]
  5. Solod, V.; Kulagin, R.; Beygelzimer, Y. A local approach to simulating bar forming in pass rolling. J. Mater. Process. Technol. 2007, 190, 23–25. [Google Scholar] [CrossRef]
  6. Hanoglu, U.; Šarler, B. Multi-pass hot rolling simulation using a meshless method. Comput. Struct. 2018, 194, 1–14. [Google Scholar] [CrossRef]
  7. Dong, X.; Cheng, J.; Hongyang, Z.; Dongying, J.; Miaoyong, Z. A new study on the growth behavior of austenite grains during heating processes. Sci. Rep. 2017, 7, 3968–3980. [Google Scholar]
  8. Sellars, C.M. Modelling microstructural development during hot rolling. Mater. Sci. Technol. 1990, 6, 1072–1081. [Google Scholar] [CrossRef]
  9. Goh, C.H.; Dachowicz, A.P.; Allen, J.K.; Mistree, F. Integrated Computational Materials Engineering (ICME) for Metals Concept Case Studies; John Willey and Sons: Hoboken, NJ, USA, 2018. [Google Scholar]
  10. Hodgson, P.D.; Gibbs, R.K. A mathematical model to predict the mechanical properties of hot rolled C-Mn and microalloyed steels. Iron Steel Inst. Jpn. Int. 1992, 32, 1329–1338. [Google Scholar] [CrossRef] [Green Version]
  11. Laasraoui, A.; Jonas, J.J. Prediction of temperature distribution, flow stress and microstructure during the multipass hot rolling of steel plate and strip. Iron Steel Inst. Jpn. Int. 1991, 31, 95–105. [Google Scholar] [CrossRef] [Green Version]
  12. Yada, H. Prediction of microstructural changes and mechanical properties in hot strip rolling. In Proceedings of the Metallurgical Society of the Canadian Institute of Mining and Metallurgy, 1st ed.; Pergamon Press: Kidlington, UK, 1988; pp. 105–119. [Google Scholar] [CrossRef]
  13. Kwon, H.C.; Lee, Y.; Kim, S.Y.; Woo, J.S. Numerical prediction of austenite grain size in round-oval-round bar rolling. Iron Steel Inst. Jpn. Int. 2003, 43, 676–683. [Google Scholar] [CrossRef] [Green Version]
  14. Sellars, C.M.; Whiteman, J.A. Recrystallization and grain growth in hot rolling. J. Met. Sci. 1979, 13, 187–194. [Google Scholar] [CrossRef]
  15. Manohar, P.A.; Lim, K.; Rollett, A.D.; Lee, Y. Computational exploration of microstructural evolution in a medium C-Mn steel and applications of rod mill. Iron Steel Inst. Jpn. Int. 2003, 43, 1421–1430. [Google Scholar] [CrossRef] [Green Version]
  16. Zhang, Y.; Zhang, H.; Wang, G.; Hu, S. Application of mathematical model for microstructure and mathematical property of hot rolled wire rods. Appl. Math. Model. 2009, 33, 1259–1269. [Google Scholar] [CrossRef]
  17. Biglou, J.; Lenard, J.G. A study of dynamic recrystallization during hot rolling of microalloyed steels. CIRP Ann. 1996, 45, 227–230. [Google Scholar] [CrossRef]
  18. Hanoglu, U.; Šarler, B. Hot rolling simulation system for steel based on advanced meshless solution. Metals 2019, 9, 788. [Google Scholar] [CrossRef] [Green Version]
  19. Zhou-Jia, F.; Li-Wen, Y.; Qiang, X.; Chein-Shan, L. A boundary collocation method for anomalous heat conduction analysis in functionally graded materials. Comput. Math. Appl. 2021, 88, 91–109. [Google Scholar]
  20. Li, Y.; Li, J.; Wen, P.H. Finite and infinite block Petrov-Galerkin method for cracks in functionally graded materials. Appl. Math. Model. 2018, 68, 306–326. [Google Scholar] [CrossRef]
  21. Vertnik, R.; Šarler, B. Meshfree local radial basis function collocation method for diffusion problems. Comput. Math. Appl. 2006, 51, 1269–1282. [Google Scholar]
  22. Intel Math Kernel Library (Intel MKL); Version 2018.0; Intel Corporation: Santa Clara, CA, USA, 2018; Available online: https://software.intel.com/oneapi/onemkl (accessed on 20 November 2019).
  23. Bai, D.Q.; Yue, S.; Maccagno, T.; Jonas, J.J. Static recrystallization of Nb and Nb-B steels under continuous cooling conditions. Iron Steel Inst. Jpn. Int. 1996, 36, 1084–1093. [Google Scholar] [CrossRef]
  24. Yang, H.; Zhou, X.; Liu, Z.; Wang, G. Determination of No-recrystallization temperature for NB-bearing Steel. J. Wuhan Univ. Technol. Mater. Sci. Ed. 2016, 31, 644–647. [Google Scholar] [CrossRef]
  25. Fang, X.; Fan, Z.; Ralph, B.; Evans, P. Influence of accelerated cooling on microstructure and mechanical properties of C–Mn steels. Mater. Sci. Technol. 2002, 18, 47–53. [Google Scholar] [CrossRef]
  26. Sellars, C.M.; Beynon, J. Microstructural development during hot rolling of titanium microalloyed steels. In International Conference on High Strength Low Alloy Steels; University of Wollongong: Wollongong, NSW, Australia, 1985; pp. 142–150. [Google Scholar]
  27. Thompson, J.F.; Soni, B.K.; Weatherill, N.P. Handbook of Grid Generation, 1st ed.; Taylor & Francis Inc.: Abingdon, UK, 1998. [Google Scholar]
  28. Hanoglu, U.; Šarler, B. Simulation of hot shape rolling of steel in continuous rolling mill by local radial basis function collocation method. Comput. Model. Eng. Sci. 2015, 109, 447–479. [Google Scholar]
  29. Fornberg, B.; Flyer, N. Fast generation of 2-D node distributions for mesh-free PDE discretizations. Comput. Math. Appl. 2015, 69, 531–544. [Google Scholar] [CrossRef]
  30. DEFORM 2D; Version 10.1; Scientific Forming Technologies Corporation: Columbus, OH, USA, 2009; Available online: http://www.deform.com (accessed on 25 May 2021).
  31. Mingard, K.P.; Roebuck, B.; Bennett, E.G.; Gee, M.G.; Nordenstrom, H.; Sweetman, G.; Chan, P. Comparison of EBSD and conventional methods of grain size measurement of hardmetals. Int. J. Refract. Met. Hard Mater. 2009, 27, 213–223. [Google Scholar] [CrossRef]
  32. Palanichamy, P.; Joseph, A.; Jayakumar, T.; Raj, B. Ultrasonic velocity measurements for estimation of grain size in austenitic stainless steel. NDT E Int. 1995, 28, 179–185. [Google Scholar] [CrossRef]
Figure 1. Scheme of slice model used in the macroscopic rolling simulations. A ( t ) is area, v ( t ) is velocity and z ( t ) is position.
Figure 1. Scheme of slice model used in the macroscopic rolling simulations. A ( t ) is area, v ( t ) is velocity and z ( t ) is position.
Materials 14 04277 g001
Figure 2. Picture of one-way macro-micro coupling.
Figure 2. Picture of one-way macro-micro coupling.
Materials 14 04277 g002
Figure 3. Scheme of the domain Ω and boundary Γ discretization with indicated local influence domains consisting of 7 node and for each central node in an influence domain is marked in blue.
Figure 3. Scheme of the domain Ω and boundary Γ discretization with indicated local influence domains consisting of 7 node and for each central node in an influence domain is marked in blue.
Materials 14 04277 g003
Figure 4. Boundary nodes during contact with groove surface and applied boundary conditions.
Figure 4. Boundary nodes during contact with groove surface and applied boundary conditions.
Materials 14 04277 g004
Figure 5. Scheme of the multiscale rolling simulation system. Micro-scale simulation steps are circled with red.
Figure 5. Scheme of the multiscale rolling simulation system. Micro-scale simulation steps are circled with red.
Materials 14 04277 g005
Figure 6. Comparison of different initial node arrangements. Regular node distribution (left), irregular node distribution (right). Av. Temp. stands for average temperature.
Figure 6. Comparison of different initial node arrangements. Regular node distribution (left), irregular node distribution (right). Av. Temp. stands for average temperature.
Materials 14 04277 g006
Figure 7. Comparison of effective strain fields after 10.5 mm of reduction. Av. Temp. stands for average temperature.
Figure 7. Comparison of effective strain fields after 10.5 mm of reduction. Av. Temp. stands for average temperature.
Materials 14 04277 g007aMaterials 14 04277 g007b
Figure 8. Irregular nodes (left) and regular nodes with mapping (right). Av. Temp. stands for average temperature.
Figure 8. Irregular nodes (left) and regular nodes with mapping (right). Av. Temp. stands for average temperature.
Materials 14 04277 g008
Figure 9. Comparison of simulated effective strain fields. Irregular nodes (a), regular nodes with mapping (b) and FEM results (c). The results of all three simulations are practically the same: (a) 70.23 mm, (b) 70.12 mm, (c) 70.15 mm. Av. Temp. stands for average temperature.
Figure 9. Comparison of simulated effective strain fields. Irregular nodes (a), regular nodes with mapping (b) and FEM results (c). The results of all three simulations are practically the same: (a) 70.23 mm, (b) 70.12 mm, (c) 70.15 mm. Av. Temp. stands for average temperature.
Materials 14 04277 g009
Figure 10. Comparison of strain fields. Meshless solution with irregular node arrangement is on top, FEM solution is shown at the bottom. The results of both calculations are almost the same: (a) 233.52 mm, (b) 233.48 mm. Av. Temp. stands for average temperature.
Figure 10. Comparison of strain fields. Meshless solution with irregular node arrangement is on top, FEM solution is shown at the bottom. The results of both calculations are almost the same: (a) 233.52 mm, (b) 233.48 mm. Av. Temp. stands for average temperature.
Materials 14 04277 g010
Figure 11. Demonstration of the roll shapes of each of the roll passes of an 8-roll rolling schedule. (a) first roll pass, (b) second roll pass, (c) third roll pass, (d) forth roll pass, (e) fifth roll pass, (f) sixth roll pass, (g) seventh roll pass, (h) eighth roll pass. Axis inc. is the length of axis increments.
Figure 11. Demonstration of the roll shapes of each of the roll passes of an 8-roll rolling schedule. (a) first roll pass, (b) second roll pass, (c) third roll pass, (d) forth roll pass, (e) fifth roll pass, (f) sixth roll pass, (g) seventh roll pass, (h) eighth roll pass. Axis inc. is the length of axis increments.
Materials 14 04277 g011aMaterials 14 04277 g011b
Figure 12. Simulation results of effective strain at the exit of each roll pass. (a) at first roll pass, (b) at second roll pass, (c) at third roll pass, (d) at forth roll pass, (e) at fifth roll pass, (f) at sixth roll pass, (g) at seventh roll pass, (h) at eighth roll pass. Axis inc. is the length of axis increments. Av. Temp. stands for average temperature.
Figure 12. Simulation results of effective strain at the exit of each roll pass. (a) at first roll pass, (b) at second roll pass, (c) at third roll pass, (d) at forth roll pass, (e) at fifth roll pass, (f) at sixth roll pass, (g) at seventh roll pass, (h) at eighth roll pass. Axis inc. is the length of axis increments. Av. Temp. stands for average temperature.
Materials 14 04277 g012aMaterials 14 04277 g012b
Figure 13. Recrystallization based on the model given by Yada [12] at the exit of each roll pass. (a) at first roll pass, (b) at second roll pass, (c) at third roll pass, (d) at forth roll pass, (e) at fifth roll pass, (f) at sixth roll pass, (g) at seventh roll pass, (h) at eighth roll pass. Possible recrystallization types are dynamic (DRX), meta-dynamic (MDRX) and static (SRX). Axis inc. is the length of axis increments and Av. Temps stands for average temperature.
Figure 13. Recrystallization based on the model given by Yada [12] at the exit of each roll pass. (a) at first roll pass, (b) at second roll pass, (c) at third roll pass, (d) at forth roll pass, (e) at fifth roll pass, (f) at sixth roll pass, (g) at seventh roll pass, (h) at eighth roll pass. Possible recrystallization types are dynamic (DRX), meta-dynamic (MDRX) and static (SRX). Axis inc. is the length of axis increments and Av. Temps stands for average temperature.
Materials 14 04277 g013aMaterials 14 04277 g013b
Figure 14. Simulation of maximum and minimum grain size based on the model by Sellars and Whiteman [14].
Figure 14. Simulation of maximum and minimum grain size based on the model by Sellars and Whiteman [14].
Materials 14 04277 g014
Figure 15. Simulation of maximum and minimum grain size based on the model by Yada [12].
Figure 15. Simulation of maximum and minimum grain size based on the model by Yada [12].
Materials 14 04277 g015
Figure 16. Simulation of maximum and minimum grain size based on the model by Manohar et al. [15].
Figure 16. Simulation of maximum and minimum grain size based on the model by Manohar et al. [15].
Materials 14 04277 g016
Figure 17. Simulation of maximum and minimum grain size based on the model by Hodgson and Gibbs [10].
Figure 17. Simulation of maximum and minimum grain size based on the model by Hodgson and Gibbs [10].
Materials 14 04277 g017
Figure 18. Ferrite transformation temperature used in the simulations based on the chemical composition of the steel.
Figure 18. Ferrite transformation temperature used in the simulations based on the chemical composition of the steel.
Materials 14 04277 g018
Figure 19. Comparison of final ferrite grain size based on the models given by Sellars and Whiteman (left) and Hodgson and Gibbs (right). Axis inc. is the length of the axis increments. Av. Temp. stands for average temperature. Vel. Inc. is the increase of the slab’s velocity towards the rolling direction. Slab. Len. Is the length of the rolled slab. LYS is the lower yields stress and UTS is the ultimate tensile strength.
Figure 19. Comparison of final ferrite grain size based on the models given by Sellars and Whiteman (left) and Hodgson and Gibbs (right). Axis inc. is the length of the axis increments. Av. Temp. stands for average temperature. Vel. Inc. is the increase of the slab’s velocity towards the rolling direction. Slab. Len. Is the length of the rolled slab. LYS is the lower yields stress and UTS is the ultimate tensile strength.
Materials 14 04277 g019
Figure 20. Comparison of estimated no-recrystallization temperature values at each rolling pass with the average and the minimum temperature values.
Figure 20. Comparison of estimated no-recrystallization temperature values at each rolling pass with the average and the minimum temperature values.
Materials 14 04277 g020
Table 1. Parameters used in the critical strain definition for carbon–manganese steel in [11,12,14] and low carbon steel in [13].
Table 1. Parameters used in the critical strain definition for carbon–manganese steel in [11,12,14] and low carbon steel in [13].
Coefficients[11][12][13][14]
A c 6.82 × 10−44.76 × 10−44.9 × 10−44.9 × 10−4
Q d 312,00066,500300,000312,000
p 0.1300.170.15
q 000.31
Table 2. Parameters used in the grain size equation for dynamic recrystallization of carbon–manganese steel.
Table 2. Parameters used in the grain size equation for dynamic recrystallization of carbon–manganese steel.
Coefficients[10][12][14][16]
A d 16,00022,600180016,000
r −0.23−0.27−0.15−0.23
Q d 312,000267,100312,000300,000
Table 3. Parameters used in the grain size equation for meta-dynamic recrystallization of carbon–manganese steel [10,14,16], and low carbon steel [10].
Table 3. Parameters used in the grain size equation for meta-dynamic recrystallization of carbon–manganese steel [10,14,16], and low carbon steel [10].
Coefficients[10][13][14][16]
A m 26,00025,000160026,000
r −0.23−0.23−0.11−0.23
Q d 312,000312,000312,000300,000
Table 4. Parameters used in the grain size prediction model when static recrystallization takes place. These parameters are obtained by Zhang [16] for high carbon steel and Sellars [8] for carbon manganese steel.
Table 4. Parameters used in the grain size prediction model when static recrystallization takes place. These parameters are obtained by Zhang [16] for high carbon steel and Sellars [8] for carbon manganese steel.
Coefficients[13][7]
A s 0.51
m −0.67−0.5
r 0.670.4
Table 5. Parameters used in the grain growth calculation when complete recrystallization takes place. These parameters were obtained by Hodgson and Gibbs [10] for three types of steel.
Table 5. Parameters used in the grain growth calculation when complete recrystallization takes place. These parameters were obtained by Hodgson and Gibbs [10] for three types of steel.
CoefficientsC-Mn-VC-Mn-TiC-Mn-Nb
n 7104.5
k s 1.45 × 1027 2.6 × 1028 4.1 × 1028
Q G G −400−437−435
Table 6. Parameters used to calculate ferrite grain size for C-Mn steel.
Table 6. Parameters used to calculate ferrite grain size for C-Mn steel.
Coefficients[10] (%C) + (%Mn)/6 > 0.35[10] (%C) + (%Mn)/6 < 0.35[26] C-Mn-Ti
α 0 −0.422.61.4
α 1 6.37−570
α 2 24.235
α 3 −5900
α 4 222222
α 5 0.0150.0150.015
Table 7. Parameters used in the rolling simulation.
Table 7. Parameters used in the rolling simulation.
General Data
Initial size of the slab60.48 × 60.48 mm
Initial temperature field of the slab1100 °C
Entry velocity towards the rolling direction760 mm/s
Coefficient of friction0.15
Initial grain size10 µm
Material Model
σ ¯ ( ε ¯ , ε ¯ ˙ , T ) = 589 ε ¯ 0.214 ε ¯ ˙ 0.2 exp ( 38 , 000 R T )   MPa , for 16MnCr5 alloyed steel
Thermal Model
Thermal conductivity29 W/mmK
Specific heat630 J/KgK
Density7450 kg/m3
Heat transfer coefficient to air20 W/m2K
Heat transfer coefficient to roll10,000 W/m2K
Roll surface temperature600 °C
Groove and Roll Data
Roll radii230 mm
Roll gaps34.5–26.5–118–22.5–112–20–107.4–18.2 (mm)
Groove geometries
(H-horizontal, V-vertical)
flat (H)—flat (H)—oval box (V)—flat (H)—flat box (V)—flat (H)—flat box (V)—flat (H)
Distance between rolling stands3 m (for the first five rolling stands), 4.5 m (after the fifth rolling stand).
Table 8. Positions of centres of each roll along the rolling direction with reference to the initial slice position.
Table 8. Positions of centres of each roll along the rolling direction with reference to the initial slice position.
Passes12345678
Position on rolling direction (mm)50035006100970013,30020,80025,30029,800
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hanoglu, U.; Šarler, B. Developments towards a Multiscale Meshless Rolling Simulation System. Materials 2021, 14, 4277. https://doi.org/10.3390/ma14154277

AMA Style

Hanoglu U, Šarler B. Developments towards a Multiscale Meshless Rolling Simulation System. Materials. 2021; 14(15):4277. https://doi.org/10.3390/ma14154277

Chicago/Turabian Style

Hanoglu, Umut, and Božidar Šarler. 2021. "Developments towards a Multiscale Meshless Rolling Simulation System" Materials 14, no. 15: 4277. https://doi.org/10.3390/ma14154277

APA Style

Hanoglu, U., & Šarler, B. (2021). Developments towards a Multiscale Meshless Rolling Simulation System. Materials, 14(15), 4277. https://doi.org/10.3390/ma14154277

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop