yade-users team mailing list archive

[raiden <question665793@xxxxxxxxxxxxxxxxxxxxx>]

Question #665793 on Yade changed:
https://answers.launchpad.net/yade/+question/665793

Description changed to:
Hello people,

I've been working on yade for a while and now I try a settlement
process. I set the particle parameters and subject them to gravitational
force and velocity to a box made from fecets. the problem is that
unbalanced force stuck at a value after hours and it goes on and on. Now
it's already more than 60 hours but still same. It is usually took 4-6
hours with same number of particles and radii but this time i changed
the young modulus. But i think it can't be the reason.


Any opinions or suggestions?

Thank you.

Edit: The code I use;

O.dt= 0.5*PWaveTimeStep()

print "Kieg eros"

for i in O.bodies:
  if i.isClumpMember:
    print i.shape.radius

def checkUnbalanced():

   if O.iter<5000: return

   if unbalancedForce()>0.01: return
   print utils.unbalancedForce()
   O.engines[2].lawDispatcher.functors[1].always_use_moment_law = True

   O.engines[2].physDispatcher.functors[1].setCohesionNow = True
   O.engines[2].physDispatcher.functors[1].setCohesionOnNewContacts = False
   O.engines[2].physDispatcher.functors[1].cohesionDisablesFriction=True
   O.engines=O.engines+[PyRunner(command='addPlotData()',iterPeriod=1)]
   # next time, do not call this function anymore, but the next one (unloadPlate) instead
   #O.run(500)
 # ----------------------------------------------------------------------------------------------------

   checker.command='stopUnloading()'

def stopUnloading():

      export.textClumps("/home/clums.txt")
      print "Save clumps"
      checker.command='addPlotData()'
      print utils.unbalancedForce()

      O.pause()

I changed YM to 5e6. When I tried before it last 4-6 hours even less. I
also decreased the timestep manually but it didn't help me to solve the
problem. I try to achive 0.01 unbalanced force, although when I try with
0.1 it also doesn't solve the problem.

-- 
You received this question notification because your team yade-users is
an answer contact for Yade.