def evaluate_buoy_forces(buoy_object):
# Set up radiation and diffraction problems
problems = [cpt.RadiationProblem(body=buoy_object, radiating_dof=degree_of_freedom, omega=omega)
for omega in omega_range]
problems += [cpt.DiffractionProblem(omega=omega, body=buoy_object, wave_direction=wave_direction)
for omega in omega_range]
# Solve each matrix problem
solver = cpt.BEMSolver(engine=cpt.HierarchicalToeplitzMatrixEngine())
results = [solver.solve(pb) for pb in sorted(problems)]
return cpt.assemble_dataset(results)
return 0.1*(-(z+1)**2 + 16)
# Generate the mesh and display it with VTK.
buoy = cpt.FloatingBody(
cpt.AxialSymmetricMesh.from_profile(shape, z_range=np.linspace(-5, 0, 30), nphi=40)
)
buoy.add_translation_dof(name="Heave")
buoy.show()
# Set up problems
omega_range = np.linspace(0.1, 5.0, 60)
problems = [cpt.RadiationProblem(body=buoy, radiating_dof='Heave', omega=omega)
for omega in omega_range]
# Solve the problems using the axial symmetry
solver = cpt.BEMSolver(engine=cpt.HierarchicalToeplitzMatrixEngine())
results = [solver.solve(pb) for pb in sorted(problems)]
dataset = capytaine.io.xarray.assemble_dataset(results)
# Plot results
import matplotlib.pyplot as plt
plt.figure()
plt.plot(
omega_range,
dataset['added_mass'].sel(radiating_dof='Heave',
influenced_dof='Heave'),
label="Added mass",
)
plt.plot(
omega_range,
dataset['radiation_damping'].sel(radiating_dof='Heave',