FiniteStrainCrystalPlasticity Class Reference

FiniteStrainCrystalPlasticity 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 <FiniteStrainCrystalPlasticity.h>

Inheritance diagram for FiniteStrainCrystalPlasticity:

Public Types
typedef DerivativeMaterialPropertyNameInterface::SymbolName	SymbolName

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

Static Public Member Functions
static InputParameters	validParams ()

Protected Member Functions
virtual void	computeQpStress ()
	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	getSlipIncrements ()
	This function updates the slip increments.
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	solveStatevar ()
	This function solves internal variables.
virtual void	postSolveStatevar ()
	This function update internal variable after solve.
virtual void	preSolveStress ()
	This function set variables for stress solve.
virtual void	solveStress ()
	This function solves for stress, updates plastic deformation gradient.
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
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

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

Definition at line 21 of file FiniteStrainCrystalPlasticity.h.

Constructor & Destructor Documentation

FiniteStrainCrystalPlasticity()

FiniteStrainCrystalPlasticity::FiniteStrainCrystalPlasticity

(

const InputParameters &

parameters

)

Definition at line 102 of file FiniteStrainCrystalPlasticity.C.

  : ComputeStressBase(parameters),
    _nss(getParam<int>("nss")),
    _gprops(getParam<std::vector<Real>>("gprops")),
    _hprops(getParam<std::vector<Real>>("hprops")),
    _flowprops(getParam<std::vector<Real>>("flowprops")),
    _slip_sys_file_name(getParam<FileName>("slip_sys_file_name")),
    _slip_sys_res_prop_file_name(getParam<FileName>("slip_sys_res_prop_file_name")),
    _slip_sys_flow_prop_file_name(getParam<FileName>("slip_sys_flow_prop_file_name")),
    _slip_sys_hard_prop_file_name(getParam<FileName>("slip_sys_hard_prop_file_name")),
    _rtol(getParam<Real>("rtol")),
    _abs_tol(getParam<Real>("abs_tol")),
    _gtol(getParam<Real>("gtol")),
    _slip_incr_tol(getParam<Real>("slip_incr_tol")),
    _maxiter(getParam<unsigned int>("maxiter")),
    _maxiterg(getParam<unsigned int>("maxitergss")),
    _num_slip_sys_flowrate_props(getParam<unsigned int>("num_slip_sys_flowrate_props")),
    _tan_mod_type(getParam<MooseEnum>("tan_mod_type")),
    _intvar_read_type(getParam<MooseEnum>("intvar_read_type")),
    _num_slip_sys_props(getParam<unsigned int>("num_slip_sys_props")),
    _gen_rndm_stress_flag(getParam<bool>("gen_random_stress_flag")),
    _input_rndm_scale_var(getParam<bool>("input_random_scaling_var")),
    _rndm_scale_var(getParam<Real>("random_scaling_var")),
    _rndm_seed(getParam<unsigned int>("random_seed")),
    _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
    _lag_e_old(
        getMaterialPropertyOld<RankTwoTensor>("lage")), // Lagrangian strain of previous increment
    _gss(declareProperty<std::vector<Real>>("gss")),    // Slip system resistances
    _gss_old(getMaterialPropertyOld<std::vector<Real>>(
        "gss")),                                  // Slip system resistances of previous increment
    _acc_slip(declareProperty<Real>("acc_slip")), // Accumulated slip
    _acc_slip_old(
        getMaterialPropertyOld<Real>("acc_slip")), // Accumulated slip of previous increment
    _update_rot(declareProperty<RankTwoTensor>(
        "update_rot")), // Rotation tensor considering material rotation and crystal orientation
    _deformation_gradient(getMaterialProperty<RankTwoTensor>("deformation_gradient")),
    _deformation_gradient_old(getMaterialPropertyOld<RankTwoTensor>("deformation_gradient")),
    _elasticity_tensor_name(_base_name + "elasticity_tensor"),
    _elasticity_tensor(getMaterialPropertyByName<RankFourTensor>(_elasticity_tensor_name)),
    _crysrot(getMaterialProperty<RankTwoTensor>("crysrot")),
    _mo(_nss * LIBMESH_DIM),
    _no(_nss * LIBMESH_DIM),
    _slip_incr(_nss),
    _tau(_nss),
    _dslipdtau(_nss),
    _s0(_nss),
    _gss_tmp(_nss),
    _gss_tmp_old(_nss),
    _dgss_dsliprate(_nss, _nss)
{
  _err_tol = false;
 
  if (_num_slip_sys_props > 0)
    _slip_sys_props.resize(_nss * _num_slip_sys_props);
 
  _pk2_tmp.zero();
  _delta_dfgrd.zero();
 
  _first_step_iter = false;
  _last_step_iter = false;
  // Initialize variables in the first iteration of substepping
  _first_substep = true;
 
  _read_from_slip_sys_file = false;
  if (_intvar_read_type == "slip_sys_file")
    _read_from_slip_sys_file = true;
 
  if (_read_from_slip_sys_file && !(_num_slip_sys_props > 0))
    mooseError("Crystal Plasticity Error: Specify number of internal variable's initial values to "
               "be read from slip system file");
 
  getSlipSystems();
 
  RankTwoTensor::initRandom(_rndm_seed);
}

Member Function Documentation

assignSlipSysRes()

void FiniteStrainCrystalPlasticity::assignSlipSysRes ( )

protectedvirtual

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 initSlipSysProps().

calc_resid_jacob()

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

protectedvirtual

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 solveStress().

calc_schmid_tensor()

void FiniteStrainCrystalPlasticity::calc_schmid_tensor ( )

protected

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 preSolveQp().

calcJacobian()

void FiniteStrainCrystalPlasticity::calcJacobian ( RankFourTensor &  jac )

protectedvirtual

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 calc_resid_jacob().

calcResidual()

void FiniteStrainCrystalPlasticity::calcResidual ( RankTwoTensor &  resid )

protectedvirtual

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 calc_resid_jacob(), and line_search_update().

calcTangentModuli()

RankFourTensor FiniteStrainCrystalPlasticity::calcTangentModuli ( )

protectedvirtual

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 postSolveQp().

computeQpElasticityTensor()

void FiniteStrainCrystalPlasticity::computeQpElasticityTensor ( )

protectedvirtual

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 ( )

protectedvirtual

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 ( )

protectedvirtual

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 calcTangentModuli().

elastoPlasticTangentModuli()

RankFourTensor FiniteStrainCrystalPlasticity::elastoPlasticTangentModuli ( )

protectedvirtual

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 calcTangentModuli().

get_current_rotation()

RankTwoTensor FiniteStrainCrystalPlasticity::get_current_rotation ( const RankTwoTensor &  a )

protected

This function perform RU decomposition to obtain the rotation tensor.

Definition at line 985 of file FiniteStrainCrystalPlasticity.C.

{
  return getMatRot(a);
}

Referenced by postSolveQp().

getFlowRateParams()

void FiniteStrainCrystalPlasticity::getFlowRateParams ( )

protectedvirtual

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 initSlipSysProps().

getHardnessParams()

void FiniteStrainCrystalPlasticity::getHardnessParams ( )

protectedvirtual

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 initSlipSysProps().

getInitSlipSysRes()

void FiniteStrainCrystalPlasticity::getInitSlipSysRes ( )

protectedvirtual

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 initSlipSysProps().

getMatRot()

RankTwoTensor FiniteStrainCrystalPlasticity::getMatRot ( const RankTwoTensor &  a )

protected

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 get_current_rotation().

getSlipIncrements()

void FiniteStrainCrystalPlasticity::getSlipIncrements ( )

protectedvirtual

This function updates the slip increments.

And derivative of slip w.r.t. resolved shear stress.

Reimplemented in FiniteStrainCPSlipRateRes.

Definition at line 961 of file FiniteStrainCrystalPlasticity.C.

{
  for (unsigned int i = 0; i < _nss; ++i)
  {
    _slip_incr(i) = _a0(i) * std::pow(std::abs(_tau(i) / _gss_tmp[i]), 1.0 / _xm(i)) *
                    std::copysign(1.0, _tau(i)) * _dt;
    if (std::abs(_slip_incr(i)) > _slip_incr_tol)
    {
      _err_tol = true;
#ifdef DEBUG
      mooseWarning("Maximum allowable slip increment exceeded ", std::abs(_slip_incr(i)));
#endif
      return;
    }
  }
 
  for (unsigned int i = 0; i < _nss; ++i)
    _dslipdtau(i) = _a0(i) / _xm(i) *
                    std::pow(std::abs(_tau(i) / _gss_tmp[i]), 1.0 / _xm(i) - 1.0) / _gss_tmp[i] *
                    _dt;
}

Referenced by calcResidual(), and FiniteStrainCPSlipRateRes::getSlipIncrements().

getSlipSystems()

void FiniteStrainCrystalPlasticity::getSlipSystems ( )

protectedvirtual

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().

initAdditionalProps()

void FiniteStrainCrystalPlasticity::initAdditionalProps ( )

protectedvirtual

This function initializes additional parameters.

Definition at line 484 of file FiniteStrainCrystalPlasticity.C.

{
}

Referenced by initQpStatefulProperties().

initQpStatefulProperties()

void FiniteStrainCrystalPlasticity::initQpStatefulProperties ( )

protectedvirtual

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 ( )

protectedvirtual

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 initQpStatefulProperties().

internalVariableUpdateNRiteration()

void FiniteStrainCrystalPlasticity::internalVariableUpdateNRiteration ( )

protected

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 solveStress().

line_search_update()

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

protected

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 solveStress().

postSolveQp()

void FiniteStrainCrystalPlasticity::postSolveQp ( )

protectedvirtual

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 computeQpStress().

postSolveStatevar()

void FiniteStrainCrystalPlasticity::postSolveStatevar ( )

protectedvirtual

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 solveQp().

postSolveStress()

void FiniteStrainCrystalPlasticity::postSolveStress ( )

protectedvirtual

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 FiniteStrainCPSlipRateRes::solveStatevar(), and solveStatevar().

preSolveQp()

void FiniteStrainCrystalPlasticity::preSolveQp ( )

protectedvirtual

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 computeQpStress().

preSolveStatevar()

void FiniteStrainCrystalPlasticity::preSolveStatevar ( )

protectedvirtual

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 solveQp().

preSolveStress()

void FiniteStrainCrystalPlasticity::preSolveStress ( )

protectedvirtual

This function set variables for stress solve.

Reimplemented in FiniteStrainCPSlipRateRes.

Definition at line 706 of file FiniteStrainCrystalPlasticity.C.

{
  if (_max_substep_iter == 1) // No substepping
  {
    _pk2_tmp = _pk2_old[_qp];
    _fp_old_inv = _fp_old[_qp].inverse();
    _fp_inv = _fp_old_inv;
    _fp_prev_inv = _fp_inv;
  }
  else
  {
    if (_first_step_iter)
    {
      _pk2_tmp = _pk2_tmp_old = _pk2_old[_qp];
      _fp_old_inv = _fp_old[_qp].inverse();
    }
    else
      _pk2_tmp = _pk2_tmp_old;
 
    _fp_inv = _fp_old_inv;
    _fp_prev_inv = _fp_inv;
  }
}

Referenced by FiniteStrainCPSlipRateRes::preSolveStress(), and solveStatevar().

readFileFlowRateParams()

void FiniteStrainCrystalPlasticity::readFileFlowRateParams ( )

protectedvirtual

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 initSlipSysProps().

readFileHardnessParams()

void FiniteStrainCrystalPlasticity::readFileHardnessParams ( )

protectedvirtual

This function read hardness parameters from file.

Definition at line 406 of file FiniteStrainCrystalPlasticity.C.

{
}

Referenced by initSlipSysProps().

readFileInitSlipSysRes()

void FiniteStrainCrystalPlasticity::readFileInitSlipSysRes ( )

protectedvirtual

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 initSlipSysProps().

solveQp()

void FiniteStrainCrystalPlasticity::solveQp ( )

protectedvirtual

This function solves stress and internal variables.

Definition at line 578 of file FiniteStrainCrystalPlasticity.C.

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

Referenced by computeQpStress().

solveStatevar()

void FiniteStrainCrystalPlasticity::solveStatevar ( )

protectedvirtual

This function solves internal variables.

Reimplemented in FiniteStrainCPSlipRateRes.

Definition at line 641 of file FiniteStrainCrystalPlasticity.C.

{
  Real gmax, gdiff;
  unsigned int iterg;
  std::vector<Real> gss_prev(_nss);
 
  gmax = 1.1 * _gtol;
  iterg = 0;
 
  while (gmax > _gtol && iterg < _maxiterg) // Check for slip system resistance update tolerance
  {
    preSolveStress();
    solveStress();
    if (_err_tol)
      return;
    postSolveStress();
 
    gss_prev = _gss_tmp;
 
    update_slip_system_resistance(); // Update slip system resistance
 
    gmax = 0.0;
    for (unsigned i = 0; i < _nss; ++i)
    {
      gdiff = std::abs(gss_prev[i] - _gss_tmp[i]); // Calculate increment size
 
      if (gdiff > gmax)
        gmax = gdiff;
    }
    iterg++;
  }
 
  if (iterg == _maxiterg)
  {
#ifdef DEBUG
    mooseWarning("FiniteStrainCrystalPLasticity: Hardness Integration error gmax", gmax, "\n");
#endif
    _err_tol = true;
  }
}

Referenced by solveQp().

solveStress()

void FiniteStrainCrystalPlasticity::solveStress ( )

protectedvirtual

This function solves for stress, updates plastic deformation gradient.

Reimplemented in FiniteStrainCPSlipRateRes.

Definition at line 731 of file FiniteStrainCrystalPlasticity.C.

{
  unsigned int iter = 0;
  RankTwoTensor resid, dpk2;
  RankFourTensor jac;
  Real rnorm, rnorm0, rnorm_prev;
 
  calc_resid_jacob(resid, jac); // Calculate stress residual
  if (_err_tol)
  {
#ifdef DEBUG
    mooseWarning(
        "FiniteStrainCrystalPLasticity: Slip increment exceeds tolerance - Element number ",
        _current_elem->id(),
        " Gauss point = ",
        _qp);
#endif
    return;
  }
 
  rnorm = resid.L2norm();
  rnorm0 = rnorm;
 
  while (rnorm > _rtol * rnorm0 && rnorm0 > _abs_tol &&
         iter < _maxiter) // Check for stress residual tolerance
  {
    dpk2 = -jac.invSymm() * resid; // Calculate stress increment
    _pk2_tmp = _pk2_tmp + dpk2;    // Update stress
    calc_resid_jacob(resid, jac);
    internalVariableUpdateNRiteration(); // update _fp_prev_inv
 
    if (_err_tol)
    {
#ifdef DEBUG
      mooseWarning(
          "FiniteStrainCrystalPLasticity: 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 && !line_search_update(rnorm_prev, dpk2))
    {
#ifdef DEBUG
      mooseWarning("FiniteStrainCrystalPLasticity: Failed with line search");
#endif
      _err_tol = true;
      return;
    }
 
    if (_use_line_search)
      rnorm = resid.L2norm();
 
    iter++;
  }
 
  if (iter >= _maxiter)
  {
#ifdef DEBUG
    mooseWarning("FiniteStrainCrystalPLasticity: Stress Integration error rmax = ", rnorm);
#endif
    _err_tol = true;
  }
}

Referenced by solveStatevar().

update_slip_system_resistance()

void FiniteStrainCrystalPlasticity::update_slip_system_resistance ( )

protectedvirtual

This function updates the slip system resistances.

Definition at line 826 of file FiniteStrainCrystalPlasticity.C.

{
  updateGss();
}

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

updateGss()

void FiniteStrainCrystalPlasticity::updateGss ( )

protectedvirtual

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 update_slip_system_resistance().

validParams()

InputParameters FiniteStrainCrystalPlasticity::validParams ( )

static

Definition at line 22 of file FiniteStrainCrystalPlasticity.C.

{
  InputParameters params = ComputeStressBase::validParams();
  params.addClassDescription("Deprecated class: please use CrystalPlasticityKalidindiUpdate and "
                             "ComputeMultipleCrystalPlasticityStress instead.  Crystal Plasticity "
                             "base class: FCC system with power law flow rule implemented");
  params.addRequiredParam<int>("nss", "Number of slip systems");
  params.addParam<std::vector<Real>>("gprops", {}, "Initial values of slip system resistances");
  params.addParam<std::vector<Real>>("hprops", {}, "Hardening properties");
  params.addParam<std::vector<Real>>("flowprops", {}, "Parameters used in slip rate equations");
  params.addRequiredParam<FileName>("slip_sys_file_name",
                                    "Name of the file containing the slip system");
  params.addParam<FileName>(
      "slip_sys_res_prop_file_name",
      "",
      "Name of the file containing the initial values of slip system resistances");
  params.addParam<FileName>(
      "slip_sys_flow_prop_file_name",
      "",
      "Name of the file containing the values of slip rate equation parameters");
  params.addParam<FileName>(
      "slip_sys_hard_prop_file_name",
      "",
      "Name of the file containing the values of hardness evolution parameters");
  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>("gtol", 1e2, "Constitutive slip system resistance residual tolerance");
  params.addParam<Real>("slip_incr_tol", 2e-2, "Maximum allowable slip in an increment");
  params.addParam<unsigned int>("maxiter", 100, "Maximum number of iterations for stress update");
  params.addParam<unsigned int>(
      "maxitergss", 100, "Maximum number of iterations for slip system resistance update");
  params.addParam<unsigned int>(
      "num_slip_sys_flowrate_props",
      2,
      "Number of flow rate properties for a slip system"); // Used for reading flow rate parameters
  params.addParam<UserObjectName>("read_prop_user_object",
                                  "The ElementReadPropertyFile "
                                  "GeneralUserObject to read element "
                                  "specific property values from file");
  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");
  MooseEnum intvar_read_options("slip_sys_file slip_sys_res_file none", "none");
  params.addParam<MooseEnum>(
      "intvar_read_type",
      intvar_read_options,
      "Read from options for initial value of internal variables: Default from .i file");
  params.addParam<unsigned int>("num_slip_sys_props",
                                0,
                                "Number of slip system specific properties provided in the file "
                                "containing slip system normals and directions");
  params.addParam<bool>(
      "gen_random_stress_flag",
      false,
      "Flag to generate random stress to perform time cutback on constitutive failure");
  params.addParam<bool>("input_random_scaling_var",
                        false,
                        "Flag to input scaling variable: _Cijkl(0,0,0,0) when false");
  params.addParam<Real>("random_scaling_var",
                        1e9,
                        "Random scaling variable: Large value can cause non-positive definiteness");
  params.addParam<unsigned int>(
      "random_seed",
      2000,
      "Random integer used to generate random stress when constitutive failure occurs");
  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");
 
  return params;
}

Referenced by FiniteStrainCPSlipRateRes::validParams().

Member Data Documentation

_a0

DenseVector<Real> FiniteStrainCrystalPlasticity::_a0

protected

Definition at line 322 of file FiniteStrainCrystalPlasticity.h.

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

_abs_tol

Real FiniteStrainCrystalPlasticity::_abs_tol

protected

Stress residual equation absolute tolerance.

Definition at line 247 of file FiniteStrainCrystalPlasticity.h.

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

_acc_slip

MaterialProperty<Real>& FiniteStrainCrystalPlasticity::_acc_slip

protected

Definition at line 307 of file FiniteStrainCrystalPlasticity.h.

Referenced by initQpStatefulProperties(), and postSolveStatevar().

_acc_slip_old

const MaterialProperty<Real>& FiniteStrainCrystalPlasticity::_acc_slip_old

protected

Definition at line 308 of file FiniteStrainCrystalPlasticity.h.

Referenced by preSolveStatevar().

_accslip_tmp

Real FiniteStrainCrystalPlasticity::_accslip_tmp

protected

Definition at line 336 of file FiniteStrainCrystalPlasticity.h.

Referenced by postSolveStatevar(), and updateGss().

_accslip_tmp_old

Real FiniteStrainCrystalPlasticity::_accslip_tmp_old

protected

Definition at line 336 of file FiniteStrainCrystalPlasticity.h.

Referenced by postSolveStatevar(), preSolveStatevar(), and 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

protected

Definition at line 317 of file FiniteStrainCrystalPlasticity.h.

Referenced by calc_schmid_tensor(), and postSolveQp().

_deformation_gradient

const MaterialProperty<RankTwoTensor>& FiniteStrainCrystalPlasticity::_deformation_gradient

protected

Definition at line 311 of file FiniteStrainCrystalPlasticity.h.

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

_deformation_gradient_old

const MaterialProperty<RankTwoTensor>& FiniteStrainCrystalPlasticity::_deformation_gradient_old

protected

Definition at line 312 of file FiniteStrainCrystalPlasticity.h.

Referenced by computeQpStress().

_delta_dfgrd

RankTwoTensor FiniteStrainCrystalPlasticity::_delta_dfgrd

protected

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 computeQpStress(), FiniteStrainCrystalPlasticity(), and preSolveQp().

_dfgrd_scale_factor

Real FiniteStrainCrystalPlasticity::_dfgrd_scale_factor

protected

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

Definition at line 351 of file FiniteStrainCrystalPlasticity.h.

Referenced by computeQpStress(), and preSolveQp().

_dfgrd_tmp

RankTwoTensor FiniteStrainCrystalPlasticity::_dfgrd_tmp

protected

Definition at line 330 of file FiniteStrainCrystalPlasticity.h.

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

_dfgrd_tmp_old

RankTwoTensor FiniteStrainCrystalPlasticity::_dfgrd_tmp_old

protected

Definition at line 349 of file FiniteStrainCrystalPlasticity.h.

Referenced by computeQpStress(), and preSolveQp().

_dgss_dsliprate

DenseMatrix<Real> FiniteStrainCrystalPlasticity::_dgss_dsliprate

protected

Definition at line 342 of file FiniteStrainCrystalPlasticity.h.

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

_dslipdtau

DenseVector<Real> FiniteStrainCrystalPlasticity::_dslipdtau

protected

Definition at line 332 of file FiniteStrainCrystalPlasticity.h.

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

_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

protected

Elasticity tensor material property.

Definition at line 316 of file FiniteStrainCrystalPlasticity.h.

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

_elasticity_tensor_name

const std::string FiniteStrainCrystalPlasticity::_elasticity_tensor_name

protected

Name of the elasticity tensor material property.

Definition at line 314 of file FiniteStrainCrystalPlasticity.h.

_err_tol

bool FiniteStrainCrystalPlasticity::_err_tol

protected

Definition at line 346 of file FiniteStrainCrystalPlasticity.h.

Referenced by calc_resid_jacob(), FiniteStrainCPSlipRateRes::calcResidJacobSlipRate(), calcResidual(), FiniteStrainCPSlipRateRes::calcResidualSlipRate(), computeQpStress(), FiniteStrainCrystalPlasticity(), FiniteStrainCPSlipRateRes::getSlipIncrements(), getSlipIncrements(), FiniteStrainCPSlipRateRes::lineSearchUpdateSlipRate(), postSolveQp(), preSolveQp(), solveQp(), FiniteStrainCPSlipRateRes::solveStatevar(), solveStatevar(), FiniteStrainCPSlipRateRes::solveStress(), and 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

protected

Definition at line 331 of file FiniteStrainCrystalPlasticity.h.

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

_first_step_iter

bool FiniteStrainCrystalPlasticity::_first_step_iter

protected

Flags to reset variables and reinitialize variables.

Definition at line 353 of file FiniteStrainCrystalPlasticity.h.

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

_first_substep

bool FiniteStrainCrystalPlasticity::_first_substep

protected

Definition at line 353 of file FiniteStrainCrystalPlasticity.h.

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

_flowprops

std::vector<Real> FiniteStrainCrystalPlasticity::_flowprops

protected

Definition at line 227 of file FiniteStrainCrystalPlasticity.h.

Referenced by getFlowRateParams().

_fp

MaterialProperty<RankTwoTensor>& FiniteStrainCrystalPlasticity::_fp

protected

Definition at line 299 of file FiniteStrainCrystalPlasticity.h.

Referenced by initQpStatefulProperties(), and postSolveStress().

_fp_inv

RankTwoTensor FiniteStrainCrystalPlasticity::_fp_inv

protected

Definition at line 331 of file FiniteStrainCrystalPlasticity.h.

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

_fp_old

const MaterialProperty<RankTwoTensor>& FiniteStrainCrystalPlasticity::_fp_old

protected

Definition at line 300 of file FiniteStrainCrystalPlasticity.h.

Referenced by preSolveStress().

_fp_old_inv

RankTwoTensor FiniteStrainCrystalPlasticity::_fp_old_inv

protected

Definition at line 331 of file FiniteStrainCrystalPlasticity.h.

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

_fp_prev_inv

RankTwoTensor FiniteStrainCrystalPlasticity::_fp_prev_inv

protected

Definition at line 331 of file FiniteStrainCrystalPlasticity.h.

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

_gen_rndm_stress_flag

bool FiniteStrainCrystalPlasticity::_gen_rndm_stress_flag

protected

Definition at line 270 of file FiniteStrainCrystalPlasticity.h.

Referenced by postSolveQp().

_gprops

std::vector<Real> FiniteStrainCrystalPlasticity::_gprops

protected

Definition at line 225 of file FiniteStrainCrystalPlasticity.h.

Referenced by getInitSlipSysRes().

_gss

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

protected

Definition at line 305 of file FiniteStrainCrystalPlasticity.h.

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

_gss_old

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

protected

Definition at line 306 of file FiniteStrainCrystalPlasticity.h.

Referenced by preSolveStatevar(), and updateGss().

_gss_tmp

std::vector<Real> FiniteStrainCrystalPlasticity::_gss_tmp

protected

Definition at line 337 of file FiniteStrainCrystalPlasticity.h.

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

_gss_tmp_old

std::vector<Real> FiniteStrainCrystalPlasticity::_gss_tmp_old

protected

Definition at line 338 of file FiniteStrainCrystalPlasticity.h.

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

_gtol

Real FiniteStrainCrystalPlasticity::_gtol

protected

Internal variable update equation tolerance.

Definition at line 249 of file FiniteStrainCrystalPlasticity.h.

Referenced by solveStatevar().

_h0

Real FiniteStrainCrystalPlasticity::_h0

protected

Definition at line 325 of file FiniteStrainCrystalPlasticity.h.

Referenced by getHardnessParams(), and updateGss().

_hprops

std::vector<Real> FiniteStrainCrystalPlasticity::_hprops

protected

Definition at line 226 of file FiniteStrainCrystalPlasticity.h.

Referenced by getHardnessParams(), and 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

protected

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

Definition at line 273 of file FiniteStrainCrystalPlasticity.h.

Referenced by postSolveQp().

_intvar_read_type

MooseEnum FiniteStrainCrystalPlasticity::_intvar_read_type

protected

Read from options for initial values of internal variables.

Definition at line 265 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity(), and initSlipSysProps().

_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(), postSolveQp(), preSolveQp(), and ComputeMultipleInelasticStressBase::updateQpStateSingleModel().

_lag_e

MaterialProperty<RankTwoTensor>& FiniteStrainCrystalPlasticity::_lag_e

protected

Definition at line 303 of file FiniteStrainCrystalPlasticity.h.

Referenced by initQpStatefulProperties(), and postSolveQp().

_lag_e_old

const MaterialProperty<RankTwoTensor>& FiniteStrainCrystalPlasticity::_lag_e_old

protected

Definition at line 304 of file FiniteStrainCrystalPlasticity.h.

_last_step_iter

bool FiniteStrainCrystalPlasticity::_last_step_iter

protected

Definition at line 353 of file FiniteStrainCrystalPlasticity.h.

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

_lsrch_max_iter

unsigned int FiniteStrainCrystalPlasticity::_lsrch_max_iter

protected

Line search bisection method maximum iteration number.

Definition at line 294 of file FiniteStrainCrystalPlasticity.h.

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

_lsrch_method

MooseEnum FiniteStrainCrystalPlasticity::_lsrch_method

protected

Definition at line 297 of file FiniteStrainCrystalPlasticity.h.

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

_lsrch_tol

Real FiniteStrainCrystalPlasticity::_lsrch_tol

protected

Line search bisection method tolerance.

Definition at line 291 of file FiniteStrainCrystalPlasticity.h.

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

_max_substep_iter

unsigned int FiniteStrainCrystalPlasticity::_max_substep_iter

protected

Maximum number of substep iterations.

Definition at line 282 of file FiniteStrainCrystalPlasticity.h.

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

_maxiter

unsigned int FiniteStrainCrystalPlasticity::_maxiter

protected

Maximum number of iterations for stress update.

Definition at line 254 of file FiniteStrainCrystalPlasticity.h.

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

_maxiterg

unsigned int FiniteStrainCrystalPlasticity::_maxiterg

protected

Maximum number of iterations for internal variable update.

Definition at line 256 of file FiniteStrainCrystalPlasticity.h.

Referenced by solveStatevar().

_mechanical_strain

const MaterialProperty<RankTwoTensor>& ComputeGeneralStressBase::_mechanical_strain

protectedinherited

Mechanical strain material property.

Definition at line 46 of file ComputeGeneralStressBase.h.

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

_min_lsrch_step

Real FiniteStrainCrystalPlasticity::_min_lsrch_step

protected

Minimum line search step size.

Definition at line 288 of file FiniteStrainCrystalPlasticity.h.

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

_mo

DenseVector<Real> FiniteStrainCrystalPlasticity::_mo

protected

Definition at line 319 of file FiniteStrainCrystalPlasticity.h.

Referenced by calc_schmid_tensor(), and getSlipSystems().

_no

DenseVector<Real> FiniteStrainCrystalPlasticity::_no

protected

Definition at line 320 of file FiniteStrainCrystalPlasticity.h.

Referenced by calc_schmid_tensor(), and getSlipSystems().

_nss

const unsigned int FiniteStrainCrystalPlasticity::_nss

protected

Number of slip system resistance.

Definition at line 223 of file FiniteStrainCrystalPlasticity.h.

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

_num_slip_sys_flowrate_props

unsigned int FiniteStrainCrystalPlasticity::_num_slip_sys_flowrate_props

protected

Number of slip system flow rate parameters.

Definition at line 259 of file FiniteStrainCrystalPlasticity.h.

Referenced by getFlowRateParams(), and readFileFlowRateParams().

_num_slip_sys_props

unsigned int FiniteStrainCrystalPlasticity::_num_slip_sys_props

protected

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(), and getSlipSystems().

_pk2

MaterialProperty<RankTwoTensor>& FiniteStrainCrystalPlasticity::_pk2

protected

Definition at line 301 of file FiniteStrainCrystalPlasticity.h.

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

_pk2_old

const MaterialProperty<RankTwoTensor>& FiniteStrainCrystalPlasticity::_pk2_old

protected

Definition at line 302 of file FiniteStrainCrystalPlasticity.h.

Referenced by preSolveStress().

_pk2_tmp

RankTwoTensor FiniteStrainCrystalPlasticity::_pk2_tmp

protected

Definition at line 335 of file FiniteStrainCrystalPlasticity.h.

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

_pk2_tmp_old

RankTwoTensor FiniteStrainCrystalPlasticity::_pk2_tmp_old

protected

Definition at line 335 of file FiniteStrainCrystalPlasticity.h.

Referenced by postSolveStress(), and preSolveStress().

_r

Real FiniteStrainCrystalPlasticity::_r

protected

Definition at line 328 of file FiniteStrainCrystalPlasticity.h.

Referenced by getHardnessParams(), and updateGss().

_read_from_slip_sys_file

bool FiniteStrainCrystalPlasticity::_read_from_slip_sys_file

protected

Definition at line 344 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity(), and getSlipSystems().

_rndm_scale_var

Real FiniteStrainCrystalPlasticity::_rndm_scale_var

protected

Scaling value.

Definition at line 276 of file FiniteStrainCrystalPlasticity.h.

Referenced by postSolveQp().

_rndm_seed

unsigned int FiniteStrainCrystalPlasticity::_rndm_seed

protected

Seed value.

Definition at line 279 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity().

_rtol

Real FiniteStrainCrystalPlasticity::_rtol

protected

Stress residual equation relative tolerance.

Definition at line 245 of file FiniteStrainCrystalPlasticity.h.

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

_s0

std::vector<RankTwoTensor> FiniteStrainCrystalPlasticity::_s0

protected

Definition at line 333 of file FiniteStrainCrystalPlasticity.h.

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

_slip_incr

DenseVector<Real> FiniteStrainCrystalPlasticity::_slip_incr

protected

Definition at line 332 of file FiniteStrainCrystalPlasticity.h.

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

_slip_incr_tol

Real FiniteStrainCrystalPlasticity::_slip_incr_tol

protected

Slip increment tolerance.

Definition at line 251 of file FiniteStrainCrystalPlasticity.h.

Referenced by getSlipIncrements().

_slip_sys_file_name

std::string FiniteStrainCrystalPlasticity::_slip_sys_file_name

protected

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

Definition at line 230 of file FiniteStrainCrystalPlasticity.h.

Referenced by getSlipSystems().

_slip_sys_flow_prop_file_name

std::string FiniteStrainCrystalPlasticity::_slip_sys_flow_prop_file_name

protected

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 initSlipSysProps(), and readFileFlowRateParams().

_slip_sys_hard_prop_file_name

std::string FiniteStrainCrystalPlasticity::_slip_sys_hard_prop_file_name

protected

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 initSlipSysProps().

_slip_sys_props

DenseVector<Real> FiniteStrainCrystalPlasticity::_slip_sys_props

protected

Definition at line 340 of file FiniteStrainCrystalPlasticity.h.

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

_slip_sys_res_prop_file_name

std::string FiniteStrainCrystalPlasticity::_slip_sys_res_prop_file_name

protected

File should contain initial values of the slip system resistances.

Definition at line 233 of file FiniteStrainCrystalPlasticity.h.

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

_tan_mod_type

MooseEnum FiniteStrainCrystalPlasticity::_tan_mod_type

protected

Type of tangent moduli calculation.

Definition at line 262 of file FiniteStrainCrystalPlasticity.h.

Referenced by calcTangentModuli().

_tau

DenseVector<Real> FiniteStrainCrystalPlasticity::_tau

protected

Definition at line 332 of file FiniteStrainCrystalPlasticity.h.

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

_tau_init

Real FiniteStrainCrystalPlasticity::_tau_init

protected

Definition at line 327 of file FiniteStrainCrystalPlasticity.h.

Referenced by getHardnessParams().

_tau_sat

Real FiniteStrainCrystalPlasticity::_tau_sat

protected

Definition at line 326 of file FiniteStrainCrystalPlasticity.h.

Referenced by getHardnessParams(), and updateGss().

_update_rot

MaterialProperty<RankTwoTensor>& FiniteStrainCrystalPlasticity::_update_rot

protected

Definition at line 309 of file FiniteStrainCrystalPlasticity.h.

Referenced by initQpStatefulProperties(), and postSolveQp().

_use_line_search

bool FiniteStrainCrystalPlasticity::_use_line_search

protected

Flag to activate line serach.

Definition at line 285 of file FiniteStrainCrystalPlasticity.h.

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

_xm

DenseVector<Real> FiniteStrainCrystalPlasticity::_xm

protected

Definition at line 323 of file FiniteStrainCrystalPlasticity.h.

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

---

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

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