yade-dev team mailing list archive

[noreply@xxxxxxxxxxxxx]

------------------------------------------------------------
revno: 3780
committer: Francois <francois.kneib@xxxxxxxxx>
timestamp: Wed 2013-12-18 15:45:51 +0100
message:
  Corrected a minor mistake in code comments : k=2*r*E
modified:
  pkg/dem/Ip2_FrictMat_FrictMat_FrictPhys.cpp


--
lp:yade
https://code.launchpad.net/~yade-pkg/yade/git-trunk

Your team Yade developers is subscribed to branch lp:yade.
To unsubscribe from this branch go to https://code.launchpad.net/~yade-pkg/yade/git-trunk/+edit-subscription

=== modified file 'pkg/dem/Ip2_FrictMat_FrictMat_FrictPhys.cpp'
--- pkg/dem/Ip2_FrictMat_FrictMat_FrictPhys.cpp	2012-06-04 21:03:29 +0000
+++ pkg/dem/Ip2_FrictMat_FrictMat_FrictPhys.cpp	2013-12-18 14:45:51 +0000
@@ -37,7 +37,7 @@
 	Ra=sphCont->refR1>0?sphCont->refR1:sphCont->refR2;
 	Rb=sphCont->refR2>0?sphCont->refR2:sphCont->refR1;
 
-	//harmonic average of the two stiffnesses when (Ri.Ei/2) is the stiffness of a contact point on sphere "i"
+	//harmonic average of the two stiffnesses when (2*Ri*Ei) is the stiffness of a contact point on sphere "i"
 	Real Kn = 2*Ea*Ra*Eb*Rb/(Ea*Ra+Eb*Rb);
 	//same for shear stiffness
 	Real Ks = 2*Ea*Ra*Va*Eb*Rb*Vb/(Ea*Ra*Va+Eb*Rb*Vb);