缺陷对薄板结构振动模态影响的数值与实验分析

1、缺陷对薄板结构振动模态影响的数值与实验分析摘要：利用有限元法和电子散斑干涉技术研究了缺陷对矩形薄铝板的振动模态的影响。通过数值模拟和实验 ?y 量同时得到了 1 块完整和 4 块带有缺陷大小与位置均不同的铝板 1 10 阶的振型图；铝板的边界条件为长对边自由，短对边固支。结果表明：缺陷的存在对薄板结构的共振振型影响较大，而对共振频率的影响较小。 共振振型随裂缝的大小、位置和方位的不同而有所变化。研究也表明，有限元计算出的模态形状与实验测量得到的结果一致性较好。 关键词： 振动测量； 模态振型； 含裂缝薄矩形板； 电子散斑干涉； 有限元法中图分类号：O329； TU33+9 文献标志码：A 文章

2、编号：10044523 （ 2017） 04056406DOI：引言板壳、杆及梁等是机械制造等领域常用结构形式12 ，其受外界冲击或振动激励易引起变形，特别是当激励频率接近其共振频率时，变形量达到最大值。其中，矩形薄板因是航空及汽车等工程设计中最为常见的结构类型之一，它的振动响应特性被广泛地研究。研究板存在裂缝缺陷时的振动行为具有重要的实用价值，目前有关含裂缝板的振动研究论文较为有限。裂缝的存在将影响板的静态和动态力学行为，这是因为板的裂缝改变其局部的刚度，使得板的静态挠曲和动态的固有频率发生变化。在有限元方法出现之前，主要用数学物理方法研究带缺陷板的振动特性。 如文献 3 研究有裂纹的矩形板

3、振动时利用格林函数表示板的挠度，进而得到第一类齐次弗雷德霍姆积分方程； Stahl 和 Keer 则利用双级数方程研究矩形板的振动和稳定性问题，最终转化为一个解第二类齐次弗雷德霍姆积分方程 4 。随着振动理论的不断完善及有限元方法的发展，对于由各向同性的材料组成，即便是复杂的结构，有限元方法均能分析和研究其振动特性 57 。振动特性测量的实验方法有传感器法8 、激光多普勒法910 、全息干涉法 11 和电子散斑干涉法等 1213 。传感器法需要将传感器直接粘贴在待测结构表面，容易改变结构局部的质量分布且只能测得离散的数据，从而影响实验结果的精确性；激光多普勒法是非接触测量，响应快、测量范围大，

4、不过它是通过快速逐点扫描进行测量，测试精度和速度依赖于测点数，且测试设备价格昂贵。时间平均的激光全息干涉和激光电子散斑干涉技术适合于薄板、壳的微米级振动变形测量，前者对于光路布置及系统隔振要求较为严苛，因此激光电子散斑干涉技术更有利于振动的测量与分析。全息干涉和电子散斑干涉得到的是稳态振动状态下位移等值线条纹图，直观地反映共振态的振型。文献 29 ，11 主要以单边固支的悬臂板为研究对象，通过理论、有限元模拟及实验方法研究缺陷对平板的振动模态、共振频率的影响，为工程应用提供有益的参考。文献14 则是利用理论模拟和电子散斑干涉技术研究了无缺陷完整平板的单边、对边及周边固支情况下的模态特性。本文采

5、用有限元法和电子散斑干涉技术对4 个含不同类型缺陷的对边固支矩形薄铝板振动特性进行研究；结果表明缺陷的位置对薄板结构的振型具有一定的影响。1 含缺陷薄板振动的数值分析含缺陷薄板结构的控制方程为Md2u （ t） dt2+Cdu（ t ）dt+Ku（ t ） =b（ t ）+f（ u）（ M ， C， K Rn×n； u， b， f Rn）（ 1）式中 u为位移矢量， M 为质量矩阵， C 为阻尼矩阵， K 为刚度矩阵，b（ t ）为与时间有关的外力，f（ u）为裂缝表面间断接触引起的非线性接触力。利用有限元软件ANSYS对铝板结构进行数值模拟，为了便于建模，假定含缺陷铝板结构为各向同

