← Back to team overview

yade-dev team mailing list archive

Re: Bugs in Cell

 

I have added

for i in range(len(O.bodies)):
	O.bodies[i].state.pos = O.cell.trsf*O.bodies[i].state.pos

which located in proper place seems to work well.

If you have only prescribed stress, then rotating stress is the easiest way, of course. But if you want to combine prescribed strain and stress in different directions, lets say prescribed strain along x axis and zero stress in y and z axis (uniaxial stress), rotating these boundary conditions is not trivial. On the other hand, saving boundary conditions and only rotating the cell solves the problem.

Thank Bruno and Vaclav for discussion :-)

Jan

______________________________________________________________
> Od: "Václav Šmilauer" <eudoxos@xxxxxxxx>
> Komu: yade-dev@xxxxxxxxxxxxxxxxxxx
> Datum: 30.06.2010 12:12
> Předmět: Re: [Yade-dev] Bugs in Cell
>
>> If you really want large and strict rigid body rotation in one step, my 
>> script is not exactly doing the job (you tried it yes?).
>The packing rotation is trivial
>https://yade-dem.org/sphinx/yade.pack.html#yade._packSpheres.SpherePack.rotate it is even easier that rotating the stress tensor, I think.
>
>Jan, I get it, of course, thanks :-)
>
>> Why not rotating the loading instead? It is equivalent and needs 0 lines 
>> of code.
>
>
>
>_______________________________________________
>Mailing list: https://launchpad.net/~yade-dev
>Post to     : yade-dev@xxxxxxxxxxxxxxxxxxx
>Unsubscribe : https://launchpad.net/~yade-dev
>More help   : https://help.launchpad.net/ListHelp
>



Follow ups

References