Detailed Description

template<typename ShapeFunction, int DisplacementDim>
class ProcessLib::LargeDeformation::LargeDeformationLocalAssembler< ShapeFunction, DisplacementDim >

Definition at line 74 of file LargeDeformationFEM.h.

#include <LargeDeformationFEM.h>

Inheritance diagram for ProcessLib::LargeDeformation::LargeDeformationLocalAssembler< ShapeFunction, DisplacementDim >:

[legend]

Collaboration diagram for ProcessLib::LargeDeformation::LargeDeformationLocalAssembler< ShapeFunction, DisplacementDim >:

[legend]

Public Types
using	ShapeMatricesType
using	NodalMatrixType = typename ShapeMatricesType::NodalMatrixType
using	GlobalDimNodalMatrixType
using	NodalVectorType = typename ShapeMatricesType::NodalVectorType
using	ShapeMatrices = typename ShapeMatricesType::ShapeMatrices
using	BMatricesType = BMatrixPolicyType<ShapeFunction, DisplacementDim>
using	BMatrixType = typename BMatricesType::BMatrixType
using	StiffnessMatrixType = typename BMatricesType::StiffnessMatrixType
using	NodalForceVectorType = typename BMatricesType::NodalForceVectorType
using	NodalDisplacementVectorType
using	GMatricesType = GMatrixPolicyType<ShapeFunction, DisplacementDim>
using	GradientVectorType = typename GMatricesType::GradientVectorType
using	GradientMatrixType = typename GMatricesType::GradientMatrixType
using	VectorTypeForFbar = typename BMatricesType::NodalForceVectorType
using	IpData

Public Member Functions
	LargeDeformationLocalAssembler (LargeDeformationLocalAssembler const &)=delete
	LargeDeformationLocalAssembler (LargeDeformationLocalAssembler &&)=delete
	LargeDeformationLocalAssembler (MeshLib::Element const &e, std::size_t const, NumLib::GenericIntegrationMethod const &integration_method, bool const is_axially_symmetric, LargeDeformationProcessData< DisplacementDim > &process_data)
void	initializeConcrete () override
double	computeOutputStrainData (double const detF0, BMatrixType const &B, GradientVectorType const &grad_u, Eigen::Ref< Eigen::VectorXd const > const &u, typename ConstitutiveRelations::OutputData< DisplacementDim > &output_data) const
ConstitutiveRelations::ConstitutiveData< DisplacementDim >	updateConstitutiveRelations (double const alpha, GradientMatrixType const &G, Eigen::Ref< Eigen::VectorXd const > const &u_prev, ParameterLib::SpatialPosition const &x_position, double const t, double const dt, typename ConstitutiveRelations::ConstitutiveSetting< DisplacementDim > &CS, MaterialPropertyLib::Medium const &medium, typename ConstitutiveRelations::StatefulData< DisplacementDim > &current_state, typename ConstitutiveRelations::StatefulDataPrev< DisplacementDim > const &prev_state, MaterialStateData< DisplacementDim > &material_state, typename ConstitutiveRelations::OutputData< DisplacementDim > &output_data) const
void	assemble (double const, double const, std::vector< double > const &, std::vector< double > const &, std::vector< double > &, std::vector< double > &, std::vector< double > &) override
std::optional< std::tuple< double, VectorTypeForFbar > >	computeFBarVariables (bool const compute_detF0_only, Eigen::VectorXd const &local_x) const
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) override
void	postTimestepConcrete (Eigen::VectorXd const &local_x, Eigen::VectorXd const &local_x_prev, double const t, double const dt, int const) override
Eigen::Map< const Eigen::RowVectorXd >	getShapeMatrix (const unsigned integration_point) const override
	Provides the shape matrix at the given integration point.
Public Member Functions inherited from ProcessLib::LargeDeformation::LargeDeformationLocalAssemblerInterface< DisplacementDim >
	LargeDeformationLocalAssemblerInterface (MeshLib::Element const &e, NumLib::GenericIntegrationMethod const &integration_method, bool const is_axially_symmetric, LargeDeformationProcessData< DisplacementDim > &process_data)
std::size_t	setIPDataInitialConditions (std::string_view name, double const *values, int const integration_order)
	Returns number of read integration points.
unsigned	getNumberOfIntegrationPoints () const
int	getMaterialID () const
std::vector< double >	getMaterialStateVariableInternalState (std::function< std::span< double >(typename MaterialLib::Solids::MechanicsBase< DisplacementDim >::MaterialStateVariables &)> const &get_values_span, int const &n_components) const
MaterialLib::Solids::MechanicsBase< DisplacementDim >::MaterialStateVariables const &	getMaterialStateVariablesAt (unsigned integration_point) const
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	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	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 int	getNumberOfVectorElementsForDeformation () const
Public Member Functions inherited from NumLib::ExtrapolatableElement
virtual	~ExtrapolatableElement ()=default

