PeridynamicsMesh.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 "MooseMesh.h"
#include "MooseApp.h"

#include "libmesh/point.h"

class PeridynamicsMesh : public MooseMesh
{
public:
  static InputParameters validParams();

  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> bonding_block_pairs,
                                  std::multimap<SubdomainID, SubdomainID> non_bonding_block_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> getBondDeformationGradientNeighbors(dof_id_type node_id,
                                                               dof_id_type neighbor_id);

  SubdomainID getNodeBlockID(dof_id_type node_id);

  void setNodeBlockID(SubdomainID id);

  Point getNodeCoord(dof_id_type node_id);

  std::vector<dof_id_type> getPDNodeIDToFEElemIDMap();

  Real getNodeVolume(dof_id_type node_id);

  Real getHorizonVolume(dof_id_type node_id);

  Real getHorizonSubsetVolume(dof_id_type node_id, dof_id_type neighbor_id);

  Real getHorizonSubsetVolumeSum(dof_id_type node_id);

  Real getHorizonSubsetVolumeFraction(dof_id_type node_id, dof_id_type neighbor_id);

  Real getNodeAverageSpacing(dof_id_type node_id);

  Real getHorizon(dof_id_type node_id);

  Real getBoundaryOffset(dof_id_type node_id);

  Real getNeighborWeight(dof_id_type node_id, dof_id_type neighbor_id);

protected:
  const Real _horizon_radius;
  const bool _has_horizon_number;
  const Real _horizon_number;
  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_average_spacing;
  std::vector<Real> & _pdnode_horizon_radius;
  std::vector<Real> & _pdnode_vol;
  std::vector<Real> & _pdnode_horizon_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>> & _pdnode_sub_vol;

  std::vector<Real> & _pdnode_sub_vol_sum;

  std::map<dof_id_type, Real> & _boundary_node_offset;

  std::vector<Real> & _pdnode_weight_normalizer;

  void createNodeHorizBasedData(std::multimap<SubdomainID, SubdomainID> bonding_block_pairs,
                                std::multimap<SubdomainID, SubdomainID> non_bonding_block_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);
};