发明名称 高低压压气机过渡流道优化设计方法
摘要 本发明公开了一种高低压压气机过渡流道优化设计方法,用于解决现有压气机过渡流道设计方法所设计的压气机过渡流道流动损失大的技术问题。技术方案是通过现有设计方法初步确定压气机过渡流道的几何形状,再采用贝塞尔曲线构建满足初始设计结果的端壁型线方程,并在此基础上进一步叠加所构造好的支板。继而随机生成机匣和轮毂的几何参数,并生成新的型线。然后用通流求解流场性能参数。最后,以支板出口处的总压损失系数为优化目标函数,并以沿轮毂的静压恢复系数梯度小于优化前压力恢复系数梯度为约束条件,采用小生境遗传算法优化直至进化到给定代数为止。由于采用贝塞尔曲线构建压气机过渡流道的型线方程,所设计的压气机过渡流道流动损失减小。
申请公布号 CN104834768A 申请公布日期 2015.08.12
申请号 CN201510174689.4 申请日期 2015.04.13
申请人 西北工业大学 发明人 吴虎;侯朝山;杨金广;刘昭威
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 西北工业大学专利中心 61204 代理人 王鲜凯
主权项 一种高低压压气机过渡流道优化设计方法,其特征在于包括以下步骤:根据给定的一维设计初始约束条件将压气机过渡流道的机匣及轮毂型线分别用两条带有5个控制点的4阶贝塞尔曲线来构建,每条曲线的参数化方程为<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><mi>x</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mn>4</mn></msubsup><msubsup><mi>C</mi><mn>4</mn><mi>i</mi></msubsup><msup><mi>t</mi><mi>i</mi></msup><msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>t</mi><mo>)</mo></mrow><mrow><mn>4</mn><mo>-</mo><mi>i</mi></mrow></msup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><msub><mi>x</mi><mi>i</mi></msub></mtd></mtr><mtr><mtd><mi>y</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mn>4</mn></msubsup><msubsup><mi>C</mi><mn>4</mn><mi>i</mi></msubsup><msup><mi>t</mi><mi>i</mi></msup><msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>t</mi><mo>)</mo></mrow><mrow><mn>4</mn><mo>-</mo><mi>i</mi></mrow></msup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><msub><mi>y</mi><mi>i</mi></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000697691270000011.GIF" wi="1452" he="167" /></maths>式中,<img file="FDA0000697691270000012.GIF" wi="78" he="98" />是二项式系数,<img file="FDA0000697691270000013.GIF" wi="524" he="97" />是Bernstein多项式;m是多项式的阶数;t是曲线的控制参数,其取值范围为0到1;x<sub>i</sub>是控制点横坐标;y<sub>i</sub>是控制点纵坐标;将已设计好的支板与上述新生成的压气机过渡流道相叠加;在适当范围内随机变化控制点坐标,利用贝塞尔曲线生成新的压气机过渡流道几何型线;对定常、无黏、轴对称的流动的径向平衡方程积分可求得沿计算站的子午速度分布;公式(2)为径向平衡方程的线性常微分形式;方程中式A(l)如公式(3)所示,B(l)如公式(4)所示;<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mfrac><msubsup><mi>dc</mi><mi>m</mi><mn>2</mn></msubsup><mi>dl</mi></mfrac><mo>+</mo><mi>A</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow><msubsup><mi>c</mi><mi>m</mi><mn>2</mn></msubsup><mo>=</mo><mi>B</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000697691270000014.GIF" wi="1236" he="127" /></maths><img file="FDA0000697691270000015.GIF" wi="1233" he="124" /><img file="FDA0000697691270000016.GIF" wi="1896" he="277" />式中,c<sub>m</sub>是子午速度;r<sub>c</sub>是流线曲率半径;<img file="FDA0000697691270000017.GIF" wi="46" he="70" />是流线倾斜角;r是计算站径向夹角;m是子午方向;ε是叶片倾斜角;W<sub>θ</sub>是相对切向速度;l是计算站方向;Ω是转子旋转角速度;I是转焓;s是熵;F<sub>θ</sub>是叶片切向力;β是相对流动角;T是静温;过渡流道中支板的流动损失的求解方程为<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msup><mi>D</mi><mo>*</mo></msup><mo>=</mo><mn>1</mn><mo>-</mo><mfrac><msub><mi>V</mi><mn>2</mn></msub><msub><mi>V</mi><mn>1</mn></msub></mfrac><mo>+</mo><mfrac><msub><mi>&Delta;V</mi><mi>&theta;</mi></msub><mrow><mn>2</mn><msub><mi>&sigma;V</mi><mn>1</mn></msub></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000697691270000018.GIF" wi="1212" he="126" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msup><mrow><mo>(</mo><mfrac><mi>&theta;</mi><mi>c</mi></mfrac><mo>)</mo></mrow><mo>*</mo></msup><mo>=</mo><mn>0.003</mn><mo>+</mo><msup><mrow><mn>0.02375</mn><mi>D</mi></mrow><mo>*</mo></msup><mo>-</mo><msup><mrow><mn>0.05</mn><mi>D</mi></mrow><mrow><mo>*</mo><mn>2</mn></mrow></msup><mo>+</mo><msup><mrow><mn>0.123</mn><mi>D</mi></mrow><mrow><mo>*</mo><mn>3</mn></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000697691270000019.GIF" wi="1539" he="121" /></maths><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msup><mi>&omega;</mi><mo>*</mo></msup><mo>=</mo><mfrac><mrow><mn>2</mn><mi>&sigma;</mi></mrow><msubsup><mrow><mi>cos</mi><mi>&beta;</mi></mrow><mn>2</mn><mo>*</mo></msubsup></mfrac><msup><mrow><mo>(</mo><mfrac><mi>&theta;</mi><mi>c</mi></mfrac><mo>)</mo></mrow><mo>*</mo></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000697691270000021.GIF" wi="1218" he="135" /></maths>式中,D<sup>*</sup>是扩散因子;V<sub>1</sub>是进口绝对速度;V<sub>2</sub>是出口绝对速度;ΔV<sub>θ</sub>是切向速度增量;σ是稠度;<img file="FDA0000697691270000022.GIF" wi="107" he="135" />是设计状态下总压损失参数;ω<sup>*</sup>是设计状态下支板出口处总压损失系数;β<sub>2</sub>是出口相对气流角;优化目标函数为min obj=ω             (9)式中,ω是支板出口处的总压损失系数;约束条件为<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mrow><mo>|</mo><mfrac><mi>dcp</mi><mi>dx</mi></mfrac><mo>|</mo></mrow><mi>max</mi></msub><mo>&lt;</mo><msub><mrow><mo>|</mo><mfrac><mi>dcp</mi><mi>dx</mi></mfrac><mo>|</mo></mrow><mrow><mi>ariginal</mi><mo>,</mo><mi>max</mi></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000697691270000023.GIF" wi="1232" he="137" /></maths>式中,cp是沿轮毂的静压恢复系数;x是过渡流道轴向坐标;采用带约束的小生境遗传算法优化求得新的型面方程,以给定进化代数为优化结束标志,达到指定进化代数则优化结束,否则生成新的控制点坐标,重复自生成新贝塞尔曲线开始的优化过程。
地址 710072 陕西省西安市友谊西路127号