Inﬂuence of Non-Metallic Inclusions on Local Deformation and Damage Behavior of Modiﬁed 16MnCrS5 Steel

: This work investigates a ferrite matrix with multiple non-metallic inclusions to evaluate their inﬂuence on the global and local deformation and damage behavior of modiﬁed 16MnCrS5 steel. For this purpose, appropriate specimens are prepared and polished. The EBSD technique is used to record local phase and orientation data, then analyze and identify the size and type of inclusions present in the material. The EBSD data are then used to run full phase crystal plasticity simulations using DAMASK-calibrated material model parameters. The qualitative and quantitative analysis of these full phase simulations provides a detailed insight into how the distribution of non-metallic inclusions within the ferrite matrix affects the local stress, strain, and damage behavior. In situ tensile tests are carried out on specially prepared miniature dog-bone-shaped notched specimens in ZEISS Gemini 450 scanning electron microscope with a Kammrath and Weiss tensile test stage. By adopting an optimized scheme, tensile tests are carried out, and local images around one large and several small MnS inclusions are taken at incremental strain values. These images are then processed using VEDDAC, a digital image correlation-based microstrain measurement tool. The damage initiation around several inclusions is recorded during the in situ tensile tests, and damage initiation, propagation, and strain localization are analyzed. The experimental results validate the simulation outcomes, providing deeper insight into the experimentally observed trends.


Introduction
The steel industry is playing an important role in providing about 45% of the raw material to the automotive sector [1]. The advancement in the steel parts used in an automobile is inevitable to competitively improve the strength to weight ratio, fuel economy, and crash safety of cars. Researchers have studied various multi-phase steels to increase the ultimate tensile strength (UTS) and percentage elongation of steels using different thermo-mechanical methods and phase combinations [2,3].
The ever-increasing demand for lightweight materials for the mobility sector has promoted the utilization of the multi-phase phenomenon in steel [4]. For instance, quenched and partitioned (Q&P) steels have shown a UTS up to 1.4 GPa with a uniform elongation of up to 20% [5]. This remarkable combination is desirable during steel formation into sheets and during individual component making with application-based mechanical strength capabilities. Small-sized hard phase inclusions added into the soft ferrite matrix during steel making act as reinforcement participants [6][7][8]. The common inclusions in the ferrite matrix-i.e., alumina (Al 2 O 3 ), cementite (Fe 3 C), manganese sulfide (MnS), and pyrite (Fe 2 S)-are very small in size,

Selection of Area for Full Phase Simulations
For the crystal structure analysis, the specimens were placed in the Gemini SEM450 from Carl Zeiss AG. The EBSD analysis was performed with an accelerating voltage of 20 kV, a working distance of 11 mm, and a specimen holder pre-tilted by 70°. The magnification was 10,000×. For the EBSD analysis, the Symmetry S2 detector from Oxford Instruments PLC was used. The analysis area was 95.4 µm × 71.8 µm, with a step size of 0.2 µm. The advanced Symmetry S2 EBSD detector utilizes a unique combination of speed (4500 pixels/second), sensitivity (CMOS technology sensor), and diffraction pattern details (1244 × 1024 pixel resolution) to efficiently produce high-quality EBSD data for the selected area, which can then be post-processed to obtain information in the desired format. Modified 16MnCrS5 steel contains inclusions that are quantitatively analyzed and recoded using EDS Ultim MaxX from Oxford Instruments PLC. The overall chemical composition of the current material is shown in Table 1 in comparison with the standard chemical composition of the alloy. It is modified 16MnCrS5 with 50% lower carbon content and 20% more sulfur. Therefore, the strength of this modified material is slightly lower; due to higher sulfur content and less carbon, it contains more MnS inclusions, with some being predominantly large (as discussed later in Sections 2.3 and 2.4). The in situ tensile tests were carried out with a ZEISS Gemini SEM 450 scanning electron microscope using a Kammrath and Weiss tensile test stage, as shown in Figure 1a,b. This module allows in situ data collection during tensile or compressive loading using SE, EBSD, and EDS in the Gemini SEM450 of Carl Zeiss AG. For the experiments, different points on the tensile stress-strain curve were targeted at a forming speed of 8 µm/s, and the specimen was held at various values of stress and strain. At the same time, analyses of the microstructure and chemistry were performed. A series of images was taken during tensile testing. All images were taken with a resolution of 2048 × 1536 pixels. In addition, all photos in the area of interest were taken with a magnification of 11k and WD (working distance) of 16.6 mm with the "InLens" module.

Selection of Area for Full Phase Simulations
For the crystal structure analysis, the specimens were placed in the Gemini SEM450 from Carl Zeiss AG. The EBSD analysis was performed with an accelerating voltage of 20 kV, a working distance of 11 mm, and a specimen holder pre-tilted by 70 • . The magnification was 10,000×. For the EBSD analysis, the Symmetry S2 detector from Oxford Instruments PLC was used. The analysis area was 95.4 µm × 71.8 µm, with a step size of 0.2 µm. The advanced Symmetry S2 EBSD detector utilizes a unique combination of speed (4500 pixels/second), sensitivity (CMOS technology sensor), and diffraction pattern details (1244 × 1024 pixel resolution) to efficiently produce high-quality EBSD data for the selected area, which can then be post-processed to obtain information in the desired format. Modified 16MnCrS5 steel contains inclusions that are quantitatively analyzed and recoded using EDS Ultim MaxX from Oxford Instruments PLC.
The multi-phase steel material Representative Volume Element (RVE) was obtained after EBSD analysis. The data were initially raw and required cleaning for unindexed points due to crystallographic noise at the grain boundaries. This noise was reduced using an intelligent algorithm implemented on the MTEX toolbox with MATLAB [40]. The procedural details of the EBSD data cleaning and the developed algorithm are published elsewhere [41]. Therefore, readers are encouraged to refer to that work for further details.
The ferrite matrix distribution is shown in Figure 2a, and the distributions of alumina (Al 2 O 3 ), cementite (Fe 3 C), and manganese sulfide (MnS) concerning the grain boundaries are displayed in Figure 2b. Most of the inclusions detected in this area of the specimens are present on the grain boundaries. This kind of microstructure with appropriately large second phase particles distributed homogeneously on the grain boundaries of adequately large ferrite grains appears to be due to the adopted manufacturing strategy.

Statistical Analysis of the EBSD Data
A quantitative comparison of the grain shape (aspect ratio) and particle size of the ferrite matrix and inclusions is presented in Figures 3 and 4, respectively. In Figure 3, the aspect ratio of 1 corresponds to almost perfectly round shapes, whereas the higher aspect ratio corresponds to the adequately elongated grain or inclusion. Considering this criterion, the aspect ratio of ferrite grains in Figure 3a predominantly lies between 1-2, except for some grains with a high degree of elongation of >4, which is the normal ferrite microstructure reported in previous literature [42,43]. On the other hand, the inclusions in Figure 3b-d are largely round (aspect ratio = 1), with some particles slightly elongated. There can be two reasons why the manufacturing process yields such inclusion shapes, as reported previously [44], or it can be due to the limitations of the measurement and magnification chosen in this work. During EBSD analysis, only these slightly large inclusions with round profiles were detected due to the chosen magnification and the set step size. In contrast, very thin and long inclusions were missed.

Statistical Analysis of the EBSD Data
A quantitative comparison of the grain shape (aspect ratio) and particle size of the ferrite matrix and inclusions is presented in Figures 3 and 4, respectively. In Figure 3, the aspect ratio of 1 corresponds to almost perfectly round shapes, whereas the higher aspect ratio corresponds to the adequately elongated grain or inclusion. Considering this criterion, the aspect ratio of ferrite grains in Figure 3a predominantly lies between 1-2, except for some grains with a high degree of elongation of >4, which is the normal ferrite microstructure reported in previous literature [42,43]. On the other hand, the inclusions in Figure 3b-d are largely round (aspect ratio = 1), with some particles slightly elongated. There can be two reasons why the manufacturing process yields such inclusion shapes, as reported previously [44], or it can be due to the limitations of the measurement and magnification chosen in this work. During EBSD analysis, only these slightly large inclusions with round profiles were detected due to the chosen magnification and the set step size. In contrast, very thin and long inclusions were missed.
The grain size in Figure 4 is measured in terms of the mean sphere diameter. A significant number of ferrite grains in the selected area for EBSD analysis are below 4 µm, as observed in Figure 4a, and some grains are larger than 6 µm. As shown in Figure 4bd, the inclusions are comparatively very small. Most inclusions are around 0.1 µm in size, whereas very few are larger than 0.2 µm. Even smaller inclusions may be present within the microstructure but were not detected due to the selected magnification and EBSD measurement step size. This lack of smaller inclusions negligibly alters the material chemistry or local phase distribution and hence is assumed to have almost no effect on the simulation results.    . Distribution of aspect ratios of (a) ferrite grains, which are relatively larger and have a wider spread; (b) cementite inclusions, which are close to 1, with few inclusions at 1.3 and 1.7; (c) MnS inclusions are also usually round, with a few with a spread of 2.8; (d) Al2O3 inclusions are also usually round with an aspect ratio of 1, with some inclusions reaching up to 2.3.  The grain size in Figure 4 is measured in terms of the mean sphere diameter. A significant number of ferrite grains in the selected area for EBSD analysis are below 4 µm, as observed in Figure 4a, and some grains are larger than 6 µm. As shown in Figure 4b-d, the inclusions are comparatively very small. Most inclusions are around 0.1 µm in size, whereas very few are larger than 0.2 µm. Even smaller inclusions may be present within the microstructure but were not detected due to the selected magnification and EBSD measurement step size. This lack of smaller inclusions negligibly alters the material chemistry or local phase distribution and hence is assumed to have almost no effect on the simulation results.
Second phase inclusions have been talked about a great deal in this section. The elemental analysis of an MnS inclusion using EDS analysis is presented in Figure 5 for reference. All the other inclusions-i.e., Al 2 O 3 and Fe 3 C-have also been identified using the same process.

Selection of Area, Tools Used, and Methodology Adopted for In Situ Strain Measurement
During the in situ tensile tests, a large, elongated sulfide inclusion (MnS) was analyzed. The inclusion, along with the chemical analysis confirming it to be MnS, is shown in Table 2. The length of the inclusion together with the cavity is about 42.5 µm. The elongated inclusion aligned with the tensile loading axis (horizontal to the shown micrograph) was expected to be prone to severe fracture and was therefore of great interest. as 0.23 µm. (d) The Al2O3 inclusions have a large spread with many having sizes close to 0.1 µm and then having a downslope spread up to 0.25 µm ER.
Second phase inclusions have been talked about a great deal in this section. The elemental analysis of an MnS inclusion using EDS analysis is presented in Figure 5 for reference. All the other inclusions-i.e., Al2O3 and Fe3C-have also been identified using the same process.

Selection of Area, Tools Used, and Methodology Adopted for In Situ Strain Measurement
During the in situ tensile tests, a large, elongated sulfide inclusion (MnS) was analyzed. The inclusion, along with the chemical analysis confirming it to be MnS, is shown in Table 2. The length of the inclusion together with the cavity is about 42.5 µm. The elongated inclusion aligned with the tensile loading axis (horizontal to the shown micrograph) was expected to be prone to severe fracture and was therefore of great interest.
The in situ tests come with their own challenge of image acquisition [35][36][37]. The images taken during the test should have a specific contrast and clear feature identification possibility [45]. To ensure that the images taken during the current tests are of acceptable quality and grayscale distribution, an analysis was performed that has not been included in the main part of the body but is presented in Appendix A. The effect of different settings and attributes is discussed and can be interesting for readers interested in performing such tests themselves. Appendix A also contains detailed information about the local and global strain measurement during experimentation.   The digital image correlation software VEDDAC from Chemnitzer Werkstoffmechanik GmbH was used to process the captured image series. When processing the images, images with strong contrast or sharpness fluctuations were removed from the calculations, and the final set contained images. Two calculation areas were selected during DIC. The first calculation area is shown in Figure 6a, which consists of the matrix around the inclusions and the cavities. The second calculation is the inclusion area and cavities around it, as shown in Figure 6b. The optimal grid sizes for a particular series of images were 9 × 9 pixels (coarse) for the matrix and 5 × 5 pixels for the inclusion (fine). More information about the local strain measurement methodology is provided in Appendix B. The in situ tests come with their own challenge of image acquisition [35][36][37]. The images taken during the test should have a specific contrast and clear feature identification possibility [45]. To ensure that the images taken during the current tests are of acceptable quality and grayscale distribution, an analysis was performed that has not been included in the main part of the body but is presented in Appendix A. The effect of different settings and attributes is discussed and can be interesting for readers interested in performing such tests themselves. Appendix A also contains detailed information about the local and global strain measurement during experimentation. The digital image correlation software VEDDAC from Chemnitzer Werkstoffmechanik GmbH was used to process the captured image series. When processing the images, images with strong contrast or sharpness fluctuations were removed from the calculations, and the final set contained images. Two calculation areas were selected during DIC. The first calculation area is shown in Figure 6a, which consists of the matrix around the inclusions and the cavities. The second calculation is the inclusion area and cavities around it, as shown in Figure 6b. The optimal grid sizes for a particular series of images were 9 × 9 pixels (coarse) for the matrix and 5 × 5 pixels for the inclusion (fine). More information about the local strain measurement methodology is provided in Appendix B. The digital image correlation software VEDDAC from Chemnitzer Werkstoffmechanik GmbH was used to process the captured image series. When processing the images, images with strong contrast or sharpness fluctuations were removed from the calculations, and the final set contained images. Two calculation areas were selected during DIC. The first calculation area is shown in Figure 6a, which consists of the matrix around the inclusions and the cavities. The second calculation is the inclusion area and cavities around it, as shown in Figure 6b. The optimal grid sizes for a particular series of images were 9 × 9 pixels (coarse) for the matrix and 5 × 5 pixels for the inclusion (fine). More information about the local strain measurement methodology is provided in Appendix B.
(a) (b) Figure 6. Calculation network for digital image correlation (a) around non-metallic inclusions (b) in non-metallic inclusions and cavity area. In this work, a grid of 9 × 9 pixels (coarse) for the matrix and 5 × 5 pixels for the inclusion (fine) was chosen to get accurate and fully resolved results. In this work, a grid of 9 × 9 pixels (coarse) for the matrix and 5 × 5 pixels for the inclusion (fine) was chosen to get accurate and fully resolved results.

Numerical Simulation Model Setup
The EBSD data and the information of each phase's chemical and crystallographic structure were used to model and run a full phase crystal plasticity-based numerical model with calibrated material models. These simulations were used to provide detailed information about the local deformation and damage behavior of this material and how different inclusions, morphologies, and distributions affect the local deformation and damage behaviors. This insight is useful to understand the data from limited in situ tests.
First developed by Hutchinson [46] and later extended by Kalidindi [47], the phenomenological model implemented in the DAMASK framework [33] calculates the plastic deformation by slip planes based on the initial and saturated slip resistance-S 0 , S s , respectively (refer to Table 3). This resistance value for slip systems = 1,2,3 . . . N slip increases from the initial to saturation value depending on the crystal structure specifications and the critical shear value on the slip plane during deformation.
This model uses a mathematical description to correlate the initial deformation gradient (F) with the resulting first Piola-Kirchoff stress tensor (P), as shown in Equations (1) and (2). The elastic part is simplified in the form of the generalized Hook's law, and for the plastic part, it is given with the help of the plastic velocity gradient given in Equation (3).
The value of the velocity plasticity gradient L p is dependent on the shear strain γ (Equation (4)), which is a function of the resolved shear stress τ α and slip-resistance on the α slip plane S α .
and α = 1, 2, 3 · · · , N slip (For BCC Ferrite N slip = 24) A non-conserved damage field ϕ is governed by the continuous release of the stored mechanical energy density from undamaged to fully damaged conditions in a region; i.e., ϕ = 1 to ϕ = 0. In the ductile damage criterion implemented in DAMASK, plastic energy dissipation at a material point is the driving force here for damage flux f ϕ. Therefore, plastic energy dissipation in the form of the following equation is used: Consequently, the minimization of the total free energy density is responsible for force to drive damage, and it is given as follows: where G = surface tension of the newly generated damage surface; l c = length scale of the diffused surface; L p = plastic velocity gradient. The plastic governing law changes to the following when damage evolution is coupled with dissipated plastic energy: Readers are encouraged to read the pioneering work by the developers of DAMASK for further details and understanding of the model [33,34,48,49]. For instance, a strategy was published recently by Qayyum et al. [50] by choosing an optimal RVE that behaves isotopically yet is small enough to produce fast results for the calibration of the material model parameters of single-phase materials. The strategy comprises a 10 × 10 × 10 size RVE containing 1000 grains, where each point represents a different grain. The RVE adopted from that published methodology is shown in Figure 7a. This RVE was used to calibrate the fitting parameters in the material model by comparing simulation outcomes with the experimental observations, and the results are shown in Figure 7b.
It is observed that the simulation results match well with the experimental observations globally. The dips in the experimental data are positions where the test was stopped and pictorial data were recorded. The adopted and calibrated material model parameters used in this study are shown in Table 3. This set of parameters was used to run full phase simulations and then analyze the local material deformation and damage behavior.
The DAMASK framework is designed to get all the inputs in the form of specifically structured text files with a complete microstructural description of RVE with geometrical details, material behavior attributes, and boundary values. Elastic and plastic phase parameters for ferrite and only elastic parameters of all the inclusions are defined to observe the micromechanical deformation of ferrite caused by stiffness degradation. Ductile damage criteria have been incorporated in the ferrite matrix as well. The RVE generated by using the EBSD data shown in Figure 2 was subjected to quasi-static tensile load with a strain rate of 10 −3 s −1 in the x-direction.
RVE containing 1000 grains, where each point represents a different grain. The RVE adopted from that published methodology is shown in Figure 7a. This RVE was used to calibrate the fitting parameters in the material model by comparing simulation outcomes with the experimental observations, and the results are shown in Figure 7b.
It is observed that the simulation results match well with the experimental observations globally. The dips in the experimental data are positions where the test was stopped and pictorial data were recorded. The adopted and calibrated material model parameters used in this study are shown in Table 3. This set of parameters was used to run full phase simulations and then analyze the local material deformation and damage behavior. The DAMASK framework is designed to get all the inputs in the form of specifically structured text files with a complete microstructural description of RVE with geometrical details, material behavior attributes, and boundary values. Elastic and plastic phase parameters for ferrite and only elastic parameters of all the inclusions are defined to observe the micromechanical deformation of ferrite caused by stiffness degradation. Ductile damage criteria have been incorporated in the ferrite matrix as well. The RVE generated by using the EBSD data shown in Figure 2 was subjected to quasi-static tensile load with a strain rate of 10 −3 s −1 in the x-direction.  Table 3. Elastic parameters of second phase inclusions obtained from the literature [10,20,23,26,28,41,51,52] and calibrated parameters of ferrite adopted from the literature [23,26,28,41,52] and adjusted by comparing with the flow curve from in situ tests. The ductile damage parameters were calibrated by comparing them with the in situ test results.

Global Results of Numerical Simulation
The global results of the simulation in comparison with experimental observations are shown in Figure 8. The global results are calculated by taking an average of the results at each solution point for each increment. The data are combined in a meaningful way to produce observable trends. The stress and damage evolution with respect to the true global strain in the selected RVE is shown here in Figure 8.

Local Results of Numerical Simulation
The local distributions of strains, stress, triaxiality, and damage at 3%, 9%, and 15.5% of global strain are shown in Figure 9. Due to the heterogeneous distribution of nonmetallic inclusions and different Schmidt factors of each ferrite grain, the local distributions of each attribute are largely heterogeneous. Therefore, the grain boundaries have been overlaid in contrasting colors to clarify the distributions of attributes within individual grains.
It is observed that at 3% global strain, the local strain distribution is relatively homogeneous within the selected RVE. As the global strain increases and reaches a maximum of 15.5%, an extreme heterogeneous distribution of strain along the grain boundaries in thick channels aligned at 45 degrees to the loading axis is observed with intermediate areas of very low local strain. On the other hand, the stress distribution is high from the beginning in the ferrite grains with a high Schmidt factor and increases significantly within those grains. The value remains similar as the global strain reaches the maximum value. Of course, the local stress and strain distribution is largely affected by the presence and distribution of non-metallic inclusions. The behavior is discussed in more detail in the next section.
Stress triaxiality is the relative degree of hydrostatic stress in each stress state, which is usually used to estimate the type of fracture that might take place in the material. In the current RVE, it is observed that the triaxiality measure in most of the grains is close to 0.1, which means that the stress (no matter how high) is predominantly hydrostatic. Therefore, the deviatoric component of stress in these grains is low and will not contribute significantly towards material flow and eventual damage. In addition, it is interesting to note that the triaxial stress remains consistent with the increase in strain and only slightly shifts to neighboring grains after the damage initiation and propagation.
Generally, it is observed that the damage initiates at the interface of the non-metallic inclusions, especially in the narrow zones locked between the cluster of inclusions. After A good correlation between global experimental and numerical simulation results is observed with less than 3% difference of predictions. The inset figures in Figure 8 from S1-S5 show damage evolution around an MnS inclusion. Inset figures E1-E4 show damage evolution around a similar-sized MnS inclusion captured during in situ testing. Interestingly, not only do the global results match well, but the local damage evolution is also comparable. These local deformation and damage trends are discussed later in this article.

Local Results of Numerical Simulation
The local distributions of strains, stress, triaxiality, and damage at 3%, 9%, and 15.5% of global strain are shown in Figure 9. Due to the heterogeneous distribution of non-metallic inclusions and different Schmidt factors of each ferrite grain, the local distributions of each attribute are largely heterogeneous. Therefore, the grain boundaries have been overlaid in contrasting colors to clarify the distributions of attributes within individual grains.
It is observed that at 3% global strain, the local strain distribution is relatively homogeneous within the selected RVE. As the global strain increases and reaches a maximum of 15.5%, an extreme heterogeneous distribution of strain along the grain boundaries in thick channels aligned at 45 degrees to the loading axis is observed with intermediate areas of very low local strain. On the other hand, the stress distribution is high from the beginning in the ferrite grains with a high Schmidt factor and increases significantly within those grains. The value remains similar as the global strain reaches the maximum value.
Of course, the local stress and strain distribution is largely affected by the presence and distribution of non-metallic inclusions. The behavior is discussed in more detail in the next section. damage accelerates. It is important to mention here that the continuum mechanics simulations are intrinsically unstable due to the loss of stiffness at several solution points. Furthermore, it takes immense computing power to calculate the corresponding equilibrium for each increment after damage initiates in the RVE. Therefore, the simulations were only computed up to 15.5% of the global strain, and the results are shown in Figure 9. Although the local results generally for the ferrite matrix are shown in Figure 9, one of the main objectives of this article is to investigate the effect of non-metallic inclusions on the material's local deformation and damage behavior. In Figure 10, the focus has specifically been shifted towards three different zones in the RVE of comparable size with different inclusion compositions and distributions. Zone-I comprises a single grain with a large MnS inclusion in the middle. Zone-II comprises several small ferrite grains with a clustered distribution of non-metallic inclusions primarily on the grain boundaries. Zone-III comprises relatively larger ferrite grains with a diffused distribution of non-metallic inclusions within the grains and on the grain boundaries. The distribution of local attributes in the three zones is selectively identified in the corresponding subplots Figure  10b-e. There is a relatively low stress concentration in zone-I and less strain due to almost nonexistent stress triaxiality, and therefore, no damage takes place in this zone. This can be due to the high Schmidt factor of this specific grain and the central position of a large hard MnS inclusion which restricts the slip planes from moving in this grain freely. Stress triaxiality is the relative degree of hydrostatic stress in each stress state, which is usually used to estimate the type of fracture that might take place in the material. In the current RVE, it is observed that the triaxiality measure in most of the grains is close to 0.1, which means that the stress (no matter how high) is predominantly hydrostatic. Therefore, the deviatoric component of stress in these grains is low and will not contribute significantly towards material flow and eventual damage. In addition, it is interesting to note that the triaxial stress remains consistent with the increase in strain and only slightly shifts to neighboring grains after the damage initiation and propagation.
Generally, it is observed that the damage initiates at the interface of the non-metallic inclusions, especially in the narrow zones locked between the cluster of inclusions. After initiations, the damage propagates at an oblique angle to the applied external load. As a result, the local voids coalesce and form larger cracks that grow faster, and material damage accelerates. It is important to mention here that the continuum mechanics simulations are intrinsically unstable due to the loss of stiffness at several solution points. Furthermore, it takes immense computing power to calculate the corresponding equilibrium for each increment after damage initiates in the RVE. Therefore, the simulations were only computed up to 15.5% of the global strain, and the results are shown in Figure 9.
Although the local results generally for the ferrite matrix are shown in Figure 9, one of the main objectives of this article is to investigate the effect of non-metallic inclusions on the material's local deformation and damage behavior. In Figure 10, the focus has specifically been shifted towards three different zones in the RVE of comparable size with different inclusion compositions and distributions. Zone-I comprises a single grain with a large MnS inclusion in the middle. Zone-II comprises several small ferrite grains with a clustered distribution of non-metallic inclusions primarily on the grain boundaries. Zone-III comprises relatively larger ferrite grains with a diffused distribution of non-metallic inclusions within the grains and on the grain boundaries. The distribution of local attributes in the three zones is selectively identified in the corresponding subplots Figure 10b-e. There is a relatively low stress concentration in zone-I and less strain due to almost nonexistent stress triaxiality, and therefore, no damage takes place in this zone. This can be due to the high Schmidt factor of this specific grain and the central position of a large hard MnS inclusion which restricts the slip planes from moving in this grain freely. Small non-metallic inclusions distributed on the grain boundaries of small ferrite grains act as high-stress concentration sites with the large plastic flow in zone-II. Consequently, the stress triaxiality in this zone is relatively high, indicating a large plastic flow and evolution of ductile damage. This damage evolution due to the coalescence of small voids that form on the inclusion/matrix interface is visible in Figure 10e. It is observed that local damage incidents join together to form microcracks that propagate at maximum shear plane oriented 45 degrees to the loading axis. This can be due to two reasons: first, due to the 2D nature of the RVE, as has been pointed out in the previous work [52]; and second, due to the high amount of slip system activation in grains and areas with maximum critical resolved shear.
Zone-III comprises relatively larger ferrite grains with a relatively broader distribution of non-metallic inclusions within the grain and grain boundaries. In zone-III, although the triaxiality is comparable with zone-II, due to the lack of clustering on nonmetallic inclusions, the stress distribution and strain distribution are lower. Hence, no damage occurs in the zone.  Small non-metallic inclusions distributed on the grain boundaries of small ferrite grains act as high-stress concentration sites with the large plastic flow in zone-II. Consequently, the stress triaxiality in this zone is relatively high, indicating a large plastic flow and evolution of ductile damage. This damage evolution due to the coalescence of small voids that form on the inclusion/matrix interface is visible in Figure 10e. It is observed that local damage incidents join together to form microcracks that propagate at maximum shear plane oriented 45 degrees to the loading axis. This can be due to two reasons: first, due to the 2D nature of the RVE, as has been pointed out in the previous work [52]; and second, due to the high amount of slip system activation in grains and areas with maximum critical resolved shear.
Zone-III comprises relatively larger ferrite grains with a relatively broader distribution of non-metallic inclusions within the grain and grain boundaries. In zone-III, although the triaxiality is comparable with zone-II, due to the lack of clustering on non-metallic inclusions, the stress distribution and strain distribution are lower. Hence, no damage occurs in the zone.
In Figure 10, a qualitative comparison of local attributes in three different zones shows how the distribution of non-metallic inclusions plays a role in affecting the local stress, strain, and damage. To quantitively compare three zones, a statistical normalized probability distribution was adopted to systematically compare the distribution of stress and strain within zone-I, zone-II, and zone-III, as shown in Figure 11. In Figure 10, a qualitative comparison of local attributes in three different zones shows how the distribution of non-metallic inclusions plays a role in affecting the local stress, strain, and damage. To quantitively compare three zones, a statistical normalized probability distribution was adopted to systematically compare the distribution of stress and strain within zone-I, zone-II, and zone-III, as shown in Figure 11. Figure 11. The normalized frequency of (a) strain and (b) stress for each zone identified in Figure  10.
To quantitatively compare the stress and strain distribution in zones I to III, a Figure 11. The normalized frequency of (a) strain and (b) stress for each zone identified in Figure 10.
To quantitatively compare the stress and strain distribution in zones I to III, a statistical distribution tool is used, and the results are presented in Figure 11. It is observed that the strain distribution in all three zones, as shown in Figure 11a, is very similar. It peaks around 0.15-0.18 and then linearly drops to a strain of 0.45. The strain in zone-I is slightly higher and that in zone-II is slightly lower due to the local orientation distribution. The stress distributions in these zones reach up to 1100 MPa.
The statistical stress distribution in the three zones is presented in Figure 11b and proves the already stated hypothesis. It is observed that the stress in zone-I and zone-III peaks around 750 MPa, with a slight bias towards the lower stress side. On the other hand, the stress distribution in zone-II, due to the close clustering of several non-metallic inclusions on the grain boundaries of small ferrite grains, peaks at 850 MPa (11% higher than the other zones), and its distribution is biased towards the higher-stress side. The distribution reaches 1250 MPa (12% higher than the other zones). This high local attribute distribution in zone-II makes it more prone to local damage initiation and propagation in high deformation regimes.

Local Results of the In Situ Tensile Test
The local results of the in situ tensile tests at 3%, 5%, and 8% of global strain are shown in Figure 12 The methodology and procedure of local strain measurement is provided in Appendix B. The complete set of results, due to its extensiveness, is provided in Appendix C of this article. Readers can refer to Appendix B and C for further details about collection and frame by frame evolution of local strain measurements.

Damage Evolution around Non-Metallic Inclusions
The MnS inclusion tracked to analyze local deformation and damage behavior in Figure 12 is an exception. Usually, the inclusion size in the material user consideration is between 1-4 µm; please refer to Figure 4 for reference. Unfortunately, such small In Figure 12, at a subsequent global strain of 3% and 400 MPa global stress, it can be observed that in the zones with increased elongation, the strain rate is much higher than in the zones where the strain rate has not increased much from the beginning of the test. The effect of the increase in elongation at the ferrite grain boundaries near the non-metallic inclusion is preserved. At this stress, two zones with a significant increase in strain can be seen in the middle of the inclusion. When comparing the entire image series, at a load of 450 MPa and global strain of 5%, the sulfide inclusion is cracked in the area of thinning. The inclusion in the left part of the crack is brighter due to its change in position, which significantly affects the strain values in this inclusion area, as the digital image correction software requires areas with contrast for accurate calculation. The breaking point itself is hidden at this point due to the topography.
The maximum values in the range of increased load are above 8.9%. It is observed that the two sulfide inclusions move further away from each other. The strain in the region increases, reaching a value of over 13-64%. Around the formed cavity of a small inclusion located on the upper right side of the elongated inclusion, the elongation near the pearlitic zone is greater than on the other side of the same grain. On the side of the ferrite grain around this cavity is a zone of increased deformation with a value of about 7-11%. There is a significant difference in elongation within and around the small cavities. In the middle cavity, the deformation is minimal and is between 2% and 3.5%, while in the left cavity, the deformation around and within the cavity is between 6.4% and 9.8%.
Due to further load increase and subsequent high strain values, the DIC program loses previously captured pixels for the calculation. By enlarging the gaps created, some pixels can significantly change their gray tone, causing the program to find a similar range of pixels elsewhere and greatly distorting the calculation values. At a maximum global strain of 8% (450 MPa), a large loss of pixels is observed (in areas beyond the specified deformation limit of dark red) because the topography changes too much, and the color range continues to change in many areas. In addition, the distance between the halves of the sulfide inclusion become significantly large. In the lower deformation zone in the matrix around the sulfide grain inclusion, the values are 43.6%. In contrast, the zone-I around the inclusion next to the perlite grain remains below the total deformation.

Damage Evolution around Non-Metallic Inclusions
The MnS inclusion tracked to analyze local deformation and damage behavior in Figure 12 is an exception. Usually, the inclusion size in the material user consideration is between 1-4 µm; please refer to Figure 4 for reference. Unfortunately, such small inclusions are hard to track and record local strain around due to high magnification, resulting in fast drifting during testing in the SEM chamber [53]. Despite this paradox, in the current work, during in situ testing, a few MnS inclusions of small size were tracked and recorded at different strains. The results are shown in Figures 13 and 14.
In Figures 13 and 14, several inclusions classified as MnS inclusions based on the EDS analysis were tracked at different external strains that were applied. The results shown here were at the beginning of the test when the inclusions were identified, at external load when the damage was initiated, and an external load when the damage started to evolve.
In Figure 13, it is observed that the damage initiates perpendicular to the line of action of the load at about 525 MPa global stress and propagates in that same direction. It is important to keep in mind that the inclusion is embedded in the matrix, and here only the surface is being observed. The surface matrix seems to be extremely deformed due to the internal propagation of the damage around the inclusion. It is also interesting to note that the damage initiates at the inclusion/matrix interface. In contrast, the inclusion remains free of any cracking, unlike the case of large inclusion, which was shown and discussed earlier in Figures 5 and 12. This is because the damage of the inclusion is a function of the size and distribution in the matrix, where small inclusions with an aspect ratio of less than 2.5 are generally not highly prone to damage. size and distribution in the matrix, where small inclusions with an aspect ratio of less than 2.5 are generally not highly prone to damage. In Figure 14, two other inclusions are tracked, again identified as MnS based on the EDS analysis shown in the same figure. Both inclusions are completely embedded in small ferrite grains and resemble the case of zone-I in Figure 10.
In Figure 14, it is observed that the damage starts to initiate at 550 MPa of external load, which is slightly higher than the damage initiation load for other inclusions but is in a similar range. At 575 MPa of external load, damage on the inclusion matrix interface grows so extensive that it is easily identifiable and is marked with the help of white arrows. In this case, it is observed that voids initiate at an oblique angle to the applied external load and not perpendicular to the load as observed earlier. This position of the damage initiation depends on several factors: i.e., orientation, size, shape, and sharpness of the inclusion edges and the place of the inclusion (inside a large grain or on the grain boundary). It would also largely depend on the orientation of the surrounding matrix In Figure 14, two other inclusions are tracked, again identified as MnS based on the EDS analysis shown in the same figure. Both inclusions are completely embedded in small ferrite grains and resemble the case of zone-I in Figure 10.
In Figure 14, it is observed that the damage starts to initiate at 550 MPa of external load, which is slightly higher than the damage initiation load for other inclusions but is in a similar range. At 575 MPa of external load, damage on the inclusion matrix interface grows so extensive that it is easily identifiable and is marked with the help of white arrows. In this case, it is observed that voids initiate at an oblique angle to the applied external load and not perpendicular to the load as observed earlier. This position of the damage initiation depends on several factors: i.e., orientation, size, shape, and sharpness of the inclusion edges and the place of the inclusion (inside a large grain or on the grain boundary). It would also largely depend on the orientation of the surrounding matrix grains. As mentioned earlier, all these inclusions have 3D geometry, and the third dimension would also affect the results observed on the surface of the specimen. grains. As mentioned earlier, all these inclusions have 3D geometry, and the third dimension would also affect the results observed on the surface of the specimen.

