MeshCut2DUserObjectBase.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 "MeshCutUserObjectBase.h"
class MeshCut2DNucleationBase;
class CrackFrontDefinition;
class MeshCut2DUserObjectBase : public MeshCutUserObjectBase
{
public:
  static InputParameters validParams();
 
  MeshCut2DUserObjectBase(const InputParameters & parameters);
 
  // initialSetup needs to be called by every derived class
  virtual void initialSetup() override final;
 
  virtual bool cutElementByGeometry(const Elem * elem,
                                    std::vector<Xfem::CutEdge> & cut_edges,
                                    std::vector<Xfem::CutNode> & cut_nodes) const override;
  virtual bool cutElementByGeometry(const Elem * elem,
                                    std::vector<Xfem::CutFace> & cut_faces) const override;
  virtual bool cutFragmentByGeometry(std::vector<std::vector<Point>> & frag_edges,
                                     std::vector<Xfem::CutEdge> & cut_edges) const override;
  virtual bool cutFragmentByGeometry(std::vector<std::vector<Point>> & frag_faces,
                                     std::vector<Xfem::CutFace> & cut_faces) const override;
  virtual unsigned int getNumberOfCrackFrontPoints() const override;
 
  virtual const std::vector<Point>
  getCrackFrontPoints(unsigned int num_crack_front_points) const override;
 
  virtual const std::vector<RealVectorValue>
  getCrackPlaneNormals(unsigned int num_crack_front_points) const override;
 
protected:
  MooseMesh & _mesh;
 
  const MeshCut2DNucleationBase * _nucleate_uo;
 
  bool _is_mesh_modified;
 
  std::vector<std::pair<dof_id_type, dof_id_type>> _original_and_current_front_node_ids;
 
  std::vector<std::pair<dof_id_type, Point>> _active_front_node_growth_vectors;
 
  CrackFrontDefinition * _crack_front_definition;
 
  virtual void findActiveBoundaryGrowth() = 0;
 
  void findOriginalCrackFrontNodes();
 
  void growFront();
 
  void addNucleatedCracksToMesh();
 
private:
  void removeNucleatedCracksTooCloseToEachOther(
      std::map<unsigned int, std::pair<RealVectorValue, RealVectorValue>> & nucleated_elems_map,
      Real nucleationRadius);
  void removeNucleatedCracksTooCloseToExistingCracks(
      std::map<unsigned int, std::pair<RealVectorValue, RealVectorValue>> & nucleated_elems_map,
      Real nucleationRadius);
};