yade-users team mailing list archive

[Andreas Plesch <andreasplesch@xxxxxxxxx>]

Here is a simple diff which hopefully changes the
friction from average to lower of the two:

orig:

                       contactPhysics->frictionAngle  
                = (fa+fb)*0.5; // FIXME - this is
actually a waste of memory space, just like initialKs
and initialKn

new:

//AP: use smaller of the two friction angles, not
average
                       if (fa<fb)
{contactPhysics->frictionAngle       = fa;}
                       else
{contactPhysics->frictionAngle     = fb;}

seems to work, Andreas


diff -u
---
yade-packages/yade-package-dem/src/Engine/EngineUnit/MacroMicroElasticRelationships/MacroMicroElasticRelationships.cpp.orig
 2005-11-14 14:20:39.000000000 -0500
+++
yade-packages/yade-package-dem/src/Engine/EngineUnit/MacroMicroElasticRelationships/MacroMicroElasticRelationships.cpp
  2005-12-21 14:18:24.000000000 -0500
@@ -92,7 +92,10 @@
                        contactPhysics->initialKn     
                 = Kn;
                        contactPhysics->initialKs     
                 = Kn*(1-alpha*Vab)/(1+Vab);
 //cerr << "Ks: " <<       contactPhysics->initialKs  
                 << endl;
-                       contactPhysics->frictionAngle 
                 = (fa+fb)*0.5; // FIXME - this is
actually a waste of memory space, just like initialKs
and initialKn
+//                     contactPhysics->frictionAngle 
                 = (fa+fb)*0.5; // FIXME - this is
actually a waste of memory space, just like initialKs
and initialKn
+//AP: use smaller of the two friction angles, not
average
+                       if (fa<fb)
{contactPhysics->frictionAngle       = fa;}
+                       else
{contactPhysics->frictionAngle     = fb;}
                       
contactPhysics->tangensOfFrictionAngle          =
std::tan(contactPhysics->frictionAngle);

                        contactPhysics->prevNormal    
                 = interactionGeometry->normal;


--- Janek Kozicki <janek_listy@xxxxx> wrote:

> > seems to indicate that the average friction is
> > calculated and actually used. Is this
> understanding
> > correct ? Would this be the correct file to change
> > that behaviour, eg. from average friction to the
> lower
> > of the two ?
> 
> yes, the average is calculated, (f1+f2)/2. If you
> want different - modify
> there.
> 
> Thanks for all your posts about XML, this week I'll
> answer them with some
> additional comments, as now I'm short on time.
> (working on lattice, yade
> docs, and laboratory experiments with concrete).
> 
> -- 
> Janek Kozicki                                       
>                  |
> _______________________________________________
> Yade-users mailing list
> Yade-users@xxxxxxxxxxxxxxxx
> http://lists.berlios.de/mailman/listinfo/yade-users
> 

_______________________________________________
Yade-users mailing list
Yade-users@xxxxxxxxxxxxxxxx
http://lists.berlios.de/mailman/listinfo/yade-users