CrystalPlasticityHCPDislocationSlipBeyerleinUpdate.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 "CrystalPlasticityStressUpdateBase.h"

class CrystalPlasticityHCPDislocationSlipBeyerleinUpdate;

class CrystalPlasticityHCPDislocationSlipBeyerleinUpdate : public CrystalPlasticityStressUpdateBase
{
public:
  static InputParameters validParams();

  CrystalPlasticityHCPDislocationSlipBeyerleinUpdate(const InputParameters & parameters);

protected:
  virtual void initQpStatefulProperties() override;
  virtual void setMaterialVectorSize() override;

  virtual void setInitialConstitutiveVariableValues() override;

  virtual void setSubstepConstitutiveVariableValues() override;

  virtual void updateSubstepConstitutiveVariableValues() override;

  virtual bool calculateSlipRate() override;

  virtual void
  calculateEquivalentSlipIncrement(RankTwoTensor & /*equivalent_slip_increment*/) override;

  virtual void calculateConstitutiveSlipDerivative(std::vector<Real> & dslip_dtau) override;

  virtual void cacheStateVariablesBeforeUpdate() override;

  virtual void calculateStateVariableEvolutionRateComponent() override;

  virtual bool updateStateVariables() override;

  void calculateForestDislocationEvolutionIncrement();

  void calculateSubstructureDensityEvolutionIncrement();

  bool calculateForestDislocationDensity();

  bool calculateSubstructureDislocationDensity();

  virtual void calculateSlipResistance() override;

  virtual void calculateGrainSizeResistance(DenseVector<Real> & lattice_resistance);

  virtual bool areConstitutiveStateVariablesConverged() override;

  bool isSubstructureDislocationDensityConverged();

  const VariableValue & _temperature;

  MaterialProperty<std::vector<Real>> & _forest_dislocation_density;
  const MaterialProperty<std::vector<Real>> & _forest_dislocation_density_old;
  MaterialProperty<std::vector<Real>> & _forest_dislocation_increment;
  MaterialProperty<std::vector<Real>> & _forest_dislocations_removed_increment;
  const Real _initial_forest_dislocation_density;

  MaterialProperty<Real> & _total_substructure_density;
  const MaterialProperty<Real> & _total_substructure_density_old;
  MaterialProperty<Real> & _total_substructure_density_increment;
  const Real _initial_substructure_density;

  const unsigned int _slip_system_modes;

  const std::vector<unsigned int> _number_slip_systems_per_mode;

  const std::vector<Real> _lattice_friction;

  const Real _reference_strain_rate;
  const Real _rate_sensitivity_exponent;

  const std::vector<Real> _burgers_vector;
  const std::vector<Real> _slip_generation_coefficient;
  const std::vector<Real> _slip_activation_energy;
  const std::vector<Real> _proportionality_factor;
  const Real _forest_interaction_coefficient;
  const Real _boltzman_constant;
  const Real _macro_applied_strain_rate;
  const Real _macro_reference_strain_rate;

  const std::vector<Real> _shear_modulus;

  const std::vector<Real> _substructure_rate_coefficient;

  const Real _substructure_hardening_coefficient;

  const std::vector<Real> _hallpetch_like_coefficient;

  const Real _grain_size;

  std::vector<Real> _previous_substep_slip_resistance;
  std::vector<Real> _previous_substep_forest_dislocations;
  Real _previous_substep_total_substructure_density;

  std::vector<Real> _slip_resistance_before_update;
  std::vector<Real> _forest_dislocations_before_update;
  Real _total_substructure_density_before_update;

  const bool _include_twinning_in_Lp;

  const MaterialProperty<Real> * const _twin_volume_fraction_total;
};