CN116882238A - 一种考虑边界扰动的拉索索力和弯曲刚度联合识别方法 - Google Patents
一种考虑边界扰动的拉索索力和弯曲刚度联合识别方法 Download PDFInfo
- Publication number
- CN116882238A CN116882238A CN202310837485.9A CN202310837485A CN116882238A CN 116882238 A CN116882238 A CN 116882238A CN 202310837485 A CN202310837485 A CN 202310837485A CN 116882238 A CN116882238 A CN 116882238A
- Authority
- CN
- China
- Prior art keywords
- cable
- boundary
- vibration
- inhaul
- force
- 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
- 238000005452 bending Methods 0.000 title claims abstract description 45
- 238000000034 method Methods 0.000 title claims abstract description 34
- 230000004044 response Effects 0.000 claims abstract description 25
- 238000012546 transfer Methods 0.000 claims abstract description 24
- 238000004873 anchoring Methods 0.000 claims abstract description 15
- 238000012360 testing method Methods 0.000 claims description 6
- 230000001133 acceleration Effects 0.000 claims description 3
- 238000005311 autocorrelation function Methods 0.000 claims description 2
- 238000005314 correlation function Methods 0.000 claims description 2
- 238000006073 displacement reaction Methods 0.000 claims description 2
- 239000000463 material Substances 0.000 claims description 2
- 239000011159 matrix material Substances 0.000 claims description 2
- YFJYSJRZDOWXDH-UHFFFAOYSA-M sodium 4-[4-(2-cyanoethynyl)benzoyl]oxy-2,3,5,6-tetrafluorobenzenesulfonate Chemical compound [Na+].[O-]S(=O)(=O)c1c(F)c(F)c(OC(=O)c2ccc(cc2)C#CC#N)c(F)c1F YFJYSJRZDOWXDH-UHFFFAOYSA-M 0.000 claims 6
- 238000012544 monitoring process Methods 0.000 abstract description 4
- 230000036541 health Effects 0.000 abstract description 3
- 238000007405 data analysis Methods 0.000 abstract description 2
- 238000007689 inspection Methods 0.000 abstract description 2
- 230000001360 synchronised effect Effects 0.000 abstract 1
- 238000004458 analytical method Methods 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 229910000831 Steel Inorganic materials 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 239000010959 steel Substances 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005316 response function Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000005481 string theory Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
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
- G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
- G01D21/00—Measuring or testing not otherwise provided for
- G01D21/02—Measuring two or more variables by means not covered by a single other subclass
-
- 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/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Architecture (AREA)
- Civil Engineering (AREA)
- Structural Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Bridges Or Land Bridges (AREA)
Abstract
本发明属于桥梁结构健康检、监测数据分析技术领域,提供了一种考虑边界扰动的拉索索力和弯曲刚度联合识别方法。本发明只需要在斜拉索梁端锚固处和索上某一位置处布置两个测点,由桥梁运营状态下随机车辆激发的索‑梁振动响应得到索‑梁传递函数,将拉索的各阶频率作为目标参数实现拉索索力和弯曲刚度的同步识别。该方法不仅考虑了拉索边界扰动的影响,而且在响应信号中消除了梁振动对索振动的影响,精确度高,可操作性强,保证了其在实际中具备更广泛的应用。
Description
技术领域
本发明属于桥梁结构健康检、监测数据分析技术领域,涉及一种考虑边界扰动的拉索索力和弯曲刚度联合识别方法。
背景技术
随着高强度钢丝的出现,拉索已成为桥梁结构中的有效张拉构件,其广泛应用于斜拉桥等索承桥梁结构中。识别索力为其施工控制、实时监测、损伤识别、异常预警和性能评估提供有效指标。因此,索力识别是近年来结构健康监测中的重要课题。
基于振动频率法的索力识别方法越来越得到工程技术人员的重视,因其便于安装、动态响应好、测试设备可重复使用等优点,已成为目前最有优势、使用范围最广的索力识别方法。通过振动模态、频响函数等动力学信息可识别拉索的索力和刚度等参数的变化,在建立索力与模态参数关系的过程中仍然面临弯曲刚度、垂度、倾角和复杂多变的边界条件等因素的挑战。
基于张紧弦理论的索力识别公式在工程测试中得到了广泛的应用,因其忽略了各种因素使得适用范围受到限制。考虑弯曲刚度的拉索固有频率解析解在各种边界条件下得到了广泛研究,Zui等人在文章“Practical formulas for estimation of cable tensionby vibration method”中给出了分段的索力识别实用公式,任伟新等人在文章“Empiricalformulas to estimate cable tension by cable fundamental frequency”中分别提出了考虑垂度和弯曲刚度的经验公式,骆兰迎等人在文章“Piecewise-fitted formula forcable force identification considering bending stiffness,sag,and inclination”中构造出考虑弯曲刚度、垂度和倾角的分段拟合公式(PFF)并给出不同情况下的使用系数。上述推导的索力识别公式在使用中忽略了弯曲刚度准确性的问题,实际的缆索是由多根钢丝并行绑扎在一起并且可以互相滑动,因此拉索真实的弯曲刚度很难计算,这将影响估算的索力的精度。此外,上述索力识别公式都是基于某一理想边界下推导而来,实际桥梁运营下车致主梁振动使得拉索的边界产生扰动,此时亟需一种可以考虑边界扰动的索力识别方法。
针对上述弯曲刚度估算不准,边界扰动影响拉索振动体系两大难题,提出了一种考虑边界扰动的拉索索力和弯曲刚度联合识别方法,通过车致主梁和拉索的振动,建立拉索梁端锚固处和索上某一位置处之间振动响应的传递关系,提出了基于索-梁传递函数的索力和弯曲刚度联合识别算法。其主要特点包括:1)利用随机车辆激励桥梁振动,获取拉索梁端锚固处和索上某一位置的振动响应,基于索-梁传递函数联合识别索力和弯曲刚度;2)所提方法考虑了拉索边界扰动的因素,识别方法避免了传统的弯曲刚度估算不准对索力识别精度的影响。
发明内容
本发明的目的是通过车致主梁和拉索的振动响应,提供了一种考虑边界扰动的拉索索力和弯曲刚度联合识别方法。
本发明的技术方案:
一种考虑边界扰动的拉索索力和弯曲刚度联合识别方法,步骤如下:
步骤一:基于车辆诱导的索-梁振动响应得到实测索-梁传递函数
桥梁运营状态下随机车辆诱导主梁和拉索振动,将两个加速度传感器分别安装在拉索梁端锚固处(即x=l)和索上某一位置处(即x=lc),通过实测振动响应得到边界扰动下的索-梁传递函数(CBTF):
式中:GXF(ω)、GFF(ω)分别为v(l,t)与v(lc,t)之间的互相关函数和自相关函数;v(l,t)为拉索在梁端锚固处即x=l的边界振动响应,v(lc,t)为索上某一位置即x=lc处的振动响应。从实测CBTF中识别拉索的第k到n阶固有频率。
步骤二:从考虑边界扰动的模型中得到理论的索-梁传递函数
(1)建立考虑边界扰动的拉索理论模型及坐标系Oxy,如图1所示;将拉索塔端边界假定为固结,将拉索梁端锚固处的振动响应vimp(t)作为已知量,得到考虑边界扰动的拉索振动方程及边界条件
v(0,t)=0
v′(0,t)=0
v(l,t)=vimp(t)
v′(l,t)=0
式中:T、l、E、I分别为索力、索长、拉索单位长度质量、材料弹性模量、截面惯性矩;v(x,t)为拉索面内振动的横向位移,v(0,t)为拉索在塔端锚固处(即x=0)的边界条件,v(l,t)为拉索在梁端锚固处(即x=l)的边界振动响应。
(2)将上述拉索振动方程及边界条件转换到频域进行求解,得到形函数在频域的解:
式中:为拉索振动形函数在频域的解的形式,ω为频域的坐标轴单位,向量e={exp(k1x)exp[k2(x-l)]exp(k3x)exp(k4x)};/>由边界条件确定,λ2=T/2EI,/>
(3)将已知边界条件代入上述形函数解的形式,得到由边界扰动引起的拉索任意一点x=lc处响应在频域的表示形式:
式中:特征矩阵
(4)定义索-梁传递函数(CBTF)为索上(即x=lc)响应与拉索梁端锚固处(即x=l)响应之比:
步骤三:利用索-梁传递函数联合识别索力和弯曲刚度
将实测CBTF作为真实值,将理论的CBTF中的索力和弯曲刚度作为待识别参数,如图2所示。建立目标函数:
式中:aj为第j阶的权重因子,默认为1。fj,test、fj,CBTF分别为测试和待确定的传递函数的第j阶频率,p代表待识别的索参数(索力和弯曲刚度);n、k分别为实测拉索第k到n阶固有频率,可测得的频率总数为n-k+1。给定索力和弯曲刚度初始值,使目标函数最小时即为拉索真实的索力和弯曲刚度。
本发明的有益效果:通过建立考虑边界扰动的拉索理论模型,可准确得到包含拉索固有特性的索-梁传递函数,从而基于实测得到拉索的索-梁传递函数,联合识别索力和弯曲刚度。该方法避免了弯曲刚度等效计算不准影响索力识别的精度,仅通过拉索梁端锚固处和索上某点响应即可实现索力识别,使得拉索索力测试更加准确高效。
附图说明
图1为本发明方法建立的考虑边界扰动的拉索理论模型。
图2为本发明方法提出的基于索-梁传递函数联合识别索力和弯曲刚度的过程示意图。
图3为本发明方法实施例中某斜拉桥模型及拉索测点布置图。
图4为本发明方法实施例中所选6根拉索的索-梁传递函数测试结果:(a)SE01;(b)SE03;(c)SE05;(d)SE07,(e)SE09;and(f)SE11。
图5为本发明的方法流程图。
具体实施方式
下面结合附图和一个数值算例来对本发明做进一步的详细说明。
本发明的识别方法分“基于车辆诱导的索-梁振动响应得到实测索-梁传递函数”、“从考虑边界扰动的模型中得到理论的索-梁传递函数”和“利用索-梁传递函数联合识别索力和弯曲刚度”三步,具体实施方式上文已经给出,接下来结合一座双塔斜拉桥的有限元模型作为数值算例对本发明方法进行可行性验证和实施说明。
如图3所示,采用ANSYS-18.2软件对某斜拉桥有限元模型进行仿真分析,主桥跨径260m,两侧跨径125.0m,主梁和桥塔采用梁单元BEAM44进行建模,主梁隔板采用MASS21进行建模。为了研究考虑抗弯刚度的斜拉索模型,采用梁单元BEAM188对斜拉索进行了建模,解除BEAM188两端的弯矩约束,使斜拉索只对主塔和主梁提供拉力。取拉索编号SE01、SE03、SE05、SE07、SE09和SE11进行研究分析。由于大跨度斜拉桥的车桥耦合和车轮荷载分布可以忽略,用集中荷载模拟车辆荷载。为模拟索的振动是由车致主梁振动引起的,本研究采用多车随机过桥的激励方法。本发明方法包括以下步骤:
提取每根拉索梁端锚固处和索上某一点的振动响应,根据步骤一得到实测的索-梁传递函数,从实测CBTF中识别拉索的第k到n阶固有频率。
根据所提考虑边界扰动的拉索理论模型,建立拉索理论的索-梁传递函数。
建立以拉索固有频率为参数的目标函数,结合所提出的用索-梁传递函数联合识别索力和弯曲刚度的算法,给定待识别的索力和弯曲刚度初始值,通过算法识别使得目标函数最小时的索力和弯曲刚度即为真实值。
通过识别结果与有限元模型成桥下的索力和弯曲刚度值对比,结果如表1所示,得出现有索力识别公式的识别误差超过了5%,本发明方法对索力和弯曲刚度的识别误差在5%以内,验证了所提方法的可行性。
从计算结果可以看出,本发明提出的考虑边界扰动的拉索索力和弯曲刚度联合识别方法可以得到更精确的识别结果。该方法避免了目前弯曲刚度等效计算不准影响索力识别的精度,并且通过两个测点即可实现索力和弯曲刚度的同时识别,使得拉索索力识别更加高效准确。
表1CBTF与现有索力识别公式的识别结果比较
注:IF为识别的索力;II为识别的弯曲刚度;RE为识别值与理论值之间的相对误差。
Claims (1)
1.一种考虑边界扰动的拉索索力和弯曲刚度联合识别方法,步骤如下:
步骤一:基于车辆诱导的索-梁振动响应得到实测索-梁传递函数
桥梁运营状态下随机车辆诱导主梁和拉索振动,将两个加速度传感器分别安装在拉索梁端锚固处即x=l和索上某一位置处即x=lc,通过实测振动响应得到边界扰动下的索-梁传递函数CBTF:
式中:GXF(ω)、GFF(ω)分别为v(l,t)与v(lc,t)之间的互相关函数和自相关函数;v(l,t)为拉索在梁端锚固处即x=l的边界振动响应,v(lc,t)为索上某一位置即x=lc处的振动响应;从实测CBTF中识别拉索的第k到n阶固有频率;
步骤二:从考虑边界扰动的模型中得到理论的索-梁传递函数
(1)建立考虑边界扰动的拉索理论模型及坐标系Oxy;将拉索塔端边界假定为固结,将拉索梁端锚固处的振动响应vimp(t)作为已知量,得到考虑边界扰动的拉索振动方程及边界条件:
v(0,t)=0
v′(0,t)=0
v(l,t)=vimp(t)
v′(l,t)=0
式中:T、l、E、I分别为索力、索长、拉索单位长度质量、材料弹性模量、截面惯性矩;v(x,t)为拉索面内振动的横向位移,v(0,t)为拉索在塔端锚固处即x=0的边界条件;
(2)将上述拉索振动方程及边界条件转换到频域进行求解,得到形函数在频域的解:
式中:为拉索振动形函数在频域的解的形式,ω为频域的坐标轴单位,向量e={exp(k1x) exp[k2(x-l)] exp(k3x) exp(k4x)};/>由边界条件确定,λ2=T/2EI,/>
(3)将已知边界条件代入上述形函数解的形式,得到由边界扰动引起的拉索任意一点x=lc处响应在频域的表示形式:
式中:特征矩阵
(4)定义索-梁传递函数CBTF为索上即x=lc响应与拉索梁端锚固处即x=l响应之比:
步骤三:利用索-梁传递函数联合识别索力和弯曲刚度
将实测CBTF作为真实值,将理论的CBTF中的索力和弯曲刚度作为待识别参数;建立目标函数:
式中:aj为第j阶的权重因子,默认为1;fj,test、fj,CBTF分别为测试和待确定的传递函数的第j阶频率;p代表待识别的索参数,包括索力和弯曲刚度;n、k分别为实测拉索第k到n阶固有频率,测得的频率总数为n-k+1;给定索力和弯曲刚度初始值,使目标函数最小时即为拉索真实的索力和弯曲刚度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310837485.9A CN116882238A (zh) | 2023-07-10 | 2023-07-10 | 一种考虑边界扰动的拉索索力和弯曲刚度联合识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310837485.9A CN116882238A (zh) | 2023-07-10 | 2023-07-10 | 一种考虑边界扰动的拉索索力和弯曲刚度联合识别方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116882238A true CN116882238A (zh) | 2023-10-13 |
Family
ID=88263752
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310837485.9A Pending CN116882238A (zh) | 2023-07-10 | 2023-07-10 | 一种考虑边界扰动的拉索索力和弯曲刚度联合识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116882238A (zh) |
-
2023
- 2023-07-10 CN CN202310837485.9A patent/CN116882238A/zh active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Gul et al. | Ambient vibration data analysis for structural identification and global condition assessment | |
Wang et al. | Identification of moving vehicle parameters using bridge responses and estimated bridge pavement roughness | |
Law et al. | Time-varying wind load identification from structural responses | |
CN110285909A (zh) | 基于同步压缩变换的索承桥梁瞬时索力计算方法 | |
CN109839441A (zh) | 一种桥梁模态参数识别方法 | |
JP3313028B2 (ja) | 張力のかかっているケーブルの曲げ剛性及び張力の測定方法 | |
CN106153176B (zh) | 一种用于对单相油浸式变压器噪声进行预估的系统及方法 | |
He et al. | FEM free damage detection of beam structures using the deflections estimated by modal flexibility matrix | |
KR100972278B1 (ko) | 구조물의 측정 진동가속도를 이용하여 동적 변위 및 속도이력을 측정하는 측정방법 및 측정장치 | |
Zhang et al. | Deflection estimation of beam structures based on the measured strain mode shape | |
Nie et al. | Output-only complete mode shape identification of bridges using a limited number of sensors | |
Benčat et al. | Bridges subjected to dynamic loading | |
Liu et al. | Bayesian inference based parametric identification of vortex-excited force using on-site measured vibration data on a long-span bridge | |
Nie et al. | Baseline-free structural damage detection using PCA-Hilbert transform with limited sensors | |
Kwon | Uncertainty of bridge flutter velocity measured at wind tunnel tests | |
Tomaszewska et al. | Study on applicability of two modal identification techniques in irrelevant cases | |
Gang et al. | Aerodynamic admittance influence on buffeting performance of suspension bridge with streamlined deck | |
Chen et al. | Experimental modal test and time-domain aerodynamic analysis of a cable-stayed bridge | |
CN116882238A (zh) | 一种考虑边界扰动的拉索索力和弯曲刚度联合识别方法 | |
CN112326787A (zh) | 一种基于专属测试车多点快速静态采集的梁式桥识别方法 | |
Liu et al. | Measurement and modelling of the vibration induced by working equipment on an offshore platform | |
Kangas et al. | Identification of cable forces on cable-stayed bridges: A novel application of the MUSIC algorithm | |
Ganguly et al. | Performance assessment of time-domain damage indicators based on output-only measurement and Poincaré map: A comparative review on nonlinear structures | |
Park et al. | Reconstruction of displacement, velocity and acceleration from measured dynamic strain for Bernoulli-beam type girders of bridges | |
Manolis et al. | Experimental evaluation of damping in beams using the acceleration generalized coordinates: A comparison of the FDD and PCA methods |
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 |