*ARCHIVED* development moved to aircraft-studio.
1
2
/home/blendux/github/UCLA_MAE_154B/example_airfoil.py,842
3
start_time = time.time()start_time15,374
4
NACA_NUM = 2412NACA_NUM18,426
5
CHORD_LENGTH = 68CHORD_LENGTH19,442
6
SEMI_SPAN = 150SEMI_SPAN20,460
7
SPAR_THICKNESS = 0.4SPAR_THICKNESS23,491
8
SKIN_THICKNESS = 0.1SKIN_THICKNESS24,512
9
AIRFOIL_MASS = 10AIRFOIL_MASS27,559
10
SPAR_MASS = 10SPAR_MASS28,577
11
STRINGER_MASS = 5STRINGER_MASS29,592
12
SPAR_CAP_AREA = 0.3SPAR_CAP_AREA32,625
13
STRINGER_AREA = 0.1STRINGER_AREA33,645
14
TOP_STRINGERS = 6TOP_STRINGERS36,688
15
BOTTOM_STRINGERS = 4BOTTOM_STRINGERS37,706
16
NOSE_TOP_STRINGERS = 3NOSE_TOP_STRINGERS38,727
17
NOSE_BOTTOM_STRINGERS = 5NOSE_BOTTOM_STRINGERS39,750
18
SAVE_PATH = '/home/blendux/github/UCLA_MAE_154B/save/'SAVE_PATH41,777
19
af = creator.Airfoil.from_dimensions(CHORD_LENGTH, SEMI_SPAN)af45,860
20
eval = evaluator.Evaluator(af)eval81,2026
21
pop = generator.Population(10)pop89,2245
22
23
/home/blendux/github/UCLA_MAE_154B/wing_scripts/my_progress.m,7747
24
Vx = 1; Vz = 1; My = 1; %test loads will be applied individuallyVx5,40
25
numStringers = numTopStringers;numStringers15,163
26
stringerGap = upperStringerGap;stringerGap16,195
27
webThickness = t_upper;webThickness17,227
28
tempStringers = topStringers;tempStringers18,251
29
dx = web(i).xStart-centroid.posX; dz = web(i).zStart-centroid.posZ;dx34,786
30
tempInt = get_int(web(i).xStart/chord,web(i).xEnd/chord,1)*chord^2; %integral of airfoil futempInt40,1142
31
triangle1 = abs( (web(i).xStart - sparCaps(1).posX)*web(i).zStart/2);triangle141,1245
32
triangle2 = abs((web(i).xEnd - sparCaps(1).posX)*web(i).zEnd/2);triangle242,1319
33
webTop = web;webTop56,2080
34
web = [];web57,2094
35
dx = web(i).xStart-centroid.posX; dz = web(i).zStart-centroid.posZ;dx67,2319
36
webRearSpar = web;webRearSpar87,3367
37
web = [];web88,3386
38
numStringers = numBottomStringers;numStringers92,3410
39
stringerGap = lowerStringerGap;stringerGap93,3445
40
webThickness = t_lower;webThickness94,3477
41
tempStringers = bottomStringers;tempStringers95,3501
42
dx = web(i).xStart-centroid.posX; dz = web(i).zStart-centroid.posZ;dx103,3817
43
tempInt = get_int(web(i).xEnd/chord,web(i).xStart/chord,0)*chord^2; %integral of airfoil futempInt118,4524
44
triangle2 = abs((web(i).xStart - sparCaps(1).posX)*web(i).zStart/2);triangle2119,4627
45
triangle1 = abs((web(i).xEnd - sparCaps(1).posX)*web(i).zEnd/2);triangle1120,4700
46
webBottom = web;webBottom136,5535
47
web = [];web137,5552
48
dx = web(i).xStart-centroid.posX; dz = web(i).zStart-centroid.posZ;dx147,5779
49
webFrontSpar = web;webFrontSpar165,6715
50
web = [];web166,6735
51
numStringers = numNoseBottomStringers;numStringers174,6781
52
stringerGap = lowerNoseStringerGap;stringerGap175,6820
53
webThickness = t_lower_front;webThickness176,6856
54
tempStringers = noseBottomStringers;tempStringers177,6886
55
dx = web(i).xStart-centroid.posX; dz = web(i).zStart-centroid.posZ;dx185,7206
56
tempInt = get_int(web(i).xEnd/chord,web(i).xStart/chord,0)*chord^2; %integral of airfoil futempInt200,7761
57
triangle1 = abs((web(i).xStart - sparCaps(2).posX)*web(i).zStart/2);triangle1201,7864
58
triangle2 = abs((web(i).xEnd - sparCaps(2).posX)*web(i).zEnd/2);triangle2202,7937
59
webLowerNose = web;webLowerNose218,8776
60
web = [];web219,8796
61
numStringers = numNoseTopStringers;numStringers222,8824
62
stringerGap = upperNoseStringerGap;stringerGap223,8860
63
webThickness = t_upper_front;webThickness224,8896
64
tempStringers = noseTopStringers;tempStringers225,8926
65
dx = web(i).xStart-centroid.posX; dz = web(i).zStart-centroid.posZ;dx233,9205
66
tempInt = get_int(web(i).xStart/chord,web(i).xEnd/chord,1)*chord^2; %integral of airfoil futempInt247,9864
67
triangle2 = abs((web(i).xStart - sparCaps(2).posX)*web(i).zStart/2);triangle2248,9967
68
triangle1 = abs((web(i).xEnd - sparCaps(2).posX)*web(i).zEnd/2);triangle1249,10040
69
webUpperNose = web;webUpperNose264,10802
70
web = [];web265,10822
71
dx = web(i).xStart-centroid.posX; dz = web(i).zStart-centroid.posZ;dx276,11050
72
webFrontSparCell2 = web;webFrontSparCell2294,11994
73
web = [];web295,12019
74
Fx = sum([webTop.qp_dx_X])+webRearSpar.qp_dx_X+ sum([webBottom.qp_dx_X])+webFrontSpar.qp_dx_X; Fx300,12064
75
Fx = Fx + sum([webLowerNose.qp_dx_X])+ sum([webUpperNose.qp_dx_X]); %cell 2Fx301,12168
76
Fz = sum([webTop.qp_dz_X])+webRearSpar.qp_dz_X+ sum([webBottom.qp_dz_X])+webFrontSpar.qp_dz_X; Fz303,12248
77
Fz = Fz + sum([webLowerNose.qp_dz_X])+ sum([webUpperNose.qp_dz_X]); %cell 2Fz304,12352
78
Fx = sum([webTop.qp_dx_Z])+webRearSpar.qp_dx_Z+ sum([webBottom.qp_dx_Z])+webFrontSpar.qp_dx_Z; Fx310,12467
79
Fx = Fx + sum([webLowerNose.qp_dx_Z])+ sum([webUpperNose.qp_dx_Z]); %cell 2Fx311,12571
80
Fz = sum([webTop.qp_dz_Z])+webRearSpar.qp_dz_Z+ sum([webBottom.qp_dz_Z])+webFrontSpar.qp_dz_Z; Fz313,12651
81
Fz = Fz + sum([webLowerNose.qp_dz_Z])+ sum([webUpperNose.qp_dz_Z]); %cell 2Fz314,12755
82
A11 = sum([webTop.dS_over_t])+webRearSpar.dS_over_t+ sum([webBottom.dS_over_t])+webFrontSpar.dS_A11323,12925
83
A22 = sum([webLowerNose.dS_over_t])+ sum([webUpperNose.dS_over_t])+webFrontSparCell2.dS_over_t;A22324,13029
84
A12 = -webFrontSpar.dS_over_t;A12325,13125
85
A21 = -webFrontSparCell2.dS_over_t;A21326,13156
86
B1_X = sum([webTop.q_dS_over_t_X])+webRearSpar.q_dS_over_t_X+ sum([webBottom.q_dS_over_t_X])+webB1_X328,13193
87
B2_X = sum([webLowerNose.q_dS_over_t_X])+ sum([webUpperNose.q_dS_over_t_X])+webFrontSparCell2.q_B2_X329,13314
88
B1_Z = sum([webTop.q_dS_over_t_Z])+webRearSpar.q_dS_over_t_Z+ sum([webBottom.q_dS_over_t_Z])+webB1_Z330,13423
89
B2_Z = sum([webLowerNose.q_dS_over_t_Z])+ sum([webUpperNose.q_dS_over_t_Z])+webFrontSparCell2.q_B2_Z331,13544
90
Amat = [A11 A12; A21 A22];Amat333,13654
91
Bmat_X = -[B1_X;B2_X];Bmat_X334,13681
92
Bmat_Z = -[B1_Z;B2_Z];Bmat_Z335,13704
93
qs_X = inv(Amat)*Bmat_X;qs_X337,13728
94
qs_Z = inv(Amat)*Bmat_Z;qs_Z338,13753
95
sum_2_a_q_X = sum([webTop.two_A_qprime_X])+webRearSpar.two_A_qprime_X+ sum([webBottom.two_A_qprisum_2_a_q_X342,13781
96
sum_2_a_q_X = sum_2_a_q_X + sum([webLowerNose.two_A_qprime_X])+ sum([webUpperNose.two_A_qprime_Xsum_2_a_q_X343,13902
97
sum_2_a_q_X = sum_2_a_q_X + 2*qs_X(1)*(sum([webTop.Area])+webRearSpar.Area+ sum([webBottom.Areasum_2_a_q_X344,14020
98
sum_2_a_q_X = sum_2_a_q_X + 2*qs_X(2)*(sum([webLowerNose.Area])+ sum([webUpperNose.Area]));sum_2_a_q_X345,14121
99
sum_2_a_q_Z = sum([webTop.two_A_qprime_Z])+webRearSpar.two_A_qprime_Z+ sum([webBottom.two_A_qprisum_2_a_q_Z347,14215
100
sum_2_a_q_Z = sum_2_a_q_Z + sum([webLowerNose.two_A_qprime_Z])+ sum([webUpperNose.two_A_qprime_Zsum_2_a_q_Z348,14336
101
sum_2_a_q_Z = sum_2_a_q_Z + 2*qs_Z(1)*(sum([webTop.Area])+webRearSpar.Area+ sum([webBottom.Areasum_2_a_q_Z349,14454
102
sum_2_a_q_Z = sum_2_a_q_Z + 2*qs_Z(2)*(sum([webLowerNose.Area])+ sum([webUpperNose.Area]));sum_2_a_q_Z350,14555
103
torque_Z = Vz*(sc.posX - 0.25*chord);torque_Z360,14875
104
torque_X = -Vx*sc.posZ;torque_X361,14913
105
Area1 = sum([webTop.Area]) + webRearSpar.Area + sum([webBottom.Area]);Area1364,14939
106
Area1_check = get_int(frontSpar,backSpar,1)*chord^2 + get_int(frontSpar,backSpar,0)*chord^2;Area1_check366,15022
107
Area2 = sum([webLowerNose.Area]) + sum([webUpperNose.Area]);Area2368,15116
108
Area2_check = get_int(0,frontSpar,1)*chord^2 + get_int(0,frontSpar,0)*chord^2;Area2_check369,15177
109
q1t_over_q2t = (A22/Area2 + webFrontSpar.dS_over_t/Area1)/(A11/Area1 + webFrontSpar.dS_over_t/Arq1t_over_q2t374,15312
110
q2t = torque_X/(2*Area1*q1t_over_q2t + 2*Area2);q2t376,15415
111
q1t = q2t*q1t_over_q2t;q1t377,15464
112
qt_X = [q1t;q2t];qt_X378,15488
113
q2t = torque_Z/(2*Area1*q1t_over_q2t + 2*Area2);q2t380,15507
114
q1t = q2t*q1t_over_q2t;q1t381,15556
115
qt_Z = [q1t;q2t];qt_Z382,15580
116
xChord = 0:.01:1;xChord403,15777
117
xChord = xChord*chord;xChord404,15795
118
upperSurface = zeros(1,length(xChord));upperSurface405,15818
119
lowerSurface = zeros(1,length(xChord));lowerSurface406,15858
120
vecX = [frontSpar*chord webTop(i).xStart webTop(i).xEnd];vecX421,16255
121
vecZ = [0 webTop(i).zStart webTop(i).zEnd];vecZ422,16316
122
vecX = [frontSpar*chord webBottom(i).xStart webBottom(i).xEnd];vecX427,16429
123
vecZ = [0 webBottom(i).zStart webBottom(i).zEnd];vecZ428,16496
124
vecX = [frontSpar*chord webUpperNose(i).xStart webUpperNose(i).xEnd];vecX433,16618
125
vecZ = [0 webUpperNose(i).zStart webUpperNose(i).zEnd];vecZ434,16691
126
vecX = [frontSpar*chord webLowerNose(i).xStart webLowerNose(i).xEnd];vecX439,16819
127
vecZ = [0 webLowerNose(i).zStart webLowerNose(i).zEnd];vecZ440,16892
128
vecX = [frontSpar*chord sparCaps(3).posX sparCaps(4).posX];vecX444,16989
129
vecZ = [0 sparCaps(3).posZ sparCaps(4).posZ];vecZ445,17052
130
sparCapSize = 18;sparCapSize449,17136
131
stringerSize = 18;stringerSize450,17154
132
133
/home/blendux/github/UCLA_MAE_154B/wing_scripts/get_z.m,481
134
function z = get_z(x,u)get_z1,0
135
M = 0.02;M9,66
136
P = 0.4;P10,76
137
T = 0.12;T11,85
138
a0 = 0.2969;a012,95
139
a1 = -0.126;a113,108
140
a2 = -0.3516;a214,121
141
a3 = 0.2843;a315,135
142
a4 = -0.1015;a416,148
143
z_camber = M/P^2*(2*P*x - x^2);z_camber19,171
144
z_camber = (M/(1-P)^2)*(1 - 2*P +2*P*x - x^2);z_camber21,212
145
z_thickness = (T/0.2)*(a0*x^.5+a1*x+a2*x^2+a3*x^3+a4*x^4);z_thickness25,301
146
z = z_camber + z_thickness;z28,369
147
z = z_camber - z_thickness;z30,406
148
149
/home/blendux/github/UCLA_MAE_154B/wing_scripts/get_dp.m,169
150
function z = get_dp(xDist,zDist,Vx,Vz,Ix,Iz,Ixz,A)get_dp1,0
151
denom = (Ix*Iz-Ixz^2);denom3,52
152
z = -A*xDist*(Ix*Vx-Ixz*Vz)/denom - A*zDist*(Iz*Vz-Ixz*Vx)/denom;z4,75
153
154
/home/blendux/github/UCLA_MAE_154B/wing_scripts/get_int.m,925
155
function z = get_int(xi,xf,u)get_int1,0
156
M = 0.02;M3,31
157
P = 0.4;P4,41
158
T = 0.12;T5,50
159
a0 = 0.2969;a06,60
160
a1 = -0.126;a17,73
161
a2 = -0.3516;a28,86
162
a3 = 0.2843;a39,100
163
a4 = -0.1015;a410,113
164
intCamb = M/P^2*(2*P*xf^2/2 - xf^3/3) - M/P^2*(2*P*xi^2/2 - xi^3/3); intCamb16,214
165
intCamb = (M/(1-P)^2)*((1 - 2*P)*xf +2*P*xf^2/2 - xf^3/3) - (M/(1-P)^2)*((1 - 2*P)*P +2*P*P^intCamb18,303
166
intCamb = intCamb + M/P^2*(2*P*P^2/2 - P^3/3) - M/P^2*(2*P*xi^2/2 - xi^3/3);intCamb19,413
167
intCamb = (M/(1-P)^2)*((1 - 2*P)*xf +2*P*xf^2/2 - xf^3/3) - (M/(1-P)^2)*((1 - 2*P)*xi +2*P*xintCamb21,499
168
intThickness = (T/0.2)*(a0*xf^1.5/1.5 + a1*xf^2/2 + a2*xf^3/3 + a3*xf^4/4 +a4*xf^5/5);intThickness27,710
169
intThickness = intThickness - (T/0.2)*(a0*xi^1.5/1.5 + a1*xi^2/2 + a2*xi^3/3 + a3*xi^4/4 +a4*xi^intThickness28,797
170
z = intCamb + intThickness; z32,964
171
z = abs(intCamb - intThickness);z34,1001
172
173
/home/blendux/github/UCLA_MAE_154B/wing_scripts/eye_beam_example.m,1004
174
length = 20; % inlength4,44
175
force = 10000; %lbsforce5,62
176
max_width = 4; % inmax_width9,105
177
min_width = 1; % inmin_width10,125
178
y_max = 4; % iny_max11,145
179
center_y = 2; % incenter_y12,161
180
M = force*length;M17,218
181
I = min_width*(2*center_y)^3/12 + 2*( max_width*(y_max-center_y)^3/12 + ...I19,237
182
sigma_max = M * y_max / I;sigma_max22,370
183
tempCoeff = force / (I * max_width);tempCoeff29,542
184
int_y_da = ((y_max+center_y)/2) * max_width*(y_max-center_y);int_y_da30,579
185
shear_1 = tempCoeff*int_y_da;shear_131,641
186
tempCoeff = force / (I * min_width);tempCoeff34,749
187
shear_2 = tempCoeff*int_y_da;shear_235,786
188
tempCoeff = force / (I * min_width);tempCoeff39,857
189
int_y_da = (center_y/2) * min_width*center_y;int_y_da40,894
190
shear_3 = shear_2+tempCoeff*int_y_da;shear_341,940
191
int_y_da_4 = force / (I * max_width)*4*(y_max^2/2 - (center_y:.1:y_max).^2/2);int_y_da_444,1027
192
int_y_da_1 = shear_2 + force / (I * min_width)*1*(center_y^2/2 - (0:.1:center_y).^2/2);int_y_da_147,1155
193
194
/home/blendux/github/UCLA_MAE_154B/wing_scripts/get_ds.m,378
195
function ds = get_ds(xi,xf,u)get_ds1,0
196
dist = 0;dist3,31
197
numSteps = 10;numSteps4,41
198
dx = (xf-xi)/numSteps;dx5,56
199
z0 = get_z(xi,u);z06,79
200
x0 = xi;x07,97
201
tempX = x0+dx;tempX9,117
202
tempZ = get_z(tempX,u);tempZ11,153
203
tempZ = 0;tempZ13,195
204
dist = dist + (dx^2+(tempZ-z0)^2)^.5;dist15,222
205
z0 = tempZ;z016,264
206
x0 = tempX;x017,280
207
208
/home/blendux/github/UCLA_MAE_154B/wing_scripts/stringersBeamExample.m,837
209
force = 8000; % lbsforce2,12
210
stringer_A = 0.5; % in^2stringer_A3,32
211
thickness = 0.04; % inthickness4,57
212
top_stringers_y = 6; % intop_stringers_y6,81
213
middle_stringers_y = 2; % inmiddle_stringers_y7,107
214
I = 2*stringer_A*top_stringers_y^2 + 2*stringer_A*middle_stringers_y^2;I9,137
215
shear_top_web = force / (I*thickness) * top_stringers_y * stringer_A;shear_top_web15,381
216
shear_middle_web = shear_top_web + (force / (I*thickness) * middle_stringers_y * stringer_A);shear_middle_web16,451
217
d_sigma = force * top_stringers_y / I; %(lbs/in^2)d_sigma36,1520
218
d_force_top = d_sigma * stringer_A;d_force_top37,1572
219
d_sigma = force * middle_stringers_y / I; %(lbs/in^2)d_sigma40,1631
220
d_force_middle = d_force_top + d_sigma*stringer_A;d_force_middle41,1686
221
check_load = 2*d_force_top*4 + d_force_middle*4;check_load44,1762
222
223
/home/blendux/github/UCLA_MAE_154B/wing_scripts/wingAnalysis_190422.m,10240
224
Vx = 1; Vz = 1; My = 1; %test loads will be applied individuallyVx8,43
225
numTopStringers = 6;numTopStringers11,124
226
numBottomStringers = 8;numBottomStringers12,145
227
numNoseTopStringers = 4;numNoseTopStringers13,169
228
numNoseBottomStringers = 4;numNoseBottomStringers14,194
229
t_upper = 0.02/12;t_upper16,223
230
t_lower = 0.02/12;t_lower17,242
231
t_upper_front = 0.02/12;t_upper_front18,261
232
t_lower_front = 0.02/12;t_lower_front19,286
233
t_frontSpar = 0.04/12;t_frontSpar20,311
234
t_rearSpar = 0.04/12;t_rearSpar21,334
235
frontSpar = 0.2;frontSpar23,357
236
backSpar = 0.7;backSpar24,374
237
chord = 5;chord25,390
238
upperStringerGap = (sparCaps(3).posX - sparCaps(1).posX)/(numTopStringers + 1);upperStringerGap42,817
239
lowerStringerGap = (sparCaps(3).posX - sparCaps(1).posX)/(numBottomStringers + 1);lowerStringerGap43,897
240
upperNoseStringerGap = (sparCaps(1).posX - 0)/(numNoseTopStringers + 1);upperNoseStringerGap44,980
241
lowerNoseStringerGap = (sparCaps(1).posX - 0)/(numNoseBottomStringers + 1);lowerNoseStringerGap45,1053
242
Ix = 0; Iz = 0; Ixz = 0;Ix98,3144
243
Ix = Ix + sparCaps(i).area*(sparCaps(i).posZ-centroid.posZ)^2;Ix101,3191
244
Iz = Iz + sparCaps(i).area*(sparCaps(i).posX-centroid.posX)^2;Iz102,3258
245
Ixz = Ixz + sparCaps(i).area*(sparCaps(i).posX-centroid.posX)*(sparCaps(i).posZ-centroid.posIxz103,3325
246
Ix = Ix + topStringers(i).area*(topStringers(i).posZ-centroid.posZ)^2;Ix108,3470
247
Iz = Iz + topStringers(i).area*(topStringers(i).posX-centroid.posX)^2;Iz109,3545
248
Ixz = Ixz + topStringers(i).area*(topStringers(i).posX-centroid.posX)*(topStringers(i).posZ-Ixz110,3620
249
Ix = Ix + bottomStringers(i).area*(bottomStringers(i).posZ-centroid.posZ)^2;Ix113,3781
250
Iz = Iz + bottomStringers(i).area*(bottomStringers(i).posX-centroid.posX)^2;Iz114,3862
251
Ixz = Ixz + bottomStringers(i).area*(bottomStringers(i).posX-centroid.posX)*(bottomStringersIxz115,3943
252
Ix = Ix + noseTopStringers(i).area*(noseTopStringers(i).posZ-centroid.posZ)^2;Ix118,4116
253
Iz = Iz + noseTopStringers(i).area*(noseTopStringers(i).posX-centroid.posX)^2;Iz119,4199
254
Ixz = Ixz + noseTopStringers(i).area*(noseTopStringers(i).posX-centroid.posX)*(noseTopStringIxz120,4282
255
Ix = Ix + noseBottomStringers(i).area*(noseBottomStringers(i).posZ-centroid.posZ)^2;Ix123,4464
256
Iz = Iz + noseBottomStringers(i).area*(noseBottomStringers(i).posX-centroid.posX)^2;Iz124,4553
257
Ixz = Ixz + noseBottomStringers(i).area*(noseBottomStringers(i).posX-centroid.posX)*(noseBotIxz125,4642
258
numStringers = numTopStringers;numStringers135,4835
259
stringerGap = upperStringerGap;stringerGap136,4867
260
webThickness = t_upper;webThickness137,4899
261
tempStringers = topStringers;tempStringers138,4923
262
dx = web(i).xStart-centroid.posX; dz = web(i).zStart-centroid.posZ;dx154,5458
263
tempInt = get_int(web(i).xStart/chord,web(i).xEnd/chord,1)*chord^2; %integral of airfoil futempInt160,5814
264
triangle1 = abs( (web(i).xStart - sparCaps(1).posX)*web(i).zStart/2);triangle1161,5917
265
triangle2 = abs((web(i).xEnd - sparCaps(1).posX)*web(i).zEnd/2);triangle2162,5991
266
webTop = web;webTop176,6752
267
web = [];web177,6766
268
dx = web(i).xStart-centroid.posX; dz = web(i).zStart-centroid.posZ;dx187,6991
269
webRearSpar = web;webRearSpar207,8039
270
web = [];web208,8058
271
numStringers = numBottomStringers;numStringers212,8082
272
stringerGap = lowerStringerGap;stringerGap213,8117
273
webThickness = t_lower;webThickness214,8149
274
tempStringers = bottomStringers;tempStringers215,8173
275
dx = web(i).xStart-centroid.posX; dz = web(i).zStart-centroid.posZ;dx223,8489
276
tempInt = get_int(web(i).xEnd/chord,web(i).xStart/chord,0)*chord^2; %integral of airfoil futempInt238,9196
277
triangle2 = abs((web(i).xStart - sparCaps(1).posX)*web(i).zStart/2);triangle2239,9299
278
triangle1 = abs((web(i).xEnd - sparCaps(1).posX)*web(i).zEnd/2);triangle1240,9372
279
webBottom = web;webBottom256,10207
280
web = [];web257,10224
281
dx = web(i).xStart-centroid.posX; dz = web(i).zStart-centroid.posZ;dx267,10451
282
webFrontSpar = web;webFrontSpar285,11387
283
web = [];web286,11407
284
numStringers = numNoseBottomStringers;numStringers294,11453
285
stringerGap = lowerNoseStringerGap;stringerGap295,11492
286
webThickness = t_lower_front;webThickness296,11528
287
tempStringers = noseBottomStringers;tempStringers297,11558
288
dx = web(i).xStart-centroid.posX; dz = web(i).zStart-centroid.posZ;dx305,11878
289
tempInt = get_int(web(i).xEnd/chord,web(i).xStart/chord,0)*chord^2; %integral of airfoil futempInt320,12433
290
triangle1 = abs((web(i).xStart - sparCaps(2).posX)*web(i).zStart/2);triangle1321,12536
291
triangle2 = abs((web(i).xEnd - sparCaps(2).posX)*web(i).zEnd/2);triangle2322,12609
292
webLowerNose = web;webLowerNose338,13448
293
web = [];web339,13468
294
numStringers = numNoseTopStringers;numStringers342,13496
295
stringerGap = upperNoseStringerGap;stringerGap343,13532
296
webThickness = t_upper_front;webThickness344,13568
297
tempStringers = noseTopStringers;tempStringers345,13598
298
dx = web(i).xStart-centroid.posX; dz = web(i).zStart-centroid.posZ;dx353,13877
299
tempInt = get_int(web(i).xStart/chord,web(i).xEnd/chord,1)*chord^2; %integral of airfoil futempInt367,14536
300
triangle2 = abs((web(i).xStart - sparCaps(2).posX)*web(i).zStart/2);triangle2368,14639
301
triangle1 = abs((web(i).xEnd - sparCaps(2).posX)*web(i).zEnd/2);triangle1369,14712
302
webUpperNose = web;webUpperNose384,15474
303
web = [];web385,15494
304
dx = web(i).xStart-centroid.posX; dz = web(i).zStart-centroid.posZ;dx396,15722
305
webFrontSparCell2 = web;webFrontSparCell2414,16666
306
web = [];web415,16691
307
Fx = sum([webTop.qp_dx_X])+webRearSpar.qp_dx_X+ sum([webBottom.qp_dx_X])+webFrontSpar.qp_dx_X; Fx420,16736
308
Fx = Fx + sum([webLowerNose.qp_dx_X])+ sum([webUpperNose.qp_dx_X]); %cell 2Fx421,16840
309
Fz = sum([webTop.qp_dz_X])+webRearSpar.qp_dz_X+ sum([webBottom.qp_dz_X])+webFrontSpar.qp_dz_X; Fz423,16920
310
Fz = Fz + sum([webLowerNose.qp_dz_X])+ sum([webUpperNose.qp_dz_X]); %cell 2Fz424,17024
311
Fx = sum([webTop.qp_dx_Z])+webRearSpar.qp_dx_Z+ sum([webBottom.qp_dx_Z])+webFrontSpar.qp_dx_Z; Fx430,17139
312
Fx = Fx + sum([webLowerNose.qp_dx_Z])+ sum([webUpperNose.qp_dx_Z]); %cell 2Fx431,17243
313
Fz = sum([webTop.qp_dz_Z])+webRearSpar.qp_dz_Z+ sum([webBottom.qp_dz_Z])+webFrontSpar.qp_dz_Z; Fz433,17323
314
Fz = Fz + sum([webLowerNose.qp_dz_Z])+ sum([webUpperNose.qp_dz_Z]); %cell 2Fz434,17427
315
A11 = sum([webTop.dS_over_t])+webRearSpar.dS_over_t+ sum([webBottom.dS_over_t])+webFrontSpar.dS_A11443,17597
316
A22 = sum([webLowerNose.dS_over_t])+ sum([webUpperNose.dS_over_t])+webFrontSparCell2.dS_over_t;A22444,17701
317
A12 = -webFrontSpar.dS_over_t;A12445,17797
318
A21 = -webFrontSparCell2.dS_over_t;A21446,17828
319
B1_X = sum([webTop.q_dS_over_t_X])+webRearSpar.q_dS_over_t_X+ sum([webBottom.q_dS_over_t_X])+webB1_X448,17865
320
B2_X = sum([webLowerNose.q_dS_over_t_X])+ sum([webUpperNose.q_dS_over_t_X])+webFrontSparCell2.q_B2_X449,17986
321
B1_Z = sum([webTop.q_dS_over_t_Z])+webRearSpar.q_dS_over_t_Z+ sum([webBottom.q_dS_over_t_Z])+webB1_Z450,18095
322
B2_Z = sum([webLowerNose.q_dS_over_t_Z])+ sum([webUpperNose.q_dS_over_t_Z])+webFrontSparCell2.q_B2_Z451,18216
323
Amat = [A11 A12; A21 A22];Amat453,18326
324
Bmat_X = -[B1_X;B2_X];Bmat_X454,18353
325
Bmat_Z = -[B1_Z;B2_Z];Bmat_Z455,18376
326
qs_X = inv(Amat)*Bmat_X;qs_X457,18400
327
qs_Z = inv(Amat)*Bmat_Z;qs_Z458,18425
328
sum_2_a_q_X = sum([webTop.two_A_qprime_X])+webRearSpar.two_A_qprime_X+ sum([webBottom.two_A_qprisum_2_a_q_X462,18453
329
sum_2_a_q_X = sum_2_a_q_X + sum([webLowerNose.two_A_qprime_X])+ sum([webUpperNose.two_A_qprime_Xsum_2_a_q_X463,18574
330
sum_2_a_q_X = sum_2_a_q_X + 2*qs_X(1)*(sum([webTop.Area])+webRearSpar.Area+ sum([webBottom.Areasum_2_a_q_X464,18692
331
sum_2_a_q_X = sum_2_a_q_X + 2*qs_X(2)*(sum([webLowerNose.Area])+ sum([webUpperNose.Area]));sum_2_a_q_X465,18793
332
sum_2_a_q_Z = sum([webTop.two_A_qprime_Z])+webRearSpar.two_A_qprime_Z+ sum([webBottom.two_A_qprisum_2_a_q_Z467,18887
333
sum_2_a_q_Z = sum_2_a_q_Z + sum([webLowerNose.two_A_qprime_Z])+ sum([webUpperNose.two_A_qprime_Zsum_2_a_q_Z468,19008
334
sum_2_a_q_Z = sum_2_a_q_Z + 2*qs_Z(1)*(sum([webTop.Area])+webRearSpar.Area+ sum([webBottom.Areasum_2_a_q_Z469,19126
335
sum_2_a_q_Z = sum_2_a_q_Z + 2*qs_Z(2)*(sum([webLowerNose.Area])+ sum([webUpperNose.Area]));sum_2_a_q_Z470,19227
336
torque_Z = Vz*(sc.posX - 0.25*chord);torque_Z480,19547
337
torque_X = -Vx*sc.posZ;torque_X481,19585
338
Area1 = sum([webTop.Area]) + webRearSpar.Area + sum([webBottom.Area]);Area1484,19611
339
Area1_check = get_int(frontSpar,backSpar,1)*chord^2 + get_int(frontSpar,backSpar,0)*chord^2;Area1_check486,19694
340
Area2 = sum([webLowerNose.Area]) + sum([webUpperNose.Area]);Area2488,19788
341
Area2_check = get_int(0,frontSpar,1)*chord^2 + get_int(0,frontSpar,0)*chord^2;Area2_check489,19849
342
q1t_over_q2t = (A22/Area2 + webFrontSpar.dS_over_t/Area1)/(A11/Area1 + webFrontSpar.dS_over_t/Arq1t_over_q2t494,19984
343
q2t = torque_X/(2*Area1*q1t_over_q2t + 2*Area2);q2t496,20087
344
q1t = q2t*q1t_over_q2t;q1t497,20136
345
qt_X = [q1t;q2t];qt_X498,20160
346
q2t = torque_Z/(2*Area1*q1t_over_q2t + 2*Area2);q2t500,20179
347
q1t = q2t*q1t_over_q2t;q1t501,20228
348
qt_Z = [q1t;q2t];qt_Z502,20252
349
xChord = 0:.01:1;xChord523,20449
350
xChord = xChord*chord;xChord524,20467
351
upperSurface = zeros(1,length(xChord));upperSurface525,20490
352
lowerSurface = zeros(1,length(xChord));lowerSurface526,20530
353
vecX = [frontSpar*chord webTop(i).xStart webTop(i).xEnd];vecX541,20927
354
vecZ = [0 webTop(i).zStart webTop(i).zEnd];vecZ542,20988
355
vecX = [frontSpar*chord webBottom(i).xStart webBottom(i).xEnd];vecX547,21101
356
vecZ = [0 webBottom(i).zStart webBottom(i).zEnd];vecZ548,21168
357
vecX = [frontSpar*chord webUpperNose(i).xStart webUpperNose(i).xEnd];vecX553,21290
358
vecZ = [0 webUpperNose(i).zStart webUpperNose(i).zEnd];vecZ554,21363
359
vecX = [frontSpar*chord webLowerNose(i).xStart webLowerNose(i).xEnd];vecX559,21491
360
vecZ = [0 webLowerNose(i).zStart webLowerNose(i).zEnd];vecZ560,21564
361
vecX = [frontSpar*chord sparCaps(3).posX sparCaps(4).posX];vecX564,21661
362
vecZ = [0 sparCaps(3).posZ sparCaps(4).posZ];vecZ565,21724
363
sparCapSize = 18;sparCapSize569,21808
364
stringerSize = 18;stringerSize570,21826
365
366
/home/blendux/github/UCLA_MAE_154B/gui.py,393
367
import tkinter as tktk20,820
368
import tkinter.ttk as ttkttk21,841
369
class MainWindow(tk.Frame):MainWindow27,963
370
def __init__(self, *args, **kwargs):__init__30,1022
371
def create_window(self):create_window64,2364
372
def new_field(parent, name):new_field72,2661
373
def set_naca(name):set_naca80,2826
374
def set_chord(name):set_chord85,2894
375
def set_semi_span(name):set_semi_span90,2957
376
377
/home/blendux/github/UCLA_MAE_154B/tools/evaluator.py,1051
378
import numpy as npnp24,959
379
import matplotlib.pyplot as pltplt26,1000
380
class Evaluator:Evaluator29,1034
381
def __init__(self, airfoil):__init__32,1130
382
def __str__(self):__str__56,1932
383
def info_print(self, round):info_print59,1991
384
def info_save(self, save_path, number):info_save74,2570
385
def get_lift_rectangular(self, lift):get_lift_rectangular92,3362
386
def get_lift_elliptical(self, L_0):get_lift_elliptical96,3507
387
def get_lift_total(self):get_lift_total103,3719
388
def get_mass_distribution(self, total_mass):get_mass_distribution108,3900
389
def get_drag(self, drag):get_drag112,4047
390
def get_centroid(self):get_centroid124,4473
391
def get_inertia_terms(self):get_inertia_terms147,5377
392
def get_dx(self, component):get_dx180,6841
393
def get_dz(self, component):get_dz183,6940
394
def get_dP(self, xDist, zDist, V_x, V_z, area):get_dP186,7039
395
def analysis(self, V_x, V_z):analysis199,7491
396
def plot_geom(evaluator):plot_geom219,8423
397
def plot_lift(evaluator):plot_lift272,10194
398
399
/home/blendux/github/UCLA_MAE_154B/tools/creator.py,1569
400
import numpy as npnp31,1142
401
import bisect as bibi33,1193
402
import matplotlib.pyplot as pltplt34,1213
403
class Airfoil:Airfoil37,1247
404
chord = 100chord51,1727
405
semi_span = 200semi_span52,1743
406
def __init__(self):__init__54,1764
407
def from_dimensions(cls, chord, semi_span):from_dimensions65,2007
408
def __str__(self):__str__75,2337
409
def add_naca(self, naca_num):add_naca78,2396
410
def get_camber(x):get_camber99,3112
411
def get_thickness(x):get_thickness113,3714
412
def get_theta(x):get_theta124,4162
413
def get_upper_coord(x):get_upper_coord133,4478
414
def get_lower_coord(x):get_lower_coord138,4663
415
def add_mass(self, mass):add_mass166,5830
416
def info_print(self, round):info_print169,5886
417
def info_save(self, save_path, number):info_save184,6468
418
class Spar(Airfoil):Spar201,7201
419
def __init__(self):__init__204,7268
420
def add_coord(self, airfoil, x_loc_percent):add_coord212,7453
421
def add_spar_caps(self, spar_cap_area):add_spar_caps239,8383
422
def add_mass(self, mass):add_mass243,8486
423
def add_webs(self, thickness):add_webs247,8576
424
class Stringer(Airfoil):Stringer258,8922
425
def __init__(self):__init__261,9001
426
def add_coord(self, airfoil,add_coord270,9214
427
def add_area(self, area):add_area330,11712
428
def add_mass(self, mass):add_mass334,11788
429
def add_webs(self, thickness):add_webs338,11899
430
def info_print(self, round):info_print348,12249
431
def plot_geom(airfoil, view: False):plot_geom354,12401
432
def main():main411,14163
433
434
/home/blendux/github/UCLA_MAE_154B/tools/generator.py,332
435
def default_airfoil():default_airfoil23,867
436
class Population(creator.Airfoil):Population45,1474
437
def __init__(self, size):__init__48,1551
438
def mutate(self, prob_mt):mutate54,1720
439
def crossover(self, prob_cx):crossover57,1824
440
def reproduce(self, prob_rp):reproduce60,1923
441
def fitness():fitness63,2034
442