This class solves the viscoplastic flow rate equations in the total form Involves 4 different types of user objects that calculates: Internal variable rates - functions of internal variables and flow rates Internal variables - functions of internal variables Strengths - functions of internal variables Flow rates - functions of strengths and PK2 stress Flow directions - functions of strengths and PK2 stress The associated derivatives from user objects are assembled and the system is solved using NR. More...

#include <FiniteStrainHyperElasticViscoPlastic.h>

Inheritance diagram for FiniteStrainHyperElasticViscoPlastic:

[legend]

Public Member Functions
	FiniteStrainHyperElasticViscoPlastic (const InputParameters &parameters)

Protected Member Functions
virtual void	initUOVariables ()
	This function initializes the properties, stateful properties and user objects The properties and stateful properties associated with user objects are only initialized here The properties have the same name as the user object name. More...

void	initNumUserObjects (const std::vector< UserObjectName > &, unsigned int &)
	This function calculates the number of each user object type. More...

template<typename T >
void	initProp (const std::vector< UserObjectName > &, unsigned int, std::vector< MaterialProperty< T > *> &)
	This function initializes properties for each user object. More...

template<typename T >
void	initPropOld (const std::vector< UserObjectName > &, unsigned int, std::vector< const MaterialProperty< T > *> &)
	This function initializes old for stateful properties associated with user object Only user objects that update internal variables have an associated old property. More...

template<typename T >
void	initUserObjects (const std::vector< UserObjectName > &, unsigned int, std::vector< const T *> &)
	This function initializes user objects. More...

virtual void	initJacobianVariables ()
	This function initialize variables required for Jacobian calculation. More...

virtual void	initQpStatefulProperties ()
	Initializes state. More...

virtual void	computeQpStress ()
	This function computes the Cauchy stress. More...

virtual void	computeQpJacobian ()
	This function computes the Jacobian. More...

virtual void	saveOldState ()
	This function saves the old stateful properties that is modified during sub stepping. More...

virtual void	preSolveQp ()
	Sets state for solve. More...

virtual bool	solveQp ()
	Solve state. More...

virtual void	postSolveQp ()
	Update state for output (Outside substepping) More...

virtual void	recoverOldState ()
	This function restores the the old stateful properties after a successful solve. More...

virtual void	preSolveFlowrate ()
	Sets state for solve (Inside substepping) More...

virtual bool	solveFlowrate ()
	Solve for flow rate and state. More...

virtual void	postSolveFlowrate ()
	Update state for output (Inside substepping) More...

virtual bool	computeFlowRateFunction ()
	Calls user objects to compute flow rates. More...

virtual bool	computeFlowDirection ()
	Calls user objects to compute flow directions. More...

virtual void	computeElasticRightCauchyGreenTensor ()
	Computes elastic Right Cauchy Green Tensor. More...

virtual void	computePK2StressAndDerivative ()
	Computes PK2 stress and derivative w.r.t elastic Right Cauchy Green Tensor. More...

virtual void	computeElasticStrain ()
	Computes elastic Lagrangian strain. More...

virtual void	computeDeeDce ()
	Computes derivative of elastic strain w.r.t elastic Right Cauchy Green Tensor. More...

virtual bool	computeFlowRateResidual ()
	Computes flow rate residual vector. More...

virtual void	computeFlowRateJacobian ()
	Computes flow rate Jacobian matrix. More...

virtual void	computeElasticPlasticDeformGrad ()
	Computes elastic and plastic deformation gradients. More...

virtual Real	computeNorm (const std::vector< Real > &)
	Computes norm of residual vector. More...

virtual void	updateFlowRate ()
	Update flow rate. More...

virtual void	computeDpk2Dfpinv ()
	Computes derivative of PK2 stress wrt inverse of plastic deformation gradient. More...

virtual bool	computeIntVarRates ()
	This function call user objects to calculate rate of internal variables. More...

virtual bool	computeIntVar ()
	This function call user objects to integrate internal variables. More...

virtual bool	computeStrength ()
	This function call user objects to compute strength. More...

virtual void	computeIntVarRateDerivatives ()
	This function call user objects to compute dintvar_rate/dintvar and dintvarrate/dflowrate. More...

virtual void	computeIntVarDerivatives ()
	This function call user objects to compute dintvar/dintvar_rate and dintvar/dflowrate. More...

void	computeStrengthDerivatives ()
	This function call user objects to compute dstrength/dintvar. More...

virtual void	computeQpProperties () override

Protected Attributes
Real	_resid_abs_tol
	Absolute tolerance for residual convergence check. More...

Real	_resid_rel_tol
	Relative tolerance for residual convergence check. More...

unsigned int	_maxiters
	Maximum number of iterations. More...

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

std::vector< UserObjectName >	_flow_rate_uo_names
	Names of flow rate user objects. More...

std::vector< UserObjectName >	_strength_uo_names
	Names of strength user objects. More...

std::vector< UserObjectName >	_int_var_uo_names
	Names of internal variable user objects. More...

std::vector< UserObjectName >	_int_var_rate_uo_names
	Names of internal variable rate user objects. More...

unsigned int	_num_flow_rate_uos
	Number of flow rate user objects. More...

unsigned int	_num_strength_uos
	Number of strength user objects. More...

unsigned int	_num_int_var_uos
	Number of internal variable user objects. More...

unsigned int	_num_int_var_rate_uos
	Number of internal variable rate user objects. More...

std::vector< const HEVPFlowRateUOBase * >	_flow_rate_uo
	Flow rate user objects. More...

std::vector< const HEVPStrengthUOBase * >	_strength_uo
	Strength user objects. More...

std::vector< const HEVPInternalVarUOBase * >	_int_var_uo
	Internal variable user objects. More...

std::vector< const HEVPInternalVarRateUOBase * >	_int_var_rate_uo
	Internal variable rate user objects. More...

std::string	_pk2_prop_name

MaterialProperty< RankTwoTensor > &	_pk2

MaterialProperty< RankTwoTensor > &	_fp

const MaterialProperty< RankTwoTensor > &	_fp_old

MaterialProperty< RankTwoTensor > &	_ce

const MaterialProperty< RankTwoTensor > &	_deformation_gradient

const MaterialProperty< RankTwoTensor > &	_deformation_gradient_old

const MaterialProperty< RankTwoTensor > &	_rotation_increment

std::vector< MaterialProperty< Real > * >	_flow_rate_prop

std::vector< MaterialProperty< Real > * >	_strength_prop

std::vector< MaterialProperty< Real > * >	_int_var_stateful_prop

std::vector< const MaterialProperty< Real > * >	_int_var_stateful_prop_old

std::vector< MaterialProperty< Real > * >	_int_var_rate_prop

std::vector< Real >	_int_var_old

RankTwoTensor	_dfgrd_tmp

RankTwoTensor	_fp_tmp_inv

RankTwoTensor	_fp_tmp_old_inv

RankTwoTensor	_fe

RankTwoTensor	_ee

RankTwoTensor	_pk2_fet

RankTwoTensor	_fe_pk2

RankFourTensor	_dpk2_dce

RankFourTensor	_dpk2_dfe

RankFourTensor	_dfe_dfpinv

RankFourTensor	_dpk2_dfpinv

RankFourTensor	_dee_dce

RankFourTensor	_dce_dfe

RankFourTensor	_dfe_df

RankFourTensor	_tan_mod

RankFourTensor	_df_dstretch_inc

std::vector< RankTwoTensor >	_flow_dirn

std::vector< RankTwoTensor >	_dflowrate_dpk2

std::vector< RankTwoTensor >	_dpk2_dflowrate

std::vector< RankTwoTensor >	_dfpinv_dflowrate

DenseVector< Real >	_dflow_rate

DenseVector< Real >	_flow_rate

DenseVector< Real >	_resid

std::vector< DenseVector< Real > >	_dintvarrate_dflowrate
	Jacobian variables. More...

std::vector< DenseVector< Real > >	_dintvar_dflowrate_tmp

DenseMatrix< Real >	_dintvarrate_dintvar

DenseMatrix< Real >	_dintvar_dintvarrate

DenseMatrix< Real >	_dintvar_dintvar

DenseMatrix< Real >	_dintvar_dflowrate

DenseMatrix< Real >	_dstrength_dintvar

DenseMatrix< Real >	_dflowrate_dstrength

DenseVector< Real >	_dintvar_dintvar_x

DenseMatrix< Real >	_jac

Real	_dt_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< RankFourTensor > &	_elasticity_tensor

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

This class solves the viscoplastic flow rate equations in the total form Involves 4 different types of user objects that calculates: Internal variable rates - functions of internal variables and flow rates Internal variables - functions of internal variables Strengths - functions of internal variables Flow rates - functions of strengths and PK2 stress Flow directions - functions of strengths and PK2 stress The associated derivatives from user objects are assembled and the system is solved using NR.

Definition at line 34 of file FiniteStrainHyperElasticViscoPlastic.h.

Constructor & Destructor Documentation

◆ FiniteStrainHyperElasticViscoPlastic()

FiniteStrainHyperElasticViscoPlastic::FiniteStrainHyperElasticViscoPlastic ( const InputParameters & parameters )

Definition at line 44 of file FiniteStrainHyperElasticViscoPlastic.C.

   : ComputeStressBase(parameters),
     _resid_abs_tol(getParam<Real>("resid_abs_tol")),
     _resid_rel_tol(getParam<Real>("resid_rel_tol")),
     _maxiters(getParam<unsigned int>("maxiters")),
     _max_substep_iter(getParam<unsigned int>("max_substep_iteration")),
     _flow_rate_uo_names(isParamValid("flow_rate_user_objects")
                             ? getParam<std::vector<UserObjectName>>("flow_rate_user_objects")
                             : std::vector<UserObjectName>(0)),
     _strength_uo_names(isParamValid("strength_user_objects")
                            ? getParam<std::vector<UserObjectName>>("strength_user_objects")
                            : std::vector<UserObjectName>(0)),
     _int_var_uo_names(isParamValid("internal_var_user_objects")
                           ? getParam<std::vector<UserObjectName>>("internal_var_user_objects")
                           : std::vector<UserObjectName>(0)),
     _int_var_rate_uo_names(
         isParamValid("internal_var_rate_user_objects")
             ? getParam<std::vector<UserObjectName>>("internal_var_rate_user_objects")
             : std::vector<UserObjectName>(0)),
     _pk2_prop_name(_base_name + "pk2"),
     _pk2(declareProperty<RankTwoTensor>(_pk2_prop_name)),
     _fp(declareProperty<RankTwoTensor>(_base_name + "fp")),
     _fp_old(getMaterialPropertyOld<RankTwoTensor>(_base_name + "fp")),
     _ce(declareProperty<RankTwoTensor>(_base_name + "ce")),
     _deformation_gradient(getMaterialProperty<RankTwoTensor>(_base_name + "deformation_gradient")),
     _deformation_gradient_old(
         getMaterialPropertyOld<RankTwoTensor>(_base_name + "deformation_gradient")),
     _rotation_increment(getMaterialProperty<RankTwoTensor>(_base_name + "rotation_increment"))
 {
   initUOVariables();
 
   initJacobianVariables();
 
   _dflow_rate.resize(_num_flow_rate_uos);
   _flow_rate.resize(_num_flow_rate_uos);
   _resid.resize(_num_flow_rate_uos);
 
   _flow_dirn.resize(_num_flow_rate_uos);
 }

Member Function Documentation

◆ computeDeeDce()

void FiniteStrainHyperElasticViscoPlastic::computeDeeDce ( )

protectedvirtual

Computes derivative of elastic strain w.r.t elastic Right Cauchy Green Tensor.

Definition at line 490 of file FiniteStrainHyperElasticViscoPlastic.C.

Referenced by initJacobianVariables().

 {
   _dee_dce.zero();
 
   for (unsigned int i = 0; i < LIBMESH_DIM; ++i)
     for (unsigned int j = 0; j < LIBMESH_DIM; ++j)
       _dee_dce(i, j, i, j) = 0.5;
 }

◆ computeDpk2Dfpinv()

void FiniteStrainHyperElasticViscoPlastic::computeDpk2Dfpinv ( )

protectedvirtual

Computes derivative of PK2 stress wrt inverse of plastic deformation gradient.

Definition at line 521 of file FiniteStrainHyperElasticViscoPlastic.C.

Referenced by computeFlowRateJacobian().

 {
   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)
         _dfe_dfpinv(i, j, k, j) = _dfgrd_tmp(i, k);
 
   _dpk2_dfpinv = _dpk2_dce * _dce_dfe * _dfe_dfpinv;
 }

◆ computeElasticPlasticDeformGrad()

void FiniteStrainHyperElasticViscoPlastic::computeElasticPlasticDeformGrad ( )

protectedvirtual

Computes elastic and plastic deformation gradients.

Definition at line 506 of file FiniteStrainHyperElasticViscoPlastic.C.

Referenced by solveFlowrate().

 {
   RankTwoTensor iden;
   iden.addIa(1.0);
 
   RankTwoTensor val;
   for (unsigned int i = 0; i < _num_flow_rate_uos; ++i)
     val += _flow_rate(i) * _flow_dirn[i] * _dt_substep;
 
   _fp_tmp_inv = _fp_tmp_old_inv * (iden - val);
   _fp_tmp_inv = std::pow(_fp_tmp_inv.det(), -1.0 / 3.0) * _fp_tmp_inv;
   _fe = _dfgrd_tmp * _fp_tmp_inv;
 }

◆ computeElasticRightCauchyGreenTensor()

void FiniteStrainHyperElasticViscoPlastic::computeElasticRightCauchyGreenTensor ( )

protectedvirtual

Computes elastic Right Cauchy Green Tensor.

Definition at line 500 of file FiniteStrainHyperElasticViscoPlastic.C.

Referenced by computeFlowRateResidual().

 {
   _ce[_qp] = _fe.transpose() * _fe;
 }

◆ computeElasticStrain()

void FiniteStrainHyperElasticViscoPlastic::computeElasticStrain ( )

protectedvirtual

Computes elastic Lagrangian strain.

Definition at line 482 of file FiniteStrainHyperElasticViscoPlastic.C.

Referenced by HyperElasticPhaseFieldIsoDamage::computePK2StressAndDerivative(), and computePK2StressAndDerivative().

 {
   RankTwoTensor iden;
   iden.addIa(1.0);
   _ee = 0.5 * (_ce[_qp] - iden);
 }

◆ computeFlowDirection()

bool FiniteStrainHyperElasticViscoPlastic::computeFlowDirection ( )

protectedvirtual

Calls user objects to compute flow directions.

Definition at line 439 of file FiniteStrainHyperElasticViscoPlastic.C.

Referenced by computeFlowRateResidual().

 {
   for (unsigned i = 0; i < _num_flow_rate_uos; ++i)
   {
     if (!_flow_rate_uo[i]->computeDirection(_qp, _flow_dirn[i]))
       return false;
   }
   return true;
 }

◆ computeFlowRateFunction()

bool FiniteStrainHyperElasticViscoPlastic::computeFlowRateFunction ( )

protectedvirtual

Calls user objects to compute flow rates.

Definition at line 450 of file FiniteStrainHyperElasticViscoPlastic.C.

Referenced by computeFlowRateResidual().

 {
   Real val = 0;
   for (unsigned i = 0; i < _num_flow_rate_uos; ++i)
   {
     if (_flow_rate_uo[i]->computeValue(_qp, val))
       _resid(i) = -val;
     else
       return false;
   }
   return true;
 }

◆ computeFlowRateJacobian()

void FiniteStrainHyperElasticViscoPlastic::computeFlowRateJacobian ( )

protectedvirtual

Computes flow rate Jacobian matrix.

Definition at line 401 of file FiniteStrainHyperElasticViscoPlastic.C.

Referenced by solveFlowrate().

 {
   computeIntVarRateDerivatives();
   computeIntVarDerivatives();
   computeStrengthDerivatives();
 
   for (unsigned int i = 0; i < _num_flow_rate_uos; ++i)
     for (unsigned int j = 0; j < _num_strength_uos; ++j)
       _flow_rate_uo[i]->computeDerivative(_qp, _strength_uo_names[j], _dflowrate_dstrength(i, j));
 
   for (unsigned int i = 0; i < _num_flow_rate_uos; ++i)
     _flow_rate_uo[i]->computeTensorDerivative(_qp, _pk2_prop_name, _dflowrate_dpk2[i]);
 
   computeDpk2Dfpinv();
 
   for (unsigned int i = 0; i < _num_flow_rate_uos; ++i)
   {
     _dfpinv_dflowrate[i] = -_fp_tmp_old_inv * _flow_dirn[i] * _dt_substep;
     _dpk2_dflowrate[i] = _dpk2_dfpinv * _dfpinv_dflowrate[i];
   }
 
   DenseMatrix<Real> dflowrate_dflowrate;
   dflowrate_dflowrate = _dflowrate_dstrength;
   dflowrate_dflowrate.right_multiply(_dstrength_dintvar);
   dflowrate_dflowrate.right_multiply(_dintvar_dflowrate);
 
   _jac.zero();
   for (unsigned int i = 0; i < _num_flow_rate_uos; ++i)
     for (unsigned int j = 0; j < _num_flow_rate_uos; ++j)
     {
       if (i == j)
         _jac(i, j) = 1;
       _jac(i, j) -= dflowrate_dflowrate(i, j);
       _jac(i, j) -= _dflowrate_dpk2[i].doubleContraction(_dpk2_dflowrate[j]);
     }
 }

◆ computeFlowRateResidual()

bool FiniteStrainHyperElasticViscoPlastic::computeFlowRateResidual ( )

protectedvirtual

Computes flow rate residual vector.

Definition at line 375 of file FiniteStrainHyperElasticViscoPlastic.C.

Referenced by solveFlowrate().

 {
   if (!computeIntVarRates())
     return false;
 
   if (!computeIntVar())
     return false;
 
   if (!computeStrength())
     return false;
 
   computeElasticRightCauchyGreenTensor();
   computePK2StressAndDerivative();
 
   if (!computeFlowRateFunction())
     return false;
 
   if (!computeFlowDirection())
     return false;
 
   _resid += _flow_rate;
 
   return true;
 }

◆ computeIntVar()

bool FiniteStrainHyperElasticViscoPlastic::computeIntVar ( )

protectedvirtual

This function call user objects to integrate internal variables.

Definition at line 597 of file FiniteStrainHyperElasticViscoPlastic.C.

Referenced by computeFlowRateResidual().

 {
   Real val = 0;
   for (unsigned int i = 0; i < _num_int_var_uos; ++i)
   {
     if (_int_var_uo[i]->computeValue(_qp, _dt_substep, val))
       (*_int_var_stateful_prop[i])[_qp] = val;
     else
       return false;
   }
   return true;
 }

◆ computeIntVarDerivatives()

void FiniteStrainHyperElasticViscoPlastic::computeIntVarDerivatives ( )

protectedvirtual

This function call user objects to compute dintvar/dintvar_rate and dintvar/dflowrate.

Definition at line 638 of file FiniteStrainHyperElasticViscoPlastic.C.

Referenced by computeFlowRateJacobian().

 {
   Real val = 0;
 
   for (unsigned int i = 0; i < _num_int_var_uos; ++i)
     for (unsigned int j = 0; j < _num_int_var_rate_uos; ++j)
     {
       _int_var_uo[i]->computeDerivative(_qp, _dt_substep, _int_var_rate_uo_names[j], val);
       _dintvar_dintvarrate(i, j) = val;
     }
 
   _dintvar_dintvar.zero();
 
   for (unsigned int i = 0; i < _num_int_var_uos; ++i)
     for (unsigned int j = 0; j < _num_int_var_uos; ++j)
     {
       if (i == j)
         _dintvar_dintvar(i, j) = 1;
       for (unsigned int k = 0; k < _num_int_var_rate_uos; ++k)
         _dintvar_dintvar(i, j) -= _dintvar_dintvarrate(i, k) * _dintvarrate_dintvar(k, j);
     }
 
   for (unsigned int i = 0; i < _num_flow_rate_uos; ++i)
     _dintvar_dintvarrate.vector_mult(_dintvar_dflowrate_tmp[i], _dintvarrate_dflowrate[i]);
 
   for (unsigned int i = 0; i < _num_flow_rate_uos; ++i)
   {
     _dintvar_dintvar_x.zero();
     _dintvar_dintvar.lu_solve(_dintvar_dflowrate_tmp[i], _dintvar_dintvar_x);
     for (unsigned int j = 0; j < _num_int_var_uos; ++j)
       _dintvar_dflowrate(j, i) = _dintvar_dintvar_x(j);
   }
 }

◆ computeIntVarRateDerivatives()

void FiniteStrainHyperElasticViscoPlastic::computeIntVarRateDerivatives ( )

protectedvirtual

This function call user objects to compute dintvar_rate/dintvar and dintvarrate/dflowrate.

Definition at line 625 of file FiniteStrainHyperElasticViscoPlastic.C.

Referenced by computeFlowRateJacobian().

 {
   Real val = 0;
 
   for (unsigned int i = 0; i < _num_int_var_rate_uos; ++i)
     for (unsigned int j = 0; j < _num_flow_rate_uos; ++j)
     {
       _int_var_rate_uo[i]->computeDerivative(_qp, _flow_rate_uo_names[j], val);
       _dintvarrate_dflowrate[j](i) = val;
     }
 }

◆ computeIntVarRates()

bool FiniteStrainHyperElasticViscoPlastic::computeIntVarRates ( )

protectedvirtual

This function call user objects to calculate rate of internal variables.

Definition at line 583 of file FiniteStrainHyperElasticViscoPlastic.C.

Referenced by computeFlowRateResidual().

 {
   Real val = 0;
   for (unsigned int i = 0; i < _num_int_var_rate_uos; ++i)
   {
     if (_int_var_rate_uo[i]->computeValue(_qp, val))
       (*_int_var_rate_prop[i])[_qp] = val;
     else
       return false;
   }
   return true;
 }

◆ computeNorm()

Real FiniteStrainHyperElasticViscoPlastic::computeNorm ( const std::vector< Real > & var )

protectedvirtual

Computes norm of residual vector.

Definition at line 532 of file FiniteStrainHyperElasticViscoPlastic.C.

Referenced by solveFlowrate().

 {
   Real val = 0.0;
   for (unsigned int i = 0; i < var.size(); ++i)
     val += Utility::pow<2>(var[i]);
   return std::sqrt(val);
 }

◆ computePK2StressAndDerivative()

void FiniteStrainHyperElasticViscoPlastic::computePK2StressAndDerivative ( )

protectedvirtual

Computes PK2 stress and derivative w.r.t elastic Right Cauchy Green Tensor.

Reimplemented in HyperElasticPhaseFieldIsoDamage.

Definition at line 464 of file FiniteStrainHyperElasticViscoPlastic.C.

Referenced by computeFlowRateResidual().

 {
   computeElasticStrain();
   _pk2[_qp] = _elasticity_tensor[_qp] * _ee;
 
   _dce_dfe.zero();
   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)
       {
         _dce_dfe(i, j, k, i) = _dce_dfe(i, j, k, i) + _fe(k, j);
         _dce_dfe(i, j, k, j) = _dce_dfe(i, j, k, j) + _fe(k, i);
       }
 
   _dpk2_dfe = _dpk2_dce * _dce_dfe;
 }

◆ computeQpJacobian()

void FiniteStrainHyperElasticViscoPlastic::computeQpJacobian ( )

protectedvirtual

This function computes the Jacobian.

Reimplemented in HyperElasticPhaseFieldIsoDamage.

Definition at line 551 of file FiniteStrainHyperElasticViscoPlastic.C.

Referenced by HyperElasticPhaseFieldIsoDamage::computeQpJacobian(), and postSolveQp().

 {
   _tan_mod = _fe.mixedProductIkJl(_fe) * _dpk2_dfe;
   _pk2_fet = _pk2[_qp] * _fe.transpose();
   _fe_pk2 = _fe * _pk2[_qp];
 
   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) += _pk2_fet(l, j);
         _tan_mod(i, j, j, l) += _fe_pk2(i, l);
       }
 
   _tan_mod /= _fe.det();
 
   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)
         _dfe_df(i, j, i, l) = _fp_tmp_inv(l, j);
 
   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)
         for (unsigned int l = 0; l < LIBMESH_DIM; ++l)
           _df_dstretch_inc(i, j, k, l) =
               _rotation_increment[_qp](i, k) * _deformation_gradient_old[_qp](l, j);
 
   _Jacobian_mult[_qp] = _tan_mod * _dfe_df * _df_dstretch_inc;
 }

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

protectedvirtual

This function computes the Cauchy stress.

Implements ComputeStressBase.

Definition at line 210 of file FiniteStrainHyperElasticViscoPlastic.C.

 {
   bool converge;
   RankTwoTensor delta_dfgrd = _deformation_gradient[_qp] - _deformation_gradient_old[_qp];
   unsigned int num_substep = 1;
   unsigned int substep_iter = 1;
 
   saveOldState();
 
   do
   {
     preSolveQp();
 
     converge = true;
     _dt_substep = _dt / num_substep;
 
     for (unsigned int istep = 0; istep < num_substep; ++istep)
     {
       _dfgrd_tmp = (istep + 1.0) * delta_dfgrd / num_substep + _deformation_gradient_old[_qp];
       if (!solveQp())
       {
         converge = false;
         substep_iter++;
         num_substep *= 2;
         break;
       }
     }
 
     if (substep_iter > _max_substep_iter)
       mooseError("Constitutive failure with substepping at quadrature point ",
                  _q_point[_qp](0),
                  " ",
                  _q_point[_qp](1),
                  " ",
                  _q_point[_qp](2));
   } while (!converge);
 
   postSolveQp();
 }

◆ computeStrength()

bool FiniteStrainHyperElasticViscoPlastic::computeStrength ( )

protectedvirtual

This function call user objects to compute strength.

Definition at line 611 of file FiniteStrainHyperElasticViscoPlastic.C.

Referenced by computeFlowRateResidual().

 {
   Real val = 0;
   for (unsigned int i = 0; i < _num_strength_uos; ++i)
   {
     if (_strength_uo[i]->computeValue(_qp, val))
       (*_strength_prop[i])[_qp] = val;
     else
       return false;
   }
   return true;
 }

◆ computeStrengthDerivatives()

void FiniteStrainHyperElasticViscoPlastic::computeStrengthDerivatives ( )

protected

This function call user objects to compute dstrength/dintvar.

Definition at line 673 of file FiniteStrainHyperElasticViscoPlastic.C.

Referenced by computeFlowRateJacobian().

 {
   Real val = 0;
 
   for (unsigned int i = 0; i < _num_strength_uos; ++i)
     for (unsigned int j = 0; j < _num_int_var_uos; ++j)
     {
       _strength_uo[i]->computeDerivative(_qp, _int_var_uo_names[j], val);
       _dstrength_dintvar(i, j) = val;
     }
 }

◆ initJacobianVariables()

void FiniteStrainHyperElasticViscoPlastic::initJacobianVariables ( )

protectedvirtual

This function initialize variables required for Jacobian calculation.

Definition at line 154 of file FiniteStrainHyperElasticViscoPlastic.C.

Referenced by FiniteStrainHyperElasticViscoPlastic().

 {
   _dintvarrate_dflowrate.resize(_num_flow_rate_uos);
   for (unsigned int i = 0; i < _num_flow_rate_uos; ++i)
     _dintvarrate_dflowrate[i].resize(_num_int_var_rate_uos);
 
   _dintvar_dflowrate_tmp.resize(_num_flow_rate_uos);
   for (unsigned int i = 0; i < _num_flow_rate_uos; ++i)
     _dintvar_dflowrate_tmp[i].resize(_num_int_var_uos);
 
   _dintvarrate_dintvar.resize(_num_int_var_rate_uos, _num_int_var_uos);
   _dintvar_dintvarrate.resize(_num_int_var_uos, _num_int_var_rate_uos);
   _dintvar_dflowrate.resize(_num_int_var_uos, _num_flow_rate_uos);
   _dintvar_dintvar.resize(_num_int_var_uos, _num_int_var_uos);
   _dstrength_dintvar.resize(_num_strength_uos, _num_int_var_uos);
   _dflowrate_dstrength.resize(_num_flow_rate_uos, _num_strength_uos);
   _dintvar_dintvar_x.resize(_num_int_var_uos);
 
   _dpk2_dflowrate.resize(_num_flow_rate_uos);
   _dflowrate_dpk2.resize(_num_flow_rate_uos);
   _dfpinv_dflowrate.resize(_num_flow_rate_uos);
 
   _jac.resize(_num_flow_rate_uos, _num_flow_rate_uos);
 
   computeDeeDce();
 }

◆ initNumUserObjects()

void FiniteStrainHyperElasticViscoPlastic::initNumUserObjects	(	const std::vector< UserObjectName > &	uo_names,
		unsigned int &	uo_num
	)

protected

This function calculates the number of each user object type.

Definition at line 109 of file FiniteStrainHyperElasticViscoPlastic.C.

Referenced by initUOVariables().

 {
   uo_num = uo_names.size();
 }

◆ initProp()

template<typename T >

void FiniteStrainHyperElasticViscoPlastic::initProp	(	const std::vector< UserObjectName > &	uo_names,
		unsigned int	uo_num,
		std::vector< MaterialProperty< T > *> &	uo_prop
	)

protected

This function initializes properties for each user object.

Definition at line 117 of file FiniteStrainHyperElasticViscoPlastic.C.

Referenced by initUOVariables().

 {
   uo_prop.resize(uo_num);
   for (unsigned int i = 0; i < uo_num; ++i)
     uo_prop[i] = &declareProperty<T>(uo_names[i]);
 }

◆ initPropOld()

template<typename T >

void FiniteStrainHyperElasticViscoPlastic::initPropOld	(	const std::vector< UserObjectName > &	uo_names,
		unsigned int	uo_num,
		std::vector< const MaterialProperty< T > *> &	uo_prop_old
	)

protected

This function initializes old for stateful properties associated with user object Only user objects that update internal variables have an associated old property.

Definition at line 128 of file FiniteStrainHyperElasticViscoPlastic.C.

Referenced by initUOVariables().

 {
   uo_prop_old.resize(uo_num);
   for (unsigned int i = 0; i < uo_num; ++i)
     uo_prop_old[i] = &getMaterialPropertyOld<T>(uo_names[i]);
 }

◆ initQpStatefulProperties()

void FiniteStrainHyperElasticViscoPlastic::initQpStatefulProperties ( )

protectedvirtual

Initializes state.

Reimplemented from ComputeStressBase.

Definition at line 182 of file FiniteStrainHyperElasticViscoPlastic.C.

 {
   _stress[_qp].zero();
 
   _fp[_qp].zero();
   _fp[_qp].addIa(1.0);
   _ce[_qp].zero();
 
   _pk2[_qp].zero();
 
   for (unsigned int i = 0; i < _num_flow_rate_uos; ++i)
     (*_flow_rate_prop[i])[_qp] = 0.0;
 
   for (unsigned int i = 0; i < _num_strength_uos; ++i)
     (*_strength_prop[i])[_qp] = 0.0;
 
   for (unsigned int i = 0; i < _num_int_var_uos; ++i)
   {
     (*_int_var_stateful_prop[i])[_qp] = 0.0;
     // TODO: remove this nasty const_cast if you can figure out how
     const_cast<MaterialProperty<Real> &>(*_int_var_stateful_prop_old[i])[_qp] = 0.0;
   }
 
   for (unsigned int i = 0; i < _num_int_var_rate_uos; ++i)
     (*_int_var_rate_prop[i])[_qp] = 0.0;
 }

◆ initUOVariables()

void FiniteStrainHyperElasticViscoPlastic::initUOVariables ( )

protectedvirtual

This function initializes the properties, stateful properties and user objects The properties and stateful properties associated with user objects are only initialized here The properties have the same name as the user object name.

Definition at line 86 of file FiniteStrainHyperElasticViscoPlastic.C.

Referenced by FiniteStrainHyperElasticViscoPlastic().

 {
   initNumUserObjects(_flow_rate_uo_names, _num_flow_rate_uos);
   initNumUserObjects(_strength_uo_names, _num_strength_uos);
   initNumUserObjects(_int_var_uo_names, _num_int_var_uos);
   initNumUserObjects(_int_var_rate_uo_names, _num_int_var_rate_uos);
 
   initProp(_flow_rate_uo_names, _num_flow_rate_uos, _flow_rate_prop);
   initProp(_strength_uo_names, _num_strength_uos, _strength_prop);
   initProp(_int_var_uo_names, _num_int_var_uos, _int_var_stateful_prop);
   initProp(_int_var_rate_uo_names, _num_int_var_rate_uos, _int_var_rate_prop);
 
   initPropOld(_int_var_uo_names, _num_int_var_uos, _int_var_stateful_prop_old);
 
   initUserObjects(_flow_rate_uo_names, _num_flow_rate_uos, _flow_rate_uo);
   initUserObjects(_strength_uo_names, _num_strength_uos, _strength_uo);
   initUserObjects(_int_var_uo_names, _num_int_var_uos, _int_var_uo);
   initUserObjects(_int_var_rate_uo_names, _num_int_var_rate_uos, _int_var_rate_uo);
 
   _int_var_old.resize(_num_int_var_uos, 0.0);
 }

◆ initUserObjects()

template<typename T >

void FiniteStrainHyperElasticViscoPlastic::initUserObjects	(	const std::vector< UserObjectName > &	uo_names,
		unsigned int	uo_num,
		std::vector< const T *> &	uo
	)

protected

This function initializes user objects.

Definition at line 140 of file FiniteStrainHyperElasticViscoPlastic.C.

Referenced by initUOVariables().

 {
   uo.resize(uo_num);
 
   if (uo_num == 0)
     mooseError("Specify atleast one user object of type", typeid(T).name());
 
   for (unsigned int i = 0; i < uo_num; ++i)
     uo[i] = &getUserObjectByName<T>(uo_names[i]);
 }

◆ postSolveFlowrate()

void FiniteStrainHyperElasticViscoPlastic::postSolveFlowrate ( )

protectedvirtual

Update state for output (Inside substepping)

Definition at line 364 of file FiniteStrainHyperElasticViscoPlastic.C.

Referenced by solveQp().

 {
   _fp_tmp_old_inv = _fp_tmp_inv;
 
   // TODO: remove this nasty const_cast if you can figure out how
   for (unsigned int i = 0; i < _num_int_var_uos; ++i)
     const_cast<MaterialProperty<Real> &>(*_int_var_stateful_prop_old[i])[_qp] =
         (*_int_var_stateful_prop[i])[_qp];
 }

◆ postSolveQp()

void FiniteStrainHyperElasticViscoPlastic::postSolveQp ( )

protectedvirtual

Update state for output (Outside substepping)

Definition at line 282 of file FiniteStrainHyperElasticViscoPlastic.C.

Referenced by computeQpStress().

 {
   recoverOldState();
 
   _stress[_qp] = _fe * _pk2[_qp] * _fe.transpose() / _fe.det();
   _fp[_qp] = _fp_tmp_inv.inverse();
 
   computeQpJacobian();
 }

◆ preSolveFlowrate()

void FiniteStrainHyperElasticViscoPlastic::preSolveFlowrate ( )

protectedvirtual

Sets state for solve (Inside substepping)

Definition at line 301 of file FiniteStrainHyperElasticViscoPlastic.C.

Referenced by solveQp().

 {
   for (unsigned int i = 0; i < _num_flow_rate_uos; ++i)
   {
     _flow_rate(i) = 0.0;
     (*_flow_rate_prop[i])[_qp] = 0.0;
   }
 
   for (unsigned int i = 0; i < _num_int_var_uos; ++i)
     (*_int_var_stateful_prop[i])[_qp] = (*_int_var_stateful_prop_old[i])[_qp];
 
   _fp_tmp_inv = _fp_tmp_old_inv;
   _fe = _dfgrd_tmp * _fp_tmp_inv;
 }

◆ preSolveQp()

void FiniteStrainHyperElasticViscoPlastic::preSolveQp ( )

protectedvirtual

Sets state for solve.

Definition at line 258 of file FiniteStrainHyperElasticViscoPlastic.C.

Referenced by computeQpStress().

 {
   _fp_tmp_old_inv = _fp_old[_qp].inverse();
 
   // TODO: remove this nasty const_cast if you can figure out how
   for (unsigned int i = 0; i < _num_int_var_uos; ++i)
     (*_int_var_stateful_prop[i])[_qp] =
         const_cast<MaterialProperty<Real> &>(*_int_var_stateful_prop_old[i])[_qp] = _int_var_old[i];
 
   _dpk2_dce = _elasticity_tensor[_qp] * _dee_dce;
 }

◆ recoverOldState()

void FiniteStrainHyperElasticViscoPlastic::recoverOldState ( )

protectedvirtual

This function restores the the old stateful properties after a successful solve.

Definition at line 293 of file FiniteStrainHyperElasticViscoPlastic.C.

Referenced by postSolveQp().

 {
   // TODO: remove this nasty const_cast if you can figure out how
   for (unsigned int i = 0; i < _num_int_var_uos; ++i)
     const_cast<MaterialProperty<Real> &>(*_int_var_stateful_prop_old[i])[_qp] = _int_var_old[i];
 }

◆ saveOldState()

void FiniteStrainHyperElasticViscoPlastic::saveOldState ( )

protectedvirtual

This function saves the old stateful properties that is modified during sub stepping.

Definition at line 251 of file FiniteStrainHyperElasticViscoPlastic.C.

Referenced by computeQpStress().

 {
   for (unsigned int i = 0; i < _num_int_var_uos; ++i)
     _int_var_old[i] = (*_int_var_stateful_prop_old[i])[_qp];
 }

◆ solveFlowrate()

bool FiniteStrainHyperElasticViscoPlastic::solveFlowrate ( )

protectedvirtual

Solve for flow rate and state.

Definition at line 317 of file FiniteStrainHyperElasticViscoPlastic.C.

Referenced by solveQp().

 {
   Real resid0, rnorm;
   unsigned int iter = 0;
 
 #ifdef DEBUG
   std::vector<Real> rnormst(_maxiters + 1), flowratest(_maxiters + 1);
 #endif
 
   if (!computeFlowRateResidual())
     return false;
 
   rnorm = computeNorm(_resid.get_values());
   resid0 = rnorm;
 
 #ifdef DEBUG
   rnormst[iter] = rnorm;
   flowratest[iter] = computeNorm(_flow_rate.get_values());
 #endif
 
   while (rnorm > _resid_abs_tol && rnorm > _resid_rel_tol * resid0 && iter < _maxiters)
   {
     computeFlowRateJacobian();
 
     updateFlowRate();
 
     computeElasticPlasticDeformGrad();
 
     if (!computeFlowRateResidual())
       return false;
 
     rnorm = computeNorm(_resid.get_values());
     iter++;
 
 #ifdef DEBUG
     rnormst[iter] = rnorm;
     flowratest[iter] = computeNorm(_flow_rate.get_values());
 #endif
   }
 
   if (iter == _maxiters && rnorm > _resid_abs_tol && rnorm > _resid_rel_tol * resid0)
     return false;
 
   return true;
 }

◆ solveQp()

bool FiniteStrainHyperElasticViscoPlastic::solveQp ( )

protectedvirtual

Solve state.

Definition at line 271 of file FiniteStrainHyperElasticViscoPlastic.C.

Referenced by computeQpStress().

 {
   preSolveFlowrate();
   if (!solveFlowrate())
     return false;
   postSolveFlowrate();
 
   return true;
 }

◆ updateFlowRate()

void FiniteStrainHyperElasticViscoPlastic::updateFlowRate ( )

protectedvirtual

Update flow rate.

Definition at line 541 of file FiniteStrainHyperElasticViscoPlastic.C.

Referenced by solveFlowrate().

 {
   _jac.lu_solve(_resid, _dflow_rate);
   _flow_rate -= _dflow_rate;
 
   for (unsigned int i = 0; i < _num_flow_rate_uos; ++i)
     (*_flow_rate_prop[i])[_qp] = _flow_rate(i);
 }

Member Data Documentation

◆ _base_name

const std::string ComputeStressBase::_base_name

protectedinherited

Definition at line 37 of file ComputeStressBase.h.

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

◆ _ce

MaterialProperty<RankTwoTensor>& FiniteStrainHyperElasticViscoPlastic::_ce

protected

Definition at line 198 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeElasticRightCauchyGreenTensor(), computeElasticStrain(), and initQpStatefulProperties().

◆ _dce_dfe

RankFourTensor FiniteStrainHyperElasticViscoPlastic::_dce_dfe

protected

Definition at line 216 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeDpk2Dfpinv(), HyperElasticPhaseFieldIsoDamage::computePK2StressAndDerivative(), computePK2StressAndDerivative(), and HyperElasticPhaseFieldIsoDamage::computeQpJacobian().

◆ _dee_dce

RankFourTensor FiniteStrainHyperElasticViscoPlastic::_dee_dce

protected

Definition at line 216 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeDeeDce(), HyperElasticPhaseFieldIsoDamage::computePK2StressAndDerivative(), HyperElasticPhaseFieldIsoDamage::computeQpJacobian(), and preSolveQp().

◆ _deformation_gradient

const MaterialProperty<RankTwoTensor>& FiniteStrainHyperElasticViscoPlastic::_deformation_gradient

protected

Definition at line 200 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeQpStress().

◆ _deformation_gradient_old

const MaterialProperty<RankTwoTensor>& FiniteStrainHyperElasticViscoPlastic::_deformation_gradient_old

protected

Definition at line 201 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeQpJacobian(), and computeQpStress().

◆ _df_dstretch_inc

RankFourTensor FiniteStrainHyperElasticViscoPlastic::_df_dstretch_inc

protected

Definition at line 217 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by HyperElasticPhaseFieldIsoDamage::computeQpJacobian(), and computeQpJacobian().

◆ _dfe_df

RankFourTensor FiniteStrainHyperElasticViscoPlastic::_dfe_df

protected

Definition at line 216 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by HyperElasticPhaseFieldIsoDamage::computeQpJacobian(), and computeQpJacobian().

◆ _dfe_dfpinv

RankFourTensor FiniteStrainHyperElasticViscoPlastic::_dfe_dfpinv

protected

Definition at line 215 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeDpk2Dfpinv().

◆ _dfgrd_tmp

RankTwoTensor FiniteStrainHyperElasticViscoPlastic::_dfgrd_tmp

protected

Definition at line 211 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeDpk2Dfpinv(), computeElasticPlasticDeformGrad(), computeQpStress(), and preSolveFlowrate().

◆ _dflow_rate

DenseVector<Real> FiniteStrainHyperElasticViscoPlastic::_dflow_rate

protected

Definition at line 224 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by FiniteStrainHyperElasticViscoPlastic(), and updateFlowRate().

◆ _dflowrate_dpk2

std::vector<RankTwoTensor> FiniteStrainHyperElasticViscoPlastic::_dflowrate_dpk2

protected

Definition at line 220 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeFlowRateJacobian(), and initJacobianVariables().

◆ _dflowrate_dstrength

DenseMatrix<Real> FiniteStrainHyperElasticViscoPlastic::_dflowrate_dstrength

protected

Definition at line 237 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeFlowRateJacobian(), and initJacobianVariables().

◆ _dfpinv_dflowrate

std::vector<RankTwoTensor> FiniteStrainHyperElasticViscoPlastic::_dfpinv_dflowrate

protected

Definition at line 222 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeFlowRateJacobian(), and initJacobianVariables().

◆ _dintvar_dflowrate

DenseMatrix<Real> FiniteStrainHyperElasticViscoPlastic::_dintvar_dflowrate

protected

Definition at line 235 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeFlowRateJacobian(), computeIntVarDerivatives(), and initJacobianVariables().

◆ _dintvar_dflowrate_tmp

std::vector<DenseVector<Real> > FiniteStrainHyperElasticViscoPlastic::_dintvar_dflowrate_tmp

protected

Definition at line 230 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeIntVarDerivatives(), and initJacobianVariables().

◆ _dintvar_dintvar

DenseMatrix<Real> FiniteStrainHyperElasticViscoPlastic::_dintvar_dintvar

protected

Definition at line 234 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeIntVarDerivatives(), and initJacobianVariables().

◆ _dintvar_dintvar_x

DenseVector<Real> FiniteStrainHyperElasticViscoPlastic::_dintvar_dintvar_x

protected

Definition at line 238 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeIntVarDerivatives(), and initJacobianVariables().

◆ _dintvar_dintvarrate

DenseMatrix<Real> FiniteStrainHyperElasticViscoPlastic::_dintvar_dintvarrate

protected

Definition at line 233 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeIntVarDerivatives(), and initJacobianVariables().

◆ _dintvarrate_dflowrate

std::vector<DenseVector<Real> > FiniteStrainHyperElasticViscoPlastic::_dintvarrate_dflowrate

protected

Jacobian variables.

Definition at line 229 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeIntVarDerivatives(), computeIntVarRateDerivatives(), and initJacobianVariables().

◆ _dintvarrate_dintvar

DenseMatrix<Real> FiniteStrainHyperElasticViscoPlastic::_dintvarrate_dintvar

protected

Definition at line 232 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeIntVarDerivatives(), and initJacobianVariables().

◆ _dpk2_dce

RankFourTensor FiniteStrainHyperElasticViscoPlastic::_dpk2_dce

protected

Definition at line 215 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeDpk2Dfpinv(), HyperElasticPhaseFieldIsoDamage::computePK2StressAndDerivative(), computePK2StressAndDerivative(), and preSolveQp().

◆ _dpk2_dfe

RankFourTensor FiniteStrainHyperElasticViscoPlastic::_dpk2_dfe

protected

Definition at line 215 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by HyperElasticPhaseFieldIsoDamage::computePK2StressAndDerivative(), computePK2StressAndDerivative(), and computeQpJacobian().

◆ _dpk2_dflowrate

std::vector<RankTwoTensor> FiniteStrainHyperElasticViscoPlastic::_dpk2_dflowrate

protected

Definition at line 221 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeFlowRateJacobian(), and initJacobianVariables().

◆ _dpk2_dfpinv

RankFourTensor FiniteStrainHyperElasticViscoPlastic::_dpk2_dfpinv

protected

Definition at line 215 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeDpk2Dfpinv(), and computeFlowRateJacobian().

◆ _dstrength_dintvar

DenseMatrix<Real> FiniteStrainHyperElasticViscoPlastic::_dstrength_dintvar

protected

Definition at line 236 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeFlowRateJacobian(), computeStrengthDerivatives(), and initJacobianVariables().

◆ _dt_substep

Real FiniteStrainHyperElasticViscoPlastic::_dt_substep

protected

Definition at line 241 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeElasticPlasticDeformGrad(), computeFlowRateJacobian(), computeIntVar(), computeIntVarDerivatives(), and computeQpStress().

◆ _ee

RankTwoTensor FiniteStrainHyperElasticViscoPlastic::_ee

protected

Definition at line 213 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by HyperElasticPhaseFieldIsoDamage::computeDamageStress(), computeElasticStrain(), HyperElasticPhaseFieldIsoDamage::computeNumStiffness(), and computePK2StressAndDerivative().

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

protectedinherited

Definition at line 44 of file ComputeStressBase.h.

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

◆ _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 FiniteStrainHyperElasticViscoPlastic::_fe

protected

Definition at line 213 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeElasticPlasticDeformGrad(), computeElasticRightCauchyGreenTensor(), HyperElasticPhaseFieldIsoDamage::computePK2StressAndDerivative(), computePK2StressAndDerivative(), HyperElasticPhaseFieldIsoDamage::computeQpJacobian(), computeQpJacobian(), postSolveQp(), and preSolveFlowrate().

◆ _fe_pk2

RankTwoTensor FiniteStrainHyperElasticViscoPlastic::_fe_pk2

protected

Definition at line 214 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeQpJacobian().

◆ _flow_dirn

std::vector<RankTwoTensor> FiniteStrainHyperElasticViscoPlastic::_flow_dirn

protected

Definition at line 219 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeElasticPlasticDeformGrad(), computeFlowDirection(), computeFlowRateJacobian(), and FiniteStrainHyperElasticViscoPlastic().

◆ _flow_rate

DenseVector<Real> FiniteStrainHyperElasticViscoPlastic::_flow_rate

protected

Definition at line 225 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeElasticPlasticDeformGrad(), computeFlowRateResidual(), FiniteStrainHyperElasticViscoPlastic(), preSolveFlowrate(), solveFlowrate(), and updateFlowRate().

◆ _flow_rate_prop

std::vector<MaterialProperty<Real> *> FiniteStrainHyperElasticViscoPlastic::_flow_rate_prop

protected

Definition at line 204 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by initQpStatefulProperties(), initUOVariables(), preSolveFlowrate(), and updateFlowRate().

◆ _flow_rate_uo

std::vector<const HEVPFlowRateUOBase *> FiniteStrainHyperElasticViscoPlastic::_flow_rate_uo

protected

Flow rate user objects.

Definition at line 186 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeFlowDirection(), computeFlowRateFunction(), computeFlowRateJacobian(), and initUOVariables().

◆ _flow_rate_uo_names

std::vector<UserObjectName> FiniteStrainHyperElasticViscoPlastic::_flow_rate_uo_names

protected

Names of flow rate user objects.

Definition at line 168 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeIntVarRateDerivatives(), and initUOVariables().

◆ _fp

MaterialProperty<RankTwoTensor>& FiniteStrainHyperElasticViscoPlastic::_fp

protected

Definition at line 196 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by initQpStatefulProperties(), and postSolveQp().

◆ _fp_old

const MaterialProperty<RankTwoTensor>& FiniteStrainHyperElasticViscoPlastic::_fp_old

protected

Definition at line 197 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by preSolveQp().

◆ _fp_tmp_inv

RankTwoTensor FiniteStrainHyperElasticViscoPlastic::_fp_tmp_inv

protected

Definition at line 212 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeElasticPlasticDeformGrad(), computeQpJacobian(), postSolveFlowrate(), postSolveQp(), and preSolveFlowrate().

◆ _fp_tmp_old_inv

RankTwoTensor FiniteStrainHyperElasticViscoPlastic::_fp_tmp_old_inv

protected

Definition at line 212 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeElasticPlasticDeformGrad(), computeFlowRateJacobian(), postSolveFlowrate(), preSolveFlowrate(), and preSolveQp().

◆ _initial_stress_fcn

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

protectedinherited

initial stress components

Definition at line 50 of file ComputeStressBase.h.

◆ _int_var_old

std::vector<Real> FiniteStrainHyperElasticViscoPlastic::_int_var_old

protected

Definition at line 209 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by initUOVariables(), preSolveQp(), recoverOldState(), and saveOldState().

◆ _int_var_rate_prop

std::vector<MaterialProperty<Real> *> FiniteStrainHyperElasticViscoPlastic::_int_var_rate_prop

protected

Definition at line 208 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeIntVarRates(), initQpStatefulProperties(), and initUOVariables().

◆ _int_var_rate_uo

std::vector<const HEVPInternalVarRateUOBase *> FiniteStrainHyperElasticViscoPlastic::_int_var_rate_uo

protected

Internal variable rate user objects.

Definition at line 192 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeIntVarRateDerivatives(), computeIntVarRates(), and initUOVariables().

◆ _int_var_rate_uo_names

std::vector<UserObjectName> FiniteStrainHyperElasticViscoPlastic::_int_var_rate_uo_names

protected

Names of internal variable rate user objects.

Definition at line 174 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeIntVarDerivatives(), and initUOVariables().

◆ _int_var_stateful_prop

std::vector<MaterialProperty<Real> *> FiniteStrainHyperElasticViscoPlastic::_int_var_stateful_prop

protected

Definition at line 206 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeIntVar(), initQpStatefulProperties(), initUOVariables(), postSolveFlowrate(), preSolveFlowrate(), and preSolveQp().

◆ _int_var_stateful_prop_old

std::vector<const MaterialProperty<Real> *> FiniteStrainHyperElasticViscoPlastic::_int_var_stateful_prop_old

protected

Definition at line 207 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by initQpStatefulProperties(), initUOVariables(), postSolveFlowrate(), preSolveFlowrate(), preSolveQp(), recoverOldState(), and saveOldState().

◆ _int_var_uo

std::vector<const HEVPInternalVarUOBase *> FiniteStrainHyperElasticViscoPlastic::_int_var_uo

protected

Internal variable user objects.

Definition at line 190 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeIntVar(), computeIntVarDerivatives(), and initUOVariables().

◆ _int_var_uo_names

std::vector<UserObjectName> FiniteStrainHyperElasticViscoPlastic::_int_var_uo_names

protected

Names of internal variable user objects.

Definition at line 172 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeStrengthDerivatives(), and initUOVariables().

◆ _jac

DenseMatrix<Real> FiniteStrainHyperElasticViscoPlastic::_jac

protected

Definition at line 239 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeFlowRateJacobian(), initJacobianVariables(), and updateFlowRate().

◆ _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.

◆ _max_substep_iter

unsigned int FiniteStrainHyperElasticViscoPlastic::_max_substep_iter

protected

Maximum number of substep iterations.

Definition at line 165 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeQpStress().

◆ _maxiters

unsigned int FiniteStrainHyperElasticViscoPlastic::_maxiters

protected

Maximum number of iterations.

Definition at line 163 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by solveFlowrate().

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

◆ _num_flow_rate_uos

unsigned int FiniteStrainHyperElasticViscoPlastic::_num_flow_rate_uos

protected

Number of flow rate user objects.

Definition at line 177 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeElasticPlasticDeformGrad(), computeFlowDirection(), computeFlowRateFunction(), computeFlowRateJacobian(), computeIntVarDerivatives(), computeIntVarRateDerivatives(), FiniteStrainHyperElasticViscoPlastic(), initJacobianVariables(), initQpStatefulProperties(), initUOVariables(), preSolveFlowrate(), and updateFlowRate().

◆ _num_int_var_rate_uos

unsigned int FiniteStrainHyperElasticViscoPlastic::_num_int_var_rate_uos

protected

Number of internal variable rate user objects.

Definition at line 183 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeIntVarDerivatives(), computeIntVarRateDerivatives(), computeIntVarRates(), initJacobianVariables(), initQpStatefulProperties(), and initUOVariables().

◆ _num_int_var_uos

unsigned int FiniteStrainHyperElasticViscoPlastic::_num_int_var_uos

protected

Number of internal variable user objects.

Definition at line 181 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeIntVar(), computeIntVarDerivatives(), computeStrengthDerivatives(), initJacobianVariables(), initQpStatefulProperties(), initUOVariables(), postSolveFlowrate(), preSolveFlowrate(), preSolveQp(), recoverOldState(), and saveOldState().

◆ _num_strength_uos

unsigned int FiniteStrainHyperElasticViscoPlastic::_num_strength_uos

protected

Number of strength user objects.

Definition at line 179 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeFlowRateJacobian(), computeStrength(), computeStrengthDerivatives(), initJacobianVariables(), initQpStatefulProperties(), and initUOVariables().

◆ _pk2

MaterialProperty<RankTwoTensor>& FiniteStrainHyperElasticViscoPlastic::_pk2

protected

Definition at line 195 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by HyperElasticPhaseFieldIsoDamage::computeNumStiffness(), HyperElasticPhaseFieldIsoDamage::computePK2StressAndDerivative(), computePK2StressAndDerivative(), computeQpJacobian(), initQpStatefulProperties(), and postSolveQp().

◆ _pk2_fet

RankTwoTensor FiniteStrainHyperElasticViscoPlastic::_pk2_fet

protected

Definition at line 214 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeQpJacobian().

◆ _pk2_prop_name

std::string FiniteStrainHyperElasticViscoPlastic::_pk2_prop_name

protected

Definition at line 194 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeFlowRateJacobian().

◆ _resid

DenseVector<Real> FiniteStrainHyperElasticViscoPlastic::_resid

protected

Definition at line 226 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeFlowRateFunction(), computeFlowRateResidual(), FiniteStrainHyperElasticViscoPlastic(), solveFlowrate(), and updateFlowRate().

◆ _resid_abs_tol

Real FiniteStrainHyperElasticViscoPlastic::_resid_abs_tol

protected

Absolute tolerance for residual convergence check.

Definition at line 159 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by solveFlowrate().

◆ _resid_rel_tol

Real FiniteStrainHyperElasticViscoPlastic::_resid_rel_tol

protected

Relative tolerance for residual convergence check.

Definition at line 161 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by solveFlowrate().

◆ _rotation_increment

const MaterialProperty<RankTwoTensor>& FiniteStrainHyperElasticViscoPlastic::_rotation_increment

protected

Definition at line 202 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeQpJacobian().

◆ _strength_prop

std::vector<MaterialProperty<Real> *> FiniteStrainHyperElasticViscoPlastic::_strength_prop

protected

Definition at line 205 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeStrength(), initQpStatefulProperties(), and initUOVariables().

◆ _strength_uo

std::vector<const HEVPStrengthUOBase *> FiniteStrainHyperElasticViscoPlastic::_strength_uo

protected

Strength user objects.

Definition at line 188 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeStrength(), computeStrengthDerivatives(), and initUOVariables().

◆ _strength_uo_names

std::vector<UserObjectName> FiniteStrainHyperElasticViscoPlastic::_strength_uo_names

protected

Names of strength user objects.

Definition at line 170 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeFlowRateJacobian(), and initUOVariables().

◆ _stress

MaterialProperty<RankTwoTensor>& ComputeStressBase::_stress

protectedinherited

Definition at line 41 of file ComputeStressBase.h.

◆ _tan_mod

RankFourTensor FiniteStrainHyperElasticViscoPlastic::_tan_mod

protected

Definition at line 217 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeQpJacobian().

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

tensor_mechanics/include/materials/FiniteStrainHyperElasticViscoPlastic.h
tensor_mechanics/src/materials/FiniteStrainHyperElasticViscoPlastic.C

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ FiniteStrainHyperElasticViscoPlastic()

Member Function Documentation

◆ computeDeeDce()

◆ computeDpk2Dfpinv()

◆ computeElasticPlasticDeformGrad()

◆ computeElasticRightCauchyGreenTensor()

◆ computeElasticStrain()

◆ computeFlowDirection()

◆ computeFlowRateFunction()

◆ computeFlowRateJacobian()

◆ computeFlowRateResidual()

◆ computeIntVar()

◆ computeIntVarDerivatives()

◆ computeIntVarRateDerivatives()

◆ computeIntVarRates()

◆ computeNorm()

◆ computePK2StressAndDerivative()

◆ computeQpJacobian()

◆ computeQpProperties()

◆ computeQpStress()

◆ computeStrength()

◆ computeStrengthDerivatives()

◆ initJacobianVariables()

◆ initNumUserObjects()

◆ initProp()

◆ initPropOld()

◆ initQpStatefulProperties()

◆ initUOVariables()

◆ initUserObjects()

◆ postSolveFlowrate()

◆ postSolveQp()

◆ preSolveFlowrate()

◆ preSolveQp()

◆ recoverOldState()

◆ saveOldState()

◆ solveFlowrate()

◆ solveQp()

◆ updateFlowRate()

Member Data Documentation

◆ _base_name

◆ _ce

◆ _dce_dfe

◆ _dee_dce

◆ _deformation_gradient

◆ _deformation_gradient_old

◆ _df_dstretch_inc

◆ _dfe_df

◆ _dfe_dfpinv

◆ _dfgrd_tmp

◆ _dflow_rate

◆ _dflowrate_dpk2

◆ _dflowrate_dstrength

◆ _dfpinv_dflowrate

◆ _dintvar_dflowrate

◆ _dintvar_dflowrate_tmp

◆ _dintvar_dintvar

◆ _dintvar_dintvar_x

◆ _dintvar_dintvarrate

◆ _dintvarrate_dflowrate

◆ _dintvarrate_dintvar

◆ _dpk2_dce

◆ _dpk2_dfe

◆ _dpk2_dflowrate

◆ _dpk2_dfpinv

◆ _dstrength_dintvar

◆ _dt_substep

◆ _ee

◆ _elastic_strain

◆ _elasticity_tensor

◆ _elasticity_tensor_name

◆ _extra_stress

◆ _fe

◆ _fe_pk2

◆ _flow_dirn

◆ _flow_rate