Notch Stress Intensity Factor (NSIF)-Based Fatigue Design to Assess Cast Steel Porosity and Related Artiﬁcially Generated Imperfections

: Shrinkage porosities and non-metallic inclusions are common manufacturing process based defects that are present within cast materials. Conventional fatigue design recommendations, such as the FKM guideline (“Forschungskuratorium Maschinenbau”), therefore propose general safety factors for the fatigue assessment of cast structures. In fact, these factors mostly lead to oversized components and do not facilitate a lightweight design process. In this work, the effect of shrinkage porosities on the fatigue strength of defect-afﬂicted large-scale specimens manufactured from the cast steel G21Mn5 is studied by means of a notch stress intensity factor-based (NSIF-based) generalized Kitagawa diagram. Additionally, the mean stress sensitivity of the material is taken into account and establishes a load stress ratio enhanced diagram. Thereby, the fatigue assessment approach is performed by utilizing the defects sizes taken either from the fracture surface of the tested specimens or from non-destructive X-ray investigations. Additionally, a numerical algorithm invoking cellular automata, which enables the generation of artiﬁcial defects, is presented. Conclusively, a comparison to the results of the experimental investigations reveals a sound agreement to the generated spatial pore geometries. To sum up, the generalized Kitagawa diagram, as well as a concept utilizing artiﬁcially generated defects, is capable of assessing the local fatigue limit of cast steel G21Mn5 components and features the mapping of imperfection grades to their corresponding fatigue strength limit.


Introduction
Cast steel components are commonly utilized in modern manufacturing processes for complexly shaped structures and demands towards higher strengths. Despite elaborated efforts in terms of numerical process simulations during component design, casting defects, such as sand inclusions, hot tears and shrinkage porosities, are manufacturing based intrinsic material properties [1]. This is caused by the material's high shrinkage during cooling and solidification. Therefore, the fatigue assessment of cast steel defects is of utmost importance referring to an operationally endurable design.
Despite the use of radiograph standards [2][3][4] for quality assurance, which allow casting defects to be classified according to their accumulation and sizes, into defect and severity grades, it is difficult to attain a conservative estimation of the local strength properties in cast steel components. A fundamental relation between casting defects, including surface and near-surface discontinuities, on the fatigue strength and the impact behaviour of cast steels has been the topic of various preliminary studies [5][6][7][8]. Furthermore, the influence of shrinkage porosity on the elastic monotonic and fatigue life by computed tomography and X-ray scanning has been studied in depth by Hardin et al. [9][10][11][12][13][14], whereas finite-element simulations were utilized for a proper localized fatigue assessment. Next, Collini et al. [15] studied the influence of casting defects on the fatigue crack initiation of ductile cast iron concluding a significant drop in fatigue strength due to defects within the material.
In addition to these simulation-based fatigue assessment approaches, fracture mechanical concepts, like Murakami's √ area concept [16], are frequently used for the assessment of casting defects, especially in aluminium alloys, which exhibit shrinkage and gas porosities with sizes up to approximately 700 µm [17]. Recent works also depicted the applicability of Murakami's concept for non-metallic inclusion in cast steel materials, which are presented in [18][19][20]. Thereby it should be noted, that according to the literature [16], the application of the fracture mechanics-based concept provides confident results only for defect sizes with √ area ≤ 1000 µm. As presented by Hardin and Beckermann [21], shrinkage porosity can feature sizes up to several millimetres, which exceed the application limit of Murakami's √ area and therefore leads to an unsatisfying design result. On the other hand, as published in recent works [22,23], the generalized Kitagawa-Takahashi diagram, which can be considered as an extended form of the linear-elastic fracture mechanical design [24], enables an auspicious assessment of shrinkage porosity afflicted large-scale cast steel specimens. Early works by Williams, Gross and Mendelson [25,26] contributed to the development of the notch stress intensity factor (NSIF) concept in terms of an analytical solution of the angle-dependent Williams eigenvalue λ at sharp V-notches, as defined in Equation (1).
sin(λ 1,2 qπ) ± λ 1,2 sin(qπ) = 0 with q = 2π − ϕ π . (1) Thereby, ϕ represents the V-notch opening angle and q the relative complementarity angle. As the NSIF concept can be understood as a generalization of the fracture mechanical approach, the notch stress intensity factor at a V-notch under mode I loading can be formulated as given in Equation (2), by considering a gross-section stress ∆σ g as uniaxial load condition. Whereat γ represents the opening-angle dependent stress singularity at the notch root and is determined as given in γ = 1 − λ. In addition, the effective fatigue notch depth is defined by the factor a e f f , which is composed of the actual notch depth a and the geometry dependent factor α γ [27][28][29].
In case of a V-notch with an opening-angle equal to 0 • , the Williams eigenvalue results in a value of 0.5 and therefore leads to the fracture mechanical base equation [30,31].
The Kitagawa-Takahashi diagram (KTD) is basically specified by two material threshold parameters: the fatigue limit range ∆σ 0 , and the long crack stress intensity factor threshold range ∆K th,lc . The KTD features two different areas that can be subdivided into a region of non-propagating cracks and, on the contrary, into a region of flaws, which exhibits a size larger than the materials intrinsic crack length a ≥ a 0 . Thereon, the intrinsic crack grows under cyclic loading if the applied stress exceeds the materials threshold value, which is often ended by structural failure [24,32].
Concerning the generalized Kitagawa diagram (GKD), the corresponding values, including the long crack notch stress intensity factor threshold range ∆K V th,lc and the intrinsic El-Haddad length for V-notches a V 0 , have to be derived from the crack-related parameters ∆K th,lc and a 0 , respectively. According to Atzori et al. [28], ∆K V th,lc can be determined by means of by Equation (3) under consideration of the opening angle dependent stress singularity γ.
Relating to the determination of the intrinsic El-Haddad length for V-notches a V 0 , the Equations (2) and (3) can be utilized, where ∆σ g is substituted by the fatigue limit ∆σ 0 . By re-arranging the resulting formula, Equation (4) is easily obtained.
Representing an enhancement of the KTD, the generalized Kitagawa diagram additionally takes the opening-angle, in terms of the stress singularity, into account. Considering a V-notch in a double infinite plate exhibiting an opening angle ϕ greater than zero and a size of 2a de f ect , the local fatigue strength increases with an increase of the opening-angle, despite keeping a constant size, as it is illustrated in Figure 1. According to Atzori et al. [33][34][35], the Kitagawa-Takahashi can be separated into the regions of non-propagating short cracks, crack closure-afflicted long cracks and the notch mechanics. Additionally to the geometrical influences of the notch, the consideration of the mean stress sensitivity is of utmost importance for service design. Thereby it is necessary to distinguish between the mean stress effect concerning the fracture mechanical behaviour and the fatigue strength of the material. Relating to the fracture mechanical literature, several methodologies, such as the approach from Elber [36] and the corresponding adaptions by Schijve [37] and Walker [38] are known to take the load stress ratio (R-ratio) into account. Mann illustrates a comparison of the referred models, which can be found in [39]. Further approaches in regard to the mean stress influence on the crack growth are given in [40][41][42].
An analytical estimation of the crack opening function, see Equation (5), was proposed by Newman [43] based on finite-element simulations and subsequently correlated between the R-ratio and the long crack threshold value [44].
In accordance with Newman, the long crack threshold ∆K th,lc diminishes continuously as the R-ratio increases due to the steadily growing separation of the crack flanks in the contact field until the effective short crack threshold ∆K th,e f f is attained. Relating to elevated load ratios, from approximately R ≥ 0.7, crack closure will not occur [43,45,46].
Thereby, ∆K th,lc,R=0 represents the long crack stress intensity factor threshold range at a tumescent load ratio of R = 0. The factor γ N is the crack opening function, A 0 is the Newman coefficient and C th is a material constant. The applicability of Newman's approach, concerning the characterization of load ratio dependent crack growth thresholds, is illustrated in [44].
Referring to the influence of the mean stress on the fatigue strength, several approaches, such as by Goodman [47], Gerber [48], and Soderberg [49] have been proposed in preliminary studies. Figure 2 illustrates these commonly applied mean stress models, as well as the FKM guideline approach [50], which introduces a mean stress dependent sensitivity factor M σ .
An overview of these and further state-of-the-art common mean stress models can be found in [51]. The FKM approach features a reduced sensitivity at high load stress ratios and is therefore characterized by a bilinear formulation. Accordingly, the R-dependent fatigue strength amplitude can be determined by Equation (6).
where the mean stress sensitivity M σ is directly related to the material's ultimate tensile strength σ UTS by the material dependent parameters a M and b M , depicted in Equation (7).
A basic differentiation of mean stress models can be conducted based on the model's dependency. Since the approaches according to Goodman, Gerber, and Soderberg are related to the material's yield strength σ YS or rather the ultimate tensile strength σ UTS , the fatigue strength amplitude at a given mean stress can be determined with diminishing experimental effort. On the contrary, further state-of-the-art models [16,38,52,53] exhibit an additional material dependent parameter, which leads to a sound consistency between model and material behaviour, but they also require more extensive material testing.
The majority of casting defects are caused by the disturbance of the local feeding condition and therefore facilitate the formation of local shrinkage porosity, which can be considered as internal notches, leading to a significantly diminished fatigue strength. Furthermore, the influence of the mean stress on the long crack threshold and the material's fatigue strength in correspondence with macroscopic casting defects has to be taken into account. Since various models highlight the effect of the mean stress on homogeneous materials, the impairment of casting flaws has to be studied comprehensively. Hence, this work aims for the accessibility of mean stress exposed macroscopic shrinkage porosities and the corresponding estimation of the local fatigue limits based on computed tomography and X-ray investigations. Furthermore, focus is laid on the comparison between real casting defects and geometrically similar, artificially generated macroscopic pores.
Hence this paper scientifically contributes to the following topics: • Evaluation of the mean stress sensitivity of both homogeneous and defect afflicted cast steel material. • Fatigue strength assessment of large-scale cast steel specimens exhibiting macroscopic shrinkage porosities based on the NSIF approach applying the generalized Kitagawa diagram. • Calculation of the fatigue strength on the basis of X-ray examinations of cast specimens. • Numerical analysis of artificially generated macroscopic defects as basis for a probabilistic fatigue assessment of imperfections with varying, but cast steel process related, geometry.

Material and Fatigue Testing
The investigated macroscopic near-defect-free specimens are manufactured from the cast steel G21Mn5, which possesses a normalizing heat treatment. The chemical composition of the analysed samples is determined by a optical emission spectroscopy and tabulated in Table 1. Table 1. Chemical composition of the investigated material in weight percent.

G21Mn5 C [%] Mn [%] Si [%] Cr [%] P [%] S [%]
Fe In order to evaluate the grain size of the investigated material, five individual sections are metallographically prepared and chemically etched. A subsequent optical analysis of the microstructure obtained a grain size range of 30 to 40 µm. To apply the FKM mean stress sensitivity approach, quasi-static tensile tests are conducted in order to obtain the material's ultimate tensile strength. Figure 3 illustrates the utilized specimen geometry.

Fatigue Test Results
Aiming for the assessment of the materials near-defect-free long life fatigue strength, high-cycle fatigue tests under axial loading with constant amplitudes are conducted. In order to cover the material's mean stress sensitivity, the load stress ratios R = −1, R = 0, and R = 0.5 are selected. The fatigue tests are performed at a resonance test rig Power Swing MAG150kN (SincoTec, Clausthal-Zellerfeld, Germany) with a testing frequency of about 90 Hz at room temperature. The experiments are executed until either ten million load cycles are achieved or a specimen burst failure occurs. Figure 5 illustrates the utilized specimen geometry for the base material fatigue tests. In terms of the fatigue strength assessment, the statistical evaluation in the finite-life region is conducted by application of the common procedure ASTM-E739 [54]. Moreover, the endurance limit region is evaluated by means of the arcsin √ P methodology [55]. Further, the second slope k 2 is set to a value of 45 as suggested in [56,57]. Table 2 exhibits the stress ranges ∆σ, normalized to the fully-reversed long life fatigue value ∆σ 0,R=−1 , at a probability of survival P S of 50%, the evaluated finite-life slopes k 1 , the number of load cycles at the transition points N T and the corresponding scatter band indices 1:T S,E of the endurance limit region. Furthermore, the results are illustrated in Figures 6-8 by means of evaluated S-N curves, which show a decrease in the long life fatigue strength due to the elevated mean stress state. For the purpose of unity, the presented fatigue limits are normalized in regard to the fatigue limit ∆σ 0,R=−1 at a load stress ratio of R = −1. The load stress ratio R = 0 exhibits a reduction of the fatigue strength of about 21% compared to the fully-reversed R = −1 fatigue limit. Accordingly, the R = 0.5 curve depicts a more significant drop of the fatigue limit by approximately 42% due to the enhanced predominant tensile mean stress state compared to the value at R = −1.   The results of the experimental tests are subsequently compared to the estimated values of the mean stress approaches in terms of a Haigh diagram, see Figure 9. Additionally, Table 3 lists the resulting deviations ∆ of the fatigue strength at R = 0. The comparison exhibits, that Soderberg's concept underestimates the empirical data by nearly 24%. Contrarily, Gerber's method leads to an non-conservative estimation of the fatigue strength. As shown in Table 3, the mean stress approach of the common engineering FKM guideline provides an applicable estimation of the material's fatigue strength at R = 0 with a quite small deviation of 4%. Therefore, the assessment of the material's mean stress sensitivity can be conducted by means of the procedure provided in [50]. It should be noted that the applied FKM-guideline mean stress model is rather conservative for high tensile mean stresses as the measured yield to tensile strength ratio is quite low. Thus, for high stress ratios, Goodman's approach utilizing a linear relation might be more appropriate.  Table 3. Deviations ∆ between the estimated stress ranges by the different mean stress models and the experimentally evaluated data at R = 0.

Model [-]
In order to apply the chosen mean stress approach on macroscopic defect afflicted cast steel components, large-scale specimens are developed by means of numerical simulations. Therefore, the casting simulation tool Magma 5 © is utilized for the design process of the representative specimen geometry. In the following, a brief summary about the design process is given, for details see [23]. At first, the specimen's dimensions are set to a diameter of D = 30 mm and a total length of L = 300 mm. The location of the intended porosity is in the centre of the specimen and characterized by a larger diameter of 52 mm and constricted by narrow sections on both sides to induce a local undersaturation. This facilitates the formation of macroscopic shrinkage porosities. The moulding geometry is depicted in Figure 10, whereas the final specimen geometry is illustrated by Figure 11.  Subsequent to the casting of the specimens, a normalizing heat treatment is conducted to ensure a consistent microstructure equivalent to the small-scale specimens. In order to evaluate the mould design and casting process, non-destructive tests on the large-scale specimen are subsequently carried out by means of X-ray and computed tomography (CT) examinations. Figures 12 and 13 illustrate the CT and X-ray examination of a macroscopic defect at centre of specimen #30. The defect exhibits spatial dimensions of 15 mm in length and about 4 mm in diameter. Furthermore, the total lengths of selected defects, both derived from CT and X-ray, are listed in Table 4.  Subsequent to the non-destructive tests, the large-scale specimens are fatigue tested with constant amplitudes under varying load stress ratios, R = 0 and R = −1. The experimental test series are carried out on a resonance test rig MOT600kN (SincoTec, Clausthal-Zellerfeld, Germany) with an average test frequency of 80 Hz at room temperature. Thereby, the investigations are executed either until burst failure occurs or ten million load cycles are achieved, identically to the small-scale test settings. Run-out specimens are re-inserted and tested again at a higher load level to provoke a burst failure. Regarding the statistical assessment of the fatigue strength, the standard [54] is applied in the finite-life and the arcsine √ P approach [55] is utilized in the long life fatigue region. The results of the high-cycle fatigue tests considering both load stress ratios are illustrated in Figures 14 and 15 in terms of filled circles and triangles. Furthermore, the re-inserted specimens are depicted by framed circles and triangles, respectively. The stresses are again normalized in regard to the fully-reversed small-scale specimen long life fatigue limit of the plain material ∆σ 0,R=−1 . Since the prior conducted non-destructive test exhibited different enclosed defects in terms of their shapes, the fatigue results are differentiated between failure type A (FT#A), for the crack-like imperfections, and failure type B (FT#B) concerning the more globular pores. Despite numerical simulations and a reproducible execution of the casting process, the observed shape and size of the macroscopic defects vary from one to each other. Accordingly, the fatigue test results exhibit an enhanced scatter band, especially in the endurance limit region. Further on, fracture surface analysis of the investigated specimens are conducted in order to evaluate the sizes and shapes of the failure initiating defects. Thereby, two fundamental different types of flaws can be identified, referred as FT#A and FT#B, which are shown as blue triangles and green circles in Figures 14 and 15. The FT#A defect is characterized by a star-shaped appearance and a rectangular orientation concerning the specimen's axis. On the contrary, the FT#B specimens exhibit a more globular shaped defect as well as an axis parallel alignment. Both defect types are depicted by representative fracture surfaces in Figure 16. However, concerning the subsequent fatigue assessment, the defect size has to be determined at first. Hence, the area of the flaw is directly measured from the fracture surface of the investigated specimens and subsequently the equivalent circle diameter (ECD), exhibiting an equal defect area, is calculated. Moreover, to evaluate the effect of defects on the mean stress sensitivity, the results of both, the small-scale and large-scale specimens, are evaluated in Figure 17. The fully-reversed fatigue limit of the defect-afflicted specimens illustrates a reduction to 42% of the plain material fatigue strength. Considering the results at the R = 0 ratio, a relative drop from the plain fatigue limit σ 0,R=0 can be determined to be about 45%. Calculation of the mean stress sensitivities based on the experimental results of the small-scale specimens, representing a near-defect-free material, the large-scale defect-afflicted specimens, and the FKM guideline approach [50], utilizing only the ultimate tensile strength as characteristic material parameter, shows severe changes in the observed value. In detail, the engineering guideline provides a slightly reduced mean stress sensitivity compared to the base material; shrinkage pore afflicted large-scale specimens exhibit a significant decrease in mean stress sensitivity, as depicted in Table 5. Hence, it can be concluded that the experimentally determined cast steel mean sensitivity is in line with the guideline, but the mean stress sensitivity ought to be affected by the shrinkage porosities.  However, a mean stress state also influences the short-and long-crack propagation behaviour of the material and needs to be considered for a latter Kitagawa-Takahashi assessment. Therefore, experimental cyclic crack growth investigations have to be performed to study the material's fracture mechanical behaviour under varying load stress ratios.

