Effect of Heat Input on Microstructure and Tensile Properties in Simulated CGHAZ of a V-Ti-N Microalloyed Weathering Steel

Abstract

The mechanical properties of a coarse-grained heat-affected zone (CGHAZ) are affected by welding thermal cycling with varied heat input (E_j), but its effect on tensile properties is rarely studied. In the present work, E_j = 15, 35, 55, 75 kJ/cm CGHAZ samples were prepared via Gleeble^TM (St. Paul, MN, USA) for a novel V-Ti-N microalloyed weathering steel. The tensile properties of CGHAZ with varied E_j were evaluated. The results indicated that mixed microstructures dominated by lath bainitic ferrite (LBF) and granular bainitic ferrite (GBF) were obtained at E_j = 15 and 35 kJ/cm, respectively, while a mixed microstructure composed of GBF, intragranular acicular ferrite (IGAF), and polygon ferrite (PF) formed at E_j = 55 and 75 kJ/cm, apart from martensite/austenite (M/A) constituents in each E_j condition. The above variation tendency in the microstructure with the increase in E_j led to coarsening of low-angle grain boundaries (LAGBs) and a decrease in dislocation density, which in turn resulted in a yield strength (YS) decrease from 480 MPa to 416 MPa. The mean equivalent diameter (MED), defined by the misorientation tolerance angles (MTAs) ranging from 2–6°, had the strongest contribution to YS due to their higher fitting coefficient of the Hall–Petch relationship. In addition, the increase in the average size (d_M/A) of M/A constituents from 0.98 μm to 1.81 μm and in their area fraction (f_M/A) from 3.11% to 4.42% enhanced the strain-hardening stress. The yield strength ratio (YR) reduced as the E_j increased, and the lower density and more uniform dislocations inside the ferrite led to a uniform elongation (uE) increase from 9.5% to 18.6%.

1. Introduction

Weathering steel (WS) is, essentially, an atmospheric-corrosion-resistant steel, usually containing no more than 0.20 wt % C, 0.25–0.55 wt % Si, and 0.8–1.5 wt % Mn, as well as one or more other alloy elements, such as Cu, Cr, Ni, with a total alloy content of 1.00–5.00 wt % [1,2,3]. WS normally possesses a stable rust layer to prevent erosion, having many advantages over laborious, costly, and high-emission coatings, and it has been increasingly used for important applications, such as uncoated bridges, buildings, and transmission towers, through the traditional arc welding process [4,5,6]. However, during the welding of WS, due to the influence of thermal cycling, the strength and toughness reduction in CGHAZ near the fusion line cannot be ignored [7,8]. Since weathering steel contains more Si [9], Cr [4,10], and other alloy elements, after subjection to the welding thermal cycle with a peak temperature higher than 1300 °C, hard brittle phase GBF and M/A constituents are enriched in CGHAZ [11,12]. When the E_j increases, the above-mentioned hard and brittle phases will coarsen and increase in quantity [13,14], further increase the reduction in toughness and strength.

Extensive investigations [15,16] have reported that the nucleation of heteronucleated ferrites on inclusions or precipitated particles can reduce the content of GBF, ferrite side plate, etc., refining the CGHAZ microstructure and improving its mechanical properties. V-Ti-N microalloyed steels with high ferritic heterogeneous nucleation capacity have received extensive attention. These utilize Ti (C, N) particles with high thermal stability to pin prior austenite grain boundaries (PAGBs) and refine them. VN and/or V (C, N) particles from the solid γ-phase reduce the mismatch with ferrite and promote the heterogeneous nucleation of ferrites [17,18]. These studies provide good microalloyed methods to improve the CGHAZ mechanical properties of weathering steel. Zhang and Hu et al. [19,20] studied the effect of E_j on the CGHAZ microstructure and impact properties of V-Ti-N microalloyed C-Mn steel; the results showed that with the increase in E_j, the ability of V-rich particles to promote heterogeneous nucleation of ferrite was enhanced, and the corresponding impact energy increased. Shi et al. [21] also obtained similar research conclusions: as the time of T_8/5 is increased from 45 s to 285 s, the CGHAZ toughness of V-Ti-N microalloyed steel significantly improved; however, the YS of CGHAZ decreased from 385 MPa to 335 MPa due to the formation of a large amount of intragranular polygon ferrite (IGPF) with the extension of T_8/5. Moreover, Wang et al. [22] investigated the CGHAZ tensile properties of V-Ti-N microalloyed steel containing 0.28 wt % Mo under varied E_j, during the process of increasing E_j from 35 kJ/cm to 65 kJ/cm and 120 kJ/cm; the YS of CGHAZ showed a trend of first increasing from 545 MPa to 610 MPa, and then decreasing to 394 MPa, respectively, the main reason being that the CGHAZ main microstructure transformed from GBF to IGAF and IGPF, respectively.

Considering that the influences of heat inputs on CGHAZ microstructure and the mechanical properties of V-Ti-N microalloyed steels with varied compositions are different, for the V-Ti-N microalloyed weathering steel, there may still be complex connections between the heat inputs and CGHAZ microstructure and mechanical properties. The effect of heat inputs on their CGHAZ impact property was reported by the author [23], because the change in CGHAZ strength can reflect its tendency towards softening and hardening, while also affecting the fatigue performance of the joint [24]. It is necessary to further evaluate the strength of CGHAZ with different E_j for V-Ti-N microalloyed weathering steel, and to the authors’ knowledge, research on heat inputs and CGHAZ strength has been scarcely reported previously.

In addition, during the continuous cooling process of the welding thermal cycle, bainitic ferrite, IGAF, IGPF, and other mixed microstructures are produced in the CGHAZ of V-Ti-N microalloyed weathering steel. It is difficult to establish relationship between complex microstructure and tensile properties [25,26]. In recent years, electron backscatter diffraction (EBSD) technology has been developed to analyze the influence of complex microstructure on tensile properties based on the variation of MTAs of different grain boundaries [25,27]. Fan et al. [28] showed that grain boundaries with MTAs ranging from 2–6° are the main boundaries controlling the tensile properties for low-alloy steel, while Zhu et al. [29] pointed out that the MED defined by MTAs smaller than 10° effectively governed the tensile properties. Olasolo et al. [30] indicated that MED defined by LAGBs less than 4° strongly contributes to tensile properties. Further, Li et al. [31] found the closest Hall–Petch relationship between the YS and the ferrite grains defined by LAGBs ranging from 2–12° in C-Mn steel. Hence, complex CGHAZ microstructure and corresponding tensile properties need further study. Moreover, the evolution of E_j will lead to differences in M/A constituents, dislocation structure, and density in CGHAZ microstructure, and its influences on strain-hardening behavior and elongation need to be elaborated.

