WARNING! rmin>rmax when running TwoPhaseFlowEngine()
Hi There,
When I run the code below, I will be given the Error when arrives to > unsat.initializ
I guessed the grains out of the box, which I asked in the code to remove that grain, but the grains were not placed outside the box.
Also, using the lower timestep but the problem remains.
Code :
## encoding: utf-8
import matplotlib; matplotlib.
import matplotlib.pyplot as plt
import pylab
from yade import pack
import numpy as np
from numpy import *
import math
## Define Parameters
num_spheres=5000
compFricDegree=30
confiningS=-12500
RhoW=1000
RhoS=2700
youngModulus=75e6
poissonRatio=0.3
## creat a packing with a specific particle size distribution (PSD)
psdSizes=
psdCumm=
psdSizesArray=
psdSizesMeter=
sp=pack.
mn,mx=Vector3(
sp.makeCloud(
sp.psd(
sp.toSimulation
## create material, which will be used as default
O.materials.
O.materials.
## create walls around the packing
walls=aabbWalls
wallIds=
O.bodies.
triax=TriaxialS
internalCompac
goal1=confiningS,
goal2=confiningS,
goal3=confiningS,
label="triax"
)
newton=
O.engines=[
ForceResetter(),
InsertionSortC
InteractionLoop(
[Ig2_
[Ip2_
[Law2_
),
GlobalStiffnes
triax,
newton
]
O.dt=.1*
O.dynDt=False
while 1:
O.run(1000,True)
unb=unbalance
e=(triax.
print ' e: ',e,' Sigma: ',-triax.
if unb<0.15 and abs(confiningS-
break
################
## Start test ###
################
triax.stressMask=2
triax.goal1=
goalTop=
triax.goal2=goalTop
triax.wall_
## Instantiate a two-phase engine
unsat=TwoPhaseF
meanDiameter=
## set boundary conditions, the drainage is controlled by decreasing Water pressure and keeping air pressure constant
unsat.bndCondIs
unsat.bndCondVa
unsat.isPhaseTr
unsat.defTolera
unsat.meshUpdat
unsat.useSolver=3
unsat.permeabil
unsat.updateTri
unsat.surfaceTe
unsat.initializ
## Drainage
for pg in arange(
unsat.
unsat.invasion()
unsat.
for b in O.bodies:
O.forces.
while 1:
O.run(
unb=
Moisture=
print ' e: ',(triax.
if unb<0.5:
break
if Moisture<18.5:
break
-------
error:
WARNING! rmin>rmax. rmin=2.10136e-05 ,rmax=1e-10
WARNING! rmin>rmax. rmin=1.75633e-05 ,rmax=1e-10
WARNING! rmin>rmax. rmin=2.02999e-05 ,rmax=1e-10
WARNING! rmin>rmax. rmin=1.50292e-05 ,rmax=1e-10
WARNING! rmin>rmax. rmin=2.16551e-05 ,rmax=1e-10
WARNING! rmin>rmax. rmin=1.67435e-05 ,rmax=1e-10
WARNING! rmin>rmax. rmin=2.45035e-05 ,rmax=1e-10
WARNING! rmin>rmax. rmin=2.10136e-05 ,rmax=1e-10
WARNING! rmin>rmax. rmin=2.27934e-05 ,rmax=1e-10
WARNING! rmin>rmax. rmin=2.02674e-05 ,rmax=1e-10
WARNING! rmin>rmax. rmin=1.51404e-05 ,rmax=1e-10
-------
Question information
- Language:
- English Edit question
- Status:
- Solved
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Solved by:
- Chao Yuan
- Solved:
- Last query:
- Last reply: