近场超声悬浮平台模态振型优化设计及可行性分析

魏 彬， BUCHER I， ATHERTON M， STOLASKI T A

(1.清华大学摩擦学重点实验室，北京 100084；2.西门子研究院(中国)，北京 100102；3.以色列理工学院机械工程学院，海法，以色列 32000；4.布鲁内尔大学工程、设计与物理科学学院，伦敦，英国 UB8 3PH)

近场超声悬浮平台模态振型优化设计及可行性分析

魏 彬1,2， BUCHER I3， ATHERTON M4， STOLASKI T A4

(1.清华大学摩擦学重点实验室，北京 100084；2.西门子研究院(中国)，北京 100102；3.以色列理工学院机械工程学院，海法，以色列 32000；4.布鲁内尔大学工程、设计与物理科学学院，伦敦，英国 UB8 3PH)

为了避免在传送过程中传送带与精密元件之间的摩擦损伤，设计了一种紧凑型非接触超声传送平台，对其可行性进行理论研究和实验验证。通过ANSYS动力学分析建立对称模型，模拟传送平台的模态和传送平台中心点的谐响应振型，结果显示，不同波数条件传送平台呈现出纯弯曲或混合波驱动状态。对所设计的紧凑型传送平台原型进行扫频实验，验证了超声传送过程的可行性以及悬浮频带和振型。理论及实验结果表明，根据不同使用条件，需要对驱动平台进行模态及承载能力的优化设计。研究成果可为超声悬浮平台模态振型和带宽的设计提供参考。

摩擦学；近场超声悬浮；谐振；模态振型；正交波形

为了满足生物工程、医疗器械、半导体等行业对精密元件清洁度、低功耗以及高稳定性生产运输环境的要求，非接触式载运传输平台逐渐进入精密制造领域。近场超声悬浮(near field acoustic levitation，简称NFAL)是近年来应用于微机电技术领域非接触式支撑与传输系统的新技术，相对于其他悬浮形式，具有清洁无污染、高精度、系统结构简单等优点。针对现有精密制造领域元件传输过程对精度保持和环境控制方面的需求，设计开发基于NFAL的无变幅机构超声传送平台新构型，并对其驱动机理和推送能力开展理论和实验研究。超声传送平台可以解决现代制造领域精密元件运输、暂存过程的精度保持及环境控制问题，具有广阔的应用前景。

超声悬浮理论起源于流体动力润滑理论中的挤压效应，最早由GROSS[1]提出，随后, SALBU[2]提出了墙式挤压模型，并假设挤压过程类似于密闭活塞情况，计算了挤压特性和承载力。自由悬浮状态挤压模型最早由BECK等[3]提出, 并通过数值差分方法进行了求解。随后KURODA等[4]提升了模型精度和算法的效率。MINIKES等[5-6]将挤压效应和声悬浮效应相结合，计算了波数和振型的关系。WEI等[7-13]和ALMURSHEDI等[14-15]建立了自由释放过程超声悬浮动力学模型，将差分法和有限元法应用于超声悬浮数值分析中。研究人员[16-20]也对超声悬浮和挤压模型的结果开展了对比研究，超声挤压膜理论得到了进一步发展。

1 挤压理论模型

直角坐标系绝热条件二维Reynolds方程的无量纲形式为[2]

(1)

可压缩瞬态Reynolds方程是二阶非线性偏微分方程，不仅得到解析解非常困难，因为是瞬态方程，所以数值解也极其耗时。为了加快求解速度，以实际情况为基础，对模型进行简化。

为了验证挤压膜的悬浮性能，假设悬浮表面与支撑表面没有横向的相对移动，即u0=0，所以Λ=0。

1)无限宽挤压膜模型

(2)

2)轴对称模型

(3)

3)二维直角坐标模型

对于振源是矩形的激励，当振源的长度、宽度差别不大时，要采用直角坐标系下的二维模型，二维模型计算复杂，计算时间较长。二维模型的控制方程无量纲形式为

(4)

固定悬浮物体的挤压膜模型的控制方程如下：

1) 膜厚的无量纲方程

(5)

2) Reynolds方程的无量纲形式

(6)

边界条件如下：

在每计算完一个周期后需检查定常周期条件是否得到满足，除了频率极低、间隙极大的特殊情况外，在计算完第1个周期后，周期边界条件是不可能得到满足的。将所得到的结果作为第2个周期的初始条件继续计算，当周期性误差小于规定的判据，就可以结束计算而把周期的结果作为问题的解答，最后可以利用无量纲压力积分求出无量纲承载力：

式中f(T)为实际承载力。

利用二阶精度无条件稳定的Crank-Nicholson形式计算此瞬态Reynolds方程(见式(3))，通过Reynolds方程来求解挤压膜的静、动态特性，按照图1所示流程编写程序。

图1 程序计算流程图Fig.1 Program flow chart

2 传送平台原型设计及计算结果分析

所设计的NFAL平台原型采用底部粘贴压电陶瓷片的形式，为了产生正交方向的悬浮波形和驱动波形，可以采用多组压电片交叉分布的形式，采用专业胶水固定，进行去应力处理，超声悬浮传送试验平台及设计原型如图2、图3所示，仿真模型如表1所示，其中材料属性如表2所示。

图2 超声悬浮传送平台设计原型Fig.2 Designed prototype of NFAL platform

图3 超声悬浮传送激振平台设计示意图Fig.3 Schematic diagram of designed NFAL platform

几何模型及材料六面体网格载荷与边界条件

表2 悬浮平台材料属性

悬浮平台电固耦合ANSYS计算结果如表3所示。结果显示出悬浮平台的3种模态振型，第1种为沿x方向的纯弯曲振型，如1阶、3阶、5阶模态；第2种为沿y方向的纯弯曲振型，如2阶、6阶、12阶模态；以及沿x，y方向具有不同波数的混合模态。

计算结果表明，纯弯曲模态振型和混合模态振型为近场超声悬浮平台设计过程中所采用的2种主要振型，其中混合模态呈现出沿x，y正交方向驻波、行波解耦驱动的模式，该模式有助于实现传送平台在法向承载的同时沿某一方向产生推动力，笔者利用横向和纵向的波数来表征混合模态的振型分布，通过S，T标识来区分某一方向的驻波与行波。该分析方法可以针对不同尺度、材料悬浮样本给出最优激振条件，实现超声悬浮正交振型的精确模拟和设计。

表3 超声悬浮平台ANSYS 模态振型

3 实验验证及讨论

采用点阵式激光扫描设备对悬浮平台的固有特性和谐响应特性进行测试，实验测试结果及误差如表4所示[21-22]。对比模拟结果(见表3)可以发现，对称模型可以很好的预计整体悬浮平台的固有特征，所有实验中的谐响应振型在理论模拟中都有所体现，其中峰值频率误差值小于15%，第3阶、第5阶、第9阶较大的误差来源于简化的边界条件。结果表明，并不是所有的谐振点都能产法相的承载能力，纯弯曲波形更适合应用于激振频率较低、振幅较大的场合，这种现象为研究人员对悬浮平台结构和模态的优化提供了一定的指导作用。

表4 点阵扫描设备实验结果

4 结 论

通过对悬浮平台理论和实验的分析，得到以下结论。

1)悬浮平台激振振型对膜厚和压力分布有很大影响，一般不能忽略。

2)对于超声悬浮传送平台，4个压电片环形布置的形式可以产生可观的承载能力，其中混合波形驱动中的行波可以使悬浮体向某一方向推进。

3)几乎所有的谐响应振型都可以通过1/4理论模型进行表达，对称的驱动波形具有更大的承载能力和悬浮稳定性。

4)悬浮平台模型受到驱动器驱动带宽的限制。对于低频悬浮，纯弯曲模态振型比较适合；相反，混合波驱动悬浮平台适合应用于较高的频率。

/References：

[1] GROSS W A. Gas Film Lubrication[M]. New York: Wiley,1962.

[2] SALBU E O J. Compressible squeeze films and squeeze bearings[J]. ASME J of Basic Eng, 1964,86(3):355-366.

[3] BECK J V,HOLIDAY W G, STRODTMAN C L. Experimental and analysis of a flat disk squeeze-film bearing including effects of supported mass motion[J]. ASME J of Lubr, 1969,91(1):138-148.

[4] KURODA S, HIRATA N. The characteristic of motion of a round plate supported on squeeze air film[J]. J of Lubr, 1984, 50(6):2727-2731.

[5] MINIKES A, BUCHER I. Levitation force induced by pressure radiation in gas squeeze films[J]. J Acoustic Soc, 2004,116(1): 217-226.

[6] MINIKES A, BUCHER I. Comparing numerical and analytical solutions for squeeze-film levitation force[J]. Journal of Fluids and Structures,2006, 22(5):713-719.

[7] WEI B, MA X Z, TANG W K. Study on characteristics of ultrasonic levitation with Piezo-ceramics exciting[J]. Piezoelectrics Acoustooptics,2011, 33(1):76-79.

[8] WANG Yanzhong, WEI Bin. Wet multi-disc brake temperature field and optimal oil supply under continuous braking condition[J]. Industrial Lubrication Tribology, 2014, 66(6): 653-661.

[9] WANG Yanzhong, WEI Bin. A linear solution for film characteristics in ultrasonic excitation condition[J]. Industry lubrication and Tribology, 2013, 34(5): 469-473.

[10]魏彬, 马希直. 挤压膜悬浮导轨的悬浮特性研究[J]. 润滑与密封, 2010,35(12): 54-58. WEI Bin, MA Xizhi. Research on the characteristic of the floating guide way with squeeze film[J]. Lubrication Engneering, 2010,35(12): 54-58.

[11]魏彬, 马希直. 考虑激振模态的挤压膜悬浮导轨特性分析[J]. 润滑与密封2010, 35(8): 32-35. WEI Bin, MA Xizhi. Research on the characteristic of the floating guide way including the mode effect[J]. Lubrication Engneering, 2010, 35(8): 32-35.

[12]魏彬，马希直.气体挤压膜悬浮平台的研究方法[J]. 江苏航空学报, 2008(sup1):145-148. WEI Bin, MA Xizhi. Research on the characteristic of the floating guide way with squeeze film[J]. Journal of Jiangsu Aviation Science, 2008(sup1):145-148.

[13]冯美君，魏彬, 佟勇. 自悬浮线性空气轴承的理论分析和轴承设计[J].军民两用技术与产品, 2012,53(1): 53-57. FENG Meijun, WEI Bin, TONG Yong. Analysis and design on self-levitation linear gas bearing[J]. Dual Use Technologies & Products, 2012,53(1): 53-57.

[14]ALMURSHEDI A, ATHERTON M, MARES C, et al. Squeeze-film levitation characteristics of plates excited by piezoelectric actuators[EB/OL]. http://bura.brunel.ac.uk/handle/2438/13663.

[15]ALMURSHEDI A, ATHERTON M, MARES C, et al. Plate actuator vibration modes for levitation[EB/OL].http://bura.brunel.ac.uk/bitstream/2438/11907/1/Fulltext.pdf.

[16]WANG C C. Bifurcation and nonlinear dynamic analysis of united gas-lubricated bearing system[J]. Computers & Mathematicswith Applications, 2012,64(5):127-132.

[17]LIU P, LI J, DING H, et al. Modeling and experimental study on near-field acoustic levitation by flexural mode[J]. IEEE Transaction Ultrasonics Ferroelectrics & Frequency Control,2009, 56(12):2679-2685.

[18]LI J, CAO W W, LIU P K, et al. Influence of gas inertia and edge effect on squeeze film in near field acoustic levitation[J].Appl Phys Lett, 2010, 96(24):96-102.

[19]XIE W J, WEI B. Parametric study of single-axisacoustic levitation[J]. Appl Phys Lett,2001, 79(6):881-884.

[20]YOSHIMOTO S, KOBAYASHI H, MIYATAKE M. Float characteristics of asqueeze-film air bearing for a linear motion guide using ultrasonic vibration[J]. Tribology International, 2007,40(3):503-510.

[21]WANG Yanzhong, WEI Bin. Mixed-modal disk gas squeeze film theoretical and experimental analysis[J]. International Journal of Modern Physics, 2013, 27(25):1-20.

[22]CHANG Xiaoni, WEI Bin, MARK A, et al. NFAL prototype design and feasibility analysis for self-levitated conveying[J]. Tribology Transitions, 2016, 59(5): 957-968.

Modal shapes optimization and feasibility analysis ofNFAL platform

WEI Bin1,2， BUCHER I3， ATHERTON M4， STOLASKI T A4

(1.State Key Laboratory of Tribology, Tsinghua University, Beijing 100084, China; 2.Corporate Technology Siemens, Beijing 100102, China; 3. Faculty of Mechanical Engineering, Israel Institute of Technology, Haifa 32000, Israel; 4. College of Engineering, Design and Physical Sciences, Brunel University, London UB8 3PH, UK)

In order to avoid friction and scratching between the conveyor and the precision components when conveying object, an compact non-contact acoustic levitation prototype is designed, and the feasibility is theoretically and experimentally verified. The symmetry model is established through kinetic analysis with ANSYS. The modal and the coupled field computation at the central point of the transfer platform are simulated. The simulation results show that pure flexural or mixed flexural wave shapes appear with different wave numbers on the platform. Sweep frequency test is conducted on the compact platform prototype. The levitation experimental results confirm the feasibility of the ultrasound transfer process, the levitation frequency range and the mode of vibration. The theoretical and experimental results show that the optimal design of the modal and the carrying capacity of the driving platform is necessary according to different conditions. The research results provide a reference for the design of the mode and bandwidth of the ultrasonic levitation platform.

tribology； near field ultrasound levitation; resonance； modal shapes； orthogonal waves

2016-10-10；

2017-04-20；责任编辑：冯 民

国家自然科学基金(51605250)；北京市自然科学基金(L161001)；博士后科学基金(2017T100071)

魏 彬(1983—)，男，河北唐山人，助理研究员，博士，主要从事超声悬浮、极端环境摩擦学等方面的研究。

E-mail:buaaweibin@126.com

1008-1542(2017)04-0313-07

10.7535/hbkd.2017yx04001

TH133.3

A

魏 彬，BUCHER I，ATHERTON M，等.近场超声悬浮平台模态振型优化设计及可行性分析[J].河北科技大学学报,2017,38(4):313-319. WEI Bin，BUCHER I，ATHERTON M，et al.Modal shapes optimization and feasibility analysis of NFAL platform[J].Journal of Hebei University of Science and Technology,2017,38(4):313-319.