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:

[legend]

Public Member Functions
	FiniteStrainUObasedCP (const InputParameters &parameters)

Static Public Member Functions
static InputParameters	validParams ()

Protected Member Functions
virtual void	computeQpStress ()
	updates the stress at a quadrature point. More...

virtual void	initQpStatefulProperties ()
	initializes the stateful properties such as stress, plastic deformation gradient, slip system resistances, etc. More...

virtual void	calcResidJacob ()
	calls the residual and jacobian functions used in the stress update algorithm. More...

virtual void	updateSlipSystemResistanceAndStateVariable ()
	updates the slip system resistances and state variables. More...

virtual void	preSolveQp ()
	set variables for stress and internal variable solve. More...

virtual void	solveQp ()
	solve stress and internal variables. More...

virtual void	postSolveQp ()
	update stress and internal variable after solve. More...

virtual void	preSolveStatevar ()
	set variables for internal variable solve. More...

virtual void	solveStatevar ()
	solve internal variables. More...

virtual void	postSolveStatevar ()
	update internal variable after solve. More...

virtual void	preSolveStress ()
	set variables for stress solve. More...

virtual void	solveStress ()
	solves for stress, updates plastic deformation gradient. More...

virtual void	postSolveStress ()
	update stress and plastic deformation gradient after solve. More...

virtual void	calcResidual ()
	calculate stress residual. More...

virtual void	calcJacobian ()
	calculate jacobian. More...

virtual void	getSlipRates ()
	updates the slip rates. More...

virtual void	calcTangentModuli ()
	calculate the tangent moduli for preconditioner. More...

virtual void	elasticTangentModuli ()
	calculate the elastic tangent moduli for preconditioner. More...

virtual void	elastoPlasticTangentModuli ()
	calculate the exact tangent moduli for preconditioner. More...

bool	lineSearchUpdate (const Real rnorm_prev, const RankTwoTensor)
	performs the line search update More...

virtual bool	isStateVariablesConverged ()
	evaluates convergence of state variables. More...

virtual void	computeQpProperties () override

Protected Attributes
std::vector< const CrystalPlasticitySlipRate * >	_uo_slip_rates
	User objects that define the slip rate. More...

std::vector< const CrystalPlasticitySlipResistance * >	_uo_slip_resistances
	User objects that define the slip resistance. More...

std::vector< const CrystalPlasticityStateVariable * >	_uo_state_vars
	User objects that define the state variable. More...

std::vector< const CrystalPlasticityStateVarRateComponent * >	_uo_state_var_evol_rate_comps
	User objects that define the state variable evolution rate component. More...

std::vector< MaterialProperty< std::vector< Real > > * >	_mat_prop_slip_rates
	Slip rates material property. More...

std::vector< MaterialProperty< std::vector< Real > > * >	_mat_prop_slip_resistances
	Slip resistance material property. More...

std::vector< MaterialProperty< std::vector< Real > > * >	_mat_prop_state_vars
	State variable material property. More...

std::vector< const MaterialProperty< std::vector< Real > > * >	_mat_prop_state_vars_old
	Old state variable material property. More...

std::vector< MaterialProperty< std::vector< Real > > * >	_mat_prop_state_var_evol_rate_comps
	State variable evolution rate component material property. More...

unsigned int	_num_uo_slip_rates
	Number of slip rate user objects. More...

unsigned int	_num_uo_slip_resistances
	Number of slip resistance user objects. More...

unsigned int	_num_uo_state_vars
	Number of state variable user objects. More...

unsigned int	_num_uo_state_var_evol_rate_comps
	Number of state variable evolution rate component user objects. More...

std::vector< std::vector< Real > >	_state_vars_old
	Local state variable. More...

std::vector< std::vector< Real > >	_state_vars_old_stored
	Local stored state variable (for sub-stepping) More...

std::vector< std::vector< Real > >	_state_vars_prev
	Local old state variable. More...

Real	_rtol
	Stress residual equation relative tolerance. More...

Real	_abs_tol
	Stress residual equation absolute tolerance. More...

Real	_stol
	Internal variable update equation tolerance. More...

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

RankTwoTensor	_resid
	Residual tensor. More...

