1. Introduction
Due to the continuously increasing requirements of our modern industrial market, focus is laid on economic and sustainable design. Lightweight engineering is often driven by topology optimization, yielding to complexly shaped components [
1]. Due to benefits in castability and comparable low density, AlSi-cast alloys are broadly used in mobility industry and transport applications [
2,
3,
4]. Besides Al-alloys, cast steel alloys play an important role as well due to improved structural load-bearing capability [
5,
6,
7]. The geometrical complexity of cast parts results in high local load factors. Moreover, strongly varying microstructural conditions may occur due to changes in local process parameters like local solidification [
8,
9,
10]. Such local microstructural conditions lead to varying shape and degree of porosity, affecting the local acting stress field significantly [
11,
12,
13]. This results in a severely limited fatigue strength of imperfective cast components compared to the strength of the ideal, defect free material [
2,
14,
15,
16,
17,
18,
19]. From the point of manufacturing quality, porosity can be subdivided into gas and shrinkage porosity according to [
20] and may be determined either destructively by fractographical inspection or non-destructively via computed tomography (CT) [
21,
22,
23].
Besides inhomogeneities, the effective mean stress state plays an important role in cyclic loading [
24,
25,
26,
27]. Although the global cyclic load predominantly results in elastic deformation, the aforementioned defects may lead to localized cyclic plasticisation and provoke therefore a change of the effective mean stress state [
28,
29]. In order to improve fatigue life of components, one can either reduce the external load, or improve the resistance of the material. If the material cannot be changed, a modification of the local strength by mechanical, chemical or thermal treatment may be helpful. For example, hardening or deep rolling of steel alloys leads to improvement of hardness and to compressive residual stresses depending on the process settings [
30,
31]. Shot peening yields an increased compressive residual stress state and hardness as well, but increases surface roughness [
26,
32,
33,
34,
35]. Moreover, welded structures can be treated by HFMI to improve fatigue strength due to local work hardening and compressive residual stresses [
36,
37]. In contrast, surface treatment by electrochemical polishing leads to a high surface quality but does not evoke residual stresses [
38]. Recent additive manufacturing technologies (3D/4D printing) feature multi-material composite design enabling tailored layer-based properties [
39]. Residual stresses can be induced by other post treatment processes as well [
40,
41,
42,
43,
44,
45,
46]. In detail, grinding and polishing by vibratory finishing as post-manufacturing process can be performed on simple fatigue specimens, but is also well applicable to complex parts like turbines blades [
47,
48,
49]. Vibratory finishing as a post manufacturing process leads to high reproducibility as well as a good surface quality with low scatter of surface roughness values as function of the process settings [
50,
51,
52]. But the fatigue design may be affected by this highly reproductive finishing processes due to the evocation of compressive residual stresses in the surface layer [
47,
53,
54] and a change of surface layer hardness due to small plastic deformations [
54,
55], depending on the material as well as the polishing medium. These effects of vibratory processing lead often to an increase of the long life fatigue strength [
56], but may also reduce fatigue strength [
57] depending on the effective residual stress state in the surface layer. Such contrary statements regarding the effects of post-manufacturing surface treatments highlight the importance of considering the effective mean stress state within the most-stressed layer properly.
Regarding nominal fatigue assessment of imperfective material, the common engineering guideline FKM (Forschungskuratorium Maschinenbau) recommends safety factors depending on the methods of non-destructive testing during quality control, but does not consider the statistical distribution of defects within the components [
58]. Surface treatment processes are incorporated as surface-layer-factor for post-treatments like hardening, deep rolling or shot peening for example, but the effective residual stress state is not taken into account in detail within this guideline. Application of such general safety factors leads to an underestimation of the fatigue strength, and therefore to a contradiction to modern light-weight design goals.
To assess the fatigue strength of notches or defects individually, the Theory of Critical Distances (TCD) by Taylor offers a widely used local methodology [
59,
60]. Within this approach, three methods are formulated, namely point, line and area methods. The point method states that the linear-elastic stress at a certain distance (critical distance) from the hot-spot is equal the fatigue strength of plain specimens
. El-Haddad et al. [
61] introduced the material intrinsic length
as a function of
and the long crack threshold
; see Equation (
1). Taylor’s critical distance is proportional to the intrinsic length of El-Haddad, as stated in [
59]. The proposed methodology was applied and validated in several studies [
62,
63,
64] with regard to notched components. Moreover, Taylor’s introduced methods were successfully applied considering the residual stress state in the surface layer of shot peened specimens [
32,
33]. Therein, sound predictions of the fatigue strength could be observed considering the residual stress state at a critical distance (point method) [
65], or by averaging the stress distribution over a distance from the surface (line method) [
32,
33,
34,
35].
In order to determine the cyclic stable stress state of ductile materials containing defects and residual stresses, a detailed knowledge about the elastic–plastic material behaviour is necessary. Due to the non-linear material response in parts loaded above a certain stress level, work hardening and subsequent mean stress relaxation occurs [
66,
67]. This implies that the cyclic loading hysteresis is shifting until stabilization of the material behaviour is reached [
68,
69]. As the distributed defects may vary severely in shape and size within the highly-stressed volume, the effective load stress due to this stabilization procedure of the hysteresis loop, the effective load–stress ratio
is introduced as a function of residual stress profile and external loading [
66,
70,
71]. The stabilization process of the residual stress profile in the surface layer has been well studied for welds [
72,
73,
74], shot peened parts [
75,
76,
77] or in the case of additive manufactured structures [
24].
In order to consider the non-linear cyclic deformation behaviour, a mathematical connection between stresses and strains is provided by the constitutive material equation [
78]. One common model to describe the elastic–plastic behaviour is provided by Chaboche [
79]. This combined hardening model includes an isotropic and kinematic part, and is based on a yield function
f, given in Equation (
2). Within this equation, the von Mises criterion is denoted by the term
, specified in Equation (
3) [
79].
In Equation (
3), the variable
describes the stress tensor and
is the kinematic back stress tensor;
and
are the related deviatoric components. The variable
k represents the initial size of the yield surface and the isotropic part
I of the combined hardening model describes the evolution of the size of the surface, as given in Equation (
4). Plastic deformation of the component takes place if the yield function reaches the limit value of zero, while
indicates a linear elastic material behaviour [
79].
The change of the size of the yield surface is described by isotropic hardening (see Equation (
4)) as a function of the accumulated plastic strain, but not referring to the position in the stress space. The constants
Q and
describe the difference between the initial and stabilized size of the asymptotic surface value and the stabilization rate as function of the accumulated plastic strain
p [
78,
80].
The translation of the yield surface in the stress space is given by the kinematic hardening model keeping the size of the yield surface constant, represented by the kinematic stress tensor
in Equation (
5) [
79]. For uniaxial loading, a non-linear approach can be expressed (see Equation (
6)) with the constants
and
[
79]. The accuracy of the model can be improved by superimposing higher order approaches
, as indicated by the variable m according to Equation (
5).
The initial conditions are represented by the parameters
and
in Equation (
6) and the value
describes the flow direction of the yield surface [
79].
and
representing the kinematic part as well as
Q and
b for the isotropic part can be experimentally determined based on low-cycle fatigue (LCF) tests using an optimisation routine, proposed in [
81].
For ductile metal structures of engineering components containing manufacturing-process-based inhomogeneities, fracture mechanics approaches feature reliable design tools to estimate the fatigue strength [
3,
82,
83,
84,
85,
86,
87,
88]. According to Murakami [
21], defects can be treated as cracks with the size
of the projected convex hull independent from their morphology. The influence of the defect’s location is introduced by means of a geometry factor
Y [
89]. Kitagawa and Takahashi [
90] introduced a threshold based methodology to describe the endurance fatigue strength as a function of the defect size, Equation (
7). Related extensions of the Kitagawa–Takahashi diagram (KTD) according to El-Haddad introduce a smooth transition between long and short crack behaviour, and Chappetti’s extension incorporates the formation of crack closure effects respectively [
61,
91]. Maierhofer et al. [
92] recommends describing the crack extension from the intrinsic short crack threshold
up to the long crack threshold
by Equation (
8) as functions of the predominant closure mechanisms
i. Regarding the in-depth influence of microstructure on crack growth parameters and the effect of the highly stressed volume, the authors refer to previous investigations [
93,
94] for Al–Si cast alloys.
with
Moreover, the threshold based tool of Kitagawa and Takahashi was applied in [
24,
95] for superimposing load and residual stresses and is quite effective for fatigue strength assessment considering residual stress measurements. Several studies conclude that residual stresses play a key role in affecting fatigue strength, and respectively crack propagation, in the high cycle fatigue (HCF) region [
36,
96,
97,
98,
99]. By superimposing the cyclically applied load and residual stress components, an effective load–stress ratio arises, which can strongly differ from the nominal load–stress ratio. This effective local load–stress ratio directly affects the local crack growth behaviour as the effective residual stress state opens or closes the crack flanks [
100,
101,
102]. To consider the influence of the load–stress ratio on the long crack threshold values, Newman [
103] introduced an engineering-feasible approach based on an initial crack threshold value
at tumescent load and a constant
, applied for steel in [
92] and for aluminium alloys in [
93]. In case of a defect free material, representing the upper border of the KTD, the effect of the load–stress ratio can be considered using the nominal stress approach of the common engineering guideline [
58]. The fatigue strength of the defect-free material
can be obtained experimentally using hot isostatic pressed (HIP) material, or otherwise estimated based on low cycle fatigue (LCF) tests [
104,
105]; see Equation (
9).
The size of the fracture-initiating defects, representing the axis of abscissae in the KTD, can be described by extreme value statistics, such as the generalized extreme value distribution (GEV) (see Equation (
10)) or the derived distributions according to Tiryakioğlu [
106]. Extreme value statistics are proved well suited for assessment of extremal defects in components, as summarized by Beretta [
107]. But one remaining challenge regarding assessment of extremal defects using computed tomography is the selection of a proper threshold value, which was investigated by Romano et al. [
108] using mean excess plots for peak over threshold (POT) sampling.
Besides the Kitagawa–Takahashi diagram, Murakami [
21] introduced the well-known
-concept as an engineering-feasible methodology for fatigue strength calculation of defect-afflicted components using the projected area of inhomogeneities. The original approach developed for steel alloys was extended by Noguchi et al. [
109] for aluminium alloys using the ratio of the Young’s modulus, as depicted in Equation (
11). The load–stress ratio is considered within the
-concept as a function of the value
, which depends on the hardness of the material. In case of defect free material, the base material strength can be calculated to
[
21].
with
This paper deals with the influence of vibratory finishing process and associated residual stress state on corresponding layer fatigue strength of three investigated cast materials, which covering steel and aluminium alloys with different manufacturing conditions. Five specimen series were taken from these materials, each of them in polished and vibratory finished condition to enable a direct comparison of the fatigue properties. After a detailed description of the experimental framework, the layer-based fatigue assessment concept, originally published by the authors in [
94] is extended with regard to the effective residual stress state within the surface layer evoked by vibratory finishing. The residual stress state in the surface layer is measured using X-ray diffraction. In order to determine the effective stable residual stress state after several load cycles, elastic–plastic simulations are set-up to analyse the material behaviour in the surface layer using a combined hardening model for each material. These results extend the probabilistic fatigue assessment methodology for aluminium alloy EN AC 46200 in T6 and HIP+T6 conditions. Subsequently, the extended methodology featuring a layer-based fatigue design with effective stress ratio dependency is validated for two cast steel alloys. Additionally, Murakami’s approach is applied to study the effect of increased surface micro-hardness by vibratory polishing. Summing up, this paper provides a detailed fatigue strength design method evaluating the effective stress ratio within the highly-stressed layers and thereby demonstrating the benefit of vibratory finishing as additional post-treatment for these cast alloys. If the fracture mechanical and fatigue strength parameters as well as imperfection distributions are known, the layer-based fatigue design method featuring the effective stress ratio
can be well applied to other materials and load-cases as well.
2. Materials and Methods
The investigated materials in this study are the cast aluminium alloy AlSi8Cu3 (EN AC 46200) in T6 and HIP+T6 heat treatment conditions and the cast steel alloys G12MnMo7-4+QT in quenched and tempered condition, and G21Mn5+N in normalized condition. The material AlSi8Cu3 was taken out of a gravity cast component which was manufactured by a core package casting process [
110]; nominal chemical composition is given in
Table 1 following the standard [
111]. After the casting process, the components were T6 heat treated to obtain a maximum in ultimate strength. Furthermore, some of the components were HIP+T6 treated to obtain a defect-free material batch. Two different positions were investigated within the components, namely position AlSi8Cu3-A possessing a secondary dendrite arm spacing of
µm in T6 and HIP+T6 condition and position AlSi8Cu3-B with
µm in T6 condition. The secondary dendrite arm spacing is not affected by the HIP processing, as the final T6 condition is identical for both manufacturing routes. Tensile tests revealed a yield strength of
MPa and tensile strength of
MPa for position AlSi8Cu3-A, and a yield strength of
MPa and tensile strength of
MPa for position AlSi8Cu3-B. It is important to note that the yield strength of
was evaluated at a strain value of 0.01% due to the brittle behaviour of position AlSi8Cu3-B. An overview of the dendritic microstructure in position AlSi8Cu3-A is given in the left column of
Figure 1 in the polished state of the metallographic sample.
A dedicated casting geometry was used for the material G21Mn5+N and G12MnMo7-4+QT to ensure a low degree of shrinkage porosity similar to the standard [
112]. The chemical composition of the material G21Mn5+N in normalized condition following the standard SEW 685 [
113] is given in
Table 2, as well as for material G12MnMo7-4+QT according to the standard SEW 520 [
114].
Quasi-static tensile tests for the material G21Mn5 revealed a upper yield strength of
MPa, a lower yield strength of
MPa, an ultimate tensile strength of
MPa and an elongation at fracture of more than 50%. An overview of the dendritic, ferrite–pearlite microstructure is given in
Figure 1 in the right column. The results of the quasistatic tensile tests of the high-strength cast steel alloy G12MnMo7-4+QT are given in a preceding work of the authors [
115]; the corresponding microstructure as a result of the heat treatment is given in the middle column in
Figure 1. The metallographic specimens of the cast steel alloys were grinded and polished, followed by etching with with 3% Nital. After the casting process, all cast parts were pre cut using water jet cutting and finally turned by a CNC machine. The final post-processing of the specimen surfaces, e.g., either manually polished or vibratory finished, is applied as final manufacturing step for all fatigue test samples.
2.1. Vibratory Finishing
It is well known that fatigue specimens have to be polished after the CNC-turning process to comply with the specifications of the standard [
116,
117]. This process is typically carried out manually, but in order to reduce the human influence for highly reproducible surface properties, vibratory finishing appears to be a useful alternative. In addition, electrochemical polishing is an alternative to manual polishing, but requires additional etching equipment and high-current power sources [
38]. Whereas polishing primarily targets a reduction of the surface roughness, vibratory finishing offers as side effect a modification of the surface layer regarding residual stress state and hardness distribution, as introduced in
Section 1.
A trough vibrator AVAtec
® TV 60 (Avatec Gmbh, Sindelfingen, Germany) was used to obtain a defined residual stress state by vibratory finishing. Therefore, different abrasive ceramic and plastic media were used in three subsequent treatment steps (strong grinding, smoothing and polishing). Additionally, a corrosion-preventing compound is necessary to absorb the abrasion of grinding media and workpieces. The specimens were grinded for several hours in the oscillating bowl for each abrasive media.
Figure 2 illustrates the obtained surface quality for the polished and vibratory finished condition of each series. Roughness parameters were evaluated with a Mahr
® MarSurf VD140 (Mahr Austria GmbH, Vienna, Austria). The arithmetic mean roughness value is
µm for the polished (POL) specimens and
µm for the vibratory finished (VF) ones, and the total height of the roughness profile is
µm for polished and
µm for vibratory finished specimens. Due to this quite high surface quality and the low difference of the roughness values, the effect of roughness on fatigue strength was purposely not taken into account in this study. For further sample designation, polished specimens are labelled with the suffix POL and vibratory finished ones with the suffix VF.
2.2. High Cycle Fatigue Tests (HCF)
The long life fatigue strength
for the investigated specimen series has been determined by uniaxial fatigue testing at room temperature following the standard [
118,
119]. The force-controlled high-cycle fatigue tests are conducted on a RUMUL
® Testonic 150 kN (Russenberger Prüfmaschinen AG, Neuhausen am Rheinfall, Switzerland) resonance test rig for the steel specimens and a RUMUL
® Microtron 20 kN (Russenberger Prüfmaschinen AG, Neuhausen am Rheinfall, Switzerland) for the aluminium specimens, imposing a global nominal load ratio of
.
Figure 3a,b depict the specimen geometries used for high cycle fatigue (HCF) testing, which are manufactured from bulk material. The blue marked test sections in
Figure 3 were finished with the respective surface treatment, either manually polished or vibratory finished. Each test terminated with the final rupture of the specimen, respectively a certain drop in frequency, or if the run-out number of
cycles is achieved. In absence of any damage, in terms of observable surface crack initiation, the run-outs were subsequently retested at high stress levels in the finite life region of the S/N-curve to study the fracture initiating defects for all fatigue samples.
Table 3 gives an overview of the high cycle fatigue tests with the related manufacturing and post-treatment coding. It should be mentioned that a part of these summarized results was obtained in preceding work by the authors [
115,
120,
121].
2.3. Low Cycle Fatigue Tests (LCF)
Low cycle fatigue experiments were conducted for the investigated cast materials to determine the elastic–plastic behaviour. The specimens were manufactured following the standard [
116,
122] featuring only a polished surface condition; see
Figure 3c for steel alloys and
Figure 3d for aluminium alloy. A servo-hydraulic test rig was used for the strain-controlled fatigue tests utilizing an extensometer for strain measurement. The specimens are tested at a triangular loading sequence featuring a constant strain rate of 1% per second at an alternate strain ratio of
.
2.4. Residual Stress Measurements
Polished and vibratory finished specimens of both steel series were subjected to residual stress measurements conducted with a X-RAYBOT from MRX-RAYS
® (MRX-RAYS, Brumath, France). This portable X-ray diffractometer was equipped with a 2 mm collimator in accordance with the standard DIN EN 15305 [
123], limiting the radius of the irradiated area to the half radius of the surface curvature. The aluminium series were measured with a Seifert XRD Charon SXL (XRD Eigenmann GmbH, Schnaittach, Germany) by Materials Center Leoben Forschung GmbH with a 0.3 mm collimator. Cr-K
radiation was used for all series with a wavelength of
Å and a
-mounting measurement setup following the standard [
124]. Furthermore, the applied measurement procedure is in accordance with the standard ASTM E2860-12 [
125]. A V-foil is placed in the secondary optical path to filter the Cr-radiation. The X-ray tube is operating at 20 kV and 1 mA and tilts in a range of ±40° in axial specimen direction (
-angle). A total number of 17 increments with an exposure time of 60 s each point is used. For the residual stress determination, the diffraction peaks from Al-alpha 311 and Fe-alpha 211 were used and the detected peak shape is fitted by means of a Pseudo-Voight function. The stress evaluation is performed based on the
method and elliptical fitting of the data points at a confidence level of 90% is applied. All residual stresses (surface and gradient in depth) were measured along the specimen axis, which represents the loading direction, at selected samples of each material batch.
To carefully remove small layers from the surface a Struers
® LectroPol-5 (Struers GmbH, Willich, Germany) was used for electrolytic polishing, exemplarily depicted as polished notch at the specimen surface in
Figure 4. A 3D-printed, flow-optimized fluid directing mask characterised by a 0.5
circular area was used at a voltage of 36 V for polishing in the test-cross section of the specimen. The electrolytic solution A2 of Struers is used at a temperature of 5 °C up to 10 °C to perform a constant material removal rate. The electrochemical erosion time was set for each increment individually, guaranteeing that the target depth was reached for subsequent residual stress measurements. The depth profile was measured by a Mahr
® MarSurf VD140 (Mahr Austria GmbH, Vienna, Austria) using the contour measuring unit, yielding a depth profile as depicted in
Figure 4.
2.5. Fracture Analysis and Computed Tomography
Subsequently to the fatigue testing, the ruptured specimens’ fracture surfaces were analysed utilizing a digital optical microscope (DOM) from KEYENCE® VHX-5000 (Keyence Corporation, Osaka, Japan) for macroscopic inspection. Local analysis of the fracture initiating defect was carried out by means of a ZEISS EVO 15 (Carl Zeiss AG, Oberkochen, Germany) scanning electron microscope (SEM), allowing a further magnified inspection of fracture origin and measurement of defects’ sizes. Non-destructive testing by means of computed tomography (CT) was realized by a Phoenix/X-ray Nanotom 180 (GE inspection technologies, Lewistown, PA, USA) using a voxel-size of just 4.5 m for series AlSi8Cu3-A-T6 and 5.5 m for AlSi8Cu3-B-T6.
2.6. Hardness Measurement
Macro hardness was measured by a ZWICK
® ZHU 2.5 (ZwickRoell GmbH & Co. KG, Ulm, Germany) test rig using a Vickers indenter. The applied testing parameters correspond to HV1 testing, according to standard ASTM E384 [
126], and were used for surface hardness measurements to differentiate between polished and vibratory finished condition. Furthermore, nano-indentation measurements were conducted at the Department of Materials Science (Montanuniversität Leoben, Leoben, Austria) for selected specimens to determine a work hardening of the surface layer by post treatment. Therefore, a Keysight Nanoindenter G200 (Keysight Technologies, Santa Rosa, CA, USA) testing unit utilizing a Berkovich indenter at a constant strain rate of 0.05
was used. The indentation depth was 150 nm for cross section in the surface layer and 2500 nm for the cross section in the bulk material.
2.7. Crack Propagation Testing
To build up the Kitagawa–Takahashi diagram and its fracture mechanical extensions, detailed information about the crack growth behaviour is necessary. Long-crack growth threshold stress intensity factors
were determined for load ratios from
up to
to obtain the Newman curve [
103]. The short crack threshold values
were determined at high load–stress ratios to minimize closure effects, as presented in in [
127,
128]. Single edge notched bending specimens were utilized. Their geometry is illustrated in
Figure 5.
After extracting the specimens from the bulk material using a CNC-milling machine, an initial notch featuring a depth of 4 mm and a width of 0.35 mm was manufactured using wire cut discharge machining (EDM). To facilitate the measurement of the initial notch geometry by the DOM, the sides of the sample were grinded and polished up to 1
m using diamond suspension. A sharp mechanical notch, which is necessary for subsequent compression pre-cracking, was produced by razor-blade polishing utilizing a 6
m diamond paste. Compressive pre-cracking [
129,
130] was conducted on a SincoTec
® Power Swing MOT100kN (SincoTec Test Systems GmbH, Clausthal-Zellerfeld, Germany) resonance test rig at a stress ratio of
. To determine the cyclic R-curve behaviour, the stepwise rising load amplitude crack growth test according to Tabernig et al. [
131] was applied. Crack propagation experiments were conducted using a RUMUL
® Cracktronic (Russenberger Prüfmaschinen AG, Neuhausen am Rheinfall, Switzerland) resonance test rig. The constant amplitude loading was started at
and stepwise increased to capture the cyclic R-curve until long-crack propagation
occurred. The crack length was measured by the direct current potential voltage drop (DCPD) method incorporating a temperature coefficient to compensate the effect of ambient temperature change.