ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim > Class Template Reference

Detailed Description

template<typename ShapeFunction, typename IntegrationMethod, int DisplacementDim>
class ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >

Definition at line 34 of file SmallDeformationLocalAssemblerMatrix.h.

#include <SmallDeformationLocalAssemblerMatrix.h>

Inheritance diagram for ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >:

Collaboration diagram for ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >:

Public Types
using	ShapeMatricesType = ShapeMatrixPolicyType< ShapeFunction, DisplacementDim >
using	NodalMatrixType = typename ShapeMatricesType::NodalMatrixType
using	NodalVectorType = typename ShapeMatricesType::NodalVectorType
using	ShapeMatrices = typename ShapeMatricesType::ShapeMatrices
using	BMatricesType = BMatrixPolicyType< ShapeFunction, DisplacementDim >
using	StiffnessMatrixType = typename BMatricesType::StiffnessMatrixType
using	NodalForceVectorType = typename BMatricesType::NodalForceVectorType
using	NodalDisplacementVectorType = typename BMatricesType::NodalForceVectorType

Public Member Functions
	SmallDeformationLocalAssemblerMatrix (SmallDeformationLocalAssemblerMatrix const &)=delete
	SmallDeformationLocalAssemblerMatrix (SmallDeformationLocalAssemblerMatrix &&)=delete
	SmallDeformationLocalAssemblerMatrix (MeshLib::Element const &e, std::size_t const local_matrix_size, bool is_axially_symmetric, unsigned const integration_order, SmallDeformationProcessData< DisplacementDim > &process_data)
void	assemble (double const, double const, std::vector< double > const &, std::vector< double > const &, std::vector< double > &, std::vector< double > &, std::vector< double > &) override
void	assembleWithJacobian (double const t, double const dt, std::vector< double > const &local_x, std::vector< double > const &, const double, const double, std::vector< double > &, std::vector< double > &, std::vector< double > &local_b_data, std::vector< double > &local_Jac_data) override
void	preTimestepConcrete (std::vector< double > const &, double const, double const) override
void	computeSecondaryVariableConcreteWithVector (double const t, Eigen::VectorXd const &local_x) override
Eigen::Map< const Eigen::RowVectorXd >	getShapeMatrix (const unsigned integration_point) const override
	Provides the shape matrix at the given integration point. More...
std::vector< double > const &	getIntPtSigmaXX (const double, std::vector< GlobalVector * > const &, std::vector< NumLib::LocalToGlobalIndexMap const * > const &, std::vector< double > &cache) const override
std::vector< double > const &	getIntPtSigmaYY (const double, std::vector< GlobalVector * > const &, std::vector< NumLib::LocalToGlobalIndexMap const * > const &, std::vector< double > &cache) const override
std::vector< double > const &	getIntPtSigmaZZ (const double, std::vector< GlobalVector * > const &, std::vector< NumLib::LocalToGlobalIndexMap const * > const &, std::vector< double > &cache) const override
std::vector< double > const &	getIntPtSigmaXY (const double, std::vector< GlobalVector * > const &, std::vector< NumLib::LocalToGlobalIndexMap const * > const &, std::vector< double > &cache) const override
std::vector< double > const &	getIntPtSigmaXZ (const double, std::vector< GlobalVector * > const &, std::vector< NumLib::LocalToGlobalIndexMap const * > const &, std::vector< double > &cache) const override
std::vector< double > const &	getIntPtSigmaYZ (const double, std::vector< GlobalVector * > const &, std::vector< NumLib::LocalToGlobalIndexMap const * > const &, std::vector< double > &cache) const override
std::vector< double > const &	getIntPtEpsilonXX (const double, std::vector< GlobalVector * > const &, std::vector< NumLib::LocalToGlobalIndexMap const * > const &, std::vector< double > &cache) const override
std::vector< double > const &	getIntPtEpsilonYY (const double, std::vector< GlobalVector * > const &, std::vector< NumLib::LocalToGlobalIndexMap const * > const &, std::vector< double > &cache) const override
std::vector< double > const &	getIntPtEpsilonZZ (const double, std::vector< GlobalVector * > const &, std::vector< NumLib::LocalToGlobalIndexMap const * > const &, std::vector< double > &cache) const override
std::vector< double > const &	getIntPtEpsilonXY (const double, std::vector< GlobalVector * > const &, std::vector< NumLib::LocalToGlobalIndexMap const * > const &, std::vector< double > &cache) const override
std::vector< double > const &	getIntPtEpsilonXZ (const double, std::vector< GlobalVector * > const &, std::vector< NumLib::LocalToGlobalIndexMap const * > const &, std::vector< double > &cache) const override
std::vector< double > const &	getIntPtEpsilonYZ (const double, std::vector< GlobalVector * > const &, std::vector< NumLib::LocalToGlobalIndexMap const * > const &, std::vector< double > &cache) const override
Public Member Functions inherited from ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerInterface
	SmallDeformationLocalAssemblerInterface ()
	SmallDeformationLocalAssemblerInterface (std::size_t n_local_size, std::vector< unsigned > dofIndex_to_localIndex)
virtual void	assembleWithJacobian (double const, double const, Eigen::VectorXd const &, Eigen::VectorXd &, Eigen::MatrixXd &)
void	computeSecondaryVariableConcrete (double const t, double const, Eigen::VectorXd const &local_x, Eigen::VectorXd const &) override
Public Member Functions inherited from ProcessLib::LocalAssemblerInterface
virtual	~LocalAssemblerInterface ()=default
void	setInitialConditions (std::size_t const mesh_item_id, NumLib::LocalToGlobalIndexMap const &dof_table, GlobalVector const &x, double const t, bool const use_monolithic_scheme, int const process_id)
virtual void	initialize (std::size_t const mesh_item_id, NumLib::LocalToGlobalIndexMap const &dof_table)
virtual void	preAssemble (double const, double const, std::vector< double > const &)
virtual void	assembleForStaggeredScheme (double const t, double const dt, Eigen::VectorXd const &local_x, Eigen::VectorXd const &local_xdot, int const process_id, std::vector< double > &local_M_data, std::vector< double > &local_K_data, std::vector< double > &local_b_data)
virtual void	assembleWithJacobianForStaggeredScheme (double const t, double const dt, Eigen::VectorXd const &local_x, Eigen::VectorXd const &local_xdot, const double dxdot_dx, const double dx_dx, int const process_id, std::vector< double > &local_M_data, std::vector< double > &local_K_data, std::vector< double > &local_b_data, std::vector< double > &local_Jac_data)
virtual void	computeSecondaryVariable (std::size_t const mesh_item_id, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_tables, double const t, double const dt, std::vector< GlobalVector * > const &x, GlobalVector const &x_dot, int const process_id)
virtual void	preTimestep (std::size_t const mesh_item_id, NumLib::LocalToGlobalIndexMap const &dof_table, GlobalVector const &x, double const t, double const delta_t)
virtual void	postTimestep (std::size_t const mesh_item_id, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_tables, std::vector< GlobalVector * > const &x, double const t, double const dt)
void	postNonLinearSolver (std::size_t const mesh_item_id, NumLib::LocalToGlobalIndexMap const &dof_table, GlobalVector const &x, GlobalVector const &xdot, double const t, double const dt, bool const use_monolithic_scheme, int const process_id)
virtual std::vector< double >	interpolateNodalValuesToIntegrationPoints (std::vector< double > const &)
virtual Eigen::Vector3d	getFlux (MathLib::Point3d const &, double const, std::vector< double > const &) const
virtual Eigen::Vector3d	getFlux (MathLib::Point3d const &, double const, std::vector< std::vector< double >> const &) const
	Fits to staggered scheme. More...
Public Member Functions inherited from NumLib::ExtrapolatableElement
virtual	~ExtrapolatableElement ()=default

Private Member Functions
std::vector< double > const &	getIntPtSigma (std::vector< double > &cache, std::size_t const component) const
std::vector< double > const &	getIntPtEpsilon (std::vector< double > &cache, std::size_t const component) const

Private Attributes
SmallDeformationProcessData< DisplacementDim > &	_process_data
std::vector< IntegrationPointDataMatrix< ShapeMatricesType, BMatricesType, DisplacementDim >, Eigen::aligned_allocator< IntegrationPointDataMatrix< ShapeMatricesType, BMatricesType, DisplacementDim > > >	_ip_data
IntegrationMethod	_integration_method
MeshLib::Element const &	_element
bool const	_is_axially_symmetric
SecondaryData< typename ShapeMatrices::ShapeType >	_secondary_data

Additional Inherited Members

Member Typedef Documentation

BMatricesType

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

using ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::BMatricesType = BMatrixPolicyType<ShapeFunction, DisplacementDim>

Definition at line 43 of file SmallDeformationLocalAssemblerMatrix.h.

NodalDisplacementVectorType

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

using ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::NodalDisplacementVectorType = typename BMatricesType::NodalForceVectorType

Definition at line 47 of file SmallDeformationLocalAssemblerMatrix.h.

NodalForceVectorType

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

using ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::NodalForceVectorType = typename BMatricesType::NodalForceVectorType

Definition at line 46 of file SmallDeformationLocalAssemblerMatrix.h.

NodalMatrixType

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

using ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::NodalMatrixType = typename ShapeMatricesType::NodalMatrixType

Definition at line 40 of file SmallDeformationLocalAssemblerMatrix.h.

NodalVectorType

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

using ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::NodalVectorType = typename ShapeMatricesType::NodalVectorType

Definition at line 41 of file SmallDeformationLocalAssemblerMatrix.h.

ShapeMatrices

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

using ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::ShapeMatrices = typename ShapeMatricesType::ShapeMatrices

Definition at line 42 of file SmallDeformationLocalAssemblerMatrix.h.

ShapeMatricesType

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

using ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::ShapeMatricesType = ShapeMatrixPolicyType<ShapeFunction, DisplacementDim>

Definition at line 38 of file SmallDeformationLocalAssemblerMatrix.h.

StiffnessMatrixType

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

using ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::StiffnessMatrixType = typename BMatricesType::StiffnessMatrixType

Definition at line 45 of file SmallDeformationLocalAssemblerMatrix.h.

Constructor & Destructor Documentation

SmallDeformationLocalAssemblerMatrix() [1/3]

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::SmallDeformationLocalAssemblerMatrix ( SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim > const &  )

delete

SmallDeformationLocalAssemblerMatrix() [2/3]

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::SmallDeformationLocalAssemblerMatrix ( SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim > &&  )

delete

SmallDeformationLocalAssemblerMatrix() [3/3]

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::SmallDeformationLocalAssemblerMatrix	(	MeshLib::Element const &	e,
		std::size_t const	local_matrix_size,
		bool	is_axially_symmetric,
		unsigned const	integration_order,
		SmallDeformationProcessData< DisplacementDim > &	process_data
	)

Definition at line 39 of file SmallDeformationLocalAssemblerMatrix-impl.h.

    : _process_data(process_data),
      _integration_method(integration_order),
      _element(e),
      _is_axially_symmetric(is_axially_symmetric)
{
    unsigned const n_integration_points =
        _integration_method.getNumberOfPoints();
 
    _ip_data.reserve(n_integration_points);
    _secondary_data.N.resize(n_integration_points);
 
    auto const shape_matrices =
        NumLib::initShapeMatrices<ShapeFunction, ShapeMatricesType,
                                  DisplacementDim>(e, is_axially_symmetric,
                                                   _integration_method);
 
    auto& solid_material = MaterialLib::Solids::selectSolidConstitutiveRelation(
        _process_data.solid_materials, _process_data.material_ids, e.getID());
 
    for (unsigned ip = 0; ip < n_integration_points; ip++)
    {
        _ip_data.emplace_back(solid_material);
        auto& ip_data = _ip_data[ip];
        auto const& sm = shape_matrices[ip];
        ip_data.N = sm.N;
        ip_data.dNdx = sm.dNdx;
        ip_data.integration_weight =
            _integration_method.getWeightedPoint(ip).getWeight() *
            sm.integralMeasure * sm.detJ;
 
        // Initialize current time step values
        static const int kelvin_vector_size =
            MathLib::KelvinVector::kelvin_vector_dimensions(DisplacementDim);
        ip_data._sigma.setZero(kelvin_vector_size);
        ip_data._eps.setZero(kelvin_vector_size);
 
        // Previous time step values are not initialized and are set later.
        ip_data._sigma_prev.resize(kelvin_vector_size);
        ip_data._eps_prev.resize(kelvin_vector_size);
 
        ip_data._C.resize(kelvin_vector_size, kelvin_vector_size);
 
        _secondary_data.N[ip] = sm.N;
    }
}

References ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::_integration_method, ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::_ip_data, ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::_process_data, ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::_secondary_data, MeshLib::Element::getID(), NumLib::initShapeMatrices(), MathLib::KelvinVector::kelvin_vector_dimensions(), ProcessLib::HeatTransportBHE::SecondaryData< ShapeMatrixType >::N, and MaterialLib::Solids::selectSolidConstitutiveRelation().

Member Function Documentation

assemble()

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

void ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::assemble ( double const  , double const  , std::vector< double > const &  , std::vector< double > const &  , std::vector< double > &  , std::vector< double > &  , std::vector< double > &    )

inlineoverridevirtual

Reimplemented from ProcessLib::LocalAssemblerInterface.

Definition at line 62 of file SmallDeformationLocalAssemblerMatrix.h.

    {
        OGS_FATAL(
            "SmallDeformationLocalAssembler: assembly without jacobian is not "
            "implemented.");
    }

References OGS_FATAL.

assembleWithJacobian()

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

void ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::assembleWithJacobian ( double const  t, double const  dt, std::vector< double > const &  local_x, std::vector< double > const &  , const double  , const double  , std::vector< double > &  , std::vector< double > &  , std::vector< double > &  local_b_data, std::vector< double > &  local_Jac_data  )

overridevirtual

Reimplemented from ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerInterface.

Definition at line 95 of file SmallDeformationLocalAssemblerMatrix-impl.h.

