Kaynağa Gözat

Monitor time, fuel directly

Petra Lamborn 6 yıl önce
ebeveyn
işleme
bdc27103a0
2 değiştirilmiş dosya ile 13 ekleme ve 12 silme
  1. 7
    3
      analysis.rmd
  2. 6
    9
      monitor.py

+ 7
- 3
analysis.rmd Dosyayı Görüntüle

16
 ```
16
 ```
17
 
17
 
18
 ```{r graph}
18
 ```{r graph}
19
-baseg <- ggplot(flight, aes(x=timepoint / 10)) + xlab("Time (seconds)") + theme_classic()
19
+baseg <- ggplot(flight, aes(x=mission_time)) + xlab("Time (seconds)") + theme_classic()
20
 alt <- baseg + geom_line(aes(y=current_altitude)) + ylab("Altitude (m)") + ggtitle("Altitude")
20
 alt <- baseg + geom_line(aes(y=current_altitude)) + ylab("Altitude (m)") + ggtitle("Altitude")
21
 vspd <- baseg + geom_line(aes(y=v_speed)) + ylab("Vertical speed (m/s)") + ggtitle("Vertical speed")
21
 vspd <- baseg + geom_line(aes(y=v_speed)) + ylab("Vertical speed (m/s)") + ggtitle("Vertical speed")
22
 hspd <- baseg + geom_line(aes(y=h_speed)) + ylab("Horizontal speed (m/s)") + ggtitle("Horizontal speed")
22
 hspd <- baseg + geom_line(aes(y=h_speed)) + ylab("Horizontal speed (m/s)") + ggtitle("Horizontal speed")
23
 mass <- baseg + geom_line(aes(y=mass)) + ylab("Vessel mass (kg)") + ggtitle("Mass")
23
 mass <- baseg + geom_line(aes(y=mass)) + ylab("Vessel mass (kg)") + ggtitle("Mass")
24
-fuel <- baseg + geom_line(aes(y=fuel)) + ylab("Fuel (kg)") + ggtitle("Fuel")
24
+ec <- baseg + geom_line(aes(y=electric_charge)) + ylab("Electric charge") + ggtitle("Electric charge")
25
+lf <- baseg + geom_line(aes(y=liquid_fuel)) + ylab("Liquid fuel (litres)") + ggtitle("Liquid fuel")
26
+ox <- baseg + geom_line(aes(y=oxidizer)) + ylab("Oxidizer (litres)") + ggtitle("Oxidizer")
25
 at <- baseg + geom_line(aes(y=available_thrust)) + ylab("Available thrust (Newtons)") + ggtitle("Available thrust")
27
 at <- baseg + geom_line(aes(y=available_thrust)) + ylab("Available thrust (Newtons)") + ggtitle("Available thrust")
26
 thrust <- baseg + geom_line(aes(y=current_thrust)) + ylab("Current thrust (Newtons)") + ggtitle("Current thrust")
28
 thrust <- baseg + geom_line(aes(y=current_thrust)) + ylab("Current thrust (Newtons)") + ggtitle("Current thrust")
27
 g_force <- baseg + geom_line(aes(y=g_force)) + ylab("G-force (m/s²)") + ggtitle("G-force")
29
 g_force <- baseg + geom_line(aes(y=g_force)) + ylab("G-force (m/s²)") + ggtitle("G-force")
32
 vspd
34
 vspd
33
 hspd
35
 hspd
34
 mass
36
 mass
35
-fuel
37
+ec
38
+lf
39
+ox
36
 at
40
 at
37
 thrust
41
 thrust
38
 g_force
42
 g_force

+ 6
- 9
monitor.py Dosyayı Görüntüle

42
 button_clicked = conn.add_stream(getattr, button, 'clicked')
42
 button_clicked = conn.add_stream(getattr, button, 'clicked')
43
 
43
 
44
 f = open(datafile, 'w')
44
 f = open(datafile, 'w')
45
-f.write('timepoint, current_altitude, v_speed, h_speed, mass, fuel, available_thrust, current_thrust, g_force\n')
46
-timepoint = 0
45
+f.write('mission_time, current_altitude, v_speed, h_speed, mass, electric_charge, liquid_fuel, oxidizer, available_thrust, current_thrust, g_force\n')
47
 
46
 
48
 vessel = conn.space_center.active_vessel
47
 vessel = conn.space_center.active_vessel
49
 body = vessel.orbit.body
48
 body = vessel.orbit.body
57
     v_speed = vessel.flight(r_frame).vertical_speed
56
     v_speed = vessel.flight(r_frame).vertical_speed
58
     h_speed = vessel.flight(r_frame).horizontal_speed
57
     h_speed = vessel.flight(r_frame).horizontal_speed
59
     c_thrust = vessel.thrust
58
     c_thrust = vessel.thrust
59
+    e_charge = vessel.resources.amount('ElectricCharge')
60
+    lf = vessel.resources.amount('LiquidFuel')
61
+    ox = vessel.resources.amount('Oxidizer')
60
     g_force = vessel.flight(r_frame).g_force
62
     g_force = vessel.flight(r_frame).g_force
61
-    f.write('{0}, {1}, {2}, {3}, {4}, {5}, {6}, {7}, {8}\n'.format(timepoint,
62
-                                                              current_altitude,
63
-                                                              v_speed, h_speed,
64
-                                                              vm, vm -
65
-                                                              dry_mass, at,
66
-                                                              c_thrust, g_force))
67
-    timepoint = timepoint + 1
63
+    timept = vessel.met
64
+    f.write('{0}, {1}, {2}, {3}, {4}, {5}, {6}, {7}, {8}, {9}, {10}\n'.format(timept, current_altitude, v_speed, h_speed, vm, e_charge, lf, ox, at, c_thrust, g_force))
68
     time.sleep(0.1)
65
     time.sleep(0.1)
69
 
66
 
70
 f.close()
67
 f.close()