Stress Gradient Effect

Stress gradient effect can be taken into consideration through either FKM guideline method or Critical Distance method.

It is supported for both shells and solid elements. For solid elements, the stress gradient effect is only available with grid point stress in fatigue analysis using results of static analysis. For solid elements, SURFSTS field on FATPARM is automatically set to GP when Stress Gradient effect is activated.

The Stress Gradient method is supported for Uniaxial and Multiaxial SN, EN and FOS Fatigue. It is not supported for Weld, Vibration, and Transient Fatigue analyses.

FKM Guideline Method

In the FKM guideline method, stress gradient effect is considered by increasing fatigue strength by a factor calculated using a rule in FKM guidelines. In OptiStruct implementation of FKM guideline method, 6 components of a stress tensor at each time step is reduced by the factor provided by FKM guidelines.

To activate Stress Gradient effect using FKM guideline method, the GRD field on FATPARM should be set to GRDFKM.

The following steps are followed to reduce stresses at the surface to take stress gradient effect into consideration.
  1. Calculate stress gradient of 6 components of a stress tensor, Δ σ ij ( t ) Δz MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaSaaaeaacq qHuoarcqaHdpWCdaWgaaWcbaGaamyAaiaadQgaaeqaaOWaaeWaaeaa caWG0baacaGLOaGaayzkaaaabaGaeuiLdqKaamOEaaaaaaa@402A@ , at each time step after linear combination of stress history. z-direction is an outward surface normal. For a solid element, the gradient is calculated by finite difference between stress at surface and stress at 1mm below the surface. The stress at 1mm below surface is an interpolated stress from grid point stresses of an element of interest. In case of 2nd order solid elements, only grid point stresses at corners are used for interpolation. For shell elements, the gradient is calculated from stresses of both layers and its thickness.
  2. Using the stress gradient obtained in Step 1, a gradient of equivalent stress in the surface normal direction, Δ σ eq ( t ) Δz MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaSaaaeaacq qHuoarcqaHdpWCdaWgaaWcbaGaamyzaiaadghaaeqaaOWaaeWaaeaa caWG0baacaGLOaGaayzkaaaabaGaeuiLdqKaamOEaaaaaaa@402D@ , is calculated in an analytical way at each time step. The equivalent stress can be either von Mises stress or absolute maximum principal stress.
  3. The related stress gradient, G σ MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4ramaaBa aaleaacqaHdpWCaeqaaaaa@38B2@ is calculated using the following normalization.(1)
    G ¯ ( t ) σ = 1 σ e q ( t ) Δ σ e q ( t ) Δ z MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGabm4rayaara WaaeWaaeaacaWG0baacaGLOaGaayzkaaWaaSbaaSqaaiabeo8aZbqa baGccqGH9aqpdaWcaaqaaiaaigdaaeaacqaHdpWCdaWgaaWcbaGaam yzaiaadghaaeqaaOWaaeWaaeaacaWG0baacaGLOaGaayzkaaaaamaa laaabaGaeuiLdqKaeq4Wdm3aaSbaaSqaaiaadwgacaWGXbaabeaakm aabmaabaGaamiDaaGaayjkaiaawMcaaaqaaiabfs5aejaadQhaaaaa aa@4DB8@
  4. Calculate the correction factor n σ ( t ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOBamaaBa aaleaacqaHdpWCaeqaaOWaaeWaaeaacaWG0baacaGLOaGaayzkaaaa aa@3B65@ . Refer to Correction Factor Calculation.
  5. Apply the correction factor n σ MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOBamaaBa aaleaacqaHdpWCaeqaaaaa@38D9@ to the surface stress tensor to obtain reduced surface stress. Apply the same n σ MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOBamaaBa aaleaacqaHdpWCaeqaaaaa@38D9@ to corresponding strain tensor to obtain reduced strain tensor when EN fatigue analysis is to be carried out with nonlinear analysis.(2)
    σ ' ij ( t )= σ ij ( t ) n σ ( t ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeq4WdmNaai 4jamaaBaaaleaacaWGPbGaamOAaaqabaGcdaqadaqaaiaadshaaiaa wIcacaGLPaaacqGH9aqpdaWcaaqaaiabeo8aZnaaBaaaleaacaWGPb GaamOAaaqabaGcdaqadaqaaiaadshaaiaawIcacaGLPaaaaeaacaWG UbWaaSbaaSqaaiabeo8aZbqabaGcdaqadaqaaiaadshaaiaawIcaca GLPaaaaaaaaa@49D6@

Correction Factor Calculation

Correction factor calculation is based on relationship between n σ MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOBamaaBa aaleaacqaHdpWCaeqaaaaa@38D9@ and G σ MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4ramaaBa aaleaacqaHdpWCaeqaaaaa@38B2@ described in the FKM guidelines.

According to FKM guidelines, the stress correction factor n σ MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOBamaaBa aaleaacqaHdpWCaeqaaaaa@38D9@ is determined by:
  • If  G ¯ σ 0.1  mm 1 n σ = 1 + G ¯ σ m m 10 ( a G 0.5 + R m b G MPa ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGceaqabeaacaqGjb GaaeOzaiaabccaceWGhbGbaebadaWgaaWcbaGaeq4Wdmhabeaakiab gsMiJkaaicdacaGGUaGaaGymaiaabccacaqGTbGaaeyBamaaCaaale qabaGaeyOeI0IaaGymaaaaaOqaaiaad6gadaWgaaWcbaGaeq4Wdmha beaakiabg2da9iaaigdacqGHRaWkceWGhbGbaebadaWgaaWcbaGaeq 4WdmhabeaakiabgwSixlaad2gacaWGTbGaeyyXICTaaGymaiaaicda daahaaWcbeqaaiabgkHiTmaabmaabaGaamyyamaaBaaameaacaWGhb aabeaaliabgkHiTiaaicdacaGGUaGaaGynaiabgUcaRmaalaaabaGa amOuamaaBaaameaacaWGTbaabeaaaSqaaiaadkgadaWgaaadbaGaam 4raaqabaWccqGHflY1caqGnbGaaeiuaiaabggaaaaacaGLOaGaayzk aaaaaaaaaa@651E@
  • If   0.1  mm 1 < G ¯ σ 1  mm 1 n σ = 1 + G ¯ σ mm 10 ( a G + R m b G MPa ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGceaqabeaacaqGjb GaaeOzaiaabccacaqGGaGaaGimaiaac6cacaaIXaGaaeiiaiaab2ga caqGTbWaaWbaaSqabeaacqGHsislcaaIXaaaaOGaeyipaWJabm4ray aaraWaaSbaaSqaaiabeo8aZbqabaGccqGHKjYOcaaIXaGaaeiiaiaa b2gacaqGTbWaaWbaaSqabeaacqGHsislcaaIXaaaaaGcbaGaamOBam aaBaaaleaacqaHdpWCaeqaaOGaeyypa0JaaGymaiabgUcaRmaakaaa baGabm4rayaaraWaaSbaaSqaaiabeo8aZbqabaGccqGHflY1caqGTb GaaeyBaaWcbeaakiabgwSixlaaigdacaaIWaWaaWbaaSqabeaacqGH sisldaqadaqaaiaadggadaWgaaadbaGaam4raaqabaWccqGHRaWkda WcaaqaaiaadkfadaWgaaadbaGaamyBaaqabaaaleaacaWGIbWaaSba aWqaaiaadEeaaeqaaSGaeyyXICTaaeytaiaabcfacaqGHbaaaaGaay jkaiaawMcaaaaaaaaa@68EB@
  • If   1  mm 1 < G ¯ σ 100  mm 1 n σ = 1 + G ¯ σ mm 4 10 ( a G + R m b G MPa ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGceaqabeaacaqGjb GaaeOzaiaabccacaqGGaGaaGymaiaabccacaqGTbGaaeyBamaaCaaa leqabaGaeyOeI0IaaGymaaaakiabgYda8iqadEeagaqeamaaBaaale aacqaHdpWCaeqaaOGaeyizImQaaGymaiaaicdacaaIWaGaaeiiaiaa b2gacaqGTbWaaWbaaSqabeaacqGHsislcaaIXaaaaaGcbaGaamOBam aaBaaaleaacqaHdpWCaeqaaOGaeyypa0JaaGymaiabgUcaRmaakeaa baGabm4rayaaraWaaSbaaSqaaiabeo8aZbqabaGccqGHflY1caqGTb GaaeyBaaWcbaGaaGinaaaakiabgwSixlaaigdacaaIWaWaaWbaaSqa beaacqGHsisldaqadaqaaiaadggadaWgaaadbaGaam4raaqabaWccq GHRaWkdaWcaaqaaiaadkfadaWgaaadbaGaamyBaaqabaaaleaacaWG IbWaaSbaaWqaaiaadEeaaeqaaSGaeyyXICTaaeytaiaabcfacaqGHb aaaaGaayjkaiaawMcaaaaaaaaa@69B1@
Table 1. Example values for Constants a G MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyyamaaBa aaleaacaWGhbaabeaaaaa@37D5@ and b G MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyyamaaBa aaleaacaWGhbaabeaaaaa@37D5@
Constants Stainless Steel Other steels GS GGG GT GG Wrought Al-Alloys Cast Al- Alloys
a G MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyyamaaBa aaleaacaWGhbaabeaaaaa@37D5@ 0.40 0.50 0.25 0.05 -0.05 -0.05 0.05 -0.05
b G MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyyamaaBa aaleaacaWGhbaabeaaaaa@37D5@ 2400 2700 2000 3200 3200 3200 850 3200
Where,
GS
Cast Steel and Heat Treatable cast steel for general purposes.
GGG
Nodular Cast Iron.
GT
Malleable Cast Iron.
GG
Cast Iron with lamellar graphite (grey cast iron).

R m MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOuamaaBa aaleaacaWGTbaabeaaaaa@37EC@ is UTS in MPa and dimension of G σ MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4ramaaBa aaleaacqaHdpWCaeqaaaaa@38B2@ is mm. OptiStruct takes care of the unit system for R m MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOuamaaBa aaleaacaWGTbaabeaaaaa@37EC@ and G σ MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4ramaaBa aaleaacqaHdpWCaeqaaaaa@38B2@ through stress units defined in MATFAT and stress unit and length unit defined in FATPARM. a G MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyyamaaBa aaleaacaWGhbaabeaaaaa@37D5@ and b G MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyyamaaBa aaleaacaWGhbaabeaaaaa@37D5@ values are user input in MATFAT after keyword STSGRD. Since the stress gradient has to be calculated in length dimension of mm, define the length units so that OptiStruct can properly locate a point that is 1mm below the surface. If G σ MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4ramaaBa aaleaacqaHdpWCaeqaaaaa@38B2@ is negative, n σ MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOBamaaBa aaleaacqaHdpWCaeqaaaaa@38D9@ is set to 1.0. If G σ MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4ramaaBa aaleaacqaHdpWCaeqaaaaa@38B2@ is greater than 100 mm-1, n σ MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOBamaaBa aaleaacqaHdpWCaeqaaaaa@38D9@ is set to 1.0 with a warning message.

User-defined Relationship

User-defined relationship between n σ MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOBamaaBa aaleaacqaHdpWCaeqaaaaa@38D9@ and G σ MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4ramaaBa aaleaacqaHdpWCaeqaaaaa@38B2@ can be specified through TABLES1 Bulk Data. Pairs of (xi,yi) = ( G σ MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4ramaaBa aaleaacqaHdpWCaeqaaaaa@38B2@ , n σ MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOBamaaBa aaleaacqaHdpWCaeqaaaaa@38D9@ ) can be defined on the TABLES1 entry. A TABLES1 that defines the relationship between n σ MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOBamaaBa aaleaacqaHdpWCaeqaaaaa@38D9@ and G σ MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4ramaaBa aaleaacqaHdpWCaeqaaaaa@38B2@ should be referenced in MATFAT after keyword STSGRD. If G σ MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4ramaaBa aaleaacqaHdpWCaeqaaaaa@38B2@ falls outside the range of xi, extrapolation behavior follows usual TABLES1 behavior. This means that n σ MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOBamaaBa aaleaacqaHdpWCaeqaaaaa@38D9@ can be lower than 1.0 when G σ MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4ramaaBa aaleaacqaHdpWCaeqaaaaa@38B2@ is negative depending on how G σ MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4ramaaBa aaleaacqaHdpWCaeqaaaaa@38B2@ is treated when being negative or greater than 100mm-1. The user-defined relationship takes precedence over the one in FKM guidelines.

Critical Distance Method

To activate Stress Gradient effect using Critical Distance method, the GRD field on FATPARM should be set to GRDCD.

Small stress concentration features or geometries with high stress gradients are less effective in fatigue than larger features or smaller gradients with the same maximum stress. A plate with a small hole, say 0.1mm, will have a much longer fatigue life than one with a large hole of 10mm even though both plates have the same stress concentration factor and maximum stress. In conventional fatigue analysis, the stress gradient effect is taken into account by using an empirical fatigue notch factor, Kf, rather than the stress concentration factor Kt. Since there is no concept of a Kt or nominal stress in a finite element model stress gradient effects are considered directly. All of the holes have the same maximum stress, three times the nominal stress.


Figure 1. Stress distribution in a plate for three different hole sizes

Figure 1 that the stresses are independent of size only at the edge of the hole and vary far from the hole. The dashed line in the figure is drawn at 0.5mm. Here the stresses increase as the size of the hole increases. Suppose crack nucleation mechanisms result in a crack with a size of 0.5mm. For the smallest hole, 0.1mm, the stress available for continued growth is only 100 MPa, the nominal stress. The same size crack is subjected to a stress of 275 MPa in the larger hole, nearly equal to the maximum stress.

For nucleation of a crack around a hole of different sizes, it is useful to think about a process zone for crack nucleation. Materials are not continuous and homogeneous on the size scale that crack nucleation mechanisms operate. The grain size of the material is a convenient way to visualize the fatigue process zone. Figure 2 shows the grain size superimposed on the stress distribution from Figure 1. What is the stress in the process zone? A simple first approximation would be to take the stress in the center of the grain. Thus, a stress of 275 MPa would be used to compute the fatigue life of a 10mm hole and a stress of 100 MPa would be used for the 0.1mm hole.


Figure 2.

The modern view of fatigue is that when a material is stressed at the fatigue limit a microcrack will form but not grow outside of the process zone. Stress gradient effects are included in the fatigue analysis in a very simple and straightforward manner. In Critical Distance method, stresses and strains at a distance L/2 (Point Method) from the surface are used rather than the surface stresses and strains. For solid elements, the stress and strain at L/2 below surface is an interpolated stress and strain from grid point stresses and strains of an element of interest. In case of 2nd order solid elements, only grid point stresses and strains at corners are used for interpolation.

The critical distance can be expressed in terms of the threshold stress intensity, Δ K T H MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeuiLdqKaam 4samaaBaaaleaacaWGubGaamisaaqabaaaaa@39FF@ , and fatigue limit range, Δ σ F L MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeuiLdqKaeq 4Wdm3aaSbaaSqaaiaadAeacaWGmbaabeaaaaa@3AE8@ , as:(3)
L = 1 π ( Δ K T H Δ σ F L ) 2 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamitaiabg2 da9maalaaabaGaaGymaaqaaiabec8aWbaadaqadaqaamaalaaabaGa euiLdqKaam4samaaBaaaleaacaWGubGaamisaaqabaaakeaacqqHuo arcqaHdpWCdaWgaaWcbaGaamOraiaadYeaaeqaaaaaaOGaayjkaiaa wMcaamaaCaaaleqabaGaaGOmaaaaaaa@45E5@
The critical distance is a unique material property. If the critical distance of the material in use is known, user can input the critical distance in MATFAT after keyword STSGRD. When you input the critical distance, it is important to define dimension of length in MATFAT as well. Computing the critical distance from the threshold stress intensity, however, is difficult because the threshold stress intensity, particularly for small microcracks, is usually unknown. Fortunately, there is a good direct correlation between the critical distance and fatigue.(4)
L= ( 7× 10 4 E Δ σ FL ) 1.92 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamitaiabg2 da9maabmaabaGaaG4naiabgEna0kaaigdacaaIWaWaaWbaaSqabeaa cqGHsislcaaI0aaaaOWaaSaaaeaacaWGfbaabaGaeuiLdqKaeq4Wdm 3aaSbaaSqaaiaadAeacaWGmbaabeaaaaaakiaawIcacaGLPaaadaah aaWcbeqaaiaaigdacaGGUaGaaGyoaiaaikdaaaaaaa@4874@

If you do not directly input the critical distance, OptiStruct uses Equation 4 to estimate the critical distance in SN fatigue analysis. Fatigue limit Δ σ F L MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeuiLdqKaeq 4Wdm3aaSbaaSqaaiaadAeacaWGmbaabeaaaaa@3AE8@ is taken after the SN curve adjustment. Dimension of L is mm.

In EN fatigue analysis, the fatigue limit Δ σ F L MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeuiLdqKaeq 4Wdm3aaSbaaSqaaiaadAeacaWGmbaabeaaaaa@3AE8@ is approximated in the following manner.(5)
Δ σ FL =2E e nc MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeuiLdqKaeq 4Wdm3aaSbaaSqaaiaadAeacaWGmbaabeaakiabg2da9iaaikdacaWG fbGaamyzamaaBaaaleaacaWGUbGaam4yaaqabaaaaa@406F@
(6)
e nc = S ' f E × N c b MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyzamaaBa aaleaacaWGUbGaam4yaaqabaGccqGH9aqpdaWcaaqaaiaadofacaGG NaWaaSbaaSqaaiaadAgaaeqaaaGcbaGaamyraaaacqGHxdaTcaWGob Waa0baaSqaaiaadogaaeaacaWGIbaaaaaa@425C@
Where,
S ' f MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4uaiaacE cadaWgaaWcbaGaamOzaaqabaaaaa@3891@
Fatigue strength coefficient.
N c MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOtamaaBa aaleaacaWGJbaabeaaaaa@37DE@
Reversal limit of endurance.
E MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyraaaa@36C1@
Young’s modulus.

If Δ σ F L MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeuiLdqKaeq 4Wdm3aaSbaaSqaaiaadAeacaWGmbaabeaaaaa@3AE8@ is 0 or the calculated L MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyraaaa@36C1@ is greater than 0.2mm, L MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyraaaa@36C1@ will be set to 0.2mm. In case of shell elements, the maximum calculated L MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyraaaa@36C1@ is thickness/4.

Input to Activate Stress Gradient Effect

Choose a method (FKM guideline or Critical Distance) to use on the GRD field after keyword STRESS in FATPARM. If FKM guideline method is chosen, the equivalent stress σ eq MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeq4Wdm3aaS baaSqaaiaadwgacaWGXbaabeaaaaa@39C6@ method to calculate stress gradient should be specified on the SCBFKM field in FATPARM. Material properties required for stress gradient effect are to be input after keyword STSGRD in MATFAT.