发明名称 变速变桨风力发电机组基于内模控制的最大功率追踪控制方法
摘要 本发明公开一种变速变桨风力发电机组基于内模控制的最大功率追踪控制方法。本发明通过对风力发电机组阻力矩响应系统进行建模得出参考模型,基于参考模型设计内模控制器并给出其解析表达式,基于内模控制方式实现低于额定风速以下轮毂与主轴低速轴连接处最优阻力矩控制,相比于传统的最大功率跟踪算法,本发明有效的降低了由于气动力矩波动、测量误差和建模误差等因素对最优阻力矩跟踪的影响,更高效的实现风轮风能吸收最大功率追踪。针对常用的单位负反馈单回路控制系统,本发明设计了基于内模控制的等效单回路单位负反馈控制系统并给出控制器解析表达式。
申请公布号 CN103362736B 申请公布日期 2015.09.02
申请号 CN201210097356.2 申请日期 2012.04.05
申请人 北京能高自动化技术股份有限公司 发明人 贾利民;刘展;雷涛
分类号 G05B13/00(2006.01)I;F03D7/00(2006.01)I 主分类号 G05B13/00(2006.01)I
代理机构 北京正理专利代理有限公司 11257 代理人 张雪梅
主权项 变速变桨风力发电机组基于内模控制的最大功率追踪控制方法,其特征是该方法考虑风力发电机组传动系统为柔性环节且存在摩擦的情况下,通过对风力发电机组传动系统进行建模,基于内模控制方式实现额定风速以下轮毂与主轴低速轴连接处最优阻力矩控制,降低了由于气动力矩波动、测量误差和建模误差因素对最优阻力矩跟踪的影响,有效的实现风轮风能吸收最大功率追踪;传动系统可以等效为低速轴和高速轴两质量块环节,按照以下方法建模:<maths num="0001" id="cmaths0001"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><mrow><mfrac><mrow><msub><mi>d&omega;</mi><mi>b</mi></msub></mrow><mrow><mi>d</mi><mi>t</mi></mrow></mfrac><mo>=</mo><mfrac><msub><mi>T</mi><mi>b</mi></msub><msub><mi>J</mi><mi>b</mi></msub></mfrac><mo>-</mo><mfrac><mi>T</mi><msub><mi>J</mi><mi>b</mi></msub></mfrac></mrow></mtd></mtr><mtr><mtd><mrow><mfrac><mrow><msub><mi>d&omega;</mi><mi>g</mi></msub></mrow><mrow><mi>d</mi><mi>t</mi></mrow></mfrac><mo>=</mo><mfrac><mrow><mi>&eta;</mi><mi>T</mi></mrow><mrow><msub><mi>NJ</mi><mi>g</mi></msub></mrow></mfrac><mo>-</mo><mfrac><msub><mi>T</mi><mi>e</mi></msub><msub><mi>J</mi><mi>g</mi></msub></mfrac></mrow></mtd></mtr><mtr><mtd><mrow><mfrac><mrow><mi>d</mi><mi>T</mi></mrow><mrow><mi>d</mi><mi>t</mi></mrow></mfrac><mo>=</mo><mi>K</mi><mrow><mo>(</mo><msub><mi>&omega;</mi><mi>b</mi></msub><mo>-</mo><mfrac><msub><mi>&omega;</mi><mi>g</mi></msub><mi>N</mi></mfrac><mo>)</mo></mrow><mo>+</mo><mi>B</mi><mrow><mo>(</mo><mfrac><mrow><msub><mi>d&omega;</mi><mi>b</mi></msub></mrow><mrow><mi>d</mi><mi>t</mi></mrow></mfrac><mo>-</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><mfrac><mrow><msub><mi>d&omega;</mi><mi>b</mi></msub></mrow><mrow><mi>d</mi><mi>t</mi></mrow></mfrac><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000750116800000011.GIF" wi="832" he="493" /></maths>式中ω<sub>b</sub>表示低速轴转速;T<sub>b</sub>表示风轮气动力矩;J<sub>b</sub>表示低速轴转动惯量;T表示轮毂与主轴低速轴连接处与旋转方向相反的阻力矩;ω<sub>g</sub>表示高速轴转速;J<sub>g</sub>表示高速轴转动惯量;T<sub>e</sub>表示发电机电磁力矩;N表示传动系统增速比;η表示传动系统传动效率;K表示传动系统等效到低速轴刚性系数;B表示传动系统等效到低速轴阻尼系数;设状态变量为<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>x</mi><mo>=</mo><mfenced open = '[' close = ']'><mtable><mtr><mtd><msub><mi>&omega;</mi><mi>b</mi></msub></mtd></mtr><mtr><mtd><msub><mi>&omega;</mi><mi>g</mi></msub></mtd></mtr><mtr><mtd><mi>T</mi></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000750116800000012.GIF" wi="219" he="232" /></maths>输入变量为<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>u</mi><mo>=</mo><mfenced open = '[' close = ']'><mtable><mtr><mtd><msub><mi>T</mi><mi>b</mi></msub></mtd></mtr><mtr><mtd><msub><mi>T</mi><mi>e</mi></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000750116800000013.GIF" wi="202" he="161" /></maths>输出为y=T,则传动系统表示为一个三阶模型:<maths num="0004" id="cmaths0004"><math><![CDATA[<mfenced open = '{' close = ''><mtable><mtr><mtd><mrow><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mo>=</mo><mfenced open = '[' close = ']'><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mrow><mo>-</mo><mfrac><mn>1</mn><msub><mi>J</mi><mi>b</mi></msub></mfrac></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mfrac><mi>N</mi><msub><mi>J</mi><mi>t</mi></msub></mfrac></mtd></mtr><mtr><mtd><mi>K</mi></mtd><mtd><mrow><mo>-</mo><mfrac><mi>K</mi><mi>N</mi></mfrac></mrow></mtd><mtd><mrow><mo>-</mo><mfrac><mi>B</mi><msub><mi>J</mi><mi>b</mi></msub></mfrac><mo>-</mo><mfrac><mi>B</mi><msub><mi>J</mi><mi>t</mi></msub></mfrac></mrow></mtd></mtr></mtable></mfenced><mi>x</mi><mo>+</mo><mfenced open = '[' close = ']'><mtable><mtr><mtd><mfrac><mn>1</mn><msub><mi>J</mi><mi>b</mi></msub></mfrac></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mo>-</mo><mfrac><mn>1</mn><msub><mi>J</mi><mi>g</mi></msub></mfrac></mrow></mtd></mtr><mtr><mtd><mfrac><mi>B</mi><msub><mi>J</mi><mi>b</mi></msub></mfrac></mtd><mtd><mfrac><mi>B</mi><mrow><msub><mi>NJ</mi><mi>g</mi></msub></mrow></mfrac></mtd></mtr></mtable></mfenced><mi>u</mi></mrow></mtd></mtr><mtr><mtd><mrow><mi>y</mi><mo>=</mo><mfenced open = '[' close = ']'><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mi>x</mi></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000750116800000021.GIF" wi="919" he="570" /></maths>其中<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>J</mi><mi>t</mi></msub><mo>=</mo><mfrac><mrow><msup><mi>N</mi><mn>2</mn></msup><msub><mi>J</mi><mi>g</mi></msub></mrow><mi>&eta;</mi></mfrac></mrow>]]></math><img file="FDA0000750116800000022.GIF" wi="234" he="155" /></maths>J<sub>t</sub>表示低速轴等效到高速轴后的转动惯量;则系统传递函数表示为:<maths num="0006" id="cmaths0006"><math><![CDATA[<mfenced open = '{' close = ''><mtable><mtr><mtd><mi>y</mi><mo>(</mo><mi>s</mi><mo>)</mo><mo>=</mo><mo>&lsqb;</mo><msub><mi>G</mi><mn>1</mn></msub><mo>(</mo><mi>s</mi><mo>)</mo><mo>,</mo><msub><mi>G</mi><mn>2</mn></msub><mo>(</mo><mi>s</mi><mo>)</mo><mo>&rsqb;</mo><mi>u</mi><mo>(</mo><mi>s</mi><mo>)</mo></mtd></mtr><mtr><mtd><msub><mi>G</mi><mn>1</mn></msub><mo>(</mo><mi>s</mi><mo>)</mo><mo>=</mo><mfrac><mrow><msub><mi>BJ</mi><mi>t</mi></msub><mi>s</mi><mo>+</mo><msub><mi>KJ</mi><mi>t</mi></msub></mrow><mrow><msub><mi>J</mi><mi>b</mi></msub><msub><mi>J</mi><mi>t</mi></msub><msup><mi>s</mi><mn>2</mn></msup><mo>+</mo><mi>B</mi><mrow><mo>(</mo><msub><mi>J</mi><mi>b</mi></msub><mo>+</mo><msub><mi>J</mi><mi>t</mi></msub><mo>)</mo></mrow><mi>s</mi><mo>+</mo><mi>K</mi><mrow><mo>(</mo><msub><mi>J</mi><mi>b</mi></msub><mo>+</mo><msub><mi>J</mi><mi>t</mi></msub><mo>)</mo></mrow></mrow></mfrac></mtd></mtr><mtr><mtd><msub><mi>G</mi><mn>2</mn></msub><mo>(</mo><mi>s</mi><mo>)</mo><mo>=</mo><mfrac><mrow><msub><mi>BJ</mi><mi>b</mi></msub><msub><mi>J</mi><mi>t</mi></msub><mi>s</mi><mo>+</mo><msub><mi>KJ</mi><mi>b</mi></msub><msub><mi>J</mi><mi>t</mi></msub></mrow><mrow><msub><mi>NJ</mi><mi>g</mi></msub><msub><mi>J</mi><mi>b</mi></msub><msub><mi>J</mi><mi>t</mi></msub><msup><mi>s</mi><mn>2</mn></msup><mo>+</mo><msub><mi>BNJ</mi><mi>g</mi></msub><mrow><mo>(</mo><msub><mi>J</mi><mi>b</mi></msub><mo>+</mo><msub><mi>J</mi><mi>t</mi></msub><mo>)</mo></mrow><mi>s</mi><mo>+</mo><msub><mi>KNJ</mi><mi>g</mi></msub><mrow><mo>(</mo><msub><mi>J</mi><mi>b</mi></msub><mo>+</mo><msub><mi>J</mi><mi>t</mi></msub><mo>)</mo></mrow></mrow></mfrac></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000750116800000023.GIF" wi="1144" he="477" /></maths>
地址 100044 北京市海淀区西直门外北京交通大学西门交大知行大厦六层
您可能感兴趣的专利