Interakce vrtulí a ramen multikoptér

Přehled problému: proč interakce vrtulí a ramen určuje výkon multikoptér

U multikoptér proudění z vrtule interaguje s blízkou strukturou – zejména s rameny a kryty motorů. Vzniká tak komplexní kombinace efektů: lokální odtržení, zesílení vírů, periodické tlaky a změny indukované rychlosti, které přímo ovlivňují tah, účinnost a hluk. Vzhledem k tomu, že typická vzdálenost hrany ramene od konce vrtule je často < 0,2–0,35R (R = poloměr vrtule), jde o zónu silné aerodynamické interakce, kde lineární aproximace selhávají a je nutná kombinace CFD simulací a validace v aerodynamickém tunelu.

Základní bezrozměrné ukazatele a definice

• Koeficient tahu: CT = T / (u03C1 n^2 D^4), kde T je tah, u03C1 hustota, n otáčky [1/s], D průměr vrtule.
• Koeficient výkonu: CP = P / (u03C1 n^3 D^5).
• Činitel indukované rychlosti: u03BB = vi / (n D) pro visení; v_i je indukovaná rychlost.
• Solídnost (solidity): u03C3 = B c / (u03C0 R) (B = počet listů, c = průměrná chord).
• Advance ratio (pokrokový poměr): J = V / (n D); pro visení je J ≈ 0, v reálném letu J ≈ 0,1–0,6.

Mechanismy interakce: jaké jevy hledat

• Clipped inflow: rameno přerušuje přítok ke kořeni listu; vzniká asymetrie zatížení a pokles lokálního C_l.
• Tip-vortex & arm impingement: špičkové víry periodicky narážejí na hrany ramene → tlakové pulzace, hlukové tóny.
• Wake re-ingestion: šíření proudu od horních vrtulí na dolní (při koaxiálních/stohovaných uspořádáních) a zpětné nasávání kolem ramen.
• Bluff-body shedding: při hranaté geometrii ramene (Re ~ 10⁵–10⁶) vzniká Kármánova ulička, která interaguje se smykovými vrstvami listů.
• Swirl & kontra-rotace: rotující proud modifikuje úhel náběhu downstream vrtule/ramene; kontra-rotace redukuje netto swirl, ale mění spektrum tlakových fluktuací.

CFD přístup: výběr modelu, domény a mřížky

• Geometrie a varianty: parametrizace ramene (profil kruh/ovál/airfoil-like), hrany (rádius/zkosení), vzdálenost od roviny rotace (z/R), vysunutí ramene od osy (r/R), kryt motoru (nacelle), ochranný kruh (duct/guard).
• Rotace listů:
  • MRF (Multiple Reference Frame): rychlá stacionární aproximace pro screening.
  • Sliding Mesh / Overset: časově přesné řešení interakce s ramenem, zachytává periodické děje (doporučené pro validaci).
• Turbulentní modely: k-ω SST pro RANS baseline; DES/DDES nebo (S)LES pro zachycení špičkových vírů a interakčního spektra. Při (S)LES cílit y⁺ ≈ 1–2 na listech a ramenech; při RANS stačí y⁺ ≈ 30–100 s wall functions, ale hrozí podmodelování odtržení.
• Doména a okrajové podmínky: hranice minimálně 5–10D od vrtule; inlet s uniformním nebo profilovaným proudem (pro dopředný let), outlet s p=0; stěny ramen a listů no-slip; rotace zadána přes RPM.
• Diskretizace a časování: pro sliding mesh zvolte krok Δθ ≈ 1–3° na krok; konvergence hodnocena na periodické stacionaritě (CT, CP) přes 3–5 otáček.
• Mřížka: hexa/prism vrstvy na listech (20–30 proudových vrstev), lokální zjemnění okolo špiček a hran ramene; celkem 30–80 mil. buněk pro detailní LES jednoho rotoru s ramenem.

Metodika nastavení: krok za krokem

1. Izolovaný rotor (baseline): vytvořte CT–RPM–P mapu v MRF; validujte proti dostupným datasheetům nebo jednoduchému ťahoměrnému stojanu.
2. Přidání ramene: nejprve stacionární MRF pro 5–10 variant trubky/profilu a polohy; vyberte top 2–3 řešení.
3. Časově přesný model: sliding mesh pro vybrané varianty; analyzujte spektrum tlaku na ramenech a listech (FFT).
4. Parametrická studie: z/R ∈ [−0,3; +0,3], r/R ∈ [0,6; 0,95], t/Ramene ∈ [0,03; 0,12], hrana = {sharp, chamfer, fillet R/t ∈ [0,2; 1,5]}.
5. Dopředný let: nastavte J = 0,2–0,4 a úhel náklonu 5–15°, sledujte posun vírů a crossflow na ramenech.

Validace v aerodynamickém tunelu: návrh experimentu

• Konfigurace: montáž rotoru a ramene v měřítku 1:1 nebo 1:2 na náboji s minimálním zasahováním držáku (sting/balance). Zajistěte stejné RPM a geometrii jako v CFD.
• Měření:
  • 6-komponentní váha (tah, krouticí moment, boční síly) → CT a CP.
  • Pitot/rake nebo PIV v rovině 0,1–0,3D pod rotory → profil indukovaných rychlostí.
  • Tlakové snímače na hranách ramene (rozteč 0,05–0,1D) → tlakové pulzace a dominantní frekvence.
  • Mikrofonové pole (akustika) → tonalita na BPF (blade-passing frequency) a jejích harmonických.
• Protokol: pro každý variant 3 opakování, 5 úrovní RPM (např. 2000–5000 rpm), logování teploty a tlaku (korekce hustoty).
• Korekce tunelu: korekce blokáže (< 5 % ideál), interference se stojanem; u otevřených testovacích sekcí sledujte recirkulaci.

Srovnání CFD a tunelových dat: metriky a tolerance

• Primární: ΔC_T a ΔC_P oproti baseline rotor-only; cíl ±3–5 % pro RANS, ±2–3 % pro LES v průměrných hodnotách.
• Spektra: porovnejte hlavní pik na BPF = B·RPM/60 a 1×–4× harmoniky; tolerance ±5 % ve frekvenci, ±2 dB v SPL pro relativní srovnání.
• Průtoky/kinematika: PIV vs. CFD – RMSE indukované rychlosti v rovině pod rotorem < 10 % maxima.
• Nejistota měření: vyhodnoťte podle ASME PTC 19.1 – kombinovaná standardní nejistota cílově < 1,5 % pro tah a < 2,5 % pro moment.

Výsledky typických studií: co se obvykle potvrdí

• Vliv vzdálenosti ramene od roviny rotace (z/R): posun ramene pod rovinu rotoru o 0,1–0,2R snižuje interakční ztráty CT o 1–3 p.b. ve srovnání s nulovou vzdáleností.
• Štíhlejší, zaoblená ramena: t/D ≤ 0,06 a fillet R/t ≥ 0,5 snižují pulzace tlaku na 0,6–0,8× oproti hranatým profilům.
• Obal před hranou ramene (leading-edge cap): redukce CP o 2–4 % při zachování CT (tj. vyšší „figure of merit“).
• Kontra-rotace sousedních vrtulí: vyrovnává swirl, zlepšuje rovnoměrnost indukce, ale může přidat sekundární píky v akustice – kompromis podle priority (účinnost vs. hluk).

Akustika: interakční hluk a jak jej modelovat

• CAA workflow: CFD (časově přesné tlaky) → akustická analýza (Ffowcs Williams–Hawkings / integrační formulace) → mapy vyzařování.
• Dominantní zdroj: rozptyl špičkových vírů na hranách ramene; sekundární odtržení na rameni při dopředném letu.
• Snížení: přiblížení ramene k „tiché zóně“ mezi listy (phase alignment), zaoblení hran, mikro-zoubky na trailing edge vrtule (serrations) – testovat podle spektra.

Doporučení pro návrh ramene a rozmístění vrtulí

• Poloha ramene: z/R = −0,1 až −0,2 (rameno mírně pod rovinou rotoru); r/R ≈ 0,75–0,85, aby se minimalizovalo impingement špičkového víru.
• Průřez ramene: eliptický nebo plytký airfoil-like (t/c ≈ 0,12–0,16), s plynulým náběhem; vyhněte se ostrým hranám a kvadratickým profilům.
• Fillet a přechody: velkorysé rádiusy na napojení na držák motoru a trup; R/t ≥ 0,5.
• Guard/duct: pokud se použije, zajistěte dostatečnou vůli (tip clearance) 1,5–3 % D a plynulý vstupní náběh; jinak hrozí pokles CT.
• Konfigurace rotace: sousední vrtule v kontra-rotaci; pro X-frame zvažte fázování listů tak, aby se BPF píky nesčítaly.

Příklad parametrického CFD porovnání (výstupní tabulky)

Checklist pro CFD–tunelový program

• Definovat baseline rotor-only (CT, CP) a testovací RPM/J.
• Parametrizovat geometrii ramene a polohu (z/R, r/R, fillet).
• Spustit MRF screening → top 3 varianty.
• Pro top varianty sliding mesh s konvergencí za 3–5 otáček.
• Navrhnout tunelový test (váha, PIV/mikrofony, tlakové snímače) s korekcemi blokáže.
• Vyhodnotit ΔCT, ΔCP, spektrum @BPF, RMSE proudění; doplnit nejistoty.
• Iterovat geometrii (fillet, leading-edge cap, z/R) a potvrdit zlepšení ≥ 1–2 p.b. v CT/CP nebo −2 dB SPL.

Typické chyby a jak se jim vyhnout

• Malá doména CFD: příliš blízké hranice falešně tlumí víry → nadhodnocení CT.
• Nedostatečný časový krok: aliasing BPF → nepřesná spektra a pulzace momentu.
• Ignorování tunelových korekcí: bez korekcí blokáže/interference jsou ΔCT v řádu procent zavádějící.
• Nezladené RPM a napájení: u elektrických pohonů se mění krouticí moment s teplotou; stabilizujte napětí a logujte teplotu vinutí.

Rozšíření: dopředný let, pozemní efekt a více rotorová interakce

• Dopředný let (J>0): boční náfuk zvyšuje asymetrii zatížení listů; rameno na návětrné straně má větší vliv – testovat v úhlech náklonu 5–15°.
• Pozemní efekt: