近年来，我国桥梁趋于大跨径、轻柔化及形式与功能的复杂化发展，因施工质量问题(如麻面、空洞、缝隙夹层)，或在长期使用过程中，由于外界因素(如车辆荷载作用、环境腐蚀效应、材料老化及突发自然灾害)的影响，引起钢筋混凝土梁桥出现开裂、断裂、缺边、掉角、分化、剥落等不同程度的损伤，其微小损伤的不断积累使桥梁结构存在较大的安全隐患，一旦发生事故，会造成巨大的人员伤亡和财产损失，给国民经济和社会发展带来不利后果.为了及时掌握桥梁运营下的健康状况，基于振动测试的损伤检测技术引起了国内外工程界和学术界广泛关注[1-2]，但受到结构参数误差、测试参数误差、模型建立误差和现场测量误差等不确定性因素的限制，而弹性波携带介质变化信息在结构中传播，使结构损伤检测成为可能.

文献[3-5]首先提出了基于Timoshenko梁理论的回传射线矩阵法(简称MRRM法)；余云燕等[6]应用该方法研究了有损伤变模量连续梁的时域瞬态响应；谢禹等[7]应用回传射线矩阵法分析有缺陷变模量简支梁的时域瞬态响应，并进行损伤诊断研究；Huang等[8-9]首先将复杂时域信号使用经验模态分解(EMD)，运用Hilbert变换构造信号的时频分布特征；也提出了许多结构损伤识别方法[10-12].但是针对结构损伤检测研究只限于单一的时域、频域或时频分析，也缺乏运用统一的动力学理论模型对结构损伤信息的时域、频域、时频域对比分析，并且结构的局部损伤用用减小单元的杨氏模量模拟，在一定程度上造成波动响应信号的延滞效应，影响损伤信息的准确诊断.本文将回传射线矩阵法应用于变截面简支梁桥的损伤检测，对比分析了有损、无损伤简支梁桥在脉冲荷载作用下速度波的时域、频域损伤信号，对已求解的时域信号运用Hilbert-Huang变换进行损伤信息的对比分析，为工程实践提供理论基础.

但是，弹一个“极弱音”很容易，困难的是弹一组连续的“极弱音”，因为一连串质地类似的“极弱音”要求我们在弹奏每一个音之前，都预备好相应的“慢动作”。如果在弹奏第一个音时，第二个音所需的手指没有开始“慢动作”，那就无法弹出下一个弱音，或只能以极慢的速度弹奏下一个弱音。

1 有损伤简支梁桥的动力学模型

基于Timoshenko梁理论模型，矩形简支梁桥局部损伤用减小横截面面积来模拟.对该简支梁桥在局部坐标系下划分单元，并进行节点编号，每个单元用其两端的节点编号表示，如图1(a)所示；在节点2上作用与水平方向成α角的方波脉冲激振力，如图1(b)所示.

  

图1 有损伤简支梁桥结构与脉冲荷载Fig.1 Structure of damaged simple beam bridge and impulsive load

对每个梁单元分别建立两个局部坐标系(x,y)JK和(x,y)KJ,在局部坐标系下，无损单元12、23、45的控制方程为

 

(1)

 

(2)

式中：u为轴向位移；v=vb+vs为简支梁桥挠度；vb由弯矩引起的挠度；vs由剪力引起的挠度；A、E、G、ρ、Iz、k′分别为无损简支梁桥的横截面面积、弹性模量、剪切模量、密度、横截面惯性矩和截面剪力系数.

2.3 2组患儿治疗前后气道阻力指标比较 2组患儿治疗后气道阻力与治疗前比较有降低（P＜0.05），治疗组治疗后低于对照组（P＜0.05）。结果见表4。

令：为纵波波速，为单元的回转半径，

302 Advances in animal models of chronic obstructive pulmonary disease

将式(1)～(2)进行无量纲化并进行Fourier变换得频率中的方程：

 

(3)

 

(4)

求解式(3)～(4)得：

 

(5)

 

(6)

 

(7)

总挠曲位移：

 

(8)

式中：a1(ω)、a2(ω)、a3(ω)为待定的入射波波幅；d1(ω)、d2(ω)、d3(ω)为待定的出射波波幅；k1、k2、k3为波数，即：k1(ω)=ω,

 

对应于波数与的比值为

 

时域中轴力F、剪力Q、弯矩Mz、转角δz、纵波Cu和横波Cv分别表示为

 

(9)

将式(9)进行Fourier变换，并通过式(5)～(8)可求得轴力、剪力、弯矩、转角、纵波和横波在频域中的表达式.

但是，当ω取某些特定频域时，(I-R)-1矩阵不存在，采用Neumann级数展开来避开奇异点，即：

式中：为各杆无量纲长度.

aJ=SJaJ+sJ.

(10)

式中：aJ、dJ表示J点入射波和出射波波幅向量；SJ、sJ为J点的局部散射矩阵和局部源矢量(外荷载列向量).

aN=PU[I-R+R2+…+RN-1+…]s=PUdN-1.

  

图2 节点3处的局部坐标系和整体坐标系Fig.2 Local and global coordinate system at joint 3

列出节点3处的位移协调及力平衡条件为

在大数据的时代背景下，计算机网络之中蕴藏了大量的信息内容，而在其开放的系统之中，这些信息的安全性大大降低，常常会出现信息的泄露与被破坏的情况。一般而言，导致信息泄露主要是不法分子通过非法手段对信息进行窃取，利用系统的故障或是病毒对用户的计算机进行攻击。

 

(11)

将式(11)代入式(9)，整理合并为

d3=S3a3+s3.

式中：

 

S3=

 

s3=(0 0 0 0 0 0)T.

Z(t)=Z(t)real+iZ(t)imag=X(t)+iH(t)=T(t)eiθ(t).

d=Sa+s.

(12)

在每个单元的局部坐标系下，入射波波幅向量和出射波波幅向量两者之间存在一个相位关系，引入传播矩阵PJK，进而得到方程式：

aJK=PJKdJK.

在频域中，对所有各节点J建立局部坐标系下的力平衡和位移协调条件，经过整理得：

将aJK，dJK分别组集到总体矩阵中，写成紧凑形式为

 

由于与d只在元素的排列顺序上不同，故引入置换矩阵U，U中只有0、1元素，得则：

a=PUd.

(13)

将式(13)代入式(12)整理得：

d=(I-R)-1s.

(14)

式中：R=SPU称为结构的回传射线矩阵；I为单位矩阵.

损伤单元34的波动方程，只须将控制方程(1)～(2)中的A、Iz改为A1、Iz1，并引入损伤程度越小，损伤程度越大，同理可求解有损单元频域中各物理量的表达式.

(I-R)-1=(I-R+R2+…+RN+…).

(15)

将式(15)中(I-R)-1取前N项代入式(12)得：

dN=[I-R+R2+…+RN+…]s.

(16)

将式(16)代入式(13)，就可求出aN，即：

只在节点2上作用有与水平方向成α角的方波脉冲激振力，故除节点2处之外，其余节点处的局部源矢量为零矢量.对于节点1、5只连接一个单元，节点2、3、4连接两个单元，故S1、S5是三阶方阵，S2、S3、S4是六阶方阵.在局部坐标系和整体坐标系下，以频域中节点3为例，如图2所示.

(17)

式中：RNs表示波经过了N次节点回传，从上式可以看出，由于局部坐标系采用了对偶坐标系，所以入射波较出射波多传播一次.

从表3可知，有损伤简支梁中速度波的第一条射线到达无量纲时间和无损伤的第一条射线相同.第二条、第四条射线为经过有损简支梁(节点3和4)的散射回传至接收点.第五条射线传播路径为2→3→2→1→2→A，纵波到达接收点A的无量纲时间为28.

求出a、d以后，可求出频域中任意一点物理量，通过Fourier逆变换(IFFT)可得相应物理量在时域中的瞬态相应，以单元上任意一点经过N次节点回传后无量纲的横波为例，其表达式为

 

