ComputeVariableIsotropicElasticityTensor defines an elasticity tensor material for isotropic materials in which the elastic constants (Young's modulus and Poisson's ratio) vary as defined by material properties. More...

#include <ComputeVariableIsotropicElasticityTensor.h>

Inheritance diagram for ComputeVariableIsotropicElasticityTensor:

[legend]

Public Member Functions
	ComputeVariableIsotropicElasticityTensor (const InputParameters &parameters)

bool	hasGuarantee (const MaterialPropertyName &prop_name, Guarantee guarantee)

Static Public Member Functions
static InputParameters	validParams ()

Protected Member Functions
virtual void	initialSetup () override

virtual void	initQpStatefulProperties () override

virtual void	computeQpElasticityTensor () override

virtual void	computeQpProperties ()

void	issueGuarantee (const MaterialPropertyName &prop_name, Guarantee guarantee)

void	revokeGuarantee (const MaterialPropertyName &prop_name, Guarantee guarantee)

Protected Attributes
const MaterialProperty< Real > &	_youngs_modulus
	Material defining the Young's Modulus. More...

const MaterialProperty< Real > &	_poissons_ratio
	Material defining the Poisson's Ratio. More...

const unsigned int	_num_args
	number of variables the moduli depend on More...

std::vector< const MaterialProperty< Real > * >	_dyoungs_modulus
	first derivatives of the Young's Modulus with respect to the args More...

std::vector< std::vector< const MaterialProperty< Real > * > >	_d2youngs_modulus
	second derivatives of the Young's Modulus with respect to the args More...

std::vector< const MaterialProperty< Real > * >	_dpoissons_ratio
	first derivatives of the Poisson's Ratio with respect to the args More...

std::vector< std::vector< const MaterialProperty< Real > * > >	_d2poissons_ratio
	second derivatives of the Poisson's Ratio with respect to the args More...

std::vector< MaterialProperty< RankFourTensor > * >	_delasticity_tensor
	first derivatives of the elasticity tensor with respect to the args More...

std::vector< std::vector< MaterialProperty< RankFourTensor > * > >	_d2elasticity_tensor
	second derivatives of the elasticity tensor with respect to the args More...

std::vector< Real >	_isotropic_elastic_constants
	Vector of elastic constants to create the elasticity tensor (member to avoid memory churn) 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

ComputeVariableIsotropicElasticityTensor defines an elasticity tensor material for isotropic materials in which the elastic constants (Young's modulus and Poisson's ratio) vary as defined by material properties.

Definition at line 24 of file ComputeVariableIsotropicElasticityTensor.h.

Constructor & Destructor Documentation

◆ ComputeVariableIsotropicElasticityTensor()

ComputeVariableIsotropicElasticityTensor::ComputeVariableIsotropicElasticityTensor ( const InputParameters & parameters )

Definition at line 31 of file ComputeVariableIsotropicElasticityTensor.C.

   : ComputeElasticityTensorBase(parameters),
     _youngs_modulus(getMaterialProperty<Real>("youngs_modulus")),
     _poissons_ratio(getMaterialProperty<Real>("poissons_ratio")),
     _num_args(coupledComponents("args")),
     _dyoungs_modulus(_num_args),
     _d2youngs_modulus(_num_args),
     _dpoissons_ratio(_num_args),
     _d2poissons_ratio(_num_args),
     _delasticity_tensor(_num_args),
     _d2elasticity_tensor(_num_args),
     _isotropic_elastic_constants(2)
 {
   // all tensors created by this class are always isotropic
   issueGuarantee(_elasticity_tensor_name, Guarantee::ISOTROPIC);
  
   // fetch prerequisite derivatives and build elasticity tensor derivatives and cross-derivatives
   for (unsigned int i = 0; i < _num_args; ++i)
   {
     const VariableName & iname = getVar("args", i)->name();
     _dyoungs_modulus[i] = &getMaterialPropertyDerivative<Real>("youngs_modulus", iname);
     _dpoissons_ratio[i] = &getMaterialPropertyDerivative<Real>("poissons_ratio", iname);
  
     _delasticity_tensor[i] =
         &declarePropertyDerivative<RankFourTensor>(_elasticity_tensor_name, iname);
  
     _d2youngs_modulus[i].resize(_num_args);
     _d2poissons_ratio[i].resize(_num_args);
     _d2elasticity_tensor[i].resize(_num_args);
  
     for (unsigned int j = i; j < _num_args; ++j)
     {
       const VariableName & jname = getVar("args", j)->name();
       _d2youngs_modulus[i][j] =
           &getMaterialPropertyDerivative<Real>("youngs_modulus", iname, jname);
       _d2poissons_ratio[i][j] =
           &getMaterialPropertyDerivative<Real>("poissons_ratio", iname, jname);
       _d2elasticity_tensor[i][j] =
           &declarePropertyDerivative<RankFourTensor>(_elasticity_tensor_name, iname, jname);
     }
   }
 }

Member Function Documentation

◆ computeQpElasticityTensor()

void ComputeVariableIsotropicElasticityTensor::computeQpElasticityTensor ( )

overrideprotectedvirtual

Implements ComputeElasticityTensorBase.

Definition at line 104 of file ComputeVariableIsotropicElasticityTensor.C.

 {
   const Real E = _youngs_modulus[_qp];
   const Real nu = _poissons_ratio[_qp];
  
   _elasticity_tensor[_qp].fillSymmetricIsotropicEandNu(E, nu);
  
   // Define derivatives of the elasticity tensor
   for (unsigned int i = 0; i < _num_args; ++i)
   {
     if (_delasticity_tensor[i])
     {
       const Real dE = (*_dyoungs_modulus[i])[_qp];
       const Real dnu = (*_dpoissons_ratio[i])[_qp];
  
       const Real dlambda = (E * dnu + dE * nu) / ((1.0 + nu) * (1.0 - 2.0 * nu)) -
                            E * nu * dnu / ((1.0 + nu) * (1.0 + nu) * (1.0 - 2.0 * nu)) +
                            2.0 * E * nu * dnu / ((1.0 + nu) * (1.0 - 2.0 * nu) * (1.0 - 2.0 * nu));
       const Real dG = dE / (2.0 * (1.0 + nu)) - 2.0 * E * dnu / (4.0 * (1.0 + nu) * (1.0 + nu));
  
       (*_delasticity_tensor[i])[_qp].fillGeneralIsotropic(dlambda, dG, 0.0);
     }
  
     for (unsigned int j = i; j < _num_args; ++j)
       if (_d2elasticity_tensor[i][j])
       {
         const Real dEi = (*_dyoungs_modulus[i])[_qp];
         const Real dnui = (*_dpoissons_ratio[i])[_qp];
  
         const Real dEj = (*_dyoungs_modulus[j])[_qp];
         const Real dnuj = (*_dpoissons_ratio[j])[_qp];
  
         const Real d2E = (*_d2youngs_modulus[i][j])[_qp];
         const Real d2nu = (*_d2poissons_ratio[i][j])[_qp];
  
         const Real d2lambda =
             1.0 / ((1.0 + nu) * (2.0 * nu - 1.0)) *
             (-E * d2nu - nu * d2E - dEi * dnuj - dEj * dnui +
              (2.0 * E * d2nu * nu + 4.0 * dnui * dnuj * E + 2.0 * dEi * dnuj * nu +
               2.0 * dEj * dnui * nu) /
                  (2.0 * nu - 1.0) -
              8.0 * dnui * dnuj * E * nu / ((2.0 * nu - 1.0) * (2.0 * nu - 1.0)) +
              (E * d2nu * nu + 2.0 * E * dnui * dnuj + dEi * dnuj * nu + dEj * dnui * nu) /
                  (nu + 1.0) -
              4.0 * E * nu * dnui * dnuj / ((1.0 + nu) * (2.0 * nu - 1.0)) -
              2.0 * E * dnui * dnuj * nu / ((nu + 1.0) * (nu + 1.0)));
         const Real d2G = 1.0 / (nu + 1.0) *
                          (0.5 * d2E - (E * d2nu + dEi * dnuj + dEj * dnui) / (2.0 * nu + 2.0) +
                           dnui * dnuj * E / ((nu + 1.0) * (nu + 1.0)));
  
         (*_d2elasticity_tensor[i][j])[_qp].fillGeneralIsotropic(d2lambda, d2G, 0.0);
       }
   }
 }

◆ 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();
 }

◆ initialSetup()

void ComputeVariableIsotropicElasticityTensor::initialSetup ( )

overrideprotectedvirtual

Definition at line 76 of file ComputeVariableIsotropicElasticityTensor.C.

 {
   validateCoupling<Real>("youngs_modulus");
   validateCoupling<Real>("poissons_ratio");
   for (unsigned int i = 0; i < _num_args; ++i)
   {
     const VariableName & iname = getVar("args", i)->name();
  
     if (!_fe_problem.isMatPropRequested(
             derivativePropertyNameFirst(_elasticity_tensor_name, iname)))
       _delasticity_tensor[i] = nullptr;
  
     for (unsigned int j = 0; j < _num_args; ++j)
     {
       const VariableName & jname = getVar("args", j)->name();
       if (!_fe_problem.isMatPropRequested(
               derivativePropertyNameSecond(_elasticity_tensor_name, iname, jname)))
         _d2elasticity_tensor[i][j] = nullptr;
     }
   }
 }

◆ initQpStatefulProperties()

void ComputeVariableIsotropicElasticityTensor::initQpStatefulProperties ( )

overrideprotectedvirtual

Definition at line 99 of file ComputeVariableIsotropicElasticityTensor.C.

100 {

101 }

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

static

Definition at line 17 of file ComputeVariableIsotropicElasticityTensor.C.

 {
   InputParameters params = ComputeElasticityTensorBase::validParams();
   params.addClassDescription("Compute an isotropic elasticity tensor for elastic constants that "
                              "change as a function of material properties");
   params.addRequiredParam<MaterialPropertyName>("youngs_modulus",
                                                 "Name of material defining the Young's Modulus");
   params.addRequiredParam<MaterialPropertyName>("poissons_ratio",
                                                 "Name of material defining the Poisson's Ratio");
   params.addRequiredCoupledVar(
       "args", "Variable dependence for the Young's Modulus and Poisson's Ratio materials");
   return params;
 }