Question #291177 on Yade changed:
https://answers.launchpad.net/yade/+question/291177
hjz posted a new comment:
Sorry, I don't put my code on.
This is the materials:
O.materials.append(FrictMat(young=5e9,poisson=0.3,frictionAngle=radians(compFricDegree),density=2600,label='spheres'))
O.materials.append(FrictMat(young=5e10,poisson=0.49,frictionAngle=0,density=0,label='walls'))
I use the method of compression to make a sample:
triax.goal1=triax.goal2=triax.goal3=-10000
while 1:
O.run(1000, True)
##the global unbalanced force on dynamic bodies, thus excluding boundaries,
which are not at equilibrium
unb=unbalancedForce()
print 'unbalanced force:',unb,' mean stress: ',triax.meanStress
if unb<0.05 and abs(-10000-triax.meanStress)/10000<0.1:
break
import sys
while triax.porosity>0.33:
# we decrease friction value and apply it to all the bodies and contacts
compFricDegree = 0.95*compFricDegree
setContactFriction(radians(compFricDegree))
print "\r Friction:",compFricDegree,"porosity:",triax.porosity
sys.stdout.flush()
# while we run steps, triax will tend to grow particles as the packing
# keeps shrinking as a consequence of decreasing friction. Consequently
# porosity will decrease
O.run(500,1)
--
You received this question notification because your team yade-users is
an answer contact for Yade.
_______________________________________________
Mailing list: https://launchpad.net/~yade-users
Post to : [email protected]
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
