problem about makeClumpCloud()

Hello,

>But it seems to don't work once the O.periodic=True is uncommented in the following codes because the box is a square box in 3d view.

Could you be more specific about what the problem/symptom is?

I would also like to kindly ask you to read through our How To Ask guide [1].

Cheers,

Robert

[1]https://www.yade-dem.org/wiki/Howtoask

Thank your reply, Robert Caulk.

I'd like to generate the initial packing of the cloud of clumps in a rectangular box whose length is 2 in the x-direction, height is 4 in the y-direction and width is 2 in the z-direction using the following codes.

*******************
O.periodic=True
sp=pack.SpherePack()
c1=pack.SpherePack([((0,0,0),.03333),((.03,0,0),.017),((0,.03,0),.017)])
sp.makeClumpCloud((0,0,0),(2,4,2),[c1],periodic=True,num=500)
sp.toSimulation()

*******************

But the problem I met is that the packing is not in a rectangular box, it is in a square box whose length is 2 in all three directions when I view it in GUI.

But the parallelepiped box I saw in the GUI of the Yade is a square box, how can I determine it is a parallelepiped measuring 2 x 4 x 2?

utils.aabbExtrema() is one way.

Thanks, Robert Caulk. I try the method or function you provide for me. Here are the codes:

*******************
from yade import utils
O.periodic=True
sp=pack.SpherePack()
c1=pack.SpherePack([((0,0,0),.03333),((.03,0,0),.017),((0,.03,0),.017)])
sp.makeClumpCloud((0,0,0),(2,4,2),[c1],periodic=True,num=500)
sp.toSimulation()
print(utils.aabbExtrema()[-1] - utils.aabbExtrema()[0]) # upper corner coordinates -
# lower corner coordinates
*******************

When I run the above codes serval times, it returns the different values of the difference between upper corner coordinates and lower corner coordinates each time.
For example, (1.63..., 3.95..., 1.80...), (1.44..., 3.39..., 1.75...) etc.. It is not equal to (2, 4, 2).

I want to know why the strange behavior occurs.

When I run your code on Yade 2018.02b it returns the correct dimensions.

Thanks, Robert Caulk. May I view the packing in Paraview through the post-process?

Thanks Robert Caulk, that solved my question.