Detailed Description

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

Definition at line 86 of file SmallDeformationFEM.h.

#include <SmallDeformationFEM.h>

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

[legend]

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

[legend]

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	BMatrixType = typename BMatricesType::BMatrixType

using	StiffnessMatrixType = typename BMatricesType::StiffnessMatrixType

using	NodalForceVectorType = typename BMatricesType::NodalForceVectorType

using	NodalDisplacementVectorType = typename BMatricesType::NodalForceVectorType

using	GMatricesType = GMatrixPolicyType< ShapeFunction, DisplacementDim >

using	GradientVectorType = typename GMatricesType::GradientVectorType

using	GradientMatrixType = typename GMatricesType::GradientMatrixType

using	IpData = IntegrationPointData< BMatricesType, ShapeMatricesType, DisplacementDim >

Public Member Functions
	SmallDeformationLocalAssembler (SmallDeformationLocalAssembler const &)=delete

	SmallDeformationLocalAssembler (SmallDeformationLocalAssembler &&)=delete

	SmallDeformationLocalAssembler (MeshLib::Element const &e, std::size_t const, bool const is_axially_symmetric, unsigned const integration_order, SmallDeformationProcessData< DisplacementDim > &process_data)

std::size_t	setIPDataInitialConditions (std::string const &name, double const *values, int const integration_order) override
	Returns number of read integration points. More...

void	initializeConcrete () override

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	postTimestepConcrete (Eigen::VectorXd const &, double const t, double const dt) override

std::vector< double > const &	getMaterialForces (std::vector< double > const &local_x, std::vector< double > &nodal_values) 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 &	getIntPtFreeEnergyDensity (const double, std::vector< GlobalVector * > const &, std::vector< NumLib::LocalToGlobalIndexMap const * > const &, std::vector< double > &cache) const override

std::size_t	setSigma (double const *values)

std::vector< double >	getSigma () const override

std::vector< double > const &	getIntPtSigma (const double, std::vector< GlobalVector * > const &, std::vector< NumLib::LocalToGlobalIndexMap const * > const &, std::vector< double > &cache) const override

std::vector< double > const &	getIntPtEpsilon (const double, std::vector< GlobalVector * > const &, std::vector< NumLib::LocalToGlobalIndexMap const * > const &, std::vector< double > &cache) const override

unsigned	getNumberOfIntegrationPoints () const override

MaterialLib::Solids::MechanicsBase< DisplacementDim >::MaterialStateVariables const &	getMaterialStateVariablesAt (unsigned integration_point) const override

void	computeSecondaryVariableConcrete (double const, double const, Eigen::VectorXd const &, 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 ProcessLib::SmallDeformation::MaterialForcesInterface
virtual	~MaterialForcesInterface ()=default

Public Member Functions inherited from NumLib::ExtrapolatableElement
virtual	~ExtrapolatableElement ()=default

Private Attributes
SmallDeformationProcessData< DisplacementDim > &	_process_data

std::vector< IpData, Eigen::aligned_allocator< IpData > >	_ip_data

IntegrationMethod	_integration_method

MeshLib::Element const &	_element

SecondaryData< typename ShapeMatrices::ShapeType >	_secondary_data

bool const	_is_axially_symmetric

Static Private Attributes
static const int	displacement_size

Member Typedef Documentation

◆ BMatricesType

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

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

Definition at line 95 of file SmallDeformationFEM.h.

◆ BMatrixType

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

using ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::BMatrixType = typename BMatricesType::BMatrixType

Definition at line 97 of file SmallDeformationFEM.h.

◆ GMatricesType

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

using ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::GMatricesType = GMatrixPolicyType<ShapeFunction, DisplacementDim>

Definition at line 103 of file SmallDeformationFEM.h.

◆ GradientMatrixType

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

using ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::GradientMatrixType = typename GMatricesType::GradientMatrixType

Definition at line 105 of file SmallDeformationFEM.h.

◆ GradientVectorType

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

using ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::GradientVectorType = typename GMatricesType::GradientVectorType

Definition at line 104 of file SmallDeformationFEM.h.

◆ IpData

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

using ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::IpData = IntegrationPointData<BMatricesType, ShapeMatricesType, DisplacementDim>

Definition at line 106 of file SmallDeformationFEM.h.

◆ NodalDisplacementVectorType

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

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

Definition at line 100 of file SmallDeformationFEM.h.

◆ NodalForceVectorType

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

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

Definition at line 99 of file SmallDeformationFEM.h.

◆ NodalMatrixType

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

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

Definition at line 92 of file SmallDeformationFEM.h.

◆ NodalVectorType

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

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

Definition at line 93 of file SmallDeformationFEM.h.

◆ ShapeMatrices

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

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

Definition at line 94 of file SmallDeformationFEM.h.

◆ ShapeMatricesType

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

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

Definition at line 90 of file SmallDeformationFEM.h.

◆ StiffnessMatrixType

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

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

Definition at line 98 of file SmallDeformationFEM.h.

Constructor & Destructor Documentation

◆ SmallDeformationLocalAssembler() [1/3]

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

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

delete

◆ SmallDeformationLocalAssembler() [2/3]

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

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

delete

◆ SmallDeformationLocalAssembler() [3/3]

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

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

inline

Definition at line 113 of file SmallDeformationFEM.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[ip].integration_weight =
                 _integration_method.getWeightedPoint(ip).getWeight() *
                 sm.integralMeasure * sm.detJ;
  
             ip_data.N = sm.N;
             ip_data.dNdx = sm.dNdx;
  
             static const int kelvin_vector_size =
                 MathLib::KelvinVector::kelvin_vector_dimensions(
                     DisplacementDim);
             // Initialize current time step values
             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);
  
             _secondary_data.N[ip] = shape_matrices[ip].N;
         }
     }

References ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::_integration_method, ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::_ip_data, ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::_process_data, ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< 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::SmallDeformation::SmallDeformationLocalAssembler< 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 229 of file SmallDeformationFEM.h.

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

References OGS_FATAL.

◆ assembleWithJacobian()

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

void ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< 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
	)

inlineoverridevirtual

Reimplemented from ProcessLib::LocalAssemblerInterface.

Definition at line 241 of file SmallDeformationFEM.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);
  
         unsigned const n_integration_points =
             _integration_method.getNumberOfPoints();
  
         MPL::VariableArray variables_prev;
         MPL::VariableArray variables;
         ParameterLib::SpatialPosition x_position;
         x_position.setElementID(_element.getID());
  
         auto const& b = _process_data.specific_body_force;
  
         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;
  
             typename ShapeMatricesType::template MatrixType<DisplacementDim,
                                                             displacement_size>
                 N_u_op = ShapeMatricesType::template MatrixType<
                     DisplacementDim,
                     displacement_size>::Zero(DisplacementDim,
                                              displacement_size);
             for (int i = 0; i < DisplacementDim; ++i)
             {
                 N_u_op
                     .template block<1, displacement_size / DisplacementDim>(
                         i, i * displacement_size / DisplacementDim)
                     .noalias() = N;
             }
  
             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_prev[static_cast<int>(MPL::Variable::stress)]
                 .emplace<
                     MathLib::KelvinVector::KelvinVectorType<DisplacementDim>>(
                     sigma_prev);
             variables_prev[static_cast<int>(MPL::Variable::mechanical_strain)]
                 .emplace<
                     MathLib::KelvinVector::KelvinVectorType<DisplacementDim>>(
                     eps_prev);
  
             double const T_ref =
                 _process_data.reference_temperature
                     ? (*_process_data.reference_temperature)(t, x_position)[0]
                     : std::numeric_limits<double>::quiet_NaN();
  
             variables_prev[static_cast<int>(MPL::Variable::temperature)]
                 .emplace<double>(T_ref);
             variables[static_cast<int>(MPL::Variable::mechanical_strain)]
                 .emplace<
                     MathLib::KelvinVector::KelvinVectorType<DisplacementDim>>(
                     eps);
             variables[static_cast<int>(MPL::Variable::temperature)]
                 .emplace<double>(T_ref);
  
             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);
  
             auto const rho = _process_data.solid_density(t, x_position)[0];
             local_b.noalias() -=
                 (B.transpose() * sigma - N_u_op.transpose() * rho * b) * w;
             local_Jac.noalias() += B.transpose() * C * B * w;
         }
     }

References ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::_element, ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::_integration_method, ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::_ip_data, ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::_is_axially_symmetric, ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::_process_data, ProcessLib::LinearBMatrix::computeBMatrix(), ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::displacement_size, MeshLib::Element::getID(), NumLib::interpolateXCoordinate(), OGS_FATAL, ParameterLib::SpatialPosition::setElementID(), ParameterLib::SpatialPosition::setIntegrationPoint(), and MaterialPropertyLib::temperature.

◆ computeSecondaryVariableConcrete()

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

void ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::computeSecondaryVariableConcrete	(	double const	,
		double const	,
		Eigen::VectorXd const &	,
		Eigen::VectorXd const &
	)

inlineoverridevirtual

Reimplemented from ProcessLib::LocalAssemblerInterface.

Definition at line 462 of file SmallDeformationFEM.h.

     {
         int const elem_id = _element.getID();
         ParameterLib::SpatialPosition x_position;
         x_position.setElementID(elem_id);
         unsigned const n_integration_points =
             _integration_method.getNumberOfPoints();
  
         auto sigma_sum = MathLib::KelvinVector::tensorToKelvin<DisplacementDim>(
             Eigen::Matrix<double, 3, 3>::Zero());
  
         for (unsigned ip = 0; ip < n_integration_points; ip++)
         {
             x_position.setIntegrationPoint(ip);
             auto const& sigma = _ip_data[ip].sigma;
             sigma_sum += sigma;
         }
  
         Eigen::Matrix<double, 3, 3, 0, 3, 3> const sigma_avg =
             MathLib::KelvinVector::kelvinVectorToTensor(sigma_sum) /
             n_integration_points;
  
         Eigen::SelfAdjointEigenSolver<Eigen::Matrix<double, 3, 3>> e_s(
             sigma_avg);
  
         Eigen::Map<Eigen::Vector3d>(
             &(*_process_data.principal_stress_values)[elem_id * 3], 3) =
             e_s.eigenvalues();
  
         auto eigen_vectors = e_s.eigenvectors();
  
         for (auto i = 0; i < 3; i++)
         {
             Eigen::Map<Eigen::Vector3d>(
                 &(*_process_data.principal_stress_vector[i])[elem_id * 3], 3) =
                 eigen_vectors.col(i);
         }
     }

References ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::_element, ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::_integration_method, ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::_ip_data, ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::_process_data, MeshLib::Element::getID(), MathLib::KelvinVector::kelvinVectorToTensor(), ParameterLib::SpatialPosition::setElementID(), and ParameterLib::SpatialPosition::setIntegrationPoint().

◆ getIntPtEpsilon()

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

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

inlineoverridevirtual

Implements ProcessLib::SmallDeformation::SmallDeformationLocalAssemblerInterface< DisplacementDim >.

Definition at line 440 of file SmallDeformationFEM.h.

     {
         return ProcessLib::getIntegrationPointKelvinVectorData<DisplacementDim>(
             _ip_data, &IpData::eps, cache);
     }

References ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::_ip_data, and ProcessLib::SmallDeformation::IntegrationPointData< BMatricesType, ShapeMatricesType, DisplacementDim >::eps.

◆ getIntPtFreeEnergyDensity()

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

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

inlineoverridevirtual

Implements ProcessLib::SmallDeformation::SmallDeformationLocalAssemblerInterface< DisplacementDim >.

Definition at line 399 of file SmallDeformationFEM.h.

     {
         cache.clear();
         cache.reserve(_ip_data.size());
  
         transform(
             cbegin(_ip_data), cend(_ip_data), back_inserter(cache),
             [](auto const& ip_data) { return ip_data.free_energy_density; });
  
         return cache;
     }

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

◆ getIntPtSigma()

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

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

inlineoverridevirtual

Implements ProcessLib::SmallDeformation::SmallDeformationLocalAssemblerInterface< DisplacementDim >.

Definition at line 430 of file SmallDeformationFEM.h.

     {
         return ProcessLib::getIntegrationPointKelvinVectorData<DisplacementDim>(
             _ip_data, &IpData::sigma, cache);
     }

References ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::_ip_data, and ProcessLib::SmallDeformation::IntegrationPointData< BMatricesType, ShapeMatricesType, DisplacementDim >::sigma.

◆ getMaterialForces()

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

std::vector<double> const& ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::getMaterialForces	(	std::vector< double > const &	local_x,
		std::vector< double > &	nodal_values
	)

inlineoverridevirtual

Implements ProcessLib::SmallDeformation::MaterialForcesInterface.

Definition at line 378 of file SmallDeformationFEM.h.

     {
         return ProcessLib::SmallDeformation::getMaterialForces<
             DisplacementDim, ShapeFunction, ShapeMatricesType,
             typename BMatricesType::NodalForceVectorType,
             NodalDisplacementVectorType, GradientVectorType,
             GradientMatrixType>(local_x, nodal_values, _integration_method,
                                 _ip_data, _element, _is_axially_symmetric);
     }

References ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::_element, ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::_integration_method, ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::_ip_data, ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::_is_axially_symmetric, and ProcessLib::SmallDeformation::getMaterialForces().

◆ getMaterialStateVariablesAt()

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

MaterialLib::Solids::MechanicsBase< DisplacementDim>::MaterialStateVariables const& ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::getMaterialStateVariablesAt ( unsigned integration_point ) const

inlineoverridevirtual

