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.
_______________________________________________
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
