# Simulation of Solidified Microstructure and Experimental Comparative Study of Twin-Roll Casting Aluminum Alloys

Authors: Teng Ma; Junting Zhang; Xiaochao Cui; Xiaosi Sun
Journal: Advances in Materials Science and Engineering
 (2017)

---

## Abstract

A coupled macro-micro mathematical model of twin-roll casting was established in the study. The continuous solidification process of Al-10Mg aluminum alloys in front of the nip point was numerically simulated by the cellular automaton method and the solidification microstructure, dendritic grain radius, and secondary dendrite arm spacing were obtained through simulation to predict the mechanical property of wedge strips. In order to verify the reliability of the simulation results, the metallographic examination and tensile tests were performed with the as-cast specimens. The results showed that the gain size, the distribution characteristics of various grain regions, dendrite arm spacing, and the yield strengths obtained from simulation were consistent with experimental results.


---

## Body

## 1. Introduction

Twin-roll casting is characterized by the short process, low energy consumption, and low manufacturing cost and widely used to produce aluminum alloy strips with good mechanical properties [1–4]. The solidification structure prepared by continuous casting affects alloy properties as well as the quality and mechanical properties of alloy strips [5–7]. Therefore, it is necessary to develop new rolling technologies and deeply explore the relationship between the casting process, solidification structure, and metal performance.The solidification process of strips can be simulated with the computer simulation technology to study nucleation and grain growth and analyze solidification structure and geometric features. The formation process of the solidification structure can be displayed directly with the simulation technology. Several mathematical models were developed for the twin-roll casting process. Some micro mathematical models were established to simulate the evolution process of the solidification microstructure in a micro region by Liu et al. [8], S. D. Chen and J. C. Chen [9], and Huang et al. [10]. However, the complete solidification process of the strip or the prediction of mechanical properties of the as-cast alloy strip was not reported. The casting-state microstructure of aluminum alloy or other alloys had been explored in the twin-roll casting experiments [11–13]. However, the solidification evolution was seldom explored.According to cellular automata theory, the mathematical model of the solidification structure of twin-roll continuous casting was established in this study. With Al-10Mg aluminum alloys as the study target, the continuous solidification process in front of the nip point was numerically simulated in Procast. The dendritic grain radius and the secondary dendrite arm spacing were obtained to predict the mechanical properties of wedge strips. Finally, the mathematic model was verified through metallographic examination and tensile tests.


## 2. Mathematic Models

### 2.1. Macroscopic Model

Macro temperature field is the basis for simulating microstructures. According to the heat transfer process between the metal liquid, casting roller, and the environment, Xiong et al. proposed the mathematical model based on unsteady heat conduction [14]:(1)ρcp∂T∂t=∂∂xλ∂T∂x+∂∂yλ∂T∂y+∂∂zλ∂T∂z-ρL∂fs∂t,where ρ is alloy density; cp is a specific heat at constant pressure; λ is thermal conductivity; L is latent heat; and fs is solid fraction.


### 2.2. Microscopic Model

#### 2.2.1. Nucleation Model

During the twin-roll continuous casting, the undercooling of the liquid phase is generally far less than the undercooling required for uniform nucleation and liquid-phase surface oxidation, broken dendrites, and heterogeneous nucleation exist in the liquid metal. Therefore, the heterogeneous nucleation model is adopted to describe the solidification process of liquid metal. At a given undercoolingΔT, the continuous nucleation model based on Gaussian distribution and proposed by Gandin and Rappaz [15] is adopted: (2)dndΔT=nmax2πΔT0exp⁡-12ΔT-ΔTNΔT02,where ΔT0 is the standard deviation; nmax is the maximum density of nuclei; and ΔTN is the mean nucleation undercooling.


#### 2.2.2. Growth Model

