CN106355008A - Method for calculating contact point during force-displacement measurement - Google Patents
Method for calculating contact point during force-displacement measurement Download PDFInfo
- Publication number
- CN106355008A CN106355008A CN201610767641.9A CN201610767641A CN106355008A CN 106355008 A CN106355008 A CN 106355008A CN 201610767641 A CN201610767641 A CN 201610767641A CN 106355008 A CN106355008 A CN 106355008A
- Authority
- CN
- China
- Prior art keywords
- contact point
- measurement data
- measurement
- computational methods
- displacement measurement
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Length Measuring Devices With Unspecified Measuring Means (AREA)
Abstract
The invention discloses a method for calculating a contact point during force-displacement measurement. The method comprises the following steps: transforming data by fast Fourier transform; after carrying out spectral analysis, determining a filtering algorithm and filtering the data; and processing the data by using an iterative algorithm to determine the final contact point. Calculating precision is high, calculating time is short, calculating errors are small, and errors of temperature excursion and system errors of nonzero value displacement calculation can be eliminated effectively.
Description
Technical field:
The invention belongs to product field of measuring technique, it is specifically related to a kind of calculating side of contact point during power-displacement measurement
Method.
Background technology:
For the measurement of elastic deformation object power-displacement, in gatherer process, measuring system will necessarily produce random noise,
And due to impacts such as system is vibrated, temperature, data has unusual fluctuations.
For above-mentioned impact, traditional method, when judging contact point, typically can go to contact as calculating with the power more than 0 value
Position, but this can bring both sides error, and one is the error of temperature drift, and two is the systematic error that nonzero value calculates displacement.That is such as
Where there is the accurate corresponding relation accurately finding out contact point power-displacement under above-mentioned disturbed condition, for the meaning of accurate measurement
Very big.
Content of the invention:
For this reason, the technical problem to be solved is in prior art to be used for elastic deformation object power-displacement
Measurement contact point determination methods using more than 0 power area as calculate contact position, thus bringing error and the non-zero of temperature drift
Value calculates the systematic error of displacement, thus proposing a kind of computational methods of contact point during power-displacement measurement.
For reaching above-mentioned purpose, technical scheme is as follows:
During a kind of power-displacement measurement, the computational methods of contact point, comprise the steps:
S1: obtain the first measurement data, described first measurement data is analyzed and filters, obtain second after denoising
Measurement data.
S2: determine drift.
S3: using iterative algorithm, described second measurement data is calculated, and determine step-size in search.
S4: judging whether the second-order differential value of the second measurement data exceedes first threshold, if exceeding, entering step s5,
Otherwise return to step s3.
S5: judging whether the second measurement data exceedes drift, if exceeding, entering step s6, otherwise return to step s3.
S6: judging whether the first differential value of the second measurement data exceedes Second Threshold, if exceeding, entering step s7,
Otherwise return to step s3；
S7: terminate, determine final contact point.
Preferred as technique scheme, described step s1 specifically includes following steps:
S11: obtain the first measurement data.
S12: using becoming scaling method, the first measurement data is entered with line translation, obtain spectrogram.
S13: the spectrogram in step s12 is analyzed.
S14: determine filtering algorithm.
S15: remove the noise in the first measurement data using described filtering algorithm, obtain the second measurement number after denoising
According to.
Preferred as technique scheme, in described step s12:
Described change scaling method adopts fast fourier transform algorithm.
Preferred as technique scheme, in described step s2:
Described drift choose not in contact with when power average.
Preferred as technique scheme, in described step s4:
Described first threshold is passed through manually to be configured.
Preferred as technique scheme, in described step s6:
Described Second Threshold passes through manually to be configured.
The beneficial effects of the present invention is: it is converted to data by fast Fourier transform, through spectrum analyses
Determine that filtering algorithm is filtered to data processing afterwards, and data is processed to determine final contact using iterative algorithm
Point, computational accuracy height, calculating time period, calculation error are little, can effectively eliminate the error of temperature drift and nonzero value calculates displacement
Systematic error.
Brief description:
The following drawings is only intended to, in doing schematic illustration and explanation to the present invention, not delimit the scope of the invention.Wherein:
Fig. 1 is the computational methods flow chart of contact point during a kind of power-displacement measurement of one embodiment of the invention；
Fig. 2 is that the filtering algorithm of one embodiment of the invention determines flow chart；
Fig. 3 is the contact point calculation flow chart of one embodiment of the invention.
Specific embodiment:
As shown in figure 1, during a kind of power-displacement measurement of the present invention contact point computational methods, comprise the steps:
S1: obtain the first measurement data, described first measurement data is analyzed and filters, obtain second after denoising
Measurement data.As shown in Fig. 2 described step s1 specifically includes following steps:
S11: obtain the first measurement data.
S12: using becoming scaling method, the first measurement data is entered with line translation, obtain spectrogram.In the present embodiment, described conversion
Algorithm adopts fast fourier transform algorithm.
S13: the spectrogram in step s12 is analyzed.
S14: determine filtering algorithm.
S15: remove the noise in the first measurement data using described filtering algorithm, obtain the second measurement number after denoising
According to.
S2: determine drift.In the present embodiment, described drift choose not in contact with when power average.
S3: using iterative algorithm, described second measurement data is calculated, and determine step-size in search.
S4: judging whether the second-order differential value of the second measurement data exceedes first threshold, if exceeding, entering step s5,
Otherwise return to step s3.In the present embodiment, described first threshold is passed through manually to be configured.
S5: judging whether the second measurement data exceedes drift, if exceeding, entering step s6, otherwise return to step s3.
S6: judging whether the first differential value of the second measurement data exceedes Second Threshold, if exceeding, entering step s7,
Otherwise return to step s3.In the present embodiment, described Second Threshold passes through manually to be configured.
S7: terminate, determine final contact point.
The operation principle of the present invention is that the concrete calculation process of contact point is as shown in Figure 3:
1st, use fast fourier transform fft analysis spectrum figure, remove the noise in data with wave filter.
2nd, be used not in contact with when power average as drift.
3rd, use iterative algorithm, determine step-size in search, by the threshold decision of second-order differential, search interval interior value more than zero
The counting of drift and first differential value exceed the counting of threshold value to determine final contact point.
The computational methods of contact point during a kind of power-displacement measurement described in the present embodiment, it passes through fast Fourier transform
Data is converted, determines that filtering algorithm is filtered to data processing after spectrum analyses, and use iterative algorithm pair
Data is processed to determine final contact point, and computational accuracy height, calculating time period, calculation error are little, can effectively disappear
Except the error of temperature drift and the systematic error of nonzero value calculating displacement.
Obviously, above-described embodiment is only intended to clearly illustrate example, and the not restriction to embodiment.Right
For those of ordinary skill in the art, can also make on the basis of the above description other multi-forms change or
Change.There is no need to be exhaustive to all of embodiment.And the obvious change thus extended out or
Change among still in the protection domain of the invention.
Claims (6)
1. during a kind of power-displacement measurement the computational methods of contact point it is characterised in that comprising the steps:
S1: obtain the first measurement data, described first measurement data is analyzed and filters, obtain the second measurement after denoising
Data；
S2: determine drift；
S3: using iterative algorithm, described second measurement data is calculated, and determine step-size in search；
S4: judging whether the second-order differential value of the second measurement data exceedes first threshold, if exceeding, entering step s5, otherwise
Return to step s3；
S5: judging whether the second measurement data exceedes drift, if exceeding, entering step s6, otherwise return to step s3；
S6: judging whether the first differential value of the second measurement data exceedes Second Threshold, if exceeding, entering step s7, otherwise
Return to step s3；
S7: terminate, determine final contact point.
2. during power-displacement measurement according to claim 1 the computational methods of contact point it is characterised in that described step s1
Specifically include following steps:
S11: obtain the first measurement data；
S12: using becoming scaling method, the first measurement data is entered with line translation, obtain spectrogram；
S13: the spectrogram in step s12 is analyzed；
S14: determine filtering algorithm；
S15: remove the noise in the first measurement data using described filtering algorithm, obtain the second measurement data after denoising.
3. during power-displacement measurement according to claim 2 the computational methods of contact point it is characterised in that described step s12
In:
Described change scaling method adopts fast fourier transform algorithm.
4. during power-displacement measurement according to claim 1 and 2 the computational methods of contact point it is characterised in that described step
In s2:
Described drift choose not in contact with when power average.
5. during power-displacement measurement according to claim 1 the computational methods of contact point it is characterised in that described step s4
In:
Described first threshold is passed through manually to be configured.
6. during power-displacement measurement according to claim 1 the computational methods of contact point it is characterised in that described step s6
In:
Described Second Threshold passes through manually to be configured.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610767641.9A CN106355008A (en) | 2016-08-30 | 2016-08-30 | Method for calculating contact point during force-displacement measurement |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610767641.9A CN106355008A (en) | 2016-08-30 | 2016-08-30 | Method for calculating contact point during force-displacement measurement |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106355008A true CN106355008A (en) | 2017-01-25 |
Family
ID=57856650
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610767641.9A Pending CN106355008A (en) | 2016-08-30 | 2016-08-30 | Method for calculating contact point during force-displacement measurement |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106355008A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112509013A (en) * | 2020-12-04 | 2021-03-16 | 杭州鲁尔物联科技有限公司 | Method for detecting precision of image displacement detection algorithm |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120108988A1 (en) * | 2010-11-03 | 2012-05-03 | Doron Moshe Ludwin | Zero-drift detection and correction in contact force measurements |
| CN103364116A (en) * | 2013-07-25 | 2013-10-23 | 安徽皖科电子工程有限公司 | Pressure measurement method of high-precision pressure sensor for time difference oscillation |
| CN103481121A (en) * | 2013-09-16 | 2014-01-01 | 同济大学 | Two-tooth helical end mill tooth radius asymmetry error detection system and method |
| CN103941756A (en) * | 2014-04-22 | 2014-07-23 | 中国人民解放军国防科学技术大学 | Non-contact active control method for vibration of micro thrust measurement system |
| CN105486454A (en) * | 2015-12-21 | 2016-04-13 | 北京航天时代激光导航技术有限责任公司 | Laser inertial measurement unit product mass center measurement method |
-
2016
- 2016-08-30 CN CN201610767641.9A patent/CN106355008A/en active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120108988A1 (en) * | 2010-11-03 | 2012-05-03 | Doron Moshe Ludwin | Zero-drift detection and correction in contact force measurements |
| CN103364116A (en) * | 2013-07-25 | 2013-10-23 | 安徽皖科电子工程有限公司 | Pressure measurement method of high-precision pressure sensor for time difference oscillation |
| CN103481121A (en) * | 2013-09-16 | 2014-01-01 | 同济大学 | Two-tooth helical end mill tooth radius asymmetry error detection system and method |
| CN103941756A (en) * | 2014-04-22 | 2014-07-23 | 中国人民解放军国防科学技术大学 | Non-contact active control method for vibration of micro thrust measurement system |
| CN105486454A (en) * | 2015-12-21 | 2016-04-13 | 北京航天时代激光导航技术有限责任公司 | Laser inertial measurement unit product mass center measurement method |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112509013A (en) * | 2020-12-04 | 2021-03-16 | 杭州鲁尔物联科技有限公司 | Method for detecting precision of image displacement detection algorithm |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Eichstädt et al. | Deconvolution filters for the analysis of dynamic measurement processes: a tutorial | |
| CN106096313B (en) | A kind of envelope Analysis Method based on unusual spectral factorization and spectrum kurtosis | |
| CN106225786B (en) | A kind of adaptive pedestrian navigation system zero-speed section detecting method | |
| RU2016116927A (en) | DEVICE, METHOD AND SYSTEM FOR PROCESSING PHYSIOLOGICAL SIGNAL | |
| ATE511128T1 (en) | STATISTICAL PROCESSING METHODS FOR DETECTING ANORMAL SITUATIONS | |
| CN105298889B (en) | Gas compressor surge detection method | |
| CN103185837B (en) | Method for measuring frequency of power system | |
| CN104880172B (en) | Measuring flatness of road surface method and device based on Kalman filtering | |
| CN103728663B (en) | A kind of Time-Frequency Analysis Method | |
| CN103913765B (en) | A kind of nucleic power spectrum Peak Search Method | |
| CN104679375A (en) | Method and device for optimizing signal-to-noise ratio parameter | |
| CN103823177B (en) | Performance of filter detection method based on THE DESIGN OF WINDOW FUNCTION and system | |
| CN102937593A (en) | Ceramic radome crack automatic detection method | |
| CN106096198B (en) | A kind of envelope Analysis Method based on variation Mode Decomposition and spectrum kurtosis | |
| CN106053069B (en) | A kind of SSD of rolling bearing, spectrum kurtosis and smooth iteration envelope Analysis Method | |
| CN115728259A (en) | Multi-component gas analysis method | |
| CN113901379A (en) | Dynamic online rapid processing method for real-time data of edge terminal | |
| CN106355008A (en) | Method for calculating contact point during force-displacement measurement | |
| CN106096199B (en) | A kind of WT of rolling bearing, spectrum kurtosis and smooth iteration envelope Analysis Method | |
| CN107656134B (en) | fault current detection method and device for adaptively filtering direct-current component | |
| CN105954031B (en) | A kind of envelope Analysis Method based on unusual spectral factorization filtering | |
| KR100722832B1 (en) | Method for adaptive notch filtering based on robust frequency estimation | |
| CN103580647B (en) | A kind of filter construction | |
| CN103490830A (en) | Denoising radio frequency spectrum peak value acquisition method based on Internet of things electric power temperature measurement equipment | |
| CN104931785A (en) | PCB characteristic impedance test method based on mode approach |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB02 | Change of applicant information | ||
| CB02 | Change of applicant information |
Address after: The Wujiang economic and Technological Development Zone West Road Wujiang District of Suzhou City, Jiangsu Province, No. 666 215200 Applicant after: Bo Seiko Polytron Technologies Inc Address before: The Wujiang economic and Technological Development Zone West Road Wujiang District of Suzhou City, Jiangsu Province, No. 666 215200 Applicant before: Suzhou Bozhong Precision Industry Technology Co., Ltd. |
|
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170125 |