Line data    Source code

       1             : //* This file is part of the MOOSE framework
       2             : //* https://mooseframework.inl.gov
       3             : //*
       4             : //* All rights reserved, see COPYRIGHT for full restrictions
       5             : //* https://github.com/idaholab/moose/blob/master/COPYRIGHT
       6             : //*
       7             : //* Licensed under LGPL 2.1, please see LICENSE for details
       8             : //* https://www.gnu.org/licenses/lgpl-2.1.html
       9             : 
      10             : #pragma once
      11             : 
      12             : #include "AuxKernel.h"
      13             : #include "CrackFrontDefinition.h"
      14             : 
      15             : /**
      16             :  * Coupled auxiliary value
      17             :  */
      18             : class DomainIntegralQFunction : public AuxKernel
      19             : {
      20             : public:
      21             :   static InputParameters validParams();
      22             : 
      23             :   /**
      24             :    * Factory constructor, takes parameters so that all derived classes can be built using the same
      25             :    * constructor.
      26             :    */
      27             :   DomainIntegralQFunction(const InputParameters & parameters);
      28             : 
      29        2496 :   virtual ~DomainIntegralQFunction() {}
      30             : 
      31             : protected:
      32             :   virtual void initialSetup();
      33             :   virtual Real computeValue();
      34             :   void projectToFrontAtPoint(Real & dist_to_front, Real & dist_along_tangent);
      35             : 
      36             : private:
      37             :   const Real _j_integral_radius_inner;
      38             :   const Real _j_integral_radius_outer;
      39             :   const CrackFrontDefinition * const _crack_front_definition;
      40             :   bool _has_crack_front_point_index;
      41             :   const unsigned int _crack_front_point_index;
      42             :   bool _treat_as_2d;
      43             :   bool _is_point_on_intersecting_boundary;
      44             : };