考虑层间变形的约束阻尼圆柱壳振动特性分析.pdf
《考虑层间变形的约束阻尼圆柱壳振动特性分析.pdf》由会员分享,可在线阅读,更多相关《考虑层间变形的约束阻尼圆柱壳振动特性分析.pdf(9页珍藏版)》请在咨信网上搜索。
1、第 43 卷第 5 期2023 年 10 月振动、测试与诊断Vol.43 No.5Oct.2023Journal of Vibration,Measurement&Diagnosis考虑层间变形的约束阻尼圆柱壳振动特性分析马宏伟1,2,陈中石1,2,孙伟1,2(1.东北大学机械工程与自动化学院 沈阳,110819)(2.东北大学航空动力装备振动及控制教育部重点实验室 沈阳,110819)摘要 在约束阻尼减振结构中,中间的黏弹性层主要依靠基体及约束层的交错运动产生剪切变形而消耗能量,因而这种层间变形关系需要被有效引入以提升建模精度,为此采用有限元法对局部贴敷约束阻尼层圆柱壳进行了动力学建模方法研
2、究。首先,基于 Donnells薄壳理论,在考虑各层之间变形关系的基础上,推导了复合单元的刚度及质量矩阵;其次,构造了一种过渡单元来连接复合壳单元和仅有基体的单层单元,实现了局部贴敷约束阻尼层圆柱壳总刚度及总质量矩阵的有效组集;然后,在考虑复合圆柱壳弹性约束边界条件的基础上,确定了复合壳的运动方程,并给出了求解振动特性的公式;最后,进行了实例研究,通过与实验和 ANSYS仿真结果的对比,验证了所提有限元模型的正确性,并进一步分析了约束层、阻尼层厚度变化对复合圆柱壳振动和阻尼特性的影响。关键词 层间变形关系;局部贴敷;约束阻尼层;圆柱壳;振动特性;有限元建模中图分类号 TH113.1;V214.
3、3;TB535引 言圆柱壳结构由于比强度、比刚度高,在航空、航天及舰船等装备中有着广泛的应用。其在机械、气动及热等载荷作用下,很容易引发过大的振动而导致结构本身或者整个装备的性能失效。众多减振方法中,在圆柱壳结构的表面贴敷约束阻尼层(通常为局部贴敷)是一种最经济、最有效的方法1,2。针对圆柱壳结构进行局部贴敷约束阻尼层减振,首先要创建可以预测复合壳结构振动特性的动力学分析模型。当前,相关研究方法主要包括解析法(或者半解析法)和有限元法。Chen 等3利用薄壳 理 论 和 唐 奈 尔穆 斯 塔 里弗 拉 索 夫(DonnellMushtari Vlasov,简称 DMV)的假设推导了局部贴敷约束
4、阻尼圆柱壳结构的一般运动微分方程,分析了贴敷厚度对复合圆柱壳振动特性的影响,这是一种典型的半解析法。基于 DMV 假设和拉格朗日方程,Zheng等4推导了局部贴敷约束阻尼层圆柱壳结构的动力学方程,分析了不同贴敷方案的约束阻尼层圆柱壳结构振动特性。基于线弹性薄壳理论,向宇等5推导了局部贴敷约束阻尼层圆柱壳的一阶常微分方程,并分析了贴敷位置对复合壳振动特性的影响。基于 Donnells 假设和线性黏弹性理论,Zheng等6构造了被动约束层阻尼圆柱壳的本构关系,利用哈密顿原理导出了约束层阻尼处理圆柱壳的运动方程。损耗因子是表征贴敷阻尼材料结构损耗能量的重要指标,张安付等7提供了一种多种阻尼处理下圆柱
5、壳损耗参数的解析换算方法。相对于解析或者半解析法,有限元法具有更广泛的适用性。Wang 等8利用有限元法将局部贴敷约束阻尼层圆柱壳等效为单层后,进行了振动和阻尼分析,并讨论了材料性能、厚度、处理区域和位置等参数对复合壳振动特性的影响。Masti 等9基于有限元法和一阶剪切变形理论,建立了一种 4 节点80自由度的复合壳单元,用于分析局部贴敷约束阻尼层圆柱壳结构的振动特性,并与实验对比验证了模型的正确性。经过约束阻尼层处理的圆柱壳部分主要由基体层、黏弹性层和约束层构成,其主要减振原理是在圆柱壳振动时,黏弹性层因基体层和约束层的交错变形出现剪切效应,导致消耗更多的能量,所以在复合圆柱壳动力学建模过
6、程中,需要考虑这种层间的变形关系。基于层间变形关系和改进的波传递法,Gao等10采用瑞利里兹法求解了贴敷约束阻尼层圆柱壳的自由振动问题。基于半解析法,Song等11在考虑层间变形的基础上,推导了约束阻尼层圆柱壳在任意边界条件下的动力学方程。基于DOI:10.16450/ki.issn.10046801.2023.05.026 中央高校基本科研业务费资助项目(N180312012)收稿日期:20210317;修回日期:20210425第 5 期马宏伟,等:考虑层间变形的约束阻尼圆柱壳振动特性分析薄壳理论,Mokhtari等12在考虑层间变形关系下,利用拉格朗日能量方程和瑞丽里兹法求解了复合圆柱壳
7、的动力学方程,并分析了厚度和半径等参数对结构振动特性的影响。利用层间变形关系和经典薄壳理论,Jin 等13采用傅里叶级数作为容许函数,推导了约束阻尼圆柱壳的动力学方程,并讨论了不同层厚和剪切参数对其振动特性的影响。上述考虑层间变形关系的约束阻尼层圆柱壳动力学建模均属于解析或者半解析法,如果这种层间变形关系没有引入到有限元建模分析中,将导致有限元分析结果与实际结果存在偏差。对于局部贴敷约束阻尼层的圆柱壳结构,既存在包含 3 层结构的复合壳单元,也存在仅包含基体的单层壳单元。在有限元法建模中,如何将这种不同类型的单元组集在一起,实现正确的单元组集非常重要。笔者采用有限元法对局部贴敷约束阻尼层圆柱壳
8、结构进行了动力学建模方法研究。在建模中考虑了层间变形关系,推导了复合单元的刚度及质量矩阵,构造了一种过渡单元来实现局部贴敷约束阻尼层圆柱壳总刚度及总质量矩阵的有效组集。1 考虑层间关系复合圆柱壳单元推导本研究的局部贴敷约束阻尼层圆柱壳模型如图 1所示,底部为弹性约束,并在自由端贴敷约束阻尼层。其中:R为复合圆柱壳的等效半径;L为圆柱壳总高度;L1为约束阻尼层贴敷的宽度;hs,hv和hc分别为基体、黏弹性层和约束层厚度。在圆柱壳中心建立柱坐标系(x,z),对应的复合圆柱壳的等效中性面上局部位移用u,v和w来表示。中性面到圆柱壳内表面的距离h014可定义为h0=12k=13E(k)(H2k-H2k
9、-1)k=13E(k)h2k (k=1,2,3)(1)其中:E(k),Hk分别为第k层材料的弹性模量和上表面到圆柱壳内表面的距离,且H0=0;k=1,2,3分别表示基体层、黏弹性层和约束层1.1层间变形关系基于Donnells薄壳理论,复合圆柱壳中基体层和约束层内任意一点位移场可表示为 uk()x,z=u0k()x,+zkx()x,vk()x,z=v0k()x,+zk()x,w()x,z=w0()x,(2)其中:u0k,v0k和w0为第k层中性面各方向位移;kx,k为各层中性面在x和方向的转角。约束阻尼层圆柱壳层间变形如图 2 所示,假设复合结构基体层和约束层中性面分别由AA 处沿周向变形了v
10、s和vc,变形至BB 处。这时,通过BB 处的CC 和DD 层面可以看出,各层间位移变化是协调连续的11。以周向变形为例,黏弹性层的位移uv和vv及转角位移vx和v可表示为 uv=12 us+uc+hs2sx-hc2cx vv=12 vs+vc+hs2s-hc2cvx=1hv uc-us-hs2sx-hc2cxv=1hv vc-vs-hs2s-hc2c(3)1.2复合单元推导复合圆柱壳单元如图 3 所示,其中Rs,Rv和Rc分别为基体、黏弹性层和约束层的几何中面半径。图 1局部贴敷约束阻尼层圆柱壳模型Fig.1Cylindrical shell model with partially att
11、ached constrained damping layer图 2约束阻尼层圆柱壳层间变形示意图Fig.2Diagram of interlayer deformation of cylindrical shells with constrained damping layers图 3复合圆柱壳单元Fig.3Shell element of composite cylindrical shell1019振 动、测试与诊断第 43 卷 该复合壳单元包含 4 个节点,设每个节点具有15个自由度,任一节点位移qi为qi=uis uisx uis vis visx vis uic uicx uicv
12、ic vicx vic wi wix wiT(4)其中:i=1,2,3,4。由层间关系描述可知,黏弹性层的位移变形可由基体层和约束层的位移变形来表达。在单元推导过程中,可以先确定基体层、约束层的位移和应变等,再引入层间变形关系确定黏弹性层的上述参数。可以假定任意方向位移,例如 z 方向位移w在曲面坐标x和方向上的位移场为w=a1+a2x+a3+a4x2+a5x+a62+a7x3+a8x2+a9x2+a103+a11x3+a12x3 (5)其中:a1,a2,a12为多项式系数。复合壳单元等效中性面上任意一点的位移场为usvsucvcw=i=14Niqi(6)其中:Ni为节点i=1,2,3,4对应
13、的形函数。Ni=NiNixNi000000000000000NiNixNi000000000000000NiNixNi000000000000000NiNixNi000000000000000NiNixNi(7)为方便构造黏弹性层的形函数,式(7)中形函数Ni也可改写为Ni=Nusi Nvsi Nuci Nvci Nwi(8)式(8)对应基体及约束层的形函数。将式(6)代入描述层间变形关系的式(3)中,并进行相应的整理,得到黏弹性层的形函数Nuvi,Nvvi,Nvxi和Nvi。复合壳单元某层任意一点应变与位移的关系为kx=0(k)x+zkx=u0kx-z 2w x2k=0(k)+zk=1Rkv
14、0k+wRk+z(1R2kv0k-1R2k2w2)kx=0(k)x+zkx=1Rkuk+vkx+z(1Rkvkx-2Rk2wx)(9)其中:kx,k和kx分别为基体、约束层的任意一点的轴向、切向及剪应变;0(k)x,0(k)和0(k)x分别为各层几何中性面上任意位置的轴向、切向及剪应变;kx,k和kx分别为各层几何中性面上任意位置轴向、切向及法向弯应变。整理可得复合壳单元各层中性面应变位移关系为0(k)x 0(k)0(k)x kx k kx=LNkqe=Bkqe (10)其 中:L为 矩 阵 微 分 算 子;各 层 形 函 数Ns=dig(Nus,Nvs,0,0,Nw),Nv=dig(Nuv,
15、Nvv,Nw),Nc=dig(0,0,Nuc,Nvc,Nw);qe为单元节点位移向量;Bk为第k层应变矩阵。黏弹性层的剪切应变为vxz=1hv(uc-us)+12hv(hs+hc+2hv)wxvz=(1hv-hc2hvhc-12Rv+hc4RvRc)vc-(1hv+hs2hvhs+12Rv+hs4RvRs)vs+(hs2hvRs+hc2hvRc+hs4RvRs-hc4RvRc+1Rv)w(11)同理可得黏弹性层的剪切应变位移关系为vxz vz=BvSqe(12)其中:BvS为黏弹性层的剪应变矩阵。黏弹性层在x和方向的剪切应力可表示为vxzvz=Qv4400Qv55 vxzvz(13)其中:为剪
16、切修正因子,一般取=5 6。Qv44=Qv55=Gv=G(1+iv)(14)其中:Gv为黏弹性材料的复剪切模量;v为黏弹性材料损耗因子。1020第 5 期马宏伟,等:考虑层间变形的约束阻尼圆柱壳振动特性分析由广义胡克定律,得到复合圆柱壳与基体及约束层相关的力与弯矩为FkxFkFkxMkxMkMkx=Ak11Ak120Bk11Bk120Ak21Ak220Bk21Bk22000Ak6600Bk66Bk11Bk120Dk11Dk120Bk21Bk220Dk21Dk22000Bk6600Dk66 0(k)x0(k)0(k)xkxkkx=Dk0(k)r 0(k)0(k)x kx k kxT(15)其中:
17、Dk为第k层弹性矩阵。黏弹性层所受的剪力为 SxzSz=Av4400Av55 vxzvz=Sv vxzvz(16)其中:Sv为剪切弹性矩阵。因此,可得各层单元刚度和质量矩阵表达式为Kek=A(Bk)DkBkdAKev=A(BvS)SvBvSdAMek=AkNkNkdA (17)其中:k为各层材料的密度;“”表示矩阵为复矩阵。复合壳单元的刚度矩阵和质量矩阵为Ke =Kes+Ke v+KecMe=Mes+Mev+Mec(18)上面推导了考虑层间变形关系下的复合壳单元的单元矩阵,而对于不考虑变形关系建模时,可假设为 3 层均满足Donnells薄壳理论,并令单元节点位移自由度为q=u ux u v
18、vx v w wx w (19)这里黏弹性层的变形不再由基体层和约束层所表征,所以不再单独考虑式(11)中黏弹性层剪切变形。各层的单元刚度矩阵和质量矩阵可表示为 Kek=A(Bk)DkBkdAMek=AkNNdA(20)2 基体单元和过渡单元推导对于局部贴敷约束阻尼层圆柱壳结构,除包含贴敷约束阻尼层区域外还包含未贴敷区域。由于未贴敷区域的基体单元不包含约束层自由度,假如均利用复合单元进行组集,会出现总刚度矩阵奇异而无法求解的情况,故需要分开考虑这 2 种单元。另外,由于有限元模型中贴敷区的复合单元和未贴敷区的基体单元的自由度不一致,直接将这 2 种单元进行组集同样存在困难。因此,本研究在贴敷区
19、和未贴敷区之间构造了一种过渡单元来连接复合单元和基体单元,形成如图4所示的局部结构单元分布。2.1基体单元构造图 4 中的未贴敷区只含有基体层,首先将复合壳单元中黏弹性层和约束层的参数设为 0,得到矩阵Kes,其维度为60 60。由于单矩阵格式较大,这里只列出单元Kes内的第 1号节点的矩阵形式Ke,1s(21)其次,消去Kes中多余的约束层自由度(式(21)中虚线框区域),形成最终的基体单元刚度矩阵Kes,该矩阵的维数为 3636,每个节点变为 9个自由度。2.2过渡单元构造图 4 中的过渡区也仅包含基体层,但此处的单元要同时与约束阻尼层单元和基体单元相连接。为确保相容性,在构造过渡单元时,
20、需要保证 5 号和6 号节点与约束阻尼单元节点自由度一致(每个节点 15个自由度),3号和 4号节点与基体自由度一致(每个节点 9 个自由度),这样过渡单元各节点就存在不同的自由度。具体构造过程为:令原始过渡单元矩阵Ket=Kes,将Ket中对应 3 号和 4 号的位置消除多余的约束层自由度,保留对应5号和6号节点位置的自由度,得到用于实际计算的过渡单元矩阵Ket,其表达式为图 4局部结构单元分布Fig.4Element distribution of local structure1021振 动、测试与诊断第 43 卷(22)最终的过渡单元刚度矩阵Ket的维数由6060缩减为4848。对于单
21、元质量矩阵进行相同的处理,可得到基体单元质量矩阵Mes和过渡单元质量矩阵Met。通过上述操作可确保局部贴敷黏弹性约束层圆柱壳结构的单元刚度和质量矩阵能进行有效组集。3 有限元方程的建立及求解3.1弹性边界条件的模拟本研究的局部贴敷约束阻尼层圆柱壳的一侧是通过螺栓固定在夹具上,因而圆柱壳受到的是一种弹性约束。参照文献 15,可将底部设置多个周向连续分布的弹簧组来模拟该弹性边界条件,周向连续非均匀分布弹性约束圆柱壳模型如图 5 所示,弹簧组对应的弹簧刚度分别为kx,k,kw,kux和kw。螺栓约束本质上是非均匀的,为模拟圆柱壳底部实际周向的非均匀性,假设连续分布的各弹簧刚度满足余弦函数分布,即ks
22、pr=Kspr(1-)cos(m)+(23)其中:spr表示各边界弹簧的方向;Kspr为弹簧单元刚度值的幅度;m为底部螺栓个数;为螺栓影响区分布因子。Kx,Kux,K,Kw,Kw和可通过实验结果并基于分解的多目标进化算法(multi objective evolutionary algorithm based on decomposition,简 称 MOEA/D)进行反推辨识得到。3.2有限元分析将复合圆柱壳单元、基体单元、过渡单元以及模拟弹性边界的弹簧单元组合在一起,可形成局部贴敷约束阻尼层圆柱壳的自由振动频域方程为(K*+Kspr-2M)X=0(24)其中:K*为复合圆柱壳整体复刚度矩阵
23、;Kspr为边界弹簧刚度矩阵;M为整体质量矩阵;X为特征向量;为自由振动角频率。对式(24)进行特征值求解,可以得到各阶复特征值(r)2=2r(1+ir),进一步可得复合圆柱壳各阶固有频率和模态损耗因子为 fr=Re(r)2 2r=Im(r)2 Re(r)2 (r=1,2,n)(25)4 实 例4.1问题描述这里以局部贴敷约束阻尼层圆柱壳为例验证本研究有限元建模方法。图 6为自由端贴敷约束阻尼层圆柱壳结构,其基体层中性面直径为 450 mm,高度为 210 mm,约束阻尼层贴敷在圆柱壳自由端,其高度为总高度的 1/3,底部由 18 个 M12 规格的螺栓将圆柱壳固定在夹具上。复合圆柱壳各层材料
24、参数如表 1所示。表 1复合圆柱壳各层材料参数Tab.1Material parameters of each layer of composite cylindrical shell材料层基体层黏弹性层约束层密度/(kg m-3)7 850.0789.52 700.0弹性模量/Pa2.1 10114.3 1067.0 1010泊松比0.300.490.30材料损耗因子0.835厚度/mm3.01.00.3图 5周向连续非均匀分布弹性约束圆柱壳模型Fig.5Cylindrical shell model with circumferential continuous nonuniform di
25、stribution of elastic constraints图 6自由端贴敷约束阻尼层圆柱壳结构Fig.6Cylindrical shell structure with constrained damping layer attached on free end1022第 5 期马宏伟,等:考虑层间变形的约束阻尼圆柱壳振动特性分析对复合圆柱壳的振动进行测试,测试结果用于校验所提的有限元模型。贴敷后圆柱壳锤击实验测试系统如图 7所示。采用 PCB模态力锤对复合壳进行激励,Polytec PDV100激光测振仪对图示激光位置进行拾振,LMS SCSDAS 便携式数据采集前端用于获取激励和响
- 配套讲稿:
如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。