Implements ProcessLib::SmallDeformation::SmallDeformationLocalAssemblerInterface< DisplacementDim >.

Definition at line 457 of file SmallDeformationFEM.h.

     {
         return *_ip_data[integration_point].material_state_variables;
     }

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

◆ getNumberOfIntegrationPoints()

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

unsigned ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::getNumberOfIntegrationPoints ( ) const

inlineoverridevirtual

Implements ProcessLib::SmallDeformation::SmallDeformationLocalAssemblerInterface< DisplacementDim >.

Definition at line 450 of file SmallDeformationFEM.h.

     {
         return _integration_method.getNumberOfPoints();
     }

References ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::_integration_method.

◆ getShapeMatrix()

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

Eigen::Map<const Eigen::RowVectorXd> ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< 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 390 of file SmallDeformationFEM.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::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::_secondary_data, and ProcessLib::HeatTransportBHE::SecondaryData< ShapeMatrixType >::N.

◆ getSigma()

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

std::vector<double> ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::getSigma ( ) const

inlineoverridevirtual

Implements ProcessLib::SmallDeformation::SmallDeformationLocalAssemblerInterface< DisplacementDim >.

Definition at line 424 of file SmallDeformationFEM.h.

     {
         return ProcessLib::getIntegrationPointKelvinVectorData<DisplacementDim>(
             _ip_data, &IpData::sigma);
     }

References ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::_ip_data, and ProcessLib::SmallDeformation::IntegrationPointData< BMatricesType, ShapeMatricesType, DisplacementDim >::sigma.

◆ initializeConcrete()

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

void ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::initializeConcrete ( )

inlineoverridevirtual

Set initial stress from parameter.

Reimplemented from ProcessLib::LocalAssemblerInterface.

Definition at line 199 of file SmallDeformationFEM.h.

     {
         unsigned const n_integration_points =
             _integration_method.getNumberOfPoints();
         for (unsigned ip = 0; ip < n_integration_points; ip++)
         {
             auto& ip_data = _ip_data[ip];
  
             if (_process_data.initial_stress != nullptr)
             {
                 ParameterLib::SpatialPosition const x_position{
                     std::nullopt, _element.getID(), ip,
                     MathLib::Point3d(
                         NumLib::interpolateCoordinates<ShapeFunction,
                                                        ShapeMatricesType>(
                             _element, ip_data.N))};
  
                 ip_data.sigma =
                     MathLib::KelvinVector::symmetricTensorToKelvinVector<
                         DisplacementDim>((*_process_data.initial_stress)(
                         std::numeric_limits<
                             double>::quiet_NaN() /* time independent */,
                         x_position));
             }
  
             ip_data.pushBackState();
         }
     }

References ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::_element, ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::_integration_method, ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::_ip_data, ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::_process_data, MeshLib::Element::getID(), NumLib::interpolateCoordinates(), and MathLib::KelvinVector::symmetricTensorToKelvinVector().

◆ postTimestepConcrete()

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

void ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::postTimestepConcrete	(	Eigen::VectorXd const &	,
		double const	t,
		double const	dt
	)

inlineoverridevirtual

Reimplemented from ProcessLib::LocalAssemblerInterface.

Definition at line 350 of file SmallDeformationFEM.h.

     {
         unsigned const n_integration_points =
             _integration_method.getNumberOfPoints();
  
         for (unsigned ip = 0; ip < n_integration_points; ip++)
         {
             _ip_data[ip].pushBackState();
         }
  
         ParameterLib::SpatialPosition x_position;
         x_position.setElementID(_element.getID());
  
         for (unsigned ip = 0; ip < n_integration_points; ip++)
         {
             x_position.setIntegrationPoint(ip);
  
             auto& ip_data = _ip_data[ip];
  
             // Update free energy density needed for material forces.
             ip_data.free_energy_density =
                 ip_data.solid_material.computeFreeEnergyDensity(
                     t, x_position, dt, ip_data.eps, ip_data.sigma,
                     *ip_data.material_state_variables);
         }
     }

References ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::_element, ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::_integration_method, ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::_ip_data, MeshLib::Element::getID(), ParameterLib::SpatialPosition::setElementID(), and ParameterLib::SpatialPosition::setIntegrationPoint().

◆ setIPDataInitialConditions()

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

std::size_t ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::setIPDataInitialConditions	(	std::string const &	name,
		double const *	values,
		int const	integration_order
	)

inlineoverridevirtual

Returns number of read integration points.

Implements ProcessLib::SmallDeformation::SmallDeformationLocalAssemblerInterface< DisplacementDim >.

