def I_x(y, n_stringers_upper, n_stringers_lower):
skins = [[upper, n_stringers_upper], [lower, n_stringers_lower]]
for n in skins:
x = np.linspace(x_values[0], x_values[1], n[1])
z = (n[0](x) - y_bar) * float(f_chord(y))
steiner = stringer_area * z**2
I_induced = np.sum(steiner)
return (A + B + C + D) * f_chord(y)**3 + I_induced
(y_upper) * n_stringers_up_3 * A_stringer +
(y_lower) * n_stringers_bm_3 * A_stringer +
y_wo_stringers * A_wo_stringers) / (A_wo_stringers + A_stringer *
(n_stringers_bm_3 + n_stringers_up_3))
# Distance from the centroid of furthest points on panels
distance_upper = 0.0633 - y_new
distance_lower = -0.0633 - y_new
distance_upper_2 = 0.0633 - y_new_2
distance_lower_2 = -0.0633 - y_new_2
distance_upper_3 = 0.0633 - y_new_3
distance_lower_3 = -0.0633 - y_new_3
d_upper = [distance_upper * f_chord(0)]
d_down = [distance_lower * f_chord(0)]
for i in nribs:
if i < 10:
d_upp = distance_upper * f_chord(i * sp)
d_d = distance_lower * f_chord(i * sp)
d_upper.append(d_upp)
d_down.append(d_d)
elif 10 <= i < 20:
d_upp = distance_upper_2 * f_chord(i * sp)
d_d = distance_lower_2 * f_chord(i * sp)
d_upper.append(d_upp)
d_down.append(d_d)
else:
d_upp = distance_upper_3 * f_chord(i * sp)
d_d = distance_lower_3 * f_chord(i * sp)
y_values = np.linspace(0, half_span, 2000)
#fig, ax = plt.subplots()
#ax.plot(y_values, shear)
#plt.show()
aoa = angle_of_attack(desired_cl)
#TORSION DIAGRAMS
momentc_4 = []
for y_value in y_values:
M = force_distribution(y_value, desired_cl,
aoa)[2] + 0.1912 * force_distribution(
y_value, desired_cl, aoa)[0] * f_chord(y_value)
momentc_4.append(M)
#Convert to a function
def conversion(moment):
M_shear_centre_f = scipy.interpolate.interp1d(y_values,
moment,
kind='cubic',
fill_value='extrapolate')
return M_shear_centre_f
M_shear_centre_f = conversion(momentc_4)
def J(y):
return I_polar * f_chord(y)**3
def lower(y):
return lng_lower * f_chord(y)
def upper(y):
return lng_upper * f_chord(y)
def aft_spar(y):
return aft * f_chord(y)
def front_spar(y):
return front * f_chord(y)