CN107727345A - 一种用于加速试验的非平稳振动信号生成方法 - Google Patents
一种用于加速试验的非平稳振动信号生成方法 Download PDFInfo
- Publication number
- CN107727345A CN107727345A CN201710844608.6A CN201710844608A CN107727345A CN 107727345 A CN107727345 A CN 107727345A CN 201710844608 A CN201710844608 A CN 201710844608A CN 107727345 A CN107727345 A CN 107727345A
- Authority
- CN
- China
- Prior art keywords
- vibration
- root
- grms
- mean
- window
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M7/00—Vibration-testing of structures; Shock-testing of structures
- G01M7/02—Vibration-testing by means of a shake table
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明涉及一种用于加速试验的非平稳振动信号生成方法。本发明包括步骤:原始振动加速度信号分类、分解成高斯振动信号和正弦振动信号、合成新的非高斯振动信号。本发明还原了道路运输过程中随机振动的非高斯振动特性及本质存在的加速度瞬变。本发明在原始振动信号分解过程中将加速度瞬变信号单独分离出来,使得重新生成的加速度信号可用于振动加速试验,本发明保证了测试振动时间压缩过程中,瞬变加速度信号的时长以及振动量级不受改变。
Description
技术领域
本发明属于振动试验领域,具体是涉及一种用于加速试验的非平稳振动信号生成方法。
背景技术
振动试验是指评定产品在预期的使用环境中抗振能力而对受振动的实物或模型进行的试验。根据施加的振动载荷的类型把振动试验分为正弦振动试验和随机振动试验两种。现有模拟随机振动的方法:将采集到的实际振动信号时域谱(加速度谱,如图1所示)通过傅里叶变换得到对应的频域谱(功率密度谱,如图2所示),即通过频域谱信号控制振动试验过程,再将所述频域谱通过逆傅立叶变换还原得到对应的时域谱(如图3所示)进行输出。该方法虽然表现出了运输振动过程的随机性,且试验操作简单,只需要对单一参数进行控制。但是该方法存在的缺点在于:输出的随机振动信号为高斯随机信号(峭度=3),而实际运输振动信号为非高斯随机信号(峭度>3);输出的是稳定振动信号,忽略了本质存在的大量的加速度瞬变,这些瞬变对产品的危害程度远大于稳态振动部分,且由于瞬变加速度的值普遍较大,其对产品的损伤机理可能包括跌落、碰撞等,因此对瞬变振动信号的忽略会致使测试结果不准确;对信号非高斯、非平稳性的忽略所造成的偏差,在加速试验下被进一步扩大。
发明内容
为了解决上述技术问题,本发明提供一种用于加速试验的非平稳振动信号生成方法。
为了实现本发明的目的,本发明采用了以下技术方案:
一种用于加速试验的非平稳振动信号生成方法,包括以下步骤:
S1、对实际振动加速度信号进行加窗处理,求得每个窗口内加速度信号均方根值Grms,其中其中N为每个窗口内的样本数量,gi为窗口内第i个加速度信号样本值;
S2、对每一个窗口的均方根值Grms进行判断并分类存储:
若该窗口的均方根值Grms大于前、后窗口均方根值之和的L倍,L>1,则认为该窗口内振动加速度信号发生瞬变并将该窗口的均方根值Grms存入瞬变数据集合中,同时将该窗口的均方根值Grms缩小至前、后窗口均方根值的平均值后存入稳态数据集合中;
若该窗口的均方根值Grms小于或等于前、后窗口均方根值之和的L倍,则认为该窗口内振动加速度信号平稳并将该窗口的均方根值Grms存入稳态数据集合中;
S3、求得所述稳态数据集合中每个均方根值Grmsjm对应的概率密度函数,所述概率密度函数其中Tjm为相应均方根值Grmsjm对应窗口的振动时长,Tj为实际振动总时长;
S4、根据公式将所述概率密度函数Pm对应的振动时长Tjm进行压缩,所述Tpm表示压缩后的测试振动时长,Grmspm表示压缩后的测试振动加速度信号均方根值,k为常数取2~5;所述其中Tp表示压缩后的测试振动总时长;
S5、以所述测试振动加速度信号均方根值Grmspm为振幅,以Tpm为测试振动时长,生成X段连续的高斯振动信号,其中X表示振动量级数且与所述均方根值Grmsjm的个数相等;
另外,将所述瞬变数据集合中每个均方根值Grmsjn对应窗口内的加速度信号转化为正弦振动信号,所述正弦振动信号gt=GrmsjnSin2πωt,其中ω为加速度信号采样频率,t∈﹝0,a﹞,a为所述窗口时长;
S6、将所述正弦振动信号与所述高斯振动信号进行叠加形成用于控制振动试验的非高斯振动信号,所述正弦振动信号与高斯振动信号叠加部分对应的均方根值Grmsjm与Grmsjn相等。
进一步,所述步骤S1中N取样本总数量的0.5~1%。
进一步,所述步骤S2中L=1.5。
进一步,所述步骤S4中Grmspm的值小于1G。
本发明的有益效果在于:
本发明还原了道路运输过程中随机振动的非高斯振动特性及本质存在的加速度瞬变。本发明在原始振动信号分解过程中将加速度瞬变信号单独分离出来,使得重新生成的加速度信号可用于振动加速试验,本发明保证了测试振动时间压缩过程中,瞬变加速度信号的时长以及振动量级不受改变。
附图说明
图1为一种原始振动信号时域谱。
图2为由图1进行傅里叶变换得到的振动信号频域谱。
图3为由2进行逆傅里叶变换得到的振动信号时域谱。
图4为一种实际振动加速度信号进行加窗处理后的结果。
图5为均方根值Grms进行判断并分类存储的流程图。
图6为一种压缩前每个均方根值Grmsjm对应的概率密度函数条形图。
图7为压缩后每个均方根值Grmspm对应的概率密度函数条形图。
图8为一种新生成的高斯振动信号波形图。
图9为一种新生成的正弦振动信号波形图。
图10为一种新生成的非高斯振动信号波形图。
图11为测试振动功率密度谱与实际振动功率密度谱比较示意图。
图12为本发明方法与传统方法分别生成的振动信号波形比较示意图。
具体实施方式
下面结合实施例对本发明技术方案做出更为具体的说明:
S1、对实际振动加速度信号进行加窗处理,求得每个窗口内加速度信号均方根值Grms,其中其中N为每个窗口内的样本数量,gi为窗口内第i个加速度信号样本值;窗口内样本数量N太少会增加计算量及数据数量,窗口内样本数量N太多可能会忽略掉存在的一些瞬变信息,窗口内样本数量N可根据采样数据量及分析需求的精度决定,一般情况下,N可取样本总数量的0.5~1%。如图4所示即为一种实际振动加速度信号进行加窗处理后的结果。
S2、对每一个窗口的均方根值Grms进行判断并分类存储,具体流程如图5所示:
若该窗口的均方根值Grms大于前、后窗口均方根值之和的L倍,L>1,则认为该窗口内振动加速度信号发生瞬变并将该窗口的均方根值Grms存入瞬变数据集合中,同时将该窗口的均方根值Grms缩小至前、后窗口均方根值的平均值后存入稳态数据集合中;
若该窗口的均方根值Grms小于或等于前、后窗口均方根值之和的L倍,则认为该窗口内振动加速度信号平稳并将该窗口的均方根值Grms存入稳态数据集合中;
S3、求得所述稳态数据集合中每个均方根值Grmsjm对应的概率密度函数,所述概率密度函数其中Tjm为相应均方根值Grmsjm对应窗口的振动时长,Tj为实际振动总时长;如图6所示即为每个均方根值Grmsjm对应的概率密度函数条形图。
S4、根据公式将所述概率密度函数Pm对应的振动时长Tjm进行压缩,所述Tpm表示压缩后的测试振动时长,Grmspm表示压缩后的测试振动加速度信号均方根值,k为常数取2~5,k值取决于包装或产品的材料、结构等;所述其中Tp表示压缩后的测试振动总时长;如图7所示即为压缩后得到的测试振动加速度信号均方根值Grmspm对应的概率密度函数条形图。
通过对每个均方根值Grmsjm对应窗口的振动时长进行压缩,从而实现对整体测试时间的压缩。经过时间压缩后的Grmspm与原Grmsjm相比具有不同的值,同时,为了确保不产生超过1G的测试振动信号,在实际操作中,应根据振动加速度信号的量级及压缩倍率合理选择压缩过程。
例如某实际采样的加速度信号,振动时长为165分钟(约10000秒),分解后其各均方根值Grmsjm对应的概率密度分布如下表所示:
表1时间压缩前的加速度均方根值的概率密度分布:
为缩短测试时间,提高试验效率,现计划采用55分钟(约3300秒)振动试验,即压缩测试时间为处理前的1/3,k取2,则
以此类推,得出3倍时间压缩后的加速度均方根值概率密度分布,如下表所示:
表2 3倍时间压缩后的加速度均方根值的概率密度分布
由于瞬变数据集合中相应的加速度均方根值代表运输过程中本质存在的由客观因素(路面崎岖、急转弯等)造成的加速度瞬变,这些瞬变对产品的危害程度远大于稳态振动部分,且由于瞬变加速度的值普遍较大,其对产品的损伤机理可能包括跌落、碰撞等。因此,在加速试验中,不应压缩该部分内容,而需要保证本质存在的加速度瞬变可以完整的呈现在振动测试中。
S5、以所述测试振动加速度信号均方根值Grmspm为振幅,以Tpm为测试振动时长,生成X段连续的高斯振动信号,其中X表示振动量级数且与所述均方根值Grmsjm的个数相等;如图8所示即为一种新生成的高斯振动信号波形图。
另外,将所述瞬变数据集合中每个均方根值Grmsjn对应窗口内的加速度信号转化为正弦振动信号,所述正弦振动信号gt=GrmsjnSin2πωt,ω为信号采样频率,t∈﹝0,a﹞,a为所述窗口时长;如图9所示即为一种新生成的正弦振动信号波形图。
S6、将所述正弦振动信号与所述高斯振动信号进行叠加并形成用于控制振动试验的非高斯振动信号,所述正弦振动信号与高斯振动信号叠加的部分对应的均方根值Grmsjm与Grmsjn相等,如图10所示即为新生成的非平稳非高斯振动信号波形图。
将本发明方法所生成的非平稳非高斯振动信号用于控制振动试验,通过对加速度传感器测得的振动试验生成的加速度信号进行傅里叶变化获得对应的振动测试功率密度谱,将所述测试振动功率密度谱与实际振动功率密度谱进行比较吻合度好,如图11所示。由此说明本发明所生成的非高斯振动信号适用于加速度振动试验,且试验结果准确,可靠性高。
Claims (4)
1.一种用于加速试验的非平稳振动信号生成方法,其特征在于,包括以下步骤:
S1、对实际振动加速度信号进行加窗处理,求得每个窗口内加速度信号均方根值Grms,其中其中N为每个窗口内的样本数量,gi为窗口内第i个加速度信号样本值;
S2、对每一个窗口的均方根值Grms进行判断并分类存储:
若该窗口的均方根值Grms大于前、后窗口均方根值之和的L倍,L>1,则认为该窗口内振动加速度信号发生瞬变并将该窗口的均方根值Grms存入瞬变数据集合中,同时将该窗口的均方根值Grms缩小至前、后窗口均方根值的平均值后存入稳态数据集合中;
若该窗口的均方根值Grms小于或等于前、后窗口均方根值之和的L倍,则认为该窗口内振动加速度信号平稳并将该窗口的均方根值Grms存入稳态数据集合中;
S3、求得所述稳态数据集合中每个均方根值Grmsjm对应的概率密度函数,所述概率密度函数其中Tjm为相应均方根值Grmsjm对应窗口的振动时长,Tj为实际振动总时长;
S4、根据公式将所述概率密度函数Pm对应的振动时长Tjm进行压缩,所述Tpm表示压缩后的测试振动时长,Grmspm表示压缩后的测试振动加速度信号均方根值,k为常数取2~5;所述其中Tp表示压缩后的测试振动总时长;
S5、以所述测试振动加速度信号均方根值Grmspm为振幅,以Tpm为测试振动时长,生成X段连续的高斯振动信号,其中X表示振动量级数且与所述均方根值Grmsjm的个数相等;
另外,将所述瞬变数据集合中每个均方根值Grmsjn对应窗口内的加速度信号转化为正弦振动信号,所述正弦振动信号gt=GrmsjnSin2πωt,其中ω为加速度信号采样频率,t∈﹝0,a﹞,a为所述窗口时长;
S6、将所述正弦振动信号与所述高斯振动信号进行叠加形成用于控制振动试验的非高斯振动信号,所述正弦振动信号与高斯振动信号叠加部分对应的均方根值Grmsjm与Grmsjn相等。
2.如权利要求1所述的用于加速试验的非平稳振动信号生成方法,其特征在于:所述步骤S1中N取样本总数量的0.5~1%。
3.如权利要求1所述的用于加速试验的非平稳振动信号生成方法,其特征在于:所述步骤S2中L=1.5。
4.如权利要求1所述的用于加速试验的非平稳振动信号生成方法,其特征在于:所述步骤S4中Grmspm的值小于1G。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710844608.6A CN107727345B (zh) | 2017-09-19 | 2017-09-19 | 一种用于加速试验的非平稳振动信号生成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710844608.6A CN107727345B (zh) | 2017-09-19 | 2017-09-19 | 一种用于加速试验的非平稳振动信号生成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107727345A true CN107727345A (zh) | 2018-02-23 |
CN107727345B CN107727345B (zh) | 2019-08-16 |
Family
ID=61206622
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710844608.6A Active CN107727345B (zh) | 2017-09-19 | 2017-09-19 | 一种用于加速试验的非平稳振动信号生成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107727345B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109490626A (zh) * | 2018-12-03 | 2019-03-19 | 中车青岛四方机车车辆股份有限公司 | 一种基于非平稳随机振动信号的标准psd获取方法及装置 |
TWI685758B (zh) * | 2018-12-12 | 2020-02-21 | 周家蓓 | 道路平整度精進加速度均方根指標演算法及其系統 |
CN111104374A (zh) * | 2019-12-25 | 2020-05-05 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种实测振动数据处理方法及系统 |
CN117191311A (zh) * | 2023-08-14 | 2023-12-08 | 暨南大学 | 物流非平稳、非高斯振动下产品的加速振动试验方法 |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111458147B (zh) * | 2020-04-09 | 2021-04-20 | 北京化工大学 | 一种基于加速度传感器的滚动轴承状态评估方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5767406A (en) * | 1996-09-30 | 1998-06-16 | Ford Motor Company | Method to specify random vibration tests for product durability validation |
JP2012233814A (ja) * | 2011-05-06 | 2012-11-29 | Microsignal Kk | 振動試験制御装置 |
CN104729677A (zh) * | 2014-12-31 | 2015-06-24 | 清华大学 | 一种非平稳噪声信号的时域数字计权方法 |
CN106248408A (zh) * | 2016-10-08 | 2016-12-21 | 中车株洲电力机车研究所有限公司 | 一种动车组变流器的振动环境长寿命试验谱的生成方法 |
-
2017
- 2017-09-19 CN CN201710844608.6A patent/CN107727345B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5767406A (en) * | 1996-09-30 | 1998-06-16 | Ford Motor Company | Method to specify random vibration tests for product durability validation |
JP2012233814A (ja) * | 2011-05-06 | 2012-11-29 | Microsignal Kk | 振動試験制御装置 |
CN104729677A (zh) * | 2014-12-31 | 2015-06-24 | 清华大学 | 一种非平稳噪声信号的时域数字计权方法 |
CN106248408A (zh) * | 2016-10-08 | 2016-12-21 | 中车株洲电力机车研究所有限公司 | 一种动车组变流器的振动环境长寿命试验谱的生成方法 |
Non-Patent Citations (2)
Title |
---|
丁康 等: "平稳和非平衡振动信号的若干处理方法及发展", 《振动工程学报》 * |
李舜酩 等: "振动信号处理方法综述", 《仪器仪表学报》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109490626A (zh) * | 2018-12-03 | 2019-03-19 | 中车青岛四方机车车辆股份有限公司 | 一种基于非平稳随机振动信号的标准psd获取方法及装置 |
CN109490626B (zh) * | 2018-12-03 | 2021-02-02 | 中车青岛四方机车车辆股份有限公司 | 一种基于非平稳随机振动信号的标准psd获取方法及装置 |
TWI685758B (zh) * | 2018-12-12 | 2020-02-21 | 周家蓓 | 道路平整度精進加速度均方根指標演算法及其系統 |
CN111104374A (zh) * | 2019-12-25 | 2020-05-05 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种实测振动数据处理方法及系统 |
CN111104374B (zh) * | 2019-12-25 | 2023-04-07 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种实测振动数据处理方法及系统 |
CN117191311A (zh) * | 2023-08-14 | 2023-12-08 | 暨南大学 | 物流非平稳、非高斯振动下产品的加速振动试验方法 |
CN117191311B (zh) * | 2023-08-14 | 2024-05-24 | 暨南大学 | 物流非平稳、非高斯振动下产品的加速振动试验方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107727345B (zh) | 2019-08-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107727345A (zh) | 一种用于加速试验的非平稳振动信号生成方法 | |
CN106840637B (zh) | 基于改进hht算法的gis机械振动信号时频分析方法 | |
CN111106783B (zh) | 一种信号制作方法、信号制作装置、振动马达及触屏设备 | |
CN103217213B (zh) | 基于响应信号时频联合分布特征的模态参数辨识方法 | |
KR900008438A (ko) | 신호처리 방법 및 음향소스 데이타 형성장치 | |
CN104462695A (zh) | 双耦合Duffing振子与变尺度相结合的微弱信号检测方法 | |
Zhou et al. | Finite element explicit dynamics simulation of motion and shedding of jujube fruits under forced vibration | |
CN106768757A (zh) | 基于变采样率采样的振动台超低频正弦扫频信号幅值识别方法 | |
CN103776901B (zh) | 基于振动响应信息的粘弹夹层结构老化状态识别方法 | |
Lancieri et al. | Spectral matching in time domain: A seismological and engineering analysis | |
CN109711355A (zh) | 一种物流包装冲击信号识别与提取方法 | |
Bass et al. | Modification of sonic boom wave forms during propagation from the source to the ground | |
Zhang et al. | Analysis of pear dynamic characteristics based on harmonic response | |
CN104568137B (zh) | 噪声测试中特定噪声源贡献率的判定方法 | |
CN105806574A (zh) | 一种基于动应力测试的振动加速模型的确定方法 | |
CN107544005A (zh) | 一种高压igbt局部放电电流时域参数的确定方法和装置 | |
CN105759156A (zh) | 一种多方向振动能量收集装置性能测试系统及方法 | |
Song et al. | Decomposition of seismic signal based on Hilbert-Huang transform | |
CN107786204A (zh) | 数模转换器参数测试系统及方法 | |
Ai et al. | Gear fault detection based on ensemble empirical mode decomposition and Hilbert-Huang transform | |
Lai et al. | The study and implementation of signal processing algorithm for digital beam position monitor | |
Sotirios et al. | Input stimulus comparison using an adaptive FPGA-based testing system | |
WO2022000656A1 (zh) | 马达系统失真的测量方法及设备、计算机可读存储介质 | |
CN106954169B (zh) | 一种声学产品的噪音测试方法及其系统 | |
Akinci et al. | Jitter Modeling for High Precision Frequency Measurements in Oscillator Circuits |
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 | ||
CB02 | Change of applicant information | ||
CB02 | Change of applicant information |
Address after: 230031 No. 888 Changjiang West Road, Shushan District, Hefei, Anhui. Applicant after: Hefei General Machinery Research Institute Co., Ltd. Address before: 230031 No. 888 Changjiang West Road, Shushan District, Hefei, Anhui. Applicant before: Hefei Universal Machinery Inst. |
|
GR01 | Patent grant | ||
GR01 | Patent grant |