mesh_name = 0.0 gid_io.InitializeMesh(mesh_name) gid_io.WriteMesh((model_part).GetMesh()) gid_io.FinalizeMesh() print(model_part) print(model_part.Properties) # the buffer size should be set up here after the mesh is read for the first time model_part.SetBufferSize(2) # importing the solver files if (fsi_structure_var.SolverType == "StaticSolver"): print(">>>>>>>>>>>>>>>>< 1yquqee ") structural_solver_static.AddDofs(model_part) solver = structural_solver_static.StaticStructuralSolver( model_part, domain_size) elif (fsi_structure_var.SolverType == "DynamicSolver"): print(" * ASRDF +******************************** FF GDFSG DG G") structural_solver_dynamic.AddDofs(model_part) solver = structural_solver_dynamic.DynamicStructuralSolver( model_part, domain_size) # creating a fluid solver object model_part.Properties[1].SetValue(CONSTITUTIVE_LAW, Isotropic3D()) model_part.Properties[2].SetValue(CONSTITUTIVE_LAW, Isotropic3D()) print("Linear elastic model selected") solver.Initialize() (solver.solver).SetEchoLevel(2)
FindNeighbours() model_part.Properties[1].SetValue(CONSTITUTIVE_LAW, Isotropic2D()) print("Linear elastic model selected") print(model_part) print(model_part.Properties) # the buffer size should be set up here after the mesh is read for the first time model_part.SetBufferSize(3) # importing the rotational dofs degrees of freedom if necessary RotationalDofs() # importing the solver files SolverType.AddDofs(model_part) solver = SolverType.StaticStructuralSolver(model_part, domain_size) tol = 1e-6 max_it = 300 restart = 150 DropTol = 1e-4 FillTol = 1e-2 FillFactor = 30 solver.structure_linear_solver = SuperLUIterativeSolver( tol, max_it, restart, DropTol, FillTol, FillFactor) CT = Kratos_Structural_Application_var.Convergence_Tolerance AT = Kratos_Structural_Application_var.Absolute_Tolerance if (Kratos_Structural_Application_var.Convergence_Criteria ==