def _create_main_search(self, key = "0"):
""" This method creates the search process that will be use during contact search
Keyword arguments:
self -- It signifies an instance of a class.
key -- The key to identify the current pair
"""
search_parameters = KM.Parameters("""{"condition_name": "", "final_string": "", "predefined_master_slave" : true, "id_name" : ""}""")
search_parameters.AddValue("type_search", self.settings["search_parameters"]["type_search"])
search_parameters.AddValue("check_gap", self.settings["search_parameters"]["check_gap"])
search_parameters.AddValue("allocation_size", self.settings["search_parameters"]["max_number_results"])
search_parameters.AddValue("bucket_size", self.settings["search_parameters"]["bucket_size"])
search_parameters.AddValue("search_factor", self.settings["search_parameters"]["search_factor"])
search_parameters.AddValue("dynamic_search", self.settings["search_parameters"]["dynamic_search"])
search_parameters.AddValue("static_check_movement", self.settings["search_parameters"]["static_check_movement"])
search_parameters.AddValue("consider_gap_threshold", self.settings["search_parameters"]["consider_gap_threshold"])
search_parameters.AddValue("debug_mode", self.settings["search_parameters"]["debug_mode"])
search_parameters["condition_name"].SetString(self._get_condition_name())
search_parameters["final_string"].SetString(self._get_final_string())
self.__assume_master_slave(key)
search_parameters["predefined_master_slave"].SetBool(self.predefined_master_slave)
search_parameters["id_name"].SetString(key)
# We compute the number of nodes of the geometry
number_nodes, number_nodes_master = self._compute_number_nodes()
# We create the search process
simple_search = self.settings["search_parameters"]["simple_search"].GetBool()
if self.dimension == 2:
if simple_search:
self.search_search[key] = CSMA.SimpleContactSearch2D2N(self.computing_model_part, search_parameters)
else:
self.search_search[key] = CSMA.AdvancedContactSearch2D2N(self.computing_model_part, search_parameters)
else:
if number_nodes == 3:
if number_nodes_master == 3:
if simple_search:
self.search_search[key] = CSMA.SimpleContactSearch3D3N(self.computing_model_part, search_parameters)
else:
self.search_search[key] = CSMA.AdvancedContactSearch3D3N(self.computing_model_part, search_parameters)
else:
if simple_search:
self.search_search[key] = CSMA.SimpleContactSearch3D3N4N(self.computing_model_part, search_parameters)
else:
self.search_search[key] = CSMA.AdvancedContactSearch3D3N4N(self.computing_model_part, search_parameters)
elif number_nodes == 4:
if number_nodes_master == 3:
if simple_search:
self.search_search[key] = CSMA.SimpleContactSearch3D4N3N(self.computing_model_part, search_parameters)
else:
self.search_search[key] = CSMA.AdvancedContactSearch3D4N3N(self.computing_model_part, search_parameters)
else:
if simple_search:
self.search_search[key] = CSMA.SimpleContactSearch3D4N(self.computing_model_part, search_parameters)
else:
self.search_search[key] = CSMA.AdvancedContactSearch3D4N(self.computing_model_part, search_parameters)
else:
raise Exception("Geometries not compatible. Check all the geometries are linear")
def __base_test_mapping(self, input_filename, num_nodes, pure_implicit):
KratosMultiphysics.Logger.GetDefaultOutput().SetSeverity(KratosMultiphysics.Logger.Severity.WARNING)
self.model = KratosMultiphysics.Model()
self.main_model_part = self.model.CreateModelPart("Structure", 2)
## Creation of the Kratos model (build sub_model_parts or submeshes)
self.StructureModel = {"Structure": self.main_model_part}
self.main_model_part.AddNodalSolutionStepVariable(KratosMultiphysics.DISPLACEMENT)
self.main_model_part.AddNodalSolutionStepVariable(KratosMultiphysics.VELOCITY)
self.main_model_part.AddNodalSolutionStepVariable(KratosMultiphysics.ACCELERATION)
self.main_model_part.AddNodalSolutionStepVariable(KratosMultiphysics.TEMPERATURE)
self.main_model_part.AddNodalSolutionStepVariable(KratosMultiphysics.NORMAL)
self.main_model_part.AddNodalSolutionStepVariable(ContactStructuralMechanicsApplication.LAGRANGE_MULTIPLIER_CONTACT_PRESSURE)
self.main_model_part.AddNodalSolutionStepVariable(ContactStructuralMechanicsApplication.WEIGHTED_GAP)
self.main_model_part.AddNodalSolutionStepVariable(KratosMultiphysics.NODAL_H)
self.main_model_part.CloneTimeStep(1.01)
KratosMultiphysics.ModelPartIO(input_filename).ReadModelPart(self.main_model_part)
KratosMultiphysics.VariableUtils().AddDof(KratosMultiphysics.DISPLACEMENT_X,self.main_model_part)
KratosMultiphysics.VariableUtils().AddDof(KratosMultiphysics.DISPLACEMENT_Y,self.main_model_part)
KratosMultiphysics.VariableUtils().AddDof(KratosMultiphysics.DISPLACEMENT_Z,self.main_model_part)
KratosMultiphysics.VariableUtils().AddDof(KratosMultiphysics.TEMPERATURE, self.main_model_part)
if (self.main_model_part.HasSubModelPart("Contact")):
interface_model_part = self.main_model_part.GetSubModelPart("Contact")
else:
interface_model_part = self.main_model_part.CreateSubModelPart("Contact")
self.mapping_model_part = self.main_model_part.GetSubModelPart("DISPLACEMENT_Displacement_Auto2")
self.model_part_slave = self.main_model_part.GetSubModelPart("Parts_Parts_Auto1")
KratosMultiphysics.VariableUtils().SetFlag(KratosMultiphysics.SLAVE, True, self.model_part_slave.Nodes)
KratosMultiphysics.VariableUtils().SetFlag(KratosMultiphysics.MASTER, False, self.model_part_slave.Nodes)
self.model_part_master = self.main_model_part.GetSubModelPart("Parts_Parts_Auto2")
KratosMultiphysics.VariableUtils().SetFlag(KratosMultiphysics.SLAVE, False, self.model_part_master.Nodes)
KratosMultiphysics.VariableUtils().SetFlag(KratosMultiphysics.MASTER, True, self.model_part_master.Nodes)
for prop in self.main_model_part.GetProperties():
prop[ContactStructuralMechanicsApplication.INTEGRATION_ORDER_CONTACT] = 3
self.main_model_part.ProcessInfo[ContactStructuralMechanicsApplication.ACTIVE_CHECK_FACTOR] = 3.0e-1
for node in self.mapping_model_part.Nodes:
node.Set(KratosMultiphysics.INTERFACE, True)
Preprocess = ContactStructuralMechanicsApplication.InterfacePreprocessCondition(self.main_model_part)
interface_parameters = KratosMultiphysics.Parameters("""{"simplify_geometry": false}""")
Preprocess.GenerateInterfacePart3D(self.mapping_model_part, interface_parameters)
# We compute NODAL_H that can be used in the search and some values computation
find_nodal_h = KratosMultiphysics.FindNodalHProcess(self.mapping_model_part)
find_nodal_h.Execute()
# We copy the conditions to the ContactSubModelPart
for cond in self.mapping_model_part.Conditions:
interface_model_part.AddCondition(cond)
del(cond)
for node in self.mapping_model_part.Nodes:
interface_model_part.AddNode(node, 0)
del(node)
# We initialize the conditions
alm_init_var = ContactStructuralMechanicsApplication.ALMFastInit(self.mapping_model_part)
alm_init_var.Execute()
search_parameters = KratosMultiphysics.Parameters("""
{
"search_factor" : 3.5,
"allocation_size" : 1000,
"check_gap" : "NoCheck",
"type_search" : "InRadius"
}
""")
if (num_nodes == 3):
contact_search = ContactStructuralMechanicsApplication.AdvancedContactSearch3D3N(self.main_model_part, search_parameters)
else:
contact_search = ContactStructuralMechanicsApplication.AdvancedContactSearch3D4N(self.main_model_part, search_parameters)
# We initialize the search utility
contact_search.CreatePointListMortar()
contact_search.InitializeMortarConditions()
contact_search.UpdateMortarConditions()
for node in self.model_part_master.Nodes:
x = node.X
y = node.Y
z = node.Z
node.SetSolutionStepValue(KratosMultiphysics.TEMPERATURE, z)
node.SetSolutionStepValue(KratosMultiphysics.DISPLACEMENT_X, x)
node.SetSolutionStepValue(KratosMultiphysics.DISPLACEMENT_Y, y)
node.SetSolutionStepValue(KratosMultiphysics.DISPLACEMENT_Z, z)
del(node)
map_parameters = KratosMultiphysics.Parameters("""
{
"echo_level" : 0,
"absolute_convergence_tolerance" : 1.0e-9,
"relative_convergence_tolerance" : 1.0e-4,
"max_number_iterations" : 10,
"integration_order" : 2
}
""")
# Copy respective conditions
for cond in self.main_model_part.Conditions:
if (cond.Is(KratosMultiphysics.SLAVE)):
self.model_part_slave.AddCondition(cond)
if (cond.Is(KratosMultiphysics.MASTER)):
self.model_part_master.AddCondition(cond)
if (pure_implicit == True):
#linear_solver = ExternalSolversApplication.SuperLUSolver()
linear_solver = KratosMultiphysics.SkylineLUFactorizationSolver()
if (num_nodes == 3):
self.mortar_mapping_double = KratosMultiphysics.SimpleMortarMapperProcess3D3NDouble(self.model_part_master, self.model_part_slave, KratosMultiphysics.TEMPERATURE, map_parameters, linear_solver)
self.mortar_mapping_vector = KratosMultiphysics.SimpleMortarMapperProcess3D3NVector(self.model_part_master, self.model_part_slave, KratosMultiphysics.DISPLACEMENT, map_parameters, linear_solver)
else:
self.mortar_mapping_double = KratosMultiphysics.SimpleMortarMapperProcess3D4NDouble(self.model_part_master, self.model_part_slave, KratosMultiphysics.TEMPERATURE, map_parameters, linear_solver)
self.mortar_mapping_vector = KratosMultiphysics.SimpleMortarMapperProcess3D4NVector(self.model_part_master, self.model_part_slave, KratosMultiphysics.DISPLACEMENT, map_parameters, linear_solver)
else:
if (num_nodes == 3):
self.mortar_mapping_double = KratosMultiphysics.SimpleMortarMapperProcess3D3NDouble(self.model_part_master, self.model_part_slave, KratosMultiphysics.TEMPERATURE, map_parameters)
self.mortar_mapping_vector = KratosMultiphysics.SimpleMortarMapperProcess3D3NVector(self.model_part_master, self.model_part_slave, KratosMultiphysics.DISPLACEMENT, map_parameters)
else:
self.mortar_mapping_double = KratosMultiphysics.SimpleMortarMapperProcess3D4NDouble(self.model_part_master, self.model_part_slave, KratosMultiphysics.TEMPERATURE, map_parameters)
self.mortar_mapping_vector = KratosMultiphysics.SimpleMortarMapperProcess3D4NVector(self.model_part_master, self.model_part_slave, KratosMultiphysics.DISPLACEMENT, map_parameters)
def _dynamic_search_tests(self, input_filename, num_nodes):
KM.Logger.GetDefaultOutput().SetSeverity(KM.Logger.Severity.WARNING)
self.model = KM.Model()
self.main_model_part = self.model.CreateModelPart("Structure", 2)
## Creation of the Kratos model (build sub_model_parts or submeshes)
self.StructureModel = {"Structure": self.main_model_part}
self.main_model_part.AddNodalSolutionStepVariable(KM.DISPLACEMENT)
self.main_model_part.AddNodalSolutionStepVariable(KM.VELOCITY)
self.main_model_part.AddNodalSolutionStepVariable(KM.ACCELERATION)
self.main_model_part.AddNodalSolutionStepVariable(
KM.VOLUME_ACCELERATION)
self.main_model_part.AddNodalSolutionStepVariable(KM.REACTION)
self.main_model_part.AddNodalSolutionStepVariable(KM.NORMAL)
self.main_model_part.AddNodalSolutionStepVariable(
CSMA.LAGRANGE_MULTIPLIER_CONTACT_PRESSURE)
self.main_model_part.AddNodalSolutionStepVariable(CSMA.WEIGHTED_GAP)
self.main_model_part.AddNodalSolutionStepVariable(KM.NODAL_H)
self.main_model_part.CloneTimeStep(1.01)
KM.ModelPartIO(input_filename).ReadModelPart(self.main_model_part)
KM.VariableUtils().AddDof(KM.DISPLACEMENT_X, KM.REACTION_X,
self.main_model_part)
KM.VariableUtils().AddDof(KM.DISPLACEMENT_Y, KM.REACTION_Y,
self.main_model_part)
KM.VariableUtils().AddDof(KM.DISPLACEMENT_Z, KM.REACTION_Z,
self.main_model_part)
KM.VariableUtils().AddDof(CSMA.LAGRANGE_MULTIPLIER_CONTACT_PRESSURE,
CSMA.WEIGHTED_GAP, self.main_model_part)
if (self.main_model_part.HasSubModelPart("Contact")):
interface_model_part = self.main_model_part.GetSubModelPart(
"Contact")
else:
interface_model_part = self.main_model_part.CreateSubModelPart(
"Contact")
self.contact_model_part = self.main_model_part.GetSubModelPart(
"DISPLACEMENT_Displacement_Auto2")
model_part_slave = self.main_model_part.GetSubModelPart(
"Parts_Parts_Auto1")
model_part_master = self.main_model_part.GetSubModelPart(
"Parts_Parts_Auto2")
KM.VariableUtils().SetFlag(KM.SLAVE, False,
self.contact_model_part.Nodes)
KM.VariableUtils().SetFlag(KM.MASTER, True,
self.contact_model_part.Nodes)
KM.VariableUtils().SetFlag(KM.SLAVE, True, model_part_slave.Nodes)
KM.VariableUtils().SetFlag(KM.MASTER, False, model_part_slave.Nodes)
for node in model_part_slave.Nodes:
# DEBUG
#node.X -= 9.81 / 32.0
#node.SetSolutionStepValue(KM.DISPLACEMENT_X, -9.81 / 32.0)
node.SetSolutionStepValue(KM.ACCELERATION_X, 1, -9.81)
del (node)
self.main_model_part.ProcessInfo[KM.STEP] = 1
self.main_model_part.ProcessInfo[KM.DELTA_TIME] = 0.5
for prop in self.main_model_part.GetProperties():
prop[CSMA.INTEGRATION_ORDER_CONTACT] = 3
self.main_model_part.ProcessInfo[CSMA.ACTIVE_CHECK_FACTOR] = 3.0e-1
KM.VariableUtils().SetFlag(KM.INTERFACE, True,
self.contact_model_part.Nodes)
pre_process = CSMA.InterfacePreprocessCondition(self.main_model_part)
interface_parameters = KM.Parameters(
"""{"simplify_geometry": false}""")
pre_process.GenerateInterfacePart3D(self.contact_model_part,
interface_parameters)
# We copy the conditions to the ContactSubModelPart
for cond in self.contact_model_part.Conditions:
interface_model_part.AddCondition(cond)
del (cond)
for node in self.contact_model_part.Nodes:
interface_model_part.AddNode(node, 0)
del (node)
# We compute NODAL_H that can be used in the search and some values computation
self.find_nodal_h = KM.FindNodalHProcess(self.contact_model_part)
self.find_nodal_h.Execute()
# We initialize the conditions
alm_init_var = CSMA.ALMFastInit(self.contact_model_part)
alm_init_var.Execute()
search_parameters = KM.Parameters("""
{
"dynamic_search" : true
}
""")
if (num_nodes == 3):
contact_search = CSMA.AdvancedContactSearch3D3N(
self.main_model_part, search_parameters)
else:
contact_search = CSMA.AdvancedContactSearch3D4N(
self.main_model_part, search_parameters)
# We initialize the search utility
contact_search.CreatePointListMortar()
contact_search.InitializeMortarConditions()
contact_search.UpdateMortarConditions()
## DEBUG
#self.__post_process()
import from_json_check_result_process
check_parameters = KM.Parameters("""
{
"check_variables" : ["NORMAL_GAP"],
"input_file_name" : "",
"model_part_name" : "Structure",
"historical_value" : false,
"time_frequency" : 0.0,
"sub_model_part_name" : "Parts_Parts_Auto1"
}
""")
check_parameters["input_file_name"].SetString(input_filename +
"_dynamic_search.json")
check = from_json_check_result_process.FromJsonCheckResultProcess(
self.StructureModel, check_parameters)
check.ExecuteInitialize()
check.ExecuteBeforeSolutionLoop()
check.ExecuteFinalizeSolutionStep()
def __base_test_integration(self, input_filename, num_nodes):
KratosMultiphysics.Logger.GetDefaultOutput().SetSeverity(
KratosMultiphysics.Logger.Severity.WARNING)
self.model = KratosMultiphysics.Model()
self.main_model_part = self.model.CreateModelPart("Structure", 2)
self.main_model_part.AddNodalSolutionStepVariable(
KratosMultiphysics.DISPLACEMENT)
self.main_model_part.AddNodalSolutionStepVariable(
KratosMultiphysics.VELOCITY)
self.main_model_part.AddNodalSolutionStepVariable(
KratosMultiphysics.ACCELERATION)
self.main_model_part.AddNodalSolutionStepVariable(
KratosMultiphysics.REACTION)
self.main_model_part.AddNodalSolutionStepVariable(
KratosMultiphysics.NORMAL)
self.main_model_part.AddNodalSolutionStepVariable(
ContactStructuralMechanicsApplication.
LAGRANGE_MULTIPLIER_CONTACT_PRESSURE)
self.main_model_part.AddNodalSolutionStepVariable(
ContactStructuralMechanicsApplication.WEIGHTED_GAP)
self.main_model_part.AddNodalSolutionStepVariable(
KratosMultiphysics.NODAL_H)
KratosMultiphysics.ModelPartIO(input_filename).ReadModelPart(
self.main_model_part)
KratosMultiphysics.VariableUtils().AddDof(
KratosMultiphysics.DISPLACEMENT_X, KratosMultiphysics.REACTION_X,
self.main_model_part)
KratosMultiphysics.VariableUtils().AddDof(
KratosMultiphysics.DISPLACEMENT_Y, KratosMultiphysics.REACTION_Y,
self.main_model_part)
KratosMultiphysics.VariableUtils().AddDof(
KratosMultiphysics.DISPLACEMENT_Z, KratosMultiphysics.REACTION_Z,
self.main_model_part)
KratosMultiphysics.VariableUtils().AddDof(
ContactStructuralMechanicsApplication.
LAGRANGE_MULTIPLIER_CONTACT_PRESSURE,
ContactStructuralMechanicsApplication.WEIGHTED_GAP,
self.main_model_part)
if (self.main_model_part.HasSubModelPart("Contact")):
interface_model_part = self.main_model_part.GetSubModelPart(
"Contact")
else:
interface_model_part = self.main_model_part.CreateSubModelPart(
"Contact")
self.contact_model_part = self.main_model_part.GetSubModelPart(
"DISPLACEMENT_Displacement_Auto2")
for node in self.contact_model_part.Nodes:
node.Set(KratosMultiphysics.SLAVE, False)
del (node)
model_part_slave = self.main_model_part.GetSubModelPart(
"Parts_Parts_Auto1")
for node in model_part_slave.Nodes:
node.Set(KratosMultiphysics.SLAVE, True)
del (node)
for prop in self.main_model_part.GetProperties():
prop[ContactStructuralMechanicsApplication.
INTEGRATION_ORDER_CONTACT] = 3
self.main_model_part.ProcessInfo[
ContactStructuralMechanicsApplication.ACTIVE_CHECK_FACTOR] = 3.0e-1
for node in self.contact_model_part.Nodes:
node.Set(KratosMultiphysics.INTERFACE, True)
Preprocess = ContactStructuralMechanicsApplication.InterfacePreprocessCondition(
self.main_model_part)
interface_parameters = KratosMultiphysics.Parameters(
"""{"simplify_geometry": false}""")
Preprocess.GenerateInterfacePart3D(self.contact_model_part,
interface_parameters)
# We copy the conditions to the ContactSubModelPart
for cond in self.contact_model_part.Conditions:
interface_model_part.AddCondition(cond)
del (cond)
for node in self.contact_model_part.Nodes:
interface_model_part.AddNode(node, 0)
del (node)
# We initialize the conditions
alm_init_var = ContactStructuralMechanicsApplication.ALMFastInit(
self.contact_model_part)
alm_init_var.Execute()
search_parameters = KratosMultiphysics.Parameters("""
{
"search_factor" : 3.5,
"allocation_size" : 1000,
"check_gap" : "NoCheck",
"type_search" : "InRadius"
}
""")
if (num_nodes == 3):
contact_search = ContactStructuralMechanicsApplication.AdvancedContactSearch3D3N(
self.main_model_part, search_parameters)
else:
contact_search = ContactStructuralMechanicsApplication.AdvancedContactSearch3D4N(
self.main_model_part, search_parameters)
# We initialize the search utility
contact_search.CreatePointListMortar()
contact_search.InitializeMortarConditions()
contact_search.UpdateMortarConditions()
if (num_nodes == 3):
## DEBUG
#print(self.main_model_part)
#self.__post_process()
self.exact_integration = KratosMultiphysics.ExactMortarIntegrationUtility3D3N(
3)
else:
## DEBUG
#print(self.main_model_part)
#self.__post_process()
self.exact_integration = KratosMultiphysics.ExactMortarIntegrationUtility3D4N(
3)
