Figure 1.
Process flow diagram; first, the multi-rotor takes off and reaches the transportation altitude, the optimization algorithm then sweeps through search-space in the rolling horizon to find an optimum policy. The multi-rotor moves with this optimum policy for this distance step. This repeats for every distance step, until the final horizon is activated. During the final horizon, at each computation, the number of stages is reduced and the final state is given.
Figure 1.
Process flow diagram; first, the multi-rotor takes off and reaches the transportation altitude, the optimization algorithm then sweeps through search-space in the rolling horizon to find an optimum policy. The multi-rotor moves with this optimum policy for this distance step. This repeats for every distance step, until the final horizon is activated. During the final horizon, at each computation, the number of stages is reduced and the final state is given.
Figure 2.
Multi-rotor schematic, with body frame and inertial frame.
Figure 2.
Multi-rotor schematic, with body frame and inertial frame.
Figure 3.
The effects of a decrease in velocity interval that results in an exponential increase in (a) number of computations O and (b) computation time.
Figure 3.
The effects of a decrease in velocity interval that results in an exponential increase in (a) number of computations O and (b) computation time.
Figure 4.
A Dynamic programming (DP) sweep with velocity interval 2 m/s and horizon length 10 m, initial velocity 2 m/s and final velocity 0 m/s. Width of the lines and intensity in redness shows a high value of combined cost function (a) ; (b) .
Figure 4.
A Dynamic programming (DP) sweep with velocity interval 2 m/s and horizon length 10 m, initial velocity 2 m/s and final velocity 0 m/s. Width of the lines and intensity in redness shows a high value of combined cost function (a) ; (b) .
Figure 5.
Description of the complete system for simulation tests of the algorithm, including the controller block, the hex-rotor block, and the optimization algorithm block.
Figure 5.
Description of the complete system for simulation tests of the algorithm, including the controller block, the hex-rotor block, and the optimization algorithm block.
Figure 6.
The variation of with respect to the pitch angle and the forward velocity is shown using the color-bar.
Figure 6.
The variation of with respect to the pitch angle and the forward velocity is shown using the color-bar.
Figure 7.
Corrected power consumption calculation based on thrust irregularity.
Figure 7.
Corrected power consumption calculation based on thrust irregularity.
Figure 8.
Flight plots for numerical simulations performed for . (a) trajectory of the UAV; (b) velocity of the UAV; (c) pitch angle of the UAV; (d) opposing drag force of the UAV; (e) energy consumption of the UAV.
Figure 8.
Flight plots for numerical simulations performed for . (a) trajectory of the UAV; (b) velocity of the UAV; (c) pitch angle of the UAV; (d) opposing drag force of the UAV; (e) energy consumption of the UAV.
Figure 9.
Software in the loop simulation architecture.
Figure 9.
Software in the loop simulation architecture.
Figure 10.
Software in the loop Gazebo DJI F550 model along with the description of the manipulator; (a) DJI-F550 model in Gazebo simulator with PX4 flight stack simulation; (b) the multi-link free joints based arm.
Figure 10.
Software in the loop Gazebo DJI F550 model along with the description of the manipulator; (a) DJI-F550 model in Gazebo simulator with PX4 flight stack simulation; (b) the multi-link free joints based arm.
Figure 11.
Publish rates and time as obtained by the ROSTOPIC tool with a sample of two consecutive messages. (a) publish rate of the velocity command topic; (b) minimum and maximum time between two commands.
Figure 11.
Publish rates and time as obtained by the ROSTOPIC tool with a sample of two consecutive messages. (a) publish rate of the velocity command topic; (b) minimum and maximum time between two commands.
Figure 12.
Flight plots for Gazebo SITL simulations performed for . (a) trajectory of the UAV; (b) velocity of the UAV; (c) pitch angle of the UAV; (d) power consumption profile; (e) energy consumption of the UAV.
Figure 12.
Flight plots for Gazebo SITL simulations performed for . (a) trajectory of the UAV; (b) velocity of the UAV; (c) pitch angle of the UAV; (d) power consumption profile; (e) energy consumption of the UAV.
Figure 13.
A description of the hardware used in the testing of the RTDP algorithm. (a) software architecture of experimental setup; (b) DJI-F550 drone used in the experiments; (c) power source and route for various components of the drone.
Figure 13.
A description of the hardware used in the testing of the RTDP algorithm. (a) software architecture of experimental setup; (b) DJI-F550 drone used in the experiments; (c) power source and route for various components of the drone.
Figure 14.
Hover power consumption in real hardware and the SITL and Simulink experiments.
Figure 14.
Hover power consumption in real hardware and the SITL and Simulink experiments.
Figure 15.
Lab scale RTDP based aerial transportation. (a) drone01 trajectory in simulation and experiment; (b) drone01 velocity profiles in experiments and simulations; (c) drone01 energy profiles in experiments and simulations; (d) drone01 trajectory in simulation and experiments also showing DP sweep with terminal costs in first horizons and DP sweep without terminal cost in the last horizon.
Figure 15.
Lab scale RTDP based aerial transportation. (a) drone01 trajectory in simulation and experiment; (b) drone01 velocity profiles in experiments and simulations; (c) drone01 energy profiles in experiments and simulations; (d) drone01 trajectory in simulation and experiments also showing DP sweep with terminal costs in first horizons and DP sweep without terminal cost in the last horizon.
Figure 16.
Lab scale RTDP based aerial transportation. (a) drone01 trajectory in experiments for λ = 0.2, 0.5, 0.7; (b) drone01 velocity profiles in experiments for λ = 0.2, 0.5, 0.7; (c) drone01 trajectory in experiments for λ = 0.2, 0.5, 0.7.
Figure 16.
Lab scale RTDP based aerial transportation. (a) drone01 trajectory in experiments for λ = 0.2, 0.5, 0.7; (b) drone01 velocity profiles in experiments for λ = 0.2, 0.5, 0.7; (c) drone01 trajectory in experiments for λ = 0.2, 0.5, 0.7.
Figure 17.
Lab scale RTDP based aerial transportation while the goal position changed at 3.5 s. (a) drone01 trajectory in simulation and experiment; (b) drone01 velocity profiles in experiment and simulations; (c) drone01 trajectory during experiment and numerical simulation, initial and final drop location are shown.
Figure 17.
Lab scale RTDP based aerial transportation while the goal position changed at 3.5 s. (a) drone01 trajectory in simulation and experiment; (b) drone01 velocity profiles in experiment and simulations; (c) drone01 trajectory during experiment and numerical simulation, initial and final drop location are shown.
Table 1.
List of all the parameters for DJI F-550 drone model.
Table 1.
List of all the parameters for DJI F-550 drone model.
Parameter | Component | Value | Parameter | Component | Value |
---|
| Controller | 0.15 | | UAV | 0.88 m |
| Controller | 0.001 | | Controller | 0.4 |
| Controller | 0.04 | | Controller | 0.4 |
| Controller | 0.15 | Max Velocity | UAV | 12 m/s |
| Controller | 0.001 | Max Altitude | UAV | 2.2 m |
| Controller | 0.04 | Min Altitude | UAV | 1 m |
| Controller | 20 | UAV mass | UAV | 3.4 Kg |
| Controller | 0 | Payload mass | Payload | 0.586 Kg |
| Controller | 8 | Arm length | UAV | 0.27 m |
| Controller | 20 | Iz | UAV | 0.05 |
| Controller | 0 | Ix | UAV | 0.037 |
| Controller | 8 | Iy | UAV | 0.037 |
Table 2.
Hellman exponents for different locations
Table 2.
Hellman exponents for different locations
Location | Hellman Exponent |
---|
Unstable air above open water surface: | 0.08 |
Neutral air above open water surface: | 0.10 |
Unstable air above flat open coast: | 0.11 |
Neutral air above flat open coast: | 0.18 |
Stable air above open water surface: | 0.27 |
Unstable air above human inhabited areas: | 0.27 |
Neutral air alcove human inhabited areas: | 0.30 |
Stable air above flat open mast: | 0.40 |
Stable air above human inhabited areas: | 0.60 |
Table 3.
Numerical simulation results with distance interval 1 m.
Table 3.
Numerical simulation results with distance interval 1 m.
Weight-age | Time (S) | Difference (%) | Energy (kJ) | Difference (%) | Velocity Interval (m/s) |
---|
| 29.73 | % | 28.2 | % | 0.1 |
| 34.39 | 0 | 23.12 | 0 | 0.1 |
| 40.446 | % | 22.39 | % | 0.1 |
Table 4.
Numerical simulation results with distance interval 2 m.
Table 4.
Numerical simulation results with distance interval 2 m.
Weight-Age | Time (S) | Difference (%) | Energy (kJ) | Difference (%) | Velocity Interval (m/s) |
---|
| 31.19 | % | 29.3 | % | 0.1 |
| 33.46 | 0 | 26.04 | 0 | 0.1 |
| 38.58 | % | 23.16 | % | 0.1 |
Table 5.
Numerical simulation results with distance interval 3.12 m.
Table 5.
Numerical simulation results with distance interval 3.12 m.
Weight-age | Time (S) | Difference (%) | Energy (kJ) | Difference (%) | Velocity Interval (m/s) |
---|
| 29 | % | 27.6 | % | 0.1 |
| 33.2 | 0 | 24.44 | 0 | 0.1 |
| 38.5 | % | 22.89 | % | 0.1 |
| 29 | % | 27.6 | % | 0.2 |
| 33.46 | 0 | 24.59 | 0 | 0.2 |
| 38.5 | % | 22.83 | % | 0.2 |
| 29.93 | % | 28.01 | % | 0.3 |
| 32.85 | 0 | 24.41 | 0 | 0.3 |
| 38.5 | % | 22.83 | % | 0.3 |
| 29.93 | % | 28.01 | % | 0.4 |
| 32.81 | 0 | 24.45 | 0 | 0.4 |
| 38.5 | % | 22.83 | % | 0.4 |
| 29.96 | % | 28.03 | % | 0.5 |
| 33.45 | 0 | 24.81 | 0 | 0.5 |
| 38.08 | % | 22.76 | % | 0.5 |
Table 6.
Numerical simulation results with distance interval 4 m.
Table 6.
Numerical simulation results with distance interval 4 m.
Weight-age | Time (S) | Difference (%) | Energy (kJ) | Difference (%) | Velocity Interval (m/s) |
---|
| 30.89 | % | 28.14 | % | 0.1 |
| 31.63 | 0 | 26.72 | 0 | 0.1 |
| 36.5 | % | 23.63 | % | 0.1 |