From a8be5ce92c2f24308970b54e311a6f9e2fe2f529 Mon Sep 17 00:00:00 2001
From: FZill <florian.zill@ufz.de>
Date: Tue, 19 Jan 2021 11:57:12 +0100
Subject: [PATCH] Revert "[HM] added permeability-strain derivative"
This reverts commit fe60a51b1a9a895ce8c88cfe9aac83a7dd42ea8c.
---
.../HydroMechanics/HydroMechanicsFEM-impl.h | 15 +--------------
1 file changed, 1 insertion(+), 14 deletions(-)
diff --git a/ProcessLib/HydroMechanics/HydroMechanicsFEM-impl.h b/ProcessLib/HydroMechanics/HydroMechanicsFEM-impl.h
index a692d09762c..061901d7b4d 100644
--- a/ProcessLib/HydroMechanics/HydroMechanicsFEM-impl.h
+++ b/ProcessLib/HydroMechanics/HydroMechanicsFEM-impl.h
@@ -175,12 +175,6 @@ void HydroMechanicsLocalAssembler<ShapeFunctionDisplacement,
pressure_size, displacement_size>::Zero(pressure_size,
displacement_size);
- typename ShapeMatricesTypeDisplacement::template MatrixType<
- pressure_size, displacement_size>
- Kpu_k = ShapeMatricesTypeDisplacement::template MatrixType<
- pressure_size, displacement_size>::Zero(pressure_size,
- displacement_size);
-
MaterialLib::Solids::MechanicsBase<DisplacementDim> const& solid_material =
*_process_data.solid_materials[0];
@@ -279,9 +273,6 @@ void HydroMechanicsLocalAssembler<ShapeFunctionDisplacement,
auto const K = MPL::formEigenTensor<DisplacementDim>(
medium->property(MPL::PropertyType::permeability)
.value(vars, x_position, t, dt));
- auto const dkde = MPL::formEigenTensor<DisplacementDim>(
- medium->property(MPL::PropertyType::permeability)
- .dValue(vars, MPL::Variable::mechanical_strain, x_position, t, dt));
auto const K_over_mu = K / mu;
@@ -330,10 +321,6 @@ void HydroMechanicsLocalAssembler<ShapeFunctionDisplacement,
//
Kpu.noalias() +=
rho_fr * alpha * N_p.transpose() * identity2.transpose() * B * w;
-
- Kpu_k.noalias() += dNdx_p.transpose() * dkde *
- (dNdx_p * p - rho_fr * b) * identity2.transpose() *
- B * rho_fr / mu;
}
// displacement equation, pressure part
local_Jac
@@ -362,7 +349,7 @@ void HydroMechanicsLocalAssembler<ShapeFunctionDisplacement,
local_Jac
.template block<pressure_size, displacement_size>(pressure_index,
displacement_index)
- .noalias() = Kpu / dt + Kpu_k;
+ .noalias() = Kpu / dt;
// pressure equation
local_rhs.template segment<pressure_size>(pressure_index).noalias() -=
--
GitLab