yade-users team mailing list archive

[xjin <question404048@xxxxxxxxxxxxxxxxxxxxx>]

New question #404048 on Yade:
https://answers.launchpad.net/yade/+question/404048

I am trying to model a triaxial test.I write the script referred to the website:
https://github.com/yade/trunk/blob/master/examples/concrete/triax.py.

When I run my script,the plot-figure(whose codes is in the end) does not show the trend curves.The valuee in the figure are only zero as time goes. Could you give me some suggestion?
The scripe is shown below.
Thanks a lot!





from yade import pack
from yade import plot
from yade import os

densitytu=float(1600)
youngtu=40e6
potu=0.3
frictionAngletu= float(atan(0.7))
globals()['densitytu']=locals()['densitytu']
globals()['youngtu']=locals()['youngtu']
globals()['potu']=locals()['potu']
globals()['frictionAngletu']=locals()['frictionAngletu']

rParticle=float(2e-3)
rRelFuzz=float(1.75/8.75)*1.5
globals()['rParticle']=locals()['rParticle']
globals()['rRelFuzz']=locals()['rRelFuzz']

vel=1 
preStress=float(-60e3)
fwire=True
globals()['vel']=locals()['vel']
globals()['preStress']=locals()['preStress']
globals()['fwire']=locals()['fwire']

L=float(25e-3)
N=36
A=float(360/N)
H=float(100e-3)
globals()['L']=locals()['L']
globals()['N']=locals()['N']
globals()['A']=locals()['A']
globals()['H']=locals()['H']

runGunplot=False
runInGui=True
globals()['runGunplot']=locals()['runGunplot']
globals()['runInGui']=locals()['runInGui']

tu=O.materials.append(FrictMat(density=densitytu,young=youngtu,poisson=potu,frictionAngle=frictionAngletu,label='tuti'))
pred=pack.inCylinder((0,0,0),(0,0,H),L)
sp=SpherePack()
sp=pack.randomDensePack(pred,spheresInCell=3000,radius=rParticle,rRelFuzz=rRelFuzz,memoizeDb='/home/dj/yade/tmp/triaxTestOnCylinder.sqlite',returnSpherePack=True)
spheres=sp.toSimulation(color=(0,0.5,0),material=tu)

globals()['spheres']=locals()['spheres']
sl1=len(O.bodies)
globals()['sl1']=locals()['sl1']

for i in range(0,N): 
O.bodies.append(facet([(0,0,0),(L*cos(i*2*pi/N),L*sin(i*2*pi/N),0),(L*cos((i+1)* 2*pi/N),L*sin((i+1)*2*pi/N),0)],wire=False,material=tu))
	
sl2=len(O.bodies)
globals()['sl2']=locals()['sl2']

ti=[]
globals()['ti']=locals()['ti']
for i in xrange(sl1,sl2):
ti=ti+[i]

for xii in ti:
O.bodies[xii].state.blockedDOFs='XYZxyz'
O.bodies[xii].state.vel=(0,0,vel)

for i in range(0,N): 
O.bodies.append(facet([(L*cos(i*2*pi/N),L*sin(i*2*pi/N),H),(0,0,H),(L*cos((i+1)* 2*pi/N),L*sin((i+1)*2*pi/N),H)],wire=False,material=tu))

sl3=len(O.bodies)
globals()['sl3']=locals()['sl3']
tii1=[]
globals()['tii1']=locals()['tii1']
for iti in xrange(sl2,sl3):
tii1=tii1+[iti]

for i in tii1:
O.bodies[i].state.blockedDOFs='XYZxyz'
O.bodies[i].state.vel=(0,0,-1*vel)

nw=36
nh=20
globals()['nw']=locals()['nw']
globals()['nh']=locals()['nh']

rCyl2=L*2*0.5/cos(pi/float(nw))

facets=[]
globals()['facets']=locals()['facets']
for r in xrange(nw):
for h in xrange(nh):
v1=Vector3(rCyl2*cos(2*pi*(r+0)/float(nw)),rCyl2*sin(2*pi*(r+0)/float(nw)),H*(h+0)/float(nh))
v2=Vector3(rCyl2*cos(2*pi*(r+1)/float(nw)),rCyl2*sin(2*pi*(r+1)/float(nw)),H*(h+0)/float(nh))
v3=Vector3(rCyl2*cos(2*pi*(r+1)/float(nw)),rCyl2*sin(2*pi*(r+1)/float(nw)),H*(h+1)/float(nh))
v4=Vector3(rCyl2*cos(2*pi*(r+0)/float(nw)),rCyl2*sin(2*pi*(r+0)/float(nw)),H*(h+1)/float(nh))
f1=facet((v1,v2,v3),color=(0,0,0.1),material=tu,wire=fwire)
f2=facet((v1,v3,v4),color=(0,0,0.1),material=tu,wire=fwire)
facets.extend((f1,f2))

O.bodies.append(facets)

sl4=len(O.bodies)
globals()['sl4']=locals()['sl4']

mass=O.bodies[0].state.mass
for f in facets:
f.state.mass=(O.bodies[0].state.mass+O.bodies[1].state.mass+O.bodies[2].state.mass+O.bodies[3].state.mass+O.bodies[4].state.mass+O.bodies[5].state.mass+O.bodies[6].state.mass+O.bodies[7].state.mass+O.bodies[8].state.mass)/float(9)
f.state.blockedDOFs='XYZz'

def addForces(): 
  for f in facets: 
 	n=f.shape.normal
	a=f.shape.area
	O.forces.addF(f.id,preStress*a*n)

globals()['addForces']=locals()['addForces']
def addForces1(): 
  for stii in ti: 
 	n=O.bodies[stii].shape.normal
	a=O.bodies[stii].shape.area
	O.forces.addF(O.bodies[stii].id,-1*preStress*a*n)

globals()['addForces1']=locals()['addForces1']
def addForces2(): 
  for tstii in tii1: 
 	n=O.bodies[tstii].shape.normal
	a=O.bodies[tstii].shape.area
	O.forces.addF(O.bodies[tstii].id,preStress*a*n)

globals()['addForces2']=locals()['addForces2']

plot.plots = {'e':('s',),}
globals()['plot']=locals()['plot']

ybz=float(max([O.bodies[s].state.pos[2] for s in spheres]))-min([O.bodies[s].state.pos[2] for s in spheres])
globals()['ybz']=locals()['ybz']

def plotAddData(): 
f1= abs((sum(O.forces.f(b)[2] for b in tii1)-sum(O.forces.f(b)[2] for b in ti))/(pi*L*L))
yb1=float(max([O.bodies[s].state.pos[2] for s in spheres])-min([O.bodies[s].state.pos[2] for s in spheres]))
yb=ybz-yb1
s=f1
e=yb/ybz
plot.addData(
                       i=O.iter,
                      s=s,
                      e=e,
                           )

globals()['plotAddData']=locals()['plotAddData']
O.dt=0.5*utils.PWaveTimeStep() 
enlargeFactor=1.5

O.engines=[
ForceResetter(),
InsertionSortCollider([
Bo1_Sphere_Aabb(aabbEnlargeFactor=enlargeFactor,label='bo1s'),
Bo1_Facet_Aabb()
]),
InteractionLoop(
[
Ig2_Sphere_Sphere_ScGeom(interactionDetectionFactor=enlargeFactor,label='ss2d3dg'),
Ig2_Facet_Sphere_ScGeom()
],
[
Ip2_FrictMat_FrictMat_FrictPhys(),
],
[
Law2_ScGeom_FrictPhys_CundallStrack()
],
),
PyRunner(iterPeriod=1,command="addForces()"),
PyRunner(iterPeriod=1,command="addForces1()"),
PyRunner(iterPeriod=1,command="addForces2()"),
NewtonIntegrator(gravity=(0,0,0),damping=0.3),
PyRunner(iterPeriod=10,command='plotAddData()'),
]
O.step()

bo1s.aabbEnlargeFactor=ss2d3dg.interactionDetectionFactor=1.0

if runInGui:
    plot.plot()
    try:
         from yade import qt
         renderer=qt.Renderer()
    except:
         pass



-- 
You received this question notification because your team yade-users is
an answer contact for Yade.