JIntegral.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://www.mooseframework.org
//*
//* 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 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<RealGradient>> * _dphi_curr_elem;
  unsigned int _qp;
  VectorPostprocessorValue & _x;
  VectorPostprocessorValue & _y;
  VectorPostprocessorValue & _z;
  VectorPostprocessorValue & _position;
  VectorPostprocessorValue & _j_integral;
};