AeroVECTOR - The Model Rocket Simulator & Tuner

This is a Model Rocket Simulator oriented towards active stabilization. It integrates the Three Degree of Freedom Equations of Motion, allowing to tune controllers used in Model Rockets. There is pre-coded controller in the file src/control.py, one can use it or run the Flight Computer's software through Software in the Loop.

https://user-images.githubusercontent.com/65097298/128649259-8a93ab7a-ccb8-4f38-8864-71a1805a3afd.mp4

Simulation Procedure

The program computes the 3DoF Equations of Motion of the Rocket and integrates their result with the trapezoidal rule. The Aerodynamic Coefficients are calculated using the same Extended Barrowman Equations that Open Rocket uses (plus some modifications). The fins' forces are computed with interpolated wind tunnel data and Diederich's Semi-Empirical Method to accurately model the behaviour for every aspect ratio and angle of attack. The program also allows for fins separated from the body. More information can be found in the technical documentation or inside src/simulation/main_simulation.py, src/aerodynamics/rocket_functions.py or src/aerodynamics/fin_aerodynamics.py

Dependencies

Mandatory

1. tkinter
2. numpy
3. scipy
4. matplotlib
5. pandas
6. vpython

Optional

7. pyserial

Without the optional dependencies, the Software in the Loop module will not work.
If someone can make an .exe that works with VPython, please let me know.

How To

The setup of the rocket is fairly simple. However, the program is not meant for designing the rocket. Open Rocket is a more comfortable option.

Creating and Opening files

To run the program, run AeroVECTOR.py. The program will open in the file tab.

One can create a new file or open an existing one. Once a file is open, a copy can be created with the Save As button.

In case of using the Python SITL module, it is recommended to save each file in a different folder since the modules can only be loaded from the SITL Modules folder located in the directory of the .txt save file, as in the examples. The SITL Modules and SITL Modules/Complementary Modules folders are created automatically when creating a new file and, in case of using the Save as command, the contents of the original modules are also copied to the new location.

Setting up the Rocket

One must fill the required parameters of the rocket. New motors can be added in the motors folder.

• Iy Liftoff/Burnout are the wet and dry pitching moment of inertia.
• Xcg Liftoff/Burnout are the position of the wet and dry Center of Gravity.
• Xt is the position of the TVC mount. If one is using fins, the program automatically calculates the force application point.
  • All distances are measured from the tip of the nose cone.
• The Servo Resolution is the minimum angle it can rotate.
• The Max Actuator Angle is the maximum angle the actuator can move (either the TVC mount or the fin).
• The Actuator Reduction is the gear ratio between the servo and the actuator.
  • The code multiplies the output of the controller times the Actuator Reduction, and then sends that output to the SERVO (not the mount). Remember that you have to multiply the output of the controller times the Actuator Reduction in you flight computer!
• The Initial Misalignment only modifies the initial angle of the TVC mount or Active Fin.
• The Servo Velocity Compensation slows down the servo according to the load, its value is 1.45 for an SG90 without load, and 2.1 with a TVC mount. The servo class found in servo_lib.py has a test method to modify this value to fit one's servo.
• The wind is positive from right to left, and the gusts follow a Gaussian distribution.
• The effective launch rod length is the length at which the rocket can pitch somewhat freely.
• Angle of the launch rod.
• The Motor Misalignment sets the motor at an angle (mainly for active fin control).
• The roughness applies to all the components of the body and fins; therefore, a little tuning might be required to get the performances right.

Please do not leave blank entries

THE SAVE BUTTON SAVES ONLY THE CURRENT TAB, BE SURE TO CLICK IT ON ALL OF THEM

Setting up the Rocket 2 - The Empire Strikes Back

Body

To draw the rocket, one must insert the point as coordinate from the nose cone tip, diameter in that point. With the Add Point button, one adds the point written in the entry. The Delete Point button deletes the point currently selected in the combo box. To modify a point, one has to select the desired point in the combo box, click the Select Point button, write the new coordinates in the entry, and at last, click the Modify Point button.

Fins

To draw the fins, one must insert the position and chord of the root and tip, separated by comma. The wingspan is measured from the root to the tip. All dimensions are in meters.
The order is the following:

Only trapezoidal fins are modelled.

After the points are written in the entries, one can either update the stabilization or control fin. Clicking the Load "" Fins button will fill the entries with the current data of the fin. The button Reset Fin sets the fin entries to a zero-area fin.

WARNING: THE USE OF CONTROL FINS DISABLES THE TVC STABILIZATION.

Examples of detached fins: Sprint - BPS.space, and our own Roll Control System Testbed

The space between the body and the fin must be considerable, if one is unsure about a fin being attached or detached, the most conservative option is the right option.
The Angle of Attack slider allows to change the AoA at which the CP (red point) is calculated. The blue point represents the CG of the rocket. One can enable and disable the fins to quickly redraw the rocket and ensure that the CG is in the correct position.

Sliders

The Set Rocket Body tab has four sliders to simulate the forces on the rocket and test changes in the configuration before simulating. For example, if you see in the simulation that your rocket looses control at 10º angle of attack, with an actuator deflection of 20º, at 3 seconds from liftoff, you can put those parameters in the sliders in the Set Rocket Body tab and make changes to make sure you have the net moment that you need (i.e. increase the control fin size to get a restaurative moment).

• Actuator deflection: The current deflection of the actuator
• Force app point: The reduction point of all the forces in the rocket, it changes color to green if the force is applied in the opposite direction.
• Time: Takes the thrust of the motor from the thrust vs time curve for the time indicated in the slider.

3D Graphics

One can activate the 3D Graphics by clicking the checkbox. IT REQUIRES VPYTHON

• Camera Shake moves the camera based on the accelerations of the rocket.
• Hide Forces hides the force arrows.
• Variable fov decreases the fov variably, maintaining the rocket of approximately the same size during the flight.
• Hide CG hides a ball that represents the CG.
• Camera type
  • "Follow" -> Follows the rocket.
  • "Fixed" -> Looks from the ground.
  • "Follow Far" -> 2D.
• Slow Motion slows the animation -> 1 = Real time, 5 = five times slower.
• Force Scale Factor scales the forces accordingly, 1 = 1 meter/Newton.
• Fov is the field of view of the camera.

Software in the Loop

Python SITL

In case of using the Python SITL module, refer to the example. One can create functions, classes, modules, etc., or modify the code at will. The objective was to make it as similar as possible to an Arduino, but there are some differences, especially with the global variables having the prefix self. Only the Python SITL module is compatible with GNSS and sample times.

On what boards can I use this software?

It was only tested on an Arduino Nano clone, so compatibility is not ensured.
Even on the Arduino, the program did not work properly with program runtimes times smaller than 5 milliseconds.

How to set up the serial communication in Python.

To use the Software in the Loop function, one must set the Port to the one in which the board is connected, and the Baudrate to the one set in the Arduino program.
One can simulate sensor noise by filling the entries with the Noise Standard Deviation of the specified sensor.

How to set up the simulation in your Arduino.

Download and include the library and create the instance with the name you want.

At the end of void setup, start the simulation.

Replace your sensor readings with Sim.getSimData() and the name of your variables.

• The Gyroscope data is in º/s.
• The Accelerometer measures the reaction force applied to it (like the real ones), the data is in g's.
• The Altimeter measures the data in meters.
• The GNSS measures the distance in the horizontal axis from the launch pad, in m.
• The GNSS measures the horizontal velocity, in m/s.

Positive values are positive in the direction of the axes!
(Refer to the Technical Documentation)

Replace your servo.write() for:

Replace servo_command for your servo variable (in º).
The parachute variable is an int, it's normally 0 and one must make it 1 when the parachute would deploy.

**REMEMBER THAT THE DATA IS IN DEGREES, G'S AND M, AND YOU HAVE TO SEND THE SERVO COMMAND IN DEGREES AND THE PARACHUTE DEPLO