资源描述
虚拟振动试验研究、验证及软件开发
0 引言
在姿控动力系统研发的过程中,振动环境试验是极其重要的环节之一,对于验证姿控动力系统的可靠性有着至关重要的作用。目前,姿控动力系统的振动环境试验一般通过振动试验台来完成。通过振动台模拟姿控动力系统在工作状态中可能遇到的外部激励情况,根据试验结果,评估姿控动力系统的抗振能力和可靠性,从而指导轨姿控动力系统的改进设计。然而,基于振动台的振动环境试验存在一些局限性:1)振动环境试验是在产品已经设计完成后进行的,且要产生一个甚至多个试验样品或模型,一旦验证产品可靠性出现问题,又要重新进行改进设计,重新制造产品,重新进行试验,试验周期长且耗费巨大。2)由于试件安装在试验台上和安装在实际产品上的耦合关系和边界条件不同,因此试验台难以完全模拟真实的激励环境,容易引起过试验或欠试验问题,存在一定的试验风险;3)振动试验受到传感器数量和安装位置的限制,能够获得的试验结果信息有限,无法完全观测到产品的各个部位的响应特性,容易漏掉可能是关键部位的响应情况。
针对采用振动台进行振动环境试验的局限性,可以采用动力学仿真手段,建立一套虚拟振动试验系统。虚拟振动试验是指通过软件环境建立的振动试验台模型和被测试件模型,在软件环境中基于振动试验台和被测试件的虚拟仿真模型,完成与实际振动试验相同或类似的振动环境试验任务。虚拟振动试验可以在物理试验之前尽可能全面的了解并及时发现问题,避免在设计后期出现设计反复和颠覆性错误。基于此虚拟振动试验系统,可以建立振动试验台模型和试件模型,并将试件模型与振动台模型进行装配,完成所需的振动试验,然后根据虚拟振动试验结果,对试件设计方案提出改进和优化。此外,可以为物理试验做好充分准备,减少对物理试验的依赖,进一步加快产品研发进程。
1虚拟振动试验平台建设
基于联合仿真技术的虚拟振动试验平台建设的基本思路:首先辨识出电动振动台的机电耦合参数模型,通过Amesim软件建立振动台系统的机电耦合模型,然后,用Virtual.lab Motion建立台体及试验件的刚柔耦合模型,最终将Virtual.lab Motion、Amesim中的模型编译为Matlab/Simulink可识别的S-function函数,以Simulink为平台建立与Amesim软件、Virtual.lab软件之间的联合仿真接口,然后通过Matlab/GUI开发的软件界面进行参数设置及调用Simulink进行联合仿真,从而实现虚拟振动试验平台建设。
1.1虚拟振动试验平台建设
根据电动振动台系统组成与工作原理,可以采用集中参数的拉格朗日方程,将振动台结构的自由度降到有限值,进行足够精度的近似分析。动圈部分的一阶频率决定了振动台的工作频率上限,可将动圈沿轴向简化为1个二自由度系统,按其结构特点和功放分为等效台面部分和等效线圈部分,它们之间通过弹簧和阻尼连接;将振动台台体看作刚体,与地面之间通过隔振气囊连接,也可以简化为单自由度弹簧质量系统。这样,振动台的机械部分便缩减为三自由度集中参数模型。振动台运动过程中所受的力是由线圈中的交变电流决定的,而电流的大小又受线圈运动速度的影响,这也就是电动振动台的机电效应,如果不考虑线圈电阻和电感随频率和温度变化的影响,施加到振动台线圈的交变电压信号主要由3部分组成,一部分为电流通过线圈电阻产生的电压,另外一部分为线圈电感产生的磁感应电动势,还有一部分为线圈在磁场内运动产生的反电动势。
不考虑振动台台体运动对系统特性的影响,可以得到电动振动台机电耦合系统的数学模型:
(1)
式中,为台面质量,为线圈质量,为线圈和台面之间的刚度,为台面支撑气囊的刚度,、分别为台面和线圈的位移,、分别为动圈的等效电阻和电感,、分别为输入的电压和线圈中的等效电流,为力常数。
对式(1)进行Laplace数学变换,且考虑到机械元件的阻尼系数,则恒压、恒流驱动下的传递函数可表示为:
(2)
式中,为系统一阶惯性圆频率;为弹性支撑圆频率;为电谐振圆频率;为动圈轴向谐振圆频率;为台面弹性支撑阻尼比;为电阻尼比;为线圈与台面之间的阻尼比。
在垂直方向振动台空台和带负载情况下获得恒压、恒流驱动下的振动系统传递函数,根据传递函数可获得式(2)中的参数,这里给出我们经常使用的T2000振动台已经辨识出的参数,如表1所示:
表1 T2000振动台垂向模型参数
机械参数
电参数
耦合参数
Kf
mt/kg
mc/kg
kc/N∙m-1
ks/N∙m-1
cc/kg∙s-1
cs/kg∙s-1
R/Ω
L/H
21.76
23.90
2.3799e+9
2.5483e+6
500
800
5.1e-3
3.7699e-6
24.31
根据辨识出来的结果,在Amesim中建立振动台系统的机电耦合模型,如图1所示。
图1电动振动台的机电耦合模型
在Amesim软件中进行仿真运算,恒压条件下传递函数如图2所示,可以看出,虚拟振动台与真实振动台幅频曲线与相频曲线吻合较好,可以较为真实的模拟振动台的动力学特性。
相位/(°)
幅值/dB
频率/Hz
图2垂向空台试验与仿真传递函数对比
1.2电动振动台刚柔耦合模型的建立
利用改进的Craig-Bampton模态综合法与多体系统动力学混合建模法在Virtual.lab Motion中建立试件与振动台的刚柔耦合模型。具体建模过程如下:
1)建立姿控动力系统的CAD模型,根据CAD模型建立试件的刚体模型;
2)使用固定副(约束所有方向上的自由度)将试件的刚体模型与振动台动圈的刚体模型连接,完成被测试件与振动台的装配;
3)根据试件的CAD模型,使用有限元方法创建被试件的柔性体模型,确立其与振动台刚体的连接关系;
4)基于模态综合法求解盒式试件的模态集,定义各阶模态阻尼,选取合适阶数的模态数据作为试件柔性体的模态基空间;
5)导入试件的模态数据,完成试件的柔性化。至此,试件与振动台的刚柔耦合模型建立完成。
1.3正弦与随机控制算法
1.3.1 正弦控制算法
正弦振动控制主要是实现对响应信号幅值的控制,进行正弦振动控制仿真时,可采用如下方式提取响应信号的幅值:
考虑正弦信号,其可以描述为如下形式:
(3)
对于一系列的采样点,可得到如下的方程:
(4)
式中,表示第1个采样时刻,表示第n个采样时刻。根据上式,利用最小二乘法即可获得与的最优估计。正弦信号的幅值即可通过下式计算得出:
(5)
对于正弦信号,每一个给定的频率,系统在此频率下的传递函数可以写为:
(6)
式中,为响应信号的幅值,为驱动信号的幅值。
对于每一个频率点记参考信号谱在当前频率下的幅值为,根据当前频率点响应信号的幅值与试验参考幅值采用如下的算法确定下一个频率点系统的反传递函数。
(7)
其中,,为系统的反传递函数,为压缩因子。
则,下一个频率点的驱动信号幅值为:
(8)
1.3.2 随机控制算法
1、随机信号时域随机化
随机振动试验是基于在频域内对功率谱的控制,但振动台推动试件都是时域的随机信号,因此首先要将功率谱密度转化到时域随机信号,而如何生成高质量的试验随机驱动信号是要解决的重要问题。
用时域随机化生成随机信号的方法包括伪随机信号生成和真随机高斯信号生成两个部分。对于一维随机信号,假设其单边自谱密度为,用Cooley-Tukey计算方法可表示为:
(9)
其中是时域信号的傅氏变换,又称傅氏谱,表示傅氏变换的模,表示采样序列的长度,是采样的时间间隔。从式(9)可得到信号的幅值谱:
(10)
由于自谱密度缺乏相位信息,直接对式(10)做逆傅里叶变换得到的时间序列不能组成伪随机信号。因此必须加入相位随机化运算,令:
(11)
其中,是服从均匀分布的随机相位。
伪随机信号是具有一定时间长度的信号在时域上的周期性延拓,其频谱是离散谱,能量集中在原始信号采样频率上,信号长度以A/D采样时间为周期,幅值分布具有尖刺,近似于高斯分布。但工程实际中需要非周期的,幅值概率密度满足高斯分布的真随机信号。为此,需对上面获得的伪随机信号再经过时域随机化处理,具体过程如图3所示。
连续谱高斯分布
延时、旋转、加窗、叠加
离散谱非高斯分布
时域随机化
真随机信号
伪随机信号
图3真随机信号生成过程示意图
2、驱动谱反馈修正
随机振动控制的目的就是使控制谱与参考谱相一致。根据线性系统的控制理论有:
(12)
式中:是输入的自功率谱;
是输出的自功率谱;
是系统的传递函数;
在控制的过程中,要对上一个循环的驱动谱进行修正,其修正公式为:
(13)
式中:
是第次修正时的驱动谱;
是第次修正时的驱动谱;
是参考的自功率谱;
是控制的自功率谱;
1.4 虚拟振动试验平台建设
闭环虚拟试验平台建设的思路为:以Matlab/Simulink作为主计算平台,Amesim与Virtual.Lab中的模型都编译为Matlab/Simulink可识别的S-function。在Matlab/Simulink发送电压驱动信号到Amesim中,通过机电耦合系统将电压信号转化为力信号加载到台面(Virtual.Lab中)上,台面受到载荷作用后会产生位移、速度和加速度信号,台面的位移、速度信号又会传送给Amesim,对动圈产生力的作用,而台面的加速度信号经过幅值提取,将幅值信号发送到Matlab/Simulink中,在Matlab/Simulink中通过幅值控制算法对其电压信号进行调整,这就形成了一个闭环虚拟试验平台系统,具体的原理图如图4所示。
图4 虚拟振动试验平台建设原理图
以Simulink为平台建立与Amesim软件、Virtual.lab软件之间的联合仿真接口,如图5所示,然后通过软件界面设置控制参数进行正弦振动试验及随机振动试验,从而实现虚拟振动试验平台建设。
图5 虚拟振动试验平台软件实现
1.5软件界面设计
软件界面是通过MATLAB/GUIDE建立,主要分为三个界面,分别为主界面、正弦振动试验界面及随机振动试验界面,如图6-图8所示。主界面的菜单栏包含有限元模型、机电耦合模型、刚柔耦合模型、联合仿真接口设置、正弦扫描试验、随机振动试验及退出菜单,点击有限元模型菜单可以打开Patran软件,进行试件的有限元建模,点击机电耦合模型可以打开Amesim软件,进行振动台机电耦合系统的建模,点击刚柔耦合模型可以打开LMS Virtual.lab软件进行刚柔耦合模型的建立,点击联合仿真接口设置可以打开simulink建立的接口界面,点击正弦扫描试验进入正弦扫描试验界面,点击随机振动试验进入随机振动试验界面,点击退出则退出主界面;在正弦振动试验界面可以设置正弦振动谱及显示正弦控制响应谱与参考谱的对比图;在随机振动试验界面可以设置随机振动谱及显示正弦控制响应谱与参考谱的对比图。通过软件可以很方便的进行虚拟正弦振动及虚拟随机振动试验。
图1 虚拟振动试验平台主界面
图7正弦扫描试验界面
图8 随机振动试验界面
2 算例验证
2.1 试验产品及试验条件
对某型动力系统进行正弦扫描试验及随机振动试验,试验条件如表2和表3所示。
表1 正弦扫描试验条件
频率范围
量级
5Hz~200Hz
0.2g
表2 随机振动试验条件
频率范围
量级
总均方根值
5Hz~20Hz
0.0065g2/Hz
0.77g
120Hz
0.0002g2/Hz
121Hz~200Hz
0.003g2/Hz
240Hz
0.0015g2/Hz
340Hz
0.00003g2/Hz
500Hz
0.00015g2/Hz
2.2 虚拟振动试验过程
1、打开Matlab软件,在命令行输入xunizhendong,运行虚拟振动试验平台软件,如图9所示。
图9 运行虚拟振动试验平台软件
2、单击虚拟振动试验平台软件左上角的“有限元模型”菜单,启动Patran软件,建立姿控动力系统的有限元模型,并进行模型修正,生成有限元模型的.bdf文件。
3、单击虚拟振动试验平台软件左上角“机电耦合模型”菜单,启动Amesim软件,建立振动台的机电耦合模型,,并将打开到simulation mode状态下,如图10所示,
图10振动台机电耦合模型
4、单击虚拟振动试验平台软件左上角“刚柔耦合模型”菜单,启动LMS Virtual.lab软件,建立振动台与动力系统的刚性模型,通过patran中生成的有限元模型.bdf文件对姿控动力系统的刚性模型进行柔性化,得到刚柔耦合模型,运行软件得到一个plantout文件。
5、单击虚拟振动试验平台软件左上角“联合仿真接口设置”菜单,启动simulink界面,在simulink界面进行软件间的接口连接,如图11所示。
图11 联合仿真接口设置
6、单击虚拟振动试验平台软件左上角“正弦扫描试验”菜单,进入正弦扫描试验界面,设置正弦扫描试验条件,点击试验,虚拟正弦扫描试验就开始运行了,如图12所示。
图12 正弦扫描试验
7、单击虚拟振动试验平台软件左上角“随机振动试验”菜单,进入随机振动试验界面,设置随机振动试验条件,点击试验,虚拟随机振动试验就开始运行了,如图13所示。
图13 随机振动试验
2.3 虚拟振动试验结果
图14给出了正弦扫描试验的控制谱与参考谱之间的对比图,图15给出了随机振动试验的控制谱与参考谱之间的对比图,从比较的结果来看,控制效果较好,与真实的控制效果比较接近.图16给出了差动贮箱与金属隔膜贮箱之间安装处的响应,最大突频点是131Hz,加速度为9.951g,图17给出了振动台上进行试验差动贮箱与金属隔膜贮箱之间安装处的响应,最大突频点是131.95Hz,加速度为10.21g,通过比较可以看出,虚拟振动试验与真实振动试验计算结果吻合较好。
图14 正弦扫描试验的控制谱与参考谱
图15 随机振动试验的控制谱与参考谱
图2 虚拟试验与真实试验正弦扫描试验结果对比
图17 虚拟试验与真实试验随机振动试验结果对比
3 结论
本文对振动台的多学科集成仿真技术进行了研究,在不同软件中建立了振动台控制系统、机电耦合系统和刚柔耦合系统,通过联合仿真技术建立了闭环正弦扫描、随机振动虚拟振动试验平台。通过算例进行验证,与物理试验结果进行对比,表明了通过虚拟试验可以有效预测试验的结果。这一技术可以弥补物理振动试验的不足,降低试验风险,节省试验成本。
µ个人信息
Individual
µ自我评价
Evaluation
n 本人是市场营销专业毕业生,有丰富的营销知识体系做基础;
n 对于市场营销方面的前沿和动向有一定的了解,善于分析和吸取经验
n 熟悉网络推广,尤其是社会化媒体方面,有独到的见解和经验
n 个性开朗,容易相处,团队荣誉感强
&技能水平
Skills
n CET-6,优秀的听说写能力
n 计算机二级,熟悉计算机各项操作
n 高级营销员,国家职业资格四级
n 极具创意的广告策划、文案写作能力,
n 善于制作PPT文档,对于配色排版有独到的见解
¹实践经历
Experience
n 大学期间任班级团支书。多次获“校优秀团员”“校优秀团干”等荣誉称号;
n 曾担任校园广播台节目部主持人,策划并主持《校园话题》、《三味书屋》等广播栏目,熟练运用Adobe Audition软件制作宣传音频;
n 多次策划并执行“广播台狂欢夜”“舞蹈大赛”“我爱记歌词”等大型校园晚会活动;
n 大二升任校园广播台网络部部长,曾管理蓝V微博账号 @北理珠校园广播台 一年
þ教育背景
Education
2008.9-2012.7 实习僧大学 市场营销(本科)
主修课程
基本会计、统计学、市场营销、国际市场营销、市场调查与预测、商业心理学、广告学、公共关系学、货币银行学、经济法、国际贸易、大学英语、经济数学、计算机应用等.
校内荣誉
þ 2009.10获国家奖学金
þ 2010.11获“三好学生称号”
þ 2010.12华南大学生创意营销大赛一等奖
þ 2011.4挑战杯创业计划大赛省级铜奖
展开阅读全文