StressBasedChemicalPotential computes chemical potential based on stress and a direction tensor Forest et. More...

#include <StressBasedChemicalPotential.h>

Inheritance diagram for StressBasedChemicalPotential:

Public Member Functions
	StressBasedChemicalPotential (const InputParameters &parameters)

Static Public Member Functions
static InputParameters	validParams ()

Protected Member Functions
virtual void	initQpStatefulProperties ()

virtual void	computeQpProperties ()

Protected Attributes
MaterialProperty< Real > &	_chemical_potential

MaterialProperty< Real > *	_dchemical_potential

const MaterialProperty< RankTwoTensor > &	_stress_old

const MaterialProperty< RealTensorValue > &	_direction_tensor

const MaterialProperty< Real > &	_prefactor

const MaterialProperty< Real > *	_dprefactor_dc

bool	_has_coupled_c

Detailed Description

StressBasedChemicalPotential computes chemical potential based on stress and a direction tensor Forest et.

al. MSMSE 2015

Definition at line 26 of file StressBasedChemicalPotential.h.

Constructor & Destructor Documentation

◆ StressBasedChemicalPotential()

StressBasedChemicalPotential::StressBasedChemicalPotential ( const InputParameters & parameters )

Definition at line 31 of file StressBasedChemicalPotential.C.

   : DerivativeMaterialInterface<Material>(parameters),
     _chemical_potential(declareProperty<Real>(getParam<MaterialPropertyName>("property_name"))),
     _stress_old(getMaterialPropertyOld<RankTwoTensor>("stress_name")),
     _direction_tensor(getMaterialProperty<RealTensorValue>("direction_tensor_name")),
     _prefactor(getMaterialProperty<Real>("prefactor_name")),
     _has_coupled_c(isCoupled("c"))
 {
   if (_has_coupled_c)
   {
     _dchemical_potential = &declarePropertyDerivative<Real>(
         getParam<MaterialPropertyName>("property_name"), getVar("c", 0)->name());
     _dprefactor_dc = &getMaterialPropertyDerivative<Real>("prefactor_name", getVar("c", 0)->name());
   }
 }

Member Function Documentation

◆ computeQpProperties()

void StressBasedChemicalPotential::computeQpProperties ( )

protectedvirtual

Definition at line 57 of file StressBasedChemicalPotential.C.

 {
   RankTwoTensor direction_tensor_rank_two = _direction_tensor[_qp];
   _chemical_potential[_qp] =
       -_stress_old[_qp].doubleContraction(direction_tensor_rank_two) * _prefactor[_qp];
  
   if (_has_coupled_c)
     (*_dchemical_potential)[_qp] =
         -_stress_old[_qp].doubleContraction(direction_tensor_rank_two) * (*_dprefactor_dc)[_qp];
 }

◆ initQpStatefulProperties()

void StressBasedChemicalPotential::initQpStatefulProperties ( )

protectedvirtual

Definition at line 48 of file StressBasedChemicalPotential.C.

 {
   _chemical_potential[_qp] = 0.0;
  
   if (_has_coupled_c)
     (*_dchemical_potential)[_qp] = 0.0;
 }

◆ validParams()

InputParameters StressBasedChemicalPotential::validParams ( )

static

Definition at line 17 of file StressBasedChemicalPotential.C.

 {
   InputParameters params = Material::validParams();
   params.addClassDescription("Chemical potential from stress");
   params.addRequiredParam<MaterialPropertyName>("property_name",
                                                 "Name of stress based chemical potential");
   params.addRequiredParam<MaterialPropertyName>("stress_name", "Name of stress property variable");
   params.addRequiredParam<MaterialPropertyName>("direction_tensor_name",
                                                 "Name of direction tensor variable");
   params.addRequiredParam<MaterialPropertyName>("prefactor_name", "Name of prefactor variable");
   params.addCoupledVar("c", "Concentration variable");
   return params;
 }