WO2021017828A1 - 一种高频驻波振幅分布的高精度测量方法 - Google Patents
一种高频驻波振幅分布的高精度测量方法 Download PDFInfo
- Publication number
- WO2021017828A1 WO2021017828A1 PCT/CN2020/102064 CN2020102064W WO2021017828A1 WO 2021017828 A1 WO2021017828 A1 WO 2021017828A1 CN 2020102064 W CN2020102064 W CN 2020102064W WO 2021017828 A1 WO2021017828 A1 WO 2021017828A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- spot
- amplitude distribution
- deflection
- standing wave
- feature
- 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.)
- Ceased
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H9/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means
Definitions
- the invention belongs to the technical field of optical engineering, and is specifically a method for measuring the amplitude distribution of a high-frequency standing wave.
- the hemispherical resonant gyroscope and other components achieve performance through the resonant standing wave generated by excitation, so the measurement of its vibration characteristics is an important guarantee for mechanism analysis and improvement of component performance.
- the measured component is plated with a metal film, and the formed electrode and the signal detector form a capacitance, so the vibration measurement can be realized by detecting the change of the capacitance characteristic [1].
- the average distance change at the electrode and the spatial distribution of the amplitude cannot be obtained.
- the purpose of the present invention is to provide a high-precision measurement method for the amplitude distribution of the standing wave of a smooth surface element, so as to facilitate the analysis of the vibration mode and vibration characteristics.
- the high-precision measurement method for the amplitude distribution of the standing wave of the smooth surface element provided by the present invention is based on the deflection technology, and the specific steps are as follows:
- (u 0 v 0 ) represents the entire
- the center of gravity coordinates of the feature spot, (u(i,j),v(i,j)) represents the coordinate of any pixel in the feature spot;
- step (4) Perform eigenvalue decomposition on the normal matrix M in step (4), where the eigenvector corresponding to the largest eigenvalue is the extension direction of the feature spot [5]; draw a straight line perpendicular to the extension direction from the center of the feature spot , The width of the characteristic spot along the straight line is the original circle spot diameter d;
- the normal deflection range of the measuring point caused by vibration can be obtained.
- the components in the x and y directions are written as s x , s y , and the z-direction component is normalized;
- the modal method is used for integral reconstruction to obtain the amplitude distribution z [6] in the whole area, that is, to solve the partial differential equation, so that the objective function is minimized:
- D x and D y respectively represent the differential matrix along the x and y directions, and the x and y components of the normal deflection can be approximated and optimized to obtain z.
- the system of the present invention has simple structure, high sensitivity and strong anti-interference ability, can measure the standing wave amplitude distribution in a local area of the component, and is of great significance for component vibration mechanism and characteristic analysis and component resonance performance guarantee.
- Figure 1 is a schematic diagram of the deflection measurement amplitude.
- Figure 2 is a schematic diagram of the formation of blurred spots in the deflection measurement.
- Figure 3 shows the binarization pattern displayed on the screen.
- Figure 4 shows the blurred image obtained by shooting.
- Figure 5 shows the recognition contour of the blurred spot.
- Figure 6 shows the extension direction of the blur spot and the calculation of the yaw distance.
- Figure 7 shows the normal x-direction yaw range.
- Figure 8 shows the normal y-direction yaw range.
- Figure 9 shows the actual solution of the standing wave amplitude distribution.
- Embodiment 1 In this embodiment, the measured component is a plane with a diameter of 70 mm, and the vibration is excited from the lower surface, and the frequency is 4500 Hz.
- the focal length of the camera is 55mm, and the size of a single pixel is 8 ⁇ m ⁇ 8 ⁇ m.
- the screen adopts ipad mini, the angle between the center point of the two and the excitation point is 45°, the distance between the center of the screen and the excitation point of the component is 172mm, and the distance between the optical center of the camera and the excitation point is 240mm, as shown in Figure 1. .
- the standing wave of the component under test causes the local normal to change, so the reflection optics will swing accordingly, and the corresponding image will form blurred spots, as shown in Figure 2.
- the round spot pattern is displayed on the screen, as shown in Figure 3.
- the camera exposure time is set to 0.05 seconds, and the collected images are shown in Figure 4.
- the Sobel operator is used to identify the boundary of the feature spot; for clarity, we only draw the identification boundary of a feature spot in the upper left corner in FIG. 5.
- the normal matrix of all pixels in the contour is subjected to eigenvalue decomposition, and the eigenvector corresponding to the largest eigenvalue is the AB direction shown in FIG. 6.
- the normal yaw range of the component can be obtained, and its z component is normalized to obtain the lateral component (s x , s y ), respectively As shown in Figure 7, Figure 8.
- the amplitude distribution obtained after reconstruction and fitting using the modal method is shown in Figure 9. Therefore, the method of the present invention can accurately obtain the amplitude distribution of the measurement area.
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Holo Graphy (AREA)
- Gyroscopes (AREA)
Abstract
Description
Claims (1)
- 一种高频驻波振幅分布精密测量方法,基于偏折技术的,其特征在于,具体步骤如下:(1)搭建偏折测量光路,投影屏幕和相机相对于被测元件左右对称放置;(2)在屏幕上显示二值化圆斑图样,经过振动元件反射后在相机成像;(3)采用Sobel微分算子对图像进行差分计算,对算子进行二值化处理,识别所采集图像中特征斑的边界,以灰度w(u,v)为权重计算斑点中心:(4)对特征斑中所有点坐标构造正规矩阵:M=∑w(i,j)δ(i,j) Tδ(i,j),其中,δ(i,j)=(u(i,j)-u 0,v(i,j)-v 0)表示特征斑中每个像素的坐标偏差,(u 0 v 0)表示整个特征斑的重心坐标,(u(i,j),v(i,j))表示特征斑中任意像素的坐标;(5)对步骤(4)中的正规矩阵M进行特征值分解,其中最大特征值对应的特征向量即为特征斑的延伸方向;从特征斑重心画出一条垂直于延伸方向的直线,特征斑沿该直线方向的宽度即为原始圆斑直径d;(6)从特征斑重心沿其延伸方向画一条直线,与特征斑边界的两个交点记为A、B;线段AB长度与圆斑直径d之差的一半即为该点极限振幅引起的成像偏差;(7)利用偏折标定的几何参数便可得到振动引起的测量点法向偏摆范围,将其x,y方向的分量分别写作s x,s y,并将z向分量归一化;(8)采用modal法进行积分重构得到全区域振幅分布z,也即解偏微分方程,使得目标函数最小:E(z)=‖zD x-s x‖ 2+||D yz-s y|| 2 (2)其中,D x、D y分别表示沿x、y方向的微分矩阵,和法向偏摆的x、y方向分量进行逼近优化,解得z;求解后振幅分布整体偏移,使得驻点振幅为零。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910685221.X | 2019-07-27 | ||
| CN201910685221.XA CN110501063B (zh) | 2019-07-27 | 2019-07-27 | 一种高频驻波振幅分布的高精度测量方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| WO2021017828A1 true WO2021017828A1 (zh) | 2021-02-04 |
Family
ID=68587576
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| PCT/CN2020/102064 Ceased WO2021017828A1 (zh) | 2019-07-27 | 2020-07-15 | 一种高频驻波振幅分布的高精度测量方法 |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN110501063B (zh) |
| WO (1) | WO2021017828A1 (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110501063B (zh) * | 2019-07-27 | 2021-06-04 | 复旦大学 | 一种高频驻波振幅分布的高精度测量方法 |
| CN113686360B (zh) * | 2021-08-27 | 2023-05-26 | 中国科学院西安光学精密机械研究所 | 一种半球谐振子驻波漂移的全局测量方法及系统 |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102679902A (zh) * | 2012-05-24 | 2012-09-19 | 天津大学 | 薄平板结构共振模态分析系统及其使用方法 |
| CN103424175A (zh) * | 2012-05-21 | 2013-12-04 | 成都星聚能科技有限公司 | 基于激光调制的桥梁振幅实时监测装置 |
| CN103438980A (zh) * | 2013-09-01 | 2013-12-11 | 北京航空航天大学 | 一种基于线阵ccd和线状红外激光的液体表面波探测方法及装置 |
| CN104776907A (zh) * | 2015-04-30 | 2015-07-15 | 杭州电子科技大学 | 一种基于多点激光散斑极值跟踪的振动检测方法 |
| CN105043524A (zh) * | 2015-04-30 | 2015-11-11 | 杭州电子科技大学 | 一种基于激光散斑离焦成像的振动检测方法 |
| JP2017003397A (ja) * | 2015-06-09 | 2017-01-05 | 国立大学法人 新潟大学 | 平面振動計測装置及び平面振動計測方法 |
| CN106768280A (zh) * | 2017-02-28 | 2017-05-31 | 北京航空航天大学 | 一种基于多波长无透镜傅里叶变换数字全息的振动检测装置 |
| CN107300413A (zh) * | 2017-07-13 | 2017-10-27 | 天津大学 | 基于条纹投影的远程低频振动测量方法 |
| CN107560709A (zh) * | 2017-08-11 | 2018-01-09 | 维沃移动通信有限公司 | 一种振幅测量方法及测量设备 |
| WO2018134932A1 (ja) * | 2017-01-19 | 2018-07-26 | 富士通株式会社 | 振幅測定装置、及び、振幅測定装置における振動測定方法 |
| CN110501063A (zh) * | 2019-07-27 | 2019-11-26 | 复旦大学 | 一种高频驻波振幅分布的高精度测量方法 |
Family Cites Families (22)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO1999029139A2 (de) * | 1997-12-03 | 1999-06-10 | Sennheiser Electronic Gmbh & Co. Kg | Optischer schallsensor, insbesondere mikrofon, mit zusätzlicher spiegelungseinrichtung |
| JP2007101394A (ja) * | 2005-10-05 | 2007-04-19 | Shinetsu Film Kk | 結晶面方位の簡易判定方法及びその装置 |
| JP4976143B2 (ja) * | 2006-09-29 | 2012-07-18 | 古河電気工業株式会社 | 衝撃振動検知装置 |
| US8570530B2 (en) * | 2009-06-03 | 2013-10-29 | Carestream Health, Inc. | Apparatus for dental surface shape and shade imaging |
| US10786158B2 (en) * | 2009-12-10 | 2020-09-29 | The Board Of Regents Of The University Of Texas System | Optoacoustic / photoacoustic / acoustic imaging system using probe beam deflection |
| JP5912712B2 (ja) * | 2012-03-21 | 2016-04-27 | スタンレー電気株式会社 | 照明用光学系 |
| CN102980746B (zh) * | 2012-11-20 | 2015-04-15 | 天津大学 | 基于数字全息的晶体光学参数测量方法 |
| CN103278935A (zh) * | 2013-05-20 | 2013-09-04 | 复旦大学 | 一种应用自由曲面分色片进行像差补偿的分色方法 |
| CN103727894A (zh) * | 2014-01-15 | 2014-04-16 | 唐山英莱科技有限公司 | 一种基于线结构光折射成像的透明体三维轮廓检测系统 |
| EP3289430B1 (en) * | 2015-04-27 | 2019-10-23 | Snap-Aid Patents Ltd. | Estimating and using relative head pose and camera field-of-view |
| CN105115908B (zh) * | 2015-08-20 | 2017-11-21 | 哈尔滨理工大学 | 金属锈斑最优光谱波段选取视觉检测装置和方法 |
| CN105387819B (zh) * | 2015-12-27 | 2017-10-13 | 河北工业大学 | 基于条纹反射法的反光物体三维形貌测量方法及装置 |
| CN105571518B (zh) * | 2016-01-19 | 2017-10-17 | 大连理工大学 | 基于折射图像偏差的三维信息视觉测量方法 |
| CN105783775B (zh) * | 2016-04-21 | 2018-05-11 | 清华大学 | 一种镜面及类镜面物体表面形貌测量装置与方法 |
| US10073026B2 (en) * | 2016-07-05 | 2018-09-11 | The United States Of America, As Represented By The Secretary Of Commerce | Optical particle sorter |
| CN107339941A (zh) * | 2017-07-21 | 2017-11-10 | 复旦大学 | 一种基于双频激光干涉原理的精确位移监测系统 |
| CN108413893B (zh) * | 2018-03-12 | 2020-06-05 | 四川大学 | 一种散斑偏折术检测平面元件面形的方法与装置 |
| CN108426539A (zh) * | 2018-03-12 | 2018-08-21 | 四川大学 | 一种基于双频正交光栅投影检测反射元件面形的方法 |
| CN108593088B (zh) * | 2018-04-25 | 2020-09-08 | 重庆医科大学 | 一种表征球形聚焦超声场的聚焦纹影系统 |
| CN108775874A (zh) * | 2018-06-04 | 2018-11-09 | 西安工业大学 | 一种应用随机相移技术的阴影三维轮廓测量装置和方法 |
| CN109186496B (zh) * | 2018-10-18 | 2020-05-01 | 淮阴师范学院 | 一种基于移动最小二乘法的三维面形测量方法 |
| CN109870129A (zh) * | 2019-03-25 | 2019-06-11 | 中国计量大学 | 一种基于相位偏折原理的晶圆表面粗糙度检测装置 |
-
2019
- 2019-07-27 CN CN201910685221.XA patent/CN110501063B/zh active Active
-
2020
- 2020-07-15 WO PCT/CN2020/102064 patent/WO2021017828A1/zh not_active Ceased
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103424175A (zh) * | 2012-05-21 | 2013-12-04 | 成都星聚能科技有限公司 | 基于激光调制的桥梁振幅实时监测装置 |
| CN102679902A (zh) * | 2012-05-24 | 2012-09-19 | 天津大学 | 薄平板结构共振模态分析系统及其使用方法 |
| CN103438980A (zh) * | 2013-09-01 | 2013-12-11 | 北京航空航天大学 | 一种基于线阵ccd和线状红外激光的液体表面波探测方法及装置 |
| CN104776907A (zh) * | 2015-04-30 | 2015-07-15 | 杭州电子科技大学 | 一种基于多点激光散斑极值跟踪的振动检测方法 |
| CN105043524A (zh) * | 2015-04-30 | 2015-11-11 | 杭州电子科技大学 | 一种基于激光散斑离焦成像的振动检测方法 |
| JP2017003397A (ja) * | 2015-06-09 | 2017-01-05 | 国立大学法人 新潟大学 | 平面振動計測装置及び平面振動計測方法 |
| WO2018134932A1 (ja) * | 2017-01-19 | 2018-07-26 | 富士通株式会社 | 振幅測定装置、及び、振幅測定装置における振動測定方法 |
| CN106768280A (zh) * | 2017-02-28 | 2017-05-31 | 北京航空航天大学 | 一种基于多波长无透镜傅里叶变换数字全息的振动检测装置 |
| CN107300413A (zh) * | 2017-07-13 | 2017-10-27 | 天津大学 | 基于条纹投影的远程低频振动测量方法 |
| CN107560709A (zh) * | 2017-08-11 | 2018-01-09 | 维沃移动通信有限公司 | 一种振幅测量方法及测量设备 |
| CN110501063A (zh) * | 2019-07-27 | 2019-11-26 | 复旦大学 | 一种高频驻波振幅分布的高精度测量方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110501063A (zh) | 2019-11-26 |
| CN110501063B (zh) | 2021-06-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113237435B (zh) | 一种高反光表面三维视觉测量系统及方法 | |
| Fu et al. | Non-contact optical dynamic measurements at different ranges: a review | |
| CN107560564B (zh) | 一种自由曲面检测方法及系统 | |
| CN109307480B (zh) | 一种透射元件多表面面形检测方法 | |
| KR20110106823A (ko) | 비구면체 측정 방법 및 장치 | |
| CN106017863A (zh) | 检测非球面的相位测量偏折方法 | |
| CN109870129A (zh) | 一种基于相位偏折原理的晶圆表面粗糙度检测装置 | |
| Akhtar et al. | An efficient automotive paint defect detection system | |
| WO2021017828A1 (zh) | 一种高频驻波振幅分布的高精度测量方法 | |
| CN120411100B (zh) | 基于相位偏折术的屏幕橘皮检测方法及系统 | |
| Guo et al. | Research on sub-pixel accuracy flange disk dimension measurement based on machine vision | |
| CN114674244A (zh) | 一种同轴正入射散斑偏折术测量方法和装置 | |
| CN112683382B (zh) | 基于单目视觉的结构三维振动测量系统及方法 | |
| JP2000111490A (ja) | 塗装面の検出装置 | |
| CN113048912A (zh) | 投影仪的标定系统及方法 | |
| Zhang et al. | High-accuracy three-dimensional shape measurement of micro solder paste and printed circuits based on digital image correlation | |
| CN119850489A (zh) | 一种高精度全域准静态应变视觉测量方法 | |
| CN218411080U (zh) | 基于机器视觉的平面平晶自动检测装置 | |
| CN205383999U (zh) | 物体平整度光学检测系统 | |
| Tay et al. | New method for measuring dynamic response of small components by fringe projection | |
| CN107869964A (zh) | 非球面面形的检测方法和装置 | |
| Liu et al. | Microscopic 3D Surface Imaging With Annular Spectrum Sampling Parallel Single‐Pixel Imaging: Resistant to Global Illumination | |
| Onishi et al. | Theoretical and experimental guideline of optimum design of defect-inspection apparatus for transparent material using phase-shift illumination approach | |
| Gao et al. | Full-field deformation measurement of 2D-DIC based on virtual camera array in small field of view | |
| KR101379677B1 (ko) | 구면파 발생 간섭계를 이용한 비구면렌즈의 편심 측정시스템 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| 121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 20847646 Country of ref document: EP Kind code of ref document: A1 |
|
| NENP | Non-entry into the national phase |
Ref country code: DE |
|
| 122 | Ep: pct application non-entry in european phase |
Ref document number: 20847646 Country of ref document: EP Kind code of ref document: A1 |
|
| 32PN | Ep: public notification in the ep bulletin as address of the adressee cannot be established |
Free format text: NOTING OF LOSS OF RIGHTS PURSUANT TO RULE 112(1) EPC (EPO FORM 1205A DATED 23.05.2023) |
|
| 122 | Ep: pct application non-entry in european phase |
Ref document number: 20847646 Country of ref document: EP Kind code of ref document: A1 |