一种振动监测方法技术

本发明专利技术公开了一种振动监测方法，包括步骤：步骤一，选取叶片的尖端为振动的参考点，使图像采集装置的镜头对准参考点，以获取振动图像；步骤二，将振动图像分解成帧，并对每一帧进行预处理，以消除噪声；步骤三，选取参考点在每一帧的像素坐标，并对像素坐标进行转化，记录参考点转换后的像素坐标，并绘制转化后的像素坐标值与时间t的关系图；步骤四，对关系图进行曲线拟合，并对拟合曲线f(t)做傅里叶变化得到F(ω)。将图像采集装置放置于整流罩，并调整角度使其对准待监测的叶片的参考点，既保证了采集图像的清晰度，亦能够实时采集图像；且采用上述方法，绘制像素坐标值与时间t的关系图，并通过拟合及傅里叶变化得到振动频率ω的组成和占比。

【技术实现步骤摘要】
一种振动监测方法
本专利技术涉及风力发电
，更具体地说，涉及一种振动监测方法。
技术介绍
叶片是风力发电机组的关键部件，占整个风力发电机制造成本的20％左右，在运行过程中，叶片会连续受到空气动力、惯性力等交变载荷的冲击，从而产生不规则的振动甚至是变形。如果振动异常或者叶片形变过于剧烈，将会导致叶片的损毁造成巨大的经济损失。因此实时监测风机叶片振动状况显得尤其重要。对风力机叶片的振动监测的研究已经引起重视，就现有技术而言，主要分为接触式法和非接触式法。前者大都采用压电陶瓷、加速度计传感器、陀螺仪传感器等放置在叶片表面进行测量；该方法的缺点是当传感器放置在叶片表面时会影响风机叶片的气动性能，并且随着风力发电机单机功率的不断提高，叶片的体积也越来越大，这种方法逐渐失去实用性。后者则主要包括激光测振和机器视觉等方法，其中以摄影测量为主的机器视觉法因设备安装方便，不会对叶片的气动性能产生影响，且在实现动态测量的基础上处理精度高可以达到亚像素级，而成为当今研究的热点。但是，就目前研究现状而言，机器视觉法尚未实际应用于风力发电机叶片的振动监测，因叶片图像采集及振动频率计算均属于研究的难点。因此，如何对风机叶片的振动进行监测，是现阶段该领域亟待解决的难题。
技术实现思路
有鉴于此，本专利技术的目的在于提供一种振动监测方法，该方法能够实现对风机叶片的振动进行监测，并获取振动频率，解决了现阶段该领域的难题。一种振动监测方法，包括步骤：步骤一，选取叶片的尖端为振动的参考点，将图像采集装置固定于风力发电机的整流罩、并使所述图像采集装置的镜头对准所述参考点，以获取所述参考点的振动图像；步骤二，将所述振动图像分解成帧，并对每一帧进行预处理，以消除噪声；步骤三，选取所述参考点在每一帧的像素坐标，并对所述像素坐标进行转化，记录所述参考点转换后的像素坐标，并绘制转化后的像素坐标值与时间t的关系图；步骤四，对所述关系图进行曲线拟合，得到拟合曲线f(t)，并对所述拟合曲线f(t)做傅里叶变化得到F(ω)，并绘制进F(ω)的曲线图，进而获得该时间段内振动频率ω的成分和占比。优选的，所述的振动监测方法，步骤三中对所述像素坐标进行转化的过程为：设所述参考点的世界坐标系为(xw，yw，zw)，图像采集装置坐标系为(xc，yc，zc)，则从世界坐标系到图像采集装置坐标系数学转化关系为：其中R为旋转矩阵，T为平移向量；T＝[txtytz]T(3)其中θ，ψ分别为世界坐标系x轴，y轴，z轴绕自身的旋转角度，tx，ty，tz分别为世界坐标系x轴，y轴，z轴相对于原坐标系平移的距离；图像采集装置坐标(X，Y)到图像坐标(xc，yc，zc)的转化关系为：其中f为图像采集装置的焦距；图像坐标(xc，yc，zc)到像素坐标(u，v)的转化关系为：式中dx，dy分别为x轴，y轴方向上的像元尺寸，(u0，v0)为图像原点坐标。优选的，所述的振动监测方法，所述步骤四中，采用正弦函数拟合，得到所述拟合曲线f(t)，f(t)＝689.3*sin(0.5517*t-1.327)+825.7*sin(2.607*t+1.243)+1283*sin(0.791*t+0.518)+721.6*sin(2.639*t+4.229)+771.3*sin(0.8907*t+3.103)+16.91*sin(4.037*t+1.5)+11.33*sin(6.941*t+2.482)+12.92*sin(4.364*t+2.082)(6)则，其中ω为振动频率；t为时间；i为虚数单位；e为自然对数函数的底数，约为2.718。优选的，所述的振动监测方法，所述参考点通过在所述叶片的尖端涂抹反光材料实现标定。优选的，所述的振动监测方法，所述步骤二中预处理的步骤包括：对所述振动图像进行灰度化，得到所述叶片的灰度图像；并对所述灰度图像进行滤波处理，滤除因外界环境变化所产生的椒盐噪声和因设备本身问题所产生的高斯噪声。优选的，所述的振动监测方法，所述图像采集装置为80帧/s的双目工业相机。优选的，所述的振动监测方法，采集所述振动图像的采集时间为10s。优选的，所述的振动监测方法，所述步骤三中，所述像素坐标值为水平方向的像素坐标值，或者为垂直方向的像素坐标值。本专利技术提出的振动监测方法，包括步骤：步骤一，选取叶片的尖端为振动的参考点，将图像采集装置固定于风力发电机的整流罩、并使图像采集装置的镜头对准参考点，以获取参考点的振动图像；步骤二，将振动图像分解成帧，并对每一帧进行预处理，以消除噪声；步骤三，选取参考点在每一帧的像素坐标，并对像素坐标进行转化，记录参考点转换后的像素坐标，并绘制转化后的像素坐标值与时间t的关系图；步骤四，对关系图进行曲线拟合，得到拟合曲线f(t)，并对拟合曲线f(t)做傅里叶变化得到F(ω)，并绘制进F(ω)的曲线图，进而获得该时间段内振动频率ω的成分和占比。风机运行时，叶片尖端的振动较为剧烈，幅度较大，将图像采集装置放置于整流罩，并调整角度使其对准待监测的叶片的参考点，既保证了采集图像的清晰度，亦能够实时采集图像；且采用上述方法，绘制像素坐标值与时间t的关系图，并通过拟合及傅里叶变化得到振动频率ω的组成和占比。因此，本专利技术提出的振动监测方法，能够实现对风机叶片的振动进行监测，并获取振动频率，解决了现阶段该领域的难题。附图说明为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。图1为本专利技术具体实施方式中风机与图像采集装置结合的示意图；图2为本专利技术具体实施方式中叶片与图像采集装置结合的放大图；图3为本专利技术具体实施方式中虚线为像素坐标值与时间t的关系图，实线为f(t)拟合曲线图；图4为本专利技术具体实施方式中傅里叶变化F(ω)的曲线图。图1-图4中：图像采集装置-1；整流罩-2；叶片-3。具体实施方式本具体实施方式的核心在于提供一种振动监测方法，该方法能够实现对风机叶片的振动进行监测，并获取振动频率，解决了现阶段该领域的难题。以下，参照附图对实施例进行说明。此外，下面所示的实施例不对权利要求所记载的
技术实现思路
起任何限定作用。另外，下面实施例所表示的构成的全部内容不限于作为权利要求所记载的专利技术的解决方案所必需的。本具体实施方式提供的振动监测方法，包括步骤：步骤一，选取叶片3的尖端为振动的参考点，将图像采集装置1固定于风力发电机的整流罩2、并使图像采集装置1的镜头对准参考点，以获取参考点的振动图像；步骤二，将振动图像分解成帧，并对每一帧进行预处理，以消除噪声；步骤三，选取参考点本文档来自技高网...

【技术保护点】
1.一种振动监测方法，其特征在于，包括步骤：/n步骤一，选取叶片(3)的尖端为振动的参考点，将图像采集装置(1)固定于风力发电机的整流罩(2)、并使所述图像采集装置(1)的镜头对准所述参考点，以获取所述参考点的振动图像；/n步骤二，将所述振动图像分解成帧，并对每一帧进行预处理，以消除噪声；/n步骤三，选取所述参考点在每一帧的像素坐标，并对所述像素坐标进行转化，记录所述参考点转换后的像素坐标，并绘制转化后的像素坐标值与时间t的关系图；/n步骤四，对所述关系图进行曲线拟合，得到拟合曲线f(t)，并对所述拟合曲线f(t)做傅里叶变化得到F(ω)，并绘制进F(ω)的曲线图，进而获得该时间段内振动频率ω的成分和占比。/n

【技术特征摘要】
1.一种振动监测方法，其特征在于，包括步骤：
步骤一，选取叶片(3)的尖端为振动的参考点，将图像采集装置(1)固定于风力发电机的整流罩(2)、并使所述图像采集装置(1)的镜头对准所述参考点，以获取所述参考点的振动图像；
步骤二，将所述振动图像分解成帧，并对每一帧进行预处理，以消除噪声；
步骤三，选取所述参考点在每一帧的像素坐标，并对所述像素坐标进行转化，记录所述参考点转换后的像素坐标，并绘制转化后的像素坐标值与时间t的关系图；
步骤四，对所述关系图进行曲线拟合，得到拟合曲线f(t)，并对所述拟合曲线f(t)做傅里叶变化得到F(ω)，并绘制进F(ω)的曲线图，进而获得该时间段内振动频率ω的成分和占比。


2.根据权利要求1所述的振动监测方法，其特征在于，步骤三中对所述像素坐标进行转化的过程为：
设所述参考点的世界坐标系为(xw，yw，zw)，图像采集装置坐标系为(xc，yc，zc)，则从世界坐标系到图像采集装置坐标系数学转化关系为：



其中R为旋转矩阵，T为平移向量；



T＝[txtytz]T(3)
其中θ，ψ分别为世界坐标系x轴，y轴，z轴绕自身的旋转角度，tx，ty，tz分别为世界坐标系x轴，y轴，z轴相对于原坐标系平移的距离；
图像采集装置坐标(X，Y)到图像坐标(xc，yc，zc)的转化关系为：



其中f为图像采集装置的焦距；
图像坐标(xc，yc，zc)到像素坐标(u，v)的转化关系为：



其中dx，dy分别为x轴，y轴方向上的像元尺寸，(u0，...

【专利技术属性】
技术研发人员：钱江波，张佳星，李振宇，张志杰，朱霄珣，姚大伟，李岩，王巧珍，
申请(专利权)人：华北电力大学保定，
类型：发明
国别省市：河北;13