B样条与曲率光滑连续性结合的风力机翼型设计方法
【技术领域】
[0001] 本发明涉及翼型设计技术领域，具体的说是一种B样条与曲率光滑连续性结合的 实现多点设计攻角情况下的风力机翼型设计方法。
【背景技术】
[0002] 目前常用风力机翼型有NREL-S系列翼型、丹麦的RIS0系列翼型、荷兰的DU系列 翼型和瑞典的FFA-W系列翼型。这些翼型廓线设计方法，无论是基于特定的翼型廓线，还是 基于函数的翼型集成数学模型表征，其出发点都是基于单点攻角（比如设计攻角为6° )情 况下，在研究其几何特性及空气动力特性的基础上进行改进，以获得性能良好的叶片翼型。 并没有考虑多点攻角情况下风力机翼型的设计。对于单点攻角情况下设计出来的翼型往往 在局部攻角范围内具有较好的气动性能，然而对于更广泛的攻角范围内其气动性能会不够 理想。所导致的结果是风力机在某一种风工况情况下，风力机捕风效率好，风工况稍微有点 变化，风力机叶片沿展向的入流速度、入流角及攻角等将发生改变，其捕风效率将会降低， 即不具备风力机工作的适用性与稳定性。而实际风轮叶片要求在更广泛的攻角范围内变化 时，叶片能够稳定运行。因此，有必要考虑多点攻角情况下翼型的设计方法与优化，不片面 追求局部攻角范围内的高气动性能，寻求在更广泛的攻角范围内其气动性能整体提高。考 虑多点攻角情况下翼型的设计与优化，其难点在于翼型气动力计算的收敛问题。即在翼型 优化过程中，在某一攻角情况下，设计翼型气动力收敛，然而在另一攻角情况下，设计翼型 气动力可能不收敛，使优化程序无法继续进行。
【发明内容】
[0003] 本发明的目的在于提供一种B样条与曲率光滑连续性结合的风力机翼型设计方 法，在建立翼型B样条函数表达式基础上，在翼型优化设计过程中，考虑翼型廓线表面曲率 光滑连续性，实现多点设计攻角情况下的翼型型线控制与参数优化。
[0004] 一种B样条与曲率光滑连续性结合的风力机翼型设计方法，包括如下步骤：
[0005] 1、翼型B样条函数设计：基于B样条曲线的翼型廓线正设计方法，采用三次均匀B样条函数；
[0006]B样条函数的表达式为：
[0008] 其中，i= 0, 1，2,…，n;P1+k为控制点；Glin(t)为η次基函数，t为结点矢量；
[0009] 式⑴为第k段η次B样条曲线段& = 0，1，"*11)，这些曲线段的全体称为11次8 样条曲线，其顶点Pji= 0, 1，···，n+m)所组成的多边形称为Β样条曲线的特征多边形；其 中，G(t)为基函数，表达式为：
[0011] 对于三次B样条曲线，其基函数表达式为：
[0013] 三次B样条函数写成矩阵的形式表示为：
[0015] 其中P。，ΡρP2,匕为4个控制点，t为B样条曲线的横坐标；本发明采用曲线通过 给定的控制点作为表达曲线的方式，其中Pc和P3为B样条曲线固定的首尾两点，PJPP2为 未知控制点；
[0016] 用式（4)来分别表示翼型上、下翼面廓线坐标，即为翼型B样条函数设计方法。为 了使翼型上、下翼面首尾两点相连且表现出光滑连续的特性，使上、下翼面B样条曲线控制 点经过首尾两个给定的点，其中翼型上翼面尾缘处端点与翼型下翼面尾缘处端点同时经过 翼型廓线坐标点（1，0)，翼型上翼面前缘处端点与翼型下翼面前缘处端点同时经过翼型廓 线固定点（〇,〇)。已知翼型上、下翼面首尾两个点，那么实际上翼型控制点只有四个，即上、 下翼面各两个。该方法只需控制四个参数点，就可变化出无穷形状的风力机翼型；
[0017] 2、翼型表面曲率光滑连续性：翼型型线曲率光滑连续性用曲率及曲率变化率来表 示；
[0020]y'为翼型纵坐标y的一阶导数；y〃为翼型纵坐标y的二阶导数；y" ^为翼型纵坐 标的三阶导数。
[0021 ] 3、多点攻角的翼型优化设计
[0022] 3.1目标函数
[0023] 在设计雷诺数为Re= 3.OX106,马赫数Ma= 0. 15的条件下，以光滑和粗糙条件 下最大升阻比作为目标函数，不以单点攻角下的升阻比最大作为目标函数，而是以攻角变 化范围在4°~9°下的升阻比加权最大作为目标函数：
[0024]f(X)=max(μμ2C，Ld) (7)
[0025] 式（7)中，μ2为运行工况在光滑与粗糙条件下的权值系数，μρμ# [0l]，
与连续攻角范围的设计攻角取值有关。
[0026] 3. 2设计变量
[0027] 选取翼型上、下翼4个控制点（翼型首尾两个固定控制点除外）共8个变量作为 翼型优化设计变量：
[0028] X=(Ρ,χ)P^,P2,x)P2,y)P',x)P'liy,P'2,x)P'2,y) ⑶
[0029] 3. 3约束条件
[0030] 为了使翼型廓线在可控制的范围内变化，将B样条曲线的控制点进行如下约束：
[0031]Xnin^X^Xnax (9)
[0032] 设计变量约束范围为：0.7彡Ρ1ιΧ< 0·9,0彡Pliy彡0· 1，0·3彡P2,x彡0.5， 0· 1 彡P2,y彡 0· 2,0· 3 彡Ρ' 1ιΧ彡 0· 5，-0· 2 彡P' 彡-ο. 1，0· 7 彡P' 2,x< 1. 0，-0· 1 彡P' 2 0. 1〇
[0033] 选取最大相对厚度为21%的三种翼型进行优化设计，设定翼型最大厚度为：
[0035] 除对最大相对厚度进行约束外，还需对翼型的最大厚度弦向位置进行约束：
[0036] 0· 24 彡L_彡 0· 35 (11)
[0037] 建立翼型曲率及曲率变化率约束不等式：
[0038] (12)
[0039] 式（12)中，Q为优化翼型第i点坐标的曲率，Cii为优化翼型第i_l点坐标的曲 率。
[0040]C，_<C，「CW_ (13)
[0041] 式（13)中，C'i为优化翼型第i点坐标的曲率变化率，C'ii为优化翼型第i-1点 坐标的曲率变化率。
[0042] 4、输出结果
[0043] 采用多目标粒子群优化程序进行求解，将该算法与RF0IL翼型分析软件耦合求解 计算翼型气动性能，对风力机翼型进行型线优化设计。
[0044] 本发明方法的有益效果在于：
[0045] 1)本发明方法考虑多点攻角情况下翼型的设计与优化，在翼型不断迭代优化过程 当中，约束翼型廓线表面曲率光滑连续性特性，可以很好解决多点攻角翼型优化设计过程 当中气动力不收敛的问题。
[0046] 2)设计出来的翼型整体气动性能具有明显的提高，从而降低叶片的弦长，减轻叶 片所需的材料；由于高的整体气动性能，使得风力机具有更广泛的适用条件。
[0047] 3)本发明方法可以推广到各种厚度的风力机翼型设计、飞机翼型设计及涡轮机叶 片型线等复杂曲线设计，具有良好的社会价值和经济效益。
【附图说明】
[0048] 图1为B样条函数控制翼型廓线。
[0049] 图2为本发明方法的流程图。
[0050] 图3为单点攻角及多点攻角优化翼型廓线。
[0051] 图4a为WQ-A210翼型与WQ-D210翼型曲率的对比曲线；
[0052] 图4b为WQ-A210翼型与WQ-D210翼型曲率变化率的对比曲线。
[0053] 图5a为WQ-A210翼型与WQ-D210翼型的升力系数的对比曲线；
[0054] 图5b为WQ-A210翼型与WQ-D210翼型的升阻比的对比曲线。
[0055] 图6a为DU93-W-210翼型与WQ-D210翼型的升力系数的对比曲线；
[0056] 图6b为DU93-W-210翼型与WQ-D210翼型的升阻比的对比曲线。
【具体实施方式】
[0057] 下面结合附图，对本发明做进一步说明。如图1-6所示，一种B样条与曲率光滑连 续性结合的风力机翼型设计方法，包括如下步骤：
[0058] 1翼型B样条函数与表面曲率光滑连续性理论
[0059] 1. 1翼型B样条函数理论
[0060] 对于风力机翼型廓线设计，以往的反设计方法是给定希望达到的压力分布以及初 始的基本翼型，通过几何和流动控制方程，逐步逼近给定的气动特性，但是这种翼型设计方 法计算量大，而且不能处理多学科优化设计问题。由于B样条曲线是依据有限个空间位置 点坐标绘制出的一条光滑曲线，再通过将复杂廓线首尾点相连就可得到封闭的曲线。因此， 本发明基于B样条曲线的翼型廓线正设计方法，采用三次均匀B样条函数，该函数只需4个 控制点即可表征一段光滑曲线，这样极大的减少了复杂曲线的控制变量，有利于风力机翼 型廓线的参数化设计。
[0061]B样条函数的一般表达方式为：
[0063] 其中，i= 0, 1，2,…，n;P1+k为控制点；Glin(t)为η次基函数，t为结点矢量；
[0064] 上式为第k段η次B样条曲线段（k= 0, 1，"·η)，这些曲线段的全体称为η次B 样条曲线，其顶APji=0，1，…，n+m)所组成的多边形称为Β样条曲线的特征多边形。其 中，G(t)为基函数，表达式为：
[0066] 由于基函数具有递推性、连续性及几何不变性等特点，使得能够较好的控制翼型 廓线变化。采用三次B样条曲线，分别通过4个顶点来控制翼型的上、下翼面，而且上、下翼 面首尾两个顶点重合。
[0067] 对于三次B样条曲线，其基函数表达式为：
[0069] 因此，三次B样条函数写成矩阵的形式表示为：
[0071] 其中P。，ΡρP2, ？3为4个控制点，t为B样条曲线的横坐标。B样条曲线由两种表 示复杂型线的方式，一种是曲线不经过给定的控制点，另外一种是曲线通过给定的控制点。 为了便于翼型廓线的优化设计，本发明采用第二种表达曲线的方式，即曲线通过给定的控 制点，其中Pc和P3为B样条曲线固定的首尾两点，PJPP2为未知控制点。
[0072] 用式（4)来分别表示翼型上、下翼面廓线坐标，即为翼型B样条函数设计方法理 论。为了使翼型上、下翼面首尾两点相连且表现出光滑连续的特性，使上、下翼面B样条曲 线控制点经过首尾两个给定的点，其中翼型上翼面尾缘处端点与翼型下翼面尾缘处端点同 时经过翼型廓线坐标点（1，0)，翼型上翼面前缘处端点与翼型下翼面前缘处端点同时经过 翼型廓线固定点（〇,〇)。已知翼型上、下翼面首尾两个点，那么实际上翼型控制点只有四个， 即上、下翼面各两个。图1即为B样条曲线控制翼型廓线图，该方法只需控制四个参数点， 就可变化出无穷形状的风力机翼型。
[0073] 1. 2翼型表面曲率光滑连续性
[0074] 相关资料已对翼型曲率光滑连续性及如何影响压力分布进行了详尽的研究。然 而，他们在翼型优化设计过程中，并未考虑翼型的曲率光滑连续性。为了解决多点设计攻角 情况下翼型气动力收敛问题，在翼型优化设计过程中，须考虑翼型表面的曲率光滑连续性。 翼型型线曲率光滑连续性通常用曲率及曲率变化率来表示：
[0077]y'为翼型纵坐标y的一阶导数；y〃为翼型纵坐标y的二阶导数；y"'为翼型纵坐 标的三阶导数。<