FiniteStrainCPSlipRateRes Class Reference

#include <FiniteStrainCPSlipRateRes.h>

Inheritance diagram for FiniteStrainCPSlipRateRes:

Public Types
typedef DerivativeMaterialPropertyNameInterface::SymbolName	SymbolName

Public Member Functions
	FiniteStrainCPSlipRateRes (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	solveStatevar ()
	This function solves internal variables.
virtual void	preSolveStress ()
	This function sets variable for internal variable solve.
virtual void	solveStress ()
	This function solves for stress, updates plastic deformation gradient.
virtual void	calcResidJacobSlipRate ()
	This function calculates residual and jacobian of slip rate.
virtual void	calcResidualSlipRate ()
	This function calculates residual of slip rate.
virtual void	calcJacobianSlipRate ()
	This function calculates jacobian of slip rate.
virtual void	getSlipIncrements ()
	This function updates the slip system resistances.
virtual void	calcDtauDsliprate ()
	This function calculates partial derivative of resolved shear stress with respect to split rate.
virtual void	calcDgssDsliprate ()
	This function calculates partial derivative of slip system resistances with respect to split rate.
void	calcUpdate ()
	This function calculates and updates the residual of slip rate.
virtual Real	calcResidNorm ()
	This function calculates the residual norm.
bool	lineSearchUpdateSlipRate (const Real, const DenseVector< Real > &)
	This function performs the line search update.
Real	calcResidDotProdUpdate (const DenseVector< Real > &)
	This function calculates the dot product of residual and update.
virtual void	computeQpStress ()
	This function updates the stress at a quadrature point.
virtual void	computeQpElasticityTensor ()
	This function updates the elasticity tensor at a quadrature point.
virtual void	initQpStatefulProperties ()
	This function initializes the stateful properties such as stress, plastic deformation gradient, slip system resistances, etc.
virtual void	calc_resid_jacob (RankTwoTensor &, RankFourTensor &)
	This function calls the residual and jacobian functions used in the stress update algorithm.
virtual void	update_slip_system_resistance ()
	This function updates the slip system resistances.
virtual void	updateGss ()
	This function updates the slip system resistances.
virtual void	getSlipSystems ()
	This function reads slip system from file - see test.
virtual void	assignSlipSysRes ()
	This function assign initial values of slip system resistances/internal variables read from getSlipSystems().
virtual void	readFileInitSlipSysRes ()
	This function read slip system resistances from file - see test.
virtual void	getInitSlipSysRes ()
	This function assign slip system resistances - see test.
virtual void	readFileFlowRateParams ()
	This function read flow rate parameters from file - see test.
virtual void	getFlowRateParams ()
	This function assign flow rate parameters - see test.
virtual void	readFileHardnessParams ()
	This function read hardness parameters from file.
virtual void	getHardnessParams ()
	This function assign flow rate parameters from .i file - see test.
virtual void	initSlipSysProps ()
	This function initializes slip system resistances.
virtual void	initAdditionalProps ()
	This function initializes additional parameters.
virtual void	preSolveQp ()
	This function set variables for stress and internal variable solve.
virtual void	solveQp ()
	This function solves stress and internal variables.
virtual void	postSolveQp ()
	This function update stress and internal variable after solve.
virtual void	preSolveStatevar ()
	This function set variables for internal variable solve.
virtual void	postSolveStatevar ()
	This function update internal variable after solve.
virtual void	postSolveStress ()
	This function update stress and plastic deformation gradient after solve.
virtual void	calcResidual (RankTwoTensor &)
	This function calculate stress residual.
virtual void	calcJacobian (RankFourTensor &)
	This function calculate jacobian.
virtual RankFourTensor	calcTangentModuli ()
	This function calculate the tangent moduli for preconditioner.
virtual RankFourTensor	elasticTangentModuli ()
	This function calculate the elastic tangent moduli for preconditioner.
virtual RankFourTensor	elastoPlasticTangentModuli ()
	This function calculate the exact tangent moduli for preconditioner.
RankTwoTensor	get_current_rotation (const RankTwoTensor &a)
	This function perform RU decomposition to obtain the rotation tensor.
RankTwoTensor	getMatRot (const RankTwoTensor &a)
	This function perform RU decomposition to obtain the rotation tensor.
void	calc_schmid_tensor ()
	This function calculate the Schmid tensor.
bool	line_search_update (const Real rnorm_prev, const RankTwoTensor)
	This function performs the line search update.
void	internalVariableUpdateNRiteration ()
	This function updates internal variables after each NewTon Raphson iteration (_fp_inv)
virtual void	computeQpProperties () override

Protected Attributes
DenseVector< Real >	_resid
DenseVector< Real >	_slip_rate
DenseVector< Real >	_dsliprate_dgss
DenseMatrix< Real >	_jacob
DenseMatrix< Real >	_dsliprate_dsliprate
const unsigned int	_nss
	Number of slip system resistance.
std::vector< Real >	_gprops
std::vector< Real >	_hprops
std::vector< Real >	_flowprops
std::string	_slip_sys_file_name
	File should contain slip plane normal and direction. See test.
std::string	_slip_sys_res_prop_file_name
	File should contain initial values of the slip system resistances.
std::string	_slip_sys_flow_prop_file_name
	File should contain values of the flow rate equation parameters.
std::string	_slip_sys_hard_prop_file_name
	The hardening parameters in this class are read from .i file. The user can override to read from file.
Real	_rtol
	Stress residual equation relative tolerance.
Real	_abs_tol
	Stress residual equation absolute tolerance.
Real	_gtol
	Internal variable update equation tolerance.
Real	_slip_incr_tol
	Slip increment tolerance.
unsigned int	_maxiter
	Maximum number of iterations for stress update.
unsigned int	_maxiterg
	Maximum number of iterations for internal variable update.
unsigned int	_num_slip_sys_flowrate_props
	Number of slip system flow rate parameters.
MooseEnum	_tan_mod_type
	Type of tangent moduli calculation.
MooseEnum	_intvar_read_type
	Read from options for initial values of internal variables.
unsigned int	_num_slip_sys_props
	Number of slip system specific properties provided in the file containing slip system normals and directions.
bool	_gen_rndm_stress_flag
bool	_input_rndm_scale_var
	Input option for scaling variable to generate random stress when convergence fails.
Real	_rndm_scale_var
	Scaling value.
unsigned int	_rndm_seed
	Seed value.
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
MaterialProperty< RankTwoTensor > &	_fp
const MaterialProperty< RankTwoTensor > &	_fp_old
MaterialProperty< RankTwoTensor > &	_pk2
const MaterialProperty< RankTwoTensor > &	_pk2_old
MaterialProperty< RankTwoTensor > &	_lag_e
const MaterialProperty< RankTwoTensor > &	_lag_e_old
MaterialProperty< std::vector< Real > > &	_gss
const MaterialProperty< std::vector< Real > > &	_gss_old
MaterialProperty< Real > &	_acc_slip
const MaterialProperty< Real > &	_acc_slip_old
MaterialProperty< RankTwoTensor > &	_update_rot
const MaterialProperty< RankTwoTensor > &	_deformation_gradient
const MaterialProperty< RankTwoTensor > &	_deformation_gradient_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 > &	_crysrot
DenseVector< Real >	_mo
DenseVector< Real >	_no
DenseVector< Real >	_a0
DenseVector< Real >	_xm
Real	_h0
Real	_tau_sat
Real	_tau_init
Real	_r
RankTwoTensor	_dfgrd_tmp
RankTwoTensor	_fe
RankTwoTensor	_fp_old_inv
RankTwoTensor	_fp_inv
RankTwoTensor	_fp_prev_inv
DenseVector< Real >	_slip_incr
DenseVector< Real >	_tau
DenseVector< Real >	_dslipdtau
std::vector< RankTwoTensor >	_s0
RankTwoTensor	_pk2_tmp
RankTwoTensor	_pk2_tmp_old
Real	_accslip_tmp
Real	_accslip_tmp_old
std::vector< Real >	_gss_tmp
std::vector< Real >	_gss_tmp_old
DenseVector< Real >	_slip_sys_props
DenseMatrix< Real >	_dgss_dsliprate
bool	_read_from_slip_sys_file
bool	_err_tol
RankTwoTensor	_delta_dfgrd
	Flag to check whether convergence is achieved.
RankTwoTensor	_dfgrd_tmp_old
Real	_dfgrd_scale_factor
	Scales the substepping increment to obtain deformation gradient at a substep iteration.
bool	_first_step_iter
	Flags to reset variables and reinitialize variables.
bool	_last_step_iter
bool	_first_substep
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

Definition at line 14 of file FiniteStrainCPSlipRateRes.h.

Constructor & Destructor Documentation

FiniteStrainCPSlipRateRes()

FiniteStrainCPSlipRateRes::FiniteStrainCPSlipRateRes

(

const InputParameters &

parameters

)

Definition at line 24 of file FiniteStrainCPSlipRateRes.C.

  : FiniteStrainCrystalPlasticity(parameters),
    _resid(_nss),
    _slip_rate(_nss),
    _dsliprate_dgss(_nss),
    _jacob(_nss, _nss),
    _dsliprate_dsliprate(_nss, _nss)
{
}

Member Function Documentation

assignSlipSysRes()

void FiniteStrainCrystalPlasticity::assignSlipSysRes ( )

protectedvirtualinherited

This function assign initial values of slip system resistances/internal variables read from getSlipSystems().

Definition at line 234 of file FiniteStrainCrystalPlasticity.C.

{
  _gss[_qp].resize(_nss);
 
  for (unsigned int i = 0; i < _nss; ++i)
    _gss[_qp][i] = _slip_sys_props(i);
}

Referenced by FiniteStrainCrystalPlasticity::initSlipSysProps().

calc_resid_jacob()

void FiniteStrainCrystalPlasticity::calc_resid_jacob ( RankTwoTensor &  resid, RankFourTensor &  jac  )

protectedvirtualinherited

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

Definition at line 881 of file FiniteStrainCrystalPlasticity.C.

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

Referenced by FiniteStrainCrystalPlasticity::solveStress().

calc_schmid_tensor()

void FiniteStrainCrystalPlasticity::calc_schmid_tensor ( )

protectedinherited

This function calculate the Schmid tensor.

Definition at line 1049 of file FiniteStrainCrystalPlasticity.C.

{
  DenseVector<Real> mo(LIBMESH_DIM * _nss), no(LIBMESH_DIM * _nss);
 
  // Update slip direction and normal with crystal orientation
  for (unsigned int i = 0; i < _nss; ++i)
  {
    for (const auto j : make_range(Moose::dim))
    {
      mo(i * LIBMESH_DIM + j) = 0.0;
      for (const auto k : make_range(Moose::dim))
        mo(i * LIBMESH_DIM + j) =
            mo(i * LIBMESH_DIM + j) + _crysrot[_qp](j, k) * _mo(i * LIBMESH_DIM + k);
    }
 
    for (const auto j : make_range(Moose::dim))
    {
      no(i * LIBMESH_DIM + j) = 0.0;
      for (const auto k : make_range(Moose::dim))
        no(i * LIBMESH_DIM + j) =
            no(i * LIBMESH_DIM + j) + _crysrot[_qp](j, k) * _no(i * LIBMESH_DIM + k);
    }
  }
 
  // Calculate Schmid tensor and resolved shear stresses
  for (unsigned int i = 0; i < _nss; ++i)
    for (const auto j : make_range(Moose::dim))
      for (const auto k : make_range(Moose::dim))
        _s0[i](j, k) = mo(i * LIBMESH_DIM + j) * no(i * LIBMESH_DIM + k);
}

Referenced by FiniteStrainCrystalPlasticity::preSolveQp().

calcDgssDsliprate()

void FiniteStrainCPSlipRateRes::calcDgssDsliprate ( )

protectedvirtual

This function calculates partial derivative of slip system resistances with respect to split rate.

Definition at line 223 of file FiniteStrainCPSlipRateRes.C.

{
  for (unsigned int i = 0; i < _nss; ++i)
    for (unsigned int j = 0; j < _nss; ++j)
      _dsliprate_dsliprate(i, j) += _dsliprate_dgss(i) * _dgss_dsliprate(i, j);
}

Referenced by calcJacobianSlipRate().

calcDtauDsliprate()

void FiniteStrainCPSlipRateRes::calcDtauDsliprate ( )

protectedvirtual

This function calculates partial derivative of resolved shear stress with respect to split rate.

Definition at line 188 of file FiniteStrainCPSlipRateRes.C.

{
  RankFourTensor dfedfpinv, deedfe, dfpinvdpk2;
  std::vector<RankTwoTensor> dtaudpk2(_nss), dfpinvdsliprate(_nss);
 
  for (unsigned int i = 0; i < _nss; ++i)
  {
    dtaudpk2[i] = _s0[i];
    dfpinvdsliprate[i] = -_fp_old_inv * _s0[i] * _dt;
  }
 
  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;
      }
 
  RankFourTensor dpk2dfpinv;
 
  dpk2dfpinv = _elasticity_tensor[_qp] * deedfe * dfedfpinv;
 
  for (unsigned int i = 0; i < _nss; ++i)
    for (unsigned int j = 0; j < _nss; ++j)
      _dsliprate_dsliprate(i, j) =
          _dslipdtau(i) * dtaudpk2[i].doubleContraction(dpk2dfpinv * dfpinvdsliprate[j]);
}

Referenced by calcJacobianSlipRate().

calcJacobian()

void FiniteStrainCrystalPlasticity::calcJacobian ( RankFourTensor &  jac )

protectedvirtualinherited

This function calculate jacobian.

Definition at line 927 of file FiniteStrainCrystalPlasticity.C.

{
  RankFourTensor dfedfpinv, deedfe, dfpinvdpk2;
 
  std::vector<RankTwoTensor> dtaudpk2(_nss), dfpinvdslip(_nss);
 
  for (unsigned int i = 0; i < _nss; ++i)
  {
    dtaudpk2[i] = _s0[i];
    dfpinvdslip[i] = -_fp_old_inv * _s0[i];
  }
 
  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 < _nss; ++i)
    dfpinvdpk2 += (dfpinvdslip[i] * _dslipdtau(i)).outerProduct(dtaudpk2[i]);
 
  jac =
      RankFourTensor::IdentityFour() - (_elasticity_tensor[_qp] * deedfe * dfedfpinv * dfpinvdpk2);
}

Referenced by FiniteStrainCrystalPlasticity::calc_resid_jacob().

calcJacobianSlipRate()

void FiniteStrainCPSlipRateRes::calcJacobianSlipRate ( )

protectedvirtual

This function calculates jacobian of slip rate.

Definition at line 171 of file FiniteStrainCPSlipRateRes.C.

{
  //_dsliprate_dsliprate not reinitialized to zero, hence order is important
  calcDtauDsliprate();
  calcDgssDsliprate();
 
  for (unsigned int i = 0; i < _nss; ++i)
    for (unsigned int j = 0; j < _nss; ++j)
    {
      _jacob(i, j) = 0.0;
      if (i == j)
        _jacob(i, j) += 1.0;
      _jacob(i, j) -= _dsliprate_dsliprate(i, j);
    }
}

Referenced by calcResidJacobSlipRate(), and solveStress().

calcResidDotProdUpdate()

Real FiniteStrainCPSlipRateRes::calcResidDotProdUpdate ( const DenseVector< Real > &  update )

protected

This function calculates the dot product of residual and update.

Definition at line 369 of file FiniteStrainCPSlipRateRes.C.

{
  Real dotprod = 0.0;
  for (unsigned int i = 0; i < _nss; ++i)
    dotprod += _resid(i) * update(i);
  return dotprod;
}

Referenced by lineSearchUpdateSlipRate().

calcResidJacobSlipRate()

void FiniteStrainCPSlipRateRes::calcResidJacobSlipRate ( )

protectedvirtual

This function calculates residual and jacobian of slip rate.

Definition at line 126 of file FiniteStrainCPSlipRateRes.C.

{
  calcResidualSlipRate();
  if (_err_tol)
    return;
  calcJacobianSlipRate();
}

Referenced by solveStress().

calcResidNorm()

Real FiniteStrainCPSlipRateRes::calcResidNorm ( )

protectedvirtual

This function calculates the residual norm.

Definition at line 261 of file FiniteStrainCPSlipRateRes.C.

{
  Real rnorm = 0.0;
  for (unsigned int i = 0; i < _nss; ++i)
    rnorm += Utility::pow<2>(_resid(i));
  rnorm = std::sqrt(rnorm) / _nss;
 
  return rnorm;
}

Referenced by lineSearchUpdateSlipRate(), and solveStress().

calcResidual()

void FiniteStrainCrystalPlasticity::calcResidual ( RankTwoTensor &  resid )

protectedvirtualinherited

This function calculate stress residual.

Definition at line 890 of file FiniteStrainCrystalPlasticity.C.

{
  RankTwoTensor iden(RankTwoTensor::initIdentity), ce, ee, ce_pk2, eqv_slip_incr, pk2_new;
 
  _fe = _dfgrd_tmp * _fp_prev_inv; // _fp_inv  ==> _fp_prev_inv
 
  ce = _fe.transpose() * _fe;
  ce_pk2 = ce * _pk2_tmp;
  ce_pk2 = ce_pk2 / _fe.det();
 
  // Calculate Schmid tensor and resolved shear stresses
  for (unsigned int i = 0; i < _nss; ++i)
    _tau(i) = ce_pk2.doubleContraction(_s0[i]);
 
  getSlipIncrements(); // Calculate dslip,dslipdtau
 
  if (_err_tol)
    return;
 
  eqv_slip_incr.zero();
  for (unsigned int i = 0; i < _nss; ++i)
    eqv_slip_incr += _s0[i] * _slip_incr(i);
 
  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_tmp - pk2_new;
}

Referenced by FiniteStrainCrystalPlasticity::calc_resid_jacob(), and FiniteStrainCrystalPlasticity::line_search_update().

calcResidualSlipRate()

void FiniteStrainCPSlipRateRes::calcResidualSlipRate ( )

protectedvirtual

This function calculates residual of slip rate.

Definition at line 135 of file FiniteStrainCPSlipRateRes.C.

{
  RankTwoTensor eqv_slip_incr, ce, ee;
  const RankTwoTensor iden(RankTwoTensor::initIdentity);
 
  _slip_incr = _slip_rate;
  _slip_incr *= _dt;
 
  for (unsigned int i = 0; i < _nss; ++i)
    eqv_slip_incr += _s0[i] * _slip_incr(i);
 
  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_tmp = _elasticity_tensor[_qp] * ee;
 
  for (unsigned int i = 0; i < _nss; ++i)
    _tau(i) = _pk2_tmp.doubleContraction(_s0[i]);
 
  update_slip_system_resistance();
  getSlipIncrements();
 
  if (_err_tol)
    return;
 
  for (unsigned int i = 0; i < _nss; ++i)
    _resid(i) = _slip_rate(i) - _slip_incr(i) / _dt;
}

Referenced by calcResidJacobSlipRate(), lineSearchUpdateSlipRate(), and solveStress().

calcTangentModuli()

RankFourTensor FiniteStrainCrystalPlasticity::calcTangentModuli ( )

protectedvirtualinherited

This function 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 1032 of file FiniteStrainCrystalPlasticity.C.

{
  RankFourTensor tan_mod;
 
  switch (_tan_mod_type)
  {
    case 0:
      tan_mod = elastoPlasticTangentModuli();
      break;
    default:
      tan_mod = elasticTangentModuli();
  }
 
  return tan_mod;
}

Referenced by FiniteStrainCrystalPlasticity::postSolveQp().

calcUpdate()

void FiniteStrainCPSlipRateRes::calcUpdate ( )

protected

This function calculates and updates the residual of slip rate.

Definition at line 247 of file FiniteStrainCPSlipRateRes.C.

{
  DenseMatrix<Real> A = _jacob;
  DenseVector<Real> r(_nss);
  DenseVector<Real> x(_nss);
 
  r = _resid;
 
  A.lu_solve(r, x);
 
  _resid = x;
}

Referenced by solveStress().

computeQpElasticityTensor()

void FiniteStrainCrystalPlasticity::computeQpElasticityTensor ( )

protectedvirtualinherited

This function updates the elasticity tensor at a quadrature point.

Presently void.

Definition at line 1027 of file FiniteStrainCrystalPlasticity.C.

{
}

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 FiniteStrainCrystalPlasticity::computeQpStress ( )

protectedvirtualinherited

This function updates the stress at a quadrature point.

Solves stress residual equation using NR.

Updates slip system resistances iteratively.

Implements ComputeGeneralStressBase.

Definition at line 493 of file FiniteStrainCrystalPlasticity.C.

{
  unsigned int substep_iter = 1; // Depth of substepping; Limited to maximum substep iteration
  unsigned int num_substep = 1;  // Calculated from substep_iter as 2^substep_iter
  Real dt_original = _dt;        // Stores original _dt; Reset at the end of solve
  _first_substep = true;         // Initialize variables at substep_iter = 1
 
  if (_max_substep_iter > 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;
    _err_tol = true; // Indicator to continue substepping
  }
 
  // Substepping loop
  while (_err_tol && _max_substep_iter > 1)
  {
    _dt = dt_original / num_substep;
 
    for (unsigned int istep = 0; istep < num_substep; ++istep)
    {
      _first_step_iter = false;
      if (istep == 0)
        _first_step_iter = true;
 
      _last_step_iter = false;
      if (istep == num_substep - 1)
        _last_step_iter = true;
 
      _dfgrd_scale_factor = (static_cast<Real>(istep) + 1) / num_substep;
 
      preSolveQp();
      solveQp();
 
      if (_err_tol)
      {
        substep_iter++;
        num_substep *= 2;
        break;
      }
    }
 
    _first_substep = false; // Prevents reinitialization
    _dt = dt_original;      // Resets dt
 
#ifdef DEBUG
    if (substep_iter > _max_substep_iter)
      mooseWarning("FiniteStrainCrystalPlasticity: Failure with substepping");
#endif
 
    if (!_err_tol || substep_iter > _max_substep_iter)
      postSolveQp(); // Evaluate variables after successful solve or indicate failure
  }
 
  // No substepping
  if (_max_substep_iter == 1)
  {
    preSolveQp();
    solveQp();
    postSolveQp();
  }
}

elasticTangentModuli()

RankFourTensor FiniteStrainCrystalPlasticity::elasticTangentModuli ( )

protectedvirtualinherited

This function calculate the elastic tangent moduli for preconditioner.

Definition at line 1121 of file FiniteStrainCrystalPlasticity.C.

{
  return _elasticity_tensor[_qp]; // update jacobian_mult
}

Referenced by FiniteStrainCrystalPlasticity::calcTangentModuli().

elastoPlasticTangentModuli()

RankFourTensor FiniteStrainCrystalPlasticity::elastoPlasticTangentModuli ( )

protectedvirtualinherited

This function calculate the exact tangent moduli for preconditioner.

Definition at line 1081 of file FiniteStrainCrystalPlasticity.C.

{
  RankFourTensor tan_mod;
  RankTwoTensor pk2fet, fepk2;
  RankFourTensor deedfe, dsigdpk2dfe;
 
  // 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_tmp * _fe.transpose();
  fepk2 = _fe * _pk2_tmp;
 
  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) = tan_mod(i, j, i, l) + pk2fet(l, j);
        tan_mod(i, j, j, l) = tan_mod(i, j, j, l) + fepk2(i, l);
      }
 
  tan_mod += dsigdpk2dfe;
 
  Real je = _fe.det();
  if (je > 0.0)
    tan_mod /= je;
 
  return tan_mod;
}

Referenced by FiniteStrainCrystalPlasticity::calcTangentModuli().

get_current_rotation()

RankTwoTensor FiniteStrainCrystalPlasticity::get_current_rotation ( const RankTwoTensor &  a )

protectedinherited

This function perform RU decomposition to obtain the rotation tensor.

Definition at line 985 of file FiniteStrainCrystalPlasticity.C.

{
  return getMatRot(a);
}

Referenced by FiniteStrainCrystalPlasticity::postSolveQp().

getFlowRateParams()

void FiniteStrainCrystalPlasticity::getFlowRateParams ( )

protectedvirtualinherited

This function assign flow rate parameters - see test.

.i input file format start_slip_sys_num, end_slip_sys_num, value1, value2

Definition at line 344 of file FiniteStrainCrystalPlasticity.C.

{
  if (_flowprops.size() <= 0)
    mooseError("FiniteStrainCrystalPLasticity: Error in reading flow rate  properties: Specify "
               "input in .i file or a slip_sys_flow_prop_file_name");
 
  _a0.resize(_nss);
  _xm.resize(_nss);
 
  unsigned int num_data_grp = 2 + _num_slip_sys_flowrate_props; // Number of data per group e.g.
                                                                // start_slip_sys, end_slip_sys,
                                                                // value1, value2, ..
 
  for (unsigned int i = 0; i < _flowprops.size() / num_data_grp; ++i)
  {
    Real vs, ve;
    unsigned int is, ie;
 
    vs = _flowprops[i * num_data_grp];
    ve = _flowprops[i * num_data_grp + 1];
 
    if (vs <= 0 || ve <= 0)
      mooseError("FiniteStrainCrystalPLasticity: Indices in flow rate parameter read must be "
                 "positive integers: is = ",
                 vs,
                 " ie = ",
                 ve);
 
    if (vs != floor(vs) || ve != floor(ve))
      mooseError("FiniteStrainCrystalPLasticity: Error in reading flow props: Values specifying "
                 "start and end number of slip system groups should be integer");
 
    is = static_cast<unsigned int>(vs);
    ie = static_cast<unsigned int>(ve);
 
    if (is > ie)
      mooseError("FiniteStrainCrystalPLasticity: Start index is = ",
                 is,
                 " should be greater than end index ie = ",
                 ie,
                 " in flow rate parameter read");
 
    for (unsigned int j = is; j <= ie; ++j)
    {
      _a0(j - 1) = _flowprops[i * num_data_grp + 2];
      _xm(j - 1) = _flowprops[i * num_data_grp + 3];
    }
  }
 
  for (unsigned int i = 0; i < _nss; ++i)
  {
    if (!(_a0(i) > 0.0 && _xm(i) > 0.0))
    {
      mooseWarning(
          "FiniteStrainCrystalPlasticity: Non-positive flow rate parameters ", _a0(i), ",", _xm(i));
      break;
    }
  }
}

Referenced by FiniteStrainCrystalPlasticity::initSlipSysProps().

getHardnessParams()

void FiniteStrainCrystalPlasticity::getHardnessParams ( )

protectedvirtualinherited

This function assign flow rate parameters from .i file - see test.

Definition at line 412 of file FiniteStrainCrystalPlasticity.C.

{
  if (_hprops.size() <= 0)
    mooseError("FiniteStrainCrystalPLasticity: Error in reading hardness properties: Specify input "
               "in .i file or a slip_sys_hard_prop_file_name");
 
  _r = _hprops[0];
  _h0 = _hprops[1];
  _tau_init = _hprops[2];
  _tau_sat = _hprops[3];
}

Referenced by FiniteStrainCrystalPlasticity::initSlipSysProps().

getInitSlipSysRes()

void FiniteStrainCrystalPlasticity::getInitSlipSysRes ( )

protectedvirtualinherited

This function assign slip system resistances - see test.

.i input file format start_slip_sys_num, end_slip_sys_num, value.

Definition at line 262 of file FiniteStrainCrystalPlasticity.C.

{
  if (_gprops.size() <= 0)
    mooseError("FiniteStrainCrystalPLasticity: Error in reading slip system resistance properties: "
               "Specify input in .i file or in slip_sys_res_prop_file or in slip_sys_file");
 
  _gss[_qp].resize(_nss, 0.0);
 
  unsigned int num_data_grp = 3; // Number of data per group e.g. start_slip_sys, end_slip_sys,
                                 // value
 
  for (unsigned int i = 0; i < _gprops.size() / num_data_grp; ++i)
  {
    Real vs, ve;
    unsigned int is, ie;
 
    vs = _gprops[i * num_data_grp];
    ve = _gprops[i * num_data_grp + 1];
 
    if (vs <= 0 || ve <= 0)
      mooseError("FiniteStrainCrystalPLasticity: Indices in gss property read must be positive "
                 "integers: is = ",
                 vs,
                 " ie = ",
                 ve);
 
    if (vs != floor(vs) || ve != floor(ve))
      mooseError("FiniteStrainCrystalPLasticity: Error in reading slip system resistances: Values "
                 "specifying start and end number of slip system groups should be integer");
 
    is = static_cast<unsigned int>(vs);
    ie = static_cast<unsigned int>(ve);
 
    if (is > ie)
      mooseError("FiniteStrainCrystalPLasticity: Start index is = ",
                 is,
                 " should be greater than end index ie = ",
                 ie,
                 " in slip system resistance property read");
 
    for (unsigned int j = is; j <= ie; ++j)
      _gss[_qp][j - 1] = _gprops[i * num_data_grp + 2];
  }
 
  for (unsigned int i = 0; i < _nss; ++i)
    if (_gss[_qp][i] <= 0.0)
      mooseError("FiniteStrainCrystalPLasticity: Value of resistance for slip system ",
                 i + 1,
                 " non positive");
}

Referenced by FiniteStrainCrystalPlasticity::initSlipSysProps().

getMatRot()

RankTwoTensor FiniteStrainCrystalPlasticity::getMatRot ( const RankTwoTensor &  a )

protectedinherited

This function perform RU decomposition to obtain the rotation tensor.

Definition at line 992 of file FiniteStrainCrystalPlasticity.C.

