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.
_______________________________________________
Mailing list: https://launchpad.net/~yade-users
Post to : [email protected]
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
