Welcome to MakerHome




We've completed our yearlong print-a-day project!
All new material is now at Hacktastic: www.mathgrrl.com


Saturday, January 25, 2014

Day 152 - Petal and tight trefoil knots

Two more trefoil conformations - as a petal knot (all crossings line up along one center stem) and a tight knot (with the smallest amount of rope length possible given its diameter):


STL file for petal: http://www.geekhaus.com/makerhome/day152_trefoil_petal_40_25_3.stl
STL file for tight: http://www.geekhaus.com/makerhome/day152_trefoil_minimumrope_25.stl
Thingiverse link: http://www.thingiverse.com/thing:234107

Settings: MakerWare .3mm/low in 20-25 minutes for the petal knot and even less for the tight knot, with the minimum-support custom slicing profile for knots described on Day 110.

Technical notes, math flavor: The beautiful petal knot conformation is from the paper Bounds on Übercrossing and Petal Numbers for Knots by Colin Adams at Williams College and his amazing team of student researchers. In this work Adams, et al made the startling decision to study "petal" projections of knots where the crossings occur all on top of each other in one place! This is a wild and crazy idea compared to the usual notion that we should look at knot projections where crossings happen nicely and one at a time. But this crazy idea bore some amazing fruit: they proved that every knot admits a petal projection, and went on to discover deep connections between the übercrossing number of knots and other traditional invariants such as crossing and unknotting numbers. I know I'm not supposed to have favorites, but of the eight trefoil knots I am printing in this collection, this is by far my favorite. It's beautiful and surprising each time you turn the knot over and the perspective shifts from one-crossing-at-a-time to all-crossings-lined-up-at-once. Wonderful! See the OpenSCAD code below if you want a (piecewise) parametrization of this knot.

Minimum-length rope knot conformations of the tiny trefoil aren't the most impressive, but for higher-crossing knots they are fascinating, and what started my journey with 3D printing last year. Jason Cantarella at the University of Georgia and Eric Rawdon at the University of St. Thomas, together with some students, have been investigating "tight" knots since their paper Knot Tightening by Constrained Gradient Descent in which they used their ridgerunner algorithm to find minimal conformations of a large number of knots. The data they made available with this paper was what helped me model my first 3D-printable knot models (for example, see Day 9Day 11Day 66, and Day 67).

Technical notes, OpenSCAD flavor:  By the time of this posting, kitwallace and the GitHub OpenSCAD community had improved this code significantly; see kitwallace's excellent post at his blog The Wallace Line. However the new code requires the new 2014.01.04 snapshot of OpenSCAD, which does not want to work on my machine at the moment. Therefore I will post the old, slow code. It works but it ain't fast. In fact, it ain't pretty either, since getting Adams' petal knot parametrization to work in OpenSCAD required creating a very clunky piecewise function, as you will see in the code below. And it just gets worse for the tight knot; I have a very inefficient piece of code that takes 50 coordinates from Cantarella's data and connects the dots with hulled pairs of spheres. I am sure a nice loop would have made quicker work of the problem, but at the time my clunky mess of code was easier to write.  I've put that code under the fold so you have to click through to see my code of shame.

// mathgrrl parametric knots
// tubify module based on tube module from kitwallace
// choose the knot you want and comment out the rest

// trefoil as petal knot
// http://arxiv.org/pdf/1311.0526v1.pdf
// based on mma algorithm from Adams, need to do z-coord in cases
// scaled to 40mm before tubifying
tempstep = .005; //.005/1 is 2/360
function g(t,a,b,u,v) = 
[ 3.5*5*cos(180*t)*sin(5*180*t), 
 3.5*5*sin(180*t)*sin(5*180*t),
 8*u*(cos(180*(((5*t-a+1)%5)-1)/2))*(cos(180*(((5*t-a+1)%5)-1)/2)) 
+ 8*v*(cos(180*(((5*t-b+1)%5)-1)/2))*(cos(180*(((5*t-b+1)%5)-1)/2)) 
];
translate([0,0,-24])
union(){
color("red")
for (t=[1: tempstep: 1.2-tempstep]) { //replacement for 0 to .2
hull() {
translate(g(t,5,1,4,1)) sphere(2.5); 
translate(g(t+tempstep,5,1,4,1)) sphere(2.5);
      }
  }
color("orange")
for (t=[.2: tempstep: .4]) { //also why here?
hull() {
translate(g(t,1,2,1,3)) sphere(2.5);
translate(g(t+tempstep,1,2,1,3)) sphere(2.5);
      }
  }
color("yellow")
for (t=[.4: tempstep: .6]) {
hull() {
translate(g(t,2,3,3,5)) sphere(2.5);
translate(g(t+tempstep,2,3,3,5)) sphere(2.5);
      }
  }
color("green")
for (t=[.6: tempstep: .8]) {
hull() {
translate(g(t,3,4,5,2)) sphere(2.5);
translate(g(t+tempstep,3,4,5,2)) sphere(2.5);
      }
  }
color("blue")
for (t=[.8: tempstep: 1]) {
hull() {
translate(g(t,4,5,2,4)) sphere(2.5);
translate(g(t+tempstep,4,5,2,4)) sphere(2.5);
      }
  }
}


// trefoil in minimal stick conformation
// http://newweb.cecm.sfu.ca/cgi-bin/KnotPlot/KnotServer/kserver?ncomp=1&ncross=3&id=1
s = 5.1;
r=2.5;
p1=[1.110299950990041*s, 1.429392328828375*s, 1.668976011298993*s];
p2=[0.7849508316994853*s, 1.390527974091808*s, 1.638400318090565*s];
p3=[0.4839768360110871*s, 1.408309308805168*s, 1.507456276424959*s];
p4=[0.2099637867533136*s, 1.526913147575392*s, 1.368746074544028*s];
p5=[-0.1044045199364073*s, 1.557072813802938*s, 1.277467938632125*s];
p6=[-0.399748259645602*s, 1.475223222916731*s, 1.157624027158315*s];
p7=[-0.6358281675566203*s, 1.307016606241804*s, 1.001722605280547*s];
p8=[-0.7851542732934931*s, 1.079193042297195*s, 0.8169408434548786*s];
p9=[-0.8311134346082066*s, 0.8212731304859858*s, 0.6176515842442668*s];
p10=[-0.7686493094139212*s, 0.5642909832121883*s, 0.4216630135406117*s];
p11=[-0.6045929612305532*s, 0.3382007063288971*s, 0.2475245715157406*s];
p12=[-0.3568306032296734*s, 0.1681613192110579*s, 0.1131468482930633*s];
p13=[-0.05215550025598645*s, 0.07145689451092788*s, 0.03453515132340339*s];
p14=[0.276462228240339*s, 0.05841653174602905*s, 0.02048291528081363*s];
p15=[0.5933655543230373*s, 0.1320624499868605*s, 0.07053994990181335*s];
p16=[0.8640791760231149*s, 0.2868919207950811*s, 0.1758864037079465*s];
p17=[1.059007552878924*s, 0.5094683811641367*s, 0.3201501562920863*s];
p18=[1.156587739687891*s, 0.7806329930439322*s, 0.4791294447142897*s];
p19=[1.146173200556639*s, 1.078483098705804*s, 0.6184647171807074*s];
p20=[1.046652555232624*s, 1.386552577187142*s, 0.6718072391049106*s];
p21=[1.004901905630199*s, 1.710284202663531*s, 0.7043995393056642*s];
p22=[0.9748500956671001*s, 2.00047580949448*s, 0.8564480524580105*s];
p23=[0.9098334192445141*s, 2.222105247942773*s, 1.090908692836742*s];
p24=[0.8019842713353851*s, 2.354321803720973*s, 1.372384217563913*s];
p25=[0.6551799860417229*s, 2.385578042488897*s, 1.665336015555774*s];
p26=[0.4805521535530748*s, 2.314511642494548*s, 1.935163624031325*s];
p27=[0.2936491326895583*s, 2.150349859052606*s, 2.150737387641385*s];
p28=[0.1127496269984707*s, 1.911868236163361*s, 2.287700493082868*s];
p29=[-0.0421605287380469*s, 1.624736696511827*s, 2.331471080472086*s];
p30=[-0.1528609563612467*s, 1.319365053676113*s, 2.27807485496556*s];
p31=[-0.2083358750431886*s, 1.029312266757455*s, 2.132647792627378*s];
p32=[-0.2054199993164743*s, 0.7873646838852222*s, 1.909466956355828*s];
p33=[-0.1485607220853743*s, 0.6216668652420217*s, 1.630782548477369*s];
p34=[-0.05040230980381472*s, 0.5531146846464222*s, 1.324177893798596*s];
p35=[0.06640401809921728*s, 0.5938663837770654*s, 1.019230285160021*s];
p36=[0.1613979533593562*s, 0.7465064519797134*s, 0.7440500547667859*s];
p37=[0.1612865861469598*s, 0.9416891449937117*s, 0.4789246463874046*s];
p38=[0.2339326806623845*s, 1.063305993816675*s, 0.1822750032128956*s];
p39=[0.4144963648776064*s, 1.167976155984675*s, -0.07214765172111942*s];
p40=[0.6734900953755177*s, 1.268817686183782*s, -0.2484598062187317*s];
p41=[0.9780770272754852*s, 1.366881693952402*s, -0.3256381806808795*s];
p42=[1.292722409633288*s, 1.458539532011356*s, -0.2948438385700362*s];
p43=[1.581835325337929*s, 1.53846510352008*s, -0.1592754396998342*s];
p44=[1.813120663011601*s, 1.601258859452873*s, 0.06637823224367628*s];
p45=[1.961278333610819*s, 1.641247203005984*s, 0.3575950904845052*s];
p46=[2.010777383524477*s, 1.652331218964934*s, 0.6828419506910011*s];
p47=[1.956284238411966*s, 1.633000534091536*s, 1.006900618530711*s];
p48=[1.802938617261106*s, 1.587763161772232*s, 1.294643490310422*s];
p49=[1.566309155931087*s, 1.525603842329356*s, 1.514867637084154*s];
p50=[1.270770899793466*s, 1.459136348798676*s, 1.643664279859707*s];
hull(){
translate(p1) sphere(r);
translate(p2) sphere(r);
}
hull(){
translate(p2) sphere(r);
translate(p3) sphere(r);
}
hull(){
translate(p3) sphere(r);
translate(p4) sphere(r);
}
hull(){
translate(p4) sphere(r);
translate(p5) sphere(r);
}
hull(){
translate(p5) sphere(r);
translate(p6) sphere(r);
}
hull(){
translate(p6) sphere(r);
translate(p7) sphere(r);
}
hull(){
translate(p7) sphere(r);
translate(p8) sphere(r);
}
hull(){
translate(p8) sphere(r);
translate(p9) sphere(r);
}
hull(){
translate(p9) sphere(r);
translate(p10) sphere(r);
}
hull(){
translate(p10) sphere(r);
translate(p11) sphere(r);
}
hull(){
translate(p11) sphere(r);
translate(p12) sphere(r);
}
hull(){
translate(p12) sphere(r);
translate(p13) sphere(r);
}
hull(){
translate(p13) sphere(r);
translate(p14) sphere(r);
}
hull(){
translate(p14) sphere(r);
translate(p15) sphere(r);
}
hull(){
translate(p15) sphere(r);
translate(p16) sphere(r);
}
hull(){
translate(p16) sphere(r);
translate(p17) sphere(r);
}
hull(){
translate(p17) sphere(r);
translate(p18) sphere(r);
}
hull(){
translate(p18) sphere(r);
translate(p19) sphere(r);
}
hull(){
translate(p19) sphere(r);
translate(p20) sphere(r);
}
hull(){
translate(p20) sphere(r);
translate(p21) sphere(r);
}
hull(){
translate(p21) sphere(r);
translate(p22) sphere(r);
}
hull(){
translate(p22) sphere(r);
translate(p23) sphere(r);
}
hull(){
translate(p23) sphere(r);
translate(p24) sphere(r);
}
hull(){
translate(p24) sphere(r);
translate(p25) sphere(r);
}
hull(){
translate(p25) sphere(r);
translate(p26) sphere(r);
}
hull(){
translate(p26) sphere(r);
translate(p27) sphere(r);
}
hull(){
translate(p27) sphere(r);
translate(p28) sphere(r);
}
hull(){
translate(p28) sphere(r);
translate(p29) sphere(r);
}
hull(){
translate(p29) sphere(r);
translate(p30) sphere(r);
}
hull(){
translate(p30) sphere(r);
translate(p31) sphere(r);
}
hull(){
translate(p31) sphere(r);
translate(p32) sphere(r);
}
hull(){
translate(p32) sphere(r);
translate(p33) sphere(r);
}
hull(){
translate(p33) sphere(r);
translate(p34) sphere(r);
}
hull(){
translate(p34) sphere(r);
translate(p35) sphere(r);
}
hull(){
translate(p35) sphere(r);
translate(p36) sphere(r);
}
hull(){
translate(p36) sphere(r);
translate(p37) sphere(r);
}
hull(){
translate(p37) sphere(r);
translate(p38) sphere(r);
}
hull(){
translate(p38) sphere(r);
translate(p39) sphere(r);
}
hull(){
translate(p39) sphere(r);
translate(p40) sphere(r);
}
hull(){
translate(p40) sphere(r);
translate(p41) sphere(r);
}
hull(){
translate(p41) sphere(r);
translate(p42) sphere(r);
}
hull(){
translate(p42) sphere(r);
translate(p43) sphere(r);
}
hull(){
translate(p43) sphere(r);
translate(p44) sphere(r);
}
hull(){
translate(p44) sphere(r);
translate(p45) sphere(r);
}
hull(){
translate(p45) sphere(r);
translate(p46) sphere(r);
}
hull(){
translate(p46) sphere(r);
translate(p47) sphere(r);
}
hull(){
translate(p47) sphere(r);
translate(p48) sphere(r);
}
hull(){
translate(p48) sphere(r);
translate(p49) sphere(r);
}
hull(){
translate(p49) sphere(r);
translate(p50) sphere(r);
}
hull(){
translate(p50) sphere(r);
translate(p1) sphere(r);
}

1 comment:

  1. you have lots of perseverance!
    can you print the Lorenz attractor come from the chaos theory.
    i think it is very special but also hard.

    ReplyDelete