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