Fatigue Crack Propagation Testing
In order to assess the material's fracture mechanical behaviour in terms of short-and long crack growth under different load stress ratios, fatigue crack propagation tests are conducted on the near-defect-free material. Therefore, single edge notched bending specimens (SENB) are manufactured from similar positions as the corresponding small-scale high cycle fatigue specimens. The utilized geometry of the investigated specimens is illustrated in Figure 18. The entire set of fracture mechanical specimens is prepared in accordance to the proposed procedure in [44]. At first, the V-shaped tip of the milled notch is razor blade polished to establish an even sharper notch at the notch root, as recommended in [58,59]. Secondly, the specimens are pre-cracked in terms of a cyclic compression-compression load at a load stress ratio of R = 20 with a maximum load level of ∆K = 18 MPa √ m until an incipient crack with a length between 100 and 200 µm is determined [60,61]. Consequently, residual tensile stresses are built-up in this plastic zone at the crack tip, provoking an emanating incipient crack exhibiting minimized crack closure effects [45]. Subsequent to the specimens preparation, cyclic crack growth tests are performed on the resonance testing machine Cracktronic © (Manufacturer: Russenberger Prüfmaschinen AG, City: Neuhausen am Rheinfall, Country: Switzerland) applying a constant amplitude load [61]. The measurement of the crack propagation is executed by the potential voltage drop methodology, as outlined in [62], in order to assess the emerging stress intensity factor during crack growth [63,64]. The results of the fracture mechanical investigations in terms of crack propagation curves for different load stress ratios are depicted in Figure 19. Further corresponding details are given in [65]. The empirical data of the R = −1 and R = 0 investigations could be tested up to the stable crack growth region. Contrarily, the R = 0.5 data region is truncated due to the upper force limit of the test rig. Nevertheless, the evaluation of the long crack stress intensity factor threshold values could be performed by utilizing the approach by Klesnil and Lukáš [66], see Equation (8).
The results illustrate the influence of the load stress ratio by means of a decreased ∆K th,lc at elevated R-ratios, which can be attributed to diminishing crack closure effects. In addition, Newman's approach, utilizing Equation (5), is applied on the experimental data, to evaluate the effect of the load stress ratio on the long crack stress intensity factor threshold values. Since no investigations under R-ratios lower than R = −1 are conducted, the long crack threshold is limited to the fully-reversed value to establish a conservative limit of the upper boundary. As depicted in Figure 20, the approach exhibits a sound agreement with the empirical data. Furthermore, Newman's approach estimates the long crack threshold at elevated R-ratios (R ≥ 0.7), which equals the effective threshold as ∆K th,e f f = 3.34 MPa √ m. The filled circles in Figure 20 correspond to the experimentally determined long crack thresholds of the crack propagation curves of Figure 19. Since, the influence of the mean stress on the fatigue limit and on the fracture mechanical threshold is determined for the investigated material, a generalized Kitagawa diagram (GKD) can be set-up to assess the fatigue strength of the porosity afflicted large-scale specimens in a unified approach.

