yade-users team mailing list archive

[hjz <question291177@xxxxxxxxxxxxxxxxxxxxx>]

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.