Tpel 2020 3013899
Tpel 2020 3013899
Abstract—Unmanned aerial vehicles are characterized by a set bushfire monitoring, agriculture, or surveillance drones [4]–[7].
of requirements, like high efficiency, resiliency, and reliability that Regardless of the mission type, the main challenges in multi-
conflict with the other main requirement of high power density
rotor UAVs design are dictated by the need to maximize payload,
aimed at minimizing the overall weight and size. This article pro-
poses a novel, modular multiphase drive for a quadrotor drone, maneuverability, and flight endurance. Payload and maneu-
realized through the integration of an axial flux permanent magnet verability are related to the size and weight of its propulsion
machine and a GaN-based power electronic converter. After an components and therefore to their power density. On the other
overview of the design process, starting from the propeller choice, side, the flight endurance depends on the stored energy and on
a brief description of the system components is presented. Focus-
the efficiency of the propulsion system, from propellers to the
ing specifically on the power electronic converter, the article then
presents a full analysis of its electrical and thermal performance. energy storage [8]. It is well known that efficiency and power
Extensive experimental tests allows to validate the predictions of density are conflicting requirements.
the design and simulation stages and demonstrated the expected As in all aeronautic applications, even if they do not carry
high power density levels. people or dangerous payload onboard, standards in terms of
Index Terms—Aerospace safety, fault tolerance, gallium nitride, safety are particularly stringent [9], [10]. So far, several studies
integrated design, motor drives, thermal analysis, thermal are reported in the literature, aiming to increase the resiliency
modeling, unmanned aerial vehicles (UAVs), variable speed drives. of these systems through appropriate control in case of failure.
For example, Vey and Lunze [11] analyze the condition of a
complete rotor loss for both a quadrotor and a hexrotor. In the
I. INTRODUCTION latter case, thanks to its redundancy, through a reconfiguration
RONES are unmanned aerial vehicles (UAVs) that can fly, of the control system, just a soft failure is experienced with a
D autonomously or remotely piloted, in open or in confined
spaces, for thousand kilometers or just for few minutes. UAVs are
decrease of maneuverability. However, for a quadrotor UAV, the
complete rotor loss can result in moderate or even catastrophic
therefore designed with very different structures and sizes, ac- failure, especially when the weight increases. In [12] and [13],
cording to the expected performances. One of the most common different approaches have been investigated to avoid crashes in
structure is the so-called multirotor UAV, commonly named also such conditions. Nevertheless, if the thrust is lost in one rotor, the
vertical take-off and landing [1], [2]. It has hoovering capabilities mission capabilities are surely compromised and the more the
and can fly in every direction, horizontally and vertically, with weight the more crucial safety aspects are. While redundancy
high maneuverability [3]. Because of their versatility, multirotor helps to prevent catastrophic failures, it is also clear that it comes
UAVs are gaining more and more interest in a variety of com- at the expense of power density.
mercial and military applications such as: delivery, ambulance, Adoption of a multiphase configuration is another approach
that can be used to increase the resiliency of electric drives since
Manuscript received December 6, 2019; revised March 18, 2020, May 24,
they are still able to operate after a phase winding or inverter leg
2020, and July 13, 2020; accepted July 25, 2020. Date of publication August 4, fault [14], [15]. While this option is given high consideration for
2020; date of current version October 30, 2020. This work was supported part more electric aircraft or fully electric aircraft applications [16],
by Infineon. Recommended for publication by Associate Editor Prof. Dian Guo
Xu. (Corresponding author: Federico Marcolini.)
the scenario is completely different in the field of multirotor
Martin Schiestl, Maurizio Incurvati, Ronald Stärz, Alejandro Secades UAVs, where the most common approach is to use a commercial
Rodríguez, and Lukas Wild are with the Department of Mechatronics, MCI power converter and machine [17]. This consideration leads to
Management Center Innsbruck, 6020 Innsbruck, Austria (e-mail: martin.
schiestl@mci.edu; maurizio.incurvati@mci.edu; ronald.staerz@mci.edu;
suggesting that the application of multiphase configurations to
alejandro.secades@mci.edu; lukas.wild@mci.edu). multirotor UAVs is still an open research field. Also in this case,
Federico Marcolini, Fabio Giulii Capponi, and Federico Caricchi are however, attention should be paid on the effects of this design
with the Department of Astronautical, Electrical and Energy Engineer-
ing, University of Rome “La Sapienza,” 00185 Roma, Italy (e-mail:
choice on the overall drive weight and size.
federico.marcolini@uniroma1.it; fabio.giuliicapponi@uniroma1.it; federico. Finally, thermal design is another critical aspect to consider,
caricchi@uniroma1.it). since not only it allows to extend the flight endurance, but it
Color versions of one or more of the figures in this article are available online
at https://ieeexplore.ieee.org.
also strongly impacts lifetime and therefore reliability of com-
Digital Object Identifier 10.1109/TPEL.2020.3013899 ponents. However, increasing reliability through proper thermal
0885-8993 © 2020 IEEE. Personal use is permitted, but republication/redistribution requires IEEE permission.
See https://www.ieee.org/publications/rights/index.html for more information.
Authorized licensed use limited to: Auckland University of Technology. Downloaded on October 29,2020 at 10:39:44 UTC from IEEE Xplore. Restrictions apply.
3160 IEEE TRANSACTIONS ON POWER ELECTRONICS, VOL. 36, NO. 3, MARCH 2021
Authorized licensed use limited to: Auckland University of Technology. Downloaded on October 29,2020 at 10:39:44 UTC from IEEE Xplore. Restrictions apply.
SCHIESTL et al.: DEVELOPMENT OF A HIGH POWER DENSITY DRIVE SYSTEM FOR UNMANNED AERIAL VEHICLES 3161
Authorized licensed use limited to: Auckland University of Technology. Downloaded on October 29,2020 at 10:39:44 UTC from IEEE Xplore. Restrictions apply.
3162 IEEE TRANSACTIONS ON POWER ELECTRONICS, VOL. 36, NO. 3, MARCH 2021
TABLE II
PERFORMANCES OF THE ELECTRIC MACHINE
Authorized licensed use limited to: Auckland University of Technology. Downloaded on October 29,2020 at 10:39:44 UTC from IEEE Xplore. Restrictions apply.
SCHIESTL et al.: DEVELOPMENT OF A HIGH POWER DENSITY DRIVE SYSTEM FOR UNMANNED AERIAL VEHICLES 3163
Fig. 9. Overall control structure of the UAV with redundant central controllers
Fig. 7. Distribution of the half-bridges of the four inverter systems. and independent local controllers.
TABLE III
POWER ELECTRONIC CONVERTER RATED PERFORMANCES
D. Powertrain Integration
Fig. 10 shows a cross section of the 3-D drawing with all
Fig. 8. (a) Bottom and (b) top views of half-bridge board module. (c) Control components assembled. Starting from the right, it is possible to
electronics board mounted on dc-link PCB. notice the shaft, where the propeller is going to be mounted;
the AFPM motor, composed of two rotors and one stator in
between; the twelve half-bridge PCB modules, mounted radially
The main characteristics of the converter are summarized in in correspondence with each coil; the dc link and control boards.
Table III, while Fig. 8 shows the prototypes of one half-bridge Mounting holes are also visible, that are used to anchor the drive
module PCB and the control board directly mounted on the dc- to the drone framework, and also the bearings that are placed on
link PCB. the front shield and the housing structure.
The whole UAV is composed of four propellers with their Correct radial positioning of the half-bridge PCB modules
integrated propulsion systems and a main body equipped with is ensured through screw terminals on the dc-link board. Each
redundant central controllers (see Fig. 9). Each local Aurix module is then pressed through a metallic clip against one of the
TC233 micro-controller implements a sensorless vector control twelve aluminum spokes, which act both as a mechanical support
scheme for each of the four inverters. The sensorless algorithm and as a thermal extraction path. In fact, the clip directly presses
is based on a back-emf estimation algorithm and used hall-effect the top side of each GaN against the aluminum (see Fig. 11).
Authorized licensed use limited to: Auckland University of Technology. Downloaded on October 29,2020 at 10:39:44 UTC from IEEE Xplore. Restrictions apply.
3164 IEEE TRANSACTIONS ON POWER ELECTRONICS, VOL. 36, NO. 3, MARCH 2021
A. Losses Estimation
A first estimation of the GaN losses, and the efficiency of the
whole power electronic converter, can be obtained using (1), (2),
[34], and (3)–(5)
Electrical insulation between the source terminal of the device 1 ma
Pc,pos_half = Rdson Ipk
2
+ cos φ (1)
and the aluminum frame is guaranteed by a thermal interface 8 3π
pad. This solution is intended to minimize the thermal resistance
1 ma
from the top side of the GaN to the ambient. The clip, therefore, Pc,neg_half = Rdson Ipk
2
− cos φ (2)
allows the mechanical frame to carry out both its mechanical 8 3π
and thermal functions. where ma is the modulation index and cosφ is the displacement
Finally, the external frame is designed to provide a totally power factor
enclosed manufacturing and is equipped with fins along the
axial direction thus taking advantage of the air flow caused by Rdson Ipk
2
Pc,tot = (3)
the propeller and therefore decreasing the thermal resistance 4
between ambient and frame. Table IV lists the size and weights 1 1
of the different components and calculates the expected power Psw = Vdc Ipk (tr + tf ) fsw + Coss Vdc2
fsw (4)
π 2
densities. The power density of the overall drive, considering
the cylindrical bulk volume, is approximately 1.43 kW/l, which 2VSD Ipk I2pk
Pdead = fsw tdead + Rdsoff + ΔI 2
is a very good result if it is compared to [33], where a value π 2
of 0.71 kW/l was reached. In addition, the breakdown of the (5)
weights shows that frame and covers are still a significant part
of the overall system. It is therefore possible to further optimize where Pc,pos_half and Pc,neg_half are the conduction losses of
the geometry of the mechanical frame to reduce weight and size the switch, respectively, during the positive and the negative half
and to increase the overall power density. cycles of the sinusoidal current, Pc,tot are the average losses of
the GaN over a fundamental cycle, Psw are the switching losses
including the term related to the output capacitance, and Pdead
E. Considerations on Vibrations are the losses occurring during the dead-time. In (5), Rdsoff is the
Due to its intrinsic modularity, the proposed concept has the differential resistance and VSD is the reverse operation voltage
potential to increase the availability and safety of the overall of the channel both at VGS = 0 V: taken together they model
system. It is however worth to point out that vibrations have to the reverse drain-source characteristic and allow to calculate the
be carefully measured in real operating conditions to evaluate conduction losses during dead-time. All parameters are defined
their effect on the overall reliability. Nevertheless, the system in Table V. The current ripple ΔI can be calculated from the
has been designed to limit the effect of vibrations basing on the dc-link voltage, Vdc , the amplitude modulation ratio, ma , and
following assumptions. the machine inductance. For the selected switching frequency
1) The motor is a slot-less, coreless, type. Therefore, cogging of 200 kHz, it is estimated to be approximately 4 A. Once the
torque and consequent vibrations are not present. current ripple is known, it is also possible to estimate the losses
Authorized licensed use limited to: Auckland University of Technology. Downloaded on October 29,2020 at 10:39:44 UTC from IEEE Xplore. Restrictions apply.
SCHIESTL et al.: DEVELOPMENT OF A HIGH POWER DENSITY DRIVE SYSTEM FOR UNMANNED AERIAL VEHICLES 3165
TABLE V
CHARACTERISTICS OF COOLGAN-E.E.S. SWITCHES
TABLE VI
MOSFET VERSUS GaN COMPARISON
Authorized licensed use limited to: Auckland University of Technology. Downloaded on October 29,2020 at 10:39:44 UTC from IEEE Xplore. Restrictions apply.
3166 IEEE TRANSACTIONS ON POWER ELECTRONICS, VOL. 36, NO. 3, MARCH 2021
Fig. 13. Arrangement of HB and main dc-link capacitors and parasitics for
LTSpice simulation of the capacitors’ current.
Fig. 14. Turn-OFF transient of a bottom switch. (a) LTSpice simulation.
(b) Measured waveform (10 ns/div 10 V/div).
TABLE VII
PARASITIC INDUCTANCES AND RESISTANCES FOR THE INVERTER BOARD
vibrations, optimized current/volume ratio, possible EMI issues,
and battery ripple current.
The amount of parallel capacitors is mainly depending on
their rms current Ic,rms which can be initially estimated through
the following equation [41]:
√ √
3 3 9
IC,rms = Iph 2ma + cos φ
2 − ma . (6)
4π π 16
loop inductance (Lloop,HB in Fig. 13). The gate inductances
(LGateHOH, LGateHOL for the high-side switch and LGateLOH, A more detailed calculation of the currents present at Cdc,HB
LGateLOL for the low-side switch) have been both calculated an- can be performed through an LTspice simulation of one half
alytically according to [36] and evaluated through simulations. bridge, including Cpar,high and Cpar,low as well as Lloop,HB , as
Table VII reports the results of the calculations and the values shown in the right end of Fig. 13. These simulation shows, in
extracted from the software. The discrepancies at the paths a worst-case operating scenario, a dc-link capacitor current of
of the low-side gate can be explained by the smaller loop, the half-bridge PCB equal to ICdc,HB,RMS = 6 A which is in
which reduces the self-inductance and is not accounted for in agreement with the initial analytical estimation of (6).
the calculations. In general, however, calculated and simulated In [19], the capacitor volume as a function of switching
valued show a good agreement, ensuring their correct estimation. frequency is shown for three different capacitor technologies:
The above-mentioned values of the parasitics were used for the comparison clearly demonstrates that ceramic capacitors
transient simulations with LTspice, to more accurately recreate allow us to minimize the required volume. Therefore MLCC type
the real switching behavior of the devices. A comparison be- devices from Murata model GCM32DC72A475ME02 were se-
tween simulation and measurements of the transient behavior of lected [42].
the switches has been carried out. An example of the turn-ON According to the device datasheet, this would result in a
transient is shown in Fig. 14. As it can be observed, the simu- temperature increase of over 35 °C if only one capacitor is used,
lation, Fig. 14(a), and the measurement, Fig. 14(b), match quite which reduces to 3 °C when using five parallel capacitors, hence
well. The oscillation frequency as well as the overshoot exhibit decreasing the individual current to IC,cap,RMS = 1.2 A (see
a good agreement, which confirms the correct evaluation of the Fig. 15). As an additional benefit, having multiple capacitors in
parasitics. A small discrepancy of 5 ns can be still observed in parallel allows us to decrease the power loop parasitic inductance
the rise time, which can be explained as a nonperfect modeling Lloop,HB , resulting in smaller overshoot and less stress on the
of the GaN or of the gate drive circuit. switches as well as in a decrease of the overall losses.
Once the parasitics are determined, the optimal sizing of the
dc-link capacitors can be addressed [37]–[39]. The number and
value of the capacitors have been selected according to [40] C. Heat Management and Thermal Simulations
taking especially into account: the nominal current of the drive, The thermal energy generated by the losses of the system
the temperature increases and therefore lifetime, the need for low needs to be dissipated by the outer structure via forced air-
inductance path in the switch layout, mechanical constraints and cooling. Fig. 16 shows the thermal circuit of one inverter board
Authorized licensed use limited to: Auckland University of Technology. Downloaded on October 29,2020 at 10:39:44 UTC from IEEE Xplore. Restrictions apply.
SCHIESTL et al.: DEVELOPMENT OF A HIGH POWER DENSITY DRIVE SYSTEM FOR UNMANNED AERIAL VEHICLES 3167
TABLE VIII
THERMAL SIMULATION PARAMETERS AND RESULTS
power, even at reduced air flow speed conditions (8.1 m/s instead
Fig. 16. (a) Thermal connection on system level. (b) Equivalent thermal circuit of the nominal value of 15 m/s, to match the experimental test
from the power sources of one inverter to the medium. conditions), both the high-side and low-side switch junction
temperatures, TJHS and TJLS , are below 60 °C.
From the simulations, the equivalent thermal resistance from
and the thermal connections at the system level as well as the
junction to ambient, RJamb , can be estimated as
location of the measurement points. The main path followed by
the heat flow goes from the switch junction over its top (RJT ), RJamb = 6.8 ◦ K/W. (7)
through the thermal interface material (RTIM )and the aluminum
structure (RALU ), till the ambient air (RALUA ). The higher The analysis also allows us to evaluate the effectiveness of the
constructive solution with no plastic cover on the GaN topside.
the airflow velocity at the outer frame, the lower the thermal
resistance between the frame and the ambient air, RALUA , will In fact, from the comparison of the values of RTIM and RJT
a ratio greater of 15 is found, clearly showing that the adopted
become. Any unbalance in losses between the two switches of
GaN packaging allows for a very efficient heat extraction.
the leg will activate also a secondary path from the junction to
its bottom (RJB ), through the board copper (RCu ) to the bottom The thermal analysis, however, needs to include also the
effects of the losses in the motor, since the aluminum frame
of the other switch. The power losses of the driver, PDriver , and
is shared with the power electronics converter. The motor power
the LDO circuitry, PLDO , are two additional power loss sources
connected from the junction to their case (RJC ). losses are transferred to the aluminum frame through RALU
and then to the ambient through the same frame-to-air thermal
All inverter sections are thermally connected to each other via
resistance of converter, RALUA . Assuming 50 W of losses in
the aluminum frame both at the inside and at the outside of the
structure. the machine (from the efficiency at rated conditions, Table II)
and considering that RALU is much smaller than RALUA , it
A finite element thermal simulation of the complete structure
can be estimated that the frame temperature will increase by
is not practicable, since it is very difficult to obtain a converging
mesh due to the small airgaps and the round surfaces. Therefore, 13 °C with respect to the value indicated in Table VIII. All other
temperatures will also increase by the same amount.
exploiting axial symmetry, only a 30° section is simulated,
Finally, a transient thermal simulation has also been per-
stretching it to yield a full orthogonal structure. The half-bridge
module PCB with the complete trace layout is exported into formed, to evaluate the time constant for the junction and yield-
ing a value of τ = 100 s.
Ansys Icepack software, to properly simulate the heat transfer
on the board.
IV. EXPERIMENTAL TESTS
Table VIII lists all important simulation parameters, together
with the temperatures of relevant spots resulting from the sim- The measurement of the converter power losses and so its
ulation. Additionally, Fig. 17 shows the 3-D temperature map efficiency and maximum output power evaluation has several
seen from two different angles. Results show that at the nominal challenges. First, the actual cooling performance of the fins on
Authorized licensed use limited to: Auckland University of Technology. Downloaded on October 29,2020 at 10:39:44 UTC from IEEE Xplore. Restrictions apply.
3168 IEEE TRANSACTIONS ON POWER ELECTRONICS, VOL. 36, NO. 3, MARCH 2021
R = Ploss −1 ΔT (8)
where Ploss is the power loss vector containing each loss com-
ponent and ΔT the temperature difference with respect to the
ambient temperature.
Once the dc-calibration procedure is completed and the ther-
mal resistance matrix is obtained, (8) is inverted to calculate
losses from temperature measurements in the desired operation
conditions.
To quantify the error introduced by this method on the evalua-
tion of efficiency and power, both type A and type B uncertainties
are accounted for. This results in an overall uncertainty on the
temperature δT, which can be calculated through the following
equation:
δT (T ) = σlin
2 + σ2
rep (T − Tamb ) + σT C−08 (T )
2 (9)
whereas σ rep and σlin are the errors due to reproducibility and Fig. 19. (a) Phase current and (b) voltage difference between phase and
negative terminal of dc-link waveforms on RL load at rated frequency. The
linearity, respectively, while σ TC-08 is the error introduced by voltage waveform is digitally filtered.
the thermocouples and the data logger.
The relative error on losses calculation, δP/P, depends on the While the electric motor is mounted in the integrated drive, it
relative uncertainty of temperature measurement at the nominal cannot be loaded due to the absence of the propeller. Therefore,
conditions, δT/T, and depends on the relative error of the thermal tests on the converter have been carried out by using four
resistance matrix, δR/R. Moreover, since the matrix R needs now independent three-phase RL loads with isolated star points. The
to be inverted, there is an error propagation phenomena, which heatsink mounted resistors and inductors are clearly visible in
is taken into account with the condition number of the thermal Fig. 18.
resistance matrix, κ(R), yielding to The inductance value was chosen to be 4.2 μH, to match
the machine inductance. The selected resistance, instead, was
δP δT δR
≤ κ (R) + . (10) 0.375 Ω. In this way, at the rated current, a voltage drop is pro-
P T R
duced that is close to the machine’s back-emf at rated speed plus
the resistive voltage drop on the machine’s resistance. As a con-
A. Measurement Setup sequence, the required phase voltage and the modulation index
Since at the present stage, the propeller prototype is not ma are very similar to the expected ones at rated conditions when
available yet, the drive is installed into a wind tunnel (see Fig. 18) the converter is feeding the motor. Fig. 19(a) shows the phase
to recreate the airflow which would be present if the propeller current waveform at 1250 Hz and 15 Arms. For the same condi-
was mounted on the device. Due to setup limitations, the airflow tions, Fig. 19(b) shows the corresponding voltage difference be-
velocity is set to 8.1 m/s, which is lower than the nominal axial tween a phase and the negative terminal of the dc-link; the wave-
velocity (15 m/s in flight mode). form is digitally lowpass filtered (fpass = 40 kHz, fstop = 80 kHz,
Authorized licensed use limited to: Auckland University of Technology. Downloaded on October 29,2020 at 10:39:44 UTC from IEEE Xplore. Restrictions apply.
SCHIESTL et al.: DEVELOPMENT OF A HIGH POWER DENSITY DRIVE SYSTEM FOR UNMANNED AERIAL VEHICLES 3169
TABLE IX
THERMAL MEASUREMENTS
Fig. 20. Efficiency and output power for the total system versus load current
rms.
Astop = −80 dB) to highlight the fundamental. SVPWM is
adopted, consistently to what explained in Section II-B; the
control loop is executed every 50 μs (20 kHz). The third colum in Table IX shows the temperature measure-
ments taken during the experimental tests at nominal operating
B. Power Loss Measurement Procedure conditions. Comparing the simulated results in Table VIII to the
measured temperatures, it can be seen that Tframe and TAlu show
Tests have been performed by placing six thermocouples a good matching, while TCuGnd is off by about 5 °K. However,
around the system: two on the ground plane of phase 1 and this result is reasonable, in consideration of the simplifications
9 inverter legs, two on the aluminum frame near to phase 5 that were made on the geometry.
and 9, one inside the system on the aluminum ribs of phase 9 Fig. 20 illustrates the efficiency and output power of the
and the last one on the dc-link. Figs. 16 and 17 show the positions whole system as a function of the rms of the phase output
of the above thermocouples. current, Irms . The error bars in the efficiency plot are the results
The dc thermal calibration is performed by imposing a posi- of the uncertainties in power calculations, (10), deriving from
tive current between the negative and the positive dc-link termi- the uncertainties introduced by the nonadiabatic semicalori-
nals to place all the devices in reverse conduction mode. In this metric approach. At the nominal value of 15 A, the output
way, it is possible to apply losses homogeneously over all the power Pout of the system is 1222 W and the total efficiency is
GaNs (without the need for any gate signal). Moreover, since 95.1% ± 0.77%. The latter results in a total system power loss of
the reverse conduction equivalent resistance Rdsoff is five times 63.1 W ± 9.83 W.
higher than Rdson , a much lower current is required. If the total losses of the drivers and LDOs are subtracted
From the temperatures and the applied power, the equivalent (the values for a single module are reported in Table VIII),
thermal resistance and, most notably, the relationship between then the measured value of the total system power losses is in
temperature and power can be calculated. Measurements at dif- good agreement with the estimated value of the GaN losses that
ferent power levels showed an almost completely linear relation was shown in Table V. Moreover, the nominal power stated in
between power losses and temperature increase, thus confirming Table III has been reached, leaving headroom for temperature
the validity of the nonadiabatic calorimetric approach, even if the and size as well as power density improvements.
calibration takes place at different power levels than the actual Finally, in terms of transient behavior, the measurements
operating conditions. Any remaining small deviation is, in any resulted in a system thermal time constant τ equal to 100 s,
case, still accounted for in (9) by the standard deviation of the thus confirming the prediction from the simulations. This result
linearity error σ lin . ensures that short overloads will have a low impact on the
junction temperature.
C. Results Summarizing, despite the complexity of the structure both
For the dc thermal calibration, an amount of power higher from a mechanical and an electrical point of view, the adopted
than the rated one is fed into the system. This results in higher analysis and simulation methods have shown to be in very good
temperatures, which decrease the relative error made onto the agreement with measurement results. The developed design
thermal resistance matrix. The reverse input current into the strategy is scalable and can be therefore used for a reliable design
switches is therefore set to 36 A and the reverse voltage of each of bigger systems.
switch is measured, resulting in a total power loss of 100.05 W
for the whole power electronic converter. The resulting calibra-
tion temperatures can be seen in the second column of Table IX. V. CONCLUSION
By comparing TCuGnd-Ph1, TCuGnd-Ph9, and TAlu-Ph10, it In this article, a high-power density integrated modular mul-
can be observed that the temperature rise of Phase 1 and Phase tiphase drive for the propulsion of a UAV was proposed and de-
9 with respect to the aluminum frame is approximately 25 °C. scribed. Detailed electrical and thermal analysis for the converter
Since the main focus of this analysis is the calculation of the was carried out. Finally, experimental tests were performed to
converter losses and since the temperature of the frame shows verify the predictions from the model. The key conclusions can
an homogenous behavior even at nominal conditions, it is easier be summarized as follows.
to take as a reference the temperature of the frame. Hence, the 1) It is possible to design an electric drive that combines at
thermal resistance matrix reduces to the scalar thermal resistance the same time high power density, high efficiency, and a
from the switch junction to the frame, R_JALU. resilient structure to ensure reliability.
Authorized licensed use limited to: Auckland University of Technology. Downloaded on October 29,2020 at 10:39:44 UTC from IEEE Xplore. Restrictions apply.
3170 IEEE TRANSACTIONS ON POWER ELECTRONICS, VOL. 36, NO. 3, MARCH 2021
2) To meet the expected requirements, an integrated design [20] R. Abebe et al., “Integrated motor drives: State of the art and future trends,”
approach has to be followed and the frame needs to fulfill IET Elect. Power Appl., vol. 10, no. 8, pp. 757–771, Sep. 2016.
[21] R. W. Prouty, Helicopter Performance, Stability, and Control. Melbourne,
both mechanical and heat extraction functions. FL, USA: Krieger, 2002.
3) Adoption of GaN semiconductors at high switching fre- [22] P. Pounds and R. Mahony, “Design principles of large quadrotors for
quency is required to allow at the same time high funda- practical applications,” in Proc. IEEE Int. Conf. Robot. Autom., May 2009,
pp. 3265–3270.
mental frequency and high power density. [23] S. Li, Y. Li, W. Choi, and B. Sarlioglu, “High-speed electric machines:
4) The predicted overall power density of 1.43 kW/l (almost Challenges and design considerations,” IEEE Trans. Transp. Electrific.,
double with respect to the most comparable literature vol. 2, no. 1, pp. 2–13, Jan. 2016.
[24] J. Roskam and C.T. Lan, Airplane Aerodynamic and Performance.
prototype [34]) is confirmed by experimental results. Lawrence, KS, USA: DARcorporation, 1997.
[25] C. N. Adkins and R. H. Liebeck, “Design of optimum propeller,” J.
REFERENCES Propulsion Power, vol. 10 no. 5, pp 676–682, Sep. 1994.
[26] A. Cavagnino, M. Lazzari, F. Profumo, and A. Tenconi, “A compari-
[1] A. C. Watts, V. G. Ambrosia, and E. A. Hinkley, “Unmanned aircraft son between the axial flux and the radial flux structures for PM syn-
systems in remote sensing and scientific research: Classification and chronous motors,” IEEE Trans. Ind. Appl., vol. 38, no. 6, pp. 1517–1524,
considerations of use,” Remote Sens., vol. 4, no. 6, pp. 1671–1692, Nov./Dec. 2002.
Jun. 2012. [27] F. Giulii Capponi, G. De Donato, and F. Caricchi, “Recent advances
[2] M. Hassanalian and A. Abdelkefi, “Classifications, applications, and in axial-flux permanent-magnet machine technology,” IEEE Trans. Ind.
design challenges of drones: A review,” Progress Aerosp. Sci., vol. 91, Appl., vol. 48, no. 6, pp. 2190–2205, Nov./Dec. 2012.
pp. 99–131, May 2017. [28] L. Alberti and N. Bianchi, “Theory and design of fractional-slot mul-
[3] R. Mahony, V. Kumar, and P. Corke, “Multirotor aerial vehicles: modeling, tilayer windings,” IEEE Trans. Ind. Appl., vol. 49, no. 2, pp. 841–849,
estimation, and control of quadrotor,” IEEE Robot. Autom. Mag., vol. 19, Mar./Apr. 2013.
no. 3, pp. 20–32, Sep. 2012. [29] E. Fornasiero, L. Alberti, N. Bianchi, and S. Bolognani, “Considerations
[4] R. D’Andrea, “Can drones deliver?,” IEEE Trans. Autom. Sci. Eng., vol. 11, on selecting fractional-slot nonoverlapped coil windings,” IEEE Trans.
no. 3, pp. 647–648, Jul. 2014. Ind. Appl., vol. 49, no. 3, pp. 1316–1324, May/Jun. 2013.
[5] N. Homainejad and C. Rizos, “Application of multiple categories of un- [30] N. Bianchi, M. Dai Prè, and L. Aliberti, Theory and Design of Fractional-
manned aircraft systems (UAS) in different airspaces for bushfire monitor- Slot PM Machines. Padova, Italy: CLEUP, 2007.
ing and response,” Remote Sens. Spatial Inf. Sci., vol. XL-1/W4, pp. 55–60, [31] M. J. Kamper, R. Wang, and F. G. Rossouw, “Analysis and performance
Aug. 2015. of axial flux permanent-magnet machine with air-cored nonoverlapping
[6] A. Momont, Ambulance Drone tudelft.nl. [Online]. Available: https: concentrated stator windings,” IEEE Trans. Ind. Appl., vol. 44, no. 5,
//www.tudelft.nl/en/ide/research/research-labs/applied-labs/ambulance- pp. 1495–1504, Sep./Oct. 2008.
drone/, Accessed on: Jul. 22, 2019. [32] B. Xia, J. Shen, P. C. Luk, and W. Fei, “Comparative study of air-cored
[7] J. Wilke, A drone program taking flight, 2019. [Online]. Available: https: axial-flux permanent-magnet machines with different stator winding con-
//blog.aboutamazon.com/transportation/a-drone-program-taking-flight, figurations,” IEEE Trans. Ind. Electron., vol. 62, no. 2, pp. 846–856,
Accessed on: Jul. 20, 2019. Feb. 2015.
[8] O. Gur and A. Rosen, “Optimizing electric propulsion systems for [33] M. Uğur and O. Keysan, “Multi-physics design optimisation of a GaN-
unmanned aerial vehicles,” J. Aircraft, vol. 46, no. 4, pp. 1340–1353, based integrated modular motor drive system,” J. Eng., vol. 2019, no. 17,
Jul. 2009. pp. 3900–3905, Jun. 2019.
[9] E. Petritoli, F. Leccese, and L. Ciani, “Reliability and maintenance analysis [34] J. W. Kolar, H. Ertl, and F. C. Zach, “Influence of the modulation method on
of unmanned aerial vehicles,” Sensors, vol. 18, no. 9, Sep. 2018, Art. the conduction and switching losses of a PWM converter system,” in Proc.
no. 3171. Conf. Rec. IEEE Ind. App. Soc. Annu. Meeting, Oct. 1990, pp. 502–511.
[10] M. Waibel, B. Keays, and F. Augugliaro, “Drone shows: Creative po- [35] A. R. Rivera-Ramos and M. Jimenez, “Analytical models for estimating
tential and best practices,” Verity Studios, white paper, 2017. [On- parasitic components in power electronics PCBs,” in Proc. 48th Midwest
line]. Available: http://www.taittowers.com/wp-content/uploads/2017/ Symp. Circuits Syst., Aug. 2005, vol. 2, pp. 1235–1238.
03/Drone-shows-Creative-potential-and-best-practices.pdf [36] H. A. Wheeler, “Transmission-line properties of a strip on a dielectric
[11] D. Vey and J. Lunze, “Structural reconfigurability analysis of multirotor sheet on a plane,” IEEE Trans. Microw. Theory Techn., vol. 25, no. 8,
UAVs after actuator failures,” in Proc. IEEE 54th Annu. Conf. Decis. pp. 631–647, Aug. 1977.
Control, Dec. 2015, pp. 5097–5104. [37] M. Ugur and O. Keysan, “Design of a GaN based integrated modular motor
[12] V. Lippiello, F. Ruggiero, and D. Serra, “Emergency landing for a quadro- drive,” in Proc. 13th Int. Conf. Elect. Mach., Sep. 2018, pp. 1471–1477.
tor in case of a propeller failure: A PID based approach,” in Proc. IEEE [38] J. Wang, Y. Li, and Y. Han, “Integrated modular motor drive design with
Int. Symp. Saf., Secur., Rescue Robot., Oct. 2014, pp. 1–7. GaN power FETs,” IEEE Trans. Ind Appl., vol. 51, no. 4, pp. 3198–3207,
[13] A. Merheb, H. Noura, and F. Bateman, “Emergency control of AR drone Jul./Aug. 2015.
quadrotor UAV suffering a total loss of one rotor,” IEEE/ASME Trans. [39] J. Van Damme, L. Verkroost, H. Vansompel, F. De Belie, and P. Sergeant,
Mechatronics, vol. 22, no. 2, pp. 961–971, Apr. 2017. “A holistic DC link architecture design method for multiphase integrated
[14] E. Levi, “Multiphase electric machines for variable-speed applica- modular motor drives,” in Proc. IEEE Int. Elect. Mach. Drives Conf., May
tions,” IEEE Trans. Ind. Electron., vol. 55, no. 5, pp. 1893–1909, 2019, pp. 1593–1598.
May 2008. [40] M. Ugur and O. Keysan, “DC link capacitor optimization for integrated
[15] L. De Lillo et al., “Multiphase power converter drive for fault-tolerant ma- modular motor drives,” in Proc. IEEE 26th Int. Symp. Ind. Electron., 2017,
chine development in aerospace applications,” IEEE Trans. Ind. Electron., vol. I, pp. 263–270.
vol. 57, no. 2, pp. 575–583, Feb. 2010. [41] J. W. Kolar and S. D. Round, “Analytical calculation of the RMS current
[16] 16W. Cao et al., “Overview of electric motor technologies used for stress on the DC-link capacitor of voltage-PWM converter systems,” IEEE
more electric aircraft (MEA),” IEEE Trans. Ind. Electron., vol. 59, no. 9, Proc., Elect. Power Appl., vol. 153, no. 4, pp. 535–543, Jul. 2006.
pp. 3523–3531, Sep. 2012. [42] Murata Manufacturing Co. Ltd, GCM32DC72A475ME02_ (1210,
[17] Q. Quan et al., “An analytical design-optimization method for elec- X7S:EIA, 4.7uF, DC100V) Reference Sheet. [Online]. Available:
tric propulsion systems of multicopter UAVs with desired hovering en- https://search.murata.co.jp/Ceramy/image/img/A01X/G101/ENG/
durance,” IEEE/ASME Trans. Mechatronics, vol. 24, no. 1, pp. 228–239, GCM32DC72A475ME02-01.pdf. Accessed on: May 7, 2020.
Feb. 2019. [43] H. Li, X. Li, Z. Zhang, J. Wang, L. Liu, and S. Bala, “A simple calorimetric
[18] N. R. Brown, T. M. Jahns, and R. D. Lorenz, “Power converter design technique for high-efficiency GaN inverter transistor loss measurement,”
for an integrated modular motor drive,” in Proc. IEEE Ind. Appl. Annu. in Proc. IEEE 5th Workshop Wide Bandgap Power Dev. Appl., Oct 2017,
Meeting, 2007, pp. 1322–1328. pp. 251–256.
[19] T. M. Jahns and H. Dai, “The past, present, and future of power electronics [44] M. Schiestl, A. Lösch, M. Incurvati, and R. Stärz, “Accurate losses multi-
integration technology in motor drives,” CPSS Trans. Power Electron. point non adiabatic calorimetric measurement technique for WBG power
Appl., vol. 2, no. 3, pp. 197–216, Sep. 2017. converters”, in Proc. PCIM Eur. Digit. Days, Jul. 2020, pp. 1330–1337.
Authorized licensed use limited to: Auckland University of Technology. Downloaded on October 29,2020 at 10:39:44 UTC from IEEE Xplore. Restrictions apply.
SCHIESTL et al.: DEVELOPMENT OF A HIGH POWER DENSITY DRIVE SYSTEM FOR UNMANNED AERIAL VEHICLES 3171
Martin Schiestl received the bachelor’s degree in Ronald Stärz received the master’s degree in exper-
mechatronics electrical engineering and the master’s imental physics from the University of Innsbruck,
degree in mechatronics smart technologies from the Innsbruck, Austria.
MCI, Innsbruck, Austria, in 2015 and 2017, respec- He has been the Head of the Hydraulic Engineering
tively. Laboratory, the University of Innsbruck, until 2008.
In April 2016, he began working with the Emerg- Since then, he is working at MCI, where he built up
ing Applications Laboratory, focusing on the topic the study programs in Mechatronics and the Emerging
of resonant wireless charging, amplifier design, and Applications Laboratory, which he leads since 2016.
drive systems. His work focuses on high-frequency power electron-
ics and on probes for fusion experiments.
Authorized licensed use limited to: Auckland University of Technology. Downloaded on October 29,2020 at 10:39:44 UTC from IEEE Xplore. Restrictions apply.