专利名称：振动台冲击响应生成方法
技术领域：
本发明涉及振动台控制技术，具体说就是一种振动台冲击响应生 成方法。
(二)
背景技术：
为了保证产品承受非重复性机械冲击的适应性，许多产品在出厂 前都需要经受冲击试验的考核，尤其是军用设备。冲击试验的目的主 要是用来确定元件、设备在使用和运输过程中经受非多次重复的机械 冲击的适应性，以及评价结构的完好性，其次还可以用于微电子器件 的内强度试验。继电器等开关类器件更容易受到冲击的影响，圆军标
65B对继电器有详细规定，当继电器按国军标360A中方法213规定 的条件A 、 B或C进行试验时，闭合触点的断开不得超过10us，断 开触点的闭合或桥接不得超过lus，不允许有机械或电气损坏现象。
电动式振动台是目前使用最广泛的一种振动测试设备。电动式振 动台的频率范围宽，小型振动台频率范围为0-10kHz;动态范围宽， 易于实现自动或手动控制，加速度波形良好，适合产生随机波，可得 到很大的加速度。虽然振动台不如冲击机产生的冲击加速度大，但由 于振动台良好的可控性，可以产生更符合标准的冲击响应。
冲击振动控制就是依据振动系统特性，针对设定的冲击波参数通 过适当的算法产生控制信号，激励电动式振动台再现设定的冲击波 形。方法213中规定的是一种复现冲击波形方法，是时域波形在振动 台上的再现。针对时域波形控制，传统方法是采用波形幅值均衡法和 传递函数均衡法，此类方法对所模拟的典型冲击脉冲时域波形进行修 正及控制，利用数字式控制方法实现对振动台控制点处的冲击响应信 号的精确控制。在控制过程中，不断修正时域波形的幅值、延迟时间 等相关参数或不断确定试验系统的传递函数，再利用规定波形的傅里 叶变换得到新的驱动波形信号，从而完成多次重复的、闭环的修正过程。此类方法是通过实时试验并修正控制信号的方法实现冲击响应， 一方面，此类方法过程繁琐并需要傅里叶变换、闭环控制等，整个控 制系统对硬件的要求被提高；另一方面，在冲击试验修正的过程中，
模拟冲击波形的一些参数，如幅值、时间延迟等将受到噪声和设备等 因素的影响发生变化，虽然冲击波形满足预期要求，但是不能保证在 冲击试验中冲击产生时，每一次的驱动信号都相同，这样会对实验结 果引入无法量化的扰动。
对振动台实现良好的控制，得到其传递函数是非常关键的，但是 目前对振动台参数的辨识研究较少，只有一些基于测量频率特性而估 计参数的研究方法，且操作繁琐，不利于工程应用。
发明内容
本发明的目的在于提供一种所生成的冲击响应具有利用系统资 源少、每次冲击下响应一致、可靠性高的振动台冲击响应生成方法。
本发明的目的是这样实现的所述的振动台冲击响应生成方法步 骤如下
一、 分析振动台的力学模型，得到其理论传递函数；根据理论传 递函数推导其离散传递函数，进而确定其零极点个数；
二、 根据辨识波形图所示形式作为振动台输入信号，通过振动台 控制驱动系统加载到振动台；
三、 加速度信号作为输出从振动台台面上的加速度传感器输出到 可存储波形的示波器，示波器同时记录振动台输入信号；
四、 将上一步输入输出数据保存至上位机，利用Matlab系统辨 识工具箱的信号处理工具处理输入输出数据；
五、 在Matlab系统辨识辨识工具箱中选择ARX模型，零极点个
数按第二步结果选择；
六、 进行系统辨识，如果辨识结果与实际波形拟合不好，改变第 四步信号处理时滤波环节的阀值，直至拟合成功(拟合度超过60% 可认定拟合成功)，拟合成功会得到振动台的离散传递函数；
七、 通Matlab的d2c命令将离散传递函数转换成连续传递函数；八、 通过理论模型形式修正系统辨识的结果，用数值代替原有参
数；
九、 将数据导入simulink参数估计，模型应用理论模型，参数初 值采用上步获得数值；
十、进行参数估计，如果估计不收敛，改变迭代规则，直至计算 收敛，最终得到接近真实模型的传递函数；
i^一、根据得到的传递函数，在simulink中进行PID仿真计算， 直至其生成符合国军标规定的波形复现法的冲击响应；
十二、将仿真结果应用到实际，开环控制即满足要求；也可以将 PID参数应用到闭环控制中，最终生成符合国军标要求的冲击响应。
本发明提供一种振动台冲击响应生成方法，该方法可以得到振动 台准确的传递函数；在得到传递函数的情况下，通过仿真可得到理想 的控制方法；本发明生成的冲击响应严格符合国军标要求，并具有利 用系统资源少，每次冲击下响应一致，可靠性高的特点。
(四)
图l为本发明的辨识波形图2为本发明的系统辨识流程图3为本发明的振动台控制驱动系统结构方框图4为本发明的电动式振动台等效参数模型图5为本发明的6ms宽半正弦冲击输入和输出响应图6为本发明的辨识出的冲击的仿真输出与实际输出图7为本发明的模型的PID控制方框图8为本发明的开环控制的输出与输入图。
具体实施方式
下面结合附图举例对本发明作进一步说明。
实施例1:结合图1、图2，本发明一种振动台冲击响应生成方 法，所述的冲击响应生成方法步骤如下
一、分析振动台的力学模型，得到其理论传递函数；根据理论传 递函数推导其离散传递函数，进而确定其零极点个数；二、 根据辨识波形图(图l)所示形式作为振动台输入信号，通
过振动台控制驱动系统加载到振动台；
三、 加速度信号作为输出从振动台台面上的加速度传感器输出到 可存储波形的示波器，示波器同时记录振动台输入信号；
四、 将上一步输入输出数据保存至上位机，利用Matlab系统辨 识工具箱的信号处理工具处理输入输出数据；
五、 在Matlab系统辨识辨识工具箱中选择ARX模型，零极点个 数按第二步结果选择；
六、 进行系统辨识，如果辨识结果与实际波形拟合不好，改变第 四步信号处理时滤波环节的阀值，直至拟合成功(拟合度超过60% 可认定拟合成功)，拟合成功会得到振动台的离散传递函数；
七、 通Matlab的d2c命令将离散传递函数转换成连续传递函数；
八、 通过理论模型形式修正系统辨识的结果，用数值代替原有参
数；
九、 将数据导入simulink参数估计，模型应用理论模型，参数初 值采用上步获得数值；
十、进行参数估计，如果估计不收敛，改变迭代规则，直至计算 收敛，最终得到接近真实模型的传递函数；
十一、根据得到的传递函数，在simulink中进行PID仿真计算， 直至其生成符合国军标规定的波形复现法的冲击响应要求；
十二、将传真结果应用到实际，开环控制即满足要求；也可以将 PID参数应用到闭环控制中，最终生成符合国军标要求的冲击响应。
实施例2，结合图l、图2、图3、图4、图5，本发明振动台冲 击响应生成方法，具体实施步骤如下
选择生成冲击的振动台是某国产小型电动式振动台。
第一步，分析振动台的力学模型，得到其理论传递函数；根据理论传
递函数推导其离散传递函数，进而确定其零极点个数；(结合图3)
□"。,+化+"e (2)
=肌=祝， (3)
0D []
F = 5￡/。=附x+cx十Ax (4)
式中^一输入电压；Z。一滤波器等效电感；C。一滤波器等效电容； i ,一滤波器等效电阻；^一线圈电流；"。一线圈电压；等效电阻； 丄一等效电感；5—气隙磁密；F —电动力；电枢和工作台质量； c—阻尼比；yt—弹性系数；x—台面位移；/一电枢等效长度；"。一 振动台等效电压；/。一驱动电流。
通常情况下可忽略模型中滤波器作用，将式(3)、 (4)代入式(2)
中，可得
W =——XH--XH--XH--X
" 说 放 说 说
经拉氏变换，可得输入电压"。与振动加速度x的传递函数G(力
刚=见-^-r (6)
w。 0) (Ly+i )(ms^ + cs+A:)+s 将式(6)转换成零极点表示形式
G(s) =-^- (7)
系统连续传函为三阶，同时含有三个极点和三个零点。为满足数 字控制的要求，将理论模型的连续传递函数转换为离散传递函数，从
而完成系统辨识。采用零阶保持器^:串接连续传函，由z^换得
式中T —采样时间；a。 = —a3ewK ; a, = a3 (e+ + es'^ ); a2 = —a3(e''r + + ; a3 = (s, - s2)(s2 - - s3);<formula>formula see original document page 8</formula>得到离散域的零、极点数，主要是目的是获得极点数，零点个数
为3，极点个数为3。
第二步，根据图1的波形图所示形式作为振动台输入信号，通过振动
台控制驱动系统加载到振动台；
第三步，加速度信号作为输出从振动台台面上的加速度传感器输出到 可存储波形的示波器，示波器同时记录振动台输入信号；(结合图4) 第四步，将上一步输入输出数据保存至上位机，利用Matlab系统辨 识工具箱的信号处理工具处理输入输出数据；
第五步，在Matlab系统辨识辨识工具箱中选择ARX模型，零点个数
为3，极点个数为3，延迟选择l (表示无延迟)；
第六步，进行系统辨识，如果辨识结果与实际波形拟合不好，改变第
四步信号处理时滤波环节的阀值，直至拟合成功(拟合度超过60%
可认定拟合成功)，拟合成功会得到振动台的离散传递函数；(结合图
<formula>formula see original document page 8</formula>
第七步，通Matlab的d2c命令将离散传递函数转换成连续传递函数；
<formula>formula see original document page 8</formula>
第八步，通过理论模型形式修正系统辨识的结果，用数值代替原有参 数；
<formula>formula see original document page 8</formula>(11)
第九步，将数据导入simulink参数估计，模型应用理论模型，参数初 值采用上步获得数值；第十步，进行参数估计，如果估计不收敛，改变迭代规则，直至计算
收敛，最终得到接近真实模型的传递函数；
卿=-n-- (13)
4.8451xl0-、3 +0.028289 +13.6235+1308.7
第H"^—步，根据得到的传递函数，在simulink中进行PID仿真计算， 直至其生成符合国军标规定的波形复现法的冲击响应要求；(结合图 6)
第十二步，将传真结果应用到实际，开环控制即满足要求；也可以将 PID参数应用到闭环控制中，最终生成符合国军标要求的冲击响应。 (结合图7)
产生一个时间为6 ms的半正弦冲击，冲击之后15ms产生了一个 向下的冲击，这是因为图8所示是加速度波形，振动台台面在冲击后 会保持静止，加速度的积分值也就是振动台台面的运动速度应为零， 在仿真时PID闭环控制下，使得这个向下的补偿冲击延迟了 15ms， 满足了国军标对冲击的要求。
权利要求
1.一种振动台冲击响应生成方法，其特征在于振动台冲击响应生成方法步骤如下一、分析振动台的力学模型，得到其理论传递函数；根据理论传递函数推导其离散传递函数，进而确定其零极点个数；二、根据辨识波形图所示形式作为振动台输入信号，通过振动台控制驱动系统加载到振动台；三、加速度信号作为输出从振动台台面上的加速度传感器输出到可存储波形的示波器，示波器同时记录振动台输入信号；四、将上一步输入输出数据保存至上位机，利用Matlab系统辨识工具箱的信号处理工具处理输入输出数据；五、在Matlab系统辨识辨识工具箱中选择ARX模型，零极点个数按第二步结果选择；六、进行系统辨识，如果辨识结果与实际波形拟合不好，改变第四步信号处理时滤波环节的阀值，直至拟合成功(拟合度超过60％可认定拟合成功)，拟合成功会得到振动台的离散传递函数；七、通过Matlab的d2c命令将离散传递函数转换成连续传递函数；八、通过理论模型形式修正系统辨识的结果，用数值代替原有参数；九、将数据导入simulink参数估计，模型应用理论模型，参数初值采用上步获得数值；十、进行参数估计，如果估计不收敛，改变迭代规则，直至计算收敛，最终得到接近真实模型的传递函数；十一、根据得到的传递函数，在simulink中进行PID仿真计算，直至其生成符合国军标规定的波形复现法的冲击响应要求；十二、将仿真结果应用到实际，开环控制即满足要求；也可以将PID参数应用到闭环控制中，最终生成符合国军标要求的冲击响应。
全文摘要
本发明的目的在于提供一种所生成的冲击响应具有利用系统资源少、每次冲击下响应一致、可靠性高的振动台冲击响应生成方法。方法包括分析振动台的力学模型，得到其理论传递函数；根据辨识波形图所示形式作为振动台输入信号，通过振动台控制驱动系统加载到振动台；进行系统辨识，根据得到的传递函数，在simulink中进行PID仿真计算，直至其生成符合国军标规定的波形复现法的冲击响应；将仿真结果应用到实际，开环控制即满足要求；也可以将PID参数应用到闭环控制中，最终生成符合国军标要求的冲击响应。本发明生成的冲击响应严格符合国军标要求，并具有利用系统资源少，每次冲击下响应一致，可靠性高的特点。
文档编号G01M7/02GK101576432SQ200910072258
公开日2009年11月11日 申请日期2009年6月12日 优先权日2009年6月12日
发明者任万滨, 康云志, 梁慧敏, 健 王, 翟国富 申请人:哈尔滨工业大学