1. Introduction
Contemporary transport engineering facilities operate at high speeds, high productivity, and high capacities to achieve the best performance. For aerospace and transport engineering, the equipment and facilities perform under high stress, which may result in elastic–plastic cyclic deformation. Overloading present particular dangers, as cyclically varying loads exceed the proportional limit of the material and cause plastic deformation and the formation of a hysteresis loop. As a result, the durability of the material decreases by hundreds or thousands of cycles.
A wise range of fatigue life prediction methods and probabilistic approaches, as well as mechanical and low cycle properties have been investigated in recent years. A considerable contribution to the calculation of probabilistic methods for mechanical and low-cycle properties was made by a series of investigators. Daunys et al. [
1,
2,
3,
4] investigated the dependences of the low-cycle durability of mechanical properties for steels of welded joints used in nuclear power plants. Ellingwood et al. [
5,
6] investigated the applicability of existing statistical data for describing the resistance of steel and reinforced concrete used in nuclear power plants. Liu et al. [
7] proposed calculating the equivalent initial flaw size (EIFS) distribution, which is very efficient for calculating the statistics of EIFS. Xiang et al. [
8] proposed a general probabilistic life prediction methodology for accurate and efficient fatigue prognosis, which is based on the inverse first-order reliability method (IFORM) to evaluate the fatigue life at an arbitrary reliability level. Bazaras et al. [
9] and Raslavičius et al. [
10] investigated the low-cycle durability of nuclear power plants’ WWER (Water–Water Energetic Reactor) steels 22k and 15Cr2MoVA. Zhu et al. [
11] developed a probabilistic methodology for low-cycle fatigue life prediction using an energy-based damage parameter under Bayes’ theorem. Fekete [
12] proposed a new low-cycle fatigue prediction model based on strain energy to account for only part of the strain energy stored in the microstructure of the material that causes fatigue damage. Strzelecki [
13] proposed the characteristics of the S–N curve using two-parameter and three-parameter Weibull distribution for fatigue limit and limited life. It was demonstrated that S–N curves can be used to determine the fatigue life for a low probability of failure when using a normal distribution. Kosturek et al. [
14] presented the results of their research on the low-cycle fatigue properties of Sc-modified AA2519-T62 extrusion. The basic mechanical properties have been established by using tensile tests and low-cycle fatigue testing has been performed on five different levels of total strain amplitude. Manouchehrynia et al. [
15] presented a mathematical model to estimate the strain-life probabilistic modelling based on the fatigue reliability prediction of an automobile coil spring under random strain loads. The obtained results demonstrated good agreement between the predicted fatigue lives of the proposed probabilistic model and the measured strain fatigue life models. Lamnauer et al. [
16] suggested the use of a probabilistic statistical model for calculating the strength of parts under cyclic fatigue loads. Statistical analysis of the samples (the average value, the corrected variance, the squared asymmetry coefficient, and the excess coefficient) was carried out according to the results of a mass experiment on the strength of samples during fatigue tests. Makhutov et al. [
17] analysed traditional engineering methods for the assessment of the lifetime characteristics of fatigue resistance. The methods used were based on deterministic parameters. The authors presented the results of experimental studies and the calculations of strength and durability for low-alloy and austenitic steels with varying mechanical properties.
Durability is one of the key criteria of structural elements. The application of appropriate probability calculation methods is important in the pursuit of extended life for in-service facilities. They also contribute to more accurate and research-based determinations of the safety values of cyclic loads at the design phase. Low-cycle strength and durability calculations based on the guaranteed mechanical characteristics rather than the standard ones retrieved from the specifications are necessary for the determination of the strength safety margin of structural elements. The strength safety margin of structural elements is, in turn, necessary for the assessment of the reliability of operation of the critical structures [
18,
19,
20].
The majority of the studies that undertake statistical assessment of low-cycle fatigue are focused on the uniaxial strain state and assessment of the durability distribution until the initiation of the fatigue crack or until the crack reaches a certain length. Currently, there are no consistent studies on the construction of probability curves for low-cycle fatigue in view of the values of the guaranteed mechanical characteristics [
21,
22,
23,
24,
25,
26].
Based on the topics discussed above, the main contributions of this paper are as follows: (1) We determine the distribution patterns of the mechanical properties, statistical parameters, and low-cycle fatigue curves; (2) we perform an analysis of the statistical assessment of cyclic elastoplastic strain diagrams and of the parameters; (3) we refine the low-cycle strength and durability calculations based on the verified values rather than standard values of key mechanical properties; and (4) we present a comparison of the low-cycle fatigue probability curves of the experimental data.
3. Identification of the Key Mechanical Properties and Correlations between Them
The standard (reference) key mechanical properties of materials were used as the input data for the calculation of the distribution under the probability methods for low-cycle strength and durability calculations [
30,
31,
32]. Hence, it was necessary to determine the preferable theoretical laws applicable to the experimental distribution functions of the key mechanical properties. The relationships between the mechanical properties simultaneously had to be determined in order to substantiate their values at a certain probability level.
Due to the very large number of results (on the order of hundreds) generated by the multiple tests, additional statistical processing—namely, statistical data series—was performed for the identification of the mechanical properties of the materials. For the histograms, the total array of the results was divided into 10 equal bins (statistical intervals). Their width was calculated using the following equation:
Exceeding 15–20 intervals would have been unreasonable due to the fact that even a very large number of results might still not ensure the accuracy of statistical characteristics. Following the division of the statistical data series into 10 intervals, the following histograms were developed (
Figure 2,
Figure 3 and
Figure 4).
Intervals equal in length were marked on the abscissa axis and the height of each interval was calculated using the following equation:
The histogram analysis showed the qualitative correspondence of the mechanical properties to the normal distribution law. Nonetheless, in order to improve the statistical assessment of the properties, statistical characteristics were calculated under the three applicable distribution laws: normal, log-normal, and Weibull distribution. The key statistical characteristics were calculated for the normal and log-normal distribution:
For the Weibull distribution [
33]:
Table 3 shows the calculated key statistical normal, log-normal, and Weibull distribution characteristics of the mechanical properties of the materials investigated.
Table 3 and the histograms (
Figure 2,
Figure 3 and
Figure 4) suggest the presence of a fairly large asymmetry of the majority of mechanical properties. The percent area reduction
ψ and percent area reduction at failure
ψu show the largest asymmetry. The coefficient of variation
V is one of the key statistical indicators. Its lowest values were obtained using the log-normal distribution. Hence, this distribution may be considered superior to normal or Weibull distribution. It should be noted that the values of both the coefficient of variation and other mechanical properties were similar in the Weibull and normal distribution. The lowest values obtained were those of the ultimate tensile stress
σu and cyclic stress
Sk of
k semicycle. The abscissa axis of
Figure 5 presents the coefficient of variation of Weibull and normal distribution laws and the ordinate axis—the coefficient of variation of log—normal distribution law.
The reliable quantitative assessment of mechanical properties is possible with a large sample available. In the case of a limited number of tests, the degree of accuracy and reliability must be provided, i.e., the confidence intervals, must be calculated (
Table 4):
For higher reliability of the probability calculations of the strength and durability of structural elements, it would be reasonable to use the calculated limit values of confidence intervals rather than the standard mechanical properties provided in the specifications [
34].
It is recommended to replace the mean value of the mechanical property with the lower endpoint of the confidence interval and the standard deviation of mechanical property with the upper endpoint of the confidence interval. Normally, when designing the facilities, standard values of mechanical properties are used. However, irrespective of the existing distribution of the properties, a considerable deviation from the true structural strength and durability is likely. In order to identify the error of the available statistical series, normalised values of mechanical properties were determined under the standard methodology and then compared to the standard values and arithmetic means of the materials under investigation.
In order to ensure a reliable safety margin, the mechanical properties were normalised from the bottom and the calculated norm values were determined using the variation series in ascending order by the variable. In this case, the calculated norm values were determined using the following equation:
The values of the statistical data density quantile
k1 depended on the sample size and materials investigated, namely, steel 15Cr2MoVa—1.41; steel C45—1.40; and aluminium alloy D16T1—1.43 [
28].
See
Table 5 for the calculations of the normalised mechanical properties performed at the required confidence level and under the procedure described above.
Along with the normalised mechanical properties,
Figure 6 presents the experimental and standard data of log-normal distribution of all the materials investigated.
Figure 6 suggests that the results of all the experimental mechanical properties are distributed linearly and this confirms the correspondence of the values to the log-normal distribution. Standard properties are mostly used when calculating the strain and durability diagram parameters: relative yield strength stress,
σ0.2; ultimate tensile stress,
σu; and relative percent area reduction,
ψ. These properties are usually provided in the reference sources.
A comparison of the experimental data, standard, and normalised properties (
σ0.2,
σu, and
ψ) is given in
Figure 6, which suggests that the values of the standard properties do not correspond to the experimental and normalised data. High probability is characteristic of the reference properties of the relative yield strength (
σ0.2) of the materials investigated: steel 15Cr2MoVa—74%; steel C45—62%; and aluminium alloy D16T1—90%. Meanwhile, the values of probability of the normalised mechanical properties
σ0.2 are considerably lower: steel 15Cr2MoVa—12%; steel C45—7.5%; and aluminium alloy D16T1—8%. Similar results were obtained for the stress (
σu) indicated in the ultimate strength standards for steel 15Cr2MoVa. The values of reference stress
σu of steel C45 and aluminium alloy D16T1 corresponded to a probability lower than 1%. Normalised
σu mechanical properties were: steel 15Cr2MoVa-25%; steel C45—105%; and aluminium alloy D16T1—8%.
The reference value of the relative percent area reduction ψ is 0.0003% for steel 15Cr2MoVa, while for steel C45, it corresponds to the experimental data, i.e., 50%. The normalised value for steel 15Cr2MoVa is 15%, for steel C45 it is 16%, and for aluminium alloy D16T1 it is 5%.
The key property defining the low-cycle durability is the relative percent area reduction ψ. Application of the comparative value ψ = 50% with 0.0003% probability for steel 15Cr2MoVa resulted in a very high durability safety factor compared to the standard safety factor. When the normalised value ψ = 76.72% with 15% probability was applied, the resulting deviation of the durability safety factor was fairly small.
For steel C45, a standard ψ value of 40–45%, corresponding to 50% probability, did not provide a sufficient confidence level. In this case, the normalised value ψ = 33.79% corresponded to a 10% probability.
The analysis of standard, normalised, and experimental results showed that the application of standard properties to the low-cycle fatigue calculations may lead to significant deviations from the actual results.
The distribution of a random event under the normal distribution law is known to be characterised by mean square deviation
s and dispersion
D. As the mean square deviation
s may be mathematically associated with the maximum and minimum values of the random measure, the ratio of the values may also be considered as the distribution characteristic:
Table 6 presents the
K values of measures
σpr,
σ0.2,
σu,
Sk,
ψ, and
ψu of the key mechanical properties.
According to
Table 6, high values of relative percent area reduction
ψu were seen for steel 15Cr2MoVa and C45. A comparison of the values of mechanical properties
σpr and
σu of the same materials showed that, for steel 15Cr2MoVa, the value of the proportionality limits
σpr of the
K coefficient was 26% higher than the ultimate strength
σu value. This was 35% for steel C45 and 16% for aluminium alloy D16T1. This was related to the considerable distribution of the properties of the proportionality limit. The diagrams of
K–
s and
K–V in
Figure 7 show the curves of the materials and mechanical properties investigated, described by the following equations [
27].
The resulting values enabled a primary assessment of the statistical properties
s,
V, and
according to the marginal values of mechanical properties usually provided in the material specifications. The resulting initial statistical characteristics may also be used for the determination of the minimum number of statistical specimens:
The analysis of calculation results of
na values (
Table 7) showed that the error Δ
a of determination of the mean value of a random measure had a considerable effect on the number of statistical tests.
A comparison of the calculated number of specimens with the values of the error Δa of determination of the mean value of a random measure (equal to 0.01–0.05) showed that the increase in Δa to 0.05 led to a 10-fold to 30-fold reduction in the number of specimens. In the same manner, the number of specimens was also affected by the reliability of normal distribution γ. The increase in its value from 0.05 to 0.1 led to a 1.5-fold increase in the number of specimens.
Single strain diagrams had a considerable effect on the low-cycle fatigue tests of the same material. The distribution of the single strain diagrams had a direct effect on the cyclic test diagrams and specimen durability [
1]. In
Figure 8 the absolute
σ–
e and relative
–
coordinates present the single strain diagrams of the limit values of the statistical series of the materials investigated.
An investigation of the tension diagrams suggested that the distribution of low-cycle test results had been affected by the type of loading and by the absolute or relative loading coordinates used.
In the case of loading with a controlled strain (
ek = constant), the effect was minor and depended on the level of loading. It could be observed in
Figure 8 that loading with controlled stress (
Sk = constant) was difficult to implement on the absolute coordinates. For steel 15Cr2MoVa, with the load being up to 400 MPa, the strain
e varied from 0.2% to 0.4%. Moreover, strain
e varied from 0.2% to 4.5%, where loading reached 450 MPa and strain
e varied from 0.2% to 11.5%. Where the loading level reached 500 MPa, strain
e varied from 0.2% to ∞.
The distribution of strain values on the relative coordinates decreased considerably. According to
Figure 8, for steel 15Cr2MoVa the value of relative strain
varied from 1 to 3 where
= 1.1. Meanwhile, where
= 1.4 the value of relative strain
ranged from 3.5% to 55%. Similar results were obtained when analysing the diagrams for steel C45 and aluminium alloy D16T1. For steel C45 (
Figure 8c), with the load being up to 300 MPa, strain
e varied from 0.2% to 6.5%. Where the load reached 500 MPa, strain
e varied from 0.2% to 8.5%, while when
σ = 650 MPa strain
e = 11.5–∞. Where the relative load was 1.1, strain
varied from 24% to 42% and where
= 1.4 strain
= 37–70% (
Figure 8d). Similar results were obtained for the aluminium alloy D16T1 (
Figure 8e,f).
4. Statistical Assessment of Low-Cycle Fatigue Curves
Probability values enabling the calculation of the theoretical low-cycle fatigue curves and their assessment from the probability perspective were determined for the mechanical properties already investigated statistically (
Table 8).
For a reliable statistical assessment of low-cycle fatigue durability properties, the paper investigated the distribution patterns and statistical parameters of the mechanical and low-cycle (strain and strength) properties of the materials with contrasting cyclic properties (hardening—aluminium alloy D16T1; softening—steel 15Cr2MoVa; and stable—steel C45).
The Coffin–Manson equation used in the strength calculations defines the dependence of durability under loading with controlled strain (
ek = constant) on the cyclic plastic strain
e0 [
35,
36]:
The modified Coffin–Manson equation was used in the present study:
In the equation α
1p <
m and
C1p <
ψ. Constants α
1p and
C1p may be determined using the mechanical properties of materials:
Manson–Langer power equations define cyclic resistance to failure and the dependence of durability between the elastoplastic strain
and number of cycles
N. Under low-cycle loading with controlled strain:
The study also employed the durability dependence presented in the design rules for the nuclear power industry, PNAE (Regularities and Norms in Nuclear Power Engineering) [
37]:
Using Equations (11), (13), and (14), low-cycle fatigue probability 1%, 10%, 30%, 50%, 70%, 90%, and 99% low-cycle durability curves were designed for steel 15Cr2MoVa, steel C45, and aluminium alloy D16T1 on the relative coordinates log
–log
Nc (
Figure 9 and
Figure 10). The relative values
of plastic strain were obtained by dividing the absolute strain values by the proportional limit strain
epr of the materials.
The comparison of experimental and theoretical curves for steel 15Cr2MoVa presented in
Figure 9 showed that the curve slope angle was similar for all cases. However, the experimental curves were slightly lower than the theoretical ones. The theoretically calculated curves fell within the experimental curve zone (
Figure 9a); however, they were hardly comparable due to the specifics of the calculation of probability constants α
1p and
C1p with a small number of cycles. With the number of cycles
N > 400, the experimental 99% probability curve corresponded to the 50% theoretical curve (
Figure 9a). The theoretical curves calculated under the PNAE rules (
Figure 9c) fell between the experimental curves. In all the calculations, the resulting arrangement was the reverse. The 99% experimental curve corresponded to the 1% theoretical curve, etc. It could be assumed that this resulted from the dependence of constant
α1p on the relative percent area reduction
ψ. According to
Table 8, the ratio of proportional limit strain
epr 99% to 1% values was 5.3:1 and for the relative percent area reduction it was 1.2:1. Moreover, the proportional limit strain
epr values were sensitive to variations in chemical composition, thermal processing technologies, surface hardening, loading conditions, and other factors of the material.
The comparison of the experimental and theoretical curves of steel C45 under loading with controlled strain (
ek = constant) in
Figure 10a suggests that, in all cases, the resulting curve slope was similar. The theoretically calculated curves were lower than the experimental ones; however, in this case, the probability arrangement of the curves was not the reverse. The 99% to 1% durability curve ratio is 7.1:1 at the relative strain amplitude
= 4% calculated by Equation (11), 7.6:1 according to Equation (12), and 10.3:1 according to Equation (12). In this case, with the relative strain amplitude
= 2%, the ratios of the durability curves were 8.4:1, 11.7:1, and 5.3:1.
Figure 10c suggests that the calculation under the PNAE rules, Equation (15), had the best correspondence with the experimental results.
The results for the D16T1 aluminium alloy under loading with controlled strain (
ek = constant) are presented in
Figure 11. A comparison of the 99% and 1% probability curves in
Figure 11a showed the clear dependence of the low-cycle durability on the strain level. The conducted analysis suggested that the 99% to 1% durability curve ratio was 37:1, where the strain amplitude
e0 = 0.3%, and 24:1 where the strain amplitude
e0 = 0.18%. The slope of the theoretical curves increased with the increase in the low-cycle failure probability. This could be related to the percent area reduction
ψ distribution (
Table 3). The
distribution band became narrower when relative coordinates were used. The 99% to 1% durability curve ratio was 3.3:1 when the strain amplitude
= 4 and 2.7:1 when the strain amplitude
= 3. The slope angles were smaller in the relative coordinate curves.
Figure 11b presents a comparison of the experimental and theoretical curves. The experimental and theoretical results differed considerably, with the theoretical curves being in the elasticity zone.