Anderson《空气动力学基础》5th读书笔记 第5记——推导二维机翼的空气动力学系数

Anderson《空气动力学基础》5th读书笔记 第5记——推导二维机翼的空气动力学系数_第1张图片 机翼的受力分析图

    我们知道,空气对一个物体产生的升力和阻力以及力矩源于作用在整个物体上的压力分布和剪切力分布,所以我们分析上图可知(取单位展长的机翼):

    对于上表面:
                       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}

我们需要对力矩的方向做一些解释,我们规定力矩方向如下图所示:

Anderson《空气动力学基础》5th读书笔记 第5记——推导二维机翼的空气动力学系数_第2张图片

 

 

因此我们可以知道机翼的总力矩为:

                 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定理

Anderson《空气动力学基础》5th读书笔记 第5记——推导二维机翼的空气动力学系数_第3张图片 二维机翼的几何关系

从上图中,我们又可以很轻松地推导出:

                     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

至此,二维机翼的空气动力学系数基本都推导完毕了,有人问干嘛把公式推这么复杂,其实答案很简单,就如我们一开始所说的空气对一个物体产生的升力和阻力以及力矩源于作用在整个物体上的压力分布和剪切力分布,记住压力分布和剪切力分布是祖先,其他所有的力都是这两个力衍生下去的,这也就是为什么我们推导公式时要把其他空气动力学系数尽量用压力和剪切力系数来表示,有根才有枝嘛!

你可能感兴趣的:(Aerodynamic)