6、性的线弹性材料，结构的形变无穷小，所以忽略了式（1）中非线性接触力f（u）的影响。本文矩形铝板结构有限元模型的物理及几何参数为： 弹性模量 E = 70 GPa，泊松比 = 0.33，密度 = 2700 kg/m3 ；试样大小为180 mm × 80 mm ×0.9mm ，边界约束条件为短对边固支长对边自由。有限元计算时选用 SHELL 181单元，并将壳单元划分为多个1 mm × 1 mm× 0.9 mm 的小单元进行计算。为了保证实验时测试试样尺寸与有限元模型一样，考虑到约束边界实际加工试样尺寸为260 mm × 80 mm ×

7、0.9 mm，具体试样示意图如图1 所示。图 1 线切割加工的（ a）完整及（ b）（ e）含预制边缘、内裂缝的薄铝板示意图Fig.1Schematic of thin Aluminum plates（ a） Intact and（b）（ e） with slits cut by EDM 第 4 期姜益军，等：缺陷对薄板结构振动模态影响的数值与实验分析振 动 工 程学 报第 30 卷 2 激光电子散斑干涉振动测试技术图 2 是本文测量薄铝板离面振动特性的系统示意图。参考图 2 的左半部分，根据干涉理论，t 时刻两束光干涉形成的散斑强度为 1415I（x，y， t ）=Io+Ir+2IoIrco

8、s（ x， y）+ （ x， y，t ） （ 2）式中 Io， Ir 分别为经试样漫反射的物光束与由漫反射面反射的参考光束的光强， （x， y）为试样静止时参考光与物光间的随机相位差； （x， y， t ）则与试样振动引起的物光光程变化有关。对于纯离面振动的情形，即面内变形或位移为 0，可知 （ x，y，t ）=4 Acos（ t ）（ 3）式中 Acos （ t ）为 t 时刻物面上点（ x， y）处的离面位移，其中 A表示离面振动的振幅， 振动圆频率， 为激光波长。事实上， CCD相机采集一幅图像需要一定的时间，因此计算机显示的图像实际上是CCD在成像周期内对光强信息的积分。如果 CCD在

9、 t 与 t+时间内曝光进行光电转换， 则其输出数字图像的灰度为g（x， y）=k t+ tIo+Ir+2IoIrcos （x， y） +（ x， y， t ）dt （ 4）式中 k 为 CCD光电转化系数， 表示曝光时间。假设相机曝光时间 是试样振动周期的整数倍（即 =2n / ， n为正整数），上式可进一步写成g（x， y）=2nk Io+Ir+2IoIrJ0（ m） ?cos（ x，y） （ 5）式中 J0（ m）为如图3 所示的第一类零阶Bessel 函数，m=4 A/ 。假设振动过程中激振器激励力的振幅微小波动，导致试样振动的振幅从 A 改变成 A+ A，则此时 CCD在另一个曝光周

10、期 内获得图像的灰度为f（ x， y） =k 0Io+Ir+2IoIrcos （ x， y） + 4（ A+ A） cos（ t） dt （ 6）令 m=4 / ，式（ 6）可改写成f（ x， y） =k 0Io+Ir+2IoIrcos （ x， y） +（ m+ m） cos（ t） dt （ 7）将 cos（m+m）cos（t ）在 m 处泰勒展开并忽略高级小量，式（ 7）简写成f（ x， y） =2nk Io+Ir+2IoIr1-（ m ）22J0（m） cos（ x， y）（ 8）式（ 5）与式（ 8）相减并取绝对值得到振幅涨落电子散斑干涉条纹灰度表达式为16G（x， y） =2 nk

11、 （m）2IoIr|J0 （ m） cos（ x，y） |（ 9）图 2 离面振动测试系统示意图Fig.2The setup for outofplane vibration testing图 3 零阶 Bessel 函数Fig.3First kind zeroorder Bessel function3 有限元及实验结果图 2 所示的测量系统中，激光器产生波长为532 nm 绿色的相干光，功率 050 mW 连续可调； CCD为德国 IDS公司的 1280× 1024 像素可编程控制的相机；信号发生器可产生 0 20000 Hz 正弦波，信号发生器、功放器及激振器均为江苏联能公司制

