def test_map_forces_conservative(self):
SetReactions(self.model_part_fluid)
if self.print_output:
mapper_test_case.VtkOutputNodesHistorical(self.model_part_fluid, KM.REACTION, "Blade_empire_mortar_Fluid_prescr_force")
self.mapper_pressure_side.InverseMap(KM.FORCE, KM.REACTION, Mapper.SWAP_SIGN | Mapper.USE_TRANSPOSE)
self.mapper_suction_side.InverseMap(KM.FORCE, KM.REACTION, Mapper.SWAP_SIGN | Mapper.USE_TRANSPOSE)
if self.print_output:
mapper_test_case.VtkOutputNodesHistorical(self.model_part_structure, KM.FORCE, "Blade_empire_mortar_Structure_mapped_force_conserv")
def test_map_displacements(self):
SetDisplacements(self.model_part_structure)
if self.print_output:
mapper_test_case.VtkOutputNodesHistorical(self.model_part_structure, KM.DISPLACEMENT, "Blade_empire_mortar_Structure_prescr_disp")
self.mapper_pressure_side.Map(KM.DISPLACEMENT, KM.MESH_DISPLACEMENT)
self.mapper_suction_side.Map(KM.DISPLACEMENT, KM.MESH_DISPLACEMENT)
if self.print_output:
mapper_test_case.VtkOutputNodesHistorical(self.model_part_fluid, KM.MESH_DISPLACEMENT, "Blade_empire_mortar_Fluid_mapped_disp")
def test_map_displacements(self):
SetDisplacements(self.model_part_structure)
if self.print_output:
mapper_test_case.VtkOutputNodesHistorical(self.model_part_structure, KM.DISPLACEMENT, "Blade_" + self.mapper_type + "_Structure_prescr_disp")
self.mapper_pressure_side.Map(KM.DISPLACEMENT, KM.MESH_DISPLACEMENT)
self.mapper_suction_side.Map(KM.DISPLACEMENT, KM.MESH_DISPLACEMENT)
if self.print_output:
mapper_test_case.VtkOutputNodesHistorical(self.model_part_fluid, KM.MESH_DISPLACEMENT, "Blade_" + self.mapper_type + "_Fluid_mapped_disp")
mapper_test_case.CheckHistoricalNonUniformValues(self.model_part_fluid, KM.MESH_DISPLACEMENT, GetFilePath(self.__GetFileName("balde_map_disp")))
def test_map_forces(self):
SetReactions(self.model_part_fluid)
if self.print_output:
mapper_test_case.VtkOutputNodesHistorical(self.model_part_fluid, KM.REACTION, "Blade_" + self.mapper_type + "_Fluid_prescr_force")
# this would be POINT_LOAD in regular StructuralMechanics (using FORCE to avoid the StructuralMechanics import)
self.mapper_pressure_side.InverseMap(KM.FORCE, KM.REACTION, KratosMapping.Mapper.SWAP_SIGN)
self.mapper_suction_side.InverseMap(KM.FORCE, KM.REACTION, KratosMapping.Mapper.SWAP_SIGN)
if self.print_output:
mapper_test_case.VtkOutputNodesHistorical(self.model_part_structure, KM.FORCE, "Blade_" + self.mapper_type + "_Structure_mapped_force")
mapper_test_case.CheckHistoricalNonUniformValues(self.model_part_structure, KM.FORCE, GetFilePath(self.__GetFileName("balde_map_force")))
def test_map_forces_conservative(self):
SetReactions(self.model_part_fluid)
if self.print_output:
mapper_test_case.VtkOutputNodesHistorical(self.model_part_fluid, KM.REACTION, "Blade_" + self.mapper_type + "_Fluid_prescr_force")
# this would be POINT_LOAD in regular StructuralMechanics (using FORCE to avoid the StructuralMechanics import)
self.mapper_pressure_side.InverseMap(KM.FORCE, KM.REACTION, KratosMapping.Mapper.SWAP_SIGN | KratosMapping.Mapper.USE_TRANSPOSE)
self.__CheckValuesSum(self.model_part_fluid.GetSubModelPart("pressure_side_tri"), self.model_part_structure.GetSubModelPart("pressure_side_quad"), KM.REACTION, KM.FORCE)
# Note: Setting the solution again because in this case some nodes are shared and hence
# would slightly influence the computation
SetReactions(self.model_part_fluid)
self.mapper_suction_side.InverseMap(KM.FORCE, KM.REACTION, KratosMapping.Mapper.SWAP_SIGN | KratosMapping.Mapper.USE_TRANSPOSE)
if self.print_output:
mapper_test_case.VtkOutputNodesHistorical(self.model_part_structure, KM.FORCE, "Blade_" + self.mapper_type + "_Structure_mapped_force_conserv")
mapper_test_case.CheckHistoricalNonUniformValues(self.model_part_structure, KM.FORCE, GetFilePath(self.__GetFileName("balde_map_force_conserv")))
self.__CheckValuesSum(self.model_part_fluid.GetSubModelPart("suction_side_tri"), self.model_part_structure.GetSubModelPart("suction_side_quad"), KM.REACTION, KM.FORCE)