FiniteStrainCrystalPlasticity.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://www.mooseframework.org
//*
//* 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"
 
class FiniteStrainCrystalPlasticity;
 
template <>
InputParameters validParams<FiniteStrainCrystalPlasticity>();
 
class FiniteStrainCrystalPlasticity : public ComputeStressBase
{
public:
  static InputParameters validParams();
 
  FiniteStrainCrystalPlasticity(const InputParameters & parameters);
 
protected:
  virtual void computeQpStress();
 
  virtual void computeQpElasticityTensor();
 
  virtual void initQpStatefulProperties();
 
  virtual void calc_resid_jacob(RankTwoTensor &, RankFourTensor &);
 
  virtual void getSlipIncrements();
 
  // Override to modify slip system resistance evolution
  virtual void update_slip_system_resistance();
 
  // Old function: Kept to avoid code break in computeQpStress
  virtual void updateGss();
 
  virtual void getSlipSystems();
 
  virtual void assignSlipSysRes();
 
  virtual void readFileInitSlipSysRes();
 
  virtual void getInitSlipSysRes();
 
  virtual void readFileFlowRateParams();
 
  virtual void getFlowRateParams();
 
  virtual void readFileHardnessParams();
 
  virtual void getHardnessParams();
 
  virtual void initSlipSysProps();
 
  virtual void initAdditionalProps();
 
  virtual void preSolveQp();
 
  virtual void solveQp();
 
  virtual void postSolveQp();
 
  virtual void preSolveStatevar();
 
  virtual void solveStatevar();
 
  virtual void postSolveStatevar();
 
  virtual void preSolveStress();
 
  virtual void solveStress();
 
  virtual void postSolveStress();
 
  virtual void calcResidual(RankTwoTensor &);
 
  virtual void calcJacobian(RankFourTensor &);
 
  virtual RankFourTensor calcTangentModuli();
 
  virtual RankFourTensor elasticTangentModuli();
 
  virtual RankFourTensor elastoPlasticTangentModuli();
 
  RankTwoTensor get_current_rotation(const RankTwoTensor & a);
 
 
  RankTwoTensor getMatRot(const RankTwoTensor & a);
 
  void calc_schmid_tensor();
 
  bool line_search_update(const Real rnorm_prev, const RankTwoTensor);
 
  void internalVariableUpdateNRiteration();
 
  const unsigned int _nss;
 
  std::vector<Real> _gprops;
  std::vector<Real> _hprops;
  std::vector<Real> _flowprops;
 
  std::string _slip_sys_file_name;
 
  std::string _slip_sys_res_prop_file_name;
 
  std::string _slip_sys_flow_prop_file_name;
 
  std::string _slip_sys_hard_prop_file_name;
 
  Real _rtol;
  Real _abs_tol;
  Real _gtol;
  Real _slip_incr_tol;
 
  unsigned int _maxiter;
  unsigned int _maxiterg;
 
  unsigned int _num_slip_sys_flowrate_props;
 
  MooseEnum _tan_mod_type;
 
  MooseEnum _intvar_read_type;
 
  unsigned int _num_slip_sys_props;
 
  bool _gen_rndm_stress_flag;
 
  bool _input_rndm_scale_var;
 
  Real _rndm_scale_var;
 
  unsigned int _rndm_seed;
 
  unsigned int _max_substep_iter;
 
  bool _use_line_search;
 
  Real _min_lsrch_step;
 
  Real _lsrch_tol;
 
  unsigned int _lsrch_max_iter;
 
  // Line search method
  MooseEnum _lsrch_method;
 
  MaterialProperty<RankTwoTensor> & _fp;
  const MaterialProperty<RankTwoTensor> & _fp_old;
  MaterialProperty<RankTwoTensor> & _pk2;
  const MaterialProperty<RankTwoTensor> & _pk2_old;
  MaterialProperty<RankTwoTensor> & _lag_e;
  const MaterialProperty<RankTwoTensor> & _lag_e_old;
  MaterialProperty<std::vector<Real>> & _gss;
  const MaterialProperty<std::vector<Real>> & _gss_old;
  MaterialProperty<Real> & _acc_slip;
  const MaterialProperty<Real> & _acc_slip_old;
  MaterialProperty<RankTwoTensor> & _update_rot;
 
  const MaterialProperty<RankTwoTensor> & _deformation_gradient;
  const MaterialProperty<RankTwoTensor> & _deformation_gradient_old;
  const std::string _elasticity_tensor_name;
  const MaterialProperty<RankFourTensor> & _elasticity_tensor;
  const MaterialProperty<RankTwoTensor> & _crysrot;
 
  DenseVector<Real> _mo;
  DenseVector<Real> _no;
 
  DenseVector<Real> _a0;
  DenseVector<Real> _xm;
 
  Real _h0;
  Real _tau_sat;
  Real _tau_init;
  Real _r;
 
  RankTwoTensor _dfgrd_tmp;
  RankTwoTensor _fe, _fp_old_inv, _fp_inv, _fp_prev_inv;
  DenseVector<Real> _slip_incr, _tau, _dslipdtau;
  std::vector<RankTwoTensor> _s0;
 
  RankTwoTensor _pk2_tmp, _pk2_tmp_old;
  Real _accslip_tmp, _accslip_tmp_old;
  std::vector<Real> _gss_tmp;
  std::vector<Real> _gss_tmp_old;
 
  DenseVector<Real> _slip_sys_props;
 
  DenseMatrix<Real> _dgss_dsliprate;
 
  bool _read_from_slip_sys_file;
 
  bool _err_tol; 
 
  RankTwoTensor _delta_dfgrd, _dfgrd_tmp_old;
  Real _dfgrd_scale_factor;
  bool _first_step_iter, _last_step_iter, _first_substep;
};