yade-users team mailing list archive
-
yade-users team
-
Mailing list archive
-
Message #06732
Re: [Question #215627]: Translation and rotation engines
Question #215627 on Yade changed:
https://answers.launchpad.net/yade/+question/215627
Status: Open => Answered
Jan Stránský proposed the following answer:
Hello,
few comments:
next time please add to the email whole error code you get, it will make it
much more easier to find the problem and to help you. In case ot
TranslationEngine and RotationEngine I would say that the problem is in
parameter ids, which should be list of integers, while in your case it is
Body instance, so write ids=[plate.id] or ids=[LW.id,RW.id] instead.
You wrote about error while using plate.state.vel. Could you provide us
with a short script producing this error?
Some more comments, not causing error but making the script nicer:
- in checkUnbalanced function, don't add GravitiyEngine. The existing one
(defined in first O.engines) should work for new particles and is placed
before NewtonIntegrator, which is correct.
- don't use "while m in range(len(O.bodies))", use "while m <
len(O.bodies)" or just "for b in O.bodies". Python does need to construct
range(len(O.bodies)) list every loop iteration.
- instead of try except pair in addfines function use "if b is not None".
cheers
Jan
2012/11/30 m.asd <question215627@xxxxxxxxxxxxxxxxxxxxx>
> New question #215627 on Yade:
> https://answers.launchpad.net/yade/+question/215627
>
> hi,
> I am a new user of yade and i am going to model a box in which some
> particles deposit. then I delete particles upper than a predefined
> elevation and add some fine particles. then I delete these fines upper than
> that elevation too. after that, I add a plate on top and define a velocity
> up to a specific load. set the velocity of the plate to zero (fix place)
> and then let sidewalls to rotate around the Y axis with a fixed angVel.
> here is my problem,
> when I use "plate.state.vel" it says that "Nonetype object has no
> attribute state"
> and when I use rotationengine and translationengine it shows an error and
> do not proceed.
> as I am not so experienced in running yade, changing the script in many
> ways made me frustrated. I will be so grateful if anyone helps me.
> furthermore, It sometimes stops running with this error
> "segmentation fault (core dumped)" what is the reason
> here is my script
>
>
> overburden=10
> from yade import pack,plot,qt
>
> idmat=O.materials.append(FrictMat(density=2700,young=1e6,poisson=0.5,frictionAngle=0))
> O.bodies.append(utils.wall(position=0,axis=2,sense=1))
> O.bodies.append(utils.wall(position=0,axis=1,sense=1))
> O.bodies.append(utils.wall(position=0.01,axis=1,sense=-1))
> O.bodies.append(utils.wall(position=0,axis=0,sense=1))
> LW=O.bodies[-1]
> O.bodies.append(utils.wall(position=0.01,axis=0,sense=-1))
> RW=O.bodies[-1]
> sp=pack.SpherePack()
> sp.makeCloud((0,0,0.001),(0.01,0.01,0.02),rMean=.0008,rRelFuzz=0)
> sp.toSimulation(color=(0,0,1)) # blue
>
> O.engines=[
> ForceResetter(),
> InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Wall_Aabb()]),
> InteractionLoop(
> [Ig2_Sphere_Sphere_ScGeom(),Ig2_Wall_Sphere_ScGeom()],
> [Ip2_FrictMat_FrictMat_FrictPhys()],
> [Law2_ScGeom_FrictPhys_CundallStrack()]
> ),
> GravityEngine(gravity=(0,0,-9.81)),
> NewtonIntegrator(damping=0.4),
>
> qt.SnapshotEngine(fileBase='3d-',iterPeriod=100000,label='snapshot'),
> PyRunner(command='unbalancedData()',iterPeriod=30000),
> PyRunner(command='checkUnbalanced()',realPeriod=1,label='checker'),
> ]
>
> O.dt=1e-6
>
> def checkUnbalanced():
> if O.iter<500000: return
> if utils.unbalancedForce()>.08: return
> for b in O.bodies:
> if isinstance(b.shape,Sphere):
> if (b.state.pos[2]+b.shape.radius)>0.007:
> O.bodies.erase(b.id)
> fines=pack.SpherePack()
>
> fines.makeCloud((0,0,0.007),(0.01,0.01,0.01),rMean=.0004,rRelFuzz=0)
> fines.toSimulation()
> O.engines=O.engines+[GravityEngine(gravity=(0,0,-9.81),mask=-1)]
> checker.command='addfines()'
>
> def addfines():
> if O.iter<650000: return
> if utils.unbalancedForce()>.08: return
> for b in O.bodies:
> if isinstance(b.shape,Sphere):
> if (b.state.pos[2]+b.shape.radius)>0.007:
> O.bodies.erase(b.id)
> posi=O.bodies[5].state.pos[2]+O.bodies[5].shape.radius
> m=5
> while m in range(len(O.bodies)):
> b=O.bodies[m]
> try:
> if isinstance(b.shape,Sphere):
> posi2=b.state.pos[2]+b.shape.radius
> if posi2>=posi: posi=posi2
> m=m+1
> except: m=m+1
> O.bodies.append(utils.wall(position=posi,axis=2,sense=-1))
> global plate
> plate=O.bodies[-1]
>
> O.engines=O.engines[:4]+[TranslationEngine(ids=plate,translationAxis=(0,0,1),velocity=0.002,label="transEngine")]+O.engines[4:]
> #here is the problem
> checker.command='shearing()'
> O.materials[0].frictionAngle=.5
> for i in O.interactions: i.phys.tangensOfFrictionAngle=tan(.5)
>
> def shearing():
> if abs(O.forces.f(plate.id)[2])<overburden: return
> transEngine.velocity=0
>
> O.engines=O.engines[:5]+[RotationEngine(ids=[LW,RW],angularVelocity=0.005,rotationAxis=(0,1,0),label="rotEngine")]+O.engines[5:]
> #and here too
> checker.command='stopShearing()'
>
> def stopShearing():
> if O.bodies[3].state.rot<(0,0.005,0): return
> O.pause()
> utils.makeVideo(snapshot.snapshots,'3dshear.mpeg',fps=10,bps=10000)
> print 'Finished'
>
> def unbalancedData():
> plot.addData(itperi=O.iter,unbalanced=utils.unbalancedForce())
>
> plot.plots={'itperi':('unbalanced',)}
> plot.plot()
>
> Gl1_Sphere.stripes=True
> qt.View()
>
> rr=yade.qt.Renderer()
> rr.shape=True
> O.saveTmp()
>
>
> --
> You received this question notification because you are a member of
> yade-users, which is an answer contact for Yade.
>
> _______________________________________________
> 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
>
--
You received this question notification because you are a member of
yade-users, which is an answer contact for Yade.
