CN103698071A - 基于监测加速度的拉索时变索力历程识别的数据驱动方法 - Google Patents
基于监测加速度的拉索时变索力历程识别的数据驱动方法 Download PDFInfo
- Publication number
- CN103698071A CN103698071A CN201310714346.3A CN201310714346A CN103698071A CN 103698071 A CN103698071 A CN 103698071A CN 201310714346 A CN201310714346 A CN 201310714346A CN 103698071 A CN103698071 A CN 103698071A
- Authority
- CN
- China
- Prior art keywords
- drag
- line
- frequency
- time
- cable
- 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.)
- Granted
Links
Images
Landscapes
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明提出了一种基于监测加速度的拉索时变索力历程识别的数据驱动方法,运用复杂度寻踪这种高效的无监督学习算法,利用拉索上布置的多通道加速度传感器的监测信息,实现了对索力时程进行实时辨识。复杂度寻踪算法可以自主地将拉索的加速度响应分解为拉索的单模态响应,进而通过极短时间的加速度信息辨识拉索的实时频率,通过张紧弦理论计算索力时程。通过结合实测索力和实测风速的斜拉桥的模拟分析以及拉索的模型试验,证明了所提出的复杂度寻踪算法能够对时变索力历程进行准确的实时识别。本发明是一种直接有效的时变索力历程辨识方法,简单易用,索力辨识精度高,时效性强且能够实现在线实时辨识,尤其适用于拉索的在线评估。
Description
技术领域
本发明涉及一种土木工程结构响应辨识的方法,具体涉及一种基于监测加速度的拉索时变索力历程识别的数据驱动方法。
背景技术
斜拉桥由于其跨度大、造型美观、施工方便等特点,是目前世界上应用最广泛的桥型。斜拉桥一般会成为一个地区的标志性建筑,而且往往是一个地区的交通枢纽,对该地区的政治经济有着至关重要的影响。
斜拉索作为斜拉桥主要的承重构件,由高强度钢丝束和PE护套组成。斜拉索在长达几十年使用期内,环境侵蚀、材料老化和荷载的长期效应、疲劳效应与突变效应等灾害因素的耦合作用下,将不可避免地导致结构和***的损伤积累和抗力衰减,从而抵抗自然灾害、甚至正常环境作用的能力下降,极端情况下将引发灾难性的突发事故。为保证桥梁结构的正常运营,最好的解决方法就是对斜拉索进行定期的检测维护,且需要依据检测结果对其进行安全评估。但是由于无损检测方法(漏磁检测、X射线检测、超声检测、基于振动法的索力检测等)的局限性,在斜拉桥实际运营过程中,只有极小部分斜拉索都得到检测。鉴于斜拉索病害的普遍性,以及斜拉索修复和更换的高昂费用,为应对斜拉索结构运营状况下进行腐蚀疲劳评估的需要,目前亟待认识斜拉索结构在服役全过程中的结构索力历程响应。
目前存在的索力测试装置有:压力环、磁通量传感器、光纤光栅智能拉索等,其中压力环和光纤光栅智能拉索能够直接测试索力历程,广泛的安装和应用于新建桥梁上。但是,这些索力监测传感器一般价格昂贵,安装复杂(只能用于新建桥梁),更主要的是传感器的耐久性比较差,这些固有缺点限制了上述传感器的大规模应用。由于更换索力监测传感器耗时耗力,价格昂贵,因而亟需发展一种省时、省力、经济的实时索力历程监测方法。
发明内容
基于以上不足之处,本发明提供一种基于多通道监测加速度辨识时变拉索索力历程的数据驱动方法,解决拉索时变索力历程识别的问题。
本发明采用如下技术方案实现:一种基于监测加速度的拉索时变索力历程识别的数据驱动方法,步骤如下:
步骤1:对待测试拉索同一平面内布设2-3个的加速度传感器,测试斜拉索在环境激励下的加速度响应,并通过一个通道的10秒内的加速度信号,辨识拉索基频及该时间段内能够激励起来的最高频率,拉索的基频用f1表示,监测加速度信号中能够辨识的拉索的最高频率用fi表示,fi≈i×f1,其中,表示拉索在索力T作用下,第i阶模态的圆频率,μ为单位索长的密度,L为斜拉索长度为L;
步骤2:设计高通滤波器,截止频率(i-1.5)f1,将拉索的多通道加速度信号中出现的最高两阶频率fi和fi-1提取出来;
步骤3:选择3秒的窗函数,对窗函数内的多通道加速度信号进行滤波预处理,预处理后加速度信号通过复杂度寻踪算法分离得到单模态响应信号,通过快速傅里叶变换或功率谱计算单模态响应信号的频率,利用张紧弦理论,通过辨识的频率计算拉索的时变索力;
步骤4:滑动窗函数,对窗函数内的加速度信号重复步骤3,辨识拉索的时变索力历程。
本发明还具有如下技术特征:
1、应通过多通道加速度信号辨识拉索的时变频率,进而辨识拉索的时变索力历程。
2、对多通道加速度信号进行高通滤波预处理,提取拉索能够被激励起的最高的2阶频率。
3、利用复杂度寻踪处理算法,将滑动窗内多模态响应分解为单模态响应,频率分辨率Δf≤0.025HZ的频率分辨率准确辨识单模态响应的频率,利用张紧弦理论计算时变索力历程。
本发明基于成熟的加速度传感器测试技术,相对于目前发展的其他类型的索力监测传感器,本发明采用的加速度传感器成熟可靠,测试精度高,价格便宜,传感器安装以及更换都很方便,从而使得本发明提出的拉索时变索力历程辨识***具有极高的可靠性和耐久性。本发明监测加速度信息,通过辨识时变频率计算索力,方法简单易用,索力辨识精度高,时效性强且能够实现在线实时辨识,本发明方法的鲁棒性及可靠性强,尤其适用于拉索的在线评估。
附图说明
图1为Benchmark桥梁的智能拉索制作构造和C8拉索位置简图;
图2为Benchmark桥梁实测风速和索力时程曲线图(测试时间为2008年1月17日1:00-2:00AM,红框内为选择的待辨识索力);
图3为Benchmark桥梁拉索加速度时程曲线图(开始时间为2008年1月17日1:40AM);
图4为C8拉索八分点通道10秒加速度功率谱图;
图5为C8拉索倍频特性图;
图6为C8拉索八分点通道和四分点通道3秒加速度滤波后的时程曲线及其功率谱图;
图7为复杂度寻踪得到的单模态响应时程曲线及其功率谱图;
图8为辨识的C8拉索时变频率图;
图9为本发明算法辨识C8拉索索力与实测索力对比图;
图10为拉索时变索力历程辨识试验装置简图;
图11为S#1测试加速度10秒功率谱图;
图12为试验拉索倍频关系图
图13为试验拉索S#1和S#3传感器3秒加速度滤波后的时程曲线及其功率谱;
图14为复杂度寻踪得到的单模态响应时程曲线及其功率谱图;
图15为辨识的试验拉索时变频率图;
图16为本发明算法辨识试验拉索索力与实测索力对比工况1图;
图17为本发明算法辨识试验拉索索力与实测索力对比工况2图。
具体实施方式
本发明的具体实施方案,通过结合实测索力和实测风速的斜拉桥的模拟分析以及拉索的模型试验进行说明。
对于小垂度的斜拉索(长度为L,截面面积为A,杨氏模量为E),其平面外的第i阶模态的横向振动运动方程可以表示为
式中,mi=μL/2表示单位索长的重量,μ表示单位索长的密度;ζi为第i阶模态阻尼比;αi=i2π2T/2L,βik=EAπ4i2k2/8L3;Fi(t)为外界激励;表示拉索在索力T作用下,第i阶模态的圆频率,ωi的单位是rad/s。的计算公式表明,拉索各阶圆频率与基频之间存在倍频关系,即ωi≈i×ω1;同时拉索的各阶模态圆频率ωi都和拉索索力T之间存在直接的映射关系。依据张紧弦理论,拉索索力T可以表示为
本发明所述的基于监测加速度的拉索时变索力历程识别的数据驱动方法的核心即是时变模态圆频率ωi的辨识算法。依据模态叠加理论,环境激励下拉索的位移向量x(t)=[x1(t),…,xn(t)]T可以表示为
式中,表示模态振型矩阵,q(t)=[q1(t),…,qn(t)]T为模态响应向量,第i阶振型(Φ矩阵的第i列)和模态响应分别表示为φi和qi(t)。上述位移向量的分解计算通过盲元分离算法实现,令表示模态响应向量估计值,表示模态振型向量估计值,yi(t)信号的复杂度定义为
本发明中, hL>>hS,和矩阵可以通过快速卷积算法计算获得。如果信号yi(t)为理想的单模态响应信号,信号复杂度函数值F(yi)最大,因此只需要对信号复杂度函数F(yi)求极值,最优响应信号yi(t)就能准确描述模态响应向量的重构值。
复杂度函数F(yi)求极值的过程可以通过经典的梯度优化算法实现,复杂度函数F(yi)对wi的微分可以表示为
上式变形为一个广义特征值问题,所以wi即为的特征向量,特征值为γi=Vi/Ui。获得wi后即可有公式(3)得到重构的单模态响应向量上述通过优化信号复杂度函数F(yi),通过监测位移响应向量x(t)获得为单模态响应q(t)=[q1(t),…,qn(t)]T的过程称为复杂度寻踪算法。对于加速度响应和速度响应可以通过同样的方式进行计算。
对于单模态响应信号q(t)=[q1(t),…,qn(t)]T,由于每个模态响应中只有控制模态的频率成分占优,所以只需要对分离出的q(t)=[q1(t),…,qn(t)]T进行快速傅里叶(FFT)变换或功率谱分析,即可求得信号的模态频率,通过公式(2)所示拉索索力与拉索频率之间的关系,求得拉索索力。
对于单模态响应信号qi(t)而言,频率分辨率Δf为选择的窗函数时间τ的倒数Δf=1/τ,存在误差项的索力计算公式为
实施例1:
Benchmark斜拉桥(图1上的C8斜拉索),由139根5mm钢丝组成,索长L=100.95m,截面面积A=2.73×10-3m2,单位长度的密度为μ=21.43kg/m。2006年,C8斜拉索在换索工程中,更换为智能拉索,能够实时监测索力历程。
步骤1:本发明算法需要利用2个或最多3个通道的加速度传感器信息辨识拉索索力,通过数值模拟计算C8斜拉索在监测风速(图2(a))和监测索力(图2(b))作用下L/8和L/4位置的加速度响应。对C8拉索L/8通道的10秒加速度时程(图3)进行功率谱分析(图4),C8拉索激励起来的最高阶频率为第30阶,同时可以看出拉索频率之间存在明显的倍频关系(图5)。
步骤2:设计高通滤波器,截止频率约等于(30-1.5)f1≈38Hz,将第29阶和30阶频率成分过滤出来;
步骤3:选择3秒的窗函数,对窗函数内的加速度曲线应用步骤2设计的滤波器进行滤波,滤波后的L/8和L/4通道的加速度时程曲线及其功率谱见图6。对滤波后的L/8和L/4通道的加速度时程曲线进行复杂度寻踪算法处理,分离出单模态响应信号,并计算其功率谱(图7),辨识单模态信号的占主导模态的频率,并利用张紧弦理论计算索力。
步骤4:逐点滑动3秒的时间窗,重复步骤3,对滑动窗内加速度信号滤波,分离单模态响应信号,通过傅里叶变换或者功率谱分析辨识频率,从而得到频率随时间变化的曲线(图8),进而得到索力时程曲线(图9)。
监测索力时程曲线与辨识索力时程曲线很好的吻合在一起,验证了本发明所提算法的准确性。
实施例2:
试验拉索长14.02m,直径为1.5cm,单位长度质量为1.33kg/m。试验拉索一段固定,另一端采用螺纹杆调节索的拉伸长度,进而调节索力的大小,同时在试验拉索端部布置测力计,监测索力时程的变化,试验拉索采用2个风机作为外界激励产生振动。
步骤1:分别在L/4、L/2和3/4L布设加速度传感器S#1-S#3测试平面外响应,采样频率为200Hz,试验装置、传感器布置以几何尺寸如图8所示。对试验拉索S#1传感器的10秒加速度时程进行功率谱分析(图11),试验拉索激励起来的最高阶频率为第14阶,同时可以看出试验拉索频率之间存在明显的倍频关系(图12);
步骤2:设计高通滤波器,截止频率约等于(14-1.5)f1≈28Hz,能够将第13阶和14阶频率成分过滤出来;
步骤3:选择3秒的窗函数,对窗函数内的加速度曲线应用步骤2设计的滤波器进行滤波,滤波后的S#1和S#3传感器的加速度时程曲线及其功率谱见图13。对滤波后的S#1和S#3传感器的加速度时程曲线进行复杂度寻踪算法处理,分离出单模态响应信号,并计算其功率谱(图14),辨识单模态信号的占主导模态的频率,并利用张紧弦理论计算索力。
步骤4:逐点滑动3秒的时间窗,重复步骤3,得到频率随时间变化的曲线(图15),进而得到索力时程曲线(图16-17)。试验拉索2个工况的监测索力时程曲线与辨识索力时程曲线很好的吻合在一起,从试验角度验证了本发明所提算法的准确性。
Claims (4)
1.一种基于监测加速度的拉索时变索力历程识别的数据驱动方法,其特征在于,步骤如下:
步骤1:对待测试拉索同一平面内布设2-3个的加速度传感器,测试斜拉索在环境激励下的加速度响应,并通过一个通道的10秒内的加速度信号,辨识拉索基频及该时间段内能够激励起来的最高频率,拉索的基频用f1表示,监测加速度信号中能够辨识的拉索的最高频率用fi表示,fi≈i×f1,其中,表示拉索在索力T作用下,第i阶模态的圆频率,μ为单位索长的密度,L为斜拉索长度为L;
步骤2:设计高通滤波器,截止频率(i-1.5)f1,将拉索的多通道加速度信号中出现的最高两阶频率fi和fi-1提取出来;
步骤3:选择3秒的窗函数,对窗函数内的多通道加速度信号进行滤波预处理,预处理后加速度信号通过复杂度寻踪算法分离得到单模态响应信号,通过快速傅里叶变换或功率谱计算单模态响应信号的频率,利用张紧弦理论,通过辨识的频率计算拉索的时变索力;
步骤4:滑动窗函数,对窗函数内的加速度信号重复步骤3,辨识拉索的时变索力历程。
2.根据权利要求1所述的基于监测加速度的拉索时变索力历程识别的数据驱动方法,其特征在于:应通过多通道加速度信号辨识拉索的时变频率,进而辨识拉索的时变索力历程。
3.根据权利要求1所述的基于监测加速度的拉索时变索力历程识别的数据驱动方法,其特征在于:对多通道加速度信号进行高通滤波预处理,提取拉索能够被激励起的最高的2阶频率。
4.根据权利要求1所述的基于监测加速度的拉索时变索力历程识别的数据驱动方法,其特征在于:利用复杂度寻踪处理算法,将滑动窗内多模态响应分解为单模态响应,频率分辨率Δf≤0.025HZ的频率分辨率准确辨识单模态响应的频率,利用张紧弦理论计算时变索力历程。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310714346.3A CN103698071B (zh) | 2013-12-12 | 2013-12-12 | 基于监测加速度的拉索时变索力历程识别的数据驱动方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310714346.3A CN103698071B (zh) | 2013-12-12 | 2013-12-12 | 基于监测加速度的拉索时变索力历程识别的数据驱动方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103698071A true CN103698071A (zh) | 2014-04-02 |
CN103698071B CN103698071B (zh) | 2015-11-25 |
Family
ID=50359676
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310714346.3A Expired - Fee Related CN103698071B (zh) | 2013-12-12 | 2013-12-12 | 基于监测加速度的拉索时变索力历程识别的数据驱动方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103698071B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104713673A (zh) * | 2015-03-11 | 2015-06-17 | 哈尔滨工业大学 | 基于扩展卡尔曼滤波器的拉索时变索力历程识别方法 |
CN108613763A (zh) * | 2018-05-11 | 2018-10-02 | 上海市建筑科学研究院 | 基于频率法的光纤光栅式拉索索力监测传感器及传感方法 |
CN111242058A (zh) * | 2020-01-16 | 2020-06-05 | 江苏东印智慧工程技术研究院有限公司 | 一种基于Android与HHT的时变索力识别方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101051226A (zh) * | 2007-02-27 | 2007-10-10 | 宁波杉工结构监测与控制工程中心有限公司 | 斜拉索索力远程在线监测方法及设备 |
KR20080047186A (ko) * | 2006-11-24 | 2008-05-28 | 부산대학교 산학협력단 | 행어 케이블의 장력 측정방법 |
CN101201282A (zh) * | 2007-12-20 | 2008-06-18 | 宁波大学 | 一种用于斜拉桥索力检测的基频识别方法 |
CN101368860A (zh) * | 2008-09-12 | 2009-02-18 | 江苏工业学院 | 频率法检测斜拉桥斜拉索索力中fft数据的修正方法 |
CN101762346A (zh) * | 2009-12-31 | 2010-06-30 | 北京市建筑工程研究院 | 一种用多频率法测量多跨钢拉索索力的方法 |
CN101936795A (zh) * | 2010-07-27 | 2011-01-05 | 同济大学 | 基于模态分析高精度的拉索索力测试方法 |
-
2013
- 2013-12-12 CN CN201310714346.3A patent/CN103698071B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20080047186A (ko) * | 2006-11-24 | 2008-05-28 | 부산대학교 산학협력단 | 행어 케이블의 장력 측정방법 |
CN101051226A (zh) * | 2007-02-27 | 2007-10-10 | 宁波杉工结构监测与控制工程中心有限公司 | 斜拉索索力远程在线监测方法及设备 |
CN101201282A (zh) * | 2007-12-20 | 2008-06-18 | 宁波大学 | 一种用于斜拉桥索力检测的基频识别方法 |
CN101368860A (zh) * | 2008-09-12 | 2009-02-18 | 江苏工业学院 | 频率法检测斜拉桥斜拉索索力中fft数据的修正方法 |
CN101762346A (zh) * | 2009-12-31 | 2010-06-30 | 北京市建筑工程研究院 | 一种用多频率法测量多跨钢拉索索力的方法 |
CN101936795A (zh) * | 2010-07-27 | 2011-01-05 | 同济大学 | 基于模态分析高精度的拉索索力测试方法 |
Non-Patent Citations (1)
Title |
---|
陈得民, 储伟伟, 王亚涛: "一种新型桥梁索力测试***及方法", 《传感器世界》, no. 3, 31 March 2013 (2013-03-31), pages 21 - 25 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104713673A (zh) * | 2015-03-11 | 2015-06-17 | 哈尔滨工业大学 | 基于扩展卡尔曼滤波器的拉索时变索力历程识别方法 |
CN108613763A (zh) * | 2018-05-11 | 2018-10-02 | 上海市建筑科学研究院 | 基于频率法的光纤光栅式拉索索力监测传感器及传感方法 |
CN108613763B (zh) * | 2018-05-11 | 2023-11-03 | 上海市建筑科学研究院 | 基于频率法的光纤光栅式拉索索力监测传感器及传感方法 |
CN111242058A (zh) * | 2020-01-16 | 2020-06-05 | 江苏东印智慧工程技术研究院有限公司 | 一种基于Android与HHT的时变索力识别方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103698071B (zh) | 2015-11-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101936795B (zh) | 基于模态分析高精度的拉索索力测试方法 | |
CN101504324B (zh) | 拉索运营状态智能判别方法及*** | |
JP5514152B2 (ja) | 構造物安全性の分析方法 | |
CN106092623B (zh) | 一种基于长标距刚度系数的桥梁损伤识别评估方法 | |
CN102353509B (zh) | 一种基于分块冲击振动测试的中小型桥梁快速检测方法 | |
CN105716814B (zh) | 一种评估桁架结构损伤的实时监测***及其方法 | |
CN104331595B (zh) | 桥梁损伤预警的移动主成分相关系数法 | |
CN108444662B (zh) | 一种基于日温度效应的桥梁损伤在线监测方法 | |
CN104198144B (zh) | 一种基于长标距光纤应变传感器的中小桥梁快速检测方法 | |
CN105241660A (zh) | 基于健康监测数据的高铁大型桥梁性能评定方法 | |
CN101762347B (zh) | 一种用半波法测量多跨钢拉索索力的方法 | |
CN112461358B (zh) | 基于车辆-桥梁***瞬时频率的桥梁模态参数识别方法 | |
CN110285909B (zh) | 基于同步压缩变换的索承桥梁瞬时索力计算方法 | |
CN101839781B (zh) | 悬索桥主缆状态快速识别方法及装置 | |
CN101266190A (zh) | 正常交通流下斜拉桥模态参数测定装置及方法 | |
CN102254068A (zh) | 大跨度桥梁抖振响应的多尺度分析方法 | |
CN103018006A (zh) | 船模波浪载荷试验用测量梁的弯矩标定装置及方法 | |
CN104502010A (zh) | 未知边界条件下索力测试方法及装置 | |
CN103698071B (zh) | 基于监测加速度的拉索时变索力历程识别的数据驱动方法 | |
CN102121858A (zh) | 部分斜拉桥拉索钢绞线张力测试方法 | |
CN112326787B (zh) | 一种基于专属测试车多点快速静态采集的梁式桥识别方法 | |
CN103913512B (zh) | 斜拉索定期检测的损伤定位*** | |
CN105912504A (zh) | 一种机械结构多源冲击载荷识别的压缩感知方法 | |
CN105823591A (zh) | 一种含复杂边界的在役拱桥短吊杆张力识别方法 | |
CN111753776B (zh) | 基于回声状态与多尺度卷积联合模型的结构损伤识别方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20151125 Termination date: 20161212 |
|
CF01 | Termination of patent right due to non-payment of annual fee |