1. 程式人生 > >Anderson《空氣動力學基礎》5th讀書筆記 第5記——推導二維機翼的空氣動力學系數

Anderson《空氣動力學基礎》5th讀書筆記 第5記——推導二維機翼的空氣動力學系數

機翼的受力分析圖

    我們知道,空氣對一個物體產生的升力和阻力以及力矩源於作用在整個物體上的壓力分佈和剪下力分佈,所以我們分析上圖可知(取單位展長的機翼):

    對於上表面:                        dN_{u}^{'} = -p_{u}ds_{u}cos\theta - \tau _{u}ds_{u}sin\theta

                       dA_{u}^{'} = -p_{u}ds_{u}cos\theta + \tau _{u}ds_{u}sin\theta

    同理對於下表面:

                       dN_{l}^{'} = p_{l}ds_{l}cos\theta - \tau _{l}ds_{l}sin\theta

                       dA_{l}^{'} = p_{l}ds_{l}cos\theta + \tau _{l}ds_{l}sin\theta 

    PS:N為法向力,A為切向力。

於是單位機翼上的總的法向力和切向力可表示為:

                        N^{'} = -\int_{LE}^{TE}(p_{u}cos\theta +\tau _{u}sin\theta )ds_{u}+\int_{LE}^{TE}(p_{l}cos\theta - \tau _{l}sin\theta )ds_{l}

                       A^{'} = \int_{LE}^{TE}(-p_{u}sin\theta +\tau _{u}cos\theta )ds_{u}+\int_{LE}^{TE}(p_{l}sin\theta + \tau _{l}cos\theta )ds_{l}  (LE代表前緣,TE代表後緣)

我們再推導機翼受到的力矩:

    上表面受到的微元力矩:

                       dM_{u}^{'} = (p_{u}cos\theta + \tau _{u}sin\theta )xds_{u}+(-p_{u}sin\theta +\tau _{u}cos\theta )yds_{u}

                       dM_{l}^{'} = (-p_{l}cos\theta + \tau _{l}sin\theta )xds_{l}+(p_{l}sin\theta +\tau _{l}cos\theta )yds_{l}

我們需要對力矩的方向做一些解釋,我們規定力矩方向如下圖所示:

因此我們可以知道機翼的總力矩為:

                 M_{LE}^{'} = \int_{LE}^{TE} [(p_{u}cos\theta + \tau _{u}sin\theta )x+(-p_{u}sin\theta +\tau _{u}cos\theta )y]ds_{u}+\int_{LE}^{TE}[ (-p_{l}cos\theta + \tau _{l}sin\theta )x+(p_{l}sin\theta +\tau _{l}cos\theta )y]ds_{l}

我們又知道空氣動力學系數為:

                       c_{n} \equiv \frac{N^{'}}{q_{\infty}c}                 c_{a} \equiv \frac{A^{'}}{q_{\infty}c}             c_{m} \equiv \frac{M^{'}}{q_{\infty}c^{2}}                 c_{p} \equiv \frac{P^{'}}{q_{\infty}}                 PS:其中q_{\infty } \equiv \frac{1}{2} \rho _{\infty }V_{\infty }^{2}

不知道上面這些無量綱係數是怎麼來的,別怕,請參考:量綱分析——白金漢PI定理

二維機翼的幾何關係

從上圖中,我們又可以很輕鬆地推匯出:

                     dx = ds cos\theta                  dy = -(ds sin\theta )            S = c(1)

有了以上這些公式我們就可以來推導c_{n}    c_{a}   c_{m_{LE}}   c_{l}    c_{d}

由於他們三者推導過程類似,所以我們只對c_{n}做詳細推導:

1.已知N^{'} = -\int_{LE}^{TE}(p_{u}cos\theta +\tau _{u}sin\theta )ds_{u}+\int_{LE}^{TE}(p_{l}cos\theta - \tau _{l}sin\theta )ds_{l}以及c_{n} \equiv \frac{N^{'}}{q_{\infty}c}

dx = ds cos\theta 、dy = -(ds sin\theta ) 、S = c(1) 、c_{p} \equiv \frac{P^{'}}{q_{\infty}}

2.  對c_{n} \equiv \frac{N^{'}}{q_{\infty}c} 進行代換得到 :c_{n} = \frac{-\int_{0}^{c}(p_{u} cos\theta + \tau _{u }sin\theta)/cos\theta dx+\int_{0}^{c}(p_{l} cos\theta - \tau _{l}sin\theta)/cos\theta dx}{q_{\infty }c}

3. 進行化簡:c_{n} = \frac{-\int_{0}^{c}(C_{p_{u}} + C_{f _{u }}sin\theta/cos\theta) dx+\int_{0}^{c}(C_{p_{l}} - C_{f _{l }}sin\theta/cos\theta) dx}{c}

4.又因為易推    \frac{sin\theta }{cos\theta} = -\frac{dy}{dx}

5.帶入3中式子,最終可得:

                       c_{n} = \frac{\int_{0}^{c}(-C_{p_{u}} +C_{p_{l}}) dx+\int_{0}^{c}( C_{f _{u }}\frac{dy_{u}}{dx} + C_{f _{l }}\frac{dy_{l}}{dx}) dx}{c}

同理,我們可以可以推的:

                        c_{a} = \frac{\int_{0}^{c}( C_{p _{u }}\frac{dy_{u}}{dx} - C_{p_{l }}\frac{dy_{l}}{dx}) dx+\int_{0}^{c}(C_{f_{u}} +C_{f_{l}}) dx}{c}

                        c_{m_{LE}} = \frac{\int_{0}^{c}(C_{p_{u}} -C_{p_{l}})x dx-\int_{0}^{c}( C_{f _{u }}\frac{dy_{u}}{dx} + C_{f _{l }}\frac{dy_{l}}{dx}) dx+\int_{0}^{c}( C_{p _{u }}\frac{dy_{u}}{dx} + C_{f_{u}}) y_{u}dx+\int_{0}^{c}( -C_{p_{l }}\frac{dy_{l}}{dx}+C_{f_{l}} ) dx}{c^{2}}

                        c_{l} = c_{n}cos\alpha - c_{a}sin\alpha

                        c_{d} = c_{n}sin\alpha + c_{a}cos\alpha

至此,二維機翼的空氣動力學系數基本都推導完畢了,有人問幹嘛把公式推這麼複雜,其實答案很簡單,就如我們一開始所說的空氣對一個物體產生的升力和阻力以及力矩源於作用在整個物體上的壓力分佈和剪下力分佈,記住壓力分佈和剪下力分佈是祖先,其他所有的力都是這兩個力衍生下去的,這也就是為什麼我們推導公式時要把其他空氣動力學系數儘量用壓力和剪下力系數來表示,有根才有枝嘛!