yade-dev team mailing list archive

[Christian Jakob <897237@xxxxxxxxxxxxxxxxxx>]

I tried another script without facets (using aabbWalls) and without changing particle radii,
it crashes after 5 or 6 steps (when neverErase=True):

5468  FATAL yade.ThreadRunner /home/me/YADE/yade/core/ThreadRunner.cpp:31 run: Exception occured:                                                                                                      
Body #3084 has velocity==NaN!

I think the problem is the changes, that Bruno implemented with
neverErase. All contacts remain active with this option. But we need a
rupture distance, at which liquid bridges break and the contact gets
inactive. Something like

if (uN<0) && (uN>rupture_distance){
	if (neverErase) {
	  phys->kn = 0.0;
	  phys->ks = 0.0;
	  #...maybe more content here...
	  phys->shearForce = phys->normalForce = Vector3r::Zero();
	}
	else {scene->interactions->requestErase(id1,id2); return;}
}
if (uN<rupture_distance) {scene->interactions->requestErase(id1,id2); return;}

-- 
You received this bug notification because you are a member of Yade
developers, which is the registrant for Yade.
https://bugs.launchpad.net/bugs/897237

Title:
  bug in neverErase?!

Status in Yet Another Dynamic Engine:
  New

Bug description:
  If I set neverErase=True Yade crashes with following output:

  O.run(3000,True)                                                                                                                                                                                             
  RuntimeError: std::exception 

  and also:

  Traceback (most recent call last):                                                                                                                                                                               
    File "/home/me/YADE/YADE2970/lib/yade-bzr2970/py/yade/qt/__init__.py", line 192, in refreshEvent                                                                                                       15:01:13
      self.refreshValues()                                                                                                                                                                                 15:01:30
    File "/home/me/YADE/YADE2970/lib/yade-bzr2970/py/yade/qt/__init__.py", line 265, in refreshValues                                                                                                      15:01:34
      s=int(t); ms=int(t*1000)%1000; us=int(t*1000000)%1000; ns=int(t*1000000000)%1000                                                                                                                             
  ValueError: cannot convert float NaN to integer  

  
  You can reproduce the bug with the attached script (similar to the one i submitted with version 2960).
  First try to set Law2_ScGeom_MindlinPhys_Mindlin(neverErase=False), this will work.
  Then set neverErase=True and yade crashes.

To manage notifications about this bug go to:
https://bugs.launchpad.net/yade/+bug/897237/+subscriptions