From be5899a9d7db1daf874e0f214fb1e02eb259d37a Mon Sep 17 00:00:00 2001
From: Wenqing Wang <wenqing.wang@ufz.de>
Date: Wed, 23 Jun 2021 12:13:42 +0200
Subject: [PATCH] [LIE] Use ElementCoordinatesMappingLocal to compute the
rotation matrix
---
ProcessLib/LIE/Common/FractureProperty.h | 25 +++++++++++--------
...llDeformationLocalAssemblerFracture-impl.h | 1 +
2 files changed, 15 insertions(+), 11 deletions(-)
diff --git a/ProcessLib/LIE/Common/FractureProperty.h b/ProcessLib/LIE/Common/FractureProperty.h
index 865f51130a9..f33125b6a91 100644
--- a/ProcessLib/LIE/Common/FractureProperty.h
+++ b/ProcessLib/LIE/Common/FractureProperty.h
@@ -10,19 +10,15 @@
#pragma once
-#include <memory>
#include <Eigen/Eigen>
-
-#include "MaterialLib/FractureModels/Permeability/Permeability.h"
+#include <memory>
#include "BranchProperty.h"
#include "JunctionProperty.h"
-#include "Utils.h"
+#include "MaterialLib/FractureModels/Permeability/Permeability.h"
+#include "MeshLib/ElementCoordinatesMappingLocal.h"
+#include "MeshLib/Elements/Element.h"
-namespace MeshLib
-{
-class Element;
-}
namespace ParameterLib
{
template <typename T>
@@ -87,9 +83,16 @@ inline void setFractureProperty(int const dim, MeshLib::Element const& e,
{
frac_prop.point_on_fracture[j] = getCenterOfGravity(e).getCoords()[j];
}
- computeNormalVector(e, dim, n);
- frac_prop.R.resize(dim, dim);
- computeRotationMatrix(e, n, dim, frac_prop.R);
+
+ const MeshLib::ElementCoordinatesMappingLocal ele_local_coord(e, dim);
+
+ // Global to local rotation matrix:
+ Eigen::MatrixXd const global2local_rotation =
+ ele_local_coord.getRotationMatrixToGlobal().transpose();
+ n = global2local_rotation.row(dim - 1);
+
+ frac_prop.R = global2local_rotation.topLeftCorner(dim, dim);
+
DBUG("Normal vector of the fracture element {:d}: [{:g}, {:g}, {:g}]",
e.getID(), n[0], n[1], n[2]);
}
diff --git a/ProcessLib/LIE/SmallDeformation/LocalAssembler/SmallDeformationLocalAssemblerFracture-impl.h b/ProcessLib/LIE/SmallDeformation/LocalAssembler/SmallDeformationLocalAssemblerFracture-impl.h
index 9290c46e489..575ef449856 100644
--- a/ProcessLib/LIE/SmallDeformation/LocalAssembler/SmallDeformationLocalAssemblerFracture-impl.h
+++ b/ProcessLib/LIE/SmallDeformation/LocalAssembler/SmallDeformationLocalAssemblerFracture-impl.h
@@ -15,6 +15,7 @@
#include "MathLib/LinAlg/Eigen/EigenMapTools.h"
#include "NumLib/Fem/InitShapeMatrices.h"
#include "ProcessLib/LIE/Common/LevelSetFunction.h"
+#include "ProcessLib/LIE/Common/Utils.h"
#include "SmallDeformationLocalAssemblerFracture.h"
namespace ProcessLib
--
GitLab