yade-users team mailing list archive

[Seungcheol Yeom <question299971@xxxxxxxxxxxxxxxxxxxxx>]

New question #299971 on Yade:
https://answers.launchpad.net/yade/+question/299971

I have tried to simulate a simple script below using both yade (1.07) and yadedaily (2016.06a-2-86aa164~trusty).
However, the results are very different.
Is this a bug?

from yade import pack,geom,plot,qt

#tag = O.tags['d.id']
#os.mkdir(tag)

r = 1.e-3 #particle radius, meter
h = 0.0 #praticle distance, meter

soilID=FrictMat(young=1.e6,poisson=0.4,frictionAngle=radians(30),density=2650.0,label='soil')
O.materials.append(soilID)

#create two sphere paticles#
O.bodies.append([
   utils.sphere(center=(0,0,-(r+h*0.5)),radius=r,material=soilID,fixed=False),
   utils.sphere(center=(0,0,r+h*0.5),radius=r,material=soilID,fixed=False)
])



#define engines#
O.engines=[
   ForceResetter(),
   InsertionSortCollider([Bo1_Sphere_Aabb()]),
   InteractionLoop(
      [Ig2_Sphere_Sphere_ScGeom()],
      [Ip2_FrictMat_FrictMat_MindlinCapillaryPhys()],
      [Law2_ScGeom_MindlinPhys_Mindlin(neverErase=True)]
   ),
   Law2_ScGeom_CapillaryPhys_Capillarity(capillaryPressure=100.e3,createDistantMeniscii=False,label='Cap'),
   NewtonIntegrator(damping=0.2,gravity=(0,0,0)),
   PyRunner(iterPeriod=1,command='check_cap()')
]

def check_cap():
	global Vol_meniscus,chk_meniscus
	i = O.interactions[0,1]
	Vol_meniscus = i.phys.vMeniscus #total volume of meniscus
	chk_meniscus = i.phys.meniscus 
	mov1 = O.bodies[0].state.pos[2]
	mov2 = O.bodies[1].state.pos[2]
	
	print Vol_meniscus,chk_meniscus
	plot.addData(i=O.iter,id1_mov=mov1,id2_mov=mov2)
#	plot.saveDataTxt(O.tags['id']+'.txt')
    
#display the force on the sphere#


plot.plots={'i':('id2_mov')}
plot.plot()
plot.live=True
plot.autozoom=True
qt.View()   
O.dt=0.5*PWaveTimeStep()
O.run(10,True)

Thanks for your help.

Seungcheol

-- 
You received this question notification because your team yade-users is
an answer contact for Yade.