CN108427827A - 一种基于加速度信号两次积分的输电线舞动轨迹测试方法 - Google Patents
一种基于加速度信号两次积分的输电线舞动轨迹测试方法 Download PDFInfo
- Publication number
- CN108427827A CN108427827A CN201810115565.2A CN201810115565A CN108427827A CN 108427827 A CN108427827 A CN 108427827A CN 201810115565 A CN201810115565 A CN 201810115565A CN 108427827 A CN108427827 A CN 108427827A
- Authority
- CN
- China
- Prior art keywords
- power transmission
- transmission line
- acceleration
- signal
- acceleration signal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
- G01B21/02—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring length, width, or thickness
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种基于加速度信号两次积分的输电线舞动轨迹测试方法。本发明包括步骤:(1)在同一根输电线上除轴向外的其他两个平面内,等距布置若干组加速度传感器,每组两个,互相垂直布置,分别测量输电线在除轴向外的其他两个平面内的振动加速度;(2)对于采集到的加速度信号进行两次积分得到位移信号;(3)用悬链线方程建立输电线的有限元模型,求出输电线模态Ф;(4)根据振动理论中的模态叠加法,输电线任意时刻的振动可以唯一地表示为各阶模态的线性组合。本发明对加速度信号进行两次积分并利用最小二乘法消除趋势项以减小初值和传感器零点漂移带来的影响,从而获得输电线的舞动轨迹。
Description
技术领域:
本发明涉及一种基于加速度信号两次积分的输电线舞动轨迹测试方法,属于输电线路检测技术领域。
背景技术:
高压输电线路长期运行于野外,由于自然、地理条件的复杂性导致其易受大风、地震、覆冰等影响,造成输电线低频、大幅度舞动。当输电线舞动时,经常会造成输电杆塔上金具损毁、输电线路相间闪络等事故,严重时还会造成输电塔倒塔等恶性事故。因此,对输电线振幅进行监测具有重要意义。
目前在电力系统中,对输电线振幅检测主要采用人工巡检、摄像头及图像识别技术、位移传感器直接测量等方法。人工巡检效率低、成本高,智能化程度低;摄像头配合图像识别技术是将摄像头安装在输电杆塔上,将拍摄的图片输入计算机进行处理,从而得到输电线的振幅,这种方法准确度较低,只能实现定性监测;位移传感器直接测量是将位移传感器布置在输电线上直接测量其位移的一种方法,但由于输电线舞动振幅较大,常规的位移传感器量程无法满足要求,而量程较大的拉绳式位移传感器又无法找到固定支点,所以位移传感器也难以准确测量输电线舞动振幅。
由于输电线舞动会造成严重事故,为了掌握输电线振动幅度,避免事故的发生,亟需一种可靠且准确的振幅测量方法。
发明内容
本发明的目的是提供一种基于加速度信号两次积分的输电线舞动轨迹测试方法,对加速度信号进行两次积分并利用最小二乘法消除趋势项以减小初值和传感器零点漂移带来的影响,从而获得输电线的舞动轨迹。
上述的目的通过以下技术方案实现:
一种基于加速度信号两次积分的输电线舞动轨迹测试方法,该方法包括如下步骤:
(1)在同一根输电线上除轴向外的其他两个平面内,等距布置若干组加速度传感器,每组两个,互相垂直布置,分别测量输电线在除轴向外的其他两个平面内的振动加速度;
(2)对于采集到的加速度信号进行两次积分得到位移信号,所述的两次积分,需要使用最小二乘法消除趋势项对原始加速度曲线以及一次积分后的速度曲线进行处理,再次积分后得到位移曲线;
(3)用悬链线方程建立输电线的有限元模型,求出输电线模态Ф;
(4)根据振动理论中的模态叠加法,模态之间相互正交线性独立,输电线任意时刻的振动可以唯一地表示为各阶模态的线性组合:
式中,xpi是描述系统运动的主坐标,只考虑输电线的前3阶模态,则i=3,通过布置3组加速度传感器即可获得任意时刻的xp1,xp2,xp3,从而获得整条输电线的形状。
所述的基于加速度信号两次积分的输电线舞动轨迹测试方法,步骤(1)中所述的加速度传感器布置有三组。
所述的基于加速度信号两次积分的输电线舞动轨迹测试方法,步骤(2)中所述对于采集到的加速度信号进行两次积分得到位移信号的具体方法是:
设初始速度和位移分别为v0和s0,加速度传感器测得的信号为a(t),则真实的加速度应为a(t)-a0,理论计算得到的速度和位移分别为:
v(t)=v0-a0t+∫a(t)dt
再采用多项式最小二乘法消除趋势项带来的误差影响,设一个多项式:
确定该多项式的各待定系数cj(j=0,1,…,m),使得该多项式与采集到的加速度信号a(t)的误差平方和最小,即
满足E有极值的条件为:
依次取E对ci求偏导,可以产生一个m+1元线性方程组:
解方程组,求出m+1个待定系数cj(j=0,1,…,m),
上述各式中,m为设定的多项式阶次,其值范围为0≤j≤m,在实际的振动信号处理中,通常取m=1~3来对采样数据进行多项式趋势项消除的处理。
本发明所产生的有益效果:
1.对于输电线低频、大振幅的舞动,采用加速度传感器进行振幅测量可以满足量程和精度的要求。
2.对加速度信号进行两次积分,使用最小二乘法消除线性项可以避免初值和零点漂移等问题带来的影响。
3.利用模态叠加法,通过有限组加速度传感器可以获得任意时刻整条输电线的舞动形状。
附图说明
图1是本发明加速度传感器安装示意图。
图2是本发明用最小二乘法消除趋势项处理前后的信号处理对比图,其中图2(a)为处理前加速度曲线图;2(b)为处理后加速度曲线;图2(c)为处理前速度曲线;图2(d)为处理后速度曲线;图2(e)为处理前位移曲线;图2(f)为处理后位移曲线。
图3是本发明用悬链线方程建立的输电线模型。
图4是本发明用有限元软件求得的输电线前3阶模态,其中图4(a)为输电线第一阶模态;图4(b)为输电线第二阶模态;图4(c)为输电线第三阶模态。
具体实施方式
下面结合具体实施方式,进一步阐明本发明,应理解下述具体实施方式仅用于说明本发明而不用于限制本发明的范围。
一种基于加速度信号两次积分的输电线舞动轨迹测试方法,该方法包括如下步骤:
(1)在同一根输电线上等距布置三组加速度传感器,每组两个,互相垂直布置,分别测量输电线在两个平面内的振幅,如图1所示。
(2)采集到的加速度信号进行两次积分得到位移信号。然而在实际的测量过程中,采集到的加速度信号不可避免地会带有测量误差和噪声干扰使得初始加速度a0形成偏移。
设初始速度和位移分别为v0和s0,加速度传感器测得的信号为a(t),则真实的加速度应为a(t)-a0,理论计算得到的速度和位移分别为:
v(t)=v0-a0t+∫a(t)dt
由上式可知,实际计算得到的速度和位移分别含有一次趋势项v0-a0t和二次趋势项s0+v0t-0.5a0t2,加速度信号随着积分次数的增加误差放大。
本发明采用多项式最小二乘法消除趋势项带来的误差影响。设一个多项式:
确定该多项式的各待定系数cj(j=0,1,…,m),使得该多项式与采集到的加速度信号a(t)的误差平方和最小,即
满足E有极值的条件为
依次取E对ci求偏导,可以产生一个m+1元线性方程组:
解方程组,求出m+1个待定系数cj(j=0,1,…,m)。上述各式中,m为设定的多项式阶次,其值范围为0≤j≤m。在实际的振动信号处理中,通常取m=1~3来对采样数据进行多项式趋势项消除的处理。用该方法处理前后的曲线对比如图2所示。
(3)用悬链线方程建立输电线的有限元模型如图3所示,同时使用有限元软件求得输电线的前三阶模态如图4所示。
(4)根据振动理论中的模态叠加法,模态之间相互正交线性独立,输电线任意时刻的振动可以唯一地表示为各阶模态的线性组合:
式中,xpi是描述系统运动的主坐标。由于输电线舞动的能量大部分集中于低阶模态,故只考虑输电线的前3阶模态,即n=3。用三组等距布置的加速度传感器可以获得任意时刻输电线在两个运动平面内三个点的运动轨迹。结合输电线的前三阶模态Ф(1),Ф(2),Ф(3),可以求出输电线运动的主坐标xp1,xp2,xp3。用主坐标和模态可以完全描述整条输电线在任意时刻的舞动轨迹。
应当指出,上述实施实例仅仅是为清楚地说明所作的举例,而并非对实施方式的限定,这里无需也无法对所有的实施方式予以穷举。本实施例中未明确的各组成部分均可用现有技术加以实现。对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (3)
1.一种基于加速度信号两次积分的输电线舞动轨迹测试方法,其特征在于,该方法包括如下步骤:
(1)在同一根输电线上除轴向外的其他两个平面内,等距布置若干组加速度传感器,每组两个,互相垂直布置,分别测量输电线在除轴向外的其他两个平面内的振动加速度;
(2)对于采集到的加速度信号进行两次积分得到位移信号,所述的两次积分,需要使用最小二乘法消除趋势项对原始加速度曲线以及一次积分后的速度曲线进行处理,再次积分后得到位移曲线;
(3)用悬链线方程建立输电线的有限元模型,求出输电线模态Ф;
(4)根据振动理论中的模态叠加法,模态之间相互正交线性独立,输电线任意时刻的振动x可以唯一地表示为各阶模态的线性组合:
式中,xpi是描述系统运动的主坐标,只考虑输电线的前3阶模态,则i=3,通过布置3组加速度传感器即可获得任意时刻的xp1,xp2,xp3,从而获得整条输电线的形状。
2.根据权利要求1所述的基于加速度信号两次积分的输电线舞动轨迹测试方法,其特征在于,步骤(1)中所述的加速度传感器布置有三组。
3.根据权利要求1或2所述的基于加速度信号两次积分的输电线舞动轨迹测试方法,其特征在于,步骤(2)中所述对于采集到的加速度信号进行两次积分得到位移信号的具体方法是:
设初始速度和位移分别为v0和s0,加速度传感器测得的信号为a(t),则真实的加速度应为a(t)-a0,理论计算得到的速度v(t)和位移s(t)分别为:
v(t)=v0-a0t+∫a(t)dt
再采用多项式最小二乘法消除趋势项带来的误差影响,设加速度多项式为:
确定的各待定系数cj(j=0,1,…,m),使得与采集到的加速度信号a(t)的误差平方和E最小,即
满足E有极值的条件为:
依次取E对ci求偏导,可以产生一个m+1元线性方程组:
解方程组,求出m+1个待定系数cj(j=0,1,…,m),
上述各式中,m为设定的多项式阶次,其值范围为0≤j≤m,在实际的振动信号处理中,通常取m=1~3来对采样数据进行多项式趋势项消除的处理。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810115565.2A CN108427827A (zh) | 2018-02-05 | 2018-02-05 | 一种基于加速度信号两次积分的输电线舞动轨迹测试方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810115565.2A CN108427827A (zh) | 2018-02-05 | 2018-02-05 | 一种基于加速度信号两次积分的输电线舞动轨迹测试方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108427827A true CN108427827A (zh) | 2018-08-21 |
Family
ID=63156635
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810115565.2A Pending CN108427827A (zh) | 2018-02-05 | 2018-02-05 | 一种基于加速度信号两次积分的输电线舞动轨迹测试方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108427827A (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS58208637A (ja) * | 1982-05-28 | 1983-12-05 | Shimadzu Corp | 送電線のギヤロツピング試験装置 |
CN102279084A (zh) * | 2011-05-03 | 2011-12-14 | 西安工程大学 | 基于微惯性测量组合的输电导线舞动定位系统及定位方法 |
CN105303055A (zh) * | 2015-11-16 | 2016-02-03 | 国家电网公司 | 风振激励作用下特高压输电线路非线性动力学分析方法 |
CN105719071A (zh) * | 2016-01-18 | 2016-06-29 | 中国电力科学研究院 | 一种强风作用下输电线路安全性能的评估方法 |
CN105740549A (zh) * | 2016-02-01 | 2016-07-06 | 西安交通大学 | 220kV输电线路塔线耦联体系风振仿真方法 |
CN107194059A (zh) * | 2017-05-19 | 2017-09-22 | 国网河南省电力公司电力科学研究院 | 一种基于覆冰偏心作用的覆冰导线舞动数值模拟方法 |
-
2018
- 2018-02-05 CN CN201810115565.2A patent/CN108427827A/zh active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS58208637A (ja) * | 1982-05-28 | 1983-12-05 | Shimadzu Corp | 送電線のギヤロツピング試験装置 |
CN102279084A (zh) * | 2011-05-03 | 2011-12-14 | 西安工程大学 | 基于微惯性测量组合的输电导线舞动定位系统及定位方法 |
CN102279084B (zh) * | 2011-05-03 | 2013-08-07 | 西安工程大学 | 基于微惯性测量组合的输电导线舞动定位方法 |
CN105303055A (zh) * | 2015-11-16 | 2016-02-03 | 国家电网公司 | 风振激励作用下特高压输电线路非线性动力学分析方法 |
CN105719071A (zh) * | 2016-01-18 | 2016-06-29 | 中国电力科学研究院 | 一种强风作用下输电线路安全性能的评估方法 |
CN105740549A (zh) * | 2016-02-01 | 2016-07-06 | 西安交通大学 | 220kV输电线路塔线耦联体系风振仿真方法 |
CN107194059A (zh) * | 2017-05-19 | 2017-09-22 | 国网河南省电力公司电力科学研究院 | 一种基于覆冰偏心作用的覆冰导线舞动数值模拟方法 |
Non-Patent Citations (9)
Title |
---|
《力学词典》编辑部: "《力学词典》", 30 June 1990, 中国大百科全书出版社 * |
NILSON BARBIERI ET AL.: "Modal Analysis of Transmission Line Cables", 《DYNAMICS OF COUPLED STRUCTURES》 * |
YU XIAN DI: "Analysis on Galloping Responses of Transmission Tower-Line System", 《APPLIED MECHANICS AND MATERIALS》 * |
周承军: "输电塔线体系动力特性及其影响因素分析", 《安徽建筑》 * |
孙月明,等: "《机械振动学 测试与分析》", 30 November 1991, 浙江大学出版社 * |
张超,等: "输电线舞动轨迹还原理论与试验研究", 《低温建筑技术》 * |
李凯: "输电线路覆冰导线舞动力学特性研究", 《中国优秀硕士学位论文全文数据库(电子期刊)工程科技II辑》 * |
陈为真: "基于时域积分的加速度信号处理", 《华中科技大学学报(自然科学版)》 * |
黄管宝,等: "输电线路导线舞动在线监测系统设计", 《南方电网技术》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109115217B (zh) | 基于电流磁场的输电线路特殊杆塔位置导线参数反演方法 | |
US9057682B2 (en) | Blade vibration measuring apparatus | |
CN107505396A (zh) | 一种结构损伤在线实时监测方法及系统 | |
CN108196166B (zh) | 基于故障波形起始点到达时差的双端行波测距方法 | |
WO2022179247A9 (zh) | 基于全索系索力测试和误差自适应分析的损伤判定方法 | |
CN113343541A (zh) | 长大跨桥梁的涡激振动预警方法、装置及终端 | |
JP2020094992A (ja) | 計測装置及び計測システム | |
WO2010032695A1 (ja) | 画像処理によるパンタグラフの鉛直加速度測定装置および測定方法 | |
CN108427827A (zh) | 一种基于加速度信号两次积分的输电线舞动轨迹测试方法 | |
JPH07301557A (ja) | ギャロッピング検出装置 | |
CN112904158B (zh) | 一种确定gis中局部放电位置的声电联合检测方法 | |
JP2007292700A (ja) | 静止誘導機器の部分放電位置特定方法 | |
JP7304783B2 (ja) | 電力系統監視装置および方法 | |
CN111998735A (zh) | 传感器可任意布阵的超声波报靶方法及装置 | |
CN111521857A (zh) | 基于tmr隧道磁阻的多导体电流测量系统 | |
CN108109394B (zh) | 基于矢量模型的单地磁车辆交通参数检测系统及方法 | |
CN106405203A (zh) | 基于最小二乘法的单相电压跌落检测方法 | |
CN103941224A (zh) | 一种次声源的定向定速方法及系统 | |
JP2000121743A (ja) | 地震動分布の評価方法 | |
CN108375415A (zh) | 一种现场环境下输电线路微风振动监测装置校准方法 | |
JP2021067496A (ja) | 波動振幅統計値計算装置 | |
KR101423152B1 (ko) | 구조물 진동 위치 추적 방법 및 시스템 | |
JP4213017B2 (ja) | 磁性体移動目標の相対位置検出方法 | |
CN117705091B (zh) | 基于大量程石英挠性加速度计的高精度姿态测量方法 | |
CN112763989B (zh) | 一种基于cdif的抖动信号分选方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180821 |