一种斜放机翼型离心风机的参数优化方法与流程

1.本发明涉及离心风机设计的技术领域，尤其涉及一种斜放机翼型离心风机的参数优化方法。

背景技术：

2.对于后倾离心风机的设计，一般先设计和轮盖型线，然后通过控制到轮盖均布的几条叶片流线来控制气流的流动分布。
3.对于离心风机，目前最常规的设计方法为二维直叶片（叶片由轮盖到叶片无变化），这种叶片结构最为简单，生成工艺容易实现，但只能控制空间中一根流线的气流流动，无法控制空间上其他流线的气流流动，导致风机效率低噪音大，已经被淘汰。
4.另外一种较为复杂的设计方法为空间扭曲的三维机翼型叶片，原理上可以在空间上控制所有流线的气流流动，设计出来的风机效率和噪音均有不错的表现，但这种叶片结构最为复杂，对生产工艺及成本要求很高，实现代价太大，所以在市场应用中并不多见。此外，目前对于二维直叶片及三元扭曲叶轮均有相关的设计软件和方法，但对于机翼不等厚斜放的离心叶型设计方法及相关设计软件不太完善，工程师在设计时一般在三维软件中手动调整，比较随意且不方便。且传统设计中，工程师人工调整叶轮通过数值cfd方法或者相关经验去制作样板测试调整，反复循环调整，效率低且资源较为浪费。
5.综上所述，现有技术中对离心风机叶型的设计一般采用二维等厚叶型，结构简单但性能较差；而三维空间扭曲叶型虽然性能优但结构复杂，导致工艺实现性较差。斜放机翼型叶片可以综合二维叶型和三维扭曲叶型的优点，在工艺较容易实现的前提下性能也可达到比较优秀的状态；但是，对于斜放机翼型叶片暂无相关设计软件可以实现，传统设计方法为工程师在三维软件中进行手动调整，设计较为繁琐，且设计参数不统一；且随着人工智能技术的推进，传统的斜放机翼型叶片设计方法无法满足计算机自动优化迭代的需求。

技术实现要素：

6.本发明的目的在于提供一种斜放机翼型离心风机的参数优化方法，实现斜放机翼型叶轮叶片的参数化自动建模，并进行计算机的自动迭代优化设计。
7.一方面，本发明实施例提供了一种斜放机翼型离心风机的参数优化方法，包括如下步骤：步骤1：结合心风机叶轮构建坐标系，获取叶轮出口半径尺寸、进口与竖直线夹角、出口与水平线夹角、叶轮高度尺寸b，并基于二次贝塞尔曲线求解离心风机叶轮的曲线坐标；步骤2：获取叶轮的轮盖进口半径尺寸r0、轮盖进口与竖直线夹角θ1、轮盖出口与水平线夹角θ2、叶轮出口宽度尺寸b，并基于二次贝塞尔曲线求解离心风机叶轮的轮盖曲线坐标；
步骤3，获取叶片进口半径尺寸r1、叶片出口半径尺寸r2、子午面叶片进出口边与竖直线夹角α，根据以上参数确定离心风机叶轮子午面叶片的进出口位置，并结合曲线求解叶片的进出口位置；步骤4，结合曲线、轮盖曲线、子午面叶片进口曲线、子午面叶片进口曲线得到叶轮子午面坐标，根据子午面坐标获取离心叶轮的子午面相关参数及待优化参数；步骤5，根据叶轮子午面叶片的进出口位置以及给定的叶片角分布求解空间叶片的坐标；步骤6，根据流线s、叶片倾斜角γ和轮盖子午流线求解轮盖流线s’，并获取离心叶轮的回转面相关参数及待寻优参数；步骤7，根据上述数据对离心叶轮进行几何造型优化；步骤8，基于遗传多目标优化算法建立优化模型；步骤9，将离心风机的相关参数输入优化模型中，对优化模型进行cfd数值预测，优化模型输出待优化参数的最优解。
8.进一步的，所述步骤1还包括如下步骤：步骤1-1，结合坐标系，基于叶轮出口半径尺寸、进口与竖直线夹角、出口与水平线夹角、叶轮高度尺寸b获取控制点坐标；步骤1-2，根据控制点坐标求解曲线任一位置处的坐标值，得到曲线坐标。
9.进一步的，所述步骤1-1中，控制点包括第一控制点p1、第二控制点p2和第三控制点p3；第一控制点p1的坐标为（0,0），第三控制点p3的坐标为（，b），第二控制点p2的坐标（，）满足：）满足：其中，、分别为第一控制点p1的横、纵坐标，、分别为第三控制点p1的横纵坐标。
10.进一步的，所述步骤1-2中，曲线的坐标满足以下公式：xh（t）=t2*p
1x
+ t*（1-t）*p
2x
+ (1-t)2*p
3x
yh（t）=t2*p
1y
+ t*（1-t）*p
2y
+ (1-t)2*p
3y
其中，t为第一控制点到第三控制点之间任意一点的相对位置；xh（t）为曲线中相对位置t处的x坐标值，即曲线的横坐标，yh（t）为曲线中相对位置t处的y坐标值，即曲线的纵坐标。
11.进一步的，所述步骤2还包括如下步骤：步骤2-1，与曲线相同坐标系，并基于轮盖进口半径尺寸r0、轮盖进口与竖直线夹角θ1、轮盖出口与水平线夹角θ2、叶轮出口宽度尺寸b获取三个控制点坐标；
步骤2-2，根据三个控制点坐标求解轮盖曲线任一位置处的坐标值，从而得到轮盖曲线坐标。
12.进一步的，所述步骤2-1中，控制点包括第四控制点p4、第五控制点p5和第六控制点p6；第四控制点p4的坐标为（,0）;第六控制点p6坐标为（r3，b-b）；第五控制点p5的坐标满足：第五控制点p5的坐标满足：其中,p
4x
、p
4y
分别为第四控制点p4的横、纵坐标，p
6x
、p
6y
分别为第四控制点p6的横、纵坐标，、分别为第五控制点p5的横、纵坐标。
13.进一步的，所述步骤2-2中，轮盖曲线的坐标满足：xs（t）=t2*p
3x
+ t*（1-t）*p
4x
+ (1-t)2*p
5xys
（t）=t2*p
3y
+ t*（1-t）*p
4y
+ (1-t)2*p
5y
其中，t∈[0,1]，0位置为起始位置p3,1位置为终止位置p5，t为0到1之间任一点的相对位置，xs（t）为轮盖曲线任一相对位置t处的x坐标值，即曲线的横坐标，ys（t）为轮盖曲线任一相对位置t处的y坐标值，即曲线的纵坐标。
[0014]
进一步的，所述步骤3还包括：步骤3-1，根据叶片进口半径尺寸r1、叶片出口半径尺寸r2、子午面叶片进出口边与竖直线夹角α，求解子午面叶片进口点leh点坐标（leh
x
, lehy）；子午面叶片进口点leh点坐标满足：leh
x
=r1leh
x
=t2*p
1x
+ t*（1-t）*p
2x
+ (1-t)2*p
3x
lehy=t2*p
1y
+ t*（1-t）*p
2y
+ (1-t)2*p
3y
其中，leh
x
和lehy分别为子午面叶片进口点leh的x坐标值和y坐标值，t∈[0,1]；步骤3-2，根据叶片进口半径尺寸r1、叶片出口半径尺寸r2、子午面叶片进出口边与竖直线夹角α，求解子午面叶片出口点teh点坐标（teh
x
, tehy）；子午面叶片出口点teh点坐标满足：teh
x
=r2teh
x =t2*p
1x
+ t*（1-t）*p
2x
+ (1-t)2*p
3x
teh
y =t2*p
1y
+ t*（1-t）*p
2y
+ (1-t)2*p
3y
其中，teh
x
和tehy分别为分别为子午面叶片出口点teh点的x坐标值；步骤3-3，通过叶片进口半径尺寸r1、叶片出口半径尺寸r2、子午面叶片进（出）口边与竖直线夹角α和上述轮盖曲线公式联立求解，得到轮盖叶片进口点的坐标值（les
x
, lesy）；子午面叶片进口直线满足
y=tanα*x+lehy‑ꢀ
tanα* leh
x
根据子午面轮盖曲线的横坐标xs（t）和纵坐标ys（t）满足：xs（t）=t2*p
3x
+ t*（1-t）*p
4x
+ (1-t)2*p
5xys
（t）=t2*p
3y
+ t*（1-t）*p
4y
+ (1-t)2*p
5y
联立以上方程组求解子午面叶片进口直线与轮盖曲线交点坐标，此点即为轮盖叶片进口点les；步骤3-4，通过叶片进口半径尺寸r1、叶片出口半径尺寸r2、子午面叶片进（出）口边与竖直线夹角α和上述轮盖曲线公式联立求解，得到轮盖叶片出口点的坐标值（tes
x
, tesy）；结合步骤3-2，子午面叶片出口直线满足：y=tanα*x+tehy‑ꢀ
tanα* teh
x
然后，根据子午面轮盖曲线公式：xs（t）=t2*p
3x
+ t*（1-t）*p
4x
+ (1-t)2*p
5xys
（t）=t2*p
3y
+ t*（1-t）*p
4y
+ (1-t)2*p
5y
联立以上方程组求解子午面叶片出口直线与轮盖曲线交点坐标，此点即为轮盖叶片出口点tes。
[0015]
进一步的，所述步骤4中，离心叶轮子午面的相关参数包括：出口半径尺寸r3、叶轮总宽尺寸b、进口与竖直线夹角θ3、出口与水平线夹角θ4、轮盖进口与竖直线夹角θ、叶轮进口半径尺寸r0、叶片进口半径尺寸r1、叶片出口半径尺寸r2、轮盖出口与水平线夹角θ2叶轮出口宽度尺寸b、子午面叶片进出口边与竖直线夹角α；离心叶轮子午面的待优化参数包括：叶轮进口半径尺寸r0、叶片进口半径尺寸r1、叶片出口半径尺寸r2、轮盖出口与水平线夹角θ2叶轮出口宽度尺寸b、子午面叶片进出口边与竖直线夹角α。
[0016]
进一步的，所述步骤4中，离心叶轮的回转面相关参数及待寻优参数包括：叶片进口轴向角β1、叶片出口轴向角β2、流线s中部位置的叶片轴向角β相对值λ以及叶片倾斜角γ。
[0017]
本发明的有益效果为：本发明提出的斜放机翼型离心风机参数优化方法，针对机翼不等厚的直叶型，在空间上叶片采用斜放处理，综合二维直叶片和三维扭曲叶片的优点，从而能在一定程度上调整空间上的气流流动，且工艺较为简单，生产成本较低；本发明提出了针对斜放机翼型叶轮叶片的参数化设计建模方法，完善了此类型叶片的设计方法，并对参数进行自动优化，实现参数化自动建模、计算及优化；本发明能够实现叶轮的参数化自动建模，并进行计算机的自动迭代优化设计。
附图说明
[0018]
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在附图中：图1为本发明实施例中提供的斜放机翼型离心风机的参数优化方法的整体流程示
意图；图2为本发明实施例中斜放机翼型叶轮叶片的子午面示意图；图3为本发明实施例中控制点的坐标示意图；图4为本发明实施例中轮盖控制点的坐标示意图；图5为本发明实施例中叶片包角与空间叶片的结构示意图；图6为本发明实施例中离心叶轮子午面和回转面的结构关系示意图；图7为本发明实施例中流线和叶型之间的空间关系示意图；图8为本发明实施例中轮盖空间叶型参数及结构示意图。
具体实施方式
[0019]
为使本发明实施例的目的、技术方案和优点更加清楚明白，下面结合附图对本发明实施例做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。
[0020]
在本说明书的描述中，所使用的“包含”、“包括”、“具有”、“含有”等，均为开放性的用语，即意指包含但不限于。参考术语“一个实施例”、“一个具体实施例”、“一些实施例”、“例如”等的描述意指结合该实施例或示例描述的具体特征、结构或者特点包含于本技术的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。各实施例中涉及的步骤顺序用于示意性说明本技术的实施，其中的步骤顺序不作限定，可根据需要作适当调整。
[0021]
本发明实施例中提供了一种斜放机翼型离心风机参数化优化方法，图1为该方法的整体流程示意图，如图1所示，该方法包括如下步骤：步骤1，获取叶轮出口半径尺寸r3、进口与竖直线夹角θ3、出口与水平线夹角θ4、叶轮高度尺寸b，并基于二次贝塞尔曲线求解离心风机叶轮的曲线坐标；具体的，步骤1的求解曲线坐标还包括如下步骤：步骤1-1，在叶轮上构建坐标系，并基于叶轮出口半径尺寸r3、进口与竖直线夹角θ3、出口与水平线夹角θ4、叶轮高度尺寸b获取三个控制点坐标；如图3所示，第一控制点p1的坐标为（0,0），即第一控制点p1为所构建坐标系的坐标原点；第二控制点p2的坐标(p
2x
,p
2y
)由第一控制点p1和第三控制点p3的坐标，以及进口与竖直线夹角θ3、出口与水平线夹角θ4求解，第二控制点p2的坐标计算式为：进口与竖直线夹角θ3、出口与水平线夹角θ4求解，第二控制点p2的坐标计算式为：其中,p
1x
为p1点的x坐标值，p
1y
为p1点的y坐标值，p
2x
为p2点的x坐标值，p
2y
为p2点的y坐标值, p
3x
为p3点的x坐标值，p
3y
为p3点的y坐标值。
[0022]
第三控制点p3的坐标由叶轮出口半径尺寸r3、叶轮高度尺寸b两个尺寸确定，其坐标为（r3，b）；步骤1-2，根据三个控制点坐标求解曲线任一位置处的坐标值，得到曲线
坐标；具体的，曲线的横纵坐标满足以下公式：xh（t）=t2*p
1x
+ t*（1-t）*p
2x
+ (1-t)2*p
3x
yh（t）=t2*p
1y
+ t*（1-t）*p
2y
+ (1-t)2*p
3y
其中，设第一控制点p1为起始位置0，第三控制点p3为终止位置1，t∈[0,1]，即t为第一控制点p1到第三控制点p3之间任一点的相对位置；xh（t）为曲线中任一相对位置t处的x坐标值，即曲线的横坐标，yh（t）为曲线中任一相对位置t处的y坐标值，即曲线的纵坐标；p
1x
为p1点的x坐标值，p
1y
为p1点的y坐标值，p
2x
为p2点的x坐标值，p
2y
为p2点的y坐标值, p
3x
为p3点的x坐标值，p
3y
为p3点的y坐标值。
[0023]
步骤2，获取叶轮的轮盖进口半径尺寸r0、轮盖进口与竖直线夹角θ1、轮盖出口与水平线夹角θ2、叶轮出口宽度尺寸b，并基于二次贝塞尔曲线求解轮盖曲线坐标；具体的，步骤2的求解轮盖曲线坐标还包括如下步骤：步骤2-1，与曲线相同坐标系，并基于轮盖进口半径尺寸r0、轮盖进口与竖直线夹角θ1、轮盖出口与水平线夹角θ2、叶轮出口宽度尺寸b获取三个控制点坐标；如图4所示，第四控制点p4的坐标为（r0,0）;第六控制点p6的坐标由叶轮出口半径尺寸r3、叶轮出口宽度尺寸b两个尺寸确定，其坐标为（r3，b-b）；第二控制点p5的坐标由第一控制点p4和第三控制点p6的坐标，以及轮盖进口与竖直线夹角θ1、轮盖出口与水平线夹角θ2求解；直线夹角θ1、轮盖出口与水平线夹角θ2求解；其中,p
4x
为p4点的x坐标值，p
4y
为p4点的y坐标值，p
5x
为p5点的x坐标值，p
5y
为p5点的y坐标值, p
6x
为p6点的x坐标值，p
6y
为p6点的y坐标值。
[0024]
步骤2-2，根据三个控制点坐标求解轮盖曲线任一位置处的坐标值，从而得到轮盖曲线坐标；具体的，轮盖曲线的横纵坐标满足以下公式：xs（t）=t2*p
3x
+ t*（1-t）*p
4x
+ (1-t)2*p
5xys
（t）=t2*p
3y
+ t*（1-t）*p
4y
+ (1-t)2*p
5y
其中，设第三控制点p3为起始位置0，第五控制点p5为终止位置1，t∈[0,1]，即t为第三控制点p3到第五控制点p5之间任一点的相对位置。
[0025]
p
3x
为p3点的x坐标值，p
3y
为p3点的y坐标值，p
4x
为p4点的x坐标值，p
4y
为p4点的y坐标值, p
5x
为p5点的x坐标值，p
5y
为p5点的y坐标值,xs（t）为轮盖曲线任一相对位置t处的x坐标值，ys（t）为轮盖曲线任一相对位置t处的y坐标值。
[0026]
步骤3，获取叶片进口半径尺寸、叶片出口半径尺寸、子午面叶片进出口边与竖直线夹角α，根据以上参数确定叶轮子午面叶片的进出口位置，并结合曲线求解叶片的进出口位置；
其中，步骤3还包括以下步骤：步骤3-1，根据叶片进口半径尺寸r1、叶片出口半径尺寸r2、子午面叶片进出口边与竖直线夹角α，求解子午面叶片进口点leh点坐标（leh
x
, lehy）；子午面叶片进口点leh点坐标满足：leh
x
=r1leh
x
=t2*p
1x
+ t*（1-t）*p
2x
+ (1-t)2*p
3x
lehy=t2*p
1y
+ t*（1-t）*p
2y
+ (1-t)2*p
3y
其中，leh
x
和lehy分别为子午面叶片进口点leh的x坐标值和y坐标值，t∈[0,1]，p
1x
为第一控制点p1点的x坐标值，p
1y
为第一控制点p1点的y坐标值，p
2x
为第二控制点p2点的x坐标值，p
2y
为第二控制点p2点的y坐标值, p
3x
为第三控制点p3点的x坐标值，p
3y
为第三控制点p3点的y坐标值；上述三元一次方程通过盛金公式可解，具体先联立前两式根据盛金公式求解三元一次方程的根t的值，在根据求解出的t值代入最后一式求解出lehy。
[0027]
步骤3-2，根据叶片进口半径尺寸r1、叶片出口半径尺寸r2、子午面叶片进出口边与竖直线夹角α，求解子午面叶片出口点teh点坐标（teh
x
, tehy）；子午面叶片出口点teh点坐标满足：teh
x
=r2teh
x =t2*p
1x
+ t*（1-t）*p
2x
+ (1-t)2*p
3x
teh
y =t2*p
1y
+ t*（1-t）*p
2y
+ (1-t)2*p
3y
其中，p
1x
为第一控制点p1点的x坐标值，p
1y
为第一控制点p1点的y坐标值，p
2x
为第二控制点p2点的x坐标值，p
2y
为第二控制点p2点的y坐标值, p
3x
为第三控制点p3点的x坐标值，p
3y
为第三控制点p3点的y坐标值，teh
x
为teh点的x坐标值，tehy为teh点的y坐标值。
[0028]
可以理解的是，子午面叶片出口点坐标的求解方式与步骤3-1中子午面叶片进口点坐标的求解方式相似，通过盛金公式可以完成求解。
[0029]
步骤3-3，通过叶片进口半径尺寸r1、叶片出口半径尺寸r2、子午面叶片进出口边与竖直线夹角α和上述轮盖曲线公式联立求解，得到轮盖叶片进口点的坐标值（les
x
, lesy）；具体的，通过叶片进口半径尺寸r1、叶片出口半径尺寸r2、子午面叶片进出口边与竖直线夹角α和上述轮盖曲线公式联立求解，采用计算机迭代求解出轮盖叶片进口点及轮盖叶片出口点坐标值。
[0030]
如图2所示，子午面叶片进（出）口边与竖直线夹角α、叶片进口点及叶片出口点坐标值并结合上述轮盖曲线通过公式联立求解，采用计算机迭代求解出通过叶片出口点坐标、倾斜角为α的直线与轮盖曲线交点坐标，以及通过叶片出口点坐标、倾斜角为α的直线与轮盖曲线交点坐标，此两个交点即为轮盖上叶片的起始和终止点，起始点和终止点分别为子午面叶片的进口位置和出口位置。
[0031]
具体的，结合步骤3-1，子午面叶片进口直线满足y=tanα*x+lehy‑ꢀ
tanα* leh
x
根据子午面轮盖曲线的横坐标xs（t）和纵坐标ys（t）满足：xs（t）=t2*p
3x
+ t*（1-t）*p
4x
+ (1-t)2*p
5xys
（t）=t2*p
3y
+ t*（1-t）*p
4y
+ (1-t)2*p
5y
联立以上方程组求解子午面叶片进口直线与轮盖曲线交点坐标，此点即为轮盖叶片进口点les。
[0032]
步骤3-4，通过叶片进口半径尺寸r1、叶片出口半径尺寸r2、子午面叶片进出口边与竖直线夹角α和上述轮盖曲线公式联立求解，得到轮盖叶片出口点的坐标值（tes
x
, tesy）；结合步骤3-2，子午面叶片出口直线满足：y=tanα*x+tehy‑ꢀ
tanα* teh
x
然后，根据子午面轮盖曲线公式：xs（t）=t2*p
3x
+ t*（1-t）*p
4x
+ (1-t)2*p
5xys
（t）=t2*p
3y
+ t*（1-t）*p
4y
+ (1-t)2*p
5y
联立以上方程组求解子午面叶片出口直线与轮盖曲线交点坐标，此点即为轮盖叶片出口点tes。
[0033]
步骤4，结合曲线、轮盖曲线、子午面叶片进口曲线、子午面叶片进口曲线得到叶轮子午面坐标，根据子午面坐标获取离心叶轮的子午面相关参数及待优化参数；具体的，叶轮子午面坐标包括步骤1~3得到的曲线、轮盖曲线、子午面叶片进口曲线、子午面叶片进口线，因此叶轮子午面坐标通过上述11个参数可解；在优化设计中，由于出口半径尺寸r3、叶轮总宽尺寸b一般为固定尺寸，且进口与竖直线夹角θ3、出口与水平线夹角θ4、轮盖进口与竖直线夹角θ1、在一定范围内对叶轮性能影响较小，因此在本发明的优化设计中以上五个参数为定值，在优化过程中不对其进行调整，具体大小由本领域技术人员根据实际情况预先进行设置。
[0034]
因此本发明在子午面上通过对叶轮性能影响较大的叶轮进口半径尺寸r0、叶片进口半径尺寸r1、叶片出口半径尺寸r2、轮盖出口与水平线夹角θ2叶轮出口宽度尺寸b、子午面叶片进出口边与竖直线夹角α共6个参数的取值范围进行寻优。
[0035]
步骤5，根据叶轮子午面叶片的进出口位置以及给定的叶片角分布求解空间叶片的坐标；步骤5-1，求解叶片包角θ；具体的，如图4所示，给定流线s上起始的两个位置处（0/1）的β值叶片进口轴向角β1、叶片出口轴向角β2，并给定流线s中部位置（0.5）的叶片轴向角β相对值λ（λ取值0~1），即可得到s中部位置β值为：β=β1+λ*（β2-β1）运用插值函数拟合得到流线s上所有点的叶片轴向角β值，再根据上图的β和叶型之间的空间关系图即可得到流线s上的叶型对应点的叶片包角θ值，叶片包角θ公式如下：其中，r——叶片上某一点半径位置尺寸d
θ
表示叶片包角的微分值；dm表示子午流线长度的微分值。
[0036]
进一步的，子午面叶片进出口所夹的轮盖两条流线即为叶片的两条子午流
线，根据子午面叶片进出口坐标、曲线和轮盖曲线可以求解出子午流线对应式，通过该式积分即可求解出子午流线对应r处的叶片包角θ。
[0037]
步骤5-2，根据叶片包角θ获取空间叶片的坐标值；如图5所示，根据叶片包角θ即可得到空间叶片任一点流线s的坐标值（sx,sy,sz）：其中：sx=r*sin(θ)，sy=r*cos(θ)，sz与子午流线上r位置处的z值相同。
[0038]
步骤6，根据流线s、叶片倾斜角γ和轮盖子午流线求解轮盖流线s’，并获取离心叶轮的回转面相关参数及待寻优参数；具体的，如图5所示，由于叶片为直叶片斜摆方式，通过斜摆方式调整空间气流流动，流线和轮盖流线在空间上是扫略关系，参数设计时规定叶片倾斜角γ，通过步骤5求解出流线s的坐标后，在流线s坐标已知情况下可根据γ和轮盖子午流线求解出轮盖流线s’。
[0039]
优选的，在叶片中线上施加一定的不等厚度分布，可以改善叶片上的气流分布，但机翼厚度分布在一定范围内对离心风机影响较小，在本发明中设计过程中，机翼厚度为定值，由本领域技术人员根据实际情况预先设置。
[0040]
进一步的，如图6所示，图6为离心叶轮子午面和回转面的结构关系示意图，子午面包括轮盖两条曲线以及叶片进出口两条直线，回转面为到轮盖中间的任一曲面，回转面设计即为在空间的回转面上对叶片的曲线进行设计。
[0041]
因此本发明在回转面上通过控制流线s上起始的两个位置处（0/1）的β值，针对叶片进口轴向角β1、叶片出口轴向角β2、流线s中部位置（0.5）的叶片轴向角β相对值λ以及叶片倾斜角γ共4个参数的取值范围进行寻优，其中λ取值为0~1。
[0042]
步骤7，根据上述数据对离心叶轮进行几何造型优化；根据上述参数求解出的型线坐标、轮盖型线坐标、流线坐标、轮盖流线坐标对离心叶轮进行3d几何造型优化，3d几何造型可以采用成熟的3d软件进行二次开发，如solidworks、ug、proe等软件，本发明采用solidworks软件。
[0043]
将型线坐标、轮盖型线坐标、流线坐标、轮盖流线坐标输入3d软件中，得到3d造型。
[0044]
通过上述坐标可实现自动3d导入造型的构建。
[0045]
步骤8，基于遗传多目标优化算法建立优化模型；其中，本发明所采用优化算法为遗传多目标优化算法；由于风机主要关注压升、流量以及效率三个参数，而在cfd数值计算中，一般边界条件会将流量或者压升给定，返求另一个物理量及效率；具体的，若给定流量q进行数值计算，则以风机压升p和效率η作为优化目标，优化决策函数y为：y=max{p，η}若给定压升p进行数值计算，则以风机流量q和效率η作为优化目标，优化决策函数为：y=max{q，η}进一步的，设定最大迭代优化次数maxstep。
[0046]
通过优化决策函数和cfd数值预测，便可通过优化算法进行打靶寻优，得到后续使用的优化模型。
[0047]
步骤9，将离心风机的相关参数输入优化模型中，对优化模型进行cfd数值预测，优化模型输出待优化参数的最优解。
[0048]
具体的，优化模型的输入为叶轮的5个定值参数以及10个待优化参数，结合步骤4和步骤六，子午面相关参数共11个参数，回转面相关参数共4个参数，子午面相关参数中5个为定值参数，即本发明预先设置好且不进行优化的参数，剩余6个参数和回转面相关4个参数作为待优化参数通过优化模型进行优化处理；优化模型的输出结果为10个待优化参数的优选值集合，包括子午面6个参数，回转面4个参数。
[0049]
通过对cfd仿真工具进行二次开发，实现自动区域离散、网格划分、计算及结果分析，预测出该模型的压升、流量及效率值，通过这三个值评估优化决策函数y值，再通过优化算法综合决策并进入下一代优化，经过多次迭代后后即可得到最优的解集。
[0050]
优选的，本发明给定流线参数返算轮盖流线，也可给定轮盖流线参数返算流线，也可给定任一流线返求轮盖流线；优选的，作为一种可选的方式，本发明给定的各参数也可以用无量纲参数替代，如r0，r1、r2
ꢀ¬
的取值可以以r3为基准，用r0/r3、r1/r3、r2/r3三个无量纲参数替代。
[0051]
优选的，作为一种可选的方式，本发明中轮盖流线不限于二次贝塞尔曲线给定，也可以为三次或者更高阶次，也可以是其他曲线或者曲线加直线的组合方式。
[0052]
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质（包括但不限于磁盘存储器、cd-rom、光学存储器等）上实施的计算机程序产品的形式。
[0053]
本发明是参照根据本发明实施例的方法、设备（系统）、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
[0054]
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
[0055]
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
[0056]
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详
细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

