ComputeElasticityTensor defines an elasticity tensor material object as a function of concentration field. More...

#include <ComputeConcentrationDependentElasticityTensor.h>

Inheritance diagram for ComputeConcentrationDependentElasticityTensor:

Public Member Functions
	ComputeConcentrationDependentElasticityTensor (const InputParameters &parameters)

bool	hasGuarantee (const MaterialPropertyName &prop_name, Guarantee guarantee)

Static Public Member Functions
static InputParameters	validParams ()

Protected Member Functions
virtual void	computeQpElasticityTensor ()

virtual void	computeQpProperties ()

void	issueGuarantee (const MaterialPropertyName &prop_name, Guarantee guarantee)

void	revokeGuarantee (const MaterialPropertyName &prop_name, Guarantee guarantee)

Protected Attributes
RankFourTensor	_Cijkl0
	Elasticity tensor for phase with zero concentration. More...

RankFourTensor	_Cijkl1
	Elasticity tensor for phase with concentration 1. More...

const VariableValue &	_c
	Concentration variable. More...

VariableName	_c_name

MaterialProperty< RankFourTensor > &	_delasticity_tensor_dc
	Derivative of elasticity tensor with respect to concentration. More...

RealVectorValue	_Euler_angles

const std::string	_base_name

std::string	_elasticity_tensor_name

MaterialProperty< RankFourTensor > &	_elasticity_tensor

MaterialProperty< Real > &	_effective_stiffness

const Function *const	_prefactor_function
	prefactor function to multiply the elasticity tensor with More...

Private Attributes
std::map< MaterialPropertyName, std::set< Guarantee > >	_guarantees

Detailed Description

ComputeElasticityTensor defines an elasticity tensor material object as a function of concentration field.

Definition at line 23 of file ComputeConcentrationDependentElasticityTensor.h.

Constructor & Destructor Documentation

◆ ComputeConcentrationDependentElasticityTensor()

ComputeConcentrationDependentElasticityTensor::ComputeConcentrationDependentElasticityTensor ( const InputParameters & parameters )

Definition at line 34 of file ComputeConcentrationDependentElasticityTensor.C.

   : ComputeRotatedElasticityTensorBase(parameters),
     _Cijkl0(getParam<std::vector<Real>>("C0_ijkl"),
             (RankFourTensor::FillMethod)(int)getParam<MooseEnum>("fill_method0")),
     _Cijkl1(getParam<std::vector<Real>>("C1_ijkl"),
             (RankFourTensor::FillMethod)(int)getParam<MooseEnum>("fill_method1")),
     _c(coupledValue("c")),
     _c_name(getVar("c", 0)->name()),
     _delasticity_tensor_dc(
         declarePropertyDerivative<RankFourTensor>(_elasticity_tensor_name, _c_name))
 {
   // Define a rotation according to Euler angle parameters
   RotationTensor R(_Euler_angles); // R type: RealTensorValue
  
   // Rotate tensors
   _Cijkl0.rotate(R);
   _Cijkl1.rotate(R);
 }

Member Function Documentation

◆ computeQpElasticityTensor()

void ComputeConcentrationDependentElasticityTensor::computeQpElasticityTensor ( )

protectedvirtual

Implements ComputeElasticityTensorBase.

Definition at line 55 of file ComputeConcentrationDependentElasticityTensor.C.

 {
   // Assign elasticity tensor at a given quad point
   _elasticity_tensor[_qp] = _Cijkl0 + (_Cijkl1 - _Cijkl0) * _c[_qp];
   // Define derivative of elasticity tensor with respect to concentration.
   _delasticity_tensor_dc[_qp] = (_Cijkl1 - _Cijkl0);
 }

◆ computeQpProperties()

void ComputeElasticityTensorBase::computeQpProperties ( )

protectedvirtualinherited

Definition at line 43 of file ComputeElasticityTensorBase.C.

 {
   _effective_stiffness[_qp] = 0; // Currently overriden by ComputeIsotropicElasticityTensor
   computeQpElasticityTensor();
  
   // Multiply by prefactor
   if (_prefactor_function)
   {
     _elasticity_tensor[_qp] *= _prefactor_function->value(_t, _q_point[_qp]);
     _effective_stiffness[_qp] *= std::sqrt(_prefactor_function->value(_t, _q_point[_qp]));
   }
 }

◆ hasGuarantee()

bool GuaranteeProvider::hasGuarantee	(	const MaterialPropertyName &	prop_name,
		Guarantee	guarantee
	)

inherited

Definition at line 16 of file GuaranteeProvider.C.

 {
   auto it = _guarantees.find(prop_name);
   if (it == _guarantees.end())
     return false;
  
   auto it2 = it->second.find(guarantee);
   return it2 != it->second.end();
 }

◆ issueGuarantee()

void GuaranteeProvider::issueGuarantee	(	const MaterialPropertyName &	prop_name,
		Guarantee	guarantee
	)

protectedinherited

Definition at line 27 of file GuaranteeProvider.C.

 {
   // intentional insertion
   _guarantees[prop_name].insert(guarantee);
 }

Referenced by ADComputeVariableIsotropicElasticityTensor< compute_stage >::ADComputeVariableIsotropicElasticityTensor(), ComputeCosseratElasticityTensor::ComputeCosseratElasticityTensor(), ComputeElasticityTensor::ComputeElasticityTensor(), ComputeIsotropicElasticityTensor::ComputeIsotropicElasticityTensor(), ComputeLayeredCosseratElasticityTensor::ComputeLayeredCosseratElasticityTensor(), ComputeVariableIsotropicElasticityTensor::ComputeVariableIsotropicElasticityTensor(), GeneralizedKelvinVoigtModel::GeneralizedKelvinVoigtModel(), and GeneralizedMaxwellModel::GeneralizedMaxwellModel().

◆ revokeGuarantee()

void GuaranteeProvider::revokeGuarantee	(	const MaterialPropertyName &	prop_name,
		Guarantee	guarantee
	)

protectedinherited

Definition at line 34 of file GuaranteeProvider.C.

 {
   auto it = _guarantees.find(prop_name);
   if (it != _guarantees.end())
     it->second.erase(guarantee);
 }

Referenced by ComputeElasticityTensorCP::ComputeElasticityTensorCP().

◆ validParams()

InputParameters ComputeConcentrationDependentElasticityTensor::validParams ( )

static

Definition at line 18 of file ComputeConcentrationDependentElasticityTensor.C.

 {
   InputParameters params = ComputeRotatedElasticityTensorBase::validParams();
   params.addClassDescription("Compute concentration dependent elasticity tensor.");
   params.addRequiredParam<std::vector<Real>>("C0_ijkl",
                                              "Stiffness tensor for zero concentration phase");
   params.addRequiredParam<std::vector<Real>>("C1_ijkl",
                                              "Stiffness tensor for phase having concentration 1.0");
   params.addParam<MooseEnum>(
       "fill_method0", RankFourTensor::fillMethodEnum() = "symmetric9", "The fill method");
   params.addParam<MooseEnum>(
       "fill_method1", RankFourTensor::fillMethodEnum() = "symmetric9", "The fill method");
   params.addRequiredCoupledVar("c", "Concentration");
   return params;
 }