New question #704035 on Yade:
https://answers.launchpad.net/yade/+question/704035
Hi all,
I am trying to continue my simulation in a rotating drum after I stop it. To do
this I save my first simulation using O.save at the end of my first script. I
also export vtu files for visualisation in Paraview. This works correctly. Then
I O.load the saved file (test.bz2) at the beginning of the second script. The
simulation goes on and gives me vtu files. However, in the vtu files some
points are Nan and I do not get a proper visualisation in Paraview (it does not
show the particles in paraview). I was wondering if you could help me with this
issue.
I have attached my two scripts here as well :
Script 1:
from yade import ymport
from yade import pack
from yade import utils, plot,wrapper
from yade import export
nCG=1
fr = 0.38
rho = 3700
En = 0.22
Et = 0.22
poi=0.3
yoM=200e6
r = 0.0005
Gamma = 0.072
Theta = 0
vB = 6.03 * 1e-9 * (nCG**3)
Eta= 0.001
CapType="Rabinovich"
Tc=0.00045
##defining material
mat=O.materials.append(
ViscElCapMat(frictionAngle=fr, density=rho, poisson=poi, young= yoM,
Vb=vB, gamma=Gamma, eta=Eta, theta=Theta, Capillar=True,
CapillarType=CapType,en=En, et=Et,tc=Tc)
)
mat2=O.materials.append(
ViscElCapMat(frictionAngle=fr, density=rho, poisson=poi, young= yoM,
Vb=vB, gamma=Gamma, eta=0, theta=90, Capillar=True, CapillarType=CapType,
en=En, et=Et,tc=Tc)
)
#defining the spheres
sp=pack.SpherePack()
sp.makeCloud((-0.035,-0.035,-0.01),(0.035,0.035,0.01),psdSizes=[0.001,0.002,0.003,0.004],
psdCumm=[0.0,0.8,0.95,1.0])
sp.toSimulation(material=mat)
Nprtcl=len(O.bodies)
print(Nprtcl)
Tt= utils.PWaveTimeStep()
O.dt = 0.5*Tt
print(O.dt)
#liquidMigration
VV=0.03
Vmin=0.0
for s in O.bodies:
if not type(s.shape)==wrapper.Sphere:
continue
s.state.Vf=VV * (4/3) * 3.14*(s.shape.radius)**3
s.state.Vmin=Vmin
##
Drum=geom.facetCylinder(material=mat2,center=(0.0,0.0,0.0),radius=0.05,height=0.03,orientation=Quaternion(Vector3(0,0,1),(pi/2.0)))
walls = O.bodies.append(Drum)
##engine
O.engines = [
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb(), Bo1_Facet_Aabb()]),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom()],
[Ip2_ViscElCapMat_ViscElCapMat_ViscElCapPhys()],
[Law2_ScGeom_ViscElCapPhys_Basic()],
),
NewtonIntegrator(gravity=[0, -9.8, 0]),
RotationEngine(ids=walls,rotationAxis=[0,0,1],rotateAroundZero=True,
zeroPoint=[0.0,0.0,0.0], angularVelocity=2.0),
LiqControl(label='lqc'),
VTKRecorder(iterPeriod=1000,recorders=['spheres','facets','colors'],fileName='test-')]
O.run( 500000,True)
O.save('test.bz2')
Script 2:
from yade import ymport
from yade import pack
from yade import utils, plot,wrapper
from yade import export
O.load('test.bz2')
Nprtcl=len(O.bodies)
print(Nprtcl)
O.run(15000,True)
###
BTW, I am using this version yade-2022-04-13.git-5bd3ade.
Thanks in advance for your respopnse!
--
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 : yade-users@lists.launchpad.net
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
