FiniteStrainCPSlipRateRes Class Reference

#include <FiniteStrainCPSlipRateRes.h>

Inheritance diagram for FiniteStrainCPSlipRateRes:

Public Member Functions
	FiniteStrainCPSlipRateRes (const InputParameters &parameters)

Protected Member Functions
virtual void	solveStatevar ()
	This function solves internal variables. More...
virtual void	preSolveStress ()
	This function sets variable for internal variable solve. More...
virtual void	solveStress ()
	This function solves for stress, updates plastic deformation gradient. More...
virtual void	calcResidJacobSlipRate ()
	This function calculates residual and jacobian of slip rate. More...
virtual void	calcResidualSlipRate ()
	This function calculates residual of slip rate. More...
virtual void	calcJacobianSlipRate ()
	This function calculates jacobian of slip rate. More...
virtual void	getSlipIncrements ()
	This function updates the slip system resistances. More...
virtual void	calcDtauDsliprate ()
	This function calculates partial derivative of resolved shear stress with respect to split rate. More...
virtual void	calcDgssDsliprate ()
	This function calculates partial derivative of slip system resistances with respect to split rate. More...
void	calcUpdate ()
	This function calculates and updates the residual of slip rate. More...
virtual Real	calcResidNorm ()
	This function calculates the residual norm. More...
bool	lineSearchUpdateSlipRate (const Real, const DenseVector< Real > &)
	This function performs the line search update. More...
Real	calcResidDotProdUpdate (const DenseVector< Real > &)
	This function calculates the dot product of residual and update. More...
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	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	postSolveStatevar ()
	This function update internal variable after solve. 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
DenseVector< Real >	_resid
DenseVector< Real >	_slip_rate
DenseVector< Real >	_dsliprate_dgss
DenseMatrix< Real >	_jacob
DenseMatrix< Real >	_dsliprate_dsliprate
const unsigned int	_nss
	Number of slip system resistance. 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 MaterialProperty< RankFourTensor > &	_elasticity_tensor
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
const std::string	_elasticity_tensor_name
const MaterialProperty< RankTwoTensor > &	_mechanical_strain
MaterialProperty< RankTwoTensor > &	_stress
MaterialProperty< RankTwoTensor > &	_elastic_strain
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

Definition at line 20 of file FiniteStrainCPSlipRateRes.h.

Constructor & Destructor Documentation

FiniteStrainCPSlipRateRes()

FiniteStrainCPSlipRateRes::FiniteStrainCPSlipRateRes

(

const InputParameters &

parameters

)

Definition at line 23 of file FiniteStrainCPSlipRateRes.C.

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

Member Function Documentation

assignSlipSysRes()

void FiniteStrainCrystalPlasticity::assignSlipSysRes ( )

protectedvirtualinherited

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

Definition at line 233 of file FiniteStrainCrystalPlasticity.C.

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

protectedvirtualinherited

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

Definition at line 881 of file FiniteStrainCrystalPlasticity.C.

Referenced by FiniteStrainCrystalPlasticity::solveStress().

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

calc_schmid_tensor()

void FiniteStrainCrystalPlasticity::calc_schmid_tensor ( )

protectedinherited

This function calculate the Schmid tensor.

Definition at line 1052 of file FiniteStrainCrystalPlasticity.C.

Referenced by FiniteStrainCrystalPlasticity::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);
}

calcDgssDsliprate()

void FiniteStrainCPSlipRateRes::calcDgssDsliprate ( )

protectedvirtual

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

Definition at line 222 of file FiniteStrainCPSlipRateRes.C.

Referenced by calcJacobianSlipRate().

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

calcDtauDsliprate()

void FiniteStrainCPSlipRateRes::calcDtauDsliprate ( )

protectedvirtual

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

Definition at line 187 of file FiniteStrainCPSlipRateRes.C.

Referenced by calcJacobianSlipRate().

{
  RankFourTensor dfedfpinv, deedfe, dfpinvdpk2;
  std::vector<RankTwoTensor> dtaudpk2(_nss), dfpinvdsliprate(_nss);

  for (unsigned int i = 0; i < _nss; ++i)
  {
    dtaudpk2[i] = _s0[i];
    dfpinvdsliprate[i] = -_fp_old_inv * _s0[i] * _dt;
  }

  for (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;
      }

  RankFourTensor dpk2dfpinv;

  dpk2dfpinv = _elasticity_tensor[_qp] * deedfe * dfedfpinv;

  for (unsigned int i = 0; i < _nss; ++i)
    for (unsigned int j = 0; j < _nss; ++j)
      _dsliprate_dsliprate(i, j) =
          _dslipdtau(i) * dtaudpk2[i].doubleContraction(dpk2dfpinv * dfpinvdsliprate[j]);
}

calcJacobian()

void FiniteStrainCrystalPlasticity::calcJacobian ( RankFourTensor &  jac )

protectedvirtualinherited

This function calculate jacobian.

Definition at line 930 of file FiniteStrainCrystalPlasticity.C.

Referenced by FiniteStrainCrystalPlasticity::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);
}

calcJacobianSlipRate()

void FiniteStrainCPSlipRateRes::calcJacobianSlipRate ( )

protectedvirtual

This function calculates jacobian of slip rate.

Definition at line 170 of file FiniteStrainCPSlipRateRes.C.

Referenced by calcResidJacobSlipRate(), and solveStress().

{
  //_dsliprate_dsliprate not reinitialized to zero, hence order is important
  calcDtauDsliprate();
  calcDgssDsliprate();

  for (unsigned int i = 0; i < _nss; ++i)
    for (unsigned int j = 0; j < _nss; ++j)
    {
      _jacob(i, j) = 0.0;
      if (i == j)
        _jacob(i, j) += 1.0;
      _jacob(i, j) -= _dsliprate_dsliprate(i, j);
    }
}

calcResidDotProdUpdate()

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

protected

This function calculates the dot product of residual and update.

Definition at line 368 of file FiniteStrainCPSlipRateRes.C.

Referenced by lineSearchUpdateSlipRate().

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

calcResidJacobSlipRate()

void FiniteStrainCPSlipRateRes::calcResidJacobSlipRate ( )

protectedvirtual

This function calculates residual and jacobian of slip rate.

Definition at line 125 of file FiniteStrainCPSlipRateRes.C.

Referenced by solveStress().

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

calcResidNorm()

Real FiniteStrainCPSlipRateRes::calcResidNorm ( )

protectedvirtual

This function calculates the residual norm.

Definition at line 260 of file FiniteStrainCPSlipRateRes.C.

Referenced by lineSearchUpdateSlipRate(), and solveStress().

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

  return rnorm;
}

calcResidual()

void FiniteStrainCrystalPlasticity::calcResidual ( RankTwoTensor &  resid )

protectedvirtualinherited

This function calculate stress residual.

Definition at line 890 of file FiniteStrainCrystalPlasticity.C.

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

{
  RankTwoTensor iden, ce, ee, ce_pk2, eqv_slip_incr, pk2_new;

  iden.zero();
  iden.addIa(1.0);

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

calcResidualSlipRate()

void FiniteStrainCPSlipRateRes::calcResidualSlipRate ( )

protectedvirtual

This function calculates residual of slip rate.

Definition at line 134 of file FiniteStrainCPSlipRateRes.C.

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

{
  RankTwoTensor eqv_slip_incr, ce, ee;
  const RankTwoTensor iden(RankTwoTensor::initIdentity);

  _slip_incr = _slip_rate;
  _slip_incr *= _dt;

  for (unsigned int i = 0; i < _nss; ++i)
    eqv_slip_incr += _s0[i] * _slip_incr(i);

  eqv_slip_incr = iden - eqv_slip_incr;

  _fp_inv = _fp_old_inv * eqv_slip_incr;
  _fe = _dfgrd_tmp * _fp_inv;

  ce = _fe.transpose() * _fe;
  ee = ce - iden;
  ee *= 0.5;

  _pk2_tmp = _elasticity_tensor[_qp] * ee;

  for (unsigned int i = 0; i < _nss; ++i)
    _tau(i) = _pk2_tmp.doubleContraction(_s0[i]);

  update_slip_system_resistance();
  getSlipIncrements();

  if (_err_tol)
    return;

  for (unsigned int i = 0; i < _nss; ++i)
    _resid(i) = _slip_rate(i) - _slip_incr(i) / _dt;
}

calcTangentModuli()

RankFourTensor FiniteStrainCrystalPlasticity::calcTangentModuli ( )

protectedvirtualinherited

This function calculate the tangent moduli for preconditioner.

Default is the elastic stiffness matrix. Exact jacobian is currently implemented. tan_mod_type can be modified to exact in .i file to turn it on.

Definition at line 1035 of file FiniteStrainCrystalPlasticity.C.

Referenced by FiniteStrainCrystalPlasticity::postSolveQp().

{
  RankFourTensor tan_mod;

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

  return tan_mod;
}

calcUpdate()

void FiniteStrainCPSlipRateRes::calcUpdate ( )

protected

This function calculates and updates the residual of slip rate.

Definition at line 246 of file FiniteStrainCPSlipRateRes.C.

Referenced by solveStress().

{
  DenseMatrix<Real> A = _jacob;
  DenseVector<Real> r(_nss);
  DenseVector<Real> x(_nss);

  r = _resid;

  A.lu_solve(r, x);

  _resid = x;
}

computeQpElasticityTensor()

void FiniteStrainCrystalPlasticity::computeQpElasticityTensor ( )

protectedvirtualinherited

This function updates the elasticity tensor at a quadrature point.

Presently void.

Definition at line 1030 of file FiniteStrainCrystalPlasticity.C.

{
}

computeQpProperties()

void ComputeStressBase::computeQpProperties ( )

overrideprotectedvirtualinherited

Definition at line 51 of file ComputeStressBase.C.

{
  computeQpStress();

  // Add in extra stress
  _stress[_qp] += _extra_stress[_qp];
}

computeQpStress()

void FiniteStrainCrystalPlasticity::computeQpStress ( )

protectedvirtualinherited

This function updates the stress at a quadrature point.

Solves stress residual equation using NR.

Updates slip system resistances iteratively.

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

protectedvirtualinherited

This function calculate the elastic tangent moduli for preconditioner.

Definition at line 1123 of file FiniteStrainCrystalPlasticity.C.

Referenced by FiniteStrainCrystalPlasticity::calcTangentModuli().

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

elastoPlasticTangentModuli()

RankFourTensor FiniteStrainCrystalPlasticity::elastoPlasticTangentModuli ( )

protectedvirtualinherited

This function calculate the exact tangent moduli for preconditioner.

Definition at line 1084 of file FiniteStrainCrystalPlasticity.C.

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

protectedinherited

This function perform RU decomposition to obtain the rotation tensor.

Definition at line 988 of file FiniteStrainCrystalPlasticity.C.

Referenced by FiniteStrainCrystalPlasticity::postSolveQp().

{
  return getMatRot(a);
}

getFlowRateParams()

void FiniteStrainCrystalPlasticity::getFlowRateParams ( )

protectedvirtualinherited

This function assign flow rate parameters - see test.

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

Definition at line 343 of file FiniteStrainCrystalPlasticity.C.

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

protectedvirtualinherited

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

Definition at line 411 of file FiniteStrainCrystalPlasticity.C.

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

protectedvirtualinherited

This function assign slip system resistances - see test.

.i input file format start_slip_sys_num, end_slip_sys_num, value.

Definition at line 261 of file FiniteStrainCrystalPlasticity.C.

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

protectedinherited

This function perform RU decomposition to obtain the rotation tensor.

Definition at line 995 of file FiniteStrainCrystalPlasticity.C.

Referenced by FiniteStrainCrystalPlasticity::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 FiniteStrainCPSlipRateRes::getSlipIncrements ( )

protectedvirtual

This function updates the slip system resistances.

Reimplemented from FiniteStrainCrystalPlasticity.

Definition at line 230 of file FiniteStrainCPSlipRateRes.C.

Referenced by calcResidualSlipRate().

{
  FiniteStrainCrystalPlasticity::getSlipIncrements();

  if (_err_tol)
    return;

  _dslipdtau *= 1.0 / _dt;

  for (unsigned int i = 0; i < _nss; ++i)
    _dsliprate_dgss(i) = -_a0(i) / _xm(i) *
                         std::pow(std::abs(_tau(i) / _gss_tmp[i]), 1.0 / _xm(i) - 1.0) * _tau(i) /
                         std::pow(_gss_tmp[i], 2.0);
}

getSlipSystems()

void FiniteStrainCrystalPlasticity::getSlipSystems ( )

protectedvirtualinherited

This function reads slip system from file - see test.

Definition at line 425 of file FiniteStrainCrystalPlasticity.C.

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

protectedvirtualinherited

This function initializes additional parameters.

Definition at line 483 of file FiniteStrainCrystalPlasticity.C.

Referenced by FiniteStrainCrystalPlasticity::initQpStatefulProperties().

{
}

initQpStatefulProperties()

void FiniteStrainCrystalPlasticity::initQpStatefulProperties ( )

protectedvirtualinherited

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

Reimplemented from ComputeStressBase.

Definition at line 188 of file FiniteStrainCrystalPlasticity.C.

{
  _stress[_qp].zero();

  _fp[_qp].zero();
  _fp[_qp].addIa(1.0);

  _pk2[_qp].zero();
  _acc_slip[_qp] = 0.0;
  _lag_e[_qp].zero();

  _update_rot[_qp].zero();
  _update_rot[_qp].addIa(1.0);

  initSlipSysProps(); // Initializes slip system related properties
  initAdditionalProps();
}

initSlipSysProps()

void FiniteStrainCrystalPlasticity::initSlipSysProps ( )

protectedvirtualinherited

This function initializes slip system resistances.

Definition at line 207 of file FiniteStrainCrystalPlasticity.C.

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

protectedinherited

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

Definition at line 1212 of file FiniteStrainCrystalPlasticity.C.

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

protectedinherited

This function performs the line search update.

Definition at line 1129 of file FiniteStrainCrystalPlasticity.C.

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

lineSearchUpdateSlipRate()

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

protected

This function performs the line search update.

Definition at line 271 of file FiniteStrainCPSlipRateRes.C.

Referenced by solveStress().

{
  if (_lsrch_method == "CUT_HALF")
  {
    Real rnorm;
    Real step = 1.0;
    do
    {
      for (unsigned int i = 0; i < update.size(); ++i)
        _slip_rate(i) += step * update(i);

      step /= 2.0;

      for (unsigned int i = 0; i < update.size(); ++i)
        _slip_rate(i) -= step * update(i);

      calcResidualSlipRate();
      if (_err_tol)
        return false;
      rnorm = calcResidNorm();
    } while (rnorm > rnorm_prev && step > _min_lsrch_step);

    if (rnorm > rnorm_prev && step <= _min_lsrch_step)
      return false;

    return true;
  }
  else if (_lsrch_method == "BISECTION")
  {
    unsigned int count = 0;
    Real step_a = 0.0;
    Real step_b = 1.0;
    Real step = 1.0;
    Real s_m = 1000.0;
    Real rnorm = 1000.0;

    Real s_b = calcResidDotProdUpdate(update);
    Real rnorm1 = calcResidNorm();

    for (unsigned int i = 0; i < update.size(); ++i)
      _slip_rate(i) += update(i);

    calcResidualSlipRate();
    Real s_a = calcResidDotProdUpdate(update);
    Real rnorm0 = calcResidNorm();

    for (unsigned int i = 0; i < update.size(); ++i)
      _slip_rate(i) -= update(i);

    if ((rnorm1 / rnorm0) < _lsrch_tol || s_a * s_b > 0)
    {
      calcResidualSlipRate();
      return true;
    }

    while ((rnorm / rnorm0) > _lsrch_tol && count < _lsrch_max_iter)
    {

      for (unsigned int i = 0; i < update.size(); ++i)
        _slip_rate(i) += step * update(i);

      step = 0.5 * (step_a + step_b);

      for (unsigned int i = 0; i < update.size(); ++i)
        _slip_rate(i) -= step * update(i);

      calcResidualSlipRate();
      s_m = calcResidDotProdUpdate(update);
      rnorm = calcResidNorm();

      if (s_m * s_a < 0.0)
      {
        step_b = step;
        s_b = s_m;
      }
      if (s_m * s_b < 0.0)
      {
        step_a = step;
        s_a = s_m;
      }
      count++;
    }

    if ((rnorm / rnorm0) < _lsrch_tol && count < _lsrch_max_iter)
      return true;

    return false;
  }
  else
  {
    mooseError("Line search meothod is not provided.");
    return false;
  }
}

postSolveQp()

void FiniteStrainCrystalPlasticity::postSolveQp ( )

protectedvirtualinherited

This function update stress and internal variable after solve.

Definition at line 587 of file FiniteStrainCrystalPlasticity.C.

Referenced by FiniteStrainCrystalPlasticity::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;
    iden.addIa(1.0);

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

protectedvirtualinherited

This function update internal variable after solve.

Definition at line 683 of file FiniteStrainCrystalPlasticity.C.

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

protectedvirtualinherited

This function update stress and plastic deformation gradient after solve.

Definition at line 802 of file FiniteStrainCrystalPlasticity.C.

Referenced by solveStatevar(), and FiniteStrainCrystalPlasticity::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 ( )

protectedvirtualinherited

This function set variables for stress and internal variable solve.

Definition at line 559 of file FiniteStrainCrystalPlasticity.C.

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

protectedvirtualinherited

This function set variables for internal variable solve.

Definition at line 621 of file FiniteStrainCrystalPlasticity.C.

Referenced by FiniteStrainCrystalPlasticity::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 FiniteStrainCPSlipRateRes::preSolveStress ( )

protectedvirtual

This function sets variable for internal variable solve.

Reimplemented from FiniteStrainCrystalPlasticity.

Definition at line 44 of file FiniteStrainCPSlipRateRes.C.

Referenced by solveStatevar().

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

readFileFlowRateParams()

void FiniteStrainCrystalPlasticity::readFileFlowRateParams ( )

protectedvirtualinherited

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

Definition at line 314 of file FiniteStrainCrystalPlasticity.C.

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

protectedvirtualinherited

This function read hardness parameters from file.

Definition at line 405 of file FiniteStrainCrystalPlasticity.C.

Referenced by FiniteStrainCrystalPlasticity::initSlipSysProps().

{
}

readFileInitSlipSysRes()

void FiniteStrainCrystalPlasticity::readFileInitSlipSysRes ( )

protectedvirtualinherited

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

Definition at line 243 of file FiniteStrainCrystalPlasticity.C.

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

protectedvirtualinherited

This function solves stress and internal variables.

Definition at line 577 of file FiniteStrainCrystalPlasticity.C.

Referenced by FiniteStrainCrystalPlasticity::computeQpStress().

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

solveStatevar()

void FiniteStrainCPSlipRateRes::solveStatevar ( )

protectedvirtual

This function solves internal variables.

Reimplemented from FiniteStrainCrystalPlasticity.

Definition at line 34 of file FiniteStrainCPSlipRateRes.C.

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

solveStress()

void FiniteStrainCPSlipRateRes::solveStress ( )

protectedvirtual

This function solves for stress, updates plastic deformation gradient.

Reimplemented from FiniteStrainCrystalPlasticity.

Definition at line 51 of file FiniteStrainCPSlipRateRes.C.

Referenced by solveStatevar().

{
  Real rnorm, rnorm0, rnorm_prev;
  unsigned int iter = 0;

#ifdef DEBUG
  std::vector<Real> rnormst(_maxiter + 1), slipratest(_maxiter + 1); // Use for Debugging
#endif

  calcResidJacobSlipRate();
  if (_err_tol)
    return;
  rnorm = calcResidNorm();
  rnorm0 = rnorm;

#ifdef DEBUG
  rnormst[iter] = rnorm;
  Real slipratemax = 0.0;
  for (unsigned int i = 0; i < _nss; ++i)
    if (std::abs(_slip_rate(i)) > slipratemax)
      slipratemax = std::abs(_slip_rate(i));
  slipratest[iter] = slipratemax;
#endif

  while (rnorm > _rtol * rnorm0 && rnorm0 > _abs_tol && iter < _maxiter)
  {
    calcUpdate();

    DenseVector<Real> update = _resid;

    _slip_rate -= update;

    calcResidualSlipRate();
    if (_err_tol)
      return;
    rnorm_prev = rnorm;
    rnorm = calcResidNorm();

    if (_use_line_search && rnorm > rnorm_prev && !lineSearchUpdateSlipRate(rnorm_prev, update))
    {
#ifdef DEBUG
      mooseWarning("FiniteStrainCrystalPLasticity: Failed with line search");
#endif
      _err_tol = true;
      return;
    }

    calcJacobianSlipRate();

    if (_use_line_search)
      rnorm = calcResidNorm();
    iter++;

#ifdef DEBUG
    slipratemax = 0.0;
    for (unsigned int i = 0; i < _nss; ++i)
      if (std::abs(_slip_rate(i)) > slipratemax)
        slipratemax = std::abs(_slip_rate(i));
    rnormst[iter] = rnorm;
    slipratest[iter] = slipratemax;
#endif
  }

  if (iter == _maxiter)
  {
#ifdef DEBUG
    mooseWarning("FiniteStrainCPSlipRateRes: NR exceeds maximum iteration ", iter, " ", rnorm);
#endif
    _err_tol = true;
    return;
  }
}

update_slip_system_resistance()

void FiniteStrainCrystalPlasticity::update_slip_system_resistance ( )

protectedvirtualinherited

This function updates the slip system resistances.

Definition at line 826 of file FiniteStrainCrystalPlasticity.C.

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

{
  updateGss();
}

updateGss()

void FiniteStrainCrystalPlasticity::updateGss ( )

protectedvirtualinherited

This function updates the slip system resistances.

Old function to update slip system resistances.

Kept to avoid code break at computeQpstress

Definition at line 836 of file FiniteStrainCrystalPlasticity.C.

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

protectedinherited

Definition at line 323 of file FiniteStrainCrystalPlasticity.h.

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

_abs_tol

Real FiniteStrainCrystalPlasticity::_abs_tol

protectedinherited

Stress residual equation absolute tolerance.

Definition at line 251 of file FiniteStrainCrystalPlasticity.h.

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

_acc_slip

MaterialProperty<Real>& FiniteStrainCrystalPlasticity::_acc_slip

protectedinherited

Definition at line 311 of file FiniteStrainCrystalPlasticity.h.

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

_acc_slip_old

const MaterialProperty<Real>& FiniteStrainCrystalPlasticity::_acc_slip_old

protectedinherited

Definition at line 312 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::preSolveStatevar().

_accslip_tmp

Real FiniteStrainCrystalPlasticity::_accslip_tmp

protectedinherited

Definition at line 337 of file FiniteStrainCrystalPlasticity.h.

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

_accslip_tmp_old

Real FiniteStrainCrystalPlasticity::_accslip_tmp_old

protectedinherited

Definition at line 337 of file FiniteStrainCrystalPlasticity.h.

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

_base_name

const std::string ComputeStressBase::_base_name

protectedinherited

Definition at line 37 of file ComputeStressBase.h.

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

_crysrot

const MaterialProperty<RankTwoTensor>& FiniteStrainCrystalPlasticity::_crysrot

protectedinherited

Definition at line 318 of file FiniteStrainCrystalPlasticity.h.

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

_deformation_gradient

const MaterialProperty<RankTwoTensor>& FiniteStrainCrystalPlasticity::_deformation_gradient

protectedinherited

Definition at line 315 of file FiniteStrainCrystalPlasticity.h.

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

_deformation_gradient_old

const MaterialProperty<RankTwoTensor>& FiniteStrainCrystalPlasticity::_deformation_gradient_old

protectedinherited

Definition at line 316 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::computeQpStress().

_delta_dfgrd

RankTwoTensor FiniteStrainCrystalPlasticity::_delta_dfgrd

protectedinherited

Flag to check whether convergence is achieved.

Used for substepping; Uniformly divides the increment in deformation gradient

Definition at line 350 of file FiniteStrainCrystalPlasticity.h.

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

_dfgrd_scale_factor

Real FiniteStrainCrystalPlasticity::_dfgrd_scale_factor

protectedinherited

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

Definition at line 352 of file FiniteStrainCrystalPlasticity.h.

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

_dfgrd_tmp

RankTwoTensor FiniteStrainCrystalPlasticity::_dfgrd_tmp

protectedinherited

Definition at line 331 of file FiniteStrainCrystalPlasticity.h.

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

_dfgrd_tmp_old

RankTwoTensor FiniteStrainCrystalPlasticity::_dfgrd_tmp_old

protectedinherited

Definition at line 350 of file FiniteStrainCrystalPlasticity.h.

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

_dgss_dsliprate

DenseMatrix<Real> FiniteStrainCrystalPlasticity::_dgss_dsliprate

protectedinherited

Definition at line 343 of file FiniteStrainCrystalPlasticity.h.

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

_dslipdtau

DenseVector<Real> FiniteStrainCrystalPlasticity::_dslipdtau

protectedinherited

Definition at line 333 of file FiniteStrainCrystalPlasticity.h.

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

_dsliprate_dgss

DenseVector<Real> FiniteStrainCPSlipRateRes::_dsliprate_dgss

protected

Definition at line 95 of file FiniteStrainCPSlipRateRes.h.

Referenced by calcDgssDsliprate(), and getSlipIncrements().

_dsliprate_dsliprate

DenseMatrix<Real> FiniteStrainCPSlipRateRes::_dsliprate_dsliprate

protected

Definition at line 97 of file FiniteStrainCPSlipRateRes.h.

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

_elastic_strain

MaterialProperty<RankTwoTensor>& ComputeStressBase::_elastic_strain

protectedinherited

Definition at line 42 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

protectedinherited

Definition at line 317 of file FiniteStrainCrystalPlasticity.h.

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

_elasticity_tensor_name

const std::string ComputeStressBase::_elasticity_tensor_name

protectedinherited

Definition at line 38 of file ComputeStressBase.h.

Referenced by ComputeFiniteStrainElasticStress::initialSetup(), ComputeSmearedCrackingStress::initialSetup(), and ComputeMultipleInelasticStress::initialSetup().

_err_tol

bool FiniteStrainCrystalPlasticity::_err_tol

protectedinherited

Definition at line 347 of file FiniteStrainCrystalPlasticity.h.

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

_extra_stress

const MaterialProperty<RankTwoTensor>& ComputeStressBase::_extra_stress

protectedinherited

Extra stress tensor.

Definition at line 47 of file ComputeStressBase.h.

Referenced by ComputeStressBase::computeQpProperties().

_fe

RankTwoTensor FiniteStrainCrystalPlasticity::_fe

protectedinherited

Definition at line 332 of file FiniteStrainCrystalPlasticity.h.

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

_first_step_iter

bool FiniteStrainCrystalPlasticity::_first_step_iter

protectedinherited

Flags to reset variables and reinitialize variables.

Definition at line 354 of file FiniteStrainCrystalPlasticity.h.

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

_first_substep

bool FiniteStrainCrystalPlasticity::_first_substep

protectedinherited

Definition at line 354 of file FiniteStrainCrystalPlasticity.h.

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

_flowprops

std::vector<Real> FiniteStrainCrystalPlasticity::_flowprops

protectedinherited

Definition at line 231 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::getFlowRateParams().

_fp

MaterialProperty<RankTwoTensor>& FiniteStrainCrystalPlasticity::_fp

protectedinherited

Definition at line 303 of file FiniteStrainCrystalPlasticity.h.

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

_fp_inv

RankTwoTensor FiniteStrainCrystalPlasticity::_fp_inv

protectedinherited

Definition at line 332 of file FiniteStrainCrystalPlasticity.h.

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

_fp_old

const MaterialProperty<RankTwoTensor>& FiniteStrainCrystalPlasticity::_fp_old

protectedinherited

Definition at line 304 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::preSolveStress().

_fp_old_inv

RankTwoTensor FiniteStrainCrystalPlasticity::_fp_old_inv

protectedinherited

Definition at line 332 of file FiniteStrainCrystalPlasticity.h.

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

_fp_prev_inv

RankTwoTensor FiniteStrainCrystalPlasticity::_fp_prev_inv

protectedinherited

Definition at line 332 of file FiniteStrainCrystalPlasticity.h.

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

_gen_rndm_stress_flag

bool FiniteStrainCrystalPlasticity::_gen_rndm_stress_flag

protectedinherited

Definition at line 274 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::postSolveQp().

_gprops

std::vector<Real> FiniteStrainCrystalPlasticity::_gprops

protectedinherited

Definition at line 229 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::getInitSlipSysRes().

_gss

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

protectedinherited

Definition at line 309 of file FiniteStrainCrystalPlasticity.h.

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

_gss_old

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

protectedinherited

Definition at line 310 of file FiniteStrainCrystalPlasticity.h.

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

_gss_tmp

std::vector<Real> FiniteStrainCrystalPlasticity::_gss_tmp

protectedinherited

Definition at line 338 of file FiniteStrainCrystalPlasticity.h.

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

_gss_tmp_old

std::vector<Real> FiniteStrainCrystalPlasticity::_gss_tmp_old

protectedinherited

Definition at line 339 of file FiniteStrainCrystalPlasticity.h.

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

_gtol

Real FiniteStrainCrystalPlasticity::_gtol

protectedinherited

Internal variable update equation tolerance.

Definition at line 253 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::solveStatevar().

_h0

Real FiniteStrainCrystalPlasticity::_h0

protectedinherited

Definition at line 326 of file FiniteStrainCrystalPlasticity.h.

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

_hprops

std::vector<Real> FiniteStrainCrystalPlasticity::_hprops

protectedinherited

Definition at line 230 of file FiniteStrainCrystalPlasticity.h.

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

_initial_stress_fcn

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

protectedinherited

initial stress components

Definition at line 50 of file ComputeStressBase.h.

_input_rndm_scale_var

bool FiniteStrainCrystalPlasticity::_input_rndm_scale_var

protectedinherited

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

Definition at line 277 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::postSolveQp().

_intvar_read_type

MooseEnum FiniteStrainCrystalPlasticity::_intvar_read_type

protectedinherited

Read from options for initial values of internal variables.

Definition at line 269 of file FiniteStrainCrystalPlasticity.h.

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

_jacob

DenseMatrix<Real> FiniteStrainCPSlipRateRes::_jacob

protected

Definition at line 96 of file FiniteStrainCPSlipRateRes.h.

Referenced by calcJacobianSlipRate(), and calcUpdate().

_Jacobian_mult

MaterialProperty<RankFourTensor>& ComputeStressBase::_Jacobian_mult

protectedinherited

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

Definition at line 53 of file ComputeStressBase.h.

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

_lag_e

MaterialProperty<RankTwoTensor>& FiniteStrainCrystalPlasticity::_lag_e

protectedinherited

Definition at line 307 of file FiniteStrainCrystalPlasticity.h.

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

_lag_e_old

const MaterialProperty<RankTwoTensor>& FiniteStrainCrystalPlasticity::_lag_e_old

protectedinherited

Definition at line 308 of file FiniteStrainCrystalPlasticity.h.

_last_step_iter

bool FiniteStrainCrystalPlasticity::_last_step_iter

protectedinherited

Definition at line 354 of file FiniteStrainCrystalPlasticity.h.

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

_lsrch_max_iter

unsigned int FiniteStrainCrystalPlasticity::_lsrch_max_iter

protectedinherited

Line search bisection method maximum iteration number.

Definition at line 298 of file FiniteStrainCrystalPlasticity.h.

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

_lsrch_method

MooseEnum FiniteStrainCrystalPlasticity::_lsrch_method

protectedinherited

Definition at line 301 of file FiniteStrainCrystalPlasticity.h.

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

_lsrch_tol

Real FiniteStrainCrystalPlasticity::_lsrch_tol

protectedinherited

Line search bisection method tolerance.

Definition at line 295 of file FiniteStrainCrystalPlasticity.h.

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

_max_substep_iter

unsigned int FiniteStrainCrystalPlasticity::_max_substep_iter

protectedinherited

Maximum number of substep iterations.

Definition at line 286 of file FiniteStrainCrystalPlasticity.h.

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

_maxiter

unsigned int FiniteStrainCrystalPlasticity::_maxiter

protectedinherited

Maximum number of iterations for stress update.

Definition at line 258 of file FiniteStrainCrystalPlasticity.h.

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

_maxiterg

unsigned int FiniteStrainCrystalPlasticity::_maxiterg

protectedinherited

Maximum number of iterations for internal variable update.

Definition at line 260 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::solveStatevar().

_mechanical_strain

const MaterialProperty<RankTwoTensor>& ComputeStressBase::_mechanical_strain

protectedinherited

Definition at line 40 of file ComputeStressBase.h.

Referenced by ComputeLinearElasticPFFractureStress::computeQpStress(), ComputeCosseratLinearElasticStress::computeQpStress(), ComputeFiniteStrainElasticStress::computeQpStress(), ComputeIsotropicLinearElasticPFFractureStress::computeQpStress(), and FiniteStrainPlasticMaterial::computeQpStress().

_min_lsrch_step

Real FiniteStrainCrystalPlasticity::_min_lsrch_step

protectedinherited

Minimum line search step size.

Definition at line 292 of file FiniteStrainCrystalPlasticity.h.

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

_mo

DenseVector<Real> FiniteStrainCrystalPlasticity::_mo

protectedinherited

Definition at line 320 of file FiniteStrainCrystalPlasticity.h.

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

_no

DenseVector<Real> FiniteStrainCrystalPlasticity::_no

protectedinherited

Definition at line 321 of file FiniteStrainCrystalPlasticity.h.

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

_nss

const unsigned int FiniteStrainCrystalPlasticity::_nss

protectedinherited

Number of slip system resistance.

Definition at line 227 of file FiniteStrainCrystalPlasticity.h.

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

_num_slip_sys_flowrate_props

unsigned int FiniteStrainCrystalPlasticity::_num_slip_sys_flowrate_props

protectedinherited

Number of slip system flow rate parameters.

Definition at line 263 of file FiniteStrainCrystalPlasticity.h.

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

_num_slip_sys_props

unsigned int FiniteStrainCrystalPlasticity::_num_slip_sys_props

protectedinherited

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

Definition at line 272 of file FiniteStrainCrystalPlasticity.h.

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

_pk2

MaterialProperty<RankTwoTensor>& FiniteStrainCrystalPlasticity::_pk2

protectedinherited

Definition at line 305 of file FiniteStrainCrystalPlasticity.h.

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

_pk2_old

const MaterialProperty<RankTwoTensor>& FiniteStrainCrystalPlasticity::_pk2_old

protectedinherited

Definition at line 306 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::preSolveStress().

_pk2_tmp

RankTwoTensor FiniteStrainCrystalPlasticity::_pk2_tmp

protectedinherited

Definition at line 336 of file FiniteStrainCrystalPlasticity.h.

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

_pk2_tmp_old

RankTwoTensor FiniteStrainCrystalPlasticity::_pk2_tmp_old

protectedinherited

Definition at line 336 of file FiniteStrainCrystalPlasticity.h.

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

_r

Real FiniteStrainCrystalPlasticity::_r

protectedinherited

Definition at line 329 of file FiniteStrainCrystalPlasticity.h.

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

_read_from_slip_sys_file

bool FiniteStrainCrystalPlasticity::_read_from_slip_sys_file

protectedinherited

Definition at line 345 of file FiniteStrainCrystalPlasticity.h.

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

_resid

DenseVector<Real> FiniteStrainCPSlipRateRes::_resid

protected

Definition at line 93 of file FiniteStrainCPSlipRateRes.h.

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

_rndm_scale_var

Real FiniteStrainCrystalPlasticity::_rndm_scale_var

protectedinherited

Scaling value.

Definition at line 280 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::postSolveQp().

_rndm_seed

unsigned int FiniteStrainCrystalPlasticity::_rndm_seed

protectedinherited

Seed value.

Definition at line 283 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::FiniteStrainCrystalPlasticity().

_rtol

Real FiniteStrainCrystalPlasticity::_rtol

protectedinherited

Stress residual equation relative tolerance.

Definition at line 249 of file FiniteStrainCrystalPlasticity.h.

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

_s0

std::vector<RankTwoTensor> FiniteStrainCrystalPlasticity::_s0

protectedinherited

Definition at line 334 of file FiniteStrainCrystalPlasticity.h.

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

_slip_incr

DenseVector<Real> FiniteStrainCrystalPlasticity::_slip_incr

protectedinherited

Definition at line 333 of file FiniteStrainCrystalPlasticity.h.

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

_slip_incr_tol

Real FiniteStrainCrystalPlasticity::_slip_incr_tol

protectedinherited

Slip increment tolerance.

Definition at line 255 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::getSlipIncrements().

_slip_rate

DenseVector<Real> FiniteStrainCPSlipRateRes::_slip_rate

protected

Definition at line 94 of file FiniteStrainCPSlipRateRes.h.

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

_slip_sys_file_name

std::string FiniteStrainCrystalPlasticity::_slip_sys_file_name

protectedinherited

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

Definition at line 234 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::getSlipSystems().

_slip_sys_flow_prop_file_name

std::string FiniteStrainCrystalPlasticity::_slip_sys_flow_prop_file_name

protectedinherited

File should contain values of the flow rate equation parameters.

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

Definition at line 243 of file FiniteStrainCrystalPlasticity.h.

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

_slip_sys_hard_prop_file_name

std::string FiniteStrainCrystalPlasticity::_slip_sys_hard_prop_file_name

protectedinherited

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

Definition at line 246 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::initSlipSysProps().

_slip_sys_props

DenseVector<Real> FiniteStrainCrystalPlasticity::_slip_sys_props

protectedinherited

Definition at line 341 of file FiniteStrainCrystalPlasticity.h.

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

_slip_sys_res_prop_file_name

std::string FiniteStrainCrystalPlasticity::_slip_sys_res_prop_file_name

protectedinherited

File should contain initial values of the slip system resistances.

Definition at line 237 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::readFileInitSlipSysRes().

_stress

MaterialProperty<RankTwoTensor>& ComputeStressBase::_stress

protectedinherited

Definition at line 41 of file ComputeStressBase.h.

Referenced by ComputeMultipleInelasticCosseratStress::computeAdmissibleState(), ComputeMultipleInelasticStress::computeAdmissibleState(), ComputeStressBase::computeQpProperties(), ComputeStrainIncrementBasedStress::computeQpStress(), ComputeLinearElasticStress::computeQpStress(), ComputeLinearElasticPFFractureStress::computeQpStress(), ComputeCosseratLinearElasticStress::computeQpStress(), ComputeDamageStress::computeQpStress(), ComputeSmearedCrackingStress::computeQpStress(), ComputeFiniteStrainElasticStress::computeQpStress(), ComputeIsotropicLinearElasticPFFractureStress::computeQpStress(), FiniteStrainPlasticMaterial::computeQpStress(), ComputeMultiPlasticityStress::computeQpStress(), ComputeLinearViscoelasticStress::computeQpStress(), ComputeMultipleInelasticStress::computeQpStressIntermediateConfiguration(), ComputeMultipleInelasticStress::finiteStrainRotation(), ComputeStressBase::initQpStatefulProperties(), FiniteStrainCrystalPlasticity::initQpStatefulProperties(), FiniteStrainUObasedCP::initQpStatefulProperties(), FiniteStrainHyperElasticViscoPlastic::initQpStatefulProperties(), ComputeMultiPlasticityStress::postReturnMap(), FiniteStrainUObasedCP::postSolveQp(), FiniteStrainHyperElasticViscoPlastic::postSolveQp(), FiniteStrainCrystalPlasticity::postSolveQp(), ComputeSmearedCrackingStress::updateCrackingStateAndStress(), ComputeMultipleInelasticStress::updateQpState(), and ComputeMultipleInelasticStress::updateQpStateSingleModel().

_tan_mod_type

MooseEnum FiniteStrainCrystalPlasticity::_tan_mod_type

protectedinherited

Type of tangent moduli calculation.

Definition at line 266 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::calcTangentModuli().

_tau

DenseVector<Real> FiniteStrainCrystalPlasticity::_tau

protectedinherited

Definition at line 333 of file FiniteStrainCrystalPlasticity.h.

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

_tau_init

Real FiniteStrainCrystalPlasticity::_tau_init

protectedinherited

Definition at line 328 of file FiniteStrainCrystalPlasticity.h.

Referenced by FiniteStrainCrystalPlasticity::getHardnessParams().

_tau_sat

Real FiniteStrainCrystalPlasticity::_tau_sat

protectedinherited

Definition at line 327 of file FiniteStrainCrystalPlasticity.h.

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

_update_rot

MaterialProperty<RankTwoTensor>& FiniteStrainCrystalPlasticity::_update_rot

protectedinherited

Definition at line 313 of file FiniteStrainCrystalPlasticity.h.

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

_use_line_search

bool FiniteStrainCrystalPlasticity::_use_line_search

protectedinherited

Flag to activate line serach.

Definition at line 289 of file FiniteStrainCrystalPlasticity.h.

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

_xm

DenseVector<Real> FiniteStrainCrystalPlasticity::_xm

protectedinherited

Definition at line 324 of file FiniteStrainCrystalPlasticity.h.

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

---

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

• tensor_mechanics/include/materials/FiniteStrainCPSlipRateRes.h
• tensor_mechanics/src/materials/FiniteStrainCPSlipRateRes.C