[发明专利]一种风力机翼型与叶片外形参数一体化设计方法

摘要

本发明公开了一种风力机翼型与叶片外形参数一体化设计方法，以三维叶片泛函集成表达式为出发点，此表达式将翼型参数表达及叶片三维形状表达集成于一体，能够很好的实现翼型与叶片并行设计；然后分别建立翼型设计模块和叶片气动外形设计模块，这两个模块之间通过气动力插值实现数据的传递机制；最后，采用MATLAB语言编制并行设计策略，制定收敛策略，若满足收敛条件则输出新翼型族及叶片气动外形数据，否则，返回开始程序。本发明突破传统的风力机叶片串行设计方法，基于叶片三维集成数学表达方程，将翼型廓线与叶片气动外形参数联立起来，通过MATLAB并行计算程序，实现翼型与叶片并行优化设计方法。

主权项

1.一种风力机翼型与叶片外形参数一体化设计方法，其特征在于，包括以下步骤：步骤1：基于三维叶片泛函集成表达式，将翼型参数表达及叶片三维形状表达集成于一体，构建三维叶片集成表达式；所述三维叶片集成表达式为：其中，x为翼型横坐标，y为翼型纵坐标，z为空间三维叶片沿展向坐标；a为1/4翼型弦长；R为叶片半径，r为翼型矢径；XM为翼型气动中心展向位置；u为叶片展向位置，u∈[0,1]；c(u)为弦长分布函数；β(u)为扭角分布函数；所述三维叶片集成表达式能表征任意可能的叶片几何形状的集成：即集成叶片的空间坐标翼型分布、叶片弦长分布、叶片扭角分布；步骤2：分别建立翼型设计模块和叶片气动外形设计模块，这两个模块之间通过气动力插值实现数据的传递机制；采用泛函集成理论作为翼型设计模块，其具体实现包括以下子步骤：步骤A1：确定翼型集成公式；式中，x为翼型横坐标，y为翼型纵坐标；r为翼型的矢径；为翼型控制方程函数；a_k、b_k为系数，θ为幅角，a为1/4翼型弦长；步骤A2：确定目标函数；以各翼型最大升阻比作为翼型设计模块目标函数，以光滑条件下升阻比最大为目标函数之一：f1(x)＝max(CL/CD)；其中CL为升力系数，CD为阻力系数；步骤A3：确定设计变量；选取翼型控制方程函数的第1到第8项系数作为优化设计的变量，确定设计变量为：X1＝(a1,b1,a2,b2,a3,b3,a4,b4)；步骤A4：确定约束条件；对翼型控制方程函数前8项系数进行约束，其约束条件如表1：表1设计变量范围对翼型最大相对厚度进行约束，其最大相对厚度分别约束在40％、35％、30％、25％、21％、18％及15％；步骤3：并行设计；其具体实现包括以下子步骤：步骤3.1：编制三维叶片集成表达主程序；步骤3.2：编制两个子程序模块，分别为翼型设计模块和叶片外形设计模块；采用MATLAB自带优化算法进行子程序优化；步骤3.3：基于动量叶素理论编制气动力插值程序，将优化中的翼型气动力插值到叶片外形中，实现气动力耦合传递机制，变化的翼型产生变化的气动力，从而又影响叶片外形；步骤3.4：开启Matlab并行计算环境：步骤3.5：在优化过程中，采用Farfor循环语句与向量化代码实现并行加速，用了parfor之后，输出参数用nargout确定；步骤4：判断翼型设计模块中目标函数与叶片设计模块中目标函数是否收敛；若否，则回转执行所述步骤1；若是，则输出新翼型族及叶片气动外形数据。