Question #681232 on Yade changed: https://answers.launchpad.net/yade/+question/681232
gaoxuesong posted a new comment: >the problem is that with rMean=rRelFuzz, the radius may be arbitrarily small (you can print the minimal value), making the PWaveTimeStep very small, too. Actually i run an example case where all the particles have the same diameter of 100 um and the time step is set as 1e7. As soon as the simulation starts, the particles explode and flow out. > the last step is mass shrinking, where the mass of all the particles is gradually decreased to the 1/1000 of the used value. Does it make sense or do you have good advice? Sorry for confusing words used. The following is an example, ### a function to reduce the mass gradually redt = 0.5 ### how long for one mass reduction retime = 1.0 ### the time when mass reduction happens mrecof = 0.1 ## the mass reduction cofficient def mass_relaxation(): global redt global retime if O.time > retime: retime = retime + redt for eb in O.bodies: if isinstance(eb.shape, Sphere): eb.state.mass = eb.state.mass*mrecof print('change mass', O.time) Then the function of mass_relaxation() is invoked by pyRunner Engine periodically. Thanks. -- 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 : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp