Seam Weld Fatigue

The seam weld fatigue implemented in SimSolid is based on linearization of stress in the conjunction with Volvo method and following that it predicts fatigue damage and life. Linearized stress decomposes a through-thickness elastic stress field into equivalent membrane and bending.

This method determines the contribution of bending to the total stress, and from this determines whether the weld is essentially stiff or flexible. The method typically requires two S-N curves. One is a bending S-N curve which is dominated by bending stress, and the other is a membrane S-N curve which dominated by membrane stress. Interpolation is made between the bending and membrane SN curves based on the degree of bending.

Predicted or likely failure (fatigue crack) locations at the weld toe are marked. These are the locations where the fatigue damage will be evaluated.


Figure 1. Fillet weld cross-section showing likely failure locations

Weld Stress Calculations

The seam weld fatigue damage calculation in SimSolid uses a structural stress at the location of interest using the stress linearization method (link to SimSolid linearized stress documentation) and the bending ratio associated with that stress. Linearized stresses are obtained in a local coordinate system along the stress linearization segment. The local coordinate system is based on the start and end points of the segment as shown in Figure 2. The X-axis of the system is along the segment from entry to exit points. The other two axes are calculated as follows:
  • If the local X-axis is not parallel to the global Y-axis:(1)
    Zlocal=Xlocal*YglobalYlocal=Zlocal*Xlocal
  • If the local X-axis is parallel to the global Y-axis:
    Local Y-axis ( Ylocal ) is negative of global-X if local-X is along positive global-Y, and vice versa.(2)
    Zlocal=Xlocal*Ylocal


Figure 2.
From the extracted stress values above, the average membrane stress tensor plus the ending bending stress tensors at the entry and exit points are calculated using numerical integration.(3)
σmi=1TT/2T/2σidx
(4)
σbiS=6T2T/2T/2σixdx
(5)
σbiE=σbiS
Where,
  • σmi is equal to the ith component of membrane stress.
  • σi is equal to the ith component of extracted stress value.
  • σbiS is equal to the ith component of bending stress at the entry.
  • σbiE is equal to the ith component of bending stress at the exit.
  • L is equal to the ith Length of the stress linearization segment.
  • x is equal to the position of a point along the segment.
The stress quantity used for the seam weld damage parameter is the sum of the membrane and bending stresses.


Figure 3.
Hence, the stress at the top and the bottom surface is derived from the below equations:(6)
σtop=σmi+σbiS
(7)
σbottom=σmiσbiS(since σbiE=σbiS)
Note: Therefore, the calculation method used in the seam weld analysis engine is equally applicable to stresses calculated from either solid or thin-shell models. This generates top and bottom surface stresses in the same form as would be generated from a thin-shell model.


Figure 4.

Bending Ratio (r)

Experiments show that two types of SN curves are required to perform seam weld fatigue analysis based on a method suggested by M. Fermér, M Andréasson, and B Frodin. Based on lab tests, two SN curves were plotted (Figure 5). The upper curve is obtained in tests where the maximum stress is dominated by bending moment and the lower curve is obtained from tests where membrane force dominates the maximum stress.


Figure 5.
The upper and lower curves are referred to as bending S-N curve and membrane S-N curve, respectively. It is recommended that membrane S-N curve should be used when membrane stress dominates in an element and a bending S-N curve should be used when bending stress dominates. Interpolation between the two curves may be carried out depending on the degree of bending dominance.(8)
r=|σbi||σbi|+|σmi|
Where,
  • σbi is the maximum bending stress equal to 6T2T/2T/2σixdx .
  • σmi is the maximum membrane stress equal to 1TT/2T/2σidx
The average bending ratio, (rAVGb) , is defined as:(9)
(rAVGb)=ni=1(r[σ2TOP]ni=1[σ2TOP]
Where,
  • σ2TOP is the square of the maximum stress at the top surface when the damage is calculated, that is, at weld toe.
  • r is the bending ratio.

It is the weighted average of the bending ratio over all points in the loading time history.

An interpolation factor (IF) is now defined as:(10)
IF=0.0 when 0.0rAVGbrTHRb
(11)
IF=rAVGbrTHRb1rTHRb when rTHRb<rAVGb1.0
The rTHRb value is defined by the bending ratio threshold in Fatigue solution settings. It is set to 0.5 by default. If average bending ratio (rAVGb ) is less than or equal to the bending ratio threshold (rTHRb ), then the membrane S-N curve is used to assess damage. If average bending ratio is greater than the bending ratio threshold, then an S-N curve that is interpolated between membrane S-N curve and the bending S-N curve is used.

Interpolation between Membrane and Bending S-N Curves



Figure 6.
The linear interpolation method, as shown in Figure 6, uses the value of the interpolation factor (IF). For the interpolated curve, the calculation is done as follows, for the Fatigue Strength coefficient value (SRI1):(12)

SRI1interpolated = SRI1membrane + (SRI1bending – SRI1membrane) ∙ IF

This defines the stress level at 1 cycle.(13)

Nc1interpolated = 10(log10Nc1membrane + (log10 Nc1bending – log10(Nc1membrane) ∙IF))

(14)

S1interpolated = S1membrane + (S1bending – S1membrane) ∙ IF

This defines the stress level at Nc1interpolated cycles. These two points define the first section of the curve up to Nc1interpolated cycles. The last section is defined by finding a third point as follows. A life value is defined being 10 times greater of the Nc1 values for the stiff and flexible curves. From these, we can calculate S2bending and S2membrane. From these, we can interpolate to get S2interpolated which defines the high cycle part of the curve.(15)

S2interpolated = S2membrane + (S2bending – S2membrane). IF

Thickness

Optionally, a thickness (size effect correction) may be applied, based on thickness t of the part. It operates as follows:

There is no effect if tTref (the reference thickness or threshold can be specified in the fatigue solution settings).

The fatigue strength is reduced by a factor of (Treft)n , where n is the thickness exponent, if t>Tref , at all lifetimes (used as a factor to up the stress).

Mean Stress Correction

FKM mean stress correction is supported for seam weld fatigue. Stress sensitivity can be defined in the fatigue solution settings dialog via the mean stress correction field. Mean stress correction for seam weld fatigue can be enabled through the seam weld fatigue solution settings dialog.

Based on FKM-Guidelines, the Haigh diagram is divided into four regimes based on the Stress ratio (R=Smin/Smax) values. The corrected value is then used to choose the SN curve for the damage and life calculation stage.

The FKM equations below illustrate the calculation of Corrected Stress Amplitude ( SAe ). The actual value of stress used in the Damage calculations is the Corrected stress Amplitude (which is 2SAe ). These equations apply for SN curves that you input.
  • Regime 1 (R>1.0): SAe=Sa(1M)
  • Regime 2 ( R0.0 ): SAe=Sa+M*Sm
  • Regime 3 (0.0<R<0.5): SAe=(1+M)Sa+(M/3)*Sm1+M/3
  • Regime 4 ( R0.5 ): SAe=3Sa(1+M)23+M
Where,
  • SAe is the stress amplitude after mean stress correction (Endurance stress).
  • Sm is the stress amplitude, and M is the mean stress sensitivity.


Figure 7.