Add contact and mpc conditions

cubitpy/cubitpy_types.py

@@ -241,6 +241,10 @@ def get_default_four_c_description(self):
             return "KINEM nonlinear EAS none ANS none THICKDIR auto"
         elif self == self.hex8_fluid or self.tet4_fluid:
             return "NA ALE"
+        elif self == self.quad4:
+            return (
+                "KINEM nonlinear EAS none THICK 1.0 STRESS_STRAIN plane_stress GP 2 2"
+            )
         else:
             raise ValueError("Got wrong element type {}!".format(self))
 
@@ -255,6 +259,9 @@ class BoundaryConditionType(Enum):
     beam_to_solid_surface_meshtying = auto()
     beam_to_solid_surface_contact = auto()
     solid_to_solid_surface_contact = auto()
+    micro_model_boundary_surface = auto()
+    periodic_rve_edge = auto()
+    line_mortar_contact = auto()
 
     # fluid
     flow_rate = auto()
@@ -322,6 +329,14 @@ def get_dat_bc_section_header(self, geometry_type):
             geometry_type == GeometryType.curve
         ):
             return "FLUID NEUMANN INFLOW LINE CONDITIONS"
+        elif self == self.line_mortar_contact and (geometry_type == GeometryType.curve):
+            return "DESIGN LINE MORTAR CONTACT CONDITIONS 2D"
+        elif self == self.periodic_rve_edge and (geometry_type == GeometryType.curve):
+            return "DESIGN LINE PERIODIC RVE 2D BOUNDARY CONDITIONS"
+        elif self == self.micro_model_boundary_surface and (
+            geometry_type == GeometryType.surface
+        ):
+            return "MICROSCALE CONDITIONS"
         else:
             raise ValueError(
                 "No implemented case for {} and {}!".format(self, geometry_type)