FiniteStrainUObasedCP Class Reference

FiniteStrainUObasedCP uses the multiplicative decomposition of deformation gradient and solves the PK2 stress residual equation at the intermediate configuration to evolve the material state. More...

#include <FiniteStrainUObasedCP.h>

Inheritance diagram for FiniteStrainUObasedCP:

Public Types
typedef DerivativeMaterialPropertyNameInterface::SymbolName	SymbolName

Public Member Functions
	FiniteStrainUObasedCP (const InputParameters &parameters)
const GenericMaterialProperty< U, is_ad > &	getDefaultMaterialProperty (const std::string &name)
const GenericMaterialProperty< U, is_ad > &	getDefaultMaterialPropertyByName (const std::string &name)
void	validateDerivativeMaterialPropertyBase (const std::string &base)
const MaterialPropertyName	derivativePropertyName (const MaterialPropertyName &base, const std::vector< SymbolName > &c) const
const MaterialPropertyName	derivativePropertyNameFirst (const MaterialPropertyName &base, const SymbolName &c1) const
const MaterialPropertyName	derivativePropertyNameSecond (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2) const
const MaterialPropertyName	derivativePropertyNameThird (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2, const SymbolName &c3) const
GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const std::vector< VariableName > &c)
GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)
GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")
GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const std::vector< VariableName > &c)
GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)
GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")
const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const std::vector< VariableName > &c)
const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)
const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")
const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, unsigned int v2, unsigned int v3=libMesh::invalid_uint)
const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, unsigned int v1, unsigned int v2=libMesh::invalid_uint, unsigned int v3=libMesh::invalid_uint)
const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const std::vector< VariableName > &c)
const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)
const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")
const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, unsigned int v2, unsigned int v3=libMesh::invalid_uint)
const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, unsigned int v1, unsigned int v2=libMesh::invalid_uint, unsigned int v3=libMesh::invalid_uint)
const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< VariableName > &c)
const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< SymbolName > &c)
const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")
const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< VariableName > &c)
const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< SymbolName > &c)
const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")
void	validateCoupling (const MaterialPropertyName &base, const std::vector< VariableName > &c, bool validate_aux=true)
void	validateCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")
void	validateCoupling (const MaterialPropertyName &base, const std::vector< VariableName > &c, bool validate_aux=true)
void	validateCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")
void	validateNonlinearCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")
void	validateNonlinearCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")
const MaterialPropertyName	propertyName (const MaterialPropertyName &base, const std::vector< SymbolName > &c) const
const MaterialPropertyName	propertyName (const MaterialPropertyName &base, const std::vector< SymbolName > &c) const
const MaterialPropertyName	propertyNameFirst (const MaterialPropertyName &base, const SymbolName &c1) const
const MaterialPropertyName	propertyNameFirst (const MaterialPropertyName &base, const SymbolName &c1) const
const MaterialPropertyName	propertyNameSecond (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2) const
const MaterialPropertyName	propertyNameSecond (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2) const
const MaterialPropertyName	propertyNameThird (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2, const SymbolName &c3) const
const MaterialPropertyName	propertyNameThird (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2, const SymbolName &c3) const

Static Public Member Functions
static InputParameters	validParams ()

Protected Member Functions
virtual void	computeQpStress ()
	updates the stress at a quadrature point.
virtual void	initQpStatefulProperties ()
	initializes the stateful properties such as stress, plastic deformation gradient, slip system resistances, etc.
virtual void	calcResidJacob ()
	calls the residual and jacobian functions used in the stress update algorithm.
virtual void	updateSlipSystemResistanceAndStateVariable ()
	updates the slip system resistances and state variables.
virtual void	preSolveQp ()
	set variables for stress and internal variable solve.
virtual void	solveQp ()
	solve stress and internal variables.
virtual void	postSolveQp ()
	update stress and internal variable after solve.
virtual void	preSolveStatevar ()
	set variables for internal variable solve.
virtual void	solveStatevar ()
	solve internal variables.
virtual void	postSolveStatevar ()
	update internal variable after solve.
virtual void	preSolveStress ()
	set variables for stress solve.
virtual void	solveStress ()
	solves for stress, updates plastic deformation gradient.
virtual void	postSolveStress ()
	update stress and plastic deformation gradient after solve.
virtual void	calcResidual ()
	calculate stress residual.
virtual void	calcJacobian ()
	calculate jacobian.
virtual void	getSlipRates ()
	updates the slip rates.
virtual void	calcTangentModuli ()
	calculate the tangent moduli for preconditioner.
virtual void	elasticTangentModuli ()
	calculate the elastic tangent moduli for preconditioner.
virtual void	elastoPlasticTangentModuli ()
	calculate the exact tangent moduli for preconditioner.
bool	lineSearchUpdate (const Real rnorm_prev, const RankTwoTensor)
	performs the line search update
virtual bool	isStateVariablesConverged ()
	evaluates convergence of state variables.
virtual void	computeQpProperties () override

Protected Attributes
std::vector< const CrystalPlasticitySlipRate * >	_uo_slip_rates
	User objects that define the slip rate.
std::vector< const CrystalPlasticitySlipResistance * >	_uo_slip_resistances
	User objects that define the slip resistance.
std::vector< const CrystalPlasticityStateVariable * >	_uo_state_vars
	User objects that define the state variable.
std::vector< const CrystalPlasticityStateVarRateComponent * >	_uo_state_var_evol_rate_comps
	User objects that define the state variable evolution rate component.
std::vector< MaterialProperty< std::vector< Real > > * >	_mat_prop_slip_rates
	Slip rates material property.
std::vector< MaterialProperty< std::vector< Real > > * >	_mat_prop_slip_resistances
	Slip resistance material property.
std::vector< MaterialProperty< std::vector< Real > > * >	_mat_prop_state_vars
	State variable material property.
std::vector< const MaterialProperty< std::vector< Real > > * >	_mat_prop_state_vars_old
	Old state variable material property.
std::vector< MaterialProperty< std::vector< Real > > * >	_mat_prop_state_var_evol_rate_comps
	State variable evolution rate component material property.
unsigned int	_num_uo_slip_rates
	Number of slip rate user objects.
unsigned int	_num_uo_slip_resistances
	Number of slip resistance user objects.
unsigned int	_num_uo_state_vars
	Number of state variable user objects.
unsigned int	_num_uo_state_var_evol_rate_comps
	Number of state variable evolution rate component user objects.
std::vector< std::vector< Real > >	_state_vars_old
	Local state variable.
std::vector< std::vector< Real > >	_state_vars_old_stored
	Local stored state variable (for sub-stepping)
std::vector< std::vector< Real > >	_state_vars_prev
	Local old state variable.
Real	_rtol
	Stress residual equation relative tolerance.
Real	_abs_tol
	Stress residual equation absolute tolerance.
Real	_stol
	Internal variable update equation tolerance.
Real	_zero_tol
	Residual tolerance when variable value is zero. Default 1e-12.
RankTwoTensor	_resid
	Residual tensor.
RankFourTensor	_jac
	Jacobian tensor.
unsigned int	_maxiter
	Maximum number of iterations for stress update.
unsigned int	_maxiterg
	Maximum number of iterations for internal variable update.
MooseEnum	_tan_mod_type
	Type of tangent moduli calculation.
unsigned int	_max_substep_iter
	Maximum number of substep iterations.
bool	_use_line_search
	Flag to activate line serach.
Real	_min_lsrch_step
	Minimum line search step size.
Real	_lsrch_tol
	Line search bisection method tolerance.
unsigned int	_lsrch_max_iter
	Line search bisection method maximum iteration number.
MooseEnum	_lsrch_method
	Line search method.
MaterialProperty< RankTwoTensor > &	_fp
const MaterialProperty< RankTwoTensor > &	_fp_old
MaterialProperty< RankTwoTensor > &	_pk2
const MaterialProperty< RankTwoTensor > &	_pk2_old
MaterialProperty< RankTwoTensor > &	_lag_e
MaterialProperty< RankTwoTensor > &	_update_rot
const MaterialProperty< RankTwoTensor > &	_update_rot_old
const std::string	_elasticity_tensor_name
	Name of the elasticity tensor material property.
const MaterialProperty< RankFourTensor > &	_elasticity_tensor
	Elasticity tensor material property.
const MaterialProperty< RankTwoTensor > &	_deformation_gradient
const MaterialProperty< RankTwoTensor > &	_deformation_gradient_old
const MaterialProperty< RankTwoTensor > &	_crysrot
	Crystal rotation.
RankTwoTensor	_dfgrd_tmp
RankTwoTensor	_fe
RankTwoTensor	_fp_old_inv
RankTwoTensor	_fp_inv
DenseVector< Real >	_tau
std::vector< MaterialProperty< std::vector< RankTwoTensor > > * >	_flow_direction
bool	_err_tol
	Flag to check whether convergence is achieved.
RankTwoTensor	_delta_dfgrd
	Used for substepping; Uniformly divides the increment in deformation gradient.
RankTwoTensor	_dfgrd_tmp_old
Real	_dfgrd_scale_factor
	Scales the substepping increment to obtain deformation gradient at a substep iteration.
Real	_substep_dt
	Current substep size.
const std::string	_base_name
	Base name prepended to all material property names to allow for multi-material systems.
const MaterialProperty< RankTwoTensor > &	_mechanical_strain
	Mechanical strain material property.
MaterialProperty< RankTwoTensor > &	_stress
	Stress material property.
MaterialProperty< RankTwoTensor > &	_elastic_strain
	Elastic strain material property.
const MaterialProperty< RankTwoTensor > &	_extra_stress
	Extra stress tensor.
std::vector< const Function * >	_initial_stress_fcn
	initial stress components
MaterialProperty< RankFourTensor > &	_Jacobian_mult
	derivative of stress w.r.t. strain (_dstress_dstrain)

Private Member Functions
bool	haveMaterialProperty (const std::string &prop_name)
std::vector< VariableName >	buildVariableVector (const VariableName &c1, const VariableName &c2, const VariableName &c3)
void	validateCouplingHelper (const MaterialPropertyName &base, const std::vector< VariableName > &c, const System &system, std::vector< VariableName > &missing)
bool	isNotObjectVariable (const VariableName &name)

Private Attributes
FEProblemBase &	_dmi_fe_problem

Detailed Description

FiniteStrainUObasedCP uses the multiplicative decomposition of deformation gradient and solves the PK2 stress residual equation at the intermediate configuration to evolve the material state.

The internal variables are updated using an interative predictor-corrector algorithm. Backward Euler integration rule is used for the rate equations.

Involves 4 different types of user objects that calculates: State variables - update state variable (derive from CrystalPlasticityStateVariable) State variable evolution compoment - individual component of state variable incremental rate (derive from CrystalPlasticityStateVariableEvolutionRateComponent) Slip resistance - calcuate slip resistances (derive from CrystalPlasticitySlipResistances) Slip rates - calcuate flow direction and slip rates (derive from CrystalPlasticitySlipRates)

Definition at line 33 of file FiniteStrainUObasedCP.h.

Constructor & Destructor Documentation

FiniteStrainUObasedCP()

FiniteStrainUObasedCP::FiniteStrainUObasedCP

(

const InputParameters &

parameters

)

Definition at line 64 of file FiniteStrainUObasedCP.C.

  : ComputeStressBase(parameters),
    _num_uo_slip_rates(parameters.get<std::vector<UserObjectName>>("uo_slip_rates").size()),
    _num_uo_slip_resistances(
        parameters.get<std::vector<UserObjectName>>("uo_slip_resistances").size()),
    _num_uo_state_vars(parameters.get<std::vector<UserObjectName>>("uo_state_vars").size()),
    _num_uo_state_var_evol_rate_comps(
        parameters.get<std::vector<UserObjectName>>("uo_state_var_evol_rate_comps").size()),
    _rtol(getParam<Real>("rtol")),
    _abs_tol(getParam<Real>("abs_tol")),
    _stol(getParam<Real>("stol")),
    _zero_tol(getParam<Real>("zero_tol")),
    _maxiter(getParam<unsigned int>("maxiter")),
    _maxiterg(getParam<unsigned int>("maxiter_state_variable")),
    _tan_mod_type(getParam<MooseEnum>("tan_mod_type")),
    _max_substep_iter(getParam<unsigned int>("maximum_substep_iteration")),
    _use_line_search(getParam<bool>("use_line_search")),
    _min_lsrch_step(getParam<Real>("min_line_search_step_size")),
    _lsrch_tol(getParam<Real>("line_search_tol")),
    _lsrch_max_iter(getParam<unsigned int>("line_search_maxiter")),
    _lsrch_method(getParam<MooseEnum>("line_search_method")),
    _fp(declareProperty<RankTwoTensor>("fp")), // Plastic deformation gradient
    _fp_old(getMaterialPropertyOld<RankTwoTensor>(
        "fp")), // Plastic deformation gradient of previous increment
    _pk2(declareProperty<RankTwoTensor>("pk2")), // 2nd Piola-Kirchoff Stress
    _pk2_old(getMaterialPropertyOld<RankTwoTensor>(
        "pk2")), // 2nd Piola Kirchoff Stress of previous increment
    _lag_e(declareProperty<RankTwoTensor>("lage")), // Lagrangian strain
    _update_rot(declareProperty<RankTwoTensor>(
        "update_rot")), // Rotation tensor considering material rotation and crystal orientation
    _update_rot_old(getMaterialPropertyOld<RankTwoTensor>("update_rot")),
    _elasticity_tensor_name(_base_name + "elasticity_tensor"),
    _elasticity_tensor(getMaterialPropertyByName<RankFourTensor>(_elasticity_tensor_name)),
    _deformation_gradient(getMaterialProperty<RankTwoTensor>("deformation_gradient")),
    _deformation_gradient_old(getMaterialPropertyOld<RankTwoTensor>("deformation_gradient")),
    _crysrot(getMaterialProperty<RankTwoTensor>("crysrot"))
{
  _err_tol = false;
 
  _delta_dfgrd.zero();
 
  // resize the material properties for each userobject
  _mat_prop_slip_rates.resize(_num_uo_slip_rates);
  _mat_prop_slip_resistances.resize(_num_uo_slip_resistances);
  _mat_prop_state_vars.resize(_num_uo_state_vars);
  _mat_prop_state_vars_old.resize(_num_uo_state_vars);
  _mat_prop_state_var_evol_rate_comps.resize(_num_uo_state_var_evol_rate_comps);
 
  // resize the flow direction
  _flow_direction.resize(_num_uo_slip_rates);
 
  // resize local state variables
  _state_vars_old.resize(_num_uo_state_vars);
  _state_vars_old_stored.resize(_num_uo_state_vars);
  _state_vars_prev.resize(_num_uo_state_vars);
 
  // resize user objects
  _uo_slip_rates.resize(_num_uo_slip_rates);
  _uo_slip_resistances.resize(_num_uo_slip_resistances);
  _uo_state_vars.resize(_num_uo_state_vars);
  _uo_state_var_evol_rate_comps.resize(_num_uo_state_var_evol_rate_comps);
 
  // assign the user objects
  for (unsigned int i = 0; i < _num_uo_slip_rates; ++i)
  {
    _uo_slip_rates[i] = &getUserObjectByName<CrystalPlasticitySlipRate>(
        parameters.get<std::vector<UserObjectName>>("uo_slip_rates")[i]);
    _mat_prop_slip_rates[i] = &declareProperty<std::vector<Real>>(
        parameters.get<std::vector<UserObjectName>>("uo_slip_rates")[i]);
    _flow_direction[i] = &declareProperty<std::vector<RankTwoTensor>>(
        parameters.get<std::vector<UserObjectName>>("uo_slip_rates")[i] + "_flow_direction");
  }
 
  for (unsigned int i = 0; i < _num_uo_slip_resistances; ++i)
  {
    _uo_slip_resistances[i] = &getUserObjectByName<CrystalPlasticitySlipResistance>(
        parameters.get<std::vector<UserObjectName>>("uo_slip_resistances")[i]);
    _mat_prop_slip_resistances[i] = &declareProperty<std::vector<Real>>(
        parameters.get<std::vector<UserObjectName>>("uo_slip_resistances")[i]);
  }
 
  for (unsigned int i = 0; i < _num_uo_state_vars; ++i)
  {
    _uo_state_vars[i] = &getUserObjectByName<CrystalPlasticityStateVariable>(
        parameters.get<std::vector<UserObjectName>>("uo_state_vars")[i]);
    _mat_prop_state_vars[i] = &declareProperty<std::vector<Real>>(
        parameters.get<std::vector<UserObjectName>>("uo_state_vars")[i]);
    _mat_prop_state_vars_old[i] = &getMaterialPropertyOld<std::vector<Real>>(
        parameters.get<std::vector<UserObjectName>>("uo_state_vars")[i]);
  }
 
  for (unsigned int i = 0; i < _num_uo_state_var_evol_rate_comps; ++i)
  {
    _uo_state_var_evol_rate_comps[i] = &getUserObjectByName<CrystalPlasticityStateVarRateComponent>(
        parameters.get<std::vector<UserObjectName>>("uo_state_var_evol_rate_comps")[i]);
    _mat_prop_state_var_evol_rate_comps[i] = &declareProperty<std::vector<Real>>(
        parameters.get<std::vector<UserObjectName>>("uo_state_var_evol_rate_comps")[i]);
  }
 
  _substep_dt = 0.0;
}

