JIntegral.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;
class JIntegral : public ElementVectorPostprocessor
{
public:
  static InputParameters validParams();
 
  JIntegral(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;
 
protected:
  Real computeQpIntegral(const std::size_t crack_front_point_index,
                         const Real scalar_q,
                         const RealVectorValue & grad_of_scalar_q);
  const CrackFrontDefinition * const _crack_front_definition;
  const enum class IntegralType { JIntegral, CIntegral, KFromJIntegral } _integral;
 
  const MaterialProperty<RealVectorValue> * const _J_thermal_term_vec;
 
  const MaterialProperty<RankTwoTensor> * _Eshelby_tensor;
 
  const MaterialProperty<RankTwoTensor> * _Eshelby_tensor_dissipation;
 
  bool _treat_as_2d;
  bool _using_mesh_cutter;
  std::vector<MooseVariableFEBase *> _fe_vars;
  const libMesh::FEType & _fe_type;
  bool _has_symmetry_plane;
  Real _poissons_ratio;
  Real _youngs_modulus;
  std::size_t _ring_index;
  const enum class QMethod { Geometry, Topology } _q_function_type;
 
  const enum class PositionType { Angle, Distance } _position_type;
 
  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;
  unsigned int _qp;
  VectorPostprocessorValue & _x;
  VectorPostprocessorValue & _y;
  VectorPostprocessorValue & _z;
  VectorPostprocessorValue & _position;
  VectorPostprocessorValue & _j_integral;
};