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.DynamicStructuralSolver(model_part, domain_size) solver.structure_linear_solver = SuperLUSolver() solver.CalculateReactionFlag = True CT = Kratos_Structural_Application_var.Convergence_Tolerance AT = Kratos_Structural_Application_var.Absolute_Tolerance if (Kratos_Structural_Application_var.Convergence_Criteria == "Displacement_Criteria"): solver.conv_criteria = DisplacementCriteria(CT, AT) elif (Kratos_Structural_Application_var.Convergence_Criteria == "Residual_Criteria"): solver.conv_criteria = ResidualCriteria(CT, AT) elif (Kratos_Structural_Application_var.Convergence_Criteria == "And_Criteria" ): Displacement = DisplacementCriteria(CT, AT)
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) Dt = fsi_structure_var.Dt nsteps = 10 gid_io.InitializeResults(mesh_name, (model_part).GetMesh()) for step in range(0, nsteps):