Re: [Yade-users] [Question #631227]: Porosity of a packing with overlapped particles
Question #631227 on Yade changed: https://answers.launchpad.net/yade/+question/631227 Status: Open => Answered Klaus Thoeni proposed the following answer: Well, are you sure you selected the right volume? Make sure you use the methods the right way and set the right parameters. Without seeing your script I can't give any further advice. Klaus -- 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
Re: [Yade-users] [Question #631227]: Porosity of a packing with overlapped particles
Question #631227 on Yade changed: https://answers.launchpad.net/yade/+question/631227 Status: Answered => Open Huihuang Xia is still having a problem: Hi Klaus, Thanks for your patience. I have selected two diverse cubic sub-volumes in my irregular packing, but they got a porosity of 1.0 Huihuang Xia -- 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
Re: [Yade-users] [Question #631227]: Porosity of a packing with overlapped particles
Question #631227 on Yade changed: https://answers.launchpad.net/yade/+question/631227 Status: Answered => Open Huihuang Xia is still having a problem: Thanks Klaus, I have tried your method, this approach works really good. Here is another question, Can I use this method to calculate the porosity of an irregular packing? Choosing a sub-volume from a irregular packing can give a result as precise as the whole volume? FW your reply, regards, Huihuang -- 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
Re: [Yade-users] [Question #631227]: Porosity of a packing with overlapped particles
Question #631227 on Yade changed: https://answers.launchpad.net/yade/+question/631227 Status: Open => Answered Klaus Thoeni proposed the following answer: Hi, you should use utils.voxelPorosity [1] instead of utils.porosity(), this works fine for overlapping particles. HTH Klaus [1] https://yade- dem.org/doc/yade.utils.html?highlight=voxelporosity#yade._utils.voxelPorosity -- 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
[Yade-users] [Question #631227]: Porosity of a packing with overlapped particles
New question #631227 on Yade: https://answers.launchpad.net/yade/+question/631227 Hi, I defined a function to calculate the porosity of a packing with lots of overlapped particles, because the build-in function "utils.porosity()" always got a negative one. This question (https://answers.launchpad.net/yade/+question/404363) gave me some useful hints. My own function has a serious error, that is my function always got a same porosity whatever the growFactor is. Here is my code segment: from yade import utils from yade import pack from yade import qt mat1 = O.materials.append(FrictMat(young=5.0e7,poisson=0.25,density=2640.0,frictionAngle=radians(18))) pred = pack.inAlignedBox((0,0,0),(0.01,0.01,0.01)) dim = pred.dim() packing = regularHexa(pred, radius=0.0005,gap=0,material=mat1) O.bodies.append(packing) # grow particles in order to get lower porosity growFactor = 1.5 utils.growParticles(growFactor) def overlappedPorosity(totVol,density): totVol = dim[0]*dim[1]*dim[2] mass = sum([b.state.mass for b in O.bodies])*pow(growFactor,-3) return (totVol - mass/density)/totVol print 'The porosity of this packing is:',overlappedPorosity(dim[0]*dim[1]*dim[2],O.materials[0].density) # 3D view and controller try: qt.Controller() qt.View() except: pass In my function, ''mass = sum([b.state.mass for b in O.bodies])*pow(growFactor,-3)'', because after the packing was generated, the radius of each particle in packing was magnified by growFactor using growParticles. Any useful advice will be highly appreciated. Thanks in advance, Huihuang Xia -- 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