Next Article in Journal
Τopology Optimization under a Single Displacement Constraint Using a Strain Energy Criterion
Next Article in Special Issue
Analytical Modeling for Mechanical Straightening Process of Case-Hardened Circular Shaft
Previous Article in Journal / Special Issue
Delamination Behavior of Highly Stretchable Soft Islands Multi-Layer Materials
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Numerically Efficient Method to Assess the Elastic–Plastic Strain Energy Density of Notched and Imperfective Cast Steel Components

Christian Doppler Laboratory for Manufacturing Process Based Component Design, Chair of Mechanical Engineering, Montanuniversität Leoben, Franz-Josef-Strasse 18, 8700 Leoben, Austria
*
Author to whom correspondence should be addressed.
Appl. Mech. 2023, 4(2), 528-566; https://doi.org/10.3390/applmech4020030
Submission received: 27 February 2023 / Revised: 29 March 2023 / Accepted: 10 April 2023 / Published: 27 April 2023
(This article belongs to the Special Issue Feature Papers in Material Mechanics)

Abstract

:
The fatigue strength of cast steel components is severely affected by manufacturing process-based bulk and surface imperfections. As these defect structures possess an arbitrary spatial shape, the utilization of local assessment methods is encouraged to design for service strength. This work applies the elastic–plastic strain energy density concept to study the fatigue strength properties of a high-strength cast steel alloy G12MnMo7-4+QT. A fatigue design limit curve is derived based on non-linear finite element analyses which merges experimental high-cycle fatigue results of unnotched and notched small-scale specimens tested at three different stress ratios into a unique narrow scatter band characterized by a scatter index of 1 : T Δ W ¯ ( t ) = 2.43 . A comparison to the linear–elastic assessment conducted in a preceding study reveals a significant improvement in prediction accuracy which is assigned to the consideration of the elastic–plastic material behaviour. In order to reduce computational effort, a novel approximation is presented which facilitates the calculation of the elastic–plastic strain energy density based on linear–elastic finite element results and Neuber’s concept. Validation of the assessment framework reveals a satisfying agreement to non-linear simulation results, showing an average root mean square deviation of only approximately eight percent in terms of total strain energy density. In order to study the effect of bulk and surface imperfections on the fatigue strength of cast steel components, defect-afflicted large-scale specimens are assessed by the presented elastic–plastic framework, yielding fatigue strength results which merge into the scatter band of the derived design limit curve. As the conducted fatigue assessment is based solely on linear–elastic two-dimensional simulations, the computational effort is substantially decreased. Within the present study, a reduction of approximately 400 times in computation time is observed. Hence, the established assessment framework presents an engineering-feasible method to evaluate the fatigue life of imperfective cast steel components based on rapid total strain energy density calculations.

1. Introduction

Due to strict environmental regulations, modern production strives for a lightweight design of mechanical components [1]. Engineering tools such as topology optimization represent feasible methods to combine application-tailored parts and economic production [2,3]. Such designs often lead to complex-shaped components possessing thin-walled cross-sections and therefore high loading conditions of their respective parts. In order to simultaneously fulfil the requirements in terms of manufacturability and fatigue strength of these complex structures, high-strength cast steel alloys are suitable engineering materials due to their favourable mouldability and high load-bearing capacities. The chemical composition of such alloys facilitates their utilization in welded hybrid constructions such as bogie frames in railway applications, consisting of complexly shaped cast components and standardized parts with various beam geometries. However, despite the latest advances in numerical simulations in the casting process [4,5,6,7,8], the formation of manufacturing process-based shrinkage imperfections within highly loaded sections cannot be fully prevented. Due to high-volume contraction, these spatial defect structures often reach dimensions of several millimetres [9,10,11,12]. Hence, the fatigue strength of the respective components is significantly affected by the presence of such macroscopic imperfections [13,14,15,16,17].
Design guidelines such as the commonly utilised guideline FKM (Forschungskuratorium Maschinenbau) [18] propose safety factors to consider the detrimental effect of defects on the fatigue strength of cast steel components. These safety factors only consider the applied methods of non-destructive testing during quality control of respective cast steel structures, fully neglecting the geometrical shape and location of inherent imperfections. Hence, the application of such safety factors often results in conservatively designed, overweight components which often repeal the advantages achieved by topology optimization. In order to take the geometrical shape of casting defects into account, the standard ASTM E446 [19] provides reference radiographs depicting representative imperfection geometries. Defect structures are classified into categories which reflect their causes of appearance, such as slag inclusions or shrinkage imperfections. Within a category, each reference radiograph represents a certain severity grade which is mainly defined based on the size and amount of the depicted imperfections. During the design process, specific combinations of defect category and severity grade are denoted as not acceptable for operation. Utilizing this defect classification system, assessment of defect-afflicted cast steel components is simply conducted by comparison of X-ray images of corresponding components to the reference radiographs. If occurring imperfections match a defect class denoted as non-acceptable, the component must be repaired or even scrapped and a new part has to be manufactured, increasing cost and time expenses significantly. A major drawback of this assessment method lies in the sole consideration of the defect geometry, fully neglecting its location in the assessed component. This circumstance implies that the existence of a certain imperfection could lead to the rejection of a manufactured component, even if it is located in a cross-section subjected to negligible stress. Furthermore, the interpretation of radiographs is a highly subjective task which depends strongly on the experience of the responsible personnel [20]. Hence, the same imperfection could be classified as either acceptable or non-acceptable, depending on the subjective interpretation of the non-destructive testing results.
In order to provide an objective quantitative assertion about the effect of imperfections on the fatigue behaviour of cast steel components, local assessment methods which consider the shape and location of inherent defect structures should be applied.
Defect structures such as shrinkage porosity networks induce a locally acting multiaxial stress state. Due to the ductile material behaviour of cast steel alloys, the resulting local loading conditions cause elastic–plastic deformation in the vicinity of each imperfection. Considering cyclic loading of notched or defect-afflicted parts, the non-linear material response results in mean stress relaxation [21,22,23,24] which alters the cyclic loading hysteresis until stabilization of the material behaviour. This stabilization process leads to the formation of a local effective stress ratio R e f f , similar to the case of residual stresses that interact with stresses induced by external loading [25,26]. In order to consider non-linear deformations during finite element analysis, the constitutive material equation provides the connection between strains and stresses [27]. Several elastic–plastic models are reported in the literature to include plasticity effects during cyclic loading such as Bauschinger’s effect or ratcheting [28,29]. A common constitutive model to describe the inelastic behaviour of metallic materials is given by Chaboche [30]. The approach presents a combined hardening model which considers isotropic as well as non-linear kinematic hardening and is therefore capable of describing the evolution of size and position of the yield surface in the space of stresses. The evaluation of the applicable constitutive model parameters is a well-researched topic and comprehensively published [31,32,33,34,35]. The application of combined hardening material models covers a broad variety of advanced research topics such as creep-fatigue analysis at elevated temperatures or uni- to multiaxial loading conditions in static and dynamic structural integrity analysis [36,37,38,39].
Dealing with the fatigue analysis of defect-afflicted ductile cast steel components, the non-linear deformation behaviour can be implemented into the assessment procedure to consider the cyclic evolution of the multiaxial stress state induced by the respective imperfections. This necessitates an assessment methodology capable of considering local plastic strains in order to accurately describe their influence on fatigue strength. The strain energy density (SED) concept presents a local fatigue assessment method that is utilised to assess arbitrarily shaped notches subjected to varying loading conditions. Fatigue failure occurs if the SED averaged over the control volume Ω reaches a critical value [40]. Herein, Ω represents a finite volume in the vicinity of a stress concentrator and is described by the material-dependent control radius R c and the notch geometry [41]. Figure 1 visualizes the definition of Ω for sharp and blunt V-notches, characterized by the notch opening angle 2 α and the notch root radius ρ [40,42].
In the case of a sharp notch, the control volume is centred in the notch root and takes the shape of a circle segment. Dealing with blunt notches as a generalized case, the control volume takes the form of a crescent shape with its centre located at a distance r 0 from the notch root, where r 0 is a function of 2 α and ρ [40,42].
In analogy to strains, the total cyclic strain energy density Δ W ( t ) can be defined as the superposition of an elastic Δ W ( e ) and plastic contribution Δ W ( p ) as given in Equation (1) [43].
Δ W ( t ) = Δ W ( e ) + Δ W ( p )
Equation (1) also holds true when dealing with SED values averaged over the control volume which are denoted by an overline accent, for example Δ W ¯ ( t ) .
Concerning linear–elastic material behaviour, analytical formulations for the calculation of the averaged SED Δ W ¯ ( e ) are reported for notched and unnotched specimens [42]. The elastic framework was successfully applied to describe the static and cyclic failure of notch-like geometries such as welded joints under uniaxial and multiaxial loading conditions [44,45,46,47,48,49,50,51,52,53,54]. More recently, the linear–elastic energy-based fatigue assessment of cast steel specimens afflicted by macroscopic shrinkage imperfections was presented, featuring the evaluation of the fatigue strength in the finite life, as well as in the endurance limit region. The calculations are based on finite element analyses of each arbitrarily shaped imperfection, as the analytical framework established for notched components cannot be directly applied in this case [10].
In order to improve the prediction accuracy during fatigue assessment, the energy dissipated by non-linear deformation will be included in the framework in terms of the plastic averaged SED Δ W ¯ ( p ) [55]. Regarding the contribution of the elastic SED Δ W ¯ ( e ) to the total SED Δ W ¯ ( t ) given in Equation (1), different approaches are reported in the literature:
  • Consideration of Δ W ¯ ( e ) regardless of its sign [56].
  • Equal contributions of Δ W ¯ ( p ) and only the positive part of Δ W ¯ ( e ) caused by tensile stresses [57,58].
  • Consideration of a weighting factor, which modifies the contribution of the positive part of Δ W ¯ ( e ) to the total strain energy density [59].
Benedetti et al. [60] summarized different approaches in a visual way, highlighting the absence of a general consensus about the relative importance of the elastic and plastic part in the calculation of the total strain energy density. In order to preserve consistency with the linear–elastic methodology reported in the literature [46,61,62,63], the elastic and plastic contributions Δ W ¯ ( e ) and Δ W ¯ ( p ) to the total SED Δ W ¯ ( t ) are equally considered following Equation (1), where Δ W ¯ ( e ) is included regardless of its sign. This approach was also referred to in a preceding study [10]. Figure 2 exemplifies the applied convention, depicting the individual parts of the SED for an illustrative stabilized hysteresis loop.
The total strain energy density concept was successfully applied to calculate the fatigue life of unnotched and notched components made of various materials under uniaxial and multiaxial loading conditions [64,65,66,67,68]. However, as an analytical calculation of Δ W ¯ ( p ) is only reported for unnotched specimens [43], the assessment procedure generally relies on numerical methods. Due to the non-linear material behaviour, computational efforts significantly rise, especially if the load history changes the hysteresis loop severely until stabilization.
Recalling established approximation concepts in notch fatigue, engineering-feasible methods by Neuber, Molski and Glinka [69,70] come to mind, which enable the approximation of the elastic–plastic stresses on the basis of linear–elastic results. Originally formulated for shear loading, Neuber [69] investigated the approximation of non-linear stresses and strains at the roots of sharp notches. A subsequent study by Seeger and Heuler confirmed the general applicability of Neuber’s concept on other loading configurations by utilizing strain measurements at the notch root [71]. The graphical interpretation of Neuber’s concept is given in Figure 3 and considers equality of the areas which are defined by the stresses and strains in the case of linear–elastic and elastic–plastic material behaviour. The respective parameters are highlighted by the superscripts ( e ) in the linear–elastic case and ( t ) in the elastic–plastic one.
Equality of the linear–elastic and elastic–plastic areas in Figure 3 is achieved if the points ( ε ( e ) , σ ( e ) ) and ( ε ( t ) , σ ( t ) ) are connected by a hyperbola [69]. This condition is expressed by Equation (2), considering the stress concentration factor K t , the nominal stress S in the net cross-section of the notched component and the Young’s Modulus E [72].
σ ( t ) · ε ( t ) = σ ( e ) · ε ( e ) = ( K t · S ) · K t · S E = K t · S 2 E
Due to the simplicity of Neuber’s concept in approximating elastic–plastic stresses and strains based on linear–elastic material behaviour, the methodology is applied in various research fields investigating local plasticity effects [73,74,75,76,77,78].
Molski and Glinka [70] presented an energy-based approach to calculate the elastic–plastic stress at the notch root. The introduced equivalent strain energy density concept assumes that the value of strain energy density in a notched body stressed by a specific load remains unaffected by the investigated material law. Hence, the strain energy density of the theoretical linear–elastic material and the real elastic–plastic one is the same [79,80]. In order to investigate the physical relationship between Neuber’s concept and the equivalent strain energy density method, Ye et al. [81,82] presented an approach based on the thermodynamic analysis of plastic deformation. The concept states that the plastic strain energy splits into heat energy and stored energy. Only the heat energy is dissipated during loading, while the remaining stored energy contributes to the local elastic–plastic stress and strain fields. The fraction of dissipated energy for one loading cycle is described by the so-called energy dissipation coefficient which constitutes the ratio of the dissipated heat energy to the plastic strain energy [82]. Furthermore, it was shown that Neuber’s concept presents a special case of the equivalent strain energy density approach, as the approach by Neuber neglects the cyclic plastic strain energy density when approximating the non-linear stresses. Hence, the concept presented in [81] yields elastic–plastic stress calculations resulting in values between the equivalent strain energy density concept as the lower limit and Neuber’s concept as the upper conservative limit. More recent studies deal with the formulation of pseudomaterials [83,84,85,86] as well as the approximation of the elastic–plastic stress and strain fields in the vicinity of crack tips [87].
However, as most fatigue studies reported in the literature focus on the assessment of notched components, the application of elastic–plastic approximation methods to arbitrarily shaped casting defects is a promising field of research. By evaluating the elastic–plastic stresses and strains induced by the complex geometry of such imperfections, the corresponding total strain energy density could be derived as given in Figure 2, considering local plasticity effects such as mean stress relaxation. Such elastic–plastic approximation methods significantly reduce computational effort and therefore feature a rapid fatigue design of defect-afflicted cast steel components.
Hence, this paper contributes to the following scientific topics:
  • Evaluation of an elastic–plastic design limit curve for the high-strength cast steel alloy G12MnMo7-4+QT.
  • Development and validation of a numerically efficient methodology to assess the elastic–plastic strain energy density based on elastic stresses and strains.
  • Application of the introduced rapid total strain energy density-based design to arbitrarily shaped, spatial bulk and surface imperfections found in cast steel components.

