← Back to team overview

yade-users team mailing list archive

Re: [Question #245072]: exciting box

 

Question #245072 on Yade changed:
https://answers.launchpad.net/yade/+question/245072

    Status: Answered => Open

masmoudimarwa is still having a problem:
HIii,
Thansk a lot

## PhysicalParameters
Young = 30e9
Poisson = 0.2
Density=1000
# Append a material
mat=O.materials.append(FrictMat(young=Young,poisson=Poisson,density=Density,frictionAngle=.5))



O.bodies.append(utils.geom.facetBox((.5,.5,.5),(.5,.5,.5),wallMask=63,material=mat))
O.bodies.append([sphere([0.5,0.5,0.5],.05,material=mat)])
O.bodies[12].state.vel=(0,0,13)
b=O.bodies[0]
bb=O.bodies[1]
b.dynamic=True
bb.dynamic=True

O.engines=[
   ForceResetter(),
   InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]),
   InteractionLoop(
      # handle sphere+sphere and facet+sphere collisions
      [Ig2_Sphere_Sphere_L3Geom(),Ig2_Facet_Sphere_L3Geom()],
      [Ip2_FrictMat_FrictMat_FrictPhys()],
      [Law2_L3Geom_FrictPhys_ElPerfPl()]
   ),
   ForceEngine(force=(0,0,10),ids=[0]),
   ForceEngine(force=(0,0,10),ids=[1]),
   GravityEngine(gravity=(0,0,-9.81)),
   NewtonIntegrator(damping=0.4),
]


O.dt=.5*utils.PWaveTimeStep()


O.saveTmp()
O.run()
But it doesn't work and i receive this message
FATAL /build/buildd/yade-0.70.0/core/ThreadRunner.cpp:31 run: Exception
occured:
Body #0 has velocity==NaN!
 Thanks a lot for help


2014-03-06 13:06 GMT+01:00 Jan Stránský <
question245072@xxxxxxxxxxxxxxxxxxxxx>:

> Your question #245072 on Yade changed:
> https://answers.launchpad.net/yade/+question/245072
>
>     Status: Open => Answered
>
> Jan Stránský proposed the following answer:
> Hello,
> the problem is that facets are non-dynamic bodies by default (therefore not
> affected by forces). Try to use some KinematicEngine as a first attempt, if
> you are successful with it, it should be possible to introduce also forces.
> cheers
> Jan
>
>
>
> 2014-03-06 12:46 GMT+01:00 masmoudimarwa <
> question245072@xxxxxxxxxxxxxxxxxxxxx>:
>
> > New question #245072 on Yade:
> > https://answers.launchpad.net/yade/+question/245072
> >
> > Hello all,
> > I just try to model a box partially filled with particle and i want to
> > excite this box with a force or with displacement that affect the
> movement
> > of particle.
> > As a first step I model a box with one particle but my model don't make
> > the box in movement .
> > I have no idea where is wrong
> > Thanks
> > Best regards
> >
> >
> O.materials.append(FrictMat(young=30e9,density=1000,poisson=.2,frictionAngle=.5))
> >
> > # create rectangular box from facets
> > O.bodies.append(utils.geom.facetBox((.5,.5,.5),(.5,.5,.5),wallMask=63))
> > O.bodies.append([sphere([0.5,0.5,0.5],.05)])
> > O.bodies[12].state.vel=(0,0,3)
> >
> >
> >
> > O.engines=[
> >    ForceResetter(),
> >    InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]),
> >    InteractionLoop(
> >       # handle sphere+sphere and facet+sphere collisions
> >       [Ig2_Sphere_Sphere_L3Geom(),Ig2_Facet_Sphere_L3Geom()],
> >       [Ip2_FrictMat_FrictMat_FrictPhys()],
> >       [Law2_L3Geom_FrictPhys_ElPerfPl()]
> >    ),
> >    GravityEngine(gravity=(0,0,-9.81)),
> >    ForceEngine(ids=[0],force=(10,0,0)),
> >    ForceEngine(force=(10,0,0),ids=[0]),
> >    ForceEngine(force=(10,0,0),ids=[1]),
> >    NewtonIntegrator(damping=0.4),
> > ]
> > O.dt=.5*utils.PWaveTimeStep()
> >
> > --
> > 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
> >
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
> https://answers.launchpad.net/yade/+question/245072/+confirm?answer_id=2
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/yade/+question/245072
>
> You received this question notification because you asked the question.
>

-- 
You received this question notification because you are a member of
yade-users, which is an answer contact for Yade.