CN109063338B - 基于完全正交化算法识别预应力桥梁现存预应力的方法 - Google Patents
基于完全正交化算法识别预应力桥梁现存预应力的方法 Download PDFInfo
- Publication number
- CN109063338B CN109063338B CN201810879043.XA CN201810879043A CN109063338B CN 109063338 B CN109063338 B CN 109063338B CN 201810879043 A CN201810879043 A CN 201810879043A CN 109063338 B CN109063338 B CN 109063338B
- Authority
- CN
- China
- Prior art keywords
- bridge
- equation
- matrix
- prestress
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000006073 displacement reaction Methods 0.000 claims abstract description 70
- 239000011159 matrix material Substances 0.000 claims description 70
- 238000013016 damping Methods 0.000 claims description 9
- 238000005452 bending Methods 0.000 claims description 3
- 230000010354 integration Effects 0.000 claims description 3
- 239000000463 material Substances 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 2
- 238000005259 measurement Methods 0.000 claims description 2
- 238000013507 mapping Methods 0.000 abstract 1
- 230000006378 damage Effects 0.000 description 3
- 230000002787 reinforcement Effects 0.000 description 2
- 208000027418 Wounds and injury Diseases 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 208000014674 injury Diseases 0.000 description 1
- 238000009659 non-destructive testing Methods 0.000 description 1
Images
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/13—Differential equations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- 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/06—Power analysis or power optimisation
-
- 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
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Operations Research (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Computing Systems (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
Abstract
Description
技术领域
本发明属于预应力桥梁现存预应力的识别技术领域,尤其涉及一种基于完全正交化算法的识别预应力桥梁现存预应力的方法。
背景技术
随着我国经济的迅猛发展,公路铁路交通量均日益增加,且在今后较长的一段时间内我国交通运输需求将继续保持快速增长趋势。货物运输中大宗货物车辆荷载重型化、车速高速化的发展虽然创造了良好的经济效益,但也对我国桥梁等交通基础设施造成了潜在的危害,因此对桥梁结构的安全性、适用性与耐久性提出更高要求。如何精确评估现役桥梁的结构安全状况,避免由于桥梁结构破坏导致的人员伤害和财产损失是桥梁工程人员迫切需要解决的工程难点。预应力桥梁现存预应力的大小是评估桥梁抗裂性能和承载能力的重要指标,采用无损检测技术快速识别预应力桥梁现存预应力具有良好的工程可行性。
发明内容
本发明的目的是提供一种基于完全正交化算法由桥梁位移实时识别预应力桥梁现存预应力的方法,识别速度快且精度较高。
为达到上述目的,本发明采用的技术方案是:一种基于完全正交化算法的识别预应力桥梁现存预应力的方法,包括以下步骤:
1)、在桥梁底面从左端依次向右布置有m个位移传感器;
3)、由测得的桥梁竖向位移计算预应力桥梁的振幅{q}N×1;
4)、建立识别桥梁现存预应力值的系统方程T=({S}T{S})-1{S}T{b};
5)、利用基于完全正交化算法,由系统方程初始解得到残差,由残差构造得到w维Krylov子空间,求得桥梁现存预应力值T。
所述步骤2)中具体包括以下步骤:
21)、建立桥梁的振动微分方程:
将桥梁考虑为一单跨简支梁,假设固定荷载P(t)位于单跨简支梁上,其距离桥梁左端距离为xp,通过在桥梁底面粘贴的位移传感器可测得距桥梁左端x位置处t时刻的动态位移为y(x,t),则其振动微分方程:
其中单跨简支梁的跨长为L,密度为ρ,截面面积为A,粘性比例阻尼为c;桥梁抗弯刚度为EI,E是材料的弹性模量,I=bh0 3/12是梁截面的惯性矩,其中b是桥梁横截面宽度,h0是桥梁横截面高度;δ(x-xp)是狄拉克函数;
22)、得到简支梁的模态振型函数的矩阵形式:
基于模态叠加原理,假设梁的第i阶模态振型函数为Yi,梁的第i阶振幅为qi(t),则动态位移y(x,t)的模态形式表示为:
把方程(2)带入方程(1)并对方程两端的每一项均乘以Yi(x),考虑模态正交条件并对方程中x值从0到L进行积分,得到新方程:
简支预应力梁的模态振型函数标准化形式表示为:
在时域内通过Newmark积分将方程(3)用矩阵形式表示为:
假定在简支桥梁底面从左端支座依次向右布置有m个位移传感器,则在模态坐标下距简支梁左端支座xm处的第m个位移传感器实测竖向位移y(xm,t)表示为:
其中,Yi为梁的第i阶模态振型函数,qi(t)为梁的第i阶振幅,方程(6)用矩阵形式表示为:
则广义坐标下基于最小二乘法的振幅伪逆解表示为:
布置在桥梁底面的第j个位移传感器在t时刻测得的竖向位移向量y(xj,t)由广义正交多项式拟合为:
其中,Nf为正交多项式函数的阶数,Gi为第i阶正交多项式,ai为第i阶正交多项式的系数;
将方程(9)带入方程(7)并写成矩阵形式得:
式中,[A],[G]分别为正交多项式的系数矩阵和正交多项式矩阵;正交多项式的系数矩阵[A]可通过方程(10)由最小二乘法得到:
所述的步骤4)具体为:
所述的方程(5)重写为:
式中,矩阵[K′]包含待识别的桥梁预应力值T,矩阵[C]包含桥梁模态阻尼ξi和模态频率ωi,矩阵[K]包含预应力桥梁的系统参数;
在时域内求解方程(12)即可得到桥梁现存预应力值T,方程(12)在时域内可表示为:
其中
则待识别预应力T由最小二乘法直接解得:
由方程(14)即可求得在t时刻预应力桥梁现存预应力值T(t),当取固定荷载P(t)在t1,t2,…,tw时刻的荷载值分别为P(t1),P(t2),…,P(tw),则方程(12)在时域内可重写为:
{S}N×wTw×1={b}N×1 (15)
其中
则待识别预应力T由最小二乘法直接解得:
T=({S}T{S})-1{S}T{b} (16)
所述的步骤5)具体包括以下步骤:
取基于桥梁动力响应识别现存预应力值系统矩阵为{S}T{S},由系统方程初始解T0得到的残差为:
r0={S}T{b}-{S}T{S}T0 (17)
取由r0构造得到的w维Krylov子空间为:
记β=‖r0‖2,则经过w步Arnoldi算法可构造一组标准正交基Vw=[v1,v2,…vw]和一个(w+1)×w的Hessenberg矩阵Hw+1,将Hessenberg矩阵Hw的最后一行元素删除即可得到大小为w×w的Hessenberg矩阵Hw,识别现存预应力值系统矩阵和标准正交基Vw以及Hessenberg矩阵Hw有下面的关系:
完全正交化算法采用正交投影理论,由标准正交基Vw和Hessenberg矩阵Hw即可在Krylov子空间解得预应力桥梁现存预应力值Tw,其表达式为:
其中,β=‖r0‖2,e1为第一个元素为1的w维单位向量,由公式(20)即可求得待识别的桥梁现存预应力值即为Tw。
本发明具有的优点是:本发明提出一种基于完全正交化算法由桥梁位移实时识别预应力桥梁现存预应力的方法,通过建立预应力桥梁现存预应力与桥梁竖向动态位移响应之间的力学映射关系,实现由桥梁位移响应识别桥梁现存预应力值。为提高识别方法的迭代效率,采用完全正交化算法以节约识别时间,最终实现实时识别桥梁现存预应力值,具有迭代收敛速度快、识别时间短等优点,具有良好的工程应用价值。
附图说明
图1是本发明的方法流程图。
具体实施方式
如图1所示,本发明公开了一种基于完全正交化算法的识别预应力桥梁现存预应力的方法,包括以下步骤:
(1)、在桥梁底面从左端依次向右布置有m个位移传感器;
21)、建立桥梁的振动微分方程:
将桥梁考虑为一单跨简支梁,假设固定荷载P(t)位于单跨简支梁上,其距离桥梁左端距离为xp,通过在桥梁底面粘贴的位移传感器可测得距桥梁左端x位置处t时刻的动态位移为y(x,t),则其振动微分方程:
其中单跨简支梁的跨长为L,密度为ρ,截面面积为A,粘性比例阻尼为c;桥梁抗弯刚度为EI,E是材料的弹性模量,I=bh0 3/12是梁截面的惯性矩,其中b是桥梁横截面宽度,h0是桥梁横截面高度;δ(x-xp)是狄拉克函数;预应力钢筋为后张无粘结,不考虑预应力钢筋偏心影响。
22)、得到简支梁的模态振型函数的矩阵形式:
基于模态叠加原理,假设梁的第i阶模态振型函数为Yi,梁的第i阶振幅为qi(t),则动态位移y(x,t)的模态形式表示为:
把方程(2)带入方程(1)并对方程两端的每一项均乘以Yi(x),考虑模态正交条件并对方程中x值从0到L进行积分,得到新方程:
简支预应力梁的模态振型函数标准化形式表示为:
在时域内通过Newmark积分将方程(3)用矩阵形式表示为:
其中,[I]是单位矩阵,{Q(t)}={q1(t),q2(t),…,qn(t)}T,[C]=diag(2ξiωi),{F(t)}={f1(t),f2(t),…,fn(t)}T,和/>分别为Q(t)的一阶导数和二阶导数。
假定在简支桥梁底面从左端支座依次向右布置有m个位移传感器,则在模态坐标下距简支梁左端支座xm处的第m个位移传感器实测竖向位移y(xm,t)表示为:
其中,Yi为梁的第i阶模态振型函数,qi(t)为梁的第i阶振幅,方程(6)用矩阵形式表示为:
则广义坐标下基于最小二乘法的振幅伪逆解表示为:
布置在桥梁底面的第j个位移传感器在t时刻测得的竖向位移向量y(xj,t)由广义正交多项式拟合为:
其中,Nf为正交多项式函数的阶数,Gi为第i阶正交多项式,ai为第i阶正交多项式的系数;
将方程(9)带入方程(7)并写成矩阵形式得:
式中,[A],[G]分别为正交多项式的系数矩阵和正交多项式矩阵;正交多项式的系数矩阵[A]可通过方程(10)由最小二乘法得到:
(4)、建立识别桥梁现存预应力值的系统方程T=({S}T{S})-1{S}T{b}。上述的方程(5)重写为:
式中,矩阵[K′]包含待识别的桥梁预应力值T,矩阵[C]包含桥梁模态阻尼ξi和模态频率ωi,矩阵[K]包含预应力桥梁的系统参数;
在时域内求解方程(12)即可得到桥梁现存预应力值T,方程(12)在时域内可表示为:
其中
则待识别预应力T由最小二乘法直接解得:
由方程(14)即可求得在t时刻预应力桥梁现存预应力值T(t),当取固定荷载P(t)在t1,t2,…,tw时刻的荷载值分别为P(t1),P(t2),…,P(tw),则方程(12)在时域内可重写为:
{S}N×wTw×1={b}N×1 (15)
其中
则待识别预应力T由最小二乘法直接解得:
T=({S}T{S})-1{S}T{b} (16)
5)、利用基于完全正交化算法求得桥梁现存预应力值T。
桥梁在固定荷载P(t)作用下,通过在桥梁底面安装位移传感器即可测得桥梁x1,x2,…xm处的竖向位移由方程(8)即可求得预应力桥梁的振幅{q}N×1,进而由方程(16)即可求得预应力桥梁在t1,t2,…,tw时刻的现存预应力值T(t1),T(t2),…,T(tw),实现由桥梁底面粘贴位移传感器识别预应力桥梁现存预应力值。在对方程(16)进行求解过程中,为避免对大型方阵求逆导致识别时间过长,识别精度过低,利用基于完全正交化算法实现高效率的识别桥梁现存预应力值。
对基于桥梁动力响应识别现存预应力值系统方程{S}T{S}T={S}T{b}进行求解,通过在Krylov子空间构造一组正交基实现完全正交化算法,最终得到计算效率更快的识别方法。
取基于桥梁动力响应识别现存预应力值系统矩阵为{S}T{S},由系统方程初始解T0得到的残差为:
r0={S}T{b}-{S}T{S}T0 (17)
取由r0构造得到的w维Krylov子空间为:
记β=‖r0‖2,则经过w步Arnoldi算法可构造一组标准正交基Vw=[v1,v2,…vw]和一个(w+1)×w的Hessenberg矩阵Hw+1,将Hessenberg矩阵Hw的最后一行元素删除即可得到大小为w×w的Hessenberg矩阵Hw,识别现存预应力值系统矩阵和标准正交基Vw以及Hessenberg矩阵Hw有下面的关系:
完全正交化算法采用正交投影理论,由标准正交基Vw和Hessenberg矩阵Hw即可在Krylov子空间解得预应力桥梁现存预应力值Tw,其表达式为:
其中,β=‖r0‖2,e1为第一个元素为1的w维单位向量,由公式(20)即可求得待识别的桥梁现存预应力值即为Tw。
Claims (2)
1.一种基于完全正交化算法识别预应力桥梁现存预应力的方法,其特征在于:包括以下步骤:
1)、在桥梁底面从左端依次向右布置有Nm个位移传感器;
3)、由测得的桥梁竖向位移计算预应力桥梁的振幅{q}N×1;N是测量位移向量中包含的振型数量;
4)、建立识别桥梁现存预应力值的系统方程;
5)、利用完全正交化算法,由系统方程初始解得到残差,由残差构造得到w维Krylov子空间,求得桥梁现存预应力值T;
所述步骤2)中具体包括以下步骤:
21)、建立桥梁的振动微分方程:
将桥梁考虑为一单跨简支梁,假设固定荷载P(t)位于单跨简支梁上,其距离桥梁左端距离为xp,通过在桥梁底面粘贴的位移传感器可测得距桥梁左端x位置处t时刻的动态位移为y(x,t),则其振动微分方程:
其中单跨简支梁的跨长为L,密度为ρ,截面面积为A,粘性比例阻尼为c;桥梁抗弯刚度为EI,E是材料的弹性模量,I=bh0 3/12是梁截面的惯性矩,其中b是桥梁横截面宽度,h0是桥梁横截面高度;δ(x-xp)是狄拉克函数;
22)、得到简支梁的模态振型函数的矩阵形式:
基于模态叠加原理,假设梁的第i阶模态振型函数为Yi,梁的第i阶振幅为qi(t),则动态位移y(x,t)的模态形式表示为:
把方程(2)带入方程(1)并对方程两端的每一项均乘以Yi(x),考虑模态正交条件并对方程中x值从0到L进行积分,得到新方程:
简支预应力梁的模态振型函数标准化形式表示为:
在时域内通过Newmark积分将方程(3)用矩阵形式表示为:
其中,[I]是单位矩阵,{Q(t)}={q1(t),q2(t),…,qn(t)}T,[C]=diag(2ξiωi),{F(t)}={f1(t),f2(t),…,fn(t)}T,和/>分别为Q(t)的一阶导数和二阶导数;ωi为模态频率;
假定在简支桥梁底面从左端支座依次向右布置有Nm个位移传感器,则在模态坐标下距简支梁左端支座xm处的第m个位移传感器实测竖向位移y(xm,t)表示为:
其中,Yi为梁的第i阶模态振型函数,qi(t)为梁的第i阶振幅,方程(6)用矩阵形式表示为:
则广义坐标下基于最小二乘法的振幅伪逆解表示为:
布置在桥梁底面的第j个位移传感器在t时刻测得的竖向位移向量y(xj,t)由广义正交多项式拟合为:
其中,Nf为正交多项式函数的阶数,Gi为第i阶正交多项式,ai为第i阶正交多项式的系数;
将方程(9)带入方程(7)并写成矩阵形式得:
式中,[A],[G]分别为正交多项式的系数矩阵和正交多项式矩阵;正交多项式的系数矩阵[A]可通过方程(10)由最小二乘法得到:
所述的步骤4)具体为:方程(5)重写为:
式中,矩阵[K′]包含待识别的桥梁预应力值T,矩阵[C]包含桥梁模态阻尼比ξi和模态频率ωi,矩阵[K]包含预应力桥梁的系统参数;
在时域内求解方程(12)即可得到桥梁现存预应力值T,方程(12)在时域内可表示为:
其中
则待识别预应力T由最小二乘法直接解得:
由方程(14)即可求得在t时刻预应力桥梁现存预应力值T(t),当取固定荷载P(t)在t1,t2,…,tw时刻的荷载值分别为P(t1),P(t2),…,P(tw),则方程(12)在时域内可重写为:
{S}N×wTw×1={b}N×1 (15)其中
则待识别预应力T由最小二乘法直接解得:
T=({S}T{S})-1{S}T{b} (16)
2.如权利要求1所述的基于完全正交化算法识别预应力桥梁现存预应力的方法,其特征在于:所述的步骤5)具体包括以下步骤:
取基于桥梁动力响应识别现存预应力值系统矩阵为{S}T{S},由系统方程初始解T0得到的残差为:
r0={S}T{b}-{S}T{S}T0 (17)
取由r0构造得到的w维Krylov子空间为:
记β=‖r0‖2,则经过w步Arnoldi算法可构造一组标准正交基Vw=[v1,v2,…vw]和一个(w+1)×w的Hessenberg矩阵Hw+1,将Hessenberg矩阵Hw+1的最后一行元素删除即可得到大小为w×w的Hessenberg矩阵Hw,识别现存预应力值系统矩阵和标准正交基Vw以及Hessenberg矩阵Hw有下面的关系:
完全正交化算法采用正交投影理论,由标准正交基Vw和Hessenberg矩阵Hw即可在Krylov子空间解得预应力桥梁现存预应力值Tw,其表达式为:
其中,β=‖r0‖2,e1为第一个元素为1的w维单位向量,由公式(20)即可求得待识别的桥梁现存预应力值即为Tw。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810879043.XA CN109063338B (zh) | 2018-08-03 | 2018-08-03 | 基于完全正交化算法识别预应力桥梁现存预应力的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810879043.XA CN109063338B (zh) | 2018-08-03 | 2018-08-03 | 基于完全正交化算法识别预应力桥梁现存预应力的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109063338A CN109063338A (zh) | 2018-12-21 |
CN109063338B true CN109063338B (zh) | 2023-06-09 |
Family
ID=64832896
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810879043.XA Active CN109063338B (zh) | 2018-08-03 | 2018-08-03 | 基于完全正交化算法识别预应力桥梁现存预应力的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109063338B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109948180B (zh) * | 2019-01-25 | 2021-05-25 | 北京航空航天大学 | 一种正交各向异性对边简支矩形薄板振动分析方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103557973A (zh) * | 2013-11-20 | 2014-02-05 | 重庆交通大学 | 在役结构预应力原位检测系统及方法 |
CN105865522A (zh) * | 2016-04-05 | 2016-08-17 | 江苏道亿智能科技有限公司 | 桥梁结构监测系统 |
CN106840336A (zh) * | 2017-01-25 | 2017-06-13 | 华北水利水电大学 | 基于qmr算法的桥面多轴移动荷载的识别方法 |
-
2018
- 2018-08-03 CN CN201810879043.XA patent/CN109063338B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103557973A (zh) * | 2013-11-20 | 2014-02-05 | 重庆交通大学 | 在役结构预应力原位检测系统及方法 |
CN105865522A (zh) * | 2016-04-05 | 2016-08-17 | 江苏道亿智能科技有限公司 | 桥梁结构监测系统 |
CN106840336A (zh) * | 2017-01-25 | 2017-06-13 | 华北水利水电大学 | 基于qmr算法的桥面多轴移动荷载的识别方法 |
Non-Patent Citations (1)
Title |
---|
齐忠强等.车辆荷载作用下预应力梁桥预应力动力识别研究.《第21届全国结构工程学术会议论文集》.2013, * |
Also Published As
Publication number | Publication date |
---|---|
CN109063338A (zh) | 2018-12-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Deng et al. | Identification of parameters of vehicles moving on bridges | |
Zheng et al. | Development of bridge influence line identification methods based on direct measurement data: A comprehensive review and comparison | |
CN104457681B (zh) | 一种基于应变模态的梁结构动挠度监测方法 | |
Deng et al. | Identification of dynamic vehicular axle loads: Demonstration by a field study | |
CN108226399A (zh) | 一种基于模态参数的张弦结构损伤组合识别方法 | |
CN104598753B (zh) | 一种基于Brakhage v方法的桥梁移动车辆荷载识别方法 | |
CN112989491B (zh) | 基于载荷应变线性叠加的应变场重构可视化方法及系统 | |
CN105930571A (zh) | 基于单位温度响应监测值的大跨钢桥有限元模型修正方法 | |
CN112949131B (zh) | 连续桥梁集群损伤诊断的概率损伤定位向量法 | |
CN110929382A (zh) | 基于弹性约束梁应变影响线曲率的桥梁损伤检测方法 | |
CN108280294B (zh) | 一种基于模态参数的索拱结构损伤组合识别方法 | |
CN104615888A (zh) | 一种基于广义最小残差方法的桥梁移动车辆荷载识别方法 | |
CN104091033A (zh) | 基于超单元结合虚拟变形法的桥梁静力有限元模型修正方法 | |
CN104122109B (zh) | 一种识别桥梁结构刚度损伤的系统 | |
CN109063338B (zh) | 基于完全正交化算法识别预应力桥梁现存预应力的方法 | |
Abramowicz et al. | Modelling and parameter identification of steel–concrete composite beams in 3D rigid finite element method | |
CN109033657A (zh) | 基于截断奇异值分解法识别预应力桥梁现存预应力的方法 | |
Li et al. | Identification of structural damping in time domain | |
CN109033656A (zh) | 基于两点步长梯度算法识别预应力桥梁现存预应力的方法 | |
CN108871221A (zh) | 一种桥梁监测装置 | |
Chen et al. | Condition assessment on thermal effects of a suspension bridge based on SHM oriented model and data | |
CN108920878A (zh) | 基于多步记忆下降法识别预应力桥梁现存预应力的方法 | |
CN110516297A (zh) | 基于长标距时频域信息的多尺度有限元模型修正方法 | |
CN114216634B (zh) | 一种浮置板轨道减振性能在线监测评估方法 | |
US20220050008A1 (en) | Method for calculating temperature-dependent mid-span vertical displacement of girder bridge |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |