From c3f93fbb7eaee4f447446fa892205c8456c0a96f Mon Sep 17 00:00:00 2001
From: Christoph Lehmann <christoph.lehmann@ufz.de>
Date: Fri, 23 Sep 2016 12:03:38 +0200
Subject: [PATCH] [PL] SmallDef init axisymmetric BMatrices
---
.../SmallDeformation/SmallDeformationFEM.h | 17 +++++++++++------
1 file changed, 11 insertions(+), 6 deletions(-)
diff --git a/ProcessLib/SmallDeformation/SmallDeformationFEM.h b/ProcessLib/SmallDeformation/SmallDeformationFEM.h
index 35405c7dea3..a86ca581c63 100644
--- a/ProcessLib/SmallDeformation/SmallDeformationFEM.h
+++ b/ProcessLib/SmallDeformation/SmallDeformationFEM.h
@@ -160,15 +160,20 @@ public:
{
_ip_data.emplace_back(*_process_data._material);
auto& ip_data = _ip_data[ip];
- ip_data._detJ = shape_matrices[ip].detJ;
- ip_data._integralMeasure = shape_matrices[ip].integralMeasure;
+ auto const& sm = shape_matrices[ip];
+ ip_data._detJ = sm.detJ;
+ ip_data._integralMeasure = sm.integralMeasure;
ip_data._b_matrices.resize(
KelvinVectorDimensions<DisplacementDim>::value,
ShapeFunction::NPOINTS * DisplacementDim);
- LinearBMatrix::computeBMatrix<
- DisplacementDim, ShapeFunction::NPOINTS,
- typename ShapeMatricesType::GlobalDimNodalMatrixType,
- BMatrixType>(shape_matrices[ip].dNdx, ip_data._b_matrices);
+
+ auto const x_coord =
+ interpolateXCoordinate<ShapeFunction, ShapeMatricesType>(e,
+ sm.N);
+ LinearBMatrix::computeBMatrix<DisplacementDim,
+ ShapeFunction::NPOINTS>(
+ sm.dNdx, ip_data._b_matrices, is_axially_symmetric, sm.N,
+ x_coord);
ip_data._sigma.resize(
KelvinVectorDimensions<DisplacementDim>::value);
--
GitLab