ComputeDilatationThermalExpansionFunctionEigenstrain computes an eigenstrain for thermal expansion from an dilatation function. More...

#include <ComputeDilatationThermalExpansionFunctionEigenstrain.h>

Inheritance diagram for ComputeDilatationThermalExpansionFunctionEigenstrainTempl< is_ad >:

[legend]

Public Types
typedef DerivativeMaterialPropertyNameInterface::SymbolName	SymbolName

Public Member Functions
	ComputeDilatationThermalExpansionFunctionEigenstrainTempl (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 >	computeDilatation (const ValueAndDerivative< is_ad > &temperature) override
	Compute the fractional linear dilatation due to thermal expansion delta L / L along with its temperature derivative. More...

virtual ValueAndDerivative< is_ad >	computeThermalStrain () override
	computeThermalStrain must be overridden in derived classes. More...

Protected Attributes
const Function &	_dilatation_function
	Dilatation function. More...

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 ComputeDilatationThermalExpansionFunctionEigenstrainTempl< is_ad >

ComputeDilatationThermalExpansionFunctionEigenstrain computes an eigenstrain for thermal expansion from an dilatation function.

Definition at line 19 of file ComputeDilatationThermalExpansionFunctionEigenstrain.h.

Constructor & Destructor Documentation

◆ ComputeDilatationThermalExpansionFunctionEigenstrainTempl()

template<bool is_ad>

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

Definition at line 32 of file ComputeDilatationThermalExpansionFunctionEigenstrain.C.

   : ComputeDilatationThermalExpansionEigenstrainBaseTempl<is_ad>(parameters),
     _dilatation_function(this->getFunction("dilatation_function"))
 {
 }

Member Function Documentation

◆ computeDilatation()

template<bool is_ad>

ValueAndDerivative< is_ad > ComputeDilatationThermalExpansionFunctionEigenstrainTempl< is_ad >::computeDilatation ( const ValueAndDerivative< is_ad > & temperature )

overrideprotectedvirtual

Compute the fractional linear dilatation due to thermal expansion delta L / L along with its temperature derivative.

Parameters

temperature current temperature

Returns: fractional linear dilatation with temperature derivative

Implements ComputeDilatationThermalExpansionEigenstrainBaseTempl< is_ad >.

Definition at line 40 of file ComputeDilatationThermalExpansionFunctionEigenstrain.C.

 {
   // we need these two branches because we cannot yet evaluate Functions with ChainedReals
   if constexpr (is_ad)
     return _dilatation_function.value(temperature);
   else
     return {_dilatation_function.value(temperature.value()),
             _dilatation_function.timeDerivative(temperature.value()) * temperature.derivatives()};
 }

◆ 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 > ComputeDilatationThermalExpansionEigenstrainBaseTempl< is_ad >::computeThermalStrain ( )

overrideprotectedvirtualinherited

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 28 of file ComputeDilatationThermalExpansionEigenstrainBase.C.

 {
   const auto stress_free = computeDilatation(this->_stress_free_temperature[_qp]);
   const auto current = computeDilatation(this->_temperature[_qp]);
 
   // in non-AD mode the T derivative of the stress_free term needs get dropped.
   // We assume _stress_free_temperature does not depend on T. In AD mode this is automatic.
   return current - (is_ad ? stress_free : stress_free.value());
 }

◆ 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 ComputeDilatationThermalExpansionFunctionEigenstrainTempl< is_ad >::validParams ( )

static

Definition at line 19 of file ComputeDilatationThermalExpansionFunctionEigenstrain.C.

 {
   InputParameters params = ComputeDilatationThermalExpansionEigenstrainBase::validParams();
   params.addClassDescription("Computes eigenstrain due to thermal expansion using a function that "
                              "describes the total dilatation as a function of temperature");
   params.addRequiredParam<FunctionName>(
       "dilatation_function",
       "Function describing the thermal dilatation as a function of temperature");
   return params;
 }