X-Git-Url: http://www.chiark.greenend.org.uk/ucgi/~ianmdlvl/git?p=reprap-play.git;a=blobdiff_plain;f=poster-tube-lid-parametric.scad.pl;h=64d4d00ebe9f5b472789f5411155e6457b95d8ec;hp=9fa0b2e1cdbd5c02e8c9ad229d36c4e812a4046a;hb=1e5609f5812a3103813be78c322b383c8a13d35b;hpb=b28141e959f59d188e7003849d8923fa951d66a4;ds=sidebyside diff --git a/poster-tube-lid-parametric.scad.pl b/poster-tube-lid-parametric.scad.pl index 9fa0b2e..64d4d00 100755 --- a/poster-tube-lid-parametric.scad.pl +++ b/poster-tube-lid-parametric.scad.pl @@ -16,20 +16,20 @@ my $slope_angle = 45 * TAU/360; my $jcurverad = 5; my $tall = 50; -my $lin_len = 5; -my $sine_size = 10; -my $sine_angle = TAU/8; +my $lin_len = 2; +my $sine_size = 5; +my $sine_angle = 1.20 * TAU/8; my $ballend_xr = $thick/2; +my $skew_slope = 0.7; + my @i_sections = qw(ball0 -6 sine0 -10 - - - lin0 2 - circle 20 + lin0 -2 + circle 40 lin1 2 sine1 10 - - ball2 6 - ); @@ -47,11 +47,13 @@ my @j_sections = qw(lin0 2 sub point ($$$$) { my ($ip,$it, $jp,$jt) = @_; - my ($i_offset, $i_outward); + my ($i_offset); - $i_outward = V( 0, - ($ip =~ m/0$/ ? -1 : +1), - 0 ); + my $i_outward = V( 0, + ($ip =~ m/0$/ ? -1 : +1), + 0 ); + + my $i_j_y_angle = 0; my $i_thickscale = 1.0; my $sine_len = $sine_size * sin($sine_angle); @@ -68,22 +70,28 @@ sub point ($$$$) { 0 ); } elsif ($ip =~ m/^sine[01]$/) { my $angle = $it * $sine_angle; - $i_offset = V( -$sine_size * sin($angle) -$lin_len, + $i_offset = V( -$lin_len -$sine_size * sin($angle), 0, +$sine_size * (1 - cos($angle)) ); + $i_j_y_angle = $angle; } elsif ($ip =~ m/^ball[02]$/) { + $i_j_y_angle = $sine_angle; my $angle = $it * TAU/4; - $i_offset = V( -$lin_len -$sine_len - -sin($angle) * $ballend_xr, + my $dx = sin($angle) * $ballend_xr; + $i_offset = V( -$lin_len -$sine_len - $dx * cos($sine_angle), 0, - +$sine_height + +$sine_height + $dx * sin($sine_angle) ); $i_thickscale = cos($angle); } else { die "$ip ?"; } + my $i_j_y_vect = V( sin($i_j_y_angle), + 0, + cos($i_j_y_angle )); + my $j_plus_th = $jp =~ m/2$/ ? $thick : 0; my $i_thick = $thick * $i_thickscale; @@ -120,11 +128,18 @@ sub point ($$$$) { die "$jp ?"; } + $j_y -= $j_qrc_y; + + if ($j_y > 0) { + $i_j_y_vect = V(0,0,1); + } + # print STDERR "@_ $j_x $j_y $i_offset $i_outward\n"; return $i_offset + $j_x * $i_outward + - V(0,0,1) * $j_y + + $i_j_y_vect * $j_y + + V(0,0,1) * $j_qrc_y + V(0,0,-$tall) ; } @@ -154,20 +169,29 @@ our @points; our %point_indices; our @triangles; +my @ipts; +my @jpts; + +my $qi; +my $qj; + sub triangle { my @pixs; foreach my $pval (@_) { my $pix = $point_indices{$pval} //= ((push @points, $pval), $#points); - return if grep { $pix eq $_ } @pixs; + if (grep { $pix eq $_ } @pixs) { + print "// elide @{ $ipts[$qi] } @{ $jpts[$qj] }\n"; + return; + } push @pixs, $pix; } - push @triangles, \@pixs; + push @triangles, [ $qi,$qj, \@pixs ]; } sub make_sheet () { - my @ipts = get_sections_ptvals(@i_sections); - my @jpts = get_sections_ptvals(@j_sections); + @ipts = get_sections_ptvals(@i_sections); + @jpts = get_sections_ptvals(@j_sections); my @sheet; foreach my $ipt (@ipts) { my @row = (); @@ -176,9 +200,9 @@ sub make_sheet () { } push @sheet, \@row; } - foreach (my $qi=0; $qi<$#ipts; $qi++) { # i direction does not wrap + foreach ($qi=0; $qi<$#ipts; $qi++) { # i direction does not wrap my $qi2 = $qi+1; - foreach (my $qj=0; $qj<@jpts; $qj++) { # j direction does wrap + foreach ($qj=0; $qj<@jpts; $qj++) { # j direction does wrap my $qj2 = ($qj+1) % @jpts; my $p0 = $sheet[$qi][$qj]; triangle($p0, $sheet[$qi2][$qj], $sheet[$qi2][$qj2]); @@ -196,7 +220,9 @@ sub write_out () { print "module ImplHeadCup(){ polyhedron(points=[\n" or die $!; print pv($_),",\n" or die $! foreach @points; print "],faces=[\n" or die $!; - print pv($_),",\n" or die $! foreach @triangles; + foreach (@triangles) { + print pv($_->[2]),", // @{ $ipts[$_->[0]] } @{ $jpts[$_->[1]] }\n" or die $!; + } print "],convexity=10); }\n" or die $!; print <