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 "GeometricCutUserObject.h"
class MeshCut2DNucleationBase;
class CrackFrontDefinition;
class MeshCut2DUserObjectBase : public GeometricCutUserObject
{
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 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;

  MeshBase & getCutterMesh() const;

protected:
  MooseMesh & _mesh;

  std::unique_ptr<MeshBase> _cutter_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);
};