def test_InverseMap_non_constant_scalar(self): SetHistoricalNonUniformSolutionScalar( self.interface_model_part_destination.Nodes, KM.TEMPERATURE) self.mapper.InverseMap(KM.PRESSURE, KM.TEMPERATURE) mapper_test_case.CheckHistoricalNonUniformValues( self.interface_model_part_origin, KM.PRESSURE, GetFilePath(self.__GetFileName("inverse_map_scalar")))
def test_InverseMap_non_constant_vector(self): SetHistoricalNonUniformSolutionVector( self.interface_model_part_destination.Nodes, KM.VELOCITY) self.mapper.InverseMap(KM.FORCE, KM.VELOCITY) mapper_test_case.CheckHistoricalNonUniformValues( self.interface_model_part_origin, KM.FORCE, GetFilePath(self.__GetFileName("inverse_map_vector")))
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)