def testSetterGetter(self):
g = np.array([0., 0., -9.81])
system = fsi.ProtChSystem()
system.setGravitationalAcceleration(g)
g0 = system.getGravitationalAcceleration()
npt.assert_almost_equal(g0, g)
mesh = fsi.ProtChMesh(system)
mesh.getChronoObject()
body = fsi.ProtChBody(system)
# position
pos = np.array([1., 2., 3.])
body.setPosition(pos)
pos0 = body.getPosition()
npt.assert_almost_equal(pos0, pos)
# inertia
inertia = np.array([[1., 4., 5.], [4., 2., 6.], [5., 6., 3.]])
inertiaXX = np.array([inertia[0, 0], inertia[1, 1], inertia[2, 2]])
inertiaXY = np.array([inertia[0, 1], inertia[0, 2], inertia[1, 2]])
body.setInertiaXX(inertiaXX)
# inertiaXX0 = body.getInertiaXX()
body.setInertiaXY(inertiaXY)
# inertiaXY0 = body.getInertiaXY()
inertia0 = body.getInertia()
# npt.assert_almost_equal(inertiaXX0, inertiaXX)
# npt.assert_almost_equal(inertiaXY0, inertiaXY)
npt.assert_almost_equal(inertia0, inertia)
# mass
mass = 10.2
body.setMass(mass)
mass0 = body.getMass()
npt.assert_almost_equal(mass0, mass)
def testHangingCableANCF(self):
g = np.array([0., 0., -9.81])
system = fsi.ProtChSystem()
system.ChSystem.Set_G_acc(chrono.ChVectorD(g[0], g[1], g[2]))
system.setTimeStep(1e-1)
# timestepper = chrono.ChTimestepperEulerImplicitLinearized()
# system.ChSystem.SetTimestepper(timestepper)
solver = chrono.ChSolverMINRES()
system.ChSystem.SetSolver(solver)
solver.SetMaxIterations(100)
solver.EnableWarmStart(True)
solver.EnableDiagonalPreconditioner(True)
solver.SetVerbose(True)
system.ChSystem.SetSolverForceTolerance(1e-10)
system.ChSystem.SetSolverMaxIterations(100)
mesh = fsi.ProtChMesh(system)
L = np.array([5.])
nb_elems = np.array([3])
d = np.array([1e-3])
rho = np.array([1000.])
E = np.array([1e10])
cable_type = b"CableANCF"
fairlead_body = fsi.ProtChBody(system)
fairlead_body.ChBody.SetBodyFixed(True)
mooring = fsi.ProtChMoorings(system=system,
mesh=mesh,
length=L,
nb_elems=nb_elems,
d=d,
rho=rho,
E=E,
beam_type=cable_type)
mooring.external_forces_manual = True # tri: should work without this line
# vertical cable
mooring.setNodesPositionFunction(lambda s: np.array([0., 0., s]),
lambda s: np.array([0., 0., 1.]))
mooring.setNodesPosition()
mooring.buildNodes()
mooring.attachBackNodeToBody(fairlead_body)
system.calculate_init()
system.calculate(0.5)
T = mooring.getTensionBack()
strain = mooring.getNodesTension(eta=1.)[-1] * np.pi * d**2 / 4 * E
T_sol = -np.ones(3) * g * rho * (np.pi * d**2 / 4. * L)
npt.assert_almost_equal(-T, T_sol)
# BOUNDARY CONDITIONS tank.BC['x+'].setNoSlip() tank.BC['x-'].setNoSlip() tank.BC['y-'].setNoSlip() tank.BC['y+'].setAtmosphere() rect.BC['x+'].setNoSlip() rect.BC['x-'].setNoSlip() rect.BC['y+'].setNoSlip() rect.BC['y-'].setNoSlip() # CHRONO system = fsi.ProtChSystem() system.ChSystem.Set_G_acc(chrono.ChVectorD(g[0], g[1], 0.)) body = fsi.ProtChBody(system=system) body.attachShape(rect) body.ChBody.SetMass(500.) body.ChBody.SetBodyFixed(True) # fixing body # OTHER PARAMS st.assembleDomain(domain) # ___ _ _ _ _ ____ _ _ _ _ # |_ _|_ __ (_) |_(_) __ _| | / ___|___ _ __ __| (_) |_(_) ___ _ __ ___ # | || '_ \| | __| |/ _` | | | | / _ \| '_ \ / _` | | __| |/ _ \| '_ \/ __| # | || | | | | |_| | (_| | | | |__| (_) | | | | (_| | | |_| | (_) | | | \__ \ # |___|_| |_|_|\__|_|\__,_|_| \____\___/|_| |_|\__,_|_|\__|_|\___/|_| |_|___/ # Initial Conditions
# translate caisson to middle of the tank caisson.translate(np.array([1*wavelength+spring_distance, water_level])) # ____ _ # / ___| |__ _ __ ___ _ __ ___ # | | | '_ \| '__/ _ \| '_ \ / _ \ # | |___| | | | | | (_) | | | | (_) | # \____|_| |_|_| \___/|_| |_|\___/ # Chrono # SYSTEM # create system system = fsi.ProtChSystem() # access chrono object chsystem = system.getChronoObject() # communicate gravity to system # can also be set with: # system.ChSystem.Set_G_acc(pychrono.ChVectorD(g[0], g[1], g[2])) system.setGravitationalAcceleration(g) # set maximum time step for system system.setTimeStep(1e-4) solver = pychrono.ChSolverMINRES() chsystem.SetSolver(solver) # BODY # create floating body
caisson.holes_ind = np.array([0]) tank.setChildShape(caisson, 0) # translate caisson to middle of the tank caisson.translate(np.array([1 * wavelength, water_level])) # ____ _ # / ___| |__ _ __ ___ _ __ ___ # | | | '_ \| '__/ _ \| '_ \ / _ \ # | |___| | | | | | (_) | | | | (_) | # \____|_| |_|_| \___/|_| |_|\___/ # Chrono # SYSTEM # create system system = fsi.ProtChSystem() # access chrono object chsystem = system.getChronoObject() # communicate gravity to system # can also be set with: # system.ChSystem.Set_G_acc(pychrono.ChVectorD(g[0], g[1], g[2])) system.setGravitationalAcceleration(g) # set maximum time step for system system.setTimeStep(1e-4) # set rate at which values are recorded system.setSampleRate(sampleRate) solver = pychrono.ChSolverMINRES() chsystem.SetSolver(solver) solver.SetMaxIterations(1000) solver.EnableWarmStart(True)
