MaterialLib::Fracture::LinearElasticIsotropic< DisplacementDim > Class Template Reference

Detailed Description

template<int DisplacementDim>
class MaterialLib::Fracture::LinearElasticIsotropic< DisplacementDim >

Definition at line 24 of file LinearElasticIsotropic.h.

#include <LinearElasticIsotropic.h>

Inheritance diagram for MaterialLib::Fracture::LinearElasticIsotropic< DisplacementDim >:

Collaboration diagram for MaterialLib::Fracture::LinearElasticIsotropic< DisplacementDim >:

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

Public Member Functions
std::unique_ptr< typename FractureModelBase< DisplacementDim >::MaterialStateVariables >	createMaterialStateVariables () override
	LinearElasticIsotropic (double const penalty_aperture_cutoff, bool const tension_cutoff, MaterialProperties material_properties)
void	computeConstitutiveRelation (double const t, ParameterLib::SpatialPosition const &x, double const aperture0, Eigen::Ref< Eigen::VectorXd const > sigma0, Eigen::Ref< Eigen::VectorXd const > w_prev, Eigen::Ref< Eigen::VectorXd const > w, Eigen::Ref< Eigen::VectorXd const > sigma_prev, Eigen::Ref< Eigen::VectorXd > sigma, Eigen::Ref< Eigen::MatrixXd > C, typename FractureModelBase< DisplacementDim >::MaterialStateVariables &material_state_variables) override
Public Member Functions inherited from MaterialLib::Fracture::FractureModelBase< DisplacementDim >
virtual	~FractureModelBase ()=default
virtual void	computeConstitutiveRelation (double const t, ParameterLib::SpatialPosition const &x, double const aperture0, Eigen::Ref< Eigen::VectorXd const > sigma0, Eigen::Ref< Eigen::VectorXd const > w_prev, Eigen::Ref< Eigen::VectorXd const > w, Eigen::Ref< Eigen::VectorXd const > sigma_prev, Eigen::Ref< Eigen::VectorXd > sigma, Eigen::Ref< Eigen::MatrixXd > C, MaterialStateVariables &material_state_variables)=0

Private Attributes
double const	_penalty_aperture_cutoff
bool const	_tension_cutoff
MaterialProperties	_mp

Constructor & Destructor Documentation

LinearElasticIsotropic()

template<int DisplacementDim>

MaterialLib::Fracture::LinearElasticIsotropic< DisplacementDim >::LinearElasticIsotropic ( double const penalty_aperture_cutoff, bool const tension_cutoff, MaterialProperties material_properties )

inlineexplicit

Definition at line 60 of file LinearElasticIsotropic.h.

        : _penalty_aperture_cutoff(penalty_aperture_cutoff),
          _tension_cutoff(tension_cutoff),
          _mp(std::move(material_properties))
    {
    }

Member Function Documentation

computeConstitutiveRelation()

template<int DisplacementDim>

void MaterialLib::Fracture::LinearElasticIsotropic< DisplacementDim >::computeConstitutiveRelation ( double const t, ParameterLib::SpatialPosition const & x, double const aperture0, Eigen::Ref< Eigen::VectorXd const > sigma0, Eigen::Ref< Eigen::VectorXd const > w_prev, Eigen::Ref< Eigen::VectorXd const > w, Eigen::Ref< Eigen::VectorXd const > sigma_prev, Eigen::Ref< Eigen::VectorXd > sigma, Eigen::Ref< Eigen::MatrixXd > C, typename FractureModelBase< DisplacementDim >::MaterialStateVariables & material_state_variables )

override

Computation of the constitutive relation for the linear elastic model.

Parameters

t	current time
x	current position in space
aperture0	initial fracture's aperture
sigma0	initial stress
w_prev	fracture displacement at previous time step
w	fracture displacement at current time step
sigma_prev	stress at previous time step
sigma	stress at current time step
C	tangent matrix for stress and fracture displacements
material_state_variables	material state variables

Definition at line 19 of file LinearElasticIsotropic.cpp.

{
    material_state_variables.reset();
 
    const int index_ns = DisplacementDim - 1;
    C.setZero();
    for (int i = 0; i < index_ns; i++)
    {
        C(i, i) = _mp.shear_stiffness(t, x)[0];
    }
 
    sigma.noalias() = C * w;
 
    double const aperture = w[index_ns] + aperture0;
 
    sigma.coeffRef(index_ns) =
        _mp.normal_stiffness(t, x)[0] * w[index_ns] *
        logPenalty(aperture0, aperture, _penalty_aperture_cutoff);
 
    C(index_ns, index_ns) =
        _mp.normal_stiffness(t, x)[0] *
        logPenaltyDerivative(aperture0, aperture, _penalty_aperture_cutoff);
 
    sigma.noalias() += sigma0;
 
    // correction for an opening fracture
    if (_tension_cutoff && sigma[index_ns] > 0)
    {
        C.setZero();
        sigma.setZero();
        material_state_variables.setTensileStress(true);
    }
}

References MaterialLib::Fracture::logPenalty(), MaterialLib::Fracture::logPenaltyDerivative(), MaterialLib::Fracture::FractureModelBase< DisplacementDim >::MaterialStateVariables::reset(), and MaterialLib::Fracture::FractureModelBase< DisplacementDim >::MaterialStateVariables::setTensileStress().

createMaterialStateVariables()

template<int DisplacementDim>

std::unique_ptr< typename FractureModelBase< DisplacementDim >::MaterialStateVariables > MaterialLib::Fracture::LinearElasticIsotropic< DisplacementDim >::createMaterialStateVariables ( )

inlineoverridevirtual

Polymorphic creator for MaterialStateVariables objects specific for a material model.

Implements MaterialLib::Fracture::FractureModelBase< DisplacementDim >.

Definition at line 52 of file LinearElasticIsotropic.h.

    {
        return std::unique_ptr<typename FractureModelBase<
            DisplacementDim>::MaterialStateVariables>{
            new MaterialStateVariables};
    }

Member Data Documentation

_mp

template<int DisplacementDim>

MaterialProperties MaterialLib::Fracture::LinearElasticIsotropic< DisplacementDim >::_mp

private

Definition at line 113 of file LinearElasticIsotropic.h.

_penalty_aperture_cutoff

template<int DisplacementDim>

double const MaterialLib::Fracture::LinearElasticIsotropic< DisplacementDim >::_penalty_aperture_cutoff

private

Compressive normal displacements above this value will not enter the computation of the normal stiffness modulus of the fracture.

Note

Setting this to the initial aperture value allows negative apertures.

Definition at line 107 of file LinearElasticIsotropic.h.

_tension_cutoff

template<int DisplacementDim>

bool const MaterialLib::Fracture::LinearElasticIsotropic< DisplacementDim >::_tension_cutoff

private

If set no resistance to open the fracture over the initial aperture is opposed.

Definition at line 111 of file LinearElasticIsotropic.h.

---

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

• MaterialLib/FractureModels/LinearElasticIsotropic.h
• MaterialLib/FractureModels/LinearElasticIsotropic.cpp