|
by
Jeff Lander This article originally
appeared in the
|
The Trials and Tribulations of Tribology
Listing 1. Code for calculating and adding in friction.
// Calculate Magnitude
of Fn
FdotN = DotProduct(&curParticle->contactN,&curParticle->f);
// Calculate Vt Velocity
Tangent to Contact Normal
VdotN = DotProduct(&curParticle->contactN,&curParticle->v);
ScaleVector(&curParticle->contactN, VdotN, &Vn);
VectorDifference(&curParticle->v, &Vn, &Vt);
NormalizeVector(&Vt);
// Get the Direction of Vt
// Multiply By Normal force magnitude and Coef of Kinetic Friction
ScaleVector(&Vt, (FdotN * m_Ckf), &Vt);
// Add into the Force
Accumulator
VectorSum(&curParticle->f,&Vt,&curParticle->f);
________________________________________________________
