CN109033656A - 基于两点步长梯度算法识别预应力桥梁现存预应力的方法 - Google Patents
基于两点步长梯度算法识别预应力桥梁现存预应力的方法 Download PDFInfo
- Publication number
- CN109033656A CN109033656A CN201810877998.1A CN201810877998A CN109033656A CN 109033656 A CN109033656 A CN 109033656A CN 201810877998 A CN201810877998 A CN 201810877998A CN 109033656 A CN109033656 A CN 109033656A
- Authority
- CN
- China
- Prior art keywords
- bridge
- equation
- matrix
- existing
- follows
- 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.)
- Withdrawn
Links
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/06—Power analysis or power optimisation
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
本发明公开了一种基于两点步长梯度算法的识别预应力桥梁现存预应力的方法,包括以下步骤:1)、在桥梁底面从左端依次向右布置有m个位移传感器;2)、测量在固定荷载P(t)作用下位移传感器测得的桥梁竖向位移3)、由测得的桥梁竖向位移计算预应力桥梁的振幅{q}N×1;4)、建立识别桥梁现存预应力值的系统方程T=({S}T{S})‑1{S}T{b};5)、利用基于两点步长梯度算法求得桥梁现存预应力值T。本发明通过建立预应力桥梁现存预应力与桥梁竖向动态位移响应之间的力学映射关系,实现由桥梁位移响应识别桥梁现存预应力值,具有识别速度快且精度较高等优点,具有良好的工程应用价值。
Description
技术领域
本发明属于预应力桥梁现存预应力的识别技术领域,尤其涉及一种基于两点步长梯度算法的识别预应力桥梁现存预应力的方法。
背景技术
近年来我国为控制车辆超载做出很多努力,但部分运输人员为了增加效益,安全意识浅薄,屡屡铤而走险,屡禁不绝的超限超重车辆对桥梁的安全性和耐久性产生不可估量的危害,这也是我国桥梁安全事故频发的主要原因。相比于发达国家的先进无损检测技术,我国桥梁结构检测技术仍处于相对落后的水平。依据我国桥梁“建养并重”的发展战略,针对目前我国已建预应力混凝土桥梁大多没有布设有效预应力值检测传感器的现状,通过测量桥梁动力响应识别预应力桥梁现存预应力值是一种可行的思路,可实现无损测量桥梁现存预应力值,进而评估桥梁结构状况和安全性能。
发明内容
本发明的目的是提供一种基于两点步长梯度算法由桥梁位移实时识别预应力桥梁现存预应力的方法,识别速度快且精度较高。
为达到上述目的,本发明采用的技术方案是:一种基于两点步长梯度算法的识别预应力桥梁现存预应力的方法,包括以下步骤:
1)、在桥梁底面从左端依次向右布置有m个位移传感器;
2)、测量在固定荷载P(t)作用下位移传感器测得的桥梁竖向位移
3)、由测得的桥梁竖向位移计算预应力桥梁的振幅{q}N×1;
4)、建立识别桥梁现存预应力值的系统方程T=({S}T{S})-1{S}T{b};
5)、利用基于两点步长梯度算法,得到系统方程初始解,得到现存预应力值的迭代解,求得桥梁现存预应力值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进行积分,得到新方程:
其中,ξi和分别是第i阶模态的阻尼比和模态质量,fi(t)=P(t)Yi(xp)是位于单跨简支梁上的模态荷载;为qi(t)的二阶导数,为qi(t)的一阶导数;
简支预应力梁的模态振型函数标准化形式表示为:
在时域内通过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)的一阶导数和二阶导数;
23)、得到桥梁竖向位移向量
假定在简支桥梁底面从左端支座依次向右布置有m个位移传感器,则在模态坐标下距简支梁左端支座xm处的第m个位移传感器实测竖向位移y(xm,t)表示为:
其中,Yi为梁的第i阶模态振型函数,qi(t)为梁的第i阶振幅,方程(6)用矩阵形式表示为:
式中,是第Nm个位移传感器测得的位移向量,N是测量位移向量中包含的振型数量;
则广义坐标下基于最小二乘法的振幅伪逆解表示为:
式中,为矩阵的逆,[Y]T为矩阵[Y]的转置;
布置在桥梁底面的第j个位移传感器在t时刻测得的竖向位移向量y(xj,t)由广义正交多项式拟合为:
其中,Nf为正交多项式函数的阶数,Gi为第i阶正交多项式,ai为第i阶正交多项式的系数;
将方程(9)带入方程(7)并写成矩阵形式得:
式中,[A],[G]分别为正交多项式的系数矩阵和正交多项式矩阵;正交多项式的系数矩阵[A]可通过方程(10)由最小二乘法得到:
将方程(11)带入方程(10)即可得到竖向位移向量
所述的步骤3)具体为:将竖向位移向量代入方程(8)即可求得预应力桥梁的振幅{q}N×1。
所述的步骤4)具体为:
所述的方程(5)重写为:
式中,矩阵[K′]包含待识别的桥梁预应力值T,矩阵[C]包含桥梁模态阻尼ξi和模态频率ωi,矩阵[K]包含预应力桥梁的系统参数;
在时域内求解方程(12)即可得到桥梁现存预应力值T,方程(12)在时域内可表示为:
其中
则待识别预应力T由最小二乘法直接解得:
式中,矩阵{B}T为矩阵{B}的转置,即为包含了方程(12)所有右端项信息的向量;
由方程(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)
式中,矩阵{S}T为矩阵{S}的转置,向量{b}N×1与类似,即为包含了方程(12)所有右端项信息的向量,可由桥梁动态参数和测量位移响应推求。
所述的步骤5)具体包括以下步骤:
取系统方程初始解T0,则待识别的现存预应力值第k+1步迭代解为:
Tk+1=Tk-Dkgk (17)
式中,gk为第k+1步迭代梯度可根据迭代收敛要求确定,两点步长梯度算法搜索方向向量D的第k步迭代解为:
Dk=αkI (18)
式中,I为对角元素为1的w行w列单位矩阵,αk为第k步迭代步长,
迭代向量sk第k步迭代解为:
sk=Tk+1-Tk (20)
迭代向量lk第k步迭代解为
lk=gk+1-gk (21)
将公式(18)-(21)代入公式(17),当Tk+1-Tk≤ε满足识别精度误差要求ε时,即可由两点步长梯度算法求得待识别的桥梁现存预应力值即为Tk+1,最终实现由桥梁动态位移响应实时识别现存预应力值。
本发明具有的优点是:本发明提出一种基于两点步长梯度算法由桥梁位移实时识别预应力桥梁现存预应力的方法,通过建立预应力桥梁现存预应力与桥梁竖向动态位移响应之间的力学映射关系,实现由桥梁位移响应识别桥梁现存预应力值。为提高预应力值的识别效率,对预应力识别系统矩阵采用两点步长梯度算法以节约识别时间,最终实现实时识别桥梁现存预应力值,具有识别速度快且精度较高等优点,具有良好的工程应用价值。
附图说明
图1是本发明的方法流程图。
具体实施方式
如图1所示,本发明公开了一种基于两点步长梯度算法的识别预应力桥梁现存预应力的方法,包括以下步骤:
(1)、在桥梁底面从左端依次向右布置有m个位移传感器;
(2)、测量在固定荷载P(t)作用下位移传感器测得的桥梁竖向位移
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进行积分,得到新方程:
其中,ξi和分别是第i阶模态的阻尼比和模态质量,fi(t)=P(t)Yi(xp)是位于单跨简支梁上的模态荷载;为qi(t)的二阶导数,为qi(t)的一阶导数;
简支预应力梁的模态振型函数标准化形式表示为:
在时域内通过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)的一阶导数和二阶导数。
23)、得到桥梁竖向位移向量
假定在简支桥梁底面从左端支座依次向右布置有m个位移传感器,则在模态坐标下距简支梁左端支座xm处的第m个位移传感器实测竖向位移y(xm,t)表示为:
其中,Yi为梁的第i阶模态振型函数,qi(t)为梁的第i阶振幅,方程(6)用矩阵形式表示为:
式中,是第Nm个位移传感器测得的位移向量,N是测量位移向量中包含的振型数量;
则广义坐标下基于最小二乘法的振幅伪逆解表示为:
式中,为矩阵的逆,[Y]T为矩阵[Y]的转置;
布置在桥梁底面的第j个位移传感器在t时刻测得的竖向位移向量y(xj,t)由广义正交多项式拟合为:
其中,Nf为正交多项式函数的阶数,Gi为第i阶正交多项式,ai为第i阶正交多项式的系数;
将方程(9)带入方程(7)并写成矩阵形式得:
式中,[A],[G]分别为正交多项式的系数矩阵和正交多项式矩阵;正交多项式的系数矩阵[A]可通过方程(10)由最小二乘法得到:
将方程(11)带入方程(10)即可得到竖向位移向量
(3)、将竖向位移向量代入方程(8)即可求得预应力桥梁的振幅{q}N×1。
(4)、建立识别桥梁现存预应力值的系统方程T=({S}T{S})-1{S}T{b}。上述的方程(5)重写为:
式中,矩阵[K′]包含待识别的桥梁预应力值T,矩阵[C]包含桥梁模态阻尼ξi和模态频率ωi,矩阵[K]包含预应力桥梁的系统参数;
在时域内求解方程(12)即可得到桥梁现存预应力值T,方程(12)在时域内可表示为:
其中
则待识别预应力T由最小二乘法直接解得:
式中,矩阵{B}T为矩阵{B}的转置,即为包含了方程(12)所有右端项信息的向量。
由方程(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)
式中,矩阵{S}T为矩阵{S}的转置,向量{b}N×1与类似,即为包含了方程(12)所有右端项信息的向量,可由桥梁动态参数和测量位移响应推求。
5)、利用基于两点步长梯度算法求得桥梁现存预应力值T。
桥梁在固定荷载P(t)作用下,通过在桥梁底面安装位移传感器即可测得桥梁x1,x2,…xm处的竖向位移由方程(8)即可求得预应力桥梁的振幅{q}N×1,进而由方程(16)即可求得预应力桥梁在t1,t2,…,tw时刻的现存预应力值T(t1),T(t2),…,T(tw),实现由桥梁底面粘贴位移传感器识别预应力桥梁现存预应力值。在对方程(16)进行求解过程中,为避免对大型方阵求解导致识别时间过长,特引入两点步长梯度算法提高现存预应力值的识别效率,实现高效率的识别桥梁现存预应力值。
两点步长梯度算法可将系统方程求解问题转化为迭代算法,取系统方程初始解T0,则待识别的现存预应力值第k+1步迭代解为:
Tk+1=Tk-Dkgk (17)
式中,gk为第k+1步迭代梯度可根据迭代收敛要求确定,两点步长梯度算法搜索方向向量D的第k步迭代解为:
Dk=αkI (18)
式中,I为对角元素为1的w行w列单位矩阵,αk为第k步迭代步长,
迭代向量sk第k步迭代解为:
sk=Tk+1-Tk (20)
迭代向量lk第k步迭代解为
lk=gk+1-gk (21)
将公式(18)-(21)代入公式(17),当Tk+1-Tk≤ε满足识别精度误差要求ε时,ε一般取一个小值,根据精度要求自行设置即可,比如取ε为0.01Tk,就是要求两次计算相对误差小于百分之一即可,即可由两点步长梯度算法求得待识别的桥梁现存预应力值即为Tk+1,最终实现由桥梁动态位移响应实时识别现存预应力值。
Claims (5)
1.一种基于两点步长梯度算法识别预应力桥梁现存预应力的方法,其特征在于:包括以下步骤:
1)、在桥梁底面从左端依次向右布置有m个位移传感器;
2)、测量在固定荷载P(t)作用下位移传感器测得的桥梁竖向位移
3)、由测得的桥梁竖向位移计算预应力桥梁的振幅{q}N×1;
4)、建立识别桥梁现存预应力值的系统方程T=({S}T{S})-1{S}T{b};
5)、利用基于两点步长梯度算法,得到系统方程初始解,得到现存预应力值的迭代解,求得桥梁现存预应力值T。
2.如权利要求1所述的基于两点步长梯度算法识别预应力桥梁现存预应力的方法,其特征在于:所述步骤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进行积分,得到新方程:
其中,ξi和分别是第i阶模态的阻尼比和模态质量,fi(t)=P(t)Yi(xp)是位于单跨简支梁上的模态荷载;为qi(t)的二阶导数,为qi(t)的一阶导数;
简支预应力梁的模态振型函数标准化形式表示为:
在时域内通过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)的一阶导数和二阶导数;
23)、得到桥梁竖向位移向量
假定在简支桥梁底面从左端支座依次向右布置有m个位移传感器,则在模态坐标下距简支梁左端支座xm处的第m个位移传感器实测竖向位移y(xm,t)表示为:
其中,Yi为梁的第i阶模态振型函数,qi(t)为梁的第i阶振幅,方程(6)用矩阵形式表示为:
式中,是第Nm个位移传感器测得的位移向量,N是测量位移向量中包含的振型数量;
则广义坐标下基于最小二乘法的振幅伪逆解表示为:
式中,为矩阵的逆,[Y]T为矩阵[Y]的转置;
布置在桥梁底面的第j个位移传感器在t时刻测得的竖向位移向量y(xj,t)由广义正交多项式拟合为:
其中,Nf为正交多项式函数的阶数,Gi为第i阶正交多项式,ai为第i阶正交多项式的系数;
将方程(9)带入方程(7)并写成矩阵形式得:
式中,[A],[G]分别为正交多项式的系数矩阵和正交多项式矩阵;正交多项式的系数矩阵[A]可通过方程(10)由最小二乘法得到:
将方程(11)带入方程(10)即可得到竖向位移向量
3.如权利要求2所述的基于两点步长梯度算法识别预应力桥梁现存预应力的方法,其特征在于:所述的步骤3)具体为:
将竖向位移向量代入方程(8)即可求得预应力桥梁的振幅{q}N×1。
4.如权利要求3所述的基于两点步长梯度算法识别预应力桥梁现存预应力的方法,其特征在于:所述的步骤4)具体为:
所述的方程(5)重写为:
式中,矩阵[K′]包含待识别的桥梁预应力值T,矩阵[C]包含桥梁模态阻尼ξi和模态频率ωi,矩阵[K]包含预应力桥梁的系统参数;
在时域内求解方程(12)即可得到桥梁现存预应力值T,方程(12)在时域内可表示为:
其中
则待识别预应力T由最小二乘法直接解得:
式中,矩阵{B}T为矩阵{B}的转置,即为包含了方程(12)所有右端项信息的向量;
由方程(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)
式中,矩阵{S}T为矩阵{S}的转置,向量{b}N×1与类似,即为包含了方程(12)所有右端项信息的向量,可由桥梁动态参数和测量位移响应推求。
5.如权利要求4所述的基于两点步长梯度算法识别预应力桥梁现存预应力的方法,其特征在于:所述的步骤5)具体包括以下步骤:
取系统方程初始解T0,则待识别的现存预应力值第k+1步迭代解为:
Tk+1=Tk-Dkgk (17)
式中,gk为第k+1步迭代梯度可根据迭代收敛要求确定,两点步长梯度算法搜索方向向量D的第k步迭代解为:
Dk=αkI (18)
式中,I为对角元素为1的w行w列单位矩阵,αk为第k步迭代步长,
迭代向量sk第k步迭代解为:
sk=Tk+1-Tk (20)
迭代向量lk第k步迭代解为
lk=gk+1-gk (21)
将公式(18)-(21)代入公式(17),当Tk+1-Tk≤ε满足识别精度误差要求ε时,即可由两点步长梯度算法求得待识别的桥梁现存预应力值即为Tk+1,最终实现由桥梁动态位移响应实时识别现存预应力值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810877998.1A CN109033656A (zh) | 2018-08-03 | 2018-08-03 | 基于两点步长梯度算法识别预应力桥梁现存预应力的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810877998.1A CN109033656A (zh) | 2018-08-03 | 2018-08-03 | 基于两点步长梯度算法识别预应力桥梁现存预应力的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109033656A true CN109033656A (zh) | 2018-12-18 |
Family
ID=64649364
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810877998.1A Withdrawn CN109033656A (zh) | 2018-08-03 | 2018-08-03 | 基于两点步长梯度算法识别预应力桥梁现存预应力的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109033656A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111783198A (zh) * | 2020-06-21 | 2020-10-16 | 西北工业大学 | 一种基于双梁模型的悬索桥精细化动力分析方法 |
CN113089495A (zh) * | 2021-04-30 | 2021-07-09 | 太原科技大学 | 预应力混凝土简支桥梁锚下有效预应力检测方法 |
-
2018
- 2018-08-03 CN CN201810877998.1A patent/CN109033656A/zh not_active Withdrawn
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111783198A (zh) * | 2020-06-21 | 2020-10-16 | 西北工业大学 | 一种基于双梁模型的悬索桥精细化动力分析方法 |
CN113089495A (zh) * | 2021-04-30 | 2021-07-09 | 太原科技大学 | 预应力混凝土简支桥梁锚下有效预应力检测方法 |
CN113089495B (zh) * | 2021-04-30 | 2022-06-28 | 太原科技大学 | 预应力混凝土简支桥梁锚下有效预应力检测方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zhang et al. | A new smart traffic monitoring method using embedded cement-based piezoelectric sensors | |
Zhu et al. | Damage detection in simply supported concrete bridge structure under moving vehicular loads | |
Rosales et al. | Crack detection in beam-like structures | |
CN110929382B (zh) | 基于弹性约束梁应变影响线曲率的桥梁损伤检测方法 | |
CN105930571A (zh) | 基于单位温度响应监测值的大跨钢桥有限元模型修正方法 | |
CN106777549B (zh) | 一种面向荷载试验的桥梁多尺度有限元模拟方法 | |
He et al. | Restoring force and dynamic loadings identification for a nonlinear chain-like structure with partially unknown excitations | |
CN109858156A (zh) | 基于车桥耦合振动的车辆与结构信息同时识别方法 | |
CN104408286B (zh) | 基于正交分解理论的平面构件变形分解与振型识别方法 | |
Lu et al. | Hybrid sensitivity matrix for damage identification in axially functionally graded beams | |
Farag | Application of smart structural system for smart sustainable cities | |
CN109033656A (zh) | 基于两点步长梯度算法识别预应力桥梁现存预应力的方法 | |
CN104091033A (zh) | 基于超单元结合虚拟变形法的桥梁静力有限元模型修正方法 | |
CN109033657A (zh) | 基于截断奇异值分解法识别预应力桥梁现存预应力的方法 | |
CN105223272A (zh) | 一种快速测空心板梁桥横向分布系数的方法及质量激振器 | |
Peralta-Braz et al. | Design optimisation of piezoelectric energy harvesters for bridge infrastructure | |
Li et al. | Identification of structural damping in time domain | |
JP2004264235A (ja) | 構造物の損傷部位検出装置、及び構造物の損傷部位検出方法 | |
CN109002638A (zh) | 基于共轭残差平方算法识别预应力桥梁现存预应力的方法 | |
Matsuda et al. | Fundamental study using new test loading scheme for steel frame subassembly with damper connection details | |
CN109063338B (zh) | 基于完全正交化算法识别预应力桥梁现存预应力的方法 | |
Choi et al. | Quasi-static responses estimation of a cable-stayed bridge from displacement data at a limited number of points | |
CN108920878A (zh) | 基于多步记忆下降法识别预应力桥梁现存预应力的方法 | |
CN109002637A (zh) | 基于加速超松弛迭代法识别预应力桥梁现存预应力的方法 | |
Björklund | Dynamic analysis of a railway bridge subjected to high speed trains |
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 | ||
WW01 | Invention patent application withdrawn after publication |
Application publication date: 20181218 |
|
WW01 | Invention patent application withdrawn after publication |