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.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"
):
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())