def Initialize(self):
self.DEM_parameters["problem_name"].SetString('benchmark' + str(benchmark_number))
#self.end_time = slt.end_time
#self.dt = slt.dt
#self.graph_print_interval = slt.graph_print_interval
super(DEMBenchamarksAnalysisStage, self).Initialize()
Logger.PrintInfo("DEM","Computing points in the curve...", 1 + self.number_of_points_in_the_graphic - self.iteration, "point(s) left to finish....",'\n')
list_of_nodes_ids = [1]
if self.nodeplotter:
os.chdir(self.main_path)
self.plotter = plot_variables.variable_plotter(self.spheres_model_part, list_of_nodes_ids, benchmark_number)
self.tang_plotter = plot_variables.tangential_force_plotter(self.spheres_model_part, list_of_nodes_ids, self.iteration)
"point(s) left to finish....", '\n')
mesh_motion = DEMFEMUtilities()
# creating a Post Utils object that executes several post-related tasks
post_utils = DEM_procedures.PostUtils(DEM_parameters,
spheres_model_part)
list_of_nodes_ids = [1]
if nodeplotter:
os.chdir(main_path)
plotter = plot_variables.variable_plotter(
spheres_model_part, list_of_nodes_ids) #Related to Benchmarks
tang_plotter = plot_variables.tangential_force_plotter(
spheres_model_part, list_of_nodes_ids,
iteration) #Related to Benchmarks
step = 0
while (time < final_time):
#
#SLS dt = spheres_model_part.ProcessInfo.GetValue(DELTA_TIME) # Possible modifications of DELTA_TIME
time = time + dt
step += 1
spheres_model_part.ProcessInfo[TIME] = time
spheres_model_part.ProcessInfo[DELTA_TIME] = dt
spheres_model_part.ProcessInfo[TIME_STEPS] = step
rigid_face_model_part.ProcessInfo[TIME] = time
procedures.ModelData(spheres_model_part, solver)
materialTest.PrintChart()
materialTest.PrepareDataForGraph()
post_utils = DEM_procedures.PostUtils(DEM_parameters, spheres_model_part)
print("Computing points in the curve...", 1 + number_of_points_in_the_graphic - iteration, "point(s) left to finish....",'\n')
list_of_nodes_ids = [1]
if nodeplotter:
os.chdir(main_path)
plotter = plot_variables.variable_plotter(spheres_model_part, list_of_nodes_ids) #Related to Benchmarks
tang_plotter = plot_variables.tangential_force_plotter(spheres_model_part, list_of_nodes_ids, iteration) #Related to Benchmarks
step = 0
##############################################################################
# MAIN LOOP #
##############################################################################
report.total_steps_expected = int(final_time / dt)
KRATOSprint(report.BeginReport(timer))
while (time < final_time):
#SLS dt = spheres_model_part.ProcessInfo.GetValue(DELTA_TIME) # Possible modifications of DELTA_TIME
time = time + dt
step += 1