FractureMechanics::CrackPropagation::FaceOrientation Struct Reference

Determine face orientation. More...

#include <users_modules/fracture_mechanics/src/CrackPropagation.hpp>

Inheritance diagram for FractureMechanics::CrackPropagation::FaceOrientation:

Collaboration diagram for FractureMechanics::CrackPropagation::FaceOrientation:

Public Member Functions
	FaceOrientation (bool crack_front, Range contact_faces, BitRefLevel bit_level)
virtual	~FaceOrientation ()
MoFEMErrorCode	getElementOrientation (MoFEM::Interface &m_field, const EntityHandle face, const BitRefLevel bit)
MoFEMErrorCode	getElementOrientation (MoFEM::Interface &m_field, const FEMethod *fe_method_ptr)
MoFEMErrorCode	getElementOrientation (MoFEM::Interface &m_field, const EntityHandle face)

Public Attributes
bool	useProjectionFromCrackFront
Range	contactFaces
BitRefLevel	defaultBitLevel
Public Attributes inherited from SurfaceSlidingConstrains::DriverElementOrientation
int	elementOrientation

Detailed Description

Determine face orientation.

Definition at line 1003 of file CrackPropagation.hpp.

Constructor & Destructor Documentation

FaceOrientation()

FractureMechanics::CrackPropagation::FaceOrientation::FaceOrientation ( bool  crack_front, Range  contact_faces, BitRefLevel  bit_level  )

inline

Definition at line 1008 of file CrackPropagation.hpp.

        : useProjectionFromCrackFront(crack_front), contactFaces(contact_faces),
          defaultBitLevel(bit_level) {}

~FaceOrientation()

virtual FractureMechanics::CrackPropagation::FaceOrientation::~FaceOrientation ( )

inlinevirtual

Definition at line 1012 of file CrackPropagation.hpp.

{}

Member Function Documentation

getElementOrientation() [1/3]

MoFEMErrorCode FractureMechanics::CrackPropagation::FaceOrientation::getElementOrientation ( MoFEM::Interface &  m_field, const EntityHandle  face  )

inlinevirtual

Reimplemented from SurfaceSlidingConstrains::DriverElementOrientation.

Definition at line 1024 of file CrackPropagation.hpp.

                                                                  {
      MoFEMFunctionBegin;
      CHKERR getElementOrientation(m_field, face, defaultBitLevel);
      MoFEMFunctionReturn(0);
    }

getElementOrientation() [2/3]

MoFEMErrorCode FractureMechanics::CrackPropagation::FaceOrientation::getElementOrientation	(	MoFEM::Interface &	m_field,
		const EntityHandle	face,
		const BitRefLevel	bit
	)

Definition at line 9737 of file CrackPropagation.cpp.

                                                                             {
  MoFEMFunctionBegin;
  elementOrientation = +1;
  Range adj_side_elems;
  CHKERR m_field.getInterface<BitRefManager>()->getAdjacencies(
      bit, &face, 1, 3, adj_side_elems, moab::Interface::INTERSECT, QUIET);
  adj_side_elems = adj_side_elems.subset_by_type(MBTET);
  if (adj_side_elems.empty()) {
    Range adj_tets_on_surface;
    BitRefLevel bit_tet_on_surface;
    bit_tet_on_surface.set(BITREFLEVEL_SIZE - 2);
    CHKERR m_field.getInterface<BitRefManager>()->getAdjacencies(
        bit_tet_on_surface, &face, 1, 3, adj_tets_on_surface,
        moab::Interface::INTERSECT, 0);
    adj_side_elems.insert(*adj_tets_on_surface.begin());
  }
  if (adj_side_elems.size() != 1) {
    adj_side_elems.clear();
    CHKERR m_field.getInterface<BitRefManager>()->getAdjacencies(
        bit, &face, 1, 3, adj_side_elems, moab::Interface::INTERSECT,
        VERY_NOISY);
    Range::iterator it = adj_side_elems.begin();
    for (; it != adj_side_elems.end(); it++) {
      Range nodes;
      CHKERR m_field.get_moab().get_connectivity(&*it, 1, nodes, true);
      PetscPrintf(PETSC_COMM_WORLD, "Connectivity %lu %lu %lu %lu\n", nodes[0],
                  nodes[1], nodes[2], nodes[3]);
    }
    int rank;
    MPI_Comm_rank(m_field.get_comm(), &rank);
    if (rank == 0) {
      EntityHandle out_meshset;
      CHKERR m_field.get_moab().create_meshset(MESHSET_SET, out_meshset);
      CHKERR m_field.get_moab().add_entities(out_meshset, adj_side_elems);
      CHKERR m_field.get_moab().add_entities(out_meshset, &face, 1);
      CHKERR m_field.get_moab().write_file("debug_error.vtk", "VTK", "",
                                           &out_meshset, 1);
      CHKERR m_field.get_moab().delete_entities(&out_meshset, 1);
    }
    SETERRQ1(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
             "Expect 1 tet but is %u", adj_side_elems.size());
  }
  EntityHandle side_elem = adj_side_elems[0];
  if (side_elem != 0) {
    int side_number, sense, offset;
    CHKERR m_field.get_moab().side_number(side_elem, face, side_number, sense,
                                          offset);
    if (sense == -1) {
      elementOrientation = -1;
    }
  }
  if (useProjectionFromCrackFront || contactFaces.contains(Range(face, face))) {
    Tag th_interface_side;
    CHKERR m_field.get_moab().tag_get_handle("INTERFACE_SIDE",
                                             th_interface_side);
    int side;
    CHKERR m_field.get_moab().tag_get_data(th_interface_side, &face, 1, &side);
    if (side == 1) {
      elementOrientation = -1;
    }
  }
  MoFEMFunctionReturn(0);
}

getElementOrientation() [3/3]

MoFEMErrorCode FractureMechanics::CrackPropagation::FaceOrientation::getElementOrientation ( MoFEM::Interface &  m_field, const FEMethod *  fe_method_ptr  )

inlinevirtual

Reimplemented from SurfaceSlidingConstrains::DriverElementOrientation.

Definition at line 1016 of file CrackPropagation.hpp.

                                                                        {
      MoFEMFunctionBegin;
      EntityHandle face = fe_method_ptr->numeredEntFiniteElementPtr->getEnt();
      BitRefLevel bit = fe_method_ptr->problemPtr->getBitRefLevel();
      CHKERR getElementOrientation(m_field, face, bit);
      MoFEMFunctionReturn(0);
    }

Member Data Documentation

contactFaces

Range FractureMechanics::CrackPropagation::FaceOrientation::contactFaces

Definition at line 1006 of file CrackPropagation.hpp.

defaultBitLevel

BitRefLevel FractureMechanics::CrackPropagation::FaceOrientation::defaultBitLevel

Definition at line 1007 of file CrackPropagation.hpp.

useProjectionFromCrackFront

bool FractureMechanics::CrackPropagation::FaceOrientation::useProjectionFromCrackFront

Definition at line 1005 of file CrackPropagation.hpp.

---

The documentation for this struct was generated from the following files:

• users_modules/fracture_mechanics/src/CrackPropagation.hpp
• users_modules/fracture_mechanics/src/impl/CrackPropagation.cpp