From 2133cc4de44ab9a5557a57a4fb458f3bd185a8d0 Mon Sep 17 00:00:00 2001 From: xipenglyu Date: Mon, 17 Jul 2023 16:31:54 +0800 Subject: [PATCH] change slender_structure_dynamics.h line 268 and 269 --- .../solid_dynamics/slender_structure_dynamics.h | 6 ++++-- 1 file changed, 4 insertions(+), 2 deletions(-) diff --git a/src/shared/particle_dynamics/solid_dynamics/slender_structure_dynamics.h b/src/shared/particle_dynamics/solid_dynamics/slender_structure_dynamics.h index d1a2db0182..b21b014bfc 100644 --- a/src/shared/particle_dynamics/solid_dynamics/slender_structure_dynamics.h +++ b/src/shared/particle_dynamics/solid_dynamics/slender_structure_dynamics.h @@ -265,8 +265,10 @@ namespace SPH Vecd local_dpseudo_n_d2t = transformation_matrix_[index_i] * dpseudo_n_d2t_[index_i]; Vecd local_dpseudo_b_n_d2t = transformation_matrix_[index_i] * dpseudo_b_n_d2t_[index_i]; - dangular_b_vel_dt_[index_i] = getRotationFromPseudoNormalForSmallDeformation_b(Vec3d(local_dpseudo_b_n_d2t), Vec3d(local_dpseudo_n_d2t), Vec3d(rotation_b_[index_i]), Vec3d(angular_b_vel_[index_i]), dt); - dangular_vel_dt_[index_i] = getRotationFromPseudoNormalForSmallDeformation(Vec3d(local_dpseudo_b_n_d2t), Vec3d(local_dpseudo_n_d2t), Vec3d(rotation_[index_i]), Vec3d(angular_vel_[index_i]), dt); + dangular_b_vel_dt_[index_i] = Vec3d(0.0, 0.0, dpseudo_b_n_d2t[0]); + + //getRotationFromPseudoNormalForSmallDeformation_b(Vec3d(local_dpseudo_b_n_d2t), Vec3d(local_dpseudo_n_d2t), Vec3d(rotation_b_[index_i]), Vec3d(angular_b_vel_[index_i]), dt); + dangular_vel_dt_[index_i] = Vec3d(0.0, dpseudo_n_d2t[0], 0.0); // getRotationFromPseudoNormalForSmallDeformation(Vec3d(local_dpseudo_b_n_d2t), Vec3d(local_dpseudo_n_d2t), Vec3d(rotation_[index_i]), Vec3d(angular_vel_[index_i]), dt); }; void update(size_t index_i, Real dt = 0.0);