https://mooseframework.inl.gov
Loading...
Searching...
No Matches
Public Types | Public Member Functions | Static Public Member Functions | Protected Member Functions | Protected Attributes | Private Member Functions | Private Attributes | List of all members
FiniteStrainHyperElasticViscoPlastic Class Reference

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 Types

typedef DerivativeMaterialPropertyNameInterface::SymbolName SymbolName
 

Public Member Functions

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

Static Public Member Functions

static InputParameters validParams ()
 

Protected Member Functions

virtual void 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.
 
void initNumUserObjects (const std::vector< UserObjectName > &, unsigned int &)
 This function calculates the number of each user object type.
 
template<typename T >
void initProp (const std::vector< UserObjectName > &, unsigned int, std::vector< MaterialProperty< T > * > &)
 This function initializes properties for each user object.
 
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.
 
template<typename T >
void initUserObjects (const std::vector< UserObjectName > &, unsigned int, std::vector< const T * > &)
 This function initializes user objects.
 
virtual void initJacobianVariables ()
 This function initialize variables required for Jacobian calculation.
 
virtual void initQpStatefulProperties ()
 Initializes state.
 
virtual void computeQpStress ()
 This function computes the Cauchy stress.
 
virtual void computeQpJacobian ()
 This function computes the Jacobian.
 
virtual void saveOldState ()
 This function saves the old stateful properties that is modified during sub stepping.
 
virtual void preSolveQp ()
 Sets state for solve.
 
virtual bool solveQp ()
 Solve state.
 
virtual void postSolveQp ()
 Update state for output (Outside substepping)
 
virtual void recoverOldState ()
 This function restores the the old stateful properties after a successful solve.
 
virtual void preSolveFlowrate ()
 Sets state for solve (Inside substepping)
 
virtual bool solveFlowrate ()
 Solve for flow rate and state.
 
virtual void postSolveFlowrate ()
 Update state for output (Inside substepping)
 
virtual bool computeFlowRateFunction ()
 Calls user objects to compute flow rates.
 
virtual bool computeFlowDirection ()
 Calls user objects to compute flow directions.
 
virtual void computeElasticRightCauchyGreenTensor ()
 Computes elastic Right Cauchy Green Tensor.
 
virtual void computePK2StressAndDerivative ()
 Computes PK2 stress and derivative w.r.t elastic Right Cauchy Green Tensor.
 
virtual void computeElasticStrain ()
 Computes elastic Lagrangian strain.
 
virtual void computeDeeDce ()
 Computes derivative of elastic strain w.r.t elastic Right Cauchy Green Tensor.
 
virtual bool computeFlowRateResidual ()
 Computes flow rate residual vector.
 
virtual void computeFlowRateJacobian ()
 Computes flow rate Jacobian matrix.
 
virtual void computeElasticPlasticDeformGrad ()
 Computes elastic and plastic deformation gradients.
 
virtual Real computeNorm (const std::vector< Real > &)
 Computes norm of residual vector.
 
virtual void updateFlowRate ()
 Update flow rate.
 
virtual void computeDpk2Dfpinv ()
 Computes derivative of PK2 stress wrt inverse of plastic deformation gradient.
 
virtual bool computeIntVarRates ()
 This function call user objects to calculate rate of internal variables.
 
virtual bool computeIntVar ()
 This function call user objects to integrate internal variables.
 
virtual bool computeStrength ()
 This function call user objects to compute strength.
 
virtual void computeIntVarRateDerivatives ()
 This function call user objects to compute dintvar_rate/dintvar and dintvarrate/dflowrate.
 
virtual void computeIntVarDerivatives ()
 This function call user objects to compute dintvar/dintvar_rate and dintvar/dflowrate.
 
void computeStrengthDerivatives ()
 This function call user objects to compute dstrength/dintvar.
 
virtual void computeQpProperties () override
 

Protected Attributes

Real _resid_abs_tol
 Absolute tolerance for residual convergence check.
 
Real _resid_rel_tol
 Relative tolerance for residual convergence check.
 
unsigned int _maxiters
 Maximum number of iterations.
 
unsigned int _max_substep_iter
 Maximum number of substep iterations.
 
std::vector< UserObjectName > _flow_rate_uo_names
 Names of flow rate user objects.
 
std::vector< UserObjectName > _strength_uo_names
 Names of strength user objects.
 
std::vector< UserObjectName > _int_var_uo_names
 Names of internal variable user objects.
 
std::vector< UserObjectName > _int_var_rate_uo_names
 Names of internal variable rate user objects.
 
unsigned int _num_flow_rate_uos
 Number of flow rate user objects.
 
unsigned int _num_strength_uos
 Number of strength user objects.
 
unsigned int _num_int_var_uos
 Number of internal variable user objects.
 
unsigned int _num_int_var_rate_uos
 Number of internal variable rate user objects.
 
std::vector< const HEVPFlowRateUOBase * > _flow_rate_uo
 Flow rate user objects.
 
std::vector< const HEVPStrengthUOBase * > _strength_uo
 Strength user objects.
 
std::vector< const HEVPInternalVarUOBase * > _int_var_uo
 Internal variable user objects.
 
std::vector< const HEVPInternalVarRateUOBase * > _int_var_rate_uo
 Internal variable rate user objects.
 
std::string _pk2_prop_name
 
MaterialProperty< RankTwoTensor > & _pk2
 
MaterialProperty< RankTwoTensor > & _fp
 
const MaterialProperty< RankTwoTensor > & _fp_old
 
MaterialProperty< RankTwoTensor > & _ce
 
const std::string _elasticity_tensor_name
 Name of the elasticity tensor material property.
 
const MaterialProperty< RankFourTensor > & _elasticity_tensor
 Elasticity tensor material property.
 
const MaterialProperty< RankTwoTensor > & _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.
 
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
 Base name prepended to all material property names to allow for multi-material systems.
 
const MaterialProperty< RankTwoTensor > & _mechanical_strain
 Mechanical strain material property.
 
MaterialProperty< RankTwoTensor > & _stress
 Stress material property.
 
MaterialProperty< RankTwoTensor > & _elastic_strain
 Elastic strain material property.
 
const MaterialProperty< RankTwoTensor > & _extra_stress
 Extra stress tensor.
 
std::vector< const Function * > _initial_stress_fcn
 initial stress components
 
MaterialProperty< RankFourTensor > & _Jacobian_mult
 derivative of stress w.r.t. strain (_dstress_dstrain)
 

Private Member Functions

bool haveMaterialProperty (const std::string &prop_name)
 
std::vector< VariableName > buildVariableVector (const VariableName &c1, const VariableName &c2, const VariableName &c3)
 
void validateCouplingHelper (const MaterialPropertyName &base, const std::vector< VariableName > &c, const System &system, std::vector< VariableName > &missing)
 
bool isNotObjectVariable (const VariableName &name)
 

Private Attributes

FEProblemBase & _dmi_fe_problem
 

Detailed Description

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 28 of file FiniteStrainHyperElasticViscoPlastic.h.

Constructor & Destructor Documentation

◆ FiniteStrainHyperElasticViscoPlastic()

FiniteStrainHyperElasticViscoPlastic::FiniteStrainHyperElasticViscoPlastic ( const InputParameters &  parameters)

Definition at line 43 of file FiniteStrainHyperElasticViscoPlastic.C.

45 : ComputeStressBase(parameters),
46 _resid_abs_tol(getParam<Real>("resid_abs_tol")),
47 _resid_rel_tol(getParam<Real>("resid_rel_tol")),
48 _maxiters(getParam<unsigned int>("maxiters")),
49 _max_substep_iter(getParam<unsigned int>("max_substep_iteration")),
50 _flow_rate_uo_names(isParamValid("flow_rate_user_objects")
51 ? getParam<std::vector<UserObjectName>>("flow_rate_user_objects")
52 : std::vector<UserObjectName>(0)),
53 _strength_uo_names(isParamValid("strength_user_objects")
54 ? getParam<std::vector<UserObjectName>>("strength_user_objects")
55 : std::vector<UserObjectName>(0)),
56 _int_var_uo_names(isParamValid("internal_var_user_objects")
57 ? getParam<std::vector<UserObjectName>>("internal_var_user_objects")
58 : std::vector<UserObjectName>(0)),
59 _int_var_rate_uo_names(
60 isParamValid("internal_var_rate_user_objects")
61 ? getParam<std::vector<UserObjectName>>("internal_var_rate_user_objects")
62 : std::vector<UserObjectName>(0)),
63 _pk2_prop_name(_base_name + "pk2"),
64 _pk2(declareProperty<RankTwoTensor>(_pk2_prop_name)),
65 _fp(declareProperty<RankTwoTensor>(_base_name + "fp")),
66 _fp_old(getMaterialPropertyOld<RankTwoTensor>(_base_name + "fp")),
67 _ce(declareProperty<RankTwoTensor>(_base_name + "ce")),
68 _elasticity_tensor_name(_base_name + "elasticity_tensor"),
69 _elasticity_tensor(getMaterialPropertyByName<RankFourTensor>(_elasticity_tensor_name)),
70 _deformation_gradient(getMaterialProperty<RankTwoTensor>(_base_name + "deformation_gradient")),
71 _deformation_gradient_old(
72 getMaterialPropertyOld<RankTwoTensor>(_base_name + "deformation_gradient")),
73 _rotation_increment(getMaterialProperty<RankTwoTensor>(_base_name + "rotation_increment"))
74{
75 initUOVariables();
76
77 initJacobianVariables();
78
79 _dflow_rate.resize(_num_flow_rate_uos);
80 _flow_rate.resize(_num_flow_rate_uos);
81 _resid.resize(_num_flow_rate_uos);
82
83 _flow_dirn.resize(_num_flow_rate_uos);
84}
ComputeGeneralStressBase::_base_name
const std::string _base_name
Base name prepended to all material property names to allow for multi-material systems.
Definition ComputeGeneralStressBase.h:43
ComputeStressBase
ComputeStressBase is the base class for stress tensors computed from MOOSE's strain calculators.
Definition ComputeStressBase.h:23
FiniteStrainHyperElasticViscoPlastic::_elasticity_tensor
const MaterialProperty< RankFourTensor > & _elasticity_tensor
Elasticity tensor material property.
Definition FiniteStrainHyperElasticViscoPlastic.h:199
FiniteStrainHyperElasticViscoPlastic::_flow_rate_uo_names
std::vector< UserObjectName > _flow_rate_uo_names
Names of flow rate user objects.
Definition FiniteStrainHyperElasticViscoPlastic.h:164
FiniteStrainHyperElasticViscoPlastic::initUOVariables
virtual void initUOVariables()
This function initializes the properties, stateful properties and user objects The properties and sta...
Definition FiniteStrainHyperElasticViscoPlastic.C:87
FiniteStrainHyperElasticViscoPlastic::_int_var_uo_names
std::vector< UserObjectName > _int_var_uo_names
Names of internal variable user objects.
Definition FiniteStrainHyperElasticViscoPlastic.h:168
FiniteStrainHyperElasticViscoPlastic::initJacobianVariables
virtual void initJacobianVariables()
This function initialize variables required for Jacobian calculation.
Definition FiniteStrainHyperElasticViscoPlastic.C:155
FiniteStrainHyperElasticViscoPlastic::_strength_uo_names
std::vector< UserObjectName > _strength_uo_names
Names of strength user objects.
Definition FiniteStrainHyperElasticViscoPlastic.h:166
FiniteStrainHyperElasticViscoPlastic::_rotation_increment
const MaterialProperty< RankTwoTensor > & _rotation_increment
Definition FiniteStrainHyperElasticViscoPlastic.h:202
FiniteStrainHyperElasticViscoPlastic::_fp_old
const MaterialProperty< RankTwoTensor > & _fp_old
Definition FiniteStrainHyperElasticViscoPlastic.h:193
FiniteStrainHyperElasticViscoPlastic::_deformation_gradient_old
const MaterialProperty< RankTwoTensor > & _deformation_gradient_old
Definition FiniteStrainHyperElasticViscoPlastic.h:201
FiniteStrainHyperElasticViscoPlastic::_num_flow_rate_uos
unsigned int _num_flow_rate_uos
Number of flow rate user objects.
Definition FiniteStrainHyperElasticViscoPlastic.h:173
FiniteStrainHyperElasticViscoPlastic::_resid_rel_tol
Real _resid_rel_tol
Relative tolerance for residual convergence check.
Definition FiniteStrainHyperElasticViscoPlastic.h:157
FiniteStrainHyperElasticViscoPlastic::_int_var_rate_uo_names
std::vector< UserObjectName > _int_var_rate_uo_names
Names of internal variable rate user objects.
Definition FiniteStrainHyperElasticViscoPlastic.h:170
FiniteStrainHyperElasticViscoPlastic::_maxiters
unsigned int _maxiters
Maximum number of iterations.
Definition FiniteStrainHyperElasticViscoPlastic.h:159
FiniteStrainHyperElasticViscoPlastic::_deformation_gradient
const MaterialProperty< RankTwoTensor > & _deformation_gradient
Definition FiniteStrainHyperElasticViscoPlastic.h:200
FiniteStrainHyperElasticViscoPlastic::_max_substep_iter
unsigned int _max_substep_iter
Maximum number of substep iterations.
Definition FiniteStrainHyperElasticViscoPlastic.h:161
FiniteStrainHyperElasticViscoPlastic::_elasticity_tensor_name
const std::string _elasticity_tensor_name
Name of the elasticity tensor material property.
Definition FiniteStrainHyperElasticViscoPlastic.h:197
FiniteStrainHyperElasticViscoPlastic::_resid_abs_tol
Real _resid_abs_tol
Absolute tolerance for residual convergence check.
Definition FiniteStrainHyperElasticViscoPlastic.h:155

Member Function Documentation

◆ computeDeeDce()

void FiniteStrainHyperElasticViscoPlastic::computeDeeDce ( )
protectedvirtual

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

Definition at line 487 of file FiniteStrainHyperElasticViscoPlastic.C.

488{
489 _dee_dce.zero();
490
491 for (const auto i : make_range(Moose::dim))
492 for (const auto j : make_range(Moose::dim))
493 _dee_dce(i, j, i, j) = 0.5;
494}
for
for(PetscInt i=0;i< nvars;++i)
dim
unsigned int dim
libMesh::make_range
IntRange< T > make_range(T beg, T end)

Referenced by initJacobianVariables().

◆ computeDpk2Dfpinv()

void FiniteStrainHyperElasticViscoPlastic::computeDpk2Dfpinv ( )
protectedvirtual

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

Definition at line 517 of file FiniteStrainHyperElasticViscoPlastic.C.

518{
519 for (const auto i : make_range(Moose::dim))
520 for (const auto j : make_range(Moose::dim))
521 for (const auto k : make_range(Moose::dim))
522 _dfe_dfpinv(i, j, k, j) = _dfgrd_tmp(i, k);
523
524 _dpk2_dfpinv = _dpk2_dce * _dce_dfe * _dfe_dfpinv;
525}

Referenced by computeFlowRateJacobian().

◆ computeElasticPlasticDeformGrad()

void FiniteStrainHyperElasticViscoPlastic::computeElasticPlasticDeformGrad ( )
protectedvirtual

Computes elastic and plastic deformation gradients.

Definition at line 503 of file FiniteStrainHyperElasticViscoPlastic.C.

504{
505 RankTwoTensor iden(RankTwoTensor::initIdentity);
506
507 RankTwoTensor val;
508 for (unsigned int i = 0; i < _num_flow_rate_uos; ++i)
509 val += _flow_rate(i) * _flow_dirn[i] * _dt_substep;
510
511 _fp_tmp_inv = _fp_tmp_old_inv * (iden - val);
512 _fp_tmp_inv = std::pow(_fp_tmp_inv.det(), -1.0 / 3.0) * _fp_tmp_inv;
513 _fe = _dfgrd_tmp * _fp_tmp_inv;
514}

Referenced by solveFlowrate().

◆ computeElasticRightCauchyGreenTensor()

void FiniteStrainHyperElasticViscoPlastic::computeElasticRightCauchyGreenTensor ( )
protectedvirtual

Computes elastic Right Cauchy Green Tensor.

Definition at line 497 of file FiniteStrainHyperElasticViscoPlastic.C.

498{
499 _ce[_qp] = _fe.transpose() * _fe;
500}
RankTwoTensorTempl::transpose
RankTwoTensorTempl< T > transpose() const

Referenced by computeFlowRateResidual().

◆ computeElasticStrain()

void FiniteStrainHyperElasticViscoPlastic::computeElasticStrain ( )
protectedvirtual

Computes elastic Lagrangian strain.

Definition at line 480 of file FiniteStrainHyperElasticViscoPlastic.C.

481{
482 RankTwoTensor iden(RankTwoTensor::initIdentity);
483 _ee = 0.5 * (_ce[_qp] - iden);
484}

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

◆ computeFlowDirection()

bool FiniteStrainHyperElasticViscoPlastic::computeFlowDirection ( )
protectedvirtual

Calls user objects to compute flow directions.

Definition at line 437 of file FiniteStrainHyperElasticViscoPlastic.C.

438{
439 for (unsigned i = 0; i < _num_flow_rate_uos; ++i)
440 {
441 if (!_flow_rate_uo[i]->computeDirection(_qp, _flow_dirn[i]))
442 return false;
443 }
444 return true;
445}
FiniteStrainHyperElasticViscoPlastic::_flow_rate_uo
std::vector< const HEVPFlowRateUOBase * > _flow_rate_uo
Flow rate user objects.
Definition FiniteStrainHyperElasticViscoPlastic.h:182

Referenced by computeFlowRateResidual().

◆ computeFlowRateFunction()

bool FiniteStrainHyperElasticViscoPlastic::computeFlowRateFunction ( )
protectedvirtual

Calls user objects to compute flow rates.

Definition at line 448 of file FiniteStrainHyperElasticViscoPlastic.C.

449{
450 Real val = 0;
451 for (unsigned i = 0; i < _num_flow_rate_uos; ++i)
452 {
453 if (_flow_rate_uo[i]->computeValue(_qp, val))
454 _resid(i) = -val;
455 else
456 return false;
457 }
458 return true;
459}
libMesh::Real
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real

Referenced by computeFlowRateResidual().

◆ computeFlowRateJacobian()

void FiniteStrainHyperElasticViscoPlastic::computeFlowRateJacobian ( )
protectedvirtual

Computes flow rate Jacobian matrix.

Definition at line 399 of file FiniteStrainHyperElasticViscoPlastic.C.

400{
401 computeIntVarRateDerivatives();
402 computeIntVarDerivatives();
403 computeStrengthDerivatives();
404
405 for (unsigned int i = 0; i < _num_flow_rate_uos; ++i)
406 for (unsigned int j = 0; j < _num_strength_uos; ++j)
407 _flow_rate_uo[i]->computeDerivative(_qp, _strength_uo_names[j], _dflowrate_dstrength(i, j));
408
409 for (unsigned int i = 0; i < _num_flow_rate_uos; ++i)
410 _flow_rate_uo[i]->computeTensorDerivative(_qp, _pk2_prop_name, _dflowrate_dpk2[i]);
411
412 computeDpk2Dfpinv();
413
414 for (unsigned int i = 0; i < _num_flow_rate_uos; ++i)
415 {
416 _dfpinv_dflowrate[i] = -_fp_tmp_old_inv * _flow_dirn[i] * _dt_substep;
417 _dpk2_dflowrate[i] = _dpk2_dfpinv * _dfpinv_dflowrate[i];
418 }
419
420 DenseMatrix<Real> dflowrate_dflowrate;
421 dflowrate_dflowrate = _dflowrate_dstrength;
422 dflowrate_dflowrate.right_multiply(_dstrength_dintvar);
423 dflowrate_dflowrate.right_multiply(_dintvar_dflowrate);
424
425 _jac.zero();
426 for (unsigned int i = 0; i < _num_flow_rate_uos; ++i)
427 for (unsigned int j = 0; j < _num_flow_rate_uos; ++j)
428 {
429 if (i == j)
430 _jac(i, j) = 1;
431 _jac(i, j) -= dflowrate_dflowrate(i, j);
432 _jac(i, j) -= _dflowrate_dpk2[i].doubleContraction(_dpk2_dflowrate[j]);
433 }
434}
FiniteStrainHyperElasticViscoPlastic::_num_strength_uos
unsigned int _num_strength_uos
Number of strength user objects.
Definition FiniteStrainHyperElasticViscoPlastic.h:175
FiniteStrainHyperElasticViscoPlastic::computeIntVarDerivatives
virtual void computeIntVarDerivatives()
This function call user objects to compute dintvar/dintvar_rate and dintvar/dflowrate.
Definition FiniteStrainHyperElasticViscoPlastic.C:635
FiniteStrainHyperElasticViscoPlastic::computeDpk2Dfpinv
virtual void computeDpk2Dfpinv()
Computes derivative of PK2 stress wrt inverse of plastic deformation gradient.
Definition FiniteStrainHyperElasticViscoPlastic.C:517
FiniteStrainHyperElasticViscoPlastic::computeIntVarRateDerivatives
virtual void computeIntVarRateDerivatives()
This function call user objects to compute dintvar_rate/dintvar and dintvarrate/dflowrate.
Definition FiniteStrainHyperElasticViscoPlastic.C:622
FiniteStrainHyperElasticViscoPlastic::computeStrengthDerivatives
void computeStrengthDerivatives()
This function call user objects to compute dstrength/dintvar.
Definition FiniteStrainHyperElasticViscoPlastic.C:670
libMesh::DenseMatrix::right_multiply
virtual void right_multiply(const DenseMatrixBase< T > &M2) override final
libMesh::DenseMatrix::zero
virtual void zero() override final

Referenced by solveFlowrate().

◆ computeFlowRateResidual()

bool FiniteStrainHyperElasticViscoPlastic::computeFlowRateResidual ( )
protectedvirtual

Computes flow rate residual vector.

Definition at line 373 of file FiniteStrainHyperElasticViscoPlastic.C.

374{
375 if (!computeIntVarRates())
376 return false;
377
378 if (!computeIntVar())
379 return false;
380
381 if (!computeStrength())
382 return false;
383
384 computeElasticRightCauchyGreenTensor();
385 computePK2StressAndDerivative();
386
387 if (!computeFlowRateFunction())
388 return false;
389
390 if (!computeFlowDirection())
391 return false;
392
393 _resid += _flow_rate;
394
395 return true;
396}
FiniteStrainHyperElasticViscoPlastic::computeElasticRightCauchyGreenTensor
virtual void computeElasticRightCauchyGreenTensor()
Computes elastic Right Cauchy Green Tensor.
Definition FiniteStrainHyperElasticViscoPlastic.C:497
FiniteStrainHyperElasticViscoPlastic::computeIntVarRates
virtual bool computeIntVarRates()
This function call user objects to calculate rate of internal variables.
Definition FiniteStrainHyperElasticViscoPlastic.C:580
FiniteStrainHyperElasticViscoPlastic::computeFlowDirection
virtual bool computeFlowDirection()
Calls user objects to compute flow directions.
Definition FiniteStrainHyperElasticViscoPlastic.C:437
FiniteStrainHyperElasticViscoPlastic::computeStrength
virtual bool computeStrength()
This function call user objects to compute strength.
Definition FiniteStrainHyperElasticViscoPlastic.C:608
FiniteStrainHyperElasticViscoPlastic::computeFlowRateFunction
virtual bool computeFlowRateFunction()
Calls user objects to compute flow rates.
Definition FiniteStrainHyperElasticViscoPlastic.C:448
FiniteStrainHyperElasticViscoPlastic::computeIntVar
virtual bool computeIntVar()
This function call user objects to integrate internal variables.
Definition FiniteStrainHyperElasticViscoPlastic.C:594
FiniteStrainHyperElasticViscoPlastic::computePK2StressAndDerivative
virtual void computePK2StressAndDerivative()
Computes PK2 stress and derivative w.r.t elastic Right Cauchy Green Tensor.
Definition FiniteStrainHyperElasticViscoPlastic.C:462

Referenced by solveFlowrate().

◆ computeIntVar()

bool FiniteStrainHyperElasticViscoPlastic::computeIntVar ( )
protectedvirtual

This function call user objects to integrate internal variables.

Definition at line 594 of file FiniteStrainHyperElasticViscoPlastic.C.

595{
596 Real val = 0;
597 for (unsigned int i = 0; i < _num_int_var_uos; ++i)
598 {
599 if (_int_var_uo[i]->computeValue(_qp, _dt_substep, val))
600 (*_int_var_stateful_prop[i])[_qp] = val;
601 else
602 return false;
603 }
604 return true;
605}
FiniteStrainHyperElasticViscoPlastic::_int_var_stateful_prop
std::vector< MaterialProperty< Real > * > _int_var_stateful_prop
Definition FiniteStrainHyperElasticViscoPlastic.h:206
FiniteStrainHyperElasticViscoPlastic::_num_int_var_uos
unsigned int _num_int_var_uos
Number of internal variable user objects.
Definition FiniteStrainHyperElasticViscoPlastic.h:177
FiniteStrainHyperElasticViscoPlastic::_int_var_uo
std::vector< const HEVPInternalVarUOBase * > _int_var_uo
Internal variable user objects.
Definition FiniteStrainHyperElasticViscoPlastic.h:186

Referenced by computeFlowRateResidual().

◆ computeIntVarDerivatives()

void FiniteStrainHyperElasticViscoPlastic::computeIntVarDerivatives ( )
protectedvirtual

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

Definition at line 635 of file FiniteStrainHyperElasticViscoPlastic.C.

636{
637 Real val = 0;
638
639 for (unsigned int i = 0; i < _num_int_var_uos; ++i)
640 for (unsigned int j = 0; j < _num_int_var_rate_uos; ++j)
641 {
642 _int_var_uo[i]->computeDerivative(_qp, _dt_substep, _int_var_rate_uo_names[j], val);
643 _dintvar_dintvarrate(i, j) = val;
644 }
645
646 _dintvar_dintvar.zero();
647
648 for (unsigned int i = 0; i < _num_int_var_uos; ++i)
649 for (unsigned int j = 0; j < _num_int_var_uos; ++j)
650 {
651 if (i == j)
652 _dintvar_dintvar(i, j) = 1;
653 for (unsigned int k = 0; k < _num_int_var_rate_uos; ++k)
654 _dintvar_dintvar(i, j) -= _dintvar_dintvarrate(i, k) * _dintvarrate_dintvar(k, j);
655 }
656
657 for (unsigned int i = 0; i < _num_flow_rate_uos; ++i)
658 _dintvar_dintvarrate.vector_mult(_dintvar_dflowrate_tmp[i], _dintvarrate_dflowrate[i]);
659
660 for (unsigned int i = 0; i < _num_flow_rate_uos; ++i)
661 {
662 _dintvar_dintvar_x.zero();
663 _dintvar_dintvar.lu_solve(_dintvar_dflowrate_tmp[i], _dintvar_dintvar_x);
664 for (unsigned int j = 0; j < _num_int_var_uos; ++j)
665 _dintvar_dflowrate(j, i) = _dintvar_dintvar_x(j);
666 }
667}
FiniteStrainHyperElasticViscoPlastic::_dintvarrate_dflowrate
std::vector< DenseVector< Real > > _dintvarrate_dflowrate
Jacobian variables.
Definition FiniteStrainHyperElasticViscoPlastic.h:229
FiniteStrainHyperElasticViscoPlastic::_num_int_var_rate_uos
unsigned int _num_int_var_rate_uos
Number of internal variable rate user objects.
Definition FiniteStrainHyperElasticViscoPlastic.h:179
FiniteStrainHyperElasticViscoPlastic::_dintvar_dflowrate_tmp
std::vector< DenseVector< Real > > _dintvar_dflowrate_tmp
Definition FiniteStrainHyperElasticViscoPlastic.h:230
libMesh::DenseMatrix::vector_mult
void vector_mult(DenseVector< T > &dest, const DenseVector< T > &arg) const
libMesh::DenseMatrix::lu_solve
void lu_solve(const DenseVector< T > &b, DenseVector< T > &x)

Referenced by computeFlowRateJacobian().

◆ computeIntVarRateDerivatives()

void FiniteStrainHyperElasticViscoPlastic::computeIntVarRateDerivatives ( )
protectedvirtual

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

Definition at line 622 of file FiniteStrainHyperElasticViscoPlastic.C.

623{
624 Real val = 0;
625
626 for (unsigned int i = 0; i < _num_int_var_rate_uos; ++i)
627 for (unsigned int j = 0; j < _num_flow_rate_uos; ++j)
628 {
629 _int_var_rate_uo[i]->computeDerivative(_qp, _flow_rate_uo_names[j], val);
630 _dintvarrate_dflowrate[j](i) = val;
631 }
632}
FiniteStrainHyperElasticViscoPlastic::_int_var_rate_uo
std::vector< const HEVPInternalVarRateUOBase * > _int_var_rate_uo
Internal variable rate user objects.
Definition FiniteStrainHyperElasticViscoPlastic.h:188

Referenced by computeFlowRateJacobian().

◆ computeIntVarRates()

bool FiniteStrainHyperElasticViscoPlastic::computeIntVarRates ( )
protectedvirtual

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

Definition at line 580 of file FiniteStrainHyperElasticViscoPlastic.C.

581{
582 Real val = 0;
583 for (unsigned int i = 0; i < _num_int_var_rate_uos; ++i)
584 {
585 if (_int_var_rate_uo[i]->computeValue(_qp, val))
586 (*_int_var_rate_prop[i])[_qp] = val;
587 else
588 return false;
589 }
590 return true;
591}
FiniteStrainHyperElasticViscoPlastic::_int_var_rate_prop
std::vector< MaterialProperty< Real > * > _int_var_rate_prop
Definition FiniteStrainHyperElasticViscoPlastic.h:208

Referenced by computeFlowRateResidual().

◆ computeNorm()

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

Computes norm of residual vector.

Definition at line 528 of file FiniteStrainHyperElasticViscoPlastic.C.

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

Referenced by solveFlowrate().

◆ computePK2StressAndDerivative()

void FiniteStrainHyperElasticViscoPlastic::computePK2StressAndDerivative ( )
protectedvirtual

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

Reimplemented in HyperElasticPhaseFieldIsoDamage.

Definition at line 462 of file FiniteStrainHyperElasticViscoPlastic.C.

463{
464 computeElasticStrain();
465 _pk2[_qp] = _elasticity_tensor[_qp] * _ee;
466
467 _dce_dfe.zero();
468 for (const auto i : make_range(Moose::dim))
469 for (const auto j : make_range(Moose::dim))
470 for (const auto k : make_range(Moose::dim))
471 {
472 _dce_dfe(i, j, k, i) = _dce_dfe(i, j, k, i) + _fe(k, j);
473 _dce_dfe(i, j, k, j) = _dce_dfe(i, j, k, j) + _fe(k, i);
474 }
475
476 _dpk2_dfe = _dpk2_dce * _dce_dfe;
477}
FiniteStrainHyperElasticViscoPlastic::computeElasticStrain
virtual void computeElasticStrain()
Computes elastic Lagrangian strain.
Definition FiniteStrainHyperElasticViscoPlastic.C:480

