Mercurial > hg > weather-synthesis-prototyping

annotate MechAeolian/newStrou.pd @ 223:3829344821ce

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Test Aeolian harp added. 6 strings with various diameters, tensions, mass and lengths. More physics to be added but not too bad a starting place. Sigmoids added for the lock in region. No hysteresis yet. No dampening factor for higher harmonics yet. Could do with an indicator on harpD1 that shows the natural freq, strouhal freq, harmonic number and if on / off.
author Rod Selfridge <r.selfridge@qmul.ac.uk>
date Sat, 28 May 2016 08:11:09 +0100
parents 45c696ea3247
children
rev   line source
r@223 1 #N canvas 179 322 718 433 10;
r@197 2 #X obj 203 48 inlet~;
r@197 3 #X obj 273 202 outlet~;
r@197 4 #X text 101 20 Reynolds number coming in;
r@197 5 #X text 211 238 Strouhal Number going out;
r@197 6 #N canvas 679 277 1026 531 eqnSelect 0;
r@197 7 #X msg 41 221 0.2684 -1.0356;
r@197 8 #X msg 148 221 0.2437 -0.8607;
r@197 9 #X msg 256 223 0.4291 -3.6735;
r@197 10 #X msg 472 223 0.2257 -0.4402;
r@197 11 #X msg 577 222 0.204 0.3364;
r@197 12 #X msg 686 224 0.1776 2.2023;
r@197 13 #X obj 365 299 unpack;
r@197 14 #X obj 75 164 < 180;
r@197 15 #X obj 181 144 > 180;
r@197 16 #X obj 227 146 < 230;
r@197 17 #X obj 289 149 > 230;
r@197 18 #X obj 335 148 < 240;
r@197 19 #X obj 398 148 > 240;
r@197 20 #X obj 444 149 < 360;
r@197 21 #X obj 505 147 > 360;
r@197 22 #X obj 551 147 < 1300;
r@197 23 #X obj 611 150 > 1300;
r@197 24 #X obj 657 150 < 5000;
r@197 25 #X obj 719 149 > 5000;
r@197 26 #X obj 765 148 < 200000;
r@197 27 #X obj 827 149 > 200000;
r@197 28 #X obj 41 196 spigot;
r@197 29 #X floatatom 286 327 8 0 0 0 - - -;
r@197 30 #X floatatom 455 330 8 0 0 0 - - -;
r@197 31 #X obj 148 198 spigot;
r@197 32 #X obj 194 171 spigot;
r@197 33 #X obj 256 200 spigot;
r@197 34 #X obj 302 173 spigot;
r@197 35 #X obj 365 200 spigot;
r@197 36 #X obj 411 173 spigot;
r@197 37 #X floatatom 546 83 8 0 0 0 - - -;
r@197 38 #X obj 472 198 spigot;
r@197 39 #X obj 518 171 spigot;
r@197 40 #X obj 578 200 spigot;
r@197 41 #X obj 624 173 spigot;
r@197 42 #X obj 686 199 spigot;
r@197 43 #X obj 732 172 spigot;
r@197 44 #X obj 426 67 snapshot~;
r@197 45 #X obj 426 37 inlet~;
r@197 46 #X obj 365 373 outlet~;
r@197 47 #X obj 398 352 outlet~;
r@197 48 #X obj 365 327 sig~;
r@197 49 #X obj 398 327 sig~;
r@197 50 #X text 539 419 This selects the values for the St* and m coefficients
r@197 51 for Fey's formula;
r@197 52 #X obj 794 201 spigot;
r@197 53 #X obj 840 174 spigot;
r@197 54 #X obj 885 149 < 1e+06;
r@197 55 #X text 784 251 Critical - Derived by author;
r@197 56 #X text 254 359 St* (I) coefficient;
r@197 57 #X text 462 361 m (J) coefficient;
r@197 58 #X msg 794 226 0.576 -175.957;
r@208 59 #X text 372 240 comment;
r@197 60 #X msg 365 223 0.2492 -0.8861;
r@197 61 #X text 370 240 Linear Int;
r@208 62 #X obj 359 43 r myMet;
r@197 63 #X connect 0 0 6 0;
r@197 64 #X connect 1 0 6 0;
r@197 65 #X connect 2 0 6 0;
r@197 66 #X connect 3 0 6 0;
r@197 67 #X connect 4 0 6 0;
r@197 68 #X connect 5 0 6 0;
r@197 69 #X connect 6 0 22 0;
r@208 70 #X connect 6 0 41 0;
r@197 71 #X connect 6 1 23 0;
r@208 72 #X connect 6 1 42 0;
r@197 73 #X connect 7 0 21 1;
r@197 74 #X connect 7 0 21 0;
r@197 75 #X connect 8 0 24 1;
r@197 76 #X connect 8 0 25 0;
r@197 77 #X connect 9 0 25 1;
r@197 78 #X connect 10 0 27 0;
r@197 79 #X connect 10 0 26 1;
r@197 80 #X connect 11 0 27 1;
r@197 81 #X connect 12 0 28 1;
r@197 82 #X connect 12 0 29 0;
r@197 83 #X connect 13 0 29 1;
r@197 84 #X connect 14 0 31 1;
r@197 85 #X connect 14 0 32 0;
r@197 86 #X connect 15 0 32 1;
r@197 87 #X connect 16 0 33 1;
r@197 88 #X connect 16 0 34 0;
r@197 89 #X connect 17 0 34 1;
r@197 90 #X connect 18 0 35 1;
r@197 91 #X connect 18 0 36 0;
r@197 92 #X connect 19 0 36 1;
r@208 93 #X connect 20 0 45 0;
r@208 94 #X connect 20 0 44 1;
r@197 95 #X connect 21 0 0 0;
r@197 96 #X connect 24 0 1 0;
r@197 97 #X connect 25 0 24 0;
r@197 98 #X connect 26 0 2 0;
r@197 99 #X connect 27 0 26 0;
r@208 100 #X connect 28 0 52 0;
r@197 101 #X connect 29 0 28 0;
r@197 102 #X connect 31 0 3 0;
r@197 103 #X connect 32 0 31 0;
r@197 104 #X connect 33 0 4 0;
r@197 105 #X connect 34 0 33 0;
r@197 106 #X connect 35 0 5 0;
r@197 107 #X connect 36 0 35 0;
r@197 108 #X connect 37 0 7 0;
r@197 109 #X connect 37 0 8 0;
r@197 110 #X connect 37 0 9 0;
r@197 111 #X connect 37 0 10 0;
r@197 112 #X connect 37 0 11 0;
r@197 113 #X connect 37 0 12 0;
r@197 114 #X connect 37 0 13 0;
r@197 115 #X connect 37 0 30 0;
r@197 116 #X connect 37 0 14 0;
r@197 117 #X connect 37 0 15 0;
r@197 118 #X connect 37 0 16 0;
r@197 119 #X connect 37 0 17 0;
r@197 120 #X connect 37 0 18 0;
r@197 121 #X connect 37 0 19 0;
r@197 122 #X connect 37 0 20 0;
r@208 123 #X connect 37 0 46 0;
r@197 124 #X connect 38 0 37 0;
r@208 125 #X connect 41 0 39 0;
r@208 126 #X connect 42 0 40 0;
r@208 127 #X connect 44 0 50 0;
r@208 128 #X connect 45 0 44 0;
r@208 129 #X connect 46 0 45 1;
r@208 130 #X connect 50 0 6 0;
r@208 131 #X connect 52 0 6 0;
r@208 132 #X connect 54 0 37 0;
r@197 133 #X restore 273 83 pd eqnSelect;
r@197 134 #X obj 203 85 sqrt~;
r@197 135 #X obj 342 119 /~;
r@197 136 #X text 135 97 sqrt(Re);
r@197 137 #X text 250 114 St*;
r@197 138 #X text 359 101 m;
r@197 139 #X text 376 152 St = St* + m / sqrt(Re);
r@197 140 #X obj 273 166 +~;
r@197 141 #X connect 0 0 4 0;
r@197 142 #X connect 0 0 5 0;
r@208 143 #X connect 4 0 11 0;
r@197 144 #X connect 4 1 6 0;
r@197 145 #X connect 5 0 6 1;
r@208 146 #X connect 6 0 11 1;
r@208 147 #X connect 11 0 1 0;