FEvoR / master

Master working branch for FEvoR
Other
2 stars 0 forks source link

Bulk strain calculation is goofy. #2

Closed jhkennedy closed 9 years ago

jhkennedy commented 9 years ago

When testing the ability of FEvoR to calculate a "good" enhancement factor, I noticed bulk strain rate is not doing what is expected.

For shear stress, it's magnitude should be larger for a single maximum than for an isotropic distribution. E.g.:

./bin/enhancementFactor 10000.0 -10.0 0.0 0.0 500.0 0.0 0.0 0.0 100.0

is a run with a perfect single maximum, T = -10 degrees, only shear stress -- Szz = 500 Pa, for a 100 year time step.

M     was: 1.2688e-16
M_iso was: 1.3922e-16
Enhancement factor was: 9.1132e-01

where the enhancement factor is

E = M / M_iso