yade-dev team mailing list archive
-
yade-dev team
-
Mailing list archive
-
Message #02998
Re: [Bug 505523] Re: Velocity Bin
Thanks Vaclav. There are too many files to upload. I have a modified Triaxial Test, Triaxial Compression Engine, TriaxialStressController, because I have used a new contact law. It is from the CohesiveFrictionalLaw but without creep and cohesion. So it is a cohesionless moment resistance law. I'm using one of the old releases. I am trying to simulate a paper by Plassiard et al. (2009) "A spherical discrete element model: calibration procedure and incremental response". They used the Weibull distribution. I have modified the SpherePack file to calculate the radius, according to their paper. So I have a group of files..
Real rMin = rMax/MaxMinRatio;
Real r = rMin + (rMax-rMin)* pow( (-Alpha*Mathr::Log(1-rnd())) ,1/Gamma);
Some particles are really small. But everytime i run a simulation, it gives me a different timestep. And sometimes it is very large, because I have very large particles in the simulation. But the calculation of the timestep should be based on the smallest particle.
Yours,
CWBoon
> Date: Sun, 10 Jan 2010 17:15:14 +0000
> From: eudoxos@xxxxxxxx
> To: yade-dev@xxxxxxxxxxxxxxxxxxx
> Subject: [Yade-dev] [Bug 505523] Re: Velocity Bin
>
> Attach your simulation. Likely it is problem on your side with
> simulation going nuts numerically.
>
> One of common causes of NaN velocity is having body with zero mass
> dynamic (such as a facet), which gives inf acceleration etc.
>
> BTW if you don't know what NaN is, check
> http://en.wikipedia.org/wiki/NaN .
>
> --
> Velocity Bin
> https://bugs.launchpad.net/bugs/505523
> You received this bug notification because you are a member of Yade
> developers, which is subscribed to Yade.
>
> Status in Yet Another Dynamic Engine: New
>
> Bug description:
> I have the following error:
>
> 614459 FATAL yade.VelocityBins /home/chiab/yadeNew/Branch3/pkg/common/DataClass/VelocityBins.cpp:77 setBins: Body #1998, velSq=nan not put in any bin!?
> 614459 [0x7f16ef013950] FATAL yade.VelocityBins null - /home/chiab/yadeNew/Branch3/pkg/common/DataClass/VelocityBins.cpp:77 setBins: Body #1998, velSq=nan not put in any bin!?
>
> Also, global stiffness timestepper seems to give me an exceedingly large step. Sometimes 1e(positive)8!
>
>
> CWBoon
>
>
>
> _______________________________________________
> Mailing list: https://launchpad.net/~yade-dev
> Post to : yade-dev@xxxxxxxxxxxxxxxxxxx
> Unsubscribe : https://launchpad.net/~yade-dev
> More help : https://help.launchpad.net/ListHelp
_________________________________________________________________
Windows 7: Find the right PC for you. Learn more.
http://windows.microsoft.com/shop
Follow ups
References