WO2021017828A1 - High-precision method for measuring high-frequency standing wave amplitude distribution - Google Patents
High-precision method for measuring high-frequency standing wave amplitude distribution 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 15
- 238000005259 measurement Methods 0.000 claims abstract description 9
- 239000011159 matrix material Substances 0.000 claims description 7
- 230000003287 optical effect Effects 0.000 claims description 5
- 230000005484 gravity Effects 0.000 claims description 4
- 238000000691 measurement method Methods 0.000 claims description 4
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- 238000003384 imaging method Methods 0.000 claims description 3
- 238000012360 testing method Methods 0.000 claims description 3
- 101100445834 Drosophila melanogaster E(z) gene Proteins 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 claims description 2
- 238000004458 analytical method Methods 0.000 abstract description 4
- 230000007246 mechanism Effects 0.000 abstract description 3
- 230000035945 sensitivity Effects 0.000 abstract description 2
- 230000010355 oscillation Effects 0.000 abstract 1
- 230000005284 excitation Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000005210 holographic interferometry Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000008569 process Effects 0.000 description 1
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.
Abstract
Description
Claims (1)
- 一种高频驻波振幅分布精密测量方法,基于偏折技术的,其特征在于,具体步骤如下:A precision measurement method for the amplitude distribution of high-frequency standing waves, based on deflection technology, is characterized in that the specific steps are as follows:(1)搭建偏折测量光路,投影屏幕和相机相对于被测元件左右对称放置;(1) Set up the deflection measurement optical path, and place the projection screen and camera symmetrically with respect to the component under test;(2)在屏幕上显示二值化圆斑图样,经过振动元件反射后在相机成像;(2) Display the binarized round spot pattern on the screen, which will be imaged on the camera after being reflected by the vibrating element;(3)采用Sobel微分算子对图像进行差分计算,对算子进行二值化处理,识别所采集图像中特征斑的边界,以灰度w(u,v)为权重计算斑点中心:(3) The Sobel differential operator is used to calculate the difference of the image, and the operator is binarized to identify the boundary of the characteristic spot in the collected image, and calculate the spot center with the gray level w(u,v) as the weight:(4)对特征斑中所有点坐标构造正规矩阵:(4) Construct a normal matrix for the coordinates of all points in the feature spot:M=∑w(i,j)δ(i,j) Tδ(i,j), 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))表示特征斑中任意像素的坐标; Among them, δ(i,j)=(u(i,j)-u 0 ,v(i,j)-v 0 ) represents the coordinate deviation of each pixel in the feature spot, (u 0 v 0 ) represents the entire feature The center of gravity coordinates of the spot, (u(i,j),v(i,j)) represents the coordinates of any pixel in the feature spot;(5)对步骤(4)中的正规矩阵M进行特征值分解,其中最大特征值对应的特征向量即为特征斑的延伸方向;从特征斑重心画出一条垂直于延伸方向的直线,特征斑沿该直线方向的宽度即为原始圆斑直径d;(5) 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; draw a straight line perpendicular to the extension direction from the center of the feature spot, the feature spot The width along the straight line is the original circle diameter d;(6)从特征斑重心沿其延伸方向画一条直线,与特征斑边界的两个交点记为A、B;线段AB长度与圆斑直径d之差的一半即为该点极限振幅引起的成像偏差;(6) Draw a straight line from the center of gravity of the characteristic spot along its extension direction, and mark the two intersection points with the boundary of the characteristic spot as A and B; the half of the difference between the length of the line segment AB and the diameter d of the spot is the imaging caused by the limit amplitude at this point deviation;(7)利用偏折标定的几何参数便可得到振动引起的测量点法向偏摆范围,将其x,y方向的分量分别写作s x,s y,并将z向分量归一化; (7) The normal deflection range of the measuring point caused by vibration can be obtained by using the geometric parameters calibrated by deflection. The components in the x and y directions are respectively written as s x and s y , and the z-direction component is normalized;(8)采用modal法进行积分重构得到全区域振幅分布z,也即解偏微分方程,使得目标函数最小:(8) Use the modal method to perform integral reconstruction to obtain the amplitude distribution z in the whole region, that is, to solve the partial differential equation, so that the objective function is minimized:E(z)=‖zD x-s x‖ 2+||D yz-s y|| 2 (2) E(z)=‖zD x -s x ‖ 2 +||D y zs y || 2 (2)其中,D x、D y分别表示沿x、y方向的微分矩阵,和法向偏摆的x、y方向分量进行逼近优化,解得z;求解后振幅分布整体偏移,使得驻点振幅为零。 Among them, 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 yaw are approximated and optimized to obtain z; after the solution, the overall amplitude distribution is shifted so that the stationary point amplitude is zero.
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910685221.XA CN110501063B (en) | 2019-07-27 | 2019-07-27 | High-precision measurement method for high-frequency standing wave amplitude distribution |
CN201910685221.X | 2019-07-27 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2021017828A1 true WO2021017828A1 (en) | 2021-02-04 |
Family
ID=68587576
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/CN2020/102064 WO2021017828A1 (en) | 2019-07-27 | 2020-07-15 | High-precision method for measuring high-frequency standing wave amplitude distribution |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN110501063B (en) |
WO (1) | WO2021017828A1 (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110501063B (en) * | 2019-07-27 | 2021-06-04 | 复旦大学 | High-precision measurement method for high-frequency standing wave amplitude distribution |
CN113686360B (en) * | 2021-08-27 | 2023-05-26 | 中国科学院西安光学精密机械研究所 | Global measurement method and system for hemispherical harmonic oscillator standing wave drift |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102679902A (en) * | 2012-05-24 | 2012-09-19 | 天津大学 | Thin flat plate structure resonance modal analysis system and using method thereof |
CN103424175A (en) * | 2012-05-21 | 2013-12-04 | 成都星聚能科技有限公司 | Bridge amplitude real-time monitoring device based on laser modulation |
CN103438980A (en) * | 2013-09-01 | 2013-12-11 | 北京航空航天大学 | Method and device for liquid surface wave detection based on linear array CCD and linear infrared laser |
CN104776907A (en) * | 2015-04-30 | 2015-07-15 | 杭州电子科技大学 | Vibration detection method based on multi-point laser speckle extreme value tracking |
CN105043524A (en) * | 2015-04-30 | 2015-11-11 | 杭州电子科技大学 | Laser speckle defocused imaging-based vibration detection method |
JP2017003397A (en) * | 2015-06-09 | 2017-01-05 | 国立大学法人 新潟大学 | Flat surface vibration measurement device and flat surface vibration measurement method |
CN106768280A (en) * | 2017-02-28 | 2017-05-31 | 北京航空航天大学 | A kind of vibration detection device based on multi-wavelength lensless fourier transform digital hologram |
CN107300413A (en) * | 2017-07-13 | 2017-10-27 | 天津大学 | Long-range low frequency vibration measurement method based on fringe projection |
CN107560709A (en) * | 2017-08-11 | 2018-01-09 | 维沃移动通信有限公司 | A kind of vibration amplitude measurement method and measuring apparatus |
WO2018134932A1 (en) * | 2017-01-19 | 2018-07-26 | 富士通株式会社 | Amplitude measuring device and vibration measuring method in amplitude measuring device |
CN110501063A (en) * | 2019-07-27 | 2019-11-26 | 复旦大学 | A kind of high-precision measuring method of high frequency standing wave distribution of amplitudes |
Family Cites Families (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1999029139A2 (en) * | 1997-12-03 | 1999-06-10 | Sennheiser Electronic Gmbh & Co. Kg | Optical acoustic sensor, especially a microphone, with an addtional reflection device |
JP2007101394A (en) * | 2005-10-05 | 2007-04-19 | Shinetsu Film Kk | Method and apparatus for simply determining azimuth of crystal plane |
JP4976143B2 (en) * | 2006-09-29 | 2012-07-18 | 古河電気工業株式会社 | Impact vibration detector |
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 (en) * | 2012-03-21 | 2016-04-27 | スタンレー電気株式会社 | Optical system for illumination |
CN102980746B (en) * | 2012-11-20 | 2015-04-15 | 天津大学 | Crystal optical parameter measurement method based on digital holography |
CN103278935A (en) * | 2013-05-20 | 2013-09-04 | 复旦大学 | Color separation method of carrying out aberration compensation by applying free-form surface color separation film |
CN103727894A (en) * | 2014-01-15 | 2014-04-16 | 唐山英莱科技有限公司 | Transparent body three-dimensional profile detection system based on line structured light refraction imaging |
EP3289430B1 (en) * | 2015-04-27 | 2019-10-23 | Snap-Aid Patents Ltd. | Estimating and using relative head pose and camera field-of-view |
CN105115908B (en) * | 2015-08-20 | 2017-11-21 | 哈尔滨理工大学 | The optimal spectral band of metal rust spots chooses vision inspection apparatus and method |
CN105387819B (en) * | 2015-12-27 | 2017-10-13 | 河北工业大学 | Reflective object 3 D measuring method and device based on streak reflex method |
CN105571518B (en) * | 2016-01-19 | 2017-10-17 | 大连理工大学 | Three-dimensional information vision measuring method based on refraction image deviation |
CN105783775B (en) * | 2016-04-21 | 2018-05-11 | 清华大学 | A kind of minute surface and class minute surface object surface appearance measuring device and method |
US10073026B2 (en) * | 2016-07-05 | 2018-09-11 | The United States Of America, As Represented By The Secretary Of Commerce | Optical particle sorter |
CN107339941A (en) * | 2017-07-21 | 2017-11-10 | 复旦大学 | A kind of precise displacement monitoring system based on double-frequency laser interference principle |
CN108426539A (en) * | 2018-03-12 | 2018-08-21 | 四川大学 | A method of detection reflecting element face shape is projected based on double frequency orthogonal grating |
CN108413893B (en) * | 2018-03-12 | 2020-06-05 | 四川大学 | Method and device for detecting surface shape of planar element by speckle deflection technique |
CN108593088B (en) * | 2018-04-25 | 2020-09-08 | 重庆医科大学 | Focusing schlieren system for representing spherical focusing ultrasonic field |
CN108775874A (en) * | 2018-06-04 | 2018-11-09 | 西安工业大学 | A kind of shaded three-dimensional contour outline measuring set and method using random phase shift technology |
CN109186496B (en) * | 2018-10-18 | 2020-05-01 | 淮阴师范学院 | Three-dimensional surface shape measuring method based on moving least square method |
CN109870129A (en) * | 2019-03-25 | 2019-06-11 | 中国计量大学 | A kind of wafer surface roughness detection device based on phase deviation principle |
-
2019
- 2019-07-27 CN CN201910685221.XA patent/CN110501063B/en active Active
-
2020
- 2020-07-15 WO PCT/CN2020/102064 patent/WO2021017828A1/en active Application Filing
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103424175A (en) * | 2012-05-21 | 2013-12-04 | 成都星聚能科技有限公司 | Bridge amplitude real-time monitoring device based on laser modulation |
CN102679902A (en) * | 2012-05-24 | 2012-09-19 | 天津大学 | Thin flat plate structure resonance modal analysis system and using method thereof |
CN103438980A (en) * | 2013-09-01 | 2013-12-11 | 北京航空航天大学 | Method and device for liquid surface wave detection based on linear array CCD and linear infrared laser |
CN104776907A (en) * | 2015-04-30 | 2015-07-15 | 杭州电子科技大学 | Vibration detection method based on multi-point laser speckle extreme value tracking |
CN105043524A (en) * | 2015-04-30 | 2015-11-11 | 杭州电子科技大学 | Laser speckle defocused imaging-based vibration detection method |
JP2017003397A (en) * | 2015-06-09 | 2017-01-05 | 国立大学法人 新潟大学 | Flat surface vibration measurement device and flat surface vibration measurement method |
WO2018134932A1 (en) * | 2017-01-19 | 2018-07-26 | 富士通株式会社 | Amplitude measuring device and vibration measuring method in amplitude measuring device |
CN106768280A (en) * | 2017-02-28 | 2017-05-31 | 北京航空航天大学 | A kind of vibration detection device based on multi-wavelength lensless fourier transform digital hologram |
CN107300413A (en) * | 2017-07-13 | 2017-10-27 | 天津大学 | Long-range low frequency vibration measurement method based on fringe projection |
CN107560709A (en) * | 2017-08-11 | 2018-01-09 | 维沃移动通信有限公司 | A kind of vibration amplitude measurement method and measuring apparatus |
CN110501063A (en) * | 2019-07-27 | 2019-11-26 | 复旦大学 | A kind of high-precision measuring method of high frequency standing wave distribution of amplitudes |
Also Published As
Publication number | Publication date |
---|---|
CN110501063A (en) | 2019-11-26 |
CN110501063B (en) | 2021-06-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Chen et al. | Modal identification of simple structures with high-speed video using motion magnification | |
WO2021017828A1 (en) | High-precision method for measuring high-frequency standing wave amplitude distribution | |
CN113237435B (en) | High-light-reflection surface three-dimensional vision measurement system and method | |
CN107144240A (en) | A kind of system and method for detecting glass panel surface defect | |
Besnard et al. | Characterization of necking phenomena in high-speed experiments by using a single camera | |
Akhtar et al. | An efficient automotive paint defect detection system | |
Guan et al. | Defect detection method for specular surfaces based on deflectometry and deep learning | |
Ortiz et al. | Location and shape measurement using a portable fringe projection system | |
Wang et al. | Phase-based motion estimation in complex environments using the illumination-invariant log-Gabor filter | |
Quan et al. | Determination of surface contour by temporal analysis of shadow moiré fringes | |
JP2000111490A (en) | Detection apparatus for coating face | |
Wu et al. | Accurate characterisation of hole size and location by projected fringe profilometry | |
Hu et al. | Microscopic fringe projection profilometry systems in Scheimpflug condition and performance comparison | |
Zhang et al. | High-accuracy three-dimensional shape measurement of micro solder paste and printed circuits based on digital image correlation | |
Yang et al. | 3D surface profile measurement of unsymmetrical microstructure using Fizeau interferometric microscope | |
CN110702505B (en) | Double-view-field video extensometer based on telecentric lens and cubic prism | |
Zhou et al. | Robust and high-precision vision system for deflection measurement of crane girder with camera shake reduction | |
CN205383999U (en) | Object roughness optical detection system | |
TWI358525B (en) | ||
CN113899528A (en) | Liquid level dynamic wave height measuring method based on 3D-DIC principle | |
CN109781153B (en) | Physical parameter estimation method and device and electronic equipment | |
Onishi et al. | Theoretical and experimental guideline of optimum design of defect-inspection apparatus for transparent material using phase-shift illumination approach | |
Zhong et al. | Three-dimensional translation vibration measurement system based on linear array sensor and composite fringe pattern | |
Tay et al. | New method for measuring dynamic response of small components by fringe projection | |
JP2003139515A (en) | Method for measuring absolute value of deformation quantity using speckle |
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 |