This class solves visco plastic model based on isotropically damaged stress The damage parameter is obtained from phase field fracture kernel Computes undamaged elastic strain energy and associated tensors used in phase field fracture kernel. More...

#include <HyperElasticPhaseFieldIsoDamage.h>

Inheritance diagram for HyperElasticPhaseFieldIsoDamage:

[legend]

Public Types
typedef DerivativeMaterialPropertyNameInterface::SymbolName	SymbolName

Public Member Functions
	HyperElasticPhaseFieldIsoDamage (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	computePK2StressAndDerivative ()
	This function computes PK2 stress.

virtual void	computeDamageStress ()
	This function computes PK2 stress modified to account for damage Computes numerical stiffness if flag is true Computes undamaged elastic strain energy and associated tensors.

virtual void	computeNumStiffness ()
	This function computes numerical stiffness.

virtual void	computeQpJacobian ()
	This function computes tensors used to construct diagonal and off-diagonal Jacobian.

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	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	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
bool	_num_stiffness
	Flag to compute numerical stiffness.

Real	_kdamage
	Small stiffness of completely damaged material point.

bool	_use_current_hist
	Use current value of history variable.

const MaterialProperty< Real > &	_l
	Material property defining crack width, declared elsewhere.

const MaterialProperty< Real > &	_gc
	Material property defining gc parameter, declared elsewhere.

Real	_zero_tol
	Used in numerical stiffness calculation to check near zero values.

Real	_zero_pert
	Perturbation value for near zero or zero strain components.

Real	_pert_val
	Perturbation value for strain components.

const VariableValue &	_c
	Compupled damage variable.

bool	_save_state
	Flag to save couple material properties.

MaterialProperty< RankTwoTensor > &	_dstress_dc

RankTwoTensor	_pk2_tmp

RankTwoTensor	_dG0_dee

RankTwoTensor	_dpk2_dc

RankFourTensor	_dpk2_dee

std::vector< RankTwoTensor >	_etens

MaterialProperty< Real > &	_F
	Elastic energy and derivatives, declared in this material.

MaterialProperty< Real > &	_dFdc

MaterialProperty< Real > &	_d2Fdc2

MaterialProperty< RankTwoTensor > &	_d2Fdcdstrain

MaterialProperty< Real > &	_hist
	History variable that prevents crack healing, declared in this material.

const MaterialProperty< Real > &	_hist_old
	Old value of history variable.

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 visco plastic model based on isotropically damaged stress The damage parameter is obtained from phase field fracture kernel Computes undamaged elastic strain energy and associated tensors used in phase field fracture kernel.

Definition at line 20 of file HyperElasticPhaseFieldIsoDamage.h.

Constructor & Destructor Documentation

◆ HyperElasticPhaseFieldIsoDamage()

HyperElasticPhaseFieldIsoDamage::HyperElasticPhaseFieldIsoDamage ( const InputParameters & parameters )

Definition at line 36 of file HyperElasticPhaseFieldIsoDamage.C.

  : FiniteStrainHyperElasticViscoPlastic(parameters),
    _num_stiffness(getParam<bool>("numerical_stiffness")),
    _kdamage(getParam<Real>("damage_stiffness")),
    _use_current_hist(getParam<bool>("use_current_history_variable")),
    _l(getMaterialProperty<Real>("l")),
    _gc(getMaterialProperty<Real>("gc_prop")),
    _zero_tol(getParam<Real>("zero_tol")),
    _zero_pert(getParam<Real>("zero_perturb")),
    _pert_val(getParam<Real>("perturbation_scale_factor")),
    _c(coupledValue("c")),
    _save_state(false),
    _dstress_dc(
        declarePropertyDerivative<RankTwoTensor>(_base_name + "stress", coupledName("c", 0))),
    _etens(LIBMESH_DIM),
    _F(declareProperty<Real>(getParam<MaterialPropertyName>("F_name"))),
    _dFdc(declarePropertyDerivative<Real>(getParam<MaterialPropertyName>("F_name"),
                                          coupledName("c", 0))),
    _d2Fdc2(declarePropertyDerivative<Real>(
        getParam<MaterialPropertyName>("F_name"), coupledName("c", 0), coupledName("c", 0))),
    _d2Fdcdstrain(declareProperty<RankTwoTensor>("d2Fdcdstrain")),
    _hist(declareProperty<Real>("hist")),
    _hist_old(getMaterialPropertyOld<Real>("hist"))
{
}

Member Function Documentation

◆ computeDamageStress()

void HyperElasticPhaseFieldIsoDamage::computeDamageStress ( )

protectedvirtual

This function computes PK2 stress modified to account for damage Computes numerical stiffness if flag is true Computes undamaged elastic strain energy and associated tensors.

Definition at line 91 of file HyperElasticPhaseFieldIsoDamage.C.

{
  Real lambda = _elasticity_tensor[_qp](0, 0, 1, 1);
  Real mu = _elasticity_tensor[_qp](0, 1, 0, 1);
 
  Real c = _c[_qp];
  Real xfac = Utility::pow<2>(1.0 - c) + _kdamage;
 
  std::vector<Real> eigval;
  RankTwoTensor evec;
  _ee.symmetricEigenvaluesEigenvectors(eigval, evec);
 
  for (const auto i : make_range(Moose::dim))
    _etens[i] = RankTwoTensor::selfOuterProduct(evec.column(i));
 
  Real etr = 0.0;
  for (const auto i : make_range(Moose::dim))
    etr += eigval[i];
 
  Real etrpos = (std::abs(etr) + etr) / 2.0;
  Real etrneg = (std::abs(etr) - etr) / 2.0;
 
  RankTwoTensor pk2pos, pk2neg;
 
  for (const auto i : make_range(Moose::dim))
  {
    pk2pos += _etens[i] * (lambda * etrpos + 2.0 * mu * (std::abs(eigval[i]) + eigval[i]) / 2.0);
    pk2neg += _etens[i] * (lambda * etrneg + 2.0 * mu * (std::abs(eigval[i]) - eigval[i]) / 2.0);
  }
 
  _pk2_tmp = pk2pos * xfac - pk2neg;
 
  if (_save_state)
  {
    std::vector<Real> epos(LIBMESH_DIM), eneg(LIBMESH_DIM);
    for (const auto i : make_range(Moose::dim))
    {
      epos[i] = (std::abs(eigval[i]) + eigval[i]) / 2.0;
      eneg[i] = (std::abs(eigval[i]) - eigval[i]) / 2.0;
    }
 
    // sum squares of epos and eneg
    Real pval(0.0), nval(0.0);
    for (const auto i : make_range(Moose::dim))
    {
      pval += epos[i] * epos[i];
      nval += eneg[i] * eneg[i];
    }
 
    // Energy with positive principal strains
    const Real G0_pos = lambda * etrpos * etrpos / 2.0 + mu * pval;
    const Real G0_neg = lambda * etrneg * etrneg / 2.0 + mu * nval;
 
    // Assign history variable and derivative
    if (G0_pos > _hist_old[_qp])
      _hist[_qp] = G0_pos;
    else
      _hist[_qp] = _hist_old[_qp];
 
    Real hist_variable = _hist_old[_qp];
    if (_use_current_hist)
      hist_variable = _hist[_qp];
 
    // Elastic free energy density
    _F[_qp] = hist_variable * xfac - G0_neg + _gc[_qp] / (2 * _l[_qp]) * c * c;
 
    // derivative of elastic free energy density wrt c
    _dFdc[_qp] = -hist_variable * 2.0 * (1.0 - c) * (1 - _kdamage) + _gc[_qp] / _l[_qp] * c;
 
    // 2nd derivative of elastic free energy density wrt c
    _d2Fdc2[_qp] = hist_variable * 2.0 * (1 - _kdamage) + _gc[_qp] / _l[_qp];
 
    _dG0_dee = pk2pos;
 
    _dpk2_dc = -pk2pos * 2.0 * (1.0 - c);
  }
}

Referenced by computeNumStiffness(), and computePK2StressAndDerivative().

◆ computeDeeDce()

void FiniteStrainHyperElasticViscoPlastic::computeDeeDce ( )

protectedvirtualinherited

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

Definition at line 487 of file FiniteStrainHyperElasticViscoPlastic.C.

{
  _dee_dce.zero();
 
  for (const auto i : make_range(Moose::dim))
    for (const auto j : make_range(Moose::dim))
      _dee_dce(i, j, i, j) = 0.5;
}

Referenced by FiniteStrainHyperElasticViscoPlastic::initJacobianVariables().

◆ computeDpk2Dfpinv()

void FiniteStrainHyperElasticViscoPlastic::computeDpk2Dfpinv ( )

protectedvirtualinherited

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

Definition at line 517 of file FiniteStrainHyperElasticViscoPlastic.C.

{
  for (const auto i : make_range(Moose::dim))
    for (const auto j : make_range(Moose::dim))
      for (const auto k : make_range(Moose::dim))
        _dfe_dfpinv(i, j, k, j) = _dfgrd_tmp(i, k);
 
  _dpk2_dfpinv = _dpk2_dce * _dce_dfe * _dfe_dfpinv;
}

Referenced by FiniteStrainHyperElasticViscoPlastic::computeFlowRateJacobian().

◆ computeElasticPlasticDeformGrad()

void FiniteStrainHyperElasticViscoPlastic::computeElasticPlasticDeformGrad ( )

protectedvirtualinherited

Computes elastic and plastic deformation gradients.

Definition at line 503 of file FiniteStrainHyperElasticViscoPlastic.C.

{
  RankTwoTensor iden(RankTwoTensor::initIdentity);
 
  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;
}

Referenced by FiniteStrainHyperElasticViscoPlastic::solveFlowrate().

◆ computeElasticRightCauchyGreenTensor()

void FiniteStrainHyperElasticViscoPlastic::computeElasticRightCauchyGreenTensor ( )

protectedvirtualinherited

Computes elastic Right Cauchy Green Tensor.

Definition at line 497 of file FiniteStrainHyperElasticViscoPlastic.C.

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

Referenced by FiniteStrainHyperElasticViscoPlastic::computeFlowRateResidual().

◆ computeElasticStrain()

void FiniteStrainHyperElasticViscoPlastic::computeElasticStrain ( )

protectedvirtualinherited

Computes elastic Lagrangian strain.

Definition at line 480 of file FiniteStrainHyperElasticViscoPlastic.C.

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

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

◆ computeFlowDirection()

bool FiniteStrainHyperElasticViscoPlastic::computeFlowDirection ( )

protectedvirtualinherited

Calls user objects to compute flow directions.

Definition at line 437 of file FiniteStrainHyperElasticViscoPlastic.C.

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

Referenced by FiniteStrainHyperElasticViscoPlastic::computeFlowRateResidual().

◆ computeFlowRateFunction()

bool FiniteStrainHyperElasticViscoPlastic::computeFlowRateFunction ( )

protectedvirtualinherited

Calls user objects to compute flow rates.

Definition at line 448 of file FiniteStrainHyperElasticViscoPlastic.C.

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

Referenced by FiniteStrainHyperElasticViscoPlastic::computeFlowRateResidual().

◆ computeFlowRateJacobian()

void FiniteStrainHyperElasticViscoPlastic::computeFlowRateJacobian ( )

protectedvirtualinherited

Computes flow rate Jacobian matrix.

Definition at line 399 of file FiniteStrainHyperElasticViscoPlastic.C.

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

Referenced by FiniteStrainHyperElasticViscoPlastic::solveFlowrate().

◆ computeFlowRateResidual()

bool FiniteStrainHyperElasticViscoPlastic::computeFlowRateResidual ( )

protectedvirtualinherited

Computes flow rate residual vector.

Definition at line 373 of file FiniteStrainHyperElasticViscoPlastic.C.

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

Referenced by FiniteStrainHyperElasticViscoPlastic::solveFlowrate().

◆ computeIntVar()

bool FiniteStrainHyperElasticViscoPlastic::computeIntVar ( )

protectedvirtualinherited

This function call user objects to integrate internal variables.

Definition at line 594 of file FiniteStrainHyperElasticViscoPlastic.C.

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

Referenced by FiniteStrainHyperElasticViscoPlastic::computeFlowRateResidual().

◆ computeIntVarDerivatives()

void FiniteStrainHyperElasticViscoPlastic::computeIntVarDerivatives ( )

protectedvirtualinherited

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

Definition at line 635 of file FiniteStrainHyperElasticViscoPlastic.C.

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

Referenced by FiniteStrainHyperElasticViscoPlastic::computeFlowRateJacobian().

◆ computeIntVarRateDerivatives()

void FiniteStrainHyperElasticViscoPlastic::computeIntVarRateDerivatives ( )

protectedvirtualinherited

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

Definition at line 622 of file FiniteStrainHyperElasticViscoPlastic.C.

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

Referenced by FiniteStrainHyperElasticViscoPlastic::computeFlowRateJacobian().

◆ computeIntVarRates()

bool FiniteStrainHyperElasticViscoPlastic::computeIntVarRates ( )

protectedvirtualinherited

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

Definition at line 580 of file FiniteStrainHyperElasticViscoPlastic.C.

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

Referenced by FiniteStrainHyperElasticViscoPlastic::computeFlowRateResidual().

◆ computeNorm()

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

protectedvirtualinherited

Computes norm of residual vector.

Definition at line 528 of file FiniteStrainHyperElasticViscoPlastic.C.

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

Referenced by FiniteStrainHyperElasticViscoPlastic::solveFlowrate().

◆ computeNumStiffness()

void HyperElasticPhaseFieldIsoDamage::computeNumStiffness ( )

protectedvirtual

This function computes numerical stiffness.

Definition at line 170 of file HyperElasticPhaseFieldIsoDamage.C.

{
  RankTwoTensor ee_tmp;
 
  for (const auto i : make_range(Moose::dim))
    for (unsigned int j = i; j < LIBMESH_DIM; ++j)
    {
      Real ee_pert = _zero_pert;
      if (std::abs(_ee(i, j)) > _zero_tol)
        ee_pert = _pert_val * std::abs(_ee(i, j));
 
      ee_tmp = _ee;
      _ee(i, j) += ee_pert;
 
      computeDamageStress();
 
      for (const auto k : make_range(Moose::dim))
        for (const auto l : make_range(Moose::dim))
        {
          _dpk2_dee(k, l, i, j) = (_pk2_tmp(k, l) - _pk2[_qp](k, l)) / ee_pert;
          _dpk2_dee(k, l, j, i) = (_pk2_tmp(k, l) - _pk2[_qp](k, l)) / ee_pert;
        }
      _ee = ee_tmp;
    }
}

Referenced by computePK2StressAndDerivative().

◆ computePK2StressAndDerivative()

void HyperElasticPhaseFieldIsoDamage::computePK2StressAndDerivative ( )

protectedvirtual

This function computes PK2 stress.

Reimplemented from FiniteStrainHyperElasticViscoPlastic.

Definition at line 63 of file HyperElasticPhaseFieldIsoDamage.C.

{
  computeElasticStrain();
 
  _save_state = true;
  computeDamageStress();
  _pk2[_qp] = _pk2_tmp;
 
  _save_state = false;
  if (_num_stiffness)
    computeNumStiffness();
 
  if (_num_stiffness)
    _dpk2_dce = _dpk2_dee * _dee_dce;
 
  _dce_dfe.zero();
  for (const auto i : make_range(Moose::dim))
    for (const auto j : make_range(Moose::dim))
      for (const auto k : make_range(Moose::dim))
      {
        _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 HyperElasticPhaseFieldIsoDamage::computeQpJacobian ( )

protectedvirtual

This function computes tensors used to construct diagonal and off-diagonal Jacobian.

Reimplemented from FiniteStrainHyperElasticViscoPlastic.

Definition at line 197 of file HyperElasticPhaseFieldIsoDamage.C.

{
  FiniteStrainHyperElasticViscoPlastic::computeQpJacobian();
 
  RankTwoTensor dG0_dce = _dee_dce.innerProductTranspose(_dG0_dee);
  RankTwoTensor dG0_dfe = _dce_dfe.innerProductTranspose(dG0_dce);
  RankTwoTensor dG0_df = _dfe_df.innerProductTranspose(dG0_dfe);
 
  // 2nd derivative wrt c and strain = 0.0 if we used the previous step's history varible
  if (_use_current_hist)
    _d2Fdcdstrain[_qp] =
        -_df_dstretch_inc.innerProductTranspose(dG0_df) * 2.0 * (1.0 - _c[_qp]) * (1 - _kdamage);
 
  usingTensorIndices(i_, j_, k_, l_);
  _dstress_dc[_qp] = _fe.times<i_, k_, j_, l_>(_fe) * _dpk2_dc;
}

◆ computeQpProperties()

void ComputeGeneralStressBase::computeQpProperties ( )

overrideprotectedvirtualinherited

Definition at line 44 of file ComputeGeneralStressBase.C.

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

◆ computeQpStress()

void FiniteStrainHyperElasticViscoPlastic::computeQpStress ( )

protectedvirtualinherited

This function computes the Cauchy stress.

Implements ComputeGeneralStressBase.

Definition at line 208 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 ( )

protectedvirtualinherited

This function call user objects to compute strength.

Definition at line 608 of file FiniteStrainHyperElasticViscoPlastic.C.

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

Referenced by FiniteStrainHyperElasticViscoPlastic::computeFlowRateResidual().

◆ computeStrengthDerivatives()

void FiniteStrainHyperElasticViscoPlastic::computeStrengthDerivatives ( )

protectedinherited

This function call user objects to compute dstrength/dintvar.

Definition at line 670 of file FiniteStrainHyperElasticViscoPlastic.C.

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

Referenced by FiniteStrainHyperElasticViscoPlastic::computeFlowRateJacobian().

◆ initJacobianVariables()

void FiniteStrainHyperElasticViscoPlastic::initJacobianVariables ( )

protectedvirtualinherited

This function initialize variables required for Jacobian calculation.

Definition at line 155 of file FiniteStrainHyperElasticViscoPlastic.C.

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

Referenced by FiniteStrainHyperElasticViscoPlastic::FiniteStrainHyperElasticViscoPlastic().

◆ initNumUserObjects()

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

protectedinherited

This function calculates the number of each user object type.

Definition at line 110 of file FiniteStrainHyperElasticViscoPlastic.C.

{
  uo_num = uo_names.size();
}

Referenced by FiniteStrainHyperElasticViscoPlastic::initUOVariables().

◆ initProp()

template<typename T >

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

protectedinherited

This function initializes properties for each user object.

Definition at line 118 of file FiniteStrainHyperElasticViscoPlastic.C.

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

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

protectedinherited

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.

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

Referenced by FiniteStrainHyperElasticViscoPlastic::initUOVariables().

◆ initQpStatefulProperties()

void FiniteStrainHyperElasticViscoPlastic::initQpStatefulProperties ( )

protectedvirtualinherited

Initializes state.

Reimplemented from ComputeGeneralStressBase.

Definition at line 183 of file FiniteStrainHyperElasticViscoPlastic.C.

{
  _stress[_qp].zero();
  _ce[_qp].zero();
  _pk2[_qp].zero();
  _fp[_qp].setToIdentity();
 
  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 ( )

protectedvirtualinherited

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.

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

Referenced by FiniteStrainHyperElasticViscoPlastic::FiniteStrainHyperElasticViscoPlastic().

◆ initUserObjects()

template<typename T >

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

protectedinherited

This function initializes user objects.

Definition at line 141 of file FiniteStrainHyperElasticViscoPlastic.C.

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

Referenced by FiniteStrainHyperElasticViscoPlastic::initUOVariables().

◆ postSolveFlowrate()

void FiniteStrainHyperElasticViscoPlastic::postSolveFlowrate ( )

protectedvirtualinherited

Update state for output (Inside substepping)

Definition at line 362 of file FiniteStrainHyperElasticViscoPlastic.C.

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

Referenced by FiniteStrainHyperElasticViscoPlastic::solveQp().

◆ postSolveQp()

void FiniteStrainHyperElasticViscoPlastic::postSolveQp ( )

protectedvirtualinherited

Update state for output (Outside substepping)

Definition at line 280 of file FiniteStrainHyperElasticViscoPlastic.C.

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

Referenced by FiniteStrainHyperElasticViscoPlastic::computeQpStress().

◆ preSolveFlowrate()

void FiniteStrainHyperElasticViscoPlastic::preSolveFlowrate ( )

protectedvirtualinherited

Sets state for solve (Inside substepping)

Definition at line 299 of file FiniteStrainHyperElasticViscoPlastic.C.

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

Referenced by FiniteStrainHyperElasticViscoPlastic::solveQp().

◆ preSolveQp()

void FiniteStrainHyperElasticViscoPlastic::preSolveQp ( )

protectedvirtualinherited

Sets state for solve.

Definition at line 256 of file FiniteStrainHyperElasticViscoPlastic.C.

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

Referenced by FiniteStrainHyperElasticViscoPlastic::computeQpStress().

◆ recoverOldState()

void FiniteStrainHyperElasticViscoPlastic::recoverOldState ( )

protectedvirtualinherited

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

Definition at line 291 of file FiniteStrainHyperElasticViscoPlastic.C.

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

Referenced by FiniteStrainHyperElasticViscoPlastic::postSolveQp().

◆ saveOldState()

void FiniteStrainHyperElasticViscoPlastic::saveOldState ( )

protectedvirtualinherited

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

Definition at line 249 of file FiniteStrainHyperElasticViscoPlastic.C.

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

Referenced by FiniteStrainHyperElasticViscoPlastic::computeQpStress().

◆ solveFlowrate()

bool FiniteStrainHyperElasticViscoPlastic::solveFlowrate ( )

protectedvirtualinherited

Solve for flow rate and state.

Definition at line 315 of file FiniteStrainHyperElasticViscoPlastic.C.

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

Referenced by FiniteStrainHyperElasticViscoPlastic::solveQp().

◆ solveQp()

bool FiniteStrainHyperElasticViscoPlastic::solveQp ( )

protectedvirtualinherited

Solve state.

Definition at line 269 of file FiniteStrainHyperElasticViscoPlastic.C.

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

Referenced by FiniteStrainHyperElasticViscoPlastic::computeQpStress().

◆ updateFlowRate()

void FiniteStrainHyperElasticViscoPlastic::updateFlowRate ( )

protectedvirtualinherited

Update flow rate.

Definition at line 537 of file FiniteStrainHyperElasticViscoPlastic.C.

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

Referenced by FiniteStrainHyperElasticViscoPlastic::solveFlowrate().

◆ validParams()

InputParameters HyperElasticPhaseFieldIsoDamage::validParams ( )

static

Definition at line 16 of file HyperElasticPhaseFieldIsoDamage.C.

{
  InputParameters params = FiniteStrainHyperElasticViscoPlastic::validParams();
  params.addParam<bool>("numerical_stiffness", false, "Flag for numerical stiffness");
  params.addParam<Real>("damage_stiffness", 1e-8, "Avoid zero after complete damage");
  params.addParam<Real>("zero_tol", 1e-12, "Tolerance for numerical zero");
  params.addParam<Real>(
      "zero_perturb", 1e-8, "Perturbation value when strain value less than numerical zero");
  params.addParam<Real>("perturbation_scale_factor", 1e-5, "Perturbation scale factor");
  params.addRequiredCoupledVar("c", "Damage variable");
  params.addParam<bool>(
      "use_current_history_variable", false, "Use the current value of the history variable.");
  params.addParam<MaterialPropertyName>(
      "F_name", "E_el", "Name of material property storing the elastic energy");
  params.addClassDescription(
      "Computes damaged stress and energy in the intermediate configuration assuming isotropy");
 
  return params;
}

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

◆ _c

const VariableValue& HyperElasticPhaseFieldIsoDamage::_c

protected

Compupled damage variable.

Definition at line 58 of file HyperElasticPhaseFieldIsoDamage.h.

Referenced by computeDamageStress(), and computeQpJacobian().

◆ _ce

MaterialProperty<RankTwoTensor>& FiniteStrainHyperElasticViscoPlastic::_ce

protectedinherited

Definition at line 194 of file FiniteStrainHyperElasticViscoPlastic.h.

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

◆ _d2Fdc2

MaterialProperty<Real>& HyperElasticPhaseFieldIsoDamage::_d2Fdc2

protected

Definition at line 71 of file HyperElasticPhaseFieldIsoDamage.h.

Referenced by computeDamageStress().

◆ _d2Fdcdstrain

MaterialProperty<RankTwoTensor>& HyperElasticPhaseFieldIsoDamage::_d2Fdcdstrain

protected

Definition at line 72 of file HyperElasticPhaseFieldIsoDamage.h.

Referenced by computeQpJacobian().

◆ _dce_dfe

RankFourTensor FiniteStrainHyperElasticViscoPlastic::_dce_dfe

protectedinherited

Definition at line 216 of file FiniteStrainHyperElasticViscoPlastic.h.

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

◆ _dee_dce

RankFourTensor FiniteStrainHyperElasticViscoPlastic::_dee_dce

protectedinherited

Definition at line 216 of file FiniteStrainHyperElasticViscoPlastic.h.

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

◆ _deformation_gradient

const MaterialProperty<RankTwoTensor>& FiniteStrainHyperElasticViscoPlastic::_deformation_gradient

protectedinherited

Definition at line 200 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by FiniteStrainHyperElasticViscoPlastic::computeQpStress().

◆ _deformation_gradient_old

const MaterialProperty<RankTwoTensor>& FiniteStrainHyperElasticViscoPlastic::_deformation_gradient_old

protectedinherited

Definition at line 201 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by FiniteStrainHyperElasticViscoPlastic::computeQpJacobian(), and FiniteStrainHyperElasticViscoPlastic::computeQpStress().

◆ _df_dstretch_inc

RankFourTensor FiniteStrainHyperElasticViscoPlastic::_df_dstretch_inc

protectedinherited

Definition at line 217 of file FiniteStrainHyperElasticViscoPlastic.h.

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

◆ _dFdc

MaterialProperty<Real>& HyperElasticPhaseFieldIsoDamage::_dFdc

protected

Definition at line 70 of file HyperElasticPhaseFieldIsoDamage.h.

Referenced by computeDamageStress().

◆ _dfe_df

RankFourTensor FiniteStrainHyperElasticViscoPlastic::_dfe_df

protectedinherited

Definition at line 216 of file FiniteStrainHyperElasticViscoPlastic.h.

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

◆ _dfe_dfpinv

RankFourTensor FiniteStrainHyperElasticViscoPlastic::_dfe_dfpinv

protectedinherited

Definition at line 215 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by FiniteStrainHyperElasticViscoPlastic::computeDpk2Dfpinv().

◆ _dfgrd_tmp

RankTwoTensor FiniteStrainHyperElasticViscoPlastic::_dfgrd_tmp

protectedinherited

Definition at line 211 of file FiniteStrainHyperElasticViscoPlastic.h.

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

◆ _dflow_rate

DenseVector<Real> FiniteStrainHyperElasticViscoPlastic::_dflow_rate

protectedinherited

Definition at line 224 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by FiniteStrainHyperElasticViscoPlastic::FiniteStrainHyperElasticViscoPlastic(), and FiniteStrainHyperElasticViscoPlastic::updateFlowRate().

◆ _dflowrate_dpk2

std::vector<RankTwoTensor> FiniteStrainHyperElasticViscoPlastic::_dflowrate_dpk2

protectedinherited

Definition at line 220 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by FiniteStrainHyperElasticViscoPlastic::computeFlowRateJacobian(), and FiniteStrainHyperElasticViscoPlastic::initJacobianVariables().

◆ _dflowrate_dstrength

DenseMatrix<Real> FiniteStrainHyperElasticViscoPlastic::_dflowrate_dstrength

protectedinherited

Definition at line 237 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by FiniteStrainHyperElasticViscoPlastic::computeFlowRateJacobian(), and FiniteStrainHyperElasticViscoPlastic::initJacobianVariables().

◆ _dfpinv_dflowrate

std::vector<RankTwoTensor> FiniteStrainHyperElasticViscoPlastic::_dfpinv_dflowrate

protectedinherited

Definition at line 222 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by FiniteStrainHyperElasticViscoPlastic::computeFlowRateJacobian(), and FiniteStrainHyperElasticViscoPlastic::initJacobianVariables().

◆ _dG0_dee

RankTwoTensor HyperElasticPhaseFieldIsoDamage::_dG0_dee

protected

Definition at line 64 of file HyperElasticPhaseFieldIsoDamage.h.

Referenced by computeDamageStress(), and computeQpJacobian().

◆ _dintvar_dflowrate

DenseMatrix<Real> FiniteStrainHyperElasticViscoPlastic::_dintvar_dflowrate

protectedinherited

Definition at line 235 of file FiniteStrainHyperElasticViscoPlastic.h.

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

◆ _dintvar_dflowrate_tmp

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

protectedinherited

Definition at line 230 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by FiniteStrainHyperElasticViscoPlastic::computeIntVarDerivatives(), and FiniteStrainHyperElasticViscoPlastic::initJacobianVariables().

◆ _dintvar_dintvar

DenseMatrix<Real> FiniteStrainHyperElasticViscoPlastic::_dintvar_dintvar

protectedinherited

Definition at line 234 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by FiniteStrainHyperElasticViscoPlastic::computeIntVarDerivatives(), and FiniteStrainHyperElasticViscoPlastic::initJacobianVariables().

◆ _dintvar_dintvar_x

DenseVector<Real> FiniteStrainHyperElasticViscoPlastic::_dintvar_dintvar_x

protectedinherited

Definition at line 238 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by FiniteStrainHyperElasticViscoPlastic::computeIntVarDerivatives(), and FiniteStrainHyperElasticViscoPlastic::initJacobianVariables().

◆ _dintvar_dintvarrate

DenseMatrix<Real> FiniteStrainHyperElasticViscoPlastic::_dintvar_dintvarrate

protectedinherited

Definition at line 233 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by FiniteStrainHyperElasticViscoPlastic::computeIntVarDerivatives(), and FiniteStrainHyperElasticViscoPlastic::initJacobianVariables().

◆ _dintvarrate_dflowrate

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

protectedinherited

Jacobian variables.

Definition at line 229 of file FiniteStrainHyperElasticViscoPlastic.h.

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

◆ _dintvarrate_dintvar

DenseMatrix<Real> FiniteStrainHyperElasticViscoPlastic::_dintvarrate_dintvar

protectedinherited

Definition at line 232 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by FiniteStrainHyperElasticViscoPlastic::computeIntVarDerivatives(), and FiniteStrainHyperElasticViscoPlastic::initJacobianVariables().

◆ _dpk2_dc

RankTwoTensor HyperElasticPhaseFieldIsoDamage::_dpk2_dc

protected

Definition at line 64 of file HyperElasticPhaseFieldIsoDamage.h.

Referenced by computeDamageStress(), and computeQpJacobian().

◆ _dpk2_dce

RankFourTensor FiniteStrainHyperElasticViscoPlastic::_dpk2_dce

protectedinherited

Definition at line 215 of file FiniteStrainHyperElasticViscoPlastic.h.

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

◆ _dpk2_dee

RankFourTensor HyperElasticPhaseFieldIsoDamage::_dpk2_dee

protected

Definition at line 65 of file HyperElasticPhaseFieldIsoDamage.h.

Referenced by computeNumStiffness(), and computePK2StressAndDerivative().

◆ _dpk2_dfe

RankFourTensor FiniteStrainHyperElasticViscoPlastic::_dpk2_dfe

protectedinherited

Definition at line 215 of file FiniteStrainHyperElasticViscoPlastic.h.

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

◆ _dpk2_dflowrate

std::vector<RankTwoTensor> FiniteStrainHyperElasticViscoPlastic::_dpk2_dflowrate

protectedinherited

Definition at line 221 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by FiniteStrainHyperElasticViscoPlastic::computeFlowRateJacobian(), and FiniteStrainHyperElasticViscoPlastic::initJacobianVariables().

◆ _dpk2_dfpinv

RankFourTensor FiniteStrainHyperElasticViscoPlastic::_dpk2_dfpinv

protectedinherited

Definition at line 215 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by FiniteStrainHyperElasticViscoPlastic::computeDpk2Dfpinv(), and FiniteStrainHyperElasticViscoPlastic::computeFlowRateJacobian().

◆ _dstrength_dintvar

DenseMatrix<Real> FiniteStrainHyperElasticViscoPlastic::_dstrength_dintvar

protectedinherited

Definition at line 236 of file FiniteStrainHyperElasticViscoPlastic.h.

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

◆ _dstress_dc

MaterialProperty<RankTwoTensor>& HyperElasticPhaseFieldIsoDamage::_dstress_dc

protected

Definition at line 62 of file HyperElasticPhaseFieldIsoDamage.h.

Referenced by computeQpJacobian().

◆ _dt_substep

Real FiniteStrainHyperElasticViscoPlastic::_dt_substep

protectedinherited

Definition at line 241 of file FiniteStrainHyperElasticViscoPlastic.h.

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

◆ _ee

RankTwoTensor FiniteStrainHyperElasticViscoPlastic::_ee

protectedinherited

Definition at line 213 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by computeDamageStress(), FiniteStrainHyperElasticViscoPlastic::computeElasticStrain(), computeNumStiffness(), and FiniteStrainHyperElasticViscoPlastic::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

protectedinherited

Elasticity tensor material property.

Definition at line 199 of file FiniteStrainHyperElasticViscoPlastic.h.

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

◆ _elasticity_tensor_name

const std::string FiniteStrainHyperElasticViscoPlastic::_elasticity_tensor_name

protectedinherited

Name of the elasticity tensor material property.

Definition at line 197 of file FiniteStrainHyperElasticViscoPlastic.h.

◆ _etens

std::vector<RankTwoTensor> HyperElasticPhaseFieldIsoDamage::_etens

protected

Definition at line 66 of file HyperElasticPhaseFieldIsoDamage.h.

Referenced by computeDamageStress().

◆ _extra_stress

const MaterialProperty<RankTwoTensor>& ComputeGeneralStressBase::_extra_stress

protectedinherited

Extra stress tensor.

Definition at line 53 of file ComputeGeneralStressBase.h.

Referenced by ComputeGeneralStressBase::computeQpProperties().

◆ _F

MaterialProperty<Real>& HyperElasticPhaseFieldIsoDamage::_F

protected

Elastic energy and derivatives, declared in this material.

Definition at line 69 of file HyperElasticPhaseFieldIsoDamage.h.

Referenced by computeDamageStress().

◆ _fe

RankTwoTensor FiniteStrainHyperElasticViscoPlastic::_fe

protectedinherited

Definition at line 213 of file FiniteStrainHyperElasticViscoPlastic.h.

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

◆ _fe_pk2

RankTwoTensor FiniteStrainHyperElasticViscoPlastic::_fe_pk2

protectedinherited

Definition at line 214 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by FiniteStrainHyperElasticViscoPlastic::computeQpJacobian().

◆ _flow_dirn

std::vector<RankTwoTensor> FiniteStrainHyperElasticViscoPlastic::_flow_dirn

protectedinherited

Definition at line 219 of file FiniteStrainHyperElasticViscoPlastic.h.

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

◆ _flow_rate

DenseVector<Real> FiniteStrainHyperElasticViscoPlastic::_flow_rate

protectedinherited

Definition at line 225 of file FiniteStrainHyperElasticViscoPlastic.h.

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

◆ _flow_rate_prop

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

protectedinherited

Definition at line 204 of file FiniteStrainHyperElasticViscoPlastic.h.

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

◆ _flow_rate_uo

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

protectedinherited

Flow rate user objects.

Definition at line 182 of file FiniteStrainHyperElasticViscoPlastic.h.

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

◆ _flow_rate_uo_names

std::vector<UserObjectName> FiniteStrainHyperElasticViscoPlastic::_flow_rate_uo_names

protectedinherited

Names of flow rate user objects.

Definition at line 164 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by FiniteStrainHyperElasticViscoPlastic::computeIntVarRateDerivatives(), and FiniteStrainHyperElasticViscoPlastic::initUOVariables().

◆ _fp

MaterialProperty<RankTwoTensor>& FiniteStrainHyperElasticViscoPlastic::_fp

protectedinherited

Definition at line 192 of file FiniteStrainHyperElasticViscoPlastic.h.

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

◆ _fp_old

const MaterialProperty<RankTwoTensor>& FiniteStrainHyperElasticViscoPlastic::_fp_old

protectedinherited

Definition at line 193 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by FiniteStrainHyperElasticViscoPlastic::preSolveQp().

◆ _fp_tmp_inv

RankTwoTensor FiniteStrainHyperElasticViscoPlastic::_fp_tmp_inv

protectedinherited

Definition at line 212 of file FiniteStrainHyperElasticViscoPlastic.h.

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

◆ _fp_tmp_old_inv

RankTwoTensor FiniteStrainHyperElasticViscoPlastic::_fp_tmp_old_inv

protectedinherited

Definition at line 212 of file FiniteStrainHyperElasticViscoPlastic.h.

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

◆ _gc

const MaterialProperty<Real>& HyperElasticPhaseFieldIsoDamage::_gc

protected

Material property defining gc parameter, declared elsewhere.

Definition at line 50 of file HyperElasticPhaseFieldIsoDamage.h.

Referenced by computeDamageStress().

◆ _hist

MaterialProperty<Real>& HyperElasticPhaseFieldIsoDamage::_hist

protected

History variable that prevents crack healing, declared in this material.

Definition at line 75 of file HyperElasticPhaseFieldIsoDamage.h.

Referenced by computeDamageStress().

◆ _hist_old

const MaterialProperty<Real>& HyperElasticPhaseFieldIsoDamage::_hist_old

protected

Old value of history variable.

Definition at line 78 of file HyperElasticPhaseFieldIsoDamage.h.

Referenced by computeDamageStress().

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

protectedinherited

Definition at line 209 of file FiniteStrainHyperElasticViscoPlastic.h.

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

◆ _int_var_rate_prop

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

protectedinherited

Definition at line 208 of file FiniteStrainHyperElasticViscoPlastic.h.

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

◆ _int_var_rate_uo

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

protectedinherited

Internal variable rate user objects.

Definition at line 188 of file FiniteStrainHyperElasticViscoPlastic.h.

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

◆ _int_var_rate_uo_names

std::vector<UserObjectName> FiniteStrainHyperElasticViscoPlastic::_int_var_rate_uo_names

protectedinherited

Names of internal variable rate user objects.

Definition at line 170 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by FiniteStrainHyperElasticViscoPlastic::computeIntVarDerivatives(), and FiniteStrainHyperElasticViscoPlastic::initUOVariables().

◆ _int_var_stateful_prop

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

protectedinherited

Definition at line 206 of file FiniteStrainHyperElasticViscoPlastic.h.

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

◆ _int_var_stateful_prop_old

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

protectedinherited

Definition at line 207 of file FiniteStrainHyperElasticViscoPlastic.h.

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

◆ _int_var_uo

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

protectedinherited

Internal variable user objects.

Definition at line 186 of file FiniteStrainHyperElasticViscoPlastic.h.

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

◆ _int_var_uo_names

std::vector<UserObjectName> FiniteStrainHyperElasticViscoPlastic::_int_var_uo_names

protectedinherited

Names of internal variable user objects.

Definition at line 168 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by FiniteStrainHyperElasticViscoPlastic::computeStrengthDerivatives(), and FiniteStrainHyperElasticViscoPlastic::initUOVariables().

◆ _jac

DenseMatrix<Real> FiniteStrainHyperElasticViscoPlastic::_jac

protectedinherited

Definition at line 239 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by FiniteStrainHyperElasticViscoPlastic::computeFlowRateJacobian(), FiniteStrainHyperElasticViscoPlastic::initJacobianVariables(), and FiniteStrainHyperElasticViscoPlastic::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.

◆ _kdamage

Real HyperElasticPhaseFieldIsoDamage::_kdamage

protected

Small stiffness of completely damaged material point.

Definition at line 44 of file HyperElasticPhaseFieldIsoDamage.h.

Referenced by computeDamageStress(), and computeQpJacobian().

◆ _l

const MaterialProperty<Real>& HyperElasticPhaseFieldIsoDamage::_l

protected

Material property defining crack width, declared elsewhere.

Definition at line 48 of file HyperElasticPhaseFieldIsoDamage.h.

Referenced by computeDamageStress().

◆ _max_substep_iter

unsigned int FiniteStrainHyperElasticViscoPlastic::_max_substep_iter

protectedinherited

Maximum number of substep iterations.

Definition at line 161 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by FiniteStrainHyperElasticViscoPlastic::computeQpStress().

◆ _maxiters

unsigned int FiniteStrainHyperElasticViscoPlastic::_maxiters

protectedinherited

Maximum number of iterations.

Definition at line 159 of file FiniteStrainHyperElasticViscoPlastic.h.

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

protectedinherited

Number of flow rate user objects.

Definition at line 173 of file FiniteStrainHyperElasticViscoPlastic.h.

◆ _num_int_var_rate_uos

unsigned int FiniteStrainHyperElasticViscoPlastic::_num_int_var_rate_uos

protectedinherited

Number of internal variable rate user objects.

Definition at line 179 of file FiniteStrainHyperElasticViscoPlastic.h.

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

◆ _num_int_var_uos

unsigned int FiniteStrainHyperElasticViscoPlastic::_num_int_var_uos

protectedinherited

Number of internal variable user objects.

Definition at line 177 of file FiniteStrainHyperElasticViscoPlastic.h.

◆ _num_stiffness

bool HyperElasticPhaseFieldIsoDamage::_num_stiffness

protected

Flag to compute numerical stiffness.

Definition at line 42 of file HyperElasticPhaseFieldIsoDamage.h.

Referenced by computePK2StressAndDerivative().

◆ _num_strength_uos

unsigned int FiniteStrainHyperElasticViscoPlastic::_num_strength_uos

protectedinherited

Number of strength user objects.

Definition at line 175 of file FiniteStrainHyperElasticViscoPlastic.h.

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

◆ _pert_val

Real HyperElasticPhaseFieldIsoDamage::_pert_val

protected

Perturbation value for strain components.

Definition at line 56 of file HyperElasticPhaseFieldIsoDamage.h.

Referenced by computeNumStiffness().

◆ _pk2

MaterialProperty<RankTwoTensor>& FiniteStrainHyperElasticViscoPlastic::_pk2

protectedinherited

Definition at line 191 of file FiniteStrainHyperElasticViscoPlastic.h.

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

◆ _pk2_fet

RankTwoTensor FiniteStrainHyperElasticViscoPlastic::_pk2_fet

protectedinherited

Definition at line 214 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by FiniteStrainHyperElasticViscoPlastic::computeQpJacobian().

◆ _pk2_prop_name

std::string FiniteStrainHyperElasticViscoPlastic::_pk2_prop_name

protectedinherited

Definition at line 190 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by FiniteStrainHyperElasticViscoPlastic::computeFlowRateJacobian().

◆ _pk2_tmp

RankTwoTensor HyperElasticPhaseFieldIsoDamage::_pk2_tmp

protected

Definition at line 64 of file HyperElasticPhaseFieldIsoDamage.h.

Referenced by computeDamageStress(), computeNumStiffness(), and computePK2StressAndDerivative().

◆ _resid

DenseVector<Real> FiniteStrainHyperElasticViscoPlastic::_resid

protectedinherited

Definition at line 226 of file FiniteStrainHyperElasticViscoPlastic.h.

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

◆ _resid_abs_tol

Real FiniteStrainHyperElasticViscoPlastic::_resid_abs_tol

protectedinherited

Absolute tolerance for residual convergence check.

Definition at line 155 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by FiniteStrainHyperElasticViscoPlastic::solveFlowrate().

◆ _resid_rel_tol

Real FiniteStrainHyperElasticViscoPlastic::_resid_rel_tol

protectedinherited

Relative tolerance for residual convergence check.

Definition at line 157 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by FiniteStrainHyperElasticViscoPlastic::solveFlowrate().

◆ _rotation_increment

const MaterialProperty<RankTwoTensor>& FiniteStrainHyperElasticViscoPlastic::_rotation_increment

protectedinherited

Definition at line 202 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by FiniteStrainHyperElasticViscoPlastic::computeQpJacobian().

◆ _save_state

bool HyperElasticPhaseFieldIsoDamage::_save_state

protected

Flag to save couple material properties.

Definition at line 60 of file HyperElasticPhaseFieldIsoDamage.h.

Referenced by computeDamageStress(), and computePK2StressAndDerivative().

◆ _strength_prop

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

protectedinherited

Definition at line 205 of file FiniteStrainHyperElasticViscoPlastic.h.

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

◆ _strength_uo

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

protectedinherited

Strength user objects.

Definition at line 184 of file FiniteStrainHyperElasticViscoPlastic.h.

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

◆ _strength_uo_names

std::vector<UserObjectName> FiniteStrainHyperElasticViscoPlastic::_strength_uo_names

protectedinherited

Names of strength user objects.

Definition at line 166 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by FiniteStrainHyperElasticViscoPlastic::computeFlowRateJacobian(), and FiniteStrainHyperElasticViscoPlastic::initUOVariables().

◆ _stress

MaterialProperty<RankTwoTensor>& ComputeGeneralStressBase::_stress

protectedinherited

Stress material property.

Definition at line 48 of file ComputeGeneralStressBase.h.

◆ _tan_mod

RankFourTensor FiniteStrainHyperElasticViscoPlastic::_tan_mod

protectedinherited

Definition at line 217 of file FiniteStrainHyperElasticViscoPlastic.h.

Referenced by FiniteStrainHyperElasticViscoPlastic::computeQpJacobian().

◆ _use_current_hist

bool HyperElasticPhaseFieldIsoDamage::_use_current_hist

protected

Use current value of history variable.

Definition at line 46 of file HyperElasticPhaseFieldIsoDamage.h.

Referenced by computeDamageStress(), and computeQpJacobian().

◆ _zero_pert

Real HyperElasticPhaseFieldIsoDamage::_zero_pert

protected

Perturbation value for near zero or zero strain components.

Definition at line 54 of file HyperElasticPhaseFieldIsoDamage.h.

Referenced by computeNumStiffness().

◆ _zero_tol

Real HyperElasticPhaseFieldIsoDamage::_zero_tol

protected

Used in numerical stiffness calculation to check near zero values.

Definition at line 52 of file HyperElasticPhaseFieldIsoDamage.h.

Referenced by computeNumStiffness().

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

solid_mechanics/include/materials/HyperElasticPhaseFieldIsoDamage.h
solid_mechanics/src/materials/HyperElasticPhaseFieldIsoDamage.C

Public Types

Public Member Functions

Static Public Member Functions

Protected Member Functions

Protected Attributes

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ HyperElasticPhaseFieldIsoDamage()

Member Function Documentation

◆ computeDamageStress()

◆ computeDeeDce()

◆ computeDpk2Dfpinv()

◆ computeElasticPlasticDeformGrad()

◆ computeElasticRightCauchyGreenTensor()

◆ computeElasticStrain()

◆ computeFlowDirection()

◆ computeFlowRateFunction()

◆ computeFlowRateJacobian()

◆ computeFlowRateResidual()

◆ computeIntVar()

◆ computeIntVarDerivatives()

◆ computeIntVarRateDerivatives()

◆ computeIntVarRates()

◆ computeNorm()

◆ computeNumStiffness()

◆ computePK2StressAndDerivative()

◆ computeQpJacobian()

◆ computeQpProperties()

◆ computeQpStress()

◆ computeStrength()

◆ computeStrengthDerivatives()

◆ initJacobianVariables()

◆ initNumUserObjects()

◆ initProp()

◆ initPropOld()

◆ initQpStatefulProperties()

◆ initUOVariables()

◆ initUserObjects()

◆ postSolveFlowrate()

◆ postSolveQp()

◆ preSolveFlowrate()

◆ preSolveQp()

◆ recoverOldState()

◆ saveOldState()

◆ solveFlowrate()

◆ solveQp()

◆ updateFlowRate()

◆ validParams()

Member Data Documentation

◆ _base_name

◆ _c

◆ _ce

◆ _d2Fdc2

◆ _d2Fdcdstrain

◆ _dce_dfe

◆ _dee_dce

◆ _deformation_gradient

◆ _deformation_gradient_old

◆ _df_dstretch_inc

◆ _dFdc

◆ _dfe_df

◆ _dfe_dfpinv

◆ _dfgrd_tmp

◆ _dflow_rate

◆ _dflowrate_dpk2

◆ _dflowrate_dstrength

◆ _dfpinv_dflowrate

◆ _dG0_dee

◆ _dintvar_dflowrate

◆ _dintvar_dflowrate_tmp

◆ _dintvar_dintvar

◆ _dintvar_dintvar_x

◆ _dintvar_dintvarrate

◆ _dintvarrate_dflowrate

◆ _dintvarrate_dintvar

◆ _dpk2_dc

◆ _dpk2_dce

◆ _dpk2_dee