Cubic Close Packing

Hello,

I have been looking at Kepler's close packing for my uni project. I am trying to see when spheres are arranged in cubic close packing and being inflated, how they would behave and if they would deform into rhombic dedocahedrons when they reach equilibrium.

http://mathworld.wolfram.com/CubicClosePacking.html

I was trying to model sphere lattice constrained within a boundary box. When inflated, they would not intersect with each other; they would stay in place; and would be malleable just enough to expand and fill in the gaps in between the spheres.

I started off with the help of this thread here(Thanks for those contributed!). As I understood, there was a bug in Kangaroo2. Solver can't handle more than one item plugged in. So I tried to understand David's Stasiuk's Script and adopted it with a few variations, please see gh file attached.

In the first 5 - I've used David Stasiuk's C# component-variable pressure (posted on June 9, 2015 at 12:25am): 'No. 4.5' being the most successful simulation so far(inflation value is kept very low so that they would not intersect);

although I realised I made some math mistake in setting the close packing grid.(could be checked by plugging voronoi3D to see if the area of the rhombic faces are regular)

No. 6-7 I tried with Kangaroo2 components.

After consulting my tutor(Andrei Jipa)'s help, I realised the following changes could be made:

- The definition posted by David on June 8, 2015 at 4:47pm with constant pressure would've worked better.

- Icosahedrons with WbCatmull(Quad divisions) would result in more even load distribution. With wbloop, vertices more concentrated at poles.

- Load in dir Z could be omitted.

Andrei has suggested to use lengths(line) in Kangaroo 2 as 'pressure' instead. And I am trying to improve the grid; and maybe try with David's constant pressure definition. I will keep you guys posted of the progress!

I am new to the parametric world, comments/advice very much appreciated! :)

Zhini