ComputeVolumetricEigenstrain Class Reference

ComputeVolumetricEigenstrain computes an eigenstrain that is defined by a set of scalar material properties that summed together define the volumetric change. More...

#include <ComputeVolumetricEigenstrain.h>

Inheritance diagram for ComputeVolumetricEigenstrain:

Public Types
typedef DerivativeMaterialPropertyNameInterface::SymbolName	SymbolName

Public Member Functions
	ComputeVolumetricEigenstrain (const InputParameters &parameters)
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 Member Functions
virtual void	initialSetup ()
virtual void	computeQpEigenstrain ()

Protected Attributes
const unsigned int	_num_args
	number of variables the material depends on More...
const std::vector< MaterialPropertyName >	_volumetric_material_names
	Names of the material properties that define volumetric change. More...
std::vector< const MaterialProperty< Real > * >	_volumetric_materials
	The material properties that define volumetric change. More...
std::vector< std::vector< const MaterialProperty< Real > * > >	_dvolumetric_materials
	first derivatives of the volumetric materials with respect to the args More...
std::vector< std::vector< std::vector< const MaterialProperty< Real > * > > >	_d2volumetric_materials
	second derivatives of the volumetric materials with respect to the args More...
std::vector< MaterialProperty< RankTwoTensor > * >	_delastic_strain
	first derivatives of the elastic strain with respect to the args More...
std::vector< std::vector< MaterialProperty< RankTwoTensor > * > >	_d2elastic_strain
	second derivatives of the elastic strain with respect to the args More...

Detailed Description

ComputeVolumetricEigenstrain computes an eigenstrain that is defined by a set of scalar material properties that summed together define the volumetric change.

This also computes the derivatives of that eigenstrain with respect to a supplied set of variable dependencies.");

Definition at line 21 of file ComputeVolumetricEigenstrain.h.

Constructor & Destructor Documentation

ComputeVolumetricEigenstrain()

ComputeVolumetricEigenstrain::ComputeVolumetricEigenstrain

(

const InputParameters &

parameters

)

Definition at line 29 of file ComputeVolumetricEigenstrain.C.

  : DerivativeMaterialInterface<ComputeEigenstrainBase>(parameters),
    _num_args(coupledComponents("args")),
    _volumetric_material_names(getParam<std::vector<MaterialPropertyName>>("volumetric_materials")),
    _volumetric_materials(_volumetric_material_names.size()),
    _dvolumetric_materials(_volumetric_material_names.size()),
    _d2volumetric_materials(_volumetric_material_names.size()),
    _delastic_strain(_num_args),
    _d2elastic_strain(_num_args)
{
  for (unsigned int i = 0; i < _volumetric_material_names.size(); ++i)
    _volumetric_materials[i] = &getMaterialProperty<Real>(_volumetric_material_names[i]);

  // fetch prerequisite derivatives and build elastic_strain derivatives and cross-derivatives
  for (unsigned int i = 0; i < _volumetric_material_names.size(); ++i)
  {
    const MaterialPropertyName & vol_matl_name = _volumetric_material_names[i];
    _dvolumetric_materials[i].resize(_num_args);
    _d2volumetric_materials[i].resize(_num_args);
    for (unsigned int j = 0; j < _num_args; ++j)
    {
      const VariableName & jname = coupledName("args", j);
      _dvolumetric_materials[i][j] = &getMaterialPropertyDerivative<Real>(vol_matl_name, jname);
      _d2volumetric_materials[i][j].resize(_num_args);

      for (unsigned int k = j; k < _num_args; ++k)
      {
        const VariableName & kname = coupledName("args", k);
        _d2volumetric_materials[i][j][k] =
            &getMaterialPropertyDerivative<Real>("prefactor", jname, kname);
      }
    }
  }

  for (unsigned int j = 0; j < _num_args; ++j)
  {
    const VariableName & jname = coupledName("args", j);
    _delastic_strain[j] =
        &declarePropertyDerivative<RankTwoTensor>(_base_name + "elastic_strain", jname);
    _d2elastic_strain[j].resize(_num_args);

    for (unsigned int k = j; k < _num_args; ++k)
    {
      const VariableName & kname = coupledName("args", k);
      _d2elastic_strain[j][k] =
          &declarePropertyDerivative<RankTwoTensor>(_base_name + "elastic_strain", jname, kname);
    }
  }
}

Member Function Documentation

computeQpEigenstrain()

void ComputeVolumetricEigenstrain::computeQpEigenstrain ( )

protectedvirtual

Definition at line 106 of file ComputeVolumetricEigenstrain.C.

{
  Real volumetric_strain = 0;
  for (unsigned int i = 0; i < _volumetric_materials.size(); ++i)
    volumetric_strain += (*_volumetric_materials[i])[_qp];

  const Real eigenstrain_comp = computeVolumetricStrainComponent(volumetric_strain);
  _eigenstrain[_qp].zero();
  _eigenstrain[_qp].addIa(eigenstrain_comp);

  // TODO: Compute derivatives of the elastic strain wrt the variables specified in args
}

initialSetup()

void ComputeVolumetricEigenstrain::initialSetup ( )

protectedvirtual

Definition at line 80 of file ComputeVolumetricEigenstrain.C.

{
  for (auto vmn : _volumetric_material_names)
    validateCoupling<Real>(vmn);

  for (unsigned int i = 0; i < _num_args; ++i)
  {
    const VariableName & iname = coupledName("args", i);
    if (_fe_problem.isMatPropRequested(
            derivativePropertyNameFirst(_base_name + "elastic_strain", iname)))
      mooseError("Derivative of elastic_strain requested, but not yet implemented");
    else
      _delastic_strain[i] = nullptr;
    for (unsigned int j = 0; j < _num_args; ++j)
    {
      const VariableName & jname = coupledName("args", j);
      if (_fe_problem.isMatPropRequested(
              derivativePropertyNameSecond(_base_name + "elastic_strain", iname, jname)))
        mooseError("Second Derivative of elastic_strain requested, but not yet implemented");
      else
        _d2elastic_strain[i][j] = nullptr;
    }
  }
}

validParams()

InputParameters ComputeVolumetricEigenstrain::validParams ( )

static

Definition at line 16 of file ComputeVolumetricEigenstrain.C.

{
  InputParameters params = ComputeEigenstrainBase::validParams();
  params.addClassDescription("Computes an eigenstrain that is defined by a set of scalar material "
                             "properties that summed together define the volumetric change.  This "
                             "also computes the derivatives of that eigenstrain with respect to a "
                             "supplied set of variable dependencies.");
  params.addRequiredParam<std::vector<MaterialPropertyName>>(
      "volumetric_materials", "List of scalar material properties defining the volumetric change");
  params.addRequiredCoupledVar("args", "variable dependencies for the volumetric_expansion");
  return params;
}

Member Data Documentation

_d2elastic_strain

std::vector<std::vector<MaterialProperty<RankTwoTensor> *> > ComputeVolumetricEigenstrain::_d2elastic_strain

protected

second derivatives of the elastic strain with respect to the args

Definition at line 48 of file ComputeVolumetricEigenstrain.h.

Referenced by ComputeVolumetricEigenstrain(), and initialSetup().

_d2volumetric_materials

std::vector<std::vector<std::vector<const MaterialProperty<Real> *> > > ComputeVolumetricEigenstrain::_d2volumetric_materials

protected

second derivatives of the volumetric materials with respect to the args

Definition at line 43 of file ComputeVolumetricEigenstrain.h.

Referenced by ComputeVolumetricEigenstrain().

_delastic_strain

std::vector<MaterialProperty<RankTwoTensor> *> ComputeVolumetricEigenstrain::_delastic_strain

protected

first derivatives of the elastic strain with respect to the args

Definition at line 46 of file ComputeVolumetricEigenstrain.h.

Referenced by ComputeVolumetricEigenstrain(), and initialSetup().

_dvolumetric_materials

std::vector<std::vector<const MaterialProperty<Real> *> > ComputeVolumetricEigenstrain::_dvolumetric_materials

protected

first derivatives of the volumetric materials with respect to the args

Definition at line 41 of file ComputeVolumetricEigenstrain.h.

Referenced by ComputeVolumetricEigenstrain().

_num_args

const unsigned int ComputeVolumetricEigenstrain::_num_args

protected

number of variables the material depends on

Definition at line 33 of file ComputeVolumetricEigenstrain.h.

Referenced by ComputeVolumetricEigenstrain(), and initialSetup().

_volumetric_material_names

const std::vector<MaterialPropertyName> ComputeVolumetricEigenstrain::_volumetric_material_names

protected

Names of the material properties that define volumetric change.

Definition at line 36 of file ComputeVolumetricEigenstrain.h.

Referenced by ComputeVolumetricEigenstrain(), and initialSetup().

_volumetric_materials

std::vector<const MaterialProperty<Real> *> ComputeVolumetricEigenstrain::_volumetric_materials

protected

The material properties that define volumetric change.

Definition at line 38 of file ComputeVolumetricEigenstrain.h.

Referenced by computeQpEigenstrain(), and ComputeVolumetricEigenstrain().

---

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

• solid_mechanics/include/materials/ComputeVolumetricEigenstrain.h
• solid_mechanics/src/materials/ComputeVolumetricEigenstrain.C