Detailed Description

template<int DisplacementDim>
class MaterialLib::Solids::LinearElasticOrthotropic< DisplacementDim >

Definition at line 22 of file LinearElasticOrthotropic.h.

#include <LinearElasticOrthotropic.h>

Inheritance diagram for MaterialLib::Solids::LinearElasticOrthotropic< DisplacementDim >:

Collaboration diagram for MaterialLib::Solids::LinearElasticOrthotropic< DisplacementDim >:

Classes
struct	EvaluatedMaterialProperties

struct	MaterialProperties
	Variables specific to the material model. More...

Public Types
using	KelvinVector = MathLib::KelvinVector::KelvinVectorType< DisplacementDim >

using	KelvinMatrix = MathLib::KelvinVector::KelvinMatrixType< DisplacementDim >

Public Types inherited from MaterialLib::Solids::MechanicsBase< DisplacementDim >
using	KelvinVector = MathLib::KelvinVector::KelvinVectorType< DisplacementDim >

using	KelvinMatrix = MathLib::KelvinVector::KelvinMatrixType< DisplacementDim >

Public Member Functions
	LinearElasticOrthotropic (MaterialProperties material_properties, std::optional< ParameterLib::CoordinateSystem > const &local_coordinate_system)

double	computeFreeEnergyDensity (double const, ParameterLib::SpatialPosition const &, double const, KelvinVector const &eps, KelvinVector const &sigma, typename MechanicsBase< DisplacementDim >::MaterialStateVariables const &) const override

std::optional< std::tuple< typename MechanicsBase< DisplacementDim >::KelvinVector, std::unique_ptr< typename MechanicsBase< DisplacementDim >::MaterialStateVariables >, typename MechanicsBase< DisplacementDim >::KelvinMatrix > >	integrateStress (MaterialPropertyLib::VariableArray const &variable_array_prev, MaterialPropertyLib::VariableArray const &variable_array, double const t, ParameterLib::SpatialPosition const &x, double const, typename MechanicsBase< DisplacementDim >::MaterialStateVariables const &material_state_variables) const override

KelvinMatrix	getElasticTensor (double const t, ParameterLib::SpatialPosition const &x, double const T) const

MaterialProperties	getMaterialProperties () const

double	getBulkModulus (double const t, ParameterLib::SpatialPosition const &x, KelvinMatrix const *const) const override

Public Member Functions inherited from MaterialLib::Solids::MechanicsBase< DisplacementDim >
virtual std::unique_ptr< MaterialStateVariables >	createMaterialStateVariables () const

virtual std::optional< std::tuple< KelvinVector, std::unique_ptr< MaterialStateVariables >, KelvinMatrix > >	integrateStress (MaterialPropertyLib::VariableArray const &variable_array_prev, MaterialPropertyLib::VariableArray const &variable_array, double const t, ParameterLib::SpatialPosition const &x, double const dt, MaterialStateVariables const &material_state_variables) const =0

virtual std::vector< InternalVariable >	getInternalVariables () const

virtual ConstitutiveModel	getConstitutiveModel () const
	Gets the type of constitutive model. More...

virtual double	getTemperatureRelatedCoefficient (double const, double const, ParameterLib::SpatialPosition const &, double const, double const) const

virtual double	computeFreeEnergyDensity (double const t, ParameterLib::SpatialPosition const &x, double const dt, KelvinVector const &eps, KelvinVector const &sigma, MaterialStateVariables const &material_state_variables) const =0

virtual	~MechanicsBase ()=default

Static Public Attributes
static int const	KelvinVectorSize

Protected Attributes
MaterialProperties	_mp

std::optional< ParameterLib::CoordinateSystem > const &	_local_coordinate_system

Member Typedef Documentation

◆ KelvinMatrix

template<int DisplacementDim>

using MaterialLib::Solids::LinearElasticOrthotropic< DisplacementDim >::KelvinMatrix = MathLib::KelvinVector::KelvinMatrixType<DisplacementDim>

Definition at line 130 of file LinearElasticOrthotropic.h.

◆ KelvinVector

template<int DisplacementDim>

using MaterialLib::Solids::LinearElasticOrthotropic< DisplacementDim >::KelvinVector = MathLib::KelvinVector::KelvinVectorType<DisplacementDim>

Definition at line 128 of file LinearElasticOrthotropic.h.

Constructor & Destructor Documentation

◆ LinearElasticOrthotropic()

template<int DisplacementDim>

MaterialLib::Solids::LinearElasticOrthotropic< DisplacementDim >::LinearElasticOrthotropic	(	MaterialProperties	material_properties,
		std::optional< ParameterLib::CoordinateSystem > const &	local_coordinate_system
	)

inline

Definition at line 133 of file LinearElasticOrthotropic.h.

         : _mp(std::move(material_properties)),
           _local_coordinate_system(local_coordinate_system)
     {
     }

Member Function Documentation

◆ computeFreeEnergyDensity()

template<int DisplacementDim>

double MaterialLib::Solids::LinearElasticOrthotropic< DisplacementDim >::computeFreeEnergyDensity	(	double const	,
		ParameterLib::SpatialPosition const &	,
		double const	,
		KelvinVector const &	eps,
		KelvinVector const &	sigma,
		typename MechanicsBase< DisplacementDim >::MaterialStateVariables const &
	)		const

inlineoverride

Definition at line 142 of file LinearElasticOrthotropic.h.

     {
         return eps.dot(sigma) / 2;
     }

◆ getBulkModulus()

template<int DisplacementDim>

double MaterialLib::Solids::LinearElasticOrthotropic< DisplacementDim >::getBulkModulus	(	double const	t,
		ParameterLib::SpatialPosition const &	x,
		KelvinMatrix const * const
	)		const

inlineoverridevirtual

Reimplemented from MaterialLib::Solids::MechanicsBase< DisplacementDim >.

Definition at line 174 of file LinearElasticOrthotropic.h.

     {
         auto const& mp = _mp.evaluate(t, x);
         auto const E = [&mp](int const i) { return mp.E(i); };
         auto const nu = [&mp](int const i, int const j) { return mp.nu(i, j); };
         // corresponds to 1/(I:S:I) --> Reuss bound.
         // Voigt bound would proceed as (I:C:I)/9. Use MFront model if you
         // prefer that.
         return E(1) * E(2) * E(3) /
                (E(1) * E(2) + E(1) * E(3) * (1 - 2 * nu(2, 3)) +
                 E(2) * E(3) * (1 - 2 * nu(1, 2) * nu(1, 3)));
     }

References MaterialLib::Solids::LinearElasticOrthotropic< DisplacementDim >::_mp, MaterialLib::Solids::LinearElasticOrthotropic< DisplacementDim >::EvaluatedMaterialProperties::E(), and MaterialLib::Solids::LinearElasticOrthotropic< DisplacementDim >::MaterialProperties::evaluate().

◆ getElasticTensor()

template<int DisplacementDim>

LinearElasticOrthotropic< DisplacementDim >::KelvinMatrix MaterialLib::Solids::LinearElasticOrthotropic< DisplacementDim >::getElasticTensor	(	double const	t,
		ParameterLib::SpatialPosition const &	x,
		double const	T
	)		const

