发明名称 可变进口导流叶片的设计方法及叶片、压气机
摘要 本发明涉及压气机技术领域,尤其涉及一种可变进口导流叶片的设计方法,包括S1:确定进口来流马赫数以及进口气流角和出口气流角,S2:确定无量纲弯角分布函数;S3:选定叶型设计点迎角,确定叶型的弯角;S4:确定叶型的最大相对厚度,前缘及尾缘厚度;S5:确定叶型的最大相对厚度位置;S6:确定厚度分布函数;S7:对称积叠获得叶身;S8:给叶身部分加上前缘和尾缘;S9:对S8获得的叶型进行数值仿真,计算其性能并判断是否合适。本发明还提供的设计方法,拓展可用进口迎角范围,并在全工况范围内降低平均总压损失系数。本发明还涉及一种可变进口导流叶片的设计方法,以及利用该方法制成的叶片。
申请公布号 CN106570213A 申请公布日期 2017.04.19
申请号 CN201610887704.4 申请日期 2016.10.11
申请人 北京航空航天大学 发明人 刘宝杰;施恒涛;于贤君
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 北京路浩知识产权代理有限公司 11002 代理人 汤财宝
主权项 一种可变进口导流叶片的设计方法,其特征在于,包括以下步骤:S1:确定进口导流叶片不同展向位置处设计点的进口来流马赫数Ma<sub>1</sub>以及设计进口气流角β<sub>1</sub>和设计出口气流角β<sub>2</sub>;S2:确定无量纲弯角分布函数φ(t),首先定义叶型无量纲坐标系,以前缘点为坐标原点,前缘点和尾缘点之间的连线为X轴,X轴表示无量纲弦长位置,t=x/C,其中C是叶型的弦长,以垂直于前缘点和尾缘点之间的连线所在的水平面的直线为Y轴,以无量纲弦长位置t表示的中弧线弯角函数:f<sub>1</sub>(t)=φ(t)·θ+θ<sub>0</sub> (0≤t≤1)式中,θ为叶型的中弧线弯角,θ<sub>0</sub>为叶型中弧线的前缘构造角,其中,θ<sub>0</sub>的正负号规定为当y/C&gt;0为负,y/C&lt;0为正,函数φ(t)是以无量纲弦长位置t表示的待定函数,函数φ(t)需要满足φ(t)=(C<sub>1</sub>t<sup>4</sup>+C<sub>2</sub>t<sup>3</sup>+C<sub>3</sub>t<sup>2</sup>+C<sub>4</sub>t) (0≤t≤1),其中,C<sub>1</sub>、C<sub>2</sub>、C<sub>3</sub>、C<sub>4</sub>需要满足C<sub>1</sub>+C<sub>2</sub>+C<sub>3</sub>+C<sub>4</sub>=1<maths num="0001"><math><![CDATA[<mrow><mfrac><mrow><msup><mi>d</mi><mn>2</mn></msup><msub><mi>f</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mrow><msup><mi>dt</mi><mn>2</mn></msup></mrow></mfrac><mo>&lt;</mo><mn>0</mn></mrow>]]></math><img file="FDA0001128856230000011.GIF" wi="292" he="143" /></maths>定义参数C<sub>FL</sub>:<maths num="0002"><math><![CDATA[<mrow><msub><mi>C</mi><mrow><mi>F</mi><mi>L</mi></mrow></msub><mo>=</mo><mfrac><mrow><mo>(</mo><msub><mi>f</mi><mn>1</mn></msub><mo>(</mo><mn>0.5</mn><mo>)</mo><mo>-</mo><msub><mi>&theta;</mi><mn>0</mn></msub><mo>)</mo></mrow><mi>&theta;</mi></mfrac></mrow>]]></math><img file="FDA0001128856230000012.GIF" wi="512" he="141" /></maths>0.60≤C<sub>FL</sub>≤0.80选定C<sub>1</sub>,C<sub>2</sub>,C<sub>3</sub>和C<sub>4</sub>,从而确定无量纲弯角分布函数φ(t);S3:定义叶栅的稠度S为叶型弦长C与栅距η之比,选定叶型在设计点的迎角i,通过下式确定叶型的弯角θ进而确定叶型的中弧线;<maths num="0003"><math><![CDATA[<mrow><mi>&delta;</mi><mo>=</mo><msub><mi>f</mi><mn>2</mn></msub><mrow><mo>(</mo><msub><mi>C</mi><mrow><mi>F</mi><mi>L</mi></mrow></msub><mo>)</mo></mrow><mfrac><mi>&theta;</mi><msup><mi>S</mi><mi>n</mi></msup></mfrac><mo>,</mo><mi>n</mi><mo>=</mo><mn>1.1</mn></mrow>]]></math><img file="FDA0001128856230000021.GIF" wi="582" he="135" /></maths>式中,δ为落后角,S为稠度,f<sub>2</sub>(C<sub>FL</sub>)为待定函数,f<sub>2</sub>(C<sub>FL</sub>)=A<sub>1</sub>C<sub>FL</sub><sup>2</sup>+A<sub>2</sub>C<sub>FL</sub>+A<sub>3</sub>其中A<sub>1</sub>,A<sub>2</sub>和A<sub>3</sub>为需要确定的系数;S4:确定最大相对厚度T、叶型的前缘厚度t<sub>LE1</sub>和尾缘厚度t<sub>TE2</sub>,对于稠度S在0.80‑1.6范围内的叶栅,最大相对厚度T:<maths num="0004"><math><![CDATA[<mrow><mi>T</mi><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>B</mi><mn>1</mn></msub><mo>&CenterDot;</mo><msub><mi>Ma</mi><mn>1</mn></msub><mo>+</mo><msub><mi>B</mi><mn>2</mn></msub><mo>,</mo></mrow></mtd><mtd><mrow><mn>0.8</mn><mo>&le;</mo><mi>S</mi><mo>&le;</mo><mn>1.1</mn></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>B</mi><mn>1</mn></msub><mo>&CenterDot;</mo><msub><mi>Ma</mi><mn>1</mn></msub><mo>+</mo><msub><mi>B</mi><mn>2</mn></msub><mo>+</mo><msub><mi>B</mi><mn>3</mn></msub><mo>&CenterDot;</mo><mrow><mo>(</mo><mi>S</mi><mo>-</mo><mn>1.1</mn><mo>)</mo></mrow><mo>,</mo></mrow></mtd><mtd><mrow><mn>1.1</mn><mo>&le;</mo><mi>S</mi><mo>&le;</mo><mn>1.6</mn></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0001128856230000022.GIF" wi="1342" he="135" /></maths>式中,B<sub>1</sub>,B<sub>2</sub>和B<sub>3</sub>为需要确定的系数;前缘厚度 t<sub>LE1</sub>=D<sub>1</sub>·T (0.05≤D<sub>1</sub>≤0.40);尾缘厚度 t<sub>TE2</sub>=D<sub>2</sub>·T (0.05≤D<sub>2</sub>≤0.40);S5:确定最大相对厚度位置P<sub>max</sub>,且P<sub>max</sub>为0.1‑0.3,对于设计点进口气流马赫数0.40≤Ma<sub>1</sub>≤0.60的可调进口导流叶片,最大相对厚度位置:<maths num="0005"><math><![CDATA[<mrow><msub><mi>P</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>E</mi><mn>1</mn></msub><mo>&CenterDot;</mo><msub><mi>Ma</mi><mn>1</mn></msub><mo>+</mo><msub><mi>E</mi><mn>2</mn></msub><mo>,</mo></mrow></mtd><mtd><mrow><mn>0.8</mn><mo>&le;</mo><mi>S</mi><mo>&le;</mo><mn>1.1</mn></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>E</mi><mn>1</mn></msub><mo>&CenterDot;</mo><msub><mi>Ma</mi><mn>1</mn></msub><mo>+</mo><msub><mi>E</mi><mn>2</mn></msub><mo>+</mo><msub><mi>E</mi><mn>3</mn></msub><mo>&CenterDot;</mo><mrow><mo>(</mo><mi>S</mi><mo>-</mo><mn>1.1</mn><mo>)</mo></mrow><mo>,</mo></mrow></mtd><mtd><mrow><mn>1.1</mn><mo>&le;</mo><mi>S</mi><mo>&le;</mo><mn>1.6</mn></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0001128856230000023.GIF" wi="1407" he="143" /></maths>式中,E<sub>1</sub>,E<sub>2</sub>和E<sub>3</sub>为需要确定的系数;S6:确定厚度分布函数f<sub>3</sub>(t),根据选定的前缘厚度、尾缘厚度以及最大相对厚度在厚度分布上的位置,确定厚度分布函数且该厚度分布函数需要满足:<maths num="0006"><math><![CDATA[<mrow><mn>5</mn><mo>&le;</mo><mfrac><mrow><mi>d</mi><mrow><mo>(</mo><msub><mi>f</mi><mn>3</mn></msub><mo>(</mo><mi>t</mi><mo>)</mo><mo>/</mo><mi>T</mi><mo>)</mo></mrow></mrow><mrow><mi>d</mi><mi>t</mi></mrow></mfrac><msub><mo>|</mo><mrow><mi>t</mi><mo>=</mo><mn>0</mn></mrow></msub><mo>&le;</mo><mn>50</mn></mrow>]]></math><img file="FDA0001128856230000024.GIF" wi="614" he="167" /></maths><maths num="0007"><math><![CDATA[<mrow><mfrac><mrow><msup><mi>d</mi><mn>2</mn></msup><msub><mi>f</mi><mn>3</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mrow><msup><mi>dt</mi><mn>2</mn></msup></mrow></mfrac><mo>&lt;</mo><mn>0</mn><mo>,</mo><mi>t</mi><mo>&Element;</mo><mo>&lsqb;</mo><mn>0</mn><mo>,</mo><mn>1</mn><mo>&rsqb;</mo></mrow>]]></math><img file="FDA0001128856230000025.GIF" wi="516" he="143" /></maths>S7:将已经获得的相对厚度分布在中弧线法向方向对称积叠获得叶型的叶身部分;S8:确定前缘及尾缘的形状,给S7获得的叶身部分加上前缘及尾缘;S9:对S8获得的基元叶型用计算流体力学软件进行数值仿真,计算其性能,根据计算结果判断设计点在整条性能曲线上的位置是否合适,如果不合适则按需求改变S3中的设计点迎角,然后重复步骤S1至S8,直至设计点在整条性能曲线上的位置是适当的。
地址 100191 北京市海淀区学院路37号