CN102230821B - 超低频振动台振级快速调整方法 - Google Patents

超低频振动台振级快速调整方法 Download PDF

Info

Publication number
CN102230821B
CN102230821B CN2011100905024A CN201110090502A CN102230821B CN 102230821 B CN102230821 B CN 102230821B CN 2011100905024 A CN2011100905024 A CN 2011100905024A CN 201110090502 A CN201110090502 A CN 201110090502A CN 102230821 B CN102230821 B CN 102230821B
Authority
CN
China
Prior art keywords
frequency
output
shaking table
signal
voltage
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
Application number
CN2011100905024A
Other languages
English (en)
Other versions
CN102230821A (zh
Inventor
何闻
王志鹏
于梅
沈润杰
马明德
贾叔仕
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang University ZJU
National Institute of Metrology
Original Assignee
Zhejiang University ZJU
National Institute of Metrology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Zhejiang University ZJU, National Institute of Metrology filed Critical Zhejiang University ZJU
Priority to CN2011100905024A priority Critical patent/CN102230821B/zh
Publication of CN102230821A publication Critical patent/CN102230821A/zh
Application granted granted Critical
Publication of CN102230821B publication Critical patent/CN102230821B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Apparatuses For Generation Of Mechanical Vibrations (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

超低频振动台振级快速调整方法,包括以下步骤:使振动台系统工作于位移反馈方式,采用保证频率变化时输出信号相位连续且幅值不变的信号发生器作为信源;使用逐次逼近法在预调整频率fk下调整振动台输出至与设定振动位移x在允许的误差范围内;使用逐步移频法将信号发生器输出等幅地、无间断相位地、从预调整频率fk逐步地递减至试验频率f,等待振动台输出稳定,调整过程结束。本发明具有使振动台在超低频工作时可以迅速精确地调整到设定振级的优点。

Description

超低频振动台振级快速调整方法
技术领域
本发明涉及一种超低频振动台振级快速调整方法。
技术背景
随着科学技术的不断发展,传感器的测试水平不断朝着极限水平发展。在航空航天、地震与海啸监测、精密加工、大型建筑的结构动力学研究、石油勘探等领域,人们对超低频测振传感器的测量频率范围不断提出新的要求,如下限频率低于0.05Hz的速度传感器以及零频响应的加速度传感器,都已经开始应用到各领域当中。这些传感器的出现对超低频计量校准系统提出了更高的要求。
常规的超低频校准是一项非常费时繁琐的工作,为确保校准精度及工作效率,缩短传感器校准过程的时间是一项非常有意义的工作。对于超低频振动校准自动控制系统,必须设计一种合理的振级调整方法,使振动台可以快速精确的达到预定振级。目前比较普遍的振级调整方法是逐次逼近法,这种方法的基本思路是:首先采集得到标准传感器输出当前振动台若干个周期的振动信号,然后根据电压有效值的计算方法得到振动台当前振级,然后以设定步长逐次递增信号发生器发出信号的幅度,直到振动台稳定在预定振级下。但是,超低频(<0.1Hz)信号的测量时间随信号周期的增大而增长,如在0.01Hz下对传感器进行校准,每改变一次信号发生器的输出信号幅值,就需要至少一个周期即100s的测量时间,假设经过5步调整可以达到预定振级,则至少需要500s的调整时间,传感器校准过程非常缓慢。
发明内容
为克服现有技术在超低频振动校准过程中测试时间长且效率较低的缺点,本发明提供了一种使振动台在超低频工作时迅速精确地达到设定振级的调整方法。
超低频振动台振级快速调整方法具体步骤如下:
1)、使振动台工作于位移反馈方式;采用保证频率变化时输出信号相位连续且幅值不变的信号发生器作为信号源,如基于直接数字合成(DDS)技术的信号发生器;
2)、获取用户设置的试验频率f和试验振动位移x,设定电压增量修正系数s(0<s<1),频率调整步数n(n>1);
3)、设置信号发生器输出频率fs为预调整频率fk、电压ui为初始uk
4)、控制信号发生器向振动台系统输入信号,待振动台输出稳定后,测取振动台的当前振动位移x0
5)、使用逐次逼近法调整振动台输出至允许的误差ε范围内:
(5.1)、计算当前比例因子 k = x 0 u i ;
(5.2)、计算信号发生器输出目标电压us
Figure BDA0000054711720000022
计算目标电压us与当前输出电压ui的电压增量Δu=us-ui;使信号发生器的输出电压 u i ⇐ u i + s × Δu ;
(5.3)、待振动台输出稳定后,测量当前振动位移x0
(5.4)、判断|x0-x|<ε是否成立,若不成立,重复执行步骤(5.1)-(5.3);若成立,则跳转至步骤6);
6)、使用逐步移频法将信号发生器的输出等幅地、无间断相位地、从fk逐步地递减至试验频率f;
7)、等待振动台稳定输出,结束振级调整。
进一步,步骤2)中,电压增量修正系数s根据试验设定,且满足0<s<1,如s=0.5;频率调整步数n也根据试验设定,且满足n>1,如n=10。
进一步,步骤3)中,fk为预调整频率,,式中,vmax为振动台输出最大速度峰值,d为设定的振动位移峰值;uk为初始电压,根据试验设定。
进一步,步骤6)中,使用逐步移频法将信号发生器的输出等幅地、无间断相位地、从fk逐步地递减至试验频率f包括以下步骤:
(6.1)、获取频率递减步长Δf=(fk-f)/n;
(6.2)、令f′=fs-Δf,其中fs为信号发生器的当前频率;
(6.3)、判断f′是否小于或等于试验频率f,若否,则将f′设为信号源频率,即
Figure BDA0000054711720000032
然后为使系统稳定,延时1s后重复执行
(6.4);若是,则将f设为信号源频率,即
Figure BDA0000054711720000033
退出移频操作至步骤7)。
本发明的技术构思是:振动台系统工作于位移反馈方式,信源采用保证频率变化时输出信号相位连续且幅值不变的信号发生器;然后使用逐次逼近法在一个相对较高的预调整频率fk下将振动台位移精确地、较快速地调整至设定位移;接着,使用逐步移频法无间断相位逐步递减信号发生器输出信号频率(幅值保持不变),直至设定频率,此时振动台的输出位移将基本保持为开始在fk下控制好的位移。
由于移频速度远远快于超低频下振级调整速度,所以可以对超低频振级进行快速调整。以在0.01Hz下对传感器进行校准为例,首先在预调整频率fk=0.1Hz下按照逐次逼近法将振动台调整到设定振动位移,假设每改变一次信号发生器的输出信号幅值,需要一个周期即10s的测量时间,并经过5步调整达到预定值,则共需50s时间;然后,将信号发生器的输出频率通过10步移频调整、花费10s时间调整到0.01Hz。总共调整时间只需60s,比技术背景中描述的常规方法调整时间缩短了近一个数量级。
本发明具有使振动台在超低频工作时迅速精确地达到设定振级的优点。
附图说明
图1是超低频振动台位移反馈控制模型。
图2是振动台运动部件简化动力学模型。
图3是振动台机电耦合动力学模型。
图4是反馈控制模型幅频特性曲线。
图5是超低频振动校准系统组成示意图。
图6是振级调整算法流程图。
具体实施方式
参照附图,进一步说明本发明:
超低频振动台振级快速调整方法具体步骤如下:
1)、使振动台工作于位移反馈方式;采用保证频率变化时输出信号相位连续且幅值不变的信号发生器作为信号源,如基于直接数字合成(DDS)技术的信号发生器;
2)、获取用户设置的试验频率f和试验振动位移x,设定电压增量修正系数s(0<s<1),频率调整步数n(n>1);
3)、设置信号发生器输出频率fs为预调整频率fk、电压ui为初始uk
4)、控制信号发生器向振动台系统输入信号,待振动台输出稳定后,测取振动台的当前振动位移x0
5)、使用逐次逼近法调整振动台输出至允许的误差ε范围内:
(5.1)、计算当前比例因子 k = x 0 u i ;
(5.2)、计算信号发生器输出目标电压us
Figure BDA0000054711720000052
计算目标电压us与当前输出电压ui的电压增量Δu=us-ui;使信号发生器的输出电压 u i ⇐ u i + s × Δu ;
(5.3)、待振动台输出稳定后,测量当前振动位移x0
(5.4)、判断|x0-x|<ε是否成立,若不成立,重复执行步骤(5.1)-(5.3);若成立,则跳转至步骤6);
6)、使用逐步移频法将信号发生器的输出等幅地、无间断相位地、从fk逐步地递减至试验频率f;
7)、等待振动台稳定输出,结束振级调整。
步骤2)中,电压增量修正系数s根据试验设定,且满足0<s<1,如s=0.5;频率调整步数n也根据试验设定,且满足n>1,如n=10。
步骤3)中,fk为预调整频率,
Figure BDA0000054711720000061
式中,vmax为振动台输出最大速度峰值,d为设定的振动位移峰值;uk为初始电压,根据试验设定。
逐步移频基本原理
振动台在超低频工作时采用如1图所示的位移反馈控制闭环调节系统。图中,功率放大器的传递函数在超低频时可以认为是一个常数,即有
G 1 ( s ) = U o ′ ( s ) U i ′ ( s ) = K - - - ( 1 )
在超低频工作条件下,振动台运动部件可简化为如图2所示单自由度机械动力学系统模型,并进而得到图3所示振动台机电耦合动力学模型,其电压-位移传递函数可表示为
G 2 ( s ) = X o ( s ) U o ′ ( s ) = Bl mLs 3 + ( mR + cL ) s 2 + [ Rc + kL + ( Bl ) 2 ] s + Rk - - - ( 2 )
式中,m为运动部件及负载总质量,k是运动部件支撑系统的刚度,c为运动部件的阻尼系数,L为动圈等效电感,R为动圈等效电阻,l为动圈等效长度。
在超低频工作时,位移传感器的传递函数H(s)可认定为常数
H ( s ) = U f ( s ) X o ( s ) = S g - - - ( 3 )
控制器采用比例控制,即
G3(s)=Kp    (4)
式中,Kp为比例系数
根据上述对图1振动台反馈控制系统各个环节的介绍,可以给出整个闭环反馈系统的传递函数
G ( s ) = X o ( s ) U i ( s ) = K p KBl mLs 3 + ( mR + cL ) s 2 + [ Rc + kL + ( Bl ) 2 ] s + Rk + BlK p KS g - - - ( 5 )
令s=jω,则有系统频响函数
G ( ω ) = K p KBl [ - ( mR + cL ) ω 2 + Rk + BlK p KS g ] + j [ - mL ω 2 + Rc + kL + ( Bl ) 2 ] ω - - - ( 6 )
当振动台在超低频直至零频的条件下工作时,公式中的ω项趋于零,通过调整控制器与功率放大器参数K、Kp,使分母中BlKpKSg?Rk,则振动台系统超低频工作时的电压-位移频响函数进一步可简化为
G ( ω ) ≈ K p KBl BlK p KS g = 1 S g - - - ( 7 )
从式中可以看出,此时系统频响函数只与位移传感器灵敏度有关,而与工作频率无关。
特别地,振动台取典型参数为:R=0.23Ω,B=0.58T,l=38m,c=3.78N/m·s-1,m=9kg,L=0.00036H,k=360N/m;功率放大器增益取K=10;控制器比例系数取Kp=1;位移传感器电压灵敏度取Sg=100V/m。将这些参数代入式(6),得
G ( ω ) = 220.4 ( - 2.0713608 ω 2 + 22122.8 ) + j ( 0.00324 ω 2 + 486.7606 ) ω - - - ( 8 )
由此得到系统幅频与相频曲线,如图4所示。从图中可以看出,当振动台处于超低频(<0.1Hz)条件下时,幅频曲线基本为一直线,表明在等电压的情况下,位移振动量基本为一常量;而相频特性接近于0°。这充分证明了式(7)近似式的有效性。通过所建立的数学模型式(8),在等电压输入条件下,取0.1Hz到接近零频的一组典型频率点值计算得到振动台输出位移,如表1所示。从表中数据可知,在等电压输入情况下,直至0.0001Hz的系统输出位移量相对于0.1Hz的输出位移量偏差极小。由此可知,采用位移反馈的振动台系统在超低频工作时输出位移量主要取决于系统的输入信号电压,而输入信号的频率对输出位移影响不大。这样,当振动台进行超低频振级调整时,先在较高的预调频率fk(如0.1Hz)下采用逐次逼近法将振动台输出位移较快速而精确地调整至设定值,然后在保持信号发生器输出幅值不变地情况下逐步无间断相位地将信号发生器输出频率递减至设定频率,则振动台输出位移将基本保持fk下幅值,从而实现超低频段时振级的快速调整。为了提高系统的稳定性,信号发生器必须保证频率变化时输出信号相位连续、幅值不变。
表1等电压超低频输入振动台输出位移变化趋势
Figure BDA0000054711720000081
下面参照附图5对本发明具体描述:系统由超低频标准振动台、功率放大器、信号发生器、多通道同步数据采集卡、位移传感器、计算机及控制软件等组成。位移传感器测量得到振动台输出位移信号。信号发生器产生一定频率和幅值的标准正弦信号,和位移传感器输出比较、通过控制器运算后输入至功率放大器,激励振动台工作,实现位移反馈工作模式。信号发生器采用DDS技术,输出幅值、频率均独立可调,保证频率变化时输出信号相位连续且幅值不变。数据采集卡测量位移传感器输出信号幅值,以作为逐次逼近法控制振级的依据;在振级调整完成后,数据采集卡同时测量被校传感器经适配器后的输出电压信号,从而计算得到被较传感器的灵敏度。信号发生器和数据采集卡均通过相应的接口与计算机通讯,并受控于计算机。
图6为振级调整算法流程图,包括如下步骤:使振动台系统工作于位移反馈方式,采用保证频率变化时输出信号相位连续且幅值不变的信号发生器作为信源;读取用户设置的试验频率f与试验振动位移x,根据试验设定电压增量修正系数s与频率调整步数n,设置信号发生器频率为预调整频率fk,初始电压为uk;使用逐次逼近法将振动台输出精确调整至设定位移值x;使用逐步移频法将信号发生器的输出等幅地、无间断相位地、从fk逐步地递减至试验频率f,从而完成超低频振级快速调整。频率增量可表示为Δf=(fk-f)/n,为保证频率变化时系统输出稳定,设置每次频率调整后等待1s。
本说明书实施例所述的内容仅仅是对发明构思的实现形式的列举,本发明的保护范围不应当被视为仅限于实施例所陈述的具体形式,本发明的保护范围也及于本领域技术人员根据本发明构思所能够想到的等同技术手段。

