GS did not converge in 20k iterations
Hello,
I used PeriodicFlowEngine to simulate Fluid-Structure Interaction, but I get ERROR :GS did not converge in 20k iterations (maybe because the reference pressure is 0?).I always faild to use PeriodicFlowEng
Thanks for your any suggestion.
#######
O.periodic=True
O.cell.
from yade import pack
num_spheres=1000# number of spheres
young=1e6
compFricDegree = 3 # initial contact friction during the confining phase
finalFricDegree = 30 # contact friction during the deviatoric loading
mn,mx=Vector3(
O.materials.
#O.materials.
#walls=
#wallIds=
sp=pack.
sp.makeCloud(
sp.toSimulation
triax=PeriTriax
mass=.2,
newton=
O.engines=[
ForceResetter(),
InsertionSortC
InteractionLoop(
[Ig2_
[Ip2_
[Law2_
),
PeriodicFlowEn
GlobalStiffnes
triax,
newton
]
while 1:
O.run(1000, True)
unb=unbalance
if unb<0.001 :
break
setContactFrict
## ______________ Oedometer section _________________
#B. Activate flow engine and set boundary conditions in order to get permeability
flow.dead=0
flow.defToleran
flow.meshUpdate
flow.permeabili
flow.viscosity=10
flow.bndCondIsP
flow.bndCondVal
flow.boundaryUs
O.dt=0.1e-3
O.dynDt=False
O.engines=
O.run(1,1)
Qin = flow.getBoundar
Qout = flow.getBoundar
permeability = abs(Qin)/1.e-4 #size is one, we compute K=V/∇H
print "Qin=",Qin," Qout=",Qout," permeability=
#C. now the oedometer test, drained at the top, impermeable at the bottom plate
flow.updateTria
newton.damping=0
#######
Question information
- Language:
- English Edit question
- Status:
- Solved
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Solved by:
- Robert Caulk
- Solved:
- Last query:
- Last reply: