ComputeThermalExpansionEigenstrainTempl< is_ad > Class Template Reference

ComputeThermalExpansionEigenstrain computes an eigenstrain for thermal expansion with a constant expansion coefficient. More...

#include <ComputeThermalExpansionEigenstrain.h>

Inheritance diagram for ComputeThermalExpansionEigenstrainTempl< is_ad >:

Public Types
typedef DerivativeMaterialPropertyNameInterface::SymbolName	SymbolName

Public Member Functions
	ComputeThermalExpansionEigenstrainTempl (const InputParameters &parameters)
virtual void	subdomainSetup () final
	resize _temperature_buffer More...
virtual void	computeProperties () final
	update _temperature_buffer More...
virtual void	computeQpEigenstrain () 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 Member Functions
virtual ValueAndDerivative< is_ad >	computeThermalStrain () override
	computeThermalStrain must be overridden in derived classes. More...

Protected Attributes
const Real &	_thermal_expansion_coeff
const std::vector< ValueAndDerivative< is_ad > > &	_temperature
	Temperature to use in the eigenstrain calculation (current value if _use_old_temperature=false, old value if _use_old_temperature=true). More...
const bool	_use_old_temperature
	lag temperature variable More...
const VariableValue &	_temperature_old
	previous time step temperature More...
MaterialProperty< RankTwoTensor > *	_deigenstrain_dT
const VariableValue &	_stress_free_temperature

Detailed Description

template<bool is_ad>
class ComputeThermalExpansionEigenstrainTempl< is_ad >

ComputeThermalExpansionEigenstrain computes an eigenstrain for thermal expansion with a constant expansion coefficient.

Definition at line 19 of file ComputeThermalExpansionEigenstrain.h.

Constructor & Destructor Documentation

ComputeThermalExpansionEigenstrainTempl()

template<bool is_ad>

ComputeThermalExpansionEigenstrainTempl< is_ad >::ComputeThermalExpansionEigenstrainTempl

(

const InputParameters &

parameters

)

Definition at line 29 of file ComputeThermalExpansionEigenstrain.C.

  : ComputeThermalExpansionEigenstrainBaseTempl<is_ad>(parameters),
    _thermal_expansion_coeff(this->template getParam<Real>("thermal_expansion_coeff"))
{
}

Member Function Documentation

computeProperties()

template<bool is_ad>

void ComputeThermalExpansionEigenstrainBaseTempl< is_ad >::computeProperties ( )

finalvirtualinherited

update _temperature_buffer

Definition at line 71 of file ComputeThermalExpansionEigenstrainBase.C.

{
  // we need to convert the temperature variable to a ChainedReal in the is_ad == false case
  for (_qp = 0; _qp < this->_qrule->n_points(); ++_qp)
    if constexpr (is_ad)
      _temperature_buffer[_qp] =
          _use_old_temperature ? _temperature_old[_qp] : _temperature_prop[_qp];
    else
    {
      if (_use_old_temperature)
        _temperature_buffer[_qp] = {_temperature_old[_qp], 0};
      else
        _temperature_buffer[_qp] = {_temperature_prop[_qp], 1};
    }

  ComputeEigenstrainBaseTempl<is_ad>::computeProperties();
}

computeQpEigenstrain()

template<bool is_ad>

void ComputeThermalExpansionEigenstrainBaseTempl< is_ad >::computeQpEigenstrain ( )

overridevirtualinherited

Definition at line 91 of file ComputeThermalExpansionEigenstrainBase.C.

{
  _eigenstrain[_qp].zero();
  const auto thermal_strain = computeThermalStrain();

  if constexpr (is_ad)
  {
    _eigenstrain[_qp].addIa(thermal_strain);
    if (_mean_thermal_expansion_coefficient)
    {
      if (_temperature[_qp] == _stress_free_temperature[_qp])
        (*_mean_thermal_expansion_coefficient)[_qp] = 0.0;
      else
        (*_mean_thermal_expansion_coefficient)[_qp] = MetaPhysicL::raw_value(
            thermal_strain / (_temperature[_qp] - _stress_free_temperature[_qp]));
    }
  }
  else
  {
    _eigenstrain[_qp].addIa(thermal_strain.value());
    if (_mean_thermal_expansion_coefficient)
    {
      if (_temperature[_qp].value() == _stress_free_temperature[_qp])
        (*_mean_thermal_expansion_coefficient)[_qp] = 0.0;
      else
        (*_mean_thermal_expansion_coefficient)[_qp] =
            thermal_strain.value() / (_temperature[_qp].value() - _stress_free_temperature[_qp]);
    }
    if (_deigenstrain_dT)
    {
      (*_deigenstrain_dT)[_qp].zero();
      if (!_use_old_temperature)
        (*_deigenstrain_dT)[_qp].addIa(thermal_strain.derivatives());
    }
  }
}

computeThermalStrain()

template<bool is_ad>

ValueAndDerivative< is_ad > ComputeThermalExpansionEigenstrainTempl< is_ad >::computeThermalStrain ( )

overrideprotectedvirtual

computeThermalStrain must be overridden in derived classes.

The return type ValueAndDerivative<is_ad> contains the value for the thermal strain and its temperature derivative. Derived classes should use _temperature[_qp] to obtain the current temperature. In the is_ad == false case that member variable is agumented and will be of the type ChainedReal. I.e. even with is_ad == false a variant of forward mode automatic differentiation will be used internally to compute the thermal strain and no manual implementation of the temperature derivative is needed.

Returns

thermal strain and its derivative, where the thermal strain is the linear thermal strain ( L / L)

Implements ComputeThermalExpansionEigenstrainBaseTempl< is_ad >.

Definition at line 38 of file ComputeThermalExpansionEigenstrain.C.

{
  return _thermal_expansion_coeff * (_temperature[_qp] - _stress_free_temperature[_qp]);
}

subdomainSetup()

template<bool is_ad>

void ComputeThermalExpansionEigenstrainBaseTempl< is_ad >::subdomainSetup ( )

finalvirtualinherited

resize _temperature_buffer

Definition at line 59 of file ComputeThermalExpansionEigenstrainBase.C.

{
  // call parent class subdomain setup, which ultimately calls Material::subdomainSetup()
  ComputeEigenstrainBaseTempl<is_ad>::subdomainSetup();

  // make sure we have enouch space to hold the augmented temperature values
  const auto nqp = this->_fe_problem.getMaxQps();
  _temperature_buffer.resize(nqp);
}

validParams()

template<bool is_ad>

InputParameters ComputeThermalExpansionEigenstrainTempl< is_ad >::validParams ( )

static

Definition at line 17 of file ComputeThermalExpansionEigenstrain.C.

{
  InputParameters params = ComputeThermalExpansionEigenstrainBaseTempl<is_ad>::validParams();
  params.addClassDescription("Computes eigenstrain due to thermal expansion "
                             "with a constant coefficient");
  params.addRequiredParam<Real>("thermal_expansion_coeff", "Thermal expansion coefficient");
  params.declareControllable("thermal_expansion_coeff");

  return params;
}

Member Data Documentation

_deigenstrain_dT

template<bool is_ad>

MaterialProperty<RankTwoTensor>* ComputeThermalExpansionEigenstrainBaseTempl< is_ad >::_deigenstrain_dT

protectedinherited

Definition at line 72 of file ComputeThermalExpansionEigenstrainBase.h.

_stress_free_temperature

template<bool is_ad>

const VariableValue& ComputeThermalExpansionEigenstrainBaseTempl< is_ad >::_stress_free_temperature

protectedinherited

Definition at line 74 of file ComputeThermalExpansionEigenstrainBase.h.

_temperature

template<bool is_ad>

const std::vector<ValueAndDerivative<is_ad> >& ComputeThermalExpansionEigenstrainBaseTempl< is_ad >::_temperature

protectedinherited

Temperature to use in the eigenstrain calculation (current value if _use_old_temperature=false, old value if _use_old_temperature=true).

We use a const reference to a private member here to prevent derived classes from accidentally overwriting any values.

Definition at line 63 of file ComputeThermalExpansionEigenstrainBase.h.

_temperature_old

template<bool is_ad>

const VariableValue& ComputeThermalExpansionEigenstrainBaseTempl< is_ad >::_temperature_old

protectedinherited

previous time step temperature

Definition at line 69 of file ComputeThermalExpansionEigenstrainBase.h.

_thermal_expansion_coeff

template<bool is_ad>

const Real& ComputeThermalExpansionEigenstrainTempl< is_ad >::_thermal_expansion_coeff

protected

Definition at line 30 of file ComputeThermalExpansionEigenstrain.h.

_use_old_temperature

template<bool is_ad>

const bool ComputeThermalExpansionEigenstrainBaseTempl< is_ad >::_use_old_temperature

protectedinherited

lag temperature variable

Definition at line 66 of file ComputeThermalExpansionEigenstrainBase.h.

Referenced by ComputeInstantaneousThermalExpansionFunctionEigenstrainTempl< is_ad >::ComputeInstantaneousThermalExpansionFunctionEigenstrainTempl(), and ComputeThermalExpansionEigenstrainBaseTempl< is_ad >::ComputeThermalExpansionEigenstrainBaseTempl().

---

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

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