Detailed Description

template<int DisplacementDim, typename ConstitutiveTraits>
struct ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >

Definition at line 27 of file LocalAssemblerInterface.h.

#include <LocalAssemblerInterface.h>

Inheritance diagram for ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >:

[legend]

Collaboration diagram for ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >:

[legend]

Public Member Functions
	LocalAssemblerInterface (MeshLib::Element const &e, NumLib::GenericIntegrationMethod const &integration_method, bool const is_axially_symmetric, ThermoRichardsMechanicsProcessData< DisplacementDim, ConstitutiveTraits > &process_data)
std::size_t	setIPDataInitialConditions (std::string_view name, double const *values, int const integration_order)
std::vector< double >	getMaterialStateVariableInternalState (std::function< std::span< double >(typename MaterialLib::Solids::MechanicsBase< DisplacementDim >::MaterialStateVariables &)> const &get_values_span, int const &n_components) const
unsigned	getNumberOfIntegrationPoints () const
int	getMaterialID () const
MaterialLib::Solids::MechanicsBase< DisplacementDim >::MaterialStateVariables const &	getMaterialStateVariablesAt (unsigned integration_point) const
void	postTimestepConcrete (Eigen::VectorXd const &, Eigen::VectorXd const &, double const, double const, int const) override
Public Member Functions inherited from ProcessLib::LocalAssemblerInterface
virtual	~LocalAssemblerInterface ()=default
virtual void	setInitialConditions (std::size_t const mesh_item_id, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_tables, std::vector< GlobalVector * > const &x, double const t, 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	assemble (double const t, double const dt, std::vector< double > const &local_x, std::vector< double > const &local_x_prev, std::vector< double > &local_M_data, std::vector< double > &local_K_data, std::vector< double > &local_b_data)
virtual void	assembleForStaggeredScheme (double const t, double const dt, Eigen::VectorXd const &local_x, Eigen::VectorXd const &local_x_prev, int const process_id, std::vector< double > &local_M_data, std::vector< double > &local_K_data, std::vector< double > &local_b_data)
virtual void	assembleWithJacobian (double const t, double const dt, std::vector< double > const &local_x, std::vector< double > const &local_x_prev, std::vector< double > &local_b_data, std::vector< double > &local_Jac_data)
virtual void	assembleWithJacobianForStaggeredScheme (double const t, double const dt, Eigen::VectorXd const &local_x, Eigen::VectorXd const &local_x_prev, int const process_id, 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_prev, 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, std::vector< GlobalVector * > const &x_prev, double const t, double const dt, int const process_id)
void	postNonLinearSolver (std::size_t const mesh_item_id, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_tables, std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &x_prev, double const t, double const dt, int const process_id)
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.
virtual std::optional< VectorSegment >	getVectorDeformationSegment () const
Public Member Functions inherited from NumLib::ExtrapolatableElement
virtual Eigen::Map< const Eigen::RowVectorXd >	getShapeMatrix (const unsigned integration_point) const =0
	Provides the shape matrix at the given integration point.
virtual	~ExtrapolatableElement ()=default

Static Public Member Functions
static auto	getReflectionDataForOutput ()

Protected Attributes
ThermoRichardsMechanicsProcessData< DisplacementDim, ConstitutiveTraits > &	process_data_
std::vector< typename ConstitutiveTraits::StatefulData >	current_states_
std::vector< typename ConstitutiveTraits::StatefulDataPrev >	prev_states_
std::vector< MaterialStateData< DisplacementDim > >	material_states_
NumLib::GenericIntegrationMethod const &	integration_method_
MeshLib::Element const &	element_
bool const	is_axially_symmetric_
ConstitutiveTraits::SolidConstitutiveRelation const &	solid_material_
std::vector< typename ConstitutiveTraits::OutputData >	output_data_

Constructor & Destructor Documentation

◆ LocalAssemblerInterface()

template<int DisplacementDim, typename ConstitutiveTraits>

ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::LocalAssemblerInterface	(	MeshLib::Element const &	e,
		NumLib::GenericIntegrationMethod const &	integration_method,
		bool const	is_axially_symmetric,
		ThermoRichardsMechanicsProcessData< DisplacementDim, ConstitutiveTraits > &	process_data )

inline

Definition at line 30 of file LocalAssemblerInterface.h.

        : process_data_(process_data),
          integration_method_(integration_method),
          element_(e),
          is_axially_symmetric_(is_axially_symmetric),
          solid_material_(MaterialLib::Solids::selectSolidConstitutiveRelation(
              process_data_.solid_materials, process_data_.material_ids,
              e.getID()))
    {
        unsigned const n_integration_points =
            integration_method_.getNumberOfPoints();
 
        current_states_.resize(n_integration_points);
        prev_states_.resize(n_integration_points);
        output_data_.resize(n_integration_points);
 
        material_states_.reserve(n_integration_points);
        for (unsigned ip = 0; ip < n_integration_points; ++ip)
        {
            material_states_.emplace_back(
                solid_material_.createMaterialStateVariables());
 
            // Set initial strain field to zero.
            std::get<StrainData<DisplacementDim>>(current_states_[ip]).eps =
                KelvinVector<DisplacementDim>::Zero();
        }
    }

References current_states_, element_, integration_method_, is_axially_symmetric_, material_states_, output_data_, prev_states_, process_data_, and solid_material_.

Referenced by ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::ThermoRichardsMechanicsLocalAssembler(), and getReflectionDataForOutput().

Member Function Documentation

◆ getMaterialID()

template<int DisplacementDim, typename ConstitutiveTraits>

int ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::getMaterialID ( ) const

inline

Definition at line 136 of file LocalAssemblerInterface.h.

    {
        return process_data_.material_ids == nullptr
                   ? 0
                   : (*process_data_.material_ids)[element_.getID()];
    }

References element_, and process_data_.

◆ getMaterialStateVariableInternalState()

template<int DisplacementDim, typename ConstitutiveTraits>

std::vector< double > ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::getMaterialStateVariableInternalState	(	std::function< std::span< double >(typename MaterialLib::Solids::MechanicsBase< DisplacementDim >::MaterialStateVariables &)> const &	get_values_span,
		int const &	n_components ) const

inline

Definition at line 118 of file LocalAssemblerInterface.h.

    {
        return ProcessLib::getIntegrationPointDataMaterialStateVariables(
            material_states_,
            &MaterialStateData<DisplacementDim>::material_state_variables,
            get_values_span, n_components);
    }

References ProcessLib::getIntegrationPointDataMaterialStateVariables(), ProcessLib::ThermoRichardsMechanics::MaterialStateData< DisplacementDim >::material_state_variables, and material_states_.

◆ getMaterialStateVariablesAt()

template<int DisplacementDim, typename ConstitutiveTraits>

MaterialLib::Solids::MechanicsBase< DisplacementDim >::MaterialStateVariables const & ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::getMaterialStateVariablesAt ( unsigned integration_point ) const

inline

Definition at line 145 of file LocalAssemblerInterface.h.

    {
        return *material_states_[integration_point].material_state_variables;
    }

References material_states_.

◆ getNumberOfIntegrationPoints()

template<int DisplacementDim, typename ConstitutiveTraits>

unsigned ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::getNumberOfIntegrationPoints ( ) const

inline

Definition at line 131 of file LocalAssemblerInterface.h.

    {
        return integration_method_.getNumberOfPoints();
    }

References integration_method_.

◆ getReflectionDataForOutput()

template<int DisplacementDim, typename ConstitutiveTraits>

auto ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::getReflectionDataForOutput ( )

inlinestatic

Definition at line 170 of file LocalAssemblerInterface.h.

    {
        using Self =
            LocalAssemblerInterface<DisplacementDim, ConstitutiveTraits>;
 
        return ProcessLib::Reflection::reflectWithoutName(
            &Self::current_states_, &Self::output_data_);
    }

References LocalAssemblerInterface(), and ProcessLib::Reflection::reflectWithoutName().

Referenced by ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsProcess< DisplacementDim, ConstitutiveTraits >::initializeConcreteProcess().

◆ postTimestepConcrete()

template<int DisplacementDim, typename ConstitutiveTraits>

void ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::postTimestepConcrete	(	Eigen::VectorXd const &	,
		Eigen::VectorXd const &	,
		double const	,
		double const	,
		int const	)

inlineoverridevirtual

Reimplemented from ProcessLib::LocalAssemblerInterface.

Definition at line 150 of file LocalAssemblerInterface.h.

    {
        unsigned const n_integration_points =
            integration_method_.getNumberOfPoints();
 
        for (unsigned ip = 0; ip < n_integration_points; ip++)
        {
            // TODO re-evaluate part of the assembly in order to be consistent?
            material_states_[ip].pushBackState();
        }
 
        for (unsigned ip = 0; ip < n_integration_points; ip++)
        {
            prev_states_[ip] = current_states_[ip];
        }
    }

References current_states_, integration_method_, material_states_, and prev_states_.

◆ setIPDataInitialConditions()

template<int DisplacementDim, typename ConstitutiveTraits>

std::size_t ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::setIPDataInitialConditions	(	std::string_view	name,
		double const *	values,
		int const	integration_order )

inline

Definition at line 63 of file LocalAssemblerInterface.h.

    {
        if (integration_order !=
            static_cast<int>(integration_method_.getIntegrationOrder()))
        {
            OGS_FATAL(
                "Setting integration point initial conditions; The integration "
                "order of the local assembler for element {:d} is different "
                "from the integration order in the initial condition.",
                element_.getID());
        }
 
        if (name == "sigma" && process_data_.initial_stress.value)
        {
            OGS_FATAL(
                "Setting initial conditions for stress from integration "
                "point data and from a parameter '{:s}' is not possible "
                "simultaneously.",
                process_data_.initial_stress.value->name);
        }
 
        // TODO (naumov) this information is runtime information and I'm not
        // sure how to put it into the reflected data structure. The
        // reflectWithName function also supports only a single return value.
        if (name.starts_with("material_state_variable_"))
        {
            name.remove_prefix(24);
 
            auto const& internal_variables =
                solid_material_.getInternalVariables();
            if (auto const iv = std::find_if(
                    begin(internal_variables), end(internal_variables),
                    [&name](auto const& iv) { return iv.name == name; });
                iv != end(internal_variables))
            {
                DBUG("Setting material state variable '{:s}'", name);
                return ProcessLib::
                    setIntegrationPointDataMaterialStateVariables(
                        values, material_states_,
                        &MaterialStateData<
                            DisplacementDim>::material_state_variables,
                        iv->reference);
            }
            return 0;
        }
 
        // TODO this logic could be pulled out of the local assembler into the
        // process. That might lead to a slightly better performance due to less
        // string comparisons.
        return ProcessLib::Reflection::reflectSetIPData<DisplacementDim>(
            name, values, current_states_);
    }

References current_states_, DBUG(), element_, integration_method_, material_states_, OGS_FATAL, process_data_, ProcessLib::Reflection::reflectSetIPData(), ProcessLib::setIntegrationPointDataMaterialStateVariables(), and solid_material_.

Member Data Documentation

◆ current_states_

template<int DisplacementDim, typename ConstitutiveTraits>

std::vector<typename ConstitutiveTraits::StatefulData> ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::current_states_

protected

Definition at line 184 of file LocalAssemblerInterface.h.

◆ element_

template<int DisplacementDim, typename ConstitutiveTraits>

MeshLib::Element const& ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::element_

protected

Definition at line 193 of file LocalAssemblerInterface.h.

◆ integration_method_

template<int DisplacementDim, typename ConstitutiveTraits>

NumLib::GenericIntegrationMethod const& ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::integration_method_

protected

Definition at line 192 of file LocalAssemblerInterface.h.

◆ is_axially_symmetric_

template<int DisplacementDim, typename ConstitutiveTraits>

bool const ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::is_axially_symmetric_

protected

Definition at line 194 of file LocalAssemblerInterface.h.

Referenced by LocalAssemblerInterface(), ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::assembleWithJacobianSingleIP(), and ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::computeSecondaryVariableConcrete().

◆ material_states_

template<int DisplacementDim, typename ConstitutiveTraits>

std::vector<MaterialStateData<DisplacementDim> > ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::material_states_

protected

Definition at line 190 of file LocalAssemblerInterface.h.

Referenced by LocalAssemblerInterface(), ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::assembleWithJacobian(), ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::computeSecondaryVariableConcrete(), getMaterialStateVariableInternalState(), getMaterialStateVariablesAt(), ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::initializeConcrete(), postTimestepConcrete(), and setIPDataInitialConditions().

◆ output_data_

template<int DisplacementDim, typename ConstitutiveTraits>

std::vector<typename ConstitutiveTraits::OutputData> ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::output_data_

protected

Definition at line 199 of file LocalAssemblerInterface.h.

Referenced by LocalAssemblerInterface(), ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::assembleWithJacobian(), and ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::computeSecondaryVariableConcrete().

◆ prev_states_

template<int DisplacementDim, typename ConstitutiveTraits>

std::vector<typename ConstitutiveTraits::StatefulDataPrev> ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::prev_states_

protected

Definition at line 186 of file LocalAssemblerInterface.h.

◆ process_data_

template<int DisplacementDim, typename ConstitutiveTraits>

ThermoRichardsMechanicsProcessData<DisplacementDim, ConstitutiveTraits>& ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::process_data_

protected

Definition at line 181 of file LocalAssemblerInterface.h.

◆ solid_material_

template<int DisplacementDim, typename ConstitutiveTraits>

ConstitutiveTraits::SolidConstitutiveRelation const& ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::solid_material_

protected

Definition at line 197 of file LocalAssemblerInterface.h.

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

ProcessLib/ThermoRichardsMechanics/CreateThermoRichardsMechanicsLocalAssemblers.h
ProcessLib/ThermoRichardsMechanics/LocalAssemblerInterface.h

Detailed Description

Public Member Functions

Static Public Member Functions

Protected Attributes

Constructor & Destructor Documentation

◆ LocalAssemblerInterface()

Member Function Documentation

◆ getMaterialID()

◆ getMaterialStateVariableInternalState()

◆ getMaterialStateVariablesAt()

◆ getNumberOfIntegrationPoints()

◆ getReflectionDataForOutput()

◆ postTimestepConcrete()

◆ setIPDataInitialConditions()

Member Data Documentation

◆ current_states_

◆ element_

◆ integration_method_

◆ is_axially_symmetric_

◆ material_states_

◆ output_data_

◆ prev_states_

◆ process_data_

◆ solid_material_