New question #685486 on Yade:
https://answers.launchpad.net/yade/+question/685486
Hello,
I run the MWE and there are warning: infinite K1. I use regularly arranged
spheres, I think maybe it is the reason.
I use Yade 20191010-2428~bf906f7~bionic1.
Thanks for your any suggestion.
############MWE#############
from yade import pack
young=1e8
compFricDegree = 30
mn,mx=Vector3(0,0,0),Vector3(3,.05,.05)
O.materials.append(FrictMat(young=young,poisson=1,frictionAngle=radians(compFricDegree),density=2500,label='spheres'))
O.materials.append(FrictMat(young=young,poisson=1,frictionAngle=0,density=0,label='walls'))
walls=aabbWalls([mn,mx],thickness=0,material='walls')
wallIds=O.bodies.append(walls)
x = 0.005;y = 0.005;z = 0.005;r = 0.005
while x < 3:
while y < 0.05:
while z < 0.05:
O.bodies.append([sphere((x,y,z),r,material='spheres')])
z+=0.01
y+=0.01;z = 0.005
x+=0.01;y=0.005;z=0.005
newton=NewtonIntegrator(damping=0.2)
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Box_Aabb()]),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom(),Ig2_Box_Sphere_ScGeom()],
[Ip2_FrictMat_FrictMat_FrictPhys()],
[Law2_ScGeom_FrictPhys_CundallStrack()],label="iloop"
),
FlowEngine(label="flow"),
GlobalStiffnessTimeStepper(active=1,timeStepUpdateInterval=100,timestepSafetyCoefficient=0.8),
newton
]
flow.defTolerance=0.3
flow.meshUpdateInterval=200
flow.useSolver=3
flow.permeabilityFactor=1
flow.viscosity=0.001
flow.bndCondIsPressure=[1,1,0,0,0,0]
flow.bndCondValue=[0,1,0,0,0,0]
flow.boundaryUseMaxMin=[0,0,0,0,0,0]
#O.dt=0.1e-20
#O.dynDt=False
--
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
