Detailed Description

template<typename ShapeFunctionDisplacement, int GlobalDim>
class ProcessLib::NormalTractionBoundaryCondition::NormalTractionBoundaryConditionLocalAssembler< ShapeFunctionDisplacement, GlobalDim >

Definition at line 51 of file NormalTractionBoundaryConditionLocalAssembler.h.

#include <NormalTractionBoundaryConditionLocalAssembler.h>

Inheritance diagram for ProcessLib::NormalTractionBoundaryCondition::NormalTractionBoundaryConditionLocalAssembler< ShapeFunctionDisplacement, GlobalDim >:

[legend]

Collaboration diagram for ProcessLib::NormalTractionBoundaryCondition::NormalTractionBoundaryConditionLocalAssembler< ShapeFunctionDisplacement, GlobalDim >:

[legend]

Public Types
using	ShapeMatricesType

using	GlobalDimVectorType = typename ShapeMatricesType::GlobalDimVectorType

using	NodalVectorType = typename ShapeMatricesType::NodalVectorType

Public Member Functions
	NormalTractionBoundaryConditionLocalAssembler (MeshLib::Element const &e, std::size_t const local_matrix_size, NumLib::GenericIntegrationMethod const &integration_method, bool const is_axially_symmetric, ParameterLib::Parameter< double > const &pressure)

void	assemble (std::size_t const id, NumLib::LocalToGlobalIndexMap const &dof_table_boundary, double const t, std::vector< GlobalVector * > const &, GlobalMatrix , GlobalVector &local_rhs, GlobalMatrix ) override

Public Member Functions inherited from ProcessLib::NormalTractionBoundaryCondition::NormalTractionBoundaryConditionLocalAssemblerInterface
virtual	~NormalTractionBoundaryConditionLocalAssemblerInterface ()=default

Public Attributes
	EIGEN_MAKE_ALIGNED_OPERATOR_NEW

Private Attributes
NumLib::GenericIntegrationMethod const &	_integration_method

ParameterLib::Parameter< double > const &	_pressure

std::vector< IntegrationPointData< ShapeMatricesType >, Eigen::aligned_allocator< IntegrationPointData< ShapeMatricesType > > >	_ip_data

ShapeMatricesType::template VectorType< displacement_size >	_local_rhs

MeshLib::Element const &	_element

Static Private Attributes
static const int	displacement_size

Member Typedef Documentation

◆ GlobalDimVectorType

template<typename ShapeFunctionDisplacement , int GlobalDim>

using ProcessLib::NormalTractionBoundaryCondition::NormalTractionBoundaryConditionLocalAssembler< ShapeFunctionDisplacement, GlobalDim >::GlobalDimVectorType = typename ShapeMatricesType::GlobalDimVectorType

Definition at line 57 of file NormalTractionBoundaryConditionLocalAssembler.h.

◆ NodalVectorType

template<typename ShapeFunctionDisplacement , int GlobalDim>

using ProcessLib::NormalTractionBoundaryCondition::NormalTractionBoundaryConditionLocalAssembler< ShapeFunctionDisplacement, GlobalDim >::NodalVectorType = typename ShapeMatricesType::NodalVectorType

Definition at line 58 of file NormalTractionBoundaryConditionLocalAssembler.h.

◆ ShapeMatricesType

template<typename ShapeFunctionDisplacement , int GlobalDim>

using ProcessLib::NormalTractionBoundaryCondition::NormalTractionBoundaryConditionLocalAssembler< ShapeFunctionDisplacement, GlobalDim >::ShapeMatricesType

Initial value:

ShapeMatrixPolicyType<ShapeFunctionDisplacement, GlobalDim>

EigenFixedShapeMatrixPolicy

Definition ShapeMatrixPolicy.h:73

Definition at line 55 of file NormalTractionBoundaryConditionLocalAssembler.h.

Constructor & Destructor Documentation

◆ NormalTractionBoundaryConditionLocalAssembler()

template<typename ShapeFunctionDisplacement , int GlobalDim>

ProcessLib::NormalTractionBoundaryCondition::NormalTractionBoundaryConditionLocalAssembler< ShapeFunctionDisplacement, GlobalDim >::NormalTractionBoundaryConditionLocalAssembler	(	MeshLib::Element const &	e,
		std::size_t const	local_matrix_size,
		NumLib::GenericIntegrationMethod const &	integration_method,
		bool const	is_axially_symmetric,
		ParameterLib::Parameter< double > const &	pressure )

inline

Definition at line 60 of file NormalTractionBoundaryConditionLocalAssembler.h.

        : _integration_method(integration_method),
          _pressure(pressure),
          _element(e)
    {
        _local_rhs.setZero(local_matrix_size);
 
        unsigned const n_integration_points =
            _integration_method.getNumberOfPoints();
 
        _ip_data.reserve(n_integration_points);
 
        auto const shape_matrices_u =
            NumLib::initShapeMatrices<ShapeFunctionDisplacement,
                                      ShapeMatricesType, GlobalDim>(
                e, is_axially_symmetric, _integration_method);
 
        GlobalDimVectorType element_normal(GlobalDim);
 
        // TODO Extend to rotated 2d meshes and line elements.
        if (e.getGeomType() == MeshLib::MeshElemType::LINE)
        {
            Eigen::Vector3d const v1 = e.getNode(1)->asEigenVector3d() -
                                       e.getNode(0)->asEigenVector3d();
            element_normal[GlobalDim - 1] = 0;
            element_normal[0] = -v1[1];
            element_normal[1] = v1[0];
            element_normal.normalize();
        }
        else
        {
            auto const n = MeshLib::FaceRule::getSurfaceNormal(e).normalized();
            for (int i = 0; i < GlobalDim; ++i)
            {
                element_normal[i] = n[i];
            }
        }
 
        for (unsigned ip = 0; ip < n_integration_points; ip++)
        {
            double const integration_weight =
                _integration_method.getWeightedPoint(ip).getWeight() *
                shape_matrices_u[ip].integralMeasure *
                shape_matrices_u[ip].detJ;
 
            _ip_data.emplace_back(shape_matrices_u[ip].N, element_normal,
                                  integration_weight);
        }
    }

References ProcessLib::NormalTractionBoundaryCondition::NormalTractionBoundaryConditionLocalAssembler< ShapeFunctionDisplacement, GlobalDim >::_integration_method, ProcessLib::NormalTractionBoundaryCondition::NormalTractionBoundaryConditionLocalAssembler< ShapeFunctionDisplacement, GlobalDim >::_ip_data, ProcessLib::NormalTractionBoundaryCondition::NormalTractionBoundaryConditionLocalAssembler< ShapeFunctionDisplacement, GlobalDim >::_local_rhs, MathLib::Point3d::asEigenVector3d(), MeshLib::Element::getGeomType(), MeshLib::Element::getNode(), NumLib::GenericIntegrationMethod::getNumberOfPoints(), MeshLib::FaceRule::getSurfaceNormal(), MathLib::WeightedPoint::getWeight(), NumLib::GenericIntegrationMethod::getWeightedPoint(), NumLib::initShapeMatrices(), and MeshLib::LINE.

Member Function Documentation

◆ assemble()

template<typename ShapeFunctionDisplacement , int GlobalDim>

void ProcessLib::NormalTractionBoundaryCondition::NormalTractionBoundaryConditionLocalAssembler< ShapeFunctionDisplacement, GlobalDim >::assemble	(	std::size_t const	id,
		NumLib::LocalToGlobalIndexMap const &	dof_table_boundary,
		double const	t,
		std::vector< GlobalVector * > const &	,
		GlobalMatrix *	,
		GlobalVector &	local_rhs,
		GlobalMatrix *	)

inlineoverridevirtual

Implements ProcessLib::NormalTractionBoundaryCondition::NormalTractionBoundaryConditionLocalAssemblerInterface.

Definition at line 115 of file NormalTractionBoundaryConditionLocalAssembler.h.

    {
        _local_rhs.setZero();
 
        unsigned const n_integration_points =
            _integration_method.getNumberOfPoints();
 
        NodalVectorType pressure =
            _pressure.getNodalValuesOnElement(_element, t);
 
        for (unsigned ip = 0; ip < n_integration_points; ip++)
        {
            auto const& w = _ip_data[ip].integration_weight;
            auto const& N = _ip_data[ip].N;
            auto const& n = _ip_data[ip].n;
 
            typename ShapeMatricesType::template MatrixType<GlobalDim,
                                                            displacement_size>
                N_u = ShapeMatricesType::template MatrixType<
                    GlobalDim, displacement_size>::Zero(GlobalDim,
                                                        displacement_size);
            for (int i = 0; i < GlobalDim; ++i)
            {
                N_u.template block<1, displacement_size / GlobalDim>(
                       i, i * displacement_size / GlobalDim)
                    .noalias() = N;
            }
 
            _local_rhs.noalias() -= n.transpose() * N_u * pressure.dot(N) * w;
        }
 
        auto const indices = NumLib::getIndices(id, dof_table_boundary);
        local_rhs.add(indices, _local_rhs);
    }