From a3ef51389f7123ac66711dfb765c0056efbc6042 Mon Sep 17 00:00:00 2001
From: FZill <florian.zill@ufz.de>
Date: Mon, 10 May 2021 14:01:09 +0200
Subject: [PATCH] [MPL] Fix strain in OrthotropicEmbeddedFracturePermeability
---
.../Properties/OrthotropicEmbeddedFracturePermeability.cpp | 6 ++++--
1 file changed, 4 insertions(+), 2 deletions(-)
diff --git a/MaterialLib/MPL/Properties/OrthotropicEmbeddedFracturePermeability.cpp b/MaterialLib/MPL/Properties/OrthotropicEmbeddedFracturePermeability.cpp
index b1470477bd7..f2c9c7392ae 100644
--- a/MaterialLib/MPL/Properties/OrthotropicEmbeddedFracturePermeability.cpp
+++ b/MaterialLib/MPL/Properties/OrthotropicEmbeddedFracturePermeability.cpp
@@ -54,8 +54,9 @@ OrthotropicEmbeddedFracturePermeability<DisplacementDim>::value(
ParameterLib::SpatialPosition const& pos, double const t,
double const /*dt*/) const
{
- auto const eps = formEigenTensor<3>(std::get<SymmetricTensor>(
- variable_array[static_cast<int>(Variable::mechanical_strain)]));
+ auto const eps = MathLib::KelvinVector::kelvinVectorToTensor(
+ std::get<MathLib::KelvinVector::KelvinVectorType<DisplacementDim>>(
+ variable_array[static_cast<int>(Variable::mechanical_strain)]));
double const k = std::get<double>(fromVector(_k(t, pos)));
double const _b0 = std::sqrt(12.0 * k);
@@ -76,6 +77,7 @@ OrthotropicEmbeddedFracturePermeability<DisplacementDim>::value(
result += H_de * (b_f / _a[i]) * ((b_f * b_f / 12.0) - k) *
(Eigen::Matrix3d::Identity() - n_i * n_i.transpose());
}
+
return result.template topLeftCorner<DisplacementDim, DisplacementDim>()
.eval();
}
--
GitLab