主权项 |
一种高低压压气机过渡流道优化设计方法,其特征在于包括以下步骤:根据给定的一维设计初始约束条件将压气机过渡流道的机匣及轮毂型线分别用两条带有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>Σ</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>Σ</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>ΔV</mi><mi>θ</mi></msub><mrow><mn>2</mn><msub><mi>σ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>θ</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>ω</mi><mo>*</mo></msup><mo>=</mo><mfrac><mrow><mn>2</mn><mi>σ</mi></mrow><msubsup><mrow><mi>cos</mi><mi>β</mi></mrow><mn>2</mn><mo>*</mo></msubsup></mfrac><msup><mrow><mo>(</mo><mfrac><mi>θ</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><</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是过渡流道轴向坐标;采用带约束的小生境遗传算法优化求得新的型面方程,以给定进化代数为优化结束标志,达到指定进化代数则优化结束,否则生成新的控制点坐标,重复自生成新贝塞尔曲线开始的优化过程。 |