Merge pull request #655 from eulovi/add-wedge-to-embeddedmesh

Add wedge6 to embeddedmesh

Author: eulovi

src/constraint_framework/4C_constraint_framework_embeddedmesh_interaction_pair_mortar.cpp

@@ -692,6 +692,7 @@ namespace CONSTRAINTS::EMBEDDEDMESH
 
   template class SurfaceToBackgroundCouplingPairMortar<t_quad4, t_hex8, t_quad4>;
   template class SurfaceToBackgroundCouplingPairMortar<t_nurbs9, t_hex8, t_nurbs9>;
+  template class SurfaceToBackgroundCouplingPairMortar<t_nurbs9, t_wedge6, t_nurbs9>;
   template class SurfaceToBackgroundCouplingPairMortar<t_quad4, t_nurbs27, t_quad4>;
   template class SurfaceToBackgroundCouplingPairMortar<t_nurbs9, t_nurbs27, t_nurbs9>;
 

src/constraint_framework/4C_constraint_framework_embeddedmesh_solid_to_solid_utils.cpp

@@ -141,6 +141,13 @@ std::shared_ptr<CONSTRAINTS::EMBEDDEDMESH::SolidInteractionPair> coupling_pair_m
               interfaceele_real, background_ele, params_ptr, cutwizard_ptr, boundary_cells);
           break;
         }
+        case Core::FE::CellType::wedge6:
+        {
+          return std::make_shared<CONSTRAINTS::EMBEDDEDMESH::SurfaceToBackgroundCouplingPairMortar<
+              GEOMETRYPAIR::t_nurbs9, GEOMETRYPAIR::t_wedge6, GEOMETRYPAIR::t_nurbs9>>(
+              interfaceele_real, background_ele, params_ptr, cutwizard_ptr, boundary_cells);
+          break;
+        }
         default:
           FOUR_C_THROW(
               "The interaction pairs with background element of type {} not yet implemented",
@@ -513,6 +520,7 @@ namespace CONSTRAINTS::EMBEDDEDMESH
   initialize_template_assemble_local_mortar_contributions(t_nurbs9, t_hex8, t_nurbs9);
   initialize_template_assemble_local_mortar_contributions(t_quad4, t_nurbs27, t_quad4);
   initialize_template_assemble_local_mortar_contributions(t_nurbs9, t_nurbs27, t_nurbs9);
+  initialize_template_assemble_local_mortar_contributions(t_nurbs9, t_wedge6, t_nurbs9);
 
 }  // namespace CONSTRAINTS::EMBEDDEDMESH
 

src/geometry_pair/4C_geometry_pair_element.hpp

@@ -41,7 +41,9 @@ namespace GEOMETRYPAIR
     //! hexahedron
     hexahedron,
     //! tetraeder
-    tetraeder
+    tetraeder,
+    //! wedge
+    wedge
   };
 
   /**
@@ -135,6 +137,12 @@ namespace GEOMETRYPAIR
     static constexpr GEOMETRYPAIR::DiscretizationTypeGeometry geometry_type_ =
         GEOMETRYPAIR::DiscretizationTypeGeometry::hexahedron;
   };
+  template <>
+  struct ElementDiscretizationToGeometryType<Core::FE::CellType::wedge6>
+  {
+    static constexpr GEOMETRYPAIR::DiscretizationTypeGeometry geometry_type_ =
+        GEOMETRYPAIR::DiscretizationTypeGeometry::wedge;
+  };
 
   /**
    * \brief Base class for the geometry pair element type.
@@ -200,6 +208,7 @@ namespace GEOMETRYPAIR
   using t_tet4 = ElementDiscretizationBase<Core::FE::CellType::tet4, 1>;
   using t_tet10 = ElementDiscretizationBase<Core::FE::CellType::tet10, 1>;
   using t_nurbs27 = ElementDiscretizationBase<Core::FE::CellType::nurbs27, 1>;
+  using t_wedge6 = ElementDiscretizationBase<Core::FE::CellType::wedge6, 1>;
 
 
   /**
@@ -265,6 +274,12 @@ namespace GEOMETRYPAIR
     static const bool value_ = true;
   };
 
+  template <>
+  struct IsLagrangeElement<t_wedge6>
+  {
+    static const bool value_ = true;
+  };
+
 
   /**
    * \brief Compile time struct to check if an element type is based on NURBS shape functions

src/geometry_pair/4C_geometry_pair_element_evaluation_functions.hpp

@@ -358,6 +358,13 @@ namespace GEOMETRYPAIR
           xi(0) + xi(1) + xi(2) < 1.0 + Constants::projection_xi_eta_tol)
         return true;
     }
+    else if (ElementType::geometry_type_ == DiscretizationTypeGeometry::wedge)
+    {
+      if (xi(0) > -Constants::projection_xi_eta_tol && xi(0) < xi_limit &&
+          xi(1) > -Constants::projection_xi_eta_tol && xi(1) < xi_limit - xi(0) &&
+          fabs(xi(2)) < xi_limit)
+        return true;
+    }
     else
     {
       FOUR_C_THROW("Wrong DiscretizationTypeGeometry given!");

src/geometry_pair/4C_geometry_pair_factory.cpp

@@ -84,6 +84,10 @@ template std::shared_ptr<GEOMETRYPAIR::GeometryPair>
 GEOMETRYPAIR::geometry_pair_line_to_volume_factory<double, GEOMETRYPAIR::t_hermite,
     GEOMETRYPAIR::t_nurbs27>(const Core::Elements::Element*, const Core::Elements::Element*,
     const std::shared_ptr<GeometryEvaluationDataBase>&);
+template std::shared_ptr<GEOMETRYPAIR::GeometryPair>
+GEOMETRYPAIR::geometry_pair_line_to_volume_factory<double, GEOMETRYPAIR::t_hermite,
+    GEOMETRYPAIR::t_wedge6>(const Core::Elements::Element*, const Core::Elements::Element*,
+    const std::shared_ptr<GeometryEvaluationDataBase>&);
 
 
 /**

src/geometry_pair/4C_geometry_pair_line_projection.cpp

@@ -548,6 +548,7 @@ namespace GEOMETRYPAIR
   initialize_template_volume_gauss_point(double, t_hermite, t_tet4);
   initialize_template_volume_gauss_point(double, t_hermite, t_tet10);
   initialize_template_volume_gauss_point(double, t_hermite, t_nurbs27);
+  initialize_template_volume_gauss_point(double, t_hermite, t_wedge6);
 
   // Define line-to-volume segmentation pairs.
 #define initialize_template_volume_segmentation(a, b, c) \
@@ -559,6 +560,7 @@ namespace GEOMETRYPAIR
   initialize_template_volume_segmentation(double, t_hermite, t_tet4);
   initialize_template_volume_segmentation(double, t_hermite, t_tet10);
   initialize_template_volume_segmentation(double, t_hermite, t_nurbs27);
+  initialize_template_volume_segmentation(double, t_hermite, t_wedge6);
 
   // Helper types for the macro initialization. The compiler has troubles inserting the templated
   // typenames into the macros.

src/geometry_pair/4C_geometry_pair_line_to_volume.cpp

@@ -318,5 +318,7 @@ template class GEOMETRYPAIR::GeometryPairLineToVolume<double, GEOMETRYPAIR::t_he
     GEOMETRYPAIR::t_tet10>;
 template class GEOMETRYPAIR::GeometryPairLineToVolume<double, GEOMETRYPAIR::t_hermite,
     GEOMETRYPAIR::t_nurbs27>;
+template class GEOMETRYPAIR::GeometryPairLineToVolume<double, GEOMETRYPAIR::t_hermite,
+    GEOMETRYPAIR::t_wedge6>;
 
 FOUR_C_NAMESPACE_CLOSE

src/geometry_pair/4C_geometry_pair_line_to_volume_gauss_point_projection.cpp

@@ -89,5 +89,7 @@ template class GEOMETRYPAIR::GeometryPairLineToVolumeGaussPointProjection<double
     GEOMETRYPAIR::t_hermite, GEOMETRYPAIR::t_tet10>;
 template class GEOMETRYPAIR::GeometryPairLineToVolumeGaussPointProjection<double,
     GEOMETRYPAIR::t_hermite, GEOMETRYPAIR::t_nurbs27>;
+template class GEOMETRYPAIR::GeometryPairLineToVolumeGaussPointProjection<double,
+    GEOMETRYPAIR::t_hermite, GEOMETRYPAIR::t_wedge6>;
 
 FOUR_C_NAMESPACE_CLOSE

src/geometry_pair/4C_geometry_pair_line_to_volume_segmentation.cpp

@@ -69,5 +69,7 @@ template class GEOMETRYPAIR::GeometryPairLineToVolumeSegmentation<double, GEOMET
     GEOMETRYPAIR::t_tet10>;
 template class GEOMETRYPAIR::GeometryPairLineToVolumeSegmentation<double, GEOMETRYPAIR::t_hermite,
     GEOMETRYPAIR::t_nurbs27>;
+template class GEOMETRYPAIR::GeometryPairLineToVolumeSegmentation<double, GEOMETRYPAIR::t_hermite,
+    GEOMETRYPAIR::t_wedge6>;
 
 FOUR_C_NAMESPACE_CLOSE