Static Public Attributes
static constexpr auto &	N_u_op

Static Private Member Functions
static constexpr auto	localDOF (std::vector< double > const &x)

Private Attributes
std::vector< IpData, Eigen::aligned_allocator< IpData > >	_ip_data
SecondaryData< typename ShapeMatrices::ShapeType >	_secondary_data

Static Private Attributes
static const int	displacement_size

Additional Inherited Members
Static Public Member Functions inherited from ProcessLib::LargeDeformation::LargeDeformationLocalAssemblerInterface< DisplacementDim >
static auto	getReflectionDataForOutput ()
Protected Attributes inherited from ProcessLib::LargeDeformation::LargeDeformationLocalAssemblerInterface< DisplacementDim >
LargeDeformationProcessData< DisplacementDim > &	process_data_
std::vector< MaterialStateData< DisplacementDim > >	material_states_
std::vector< typename ConstitutiveRelations::StatefulData< DisplacementDim > >	current_states_
std::vector< typename ConstitutiveRelations::StatefulDataPrev< DisplacementDim > >	prev_states_
std::vector< typename ConstitutiveRelations::OutputData< DisplacementDim > >	output_data_
NumLib::GenericIntegrationMethod const &	integration_method_
MeshLib::Element const &	element_
bool const	is_axially_symmetric_
ConstitutiveRelations::SolidConstitutiveRelation< DisplacementDim > const &	solid_material_

Member Typedef Documentation

◆ BMatricesType

template<typename ShapeFunction, int DisplacementDim>

using ProcessLib::LargeDeformation::LargeDeformationLocalAssembler< ShapeFunction, DisplacementDim >::BMatricesType = BMatrixPolicyType<ShapeFunction, DisplacementDim>

Definition at line 85 of file LargeDeformationFEM.h.

◆ BMatrixType

template<typename ShapeFunction, int DisplacementDim>

using ProcessLib::LargeDeformation::LargeDeformationLocalAssembler< ShapeFunction, DisplacementDim >::BMatrixType = typename BMatricesType::BMatrixType

Definition at line 87 of file LargeDeformationFEM.h.

◆ GlobalDimNodalMatrixType

template<typename ShapeFunction, int DisplacementDim>

using ProcessLib::LargeDeformation::LargeDeformationLocalAssembler< ShapeFunction, DisplacementDim >::GlobalDimNodalMatrixType

Initial value:

typename ShapeMatricesType::GlobalDimNodalMatrixType

EigenFixedShapeMatrixPolicy::GlobalDimNodalMatrixType

MatrixType< GlobalDim, ShapeFunction::NPOINTS > GlobalDimNodalMatrixType

Definition ShapeMatrixPolicy.h:75

Definition at line 81 of file LargeDeformationFEM.h.

◆ GMatricesType

template<typename ShapeFunction, int DisplacementDim>

using ProcessLib::LargeDeformation::LargeDeformationLocalAssembler< ShapeFunction, DisplacementDim >::GMatricesType = GMatrixPolicyType<ShapeFunction, DisplacementDim>

Definition at line 93 of file LargeDeformationFEM.h.

◆ GradientMatrixType

template<typename ShapeFunction, int DisplacementDim>

using ProcessLib::LargeDeformation::LargeDeformationLocalAssembler< ShapeFunction, DisplacementDim >::GradientMatrixType = typename GMatricesType::GradientMatrixType

Definition at line 95 of file LargeDeformationFEM.h.

◆ GradientVectorType

template<typename ShapeFunction, int DisplacementDim>

using ProcessLib::LargeDeformation::LargeDeformationLocalAssembler< ShapeFunction, DisplacementDim >::GradientVectorType = typename GMatricesType::GradientVectorType

Definition at line 94 of file LargeDeformationFEM.h.

◆ IpData

template<typename ShapeFunction, int DisplacementDim>

using ProcessLib::LargeDeformation::LargeDeformationLocalAssembler< ShapeFunction, DisplacementDim >::IpData

Initial value:

IntegrationPointData<BMatricesType, ShapeMatricesType, DisplacementDim>

ProcessLib::HMPhaseField::IntegrationPointData

Definition HMPhaseFieldFEM.h:37

Definition at line 99 of file LargeDeformationFEM.h.

◆ NodalDisplacementVectorType

template<typename ShapeFunction, int DisplacementDim>

using ProcessLib::LargeDeformation::LargeDeformationLocalAssembler< ShapeFunction, DisplacementDim >::NodalDisplacementVectorType

Initial value:

typename BMatricesType::NodalForceVectorType

ProcessLib::BMatrixPolicyType::NodalForceVectorType

VectorType< _number_of_dof > NodalForceVectorType

Rhs residual.

Definition BMatrixPolicy.h:37

Definition at line 90 of file LargeDeformationFEM.h.

◆ NodalForceVectorType

template<typename ShapeFunction, int DisplacementDim>

