MeshCut2DNucleationBase.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 "ElementUserObject.h"

class XFEM;

class MeshCut2DNucleationBase : public ElementUserObject
{
public:
  static InputParameters validParams();

  MeshCut2DNucleationBase(const InputParameters & parameters);

  virtual void initialize() override;
  virtual void execute() override;
  virtual void threadJoin(const UserObject & y) override;
  virtual void finalize() override;
  std::map<unsigned int, std::pair<RealVectorValue, RealVectorValue>> getNucleatedElemsMap() const
  {
    return _nucleated_elems;
  }
  Real getNucleationRadius() const { return _nucleation_radius; }

protected:
  virtual bool doesElementCrack(std::pair<RealVectorValue, RealVectorValue> & cutterElemNodes) = 0;

private:
  MooseMesh & _mesh;
  // New cracks can only be nucleated if they are at least this far from another crack
  Real _nucleation_radius;
  std::shared_ptr<XFEM> _xfem;
  // Boundaries where cracks can nucleate
  std::vector<BoundaryID> _initiation_boundary_ids;
  // map with key for element id and value is a pair containing the node points for creating a
  // nucleated element on the xfem cutter mesh.
  std::map<unsigned int, std::pair<RealVectorValue, RealVectorValue>> _nucleated_elems;
};