发明名称 一种抑制塔架振动方法与系统
摘要 本发明公开了一种抑制塔架振动方法与系统,应用所述抑制塔架振动方法与系统,通过对前后加速度消差,计算补偿变桨角度系数得到需要调整的补偿变桨角度进行变桨控制,在不改变机组运行状态的基础上,可以准确、稳定地通过调整补偿变桨角度抑制塔架前后振动。
申请公布号 CN104533717B 申请公布日期 2017.05.03
申请号 CN201410765055.1 申请日期 2014.12.12
申请人 中船重工(重庆)海装风电设备有限公司 发明人 杨微;陶友传;陶芬;刘杰;韩花丽
分类号 F03D7/00(2006.01)I 主分类号 F03D7/00(2006.01)I
代理机构 北京集佳知识产权代理有限公司 11227 代理人 罗满
主权项 一种抑制塔架振动方法,其特征在于,应用于风电机组,所述方法包括:测量机舱前后加速度;消除测量到的所述加速度的数据零点漂移和含有长期趋势向的干扰信号得到消差后的加速度;对所述消差后的加速度进行积分得到前后振动速度;根据所述前后振动速度截取得到塔架一阶固有频率对应速度;以上一周期变桨补偿角度<img file="FDA0001232389840000015.GIF" wi="123" he="63" />作为输入,其绝对值与变桨补偿最大限制值θ<sub>comp_max</sub>之差比上所述变桨补偿最大限制值θ<sub>comp_max</sub>得到变桨补偿角度余量系数K<sub>re</sub>,其关系式为:<maths num="0001"><math><![CDATA[<mrow><msub><mi>K</mi><mrow><mi>r</mi><mi>e</mi></mrow></msub><mo>=</mo><mrow><mo>(</mo><msub><mi>&theta;</mi><mrow><mi>c</mi><mi>o</mi><mi>m</mi><mi>p</mi><mo>_</mo><mi>max</mi></mrow></msub><mo>-</mo><mo>|</mo><msubsup><mi>&theta;</mi><mrow><mi>c</mi><mi>o</mi><mi>m</mi><mi>p</mi></mrow><mo>*</mo></msubsup><mo>|</mo><mo>)</mo></mrow><mo>/</mo><msub><mi>&theta;</mi><mrow><mi>c</mi><mi>o</mi><mi>m</mi><mi>p</mi><mo>_</mo><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub></mrow>]]></math><img file="FDA0001232389840000011.GIF" wi="888" he="82" /></maths>其中,所述变桨补偿最大限制值θ<sub>comp_max</sub>取值根据风电机组振动幅值确定;依据风电机组测量的变桨速度V<sub>pitchspeed</sub>和变桨角度θ<sub>PI</sub>计算最大变桨速度限制余量系数K<sub>v</sub>和变桨角度余量系数K<sub>θ</sub>,其关系式为: K<sub>v</sub>=(V<sub>pitchspeed_max</sub>‑|V<sub>pitchspeed</sub>|‑V<sub>pitchspeed_re</sub>)/(V<sub>pitchspeed_max</sub>‑V<sub>pitchspeed_re</sub>)其中,所述变桨速度V<sub>pitchspeed_max</sub>为变桨系统允许在正常运行过程中的最大变桨速度,V<sub>pitchspeed_re</sub>是变桨速度冗余量;计算变桨角度余量系数K<sub>θ</sub>、最大和最小变桨角度,其关系式为:θ<sub>PI_max</sub>=θ<sub>max</sub>‑θ<sub>PI</sub>‑θ<sub>re</sub>θ<sub>PI_min</sub>=θ<sub>PI</sub>‑θ<sub>min</sub>‑θ<sub>re</sub>K<sub>θ</sub>=min(θ<sub>PI_max</sub>,θ<sub>PI_min</sub>)/θ<sub>comp_max</sub>K<sub>v_θ</sub>=min(K<sub>θ</sub>,K<sub>v</sub>)K<sub>v_θ</sub>=0,if K<sub>θ</sub><0 or K<sub>v</sub><0其中,θ<sub>max</sub>为变桨系统允许在正常运行过程中的最大变桨角度,θ<sub>min</sub>为变桨系统允许在正常运行过程中的最小变桨角度,θ<sub>re</sub>是变桨角度冗余量;K<sub>v_θ</sub>为计算变桨需求余量系数,K<sub>re</sub>为变桨补偿角度余量系数;依据所述计算变桨需求余量系数K<sub>v_θ</sub>和所述变桨补偿角度余量系数K<sub>re</sub>计算得到变桨余量总体尺度系数K<sub>c</sub>,关系式如下:<maths num="0002"><math><![CDATA[<mrow><msub><mi>K</mi><mi>c</mi></msub><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>K</mi><mrow><mi>v</mi><mo>_</mo><mi>&theta;</mi></mrow></msub><mo>,</mo></mrow></mtd><mtd><mrow><mi>i</mi><mi>f</mi></mrow></mtd><mtd><mrow><msub><mi>K</mi><mrow><mi>r</mi><mi>e</mi></mrow></msub><mo>&gt;</mo><msub><mi>K</mi><mrow><mi>v</mi><mo>_</mo><mi>&theta;</mi></mrow></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>K</mi><mrow><mi>r</mi><mi>e</mi></mrow></msub><mo>,</mo></mrow></mtd><mtd><mrow><mi>i</mi><mi>f</mi></mrow></mtd><mtd><mrow><msub><mi>K</mi><mrow><mi>r</mi><mi>e</mi></mrow></msub><mo>&le;</mo><msub><mi>K</mi><mrow><mi>v</mi><mi>&theta;</mi></mrow></msub></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0001232389840000012.GIF" wi="623" he="127" /></maths>结合关系式<img file="FDA0001232389840000013.GIF" wi="572" he="63" />及系统阻尼比关系式<img file="FDA0001232389840000014.GIF" wi="185" he="92" />得到变桨系统附加阻尼比ξ<sub>f</sub>与补偿角度Δθ之间的关系为:<maths num="0003"><math><![CDATA[<mrow><mi>&Delta;</mi><mi>&theta;</mi><mo>=</mo><mfrac><mrow><mo>-</mo><mn>2</mn><msub><mi>&omega;M&xi;</mi><mi>f</mi></msub><mover><mi>x</mi><mo>&CenterDot;</mo></mover></mrow><mrow><msub><mi>&Delta;F</mi><mi>a</mi></msub><mo>/</mo><mi>&Delta;</mi><mi>&theta;</mi></mrow></mfrac></mrow>]]></math><img file="FDA0001232389840000021.GIF" wi="305" he="115" /></maths>其中:M为机组质量,D为阻尼,K为等效刚度,x为塔架顶部前后位移,<img file="FDA0001232389840000022.GIF" wi="32" he="51" />为塔架顶部前后移动速度,<img file="FDA0001232389840000023.GIF" wi="35" he="43" />为塔架顶部前后移动加速度,F<sub>a</sub>为作用在塔架顶部的水平气动推力;变桨系统附加阻尼比ξ<sub>f</sub>根据塔架振动幅度确定,据此得到非线性变桨补偿放大系数<img file="FDA0001232389840000024.GIF" wi="312" he="105" />依据所述变桨余量总体尺度系数K<sub>c</sub>和所述非线性变桨补偿放大系数K<sub>un</sub>后可得到补偿系数K,其关系式为:K=K<sub>c</sub>×K<sub>un</sub>;将所述补偿系数与所述截取得到的塔架一阶固有频率对应速度相乘得到抑制塔架前后振动的目标变桨补偿角度;计算所述目标变桨补偿角度与所述变桨角度的角度和;根据所述角度和进行变桨控制。
地址 401122 重庆市北部新区经开园金渝大道30号