New question #687702 on Yade:
https://answers.launchpad.net/yade/+question/687702
Hi there,
I am changing the damping force formula in the source code, and I want the
damping force to be proportional to incidentVn**1.5 instead of incidentVn. Am I
correct as follows?
Thanks!
double cn = Cn_crit*phys->viscousDamping; // Damping normal coefficient
/* Add normal viscous component if damping is included */
Vector3r incidentVn_Test=incidentVn;
double sign_1=Sgn(incidentVn_Test[0]);
double sign_2=Sgn(incidentVn_Test[1]);
double sign_3=Sgn(incidentVn_Test[2]);
#Here Sgn is a function to return the sign (1 or minus 1)
incidentVn_Test[0]=sign_1*pow(fabs(incidentVn_Test[0]),1.5);
incidentVn_Test[1]=sign_2*pow(fabs(incidentVn_Test[1]),1.5);
incidentVn_Test[2]=sign_3*pow(fabs(incidentVn_Test[2]),1.5);
#Here I change the component of the Vector3r respectively
phys->normalViscous = cn*incidentVn_Test;
/* Add normal viscous component if damping is included */
phys->normalForce -= phys->normalViscous;
--
You received this question notification because your team yade-users is
an answer contact for Yade.
_______________________________________________
Mailing list: https://launchpad.net/~yade-users
Post to : [email protected]
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