Discussion
The global and local deformation and damage behavior of modified 16MnCrS5 steel (composition in Table 1) is analyzed. Before the in situ testing, the specimens were prepared using a special grinding and polishing technique (presented in detail in Appendix A), and detailed EBSD data were collected to understand the morphology of the matrix and the second phase inclusions. The outcome of the statistical analysis of the collected microstructural data is presented in Figures 3 and 4. It is evident that, generally, the non-metallic inclusions are quite small, and the MnS inclusions are slightly larger and more visible at the set magnification. Although the MnS inclusions were the focus of this work, other inclusions are also expected to behave similarly.
The parameters of the crystal plasticity model were calibrated and then incorporated in the full phase model employing the recorded and cleaned EBSD data. It is observed that the global results of the simulation match well with the experimental observations. Therefore, researchers can successfully use the identified parameters to model this material using DAMASK in the future. However, one can argue that this match of global results is because a similar data set was used to calibrate these parameters. To answer that

Discussion
The global and local deformation and damage behavior of modified 16MnCrS5 steel (composition in Table 1) is analyzed. Before the in situ testing, the specimens were prepared using a special grinding and polishing technique (presented in detail in Appendix A), and detailed EBSD data were collected to understand the morphology of the matrix and the second phase inclusions. The outcome of the statistical analysis of the collected microstructural data is presented in Figures 3 and 4. It is evident that, generally, the non-metallic inclusions are quite small, and the MnS inclusions are slightly larger and more visible at the set magnification. Although the MnS inclusions were the focus of this work, other inclusions are also expected to behave similarly.
The parameters of the crystal plasticity model were calibrated and then incorporated in the full phase model employing the recorded and cleaned EBSD data. It is observed that the global results of the simulation match well with the experimental observations. Therefore, researchers can successfully use the identified parameters to model this material using DAMASK in the future. However, one can argue that this match of global results is because a similar data set was used to calibrate these parameters. To answer that critique, the damage pattern around a comparable inclusion in experiments and simulations was made at different strains, and a remarkable similarity in the trends was observed.
Apart from the global results, the local simulation results provide great insight into how the inclusion size, position within the matrix, and their distribution affect the local stress, strain, damage, and triaxiality. The trends of local strain distribution match well with the previously published data [28,54]. A qualitative and quantitative comparison of the three different zones was presented in Section 4.2. All three zones are intrinsically different based on the composition, inclusion size, and distribution. In addition, there is a clear difference in the damage behavior in these zones, and local stress distribution is observed. These local attributes account for the different local material behavior observed in the experimental part of this work.
Since the previous works have shown that the second phase inclusions are major players in defining the formability of the material [15,41,55], a special focus was given to them in this work. Previous work [51][52][53][54][55][56] shows that the morphology and distribution of the second phase inclusions dictate the formability limit and damage degradation in a material. In this work, in situ tests were carried out focusing on the strain evolution and damage evolution around and within these second phase inclusions. The second phase inclusion size distribution is observed to be close to 1 µm with a few larger inclusions in the specimen. The local strain around and within one of these exceptionally large MnS inclusions is studied in detail. MnS inclusions, due to their larger size, are easily identifiable at the magnification range selected in this study, and therefore they appear predominantly in this work. However, it is assumed that other inclusions behave similarly. This is also verified by the simulation results presented in Section 4.2, where all the inclusions behave almost indifferently under applied external load.
The results of the local strain measurement from the in situ tests are presented in Section 4.3. It is observed that the strain in the matrix is higher at some points and lower at some points, depending on the orientation of the corresponding ferrite grain. The inclusion is observed to have a low strain distribution until a brittle fracture appears, which damages the whole microstructure, and it becomes hard to track the local strain in the material further. Small MnS inclusions were also tracked in this work at increasing external loads. For these smaller inclusions, it is observed that matrix/inclusion interface decohesion starts to take place and grow at high strain regimes of >550 MPa. These findings match well with the previously published work of other researchers [4,6,19,44,56].
Although similar research with in situ tests and crystal plasticity simulations has been previously carried out by some researchers before [21,27,28,57,58], it was either pure simulation or pure experimental work. No correlation using both methods has been presented before concerning how the inclusion distribution, size, and morphology affect the material deformation and degradation. However, in this work, the findings of the developed crystal plasticity-based numerical simulation model have been validated. Researchers and industry can now use this tool to analyze and optimize the non-metallic inclusion size, distribution, and morphology to attain the desired material formability.

Conclusions
The global and local deformation and damage behavior of modified 16MnCrS5 steel using in situ experimental and crystal plasticity-based numerical simulation model is analyzed. Recorded EBSD data were used to evaluate the local microstructure of the material. The data were then adopted to run full phase crystal plasticity simulations in DAMASK with calibrated material model parameters. In situ tensile tests were carried out on specially prepared specimens with incremental data collection for different inclusions distributed within the ferrite matrix. A images were processed using digital image correlation-based tool to obtain microstrain measurements. The damage initiation, propagation, and strain localization around a large MnS inclusion were analyzed in detail. Other smaller inclusions tracked in this study helped in understanding the local material behavior. The conclusions of this extensive study are as follows: 1. The non-metallic inclusions are heterogeneously sized and heterogeneously distributed within the ferrite matrix. The inclusions are usually small (size~2 µm) and elliptical (aspect ratio <~2), with the exceptions of some extremely large (>10 µm up to 50 µm) and elongated (aspect ratio > 5) inclusions that are also present in the matrix. These large and elongated inclusions play a critical role in defining the limiting formability of the steel under consideration. 2. The 2D full phase simulation model developed in the current work provides accurate information about the material's local damage initiation and propagation under consideration. Although the same areas were not compared quantitatively, the simulation results match the experimental observations of global stress-strain response and the local damage initiation around the inclusions. Hence, the model can be used by engineers and researchers for further material engineering and optimization with confidence. 3. The local stress and strain largely depend on the local composition and distribution of non-metallic inclusions and the size and orientation of the neighboring ferrite grains. The local stress in highly clustered zones is~12% higher than the other material zones where the inclusions are more dispersed within the matrix. 4. The damage initiation and propagation also depend on the inclusion size and position.
If the inclusion is very large, brittle fracture occurs at relatively lower applied external stress (~450 MPa), which results in the fast damage initiation and propagation in the matrix. For most small and relatively elliptical inclusions, the damage initiates on the matrix/inclusion interface at relatively high-stress regimes (~550 MPa). It propagates at an oblique angle to the applied load. These smaller and relatively elliptical inclusions are also less prone to brittle cracking and strengthen the matrix during deformation. 5. An adequate material manufacturing methodology should be employed for the material class under consideration, resulting in small, elliptical, homogeneously distributed inclusions within the ferrite matrix. This would result in a material with better formability and higher damage resistance. resulting depressions ( Figure A2). The available DIC systems, such as VEDDAC, GOM, and VIC-2D, have therefore not made a flat closed evaluation possible. Therefore, the procedure with balls was not pursued in further investigations. The focus was placed on a special etching of the specimen surface with the application of the "InLens" module. The method for the deposition of SiO 2 nanoparticles on the specimen surface was as follows: 1. The surface of the specimen was metallographically prepared for colloidal silicon dioxide (OPS); 2. On a clean polished dip plate, 10-20 drops of OPS were placed along the diameters; 3. The specimen was pressed onto the cloth with 1-2 MPa pressure and turned for 1 s at 2 rpm; 4. With water flushing and 100 RPS plate speed, the specimen was rotated on the cloth at 1-2 MPa pressure at 2 RPS for 3 s; 5. The specimen surface was rinsed with ethanol and dried with a blower.
etching, it was difficult to achieve an even distribution on the surface, as all the beads accumulated in the resulting depressions ( Figure A2). The available DIC systems, such as VEDDAC, GOM, and VIC-2D, have therefore not made a flat closed evaluation possible. Therefore, the procedure with balls was not pursued in further investigations. The focus was placed on a special etching of the specimen surface with the application of the "InLens" module.
The method for the deposition of SiO2 nanoparticles on the specimen surface was as follows: 1. The surface of the specimen was metallographically prepared for colloidal silicon dioxide (OPS); 2. On a clean polished dip plate, 10-20 drops of OPS were placed along the diameters; 3. The specimen was pressed onto the cloth with 1-2 MPa pressure and turned for 1 s at 2 rpm; 4. With water flushing and 100 RPS plate speed, the specimen was rotated on the cloth at 1-2 MPa pressure at 2 RPS for 3 s; 5. The specimen surface was rinsed with ethanol and dried with a blower. As part of this work, various tests were also carried out to create contrasts on the surface of the specimens. To achieve the desired contrast, different etching times of 2 to 15 s were used ( Figure A3). As part of this work, various tests were also carried out to create contrasts on the surface of the specimens. To achieve the desired contrast, different etching times of 2 to 15 s were used ( Figure A3).

Appendix B. Local Mises Strain Measurement
The von Mises strain at the meso plane was calculated using a scale in the areas studied ( Figure A6). On each surface of the undeformed specimen, three points were selected so that one point was common; they formed two identical straights, and there was an angle of 90° between these straights. As a result of the tests, the meso deformation in the long inclusion in the tensile specimen was calculated. The results of meso-scale deformation are presented in Table  A1.

Appendix B. Local Mises Strain Measurement
The von Mises strain at the meso plane was calculated using a scale in the areas studied ( Figure A6). On each surface of the undeformed specimen, three points were selected so that one point was common; they formed two identical straights, and there was an angle of 90 • between these straights.

Appendix B. Local Mises Strain Measurement
The von Mises strain at the meso plane was calculated using a sca studied ( Figure A6). On each surface of the undeformed specimen, thre selected so that one point was common; they formed two identical straig was an angle of 90° between these straights. The strain was calculated using the following formula: As a result of the tests, the meso deformation in the long inclusion in the tensile specimen was calculated. The results of meso-scale deformation are presented in Table A1. The tests yielded the stress-strain diagram shown in Figure A7. The diagram shows that the range of plastic deformation started at 300 MPa. The maximum tensile strength achieved was 462 Mpa, and the breakage of the specimen began at a stretch of 7.9% Table A1. Change in the length of the account of the selected area and angle at the meso plane for the tensile specimen. The tests yielded the stress-strain diagram shown in Figure A7. The diagram shows that the range of plastic deformation started at 300 MPa. The maximum tensile strength achieved was 462 Mpa, and the breakage of the specimen began at a stretch of 7.9% Based on the data obtained, diagrams were created to compare the meso deformation and micro deformation in different areas of the study with the names of the selected areas ( Figure A8). A diagram has also been created to compare macro and micro deformation ( Figure A9). Based on the data obtained, diagrams were created to compare the meso deformation and micro deformation in different areas of the study with the names of the selected areas ( Figure A8). A diagram has also been created to compare macro and micro deformation ( Figure A9). The diagram in Figure A9 shows that the intensity of micro deformation is significantly higher compared to macro deformation at the points far from the pearlite grains. Figure A9. Comparison of meso deformation and micro deformation.

Appendix C. Local Strain Measurement Overlaid on SEM Micrographs
The frame by frame evolution of local strain distribution around (left) and within (right) the MnS inclusion is presented in Figure A10 from 1.6% to 7.75% global true strain. The diagram in Figure A9 shows that the intensity of micro deformation is significantly higher compared to macro deformation at the points far from the pearlite grains. The diagram in Figure A9 shows that the intensity of micro deformation is significantly higher compared to macro deformation at the points far from the pearlite grains.