Skip to content
Snippets Groups Projects
Commit 49638df7 authored by Oliver Sander's avatar Oliver Sander Committed by sander
Browse files

[bugfix] Add missing skew part in finite-strain membrane energy

This bug has now influence as long as mu_c is zero, but for
non-zero mu_c there was one term missing.

[[Imported from SVN: r9865]]
parent 929e8224
No related branches found
No related tags found
No related merge requests found
......@@ -204,7 +204,7 @@ public:
RT detU = U.determinant();
return mu_ * sym(UMinus1).frobenius_norm2()
return mu_ * sym(UMinus1).frobenius_norm2() + mu_c_ * skew(UMinus1).frobenius_norm2()
+ (mu_*lambda_)/(2*mu_ + lambda_) * 0.5 * ((detU-1)*(detU-1) + (1.0/detU -1)*(1.0/detU -1));
}
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment