CN104298870B - 一种移动荷载下简支梁损伤和移动力同时识别方法 - Google Patents
一种移动荷载下简支梁损伤和移动力同时识别方法 Download PDFInfo
- Publication number
- CN104298870B CN104298870B CN201410527720.3A CN201410527720A CN104298870B CN 104298870 B CN104298870 B CN 104298870B CN 201410527720 A CN201410527720 A CN 201410527720A CN 104298870 B CN104298870 B CN 104298870B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- msup
- mfrac
- mtd
- 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
Landscapes
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本发明提出了一种移动荷载下简支梁损伤和移动力同时识别方法:测量结构在损伤状态的加速度动力响应,以结构的加速度响应作为结构损伤和移动力识别的动力指标;根据结构的数值模型,将未知移动荷载采用切比雪夫多项式表示,并利用梁单元形函数概念将移动荷载等效成单元节点力,分别计算结构加速度响应对损伤和移动力参数的灵敏度矩阵;结合敏感性分析结果,以结构测量加速度响应和计算所得的有限元加速度响应的差值作为目标函数,采用一阶泰勒展开的识别方程逐步迭代同时识别出简支梁的损伤和移动力。通过本发明方法,能够仅仅通过结构所测的少数加速度响应来实现结构损伤和移动力的同时识别,具有较强的实用性。
Description
技术领域
本发明属于结构健康监测技术领域,具体涉及一种移动荷载下简支梁损伤和移动力同时识别方法。
背景技术
结构的损伤识别和外部荷载识别是结构健康监测体系的两个重要方面,在结构安全维护和安全评定中具有重要作用。由于结构的振动特性(如时程响应、频响函数、固有频率、振型等)是结构物理参数(如刚度、质量、阻尼等)的函数,结构损伤即意味着结构参数的改变,而物理参数的改变必然引起结构振动特性的改变。目前基于结构振动信息的损伤识别方法,主要分为频域法和时域法。频域法是指利用实测结构的模态参数(固有频率和振型),通过求解结构动力特征值的反问题识别结构的物理参数的变化。时域法是指通过安装在结构上的传感设备,对结构的振动进行实时监测,可以获得结构不同位置的动力响应,对损伤前后结构的动力响应变化进行分析及处理,就有可能获得结构物理参数的变化情况。
动态载荷识别一般是在预先已知结构的参数信息的情况下,利用测量结构响应(位移、速度、加速度及频响函数等),识别出无法测量的未知输入力。频域法和时域法是动态荷载识别的两类经典方法。频域法利用频响函数及模态参数在频率域中进行识别;时域法采用模态坐标变换,将运动方程转化为非耦合的型式,并假定微小时间间隔内动态荷载为阶跃函数,对非耦合方程进行求解,从而利用响应历程获得动态荷载。频域识别法已形成了较完善的理论及计算方法,但是该方法要求采样信号有一定的数据长度、一般适应于平稳振动信号,且高阶模态截断会带来模型误差,必然会影响动荷载的识别结果。由于时域法可直接得到荷载时间历程,无需时频转化,不存在中间过程的误差,所得结果比较直观,便于工程实际应用,且时域荷载识别法是处理非平稳载荷的根本途径,因此荷载识别方法多数基于时域内的结构振动响应进行识别。
在基于结构振动信息的损伤识别过程中,一般需要测量或已知结构外部荷载的时间历程,然而在实际应用中,结构的未知荷载和未知损伤常常是共存的,特别是在桥梁结构中,车辆的移动荷载通常难以准确测量,因此在未知移动荷载下受损结构的损伤识别和移动力识别情况同时存在。
李杰和陈隽首次将这种在时域中识别结构参数同时反演未知激励的反问题归结为动力复合反演问题,并提出全量补偿法、分组归一化统计平均法等多种复合反演算法,较好地解决了部分输入未知条件下的结构参数识别及荷载反演问题[未知输入条件下的结构物理参数识别研究,计算力学学报,1999,16(1):32-40]。Lu等[Identification of systemparameters and Input force from output only,Mechanical Systems and SignalProcessing,2007,21:2099-2111]基于结构响应的灵敏度利用两步反演法同时识别系统参数和输入激励,未知荷载由一系列正弦函数加一个常数逼近组成。Zhang等[A probabilistic damage identification approach for structures with uncertainties underunknown input,Mechanical Systems and Signal Processing,2011,25:1126-1145]利用切比雪夫多项式表示未知荷载,把识别荷载等效为识别多项式系数,识别中能够对损伤参数和荷载参数二者同时进行迭代修正。Zhang等[Identification of coexistent loadand damage,Structural and Multidisciplinary Optimization,2010,41(2):243-253]基于虚拟变形方法(VDM),用虚拟变形等效结构的损伤,继而借助未损伤结构通过实测损伤结构响应识别荷载和损伤,包括损伤类型和大小,该方法要求传感器数目不小于未知荷载数目和虚拟变形数目的总和。以上方法都能解决未知荷载和结构损伤的同时识别,但仅局限于固定位置的动态荷载识别。
发明内容
针对实际桥梁结构中车辆的移动荷载难以准确测量,移动车辆荷载同时也会对桥梁结构产生损伤影响,因此桥梁结构移动荷载和损伤的同时识别问题尤为重要,本发明采用时域内动力响应灵敏度的模型修正方法,将未知移动荷载采用切比雪夫多项式表示,结合梁单元形函数概念将移动荷载等效成单元节点力,仅利用少数测量的动力响应同时识别出简支梁的损伤状况和移动荷载的时程曲线。
本发明提供一种用于同时识别简支梁损伤和移动力的方法,解决在未知移动荷载作用下简支梁的损伤状况以及移动力识别的快捷计算方法。
本发明提供了一种移动荷载下简支梁损伤和移动力同时识别方法,其步骤包括:
(1)测量结构损伤状态下的加速度响应YA,假定结构的初始损伤值和移动力的正交系数分别为α0和c0;
(2)在状态空间域内表示移动荷载下结构的运动方程,移动力采用切比雪夫多项式表示为在移动荷载表达式Fδ(l-vt)中,F表示以一恒定速度v运动的时变荷载,l表示结构的某一位置点,vt表示移动荷载在t时刻的位置,δ(·)表示狄拉克函数。结合梁单元形函数概念将任一时刻的移动力等效成单元节点力,通过离散化的运动方程计算在状态空间域内的马尔科夫系数矩阵HL,再计算状态空间域内的结构加速度响应YA:
(2.1)状态空间域内移动荷载下结构的运动方程表示为:
其中R=[Rd-RaM-1K Rv-RaM-1C],K*为系统矩阵,B*为输入矩阵,Y表示输出矩阵,Ra,Rv和Rd分别表示测量加速度,速度和位移的映射矩阵,ck和Tk分别表示移动荷载的正交系数和切比雪夫正交式。
(2.2)将状态空间域内连续运动方程离散化,计算在状态空间域内的马尔科夫系数矩阵HL,离散化的状态空间方程可表示为下式
其中A=exp(K*h),B=K*-1(A-I)B*.N表示采样点数量,h表示时间间隔。结合以上两式,输出矩阵可用下式表示:
令:
其中,HL表示状态空间域内的马尔科夫系数矩阵,HL中各元素Hm(m=0,1…N-1)表示在单位脉冲荷载下离散结构的响应,Hm的计算式如下:
(2.3)根据步骤(2.2)中计算结果HL,结构的输出矩阵可简写为:
其中YA表示在状态空间域内结构的加速度响应。
(3)利用步骤(2)中计算所得状态空间域内结构的加速度响应,进行结构加速度响应对单元刚度损伤参数和移动力正交参数的一阶灵敏度矩阵和
其中,和分别表示结构加速度响应对单元刚度损伤参数和移动力正交参数的一阶灵敏度矩阵,表示结构马尔科夫系数关于单元刚度损伤参数的一阶偏导,具体计算如下:
其中和的计算如下:
(4)将步骤(1)中的测量加速度响应YE和步骤(2)中计算的结构加速度响应值YA的差值作为目标函数,表示如下:
ΔY=YE-YA
(5)根据步骤(3)计算所得的动力加速度响应灵敏度矩阵和采用基于动力响应灵敏度模型修正的识别方程计算结构的损伤α和移动力的正交系数c。
(6)重复步骤(2)-(5),进行下一次循环的计算,可以得到修正后的结构损伤值和移动力的正交系数,直至目标函数达到设定允许值停止计算。
与现有技术相比,本发明具有以下有益效果:传统的结构损伤识别理论都是建立在结构输入和输出均已知的基础上的,但是对于实际桥梁结构而言,其移动车载输入信息往往很难测得或者很难被准确地测试,本发明仅利用少数几个测量的动力响应可同时识别出结构的损伤和移动外荷载的时程曲线,实现未知移动荷载下结构损伤和力的同时识别。
附图说明
图1为移动荷载下简支梁的有限元模型示意图;
图2表示测量加速度AY(6)和AY(10)示意图;
图3表示测量加速度对第六个单元刚度参数α6的一阶导数;其中:
图3(a)表示测量加速度AY(6)对第六个单元刚度参数α6的一阶导数;
图3(b)表示测量加速度AY(10)对第六个单元刚度参数α6的一阶导数;
图4表示测量加速度对移动力第一个正交参数c1的一阶导数;其中:
图4(a)表示测量加速度AY(6)对移动力第一个正交参数c1的一阶导数;
图4(b)表示测量加速度AY(10)对移动力第一个正交参数c1的一阶导数;
图5表示测量加速度对移动力第二个正交参数c2的一阶导数;其中:
图5(a)表示测量加速度AY(6)对移动力第二个正交参数c2的一阶导数;
图5(b)表示测量加速度AY(10)对移动力第二个正交参数c2的一阶导数;
图6表示简支梁的损伤识别结果图;
图7表示简支梁未知移动力的识别结果;其中:
图7(a)表示无测量噪音情况下移动力的识别结果;
图7(b)表示5%测量噪音情况下移动力的识别结果。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
以图1中移动荷载下的简支梁为研究对象,来阐述结构损伤和移动力荷载同时识别的过程。简支梁全长12m,弹性模量210GPa,密度为7800kg/m3,泊松比为0.3。梁横截面为0.5m×0.5m。瑞尼阻尼系数a和b分别为0.8180和7.823×10-5,假设未知移动荷载F=-12000(sin(24t)+0.5sin(12t))N以v=12m/s的速度从梁左端行驶到右端。假设简支梁6m处刚度出现15%的损伤,采样频率为120Hz,测量加速度为AY(6)和AY(10),其中AY(6)表示测量节点6的Y方向的加速度,移动力的切比雪夫正交多项式的阶数为15。
图2表示测量加速度AY(6)和AY(10)示意图,通过步骤(3)可计算测量加速度响应对损伤和移动力参数的一阶灵敏度矩阵,图3(a)和(b)分别表示测量加速度AY(6)和AY(10)对第六个单元刚度参数α6的一阶导数,图4(a)和(b)分别表示测量加速度AY(6)和AY(10)对移动力第一个正交参数c1的一阶导数,图5(a)和(b)分别表示测量加速度AY(6)和AY(10)对移动力第两个正交参数c2的一阶导数。为了验证所计算的动力响应灵敏度矩阵正确性,采用Newmark法(Lu and Law,Features of dynamic response sensitivity and itsapplication in damage detection,Journal of Sound and Vibration,2007,303:305-329)进行对比参考。由图3-5可以看出,采用提出方法计算的加速度响应对损伤和移动力参数的一阶灵敏度矩阵与Newmark法计算结果保持一致,验证了其正确性。
根据计算的测量加速度响应对损伤和移动力参数的一阶灵敏度矩阵,通过步骤(4)-(6)同时识别出简支梁的损伤状况和移动力。图6表示简支梁在无噪音和5%噪音下损伤结果示意图,由图6可以看出,在无噪音情况下,简支梁的损伤位置和损伤程度都能准确识别出来,相对误差值为1.64%;在5%噪音情况下,简支梁的损伤位置和损伤程度也能识别出来,但其相对误差值增大,在简支梁6m位置时的刚度损伤为11.84%。当考虑测量噪音时,损伤结果的相对误差会增大,说明损伤结果对测量噪音较为敏感。图7所示为简支梁的移动力时程曲线识别结果,其中图7(a)和(b)分别表示简支梁在无测量噪声和5%测量噪声情况下移动力时程曲线识别结果。从图7可以看出,在无噪音和5%测量噪音情况下,识别出的移动力时程曲线与实际移动荷载都很接近,其相对误差值分别为8%和10.97%。误差来源主要来自切比雪夫多项式的阶数,采样时间以及采样频率等。
由图6和图7可以看出,采用本发明,可以在未知移动力状况下同时准确识别简支梁的损伤状况和移动力的时程曲线。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (5)
1.一种移动荷载下简支梁损伤和移动力同时识别方法,其特征在于,所述方法包括如下步骤:
(1)测量结构损伤状态下的加速度响应YE,假定结构的初始损伤值为α0和初始移动力正交系数为c0;
(2)在状态空间域内表示移动荷载下结构的运动方程,移动力采用切比雪夫多项式表示为在移动荷载表达式Fδ(l-vt)中,F表示以一恒定速度v运动的时变荷载,ck和Tk分别表示移动力正交系数和切比雪夫正交式,l表示结构的某一位置点,vt表示移动荷载在t时刻的位置,δ(·)表示狄拉克函数,结合梁单元形函数概念将任一时刻的移动力等效成单元节点力,通过离散化的运动方程计算在状态空间域内的马尔科夫系数矩阵HL,再计算状态空间域内的结构加速度响应YA;
(3)利用步骤(2)中计算所得状态空间域内结构的加速度响应YA,计算结构加速度响应对单元刚度损伤参数和移动力正交参数的一阶灵敏度矩阵
(4)将步骤(1)中的测量加速度响应YE和步骤(2)中计算的结构加速度响应值YA的差值作为目标函数,表示如下:
ΔY=YE-YA
(5)根据步骤(3)计算所得的结构加速度响应对单元损伤参数和移动力正交参数的一阶灵敏度矩阵采用基于动力响应灵敏度模型修正的识别方程计算结构的损伤α和移动力的正交系数c;
(6)重复步骤(2)-(5),进行下一次循环的计算,可以得到修正后的结构损伤值和移动力的正交系数,直至目标函数达到设定允许值停止计算。
2.如权利要求1所述的移动荷载下简支梁损伤和移动力同时识别方法,其特征在于,所述步骤(2)中状态空间域内移动荷载下结构的运动方程表示为:
<mrow>
<mover>
<mi>X</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<msup>
<mi>K</mi>
<mo>*</mo>
</msup>
<mi>X</mi>
<mo>+</mo>
<msup>
<mi>B</mi>
<mo>*</mo>
</msup>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>Q</mi>
</munderover>
<msub>
<mi>c</mi>
<mi>k</mi>
</msub>
<msub>
<mi>T</mi>
<mi>k</mi>
</msub>
<mi>&delta;</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>-</mo>
<mi>v</mi>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>Y</mi>
<mo>=</mo>
<mi>R</mi>
<mi>X</mi>
<mo>+</mo>
<mi>D</mi>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>Q</mi>
</munderover>
<msub>
<mi>c</mi>
<mi>k</mi>
</msub>
<msub>
<mi>T</mi>
<mi>k</mi>
</msub>
<mi>&delta;</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>-</mo>
<mi>v</mi>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
其中R=[Rd-RaM-1K Rv-RaM-1C],K*为系统矩阵,B*为输入矩阵,Y表示输出矩阵,Ra,Rv和Rd分别表示测量加速度,速度和位移的映射矩阵。
3.如权利要求2所述的移动荷载下简支梁损伤和移动力同时识别方法,其特征在于,所述步骤(2)中通过离散化的运动方程计算在状态空间域内的马尔科夫系数矩阵HL具体为:
<mrow>
<mi>X</mi>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>A</mi>
<mi>X</mi>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>B</mi>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>Q</mi>
</munderover>
<msub>
<mi>c</mi>
<mi>k</mi>
</msub>
<msub>
<mi>T</mi>
<mi>k</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mi>&delta;</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>-</mo>
<mi>v</mi>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>Y</mi>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>R</mi>
<mi>X</mi>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>D</mi>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>Q</mi>
</munderover>
<msub>
<mi>c</mi>
<mi>k</mi>
</msub>
<msub>
<mi>T</mi>
<mi>k</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mi>&delta;</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>-</mo>
<mi>v</mi>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mi>N</mi>
</mrow>
其中A=exp(K*h),B=K*-1(A-I)B*,N表示采样点数量,h表示时间间隔,结合以上两式,输出矩阵用下式表示:
令:
其中,HL表示状态空间域内的马尔科夫系数矩阵,HL中各元素Hm表示在单位脉冲荷载下离散结构的响应,m=0,1…N-1,Hm的计算式如下:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>H</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mi>D</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>H</mi>
<mi>m</mi>
</msub>
<mo>=</mo>
<msup>
<mi>RA</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msup>
<mi>B</mi>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mo>...</mo>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>.</mo>
</mrow>
4.如权利要求3所述的移动荷载下简支梁损伤和移动力同时识别方法,其特征在于,所述计算状态空间域内的结构加速度响应YA具体为:
根据步骤(2)中计算结果HL,结构的输出矩阵简写为:
<mrow>
<msup>
<mi>Y</mi>
<mi>A</mi>
</msup>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>Q</mi>
</munderover>
<msub>
<mi>c</mi>
<mi>k</mi>
</msub>
<msub>
<mi>H</mi>
<mi>L</mi>
</msub>
<msub>
<mi>T</mi>
<mi>k</mi>
</msub>
<mi>&delta;</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>-</mo>
<mi>v</mi>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
其中YA表示在状态空间域内结构的加速度响应。
5.如权利要求1至4任一项所述的移动荷载下简支梁损伤和移动力同时识别方法,其特征在于,所述步骤(3)具体根据下式进行:
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>Y</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mrow>
<mo>(</mo>
<mrow>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>Q</mi>
</munderover>
<msub>
<mi>c</mi>
<mi>k</mi>
</msub>
<msub>
<mi>H</mi>
<mi>L</mi>
</msub>
<msub>
<mi>T</mi>
<mi>k</mi>
</msub>
<mi>&delta;</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>l</mi>
<mo>-</mo>
<mi>v</mi>
<mi>t</mi>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>H</mi>
<mi>L</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>Q</mi>
</munderover>
<msub>
<mi>c</mi>
<mi>k</mi>
</msub>
<msub>
<mi>T</mi>
<mi>k</mi>
</msub>
<mi>&delta;</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>l</mi>
<mo>-</mo>
<mi>v</mi>
<mi>t</mi>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>Y</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>c</mi>
<mi>k</mi>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mrow>
<mo>(</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>Q</mi>
</munderover>
<msub>
<mi>c</mi>
<mi>k</mi>
</msub>
<msub>
<mi>H</mi>
<mi>L</mi>
</msub>
<msub>
<mi>T</mi>
<mi>k</mi>
</msub>
<mi>&delta;</mi>
<mo>(</mo>
<mrow>
<mi>l</mi>
<mo>-</mo>
<mi>v</mi>
<mi>t</mi>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>c</mi>
<mi>k</mi>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<msub>
<mi>H</mi>
<mi>L</mi>
</msub>
<msub>
<mi>T</mi>
<mi>k</mi>
</msub>
<mi>&delta;</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>-</mo>
<mi>v</mi>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
其中,和分别表示结构加速度响应对单元刚度损伤参数和移动力正交参数的一阶灵敏度矩阵,表示结构马尔科夫系数关于单元刚度损伤参数的一阶偏导,具体计算如下:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>H</mi>
<mn>0</mn>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<mn>0</mn>
<mo>,</mo>
<mi>m</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>H</mi>
<mi>m</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mrow>
<mo>(</mo>
<mrow>
<msup>
<mi>RA</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
</mrow>
</msup>
<mi>B</mi>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>R</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<msup>
<mi>A</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
</mrow>
</msup>
<mi>B</mi>
<mo>+</mo>
<mi>R</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mrow>
<mo>(</mo>
<msup>
<mi>A</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
</mrow>
</msup>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<mi>B</mi>
<mo>+</mo>
<msup>
<mi>RA</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
</mrow>
</msup>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>B</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<mo>,</mo>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中和的计算如下:
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>R</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>R</mi>
<mi>a</mi>
</msub>
<msup>
<mi>M</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>K</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>R</mi>
<mi>a</mi>
</msub>
<msup>
<mi>M</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>C</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>R</mi>
<mi>a</mi>
</msub>
<msup>
<mi>M</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>K</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>bR</mi>
<mi>a</mi>
</msub>
<msup>
<mi>M</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>K</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mrow>
<mo>(</mo>
<msup>
<mi>A</mi>
<mrow>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mrow>
<mo>(</mo>
<msup>
<mi>A</mi>
<mrow>
<mi>m</mi>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msup>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<mi>A</mi>
<mo>+</mo>
<msup>
<mi>A</mi>
<mrow>
<mi>m</mi>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msup>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>A</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
</mrow>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>A</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<mi>A</mi>
<mi>h</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<msup>
<mi>K</mi>
<mo>*</mo>
</msup>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
</mrow>
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>B</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mo>&part;</mo>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<msup>
<mi>K</mi>
<mrow>
<mo>*</mo>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mrow>
<mi>A</mi>
<mo>-</mo>
<mi>I</mi>
</mrow>
<mo>)</mo>
</mrow>
<msup>
<mi>B</mi>
<mo>*</mo>
</msup>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msup>
<mi>K</mi>
<mrow>
<mo>*</mo>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<mrow>
<mi>A</mi>
<mo>-</mo>
<mi>I</mi>
</mrow>
<mo>)</mo>
</mrow>
<msup>
<mi>B</mi>
<mo>*</mo>
</msup>
<mo>+</mo>
<msup>
<mi>K</mi>
<mrow>
<mo>*</mo>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>A</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<msup>
<mi>B</mi>
<mo>*</mo>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mrow>
<mo>-</mo>
<msup>
<mi>K</mi>
<mrow>
<mo>*</mo>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mfrac>
<mrow>
<mo>&part;</mo>
<msup>
<mi>K</mi>
<mo>*</mo>
</msup>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<msup>
<mi>K</mi>
<mrow>
<mo>*</mo>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mrow>
<mi>A</mi>
<mo>-</mo>
<mi>I</mi>
</mrow>
<mo>)</mo>
</mrow>
<msup>
<mi>B</mi>
<mo>*</mo>
</msup>
<mo>+</mo>
<msup>
<mi>K</mi>
<mrow>
<mo>*</mo>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mi>A</mi>
<mi>h</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<msup>
<mi>K</mi>
<mo>*</mo>
</msup>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<msup>
<mi>B</mi>
<mo>*</mo>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<msup>
<mi>K</mi>
<mo>*</mo>
</msup>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msup>
<mi>M</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>K</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msup>
<mi>M</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>C</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msup>
<mi>M</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>K</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msup>
<mi>bM</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>K</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>.</mo>
</mrow>
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410527720.3A CN104298870B (zh) | 2014-10-08 | 2014-10-08 | 一种移动荷载下简支梁损伤和移动力同时识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410527720.3A CN104298870B (zh) | 2014-10-08 | 2014-10-08 | 一种移动荷载下简支梁损伤和移动力同时识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104298870A CN104298870A (zh) | 2015-01-21 |
CN104298870B true CN104298870B (zh) | 2018-01-26 |
Family
ID=52318593
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410527720.3A Active CN104298870B (zh) | 2014-10-08 | 2014-10-08 | 一种移动荷载下简支梁损伤和移动力同时识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104298870B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104598931B (zh) * | 2015-02-06 | 2016-10-05 | 华北水利水电大学 | 基于Kaczmarz代数迭代重建方法的桥梁移动车辆荷载识别方法 |
CN105675493A (zh) * | 2016-03-10 | 2016-06-15 | 吉林建筑大学 | 移动荷载作用下三跨连续梁桥的损伤识别方法 |
CN106202789B (zh) * | 2016-07-20 | 2019-04-05 | 暨南大学 | 一种基于l1/2正则化的桥梁移动荷载与损伤协同识别方法 |
CN106874583A (zh) * | 2017-01-25 | 2017-06-20 | 华北水利水电大学 | 基于对角松弛正交投影迭代算法的桥面多轴移动荷载的识别方法 |
CN107633132B (zh) * | 2017-09-18 | 2020-11-06 | 湖南大学 | 空间分布动态载荷的等效识别方法 |
CN107704719A (zh) * | 2017-11-27 | 2018-02-16 | 中国铁道科学研究院铁道建筑研究所 | 一种基于灵敏度的桥梁整体损伤识别方法 |
CN110017929B (zh) * | 2019-03-12 | 2020-12-08 | 华中科技大学 | 基于子结构灵敏度分析的船撞桥荷载与损伤同步识别方法 |
CN112964165B (zh) * | 2021-02-08 | 2021-12-14 | 合肥工业大学 | 一种移动荷载作用下的桥梁位移重构方法 |
CN113094928B (zh) * | 2021-04-26 | 2022-11-04 | 东北农业大学 | 一种人行单步荷载模型构建方法 |
CN113591233B (zh) * | 2021-06-03 | 2024-04-30 | 南京林业大学 | 一种用于求解柔性梁系统动力学方程的分裂迭代方法 |
CN117610307B (zh) * | 2023-12-15 | 2024-05-17 | 大连海事大学 | 一种移动质量作用下简支梁的数字孪生构建方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102080955A (zh) * | 2010-12-02 | 2011-06-01 | 中山大学 | 一种用于桥梁安全检查的移动荷载定位方法 |
CN102620900A (zh) * | 2012-04-01 | 2012-08-01 | 长安大学 | 基于动载试验检测桥梁冲击系数的方法 |
CN103674598A (zh) * | 2013-12-18 | 2014-03-26 | 西南交通大学 | 模拟移动荷载的轨道结构疲劳试验方法及装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103063451B (zh) * | 2012-12-26 | 2014-07-23 | 浙江大学 | 轨道交通列车整车移动荷载模拟加载方法及装置 |
-
2014
- 2014-10-08 CN CN201410527720.3A patent/CN104298870B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102080955A (zh) * | 2010-12-02 | 2011-06-01 | 中山大学 | 一种用于桥梁安全检查的移动荷载定位方法 |
CN102620900A (zh) * | 2012-04-01 | 2012-08-01 | 长安大学 | 基于动载试验检测桥梁冲击系数的方法 |
CN103674598A (zh) * | 2013-12-18 | 2014-03-26 | 西南交通大学 | 模拟移动荷载的轨道结构疲劳试验方法及装置 |
Non-Patent Citations (2)
Title |
---|
移动荷载作用下桥梁动态响应的波谱单元法;张俊兵 等;《土木工程与管理学报》;20110915;第28卷(第3期);第206-211页 * |
高速移动荷载下轨道系统振动模拟的数值算法;谢伟平 等;《武汉理工大学学报》;20011225;第23卷(第12期);第56-59页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104298870A (zh) | 2015-01-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104298870B (zh) | 一种移动荷载下简支梁损伤和移动力同时识别方法 | |
Oshima et al. | Damage assessment of a bridge based on mode shapes estimated by responses of passing vehicles | |
Shu et al. | The application of a damage detection method using Artificial Neural Network and train-induced vibrations on a simplified railway bridge model | |
Kim et al. | Modal parameter identification of short span bridges under a moving vehicle by means of multivariate AR model | |
Xu et al. | Energy damage detection strategy based on acceleration responses for long-span bridge structures | |
Sung et al. | Damage detection for beam-like structures using the normalized curvature of a uniform load surface | |
Hernandez-Garcia et al. | An experimental investigation of change detection in uncertain chain-like systems | |
Keenahan et al. | Determination of road profile using multiple passing vehicle measurements | |
Zhu et al. | A rapid structural damage detection method using integrated ANFIS and interval modeling technique | |
Park et al. | Sequential damage detection approaches for beams using time-modal features and artificial neural networks | |
CN106202789A (zh) | 一种基于l1/2正则化的梁桥移动荷载与损伤协同识别方法 | |
Ni et al. | Structural damage detection using auto/cross-correlation functions under multiple unknown excitations | |
CN103884776B (zh) | 一种提高随机损伤定位向量法监测结果准确性的方法 | |
CN104598753A (zh) | 一种基于Brakhage v方法的桥梁移动车辆荷载识别方法 | |
JP2017044614A (ja) | 構造物の伝達関数の高精度評価・地震応答予測・劣化診断システム及びその方法 | |
Jayawardhana et al. | Statistical damage sensitive feature for structural damage detection using AR model coefficients | |
George et al. | Damage detection in railway bridges under moving train load | |
Daneshvar et al. | Structural health monitoring using high-dimensional features from time series modeling by innovative hybrid distance-based methods | |
Mahmood et al. | A multi‐reference‐based mode selection approach for the implementation of NExT–ERA in modal‐based damage detection | |
Damadipour et al. | Seismic modal identification using a new approach based on weighted transmissibility | |
CN103530275A (zh) | 基于振动传递率函数主成分置信度的结构损伤预警方法 | |
KR20150004127A (ko) | 가속도 및 변형률 측정값에 기반한 구조물의 다 지점 변위 추정 방법 | |
Golmohamadi et al. | Damage diagnosis in bridges using wavelet | |
Gul et al. | A modified time series analysis for identification, localization, and quantification of damage | |
Zhan et al. | A double-pass method for bridge assessment considering surface roughness using normalized contact point responses |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |