Question #669990 on Yade changed:
https://answers.launchpad.net/yade/+question/669990
Status: Answered => Open
rhaven is still having a problem:
Hi Bruno,
thank you for the reply. I try with other packings but cant get around the
Vh=NULL error
for example the code below with the following input file
https://pastebin.com/8a2jsCqV
I understand that the error isnt necessarily coming from saveVTK, should
I open a new question instead? Is there another way to access the
triangulation data? My goal is actually to be able to read the edges and
vertices of the regular triangulation.
code:
O.reset()
sphs = ymport.text('clump' + str(clumpNum) + '.txt')
O.bodies.append(sphs)
print('Number of particles: ', len(sphs))
##########
raw_input('Before FlowEngine saveVtk')
O.engines = [
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Box_Aabb()]),
#InteractionLoop(
#[Ig2_Sphere_Sphere_ScGeom(),Ig2_Box_Sphere_ScGeom()],
##triax only
# [Ig2_Sphere_Sphere_ScGeom()],
# [Ip2_FrictMat_FrictMat_FrictPhys()],
# [Law2_ScGeom_FrictPhys_CundallStrack()]
#),
FlowEngine(dead=1,label="flow"),#introduced as a dead engine for the
moment, see 2nd section
GlobalStiffnessTimeStepper(active=1,timeStepUpdateInterval=100,timestepSafetyCoefficient=0.8),
#PeriTriaxController(dynCell=True,mass=0.2,maxUnbalanced=0.01,relStressTol=0.02,goal=(-1e6,-1e6,-1e6),stressMask=3,globUpdate=5,maxStrainRate=(1e10,1e10,1e10),doneHook='triaxDone()',label='compressor'),
#NewtonIntegrator(damping=.6), ##Triaxonly
NewtonIntegrator(damping=.2),
]
factor = .5
O.dt = factor * PWaveTimeStep()
#B. Activate flow engine and set boundary conditions in order to get
permeability
#TW=TesselationWrapper()
#TW.triangulate()
#TW.computeVolumes()
flow.dead=0
flow.defTolerance=0.3
flow.meshUpdateInterval=200
flow.useSolver=3
flow.permeabilityFactor=1
flow.viscosity=10
flow.bndCondIsPressure=[0,0,1,1,0,0]
flow.bndCondValue=[0,0,1,0,0,0]
flow.boundaryUseMaxMin=[0,0,0,0,0,0]
O.dt=0.1e-3
O.dynDt=False
flow.emulateAction()
flow.saveVtk()
raw_input('FlowEngine saveVtk END')
--
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