技术特征：

1.一种斜放机翼型离心风机的参数优化方法，其特征在于，包括如下步骤：步骤1：构建离心风机叶轮构建坐标系，获取叶轮出口半径尺寸、进口与竖直线夹角、出口与水平线夹角和叶轮高度尺寸b，结合以上参数并基于二次贝塞尔曲线求解离心风机叶轮的曲线坐标；步骤2：获取叶轮的轮盖进口半径尺寸、轮盖进口与竖直线夹角、轮盖出口与水平线夹角和叶轮出口宽度尺寸b，结合以上参数并基于二次贝塞尔曲线求解离心风机叶轮的轮盖曲线坐标；步骤3，获取叶片进口半径尺寸、叶片出口半径尺寸、子午面叶片进出口边与竖直线夹角α，根据以上参数确定离心风机叶轮子午面叶片的进出口位置，并结合曲线求解叶片的进出口位置；步骤4，结合曲线、轮盖曲线、子午面叶片进口曲线、子午面叶片进口曲线得到叶轮子午面坐标，根据子午面坐标获取离心叶轮的子午面相关参数及待优化参数；步骤5，根据叶轮子午面叶片的进出口位置以及给定的叶片角分布求解空间叶片的坐标；步骤6，根据流线s、叶片倾斜角γ和轮盖子午流线求解轮盖流线s’，并获取离心叶轮的回转面相关参数及待寻优参数；步骤7，根据上述数据对离心叶轮进行几何造型优化；步骤8，基于遗传多目标优化算法建立优化模型；步骤9，将离心风机的相关参数输入优化模型中，对优化模型进行cfd数值预测，优化模型输出待优化参数的最优解。2.如权利要求1所述的斜放机翼型离心风机的参数优化方法，其特征在于，所述步骤1还包括如下步骤：步骤1-1，结合坐标系，基于叶轮出口半径尺寸、进口与竖直线夹角、出口与水平线夹角、叶轮高度尺寸b获取控制点坐标；步骤1-2，根据控制点坐标求解曲线任一位置处的坐标值，得到曲线坐标。3.如权利要求2所述的斜放机翼型离心风机的参数优化方法，其特征在于，所述步骤1-1中，控制点包括第一控制点p1、第二控制点p2和第三控制点p3；第一控制点p1的坐标为（0,0），第三控制点p3的坐标为（，b），第二控制点p2的坐标（，）满足：）满足：其中，、分别为第一控制点p1的横、纵坐标，、分别为第三控制点p1的横纵
坐标。4.如权利要求3所述的斜放机翼型离心风机的参数优化方法，其特征在于，所述步骤1-2中，曲线的坐标满足以下公式：x
h
（t）=t2*p
1x
+ t*（1-t）*p
2x
+ (1-t)2*p
3x
y
h
（t）=t2*p
1y
+ t*（1-t）*p
2y
+ (1-t)2*p
3y
其中，t为第一控制点到第三控制点之间任意一点的相对位置，x
h
（t）为曲线中相对位置t处的x坐标值，即曲线的横坐标，y
h
（t）为曲线中相对位置t处的y坐标值，即曲线的纵坐标。5.如权利要求4所述的斜放机翼型离心风机的参数优化方法，其特征在于，所述步骤2还包括如下步骤：步骤2-1，与曲线相同坐标系，并基于轮盖进口半径尺寸r0、轮盖进口与竖直线夹角θ1、轮盖出口与水平线夹角θ2、叶轮出口宽度尺寸b获取三个控制点坐标；步骤2-2，根据三个控制点坐标求解轮盖曲线任一位置处的坐标值，从而得到轮盖曲线坐标。6.如权利要求5所述的斜放机翼型离心风机的参数优化方法，其特征在于，所述步骤2-1中，控制点包括第四控制点p4、第五控制点p5和第六控制点p6；第四控制点p4的坐标为（,0）;第六控制点p6坐标为（，b-b）；第五控制点p5的坐标满足：第五控制点p5的坐标满足：其中,p
4x
、p
4y
分别为第四控制点p4的横、纵坐标，p
6x
、p
6y
分别为第四控制点p6的横、纵坐标，、分别为第五控制点p5的横、纵坐标。7.如权利要求6所述的斜放机翼型离心风机的参数优化方法，其特征在于，所述步骤2-2中，轮盖曲线的坐标满足：x
s
（t）=t2*p
3x
+ t*（1-t）*p
4x
+ (1-t)2*p
5x
y
s
（t）=t2*p
3y
+ t*（1-t）*p
4y
+ (1-t)2*p
5y
其中，t∈[0,1]，0位置为起始位置p3,1位置为终止位置p5，t为0到1之间任一点的相对位置，x
s
（t）为轮盖曲线任一相对位置t处的x坐标值，即曲线的横坐标，y
s
（t）为轮盖曲线任一相对位置t处的y坐标值，即曲线的纵坐标。8.如权利要求7所述的斜放机翼型离心风机的参数优化方法，其特征在于，所述步骤3还包括：步骤3-1，根据叶片进口半径尺寸r1、叶片出口半径尺寸r2、子午面叶片进出口边与竖直线夹角α，求解子午面叶片进口点leh点坐标（leh
x
, leh
y
）；子午面叶片进口点leh点坐标满足：leh
x
=r1leh
x
=t2*p
1x
+ t*（1-t）*p
2x
+ (1-t)2*p
3x
leh
y
=t2*p
1y
+ t*（1-t）*p
2y
+ (1-t)2*p
3y
其中，leh
x
和leh
y
分别为子午面叶片进口点leh的x坐标值和y坐标值，t∈[0,1]；步骤3-2，根据叶片进口半径尺寸、叶片出口半径尺寸、子午面叶片进出口边与竖直线夹角α，求解子午面叶片出口点teh点坐标（teh
x
, teh
y
）；子午面叶片出口点teh点坐标满足：teh
x
=teh
x =t2*p
1x
+ t*（1-t）*p
2x
+ (1-t)2*p
3x
teh
y =t2*p
1y
+ t*（1-t）*p
2y
+ (1-t)2*p
3y
其中，teh
x
和teh
y
分别为分别为子午面叶片出口点teh点的x坐标值；步骤3-3，通过叶片进口半径尺寸、叶片出口半径尺寸、子午面叶片进出口边与竖直线夹角α和上述轮盖曲线公式联立求解，得到轮盖叶片进口点的坐标值（les
x
, les
y
）；子午面叶片进口直线满足y=tanα*x+leh
y
‑ꢀ
tanα* leh
x
根据子午面轮盖曲线的横坐标x
s
（t）和纵坐标y
s
（t）满足：x
s
（t）=t2*p
3x
+ t*（1-t）*p
4x
+ (1-t)2*p
5x
y
s
（t）=t2*p
3y
+ t*（1-t）*p
4y
+ (1-t)2*p
5y
联立以上方程组求解子午面叶片进口直线与轮盖曲线交点坐标，此点即为轮盖叶片进口点les；步骤3-4，通过叶片进口半径尺寸、叶片出口半径尺寸、子午面叶片进出口边与竖直线夹角α和上述轮盖曲线公式联立求解，得到轮盖叶片出口点的坐标值（tes
x
, tes
y
）；结合步骤3-2，子午面叶片出口直线满足：y=tanα*x+teh
y
‑ꢀ
tanα* teh
x
然后，根据子午面轮盖曲线公式：x
s
（t）=t2*p
3x
+ t*（1-t）*p
4x
+ (1-t)2*p
5x
y
s
（t）=t2*p
3y
+ t*（1-t）*p
4y
+ (1-t)2*p
5y
联立以上方程组求解子午面叶片出口直线与轮盖曲线交点坐标，此点即为轮盖叶片出口点tes。9.如权利要求8所述的斜放机翼型离心风机的参数优化方法，其特征在于，所述步骤4中，离心叶轮子午面的相关参数包括：出口半径尺寸、叶轮总宽尺寸b、进口与竖直线夹角θ3、出口与水平线夹角θ4、轮盖进口与竖直线夹角θ、叶轮进口半径尺寸、叶片进口半径尺寸、叶片出口半径尺寸、轮盖出口与水平线夹角θ2叶轮出口宽度尺寸b、子午面叶片进出口边与竖直线夹角α；离心叶轮子午面的待优化参数包括：叶轮进口半径尺寸、叶片进口半径尺寸、叶片出口半径尺寸、轮盖出口与水平线夹角θ2叶轮出口宽度尺寸b、子午面叶片进出
口边与竖直线夹角α。10.如权利要求8所述的斜放机翼型离心风机的参数优化方法，其特征在于，所述步骤4中，离心叶轮的回转面相关参数及待寻优参数包括：叶片进口轴向角β1、叶片出口轴向角β2、流线s中部位置的叶片轴向角β相对值λ以及叶片倾斜角γ。

技术总结

本发明提供了一种斜放机翼型离心风机的参数优化方法，包括如下步骤：基于离心风机结构构建坐标系，获取离心风机结构的相关参数及其关系，得到叶轮子午面相关参数及待优化参数、以及离心叶轮的回转面相关参数及待寻优参数，构建优化模型，将离心风机的相关参数输入优化模型中，对优化模型进行CFD数值预测，优化模型输出待优化参数的最优解。本发明提出了针对斜放机翼型叶轮叶片的参数化设计建模方法，完善了此类型叶片的设计方法，并对参数进行自动优化，实现参数化自动建模、计算及优化；能够实现叶轮的参数化自动建模，并进行计算机的自动迭代优化设计。动迭代优化设计。动迭代优化设计。