FiniteStrainUObasedCP.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"
 
#include "CrystalPlasticitySlipRate.h"
#include "CrystalPlasticitySlipResistance.h"
#include "CrystalPlasticityStateVariable.h"
#include "CrystalPlasticityStateVarRateComponent.h"
 
class FiniteStrainUObasedCP;
 
template <>
InputParameters validParams<FiniteStrainUObasedCP>();
 
class FiniteStrainUObasedCP : public ComputeStressBase
{
public:
  static InputParameters validParams();
 
  FiniteStrainUObasedCP(const InputParameters & parameters);
 
protected:
  virtual void computeQpStress();
 
  virtual void initQpStatefulProperties();
 
  virtual void calcResidJacob();
 
  virtual void updateSlipSystemResistanceAndStateVariable();
 
  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();
 
  virtual void calcJacobian();
 
  virtual void getSlipRates();
 
  virtual void calcTangentModuli();
 
  virtual void elasticTangentModuli();
 
  virtual void elastoPlasticTangentModuli();
 
  bool lineSearchUpdate(const Real rnorm_prev, const RankTwoTensor);
 
  virtual bool isStateVariablesConverged();
 
  std::vector<const CrystalPlasticitySlipRate *> _uo_slip_rates;
 
  std::vector<const CrystalPlasticitySlipResistance *> _uo_slip_resistances;
 
  std::vector<const CrystalPlasticityStateVariable *> _uo_state_vars;
 
  std::vector<const CrystalPlasticityStateVarRateComponent *> _uo_state_var_evol_rate_comps;
 
  std::vector<MaterialProperty<std::vector<Real>> *> _mat_prop_slip_rates;
 
  std::vector<MaterialProperty<std::vector<Real>> *> _mat_prop_slip_resistances;
 
  std::vector<MaterialProperty<std::vector<Real>> *> _mat_prop_state_vars;
 
  std::vector<const MaterialProperty<std::vector<Real>> *> _mat_prop_state_vars_old;
 
  std::vector<MaterialProperty<std::vector<Real>> *> _mat_prop_state_var_evol_rate_comps;
 
  unsigned int _num_uo_slip_rates;
 
  unsigned int _num_uo_slip_resistances;
 
  unsigned int _num_uo_state_vars;
 
  unsigned int _num_uo_state_var_evol_rate_comps;
 
  std::vector<std::vector<Real>> _state_vars_old;
 
  std::vector<std::vector<Real>> _state_vars_old_stored;
 
  std::vector<std::vector<Real>> _state_vars_prev;
 
  Real _rtol;
  Real _abs_tol;
  Real _stol;
  Real _zero_tol;
 
  RankTwoTensor _resid;
 
  RankFourTensor _jac;
 
  unsigned int _maxiter;
  unsigned int _maxiterg;
 
  MooseEnum _tan_mod_type;
 
  unsigned int _max_substep_iter;
 
  bool _use_line_search;
 
  Real _min_lsrch_step;
 
  Real _lsrch_tol;
 
  unsigned int _lsrch_max_iter;
 
  MooseEnum _lsrch_method;
 
  MaterialProperty<RankTwoTensor> & _fp;
  const MaterialProperty<RankTwoTensor> & _fp_old;
  MaterialProperty<RankTwoTensor> & _pk2;
  const MaterialProperty<RankTwoTensor> & _pk2_old;
  MaterialProperty<RankTwoTensor> & _lag_e;
  MaterialProperty<RankTwoTensor> & _update_rot;
  const MaterialProperty<RankTwoTensor> & _update_rot_old;
 
  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> & _crysrot;
 
  RankTwoTensor _dfgrd_tmp;
  RankTwoTensor _fe, _fp_old_inv, _fp_inv;
  DenseVector<Real> _tau;
  std::vector<MaterialProperty<std::vector<RankTwoTensor>> *> _flow_direction;
 
  bool _err_tol;
 
  RankTwoTensor _delta_dfgrd, _dfgrd_tmp_old;
  Real _dfgrd_scale_factor;
};