12、造。试样由贴在其后表面薄压电陶瓷激励，陶瓷直径10 mm ，厚 0.4 mm 。为了对比， 图 4 给出有限元法模拟计算及实验测量获得的图 1（ a）所示的无缺陷完整、长对边固支短对边自由的薄铝板试样 1 的离面 110 阶振型云图和散斑干涉条纹图。由图 3 及式（ 9）知图 4（A）（ J）散斑干涉条纹图中最亮的条纹对应于振动节线（驻点线）。由于铝板固支部分无法达到完全刚性夹持，导致对应于图4（ B）与图 4（F）的一阶纯扭及二阶纯扭的节线位置与有限元计算结果有一定偏差。图 5 是模拟与测量得到的图1（ b）所示试样2 的 1 10阶振型图。计算模型和测试试样中预制边缘裂缝距下固定边90 m

13、m 处，缝宽 1 mm ，长 32 mm 。由图 5（a）（ c），（e）与（ f）可以看出，对于居中对称的边缘裂缝，其对纯弯及纯扭振动的振型影响很小；但对弯扭组合振动的振型影响较大。模拟与测量得到的图1（c）所示试样 1 10 阶振型图如图 6 所示。计算模型和测试试样中预制边缘裂缝距下固定边115 mm 处，缝宽 1 mm ，长 30 mm 。由图 6（ a）（ b）看出，对于非居中的边缘裂缝，除了对一阶纯弯及一阶纯扭振动的振型影响较小外；对其他高阶的弯、扭及弯扭组合振动的振型影响都较大。图 7 是对应于图 1（ d）所示试样的 1 10 阶振型图。 该试样中距下固定边 60 mm 处有一

14、与其平行的预制内部裂缝，裂缝长 50 mm ，宽 1 mm。图 7（a）（ d）、图 7（A）（ D）以及图 7（f）（ g）表明，该内部裂缝对低阶弯、扭及弯扭振动的振型影响不大。含内部斜裂缝的薄铝板110 阶振动振型图如图8 所示。该试样的裂缝以对称轴上?c 为中心对称，长 50 mm ，宽 1 mm ，与对称轴成60°。模拟与实验结果表明，该斜裂缝仅对前4阶振型影响较小。图8（ H）和图 8（ I）的 8，9 阶振型，与图 8（ h）和图 8（ i）有限元模拟得到振型位次互换，即出现了模态的跃迁。表 1 为实验中 5 个试样 1 10 阶共振频率的测量和有限元模拟的结果。4 讨论

15、与结论由于实验中试样夹持边界条件不是理想的固支约束，导致测得纯弯模态，特别是一阶纯弯的共振频率与有限元模拟值相差较大（参见表1）。模拟和实验结果表明，缺陷对薄板各阶的模态和共振频率都会产生影响，缺陷位置的不同对振型影响也不同，缺陷对共振频率的影响相对较小。对比图 5（ f），（F）及图 6（ g），（ G）可知，缺陷的存在除了上述影响外，缺陷位置和大小的不同有可能导致薄板振动模态的跃迁现象产生，这在文献6 有限元图4 完整薄铝板1 10 阶离面振动振型：（a）（ j）有限元计算； （ A）（ J）实验测试结果Fig.4Mode shapes of first 10 resonant mode

16、of CFCFAluminunm plate from ：（ a）（ j） FEM simulation and （ A）（ J） experiment measurement ， respectively 图 5（a）（j）和（ A）（ J）分别对应于含边缘居中裂缝薄板1 10 阶离面振动振型的模拟和实验结果Fig.5Mode shapes of first 10 resonant mode of Aluminunmplate with central edge slit results from（ a）（ j） FEM and（ A）（ J） experiment methods ， re

17、spectively 图 6（ a）（j）和（ A）（ J）分别为含边缘非居中裂缝薄板1 10 阶离面振动振型的模拟和实验结果Fig.6Mode shapes offirst 10 resonant mode of Aluminunm plate with noncentral edge slit results from （ a）（ j） FEM and （ A）（ J） experiment methods ， respectively?D7 含内部裂缝薄板 1 10 阶离面振动振型图： （ a）（ j）为计算值；（ A）（ J）对应于实验测量值Fig.7Mode shapes of fi

18、rst 10 resonant mode of Aluminunmplate with inner horizontal slit results from（a）（ j） FEM and（ A）（ J） experiment methods ， respectively 图 8 内部斜裂缝薄板 1 10 阶离面振动振型图： （ a）（ j）为计算值；（ A）（ J）对应于实验测量值Fig.8Mode shapes of first 10 resonant mode of Aluminunmplate with inner skewed slit results from（ a）（ j） FEM

