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