基于虚拟激励法的多激励振动试验数值分析.pdf
《基于虚拟激励法的多激励振动试验数值分析.pdf》由会员分享,可在线阅读,更多相关《基于虚拟激励法的多激励振动试验数值分析.pdf(8页珍藏版)》请在咨信网上搜索。
1、第 6 卷 第 2 期2023 年 6 月空天防御AIR&SPACE DEFENSEVol.6,No.2Jun.,2023基于虚拟激励法的多激励振动试验数值分析王肇喜,翟师慧,赵凡,王者蓝,谢夏阳(上海航天精密机械研究所,上海201600)摘要:针对细长型飞行器双台振动试验开展仿真计算与控制效果分析。根据模态叠加原理,建立细长型飞行器多输入多输出(multi-input multi-output,MIMO)随机振动试验中的动力学模型;采用虚拟激励法和逆虚拟激励法计算方法,分析细长型飞行器MIMO随机振动试验中不同控制测点数量、控制谱矩阵设置的可实现性及影响效果,并以其模拟件为研究对象,完成其两
2、点、多点控制试验方式下不同控制谱矩阵设置的振动试验的数值模拟。经分析,虚拟激励法大幅提升了振动仿真计算精度,且双台振动系统在不同控制点组合下振动响应差异性加大。关键词:振动试验;虚拟激励法;多输入多输出(MIMO);数值模拟中图分类号:TJ761 文献标志码:A 文章编号:2096-4641(2023)02-0069-08Numerical Simulation of the Multi-Input Multi-Output Random Vibration Tests Based on Pseudo Excitation MethodWANG Zhaoxi,ZHAI Shihui,ZHAO
3、Fan,WANG Zhelan,XIE Xiayang(Shanghai Spaceflight Precision Machinery Institute,Shanghai 201600,China)Abstract:In this paper,the numerical simulation and control effect of double-vibration test for the slender aerocraft were studied.Using the Mode-superposition method,a dynamic model of the slender a
4、erocraft under multi-input multi-output random vibration tests was proposed in this paper.The pseudo excitation method and inverse pseudo excitation method were adopted to study the influence of the reference point number and the reference spectral density matrix.The numerical simulation of the Mult
5、i-Input Multi-Output random vibration tests of a simulator of the slender aerocraft under different reference point numbers and reference spectral density matrices was acquired.Finally,the numerical simulation results were promoted through the pseudo excitation method,and the vibration responses of
6、the various control point combinations were different for the double-vibration system.Keywords:vibration test;pseudo excitation method;MIMO(multi-input multi-output);numerical simulation0引言振动环境试验是目前航天产品可靠性和环境适应性最主要的检验手段,用于发现产品设计上的缺陷、工艺上的不成熟及可靠性的薄弱环节,为产品设计及工艺上的改进与优化指明方向1。细长型飞行器的振动环境试验为工程中的一大难题,由于此类飞行
7、器外形细长,振动试验中其超出振动台面的悬臂部位过多,振动响应悬臂端低频放大、高频衰减,使试验室状态下飞行器的振动环境较之飞行状态有较大失真,飞行器各部位的过试验、欠试验现象异常严重。多输入多输出(multi-input multi-output,MIMO)振动试验技术采用两(多)个振动台提供激振力,使试件不同部位产生不同的振动响应。这种试验方式能够更加准确地模拟细长型飞行器各种飞行状态下的振动环境,有效地避免传统单振动台试验方法中振动应力集中的现象,为力学环境试验技术的重大进步和收稿日期:2022-12-14;修订日期:2023-04-21作者简介:王肇喜(1986),男,硕士,高级工程师,主
8、要研究方向为振动、冲击试验理论与技术研究。空天防御第 6 卷未来发展趋势1。目前,国内少数航空、航天单位已经开始应用MIMO振动试验方法。但由于我国现有的试验标准对该类型试验的描述较为笼统,并且国内对互谱的作用、双振动台模拟互谱的可实现性及其效果的研究不够透彻1-2,加之此项技术的军事背景,使得国外相关的技术资料难以获得,因此,我国工程人员在进行MIMO振动试验时,面对控制参数的选取和试验条件的制定等问题,大多凭借工程经验,采用试验中不断变换、试探、摸索的方法,这将会走很多弯路,造成大量人力、物力、财力的浪费,而最终的试验效果也难以达到最优。本文基于虚拟激励法、逆虚拟激励法原理3-5,研究MI
9、MO随机振动试验中控制谱矩阵设置的可实现性,并对某一细长型飞行器模拟件在不同控制点数量、控制参数设置情况下的振动试验进行数值模拟,分析控制点数量及控制互功率谱对双台振动试验的影响。应用此方法,在正式试验前对各种试验方案进行预试验,使试验人员提前了解振动试验中不同控制点位置、参数设置所产生的效果,提前知晓试件在试验过程中各部位振动响应情况,提前评判所选用的试验仪器设备能否顺利完成试验,避免在试验条件定制过程中的盲目性,由此可降低大型复杂试件振动试验的风险,获取最优的试验设计方案。1随机振动问题的分析方法1.1随机振动基本理论对于多输入多输出振动系统,其输入、输出谱密度矩阵为5-6Sxx()=Sx
10、1x1()Sx2x1()Sxnx1()Sx1x2()Sx2x2()Sxnx2()Sx1xn()Sx2xn()Sxnxn()(1)式中:为圆频率;n为系统的输入或输出个数。对角线元素为各点的自功率谱,它表述随机振动在频域的能量分布,非对角线元素为两点间的互功率谱,它含有相干系数和相位这两个参数的信息,表达式为Sxixj=SxixiSxjxjxixjejxixj(2)式中:xixj为相干系数;xixj为相位差;i,j为公式索引变量。激励与响应的功率谱密度矩阵的关系式为6-7Syy()=-H()Sxx()HT()(3)式中:-H()为传递函数矩阵H()的共轭;HT()为传递函数矩阵的转置;Sxx()
11、为系统的输出功率谱密度矩阵;Syy()为系统的输入功率谱密度矩阵。1.2虚拟激励法理论直接使用式(3)计算MIMO随机振动问题,计算量太大,因此林家浩等8提出虚拟激励法,这种方法在不降低计算精度的前提下,大大降低了求解该类问题的计算量。虚拟激励法3-4基本手段为矩阵的对角化变换,MIMO 系统的输入、输出功率谱密度矩阵Sxx()为Hermite矩阵,根据矩阵理论,Hermite矩阵一定酉相似于对角形矩阵,且其特征值为实数。则功率谱密度矩阵可转化为Sxx()=U1()U2()Un()1()2()k()00nnU*1()U*2()U*n()(4)式中:1()k()为矩阵Sxx()的特征值;k为矩阵
12、Sxx()的秩;U1()Un()为酉矩阵的列向量;Ui*()为Ui()的共轭转置。式(4)可转化为Sxx()=i=1kUi()i()U*i()(5)令-xi()=i()Ui(),则Sxx()=i=1k-xi()xTi()(6)1)各激振力相互之间完全相干时,Sxx()的秩k=1,式(6)可改写为Sxx()=-x()xT()(7)70第 2 期王肇喜,等:基于虚拟激励法的多激励振动试验数值分析2)各激振力相互之间完全不相干时,Sxx()为对角线矩阵,式(6)可改写为Sxx()=-x1()xT1()+-x2()xT2()+-xn()xTn()(8)3)激振力部分相干时,Sxx()的秩为k(kn),
13、式(6)可改写为Sxx()=-x1()xT1()+-x2()xT2()+-xk()xTk()(9)由于xi()为列向量,且相互正交,则可将MIMO随机振动系统的输入项视为几个相互正交的虚拟简谐激励力xi()eit叠加的效果。将式(6)代入式(3)得:Syy()=-H()i=1k-xi()xTi()HT()=i=1k-H()-xi()H()xi()T(10)则Syy()=i=1k-yi()yTi()(11)其中yi()=H()xi()(12)同理,根据式(11),可将MIMO系统的振动响应视为虚拟简谐激励xi()eit产生的虚拟简谐响应yi()eit叠加的效果,其虚拟激励与虚拟响应的关系见式(1
14、2)。应用虚拟激励法,求解 MIMO 系统随机振动响应的基本流程为:将输入的谱密度矩阵分解为式(6)的形式,构造虚拟激励xi();由式(12)得出虚拟响应yi();由式(11)输出功率谱密度矩阵。1.3逆虚拟激励法理论逆虚拟激励法5多用于随机振动载荷识别问题领域,即已知结构的随机振动响应和结构的固有特性(频响函数矩阵),反推结构所受的激振力。逆虚拟激励法的计算流程为:将Syy()改写为式(11)的形式,构造虚拟响应yi();由式(12)构造方程,求解xi();由式(6)得出激振力的功率谱密度矩阵。2双振动台试验模型分析两个振动台组成的双台振动试验系统如图 1 所示,其中加速度传感器、电荷放大器
15、、计算机、功率放大器、振动台构成一个闭环反馈控制系统。试验过程中闭环控制系统不断反馈修正,调整两振动台输出的激振力,使试件上控制测点的振动响应功率谱与预设值相等。不同厂家的振动控制系统,其控制器的控制迭代算法各有不同,但经反馈修正,最终稳定时,两个振动台输出的激振力都使控制点实际振动响应与预设值相等或者偏差最小。根据此原理,假设控制器性能理想,对细长型飞行器双振动台试验建立力学模型,并进行分析。图1双台振动试验系统Fig.1The system of double vibration test 71空天防御第 6 卷2.1双振动台试验分析模型的建立由于球头(见图1)的作用,将试件与振动台连接处
16、的转动自由度释放,且弹簧绳的拉力与试件的重力相抵消,则试件的边界条件可视为自由状态,其受力模型如图2所示,图中fi、fj分别为i、j位置点所受的力。该受力模型的动力学方程为9-10Mx +Cx +Kx=F(13)式中:M、C、K、F分别为飞行器的质量矩阵、阻尼矩阵、刚度矩阵、所受的力矩阵。式(13)经傅里叶变换表示为(-2M+jC+K)X()=F()(14)即X()=H()F()(15)其中,F()为F()=00fi()00fj()00T(16)式中:fi()、fj()分别为图 2 中 i、j 位置点所受的力谱。结构的整个传递函数矩阵,可由模态测试的结果,根据模态叠加法11-12近似得到:H(
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 虚拟 激励 振动 试验 数值 分析
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。