LLNL / serac

Serac is a high order nonlinear thermomechanical simulation code
BSD 3-Clause "New" or "Revised" License
173 stars 30 forks source link

Fix error in conversion of Mandel stress to Cauchy stress #1166

Closed btalamini closed 6 days ago

btalamini commented 1 week ago

One-line fix. Actually, "fix" is not quite the right word, since the results obtained before were numerically correct by a coincidence. The conversion from the Mandel stress (conjugate to the elastic log strain) to the Cauchy stress was incorrect for the general case. However, for the special case of isotropic behavior (as in the implemented model), the expression in the code would work because of some matrix commuting and cancellation. Sorry, I should have checked this more closely when I implemented the model.