In the present work, a novel V-Ti-N microalloyed weathering steel was designed and prepared. The steel was subjected to a CGHAZ welding thermal cycle with E_j = 15, 35, 55, 75 kJ/cm on a Gleeble-3500, and then the CGHAZ microstructure under each heat input condition was characterized and its tensile properties were tested. Moreover, the CGHAZ microstructure was further analyzed and quantitatively characterized by BESD and TEM, the relationship between the complex CGHAZ microstructure and tensile properties was evaluated, and the most effective structural unit for controlling YS was determined. In addition, the influence of CGHAZ microstructure on strain-hardening ability and elongation under different heat input conditions had also been studied. The research in this article will further clarify the influence of heat input on the microstructure and tensile properties of the welding heat-affected zone and understand the tensile fracture behavior with multiphase microstructure.

2. Materials and Experimental Method

The steel used in this article is a novel V-Ti-N microalloyed weathering steel. The specific chemical composition is shown in Table 1, and the mechanical properties of the test steel are shown in Table 2; this indicates that this novel weathering steel has the advantages of high strength, high toughness, and low yield ratio.

Round bars with sizes φ15 × 80 (rolling direction) mm² were machined from the test steel. The simulated CGHAZ with E_j = 15, 35, 55, and 75 kJ/cm were conducted in a Gleeble-3500 system (Figure 1a). The samples were heated to 1320 °C at 100 °C/s and cooled based on the thermal cycle curve developed via the Rykalin-2D mode; the model formula is shown in Equation (1), the T_8/5 time corresponding to E_j = 15, 35, 55, and 75 kJ/cm is 18.9 s, 50.6 s, 112.3 s, and 232.5 s, respectively, and the thermal cycle curves are shown in Figure 1b. Taking samples from a 2 mm area on both sides of the thermocouple of the round bars, the CGHAZ microstructure and tensile properties of each E_j sample was detected, as shown in Figure 1c. Tensile properties are tested on an MTS servo hydraulic testing machine at room temperature with a strain of 0.25 mm/min, and the specimens obtained for each E_j were tested three times. The YS was determined by the 0.2% offset flow stress. In addition, a Vickers hardness tester (HV-10Z) was used to test the hardness of different heat input samples, and the average value was measured 5 times under each condition.

The microstructure was characterized via optical microscopy (OM; Axiover-200MA T, ZEISS, Jena, Germany), scanning electron microscopy (SEM; Hitachi-SU5000, Hitachi, Tokyo, Japan), and transmission electron microscopy (TEM; Talos-F200X, FEI, Hillsboro, USA) for each E_j sample. OM and SEM specimens were polished and etched with 4% nitrate alcohol. TEM samples for observing precipitated particles at different E_j were prepared by carbon extraction replication. The film samples of TEM were obtained by sanding to 30 μm and subsequently electrolytic polishing, and the polishing occurred in a solution of 5% perchloric acid in ethyl alcohol with v = 25 V, i = 55–65 mA, and T = 25 °C. The OM samples were re-polished and then etched with LePera’s reagent. The M/A constituents were corroded into bright white and subsequently observed by OM, then Image-Pro Plus (IPP) software (Media Cybernetics, Austin, USA) was used to evaluate the size and area fraction of M/A constituents for each E_j through more than 1000 M/A constituents in 15 fields of view. Samples from each E_j were ground and then electrolytically polished in 90% methanol solution of perchlorate with v = 18 V, I = 1 A, and t = 30 s, and evaluated by SEM equipped with a TSL electron backscatter diffraction (EBSD); the scanning step was 0.20 μm, and the MED of grains at different MTAs was determined by EBSD. In addition, an XRD (Rigaku D/max-2500 PC, Rigaku, Tokyo, Japan) with Cu K-α radiation was employed to determine the dislocation density for each E_j sample, and the specimens were scanned at a step width of 0.02° and counting time of 2 s over the range of 40°–105°, with three tests per sample.

(1)$E_j=\sqrt{\frac{4\mathsf{\pi }kc\rho \Delta \mathrm{t}}{1/{\left(T_2 -T_0\right)}^2 - 1/{\left(T_1 - T_0\right)}^2}}\times \mathsf{\delta }$

where E_j is heat input; c is the specific heat capacity, 0.586 J/(g·°C); ρ is the density of iron, 7.800 g/cm³; k is the heat conductivity coefficient, 0.360 J/(cm·s·°C); δ is the thickness of plating, 0.24 cm; T₁ and T₂ are the defined cooling time temperatures T₁ = 800 °C, T₂ = 500 °C, respectively; T₀ is the preheating temperature, 25 °C; Δ_t is temperature cooling time from T₁ to T₂.

3. Experimental Result

3.1. Microstructural Characteristics

OM observations of simulated CGHAZ samples under different E_j are displayed in Figure 2; there are obvious differences in the morphology of the microstructures. When the E_j = 15 and 35 kJ/cm, the prior austenite grain boundaries (PAGBs) can be observed, and the PAGBs may be further divided by lath bainitic ferrite (LBF) and granular bainitic ferrite (GBF). However, the proportion of LBF is higher at E_j = 15 kJ/cm, while a microstructure of GBF is dominant at E_j = 35 kJ/cm. As the E_j further increased to 55 and 75 kJ/cm with two large heat inputs, PAGBs were almost covered by growth of ferrite and the CGHAZ microstructures were mainly intragranular acicular ferrite (IGAF) and polygon ferrite (PF). These types of microstructures occupy a higher proportion at E_j = 75 kJ/cm. In addition, developed GBF can also be found in CGHAZ microstructure of two large heat inputs.

Moreover, Figure 3 characterizes the variation in M/A constituents of simulated CGHAZ under different E_j; the M/A constituents appear bright white after being corroded by LePera’s reagent. With the increase in E_j, the macroscopic morphology of M/A constituents gradually transforms from elongated to blocky. Meanwhile, according to the statistics of the IPP software, the average size (d_M/A) and area fraction (f_M/A) of M/A constituents increase with the increase in E_j. Specifically, d_M/A increases from 0.98 μm to 1.13 μm, 1.68 μm, and 1.81 μm, and f_M/A increases from 3.11% to 3.56%, 4.21%, and 4.42%, respectively.

The effect of E_j on CGHAZ microstructure is further characterized by TEM. With the increase in E_j from 15 kJ/cm to 75 kJ/cm, the ferrite morphology changes from parallel lath bundles to blocky, and the size of ferrite coarsens significantly, with the results shown in Figure 4a–d. The black island structures can be found in each E_j sample, and their size increases with the increase in E_j. The island structure is identified as M/A constituent by light and dark fields and selected area electron diffraction (SAED) pattern (Figure 4e–e2). Moreover, the dislocation structures inside the ferrite are also characterized; when the E_j = 15 and 35 kJ/cm, there are high-density dislocation structures in the ferrite. In the partial enlarged view of Figure 4a,b, dislocation tangles (DTs) occur around the elongated M/A constituents (Figure 4a1,b1). With a further increase in E_j to 55 kJ/cm, dislocation density decreases significantly compared to E_j = 15 and 35 kJ/cm, but compared to the relatively uniformly distributed dislocation lines (UDDLs) at E_j = 75 kJ/cm (Figure 4d1), the internal dislocation lines (DLs) are densely distributed (Figure 4c1).

In addition, XRD patterns of samples that are subjected to different E_j are displayed in Figure 5a. The dislocation density can be calculated by Equation (2) according to Ref. [32]; the relevant calculation results are shown in Figure 5b. With the increase in E_j from 15 kJ/cm to 75 kJ/cm, the dislocation density decreases from 2.20 × 10¹⁴ m⁻² to 0.73 × 10¹⁴ m⁻².

(2)$\rho =k{\epsilon }^2/{Fb}^2$

where ρ is the dislocation density, k is the Burgers vector of dislocations, ɛ is microstrain, (b) along <111> for body-centered cubic metals; the value of F is assumed to be 1.

Moreover, TEM characterizations of CGHAZ carbon extraction replica samples are shown in Figure 6. When the E_j = 15 and 35 kJ/cm (Figure 6a,b), the precipitated particles are mainly square particles. When the E_j increases to 55 and 75 kJ/cm (Figure 6c,d), in addition to square particles, there are also round particles in the microstructure, and the density and size of precipitated particles is enhanced with the increase in E_j. Further, the typical precipitated particles under each E_j condition are analyzed by energy dispersive spectroscopy (EDS), as shown in Figure 6I–IV. The EDS results indicate that the square particles are Ti (C, N) particles at low E_j = 15 and 35 kJ/cm (Figure 6I,II). However, when the E_j increases to 55 and 75 kJ/cm, based on the results of EDS analysis, the circular precipitated particles appearing in the sample are V (C, N) particles (Figure 6III), and the chemical composition of the square particles transforms into V-rich (Ti, V) (C, N) particles (Figure 6IV); this is also consistent with the research results of Zhang and Hu et al. [19,33]. In addition, 20 square particles and round particles under each heat input condition are selected respectively, and their EDS data were further processed and the atomic fraction of elements were statistically analyzed; the results are shown in Table 3. The results indicate that with the increase in E_j, the atomic fraction of V in square and circular precipitates is increased.

The EBSD inverse pole figure maps and the corresponding kernel average misorientation (KAM) maps of the welding thermal cycle at E_j = 15–75 kJ/cm are shown in Figure 7a–d and Figure 7e–h, respectively. The high-angle grain boundaries (HAGBs) are defined by misorientation tolerance angles (MTAs) of boundaries higher than 15° and marked by black lines, while the low-angle grain boundaries (LAGBs) are marked by white lines, which are defined by MTAs of boundaries ranging from 2°–15° [34]. The mean equivalent diameter (MED) of grains with varied MTA is evaluated, with the results counted in Table 4. With the increase in E_j, coarsening occurs in the MED defined by MTAs of 2°, 4°, 6°, 8°, 10° and 12°. Moreover, the KAM maps can be employed to estimate the degree of plastic deformation and the density of geometrically necessary dislocations (GNDs) [35,36]. The average KAM value (KAM_ave) decreases from 0.83° to 0.55° with the increase in E_j; this trend is also consistent with the results of TEM and XRD observations, where an increase in E_j leads to a decrease in dislocation density.

3.2. Mechanical Properties of CGHAZ

The tensile properties and Vickers hardness of CGHAZ samples with different welding heat inputs were tested, and the relevant results are shown in Table 5, and the typical CGHAZ stress-strain curves and detailed tensile properties of each E_j sample are shown in Figure 8a and Figure 8b, respectively. With the increase in E_j, the hardness of the CGHAZ reduced from 238.4 HV10 to 201HV10, and the YS and TS decreased from 480 MPa to 416 MPa and from 615 MPa to 574 MPa, respectively, while the YR decreased monotonically from 0.78 to 0.72. In addition, the uE increased significantly from 9.5% to 20.1%.

4. Discussion

4.1. Effect of E_j on CGHAZ Microstructure

With the increase in E_j from 15 kJ/cm to 35 kJ/cm and 55/75 kJ/cm, the main CGHAZ microstructure changes from LBF to GBF, and further changes to a mixture of PF and IGAF. The evolutions in microstructure are related to the difference in the cooling process caused by varied E_j. When the E_j = 15 kJ/cm to 35 kJ/cm, the high-temperature residence time is shorter and the cooling rate is faster, which inhibits the diffusion of C atoms, and the supercooling degree is larger and hence promotes the transformation of bainitic ferrite [37]. The PAGBs are known to be the nucleation sites of bainitic ferrite due to the irregular atomic arrangement and lattice distortion [11]. The γ (austenite)→LBF transition tends to require a greater degree of supercooling than the γ→GBF transition [8]; therefore, with the increase in E_j, the nucleation of LBF is inhibited, resulting in a large amount of GBF at E_j = 35 kJ/cm, while the microstructure is mainly LBF at E_j = 15 kJ/cm.

