-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathkoefs_estimation.py
116 lines (114 loc) · 8.64 KB
/
koefs_estimation.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
import numpy as np
# http://www.aerodynamics4students.com/propulsion/propel2.php?dia=0.1778&B=3&pitch=0.1016&taper=0.4&tc=0.01&angle=0.0&option=Solve
stats = np.array([
[0.32138839787884, 4.8498935547717E+49, -8.6921212564081E+48, -0.28540183230037, 0 ],
[0.48208259681825, 2.4232212054258E+47, -4.3429680076044E+46, -0.42810274845055, 0 ],
[0.64277679575767, 1.4660673756184E+49, 4.9095927261276E+47, 3.0548422230328, 0 ],
[0.80347099469709, 5.8871698527809E+49, 1.9715060008411E+48, 3.8185527787909, 0 ],
[0.96416519363651, 1.2519035812791E+50, 4.1923971699244E+48, 4.5822633345491, 0 ],
[1.1248593925759, 2.108868096508E+50, 7.0622153109509E+48, 5.3459738903073, 0 ],
[1.2855535915153, 3.149570782219E+50, 1.0547339133226E+49, 6.1096844460655, 0 ],
[1.4462477904548, 4.3695666840932E+50, 1.4632883293992E+49, 6.8733950018237, 0 ],
[1.6069419893942, 5.7665294903116E+50, 1.9311057398499E+49, 7.6371055575819, 0 ],
[1.7676361883336, 7.3391104106035E+50, 2.4577344593693E+49, 8.4008161133401, 0 ],
[1.928330387273, 9.0864743480093E+50, 3.0428948291898E+49, 9.1645266690983, 0 ],
[2.0890245862124, 1.1008079370664E+51, 3.6864053661961E+49, 9.9282372248565, 0 ],
[2.2497187851519, 1.3103560659682E+51, 4.3881438991852E+49, 10.691947780615, 0 ],
[2.4104129840913, 1.5372665101776E+51, 5.1480256651258E+49, 11.455658336373, 0 ],
[2.5711071830307, 1.7815212524474E+51, 5.9659903275371E+49, 12.219368892131, 0 ],
[2.7318013819701, 2.0431071793688E+51, 6.8419939720005E+49, 12.983079447889, 0 ],
[2.8924955809095, 2.3220145584407E+51, 7.7760040061418E+49, 13.746790003647, 0 ],
[3.0531897798489, 2.6182360401538E+51, 8.7679958177923E+49, 14.510500559406, 0 ],
[3.2138839787884, 2.9317659878798E+51, 9.8179505309101E+49, 15.274211115164, 0 ],
[3.3745781777278, 3.2626000180093E+51, 1.0925853465585E+50, 16.037921670922, 0 ],
[3.5352723766672, 3.6107346781724E+51, 1.209169306046E+50, 16.80163222668, 0 ],
[3.6959665756066, 3.9761672180977E+51, 1.3315460105375E+50, 17.565342782438, 0 ],
[3.856660774546, 4.3588954238165E+51, 1.4597147186154E+50, 18.329053338197, 0 ],
[4.0173549734855, 4.758917495947E+51, 1.5936748276993E+50, 19.092763893955, 0 ],
[4.1780491724249, 5.1762319591374E+51, 1.7334258437209E+50, 19.856474449713, 0 ],
[4.3387433713643, 5.6108375938659E+51, 1.8789673582844E+50, 20.620185005471, 0 ],
[4.4994375703037, 6.0627333844876E+51, 2.030299031269E+50, 21.383895561229, 0 ],
[4.6601317692431, 6.5319184792477E+51, 2.1874205774373E+50, 22.147606116988, 0 ],
[4.8208259681825, 7.0183921592002E+51, 2.3503317560276E+50, 22.911316672746, 0 ],
[4.981520167122, 7.5221538138308E+51, 2.5190323625896E+50, 23.675027228504, 0 ],
[5.1422143660614, 8.0432029217764E+51, 2.6935222225284E+50, 24.438737784262, 0 ],
[5.3029085650008, 8.5815390354539E+51, 2.8738011859565E+50, 25.20244834002, 0 ],
[5.4636027639402, 9.1371617687144E+51, 3.0598691235598E+50, 25.966158895778, 0 ],
[5.6242969628796, 9.7100707868618E+51, 3.2517259232545E+50, 26.729869451537, 0 ],
[5.7849911618191, 1.0300265798523E+52, 3.449371487465E+50, 27.493580007295, 0 ],
[5.9456853607585, 1.0907746548995E+52, 3.6528057308961E+50, 28.257290563053, 0 ],
[6.1063795596979, 1.1532512814758E+52, 3.8620285786952E+50, 29.021001118811, 0 ],
[6.2670737586373, 1.2174564398938E+52, 4.0770399649323E+50, 29.784711674569, 0 ],
[6.4277679575767, 1.283390112752E+52, 4.2978398313333E+50, 30.548422230328, 0 ],
[6.5884621565161, 1.3510522846186E+52, 4.5244281262201E+50, 31.312132786086, 0 ],
[6.7491563554556, 1.4204429417658E+52, 4.7568048036202E+50, 32.075843341844, 0 ],
[6.909850554395, 1.4915620719444E+52, 4.9949698225138E+50, 32.839553897602, 0 ],
[7.0705447533344, 1.5644096641941E+52, 5.2389231461964E+50, 33.60326445336, 0 ],
[7.2312389522738, 1.6389857086807E+52, 5.4886647417355E+50, 34.366975009119, 0 ],
[7.3919331512132, 1.715290196558E+52, 5.7441945795067E+50, 35.130685564877, 0 ],
[7.5526273501527, 1.7933231198485E+52, 6.0055126327946E+50, 35.894396120635, 0 ],
[7.7133215490921, 1.8730844713414E+52, 6.2726188774518E+50, 36.658106676393, 0 ],
[7.8740157480315, 1.954574244504E+52, 6.5455132916011E+50, 37.421817232151, 0 ],
[8.0347099469709, 2.0377924334051E+52, 6.8241958553802E+50, 38.185527787909, 0 ],
[8.1954041459103, 2.1227390326483E+52, 7.1086665507178E+50, 38.949238343668, 0 ],
[8.3560983448497, 2.2094140373143E+52, 7.3989253611396E+50, 39.712948899426, 0 ],
[8.5167925437892, 2.2978174429101E+52, 7.6949722715996E+50, 40.476659455184, 0 ],
[8.6774867427286, 2.3879492453243E+52, 7.9968072683294E+50, 41.240370010942, 0 ],
[8.838180941668, 2.4798094407881E+52, 8.304430338708E+50, 42.0040805667, 0 ],
[8.9988751406074, 2.573398025841E+52, 8.6178414711468E+50, 42.767791122459, 0 ],
[9.1595693395468, 2.6687149973001E+52, 8.9370406549872E+50, 43.531501678217, 0 ],
[9.3202635384863, 2.765760352234E+52, 9.2620278804118E+50, 44.295212233975, 0 ],
[9.4809577374257, 2.8645340879375E+52, 9.592803138362E+50, 45.058922789733, 0 ],
[9.6416519363651, 2.9650362019118E+52, 9.9293664204693E+50, 45.822633345491, 0 ],
[9.8023461353045, 3.0672666918445E+52, 1.027171771899E+51, 46.58634390125, 0 ],
[9.9630403342439, 3.1712255555932E+52, 1.061985702675E+51, 47.350054457008, 0 ],
[10.123734533183, 3.2769127911708E+52, 1.0973784337094E+51, 48.113765012766, 0 ],
[10.284428732123, 3.3843283967309E+52, 1.133349964384E+51, 48.877475568524, 0 ],
[10.445122931062, 3.4934723705564E+52, 1.1699002941237E+51, 49.641186124282, 0 ],
[10.605817130002, 3.6043447110486E+52, 1.2070294223931E+51, 50.40489668004, 0 ],
[10.766511328941, 3.7169454167166E+52, 1.244737348693E+51, 51.168607235799, 0 ],
[10.92720552788, 3.8312744861695E+52, 1.2830240725575E+51, 51.932317791557, 0 ],
[11.08789972682, 3.9473319181074E+52, 1.3218895935513E+51, 52.696028347315, 0 ],
[11.248593925759, 4.0651177113147E+52, 1.3613339112675E+51, 53.459738903073, 0 ],
[11.409288124699, 4.1846318646527E+52, 1.4013570253246E+51, 54.223449458831, 0 ],
[11.569982323638, 4.305874377055E+52, 1.4419589353655E+51, 54.98716001459, 0 ],
[11.730676522578, 4.4288452475207E+52, 1.4831396410551E+51, 55.750870570348, 0 ],
[11.891370721517, 4.5535444751103E+52, 1.5248991420786E+51, 56.514581126106, 0 ],
[12.052064920456, 4.6799720589403E+52, 1.5672374381403E+51, 57.278291681864, 0 ],
[12.212759119396, 4.8081279981801E+52, 1.610154528962E+51, 58.042002237622, 0 ],
[12.373453318335, 4.9380122920473E+52, 1.6536504142818E+51, 58.805712793381, 0 ],
[12.534147517275, 5.0696249398052E+52, 1.6977250938527E+51, 59.569423349139, 0 ],
[12.694841716214, 5.2029659407585E+52, 1.7423785674422E+51, 60.333133904897, 0 ],
[12.855535915153, 5.3380352942513E+52, 1.7876108348304E+51, 61.096844460655, 0 ],
[13.016230114093, 5.4748329996638E+52, 1.8334218958098E+51, 61.860555016413, 0 ],
[13.176924313032, 5.6133590564102E+52, 1.8798117501843E+51, 62.624265572172, 0 ],
[13.337618511972, 5.7536134639368E+52, 1.9267803977684E+51, 63.38797612793, 0 ],
[13.498312710911, 5.8955962217189E+52, 1.9743278383865E+51, 64.151686683688, 0 ],
[13.659006909851, 6.0393073292598E+52, 2.0224540718721E+51, 64.915397239446, 0 ],
[13.81970110879, 6.1847467860884E+52, 2.0711590980675E+51, 65.679107795204, 0 ],
[13.980395307729, 6.3319145917577E+52, 2.1204429168231E+51, 66.442818350963, 0 ],
[14.141089506669, 6.4808107458443E+52, 2.1703055279969E+51, 67.206528906721, 0 ],
[14.301783705608, 6.631435247945E+52, 2.2207469314542E+51, 67.970239462479, 0 ],
[14.462477904548, 6.7837880976771E+52, 2.2717671270665E+51, 68.733950018237, 0 ],
[14.623172103487, 6.9378692946764E+52, 2.3233661147121E+51, 69.497660573995, 0 ],
[14.783866302426, 7.0936788385973E+52, 2.3755438942751E+51, 70.261371129753, 0 ],
[14.944560501366, 7.2512167291091E+52, 2.4283004656448E+51, 71.025081685512, 0 ],
[15.105254700305, 7.4104829658982E+52, 2.4816358287162E+51, 71.78879224127, 0 ],
[15.265948899245, 7.5714775486654E+52, 2.5355499833891E+51, 72.552502797028, 0 ],
[15.426643098184, 7.7342004771249E+52, 2.5900429295678E+51, 73.316213352786, 0 ],
[15.587337297124, 7.8986517510043E+52, 2.645114667161E+51, 74.079923908544, 0 ],
[15.748031496063, 8.0648313700434E+52, 2.7007651960817E+51, 74.843634464303, 0 ],
[15.908725695002, 8.2327393339939E+52, 2.7569945162466E+51, 75.607345020061, 0],
[ 16.069419893942, 8.402375642619E+52, 2.8138026275765E+51, 76.371055575819, 0],
])
J = stats[:, 0]
T = stats[:, 1]
Q = stats[:, 2]
Eff = stats[:, 3]
n = 2100/60
rho = 1.225
D = Eff*Q/(J*T)
Ct = T/(rho*(n**2)*(D**4))
Cq = Q/(rho*(n**2)*(D**5))
# print(Eff*Q/(J*T))
print(Ct)