ComputeMultipleCrystalPlasticityStress.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 "ComputeFiniteStrainElasticStress.h"
 
#include "CrystalPlasticityStressUpdateBase.h"
#include "ComputeCrystalPlasticityEigenstrainBase.h"
 
#include "RankTwoTensor.h"
#include "RankFourTensor.h"
 
class ComputeMultipleCrystalPlasticityStress : public ComputeFiniteStrainElasticStress
{
public:
  static InputParameters validParams();
 
  ComputeMultipleCrystalPlasticityStress(const InputParameters & parameters);
 
  virtual void initialSetup() override;
 
protected:
  virtual void computeQpStress() override;
 
  virtual void updateStress(RankTwoTensor & cauchy_stress, RankFourTensor & jacobian_mult);
 
  virtual void initQpStatefulProperties() override;
 
  void calculateResidualAndJacobian();
 
  void preSolveQp();
 
  void solveQp();
 
  void postSolveQp(RankTwoTensor & stress_new, RankFourTensor & jacobian_mult);
 
  void solveStateVariables();
 
  void solveStress();
 
  void calculateResidual();
 
  void calculateJacobian();
 
  void calcTangentModuli(RankFourTensor & jacobian_mult);
  void elasticTangentModuli(RankFourTensor & jacobian_mult);
  void elastoPlasticTangentModuli(RankFourTensor & jacobian_mult);
 
  bool lineSearchUpdate(const Real & rnorm_prev, const RankTwoTensor & dpk2);
 
  void calculateEigenstrainDeformationGrad();
 
  const unsigned _num_models;
 
  std::vector<CrystalPlasticityStressUpdateBase *> _models;
 
  const unsigned _num_eigenstrains;
 
  std::vector<ComputeCrystalPlasticityEigenstrainBase *> _eigenstrains;
 
  const std::string _base_name;
 
  const MaterialProperty<RankFourTensor> & _elasticity_tensor;
 
  Real _rtol;
  Real _abs_tol;
 
  RankTwoTensor _residual_tensor;
  RankFourTensor _jacobian;
 
  unsigned int _maxiter;
  unsigned int _maxiterg;
 
  const enum class TangentModuliType { EXACT, NONE } _tan_mod_type;
 
  unsigned int _max_substep_iter;
 
  Real _substep_dt;
 
  bool _use_line_search;
 
  Real _min_line_search_step_size;
 
  Real _line_search_tolerance;
 
  unsigned int _line_search_max_iterations;
 
  const enum class LineSearchMethod { CutHalf, Bisection } _line_search_method;
 
  MaterialProperty<RankTwoTensor> & _plastic_deformation_gradient;
  const MaterialProperty<RankTwoTensor> & _plastic_deformation_gradient_old;
 
  MaterialProperty<RankTwoTensor> * _eigenstrain_deformation_gradient;
  const MaterialProperty<RankTwoTensor> * _eigenstrain_deformation_gradient_old;
 
  const MaterialProperty<RankTwoTensor> & _deformation_gradient;
  const MaterialProperty<RankTwoTensor> & _deformation_gradient_old;
 
  MaterialProperty<RankTwoTensor> & _pk2;
  const MaterialProperty<RankTwoTensor> & _pk2_old;
 
  MaterialProperty<RankTwoTensor> & _total_lagrangian_strain;
 
  MaterialProperty<RankTwoTensor> & _updated_rotation;
  const MaterialProperty<RankTwoTensor> & _updated_rotation_old;
 
  MaterialProperty<Real> & _misorientation;
 
  const MaterialProperty<RankTwoTensor> & _crysrot;
 
  RankTwoTensor _temporary_deformation_gradient;
  RankTwoTensor _elastic_deformation_gradient;
  RankTwoTensor _inverse_plastic_deformation_grad;
  RankTwoTensor _inverse_plastic_deformation_grad_old;
  RankTwoTensor _inverse_eigenstrain_deformation_grad;
 
  const bool _print_convergence_message;
 
  bool _convergence_failed;
 
  RankTwoTensor _delta_deformation_gradient;
  RankTwoTensor _temporary_deformation_gradient_old;
 
  Real _dfgrd_scale_factor;
};