19、 and（ A）（ J） experiment methods ， respectively表 1 测试试样的1 10 阶共振频率（ Hz）Tab.1First 10 resonant frequencies （ Hz） of the testedspecimens from FEM and experimental methods ，respectively阶数试样 1 试样 2 试样 3 试样 4 试样 5 计算值实验值计算值实验值计算值实验值计算值实验值计算值实验值.8 1309 分析含缺陷悬臂薄板振动时也有类似现象。从弹性力学及振动理论观点来说，薄板的局部损伤导致结构局部刚度产生变化，

20、引起非线性响应导致模态跃迁。本文研究表明，在一定假设前提下，有限元模拟计算的结果必然与实验的结果有一定的差异。因此，在计算机模拟的基础上开展实验研究是很有必要的。参考文献：1Dimarogonas A D. Vibration of cracked structures ： astate of the art reviewJ. Engineering Fracture Mechanics ，1996， 55（ 5）： 831 857.2Qiu Z C， Wu H X， Ye C D. Acceleration sensors based modal identification and act

21、ive vibration control of flexiblesmart cantilever plateJ. Aerospace Science and Technology ，2009， 13（ 6）： 277 290.3Lynn P P， Kumbasar N. Free vibration of thin rectangular plates having narrow cracks with simply supported edgesC. Proceedings of 10th Midwestern Mechanics Conference， 1967： 911 928.4St

22、ahl B， Keer L M. Vibration and stability of cracked rectangular platesJ. International Journal of Solids and Structures ， 1972， 8， 69 91.5Wu G Y， Shih Y S. Dynamic instability of rectangular plate with an edge crackJ. Computers and Structures ， 2005，84： 110.6Saito A， Castanier M P， Pierre C. Estimat

23、ion and veering analysis of nonlinear resonant frequencies of cracked platesJ. Journal of Sound and Vibration ， 2009， 326： 725739. 7Huang C S， Leissa A W， Li R S. Accurate vibration analysis of thick ， cracked rectangular platesJ. Journal of Sound and Vibration ， 2011， 330：2079 2093.8Yang Z B， Chen

24、X F， Yu J， et al. A damageidentification approach for plate structures based on frequencymeasurementsJ. Nondestructive Testing and Evaluation ，2013， 28（ 4）： 321 341.9Sundaresan M J， Ghoshal A， Li J， et al. Experimental damage detection on a wing panel using vibration deflection shapesJ. Structural H

25、ealth Monitoring ， 2003，2（ 3）：243 256.10Qiao P Z， Lu K， Lestari W， et al. Curvature mode shapebased damage detection in composite laminated platesJ. Composite Structures ， 2007，80： 409428.11Maruyama K ， Ichinomiya O. Experimental study of free vibration of clamped rectangular plates with straight na

26、rrow slitsJ. JSME Int. J. Ser. III， 1989，32： 187193.12Ma C C， Huang C H. Experimental and numerical analysis of vibrating cracked plates at resonant frequenciesJ. Experimental Mechanics ， 2001， 41（ 1）： 8 18.13 季业 .电子散斑干涉测量技术在薄板结构模态分析中的研究 D. 天津：天津大学， 2012.Ji Ye. The study on model analysis of thin pl

27、ate based onElectronic Speckle Pattern InterferometryD. Tianjin ： Tianjin University ， 2012.14Butters J N， Leendertz J A. Speckle pattern andholographic techniques in engineering metrologyJ. Optics andLaser Technology， 1971，3（ 1）：26 30. 15Creath K. Slettemoen G A. Vibration observationtechniques for

28、 digital specklepattern interferometryJ. Journal of the Optical Society of America A ， 1985， 2（10）： 1629 1636.16Wang W C， Hwang C H， Lin S Y. Vibrationmeasurement by the timeaveraged electronic speckle pattern interferometry methodsJ. Applied Optics ， 1996， 35（ 22）： 4502 4509.Investigation on vibration characteristics of thin plates withdefect by numerical and experimental methodsJIANG Yijun ，JIANG Shan， TAO Nan， YANG Fujun， HE Xiaoyuan（ Department of Engineering Mechanics ， SoutheastUniversity ， Nanjing 210096 ， China） Abstract： Vibration characteristics of thin rectan