实际做结构动力测试这些年，最常被问到的问题就是：“测到的时变位移到底准不准？”“频响函数算出来怎么跟书上对不上？”说实话，这两个问题背后藏着的是一整套关于测试方案、信号处理和参数识别的门道。结构在风、车辆、地震这类环境荷载作用下，它的响应是随时在变的——我们把这叫时变响应。而要摸清结构自身的动态特性，就必须借助频响函数这把标尺。这篇文章我就结合自己动手做过的测试项目，把结构响应时变位移的测量与分析、频响函数估计的完整链条拆开讲一遍，包括中间踩过的坑和总结出来的经验。
1. 项目核心概念拆解：时变位移与频响函数的关系
1.1 时变位移到底在测什么
先说说时变位移。所谓时变位移，就是结构上某一点的位置随时间变化的量，通常用绝对位移或者相对位移表示。我们在桥面上、楼板上、大型设备基础上安装传感器，测到的物理量有可能是加速度、速度，也有可能是位移。但工程上真正关心的往往是位移——因为位移直接对应变形，变形大了就可能开裂、屈曲甚至倒塌。
我见过很多刚入行的工程师，拿到一支加速度计就只认加速度数据，忽略了积分位移的用途。实际上，加速度时程经过一次积分得到速度，经过二次积分得到位移，但这里有个大坑：低频漂移。现场测试时传感器自身的零点漂移、温度变化、线缆扰动都会在低频段产生虚假成分，而积分运算恰恰会把低频成分无限放大。结果就是位移时程像一条歪歪扭扭的龙，完全看不出真实的结构响应形状。这也是为什么业内常说“加速度计测位移，靠的不是硬件，是信号处理”。
时变位移的典型应用场景包括：大跨度桥梁在车辆荷载下的竖向挠度跟踪、高层建筑在脉动风作用下的层间位移角评估、以及大型旋转机械基础在启停机过程中的瞬态位移监测。在这些场景里，时变位移不仅是结构性能的直接反映，更是后续模态分析和安全评估的原始依据。
1.2 频响函数就是结构的指纹
频响函数，简称FRF，说白了就是输出与输入之间的频率域比值关系。对于一个线性时不变系统，如果在某一频率点以单位力激励结构，在响应点测到的响应幅值和相位，就是该频率下的频响函数值。把全频段扫描一遍，得到的一条曲线，就是这个结构的“指纹”。
这条指纹包含了结构的固有频率、阻尼比和振型信息。为什么叫指纹呢？因为结构一旦建成，它的几何尺寸、材料特性和边界约束就都定了，其固有频率和振型就基本定型了。哪怕结构出现裂缝、刚度退化，频响函数也会随之改变——峰值频率往左移、峰形变宽。健康监测、损伤识别，很多都是从频响函数的变化里找线索。
这里我需要强调一个概念，频响函数分为原点和跨点两种。原点FRF指的是激励点和响应点重合的情况，跨点FRF则对应不同位置之间的传递关系。原点FRF通常用于确定模态参数时辅助识别，而跨点FRF是振型归一化的基础。测量中两者都需要，但在传感器布置和锤击点选择上要有所侧重。
1.3 二者之间如何打通
你可能要问，时变位移和频响函数看起来是两套东西，一个在时域，一个在频域，怎么建立起联系？
答案就是傅里叶变换这座桥。时域里测得的位移时程，经过傅里叶变换得到频谱；激励力时程同样变换到频域。两者相除，就得到了频响函数。而频响函数经过逆傅里叶变换，得到脉冲响应函数——这又回到了时域。换句话说，时变位移是“原料”，频响函数是“产品”，中间的关键工艺是模态参数识别，包括峰值拾取、曲线拟合等步骤。
在实际项目中，我常用的策略是两条腿走路：先时域分析脉冲响应或自由衰减信号，初步判断结构的固有频率；再频域计算FRF，精确拟合阻尼比和振型。两条路交叉验证，结果才靠谱。单一依赖某个方法，碰到密集模态或者大阻尼结构的时候，很容易出偏差。
2. 实测方案的选型与准备：从传感器到采集仪的全链路考虑
2.1 传感器的类型选择与测点布置
测时变位移和FRF，传感器怎么选，是个真问题。
首选方案当然是位移传感器，比如拉线式位移计、电涡流位移传感器和激光位移传感器。拉线式安装方便、量程大，但有机械接触，对被测结构有一定附加质量影响；电涡流适合金属结构，非接触、精度高，但量程一般较小；激光位移传感器精度高、非接触，可是价格昂贵，现场对光条件要求也苛刻。
如果你的目标是同时获得频响函数，我更推荐以加速度计为主、位移计为辅的组合方案。加速度计频响宽、安装方便，重量轻对结构的附加质量影响小，适合做模态测试；位移计用来校核低频位移的准确性，弥补加速度积分在低频段的先天不足。
测点布置方面，有一条黄金法则：测点必须避开节点。所谓节点，就是某一阶振型中位移为零的位置。比如一个简支梁的一阶振型节点在支座处，二阶振型节点在1/4和3/4跨处。如果你把传感器恰好放在节点上，那这阶模态就完全激发不出来，FRF上对应频率根本没有峰。
那怎么提前知道节点在哪？结构越复杂，节点位置越难靠直觉判断。我的建议是：先做一次密集测点的预扫描，或者用有限元模型算出振型大致位置，再优化测点。对于大跨度结构，测点间距宜取跨度的1/20到1/10，先密后疏。
2.2 激励方式的选择：锤击、激振器还是环境激励
频响函数必须有输入和输出，这意味着你必须知道力长什么样。激励方式基本分为三类：锤击激励、激振器激励和环境激励。
锤击激励是最经济快捷的方式。一把力锤配上不同硬度的锤头，软锤头激励低频丰富，硬锤头激励高频丰富。测小型结构、机械部件，我基本都用锤击。锤击有个好处是激励信号是瞬态的，测试时间短，结构处于近似线性时不变状态，FRF质量高。但锤击能量有限，大型桥梁、高层建筑根本激发不起来。
激振器激励适合中大型结构。通过功放给激振器输入随机信号或者正弦扫频信号，把力持续地注入结构。激振器能提供稳定可控的激励信号，FRF的信噪比高，还能做多输入多输出测试，识别密集模态。但激振器的安装是个技术活——顶杆与结构之间必须紧密接触，不能有间隙，否则信号的相位就乱了，后期拟合模态参数时误差很大。
环境激励则是既不锤也不振，直接用脉动、车流、风这类自然荷载作为激励源。这种方法不需要专门的激励设备，适合大型结构运行状态下的实时监测。但环境激励下输入力无法直接测量，只能通过输出信号之间互相参考，采用频域分解法或随机子空间法来识别模态参数。严格来说这得到的是“输出-only”频响函数或相关函数，跟传统FRF不完全等价，但对大型结构而言，这是最现实、成本最低的方案。
2.3 采集系统的关键参数配置
数据采集系统是整套测试的中枢。配置采集参数时，有几个参数必须反复确认：
采样率：根据奈奎斯特定理，采样率至少要达到感兴趣最高频率的2倍，工程上通常留20%~30%的余量，即取最高分析频率的2.5倍左右。比如关心到100Hz，采样率建议256Hz以上，实际操作中我常用512Hz。
采样时长：决定了频率分辨率。频率分辨率等于采样率的倒数乘上采样点数。想要分辨两个频率相差0.1Hz的相邻模态，就需要至少10秒的记录长度。这是一笔非常容易忽略的账。
量程设置：提前估计结构响应的最大幅值，确保信号不削顶也不淹没在噪声里。一个简单的方法：先以较大量程做预采，观察信号峰值，再调整为峰值的1.5至2倍作为正式量程。
抗混叠滤波：所有通道必须开启抗混叠滤波器，否则高频成分折叠到低频段，FRF曲线会出现根本不存在的假峰。这个我在早期测试中吃过亏，后面还会细说。
另外，通道间的相位一致性必须校准。多通道采集时，如果各通道之间存在不一致的相位延迟，识别出来的振型会存在系统性偏差，尤其在跨点FRF计算时影响显著。
3. 信号处理与FRF计算实操：从原始数据到可信频响函数
3.1 时变位移数据的预处理流程
拿到位移或者加速度原始数据，第一件事不是算FRF，而是做预处理。
去趋势：时域数据存在线性趋势项，多半是传感器零点漂移或者线缆摆动引起的。用最小二乘法拟合一条直线并从原始信号中减去，这一步在做积分之前尤其关键。
滤波：根据目标频段选择带通滤波。比如关心0.5Hz到50Hz的结构响应，就用带通滤波器滤掉低频漂移和高频噪声。这里我提一个细节：数字滤波器有相位延迟，处理位移数据时要用零相位滤波，即正向滤波一次再反向滤波一次，保证滤波后的信号不产生相位失真。如果不做零相位处理，后面积分出来的位移相位乱套，FRF的相位谱就废了。
积分：加速度积分到位移，时域积分存在初值问题和漂移问题，频域积分则更稳定。常见做法是将信号变换到频域，对目标频段外的分量清零，然后对谱除以(jω)²完成二次积分，再逆变换回时域。这种做法本质上是频域带通积分，避免了时域积分的累积误差。
我实际做桥梁测试的时候，遇到过一组加速度信号，积分前看波形挺好，积分后位移居然整体飘了几十毫米——最后查下来就是因为原始信号里混入了0.1Hz附近的低频干扰。用高通滤波把0.5Hz以下清掉后再积分，位移波形才恢复正常。
3.2 FRF估计方法：H1与H2的取舍
频响函数估计不是简单地输出谱除以输入谱。实测信号里输入和输出都有噪声，直接除法得到的FRF方差很大。工程上常用H1和H2两种估计器。
H1估计假设噪声主要在输出端，估计公式为输入自功率谱除以输入输出互功率谱。H1能有效抑制输出噪声的影响，在共振峰附近估计值偏小，通常被认为是保守估计。
H2估计假设噪声主要在输入端，用输出自功率谱除以输入输出互功率谱。在反共振点附近，H2估计结果比较可靠。
实际操作中，我是这么用的：如果激励力信号干净（比如力锤信号信噪比高），优先用H1；如果响应信号很干净但力信号噪声大（比如激振器顶杆连接松动），用H2。两者都算出来，合在一起看，如果H1和H2严重不一致，说明系统存在非线性或者信号质量太差，这时FRF本身就不可信。
还要提一下相干函数。相干函数是FRF可靠性的一个晴雨表，取值0到1。相干接近1，说明输出完全由输入引起，FRF可信；相干掉到0.5以下，说明有外界噪声或者其他激励源混入，要重点排查。每次测试完，我必看相干曲线，这个习惯帮我挡掉了不少无效数据。
3.3 窗函数与平均次数设置
锤击激励测试时，力锤信号是瞬态的，响应信号是衰减的。通常对力信号加力窗，保留主脉冲附近的信息，尾部噪声清零；对响应信号加指数窗，强迫响应在窗尾衰减到零，避免截断泄漏。
这里有一个参数需要探索：指数窗的衰减系数。衰减太快，相当于人为增加了阻尼，识别出的阻尼比偏大；衰减太慢，截断泄漏没压住，谱峰附近出现涟漪。我处理钢结构时，指数窗时间常数一般取结构衰减到1%所需时间的1到1.5倍，然后通过试算对比不同衰减系数下的FRF峰值稳定性来确定。
平均次数方面，锤击测试每测点至少锤击5次取平均，激振器随机激励至少平均50次，环境激励则看数据长度，能取多少段取多少段，但每段之间要保证有重叠率（一般50%到75%），提高谱估计的平滑性。
此外，务必检查每次锤击的力幅一致性。如果第一锤500N，第二锤800N，第三锤只有200N，那平均后的FRF会被大力锤主导，相当于给数据加了权，会引入人为误差。我在测试流程中会设置触发阈值和拒绝逻辑，力峰值超过设定范围的数据点自动剔除。
4. 模态参数识别与结果验证：把FRF变成工程结论
4.1 峰值拾取法与曲线拟合法的实战对比
有了高质量的FRF，接下来要从中提取固有频率、阻尼比和振型。
峰值拾取法最简单直观：在FRF幅值曲线上找出峰值对应的频率，那就是固有频率；用半功率带宽法算阻尼比——找到峰值的0.707倍处两个频率点，它们的差除以两倍中心频率就是阻尼比。这个方法速度快，适合模态稀疏、阻尼不大的结构。但频率分辨率不足时，峰值位置和半功率带宽都不准，识别结果误差大。
曲线拟合法更精细。用一个有理分式模型去拟合实测FRF，通过优化算法求出极点、留数，进而得到模态参数。这种方法可以处理密集模态和较大阻尼的情况，精度高，但需要经验去设定拟合频段和初值。我常用的做法是先用峰值拾取法获得初始估计，再用曲线拟合法精修，人机结合效率最高。
实操中还有一种多自由度时域方法——ITD法或者STD法，直接用时域自由响应信号识别模态参数，不需要FRF，适合环境激励下的运行模态分析。这类方法的好处是不需要人工选择频段、不容易遗漏耦合模态，但对噪声敏感，需要配合奇异值分解来降秩。
4.2 模态验证的三个常规指标
模态参数识别出来，不能直接拿去写报告。三个常规验证指标必须检查：
稳态图：这是模态分析软件里最常用的一眼验证工具。在不同模型阶数下，频率和阻尼估计值稳定的点连成一条竖线，说明这个模态是真实的；到处飘忽不定的点，多半是噪声造成的虚假模态。
MAC矩阵：模态置信准则，用来验证不同模态之间振型的正交性。两个不同模态的MAC值应接近0，对角线上的MAC值接近1。如果两个模态MAC很高，说明它们其实是同一个模态，或者是测试点布置不够导致振型无法区分。
模态参与因子：反映各阶模态在实测响应中的贡献大小。某阶模态参与因子过低，说明激励能量没分布到这一阶上，识别的参数不可靠。要看激励点在振型上的位置——如果激励点恰好在某阶振型的节点附近，这阶模态就激不出来了。
我在一个斜拉桥模型的测试里，就遇到过一阶对称竖弯和正对称侧弯频率靠得很近的情况。峰值拾取法只看幅值谱，两者峰完全重叠，误判成一个模态。后来用曲线拟合法加上MAC矩阵验证，才把它们分开。这种经验告诉我：模态参数识别永远要留一手验证手段。
4.3 时变位移与FRF结合评估结构状态
做模态参数识别不是为了发论文，而是为了判断结构状态。一个完整的评估流程应该是：
第一，对比实测FRF与初始健康状态下的FRF。如果固有频率下降超过5%，多半意味着刚度退化，典型原因是裂缝开展或支座松动。如果阻尼比明显上升，说明结构内部摩擦增大，可能伴随损伤。
第二，用模态振型结合时变位移判断损伤位置。某阶振型在某个测点出现明显趋势异常，说明该区域存在局部变形集中。比如梁体某段实测位移比理论预测大，且该段FRF幅值谱出现新峰，就要重点检查该区段是否存在截面削弱。
第三，用FRF曲线的重复性测试评估结构时变性。同一个测点，在不同温度、不同荷载状态下测FRF，观察峰频和峰型的漂移规律。如果是混凝土结构，温度引起的频率变化范围通常在1%到2%以内，超出这个范围就要警惕材料劣化。
实测中要注意的是，时变位移和FRF识别结果只能作为评估依据之一，绝对不能替代定期外观检查和材料检测。它们告诉你“结构状态可能变了”，但要回答“哪里变了、为什么变”，还需要结合检测数据综合判断。
5. 现场常见问题与排查技巧实录
5.1 数据漂移与低频失真问题
这是时变位移测试的第一号杀手。现象是位移时程整体线性爬升或缓慢波动，甚至出现完全不合理的单向偏移。排查思路按顺序执行：
检查传感器是否固定牢靠，尤其是磁吸座和胶粘方式。磁吸座在钢结构上看似结实，受到轻微冲击时会有微观滑移，导致基线偏移。检查线缆是否存在摆动，线缆在风中摆动产生的运动伪迹，在低频段体现特别明显。检查采集仪的耦合方式，如果用了交流耦合而非直流耦合，低频成分会被滤掉，位移时程峰值会严重衰减。
正常情况下，温度变化引起的结构位移速率应该与气温变化速率同步。如果位移信号中含有明显的与结构振动频率无关的缓慢分量，就要考虑是不是地脉动或者附近施工振动引入的低频干扰。
5.2 频响函数曲线毛刺过多
FRF曲线像刺猬一样到处是毛刺，频谱上到处是假峰，这种情况多是信噪比不足导致的。
先看力信号的质量。锤击测试中，力信号如果触发阈值太低，力通道记录到大量无效的微小锤击，计算FRF时噪声占了主导。对策是提高触发阈值，或者把力窗变窄，只保留主脉冲。
再看响应信号。如果响应的信号幅值只有采集量程的1%以下，量化噪声就非常严重，FRF上每个频点都带着大量噪声。对策是调整量程，让响应幅值占量程的30%到80%。
还有一个常见原因是锤击点位置漂移。手动锤击时如果每次敲在结构上的位置略有偏差，激励点与响应点之间的传递路径每次都不一样，平均出来的FRF自然毛刺丛生。对策是加装锤击定位导向套，或者用自动锤击装置。
5.3 模态丢失与振型错乱的原因排查
结构明明存在某阶模态，FRF上却没有对应峰值，或者振型形状完全不对，这通常有三种原因：
激励点在该阶模态的节点附近。就像前面说的，节点处振幅为零，你在这里输入能量，该阶模态根本起不来。对策是通过有限元预分析或预扫描判断节点位置，把激励点挪到远离节点的位置。
传感器安装位置刚度不足。传感器安装底座用胶水粘在一个薄板上，薄板本身局部振动模态可能比结构模态更强，实测FRF上全是局部模态，结构模态反倒被淹没。对策是增加底座刚度，确保安装面与被测结构融为一体。
响应信号相位翻转。这种情况常见于传感器安装方向不一致。加速度计有方向标记，如果装反了，振型符号全反，MAC矩阵对角线元素甚至会变成负数。对策是安装前逐通道校准方向，并在测试记录表上标注每个测点的正方向。
5.4 测试数据的快速自检清单
结合多次项目经验，我给自己定了一套快速自检清单，每次测试收工前逐项打勾：
所有通道的时域波形是否存在明显的削顶或饱和。所有通道的频谱在关心频段内是否有明显的谱峰，谱峰数量是否与结构阶数大致对应。相干函数在主要峰值频率附近是否大于0.9。FRF的峰值频率与有限元预估值的差异是否在合理范围（通常10%以内）。锤击测试中，不同次锤击的FRF峰值变化是否在3%以内。如果以上任何一项不满足，宁可重新测试，也不要拿有问题的数据进入分析流程。
这套清单看起来简单，但真金白银换来的教训是：现场多花半小时做自检，省下的是后面数据分析阶段三天的返工时间。早期我试过直接把数据拷走，回办公室发现半批数据不能用，再跑一趟现场补测，费力费钱还丢面子。
6. 个人体会与后续扩展建议
做结构响应测试与频响函数分析这么多年，我最大的感触是：高级的分析方法固然重要，但扎实的测试基本功才是决定成败的关键。一套配置合理的测试方案，一次严格的数据自检，比任何花哨的识别算法都更能保障结论的可靠性。很多初学者一上来就想用复杂的时域算法识别模态，结果输入数据一塌糊涂，再先进的算法也只是在垃圾数据上精雕细琢。
对于想在这个方向深入的朋友，我建议按这样的路径扩展：
把单点激励扩展为多点激励。多参考点锤击或者多激振器激励，能有效激发密集模态和扭转模态，识别结果更全面。把稳态测试扩展为跟踪测试。在结构上长期布置传感器，定时采集环境激励数据，用频响函数峰值频率的变化趋势做长期健康监测。把线性假设扩展到非线性识别。当结构出现明显非线性行为时，比如裂缝开闭、松动接触，常规FRF已经失效，需要引入希尔伯特变换或非线性模态分析的方法。把单学科扩展到交叉融合。结构响应数据结合有限元模型修正，构建数字孪生体，实时评估结构的安全状态。
最后分享一个小技巧：对大型结构做环境激励测试时，不要只看单一测点的自功率谱，要把多个测点的互功率谱相位信息结合起来看。互功率谱的相位在结构模态频率处会发生规律性翻转，比自功率谱的峰值更稳定、更容易识别真实模态，尤其在模态密集的情况下效果显著。这个小技巧帮我处理过好几个棘手的密集模态项目，希望对你也管用。
