Triaxial compression test displacement field
I tried to make the displacement field of the second stage of triaxial compression. I have this answer in the following question, but I have not implemented this method. Here is my code, please tell me where the error is.
https:/
from yade import pack,plot,qt,export
nRead=readParam
)
from yade.params import table
num_spheres=
key=table.key
targetPorosity = 0.42
compFricDegree = table.compFricD
finalFricDegree = 30
rate=-0.02
damp=0.2
stabilityThresh
young=5e8
mn,mx=Vector3(
weiya=100000
poisson=0.1
O.materials.
O.materials.
walls=aabbWalls
wallIds=
psdSizes,psdCumm =[0.19,
sp=pack.
sp.makeCloud(
O.bodies.
triax=TriaxialS
thickness = 0,
stressMask = 7,
)
newton=
O.engines=[
),
triax,
newton,
]
yade.qt.
triax.goal1=
while 1:
O.run(1000, True)
unb=unbalance
print 'unbalanced force:',unb,' mean stress: ',triax.meanStress
if unb<stabilityTh
break
print "### Isotropic state saved ###"
import sys
while triax.porosity>
compFricDegree = 0.95*compFricDegree
setContactFric
print "\r Friction: ",compFricDegree," porosity:
sys.stdout.flush()
O.run(500,1)
print 'Box Volume: ', triax.boxVolume
print "### Compacted state saved ###"
triax.internalC
setContactFrict
triax.stressMask = 5
triax.goal1=-weiya
triax.goal2=rate
triax.goal3=-weiya
newton.damping=0.1
O.saveTmp()
vtk = export.
vtk.exportSpher
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: