实测桥梁加速度数字积分得速度,位移,提供原始数据
本帖最后由 westrongmc 于 2012-12-28 00:40 编辑对于南京某桥梁，现场实测的加速度信号进行一次数字积分后得到速度信号，结果如下：
由上自下依次为：
加速度信号--积分得到的速度信号--前二者重叠显示。
拉开看局部细节：
进行二次数字积分后得到位移信号，结果如下：
由上自下依次为：
原始加速度信号（蓝）与 二次积分后的位移信号（红）重叠显示；拉开看其局部细节。
处理的基本思路：
利用FIR数字算法进行一次积分，二次积分处理。
其中二次积分时，对信号进行了0.3Hz的高通滤波处理。
网上搜了一张桥梁照片供参考。
如大家有兴趣，可下载附件中的原始加速度信号进行处理。
附件数据格式为ASCII码，采样频率Fs=10.245464Hz，共12288个样本点，单位：mm/s^2。
不妨把你的结果展示一下，让大家看看效果。
我觉得你的加速度二次积分后的位移信号幅值好像不对。 本帖最后由 impulse 于 2012-12-28 10:57 编辑
分别采用NI的声学振动工具包和自己编制的C++算法，结果如下：
1、一次积分
2、二次积分
3、加速度、速度、位移频谱图
能否把积分结果与原始数据重叠显示的波形拉开？看看细节。我还关心相位信息，呵呵。 本帖最后由 impulse 于 2012-12-28 11:09 编辑
westrongmc 发表于 2012-12-28 10:36 static/image/common/back.gif
能否把积分结果与原始数据重叠显示的波形拉开？看看细节。我还关心相位信息，呵呵。
NI没有采用零相位滤波，它是用来做实时分析的，相位肯定有点误差。以下结果位移放大了十倍。
1、加速度、速度信号局部放大比较
2、加速度、位移信号局部放大比较
这是南京大胜关长江大桥啊，不知道测点在该桥的什么位置。 zswseu 发表于 2012-12-28 12:01 static/image/common/back.gif
这是南京大胜关长江大桥啊，不知道测点在该桥的什么位置。
1/4跨位置。 从相位上来看，楼主的处理效果很好，不知道断点效果如何，能否看一下0时刻的对比结果。 本帖最后由 westrongmc 于 2013-1-1 23:35 编辑
前50秒数据，如下图所示。
后50秒数据，如下图所示。
自己写的C++程序积分得的位移结果跟硬件积分得的位移结果一样吗？看纵坐标的值貌似有点不同 impulse 发表于 2012-12-28 09:37 static/image/common/back.gif
分别采用NI的声学振动工具包和自己编制的C++算法，结果如下：
1、一次积分
从频谱上看好像做滤波了，是吗？ liufengsjzri 发表于 2013-1-6 22:50 static/image/common/back.gif
从频谱上看好像做滤波了，是吗？
积分的过程需要有高通滤波 {:{39}:} 请问高通滤波是怎么定的频率？0.5Hz?我用PCB的加速度计， 用NI的声学工具SVL Integration.vi测量， 按加速度计的要求设置了0.5Hz的高通滤波，在没有振动的时候位移积分有明显的底噪，位移积分最大差不多有0.1mm. 底噪感觉还是周期波， 频率差不多2Hz。 如果把高通滤波设高， 5Hz就基本上没有了， 不知道你在测试没有振动的时候， 底噪位移有多大？谢谢