yade-dev team mailing list archive

Thread
Date

[Branch ~yade-dev/yade/trunk] Rev 2803: Fix small bug in HM (only in case of viscous damping).

To: Yade developers <yade-dev@xxxxxxxxxxxxxxxxxxx>
From: noreply@xxxxxxxxxxxxx
Date: Tue, 05 Apr 2011 16:24:24 -0000
Reply-to: noreply@xxxxxxxxxxxxx
Sender: bounces@xxxxxxxxxxxxx

------------------------------------------------------------
revno: 2803
committer: Chiara Modenese <c.modenese@xxxxxxxxx>
branch nick: yade
timestamp: Tue 2011-04-05 17:19:20 +0000
message:
  Fix small bug in HM (only in case of viscous damping).
modified:
  pkg/dem/HertzMindlin.cpp


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

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

=== modified file 'pkg/dem/HertzMindlin.cpp'
--- pkg/dem/HertzMindlin.cpp	2011-02-08 16:37:12 +0000
+++ pkg/dem/HertzMindlin.cpp	2011-04-05 17:19:20 +0000
@@ -369,6 +369,7 @@
 	/********************/
 	
 	phys->shearViscous=Vector3r::Zero(); // reset so that during sliding, the previous values is not there
+	Fn = phys->normalForce.norm();
 	if (!includeAdhesion) {
 		Real maxFs = Fn*phys->tangensOfFrictionAngle;
 		if (shearElastic.squaredNorm() > maxFs*maxFs){