12、平面流动的面板法及有限翼展机翼的三维流动分析

平面流动的面板法及有限翼展机翼的三维流动分析

1. 平面流动面板法基础

在分析平面流动时，面板法是一种有效的工具。首先，我们来介绍一些关键的函数和概念。

1.1 NACA 4 位数机翼截面坐标生成

naca4z 函数用于返回描述 NACA 4 位数机翼截面的点的向量，代码如下：

function z = naca4z (pts, mptt)
    [m, p, t] = naca4pars (mptt);
    chi = linspace (0, 2*pi, pts);
    x = (1 + cos (chi)) / 2;
    xte = fsolve (@naca4thickness, 1.0);
    nte = max (2, ceil ((xte - x(1)) / (x(1) - x(2))));
    xx = linspace (xte, x(1), nte);
    x = [xx, x(2:(end-1)), fliplr(xx)];
    chi = [zeros(1, nte-1), chi, 2*pi*ones(1, nte - 1)];
    z = (complex (x, naca4meanline (x, m, p)’) + ...
         ((2 * (chi < pi) - 1) * 1i * t .* ...
          naca4thickness (x) .* ...
          exp (1i * atan (naca4slope (x, m, p))) ) );
<