VTK
hi, I m new on YADE. I have problems with vtk recorder. I m able to export the particles positions of my deposition and the facets that I have, and I am able to generate vtp and vtu files but I am not able to set the conditions in paraview such that I can get a variation of velocities of sphere during collapsing of the slope. Also while executing this script I am getting a statement like
"!joint cohesion total degradation! | iteration= 2000"
So does it have some reference why all the sphere are not falling down during collapsing of the slope above the persistent plane
from yade import plot, pack,utils,ymport
#### controling parameters
packing=
smoothContact=True
jointFrict=
jointDil=radians(0)
output=
maxIter=10000
#### Import of the sphere assembly
def sphereMat(): return JCFpmMat(
O.bodies.
## preprocessing to get dimensions of the packing
dim=utils.
dim=utils.
xinf=dim[0][0]
xsup=dim[1][0]
X=xsup-xinf
yinf=dim[0][1]
ysup=dim[1][1]
Y=ysup-yinf
zinf=dim[0][2]
zsup=dim[1][2]
Z=zsup-zinf
## preprocessing to get spheres dimensions
R=0
Rmax=0
numSpheres=0.
for o in O.bodies:
if isinstance(
numSpheres+=1
R+=o.
if o.shape.
Rmax=
Rmean=R/numSpheres
#### Identification of the spheres on joint (some DIY here!) -> work to do on import function textExt to directly load material properties from the ascii file
inFile=
for line in inFile:
if '#' in line : continue
id = int(line.
onJ = int(line.
nj = int(line.
j11 = float(line.
j12 = float(line.
j13 = float(line.
j21 = float(line.
j22 = float(line.
j23 = float(line.
j31 = float(line.
j32 = float(line.
j33 = float(line.
O.bodies[
O.bodies[
O.bodies[
O.bodies[
O.bodies[
inFile.close
#### Boundary conditions
e=2*Rmean
Xmax=0
Ymax=0
baseBodies=[]
for o in O.bodies:
if isinstance(
o.
## to fix boundary particles on ground
if o.state.
o.state.
baseBodies.
o.shape.
## to identify indicator on top
if o.state.
refPoint=o.id
p0=o.state.pos[1]
baseBodies=
O.bodies[
#### Engines definition
interactionRadi
O.engines=[
ForceResetter(),
InsertionSortC
InteractionLoop(
[Ig2_
[Ip2_
[Law2_
),
GlobalStiffnes
PyRunner(
VTKRecorder(
PyRunner(
NewtonIntegrat
]
#### dataCollector
plot.plots=
def dataCollector():
R=O.bodies[
plot.addData(
plot.saveDataT
#### joint strength degradation
stableIter=2000
stableVel=0.001
degrade=True
def jointStrengthDe
global degrade
if degrade and O.iter>=stableIter and abs(O.bodies[
print '!joint cohesion total degradation!', ' | iteration=', O.iter
degrade=False
for i in O.interactions:
if i.phys.isOnJoint :
if i.phys.isCohesive:
i.phys.
i.phys.FnMax=0.
i.phys.FsMax=0.
#### YADE windows
from yade import qt
v=qt.Controller()
v=qt.View()
O.dt=0.
#### set cohesive links with interaction radius>=1
O.step();
#### initializes now the interaction detection factor to strictly 1
ss2d3dg.
is2aabb.
#### if you want to avoid contact detection (Lattice like)
#O.engines=
#### RUN!!!
O.dt=-0.
O.run(maxIter)
plot.plot()
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:
This question was reopened
- by shiv