# Sine Sweep Fatigue Analysis

The study of fatigue life of structures under Sinusoidal Loading.

Frequency Response Stress results from the Frequency Response Analysis are used to calculate the frequency-dependent stress range, which are subsequently used to calculate Fatigue damage based on sweep rate.

## Input

Sine sweep is a well-known method of vibration testing which is conducted to evaluate the vibration behavior of a specimen.

- The Frequency Response subcase provides the input sinusoidal vibration profile that is required to conduct sine sweep fatigue analysis.
- The corresponding sine
sweep parameters, i.e. the
`SR`(sweep rate) and`SRUNIT`(sweep rate unit) fields, are specified with each fatigue event created. - The frequencies at which the sine sweep fatigue calculation are performed are specified in SN or EN uniaxial contexts.
- The static subcase can be added to the event in the load map for mean stress correction.
- The number of sweeps are to be specified alongside the event (default = 1 Sweep). The number of sweeps will be considered as time ‘T’ seconds if only one Frequency is considered, or if the SR = 0.0.

## Stress-Life

Uses the frequency-based stress range output from Frequency Response Analysis to compile a Damage versus Frequency curve.

The key for Sine Sweep Fatigue Analysis is to find the number of cycles that occur per frequency. Subsequently, the amount of Damage per cycle at a particular frequency can be calculated. This is multiplied by the number of cycles at that frequency to find the total damage per frequency. This is then integrated over the entire frequency range to find the damage per sweep. Fatigue damage for multiple sweeps is just a multiplication of the damage per sweep times the number of sweeps.

### Sweep Rate in Hertz Per Second (SRUNIT=HZPS)

Sweep rate in Hertz per second implies the following relationship between frequency and time, for a given sweep rate $\upsilon $ .

- $f$
- Frequency value of interest .
- ${f}_{0}$
- Initial frequency value of loading frequency range from Frequency Response subcase.
- $\upsilon $
- Sweep rate in Hertz per second (
`SR`field in the load map). - $t$
- Time at which the stress range $\text{\Delta}{\sigma}_{f}$ is applied with a frequency $f$ in sine-sweep fatigue.

The stress range versus frequency is output from the Frequency response analysis. This is used to identify the stress range for a particular frequency. To calculate the number of cycles of the stress range, $\text{\Delta}{\sigma}_{f}$ , at a particular frequency, the following formula is used.

- $f$
- Frequency value of interest .
- $\upsilon $
- Sweep rate in Hertz per second (
`SR`field in the load map).

You can see that the number of cycles of a stress range $\text{\Delta}{\sigma}_{f}$ at a particular frequency $f$ is equal to ${n}_{f}^{cyc}$ . This can be also defined as the rate of change of number of cycles with respect to frequency, at an arbitrary frequency.

The damage value per cycle ${D}^{cyc}\left(f\right)$ at frequency $f$ , is calculated based on the stress range $\text{\Delta}{\sigma}_{f}$ from the SN curve.

- ${D}^{cyc}\left(f\right)$
- Damage per cycle at frequency $f$
- ${D}_{SN}^{\text{\Delta}{\sigma}_{f}}$
- Damage value corresponding to stress range $\text{\Delta}{\sigma}_{f}$ from the SN curve.

Based on the damage per cycle at a frequency, ${D}^{cyc}\left(f\right)$ , a curve of Damage versus Frequency is constructed. This is constructed by first calculating the Damage per frequency:

Where, $D\left(f\right)$ is the damage per frequency.

Based on the Miner’s Damage accumulation rule, the total damage for a single sine sweep is calculated by integrating for the area under the Damage versus frequency curve.

Where, ${D}_{Sweep}$ is the total damage per sweep.

The total damage is calculated by multiplying the damage for each sweep by the number of sweeps.

Where, ${N}_{Sweep}$ is the number of sweeps.

If `SR`=0.0 (no sweep) or only one loading frequency is defined in the
frequency response subcase, then the damage value is calculated for that particular
frequency. In the case of `SR`=0.0, the first loading frequency in the
frequency response subcase is used. For the case of a single loading frequency input on
frequency response subcase, that particular frequency is chosen. In such cases, the total
damage is the damage per frequency multiplied by the total time `T#`
defined for the event.

### Sweep Rate in Octaves Per Minute (SRUNIT=OCTPM)

If the Sweep rate is defined in Octaves per minute, the sweep rate calculation involves a different equation compared to that of Hertz per second.

Sweep rate in Octaves per minute implies the following relationship between frequency and time, for a given sweep rate $\upsilon $ .

- $f$
- Frequency value of interest .
- ${f}_{0}$
- Initial frequency value of loading frequency range from Frequency Response subcase.
- $\upsilon $
- Sweep rate in Octvates per minute (
`SR`field in the load map). - $t$
- Time at which the stress range $\text{\Delta}{\sigma}_{f}$ is applied with a frequency $f$ in sine-sweep fatigue.

Applying natural logarithm function on both sides and taking derivative with respect to frequency, the equation reduces to:

The stress range versus frequency is output from the Frequency response analysis. This is used to identify the stress range for a particular frequency. To calculate the number of cycles of the stress range, $\text{\Delta}{\sigma}_{f}$ , at a particular frequency, the following formula is used.

- $f$
- Frequency value of interest .
- $\upsilon $
- Sweep rate in Octvates per minute (
`SR`field in the load map).

You can see that the number of cycles of a stress range $\text{\Delta}{\sigma}_{f}$ at a particular frequency $f$ is equal to ${n}_{f}^{cyc}$ (which can also be defined as the number of cycles per unit frequency, at an arbitrary frequency). Also, from the formula, for the case of Octaves per minute, the number of cycles of each frequency is independent of the frequency. Therefore, all frequencies have the same number of cycles.

The damage value per cycle ${D}^{cyc}\left(f\right)$ at frequency $f$ , is calculated based on the stress range $\text{\Delta}{\sigma}_{f}$ from the SN curve.

- ${D}^{cyc}\left(f\right)$
- Damage per cycle at frequency $f$
- ${D}_{SN}^{\text{\Delta}{\sigma}_{f}}$
- Damage value corresponding to stress range $\text{\Delta}{\sigma}_{f}$ from the SN curve.

Based on the damage per cycle at a frequency, ${D}^{cyc}\left(f\right)$ , a curve of Damage versus Frequency is constructed (Figure 2). This is constructed by first calculating the Damage per frequency.

Where, $D\left(f\right)$ is the damage per frequency.

Based on the Miner’s Damage accumulation rule, the total damage for a single sine sweep is calculated by integrating for the area under the Damage versus frequency curve.

Where, ${D}_{Sweep}$ is the total damage per sweep.

The total damage is calculated by multiplying the damage for each sweep by the number of sweeps.

Where, ${N}_{Sweep}$ is the number of sweeps.

If `SR`=0.0 (no sweep) or only one loading frequency is defined in the
frequency response subcase, then the damage value is calculated for that particular
frequency. In the case of `SR`=0.0, the first loading frequency in the
frequency response subcase is used. For the case of a single loading frequency input on
frequency response subcase, that particular frequency is chosen. In such cases, the total
damage is the damage per frequency multiplied by the total time `T#`
defined for the event.

## Strain-Life

Used for vibration based fatigue is calculated based on the damage caused by the highest stress-range.

The process for Strain-life fatigue is similar to that of stress-life, wherein, the characterization of number of cycles versus the stress-range is utilized to calculate the damage due to sine-sweep.

Since loading and unloading path information is lost in vibration analysis, the strain-life approach in vibration fatigue is a conservative extension of the stress-range based data.

Using this updated stress strain cycle data, the rest of the sine-sweep procedure continues in the same way as in stress-life sine sweep fatigue analysis.

The SWT or Morrow damage models are used to incorporate the effect of updated stress information (as a result of hanging procedure explained in the previous section). The updated mean stress for the new inner loops is used for Morrow and the updated maximum stress for the new inner loops is used for SWT damage models.

The damage value per cycle ${D}^{cyc}\left(f\right)$ at frequency $f$ , is calculated based on the strain values and updated mean/maximum stress range from the EN curve. For further information, refer to Mean Stress Correction.

- ${D}^{cyc}\left(f\right)$
- Damage per cycle at frequency $f$
- ${D}_{EN}^{\text{\Delta}{\sigma}_{f}}$
- Damage value corresponding to strain and updated mean stress from the EN curve.

Based on the damage per cycle at a frequency, ${D}^{cyc}\left(f\right)$ , a curve of Damage versus Frequency is constructed (Figure 2). This is constructed by first calculating the Damage per frequency.

- $D\left(f\right)$
- Damage per frequency.

Based on the Miner’s damage accumulation rule, the total damage for a single sine sweep is calculated by integrating for the area under the Damage versus Frequency curve.

- ${D}_{Sweep}$
- Total damage per sweep.

If `SR`=0.0 (no sweep) or only one loading frequency is defined in the
frequency response subcase, then the damage value is calculated for that particular
frequency. In the case of `SR`=0.0, the first loading frequency in the
frequency response subcase is used. For the case of a single loading frequency input on
frequency response subcase, that particular frequency is chosen. In such cases, the total
damage is the damage per frequency multiplied by the total time `T#`
defined for the event.

## Current Support of Sine Sweep Fatigue Analysis

- Only von Mises is supported for the stress combination.
- Only uniaxial fatigue analysis is supported.
- SN and EN for solid and shell elements are supported.
- Only one FRF subcase is allowed per sine sweep event.