New question #466297 on Yade:
https://answers.launchpad.net/yade/+question/466297
Hi,
I'm new to Yade and would like to use the PeriodicFlowEngine. I try to simulate
a bi-periodic packing. The packing is supported by a bottom plane and another
plane is located at the top. The impermeable top plane should move downwards
until a certain maximum load is reached. Furthermore, the top plane moves
sidewards to apply a shearing force onto the packing.
I get the 'periodicity is broken' error and the triangulation cannot be build:
'buildTriangulation: problem here, flow.errorCode>0'.
I'm using Yade 2016 on Ubuntu 16.04.
Does anyone know what I did wrong?
Thanks,
Matari
# -*- coding: utf-8 -*-
from yade import pack, plot
O.periodic=True
radius=0.3
rRelFuzz=0.1
numParticles=2000
numCohParticles=100
porosity=0.3
height=40*radius
length=25*radius
width=25*radius
youngSpheres=1e6
poissonSpheres=0.5
frictDegreeSpheres=30
densitySpheres=2000
youngWalls=1e6
poissonWalls=0.5
frictDegreeWalls=30
densityWalls=7000
thickness=0.01
maxLoad=1000.*length*width
coh=1000000.
gap=0.01*length
O.cell.hSize=Matrix3(length,0.,0.,
0.,3*height,0.,
0.,0.,width)
O.materials.append(CohFrictMat(young=youngSpheres,poisson=poissonSpheres,frictionAngle=radians(frictDegreeSpheres),density=densitySpheres,isCohesive=False,label='spheres'))
O.materials.append(CohFrictMat(young=youngSpheres,poisson=poissonSpheres,frictionAngle=radians(frictDegreeSpheres),density=densitySpheres,fragile=False,normalCohesion=coh,shearCohesion=poissonSpheres*coh,label='cohSpheres'))
O.materials.append(CohFrictMat(young=youngWalls,poisson=poissonWalls,frictionAngle=radians(frictDegreeWalls),density=densityWalls,fragile=False,normalCohesion=coh,shearCohesion=poissonWalls*coh,label='walls'))
sp1=pack.SpherePack()
sp1.makeCloud((0.,height+1.1*radius,0.),(length,2.*height,width),radius,rRelFuzz,numParticles,periodic=True,porosity=porosity)
O.bodies.append([utils.sphere(s[0],s[1],color=(1,1,1),material='spheres') for s
in sp1])
pred=pack.inAlignedBox(Vector3(0.,2.0*height,0.),Vector3(length,2.1*height,width))
sp2=pack.regularHexa(pred,radius=radius,gap=gap,material='cohSpheres')
O.bodies.append(sp2)
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Box_Aabb()],allowBiggerThanPeriod=True),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom6D(),Ig2_Box_Sphere_ScGeom6D()],
[Ip2_CohFrictMat_CohFrictMat_CohFrictPhys(setCohesionNow=False,setCohesionOnNewContacts=True)],
[Law2_ScGeom6D_CohFrictPhys_CohesionMoment()]
),
GlobalStiffnessTimeStepper(active=1,timeStepUpdateInterval=100,timestepSafetyCoefficient=0.8),
NewtonIntegrator(gravity=(0,0,0),damping=0.2),
PeriodicFlowEngine(dead=1,label="flow"),
PyRunner(command='loading()',iterPeriod=100,label='stop')
]
bottomWall = utils.box(center=(length/2.0,height-thickness/2.0,width/2.0),
extents=(length*1000.0,thickness/2.0,width*1000.0)
,fixed=True,wire=False,material='walls')
topWall = utils.box(center=(length/2.0,2.1*height-thickness/2.0,width/2.0),
extents=(length*1000.0,thickness/2.0,width*1000.0)
,fixed=True,wire=False,material='walls')
O.bodies.append([bottomWall,topWall])
def loading():
topWall.state.vel=(1,0,0)
if abs(O.forces.f(topWall.id)[1])<maxLoad:
topWall.state.vel=(1,-0.1,0)
O.engines=O.engines+[PyRunner(command='addPlotData()',iterPeriod=100)]
def addPlotData():
force=O.forces.f(topWall.id)[1]
plot.addData(force=force,i=O.iter)
plot.plots={'i':('force',)}
plot.plot()
flow.dead=0
flow.defTolerance=0.3
flow.meshUpdateInterval=200
flow.useSolver=3
flow.permeabilityFactor=1
flow.eps=0.04*radius
flow.viscosity=1
flow.bndCondIsPressure=[1,1,0,1,1,1]
flow.bndCondValue=[0,0,0,0,0,0]
O.dt=0.1e-3
O.dynDt=True
--
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