Definition at line 169 of file SmallDeformationFEM.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_ip")
         {
             if (_process_data.initial_stress != nullptr)
             {
                 OGS_FATAL(
                     "Setting initial conditions for stress from integration "
                     "point data and from a parameter '{:s}' is not possible "
                     "simultaneously.",
                     _process_data.initial_stress->name);
             }
             return setSigma(values);
         }
  
         return 0;
     }

References ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::_element, ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::_integration_method, ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::_process_data, MeshLib::Element::getID(), MaterialPropertyLib::name, OGS_FATAL, and ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::setSigma().

◆ setSigma()

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

std::size_t ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::setSigma ( double const * values )

inline

Definition at line 415 of file SmallDeformationFEM.h.

     {
         return ProcessLib::setIntegrationPointKelvinVectorData<DisplacementDim>(
             values, _ip_data, &IpData::sigma);
     }

References ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::_ip_data, and ProcessLib::SmallDeformation::IntegrationPointData< BMatricesType, ShapeMatricesType, DisplacementDim >::sigma.

Referenced by ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::setIPDataInitialConditions().

Member Data Documentation

◆ _element

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

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

private

Definition at line 509 of file SmallDeformationFEM.h.

Referenced by ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::assembleWithJacobian(), ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::computeSecondaryVariableConcrete(), ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::getMaterialForces(), ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::initializeConcrete(), ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::postTimestepConcrete(), and ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::setIPDataInitialConditions().

◆ _integration_method

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

IntegrationMethod ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::_integration_method

private

Definition at line 508 of file SmallDeformationFEM.h.

Referenced by ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::SmallDeformationLocalAssembler(), ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::assembleWithJacobian(), ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::computeSecondaryVariableConcrete(), ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::getMaterialForces(), ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::getNumberOfIntegrationPoints(), ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::initializeConcrete(), ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::postTimestepConcrete(), and ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::setIPDataInitialConditions().

◆ _ip_data

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

std::vector<IpData, Eigen::aligned_allocator<IpData> > ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::_ip_data

private

Definition at line 506 of file SmallDeformationFEM.h.

◆ _is_axially_symmetric

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

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

private

Definition at line 511 of file SmallDeformationFEM.h.

Referenced by ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::assembleWithJacobian(), and ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::getMaterialForces().

◆ _process_data

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

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

private

Definition at line 504 of file SmallDeformationFEM.h.

Referenced by ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::SmallDeformationLocalAssembler(), ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::assembleWithJacobian(), ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::computeSecondaryVariableConcrete(), ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::initializeConcrete(), and ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::setIPDataInitialConditions().

◆ _secondary_data

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

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

private

Definition at line 510 of file SmallDeformationFEM.h.

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

◆ displacement_size

template<typename ShapeFunction , typename IntegrationMethod , int DisplacementDim>

const int ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::displacement_size

staticprivate

Initial value:

=

ShapeFunction::NPOINTS * DisplacementDim

Definition at line 513 of file SmallDeformationFEM.h.

Referenced by ProcessLib::SmallDeformation::SmallDeformationLocalAssembler< ShapeFunction, IntegrationMethod, DisplacementDim >::assembleWithJacobian().

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

ProcessLib/SmallDeformation/SmallDeformationFEM.h

Detailed Description

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

Public Types

Public Member Functions

Private Attributes

Static Private Attributes

Member Typedef Documentation

◆ BMatricesType

◆ BMatrixType

◆ GMatricesType

◆ GradientMatrixType

◆ GradientVectorType

◆ IpData

◆ NodalDisplacementVectorType

◆ NodalForceVectorType

◆ NodalMatrixType

◆ NodalVectorType

◆ ShapeMatrices

◆ ShapeMatricesType

◆ StiffnessMatrixType

Constructor & Destructor Documentation

◆ SmallDeformationLocalAssembler() [1/3]

◆ SmallDeformationLocalAssembler() [2/3]

◆ SmallDeformationLocalAssembler() [3/3]

Member Function Documentation

◆ assemble()

◆ assembleWithJacobian()

◆ computeSecondaryVariableConcrete()

◆ getIntPtEpsilon()

◆ getIntPtFreeEnergyDensity()

◆ getIntPtSigma()

◆ getMaterialForces()

◆ getMaterialStateVariablesAt()

◆ getNumberOfIntegrationPoints()

◆ getShapeMatrix()

◆ getSigma()

◆ initializeConcrete()

◆ postTimestepConcrete()

◆ setIPDataInitialConditions()

◆ setSigma()

Member Data Documentation

◆ _element

◆ _integration_method

◆ _ip_data

◆ _is_axially_symmetric

◆ _process_data

◆ _secondary_data

◆ displacement_size

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