Fatigue Damage Accumulation Modeling of Metals Alloys under High Amplitude Loading at Elevated Temperatures

: This article presents an approach related to the modeling of the fatigue life of constructional metal alloys working under elevated temperature conditions and in the high-amplitude load range. The article reviews the fatigue damage accumulation criteria that makes it possible to determine the number of loading cycles until damage occurs. Results of experimental tests conducted on various technical metal alloys made it possible to develop a fatigue damage accumulation model for the LCF (Low Cycle Fatigue) range. In modeling, the material’s damage state variable was deﬁned, and the damage accumulation law was formulated incrementally so as to enable the analysis of the inﬂuence of loading history on the material’s fatigue life. In the proposed model, the increment of the damage state variable was made dependent on the increment of plastic strain, on the tensile stress value in the sample, and also on the actual value of the damage state variable. The model was veriﬁed on the basis of data obtained from experiments in the ﬁeld of uniaxial and multiaxial loads. Samples made of EN AW 2024T3 aluminum alloy were used for this purpose.


Introduction
During operation, machine parts are frequently subjected to fatigue loading at elevated temperatures. Such operating conditions have an adverse effect on the materials' behavior. An elevated temperature accelerates wear processes and contributes to increasing the rate of fatigue crack propagation [1,2]. Aluminum alloys, e.g., EN AW 2024T3, are among the materials used to make constructions capable of working at an elevated temperature. It is widely used in the motorization and aviation industries, which place strong emphasis on ensuring high strength and light weight at the same time [3]. These features make it suitable for the manufacture of parts working in close proximity to combustion and jet engines where elevated temperature zones are present. In the work of Szusta and Seweryn [4,5], the influence of elevated temperature on the fatigue life of this material was tested. The authors demonstrated that strength parameters decreased as the test temperature increased. Based on the tests conducted, the authors determined the parameters of the Manson-Coffin fatigue life curve as a function of temperature: where N f is the number of cycles until failure; E(T) is Young's modulus determined at temperature T; curve for the analyzed temperature T; and b(T) and c(T) are, respectively, exponents of the elastic and plastic fatigue life curve for temperature T. Furthermore, based on observations of the material's cracking mechanisms at an elevated temperature, the authors proposed a semi-empirical model for estimating the material's fatigue life. In this model, it was accepted that the material's fatigue life drops as the test temperature (isothermal) increases. The lower the material's melting temperature, the more substantial this drop is. Moreover, it was assumed that there exists a limit temperature T m at which the given material no longer has any immediate strength. Therefore, knowing the value of this temperature and the function defining the reduction in fatigue life as temperature increases, calculations of the number of cycles until crack initiation can be performed with engineering accuracy for a given temperature level and strain amplitude. For uniaxial tensile-compressive loads, the function of the temperature's influence on the material's fatigue life has been defined as follows: where A and B are material constants determined in periodical tensile-compression tests at temperature T.
In this model, the material constant B is a function of the material's recrystallization temperature T m (B = −AT m ). The value of parameter A is calculated on the basis of the experimentally determined number of cycles to failure (obtained for a given load course at room temperature) related to the value of temperature difference: the material recrystallization and room temperature. Parameter A is determined from the following equation: where N f RT, ∆ε eq 2 is determined from Manson-Coffin-Basquiun equation at room temperature RT for the set total strain amplitude ∆ε eq 2 . Modified cast irons and cast steels, which are used in the manufacture of brake disks, are another group of materials that can be exploited at elevated temperatures. In the works of Samec et al. [6], Peve et al. [7], and Li et al. [8], the authors analyzed the fatigue life of the following materials: EN-GJS-500-7, EN-GJL-250, and Cr-Mo-V. Material cracking mechanisms were studied within the LCF range at temperature up to 700 • C, and the influence of temperature on the material's physical properties was defined.
Alloy steels are another group of materials used to make constructions which are exposed to the action of elevated temperatures. Here, some examples include 8Cr-2WvTa, (RAFM) JLF-1, or 617M steels, which are used to build thermoreactors and machine parts exposed to the action of high temperatures. In the works of Ishii et al. [9], Mariappan et al. [10], and Shankar et al. [11], the characteristics of these steels' cyclic properties were determined as a function of temperature. Degradation mechanisms of materials subjected to the action of fatigue loads and temperatures were analyzed. In the article by Ishii et al. [9], an approach to modeling fatigue life of the material was also presented, accounting for the influence of creep on the material's fatigue life: where ε eq is the equivalent strain, ε e is the elastic strain, ε pl is the plastic strain, and ε cr is the creep-induced strain. Works by Nagode and Zingsheim [12] and Nagode and Hack [13] present the results of experimental tests performed on steels for work at elevated temperatures: 10 CrMo 9 10 and X22CrMoV121. Fatigue damage accumulation models of the material under non-isothermal conditions at elevated temperatures were developed on the basis of these results. These models were developed for estimation of the fatigue life of materials which undergo changes of the load amplitude and temperature transient over the course of its exploitation. These models do not account for the influence of creep or hardening of the material on its fatigue life. The authors presented a stress and strain approach in which it is assumed that the instant working temperature will be accounted for. In models, stabilized hysteresis loops were described by means of the Ramberg-Osgood equation. The damage parameter P SWT was identified similarly as in the SWT (Smith-Watson-Topper) model: where T e is the actual temperature; σ f (T e ) and b(T e ) are, respectively, the coefficient and exponent of the fatigue life curve at the tested temperature; K (T e ) and n (T e ) are, respectively, the coefficient and exponent of the cyclic hardening curve at the tested temperature; and ε f (T e ) and c(T e ) are, respectively, the coefficient and exponent of the plastic fatigue life curve. Elements of gas turbines, e.g., rotor blades working under extremely difficult conditions, are also made from Inconel 718 alloy. Test results for this material was given in the work of Schlesinger et al. [14]. Here, a dispersion-hardening nickel-chromium alloy was analyzed with respect to fatigue loading at temperatures up to 650 • C. Similarly, in a different study [15], the same authors studied the Inconel 718 alloy. They proposed a model for estimation of the fatigue life of material working at elevated temperatures. The Manson-Coffin equation was used during modeling, in which parameters of the fatigue life curve were made functionally dependent on the structure and grain size in the given material G, test temperature T, and strength parameter S: Materials adapted for manufacturing tools like casting molds and forging dies play an important role among materials intended for work at elevated temperature. In the study by Tunthawiroon et al. [16], the single-phase Co-29Cr-6Mo steel alloy, intended for aluminum casting molds, was tested. The study included analysis of the alloy's cracking mechanisms at temperatures up to 700 • C. A simple fatigue damage accumulation model at elevated temperatures was also proposed based on the Arrhenius equation: where N f is the number of cycles until crack initiation, T is the test temperature, R and A are material constants, and Q is the activation energy. Similarly, the works of Gopinath et al. [17] and He et al. [18] present the results of experimental tests conducted on 720Li and HAYNES HR-120 superalloys within the low-cycle fatigue range at temperatures up to 980 • C. The material was tested under uniaxial tensile-compressive loads. The mechanisms accompanying the cracking of these materials were defined.
ACI HH50 austenitic stainless steels can also work at elevated temperatures. The results given in the work by Kim and Jang [19] can serve as evidence of this. Due to its high resistance to elevated temperatures, high fatigue strength, and resistance to pitting and corrosion, the material in question finds applications in parts of combustion engines and power unit structures of nuclear power plants. The paper investigates the material's degradation mechanisms and presents a model estimating fatigue life in the LCF range. The model assumes that fatigue crack initiation will occur when damage reaches a limit value given by the following function: where α is the material constant associated with dislocation movement over the course of the material loading process, ν is the frequency of load change during the test, .
ε 0 and C are material constants independent of test temperature, Q is activation energy, R is the gas constant, and T is the test temperature.
The non-linear creep damage accumulation model under uniaxial loads at elevated temperatures (700 • C) that was also developed on the basis of tests performed on stainless steel X-8-CrNiMoNb-16-16 was presented by Pavlou 2001 [20]. The model was developed and verified by other researchers. In the work of Grell et al. [21], attempts have been made to estimate durability under the conditions of uniaxial isothermal loads (constant and incremental) at elevated temperatures using the dependencies proposed by Pavlou 2001 [20]. The work demonstrates greater accuracy in predicting durability with the use of non-linear models of the accumulation of defects in relation to linear models.
Titanium alloys are the next group of materials used to manufacture components working at elevated temperatures. One study [22] analyzes the influence of different methods of the surface hardening of the Ti-6Al-4V titanium alloy on the fatigue life of this material under uniaxial tensile-compressive loading conditions at temperatures up to 555 • C. The influence of temperature on the evolution of the cyclic properties of the TNB-V2 titanium alloy is also presented in the work [23]. Specimens of the material were subjected to constant-amplitude strains within the low-cycle loading range at temperatures within the range of 550-850 • C. Due to its good mechanical properties at elevated temperatures, this alloy finds applications in motorization and aviation and in engine parts. The authors tested the material within the entire safe temperature range, and determined that, in the case of this alloy, the Manson-Coffin fatigue life model does not allow for the prediction of the material's fatigue life at elevated temperatures due to the length of the phase between fatigue crack initiation and its propagation.
The fatigue damage accumulation model proposed in this paper was designed to predict the fatigue life of the material operated at elevated temperatures. The concept of elevated temperature can be defined here, according to the works by Chen et al. [24], as the temperature corresponding to two-thirds of the melting point of the bulk material T c . In most cases, as temperature increases, a material's strength and fatigue properties are reduced in a predictable manner until a certain threshold temperature value is reached [25]. Under these conditions, the operation of technical machinery may be safe when the material's response to loading conditions and the working environment are known. However, if the conventional upper threshold of elevated temperature is crossed even slightly, the material's properties cannot be predicted directly. The material's strength suddenly drops (in a non-linear manner), and use of machinery under such conditions may pose a safety risk.
Uncontrolled loss of durability may lead to costly failures, standstills on process lines, and sometimes even to the death of those nearby when a structural component cracks. This is why information about how long and under what loading and temperature conditions a structural component can work safely is important. The fatigue damage accumulation models of materials working under low-cycle loading conditions at elevated temperatures presented in the literature are modified functions formulated to estimate the material's fatigue life at room temperature. The equations mentioned earlier can serve as an example: Manson-Coffin-Basquin, Ramberg-Osgood, and Smith-Watson-Topper. In these models, material parameters are determined independently for each of the analyzed temperatures. The process of calculating these parameters is not complicated in itself; however, preparing the data for determining them is very expensive and time-consuming and requires specialized testing apparatus. These models are dedicated for specific loading cases and metal alloys. Due to the low number of functioning fatigue criteria in this scope, there is a justified need to create new models that will provide a more accurate description of the material cracking process, thus enabling calculation of fatigue life with greater accuracy. That is also why it this paper set out to create a new fatigue damage accumulation model that will make it possible to determine the number of cycles until failure of a material working under conditions of fatigue loads of constant amplitude at elevated temperatures with engineering accuracy.

Fatigue Life of EN AW-2024T3 Aluminum Alloy at Elevated Temperatures
To develop a fatigue damage accumulation model, it is indispensable to perform experimental tests. They constitute the physical basis of formulated numerical dependencies. The fatigue damage accumulation model proposed in this paper was created based on an experiment in which EN AW-2024T3 aluminum alloy samples were subjected to the action of high-amplitude fatigue loading in the form of uniaxial tension-compression and complex loading (tension-compression with simultaneous torsion) at elevated temperatures. More information about the experiment can be found in other works [4,5].
Experimental working under low-cycle loading conditions at elevated temperatures presented in the literature are modified functions formulated to estimate the material's fatigue life at room temperature. The equations mentioned earlier can serve as an example: Manson-Coffin-Basquin, Ramberg-Osgood, and Smith-Watson-Topper. In these models, material parameters are determined independently for each of the analyzed temperatures. The process of calculating these parameters is not complicated in itself; however, preparing the data for determining them is very expensive and time-consuming and requires specialized testing apparatus. These models are dedicated for specific loading cases and metal alloys. Due to the low number of functioning fatigue criteria in this scope, there is a justified need to create new models that will provide a more accurate description of the material cracking process, thus enabling calculation of fatigue life with greater accuracy. That is also why it this paper set out to create a new fatigue damage accumulation model that will make it possible to determine the number of cycles until failure of a material working under conditions of fatigue loads of constant amplitude at elevated temperatures with engineering accuracy.

Fatigue Life of EN AW-2024T3 Aluminum Alloy at Elevated Temperatures
To develop a fatigue damage accumulation model, it is indispensable to perform experimental tests. They constitute the physical basis of formulated numerical dependencies. The fatigue damage accumulation model proposed in this paper was created based on an experiment in which EN AW-2024T3 aluminum alloy samples were subjected to the action of high-amplitude fatigue loading in the form of uniaxial tension-compression and complex loading (tension-compression with simultaneous torsion) at elevated temperatures. More information about the experiment can be found in other works [4,5].
Experimental The scope of performed tests covered, first of all, monotonic tests determining the influence of elevated temperatures on the mechanical properties of the EN AW-2024 T3 aerospace aluminum alloy. Table 1 presents the obtained results. The scope of performed tests covered, first of all, monotonic tests determining the influence of elevated temperatures on the mechanical properties of the EN AW-2024 T3 aerospace aluminum alloy. Table 1 presents the obtained results. Next, fatigue tests were carried out. Samples were subjected to uniaxial loading, periodically varying (tension-compression), until their failure at different set temperatures. Tests were performed at temperatures of 20, 100, 200, and 300 • C and at a constant value of the strain change range ε, frequency f = 1 Hz, and R ε = −1. The following values of control variable ε were applied: 0.015, 0.01, 0.0095, 0.008, and 0.006. The tests were performed on three samples for every value of the control variable's range and every test temperature.
Based on the performed tests, parameters describing the process of cyclic deformation of the material (as expressed by the Ramberg-Osgood equation) at elevated temperatures were determined by the following equation: where σ a is the amplitude of normal stress induced by the action of periodic axial force, K(T) and n(T) are the coefficient and exponent of Ramberg-Osgood cyclic strain curve, which are dependent on temperature. Table 2 presents the calculated parameters of the cyclic strain curve. In the next stage, parameters of the Manson-Coffin-Basquin strain-based fatigue life curve were determined. Its parameters, calculated for various temperatures, are given in Table 3. Table 3. Coefficients and factors of the fatigue curve for the EN AW-2024T3 alloy at different temperature values [4].
Next, samples were subjected to the action of cyclic loading with torque moment and tensile or compressive force (during loading of samples, a constant quotient of maximum longitudinal and shear strain was assumed over the course of a loading cycle: ε a /γ a = √ 3). In all cases, loading histories formed oriented, multi-segment loops. Their configuration and parameters are given in Figure 2. Just as in the case of uniaxial tests, the loading process was controlled by means of increments of strain tensor components. Tubular specimens (Figure 1b) were used in the experiment, and they were subjected to proportional and non-proportional tension-torsion cyclic loading with cycle asymmetry factor R = 0 (εmin/εmax; γmin/γmax). Fatigue tests were conducted for different combinations of axial and shear strains as well as for four temperature values: 20, 100, 200, and 300 °C . The specimen's kinematic input was controlled by means of the MTS 632.68F-08 (Eden Prairie, MN, USA) biaxial extensometer, and increments of strain components (linear and shear), averaged over a 25 mm-long segment of the measuring base, were used for this purpose. Three tests were conducted for every load with frequency of load changes f = 1 Hz.
The crack initiation moment was determined on the basis of the analysis of the maximum force value in the loading cycle. A 10% reduction of the force value relative to its maximum value at a given strain level was accepted as the criterion for crack initiation.
A compilation of the results of experimental tests (number of loading cycles until crack initiation) is given in Table 4.
The registered hysteresis loops for selected multiaxial loading configurations (individual loading components) have been presented in Figure 3. Tubular specimens (Figure 1b) were used in the experiment, and they were subjected to proportional and non-proportional tension-torsion cyclic loading with cycle asymmetry factor R = 0 (ε min /ε max ; γ min /γ max ). Fatigue tests were conducted for different combinations of axial and shear strains as well as for four temperature values: 20, 100, 200, and 300 • C. The specimen's kinematic input was controlled by means of the MTS 632.68F-08 (Eden Prairie, MN, USA) biaxial extensometer, and increments of strain components (linear and shear), averaged over a 25 mm-long segment of the measuring base, were used for this purpose. Three tests were conducted for every load with frequency of load changes f = 1 Hz.
The crack initiation moment was determined on the basis of the analysis of the maximum force value in the loading cycle. A 10% reduction of the force value relative to its maximum value at a given strain level was accepted as the criterion for crack initiation.
A compilation of the results of experimental tests (number of loading cycles until crack initiation) is given in Table 4.
The registered hysteresis loops for selected multiaxial loading configurations (individual loading components) have been presented in Figure 3.

Modeling of Damage Accumulation and Fracture of Material at Elevated Temperatures
Preliminary tests based on the cyclic loading of a flat sample and observation of the amount of energy dissipating from the material due to the action of plastic strains were performed in order to investigate the fatigue damage accumulation process. The energy released as a result of fatigue loading was in the form of heat and was observed through a thermal imaging camera. Tests were performed at room temperature on the stand shown in Figure 4.

Modeling of Damage Accumulation and Fracture of Material at Elevated Temperatures
Preliminary tests based on the cyclic loading of a flat sample and observation of the amount of energy dissipating from the material due to the action of plastic strains were performed in order to investigate the fatigue damage accumulation process. The energy released as a result of fatigue loading was in the form of heat and was observed through a thermal imaging camera. Tests were performed at room temperature on the stand shown in Figure 4. During the tests, changes in the deformation fields (Aramis 4M-DIC system) and temperature fields (Cedip Titanium-IRT system) were recorded simultaneously. On their basis, areas of plastic deformation were identified and the distribution of plastic strain and energy converted into heat was determined. As can be seen in Figure 5, the temperature distribution corresponds to the strain distribution. In connection with the above, it can be assumed that the change in the temperature of the object caused by the action of external load is an indicator of damage generated during the loading process.
An example of temperature and strain fields is presented in Figure 5. Those are the fields corresponding to the localization of the plastic strain. This localization was manifested by the nonuniform temperature field and by the heterogeneity of the strain field on the specimen surface. The evolution of the temperature field due to the plastic deformation was used to determine the evolution of the field of heat dissipated by the specimen. The distribution of the heat and the plastic work on the surface of the deformed specimen was used to determine the failure mechanism.
The amount of heat emitted over the course of the tests can be related to the number of faults (damage) appearing in the material as a result of the action of external load. As the number of damages in the material increases, the object's temperature increases. Figure 6 presents the function of the maximum sample temperature at the crack initiation site as a function of test time (number of loading cycles). During the tests, changes in the deformation fields (Aramis 4M-DIC system) and temperature fields (Cedip Titanium-IRT system) were recorded simultaneously. On their basis, areas of plastic deformation were identified and the distribution of plastic strain and energy converted into heat was determined. As can be seen in Figure 5, the temperature distribution corresponds to the strain distribution. In connection with the above, it can be assumed that the change in the temperature of the object caused by the action of external load is an indicator of damage generated during the loading process.
An example of temperature and strain fields is presented in Figure 5. Those are the fields corresponding to the localization of the plastic strain. This localization was manifested by the non-uniform temperature field and by the heterogeneity of the strain field on the specimen surface. The evolution of the temperature field due to the plastic deformation was used to determine the evolution of the field of heat dissipated by the specimen. The distribution of the heat and the plastic work on the surface of the deformed specimen was used to determine the failure mechanism.
The amount of heat emitted over the course of the tests can be related to the number of faults (damage) appearing in the material as a result of the action of external load. As the number of damages in the material increases, the object's temperature increases. Figure 6   As can be seen in the initial phase of the loading process, damage accumulated according to the linear hypothesis of damage summation. At a certain time in the loading process, there was a sharp spike in the amount of damage, observed as an increase in released thermal energy. Damage then increased exponentially. In the tensile half-cycle of loading, the object's temperature increased as a result of growing damage. Meanwhile, the compressive half-cycle slowed down the damage accumulation process, resulting in a drop of the sample's temperature.
Based on observations of fatigue cracking planes obtained over the course of experimental tests, it was noted that slip lines form first on the samples' polished surfaces under the influence of cyclic loadings at elevated temperature, even with values less than that of fatigue strength, and they transform into fatigue slip bands (numerous parallel faults) as the number of loading cycles grew (Figure 7).  As can be seen in the initial phase of the loading process, damage accumulated according to the linear hypothesis of damage summation. At a certain time in the loading process, there was a sharp spike in the amount of damage, observed as an increase in released thermal energy. Damage then increased exponentially. In the tensile half-cycle of loading, the object's temperature increased as a result of growing damage. Meanwhile, the compressive half-cycle slowed down the damage accumulation process, resulting in a drop of the sample's temperature.
Based on observations of fatigue cracking planes obtained over the course of experimental tests, it was noted that slip lines form first on the samples' polished surfaces under the influence of cyclic loadings at elevated temperature, even with values less than that of fatigue strength, and they transform into fatigue slip bands (numerous parallel faults) as the number of loading cycles grew (Figure 7). As can be seen in the initial phase of the loading process, damage accumulated according to the linear hypothesis of damage summation. At a certain time in the loading process, there was a sharp spike in the amount of damage, observed as an increase in released thermal energy. Damage then increased exponentially. In the tensile half-cycle of loading, the object's temperature increased as a result of growing damage. Meanwhile, the compressive half-cycle slowed down the damage accumulation process, resulting in a drop of the sample's temperature.
Based on observations of fatigue cracking planes obtained over the course of experimental tests, it was noted that slip lines form first on the samples' polished surfaces under the influence of cyclic loadings at elevated temperature, even with values less than that of fatigue strength, and they transform into fatigue slip bands (numerous parallel faults) as the number of loading cycles grew (Figure 7). Similarly, as in the case of using room temperature [26], a fatigue damage accumulation criterion can be proposed based on the assumption that the material's damage state mainly depends on variable slips on physical planes and the normal stresses on these planes.
Considering the above, this paper proposes a damage accumulation model intended for the analysis of the fatigue life of constructional elements working under uniaxial and multiaxial loading conditions at elevated temperatures. This model assumes that permanent slip bands are the site of crack nucleation and defect formation [26], where crack initiation is the result of the accumulation of this damage. In this case, damage is induced by the action of the stress vector's normal component. Here, it is accepted that tensile stress in the loading cycle is responsible for generating new faults and accelerating the development of existing faults (temperature increase in loading cycle; Figure 6), while compressive stress retards the damage accumulation process (temperature drop in loading cycle; Figure 6). This fact was accounted for in the definition of the damage state variable dωn.
The presented law of damage accumulation induced by plastic strains under elevated temperature conditions was linked to the stress-based damage accumulation function Ψp, the value of which is dependent on the value of normal stresses σn, damage state variable ωn on the physical plane, and temperature T. The increment of the damage state variable on the physical plane dωn, as caused by the development of plastic strains at elevated temperature T, was made dependent on the increment of plastic shear strains dγn p on the same plane [27], namely: (10) where Ap(T) is the material variable describing the evolution of the material's plastic properties depending on the actual stress state and temperature. In numerical simulations, the following values of parameter Ap were accepted: Ap(20 °C ) = 0.55; Ap(100 °C ) = 1.4; Ap(200 °C ) = 1.8; and Ap(300 °C ) = 1.95.
The damage accumulation mechanism described by this law is presented in Figure 8. Similarly, as in the case of using room temperature [26], a fatigue damage accumulation criterion can be proposed based on the assumption that the material's damage state mainly depends on variable slips on physical planes and the normal stresses on these planes.
Considering the above, this paper proposes a damage accumulation model intended for the analysis of the fatigue life of constructional elements working under uniaxial and multiaxial loading conditions at elevated temperatures. This model assumes that permanent slip bands are the site of crack nucleation and defect formation [26], where crack initiation is the result of the accumulation of this damage. In this case, damage is induced by the action of the stress vector's normal component. Here, it is accepted that tensile stress in the loading cycle is responsible for generating new faults and accelerating the development of existing faults (temperature increase in loading cycle; Figure 6), while compressive stress retards the damage accumulation process (temperature drop in loading cycle; Figure 6). This fact was accounted for in the definition of the damage state variable dω n .
The presented law of damage accumulation induced by plastic strains under elevated temperature conditions was linked to the stress-based damage accumulation function Ψ p , the value of which is dependent on the value of normal stresses σ n , damage state variable ω n on the physical plane, and temperature T. The increment of the damage state variable on the physical plane dω n , as caused by the development of plastic strains at elevated temperature T, was made dependent on the increment of plastic shear strains dγ n p on the same plane [27], namely: where A p (T) is the material variable describing the evolution of the material's plastic properties depending on the actual stress state and temperature. In numerical simulations, the following values of parameter A p were accepted: A p (20 • C) = 0.55; A p (100 • C) = 1.4; A p (200 • C) = 1.8; and A p (300 • C) = 1.95. The damage accumulation mechanism described by this law is presented in Figure 8. Similar to the dependencies presented in the work [26], the damage accumulation function Ψp can be proposed in the following form: where Rσ is the stress function of cracking given by the condition of maximum normal stresses, which is dependent on the test temperature:  Table 3).
The crack initiation criterion can be written in the form of the condition of maximum normal stresses (linked to the physical plane), namely: where Rfσ is the stress cracking coefficient dependent on the state of stress and damage in the material as well as on the test temperature T, and σc is the actual value of normal failure stresses for the material at the given temperature, which is dependent on the damage state variable and temperature, namely: where σc0(T) is the critical stress for undamaged material, which is dependent on temperature. The value of parameter σc0(T) can be equated with the value of the coefficient of the Manson-Coffin fatigue life curve σf′(T). Critical stress values according to Table 3 were used to calculate the accumulation of fatigue damage by means of the proposed model. An additional cracking condition should also be introduced, in which it is accepted that crack initiation will occur when the damage state variable induced by plastic strains on any physical plane reaches critical value, that is: At that moment, the material has no strength (σc = 0). In the presented model, stress and strain values were determined independently for each test temperature by using the generalized Hooke's law, the Huber-von Mises yield criterion, the gradient flow law, and Mroz-Garud multiple-surface material hardening model associated with the yield criterion. They are an integral part of the numerical model of damage accumulation, making it possible to determine stress and strain state in any loading state [28].
When estimating the material's fatigue life, the principle of additivity of components of elastic and plastic strain increment tensors was adopted. Increments of elastic strain were determined by Similar to the dependencies presented in the work [26], the damage accumulation function Ψ p can be proposed in the following form: where R σ is the stress function of cracking given by the condition of maximum normal stresses, which is dependent on the test temperature: c(T) is the exponent present in the Manson-Coffin-Basquin equation, which is also dependent on temperature (values of parameter c(T) used in numerical simulations are given in Table 3).
The crack initiation criterion can be written in the form of the condition of maximum normal stresses (linked to the physical plane), namely: where R fσ is the stress cracking coefficient dependent on the state of stress and damage in the material as well as on the test temperature T, and σ c is the actual value of normal failure stresses for the material at the given temperature, which is dependent on the damage state variable and temperature, namely: where σ c0 (T) is the critical stress for undamaged material, which is dependent on temperature. The value of parameter σ c0 (T) can be equated with the value of the coefficient of the Manson-Coffin fatigue life curve σ f (T). Critical stress values according to Table 3 were used to calculate the accumulation of fatigue damage by means of the proposed model. An additional cracking condition should also be introduced, in which it is accepted that crack initiation will occur when the damage state variable induced by plastic strains on any physical plane reaches critical value, that is: max At that moment, the material has no strength (σ c = 0). In the presented model, stress and strain values were determined independently for each test temperature by using the generalized Hooke's law, the Huber-von Mises yield criterion, the gradient flow law, and Mroz-Garud multiple-surface material hardening model associated with the yield criterion. They are an integral part of the numerical model of damage accumulation, making it possible to determine stress and strain state in any loading state [28].
When estimating the material's fatigue life, the principle of additivity of components of elastic and plastic strain increment tensors was adopted. Increments of elastic strain were determined by means of the generalized Hooke's law, and increments of plastic strain were calculated by means of the gradient flow law, associated with the yield surface, namely (T = const): where both the proportionality coefficient dλ and the elasticity constants E and ν are dependent on temperature. In numerical simulations, the values of parameters determined in monotonic tensile tests were used according to Table 1. The material's hardening curve was plotted according to the parameters given in Table 2. Yield surface f = 0 was determined using the Huber-von Mises criterion: where R(T) defines the temperature-dependent size of the yield surface.

Results
Simulations of the fatigue life of the aerospace aluminum alloy within the range of uniaxial loads (tension-compression) at elevated temperatures were conducted for the modified Manson-Coffin-Basquin model (Equation (1)), semi-empirical model S-S_1 (Equation (2)), and the proposed approach S-S_2 (Equation (10)). Damage accumulation simulations were carried out for four different temperatures and five strain levels. The results of the simulations are provided in Table 5.  Figure 9a presents the evolution of the damage state variable over the course of the fatigue loading process, plotted based on Equation (10) for the analyzed levels of temperature and strain ε a = 0.015. As can be seen in its evolution, the fatigue damage accumulation process progresses in a near-linear manner in the first loading stage until it reaches a certain limit state, and after this limit state is crossed, there is a sudden spike in damage. This was confirmed by thermograms (Figure 9b), which show a sudden temperature increase in the final phase of periodic loading. The characters of the fatigue damage accumulation process and of the evolution of the sample's temperature were similar over the course of cyclic loading. It can be surmised that the object's temperature may be an indicator of the number of damage accumulated over the course of the material loading process.  Figure 9a presents the evolution of the damage state variable over the course of the fatigue loading process, plotted based on Equation (10) for the analyzed levels of temperature and strain εa = 0.015. As can be seen in its evolution, the fatigue damage accumulation process progresses in a nearlinear manner in the first loading stage until it reaches a certain limit state, and after this limit state is crossed, there is a sudden spike in damage. This was confirmed by thermograms (Figure 9b), which show a sudden temperature increase in the final phase of periodic loading. The characters of the fatigue damage accumulation process and of the evolution of the sample's temperature were similar over the course of cyclic loading. It can be surmised that the object's temperature may be an indicator of the number of damage accumulated over the course of the material loading process. Temperature influences the rate at which the material subjected to the action of cyclic loads reaches the state of "damage saturation". Samples loaded at 300 °C reached the limit state after working 0.26 of the number of cycles until failure, while at room temperature, this state was reached after 0.38 of the number of cycles until failure. This state can be explained by the fact that movement of dislocations is facilitated as temperature increases, since they do not encounter obstacles that could slow them down on their path due to the elevated mobility of atoms. After crossing the "limit" state, further damage accumulation progresses exponentially. Microcracks appear in the material, leading to the material's decohesion as they are joined.  Temperature influences the rate at which the material subjected to the action of cyclic loads reaches the state of "damage saturation". Samples loaded at 300 • C reached the limit state after working 0.26 of the number of cycles until failure, while at room temperature, this state was reached after 0.38 of the number of cycles until failure. This state can be explained by the fact that movement of dislocations is facilitated as temperature increases, since they do not encounter obstacles that could slow them down on their path due to the elevated mobility of atoms. After crossing the "limit" state, further damage accumulation progresses exponentially. Microcracks appear in the material, leading to the material's decohesion as they are joined. Figure 10 presents a comparison of numerical calculations of fatigue life (for the modified Manson-Coffin-Basquin model (Equation (1)), semi-empirical model S-S_1 (Equation (2)), and the proposed approach S-S_2 (Equation (10) Next, the proposed fatigue damage accumulation model was also verified for the case of multiaxial loads. The results of original experimental tests performed on the EN AW-2024T3 aerospace aluminum alloy were used for this purpose [5]. Figure 11 presents normalized evolutions of the damage state variable, plotted for load paths RS0 and RS90 at the analyzed temperature. A comparison between the results of experimental tests and of numerical simulations using the proposed SS2 approach (Equation (10)) and the SS1 semi-empirical model (Equation (2)) described in Reference [5] is given in Figure 12 and in Table 6. Next, the proposed fatigue damage accumulation model was also verified for the case of multiaxial loads. The results of original experimental tests performed on the EN AW-2024T3 aerospace aluminum alloy were used for this purpose [5]. Figure 11 presents normalized evolutions of the damage state variable, plotted for load paths RS0 and RS90 at the analyzed temperature. A comparison between the results of experimental tests and of numerical simulations using the proposed SS2 approach (Equation (10)) and the SS1 semi-empirical model (Equation (2)) described in Reference [5] is given in Figure 12 and in Table 6.          Table 6. Comparison of the results of numerical simulations using the SS1 model (Szusta Seweryn 2015) and SS2 model (Equation (10)) and experimental tests of the number of cycles N f for different loading paths (EN AW-2024 T3).

Load Path Calculations N f (cycle)
Exp. N f