[go: up one dir, main page]

CN116644644B - Restoring force prediction method based on eccentric and inclined states of suspended object - Google Patents

Restoring force prediction method based on eccentric and inclined states of suspended object Download PDF

Info

Publication number
CN116644644B
CN116644644B CN202310872125.2A CN202310872125A CN116644644B CN 116644644 B CN116644644 B CN 116644644B CN 202310872125 A CN202310872125 A CN 202310872125A CN 116644644 B CN116644644 B CN 116644644B
Authority
CN
China
Prior art keywords
suspended object
restoring force
air film
suspended
eccentricity
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202310872125.2A
Other languages
Chinese (zh)
Other versions
CN116644644A (en
Inventor
刘远远
李鑫
郑妍
孙先涛
陈文杰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Anhui University
Original Assignee
Anhui University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Anhui University filed Critical Anhui University
Priority to CN202310872125.2A priority Critical patent/CN116644644B/en
Publication of CN116644644A publication Critical patent/CN116644644A/en
Application granted granted Critical
Publication of CN116644644B publication Critical patent/CN116644644B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01LMEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
    • G01L11/00Measuring steady or quasi-steady pressure of a fluid or a fluent solid material by means not provided for in group G01L7/00 or G01L9/00
    • G01L11/004Measuring steady or quasi-steady pressure of a fluid or a fluent solid material by means not provided for in group G01L7/00 or G01L9/00 by the use of counterbalancing forces
    • G01L11/006Measuring steady or quasi-steady pressure of a fluid or a fluent solid material by means not provided for in group G01L7/00 or G01L9/00 by the use of counterbalancing forces hydraulic or pneumatic counterbalancing forces
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/28Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/08Fluids
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Algebra (AREA)
  • Geometry (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Data Mining & Analysis (AREA)
  • Fluid Mechanics (AREA)
  • Computing Systems (AREA)
  • Operations Research (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Extrusion Moulding Of Plastics Or The Like (AREA)

Abstract

The invention belongs to the technical field of non-contact suspension, and particularly relates to a restoring force prediction method based on eccentric and inclined states of a suspended object, which comprises the steps of obtaining resonant frequency and amplitude distribution of an ultrasonic transducer; establishing a restoring force prediction theoretical model taking eccentricity and inclination of a suspended object into consideration; deducing a fluid motion control equation in the extrusion film based on a gas lubrication theory; obtaining the inclination angle and the restoring force of the suspended object by using a finite difference method and a spline interpolation method; and establishing a mathematical relationship between the restoring force and the eccentric distance, and verifying with an experimental result. The invention discloses a mechanism for generating restoring force under the near-field acoustic levitation eccentricity from the angle of extrusion fluid movement, and provides a restoring force prediction method which simultaneously considers the levitation object eccentricity and inclination, so that the magnitude of the restoring force when the levitation object is eccentric can be accurately predicted, and a theoretical basis is provided for the application of the near-field acoustic levitation in the field of non-contact precise levitation.

Description

Restoring force prediction method based on eccentric and inclined states of suspended object
Technical Field
The invention belongs to the technical field of non-contact suspension, and particularly relates to a restoring force prediction method based on an eccentric and inclined state of a suspended object.
Background
The principle of near-field ultrasonic levitation is to utilize vibration at ultrasonic frequencies to generate a nonlinear extrusion flow field. Under the excitation of sine alternating voltage, piezoelectric ceramics contained in the ultrasonic transducer converts electric energy into mechanical energy by utilizing inverse piezoelectric effect, so that the ultrasonic transducer periodically vibrates at a certain frequency and amplitude. The ultrasonic vibration of the radiating disk in the ultrasonic transducer can generate an extrusion flow field between the radiating disk and a suspended object. In the extrusion flow field, the average air pressure value in one period is larger than the ambient air pressure value, so that levitation force is generated to realize non-contact levitation.
Because the near-field ultrasonic levitation has the advantages of compact structure, easy control, no physical contact and the like, the positioning and transmission system based on the near-field ultrasonic levitation can be applied to levitation and wafer transfer in semiconductor manufacturing. In addition, the extrusion air film bearing based on near field ultrasonic suspension has the advantages of no friction, low maintenance cost, high speed and the like.
For near-field ultrasonic levitation, the flow field distribution generated by vibration is symmetrical in an ideal situation, so that the levitated object can be stabilized in an equilibrium position. This symmetry breaks when the levitated object is eccentric with respect to the radiating disk. This asymmetric flow field distribution produces a restoring force to urge the suspended objects back to an equilibrium position due to the presence of gas viscous forces.
The mechanism and prediction of the generation of the restoring force is critical for applications of near-field acoustic levitation in the field of precision levitation. In the aspect of theoretical research, some scholars propose a spring vibrator model based on the near boundary acoustic flow theory to estimate the stability and restoring force of a suspended object. The surface restoring force has a large relationship with the amplitude distribution of the radiating disc, the mass of the suspended object and the offset position as a result of the study. However, due to its linear assumption, the spring vibrator model can only handle small amplitude oscillations. Other models for studying restoring forces are based on time-averaged local potential energy methods or gas lubrication theory. It is not to be neglected that the deflection and tilting of the suspended mass are symbiotic, so that the model is not accurate enough in either case of analysis alone, especially for larger deflection displacements, the influence of the tilting on the gas film pressure distribution and the magnitude of the restoring force is greater. Up to now, no model is established that takes both cases into consideration at the same time. Therefore, establishing a mathematical model that takes into account the offset and tilt of the levitated object is of great importance for accurately predicting the stability of the levitated system.
Disclosure of Invention
Aiming at the defects and improvement demands of the prior art, the invention provides a restoring force prediction method which simultaneously considers the eccentricity and the inclination of a suspended object. The method aims to effectively predict the restoring force generated by the system when the suspended object deflects, and provide theoretical guidance for parameter optimization for improving the stability of the system and precise motion control of the suspended object.
The invention realizes the above purpose through the following technical scheme:
the restoring force prediction method based on the eccentric and inclined states of the suspended object is applied to predicting the restoring force of an extrusion air film generated in near-field ultrasonic suspension established by an ultrasonic transducer on the suspended object, and comprises the following steps:
s1: obtaining the resonant frequency of an ultrasonic transducer, and carrying out measurement experiments according to the resonant frequency to obtain the amplitude distribution of a radiation disc in the ultrasonic transducer;
s2: constructing a restoring force prediction theoretical model based on the amplitude distribution and the eccentricity of the suspended object in the eccentric and inclined states, and obtaining a dimensionless expression of the extrusion air film under the prediction theoretical model;
s3: constructing a fluid motion control equation in the extrusion gas film based on a gas lubrication theory and a dimensionless expression of the extrusion gas film;
s4: introducing an eight-point discrete method to solve the fluid motion control equation to obtain dimensionless pressure distribution in the extrusion air film, constructing a moment balance equation and a suspension force balance equation when a suspended object is suspended above the extrusion air film based on the pressure distribution, and processing the moment balance equation and the suspension force balance equation by adopting a finite difference method and a spline interpolation method to obtain an inclination angle and a restoring force expression of the suspended object;
s5: and calculating the restoring force of the suspended object under different eccentricities by adopting the restoring force expression, generating a nonlinear relation between the restoring force and the corresponding eccentricity, using the actual restoring force and the corresponding eccentricity, which are tested and measured by the suspended object under different eccentricities, as a comparison group, and comparing the nonlinear relation with the comparison group to obtain the prediction precision.
As a further optimization scheme of the present invention, in step S1, the steps include:
s1.1: measuring the resonant frequency f of the ultrasonic transducer by adopting a precise impedance analyzer;
s1.2: setting an ultrasonic generator to output a sinusoidal voltage signal with the same resonant frequency f to the ultrasonic transducer;
s1.3: and measuring the amplitude distribution V (r) of the radiation disc of the ultrasonic transducer by adopting a scanning laser vibration meter.
As a further optimization scheme of the invention, in the step S2, when the eccentricity of the suspended object relative to the radiation disk is e, the overlapped part of the suspended object and the radiation disk is provided with a squeezing domain omega s The rest is set as a non-extrusion domain omega n In order to obtain an expression solving the intra-domain squeeze film, it is assumed that the suspended object has grooves of a size not to squeeze domain Ω n The dimensionless air film thickness expression is as follows:
wherein r and θ are coordinates of any point in the solution domain, α is an inclination angle of the floating object, h a For the distance between the center of the suspended object and the center of the radiation disk in the vertical direction, T=2pi ft is dimensionless time, h g For suspending the depth of the slot in the object.
As a further optimization scheme of the present invention, in step S3, the equation of fluid motion control in the squeeze film is:
in the formula, σ=12 μ a ωL 2 /p a h a 2 To characterize the number of extrusions, the pressure, μ, generated by the extrusion effect a Is the dynamic viscosity coefficient of air, Λ s =6μ a L 2 /p a h a 2 Sum lambda a =-ρ s L 2 /2p a Representing the degree of influence of the speed and acceleration of the relative movement of the suspended object on the pressure distribution, p=p/P a ,H=h/h a r=r/L and u=u/L are the dimensionless gas film pressure, gas film thickness, radial position and relative motion displacement, respectively, wherein ρ s To squeeze the density distribution of the film, p a For the ambient pressure, L is the radius of the radiating disk and the suspended object, p and h represent the actual air pressure and thickness in the squeeze film respectively,and->The speed and acceleration of the relative motion of the levitated object, respectively.
In step S4, based on the discontinuous thickness of the extruded air film caused by the offset of the suspended object, eight-point discrete method is introduced to solve the formula (2), so as to obtain the dimensionless air film pressure P in the extruded air film, the inclination caused by the offset satisfies the moment balance of the air film force distribution relative to the center of the suspended object, and the moment balance equation is:
in addition, the balance of the gravity of the suspended object and the levitation force generated by the squeeze film is expressed as:
wherein m and g are the mass of the suspended object and the local gravitational acceleration, respectively, and the position of the suspended object satisfying the formulas (3) and (4) is defined as a local equilibrium position, and the restoring force generated by the squeeze film is expressed as:
in the formula, the shear stress tau And τ zr Gas flow from the θ and r directions, θ 1 And obtaining the inclination angle and the restoring force of the suspended object by using a spline interpolation method for the angle of the connecting line of the micro-element rdrd theta and the center of the suspended object.
The invention has the beneficial effects that:
(1) In the invention, a theoretical calculation model which simultaneously considers the eccentricity and the inclination of the suspended object is established, and the model is closer to the actual motion condition, so that the method can accurately predict the return force values under different eccentricities.
(2) The prediction method provided by the invention realizes the prediction of restoring force by calculating the actual flow field distribution in the extrusion gas film, so that the influence mechanism of different designs and operation parameters on the stability of the suspension system can be analyzed theoretically, and the prediction method is further used for optimizing the design operation parameters and guiding the actual production.
(3) The prediction method provided by the invention can provide a theoretical basis for realizing flow field distribution control of accurate positioning.
Drawings
FIG. 1 is a flow chart of a near-field ultrasonic levitation restoring force prediction method in the present invention;
FIG. 2 is a schematic diagram of the structure of an ultrasonic transducer of the present invention;
FIG. 3 is a schematic diagram of near-field acoustic levitation in the present invention taking both levitating object eccentricity and tilt into account;
FIG. 4 is a schematic illustration of the structure of an imaginary suspension tray with slots in an embodiment of the invention;
FIG. 5 is a flow chart of the restoring force calculation in the present invention;
fig. 6 is a graph comparing the calculation results with the experimental results of the prediction method proposed in the present invention.
In fig. 2: 1. a radiating plate; 2. a horn; 3. a front cover plate; 4. a piezoelectric ceramic sheet; 5. a back cover plate; 6. an ultrasonic generator; 7. and (5) conducting wires.
Detailed Description
The invention will now be described in further detail with reference to the accompanying drawings, wherein it is to be understood that the following detailed description is for the purpose of illustration only and is not to be construed as limiting the scope of the invention, as various insubstantial modifications and adaptations of the invention to those skilled in the art may be made in light of the foregoing disclosure.
Example 1
As shown in fig. 1, the present application provides a restoring force prediction method based on eccentric and oblique states of a suspended object, where the method is applied to restoring force prediction of an extrusion air film generated in near-field ultrasonic suspension established by an ultrasonic transducer on the suspended object, and includes the following steps:
s1: obtaining the resonant frequency of an ultrasonic transducer, and carrying out measurement experiments according to the resonant frequency to obtain the amplitude distribution of a radiation disc in the ultrasonic transducer;
s2: constructing a restoring force prediction theoretical model based on amplitude distribution and the eccentricity of a suspended object in eccentric and inclined states, and obtaining a dimensionless expression of an extrusion air film under the prediction theoretical model;
s3: based on a gas lubrication theory and a dimensionless expression of the extrusion gas film, constructing a fluid motion control equation in the extrusion gas film;
s4: introducing an eight-point discrete method to solve a fluid motion control equation to obtain dimensionless pressure distribution in an extrusion air film, constructing a moment balance equation and a suspension force balance equation when a suspended object is suspended above the extrusion air film based on the pressure distribution, and processing the moment balance equation and the suspension force balance equation by adopting a finite difference method and a spline interpolation method to obtain an inclination angle and a restoring force expression of the suspended object;
s5: and calculating the restoring forces of a plurality of groups of different eccentricities by adopting a restoring force expression, generating nonlinear relations between the plurality of groups of restoring forces and corresponding eccentricities, using the actual restoring forces and the corresponding eccentricities which are tested under the different eccentricities as a comparison group, and comparing the nonlinear relations with the comparison group to obtain the prediction precision.
Further, in step S1, the steps include:
s1.1: measuring the resonant frequency f of the ultrasonic transducer by adopting a precise impedance analyzer;
s1.2: setting an ultrasonic generator to output a sine voltage signal with the same resonant frequency f to an ultrasonic transducer;
s1.3: the amplitude distribution V (r) of the radiating disk of the ultrasonic transducer is measured using a scanning laser vibrometer.
In this embodiment, as shown in fig. 2, the ultrasonic transducer mainly includes a piezoelectric ceramic plate 4, a back cover plate 5, a front cover plate 3, a horn 2, and a radiation disk 1. The resonant frequency f of the ultrasonic transducer is measured using a precision impedance analyzer. The ultrasonic generator 6 is arranged to output a sine voltage signal with the same frequency f, the sine voltage signal is applied to the piezoelectric ceramic plate 4 through a lead 7, the piezoelectric ceramic plate 4 generates vibration with the same frequency due to the inverse piezoelectric effect, the amplitude of the vibration is amplified by the amplitude transformer 2 and is transmitted to the radiation disc 1, and the amplitude distribution V (r) of the radiation disc is measured by the scanning laser vibration meter.
Further, as shown in FIG. 3, in step S2, the levitated object is levitated freely on the radiation tray 1, centers of the radiation tray 1 and the levitated object are denoted by O and O, respectively r When the offset distance of the levitated object relative to the radiation tray 1 is e, the levitated object is tilted by an angle α, and since the angle of tilt is very small, the projected area of the levitated object on the radiation tray is approximately equal to the area of the levitated object. Since the area of the radiating disk serves as a solving domain, it can be divided into a squeeze domain Ω s (portion where the levitated object and the radiating disk overlap) and non-extrusion domain Ω n (the remainder); to allow the film thickness to be expressed across the solution domain, a suspended object with grooves as shown in FIG. 4 was introduced. The size of the groove is the same as the size of the non-extrusion domain, and the depth of the groove is defined as h g . Then, the non-dimensional expression of the squeeze film H is:
wherein r and θ are coordinates of any point in the solution domain, α is an inclination angle of the floating object, h a For the distance between the center of the suspended object and the center of the radiation disk in the vertical direction, T=2pi ft is dimensionless time, h g For suspending the depth of the slot in the object.
In the present embodiment, since the inclination of the suspended object is considered, the squeeze film H satisfies at a minimum value near the center OWhere r=r/L is the dimensionless radial position and L is the radius of the radiating disc and the suspended object. Therefore, in this position the squeezing effect is most pronounced, with the air pressure variation satisfying +.>Wherein p=p/P a Is the pressure of a dimensionless air film, p a Is at ambient atmospheric pressure. Further, as shown in fig. 3, at the boundary of the squeeze film, the air pressure value should satisfy P (r=l) =1.
Further, based on the gas lubrication theory, in step S3, when the suspended object is considered to generate motion due to the restoring force, the fluid motion control equation in the squeeze film is as follows:
in the formula, σ=12 μ a ωL 2 /p a h a 2 To characterize the number of extrusions, the pressure, μ, generated by the extrusion effect a Is the dynamic viscosity coefficient of air, Λ s =6μ a L 2 /p a h a 2 Sum lambda a =-ρ s L 2 /2p a Representing the degree of influence of the velocity and acceleration of the suspended object movement on the pressure distribution, p=p/P a ,H=h/h a r=r/L and u=u/L are each a dimensionless gasFilm pressure, film thickness, radial position and relative motion displacement, where ρ s To squeeze the density distribution of the film, p a For the ambient pressure, L is the radius of the radiating disk and the suspended object, p and h represent the actual air pressure and thickness in the squeeze film respectively,and->The speed and acceleration of the relative motion of the levitated object, respectively.
Further, in step S4, based on the discontinuous thickness of the extruded air film caused by the offset of the suspended object, an eight-point discrete method and a newton-lavson method are introduced to solve the formula (2), so as to obtain the dimensionless air film pressure P in the extruded air film, the inclination caused by the offset satisfies the moment balance of the air film force distribution relative to the center of the suspended object, and the moment balance equation is:
in addition, the balance of the gravity of the suspended object and the levitation force generated by the squeeze film is expressed as:
wherein m and g are the mass of the suspended object and the local gravitational acceleration, respectively, and the position of the suspended object satisfying the formulas (3) and (4) is defined as a local equilibrium position, and the restoring force generated by the squeeze film is expressed as:
in the formula, the shear stress tau And τ zr Gas flow from squeeze film θ and r directions, θ 1 Is connected with the center of the suspended object by the infinitesimal rdrdrd thetaAngle of line.
Since the mass of the suspended object is measured by a precision balance instrument, the magnitude of the suspended force is known, and the suspended height and the inclination angle are unknown, the magnitude of the restoring force of the suspended object under the inclination condition can be obtained by using a two-time spline interpolation method as shown in fig. 5.
By varying the eccentricity e, a mathematical relationship between the restoring force and the eccentric distance can be established, and as shown in fig. 6, a nonlinear increase in restoring force can be found with an increase in eccentricity. When the eccentricity is smaller, the calculated result of the prediction method provided by the invention is relatively close to the calculated result without considering the inclination because the inclination angle is smaller at the moment. When the eccentricity is larger, the inclination angle is larger, and compared with the calculation result without considering inclination, the calculation result of the prediction method provided by the invention is more approximate to the experimental result. The closeness of the calculation result to the experimental result is described by R-square (determination coefficient), the value of which is approximately close to 1, which means that the higher the closeness is. The R-square value of the calculation result and the experimental result proposed by the invention is 0.9374, and the R-square value of the calculation result and the experimental result under no consideration of inclination is 0.8769. Therefore, the method realizes the prediction of restoring force by calculating the actual flow field distribution in the extrusion gas film, so that the influence mechanism of different designs and operation parameters on the stability of the suspension system can be analyzed theoretically, and the method is further used for optimizing the design operation parameters and guiding the actual production. In addition, the method can provide a theoretical basis for realizing flow field distribution control of accurate positioning.
The above embodiments are only for illustrating the technical solution of the present application, and are not limiting thereof; although the present application has been described in detail with reference to the foregoing embodiments, one of ordinary skill in the art will appreciate that: the technical scheme described in the foregoing embodiments can be modified or some technical features thereof can be replaced by equivalents; such modifications and substitutions do not drive the essence of the corresponding technical solutions to depart from the spirit and scope of the technical solutions of the embodiments of the present application.

Claims (3)

1.一种基于悬浮物体偏心和倾斜状态下的回复力预测方法,所述方法应用于超声换能器建立的近场超声悬浮中产生的挤压气膜对悬浮物体的回复力预测,其特征在于,包括如下步骤:1. A method for predicting the restoring force of a suspended object under eccentric and tilted conditions, the method being applied to predicting the restoring force of a suspended object by a squeezed air film generated in a near-field ultrasonic suspension established by an ultrasonic transducer, and characterized in that it comprises the following steps: S1:获取超声换能器的谐振频率,依据所述谐振频率进行测量实验得出超声换能器中辐射盘的振幅分布;S1: Obtaining the resonant frequency of the ultrasonic transducer, and performing a measurement experiment based on the resonant frequency to obtain the amplitude distribution of the radiation disk in the ultrasonic transducer; S2:基于所述振幅分布和悬浮物体在偏心与倾斜状态下的偏心距构建回复力预测理论模型,并得出该预测理论模型下的挤压气膜的无量纲表达式;S2: constructing a restoring force prediction theoretical model based on the amplitude distribution and the eccentricity of the suspended object in the eccentric and tilted state, and deriving a dimensionless expression of the squeezed air film under the prediction theoretical model; S3:基于气体润滑理论和所述挤压气膜的无量纲表达式,构建挤压气膜内流体运动控制方程;S3: Based on the gas lubrication theory and the dimensionless expression of the squeeze gas film, the control equation of the fluid motion in the squeeze gas film is constructed; S4:引入八点离散法对所述流体运动控制方程进行求解以获得挤压气膜内无量纲的压力分布,再基于所述压力分布构建悬浮物体位于挤压气膜上方悬浮时的力矩平衡方程和悬浮力平衡方程,采用有限差分法和样条插值法对所述力矩平衡方程及悬浮力平衡方程进行处理,以得出悬浮物体的倾斜角和回复力表达式;S4: Introducing an eight-point discretization method to solve the fluid motion control equation to obtain a dimensionless pressure distribution in the squeeze air film, and then constructing a moment balance equation and a suspension force balance equation when the suspended object is suspended above the squeeze air film based on the pressure distribution, and using a finite difference method and a spline interpolation method to process the moment balance equation and the suspension force balance equation to obtain an expression for the inclination angle and restoring force of the suspended object; S5:采用所述回复力表达式计算悬浮物体在不同偏心距下的回复力,并生成回复力和对应偏心距的非线性关系,再将悬浮物体在不同偏心距下试验测得的实际回复力与对应偏心距作为对照组,将所述非线性关系与所述对照组进行比较以得到预测精度;S5: using the restoring force expression to calculate the restoring force of the suspended object at different eccentricities, and generating a nonlinear relationship between the restoring force and the corresponding eccentricity, then using the actual restoring force and the corresponding eccentricity measured by the suspended object at different eccentricities as a control group, and comparing the nonlinear relationship with the control group to obtain the prediction accuracy; 步骤S1中,步骤包括:In step S1, the steps include: S1.1:采用精密阻抗分析仪对超声换能器的谐振频率f进行测量;S1.1: Use a precision impedance analyzer to measure the resonant frequency f of the ultrasonic transducer; S1.2:设置超声发生器输出与所述谐振频率f相同的正弦电压信号至所述超声换能器中;S1.2: Setting the ultrasonic generator to output a sinusoidal voltage signal having the same resonant frequency f to the ultrasonic transducer; S1.3:采用扫描激光测振仪对所述超声换能器的辐射盘的振幅分布V(r)进行测量;S1.3: using a scanning laser vibrometer to measure the amplitude distribution V(r) of the radiation disk of the ultrasonic transducer; 步骤S2中,悬浮物体相对于辐射盘的偏心距为e时,悬浮物体与辐射盘重叠部分设有挤压域 ,其余部分设为非挤压域/>,为了获得求解域内挤压气膜的表达式,假设所述悬浮物体带有槽,槽的大小与非挤压域/>的大小相同,则无量纲气膜厚度H表达式为: In step S2, when the eccentricity of the suspended object relative to the radiation disk is e , the overlapping part of the suspended object and the radiation disk is provided with an extrusion zone , the rest is set as non-extrusion domain/> In order to obtain the expression for squeezing the air film in the solution domain, it is assumed that the suspended object has a groove, and the size of the groove is similar to that of the non-squeezed domain/> The size of is the same, then the expression of dimensionless gas film thickness H is: 式中,V(r)为辐射盘的振幅分布,rθ为求解域内任意一点的坐标,α为悬浮物体的倾斜角,为悬浮物体中心与辐射盘中心在竖直方向的距离,/>为无量纲的时间,/>悬浮物体中槽的深度。Where V(r) is the amplitude distribution of the radiation disk, r and θ are the coordinates of any point in the solution domain, α is the inclination angle of the suspended object, is the vertical distance between the center of the suspended object and the center of the radiation disk, /> is the dimensionless time, /> The depth of the groove in the suspended object. 2.根据权利要求1所述的一种基于悬浮物体偏心和倾斜状态下的回复力预测方法,其特征在于:步骤S3中,挤压气膜内流体运动控制方程为: 2. The method for predicting the restoring force of a suspended object under eccentric and tilted conditions according to claim 1, characterized in that: in step S3, the control equation for the motion of the fluid in the squeezed air film is: 式中,为挤压数,表征由挤压效应产生的压力,/>为空气的动力粘度系数,/>分别代表悬浮物体相对运动的速度和加速度对压力分布影响的程度,/>, />,/>和/>分别为无量纲气膜压力、气膜厚度、径向位置和相对运动位移,其中,/>为挤压气膜的密度分布,/>为环境压力,L为辐射盘和悬浮物体的半径,ph分别代表挤压气膜内实际气压和厚度,/>和/>分别为悬浮物体相对运动的速度和加速度。 In the formula, is the extrusion number, which characterizes the pressure generated by the extrusion effect,/> is the dynamic viscosity coefficient of air, /> They represent the influence of the relative motion speed and acceleration of the suspended object on the pressure distribution, respectively. , /> ,/> and/> are dimensionless air film pressure, air film thickness, radial position and relative motion displacement, respectively, where,/> is the density distribution of the squeezed air film, /> is the ambient pressure, L is the radius of the radiation disk and the suspended object, p and h represent the actual pressure and thickness of the squeezed air film, respectively. and/> are the velocity and acceleration of the relative motion of the suspended object respectively. 3.根据权利要求2所述的一种基于悬浮物体偏心和倾斜状态下的回复力预测方法,其特征在于:步骤S4中,基于悬浮物体偏移导致的挤压气膜厚度不连续,引入八点离散法来对公式(2)进行求解,进而获得挤压气膜内无量纲气膜压力P,由偏移导致的倾斜满足气膜力分布相对于悬浮物体中心的力矩平衡,力矩平衡方程为: 3. A restoring force prediction method based on the eccentricity and tilt state of a suspended object according to claim 2, characterized in that: in step S4, based on the discontinuity of the squeezed air film thickness caused by the displacement of the suspended object, an eight-point discretization method is introduced to solve formula (2), thereby obtaining the dimensionless air film pressure P in the squeezed air film, and the tilt caused by the displacement satisfies the moment balance of the air film force distribution relative to the center of the suspended object, and the moment balance equation is: 此外,悬浮物体的重力和挤压气膜产生的悬浮力平衡,表示为: In addition, the gravity of the suspended object and the suspension force generated by the squeezed air film are balanced, which can be expressed as: 其中mg分别为悬浮物体的质量和当地的重力加速度,此时满足公式(3)和(4)的悬浮物体位置定义为局部平衡位置,挤压气膜产生的回复力表示为: Where m and g are the mass of the suspended object and the local gravitational acceleration, respectively. The position of the suspended object that satisfies formulas (3) and (4) is defined as the local equilibrium position. The restoring force generated by squeezing the air film is expressed as: 式中,剪切应力和/>源自挤压气膜θr方向的气体流动,/>为微元/>与悬浮物体中心连线的角度,利用样条插值法,以获得悬浮物体的倾斜角和回复力大小。In the formula, shear stress and/> Gas flow in the θ and r directions originating from the squeeze film,/> For microelement/> The angle of the line connecting the center of the suspended object is calculated using the spline interpolation method to obtain the inclination angle and the magnitude of the restoring force of the suspended object.
CN202310872125.2A 2023-07-17 2023-07-17 Restoring force prediction method based on eccentric and inclined states of suspended object Active CN116644644B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310872125.2A CN116644644B (en) 2023-07-17 2023-07-17 Restoring force prediction method based on eccentric and inclined states of suspended object

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310872125.2A CN116644644B (en) 2023-07-17 2023-07-17 Restoring force prediction method based on eccentric and inclined states of suspended object

Publications (2)

Publication Number Publication Date
CN116644644A CN116644644A (en) 2023-08-25
CN116644644B true CN116644644B (en) 2024-04-16

Family

ID=87619199

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310872125.2A Active CN116644644B (en) 2023-07-17 2023-07-17 Restoring force prediction method based on eccentric and inclined states of suspended object

Country Status (1)

Country Link
CN (1) CN116644644B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN118569038B (en) * 2024-06-05 2025-02-28 安徽大学 A method for predicting and controlling the motion frequency of an ultrasonically suspended object under eccentricity
CN119167459A (en) * 2024-09-14 2024-12-20 安徽大学 A method for predicting the levitation accuracy of near-field ultrasonic levitation objects

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008227480A (en) * 2007-02-14 2008-09-25 Institute Of National Colleges Of Technology Japan Noncontact holding device and noncontact conveying device
WO2011055482A1 (en) * 2009-11-06 2011-05-12 村田機械株式会社 Contact-less holding device, transfer device, and contact-less holding method
CN102567563A (en) * 2011-10-18 2012-07-11 上海交通大学 Accurately quantitative calculation method for near-field ultrasonic levitation force
WO2015114634A1 (en) * 2014-02-02 2015-08-06 Technion Research & Development Foundation Limited. Method and system for non-contact levitation
CN109780056A (en) * 2019-01-15 2019-05-21 浙江大学 Ultrasonic suspended radial bearing capable of self-adaptively stabilizing the position of the inner ring of the bearing
CN112069710A (en) * 2020-09-07 2020-12-11 电子科技大学 Prediction method for self-excited vibration of gas static pressure main shaft
GB202107416D0 (en) * 2021-05-25 2021-07-07 Alker Graham Prediction of far-field radiated acoustic fields in a fluid using near-field acoustic measurements
CN113569450A (en) * 2021-07-15 2021-10-29 东南大学 Method for estimating and controlling suspension and residence of liquid drops
CN116317684A (en) * 2023-01-31 2023-06-23 安徽大学 Non-contact motor based on standing wave type near field ultrasonic levitation and static pressure gas levitation coupling

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008227480A (en) * 2007-02-14 2008-09-25 Institute Of National Colleges Of Technology Japan Noncontact holding device and noncontact conveying device
WO2011055482A1 (en) * 2009-11-06 2011-05-12 村田機械株式会社 Contact-less holding device, transfer device, and contact-less holding method
CN102567563A (en) * 2011-10-18 2012-07-11 上海交通大学 Accurately quantitative calculation method for near-field ultrasonic levitation force
WO2015114634A1 (en) * 2014-02-02 2015-08-06 Technion Research & Development Foundation Limited. Method and system for non-contact levitation
CN109780056A (en) * 2019-01-15 2019-05-21 浙江大学 Ultrasonic suspended radial bearing capable of self-adaptively stabilizing the position of the inner ring of the bearing
CN112069710A (en) * 2020-09-07 2020-12-11 电子科技大学 Prediction method for self-excited vibration of gas static pressure main shaft
GB202107416D0 (en) * 2021-05-25 2021-07-07 Alker Graham Prediction of far-field radiated acoustic fields in a fluid using near-field acoustic measurements
CN113569450A (en) * 2021-07-15 2021-10-29 东南大学 Method for estimating and controlling suspension and residence of liquid drops
CN116317684A (en) * 2023-01-31 2023-06-23 安徽大学 Non-contact motor based on standing wave type near field ultrasonic levitation and static pressure gas levitation coupling

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Self-floating and self-rotating non-contact ultrasonic motor with single active vibrator;Minghui Shi,Shujie Chen;Tribology International;20230228;第180卷;1-13 *
基于近场超声挤压的非接触支撑与传输系统理论研究;刘远远;硕士电子期刊;20180715(2018年第07期);C029-3 *

Also Published As

Publication number Publication date
CN116644644A (en) 2023-08-25

Similar Documents

Publication Publication Date Title
CN116644644B (en) Restoring force prediction method based on eccentric and inclined states of suspended object
Pandey et al. Effect of flexural modes on squeeze film damping in MEMS cantilever resonators
Cattafesta III et al. Development of piezoelectric actuators for active flow control
WO2013111608A1 (en) Viscoelasticity measurement method and viscoelasticity measurement device
Go et al. Heat transfer enhancement using flow-induced vibration of a microfin array
Liu et al. Stabilizing near-field acoustic levitation: Investigation of non-linear restoring force generated by asymmetric gas squeeze film
Feng et al. A novel squeeze-film air bearing with flexure pivot-tilting pads: numerical analysis and measurement
Feng et al. Numerical analysis of the transportation characteristics of a self-running sliding stage based on near-field acoustic levitation
Beltman et al. Air loads on a rigid plate oscillating normal to a fixed surface
Wang et al. Stiffness modeling for near field acoustic levitation bearings
JPH11211705A (en) Method and apparatus for separately measuring density and viscosity of liquid
Chandrasekaran et al. Dynamic calibration technique for thermal shear-stress sensors with mean flow
Higashino et al. Self-excited vibrational cantilever-type viscometer driven by piezo-actuator
Stolarski Acoustic levitation-A novel alternative to traditional lubrication of contacting surfaces
Koyama et al. An ultrasonic air pump using an acoustic traveling wave along a small air gap
JP5831904B2 (en) Viscoelasticity measuring method and viscoelasticity measuring device
Almurshedi et al. Modelling influence of Poisson's contraction on squeeze film levitation of planar objects
Koyama et al. A self-running standing wave-type bidirectional slider for the ultrasonically levitated thin linear stage
CN114113315B (en) A method for measuring gas viscosity based on acoustic resonance energy dissipation
Hosaka et al. Coupled vibration of microcantilever array induced by airflow force
CN118569038B (en) A method for predicting and controlling the motion frequency of an ultrasonically suspended object under eccentricity
Mahajan et al. Experimental and analytical investigation of a dynamic gas squeeze film bearing including asperity contact effects
Ma et al. Influence of the disturbance frequency on the dynamic characteristics of an ultrasonic gas squeeze film
Stolarski et al. Acoustic journal bearing with changeable geometry and built-in flexibility
Hu et al. Characterization of the acoustic field generated by a horn shaped ultrasonic transducer

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant