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