EshelbyTensor Class Reference

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

#include <EshelbyTensor.h>

Inheritance diagram for EshelbyTensor:

Public Member Functions
	EshelbyTensor (const InputParameters &parameters)
virtual void	initQpStatefulProperties () override
virtual void	computeQpProperties () override

Static Public Member Functions
static InputParameters	validParams ()

Protected Attributes
const std::string	_base_name
const MaterialProperty< Real > &	_sed
MaterialProperty< RankTwoTensor > &	_eshelby_tensor
const MaterialProperty< RankTwoTensor > &	_stress
const MaterialProperty< RankTwoTensor > &	_stress_old
std::vector< const VariableGradient * >	_grad_disp
MaterialProperty< RealVectorValue > &	_J_thermal_term_vec
const VariableGradient &	_grad_temp
const bool	_has_temp
const MaterialProperty< RankTwoTensor > *	_total_deigenstrain_dT

Detailed Description

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

Definition at line 27 of file EshelbyTensor.h.

Constructor & Destructor Documentation

EshelbyTensor()

EshelbyTensor::EshelbyTensor

(

const InputParameters &

parameters

)

Definition at line 36 of file EshelbyTensor.C.

  : DerivativeMaterialInterface<Material>(parameters),
    _base_name(isParamValid("base_name") ? getParam<std::string>("base_name") + "_" : ""),
    _sed(getMaterialPropertyByName<Real>(_base_name + "strain_energy_density")),
    _eshelby_tensor(declareProperty<RankTwoTensor>(_base_name + "Eshelby_tensor")),
    _stress(getMaterialProperty<RankTwoTensor>(_base_name + "stress")),
    _stress_old(getMaterialPropertyOld<RankTwoTensor>(_base_name + "stress")),
    _grad_disp(3),
    _J_thermal_term_vec(declareProperty<RealVectorValue>("J_thermal_term_vec")),
    _grad_temp(coupledGradient("temperature")),
    _has_temp(isCoupled("temperature")),
    _total_deigenstrain_dT(hasMaterialProperty<RankTwoTensor>("total_deigenstrain_dT")
                               ? &getMaterialProperty<RankTwoTensor>("total_deigenstrain_dT")
                               : nullptr)
{
  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;
 
  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.");
}

Member Function Documentation

computeQpProperties()

void EshelbyTensor::computeQpProperties ( )

overridevirtual

Definition at line 78 of file EshelbyTensor.C.

{
  RankTwoTensor F((*_grad_disp[0])[_qp],
                  (*_grad_disp[1])[_qp],
                  (*_grad_disp[2])[_qp]); // Deformation gradient
 
  RankTwoTensor H(F);
  F.addIa(1.0);
  Real detF = F.det();
  RankTwoTensor FinvT(F.inverse().transpose());
 
  // 1st Piola-Kirchoff Stress (P):
  RankTwoTensor P = detF * _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;
 
  if (_has_temp)
  {
    Real sigma_alpha = _stress[_qp].doubleContraction((*_total_deigenstrain_dT)[_qp]);
    _J_thermal_term_vec[_qp] = sigma_alpha * _grad_temp[_qp];
  }
  else
    _J_thermal_term_vec[_qp].zero();
}

initQpStatefulProperties()

void EshelbyTensor::initQpStatefulProperties ( )

overridevirtual

Definition at line 73 of file EshelbyTensor.C.

{
}

validParams()

InputParameters EshelbyTensor::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-Kirchoff 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.addCoupledVar("temperature", "Coupled temperature");
  return params;
}

Member Data Documentation

_base_name

const std::string EshelbyTensor::_base_name

protected

Definition at line 38 of file EshelbyTensor.h.

_eshelby_tensor

MaterialProperty<RankTwoTensor>& EshelbyTensor::_eshelby_tensor

protected

Definition at line 41 of file EshelbyTensor.h.

Referenced by computeQpProperties().

_grad_disp

std::vector<const VariableGradient *> EshelbyTensor::_grad_disp

protected

Definition at line 44 of file EshelbyTensor.h.

Referenced by computeQpProperties(), and EshelbyTensor().

_grad_temp

const VariableGradient& EshelbyTensor::_grad_temp

protected

Definition at line 47 of file EshelbyTensor.h.

Referenced by computeQpProperties().

_has_temp

const bool EshelbyTensor::_has_temp

protected

Definition at line 48 of file EshelbyTensor.h.

Referenced by computeQpProperties(), and EshelbyTensor().

_J_thermal_term_vec

MaterialProperty<RealVectorValue>& EshelbyTensor::_J_thermal_term_vec

protected

Definition at line 46 of file EshelbyTensor.h.

Referenced by computeQpProperties().

_sed

const MaterialProperty<Real>& EshelbyTensor::_sed

protected

Definition at line 40 of file EshelbyTensor.h.

Referenced by computeQpProperties().

_stress

const MaterialProperty<RankTwoTensor>& EshelbyTensor::_stress

protected

Definition at line 42 of file EshelbyTensor.h.

Referenced by computeQpProperties().

_stress_old

const MaterialProperty<RankTwoTensor>& EshelbyTensor::_stress_old

protected

Definition at line 43 of file EshelbyTensor.h.

_total_deigenstrain_dT

const MaterialProperty<RankTwoTensor>* EshelbyTensor::_total_deigenstrain_dT

protected

Definition at line 49 of file EshelbyTensor.h.

Referenced by computeQpProperties(), and EshelbyTensor().

---

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

• tensor_mechanics/include/materials/EshelbyTensor.h
• tensor_mechanics/src/materials/EshelbyTensor.C