Yeoh (/MAT/LAW94)

LAW94 is a hyperelastic material model that can be used to describe incompressible materials.

The strain energy density function of LAW94 only depends on the first strain invariant and is computed as:(1)
W= i=1 3 [ C i0 ( I ¯ 1 3 ) i W( I ¯ 1 ) + 1 D i ( J1 ) 2i U(J) ] MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4vaiabg2 da9maaqahabaWaamWaaeaadaagaaqaaiaadoeadaWgaaWcbaGaamyA aiaaicdaaeqaaOWaaeWaaeaaceWGjbGbaebadaWgaaWcbaGaaGymaa qabaGccqGHsislcaaIZaaacaGLOaGaayzkaaWaaWbaaSqabeaacaWG PbaaaaqaaiaadEfacaGGOaGabmysayaaraWaaSbaaWqaaiaaigdaae qaaSGaaiykaaGccaGL44pacqGHRaWkdaagaaqaamaalaaabaGaaGym aaqaaiaadseadaWgaaWcbaGaamyAaaqabaaaaOWaaeWaaeaacaWGkb GaeyOeI0IaaGymaaGaayjkaiaawMcaamaaCaaaleqabaGaaGOmaiaa dMgaaaaabaGaamyvaiaacIcacaWGkbGaaiykaaGccaGL44paaiaawU facaGLDbaaaSqaaiaadMgacqGH9aqpcaaIXaaabaGaaG4maaqdcqGH ris5aaaa@5CAE@
Where,
I ¯ 1 = λ ¯ 1 2 + λ ¯ 2 2 + λ ¯ 3 2 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGabmysayaara WaaSbaaSqaaiaaigdaaeqaaOGaeyypa0Jafq4UdWMbaebadaqhaaWc baGaaGymaaqaaiaaikdaaaGccqGHRaWkcuaH7oaBgaqeamaaDaaale aacaaIYaaabaGaaGOmaaaakiabgUcaRiqbeU7aSzaaraWaa0baaSqa aiaaiodaaeaacaaIYaaaaaaa@44FF@
First strain invariant
λ ¯ i = J 1 3 λ i MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGafq4UdWMbae badaWgaaWcbaGaamyAaaqabaGccqGH9aqpcaWGkbWaaWbaaSqabeaa cqGHsisldaWcaaqaaiaaigdaaeaacaaIZaaaaaaakiabeU7aSnaaBa aaleaacaWGPbaabeaaaaa@4036@
Deviatoric stretch
The Cauchy stress is:(2)
σ i = λ i J W λ i MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeq4Wdm3aaS baaSqaaiaadMgaaeqaaOGaeyypa0ZaaSaaaeaacqaH7oaBdaWgaaWc baGaamyAaaqabaaakeaacaWGkbaaamaalaaabaGaeyOaIyRaam4vaa qaaiabgkGi2kabeU7aSnaaBaaaleaacaWGPbaabeaaaaaaaa@4421@

Material Parameters

For incompressible materials with i MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyAaaaa@36E4@ =1 only and D 1 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiramaaBa aaleaacaaIXaaabeaaaaa@37A7@ are input and the Yeoh model is reduced to a Neo-Hookean model.
C 10 , C 20 , C 30 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4qamaaBa aaleaacaaIXaGaaGimaaqabaGccaGGSaGaam4qamaaBaaaleaacaaI YaGaaGimaaqabaGccaGGSaGaam4qamaaBaaaleaacaaIZaGaaGimaa qabaaaaa@3EA8@
Material constants specify the deviatoric part (shape change) of the material
D 1 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiramaaBa aaleaacaaIXaaabeaaaaa@37A7@ , D 2 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiramaaBa aaleaacaaIXaaabeaaaaa@37A7@ , D 3 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiramaaBa aaleaacaaIXaaabeaaaaa@37A7@
Parameters specify the volumetric change of the material

These six material constants need to be calculated by curve fitting material test data. RD-E: 5600 Hyperelastic Material with Curve Input includes a Yeoh fitting Compose script for uniaxial test data. The Yeoh material model has been shown to model all deformation models, even if the curve fit was obtained using only uniaxial test data.

The initial shear modulus and the bulk modulus are:(3)
μ=2 C 10 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqiVd0Maey ypa0JaaGOmaiabgwSixlaadoeadaWgaaWcbaGaaGymaiaaicdaaeqa aaaa@3E22@
and(4)
K= 2 D 1 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4saiabg2 da9maalaaabaGaaGOmaaqaaiaadseadaWgaaWcbaGaaGymaaqabaaa aaaa@3A49@

Poisson's Ratio and Material Incompressibility

LAW94 is available only as an incompressible material model.

If D 1 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiramaaBa aaleaacaaIXaaabeaaaaa@37A7@ = 0, an incompressible material is considered where, ν = 0.495 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqyVd4Maey ypa0JaaGimaiaac6cacaaI0aGaaGyoaiaaiwdaaaa@3C61@ and D 1 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiramaaBa aaleaacaaIXaaabeaaaaa@37A7@ is calculated as:(5)
D 1 = 3 ( 1 2 v ) μ ( 1 + v ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiramaaBa aaleaacaaIXaaabeaakiabg2da9maalaaabaGaaG4maiaacIcacaaI XaGaeyOeI0IaaGOmaiaadAhacaGGPaaabaGaeqiVd0Maaiikaiaaig dacqGHRaWkcaWG2bGaaiykaaaaaaa@43E2@

References

1 Yeoh, O. H. "Some forms of the strain energy function for rubber." Rubber Chemistry and technology 66, no. 5 (1993): 754-771