correct forces sign
@@ -776,15 +776,15 @@ class FloatingForces(generator_interface.BaseGenerator):
if self.settings['method_matrices_freq'] == 'constant':
hd_f_qdot_g -= np.dot(self.hd_damping, self.qdot[data.ts, :])
- hd_f_qdotdot_g = np.dot(self.hd_added_mass, self.qdotdot[data.ts, :])
+ hd_f_qdotdot_g = -np.dot(self.hd_added_mass, self.qdotdot[data.ts, :])
equiv_hd_added_mass = self.hd_added_mass
elif self.settings['method_matrices_freq'] in ['interp_matrices', 'rational_function']:
- hd_f_qdot_g += response_freq_dep_matrix(self.hd_damping, self.ab_freq_rads, self.qdot, data.ts, self.settings['dt'])
- hd_f_qdotdot_g = response_freq_dep_matrix(self.hd_added_mass, self.ab_freq_rads, self.qdotdot, data.ts, self.settings['dt'])
+ hd_f_qdot_g -= response_freq_dep_matrix(self.hd_damping, self.ab_freq_rads, self.qdot, data.ts, self.settings['dt'])
+ hd_f_qdotdot_g = -response_freq_dep_matrix(self.hd_added_mass, self.ab_freq_rads, self.qdotdot, data.ts, self.settings['dt'])
# Compute the equivalent added mass matrix
- equiv_hd_added_mass = compute_equiv_hd_added_mass(hd_f_qdotdot_g, self.qdotdot[data.ts, :])
+ equiv_hd_added_mass = compute_equiv_hd_added_mass(-hd_f_qdotdot_g, self.qdotdot[data.ts, :])
# Update added mass in the beam
if self.added_mass_in_mass_matrix:
