openscad/dotSCAD
1
0
Fork 0
mirror of https://github.com/JustinSDK/dotSCAD.git synced 2025-08-22 06:13:10 +02:00
Code Issues Releases Wiki Activity

use circle_packing

Browse Source
This commit is contained in:
Justin Lin
2022-05-31 14:36:23 +08:00
parent 8c3a7d7f7a
commit b80ea25688
1 changed files with 15 additions and 54 deletions
Download Patch File Download Diff File Expand all files Collapse all files

69
examples/packing_circles.scad
View File

@@ -1,66 +1,27 @@
use <util/rand.scad>;
use <experimental/circle_packing.scad>;
size = [500, 250];
min_radius = 5;
max_radius = 50;
total_circles = 100;
size = [200, 100];
density = 2;
min_radius = 1;
point_numbers = 100;
function _packing_circles_out_of_range(size, c) =
c[0] + c[2] >= size[0] ||
c[0] - c[2] <= 0 ||
c[1] + c[2] >= size[1] ||
c[1] - c[2] <= 0;
points = [
for(i = [0:point_numbers - 1])
[rand(0, size.x), rand(0, size.y)]
];
function _packing_circles_overlapping(circles, c, i = 0) =
i == len(circles) ? false :
let(
x = c[0] - circles[i][0],
y = c[1] - circles[i][1],
a = c[2] + circles[i][2],
collision = a ^ 2 >= x ^2 + y ^ 2
)
collision || _packing_circles_overlapping(circles, c, i + 1);
function _packing_circles_packable(size, circles, c) =
!_packing_circles_overlapping(circles, c) &&
!_packing_circles_out_of_range(size, c);
function _packing_circles_new_min_circle(size, min_radius, attempts, circles, i = 0) =
i == attempts ? [] :
let(c = [rand() * size[0], rand() * size[1], min_radius])
_packing_circles_packable(size, circles, c) ? c :
_packing_circles_new_min_circle(size, min_radius, attempts, circles, i + 1);
function _packing_circles_increase_radius(size, circles, c, max_radius) =
c[2] == max_radius || !_packing_circles_packable(size, circles, c) ?
[c[0], c[1], c[2] - 1] :
_packing_circles_increase_radius(size, circles, [c[0], c[1], c[2] + 1], max_radius);
function _packing_circles_new_circle(size, min_radius, max_radius, attempts, circles) =
let(c = _packing_circles_new_min_circle(size, min_radius, attempts, circles))
c == [] ? [] : _packing_circles_increase_radius(size, circles, c, max_radius);
function _packing_circles(size, min_radius, max_radius, total_circles, attempts, circles = [], i = 0) =
i == total_circles ? circles :
let(c = _packing_circles_new_circle(size, min_radius, max_radius, attempts, circles))
c == [] ? _packing_circles(size, min_radius, max_radius, total_circles, attempts, circles) :
_packing_circles(size, min_radius, max_radius, total_circles, attempts, [each circles, c], i + 1);
function packing_circles(size, min_radius, max_radius, total_circles, attempts = 100) =
_packing_circles(is_num(size) ? [size, size] : size, min_radius, max_radius, total_circles, attempts);
circles = packing_circles(size, min_radius, max_radius, total_circles);
mr = max([for(c = circles) c[2]]);
circles = circle_packing(points, density, min_radius);
mr = max([for(c = circles) c[1]]);
translate([0, 0, mr])
for(c = circles) {
translate([c[0], c[1]])
sphere(c[2], $fn = 48);
translate(c[0])
sphere(c[1], $fn = 48);
}
for(c = circles) {
translate([c[0], c[1]])
translate(c[0])
linear_extrude(mr)
circle(c[2]/ 3, $fn = 48);
circle(c[1]/ 3, $fn = 48);
}
linear_extrude(1) square(size);