Abstract
The influence of different cumulative rolling processes on the distribution of thermal parameters of sheet was simulated using DEFORM commercial finite element software. The dynamic globularization process of TA15 titanium alloy after thermal compression and cumulative rolling was simulated by the cellular automata. The result shows that during the accumulative roll bonding (ARB) process, the microstructure and properties of TA15 titanium alloy can be effectively improved by reducing the deformation of single pass and rolling speed as well as maintaining a certain temperature. A reasonable dynamic globularization cellular automata (CA) model was imported into Deform-3D commercial finite element simulation software and the microstructure evolution was simulated in the process of thermal compression and accumulative rolling-bonding.
Science Press
Titanium alloys possess high specific strength, thermal stability and corrosion resistance, and have been increasingly used as key structural parts in aerospace, automotive, and biomedical application
In the past 30 years, the severe plastic deformation (SPD) technology has been investigated over a wide range. The typical SPD technologies mainly include high pressure torsion (HPT), multi-directional compression/forging (MF), equal channel angular pressing (ECAP), torsion extrusion (TS), continuous shear deformation (CP) and accumulative roll bonding (ARB
In recent years, a lot of researches about the preparation process of ARB technology have been studied. Based on the theoretical derivation of equivalent strain, Xiang et a
Apart from experimental analyses, a powerful research tool for microstructure modeling, finite element model, has been developed, but it has less application in ARB simulatio
Hence, the ARB finite element model of TA15 alloy sheet was established in terms of the constitutive relation of thermo-mechanical coupling and accurate boundary conditions. Moreover, the distribution of thermal parameters under different process conditions was analyzed by numerical simulation, and the microstructure was predicted and verified by the cellular automata (CA
The material used was TA15 alloy with a composition (wt%) of 6.46Al-2.12Zr-1.12Mo-1.21V. It was received as the twice hot-rolled sheet with the size of 3000 mm×3000 mm×14.5 mm. The specimen was first heated to 1010 °C for 2 h, then cooled to 925 °C for 2 h, and finally water cooled to room temperature. The water cooling after heating was to ensure that the obtained microstructure at high temperatures is retained. The phase transformation temperature was measured at 995 °C by metallographic observation method in laser con-focal high-temperature scanning microscope (Lasertech Inc., Yonekura, Japan).

Fig.1 Microstructure of TA15 sheet after heat treatment
Cylindrical specimens with Φ8 mm×12 mm were cut from the as-received and heat treated material by wire-electrode cutting. The isothermal compression was carried out on a Gleeble-3800 machine (Dynamic Systems Inc., North Carolina, America) at 800, 850, 900, 930, and 960 °C, and strain rates of 0.001, 0.1, 1, and 10

The thermo-elasto-plastic coupling equations for ARB are established based on the virtual work principle and large deformation theory, with the following form:
(1) |
where is elastoplastic stiffness matrix, is geometric stiffness matrix, is node velocity, is change rate relative to temperature and equivalent plastic strain rate, is strain matrix, and is rate of external load change.
Coulomb friction assumed that the pressure of the contact surface is directly proportional to the friction, expressed by
(2) |
where is the coulomb friction coefficient, is the shear yield limit of material, and is the normal stress on the friction surface.
Arctangent friction model of relative sliding velocity can be expressed by:
(3) |
where is the arctangent friction factor, K is the shear yield strength of material, is the relative speed, is the absolute value of rolling speed and is tuning parameter, generally 1
The flow stress of the treated sheet was imported into the material library of Deform-3D software. In traditional rolling process, the friction is generally defined as shear friction when the sheet contacts with roller and push shee

The schematic diagram of the deformation area during accumulative rolling process is summarized in Fig.4. There-fore, it is difficult to investigate the deformation of the same section based on time because of different rolling direction (RD) positions. In general, 1/2 of the sum of biting angle and neutral angle can be selected as a reference. Table 2 lists the process parameters of ARB under different thickness reduction.

Table 2 Process parameters of ARB
Cellular automata (CA) model is a grid dynamic model with discrete space-time situation and local spatial interaction, which particularly embody practical causality. In recent years, it has become an important method to predict the evolution of structure through the flexible adaptation to space-time scale.
The dynamic globularization can be divided into two stages: fracture nucleation and growth of lamellar. In order to simplify the simulation process, the dynamic globularization model adopts the following assumptions: the dislocation density in matrix can conform to Kocks-Mecking equatio

Fig.5 CA model block diagram
The polycrystal plasticity FE model is established based on Hill's quadratic anisotropy yield functio
(4) |
where F, G, H, L, M and N are the anisotropic constants.
The dynamic recrystallization adopts the nucleation model proposed by Ding and Gu
(5) |
where and m are material constants, Qact is activation energy, is nucleation rate, is strain rate, is temperature, and is Stefan-Boltzmann constant.
In order to simplify the CA microstructure model constructed, this work adopts the following assumptions: (1) the rotation of crystal grains is ignored; (2) there is no texture in the deformed microstructure, and the influence of texture on dynamic recrystallization in the actual deformed structure is ignored; (3) only the dislocations, crystal grains, grain boundaries and second phase particles existing in the crystal are considered, and the defects such as vacancies, interstitial atoms, stacking faults, twins are ignored; (4) only the dislocation density changes caused by dynamic recovery are considered, and its impact on microstructure changes is ignored.
The region was divided into a quadrilateral cell space of 200 μm×200 μ
The distribution of thermal parameters from the normal direction (ND) to RD can directly reflect the deformation process of the accumulative rolling process. The influence rules of the process parameters on stress, strain, strain rate and temperature field were established using the accurate finite element model to simulate the cumulative rolling process of TA15, and thus the measures to improve the homogeneity of microstructure can be obtained in the cumulative rolling process of TA15.

The change of the equivalent stress of A1 section in ND direction is illustrated in Fig.6. It can be seen that the surface stress is lower than the central stress, and the central stress increases with the increase of deformation. However, as for the stress distribution of traditional rolling process, the surface stress is generally higher than the center stress. The heat generated by surface friction is difficult to be transmitted to the interior through the way of heat conduction which is due to poor thermal conductivity of 6.7 W/(m·K) for titanium alloy. Moreover, the high strain rate generated by the accumulated rolling process makes the surface generate adiabatic temperature rise, which will further lead to the rapid rise of temperature and the reduction of equivalent stress. With the increase of deformation, the contact length between rolling piece and roll keeps rising, and the effect of rising temperature is also expanded. For the central stress, the higher deformation corresponds to the larger deformation area. Therefore, the central stress reaches a higher value because of the larger rolling pressure and average unit pressure.

Fig.7 shows the change of the equivalent stress of the central node with the rolling time at different deformation conditions. In Fig.7, the signs of t1, t2, t3, t4 and t5 correspond to the rolling time of the cross section under five kinds of deformation. From Fig.7, the geometric parameters (length/width) of the rolling deformation zone and the rolling time of the section in the rolling deformation zone all increase with the increase of deformation. During the period of stable rolling in the deformation zone, the larger the deformation, the higher the equivalent stress. In addition, the equivalent stress becomes higher when it is closer to the narrowest part of the wedge-shaped deformation zone with the rolling processing. The equivalent stress fluctuates significantly in two sections of entering and leaving the deformation zone, which is deter-mined by two unstable stages of biting in and throwing out.

Fig.8 represents the variation of equivalent plastic strain along the direction of ND under different thickness reduction. The distribution of the total equivalent effect and the variation of strain in the direction of thickness all present a rising trend as the deformation increases. The growth rate curves of strain show the growth laws that the deformation is characterized by an uneven increase. Because the simulation result is obtained under the single pass deformation, the conclusion can be derived that the strain gradient in thickness direction is affected by the deformation passes. Therefore, excessive deformation of the rolling pass should be avoided in the design of cumulative rolling to reduce the strain gradient.


Fig.9 shows the distribution of the equivalent plastic strain in centre and surface changing with time under the conditions of different thickness reduction. It can be seen that the equivalent plastic strain of both surface layer and central layer increases with the rolling time. With the increase of thickness reduction, the difference of the strain between centre and surface enlarges step by step, and the growth of centre is larger than that of surface, which leads to the uneven deformation in the thickness direction.
Fig.10 gives the distribution of shear strain of the RD-ND surface along ND direction in the section. It is obvious that the shear strain displays the similar behavior with the increase of deformation.
It should also be noted that a significant difference of the direction of shear strain occurs in the two sides of the intersection (at a thickness of 6 mm). This particular distri-bution is closely related to the characteristics of cumulative rolling. According to previous studies, friction force is the driving force of cumulative rolling and it belongs to the surface shear stress. The effect of friction force is attributed to the excellent ductility of rolled piece surface, but there is a spread gradient between single sheet and multi-layer sheet in the thickness direction. Therefore, the transformation of shear strain exists in the intersection. As clearly observed from the influence of deformation on thermal parameters, small deformation with the single pass can effectively improve the non-uniformity of distribution.
The effect of rolling speed on equivalent stress along ND direction is presented in Fig.11. It shows that the equivalent stress increases along with the increase of rolling speed and the growth of centre is significantly higher than that of surface.
Fig.11 Equivalent stress along ND direction under different rolling speeds


According to the distribution of temperature field shown in Fig.12, high temperature makes the metal flow and softening can be observed. Furthermore, the strain rate from surface to centre presents the evolution law of decreasing first and then increasing, which is influenced by the temperature change (Fig.13). In terms of the thermo-mechanical coupling and constitutive law, TA15 alloy has a high sensitivity to strain rate and deformation temperature. While the stress in the surface produced with high strain rate is offset by the high temperature, which directly results in the stress to stop growing or even to decrease. According to the influence of rolling speed on thermal parameters, temperature distribution from the centre to the surface can be improved by lowering rolling speed.

Fig.14 represents the distribution of equivalent stress and equivalent strain along the ND direction at different rolling temperatures. The equivalent stress of the section decreases as the rolling temperature increases for the same deformation area. However, the distribution of strain exhibits an opposite trend, i.e. it is lower in center while higher in surface. In addition, the strain range increases simultaneously with the increase of rolling temperature and deformation of centre and surface.


It can be seen from Fig.15 that the temperature rise of the sheet shows a trend of being lower in center and being higher in surface, which also indicates that the temperature rise of the surface is mainly derived from the high strain rate caused by high-speed cumulative rolling, thus leading to the “skin effect” in the superficial area. From the influence law of rolling temperature on the thermal parameters, the medium rolling temperature can effectively improve the strain and temperature range of the accumulative roll-bonding sheets, so as to obtain the optimal structure property.
According to the comparison from the above results, the grain size of 9.2, 13.6, and 18.9 μm can be obtained before the hot deformation at 800, 850, and 900 °C, respectively. Based on the assumption and simplification from the dynamic globularization above, the change of microstructure of TA15 alloy under different conditions (temperature of 800, 850, and 900 °C; strain rate of 0.01, 0.1, 1, and 10
From the comparison between the simulation results and experiment results in Fig.16, the results of CA simulation have a high advantage in accuracy in terms of grain morphology and volume fraction of recrystallization. On condition that the macro data of rolling stress, strain and temperature field can be extracted, the simulation analysis of microstructure of rolled sheet can also be achieved. The microstructure under optical microscope (OM) is similar to that under CA simulation, in which the proportion of primary phase is less, and the overall deformation direction of matrix is perpendicular to compression reverse. Fig.16c shows that equiaxed grains with an average size of about 0.86 μm can be observed in the area with large deformation. This prediction model is reliable when we combine the microstructure of characteristic points in Fig.16b with the analysis of Fig.16a.
Finally, it can be seen from



Fig.17 Microstructure comparison under 930 °C/1
1) The thermal parameters and the microstructure evolution in ARB process of TA15 alloy are predicted using finite element simulation software of Deform-3D, as implemented by the combination of thermo-mechanical FEM and cellular automata model.
2) A small thickness reduction in single pass ARB can effectively improve the nonuniformity of shear strain distribution, and the shear strain displays the similar behavior with the increase of deformation. The temperature distribution from the centre to the surface can be improved by lowering rolling speed.
3) The microstructure and properties of TA15 alloys can
be effectively improved by reducing the thickness reduction (
Reference
Yi J, Wang X B, Jiao L et al. Thin-walled Structures[J], 2019,144: 106 329 [Baidu Scholar]
Xiang J F, Yi J. CIRP Journal of Manufacturing Science and Technology[J], 2021, 32: 356 [Baidu Scholar]
Fan X G, Yang H, Gao P F. Journal of Materials Science[J], 2011, 46: 6018 [Baidu Scholar]
Zhou Xiaolong, Tao Qiying, Zhou Yunhong. Rare Metal Materials and Engineering[J], 2017, 46: 942 [Baidu Scholar]
Xiang J, Pang S, Xie L J. The International Journal of Advanced Manufacturing Technology[J], 2018, 97: 1673 [Baidu Scholar]
Saito Y, Utsunomiya H, Tsuji N. Acta Materialia[J], 1999, 47: 579 [Baidu Scholar]
Xiang J, Pang S, Xie L J. The International Journal of Advanced Manufacturing Technology[J], 2018, 98: 1237 [Baidu Scholar]
Xing Z P, Kang S B, Kim H W. Metallurgical and Materials Transactions A[J], 2002, 33: 1521 [Baidu Scholar]
Kwan C, Wang Z, Kang S. Materials Science and Engineering A[J], 2008, 480: 148 [Baidu Scholar]
Tsuji N, Saito Y, Lee S H. Advanced Engineering Materials[J], 2003, 5: 338 [Baidu Scholar]
Lee S H, Saito Y, Tsuji N. Scripta Materialia[J], 2002, 46: 281 [Baidu Scholar]
Xiang J F, Xie L J, Gao F. Chinese Journal of Aeronautics[J], 2020, 34(8): 218 [Baidu Scholar]
Wang H, Lu C, Tieu K. Crystals[J], 2019, 9: 119 [Baidu Scholar]
Wang H, Su L, Yu H et al. Materials Science and Engineering A[J], 2018, 726: 93 [Baidu Scholar]
Khatami R, Fattah-Alhosseini A, Mazaheri Y. The International Journal of Advanced Manufacturing Technology[J], 2017, 93: 681 [Baidu Scholar]
Xiang J, Xie L, Gao F et al. Ceramics International[J], 2018, 44: 5341 [Baidu Scholar]
Yi J, Xiang J, Yi F et al. The International Journal of Advanced Manufacturing Technology[J], 2018, 97: 1673 [Baidu Scholar]
Xiang J, Xie L, Gao F et al. Ceramics International[J], 2019, 44: 117 65 [Baidu Scholar]
Xiang J, Pang S Q, Xie L J. Materials[J], 2018, 11: 252 [Baidu Scholar]
Mecking H, Kocks U F. Acta Metallurgica[J], 1981, 29: 1865 [Baidu Scholar]
Ding R, Guo Z X. Acta Materialia[J], 2001, 49: 3163 [Baidu Scholar]