ComputeElasticityTensorCP defines an elasticity tensor material object for crystal plasticity. More...

#include <ComputeElasticityTensorCP.h>

Inheritance diagram for ComputeElasticityTensorCP:

Public Member Functions
	ComputeElasticityTensorCP (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	assignEulerAngles ()

virtual void	computeQpProperties ()

void	issueGuarantee (const MaterialPropertyName &prop_name, Guarantee guarantee)

void	revokeGuarantee (const MaterialPropertyName &prop_name, Guarantee guarantee)

Protected Attributes
const ElementPropertyReadFile *	_read_prop_user_object
	Element property read user object Presently used to read Euler angles - see test. More...

MaterialProperty< RealVectorValue > &	_Euler_angles_mat_prop

MaterialProperty< RankTwoTensor > &	_crysrot
	Crystal Rotation Matrix. More...

RotationTensor	_R
	Rotation matrix. More...

RankFourTensor	_Cijkl
	Individual material information. 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

ComputeElasticityTensorCP defines an elasticity tensor material object for crystal plasticity.

Definition at line 25 of file ComputeElasticityTensorCP.h.

Constructor & Destructor Documentation

◆ ComputeElasticityTensorCP()

ComputeElasticityTensorCP::ComputeElasticityTensorCP ( const InputParameters & parameters )

Definition at line 29 of file ComputeElasticityTensorCP.C.

   : ComputeElasticityTensor(parameters),
     _read_prop_user_object(isParamValid("read_prop_user_object")
                                ? &getUserObject<ElementPropertyReadFile>("read_prop_user_object")
                                : nullptr),
     _Euler_angles_mat_prop(declareProperty<RealVectorValue>("Euler_angles")),
     _crysrot(declareProperty<RankTwoTensor>("crysrot")),
     _R(_Euler_angles)
 {
   // the base class guarantees constant in time, but in this derived class the
   // tensor will rotate over time once plastic deformation sets in
   revokeGuarantee(_elasticity_tensor_name, Guarantee::CONSTANT_IN_TIME);
  
   // the base class performs a passive rotation, but the crystal plasticity
   // materials use active rotation: recover unrotated _Cijkl here
   _Cijkl.rotate(_R.transpose());
 }

Member Function Documentation

◆ assignEulerAngles()

void ComputeElasticityTensorCP::assignEulerAngles ( )

protectedvirtual

Definition at line 48 of file ComputeElasticityTensorCP.C.

 {
   if (_read_prop_user_object)
   {
     _Euler_angles_mat_prop[_qp](0) = _read_prop_user_object->getData(_current_elem, 0);
     _Euler_angles_mat_prop[_qp](1) = _read_prop_user_object->getData(_current_elem, 1);
     _Euler_angles_mat_prop[_qp](2) = _read_prop_user_object->getData(_current_elem, 2);
   }
   else
     _Euler_angles_mat_prop[_qp] = _Euler_angles;
 }

Referenced by computeQpElasticityTensor().

◆ computeQpElasticityTensor()

void ComputeElasticityTensorCP::computeQpElasticityTensor ( )

protectedvirtual

Reimplemented from ComputeElasticityTensor.

Definition at line 61 of file ComputeElasticityTensorCP.C.

 {
   // Properties assigned at the beginning of every call to material calculation
   assignEulerAngles();
  
   _R.update(_Euler_angles_mat_prop[_qp]);
  
   _crysrot[_qp] = _R.transpose();
   _elasticity_tensor[_qp] = _Cijkl;
   _elasticity_tensor[_qp].rotate(_crysrot[_qp]);
 }

◆ 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().

◆ validParams()

InputParameters ComputeElasticityTensorCP::validParams ( )

static

Definition at line 18 of file ComputeElasticityTensorCP.C.

 {
   InputParameters params = ComputeElasticityTensor::validParams();
   params.addClassDescription("Compute an elasticity tensor for crystal plasticity.");
   params.addParam<UserObjectName>("read_prop_user_object",
                                   "The ElementReadPropertyFile "
                                   "GeneralUserObject to read element "
                                   "specific property values from file");
   return params;
 }