← Back to team overview

yade-users team mailing list archive

[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.