microstrain field
Hi all,
I am trying to get the microstrain field. I did as the document said[1]. I want to use PyRunner command to export a series of .vtk files. It seems that I can just get one VTK file. is there someone can help me to check the ### get the strain field part ### in the MWE.
from yade import pack, plot
from yade import export
from yade import utils
compFricDegree=45.0
finalFricDegree
rate=-15
damp=0.7
target_strain = 0.15
stabilityThresh
young=4e8
mn,mx=Vector3(
O.materials.
O.materials.
walls=aabbWalls
wallIds=
sp=pack.
psdSizes=
psdCumm=
sp.makeCloud(
sp.toSimulation
O.dt=.5*
triax=TriaxialS
maxMultiplier=
finalMaxMultip
thickness=0,
stressMask=7,
internalCompac
)
newton=
O.engines=[
ForceResetter(),
InsertionSortC
InteractionLoop(
[Ig2_
[Ip2_
[Law2_
),
GlobalStiffnes
triax,
newton,
]
triax.goal1=
while 1:
O.run(3000,True)
unb=unbalanced
print 'unbalanced force: ',unb,'mean stress: ',triax.meanStress
if unb<stabilityTh
break
triax.wall_
triax.wall_
triax.wall_
triax.wall_
triax.wall_
triax.wall_
setContactFrict
re11=-triax.
re22=-triax.
re33=-triax.
triax.stressMask = 5
newton=
triax.goal2=rate
triax.strainRat
triax.strainRat
def stop_loading():
if -triax.
O.pause()
#######
TW=TesselationW
TW.triangulate()
TW.computeVolumes()
TW.setState(0)
#step = 0
#while -triax.
#step = step+500
#O.run(400,True)
#TW.setState(1)
#TW.defToVtk(
O.run(100,True)
TW.setState(1)
def strain_export():
#O.run(100,True)
TW.setState(1)
TW.defToVtk(
#######
O.engines=[
ForceResetter(),
InsertionSortC
InteractionLoop(
[Ig2_
[Ip2_
[Law2_
),
GlobalStiffnes
triax,
newton,
PyRunner(
PyRunner(
]
O.run(50000,True)
I found that there is a question is almost the same as my question[2], they use a while function to solve this problem. But I want to know how can we solve this problem by using pyrunner function. thanks.
[1]https:/
[2]https:/
best,
yong
Question information
- Language:
- English Edit question
- Status:
- Solved
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Solved by:
- Bruno Chareyre
- Solved:
- Last query:
- Last reply: