yade-users team mailing list archive

[mrhappy <question689293@xxxxxxxxxxxxxxxxxxxxx>]

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

Hi,

I have some experience with YADE, but this is the first time I am using PFV. I want to experiment with DEM-PFV. I am aware that there is an example by the name Odometer Example. I have 4 spheres. Each sphere is connected to its neighboring sphere. I want to apply fluid pressure at the center of my model (0,0,0) and have the spheres spread outwards due to the pressure. I have made an attempt, but I have no clue what I am doing (the PFV part). I get an error:
0 : Vh==NULL!! id=0 Point=0.5 0.5 0 rad=0.5
1 : Vh==NULL!! id=1 Point=-0.5 0.5 0 rad=0.5
2 : Vh==NULL!! id=2 Point=-0.5 -0.5 0 rad=0.5
3 : Vh==NULL!! id=3 Point=0.5 -0.5 0 rad=0.5
AREA <= 0!!
AREA <= 0!!
AREA <= 0!!
AREA <= 0!!
AREA <= 0!!
AREA <= 0!!
AREA <= 0!!
AREA <= 0!!
I will attach my code. Any help on this would be great. 

Thank you soo much!


# MATERIAL PROPERTIES
intR=0.1 # allows near neighbour interaction (can be adjusted for every packing)
DENS=2500 # Density
YOUNG=1800
FRICT=7
ALPHA=0.1
TENS=100
COH=1000
iterMax = 1

O.materials.append(JCFpmMat(type=1,density=DENS,young=YOUNG,poisson=ALPHA,frictionAngle=radians(FRICT),tensileStrength=TENS,cohesion=COH,label='mat1'))

# PARTICLES
O.bodies.append([
 	sphere(center=( 0.5, 0.5, 0.0),radius=0.5,material='mat1'),
 	sphere(center=(-0.5, 0.5, 0.0),radius=0.5,material='mat1'),
	sphere(center=(-0.5,-0.5, 0.0),radius=0.5,material='mat1'),
	sphere(center=( 0.5,-0.5, 0.0),radius=0.5,material='mat1'),
])


# SIMULATION LOOP
O.engines=[
     	ForceResetter(),
     	InsertionSortCollider([Bo1_Sphere_Aabb()]),
  	InteractionLoop(
    		[Ig2_Sphere_Sphere_ScGeom(interactionDetectionFactor=intR)],
   		[Ip2_JCFpmMat_JCFpmMat_JCFpmPhys(cohesiveTresholdIteration=1)],
   		[Law2_ScGeom_JCFpmPhys_JointedCohesiveFrictionalPM()]
 	 ),
     	GlobalStiffnessTimeStepper(active=1,timeStepUpdateInterval=10,timestepSafetyCoefficient=0.5, defaultDt=utils.PWaveTimeStep()),
 	NewtonIntegrator(damping=0.1),
	FlowEngine(dead=0,label="flow"),
]

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]
flow.bndCondIsPressure=[0,0,0,1,0,0]
flow.bndCondValue=[0,0,0,0,0,0]
flow.updateTriangulation=True


i=createInteraction(0,1)
i=createInteraction(1,2)
i=createInteraction(2,3)
i=createInteraction(3,0)
i=createInteraction(0,2)
i=createInteraction(3,1)

#PLOTTING

# SAVE SIMULATION


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