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:

[legend]

Public Member Functions
	FiniteStrainCrystalPlasticity (const InputParameters &parameters)

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

virtual void	computeQpElasticityTensor ()
	This function updates the elasticity tensor at a quadrature point. More...

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

virtual void	calc_resid_jacob (RankTwoTensor &, RankFourTensor &)
	This function calls the residual and jacobian functions used in the stress update algorithm. More...

virtual void	getSlipIncrements ()
	This function updates the slip increments. More...

virtual void	update_slip_system_resistance ()
	This function updates the slip system resistances. More...

virtual void	updateGss ()
	This function updates the slip system resistances. More...

virtual void	getSlipSystems ()
	This function reads slip system from file - see test. More...

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

virtual void	readFileInitSlipSysRes ()
	This function read slip system resistances from file - see test. More...

virtual void	getInitSlipSysRes ()
	This function assign slip system resistances - see test. More...

virtual void	readFileFlowRateParams ()
	This function read flow rate parameters from file - see test. More...

virtual void	getFlowRateParams ()
	This function assign flow rate parameters - see test. More...

virtual void	readFileHardnessParams ()
	This function read hardness parameters from file. More...

virtual void	getHardnessParams ()
	This function assign flow rate parameters from .i file - see test. More...

virtual void	initSlipSysProps ()
	This function initializes slip system resistances. More...

virtual void	initAdditionalProps ()
	This function initializes additional parameters. More...

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

virtual void	solveQp ()
	This function solves stress and internal variables. More...

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

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

virtual void	solveStatevar ()
	This function solves internal variables. More...

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

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

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

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

virtual void	calcResidual (RankTwoTensor &)
	This function calculate stress residual. More...

virtual void	calcJacobian (RankFourTensor &)
	This function calculate jacobian. More...

virtual RankFourTensor	calcTangentModuli ()
	This function calculate the tangent moduli for preconditioner. More...

virtual RankFourTensor	elasticTangentModuli ()
	This function calculate the elastic tangent moduli for preconditioner. More...

virtual RankFourTensor	elastoPlasticTangentModuli ()
	This function calculate the exact tangent moduli for preconditioner. More...

RankTwoTensor	get_current_rotation (const RankTwoTensor &a)
	This function perform RU decomposition to obtain the rotation tensor. More...

RankTwoTensor	getMatRot (const RankTwoTensor &a)
	This function perform RU decomposition to obtain the rotation tensor. More...

void	calc_schmid_tensor ()
	This function calculate the Schmid tensor. More...

bool	line_search_update (const Real rnorm_prev, const RankTwoTensor)
	This function performs the line search update. More...

void	internalVariableUpdateNRiteration ()
	This function updates internal variables after each NewTon Raphson iteration (_fp_inv) More...

virtual void	computeQpProperties () override

Protected Attributes
const unsigned int	_nss
	Number of slip system resistance. More...

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

std::string	_slip_sys_res_prop_file_name
	File should contain initial values of the slip system resistances. More...

std::string	_slip_sys_flow_prop_file_name
	File should contain values of the flow rate equation parameters. More...

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

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

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

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

Real	_slip_incr_tol
	Slip increment tolerance. More...

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

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

unsigned int	_num_slip_sys_flowrate_props
	Number of slip system flow rate parameters. More...

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

MooseEnum	_intvar_read_type
	Read from options for initial values of internal variables. More...

unsigned int	_num_slip_sys_props
	Number of slip system specific properties provided in the file containing slip system normals and directions. More...

bool	_gen_rndm_stress_flag

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

Real	_rndm_scale_var
	Scaling value. More...

unsigned int	_rndm_seed
	Seed value. More...

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

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

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

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

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

MooseEnum	_lsrch_method

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

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

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

RankTwoTensor	_dfgrd_tmp_old

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

bool	_first_step_iter
	Flags to reset variables and reinitialize variables. More...

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

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

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

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

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

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

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

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 27 of file FiniteStrainCrystalPlasticity.h.

Constructor & Destructor Documentation

◆ FiniteStrainCrystalPlasticity()

FiniteStrainCrystalPlasticity::FiniteStrainCrystalPlasticity ( const InputParameters & parameters )

Definition at line 100 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 alip 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 232 of file FiniteStrainCrystalPlasticity.C.

Referenced by initSlipSysProps().

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

◆ 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 879 of file FiniteStrainCrystalPlasticity.C.

Referenced by solveStress().

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

◆ calc_schmid_tensor()

void FiniteStrainCrystalPlasticity::calc_schmid_tensor ( )

protected

This function calculate the Schmid tensor.

Definition at line 1047 of file FiniteStrainCrystalPlasticity.C.

Referenced by preSolveQp().

 {
   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 (unsigned int j = 0; j < LIBMESH_DIM; ++j)
     {
       mo(i * LIBMESH_DIM + j) = 0.0;
       for (unsigned int k = 0; k < LIBMESH_DIM; ++k)
         mo(i * LIBMESH_DIM + j) =
             mo(i * LIBMESH_DIM + j) + _crysrot[_qp](j, k) * _mo(i * LIBMESH_DIM + k);
     }
 
     for (unsigned int j = 0; j < LIBMESH_DIM; ++j)
     {
       no(i * LIBMESH_DIM + j) = 0.0;
       for (unsigned int k = 0; k < LIBMESH_DIM; ++k)
         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 (unsigned int j = 0; j < LIBMESH_DIM; ++j)
       for (unsigned int k = 0; k < LIBMESH_DIM; ++k)
         _s0[i](j, k) = mo(i * LIBMESH_DIM + j) * no(i * LIBMESH_DIM + k);
 }

◆ calcJacobian()

void FiniteStrainCrystalPlasticity::calcJacobian ( RankFourTensor & jac )

protectedvirtual

This function calculate jacobian.

Definition at line 925 of file FiniteStrainCrystalPlasticity.C.

Referenced by calc_resid_jacob().

 {
   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 (unsigned int i = 0; i < LIBMESH_DIM; ++i)
     for (unsigned int j = 0; j < LIBMESH_DIM; ++j)
       for (unsigned int k = 0; k < LIBMESH_DIM; ++k)
         dfedfpinv(i, j, k, j) = _dfgrd_tmp(i, k);
 
   for (unsigned int i = 0; i < LIBMESH_DIM; ++i)
     for (unsigned int j = 0; j < LIBMESH_DIM; ++j)
       for (unsigned int k = 0; k < LIBMESH_DIM; ++k)
       {
         deedfe(i, j, k, i) = deedfe(i, j, k, i) + _fe(k, j) * 0.5;
         deedfe(i, j, k, j) = deedfe(i, j, k, j) + _fe(k, i) * 0.5;
       }
 
   for (unsigned int i = 0; i < _nss; ++i)
     dfpinvdpk2 += (dfpinvdslip[i] * _dslipdtau(i)).outerProduct(dtaudpk2[i]);
 
   jac =
       RankFourTensor::IdentityFour() - (_elasticity_tensor[_qp] * deedfe * dfedfpinv * dfpinvdpk2);
 }

◆ calcResidual()

void FiniteStrainCrystalPlasticity::calcResidual ( RankTwoTensor & resid )

protectedvirtual

This function calculate stress residual.

Definition at line 888 of file FiniteStrainCrystalPlasticity.C.

Referenced by calc_resid_jacob(), and line_search_update().

 {
   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;
 }

◆ 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 1030 of file FiniteStrainCrystalPlasticity.C.

Referenced by postSolveQp().

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

◆ computeQpElasticityTensor()

void FiniteStrainCrystalPlasticity::computeQpElasticityTensor ( )

protectedvirtual

This function updates the elasticity tensor at a quadrature point.

Presently void.

Definition at line 1025 of file FiniteStrainCrystalPlasticity.C.

1026 {

1027 }

◆ computeQpProperties()

void ComputeStressBase::computeQpProperties ( )

overrideprotectedvirtualinherited

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

Definition at line 491 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 1118 of file FiniteStrainCrystalPlasticity.C.

Referenced by calcTangentModuli().

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

◆ elastoPlasticTangentModuli()

RankFourTensor FiniteStrainCrystalPlasticity::elastoPlasticTangentModuli ( )

protectedvirtual

This function calculate the exact tangent moduli for preconditioner.

Definition at line 1079 of file FiniteStrainCrystalPlasticity.C.

Referenced by calcTangentModuli().

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

◆ 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 983 of file FiniteStrainCrystalPlasticity.C.

Referenced by postSolveQp().

 {
   return getMatRot(a);
 }

◆ 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 342 of file FiniteStrainCrystalPlasticity.C.

Referenced by initSlipSysProps().

 {
   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;
     }
   }
 }

◆ getHardnessParams()

void FiniteStrainCrystalPlasticity::getHardnessParams ( )

protectedvirtual

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

Definition at line 410 of file FiniteStrainCrystalPlasticity.C.

Referenced by initSlipSysProps().

 {
   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];
 }

◆ 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 260 of file FiniteStrainCrystalPlasticity.C.

Referenced by initSlipSysProps().

 {
   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");
 }

◆ getMatRot()

RankTwoTensor FiniteStrainCrystalPlasticity::getMatRot ( const RankTwoTensor & a )

protected

This function perform RU decomposition to obtain the rotation tensor.

Definition at line 990 of file FiniteStrainCrystalPlasticity.C.

Referenced by get_current_rotation().

 {
   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 (unsigned int i = 0; i < LIBMESH_DIM; ++i)
     for (unsigned int j = 0; j < LIBMESH_DIM; ++j)
       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 (unsigned int i = 0; i < LIBMESH_DIM; ++i)
     diag(i, i) = std::sqrt(w[i]);
 
   for (unsigned int i = 0; i < LIBMESH_DIM; ++i)
     for (unsigned int j = 0; j < LIBMESH_DIM; ++j)
       evec(i, j) = cmat[i][j];
 
   rot = a * ((evec.transpose() * diag * evec).inverse());
 
   return rot;
 }

◆ 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 959 of file FiniteStrainCrystalPlasticity.C.

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

 {
   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;
 }

◆ getSlipSystems()

void FiniteStrainCrystalPlasticity::getSlipSystems ( )

protectedvirtual

This function reads slip system from file - see test.

Definition at line 424 of file FiniteStrainCrystalPlasticity.C.

Referenced by FiniteStrainCrystalPlasticity().

 {
   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 (unsigned int j = 0; j < LIBMESH_DIM; ++j)
       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 (unsigned int j = 0; j < LIBMESH_DIM; ++j)
       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 (unsigned int j = 0; j < LIBMESH_DIM; ++j)
       _mo(i * LIBMESH_DIM + j) = vec[j] / mag;
 
     mag = 0.0;
     for (unsigned int j = 0; j < LIBMESH_DIM; ++j)
       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();
 }

◆ initAdditionalProps()

void FiniteStrainCrystalPlasticity::initAdditionalProps ( )

protectedvirtual

This function initializes additional parameters.

Definition at line 482 of file FiniteStrainCrystalPlasticity.C.

Referenced by initQpStatefulProperties().

483 {

484 }

◆ initQpStatefulProperties()

void FiniteStrainCrystalPlasticity::initQpStatefulProperties ( )

protectedvirtual

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

Reimplemented from ComputeStressBase.

Definition at line 189 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 206 of file FiniteStrainCrystalPlasticity.C.

Referenced by initQpStatefulProperties().

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

◆ internalVariableUpdateNRiteration()

void FiniteStrainCrystalPlasticity::internalVariableUpdateNRiteration ( )

protected

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

Definition at line 1207 of file FiniteStrainCrystalPlasticity.C.

Referenced by solveStress().

 {
   _fp_prev_inv = _fp_inv; // update _fp_prev_inv
 }

◆ 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 1124 of file FiniteStrainCrystalPlasticity.C.

Referenced by solveStress().

 {
   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;
   }
 }

◆ postSolveQp()

void FiniteStrainCrystalPlasticity::postSolveQp ( )

protectedvirtual

This function update stress and internal variable after solve.

Definition at line 586 of file FiniteStrainCrystalPlasticity.C.

Referenced by computeQpStress().

 {
   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];
   }
 }

◆ postSolveStatevar()

void FiniteStrainCrystalPlasticity::postSolveStatevar ( )

protectedvirtual

This function update internal variable after solve.

Definition at line 681 of file FiniteStrainCrystalPlasticity.C.

Referenced by solveQp().

 {
   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;
     }
   }
 }

◆ postSolveStress()

void FiniteStrainCrystalPlasticity::postSolveStress ( )

protectedvirtual

This function update stress and plastic deformation gradient after solve.

Definition at line 800 of file FiniteStrainCrystalPlasticity.C.

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

 {
   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;
     }
   }
 }

◆ preSolveQp()

void FiniteStrainCrystalPlasticity::preSolveQp ( )

protectedvirtual

This function set variables for stress and internal variable solve.

Definition at line 558 of file FiniteStrainCrystalPlasticity.C.

Referenced by computeQpStress().

 {
   // 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;
 }

◆ preSolveStatevar()

void FiniteStrainCrystalPlasticity::preSolveStatevar ( )

protectedvirtual

This function set variables for internal variable solve.

Definition at line 619 of file FiniteStrainCrystalPlasticity.C.

Referenced by solveQp().

 {
   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;
   }
 }

◆ preSolveStress()

void FiniteStrainCrystalPlasticity::preSolveStress ( )

protectedvirtual

This function set variables for stress solve.

Reimplemented in FiniteStrainCPSlipRateRes.

Definition at line 704 of file FiniteStrainCrystalPlasticity.C.

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

 {
   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;
   }
 }

◆ readFileFlowRateParams()

void FiniteStrainCrystalPlasticity::readFileFlowRateParams ( )

protectedvirtual

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

Definition at line 313 of file FiniteStrainCrystalPlasticity.C.

Referenced by initSlipSysProps().

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

◆ readFileHardnessParams()

void FiniteStrainCrystalPlasticity::readFileHardnessParams ( )

protectedvirtual

This function read hardness parameters from file.

Definition at line 404 of file FiniteStrainCrystalPlasticity.C.

Referenced by initSlipSysProps().

405 {

406 }

◆ readFileInitSlipSysRes()

void FiniteStrainCrystalPlasticity::readFileInitSlipSysRes ( )

protectedvirtual

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

Definition at line 242 of file FiniteStrainCrystalPlasticity.C.

Referenced by initSlipSysProps().

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

◆ solveQp()

void FiniteStrainCrystalPlasticity::solveQp ( )

protectedvirtual

This function solves stress and internal variables.

Definition at line 576 of file FiniteStrainCrystalPlasticity.C.

Referenced by computeQpStress().

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

◆ solveStatevar()

void FiniteStrainCrystalPlasticity::solveStatevar ( )

protectedvirtual

This function solves internal variables.

Reimplemented in FiniteStrainCPSlipRateRes.

Definition at line 639 of file FiniteStrainCrystalPlasticity.C.

Referenced by solveQp().

 {
   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;
   }
 }

◆ solveStress()

void FiniteStrainCrystalPlasticity::solveStress ( )

protectedvirtual

This function solves for stress, updates plastic deformation gradient.

Reimplemented in FiniteStrainCPSlipRateRes.

Definition at line 729 of file FiniteStrainCrystalPlasticity.C.

Referenced by solveStatevar().

 {
   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;
   }
 }

◆ update_slip_system_resistance()

void FiniteStrainCrystalPlasticity::update_slip_system_resistance ( )

protectedvirtual

This function updates the slip system resistances.

Definition at line 824 of file FiniteStrainCrystalPlasticity.C.

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

 {
   updateGss();
 }

◆ 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 834 of file FiniteStrainCrystalPlasticity.C.

Referenced by update_slip_system_resistance().

 {
   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;
     }
   }
 }

Member Data Documentation

◆ _a0

DenseVector<Real> FiniteStrainCrystalPlasticity::_a0

protected

Definition at line 326 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 251 of file FiniteStrainCrystalPlasticity.h.

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

◆ _acc_slip

MaterialProperty<Real>& FiniteStrainCrystalPlasticity::_acc_slip

protected

Definition at line 311 of file FiniteStrainCrystalPlasticity.h.

Referenced by initQpStatefulProperties(), and postSolveStatevar().

◆ _acc_slip_old

const MaterialProperty<Real>& FiniteStrainCrystalPlasticity::_acc_slip_old

protected

Definition at line 312 of file FiniteStrainCrystalPlasticity.h.

Referenced by preSolveStatevar().

◆ _accslip_tmp

Real FiniteStrainCrystalPlasticity::_accslip_tmp

protected

Definition at line 340 of file FiniteStrainCrystalPlasticity.h.

Referenced by postSolveStatevar(), and updateGss().

◆ _accslip_tmp_old

Real FiniteStrainCrystalPlasticity::_accslip_tmp_old

protected

Definition at line 340 of file FiniteStrainCrystalPlasticity.h.

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

◆ _base_name

const std::string ComputeStressBase::_base_name

protectedinherited

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

Definition at line 44 of file ComputeStressBase.h.

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

◆ _crysrot

const MaterialProperty<RankTwoTensor>& FiniteStrainCrystalPlasticity::_crysrot

protected

Definition at line 321 of file FiniteStrainCrystalPlasticity.h.

Referenced by calc_schmid_tensor(), and postSolveQp().

◆ _deformation_gradient

const MaterialProperty<RankTwoTensor>& FiniteStrainCrystalPlasticity::_deformation_gradient

protected

Definition at line 315 of file FiniteStrainCrystalPlasticity.h.

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

◆ _deformation_gradient_old

const MaterialProperty<RankTwoTensor>& FiniteStrainCrystalPlasticity::_deformation_gradient_old

protected

Definition at line 316 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 353 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 355 of file FiniteStrainCrystalPlasticity.h.

Referenced by computeQpStress(), and preSolveQp().

◆ _dfgrd_tmp

RankTwoTensor FiniteStrainCrystalPlasticity::_dfgrd_tmp

protected

Definition at line 334 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 353 of file FiniteStrainCrystalPlasticity.h.

Referenced by computeQpStress(), and preSolveQp().

◆ _dgss_dsliprate

DenseMatrix<Real> FiniteStrainCrystalPlasticity::_dgss_dsliprate

protected

Definition at line 346 of file FiniteStrainCrystalPlasticity.h.

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

◆ _dslipdtau

DenseVector<Real> FiniteStrainCrystalPlasticity::_dslipdtau

protected

Definition at line 336 of file FiniteStrainCrystalPlasticity.h.

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

◆ _elastic_strain

MaterialProperty<RankTwoTensor>& ComputeStressBase::_elastic_strain

protectedinherited

Elastic strain material property.

Definition at line 51 of file ComputeStressBase.h.

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

◆ _elasticity_tensor

const MaterialProperty<RankFourTensor>& FiniteStrainCrystalPlasticity::_elasticity_tensor

protected

Elasticity tensor material property.

Definition at line 320 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 318 of file FiniteStrainCrystalPlasticity.h.

◆ _err_tol

bool FiniteStrainCrystalPlasticity::_err_tol

protected

Definition at line 350 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>& ComputeStressBase::_extra_stress

protectedinherited

Extra stress tensor.

Definition at line 54 of file ComputeStressBase.h.

Referenced by ComputeStressBase::computeQpProperties().

◆ _fe

RankTwoTensor FiniteStrainCrystalPlasticity::_fe

protected

Definition at line 335 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 357 of file FiniteStrainCrystalPlasticity.h.

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

◆ _first_substep

bool FiniteStrainCrystalPlasticity::_first_substep

protected

Definition at line 357 of file FiniteStrainCrystalPlasticity.h.

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

◆ _flowprops

std::vector<Real> FiniteStrainCrystalPlasticity::_flowprops

protected

Definition at line 231 of file FiniteStrainCrystalPlasticity.h.

Referenced by getFlowRateParams().

◆ _fp

MaterialProperty<RankTwoTensor>& FiniteStrainCrystalPlasticity::_fp

protected

Definition at line 303 of file FiniteStrainCrystalPlasticity.h.

Referenced by initQpStatefulProperties(), and postSolveStress().

◆ _fp_inv

RankTwoTensor FiniteStrainCrystalPlasticity::_fp_inv

protected

Definition at line 335 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 304 of file FiniteStrainCrystalPlasticity.h.

Referenced by preSolveStress().

◆ _fp_old_inv

RankTwoTensor FiniteStrainCrystalPlasticity::_fp_old_inv

protected

Definition at line 335 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 335 of file FiniteStrainCrystalPlasticity.h.

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

◆ _gen_rndm_stress_flag

bool FiniteStrainCrystalPlasticity::_gen_rndm_stress_flag

protected

Definition at line 274 of file FiniteStrainCrystalPlasticity.h.

Referenced by postSolveQp().

◆ _gprops

std::vector<Real> FiniteStrainCrystalPlasticity::_gprops

protected

Definition at line 229 of file FiniteStrainCrystalPlasticity.h.

Referenced by getInitSlipSysRes().

◆ _gss

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

protected

Definition at line 309 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 310 of file FiniteStrainCrystalPlasticity.h.

Referenced by preSolveStatevar(), and updateGss().

◆ _gss_tmp

std::vector<Real> FiniteStrainCrystalPlasticity::_gss_tmp

protected

Definition at line 341 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 342 of file FiniteStrainCrystalPlasticity.h.

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

◆ _gtol

Real FiniteStrainCrystalPlasticity::_gtol

protected

Internal variable update equation tolerance.

Definition at line 253 of file FiniteStrainCrystalPlasticity.h.

Referenced by solveStatevar().

◆ _h0

Real FiniteStrainCrystalPlasticity::_h0

protected

Definition at line 329 of file FiniteStrainCrystalPlasticity.h.

Referenced by getHardnessParams(), and updateGss().

◆ _hprops

std::vector<Real> FiniteStrainCrystalPlasticity::_hprops

protected

Definition at line 230 of file FiniteStrainCrystalPlasticity.h.

Referenced by getHardnessParams(), and updateGss().

◆ _initial_stress_fcn

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

protectedinherited

initial stress components

Definition at line 57 of file ComputeStressBase.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 277 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 269 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity(), and initSlipSysProps().

◆ _Jacobian_mult

MaterialProperty<RankFourTensor>& ComputeStressBase::_Jacobian_mult

protectedinherited

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

Definition at line 60 of file ComputeStressBase.h.

◆ _lag_e

MaterialProperty<RankTwoTensor>& FiniteStrainCrystalPlasticity::_lag_e

protected

Definition at line 307 of file FiniteStrainCrystalPlasticity.h.

Referenced by initQpStatefulProperties(), and postSolveQp().

◆ _lag_e_old

const MaterialProperty<RankTwoTensor>& FiniteStrainCrystalPlasticity::_lag_e_old

protected

Definition at line 308 of file FiniteStrainCrystalPlasticity.h.

◆ _last_step_iter

bool FiniteStrainCrystalPlasticity::_last_step_iter

protected

Definition at line 357 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 298 of file FiniteStrainCrystalPlasticity.h.

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

◆ _lsrch_method

MooseEnum FiniteStrainCrystalPlasticity::_lsrch_method

protected

Definition at line 301 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 295 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 286 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 258 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 260 of file FiniteStrainCrystalPlasticity.h.

Referenced by solveStatevar().

◆ _mechanical_strain

const MaterialProperty<RankTwoTensor>& ComputeStressBase::_mechanical_strain

protectedinherited

Mechanical strain material property.

Definition at line 47 of file ComputeStressBase.h.

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

◆ _min_lsrch_step

Real FiniteStrainCrystalPlasticity::_min_lsrch_step

protected

Minimum line search step size.

Definition at line 292 of file FiniteStrainCrystalPlasticity.h.

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

◆ _mo

DenseVector<Real> FiniteStrainCrystalPlasticity::_mo

protected

Definition at line 323 of file FiniteStrainCrystalPlasticity.h.

Referenced by calc_schmid_tensor(), and getSlipSystems().

◆ _no

DenseVector<Real> FiniteStrainCrystalPlasticity::_no

protected

Definition at line 324 of file FiniteStrainCrystalPlasticity.h.

Referenced by calc_schmid_tensor(), and getSlipSystems().

◆ _nss

const unsigned int FiniteStrainCrystalPlasticity::_nss

protected

◆ _num_slip_sys_flowrate_props

unsigned int FiniteStrainCrystalPlasticity::_num_slip_sys_flowrate_props

protected

Number of slip system flow rate parameters.

Definition at line 263 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 272 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity(), and getSlipSystems().

◆ _pk2

MaterialProperty<RankTwoTensor>& FiniteStrainCrystalPlasticity::_pk2

protected

Definition at line 305 of file FiniteStrainCrystalPlasticity.h.

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

◆ _pk2_old

const MaterialProperty<RankTwoTensor>& FiniteStrainCrystalPlasticity::_pk2_old

protected

Definition at line 306 of file FiniteStrainCrystalPlasticity.h.

Referenced by preSolveStress().

◆ _pk2_tmp

RankTwoTensor FiniteStrainCrystalPlasticity::_pk2_tmp

protected

Definition at line 339 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 339 of file FiniteStrainCrystalPlasticity.h.

Referenced by postSolveStress(), and preSolveStress().

◆ _r

Real FiniteStrainCrystalPlasticity::_r

protected

Definition at line 332 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 348 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity(), and getSlipSystems().

◆ _rndm_scale_var

Real FiniteStrainCrystalPlasticity::_rndm_scale_var

protected

Scaling value.

Definition at line 280 of file FiniteStrainCrystalPlasticity.h.

Referenced by postSolveQp().

◆ _rndm_seed

unsigned int FiniteStrainCrystalPlasticity::_rndm_seed

protected

Seed value.

Definition at line 283 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity().

◆ _rtol

Real FiniteStrainCrystalPlasticity::_rtol

protected

Stress residual equation relative tolerance.

Definition at line 249 of file FiniteStrainCrystalPlasticity.h.

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

◆ _s0

std::vector<RankTwoTensor> FiniteStrainCrystalPlasticity::_s0

protected

Definition at line 337 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 336 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 255 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 234 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 243 of file FiniteStrainCrystalPlasticity.h.

Referenced by initSlipSysProps(), and readFileFlowRateParams().