time for reaching a static equilibrium
Hi,
I want to ask a question I have encountered, every time I change box size, it becomes very difficult to reach equilibrium. for example,
MWE:
from yade import pack
import math
nRead=utils.
young = 0.858e9,
poisson = 0.3,
num_spheres = 10000,
confiningS = 100e3,
rate = -0.02,
damp = 0.2,
unknownOk = True
)
from yade.params import table
mn,mx=Vector3(
young = table.young
poisson = table.poisson
compFricDegree = table.compFricD
finalFricDegree = table.finalFric
num_spheres = table.num_spheres
confiningS = table.confiningS
isoconfining = table.isoconfining
damp = table.damp
rate = table.rate
AxialStrainLimi
key='_cf-
cd = 1 # cd=1 - consolidate drained cd=0- undrained
porosity = 0.72
stabilityThreshold = 0.001
thick = 0.01
d_limit = 0.15 # (d_limit/dmin = 0.25)
roll_stiff = 0.01
targetPorosity = 0.43
O.materials.
young=young,
poisson=
density=2650,
frictionAng
label=
)
O.materials.
young=young/10,
poisson=
frictionAng
density=0,
label='walls')
)
# create walls around the packing
walls=aabbWalls
wallIds=
sp=pack.
sp.makeCloud(
O.bodies.
triax=TriaxialS
## TriaxialStressC
maxMultipli
finalMaxMul
thickness = thick,
stressMask = 7, # 1*1 + 1*2 + 1*3 = 7 means x,y,z are on
internalCom
)
newton = NewtonIntegrato
O.engines=[
ForceResett
InsertionSo
Interaction
),
GlobalStiffnes
triax,
TriaxialSta
newton,
]
triax.goal1=
timing.reset()
while 1:
O.run(1000, True)
unb=
print 'unbalanced force:',unb,' mean stress: ',triax.meanStress
if unb < stabilityThreshold and abs(-isoconfini
break
O.save(
print "### Initial state saved ###"
Thank you in advance!
Question information
- Language:
- English Edit question
- Status:
- Solved
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Solved by:
- Cloud
- Solved:
- Last query:
- Last reply: