CN107766630A - 高速铁路桥梁支座在设计使用寿命内累积磨损的评估方法 - Google Patents
高速铁路桥梁支座在设计使用寿命内累积磨损的评估方法 Download PDFInfo
- Publication number
- CN107766630A CN107766630A CN201710932238.1A CN201710932238A CN107766630A CN 107766630 A CN107766630 A CN 107766630A CN 201710932238 A CN201710932238 A CN 201710932238A CN 107766630 A CN107766630 A CN 107766630A
- Authority
- CN
- China
- Prior art keywords
- value
- cumulative
- longitudinal
- displacement
- dynamic displacement
- 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
- 230000001186 cumulative effect Effects 0.000 title claims abstract description 192
- 238000013461 design Methods 0.000 title claims abstract description 32
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000006073 displacement reaction Methods 0.000 claims abstract description 258
- 238000005315 distribution function Methods 0.000 claims abstract description 77
- 238000012544 monitoring process Methods 0.000 claims abstract description 71
- 230000003203 everyday effect Effects 0.000 claims abstract description 23
- 238000004088 simulation Methods 0.000 claims description 38
- 238000000354 decomposition reaction Methods 0.000 claims description 25
- 238000005070 sampling Methods 0.000 claims description 19
- 238000011156 evaluation Methods 0.000 claims description 16
- 238000012614 Monte-Carlo sampling Methods 0.000 claims description 15
- 238000007619 statistical method Methods 0.000 claims description 11
- 238000012545 processing Methods 0.000 claims description 10
- 238000005516 engineering process Methods 0.000 claims description 5
- 230000002354 daily effect Effects 0.000 claims description 3
- 230000007613 environmental effect Effects 0.000 description 7
- 238000004458 analytical method Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 230000036541 health Effects 0.000 description 3
- 238000012423 maintenance Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000003628 erosive effect Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000004643 material aging Methods 0.000 description 1
- 238000012067 mathematical method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/04—Ageing analysis or optimisation against ageing
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Theoretical Computer Science (AREA)
- Civil Engineering (AREA)
- Structural Engineering (AREA)
- Computational Mathematics (AREA)
- Architecture (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明公开了一种高速铁路桥梁支座在设计使用寿命内累积磨损的评估方法,步骤10)对支座纵向位移进行监测,采集纵向位移监测数据;步骤20)提取纵向动位移值;步骤30)获取纵向动位移累积值;步骤40)将每天的纵向动位移累积值划分为两类;步骤50)将两类纵向动位移累积值视为随机变量,并确定纵向动位移累积值威布尔分布函数;步骤60)确定支座在设计使用寿命内的支座累积行程模拟值;步骤70)确定支座累积行程模拟值的广义极值分布函数;步骤80)计算支座累积行程模拟值超过支座磨损上限值的失效概率,比较失效概率与预设概率阈值。本发明可更准确地对桥梁支座在设计使用寿命内的累积磨损值进行安全评估。
Description
技术领域
本发明属于铁路桥梁结构健康监测与安全评估领域,具体来说,涉及一种高速铁路桥梁支座在设计使用寿命内累积磨损的评估方法。
背景技术
目前,中国高速铁路网飞速发展。高速铁路桥梁多具有刚度较大、整体性好、减振性能和横向位移控制要求高等特点,高速铁路桥梁在服役期间的运营安全问题备受瞩目。随着服役时间增长,高速铁路桥梁结构会受到环境侵蚀、材料老化、列车荷载等多重因素影响,其安全性和耐久性会不可避免地发生退化,这必然影响到高速铁路桥梁结构的安全运营。与此同时,高速铁路桥梁结构中的一些关键性构件在复杂环境和各类荷载的共同作用下,很容易产生损伤。因此,针对关键性结构构件进行实时评估与安全预警显得尤为重要。支座是大跨铁路桥梁结构中最重要的构件之一,在实际服役运营中,支座会出现支座变形和磨损、橡胶圈液化和承载力下降等严重问题,将会直接影响到铁路桥梁结构的使用寿命,并危及桥梁结构的整体安全性。因此,对高速铁路桥梁支座在设计使用寿命内的累积磨损程度进行安全评估具有重要意义。
近年来桥梁健康监测系统已成为桥梁管理及养护的重要辅助手段。基于桥梁健康监测手段,可以直接监测到高速铁路桥梁支座的纵向动位移累积行程。但目前尚未提出利用纵向动位移累积行程监测数据对支座在服役期内的磨损程度进行安全评估的方法。环境荷载引起的支座纵向动位移与列车荷载引起的支座纵向动位移具有完全不同的随机特性。目前,研究尚未将这两种随机特性进行区别对待。此外,监测数据往往不能在支座服役期内持续采集,达不到对支座在服役期内的磨损程度进行安全评估所需的足够样本数量要求。
发明内容
技术问题:本发明所要解决的技术问题是:提供一种高速铁路桥梁支座在设计使用寿命内累积磨损的评估方法,能够考虑支座纵向动位移的双重随机特性,并能有效解决纵向动位移监测数据量过少的问题,使得支座在设计使用寿命内的累积磨损的评估更加准确。
技术方案:为解决上述技术问题,本发明采用以下的技术方案:
一种高速铁路桥梁支座在设计使用寿命内累积磨损的评估方法,包括如下步骤:
步骤10):在桥梁主梁的支座上安装纵向位移监测传感器,对支座纵向位移进行监测,并采集纵向位移监测数据,其中,第t时刻支座纵向位移值记为v(t),t表示时间,且t=1/f、2/f、…、L/f,f表示采样频率,L表示采样总数,L=86400df,d表示采样天数;
步骤20):利用小波包分解法提取纵向位移监测数据中的纵向动位移值;
步骤30):对所述纵向动位移值作一阶差分处理,并取绝对值累加,得到纵向动位移累积值;
步骤40):将每天的纵向动位移累积值划分为两类,构建第一纵向动位移累积值序列和第二纵向动位移累积值序列;
步骤50):将两类纵向动位移累积值视为随机变量,并确定纵向动位移累积值威布尔分布函数;
步骤60):利用蒙特卡洛抽样方法确定支座在设计使用寿命内的支座累积行程模拟值;
步骤70):将支座累积行程模拟值视为随机变量,并确定支座累积行程模拟值的广义极值分布函数;
步骤80):计算支座累积行程模拟值超过支座磨损上限值的失效概率,比较失效概率与预设概率阈值,如果失效概率大于或等于预设概率阈值,则支座达到磨损上限,更换支座。
作为优选例,所述步骤20)包括:利用小波包分解技术,对支座纵向位移监测数据进行8层尺度的小波包分解,得到32个小波包分解系数,剔除第一个小波包分解系数后,对剩余小波包分解系数重构并叠加,得到支座纵向动位移值。
作为优选例,所述步骤30)具体包括:将每天凌晨0时至6时的纵向动位移监测数据划分为n1份,其中第i天中第j份纵向动位移监测数据采用vi,j表示,其中i=1、2、…、d,j=1、2、…、n1,对每份纵向动位移值vi,j作一阶差分处理,并取绝对值累加,得到第i天中第j份纵向动位移监测数据的动位移累积值Ci,j,如式(1)所示:
式中,Ni,j表示第i天中第j份纵向动位移监测数据中的数据总数,vi,j,p+1表示第i天中第j份纵向动位移监测数据中的第p+1个值,vi,j,p表示第i天中第j份纵向动位移监测数据中的第p个值;
将每天凌晨6时至夜晚24时的纵向动位移监测数据划分为n2份,其中第m天中第n份纵向动位移监测数据采用vm,n表示,其中m=1、2、…、d,n=1、2、…、n2,对每份纵向动位移值vm,n作一阶差分处理,并取绝对值累加,得到第m天中第n份纵向动位移监测数据的动位移累积值Dm,n,如式(2)所示:
式中,Nm,n表示第m天中第n份纵向动位移监测数据中的数据总数,vm,n,q+1表示第m天中第n份纵向动位移监测数据中的第q+1个值,vm,n,q表示第m天中第n份纵向动位移监测数据中的第q个值。
作为优选例,所述步骤40)具体包括:将所有天数在每天凌晨0时至6时的纵向动位移累积值组合构成第一纵向动位移累积值序列C,C=[C1,1、C1,2、…、C1,n1、C2,1、C2,2、…、C2,n1、…、Cd,1、Cd,2、…、Cd,n1];
将所有天数在每天凌晨6时至24时的纵向动位移累积值组合构成第二纵向动位移累积值序列D,D=[D1,1、D1,2、…、D1,n2、D2,1、D2,2、…、D2,n2、…、Dd,1、Dd,2、…、Dd,n2]。
作为优选例,所述步骤50)具体包括:将第一纵向动位移累积值序列C中的纵向动位移累积值视为统计变量,利用统计学方法计算统计变量的概率密度值,再利用威布尔分布函数对统计变量概率密度值进行拟合,得到威布尔分布函数的第一形状参数k1和第一比例参数λ1的拟合值,从而确定第一纵向动位移累积值序列中纵向动位移累积值的威布尔分布函数,如式(3)所示:
式中,g(C)表示第一纵向动位移累积值序列中纵向动位移累积值的威布尔分布函数,k1表示威布尔分布函数的第一形状参数,λ1表示威布尔分布函数的第一比例参数,e表示指数函数的底,e=2.71828;
将第二纵向动位移累积值序列D中的纵向动位移累积值视为统计变量,利用统计学方法计算统计变量概率密度值,在利用威布尔分布函数对统计变量概率密度值进行拟合,得到威布尔分布函数的第二形状参数k2和第二比例参数λ2的拟合值,从而确定第二纵向动位移累积值序列中纵向动位移累积值的威布尔分布函数,如式(4)所示:
式中,g(D)表示第二纵向动位移累积值序列D中纵向动位移累积值的威布尔分布函数,k2表示威布尔分布函数的第二形状参数,λ2表示威布尔分布的第二比例参数。
作为优选例,所述步骤60)具体包括:设支座的设计使用寿命为R年,利用第一纵向动位移累积值序列中纵向动位移累积值的威布尔分布函数g(C)对支座在R年内每天凌晨0时至6时的所有纵向动位移累积值进行蒙特卡洛抽样模拟,模拟总数NC=365Rn1;利用第二纵向动位移累积值序列D中纵向动位移累积值的威布尔分布函数g(D)对支座在R年内每天凌晨6时至24时的所有纵向动位移累积值进行蒙特卡洛抽样模拟,模拟总数ND=365Rn2;
在R年内共得到31536000Rf个模拟的纵向动位移累积值,即NC+ND=31536000Rf,将31536000Rf个纵向动位移累积模拟值进行求和,得到支座在设计使用寿命内的累积行程模拟值。
作为优选例,所述步骤70)具体包括:将步骤50)和步骤60)重复操作100次,得到100个累积行程模拟值,将累积行程模拟值视为统计变量,并采用H表示;利用统计学方法计算所述100个累积行程模拟值的广义极值分布值,利用广义极值分布函数对累积行程模拟值广义极值分布值进行拟合,得到广义极值分布函数的形状参数r3、比例参数a3和位置参数b3的拟合值,从而确定支座累积行程模拟值广义极值分布函数,如式(5)所示:
式中,t(H)表示累积行程模拟值的广义极值分布函数,r3表示广义极值分布函数的形状参数,a3表示广义极值分布函数的比例参数,b3表示广义极值分布函数的位置参数。
作为优选例,所述步骤80)具体包括:设支座磨损上限值为[H],利用式(6)计算支座累积行程模拟值超过支座磨损上限值的失效概率β:
其中,d表示采样天数;
若β大于或等于预设概率阈值,则认为支座已经达到磨损上限,需要及时更换。
作为优选例,所述预设概率阈值为5%。
作为优选例,所述步骤10)中,采样频率f≥1Hz,采样天数d≥1;所述步骤20)中,n1的取值在30~360之间,n2的取值在90~1080之间。
有益效果:与现有技术相比,本发明具有以下优点:本发明实施例的方法评估准确性更高。环境荷载引起的支座纵向动位移与列车荷载引起的支座纵向动位移具有完全不同的随机特性。本发明实施例步骤30)将每天内的纵向动位移累积值划分为两类:由随机环境引起的纵向动位移累积值、随机环境与列车荷载共同引起的纵向动位移累积值,区别这两种随机特性,从而更准确地考虑了支座纵向动位移的两种随机特性,以提高评估准确性。另外,现有技术中,不能在支座服役期内持续采集纵向动位移监测数据,从而达不到对支座在服役期内的磨损程度进行安全评估所需的足够样本数量要求。本发明利用步骤40)首先确定了纵向动位移累积值的概率密度函数,然后利用步骤50)通过蒙特卡洛抽样模拟得到了支座在设计使用寿命内的支座累积行程模拟值,从而满足了样本数量要求,使得评估准确性更高。
附图说明
图1为本发明实施例中对支座采集的纵向位移监测数据示意图;
图2为本发明实施例中对支座采集的纵向动位移监测数据示意图;
图3为本发明实施例中纵向动位移累积值序列C的变化趋势;
图4为本发明实施例中纵向动位移累积值序列D的变化趋势;
图5为本发明实施例中威布尔分布函数g(C)的拟合曲线;
图6为本发明实施例中威布尔分布函数g(D)的拟合曲线;
图7为本发明实施例中100个累积行程模拟值的变化趋势;
图8为本发明实施例中广义极值分布函数t(H)的拟合曲线。
具体实施方式
下面将参照附图,对本发明的技术方案进行详细的说明。
本发明实施例的一种高速铁路桥梁支座在设计使用寿命内累积磨损的评估方法,包括如下步骤:
步骤10):在桥梁主梁的支座上安装纵向位移监测传感器,对支座纵向位移进行监测,并采集纵向位移监测数据,其中,第t时刻支座纵向位移值记为v(t),t表示时间,且t=1/f、2/f、…、L/f,f表示采样频率,L表示采样总数,L=86400df,d表示采样天数。
步骤10)中,采样频率f≥1Hz,采样天数d≥1;所述步骤20)中,n1的取值在30~360之间,n2的取值在90~1080之间。
步骤20):利用小波包分解法提取纵向位移监测数据中的纵向动位移值。
步骤20)具体包括:利用小波包分解技术,对支座纵向位移监测数据进行8层尺度的小波包分解,得到32个小波包分解系数,剔除第一个小波包分解系数后,对剩余小波包分解系数重构并叠加,得到支座纵向动位移值。小波包分解法是一种精细分析方法,能够根据被分析信号的特征,自适应地选取最佳基函数,使之与信号相匹配,以提高信号的分析能力。
步骤30):对所述纵向动位移值作一阶差分处理,并取绝对值累加,得到纵向动位移累积值。
步骤30)具体包括:将每天凌晨0时至6时的纵向动位移监测数据划分为n1份,其中第i天中第j份纵向动位移监测数据采用vi,j表示,其中i=1、2、…、d,j=1、2、…、n1。作为优选,n1的取值在30~360之间,可使式(3)中威布尔分布函数的拟合效果更好。对每份纵向动位移值vi,j作一阶差分处理,并取绝对值累加,得到第i天中第j份纵向动位移监测数据的动位移累积值Ci,j,如式(1)所示:
式中,Ni,j表示第i天中第j份纵向动位移监测数据中的数据总数,vi,j,p+1表示第i天中第j份纵向动位移监测数据中的第p+1个值,vi,j,p表示第i天中第j份纵向动位移监测数据中的第p个值;
将每天凌晨6时至夜晚24时的纵向动位移监测数据划分为n2份,其中第m天中第n份纵向动位移监测数据采用vm,n表示,其中m=1、2、…、d,n=1、2、…、n2。作为优选,n2的取值在90~1080之间,可使式(4)中威布尔分布函数的拟合效果更好。对每份纵向动位移值vm,n作一阶差分处理,并取绝对值累加,得到第m天中第n份纵向动位移监测数据的动位移累积值Dm,n,如式(2)所示:
式中,Nm,n表示第m天中第n份纵向动位移监测数据中的数据总数,vm,n,q+1表示第m天中第n份纵向动位移监测数据中的第q+1个值,vm,n,q表示第m天中第n份纵向动位移监测数据中的第q个值。
步骤40):将每天的纵向动位移累积值划分为两类,构建第一纵向动位移累积值序列和第二纵向动位移累积值序列。将每天内的纵向动位移累积值划分为两类:由随机环境引起的纵向动位移累积值,以及随机环境与列车荷载共同引起的纵向动位移累积值。在每天凌晨0时至6时的纵向动位移累积值主要由随机环境荷载引起,作为第一纵向动位移累积值序列。在每天凌晨6时至夜晚24时的纵向动位移累积值主要由随机环境荷载和列车荷载共同引起,作为第二纵向动位移累积值序列。
步骤40)具体包括:将所有天数在每天凌晨0时至6时的纵向动位移累积值组合构成第一纵向动位移累积值序列C,C=[C1,1、C1,2、…、C1,n1、C2,1、C2,2、…、C2,n1、…、Cd,1、Cd,2、…、Cd,n1];
将所有天数在每天凌晨6时至24时的纵向动位移累积值组合构成第二纵向动位移累积值序列D,D=[D1,1、D1,2、…、D1,n2、D2,1、D2,2、…、D2,n2、…、Dd,1、Dd,2、…、Dd,n2]。
步骤50):将两类纵向动位移累积值视为随机变量,并确定纵向动位移累积值威布尔分布函数。
步骤50)具体包括:将第一纵向动位移累积值序列C中的纵向动位移累积值视为统计变量,利用统计学方法计算统计变量的概率密度值。在统计学中,概率密度分析是在统计学中最常用的一种概率密度分析方法。概率密度等于一段区间的概率除以该段区间的长度。它的值是非负的,概率密度对区间的积分就是事件在这个区间发生的概率,概率密度区间中所有面积的和为1。在计算统计变量的概率密度值之后,利用威布尔分布函数对统计变量概率密度值进行拟合,得到威布尔分布函数的第一形状参数k1和第一比例参数λ1的拟合值,从而确定第一纵向动位移累积值序列中纵向动位移累积值的威布尔分布函数,如式(3)所示:
式中,g(C)表示第一纵向动位移累积值序列中纵向动位移累积值的威布尔分布函数,k1表示威布尔分布函数的第一形状参数,λ1表示威布尔分布函数的第一比例参数,e表示指数函数的底,e=2.71828。
将实测第一纵向动位移累积值序列C中的纵向动位移累积值及其概率密度值代入式(3)中,并进行最小二乘拟合,便可确定g(C)中参数k1和λ1的具体取值。最小二乘法是函数拟合中最常用的一种方法,是一种数学优化技术。它通过最小化误差的平方和寻找数据的最佳函数匹配。
将第二纵向动位移累积值序列D中的纵向动位移累积值视为统计变量,利用统计学方法计算统计变量概率密度值,在利用威布尔分布函数对统计变量概率密度值进行拟合,得到威布尔分布函数的第二形状参数k2和第二比例参数λ2的拟合值,从而确定第二纵向动位移累积值序列中纵向动位移累积值的威布尔分布函数,如式(4)所示:
式中,g(D)表示第二纵向动位移累积值序列D中纵向动位移累积值的威布尔分布函数,k2表示威布尔分布函数的第二形状参数,λ2表示威布尔分布的第二比例参数。
将实测第二纵向动位移累积值序列D中的纵向动位移累积值及其概率密度值代入式(4)中并进行最小二乘拟合,便可确定g(D)中参数k2和λ2的具体取值。
步骤60):利用蒙特卡洛抽样方法确定支座在设计使用寿命内的支座累积行程模拟值。
步骤60)具体包括:设支座的设计使用寿命为R年,利用第一纵向动位移累积值序列中纵向动位移累积值的威布尔分布函数g(C)对支座在R年内每天凌晨0时至6时的所有纵向动位移累积值进行蒙特卡洛抽样模拟,模拟总数NC=365Rn1。其中,蒙特卡罗抽样模拟方法又称统计模拟法、随机抽样技术,是使用随机数将所求解的问题同概率模型相联系,用电子计算机实现统计模拟或抽样,以获得问题的近似解。利用第二纵向动位移累积值序列D中纵向动位移累积值的威布尔分布函数g(D)对支座在R年内每天凌晨6时至24时的所有纵向动位移累积值进行蒙特卡洛抽样模拟,模拟总数ND=365Rn2;
在R年内共得到31536000Rf个模拟的纵向动位移累积值,即NC+ND=31536000Rf,将31536000Rf个纵向动位移累积模拟值进行求和,得到支座在设计使用寿命内的累积行程模拟值。
步骤70):将支座累积行程模拟值视为随机变量,并确定支座累积行程模拟值的广义极值分布函数。
步骤70)具体包括:将步骤50)和步骤60)重复操作100次,得到100个累积行程模拟值,将累积行程模拟值视为统计变量,并采用H表示;利用统计学方法计算所述100个累积行程模拟值的广义极值分布值,利用广义极值分布函数对累积行程模拟值广义极值分布值进行拟合,得到广义极值分布函数的形状参数r3、比例参数a3和位置参数b3的拟合值,从而确定支座累积行程模拟值广义极值分布函数,如式(5)所示:
式中,t(H)表示累积行程模拟值的广义极值分布函数,r3表示广义极值分布函数的形状参数,a3表示广义极值分布函数的比例参数,b3表示广义极值分布函数的位置参数。
将100个累积行程模拟值及其概率密度值代入式(5)中,并进行最小二乘拟合,便可确定t(H)中参数r3、a3和b3的具体取值。
步骤80):计算支座累积行程模拟值超过支座磨损上限值的失效概率,比较失效概率与预设概率阈值,如果失效概率大于或等于预设概率阈值,则支座达到磨损上限,更换支座。
步骤80)具体包括:设支座磨损上限值为[H],利用式(6)计算支座累积行程模拟值超过支座磨损上限值的失效概率β:
其中,d表示采样天数;
若β大于或等于预设概率阈值,则认为支座已经达到磨损上限,需要及时更换。作为优选,所述预设概率阈值为5%。
上述实施例的高速铁路桥梁支座在设计使用寿命内累积磨损的评估方法,支座磨损程度可以通过支座纵向动位移的累积行程来反映,因此本发明通过计算支座在服役期内的纵向动位移累积行程并与支座磨损上限值进行对比,来判断支座磨损程度。在计算支座纵向动位移累积行程时,现有方法不区分在0点至6点由随机环境荷载引起的纵向动位移和在6点至24点由随机环境荷载、列车荷载共同引起的纵向动位移。事实表明,这两种纵向动位移的概率密度特性是不一样的。本发明提出的方法充分考虑了这两种纵向动位移的差异性,并分别建立这两类纵向动位移的概率密度函数,使得评估结果更加准确。
步骤(1):在主梁上的支座位置安装纵向位移监测传感器,对支座纵向位移进行长期监测和数据采集;步骤(2):利用小波包分解法,提取出纵向位移监测数据中的动位移值;步骤(3):对基本时距内所有纵向动位移值的一阶差分绝对值累加得到纵向动位移累积值;步骤(4):将每天内的纵向动位移累积值划分为两类:由随机环境引起的纵向动位移累积值、随机环境与列车荷载共同引起的纵向动位移累积值;步骤(5):分别将两类纵向动位移累积值视为随机变量并确定其概率密度函数;步骤(6):利用蒙特卡洛抽样方法确定支座在设计使用寿命内的支座累积行程模拟值;步骤(7):将支座累积行程模拟值视为随机变量并确定其概率密度函数;步骤(8):计算支座累积行程模拟值超过支座磨损上限值的失效概率并进行安全评估。
下面例举一具体实施例:以南京大胜关长江大桥的支座纵向位移监测数据为例,说明本发明的具体实施过程以及评估效果。
步骤10)在南京大胜关长江大桥的支座位置安装纵向位移监测传感器,对支座纵向位移进行监测,并采集纵向位移监测数据,采样频率f为1Hz,采样天数d=1;支座纵向位移数据的采样总个数为L=86400。实际采集纵向位移监测数据,如图1所示。
步骤20)利用小波包分解技术,对支座纵向位移监测数据进行8层尺度的小波包分解,得到32个小波包分解系数,剔除第一个小波包分解系数后,对剩余小波包分解系数重构并叠加,得到支座纵向动位移值,如图2所示。
步骤30)将每天凌晨0时至6时的纵向动位移监测数据划分为200份,并对每份纵向动位移值作一阶差分处理,并取绝对值累加,得到第i天中第j份纵向动位移监测数据的动位移累积值Ci,j:
将每天凌晨6时至夜晚24时的纵向动位移监测数据划分为600份,并对每份纵向动位移值作一阶差分处理,并取绝对值累加,得到第m天中第n份纵向动位移监测数据的动位移累积值Dm,n:
步骤40)将凌晨0时至6时的纵向动位移累积值组合构成第一纵向动位移累积值序列C,如图3所示,C=[C1,1,C1,2,…,C1,200];将凌晨6时至24时的纵向动位移累积值组合构成第二纵向动位移累积值序列D,如图4所示,D=[D1,1,D1,2,…,D1,600]。
步骤50)将第一纵向动位移累积值序列C中的纵向动位移累积值视为统计变量,利用统计学方法计算统计变量的概率密度值,再利用威布尔分布函数g(C)对统计变量概率密度值进行拟合,如图5所示,g(C)的表达式如下:
将第一纵向动位移累积值序列C中的纵向动位移累积值及其概率密度值代入上式中,并进行最小二乘拟合,便可确定g(C)中参数k1和λ1的具体取值分别为:k1=3.03406,λ1=2.6409。得到威布尔分布函数的第一形状参数k1和第一比例参数λ1的拟合值,从而确定第一纵向动位移累积值序列中纵向动位移累积值的威布尔分布函数。
将第二纵向动位移累积值序列D中的纵向位移累积值视为统计变量,利用统计学方法计算其概率密度值,利用威布尔分布函数g(D)对其概率密度值进行拟合,如图6所示,其中g(D)的表达式如下:
将第二纵向动位移累积值序列D中的纵向动位移累积值及其概率密度值代入上式中,并进行最小二乘拟合,便可确定g(D)中参数k2和λ2的具体取值分别为:k2=2.50608,λ2=6.13737;
步骤60)设支座的设计使用寿命为15年。利用第一纵向动位移累积值序列中纵向动位移累积值的威布尔分布函数g(C)对支座在15年内每天凌晨0时至6时的所有纵向动位移累积值进行蒙特卡洛抽样模拟,模拟总个数NC=1095000;此外,利用第二纵向动位移累积值序列中纵向动位移累积值的威布尔分布函数g(D)对支座在15年内每天凌晨6时至24时的所有纵向动位移累积值进行蒙特卡洛抽样模拟,模拟总个数ND=3285000;在完成蒙特卡洛抽样模拟后,在15年内共得到4380000个模拟的纵向动位移累积值,即NC+ND=4380000,将这4380000个纵向动位移累积模拟值进行求和,得到支座在设计使用寿命内的累积行程模拟值;
步骤70)将步骤50)和步骤60)重复操作100次,得到100个累积行程模拟值,如图7所示。将累积行程模拟值视为统计变量并采用H表示,利用统计学方法计算这100个累积行程模拟值的广义极值分布值,在此基础上利用广义极值分布函数t(H)对广义极值分布值行拟合,如图8所示,其中t(H)的表达式为:
将100个累积行程模拟值及其广义极值分布值代入式中,并进行最小二乘拟合,得到r3=-0.200867,a3=0.004297,b3=20.47。
步骤80)支座磨损上限值在制作出厂时标定为[H]=25km,利用下式计算支座累积行程模拟值超过支座磨损上限值的失效概率β:
预设概率阈值为5%。因为β=0<5%,则认为支座尚未达到磨损上限,无需更换。
上述实施例的评估方法,可以更加真实、准确和全面地反映出高速铁路桥梁支座在设计使用寿命内的实际磨损程度,从而能够实现对磨损支座进行及时地保养、维护和更换。同时,上述实施例的评估方法完全可在高速铁路桥梁支座的累积磨损安全评估领域得到广泛推广和应用。
以上实施例仅仅是对本发明方案的进一步具体说明,在阅读了本发明实施例之后,本领域普通技术人员对本发明的各种等同形式的修改和替换均属于本发明申请权利要求所限定的保护的范围。
Claims (10)
1.一种高速铁路桥梁支座在设计使用寿命内累积磨损的评估方法,其特征在于,该方法包括如下步骤:
步骤10):在桥梁主梁的支座上安装纵向位移监测传感器,对支座纵向位移进行监测,并采集纵向位移监测数据,其中,第t时刻支座纵向位移值记为v(t),t表示时间,且t=1/f、2/f、…、L/f,f表示采样频率,L表示采样总数,L=86400df,d表示采样天数;
步骤20):利用小波包分解法提取纵向位移监测数据中的纵向动位移值;
步骤30):对所述纵向动位移值作一阶差分处理,并取绝对值累加,得到纵向动位移累积值;
步骤40):将每天的纵向动位移累积值划分为两类,构建第一纵向动位移累积值序列和第二纵向动位移累积值序列;
步骤50):将两类纵向动位移累积值视为随机变量,并确定纵向动位移累积值威布尔分布函数;
步骤60):利用蒙特卡洛抽样方法确定支座在设计使用寿命内的支座累积行程模拟值;
步骤70):将支座累积行程模拟值视为随机变量,并确定支座累积行程模拟值的广义极值分布函数;
步骤80):计算支座累积行程模拟值超过支座磨损上限值的失效概率,比较失效概率与预设概率阈值,如果失效概率大于或等于预设概率阈值,则支座达到磨损上限,更换支座。
2.按照权利要求1所述的高速铁路桥梁支座在设计使用寿命内累积磨损的评估方法,其特征在于,所述步骤20)包括:利用小波包分解技术,对支座纵向位移监测数据进行8层尺度的小波包分解,得到32个小波包分解系数,剔除第一个小波包分解系数后,对剩余小波包分解系数重构并叠加,得到支座纵向动位移值。
3.按照权利要求1所述的高速铁路桥梁支座在设计使用寿命内累积磨损的评估方法,其特征在于,所述步骤30)具体包括:将每天凌晨0时至6时的纵向动位移监测数据划分为n1份,其中第i天中第j份纵向动位移监测数据采用vi,j表示,其中i=1、2、…、d,j=1、2、…、n1,对每份纵向动位移值vi,j作一阶差分处理,并取绝对值累加,得到第i天中第j份纵向动位移监测数据的动位移累积值Ci,j,如式(1)所示:
式中,Ni,j表示第i天中第j份纵向动位移监测数据中的数据总数,vi,j,p+1表示第i天中第j份纵向动位移监测数据中的第p+1个值,vi,j,p表示第i天中第j份纵向动位移监测数据中的第p个值;
将每天凌晨6时至夜晚24时的纵向动位移监测数据划分为n2份,其中第m天中第n份纵向动位移监测数据采用vm,n表示,其中m=1、2、…、d,n=1、2、…、n2,对每份纵向动位移值vm,n作一阶差分处理,并取绝对值累加,得到第m天中第n份纵向动位移监测数据的动位移累积值Dm,n,如式(2)所示:
式中,Nm,n表示第m天中第n份纵向动位移监测数据中的数据总数,vm,n,q+1表示第m天中第n份纵向动位移监测数据中的第q+1个值,vm,n,q表示第m天中第n份纵向动位移监测数据中的第q个值。
4.按照权利要求1所述的高速铁路桥梁支座在设计使用寿命内累积磨损的评估方法,其特征在于,所述步骤40)具体包括:将所有天数在每天凌晨0时至6时的纵向动位移累积值组合构成第一纵向动位移累积值序列C,C=[C1,1、C1,2、…、C1,n1、C2,1、C2,2、…、C2,n1、…、Cd,1、Cd,2、…、Cd,n1];
将所有天数在每天凌晨6时至24时的纵向动位移累积值组合构成第二纵向动位移累积值序列D,D=[D1,1、D1,2、…、D1,n2、D2,1、D2,2、…、D2,n2、…、Dd,1、Dd,2、…、Dd,n2]。
5.按照权利要求1所述的高速铁路桥梁支座在设计使用寿命内累积磨损的评估方法,其特征在于,所述步骤50)具体包括:将第一纵向动位移累积值序列C中的纵向动位移累积值视为统计变量,利用统计学方法计算统计变量的概率密度值,再利用威布尔分布函数对统计变量概率密度值进行拟合,得到威布尔分布函数的第一形状参数k1和第一比例参数λ1的拟合值,从而确定第一纵向动位移累积值序列中纵向动位移累积值的威布尔分布函数,如式(3)所示:
式中,g(C)表示第一纵向动位移累积值序列中纵向动位移累积值的威布尔分布函数,k1表示威布尔分布函数的第一形状参数,λ1表示威布尔分布函数的第一比例参数,e表示指数函数的底,e=2.71828;
将第二纵向动位移累积值序列D中的纵向动位移累积值视为统计变量,利用统计学方法计算统计变量概率密度值,再利用威布尔分布函数对统计变量概率密度值进行拟合,得到威布尔分布函数的第二形状参数k2和第二比例参数λ2的拟合值,从而确定第二纵向动位移累积值序列中纵向动位移累积值的威布尔分布函数,如式(4)所示:
式中,g(D)表示第二纵向动位移累积值序列D中纵向动位移累积值的威布尔分布函数,k2表示威布尔分布函数的第二形状参数,λ2表示威布尔分布的第二比例参数。
6.按照权利要求1所述的高速铁路桥梁支座在设计使用寿命内累积磨损的评估方法,其特征在于,所述步骤60)具体包括:设支座的设计使用寿命为R年,利用第一纵向动位移累积值序列中纵向动位移累积值的威布尔分布函数g(C)对支座在R年内每天凌晨0时至6时的所有纵向动位移累积值进行蒙特卡洛抽样模拟,模拟总数NC=365Rn1;利用第二纵向动位移累积值序列中纵向动位移累积值的威布尔分布函数g(D)对支座在R年内每天凌晨6时至24时的所有纵向动位移累积值进行蒙特卡洛抽样模拟,模拟总数ND=365Rn2;
在R年内共得到31536000Rf个模拟的纵向动位移累积值,即NC+ND=31536000Rf,将31536000Rf个纵向动位移累积模拟值进行求和,得到支座在设计使用寿命内的累积行程模拟值。
7.按照权利要求1所述的高速铁路桥梁支座在设计使用寿命内累积磨损的评估方法,其特征在于,所述步骤70)具体包括:将步骤50)和步骤60)重复操作100次,得到100个累积行程模拟值,将累积行程模拟值视为统计变量,并采用H表示;利用统计学方法计算所述100个累积行程模拟值的广义极值分布值,利用广义极值分布函数对累积行程模拟值广义极值分布值进行拟合,得到广义极值分布函数的形状参数r3、比例参数a3和位置参数b3的拟合值,从而确定支座累积行程模拟值广义极值分布函数,如式(5)所示:
式中,t(H)表示累积行程模拟值的广义极值分布函数,r3表示广义极值分布函数的形状参数,a3表示广义极值分布函数的比例参数,b3表示广义极值分布函数的位置参数。
8.按照权利要求1所述的高速铁路桥梁支座在设计使用寿命内累积磨损的评估方法,其特征在于,所述步骤80)具体包括:设支座磨损上限值为[H],利用式(6)计算支座累积行程模拟值超过支座磨损上限值的失效概率β:
其中,d表示采样天数;
若β大于或等于预设概率阈值,则认为支座已经达到磨损上限,需要及时更换。
9.按照权利要求8所述的高速铁路桥梁支座在设计使用寿命内累积磨损的评估方法,其特征在于,所述预设概率阈值为5%。
10.按照权利要求1所述的高速铁路桥梁支座在设计使用寿命内累积磨损的评估方法,其特征在于,所述步骤10)中,采样频率f≥1Hz,采样天数d≥1;所述步骤20)中,n1的取值在30~360之间,n2的取值在90~1080之间。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710932238.1A CN107766630B (zh) | 2017-10-10 | 2017-10-10 | 高速铁路桥梁支座在设计使用寿命内累积磨损的评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710932238.1A CN107766630B (zh) | 2017-10-10 | 2017-10-10 | 高速铁路桥梁支座在设计使用寿命内累积磨损的评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107766630A true CN107766630A (zh) | 2018-03-06 |
CN107766630B CN107766630B (zh) | 2019-01-25 |
Family
ID=61267016
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710932238.1A Active CN107766630B (zh) | 2017-10-10 | 2017-10-10 | 高速铁路桥梁支座在设计使用寿命内累积磨损的评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107766630B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108897917A (zh) * | 2018-05-31 | 2018-11-27 | 南京东南建筑机电抗震研究院有限公司 | 一种用于高速铁路桥梁支座易损性评估的方法 |
CN114756946A (zh) * | 2022-05-13 | 2022-07-15 | 东南大学 | 一种基于短期检测数据的桥梁梁端长期累积位移估计方法 |
CN116956664A (zh) * | 2023-06-26 | 2023-10-27 | 中交公路长大桥建设国家工程研究中心有限公司 | 基于监测数据的桥梁支座滑板磨损及寿命预测方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005214828A (ja) * | 2004-01-30 | 2005-08-11 | Nsk Ltd | 故障判定手段を有する変位測定装置および変位測定装置の故障検出方法 |
CN201567528U (zh) * | 2009-11-27 | 2010-09-01 | 成都市虹筑路桥机械有限公司 | 桥梁型钢伸缩装置转动结构 |
CN103440404A (zh) * | 2013-08-07 | 2013-12-11 | 东南大学 | 基于横向风载效应的桥梁加劲梁抗侧力性能退化报警方法 |
CN103900785A (zh) * | 2014-04-14 | 2014-07-02 | 东南大学 | 一种确定大跨桥梁结构主梁横向动位移的方法 |
CN104392148A (zh) * | 2014-12-15 | 2015-03-04 | 重庆交通大学 | 大跨度轨道专用斜拉桥预拱度设置方法 |
-
2017
- 2017-10-10 CN CN201710932238.1A patent/CN107766630B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005214828A (ja) * | 2004-01-30 | 2005-08-11 | Nsk Ltd | 故障判定手段を有する変位測定装置および変位測定装置の故障検出方法 |
CN201567528U (zh) * | 2009-11-27 | 2010-09-01 | 成都市虹筑路桥机械有限公司 | 桥梁型钢伸缩装置转动结构 |
CN103440404A (zh) * | 2013-08-07 | 2013-12-11 | 东南大学 | 基于横向风载效应的桥梁加劲梁抗侧力性能退化报警方法 |
CN103900785A (zh) * | 2014-04-14 | 2014-07-02 | 东南大学 | 一种确定大跨桥梁结构主梁横向动位移的方法 |
CN104392148A (zh) * | 2014-12-15 | 2015-03-04 | 重庆交通大学 | 大跨度轨道专用斜拉桥预拱度设置方法 |
Non-Patent Citations (4)
Title |
---|
QIHUAN HUANG等: "MONITORING AND EVALUATION OF A LONG-SPAN RAIWAY BRIDGE USING SENTINEL-1 DATA", 《ISPRS ANNALS OF THE PHOTOGRAMMETRY, REMOTE SENSING AND SPATIAL INFORMATION SCIENCES》 * |
WANG G X等: "Detection and location of the degraded bearings based on monitoring the longitudinal expansion performance of the main girder of the Dashengguan Yangtze Bridge ", 《JOURNAL OF PERFORMANCE OF CONSTRUCTED FACILITIES》 * |
周凯: "基于静态监测数据的桥梁结构状态评估和预警方法研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
黄其欢 等: "基于InSAR的南京大胜关大桥纵向位移监测与分析", 《东南大学学报》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108897917A (zh) * | 2018-05-31 | 2018-11-27 | 南京东南建筑机电抗震研究院有限公司 | 一种用于高速铁路桥梁支座易损性评估的方法 |
CN114756946A (zh) * | 2022-05-13 | 2022-07-15 | 东南大学 | 一种基于短期检测数据的桥梁梁端长期累积位移估计方法 |
CN116956664A (zh) * | 2023-06-26 | 2023-10-27 | 中交公路长大桥建设国家工程研究中心有限公司 | 基于监测数据的桥梁支座滑板磨损及寿命预测方法及系统 |
CN116956664B (zh) * | 2023-06-26 | 2025-04-15 | 中交公路长大桥建设国家工程研究中心有限公司 | 基于监测数据的桥梁支座滑板磨损及寿命预测方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN107766630B (zh) | 2019-01-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109472110B (zh) | 一种基于lstm网络和arima模型的航空发动机剩余使用寿命预测方法 | |
CN111191191B (zh) | 一种精准预测混凝土坝变形效应的组合模型的构建方法 | |
Lin et al. | An approach to signal processing and condition-based maintenance for gearboxes subject to tooth failure | |
CN112083244B (zh) | 综合化航空电子设备故障智能诊断系统 | |
CN103278343B (zh) | 一种结构关键疲劳构件的识别方法 | |
Cheung et al. | The application of statistical pattern recognition methods for damage detection to fielddata | |
Ding et al. | Assessment of bridge expansion joints using long-term displacement measurement under changing environmental conditions | |
CN107871161A (zh) | 一种基于神经网络的桥梁整体损伤预警方法 | |
CN103293014A (zh) | 一种桥梁疲劳损伤状态与剩余寿命的评估方法 | |
CN110175369A (zh) | 一种基于二维卷积神经网络的齿轮剩余寿命预测方法 | |
CN107766630A (zh) | 高速铁路桥梁支座在设计使用寿命内累积磨损的评估方法 | |
CN115408755B (zh) | 一种考虑时变效应的组合梁桥动力疲劳可靠性评估方法 | |
CN108170925A (zh) | 一种基于arma模型的桥梁损伤预警方法 | |
CN101799674A (zh) | 一种数控装备服役状态分析方法 | |
CN107480731A (zh) | 一种火电厂汽动给水泵组故障特征的早期识别方法 | |
Li et al. | Data-driven structural condition assessment for high-speed railway bridges using multi-band FIR filtering and clustering | |
CN102937534A (zh) | 一种基于组合预测模型对机电设备的故障预测方法 | |
CN113673010A (zh) | 一种基于监测数据的钢箱梁评估方法与系统 | |
CN108304348A (zh) | 一种基于二元维纳过程的轴承剩余寿命预测方法 | |
CN113204850A (zh) | 一种桥梁挠度监测中的温度效应分离方法 | |
CN114757365A (zh) | 一种基于深度学习的高速铁路路基沉降预测及预警方法 | |
Lee et al. | Response pattern analysis‐based structural health monitoring of cable‐stayed bridges | |
CN118071136A (zh) | 一种水库大坝风险评估方法 | |
CN117972842A (zh) | 一种基于Prophet+LSTM模型的桥梁应变应力数据预测方法 | |
CN113673145A (zh) | 基于pca-lstm神经网络的地表沉降预测方法 |
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 | ||
CB03 | Change of inventor or designer information | ||
CB03 | Change of inventor or designer information |
Inventor after: Wang Gaoxin Inventor after: Ye Jihong Inventor after: Dong Zhaoxing Inventor after: Liu Zhijun Inventor before: Wang Gaoxin Inventor before: Ye Jihong Inventor before: Dong Zhaoxing |
|
GR01 | Patent grant | ||
GR01 | Patent grant |