using ProcessLib::LargeDeformation::LargeDeformationLocalAssembler< ShapeFunction, DisplacementDim >::NodalForceVectorType = typename BMatricesType::NodalForceVectorType

Definition at line 89 of file LargeDeformationFEM.h.

◆ NodalMatrixType

template<typename ShapeFunction, int DisplacementDim>

using ProcessLib::LargeDeformation::LargeDeformationLocalAssembler< ShapeFunction, DisplacementDim >::NodalMatrixType = typename ShapeMatricesType::NodalMatrixType

Definition at line 80 of file LargeDeformationFEM.h.

◆ NodalVectorType

template<typename ShapeFunction, int DisplacementDim>

using ProcessLib::LargeDeformation::LargeDeformationLocalAssembler< ShapeFunction, DisplacementDim >::NodalVectorType = typename ShapeMatricesType::NodalVectorType

Definition at line 83 of file LargeDeformationFEM.h.

◆ ShapeMatrices

template<typename ShapeFunction, int DisplacementDim>

using ProcessLib::LargeDeformation::LargeDeformationLocalAssembler< ShapeFunction, DisplacementDim >::ShapeMatrices = typename ShapeMatricesType::ShapeMatrices

Definition at line 84 of file LargeDeformationFEM.h.

◆ ShapeMatricesType

template<typename ShapeFunction, int DisplacementDim>

using ProcessLib::LargeDeformation::LargeDeformationLocalAssembler< ShapeFunction, DisplacementDim >::ShapeMatricesType

Initial value:

ShapeMatrixPolicyType<ShapeFunction, DisplacementDim>

ShapeMatrixPolicyType

EigenFixedShapeMatrixPolicy< ShapeFunction, GlobalDim > ShapeMatrixPolicyType

Definition ShapeMatrixPolicy.h:134

Definition at line 78 of file LargeDeformationFEM.h.

◆ StiffnessMatrixType

template<typename ShapeFunction, int DisplacementDim>

using ProcessLib::LargeDeformation::LargeDeformationLocalAssembler< ShapeFunction, DisplacementDim >::StiffnessMatrixType = typename BMatricesType::StiffnessMatrixType

Definition at line 88 of file LargeDeformationFEM.h.

◆ VectorTypeForFbar

template<typename ShapeFunction, int DisplacementDim>

using ProcessLib::LargeDeformation::LargeDeformationLocalAssembler< ShapeFunction, DisplacementDim >::VectorTypeForFbar = typename BMatricesType::NodalForceVectorType

Definition at line 97 of file LargeDeformationFEM.h.

Constructor & Destructor Documentation

◆ LargeDeformationLocalAssembler() [1/3]

template<typename ShapeFunction, int DisplacementDim>

ProcessLib::LargeDeformation::LargeDeformationLocalAssembler< ShapeFunction, DisplacementDim >::LargeDeformationLocalAssembler ( LargeDeformationLocalAssembler< ShapeFunction, DisplacementDim > const & )

delete

References LargeDeformationLocalAssembler().

Referenced by LargeDeformationLocalAssembler(), and LargeDeformationLocalAssembler().

◆ LargeDeformationLocalAssembler() [2/3]

template<typename ShapeFunction, int DisplacementDim>

ProcessLib::LargeDeformation::LargeDeformationLocalAssembler< ShapeFunction, DisplacementDim >::LargeDeformationLocalAssembler ( LargeDeformationLocalAssembler< ShapeFunction, DisplacementDim > && )

delete

References LargeDeformationLocalAssembler().

◆ LargeDeformationLocalAssembler() [3/3]

template<typename ShapeFunction, int DisplacementDim>

ProcessLib::LargeDeformation::LargeDeformationLocalAssembler< ShapeFunction, DisplacementDim >::LargeDeformationLocalAssembler	(	MeshLib::Element const &	e,
		std::size_t const	,
		NumLib::GenericIntegrationMethod const &	integration_method,
		bool const	is_axially_symmetric,
		LargeDeformationProcessData< DisplacementDim > &	process_data )

inline

Definition at line 109 of file LargeDeformationFEM.h.

        : LargeDeformationLocalAssemblerInterface<DisplacementDim>(
              e, integration_method, is_axially_symmetric, process_data)
    {
        unsigned const n_integration_points =
            this->integration_method_.getNumberOfPoints();
 
        _ip_data.resize(n_integration_points);
        _secondary_data.N.resize(n_integration_points);
 
        auto const shape_matrices =
            NumLib::initShapeMatrices<ShapeFunction, ShapeMatricesType,
                                      DisplacementDim>(
                e, is_axially_symmetric, this->integration_method_);
 
        for (unsigned ip = 0; ip < n_integration_points; ip++)
        {
            auto& ip_data = _ip_data[ip];
            auto const& sm = shape_matrices[ip];
            _ip_data[ip].integration_weight =
                this->integration_method_.getWeightedPoint(ip).getWeight() *
                sm.integralMeasure * sm.detJ;
 
            ip_data.N_u = sm.N;
            ip_data.dNdx_u = sm.dNdx;
 
            _secondary_data.N[ip] = shape_matrices[ip].N;
        }
    }

References ProcessLib::LargeDeformation::LargeDeformationLocalAssemblerInterface< DisplacementDim >::LargeDeformationLocalAssemblerInterface(), _ip_data, _secondary_data, NumLib::initShapeMatrices(), and ProcessLib::LargeDeformation::LargeDeformationLocalAssemblerInterface< DisplacementDim >::integration_method_.

Member Function Documentation

◆ assemble()

template<typename ShapeFunction, int DisplacementDim>

void ProcessLib::LargeDeformation::LargeDeformationLocalAssembler< ShapeFunction, 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 280 of file LargeDeformationFEM.h.

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

References OGS_FATAL.

◆ assembleWithJacobian()

template<typename ShapeFunction, int DisplacementDim>

void ProcessLib::LargeDeformation::LargeDeformationLocalAssembler< ShapeFunction, DisplacementDim >::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 )

inlineoverridevirtual

Reimplemented from ProcessLib::LocalAssemblerInterface.

Definition at line 319 of file LargeDeformationFEM.h.

    {
        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);
 
        auto [u] = localDOF(local_x);
        auto [u_prev] = localDOF(local_x_prev);
 
        unsigned const n_integration_points =
            this->integration_method_.getNumberOfPoints();
 
        typename ConstitutiveRelations::ConstitutiveSetting<DisplacementDim>
            constitutive_setting;
        auto const& medium =
            *this->process_data_.media_map.getMedium(this->element_.getID());
 
        auto const f_bar_variables =
            computeFBarVariables(/*compute_detF0_only?*/ false, u);
 
        for (unsigned ip = 0; ip < n_integration_points; ip++)
        {
            auto const& w = _ip_data[ip].integration_weight;
            auto const& N = _ip_data[ip].N_u;
            auto const& dNdx = _ip_data[ip].dNdx_u;
 
            ParameterLib::SpatialPosition const x_position{
                std::nullopt, this->element_.getID(),
                MathLib::Point3d(
                    NumLib::interpolateCoordinates<ShapeFunction,
                                                   ShapeMatricesType>(
                        this->element_, N))};
 
            auto const x_coord = x_position.getCoordinates().value()[0];
 
            // For the 2D case the 33-component is needed (and the four entries
            // of the non-symmetric matrix); In 3d there are nine entries.
            GradientMatrixType G(DisplacementDim * DisplacementDim +
                                     (DisplacementDim == 2 ? 1 : 0),
                                 ShapeFunction::NPOINTS * DisplacementDim);
            Deformation::computeGMatrix<DisplacementDim,
                                        ShapeFunction::NPOINTS>(
                dNdx, G, this->is_axially_symmetric_, N, x_coord);
 
            GradientVectorType const grad_u = G * u;
 
            auto const B_0 = LinearBMatrix::computeBMatrix<
                DisplacementDim, ShapeFunction::NPOINTS,
                typename BMatricesType::BMatrixType>(
                dNdx, N, x_coord, this->is_axially_symmetric_);
 
            auto const B_N = NonLinearBMatrix::computeBMatrix<
                DisplacementDim, ShapeFunction::NPOINTS,
                typename BMatricesType::BMatrixType>(
                dNdx, N, grad_u, x_coord, this->is_axially_symmetric_);
 
            double const detF0 =
                f_bar_variables ? std::get<double>(*f_bar_variables) : 0.0;
            double const alpha = computeOutputStrainData(
                detF0, B_0 + 0.5 * B_N, grad_u, u, this->output_data_[ip]);
 
            auto const CD = updateConstitutiveRelations(
                alpha, G, u_prev, x_position, t, dt, constitutive_setting,
                medium, this->current_states_[ip], this->prev_states_[ip],
                this->material_states_[ip], this->output_data_[ip]);
 
            BMatrixType B = B_0 + B_N;
 
            auto const& sigma =
                std::get<StressData<DisplacementDim>>(this->current_states_[ip])
                    .sigma;
            auto const& b = *std::get<VolumetricBodyForce<DisplacementDim>>(CD);
            auto const& C =
                std::get<ConstitutiveRelations::SolidMechanicsDataStateless<
                    DisplacementDim>>(CD)
                    .stiffness_tensor;
 
            auto const sigma_geom =
                computeSigmaGeom<DisplacementDim, ShapeMatricesType>(
                    MathLib::KelvinVector::kelvinVectorToTensor(sigma));
 
            if (!f_bar_variables)
            {
                local_b.noalias() -=
                    (B.transpose() * sigma - N_u_op(N).transpose() * b) * w;
 
                local_Jac.noalias() += G.transpose() * sigma_geom * G * w;
 
                local_Jac.noalias() += B.transpose() * C * B * w;
            }
            else  // with F bar
            {
                double const alpha_p2 = alpha * alpha;
                local_Jac.noalias() +=
                    alpha_p2 * G.transpose() * sigma_geom * G * w;
 
                auto const q0 = std::get<VectorTypeForFbar>(*f_bar_variables);
 
                auto const& F =
                    std::get<DeformationGradientData<DisplacementDim>>(
                        this->output_data_[ip])
                        .deformation_gradient;
                VectorTypeForFbar q = NonLinearFbar::computeQVector<
                    DisplacementDim, ShapeFunction::NPOINTS, VectorTypeForFbar,
                    GradientVectorType,
                    typename ShapeMatricesType::GlobalDimNodalMatrixType>(
                    dNdx, F, this->is_axially_symmetric_);
 
                auto const q0_q = q0 - q;
 
                auto const S_q = sigma * (q0_q).transpose();
 
                local_Jac.noalias() +=
                    2.0 * alpha_p2 * S_q.transpose() * B * w / DisplacementDim;
 
                auto const twoEplsI =
                    NonLinearFbar::compute2EPlusI<DisplacementDim>(
                        alpha_p2,
                        std::get<StrainData<DisplacementDim>>(
                            this->output_data_[ip])
                            .eps,
                        this->is_axially_symmetric_);
 
                // B bar:
                BMatrixType const Bbar =
                    alpha_p2 *
                    (B + twoEplsI * (q0_q).transpose() / DisplacementDim);
 
                local_Jac.noalias() +=
                    2.0 * Bbar.transpose() * S_q * w / DisplacementDim;
 
                local_Jac.noalias() += Bbar.transpose() * C * Bbar * w;
 
                local_Jac.noalias() -=
                    alpha_p2 * sigma.dot(twoEplsI) *
                    (q0 * q0.transpose() - q * q.transpose()) * w /
                    DisplacementDim;
 
                local_b.noalias() -=
                    (Bbar.transpose() * sigma - N_u_op(N).transpose() * b) * w;
            }
        }
    }

◆ computeFBarVariables()

template<typename ShapeFunction, int DisplacementDim>

std::optional< std::tuple< double, VectorTypeForFbar > > ProcessLib::LargeDeformation::LargeDeformationLocalAssembler< ShapeFunction, DisplacementDim >::computeFBarVariables	(	bool const	compute_detF0_only,
		Eigen::VectorXd const &	local_x ) const

inline

Definition at line 292 of file LargeDeformationFEM.h.

    {
        if (this->process_data_.bar_det_f_type ==
            NonLinearFbar::BarDetFType::NONE)
        {
            return std::nullopt;
        }
 
        if (this->process_data_.bar_det_f_type ==
            NonLinearFbar::BarDetFType::ELEMENT_AVERAGE)
        {
            return NonLinearFbar::computeFBarInitialVariablesAverage<
                DisplacementDim, GradientVectorType, VectorTypeForFbar,
                NodalVectorType, ShapeFunction, ShapeMatricesType, IpData>(
                _ip_data, compute_detF0_only, local_x,
                this->integration_method_, this->element_,
                this->is_axially_symmetric_);
        }
 
        return NonLinearFbar::computeFBarInitialVariablesElementCenter<
            DisplacementDim, GradientVectorType, GradientMatrixType,
            VectorTypeForFbar, ShapeFunction, ShapeMatricesType>(
            compute_detF0_only, local_x, this->element_,
            this->is_axially_symmetric_);
    }

References _ip_data, ProcessLib::NonLinearFbar::computeFBarInitialVariablesAverage(), ProcessLib::NonLinearFbar::computeFBarInitialVariablesElementCenter(), ProcessLib::LargeDeformation::LargeDeformationLocalAssemblerInterface< DisplacementDim >::element_, ProcessLib::NonLinearFbar::ELEMENT_AVERAGE, ProcessLib::LargeDeformation::LargeDeformationLocalAssemblerInterface< DisplacementDim >::integration_method_, ProcessLib::LargeDeformation::LargeDeformationLocalAssemblerInterface< DisplacementDim >::is_axially_symmetric_, ProcessLib::NonLinearFbar::NONE, and ProcessLib::LargeDeformation::LargeDeformationLocalAssemblerInterface< DisplacementDim >::process_data_.

Referenced by assembleWithJacobian(), and postTimestepConcrete().

◆ computeOutputStrainData()

template<typename ShapeFunction, int DisplacementDim>

double ProcessLib::LargeDeformation::LargeDeformationLocalAssembler< ShapeFunction, DisplacementDim >::computeOutputStrainData	(	double const	detF0,
		BMatrixType const &	B,
		GradientVectorType const &	grad_u,
		Eigen::Ref< Eigen::VectorXd const > const &	u,
		typename ConstitutiveRelations::OutputData< DisplacementDim > &	output_data ) const

inline

Definition at line 185 of file LargeDeformationFEM.h.

    {
        // Note: Here B=B_{linear}+B_{nonlinear}/2  For Green-Lagrange strain.
        auto& eps_data = std::get<StrainData<DisplacementDim>>(output_data);
        auto& deformation_gradient_data =
            std::get<DeformationGradientData<DisplacementDim>>(output_data);
 
        eps_data.eps = B * u;
        deformation_gradient_data.deformation_gradient =
            grad_u + MathLib::VectorizedTensor::identity<DisplacementDim>();
        deformation_gradient_data.volume_ratio =
            MathLib::VectorizedTensor::determinant(
                deformation_gradient_data.deformation_gradient);
 
        if (this->process_data_.bar_det_f_type ==
            NonLinearFbar::BarDetFType::NONE)
        {
            return 1.0;
        }
 
        double const detF = deformation_gradient_data.volume_ratio;
        double const J_ratio = detF0 / detF;
 
        if (J_ratio < .0)
        {
            OGS_FATAL(
                "det(F0)/det(F) is negative with det(F0) = {:g}, and det(F) = "
                "{:g} ",
                detF0, detF);
        }
 
        double const alpha =
            DisplacementDim == 3 ? std::cbrt(J_ratio) : std::sqrt(J_ratio);
 
        deformation_gradient_data.deformation_gradient
            .template segment<DisplacementDim * DisplacementDim>(0) *= alpha;
        deformation_gradient_data.volume_ratio *=
            std::pow(alpha, DisplacementDim);
 
        double const alpha_p2 = alpha * alpha;
        eps_data.eps = alpha_p2 * eps_data.eps +
                       0.5 * (alpha_p2 - 1) *
                           NonLinearFbar::identityForF<DisplacementDim>(
                               this->is_axially_symmetric_);
        return alpha;
    }

References MathLib::VectorizedTensor::determinant(), MathLib::VectorizedTensor::identity(), ProcessLib::NonLinearFbar::identityForF(), ProcessLib::LargeDeformation::LargeDeformationLocalAssemblerInterface< DisplacementDim >::is_axially_symmetric_, ProcessLib::NonLinearFbar::NONE, OGS_FATAL, and ProcessLib::LargeDeformation::LargeDeformationLocalAssemblerInterface< DisplacementDim >::process_data_.

Referenced by assembleWithJacobian(), and postTimestepConcrete().

◆ getShapeMatrix()

template<typename ShapeFunction, int DisplacementDim>

Eigen::Map< const Eigen::RowVectorXd > ProcessLib::LargeDeformation::LargeDeformationLocalAssembler< ShapeFunction, DisplacementDim >::getShapeMatrix ( const unsigned integration_point ) const

inlineoverridevirtual

Provides the shape matrix at the given integration point.

Implements NumLib::ExtrapolatableElement.

Definition at line 544 of file LargeDeformationFEM.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 _secondary_data.

◆ initializeConcrete()

template<typename ShapeFunction, int DisplacementDim>

void ProcessLib::LargeDeformation::LargeDeformationLocalAssembler< ShapeFunction, DisplacementDim >::initializeConcrete ( )

inlineoverridevirtual

Set initial stress from parameter.

Reimplemented from ProcessLib::LocalAssemblerInterface.

Definition at line 144 of file LargeDeformationFEM.h.

    {
        unsigned const n_integration_points =
            this->integration_method_.getNumberOfPoints();
        for (unsigned ip = 0; ip < n_integration_points; ip++)
        {
            auto const& ip_data = _ip_data[ip];
 
            ParameterLib::SpatialPosition const x_position{
                std::nullopt, this->element_.getID(),
                MathLib::Point3d(
                    NumLib::interpolateCoordinates<ShapeFunction,
                                                   ShapeMatricesType>(
                        this->element_, ip_data.N_u))};
            if (this->process_data_.initial_stress != nullptr)
            {
                std::get<StressData<DisplacementDim>>(this->current_states_[ip])
                    .sigma.noalias() =
                    MathLib::KelvinVector::symmetricTensorToKelvinVector<
                        DisplacementDim>((*this->process_data_.initial_stress)(
                        std::numeric_limits<
                            double>::quiet_NaN() /* time independent */,
                        x_position));
            }
 
            double const t = 0;  // TODO (naumov) pass t from top
            auto& material_state = this->material_states_[ip];
            this->solid_material_.initializeInternalStateVariables(
                t, x_position, *material_state.material_state_variables);
 
            typename ConstitutiveRelations::ConstitutiveSetting<DisplacementDim>
                constitutive_setting;
            constitutive_setting.init();
 
            material_state.pushBackState();
            this->prev_states_[ip] = this->current_states_[ip];
        }
    }

◆ localDOF()

template<typename ShapeFunction, int DisplacementDim>

constexpr auto ProcessLib::LargeDeformation::LargeDeformationLocalAssembler< ShapeFunction, DisplacementDim >::localDOF ( std::vector< double > const & x )

inlinestaticconstexprprivate

Definition at line 554 of file LargeDeformationFEM.h.

    {
        return NumLib::localDOF<
            NumLib::Vectorial<ShapeFunction, DisplacementDim>>(x);
    }

References NumLib::localDOF().

Referenced by assembleWithJacobian().

◆ postTimestepConcrete()

template<typename ShapeFunction, int DisplacementDim>

void ProcessLib::LargeDeformation::LargeDeformationLocalAssembler< ShapeFunction, DisplacementDim >::postTimestepConcrete	(	Eigen::VectorXd const &	local_x,
		Eigen::VectorXd const &	local_x_prev,
		double const	t,
		double const	dt,
		int const	)

inlineoverridevirtual

Reimplemented from ProcessLib::LocalAssemblerInterface.

Definition at line 471 of file LargeDeformationFEM.h.

    {
        unsigned const n_integration_points =
            this->integration_method_.getNumberOfPoints();
 
        typename ConstitutiveRelations::ConstitutiveSetting<DisplacementDim>
            constitutive_setting;
 
        auto& medium =
            *this->process_data_.media_map.getMedium(this->element_.getID());
 
        auto const f_bar_variables =
            computeFBarVariables(/*compute_detF0_only?*/ true, local_x);
 
        for (unsigned ip = 0; ip < n_integration_points; ip++)
        {
            auto const& N = _ip_data[ip].N_u;
            auto const& dNdx = _ip_data[ip].dNdx_u;
 
            ParameterLib::SpatialPosition const x_position{
                std::nullopt, this->element_.getID(),
                MathLib::Point3d(
                    NumLib::interpolateCoordinates<ShapeFunction,
                                                   ShapeMatricesType>(
                        this->element_, N))};
 
            auto const x_coord = x_position.getCoordinates().value()[0];
 
            // For the 2D case the 33-component is needed (and the four entries
            // of the non-symmetric matrix); In 3d there are nine entries.
            GradientMatrixType G(DisplacementDim * DisplacementDim +
                                     (DisplacementDim == 2 ? 1 : 0),
                                 ShapeFunction::NPOINTS * DisplacementDim);
            Deformation::computeGMatrix<DisplacementDim,
                                        ShapeFunction::NPOINTS>(
                dNdx, G, this->is_axially_symmetric_, N, x_coord);
 
            GradientVectorType const grad_u = G * local_x;
 
            // B for Green-Lagrange strain.
            auto B = LinearBMatrix::computeBMatrix<
                DisplacementDim, ShapeFunction::NPOINTS,
                typename BMatricesType::BMatrixType>(
                dNdx, N, x_coord, this->is_axially_symmetric_);
 
            B.noalias() += 0.5 * NonLinearBMatrix::computeBMatrix<
                                     DisplacementDim, ShapeFunction::NPOINTS,
                                     typename BMatricesType::BMatrixType>(
                                     dNdx, N, grad_u, x_coord,
                                     this->is_axially_symmetric_);
 
            double const detF0 =
                f_bar_variables ? std::get<double>(*f_bar_variables) : 0.0;
            double const alpha = computeOutputStrainData(
                detF0, B, grad_u, local_x, this->output_data_[ip]);
 
            updateConstitutiveRelations(
                alpha, G, local_x_prev, x_position, t, dt, constitutive_setting,
                medium, this->current_states_[ip], this->prev_states_[ip],
                this->material_states_[ip], this->output_data_[ip]);
 
            this->material_states_[ip].pushBackState();
        }
 
        for (unsigned ip = 0; ip < n_integration_points; ip++)
        {
            this->prev_states_[ip] = this->current_states_[ip];
        }
    }

◆ updateConstitutiveRelations()

template<typename ShapeFunction, int DisplacementDim>

ConstitutiveRelations::ConstitutiveData< DisplacementDim > ProcessLib::LargeDeformation::LargeDeformationLocalAssembler< ShapeFunction, DisplacementDim >::updateConstitutiveRelations	(	double const	alpha,
		GradientMatrixType const &	G,
		Eigen::Ref< Eigen::VectorXd const > const &	u_prev,
		ParameterLib::SpatialPosition const &	x_position,
		double const	t,
		double const	dt,
		typename ConstitutiveRelations::ConstitutiveSetting< DisplacementDim > &	CS,
		MaterialPropertyLib::Medium const &	medium,
		typename ConstitutiveRelations::StatefulData< DisplacementDim > &	current_state,
		typename ConstitutiveRelations::StatefulDataPrev< DisplacementDim > const &	prev_state,
		MaterialStateData< DisplacementDim > &	material_state,
		typename ConstitutiveRelations::OutputData< DisplacementDim > &	output_data ) const

inline

Definition at line 240 of file LargeDeformationFEM.h.

    {
        double const T_ref =
            this->process_data_.reference_temperature
                ? (*this->process_data_.reference_temperature)(t, x_position)[0]
                : std::numeric_limits<double>::quiet_NaN();
 
        auto models =
            ConstitutiveRelations::createConstitutiveModels<DisplacementDim>(
                this->process_data_, this->solid_material_);
        typename ConstitutiveRelations::ConstitutiveTempData<DisplacementDim>
            tmp;
        typename ConstitutiveRelations::ConstitutiveData<DisplacementDim> CD;
 
        CS.eval(
            models, t, dt, x_position,  //
            medium,                     //
            T_ref,                      //
            std::get<DeformationGradientData<DisplacementDim>>(output_data),
            alpha * (G * u_prev +
                     MathLib::VectorizedTensor::identity<DisplacementDim>()),
            current_state, prev_state, material_state, tmp, CD);
 
        return CD;
    }

References ProcessLib::LargeDeformation::ConstitutiveRelations::createConstitutiveModels(), ProcessLib::LargeDeformation::ConstitutiveRelations::ConstitutiveSetting< DisplacementDim >::eval(), MathLib::VectorizedTensor::identity(), ProcessLib::LargeDeformation::LargeDeformationLocalAssemblerInterface< DisplacementDim >::process_data_, and ProcessLib::LargeDeformation::LargeDeformationLocalAssemblerInterface< DisplacementDim >::solid_material_.

Referenced by assembleWithJacobian(), and postTimestepConcrete().

Member Data Documentation

◆ _ip_data

template<typename ShapeFunction, int DisplacementDim>

std::vector<IpData, Eigen::aligned_allocator<IpData> > ProcessLib::LargeDeformation::LargeDeformationLocalAssembler< ShapeFunction, DisplacementDim >::_ip_data

private

Definition at line 561 of file LargeDeformationFEM.h.

Referenced by LargeDeformationLocalAssembler(), assembleWithJacobian(), computeFBarVariables(), initializeConcrete(), and postTimestepConcrete().

◆ _secondary_data

template<typename ShapeFunction, int DisplacementDim>

SecondaryData<typename ShapeMatrices::ShapeType> ProcessLib::LargeDeformation::LargeDeformationLocalAssembler< ShapeFunction, DisplacementDim >::_secondary_data

private

Definition at line 563 of file LargeDeformationFEM.h.

Referenced by LargeDeformationLocalAssembler(), and getShapeMatrix().

◆ displacement_size

template<typename ShapeFunction, int DisplacementDim>

const int ProcessLib::LargeDeformation::LargeDeformationLocalAssembler< ShapeFunction, DisplacementDim >::displacement_size

staticprivate

Initial value:

=

ShapeFunction::NPOINTS * DisplacementDim

Definition at line 565 of file LargeDeformationFEM.h.

◆ N_u_op

template<typename ShapeFunction, int DisplacementDim>

auto& ProcessLib::LargeDeformation::LargeDeformationLocalAssembler< ShapeFunction, DisplacementDim >::N_u_op

staticconstexpr

Initial value:

= MathLib::eigenBlockMatrixView<

DisplacementDim, typename ShapeMatricesType::NodalRowVectorType>

MathLib::eigenBlockMatrixView

constexpr Eigen::CwiseNullaryOp< EigenBlockMatrixViewFunctor< D, M >, typename EigenBlockMatrixViewFunctor< D, M >::Matrix > eigenBlockMatrixView(const Eigen::MatrixBase< M > &matrix)

Definition EigenBlockMatrixView.h:43

EigenFixedShapeMatrixPolicy::NodalRowVectorType

RowVectorType< ShapeFunction::NPOINTS > NodalRowVectorType

Definition ShapeMatrixPolicy.h:71

Definition at line 102 of file LargeDeformationFEM.h.

Referenced by assembleWithJacobian().

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

LargeDeformationFEM.h

Detailed Description

Public Types

Public Member Functions

Static Public Attributes

Static Private Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Member Typedef Documentation

◆ BMatricesType

◆ BMatrixType

◆ GlobalDimNodalMatrixType

◆ GMatricesType

◆ GradientMatrixType

◆ GradientVectorType

◆ IpData

◆ NodalDisplacementVectorType

◆ NodalForceVectorType

◆ NodalMatrixType

◆ NodalVectorType

◆ ShapeMatrices

◆ ShapeMatricesType

◆ StiffnessMatrixType

◆ VectorTypeForFbar

Constructor & Destructor Documentation

◆ LargeDeformationLocalAssembler() [1/3]

◆ LargeDeformationLocalAssembler() [2/3]

◆ LargeDeformationLocalAssembler() [3/3]

Member Function Documentation

◆ assemble()

◆ assembleWithJacobian()

◆ computeFBarVariables()

◆ computeOutputStrainData()

◆ getShapeMatrix()

◆ initializeConcrete()

◆ localDOF()

◆ postTimestepConcrete()

◆ updateConstitutiveRelations()

Member Data Documentation

◆ _ip_data

◆ _secondary_data

◆ displacement_size

◆ N_u_op