Application of the Generalized Kitagawa Diagram
Firstly, in order to set up the GKD, the characteristic values ∆K V th,lc and a V 0 must be determined utilizing Equations (3) and (4). In order to take the mean stress sensitivity into account, the R-ratio dependencies according to the FKM guideline [50] and Newman's approach [43] are used.
Secondly, considering the assessment of the macroscopic shrinkage porosities, a relation between the defect size, referring to the equivalent circle diameter (ECD), and the effective defect size a e f f has to be established. Defect types, as depicted in Figure 16, are difficult to assess concerning the effect on the local fracture mechanical stress intensity factors due to their partially irregular shapes. The literature offers a wide variety of studied simple geometrical shapes and the corresponding fracture mechanical geometry factor solutions [67,68]. Referring to an engineering-feasible substitutional shape, a correlation between the arbitrarily-shaped defects and a corresponding penny-shaped crack has been studied in [23] by means of numerical crack growth simulations. In the following, the assessment approach is briefly reviewed. Details about the crack propagation rates, corresponding simulations and the results were given in [23].
The numerical analysis aims to link the arbitrary flaw ECD and the initial diameter of the penny-shaped crack in terms of a similar crack propagation lifetime. Therefore, several defect shapes of both defect types, FT#A and FT#B, are imported from the fracture surfaces into three-dimensional finite-element models, whereas the defects are simplified and integrated as two-dimensional cracks, considering an already initiated crack. The applied stresses in the simulation are in correspondence with the applied stresses of the experimental tests. The calculation of the local stress intensity factor ranges ∆K is conducted by the finite-element simulation tool Abaqus © (Version: 6.14). The subsequent determination of the local crack propagation rate da/dN around the perimeter is executed by means of the modified NASGRO equation [44]. The incremental simulations feature a crack growth of ∆a = 5 µm per step. Furthermore, the self-iterating crack growth model is stopped if the critical stress intensity range ∆K IC is detected at any point of the crack front perimeter.
A subsequent comparison of the numerical results exhibits a similar crack propagation rate between an FT#A defect and an ECD/4 penny-shaped crack. But for the FT#B defects, a penny-shaped crack of size ECD/2 results in a similar lifetime [23].
According to these findings, the defect ECDs, which are derived from the fracture surface analysis, are recalculated to obtain the effective defect sizes. Depending on the shape and orientation of the pores, the a e f f values are either determined by the factor ECD/2 or ECD/4 [23]. The results of the fatigue assessment are illustrated in Figures 21 and 22, where the effective defect sizes are obtained by fracture surface analysis. Thereby, the experimental load level of each specimen is plotted over the corresponding effective defect size. The aforementioned defect types are distinct by means of blue and green markers for FT#A and FT#B defects, respectively. Furthermore, tests that resulted in premature burst failure are depicted as '×' and run-outs are illustrated as '•'. The blue and green lines represent the long-life fatigue limits at ten million load cycles for V-notches with 45 • and 135 • opening-angles and are calculated in accordance with Equation (2). Additionally, it should be noted that the sketched fatigue strength design curves correspond to a probability of survival of 50%.  The calculated fracture mechanics results show a sound agreement for both opening angles. Thereby, a valid relation between the FT#A and the ϕ = 45 • line and between FT#B and the ϕ = 135 • line can be detected for both R-ratios. Consequently, a link between the spatial orientation of the defects and the opening-angle dependent GKD design curves is present in this context. The FT#A defects, which are rectangularly aligned regarding the specimen's axis and additionally exhibit a shallow depth, depict a crack similar geometry. The determined stress singularity of a 45 • V-notch results in γ = 0.495, calculated by Equation (1). As a comparison, the stress singularity of a crack results in a value of γ = 0.5. Hence, the assessment of shrinkage pores, which exhibit a crack like orientation and shape, can be satisfactorily conducted by a 45 • NSIF approximation based on these results.
Since a fracture surface based defect assessment is not feasible in an industrial manufacturing environment, the sizes of the defects are additionally non-destructively determined from X-ray results on specimens before cyclic testing. Therefore, the radiographic films from two rectangular directions are utilized for the measurement of the defect size. The average flaw size values are further utilized for the generalized Kitagawa diagram, as illustrated in Figure 23. Compared to Figure 22, which illustrates the GKD where the defect sizes are obtained by fracture surface analysis, a similar distribution concerning FT#A and FT#B failures and run-outs can be determined. Nevertheless, a reduced quantity of data points is depicted in the X-ray based diagram, due to a limited detectability of the porosities during the non-destructive testing. Hence, an X-ray based assessment enables equivalent results compared to a fracture surface based evaluation; however, the probability and reliability of detection represents the key aspects in this concern for engineering applicability. Due to the estimations based on Newman's approach and the fatigue limit model of the FKM guideline, the GKD can be enhanced towards the load stress ratio as the diagram's third axis. This enables a holistic defect assessment considering the spatial orientation, defect size and the present mean stress, as depicted by Figures 24 and 25.  Based on the overall shape, the derived GKD can be separated into two regions, represented by the red line as a common border. The area exhibiting defect sizes a e f f ≤ a V 0 is solely influenced by the Haigh diagram given by the FKM guideline. Contrarily, the area based on notch fracture mechanics, a e f f ≥ a V 0 , is additionally affected by Newman's approach. Therefore, the overall shape of the GKD surface is dependent on the applied load stress ratio, except the slope towards the defect size, which is itself entirely influenced by the stress singularity γ, hence the notch-opening angle. The red transition line exhibits a slight decrease of the a V 0 in the range −1 ≤ R ≤ 0. However, if the load stress ratio exceeds a value of R = 0, the transition knee point increases significantly to higher effective defect sizes. Furthermore, a considerable drop of the fatigue limit is obtained at load stress ratios R ≥ 0 due to elastic yield limitations resulting from the Haigh diagram.
Since the build-up of a probabilistic fatigue database, which consists of a large bulk of defects that exhibit varying shapes and spatial orientations, can substitute a high quantity of experimental effort, a numerical algorithm is further developed in order to generate such artificial porosities. Thereby, it is necessary to partially specify and limit the global shape of the flaws to keep the investigated geometries related to the investigated manufacturing process itself, but to also leave space for randomly artificially generated defects.