Twin-roll continuous casting is a near-rapid solidification process. Therefore, the modified KGT model is adopted to describe the constrained dendrite tip growth. The dynamic effect in the dendrite growth is not significant [16], so compositional undercooling and curvature undercooling other than the kinetic degree of the undercooling of dendrite tip are considered in the model. Assuming that the equilibrium distribution coefficient of alloy and diffusion coefficients of solute remain unchanged, the modified KGT model is given as [8](3)ΔT=ΔTc+ΔTR,ΔTc=mc1-11-Ivpe1-k,ΔTR=2ΓR,where c is the solute concentration at the liquid interface; k is the solute partition coefficient; (pe) is the Ivantsov function; and Γ is the Gibbs-Thomson coefficient. According to the stability principle of the dendrite tip, the growth velocity of dendrite tip v and the growing radius of dendrite tip R can be calculated as follows:(4)R=2πΓmGcζc-G,pe=VR2DL,v2π2Γpe2DL2+vmLC01-K0DL1-1-K0Ivpe+G=0.The secondary dendrite spacing is an important parameter to describe characteristics of solidification structure and can reflect the micro segregation and other casting defects. The secondary dendrite spacing,λ2, is given as [17](5)λ2=5.5Mtf1/3,where(6)tf=ΔT0GLR,M=TmΓDLln⁡Cf/C0mLk-1Cf/C0,where Cf is the final composition of liquid metal of dendrite root. The growth velocity of columnar crystals and dendrite tip as well as the dendrite spacing in the process of twin-roll continuous casting can be expressed by the modified KGT model.


## 2.1. Macroscopic Model

Macro temperature field is the basis for simulating microstructures. According to the heat transfer process between the metal liquid, casting roller, and the environment, Xiong et al. proposed the mathematical model based on unsteady heat conduction [14]:(1)ρcp∂T∂t=∂∂xλ∂T∂x+∂∂yλ∂T∂y+∂∂zλ∂T∂z-ρL∂fs∂t,where ρ is alloy density; cp is a specific heat at constant pressure; λ is thermal conductivity; L is latent heat; and fs is solid fraction.


## 2.2. Microscopic Model

### 2.2.1. Nucleation Model

During the twin-roll continuous casting, the undercooling of the liquid phase is generally far less than the undercooling required for uniform nucleation and liquid-phase surface oxidation, broken dendrites, and heterogeneous nucleation exist in the liquid metal. Therefore, the heterogeneous nucleation model is adopted to describe the solidification process of liquid metal. At a given undercoolingΔT, the continuous nucleation model based on Gaussian distribution and proposed by Gandin and Rappaz [15] is adopted: (2)dndΔT=nmax2πΔT0exp⁡-12ΔT-ΔTNΔT02,where ΔT0 is the standard deviation; nmax is the maximum density of nuclei; and ΔTN is the mean nucleation undercooling.


### 2.2.2. Growth Model

Twin-roll continuous casting is a near-rapid solidification process. Therefore, the modified KGT model is adopted to describe the constrained dendrite tip growth. The dynamic effect in the dendrite growth is not significant [16], so compositional undercooling and curvature undercooling other than the kinetic degree of the undercooling of dendrite tip are considered in the model. Assuming that the equilibrium distribution coefficient of alloy and diffusion coefficients of solute remain unchanged, the modified KGT model is given as [8](3)ΔT=ΔTc+ΔTR,ΔTc=mc1-11-Ivpe1-k,ΔTR=2ΓR,where c is the solute concentration at the liquid interface; k is the solute partition coefficient; (pe) is the Ivantsov function; and Γ is the Gibbs-Thomson coefficient. According to the stability principle of the dendrite tip, the growth velocity of dendrite tip v and the growing radius of dendrite tip R can be calculated as follows:(4)R=2πΓmGcζc-G,pe=VR2DL,v2π2Γpe2DL2+vmLC01-K0DL1-1-K0Ivpe+G=0.The secondary dendrite spacing is an important parameter to describe characteristics of solidification structure and can reflect the micro segregation and other casting defects. The secondary dendrite spacing,λ2, is given as [17](5)λ2=5.5Mtf1/3,where(6)tf=ΔT0GLR,M=TmΓDLln⁡Cf/C0mLk-1Cf/C0,where Cf is the final composition of liquid metal of dendrite root. The growth velocity of columnar crystals and dendrite tip as well as the dendrite spacing in the process of twin-roll continuous casting can be expressed by the modified KGT model.


## 2.2.1. Nucleation Model

During the twin-roll continuous casting, the undercooling of the liquid phase is generally far less than the undercooling required for uniform nucleation and liquid-phase surface oxidation, broken dendrites, and heterogeneous nucleation exist in the liquid metal. Therefore, the heterogeneous nucleation model is adopted to describe the solidification process of liquid metal. At a given undercoolingΔT, the continuous nucleation model based on Gaussian distribution and proposed by Gandin and Rappaz [15] is adopted: (2)dndΔT=nmax2πΔT0exp⁡-12ΔT-ΔTNΔT02,where ΔT0 is the standard deviation; nmax is the maximum density of nuclei; and ΔTN is the mean nucleation undercooling.


