Fixed small bug in evaluation of elastic tangent in cubic-CLE and ort…

…ho-CLE materials
febiosoftware · Oct 3, 2024 · ce496e4 · ce496e4
1 parent b27cc08
commit ce496e4
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/FEBioMech/FECubicCLE.cpp b/FEBioMech/FECubicCLE.cpp
@@ -181,7 +181,7 @@ tens4ds FECubicCLE::Tangent(FEMaterialPoint& mp)
  for (i=0; i<3; ++i) {
  c += dyad4s(A[i], B)*mu[i];
  for (j=0; j<3; ++j) {
- c += dyad1s(A[i],A[j])*lam[i][j];
+ c += dyad1s(A[i],A[j])*lam[i][j]/2;
  }
  }
  c /= J;

diff --git a/FEBioMech/FEOrthotropicCLE.cpp b/FEBioMech/FEOrthotropicCLE.cpp
@@ -188,7 +188,7 @@ tens4ds FEOrthotropicCLE::Tangent(FEMaterialPoint& mp)
  for (i=0; i<3; ++i) {
  c += dyad4s(A[i], B)*mu[i];
  for (j=0; j<3; ++j) {
- c += dyad1s(A[i],A[j])*lam[i][j];
+ c += dyad1s(A[i],A[j])*lam[i][j]/2;
  }
  }
  c /= J;