In addition, it is noteworthy that uniform nucleation of ferrite and heterogeneous nucleation of ferrite are very difficult to detect under small heat inputs, which is related to limited carbon diffusion and high undercooling in favor of bainite ferrite transformation, as well as to heteronuclear particles. Although the nucleation temperatures of GBF and IGAF are both between 400–600 °C according to Ref. [38], at small heat inputs, the precipitated particles that reach the critical nucleation size (particles with an equivalent circular diameter higher than 100 nm have nucleation ability according to Refs. [38,39]) are Ti (C, N) particles (Figure 6a,b). Previous studies [40,41] have reported that the ability of Ti (C, N) particles as heterogeneous nucleation sites is weak; as shown in Figure 9a, Ti (C, N) particles fail to promote heterogeneous nucleation of ferrite. The competitive nucleation and growth of GBF and IGAF under small heat inputs are mainly dominated by GBF, while IGAF is inhibited. Furthermore, increasing the E_j to 55/75 kJ/cm, the high-temperature residence time is extended, the corresponding undercooling decreases, and the diffusion of carbon atoms is promoted, all of which are beneficial for γ→α (ferrite) phase transition to occur at a higher temperature. Meanwhile, the precipitated particles transform from Ti (C, N) to V-rich (Ti, V) (C, N) particles (Figure 6c,d). Compared with Ti (C, N) particles, the mismatch between V-rich (Ti, V) (C,N) particles and ferrite is lower, which can significantly promote the heterogeneous nucleation of ferrite [42], as shown in Figure 9b; the heterogeneous nucleation of the ferrite on the V-rich (Ti, V) (C, N) particles is evident. Accordingly, the nucleation at PAGBs is mainly GBPF, while the heterogeneous nucleation of IGAF and IGPF occurs in the grain; thereby high proportions of PF and IGAF are formed in the CGHAZ when E_j = 55 and 75 kJ/cm.

In fact, the γ→α phase transformation is accompanied by the diffusion of C atoms from α to γ due to the higher solubility of C atoms at γ, and the C-rich retained austenite (γ’) forms during this period. γ’ can be stable at relatively high temperature because of the high C content, and then transform into M/A constituents in the air cooling process [43]. With the increase in E_j from 15 kJ/cm to 75 kJ/cm, the high-temperature residence time gradually prolongs, and the diffusion of C becomes more sufficient. Therefore, with the increase in E_j, the mean size and area fraction of M/A constituents increase.

4.2. Effect of E_j on CGHAZ Tensile Properties

With the increase in E_j, the grain size, the dislocation density, and the size and volume fraction of precipitated particles in CGHAZ will all change, and these will have an effect on the strength. According to previous research [28,44], the YS of low-alloy steel is generally determined by the cumulative contributions of various strengthening mechanisms, and the main strengthening modes are shown in Equation (3):

(3)${\mathsf{\sigma }}_{\mathrm{y}}={\mathsf{\sigma }}_0+{\sigma }_{ss}+{\sigma }_d+{\sigma }_{\rho }+{\sigma }_{pp}+{\sigma }_{M/A}$

Here, σ_y is the YS, σ₀ represents the Peierls–Nabarro stress or lattice friction stress, σ_ss represents the solid solution strengthening, σ_d represents the boundary strengthening, σ_ρ represents the dislocation strengthening, σ_pp represents the precipitation strengthening, and σ_M/A represents the M/A constituent strengthening. Among these strengthening mechanisms, the boundary strengthening exerted by LAGBs and the dislocation strengthening exerted by dislocation structures in the microstructure contribute the greatest [28,45]. With the increase in E_j, the MED defined by LAGBs are coarsened (Figure 7 and Table 4), while the dislocation density decreases (Figure 4 and Figure 5), resulting in a monotonic decrease in YS of CGHAZ.

Typically, the LAGBs defined by MTAs ranging from 2° to 15° have the most effective contribution to boundary strengthening [22], the HAGBs (θ ≥ 15°) have the most significant contribution to blocking crack propagation [21,46], and the remaining boundaries (θ < 2°) contribute to dislocation strengthening [47]. YS and grain size can be described by the Hall–Petch equation, as shown in Equation (4).

(4)${\sigma }_y={\sigma }_0*+{ \mathrm{k}}_{\mathrm{hp}}{d}^{-1/2}$

Here, σ_y is the YS, σ₀* is the is the sum of other strengthening methods except the boundary strengthening, k_hp is the structural constant, and d is the grain size corresponding to MED defined by varied MTAs. It is worth noting that the microstructure is complex; under different E_j conditions, there are obvious differences in the size, morphology and distribution of CGHAZ microstructures. In this paper, with the change in E_j, the YS have been plotted as a function of the MED^−1/2 with MTAs ranging from 2° to 12°, with the results shown in Figure 10. The results indicate that the MED_{MTA≥2–6°}^−1/2 defined by MTAs ranging from 2° to 6° showed a fairly strong linear fit with YS, with the correlation factor of 0.99–0.96 (Figure 10a–c). However, Figure 10 d–f show the correlation factor was reduced to 0.82–0.76 for the YS and MED_{MTA≥8–12°} defined by MTAs ranging from 8° to 12°. Accordingly, the MED^−1/2 defined by MTAs ranging from 2° to 6° had a higher correlation factor compared with MED^−1/2 defined by MTAs ranging from 8° to 12°. This result made it clear that that MED defined by MTA ranging from 2° to 6°of LAGBs is the most efficient microstructural unit to controlling YS; this result is also consistent with that reported in related literature [28].

Furthermore, with the increase in E_j from 15 kJ/cm to 75 kJ/cm, the YR monotonically decreases, which means that an increase in E_j leads to an increase in strain-hardening ability. This phenomenon is mainly related to the area fraction and size of the M/A constituents. M/A constituents exist in ferrite matrixes as hard phases, and their hardnesses are higher than those of ferrite matrixes, and M/A constituents can be plastically deformed by two-stage or three-stage work hardening. This deformation mechanism is mainly controlled by high-strength local plastic deformation; the local internal stress applied from the M-A constituent hard phase acts on the ductile ferrite phase, thereby improving the strain-hardening ability [48,49]. Accordingly, with the increase in E_j from 15 to 75 kJ/cm, the d_M/A increases from 0.98 μm to 1.81 μm, and f_M/A increases from 3.11% to 4.42%, leading to an increase in the strain-hardening ability, and the YR monotonically decreases.

Furthermore, the uE increases significantly with the increase in E_j from 15 to 75 kJ/cm. The uE is closely related to slip of dislocation: once dislocations slip over long distances without encountering barriers, high uE is generated [50,51]. When the E_j = 15 and 35 kJ/cm, there are many DTs and ferrite lath boundaries in the microstructure, and the dislocation density is high (Figure 4a,b and Figure 7e,f). The dislocations will be significantly hindered after short-distance movement during the stretching process. While E_j = 55 and 75 kJ/cm, the proportion of PF in microstructure increases, relatively uniform dislocations are distributed in the microstructure, and the dislocation density is low (Figure 4c,d and Figure 7g,h), these are beneficial for long-distance slip of dislocations and obtaining higher uE. In addition, the corresponding characterization results of tensile fracture are shown in Figure 11. The tensile fractures under low E_j and large E_j are composed of two parts: fiber zone and shear-lip zone: the fiber zone is mainly composed of dimples, but deeper and smaller dimples appear at E_j = 75 kJ/cm. Further, although microvoids can be found in all fiber regions, the number and the size of microvoids are larger at E_j = 15 kJ/cm. Compared to deeper dimples, the appearance of large-sized microvoids indicates that microcracks are more likely to occur during the stretching process [52], which is not conducive to obtaining high uE. Accordingly, the characterizations of tensile fracture further demonstrate that increasing E_j can improve uE.

5. Conclusions

The effects of heat input on CGHAZ microstructure and tensile properties were studied for a novel V-N-Ti microalloyed weathering steel, and the following main conclusions can be drawn:

• With the increase in E_j from 15 kJ/cm to 35 kJ/cm, and 55/75 kJ/cm, the dominant CGHAZ microstructure transforms from LBF to GBF, and then to a mixture of IGAF+PF. The size and area fraction of M/A constituents increase monotonically with the increase in E_j, the d_M/A increases from 0.98 μm to 1.81 μm, and the f_M/A increases from 3.11% to 4.42%.

• With the increase in E_j, the increase in MED defined by the LAGBs and a reduction in dislocation density led to the decrease in YS. And the MED with MTAs ranging from 2°–6° strongly controlled the YS due to their higher fitting coefficient of the Hall–Petch relationship.

• With the increase in E_j, the decrease in YR was mainly due to the increase in the size and area fraction of M/A constituents, which improved the strain-hardening ability. Meanwhile, the lower density and more uniform dislocation structures improved the uE.

Author Contributions

Q.W. (Qingfeng Wang) and B.H. conceived and designed the experiments; Q.W. (Qiuming Wang) and B.H. performed the experiments and analyzed the data; Q.W. (Qingfeng Wang) contributed the experimental material and analysis tools; and Q.W. (Qingfeng Wang) and B.H. wrote the paper. All authors have read and agreed to the published version of the manuscript.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

Acronyms	Full Names
WS	Weathering steel
CGHAZ	Coarse-grained heat-affected zone
Ej	Heat input
YS	Yield strength
TS	Tensile strength
YR	Yield strength ratio
uE	Uniform elongation
T 8/5	The time of cooling from 800 °C to 500 °C during the welding thermal cycle cooling stage
LAGBs	Low-angle grain boundaries
HAGBs	High-angle grain boundaries
MED	Mean equivalent diameter
MTAs	Misorientation tolerance angles
IPP software	Image-Pro Plus software
EBSD	Electron backscatter diffraction
d M/A	The average size of martensite/austenite constituents
f M/A	Area fraction of martensite/austenite constituents
SAED	Selected area electron diffraction
LBF	Lath bainitic ferrite
GBF	Granular bainitic ferrite
PF	Polygon ferrite
IGAF	Intragranular acicular ferrite
IGPF	Intragranular polygon ferrite
M/A constituents	Martensite/austenite constituents
PAGBs	Prior austenite grain boundaries
DTs	Dislocation tangles
UDDLs	Uniformly distributed dislocation lines
DLs	Dislocation lines
GNDs	Geometrically necessary dislocations
KAM	Kernel average misorientation

Footnotes

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Figures and Tables

Figure 1. Welding thermal simulation diagram of the round bars: (a), the welding thermal cycle curves with different heat input (b) [23], microstructure observation and micro-tensile specimens taken from the thermocouple region (c).

Figure 2. Optical micrographs of the simulated CGHAZ at Ej = 15 kJ/cm (a), 35 kJ/cm (b), 55 kJ/cm (c), and 75 kJ/cm (d). PF—intragranular polygon ferrite (IGPF) and/or grain boundary polygonal ferrite (GBPF).

Figure 3. The observation of M/A constituents in CGHAZ after LePera’s agent corrosion, Ej = 15 kJ/cm (a), 35 kJ/cm (b), 55 kJ/cm (c), and 75 kJ/cm (d).

Figure 4. Observation of CGHAZ ferrite and M/A constituents via TEM at Ej = 15 kJ/cm (a,a1), 35 kJ/cm (b,b1), 55 kJ/cm (c,c1), and 75 kJ/cm (d,d1). The bright field (e) and dark field (e1) of M/A constituent at Ej = 15 kJ/cm and SAED pattern of M/A constituent (e2).

Figure 5. XRD patterns (a) and dislocation density of specimens after heat cycle with different Ej (b).

Figure 6. TEM observation results of precipitated particles at Ej = 15 kJ/cm (a), 35 kJ/cm (b), 55 kJ/cm (c), 75 kJ/cm (d), the inset images (I–IV) in (a–d) show the TEM-EDS of particles.

Figure 7. The inverse pole figure maps (a–d) and the corresponding KAM maps (e–h) of CGHAZ, Ej = 15 kJ/cm (a,e), 35 kJ/cm (b,f), 55 kJ/cm (c,g), and 75 kJ/cm (d,h), and the inset images in (e–h) show the KAMave value of each Ej sample.

Figure 8. The typical stress-strain curves for each Ej sample (a) and statistical chart of CGHAZ tensile performance testing results (b).

Figure 9. Observation of precipitated particles in CGHAZ microstructures at Ej = 35 kJ/cm (a) and the corresponding EDS analysis of particles (I); observation of precipitated particles in CGHAZ microstructures at Ej = 75 kJ/cm (b), and the corresponding EDS analysis of particles (II).

Figure 10. The YS as a function of the MED−1/2 that defined by MTA of 2° (a); the YS as a function of the MED−1/2 that defined by MTA of 4° (b); the YS as a function of the MED−1/2 that defined by MTA of 6° (c), the YS as a function of the MED−1/2 that defined by MTA of 8° (d), the YS as a function of the MED−1/2 that defined by MTA of 10° (e), and the YS as a function of the MED−1/2 that defined by MTA of 12° (f).

Figure 11. Typical observations of tensile fracture morphology at Ej = 15 kJ/cm (a–a2) and Ej = 75 kJ/cm (b–b2).

Table 1

Chemical composition of experimental steel (wt. %) [23].

C	Si	Mn	P	S	Cr	Cu	Ni	Mo	Ti	V	N	Fe
0.046	0.19	1.30	0.009	0.005	0.45	0.32	0.37	0.05	0.015	0.024	0.0073	Bal.

Table 2

Mechanical properties of experimental steel [23].

YSReL/MPa	TSRm/MPa	YRReL/Rm	ElongationA/%	Energy Absorbed at−40 °C AkV/J	Micro-Hardness/HV10
489	594	0.82	26.5	238.4	226.0

Table 3

The atomic fraction of each chemical element in particles (wt %).

Ej kJ/cm	Type of Element in Square Particles					Type of Element in Circular Particles
	Ti	V	C	N	Fe	V	C	N	Fe
15	23.15%	0%	41.81%	24.15%	10.89%	/	/	/	/
35	23.19%	0.05%	42.43%	24.89%	10.44%	0.02%	60.15%	5.89%	33.94%
55	24.55%	0.95%	43.05%	26.15%	5.30%	0.88%	61.89%	10.25%	26.98%
75	25.11%	1.60%	43.10%	27.33%	2.86%	1.45%	62.82%	12.63%	23.10%

Table 4

MED of grains with boundaries at MTAs from 2° to 15°.

EjkJ/cm	MEDMTA≥2°/μm	MEDMTA≥4°/μm	MEDMTA≥6°/μm	MEDMTA≥8°/μm	MEDMTA≥10°/μm	MEDMTA≥12°/μm	MEDMTA≥15°/μm
15	2.75	3.15	3.45	4.15	5.21	5.47	6.36
35	3.16	3.6	4.13	4.96	5.62	5.76	6.73
55	4.03	4.54	5.12	5.67	5.79	5.84	6.15
75	4.69	5.56	5.68	5.76	5.87	5.91	5.96

MED—Mean equivalent diameter; MTAs—Misorientation tolerance angles; MED_MTA≥θ—The MED of grains with MTAs higher than θ.

Table 5

Summary of the tensile results at different E_j.

Heat Input/kJ·cm−1	YS/MPa	TS/MPa	YR	uE/%	Micro-Hardness/HV10
15	480 ± 5	615 ± 4	0.78	9.5 ± 0.5	238.4 ± 3.5
35	461 ± 6	605 ± 3	0.76	13.2 ± 0.7	222 ± 5.2
55	437 ± 6	585 ± 3	0.75	16.4 ± 0.7	210 ± 4.8
75	416 ± 5	574 ± 4	0.72	18.6 ± 0.6	201 ± 3.5

References

1. Wang, Z.; Zhang, X.; Cheng, L. Role of inclusion and microstructure on corrosion initiation and propagation of weathering steels in marine environment. J. Mater. Res. Technol.; 2021; 10, pp. 306-321. [DOI: https://dx.doi.org/10.1016/j.jmrt.2020.11.096]

2. Travassos, S.; Almeida, M.B.D.; Tomachuk, C. Non-destructive thickness measurement as a tool to evaluate the evolution of patina layer formed on weathering steel exposed to the atmosphere. J. Mater. Res. Technol.; 2020; 9, pp. 687-699. [DOI: https://dx.doi.org/10.1016/j.jmrt.2019.11.010]

3. Yamashita, M.; Nagano, H.; Misawa, T. Structure of protective rust layers formed on weathering steels by long-term exposure in the industrial atmospheres of Japan and North America. ISIJ Int.; 1998; 38, pp. 285-290. [DOI: https://dx.doi.org/10.2355/isijinternational.38.285]

4. Wang, D.; Zhong, Q.; Yang, J. Effects of Cr and Ni on the microstructure and corrosion resistance of high-strength low alloy steel. J. Mater. Res. Technol.; 2023; 23, pp. 36-52. [DOI: https://dx.doi.org/10.1016/j.jmrt.2022.12.191]

5. Morcillo, M.; Diaz, I.; Chico, B.; Cano, H. Weathering steels: From empirical development to scientific design. A review. Corros. Sci.; 2014; 83, pp. 6-31. [DOI: https://dx.doi.org/10.1016/j.corsci.2014.03.006]

6. Prasad, S.N.; Mediratta, S.R.; Sarma, D.S. Influence of austenitisation temperature on the structure and properties of weather resistant steels. Mater. Sci. Eng. A; 2003; 358, pp. 288-297. [DOI: https://dx.doi.org/10.1016/S0921-5093(03)00303-4]

7. Amer, A.E.; Koo, M.Y.; Lee, K.H. Effect of welding heat input on microstructure and mechanical properties of simulated HAZ in Cu containing microalloyed steel. J. Mater. Sci.; 2010; 45, pp. 1248-1254. [DOI: https://dx.doi.org/10.1007/s10853-009-4074-7]

8. Yang, X.; Di, X.; Liu, X.; Wang, D.; Li, C. Effects of heat input on microstructure and fracture toughness of simulated coarse-grained heat affected zone for HSLA steels. Mater. Charact.; 2019; 155, 109818. [DOI: https://dx.doi.org/10.1016/j.matchar.2019.109818]

9. Cao, R.; Chan, Z.S.; Yuan, J.J. The effects of Silicon and Copper on microstructures, tensile and Charpy properties of weld metals by refined X120 wire. Mater. Sci. Eng. A; 2018; 718, pp. 350-362. [DOI: https://dx.doi.org/10.1016/j.msea.2018.01.080]

10. Jorge, J.; Souza, L.; Rebello, J. The effect of chromium on the microstructure/toughness relationship of C–Mn weld metal deposits. Mater. Charact.; 2001; 47, pp. 195-205. [DOI: https://dx.doi.org/10.1016/S1044-5803(01)00168-1]

11. Zhou, P.; Wang, B.; Liang, W. Effect of welding heat input on grain boundary evolution and toughness properties in CGHAZ of X90 pipeline steel. Mater. Sci. Eng. A; 2018; 722, pp. 112-121. [DOI: https://dx.doi.org/10.1016/j.msea.2018.03.029]

12. Bao, L.; Huang, Y.; Liu, F. Investigation on microstructure and impact toughness of simulated heat affected zone of high strength low alloy steels by laser-arc hybrid welding. Mater. Res. Express; 2022; 9, 096519. [DOI: https://dx.doi.org/10.1088/2053-1591/ac93e9]

13. Wang, J.; Shen, Y.F.; Xue, W.Y. The significant impact of introducing nanosize precipitates and decreased effective grain size on retention of high toughness of simulated heat affected zone (HAZ). Mater. Sci. Eng. A; 2021; 803, 140484. [DOI: https://dx.doi.org/10.1016/j.msea.2020.140484]

14. Zhu, Z.; Liu, Y.; Gou, G. Effect of heat input on interfacial characterization of the butter joint of hot-rolling CP-Ti/Q235 bimetallic sheets by Laser+ CMT. Sci. Rep.; 2021; 11, 10020. [DOI: https://dx.doi.org/10.1038/s41598-021-89343-9]

15. Jia, X.; Chen, Y.; Li, H. Study on the mechanism of AF nucleation induced by complex oxide inclusions after LF refining in oxide metallurgical steel. Mater. Charact.; 2023; 204, 113239. [DOI: https://dx.doi.org/10.1016/j.matchar.2023.113239]

16. Qi, H.; Pang, Q.; Li, W. In Situ observation of CGHAZ microstructure evolution of hot-rolled ship plate steel containing 0.012 yttrium under high-heat-input welding. Steel Res. Int.; 2023; 98, 2200771. [DOI: https://dx.doi.org/10.1002/srin.202200771]

17. Fadel, A.; Glisic, D.; Radovic, N.; Drobnjak, D. Intragranular ferrite morphologies in medium carbon vanadium-micro alloyed steel. J. Min. Metall. B; 2013; 49, pp. 237-244. [DOI: https://dx.doi.org/10.2298/JMMB120820001F]

18. Medina, S.F.; Gómez, M.; Rancel, L. Grain refinement by intragranular nucleation of ferrite in a high nitrogen content vanadium microalloyed steel. Scr. Mater.; 2008; 58, pp. 1110-1113. [DOI: https://dx.doi.org/10.1016/j.scriptamat.2008.02.004]

19. Zhang, J.; Xin, W.; Luo, G. Effect of welding heat input on microstructural evolution, precipitation behavior and resultant properties of the simulated CGHAZ in high-N V-alloyed steel. Mater. Charact.; 2020; 162, 110201. [DOI: https://dx.doi.org/10.1016/j.matchar.2020.110201]

20. Hu, J.; Du, L.-X.; Wang, J.-J.; Xie, H.; Gao, C.-R.; Misra, R.D.K. High toughness in the intercritically reheated coarse-grained (ICRCG) heat-affected zone (HAZ) of low carbon microalloyed steel. Mater. Sci. Eng. A; 2014; 590, pp. 323-328. [DOI: https://dx.doi.org/10.1016/j.msea.2013.10.062]

21. Shi, Z.; Yang, C.; Wang, R.; Su, H.; Chai, F. Effect of nitrogen on the microstructures and mechanical properties in simulated CGHAZ of vanadium microalloyed steel varied with different heat inputs. Mater. Sci. Eng. A; 2016; 651, 184. [DOI: https://dx.doi.org/10.1016/j.msea.2015.10.117]

22. Wang, L.; Fan, H.; Shi, G. Effect of ferritic morphology on yield strength of CGHAZ in a low carbon Mo-V-N-Ti-B steel. Metals; 2021; 11, 1863. [DOI: https://dx.doi.org/10.3390/met11111863]

23. Hu, B.; Shi, G.; Wang, Q. Elucidating the heat input on CGHAZ microstructure and its irregular effect on impact toughness for a novel V–N microalloying weathering steel. J. Mater. Res. Technol.; 2023; 25, pp. 5888-5906. [DOI: https://dx.doi.org/10.1016/j.jmrt.2023.07.086]

24. Wang, B.; Chen, N.; Cai, Y. Effect of crystallographic features on low-temperature fatigue ductile-to-brittle transition for simulated coarse-grained heat-affected zone of bainite steel weld. Int. J. Fatigue; 2023; 170, 107523. [DOI: https://dx.doi.org/10.1016/j.ijfatigue.2023.107523]

25. Zhang, X.; Gao, H.; Zhang, X. Effect of volume fraction of bainite on microstructure and mechanical properties of X80 pipeline steel with excellent deformability. Mater. Sci. Eng. A; 2012; 531, pp. 84-90. [DOI: https://dx.doi.org/10.1016/j.msea.2011.10.035]

26. Fu, Z.; Yang, B.; Shan, M. Hydrogen embrittlement behavior of SUS301L-MT stainless steel laser-arc hybrid welded joint localized zones. Corros. Sci.; 2020; 164, 108337. [DOI: https://dx.doi.org/10.1016/j.corsci.2019.108337]

27. Guo, A.; Misra, R.D.K.; Liu, J. An analysis of the microstructure of the heat-affected zone of an ultra-low carbon and niobium-bearing acicular ferrite steel using EBSD and its relationship to mechanical properties. Mater. Sci. Eng. A; 2010; 527, pp. 6440-6448. [DOI: https://dx.doi.org/10.1016/j.msea.2010.06.092]

28. Fan, L.; Zhou, D.; Wang, T. Tensile properties of an acicular ferrite and martensite/austenite constituent steel with varying cooling rates. Mater. Sci. Eng. A; 2014; 590, pp. 224-231. [DOI: https://dx.doi.org/10.1016/j.msea.2013.10.037]

29. Zhu, K.; Bouaziz, O.; Oberbillig, C. An approach to define the effective lath size controlling yield strength of bainite. Mater. Sci. Eng. A; 2010; 527, pp. 6614-6619. [DOI: https://dx.doi.org/10.1016/j.msea.2010.06.061]

30. Olasolo, M.; Uranga, P.; Rodriguez-Ibabe, J.M. Effect of austenite microstructure and cooling rate on transformation characteristics in a low carbon Nb-V microalloyed steel. Mater. Sci. Eng. A; 2011; 528, pp. 2559-2569. [DOI: https://dx.doi.org/10.1016/j.msea.2010.11.078]

31. Li, C.; Duan, R.; Fu, W. Improvement of mechanical properties for low carbon ultra-high strength steel strengthened by Cu-rich multistructured precipitation via modification to bainite. Mater. Sci. Eng. A; 2021; 817, 141337. [DOI: https://dx.doi.org/10.1016/j.msea.2021.141337]

32. Williamson, G.K.; Smallman, R.E.J.P.M., III. Dislocation densities in some annealed and cold-worked metals from measurements on the X-ray debye-scherrer spectrum. Philos. Mag. A; 1956; 1, pp. 34-46. [DOI: https://dx.doi.org/10.1080/14786435608238074]

33. Hu, J.; Du, X.; Wang, J.; Gao, R. Effect of welding heat input on microstructures and toughness in simulated CGHAZ of V-N high strength steel. Mater. Sci. Eng. A; 2013; 577, pp. 161-168. [DOI: https://dx.doi.org/10.1016/j.msea.2013.04.044]

34. Wang, X.L.; Wang, Z.Q.; Ma, X.P. Analysis of impact toughness scatter in simulated coarse-grained HAZ of E550 grade offshore engineering steel from the aspect of crystallographic structure. Mater. Charact.; 2018; 140, pp. 312-319. [DOI: https://dx.doi.org/10.1016/j.matchar.2018.03.037]

35. Zhang, S.; Yang, M.; Yuan, F. Extraordinary fracture toughness in nickel induced by heterogeneous grain structure. Mater. Sci. Eng. A; 2022; 830, 142313. [DOI: https://dx.doi.org/10.1016/j.msea.2021.142313]

36. Wang, J.; Pan, Z.; Wang, Y. Evolution of crystallographic orientation, precipitation, phase transformation and mechanical properties realized by enhancing deposition current for dual-wire arc additive manufactured Ni-rich NiTi alloy. Addit. Manuf.; 2020; 34, 101240. [DOI: https://dx.doi.org/10.1016/j.addma.2020.101240]

37. Moeinifar, S.; Kokabi, A.H.; Hosseini, H. Influence of peak temperature during simulation and real thermal cycles on microstructure and fracture properties of the reheated zones. Mater. Des.; 2010; 31, pp. 2948-2955. [DOI: https://dx.doi.org/10.1016/j.matdes.2009.12.023]

38. Baker, T.N. Microalloyed steels. Ironmak. Steelmak.; 2016; 43, pp. 264-307. [DOI: https://dx.doi.org/10.1179/1743281215Y.0000000063]

39. Tomita, Y.; Saito, N.; Tsuzuki, T. Improvement in HAZ toughness of steel by TiN-MnS addition. ISIJ Int.; 1994; 34, pp. 829-835. [DOI: https://dx.doi.org/10.2355/isijinternational.34.829]

40. Yang, L.; Peng, D.; Zhou, J. Effects of TiC on the microstructure and formation of acicular ferrite in ferritic stainless steel. Metall. Mater.; 2019; 26, pp. 1385-1395.

41. Gregg, J.M.; Bhadeshia, H. Solid-state nucleation of acicular ferrite on minerals added to molten steel. Acta Mater.; 1997; 45, pp. 739-748. [DOI: https://dx.doi.org/10.1016/S1359-6454(96)00187-5]

42. Ishikawa, F.; Takahashi, T.; Ochi, T.J.M. Intragranular ferrite nucleation in medium-carbon vanadium steels. Metall. Mater. Trans. A; 1994; 25, pp. 929-936. [DOI: https://dx.doi.org/10.1007/BF02652268]

43. Tan, X.; Lu, W.; Guo, N. Effect of tempering and partitioning (T&P) treatment on microstructure and mechanical properties of a low-carbon low-alloy quenched and dynamically partitioned (Q-DP) steel. Mater. Sci. Eng. A; 2023; 872, 144968.

44. Sun, X.; Wang, Y.; Sun, D. Microstructure and mechanical properties of Si-rich H13 steel processed via austempering and tempering. Mater. Sci. Eng. A; 2022; 834, 142616. [DOI: https://dx.doi.org/10.1016/j.msea.2022.142616]

45. Iza-Mendia, A.; Gutierrez, I. Generalization of the existing relations between microstructure and yield stress from ferrite-pearlite to high strength steels. Mater. Sci. Eng. A; 2013; 561, pp. 40-51. [DOI: https://dx.doi.org/10.1016/j.msea.2012.10.012]

46. Chen, J.; Shuai, T.; Liu, Z.Y. Microstructural characteristics with various cooling paths and the mechanism of embrittlement and toughening in low-carbon high performance bridge steel. Mater. Sci. Eng. A; 2013; 559, pp. 241-249. [DOI: https://dx.doi.org/10.1016/j.msea.2012.08.091]

47. Gutierrez, I. Effect of microstructure on the impact toughness of Nb-microalloyed steel: Generalisation of existing relations from ferrite-pearlite to high strength microstructures. Mater. Sci. Eng. A; 2013; 571, pp. 57-67. [DOI: https://dx.doi.org/10.1016/j.msea.2013.02.006]

48. Hüper, T.; Endo, S.; Ishikawa, N. Effect of volume fraction of constituent phases on the stress-strain relationship of dual phase steels. ISIJ Int.; 2007; 39, pp. 288-294. [DOI: https://dx.doi.org/10.2355/isijinternational.39.288]

49. Li, Z.; Chen, Y.; Yuan, H. Microstructure and properties of low alloy Cr-Mo steel processed with varied isothermal temperatures. Mater. Today Commun.; 2023; 35, 106158. [DOI: https://dx.doi.org/10.1016/j.mtcomm.2023.106158]

50. Sung, H.K.; Lee, D.H.; Lee, S. Correlation between microstructures and tensile properties of strain-based API X60 pipeline steels. Metall. Mater. Trans. A; 2016; 47, pp. 2726-2738. [DOI: https://dx.doi.org/10.1007/s11661-016-3453-3]

51. Bandhu, D.; Vora, J.J.; Das, S. Experimental study on application of gas metal arc welding based regulated metal deposition technique for low alloy steel. Mater. Manuf. Process.; 2022; 37, pp. 1727-1745. [DOI: https://dx.doi.org/10.1080/10426914.2022.2049298]

52. Noell, P.; Carroll, J.; Hattar, K. Do voids nucleate at grain boundaries during ductile rupture?. Acta Mater.; 2017; 137, pp. 103-114. [DOI: https://dx.doi.org/10.1016/j.actamat.2017.07.004]