gid_io.ReadModelPart(fluid_model_part)
print(fluid_model_part)
print("fluid model read correctly")
# reading the structural part
data_io = DatafileIO(input_file_name + str("_structure"))
data_io.ReadModelPart(structure_model_part)
print(structure_model_part)
print("structural model read correctly")
# setting up the buffer size: SHOULD BE DONE AFTER READING!!!
fluid_model_part.SetBufferSize(3)
structure_model_part.SetBufferSize(3)
# adding dofs
incompressible_fluid_solver.AddDofs(fluid_model_part)
mesh_solver.AddDofs(fluid_model_part)
structural_solver_dynamic.AddDofs(structure_model_part)
# assigning the fluid properties
density = 1.18
viscosity = 1.82e-5 / density
for node in fluid_model_part.Nodes:
node.SetSolutionStepValue(VISCOSITY, 0, viscosity)
node.SetSolutionStepValue(DENSITY, 0, density)
# creating the solvers
# fluid solver
fluid_solver = incompressible_fluid_solver.IncompressibleFluidSolver(fluid_model_part, domain_size)
def AddDofs(fluid_model_part, structure_model_part):
import incompressible_fluid_solver
incompressible_fluid_solver.AddDofs(fluid_model_part)