Termální management: Modelování a validace rozptylu tepla v pohonných systémech UAV

Manažment tepla v pohonných systémech UAV: propojení modelování a experimentu

Růst hustoty výkonu v pohonech UAV (motory, regulátory otáček, baterie, převodovka) zvyšuje tepelné zatížení a riziko degradace spolehlivosti. Cílem tohoto článku je představit praktický rámec řízení tepla založený na kombinaci víceúrovňového modelování a experimentálního ověření. Text se zaměřuje na multirotorová i pevná křídla a zahrnuje elektrické i hybridní konfigurace.

Typické zdroje tepla v pohonu UAV

  • Elektrický motor (BLDC/PMSM): měděné ztráty I2R, železné ztráty (hysterézní a vířivé), mechanické ztráty ložisek a ventilace.
  • Regulátor (ESC): spínací ztráty polovodičů, vodivé ztráty v MOSFET/IGBT, ztráty na šuntech a v driveru.
  • Baterie: ohmické ztráty v článcích a kolektorech, teplo z chemických procesů při vysokém C-rate, riziko tepelného úniku.
  • Převodovka a ložiska: viskózní a třecí ztráty, zvláště při dlouhodobém letu s proměnlivým zatížením.

Architektura chlazení: od zdroje k prostředí

Tepelná cesta obvykle zahrnuje: tvorba tepla → vedení (kondukce) v obalech a podložkách → rozvod (heat spreader) → předání do proudícího vzduchu (konvekce) → vyzařování. Klíčové je minimalizovat tepelné odpory v sérii a maximalizovat konvekční přenos při typických letových rychlostech 5–40 m·s−1.

Materiály a rozhraní

Materiál Teplotní vodivost λ [W·m−1·K−1] Poznámka
Hliník (6061/7075) 150–170 lehký, vhodný na chladiče a rámy
Měděné vložky 380–400 heat spreader, lokální vložky pod MOSFET
Grafit/pyrolytický grafit 200–1500 (anizotropní) ultra lehké rozváděče tepla ve fólii
TIM (pasty, podložky) 2–10 kritická je tloušťka a tlak spoje
Kompozity CFRP 5–20 (závislé na orientaci) dbejte na izolující pryskyřice a lokální přehřívání

Modelování: víceúrovňový přístup

Doporučuje se kombinovat rychlé náhradní modely pro předběžný návrh s detailními simulacemi pro finální doladění. Následují čtyři úrovně věrnosti:

  1. Lumped-parameter (RC) modely: tepelné obvody s tepelnými kapacitami C a odporem Rth. Rovnice pro uzel: C·dT/dt = Pztráty − (T − Tamb)/Rth.
  2. 1D/2D analytické aproximace: přenos tepla v deskách a žebrech, výpočet optimální výšky/rozestupu žeber při daném Reynolds a Nusseltově korelaci.
  3. CFD konvekce: simulace vnějšího omývání trupu a nasávačů, interakce proudění vrtule s chladičem, tlakové ztráty kanálů.
  4. 3D FEA (kondukce a spoje): lokální hot-spoty v ESC a vinutí, návrh heat spreaderů, validace kontaktů TIM.

Ztráty a tepelné zdroje: odhad velikosti

  • Motor: PCu = I2·R(T), kde R(T) ≈ R20[1 + α(T−20 °C)], α≈0,004 K−1. Železné ztráty rostou přibližně s n·B2.
  • ESC: Pcond ≈ I2·RDS(on), Psw ≈ 0,5·V·I·(tr+tf)·fsw·Nfet.
  • Baterie: Pbat ≈ I2·Rint; při vysokých C-rate se uplatní i entropické teplo (značka ± podle stavu nabití).

Konvekce a aerotermika: rychlé odhady

Nusseltovy korelace pro proudění vzduchu poskytují závěry pro návrh:

  • Vnější omývání desky (turbulentní): NuL ≈ 0,037·ReL0,8·Pr1/3.
  • Kanál obdélníkového průřezu: použijte hydraulický průměr Dh, Nu ≈ f(Re, Pr, L/Dh); tlaková ztráta Δp ≈ f·(L/Dh)·(ρv2/2).
  • Vliv vrtule: indukované zrychlení proudu zvyšuje přenos tepla h, ale snižuje statický tlak; optimalizujte umístění chladiče mimo odtokovou separaci.

Návrh chladiče a kanálů: metodika

  1. Stanovte Pteplo (worst-case, 95. percentil zatížení) a povolený ΔT = Tmax − Tamb.
  2. Určete potřebný tepelný odpor Rth,celk ≤ ΔT/Pteplo a rozložte jej na Rspoj + Rrozvod + Rkonvekce.
  3. Navrhněte žebra: poměr výšky k tloušťce 10–20, rozestup odpovídající cílovému Re, omezte tlakové ztráty pro malá dmychadla.
  4. Minimalizujte tloušťku TIM (< 100–200 µm), zajistěte rovinnost a upínací tlak.
  5. Validujte proudění: zda je přívodní otvor ≥ 1,2× výstupního a zda nedochází k recirkulaci.

Pokročilé prvky: heat pipes, PCM a kapalné chlazení

  • Heat pipes / vapor chambers: vysoká efektivní vodivost pro rozložení hot-spotů v ESC a pod články baterie.
  • Fázově měnitelné materiály (PCM): stabilizace teplotních špiček při náběhu zatížení; důležitá je vnesená hmotnost a cyklická životnost.
  • Uzavřené kapalné okruhy: pro hybridní/VTOL s vysokým trvalým výkonem; penalizace hmotnosti a riziko úniku vyžadují FMEA.

Modelování baterií: teplo a bezpečnost

Termoelektrický model článku spojuje elektrický ekvivalent (R0−R1||C1) s parciální diferenciální rovnicí vedení tepla. Pro pack platí síť uzlů (články, mezivrstva, separátory) s anizotropním vedením. Bezpečnostní hranice: Tprovoz 0–60 °C (typicky), Tmax článku < 55 °C při dlouhodobém zatížení, ΔT mezi články < 5 °C.

Experimentální měření: senzory a metodika

  • Teplotní senzory: termočlánky typu K (izolované), RTD Pt100/1000 pro referenci, NTC na PCB ESC.
  • IR termografie: kalibrace emisivity, matné černé nátěry, korekce odrazů v průtokových tunelech.
  • Proudění a tlak: Pitotovy trubice, anemometry, diferenční tlakoměry pro kanály.
  • Elektrické veličiny: přesná metrologie V, I, ripple, fsw, záznam telemetrie z ESC/baterie.
  • Datalogging: ≥10 Hz pro teploty, ≥1 kHz pro elektrické veličiny; synchronizace časové základny.

Testy v laboratoři a v terénu

  1. Stojan motoru: řízený krouticí moment/otáčky, simulace vrtule, mapy účinnosti a teplot vinutí/magnetů.
  2. ESC bench: stupňované zatížení, měření teplot polovodičů, validace chladiče a proudění v kanálu.
  3. Pack cycling: cykly nabíjení/vybíjení při různých C-rate a proudění, sledování ΔT mezi články.
  4. Letové testy: profil „stoupání–křížení–hovor–klesání“, větrné podmínky, sběr letové telemetrie a teplot na palubě.

Plán experimentu (DoE) a metriky

  • Faktory: rychlost proudění, zatížení [% výkonu], konfigurace žeber, poloha ESC, typ TIM.
  • Odezvy: Tmax, Rth,celk, ΔT mezi články, tepelná časová konstanta τ, účinnost η.
  • Randomizace a replikace: minimalizace driftu podmínek.
  • Analýza: ANOVA, plocha odezvy pro optimalizaci hmotnost ↔ ΔT.

Nejistota měření a kalibrace

Kombinovaná nejistota uc se stanovuje pomocí odmocniny ze součtu čtverců dílčích nejistot (senzory, proudění, emisivita). Kalibrujte termočlánky ve vodních/kombinovaných lázních, validujte IR kameru černým tělesem. Při reportování uvádějte interval spolehlivosti (například k = 2) a časové průměry vs. špičky.

Průběh případu: ESC v kanálu trupu

  1. Hypotéza: ESC překračuje 95 °C při stoupání nad 80 % výkonu.
  2. RC model: odhad Rth z ESC do vzduchu 2,5 K·W−1, C = 45 J·K−1 → τ ≈ 113 s.
  3. CFD: recirkulace v zadní části kapsy; návrh deflektoru a zvýšení výstupu o 20 %.
  4. Prototyp: přidána vapor chamber + TIM 100 µm, žebra 12 mm, rozestup 4 mm.
  5. Test: Tmax klesla na 82 °C, Δp kanálu +12 Pa (akceptovatelné), stabilita při delším stoupání.

Integrace s aerodynamikou vrtule

Vrtulový slipstream zvyšuje lokální h a současně může měnit směr proudění v kanálech. Umisťujte chladiče mimo silné gradienty rychlosti a vyhněte se tečení na trailing-edge, kde vzniká separace. Pro VTOL konverze validujte režimy hover vs. forward flight – tepelné podmínky se mohou zásadně lišit.

Hybridní pohony: motor/generátor a spalovací zdroj

  • Výfukové teplo: stínění baterií a ESC před vyzařováním, keramické fólie, ventilační clony.
  • Kapalné okruhy: spárování s chladičem oleje/kapaliny generátoru; potřeba nízkohmotnostního čerpadla a redundantních cest.
  • Vibrace: vliv na kontaktní odpor TIM a uvolnění šroubových spojů; používejte pojistné prvky.

Bezpečnost a FMEA pro tepelné události

  • Definujte kritické teploty Tkrit pro motor, ESC a baterie a implementujte softwarové limitery.
  • Detekujte rychlost nárůstu teploty dT/dt jako včasný varovný signál.
  • Navrhněte bezpečné větrání a oddělení prostoru baterie, tepelné pojistky.
  • Simulujte selhání ventilace (zanesení filtru, zakrytí otvoru) a ověřte čas do dosažení Tkrit.

Digitální dvojče řízení tepelného managementu

Propojením RC modelů, CFD/FEA a telemetrie z letu vzniká digitální dvojče, které umožňuje prediktivní řízení: adaptivní limit výkonu, předehřev baterií při nízké teplotě, dynamické nasměrování proudění (aktivní klapky) či plánování profilu letu s ohledem na tepelné rezervy.

Validace modelu: kroky a kritéria

  1. Identifikace parametrů: fit Rth a C z krokové odezvy výkonu → teplota.
  2. <