CN102251983A - 一种抗空蚀离心泵叶轮优化设计方法 - Google Patents
一种抗空蚀离心泵叶轮优化设计方法 Download PDFInfo
- Publication number
- CN102251983A CN102251983A CN2011102025245A CN201110202524A CN102251983A CN 102251983 A CN102251983 A CN 102251983A CN 2011102025245 A CN2011102025245 A CN 2011102025245A CN 201110202524 A CN201110202524 A CN 201110202524A CN 102251983 A CN102251983 A CN 102251983A
- Authority
- CN
- China
- Prior art keywords
- dimensional airfoil
- bone line
- dimensional
- airfoil
- front shroud
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F04—POSITIVE - DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS FOR LIQUIDS OR ELASTIC FLUIDS
- F04D—NON-POSITIVE-DISPLACEMENT PUMPS
- F04D29/00—Details, component parts, or accessories
- F04D29/18—Rotors
- F04D29/22—Rotors specially for centrifugal pumps
- F04D29/2205—Conventional flow pattern
- F04D29/2216—Shape, geometry
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F04—POSITIVE - DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS FOR LIQUIDS OR ELASTIC FLUIDS
- F04D—NON-POSITIVE-DISPLACEMENT PUMPS
- F04D29/00—Details, component parts, or accessories
- F04D29/66—Combating cavitation, whirls, noise, vibration or the like; Balancing
- F04D29/669—Combating cavitation, whirls, noise, vibration or the like; Balancing especially adapted for liquid pumps
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Structures Of Non-Positive Displacement Pumps (AREA)
Abstract
本发明一种抗空蚀离心泵叶轮优化设计方法,获得原始叶轮的周向XY平面上的前盖板二维翼型骨线、中间二维翼型骨线以及后盖板二维翼型骨线,并参数化后,构建叶轮优化参数,构建变化后的骨线,再将其三维化得到的变化后的前盖板三维翼型型线、中间三维翼型型线以及后盖板三维翼型型线,输入到几何造型软件中,形成变化后的叶轮三维形状,在单流道网格划分的基础上进行水利性能测算,以多目标遗传算法NSGA-II遗传算法为优化工具,对离心泵叶轮参数进行多目标优化设计,选取水利性能最佳的一组叶轮参数作为优化后的叶轮参数;本发明方法优化参数构建科学,优化过程自动化程度高,优化后离心泵叶轮叶片效率,空化空蚀性能优良。
Description
技术领域
本发明属于离心泵叶轮制造技术领域,具体涉及一种抗空蚀离心泵叶轮优化设计方法。
背景技术
离心泵是使用最广泛的通用机械之一,在电力工业领域、化学工业领域等都具有广泛的应用。离心泵的作用是根据离心力原理,使叶轮叶片高速旋转从而带动水转动并将水甩出,达到输送水的目的。
我国现有的离心泵普遍效率较低,空蚀性能差,其主要原因是优化设计方法落后。离心泵内部的流动是复杂的全三维粘性流动,内部存在二次流以及分离流动在内的复杂非定常流动,而现在大部分离心泵叶轮的优化设计都是采用一元或二元理论,尽管有的一些厂家和研究所采用的是全三元的设计方法,但所采用的全三元的设计方法是基于无粘假设建立的,无法得出准确的符合离心泵叶轮内部流动特性的叶轮形状。
发明内容
本发明的目的是提供一种抗空蚀离心泵叶轮优化设计方法,优化参数构建科学,优化过程自动化程度高,优化后离心泵叶轮叶片效率,空化空蚀性能优良。
本发明所采用的技术方案是一种抗空蚀离心泵叶轮优化设计方法,包括以下步骤:
步骤1、得到原始叶轮的周向XY平面上的前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线Gh;
步骤2、采用四点三次贝塞尔曲线分别对上述前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线Gh进行参数化,以得到各骨线的四个控制点P0,P1,P2,P3的在周向XY平面上的坐标,其中,控制点P0为首点,控制点P3为末点,控制点P1和P2为中间点;
步骤3、构建叶轮优化参数:
步骤3.1、构建前盖板二维翼型骨线Gs的挠度参数C1和C2:
步骤3.2、依照步骤3.1得到中间二维翼型骨线Gm的挠度参数C3和C4,以及后盖板二维翼型骨线Gh的挠度参数C5和C6;
步骤3.3、定义前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线Gh各自控制点P0和P3的变化量:
将前盖板二维翼型骨线Gs的控制点P0的柱坐标为(rs,θs),中间二维翼型骨线Gm的控制点P0的柱坐标为(rm,θm),后盖板二维翼型骨线Gh的控制点P0的柱坐标为(rh,θh);定义各控制点P0的r坐标的优化变化量为Δrs,Δrm,Δrh,定义各控制点P0的θ的坐标变化量为Δθs,Δθm,Δθh;
所述前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线Gh的控制点P3点均重合,且该P3点的柱坐标为(re,θe),定义该P3点的θ的坐标变化量为Δθe;
步骤4、构建变化后的前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线Gh:
依照上述方法,变化C3和C4的取值,得到中间二维翼型骨线Gm变化后的控制点P1和P2在周向XY平面上的坐标;变化C5和C6的取值,得到后盖板二维翼型骨线Gh变化后的控制点P1和P2在周向XY平面上的坐标;
步骤4.2、变化Δrs,Δrm,Δrh的取值,以及Δθs,Δθm,Δθh的取值,分别得到前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线Gh变化后的控制点P0的柱坐标值,根据柱坐标系与平面直角坐标系的转换关系,分别得到前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线Gh变化后的控制点P0在周向XY平面上的坐标;
步骤4.3、变化Δθe的取值,分别得到前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线变化后的Gh的控制点P3点的柱坐标值,根据柱坐标系与平面直角坐标系的转换关系,得到前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线Gh变化后的控制点P3点在周向XY平面上的坐标;
步骤4.4、根据变化后的前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线Gh变化后的四个控制点P0,P1,P2,P3的坐标,构建得到在周向XY平面上的变化后的前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线Gh;
步骤5、将步骤4得到变化后的前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线Gh,按照翼型厚度分布规律进行加厚,分别得到变化后的前盖板二维翼型型线中间二维翼型型线以及后盖板二维翼型型线
步骤6、求得步骤5得到的变化后的前盖板二维翼型型线中间二维翼型型线以及后盖板二维翼型型线的多个点的柱坐标,再依据轮三维翼型的轴坐标在轴面投影图上的形成的z=f(r)函数,分别计算变化后的前盖板二维翼型型线中间二维翼型型线以及后盖板二维翼型型线各点的Z轴坐标,此时,得到变化后的前盖板三维翼型型线、中间三维翼型型线以及后盖板三维翼型型线;
步骤7、将步骤6得到的变化后的前盖板三维翼型型线、中间三维翼型型线以及后盖板三维翼型型线,输入到几何造型软件中,形成变化后的叶轮三维形状;
步骤8、将步骤7得到的变化后的叶轮三维形状进行单流道网格划分,并进行水利性能测算,将测算得到的性能数据在全三维CFD软件中,生成全三维CFD粘性计算结果作为评价指标进行评价;
步骤9、以多目标遗传算法NSGA-II遗传算法为优化工具,选取多组变化参数C1,C2,C3,C4,C5,C6,Δrs,Δrm,Δrh,Δθs,Δθm,Δθh,Δθe,并重复步骤4至步骤8,得到对应的多个变化后的叶轮的水利性能参数,以叶轮叶片的效率和叶轮叶片上的最低静压值作为优化目标,对离心泵叶轮形状参数进行多目标优化设计,最终选取水利性能最佳的叶轮作为优化后的叶轮。
其中,步骤1的具体步骤为:由原始叶轮的轴面投影图得到前盖板流线ls、中间流线lm以及后盖板流线lh,将前盖板流线ls、中间流线lm以及后盖板流线lh分别绕叶轮旋转轴旋转得到前盖板流面、中间流面以及后盖板流面;原始叶轮与得到的前盖板流面、中间流面和后盖板流面相交得到相对应的三维翼型,再投影到周向平面上,即得到前盖板二维翼型型线中间二维翼型型线以及后盖板二维翼型型线沿所述前盖板二维翼型型线作一系列内切圆,连接内切圆圆心得到的光滑曲线为前盖板二维翼型骨线Gs;沿所述中间二维翼型型线作一系列内切圆,连接内切圆圆心得到的光滑曲线为中间二维翼型骨线Gm;沿所述后盖板二维翼型型线作一系列内切圆,连接内切圆圆心得到的光滑曲线为后盖板二维翼型骨线Gh。
其中,步骤2中,根据所述前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线Gh上的已知的多个点坐标,利用以下公式构成方程组,求各骨线对应的四个控制点P0,P1,P2,P3的坐标:
x(t)=(1-t)3x0+3t(1-t)2x1+3t2(1-t)x2+t3x3,
y(t)=(1-t)3y0+3t(1-t)2y1+3t2(1-t)y2+t3y3,
其中,t为四点三次贝塞尔曲线的参数,t∈[0,1],(x(t),y(t))为该四点三次贝塞尔曲线上任一点的坐标值,(x0,y0)为对应的控制点P0的坐标值,(x1,y1)为对应的控制点P1的坐标值,(x2,y2)为对应的控制点P2的坐标值(x3,y3)为对应的控制点P3的坐标值。
其中,所述步骤3.3中,定义前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线Gh的控制点P0的θ坐标的变化量:Δθs=Δθm=Δθh=Δθ。
其中,步骤5加厚的具体方法为:以变化后的前盖板二维翼型骨线GS的多个点为圆心,以已知的叶轮翼型厚度分布值为半径作一系列圆,连接所作圆的外部轮廓线得到的光滑曲线为变化后的前盖板二维翼型型线依照上述方法,根据变化后的中间二维翼型骨线Gm和后盖板二维翼型骨线Gh,分别得到变化后的中间二维翼型型线和后盖板二维翼型型线
其中,步骤6的具体步骤为:
由原始叶轮的轴面投影图得到前盖板流线ls、中间流线lm以及后盖板流线lh;
在前盖板流线ls上取均匀分布的多个点,并得到该多个取点的柱坐标,根据三次样条差值函数拟合得出z=f(r)函数,将变化后的前盖板二维翼型型线上各点柱坐标的r坐标代入z=f(r)函数,求得变化后的前盖板二维翼型型线各点的Z轴坐标;
本发明的有益效果是,优化参数构建合理,采用NSGA-II遗传算法为优化工具,采用本优化方法设计出来的离心泵叶轮具有效率高,空化空蚀性能好的特点,相比于现在依靠设计人员经验来修改叶片形状来进行优化的优化效果好。另一方面给,由于在优化设计过程中不需要依赖,全部过程由计算机自动完成,节省了人力物力。
附图说明
图1为原始叶轮的前盖板流线ls、中间流线lm以及后盖板流线lh的示意图;
图3为对图2中各二维翼型型线进行骨化后得到的前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线Gh;
图4为四点三次贝塞尔曲线中骨线各控制点的位置分布示意图;
图5为变化后的叶轮三维形状进行单流道网格划分的示意图。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
本发明一种抗空蚀离心泵叶轮优化设计方法,将离心泵叶轮叶片三维叶片形状采用Bézier曲线(即贝塞尔曲线)参数化技术进行参数化,将前盖板流面、中间流面以及后盖板流面分别与叶片相交得到的三维翼型投影到周向XY平面上,然后采用三次Bézier曲线将周向XY平面上的翼型骨线参数化,通过控制翼型骨线的变化来控制平面上翼型的变化。周向XY平面上的翼型变化后,然后通过柱坐标系下三个流面的r坐标和z坐标的函数关系z=f(r),可将平面翼型变换为三维翼型。然后,以叶轮叶片的全三维CFD粘性计算结果为评价指标,以NSGA-II遗传算法为优化工具,以叶轮叶片的效率和叶片上的最低静压值作为优化目标对离心泵叶轮叶片形状参数进行多目标优化设计,从而得出效率和空化空蚀性能均优的离心泵叶轮。
本发明方法包括以下步骤:
步骤1、得到原始叶轮的周向XY平面上的前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线Gh。
其中,步骤1的具体步骤为:由原始叶轮的轴面投影图得到前盖板流线ls、中间流线lm以及后盖板流线lh,将前盖板流线ls、中间流线lm以及后盖板流线lh分别绕叶轮旋转轴旋转得到前盖板流面、中间流面以及后盖板流面;原始叶轮与得到的前盖板流面、中间流面和后盖板流面相交得到相对应的三维翼型,将得到的三个三维翼型投影到周向平面上即得到前盖板二维翼型型线中间二维翼型型线以及后盖板二维翼型型线沿所述前盖板二维翼型型线作一系列内切圆,连接内切圆圆心得到的光滑曲线为前盖板二维翼型骨线Gs;沿所述中间二维翼型型线作一系列内切圆,连接内切圆圆心得到的光滑曲线为中间二维翼型骨线Gm;沿所述后盖板二维翼型型线作一系列内切圆,连接内切圆圆心得到的光滑曲线为后盖板二维翼型骨线Gh。
步骤2、采用四点三次贝塞尔曲线分别对上述前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线Gh进行参数化,以得到各骨线的四个控制点P0,P1,P2,P3的在周向XY平面上的坐标,其中控制点P0为首点,控制点P3为末点,控制点P1和P2为中间点。
Bezier曲线表达式如下: 0≤t≤1,其中,Pi(i=0,1,2,L L,n)为折线多边形的顶点;Bi,n(t)称为伯恩斯坦基函数, (i=0,1,2,L L,n)。本发明方法选用四点三次贝塞尔曲线。
其中,根据所述前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线Gh上的已知的多个点坐标,利用以下公式构成方程组,求各骨线对应的四个控制点P0,P1,P2,P3的坐标:
x(t)=(1-t)3x0+3t(1-t)2x1+3t2(1-t)x2+t3x3,
y(t)=(1-t)3y0+3t(1-t)2y1+3t2(1-t)y2+t3y3,
其中,t为四点三次贝塞尔曲线的参数,t∈[0,1],(x(t),y(t))为该四点三次贝塞尔曲线上任一点的坐标值,(x0,y0)为对应的控制点P0的坐标值,(x1,y1)为对应的控制点P1的坐标值,(x2,y2)为对应的控制点P2的坐标值(x3,y3)为对应的控制点P3的坐标值。
t=0时,贝塞尔曲线与特征多边形的第一条边P0P1相切;在末点处,即t=1的时候,贝塞尔曲线与特征多边形的最后一条边P2P3相切。根据上述原理,将周向XY平面的二维翼型骨线用Bézier曲线参数化方法进行参数化。
步骤3、构建叶轮优化参数:
叶轮优化参数的构建是本发明方法的入口,它在整个优化过程中作用非常重要。
步骤3.1、构建前盖板二维翼型骨线Gs的挠度参数C1和C2:
为了减少设计变量,在进行骨线参数化时应尽可能用最少的设计参数来表达骨线的形状。连接前盖板二维翼型骨线Gs的控制点P0和P1,以及控制点P2和P3,延长直线P0P1和直线P2P3并交于Q点,令: 则P1和P2点坐标满足: 其中,c1∈[0,1],c2∈[0,1]。因此,变化C1和C2的取值,就能控制骨线中间点P1和P2在直线P0Q和直线P1Q上的位置,即中间点P1和P2在周向XY平面上的坐标,也就能控制前盖板二维翼型骨线Gs的挠度。
步骤3.2、依照步骤3.1得到中间二维翼型骨线Gm的挠度参数C3和C4,以及后盖板二维翼型骨线Gh的挠度参数C5和C6;
步骤3.3、定义前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线Gh各自控制点P0和P3的变化量。骨线首点P0和骨线末点P3的参数设定,是为了通过控制该两个点参数的变化,来控制叶轮叶片进出口边的变化。
将前盖板二维翼型骨线Gs的控制点P0的柱坐标为(rs,θs),中间二维翼型骨线Gm的控制点P0的柱坐标为(rm,θm),后盖板二维翼型骨线Gh的控制点P0的柱坐标为(rh,θh);定义各控制点P0的r坐标的优化变化量为Δrs,Δrm,Δrh,定义各控制点P0的θ的坐标变化量为Δθs,Δθm,Δθh;
前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线Gh的控制点P3点均重合,且该P3点的柱坐标为(re,θe),定义该P3点的θ的坐标变化量为Δθe。
在不影响优化效果的前提下,为了减少系统的运算量,步骤3.3中,定义前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线Gh的控制点P0的θ坐标的变化量:Δθs=Δθm=Δθh=Δθ。
步骤4、构建变化后的前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线Gh:
依照上述方法,变化C3和C4的取值,得到中间二维翼型骨线Gm变化后的控制点P1和P2在周向XY平面上的坐标;变化C5和C6的取值,得到后盖板二维翼型骨线Gh变化后的控制点P1和P2在周向XY平面上的坐标。
步骤4.2、变化Δrs,Δrm,Δrh的取值,以及Δθs,Δθm,Δθh的取值,分别得到前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线Gh变化后的控制点P0的柱坐标值,根据柱坐标系与平面直角坐标系的转换关系,分别得到前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线Gh变化后的控制点P0在周向XY平面上的坐标。
步骤4.3、变化Δθe的取值,分别得到前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线Gh变化后的的控制点P3点的柱坐标值,根据柱坐标系与平面直角坐标系的转换关系,得到前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线Gh变化后的控制点P3点在周向XY平面上的坐标。
步骤4.4、根据变化后的前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线Gh变化后的四个控制点P0,P1,P2,P3的坐标,构建得到在周向XY平面上的变化后的前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线Gh。
步骤5、将步骤4得到变化后的前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线Gh,按照翼型厚度分布规律进行加厚,分别得到变化后的前盖板二维翼型型线中间二维翼型型线以及后盖板二维翼型型线
其中,对骨线进行加厚的具体方法为:以变化后的前盖板二维翼型骨线Gs的多个点为圆心,以已知的叶轮翼型厚度分布值为半径作一系列圆,连接所作圆的外部轮廓线得到的光滑曲线为变化后的前盖板二维翼型型线依照上述方法,根据变化后的中间二维翼型骨线Gm和后盖板二维翼型骨线Gh,分别得到变化后的中间二维翼型型线和后盖板二维翼型型线
步骤6、求得步骤5得到的变化后的前盖板二维翼型型线中间二维翼型型线以及后盖板二维翼型型线的多个点的柱坐标,再依据轮三维翼型的轴坐标在轴面投影图上的形成的z=f(r)函数,分别计算变化后的前盖板二维翼型型线中间二维翼型型线以及后盖板二维翼型型线各点的Z轴坐标,此时,得到变化后的前盖板三维翼型型线、中间三维翼型型线以及后盖板三维翼型型线。
其中,具体步骤为:
由原始叶轮的轴面投影图得到前盖板流线ls、中间流线lm以及后盖板流线lh;
在前盖板流线ls上取均匀分布的多个点,并得到该多个取点的柱坐标,根据三次样条差值函数拟合得出z=f(r)函数,将变化后的前盖板二维翼型型线上各点柱坐标的r坐标代入z=f(r)函数,求得变化后的前盖板二维翼型型线各点的Z轴坐标;
步骤7、将步骤6得到的变化后的前盖板三维翼型型线、中间三维翼型型线以及后盖板三维翼型型线,输入到几何造型软件中,形成变化后的叶轮三维形状。
步骤8、将步骤7得到的变化后的叶轮三维形状进行单流道网格划分,并进行水利性能测算,将测算得到的性能数据在全三维CFD软件中,生成全三维CFD粘性计算结果作为评价指标进行评价。
网格生成应以能满足解决实际问题的需求为目标。离心泵叶轮内部的流道是一个非常复杂的几何形状,为了节约计算机资源,在本方法中对叶轮的流道采用块结构化网格进行划分。同时为了节约单个优化点的计算时间,在对叶轮性能进行评估时只选取了单个叶轮流道进行计算,因此块结构化网格划分也只对叶轮单流道进行。
此步骤的优化过程中,对变化后的叶轮三维形状的全三维粘性CFD计算采用商业CFD程序CFX11.0进行,过程采用批处理加载CFX宏命令实现前处理快速加载、求解器后台运行和后处理自动成叶片性能数据并输出。
步骤9、以多目标遗传算法NSGA-II遗传算法为优化工具,选取多组变化参数C1,C2,C3,C4,C5,C6,Δrs,Δrm,Δrh,Δθs,Δθm,Δθh,Δθe,并重复步骤4至步骤8,得到对应的多个变化后的叶轮的水利性能参数,以叶轮叶片的效率和叶轮叶片上的最低静压值作为优化目标,对离心泵叶轮形状参数进行多目标优化设计,最终选取水力性能最佳的叶轮作为优化后的叶轮。
以智能算法NSGA-II遗传算法代替人工来完成对叶片性能优劣的判定以及对叶片形状的修改,这相比传统的优化设计方法,把设计者从不断的重复劳动中解放了出来,而且采用NSGA-II遗传算法来对这种高度非线性、多峰值的优化问题来进行寻优,鲁棒性更强,也更不易陷入局部最优的情况。
叶轮效率值代表着叶轮的能量性能,而叶片表面上的最低静压值代表着离心泵叶轮的空蚀性能,提高叶片表面上的最低静压值可以改善叶轮的空蚀性能。通过这两个优化目标,即能控制叶轮的能量和空蚀性能。通过NSGA-II遗传算法来优化提高这两个目标函数值,即可达到提高离心泵叶轮能量性能和空蚀性能的效果。
具体实施例
采用本发明方法对某原始叶轮外径为720mm,叶片数为9,转速为960r/min,流量设计为440L/s的离心泵叶轮进行了优化设计。
步骤1、得到原始叶轮的周向XY平面上的前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线Gh。其中,如图1所示,为原始叶轮的轴面投影图得到前盖板流线ls、中间流线lm以及后盖板流线lh。如图2所示,为原始叶轮的前盖板二维翼型型线中间二维翼型型线以及后盖板二维翼型型线如图3所示,为对图2中各二维翼型型线进行骨化后得到的前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线Gh。
步骤2、采用四点三次贝塞尔曲线分别对上述前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线Gh进行参数化,以得到各骨线的四个控制点P0,P1,P2,P3的在周向XY平面上的坐标,如图4所示,控制点P0为首点,控制点P3为末点,控制点P1和P2为中间点。
步骤3、构建叶轮优化参数:
构建前盖板二维翼型骨线Gs的挠度参数C1和C2,中间二维翼型骨线Gm的挠度参数C3和C4,以及后盖板二维翼型骨线Gh的挠度参数C5和C6。定义前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线Gh各自控制点P0和P3的变化量:定义三个控制点P0的r坐标的优化变化量为Δrs,Δrm,Δrh,定义三个控制点P0的θ坐标的变化量Δθs=Δθm=Δθh=Δθ,定义重合P3点的θ的坐标变化量为Δθe。
步骤4、构建变化后的前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线Gh。
步骤5、将步骤4得到变化后的前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线Gh,按照翼型厚度分布规律进行加厚,分别得到变化后的前盖板二维翼型型线中间二维翼型型线以及后盖板二维翼型型线
其中,步骤6的具体步骤为:
由原始叶轮的轴面投影图得到前盖板流线ls、中间流线lm以及后盖板流线lh;
在前盖板流线ls上取均匀分布的100个点,并得到该100个取点的柱坐标,根据三次样条差值拟合得出z=f(r)函数,将变化后的前盖板二维翼型型线上各点柱坐标的r坐标代入z=f(r)函数,求得变化后的前盖板二维翼型型线各点的Z轴坐标。
步骤7、将步骤6得到的变化后的前盖板三维翼型型线、中间三维翼型型线以及后盖板三维翼型型线,输入到几何造型软件中,形成变化后的叶轮三维形状.
步骤8、如图5所示,将步骤7得到的变化后的叶轮三维形状进行单流道网格划分,并进行水利性能测算,将测算得到的性能数据在全三维CFD软件中,生成全三维CFD粘性计算结果作为评价指标进行评价。
步骤9、以多目标遗传算法NSGA-II遗传算法为优化工具,选取多组变化参数C1,C2,C3,C4,C5,C6,Δrs,Δrm,Δrh,Δθ,Δθe,并重复步骤4至步骤8,得到对应的多个变化后的叶轮的水力性能参数,以叶轮叶片的效率和叶轮叶片上的最低静压值作为优化目标,对离心泵叶轮形状参数进行多目标优化设计,最终选取水利性能最佳的叶轮作为优化后的叶轮。
考虑到该原始叶轮的实际运行条件以及其主要用途,在优化时选取1.1Qr(Qr为设计流量)工况点进行优化,
本实施例中,采用的NSGA-II算法的变异概率为0.8,种群数为40,遗传代数取为35代。为了使叶片在优化过程中避免出现非常规形状,设定C1,C2,C3,C4,C5,C6,Δrs,Δrm,Δrh,Δθ,Δθe各参数约束如下表所示:
参数 | 上限 | 下限 |
C1(无量纲数) | 0.45 | 0.7 |
C2(无量纲数) | 0.3 | 0.7 |
C3(无量纲数) | 0.4 | 0.7 |
C4(无量纲数) | 0.4 | 0.7 |
C5(无量纲数) | 0.35 | 0.7 |
C6(无量纲数) | 0.25 | 0.7 |
Δrm(mm) | -5.0 | 5.0 |
Δrh(mm) | -8.0 | 2.0 |
Δrs(mm) | 0 | 5.0 |
Δθ(°) | 0 | 5.0 |
Δθe(°) | 0 | 5.0 |
如下表所示为,优化后的各参数与优化值的对比表:
如下表所示为,优化前后叶轮叶片的效率对比表。
0.7Q | 0.8Q | 0.9Q | Q | 1.1Q | 1.2Q | 1.3Q | 1.4Q | |
流量(L/s) | 308 | 352 | 396 | 440 | 484 | 528 | 572 | 616 |
优化前η% | 0.961 | 0.962 | 0.959 | 0.953 | 0.944 | 0.931 | 0.913 | 0.892 |
优化后η% | 0.957 | 0.96 | 0.961 | 0.96 | 0.957 | 0.951 | 0.943 | 0.932 |
如下表所示为,优化前后叶轮叶片上的最低静压值对比表。
如下表所示为,优化前后叶轮叶片扬程对比表。
0.7Q | 0.8Q | 0.9Q | Q | 1.1Q | 1.2Q | 1.3Q | 1.4Q | |
流量(L/s) | 308 | 352 | 396 | 440 | 484 | 528 | 572 | 616 |
优化前H(m) | 91.94 | 88.76 | 85.86 | 82.96 | 79.4 | 75.1 | 71.1 | 67.7 |
优化后H(m) | 91.74 | 88.69 | 85.48 | 82.12 | 78.74 | 76.4 | 73.5 | 69.1 |
由此可以看出,采用本发明方法优化后的离心泵叶轮在设计流量工况到大流量工况的工作区间内,效率及空化性能都有了较大的提高,扬程变化很小。本实施例中的离心泵是用于抽水灌溉,因此其工作时间大部分会处于设计流量及大流量工作区间,工作在小流量工况的时间很少,而本发明方法所采用的优化方案使得该离心泵在其运行区间内的能量性能和空化性能得到较大的提升,满足了实际工作的需要。
Claims (6)
1.一种抗空蚀离心泵叶轮优化设计方法,其特征在于,包括以下步骤:
步骤1、得到原始叶轮的周向XY平面上的前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线Gh;
步骤2、采用四点三次贝塞尔曲线分别对上述前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线Gh进行参数化,以得到各骨线的四个控制点P0,P1,P2,P3在周向XY平面上的坐标,其中,控制点P0为首点,控制点P3为末点,控制点P1和P2为中间点;
步骤3、构建叶轮优化参数:
步骤3.1、构建前盖板二维翼型骨线Gs的挠度参数C1和C2:
步骤3.2、依照步骤3.1得到中间二维翼型骨线Gm的挠度参数C3和C4,以及后盖板二维翼型骨线Gh的挠度参数C5和C6;
步骤3.3、定义前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线Gh各自控制点P0和P3的变化量:
将前盖板二维翼型骨线Gs的控制点P0的柱坐标为(rs,θs),中间二维翼型骨线Gm的控制点P0的柱坐标为(rm,θm),后盖板二维翼型骨线Gh的控制点P0的柱坐标为(rh,θh);定义各控制点P0的r坐标的优化变化量为Δrs,Δrm,Δrh,定义各控制点P0的θ的坐标变化量为Δθs,Δθm,Δθh;
所述前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线Gh的控制点P3点均重合,且该P3点的柱坐标为(re,θe),定义该P3点的θ的坐标变化量为Δθe;
步骤4、构建变化后的前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线Gh:
依照上述方法,变化C3和C4的取值,得到中间二维翼型骨线Gm变化后的控制点P1和P2在周向XY平面上的坐标;变化C5和C6的取值,得到后盖板二维翼型骨线Gh变化后的控制点P1和P2在周向XY平面上的坐标;
步骤4.2、变化Δrs,Δrm,Δrh的取值,以及Δθs,Δθm,Δθh的取值,分别得到前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线Gh变化后的控制点P0的柱坐标值,根据柱坐标系与平面直角坐标系的转换关系,分别得到前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线Gh变化后的控制点P0在周向XY平面上的坐标;
步骤4.3、变化Δθe的取值,分别得到前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线变化后的Gh的控制点P3点的柱坐标值,根据柱坐标系与平面直角坐标系的转换关系,得到前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线Gh变化后的控制点P3点在周向XY平面上的坐标;
步骤4.4、根据变化后的前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线Gh变化后的四个控制点P0,P1,P2,P3的坐标,构建得到在周向XY平面上的变化后的前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线Gh;
步骤5、将步骤4得到变化后的前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线Gh,按照翼型厚度分布规律进行加厚,分别得到变化后的前盖板二维翼型型线中间二维翼型型线以及后盖板二维翼型型线
步骤6、求得步骤5得到的变化后的前盖板二维翼型型线中间二维翼型型线以及后盖板二维翼型型线的多个点的柱坐标,再依据轮三维翼型的轴坐标在轴面投影图上的形成的z=f(r)函数,分别计算变化后的前盖板二维翼型型线中间二维翼型型线以及后盖板二维翼型型线各点的Z轴坐标,此时,得到变化后的前盖板三维翼型型线、中间三维翼型型线以及后盖板三维翼型型线;
步骤7、将步骤6得到的变化后的前盖板三维翼型型线、中间三维翼型型线以及后盖板三维翼型型线,输入到几何造型软件中,形成变化后的叶轮三维形状;
步骤8、将步骤7得到的变化后的叶轮三维形状进行单流道网格划分,并进行水利性能测算,将测算得到的性能数据在全三维CFD软件中,生成全三维CFD粘性计算结果作为评价指标进行评价;
步骤9、以多目标遗传算法NSGA-II遗传算法为优化工具,选取多组变化参数C1,C2,C3,C4,C5,C6,Δrs,Δrm,Δrh,Δθs,Δθm,Δθh,Δθe,并重复步骤4至步骤8,得到对应的多个变化后的叶轮的水利性能参数,以叶轮叶片的效率和叶轮叶片上的最低静压值作为优化目标,对离心泵叶轮参数进行多目标优化设计,选取水利性能最佳的一组叶轮参数作为优化后的叶轮参数。
2.按照权利要求1所述抗空蚀离心泵叶轮优化设计方法,其特征在于,其中,步骤1的具体步骤为:由原始叶轮的轴面投影图得到前盖板流线ls、中间流线lm以及后盖板流线lh,将前盖板流线ls、中间流线lm以及后盖板流线lh分别绕叶轮旋转轴旋转得到前盖板流面、中间流面以及后盖板流面;原始叶轮与得到的前盖板流面、中间流面和后盖板流面相交得到相对应的三维翼型,再投影到周向平面上,即得到前盖板二维翼型型线中间二维翼型型线以及后盖板二维翼型型线沿所述前盖板二维翼型型线作一系列内切圆,连接内切圆圆心得到的光滑曲线为前盖板二维翼型骨线Gs;沿所述中间二维翼型型线作一系列内切圆,连接内切圆圆心得到的光滑曲线为中间二维翼型骨线Gm;沿所述后盖板二维翼型型线作一系列内切圆,连接内切圆圆心得到的光滑曲线为后盖板二维翼型骨线Gh。
3.按照权利要求1所述抗空蚀离心泵叶轮优化设计方法,其特征在于,其中,步骤2中,根据所述前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线Gh上的已知的多个点坐标,利用以下公式构成方程组,求各骨线对应的四个控制点P0,P1,P2,P3的坐标:
x(t)=(1-t)3x0+3t(1-t)2x1+3t2(1-t)x2+t3x3,
y(t)=(1-t)3y0+3t(1-t)2y1+3t2(1-t)y2+t3y3,
其中,t为四点三次贝塞尔曲线的参数,t∈[0,1],(x(t),y(t))为该四点三次贝塞尔曲线上任一点的坐标值,(x0,y0)为对应的控制点P0的坐标值,(x1,y1)为对应的控制点P1的坐标值,(x2,y2)为对应的控制点P2的坐标值(x3,y3)为对应的控制点P3的坐标值。
4.按照权利要求1所述抗空蚀离心泵叶轮优化设计方法,其特征在于,其中,所述步骤3.3中,定义前盖板二维翼型骨线Gs、中间二维翼型骨线Gm以及后盖板二维翼型骨线Gh的控制点P0的θ坐标的变化量:Δθs=Δθm=Δθh=Δθ。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011102025245A CN102251983B (zh) | 2011-07-20 | 2011-07-20 | 一种抗空蚀离心泵叶轮优化设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011102025245A CN102251983B (zh) | 2011-07-20 | 2011-07-20 | 一种抗空蚀离心泵叶轮优化设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102251983A true CN102251983A (zh) | 2011-11-23 |
CN102251983B CN102251983B (zh) | 2013-03-13 |
Family
ID=44979477
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2011102025245A Expired - Fee Related CN102251983B (zh) | 2011-07-20 | 2011-07-20 | 一种抗空蚀离心泵叶轮优化设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102251983B (zh) |
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102750414A (zh) * | 2012-06-25 | 2012-10-24 | 江苏大学 | 一种验证离心泵网格质量与计算精度关系的方法 |
CN103016396A (zh) * | 2012-12-11 | 2013-04-03 | 江苏大学 | 一种靠汽蚀控制最大流量的离心泵水力设计方法 |
CN103452900A (zh) * | 2013-09-18 | 2013-12-18 | 成都市永益泵业有限公司 | 一种可提高泵效率的离心叶轮几何参数优化设计方法 |
CN103473422A (zh) * | 2013-09-18 | 2013-12-25 | 成都市永益泵业有限公司 | 基于奇点分布法的轴流叶轮翼型设计方法 |
CN103527514A (zh) * | 2013-10-29 | 2014-01-22 | 中国北车集团大连机车研究所有限公司 | 机车用离心通风机叶轮 |
CN104613157A (zh) * | 2015-01-26 | 2015-05-13 | 北京理工大学 | 基于贝塞尔曲线的液力变矩器二维叶片型线构造方法 |
CN105240310A (zh) * | 2015-12-03 | 2016-01-13 | 浙江工业大学之江学院 | 一种双叶片无堵塞离心叶轮水力设计方法 |
CN106227967A (zh) * | 2016-08-01 | 2016-12-14 | 杭州汽轮机股份有限公司 | 工业汽轮机低压级组叶片型线优化方法 |
RU2613545C1 (ru) * | 2015-12-02 | 2017-03-17 | федеральное государственное бюджетное образовательное учреждение высшего образования "Национальный исследовательский университет "МЭИ" | Реактивное рабочее колесо центробежного насоса |
CN106640210A (zh) * | 2016-09-14 | 2017-05-10 | 西安交通大学苏州研究院 | 一种径‑轴流式透平膨胀机叶轮叶片三维型线的设计方法 |
CN107066686A (zh) * | 2017-02-22 | 2017-08-18 | 江苏大学 | 一种基于遗传算法的轴流泵叶轮水力优化设计方法 |
CN107529643A (zh) * | 2017-09-15 | 2018-01-02 | 西安理工大学 | 一种兼顾水力、强度性能的混流式转轮多工况优化方法 |
CN108446452A (zh) * | 2018-02-27 | 2018-08-24 | 江苏大学 | 一种混流泵叶轮鲁棒优化设计方法 |
CN111734674A (zh) * | 2020-04-26 | 2020-10-02 | 上海凯泉泵业(集团)有限公司 | 一种基于遗传算法的离心泵多工况节能优化方法 |
CN111852937A (zh) * | 2020-06-02 | 2020-10-30 | 江苏大学 | 一种离心空化泵 |
CN112214850A (zh) * | 2020-09-29 | 2021-01-12 | 西安理工大学 | 一种双向全贯流式转轮叶片几何形状的参数化控制方法 |
CN112528427A (zh) * | 2020-11-18 | 2021-03-19 | 江苏徐工工程机械研究院有限公司 | 一种叶轮机械叶片骨线设计方法 |
CN112541234A (zh) * | 2020-11-13 | 2021-03-23 | 江苏徐工工程机械研究院有限公司 | 一种叶轮机械非等厚叶片设计方法 |
CN113221264A (zh) * | 2021-04-21 | 2021-08-06 | 山东双轮股份有限公司 | 海水淡化泵流道式导叶结构设计优化方法 |
CN115994394A (zh) * | 2023-02-21 | 2023-04-21 | 上海中韩杜科泵业制造有限公司 | 离心泵叶轮造型方法、装置以及设备 |
CN117539168A (zh) * | 2024-01-09 | 2024-02-09 | 江西江投能源技术研究有限公司 | 一种基于半实物仿真的水轮机空化诊断系统及方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6796770B2 (en) * | 2002-11-06 | 2004-09-28 | Spx Corporation | Impeller and method using solid free form fabrication |
CN101520052A (zh) * | 2008-10-28 | 2009-09-02 | 兰州理工大学 | 一种离心泵叶轮的逆向设计方法 |
-
2011
- 2011-07-20 CN CN2011102025245A patent/CN102251983B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6796770B2 (en) * | 2002-11-06 | 2004-09-28 | Spx Corporation | Impeller and method using solid free form fabrication |
CN101520052A (zh) * | 2008-10-28 | 2009-09-02 | 兰州理工大学 | 一种离心泵叶轮的逆向设计方法 |
Non-Patent Citations (2)
Title |
---|
刘小民等: "采用遗传算法的离心叶轮多目标自动优化设计", 《西安交通大学学报》 * |
郭涛等: "基于遗传算法的离心泵叶轮参数化造型及优化设计", 《排灌机械工程学报》 * |
Cited By (33)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102750414B (zh) * | 2012-06-25 | 2015-12-09 | 江苏大学 | 一种验证离心泵网格质量与计算精度关系的方法 |
CN102750414A (zh) * | 2012-06-25 | 2012-10-24 | 江苏大学 | 一种验证离心泵网格质量与计算精度关系的方法 |
CN103016396A (zh) * | 2012-12-11 | 2013-04-03 | 江苏大学 | 一种靠汽蚀控制最大流量的离心泵水力设计方法 |
CN103452900A (zh) * | 2013-09-18 | 2013-12-18 | 成都市永益泵业有限公司 | 一种可提高泵效率的离心叶轮几何参数优化设计方法 |
CN103473422A (zh) * | 2013-09-18 | 2013-12-25 | 成都市永益泵业有限公司 | 基于奇点分布法的轴流叶轮翼型设计方法 |
CN103452900B (zh) * | 2013-09-18 | 2015-06-10 | 成都市永益泵业有限公司 | 一种可提高泵效率的离心叶轮几何参数优化设计方法 |
CN103527514A (zh) * | 2013-10-29 | 2014-01-22 | 中国北车集团大连机车研究所有限公司 | 机车用离心通风机叶轮 |
CN104613157A (zh) * | 2015-01-26 | 2015-05-13 | 北京理工大学 | 基于贝塞尔曲线的液力变矩器二维叶片型线构造方法 |
RU2613545C1 (ru) * | 2015-12-02 | 2017-03-17 | федеральное государственное бюджетное образовательное учреждение высшего образования "Национальный исследовательский университет "МЭИ" | Реактивное рабочее колесо центробежного насоса |
CN105240310B (zh) * | 2015-12-03 | 2017-09-26 | 浙江工业大学之江学院 | 一种双叶片无堵塞离心叶轮水力设计方法 |
CN105240310A (zh) * | 2015-12-03 | 2016-01-13 | 浙江工业大学之江学院 | 一种双叶片无堵塞离心叶轮水力设计方法 |
CN106227967A (zh) * | 2016-08-01 | 2016-12-14 | 杭州汽轮机股份有限公司 | 工业汽轮机低压级组叶片型线优化方法 |
CN106640210B (zh) * | 2016-09-14 | 2019-03-05 | 西安交通大学苏州研究院 | 一种径-轴流式透平膨胀机叶轮叶片三维型线的设计方法 |
CN106640210A (zh) * | 2016-09-14 | 2017-05-10 | 西安交通大学苏州研究院 | 一种径‑轴流式透平膨胀机叶轮叶片三维型线的设计方法 |
CN107066686B (zh) * | 2017-02-22 | 2020-09-25 | 江苏大学 | 一种基于遗传算法的轴流泵叶轮水力优化设计方法 |
CN107066686A (zh) * | 2017-02-22 | 2017-08-18 | 江苏大学 | 一种基于遗传算法的轴流泵叶轮水力优化设计方法 |
CN107529643A (zh) * | 2017-09-15 | 2018-01-02 | 西安理工大学 | 一种兼顾水力、强度性能的混流式转轮多工况优化方法 |
CN108446452A (zh) * | 2018-02-27 | 2018-08-24 | 江苏大学 | 一种混流泵叶轮鲁棒优化设计方法 |
CN108446452B (zh) * | 2018-02-27 | 2019-08-02 | 江苏大学 | 一种混流泵叶轮鲁棒优化设计方法 |
CN111734674A (zh) * | 2020-04-26 | 2020-10-02 | 上海凯泉泵业(集团)有限公司 | 一种基于遗传算法的离心泵多工况节能优化方法 |
CN111852937A (zh) * | 2020-06-02 | 2020-10-30 | 江苏大学 | 一种离心空化泵 |
CN111852937B (zh) * | 2020-06-02 | 2022-01-11 | 江苏大学 | 一种离心空化泵 |
CN112214850B (zh) * | 2020-09-29 | 2021-12-17 | 西安理工大学 | 一种双向全贯流式转轮叶片几何形状的参数化控制方法 |
CN112214850A (zh) * | 2020-09-29 | 2021-01-12 | 西安理工大学 | 一种双向全贯流式转轮叶片几何形状的参数化控制方法 |
CN112541234A (zh) * | 2020-11-13 | 2021-03-23 | 江苏徐工工程机械研究院有限公司 | 一种叶轮机械非等厚叶片设计方法 |
CN112541234B (zh) * | 2020-11-13 | 2023-10-03 | 江苏徐工工程机械研究院有限公司 | 一种叶轮机械非等厚叶片设计方法 |
CN112528427A (zh) * | 2020-11-18 | 2021-03-19 | 江苏徐工工程机械研究院有限公司 | 一种叶轮机械叶片骨线设计方法 |
CN112528427B (zh) * | 2020-11-18 | 2023-05-05 | 江苏徐工工程机械研究院有限公司 | 一种叶轮机械叶片骨线设计方法 |
CN113221264A (zh) * | 2021-04-21 | 2021-08-06 | 山东双轮股份有限公司 | 海水淡化泵流道式导叶结构设计优化方法 |
CN115994394A (zh) * | 2023-02-21 | 2023-04-21 | 上海中韩杜科泵业制造有限公司 | 离心泵叶轮造型方法、装置以及设备 |
CN115994394B (zh) * | 2023-02-21 | 2023-11-17 | 上海中韩杜科泵业制造有限公司 | 离心泵叶轮造型方法、装置以及设备 |
CN117539168A (zh) * | 2024-01-09 | 2024-02-09 | 江西江投能源技术研究有限公司 | 一种基于半实物仿真的水轮机空化诊断系统及方法 |
CN117539168B (zh) * | 2024-01-09 | 2024-03-26 | 江西江投能源技术研究有限公司 | 一种基于半实物仿真的水轮机空化诊断系统及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102251983B (zh) | 2013-03-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102251983B (zh) | 一种抗空蚀离心泵叶轮优化设计方法 | |
CN101520052B (zh) | 一种离心泵叶轮的逆向设计方法 | |
CN102608914B (zh) | 径流式液力透平优化设计方法 | |
CN108153998B (zh) | 离心鼓风机叶轮的全三维气动优化设计方法 | |
CN104166752B (zh) | 液力变矩器全流道瞬态数值模拟计算方法 | |
CN106650105A (zh) | 一种混流泵叶轮的设计方法 | |
CN105201916A (zh) | 一种空间导叶离心泵水力设计方法 | |
CN106446324B (zh) | 大型工业汽轮机末级扭叶片设计方法 | |
CN109858078A (zh) | 一种双叶片离心叶轮自寻优设计方法 | |
CN109598081B (zh) | 基于数据降维及多二维流面的径流式透平气动优化方法 | |
CN105805043A (zh) | 一种具有长短叶片特征的不可调轴流泵叶轮的设计方法 | |
CN110617238B (zh) | 一种离心泵叶轮的优化设计方法 | |
CN115017843A (zh) | 一种离心压缩机气动性能优化设计方法 | |
CN106949095A (zh) | 低压轴流式通风机叶轮叶片的优化方法 | |
CN107908914B (zh) | 离心压缩机闭式叶轮可加工性判断及中间截面的计算方法 | |
CN105971931A (zh) | 一种离心式叶轮分流叶片的设计方法 | |
CN112784370A (zh) | 一种多级离心泵的空间导叶的设计方法 | |
CN110245408B (zh) | 一种汽轮机单圆弧压力面叶型设计方法 | |
CN112487570A (zh) | 基于自由变形技术的离心压缩机分流叶片形状优化方法 | |
CN106202711B (zh) | 液力变矩器流道模型的参数化方法 | |
CN108446452B (zh) | 一种混流泵叶轮鲁棒优化设计方法 | |
CN104675713B (zh) | 一种基于数据样本的离心泵无过载设计方法 | |
CN109236726A (zh) | 一种高比转速轴流泵叶轮出口角和厚度设计方法 | |
CN103321957A (zh) | 轴流通风机叶轮叶片的优化设计方法 | |
CN107013490B (zh) | 一种低压轴流式通风机叶轮叶片的优化方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130313 Termination date: 20160720 |
|
CF01 | Termination of patent right due to non-payment of annual fee |