发明名称 一种建立发电机轴系多质量块变阻尼模型的方法及装置
摘要 一种建立发电机轴系多质量块变阻尼模型的方法及装置,属于电力系统实时数字仿真领域,以解决在对现有多质量块模型的仿真计算过程中不能完全反映实际情况的问题。该方法包括确定建立弹性多质量块阻尼模型需要的参数;根据所述参数建立的任意一个质量块旋转运动方程,获得第i个质量块的旋转运动方程;通过梯形方法对第i个质量块的旋转运动方程进行离散化;根据离散化的第i个质量块的旋转运动方程,得到第阶质量块的发电机轴系转子运动方程。本发明具有能够在对现有多质量块模型的仿真计算过程中反映实际情况的优点。
申请公布号 CN101515719B 申请公布日期 2011.12.14
申请号 CN200910081098.7 申请日期 2009.04.02
申请人 华北电力大学 发明人 张鹏;毕天姝;薛安成;杨奇逊
分类号 H02J3/00(2006.01)I;G06F17/50(2006.01)I 主分类号 H02J3/00(2006.01)I
代理机构 北京凯特来知识产权代理有限公司 11260 代理人 郑立明;孟丽娟
主权项 1.一种建立发电机轴系多质量块变阻尼模型的方法,其特征在于,包括:确定建立发电机轴系多质量块变阻尼模型需要的参数;根据所述参数建立的任意一个质量块旋转运动方程,获得第i个质量块的旋转运动方程,其中i表示自然数,所述任意一个质量块旋转运动方程为:<maths num="0001"><![CDATA[<math><mrow><mn>2</mn><mi>H</mi><mover><mi>&omega;</mi><mo>&CenterDot;</mo></mover><mo>=</mo><mi>&Delta;T</mi></mrow></math>]]></maths>其中ω表示旋转角速度,<img file="FDA0000082126880000012.GIF" wi="42" he="73" />表示旋转角加速度,H为质量块围绕旋转轴的惯性时间常数,ΔT为不平衡转矩;通过梯形方法对第i个质量块的旋转运动方程进行离散化;根据离散化的第i个质量块的旋转运动方程,得到第i个质量块的发电机轴系转子运动方程;其中,所述第i个质量块的旋转运动方程为:<maths num="0002"><![CDATA[<math><mrow><mn>2</mn><mi>H</mi><mover><mi>&omega;</mi><mo>&CenterDot;</mo></mover><mo>=</mo><msub><mi>K</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>i</mi></mrow></msub><mrow><mo>(</mo><msub><mi>&Delta;&theta;</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>&Delta;&theta;</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>&Delta;T</mi><mi>i</mi></msub><mo>-</mo><msub><mi>K</mi><mrow><mi>i</mi><mo>,</mo><mi>i</mi><mo>+</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><msub><mi>&Delta;&theta;</mi><mi>i</mi></msub><mo>-</mo><msub><mi>&Delta;&theta;</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mo>-</mo><mi>D&Delta;</mi><msub><mi>&omega;</mi><mi>i</mi></msub></mrow></math>]]></maths>式中的θ<sub>i</sub>表示第i个质量块的扭转角度,θ<sub>i+1</sub>表示第(i+1)个质量块的扭转角度,θ<sub>i-1</sub>表示第(i-1)个质量块的扭转角度,ω<sub>i</sub>表示第i个质量块的旋转角速度,K<sub>i,i+1</sub>表示第i个质量块和第(i+1)个质量块间轴段的刚度,K<sub>i-1,i</sub>表示第(i-1)个质量块和第i个质量块间轴段的刚度,ΔT<sub>i</sub>为第i个质量块的不平衡转矩,K<sub>i</sub>表示第i个轴段刚度,并且K<sub>i-1,i</sub>|<sub>i=1</sub>=0,K<sub>i,i+1</sub>|<sub>i=n</sub>=0,i表示自然数;将第i个质量块的旋转运动方程离散化得到离散化的第i个质量块的旋转运动方程为:<maths num="0003"><![CDATA[<math><mrow><mn>2</mn><mi>H</mi><mfrac><mrow><mi>&omega;</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mi>&omega;</mi><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mi>&Delta;t</mi><mo>)</mo></mrow></mrow><mi>&Delta;t</mi></mfrac><mo>+</mo><mi>D</mi><mfrac><mrow><mi>&omega;</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mi>&omega;</mi><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mi>&Delta;t</mi><mo>)</mo></mrow></mrow><mn>2</mn></mfrac><mo>+</mo><mi>K</mi><mfrac><mrow><mi>&theta;</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mi>&theta;</mi><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mi>&Delta;t</mi><mo>)</mo></mrow></mrow><mn>2</mn></mfrac><mo>=</mo><mfrac><mrow><mi>T</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mi>T</mi><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mi>&Delta;t</mi><mo>)</mo></mrow></mrow><mn>2</mn></mfrac></mrow></math>]]></maths>式中的K为常数,表示轴刚度,转矩T和阻尼D作为已知量,ω(t)和θ(t)为未知量;将离散化的第i个质量块的旋转运动方程和t时刻质量块的扭转角度θ<sub>i</sub>的表达式带入所述的第i个质量块的旋转运动方程,消去变量θ(t),得到第i个质量块的发电机轴系转子运动方程为:<img file="FDA0000082126880000021.GIF" wi="1710" he="84" />式中的[ω<sub>i</sub>(t)]表示第i个质量块转速的本次迭代结果,[T<sub>i</sub>(t)]表示用于本次迭代的作用于第i个质量块的转矩,[ω<sub>i</sub>(t-Δt)]表示第i个质量块转速值的上次迭代结果,[B]表示参数B形成的矩阵,[2K]表示刚度参数2K矩阵,[θ<sub>i</sub>(t-Δt)]表示第i个质量块旋转角度的上次迭代结果,[A<sup>-1</sup>]中A的表达式、参数B以及t时刻质量块的扭转角度θ<sub>i</sub>的表达式分别为:<maths num="0004"><![CDATA[<math><mrow><mi>A</mi><mo>=</mo><mo>[</mo><mfrac><mn>2</mn><mi>&Delta;t</mi></mfrac><mn>2</mn><mi>H</mi><mo>+</mo><mi>D</mi><mo>+</mo><mfrac><mi>&Delta;t</mi><mn>2</mn></mfrac><mi>K</mi><mo>]</mo></mrow></math>]]></maths><maths num="0005"><![CDATA[<math><mrow><mi>B</mi><mo>=</mo><mo>[</mo><mfrac><mn>2</mn><mi>&Delta;t</mi></mfrac><mn>2</mn><mi>H</mi><mo>-</mo><mi>D</mi><mo>+</mo><mfrac><mi>&Delta;t</mi><mn>2</mn></mfrac><mi>K</mi><mo>]</mo></mrow></math>]]></maths><maths num="0006"><![CDATA[<math><mrow><msub><mi>&theta;</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mi>&Delta;t</mi><mn>2</mn></mfrac><mo>[</mo><msub><mi>&omega;</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>&omega;</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mi>&Delta;t</mi><mo>)</mo></mrow><mo>]</mo><mo>+</mo><msub><mi>&theta;</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mi>&Delta;t</mi><mo>)</mo></mrow><mo>.</mo></mrow></math>]]></maths>
地址 102206 北京市昌平区北农路2号