/ANIM/SHELL/DAMA

Engine Keyword Generates animation files containing a damage value for each layer of a shell element. The damage value is the maximum of damage over time of all failure criteria associated to the material.

Format

/ANIM/SHELL/DAMA/Keyword4/Keyword5

Definition

Field Contents SI Unit Example
Keyword4 Output location:
N
Layer number N
UPPER
Upper layer
LOWER
Lower layer (number of integration point is equal to one)
ALL
All the layers
Keyword5 Result type or location:
NIP
Number of integration points through layer N 5
UPPER
Upper integration point (through NIP) of the upper layer N
LOWER
Lower integration point (through NIP) of the lower layer N
ALL
All the integration points (NIP) through the layer N

Comments

  1. In case of BATOZ shell element, the value is the average value of the Gauss points in the layer.

    This applies only to failure criteria defined with Failure Models keywords, failure defined inside material laws (/MAT/LAW52 (GURSON)) is currently not available.

  2. The damage value, D is the maximum value of all failure criteria of the model.

    The damage value, which is displayed, is the maximum damage value over time.

    In case of composite shell element property TYPE51, the damage value for one layer is computed as the maximum value over the all integration points of its layer.

    D ( T ) = max t T ( D ( t ) )

  3. The damage value, D is 0 ≤ D ≤ 1. The status for fracture is:

    Free, if 0 ≤ D < 1

    Failure, if D = 1

  4. D is computed for every failure criteria as follows:
    • Strain based failure (/FAIL/BIQUAD)
      D = Δ ε p ε f
    • Johnson-Cook failure (/FAIL/JOHNSON)
      D = Δ ɛ p ɛ f
    • Tuler-Butcher failure (/FAIL/TBUTCHER)
      D = 0 t ( σ σ r ) λ d t K
    • Wilkins failure (/FAIL/WILKINS)
      D = W 1 W 2 d ε p D f
    • FLD failure (/FAIL/FLD)
      D = M a x t i m e ( ε m a j o r ε lim )
    Where,
    ε lim
    Major strain as limit from FLD function
    Figure 1.
    • BAO-XUE-Wierzbicki failure (/FAIL/WIERZBICKI)
      D = Δ ε p ε ¯ f
    • Strain failure model (/FAIL/TENSSTRAIN)
      D = Max time ( ɛ 1 ɛ t 1 ɛ t 2 ɛ t 1 )
    • Strain failure model (/FAIL/FABRIC)
      D = M a x ( ε 1 ε f 1 ε r 1 ε f 1 , ε 2 ε f 2 ε r 2 ε f 2 ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGebGaey ypa0JaamytaiaadggacaWG4bWaaeWaaeaadaWcaaqaaiabew7aLnaa BaaaleaacaaIXaaabeaakiabgkHiTiabew7aLnaaBaaaleaacaWGMb GaaGymaaqabaaakeaacqaH1oqzdaWgaaWcbaGaamOCaiaaigdaaeqa aOGaeyOeI0IaeqyTdu2aaSbaaSqaaiaadAgacaaIXaaabeaaaaGcca GGSaWaaSaaaeaacqaH1oqzdaWgaaWcbaGaaGOmaaqabaGccqGHsisl cqaH1oqzdaWgaaWcbaGaamOzaiaaikdaaeqaaaGcbaGaeqyTdu2aaS baaSqaaiaadkhacaaIYaaabeaakiabgkHiTiabew7aLnaaBaaaleaa caWGMbGaaGOmaaqabaaaaaGccaGLOaGaayzkaaaaaa@5C57@
    • Energy density failure model (/FAIL/ENERGY):
      D = Max time ( E E 1 E 2 E 1 )
    • Chang failure (/FAIL/CHANG)

      The maximum damage for different failure mode:

      D = Max ( e f 2 , e c 2 , e m 2 , e d 2 )

    • Hashin Composite failure (/FAIL/HASHIN)

      The maximum damage for different failure mode

      For uni-directional lamina model:

      D = M a x ( F 1 , F 2 , F 3 , F 4 , F 5 ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiraiabg2 da9iaad2eacaWGHbGaamiEamaabmaabaGaamOramaaBaaaleaacaaI XaaabeaakiaacYcacaWGgbWaaSbaaSqaaiaaikdaaeqaaOGaaiilai aadAeadaWgaaWcbaGaaG4maaqabaGccaGGSaGaamOramaaBaaaleaa caaI0aaabeaakiaacYcacaWGgbWaaSbaaSqaaiaaiwdaaeqaaaGcca GLOaGaayzkaaaaaa@4779@

      For fabric lamina model:

      D = M a x ( F 1 , F 2 , F 3 , F 4 , F 5 , F 6 , F 7 ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiraiabg2 da9iaad2eacaWGHbGaamiEamaabmaabaGaamOramaaBaaaleaacaaI XaaabeaakiaacYcacaWGgbWaaSbaaSqaaiaaikdaaeqaaOGaaiilai aadAeadaWgaaWcbaGaaG4maaqabaGccaGGSaGaamOramaaBaaaleaa caaI0aaabeaakiaacYcacaWGgbWaaSbaaSqaaiaaiwdaaeqaaOGaai ilaiaadAeadaWgaaWcbaGaaGOnaaqabaGccaGGSaGaamOramaaBaaa leaacaaI3aaabeaaaOGaayjkaiaawMcaaaaa@4C5C@

    • Puck Composite failure (/FAIL/PUCK)

      The maximum damage for different failure mode:

      D = Max ( e f ( tensile ) , e f ( compression ) , e f ( ModeA ) , e f ( ModeB ) , e f ( ModeC ) )

    • Ladeveze failure (/FAIL/LAD_DAMA)

      Only used for Plyxfem in /PROP/TYPE17:

      d ˙ = k a [ 1 exp ( a w ( Y ) d ) ] MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGabmizayaaca Gaeyypa0ZaaSaaaeaacaWGRbaabaGaamyyaaaadaWadaqaaiaaigda cqGHsislciGGLbGaaiiEaiaacchadaqadaqaaiabgkHiTiaadggada aadaqaaiaadEhadaqadaqaaiaadMfaaiaawIcacaGLPaaacqGHsisl caWGKbaacaGLPmIaayPkJaaacaGLOaGaayzkaaaacaGLBbGaayzxaa aaaa@4AAF@

    • Strain Failure Model with dependence on Lode angle (/FAIL/TAB1)
      D = Δ D D c r i t
    • NXT Failure Model (/FAIL/NXT)
      D = λ f 2
    • Strain based Ductile Failure Model (/FAIL/HC_DSSE)
      D = Δ ε ¯ p ε ¯ H C p r ( η ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiaadseacqGH9a qpdaaeabqaamaalaaabaGaeuiLdqKafqyTduMbaebadaWgaaWcbaGa amiCaaqabaaakeaacuaH1oqzgaqeamaaDaaaleaacaWGibGaam4qaa qaaiaadchacaWGYbaaaOWaaeWaaeaacqaH3oaAaiaawIcacaGLPaaa aaaaleqabeqdcqGHris5aaaa@46DB@
  5. This option is available only if Keyword4 = N (for each defined layer N allows to compute shell damage value within each integration point NIP of layer N). This option is only available for /PROP/TYPE51 and PCOMPP (/PROP/TYPE51 or /PROP/PCOMPP with multiple integration points through each layer).

    Nevertheless, Keyword4=ALL combined with Keyword5=ALL is allowed.