Using the Fourier Methods for Cycle Counting of Bimodal Stress Histories with Variable in Time Amplitudes of Components

The application of the Fourier methods to perform cycle identification and description for different cases of bimodal stress histories are presented and discussed in the paper. The direct spectral method and the modified direct spectral method, previously co-proposed by the author, together with the aspects of the use of Fourier methods discussed in this article, offer a unique alternative to the methods known in the literature for identifying and counting various types of bimodal stress variations in time (with constant or variable amplitudes; completely reversed, reversal, or pulsed type). The application of the Fast Fourier Transform (FFT) and the Short-Time Fourier Transform (STFT) is discussed. The method is useful, especially in cases when it is known that due to the form of work of the mechanical system the two existing components of vibrations with various frequencies can be identified, of which the one with higher frequency has a vibration amplitude lower than a component with a lower frequency and, above all, a variable in time.


Introduction
The vibration of engineering structures often produces a few-modal type of stress history. The most common is the dominant one frequency history (one-modal case). All methods of fatigue cycle-counting are referred to this case. The second type of considered cases is the bimodal one. This type can be observed, for example, in a ship's hull when it is forced to vibrate by the propulsion system with one dominant frequency and by the impact of sea waves [1], or in the suspension of the car [2]. In previous research, the case of the tri-modal type has been separated in fatigue analysis [3]. An example of the application is a new tri-modal model under stochastic loads by C. Han, X. Qu, S. Ding, and Y. Ma [4].
This article relates to the case of bimodal vibrations generating a bimodal history of stress variability in structure. The literature presents the results of scientific research on the issues of fatigue that analyze the bimodal process. S. Sakai and H. Okamura [5] analyzed a case with two dominant frequencies, clearly spaced from each other. It was a generalization of previous work on the problem in which a narrow frequency band was associated with one frequency [6,7]. The disadvantage of this proposed approach was the lack of consideration of the increase in the stress amplitude of the primary cycle with low frequency due to the occurrence of a component with a higher frequency. Therefore, T.-T. Fu and D. Cebon proposed their own concept of the counting of stress cycles that was based on the idea of taking into account the superposition of the amplitudes of harmonic components with low and high frequencies. Based on the above-mentioned methods, D. Benasciutti and R. Tovo [8] proposed the so-called modified Fu-Cebon method. The other attempt dedicated to the case of the superposition of deterministic-excited vibrations with impulse-excited ones whilst considering the damping effect was proposed by G. Jiao and T. Moan [9]. C. Han, Y. Ma, X. Qu, and M. Yang [10] were focused on the above-mentioned three methods: Jiao-Moan, Fu-Cebo,n and modified Fu-Cebon. It has been pointed out that these bimodal methods are based on the knowledge of the frequency characteristics of stress and not the variation of stress in the time domain. Moreover, in the article [10] an analytical solution for predicting the vibration-fatigue-life in bimodal spectra is developed based on random attempts.
M.S. Kozień and D. Smolarski formulated their own method for the cycle-counting of bimodal stress of completely reversed type, based on their frequency characteristics. The method was formulated at the beginning for uniaxial stress [11,12] and later for multiaxial stress [13]. In the article [14], P.J. Romanowicz, D. Smolarski, and M.S. Kozień formulated and applied the modified direct spectral method. The idea was to apply the correction factor taking into account the influence of the mean compressive stresses. The method was presented with the example of a thrust roller bearing, using the results of numerical simulations of stress variation in time and based on analysis of the time domain. It was concluded that the proposed method (in contrast to multiaxial high-cycle fatigue criteria) included both kinds of stress waveforms in the fatigue analyses and could be applied for the identification of stress components and determination of fatigue life. It should be noted that, as far as the analysis presented in this article is concerned, the aim of which is to show the application of Fourier methods for identifying and describing cycles, the approaches to the direct spectral method [11][12][13] and the modified direct spectral method [14] presented in these articles are no different.
In engineering practice, the form of stress variability (stress tensor) is usually known. It can be ascertained directly from measurements carried out on the real object (for the structure or sample element), from the analytical solution of the problem, or as the result of computer simulations carried out based on the structure calculation model used. To identify the amplitude-frequency characteristics, it is common to use Fourier analysis methods of signal, in particular Fourier Series (FS) and Fast Fourier Transform (FFT) for stationary signals and Short Time Fourier Transform (STFT) for non-stationary ones. This article discusses the problem of using the FFT and STFT functions to properly identify the amplitude values of bimodal components, considering on the one hand the stress functions with constant or variable amplitudes, and on the other hand the components of the completely reversed type, reversal, as well as the pulsed type. While these types of courses are often found in practice, in particular the analysis of the variable amplitude course is not discussed in literature. In these cases, it is suggested to use one of the methods of counting irregular stress cycles, most often the rain-flow method. In the case of the bimodal cases, these methods seem to be unnatural. Another possible approach is the application of dedicated spectral methods. However, also in this case, the variability of the amplitude of stress in time is difficult to identify, especially considering that in the fatigue analysis, the averaged amplitude of stress variability is not enough to analyze. The key role is played by a maximal range of variability, and the identification of such cycles should be focused upon above all. The direct spectral method and the modified direct spectral method, previously co-proposed by the author, together with the aspects of the use of Fourier methods discussed in this article, provides a unique alternative to the methods known in the literature for identifying and counting various types of bimodal stress variations in time (with constant or variable amplitudes; completely reversed, reversal, or pulsed type).
In the fatigue life analysis of the structure, the identified stress cycles of various types, usually characterized by their number, average stress value, and stress amplitude value, are data for further calculations. Different methods of cumulative damage are used for this purpose. Commonly used in engineering applications is the Palmgren-Miner linear cumulative damage rule [15,16]. The method is simple in its principle but has its limitations in application. Therefore, different non-linear approaches are currently used, e.g., the double linear damage rule; Marco and Starkey; Subramanyan, Hashin, and Rotem; Corten and Dolon; Freudenthal-Heller; Serensen; Subramanayan; Bui-Quoc [17]; and Liou [7]. This aspect of fatigue analysis is not the purpose of this article.
The research question posed in this article is as follows: to what extent is it possible and effective to use Fourier's methods of function analysis, in combination with the direct and modified direct spectral method, to count stress cycles with variable amplitudes in time, of the bimodal type, whose courses are irregular in time.
The paper is divided into five sections. A state of the art introduction to the bimodal stress fatigue problem is given in Section 1. The theoretical background and interpretation of the application of the Fourier methods (FS, FFT, STFT) is discussed in Section 2. The theoretical background of direct spectral methods is described in Section 3. Cycle identification using FFT and STFT methods for three specially selected examples of stress variation over time, involving a complex combination of types of variation (reversal, pulsed) and variation of amplitude of the higher frequency component is thoroughly discussed in Section 4. The two examples of stress variation are specially generated for analysis, and the third refers to the stresses generated in a ball bearing. The results obtained using the rainflow cycle counting method, commonly used for irregular stress histories, and the direct spectral method, using FFT or STFT for identification process, were compared. Finally, the paper is concluded in Section 5.

Identification of Pulsating Waveform Based on Fourier Analysis
The commonly used method for identifying the frequency characteristic of the function is the so-called Fourier analysis. For stationary functions (signals in Theory of signals), it contains Fourier Series (FS) and Fourier Transform (FT). In computer realization commonly used in algorithm of FT is the Fast Fourier Transform (FFT). For non-stationary ones, the Short Time Fourier Transform (STFT) is used. Let us define function x(t) in the range t∈[0,T], assuming the fulfillment of the Dirichlet conditions. It can be shown in the form of infinite series (1), called as Fourier Series (FS) representation of x(t), where: Each function x(t), for which the integral of its absolute value over the real space R is finite, can be transformed into the frequency value using the Fourier Transform (2), where I is the imaginary unit. The FT Fourier Transform (FT) is defined in the complex space. Usually, the Fast Fourier Transform (FFT) algorithm is used in this case.
When the function is non-stationary, the generalization of the FT (2) is attained by inserting the specially defined window function g(t) centered for time τ running along time domain (independent time variable), defining the STFT function of x(t) (3).
The FS and especially FT in FFT realization is commonly used in identifying the frequency form of the signals in the dynamic analysis of structures for stationary functions and STFT for non-stationary ones. In particular, they can be applied to the identification or verification of the bimodal form of the signal.
In fatigue analysis, the most commonly analyzed stress variation in time domain has the form of a completely reversed form. Sometimes, however, the pulsating form is recognized. It should be noted that the spectral forms of the two types of realization have completely different forms in Fourier analyses. Therefore, it is very important to properly interpret the Fourier spectras of stress signals, because the direct spectral method is based on reconstruction of the two harmonic components based on the spectrum of the stress signal. Let us consider the two test signals-the harmonic sine one modeling the completely reversed stress (4) and the signal in the form of absolute value of harmonic sine modelling the pulsating (zero-to-tension or zero-to compression) stress (5).
x(t) = A|sin(ωt)| The time-history of the harmonic function (4) and module of its FFT function is shown in Figure 1. The second pair of graphs depicting two cases of analysis of the signal (5) is shown in Figures 2 and 3. In Figure 3 the signal in the form (5) shifted vertically after the procedure subtract mean was applied. The subtract mean procedure consisted of shifting the graph of the function (function course) in the vertical direction (change of function value) in such a way that, finally, the average value of the function in the analyzed interval was equal to zero. interpret the Fourier spectras of stress signals, because the direct spectral method is based on reconstruction of the two harmonic components based on the spectrum of the stress signal. Let us consider the two test signals-the harmonic sine one modeling the completely reversed stress (4) and the signal in the form of absolute value of harmonic sine modelling the pulsating (zero-to-tension or zero-to compression) stress (5).
The time-history of the harmonic function (4) and module of its FFT function is shown in Figure 1. The second pair of graphs depicting two cases of analysis of the signal (5) is shown in Figures 2 and 3. In Figure 3 the signal in the form (5) shifted vertically after the procedure subtract mean was applied. The subtract mean procedure consisted of shifting the graph of the function (function course) in the vertical direction (change of function value) in such a way that, finally, the average value of the function in the analyzed interval was equal to zero.    interpret the Fourier spectras of stress signals, because the direct spectral method is based on reconstruction of the two harmonic components based on the spectrum of the stress signal. Let us consider the two test signals-the harmonic sine one modeling the completely reversed stress (4) and the signal in the form of absolute value of harmonic sine modelling the pulsating (zero-to-tension or zero-to compression) stress (5).
The time-history of the harmonic function (4) and module of its FFT function is shown in Figure 1. The second pair of graphs depicting two cases of analysis of the signal (5) is shown in Figures 2 and 3. In Figure 3 the signal in the form (5) shifted vertically after the procedure subtract mean was applied. The subtract mean procedure consisted of shifting the graph of the function (function course) in the vertical direction (change of function value) in such a way that, finally, the average value of the function in the analyzed interval was equal to zero.   on reconstruction of the two harmonic components based on the spectrum of the stress signal. Let us consider the two test signals-the harmonic sine one modeling the completely reversed stress (4) and the signal in the form of absolute value of harmonic sine modelling the pulsating (zero-to-tension or zero-to compression) stress (5).
The time-history of the harmonic function (4) and module of its FFT function is shown in Figure 1. The second pair of graphs depicting two cases of analysis of the signal (5) is shown in Figures 2 and 3. In Figure 3 the signal in the form (5) shifted vertically after the procedure subtract mean was applied. The subtract mean procedure consisted of shifting the graph of the function (function course) in the vertical direction (change of function value) in such a way that, finally, the average value of the function in the analyzed interval was equal to zero.    The following remarks relate to the corresponding forms of the FS representations: • For the signal (4) the values of coefficients are: a 0 = 0, a 1 = 1, a n = 0, n = 2, . . . , +∞, b n = 0, n = 1, . . . , +∞.
• For the signal (5) the values of coefficients are: The most important differences in the FFT spectras of signals (4) and (5) are: • The harmonic function (4) has only the single strip for the given frequency of the signal for the function; • The absolute value of harmonic function (5) has a series of strips with decreasing amplitudes with non-zero value for zero-frequency and multiplied (especially doubled) basic frequency.

Formulation of the Direct Spectral Method of Cycle Counting
The direct spectral method and the modified direct spectral method are dedicated to cases where the variability of bimodal stresses is the source of structure vibrations. The methods are based on the following assumptions:

•
The primary cycle is of a low-frequency character with a frequency ω 1 and a larger constant amplitude A 1 than the higher-frequency component with a frequency ω 2 and the time-varying amplitude A 2 (t).

•
In general, the frequency ω 2 is much greater than the frequency ω 1 and the amplitude A 1 is much greater than the maximum value of the amplitude A 2 (t).

•
The analysis is performed for the stress block period T B , calculated based on the values of T 1 and T 2 , where T B is the smallest total multiple of the T 1 period for which the ratio T B /T 1 is the approximate integer. • Generally, the variation of stress for low-frequency cycles has harmonic shape; if realistic forms is other (identified in computer simulation or measurement) the dedicated analytical function can be used for an approximation of this shape. • High-frequency components varying with frequency ω 2 generates the so-called secondary cycles.

•
The difference between the direct spectral method and the modified direct spectral method is the formula (the applied theory) that is used to consider the average stress value and determine finally the equivalent completely reversed stress. The direct spectral method was applied in the analyses discussed in [11][12][13] and the modified direct spectral method is presented in [14].

•
The difference in the use of direct spectral methods [11][12][13][14] and the spectral methods known in the literature [1][2][3][4][5][6][7][8][9][10], especially [7], lies in the description. The formulation of the direct spectral methods is deterministic, and the spectral methods are random. It seems, therefore, that direct methods are more natural in the analysis of bimodal waveforms. However, the key here is to identify the amplitudes and average stress values of the main cycle and secondary cycles. Fourier methods (FFT, STFT) can be used for this purpose.

•
The algorithm of direct spectral methods was successively developed and presented in articles [11][12][13][14]. In particular, article [14] contains a detailed algorithm for the last modified direst spectral method.

General Remarks
To test the applicability and limitations of the Fourier analysis for quantitative identification of the parameters of the identified stress cycles, an analysis was performed for three cases of bimodal stress cycles (time signals). These three cases are as follows: • Analytical (generated based on the analytical formula). Superposition of the absolute value of sine with constant amplitude for frequency ω 1 (pulsating form, zero-totension) and sine with constant amplitude for frequency ω 2 (reversal stress).
• Analytical (generated based on the analytical formula). Superposition of the absolute value of sine with constant amplitude for frequency ω 1 (pulsating form, zero-totension) and sine with amplitude linearly varying in time for frequency ω 2 . For the first half of period T 1 the amplitude grew to maximum value (A 2 ) max = A 2 (T 1 /2), and for the second half of period T 1, it fell to zero. • Numerical simulation based on solution of the realistic problem of stresses existing in rolling bearings due to repeated contact between rollers and rings. The signal was one of the components of stress tensor analysed and discussed by P.J. Romanowicz, D. Smolarski, and M.S. Kozień [14]. Due to the nature of the contact, the component with frequency ω 1 had a zero-to-tension form, but the variation in time did not have a sine form, and had to be approximated in a suitable form. Moreover, the component with frequency ω 2 ω 2 had a sine form with various in time amplitude that grew from zero to maximum value in the middle of the period T 1 T 1 , before falling to zero. The variability of amplitude A 2 in time depended on the rolling process of the ball in contact with the bearing raceway and was not defined in analytical form.
For all cases, the values of the average stress and stress amplitude were compared with those obtained using the rainflow method [18], commonly used in fatigue analyses. The VibrationData Toolbox ver. 13.3 by Tom Irvine in Partnership with ENDAQ was used for the analysis, in which the procedure was carried out in accordance with the standard ASTM E 1049-85. The values of the stress amplitude and mean stress were determined using the Fourier methods for the identified stress cycles and then the direct spectral method (called FOURIER in Tables); those obtained using the rainflow method (called RAINFLOW in Tables) were then compared. For the rainflow method, the following values were given: the average (Ave.) and maximum (Max.) value of the stress amplitude, as well as the minimum (Min.) and maximum (Max.) values of the mean stress value.

Analytical Simulation-Stationary Subcycles
The first example was devoted to the stress variation defined in the analytical form (8). The primary cycle had the pulsating form of absolute value of sine with amplitude of 800 MPa and frequency of 0.5 Hz (period equal to T 1 = 2 s). Due to the action of the absolute value type function, the final frequency of the signal was equal to 1 Hz with period equal to T B = 1 s (see Figure 4-left). The secondary cycles were reversal stress with amplitude of 200 MPa and frequency of 12.5 Hz. The visualization of the signal σ(t), named as sigma(t) on the plot, is shown in Figure 4 (left). The graphical form of FFT function of the function (6), named as SIGMA(f) on the plot, is shown in Figure 4 (right).
value of sine with constant amplitude for frequency ω1 (pulsating form, zero-to-tension) and sine with amplitude linearly varying in time for frequency ω2. For the first half of period T1 the amplitude grew to maximum value (A2)max = A2(T1/2), and for the second half of period T1, it fell to zero.  Numerical simulation based on solution of the realistic problem of stresses existing in rolling bearings due to repeated contact between rollers and rings. The signal was one of the components of stress tensor analysed and discussed by P.J. Romanowicz, D. Smolarski, and M.S. Kozień [14]. Due to the nature of the contact, the component with frequency ω1 had a zero-to-tension form, but the variation in time did not have a sine form, and had to be approximated in a suitable form. Moreover, the component with frequency ω2 ω2 had a sine form with various in time amplitude that grew from zero to maximum value in the middle of the period T1 T1, before falling to zero. The variability of amplitude A2 in time depended on the rolling process of the ball in contact with the bearing raceway and was not defined in analytical form.
For all cases, the values of the average stress and stress amplitude were compared with those obtained using the rainflow method [18], commonly used in fatigue analyses. The VibrationData Toolbox ver. 13.3 by Tom Irvine in Partnership with ENDAQ was used for the analysis, in which the procedure was carried out in accordance with the standard ASTM E 1049-85. The values of the stress amplitude and mean stress were determined using the Fourier methods for the identified stress cycles and then the direct spectral method (called FOURIER in Tables); those obtained using the rainflow method (called RAINFLOW in Tables) were then compared. For the rainflow method, the following values were given: the average (Ave.) and maximum (Max.) value of the stress amplitude, as well as the minimum (Min.) and maximum (Max.) values of the mean stress value.

Analytical Simulation-Stationary Subcycles
The first example was devoted to the stress variation defined in the analytical form (8). The primary cycle had the pulsating form of absolute value of sine with amplitude of 800 MPa and frequency of 0.5 Hz (period equal to T1 = 2 s). Due to the action of the absolute value type function, the final frequency of the signal was equal to 1 Hz with period equal to TB = 1 s (see Figure 4-left). The secondary cycles were reversal stress with amplitude of 200 MPa and frequency of 12.5 Hz. The visualization of the signal (t), named as sigma(t) on the plot, is shown in Figure 4 (left). The graphical form of FFT function of the function (6), named as SIGMA(f) on the plot, is shown in Figure 4    Identified cycles based on the applied spectral method are summarized in Table 1 (RAINFLOW denotes results obtained using the rainflow method and FOURIER those obtained by application of the Fourier transform and the direct spectral method). It should Identified cycles based on the applied spectral method are summarized in Table 1 (RAINFLOW denotes results obtained using the rainflow method and FOURIER those obtained by application of the Fourier transform and the direct spectral method). It should be emphasized that in the analyzed example, the stress block period was exactly equal to T B = 2 s. Comparison of different identified parameters of mean values and amplitudes of cycles, based on the rainflow method and application of the direct spectral method for stress variation defined by (6) is summarized in Table 1.   Figure 2 (right). The frequency of the second component was correctly identified as 12.5 Hz. However, the amplitude value read was 199 MPa and was very close to the value resulting from the analytical form of a function of 200 MPa. In this case, the value of the amplitude resulting from the FFT analysis was not multiplied by a scaling factor, because its source was a sine wave. This interpretation also resulted from the discussion presented in Section 2 referring to the FFT function for the sine wave and can be identified in Figure 1 (right). From the analysis of the case under consideration, it can be seen that for the correct identification of amplitude-frequency, it is necessary to reference quality to the course of the function σ(t) in the time domain.
Both methods identify approximately the same number of cycles (rainflow-26.5, Fourier-27). The Fourier method gave slightly higher stress amplitudes for the secondary cycles. On the other hand, the rainflow method additionally identified a stress half-cycle with a relatively high amplitude, which could not be identified on the time course. This approach is on the safe side and is conservative. On the other hand, the stress amplitudes identified for the two twin primary cycles had slightly higher values in the rainflow method. These results were not justified considering the time course of the function.
The identified values seem to be at a good level of agreement. The identification accuracy can also be increased by using the STFT method.

Analytical Simulation-Non-Stationary Subcycles
The second case was devoted to the stress variation whose envelope was, as in the previous example, the absolute value of sine with amplitude of 800 MPa and frequency of 0.5 Hz. The secondary cycles were reverse stress of sinus type with linearly variable amplitude starting with zero values for zero time, a maximal value of 200 MPa for midperiod time, zero value at the end of the period, and frequency of 12.5 Hz. Analytical form of the signal was defined by formula (7) where the weighting function was given by (8). The visualization of the function σ(t), named as sigma(t) is shown in Figure 5 (left).
EER REVIEW 9 of 14  Identified cycles based on the applied spectral method are summarized in Table 2  Because the second component of the signal with frequency ω 2 had amplitude varying in time, it was reflected in the transformation of the FFT function (see Figure 5-right). The identified frequency value (1 Hz) and amplitude 399 MPa (corresponding in reality to 799 MPa) was the same as in the case under consideration in Section 4.2. This component had constant amplitude. On the other hand, the identified frequency for the second component was located around 12.5 Hz, but the identified value of the amplitude was about 100 MPa and so was lower than the maximum value of 200 MPa. This is because in the entire analyzed period, the amplitude of this ingredient changed over time, and the obtained appointment in FFT analysis was a certain averaging in the entire period of the course analysis. This means that, for this case, the FFT function could not be used to identify the amplitude of the variability of stress amplitudes, especially given the fact that in fatigue analysis the identification of the range of the highest variable stress was of key importance. In this case, with the supported identification of the amplitude of the second component as a function of time, it was necessary to use the STFT function. The graphical form of STFT function is shown in Figure 6.  Identified cycles based on the applied spectral method are summariz (RAINFLOW denotes results obtained using the rainflow method and FO obtained by application of the Fourier transform and the direct spectral stress block period was exactly equal to TB = 2 s, as in the case analysed Comparison of different identified parameters of mean values and amplitu based on the rainflow method and application of the direct spectral method iation defined by (7), are summarized in Table 2. The number of identifi slightly higher for the Fourier method (by four cycles), but the amplitude v cycles were small and insignificant in fatigue analyses for a sample made o steel. The identified stress amplitude values for primary cycles were very The stress amplitude values were generally at the same level for identifi cycles with smaller amplitudes, the Fourier method provided more conserv   Identified cycles based on the applied spectral method are summarized in Table 3 (RAINFLOW denotes results obtained using the rainflow method and FOURIER those obtained by application of the Fourier transform and the direct spectral method). The stress block period was exactly equal to T B = 2 s, as in the case analysed in Section 4.2. Comparison of different identified parameters of mean values and amplitudes of cycles, based on the rainflow method and application of the direct spectral method for stress variation defined by (7), are summarized in Table 3. The number of identified cycles was slightly higher for the Fourier method (by four cycles), but the amplitude values of these cycles were small and insignificant in fatigue analyses for a sample made of aluminum or steel. The identified stress amplitude values for primary cycles were very close in value. The stress amplitude values were generally at the same level for identified cycles. For cycles with smaller amplitudes, the Fourier method provided more conservative results. Table 2. Comparison of identified mean values and amplitudes of cycles based on rainflow and Fourier methods for stress variation defined by (9).

Cycle Type
No.  The identified values seem to be at a good level of agreement. Accuracy of the identification can also be increased by using the STFT method with different parameters. This is possible by using a shorter time discretization of the time series of the stress course, which is built for the needs of discrete STFT analysis.

Stress Component of Rolling Bearing
The third example to show the possibility of using Fourier methods to identify irregular (with varying amplitudes in time) stress cycles refers to the stresses generated in the bearing balls during their operation. The issue of estimation of emerging stresses was discussed in articles [19,20]. However, in article [14] the possibility of using the modified spectral method for this purpose was shown, and its advantages over other commonly used identification methods were also shown, in particular considering compressive stresses. The analyses concerning the application of the direct method presented in article [14] referred to the comparison of methods and were based on the analysis of stress courses over time, obtained based on computer simulations. In the presented analysis concerning the use of Fourier methods in the implementation of the direct spectral method, the time course of the most difficult in the analysis case of the variability of the stress tensor component σ z in the bearing balls was used.    For the analyzed case of the stress component σz caused by the press on the bearing raceway, the stresses are pulsating type in nature, but the sha with frequency 1 was not sine-shaped (see Figure 7-left). This shape w exponent-like function. Therefore, for better identification of the mean str ary cycles in the direct spectral method, the exponential function was chos imation.
For the analyzed case TB = T1 = 0.8 s, the analytical form identified by t tral method stress variation in time is given by Formula (9). The formula is range of time ∈ [0, ], a Concat function denotes the concatenation of t nential functions along time and d is a parameter of the exponential functio proper shape of this envelope function (for the analysed case A1 = 800 MPa ( ) Identified cycles based on the applied spectral method are summari Comparison of identified mean values and amplitudes of cycles, based o method and the application of the Fourier method is summarized in Table   Table 3  For the analyzed case of the stress component σ z caused by the pressure of the ball on the bearing raceway, the stresses are pulsating type in nature, but the shape of variation with frequency ω1 was not sine-shaped (see Figure 7-left). This shape was more like an exponent-like function. Therefore, for better identification of the mean stress for secondary cycles in the direct spectral method, the exponential function was chosen for approximation.
For the analyzed case T B = T 1 = 0.8 s, the analytical form identified by the direct spectral method stress variation in time is given by Formula (9). The formula is defined in the range of time t ∈ [0, T 1 ], a Concat function denotes the concatenation of two used exponential functions along time and d is a parameter of the exponential function to define the proper shape of this envelope function (for the analysed case A 1 = 800 MPa, d = 0.06).
Identified cycles based on the applied spectral method are summarized in Table 4. Comparison of identified mean values and amplitudes of cycles, based on the rainflow method and the application of the Fourier method is summarized in Table 4.
The number of identified stress cycles using both methods was comparable, especially bearing in mind the fact that starting from about the 20th cycles, there were small stress amplitudes for aluminum or steel.
The values of the identified secondary stress cycles had higher values than the corresponding cycles identified by the rainflow method. The results obtained in this case using the Fourier method, combined with the direct spectral method, gave more reliable results. The value of the identified stress amplitude for the primary cycle for the rainflow method was clearly too high, bearing in mind the zero-to-compression (pulsating) nature of the course. This value obtained by the Fourier method, together with the use of the direct spectral method, was correctly identified.
Overall, the identification of cycles using the Fourier method, combined with the direct spectral method, gave better results than those obtained using the rainflow method.

Conclusions
The application of the Fourier methods for cycle identification and description for different cases of bimodal irregular stress histories was discussed in the paper. Based on this study, the following conclusions can be formulated:

1.
The direct spectral method and the modified direct spectral method, previously coproposed by the author, together with the aspects of the use of Fourier methods discussed in this article, provides a unique alternative to the methods known in the literature for identifying and counting various types of bimodal stress variation in time (with constant or variable amplitudes; completely reversed, reversal, or pulsed type).

2.
The direct spectral method [10][11][12][13] and the modified direct spectral method [14] method are useful, especially in cases when it is known that due to the form of work of the mechanical system, the two existing components of vibrations with various frequencies can be identified, of which the one with the higher frequency has a vibration amplitude lower than a component with a lower frequency and, above all, a variable in time. In these cases, the application of the spectral methods or the rainflow method for cycle-counting seems to be unnatural.

3.
The Fourier-based identification method, applied together with the direct spectral method for bimodal waveforms, is a more natural method. For irregular waveforms with variable amplitude, it can give better results than other methods. In this case, it is worth using this method at least as a checking method in relation to other methods.

4.
Frequency analysis is often carried out using Fourier methods. The article thoroughly discusses the way the results are used and interpreted using the FFT and STFT methods. For frequency analysis of the stress with variable in time amplitude of the component, the Wavelet Transform (WT) can be applied too, as an alternative approach. 5.
When using Fourier methods in a discreet formulation, as is the case in computer analyses, the proper selection of the parameters of the analysis is crucial. For FFT analyses, they are sampling time ∆t and associating it with discretization in frequency ∆f. For STFT analyses, they are sampling time ∆t, number of segments, number of samples per segment, time duration of the segment, and discretization in frequency ∆f.