[发明专利]基于奇点分布法的轴流叶轮翼型设计方法

权利要求书

1.基于奇点分布法的轴流叶轮翼型设计方法，其特征在于：它包括以下步骤：

S1：计算总体及各流面的基本几何和流动参数；

S2：计算各流面基本翼型骨线节点上各分速度及骨线节点位置；

S3：循环完成所有流面上的骨线运算，程序结束。

2.根据权利要求1所述的基于奇点分布法的轴流叶轮翼型设计方法，其特征在于：所述的步骤S1包括以下子步骤：

S101：确定计算流面个数；

S102：循环计算各流面半径，计算各流面叶栅稠密度l/t、栅距t及弦长l；

S103：根据用户给定的泵在设计点的性能参数，循环计算各计算流面上均匀流大小W∞和流动方向角β∞，进一步计算各流面上均匀流的水平及铅垂分量W∞u、W∞z，

S104：确定各流面的冲角α后，循环计算各流面上翼型安放角β_e；

S105：循环计算各流面上绕单个翼型的环量Γ_B；

S106：以插值方法获取反映孤立翼型和栅中翼型性能差异的系数l_n的插值方程；

S107：根据各流面上已知的翼型安放角β_e和叶栅稠密度l/t确定流面的A₀值，并由此及Γ_B值循环计算各流面的系数A₁。

3.根据权利要求1所述的基于奇点分布法的轴流叶轮翼型设计方法，其特征在于：所述的步骤S2包括一个基本翼型骨线节点上各分速度的计算步骤和一个基本翼型骨线节点位置计算步骤：

基本翼型骨线节点上各分速度的计算步骤包括：

S211：将基本翼型骨线节点各点的速度分解为三项：

（1）由基本翼型骨线上的连续分布涡产生的诱导速度V_1u,V_1z；

（2）由除开基本翼型的其余翼型的连续分布涡产生的诱导速度V_2u,V_2z；

（3）均匀流的速度分量W∞u，W∞z；

S212：将一流面的长l翼型骨线等分成6段，在各节点上，由于基本翼型上的连续分布涡产生的诱导速度的两个分量V_1u,V_1z，由下式计算，V1u=-Vysinβe=-(A02+A1s0/l)sinβeV1z=Vycosβe=(A02+A1s0/l)cosβe;]]>

S213：计算各分点S₀上的a(s_o,-3l/6),a(s_o,-2l/6),......,a(s_o,3l/6)值及b值，计算流面翼型骨线各分点上的V_2u,V_2z：

将一流面上的7个分点从下向上编号，编号分别为0，1，……6，在程序中形成一个7×7的方矩数组a(7,7)，数组元素a(m,k)(m,k＝0,1,2,......6)表示编号为m的点涡在编号为k的计算点产生的诱导速度的水平分量；

该方阵关于主对角线反对称，且平行于主对线的直线上的元素相等；

正弦函数和双曲正弦函数都是奇函数，余弦函数和双曲余弦函数都是偶函数，由函数的奇偶性，可以看出a(i,j)＝-a(j,i)，其中i∈[0，6]，j∈[0，6]，每对关于主对角线对称的元素都有这一特性；

首先计算主对角线右上方的全部元素，以u1,z1表示流面骨线上任意两个相邻计算点的横坐标与纵坐标之差，在同一流面上这是两个常数；在变量p一定时（p=0，1，2，……6），q从p+1到6循环，计算p，q代表的两点的横、纵坐标之差2π(u₀-u)/t,2π(z₀-z)/t，进一步计算sh[2π(z₀-z)/t],ch[2π(z₀-z)/t],cos[2π(u₀-u)/t]，将这些值代入式a(s0,s)=-12·sh[2πt(z0-z)]ch(2πt(z0-z))-cos(2πt(u0-u))+t2πz0-z(z0-z)2+(u0-u)2]]>中，即可获得方阵右上角全部a(s₀.s)值；由方阵元素的负对称性，形成主对角线元素左下方全部元素；主对线上n+1个元素全部置0；

在利用式：

V2U=12560·πlt334a(x0,-3l6)A0+630a(x0,-2l6)A0+210a(x0,-2l6)A1-180a(x0,-l6)A0-120a(x0,-l6)A1+460a(x0,0)A0+460a(x0,0)A1-90a(x0,l6)A0-120a(x0,l6)A0+126a(x0,2l6)A0+210a(x0,2l6)A1]]>计算一个流面骨线上的V_2u时，调用方阵元素即可得到a(s_o,-3l/6),a(s_o,-2l/6),......,a(s_o,3l/6)值，以同样方法形成b方程的全部元素；

基本翼型骨线节点位置计算步骤为：在确定了骨线上各分点的全部速度分量后，以公式确定各分点上的β角，并得到最终骨线上各分点位置的坐标。

4.根据权利要求1所述的基于奇点分布法的轴流叶轮翼型设计方法，其特征在于：所述的翼型骨线上的旋涡密度函数为：或r(s)=A01-s/(l/2)1+s/(l/2)+A11-[s/(l/2)]2;]]>附加项为：f(θ)=A0tgθ2,]]>或f(s)=A01-s/(l/2)1+s/(l/2),]]>其中，l表示计算翼型的弦长，s表示弦上计算节点的曲线坐标，A₀与A₁是密度函数r(s)的两个系数，变换量θ由cosθ＝s/(l/2)定义。