Scaling1 plot update; within-isoline

avikde commented 4 years ago

After #133

Do we need to vary Rpow while doing this?

avikde commented 4 years ago

First naively redo it after #138

Then the best thing will be to create plots that show divergence from naive scaling. scaling1_138.pdf scaling1_138.zip

avikde commented 4 years ago

Scaling2 - same lift, just change dt weight

get force, power with different trends. can produce a plot with the ideal design in each case

This is with phi fixed at 120

scaling2_300al.zip
Trying to get something for the HalfBee1 design
scaling2_220al.zip
-> Need to lower T

Lower phi to 90->80

scaling2_220al_90.zip

Opt Φ=90, Rpow=0.0, minal=220, τ2/1 lim=2 => -2, [21.722 2.496 1.129 5.265 3.79 86.888 0.083], fHz=149.9, al[mg]=226.9, u∞=126.6, FD∞=81.0, pow=17.7, J=118.3, AR=4.0, x=29.3
Opt Φ=90, Rpow=0.0, minal=220, τ2/1 lim=2 => -2, [20.782 2.515 1.084 5.208 3.177 83.409 0.08], fHz=156.2, al[mg]=227.2, u∞=125.7, FD∞=79.2, pow=18.1, J=132.7, AR=4.0, x=28.7
Opt Φ=90, Rpow=0.0, minal=220, τ2/1 lim=2 => 0, [20.013 2.491 1.04 5.155 3.941 80.05 0.077], fHz=162.7, al[mg]=224.5, u∞=123.7, FD∞=75.0, pow=18.3, J=147.5, AR=4.0, x=28.1
Opt Φ=90, Rpow=0.0, minal=220, τ2/1 lim=2 => 0, [19.017 2.494 1.001 5.078 3.837 77.044 0.074], fHz=169.1, al[mg]=226.2, u∞=126.9, FD∞=73.7, pow=18.7, J=161.9, AR=4.1, x=27.8
Opt Φ=90, Rpow=0.0, minal=220, τ2/1 lim=2 => -2, [17.787 2.507 0.963 4.965 3.399 74.124 0.071], fHz=175.7, al[mg]=231.1, 
u∞=131.3, FD∞=73.9, pow=19.5, J=175.0, AR=4.2, x=27.6
Opt Φ=90, Rpow=0.0, minal=220, τ2/1 lim=2 => 0, [16.975 2.508 0.934 4.879 3.374 71.936 0.069], fHz=181.1, al[mg]=233.8, u∞=134.8, FD∞=73.5, pow=20.0, J=187.5, AR=4.2, x=27.4
Opt Φ=90, Rpow=0.0, minal=220, τ2/1 lim=2 => 0, [16.122 2.512 0.906 4.782 3.251 69.763 0.067], fHz=186.7, al[mg]=237.5, u∞=138.3, FD∞=73.7, pow=20.6, J=199.0, AR=4.3, x=27.3
Opt Φ=90, Rpow=0.0, minal=220, τ2/1 lim=2 => 0, [15.456 2.518 0.883 4.7 3.046 67.979 0.065], fHz=191.6, al[mg]=240.3, u∞=140.8, FD∞=73.8, pow=21.1, J=210.0, AR=4.4, x=27.2
Opt Φ=90, Rpow=0.0, minal=220, τ2/1 lim=2 => -2, [14.878 2.52 0.863 4.64 2.986 66.474 0.064], fHz=196.0, al[mg]=243.2, u∞=143.7, FD∞=74.1, pow=21.6, J=220.6, AR=4.5, x=27.1
Opt Φ=90, Rpow=0.0, minal=220, τ2/1 lim=2 => 0, [14.396 2.523 0.846 4.58 2.89 65.116 0.062], fHz=200.0, al[mg]=245.0, u∞=145.7, FD∞=74.0, pow=22.0, J=230.7, AR=4.5, x=27.0

scaling2_180al_90.zip

Opt Φ=90, Qdt=1000.0, minal=180, τ2/1 lim=2 => -2, [19.352 2.462 1.019 4.281 4.924 78.48 0.083], fHz=150.1, al[mg]=187.7, u∞=103.1, FD∞=64.0, pow=14.3, J=86.4, AR=4.1, x=28.0
Opt Φ=90, Qdt=25750.0, minal=180, τ2/1 lim=2 => -2, [14.6 2.462 0.839 3.918 4.922 64.596 0.069], fHz=182.4, al[mg]=199.1, u∞=98.0, FD∞=60.6, pow=16.5, J=152.1, AR=4.4, x=26.6
Opt Φ=90, Qdt=50500.0, minal=180, τ2/1 lim=2 => 0, [12.168 2.462 0.75 3.659 4.924 57.758 0.061], fHz=204.0, al[mg]=209.2, u∞=109.6, FD∞=60.5, pow=18.1, J=200.2, AR=4.7, x=26.0
Opt Φ=90, Qdt=75250.0, minal=180, τ2/1 lim=2 => 0, [10.872 2.48 0.698 3.518 4.297 53.775 0.057], fHz=219.1, al[mg]=214.0, u∞=115.4, FD∞=59.5, pow=19.2, J=241.2, AR=4.9, x=25.6
Opt Φ=90, Qdt=100000.0, minal=180, τ2/1 lim=2 => 0, [9.917 2.527 0.659 3.402 2.759 50.749 0.054], fHz=232.2, al[mg]=217.9, u∞=117.9, FD∞=59.2, pow=20.1, J=277.4, AR=5.1, x=25.

Choosing ret2 = @time opt1(m, ret1["traj"], ret1["param"], 1, 180; Φ=90, Qdt=5e4) Opt Φ=90, Qdt=50000.0, minal=180, τ2/1 lim=2.0 => 0, [12.238 2.462 0.752 3.666 4.924 57.955 0.061], fHz=203.3, al[mg]=208.9, u∞=109.2, FD∞=60.5, pow=18.1, J=199.4, AR=4.7, x=26.0