PeridynamicsMesh.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 "MooseMesh.h"
#include "MooseApp.h"
 
#include "libmesh/point.h"
 
class PeridynamicsMesh;
 
template <>
InputParameters validParams<PeridynamicsMesh>();
 
class PeridynamicsMesh : public MooseMesh
{
public:
  PeridynamicsMesh(const InputParameters & parameters);
  PeridynamicsMesh(const PeridynamicsMesh & /* other_mesh */) = default;
 
  PeridynamicsMesh & operator=(const PeridynamicsMesh & other_mesh) = delete;
  virtual std::unique_ptr<MooseMesh> safeClone() const override;
  virtual void buildMesh() override;
  virtual unsigned int dimension() const override;
 
  dof_id_type nPDNodes() const;
 
  dof_id_type nPDBonds() const;
 
  void
  createPeridynamicsMeshData(MeshBase & fe_mesh,
                             std::set<dof_id_type> converted_elem_id,
                             std::multimap<SubdomainID, SubdomainID> connect_block_id_pairs,
                             std::multimap<SubdomainID, SubdomainID> non_connect_block_id_pairs);
 
  std::vector<dof_id_type> getNeighbors(dof_id_type node_id);
 
  dof_id_type getNeighborIndex(dof_id_type node_i, dof_id_type node_j);
 
  std::vector<dof_id_type> getBonds(dof_id_type node_id);
 
  std::vector<dof_id_type> getDefGradNeighbors(dof_id_type node_id, unsigned int neighbor_id);
 
  SubdomainID getNodeBlockID(dof_id_type node_id);
 
  void setNodeBlockID(SubdomainID id);
 
  Point getPDNodeCoord(dof_id_type node_id);
 
  std::vector<dof_id_type> getPDNodeIDToFiniteElemIDMap();
 
  Real getPDNodeVolume(dof_id_type node_id);
 
  Real getHorizVolume(dof_id_type node_id);
 
  Real getDefGradVolFraction(dof_id_type node_id, dof_id_type neighbor_id);
 
  Real getNodeAvgSpacing(dof_id_type node_id);
 
  Real getHorizon(dof_id_type node_id);
 
  Real getBoundaryOffset(dof_id_type node_id);
 
protected:
  const Real _horiz_rad;
  const bool _has_horiz_num;
  const Real _horiz_num;
  const Real _bah_ratio;
 
  const bool _has_cracks;
  std::vector<Point> _cracks_start;
  std::vector<Point> _cracks_end;
  std::vector<Real> _cracks_width;
 
  unsigned int & _dim;
 
  unsigned int & _n_pdnodes;
 
  unsigned int & _n_pdbonds;
 
  std::vector<Point> _pdnode_coord;
  std::vector<Real> & _pdnode_avg_spacing;
  std::vector<Real> & _pdnode_horiz_rad;
  std::vector<Real> & _pdnode_vol;
  std::vector<Real> & _pdnode_horiz_vol;
  std::vector<SubdomainID> & _pdnode_blockID;
  std::vector<dof_id_type> & _pdnode_elemID;
 
  std::vector<std::vector<dof_id_type>> & _pdnode_neighbors;
 
  std::vector<std::vector<dof_id_type>> & _pdnode_bonds;
 
  std::vector<std::vector<std::vector<dof_id_type>>> & _dg_neighbors;
 
  std::vector<std::vector<Real>> & _dg_vol_frac;
 
  std::map<dof_id_type, Real> & _boundary_node_offset;
 
  void createNodeHorizBasedData(std::multimap<SubdomainID, SubdomainID> connect_block_id_pairs,
                                std::multimap<SubdomainID, SubdomainID> non_connect_block_id_pairs);
 
  bool checkInterface(SubdomainID pdnode_blockID_i,
                      SubdomainID pdnode_blockID_j,
                      std::multimap<SubdomainID, SubdomainID> blockID_pairs);
 
  void createNeighborHorizonBasedData();
 
  bool
  checkPointInsideRectangle(Point point, Point rec_p1, Point rec_p2, Real rec_height, Real tol = 0);
 
  bool checkCrackIntersectBond(
      Point crack_p1, Point crack_p2, Real crack_width, Point bond_p1, Point bond_p2);
 
  bool checkSegmentIntersectSegment(Point seg1_p1, Point seg1_p2, Point seg2_p1, Point seg2_p2);
};