CN103837358B - The method for early warning of the overall lateral resistance behavior exception of long-span bridges - Google Patents
The method for early warning of the overall lateral resistance behavior exception of long-span bridges Download PDFInfo
- Publication number
- CN103837358B CN103837358B CN201410054835.5A CN201410054835A CN103837358B CN 103837358 B CN103837358 B CN 103837358B CN 201410054835 A CN201410054835 A CN 201410054835A CN 103837358 B CN103837358 B CN 103837358B
- Authority
- CN
- China
- Prior art keywords
- static
- sequence
- wind speed
- displacement
- month
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 22
- 230000003068 static effect Effects 0.000 claims abstract description 133
- 238000006073 displacement reaction Methods 0.000 claims abstract description 81
- 230000000694 effects Effects 0.000 claims abstract description 14
- 239000013598 vector Substances 0.000 claims abstract description 14
- 230000005856 abnormality Effects 0.000 claims abstract description 10
- 238000004458 analytical method Methods 0.000 claims abstract description 10
- 238000012544 monitoring process Methods 0.000 claims abstract description 10
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 9
- 230000002159 abnormal effect Effects 0.000 claims abstract description 7
- 238000012545 processing Methods 0.000 claims abstract description 5
- 238000012360 testing method Methods 0.000 claims description 11
- 230000008859 change Effects 0.000 claims description 9
- 238000001744 unit root test Methods 0.000 claims description 9
- 238000007689 inspection Methods 0.000 claims description 7
- 238000010219 correlation analysis Methods 0.000 claims description 3
- 238000012423 maintenance Methods 0.000 claims description 3
- 238000011160 research Methods 0.000 abstract description 3
- 230000015556 catabolic process Effects 0.000 description 4
- 238000006731 degradation reaction Methods 0.000 description 4
- 238000013461 design Methods 0.000 description 3
- 230000007774 longterm Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 229910000831 Steel Inorganic materials 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000000780 augmented Dickey–Fuller test Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
大跨桥梁结构整体抗侧力性能异常的预警方法包含如下步骤:(1)对主梁跨中部位的三维风场和GPS位移效应进行数据采集;(2)对采集数据进行矢量分解和均值处理,得到横桥向静风速序列和静位移序列;(3)对静风速序列和静位移序列按月份划分,并利用小波包互相关系数确定两者在同月份的主相关序列;(4)依次对各月主相关序列进行傅里叶级数拟合,利用拟合参数值的月变化特征对大跨桥梁结构的整体抗侧力性能进行异常预警。本发明提出的一种基于实测风载效应的大跨桥梁结构整体抗侧力性能异常预警方法,弥补了大跨桥梁结构在整体抗侧力性能异常预警方面的研究空白,可为大跨桥梁结构整体抗侧力性能的监测和分析工作提供重要参考。
The early warning method for the abnormality of the overall lateral force resistance of long-span bridge structures includes the following steps: (1) Collect data on the three-dimensional wind field and GPS displacement effect at the mid-span of the main beam; (2) Perform vector decomposition and mean value processing on the collected data , to obtain the static wind speed sequence and static displacement sequence of the cross bridge; (3) divide the static wind speed sequence and static displacement sequence by month, and use the wavelet packet cross-correlation coefficient to determine the main correlation sequence of the two in the same month; (4) sequentially Fourier series fitting is carried out on the main correlation sequence of each month, and the abnormal early warning of the overall lateral force resistance performance of the long-span bridge structure is carried out by using the monthly variation characteristics of the fitting parameter values. The present invention proposes an early warning method for the abnormality of the overall lateral force performance of long-span bridge structures based on the measured wind load effect, which makes up for the research gap in the early warning of the overall lateral force performance of long-span bridge structures, and can be used for long-span bridge structures. It provides an important reference for the monitoring and analysis of the overall lateral force resistance performance.
Description
技术领域 technical field
本发明涉及一种基于实测风载效应的大跨桥梁结构整体抗侧力性能异常预警方法。 The invention relates to an early warning method for abnormality of the overall lateral force resistance performance of a long-span bridge structure based on measured wind load effects.
背景技术 Background technique
风载作用下大跨桥梁结构的整体抗侧力性能直接关系到整个桥梁结构在运营期间的正常使用和安全性能。在100年设计服役期内,桥梁构件由于长期受到气候、环境、荷载等外界因素影响,其结构材料会被逐渐腐蚀、老化和形成损伤累积,使得桥梁结构的整体抗侧力性能逐渐发生退化。然而,目前针对桥梁结构整体抗侧力性能退化的实时监测和分析工作甚少,普遍将桥梁结构的整体抗侧力性能始终视为运营初期的完好状态,而且桥梁结构重要构件的侧向抗力设计也很少会考虑到“抗侧力性能退化”这一影响因素。可见,工程界对于大跨桥梁结构在服役期间的整体抗侧力性能退化行为缺乏足够认识,有必要深入研究大跨桥梁结构的整体抗侧力性能退化分析方法。由于大跨桥梁结构在横桥向主要承受风荷载作用,由风荷载引起的主梁侧向位移响应大小反映着整个桥梁结构的抗侧力性能,这一点为开展大跨桥梁结构的整体抗侧力性能分析工作提供了契机。 The overall lateral force resistance of long-span bridge structures under wind loads is directly related to the normal use and safety performance of the entire bridge structure during operation. During the 100-year design service period, due to the long-term influence of external factors such as climate, environment, and load, the structural materials of bridge components will be gradually corroded, aged, and accumulated damage, which will gradually degrade the overall lateral force resistance of the bridge structure. However, at present, there is little real-time monitoring and analysis of the degradation of the overall lateral force resistance of bridge structures. The overall lateral force resistance of bridge structures is generally considered to be in good condition at the initial stage of operation. The influencing factor of "degradation of lateral force resistance" is rarely taken into account. It can be seen that the engineering community lacks sufficient understanding of the overall lateral force performance degradation behavior of long-span bridge structures during service, and it is necessary to study the overall lateral force performance degradation analysis methods of long-span bridge structures. Since the long-span bridge structure mainly bears the wind load in the transverse bridge direction, the lateral displacement response of the main girder caused by the wind load reflects the lateral force resistance performance of the entire bridge structure. Force performance analysis work provides an opportunity.
鉴于此,本发明提出一种基于实测风载效应的大跨桥梁结构整体抗侧力性能异常预警方法。 In view of this, the present invention proposes an early warning method for abnormality of the overall lateral force resistance performance of a long-span bridge structure based on the measured wind load effect.
发明内容 Contents of the invention
技术问题:本发明针对现有技术中关于大跨桥梁结构在整体抗侧力性能评价方面的研究空白,提出一种基于实测风载效应的大跨桥梁结构整体抗侧力性能异常预警方法。 Technical problem: The present invention aims at the research gap in the prior art on the evaluation of the overall lateral force resistance performance of long-span bridge structures, and proposes an early warning method for the abnormality of the overall lateral force resistance performance of long-span bridge structures based on the measured wind load effect.
技术方案:为解决上述技术问题,本发明的一种基于实测风载效应的大跨桥梁结构整体抗侧力性能异常预警方法采用如下技术方案: Technical solution: In order to solve the above-mentioned technical problems, the present invention adopts the following technical solution for the abnormal early warning method of the overall lateral force resistance performance of long-span bridge structures based on the measured wind load effect:
本发明提出的一种基于实测风载效应的大跨桥梁结构整体抗侧力性能异常预警方法,该方法具体包括如下步骤: The present invention proposes a method for early warning of abnormality of the overall lateral force resistance performance of long-span bridge structures based on measured wind load effects. The method specifically includes the following steps:
步骤(1):对主梁跨中部位的三维风场和GPS位移效应进行数据采集: Step (1): Collect data on the 3D wind field and GPS displacement effect at the mid-span of the main girder:
在大跨桥梁的主梁跨中处安装三维超声风速仪和GPS位移监测站,对此处风向量v(t)及位移向量u(t)进行实时监测并以时间序列存储,其中v(t)=[vr(t),α(t),β(t)],u(t)=[ux(t),uy(t),uz(t)],vr(t),α(t),β(t)分别为绝对风速、风攻角与风向角,ux(t),uy(t),uz(t)分别为GPS坐标系下的三轴方向位移,t表示时间,t=1,2,...,L,单位为秒,L表示时间长度; A three-dimensional ultrasonic anemometer and a GPS displacement monitoring station are installed at the middle of the main girder span of the long-span bridge, where the wind vector v(t) and displacement vector u(t) are monitored in real time and stored in time series, where v(t )=[v r (t), α(t), β(t)], u(t)=[u x (t), u y (t), u z (t)], v r (t) , α(t), β(t) are the absolute wind speed, wind attack angle and wind direction angle respectively, u x (t), u y (t), u z (t) are the three-axis direction displacement in the GPS coordinate system , t represents time, t=1,2,...,L, the unit is second, L represents the length of time;
步骤(2):对采集数据进行矢量分解和均值处理,得到横桥向静风速序列和静位移序列: Step (2): Perform vector decomposition and mean value processing on the collected data to obtain the static wind speed sequence and static displacement sequence of the cross bridge:
利用以下两式将时间序列v(t)、u(t)进行矢量分解,得到横桥向风速时程vh(t)和位移时程ur(t): Use the following two formulas to decompose the time series v(t) and u(t) into vectors to obtain the time history of wind speed v h (t) and the time history of displacement u r (t) across the bridge:
vh(t)=vr(t)·cos(α(t))·sin(β(t)) v h (t) = v r (t)·cos(α(t))·sin(β(t))
ur(t)=ux(t)·sin(γ)-uy(t)·cos(γ) u r (t)=u x (t)·sin(γ)-u y (t)·cos(γ)
式中γ表示GPS坐标系中的x轴与主梁纵向轴线的夹角;后将L划分为n个10min时间段,并利用下式计算每个时间段内对应vh(t)和ur(t)的平均值,得到静风速序列vm(k)和静位移序列um(k): In the formula, γ represents the angle between the x-axis in the GPS coordinate system and the longitudinal axis of the main girder; then divide L into n 10-min time periods, and use the following formula to calculate the corresponding v h (t) and u r in each time period (t) to get the static wind speed sequence v m (k) and the static displacement sequence u m (k):
式中k=1,2,…,n; In the formula, k=1,2,...,n;
步骤(3):对静风速序列和静位移序列按月份划分,并利用小波包互相关系数确定两者在同月份的主相关序列: Step (3): Divide the static wind speed sequence and static displacement sequence by month, and use the wavelet packet cross-correlation coefficient to determine the main correlation sequence of the two in the same month:
①对静风速序列和静位移序列按月份进行划分,以第q个月份的划分结果为分析例,将此月静风速值按递增排列并等间距划分为p段,其中静风速值排列前后的位置变化记为R(n1,n2),R(n1,n2)具体表示第n1个静风速值在排列后位于第n2个位置;此外按照R(n1,n2)的排列规则将此月静位移值排列并同样均分为p段; ① Divide the static wind speed sequence and the static displacement sequence by month. Taking the division result of the qth month as an analysis example, the monthly static wind speed values are arranged in increasing order and divided into p segments at equal intervals. Among them, the static wind speed values before and after the arrangement are The position change is recorded as R(n 1 ,n 2 ), and R(n 1 ,n 2 ) specifically means that the n 1st static wind speed value is located at the n 2nd position after arrangement; in addition, according to R(n 1 ,n 2 ) Arrange the monthly static displacement values according to the arrangement rules and equally divide them into p segments;
②以第s段静风速值和静位移值为例,对此段静风速值和静位移值进行第4尺度小波包分解,两者均得到按结点位置排列的16个小波包系数,利用下式对相同结点位置小波包系数的重构序列进行逐一互相关分析: ② Taking the static wind speed value and static displacement value of the s segment as an example, the fourth-scale wavelet packet decomposition is performed on the static wind speed value and static displacement value of this segment, and 16 wavelet packet coefficients arranged according to the node positions are obtained for both of them. The following formula performs cross-correlation analysis one by one on the reconstructed sequence of wavelet packet coefficients at the same node position:
式中,rvu(g)表示第g个小波包系数的重构静风速序列和静位移序列之间的互相关系数,分别为第g个小波包系数的第ng个重构静风速值和静位移值,nt(g)为第g个小波包系数的重构静风速值总个数,分别为第g个小波包系数的重构静风速序列和静位移序列的的均值,g=1,2,...,16; In the formula, r vu (g) represents the cross-correlation coefficient between the reconstructed static wind speed sequence and static displacement sequence of the gth wavelet packet coefficient, are respectively the n gth reconstructed static wind velocity value and static displacement value of the gth wavelet packet coefficient, n t (g) is the total number of reconstructed static wind velocity values of the gth wavelet packet coefficient, are the mean values of the reconstructed static wind speed sequence and static displacement sequence of the gth wavelet packet coefficient respectively, g=1,2,...,16;
③从16个小波包系数中剔除掉互相关系数绝对值小于0.9对应的小波包系数,后对剩余小波包系数的重构序列叠加得到第s段静风速值和静位移值的重构序列,按此方法将得到的p段静风速值和静位移值的重构序列重新组合成第q个月份的静风速和静位移序列,以此作为静风速和静位移在此月的主相关序列,后对各月遍历得到所有月份静风速和静位移之间的主相关序列; ③ Eliminate the wavelet packet coefficients corresponding to the absolute value of the cross-correlation coefficient less than 0.9 from the 16 wavelet packet coefficients, and then superimpose the reconstruction sequence of the remaining wavelet packet coefficients to obtain the reconstruction sequence of the static wind speed value and static displacement value of the s segment, According to this method, the reconstructed sequence of the static wind speed and static displacement of the p segment is recombined into the static wind speed and static displacement sequence of the qth month, which is used as the main correlation sequence of the static wind speed and static displacement in this month, Afterwards, the main correlation sequence between the static wind speed and static displacement of all months is obtained by traversing each month;
步骤(4):依次对各月主相关序列进行傅里叶级数拟合,利用拟合参数值的月变化特征对大跨桥梁结构的整体抗侧力性能进行异常预警: Step (4): Carry out Fourier series fitting on the main correlation sequence of each month in turn, and use the monthly variation characteristics of the fitted parameter values to give an abnormal early warning of the overall lateral force resistance performance of the long-span bridge structure:
利用下式所示的二阶傅里叶级数,依次对各月主相关序列进行最小二乘拟合并确定各月参数估计值: Using the second-order Fourier series shown in the following formula, the least squares fitting is performed on the main correlation series of each month in turn and the parameter estimates of each month are determined:
式中,vm表示第m个月主相关序列中的静风速值,um表示第m个月主相关序列中的静位移值,ae(m)、bx(m)和w(m)分别为第m个月傅里叶级数的参数估计值。各参数估计值的月变化平稳特性直接关系到大跨桥梁结构的整体抗侧力性能,因此分别对其进行ADF单位根检验。(ADF单位根检验是判断时间序列平稳性的一种统计分析方法,又称为增广迪基-福勒检验,具体可通过调用MATLAB数学软件中的函数命令adftest得到检验结果,若检验结果不拒绝存在一个单位根的原假设,则表明时间序列具有非平稳性;若检验结果拒绝存在一个单位根的原假设,则表明时间序列具有平稳性。)对于各参数估计值的月变化序列的ADF检验结果可分为以下四种情况: In the formula, v m represents the static wind speed value in the mth month main correlation sequence, u m represents the static displacement value in the mth month main correlation sequence, a e (m), b x (m) and w (m ) are the parameter estimates of the Fourier series for the mth month, respectively. The monthly variation of the estimated value of each parameter is directly related to the overall lateral force resistance performance of the long-span bridge structure, so the ADF unit root test is carried out separately. (ADF unit root test is a statistical analysis method for judging the stationarity of time series, also known as Augmented Dickey-Fuller test, the test result can be obtained by calling the function command adftest in the MATLAB mathematical software, if the test result is not Rejecting the null hypothesis that there is a unit root indicates that the time series is non-stationary; if the test result rejects the null hypothesis that there is a unit root, it indicates that the time series is stationary.) For the ADF of the monthly change series of each parameter estimate The test results can be divided into the following four situations:
①若每个参数估计值的检验结果均拒绝存在一个单位根的原假设,则大跨桥梁结构的整体抗侧力性能处于良好状态;②若存在1至2个参数估计值的检验结果不拒绝存在一个单位根的原假设,则对大跨桥梁结构的整体抗侧力性能进行黄色预警,并密切跟踪各参数估计值月变化特征的ADF单位根检验结果;③若存在3至4个参数估计值的检验结果不拒绝存在一个单位根的原假设,则对大跨桥梁结构的整体抗侧力性能进行橙色预警,并派桥梁检修人员对桥梁结构关键构件和部位进行现场检测,根据检测结果采取应对措施;④若存在4个以上参数估计值的检验结果不拒绝存在一个单位根的原假设,则对大跨桥梁结构的整体抗侧力性能进行红色预警,并派桥梁专业人员对桥梁结构的整体抗侧力性能进行安全评估和决策。 ①If the test results of each parameter estimate reject the null hypothesis that there is a unit root, the overall lateral force resistance of the long-span bridge structure is in a good state; ②If there are 1 to 2 parameter estimates, the test results do not reject If there is a null hypothesis of a unit root, a yellow warning will be given to the overall lateral force resistance performance of the long-span bridge structure, and the ADF unit root test results of the monthly variation characteristics of the estimated values of each parameter will be closely tracked; ③ If there are 3 to 4 parameter estimates If the test result of the value does not reject the null hypothesis that there is a unit root, an orange warning will be given to the overall lateral force resistance of the long-span bridge structure, and bridge maintenance personnel will be sent to conduct on-site inspections on the key components and parts of the bridge structure. Countermeasures; ④ If the test results of more than 4 parameter estimates do not reject the null hypothesis that there is a unit root, a red warning will be given to the overall lateral force resistance of the long-span bridge structure, and bridge professionals will be sent to inspect the bridge structure. The overall lateral force resistance performance is used for safety assessment and decision-making.
有益效果:本发明提出的一种基于实测风载效应的大跨桥梁结构整体抗侧力性能异常预警方法,利用小波包互相关系数分析法提取出横桥向静风速与静位移之间的主相关序列,并利用傅里叶级数参数拟合值的月变化特性对大跨桥梁结构的整体抗侧力性能进行异常预警,弥补了大跨桥梁结构在整体抗侧力性能异常预警方面的研究空白,可为大跨桥梁结构整体抗侧力性能的监测和分析工作提供重要参考。 Beneficial effects: The present invention proposes a method for early warning of anomalies in overall lateral force resistance of long-span bridge structures based on measured wind load effects. Correlation sequence, and use the monthly variation characteristics of the Fourier series parameter fitting value to carry out abnormal early warning of the overall lateral force resistance performance of long-span bridge structures, making up for the research on the abnormal early warning of overall lateral force resistance performance of long-span bridge structures It can provide an important reference for the monitoring and analysis of the overall lateral force performance of long-span bridge structures.
附图说明 Description of drawings
图1为本发明实施例苏通大桥主梁侧向风载效应监测点布置(单位:m); Fig. 1 is the arrangement of monitoring points for lateral wind load effect of main girder of Sutong Bridge according to the embodiment of the present invention (unit: m);
图2为本发明实施例横桥向静风速序列; Fig. 2 is the static wind speed sequence of the transverse bridge of the embodiment of the present invention;
图3为本发明实施例横桥向静位移序列; Fig. 3 is the static displacement sequence of the horizontal bridge according to the embodiment of the present invention;
图4为本发明实施例月静风速值与月静位移值之间的11段相关性散点图; Fig. 4 is 11 segment correlation scatter plots between monthly static wind speed value and monthly static displacement value in the embodiment of the present invention;
图5为本发明实施例月静风速与月静位移的11段主相关序列; Fig. 5 is 11 sections main correlation sequences of monthly static wind speed and monthly static displacement of the embodiment of the present invention;
图6是傅里叶级数拟合参数值的月变化特征中参数a0、a1和a2变化时程; Fig. 6 is the change time course of parameters a 0 , a 1 and a 2 in the monthly variation characteristics of Fourier series fitting parameter values;
图7是傅里叶级数拟合参数值的月变化特征中参数b1和b2变化时程; Fig. 7 is parameter b 1 and b 2 change time course in the month change characteristic of Fourier series fitting parameter value;
图8是傅里叶级数拟合参数值的月变化特征中参数w变化时程。 Fig. 8 is the change time course of parameter w in the monthly change feature of Fourier series fitting parameter value.
具体实施方式 detailed description
下面结合附图和具体实施例,进一步阐明本发明,应理解这些实施例仅用于说明本发明而不用于限制本发明的范围,在阅读了本发明之后,本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。 Below in conjunction with accompanying drawing and specific embodiment, further illustrate the present invention, should be understood that these embodiments are only for illustrating the present invention and are not intended to limit the scope of the present invention, after having read the present invention, those skilled in the art will understand various aspects of the present invention Modifications in equivalent forms all fall within the scope defined by the appended claims of this application.
本发明的一种基于实测风载效应的大跨桥梁结构整体抗侧力性能异常预警方法,该方法具体包括如下步骤: A method for early warning of abnormality of the overall lateral force resistance performance of a long-span bridge structure based on the measured wind load effect of the present invention, the method specifically includes the following steps:
步骤(1):对主梁跨中部位的三维风场和GPS位移效应进行数据采集: Step (1): Collect data on the 3D wind field and GPS displacement effect at the mid-span of the main girder:
在大跨桥梁的主梁跨中处安装三维超声风速仪和GPS位移监测站,对此处风向量v(t)及位移向量u(t)进行实时监测并以时间序列存储,其中v(t)=[vr(t),α(t),β(t)],u(t)=[ux(t),uy(t),uz(t)],vr(t),α(t),β(t)分别为绝对风速、风攻角与风向角,ux(t),uy(t),uz(t)分别为GPS坐标系下的三轴方向位移,t表示时间,t=1,2,...,L,单位为秒,L表示时间长度; A three-dimensional ultrasonic anemometer and a GPS displacement monitoring station are installed at the middle of the main girder span of the long-span bridge, where the wind vector v(t) and displacement vector u(t) are monitored in real time and stored in time series, where v(t )=[v r (t), α(t), β(t)], u(t)=[u x (t), u y (t), u z (t)], v r (t) , α(t), β(t) are the absolute wind speed, wind attack angle and wind direction angle respectively, u x (t), u y (t), u z (t) are the three-axis direction displacement in the GPS coordinate system , t represents time, t=1,2,...,L, the unit is second, L represents the length of time;
步骤(2):对采集数据进行矢量分解和均值处理,得到横桥向静风速序列和静位移序列: Step (2): Perform vector decomposition and mean value processing on the collected data to obtain the static wind speed sequence and static displacement sequence of the cross bridge:
利用以下两式将时间序列v(t)、u(t)进行矢量分解,得到横桥向风速时程vh(t)和位移时程ur(t): Use the following two formulas to decompose the time series v(t) and u(t) into vectors to obtain the time history of wind speed v h (t) and the time history of displacement u r (t) across the bridge:
vh(t)=vr(t)·cos(α(t))·sin(β(t)) v h (t) = v r (t)·cos(α(t))·sin(β(t))
ur(t)=ux(t)·sin(γ)-uy(t)·cos(γ) u r (t)=u x (t)·sin(γ)-u y (t)·cos(γ)
式中γ表示GPS坐标系中的x轴与主梁纵向轴线的夹角;后将L划分为n个10min时间段,并利用下式计算每个时间段内对应vh(t)和ur(t)的平均值,得到静风速序列vm(k)和静位移序列um(k): In the formula, γ represents the angle between the x-axis in the GPS coordinate system and the longitudinal axis of the main girder; then divide L into n 10-min time periods, and use the following formula to calculate the corresponding v h (t) and u r in each time period (t) to get the static wind speed sequence v m (k) and the static displacement sequence u m (k):
式中k=1,2,…,n; In the formula, k=1,2,...,n;
步骤(3):对静风速序列和静位移序列按月份划分,并利用小波包互相关系数确定两者在同月份的主相关序列: Step (3): Divide the static wind speed sequence and static displacement sequence by month, and use the wavelet packet cross-correlation coefficient to determine the main correlation sequence of the two in the same month:
①对静风速序列和静位移序列按月份进行划分,以第q个月份的划分结果为分析例,将此月静风速值按递增排列并等间距划分为p段,其中静风速值排列前后的位置变化记为R(n1,n2),R(n1,n2)具体表示第n1个静风速值在排列后位于第n2个位置;此外按照R(n1,n2)的排列规则将此月静位移值排列并同样均分为p段; ① Divide the static wind speed sequence and the static displacement sequence by month. Taking the division result of the qth month as an analysis example, the monthly static wind speed values are arranged in increasing order and divided into p segments at equal intervals. Among them, the static wind speed values before and after the arrangement are The position change is recorded as R(n 1 ,n 2 ), and R(n 1 ,n 2 ) specifically means that the n 1st static wind speed value is located at the n 2nd position after arrangement; in addition, according to R(n 1 ,n 2 ) Arrange the monthly static displacement values according to the arrangement rules and equally divide them into p segments;
②以第s段静风速值和静位移值为例,对此段静风速值和静位移值进行第4尺度小波包分解,两者均得到按结点位置排列的16个小波包系数,利用下式对相同结点位置小波包系数的重构序列进行逐一互相关分析: ② Taking the static wind speed value and static displacement value of the s segment as an example, the fourth-scale wavelet packet decomposition is performed on the static wind speed value and static displacement value of this segment, and 16 wavelet packet coefficients arranged according to the node positions are obtained for both of them. The following formula performs cross-correlation analysis one by one for the reconstructed sequence of wavelet packet coefficients at the same node position:
式中,rvu(g)表示第g个小波包系数的重构静风速序列和静位移序列之间的互相关系数,分别为第g个小波包系数的第ng个重构静风速值和静位移值,nt(g)为第g个小波包系数的重构静风速值总个数,分别为第g个小波包系数的重构静风速序列和静位移序列的的均值,g=1,2,...,16; In the formula, r vu (g) represents the cross-correlation coefficient between the reconstructed static wind speed sequence and static displacement sequence of the gth wavelet packet coefficient, are respectively the n gth reconstructed static wind velocity value and static displacement value of the gth wavelet packet coefficient, n t (g) is the total number of reconstructed static wind velocity values of the gth wavelet packet coefficient, are the mean values of the reconstructed static wind speed sequence and static displacement sequence of the gth wavelet packet coefficient respectively, g=1,2,...,16;
③从16个小波包系数中剔除掉互相关系数绝对值小于0.9对应的小波包系数,后对剩余小波包系数的重构序列叠加得到第s段静风速值和静位移值的重构序列,按此方法将得到的p段静风速值和静位移值的重构序列重新组合成第q个月份的静风速和静位移序列,以此作为静风速和静位移在此月的主相关序列,后对各月遍历得到所有月份静风速和静位移之间的主相关序列; ③ Eliminate the wavelet packet coefficients corresponding to the absolute value of the cross-correlation coefficient less than 0.9 from the 16 wavelet packet coefficients, and then superimpose the reconstruction sequence of the remaining wavelet packet coefficients to obtain the reconstruction sequence of the static wind speed value and static displacement value of the s segment, According to this method, the reconstructed sequence of the static wind speed and static displacement of the p segment is recombined into the static wind speed and static displacement sequence of the qth month, which is used as the main correlation sequence of the static wind speed and static displacement in this month, Afterwards, the main correlation sequence between the static wind speed and static displacement of all months is obtained by traversing each month;
步骤(4):依次对各月主相关序列进行傅里叶级数拟合,利用拟合参数值的月变化特征对大跨桥梁结构的整体抗侧力性能进行异常预警: Step (4): Carry out Fourier series fitting on the main correlation sequence of each month in turn, and use the monthly variation characteristics of the fitted parameter values to give an abnormal early warning of the overall lateral force resistance performance of the long-span bridge structure:
利用下式所示的二阶傅里叶级数,依次对各月主相关序列进行最小二乘拟合并确定各月参数估计值: Using the second-order Fourier series shown in the following formula, the least squares fitting is performed on the main correlation series of each month in turn and the parameter estimates of each month are determined:
式中,vm表示第m个月主相关序列中的静风速值,um表示第m个月主相关序列中的静位移值,ae(m)、bx(m)和w(m)分别为第m个月傅里叶级数的参数估计值;分别对各参数估计值的月变化特性进行ADF单位根检验,若每个参数估计值的检验结果均拒绝存在一个单位根的原假设,则大跨桥梁结构的整体抗侧力性能处于良好状态;若存在1至2个参数估计值的检验结果不拒绝存在一个单位根的原假设,则对大跨桥梁结构的整体抗侧力性能进行黄色预警,并密切跟踪各参数估计值月变化特征的ADF单位根检验结果;若存在3至4个参数估计值的检验结果不拒绝存在一个单位根的原假设,则对大跨桥梁结构的整体抗侧力性能进行橙色预警,并派桥梁检修人员对桥梁结构关键构件和部位进行现场检测,根据检测结果采取应对措施;若存在4个以上参数估计值的检验结果不拒绝存在一个单位根的原假设,则对大跨桥梁结构的整体抗侧力性能进行红色预警,并派桥梁专业人员对桥梁结构的整体抗侧力性能进行安全评估和决策; In the formula, v m represents the static wind speed value in the mth month main correlation sequence, u m represents the static displacement value in the mth month main correlation sequence, a e (m), b x (m) and w (m ) are the parameter estimates of the Fourier series in the mth month; the ADF unit root test is performed on the monthly variation characteristics of each parameter estimate, and if the test results of each parameter estimate reject the existence of a unit root Assuming that , the overall lateral force resistance of long-span bridge structures is in good condition; A yellow warning will be issued for the performance, and the ADF unit root test results of the monthly variation characteristics of the estimated values of each parameter will be closely tracked; The overall lateral force resistance performance will be given an orange warning, and bridge maintenance personnel will be sent to conduct on-site inspections of key components and parts of the bridge structure, and countermeasures will be taken according to the inspection results; if there are more than 4 parameter estimates, the inspection results do not reject the existence of a unit root The null hypothesis of the long-span bridge structure is given a red warning for the overall lateral force resistance performance, and bridge professionals are sent to conduct safety assessment and decision-making on the overall lateral force resistance performance of the bridge structure;
实施例1 Example 1
下面以苏通大桥为分析对象,说明本发明的具体实施过程: Take Sutong Bridge as the object of analysis below to illustrate the concrete implementation process of the present invention:
(1)苏通大桥是连接南通与镇江两市的跨长江大桥,采用双塔双索面斜拉桥结构体系,其中主梁构件采用流线型扁平钢箱梁形式,主跨部分纵向设计尺寸达到1088m,这一设计尺寸使得主梁构件在横桥向风载作用下,跨中部位会出现明显的侧向位移效应。基于桥梁结构健康监测系统,对主梁跨中部位的三维风场和GPS位移响应进行长期监测和数据采集,具体监测仪器布置如图1所示,仪器采样频率均设定为1Hz; (1) Sutong Bridge is a bridge spanning the Yangtze River connecting Nantong and Zhenjiang. It adopts a double-tower double-cable-plane cable-stayed bridge structure system, in which the main girder is in the form of a streamlined flat steel box girder, and the longitudinal design dimension of the main span reaches 1088m , this design dimension makes the main girder members have obvious lateral displacement effect at the mid-span under the action of wind load in the transverse bridge direction. Based on the bridge structural health monitoring system, long-term monitoring and data collection are carried out for the three-dimensional wind field and GPS displacement response at the mid-span of the main girder. The specific monitoring instrument layout is shown in Figure 1, and the sampling frequency of the instruments is set at 1 Hz;
(2)基于步骤2)对采集数据进行矢量分解和均值处理,得到横桥向静风速序列和静位移序列分别如图2和图3所示(以2012年8月1日至8月10日为例); (2) Based on step 2), vector decomposition and mean value processing are performed on the collected data, and the static wind speed sequence and static displacement sequence of the cross bridge are obtained as shown in Figure 2 and Figure 3 respectively (taken from August 1 to August 10, 2012 example);
(3)基于步骤3)对对静风速序列和静位移序列按月份进行划分,以8月份的划分结果为分析例,将此月静风速值按递增排列并等间距划分为11段,此外按照R(n1,n2)的排列规则将此月静位移值排列并同样均分为11段,11段划分结果采用月静风速值与月静位移值之间的相关性散点图表示如图4所示; (3) Based on step 3), the static wind speed series and static displacement series are divided by month. Taking the division results in August as an analysis example, the monthly static wind speed values are arranged in increasing order and divided into 11 segments at equal intervals. In addition, according to The arrangement rule of R(n 1 ,n 2 ) arranges the monthly static displacement values and equally divides them into 11 segments. The results of the division of 11 segments are represented by the correlation scatter diagram between the monthly static wind speed value and the monthly static displacement value as follows: As shown in Figure 4;
(4)基于步骤3)分别对每段静风速值和静位移值进行第4尺度小波包分解,每段静风速值和静位移值均得到按结点位置排列的16个小波包系数,对相同结点位置小波包系数的重构序列进行逐一互相关分析,剔除掉互相关系数绝对值小于0.9的小波包系数,并对剩余小波包系数重构得到每段静风速值和静位移值的重构序列,将11段重构序列组合得到静风速和静位移在此月的主相关序列如图5所示; (4) Based on step 3), the fourth-scale wavelet packet decomposition is performed on each static wind speed value and static displacement value, and each static wind speed value and static displacement value can obtain 16 wavelet packet coefficients arranged according to the node position. The reconstruction sequence of wavelet packet coefficients at the same node position is analyzed one by one by cross-correlation, and the wavelet packet coefficients whose absolute value of cross-correlation coefficient is less than 0.9 are eliminated, and the remaining wavelet packet coefficients are reconstructed to obtain the static wind speed value and static displacement value of each segment. Reconstruction sequence, combining 11 segments of reconstruction sequence to obtain the main correlation sequence of static wind speed and static displacement in this month, as shown in Figure 5;
(5)基于步骤4)对各月主相关序列进行二阶傅里叶级数拟合,各拟合参数值的月变化特征如图6-8所示,对各参数估计值的月变化特性进行ADF单位根检验,存在2个参数估计值的检验结果不拒绝存在一个单位根的原假设,则对大跨桥梁结构的整体抗侧力性能进行黄色预警,并密切跟踪各参数估计值月变化特征的ADF单位根检验结果。 (5) Based on step 4), the second-order Fourier series fitting is performed on each monthly main correlation sequence. The monthly variation characteristics of each fitting parameter value are shown in Figure 6-8. The monthly variation characteristics of each parameter estimated value Carry out the ADF unit root test, and if the test results of two parameter estimates do not reject the null hypothesis of the existence of a unit root, a yellow warning will be given to the overall lateral force resistance of the long-span bridge structure, and the monthly changes in the estimated values of each parameter will be closely tracked The ADF unit root test results of the features.
Claims (3)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410054835.5A CN103837358B (en) | 2014-02-18 | 2014-02-18 | The method for early warning of the overall lateral resistance behavior exception of long-span bridges |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410054835.5A CN103837358B (en) | 2014-02-18 | 2014-02-18 | The method for early warning of the overall lateral resistance behavior exception of long-span bridges |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103837358A CN103837358A (en) | 2014-06-04 |
CN103837358B true CN103837358B (en) | 2016-03-02 |
Family
ID=50801076
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410054835.5A Expired - Fee Related CN103837358B (en) | 2014-02-18 | 2014-02-18 | The method for early warning of the overall lateral resistance behavior exception of long-span bridges |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103837358B (en) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104048843B (en) * | 2014-06-13 | 2017-10-10 | 东南大学 | Loads of Long-span Bridges steel box-girder damage alarming method based on GPS displacement monitorings |
CN104122111B (en) * | 2014-08-12 | 2016-11-23 | 招商局重庆交通科研设计院有限公司 | A kind of three class comprehensive safety method for early warning of bridge structure |
CN104820234B (en) * | 2015-04-13 | 2017-04-12 | 东南大学 | Quick identification method for earthquake load of bridge |
CN106568557B (en) * | 2016-10-31 | 2018-08-21 | 东南大学 | A kind of safe early warning method of high-speed railway bridge vehicle-bridge vibration performance |
CN106570299A (en) * | 2016-11-14 | 2017-04-19 | 东南大学 | Method for determining vehicle-bridge resonance performance curves of high-speed railway steel truss arch bridges |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102565194A (en) * | 2012-02-09 | 2012-07-11 | 东南大学 | Method for carrying out early warning on damage to steel box girder of long span bridge in operation state |
CN102567630A (en) * | 2011-12-20 | 2012-07-11 | 东南大学 | Method for determining wind-induced vibrating response of long-span bridge structure |
CN103440404A (en) * | 2013-08-07 | 2013-12-11 | 东南大学 | Lateral force resisting performance degradation alarm method for bridge stiffening girder based on transverse wind load effect |
CN103530521A (en) * | 2013-10-22 | 2014-01-22 | 东南大学 | Sunlight temperature time interval simulation method based on Fourier series and ARMA model |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3212443B2 (en) * | 1994-05-20 | 2001-09-25 | 三菱重工業株式会社 | Seismic strong wind observation method for structures |
US9581570B2 (en) * | 2011-02-10 | 2017-02-28 | University Of South Carolina | Determination of the remaining life of a structural system based on acoustic emission signals |
-
2014
- 2014-02-18 CN CN201410054835.5A patent/CN103837358B/en not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102567630A (en) * | 2011-12-20 | 2012-07-11 | 东南大学 | Method for determining wind-induced vibrating response of long-span bridge structure |
CN102565194A (en) * | 2012-02-09 | 2012-07-11 | 东南大学 | Method for carrying out early warning on damage to steel box girder of long span bridge in operation state |
CN103440404A (en) * | 2013-08-07 | 2013-12-11 | 东南大学 | Lateral force resisting performance degradation alarm method for bridge stiffening girder based on transverse wind load effect |
CN103530521A (en) * | 2013-10-22 | 2014-01-22 | 东南大学 | Sunlight temperature time interval simulation method based on Fourier series and ARMA model |
Non-Patent Citations (1)
Title |
---|
"桥梁健康监测海量数据分析与评估——"结构健康监测"研究进展";李爱群 等;《中国科学:技术科学》;20120820;第42卷(第8期);第972-984页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103837358A (en) | 2014-06-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106197910B (en) | Bridge detection method and detection system based on vehicle-bridge coupling vibration analysis | |
CN103837358B (en) | The method for early warning of the overall lateral resistance behavior exception of long-span bridges | |
Wahab et al. | Damage detection in bridges using modal curvatures: application to a real damage scenario | |
CN105760934A (en) | Bridge abnormity monitoring restoration method based on wavelet and BP neural network | |
CN105241660A (en) | High-speed rail large-scale bridge performance evaluation method based on health monitoring data | |
CN106897543A (en) | The girder construction damnification recognition method of On Modal Flexibility Curvature matrix norm | |
CN103984875A (en) | Bridge structure damage diagnosis method based on accumulated damage factor sequence under complex environment | |
CN104268658A (en) | Bridge structure safety monitoring data prediction method | |
Martín-Sanz et al. | Monitoring of the UHPFRC strengthened Chillon viaduct under environmental and operational variability | |
Alves et al. | A fast and efficient feature extraction methodology for structural damage localization based on raw acceleration measurements | |
CN104915534A (en) | Deformation analysis and decision-making method of electric power tower based on sequence learning | |
CN112816352A (en) | Engineering structure damage identification method and device, computer equipment and storage medium | |
Chen et al. | Damage detection of long-span bridges using stress influence lines incorporated control charts | |
Amanollah et al. | Damage detection of structures based on wavelet analysis using improved AlexNet | |
KR101579732B1 (en) | A method for novel health monitoring scheme for smart concrete structures | |
Ahmadi et al. | New damage index based on least squares distance for damage diagnosis in steel girder of bridge's deck | |
CN112326787A (en) | Beam bridge identification method based on multipoint rapid static acquisition of exclusive test car | |
Morichika et al. | Estimation of displacement response in steel plate girder bridge using a single MEMS accelerometer | |
CN102998133A (en) | Energy damage identification method based on quasi-distributed acceleration data | |
CN103913512A (en) | Damage positioning system and damage positioning method for stay cable periodic detection | |
Singh et al. | Recent trends in operation modal analysis techniques and its application on a steel truss bridge | |
CN106446384A (en) | A damage identification method for the main girder structure of bridge crane | |
CN110008520A (en) | Structural Damage Identification based on dynamic respond covariance parameter and Bayesian Fusion | |
CN110657882B (en) | A real-time monitoring method of bridge safety state using single measuring point response | |
CN103911958B (en) | The damage reason location system of suspension bridge and arch bridge suspender periodic detection and method thereof |
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: 20160302 Termination date: 20190218 |
|
CF01 | Termination of patent right due to non-payment of annual fee |