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)

Static Public Member Functions
static InputParameters	validParams ()

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< const Function * >	_initial_stress_fcn
	initial stress components More...

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

Detailed Description

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

Constructor & Destructor Documentation

◆ FiniteStrainCrystalPlasticity()

FiniteStrainCrystalPlasticity::FiniteStrainCrystalPlasticity ( const InputParameters & parameters )

Definition at line 101 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 233 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 880 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 1048 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 (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);
 }

Referenced by preSolveQp().

◆ calcJacobian()

void FiniteStrainCrystalPlasticity::calcJacobian ( RankFourTensor & jac )

protectedvirtual

This function calculate jacobian.

Definition at line 926 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 (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);
 }

Referenced by calc_resid_jacob().

◆ calcResidual()

void FiniteStrainCrystalPlasticity::calcResidual ( RankTwoTensor & resid )

protectedvirtual

This function calculate stress residual.

Definition at line 889 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 1031 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 1026 of file FiniteStrainCrystalPlasticity.C.

1027 {

1028 }

◆ computeQpProperties()

void ComputeStressBase::computeQpProperties ( )

overrideprotectedvirtualinherited

Definition at line 50 of file ComputeStressBase.C.

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

◆ computeQpStress()

void 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 492 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 1119 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 1080 of file FiniteStrainCrystalPlasticity.C.

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

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 984 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 343 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 411 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 261 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 991 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 (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;
 }

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

Referenced by FiniteStrainCrystalPlasticity().

◆ initAdditionalProps()

void FiniteStrainCrystalPlasticity::initAdditionalProps ( )

protectedvirtual

This function initializes additional parameters.

Definition at line 483 of file FiniteStrainCrystalPlasticity.C.

484 {

485 }

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

Definition at line 190 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 207 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 1208 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 1125 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 587 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 682 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 801 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 559 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 620 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 705 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 314 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 405 of file FiniteStrainCrystalPlasticity.C.

406 {

407 }

Referenced by initSlipSysProps().

◆ readFileInitSlipSysRes()

void FiniteStrainCrystalPlasticity::readFileInitSlipSysRes ( )

protectedvirtual

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

Definition at line 243 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 577 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 640 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 730 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 825 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 835 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(
       "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 327 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 252 of file FiniteStrainCrystalPlasticity.h.

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

◆ _acc_slip

MaterialProperty<Real>& FiniteStrainCrystalPlasticity::_acc_slip

protected

Definition at line 312 of file FiniteStrainCrystalPlasticity.h.

Referenced by initQpStatefulProperties(), and postSolveStatevar().

◆ _acc_slip_old

const MaterialProperty<Real>& FiniteStrainCrystalPlasticity::_acc_slip_old

protected

Definition at line 313 of file FiniteStrainCrystalPlasticity.h.

Referenced by preSolveStatevar().

◆ _accslip_tmp

Real FiniteStrainCrystalPlasticity::_accslip_tmp

protected

Definition at line 341 of file FiniteStrainCrystalPlasticity.h.

Referenced by postSolveStatevar(), and updateGss().

◆ _accslip_tmp_old

Real FiniteStrainCrystalPlasticity::_accslip_tmp_old

protected

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

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

◆ _crysrot

const MaterialProperty<RankTwoTensor>& FiniteStrainCrystalPlasticity::_crysrot

protected

Definition at line 322 of file FiniteStrainCrystalPlasticity.h.

Referenced by calc_schmid_tensor(), and postSolveQp().

◆ _deformation_gradient

const MaterialProperty<RankTwoTensor>& FiniteStrainCrystalPlasticity::_deformation_gradient

protected

Definition at line 316 of file FiniteStrainCrystalPlasticity.h.

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

◆ _deformation_gradient_old

const MaterialProperty<RankTwoTensor>& FiniteStrainCrystalPlasticity::_deformation_gradient_old

protected

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

Referenced by computeQpStress(), and preSolveQp().

◆ _dfgrd_tmp

RankTwoTensor FiniteStrainCrystalPlasticity::_dfgrd_tmp

protected

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

Referenced by computeQpStress(), and preSolveQp().

◆ _dgss_dsliprate

DenseMatrix<Real> FiniteStrainCrystalPlasticity::_dgss_dsliprate

protected

Definition at line 347 of file FiniteStrainCrystalPlasticity.h.

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

◆ _dslipdtau

DenseVector<Real> FiniteStrainCrystalPlasticity::_dslipdtau

protected

Definition at line 337 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 52 of file ComputeStressBase.h.

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

◆ _elasticity_tensor

const MaterialProperty<RankFourTensor>& FiniteStrainCrystalPlasticity::_elasticity_tensor

protected

Elasticity tensor material property.

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

◆ _err_tol

bool FiniteStrainCrystalPlasticity::_err_tol

protected

Definition at line 351 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 55 of file ComputeStressBase.h.

Referenced by ComputeStressBase::computeQpProperties().

◆ _fe

RankTwoTensor FiniteStrainCrystalPlasticity::_fe

protected

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

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

◆ _first_substep

bool FiniteStrainCrystalPlasticity::_first_substep

protected

Definition at line 358 of file FiniteStrainCrystalPlasticity.h.

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

◆ _flowprops

std::vector<Real> FiniteStrainCrystalPlasticity::_flowprops

protected

Definition at line 232 of file FiniteStrainCrystalPlasticity.h.

Referenced by getFlowRateParams().

◆ _fp

MaterialProperty<RankTwoTensor>& FiniteStrainCrystalPlasticity::_fp

protected

Definition at line 304 of file FiniteStrainCrystalPlasticity.h.

Referenced by initQpStatefulProperties(), and postSolveStress().

◆ _fp_inv

RankTwoTensor FiniteStrainCrystalPlasticity::_fp_inv

protected

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

Referenced by preSolveStress().

◆ _fp_old_inv

RankTwoTensor FiniteStrainCrystalPlasticity::_fp_old_inv

protected

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

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

◆ _gen_rndm_stress_flag

bool FiniteStrainCrystalPlasticity::_gen_rndm_stress_flag

protected

Definition at line 275 of file FiniteStrainCrystalPlasticity.h.

Referenced by postSolveQp().

◆ _gprops

std::vector<Real> FiniteStrainCrystalPlasticity::_gprops

protected

Definition at line 230 of file FiniteStrainCrystalPlasticity.h.

Referenced by getInitSlipSysRes().

◆ _gss

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

protected

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

Referenced by preSolveStatevar(), and updateGss().

◆ _gss_tmp

std::vector<Real> FiniteStrainCrystalPlasticity::_gss_tmp

protected

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

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

◆ _gtol

Real FiniteStrainCrystalPlasticity::_gtol

protected

Internal variable update equation tolerance.

Definition at line 254 of file FiniteStrainCrystalPlasticity.h.

Referenced by solveStatevar().

◆ _h0

Real FiniteStrainCrystalPlasticity::_h0

protected

Definition at line 330 of file FiniteStrainCrystalPlasticity.h.

Referenced by getHardnessParams(), and updateGss().

◆ _hprops

std::vector<Real> FiniteStrainCrystalPlasticity::_hprops

protected

Definition at line 231 of file FiniteStrainCrystalPlasticity.h.

Referenced by getHardnessParams(), and updateGss().

◆ _initial_stress_fcn

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

protectedinherited

initial stress components

Definition at line 58 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 278 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 270 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 61 of file ComputeStressBase.h.

◆ _lag_e

MaterialProperty<RankTwoTensor>& FiniteStrainCrystalPlasticity::_lag_e

protected

Definition at line 308 of file FiniteStrainCrystalPlasticity.h.

Referenced by initQpStatefulProperties(), and postSolveQp().

◆ _lag_e_old

const MaterialProperty<RankTwoTensor>& FiniteStrainCrystalPlasticity::_lag_e_old

protected

Definition at line 309 of file FiniteStrainCrystalPlasticity.h.

◆ _last_step_iter

bool FiniteStrainCrystalPlasticity::_last_step_iter

protected

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

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

◆ _lsrch_method

MooseEnum FiniteStrainCrystalPlasticity::_lsrch_method

protected

Definition at line 302 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 296 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 287 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 259 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 261 of file FiniteStrainCrystalPlasticity.h.

Referenced by solveStatevar().

◆ _mechanical_strain

const MaterialProperty<RankTwoTensor>& ComputeStressBase::_mechanical_strain

protectedinherited

Mechanical strain material property.

Definition at line 48 of file ComputeStressBase.h.

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

◆ _min_lsrch_step

Real FiniteStrainCrystalPlasticity::_min_lsrch_step

protected

Minimum line search step size.

Definition at line 293 of file FiniteStrainCrystalPlasticity.h.

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

◆ _mo

DenseVector<Real> FiniteStrainCrystalPlasticity::_mo

protected

Definition at line 324 of file FiniteStrainCrystalPlasticity.h.

Referenced by calc_schmid_tensor(), and getSlipSystems().

◆ _no

DenseVector<Real> FiniteStrainCrystalPlasticity::_no

protected

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

Referenced by FiniteStrainCrystalPlasticity(), and getSlipSystems().

◆ _pk2

MaterialProperty<RankTwoTensor>& FiniteStrainCrystalPlasticity::_pk2

protected

Definition at line 306 of file FiniteStrainCrystalPlasticity.h.

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

◆ _pk2_old

const MaterialProperty<RankTwoTensor>& FiniteStrainCrystalPlasticity::_pk2_old

protected

Definition at line 307 of file FiniteStrainCrystalPlasticity.h.

Referenced by preSolveStress().

◆ _pk2_tmp

RankTwoTensor FiniteStrainCrystalPlasticity::_pk2_tmp

protected

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

Referenced by postSolveStress(), and preSolveStress().

◆ _r

Real FiniteStrainCrystalPlasticity::_r

protected

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

Referenced by FiniteStrainCrystalPlasticity(), and getSlipSystems().

◆ _rndm_scale_var

Real FiniteStrainCrystalPlasticity::_rndm_scale_var

protected

Scaling value.

Definition at line 281 of file FiniteStrainCrystalPlasticity.h.

Referenced by postSolveQp().

◆ _rndm_seed

unsigned int FiniteStrainCrystalPlasticity::_rndm_seed

protected

Seed value.

Definition at line 284 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity().

◆ _rtol

Real FiniteStrainCrystalPlasticity::_rtol

protected

Stress residual equation relative tolerance.

Definition at line 250 of file FiniteStrainCrystalPlasticity.h.

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

◆ _s0

std::vector<RankTwoTensor> FiniteStrainCrystalPlasticity::_s0

protected

Definition at line 338 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 337 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 256 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 235 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 244 of file FiniteStrainCrystalPlasticity.h.

Referenced by initSlipSysProps(), and readFileFlowRateParams().