yade-users team mailing list archive

[Jerome Duriez <duriez@xxxxxxxxxxxxxxx>]

Unfortunately I know very very few python, so I do not know how to do 
this in this framework. With C++ instructions, (inside the code) it is 
like that

dynamic_cast<Force*>(*Metabody-*>physicalActions->find(*id_facet*,actionForce->getClassIndex()) 
. get() )->force;
(I see clearly what can be the advantages of python by writing this...)

   An intermediate mean for you is (until someone can answer you more 
precisely) to add an Engine "ForceRecorder" (source in 
pkg/dem/Engine/StandAloneEngine/), defined with the id of the facet, to 
o.engines. This will store the force you want in an output file you can 
later analyse with the software you want.

Jerome

Anton Gladky a écrit :
> How can I get forces affected on 1 facet from spheres: I need it for plot?
>
> Thank you
>
> PS I promise, that I will place all answers to instructions or FAQ.
> ______________________________
> [ENG] Best Regards
> [GER] Mit freundlichen Grüßen
> [RUS] С наилучшими пожеланиями
> [UKR] З найкращими побажаннями
>
> Anton Gladkyy
>
>
> 2009/6/26 Anton Gladky <gladky.anton@xxxxxxxxx 
> <mailto:gladky.anton@xxxxxxxxx>>
>
>     Amazing, but it works!
>     Thank you!
>
>     ______________________________
>     [ENG] Best Regards
>     [GER] Mit freundlichen Grüßen
>     [RUS] С наилучшими пожеланиями
>     [UKR] З найкращими побажаннями
>
>     Anton Gladkyy
>
>
>     2009/6/26 Václav Šmilauer <eudoxos@xxxxxxxx
>     <mailto:eudoxos@xxxxxxxx>>
>
>
>             pp={'se3':[center[0],center[1],center[2],1,0,0,0],'refSe3':[center[0],center[1],center[2],1,0,0,0],'young':young,'poisson':poisson,'frictionAngle':frictionAngle,'inertia':[0,0,0],
>             'velocity':[velocity[0],velocity[1],velocity[2]]}
>
>         You have to use TranslationEngine and subscribe facets to it.
>         That's what I do. See e.g. scripts/test/Dem3DofGeom.py for an
>         example. The sphere is moved that way. Facet ids are returned
>         by O.bodies.append.
>
>         ids=O.bodies.append([utils.facet(...),utils.facet(...)])
>
>         print ids ## [0,1,2,3,4]
>         O.engines=[...,TranslationEngine(subscribedBodies=ids,translationAxis=[0,0,1],velocity=1),...]
>
>         HTH, Vaclav
>
>
>
>         _______________________________________________
>         Mailing list: https://launchpad.net/~yade-users
>         <https://launchpad.net/%7Eyade-users>
>         Post to     : yade-users@xxxxxxxxxxxxxxxxxxx
>         <mailto:yade-users@xxxxxxxxxxxxxxxxxxx>
>         Unsubscribe : https://launchpad.net/~yade-users
>         <https://launchpad.net/%7Eyade-users>
>         More help   : https://help.launchpad.net/ListHelp
>
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> 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
>