2. Materials and Methods

In order to establish the elastic–plastic assessment framework, the high-strength cast steel alloy G12MnMo7-4+QT in quenched and tempered heat treatment condition was chosen. This material is a promising alternative to weldable low-strength alloys such as G21Mn5+N to facilitate improved lightweight design. This holds especially true for large structures subjected to a high number of load cycles, as for example in railway components. The chemical composition of G12MnMo7-4+QT is given in the standard SEW 520 [88]. The high-cycle fatigue behaviour of the high-strength cast steel alloy was reported by axisymmetric, unnotched plain and notched small-scale specimens in a previous study [10]. Concerning the notched samples, V-notch opening angles of 2 α = 45 and 2 α = 135 were considered, where both geometries possess identical notch depths of 3 mm and notch root radii of ρ = 0.1 mm. The investigated stress ranges are referred to the cross-section of the V-notched specimens with a diameter of 16 mm. Concerning the plain specimens, stress ranges are referred to the diameter of 10 mm. Figure 4 presents the specimen geometries investigated in [10].
The high-cycle fatigue experiments were conducted on a RUMUL® Testronic 150 kN (Russenberger Prüfmaschinen AG, Neuhausen am Rheinfall, Switzerland) resonance test rig. The experimental setup is depicted in Figure 5.
All experimental investigations concerned uniaxial tension-compression loading at room temperature. Detailed information about the utilised number of small-scale specimens and the considered stress ratios are given in [10]. Additionally, this preceding paper includes high-cycle fatigue test results of large-scale specimens affected by macroscopic bulk shrinkage imperfections. It should be noted that only a linear–elastic energy-based approach was applied for fatigue assessment in this preceding study [10].
The present study reassessed these small- and large-scale specimen test results by means of an elastic–plastic energy-based approach. In order to characterize the non-linear material behaviour, low-cycle fatigue experiments were carried out. Furthermore, additional high-cycle fatigue experiments utilizing bulk defect-afflicted large-scale specimens and a stress ratio of R = 0.5 were conducted in the present study to investigate the influence of high mean stresses. The current experimental investigations also involved high-cycle fatigue experiments of surface defect-afflicted large-scale specimens in order to further enhance the experimental database of real-life casting imperfections. This extended dataset is used to demonstrate the uniform applicability of the formulated assessment framework. All experimental investigations were conducted at room temperature.

2.1. Low-Cycle Fatigue Investigations

The specimen geometry utilised in the conducted low-cycle fatigue experiments corresponds to the standard DIN EN 3988 [89] and is depicted in Figure 6.
The strain-controlled experiments were conducted on a servo-hydraulic test rig, utilizing an extensometer attached to the specimen for strain measurement. A triangular loading sequence was applied, featuring a constant strain rate of 1 % per second. The specimens were tested at an alternate strain ratio of R ε = 1 and subjected to total strain amplitudes between 0.18 % and 1 % . Figure 7 depicts the experimental setup utilised for the conducted low-cycle fatigue investigations.

2.2. High-Cycle Fatigue Investigations of Bulk Defect-Afflicted Specimens

To extend the existing experimental database reported in [10] towards higher stress ratios, additional high-cycle fatigue experiments of bulk defect-afflicted large-scale specimens made of G12MnMo7-4+QT were conducted. The utilised specimen geometry was identical to the one investigated in [10] and is presented in Figure 8. It should be noted that although the specimen geometry was the same, each sample was affected by an arbitrarily shaped bulk imperfection possessing a unique geometry.
X-ray investigations in two perpendicular planes, denoted as X- and Y-view, were carried out prior to testing in order to characterize the spatial shape of each individual bulk defect. The obtained radiographs were evaluated in terms of the calculation of geometrical shape parameters which are given in Table 1.
The experiments were carried out on a SincoTec® MOT600kN (SincoTec Test Systems GmbH, Clausthal-Zellerfeld, Germany) resonance test rig under constant amplitude and uniaxial tension-compression loading. Six specimens were tested utilizing a stress ratio of R = 0.5 and stress ranges between 100 MPa and 175 MPa at a testing frequency of approximately 75 Hz. The test strategy introduced in preceding studies [10,11,17] was followed for the experimental investigations. This strategy involved testing of each individual specimen at up to two target load levels. The first experiment was conducted at a stress range which was expected to result in a run-out of the specimen as the experiment reached ten million load cycles. In the following experiment, the load was increased in order to provoke burst failure of the re-inserted specimen. Hence, the testing procedure provided an experimentally derived assertion about the limiting impact each individual defect had on the fatigue strength. All investigated load levels refer to the cross-section of the specimens possessing a diameter of 30 mm. Figure 9 presents the experimental setup utilised for the conducted high-cycle fatigue investigations.
Subsequent to testing, the fracture surfaces were investigated using a KEYENCE® VHX-5000 (Keyence Corporation, Osaka, Japan) digital microscope. The conducted analysis revealed the bulk defects as the primary sites of fatigue crack initiation. Additionally, the fracture surfaces were geometrically evaluated in terms of shape parameters. Table 2 depicts the evaluated results and exemplifies the geometrical variety of investigated imperfections. As these defects were assessed in terms of fatigue strength by the framework which will be introduced in the following sections, the presented shape parameter results serve as a reference for possible future work.
In order to graphically demonstrate the variety in geometrical appearance of defect structures inherent in cast steel components, Figure 10 depicts the radiographs and fracture surfaces of three representative investigated bulk imperfections.

2.3. High-Cycle Fatigue Investigations of Surface Defect-Afflicted Specimens

In order to investigate the influence of macroscopic surface imperfections on the fatigue strength of the cast steel alloy G12MnMo7-4+QT, high-cycle fatigue investigations were conducted. The experiments utilize surface defect-afflicted large-scale specimens possessing a cross-section diameter of 32.5 mm, as depicted in Figure 11.
In order to ensure the surface defect as the site of fatigue crack initiation, a shape optimization study using Tosca® was conducted. The clamping conditions of the utilised test rig were modelled and the region around the clamping area was optimized aiming for a minimum stress concentration. Figure 12 depicts the result of the numerical study, revealing a maximum stress concentration factor of K t , m a x = 1.115 in the optimized specimen design.
The specimens were manufactured by the shell mould casting process. In order to achieve an arbitrary surface defect, a sand particle of several millimetres in size was inserted into the mould. During casting, the particle floated on top of the liquid metal due to the difference in density, featuring the intended macroscopic as-cast surface notch. Subsequent to casting, the specimens were machined in the clamping areas to achieve the geometry depicted in Figure 11. Non-destructive testing in terms of X-ray examinations was conducted to prove that no other sand particles were enclosed in the specimen cross-section which could act as possible sites of fatigue crack initiation. Furthermore, the spatial surface defect geometry of each specimen was measured non-destructively by a digital microscope. A three-dimensional image stack was captured for each defect, acquiring its spatial surface which served as an input parameter for the fatigue assessment procedure outlined in the following sections. Table 3 summarizes the minimum, maximum and mean values, as well as the standard deviation of the derived shape parameters, whereat 17 specimens were investigated in the present study.
Subsequent to the geometrical characterization of the defect structures, high-cycle fatigue experiments were conducted on the resonance test rig depicted in Figure 9, applying constant amplitude and uniaxial tension-compression loading. In order to investigate the effect of mean stress on the fatigue strength, the specimens were tested at stress ratios of R = 1 , R = 0 and R = 0.5 . The same test strategy involving subsequent testing of individual specimens at up to two load levels, as described in Section 2.2, was applied to surface imperfections samples. The conducted experiments featured stress ranges between 50 MPa and 325 MPa at an average testing frequency of approximately 84 Hz. It should be noted that the investigated load levels refer to the cross-section of the specimens with a diameter of 32.5 mm, as the definition of a net cross-section is complicated due to the arbitrary spatial shape of the surface imperfections.
After completion of the experiments, fracture surface analysis was carried out utilizing the digital microscope. The analysis confirmed the macroscopic surface imperfections as the primary site of fatigue crack initiation. Analogous to the optical assessment prior to the conducted high-cycle fatigue experiments, the geometrical shape parameters of the defects were evaluated on the fracture surfaces. Table 4 presents the derived minimum, maximum and mean values, as well as the standard deviations, providing further information about the spatial geometry of the investigated surface defects.
Figure 13 depicts three representative surface defect-afflicted specimens prior to testing, as well as the respective fracture surfaces investigated after burst failure.

3. Results and Discussion

As the experimental database covering low- and high-cycle investigations is concluded, evaluation of the non-linear material behaviour as well as the definition, validation and application of the energy-based fatigue assessment framework is studied next.

3.1. Characterization of the Low-Cycle Fatigue Strength

Evaluation of the experimental results obtained by the low-cycle investigations presented in Section 2.1 is conducted according to DIN EN 3988 [89], identifying the number of cycles until failure N f by a decrease in the measured stresses of 10 % subsequent to stabilization of the material behaviour. Furthermore, the elastic–plastic strains at half-life N f / 2 are evaluated and a strain-based fatigue assessment curve is derived according to Basquin [90], Manson [91] and Coffin [92]. Figure 14 depicts the evaluated Manson–Coffin curve.
The total strain amplitude ε a ( t ) is given by Equation (3) as a superposition of the elastic and plastic strains [90,91,92].
ε a ( t ) = ε a ( e ) + ε a ( p ) = σ f E · ( 2 N f ) b + ε f · ( 2 N f ) c
The elastic strain amplitude ε a ( e ) is described by the fatigue strength coefficient σ f , the stabilized cyclic Young’s Modulus E and Basquin’s fatigue strength exponent b. The plastic strain amplitude ε a ( p ) is evaluated utilizing the fatigue ductility coefficient ε f and Manson–Coffin’s fatigue ductility exponent c [90,92,93].
To describe ε a ( t ) in the stabilized material condition in terms of the acting stress σ a , the formulation according to Ramberg–Osgood given in Equation (4) is utilised [94].
ε a ( t ) = ε a ( e ) + ε a ( p ) = σ a E + σ a K 1 n
The cyclic strength coefficient K and the cyclic strain hardening exponent n featured in Equation (4) are related to the parameters of the Manson–Coffin curve, following Equations (5) and (6) [95,96].
K = σ f · ( ε f ) n
n = b c
Equation (4) can be visualized as a cyclic stress–strain curve, utilizing the experimental data in the stabilized material condition; see Figure 15. Herein, cyclic softening of the investigated cast steel alloy G12MnMo7-4+QT can be perceived clearly, as the initial yield stress determined by preceding quasi-static tensile tests possesses a value of σ Y S = 549 MPa [10].
Table 5 summarizes the parameters of the cyclic stress–strain curve.

3.2. Parametrization of Elastic–Plastic Material Model

The low-cycle fatigue experiments reported in Section 2.1 are conducted utilizing the plain specimen geometry depicted in Figure 6. Hence, the description of the elastic–plastic strain as a function of stress according to Equation (4) is only valid for unnotched specimens, where the deformation behaviour is uniform in the loaded cross-section. As notched specimens exhibit a stress gradient emanating from the notch root, local plastic deformation can occur in the highly-stressed vicinity of the notch root, while elastic deformation is predominant in regions further away from it. Dealing with cyclic loading, this local plastic deformation results in effects such as mean stress relaxation in the notch vicinity, leading to a complex evolution of the acting stress field which cannot be described in an analytical manner. Thus, numerical analyses are necessary utilizing a proper material model capable of describing the inelastic material behaviour.
In the present study, a combined hardening model considering isotropic and non-linear kinematic hardening is chosen to describe the non-linear material behaviour. The model is based on the definition of a yield function f given in Equation (7) [30].
f = J σ X I k
The term J σ X corresponds to the von Mises criterion and is defined by Equation (8) [30].
J σ X = 3 2 σ X : σ X
Herein, the parameter σ represents the stress tensor and X corresponds to the kinematic back stress tensor. The parameters σ and X denote the deviatoric parts of σ and X , respectively.
The yield function f reported in Equation (7) represents a surface in the space of stresses where plastic deformation of the stressed component takes place if the condition f = 0 is fulfilled. In the case of f < 0 , the material flow is of a linear–elastic type. The initial size of f at zero plastic strain is described by the parameter k. During cyclic plastic loading, the dislocation structure inherent in the material is altered, resulting in a change of the size and the position of the yield surface f. The evolution of the surface size is described by the isotropic part I of the combined hardening approach according to Equation (9) [97].
I = Q · 1 e b s p
The size of the yield surface changes as a function of the accumulated plastic strain p, where the parameter Q denotes the difference between the initial and stabilized size of the yield surface and b s represents the stabilization rate under constant strain amplitude loading [28,30].
The position and translation of the yield surface in the space of stresses are governed by the kinematic part of the combined hardening model [97]. In order to improve the description accuracy of the non-linear material model, the kinematic back stress tensor X is defined according to Equation (10) as a superposition of individual contributions X i [30].
X = i = 1 n X i
Each part X i follows a differential equation [98] which can be integrated to give Equation (11) for the uniaxial load case [28,30].
X i = ψ C i γ i + X 0 , i ψ C i γ i · e ψ γ i ε ( p ) ε 0 ( p )
The parameters X 0 , i and ε 0 ( p ) represent the initial conditions, where ε ( p ) denotes the plastic strain. The direction of the flow is described by ψ which takes values of ψ = ± 1 .
In the present study, the parameters of the combined hardening model, including the kinematic material parameters C i and γ i , as well as Q and b s of the isotropic part, are evaluated utilizing an optimization routine [31] which was developed at the Chair of Mechanical Engineering at the Montanuniversität Leoben. This routine was comprehensively validated for various materials under ambient and elevated temperatures [99,100,101,102,103]. In the present study, the material parameters featured in Equation (9) and (11) are determined based on the experimental data reported in Section 2.1. Concerning the non-linear kinematic part of the combined hardening model, two back stress components are considered. Table 6 summarizes the parameters of the elastic–plastic material model.
As the material parameters are determined, validation of the constitutional model is carried out by a non-linear finite element analysis in Abaqus® 2021. A single three-dimensional element of unity edge length is utilised to evaluate the stabilized hysteresis loop at a total strain amplitude of ε a ( t ) = 0.5 % . The element is of the type C3D20 and features quadratic shape functions. Figure 16 depicts the finite element model, as well as the numerically evaluated hysteresis loop for the stabilized material condition in contrast to the experimental results. In this case, a stabilization of the hysteresis loop has been obtained after approximately 1100 cycles.
The comparison of numerical and experimental results given in Figure 16 reveals a satisfying representation of the stabilized cyclic material behaviour by the established combined hardening model.

3.3. Evaluation of the Elastic–Plastic Design Limit Curve

As the non-linear material behaviour of the investigated cast steel alloy G12MnMo7-4+QT is characterized, the elastic–plastic design limit curve can be assessed by means of an energy-based approach. This limit curve is further utilised to calculate the fatigue strength of defect-afflicted components. In order to derive the elastic–plastic design limit curve, the experimental high-cycle fatigue results of small-scale specimens, reported in a preceding study [10], are reassessed.
The elastic–plastic design limit curve is evaluated based on the averaged total strain energy density Δ W ¯ ( t ) . Therefore, finite element analyses are conducted utilizing Abaqus® 2021. The non-linear material behaviour is considered by implementation of the combined hardening model presented in Section 3.2. Axisymmetric models are utilised to evaluate Δ W ¯ ( t ) of the respective plain and notched specimens, featuring elements of type CAX6/CAX8 with quadratic shape functions. The control volume is modelled by partitioning of the notched geometries, considering a control radius of R c = 0.083 mm, as evaluated in [10] for G12MnMo7-4+QT. Herein, the evaluation of R c was conducted following the linear–elastic framework given in [40,41], considering the fatigue strength and long-crack growth threshold stress intensity factor range which correspond to the fatigue limit at ten million load cycles. Three different stress ratios R have been observed, where R c is defined as the arithmetic mean value. In order to ensure appropriate element sizing within the control volume Ω , the results of a mesh sensitivity study conducted in a preceding study [104] are considered. The SED within Ω is rather insensitive to the element size. Hence, control volumes consisting of a minimum of three elements yield appropriate SED values without sacrificing accuracy, as also stated in [105]. The position r 0 of the control volume center point is derived according to [42] for the investigated notch opening angles of 2 α = 45 and 2 α = 135 . The numerical models are stressed under uniaxial tension-compression loading, covering the tested load amplitudes and stress ratios of the high-cycle fatigue experiments conducted in [10]. A cyclic load is applied during numerical analysis until stabilization of the material behaviour is achieved. Stabilization is assumed if the iterative change in stresses of each cycle is less than approximately one percent. This was obtained after approximately one hundred load cycles for the conducted simulations. The cyclic load in Abaqus® 2021 is realized by the definition of a time-dependent multiplier m R , altering the currently acting load which corresponds to the investigated stress range Δ σ . Figure 17 exemplarily depicts the finite element model of the investigated V-notched specimen geometry possessing a notch opening angle of 2 α = 45 . The presented stress result refers to the stabilized material condition at a stress range of Δ σ = 120 MPa and a stress ratio of R = 1 . The specimen is depicted for the maximum value of the numerical cyclic load. Furthermore, the course of the load multiplier m R utilised for the cyclic simulations is depicted for all investigated load ratios for the first two cycles.
Finally, the elastic and plastic strain energy contributions are evaluated from the control volume elements and the averaged total strain energy density in the stabilized material condition is calculated. Figure 18 summarizes the elastic–plastic fatigue assessment results.
Concerning the long-life fatigue strength, the presented results are statistically evaluated according to the a r c s i n P approach [106]. The finite-life fatigue regime is assessed following ASTM E739 [107], including the statistical evaluation of the fatigue limit curves for a probability of survival of P S = 50 % , as well as for values of P S = 10 % and P S = 90 % . Furthermore, the limit curve for P S = 97.5 % is depicted, which serves as an assessment basis in fatigue design, according to the widely used guideline FKM (Forschungskuratorium Maschinenbau) [18].
As depicted in Figure 18, the elastic–plastic fatigue assessment based on the total strain energy density summarizes the fatigue results of three different specimen geometries tested at three different stress ratios into a unique narrow scatter band. The scatter index which is defined by the limit curves of P S = 10 % and P S = 90 % possesses a value of 1 : T Δ W ¯ ( t ) = 2.43 . Comparing to the scatter index of 1 : T Δ W ¯ ( e ) = 4.32 as reported in [10] utilizing a linear–elastic SED approach, the consideration of the non-linear material behaviour by the total SED provides a significant improvement in prediction accuracy. As stated in [10], the comparably high value of 1 : T Δ W ¯ ( e ) is primarily assigned to the assessment results at a stress ratio of R = 0.5 . Due to the high maximum loads, inelastic material deformation and related mean stress relaxation occur, which cannot be considered by the linear–elastic framework. This leads to an overestimation of the local inherent deformation energy and therefore an increased scatter index when gathering the fatigue assessment results obtained at different stress ratios. The overestimation due to the linear–elastic formulation holds especially true when dealing with high provoked stress concentrations, as in the case of sharply notched components. To illustrate the localized plastic deformation effect and its associated mean stress relaxation, a comparison between the results of the linear–elastic and the elastic–plastic SED concept is presented in Figure 19. Both concepts are utilised to assess V-notched specimens with a notch opening angle of 2 α = 45 and a notch root radius of ρ = 0.1 mm. The linear–elastic SED Δ W ¯ ( e ) is calculated analytically utilizing the framework given in [42]. To consider the stress ratio in the calculation of Δ W ¯ ( e ) , the correction factor c w according to [61] is implemented. The total SED Δ W ¯ ( t ) is evaluated numerically utilizing the set up combined hardening model. The derived SED values are depicted over the linear–elastic maximum stress as a reference.
Reflecting the data points given in Figure 19, the effect of cyclic mean stress relaxation can be observed. As the theoretical linear–elastic maximum stress increases, the total SED and the purely linear–elastic SED diverge from each other due to occurring plastic deformation in the vicinity of the notch root. In order to describe the effect of stress ratio, the behaviour of a common load level Δ σ i is studied first. These calculations are highlighted by an additional bounding box in Figure 19. Δ σ i is chosen such that the deformation is mainly linear–elastic in the case of R = 1 . Hence, there is no significant difference between Δ W ¯ ( e ) and Δ W ¯ ( t ) . As the stress ratio and therefore the maximum stress increases, the mismatch between the elastic–plastic and the purely linear–elastic solution becomes more pronounced. This behaviour is explained by the occurring mean stress relaxation due to plastic deformation and the associated decrease in the locally acting stress ratio.
Second, the combined effect of stress range and stress ratio is investigated. Examining the depicted results in Figure 19, one can observe opposite trends of Δ W ¯ ( t ) compared to Δ W ¯ ( e ) as the stress range and stress ratio increase; while Δ W ¯ ( t ) takes higher values than Δ W ¯ ( e ) with an increasing stress range in the case of R = 1 , the calculated results of Δ W ¯ ( t ) decrease compared to Δ W ¯ ( e ) in the case of higher stress ratios. This is reasoned by the cyclic plastic deformation which contributes a significant amount to the total SED in the case of R = 1 . Considering R = 0 and R = 0.5 , the bearable stress ranges are relatively small compared to R = 1 , resulting in a mainly elastic cyclic deformation. Hence, the contribution of the cyclic plastic SED to Δ W ¯ ( t ) cannot compensate for the dominant reduction of the cyclic elastic SED term due to mean stress relaxation. This results in an overall decrease of Δ W ¯ ( t ) compared to Δ W ¯ ( e ) .
As demonstrated, the implementation of the non-linear material behaviour severely affects the local stress field and therefore the fatigue behaviour of cast steel components. However, the elastic–plastic finite element analyses necessary to evaluate the total strain energy density increase computational efforts significantly. In order to keep time and cost expenses at a minimum, an alternative evaluation methodology based on the acting linear–elastic stress and strain field is studied next.

3.4. A Numerically Efficient Approximation Approach to Evaluate the Elastic–Plastic Strain Energy Density

In order to derive an engineering-feasible methodology for the approximation of the averaged elastic–plastic strain energy density Δ W ¯ ( t ) at complexly shaped casting imperfections, the V-notched specimens are investigated firstly. The respective specimen geometries correspond to the ones utilised for evaluation of the elastic–plastic design limit curve presented in Section 3.3. As Δ W ¯ ( t ) is defined by the stabilized hysteresis loop, see Figure 2, the non-linear stress and strain fields have to be calculated within the finite control volume Ω in order to derive Δ W ¯ ( t ) . Due to the general multiaxiality of the local acting stress state, the presented methodology is formulated to include all components of the stress tensor σ and strain tensor ε given in Equations (12) and (13).
σ = σ 11 σ 12 σ 13 σ 21 σ 22 σ 23 σ 31 σ 32 σ 33
ε = ε 11 ε 12 ε 13 ε 21 ε 22 ε 23 ε 31 ε 32 ε 33
In the following, the components of σ and ε are denoted as σ i j and ε i j , where i , j { 1 , 2 , 3 } . For the sake of completeness, it should be noted that σ i j = σ j i and ε i j = ε j i , due to the symmetry of σ and ε .
The definition of the applied approximation method for a component loaded under an arbitrary stress range Δ σ and a stress ratio of R is demonstrated for V-notched specimens, possessing a notch opening angle of 2 α = 45 and a notch root radius of ρ = 0.1 mm. In the present study, the conservative concept according to Neuber [69] is chosen to evaluate the elastic–plastic stresses and strains in the control volume Ω . This necessitates the evaluation of the linear–elastic stress tensor σ k for each finite element C k within Ω . Therefore, a linear–elastic finite element analysis is carried out, utilizing a stress of 1 MPa in the cross-section with a diameter of 16 mm. Figure 20 schematically depicts the numerical model and the control volume of the investigated V-notched specimen geometry.
The stress components σ i j , k are evaluated in the centroid of the respective element C k . Due to the utilization of a stress of 1 MPa in the specimen cross-section, σ i j , k constitutes the elemental stress concentration factors for the individual directions in the space of stresses. In order to approximate the elastic–plastic stresses and strains for cyclic loading of the specimen, the cyclic stress–strain curve, as evaluated in Section 3.1, is utilised with some additional considerations. The low-cycle fatigue experiments conducted to determine the cyclic stress–strain curve utilised plain specimens loaded under uniaxial tension-compression loading. This implies the existence of a uniaxial stress state in the specimen cross-section. Hence, the cyclic stress–strain curve described by Equation (4) is only valid for a uniaxial stress state. As the evaluated stress tensor of element C k possesses in general more than one non-zero component, an equivalent stress has to be calculated in order to apply Neuber’s approximation concept. Therefore, the von Mises yield criterion [108] is applied, which is based on the second invariant J 2 of the deviatoric part of the stress tensor σ k . The derived equivalent stress represents the stress concentration factor K t , k ( u n i ) , which is calculated according to Equation (14). The superscript ( u n i ) denotes that the respective parameter is evaluated for the equivalent uniaxial stress state.
K t , k ( u n i ) = 3 · J 2 = 1 2 σ 11 , k σ 22 , k 2 + σ 22 , k σ 33 , k 2 + σ 33 , k σ 11 , k 2 + 6 σ 23 , k 2 + σ 31 , k 2 + σ 12 , k 2
Utilizing K t , k ( u n i ) and the cyclic stress–strain curve given by Equation (4), the elastic–plastic stresses and strains within C k can be approximated. At first, the uniaxial maximum stress σ ˜ m a x , k ( u n i ) during cyclic loading is calculated by minimizing the objective function given in Equation (15). The tilde accent denotes elastic–plastic parameters evaluated by the approximation framework.
1 E σ ˜ m a x , k ( u n i ) 2 + σ ˜ m a x , k ( u n i ) σ ˜ m a x , k ( u n i ) K 1 n 1 E Δ σ · K t , k ( u n i ) 1 R 2 = 0
To derive the maximum stress values of the components of the elastic–plastic stress tensor σ ˜ k from the approximated uniaxial maximum stress σ ˜ m a x , k ( u n i ) , an additional assumption has to be made. Hereby, the ratio between the linear–elastic components σ i j , k and the equivalent stress K t , k ( u n i ) is studied and it is assumed that this ratio remains valid in the case of non-linear material behaviour. Thus, the elastic–plastic maximum values σ ˜ i j , m a x , k of the components of σ ˜ k can be computed according to Equation (16).
σ ˜ i j , m a x , k = σ i j , k · σ ˜ m a x , k ( u n i ) K t , k ( u n i )
The values of σ ˜ i j , m a x , k constitute the upper reversal points of the approximated hysteresis loops during cyclic loading. In order to define also the lower reversal points, the elastic–plastic stress ranges Δ σ ˜ i j , k have to be calculated. As the investigated cast steel alloy G12MnMo7-4+QT exhibits Masing behaviour [109], the cyclic strains can be described as a function of the stress ranges according to Equation (17).
Δ ε ( t ) = Δ ε ( e ) + Δ ε ( p ) = Δ σ E + 2 Δ σ 2 · K 1 n
It should be noted that Equation (17) results from Equation (4) by substituting the amplitude terms for half of the respective ranges, for example, σ a = Δ σ / 2 . Thus, Equation (17) is also only valid for the uniaxial stress state. Utilizing K t , k ( u n i ) , the equivalent uniaxial stress range Δ σ ˜ k ( u n i ) of the finite element C k within Ω is calculated by minimizing the objective function given in Equation (18).
1 E Δ σ ˜ k ( u n i ) 2 + 2 · Δ σ ˜ k ( u n i ) Δ σ ˜ k ( u n i ) 2 · K 1 n 1 E Δ σ · K t , k ( u n i ) 2 = 0
In order to evaluate the elastic–plastic stress ranges Δ σ ˜ i j , k , the relationship reported in Equation (16) is utilised in terms of cyclic stresses according to Equation (19).
Δ σ ˜ i j , k = σ i j , k · Δ σ ˜ k ( u n i ) K t , k ( u n i )
Subsequently, the minimum stresses σ ˜ i j , m i n , k can be calculated following Equation (20), constituting the lower reversal points of the approximated cyclic hysteresis loops.
σ ˜ i j , m i n , k = σ ˜ i j , m a x , k Δ σ ˜ i j , k
As the cyclic hysteresis loops in element C k are characterized in terms of the acting stress components, the corresponding strain components of the elastic–plastic strain tensor ε ˜ k have to be evaluated. Utilizing Hooke’s law, the elastic strain ranges are computed according to Equation (21), where ν represents Poisson’s ratio and δ i j represents the Kronecker delta.
Δ ε ˜ i j , k ( e ) = 1 + ν E Δ σ ˜ i j , k ν E Δ σ ˜ 11 , k + Δ σ ˜ 22 , k + Δ σ ˜ 33 , k δ i j
While the elastic strain components in the multiaxial stress state can be analytically evaluated following Equation (21), the calculation of the plastic strain contributions Δ ε ˜ i j , k ( p ) needs some further consideration. The evaluation of the plastic strain ranges following Equation (17) is only valid for the uniaxial stress state, as the underlying experiments were conducted utilizing plain specimens subjected to uniaxial tension-compression loading. Hence, a direct calculation of the individual components Δ ε ˜ i j , k ( p ) of ε ˜ k utilizing Equation (17) is impermissible. Similar to Equations (16) and (19), an assumption based on the behaviour in the uniaxial stress state has to be made to approximate the plastic contributions Δ ε ˜ i j , k ( p ) . Therefore, the ratio between the plastic and elastic strain for the equivalent stress Δ σ ˜ k ( u n i ) is calculated based on Equation (17) and it is assumed that this ratio is valid also in the case of a multiaxial stress state. Hence, the plastic strain ranges Δ ε ˜ i j , k ( p ) are approximated according to Equation (22).
Δ ε ˜ i j , k ( p ) = Δ ε ˜ i j , k ( e ) · E · 2 n 1 K 1 n · Δ σ ˜ k ( u n i ) 1 n n
Finally, the total strain ranges Δ ε ˜ i j , k ( t ) are calculated as a superposition of the elastic and plastic contributions, as given in Equation (23).
Δ ε ˜ i j , k ( t ) = Δ ε ˜ i j , k ( e ) + Δ ε ˜ i j , k ( p )
As the stress and strain parameters describing the elastic–plastic hysteresis loops during cyclic loading are evaluated for the elements within the control volume Ω , the total averaged strain energy density Δ W ¯ ( t ) is calculated. The evaluation is carried out for each component of σ ˜ k . After computation of the individual contributions, a superposition is carried out to derive Δ W ¯ ( t ) . The general evaluation procedure is exemplified for element C k within Ω .
As illustrated in Figure 2, the cyclic stabilized hysteresis loop defines the elastic and plastic SED during cyclic loading. Hence, the respective contributions to the total SED can be evaluated simply by calculating the corresponding areas in the stress–strain diagram. In the uniaxial load case, the ascending and descending branches of the hysteresis loop are described by Equation (17). When dealing with a local multiaxial stress state, simply inserting the stress range component Δ σ ˜ i j , k into Equation (17) will therefore not result in the calculated total strain range Δ ε ˜ i j , k ( t ) . In order to derive a valid description of Δ ε ˜ i j , k ( t ) when dealing with individual stress components Δ σ ˜ i j , k , some other considerations have to be made. Stress and strain in a specific volume element are connected through the stiffness of the respective body. Hence, modulation of the stiffness leads to a change in strain when the acting stress parameter is kept constant. In terms of cyclic loading, the stiffness in the stress–strain relationship according to Equation (17) is described by E and K . The influence of the two parameters on the shape of the cyclic stabilized hysteresis loop is depicted in Figure 21.
Concerning Figure 21, the hysteresis loop depicted in blue constitutes the reference state. Varying E , as shown in Figure 21a, alters the elastic SED Δ W ( e ) , while the plastic SED Δ W ( p ) is unaffected. Then again, modulating K changes Δ W ( p ) , while Δ W ( e ) remains virtually unchanged, as depicted in Figure 21b.
The difference in strain response between the uniaxial and multiaxial stress state can be interpreted as some kind of constraint, which can be considered by a local change of stiffness in the volume element. In order to quantitatively describe this variation of stiffness, the modulation factors ξ i j , k ( e ) and ξ i j , k ( p ) are defined according to Equations (24) and (25). Both parameters are based on the comparison between the respective strain component Δ ε ˜ i j , k ( e ) or Δ ε ˜ i j , k ( p ) in the present multiaxial stress state and the corresponding hypothetical uniaxial strain component.
ξ i j , k ( e ) = Δ σ ˜ i j , k E · Δ ε ˜ i j , k ( e )
ξ i j , k ( p ) = Δ σ ˜ i j , k 2 1 n · K · Δ ε ˜ i j , k ( p ) n
Utilizing ξ i j , k ( e ) and ξ i j , k ( p ) , the stress–strain relationship given in Equation (26) is defined, which yields the total strain Δ ε ˜ i j , k ( t ) in the multiaxial stress state consistent with Equation (23).
Δ ε ˜ i j , k ( t ) = Δ ε ˜ i j , k ( e ) + Δ ε ˜ i j , k ( p ) = Δ σ ˜ i j , k ξ i j , k ( e ) · E + 2 Δ σ ˜ i j , k 2 · ξ i j , k ( p ) · K 1 n
As the stress–strain relation for the individual components of the multiaxial stress state is now defined according to Equation (26), the elastic and plastic SED contributions can be evaluated. At first, the calculation of the respective SED parts is carried out for a stress ratio of R = 0 , independent of the effective stress ratio R e f f , i j , k . Combining Equation (26) with trivial geometrical observations based on Figure 2, the total SED Δ W i j , k , R = 0 ( t ) is evaluated according to Equation (27).
Δ W i j , k , R = 0 ( t ) = Δ σ ˜ i j , k · Δ ε ˜ i j , k ( t ) 0 Δ σ ˜ i j , k σ ^ ξ i j , k ( e ) · E + 2 σ ^ 2 · ξ i j , k ( p ) · K 1 n d σ ^
Subsequently, the plastic SED Δ W i j , k , R = 0 ( p ) is calculated following Equation (28).
Δ W i j , k , R = 0 ( p ) = Δ σ ˜ i j , k · Δ ε ˜ i j , k ( t ) 2 0 Δ σ ˜ i j , k σ ^ ξ i j , k ( e ) · E + 2 σ ^ 2 · ξ i j , k ( p ) · K 1 n d σ ^
The elastic SED Δ W i j , k , R = 0 ( e ) represents the difference between the total and the plastic SED, as given by Equation (29).
Δ W i j , k , R = 0 ( e ) = Δ W i j , k , R = 0 ( t ) Δ W i j , k , R = 0 ( p ) = 0 Δ σ ˜ i j , k σ ^ ξ i j , k ( e ) · E + 2 σ ^ 2 · ξ i j , k ( p ) · K 1 n d σ ^
Finally, the SED components are calculated for the effective stress ratio R e f f , i j , k , which is induced by the mean stress relaxation during cyclic loading and is given by Equation (30).
R e f f , i j , k = σ ˜ i j , m i n , k σ ˜ i j , m a x , k
As the plastic SED relies solely on the shape of the cyclic stabilized hysteresis loop, it follows that Δ W i j , k ( p ) at R e f f , i j , k equals Δ W i j , k , R = 0 ( p ) given by Equation (28).
The elastic SED Δ W i j , k ( e ) depends on R e f f , i j , k and is evaluated according to Equation (31).
Δ W i j , k ( e ) = Δ W i j , k , R = 0 ( e ) + σ ˜ i j , m i n , k · Δ ε ˜ i j , k ( t ) 0 R e f f , i j , k 1 0 | σ ˜ i j , m i n , k | σ ^ ξ i j , k ( e ) · E + 2 σ ^ 2 · ξ i j , k ( p ) · K 1 n d σ ^ + 0 σ ˜ i j , m a x , k σ ^ ξ i j , k ( e ) · E + 2 σ ^ 2 · ξ i j , k ( p ) · K 1 n d σ ^ R e f f , i j , k < 0
As the elastic and plastic contributions are evaluated, the total SED Δ W i j , k ( t ) is calculated as their superposition following Equation (32).
Δ W i j , k ( t ) = Δ W i j , k ( e ) + Δ W i j , k ( p )
Figure 22 summarizes the individual steps of the approximation methodology in a graphical way.
The last step in the calculation procedure for the element C k constitutes the evaluation of the strain energy. Utilizing the element volume V E l e m e n t , k derived from the finite element model, the elastic, plastic and total strain energy contributions are evaluated following Equations (33)–(35).
Δ E i j , k ( e ) = Δ W i j , k ( e ) · V E l e m e n t , k
Δ E i j , k ( p ) = Δ W i j , k ( p ) · V E l e m e n t , k
Δ E i j , k ( t ) = Δ W i j , k ( t ) · V E l e m e n t , k
The approximation framework according to Equations (14)–(16) and (18)–(35) is utilised to calculate the strain energy terms for each of the n finite elements within Ω . In order to calculate the SED terms averaged over the whole control volume, the individual strain energy contributions of each element are summed and finally divided by the sum of the element volumes within Ω . To preserve consistency with the finite element simulations utilised to validate the approximation framework, it should be noted that each strain energy term Δ E i j , k ( e ) or Δ E i j , k ( p ) contributes either positively or negatively to the total sum. The respective sign depends on the ones of the stress and strain components and is implemented based on their product. As the resulting plastic strain energy is dissipated, the resulting sum of the individual contributions is defined as positive by convention. Finally, Equations (36)–(38) yield the averaged SED terms considering all n elements within Ω , where s g n represents the signum function.
Δ W ¯ ( e ) = k = 1 n i = 1 3 j = 1 3 s g n Δ σ ˜ i j , k · Δ ε ˜ i j , k ( e ) · Δ E i j , k ( e ) k = 1 n V E l e m e n t , k
Δ W ¯ ( p ) = k = 1 n i = 1 3 j = 1 3 s g n Δ σ ˜ i j , k · Δ ε ˜ i j , k ( p ) · Δ E i j , k ( p ) k = 1 n V E l e m e n t , k
Δ W ¯ ( t ) = Δ W ¯ ( e ) + Δ W ¯ ( p )
Equation (38) concludes the definition of the introduced elastic–plastic approximation framework. In order to prove its capabilities for the elastic–plastic fatigue assessment of notched and defect-afflicted cast steel components, the following section validates the methodology.

3.5. Validation of the Elastic–Plastic Approximation Framework

The obtained results of the approximation framework defined in Section 3.4 are validated utilizing non-linear finite element analyses. The validation procedure covers the non-linear simulations of the V-notched fatigue samples, which served as an experimental basis for the evaluation of the elastic–plastic design limit curve presented in Section 3.3. Both specimen geometries possessing notch opening angles of 2 α = 45 and 2 α = 135 are assessed by the presented elastic–plastic approximation framework, covering all experimentally tested stress ranges Δ σ and stress ratios R. The elastic and plastic SED terms Δ W ¯ ( e ) and Δ W ¯ ( p ) , as well as the total SED Δ W ¯ ( t ) , are evaluated and compared to the ones obtained by the finite element analyses. Figure 23 depicts a comparison of the non-linear simulation results and those obtained by the elastic–plastic approximation framework for both specimen geometries and a stress ratio of R = 1 .
Reflecting the results depicted in Figure 23, a sound accordance of the results obtained by elastic–plastic finite element analyses Δ W ¯ F E M and the ones derived by the presented approximation framework Δ W ¯ can be perceived. Both the elastic and plastic SED components show only a small deviation from the simulation results, therefore leading to a quite satisfying approximation of the total SED. Comparing (a) and (b) in Figure 23, the effect of the notch geometry can be perceived clearly. Considering loading of both notch geometries under an identical stress range Δ σ , the elastic and plastic SED contributions possess higher values at a notch opening angle of 2 α = 45 than at 2 α = 135 . To further investigate the capabilities of the introduced assessment method, the stabilized hysteresis loops evaluated by non-linear simulation and the presented approximation procedure are compared to each other. Figure 24 depicts respective hysteresis loops for the stress component σ 22 parallel to the specimen axis which yields the main contribution to Δ W ¯ ( t ) .
Table 7 summarizes the comparison between elastic–plastic simulations and the approximation framework for the investigated notch geometries and stress ratios in terms of the unsigned and signed deviations Δ u n s i g n e d and Δ s i g n e d , as well as the root mean square error Δ R M S ; see Equations (39)–(41).
Δ u n s i g n e d = 1 n · i = 1 n Δ W ¯ i ( t ) Δ W ¯ F E M , i ( t ) Δ W ¯ F E M , i ( t ) · 100 %
Δ s i g n e d = 1 n · i = 1 n Δ W ¯ i ( t ) Δ W ¯ F E M , i ( t ) Δ W ¯ F E M , i ( t ) · 100 %
Δ R M S = 1 n · i = 1 n Δ W ¯ i ( t ) Δ W ¯ F E M , i ( t ) Δ W ¯ F E M , i ( t ) · 100 % 2
The depicted values in Table 7 confirm the sound applicability of the presented approximation methodology as a maximum root mean square deviation of only 12.0 % is achieved. To put this value into perspective, it should be noted that this energy-based deviation translates to 5.8 % in terms of stress when considering a hypothetical linear–elastic material behaviour. In order to demonstrate the difference in numerical efficiency, Figure 25 depicts a comparison of the computation time needed for total strain energy density calculation using non-linear finite element analyses and the elastic–plastic approximation framework. The depicted values summarize the evaluation effort of all experimentally tested stress ranges Δ σ .
Considering Figure 25, the enhanced numerical efficiency obtained by the presented elastic–plastic approximation framework can be perceived clearly. Regarding the non-linear finite element analyses, the evaluation of the cyclic stabilized material behaviour constitutes the main contribution to the observed computation times. Due to the occurring elastic–plastic deformations, the investigated loads have to be applied in small increments in order to achieve an equilibrium condition in the implicit simulations. As the presented approximation procedure is solely based upon a quasi-static linear–elastic finite element analysis, no incremental load stepping has to be applied, which significantly reduces numerical efforts and therefore computation time. Hence, the application of the established approximation method enables a calculation of the total strain energy density which is on average approximately 400 times faster compared to non-linear finite element analyses.
Finally, Figure 26 depicts the approximation results for the notched specimens with 2 α = 45 and 2 α = 135 tested at R = 1 , R = 0 and R = 0.5 , which merge well into the previously evaluated scatter band; compare Figure 18 as a reference.

3.6. Calculation of the Elastic–Plastic Strain Energy Density of Bulk and Surface Defects

As the established approximation framework is validated, the methodology is applied to assess the fatigue strength of large-scale specimens affected by bulk and surface imperfections.
Starting with bulk defect-afflicted cast steel specimens, the consistency of the presented elastic–plastic approximation method in the linear–elastic fatigue regime is investigated at first. Therefore, the numerical results of the large-scale specimens given in [10] are reassessed utilizing Abaqus® 2021 and plane strain elements of the type CPE6/CPE8. Twenty defect contours are considered, including ten specimens tested at R = 1 and R = 0 , respectively. Furthermore, the selected specimens tested at R = 0.5 which are reported in Section 2.2 are also included. A stress range of Δ σ = 1 MPa acting in the cross-section with a diameter of 30 mm is utilised in order to evaluate the elastic SED Δ W ¯ D E F , 97.5 % ( e ) according to [10]. The evaluation of Δ W ¯ D E F , 97.5 % ( e ) includes the partitioning of individual control volumes along the respective defect contour which corresponds to Method 1 reported in a preceding study [104]. Based on the investigated control volumes, the fatigue parameter Δ W ¯ D E F , 97.5 % ( e ) is evaluated, which is defined as the averaged SED value that exceeds 97.5 % of all derived values along the defect perimeter [10]. Subsequent to the evaluation of Δ W ¯ D E F , 97.5 % ( e ) , the derived linear–elastic stress fields in the investigated control volumes along the defect perimeter are assessed by the methodology introduced in Section 3.4 and evaluated as a 97.5 % elastic–plastic SED exceedance limit Δ W ¯ D E F , 97.5 % ( t ) . Figure 27 depicts a comparison of Δ W ¯ D E F , 97.5 % ( e ) and Δ W ¯ D E F , 97.5 % ( t ) for a gross stress of Δ σ = 1 MPa.
The consistency between the approximation framework and the fatigue assessment method introduced in [10] in the linear–elastic fatigue regime can be perceived clearly. Nearly all evaluated data points are on the diagonal line, indicating a match of the calculated SED values. The derived statistical deviations possess values of Δ u n s i g n e d = 0.56 % , Δ s i g n e d = 0.19 % and Δ R M S = 1.24 % , confirming the sound applicability of the elastic–plastic approximation framework in the linear–elastic fatigue regime.
Subsequent to the consistency check, the elastic–plastic approximation framework is utilised for the elastic–plastic fatigue assessment of the bulk defect-afflicted large-scale specimens as reported in Section 2.2 and [10]. This includes the evaluation of the linear–elastic stress fields in the control volumes along the defect contour. Two-dimensional finite element models subjected to a gross stress of Δ σ = 1 MPa are used, featuring plane strain elements of type CPE6/CPE8. Considering the stress ranges investigated in the respective experiments, the total SED Δ W ¯ D E F , 97.5 % ( t ) is evaluated at increased stress levels of hundreds of Megapascals of gross stress under uniaxial tension, utilizing an exceedance limit of 97.5 % along the defect contour in accordance with [10]. It should be noted that each defect is assessed in terms of the two radiographs captured prior to testing. Hence, the energy-based fatigue assessment method yields two values of Δ W ¯ D E F , 97.5 % ( t ) for each individual defect.
Concerning the surface defect-afflicted specimens presented in Section 2.3, the energy-based fatigue assessment considers the three-dimensional image stacks of the individual defect structures, which were captured prior to testing. The image stack contains the spatial coordinates of the defect contour which are processed for projected edge detection. The scripted procedure includes the alignment of the defect structure on the specimen. The three-dimensional point cloud is utilised to derive a virtual radiograph by projecting the volumetric surface defect onto a plane which is parallel to the specimen axis, as shown in Figure 28.
Based on this planar edge projection, a linear–elastic plane strain finite element model is set up. Analogous to the assessment of the bulk defect-afflicted specimens, the finite element submodeling technique is utilised to derive the linear–elastic stress tensors of each element within the control volume Ω i . Subsequently, Ω i is assessed following the elastic–plastic approximation framework presented in Section 3.4. As the total SED is evaluated for Ω i , the numerical procedure is repeated until the whole defect contour has been assessed. As in the case of the bulk defect-afflicted specimens, the fatigue effective total SED Δ W ¯ D E F , 97.5 % ( t ) is derived by considering an exceedance limit of 97.5 % . Figure 29 exemplarily depicts the implemented finite element submodeling technique.
As the investigated bulk and surface defects are assessed by the elastic–plastic approximation framework, the evaluated data points Δ W ¯ D E F , 97.5 % ( t ) are plotted as a function of the experimentally derived numbers of load cycles. Figure 30 depicts the original scatter band, as evaluated in Section 3.3, as well as the fatigue assessment results of the defect-afflicted large-scale specimens.
Concerning Figure 30, the elastic–plastic fatigue assessment of bulk and surface defect-afflicted cast steel specimens yields satisfying results, as nearly all data points fall within the scatter band derived by small-scale specimen experiments. This highlights the ability of the total strain energy density concept to assess arbitrarily notched geometrical features. Bearing in mind that each data point in Figure 30 represents only a planar projection of a complexly shaped spatial imperfection, the depicted results also prove the presented approximation framework as an engineering-feasible fatigue assessment methodology. Focusing on the data points at a stress ratio of R = 0.5 , the implementation of the mean stress relaxation and the resulting effective stress ratio is emphasized, which constitutes a significant improvement in prediction accuracy compared to the linear–elastic SED concept. As associated numerical efforts are substantially reduced compared to non-linear finite element analyses, the elastic–plastic approximation method represents a promising tool to assess arbitrarily shaped defect geometries. In order to minimize computational efforts for complex defect geometries or imperfection networks even further, the implementation of a free-mesh approach as reported in [104] is encouraged for future work. Moreover, the presented elastic–plastic approximation framework can be expanded from planar case studies towards spatial fatigue assessment as well.

4. Conclusions

The present paper focuses on the elastic–plastic energy-based fatigue assessment of notched and imperfective cast steel components. Based on the conducted experimental and numerical investigations, the following scientific conclusions are drawn:
  • A design limit curve for the high-strength cast steel alloy G12MnMo7-4+QT was evaluated based on the elastic–plastic strain energy density concept and non-linear finite element analyses. The evaluated design limit curve considers the results of high-cycle fatigue experiments of three different specimen geometries tested at stress ratios of R = 1 , R = 0 and R = 0.5 . Implementation of the elastic–plastic material behaviour results in a significant improvement of prediction accuracy compared to the linear–elastic assessment conducted in a previous study, as the considered experimental results merge into a unique narrow scatter band which is characterized by a scatter index of 1 : T Δ W ¯ ( t ) = 2.43 .
  • The effect of mean stress relaxation results in the formation of an effective stress ratio that is smaller than the stress ratio of the applied load and generally tends towards R = 1 . For cyclic loading at stress ratios of R 0 , this leads to a decrease of the total SED compared to the results of the linear–elastic SED concept.
  • A novel elastic–plastic approximation approach to evaluate the total strain energy density was presented. The framework facilitates the rapid assessment of the elastic–plastic stress and strain tensors in the cyclic stabilized material state, invoking only linear–elastic finite element analysis and low-cycle fatigue properties. Hence, numerical efficiency is significantly increased compared to time-consuming non-linear finite element analyses.
  • Validation of the elastic–plastic approximation framework utilizing comprehensive elastic–plastic finite element analyses of two V-notched small-scale specimen geometries proved a sound agreement of obtained approximation and non-linear simulation results, showing an average root mean square deviation of only approximately eight percent.
  • The fatigue assessment of bulk and surface defect-afflicted large-scale specimens utilizing the introduced elastic–plastic approximation method results in quite satisfying fatigue life calculations. A total of 56 defect-afflicted large-scale specimens tested at stress ratios of R = 1 , R = 0 and R = 0.5 were investigated by the presented framework, yielding calculated fatigue life results which merge into the scatter band of the derived high-strength cast steel design fatigue curve.
  • As the elastic–plastic approximation framework is solely based on linear–elastic simulations currently utilizing planar projections of complexly shaped spatial imperfections, cost and time expenses are significantly decreased compared to three-dimensional non-linear finite element simulations. Hence, the presented approximation methodology presents an engineering-feasible design tool to assess the fatigue strength of bulk and surface defect-afflicted cast steel components.

Author Contributions

Conceptualization, M.H. and M.S.; methodology, M.H. and M.O.; software, M.H. and M.O.; validation, M.H., M.S. and M.O.; formal analysis, M.H. and M.O.; investigation, M.H. and M.O.; resources, M.S.; data curation, M.H. and M.O.; writing—original draft preparation, M.H.; writing—review and editing, M.H., M.O. and M.S.; visualization, M.H.; supervision, M.S.; project administration, M.S.; funding acquisition, M.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Austrian Federal Ministry for Digital and Economic Affairs and the National Foundation for Research, Technology, and Development.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

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 funders had no roles 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.

Abbreviations

The following abbreviations are used in this manuscript:
1 : T Δ W ¯ ( e ) Scatter index in terms of elastic strain energy density
1 : T Δ W ¯ ( t ) Scatter index in terms of total strain energy density
2 α Notch opening angle
γ i Kinematic material parameter
δ i j Kronecker delta
ε ˜ Elastic–plastic strain tensor
ε ( e ) Strain considering linear–elastic material behaviour
ε ( p ) Plastic strain
ε ( t ) Strain considering elastic–plastic material behaviour
ε a Strain amplitude
ε a ( e ) Elastic strain amplitude
ε a ( p ) Plastic strain amplitude
ε a ( t ) Total strain amplitude
ε f Fatigue ductility coefficient
ε i j Component of the linear–elastic strain tensor
ε Strain tensor
ν Poisson’s ratio
ξ i j ( e ) Elastic modulation factor assigned to component of the stress tensor
ξ i j ( p ) Plastic modulation factor assigned to component of the stress tensor
ρ Notch root radius
σ Stress tensor
σ Deviatoric part of the stress tensor
σ ( e ) Stress considering linear–elastic material behaviour
σ ( t ) Stress considering elastic–plastic material behaviour
σ a Stress amplitude
σ f Fatigue strength coefficient
σ i j Component of the linear–elastic stress tensor
σ m a x Maximum stress
σ m a x ( e ) Elastic maximum stress
σ m i n Minimum stress
σ Y S Yield strength
σ ^ Integration variable
σ ˜ Elastic–plastic stress tensor
σ ˜ i j , m a x Maximum value of the component of the elastic–plastic stress tensor
σ ˜ i j , m i n Minimum value of the component of the elastic–plastic stress tensor
σ ˜ m a x ( u n i ) Maximum stress in the equivalent uniaxial stress state
ψ Parameter to describe the direction of the flow
Δ R M S Root mean square error
Δ s i g n e d Signed deviation
Δ u n s i g n e d Unsigned deviation
Δ ε ( e ) Elastic strain range
Δ ε ( p ) Plastic strain range
Δ ε ( t ) Total strain range
Δ ε ˜ i j ( e ) Elastic range of the component of the elastic–plastic strain tensor
Δ ε ˜ i j ( p ) Plastic range of the component of the elastic–plastic strain tensor
Δ ε ˜ i j ( t ) Total range of the component of the elastic–plastic strain tensor
Δ σ Stress range
Δ σ i Stress range to be assessed
Δ σ ˜ ( u n i ) Equivalent uniaxial stress range
Δ σ ˜ i j Range of the component of the elastic–plastic stress tensor
Δ E i j ( e ) Elastic strain energy component
Δ E i j ( p ) Plastic strain energy component
Δ E i j ( t ) Total strain energy component
Δ W ( e ) Elastic strain energy density range
Δ W ( p ) Plastic strain energy density range
Δ W ( t ) Total strain energy density range
Δ W ¯ ( e ) Average elastic strain energy density range
Δ W ¯ ( p ) Average plastic strain energy density range
Δ W ¯ ( t ) Average total strain energy density range
Δ W ¯ D E F , 97.5 % ( e ) Average elastic strain energy density range for an exceedance limit of 97.5 %
Δ W ¯ D E F , 97.5 % ( t ) Average total strain energy density range for an exceedance limit of 97.5 %
Δ W i j ( e ) Elastic strain energy density range component
Δ W i j ( p ) Plastic strain energy density range component
Δ W i j ( t ) Total strain energy density range component
Δ W i j , R = 0 ( e ) Elastic strain energy density range component for R = 0
Δ W i j , R = 0 ( p ) Plastic strain energy density range component for R = 0
Δ W i j , R = 0 ( t ) Total strain energy density range component for R = 0
Δ W ¯ F E M ( e ) Average elastic strain energy density range derived by non-linear simulations
Δ W ¯ F E M ( p ) Average plastic strain energy density range derived by non-linear simulations
Δ W ¯ F E M ( t ) Average total strain energy density range derived by non-linear simulations
Ω Control volume
Ω i Control volume at the node P i
bFatigue strength exponent
b s Stabilization rate
cFatigue ductility exponent
fYield function
kInitial size of the yield function
m R Load multiplier for numerical cyclic load
n Cyclic strain hardening exponent
pAccumulated plastic strain
r 0 Distance between V-notch apex and origin of coordinate system
C i Kinematic material parameter
C k Element within the control volume
EYoung’s Modulus
E Stabilized cyclic Young’s Modulus
IIsotropic part of the combined hardening model
J ( σ X ) von Mises criterion
K Cyclic strength coefficient
K t Stress concentration factor
K t , m a x Maximum stress concentration factor
K t ( u n i ) Stress concentration factor in the equivalent uniaxial stress state
NNumber of cycles
N f Number of cycles to failure
P i Node on the defect contour in the finite element model
P S Probability of survival
QDifference between initial and stabilized size of yield surface
RStress ratio
R ε Strain ratio
R c Control radius
R e f f Effective stress ratio
R e f f , i j Effective stress ratio assigned to component of the stress tensor
SNominal stress in the net cross-section of a notched component
V E l e m e n t Volume of finite element
X Kinematic back stress tensor
X Deviatoric part of the kinematic back stress tensor
FKMForschungskuratorium Maschinenbau
SEDStrain energy density

References

  1. Fiedler, K.; Rolfe, B.F.; de Souza, T. Integrated Shape and Topology Optimization—Applications in Automotive Design and Manufacturing. SAE Int. J. Mater. Manuf. 2017, 10, 385–394. [Google Scholar] [CrossRef]
  2. Du, W.; Wang, H.; Zhu, L.; Zhao, Y.; Wang, Y.; Hao, R.; Yang, M. Innovative Joint for Cable Dome Structure Based on Topology Optimization and Additive Manufacturing. Materials 2021, 14, 5158. [Google Scholar] [CrossRef] [PubMed]
  3. Wang, L.; Du, W.; He, P.; Yang, M. Topology Optimization and 3D Printing of Three-Branch Joints in Treelike Structures. J. Struct. Eng. 2020, 146, 04019167. [Google Scholar] [CrossRef]
  4. Carlson, K.D.; Beckermann, C. Prediction of Shrinkage Pore Volume Fraction Using a Dimensionless Niyama Criterion. Metall. Mater. Trans. A 2009, 40, 163–175. [Google Scholar] [CrossRef]
  5. Futáš, P.; Pribulová, A.; Fedorko, G.; Molnár, V.; Junáková, A.; Laskovský, V. Failure analysis of a railway brake disc with the use of casting process simulation. Eng. Fail. Anal. 2019, 95, 226–238. [Google Scholar] [CrossRef]
  6. Jolly, M.; Katgerman, L. Modelling of defects in aluminium cast products. Prog. Mater. Sci. 2022, 123, 100824. [Google Scholar] [CrossRef]
  7. Kwon, H.J.; Kwon, H.K. Computer aided engineering (CAE) simulation for the design optimization of gate system on high pressure die casting (HPDC) process. Robot. Comput.-Integr. Manuf. 2019, 55, 147–153. [Google Scholar] [CrossRef]
  8. Patil, P.; Nalawade, R.; Balachandran, G.; Balasubramanian, V. Analysis of solidification behaviour of low alloy steel ingot casting—Simulation and experimental validation. Ironmak. Steelmak. 2015, 42, 512–524. [Google Scholar] [CrossRef]
  9. Meimandi, S.; Vanderesse, N.; Thibault, D.; Bocher, P.; Viens, M. Macro-defects characterization in cast CA-6NM martensitic stainless steel. Mater. Charact. 2017, 124, 31–39. [Google Scholar] [CrossRef]
  10. Horvath, M.; Oberreiter, M.; Stoschka, M. Energy-Based Fatigue Assessment of Defect-Afflicted Cast Steel Components by Means of a Linear–elastic Approach. Appl. Sci. 2023, 13, 3768. [Google Scholar] [CrossRef]
  11. Horvath, M.; Stoschka, M.; Fladischer, S. Fatigue strength study based on geometric shape of bulk defects in cast steel. Int. J. Fatigue 2022, 163, 107082. [Google Scholar] [CrossRef]
  12. Hardin, R.A.; Beckermann, C. Prediction of the Fatigue Life of Cast Steel Containing Shrinkage Porosity. Metall. Mater. Trans. A 2009, 40, 581–597. [Google Scholar] [CrossRef]
  13. Hardin, R.A.; Beckermann, C. Effect of Porosity on Deformation, Damage, and Fracture of Cast Steel. Metall. Mater. Trans. A 2013, 44, 5316–5332. [Google Scholar] [CrossRef]
  14. Makino, T.; Shimokawa, Y.; Yamamoto, M. Fatigue Property and Design Criterion of Cast Steel for Railway Bogie Frames. Mater. Trans. 2019, 60, 950–958. [Google Scholar] [CrossRef]
  15. Nadot, Y. Fatigue from Defect: Influence of Size, Type, Position, Morphology and Loading. Int. J. Fatigue 2022, 154, 106531. [Google Scholar] [CrossRef]
  16. Nagel, S.; Rauber, C.; Veile, I.; Knödel, P.; Ummenhofer, T. Influence of internal imperfections on the fatigue resistance of cast steel—Testing methodology. MATEC Web Conf. 2019, 300, 09001. [Google Scholar] [CrossRef]
  17. Schuscha, M.; Horvath, M.; Leitner, M.; Stoschka, M. Notch Stress Intensity Factor (NSIF)-Based Fatigue Design to Assess Cast Steel Porosity and Related Artificially Generated Imperfections. Metals 2019, 9, 1097. [Google Scholar] [CrossRef]
  18. Rennert, R.; Kullig, E.; Vormwald, M. Analytical Strength Assessment of Components Made of Steel, Cast Iron and Aluminium Materials in Mechanical Engineering: FKM Guideline, 6th revised ed.; VDMA Verlag: Frankfurt am Main, Germany, 2012. [Google Scholar]
  19. ASTM E446; Standard Reference Radiographs for Steel Castings Up to 2 in. (50.8 mm) in Thickness. ASTM International: West Conshohocken, PA, USA, 2020. [CrossRef]
  20. Carlson, K.; Ou, S.; Hardin, R.; Beckermann, C. Analysis of ASTM X-ray shrinkage rating for steel castings. In Proceedings of the 54th SFSA Technical and Operating Conference, Chicago, IL, USA, 3 November 2000; Steel Founder’s Society of America: Crystal Lake, IL, USA, 2000. [Google Scholar]
  21. Morrow, J.; Sinclair, G.M. Cycle-Dependent Stress Relaxation. In Symposium on Basic Mechanisms of Fatigue; American Society for Testing and Materials, Ed.; ASTM International: West Conshohocken, PA, USA, 1959; pp. 83–99. [Google Scholar] [CrossRef]
  22. Hao, H.; Ye, D.; Chen, Y.; Feng, M.; Liu, J. A study on the mean stress relaxation behavior of 2124-T851 aluminum alloy during low-cycle fatigue at different strain ratios. Mater. Des. 2015, 67, 272–279. [Google Scholar] [CrossRef]
  23. Arcari, A.; de Vita, R.; Dowling, N.E. Mean stress relaxation during cyclic straining of high strength aluminum alloys. Int. J. Fatigue 2009, 31, 1742–1750. [Google Scholar] [CrossRef]
  24. Paul, S.K.; Stanford, N.; Taylor, A.; Hilditch, T. The effect of low cycle fatigue, ratcheting and mean stress relaxation on stress–strain response and microstructural development in a dual phase steel. Int. J. Fatigue 2015, 80, 341–348. [Google Scholar] [CrossRef]
  25. Landgraf, R.W.; Chernenkoff, R.A. Residual stress effects on fatigue of surface processed steels. In Analytical and Experimental Methods for Residual Stress Effects in Fatigue, ASTM STP 1004; ASTM International: West Conshohocken, PA, USA, 1988; pp. 1–12. [Google Scholar]
  26. Lee, Y. Effects of redistributing residual stress on the fatigue behavior of SS330 weldment. Int. J. Fatigue 1998, 20, 565–573. [Google Scholar] [CrossRef]
  27. Chaboche, J.L. A review of some plasticity and viscoplasticity constitutive theories. Int. J. Plast. 2008, 24, 1642–1693. [Google Scholar] [CrossRef]
  28. Halama, R.; Sedlak, J.; Sofer, M. Phenomenological Modelling of Cyclic Plasticity. In Numerical Modelling; Miidla, P., Ed.; InTech: London, UK, 2012. [Google Scholar] [CrossRef]
  29. Bertini, L.; Le Bone, L.; Santus, C.; Chiesi, F.; Tognarelli, L. High Load Ratio Fatigue Strength and Mean Stress Evolution of Quenched and Tempered 42CrMo4 Steel. J. Mater. Eng. Perform. 2017, 26, 3784–3793. [Google Scholar] [CrossRef]
  30. Chaboche, J.L. Time-independent constitutive theories for cyclic plasticity. Int. J. Plast. 1986, 2, 149–188. [Google Scholar] [CrossRef]
  31. Seisenbacher, B.; Winter, G.; Grün, F. Improved Approach to Determine the Material Parameters for a Combined Hardening Model. Mater. Sci. Appl. 2018, 9, 357–367. [Google Scholar] [CrossRef]
  32. Koo, S.; Han, J.; Marimuthu, K.P.; Lee, H. Determination of Chaboche combined hardening parameters with dual backstress for ratcheting evaluation of AISI 52100 bearing steel. Int. J. Fatigue 2019, 122, 152–163. [Google Scholar] [CrossRef]
  33. Santus, C.; Grossi, T.; Romanelli, L.; Pedranz, M.; Benedetti, M. A computationally fast and accurate procedure for the identification of the Chaboche isotropic-kinematic hardening model parameters based on strain-controlled cycles and asymptotic ratcheting rate. Int. J. Plast. 2023, 160, 103503. [Google Scholar] [CrossRef]
  34. Liu, S.; Liang, G.; Yang, Y. A strategy to fast determine Chaboche elasto-plastic model parameters by considering ratcheting. Int. J. Press. Vessel. Pip. 2019, 172, 251–260. [Google Scholar] [CrossRef]
  35. Li, X.; Sun, Q.; Zhang, X. Identification of combined hardening model parameters by considering pre-tensile stress for 2024-T3 aluminum alloy. Mater. Today Commun. 2021, 26, 102065. [Google Scholar] [CrossRef]
  36. Hemmesi, K.; Mallet, P.; Farajian, M. Numerical evaluation of surface welding residual stress behavior under multiaxial mechanical loading and experimental validations. Int. J. Mech. Sci. 2020, 168, 105127. [Google Scholar] [CrossRef]
  37. Zhuang, X.; Ma, Y.; Zhao, Z. Fracture prediction under nonproportional loadings by considering combined hardening and fatigue-rule-based damage accumulation. Int. J. Mech. Sci. 2019, 150, 51–65. [Google Scholar] [CrossRef]
  38. Pandey, V.B.; Singh, I.V.; Mishra, B.K. A Strain-based continuum damage model for low cycle fatigue under different strain ratios. Eng. Fract. Mech. 2021, 242, 107479. [Google Scholar] [CrossRef]
  39. Mao, J.; Yang, J.; Zhu, J.; Ding, Z.; Zhong, F.; Wang, D. Experimental and theoretical research on creep–fatigue behaviors of 316L steel with and without 650 °C thermal aging. Fatigue Fract. Eng. Mater. Struct. 2022, 45, 1179–1198. [Google Scholar] [CrossRef]
  40. Lazzarin, P.; Zambardi, R. A finite-volume-energy based approach to predict the static and fatigue behavior of components with sharp V-shaped notches. Int. J. Fract. 2001, 112, 275–298. [Google Scholar] [CrossRef]
  41. Yosibash, Z.; Bussiba, A.; Gilad, I. Failure criteria for brittle elastic materials. Int. J. Fract. 2004, 125, 307–333. [Google Scholar] [CrossRef]
  42. Lazzarin, P.; Berto, F. Some Expressions for the Strain Energy in a Finite Volume Surrounding the Root of Blunt V-notches. Int. J. Fract. 2005, 135, 161–185. [Google Scholar] [CrossRef]
  43. Ellyin, F.; Kujawski, D. Plastic Strain Energy in Fatigue Failure. J. Press. Vessel Technol. 1984, 106, 342. [Google Scholar] [CrossRef]
  44. Gómez, F.J.; Elices, M.; Berto, F.; Lazzarin, P. Local strain energy to assess the static failure of U-notches in plates under mixed mode loading. Int. J. Fract. 2007, 145, 29–45. [Google Scholar] [CrossRef]
  45. Lazzarin, P.; Berto, F.; Elices, M.; Gómez, J. Brittle failures from U- and V-notches in mode I and mixed, I + II, mode: A synthesis based on the strain energy density averaged on finite-size volumes. Fatigue Fract. Eng. Mater. Struct. 2009, 32, 671–684. [Google Scholar] [CrossRef]
  46. Lazzarin, P.; Livieri, P.; Berto, F.; Zappalorto, M. Local strain energy density and fatigue strength of welded joints under uniaxial and multiaxial loading. Eng. Fract. Mech. 2008, 75, 1875–1889. [Google Scholar] [CrossRef]
  47. Berto, F.; Campagnolo, A.; Meneghetti, G.; Tanaka, K. Averaged strain energy density-based synthesis of crack initiation life in notched steel bars under torsional fatigue. Frat. Integrita Strutt. 2016, 10, 215–223. [Google Scholar] [CrossRef]
  48. Berto, F.; Lazzarin, P. A review of the volume-based strain energy density approach applied to V-notches and welded structures. Theor. Appl. Fract. Mech. 2009, 52, 183–194. [Google Scholar] [CrossRef]
  49. Gallo, P.; Berto, F. Advanced Materials for Applications at High Temperature: Fatigue Assessment by Means of Local Strain Energy Density. Adv. Eng. Mater. 2016, 18, 2010–2017. [Google Scholar] [CrossRef]
  50. Livieri, P.; Lazzarin, P. Fatigue strength of steel and aluminium welded joints based on generalised stress intensity factors and local strain energy values. Int. J. Fract. 2005, 133, 247–276. [Google Scholar] [CrossRef]
  51. Razavi, S.; Ferro, P.; Berto, F.; Torgersen, J. Fatigue strength of blunt V-notched specimens produced by selective laser melting of Ti-6Al-4V. Theor. Appl. Fract. Mech. 2017, 97, 376–384. [Google Scholar] [CrossRef]
  52. Ding, K.; Yang, Y.; Wang, Z.; Zhang, T.; Guo, W. Relationship between local strain energy density and fatigue life of riveted Al-Li alloy plate. Theor. Appl. Fract. Mech. 2022, in press. [Google Scholar] [CrossRef]
  53. Foti, P.; Santonocito, D.; Risitano, G.; Berto, F. Fatigue assessment of cruciform joints: Comparison between Strain Energy Density predictions and current standards and recommendations. Eng. Struct. 2021, 230, 111708. [Google Scholar] [CrossRef]
  54. Foti, P.; Berto, F. Fatigue assessment of high strength welded joints through the strain energy density method. Fatigue Fract. Eng. Mater. Struct. 2020, 43, 2694–2702. [Google Scholar] [CrossRef]
  55. Schuscha, M.; Leitner, M.; Stoschka, M.; Meneghetti, G. Local strain energy density approach to assess the fatigue strength of sharp and blunt V-notches in cast steel. Int. J. Fatigue 2020, 132, 105334. [Google Scholar] [CrossRef]
  56. Martins, R.F.; Branco, R.; Long, X. Fatigue Life Assessment in Bainitic Steels Based on The Cumulative Strain Energy Density. Appl. Sci. 2020, 10, 7774. [Google Scholar] [CrossRef]
  57. Ellyin, F.; Golos, K. Multiaxial Fatigue Damage Criterion. J. Eng. Mater. Technol. 1988, 110, 63–68. [Google Scholar] [CrossRef]
  58. Golos, K.; Ellyin, F. Total Strain Energy Density as a Fatigue Damage Parameter. In Advances in Fatigue Science and Technology; Branco, C.M., Rosa, L.G., Eds.; Springer: Dordrecht, The Netherlands, 1989; pp. 849–858. [Google Scholar] [CrossRef]
  59. Dallmeier, J.; Denk, J.; Huber, O.; Saage, H.; Eigenfeld, K. Deformation Behavior and Fatigue Analysis of Magnesium Wrought Alloys under Variable Amplitude Loading. Mater. Today Proc. 2015, 2, S119–S124. [Google Scholar] [CrossRef]
  60. Benedetti, M.; Berto, F.; Le Bone, L.; Santus, C. A novel Strain-Energy-Density based fatigue criterion accounting for mean stress and plasticity effects on the medium-to-high-cycle uniaxial fatigue strength of plain and notched components. Int. J. Fatigue 2020, 133, 105397. [Google Scholar] [CrossRef]
  61. Lazzarin, P.; Sonsino, C.M.; Zambardi, R. A notch stress intensity approach to assess the multiaxial fatigue strength of welded tube-to-flange joints subjected to combined loadings. Fatigue Fract. Eng. Mater. Struct. 2004, 27, 127–140. [Google Scholar] [CrossRef]
  62. Berto, F.; Lazzarin, P.; Yates, J.R. Multiaxial fatigue of V-notched steel specimens: A non-conventional application of the local energy method. Fatigue Fract. Eng. Mater. Struct. 2011, 34, 921–943. [Google Scholar] [CrossRef]
  63. Berto, F.; Campagnolo, A.; Lazzarin, P. Fatigue strength of severely notched specimens made of Ti-6Al-4V under multiaxial loading. Fatigue Fract. Eng. Mater. Struct. 2015, 38, 503–517. [Google Scholar] [CrossRef]
  64. Li, X.K.; Chen, S.; Zhu, S.P.; Liao, D.; Gao, J.W. Probabilistic fatigue life prediction of notched components using strain energy density approach. Eng. Fail. Anal. 2021, 124, 105375. [Google Scholar] [CrossRef]
  65. Branco, R.; Prates, P.; Costa, J.D.; Cruces, A.; Lopez-Crespo, P.; Berto, F. On the applicability of the cumulative strain energy density for notch fatigue analysis under multiaxial loading. Theor. Appl. Fract. Mech. 2022, 120, 103405. [Google Scholar] [CrossRef]
  66. Branco, R.; Martins, R.F.; Correia, J.; Marciniak, Z.; Macek, W.; Jesus, J. On the use of the cumulative strain energy density for fatigue life assessment in advanced high-strength steels. Int. J. Fatigue 2022, 164, 107121. [Google Scholar] [CrossRef]
  67. Ren, Z.; Qin, X.; Zhang, Q.; Sun, Y. Multiaxial fatigue life prediction model based on an improved strain energy density criterion. J. Press. Vessel. Pip. 2022, 199, 104724. [Google Scholar] [CrossRef]
  68. Branco, R.; Prates, P.A.; Costa, J.D.; Borrego, L.P.; Berto, F.; Kotousov, A.; Antunes, F.V. Rapid assessment of multiaxial fatigue lifetime in notched components using an averaged strain energy density approach. Int. J. Fatigue 2019, 124, 89–98. [Google Scholar] [CrossRef]
  69. Neuber, H. Theory of Stress Concentration for Shear-Strained Prismatical Bodies With Arbitrary Nonlinear Stress-Strain Law. J. Appl. Mech. 1961, 28, 544–550. [Google Scholar] [CrossRef]
  70. Molski, K.; Glinka, G. A method of elastic–plastic stress and strain calculation at a notch root. Mater. Sci. Eng. 1981, 50, 93–100. [Google Scholar] [CrossRef]
  71. Seeger, T.; Heuler, P. Generalized Application of Neuber’s Rule. J. Test. Eval. 1980, 8, 6. [Google Scholar] [CrossRef]
  72. Liao, D.; Zhu, S.P.; Correia, J.A.; de Jesus, A.M.; Berto, F. Recent advances on notch effects in metal fatigue: A review. Fatigue Fract. Eng. Mater. Struct. 2020, 43, 637–659. [Google Scholar] [CrossRef]
  73. Visvanatha, S.K.; Straznicky, P.V.; Hewitt, R.L. Influence of strain estimation methods on life predictions using the local strain approach. Int. J. Fatigue 2000, 22, 675–681. [Google Scholar] [CrossRef]
  74. Gu, H.H.; Wang, R.Z.; Tang, M.J.; Zhang, X.C.; Tu, S.T. Creep-fatigue reliability assessment for high-temperature components fusing on-line monitoring data and physics-of-failure by engineering damage mechanics approach. Int. J. Fatigue 2023, 169, 107481. [Google Scholar] [CrossRef]
  75. Hatami, F.; Varvani-Farahani, A. Accumulation of Plastic Strain at Notch Root of Steel Specimens Undergoing Asymmetric Fatigue Cycles: Analysis and Simulation. Materials 2023, 16, 2153. [Google Scholar] [CrossRef] [PubMed]
  76. Jena, S.K.; Arora, P.; Gupta, S.K.; Chattopadhyay, J. Validation of notch stress estimation schemes for low C–Mn steel. Fatigue Fract. Eng. Mater. Struct. 2022, 45, 1145–1164. [Google Scholar] [CrossRef]
  77. Seiler, M.; Linse, T.; Hantschke, P.; Kästner, M. An efficient phase-field model for fatigue fracture in ductile materials. Eng. Fract. Mech. 2020, 224, 106807. [Google Scholar] [CrossRef]
  78. Yadegari, P.; Schlitzer, T.; Vormwald, M. Approximation of elastic–plastic local strains in surface-hardened notched components. MATEC Web Conf. 2021, 349, 04008. [Google Scholar] [CrossRef]
  79. Glinka, G. Energy density approach to calculation of inelastic strain-stress near notches and cracks. Eng. Fract. Mech. 1985, 22, 485–508. [Google Scholar] [CrossRef]
  80. Singh, M.N.K.; Glinka, G.; Dubey, R.N. Elastic–plastic stress–strain calculation in notched bodies subjected to non-proportional loading. Int. J. Fract. 1996, 76, 39–60. [Google Scholar] [CrossRef]
  81. Ye, D.; Matsuoka, S.; Suzuki, N.; Maeda, Y. Further investigation of Neuber’s rule and the equivalent strain energy density (ESED) method. Int. J. Fatigue 2004, 26, 447–455. [Google Scholar] [CrossRef]
  82. Ye, D.; Hertel, O.; Vormwald, M. A unified expression of elastic–plastic notch stress–strain calculation in bodies subjected to multiaxial cyclic loading. Int. J. Solids Struct. 2008, 45, 6177–6189. [Google Scholar] [CrossRef]
  83. Li, J.; Zhang, Z.p.; Li, C.w. A Coupled Armstrong-Frederick Type Plasticity Correction Methodology for Calculating Multiaxial Notch Stresses and Strains. J. Fail. Anal. Prev. 2017, 17, 706–716. [Google Scholar] [CrossRef]
  84. Tao, Z.Q.; Shang, D.G.; Sun, Y.J. New pseudo stress correction method for estimating local strains at notch under multiaxial cyclic loading. Int. J. Fatigue 2017, 103, 280–293. [Google Scholar] [CrossRef]
  85. Li, D.H.; Shang, D.G.; Xue, L.; Li, L.J.; Wang, L.W.; Cui, J. Notch stress–strain estimation method based on pseudo stress correction under multiaxial thermo-mechanical cyclic loading. Int. J. Solids Struct. 2020, 199, 144–157. [Google Scholar] [CrossRef]
  86. Lutovinov, M.; Halama, R.; Papuga, J.; Bartošák, M.; Kuželka, J.; Růžička, M. An Approximate Method for Calculating Elastic–plastic Stress and Strain on Notched Specimens. Materials 2022, 15, 1432. [Google Scholar] [CrossRef]
  87. Ince, A.; Glinka, G. Approximation modeling framework for elastic–plastic stress–strain fields near cracks with a small finite crack tip radius. Theor. Appl. Fract. Mech. 2022, 121, 103452. [Google Scholar] [CrossRef]
  88. SEW 520; Hochfester Stahlguß mit Guter Schweißeignung. Verlag Stahleisen GmbH: Düsseldorf, Germany, 1996.
  89. DIN EN 3988; Test Methods for Metallic Materials—Constant Amplitude Strain-Controlled Low Cycle Fatigue Testing. German Institute for Standardisation: Berlin, Germany, 1998.
  90. Basquin, O. The exponential law of endurance tests. Am. Soc. Test. Mater. Proc. 1910, 10, 625–630. [Google Scholar]
  91. Manson, S.S. Behavior of Materials under Conditions of Thermal Stress; National Advisory Committee for Aeronautics: Cleveland, OH, USA, 1954.
  92. Coffin, L.F. A Study of the Effects of Cyclic Thermal Stresses on a Ductile Metal. Trans. ASME 1954, 76, 931–950. [Google Scholar] [CrossRef]
  93. Manson, S.S. Fatigue: A complex subject—Some simple approximations. Exp. Mech. 1965, 5, 193–226. [Google Scholar] [CrossRef]
  94. Ramberg, W.; Osgood, W.R. Description of Stress-Strain Curves by Three Parameters; NASA: Washington, DC, USA, 1943. [Google Scholar]
  95. Lukáš, P.; Klesnil, M. Cyclic stress–strain response and fatigue life of metals in low amplitude region. Mater. Sci. Eng. 1973, 11, 345–356. [Google Scholar] [CrossRef]
  96. Nieslony, A.; Dsoki, C.; Kaufmann, H.; Krug, P. New method for evaluation of the Manson–Coffin–Basquin and Ramberg–Osgood equations with respect to compatibility. Int. J. Fatigue 2008, 30, 1967–1977. [Google Scholar] [CrossRef]
  97. Chaboche, J.L.; Van, K.D.; Cordier, G. Modelization of the Strain Memory Effect on the Cyclic Hardening of 316 Stainless Steel; NC State University: Raleigh, NC, USA, 1979. [Google Scholar]
  98. Armstrong, P.J.; Frederick, C.O. A Mathematical Representation of the Multiaxial Bauschinger Effect; G.E.G.B. Report RD/B/N 731; Berkeley Nuclear Laboratories: Berkeley, CA, USA, 1966. [Google Scholar]
  99. Himmelbauer, F.; Seisenbacher, B.T.; Winter, G.; Grün, F. (Eds.) Fatigue Life Assessment and Modelling of the Cyclic Material Behaviour of a High Strength Martensitic Steel; Deutscher Verband für Materialforschung und -Prüfung e.V.: Berlin, Germany, 2022. [Google Scholar] [CrossRef]
  100. Himmelbauer, F.; Winter, G.; Grün, F.; Kiesling, C. On the Transferability of Fatigue and Cyclic Deformation Data to 100 µm Thin Structures. Metals 2022, 12, 1524. [Google Scholar] [CrossRef]
  101. Himmelbauer, F.; Winter, G.; Seisenbacher, B.; Grün, F.; Kiesling, C. Very high cycle fatigue assessment at elevated temperature of 100 µm thin structures made of high-strength steel X5CrNiCuNb16-4. J. Mater. Res. Technol. 2022, 21, 1811–1829. [Google Scholar] [CrossRef]
  102. Seisenbacher, B.; Winter, G.; Grün, F. Investigation and simulation of the age-dependent behaviour of an aluminium alloy under thermomechanical loading. Mech. Mater. 2020, 151, 103612. [Google Scholar] [CrossRef]
  103. Seisenbacher, B.; Winter, G.; Grün, F. Modelling the effect of ageing on the yield strength of an aluminium alloy under cyclic loading at different ageing temperatures and test temperatures. Int. J. Fatigue 2020, 137, 105635. [Google Scholar] [CrossRef]
  104. Stoschka, M.; Horvath, M.; Fladischer, S.; Oberreiter, M. Study of Local Fatigue Methods (TCD, N-SIF, and ESED) on Notches and Defects Related to Numerical Efficiency. Appl. Sci. 2023, 13, 2247. [Google Scholar] [CrossRef]
  105. Lazzarin, P.; Berto, F. Control volumes and strain energy density under small and large scale yielding due to tension and torsion loading. Fatigue Fract. Eng. Mater. Struct. 2008, 31, 95–107. [Google Scholar] [CrossRef]
  106. Dengel, D. Die arc sin √P-Transformation—Ein einfaches Verfahren zur grafischen und rechnerischen Auswertung geplanter Wöhlerversuche. Mater. Werkst. 1975, 6, 253–261. [Google Scholar] [CrossRef]
  107. ASTM E739; Standard Practice for Statistical Analysis of Linear or Linearized Stress-Life (S-N) and Strain-Life (e-N) Fatigue Data. ASTM International: West Conshohocken, PA, USA. [CrossRef]
  108. Von Mises, R. Mechanik der festen Körper im plastisch-deformablen Zustand. Nachrichten Ges. Wiss. Göttingen-Math-Phys. Kl. 1913, 1913, 582–592. [Google Scholar]
  109. Masing, G. Eigenspannungen und Verfestigung beim Messing. In Proceedings of the 2nd International Congress of Applied Mechanics, Zürich, Switzerland, 12–17 September 1926; pp. 332–335. [Google Scholar]