## 2.2.2. Growth Model

Twin-roll continuous casting is a near-rapid solidification process. Therefore, the modified KGT model is adopted to describe the constrained dendrite tip growth. The dynamic effect in the dendrite growth is not significant [16], so compositional undercooling and curvature undercooling other than the kinetic degree of the undercooling of dendrite tip are considered in the model. Assuming that the equilibrium distribution coefficient of alloy and diffusion coefficients of solute remain unchanged, the modified KGT model is given as [8](3)ΔT=ΔTc+ΔTR,ΔTc=mc1-11-Ivpe1-k,ΔTR=2ΓR,where c is the solute concentration at the liquid interface; k is the solute partition coefficient; (pe) is the Ivantsov function; and Γ is the Gibbs-Thomson coefficient. According to the stability principle of the dendrite tip, the growth velocity of dendrite tip v and the growing radius of dendrite tip R can be calculated as follows:(4)R=2πΓmGcζc-G,pe=VR2DL,v2π2Γpe2DL2+vmLC01-K0DL1-1-K0Ivpe+G=0.The secondary dendrite spacing is an important parameter to describe characteristics of solidification structure and can reflect the micro segregation and other casting defects. The secondary dendrite spacing,λ2, is given as [17](5)λ2=5.5Mtf1/3,where(6)tf=ΔT0GLR,M=TmΓDLln⁡Cf/C0mLk-1Cf/C0,where Cf is the final composition of liquid metal of dendrite root. The growth velocity of columnar crystals and dendrite tip as well as the dendrite spacing in the process of twin-roll continuous casting can be expressed by the modified KGT model.


## 3. Experiment and Model

### 3.1. Experimental

The alloy material used in the study is Al-10Mg aluminum alloy (Table1). Aluminum ingot and magnesium ingot were melted in an airtight furnace according to the proportion of Al-10Mg aluminum alloy and inert gases were injected into the furnace to protect the melting process. The casting roller size is Ø400 mm/Ø300 mm; the roll gap width is 2–4 mm; the gap between nozzle and roll is 0.10~0.15 mm. The chamber was preheated to 350~400°C and the temperature was maintained. Then the twin-roll caster was started and molten metal was poured into the chamber at the pouring temperature of 690°C. The casting speed of the twin-roll caster was 1 m/min.Table 1
Chemical composition of Al-10Mg alloy.
Al
Mg
Si
Fe
Cu
Zn
Other components (Mn + Ni)
89.70
9.90
0.11
0.18
0.01
0.03
0.01When the aluminum alloy strip went through the rolling area, plastic deformation occurred and primary crystal shape was modified. In order to get initial solidification, the wedge strip was taken out after casting and processed to the length of 15 mm (Figure1). After polishing, the strip was corroded in 50% HF. The rolled plane and the longitudinal section of the sample were observed by the metallographic microscope (LeiCADM2700M).Figure 1
Twin-roll cast Al-10Mg alloy wedge-shaped plate specimen: (a) wedge-shaped plate; (b) mounting specimen.
(a)
(b)Three groups of standard samples (10 mm × 4 mm × 3 mm) were taken out along three directions: the horizontal direction, rolling direction, and the angle of 45°. The tensile tests were carried out with universal testing machine to obtain mechanical properties and stress-strain curve of Al-10Mg (strain rateε=5×10-4s-1).


### 3.2. Geometrical Model

In the process of twin-roll continuous casting, the continuous nucleation of the molten metal occurred at the meniscus of the roller when the rollers were reversely rotated. The cast-rolled strip was thickened constantly in the growth of dendrite. When the strip was close to the nip point, the solidification process was completed and the initial structure was formed. Therefore, the continuous solidification process in front of the nip point is numerically simulated in the study. The calculation area is shown in Figure2 and the thermal physical properties of Al-10Mg alloy are shown in Table 2 [18]. The microstructure was analyzed and the mechanical properties were predicted. The simulation parameters were consistent with the production parameters. The thickness of the strip is 3 mm. Figure 2 shows the schematic diagram of the twin-roll casting process. Under the conditions of continuous casting, the heat transfer in the roll barrel direction is far less than that in the other two directions. Therefore, the changes in the grain morphology mainly occurred in the longitudinal section and the rolled plane of the sample, but the grain morphology was the same in the strip width direction. Therefore, in order to simplify the model and reduce the calculation load, a certain section along the roll barrel direction was selected for analysis in the model under the assumption of adiabatic side seal. The strip width is 5 mm in the simplified model.Table 2
Thermal physical properties of Al-10Mg alloy.
Melting temperature (°C)
Specific heat (J·kgK−1)
Density of liquid (kg·m−3)
Thermal conductivity of liquid (W·mk−1)
607
1088.6
2.55
108.9Figure 2
Diagram of the twin-roll casting process.


## 3.1. Experimental

The alloy material used in the study is Al-10Mg aluminum alloy (Table1). Aluminum ingot and magnesium ingot were melted in an airtight furnace according to the proportion of Al-10Mg aluminum alloy and inert gases were injected into the furnace to protect the melting process. The casting roller size is Ø400 mm/Ø300 mm; the roll gap width is 2–4 mm; the gap between nozzle and roll is 0.10~0.15 mm. The chamber was preheated to 350~400°C and the temperature was maintained. Then the twin-roll caster was started and molten metal was poured into the chamber at the pouring temperature of 690°C. The casting speed of the twin-roll caster was 1 m/min.Table 1
Chemical composition of Al-10Mg alloy.
Al
Mg
Si
Fe
Cu
Zn
Other components (Mn + Ni)
89.70
9.90
0.11
0.18
0.01
0.03
0.01When the aluminum alloy strip went through the rolling area, plastic deformation occurred and primary crystal shape was modified. In order to get initial solidification, the wedge strip was taken out after casting and processed to the length of 15 mm (Figure1). After polishing, the strip was corroded in 50% HF. The rolled plane and the longitudinal section of the sample were observed by the metallographic microscope (LeiCADM2700M).Figure 1
Twin-roll cast Al-10Mg alloy wedge-shaped plate specimen: (a) wedge-shaped plate; (b) mounting specimen.
(a)
(b)Three groups of standard samples (10 mm × 4 mm × 3 mm) were taken out along three directions: the horizontal direction, rolling direction, and the angle of 45°. The tensile tests were carried out with universal testing machine to obtain mechanical properties and stress-strain curve of Al-10Mg (strain rateε=5×10-4s-1).


## 3.2. Geometrical Model

In the process of twin-roll continuous casting, the continuous nucleation of the molten metal occurred at the meniscus of the roller when the rollers were reversely rotated. The cast-rolled strip was thickened constantly in the growth of dendrite. When the strip was close to the nip point, the solidification process was completed and the initial structure was formed. Therefore, the continuous solidification process in front of the nip point is numerically simulated in the study. The calculation area is shown in Figure2 and the thermal physical properties of Al-10Mg alloy are shown in Table 2 [18]. The microstructure was analyzed and the mechanical properties were predicted. The simulation parameters were consistent with the production parameters. The thickness of the strip is 3 mm. Figure 2 shows the schematic diagram of the twin-roll casting process. Under the conditions of continuous casting, the heat transfer in the roll barrel direction is far less than that in the other two directions. Therefore, the changes in the grain morphology mainly occurred in the longitudinal section and the rolled plane of the sample, but the grain morphology was the same in the strip width direction. Therefore, in order to simplify the model and reduce the calculation load, a certain section along the roll barrel direction was selected for analysis in the model under the assumption of adiabatic side seal. The strip width is 5 mm in the simplified model.Table 2
Thermal physical properties of Al-10Mg alloy.
Melting temperature (°C)
Specific heat (J·kgK−1)
Density of liquid (kg·m−3)
Thermal conductivity of liquid (W·mk−1)
607
1088.6
2.55
108.9Figure 2
Diagram of the twin-roll casting process.


## 4. Results and Discussion