Definition at line 57 of file LinearElasticOrthotropic.cpp.

 {
     using namespace MathLib::KelvinVector;
  
     auto const& mp = _mp.evaluate(t, x);
     auto const E = [&mp](int const i) { return mp.E(i); };
     auto const G = [&mp](int const i, int const j) { return mp.G(i, j); };
     auto const nu = [&mp](int const i, int const j) { return mp.nu(i, j); };
  
     KelvinMatrixType<3> S_ortho = KelvinMatrixType<3>::Zero();
     // clang-format off
     S_ortho.template topLeftCorner<3, 3>() <<
                1. / E(1), -nu(2, 1) / E(2), -nu(3, 1) / E(3),
         -nu(1, 2) / E(1),        1. / E(2), -nu(3, 2) / E(3),
         -nu(1, 3) / E(1), -nu(2, 3) / E(2),        1. / E(3);
  
     S_ortho.template bottomRightCorner<3, 3>().diagonal() <<
         1. / (2 * G(1, 2)),
         1. / (2 * G(2, 3)),
         1. / (2 * G(1, 3));
     // clang-format on
  
     KelvinMatrixType<3> const C_ortho = S_ortho.inverse();
     auto const Q = [this, &x]() -> KelvinMatrixType<3>
     {
         if (!_local_coordinate_system)
         {
             return KelvinMatrixType<3>::Identity();
         }
         Eigen::Matrix3d R = Eigen::Matrix3d::Identity();
         R.template topLeftCorner<DisplacementDim, DisplacementDim>().noalias() =
             _local_coordinate_system->transformation<DisplacementDim>(x);
         return fourthOrderRotationMatrix(R);
     }();
  
     // Rotate the forth-order tenser in Kelvin mapping with Q*C_ortho*Q^T and
     // return the top left corner block of size 4x4 for two-dimensional case or
     // the full 6x6 matrix is returned in the three-dimensional case.
     return (Q * C_ortho * Q.transpose())
         .template topLeftCorner<kelvin_vector_dimensions(DisplacementDim),
                                 kelvin_vector_dimensions(DisplacementDim)>();
 }

References MathLib::KelvinVector::fourthOrderRotationMatrix(), and MathLib::KelvinVector::kelvin_vector_dimensions().

◆ getMaterialProperties()

template<int DisplacementDim>

MaterialProperties MaterialLib::Solids::LinearElasticOrthotropic< DisplacementDim >::getMaterialProperties ( ) const

inline

Definition at line 172 of file LinearElasticOrthotropic.h.

172 { return _mp; }

References MaterialLib::Solids::LinearElasticOrthotropic< DisplacementDim >::_mp.

◆ integrateStress()

template<int DisplacementDim>

std::optional< std::tuple< typename MechanicsBase< DisplacementDim >::KelvinVector, std::unique_ptr< typename MechanicsBase< DisplacementDim >::MaterialStateVariables >, typename MechanicsBase< DisplacementDim >::KelvinMatrix > > MaterialLib::Solids::LinearElasticOrthotropic< DisplacementDim >::integrateStress	(	MaterialPropertyLib::VariableArray const &	variable_array_prev,
		MaterialPropertyLib::VariableArray const &	variable_array,
		double const	t,
		ParameterLib::SpatialPosition const &	x,
		double const	,
		typename MechanicsBase< DisplacementDim >::MaterialStateVariables const &	material_state_variables
	)		const

override

Definition at line 27 of file LinearElasticOrthotropic.cpp.

 {
     auto const& eps_m = std::get<MPL::SymmetricTensor<DisplacementDim>>(
         variable_array[static_cast<int>(MPL::Variable::mechanical_strain)]);
     auto const& eps_m_prev = std::get<MPL::SymmetricTensor<DisplacementDim>>(
         variable_array_prev[static_cast<int>(
             MPL::Variable::mechanical_strain)]);
     auto const& sigma_prev = std::get<MPL::SymmetricTensor<DisplacementDim>>(
         variable_array_prev[static_cast<int>(MPL::Variable::stress)]);
     auto const& T = std::get<double>(
         variable_array_prev[static_cast<int>(MPL::Variable::temperature)]);
  
     KelvinMatrix C = getElasticTensor(t, x, T);
  
     KelvinVector sigma = sigma_prev + C * (eps_m - eps_m_prev);
  
     return {std::make_tuple(
         sigma,
         std::make_unique<
             typename MechanicsBase<DisplacementDim>::MaterialStateVariables>(),
         C)};
 }