Figure 1. Definition of the control volume Ω : (a) Sharp V-Notch. (b) Blunt V-Notch.
Figure 1. Definition of the control volume Ω : (a) Sharp V-Notch. (b) Blunt V-Notch.
Applmech 04 00030 g001
Figure 2. Definition of the elastic and plastic strain energy density on a stabilized hysteresis loop.
Figure 2. Definition of the elastic and plastic strain energy density on a stabilized hysteresis loop.
Applmech 04 00030 g002
Figure 3. Graphical interpretation of Neuber’s concept.
Figure 3. Graphical interpretation of Neuber’s concept.
Applmech 04 00030 g003
Figure 4. Plain and notched specimen geometries investigated in a preceding study [10].
Figure 4. Plain and notched specimen geometries investigated in a preceding study [10].
Applmech 04 00030 g004
Figure 5. Experimental setup of small-scale specimen high-cycle fatigue investigations conducted in a preceding study [10].
Figure 5. Experimental setup of small-scale specimen high-cycle fatigue investigations conducted in a preceding study [10].
Applmech 04 00030 g005
Figure 6. Specimen geometry utilised in low-cycle fatigue experiments.
Figure 6. Specimen geometry utilised in low-cycle fatigue experiments.
Applmech 04 00030 g006
Figure 7. Experimental setup of low-cycle fatigue investigations.
Figure 7. Experimental setup of low-cycle fatigue investigations.
Applmech 04 00030 g007
Figure 8. Geometry of bulk defect-afflicted large-scale specimens utilised in high-cycle fatigue experiments.
Figure 8. Geometry of bulk defect-afflicted large-scale specimens utilised in high-cycle fatigue experiments.
Applmech 04 00030 g008
Figure 9. Experimental setup utilised for high-cycle fatigue investigations of defect-afflicted large-scale specimens.
Figure 9. Experimental setup utilised for high-cycle fatigue investigations of defect-afflicted large-scale specimens.
Applmech 04 00030 g009
Figure 10. Radiographs and fracture surfaces of representative bulk defect-afflicted large-scale specimens tested at R = 0.5 : (a) Specimen #33. (b) Specimen #35. (c) Specimen #36.
Figure 10. Radiographs and fracture surfaces of representative bulk defect-afflicted large-scale specimens tested at R = 0.5 : (a) Specimen #33. (b) Specimen #35. (c) Specimen #36.
Applmech 04 00030 g010
Figure 11. Geometry of surface defect-afflicted large-scale specimens utilised in high-cycle fatigue experiments.
Figure 11. Geometry of surface defect-afflicted large-scale specimens utilised in high-cycle fatigue experiments.
Applmech 04 00030 g011
Figure 12. Resulting stress distribution around the clamping area of the optimized specimen design.
Figure 12. Resulting stress distribution around the clamping area of the optimized specimen design.
Applmech 04 00030 g012
Figure 13. Surface defects and fracture surfaces of investigated large-scale specimens: (a) Specimen #02 tested at R = 1 . (b) Specimen #09 tested at R = 0 . (c) Specimen #17 tested at R = 0.5 .
Figure 13. Surface defects and fracture surfaces of investigated large-scale specimens: (a) Specimen #02 tested at R = 1 . (b) Specimen #09 tested at R = 0 . (c) Specimen #17 tested at R = 0.5 .
Applmech 04 00030 g013
Figure 14. Manson–Coffin curve for R ε = 1 .
Figure 14. Manson–Coffin curve for R ε = 1 .
Applmech 04 00030 g014
Figure 15. Cyclic stress–strain curve for the stabilized material behaviour.
Figure 15. Cyclic stress–strain curve for the stabilized material behaviour.
Applmech 04 00030 g015
Figure 16. Validation of the evaluated combined hardening model: (a) Single element test model and loading sequence. (b) Comparison of numerically and experimentally derived stabilized hysteresis loop for ε a ( t ) = 0.5 % .
Figure 16. Validation of the evaluated combined hardening model: (a) Single element test model and loading sequence. (b) Comparison of numerically and experimentally derived stabilized hysteresis loop for ε a ( t ) = 0.5 % .
Applmech 04 00030 g016
Figure 17. Finite element analysis strategy for the assessment of investigated small-scale specimens: (a) Numerical model and boundary conditions. (b) Course of the load multiplier m R utilised to apply numerical cyclic loading.
Figure 17. Finite element analysis strategy for the assessment of investigated small-scale specimens: (a) Numerical model and boundary conditions. (b) Course of the load multiplier m R utilised to apply numerical cyclic loading.
Applmech 04 00030 g017
Figure 18. Elastic–plastic strain energy density-based fatigue assessment of small-scale specimens.
Figure 18. Elastic–plastic strain energy density-based fatigue assessment of small-scale specimens.
Applmech 04 00030 g018
Figure 19. Comparison of linear–elastic and elastic–plastic strain energy density for a sharp V-notched specimen geometry.
Figure 19. Comparison of linear–elastic and elastic–plastic strain energy density for a sharp V-notched specimen geometry.
Applmech 04 00030 g019
Figure 20. Evaluation of the linear–elastic stress tensor for each element within the control volume.
Figure 20. Evaluation of the linear–elastic stress tensor for each element within the control volume.
Applmech 04 00030 g020
Figure 21. Influence of the parameters E and K on the shape of the cyclic stabilized hysteresis loop: (a) Influence of E , Δ σ and K are kept constant. (b) Influence of K , Δ σ and E are kept constant.
Figure 21. Influence of the parameters E and K on the shape of the cyclic stabilized hysteresis loop: (a) Influence of E , Δ σ and K are kept constant. (b) Influence of K , Δ σ and E are kept constant.
Applmech 04 00030 g021
Figure 22. Evaluation of Δ W i j , k ( t ) for the finite element C k by the introduced approximation procedure.
Figure 22. Evaluation of Δ W i j , k ( t ) for the finite element C k by the introduced approximation procedure.
Applmech 04 00030 g022
Figure 23. Comparison of assessment results obtained by non-linear finite element analyses and the presented approximation methodology for R = 1 : (a) V-Notch with 2 α = 45 . (b) V-Notch with 2 α = 135 .
Figure 23. Comparison of assessment results obtained by non-linear finite element analyses and the presented approximation methodology for R = 1 : (a) V-Notch with 2 α = 45 . (b) V-Notch with 2 α = 135 .
Applmech 04 00030 g023
Figure 24. Comparison of stabilized hysteresis loops evaluated by elastic–plastic simulations and the presented approximation methodology.
Figure 24. Comparison of stabilized hysteresis loops evaluated by elastic–plastic simulations and the presented approximation methodology.
Applmech 04 00030 g024
Figure 25. Comparison of computation time to evaluate the total strain energy density based on non-linear finite element analyses and the presented elastic–plastic approximation framework: (a) V-Notch with 2 α = 45 . (b) V-Notch with 2 α = 135 .
Figure 25. Comparison of computation time to evaluate the total strain energy density based on non-linear finite element analyses and the presented elastic–plastic approximation framework: (a) V-Notch with 2 α = 45 . (b) V-Notch with 2 α = 135 .
Applmech 04 00030 g025
Figure 26. Fatigue assessment results of V-notched small-scale specimens based on the presented approximation framework.
Figure 26. Fatigue assessment results of V-notched small-scale specimens based on the presented approximation framework.
Applmech 04 00030 g026
Figure 27. Consistency check of the established approximation method in the linear–elastic fatigue regime.
Figure 27. Consistency check of the established approximation method in the linear–elastic fatigue regime.
Applmech 04 00030 g027
Figure 28. Schematic representation of the evaluation of the planar defect projection.
Figure 28. Schematic representation of the evaluation of the planar defect projection.
Applmech 04 00030 g028
Figure 29. Finite element analysis strategy to assess surface defect-afflicted large-scale specimens.
Figure 29. Finite element analysis strategy to assess surface defect-afflicted large-scale specimens.
Applmech 04 00030 g029
Figure 30. Elastic–plastic fatigue assessment of defect-afflicted large-scale specimens based on the presented approximation framework.
Figure 30. Elastic–plastic fatigue assessment of defect-afflicted large-scale specimens based on the presented approximation framework.
Applmech 04 00030 g030
Table 1. Shape parameters of investigated bulk defects evaluated based upon radiographs.
Table 1. Shape parameters of investigated bulk defects evaluated based upon radiographs.
ParameterMinimumMaximumMeanStandard Deviation
Area [ mm 2 ]89.99160.98132.5422.10
Feret Max. [mm]25.5549.5336.518.20
Feret Min. [mm]6.629.968.070.87
Aspect ratio [−]3.057.474.601.32
Circularity [−]0.080.320.180.08
Convexity [−]0.470.720.610.08
Elongation [−]6.9718.0410.283.74
Table 2. Shape parameters of investigated bulk defects evaluated on the fracture surfaces of respective specimens.
Table 2. Shape parameters of investigated bulk defects evaluated on the fracture surfaces of respective specimens.
ParameterMinimumMaximumMeanStandard Deviation
Area [ mm 2 ]22.4163.4435.7215.93
Feret Max. [mm]7.5313.029.472.14
Feret Min. [mm]4.588.506.411.26
Aspect ratio [−]1.071.971.480.36
Circularity [−]0.110.300.200.07
Convexity [−]0.640.820.760.07
Elongation [−]1.893.772.660.63
Table 3. Shape parameters of investigated surface defects evaluated on the shell surfaces of respective specimens.
Table 3. Shape parameters of investigated surface defects evaluated on the shell surfaces of respective specimens.
ParameterMinimumMaximumMeanStandard Deviation
Area [ mm 2 ]6.82114.8959.7123.36
Feret Max. [mm]3.5216.9311.032.93
Feret Min. [mm]2.849.977.671.70
Aspect ratio [−]1.061.871.440.27
Circularity [−]0.460.770.630.08
Convexity [−]0.830.960.910.03
Elongation [−]1.612.752.120.33
Table 4. Shape parameters of investigated surface defects evaluated on the fracture surfaces of respective specimens.
Table 4. Shape parameters of investigated surface defects evaluated on the fracture surfaces of respective specimens.
ParameterMinimumMaximumMeanStandard Deviation
Area [ mm 2 ]2.9635.7222.037.16
Feret Max. [mm]4.2415.2910.012.60
Feret Min. [mm]1.265.173.680.93
Aspect ratio [−]1.345.162.721.05
Circularity [−]0.340.620.490.10
Convexity [−]0.760.950.890.06
Elongation [−]2.687.244.911.61
Table 5. Parameters of the cyclic stress–strain curve.
Table 5. Parameters of the cyclic stress–strain curve.
E [MPa] K [MPa] n [ ] b [ ] c [ ]
203,981892.560.0901 0.069 0.766
Table 6. Parameters of the combined hardening material model.
Table 6. Parameters of the combined hardening material model.
k [MPa] Q [MPa] b S [ ] C 1 [MPa] γ 1 [ ] C 2 [MPa] γ 2 [ ]
470 290.5 150.2899,4368996133,620522.7
Table 7. Average deviations obtained by the applied assessment methodologies.
Table 7. Average deviations obtained by the applied assessment methodologies.
2 α [ ] R [ ] Δ unsigned [ % ] Δ signed [ % ] Δ RMS [ % ]
1 3.9+0.34.5
4506.4 5.3 6.9
0.510.7 5.0 11.7
1 2.1+1.72.9
135010.9 10.6 12.0
0.56.9+5.910.7
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Horvath, M.; Oberreiter, M.; Stoschka, M. A Numerically Efficient Method to Assess the Elastic–Plastic Strain Energy Density of Notched and Imperfective Cast Steel Components. Appl. Mech. 2023, 4, 528-566. https://doi.org/10.3390/applmech4020030

AMA Style

Horvath M, Oberreiter M, Stoschka M. A Numerically Efficient Method to Assess the Elastic–Plastic Strain Energy Density of Notched and Imperfective Cast Steel Components. Applied Mechanics. 2023; 4(2):528-566. https://doi.org/10.3390/applmech4020030

Chicago/Turabian Style

Horvath, Michael, Matthias Oberreiter, and Michael Stoschka. 2023. "A Numerically Efficient Method to Assess the Elastic–Plastic Strain Energy Density of Notched and Imperfective Cast Steel Components" Applied Mechanics 4, no. 2: 528-566. https://doi.org/10.3390/applmech4020030

APA Style

Horvath, M., Oberreiter, M., & Stoschka, M. (2023). A Numerically Efficient Method to Assess the Elastic–Plastic Strain Energy Density of Notched and Imperfective Cast Steel Components. Applied Mechanics, 4(2), 528-566. https://doi.org/10.3390/applmech4020030

Article Metrics

Back to TopTop