Member Function Documentation

calcJacobian()

void FiniteStrainUObasedCP::calcJacobian ( )

protectedvirtual

calculate jacobian.

Definition at line 541 of file FiniteStrainUObasedCP.C.

{
  RankFourTensor dfedfpinv, deedfe, dfpinvdpk2;
 
  for (const auto i : make_range(Moose::dim))
    for (const auto j : make_range(Moose::dim))
      for (const auto k : make_range(Moose::dim))
        dfedfpinv(i, j, k, j) = _dfgrd_tmp(i, k);
 
  for (const auto i : make_range(Moose::dim))
    for (const auto j : make_range(Moose::dim))
      for (const auto k : make_range(Moose::dim))
      {
        deedfe(i, j, k, i) = deedfe(i, j, k, i) + _fe(k, j) * 0.5;
        deedfe(i, j, k, j) = deedfe(i, j, k, j) + _fe(k, i) * 0.5;
      }
 
  for (unsigned int i = 0; i < _num_uo_slip_rates; ++i)
  {
    unsigned int nss = _uo_slip_rates[i]->variableSize();
    std::vector<RankTwoTensor> dtaudpk2(nss), dfpinvdslip(nss);
    std::vector<Real> dslipdtau;
    dslipdtau.resize(nss);
    _uo_slip_rates[i]->calcSlipRateDerivative(_qp, _substep_dt, dslipdtau);
    for (unsigned int j = 0; j < nss; j++)
    {
      dtaudpk2[j] = (*_flow_direction[i])[_qp][j];
      dfpinvdslip[j] = -_fp_old_inv * (*_flow_direction[i])[_qp][j];
      dfpinvdpk2 += (dfpinvdslip[j] * dslipdtau[j] * _substep_dt).outerProduct(dtaudpk2[j]);
    }
  }
  _jac =
      RankFourTensor::IdentityFour() - (_elasticity_tensor[_qp] * deedfe * dfedfpinv * dfpinvdpk2);
}

Referenced by calcResidJacob().

calcResidJacob()

void FiniteStrainUObasedCP::calcResidJacob ( )

protectedvirtual

calls the residual and jacobian functions used in the stress update algorithm.

Definition at line 492 of file FiniteStrainUObasedCP.C.

{
  calcResidual();
  if (_err_tol)
    return;
  calcJacobian();
}

Referenced by solveStress().

calcResidual()

void FiniteStrainUObasedCP::calcResidual ( )

protectedvirtual

calculate stress residual.

Definition at line 514 of file FiniteStrainUObasedCP.C.

{
  RankTwoTensor iden(RankTwoTensor::initIdentity), ce, ee, ce_pk2, eqv_slip_incr, pk2_new;
 
  getSlipRates();
  if (_err_tol)
    return;
 
  for (unsigned int i = 0; i < _num_uo_slip_rates; ++i)
    for (unsigned int j = 0; j < _uo_slip_rates[i]->variableSize(); ++j)
      eqv_slip_incr +=
          (*_flow_direction[i])[_qp][j] * (*_mat_prop_slip_rates[i])[_qp][j] * _substep_dt;
 
  eqv_slip_incr = iden - eqv_slip_incr;
  _fp_inv = _fp_old_inv * eqv_slip_incr;
  _fe = _dfgrd_tmp * _fp_inv;
 
  ce = _fe.transpose() * _fe;
  ee = ce - iden;
  ee *= 0.5;
 
  pk2_new = _elasticity_tensor[_qp] * ee;
 
  _resid = _pk2[_qp] - pk2_new;
}

Referenced by calcResidJacob(), and lineSearchUpdate().

calcTangentModuli()

void FiniteStrainUObasedCP::calcTangentModuli ( )

protectedvirtual

calculate the tangent moduli for preconditioner.

Default is the elastic stiffness matrix. Exact jacobian is currently implemented. tan_mod_type can be modified to exact in .i file to turn it on.

Definition at line 577 of file FiniteStrainUObasedCP.C.

{
  switch (_tan_mod_type)
  {
    case 0:
      elastoPlasticTangentModuli();
      break;
    default:
      elasticTangentModuli();
  }
}

Referenced by postSolveQp().

computeQpProperties()

void ComputeGeneralStressBase::computeQpProperties ( )

overrideprotectedvirtualinherited

Definition at line 44 of file ComputeGeneralStressBase.C.

{
  computeQpStress();
 
  // Add in extra stress
  _stress[_qp] += _extra_stress[_qp];
}

computeQpStress()

void FiniteStrainUObasedCP::computeQpStress ( )

protectedvirtual

updates the stress at a quadrature point.

Solves stress residual equation using NR.

Updates slip system resistances iteratively.

Implements ComputeGeneralStressBase.

Definition at line 207 of file FiniteStrainUObasedCP.C.

{
  // Userobject based crystal plasticity does not support face/boundary material property
  // calculation.
  if (isBoundaryMaterial())
    return;
  // Depth of substepping; Limited to maximum substep iteration
  unsigned int substep_iter = 1;
  // Calculated from substep_iter as 2^substep_iter
  unsigned int num_substep = 1;
 
  _dfgrd_tmp_old = _deformation_gradient_old[_qp];
  if (_dfgrd_tmp_old.det() == 0)
    _dfgrd_tmp_old.addIa(1.0);
 
  _delta_dfgrd = _deformation_gradient[_qp] - _dfgrd_tmp_old;
 
  // Saves the old stateful properties that are modified during sub stepping
  for (unsigned int i = 0; i < _num_uo_state_vars; ++i)
    _state_vars_old[i] = (*_mat_prop_state_vars_old[i])[_qp];
 
  for (unsigned int i = 0; i < _num_uo_slip_rates; ++i)
    _uo_slip_rates[i]->calcFlowDirection(_qp, (*_flow_direction[i])[_qp]);
 
  do
  {
    _err_tol = false;
 
    preSolveQp();
 
    _substep_dt = _dt / num_substep;
 
    for (unsigned int istep = 0; istep < num_substep; ++istep)
    {
      _dfgrd_tmp = (static_cast<Real>(istep) + 1) / num_substep * _delta_dfgrd + _dfgrd_tmp_old;
 
      solveQp();
 
      if (_err_tol)
      {
        substep_iter++;
        num_substep *= 2;
        break;
      }
    }
 
    if (substep_iter > _max_substep_iter && _err_tol)
      throw MooseException("FiniteStrainUObasedCP: Constitutive failure.");
  } while (_err_tol);
 
  postSolveQp();
}

elasticTangentModuli()

void FiniteStrainUObasedCP::elasticTangentModuli ( )

protectedvirtual

calculate the elastic tangent moduli for preconditioner.

Definition at line 634 of file FiniteStrainUObasedCP.C.

{
  // update jacobian_mult
  _Jacobian_mult[_qp] = _elasticity_tensor[_qp];
}

Referenced by calcTangentModuli().

elastoPlasticTangentModuli()

void FiniteStrainUObasedCP::elastoPlasticTangentModuli ( )

protectedvirtual

calculate the exact tangent moduli for preconditioner.

Definition at line 590 of file FiniteStrainUObasedCP.C.

{
  RankFourTensor tan_mod;
  RankTwoTensor pk2fet, fepk2;
  RankFourTensor deedfe, dsigdpk2dfe, dfedf;
 
  // Fill in the matrix stiffness material property
  for (const auto i : make_range(Moose::dim))
    for (const auto j : make_range(Moose::dim))
      for (const auto k : make_range(Moose::dim))
      {
        deedfe(i, j, k, i) = deedfe(i, j, k, i) + _fe(k, j) * 0.5;
        deedfe(i, j, k, j) = deedfe(i, j, k, j) + _fe(k, i) * 0.5;
      }
 
  usingTensorIndices(i_, j_, k_, l_);
  dsigdpk2dfe = _fe.times<i_, k_, j_, l_>(_fe) * _elasticity_tensor[_qp] * deedfe;
 
  pk2fet = _pk2[_qp] * _fe.transpose();
  fepk2 = _fe * _pk2[_qp];
 
  for (const auto i : make_range(Moose::dim))
    for (const auto j : make_range(Moose::dim))
      for (const auto l : make_range(Moose::dim))
      {
        tan_mod(i, j, i, l) += pk2fet(l, j);
        tan_mod(i, j, j, l) += fepk2(i, l);
      }
 
  tan_mod += dsigdpk2dfe;
 
  Real je = _fe.det();
  if (je > 0.0)
    tan_mod /= je;
 
  for (const auto i : make_range(Moose::dim))
    for (const auto j : make_range(Moose::dim))
      for (const auto l : make_range(Moose::dim))
        dfedf(i, j, i, l) = _fp_inv(l, j);
 
  _Jacobian_mult[_qp] = tan_mod * dfedf;
}

Referenced by calcTangentModuli().

getSlipRates()

void FiniteStrainUObasedCP::getSlipRates ( )

protectedvirtual

updates the slip rates.

Definition at line 501 of file FiniteStrainUObasedCP.C.

{
  for (unsigned int i = 0; i < _num_uo_slip_rates; ++i)
  {
    if (!_uo_slip_rates[i]->calcSlipRate(_qp, _substep_dt, (*_mat_prop_slip_rates[i])[_qp]))
    {
      _err_tol = true;
      return;
    }
  }
}

Referenced by calcResidual().

initQpStatefulProperties()

void FiniteStrainUObasedCP::initQpStatefulProperties ( )

protectedvirtual

initializes the stateful properties such as stress, plastic deformation gradient, slip system resistances, etc.

Reimplemented from ComputeGeneralStressBase.

Definition at line 167 of file FiniteStrainUObasedCP.C.

{
  for (unsigned int i = 0; i < _num_uo_slip_rates; ++i)
  {
    (*_mat_prop_slip_rates[i])[_qp].resize(_uo_slip_rates[i]->variableSize());
    (*_flow_direction[i])[_qp].resize(_uo_slip_rates[i]->variableSize());
  }
 
  for (unsigned int i = 0; i < _num_uo_slip_resistances; ++i)
    (*_mat_prop_slip_resistances[i])[_qp].resize(_uo_slip_resistances[i]->variableSize());
 
  for (unsigned int i = 0; i < _num_uo_state_vars; ++i)
  {
    (*_mat_prop_state_vars[i])[_qp].resize(_uo_state_vars[i]->variableSize());
    _state_vars_old[i].resize(_uo_state_vars[i]->variableSize());
    _state_vars_old_stored[i].resize(_uo_state_vars[i]->variableSize());
    _state_vars_prev[i].resize(_uo_state_vars[i]->variableSize());
  }
 
  for (unsigned int i = 0; i < _num_uo_state_var_evol_rate_comps; ++i)
    (*_mat_prop_state_var_evol_rate_comps[i])[_qp].resize(
        _uo_state_var_evol_rate_comps[i]->variableSize());
 
  _stress[_qp].zero();
  _pk2[_qp].zero();
  _lag_e[_qp].zero();
 
  _fp[_qp].setToIdentity();
  _update_rot[_qp].setToIdentity();
 
  for (unsigned int i = 0; i < _num_uo_state_vars; ++i)
    // Initializes slip system related properties
    _uo_state_vars[i]->initSlipSysProps((*_mat_prop_state_vars[i])[_qp], _q_point[_qp]);
}

isStateVariablesConverged()

bool FiniteStrainUObasedCP::isStateVariablesConverged ( )

protectedvirtual

evaluates convergence of state variables.

Definition at line 349 of file FiniteStrainUObasedCP.C.

{
  Real diff;
 
  for (unsigned int i = 0; i < _num_uo_state_vars; ++i)
  {
    unsigned int n = (*_mat_prop_state_vars[i])[_qp].size();
    for (unsigned j = 0; j < n; j++)
    {
      diff = std::abs((*_mat_prop_state_vars[i])[_qp][j] - _state_vars_prev[i][j]);
 
      if (std::abs(_state_vars_old_stored[i][j]) < _zero_tol && diff > _zero_tol)
        return true;
      if (std::abs(_state_vars_old_stored[i][j]) > _zero_tol &&
          diff > _stol * std::abs(_state_vars_old_stored[i][j]))
        return true;
    }
  }
  return false;
}

Referenced by solveStatevar().

lineSearchUpdate()

bool FiniteStrainUObasedCP::lineSearchUpdate ( const Real  rnorm_prev, const RankTwoTensor  dpk2  )

protected

performs the line search update

Definition at line 641 of file FiniteStrainUObasedCP.C.

{
  switch (_lsrch_method)
  {
    case 0: // CUT_HALF
    {
      Real rnorm;
      Real step = 1.0;
 
      do
      {
        _pk2[_qp] = _pk2[_qp] - step * dpk2;
        step /= 2.0;
        _pk2[_qp] = _pk2[_qp] + step * dpk2;
 
        calcResidual();
        rnorm = _resid.L2norm();
      } while (rnorm > rnorm_prev && step > _min_lsrch_step);
 
      // has norm improved or is the step still above minumum search step size?
      return (rnorm <= rnorm_prev || step > _min_lsrch_step);
    }
 
    case 1: // BISECTION
    {
      unsigned int count = 0;
      Real step_a = 0.0;
      Real step_b = 1.0;
      Real step = 1.0;
      Real s_m = 1000.0;
      Real rnorm = 1000.0;
 
      calcResidual();
      Real s_b = _resid.doubleContraction(dpk2);
      Real rnorm1 = _resid.L2norm();
      _pk2[_qp] = _pk2[_qp] - dpk2;
      calcResidual();
      Real s_a = _resid.doubleContraction(dpk2);
      Real rnorm0 = _resid.L2norm();
      _pk2[_qp] = _pk2[_qp] + dpk2;
 
      if ((rnorm1 / rnorm0) < _lsrch_tol || s_a * s_b > 0)
      {
        calcResidual();
        return true;
      }
 
      while ((rnorm / rnorm0) > _lsrch_tol && count < _lsrch_max_iter)
      {
        _pk2[_qp] = _pk2[_qp] - step * dpk2;
        step = 0.5 * (step_b + step_a);
        _pk2[_qp] = _pk2[_qp] + step * dpk2;
        calcResidual();
        s_m = _resid.doubleContraction(dpk2);
        rnorm = _resid.L2norm();
 
        if (s_m * s_a < 0.0)
        {
          step_b = step;
          s_b = s_m;
        }
        if (s_m * s_b < 0.0)
        {
          step_a = step;
          s_a = s_m;
        }
        count++;
      }
 
      // below tolerance and max iterations?
      return ((rnorm / rnorm0) < _lsrch_tol && count < _lsrch_max_iter);
    }
 
    default:
      mooseError("Line search method is not provided.");
  }
}

Referenced by solveStress().

postSolveQp()

void FiniteStrainUObasedCP::postSolveQp ( )

protectedvirtual

update stress and internal variable after solve.

Definition at line 282 of file FiniteStrainUObasedCP.C.

{
  _stress[_qp] = _fe * _pk2[_qp] * _fe.transpose() / _fe.det();
 
  // Calculate jacobian for preconditioner
  calcTangentModuli();
 
  RankTwoTensor iden(RankTwoTensor::initIdentity);
 
  _lag_e[_qp] = _deformation_gradient[_qp].transpose() * _deformation_gradient[_qp] - iden;
  _lag_e[_qp] = _lag_e[_qp] * 0.5;
 
  RankTwoTensor rot;
  // Calculate material rotation
  _deformation_gradient[_qp].getRUDecompositionRotation(rot);
  _update_rot[_qp] = rot * _crysrot[_qp];
}

Referenced by computeQpStress().

postSolveStatevar()

void FiniteStrainUObasedCP::postSolveStatevar ( )

protectedvirtual

update internal variable after solve.

Definition at line 371 of file FiniteStrainUObasedCP.C.

{
  for (unsigned int i = 0; i < _num_uo_state_vars; ++i)
    _state_vars_old_stored[i] = (*_mat_prop_state_vars[i])[_qp];
 
  _fp_old_inv = _fp_inv;
}

Referenced by solveQp().

postSolveStress()

void FiniteStrainUObasedCP::postSolveStress ( )

protectedvirtual

update stress and plastic deformation gradient after solve.

Definition at line 464 of file FiniteStrainUObasedCP.C.

{
  _fp[_qp] = _fp_inv.inverse();
}

Referenced by solveStatevar().

preSolveQp()

void FiniteStrainUObasedCP::preSolveQp ( )

protectedvirtual

set variables for stress and internal variable solve.

Definition at line 261 of file FiniteStrainUObasedCP.C.

{
  for (unsigned int i = 0; i < _num_uo_state_vars; ++i)
    (*_mat_prop_state_vars[i])[_qp] = _state_vars_old_stored[i] = _state_vars_old[i];
 
  _pk2[_qp] = _pk2_old[_qp];
  _fp_old_inv = _fp_old[_qp].inverse();
}

Referenced by computeQpStress().

preSolveStatevar()

void FiniteStrainUObasedCP::preSolveStatevar ( )

protectedvirtual

set variables for internal variable solve.

Definition at line 301 of file FiniteStrainUObasedCP.C.

{
  for (unsigned int i = 0; i < _num_uo_state_vars; ++i)
    (*_mat_prop_state_vars[i])[_qp] = _state_vars_old_stored[i];
 
  for (unsigned int i = 0; i < _num_uo_slip_resistances; ++i)
    _uo_slip_resistances[i]->calcSlipResistance(_qp, (*_mat_prop_slip_resistances[i])[_qp]);
 
  _fp_inv = _fp_old_inv;
}

Referenced by solveQp().

preSolveStress()

void FiniteStrainUObasedCP::preSolveStress ( )

protectedvirtual

set variables for stress solve.

Definition at line 380 of file FiniteStrainUObasedCP.C.

{
}

Referenced by solveStatevar().

solveQp()

void FiniteStrainUObasedCP::solveQp ( )

protectedvirtual

solve stress and internal variables.

Definition at line 271 of file FiniteStrainUObasedCP.C.

{
  preSolveStatevar();
  solveStatevar();
  if (_err_tol)
    return;
 
  postSolveStatevar();
}

Referenced by computeQpStress().

solveStatevar()

void FiniteStrainUObasedCP::solveStatevar ( )

protectedvirtual

solve internal variables.

Definition at line 313 of file FiniteStrainUObasedCP.C.

{
  unsigned int iterg;
  bool iter_flag = true;
 
  iterg = 0;
  // Check for slip system resistance update tolerance
  while (iter_flag && iterg < _maxiterg)
  {
    preSolveStress();
    solveStress();
    if (_err_tol)
      return;
 
    postSolveStress();
 
    // Update slip system resistance and state variable
    updateSlipSystemResistanceAndStateVariable();
 
    if (_err_tol)
      return;
 
    iter_flag = isStateVariablesConverged();
    iterg++;
  }
 
  if (iterg == _maxiterg)
  {
#ifdef DEBUG
    mooseWarning("FiniteStrainUObasedCP: Hardness Integration error\n");
#endif
    _err_tol = true;
  }
}

Referenced by solveQp().

solveStress()

void FiniteStrainUObasedCP::solveStress ( )

protectedvirtual

solves for stress, updates plastic deformation gradient.

Definition at line 385 of file FiniteStrainUObasedCP.C.

{
  unsigned int iter = 0;
  RankTwoTensor dpk2;
  Real rnorm, rnorm0, rnorm_prev;
 
  // Calculate stress residual
  calcResidJacob();
  if (_err_tol)
  {
#ifdef DEBUG
    mooseWarning("FiniteStrainUObasedCP: Slip increment exceeds tolerance - Element number ",
                 _current_elem->id(),
                 " Gauss point = ",
                 _qp);
#endif
    return;
  }
 
  rnorm = _resid.L2norm();
  rnorm0 = rnorm;
 
  // Check for stress residual tolerance
  while (rnorm > _rtol * rnorm0 && rnorm0 > _abs_tol && iter < _maxiter)
  {
    // Calculate stress increment
    dpk2 = -_jac.invSymm() * _resid;
    _pk2[_qp] = _pk2[_qp] + dpk2;
 
    calcResidJacob();
 
    if (_err_tol)
    {
#ifdef DEBUG
      mooseWarning("FiniteStrainUObasedCP: Slip increment exceeds tolerance - Element number ",
                   _current_elem->id(),
                   " Gauss point = ",
                   _qp);
#endif
      return;
    }
 
    rnorm_prev = rnorm;
    rnorm = _resid.L2norm();
 
    if (_use_line_search && rnorm > rnorm_prev && !lineSearchUpdate(rnorm_prev, dpk2))
    {
#ifdef DEBUG
      mooseWarning("FiniteStrainUObasedCP: Failed with line search");
#endif
      _err_tol = true;
      return;
    }
 
    if (_use_line_search)
      rnorm = _resid.L2norm();
 
    iter++;
  }
 
  if (iter >= _maxiter)
  {
#ifdef DEBUG
    mooseWarning("FiniteStrainUObasedCP: Stress Integration error rmax = ",
                 rnorm,
                 " and the tolerance is ",
                 _rtol * rnorm0,
                 " when the rnorm0 value is ",
                 rnorm0,
                 " for element ",
                 _current_elem->id(),
                 " and qp ",
                 _qp);
#endif
    _err_tol = true;
  }
}

Referenced by solveStatevar().

updateSlipSystemResistanceAndStateVariable()

void FiniteStrainUObasedCP::updateSlipSystemResistanceAndStateVariable ( )

protectedvirtual

updates the slip system resistances and state variables.

override to modify slip system resistance and state variable evolution.

Definition at line 470 of file FiniteStrainUObasedCP.C.

{
  for (unsigned int i = 0; i < _num_uo_state_vars; ++i)
    _state_vars_prev[i] = (*_mat_prop_state_vars[i])[_qp];
 
  for (unsigned int i = 0; i < _num_uo_state_var_evol_rate_comps; ++i)
    _uo_state_var_evol_rate_comps[i]->calcStateVariableEvolutionRateComponent(
        _qp, (*_mat_prop_state_var_evol_rate_comps[i])[_qp]);
 
  for (unsigned int i = 0; i < _num_uo_state_vars; ++i)
  {
    if (!_uo_state_vars[i]->updateStateVariable(
            _qp, _substep_dt, (*_mat_prop_state_vars[i])[_qp], _state_vars_old_stored[i]))
      _err_tol = true;
  }
 
  for (unsigned int i = 0; i < _num_uo_slip_resistances; ++i)
    _uo_slip_resistances[i]->calcSlipResistance(_qp, (*_mat_prop_slip_resistances[i])[_qp]);
}

Referenced by solveStatevar().

validParams()

InputParameters FiniteStrainUObasedCP::validParams ( )

static

Definition at line 21 of file FiniteStrainUObasedCP.C.

{
  InputParameters params = ComputeStressBase::validParams();
  params.addClassDescription("UserObject based Crystal Plasticity system.");
  params.addParam<Real>("rtol", 1e-6, "Constitutive stress residue relative tolerance");
  params.addParam<Real>("abs_tol", 1e-6, "Constitutive stress residue absolute tolerance");
  params.addParam<Real>(
      "stol", 1e-2, "Constitutive slip system resistance relative residual tolerance");
  params.addParam<Real>(
      "zero_tol", 1e-12, "Tolerance for residual check when variable value is zero");
  params.addParam<unsigned int>("maxiter", 100, "Maximum number of iterations for stress update");
  params.addParam<unsigned int>(
      "maxiter_state_variable", 100, "Maximum number of iterations for state variable update");
  MooseEnum tan_mod_options("exact none", "none"); // Type of read
  params.addParam<MooseEnum>("tan_mod_type",
                             tan_mod_options,
                             "Type of tangent moduli for preconditioner: default elastic");
  params.addParam<unsigned int>(
      "maximum_substep_iteration", 1, "Maximum number of substep iteration");
  params.addParam<bool>("use_line_search", false, "Use line search in constitutive update");
  params.addParam<Real>("min_line_search_step_size", 0.01, "Minimum line search step size");
  params.addParam<Real>("line_search_tol", 0.5, "Line search bisection method tolerance");
  params.addParam<unsigned int>(
      "line_search_maxiter", 20, "Line search bisection method maximum number of iteration");
  MooseEnum line_search_method("CUT_HALF BISECTION", "CUT_HALF");
  params.addParam<MooseEnum>(
      "line_search_method", line_search_method, "The method used in line search");
  params.addRequiredParam<std::vector<UserObjectName>>(
      "uo_slip_rates",
      "List of names of user objects that define the slip rates for this material.");
  params.addRequiredParam<std::vector<UserObjectName>>(
      "uo_slip_resistances",
      "List of names of user objects that define the slip resistances for this material.");
  params.addRequiredParam<std::vector<UserObjectName>>(
      "uo_state_vars",
      "List of names of user objects that define the state variable for this material.");
  params.addRequiredParam<std::vector<UserObjectName>>(
      "uo_state_var_evol_rate_comps",
      "List of names of user objects that define the state "
      "variable evolution rate components for this material.");
  return params;
}

Member Data Documentation

_abs_tol

Real FiniteStrainUObasedCP::_abs_tol

protected

Stress residual equation absolute tolerance.

Definition at line 203 of file FiniteStrainUObasedCP.h.

Referenced by solveStress().

_base_name

const std::string ComputeGeneralStressBase::_base_name

protectedinherited

Base name prepended to all material property names to allow for multi-material systems.

Definition at line 43 of file ComputeGeneralStressBase.h.

Referenced by ComputeCosseratLinearElasticStress::initialSetup(), and ComputeLinearElasticStress::initialSetup().

_crysrot

const MaterialProperty<RankTwoTensor>& FiniteStrainUObasedCP::_crysrot

protected

Crystal rotation.

Definition at line 257 of file FiniteStrainUObasedCP.h.

Referenced by postSolveQp().

_deformation_gradient

const MaterialProperty<RankTwoTensor>& FiniteStrainUObasedCP::_deformation_gradient

protected

Definition at line 253 of file FiniteStrainUObasedCP.h.

Referenced by computeQpStress(), and postSolveQp().

_deformation_gradient_old

const MaterialProperty<RankTwoTensor>& FiniteStrainUObasedCP::_deformation_gradient_old

protected

Definition at line 254 of file FiniteStrainUObasedCP.h.

Referenced by computeQpStress().

_delta_dfgrd

RankTwoTensor FiniteStrainUObasedCP::_delta_dfgrd

protected

Used for substepping; Uniformly divides the increment in deformation gradient.

Definition at line 268 of file FiniteStrainUObasedCP.h.

Referenced by computeQpStress(), and FiniteStrainUObasedCP().

_dfgrd_scale_factor

Real FiniteStrainUObasedCP::_dfgrd_scale_factor

protected

Scales the substepping increment to obtain deformation gradient at a substep iteration.

Definition at line 270 of file FiniteStrainUObasedCP.h.

_dfgrd_tmp

RankTwoTensor FiniteStrainUObasedCP::_dfgrd_tmp

protected

Definition at line 259 of file FiniteStrainUObasedCP.h.

Referenced by calcJacobian(), calcResidual(), and computeQpStress().

_dfgrd_tmp_old

RankTwoTensor FiniteStrainUObasedCP::_dfgrd_tmp_old

protected

Definition at line 268 of file FiniteStrainUObasedCP.h.

Referenced by computeQpStress().

_elastic_strain

MaterialProperty<RankTwoTensor>& ComputeGeneralStressBase::_elastic_strain

protectedinherited

Elastic strain material property.

Definition at line 50 of file ComputeGeneralStressBase.h.

Referenced by ComputeSmearedCrackingStress::computeCrackStrainAndOrientation(), ComputeMultiPlasticityStress::computeQpStress(), FiniteStrainPlasticMaterial::computeQpStress(), ComputeCosseratLinearElasticStress::computeQpStress(), ComputeFiniteStrainElasticStress::computeQpStress(), ComputeLinearElasticStress::computeQpStress(), ComputeLinearViscoelasticStress::computeQpStress(), ComputeSmearedCrackingStress::computeQpStress(), ComputeMultipleInelasticStressBase::computeQpStressIntermediateConfiguration(), ComputeMultipleInelasticStressBase::finiteStrainRotation(), and ComputeGeneralStressBase::initQpStatefulProperties().

_elasticity_tensor

const MaterialProperty<RankFourTensor>& FiniteStrainUObasedCP::_elasticity_tensor

protected

Elasticity tensor material property.

Definition at line 252 of file FiniteStrainUObasedCP.h.

Referenced by calcJacobian(), calcResidual(), elasticTangentModuli(), and elastoPlasticTangentModuli().

_elasticity_tensor_name

const std::string FiniteStrainUObasedCP::_elasticity_tensor_name

protected

Name of the elasticity tensor material property.

Definition at line 250 of file FiniteStrainUObasedCP.h.

_err_tol

bool FiniteStrainUObasedCP::_err_tol

protected

Flag to check whether convergence is achieved.

Definition at line 265 of file FiniteStrainUObasedCP.h.

Referenced by calcResidJacob(), calcResidual(), computeQpStress(), FiniteStrainUObasedCP(), getSlipRates(), solveQp(), solveStatevar(), solveStress(), and updateSlipSystemResistanceAndStateVariable().

_extra_stress

const MaterialProperty<RankTwoTensor>& ComputeGeneralStressBase::_extra_stress

protectedinherited

Extra stress tensor.

Definition at line 53 of file ComputeGeneralStressBase.h.

Referenced by ComputeGeneralStressBase::computeQpProperties().

_fe

RankTwoTensor FiniteStrainUObasedCP::_fe

protected

Definition at line 260 of file FiniteStrainUObasedCP.h.

Referenced by calcJacobian(), calcResidual(), elastoPlasticTangentModuli(), and postSolveQp().

_flow_direction

std::vector<MaterialProperty<std::vector<RankTwoTensor> > *> FiniteStrainUObasedCP::_flow_direction

protected

Definition at line 262 of file FiniteStrainUObasedCP.h.

Referenced by calcJacobian(), calcResidual(), computeQpStress(), FiniteStrainUObasedCP(), and initQpStatefulProperties().

_fp

MaterialProperty<RankTwoTensor>& FiniteStrainUObasedCP::_fp

protected

Definition at line 241 of file FiniteStrainUObasedCP.h.

Referenced by initQpStatefulProperties(), and postSolveStress().

_fp_inv

RankTwoTensor FiniteStrainUObasedCP::_fp_inv

protected

Definition at line 260 of file FiniteStrainUObasedCP.h.

Referenced by calcResidual(), elastoPlasticTangentModuli(), postSolveStatevar(), postSolveStress(), and preSolveStatevar().

_fp_old

const MaterialProperty<RankTwoTensor>& FiniteStrainUObasedCP::_fp_old

protected

Definition at line 242 of file FiniteStrainUObasedCP.h.

Referenced by preSolveQp().

_fp_old_inv

RankTwoTensor FiniteStrainUObasedCP::_fp_old_inv

protected

Definition at line 260 of file FiniteStrainUObasedCP.h.

Referenced by calcJacobian(), calcResidual(), postSolveStatevar(), preSolveQp(), and preSolveStatevar().

_initial_stress_fcn

std::vector<const Function *> ComputeGeneralStressBase::_initial_stress_fcn

protectedinherited

initial stress components

Definition at line 56 of file ComputeGeneralStressBase.h.

_jac

RankFourTensor FiniteStrainUObasedCP::_jac

protected

Jacobian tensor.

Definition at line 213 of file FiniteStrainUObasedCP.h.

Referenced by calcJacobian(), and solveStress().

_Jacobian_mult

MaterialProperty<RankFourTensor>& ComputeGeneralStressBase::_Jacobian_mult

protectedinherited

derivative of stress w.r.t. strain (_dstress_dstrain)

Definition at line 59 of file ComputeGeneralStressBase.h.

Referenced by ComputeStrainIncrementBasedStress::computeQpJacobian(), FiniteStrainHyperElasticViscoPlastic::computeQpJacobian(), ComputeMultipleInelasticStressBase::computeQpJacobianMult(), ComputeMultipleInelasticCosseratStress::computeQpJacobianMult(), ComputeMultiPlasticityStress::computeQpStress(), FiniteStrainPlasticMaterial::computeQpStress(), ComputeCosseratLinearElasticStress::computeQpStress(), ComputeFiniteStrainElasticStress::computeQpStress(), ComputeLinearElasticPFFractureStress::computeQpStress(), ComputeLinearElasticStress::computeQpStress(), ComputeLinearViscoelasticStress::computeQpStress(), ComputeMultipleInelasticStressBase::computeQpStress(), ComputeSmearedCrackingStress::computeQpStress(), ComputeMultipleCrystalPlasticityStress::computeQpStress(), ComputeMultipleInelasticStressBase::computeQpStressIntermediateConfiguration(), ComputeLinearElasticPFFractureStress::computeStrainSpectral(), ComputeLinearElasticPFFractureStress::computeStrainVolDev(), ComputeLinearElasticPFFractureStress::computeStressSpectral(), elasticTangentModuli(), elastoPlasticTangentModuli(), ComputeMultipleInelasticStressBase::finiteStrainRotation(), ComputeMultiPlasticityStress::postReturnMap(), FiniteStrainCrystalPlasticity::postSolveQp(), FiniteStrainCrystalPlasticity::preSolveQp(), and ComputeMultipleInelasticStressBase::updateQpStateSingleModel().

_lag_e

MaterialProperty<RankTwoTensor>& FiniteStrainUObasedCP::_lag_e

protected

Definition at line 245 of file FiniteStrainUObasedCP.h.

Referenced by initQpStatefulProperties(), and postSolveQp().

_lsrch_max_iter

unsigned int FiniteStrainUObasedCP::_lsrch_max_iter

protected

Line search bisection method maximum iteration number.

Definition at line 236 of file FiniteStrainUObasedCP.h.

Referenced by lineSearchUpdate().

_lsrch_method

MooseEnum FiniteStrainUObasedCP::_lsrch_method

protected

Line search method.

Definition at line 239 of file FiniteStrainUObasedCP.h.

Referenced by lineSearchUpdate().

_lsrch_tol

Real FiniteStrainUObasedCP::_lsrch_tol

protected

Line search bisection method tolerance.

Definition at line 233 of file FiniteStrainUObasedCP.h.

Referenced by lineSearchUpdate().

_mat_prop_slip_rates

std::vector<MaterialProperty<std::vector<Real> > *> FiniteStrainUObasedCP::_mat_prop_slip_rates

protected

Slip rates material property.

Definition at line 165 of file FiniteStrainUObasedCP.h.

Referenced by calcResidual(), FiniteStrainUObasedCP(), getSlipRates(), and initQpStatefulProperties().

_mat_prop_slip_resistances

std::vector<MaterialProperty<std::vector<Real> > *> FiniteStrainUObasedCP::_mat_prop_slip_resistances

protected

Slip resistance material property.

Definition at line 168 of file FiniteStrainUObasedCP.h.

Referenced by FiniteStrainUObasedCP(), initQpStatefulProperties(), preSolveStatevar(), and updateSlipSystemResistanceAndStateVariable().

_mat_prop_state_var_evol_rate_comps

std::vector<MaterialProperty<std::vector<Real> > *> FiniteStrainUObasedCP::_mat_prop_state_var_evol_rate_comps

protected

State variable evolution rate component material property.

Definition at line 177 of file FiniteStrainUObasedCP.h.

Referenced by FiniteStrainUObasedCP(), initQpStatefulProperties(), and updateSlipSystemResistanceAndStateVariable().

_mat_prop_state_vars

std::vector<MaterialProperty<std::vector<Real> > *> FiniteStrainUObasedCP::_mat_prop_state_vars

protected

State variable material property.

Definition at line 171 of file FiniteStrainUObasedCP.h.

Referenced by FiniteStrainUObasedCP(), initQpStatefulProperties(), isStateVariablesConverged(), postSolveStatevar(), preSolveQp(), preSolveStatevar(), and updateSlipSystemResistanceAndStateVariable().

_mat_prop_state_vars_old

std::vector<const MaterialProperty<std::vector<Real> > *> FiniteStrainUObasedCP::_mat_prop_state_vars_old

protected

Old state variable material property.

Definition at line 174 of file FiniteStrainUObasedCP.h.

Referenced by computeQpStress(), and FiniteStrainUObasedCP().

_max_substep_iter

unsigned int FiniteStrainUObasedCP::_max_substep_iter

protected

Maximum number of substep iterations.

Definition at line 224 of file FiniteStrainUObasedCP.h.

Referenced by computeQpStress().

_maxiter

unsigned int FiniteStrainUObasedCP::_maxiter

protected

Maximum number of iterations for stress update.

Definition at line 216 of file FiniteStrainUObasedCP.h.

Referenced by solveStress().

_maxiterg

unsigned int FiniteStrainUObasedCP::_maxiterg

protected

Maximum number of iterations for internal variable update.

Definition at line 218 of file FiniteStrainUObasedCP.h.

Referenced by solveStatevar().

_mechanical_strain

const MaterialProperty<RankTwoTensor>& ComputeGeneralStressBase::_mechanical_strain

protectedinherited

Mechanical strain material property.

Definition at line 46 of file ComputeGeneralStressBase.h.

Referenced by ComputeStrainIncrementBasedStress::computeQpStress(), FiniteStrainPlasticMaterial::computeQpStress(), ComputeCosseratLinearElasticStress::computeQpStress(), ComputeFiniteStrainElasticStress::computeQpStress(), ComputeLinearElasticPFFractureStress::computeQpStress(), ComputeLinearElasticStress::computeQpStress(), ComputeLinearViscoelasticStress::computeQpStress(), ComputeLinearElasticPFFractureStress::computeStrainSpectral(), ComputeLinearElasticPFFractureStress::computeStrainVolDev(), and ComputeLinearElasticPFFractureStress::computeStressSpectral().

_min_lsrch_step

Real FiniteStrainUObasedCP::_min_lsrch_step

protected

Minimum line search step size.

Definition at line 230 of file FiniteStrainUObasedCP.h.

Referenced by lineSearchUpdate().

_num_uo_slip_rates

unsigned int FiniteStrainUObasedCP::_num_uo_slip_rates

protected

Number of slip rate user objects.

Definition at line 180 of file FiniteStrainUObasedCP.h.

Referenced by calcJacobian(), calcResidual(), computeQpStress(), FiniteStrainUObasedCP(), getSlipRates(), and initQpStatefulProperties().

_num_uo_slip_resistances

unsigned int FiniteStrainUObasedCP::_num_uo_slip_resistances

protected

Number of slip resistance user objects.

Definition at line 183 of file FiniteStrainUObasedCP.h.

Referenced by FiniteStrainUObasedCP(), initQpStatefulProperties(), preSolveStatevar(), and updateSlipSystemResistanceAndStateVariable().

_num_uo_state_var_evol_rate_comps

unsigned int FiniteStrainUObasedCP::_num_uo_state_var_evol_rate_comps

protected

Number of state variable evolution rate component user objects.

Definition at line 189 of file FiniteStrainUObasedCP.h.

Referenced by FiniteStrainUObasedCP(), initQpStatefulProperties(), and updateSlipSystemResistanceAndStateVariable().

_num_uo_state_vars

unsigned int FiniteStrainUObasedCP::_num_uo_state_vars

protected

Number of state variable user objects.

Definition at line 186 of file FiniteStrainUObasedCP.h.

Referenced by computeQpStress(), FiniteStrainUObasedCP(), initQpStatefulProperties(), isStateVariablesConverged(), postSolveStatevar(), preSolveQp(), preSolveStatevar(), and updateSlipSystemResistanceAndStateVariable().

_pk2

MaterialProperty<RankTwoTensor>& FiniteStrainUObasedCP::_pk2

protected

Definition at line 243 of file FiniteStrainUObasedCP.h.

Referenced by calcResidual(), elastoPlasticTangentModuli(), initQpStatefulProperties(), lineSearchUpdate(), postSolveQp(), preSolveQp(), and solveStress().

_pk2_old

const MaterialProperty<RankTwoTensor>& FiniteStrainUObasedCP::_pk2_old

protected

Definition at line 244 of file FiniteStrainUObasedCP.h.

Referenced by preSolveQp().

_resid

RankTwoTensor FiniteStrainUObasedCP::_resid

protected

Residual tensor.

Definition at line 210 of file FiniteStrainUObasedCP.h.

Referenced by calcResidual(), lineSearchUpdate(), and solveStress().

_rtol

Real FiniteStrainUObasedCP::_rtol

protected

Stress residual equation relative tolerance.

Definition at line 201 of file FiniteStrainUObasedCP.h.

Referenced by solveStress().

_state_vars_old

std::vector<std::vector<Real> > FiniteStrainUObasedCP::_state_vars_old

protected

Local state variable.

Definition at line 192 of file FiniteStrainUObasedCP.h.

Referenced by computeQpStress(), FiniteStrainUObasedCP(), initQpStatefulProperties(), and preSolveQp().

_state_vars_old_stored

std::vector<std::vector<Real> > FiniteStrainUObasedCP::_state_vars_old_stored

protected

Local stored state variable (for sub-stepping)

Definition at line 195 of file FiniteStrainUObasedCP.h.

Referenced by FiniteStrainUObasedCP(), initQpStatefulProperties(), isStateVariablesConverged(), postSolveStatevar(), preSolveQp(), preSolveStatevar(), and updateSlipSystemResistanceAndStateVariable().

_state_vars_prev

std::vector<std::vector<Real> > FiniteStrainUObasedCP::_state_vars_prev

protected

Local old state variable.

Definition at line 198 of file FiniteStrainUObasedCP.h.

Referenced by FiniteStrainUObasedCP(), initQpStatefulProperties(), isStateVariablesConverged(), and updateSlipSystemResistanceAndStateVariable().

_stol

Real FiniteStrainUObasedCP::_stol

protected

Internal variable update equation tolerance.

Definition at line 205 of file FiniteStrainUObasedCP.h.

Referenced by isStateVariablesConverged().

_stress

MaterialProperty<RankTwoTensor>& ComputeGeneralStressBase::_stress

protectedinherited

Stress material property.

Definition at line 48 of file ComputeGeneralStressBase.h.

Referenced by ComputeMultipleInelasticStressBase::computeAdmissibleState(), ComputeMultipleInelasticCosseratStress::computeAdmissibleState(), ComputeGeneralStressBase::computeQpProperties(), ComputeMultiPlasticityStress::computeQpStress(), ComputeStrainIncrementBasedStress::computeQpStress(), FiniteStrainPlasticMaterial::computeQpStress(), AbaqusUMATStress::computeQpStress(), ComputeCosseratLinearElasticStress::computeQpStress(), ComputeFiniteStrainElasticStress::computeQpStress(), ComputeLinearElasticPFFractureStress::computeQpStress(), ComputeLinearElasticStress::computeQpStress(), ComputeLinearViscoelasticStress::computeQpStress(), ComputeMultipleInelasticStressBase::computeQpStress(), ComputeSmearedCrackingStress::computeQpStress(), ComputeMultipleCrystalPlasticityStress::computeQpStress(), ComputeMultipleInelasticStressBase::computeQpStressIntermediateConfiguration(), ComputeLinearElasticPFFractureStress::computeStrainSpectral(), ComputeLinearElasticPFFractureStress::computeStrainVolDev(), ComputeLinearElasticPFFractureStress::computeStressSpectral(), ComputeCreepPlasticityStress::computeTangentOperators(), ComputeMultipleInelasticStressBase::finiteStrainRotation(), FiniteStrainCrystalPlasticity::initQpStatefulProperties(), initQpStatefulProperties(), FiniteStrainHyperElasticViscoPlastic::initQpStatefulProperties(), ComputeGeneralStressBase::initQpStatefulProperties(), ComputeMultiPlasticityStress::postReturnMap(), FiniteStrainCrystalPlasticity::postSolveQp(), postSolveQp(), FiniteStrainHyperElasticViscoPlastic::postSolveQp(), ComputeSmearedCrackingStress::updateCrackingStateAndStress(), ComputeCreepPlasticityStress::updateQpState(), ComputeMultipleInelasticStress::updateQpState(), and ComputeMultipleInelasticStressBase::updateQpStateSingleModel().

_substep_dt

Real FiniteStrainUObasedCP::_substep_dt

protected

Current substep size.

Definition at line 272 of file FiniteStrainUObasedCP.h.

Referenced by calcJacobian(), calcResidual(), computeQpStress(), FiniteStrainUObasedCP(), getSlipRates(), and updateSlipSystemResistanceAndStateVariable().

_tan_mod_type

MooseEnum FiniteStrainUObasedCP::_tan_mod_type

protected

Type of tangent moduli calculation.

Definition at line 221 of file FiniteStrainUObasedCP.h.

Referenced by calcTangentModuli().

_tau

DenseVector<Real> FiniteStrainUObasedCP::_tau

protected

Definition at line 261 of file FiniteStrainUObasedCP.h.

_uo_slip_rates

std::vector<const CrystalPlasticitySlipRate *> FiniteStrainUObasedCP::_uo_slip_rates

protected

User objects that define the slip rate.

Definition at line 153 of file FiniteStrainUObasedCP.h.

Referenced by calcJacobian(), calcResidual(), computeQpStress(), FiniteStrainUObasedCP(), getSlipRates(), and initQpStatefulProperties().

_uo_slip_resistances

std::vector<const CrystalPlasticitySlipResistance *> FiniteStrainUObasedCP::_uo_slip_resistances

protected

User objects that define the slip resistance.

Definition at line 156 of file FiniteStrainUObasedCP.h.

Referenced by FiniteStrainUObasedCP(), initQpStatefulProperties(), preSolveStatevar(), and updateSlipSystemResistanceAndStateVariable().

_uo_state_var_evol_rate_comps

std::vector<const CrystalPlasticityStateVarRateComponent *> FiniteStrainUObasedCP::_uo_state_var_evol_rate_comps

protected

User objects that define the state variable evolution rate component.

Definition at line 162 of file FiniteStrainUObasedCP.h.

Referenced by FiniteStrainUObasedCP(), initQpStatefulProperties(), and updateSlipSystemResistanceAndStateVariable().

_uo_state_vars

std::vector<const CrystalPlasticityStateVariable *> FiniteStrainUObasedCP::_uo_state_vars

protected

User objects that define the state variable.

Definition at line 159 of file FiniteStrainUObasedCP.h.

Referenced by FiniteStrainUObasedCP(), initQpStatefulProperties(), and updateSlipSystemResistanceAndStateVariable().

_update_rot

MaterialProperty<RankTwoTensor>& FiniteStrainUObasedCP::_update_rot

protected

Definition at line 246 of file FiniteStrainUObasedCP.h.

Referenced by initQpStatefulProperties(), and postSolveQp().

_update_rot_old

const MaterialProperty<RankTwoTensor>& FiniteStrainUObasedCP::_update_rot_old

protected

Definition at line 247 of file FiniteStrainUObasedCP.h.

_use_line_search

bool FiniteStrainUObasedCP::_use_line_search

protected

Flag to activate line serach.

Definition at line 227 of file FiniteStrainUObasedCP.h.

Referenced by solveStress().

_zero_tol

Real FiniteStrainUObasedCP::_zero_tol

protected

Residual tolerance when variable value is zero. Default 1e-12.

Definition at line 207 of file FiniteStrainUObasedCP.h.

Referenced by isStateVariablesConverged().

---

The documentation for this class was generated from the following files:

• solid_mechanics/include/materials/crystal_plasticity/FiniteStrainUObasedCP.h
• solid_mechanics/src/materials/crystal_plasticity/FiniteStrainUObasedCP.C