InteractionIntegral.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 "ElementVectorPostprocessor.h"
 
// Forward Declarations
class CrackFrontDefinition;
 
template <bool is_ad>
class InteractionIntegralTempl : public ElementVectorPostprocessor
{
public:
  static InputParameters validParams();
 
  InteractionIntegralTempl(const InputParameters & parameters);
 
  virtual void initialSetup() override;
  virtual void initialize() override;
  virtual void execute() override;
  virtual void finalize() override;
  virtual void threadJoin(const UserObject & y) override;
 
  static MooseEnum qFunctionType();
 
  static MooseEnum sifModeType();
 
protected:
  Real computeQpIntegral(const std::size_t crack_front_point_index,
                         const Real scalar_q,
                         const RealVectorValue & grad_of_scalar_q);
 
  void computeAuxFields(RankTwoTensor & aux_stress,
                        RankTwoTensor & grad_disp,
                        RankTwoTensor & aux_strain,
                        RankTwoTensor & aux_disp);
  void computeTFields(RankTwoTensor & aux_stress, RankTwoTensor & grad_disp);
 
  std::size_t _ndisp;
  const CrackFrontDefinition * const _crack_front_definition;
  bool _treat_as_2d;
  bool _using_mesh_cutter;
  const GenericMaterialProperty<RankTwoTensor, is_ad> & _stress;
  const GenericMaterialProperty<RankTwoTensor, is_ad> & _strain;
  std::vector<MooseVariableFEBase *> _fe_vars;
  const libMesh::FEType & _fe_type;
  std::vector<const VariableValue *> _disp;
  std::vector<const VariableGradient *> _grad_disp;
  const bool _has_temp;
  const VariableGradient & _grad_temp;
  const MaterialProperty<Real> * _functionally_graded_youngs_modulus_crack_dir_gradient;
  const MaterialProperty<Real> * _functionally_graded_youngs_modulus;
  Real _K_factor;
  bool _has_symmetry_plane;
  Real _poissons_ratio;
  Real _youngs_modulus;
  const bool _fgm_crack;
  std::size_t _ring_index;
  const GenericMaterialProperty<RankTwoTensor, is_ad> * const _total_deigenstrain_dT;
  bool _has_additional_eigenstrain;
  const VariableGradient * _additional_eigenstrain_gradient_00;
  const VariableGradient * _additional_eigenstrain_gradient_01;
  const VariableGradient * _additional_eigenstrain_gradient_11;
  const VariableGradient * _additional_eigenstrain_gradient_22;
  const VariableGradient * _additional_eigenstrain_gradient_02;
  const VariableGradient * _additional_eigenstrain_gradient_12;
  std::vector<Real> _q_curr_elem;
  const std::vector<std::vector<Real>> * _phi_curr_elem;
  const std::vector<std::vector<libMesh::RealGradient>> * _dphi_curr_elem;
  Real _kappa;
  Real _shear_modulus;
  Real _r;
  Real _theta;
  unsigned int _qp;
 
  const enum class QMethod { Geometry, Topology } _q_function_type;
 
  const enum class PositionType { Angle, Distance } _position_type;
 
  const enum class SifMethod { KI, KII, KIII, T } _sif_mode;
 
  VectorPostprocessorValue & _x;
  VectorPostprocessorValue & _y;
  VectorPostprocessorValue & _z;
  VectorPostprocessorValue & _position;
  VectorPostprocessorValue & _interaction_integral;
 
  const GenericMaterialProperty<RankThreeTensor, is_ad> * _eigenstrain_gradient;
  const MaterialProperty<RealVectorValue> * _body_force;
};
 
typedef InteractionIntegralTempl<false> InteractionIntegral;
typedef InteractionIntegralTempl<true> ADInteractionIntegral;