附加自由阻尼板高频振动响应的能量流模型.pdf
《附加自由阻尼板高频振动响应的能量流模型.pdf》由会员分享,可在线阅读,更多相关《附加自由阻尼板高频振动响应的能量流模型.pdf(10页珍藏版)》请在文库网上搜索。
1、振动与冲击第42 卷第14期JOURNAL OFVIBRATIONAND SHOCKVol.42 No.14 2023附加自由阻尼板高频振动响应的能量流模型江旭东,韩月东,滕晓艳(1.哈尔滨理工大学机械动力工程学院,哈尔滨150 0 8 0;2.哈尔滨工程大学机电工程学院,哈尔滨150 0 0 1)摘要:为了预示自由阻尼结构(free layer damping,FLD)的高频振动响应,将复刚度法与能量流分析(energy flowanalysis,EFA)理论相结合,推导了大阻尼条件下高频振动自由阻尼板的能量密度方程。分析了阻尼交界面处的能量传递关系,利用耦合结构的能量流分析方法推导了局部附
2、加阻尼处理的薄板结构的能量密度方程。基于上述能量密度方程,构建了完全/局部自由阻尼板的能量有限元法(energyfinite elementmethod,EFEM)模型,求解了阻尼层合结构的高频能量流响应。通过对比模态解析解与能量有限元模型的数值解,验证了所建立的自由阻尼板高频振动能量流模型的有效性。关键词:能量流分析;自由阻尼结构;能量有限元法(EFEM);高频振动中图分类号:TB123Energy flow model for a high-frequency vibrating plate with free-layer damping treatment(1.School of Mec
3、hanical and Power Engineering,Harbin Science and Technology University,Harbin 150080,China;2.College of Mechanical and Electrical Engineering,Harbin Engineering University,Harbin 150001,China)Abstract:To predict the high-frequency vibrational responses of a plate with free layer damping(FLD)treatmen
4、t,the energy density equation considering high damping effect was derived by using the energy flow analysis(EFA)methodcombined with the complex stiffness method.The wave transmission analysis was performed for clarifying the powertransfer relation at discontinuous interfaces of the FLD.Then the ener
5、gy density equation of a plate with partial FLDtreatment was also established based on the EFA of coupled plates.The energy finite element method(EFEM)wasemployed to solve the derived energy density equations of the plates with full or partial FLD treatment.The simulatedresults obtained by the EFEM
6、were compared with those obtained by the modal superposition method.The resultsdemonstrate that the proposed EFA model can predict the energy density of the high-frequency vibrating plate with FLDtreatment with reasonable accuracy.Key words:energy flow analysis;free damping structure;energy finite e
7、lement method;high-frequency vibration文献标志码:AJIANG Xudong,HAN Yuedong,TENG Xiaoyan?2D0I:10.13465/ki.jvs.2023.14.038高速飞行器通常工作于高频、高声压级的极端噪声环境中,诱发的结构振动将影响舱室内精密仪器的正常运行,甚至造成关键电子设备的损伤或失效1-3。阻尼减振技术通过敷设于大尺寸薄壁构件的黏弹性材料层,可以实现振动与噪声的有效控制,因而应用于高速飞行器舱壁结构的减振降噪中4-5。由此,在概念设计阶段对附加阻尼薄壁结构动力学特性的预测分析,实现基于结构声振性能的优化设计,具有重要的
8、理论与工程价值。在空天结构领域,板壳结构广泛应用于高速飞行基金项目:国家自然科学基金(52 17 550 2;5150 50 96);黑龙江省自然科学基金联合引导项目(LH2020E064)收稿日期:2 0 2 2-0 5-11修改稿收到日期:2 0 2 2-0 9-2 1第一作者江旭东男,博士,副教授,197 7 年生器的机身和舱室等轻质薄壁结构的动力学建模分析。夏小均等6 基于波函数法理论,推导了自由阻尼薄板振动分析模型以及包含自由阻尼结构与声腔的三维耦合模型的建模方法。Liu 等7 考虑了边界的弹性与阻尼特性,研究了任意边界条件下的薄板弯曲振动问题。高速飞行器一般工作于气动流和外界噪声作
9、用的高频声振环境,力学行为上表现为短波高频振动特性。上述以位移为基本变量的确定性分析方法由于计算量以及参数敏感性等原因难以适用于高频动力学分析8 空天结构的高频振动分析一般采用以子系统平均声振能量为基本变量的统计能量法(statistical energymethod,SEM)。Pa v a n 等9分析了湍流激励下飞机舱室内的高频声场分布。Raef等10 预测了锥壳封闭声腔的声传递损失。Poblet-Puigl将谱单元方法与统计能量法相结合,建立了L型折板的传递损失分析模型;322陈强等12-13研究了热环境下长宽比对L型折板统计能量分析参数的影响。胡婉璐等14耦合梁结构间的耦合损耗因子。G
10、ao 等15基于FE-SEA 混合建模技术,提出了板梁耦合结构中频振动问题的设计优化方法。张婧等16 基于瞬态SEA方法,预测了线耦合平板的高频振动响应。但是,统计能量法仅能给出子系统的平均振动响应,无法预测子系统内部任意空间位置的精确响应,因而统计能量法难以精准地预示附加阻尼结构的高频振动响应。与SEA不同,能量流分析方法以能量密度为基本变量,可以精确地预测振动能量的空间分布。Zhang等17 和Wang等18 考虑了热应力和材料参数的温度相关性,建立了均匀温度场下梁和板结构的能量流模型。上述研究19忽略了热应力对势能密度的影响,无法精确预示预应力结构的高频振动响应。Yan20首次利用能量有
11、限元法(energy finite element method,振动与冲击1完全自由阻尼薄板的能量流模型1.1自由阻尼板的等效处理如图1所示为简谐激励作用下的自由阻尼板,板长和板宽分别为l和l阻尼层厚度为h2,基板厚度为hi简谐点激励作用于板(xo,y o)处。依据Kirchhoff 经典薄板理论,横向弯曲振动自由阻尼板在,y,z方向的位移场表示为u(x,y,z)=-(z-zo)awo(x,y)/axv(x,y,z)=-(z-zo)awo(x,y)/ayw(x,y)=wo(x,y)式中,z。为物理中面的坐标。自由阻尼层、22023年第42 卷(1)(2)(3)h2FoeJo(Xo.o)hlE
12、FEM)研究了复合材料层合板结构的高频振动数值预报方法。Liu等2 1-2 建立了功能梯度梁结构的能量流模型,分析了热梯度环境下功能梯度梁的高频振动响应。Kim等2 3基于等效流体模型,提出了压缩波在高频振动多孔介质结构中传播的能量流分析方法。Cho24研究了应力波在结构连接处的波型转化与传递特性,建立了组合结构的能量有限元模型。Liu等2 5利用耦合损耗因子描述混响波场在耦合结构间的功率耦合关系,提出了改进的耦合结构问题的能量有限元法。Park26构建了耦合Mindlin板面内、面外振动的能量流模型。Liu等2 7 融合能量有限元法和水平集法,提出了板梁耦合结构的低噪声拓扑优化方法。Yeo等
13、2 8 建立了压力波在高频振动三维弹性体中传播的能量流模型。Zeng等2 9预测了机械激励与轨道噪声激励联合作用下的高速列车车内噪声水平。陈兆林等30 和Zhao等31混合能量有限元法和虚拟模态综合法,提出了高频冲击载荷作用下结构瞬态响应预报方法。Chen等32 研究了超声速气流环境下薄板结构的高频振动问题。上述研究成果均基于结构的小阻尼假设,在能量密度控制方程的推导中对波数实施了一阶近似估计,对于预测具有大阻尼损耗因子结构的高频振动响应将产生显著的误差,Han 等33-35已通过一维大阻尼系统的能量流分析验证了上述观点。自由阻尼板结构属于具有大阻尼特征的层合板结构,基于小阻尼假设的能量流分析
14、方法不再适用于高频振动响应的预测。本文将复刚度理论与能量流分析理论相结合,推导大阻尼层合板结构的能量密度方程,建立完全/局部自由阻尼板的能量有限元模型。以模态解析解作为基准解,验证所提出的理论模型对于预测附加自由阻尼结构高频振动响应的有效性。XZ3基板图1简谐激励作用的自由阻尼板Fig.1A plate with full free layer damping treatmentunder a concentrated harmonic load由式(1)式(3),自由阻尼板的应变场表示为8=(z-Zo)K式中,=kKKyT 为曲率矩阵,其中 K=-w/ax*,k,=-a*w/ay2,y=-2
15、aw/axay。假设自由阻尼板中基板和阻尼层的应力矩阵为g=o 0(h=1为基板,h=2为自由阻尼层),则小变形条件下结构的应力应变的关系为d=Q10Q=三10,三00(1-vk)/2式中:E,分别为各层的弹性模量和泊松比;Q为各层的弹性矩阵。自由阻尼板的弯曲应变能等于基板和附加自由阻尼层的弯曲应变能之和,即eD.=22三(z-zo)dzzk+12zkD(1-v)三(z-zo)dzeqvJzk+1式中:De和Deqy分别为自由阻尼板的等效弯曲刚度和等效扭曲刚度;h,和h分别为基板和阻尼层的厚度。层合阻尼板的物理中面坐标可以由x轴方向的力物理中面E2(6)2xydS(7)(8)z(9)(4)(5
16、)第14期平衡得到三2(z-zo)dz+三l由此,自由阻尼板中面的位置zo为三?h2+2三 hhi+三 hiZo=2三2 h2+2三hl为了反映自由阻尼板的刚度的阻尼特性,将式(8)的等效复刚度表示为D=+三*(1+jmi)(z-zo)dzk=1式中,为各层板的材料阻尼因子。根据式(12),自由阻尼板的等效阻尼损耗因子为2Im(De)meq=Re(Dt)另外,依据式(7),自由阻尼板的等效泊松比为U=1Deq1.2能量密度控制方程通过哈密尔顿原理,自由阻尼板的弯曲振动控制方程为DaVw+m.aw=F。(x -x o,y -y o)e j o l (15)at?式中:me=pih,+pzhz(p
17、1,p2分别为基板和自由阻尼层的密度)为自由阻尼板的等效表面质量;F。,=2 f(f 为激振频率)为简谐激励的幅值和激振圆频率,(-o)(-x)为二维狄拉克函数。自由阻尼板振动控制方程式(15)的稳态振动解为W(x,y,t)=Ae-i(k+ky),ejol(16)式中,k,k,分别为板在x,y方向上的波数分量,其中k=kicos,ky=krsin,kr为弯曲总波数,为弯曲波的散射角。将式(16)代入控制方程式(15)中,得到结构的频散方程(17)由式(17),弯曲波数表示为44m.medeq0kf=Deq根据波动理论,式(18)中的前两项对应于空间振荡的远场波,后两项对应于指数衰减的近场波。在
18、高频振动状态下,弯曲波的波长远小于板的面内尺寸,因而能量流分析中一般忽略近场波带来的影响。在小阻尼假设下(neg1),弯曲波数的一阶近似为4meqkr=Re(De)在混响波场假设下,弯曲波在各个可能的方向上江旭东等:附加自由阻尼板高频振动响应的能量流模型(z-zo)dz=0(10)1krcos odp(11)12TJkrsin do根据式(18)、式(2 0)和式(2 1),式(16)重新表示为Z(12)wr(x,y,t)=(Afe-ix-iy+Bejkx-jk,y+z:+1*=(z-z0)dz2=(z-zo)dz+k=12Deqy2=022323传播,因而总波数、波数分量在散射角0,2 内需
19、要均值化处理,则有(20)(21)(22)式中,Ar,Br,Cr,D,为各远场行波的幅值。(13)另外,一般金属结构的阻尼因子均较小,因而现有的能量流理论基于小阻尼假设,通过波数对于阻尼因子的一阶近似,推导结构的能量密度方程。然而,黏弹阻尼层的敷设导致附加阻尼结构的损耗因子一般远大(14)于基板的阻尼,对波数进行小阻尼近似将带来显著的系统误差。因此,为了精确预示自由阻尼板结构的高频振动响应,应该对现有的传统能量流理论进行大阻尼修正,通过精确的波数表达式(18)进行结构能量密度控制方程的推导。薄板结构在一个振动周期内时域平均的能量密度e)由动能能量密度 ek和弯曲变形势能能量密度组成。e)=ek
20、)+ep)薄板结构时域平均的动能能量密度、弯曲变形势能能量密度分别表示为mee4dWaw2+2axdw(w)dw+2(1y式中:()*为共轭运算;Re()为复数取实部函数。类似于与能量密度,时域平均的能量强度 I.和2I,分别表示为(18)1.)=ReD-Q2根据薄板弯曲振动理论,可以得到弯矩M,M,,扭矩My,剪力Qx,Q y。将式(2 2)给出的远场解分别代人式(2 3)式(2 7),经过一个波长范围内的局部空间平均以消除行波之间干涉项引起的功率流,可以获得时neq(19)(23)w(24)att2ved.x()axdy)(25)ww+M.+Mdxdtdyt+Mw+Mw1(27)yxdyd
21、tdxdt间空间平均处理后的能量密度(J/m)与能量强度(J/s)。W+2(26)324振动与冲击ReiDaRit+2uhi,+hgi+=2mmih+2h2,+kgkk.22023年第42 卷2(1-ag)k,k,k k,/入2(1-g)h,k;h,+mg).+)(2 8)采用无限导纳计算平均输人功率元),则有12wRelDi+kig+his+vegkgh:+(1-ve)h,k,k,J1-+-+)(29)1wRelDnn+hh,+vaii,+2hgh,+(1-e)k,k,k1.-+)式中:)=A)-*B)+*C)-+D)+,其中系数 A)*满足式(31),其他系数类同。k.k,A)=x2Fy2
22、式中,kx,,和k,分别为k和k,的虚部。对比式(2 8)与式(2 9)和式(30),描述能量密度与能量强度关系的能量传递方程为ae)xoL+h+Rk+3k:+2(1-0(33)B:ky,;入入=k,1+20ghg+hg+2(1-0g)hkhky+meqa/De(35)根据能量守恒定律,输入功率一部分以能量强度梯度的形式存储于结构内部,另一部分被耗散掉,因此功率平衡方程可以表示为=+V.i)式中:V.i)=a /a x+/a y 为函数的散度运算;元),元s)分别为平均输入与耗散功率。薄板结构的平均耗散功率与频率、结构阻尼因子和势能密度相关,即能量耗散方程为=2m0(e,式中:n为结构损耗因子
23、;e为结构总变形势能密度。在小阻尼假设下,势能密度与动能密度相等,能量耗散方程为变为=no(e)将式(2 5)和式(2 8)代人式(37),自由阻尼板的平均耗散功率表示为=e)(40)元)=2Re8meq对于尺寸有限的结构,其确切的阻抗取决于边界条件和加载的位置,因此很难直接计算,除非通过有限元法、解析法或其他传统方法求解运动方程。然而,具有(30)无限维度的结构的阻抗在高频下非常接近有限结构的阻抗3-3,因此,在EFA和EFEM中,总是应用无限结构的阻抗来逼近有限结构的精确阻抗。由此,将式(31)、式(37)和式(39)代人式(35),自由阻尼板的能量密度控制方程为e)e+(e)=8(-0,
24、y-yo)(42)dy式(42)为大阻尼条件下附加自由阻尼版高频振动的能量密度控制方程,与传统的薄板结构的能量密度控制方程相比,考虑了附加阻尼层对薄板振动的影响,采用了精确的波数表达式,无需小阻尼假设,因此能更3e)ay(41)(32)好的预示附加阻尼板结构的高频振动响应。采用有限元法求解自由阻尼板的能量密度控制方程,则有Ke=Fe+Q2-2-1-22K=-ab6a-111-1-221-1-21Ba166-1-22一2-112(36)Fe式中:K为单元刚度矩阵;F为单元输入功率矩阵;N为4节点双线性四边形单元的形函数矩阵;Q为单元边界能量通量矩阵;e为单元节点能量密度矩阵;,b(37)为单元边
25、长尺寸的一半。2局部自由阻尼薄板的能量流模型工程实际中,附加阻尼层往往敷设于振动剧烈区(38)域。因此,本章构建局部自由阻尼板的能量流响应模型。如图2 所示为局部自由阻尼板结构,沿阻尼边界将整个自由阻尼板分成子域和子域。当弯曲波人(39)射于附加阻尼层交界面处时,不同的波数将引起弯曲(43)171-12-222-2-1yab2421(44)+1912422124-Nin)o(x-xo,y-yo)do(45)2eQ=f.NTqdre212(46)(56)第14期波的反射和透射现象。因此,可以将局部自由阻尼板等效为具有不同物理参数的耦合板结构。自由阻尼层基板IFoeio(xoyo)X图2 简谐激励
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 附加 自由 阻尼 高频 振动 响应 能量 模型