Figures3 and 4, respectively, show the simulated and experimentally obtained microstructures of rolled plane and longitudinal section, which are 1 mm away from the roll in the continuous casting process. Figure 5 shows the schematic diagram of sample selection for microstructure and texture analysis. On the surface of rolled plane, the size of grains is uniform and equiaxed (Figure 3(a)). The grains near the casting zone are relatively larger (the left side of Figure 3(a)), whereas the grains near the rolling zone are relatively smaller (the right side of Figure 3(a)). In the simulation results, the grain size is 3.2 × 10−8 m2–4.9 × 10−7 m2, whereas metallographic test results indicate that the grain size is 8.0 × 10−8 m2–8.1 × 10−7 m2. The above results prove that the calculated results are well consistent with experimental results. The simulated grain size is a little smaller than the result by metallographic test because of the assumption made in the theoretical models. In addition, the solidification process is a complicated physical process and it is difficult to exactly determine the changing rule of thermal properties with temperature. In Figure 3(b), the solidification structure of longitudinal section from the roll surface into the strip center is, respectively, composed of the fine grain zone, columnar grain zone, and equiaxed crystal zone. The predicted results including the grain shape and grain size are consistent with the results of metallographic tests.Figure 3
Solidification simulation of twin-roll cast Al-10Mg alloy wedge-shaped plate: (a) rolled plane; (b) longitudinal section.
(a)
(b)Figure 4
Solidification simulation of the cross section and lengthwise section of twin-roll cast Al-10Mg alloy wedge-shaped plate: (a) rolled plane; (b) longitudinal section.
(a)
(b)Figure 5
Schematic illustration of sample selection for microstructure and texture analysis.Figure6(a) shows the fraction solid distribution of molten metal at the steady state. Figures 6(b) and 6(c), respectively, show primary dendritic radius distribution and secondary dendrite arm spacing distribution. The primary dendritic radius and the secondary dendrite arm spacing are important parameters for estimating the solidification structure and can reflect various roll casting defects. The solidification process of the solution reaches a stable state after the casting process becomes stable. The primary dendritic radius in the front of columnar grain growth zone reflects the growth condition of columnar grain growth. The smaller radius indicates the more obvious growth tendency. Figure 6(b) shows that the primary dendritic radius of strip surface is smaller than that of strip canter. As shown in Figure 6(c), the secondary dendrite arm spacing is gradually decreased from the outer zone to the inner zone along the thickness direction and is gradually increased from the outer zone to the inner zone along the rolling direction with the solidification. The largest secondary dendrite arm spacing occurs near the nip point and is about 19.55 μm. The measured result (23 μm) in Figure 4(a) is basically consistent with the simulation result, indicating that the model is applicable to predict the dendrite arm spacing.Figure 6
Solid fraction graph, dendritic grain radius graph, and secondary dendrite arm spacing graph of twin-roll cast Al-10Mg alloy wedge-shaped plate: (a) solid fraction graph; (b) dendritic grain radius graph; (c) secondary dendrite arm spacing graph.
(a)
(b)
(c)Roll casting solidification is a sub-rapid solidification process. Compared with the traditional casting structures, the strip structure possesses the smaller grain size, lower segregation degree, and excellent mechanical properties. In this paper, the mechanical properties of the strip were predicted based on the solidification structure obtained by the mathematic model of the casting process. The predicted yield strength of twin-roll cast Al-10Mg strip was 158 Mpa, which was a little higher than the experimental result, 142 Mpa (Figure7), due to the defects such as segregation, gas cavity, and the heterogeneity of grain size. The experimental results of uniaxial tensile mechanical properties of Al-10Mg aluminum alloy casting strip at room temperature are shown in Table 3. The corresponding engineering stress-strain curve is shown in Figure 8. According to the above results, the Al-10Mg casting strip shows the good mechanical properties, obvious yield phenomenon, and higher elongation.Table 3
Tensile properties of twin-roll cast Al-10Mg plate.
Orientation
σ
0.2
/
M
P
a
σ
b
/
M
P
a
δ
0°
142.9
207.2
8.6
90°
121.1
183.3
7.4
45°
108.6
159.4
5.3Figure 7
Simulation results of the yield strength of Al-10Mg alloy plate.Figure 8
Engineering stress-strain curves of twin-roll cast Al-10Mg alloy plate at room temperature.


## 5. Conclusions

(
1
) A mathematical model of the solidification process on twin-roll continuous casting thin strip was developed based on the cellular automata (CA) theory. The continuous solidification process in front of the nip point was numerically simulated by Procast. The solidification microstructure, the grains size, the dendritic grain radius, and the secondary dendrite arm spacing were obtained to predict the mechanical properties of the wedge strip.(
2
) The simulation results were verified by metallographic examination and tensile tests. The simulation results of the microstructure on the wedge strip rolled plane and cross section, which was 1 mm away from the roll, were consistent with the experimental results. The simulation results of dendritic grain radius and the secondary dendrite arm spacing were, respectively, 1133 μm and 16 μm, which were consistent with the experimental results. The predicted yield strength of the aluminum alloy strip was 158 MPa, which was consistent with the experimental results (142.9 MPa).(
3
) The solidification process of twin-roll cast alloy was described accurately with the established mathematical model. The model can be used to predict the solidification structure and mechanical properties of the strip under different process parameters and solidification parameters.



---

*Source: 9304038-2017-10-02.xml*