CN110186632B - 一种基于三角形下凹的基础激励振动条件的振动试验方法 - Google Patents
一种基于三角形下凹的基础激励振动条件的振动试验方法 Download PDFInfo
- Publication number
- CN110186632B CN110186632B CN201910428714.5A CN201910428714A CN110186632B CN 110186632 B CN110186632 B CN 110186632B CN 201910428714 A CN201910428714 A CN 201910428714A CN 110186632 B CN110186632 B CN 110186632B
- Authority
- CN
- China
- Prior art keywords
- beta
- vibration
- points
- excitation
- normalized
- 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
Images
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
- G01M7/025—Measuring arrangements
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明公开了一种基于三角下凹基础激励振动条件的振动试验方法,利用单自由度系统传递函数可快速求得响应限制对应的频点β1、β2,利用频点β1、β2及β0构造出二次曲线,进而构造出基础激励的三角形下凹方案,在本发明的三角形下凹方法下,产品响应在响应限制附近有微小浮动,既可保证产品响应不会太大以致造成产品破坏,又可保证产品响应不会太小以致造成“欠试验”的情况。
Description
技术领域
本发明属于结构动力学及振动试验技术领域,具体涉及一种基于三角形下凹基础激励振动条件的振动试验方法。
背景技术
工业产品及航天产品(尤其是军品)在产品设计及验证时通常都需要做振动试验,用于考核产品是否能够承受生命周期内的振动载荷。典型的产品做振动试验时的状态如附图1所示:M0为振动台,m为产品;X0为M0的位移,即基础激励振动条件;X为m的位移,即产品的响应。
在正弦振动试验或随机振动试验中,在产品共振区附近,极易出现产品的响应X大于响应限制Xlim的情况。当产品的响应大于响应限制时,产品会出现破坏。如附图2所示,在频段β1~β2内(β为归一化频率,β=ω0/ω,其中ω0为激励X0的圆频率,ω为产品的圆频率),产品的响应超出响应限制Xlim (响应限制Xlim是产品的固有属性,产品响应超过此数值时,产品易发生破坏)。为防止产品发生破坏,需要修改基础激励振动条件X0,即对基础激励振动条件进行下凹。
上述下凹方法为矩形下凹。下凹之后的产品响应如附图3所示。由附图3 可知,在β=β1、β≈1、β=β2处,产品最大响应等于响应限制Xlim。但是在频段β1~1和频段1~β2内,产品响应较小,远低于响应限制Xlim,这必然导致在频段β1~1和频段1~β2内,产品性能没有被充分考核,振动试验存在“欠试验”的风险。
发明内容
有鉴于此,本发明的目的是提供一种基于三角形下凹基础激励振动条件的振动试验方法,可以保证产品在振动试验中的最大响应不大于响应限制Xlim,又能够保证在下凹频段内不会出现因条件下凹太大而出现产品响应太小以致“欠试验”的情况。
一种振动试验方法,包括如下步骤:
步骤1、建立正弦振动系统传递函数:
其中,X0为振动台M0的位移,X为试验产品m的位移;j为复数单位,β为归一化频率;ζ为试验产品的阻尼比;
步骤2、设正弦振动试验中位移X的最大值为Xmax,响应限制为Xlim;
步骤6、在基础激励条件坐标系中,过(β1,1)和(β0,a)两点并以点(β0,a)作为极值点,构造一条二次曲线;
步骤7、令β3=β1+λsin×(β0-β1),求出步骤6构造的二次曲线在β3处的值 x3;
步骤9、过(β2,1)和(β0,a)两点并以点(β0,a)作为极值点,构造一条二次曲线;定义β4=β2-λsin×(β2-β0),求出本步骤构造的二次曲线在β4 处的值x4;其中λsin=0.6;
步骤12、采用步骤11获得的基础振动激励X0NEW进行振动试验。
一种振动试验方法,包括如下步骤:
步骤1、建立随机振动系统传递函数:
其中,j为复数单位,β为归一化频率;ζ为试验产品的阻尼比;
步骤2、设正弦振动试验中试验产品的位移X的最大值为Xmax,响应限制为 Xlim;
步骤6、在基础激励条件坐标系中,过(β1,1)和(β0,a)两点并以点 (β0,a)作为极值点,构造一条二次曲线;
步骤7、令β3=β1+λrand×(β0-β1),求出步骤6构造的二次曲线在β3处的值 x3;λrand=0.5;
步骤9、过(β2,1)和(β0,a)两点并以点(β0,a)作为极值点,构造一条二次曲线;定义β4=β2-λrand×(β2-β0),求出本步骤构造的二次曲线在β4 处的值x4;
步骤12、采用步骤11获得的基础振动激励X0NEW进行振动试验。
本发明具有如下有益效果:
本发明提出的一种基于三角下凹基础激励振动条件的振动试验方法,利用单自由度系统传递函数可快速求得响应限制对应的频点β1、β2,利用频点β1、β2及β0构造出二次曲线,进而构造出基础激励的三角形下凹方案,在本发明的三角形下凹方法下,产品响应在响应限制附近有微小浮动,既可保证产品响应不会太大以致造成产品破坏,又可保证产品响应不会太小以致造成“欠试验”的情况。
附图说明
图1为典型的振动试验状态图;
图2为振动试验中产品响应超出响应限制示意图;
图3为传统的矩形下凹之后产品响应示意图;
图4为无量纲归一化响应与下凹参数示意图;
图5为无量纲归一化基础激励的三角形下凹示意图;
图6为正弦振动中,传统矩形下凹与三角形下凹得到的响应对比示意图。其中下凹参数a=0.6、产品阻尼比ζ=0.05;
图7为正弦振动中,传统矩形下凹与三角形下凹得到的响应对比示意图。其中下凹参数a从0.3到0.95变化,产品阻尼比ζ=0.05;
图8为正弦振动中,传统矩形下凹与三角形下凹得到的响应对比示意图。其中下凹参数a=0.6,产品阻尼比ζ从0.01到0.1变化;
图9为随机振动中,传统矩形下凹与三角形下凹得到的响应对比示意图。其中下凹参数a=0.6、产品阻尼比ζ=0.05;
图10为随机振动中,传统矩形下凹与三角形下凹得到的响应对比示意图。其中下凹参数a从0.3到0.95变化,产品阻尼比ζ=0.05;
图11为随机振动中,传统矩形下凹与三角形下凹得到的响应对比示意图。其中下凹参数a=0.6,产品阻尼比ζ从0.01到0.1变化。
具体实施方式
下面结合附图并举实施例,对本发明进行详细描述。
步骤1、建立系统传递函数。不失一般性,可将试验产品简化为附图1所示的单自由度振动系统,M0为振动台,m为试验产品;X0为M0的位移,X为m的位移;系统的各个参数可在试验之前进行预估。
在正弦振动中,系统传递函数H(ω)由X0和X的幅值比值确定:有如下关系:
其中,j为复数单位,β为归一化频率,β=ω0/ω,其中ω0为激励X0的圆频率,ω为试验产品的圆频率;ζ为试验产品的阻尼比;
在随机振动中,系统传递函数HPSD(ω)由X0和X的功率谱密度(PSD谱)比值表示,有如下关系:
上述公式(1)和公式(2)为结构动力学领域中的基本公式,在此不做过多解释。
步骤2、设正弦振动试验中X的最大值为Xmax,(结构动力学一般知识, Xmax≈X0/2ζ)响应限制为Xlim;
步骤3、定义下凹参数其中,a的值根据实验中响应限制Xlim确定;定义X的无量纲归一化响应则的最大值为1,取最大值时对应的频点为β0;其中,由结构动力学一般知识可知,β0实际值比1略小,跟系统的阻尼比ζ有关,即β0≈1;设时,对应的频点为β1、β2,如附图4所示(本图中以下凹参数a的值等于0.6为例,下同);
附表1:正弦振动参数
步骤7、令β3=β1+λsin×(β0-β1),求出步骤6构造的二次曲线在β3处的值 x3,其中λsin=0.6;附表1中已列出不同ζ、不同a对应的β3、x3的值,可直接使用;
步骤9、同理,定义β4=β2-λsin×(β2-β0),使用步骤6、步骤7相同方法计算得到β4处的值x4,其中λsin=0.6;附表1中已列出不同ζ、不同a对应的β4、 x4的值,可直接使用;
步骤12、采用步骤11获得的基础振动激励X0-NEW进行振动试验。
需要注说明的是,对于随机振动试验,基础振动激励条件的下凹可同样参照步骤2~步骤13;不同之处是使用公式(2)求解频点β1、β2的值,且β3、β4 的定义如下:
β3=β1+λrand×(β0-β1)
β4=β2-λrand×(β2-β0)
其中,λrand=0.5。附表2中已列出不同ζ、不同a对应的β1、β2、β3、β4、 x3、x4的值,可直接使用;
附表2:随机振动参数
本发明使用的前提条件是可将试验产品简化为图1所示的系统。对于复杂的多自由度系统,在低频段,可将系统的前几阶(一般取前3-4阶)主模态分别简化为图1所示的振动系统,然后在每阶主模态附近分别使用本发明进行条件下凹。
图6表示的是正弦振动试验中,当下凹参数a=0.6、产品阻尼比ζ=0.05时使用传统矩形下凹方法和使用本发明的三角形下凹方法后,计算得到的产品无量纲归一化响应。由图6可知,在三角形下凹方法下,产品响应在响应限制附近有微小浮动,既可保证产品响应不会太大以致造成产品破坏,又可保证产品响应不会太小以致造成“欠试验”的情况。
图7表示在正弦振动试验中,对于同一个产品(产品不变则阻尼比ζ不变),变化下凹参数a,传统矩形下凹与三角形下凹得到的响应对比示意图。图8为正弦振动试验中,对不同产品(阻尼比不同),使用相同的下凹参数a,传统矩形下凹与三角形下凹得到的响应对比示意图。两图都能得到如图6一样的结论。
图9~图11是图6~图8分别对应的随机振动的情况,也能得到如图6相同的结论。
本发明提出的一种基于传函的基础激励振动条件的三角形下凹方法,利用系统传递函数可快速求得响应限制对应的频点β1、β2,利用频点β1、β2及β0构造出二次曲线,进而构造出基础激励的三角形下凹方案。在本发明的三角形下凹方法下,产品响应在响应限制附近有微小浮动,既可保证产品响应不会太大以致造成产品破坏,又可保证产品响应不会太小以致造成“欠试验”的情况。
综上所述,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (2)
1.一种振动试验方法,其特征在于,包括如下步骤:
步骤1、建立正弦振动系统传递函数:
其中,X0为振动台M0的位移,X为试验产品m的位移;j为复数单位,β为归一化频率;ζ为试验产品的阻尼比;
步骤2、设正弦振动试验中位移X的最大值为Xmax,响应限制为Xlim;
步骤6、在基础激励条件坐标系中,过(β1,1)和(β0,a)两点并以点(β0,a)作为极值点,构造一条二次曲线;
步骤7、令β3=β1+λsin×(β0-β1),求出步骤6构造的二次曲线在β3处的值x3;
步骤9、过(β2,1)和(β0,a)两点并以点(β0,a)作为极值点,构造一条二次曲线;定义β4=β2-λsin×(β2-β0),求出本步骤构造的二次曲线在β4处的值x4;其中λsin=0.6;
步骤12、采用步骤11获得的基础振动激励X0NEW进行振动试验。
2.一种振动试验方法,其特征在于,包括如下步骤:
步骤1、建立随机振动系统传递函数:
其中,j为复数单位,β为归一化频率;ζ为试验产品的阻尼比;
步骤2、设随机振动试验中试验产品的位移X的最大值为Xmax,响应限制为Xlim;
步骤6、在基础激励条件坐标系中,过(β1,1)和(β0,a)两点并以点(β0,a)作为极值点,构造一条二次曲线;
步骤7、令β3=β1+λrand×(β0-β1),求出步骤6构造的二次曲线在β3处的值x3;λrand=0.5;
步骤9、过(β2,1)和(β0,a)两点并以点(β0,a)作为极值点,构造一条二次曲线;定义β4=β2-λrand×(β2-β0),求出本步骤构造的二次曲线在β4处的值x4;
步骤12、采用步骤11获得的基础振动激励X0NEW进行振动试验。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910428714.5A CN110186632B (zh) | 2019-05-22 | 2019-05-22 | 一种基于三角形下凹的基础激励振动条件的振动试验方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910428714.5A CN110186632B (zh) | 2019-05-22 | 2019-05-22 | 一种基于三角形下凹的基础激励振动条件的振动试验方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110186632A CN110186632A (zh) | 2019-08-30 |
CN110186632B true CN110186632B (zh) | 2020-11-20 |
Family
ID=67717344
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910428714.5A Active CN110186632B (zh) | 2019-05-22 | 2019-05-22 | 一种基于三角形下凹的基础激励振动条件的振动试验方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110186632B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002202217A (ja) * | 2001-01-04 | 2002-07-19 | Nec Corp | 振動試験機の制御装置 |
CN106596005A (zh) * | 2016-07-15 | 2017-04-26 | 北京卫星环境工程研究所 | 振动控制系统中主动下凹力学条件制定模块 |
CN107782520A (zh) * | 2017-09-22 | 2018-03-09 | 上海卫星工程研究所 | 微小卫星随机振动试验下凹控制方法 |
CN109000868A (zh) * | 2018-05-07 | 2018-12-14 | 北京空间飞行器总体设计部 | 一种航天器正弦振动试验下凹条件制定方法 |
-
2019
- 2019-05-22 CN CN201910428714.5A patent/CN110186632B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002202217A (ja) * | 2001-01-04 | 2002-07-19 | Nec Corp | 振動試験機の制御装置 |
CN106596005A (zh) * | 2016-07-15 | 2017-04-26 | 北京卫星环境工程研究所 | 振动控制系统中主动下凹力学条件制定模块 |
CN107782520A (zh) * | 2017-09-22 | 2018-03-09 | 上海卫星工程研究所 | 微小卫星随机振动试验下凹控制方法 |
CN109000868A (zh) * | 2018-05-07 | 2018-12-14 | 北京空间飞行器总体设计部 | 一种航天器正弦振动试验下凹条件制定方法 |
Non-Patent Citations (1)
Title |
---|
卫星随机振动带谷试验方法探讨;邱汉平等;《强度与环境》;20111231;第36卷(第6期);45-49 * |
Also Published As
Publication number | Publication date |
---|---|
CN110186632A (zh) | 2019-08-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Cunefare et al. | Energy harvesting from hydraulic pressure fluctuations | |
US7911205B2 (en) | Electromagnetic resonance frequency inspection systems and methods | |
Song et al. | Temporal kurtosis of dynamic pressure signal as a quantitative measure of combustion instability | |
US10371772B2 (en) | Magnetic resonance imaging apparatus | |
JP2011158444A (ja) | 二次電池の残存容量検出方法および装置 | |
US9182379B2 (en) | Method of carrying out a vibratory fatigue test of a mechanical part | |
CN109632229A (zh) | 共振疲劳测试方法、装置及工程疲劳测试平台 | |
CN110849568B (zh) | 一种结构疲劳寿命的试验方法 | |
CN110186632B (zh) | 一种基于三角形下凹的基础激励振动条件的振动试验方法 | |
CN103018027B (zh) | 一种激发飞机刹车阀故障的方法 | |
CN107144720B (zh) | 循环载荷下高温超导带材临界电流测试装置及其测试方法 | |
CN106815437B (zh) | 变压器稳态工况下油箱振动敏感区域确定方法及装置 | |
CN111881564A (zh) | 一种关于机械结构变幅疲劳寿命预测方法 | |
US9835034B2 (en) | Method for detuning a rotor-blade cascade | |
CN105092080B (zh) | 一种应用在声表面波无线无源测温系统中的新型扫频算法 | |
Holzinger et al. | Self-excited blade vibration experimentally investigated in transonic compressors: Acoustic resonance | |
CN110375989A (zh) | 一种柴油机燃烧噪声检测系统 | |
Nelson et al. | Three dimensional experimental characterization of a NiMnGa alloy | |
Yang et al. | Experimental investigation on vibration isolation effectiveness of two sets of underwater vehicle dynamic systems | |
CN109992835A (zh) | 一种基于Fatemi-Socie准则的增压器叶片疲劳寿命评估方法 | |
CN104281757B (zh) | 抗振晶振的机电集成设计方法 | |
Kotkowski et al. | Diagnostics of turbine blades, based on estimation of frequency response function | |
Berruti et al. | A test rig for non-contact travelling wave excitation of a bladed disk with underplatform dampers | |
JP6956597B2 (ja) | 磁気共鳴イメージング装置 | |
Le et al. | VIBRATION CHARACTERISTIC CHANGES IN DAMAGED STRUCTURE WITH REDUCTION OF CROSS-SECTION |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |