12 #include "libmesh/replicated_mesh.h" 23 return std::get<0>(s) < i;
27 return i < std::get<0>(s);
39 const std::vector<libMesh::BoundaryInfo::BCTuple> & bdry_side_list,
41 std::vector<dof_id_type> & converted_elems_ids);
51 const std::vector<libMesh::BoundaryInfo::BCTuple> & bdry_side_list,
53 std::vector<dof_id_type> & converted_elems_ids);
63 const std::vector<libMesh::BoundaryInfo::BCTuple> & bdry_side_list,
65 std::vector<dof_id_type> & converted_elems_ids);
77 const unsigned int sec_min_pos,
78 std::vector<unsigned int> & face_rotation,
79 std::vector<unsigned int> & node_rotation);
100 std::vector<std::vector<unsigned int>> & rotated_tet_face_indices,
101 std::vector<std::vector<const Node *>> & tet_nodes_list);
129 std::vector<std::vector<unsigned int>>
131 std::vector<std::vector<unsigned int>> & tet_face_indices);
140 std::vector<unsigned int> & face_rotation,
141 std::vector<unsigned int> & node_rotation);
155 std::vector<std::vector<unsigned int>> & rotated_tet_face_indices,
156 std::vector<std::vector<const Node *>> & tet_nodes_list);
167 std::vector<std::vector<unsigned int>>
169 std::vector<std::vector<unsigned int>> & tet_face_indices);
180 std::vector<unsigned int> & face_rotation,
181 std::vector<unsigned int> & node_rotation);
195 std::vector<std::vector<unsigned int>> & rotated_tet_face_indices,
196 std::vector<std::vector<const Node *>> & tet_nodes_list);
211 const std::vector<std::pair<dof_id_type, bool>> & elems_to_process,
212 std::vector<dof_id_type> & converted_elems_ids_to_track,
214 const bool delete_block_to_remove);
232 const unsigned short n_elem_sides,
233 std::vector<std::vector<boundary_id_type>> & elem_side_list);
246 const std::vector<unsigned int> & side_indices,
247 const std::vector<std::vector<boundary_id_type>> & elem_side_info,
260 const std::vector<unsigned int> & side_indices,
261 const std::vector<std::vector<boundary_id_type>> & elem_side_info,
275 const unsigned int & side_index,
276 const Node * new_node,
277 const std::vector<boundary_id_type> & side_info,
291 const unsigned int & side_index,
292 const Node * new_node,
293 const std::vector<boundary_id_type> & side_info,
306 const std::vector<unsigned int> & side_indices,
307 const std::vector<std::vector<boundary_id_type>> & elem_side_info,
320 const unsigned int & side_index,
321 const Node * new_node,
322 const std::vector<boundary_id_type> & side_info,
337 const unsigned int & side_index,
338 const Node * new_node,
339 const std::vector<boundary_id_type> & side_info,
351 const std::vector<std::vector<boundary_id_type>> & elem_side_info,
373 const std::vector<BoundaryName> & boundary_names,
374 const unsigned int conversion_element_layer_number,
375 const bool external_boundaries_checking);
388 ReplicatedMesh & mesh,
389 const std::set<subdomain_id_type> & original_subdomain_ids,
391 const SubdomainName & tet_suffix,
392 const SubdomainName & pyramid_suffix);
std::tuple< dof_id_type, unsigned short int, boundary_id_type > BCTuple
void convertPyramid5Elem(ReplicatedMesh &mesh, const dof_id_type &elem_id, const std::vector< std::vector< boundary_id_type >> &elem_side_info, const SubdomainID &subdomain_id_shift_base)
Convert a PYRAMID5 element to elements with TRI3 surfaces on the original QUAD4 side.
void pyramidElemSplitter(ReplicatedMesh &mesh, const std::vector< libMesh::BoundaryInfo::BCTuple > &bdry_side_list, const dof_id_type elem_id, std::vector< dof_id_type > &converted_elems_ids)
Split a PYRAMID5 element into two TET4 elements.
std::vector< bool > quadFaceDiagonalDirectionsHex(const std::vector< const Node *> &hex_nodes)
For a HEX8 element, determine the direction of the diagonal line of each face that involves the node ...
void nodeRotationHEX8(const unsigned int min_id_index, const unsigned int sec_min_pos, std::vector< unsigned int > &face_rotation, std::vector< unsigned int > &node_rotation)
Rotate a HEX8 element's nodes to ensure that the node with the minimum id is the first node; and the ...
void createUnitTet4FromHex8(ReplicatedMesh &mesh, const dof_id_type &elem_id, const unsigned int &side_index, const Node *new_node, const std::vector< boundary_id_type > &side_info, const SubdomainID &subdomain_id_shift_base)
Create two TET4 elements based on a side and the centroid of the HEX8 element.
void createUnitPyramid5FromPrism6(ReplicatedMesh &mesh, const dof_id_type &elem_id, const unsigned int &side_index, const Node *new_node, const std::vector< boundary_id_type > &side_info, const SubdomainID &subdomain_id_shift_base)
Create a PYRAMID5 element opposite the sides converted to tets and the centroid of the PRISM6 element...
bool quadFaceDiagonalDirection(const std::vector< const Node *> &quad_nodes)
For a QUAD4 element, determine the direction of the diagonal line that involves the node with the min...
void hexNodesToTetNodesDeterminer(std::vector< const Node *> &hex_nodes, std::vector< std::vector< unsigned int >> &rotated_tet_face_indices, std::vector< std::vector< const Node *>> &tet_nodes_list)
For a vector of rotated nodes that can form a HEX8 element, create a vector of four-node sets that ca...
void elementBoundaryInfoCollector(const std::vector< libMesh::BoundaryInfo::BCTuple > &bdry_side_list, const dof_id_type elem_id, const unsigned short n_elem_sides, std::vector< std::vector< boundary_id_type >> &elem_side_list)
Collect the boundary information of the given element in a mesh.
std::vector< unsigned int > neighborNodeIndicesHEX8(unsigned int min_id_index)
Calculate the indices (within the element nodes) of the three neighboring nodes of a node in a HEX8 e...
void retainEEID(ReplicatedMesh &mesh, const dof_id_type &elem_id, Elem *new_elem_ptr)
Retain the extra integer of the original element in a new element.
std::vector< std::vector< unsigned int > > tetNodesForPrism(const bool diagonal_direction, std::vector< std::vector< unsigned int >> &tet_face_indices)
Creates sets of four nodes indices that can form TET4 elements to replace the original PRISM6 element...
void nodeRotationPRISM6(unsigned int min_id_index, std::vector< unsigned int > &face_rotation, std::vector< unsigned int > &node_rotation)
Rotate a PRISM6 element nodes to ensure that the node with the minimum id is the first node...
void createUnitTet4FromPrism6(ReplicatedMesh &mesh, const dof_id_type &elem_id, const unsigned int &side_index, const Node *new_node, const std::vector< boundary_id_type > &side_info, const SubdomainID &subdomain_id_shift_base)
Create one or two TET4 elements based on a side and the centroid of the PRISM6 element.
void nodeRotationPYRAMID5(unsigned int min_id_index, std::vector< unsigned int > &face_rotation, std::vector< unsigned int > &node_rotation)
Rotate a PYRAMID5 element nodes to ensure that the node with the minimum id is the first node for the...
void convert3DMeshToAllTet4(ReplicatedMesh &mesh, const std::vector< std::pair< dof_id_type, bool >> &elems_to_process, std::vector< dof_id_type > &converted_elems_ids_to_track, const subdomain_id_type block_id_to_remove, const bool delete_block_to_remove)
Convert all the elements in a 3D mesh, consisting of only linear elements, into TET4 elements...
bool operator()(const libMesh::BoundaryInfo::BCTuple &s, dof_id_type i) const
void convertElem(ReplicatedMesh &mesh, const dof_id_type &elem_id, const std::vector< unsigned int > &side_indices, const std::vector< std::vector< boundary_id_type >> &elem_side_info, const SubdomainID &subdomain_id_shift_base)
Convert the element to an element with TRI3 side-elements on the user-specified sides by modifying th...
void prismElemSplitter(ReplicatedMesh &mesh, const std::vector< libMesh::BoundaryInfo::BCTuple > &bdry_side_list, const dof_id_type elem_id, std::vector< dof_id_type > &converted_elems_ids)
Split a PRISM6 element into three TET4 elements.
void assignConvertedElementsSubdomainNameSuffix(ReplicatedMesh &mesh, const std::set< subdomain_id_type > &original_subdomain_ids, const subdomain_id_type sid_shift_base, const SubdomainName &tet_suffix, const SubdomainName &pyramid_suffix)
Assign a subdomain name suffix to the converted elements created during transition layer generation...
bool operator()(dof_id_type i, const libMesh::BoundaryInfo::BCTuple &s) const
void pyramidNodesToTetNodesDeterminer(std::vector< const Node *> &pyramid_nodes, std::vector< std::vector< unsigned int >> &rotated_tet_face_indices, std::vector< std::vector< const Node *>> &tet_nodes_list)
For a rotated nodes that can form a PYRAMID5 element, create a series of four-node set that can form ...
void hexElemSplitter(ReplicatedMesh &mesh, const std::vector< libMesh::BoundaryInfo::BCTuple > &bdry_side_list, const dof_id_type elem_id, std::vector< dof_id_type > &converted_elems_ids)
Split a HEX8 element into six TET4 elements.
void convertPrism6Elem(ReplicatedMesh &mesh, const dof_id_type &elem_id, const std::vector< unsigned int > &side_indices, const std::vector< std::vector< boundary_id_type >> &elem_side_info, const SubdomainID &subdomain_id_shift_base)
Convert a PRISM6 element to elements with TRI3 surfaces on the given original QUAD4 side(s)...
void convertHex8Elem(ReplicatedMesh &mesh, const dof_id_type &elem_id, const std::vector< unsigned int > &side_indices, const std::vector< std::vector< boundary_id_type >> &elem_side_info, const SubdomainID &subdomain_id_shift_base)
Convert a HEX8 element to elements with TRI3 surfaces on the given original QUAD4 side(s)...
std::vector< std::vector< unsigned int > > tetNodesForHex(const std::vector< bool > diagonal_directions, std::vector< std::vector< unsigned int >> &tet_face_indices)
Creates sets of four nodes indices that can form TET4 elements to replace the original HEX8 element...
void prismNodesToTetNodesDeterminer(std::vector< const Node *> &prism_nodes, std::vector< std::vector< unsigned int >> &rotated_tet_face_indices, std::vector< std::vector< const Node *>> &tet_nodes_list)
For a rotated nodes that can form a PRISM6 element, create a series of four-node set that can form TE...
void createUnitPyramid5FromHex8(ReplicatedMesh &mesh, const dof_id_type &elem_id, const unsigned int &side_index, const Node *new_node, const std::vector< boundary_id_type > &side_info, const SubdomainID &subdomain_id_shift_base)
Create one PYRAMID5 element based on a side and the centroid of the HEX8 element. ...
void transitionLayerGenerator(ReplicatedMesh &mesh, const std::vector< BoundaryName > &boundary_names, const unsigned int conversion_element_layer_number, const bool external_boundaries_checking)
Generate a transition layer of elements with TRI3 surfaces on the given boundaries.