← Back to team overview

yade-users team mailing list archive

Re: Contact forces


Hi Nejib,

Presently, I think plotting forces against orientations has not been done.  YADE is able to plot the histogram (number) of contact orientations.  https://yade-dem.org/index.php/ScriptPlotContactDirectionDistribution   You may have to code it yourself.  In principle, it would be same as the one above.  Instead of counting the number of contacts in that direction, you can sum the forces and average it.  

Alternatively, if you wish to store the data into a textfile, you can look at the TriaxialStateRecorder.cpp/hpp.  Also do look at TriaxialCompressionEngine.cpp/hpp.  To loop over the interactions, you can use:
FOREACH(const shared_ptr<Interaction>& I, *ncb->interactions){  }

The contact normal should be stored in ScGeom or Dem3Dof depending on which contact laws you are using.   I'm quite new to YADE, and this is all that I can suggest.



> From: nejib.hadda@xxxxxxxxxxx
> To: yade-users@xxxxxxxxxxxxxxxxxxx
> Date: Tue, 15 Dec 2009 15:20:37 +0100
> Subject: [Yade-users] Contact forces
> Hi evrybody,
> I want to draw mean contact forces against contact orientations for a
> triaxial test. I take a look into the contactstresses files, but there
> is no informations about the contact directions. 
> _______________________________________________
> Mailing list: https://launchpad.net/~yade-users
> Post to     : yade-users@xxxxxxxxxxxxxxxxxxx
> Unsubscribe : https://launchpad.net/~yade-users
> More help   : https://help.launchpad.net/ListHelp
New Windows 7: Simplify what you do everyday. Find the right PC for you.

Follow ups