some question in triaxial test modeling
Hello everyone:
the following is my script:
*******
from yade import pack,utils,plot
from yade import utils
matP = FrictMat()
matP.density=
matP.young = 3e9
matP.poisson = 0.21
matP.frictionAngle = 0.65
matW = FrictMat()
matW.density = 1950#kg/m^3
matW.young = 1.2e9
matW.poisson = 0.33
matW.frictionAngle = 0.23
O.materials.
O.materials.
#cylinder in triaxial
#cyl=yade.
#O.bodies.
#O.bodies.erase()
#wall1=
#wall2=
#O.bodies.
#pred=pack.
#sp=SpherePack()
#sp=pack.
#spheres=
#former context is concern about how to build the model in cylinder,have some question.
wall1=utils.
wall2=utils.
wall3=utils.
wall4=utils.
wall5=utils.
O.bodies.
#pre=pack.
#spheres=
#O.bodies.
sp=pack.
sp.makeCloud(
sp.toSimulation()
#pre=pack.
#spheres=
#O.bodies.
sp=pack.
sp.makeCloud(
sp.toSimulation()
lpm=O.bodies[-1].id
O.bodies[
global inipos# global var, attention
inipos=
O.bodies.erase(lpm)
wall6= utils.wall(
O.bodies.
lpm=O.bodies[
inipos=
from yade import timing
yade.timing.reset() #Zero all timing data
O.engines=[
ForceResetter(),
InsertionSortC
InteractionLoop(
[Ig2_
[Ip2_
[Law2_
),
NewtonIntegrat
PyRunner(
PyRunner(
]
def Step():
from yade import utils, plot
global stepnum, inipos, deltpos, n
lpm=O.
platepos=
plateforce=
O.bodies[
O.bodies[
platepos2=
wallx1pos=
wallx2pos=
wally1pos=
wally2pos=
areaplate=
wallx1force=
wallx2force=
wally1force=
wally2force=
fmw=[abs(
G=(0.5*
wav=G*(fmw-500)
O.bodies[
O.bodies[
O.bodies[
O.bodies[
delvol=
plot.addData(
plot.saveDataT
if stepnum==1:
O.bodies[
if plateforce< 4500: return
stepnum=stepnum+1
elif stepnum==2:
O.bodies[
delconf=
if abs > 100: return
export.
print "wall6pos=
O.pause()
def Save():
global n
O.save(
n=n+1
#plot.plots=
plot.plots=
plot.plot()
#plot.saveDataT
from yade import qt
O.dt = 0.0005
V = qt.View()
V.screenSize = (550,450)
V.sceneRadius = 20
V.eyePosition = (0.7,0.5,0.1)
V.upVector = (0,0,1)
V.lookAt = (0.15,0.15,0.1)
O.run()
*******
However, when i tested it in yadedaily it showed some problems in below:
FATAL /data/trunk/
FATAL /data/trunk/
/data/trunk/
Besides that, i found that the walls in lateral are different from what i designed in line 35-41 after i added particles in the model.
So, i would be greatful if you can help me find out the causes of these problems.
Thanks a lot!
xue
Question information
- Language:
- English Edit question
- Status:
- Expired
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Last query:
- Last reply: