yade-users team mailing list archive

[mrhappy <question689293@xxxxxxxxxxxxxxxxxxxxx>]

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

mrhappy gave more information on the question:
Hi guys,

I have played around with my model and I have added walls to my model. I
have read the following thesis [1].  Are there any other reading
material that would be usefull? With my new code (below), I get he
following error:

infinite K1!
infinite K1!
infinite K1!
infinite K1!
infinite K1!
infinite K1!
infinite K1!
infinite K1!

[1] A pore-scale coupled hydromechanical model for biphasic granular
media.

As always, thank you for your help.

My new code:
# 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'))

mn,mx=Vector3(-1,-1,-0.5),Vector3(1,1,0.5) # corners of the initial packing
walls=aabbWalls([mn,mx],thickness=0,material='mat1')
wallIds=O.bodies.append(walls)


# 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(6,7)
i=createInteraction(7,8)
i=createInteraction(8,9)
i=createInteraction(9,6)
i=createInteraction(6,8)
i=createInteraction(7,9)

#PLOTTING

# SAVE SIMULATION

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