FiniteStrainHyperElasticViscoPlastic.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html
 
#pragma once
 
#include "ComputeStressBase.h"
#include "HEVPFlowRateUOBase.h"
#include "HEVPStrengthUOBase.h"
#include "HEVPInternalVarUOBase.h"
#include "HEVPInternalVarRateUOBase.h"
 
class FiniteStrainHyperElasticViscoPlastic : public ComputeStressBase
{
public:
  static InputParameters validParams();
 
  FiniteStrainHyperElasticViscoPlastic(const InputParameters & parameters);
 
protected:
  virtual void initUOVariables();
 
  void initNumUserObjects(const std::vector<UserObjectName> &, unsigned int &);
 
  template <typename T>
  void
  initProp(const std::vector<UserObjectName> &, unsigned int, std::vector<MaterialProperty<T> *> &);
 
  template <typename T>
  void initPropOld(const std::vector<UserObjectName> &,
                   unsigned int,
                   std::vector<const MaterialProperty<T> *> &);
 
  template <typename T>
  void initUserObjects(const std::vector<UserObjectName> &, unsigned int, std::vector<const T *> &);
 
  virtual void initJacobianVariables();
 
  virtual void initQpStatefulProperties();
 
  virtual void computeQpStress();
 
  virtual void computeQpJacobian();
 
  virtual void saveOldState();
 
  virtual void preSolveQp();
 
  virtual bool solveQp();
 
  virtual void postSolveQp();
 
  virtual void recoverOldState();
 
  virtual void preSolveFlowrate();
 
  virtual bool solveFlowrate();
 
  virtual void postSolveFlowrate();
 
  virtual bool computeFlowRateFunction();
 
  virtual bool computeFlowDirection();
 
  virtual void computeElasticRightCauchyGreenTensor();
 
  virtual void computePK2StressAndDerivative();
 
  virtual void computeElasticStrain();
 
  virtual void computeDeeDce();
 
  virtual bool computeFlowRateResidual();
 
  virtual void computeFlowRateJacobian();
 
  virtual void computeElasticPlasticDeformGrad();
 
  virtual Real computeNorm(const std::vector<Real> &);
 
  virtual void updateFlowRate();
 
  virtual void computeDpk2Dfpinv();
 
  virtual bool computeIntVarRates();
 
  virtual bool computeIntVar();
 
  virtual bool computeStrength();
 
  virtual void computeIntVarRateDerivatives();
 
  virtual void computeIntVarDerivatives();
 
  void computeStrengthDerivatives();
 
  Real _resid_abs_tol;
  Real _resid_rel_tol;
  unsigned int _maxiters;
  unsigned int _max_substep_iter;
 
  std::vector<UserObjectName> _flow_rate_uo_names;
  std::vector<UserObjectName> _strength_uo_names;
  std::vector<UserObjectName> _int_var_uo_names;
  std::vector<UserObjectName> _int_var_rate_uo_names;
 
  unsigned int _num_flow_rate_uos;
  unsigned int _num_strength_uos;
  unsigned int _num_int_var_uos;
  unsigned int _num_int_var_rate_uos;
 
  std::vector<const HEVPFlowRateUOBase *> _flow_rate_uo;
  std::vector<const HEVPStrengthUOBase *> _strength_uo;
  std::vector<const HEVPInternalVarUOBase *> _int_var_uo;
  std::vector<const HEVPInternalVarRateUOBase *> _int_var_rate_uo;
 
  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;
  const MaterialProperty<RankFourTensor> & _elasticity_tensor;
  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, _fp_tmp_old_inv;
  RankTwoTensor _fe, _ee;
  RankTwoTensor _pk2_fet, _fe_pk2;
  RankFourTensor _dpk2_dce, _dpk2_dfe, _dfe_dfpinv, _dpk2_dfpinv;
  RankFourTensor _dee_dce, _dce_dfe, _dfe_df;
  RankFourTensor _tan_mod, _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;
  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;
};