{
  RankTwoTensor rot;
  RankTwoTensor c, diag, evec;
  PetscScalar cmat[LIBMESH_DIM][LIBMESH_DIM], work[10];
  PetscReal w[LIBMESH_DIM];
  PetscBLASInt nd = LIBMESH_DIM, lwork = 10, info;
 
  c = a.transpose() * a;
 
  for (const auto i : make_range(Moose::dim))
    for (const auto j : make_range(Moose::dim))
      cmat[i][j] = c(i, j);
 
  LAPACKsyev_("V", "U", &nd, &cmat[0][0], &nd, w, work, &lwork, &info);
 
  if (info != 0)
    mooseError("FiniteStrainCrystalPLasticity: DSYEV function call in getMatRot function failed");
 
  diag.zero();
 
  for (const auto i : make_range(Moose::dim))
    diag(i, i) = std::sqrt(w[i]);
 
  for (const auto i : make_range(Moose::dim))
    for (const auto j : make_range(Moose::dim))
      evec(i, j) = cmat[i][j];
 
  rot = a * ((evec.transpose() * diag * evec).inverse());
 
  return rot;
}

Referenced by FiniteStrainCrystalPlasticity::get_current_rotation().

getSlipIncrements()

void FiniteStrainCPSlipRateRes::getSlipIncrements ( )

protectedvirtual

This function updates the slip system resistances.

Reimplemented from FiniteStrainCrystalPlasticity.

Definition at line 231 of file FiniteStrainCPSlipRateRes.C.

{
  FiniteStrainCrystalPlasticity::getSlipIncrements();
 
  if (_err_tol)
    return;
 
  _dslipdtau *= 1.0 / _dt;
 
  for (unsigned int i = 0; i < _nss; ++i)
    _dsliprate_dgss(i) = -_a0(i) / _xm(i) *
                         std::pow(std::abs(_tau(i) / _gss_tmp[i]), 1.0 / _xm(i) - 1.0) * _tau(i) /
                         std::pow(_gss_tmp[i], 2.0);
}

Referenced by calcResidualSlipRate().

getSlipSystems()

void FiniteStrainCrystalPlasticity::getSlipSystems ( )

protectedvirtualinherited

This function reads slip system from file - see test.

Definition at line 426 of file FiniteStrainCrystalPlasticity.C.

{
  Real vec[LIBMESH_DIM];
  std::ifstream fileslipsys;
 
  MooseUtils::checkFileReadable(_slip_sys_file_name);
 
  fileslipsys.open(_slip_sys_file_name.c_str());
 
  for (unsigned int i = 0; i < _nss; ++i)
  {
    // Read the slip normal
    for (const auto j : make_range(Moose::dim))
      if (!(fileslipsys >> vec[j]))
        mooseError("Crystal Plasticity Error: Premature end of file reading slip system file \n");
 
    // Normalize the vectors
    Real mag;
    mag = Utility::pow<2>(vec[0]) + Utility::pow<2>(vec[1]) + Utility::pow<2>(vec[2]);
    mag = std::sqrt(mag);
 
    for (unsigned j = 0; j < LIBMESH_DIM; ++j)
      _no(i * LIBMESH_DIM + j) = vec[j] / mag;
 
    // Read the slip direction
    for (const auto j : make_range(Moose::dim))
      if (!(fileslipsys >> vec[j]))
        mooseError("Crystal Plasticity Error: Premature end of file reading slip system file \n");
 
    // Normalize the vectors
    mag = Utility::pow<2>(vec[0]) + Utility::pow<2>(vec[1]) + Utility::pow<2>(vec[2]);
    mag = std::sqrt(mag);
 
    for (const auto j : make_range(Moose::dim))
      _mo(i * LIBMESH_DIM + j) = vec[j] / mag;
 
    mag = 0.0;
    for (const auto j : make_range(Moose::dim))
      mag += _mo(i * LIBMESH_DIM + j) * _no(i * LIBMESH_DIM + j);
 
    if (std::abs(mag) > 1e-8)
      mooseError(
          "Crystal Plasicity Error: Slip direction and normal not orthonormal, System number = ",
          i,
          "\n");
 
    if (_read_from_slip_sys_file)
      for (unsigned int j = 0; j < _num_slip_sys_props; ++j)
        if (!(fileslipsys >> _slip_sys_props(i * _num_slip_sys_props + j)))
          mooseError("Crystal Plasticity Error: Premature end of file reading slip system file - "
                     "check in slip system file read input options/values\n");
  }
 
  fileslipsys.close();
}

Referenced by FiniteStrainCrystalPlasticity::FiniteStrainCrystalPlasticity().

initAdditionalProps()

void FiniteStrainCrystalPlasticity::initAdditionalProps ( )

protectedvirtualinherited

This function initializes additional parameters.

Definition at line 484 of file FiniteStrainCrystalPlasticity.C.

{
}

Referenced by FiniteStrainCrystalPlasticity::initQpStatefulProperties().

initQpStatefulProperties()

void FiniteStrainCrystalPlasticity::initQpStatefulProperties ( )

protectedvirtualinherited

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

Reimplemented from ComputeGeneralStressBase.

Definition at line 191 of file FiniteStrainCrystalPlasticity.C.

{
  _stress[_qp].zero();
 
  _fp[_qp].setToIdentity();
 
  _pk2[_qp].zero();
  _acc_slip[_qp] = 0.0;
  _lag_e[_qp].zero();
 
  _update_rot[_qp].setToIdentity();
 
  initSlipSysProps(); // Initializes slip system related properties
  initAdditionalProps();
}

initSlipSysProps()

void FiniteStrainCrystalPlasticity::initSlipSysProps ( )

protectedvirtualinherited

This function initializes slip system resistances.

Definition at line 208 of file FiniteStrainCrystalPlasticity.C.

{
  switch (_intvar_read_type)
  {
    case 0:
      assignSlipSysRes();
      break;
    case 1:
      readFileInitSlipSysRes();
      break;
    default:
      getInitSlipSysRes();
  }
 
  if (_slip_sys_flow_prop_file_name.length() != 0)
    readFileFlowRateParams();
  else
    getFlowRateParams();
 
  if (_slip_sys_hard_prop_file_name.length() != 0)
    readFileHardnessParams();
  else
    getHardnessParams();
}

Referenced by FiniteStrainCrystalPlasticity::initQpStatefulProperties().

internalVariableUpdateNRiteration()

void FiniteStrainCrystalPlasticity::internalVariableUpdateNRiteration ( )

protectedinherited

This function updates internal variables after each NewTon Raphson iteration (_fp_inv)

Definition at line 1210 of file FiniteStrainCrystalPlasticity.C.

{
  _fp_prev_inv = _fp_inv; // update _fp_prev_inv
}

Referenced by FiniteStrainCrystalPlasticity::solveStress().

line_search_update()

bool FiniteStrainCrystalPlasticity::line_search_update ( const Real  rnorm_prev, const RankTwoTensor  dpk2  )

protectedinherited

This function performs the line search update.

Definition at line 1127 of file FiniteStrainCrystalPlasticity.C.

{
  if (_lsrch_method == "CUT_HALF")
  {
    Real rnorm;
    RankTwoTensor resid;
    Real step = 1.0;
 
    do
    {
      _pk2_tmp = _pk2_tmp - step * dpk2;
      step /= 2.0;
      _pk2_tmp = _pk2_tmp + step * dpk2;
 
      calcResidual(resid);
      rnorm = resid.L2norm();
    } while (rnorm > rnorm_prev && step > _min_lsrch_step);
 
    if (rnorm > rnorm_prev && step <= _min_lsrch_step)
      return false;
 
    return true;
  }
  else if (_lsrch_method == "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;
 
    RankTwoTensor resid;
    calcResidual(resid);
    Real s_b = resid.doubleContraction(dpk2);
    Real rnorm1 = resid.L2norm();
    _pk2_tmp = _pk2_tmp - dpk2;
    calcResidual(resid);
    Real s_a = resid.doubleContraction(dpk2);
    Real rnorm0 = resid.L2norm();
    _pk2_tmp = _pk2_tmp + dpk2;
 
    if ((rnorm1 / rnorm0) < _lsrch_tol || s_a * s_b > 0)
    {
      calcResidual(resid);
      return true;
    }
 
    while ((rnorm / rnorm0) > _lsrch_tol && count < _lsrch_max_iter)
    {
      _pk2_tmp = _pk2_tmp - step * dpk2;
      step = 0.5 * (step_b + step_a);
      _pk2_tmp = _pk2_tmp + step * dpk2;
      calcResidual(resid);
      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++;
    }
 
    if ((rnorm / rnorm0) < _lsrch_tol && count < _lsrch_max_iter)
      return true;
 
    return false;
  }
  else
  {
    mooseError("Line search meothod is not provided.");
    return false;
  }
}

Referenced by FiniteStrainCrystalPlasticity::solveStress().

lineSearchUpdateSlipRate()

bool FiniteStrainCPSlipRateRes::lineSearchUpdateSlipRate ( const Real  rnorm_prev, const DenseVector< Real > &  update  )

protected

This function performs the line search update.

Definition at line 272 of file FiniteStrainCPSlipRateRes.C.

{
  if (_lsrch_method == "CUT_HALF")
  {
    Real rnorm;
    Real step = 1.0;
    do
    {
      for (unsigned int i = 0; i < update.size(); ++i)
        _slip_rate(i) += step * update(i);
 
      step /= 2.0;
 
      for (unsigned int i = 0; i < update.size(); ++i)
        _slip_rate(i) -= step * update(i);
 
      calcResidualSlipRate();
      if (_err_tol)
        return false;
      rnorm = calcResidNorm();
    } while (rnorm > rnorm_prev && step > _min_lsrch_step);
 
    if (rnorm > rnorm_prev && step <= _min_lsrch_step)
      return false;
 
    return true;
  }
  else if (_lsrch_method == "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;
 
    Real s_b = calcResidDotProdUpdate(update);
    Real rnorm1 = calcResidNorm();
 
    for (unsigned int i = 0; i < update.size(); ++i)
      _slip_rate(i) += update(i);
 
    calcResidualSlipRate();
    Real s_a = calcResidDotProdUpdate(update);
    Real rnorm0 = calcResidNorm();
 
    for (unsigned int i = 0; i < update.size(); ++i)
      _slip_rate(i) -= update(i);
 
    if ((rnorm1 / rnorm0) < _lsrch_tol || s_a * s_b > 0)
    {
      calcResidualSlipRate();
      return true;
    }
 
    while ((rnorm / rnorm0) > _lsrch_tol && count < _lsrch_max_iter)
    {
 
      for (unsigned int i = 0; i < update.size(); ++i)
        _slip_rate(i) += step * update(i);
 
      step = 0.5 * (step_a + step_b);
 
      for (unsigned int i = 0; i < update.size(); ++i)
        _slip_rate(i) -= step * update(i);
 
      calcResidualSlipRate();
      s_m = calcResidDotProdUpdate(update);
      rnorm = calcResidNorm();
 
      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++;
    }
 
    if ((rnorm / rnorm0) < _lsrch_tol && count < _lsrch_max_iter)
      return true;
 
    return false;
  }
  else
  {
    mooseError("Line search meothod is not provided.");
    return false;
  }
}

Referenced by solveStress().

postSolveQp()

void FiniteStrainCrystalPlasticity::postSolveQp ( )

protectedvirtualinherited

This function update stress and internal variable after solve.

Definition at line 588 of file FiniteStrainCrystalPlasticity.C.

{
  if (_err_tol)
  {
    _err_tol = false;
    if (_gen_rndm_stress_flag)
    {
      if (!_input_rndm_scale_var)
        _rndm_scale_var = _elasticity_tensor[_qp](0, 0, 0, 0);
 
      _stress[_qp] = RankTwoTensor::genRandomSymmTensor(_rndm_scale_var, 1.0);
    }
    else
      mooseError("FiniteStrainCrystalPlasticity: Constitutive failure");
  }
  else
  {
    _stress[_qp] = _fe * _pk2[_qp] * _fe.transpose() / _fe.det();
 
    _Jacobian_mult[_qp] += calcTangentModuli(); // Calculate jacobian for preconditioner
 
    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;
    rot = get_current_rotation(_deformation_gradient[_qp]); // Calculate material rotation
    _update_rot[_qp] = rot * _crysrot[_qp];
  }
}

Referenced by FiniteStrainCrystalPlasticity::computeQpStress().

postSolveStatevar()

void FiniteStrainCrystalPlasticity::postSolveStatevar ( )

protectedvirtualinherited

This function update internal variable after solve.

Definition at line 683 of file FiniteStrainCrystalPlasticity.C.

{
  if (_max_substep_iter == 1) // No substepping
  {
    _gss[_qp] = _gss_tmp;
    _acc_slip[_qp] = _accslip_tmp;
  }
  else
  {
    if (_last_step_iter)
    {
      _gss[_qp] = _gss_tmp;
      _acc_slip[_qp] = _accslip_tmp;
    }
    else
    {
      _gss_tmp_old = _gss_tmp;
      _accslip_tmp_old = _accslip_tmp;
    }
  }
}

Referenced by FiniteStrainCrystalPlasticity::solveQp().

postSolveStress()

void FiniteStrainCrystalPlasticity::postSolveStress ( )

protectedvirtualinherited

This function update stress and plastic deformation gradient after solve.

Definition at line 802 of file FiniteStrainCrystalPlasticity.C.

{
  if (_max_substep_iter == 1) // No substepping
  {
    _fp[_qp] = _fp_inv.inverse();
    _pk2[_qp] = _pk2_tmp;
  }
  else
  {
    if (_last_step_iter)
    {
      _fp[_qp] = _fp_inv.inverse();
      _pk2[_qp] = _pk2_tmp;
    }
    else
    {
      _fp_old_inv = _fp_inv;
      _pk2_tmp_old = _pk2_tmp;
    }
  }
}

Referenced by solveStatevar(), and FiniteStrainCrystalPlasticity::solveStatevar().

preSolveQp()

void FiniteStrainCrystalPlasticity::preSolveQp ( )

protectedvirtualinherited

This function set variables for stress and internal variable solve.

Definition at line 560 of file FiniteStrainCrystalPlasticity.C.

{
  // Initialize variable
  if (_first_substep)
  {
    _Jacobian_mult[_qp].zero(); // Initializes jacobian for preconditioner
    calc_schmid_tensor();
  }
 
  if (_max_substep_iter == 1)
    _dfgrd_tmp = _deformation_gradient[_qp]; // Without substepping
  else
    _dfgrd_tmp = _dfgrd_scale_factor * _delta_dfgrd + _dfgrd_tmp_old;
 
  _err_tol = false;
}

Referenced by FiniteStrainCrystalPlasticity::computeQpStress().

preSolveStatevar()

void FiniteStrainCrystalPlasticity::preSolveStatevar ( )

protectedvirtualinherited

This function set variables for internal variable solve.

Definition at line 621 of file FiniteStrainCrystalPlasticity.C.

{
  if (_max_substep_iter == 1) // No substepping
  {
    _gss_tmp = _gss_old[_qp];
    _accslip_tmp_old = _acc_slip_old[_qp];
  }
  else
  {
    if (_first_step_iter)
    {
      _gss_tmp = _gss_tmp_old = _gss_old[_qp];
      _accslip_tmp_old = _acc_slip_old[_qp];
    }
    else
      _gss_tmp = _gss_tmp_old;
  }
}

Referenced by FiniteStrainCrystalPlasticity::solveQp().

preSolveStress()

void FiniteStrainCPSlipRateRes::preSolveStress ( )

protectedvirtual

This function sets variable for internal variable solve.

Reimplemented from FiniteStrainCrystalPlasticity.

Definition at line 45 of file FiniteStrainCPSlipRateRes.C.

{
  FiniteStrainCrystalPlasticity::preSolveStress();
  _slip_rate.zero();
}

Referenced by solveStatevar().

readFileFlowRateParams()

void FiniteStrainCrystalPlasticity::readFileFlowRateParams ( )

protectedvirtualinherited

This function read flow rate parameters from file - see test.

Definition at line 315 of file FiniteStrainCrystalPlasticity.C.

{
  _a0.resize(_nss);
  _xm.resize(_nss);
 
  MooseUtils::checkFileReadable(_slip_sys_flow_prop_file_name);
 
  std::ifstream file;
  file.open(_slip_sys_flow_prop_file_name.c_str());
 
  std::vector<Real> vec;
  vec.resize(_num_slip_sys_flowrate_props);
 
  for (unsigned int i = 0; i < _nss; ++i)
  {
    for (unsigned int j = 0; j < _num_slip_sys_flowrate_props; ++j)
      if (!(file >> vec[j]))
        mooseError(
            "Error FiniteStrainCrystalPlasticity: Premature end of slip_sys_flow_rate_param file");
 
    _a0(i) = vec[0];
    _xm(i) = vec[1];
  }
 
  file.close();
}

Referenced by FiniteStrainCrystalPlasticity::initSlipSysProps().

readFileHardnessParams()

void FiniteStrainCrystalPlasticity::readFileHardnessParams ( )

protectedvirtualinherited

This function read hardness parameters from file.

Definition at line 406 of file FiniteStrainCrystalPlasticity.C.

{
}

Referenced by FiniteStrainCrystalPlasticity::initSlipSysProps().

readFileInitSlipSysRes()

void FiniteStrainCrystalPlasticity::readFileInitSlipSysRes ( )

protectedvirtualinherited

This function read slip system resistances from file - see test.

Definition at line 244 of file FiniteStrainCrystalPlasticity.C.

{
  _gss[_qp].resize(_nss);
 
  MooseUtils::checkFileReadable(_slip_sys_res_prop_file_name);
 
  std::ifstream file;
  file.open(_slip_sys_res_prop_file_name.c_str());
 
  for (unsigned int i = 0; i < _nss; ++i)
    if (!(file >> _gss[_qp][i]))
      mooseError("Error FiniteStrainCrystalPlasticity: Premature end of slip_sys_res_prop file");
 
  file.close();
}

Referenced by FiniteStrainCrystalPlasticity::initSlipSysProps().

solveQp()

void FiniteStrainCrystalPlasticity::solveQp ( )

protectedvirtualinherited

This function solves stress and internal variables.

Definition at line 578 of file FiniteStrainCrystalPlasticity.C.

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

Referenced by FiniteStrainCrystalPlasticity::computeQpStress().

solveStatevar()

void FiniteStrainCPSlipRateRes::solveStatevar ( )

protectedvirtual

This function solves internal variables.

Reimplemented from FiniteStrainCrystalPlasticity.

Definition at line 35 of file FiniteStrainCPSlipRateRes.C.

{
  preSolveStress();
  solveStress();
  if (_err_tol)
    return;
  postSolveStress();
}

solveStress()

void FiniteStrainCPSlipRateRes::solveStress ( )

protectedvirtual

This function solves for stress, updates plastic deformation gradient.

Reimplemented from FiniteStrainCrystalPlasticity.

Definition at line 52 of file FiniteStrainCPSlipRateRes.C.

{
  Real rnorm, rnorm0, rnorm_prev;
  unsigned int iter = 0;
 
#ifdef DEBUG
  std::vector<Real> rnormst(_maxiter + 1), slipratest(_maxiter + 1); // Use for Debugging
#endif
 
  calcResidJacobSlipRate();
  if (_err_tol)
    return;
  rnorm = calcResidNorm();
  rnorm0 = rnorm;
 
#ifdef DEBUG
  rnormst[iter] = rnorm;
  Real slipratemax = 0.0;
  for (unsigned int i = 0; i < _nss; ++i)
    if (std::abs(_slip_rate(i)) > slipratemax)
      slipratemax = std::abs(_slip_rate(i));
  slipratest[iter] = slipratemax;
#endif
 
  while (rnorm > _rtol * rnorm0 && rnorm0 > _abs_tol && iter < _maxiter)
  {
    calcUpdate();
 
    DenseVector<Real> update = _resid;
 
    _slip_rate -= update;
 
    calcResidualSlipRate();
    if (_err_tol)
      return;
    rnorm_prev = rnorm;
    rnorm = calcResidNorm();
 
    if (_use_line_search && rnorm > rnorm_prev && !lineSearchUpdateSlipRate(rnorm_prev, update))
    {
#ifdef DEBUG
      mooseWarning("FiniteStrainCrystalPLasticity: Failed with line search");
#endif
      _err_tol = true;
      return;
    }
 
    calcJacobianSlipRate();
 
    if (_use_line_search)
      rnorm = calcResidNorm();
    iter++;
 
#ifdef DEBUG
    slipratemax = 0.0;
    for (unsigned int i = 0; i < _nss; ++i)
      if (std::abs(_slip_rate(i)) > slipratemax)
        slipratemax = std::abs(_slip_rate(i));
    rnormst[iter] = rnorm;
    slipratest[iter] = slipratemax;
#endif
  }
 
  if (iter == _maxiter)
  {
#ifdef DEBUG
    mooseWarning("FiniteStrainCPSlipRateRes: NR exceeds maximum iteration ", iter, " ", rnorm);
#endif
    _err_tol = true;
    return;
  }
}

Referenced by solveStatevar().

update_slip_system_resistance()

void FiniteStrainCrystalPlasticity::update_slip_system_resistance ( )

protectedvirtualinherited

This function updates the slip system resistances.

Definition at line 826 of file FiniteStrainCrystalPlasticity.C.

{
  updateGss();
}

Referenced by calcResidualSlipRate(), and FiniteStrainCrystalPlasticity::solveStatevar().

updateGss()

void FiniteStrainCrystalPlasticity::updateGss ( )

protectedvirtualinherited

This function updates the slip system resistances.

Old function to update slip system resistances.

Kept to avoid code break at computeQpstress

Definition at line 836 of file FiniteStrainCrystalPlasticity.C.

{
  DenseVector<Real> hb(_nss);
  Real qab;
 
  Real a = _hprops[4]; // Kalidindi
 
  _accslip_tmp = _accslip_tmp_old;
  for (unsigned int i = 0; i < _nss; ++i)
    _accslip_tmp += std::abs(_slip_incr(i));
 
  // Real val = std::cosh(_h0 * _accslip_tmp / (_tau_sat - _tau_init)); // Karthik
  // val = _h0 * std::pow(1.0/val,2.0); // Kalidindi
 
  for (unsigned int i = 0; i < _nss; ++i)
    // hb(i)=val;
    hb(i) = _h0 * std::pow(std::abs(1.0 - _gss_tmp[i] / _tau_sat), a) *
            std::copysign(1.0, 1.0 - _gss_tmp[i] / _tau_sat);
 
  for (unsigned int i = 0; i < _nss; ++i)
  {
    if (_max_substep_iter == 1) // No substepping
      _gss_tmp[i] = _gss_old[_qp][i];
    else
      _gss_tmp[i] = _gss_tmp_old[i];
 
    for (unsigned int j = 0; j < _nss; ++j)
    {
      unsigned int iplane, jplane;
      iplane = i / 3;
      jplane = j / 3;
 
      if (iplane == jplane) // Kalidindi
        qab = 1.0;
      else
        qab = _r;
 
      _gss_tmp[i] += qab * hb(j) * std::abs(_slip_incr(j));
      _dgss_dsliprate(i, j) = qab * hb(j) * std::copysign(1.0, _slip_incr(j)) * _dt;
    }
  }
}

Referenced by FiniteStrainCrystalPlasticity::update_slip_system_resistance().

validParams()

InputParameters FiniteStrainCPSlipRateRes::validParams ( )

static

Definition at line 16 of file FiniteStrainCPSlipRateRes.C.

{
  InputParameters params = FiniteStrainCrystalPlasticity::validParams();
  params.addClassDescription("Deprecated class: please use CrystalPlasticityKalidindiUpdate and "
                             "ComputeMultipleCrystalPlasticityStress instead.");
  return params;
}

Member Data Documentation

_a0

DenseVector<Real> FiniteStrainCrystalPlasticity::_a0

protectedinherited

Definition at line 322 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::getFlowRateParams(), getSlipIncrements(), FiniteStrainCrystalPlasticity::getSlipIncrements(), and FiniteStrainCrystalPlasticity::readFileFlowRateParams().

_abs_tol

Real FiniteStrainCrystalPlasticity::_abs_tol

protectedinherited

Stress residual equation absolute tolerance.

Definition at line 247 of file FiniteStrainCrystalPlasticity.h.

Referenced by solveStress(), and FiniteStrainCrystalPlasticity::solveStress().

_acc_slip

MaterialProperty<Real>& FiniteStrainCrystalPlasticity::_acc_slip

protectedinherited

Definition at line 307 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::initQpStatefulProperties(), and FiniteStrainCrystalPlasticity::postSolveStatevar().

_acc_slip_old

const MaterialProperty<Real>& FiniteStrainCrystalPlasticity::_acc_slip_old

protectedinherited

Definition at line 308 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::preSolveStatevar().

_accslip_tmp

Real FiniteStrainCrystalPlasticity::_accslip_tmp

protectedinherited

Definition at line 336 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::postSolveStatevar(), and FiniteStrainCrystalPlasticity::updateGss().

_accslip_tmp_old

Real FiniteStrainCrystalPlasticity::_accslip_tmp_old

protectedinherited

Definition at line 336 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::postSolveStatevar(), FiniteStrainCrystalPlasticity::preSolveStatevar(), and FiniteStrainCrystalPlasticity::updateGss().

_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>& FiniteStrainCrystalPlasticity::_crysrot

protectedinherited

Definition at line 317 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::calc_schmid_tensor(), and FiniteStrainCrystalPlasticity::postSolveQp().

_deformation_gradient

const MaterialProperty<RankTwoTensor>& FiniteStrainCrystalPlasticity::_deformation_gradient

protectedinherited

Definition at line 311 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::computeQpStress(), FiniteStrainCrystalPlasticity::postSolveQp(), and FiniteStrainCrystalPlasticity::preSolveQp().

_deformation_gradient_old

const MaterialProperty<RankTwoTensor>& FiniteStrainCrystalPlasticity::_deformation_gradient_old

protectedinherited

Definition at line 312 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::computeQpStress().

_delta_dfgrd

RankTwoTensor FiniteStrainCrystalPlasticity::_delta_dfgrd

protectedinherited

Flag to check whether convergence is achieved.

Used for substepping; Uniformly divides the increment in deformation gradient

Definition at line 349 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::computeQpStress(), FiniteStrainCrystalPlasticity::FiniteStrainCrystalPlasticity(), and FiniteStrainCrystalPlasticity::preSolveQp().

_dfgrd_scale_factor

Real FiniteStrainCrystalPlasticity::_dfgrd_scale_factor

protectedinherited

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

Definition at line 351 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::computeQpStress(), and FiniteStrainCrystalPlasticity::preSolveQp().

_dfgrd_tmp

RankTwoTensor FiniteStrainCrystalPlasticity::_dfgrd_tmp

protectedinherited

Definition at line 330 of file FiniteStrainCrystalPlasticity.h.

Referenced by calcDtauDsliprate(), FiniteStrainCrystalPlasticity::calcJacobian(), FiniteStrainCrystalPlasticity::calcResidual(), calcResidualSlipRate(), and FiniteStrainCrystalPlasticity::preSolveQp().

_dfgrd_tmp_old

RankTwoTensor FiniteStrainCrystalPlasticity::_dfgrd_tmp_old

protectedinherited

Definition at line 349 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::computeQpStress(), and FiniteStrainCrystalPlasticity::preSolveQp().

_dgss_dsliprate

DenseMatrix<Real> FiniteStrainCrystalPlasticity::_dgss_dsliprate

protectedinherited

Definition at line 342 of file FiniteStrainCrystalPlasticity.h.

Referenced by calcDgssDsliprate(), and FiniteStrainCrystalPlasticity::updateGss().

_dslipdtau

DenseVector<Real> FiniteStrainCrystalPlasticity::_dslipdtau

protectedinherited

Definition at line 332 of file FiniteStrainCrystalPlasticity.h.

Referenced by calcDtauDsliprate(), FiniteStrainCrystalPlasticity::calcJacobian(), getSlipIncrements(), and FiniteStrainCrystalPlasticity::getSlipIncrements().

_dsliprate_dgss

DenseVector<Real> FiniteStrainCPSlipRateRes::_dsliprate_dgss

protected

Definition at line 91 of file FiniteStrainCPSlipRateRes.h.

Referenced by calcDgssDsliprate(), and getSlipIncrements().

_dsliprate_dsliprate

DenseMatrix<Real> FiniteStrainCPSlipRateRes::_dsliprate_dsliprate

protected

Definition at line 93 of file FiniteStrainCPSlipRateRes.h.

Referenced by calcDgssDsliprate(), calcDtauDsliprate(), and calcJacobianSlipRate().

_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>& FiniteStrainCrystalPlasticity::_elasticity_tensor

protectedinherited

Elasticity tensor material property.

Definition at line 316 of file FiniteStrainCrystalPlasticity.h.

Referenced by calcDtauDsliprate(), FiniteStrainCrystalPlasticity::calcJacobian(), FiniteStrainCrystalPlasticity::calcResidual(), calcResidualSlipRate(), FiniteStrainCrystalPlasticity::elasticTangentModuli(), FiniteStrainCrystalPlasticity::elastoPlasticTangentModuli(), and FiniteStrainCrystalPlasticity::postSolveQp().

_elasticity_tensor_name

const std::string FiniteStrainCrystalPlasticity::_elasticity_tensor_name

protectedinherited

Name of the elasticity tensor material property.

Definition at line 314 of file FiniteStrainCrystalPlasticity.h.

_err_tol

bool FiniteStrainCrystalPlasticity::_err_tol

protectedinherited

Definition at line 346 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::calc_resid_jacob(), calcResidJacobSlipRate(), FiniteStrainCrystalPlasticity::calcResidual(), calcResidualSlipRate(), FiniteStrainCrystalPlasticity::computeQpStress(), FiniteStrainCrystalPlasticity::FiniteStrainCrystalPlasticity(), getSlipIncrements(), FiniteStrainCrystalPlasticity::getSlipIncrements(), lineSearchUpdateSlipRate(), FiniteStrainCrystalPlasticity::postSolveQp(), FiniteStrainCrystalPlasticity::preSolveQp(), FiniteStrainCrystalPlasticity::solveQp(), solveStatevar(), FiniteStrainCrystalPlasticity::solveStatevar(), solveStress(), and FiniteStrainCrystalPlasticity::solveStress().

_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 FiniteStrainCrystalPlasticity::_fe

protectedinherited

Definition at line 331 of file FiniteStrainCrystalPlasticity.h.

Referenced by calcDtauDsliprate(), FiniteStrainCrystalPlasticity::calcJacobian(), FiniteStrainCrystalPlasticity::calcResidual(), calcResidualSlipRate(), FiniteStrainCrystalPlasticity::elastoPlasticTangentModuli(), and FiniteStrainCrystalPlasticity::postSolveQp().

_first_step_iter

bool FiniteStrainCrystalPlasticity::_first_step_iter

protectedinherited

Flags to reset variables and reinitialize variables.

Definition at line 353 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::computeQpStress(), FiniteStrainCrystalPlasticity::FiniteStrainCrystalPlasticity(), FiniteStrainCrystalPlasticity::preSolveStatevar(), and FiniteStrainCrystalPlasticity::preSolveStress().

_first_substep

bool FiniteStrainCrystalPlasticity::_first_substep

protectedinherited

Definition at line 353 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::computeQpStress(), FiniteStrainCrystalPlasticity::FiniteStrainCrystalPlasticity(), and FiniteStrainCrystalPlasticity::preSolveQp().

_flowprops

std::vector<Real> FiniteStrainCrystalPlasticity::_flowprops

protectedinherited

Definition at line 227 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::getFlowRateParams().

_fp

MaterialProperty<RankTwoTensor>& FiniteStrainCrystalPlasticity::_fp

protectedinherited

Definition at line 299 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::initQpStatefulProperties(), and FiniteStrainCrystalPlasticity::postSolveStress().

_fp_inv

RankTwoTensor FiniteStrainCrystalPlasticity::_fp_inv

protectedinherited

Definition at line 331 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::calcResidual(), calcResidualSlipRate(), FiniteStrainCrystalPlasticity::internalVariableUpdateNRiteration(), FiniteStrainCrystalPlasticity::postSolveStress(), and FiniteStrainCrystalPlasticity::preSolveStress().

_fp_old

const MaterialProperty<RankTwoTensor>& FiniteStrainCrystalPlasticity::_fp_old

protectedinherited

Definition at line 300 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::preSolveStress().

_fp_old_inv

RankTwoTensor FiniteStrainCrystalPlasticity::_fp_old_inv

protectedinherited

Definition at line 331 of file FiniteStrainCrystalPlasticity.h.

Referenced by calcDtauDsliprate(), FiniteStrainCrystalPlasticity::calcJacobian(), FiniteStrainCrystalPlasticity::calcResidual(), calcResidualSlipRate(), FiniteStrainCrystalPlasticity::postSolveStress(), and FiniteStrainCrystalPlasticity::preSolveStress().

_fp_prev_inv

RankTwoTensor FiniteStrainCrystalPlasticity::_fp_prev_inv

protectedinherited

Definition at line 331 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::calcResidual(), FiniteStrainCrystalPlasticity::internalVariableUpdateNRiteration(), and FiniteStrainCrystalPlasticity::preSolveStress().

_gen_rndm_stress_flag

bool FiniteStrainCrystalPlasticity::_gen_rndm_stress_flag

protectedinherited

Definition at line 270 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::postSolveQp().

_gprops

std::vector<Real> FiniteStrainCrystalPlasticity::_gprops

protectedinherited

Definition at line 225 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::getInitSlipSysRes().

_gss

MaterialProperty<std::vector<Real> >& FiniteStrainCrystalPlasticity::_gss

protectedinherited

Definition at line 305 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::assignSlipSysRes(), FiniteStrainCrystalPlasticity::getInitSlipSysRes(), FiniteStrainCrystalPlasticity::postSolveStatevar(), and FiniteStrainCrystalPlasticity::readFileInitSlipSysRes().

_gss_old

const MaterialProperty<std::vector<Real> >& FiniteStrainCrystalPlasticity::_gss_old

protectedinherited

Definition at line 306 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::preSolveStatevar(), and FiniteStrainCrystalPlasticity::updateGss().

_gss_tmp

std::vector<Real> FiniteStrainCrystalPlasticity::_gss_tmp

protectedinherited

Definition at line 337 of file FiniteStrainCrystalPlasticity.h.

Referenced by getSlipIncrements(), FiniteStrainCrystalPlasticity::getSlipIncrements(), FiniteStrainCrystalPlasticity::postSolveStatevar(), FiniteStrainCrystalPlasticity::preSolveStatevar(), FiniteStrainCrystalPlasticity::solveStatevar(), and FiniteStrainCrystalPlasticity::updateGss().

_gss_tmp_old

std::vector<Real> FiniteStrainCrystalPlasticity::_gss_tmp_old

protectedinherited

Definition at line 338 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::postSolveStatevar(), FiniteStrainCrystalPlasticity::preSolveStatevar(), and FiniteStrainCrystalPlasticity::updateGss().

_gtol

Real FiniteStrainCrystalPlasticity::_gtol

protectedinherited

Internal variable update equation tolerance.

Definition at line 249 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::solveStatevar().

_h0

Real FiniteStrainCrystalPlasticity::_h0

protectedinherited

Definition at line 325 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::getHardnessParams(), and FiniteStrainCrystalPlasticity::updateGss().

_hprops

std::vector<Real> FiniteStrainCrystalPlasticity::_hprops

protectedinherited

Definition at line 226 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::getHardnessParams(), and FiniteStrainCrystalPlasticity::updateGss().

_initial_stress_fcn

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

protectedinherited

initial stress components

Definition at line 56 of file ComputeGeneralStressBase.h.

_input_rndm_scale_var

bool FiniteStrainCrystalPlasticity::_input_rndm_scale_var

protectedinherited

Input option for scaling variable to generate random stress when convergence fails.

Definition at line 273 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::postSolveQp().

_intvar_read_type

MooseEnum FiniteStrainCrystalPlasticity::_intvar_read_type

protectedinherited

Read from options for initial values of internal variables.

Definition at line 265 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::FiniteStrainCrystalPlasticity(), and FiniteStrainCrystalPlasticity::initSlipSysProps().

_jacob

DenseMatrix<Real> FiniteStrainCPSlipRateRes::_jacob

protected

Definition at line 92 of file FiniteStrainCPSlipRateRes.h.

Referenced by calcJacobianSlipRate(), and calcUpdate().

_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(), FiniteStrainUObasedCP::elasticTangentModuli(), FiniteStrainUObasedCP::elastoPlasticTangentModuli(), ComputeMultipleInelasticStressBase::finiteStrainRotation(), ComputeMultiPlasticityStress::postReturnMap(), FiniteStrainCrystalPlasticity::postSolveQp(), FiniteStrainCrystalPlasticity::preSolveQp(), and ComputeMultipleInelasticStressBase::updateQpStateSingleModel().

_lag_e

MaterialProperty<RankTwoTensor>& FiniteStrainCrystalPlasticity::_lag_e

protectedinherited

Definition at line 303 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::initQpStatefulProperties(), and FiniteStrainCrystalPlasticity::postSolveQp().

_lag_e_old

const MaterialProperty<RankTwoTensor>& FiniteStrainCrystalPlasticity::_lag_e_old

protectedinherited

Definition at line 304 of file FiniteStrainCrystalPlasticity.h.

_last_step_iter

bool FiniteStrainCrystalPlasticity::_last_step_iter

protectedinherited

Definition at line 353 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::computeQpStress(), FiniteStrainCrystalPlasticity::FiniteStrainCrystalPlasticity(), FiniteStrainCrystalPlasticity::postSolveStatevar(), and FiniteStrainCrystalPlasticity::postSolveStress().

_lsrch_max_iter

unsigned int FiniteStrainCrystalPlasticity::_lsrch_max_iter

protectedinherited

Line search bisection method maximum iteration number.

Definition at line 294 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::line_search_update(), and lineSearchUpdateSlipRate().

_lsrch_method

MooseEnum FiniteStrainCrystalPlasticity::_lsrch_method

protectedinherited

Definition at line 297 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::line_search_update(), and lineSearchUpdateSlipRate().

_lsrch_tol

Real FiniteStrainCrystalPlasticity::_lsrch_tol

protectedinherited

Line search bisection method tolerance.

Definition at line 291 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::line_search_update(), and lineSearchUpdateSlipRate().

_max_substep_iter

unsigned int FiniteStrainCrystalPlasticity::_max_substep_iter

protectedinherited

Maximum number of substep iterations.

Definition at line 282 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::computeQpStress(), FiniteStrainCrystalPlasticity::postSolveStatevar(), FiniteStrainCrystalPlasticity::postSolveStress(), FiniteStrainCrystalPlasticity::preSolveQp(), FiniteStrainCrystalPlasticity::preSolveStatevar(), FiniteStrainCrystalPlasticity::preSolveStress(), and FiniteStrainCrystalPlasticity::updateGss().

_maxiter

unsigned int FiniteStrainCrystalPlasticity::_maxiter

protectedinherited

Maximum number of iterations for stress update.

Definition at line 254 of file FiniteStrainCrystalPlasticity.h.

Referenced by solveStress(), and FiniteStrainCrystalPlasticity::solveStress().

_maxiterg

unsigned int FiniteStrainCrystalPlasticity::_maxiterg

protectedinherited

Maximum number of iterations for internal variable update.

Definition at line 256 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::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 FiniteStrainCrystalPlasticity::_min_lsrch_step

protectedinherited

Minimum line search step size.

Definition at line 288 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::line_search_update(), and lineSearchUpdateSlipRate().

_mo

DenseVector<Real> FiniteStrainCrystalPlasticity::_mo

protectedinherited

Definition at line 319 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::calc_schmid_tensor(), and FiniteStrainCrystalPlasticity::getSlipSystems().

_no

DenseVector<Real> FiniteStrainCrystalPlasticity::_no

protectedinherited

Definition at line 320 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::calc_schmid_tensor(), and FiniteStrainCrystalPlasticity::getSlipSystems().

_nss

const unsigned int FiniteStrainCrystalPlasticity::_nss

protectedinherited

Number of slip system resistance.

Definition at line 223 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::assignSlipSysRes(), FiniteStrainCrystalPlasticity::calc_schmid_tensor(), calcDgssDsliprate(), calcDtauDsliprate(), FiniteStrainCrystalPlasticity::calcJacobian(), calcJacobianSlipRate(), calcResidDotProdUpdate(), calcResidNorm(), FiniteStrainCrystalPlasticity::calcResidual(), calcResidualSlipRate(), calcUpdate(), FiniteStrainCrystalPlasticity::FiniteStrainCrystalPlasticity(), FiniteStrainCrystalPlasticity::getFlowRateParams(), FiniteStrainCrystalPlasticity::getInitSlipSysRes(), getSlipIncrements(), FiniteStrainCrystalPlasticity::getSlipIncrements(), FiniteStrainCrystalPlasticity::getSlipSystems(), FiniteStrainCrystalPlasticity::readFileFlowRateParams(), FiniteStrainCrystalPlasticity::readFileInitSlipSysRes(), FiniteStrainCrystalPlasticity::solveStatevar(), solveStress(), and FiniteStrainCrystalPlasticity::updateGss().

_num_slip_sys_flowrate_props

unsigned int FiniteStrainCrystalPlasticity::_num_slip_sys_flowrate_props

protectedinherited

Number of slip system flow rate parameters.

Definition at line 259 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::getFlowRateParams(), and FiniteStrainCrystalPlasticity::readFileFlowRateParams().

_num_slip_sys_props

unsigned int FiniteStrainCrystalPlasticity::_num_slip_sys_props

protectedinherited

Number of slip system specific properties provided in the file containing slip system normals and directions.

Definition at line 268 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::FiniteStrainCrystalPlasticity(), and FiniteStrainCrystalPlasticity::getSlipSystems().

_pk2

MaterialProperty<RankTwoTensor>& FiniteStrainCrystalPlasticity::_pk2

protectedinherited

Definition at line 301 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::initQpStatefulProperties(), FiniteStrainCrystalPlasticity::postSolveQp(), and FiniteStrainCrystalPlasticity::postSolveStress().

_pk2_old

const MaterialProperty<RankTwoTensor>& FiniteStrainCrystalPlasticity::_pk2_old

protectedinherited

Definition at line 302 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::preSolveStress().

_pk2_tmp

RankTwoTensor FiniteStrainCrystalPlasticity::_pk2_tmp

protectedinherited

Definition at line 335 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::calcResidual(), calcResidualSlipRate(), FiniteStrainCrystalPlasticity::elastoPlasticTangentModuli(), FiniteStrainCrystalPlasticity::FiniteStrainCrystalPlasticity(), FiniteStrainCrystalPlasticity::line_search_update(), FiniteStrainCrystalPlasticity::postSolveStress(), FiniteStrainCrystalPlasticity::preSolveStress(), and FiniteStrainCrystalPlasticity::solveStress().

_pk2_tmp_old

RankTwoTensor FiniteStrainCrystalPlasticity::_pk2_tmp_old

protectedinherited

Definition at line 335 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::postSolveStress(), and FiniteStrainCrystalPlasticity::preSolveStress().

_r

Real FiniteStrainCrystalPlasticity::_r

protectedinherited

Definition at line 328 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::getHardnessParams(), and FiniteStrainCrystalPlasticity::updateGss().

_read_from_slip_sys_file

bool FiniteStrainCrystalPlasticity::_read_from_slip_sys_file

protectedinherited

Definition at line 344 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::FiniteStrainCrystalPlasticity(), and FiniteStrainCrystalPlasticity::getSlipSystems().

_resid

DenseVector<Real> FiniteStrainCPSlipRateRes::_resid

protected

Definition at line 89 of file FiniteStrainCPSlipRateRes.h.

Referenced by calcResidDotProdUpdate(), calcResidNorm(), calcResidualSlipRate(), calcUpdate(), and solveStress().

_rndm_scale_var

Real FiniteStrainCrystalPlasticity::_rndm_scale_var

protectedinherited

Scaling value.

Definition at line 276 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::postSolveQp().

_rndm_seed

unsigned int FiniteStrainCrystalPlasticity::_rndm_seed

protectedinherited

Seed value.

Definition at line 279 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::FiniteStrainCrystalPlasticity().

_rtol

Real FiniteStrainCrystalPlasticity::_rtol

protectedinherited

Stress residual equation relative tolerance.

Definition at line 245 of file FiniteStrainCrystalPlasticity.h.

Referenced by solveStress(), and FiniteStrainCrystalPlasticity::solveStress().

_s0

std::vector<RankTwoTensor> FiniteStrainCrystalPlasticity::_s0

protectedinherited

Definition at line 333 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::calc_schmid_tensor(), calcDtauDsliprate(), FiniteStrainCrystalPlasticity::calcJacobian(), FiniteStrainCrystalPlasticity::calcResidual(), and calcResidualSlipRate().

_slip_incr

DenseVector<Real> FiniteStrainCrystalPlasticity::_slip_incr

protectedinherited

Definition at line 332 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::calcResidual(), calcResidualSlipRate(), FiniteStrainCrystalPlasticity::getSlipIncrements(), and FiniteStrainCrystalPlasticity::updateGss().

_slip_incr_tol

Real FiniteStrainCrystalPlasticity::_slip_incr_tol

protectedinherited

Slip increment tolerance.

Definition at line 251 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::getSlipIncrements().

_slip_rate

DenseVector<Real> FiniteStrainCPSlipRateRes::_slip_rate

protected

Definition at line 90 of file FiniteStrainCPSlipRateRes.h.

Referenced by calcResidualSlipRate(), lineSearchUpdateSlipRate(), preSolveStress(), and solveStress().

_slip_sys_file_name

std::string FiniteStrainCrystalPlasticity::_slip_sys_file_name

protectedinherited

File should contain slip plane normal and direction. See test.

Definition at line 230 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::getSlipSystems().

_slip_sys_flow_prop_file_name

std::string FiniteStrainCrystalPlasticity::_slip_sys_flow_prop_file_name

protectedinherited

File should contain values of the flow rate equation parameters.

Values for every slip system must be provided. Should have the same order of slip systens as in slip_sys_file. See test. The option of reading all the properties from .i is still present.

Definition at line 239 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::initSlipSysProps(), and FiniteStrainCrystalPlasticity::readFileFlowRateParams().

_slip_sys_hard_prop_file_name

std::string FiniteStrainCrystalPlasticity::_slip_sys_hard_prop_file_name

protectedinherited

The hardening parameters in this class are read from .i file. The user can override to read from file.

Definition at line 242 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::initSlipSysProps().

_slip_sys_props

DenseVector<Real> FiniteStrainCrystalPlasticity::_slip_sys_props

protectedinherited

Definition at line 340 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::assignSlipSysRes(), FiniteStrainCrystalPlasticity::FiniteStrainCrystalPlasticity(), and FiniteStrainCrystalPlasticity::getSlipSystems().

_slip_sys_res_prop_file_name

std::string FiniteStrainCrystalPlasticity::_slip_sys_res_prop_file_name

protectedinherited

File should contain initial values of the slip system resistances.

Definition at line 233 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::readFileInitSlipSysRes().

_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(), FiniteStrainUObasedCP::initQpStatefulProperties(), FiniteStrainHyperElasticViscoPlastic::initQpStatefulProperties(), ComputeGeneralStressBase::initQpStatefulProperties(), ComputeMultiPlasticityStress::postReturnMap(), FiniteStrainCrystalPlasticity::postSolveQp(), FiniteStrainUObasedCP::postSolveQp(), FiniteStrainHyperElasticViscoPlastic::postSolveQp(), ComputeSmearedCrackingStress::updateCrackingStateAndStress(), ComputeCreepPlasticityStress::updateQpState(), ComputeMultipleInelasticStress::updateQpState(), and ComputeMultipleInelasticStressBase::updateQpStateSingleModel().

_tan_mod_type

MooseEnum FiniteStrainCrystalPlasticity::_tan_mod_type

protectedinherited

Type of tangent moduli calculation.

Definition at line 262 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::calcTangentModuli().

_tau

DenseVector<Real> FiniteStrainCrystalPlasticity::_tau

protectedinherited

Definition at line 332 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::calcResidual(), calcResidualSlipRate(), getSlipIncrements(), and FiniteStrainCrystalPlasticity::getSlipIncrements().

_tau_init

Real FiniteStrainCrystalPlasticity::_tau_init

protectedinherited

Definition at line 327 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::getHardnessParams().

_tau_sat

Real FiniteStrainCrystalPlasticity::_tau_sat

protectedinherited

Definition at line 326 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::getHardnessParams(), and FiniteStrainCrystalPlasticity::updateGss().

_update_rot

MaterialProperty<RankTwoTensor>& FiniteStrainCrystalPlasticity::_update_rot

protectedinherited

Definition at line 309 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::initQpStatefulProperties(), and FiniteStrainCrystalPlasticity::postSolveQp().

_use_line_search

bool FiniteStrainCrystalPlasticity::_use_line_search

protectedinherited

Flag to activate line serach.

Definition at line 285 of file FiniteStrainCrystalPlasticity.h.

Referenced by solveStress(), and FiniteStrainCrystalPlasticity::solveStress().

_xm

DenseVector<Real> FiniteStrainCrystalPlasticity::_xm

protectedinherited

Definition at line 323 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::getFlowRateParams(), getSlipIncrements(), FiniteStrainCrystalPlasticity::getSlipIncrements(), and FiniteStrainCrystalPlasticity::readFileFlowRateParams().

---

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

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