yade-users team mailing list archive

[rhaven <question669990@xxxxxxxxxxxxxxxxxxxxx>]

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.