yade-users team mailing list archive
-
yade-users team
-
Mailing list archive
-
Message #07930
Re: [Question #232941]: Capillary pressure between two particles
Question #232941 on Yade changed:
https://answers.launchpad.net/yade/+question/232941
Seungcheol Yeom posted a new comment:
Hello all,
First of all, thanks for all your response. I have learned a lot from you.
By the way, I have modified the script as shown below:
from yade import plot,qt
r = 1e-4 #particle radius
h = 1e-5 #praticle distance
#create two sphere paticles#
O.bodies.append([
utils.sphere(center=(0,0,0),radius=r,fixed=False),
utils.sphere((0,0,2*r+h),r)
])
#define engines#
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb()]),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom()],
[Ip2_FrictMat_FrictMat_CapillaryPhys()],
[Law2_ScGeom_FrictPhys_CundallStrack(neverErase=True)]
),
Law2_ScGeom_CapillaryPhys_Capillarity(capillaryPressure=10000,label='Cap'),
NewtonIntegrator(damping=0.4,gravity=(0,0,0)),
PyRunner(iterPeriod=100,command='addPlotData()')
]
Cap.createDistantMeniscii=True
O.run(1,True)
Cap.createDistantMeniscii=False
def addPlotData():
plot.addData(i=O.iter,b0displ=O.bodies[0].state.displ(),b1displ=O.bodies[1].state.displ())
plot.plots={'i':('b0displ')}
plot.plot()
plot.live=True
plot.autozoom=True
O.dt=0.5*PWaveTimeStep()
qt.View()
O.saveTmp()
When I ran it, I am not able to see any displacements based on the plot. I hoped I can see the displacement as a function of the number of iteration but it seems something went wrong.
Am I making something wrong?
Thanks for your help!
Seungcheol
--
You received this question notification because you are a member of
yade-users, which is an answer contact for Yade.