发明名称 用于燃气轮机拉杆式转子扭振模态的计算验证方法
摘要 本发明公开了一种用于燃气轮机拉杆式转子扭振模态的计算验证方法,(1)建立拉杆转子接触段的传递矩阵方程,然后在该方程中引入接触段抗扭刚度修正系数;(2)计算上述方程中拉杆转子接触段抗扭刚度修正系数;(3)利用上述拉杆转子接触段的传递矩阵方程和接触段抗扭刚度修正系数值,通过传递矩阵法计算转子的模态频率和模态振型,计算中考虑了转子接触段接触效应;(4)通过扭振模态实验测试实验拉杆转子在不同预紧力下的扭振模态频率,并与上述计算结果比较。本发明对拉杆式转子的扭振模态的计算能有效地考虑接触效应的影响,可以更为准确地计算拉杆式转子的扭振模态频率,这为合理地设计该类型转子的扭转振动特性提供了可靠依据。
申请公布号 CN102222123B 申请公布日期 2013.04.17
申请号 CN201110062849.8 申请日期 2011.03.16
申请人 西安交通大学 发明人 袁奇;高进;李浦
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 西安通大专利代理有限责任公司 61200 代理人 朱海临
主权项 1.一种用于燃气轮机拉杆式转子扭振模态的计算验证方法,其特征在于,包括下述步骤:1)建立拉杆转子接触段的传递矩阵方程对于接触段长度为L<sub>s</sub>,极惯量矩为I<sub>p</sub>,直径转动惯量为I的拉杆式转子,接触段一端的角位移和扭矩分别为Φ<sub>i</sub>和T<sub>i</sub>,另一端为Φ<sub>i+1</sub>和T<sub>i+1</sub>,考虑接触效应对转子接触段抗扭刚度的影响,其传递矩阵方程为:<maths num="0001"><![CDATA[<math><mrow><msub><mfenced open='{' close='}'><mtable><mtr><mtd><mi>&Phi;</mi></mtd></mtr><mtr><mtd><mi>T</mi></mtd></mtr></mtable></mfenced><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><msub><mfenced open='[' close=']'><mtable><mtr><mtd><mi>cos</mi><mfrac><mi>&omega;</mi><mrow><msqrt><mi>&eta;</mi></msqrt><msub><mi>a</mi><mi>T</mi></msub></mrow></mfrac><msub><mi>L</mi><mi>s</mi></msub></mtd><mtd><mfrac><mn>1</mn><mrow><mi>&eta;</mi><msub><mi>GI</mi><mi>p</mi></msub><mfrac><mi>&omega;</mi><mrow><msqrt><mi>&eta;</mi></msqrt><msub><mi>a</mi><mi>T</mi></msub></mrow></mfrac></mrow></mfrac><mi>sin</mi><mfrac><mi>&omega;</mi><mrow><msqrt><mi>&eta;</mi></msqrt><msub><mi>a</mi><mi>T</mi></msub></mrow></mfrac><msub><mi>L</mi><mi>s</mi></msub></mtd></mtr><mtr><mtd><mo>-</mo><mi>&eta;</mi><msub><mi>GI</mi><mi>p</mi></msub><mfrac><mi>&omega;</mi><mrow><msqrt><mi>&eta;</mi></msqrt><msub><mi>a</mi><mi>T</mi></msub></mrow></mfrac><mi>sin</mi><mfrac><mi>&omega;</mi><mrow><msqrt><mi>&eta;</mi></msqrt><msub><mi>a</mi><mi>T</mi></msub></mrow></mfrac><msub><mi>L</mi><mi>s</mi></msub></mtd><mtd><mi>cos</mi><mfrac><mi>&omega;</mi><mrow><msqrt><mi>&eta;</mi></msqrt><msub><mi>a</mi><mi>T</mi></msub></mrow></mfrac><msub><mi>L</mi><mi>s</mi></msub></mtd></mtr></mtable></mfenced><mi>i</mi></msub><msub><mfenced open='{' close='}'><mtable><mtr><mtd><mi>&Phi;</mi></mtd></mtr><mtr><mtd><mi>T</mi></mtd></mtr></mtable></mfenced><mi>i</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中:ω为自由振动频率,G为拉杆式转子轴段材料的剪切模量,<img file="FDA00002523218900012.GIF" wi="322" he="69" />η为拉杆式转子接触段抗扭刚度修正系数,其表示考虑接触效应的接触段和相应的无接触面连续轴段的抗扭刚度之比;2)计算拉杆式转子接触段抗扭刚度修正系数η<maths num="0002"><![CDATA[<math><mrow><mi>&eta;</mi><mo>=</mo><mfrac><mn>1</mn><mrow><mrow><mo>(</mo><mn>1</mn><mo>/</mo><msub><mi>&eta;</mi><mi>Lc</mi></msub><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>&zeta;</mi><mo>+</mo><mn>1</mn></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>21</mn><mo>)</mo></mrow></mrow></math>]]></maths>式(21)中:ζ=L<sub>c</sub>/L<sub>s</sub>,L<sub>c</sub>为接触层长度,η<sub>Lc</sub>为接触层抗扭刚度修正系数,由下式表示:<maths num="0003"><![CDATA[<math><mrow><msub><mi>&eta;</mi><mi>Lc</mi></msub><mo>=</mo><mfrac><msub><mi>K</mi><mi>t</mi></msub><msub><mi>K</mi><mn>0</mn></msub></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>20</mn><mo>)</mo></mrow></mrow></math>]]></maths>式(20)中:K<sub>0</sub>=GI<sub>p</sub>L<sup>-1</sup>为与接触层同样长度的转子连续轴段的抗扭刚度,其中L为连续轴段的长度,K<sub>t</sub>为接触层的抗扭刚度,由下式表示:<maths num="0004"><![CDATA[<math><mrow><msub><mi>K</mi><mi>t</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>G</mi><mi>T</mi></msub><msub><mi>I</mi><mi>p</mi></msub></mrow><msub><mi>L</mi><mi>c</mi></msub></mfrac><mo>=</mo><mfrac><mrow><mn>2.6256</mn><msup><mi>&beta;</mi><mn>0.1864</mn></msup><msubsup><mi>P</mi><mn>0</mn><mn>0.8136</mn></msubsup><msubsup><mi>I</mi><mi>p</mi><mn>2</mn></msubsup></mrow><mrow><msub><mi>&sigma;</mi><mi>m</mi></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msub><mi>v</mi><mn>0</mn></msub><mo>)</mo></mrow><mrow><mo>(</mo><mn>2</mn><mo>-</mo><msub><mi>v</mi><mn>0</mn></msub><mo>)</mo></mrow><msub><mi>I</mi><mi>T</mi></msub></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>18</mn><mo>)</mo></mrow></mrow></math>]]></maths>式(18)中:G<sub>T</sub>为接触层等效剪切模量,P<sub>0</sub>为接触面名义压力,σ<sub>m</sub>为微凸体高度分布标准差,ν<sub>0</sub>为接触层材料泊松比,<img file="FDA00002523218900016.GIF" wi="377" he="67" />其中m为接触面微凸体总数,A为接触面积,E为当量弹性模量,R<sub>s</sub>为接触面的微凸体平均曲率半径,式(18)中,I<sub>T</sub>由下式表示:<maths num="0005"><![CDATA[<math><mrow><msub><mi>I</mi><mi>T</mi></msub><mo>=</mo><munder><mrow><mo>&Integral;</mo><mo>&Integral;</mo></mrow><mi>A</mi></munder><mfrac><msup><mi>&rho;</mi><mn>2</mn></msup><msup><mi>&xi;</mi><mfrac><mn>1</mn><mn>3</mn></mfrac></msup></mfrac><mi>dA</mi><mo>=</mo><munder><mrow><mo>&Integral;</mo><mo>&Integral;</mo></mrow><mi>A</mi></munder><mfrac><msup><mi>&rho;</mi><mn>2</mn></msup><msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mfrac><mrow><msub><mi>M</mi><mi>T</mi></msub><mi>&rho;</mi></mrow><mrow><msub><mi>&mu;I</mi><mi>p</mi></msub><msub><mi>P</mi><mn>0</mn></msub></mrow></mfrac><mo>)</mo></mrow><mfrac><mn>1</mn><mn>3</mn></mfrac></msup></mfrac><mi>dA</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>19</mn><mo>)</mo></mrow></mrow></math>]]></maths>式(19)中:μ为接触面最大静摩擦系数,MT为接触段所受扭矩,ρ为接触面半径;3)计算拉杆转子在各预紧力下的扭振模态参数利用上述步骤1)~2)得到的接触段传递矩阵方程和接触段抗扭刚度修正系数η,采用传递矩阵法计算得到实验拉杆转子在各预紧力下的扭振模态参数,即模态频率和模态振型;4)对接触段抗扭刚度修正系数η进行实验验证通过扭振模态实验测得拉杆转子在各预紧力下的扭振模态频率,然后与步骤3)计算得到的拉杆转子在各对应预紧力下的扭振模态频率对比,分析计算和实验的结果随预紧力的变化规律是否一致,则可验证接触段抗扭刚度修正系数的计算方法是否准确。
地址 710049 陕西省西安市咸宁西路28号