{
    assert(_element.getDimension() == DisplacementDim);
 
    auto const local_matrix_size = local_x.size();
 
    auto local_Jac = MathLib::createZeroedMatrix<StiffnessMatrixType>(
        local_Jac_data, local_matrix_size, local_matrix_size);
 
    auto local_b = MathLib::createZeroedVector<NodalDisplacementVectorType>(
        local_b_data, local_matrix_size);
 
    unsigned const n_integration_points =
        _integration_method.getNumberOfPoints();
 
    MPL::VariableArray variables;
    MPL::VariableArray variables_prev;
    ParameterLib::SpatialPosition x_position;
    x_position.setElementID(_element.getID());
 
    for (unsigned ip = 0; ip < n_integration_points; ip++)
    {
        x_position.setIntegrationPoint(ip);
        auto const& w = _ip_data[ip].integration_weight;
 
        auto const& N = _ip_data[ip].N;
        auto const& dNdx = _ip_data[ip].dNdx;
        auto const x_coord =
            NumLib::interpolateXCoordinate<ShapeFunction, ShapeMatricesType>(
                _element, N);
        auto const B =
            LinearBMatrix::computeBMatrix<DisplacementDim,
                                          ShapeFunction::NPOINTS,
                                          typename BMatricesType::BMatrixType>(
                dNdx, N, x_coord, _is_axially_symmetric);
 
        auto const& eps_prev = _ip_data[ip]._eps_prev;
        auto const& sigma_prev = _ip_data[ip]._sigma_prev;
 
        auto& eps = _ip_data[ip]._eps;
        auto& sigma = _ip_data[ip]._sigma;
        auto& state = _ip_data[ip]._material_state_variables;
 
        eps.noalias() =
            B * Eigen::Map<typename BMatricesType::NodalForceVectorType const>(
                    local_x.data(), ShapeFunction::NPOINTS * DisplacementDim);
 
        variables[static_cast<int>(
                      MaterialPropertyLib::Variable::mechanical_strain)]
            .emplace<MathLib::KelvinVector::KelvinVectorType<DisplacementDim>>(
                eps);
 
        variables_prev[static_cast<int>(MaterialPropertyLib::Variable::stress)]
            .emplace<MathLib::KelvinVector::KelvinVectorType<DisplacementDim>>(
                sigma_prev);
        variables_prev[static_cast<int>(
                           MaterialPropertyLib::Variable::mechanical_strain)]
            .emplace<MathLib::KelvinVector::KelvinVectorType<DisplacementDim>>(
                eps_prev);
        variables_prev[static_cast<int>(
                           MaterialPropertyLib::Variable::temperature)]
            .emplace<double>(_process_data._reference_temperature);
 
        auto&& solution = _ip_data[ip]._solid_material.integrateStress(
            variables_prev, variables, t, x_position, dt, *state);
 
        if (!solution)
        {
            OGS_FATAL("Computation of local constitutive relation failed.");
        }
 
        MathLib::KelvinVector::KelvinMatrixType<DisplacementDim> C;
        std::tie(sigma, state, C) = std::move(*solution);
 
        local_b.noalias() -= B.transpose() * sigma * w;
        local_Jac.noalias() += B.transpose() * C * B * w;
    }
}

References ProcessLib::LinearBMatrix::computeBMatrix(), MaterialPropertyLib::mechanical_strain, OGS_FATAL, ParameterLib::SpatialPosition::setElementID(), ParameterLib::SpatialPosition::setIntegrationPoint(), MaterialPropertyLib::stress, and MaterialPropertyLib::temperature.

computeSecondaryVariableConcreteWithVector()

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

void ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::computeSecondaryVariableConcreteWithVector ( double const  t, Eigen::VectorXd const &  local_x  )

overridevirtual

Implements ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerInterface.

Definition at line 185 of file SmallDeformationLocalAssemblerMatrix-impl.h.

{
    // Compute average value per element
    const int n = DisplacementDim == 2 ? 4 : 6;
    Eigen::VectorXd ele_stress = Eigen::VectorXd::Zero(n);
    Eigen::VectorXd ele_strain = Eigen::VectorXd::Zero(n);
 
    unsigned const n_integration_points =
        _integration_method.getNumberOfPoints();
    for (unsigned ip = 0; ip < n_integration_points; ip++)
    {
        auto& ip_data = _ip_data[ip];
 
        ele_stress += ip_data._sigma;
        ele_strain += ip_data._eps;
    }
    ele_stress /= n_integration_points;
    ele_strain /= n_integration_points;
 
    (*_process_data._mesh_prop_stress_xx)[_element.getID()] = ele_stress[0];
    (*_process_data._mesh_prop_stress_yy)[_element.getID()] = ele_stress[1];
    (*_process_data._mesh_prop_stress_zz)[_element.getID()] = ele_stress[2];
    (*_process_data._mesh_prop_stress_xy)[_element.getID()] = ele_stress[3];
    if (DisplacementDim == 3)
    {
        (*_process_data._mesh_prop_stress_yz)[_element.getID()] = ele_stress[4];
        (*_process_data._mesh_prop_stress_xz)[_element.getID()] = ele_stress[5];
    }
 
    (*_process_data._mesh_prop_strain_xx)[_element.getID()] = ele_strain[0];
    (*_process_data._mesh_prop_strain_yy)[_element.getID()] = ele_strain[1];
    (*_process_data._mesh_prop_strain_zz)[_element.getID()] = ele_strain[2];
    (*_process_data._mesh_prop_strain_xy)[_element.getID()] = ele_strain[3];
    if (DisplacementDim == 3)
    {
        (*_process_data._mesh_prop_strain_yz)[_element.getID()] = ele_strain[4];
        (*_process_data._mesh_prop_strain_xz)[_element.getID()] = ele_strain[5];
    }
}

getIntPtEpsilon()

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

std::vector<double> const& ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::getIntPtEpsilon ( std::vector< double > &  cache, std::size_t const  component  ) const

inlineprivate

Definition at line 242 of file SmallDeformationLocalAssemblerMatrix.h.

    {
        cache.clear();
        cache.reserve(_ip_data.size());
 
        for (auto const& ip_data : _ip_data)
        {
            if (component < 3)
            {  // xx, yy, zz components
                cache.push_back(ip_data._eps[component]);
            }
            else
            {  // mixed xy, yz, xz components
                cache.push_back(ip_data._eps[component] / std::sqrt(2));
            }
        }
 
        return cache;
    }

References ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::_ip_data.

Referenced by ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::getIntPtEpsilonXX(), ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::getIntPtEpsilonXY(), ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::getIntPtEpsilonXZ(), ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::getIntPtEpsilonYY(), ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::getIntPtEpsilonYZ(), and ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::getIntPtEpsilonZZ().

getIntPtEpsilonXX()

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

std::vector<double> const& ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::getIntPtEpsilonXX ( const double  , std::vector< GlobalVector * > const &  , std::vector< NumLib::LocalToGlobalIndexMap const * > const &  , std::vector< double > &  cache  ) const

inlineoverridevirtual

Implements ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerInterface.

Definition at line 164 of file SmallDeformationLocalAssemblerMatrix.h.

    {
        return getIntPtEpsilon(cache, 0);
    }

References ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::getIntPtEpsilon().

getIntPtEpsilonXY()

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

std::vector<double> const& ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::getIntPtEpsilonXY ( const double  , std::vector< GlobalVector * > const &  , std::vector< NumLib::LocalToGlobalIndexMap const * > const &  , std::vector< double > &  cache  ) const

inlineoverridevirtual

Implements ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerInterface.

Definition at line 191 of file SmallDeformationLocalAssemblerMatrix.h.

    {
        return getIntPtEpsilon(cache, 3);
    }

References ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::getIntPtEpsilon().

getIntPtEpsilonXZ()

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

std::vector<double> const& ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::getIntPtEpsilonXZ ( const double  , std::vector< GlobalVector * > const &  , std::vector< NumLib::LocalToGlobalIndexMap const * > const &  , std::vector< double > &  cache  ) const

inlineoverridevirtual

Implements ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerInterface.

Definition at line 200 of file SmallDeformationLocalAssemblerMatrix.h.

    {
        assert(DisplacementDim == 3);
        return getIntPtEpsilon(cache, 4);
    }

References ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::getIntPtEpsilon().

getIntPtEpsilonYY()

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

std::vector<double> const& ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::getIntPtEpsilonYY ( const double  , std::vector< GlobalVector * > const &  , std::vector< NumLib::LocalToGlobalIndexMap const * > const &  , std::vector< double > &  cache  ) const

inlineoverridevirtual

Implements ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerInterface.

Definition at line 173 of file SmallDeformationLocalAssemblerMatrix.h.

    {
        return getIntPtEpsilon(cache, 1);
    }

References ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::getIntPtEpsilon().

getIntPtEpsilonYZ()

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

std::vector<double> const& ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::getIntPtEpsilonYZ ( const double  , std::vector< GlobalVector * > const &  , std::vector< NumLib::LocalToGlobalIndexMap const * > const &  , std::vector< double > &  cache  ) const

inlineoverridevirtual

Implements ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerInterface.

Definition at line 210 of file SmallDeformationLocalAssemblerMatrix.h.

    {
        assert(DisplacementDim == 3);
        return getIntPtEpsilon(cache, 5);
    }

References ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::getIntPtEpsilon().

getIntPtEpsilonZZ()

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

std::vector<double> const& ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::getIntPtEpsilonZZ ( const double  , std::vector< GlobalVector * > const &  , std::vector< NumLib::LocalToGlobalIndexMap const * > const &  , std::vector< double > &  cache  ) const

inlineoverridevirtual

Implements ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerInterface.

Definition at line 182 of file SmallDeformationLocalAssemblerMatrix.h.

    {
        return getIntPtEpsilon(cache, 2);
    }

References ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::getIntPtEpsilon().

getIntPtSigma()

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

std::vector<double> const& ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::getIntPtSigma ( std::vector< double > &  cache, std::size_t const  component  ) const

inlineprivate

Definition at line 221 of file SmallDeformationLocalAssemblerMatrix.h.

    {
        cache.clear();
        cache.reserve(_ip_data.size());
 
        for (auto const& ip_data : _ip_data)
        {
            if (component < 3)
            {  // xx, yy, zz components
                cache.push_back(ip_data._sigma[component]);
            }
            else
            {  // mixed xy, yz, xz components
                cache.push_back(ip_data._sigma[component] / std::sqrt(2));
            }
        }
 
        return cache;
    }

References ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::_ip_data.

Referenced by ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::getIntPtSigmaXX(), ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::getIntPtSigmaXY(), ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::getIntPtSigmaXZ(), ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::getIntPtSigmaYY(), ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::getIntPtSigmaYZ(), and ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::getIntPtSigmaZZ().

getIntPtSigmaXX()

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

std::vector<double> const& ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::getIntPtSigmaXX ( const double  , std::vector< GlobalVector * > const &  , std::vector< NumLib::LocalToGlobalIndexMap const * > const &  , std::vector< double > &  cache  ) const

inlineoverridevirtual

Implements ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerInterface.

Definition at line 108 of file SmallDeformationLocalAssemblerMatrix.h.

    {
        return getIntPtSigma(cache, 0);
    }

References ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::getIntPtSigma().

getIntPtSigmaXY()

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

std::vector<double> const& ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::getIntPtSigmaXY ( const double  , std::vector< GlobalVector * > const &  , std::vector< NumLib::LocalToGlobalIndexMap const * > const &  , std::vector< double > &  cache  ) const

inlineoverridevirtual

Implements ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerInterface.

Definition at line 135 of file SmallDeformationLocalAssemblerMatrix.h.

    {
        return getIntPtSigma(cache, 3);
    }

References ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::getIntPtSigma().

getIntPtSigmaXZ()

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

std::vector<double> const& ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::getIntPtSigmaXZ ( const double  , std::vector< GlobalVector * > const &  , std::vector< NumLib::LocalToGlobalIndexMap const * > const &  , std::vector< double > &  cache  ) const

inlineoverridevirtual

Implements ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerInterface.

Definition at line 144 of file SmallDeformationLocalAssemblerMatrix.h.

    {
        assert(DisplacementDim == 3);
        return getIntPtSigma(cache, 4);
    }

References ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::getIntPtSigma().

getIntPtSigmaYY()

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

std::vector<double> const& ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::getIntPtSigmaYY ( const double  , std::vector< GlobalVector * > const &  , std::vector< NumLib::LocalToGlobalIndexMap const * > const &  , std::vector< double > &  cache  ) const

inlineoverridevirtual

Implements ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerInterface.

Definition at line 117 of file SmallDeformationLocalAssemblerMatrix.h.

    {
        return getIntPtSigma(cache, 1);
    }

References ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::getIntPtSigma().

getIntPtSigmaYZ()

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

std::vector<double> const& ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::getIntPtSigmaYZ ( const double  , std::vector< GlobalVector * > const &  , std::vector< NumLib::LocalToGlobalIndexMap const * > const &  , std::vector< double > &  cache  ) const

inlineoverridevirtual

Implements ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerInterface.

Definition at line 154 of file SmallDeformationLocalAssemblerMatrix.h.

    {
        assert(DisplacementDim == 3);
        return getIntPtSigma(cache, 5);
    }

References ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::getIntPtSigma().

getIntPtSigmaZZ()

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

std::vector<double> const& ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::getIntPtSigmaZZ ( const double  , std::vector< GlobalVector * > const &  , std::vector< NumLib::LocalToGlobalIndexMap const * > const &  , std::vector< double > &  cache  ) const

inlineoverridevirtual

Implements ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerInterface.

Definition at line 126 of file SmallDeformationLocalAssemblerMatrix.h.

    {
        return getIntPtSigma(cache, 2);
    }

References ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::getIntPtSigma().

getShapeMatrix()

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

Eigen::Map<const Eigen::RowVectorXd> ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::getShapeMatrix ( const unsigned  integration_point ) const

inlineoverridevirtual

Provides the shape matrix at the given integration point.

Implements NumLib::ExtrapolatableElement.

Definition at line 99 of file SmallDeformationLocalAssemblerMatrix.h.

    {
        auto const& N = _secondary_data.N[integration_point];
 
        // assumes N is stored contiguously in memory
        return Eigen::Map<const Eigen::RowVectorXd>(N.data(), N.size());
    }

References ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::_secondary_data, and ProcessLib::HeatTransportBHE::SecondaryData< ShapeMatrixType >::N.

preTimestepConcrete()

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

void ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::preTimestepConcrete ( std::vector< double > const &  , double const  , double const    )

inlineoverridevirtual

Reimplemented from ProcessLib::LocalAssemblerInterface.

Definition at line 83 of file SmallDeformationLocalAssemblerMatrix.h.

    {
        unsigned const n_integration_points =
            _integration_method.getNumberOfPoints();
 
        for (unsigned ip = 0; ip < n_integration_points; ip++)
        {
            _ip_data[ip].pushBackState();
        }
    }

References ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::_integration_method, and ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::_ip_data.

Member Data Documentation

_element

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

MeshLib::Element const& ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::_element

private

Definition at line 272 of file SmallDeformationLocalAssemblerMatrix.h.

_integration_method

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

IntegrationMethod ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::_integration_method

private

Definition at line 271 of file SmallDeformationLocalAssemblerMatrix.h.

Referenced by ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::SmallDeformationLocalAssemblerMatrix(), and ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::preTimestepConcrete().

_ip_data

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

std::vector<IntegrationPointDataMatrix<ShapeMatricesType, BMatricesType, DisplacementDim>, Eigen::aligned_allocator<IntegrationPointDataMatrix< ShapeMatricesType, BMatricesType, DisplacementDim> > > ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::_ip_data

private

Definition at line 269 of file SmallDeformationLocalAssemblerMatrix.h.

Referenced by ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::SmallDeformationLocalAssemblerMatrix(), ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::getIntPtEpsilon(), ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::getIntPtSigma(), and ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::preTimestepConcrete().

_is_axially_symmetric

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

bool const ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::_is_axially_symmetric

private

Definition at line 273 of file SmallDeformationLocalAssemblerMatrix.h.

_process_data

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

SmallDeformationProcessData<DisplacementDim>& ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::_process_data

private

Definition at line 263 of file SmallDeformationLocalAssemblerMatrix.h.

Referenced by ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::SmallDeformationLocalAssemblerMatrix().

_secondary_data

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

SecondaryData<typename ShapeMatrices::ShapeType> ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::_secondary_data

private

Definition at line 274 of file SmallDeformationLocalAssemblerMatrix.h.

Referenced by ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::SmallDeformationLocalAssemblerMatrix(), and ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrix< ShapeFunction, IntegrationMethod, DisplacementDim >::getShapeMatrix().

---

The documentation for this class was generated from the following files:

• ProcessLib/LIE/SmallDeformation/LocalAssembler/SmallDeformationLocalAssemblerMatrix.h
• ProcessLib/LIE/SmallDeformation/LocalAssembler/SmallDeformationLocalAssemblerMatrix-impl.h