yade-users team mailing list archive

[onyourself <question691171@xxxxxxxxxxxxxxxxxxxxx>]

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

    Status: Answered => Open

onyourself is still having a problem:
>my randomOrientation most likely is not uniform, just "some randomness"

that works for me too(i didn't figure out way to achieve ‘random
directions and uniform distribution’ so…)

>cylinders (boxes) replaces spheres, so no sphere is generated within
spherical space (not cubic).

replaces?  i don't get it…the code below makes cloud of 40 spheres and
30 spheres first(not adding to simulation),then when i<40, add every box
to bodies,(where are those 40 spheres?),then add 30 spheres to bodies?
Why not adding all 70 spheres to bodies?

>>for i,(pos,radius) in enumerate(sp):
>>O.bodies.append(sphere(pos,radius))


#########
def randomOrientation(): # most likely not very well random, but for this example it is ok
   from random import random as r
   return Quaternion((r(),r(),r()),2*pi*r())

# create loose spherical packing using makeCloud
sp = yade.pack.SpherePack()
mi,ma = (0,0,0),(20,20,20)
nCyls,nSphs = 40,30
cylinder(length=0.012,radius=0.0008,color=[0,0,1],)
sp.makeCloud(mi,ma,rMean=2,rRelFuzz=.5,num=nCyls) # makeCloud for cylinders
sp.makeCloud(mi,ma,rMean=1,rRelFuzz=.5,num=nSphs) # makeCloud for spheres
for i,(pos,radius) in enumerate(sp):
   if i < nCyls: # add cylinder
      ori = randomOrientation()
      extents = (radius,.1*radius,.1*radius)
      O.bodies.append(box(pos,extents,ori))
   else: # add sphere
      O.bodies.append(sphere(pos,radius))
##########

1)i need to simulate fiber so i changed a bit then it's not working, cylinders are out of space requested before. 
2)i set orientation=Quaternion ((r(),r(),r()),r()) but it turns out cylinders' orientation are not random.

my code is below.

#########

from random import random as r

def randomOrientation(): # most likely not very well random, but for this example it is ok
   from random import random as r
   return Quaternion((r(),r(),r()),r())

# create loose spherical packing using makeCloud
sp = yade.pack.SpherePack()
mi,ma = (0,0,0),(0.08,0.08,0.08)
nCyls,nSphs = 1000,100000
sp.makeCloud(mi,ma,rMean=0.012,rRelFuzz=.003,num=nCyls) # makeCloud for cylinders
sp.makeCloud(mi,ma,rMean=0.0011,rRelFuzz=.0005,num=nSphs) # makeCloud for spheres
for i,(pos,radius) in enumerate(sp):
    if i < nCyls: # add cylinder
      O.bodies.append(geom.facetCylinder(center=(r(),r(),r()),height=0.012,radius=0.0015,orientation=Quaternion(),segmentsNumber=10,wallMask=7))
    else: # add sphere
      #O.bodies.append(sphere())    
      O.bodies.append(sphere(pos,radius))    

##########

thank u soooo much
xxxe

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