yade-users team mailing list archive
-
yade-users team
-
Mailing list archive
-
Message #17997
[Question #672525]: Free fall simulation did not have good results
New question #672525 on Yade:
https://answers.launchpad.net/yade/+question/672525
Hello, I am a beginner of Yade. I am trying to couple the dem with LBM (Using Bouzidi boundary condition, some work like Franck Lominé's model). But when I test the free fall simulation, the Yade did not show a good result to me. I'm a little confused.
The following is the code I typed in yade:
---------------------------------------------------------------------------------------------------------------------
s=utils.sphere(center=[0,0,0],radius=1)
O.bodies.append(s)
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Wall_Aabb()]),
InteractionLoop(
[Ig2_Sphere_Sphere_L3Geom_Inc(),Ig2_Wall_Sphere_L3Geom_Inc()],
[Ip2_FrictMat_FrictMat_FrictPhys()],
[Law2_L3Geom_FrictPhys_ElPerfPl()]
),
NewtonIntegrator(damping=0,gravity=(0,0.,-9.81),label="ENewton"),
]
O.dt=1
O.bodies[0].state.vel #Result: Vector3(0,0,0)
O.bodies[0].state.pos #Result: Vector3(0,0,0)
O.step() #1st iter
O.bodies[0].state.vel #Result: Vector3(0,0,-9.81)
O.bodies[0].state.pos #Result: Vector3(0,0,-9.81) this doesn't make sense, I think it should be (0,0,-4.905)
O.step() #2nd iter
O.bodies[0].state.vel #Result: Vector3(0,0,-19.62)
O.bodies[0].state.pos #Result: Vector3(0,0,-29.43) this doesn't make sense, I think it should be (0,0,-19.62)
O.step() #3rd iter
O.bodies[0].state.vel #Result: Vector3(0,0,-29.43)
O.bodies[0].state.pos #Result: Vector3(0,0,-58.86) this doesn't make sense, I think it should be (0,0,-44.145)
---------------------------------------------------------------------------------------------
Thanks in advance, I appreciate your help.
Xiaohu
--
You received this question notification because your team yade-users is
an answer contact for Yade.