-
Notifications
You must be signed in to change notification settings - Fork 0
/
classes.py
108 lines (92 loc) · 3.47 KB
/
classes.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
class Aircraft:
mission = None
systems = {}
f_length = 50
height = 10.98
def CG(self):
total_weight = 0
numerator = 0
for name, system in self.systems.items():
total_weight += system.weight
numerator += system.weight * system.loc
return numerator / total_weight
def z_CG(self):
total_weight = 0
numerator = 0
for name, system in self.systems.items():
total_weight += system.weight
numerator += system.weight * (system.z_loc / self.height)
return numerator / total_weight
def Ixx(self):
"""
Calculate Ixx of aircraft with Class II Method (Roskam pg.121)
divide system.weight by 32.2 to convert lbf to slug
system.z_loc in ft but self.z_CG() in % aircraft height,
therefore convert self.z_CG() to ft by multiplying by self.height
:returns ixx[float] units: slug*ft^2
"""
ixx = 0
for name, system in self.systems.items():
ixx += (system.weight / 32.2) * (system.z_loc - self.z_CG() * self.height) ** 2
return ixx
def Iyy(self):
height = 10.98
f_length = 50
iyy = 0
for name, system in self.systems.items():
iyy += (system.weight / 32.2) * ((system.loc * f_length - self.CG() * f_length) ** 2 + (
system.z_loc - self.z_CG() * height) ** 2)
return iyy
def Izz(self):
height = 10.98
izz = 0
for name, system in self.systems.items():
izz += (system.weight / 32.2) * (system.loc * self.f_length - self.CG() * self.f_length) ** 2
return izz
def Izx(self):
izx = 0
for name, system in self.systems.items():
izx += (system.weight / 32.2) * (system.loc * self.f_length - self.CG() * self.f_length) + (
system.z_loc - self.z_CG() * self.height)
return izx
def CG_no_HT(self):
total_weight = 0
numerator = 0
for name, system in self.systems.items():
if name != 'HT':
total_weight += system.weight
numerator += system.weight * system.loc
return numerator / total_weight
def W_total(self):
total_weight = 0
for name, system in self.systems.items():
total_weight += system.weight
return total_weight
def W_e(self):
total_weight = 0
for name, system in self.systems.items():
if name not in ['Fuel wing + drop tank', 'Fuel fuselage', 'Payload wing', 'Gun', 'Payload Internal Bay',
'Pilot']:
total_weight += system.weight
return total_weight
def W_no_HT(self):
total_weight = 0
for name, system in self.systems.items():
if name != 'HT':
total_weight += system.weight
return total_weight
def W_f(self):
return self.systems["Fuel wing + drop tank"].weight + self.systems["Fuel fuselage"].weight
def FDGW(self):
return self.W_total() - self.W_f()
def display_systems(self):
for name, system in self.systems.items():
print(name, system.weight, system.loc)
class System:
reduction_factor = 0.9
def __init__(self, name, weight, loc, system_group, z_loc=0):
self.name = name
self.weight = weight * self.reduction_factor
self.loc = loc
self.z_loc = z_loc
system_group[self.name] = self