椭圆穿孔板超材料的带隙特性计算分析.pdf
《椭圆穿孔板超材料的带隙特性计算分析.pdf》由会员分享,可在线阅读,更多相关《椭圆穿孔板超材料的带隙特性计算分析.pdf(9页珍藏版)》请在咨信网上搜索。
1、Modeling and Simulation 建模与仿真建模与仿真,2023,12(4),3327-3335 Published Online July 2023 in Hans.https:/www.hanspub.org/journal/mos https:/doi.org/10.12677/mos.2023.124305 文章引用文章引用:张祖坚,郭辉,袁涛,孙裴,黄双,马西沛.椭圆穿孔板超材料的带隙特性计算分析J.建模与仿真,2023,12(4):3327-3335.DOI:10.12677/mos.2023.124305 椭圆穿孔板椭圆穿孔板超材料的超材料的带隙特性带隙特性计算分析
2、计算分析 张祖坚,郭张祖坚,郭 辉辉*,袁,袁 涛,孙涛,孙 裴,黄裴,黄 双,马西沛双,马西沛 上海工程技术大学机械与汽车工程学院,上海 收稿日期:2023年3月29日;录用日期:2023年6月27日;发布日期:2023年7月5日 摘摘 要要 为研究二维多孔结构的带隙特性,本文在均质方形板基础上设计一种由三组椭圆孔组成的超材料结构,为研究二维多孔结构的带隙特性,本文在均质方形板基础上设计一种由三组椭圆孔组成的超材料结构,建立椭圆穿孔板有限元模型,利用限元法研究结构在自由状态和建立椭圆穿孔板有限元模型,利用限元法研究结构在自由状态和1 mm指定位移时的带隙特征,分析椭指定位移时的带隙特征,分析
3、椭圆穿孔板的带隙特性、传输特性曲线和带隙边界的振型。结果表明,合理改变结构的边界条件,不仅可圆穿孔板的带隙特性、传输特性曲线和带隙边界的振型。结果表明,合理改变结构的边界条件,不仅可以拓宽带隙宽度,还可将归一化频率以拓宽带隙宽度,还可将归一化频率00.5范围内范围内方向带隙转化为完全带隙,能更有效地衰减和抑制弹方向带隙转化为完全带隙,能更有效地衰减和抑制弹性波传播。性波传播。关键词关键词 边界条件边界条件,超材料超材料,带隙带隙,椭圆穿孔板椭圆穿孔板,有限元法有限元法 Calculation and Analysis of Band Gap Characteristics of Ellipti
4、cal Perforated Plate Metamaterials Zujian Zhang,Hui Guo*,Tao Yuan,Pei Sun,Shuang Huang,Xipei Ma School of Mechanical and Automotive Engineering,Shanghai University of Engineering Science,Shanghai Received:Mar.29th,2023;accepted:Jun.27th,2023;published:Jul.5th,2023 Abstract In order to study the band
5、gap characteristics of two-dimensional porous structures,a metama-terial structure composed of three groups of elliptical holes is designed on the basis of homoge-neous square plates,a finite element model of elliptical perforated plates is established,and the *通讯作者。张祖坚 等 DOI:10.12677/mos.2023.12430
6、5 3328 建模与仿真 bandgap characteristics of the structure in the free state and 1mm specified displacement are stu-died by the finite element method.The bandgap characteristics,transmission loss curves and modes of the bandgap boundary of elliptical perforated panel were analyzed.The results show that r
7、easonable change of the boundary conditions of the structure can not only expand the width of the wide bandgap,but also convert the directional band gap in the normalized frequency range of 00.5 into a complete band gap,which can more effectively attenuate and suppress the propa-gation of elastic wa
8、ves.Keywords Boundary Conditions,Metamaterials,Bandgap,Elliptical Perforated-Panel,FEM Copyright 2023 by author(s)and Hans Publishers Inc.This work is licensed under the Creative Commons Attribution International License(CC BY 4.0).http:/creativecommons.org/licenses/by/4.0/1.引言引言 目前,与多孔结构相结合的机械超材料具有
9、典型的周期穿孔结构,正广泛应用于吸声降噪1 2、超分辨率成像3、水声探测4、抵抗抗冲击5等方面。轻质多孔超材料主要有蜂窝结构、手性和反手性结构,以及穿孔板类结构等。因多孔结具有质地轻、强度高等良好的静力学特性,越来越多研究者将多孔结构应用在减振降噪方面,并取得良好的效果。Gupta等6在对周期性梁和板振动的研究过程中,解释了周期结构中弹性波带隙产生的原因,发现了带隙边界和原胞与模态频率之间的关系。Wu 等7将微穿孔板和卷曲的空间结构组合一起,实现了结构厚度为 50 mm 时,在归一化频率后的带宽达到 0.82。Xiao 等8在薄板基础上附加周期局域振子,构建的板状超材料可大幅提高隔声量。在此基
10、础上,Song 等9进一步引入蜂窝层级板结构,将软橡胶和铝制作的局域振子附加在蜂窝夹层板上,在中高频段获得比同等质量的蜂窝板结构更好的隔声效果。目前对多孔结构的应用研究除了在吸声方面,许多研究人员对多孔结构的减振性能同样作了大量工作。如 Sun 等10 11通过用改进的遗传算法对方形框架结构中材料的分布进行了优化,可有效拓宽和调控低频带隙。Wang 等12 13对方形、“十十”字形和“X”形穿孔板的内部孔型排列方式对结构带隙特性的变化作了深入研究,发现旋转对称的孔型结构更容易打开带隙。Wu 等14在以螺旋孔为基体的原胞的两侧附加钢柱,在 250 Hz 以下产生了较宽的低频带隙。然而,目前对多孔
11、结构的研究除了对经典的方形、圆形等孔型结构研究外,因带隙频率高和带宽较窄,因此有必要对穿孔板结构进行设计研究。本文在第一层级穿孔板结构15 16基础上,引入第二、第三组椭圆孔设计建立复合椭圆穿孔板模型,计算了椭圆穿孔板在自由状态和边界受到指定位移作用下的带隙特性。结果表明,结构边界条件改变可以增大带隙宽度,并且在一定范围内打开更多的完全带隙,在 0.5 范围内打开频率更低的完全带隙。2.结构模型的建立结构模型的建立 本文所述椭圆穿孔板结构是在均质方形板的基础上设计的,结构内部共分为 3 组大小不同的椭圆孔。假设 oxy,oxy是两个坐标系,o 和 o分别为原点,且 o在 oxy 中坐标为(x0
12、,y0),由 x 轴到 x轴的角度为/4,图1中P点坐标(x1,y1),则通过式(1)可得出变换后的P坐标(x1,y1)。不同椭圆孔的几何参数通过公式(1)(3)进行定义。Open AccessOpen Access张祖坚 等 DOI:10.12677/mos.2023.124305 3329 建模与仿真 110110cossinsincosxxyxyxyy=+=+(1)()215iAai=(2)25iiAb=(3)其中,x 和 y 分别表示,坐标变化后椭圆上点的横、纵坐标,表示椭圆孔绕 y 轴旋转的角度;A 表示原胞尺寸,ai表示第 i(i 3)组椭圆孔的长半轴大小,bi表示第 i 组椭圆孔
13、的短半轴大小。如图 1 所示,根据公式(1)将坐标系中心的椭圆变换得到三组不同大小的椭圆孔。第一、二组均由 4个椭圆孔组成,第三组由 8 个椭圆孔组成,如图 2 所示。其中,第 1 组在方形板四个角的位置处,该组的 4 个椭圆长半轴均与各自所在卦限的角平分线垂直,第 2 组和第 3 组椭圆长半轴则均与坐标轴垂直。Figure 1.Schematic diagram of elliptical hole position transformation 图图 1.椭圆孔位置变换示意图 该结构选用的基体材料为尼龙,材料的弹性模量、密度和泊松比分别为 E0=1.8 109(Pa)、0=1010(kg/
14、m)、0=0.38。同时,在结构中间填充半径 r=8mm 的铜质圆柱体,其材料参数分别为 E0=1.8 109(Pa)、0=1010(kg/m)、0=0.38。Figure 2.Elliptical perforated-plate model 图图 2.椭圆穿孔板模型 根据公式(2)与公式(3)建立的椭圆穿孔板模型的相关几何结构参数如图 2 所示。根据椭圆孔几何参数的定义公式,原胞的几何结构参数数值如表 1 所示。张祖坚 等 DOI:10.12677/mos.2023.124305 3330 建模与仿真 Table 1.Geometric parameters value of unit c
15、ell 表表 1.原胞结构几何参数数值 参数 值/mm 说明 A 40 晶格尺寸 h 8 穿孔板厚度 a1 8 第一组椭圆孔长半轴 b1 4 第一组椭圆孔短半轴 a2 6 第二组椭圆孔长半轴 b2 2 第二组椭圆孔短半轴 a3 4 第组椭圆孔短半轴 b3 1 第三组椭圆孔短半轴 R 8 填充体半径 3.带隙机理及布里渊区带隙机理及布里渊区 根据弹性波在固体中传播的波动方程,利用有限元软件 COMSOL 对弹性波在椭圆穿孔板结构中的传播过程进行仿真,通过弹性波理论求解特征方程,得到结构的本征频率。弹性波在 xoy 平面中的传播时,其波动方程如公式(4)所示:()()()()()()()2,xyx
16、yxytttuuu=+rrrrrrr (4)其中,为介质的密度,u(x,y)和 r=(x,y)分别表示位移矢量、位置矢量。此外,因结构周期性可将位移表示为公式(5):()edSidS=k r ruu (5)式中,k 为波矢,us是与各材料参数具有相同周期的函数。当波矢 k 沿着不可约 Brillouin 区(阴影部分)边界扫掠,如图 3 所示。由于弹性波受到内部周期性结构作用形成特殊的色散关系,可绘制出频率 f 和波矢 k 的变化曲线,即能带结构。Figure 3.Irreducible Brillouin zone 图图 3.不可约布里渊区 基于弹性波理论和有限元法计算出各层级的能带结构及相
17、应振型,考虑到精确度和计算效率,本文研究前 20 阶振型。为便于观察带隙频率与材料参数、结构参数之间的关系,将能带结构图的纵坐标数值通过公式(6)进行归一化处理:tfac=(6)张祖坚 等 DOI:10.12677/mos.2023.124305 3331 建模与仿真 其中,为归一化频率,f 为频率,a 为晶格常数,ct为横波波速。4.结果与结果与分析分析 通过位移法研究该椭圆穿孔板结构 1 6 单元的传输损耗模型,如图 4 所示。根据第 2 节中带隙产生的基本理论,计算得到的椭圆穿孔板能带结构如图 5(a)所示,图 5(a)中灰色阴影部分表示完全带隙,黑色阴影部分表示方向带隙。Figure
18、4.Transmission loss model before boundary conditions change 图图 4.边界条件改变前的传输损耗模型 计算得到传输损耗曲线如图 5(b)所示。从图 5(a)中可知,该结构在前 20 条能带打开了 2 个完全带隙,且在 M、X 和 XM 方向上各打开了 2 个方向带隙,其中第一个完全带隙同时也是带宽最大的带隙,范围为0.580.76,分布在第6、7能带之间;第二个完全带隙分布于第12、13能带之间,范围是1.021.17;此外,在第 3 和第 4 能带之间形成了 XM 方向带隙,带宽范围为 0.370.40。(a)(b)Figure 5.
19、BGs characteristics and transmission curve(a)Energy band structure;(b)Transmission loss 图图 5.带隙特性和传输特性曲线(a)能带结构;(b)传输损耗 从图 5(b)可知,计算得到的传输损耗系数在 2 个完全带隙及 6 个方向带隙范围均有不同程度的下降,且均在150 dB 以下,传输损耗系数的大幅衰减表明弹性波在相关带隙范围内可被有效衰减。为进一步阐述带隙打开的机理,将第一个完全带隙和第二个完全带隙上、下边界处的模态作为研究对象,如图 6 所示。从图 5(a)知第一带隙下边界的能带呈平直状态,这是结构产生局
20、域共振所导致的,结合图 6 的 B1可知,弹性波频率等于结构第 6 阶特征频率时,该结构中心几乎保持静止状态,四个角部分等效为局域振 张祖坚 等 DOI:10.12677/mos.2023.124305 3332 建模与仿真 Figure 6.The modes of BGs boundary 图图 6.带隙边界处的振型 子,激发局域振子产生局域共振,且左右两边振子运动方向相反,弹性波传播过程中的能量在该区域被耗散,使得弹性波不能继续向下一个单元继续传播,带隙开始打开。图 6 中的 B2显示,当弹性波频率大于结构第 7 阶特征频率时,该结构四周边界的中心沿着逆时针方向剧烈振动,第二组和三组椭圆
- 配套讲稿:
如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。