MaterialLib::Solids::LinearElasticOrthotropic< DisplacementDim > Class Template Reference

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
using	KelvinMatrix
Public Types inherited from MaterialLib::Solids::MechanicsBase< DisplacementDim >
using	KelvinVector
using	KelvinMatrix

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 void	initializeInternalStateVariables (double const, ParameterLib::SpatialPosition const &, typename MechanicsBase< DisplacementDim >::MaterialStateVariables &) 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.
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

Initial value:

 
        MathLib::KelvinVector::KelvinMatrixType<DisplacementDim>

Definition at line 130 of file LinearElasticOrthotropic.h.

KelvinVector

template<int DisplacementDim>

using MaterialLib::Solids::LinearElasticOrthotropic< DisplacementDim >::KelvinVector

Initial value:

 
        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 55 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)>();
}

Referenced by MaterialLib::Solids::LinearElasticOrthotropic< DisplacementDim >::integrateStress().

getMaterialProperties()

template<int DisplacementDim>

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

inline

Definition at line 172 of file LinearElasticOrthotropic.h.

{ 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.mechanical_strain);
    auto const& eps_m_prev = std::get<MPL::SymmetricTensor<DisplacementDim>>(
        variable_array_prev.mechanical_strain);
    auto const& sigma_prev = std::get<MPL::SymmetricTensor<DisplacementDim>>(
        variable_array_prev.stress);
    auto const& T = variable_array_prev.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)};
}

References MaterialLib::Solids::LinearElasticOrthotropic< DisplacementDim >::getElasticTensor(), MaterialPropertyLib::VariableArray::mechanical_strain, MaterialPropertyLib::VariableArray::stress, and MaterialPropertyLib::VariableArray::temperature.

Member Data Documentation

_local_coordinate_system

template<int DisplacementDim>

std::optional<ParameterLib::CoordinateSystem> const& MaterialLib::Solids::LinearElasticOrthotropic< DisplacementDim >::_local_coordinate_system

protected

Definition at line 192 of file LinearElasticOrthotropic.h.

_mp

template<int DisplacementDim>

MaterialProperties MaterialLib::Solids::LinearElasticOrthotropic< DisplacementDim >::_mp

protected

Definition at line 190 of file LinearElasticOrthotropic.h.

Referenced by MaterialLib::Solids::LinearElasticOrthotropic< DisplacementDim >::getBulkModulus(), and MaterialLib::Solids::LinearElasticOrthotropic< DisplacementDim >::getMaterialProperties().

KelvinVectorSize

template<int DisplacementDim>

int const MaterialLib::Solids::LinearElasticOrthotropic< DisplacementDim >::KelvinVectorSize

static

Initial value:

=
        MathLib::KelvinVector::kelvin_vector_dimensions(DisplacementDim)

Definition at line 126 of file LinearElasticOrthotropic.h.

---

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

• MaterialLib/SolidModels/LinearElasticOrthotropic.h
• MaterialLib/SolidModels/LinearElasticOrthotropic.cpp