Describe the bug
A mistake in applying friction constraint at nodes. When applying friction constraint at 3D condition, if the friction constraint belongs to kinetic friction and current tangential velocity (vel_net_t) is smaller than velocity redution due to friction (vel_fricion), the new tangential velocity should be updated to zero. But it doesn't.
Mistake should be from node.tcc line 487 to 488
if (vel_net_t <= vel_fricion) {acc(dir_t0) = -vel(dir_t0); // To set particle velocity to zeroacc(dir_t1) = -vel(dir_t1);}
Expectedif (vel_net_t <= vel_fricion) {acc(dir_t0) = -vel(dir_t0 /dt); // To set particle velocity to zeroacc(dir_t1) = -vel(dir_t1) /dt;}
Describe the bug A mistake in applying friction constraint at nodes. When applying friction constraint at 3D condition, if the friction constraint belongs to kinetic friction and current tangential velocity (vel_net_t) is smaller than velocity redution due to friction (vel_fricion), the new tangential velocity should be updated to zero. But it doesn't.
Mistake should be from node.tcc line 487 to 488
if (vel_net_t <= vel_fricion) {
acc(dir_t0) = -vel(dir_t0); // To set particle velocity to zero
acc(dir_t1) = -vel(dir_t1);
}
Expected
if (vel_net_t <= vel_fricion) {
acc(dir_t0) = -vel(dir_t0 /dt); // To set particle velocity to zero
acc(dir_t1) = -vel(dir_t1) /dt;
}