RankFourTensor	_jac
	Jacobian tensor. More...

unsigned int	_maxiter
	Maximum number of iterations for stress update. More...

unsigned int	_maxiterg
	Maximum number of iterations for internal variable update. More...

MooseEnum	_tan_mod_type
	Type of tangent moduli calculation. More...

unsigned int	_max_substep_iter
	Maximum number of substep iterations. More...

bool	_use_line_search
	Flag to activate line serach. More...

Real	_min_lsrch_step
	Minimum line search step size. More...

Real	_lsrch_tol
	Line search bisection method tolerance. More...

unsigned int	_lsrch_max_iter
	Line search bisection method maximum iteration number. More...

MooseEnum	_lsrch_method
	Line search method. More...

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. More...

const MaterialProperty< RankFourTensor > &	_elasticity_tensor
	Elasticity tensor material property. More...

const MaterialProperty< RankTwoTensor > &	_deformation_gradient

const MaterialProperty< RankTwoTensor > &	_deformation_gradient_old

const MaterialProperty< RankTwoTensor > &	_crysrot
	Crystal rotation. More...

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. More...

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

RankTwoTensor	_dfgrd_tmp_old

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

const std::string	_base_name
	Base name prepended to all material property names to allow for multi-material systems. More...

const MaterialProperty< RankTwoTensor > &	_mechanical_strain
	Mechanical strain material property. More...

MaterialProperty< RankTwoTensor > &	_stress
	Stress material property. More...

MaterialProperty< RankTwoTensor > &	_elastic_strain
	Elastic strain material property. More...

const MaterialProperty< RankTwoTensor > &	_extra_stress
	Extra stress tensor. More...

std::vector< const Function * >	_initial_stress_fcn
	initial stress components More...

MaterialProperty< RankFourTensor > &	_Jacobian_mult
	derivative of stress w.r.t. strain (_dstress_dstrain) More...

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 38 of file FiniteStrainUObasedCP.h.

Constructor & Destructor Documentation

◆ FiniteStrainUObasedCP()

FiniteStrainUObasedCP::FiniteStrainUObasedCP ( const InputParameters & parameters )

Definition at line 66 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]);
   }
 }

Member Function Documentation

◆ calcJacobian()

void FiniteStrainUObasedCP::calcJacobian ( )

protectedvirtual

calculate jacobian.

Definition at line 531 of file FiniteStrainUObasedCP.C.

 {
   RankFourTensor dfedfpinv, deedfe, dfpinvdpk2;
  
   for (unsigned int i = 0; i < LIBMESH_DIM; ++i)
     for (unsigned int j = 0; j < LIBMESH_DIM; ++j)
       for (unsigned int k = 0; k < LIBMESH_DIM; ++k)
         dfedfpinv(i, j, k, j) = _dfgrd_tmp(i, k);
  
   for (unsigned int i = 0; i < LIBMESH_DIM; ++i)
     for (unsigned int j = 0; j < LIBMESH_DIM; ++j)
       for (unsigned int k = 0; k < LIBMESH_DIM; ++k)
       {
         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, _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];
     }
  
     for (unsigned int j = 0; j < nss; j++)
       dfpinvdpk2 += (dfpinvdslip[j] * dslipdtau[j] * _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 483 of file FiniteStrainUObasedCP.C.

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

Referenced by solveStress().

◆ calcResidual()

void FiniteStrainUObasedCP::calcResidual ( )

protectedvirtual

calculate stress residual.

Definition at line 505 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] * _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 569 of file FiniteStrainUObasedCP.C.

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

Referenced by postSolveQp().

◆ computeQpProperties()

void ComputeStressBase::computeQpProperties ( )

overrideprotectedvirtualinherited

Definition at line 50 of file ComputeStressBase.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 ComputeStressBase.

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;
   // Store original _dt; Reset at the end of solve
   Real dt_original = _dt;
  
   _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 is 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();
  
     _dt = dt_original / 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);
  
   _dt = dt_original;
  
   postSolveQp();
 }

◆ elasticTangentModuli()

void FiniteStrainUObasedCP::elasticTangentModuli ( )

protectedvirtual

calculate the elastic tangent moduli for preconditioner.

Definition at line 625 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 582 of file FiniteStrainUObasedCP.C.

 {
   RankFourTensor tan_mod;
   RankTwoTensor pk2fet, fepk2;
   RankFourTensor deedfe, dsigdpk2dfe, dfedf;
  
   // Fill in the matrix stiffness material property
   for (unsigned int i = 0; i < LIBMESH_DIM; ++i)
     for (unsigned int j = 0; j < LIBMESH_DIM; ++j)
       for (unsigned int k = 0; k < LIBMESH_DIM; ++k)
       {
         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;
       }
  
   dsigdpk2dfe = _fe.mixedProductIkJl(_fe) * _elasticity_tensor[_qp] * deedfe;
  
   pk2fet = _pk2[_qp] * _fe.transpose();
   fepk2 = _fe * _pk2[_qp];
  
   for (unsigned int i = 0; i < LIBMESH_DIM; ++i)
     for (unsigned int j = 0; j < LIBMESH_DIM; ++j)
       for (unsigned int l = 0; l < LIBMESH_DIM; ++l)
       {
         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 (unsigned int i = 0; i < LIBMESH_DIM; ++i)
     for (unsigned int j = 0; j < LIBMESH_DIM; ++j)
       for (unsigned int l = 0; l < LIBMESH_DIM; ++l)
         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 492 of file FiniteStrainUObasedCP.C.

 {
   for (unsigned int i = 0; i < _num_uo_slip_rates; ++i)
   {
     if (!_uo_slip_rates[i]->calcSlipRate(_qp, _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 ComputeStressBase.

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 350 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]); // Calculate increment size
       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 632 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 284 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 372 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 455 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 264 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 303 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 381 of file FiniteStrainUObasedCP.C.

382 {

383 }

Referenced by solveStatevar().

◆ solveQp()

void FiniteStrainUObasedCP::solveQp ( )

protectedvirtual

solve stress and internal variables.

Definition at line 274 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 315 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 386 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);
 #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 461 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, _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 23 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 208 of file FiniteStrainUObasedCP.h.

Referenced by solveStress().

◆ _base_name

const std::string ComputeStressBase::_base_name

protectedinherited

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

Definition at line 45 of file ComputeStressBase.h.

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

◆ _crysrot

const MaterialProperty<RankTwoTensor>& FiniteStrainUObasedCP::_crysrot

protected

Crystal rotation.

Definition at line 262 of file FiniteStrainUObasedCP.h.

Referenced by postSolveQp().

◆ _deformation_gradient

const MaterialProperty<RankTwoTensor>& FiniteStrainUObasedCP::_deformation_gradient

protected

Definition at line 258 of file FiniteStrainUObasedCP.h.

Referenced by computeQpStress(), and postSolveQp().

◆ _deformation_gradient_old

const MaterialProperty<RankTwoTensor>& FiniteStrainUObasedCP::_deformation_gradient_old

protected

Definition at line 259 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 273 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 275 of file FiniteStrainUObasedCP.h.

◆ _dfgrd_tmp

RankTwoTensor FiniteStrainUObasedCP::_dfgrd_tmp

protected

Definition at line 264 of file FiniteStrainUObasedCP.h.

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

◆ _dfgrd_tmp_old

RankTwoTensor FiniteStrainUObasedCP::_dfgrd_tmp_old

protected

Definition at line 273 of file FiniteStrainUObasedCP.h.

Referenced by computeQpStress().

◆ _elastic_strain

MaterialProperty<RankTwoTensor>& ComputeStressBase::_elastic_strain

protectedinherited

Elastic strain material property.

Definition at line 52 of file ComputeStressBase.h.

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

◆ _elasticity_tensor

const MaterialProperty<RankFourTensor>& FiniteStrainUObasedCP::_elasticity_tensor

protected

Elasticity tensor material property.

Definition at line 257 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 255 of file FiniteStrainUObasedCP.h.

◆ _err_tol

bool FiniteStrainUObasedCP::_err_tol

protected

Flag to check whether convergence is achieved.

Definition at line 270 of file FiniteStrainUObasedCP.h.

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

◆ _extra_stress

const MaterialProperty<RankTwoTensor>& ComputeStressBase::_extra_stress

protectedinherited

Extra stress tensor.

Definition at line 55 of file ComputeStressBase.h.

Referenced by ComputeStressBase::computeQpProperties().

◆ _fe

RankTwoTensor FiniteStrainUObasedCP::_fe

protected

Definition at line 265 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 267 of file FiniteStrainUObasedCP.h.

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

◆ _fp

MaterialProperty<RankTwoTensor>& FiniteStrainUObasedCP::_fp

protected

Definition at line 246 of file FiniteStrainUObasedCP.h.

Referenced by initQpStatefulProperties(), and postSolveStress().

◆ _fp_inv

RankTwoTensor FiniteStrainUObasedCP::_fp_inv

protected

Definition at line 265 of file FiniteStrainUObasedCP.h.

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

◆ _fp_old

const MaterialProperty<RankTwoTensor>& FiniteStrainUObasedCP::_fp_old

protected

Definition at line 247 of file FiniteStrainUObasedCP.h.

Referenced by preSolveQp().

◆ _fp_old_inv

RankTwoTensor FiniteStrainUObasedCP::_fp_old_inv

protected

Definition at line 265 of file FiniteStrainUObasedCP.h.

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

◆ _initial_stress_fcn

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

protectedinherited

initial stress components

Definition at line 58 of file ComputeStressBase.h.

◆ _jac

RankFourTensor FiniteStrainUObasedCP::_jac

protected

Jacobian tensor.

Definition at line 218 of file FiniteStrainUObasedCP.h.

Referenced by calcJacobian(), and solveStress().

◆ _Jacobian_mult

MaterialProperty<RankFourTensor>& ComputeStressBase::_Jacobian_mult

protectedinherited

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

Definition at line 61 of file ComputeStressBase.h.

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

◆ _lag_e

MaterialProperty<RankTwoTensor>& FiniteStrainUObasedCP::_lag_e

protected

Definition at line 250 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 241 of file FiniteStrainUObasedCP.h.

Referenced by lineSearchUpdate().

◆ _lsrch_method

MooseEnum FiniteStrainUObasedCP::_lsrch_method

protected

Line search method.

Definition at line 244 of file FiniteStrainUObasedCP.h.

Referenced by lineSearchUpdate().

◆ _lsrch_tol

Real FiniteStrainUObasedCP::_lsrch_tol

protected

Line search bisection method tolerance.

Definition at line 238 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 170 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 173 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 182 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 176 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 179 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 229 of file FiniteStrainUObasedCP.h.

Referenced by computeQpStress().

◆ _maxiter

unsigned int FiniteStrainUObasedCP::_maxiter

protected

Maximum number of iterations for stress update.

Definition at line 221 of file FiniteStrainUObasedCP.h.

Referenced by solveStress().

◆ _maxiterg

unsigned int FiniteStrainUObasedCP::_maxiterg

protected

Maximum number of iterations for internal variable update.

Definition at line 223 of file FiniteStrainUObasedCP.h.

Referenced by solveStatevar().

◆ _mechanical_strain

const MaterialProperty<RankTwoTensor>& ComputeStressBase::_mechanical_strain

protectedinherited

Mechanical strain material property.

Definition at line 48 of file ComputeStressBase.h.

Referenced by ComputeLinearElasticStress::computeQpStress(), ComputeFiniteStrainElasticStress::computeQpStress(), ComputeCosseratLinearElasticStress::computeQpStress(), FiniteStrainPlasticMaterial::computeQpStress(), ComputeLinearElasticPFFractureStress::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 235 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 185 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 188 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 194 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 191 of file FiniteStrainUObasedCP.h.

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

◆ _pk2

MaterialProperty<RankTwoTensor>& FiniteStrainUObasedCP::_pk2

protected

Definition at line 248 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 249 of file FiniteStrainUObasedCP.h.

Referenced by preSolveQp().

◆ _resid

RankTwoTensor FiniteStrainUObasedCP::_resid

protected

Residual tensor.

Definition at line 215 of file FiniteStrainUObasedCP.h.

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

◆ _rtol

Real FiniteStrainUObasedCP::_rtol

protected

Stress residual equation relative tolerance.

Definition at line 206 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 197 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 200 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 203 of file FiniteStrainUObasedCP.h.

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

◆ _stol

Real FiniteStrainUObasedCP::_stol

protected

Internal variable update equation tolerance.

Definition at line 210 of file FiniteStrainUObasedCP.h.

Referenced by isStateVariablesConverged().

◆ _stress

MaterialProperty<RankTwoTensor>& ComputeStressBase::_stress

protectedinherited

Stress material property.

Definition at line 50 of file ComputeStressBase.h.

◆ _tan_mod_type

MooseEnum FiniteStrainUObasedCP::_tan_mod_type

protected

Type of tangent moduli calculation.

Definition at line 226 of file FiniteStrainUObasedCP.h.

Referenced by calcTangentModuli().

◆ _tau

DenseVector<Real> FiniteStrainUObasedCP::_tau

protected

Definition at line 266 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 158 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 161 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 167 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 164 of file FiniteStrainUObasedCP.h.

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

◆ _update_rot

MaterialProperty<RankTwoTensor>& FiniteStrainUObasedCP::_update_rot

protected

Definition at line 251 of file FiniteStrainUObasedCP.h.

Referenced by initQpStatefulProperties(), and postSolveQp().

◆ _update_rot_old

const MaterialProperty<RankTwoTensor>& FiniteStrainUObasedCP::_update_rot_old

protected

Definition at line 252 of file FiniteStrainUObasedCP.h.

◆ _use_line_search

bool FiniteStrainUObasedCP::_use_line_search

protected

Flag to activate line serach.

Definition at line 232 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 212 of file FiniteStrainUObasedCP.h.

Referenced by isStateVariablesConverged().

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

tensor_mechanics/include/materials/FiniteStrainUObasedCP.h
tensor_mechanics/src/materials/FiniteStrainUObasedCP.C

Public Member Functions

Static Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ FiniteStrainUObasedCP()

Member Function Documentation

◆ calcJacobian()

◆ calcResidJacob()

◆ calcResidual()

◆ calcTangentModuli()

◆ computeQpProperties()

◆ computeQpStress()

◆ elasticTangentModuli()

◆ elastoPlasticTangentModuli()

◆ getSlipRates()

◆ initQpStatefulProperties()

◆ isStateVariablesConverged()

◆ lineSearchUpdate()

◆ postSolveQp()

◆ postSolveStatevar()

◆ postSolveStress()

◆ preSolveQp()

◆ preSolveStatevar()

◆ preSolveStress()

◆ solveQp()

◆ solveStatevar()

◆ solveStress()

◆ updateSlipSystemResistanceAndStateVariable()

◆ validParams()

Member Data Documentation

◆ _abs_tol

◆ _base_name

◆ _crysrot

◆ _deformation_gradient

◆ _deformation_gradient_old

◆ _delta_dfgrd

◆ _dfgrd_scale_factor

◆ _dfgrd_tmp

◆ _dfgrd_tmp_old

◆ _elastic_strain

◆ _elasticity_tensor

◆ _elasticity_tensor_name

◆ _err_tol

◆ _extra_stress

◆ _fe

◆ _flow_direction

◆ _fp

◆ _fp_inv

◆ _fp_old

◆ _fp_old_inv

◆ _initial_stress_fcn

◆ _jac

◆ _Jacobian_mult

◆ _lag_e

◆ _lsrch_max_iter

◆ _lsrch_method

◆ _lsrch_tol

◆ _mat_prop_slip_rates

◆ _mat_prop_slip_resistances

◆ _mat_prop_state_var_evol_rate_comps

◆ _mat_prop_state_vars

◆ _mat_prop_state_vars_old

◆ _max_substep_iter

◆ _maxiter

◆ _maxiterg

◆ _mechanical_strain

◆ _min_lsrch_step

◆ _num_uo_slip_rates

◆ _num_uo_slip_resistances

◆ _num_uo_state_var_evol_rate_comps

◆ _num_uo_state_vars

◆ _pk2

◆ _pk2_old

◆ _resid

◆ _rtol

◆ _state_vars_old

◆ _state_vars_old_stored

◆ _state_vars_prev