Insect-Inspired-Micro-Aerial-Vehicles

昆蟲啟發的微型飛行器:撲翼空氣動力學的非定常升力機制、壓電致動器設計與低雷諾數自主飛行控制。

FLAPPING FREQ: 120 Hz · Re ~ 10³ · LIFT 1.5× BODY WEIGHT
WINGSPAN3–15 cm昆蟲級MAV的典型翼展範圍
MASS< 5 g含致動器、電池與控制電路的總起飛質量
FLAP FREQ30–200 Hz從蜻蜓到蜜蜂的撲翼頻率範圍
REYNOLDS10²–10⁴低雷諾數流體力學的典型工作區間

當物理定律拒絕無人機縮小

將一架四軸無人機從30公分縮小到3公分並不是簡單的等比縮放——這根本行不通。當飛行器的尺寸縮小一個數量級時,翼面積(升力來源)以平方衰減,而體積(質量來源)以立方衰減。與此同時,空氣的性質也發生了變化:雷諾數(Re)從幾十萬跌落到幾十到幾百的範圍,空氣變得「黏稠」得像蜂蜜——流體慣性力相對於黏滯力的比值急劇減少,傳統的固定翼升力機制徹底失效。

然而,大自然早在數億年前就解決了這個問題。一隻蜜蜂在Re ≈ 2,000的條件下飛行——這個數字遠低於任何人造飛行器——使用的正是撲翼飛行的非定常空氣動力學機制:前緣渦(Leading Edge Vortex, LEV)在每次撲翼的下拍過程中在翼面上方形成一個穩定的低壓區,提供遠超定常理論預測的升力。

前緣渦:昆蟲的升力秘密

當一個薄翼以高攻角(通常在30°至45°之間)進行撲動時,氣流在前緣分離,形成一個強烈的旋渦結構——前緣渦(LEV)。在常規的定常流體力學中,這樣的分離渦應該會迅速增長並脫落,導致失速。但在昆蟲尺度,翼尖渦(tip vortex)與展向流(spanwise flow)共同作用將LEV穩定地錨定在翼面上方,使其在下拍過程中保持穩定並持續提供升力增強。

此外,昆蟲還利用撲翼循環中的其他非定常機制:上拍剛結束時翼的快速翻轉(clap-and-fling)可以在兩個翼面之間捕捉並加速空氣;尾流捕捉(wake capture)則使翼在穿過前一次撲動留下的尾流時獲得額外的能量回收。這些機制的協同作用使昆蟲在Re < 100的極低雷諾數下仍能實現升力/體重比大於1.5的優異性能。

Insect-inspired micro drone
Fig 1. 昆蟲撲翼飛行中的前緣渦視覺化示意圖:LEV在下拍過程中穩定錨定於翼面上方Source: Unsplash

壓電致動器:人工肌肉的工程實現

驅動一個3公分翼展的撲翼需要什麼樣的致動器?電動馬達在縮小到毫克級時效率急劇下降(電磁力與體積的關係不利於微型化),而壓電陶瓷(piezoelectric ceramics)——如PZT(鋯鈦酸鉛)——在微型尺度展現出優異的功率密度與響應頻率。一片50 mg的壓電雙晶片(bimorph)可以在電場驅動下產生亞毫米級的彎曲位移,並在100 Hz以上的頻率下穩定振動數百萬次而不疲勞。

壓電致動器的主要限制是其位移幅度較小(通常小於總長度的1%),需要透過精巧的位移放大機構(如四連桿機構或共振放大架構)來將微小的壓電應變轉化為足夠的翼尖振幅。目前最先進的昆蟲級MAV——如哈佛大學的RoboBee——使用兩個獨立控制的壓電致動器分別驅動左右翼,實現了基本的俯仰與偏航姿態控制。

撲翼空氣動力學模擬

以下Python程式模擬了簡化的準穩態撲翼升力計算。

FlappingWingAerodynamics.pyPYTHON 3.11
import numpy as np

class FlappingWingAerodynamics:
    """昆蟲級撲翼的準穩態空氣動力學模型"""

    def __init__(self, wing_length=0.015, chord=0.005, freq=120, amplitude=60):
        self.R = wing_length   # 翼長 (m)
        self.c = chord         # 弦長 (m)
        self.f = freq          # 撲翼頻率 (Hz)
        self.phi = np.radians(amplitude)  # 撲翼振幅 (rad)
        self.rho = 1.225       # 空氣密度 (kg/m³)
        self.mu = 1.81e-5     # 空氣動力黏度 (Pa·s)

    def reynolds(self):
        """計算平均雷諾數"""
        tip_vel = 2 * self.phi * self.R * self.f  # 翼尖平均速度
        return self.rho * tip_vel * self.c / self.mu

    def quasi_steady_lift(self, alpha=40):
        """準穩態升力計算 (包含前緣渦增強因子)"""
        alpha_rad = np.radians(alpha)
        # 前緣渦增強係數 (實驗擬合)
        C_L_LEV = 1.8 * np.sin(2 * alpha_rad)
        # 第二力矩面積
        S2 = self.c * self.R ** 2 / 2
        omega = 2 * self.phi * self.f  # 平均角速度
        return 0.5 * self.rho * C_L_LEV * S2 * omega ** 2

    def lift_to_weight(self, mass, alpha=40):
        """計算升力/體重比"""
        lift = self.quasi_steady_lift(alpha)
        return lift / (mass * 9.81)

mav = FlappingWingAerodynamics(wing_length=0.015, chord=0.005, freq=120)
print(f"雷諾數: {mav.reynolds():.0f}")
print(f"升力/體重: {mav.lift_to_weight(0.0001):.2f}x")

微型化的下一站:自主感知

讓昆蟲級MAV實現真正的自主飛行需要在機載感知方面取得突破——這台不到5克的飛行器需要同時整合視覺里程計、姿態估計與避障功能,而所有這些運算的功耗必須被壓縮在數十毫瓦以內。仿生視覺系統——如基於事件相機(event camera)的光流估計——正在為這一目標提供低延遲、低功耗的感知方案。當這些技術成熟時,昆蟲級MAV將能以群體形式執行大面積的環境監測、搜索救援與精準農業任務。

免責聲明:本文所述氣動模型為準穩態近似,未考慮尾流-翼面互動、翼變形與湍流擾動等非定常效應。