基于Lamb波高阶模态厚度共振的薄片厚度表征_孙铭.pdf
《基于Lamb波高阶模态厚度共振的薄片厚度表征_孙铭.pdf》由会员分享,可在线阅读,更多相关《基于Lamb波高阶模态厚度共振的薄片厚度表征_孙铭.pdf(8页珍藏版)》请在咨信网上搜索。
1、第 41 卷第 6 期2022 年 12 月Vol.41,No.6Dec.,2022声学技术Technical Acoustics基于Lamb波高阶模态厚度共振的薄片厚度表征孙 铭1,2,袁 玲1,孙凯华2,沈中华1,钟 华2(1.南京理工大学理学院,江苏南京 210094;2.中国工程物理研究院机械制造工艺研究所,四川绵阳 621900)摘要:爆破片是气路控制系统阀门中的关键元件,其厚度是影响压力开关阈值的关键参数,对其厚度的精确表征具有重要意义。为此,研究了薄片的Lamb波高阶模态的厚度共振响应与薄片厚度之间的关系。对200 m和50 m两种不同厚度的不锈钢薄片的厚度共振模式进行了数值模拟
2、,并通过厚度共振频率和厚度之间关系的理论公式计算其厚度。建立了一个使用超声纵波直探头激发、探测不锈钢薄片共振模式的实验系统,得到两种不锈钢薄片的厚度共振频率,通过理论公式计算厚度并与点探头测厚仪的测量结果对比。利用模拟及实验得到的共振峰频率所计算的厚度与实际的厚度测量结果均吻合较好。利用Lamb波在截止频率处可以产生多模态的厚度共振,这些模态均可用于薄片的高精度厚度表征。关键词:Lamb波;厚度共振;数值模拟;超声探头;不锈钢薄片中图分类号:O426 文献标志码:A 文章编号:1000-3630(2022)-06-0838-08Characterization of thin sheet th
3、ickness by high-order mode thickness resonance of Lamb waveSUN Ming1,2,YUAN Ling1,SUN Kaihua2,SHEN Zhonghua1,ZHONG Hua2(1.College of science,Nanjing University of Science and Technology,Nanjing 210094,Jiangsu,China;2.Institute of Mechanical Manufacturing Technology,China Academy of Engineering Physi
4、cs,Mianyang 621900,Sichuan,China)Abstract:The rupture disc is a key element in the valve of the gas circuit control system,and its thickness is a key parameter that affects the pressure switch threshold.Therefore,for the accurate characterization of its thickness,the relationship between the thickne
5、ss resonance response of the Lamb wave in the thin sheet and the thickness of the thin sheet is studied.Numerical simulations of thickness resonance modes in two stainless steel sheets with different thicknesses of 200 m and 50 m are conducted,and their thickness values are calculated by the theoret
6、ical formulas that expresses the relationship between the thickness and the corresponding thickness resonance frequency.An experimental system is established to excite and detect the resonance mode of the stainless steel sheet with a straight ultrasonic longitudinal wave probe,and the thickness reso
7、nance frequencies of two stainless steel sheets are obtained.The thicknesses are also calculated by the theoretical formula and compared with the measurement values of the point probe thickness gauge.The results show that the thicknesses of the two stainless steel sheets obtained by simulation and e
8、xperiment are in good agreement with the actual measured values.So,Lamb waves can be used to generate multi-mode thickness resonance at the cut-off frequency and these modes can be used for high-precision thickness characterization of thin sheet.Key words:Lamb wave;thickness resonance;numerical simu
9、lation;ultrasonic probe;stainless steel slice0引 言在石油、化工、药品及食品等行业中,具有高温高压性质的产品通常存放于密闭的容器中,而容器一旦超压将会发生爆炸,后果不堪设想。爆破片以泄压面积大、无泄漏等优点,已成为应用最广泛的超压泄放装置。爆破片的爆破压强与厚度有关,故对金属爆破片厚度的精准测量对爆破片的设计、制造及使用具有重要指导意义1。超声测厚技术由于其轻便、测量速度快、对人体无害、灵敏度高等特点,广泛用于物体厚度的测量2。超声测厚有许多种方法,常用方法为体波法和Lamb波法等。体波法可以分为脉冲回波法和共振法等。脉冲回波法是通过测量两相邻回波
10、的传递时间,结合已知声速来测量厚度的方法。Kruger等3利用脉冲回引用格式:孙铭,袁玲,孙凯华,等.基于Lamb波高阶模态厚度共振的薄片厚度表征J.声学技术,2022,41(6):838-845.SUN Ming,YUAN Ling,SUN Kaihua,et al.Characterization of thin sheet thickness by high-order mode thickness resonance of Lamb waveJ.Technical Acous-tics,2022,41(6):838-845.DOI:10.16300/ki.1000-3630.2022.
11、06.008收稿日期:2021-04-26;修回日期:2021-06-08基金项目:国家自然科学基金资助项目(61975080、61801451)。作者简介:孙铭(1997),男,吉林白城人,硕士研究生,研究方向为激光与物质相互作用。通信作者:孙凯华,E-mail:第 6 期孙铭等:基于Lamb波高阶模态厚度共振的薄片厚度表征波法在1 250测量厚度达100 mm的钢样品时,证明了此方法的可行性。Takimoto等4采用16个频率为5 MHz的超声换能器对管道壁进行检测,得到了厚度随时间与温度变化的图像。共振法是材料厚度为半波长/2的整数倍时,反射波与入射波会相互干涉形成驻波,测得两个相邻的共
12、振频率即可计算材料的厚度。Dixon等5利用横波的体波共振模式,测量了220 m铝板上的厚度约为11 m聚合物涂层的厚度。Shekhawat等6开发了一种无损扫描近场厚度共振声显微镜,激发了样品厚度共振模式,观察其相位的变化以提供表面的纳米级特征。Terence等7介绍了一种新的厚度测量技术,使用波长相对较长的空气耦合的超声波作用于无孔膜,估计一些金属和非金属样板材料的厚度,通过实验验证了该技术。厚度的测量精度通常为其标称值的6%以内,该技术的精度范围约为10100 m。Lamb波法可分为Lamb波频散法和Lamb波的零群速度共振法等。Lamb波频散法是利用Lamb波频散曲线与厚度的相关性,通
13、过对超声信号进行时频处理,提取各模态频散曲线从而得到厚度。许凯亮等8准确地提取了 3、4、5 mm 厚铝板的宽带Lamb波频散曲线,利用不同时频分析方法对铝板的厚度进行估计,厚度的估计值与板厚的实际测量值的误差分别为0,6.25%和5%。Ma等9利用信号在双层介质中的Lamb波传播特性来测量铝板表面附着污垢的厚度,对于 2 mm 厚铝板可计算 0.53 mm污垢层厚度。零群速度法的原理为Lamb波在群速度为零处的能量不传播,在频域会产生一个尖锐明显的共振峰信号,利用该信号对厚度的敏感性可以求得厚度。Gruensteidl等10提出了一种同时确定各向同性板厚度及纵向和剪切弹性波速度的方法,对厚度
14、测量的误差为1.7%,纵波和横波速度的测量误差为1.7%。Prada等11通过实验研究了零群速度模式的共振效应、后退波的传播等特性,同时讨论了各向同性板中零群速度共振模式存在所需的条件。由于爆破片的厚度较薄,利用脉冲回波法进行厚度测量时,需要高频探头进行激发同时需要带宽较宽的探头进行探测。Lamb波频散法需要对超声信号进行时频分析处理12,复杂的频散曲线易引入不必要的误差,而基于Lamb波的零群速度法在实验中不易观察到明显的零群速度共振现象。因此,本文以不锈钢薄片为研究对象,对体波法和Lamb波法进行了结合,基于Lamb波高阶模态共振理论,利用Lamb波在截止频率处产生的多模态厚度共振,采用有
15、限无探头激发实验的方法,对金属薄片的厚度进行表征。1Lamb波厚度共振的理论研究Lamb波是在自由板中产生的平面应变波,随着波的入射角和频率的改变,会产生不同的模态结构13。Lamb波具有衰减小、传播距离远、信号强、多模态等特点,可广泛用于无损检测14和评估15-19。在厚度d=2h的板界面处,对Lamb波的传播方程和边界条件进行求解,可得到对称模式和反对称两种模式,分别如式(1)、(2)所示:tan(qh)tan(ph)+4k2pq(q2-k2)2=0(1)tan(qh)tan(ph)+(q2-k2)24k2pq=0(2)其中:p=2c2L-k2,q=2c2T-k2,k=cp,=2f,k为波
16、数,cL为纵波速度,cT为横波速度,cp为相速度,f为Lamb波频率,为角频率。方程(1)、(2)的所有解都可以通过数值计算求得,并以色散曲线的形式表示,这些色散曲线是材料特性的独特表征,是利用板厚d来归一化波数与频率的色散关系。例如,不锈钢316L金属薄片(cL=5.678 mm s-1,cT=3.067 mm s-1)的色散曲线如图1所示,其力学参数如表1所示。图1中,对称模式(实线)为Sn,反对称模式(虚线)为An,fd表示频厚积,kd为波数与厚度的乘积。当波数k=0时,厚度共振响应就会出现,厚度共振是空穴共振,类似于光学中的法布里-珀罗共图1 不锈钢316L材料的色散曲线Fig.1 D
17、ispersion curves of 316L stainless steel表1不锈钢316L的力学参数Table 1Mechanical parameters of 316L stainless steel杨氏模量/GPa194泊松比0.294密度/(gcm-3)7.97/GPa106.98/GPa74.968392022 年声学技术振,这些厚度共振模式与纵向或横向体波相关联20。由于前两个模式A0和S0不存在截止频率(当k=0时,=0),故也不存在对应厚度共振模式,除A0和S0外的所有高阶模式都具厚度共振模式。根据Mindlin21提出的分类,可利用对称性和阶数的奇偶性分解厚度共振模式
18、,对于对称模态,偶数阶模态S2n和奇数阶模态S2m+1的共振频率分别如式(3)、(4)所示:fc=ncTd(3)fc=(m+12)cLd(4)对于反对称模态偶数阶模态A2n和奇数阶模态A2m+1的共振频率分别如式(5)、(6)所示:fc=ncLd(5)fc=(m+12)cTd(6)式中:n为正整数,m为自然数。同时,厚度共振又可分为厚度剪切共振(模式S2n或A2m+1)或厚度拉伸共振(模式S2m+1或A2n),则可将上面两式进行合并,对于厚度拉伸共振:fc=n2dcL(7)对于厚度剪切共振:fc=m2dcT(8)其中:m和n均为正整数。利用不同厚度共振的共振峰频率值,可通过下式求得厚度值:de
19、=nc2fn(9)正整数n表示不同的共振模式,fn为共振模式对应的共振峰频率测量值,de为厚度的计算值,c为超声波的波速,当共振模式为厚度拉伸共振时,c为纵波速度;当共振模式为厚度剪切共振时,c为横波速度。可以看到,式(9)与体波的厚度共振的公式是一致的,这也说明了Lamb波高阶模态的厚度共振实质上是多模态的体波共振模式。由公式(3)(6)可知,在板厚固定的情况下,截止频率仅取决于体波的速度。根据Mindlin18提出的分类,对称模式与反对称模式的排列顺序仅取决于截止频率的大小。通过图1的色散曲线,提取了不锈钢316L材料金属薄片的各厚度共振模式对应频厚积fd,结果如表2所示。2Lamb波厚度
20、共振的数值模拟为更直观看到Lamb波的厚度共振模式,利用有限元方法对Lamb波厚度共振进行了仿真模拟。仿真的对象为不锈钢316L材料的薄片,其力学参数已在表1中给出。模拟选用固体力学模块,为方便计算,将位移加载情况等效为二维平面应变模型,在几何模块中构建矩形不锈钢薄片如图2所示,模型长度为16 mm,厚度分别为200 m和50 m。为研究厚度剪切共振与厚度拉伸共振模式,在图2材料的下方添加1 mm宽度的指定位移:u=Acos2f(t-t0)e-(t-t0tw)2e-(xr0)2(10)式中:振幅A为100 nm,f为中心频率,对200 m和50 m的不同厚度,为满足频域宽度需求,f分别取20
21、MHz和50 MHz,t0为时延,tw为脉冲宽度,r0为空间尺寸半径。为计算出超声波,必须选取合适的时间步长t,在保证求解精度的情况下,时间步长一般通过式(11)来计算,同时为了在计算中保证足够的空间分辨率,利用式(12)来计算单元格的长度:t=120fmax(11)Le=min10(12)式中:fmax是所期望的最高频率,min为最小波长,Le为单元的长度。经过计算,对于200 m厚的薄片,取时间步长为1 ns,最大单元的长度为10 m;对于50 m厚的薄片,取时间步长为0.1 ns,最大单元的长度为10 m。为减小纵波或横波在模型边界处发生相互转换表2不锈钢316L薄片对称和非对称厚度共振
22、模式对应的频厚积Table 2The corresponding frequency-thickness products of the symmetric and asymmetric thickness resonance modes of 316 L stainless steel slice模式对称模式模态S0S1S2S4S3S6S8S5fd/(MHz mm)2.843.076.138.519.2012.2614.19模式反对称模式模态A0A1A3A2A5A7A4A9fd/(MHz mm)1.534.605.687.6610.7311.3613.80图2 仿真模型示意图Fig.2 Si
23、mulation model diagram840第 6 期孙铭等:基于Lamb波高阶模态厚度共振的薄片厚度表征对模拟结果造成的影响,对于厚度拉伸共振的模拟,只加载垂直方向的指定纵波位移并提取对心位置的离面位移数据;对于厚度剪切共振的模拟,只加载水平方向的指定横波位移并提取对心位置的面内位移数据。对200 m厚的薄片两种加载方式的对心位置超声信号如图3(a)、3(b)所示。利用傅里叶变换得到两时域信号对应的频谱图,归一化后的结果分别如图4(a)、图4(b)所示,在两图中可清晰看到多个厚度共振峰,各厚度共振峰对应的共振模式已在图中标出。将两图中的各共振峰频率值代入式(9),计算得到的各模式对应的
24、厚度值及与给定厚度值的相对误差已在表3、4中给出。除A6共振模式外,各共振模式计算得到的厚度的相对误差均小于1%,利用模拟得到的数据计算结果与模拟标准值十分吻合。同理,对于50 m厚薄片,两种加载方式得到的超声信号对应的频谱图,归一化后的结果如图5(a)、5(b)所示。采用相同的计算方法得到的厚度值及相对误差表5、6所示。相比于200 m厚薄片的计算结果,50 m薄片的计算误差均有所增大,但利用前三种共振模态计算得到的厚度相对误差仍小于1.5%。对Lamb波厚度共振的模拟结果总结如下:对于200 m厚不锈钢薄片的厚度拉伸共振和厚度剪切共振的模拟,除 A6厚度共振模式相对误差为1.57%外,误差
25、值均在0.7%以下,分析A6模式的相对误差增大的原因为加载位移中心频率为 20 MHz,A6共振模式对的应共振频率在 80 MHz 以上,加载的频域振幅在高频处会减弱而导致偏差,图4模拟厚度为200 m薄片的厚度拉伸和剪切共振频谱图Fig.4Spectra maps of the thickness stretching and shear reso-nances of the simulated 200 m thick slice表3200 m薄片各厚度拉伸共振模式对应厚度及相对误差Table 3Thickness and its relative error of 200 m thick
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 Lamb 波高阶模态 厚度 共振 薄片 表征 孙铭
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。