1. Introduction
Extra-thick steel plates are primarily utilized in offshore engineering, the military industry, nuclear power, and high-rise buildings, among other critical technical equipment manufacturing fields [
1,
2]. It is a hot issue that the production of extra-thick steel plates using the ingot casting process results in a comparatively low metal yield. Given this situation, the continuous casting (CC) process is considered a more efficient method that contributes to energy saving and carbon reduction [
3]. However, the CC process has numerous influencing factors, and the accessibility of the appropriate process parameters increases with the thickness of the slab, so the extra-thick slabs are prone to internal defects such as central segregation and porosity shrinkage [
4,
5]. The solidification structure characteristics determine the formation of the internal defects of slabs, the adequate improvement of which is a pivotal factor for the production of ultra-thick steel plates [
6]. Furthermore, the wider equiaxed crystal zone in the center of the slab protects against the formation of mini-ingots, leading to severe V-shaped segregation and porosity shrinkage [
7]. To some extent, the fluctuation of the alloy elements content, a non-negligible effect on the final solidification structure, also causes variations in the solidus temperature, liquidus temperature, and dendrite growth parameters [
8]. Meanwhile, the tiny shrinkage holes in the slab are not wholly compressed in the subsequent rolling process, as the compaction ratio is limited. Also, hardly will the macroscopic chemical composition distribution maintain uniformity during the heat treatment process [
9]. Therefore, it is essential for the realization of high-quality extra-thick slabs to explore the influence of CC process parameters and alloying element content on the solidification structure characteristics.
To date, numerous studies have been conducted on the effect of CC process parameters (e.g., superheat, casting speed, and specific water flow) on the solidification structure of slabs across a range of thicknesses [
10,
11,
12]. The heat transfer rate in the primary and secondary cooling zones, the main factor influencing the local temperature gradient and solidification rate, is mainly determined by the CC process parameters during slab solidification. Additionally, Hunt et al. [
13] concluded that the criterion for the columnar to equiaxed transition (CET) depends on both the temperature gradient and the solidification rate, which determines the distribution characteristics of the central equiaxial and columnar crystal zones in the billet. Unlike columnar crystals, equiaxed crystals have no specific direction in the growth process. Thus, the residual steel enriched with solute elements can be distributed relatively uniformly in the mushy zone, alleviating the central segregation of the slabs [
14]. Zhang et al. [
15] developed a macro-segregation model and a solidification structure model for 300 mm thick slabs, and the results showed that a lower superheat can reduce grain size, increase ECR, and lessen central segregation. Furthermore, Chen et al. [
16] studied the microstructure of 400 mm extra-thick slabs by numerical simulation, and believed that the center ECR and the average grain size are little changed along with the change in casting speed and secondary cooling compared to superheat. In contrast to extra-thick slabs, however, the cooling intensity in the secondary cooling zone has a larger impact on the fluctuation of the solidification structure of the slab with common thickness. As for slabs with a thickness of 200 mm, Zhang et al. [
17] found that reducing the specific water flow is conducive to enhancing the ECR of 24Mn steel slabs through a numerical–physical coupled simulation method. Sheng et al. [
18] observed that the CET position near the corners of a 250 mm thick slab progressively moved towards the center of the slab and the ECR decreased due to the cooling intensity of the slab corners declining.
Also, the variations in the alloy element content have an essential influence on the thermophysical properties of the steel grade, which determines the solidification structure development of the slab with defined CC process parameters [
8,
19,
20]. Cai et al. [
21] concluded that the ECR initially increased and then decreased when the carbon (C) content was improved from 0.09% to 0.53%, and the ECR is at a peak at a C content of 0.3%. Moreover, the dendrite arm spacing is a fundamental feature parameter for characterizing the compactness degree of the solidification structure. Zhang et al. [
22], through a combination of numerical simulations and experiments, found that decreasing the silicon (Si) and C content can reduce the spacing of secondary dendrite arms (SDAS) in slabs with a thickness of 295 mm, promoting the improvement of central segregation. Likewise, Li et al. [
23] analyzed the relationships between Si, Cr, manganese (Mn) contents, and solidification structure based on the numerical simulations method for a 72 mm thickness slab. Yuan et al. [
24] discovered that higher Mn and C contents would suffer from coarse solidification structure and carbide, which made it difficult to find suitable parameters for a continuous casting process, using various experimental methods for Mn13 steel slab with a 230 mm thickness. Gao et al. [
19] developed the solidification structure model and SDAS model for a 370 mm extra-thick slab. It was found that the increase in Nb and V contents can refine the grain, and SDAS has enormous fluctuation with the elevation of Ti element contents. However, there are limited reports related to the influence of alloying element content on the solidification structure.
Analyzing the solidification structure evolution of billets through numerical simulation can save a huge amount of resources and time compared to traditional metallurgical analysis methods [
25]. Current simulation methods for solidification structures include stochastic methods, deterministic methods, and phase-field methods [
26]. Furthermore, the cellular automaton (CA) method is one of the stochastic methods, which has the physical background mechanism of dendrite tip growth kinetics and the characteristics of high simulation accuracy [
27]. Combining it with the finite element (FE) model of macroscopic heat transfer contributes to handling nonlinearity and complexity problems [
28]. Utilizing the cellular automaton–finite element (CAFÉ) model, Lu et al. [
29] determined the optimal combination for temperature distribution and ECR within the specified parameter range of thin slab by cross-scale and full-process calculations. A unique CAFÉ analysis approach was suggested by Tan et al. [
30] and is appropriate for the solidification structure simulation of 38CrMoAl large round bloom considering mold and final electromagnetic stirring. Nonetheless, the solidification structures simulated by the CAFÉ model are primarily billets, thin slabs, and round blooms, while there are relatively few simulations of ultra-thick slabs over 400 mm.
In addition, it is crucial to emphasize that the slab thickness is closely related to solidification structure characteristics. The bigger section size will give rise to slower solidification in the central region and longer local solidification times, which is more severe in comparison with the internal quality of the thin slabs [
31]. Zhang et al. [
32] compared the relationship between the solidification structure and central segregation of slabs with thicknesses ranging from 245 mm to 440 mm by numerical simulation. It was observed that the central segregation region and ECR are expanded with the increase in billet thickness. In addition, Xu et al. [
33] analyzed the microstructure evolution with 420 mm extra-thick slab under different casting speeds, superheats, and specific water flow, and the findings provide strategic guidance for optimizing the CC process and alleviating the central segregation of extra-thick slabs. Currently, the maximum thickness of extra-thick slabs produced by the vertical curved continuous caster is up to 475 mm [
34]. There is no doubt that the occurrence of internal defects in the slab obviously rises with the thickness of the casting slab. However, the study of the solidification behavior and solidification structure of extra-thick slabs with a thickness of 475 mm has not yet been reported.
In the present study, models for the solidification heat transfer and the cellular automaton–finite element of high-strength structural steel Q355 were established, with a section size of 475 mm × 2000 mm. The accuracy of these models was verified through nail penetration experiments, surface temperature measurements, and acid-etching tests. Moreover, the influence rules of various casting speeds, specific water flow, and superheats on the solidification behavior of extra-thick slabs were investigated in detail. Finally, the solidification structure characteristics of extra-thick slabs with different CC process parameters and alloy element content were calculated and analyzed. The outcomes from this study are expected to shed light on a better understanding of the changing characteristics of the solidification structure in 475 mm ultra-thick slabs, leading to a theoretical basis to achieve higher quality extra-thick steel plates.