------------------------------------------------------------ revno: 3798 committer: Klaus Thoeni <[email protected]> timestamp: Wed 2014-01-22 16:07:34 +1100 message: Correct shear stiffness averaging for Ip2_FrictMat_FrictMat_ViscoFrictPhys 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 2014-01-08 14:36:59 +0000 +++ pkg/dem/Ip2_FrictMat_FrictMat_FrictPhys.cpp 2014-01-22 05:07:34 +0000 @@ -91,7 +91,7 @@ //harmonic average of the two stiffnesses when (Ri.Ei/2) 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*Va); + Real Ks = 2*Ea*Ra*Va*Eb*Rb*Vb/(Ea*Ra*Va+Eb*Rb*Vb); Real frictionAngle = (!frictAngle) ? std::min(mat1->frictionAngle,mat2->frictionAngle) : (*frictAngle)(mat1->id,mat2->id,mat1->frictionAngle,mat2->frictionAngle); contactPhysics->tangensOfFrictionAngle = std::tan(frictionAngle);
_______________________________________________ Mailing list: https://launchpad.net/~yade-dev Post to : [email protected] Unsubscribe : https://launchpad.net/~yade-dev More help : https://help.launchpad.net/ListHelp