Assessment of Artificially Generated Defects
In order to establish artificial casting defects numerically, the algorithm of cellular automata (CA) is utilized. As already proposed in recent publications, the application of CA leads to satisfying results concerning the modelling of porosity formation and pore growth [69][70][71]. The presented results within this paper demonstrate the basic applicability of the developed tool, further details are provided in a preceding scientific work [72].
At first, a 50mm × 50 mm square grid, which acts as the plane basis of the porosity is generated. To construct a third dimension, 50 entities of this grid are equidistantly stacked on top of each other. Each cell of the generated cube can exhibit a state of either zero or one. For an initial set of the generated array, a random distribution of the available states is assigned. Thereby, a limiting hull curve is prescribed to constrain the basic shape of the porosity. Inside of the hull curve, the cells are able to be assigned either zero or one. Contrarily, outside the limiting curve, the individual cells are allocated with zero.
Secondly, the algorithm of the cellular automaton is applied under consideration of the three-dimensional Moore neighbourhood [73]. In order to avoid an excessive smoothing of the surface, the routine is applied in two subsequent runs.
Finally, the resulting array is converted into a point cloud, representing all cells that have the value of one assigned. Thereafter, the surface points are triangulated and smoothed. Figure 26 exemplarily illustrates a generated defect, which represents quite well the observed shrinkage porosity of FT#B within the investigated cast steel large-scale specimens. Moreover, the triangulated pore is subsequently imported into a finite-element model exhibiting the large-scale specimen's geometry. The model is meshed by means of quadratic tetrahedral elements and loaded with a uniform tensile stress along the model's z-axis.
In order to assess the numerically generated pores on the basis of the generalized Kitagawa diagram, the defect's ECD and an effective notch-opening angle has to be determined. The equivalent circle diameter is determined during the numerical generation of the pore by means of the projected defect area in the x-y plane. Relating to the effective notch-opening angle, an assessment of the vicinity around the maximum stress has to be conducted. Therefore, the nodal maximum principal stress and the corresponding normal surface vector has to obtained for the whole modelled defect surface. Subsequently, the defect is intersected with a plane that is spanned by the normal surface vector and a parallel vector, to the centre axis of the specimen. The obtained intersection curve of the pore, especially the section around the node of the maximum principal stress, is utilized for the assessment of an effective notch-opening angle. Hence, an empirical study of computed tomography measured imperfections fatigue data yielded that the length of the used secants should correspond to a value of ECD/12 for a proper assessment of the equivalent notch-opening angle ϕ as depicted in Figure 27. The resulting angle between both secants is subsequently evaluated and further utilized as the equivalent opening angle for the determination of the local NSIF value, applicable both for measured and artificially generated defects. Additionally, it should be mentioned that Neuber's rule may also be applicable once the defects and notch radii attain a certain size. MAX. In order to examine the sensitivity of the proposed approach, the numerical algorithm is utilized to generate several macro-geometry related porosities for each individually measured defect. Hereby, attention is paid to the occurrence of similar geometrical characteristics, as found at the computed tomography scans of the large-scale specimens. Further on, finite-element simulations, with a unified applied stress, are conducted in order to assess the local NSIF values and the corresponding fatigue limits, according to Equation (9). The resulting data is subsequently averaged and compared to the experimental test data as sensitivity qualifier, see Equation (10).
For illustration, three different porosities from the CT measurements are selected. Correspondingly, artificially numerical pores, which possess a similar overall shape and ECD, are generated. The selected casting defects and a corresponding numerically generated defect are sketched in Figure A1.
The results of the finite-element simulation and the subsequent calculation of the estimated fatigue limits as well as the equivalent circle diameters are summarized in Table 6. Referring to the results of the specimens #31 and #32, the estimated mean fatigue limits are both below the applied experimental stress ranges. In the case of specimen #31, the applied load was endured for about 9.4 × 10 6 load cycles, whereby ten million load cycles are regarded as run-out. In terms of specimen #32, the test was aborted due to burst failure at about 6 × 10 6 load cycles. Considering both estimated stress range limits, specimen #31 exhibits a smaller deviation to the applied stress range by 6 MPa than specimen #32 with 23 MPa. Concerning the third illustrated porosity, specimen #34, the experiment resulted in a run-out at ∆σ exp = 150 MPa, and further, to a burst failure at 175 MPa at about 1.3 × 10 6 load cycles. The accordingly generated defects and the accompanying fatigue assessment obtained an estimated fatigue stress range of 162 MPa, which is located between both experiments. Summing up, the numerical algorithm and the subsequent fatigue assessment deliver an averaged conservative deviation of 7.8%. Moreover, in [72] an extended study of ten specimens and their corresponding defects is performed, which results an average conservative deviation of 6%.
Furthermore, a mesh sensitivity study is conducted in order to evaluate the fatigue limit dependency on the element size. Therefore, the artificially generated porosity of specimen #31 is numerically evaluated by means of constant boundary conditions and an identical flaw shape, but exhibiting four different mesh sizes, including the average element sizes 500 µm, 250 µm, 167 µm and 125 µm as depicted in Figure 28. The results of the subsequent assessment is listed in Table 7.  Thereby it can be stated that the node of maximum principal stress is almost unchanged and thus independent from the mesh size. Contrarily, the refined mesh exhibits an increasing stress concentration factor, but the evaluated NSIF-value is only slightly affected. Referring to the evaluated fatigue limits, which are normalized to the result of the 500 µm model, the refined mesh study exhibits only small deviations. In detail, models with mesh seed of 250 µm, 167 µm and 125 µm lead to variances of 4%, 3% and 1%, respectively.
Concluding, an increasing mesh refinement results in a more accurate stress concentration factor. Furthermore, the study shows that the estimated fatigue limits are nearly unchanged. Additionally, it should be noted that a mesh refinement of the analysed imperfections lead to a significant increase of up to 400 times in terms of in calculation time.
Considering the results and a comparison of the experiments and numerically generated artificial pores, the NSIF-based approach offers sound applicability as an assessment tool. Furthermore, the utilization of artificially generated porosities provides a basic database for the probabilistic assessment of non-destructively detected casting flaws, which will be statistically covered in more detail in ongoing research work.

Conclusions
Based on the conducted investigations and presented results in this paper, the following scientific conclusions can be drawn: • The mean stress sensitivity approach of the FKM guideline provides a basically matching fatigue assessment relating to the cast G21Mn5 material. Referring to the defect-free base material on small-scale specimens, the guideline assesses a 16% reduced mean stress sensitivity. On the contrary, considering the large-scale defect-afflicted specimens, the FKM guideline estimates a 24% increased mean stress sensitivity. • The NSIF-based generalized Kitagawa diagram (GKD) represents an engineering feasible tool for the assessment of shrinkage porosities in cast steel for different types of defect shapes and orientations. • Furthermore, the fatigue assessment of shrinkage defects, by applying the GKD based on defect sizes taken from X-ray measurements, exhibits a satisfying agreement regarding the experimental results. It should be noted that it is essential to ensure the detectability of casting defects at X-ray investigations to facilitate a conservative fatigue design. Hence, it is recommendable to assume a minimum defect size according to the invoked NDT technique for engineering applications. • The enhancement of the generalized Kitagawa diagram towards the applied load stress ratio facilitates a unified fatigue assessment approach, which takes both the defect parameters and the loading conditions into account. • The presented numerical cellular automata algorithm, which enables the generation of artificial defects, shows a sound consistency with referenced experimental data leading to a conservative deviation of 7.8% and therefore offers ongoing research the set-up of a comprehensive defect database, consisting of different shapes, orientations, and sizes and even their corresponding estimated fatigue limits. Funding: The financial support by the Austrian Federal Ministry for Digital and Economic Affairs and the National Foundation for Research, Technology and Development is gratefully acknowledged.

Conflicts of Interest:
The authors declare no conflict of interest. The founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.