yade-users team mailing list archive
-
yade-users team
-
Mailing list archive
-
Message #09127
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.