shear band
How do I export a VTK file of a cut strip in a triaxial compression test? I have seen images of shear bands in the URL [1], but I don't know how to form such VTK files.
I applied the following code to generate the VTK file. But it does not show the perfect cut strip image as in [1]. Where am I going wrong, please?
[1]. https:/
#######
# -*- 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(
]
#def doExport():
#vtk.
O.run(100,True)
TW=TesselationW
TW.triangulate()
TW.computeVolumes()
TW.volume(10)
TW.setState(0)
O.run(20000,True)
TW.setState(1)
TW.defToVtk(
triax.strainRat
O.engines=
def stop():
if O.iter>40000 :
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()
#######
Question information
- Language:
- English Edit question
- Status:
- Answered
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Last query:
- Last reply:
Can you help with this problem?
Provide an answer of your own, or ask Xifan Li for more information if necessary.