Exemplo n.º 1
0
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)
Exemplo n.º 2
0
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 ==