(18)

对求出频率中的任意物理量可以分解成实部和虚部，以横波为例，用极坐标形式表示成振幅及相位φ(ω)为

使得：yTα=0，其中a为拉格朗日乘子(0≤αi≤C)；e表示数值全为1的向量；Q是l乘l的半正定矩阵，Qij≡yiyjK(xi,xj)，K(xi,xj)≡ϕ(xi)Tϕ(xj)，K为核函数，ϕ是训练向量低xi维空间到高维空间的映射函数。最终得到二分类训练决策模型：

 

(19)

 

(20)

2 Hilbert-Huang变换方法

Hilbert-Huang变换(简称HHT变换)，HHT变换的信号处理方法是由Huang等[8-9]于1998年提出，首先是使用经验模态分解(Empirical Mode Decomposition，简称EMD)，将给定的非线性、非平稳信号自适应地分解为一系列若干个表征信号特征时间尺度的固有模态函数(intrinsic mode function，简称IMF).分解过程通过一种称为“筛选”[12]的计算步骤实现，即所谓的自适应分解基，分解出n的个IMF(B1(t)，B2(t)，…，Bn(t))包含了原信号从高频段到低频段的成分，原始信号将是一组IMF分量与残余项rn(t)的总和，即：

过去矿业用地租赁集体土地的情况较为频繁，多为农用地或者未利用地，其实际用途与二调现状用途并不一致，难以与现有政策衔接，造成后续使用困难。因此建议办理建设用地转用手续，但是具体操作中，根据其与二调现状用途的出入情况，进行差别化处理：一是对于二调已确认为建设用地的，按照建设用地管理，需要处罚的，依法处罚到位后，为用地单位补办建设用地手续；对于二调确定为非建设用地的，实际上因矿山开采造成损毁不可恢复的，经年度变更调查、核减耕地保有量、补划基本农田后，办理建设用地转用手续。

 

(21)

对给定的连续原始信号X(t)，Hilbert变换H(t)是将信号X(t)与脉冲响应Y(t)=1/πt做卷积，即：

 

(22)

对应于X(t)的解析信号Z(t)为

通过对测定指标的统计分析、相关性分析和因子评价，确定了鲜切黄瓜加工适宜性的评价指标为口感、果肉厚度、TSSC、果肉硬度、果皮硬度、果肉色差a*和L*值；通过层次分析确定各主评价指标的权重分别为35.45%、23.66%、15.82%、10.50%、6.92%、4.59%和3.07%；应用灰色关联分析对19个黄瓜样品材料的鲜切加工适宜性进行了排序，结合各样品的田间栽培性状优选出较适宜的鲜切加工样品材料为4#和2#，表明该评价体系能够较合理的指导黄瓜田间育种，筛选适宜的黄瓜鲜切加工品种。

还有一个决定性的因素，跟讲义夹里面的内容有关。三年多，我父亲留下的账页有三四十张，每行记着日期、村名、巷名和东家姓名。正反两面都是，密密麻麻的。李锦文翻着账页，眼睛红了，带着哭腔问我：晓得叫你将来一一去拜访人家是什么意思吗？我摇头。李打油把他那段打油诗再背一遍，就是“腋下大夹子”那首。再问，我还是摇头。他急了，他说你设身处地想想，每天跟在猪牯后面处处去、为了三块钱家家候，那是什么感觉？

将所有节点的a、d、S、s组集到整体矩阵中，得到整体坐标系下的整体散射矩阵和整体源矢量，即：

(23)

式中：T(t)、θ(t)分别为原始信号的瞬时振幅和瞬时相位，

由瞬时相位的导数可得瞬时频率的表达式为

Φ(t)=dθ(t)/dt.

(24)

瞬时频率Φ(t)是时间t的函数，对每一个IMF做Hilbert变换，求出每一个IMF所对应的瞬时频率，即：

 

(25)

H(t)既是时间的函数，又是瞬时频率的函数，而瞬时频率也是时间的函数，取上式的实部，得到原始信号X(t)的Hilbert谱，反映信号幅值随时间和频率的变化情况，记为H(ω,t).将Hilbert谱对时间积分，得到Hilbert边际谱，即：

h(ω)=H(ω,t)dt.

(26)

3 基于MRRM与HHT法有损简支梁桥检测的对比分析

如图1(a)所示的有损简支梁桥，设节点2为脉冲荷载作用点，A点为信号接收点，距离脉冲作用点的无量纲长度为4，损伤位置分别为3和4节点，节点编号如图1(a)所示.无损单元和有损单元的横截面面积、惯性矩分别为A、Iz和A1、Iz1，损伤程度m=A1/A=0.5，杨氏模量与修正剪切模量之比n=2.99，瞬态响应展开项级数N=8，脉冲荷载为方波脉冲，无量纲化脉冲持续时间如图1(b)所示.本文的计算数据采用无量纲方法，简支梁桥无量纲总长度无损各单元的无量纲长度分别为有损单元的无量纲长度为以杆件12单元为例，其纵波传播无量纲时间横波传播的无量纲时间为其中:为纵波波速，为横波波速，其余杆件的求法类似，各杆件的无量纲长度及纵波、横波传播的无量纲时间如表1所列.

通常情况下，不同的地区有着各自不同的土质，不同土质的稳定性也有着很大的差距，例如，陕西地区就属于黄土土质，又因为大多处于特殊的地质构造处，所以极其容易发生地面沉降的事故。另外，陕西地区的湿陷性黄土地基也较为严重，对地基处理不当是造成地基基础下沉的主要原因[3]。

 

表1 无量纲长度及速度波传播的无量纲时间

 

Tab.1 Dimensionless length and dimensionless time for the movementof axial velocity waves and transverse velocity waves

  

杆件名称杆件无量纲长度纵波传播无量纲时间横波传播无量纲时间杆12446．928杆238813．856杆34446．928杆45446．928杆2A446．928杆3A446．928

3.1 基于MRRM法对简支梁桥有损伤和无损伤纵波、横波的时域对比分析

对于无损伤简支梁桥，由脉冲作用点(图1中的节点2)出射的波向节点1和接收点A传播，第一条射线传播路径为2→A，纵波到达接收点A的无量纲时间为4，其首波到达时间也很准确，与预测值一致.第二条射线传播路径为2→5→A，纵波到达接收点A的无量纲时间为28，由于经过铰节点5的反射，反射后波的脉冲方向改变一次，波形不变.第三条射线传播路径为2→1→2→A，纵波到达接收点A的无量纲时间为12，也经过了铰节点1的反射，脉冲方向随之改变，波形不变，反映了节点约束条件.第四条射线传播路径为2→1→2→5→A，纵波到达接收点的无量纲时间为36.第五条射线传播路径为2→5→2→1→2→A，纵波到达接收点A的无量纲时间为44.无损伤简支梁中纵波、横波的回传次数、传播路径及到达时间如表2所列.

 

表2 速度波的回传次数、传播路径及到达时间

 

Tab.2 Times of return, movement path and arrival time ofaxial velocity waves and transverse velocity waves

  

射线传播路径纵波传播无量纲时间t′i横波传播无量纲时间t′1i第一条2→A46．928第二条2→5→A2848．497第三条2→1→2→A1220．785第四条2→1→2→5→A3662．354第五条2→5→2→1→2→A4476．210

表3为有损伤简支梁中纵波、横波的回传次数、传播路径及到达时间.因表3中所列的前五条回传射线可以准确定量判定缺陷，故后面的传播路线不再罗列.

 

表3 有损伤简支梁中速度波的回传次数、传播路径及到达时间

 

Tab.3 Times of return, movement path and arrival time ofaxial velocity waves and transverse velocity wavesin the damaged simple beam bridges

  

射线传播路径纵波传播无量纲时间ti横波传播无量纲时间t1i第一条2→A46．928第二条2→3→A1220．785第三条2→1→2→A1220．785第四条2→1→2→3→A2034．641(2→3→4→3→A)2034．641第五条2→3→2→1→2→A2848．497

本设计方案最大的创新之处在于身份识别系统基于以太坊智能合约技术，立足于传统互联网行业的现状，解决了中心化管理的麻烦，以及用户信息容易被篡改、被黑客盗用，中心服务器被攻击等事关国家信息安全痛点的问题。

图3为在方波脉冲激振力作用下，有损伤简支梁桥和无损伤简支梁桥中纵波的对比.根据表2和表3的分析，若没有损伤存在，A点纵波在无量纲时间12～28无其他射线到达且波的形状保持不变，在无量纲时间12时纵波波幅减小(脉冲持续时间不变)，且无量纲时间20时又有一条射线到达，说明简支梁桥存在损伤.在无量纲时间12时纵波传播路径为2→3→A(2→1→2→A)，经过损伤处(图1中节点3)的反射和透射，使无量纲纵波波幅减小.在无量纲时间20时这条射线的传播路径为2→1→2→3→A(2→3→4→3→A)，这是由于简支梁桥存在损伤，纵波经过损伤处(图1中节点3和4)的散射，纵波波幅减小，由于没有经过铰节点，所以波的脉冲方向与无量纲时间12时的脉冲方向相同.对于结构距接收点A的损伤位置，可根据表3中第一条射线2→A和第二条射线2→3→A确定，设简支梁桥的横截面面积A为0.2 m2，即：根据表3中第二条射线(2→3→A)和第四条射线(2→3→4→3→A)确定损伤区域，即：

  

图3 有无损伤简支梁桥在接收点A的纵波比较Fig.3 Comparison of axial velocity waves at receiving spot A between damaged and non-damaged simple beam bridges

图4为有损伤和无损伤简支梁桥在接收点A的横波比较结果.

另外，从立法模式上讲，“我国关于偷税罪的规定实际上采取的是一种混合式立法模式，即在税法如《税收征管法》里规定涉及偷税罪的有关行为和表现，但真正明确罪刑关系是在刑法典中。因此，我国偷税罪立法的改革完善还要取决于立法模式的改变。从世界发展趋势来看，随着各国经济的发展、税制得复杂化、税种的不断增加，把税收犯罪规定于各种税收管理法规之下是一种理想选择。”或许，该选择是目前一些经济文化发达国家税法发展的主要方向，但是对于我国现今经济文化与水平来说，或许不是最佳选择，我国仍需如当年一样“摸着石头过河”来探索最适合国情的立法方向。

  

图4 有无损伤简支梁桥在接收点A的横波比较Fig.4 Comparison of transverse velocity waves at receiving spot A between damaged and non-damaged simple beam bridges

从图4可以看出，在横波到达接收点A的无量纲时间为6.928时，其传播路径为2→A，由于没有经过损伤区域，两条曲线完全重合，之后两条曲线在无量纲时间20.8时开始分离，横波传播路径为2→3→A，这是由于简支梁桥中的横波在损伤位置(节点3)处发生节点反射和透射，使无量纲横波波幅明显减小.根据表3中横波传播的第一条射线2→A和第二条射线2→3→A，可以确定简支梁桥存在损伤，因横波在结构中横向传播，不能通过横波波速与无量纲时间的乘积来确定损伤位置及损伤区域.由于横波的每条射线到达接收点A的信号没有纵波强烈，且首波不是方波，所以横波用于结构检测时存在局限性.以上通过对纵波和横波的对比分析发现，纵波用于结构检测更具有优越性，在下面的讨论中以纵波的时频、频域进行对比分析.

3.2 基于MRRM法对有损伤和无损伤简支梁桥纵波波谱的频域对比分析

图5为有损伤和无损伤简支梁桥横波在接收点A的相位谱频域对比结果.从图5中可以看到，在频域里，有损伤纵波的相位谱变化信号不强烈，简支梁桥损伤很难判别.

  

图5 有无损伤简支梁桥纵波在接收点A的相位谱比较Fig.5 Phase spectrum comparison of axial velocity waves at receiving spot A between damaged and non- damaged simple beam bridges

图6为有损伤和无损伤简支梁桥横波在接收点A的振幅谱对比结果.

  

图6 有无损伤简支梁桥纵波在接收点A的振幅谱比较Fig.6 Amplitude spectrum comparison of axial velocity waves at receiving spot A between damaged and non-damaged simple beam bridges

图6体现了一定的规律性，从频域里可以看到，有损伤简支梁桥横波在接收点A的振幅谱波峰值比无损伤的低，并且在对应的两峰值之间，有损伤简支梁桥横波的振幅谱有明显的变化结果，可以判别简支梁桥存在损伤，但是对损伤位置和损伤区域很难判断.

3.3 用HHT法对有损伤简支梁桥纵波的时频分析

图7为有损伤简支梁桥纵波在接收点A的IMF2瞬时频率.接收点A在图7中归一化时间和表3中的无量纲时间一致；在方波脉冲荷载作用下，首波到达接收点A的归一化时间为4，对应的瞬时频率为20 Hz；在归一化时间12时，瞬时频率发生突变，经过突变以后，对应的瞬时频率为14 Hz，由此可以判定简支梁桥存在损伤，波传播经过损伤位置(图1中的节点3)时，发生了散射，在频域里表现为瞬时频率的减小；在量纲一时间20、28时，是由于波在传播过程中相遇、叠加，使瞬时频率幅值变化很小；在量纲一时间36时，瞬时频率减小至12 Hz，根据量纲一时间12、36，可确定简支梁桥损伤的位置和损伤区域.

  

图7 有损伤简支梁桥纵波在接收点A的IMF2瞬时频率Fig.7 The IMF2 instantaneous frequency of axial velocity waves at receiving spot A in damaged simple beam bridges

图8为有损伤简支梁桥纵波在接收点A的Hilbert边际谱.对有损伤简支梁桥在接收点A的纵波进行了Hilbert变换，得到响应信号的Hilbert边际谱，而Hilbert边际谱中某一频率f处能量的存在，就表示一定有这样的振动波在局部出现过.

在图2中比较了ESD和MSD的计算复杂度与Rx-SD，Tx-SD的计算复杂度.图中的复杂度定义为Cre(%)=100×(CML-CSD)/CML[9]，表示SD相对于ML检测减小的相关复杂度.图2(a)中可以看出，在相同的频谱效率(频谱效率为6 bps/Hz)下，Nt>Nr时，改进算法比ML检测和经典的SD检测具有更低的计算复杂度，尤其是在信噪比较低的情况下优势更加明显，在信噪比为3dB时，ESD和MSD的相对复杂度降低了33%.图2(b)中，Nt≤Nr时，与Tx-SD进行比较，MSD的相关复杂度降低了20%以上.此外，在相同条件下ESD与MSD相比，ESD能够降低约2%的相关复杂度.

图8中表现出来的频率比较丰富，可以清晰的看到信号能量在高频和低频时的分布情况，当简支梁桥受到脉冲荷载作用时，结构中的高频部分被激起，且随时间变化.如图8中归一化时间为4时，接收点A接收到瞬态波动能量信号；随着时间的变化，在归一化时间12时，有损伤简支梁桥横波在接收点A的高频变化很低，是因为在损伤处(图1中节点3)，发生了波的反射和透射.在归一化时间20、28、36、44时，使信号能量减小的过程，是因为波在往返传播过程中，在接收点A处形成了波的叠加，并且波传播至缺损处(图1中的节点3、4)，使波不断的反射和透射，使能量逐渐衰减的过程.在归一化时间52时，被激起的高频部分靠近结构的低频部分，这个直观的图形是Fourier谱无法展示的，Fourier谱只能反映实际频率成分的大致情况，而Hilbert边际谱能够真实反映频率在信号中是否存在.对于结构损伤检测，时频分析比时域分析或频域分析，更具有直观性、优越性和不可替代性.图8展示的归一化时间与表3中纵波传播无量纲时间一致，根据表3和图8中波的传播路径和归一化时间，可容易确定简支梁桥存在损伤、损伤的位置和损伤区域.

  

图8 有损伤简支梁桥纵波在接收点A的Hilbert边际谱Fig.8 The Hilbert marginal spectrum of axial velocity waves at receiving spot A in damaged simple beam bridges

4 结论

本文基于回传射线矩阵法和HHT变换，研究了有损伤简支梁桥在方波脉冲荷载作用下的瞬态波动响应，对接收点速度波的时域、频域和时频进行了对比分析，得出以下结论：

1) 应用回传射线矩阵法进行时域分析，纵波可判定简支梁桥存在损伤，确定损伤位置、损伤区域；横波可判定简支梁桥存在损伤，不能确定损伤位置及损伤区域.

2) 应用回传射线矩阵法进行频域分析，纵波的振幅谱，功率谱可判定简支梁桥存在损伤，不能确定损伤位置、损伤区域；相位谱不能用于结构检测.

3) 应用HHT法进行时频分析，纵波可以清晰、直观的确定简支梁桥存在损伤、损伤位置及损伤区域.

二是秉承山水林田湖草“生命共同体”理念，树立“大土地观”，打造土地质量地质调查的升级版。要充分发挥地球化学测量多目标、高精度和信息化程度高的技术优势，不断拓展调查的深度和广度，不断丰富土地质量的内涵，不断创新成果的表达方式。在问题和需求的引领下，着力做好知识储备、技术储备、人才储备，快速提升服务支撑能力，以满足自然资源管理对科技工作的需求。

4) 应用回传射线矩阵法和HHT变换，检测方法简单，且易于编程，结果准确，只需得到方波脉冲荷载作用下接收点的速度波响应，无需结构的模态参数分析，对于多处损伤的复杂结构的损伤检测，仅增加程序计算时间.

参考文献：

[1] 谭林.基于动力指纹的结构损伤识别可靠度方法研究[D].广州:华南理工大学,2010.

[2] 王名月,缪炳荣,李旭娟,等.基于转角模态和小波神经网络的连续梁损伤识别研究[J].力学季刊,2016,37(4):686-691.

[3] PAO Y H,KEH D C,HOWARD S M.Dynamic response and wave propagation in plane trusses and frames[J].AIAA.Journal,1999,37(5):594-603.

[4] HOWARD S M,PAO Y H.Analysis and experiments on stress waves in planar trusses[J].Journal of Engineering Mechanics(ASCE),1998,124(8):884-891.

[5] CHEN J F,PAO Y H.Effects of causality and joint conditions on method of reverberation-ray matrix[J].AIAA Journa1,2003,41(6):1138-1142.

[6] 余云燕.有损伤连续梁的瞬态响应[J].爆炸与冲击,2009,29(1):61-66.

[7] 谢禹,余云燕.有损伤简支桥梁的波动分析[J].兰州交通大学学报,2006,25(6):74-77.

[8] HUANG N E,ZHENG S,LONG S R,et al.The empirical mode decomposition and Hilbert spectrum for nonlinear and non-stationary time series analysis[J].Proceedings of the Royal Society of London, Series A,1998,454:903-995.

[9] HUANG N E,SHEN Z,LONG S R.A new view of nonlinear water waves-hilbert spectrum[J].Annual Review of Fluid Mechanics,1999,31:417-457.

[10] 杨智春,张慕宇,丁燕,等.Hilbert边际谱在框架结构损伤检测中的应用[J].振动与冲击,2010,26(11):250-255.

[11] 任宜春,翁璞.基于改进 Hilbert-Huang变换的结构损伤识别方法研究[J].振动与冲击,2015,34(18):195-199.

[12] 冯志鹏,秦嗣峰.基于Hilbert振动分解和高阶能量算子的行星齿轮箱故障诊断研究[J].振动与冲击,2016,35(5):47-54.