EshelbyTensor defines a strain increment and rotation increment, for finite strains. More...

#include <EshelbyTensor.h>

Inheritance diagram for EshelbyTensorTempl< is_ad >:

Public Types
typedef DerivativeMaterialPropertyNameInterface::SymbolName	SymbolName

Public Member Functions
	EshelbyTensorTempl (const InputParameters &parameters)

virtual void	initialSetup () override

virtual void	initQpStatefulProperties () override

virtual void	computeQpProperties () override

const GenericMaterialProperty< U, is_ad > &	getDefaultMaterialProperty (const std::string &name)

const GenericMaterialProperty< U, is_ad > &	getDefaultMaterialPropertyByName (const std::string &name)

void	validateDerivativeMaterialPropertyBase (const std::string &base)

const MaterialPropertyName	derivativePropertyName (const MaterialPropertyName &base, const std::vector< SymbolName > &c) const

const MaterialPropertyName	derivativePropertyNameFirst (const MaterialPropertyName &base, const SymbolName &c1) const

const MaterialPropertyName	derivativePropertyNameSecond (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2) const

const MaterialPropertyName	derivativePropertyNameThird (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2, const SymbolName &c3) const

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const std::vector< VariableName > &c)

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const std::vector< VariableName > &c)

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const std::vector< VariableName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, unsigned int v2, unsigned int v3=libMesh::invalid_uint)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, unsigned int v1, unsigned int v2=libMesh::invalid_uint, unsigned int v3=libMesh::invalid_uint)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const std::vector< VariableName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, unsigned int v2, unsigned int v3=libMesh::invalid_uint)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, unsigned int v1, unsigned int v2=libMesh::invalid_uint, unsigned int v3=libMesh::invalid_uint)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< VariableName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< SymbolName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< VariableName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< SymbolName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

void	validateCoupling (const MaterialPropertyName &base, const std::vector< VariableName > &c, bool validate_aux=true)

void	validateCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")

void	validateCoupling (const MaterialPropertyName &base, const std::vector< VariableName > &c, bool validate_aux=true)

void	validateCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")

void	validateNonlinearCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")

void	validateNonlinearCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")

const MaterialPropertyName	propertyName (const MaterialPropertyName &base, const std::vector< SymbolName > &c) const

const MaterialPropertyName	propertyName (const MaterialPropertyName &base, const std::vector< SymbolName > &c) const

const MaterialPropertyName	propertyNameFirst (const MaterialPropertyName &base, const SymbolName &c1) const

const MaterialPropertyName	propertyNameFirst (const MaterialPropertyName &base, const SymbolName &c1) const

const MaterialPropertyName	propertyNameSecond (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2) const

const MaterialPropertyName	propertyNameSecond (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2) const

const MaterialPropertyName	propertyNameThird (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2, const SymbolName &c3) const

const MaterialPropertyName	propertyNameThird (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2, const SymbolName &c3) const

Static Public Member Functions
static InputParameters	validParams ()

Protected Attributes
const std::string	_base_name
	Base name of the material system.

const bool	_compute_dissipation
	Whether to also compute Eshelby tensor's dissipation for C(t) integral.

const MaterialProperty< Real > &	_sed

const MaterialProperty< Real > *	_serd

MaterialProperty< RankTwoTensor > &	_eshelby_tensor

MaterialProperty< RankTwoTensor > *	_eshelby_tensor_dissipation

const GenericMaterialProperty< RankTwoTensor, is_ad > &	_stress
	The stress tensor.

const MaterialProperty< RankTwoTensor > &	_stress_old
	The old stress tensor.

std::vector< const VariableGradient * >	_grad_disp

std::vector< const VariableGradient * >	_grad_disp_old

MaterialProperty< RealVectorValue > &	_J_thermal_term_vec

const VariableGradient &	_grad_temp

const bool	_has_temp

const OptionalMaterialProperty< RankTwoTensor > &	_total_deigenstrain_dT

Private Member Functions
bool	haveMaterialProperty (const std::string &prop_name)

std::vector< VariableName >	buildVariableVector (const VariableName &c1, const VariableName &c2, const VariableName &c3)

void	validateCouplingHelper (const MaterialPropertyName &base, const std::vector< VariableName > &c, const System &system, std::vector< VariableName > &missing)

bool	isNotObjectVariable (const VariableName &name)

Private Attributes
FEProblemBase &	_dmi_fe_problem

Detailed Description

template<bool is_ad>
class EshelbyTensorTempl< is_ad >

EshelbyTensor defines a strain increment and rotation increment, for finite strains.

Definition at line 21 of file EshelbyTensor.h.

Constructor & Destructor Documentation

◆ EshelbyTensorTempl()

template<bool is_ad>

EshelbyTensorTempl< is_ad >::EshelbyTensorTempl ( const InputParameters & parameters )

Definition at line 42 of file EshelbyTensor.C.

  : DerivativeMaterialInterface<Material>(parameters),
    _base_name(isParamValid("base_name") ? getParam<std::string>("base_name") + "_" : ""),
    _compute_dissipation(getParam<bool>("compute_dissipation")),
    _sed(getMaterialPropertyByName<Real>(_base_name + "strain_energy_density")),
    _serd(_compute_dissipation
              ? &getMaterialPropertyByName<Real>(_base_name + "strain_energy_rate_density")
              : nullptr),
    _eshelby_tensor(declareProperty<RankTwoTensor>(_base_name + "Eshelby_tensor")),
    _eshelby_tensor_dissipation(
        _compute_dissipation
            ? &declareProperty<RankTwoTensor>(_base_name + "Eshelby_tensor_dissipation")
            : nullptr),
    _stress(getGenericMaterialProperty<RankTwoTensor, is_ad>(_base_name + "stress")),
    _stress_old(getMaterialPropertyOld<RankTwoTensor>(_base_name + "stress")),
    _grad_disp(3),
    _grad_disp_old(3),
    _J_thermal_term_vec(declareProperty<RealVectorValue>("J_thermal_term_vec")),
    _grad_temp(coupledGradient("temperature")),
    _has_temp(isCoupled("temperature")),
    _total_deigenstrain_dT(getOptionalMaterialProperty<RankTwoTensor>("total_deigenstrain_dT"))
{
  unsigned int ndisp = coupledComponents("displacements");
 
  // Checking for consistency between mesh size and length of the provided displacements vector
  if (ndisp != _mesh.dimension())
    mooseError(
        "The number of variables supplied in 'displacements' must match the mesh dimension.");
 
  // fetch coupled gradients
  for (unsigned int i = 0; i < ndisp; ++i)
    _grad_disp[i] = &coupledGradient("displacements", i);
 
  // set unused dimensions to zero
  for (unsigned i = ndisp; i < 3; ++i)
    _grad_disp[i] = &_grad_zero;
 
  // Need previous step's displacements to compute deformation gradient time rate
  if (_compute_dissipation)
  {
    // fetch coupled gradients previous step
    for (unsigned int i = 0; i < ndisp; ++i)
      _grad_disp_old[i] = &coupledGradientOld("displacements", i);
 
    // set unused dimensions to zero
    for (unsigned i = ndisp; i < 3; ++i)
      _grad_disp_old[i] = &_grad_zero;
  }
}

Member Function Documentation

◆ computeQpProperties()

template<bool is_ad>

void EshelbyTensorTempl< is_ad >::computeQpProperties ( )

overridevirtual

Definition at line 110 of file EshelbyTensor.C.

{
  // Deformation gradient
  auto F = RankTwoTensor::initializeFromRows(
      (*_grad_disp[0])[_qp], (*_grad_disp[1])[_qp], (*_grad_disp[2])[_qp]);
 
  RankTwoTensor H(F);
  F.addIa(1.0);
  Real detF = F.det();
  RankTwoTensor FinvT(F.inverse().transpose());
 
  // 1st Piola-Kirchoff Stress (P):
  RankTwoTensor P = detF * MetaPhysicL::raw_value(_stress[_qp]) * FinvT;
 
  // HTP = H^T * P = H^T * detF * sigma * FinvT;
  RankTwoTensor HTP = H.transpose() * P;
 
  RankTwoTensor WI = RankTwoTensor(RankTwoTensor::initIdentity);
  WI *= (_sed[_qp] * detF);
 
  _eshelby_tensor[_qp] = WI - HTP;
 
  // Compute deformation gradient rate
  if (_compute_dissipation == true)
  {
    auto H_old = RankTwoTensor::initializeFromRows(
        (*_grad_disp_old[0])[_qp], (*_grad_disp_old[1])[_qp], (*_grad_disp_old[2])[_qp]);
 
    RankTwoTensor Wdot = RankTwoTensor(RankTwoTensor::initIdentity);
    Wdot *= ((*_serd)[_qp] * detF);
 
    // F_dot = (F - F_old)/dt
    RankTwoTensor F_dot = (H - H_old) / _dt;
 
    // FdotTP = Fdot^T * P = Fdot^T * detF * sigma * FinvT;
    RankTwoTensor FdotTP = F_dot.transpose() * P;
 
    (*_eshelby_tensor_dissipation)[_qp] = Wdot - FdotTP;
  }
 
  if (_has_temp)
  {
    const Real sigma_alpha =
        MetaPhysicL::raw_value(_stress[_qp]).doubleContraction(_total_deigenstrain_dT[_qp]);
    _J_thermal_term_vec[_qp] = sigma_alpha * _grad_temp[_qp];
  }
  else
    _J_thermal_term_vec[_qp].zero();
}

◆ initialSetup()

template<bool is_ad>

void EshelbyTensorTempl< is_ad >::initialSetup ( )

overridevirtual

Definition at line 94 of file EshelbyTensor.C.

{
  if (_has_temp && !_total_deigenstrain_dT)
    mooseError("EshelbyTensor Error: To include thermal strain term in Fracture integral "
               "calculation, must both couple temperature in DomainIntegral block and compute "
               "total_deigenstrain_dT using ThermalFractureIntegral material model.");
}

◆ initQpStatefulProperties()

template<bool is_ad>

void EshelbyTensorTempl< is_ad >::initQpStatefulProperties ( )

overridevirtual

Definition at line 104 of file EshelbyTensor.C.

105{

106}

◆ validParams()

template<bool is_ad>

InputParameters EshelbyTensorTempl< is_ad >::validParams ( )

static

Definition at line 19 of file EshelbyTensor.C.

{
  InputParameters params = Material::validParams();
  params.addClassDescription("Computes the Eshelby tensor as a function of "
                             "strain energy density and the first "
                             "Piola-Kirchhoff stress");
  params.addRequiredCoupledVar(
      "displacements",
      "The displacements appropriate for the simulation geometry and coordinate system");
  params.addParam<std::string>("base_name",
                               "Optional parameter that allows the user to define "
                               "multiple mechanics material systems on the same "
                               "block, i.e. for multiple phases");
  params.addParam<bool>(
      "compute_dissipation",
      false,
      "Whether to compute Eshelby tensor's dissipation (or rate of change). This tensor"
      "yields the increase in dissipation per unit crack advanced");
  params.addCoupledVar("temperature", "Coupled temperature");
  return params;
}

Member Data Documentation

◆ _base_name

template<bool is_ad>

const std::string EshelbyTensorTempl< is_ad >::_base_name

protected

Base name of the material system.

Definition at line 34 of file EshelbyTensor.h.

◆ _compute_dissipation

template<bool is_ad>

const bool EshelbyTensorTempl< is_ad >::_compute_dissipation

protected

Whether to also compute Eshelby tensor's dissipation for C(t) integral.

Definition at line 37 of file EshelbyTensor.h.

Referenced by EshelbyTensorTempl< is_ad >::EshelbyTensorTempl().

◆ _eshelby_tensor

template<bool is_ad>

MaterialProperty<RankTwoTensor>& EshelbyTensorTempl< is_ad >::_eshelby_tensor

protected

Definition at line 42 of file EshelbyTensor.h.

◆ _eshelby_tensor_dissipation

template<bool is_ad>

MaterialProperty<RankTwoTensor>* EshelbyTensorTempl< is_ad >::_eshelby_tensor_dissipation

protected

Definition at line 43 of file EshelbyTensor.h.

◆ _grad_disp

template<bool is_ad>

std::vector<const VariableGradient *> EshelbyTensorTempl< is_ad >::_grad_disp

protected

Definition at line 51 of file EshelbyTensor.h.

Referenced by EshelbyTensorTempl< is_ad >::EshelbyTensorTempl().

◆ _grad_disp_old

template<bool is_ad>

std::vector<const VariableGradient *> EshelbyTensorTempl< is_ad >::_grad_disp_old

protected

Definition at line 52 of file EshelbyTensor.h.

Referenced by EshelbyTensorTempl< is_ad >::EshelbyTensorTempl().

◆ _grad_temp

template<bool is_ad>

const VariableGradient& EshelbyTensorTempl< is_ad >::_grad_temp

protected

Definition at line 55 of file EshelbyTensor.h.

◆ _has_temp

template<bool is_ad>

const bool EshelbyTensorTempl< is_ad >::_has_temp

protected

Definition at line 56 of file EshelbyTensor.h.

◆ _J_thermal_term_vec

template<bool is_ad>

MaterialProperty<RealVectorValue>& EshelbyTensorTempl< is_ad >::_J_thermal_term_vec

protected

Definition at line 54 of file EshelbyTensor.h.

◆ _sed

template<bool is_ad>

const MaterialProperty<Real>& EshelbyTensorTempl< is_ad >::_sed

protected

Definition at line 39 of file EshelbyTensor.h.

◆ _serd

template<bool is_ad>

const MaterialProperty<Real>* EshelbyTensorTempl< is_ad >::_serd

protected

Definition at line 40 of file EshelbyTensor.h.

◆ _stress

template<bool is_ad>

const GenericMaterialProperty<RankTwoTensor, is_ad>& EshelbyTensorTempl< is_ad >::_stress

protected

The stress tensor.

Definition at line 46 of file EshelbyTensor.h.

◆ _stress_old

template<bool is_ad>

const MaterialProperty<RankTwoTensor>& EshelbyTensorTempl< is_ad >::_stress_old

protected

The old stress tensor.

Definition at line 49 of file EshelbyTensor.h.

◆ _total_deigenstrain_dT

template<bool is_ad>

const OptionalMaterialProperty<RankTwoTensor>& EshelbyTensorTempl< is_ad >::_total_deigenstrain_dT

protected

Definition at line 57 of file EshelbyTensor.h.

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

solid_mechanics/include/materials/EshelbyTensor.h
solid_mechanics/src/materials/EshelbyTensor.C

Public Types

Public Member Functions

Static Public Member Functions

Protected Attributes

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ EshelbyTensorTempl()

Member Function Documentation

◆ computeQpProperties()

◆ initialSetup()

◆ initQpStatefulProperties()

◆ validParams()

Member Data Documentation

◆ _base_name

◆ _compute_dissipation

◆ _eshelby_tensor

◆ _eshelby_tensor_dissipation

◆ _grad_disp

◆ _grad_disp_old

◆ _grad_temp

◆ _has_temp

◆ _J_thermal_term_vec

◆ _sed

◆ _serd

◆ _stress

◆ _stress_old

◆ _total_deigenstrain_dT