Referenced by computeFlowRateResidual().

◆ computeQpJacobian()

void FiniteStrainHyperElasticViscoPlastic::computeQpJacobian ( )
protectedvirtual

This function computes the Jacobian.

Reimplemented in HyperElasticPhaseFieldIsoDamage.

Definition at line 547 of file FiniteStrainHyperElasticViscoPlastic.C.

548{
549 usingTensorIndices(i_, j_, k_, l_);
550 _tan_mod = _fe.times<i_, k_, j_, l_>(_fe) * _dpk2_dfe;
551 _pk2_fet = _pk2[_qp] * _fe.transpose();
552 _fe_pk2 = _fe * _pk2[_qp];
553
554 for (const auto i : make_range(Moose::dim))
555 for (const auto j : make_range(Moose::dim))
556 for (const auto l : make_range(Moose::dim))
557 {
558 _tan_mod(i, j, i, l) += _pk2_fet(l, j);
559 _tan_mod(i, j, j, l) += _fe_pk2(i, l);
560 }
561
562 _tan_mod /= _fe.det();
563
564 for (const auto i : make_range(Moose::dim))
565 for (const auto j : make_range(Moose::dim))
566 for (const auto l : make_range(Moose::dim))
567 _dfe_df(i, j, i, l) = _fp_tmp_inv(l, j);
568
569 for (const auto i : make_range(Moose::dim))
570 for (const auto j : make_range(Moose::dim))
571 for (const auto k : make_range(Moose::dim))
572 for (const auto l : make_range(Moose::dim))
573 _df_dstretch_inc(i, j, k, l) =
574 _rotation_increment[_qp](i, k) * _deformation_gradient_old[_qp](l, j);
575
576 _Jacobian_mult[_qp] = _tan_mod * _dfe_df * _df_dstretch_inc;
577}
ComputeGeneralStressBase::_Jacobian_mult
MaterialProperty< RankFourTensor > & _Jacobian_mult
derivative of stress w.r.t. strain (_dstress_dstrain)
Definition ComputeGeneralStressBase.h:59
RankTwoTensorTempl::times
RankFourTensorTempl< T > times(const RankTwoTensorTempl< T > &b) const

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

◆ computeQpProperties()

void ComputeGeneralStressBase::computeQpProperties ( )
overrideprotectedvirtualinherited

Definition at line 44 of file ComputeGeneralStressBase.C.

45{
46 computeQpStress();
47
48 // Add in extra stress
49 _stress[_qp] += _extra_stress[_qp];
50}
ComputeGeneralStressBase::_extra_stress
const MaterialProperty< RankTwoTensor > & _extra_stress
Extra stress tensor.
Definition ComputeGeneralStressBase.h:53
ComputeGeneralStressBase::computeQpStress
virtual void computeQpStress()=0
Compute the stress and store it in the _stress material property for the current quadrature point.
ComputeGeneralStressBase::_stress
MaterialProperty< RankTwoTensor > & _stress
Stress material property.
Definition ComputeGeneralStressBase.h:48

◆ computeQpStress()

void FiniteStrainHyperElasticViscoPlastic::computeQpStress ( )
protectedvirtual

This function computes the Cauchy stress.

Implements ComputeGeneralStressBase.

Definition at line 208 of file FiniteStrainHyperElasticViscoPlastic.C.

209{
210 bool converge;
211 RankTwoTensor delta_dfgrd = _deformation_gradient[_qp] - _deformation_gradient_old[_qp];
212 unsigned int num_substep = 1;
213 unsigned int substep_iter = 1;
214
215 saveOldState();
216
217 do
218 {
219 preSolveQp();
220
221 converge = true;
222 _dt_substep = _dt / num_substep;
223
224 for (unsigned int istep = 0; istep < num_substep; ++istep)
225 {
226 _dfgrd_tmp = (istep + 1.0) * delta_dfgrd / num_substep + _deformation_gradient_old[_qp];
227 if (!solveQp())
228 {
229 converge = false;
230 substep_iter++;
231 num_substep *= 2;
232 break;
233 }
234 }
235
236 if (substep_iter > _max_substep_iter)
237 mooseError("Constitutive failure with substepping at quadrature point ",
238 _q_point[_qp](0),
239 " ",
240 _q_point[_qp](1),
241 " ",
242 _q_point[_qp](2));
243 } while (!converge);
244
245 postSolveQp();
246}
mooseError
void mooseError(Args &&... args)
FiniteStrainHyperElasticViscoPlastic::postSolveQp
virtual void postSolveQp()
Update state for output (Outside substepping)
Definition FiniteStrainHyperElasticViscoPlastic.C:280
FiniteStrainHyperElasticViscoPlastic::saveOldState
virtual void saveOldState()
This function saves the old stateful properties that is modified during sub stepping.
Definition FiniteStrainHyperElasticViscoPlastic.C:249

◆ computeStrength()

bool FiniteStrainHyperElasticViscoPlastic::computeStrength ( )
protectedvirtual

This function call user objects to compute strength.

Definition at line 608 of file FiniteStrainHyperElasticViscoPlastic.C.

609{
610 Real val = 0;
611 for (unsigned int i = 0; i < _num_strength_uos; ++i)
612 {
613 if (_strength_uo[i]->computeValue(_qp, val))
614 (*_strength_prop[i])[_qp] = val;
615 else
616 return false;
617 }
618 return true;
619}
FiniteStrainHyperElasticViscoPlastic::_strength_prop
std::vector< MaterialProperty< Real > * > _strength_prop
Definition FiniteStrainHyperElasticViscoPlastic.h:205
FiniteStrainHyperElasticViscoPlastic::_strength_uo
std::vector< const HEVPStrengthUOBase * > _strength_uo
Strength user objects.
Definition FiniteStrainHyperElasticViscoPlastic.h:184

Referenced by computeFlowRateResidual().

◆ computeStrengthDerivatives()

void FiniteStrainHyperElasticViscoPlastic::computeStrengthDerivatives ( )
protected

This function call user objects to compute dstrength/dintvar.

Definition at line 670 of file FiniteStrainHyperElasticViscoPlastic.C.

671{
672 Real val = 0;
673
674 for (unsigned int i = 0; i < _num_strength_uos; ++i)
675 for (unsigned int j = 0; j < _num_int_var_uos; ++j)
676 {
677 _strength_uo[i]->computeDerivative(_qp, _int_var_uo_names[j], val);
678 _dstrength_dintvar(i, j) = val;
679 }
680}

Referenced by computeFlowRateJacobian().

◆ initJacobianVariables()

void FiniteStrainHyperElasticViscoPlastic::initJacobianVariables ( )
protectedvirtual

This function initialize variables required for Jacobian calculation.

Definition at line 155 of file FiniteStrainHyperElasticViscoPlastic.C.

156{
157 _dintvarrate_dflowrate.resize(_num_flow_rate_uos);
158 for (unsigned int i = 0; i < _num_flow_rate_uos; ++i)
159 _dintvarrate_dflowrate[i].resize(_num_int_var_rate_uos);
160
161 _dintvar_dflowrate_tmp.resize(_num_flow_rate_uos);
162 for (unsigned int i = 0; i < _num_flow_rate_uos; ++i)
163 _dintvar_dflowrate_tmp[i].resize(_num_int_var_uos);
164
165 _dintvarrate_dintvar.resize(_num_int_var_rate_uos, _num_int_var_uos);
166 _dintvar_dintvarrate.resize(_num_int_var_uos, _num_int_var_rate_uos);
167 _dintvar_dflowrate.resize(_num_int_var_uos, _num_flow_rate_uos);
168 _dintvar_dintvar.resize(_num_int_var_uos, _num_int_var_uos);
169 _dstrength_dintvar.resize(_num_strength_uos, _num_int_var_uos);
170 _dflowrate_dstrength.resize(_num_flow_rate_uos, _num_strength_uos);
171 _dintvar_dintvar_x.resize(_num_int_var_uos);
172
173 _dpk2_dflowrate.resize(_num_flow_rate_uos);
174 _dflowrate_dpk2.resize(_num_flow_rate_uos);
175 _dfpinv_dflowrate.resize(_num_flow_rate_uos);
176
177 _jac.resize(_num_flow_rate_uos, _num_flow_rate_uos);
178
179 computeDeeDce();
180}
FiniteStrainHyperElasticViscoPlastic::computeDeeDce
virtual void computeDeeDce()
Computes derivative of elastic strain w.r.t elastic Right Cauchy Green Tensor.
Definition FiniteStrainHyperElasticViscoPlastic.C:487
libMesh::DenseMatrix::resize
void resize(const unsigned int new_m, const unsigned int new_n)

Referenced by FiniteStrainHyperElasticViscoPlastic().

◆ 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 110 of file FiniteStrainHyperElasticViscoPlastic.C.

112{
113 uo_num = uo_names.size();
114}

Referenced by initUOVariables().

◆ 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 118 of file FiniteStrainHyperElasticViscoPlastic.C.

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

Referenced by initUOVariables().

◆ 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 129 of file FiniteStrainHyperElasticViscoPlastic.C.

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

Referenced by initUOVariables().

◆ initQpStatefulProperties()

void FiniteStrainHyperElasticViscoPlastic::initQpStatefulProperties ( )
protectedvirtual

Initializes state.

Reimplemented from ComputeGeneralStressBase.

Definition at line 183 of file FiniteStrainHyperElasticViscoPlastic.C.

184{
185 _stress[_qp].zero();
186 _ce[_qp].zero();
187 _pk2[_qp].zero();
188 _fp[_qp].setToIdentity();
189
190 for (unsigned int i = 0; i < _num_flow_rate_uos; ++i)
191 (*_flow_rate_prop[i])[_qp] = 0.0;
192
193 for (unsigned int i = 0; i < _num_strength_uos; ++i)
194 (*_strength_prop[i])[_qp] = 0.0;
195
196 for (unsigned int i = 0; i < _num_int_var_uos; ++i)
197 {
198 (*_int_var_stateful_prop[i])[_qp] = 0.0;
199 // TODO: remove this nasty const_cast if you can figure out how
200 const_cast<MaterialProperty<Real> &>(*_int_var_stateful_prop_old[i])[_qp] = 0.0;
201 }
202
203 for (unsigned int i = 0; i < _num_int_var_rate_uos; ++i)
204 (*_int_var_rate_prop[i])[_qp] = 0.0;
205}
FiniteStrainHyperElasticViscoPlastic::_flow_rate_prop
std::vector< MaterialProperty< Real > * > _flow_rate_prop
Definition FiniteStrainHyperElasticViscoPlastic.h:204
FiniteStrainHyperElasticViscoPlastic::_int_var_stateful_prop_old
std::vector< const MaterialProperty< Real > * > _int_var_stateful_prop_old
Definition FiniteStrainHyperElasticViscoPlastic.h:207

◆ 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 87 of file FiniteStrainHyperElasticViscoPlastic.C.

88{
89 initNumUserObjects(_flow_rate_uo_names, _num_flow_rate_uos);
90 initNumUserObjects(_strength_uo_names, _num_strength_uos);
91 initNumUserObjects(_int_var_uo_names, _num_int_var_uos);
92 initNumUserObjects(_int_var_rate_uo_names, _num_int_var_rate_uos);
93
94 initProp(_flow_rate_uo_names, _num_flow_rate_uos, _flow_rate_prop);
95 initProp(_strength_uo_names, _num_strength_uos, _strength_prop);
96 initProp(_int_var_uo_names, _num_int_var_uos, _int_var_stateful_prop);
97 initProp(_int_var_rate_uo_names, _num_int_var_rate_uos, _int_var_rate_prop);
98
99 initPropOld(_int_var_uo_names, _num_int_var_uos, _int_var_stateful_prop_old);
100
101 initUserObjects(_flow_rate_uo_names, _num_flow_rate_uos, _flow_rate_uo);
102 initUserObjects(_strength_uo_names, _num_strength_uos, _strength_uo);
103 initUserObjects(_int_var_uo_names, _num_int_var_uos, _int_var_uo);
104 initUserObjects(_int_var_rate_uo_names, _num_int_var_rate_uos, _int_var_rate_uo);
105
106 _int_var_old.resize(_num_int_var_uos, 0.0);
107}
FiniteStrainHyperElasticViscoPlastic::initNumUserObjects
void initNumUserObjects(const std::vector< UserObjectName > &, unsigned int &)
This function calculates the number of each user object type.
Definition FiniteStrainHyperElasticViscoPlastic.C:110
FiniteStrainHyperElasticViscoPlastic::initUserObjects
void initUserObjects(const std::vector< UserObjectName > &, unsigned int, std::vector< const T * > &)
This function initializes user objects.
Definition FiniteStrainHyperElasticViscoPlastic.C:141
FiniteStrainHyperElasticViscoPlastic::initProp
void initProp(const std::vector< UserObjectName > &, unsigned int, std::vector< MaterialProperty< T > * > &)
This function initializes properties for each user object.
Definition FiniteStrainHyperElasticViscoPlastic.C:118
FiniteStrainHyperElasticViscoPlastic::initPropOld
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 t...
Definition FiniteStrainHyperElasticViscoPlastic.C:129

Referenced by FiniteStrainHyperElasticViscoPlastic().

◆ 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 141 of file FiniteStrainHyperElasticViscoPlastic.C.

144{
145 uo.resize(uo_num);
146
147 if (uo_num == 0)
148 mooseError("Specify atleast one user object of type", typeid(T).name());
149
150 for (unsigned int i = 0; i < uo_num; ++i)
151 uo[i] = &getUserObjectByName<T>(uo_names[i]);
152}
T
const double T
Definition GasLiquidMassTransferTest.C:19
name
const std::string name
Definition Setup.h:21

Referenced by initUOVariables().

◆ postSolveFlowrate()

void FiniteStrainHyperElasticViscoPlastic::postSolveFlowrate ( )
protectedvirtual

Update state for output (Inside substepping)

Definition at line 362 of file FiniteStrainHyperElasticViscoPlastic.C.

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

Referenced by solveQp().

◆ postSolveQp()

void FiniteStrainHyperElasticViscoPlastic::postSolveQp ( )
protectedvirtual

Update state for output (Outside substepping)

Definition at line 280 of file FiniteStrainHyperElasticViscoPlastic.C.

281{
282 recoverOldState();
283
284 _stress[_qp] = _fe * _pk2[_qp] * _fe.transpose() / _fe.det();
285 _fp[_qp] = _fp_tmp_inv.inverse();
286
287 computeQpJacobian();
288}
FiniteStrainHyperElasticViscoPlastic::recoverOldState
virtual void recoverOldState()
This function restores the the old stateful properties after a successful solve.
Definition FiniteStrainHyperElasticViscoPlastic.C:291
FiniteStrainHyperElasticViscoPlastic::computeQpJacobian
virtual void computeQpJacobian()
This function computes the Jacobian.
Definition FiniteStrainHyperElasticViscoPlastic.C:547
RankTwoTensorTempl::inverse
RankTwoTensorTempl< T > inverse() const

Referenced by computeQpStress().

◆ preSolveFlowrate()

void FiniteStrainHyperElasticViscoPlastic::preSolveFlowrate ( )
protectedvirtual

Sets state for solve (Inside substepping)

Definition at line 299 of file FiniteStrainHyperElasticViscoPlastic.C.

300{
301 for (unsigned int i = 0; i < _num_flow_rate_uos; ++i)
302 {
303 _flow_rate(i) = 0.0;
304 (*_flow_rate_prop[i])[_qp] = 0.0;
305 }
306
307 for (unsigned int i = 0; i < _num_int_var_uos; ++i)
308 (*_int_var_stateful_prop[i])[_qp] = (*_int_var_stateful_prop_old[i])[_qp];
309
310 _fp_tmp_inv = _fp_tmp_old_inv;
311 _fe = _dfgrd_tmp * _fp_tmp_inv;
312}

Referenced by solveQp().

◆ preSolveQp()

void FiniteStrainHyperElasticViscoPlastic::preSolveQp ( )
protectedvirtual

Sets state for solve.

Definition at line 256 of file FiniteStrainHyperElasticViscoPlastic.C.

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

Referenced by computeQpStress().

◆ recoverOldState()

void FiniteStrainHyperElasticViscoPlastic::recoverOldState ( )
protectedvirtual

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

Definition at line 291 of file FiniteStrainHyperElasticViscoPlastic.C.

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

Referenced by postSolveQp().

◆ saveOldState()

void FiniteStrainHyperElasticViscoPlastic::saveOldState ( )
protectedvirtual

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

Definition at line 249 of file FiniteStrainHyperElasticViscoPlastic.C.

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

Referenced by computeQpStress().

◆ solveFlowrate()

bool FiniteStrainHyperElasticViscoPlastic::solveFlowrate ( )
protectedvirtual

Solve for flow rate and state.

Definition at line 315 of file FiniteStrainHyperElasticViscoPlastic.C.

316{
317 Real resid0, rnorm;
318 unsigned int iter = 0;
319
320#ifdef DEBUG
321 std::vector<Real> rnormst(_maxiters + 1), flowratest(_maxiters + 1);
322#endif
323
324 if (!computeFlowRateResidual())
325 return false;
326
327 rnorm = computeNorm(_resid.get_values());
328 resid0 = rnorm;
329
330#ifdef DEBUG
331 rnormst[iter] = rnorm;
332 flowratest[iter] = computeNorm(_flow_rate.get_values());
333#endif
334
335 while (rnorm > _resid_abs_tol && rnorm > _resid_rel_tol * resid0 && iter < _maxiters)
336 {
337 computeFlowRateJacobian();
338
339 updateFlowRate();
340
341 computeElasticPlasticDeformGrad();
342
343 if (!computeFlowRateResidual())
344 return false;
345
346 rnorm = computeNorm(_resid.get_values());
347 iter++;
348
349#ifdef DEBUG
350 rnormst[iter] = rnorm;
351 flowratest[iter] = computeNorm(_flow_rate.get_values());
352#endif
353 }
354
355 if (iter == _maxiters && rnorm > _resid_abs_tol && rnorm > _resid_rel_tol * resid0)
356 return false;
357
358 return true;
359}
FiniteStrainHyperElasticViscoPlastic::computeElasticPlasticDeformGrad
virtual void computeElasticPlasticDeformGrad()
Computes elastic and plastic deformation gradients.
Definition FiniteStrainHyperElasticViscoPlastic.C:503
FiniteStrainHyperElasticViscoPlastic::computeFlowRateResidual
virtual bool computeFlowRateResidual()
Computes flow rate residual vector.
Definition FiniteStrainHyperElasticViscoPlastic.C:373
FiniteStrainHyperElasticViscoPlastic::computeFlowRateJacobian
virtual void computeFlowRateJacobian()
Computes flow rate Jacobian matrix.
Definition FiniteStrainHyperElasticViscoPlastic.C:399
FiniteStrainHyperElasticViscoPlastic::computeNorm
virtual Real computeNorm(const std::vector< Real > &)
Computes norm of residual vector.
Definition FiniteStrainHyperElasticViscoPlastic.C:528

Referenced by solveQp().

◆ solveQp()

bool FiniteStrainHyperElasticViscoPlastic::solveQp ( )
protectedvirtual

Solve state.

Definition at line 269 of file FiniteStrainHyperElasticViscoPlastic.C.

270{
271 preSolveFlowrate();
272 if (!solveFlowrate())
273 return false;
274 postSolveFlowrate();
275
276 return true;
277}
FiniteStrainHyperElasticViscoPlastic::preSolveFlowrate
virtual void preSolveFlowrate()
Sets state for solve (Inside substepping)
Definition FiniteStrainHyperElasticViscoPlastic.C:299
FiniteStrainHyperElasticViscoPlastic::postSolveFlowrate
virtual void postSolveFlowrate()
Update state for output (Inside substepping)
Definition FiniteStrainHyperElasticViscoPlastic.C:362
FiniteStrainHyperElasticViscoPlastic::solveFlowrate
virtual bool solveFlowrate()
Solve for flow rate and state.
Definition FiniteStrainHyperElasticViscoPlastic.C:315

Referenced by computeQpStress().

◆ updateFlowRate()

void FiniteStrainHyperElasticViscoPlastic::updateFlowRate ( )
protectedvirtual

Update flow rate.

Definition at line 537 of file FiniteStrainHyperElasticViscoPlastic.C.

538{
539 _jac.lu_solve(_resid, _dflow_rate);
540 _flow_rate -= _dflow_rate;
541
542 for (unsigned int i = 0; i < _num_flow_rate_uos; ++i)
543 (*_flow_rate_prop[i])[_qp] = _flow_rate(i);
544}

Referenced by solveFlowrate().

◆ validParams()

InputParameters FiniteStrainHyperElasticViscoPlastic::validParams ( )
static

Definition at line 16 of file FiniteStrainHyperElasticViscoPlastic.C.

17{
18 InputParameters params = ComputeStressBase::validParams();
19 params.addParam<Real>(
20 "resid_abs_tol", 1e-10, "Absolute Tolerance for flow rate residual equation");
21 params.addParam<Real>(
22 "resid_rel_tol", 1e-6, "Relative Tolerance for flow rate residual equation");
23 params.addParam<unsigned int>("maxiters", 50, "Maximum iteration for flow rate update");
24 params.addParam<unsigned int>("max_substep_iteration", 1, "Maximum number of substep iteration");
25 params.addParam<std::vector<UserObjectName>>(
26 "flow_rate_user_objects",
27 "List of User object names that computes flow rate and derivatives");
28 params.addParam<std::vector<UserObjectName>>(
29 "strength_user_objects",
30 "List of User object names that computes strength variables and derivatives");
31 params.addParam<std::vector<UserObjectName>>(
32 "internal_var_user_objects",
33 "List of User object names that integrates internal variables and computes derivatives");
34 params.addParam<std::vector<UserObjectName>>(
35 "internal_var_rate_user_objects",
36 "List of User object names that computes internal variable rates and derivatives");
37 params.addClassDescription("Material class for hyper-elastic viscoplatic flow: Can handle "
38 "multiple flow models defined by flowratemodel type user objects");
39
40 return params;
41}
ComputeStressBase::validParams
static InputParameters validParams()
Definition ComputeStressBase.C:15
InputParameters::addParam
void addParam(const std::string &name, const std::initializer_list< typename T::value_type > &value, const std::string &doc_string)
InputParameters::addClassDescription
void addClassDescription(const std::string &doc_string)

Referenced by HyperElasticPhaseFieldIsoDamage::validParams().

Member Data Documentation

◆ _base_name

const std::string ComputeGeneralStressBase::_base_name
protectedinherited

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

Definition at line 43 of file ComputeGeneralStressBase.h.

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

◆ _ce

MaterialProperty<RankTwoTensor>& FiniteStrainHyperElasticViscoPlastic::_ce
protected

Definition at line 194 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(), computePK2StressAndDerivative(), HyperElasticPhaseFieldIsoDamage::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 computeQpJacobian(), and HyperElasticPhaseFieldIsoDamage::computeQpJacobian().

◆ _dfe_df

RankFourTensor FiniteStrainHyperElasticViscoPlastic::_dfe_df
protected

Definition at line 216 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeQpJacobian(), and HyperElasticPhaseFieldIsoDamage::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(), computePK2StressAndDerivative(), HyperElasticPhaseFieldIsoDamage::computePK2StressAndDerivative(), and preSolveQp().

◆ _dpk2_dfe

RankFourTensor FiniteStrainHyperElasticViscoPlastic::_dpk2_dfe
protected

Definition at line 215 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computePK2StressAndDerivative(), HyperElasticPhaseFieldIsoDamage::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>& ComputeGeneralStressBase::_elastic_strain
protectedinherited

Elastic strain material property.

Definition at line 50 of file ComputeGeneralStressBase.h.

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

◆ _elasticity_tensor

const MaterialProperty<RankFourTensor>& FiniteStrainHyperElasticViscoPlastic::_elasticity_tensor
protected

Elasticity tensor material property.

Definition at line 199 of file FiniteStrainHyperElasticViscoPlastic.h.

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

◆ _elasticity_tensor_name

const std::string FiniteStrainHyperElasticViscoPlastic::_elasticity_tensor_name
protected

Name of the elasticity tensor material property.

Definition at line 197 of file FiniteStrainHyperElasticViscoPlastic.h.

◆ _extra_stress

const MaterialProperty<RankTwoTensor>& ComputeGeneralStressBase::_extra_stress
protectedinherited

Extra stress tensor.

Definition at line 53 of file ComputeGeneralStressBase.h.

Referenced by ComputeGeneralStressBase::computeQpProperties().

◆ _fe

RankTwoTensor FiniteStrainHyperElasticViscoPlastic::_fe
protected

Definition at line 213 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeElasticPlasticDeformGrad(), computeElasticRightCauchyGreenTensor(), computePK2StressAndDerivative(), HyperElasticPhaseFieldIsoDamage::computePK2StressAndDerivative(), computeQpJacobian(), HyperElasticPhaseFieldIsoDamage::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 182 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 164 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeIntVarRateDerivatives(), and initUOVariables().

◆ _fp

MaterialProperty<RankTwoTensor>& FiniteStrainHyperElasticViscoPlastic::_fp
protected

Definition at line 192 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by initQpStatefulProperties(), and postSolveQp().

◆ _fp_old

const MaterialProperty<RankTwoTensor>& FiniteStrainHyperElasticViscoPlastic::_fp_old
protected

Definition at line 193 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<const Function *> ComputeGeneralStressBase::_initial_stress_fcn
protectedinherited

initial stress components

Definition at line 56 of file ComputeGeneralStressBase.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 188 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 170 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 186 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 168 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>& ComputeGeneralStressBase::_Jacobian_mult
protectedinherited

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

Definition at line 59 of file ComputeGeneralStressBase.h.

Referenced by ComputeStrainIncrementBasedStress::computeQpJacobian(), computeQpJacobian(), ComputeMultipleInelasticStressBase::computeQpJacobianMult(), ComputeMultipleInelasticCosseratStress::computeQpJacobianMult(), ComputeMultiPlasticityStress::computeQpStress(), FiniteStrainPlasticMaterial::computeQpStress(), ComputeCosseratLinearElasticStress::computeQpStress(), ComputeFiniteStrainElasticStress::computeQpStress(), ComputeLinearElasticPFFractureStress::computeQpStress(), ComputeLinearElasticStress::computeQpStress(), ComputeLinearViscoelasticStress::computeQpStress(), ComputeMultipleInelasticStressBase::computeQpStress(), ComputeSmearedCrackingStress::computeQpStress(), ComputeMultipleCrystalPlasticityStress::computeQpStress(), ComputeMultipleInelasticStressBase::computeQpStressIntermediateConfiguration(), ComputeLinearElasticPFFractureStress::computeStrainSpectral(), ComputeLinearElasticPFFractureStress::computeStrainVolDev(), ComputeLinearElasticPFFractureStress::computeStressSpectral(), FiniteStrainUObasedCP::elasticTangentModuli(), FiniteStrainUObasedCP::elastoPlasticTangentModuli(), ComputeMultipleInelasticStressBase::finiteStrainRotation(), ComputeMultiPlasticityStress::postReturnMap(), FiniteStrainCrystalPlasticity::postSolveQp(), FiniteStrainCrystalPlasticity::preSolveQp(), and ComputeMultipleInelasticStressBase::updateQpStateSingleModel().

◆ _max_substep_iter

unsigned int FiniteStrainHyperElasticViscoPlastic::_max_substep_iter
protected

Maximum number of substep iterations.

Definition at line 161 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeQpStress().

◆ _maxiters

unsigned int FiniteStrainHyperElasticViscoPlastic::_maxiters
protected

Maximum number of iterations.

Definition at line 159 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by solveFlowrate().

◆ _mechanical_strain

const MaterialProperty<RankTwoTensor>& ComputeGeneralStressBase::_mechanical_strain
protectedinherited

Mechanical strain material property.

Definition at line 46 of file ComputeGeneralStressBase.h.

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

◆ _num_flow_rate_uos

unsigned int FiniteStrainHyperElasticViscoPlastic::_num_flow_rate_uos
protected

Number of flow rate user objects.

Definition at line 173 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 179 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 177 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 175 of file FiniteStrainHyperElasticViscoPlastic.h.

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

◆ _pk2

MaterialProperty<RankTwoTensor>& FiniteStrainHyperElasticViscoPlastic::_pk2
protected

Definition at line 191 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by HyperElasticPhaseFieldIsoDamage::computeNumStiffness(), computePK2StressAndDerivative(), HyperElasticPhaseFieldIsoDamage::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 190 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 155 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 157 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 184 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 166 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeFlowRateJacobian(), and initUOVariables().

◆ _stress

MaterialProperty<RankTwoTensor>& ComputeGeneralStressBase::_stress
protectedinherited

Stress material property.

Definition at line 48 of file ComputeGeneralStressBase.h.

Referenced by ComputeMultipleInelasticStressBase::computeAdmissibleState(), ComputeMultipleInelasticCosseratStress::computeAdmissibleState(), ComputeGeneralStressBase::computeQpProperties(), ComputeMultiPlasticityStress::computeQpStress(), ComputeStrainIncrementBasedStress::computeQpStress(), FiniteStrainPlasticMaterial::computeQpStress(), AbaqusUMATStress::computeQpStress(), ComputeCosseratLinearElasticStress::computeQpStress(), ComputeFiniteStrainElasticStress::computeQpStress(), ComputeLinearElasticPFFractureStress::computeQpStress(), ComputeLinearElasticStress::computeQpStress(), ComputeLinearViscoelasticStress::computeQpStress(), ComputeMultipleInelasticStressBase::computeQpStress(), ComputeSmearedCrackingStress::computeQpStress(), ComputeMultipleCrystalPlasticityStress::computeQpStress(), ComputeMultipleInelasticStressBase::computeQpStressIntermediateConfiguration(), ComputeLinearElasticPFFractureStress::computeStrainSpectral(), ComputeLinearElasticPFFractureStress::computeStrainVolDev(), ComputeLinearElasticPFFractureStress::computeStressSpectral(), ComputeCreepPlasticityStress::computeTangentOperators(), ComputeMultipleInelasticStressBase::finiteStrainRotation(), FiniteStrainCrystalPlasticity::initQpStatefulProperties(), FiniteStrainUObasedCP::initQpStatefulProperties(), initQpStatefulProperties(), ComputeGeneralStressBase::initQpStatefulProperties(), ComputeMultiPlasticityStress::postReturnMap(), FiniteStrainCrystalPlasticity::postSolveQp(), FiniteStrainUObasedCP::postSolveQp(), postSolveQp(), ComputeSmearedCrackingStress::updateCrackingStateAndStress(), ComputeCreepPlasticityStress::updateQpState(), ComputeMultipleInelasticStress::updateQpState(), and ComputeMultipleInelasticStressBase::updateQpStateSingleModel().

◆ _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:
Generated on Fri Aug 21 2026 13:42:44 for https://mooseframework.inl.gov by doxygen 1.9.8