CN109194050A - 一种基于中心复合设计的转子结构优化方法 - Google Patents
一种基于中心复合设计的转子结构优化方法 Download PDFInfo
- Publication number
- CN109194050A CN109194050A CN201811032215.6A CN201811032215A CN109194050A CN 109194050 A CN109194050 A CN 109194050A CN 201811032215 A CN201811032215 A CN 201811032215A CN 109194050 A CN109194050 A CN 109194050A
- Authority
- CN
- China
- Prior art keywords
- design variable
- rotor structure
- central composite
- design
- model
- 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
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02K—DYNAMO-ELECTRIC MACHINES
- H02K15/00—Methods or apparatus specially adapted for manufacturing, assembling, maintaining or repairing of dynamo-electric machines
- H02K15/02—Methods or apparatus specially adapted for manufacturing, assembling, maintaining or repairing of dynamo-electric machines of stator or rotor bodies
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/06—Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Power Engineering (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Manufacturing & Machinery (AREA)
- Turbine Rotor Nozzle Sealing (AREA)
Abstract
本发明公开了一种基于中心复合设计的转子结构优化方法,包括步骤:1、确定转子结构优化的设计变量,确定设计变量的范围与约束条件;2、采用中心复合设计模型生成所需要的N组设计变量,对这N组样本进行参数化有限元计算,得到N组样本所对应的最大等效应力值作为目标函数;3、采用最小二乘法拟合设计变量与目标函数的二阶多项式回归模型;4、对所拟合的回归模型进行优化,给定初始点、初始步长、加速系数、收缩系数和精度,设置约束条件,由表决判断轴向探测移动和模式移动方式,直至满足收敛条件,完成对转子结构的优化。本发明能够快速、准确地解决转子结构优化的问题。
Description
技术领域
本发明属于汽轮机转子技术领域,特别涉及一种基于中心复合设计的转子结构优化方法。
背景技术
随着高参数大容量的火电机组的普遍采用,汽轮机机组的安全性和可靠性问题也日益突出。对于一直处于高温高压恶劣环境的汽轮机转子,其承受着离心力、气流力、交变热应力等多种载荷的共同作用。转子上的圆角、轴肩以及凹槽等形状突变部位将产生较大的热应力集中,这些部位的热应力有可能会超过材料对应温度下的屈服极限,产生较大的塑性变形,经过一定周次的交变热应力循环作用,在某些应力集中区域就会产生疲劳初始裂纹,而转子上的裂纹不易修复,在此后的循环载荷下,裂纹将进一步扩展直到转子发生断裂,从而造成巨大的经济损失,所以有必要改进转子结构以降低其应力水平,提高机组运行安全性。
将有限元法和优化算法相结合是目前主要采用的数值计算优化方法。通过有限元软件完成对转子热结构耦合分析模型的计算并提取其应力结果,然后使用优化算法通过“分析——评估——修正”的方式完成对最优转子结构的搜索,以达到降低应力水平的目的。这种优化方法具有优化结果精度高的优点,但由于热结构耦合分析模型的有限元网格数量巨大,优化过程需要耗费大量时间和计算机资源。
发明内容
本发明的目的是提供一种基于中心复合设计的转子结构优化方法,该方法通过构造目标函数与优化参数的回归模型,再对回归模型进行优化,获取最优设计的转子结构,因此,本发明能够快速、准确地解决转子结构优化的问题。
本发明采用如下技术方案来实现的:
一种基于中心复合设计的转子结构优化方法,该方法包括:首先确定转子结构优化过程中的设计变量,基于中心复合设计方法生成所需要的N组设计变量,参数化建立相应的有限元模型并计算获得目标函数值,再采用最小二乘法拟合设计变量与目标函数的回归模型,最后利用变步长的模式搜索法对回归模型进行优化,获得回归模型的最优组合。
本发明进一步的改进在于,具体包括以下步骤:
第一步,确定转子结构优化的设计变量,确定设计变量的范围与约束条件;
第二步,采用中心复合设计模型生成所需要的N组设计变量,对这N组样本进行参数化有限元计算,得到N组样本所对应的最大等效应力值作为目标函数;
第三步,采用最小二乘法拟合设计变量与目标函数的二阶多项式回归模型;
第四步,对所拟合的回归模型进行优化,给定初始点、初始步长、加速系数、收缩系数、精度,设置约束条件,由表决判断轴向探测移动和模式移动方式,直至满足收敛条件,完成对转子结构的优化。
本发明进一步的改进在于,步骤一中,选定其需要优化的结构,通过n个参数来确定其具体的几何形状,即x=(x1,x2,...,xn)为该转子的几何形状参数;同时,为了避免结构内部发生相互干涉,防止在建立几何面等操作时发生错误,确定设计变量的范围与约束条件,即Li≤xi≤Ui,i=1,2,...,n,Li与Ui分别表示变量参数xi上下极限值。
本发明进一步的改进在于,步骤二中,采用中心复合设计模型生成N组设计变量,N由中心复合设计模型中参数数量n决定,n由步骤一所确定;同时,对以上N组设计变量所确定的结构实现建模、载荷和边界条件的加载、求解计算和后处理分析过程,提取后处理结果中的最大等效应力值作为目标函数。
本发明进一步的改进在于,步骤三中,依据步骤二所得的N组设计变量与目标函数值,选用基于泰勒展开式的二阶多项式模型,对设计变量与目标函数进行多元回归拟合,得到表征最大等效应力值与n维设计变量非线性关系的回归模型;
回归模型形式如:其中f(x)为目标函数,即最大等效应力值,式中包括常数项β0、线性项线性交叉项和二次项βj、βij、βjj分别为线性项、线性交叉项和二次项系数;
采用最小二乘法对以上回归模型进行拟合。
本发明进一步的改进在于,步骤四中,具体包括:
401)由原转子结构模型,测定其几何参数以给定设计变量的初始值x(0)=(x1,x2,...,xn),给定初始步长δ(0)=(δ1,δ2,...,δn)T、加速系数γ、收缩系数θ和收敛精度ε;
402)优化的目标函数即为步骤三所拟合的二阶多项式f(x),x∈Rn,模式搜索的坐标方向n维空间的基矢向量e1=(1,0,...,0)T,e2=(0,1,...,0)T,...,en=(0,0,...,1)T,同时,令i=0,取y=x(i),令j=1;
403)做平行于单位矢量ej的探测移动循环:
正向探测移动:若则令j=j+1,否则做负向探测移动;
负向探测移动:若则令j=j+1,否则令y=y,j=j+1;
404)令x(i+1)=y,若f(x(i+1))<f(x(i)),表明模式移动成功,则对x(i+1)沿加速方向做模式移动,令y=x(i+1)+γ(x(i+1)-x(i)),δ(i+1)=θδ(i),i=i+1,进入循环403);若f(x(i+1))≥f(x(i)),表明模式移动失败,则进入循环405);
判断优化迭代是否收敛,若|δ(i)|<ε,则停止迭代,输出x(i),即为所求设计变量;否则,若x(i+1)=x(i),令y=x(i+1),δ(i+1)=θδ(i),i=i+1,进入循环403),若x(i+1)≠x(i),令y=x(i+1),δ(i+1)=δ(i),i=i+1,进入循环403)。
本发明具有如下有益的技术效果:
本发明所建立的转子结构优化方法中,采用中心复合设计仅通过很少的有限元计算便可以拟合出精度很高的目标函数与设计变量的复杂非线性回归模型,可以节省大量计算机资源,在此基础上对拟合的回归模型采用变步长的模式搜索法进行优化计算,可以提升优化效果。
附图说明
图1为本发明方法的总流程图;
图2为某示例转子结构图;其中图2(a)为转子整体结构;图2(b)为待优化局部结构;
图3为模式搜索算法的流程图。
具体实施方式
下面将结合附图详细说明本发明的实施方式。
请参阅图1所示,本发明提供的一种基于中心复合设计的转子结构优化方法,包括以下步骤:
一、确定转子结构优化的设计变量,确定设计变量的取值范围。
对于某种确定类型的转子,选定其需要优化的结构,往往需要n个参数来确定其具体的几何形状;同时,为了避免结构内部发生相互干涉,防止建立几何面等操作发生错误,需要确定设计变量的范围与约束条件。
以图2中的转子为例,由该转子的整体结构图和待优化局部结构图可知,该凹槽结构由n=5个参数来确定其几何形状:H1,H2为高度参数,控制着凹槽的深度;R1,R2,R3为半径参数,决定凹槽底部圆弧过渡的大小,结构的曲率半径变化显著影响其应力水平。转子凹槽具有一定深度,因此调节级右端面的H2不能低于左端面相应段的高度。同时,为了避免结构内部发生相互干涉,防止建立几何面等操作发生错误,各几何点的相对位置需要维持不变,确定的变量范围如表1所示。
表1设计变量的取值范围
一般的对于设计变量为x=(x1,x2,...,xn),其每个特征取值范围记为Li≤xi≤Ui,i=1,2,...,n,Li与Ui分别表示特征xi的上下极限值。
二、采用中心复合设计模型生成所需要的N组设计变量,建立这N组样本进行参数化有限元计算,得到N组样本所对应的最大等效应力值。
采用二水平全因子的中心复合表面设计,设计点由立方点、轴向点和中心点组成,由每个特征的三个取值水平(-1,0,1)所决定,即特征取值范围的最小值、中值和最大值。立方点由二水平对应的-1和+1点构成,共2n组;轴向点由一个坐标为-1或+1其余坐标为0的点组成,共2n个;中心点的坐标全为0,所以以此方式生成的设计变量组共2n+2n+1个。
以图2的转子凹槽结构为例,存在5个设计变量,各因素的取值水平如下表2所示。
表2二水平全因子中心复合设计下的变量取值
由图2的转子凹槽结构存在5个设计变量,所以采用以上设计方法共生成设计变量组N=2n+2n+1=43,安排如下表3所示。
表3二水平全因子中心复合设计的设计变量安排
采用以上二水平全因子的中心复合设计安排实验因素,对以上N组设计变量所确定的结构实现建模、载荷和边界条件的加载、求解计算和后处理等分析过程,提取后处理结果中的最大等效应力值即为目标函数,记为σ=(σ1,σ2,...,σN)T,N=43。
三、采用最小二乘法拟合设计变量与目标函数的二阶多项式回归模型。
回归模型形式如:其中f(x)为目标函数,即最大等效应力值,式中包括常数项β0、线性项线性交叉项和二次项βj、βij、βjj分别为线性项、线性交叉项和二次项系数。
采用最小二乘法对以上回归模型进行拟合,以图2所示的转子凹槽结构为例:
令x=(x1,x2,x3,x4,x5)分别为5个设计变量H1,H2,R1,R2,R3;
使用α0,α1,...,αm分别表示线性项、线性交叉项和二次项系数β0,βj(j=1,2,...,5),βij(i<j,j≤5),βjj(j=1,2,...,5),得α=(α0,α1,...,αm),m=20;
同理,使用v0,v1,...,vm分别表示1,xj(j=1,2,...,5),xixj(i<j,j≤5),xj 2(j=1,2,...,5),得v=(v0,v1,...,vm),m=20,则N组设计变量所对应的1,xkj(j=1,2,...,5),xkixkj(i<j,j≤5),xkj 2(j=1,2,...,5)记为vk=(vk0,vk1,...,vkm),k=1,2,...,N。
那么,可以改写为用最小二乘法来拟合就是要设法求出多项式f(x)的系数α0,α1,...,αm,使得拟合的多项式的值与已知步骤二中所求目标函数即最大应力值σ=(σ1,σ2,...,σN)T之差的平方和最小。
根据微积分中的极值原理可知所求Δ对α0,α1,...,αm的一阶偏导数为0,即:
化简得将代入,得
令则式(2)改写为
其中cij与dj由步骤二中所得设计变量和目标函数所确定,求解线性方程组(3)即可得α0,α1,...,αm,从而得到回归模型
四、优化回归模型,给定初始点、初始步长、加速系数、收缩系数、精度,设置约束条件,由完全表决判断轴向探测移动方向和模式移动,直至满足收敛条件,完成对转子结构的优化。
401)由原转子模型,测定其几何参数以给定设计变量的初始值x(0)=(x1,x2,...,xn),如表一所示。给定初始步长δ(0)=(δ1,δ2,...,δn)T为1/10的变量区间长度,加速系数γ∈[1,2],收缩系数θ∈[0.1,0.5],以及收敛精度ε=1.0×10-4。
402)优化的目标函数即为步骤三所拟合的二阶多项式f(x),x∈Rn,模式搜索的方向为n维空间的基矢向量e1=(1,0,...,0)T,e2=(0,1,...,0)T,...,en=(0,0,...,1)T。同时,令i=0,取y=x(i),令j=1。
403)做平行于单位矢量ej的探测移动循环;
正向探测移动:若则令j=j+1,否则做负向探测移动;
负向探测移动:若则令j=j+1,否则令y=y,j=j+1;
404)令x(i+1)=y,若f(x(i+1))<f(x(i)),表明模式移动成功,则对x(i+1)沿加速方向做模式移动,令y=x(i+1)+γ(x(i+1)-x(i)),δ(i+1)=θδ(i),i=i+1,进入循环403)。若f(x(i+1))≥f(x(i)),表明模式移动失败,则进入循环405)。
405)判断优化迭代是否收敛,若|δ(i)|<ε,则停止迭代,输出x(i),即为所求设计变量;否则,若x(i+1)=x(i),令y=x(i+1),δ(i+1)=θδ(i),i=i+1,进入循环403),若x(i+1)≠x(i),令y=x(i+1),δ(i+1)=δ(i),i=i+1,进入循环403)。
Claims (6)
1.一种基于中心复合设计的转子结构优化方法,其特征在于,该方法包括:首先确定转子结构优化过程中的设计变量,基于中心复合设计方法生成所需要的N组设计变量,参数化建立相应的有限元模型并计算获得目标函数值,再采用最小二乘法拟合设计变量与目标函数的回归模型,最后利用变步长的模式搜索法对回归模型进行优化,获得回归模型的最优组合。
2.根据权利要求1所述的一种基于中心复合设计的转子结构优化方法,其特征在于,具体包括以下步骤:
第一步,确定转子结构优化的设计变量,确定设计变量的范围与约束条件;
第二步,采用中心复合设计模型生成所需要的N组设计变量,对这N组样本进行参数化有限元计算,得到N组样本所对应的最大等效应力值作为目标函数;
第三步,采用最小二乘法拟合设计变量与目标函数的二阶多项式回归模型;
第四步,对所拟合的回归模型进行优化,给定初始点、初始步长、加速系数、收缩系数、精度,设置约束条件,由表决判断轴向探测移动和模式移动方式,直至满足收敛条件,完成对转子结构的优化。
3.根据权利要求2所述的一种基于中心复合设计的转子结构优化方法,其特征在于,步骤一中,选定其需要优化的结构,通过n个参数来确定其具体的几何形状,即x=(x1,x2,...,xn)为该转子的几何形状参数;同时,为了避免结构内部发生相互干涉,防止在建立几何面等操作时发生错误,确定设计变量的范围与约束条件,即Li≤xi≤Ui,i=1,2,...,n,Li与Ui分别表示变量参数xi上下极限值。
4.根据权利要求3所述的一种基于中心复合设计的转子结构优化方法,其特征在于,步骤二中,采用中心复合设计模型生成N组设计变量,N由中心复合设计模型中参数数量n决定,n由步骤一所确定;同时,对以上N组设计变量所确定的结构实现建模、载荷和边界条件的加载、求解计算和后处理分析过程,提取后处理结果中的最大等效应力值作为目标函数。
5.根据权利要求4所述的一种基于中心复合设计的转子结构优化方法,其特征在于,步骤三中,依据步骤二所得的N组设计变量与目标函数值,选用基于泰勒展开式的二阶多项式模型,对设计变量与目标函数进行多元回归拟合,得到表征最大等效应力值与n维设计变量非线性关系的回归模型;
回归模型形式如:其中f(x)为目标函数,即最大等效应力值,式中包括常数项β0、线性项线性交叉项和二次项βj、βij、βjj分别为线性项、线性交叉项和二次项系数;
采用最小二乘法对以上回归模型进行拟合。
6.根据权利要求5所述的一种基于中心复合设计的转子结构优化方法,其特征在于,步骤四中,具体包括:
401)由原转子结构模型,测定其几何参数以给定设计变量的初始值x(0)=(x1,x2,...,xn),给定初始步长δ(0)=(δ1,δ2,...,δn)T、加速系数γ、收缩系数θ和收敛精度ε;
402)优化的目标函数即为步骤三所拟合的二阶多项式f(x),x∈Rn,模式搜索的坐标方向n维空间的基矢向量e1=(1,0,...,0)T,e2=(0,1,...,0)T,...,en=(0,0,...,1)T,同时,令i=0,取y=x(i),令j=1;
403)做平行于单位矢量ej的探测移动循环:
正向探测移动:若则令否则做负向探测移动;
负向探测移动:若则令否则令y=y,j=j+1;
404)令x(i+1)=y,若f(x(i+1))<f(x(i)),表明模式移动成功,则对x(i+1)沿加速方向做模式移动,令y=x(i+1)+γ(x(i+1)-x(i)),δ(i+1)=θδ(i),i=i+1,进入循环403);若f(x(i+1))≥f(x(i)),表明模式移动失败,则进入循环405);
判断优化迭代是否收敛,若|δ(i)|<ε,则停止迭代,输出x(i),即为所求设计变量;否则,若x(i+1)=x(i),令y=x(i+1),δ(i+1)=θδ(i),i=i+1,进入循环403),若x(i+1)≠x(i),令y=x(i+1),δ(i+1)=δ(i),i=i+1,进入循环403)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811032215.6A CN109194050B (zh) | 2018-09-05 | 2018-09-05 | 一种基于中心复合设计的转子结构优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811032215.6A CN109194050B (zh) | 2018-09-05 | 2018-09-05 | 一种基于中心复合设计的转子结构优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109194050A true CN109194050A (zh) | 2019-01-11 |
CN109194050B CN109194050B (zh) | 2020-10-27 |
Family
ID=64914706
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811032215.6A Active CN109194050B (zh) | 2018-09-05 | 2018-09-05 | 一种基于中心复合设计的转子结构优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109194050B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111478466A (zh) * | 2020-04-14 | 2020-07-31 | 合肥工业大学 | 一种同步磁阻电机转子优化设计方法 |
CN112926129A (zh) * | 2021-04-12 | 2021-06-08 | 交通运输部天津水运工程科学研究所 | 一种应用于港口石化钢制储罐抗风圈优化设计方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103324798A (zh) * | 2013-06-25 | 2013-09-25 | 福州大学 | 基于区间响应面模型的随机模型修正方法 |
KR101600001B1 (ko) * | 2014-09-22 | 2016-03-14 | 주식회사 져스텍 | 모터의 최적 설계 방법 |
CN105846560A (zh) * | 2015-01-30 | 2016-08-10 | 尼得科Sr驱动有限公司 | 用于电机的定子 |
CN105956254A (zh) * | 2016-04-27 | 2016-09-21 | 西安交通大学 | 一种用于转子系统动力学建模的刚体单元法 |
CN108280295A (zh) * | 2018-01-23 | 2018-07-13 | 北京航空航天大学 | 一种基于几何特征驱动的多级轮盘转子结构交互式设计方法 |
-
2018
- 2018-09-05 CN CN201811032215.6A patent/CN109194050B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103324798A (zh) * | 2013-06-25 | 2013-09-25 | 福州大学 | 基于区间响应面模型的随机模型修正方法 |
KR101600001B1 (ko) * | 2014-09-22 | 2016-03-14 | 주식회사 져스텍 | 모터의 최적 설계 방법 |
CN105846560A (zh) * | 2015-01-30 | 2016-08-10 | 尼得科Sr驱动有限公司 | 用于电机的定子 |
CN105956254A (zh) * | 2016-04-27 | 2016-09-21 | 西安交通大学 | 一种用于转子系统动力学建模的刚体单元法 |
CN108280295A (zh) * | 2018-01-23 | 2018-07-13 | 北京航空航天大学 | 一种基于几何特征驱动的多级轮盘转子结构交互式设计方法 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111478466A (zh) * | 2020-04-14 | 2020-07-31 | 合肥工业大学 | 一种同步磁阻电机转子优化设计方法 |
CN112926129A (zh) * | 2021-04-12 | 2021-06-08 | 交通运输部天津水运工程科学研究所 | 一种应用于港口石化钢制储罐抗风圈优化设计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109194050B (zh) | 2020-10-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Mansour | A developed algorithm for simulation of blades to reduce the measurement points and time on coordinate measuring machine (CMM) | |
CN106446343B (zh) | 一种径流叶轮直纹叶片参数化型线的自动提取方法 | |
CN105528503B (zh) | 一种基于结构分解的大型构件动态优化设计方法 | |
CN104008257B (zh) | 一种针对具有复杂曲面的构件的逆向设计方法 | |
CN103970965A (zh) | 燃气涡轮发动机加速寿命试验试车方法 | |
CN101672637A (zh) | 一种复杂曲面的数字化检测方法 | |
CN105446264A (zh) | 基于特征的机床精度优化设计方法 | |
CN104133417A (zh) | 叶片式流体机械数控加工切削力的快速预测方法 | |
CN108334709A (zh) | 基于知识库数据统一管理的透平机械cea集成平台 | |
CN109194050A (zh) | 一种基于中心复合设计的转子结构优化方法 | |
CN106489068B (zh) | 测量值分析装置以及测量值分析方法 | |
CN104153981A (zh) | 一种用于电站循环水泵启停时背压数据的预估方法 | |
CN111504191A (zh) | 一种基于三维激光扫描的航空零件自动快速测量方法 | |
US20200103846A1 (en) | Blade sentencing | |
CN106682292A (zh) | 一种降维模拟退火算法的叶根结构优化方法 | |
CN101839149A (zh) | 叶轮机械叶片回转截面型线计算方法 | |
CN108536932B (zh) | 基于互扭约束条件下的航空叶片积叠轴垂直度计算方法 | |
CN108628956B (zh) | 产品创新设计的方法及系统 | |
CN108574290B (zh) | 受迫振荡的振荡源定位方法、装置、终端及可读存储介质 | |
Derakhshan et al. | Optimal design of axial hydro turbine for micro hydropower plants | |
Safari et al. | A novel combination of adaptive tools for turbomachinery airfoil shape optimization using a real-coded genetic algorithm | |
CN104180789A (zh) | 基于图形匹配算法的叶片检测方法 | |
CN105005210B (zh) | 机电一体化仿真系统及使用其的方法 | |
Zhang et al. | The dynamic position compensation method of the online laser detection system for the ring workpiece | |
CN106339543A (zh) | 水轮机过流部件形貌修复模型构建方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |