def checkResults(self):
with UnitTest.WorkFolderScope(self.work_folder, __file__):
if self.print_reference_values:
json_output_settings = KratosMultiphysics.Parameters(r'''{
"output_variables": ["TEMPERATURE"],
"output_file_name": "",
"model_part_name": "MainModelPart",
"time_frequency": 0.0
}''')
json_output_settings["output_file_name"].SetString(GetFilePath("BFECCConvectionTest/" + self.reference_file + "_results.json"))
json_output_process = JsonOutputProcess(self.model, json_output_settings)
json_output_process.ExecuteInitialize()
json_output_process.ExecuteBeforeSolutionLoop()
json_output_process.ExecuteFinalizeSolutionStep()
else:
json_check_parameters = KratosMultiphysics.Parameters(r'''{
"check_variables" : ["TEMPERATURE"],
"input_file_name" : "",
"model_part_name" : "MainModelPart",
"time_frequency" : 0.0
}''')
json_check_parameters["input_file_name"].SetString(GetFilePath("BFECCConvectionTest/" + self.reference_file + "_results.json"))
json_check_process = FromJsonCheckResultProcess(self.model, json_check_parameters)
json_check_process.ExecuteInitialize()
json_check_process.ExecuteBeforeSolutionLoop()
json_check_process.ExecuteFinalizeSolutionStep()
def RunJsonOutputProcess(model):
parameters = KM.Parameters("""
{
"output_variables" : ["SHAPE_CHANGE_X","SHAPE_CHANGE_Y","SHAPE_CHANGE_Z"],
"output_file_name" : "shape_change_results.json",
"model_part_name" : "structure.Parts_structure",
"time_frequency" : 0.0
}
""")
process = JsonOutputProcess(model, parameters)
process.ExecuteInitialize()
process.ExecuteBeforeSolutionLoop()
process.ExecuteFinalizeSolutionStep()
def _create_reference_solution(current_model):
# The following is used to create the solution database
from KratosMultiphysics.json_output_process import JsonOutputProcess
out_parameters = KratosMultiphysics.Parameters("""
{
"gauss_points_output_variables": ["VON_MISES_STRESS"],
"output_file_name" : "",
"time_frequency" : 0.01,
"model_part_name" : "solid_part",
"sub_model_part_name" : "Body"
}
""")
out_parameters["output_file_name"].SetString(GetFilePath("test_perfect_plasticity_implementation_verification_reference.json"))
out = JsonOutputProcess(current_model, out_parameters)
out.ExecuteInitialize()
out.ExecuteBeforeSolutionLoop()
return out
class TestCookMembrane(KratosUnittest.TestCase):
def setUp(self):
self.print_output = False
self.print_results = False
def test_cook_membrane_2d(self):
results_filename = "cook_membrane_test/cook_membrane_results.json"
parameters_filename = "cook_membrane_test/cook_membrane_parameters.json"
with open(parameters_filename, 'r') as parameter_file:
parameters = KratosMultiphysics.Parameters(parameter_file.read())
model = KratosMultiphysics.Model()
simulation = StructuralMechanicsAnalysis(model, parameters)
simulation.Run()
# self._check_results(model_part, A, b)
if self.print_results:
self.__print_results(model, results_filename)
if self.print_output:
self.__post_process(
model.GetModelPart(parameters["solver_settings"]
["model_part_name"].GetString()))
self.__check_results(model, results_filename)
def test_cook_membrane_incompressible_2d(self):
results_filename = "cook_membrane_test/cook_membrane_incompressible_results.json"
parameters_filename = "cook_membrane_test/cook_membrane_parameters.json"
with open(parameters_filename, 'r') as parameter_file:
parameters = KratosMultiphysics.Parameters(parameter_file.read())
parameters["solver_settings"]["material_import_settings"][
"materials_filename"].SetString(
"cook_membrane_test/cook_membrane_incompressible_materials.json"
)
model = KratosMultiphysics.Model()
simulation = StructuralMechanicsAnalysis(model, parameters)
simulation.Run()
# self._check_results(model_part, A, b)
if self.print_results:
self.__print_results(model, results_filename)
if self.print_output:
self.__post_process(
model.GetModelPart(parameters["solver_settings"]
["model_part_name"].GetString()))
self.__check_results(model, results_filename)
def __print_results(self, model, results_filename):
json_output_settings = KratosMultiphysics.Parameters(r"""
{
"output_variables": ["DISPLACEMENT_X","DISPLACEMENT_Y","VOLUMETRIC_STRAIN"],
"output_file_name": "",
"time_frequency": 0.00,
"model_part_name": "cook_membrane.Parts_ResultsCheck"
}""")
json_output_settings["output_file_name"].SetString(results_filename)
self.json_output = JsonOutputProcess(model, json_output_settings)
self.json_output.ExecuteInitialize()
self.json_output.ExecuteBeforeSolutionLoop()
self.json_output.ExecuteFinalizeSolutionStep()
def __check_results(self, model, results_filename):
from_json_check_result_settings = KratosMultiphysics.Parameters(r"""
{
"check_variables": ["DISPLACEMENT_X","DISPLACEMENT_Y","VOLUMETRIC_STRAIN"],
"input_file_name": "",
"model_part_name": "cook_membrane.Parts_ResultsCheck"
}""")
from_json_check_result_settings["input_file_name"].SetString(
results_filename)
self.from_json_check_result = FromJsonCheckResultProcess(
model, from_json_check_result_settings)
self.from_json_check_result.ExecuteInitialize()
self.from_json_check_result.ExecuteFinalizeSolutionStep()
def __post_process(self, main_model_part, post_type="gid"):
if post_type == "gid":
self.gid_output = GiDOutputProcess(
main_model_part, main_model_part.Name,
KratosMultiphysics.Parameters(r"""
{
"result_file_configuration" : {
"gidpost_flags": {
"GiDPostMode": "GiD_PostBinary",
"WriteDeformedMeshFlag": "WriteUndeformed",
"WriteConditionsFlag": "WriteConditions",
"MultiFileFlag": "SingleFile"
},
"nodal_results" : ["DISPLACEMENT","VOLUMETRIC_STRAIN"],
"gauss_point_results" : []
}
}"""))
self.gid_output.ExecuteInitialize()
self.gid_output.ExecuteBeforeSolutionLoop()
self.gid_output.ExecuteInitializeSolutionStep()
self.gid_output.PrintOutput()
self.gid_output.ExecuteFinalizeSolutionStep()
self.gid_output.ExecuteFinalize()
elif post_type == "vtk":
vtk_output_parameters = KratosMultiphysics.Parameters(r"""
{
"model_part_name": "",
"extrapolate_gauss_points": false,
"nodal_solution_step_data_variables" : ["DISPLACEMENT","VOLUMETRIC_STRAIN"],
"gauss_point_variables": []
}""")
vtk_output_parameters["model_part_name"].SetString(
main_model_part.Name)
self.vtk_output_process = VtkOutputProcess(
main_model_part.GetModel(), vtk_output_parameters)
self.vtk_output_process.ExecuteInitialize()
self.vtk_output_process.ExecuteBeforeSolutionLoop()
self.vtk_output_process.ExecuteInitializeSolutionStep()
self.vtk_output_process.PrintOutput()
self.vtk_output_process.ExecuteFinalizeSolutionStep()
self.vtk_output_process.ExecuteFinalize()