一种可提高核主泵惰转时间的叶轮的优化设计方法
【专利摘要】本发明涉及一种可提高核主泵惰转时间的叶轮的优化设计方法，首先建立可以通过优化叶轮主要参数来提高惰转工况运转效率的主目标函数，然后将主要参数分为三组分别作为设计变量，在初始约束条件内，对这三个分目标函数分别利用复合形法对惰转效率优化，分别得到使效率最高的叶轮主要几何参数的最优点，然后在最优点附近小范围确定一个新的约束，在缩小的新约束内用可行方向法对主目标函数优化求解，得到最终优化结果。本发明对核主泵叶轮的设计，能够改善叶轮惰转运行效率偏低问题，延长惰转时间，提高核安全性。同时，本方法创新运用分目标缩小约束范围和主目标寻找最优点的分步优化设计方法，提高了优化的速度，保证了优化结果有效性。
【专利说明】
-种可提高核主累惰转时间的叶轮的优化设计方法
技术领域
[0001] 本发明所设及的是核主累惰转效率分析及叶轮机械优化设计领域，具体设及采用 现代机械优化设计方法对累叶轮参数进行分步优化的设计方法。
【背景技术】
[0002] 核电作为重要的清洁能源，已成为多数发达国家的主要电力来源。核主累全称核 反应堆冷却剂主循环累，是核电站的"屯、脏"，其功能是驱动核岛内高放射性高溫高压水循 环，将反应堆忍核裂变的热能传递给蒸汽发生器产生蒸汽，推动汽轮机发电。在压水堆核电 站中，一回路循环的工质水与核反应直接接触，吸收和传递大量的热量，使工质工作溫度、 工作压力很高，工作流量很大，还具有放射性，因此与之匹配的驱动冷却剂工质循环的核主 累具有高溫高压大流量的特征，运一特征就是核主累的设计难点之一。另一方面，由于核反 应堆的热量是由核主累驱动一回路工质水循环带走，所W核主累是核电站中最重要的累， 也是核岛内唯一高速运转的装备，无论是对安全性还是对性能，它的要求都比普通水累高 很多。
[0003] 核主累若频繁的停机会给电厂造成巨大的经济损失，核主累的突发性故障甚至可 能带来不可估量的灾难，因此，对核主累安全系数的要求非常高，一般要求其能够安全无故 障运行年。在全厂断电情况下，核主累依靠惯性继续惰转，保持一回路一定的流量，维持堆 忍冷却，保证反应堆不会达到泡核沸腾状态，从而保障核安全。因此在优化设计时，要求核 主累在断电惰转的工况下运转效率高，从而延长惰转时间，尽量降低核主累因失去外电源 时的流量突变，W确保全厂断电后堆忍安全。而在核主累中，叶轮是最重要的过流部件，叶 轮设计的合理与否直接影响累的各项性能，因此对叶轮进行优化对提高惰转时间有一定的 指导意义。
[0004] 经检索，与本发明相关的专利有：一种核主累叶轮设计方法（公开号： CN102691671A)，设及一种带核主累叶轮的设计方法，其设计优化叶轮的叶片厚度和包角， 但是仅从叶轮叶片厚度进行优化，可能会带来其他惰转运行问题;一种失水事故下核主累 极大效率水力设计方法(CN104595232A)，提供了一种水力设计方法确定叶轮外径、出口宽 度、出口安放角、叶片数和喉部面积，使核主累在失水事故下设计的额定流量点附近有最大 效率点，提高了核主累在失水事故下的效率，同时提高了核主累的抗气蚀性能和运行可靠 性，但此方法依靠设计者的经验，缺乏经验的情况下结果会存在很大偏差。
【发明内容】
[0005] 针对W上存在问题，本发明为核主累惰转优化设计提供一种更加简单、系统的方 法和思路，W优化叶轮的性能参数为目标提高核主累惰转时的运转效率，获得较长的惰转 时间，保证核主累的断电情况下安全高效地运行。
[0006] 本发明的技术方案是：
[0007] -种可提高核主累惰转时间的叶轮的优化设计方法，具体包括如下步骤：
[000引 SI:首先寻找一个可W通过优化叶轮主要参数:叶轮进口直径Do、叶轮出口直径化、 叶轮出口宽度b2、叶片进口角01、叶片出口安放角阮、叶片包角巧、叶片数Z，来提高惰转工况 下运转效率的主目标函数；
[0009] S2:然后将主要几何参数分为=组设计变量，并建立=个分目标函数，再在初始约 束条件内利用复合形法对惰转效率分目标函数优化，分别得到=个分目标函数的最优解，
[0010] S3:最后在最优解各叶轮参数附近小范围确定一个新的约束条件，在新约束条件 内用可行方向法对主目标函数优化求解，得到最终优化结果。
[0011] 步骤Sl中，所述的设计变量为：
[001^ 主目标函数为：
[0013]
[0014] 式中，Do-叶轮进口直径，单位米；
[001引 D2-叶轮出日直径，单位米；
[0016] b2-叶轮进口宽度，单位米；
[0017]执-叶轮进日角，单位胺
[001引 叶轮出口安放角，单位度；
[0019] <?^-叶片包角，单位度；
[0020] Z-叶片数。
[0021] 步骤S2中，所述的=组分设计变量为：
[0022] Xi=[Dg、D2]t;
[0023] X2=[b2、01]T;
[0024]
[0025]
[0026]
[0027]
[002引
[0029] 步骤S2中，所述初始约束条件内的相关设计变量先参照如下设计方法计算：
[0030] 1)叶轮入口直径Do在兼顾效率和汽蚀的情况下满足，
[0031]
[0032] 式中，Q-流量，单位立方米每秒;n-转速，单位立转每分；
[0033] 2)叶轮出口直径化在考虑修正系数心时的约束条件满足，
[0034]
[0035] 式中，ns-比转速；
[0036] 3)叶轮出口宽度b2由经验及相应推导公式得，其约束条件满足，
[0037]
b2取中间值；
[0038] 4)叶片进口角执由经验公式及累的参数得，其约束条件满足，
[0039]
[0040] 式中，累的容积效率；
[0041] F-累的过水断面积，单位平方米；
[0042] (6-累的叶片排挤系数；
[0043] D-累的进口处计算点直径，单位米；
[0044] m-经验系数；
[0045] A执-累的进口冲角；
[0046] 5)叶片出口角&由经验得，其约束条件满足，22°
[0047] 6)确定叶片包角巧的通常取值为90 <v>若J20,取梦= 100;
[004引 7)叶片数的选择为4《2《7;
[0049] 上述屯项计算完成后，所述的初始约束条件得到为：
[0050] 分目标函数f (Xi)的初始约束条件为[0.9Do，1. IDo]、[0.9D2，1. 1化];
[0051] 分目标函数f(X2)的初始约束条件为[0.抓2，1.化2]、[0.9执，1.1&];
[0化2] 分目标函数f(X3)的初始约束条件为[22°，30°]、[0,辦)，：1.1础[4,7]。
[0053] 步骤S2中，所述的运用复合形法对惰转效率进行优化，即对分目标函数进行优化 求解，具体包括如下步骤：
[0054] (1)选择复合形顶点数k，一般取n+l《k《化，在可行域内构成只有k个顶点的初始 复合形；
[0055] (2)计算复合形k个顶点的目标函数值，选出其中最大值，即 [0化6]最坏点X 邮，F(X 化)）=max{F(X(j))J = l，2r..，k};
[0化7] 好次点乂似，尸村似）=111曰又巧村〇'）^' = 1，2，...，4;但_]'声田；
[0化引最好点 X化)，F(X化)）=min{F(X(j))，j = l，2，...k};
[0059] ( 3 )计算除最坏点X ( H )外其余k - 1个顶点的中屯、点X ( S )，即
检验中屯、点是否在可行域内。如果在可行域 内，则继续执行第(4)步，否则转到第(5)步；
[0060] (4)若X(s>点在可行域内，则X(H哺X(s>的连线方向上取映射点X ("，X("=xW+a(X W-X?)，式中，a-映射系数，一般取a = 1.3;
[0061 ]若XW越出了可行域，则需将其退回，即将映射系数a减半，重新计算XW ;如果还不 满足可行性，继续将a减半，直到映射点成为可行点为止；
[0062] (5)若点不在可行域内，此时可行域为非凸集;按上述第(4)步计算的映射点不 可能是可行点，此时利用中屯、点X(s哺最好点X(U重新确立一个区间，在此区间内重新随机 产生k个顶点构成复合形，新的区间其边界值为，
[0063] 若Xi(L)<xi(s) (i = I，2，…，n)，则巧
[0064] 若xi(L)〉xi(s) (i = 1，2，…，n)，则巧
[00化]重新构成复合形再重复第(2)、（3)步，直至成为可行点为止。
[0066] (6)计算映射点目标函数值FU(K));
[0067] 若f(xW)<f(x?)，则用映射点代替最坏点xW，构成新复合形，完成一次迭代 计算，转向第(2)步，否则计息下一步。
[006引（7)若F(XW)^(XW)，则将映射系数a减半，重新计算映射点；
[0069] 如果新的映射点X("既为可行点，又满足F(XWkfuW),则WX("代替X?，完成本 次迭代;否则继续将a减半，直到当a值小于预先给定的一个很小数U例如C = ICT5)时;若目 标函数仍无改进，则转向第(4)步，但此时改用次点来代替前次的最坏点X?进行映射；
[0070] (8)终止准则：反复执行上述迭代过程，复合形逐渐变小且向最优点逼近，直到满 足
时迭代计算可W结束;此时复合形中目标函数值最小的顶点 即为最优解，其中，为复合形所有顶点的点集中屯、，巧
[0071] 在分目标函数求解得到最优解 Xl*=[D0(l)*、D2(l)*]T、X2*=[b2(2)*、ei(2)*]T、
[0072] 步骤S3中，所述的W最优解的各叶轮参数为中屯、，在小范围内确定新的约束条
[0073] 步骤S3中，主目标函数在新的约束条件下运用ZoutendUk方法改进后的Topkis- Veinott可行方向法优化求解，具体步骤如下：
[0074] (1)给定初始点XW、允许误差e，确定初始可行方向SW，设置k = l;
[00对 （2)确定标号集E(xW) = {j|g^xW) = 0，l《j《m}是否为空集，若为空集转向第 (3)步，否则转向第(4)步，其中m为约束条件数；
[0076] (3)判断||V/(X…)1^，若成立则直接得到最优解X* = X…并输出结果，否则令
索后转向第(5)步；
[0077] (4)求解W下线型规戈
，并判断Zk = O是否成 立，若成立则直接得到最优解x*=x?并输出结果，否则转向第巧)步；
[007引（5)求解一维捜索问题. + /j' \得到Ak，其中
[义to < 乂 < A胃
[0079] 入。3义=3啡{入|的（乂化)+入5化）《〇，11=1，2，...，111};
[0080]
然后转向第(2)步；
[0081] 在新的约束条件下对主目标函数进行上述优化设计后，对求解的最优解进行圆整 验证后，可W得到满足一种可提高核主累惰转时间的叶轮的优化设计方法的最终优化结果
[0082] 本发明的有益效果在于：
[0083] (1)本发明提出的一种可提高核主累惰转时间的叶轮的优化设计方法，对核主累 叶轮的设计制造具有一定指导意义，能够改善叶轮运行过程中效率偏低的问题，延长核主 累全厂断电状况下的惰转时间，提高了核主累的核安全性。
[0084] (2)本方法创新运用了分目标缩小约束范围和主目标寻找最优点的分步优化的设 计方法，不仅提高了优化的速度，同时保证了优化结果的有效性。
【附图说明】
[0085] 图1是复合形法程序框图。
[0086] 图2是可行方向法程序框图。
【具体实施方式】
[0087] 下面结合附图和具体技术方案对本发明进一步说明。
[0088] -种可提高核主累惰转时间的叶轮的优化设计方法，建立通过优化叶轮主要参数 (叶轮进口直径Do、叶轮出口直径〇2、叶轮出口宽度b2、叶片进口角01、叶片出口安放角02、叶 片包角譽、叶片数Z)来提高惰转工况下运转效率的主目标函数，具体过程为：
[0089] 设计变量为:文=[D。、D:、b,、知/，':、巧、到T.，主目标函数为：
[0090]
[0091 ] 式中，Do-叶轮进口直径，单位米；
[0092] D2-叶轮出口直径，单位米；
[0093] b2-叶轮进口宽度，单位米；
[0094] 执-叶轮进口角，单位度；
[0095] 叶轮出口安放角，单位度；
[0096] 巧-叶片包角，单位度；
[0097] Z-叶片数；
[0098] 然后将主要设计变量分为=组分设计变量为：
[0099] Xi= [Dg、D2]t;
[0100] X2=[b2、&]T;
[0101]
[0103]
[0102] 构建的S个分目标函数为：
[0104]
[0105]
[0106] 需要确定求解优化过程中必须满足的最初约束条件，其具体实施过程中必须结合 相关的设计参数和经验：
[0107] 1)叶轮入口直径Do在兼顾效率和汽蚀的情况下满足，
[010 引
[0109] 式中，Q-流量，单位立方米每秒；
[0110] n-转速，单位立转每分；
[0111] 2)叶轮出口直径化在考虑修正系数Ka,时的约束条件满足，
[0112]
[0113] 式中，ns-比转速；
[0114] 3)叶轮出口宽度b2由经验及相应推导公式得，其约束条件满足，
[0115]
，b2取中间值；
[0116] 4)叶片进口角执由经验公式及累的参数得，其约束条件满足，
[0117]
[011引式中，累的容积效率；
[0119] F-累的过水断面积，单位平方米；
[0120] (6-累的叶片排挤系数；
[0121] D-累的进口处计算点直径，单位米；
[0122] m-经验系数；
[0123] A执-累的进口冲角；
[0124] 5)叶片出口角&由经验得，其约束条件满足，22°
[012引6)确定叶片包角巧的通常取值为9炉《巧含120°戚巧=100。，
[0126] 7)叶片数的选择为4《2《7;
[0127] 上述屯项计算完成后，得到的所述的初始约束条件为：
[01測分目标函数f (Xi)的初始约束条件为[0.9D0，1. IDo]、[0.9D2，1. 1化];
[0129] 分目标函数f (X2)的初始约束条件为[0.抓2，1.化2]、[0.9执，1.101];
[0130] 分目标函数f (X3)的初始约束条件为[22。，30。]、[0.9;。J .岭]、|[4，7]。
[0131] 在初始约束条件内对分目标函数运用复合形法进行求解计算，其中复合形法就是 在n维设计空间的可行域内，对复合形（即n维空间内由n+l^k^2n个顶点所构成的多面体） 各顶点的目标函数逐一进行比较，不断地去掉最坏点(对于求极小值问题，即目标函数的最 大点），代之W既能使目标函数值有所下降，又满足所有约束条件的新点，逐步调向最优点。
[0132] 如附图1所示的复合形法程序框图，具体步骤如下：
[0133] (1)选择复合形顶点数k，一般取n+l《k《化，在可行域内构成只有k个顶点的初始 复合形。
[0134] (2)计算复合形k个顶点的目标函数值，选出其中最大值，即 [013 引最坏点X 邮，F(X 化)）=max{F(X(j))J = l，2r..，k};
[0136] 好次点乂似，尸村("）=111曰义巧村〇'）^' = 1，2，...，4;但_]'声田；
[0137] 最好点 X(l)，F(X 化)）=min{F(X(j))J = l，2，.-k}。
[013引（3 )计算除最坏点X ( H >外其余k - 1个顶点的中屯、点X ( S )，即
，检验中屯、点是否在可行域内。如果在可行域 内，则继续执行第(4)步，否则转到第(5)步。
[0139] (4)若X(s>点在可行域内，则X(H哺X (s>的连线方向上取映射点X("，X("=xW+a(X (引-X邮），
[0140] 式中a-映射系数，一般取a = l.3。
[0141] 若XW越出了可行域，则需将其退回，即将映射系数a减半，重新计算XW;如果还不 满足可行性，继续将a减半，直到映射点成为可行点为止。
[0142] (5)若点不在可行域内，此时可行域为非凸集。按上述第(4)步计算的映射点不 可能是可行点，此时利用中屯、点X(s哺最好点X(U重新确立一个区间，在此区间内重新随机 产生k个顶点构成复合形。新的区间其边界值为，
[0143] 若Xi(L)<xi(s) (i = 1，2，...，n)，则耳
[0144] 若Xi(L)〉xi(s) (i = 1，2，...，n)，则耳
[0145] 重新构成复合形再重复第(2)、（3)步，直至X?成为可行点为止。
[0146] (6)计算映射点目标函数值F(XW)。若F(XWkfuW)，则用映射点X("代替最坏点 X?，构成新复合形，完成一次迭代计算，转向第(2)步，否则计息下一步。
[0147] (7)若F(XW)WX^),则将映射系数a减半，重新计算映射点。如果新的映射点 既为可行点，又满足F(xW)<F(xW)，则WxW代替X?，完成本次迭代。否则继续将a减半，直 到当a值小于预先给定的一个很小数C(例如C = ICT5)时。若目标函数仍无改进，则转向第(4) 步，但此时改用次点来代替前次的最坏点X?进行映射。
[0148] (8)终止准则：反复执行上述迭代过程，复合形逐渐变小且向最优点逼近，直到满 足
时迭代计算可W结束。此时复合形中目标函数值最小的顶 点即为最优解，其中，为复合形所有顶点的点集中屯、，即
[0149] 分目标函数求解得到最优解为：
[0150] Xl*= [Do(i)*、D2 山叩、拉 *=化 2(2)*、&(2)叩
[0151] 我们W最优解各导叶参数为中屯、，在小范围内确定新的约束条件如下：
[0152] 如附图2所示的可行方向法程序框图，具体步骤如下：
[0153] (1)给定初始点X(i>、允许误差￡，确定初始可行方向S(i>，设置k = l。
[0154] (2)确定标号集E(X(k)) = {j|gj(X(k))=〇，l《j《m}是否为空集，
[0155] 若为空集转向第(3)步，否则转向第(4)步，其中m为约束条件数。
[0156] (3)判断||W(^")|^€，若成立则直接得到最优解X* = XW并输出结果，否则令 沪"=-Y化护"），然后转向第(5)步。
[0157] (4)求解W下线型规划
E巧义">)，并判断Zk = O是否成立， 若
[015引成立则直接得到最优解X* = X?并输出结果，否则转向第巧)步。
[0159] (5)求解一维捜索问是
，得到入k，其中
[0160] Ajiax = sup{A I 扣（X(k)+入S(k)《〇 ,U = 1,2,... ,m}。
[0161] (6)令
，然后转向第(2)步。
[0162] 在新的约束条件下对主目标函数进行上述优化设计后，对求解的最优解进行圆整 验证后，可W得到满足一种可提高核主累惰转时间的叶轮的优化设计方法的最终优化结果
[0163] 本发明不限于上述实施例，也包含本发明构思范围内其它实施例和变形例。
【主权项】
1. 一种可提高核主栗惰转时间的叶轮的优化设计方法，其特征在于，包括如下步骤： Sl:首先寻找一个可以通过优化叶轮主要参数:叶轮进口直径D〇、叶轮出口直径D2、叶轮 出口宽度b2、叶片进口角扮、叶片出口安放角β2、叶片包角辦、叶片数Z，来提高惰转工况下运 转效率的主目标函数； S2:然后将主要几何参数分为三组设计变量，并建立三个分目标函数，再在初始约束条 件内利用复合形法对惰转效率分目标函数优化，分别得到三个分目标函数的最优解， S3:最后在最优解各叶轮参数附近小范围确定一个新的约束条件，在新约束条件内用 可行方向法对主目标函数优化求解，得到最终优化结果。2. 根据权利要求1所述的一种可提高核主栗惰转时间的叶轮的优化设计方法，其特征 在于，步骤Sl中，所述的设计变量为：Y = [D(,、A、b2、/?、爲、A 主目标函数为：式中，Dq-叶轮进口直径，单位米； D2-叶轮出口直径，单位米； b2_叶轮进口宽度，单位米； -叶轮进口角，单位度； 叶轮出口安放角，单位度； 史-叶片包角，单位度； Z-叶片数。3. 根据权利要求1所述的一种可提高核主栗惰转时间的叶轮的优化设计方法，其特征 在于，步骤S2中，所述的三组分设计变量为：4. 根据权利要求1所述的一种可提高核主栗惰转时间的叶轮的优化设计方法，其特征 在于，步骤S2中，所述的初始约束条件为： 分目标函数f (X1)的初始约束条件为[0.9D〇，I. IDo]、[0.9D2，I. ID2]; 分目标函数f (X2)的初始约束条件为[0.9b2，I. Ib2]、[0.9?，1. 1?]; 分目标函数f (X3)的初始约束条件为5. 根据权利要求4所述的一种可提高核主栗惰转时间的叶轮的优化设计方法，其特征 在于，步骤S2中，所述初始约束条件内的相关设计变量先参照如下设计方法计算： 1) 叶轮入口直径Do在兼顾效率和汽蚀的情况下满足，式中，Q-流量，单位立方米每秒;η-转速，单位立转每分； 2) 叶轮出口直径02在考虑修正系数Ia时的约束条件满足，式中，ns-比转速； 3) 叶轮出口宽度b2由经验及相应推导公式得，其约束条件满足，4) 叶片进口角扮由经验公式及栗的参数得，其约束条件满足，式中，栗的容积效率； F-栗的过水断面积，单位平方米； Φ-栗的叶片排挤系数； D-栗的进口处计算点直径，单位米； m-经验系数； A βι_栗的进口冲角； 5) 叶片出口角β2由经验得，其约束条件满足，22°彡β2彡30° ; 6) 确定叶片包角P的通常取值为90 120，.取P= 100 ; 7) 叶片数的选择为4<Ζ<7。6.根据权利要求1所述的一种可提高核主栗惰转时间的叶轮的优化设计方法，其特征 在于，步骤S2中，所述的运用复合形法对惰转效率进行优化，即对分目标函数进行优化求 解，具体包括如下步骤： (1) 选择复合形顶点数k，一般取η+1彡k彡2η，在可行域内构成只有k个顶点的初始复合 形； (2) 计算复合形k个顶点的目标函数值，选出其中最大值，BP 最坏点 X(H)，F(X(H))=maX{F(X(j))，j = l，2r"，k}; 好次点 X(G)，F(X(G))=maX{F(X(j)，j = l，2，."，kdSj#H}; 最好点 X(L)，F(X(L))=min{F(X(j))，j = l，2，...k}; ⑶计算除最坏点X(H)外其余k-1个顶点的中心点X(S)，B但/WH)，检验中心点X(s)是否在可行域内。如果在可行域内，则继续执行第(4)步，否则转到 第(5)步； (4)若X(s)点在可行域内，则X(H)和X(s)的连线方向上取映射点X (R)，X(R)=X(s)+a(X(s)-X (11))，式中，〇-映射系数，一般取€[ = 1.3; 若X(R)越出了可行域，则需将其退回，即将映射系数α减半，重新计算X(R);如果还不满足 可行性，继续将α减半，直到映射点X(R)成为可行点为止； (5) 若X(s)点不在可行域内，此时可行域为非凸集;按上述第(4)步计算的映射点不可能 是可行点，此时利用中心点X (S)和最好点Χα)重新确立一个区间，在此区间内重新随机产生k 个顶点构成复合形，新的区间其边界值为，重新构成复合形再重复第(2)、（3)步，直至X(s)成为可行点为止。 (6) 计算映射点目标函数值F(Xw); 若F(X(R))<F(X(H))，则用映射点X(R)代替最坏点X (H)，构成新复合形，完成一次迭代计 算，转向第(2)步，否则计息下一步。 (7) 若F(X(R)) >F(X(H))，则将映射系数α减半，重新计算映射点； 如果新的映射点X(R)既为可行点，又满足F(X(R))<F(X(H))，则以X (R)代替Χ(Η)，完成本次 迭代;否则继续将α减半，直到当α值小于预先给定的一个很小数ξ(例如ξ = 10-5)时;若目标 函数仍无改进，则转向第(4)步，但此时改用次点X(e)来代替前次的最坏点Χ (Η)进行映射； (8) 终止准则：反复执行上述迭代过程，复合形逐渐变小且向最优点逼近，直到满足时迭代计算可以结束;此时复合形中目标函数值最小的顶 点即为最优解，其中，X(e)为复合形所有顶点的点集中心，B丨 ?.尸在分目标函数求解得到最优_7. 根据权利要求1所述的一种可提高核主栗惰转时间的叶轮的优化设计方法，其特征 在于，步骤S3中，所述的以最优解的各叶轮参数为中心，在小范围内确定新的约束条件如 下.8. 根据权利要求1所述的一种可提高核主栗惰转时间的叶轮的优化设计方法，其特征 在于，步骤S3中，主目标函数在新的约束条件下运用Zoutendi jk方法改进后的TopkiS-Veinott可行方向法优化求解，具体步骤如下： (1) 给定初始点X(1)、允许误差ε，确定初始可行方向S(1)，设置k= 1; (2) 确定标号集E(X(k)) = {j|gj(X(k))=0，l<j<m}是否为空集，若为空集转向第⑶步， 否则转向第(4)步，其中m为约束条件数； (3) 判断|v/(xw)|以，若成立则直接得到最优解X* = X(k)并输出结果，否则令瘩后转向第(5)步； (4) 求解以下线型规j，并判断21{ = 0是否成立，若 成立则直接得到最优解X*=X(k)并输出结果，否则转向第(5)步；在新的约束条件下对主目标函数进行上述优化设计后，对求解的最优解进行圆整验证 后，可以得到满足一种可提高核主栗惰转时间的叶轮的优化设计方法的最终优化结果
【文档编号】F04D29/18GK106015082SQ201610318493
【公开日】2016年10月12日
【申请日】2016年5月13日
【发明人】朱荣生, 蔡峥, 王秀礼, 卢永刚, 付强, 康俊鋆
【申请人】江苏大学