Compute an evolving elasticity tensor coupled to a grain growth phase field model. More...

#include <ComputePolycrystalElasticityTensor.h>

Inheritance diagram for ComputePolycrystalElasticityTensor:

Public Member Functions
	ComputePolycrystalElasticityTensor (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
Real	_length_scale

Real	_pressure_scale

const GrainDataTracker< RankFourTensor > &	_grain_tracker
	Grain tracker object. More...

const unsigned int	_op_num
	Number of order parameters. More...

std::vector< const VariableValue * >	_vals
	Order parameters. More...

std::vector< MaterialProperty< RankFourTensor > * >	_D_elastic_tensor
	vector of elasticity tensor material properties More...

const Real	_JtoeV
	Conversion factor from J to eV. More...

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

Compute an evolving elasticity tensor coupled to a grain growth phase field model.

Definition at line 25 of file ComputePolycrystalElasticityTensor.h.

Constructor & Destructor Documentation

◆ ComputePolycrystalElasticityTensor()

ComputePolycrystalElasticityTensor::ComputePolycrystalElasticityTensor ( const InputParameters & parameters )

Definition at line 31 of file ComputePolycrystalElasticityTensor.C.

   : ComputeElasticityTensorBase(parameters),
     _length_scale(getParam<Real>("length_scale")),
     _pressure_scale(getParam<Real>("pressure_scale")),
     _grain_tracker(getUserObject<GrainDataTracker<RankFourTensor>>("grain_tracker")),
     _op_num(coupledComponents("v")),
     _vals(_op_num),
     _D_elastic_tensor(_op_num),
     _JtoeV(6.24150974e18)
 {
   // Loop over variables (ops)
   for (MooseIndex(_op_num) op_index = 0; op_index < _op_num; ++op_index)
   {
     // Initialize variables
     _vals[op_index] = &coupledValue("v", op_index);
  
     // declare elasticity tensor derivative properties
     _D_elastic_tensor[op_index] = &declarePropertyDerivative<RankFourTensor>(
         _elasticity_tensor_name, getVar("v", op_index)->name());
   }
 }

Member Function Documentation

◆ computeQpElasticityTensor()

void ComputePolycrystalElasticityTensor::computeQpElasticityTensor ( )

protectedvirtual

Implements ComputeElasticityTensorBase.

Definition at line 55 of file ComputePolycrystalElasticityTensor.C.

 {
   // Get list of active order parameters from grain tracker
   const auto & op_to_grains = _grain_tracker.getVarToFeatureVector(_current_elem->id());
  
   // Calculate elasticity tensor
   _elasticity_tensor[_qp].zero();
   Real sum_h = 0.0;
   for (MooseIndex(op_to_grains) op_index = 0; op_index < op_to_grains.size(); ++op_index)
   {
     auto grain_id = op_to_grains[op_index];
     if (grain_id == FeatureFloodCount::invalid_id)
       continue;
  
     // Interpolation factor for elasticity tensors
     Real h = (1.0 + std::sin(libMesh::pi * ((*_vals[op_index])[_qp] - 0.5))) / 2.0;
  
     // Sum all rotated elasticity tensors
     _elasticity_tensor[_qp] += _grain_tracker.getData(grain_id) * h;
     sum_h += h;
   }
  
   const Real tol = 1.0e-10;
   sum_h = std::max(sum_h, tol);
   _elasticity_tensor[_qp] /= sum_h;
  
   // Calculate elasticity tensor derivative: Cderiv = dhdopi/sum_h * (Cop - _Cijkl)
   for (MooseIndex(_op_num) op_index = 0; op_index < _op_num; ++op_index)
     (*_D_elastic_tensor[op_index])[_qp].zero();
  
   for (MooseIndex(op_to_grains) op_index = 0; op_index < op_to_grains.size(); ++op_index)
   {
     auto grain_id = op_to_grains[op_index];
     if (grain_id == FeatureFloodCount::invalid_id)
       continue;
  
     Real dhdopi = libMesh::pi * std::cos(libMesh::pi * ((*_vals[op_index])[_qp] - 0.5)) / 2.0;
     RankFourTensor & C_deriv = (*_D_elastic_tensor[op_index])[_qp];
  
     C_deriv = (_grain_tracker.getData(grain_id) - _elasticity_tensor[_qp]) * dhdopi / sum_h;
  
     // Convert from XPa to eV/(xm)^3, where X is pressure scale and x is length scale;
     C_deriv *= _JtoeV * (_length_scale * _length_scale * _length_scale) * _pressure_scale;
   }
 }

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

staticinherited

Definition at line 16 of file ComputeElasticityTensorBase.C.

 {
   InputParameters params = Material::validParams();
   params.addParam<FunctionName>(
       "elasticity_tensor_prefactor",
       "Optional function to use as a scalar prefactor on the elasticity tensor.");
   params.addParam<std::string>("base_name",
                                "Optional parameter that allows the user to define "
                                "multiple mechanics material systems on the same "
                                "block, i.e. for multiple phases");
   return params;
 }