no vtk file exported
I am simulating a two-dimensional biaxial compression experiment and I would like to use paraview for post-image optimisation (e.g. to form shear bands), but I am not sure why I am always unable to generate the appropriate VTK files. Here is my code.
*******
# -*- encoding=utf-8 -*-
# the script demonstrates a simple case of triaxial simulation using TriaxialCompres
from yade import pack
sp=pack.
## corners of the initial packing
mn,mx=Vector3(
## box between mn and mx, avg radius ± ½(20%), 2k spheres
sp.makeCloud(
## create material #0, which will be used as default
O.materials.
O.materials.
## copy spheres from the packing into the scene
## use default material, don't care about that for now
O.bodies.
## create walls around the packing
walls=aabbWalls
wallIds=
triax=TriaxialC
wall_bottom_
wall_top_
wall_left_
wall_right_
wall_back_
wall_front_
internalCompac
## define the rest of triax params here
## see in pkg/dem/
## which are assigned in the c++ preprocessor actually
sigmaIsoCompac
sigmaLateralCo
max_vel=10,
label="triax"
)
O.engines=[
ForceResetter(),
InsertionSortC
InteractionLoop(
[Ig2_
[Ip2_
[Law2_
),
GlobalStiffnes
triax,
TriaxialStateR
NewtonIntegrat
VTKRecorder(
]
triax.strainRat
O.engines=
def stop():
if O.iter>50000 :
O.pause()
from yade import plot
def history():
plot.addData(
s11=-
s22=-
i=O.iter)
plot.plots={'e11': ('s11',),'e22': ('s22',)}
O.saveTmp()
plot.plot()
*******
Any reply would be helpful.
Thank you.
Question information
- Language:
- English Edit question
- Status:
- Solved
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Solved by:
- Karol Brzezinski
- Solved:
- Last query:
- Last reply: