yade-users team mailing list archive
-
yade-users team
-
Mailing list archive
-
Message #01286
Re: Probleme avec un mur à déplacer sous Python
-
To:
Jean-Francois Jerier <jerier@xxxxxxxxxxxxxxx>
-
From:
Václav Šmilauer <eudoxos@xxxxxxxx>
-
Date:
Wed, 29 Oct 2008 19:34:42 +0100
-
Cc:
yade-users@xxxxxxxxxxxxxxxxxxx
-
In-reply-to:
<49089B77.4070308@geo.hmg.inpg.fr>
-
User-agent:
Mozilla/5.0 (X11; U; Linux x86_64; en-US; rv:1.8.1.17) Gecko/20080925 Lightning/0.8 Thunderbird/2.0.0.17 Mnenhy/0.7.5.666
> Je ne comprends pas pourquoi mon mur du dessus qui a une vitesse nulle
> se deplace à l'infini.
> Pourtant dans mon fichier j'ai bien mis:
>
> o.bodies[41]['isDynamic']=True; #Premiere facette
> o.bodies[41].phys['velocity']=[0,0,0]; #Deuxième facette
> o.bodies[42]['isDynamic']=True; #Premiere facette
> o.bodies[42].phys['velocity']=[0,0,0] #Deuxième facette
>
>
> DeusExMachina('TranslationEngine',{'subscribedBodies':[41],'velocity':0,'translationAxis':[0,0,-1]}),
> #Premiere facette
>
> DeusExMachina('TranslationEngine',{'subscribedBodies':[42],'velocity':0,'translationAxis':[0,0,-1]}),
> #Deuxième facette
<footnote>
You can do this in a more efficient way:
DeusExMachina('TranslationEngine',{'subscribedBodies':[41,42],'velocity':0,'translationAxis':[0,0,-1]}),
</footnote>
Well, it actually doesn't make much sense to have the body dynamic and
prescribe translation at the same time, does it? The problem of
displacement being infinity is that facets have zero mass (see
o.bodies[41].phys['mass']) and that doesn't do any good to the integrator.
I think that for the wall that does the compression, it would be perhaps
better to use regular box instead of facets; if you have it with
wire=True, you will even still see what is goig on down there.
I am not sure what you're trying to achieve? A box compressing spheres
from above, ok, but if it has contact speed, it will explode at one
point, since no force will stop that box. Wouldn't it be better to
actually prescribe force? But force will be constant, if you want to
change at runtime, you would need something more sophisticated, like
TriaxialTest, but only 1d.
HTH, Vaclav