Claims (2)

1.超低频振动台振级快速调整方法,包括以下步骤:
1)、使振动台工作于位移反馈方式;采用保证频率变化时输出信号相位连续且幅值不变的信号发生器作为信号源;
2)、获取用户设置的试验频率f和试验振动位移x,设定电压增量修正系数s(0<s<1),频率调整步数n(n>1);
3)、设置信号发生器输出频率fs为预调整频率fk、电压ui为初始电压uk
4)、控制信号发生器向振动台系统输入信号,待振动台输出稳定后,测取振动台的当前振动位移x0
5)、使用逐次逼近法调整振动台输出至允许的误差ε范围内:
(5.1)、计算当前比例因子 
Figure FDA00001909944300011
(5.2)、计算信号发生器输出目标电压us, 
Figure FDA00001909944300012
计算目标电压us与当前输出电压ui的电压增量Δu=us-ui;使信号发生器的输出电压 
Figure FDA00001909944300013
(5.3)、待振动台输出稳定后,测量当前振动位移x0
(5.4)、判断|x0-x|<ε是否成立,若不成立,重复执行步骤(5.1)-(5.3);若成立,则跳转至步骤6);
6)、使用逐步移频法将信号发生器的输出等幅地、无间断相位地、从fk逐步地递减至试验频率f;
(6.1)、获取频率递减步长Δf=(fk-f)/n; 
(6.2)、令f′=fs-Δf,其中fs为信号发生器的当前频率;
(6.3)、判断f′是否小于或等于试验频率f,若否,则将f′设为信号源频率,即 
Figure FDA00001909944300021
然后为使系统稳定,延时1s后重复执行(6.2);若是,则将f设为信号源频率,即 
Figure FDA00001909944300022
退出移频操作至步骤7);
7)、等待振动台稳定输出,结束振级调整。
2.如权利要求1所述的超低频振动台振级快速调整方法,其特征在于:步骤2)中,电压增量修正系数s根据试验设定,且满足0<s<1;频率调整步数n也根据试验设定,且满足n>1。 
CN2011100905024A 2011-04-11 2011-04-11 超低频振动台振级快速调整方法 Active CN102230821B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2011100905024A CN102230821B (zh) 2011-04-11 2011-04-11 超低频振动台振级快速调整方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2011100905024A CN102230821B (zh) 2011-04-11 2011-04-11 超低频振动台振级快速调整方法

Publications (2)

Publication Number Publication Date
CN102230821A CN102230821A (zh) 2011-11-02
CN102230821B true CN102230821B (zh) 2012-11-07

Family

ID=44843409

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2011100905024A Active CN102230821B (zh) 2011-04-11 2011-04-11 超低频振动台振级快速调整方法

Country Status (1)

Country Link
CN (1) CN102230821B (zh)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102661782B (zh) * 2012-05-14 2013-10-23 中国航空工业集团公司北京长城计量测试技术研究所 一种超低频振动参数量值的快速测量方法
US9377375B2 (en) * 2012-05-16 2016-06-28 Venturedyne, Ltd. Repetitive shock vibration testing system and method
CN103822703B (zh) * 2014-03-19 2017-04-19 浙江大学 一种超低频水平向振动台导轨不平顺动态补偿方法
CN106289693B (zh) * 2016-08-17 2018-08-14 西安交通大学 一种液压振动系统的低频拓展控制方法
CN108426690A (zh) * 2018-01-22 2018-08-21 浙江大学 基于交叠残周期采样的振级快速调整方法
CN108444591A (zh) * 2018-01-22 2018-08-24 浙江大学 残周期采样确定振动台振级的方法
CN108444590A (zh) * 2018-01-22 2018-08-24 浙江大学 基于交叠周期采样的振级快速调整方法
CN108827457A (zh) * 2018-06-12 2018-11-16 太原理工大学 基于传递函数的振动幅值迭代控制方法
CN111256927B (zh) * 2019-12-31 2022-03-25 杭州亿恒科技有限公司 自适应调节双闭环时间正弦振动控制方法
CN113932917B (zh) * 2021-11-16 2022-06-10 厦门乃尔电子有限公司 一种振动传感器校准系统输出信号控制方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2223963A5 (en) * 1973-03-05 1974-10-25 Commissariat Energie Atomique Ultra-low frequency vibration generator - has opposed linear motors with common armature and variable three-phase current
CN101165460A (zh) * 2006-10-20 2008-04-23 上海瑞视仪表电子有限公司 电涡流位移/振动传感器系统校验仪

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2223963A5 (en) * 1973-03-05 1974-10-25 Commissariat Energie Atomique Ultra-low frequency vibration generator - has opposed linear motors with common armature and variable three-phase current
CN101165460A (zh) * 2006-10-20 2008-04-23 上海瑞视仪表电子有限公司 电涡流位移/振动传感器系统校验仪

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
Asan Gani等.A labview based data acquisition system for vibration monitoring and analysis.《2002 Student conference on research and development proceedings会议论文集》.2002,全文. *
H.I. Abu-mulaweh等.The use of labview and data acquisition unit to monitor and control a bench-top air-to-water hear pump.《Computer applications in engineering education》.2008,第16卷(第2期),全文. *
与梅.低频振动传感器幅值和相移测量的不确定度.《振动与冲击》.2009,第28卷(第4期),全文. *
于梅.低频超低频振动计量技术的研究与展望.《振动与冲击》.2007,第26卷(第11期),全文. *
于梅.精密仪器环境振动测量和评价方法的研究.《振动与冲击》.2010,第29卷(第8期),全文. *

Also Published As

Publication number Publication date
CN102230821A (zh) 2011-11-02

Similar Documents

Publication Publication Date Title
CN102230821B (zh) 超低频振动台振级快速调整方法
CN103255755B (zh) 一种快速实时评价土石料填筑压实质量的无损方法及其评价装置
JP5340122B2 (ja) 共振測定システムの作動方法および共振測定システム
CN104035329B (zh) 一种振动自动化校准方法
CN101666730B (zh) 一种使用测试空气弹簧振动特性设备的测试方法
CN101685026B (zh) 标定倾角传感器敏感轴零位输出值的方法及获取其参数的方法
CN101878417A (zh) 发动机台架系统控制系统
CN101382425A (zh) 一种微机械陀螺自激驱动与解调装置
CN105573256A (zh) 增益自动调整支援装置
CN108645583A (zh) 一种装配结合面法向接触阻尼高精度检测装置及方法
CN106052667A (zh) 振动陀螺仪中谐振器和科里奥利轴控制的系统、装置、方法
CN105547274A (zh) 一种用于机抖激光陀螺的主动式抖动消振控制方法
CN108388109A (zh) 回收机控制系统及方法
CN102818578A (zh) 用于启动期间的增强速率估计的陀螺仪动态电动机振幅补偿
CN106918438B (zh) 一种多分量力及力矩的测量方法及系统
CN105571576A (zh) 一种mems陀螺模态匹配电压自动测试方法
CN103884870A (zh) 提高加速度计校准精度的方法和装置
CN102175394B (zh) 刚性转子软支承动不平衡测试中的永久标定方法
CN106354016A (zh) 一种金属管夯实高度的控制方法
US20160124097A1 (en) Seismic shaker
CN104597753A (zh) 一种同步碎石封层车沥青碎石撒布智能控制方法及装置
US10281269B2 (en) Method and device for determining a height of a settled bed in a mixture in a loading space
CN103575469A (zh) 飞机可移动控制表面的惯性性质的测量
CN102648392B (zh) 通过转动惯性设备确定航向的方法
CN202770590U (zh) 加工机的动平衡刀把及其动平衡控制系统

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant