发明名称 一种针对大型风电塔筒结构应力的在线监测方法
摘要 本发明公开了一种针对大型风电塔筒结构应力的在线监测方法,属于结构健康监测领域。由于风电塔筒尺寸庞大,所处环境恶劣,传统的使用应变片测量结构应力的方法所需应变片数量较多,而且应变片不能长期处于恶劣的环境之中,不适合作为风电塔筒结构的应力监测方法。本发明根据风电塔筒晃动特点,将晃动分解为晃动平衡位置和围绕平衡位置的晃动幅度,分别对应交变应力中的平均应力和应力幅,并建立了晃动与应力之间的关系,从而可以通过监测风电塔筒的晃动来反映风电塔筒的结构应力,解决了应变片测应力的方法受限于风电塔筒的结构尺寸以及工作环境的问题。此外,本发明方便实用,适合大型风电塔筒结构的在线应力监测。
申请公布号 CN104122013B 申请公布日期 2016.06.08
申请号 CN201410335944.4 申请日期 2014.07.15
申请人 西安交通大学 发明人 侯成刚;杜杰
分类号 G01L1/08(2006.01)I;G01P15/08(2006.01)I 主分类号 G01L1/08(2006.01)I
代理机构 西安通大专利代理有限责任公司 61200 代理人 陆万寿
主权项 一种针对大型风电塔筒结构应力的在线监测方法,其特征在于,该方法通过建立风电塔筒晃动与结构应力之间的关系,通过监测塔筒的晃动来进行塔筒结构应力监测;具体包括以下步骤:1)采用安装于风电塔筒机舱的电容式三轴MEMS加速度传感器实时同步采集塔筒的结构振动加速度信号;2)对步骤1)采集的结构振动加速度信号进行处理,分别提取其中的重力加速度分量和运动加速度分量;3)对步骤2)中所提取的重力加速度分量,根据三轴MEMS加速度传感器倾斜测量原理计算得到风电塔筒顶部倾角,再由风电塔筒顶部倾角与其挠度几何关系,得到风电塔筒晃动的平衡位置;对步骤2)中所提取的运动加速度分量,经过二次积分,得到风电塔筒晃动的幅度;4)根据风电塔筒最大应力与风电塔筒晃动的平衡位置及幅度的关系,将晃动平衡位置转化为风电塔筒最大应力点处的平均应力,晃动幅度转化为风电塔筒最大应力点处的应力幅,进而表征风电塔筒晃动过程中的交变应力,实现风电塔筒结构的在线应力监测;风电塔筒的晃动平衡位置和围绕平衡位置的晃动幅度的计算方法为:对于加速度信号中的直流分量,根据三轴MEMS加速度传感器倾角测量原理计算得到晃动处于平衡位置时的弯曲倾角,再根据风电塔筒顶部挠度与弯曲倾角之间的几何关系确定晃动平衡位置;对于运动加速度,经过二次积分,得到风电塔筒的晃动幅度;其中,风电塔筒晃动平衡位置反映平均应力,围绕平衡位置的晃动幅度反映应力幅,根据平均应力和应力幅表征结构的交变应力;风电塔筒的晃动位移S(t)的计算公式如下:S(t)=∫[∫a(t)dt]dt          (1)式中:a为风电塔筒围绕平衡位置晃动的运动加速度,m·s<sup>‑2</sup>;t为风电塔筒围绕平衡位置晃动的采样时间,s;其中,风电塔筒的晃动幅度y<sub>A</sub>取晃动位移S(t)的最大值;风电塔筒晃动处于平衡位置时的弯曲倾角包括风电塔筒塔顶X方向倾角θ<sub>1</sub>和风电塔筒塔顶Y方向倾角θ<sub>2</sub>,其计算公式分别如下:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>&theta;</mi><mn>1</mn></msub><mo>=</mo><mi>a</mi><mi>r</mi><mi>c</mi><mi>t</mi><mi>a</mi><mi>n</mi><mfrac><msub><mi>a</mi><mi>x</mi></msub><msqrt><mrow><msubsup><mi>a</mi><mi>y</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>a</mi><mi>z</mi><mn>2</mn></msubsup></mrow></msqrt></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000913904960000021.GIF" wi="1190" he="149" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>&theta;</mi><mn>2</mn></msub><mo>=</mo><mi>a</mi><mi>r</mi><mi>c</mi><mi>t</mi><mi>a</mi><mi>n</mi><mfrac><msub><mi>a</mi><mi>y</mi></msub><msqrt><mrow><msubsup><mi>a</mi><mi>x</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>a</mi><mi>z</mi><mn>2</mn></msubsup></mrow></msqrt></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000913904960000022.GIF" wi="1198" he="151" /></maths>式中:a<sub>x</sub>、a<sub>y</sub>、a<sub>z</sub>分别为风电塔筒坐标系的X、Y、Z三个方向的重力加速度分量,m·s<sup>‑2</sup>;且有,<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msqrt><mrow><msubsup><mi>a</mi><mi>x</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>a</mi><mi>y</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>a</mi><mi>z</mi><mn>2</mn></msubsup></mrow></msqrt><mo>=</mo><mi>g</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000913904960000023.GIF" wi="1189" he="102" /></maths>其中:g为重力加速度,m·s<sup>‑2</sup>;风电塔筒塔顶挠度和弯曲倾角的几何关系式如下:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>y</mi><mn>1</mn></msub><mo>=</mo><msub><mi>&theta;</mi><mn>1</mn></msub><mo>&CenterDot;</mo><mfrac><mrow><mn>2</mn><mi>l</mi></mrow><mn>3</mn></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000913904960000024.GIF" wi="1172" he="126" /></maths><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>y</mi><mn>2</mn></msub><mo>=</mo><msub><mi>&theta;</mi><mn>2</mn></msub><mo>&CenterDot;</mo><mfrac><mrow><mn>2</mn><mi>l</mi></mrow><mn>3</mn></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000913904960000025.GIF" wi="1155" he="119" /></maths>式中:y<sub>1</sub>为风电塔筒塔顶X方向挠度,m;y<sub>2</sub>为风电塔筒塔顶Y方向挠度,m;l为风电塔筒高度/m;风电塔筒位于平衡位置处的塔顶挠度y<sub>B</sub>的计算公式如下:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>y</mi><mi>B</mi></msub><mo>=</mo><msqrt><mrow><msubsup><mi>y</mi><mn>1</mn><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>y</mi><mn>2</mn><mn>2</mn></msubsup></mrow></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>12</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000913904960000026.GIF" wi="1189" he="95" /></maths>将晃动平衡位置转化为风电塔筒最大应力点处的平均应力σ<sub>m</sub>,其计算公式如下:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>&sigma;</mi><mi>m</mi></msub><mo>=</mo><mo>-</mo><msub><mi>y</mi><mi>B</mi></msub><mo>&CenterDot;</mo><mfrac><mrow><mn>3</mn><mi>E</mi><mi>I</mi></mrow><mrow><msub><mi>W</mi><mi>z</mi></msub><msup><mi>l</mi><mn>2</mn></msup></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>13</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000913904960000027.GIF" wi="1182" he="142" /></maths>式中:E为材料弹性模量,Pa;I为截面惯性矩,m<sup>4</sup>;将晃动幅度转化为风电塔筒最大应力点处的应力幅,其计算公式如下:<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msub><mi>&sigma;</mi><mi>a</mi></msub><mo>=</mo><mo>-</mo><msub><mi>y</mi><mi>A</mi></msub><mo>&CenterDot;</mo><mfrac><mrow><mn>3</mn><mi>E</mi><mi>I</mi></mrow><mrow><msub><mi>W</mi><mi>z</mi></msub><msup><mi>l</mi><mn>2</mn></msup></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>14</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000913904960000031.GIF" wi="1174" he="134" /></maths>至此,就建立了风电塔筒最大应力与塔顶晃动位移的关系,通过监测塔筒的晃动来进行塔筒结构应力监测。
地址 710049 陕西省西安市咸宁西路28号