From 26a09dcbbe4b58902b473aaee8acddcb55cd942f Mon Sep 17 00:00:00 2001
From: Zheng <ichbinzty@gmail.com>
Date: Thu, 19 Jan 2017 00:50:42 +0100
Subject: [PATCH] add the integralMeasure term in HM process to fix the
axsymmetric problem
---
ProcessLib/HydroMechanics/HydroMechanicsFEM.h | 3 ++-
1 file changed, 2 insertions(+), 1 deletion(-)
diff --git a/ProcessLib/HydroMechanics/HydroMechanicsFEM.h b/ProcessLib/HydroMechanics/HydroMechanicsFEM.h
index c39a47baf5d..6275c1d5c89 100644
--- a/ProcessLib/HydroMechanics/HydroMechanicsFEM.h
+++ b/ProcessLib/HydroMechanics/HydroMechanicsFEM.h
@@ -144,8 +144,9 @@ public:
// displacement (subscript u)
_ip_data.emplace_back(*_process_data.material);
auto& ip_data = _ip_data[ip];
+ auto const& sm = shape_matrices_u[ip];
_ip_data[ip].integration_weight =
- _integration_method.getWeightedPoint(ip).getWeight() *
+ _integration_method.getWeightedPoint(ip).getWeight() * sm.integralMeasure *
shape_matrices_u[ip].detJ;
ip_data.b_matrices.resize(
kelvin_vector_size,
--
GitLab