CN105277274A - Vibrating rod working state determination method based on noise signal analysis - Google Patents
Vibrating rod working state determination method based on noise signal analysis Download PDFInfo
- Publication number
- CN105277274A CN105277274A CN201510616361.3A CN201510616361A CN105277274A CN 105277274 A CN105277274 A CN 105277274A CN 201510616361 A CN201510616361 A CN 201510616361A CN 105277274 A CN105277274 A CN 105277274A
- Authority
- CN
- China
- Prior art keywords
- frequency
- noise signal
- vibrating spear
- signal
- threshold values
- 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
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
The invention discloses a vibrating rod working state determination method based on noise signal analysis. The method comprises steps of 1, performing noise signal collection, wherein noise signals are timely collected via a pickup when a vibrating rod works; 2, performing time frequency transformation, wherein the collected noise signals are transformed into other frequency spectrums via FFT (Fast Fourier Transform); 3, selecting three feature frequency points in the obtained frequency spectrums via transformation, and performing energy statistics of the three feature points; and 4, obtaining the working state of the vibrating rod through a percentage threshold values and an energy statistics result. The invention can timely and effectively determine the working state of the vibrating rod, and provides an effective visualized guidance index for concrete construction.
Description
Technical field
The invention belongs to digital signal processing technique field, particularly relate to a kind of vibrating spear condition judgement method based on noise signal analysis.
Background technology
Along with interchangeable manufacturing degree is more and more higher, how to ensure that the quality of operation is the problem that must consider, in concrete construction, the reliable test of effect of vibrating is built to fresh concrete and the overall evaluation is one of important step of concrete Quality Control.Vibrating quality control in current working-yard, is control the insertion position of vibrating spear, spacing etc. of vibrating by the mode of experience, thus ensure that concrete vibrating is closely knit.In practical operation, when workman moves and to vibrate, the accurate assurance of the position/degree of depth/vibrating time cannot be accomplished, mostly often be experiential operating, by human factor and operating infulence very large; Once occur that owing shakes/cross and shake/leak shaking, vibrating can not meet construction quality and require, leave mass defect and cannot know in time.This has become concreting and to have vibrated the common fault of quality, but there is no accurately effective technical solution at present both at home and abroad.
In traditional construction quality method controls, highly rely on consciousness and the construction experience of workman, extensive form of construction work can cause various hidden dangers in project, and the later stage is once pinpoint the problems, remedial measures will be difficult to implement, and be difficult to meet rapid construction and the information-based requirement of quality control on construction.
Summary of the invention
A kind of vibrating spear duty decision method based on noise signal analysis is the invention provides for the deficiencies in the prior art, effectively can carry out real-time judgement to the duty of vibrating spear, the present invention adopts following technical scheme: a kind of vibrating spear condition judgement method based on noise signal analysis, comprises the steps:
Step 1, noise signal collection, noise signal when being worked by acoustic pickup Real-time Collection vibrating spear;
Step 2, time-frequency conversion, obtain its frequency spectrum by the noise signal collected through FFT conversion;
Step 3, select 3 characteristic frequency points converting in the signal spectrum that obtains, energy carries out to it relatively strong
Degree statistics;
Step 4, by number percent threshold values and and the comparison of energy relative intensity statistics obtain the duty of vibrating spear.
Further, the noise signal passed through when acoustic pickup Real-time Collection vibrating spear works in described step 1
Time sample frequency must meet formula (1),
F
s＞2·f
max(1)
Wherein, f
maxfor being sampled the maximal value of signal frequency, F
sfor sample frequency.
Further, in described step 2, the noise signal collected is obtained being specially in its frequency spectrum through FFT conversion: adopt the butterfly conversion in FFT conversion to carry out time-frequency conversion according to formula (2) to the noise collected and obtain its frequency spectrum,
Wherein, wherein x (n) represents the discrete digital signal sequence of input, and W is twiddle factor, and X (k) is the relative amplitude of N number of discrete point in frequency of list entries x (n) correspondence.
Further, converting selection 3 characteristic frequency points in the signal spectrum obtained described in described step 3, energy relative intensity statistics carried out to it and is specially:
Using 3 characteristic frequency points as centre frequency, calculate [middle-shift1, middle+shift1] in frequency range with [middle-shift2, middle+shift2] FFT of signal converts the ratio of the signal amplitude summation obtained in frequency range, wherein frequency centered by middle, shif1 and shift2 is different frequency offset of signal.
Further, the duty obtaining vibrating spear by number percent threshold values and energy relative intensity statistics described in described step 4 is specially: for 3 characteristic frequency points, a frequency threshold values is selected respectively from the signal spectrum obtained, the energy relative intensity statistical value of each characteristic frequency point is compared with the threshold values of its correspondence respectively, if the energy relative intensity statistical value of this characteristic frequency point is greater than corresponding threshold values, the comparative result of this characteristic frequency point is "True", and being less than is then "false"; When having two or more characteristic frequency point and the threshold values comparative result corresponding to it is "True", judging that vibrating spear has and carrying, otherwise determining that it is no-load.
Beneficial effect of the present invention: effectively can carry out real-time judgement to the duty of vibrating spear, reduce the dependence to artificial experience.
Accompanying drawing explanation
Fig. 1 is the schematic diagram of butterfly conversion;
Fig. 2 is signal frequency spectrum analysis figure of the present invention;
Fig. 3 is process flow diagram of the present invention.
Embodiment
In order to make those skilled in the art person understand technical scheme of the present invention better, and enable above-mentioned purpose of the present invention, feature and advantage become apparent more, below in conjunction with accompanying drawing, the present invention is further described.
As shown in Figure 1-Figure 3, the vibrating spear condition judgement method based on noise signal analysis of the present invention's design, design object is in order in real time, judges two kinds of duties of vibrating spear exactly: unloaded and have and carry.Based on the DSP development board of TMS320VC5402, its specific implementation process is as follows:
1, Initialize installation: close house dog, serial ports parameter, AD conversion correlation parameter (sample frequency, sampling channel etc.) are set.
2, data acquisition: by acoustic pickup, gathers the data of a period of time and (hyperchannel delays by McBSP
Rush serial port) transfer to process chip and store, wherein, sample frequency must meet formula (1),
F
s＞2·f
max(1)
Wherein, f
maxfor being sampled the maximal value of signal frequency, F
sfor sample frequency.
3, FFT (fastFouriertransform, Fourier transform) conversion.Fast Fourier Transform (FFT) is carried out to data temporary in buffer memory, obtains its frequency domain information.
The butterfly conversion be specially in adopting FFT to convert is carried out time-frequency conversion according to formula (2) to the noise collected and is obtained its frequency spectrum,
Wherein, wherein x (n) represents the discrete digital signal sequence of input, and W is twiddle factor, and X (k) is the relative amplitude of N number of discrete point in frequency of list entries x (n) correspondence.
4, energy distribution statistics: choose a characteristic frequency point from the frequency spectrum of the signal obtained, using the spectrum mean of two kinds of states of vibrating as threshold values.Therefore, use the encircled energy (referring to the amplitude of signal unit bandwidth) of signal near centre frequency as basis for estimation, namely by comparing [middle-shift1, middle+shift1] FFT of signal converts the signal amplitude summation that obtains at [middle-shift2 in frequency range, middle+shift2] between ratio, wherein frequency centered by middle, shif1 and shift2 is different frequency offset of signal.Using 3 characteristic frequency points as centre frequency, calculate [middle-shift1, middle+shift1] in frequency range with [middle-shift2, middle+shift2] FFT of signal converts the ratio of the signal amplitude summation obtained in frequency range, select 3 characteristic frequency point 200Hz, 400Hz, 600Hz, the result that current FFT converts is added up according to following rule: during middle=200Hz, shif1=50Hz, shift2=100Hz, remove statistical frequency [150 after maximal value, 250] amplitude of the signal between after FFT converts and [100, 300] ratio between, during middle=400Hz, the amplitude of signal after FFT converts between shif1=20Hz, shift2=50Hz statistical frequency [380,420] and the ratio between [350,450], during middle=600Hz, the amplitude of signal after FFT converts between shif1=20Hz, shift2=50Hz statistical frequency [580,620] and the ratio between [550,650].The actual measurement statistics of encircled energy is respectively as shown in table 1,2,3.
Statistics when table 1 centre frequency is 200Hz
Statistics when table 2 centre frequency is 400Hz
Statistics when table 3 centre frequency is 600Hz
5, state judges.According to the encircled energy number percent obtained in step 4, under can finding out three kinds of different center frequency point (200Hz, 400Hz, 600Hz), the percentage range of respective encircled energy has obvious boundary.Therefore for three kinds of states, threshold value can be chosen for 65%, 55%, 60% respectively, the result of comparative statistics with the relation of upper threshold value, being greater than is then "True", and being less than is then "false".When having two "True" in comparative result at least, be namely judged to be " have and carry "; Otherwise, be " zero load ".
By the method that this patent designs, for the differentiation of three kinds of states, the accuracy of about 90% in actual measurement, can be reached.The method had under effectively compensate for and shaking carries the higher situation of state False Rate, for visualizing monitor concrete construction quality provides the foundation effectively quantizing to judge.
More than show and describe ultimate principle of the present invention, principal character and advantage.The technician of the industry should understand, and above-described embodiment does not limit the present invention in any form, the technical scheme that the mode that all employings are equal to replacement or equivalent transformation obtains, and all drops in protection scope of the present invention.
Claims (5)
1., based on the vibrating spear duty decision method of noise signal analysis, it is characterized in that, comprise the steps:
Step 1, noise signal collection, noise signal when being worked by acoustic pickup Real-time Collection vibrating spear;
Step 2, time-frequency conversion, obtain its frequency spectrum by the noise signal collected through FFT conversion;
Step 3, select 3 characteristic frequency points converting in the signal spectrum that obtains, energy relative intensity statistics is carried out to it;
Step 4, by number percent threshold values and and the comparison of energy relative intensity statistics obtain the duty of vibrating spear.
2. the vibrating spear duty decision method based on noise signal analysis according to claim 1, is characterized in that:
During noise signal during being worked by acoustic pickup Real-time Collection vibrating spear in described step 1, sample frequency must meet formula (1),
F
s＞2·f
max(1)
Wherein, f
maxfor being sampled the maximal value of signal frequency, F
sfor sample frequency.
3. the vibrating spear duty decision method based on noise signal analysis according to claim 1, it is characterized in that, in described step 2, the noise signal collected is obtained being specially in its frequency spectrum through FFT conversion: adopt the butterfly conversion in FFT conversion to carry out time-frequency conversion according to formula (2) to the noise collected and obtain its frequency spectrum
Wherein, wherein x (n) represents the discrete digital signal sequence of input, and W is twiddle factor, and X (k) is the relative amplitude of N number of discrete point in frequency of list entries x (n) correspondence.
4. the vibrating spear duty decision method based on noise signal analysis according to claim 2, is characterized in that, is converting selection 3 characteristic frequency points in the signal spectrum obtained, carry out energy relative intensity statistics be specially it described in described step 3:
Using 3 characteristic frequency points as centre frequency, calculate [middle-shift1, middle+shift1] in frequency range with [middle-shift2, middle+shift2] FFT of signal converts the ratio of the signal amplitude summation obtained in frequency range, wherein frequency centered by middle, shif1 and shift2 is different frequency offset of signal.
5. the vibrating spear duty decision method based on noise signal analysis according to claim 1, it is characterized in that: the duty obtaining vibrating spear by number percent threshold values and energy relative intensity statistics described in described step 4 is specially: for 3 characteristic frequency points, a frequency threshold values is selected respectively from the signal spectrum obtained, the energy relative intensity statistical value of each characteristic frequency point is compared with the threshold values of its correspondence respectively, if the energy relative intensity statistical value of this characteristic frequency point is greater than corresponding threshold values, the comparative result of this characteristic frequency point is "True", be less than is then "false", when having two or more characteristic frequency point and the threshold values comparative result corresponding to it is "True", judging that vibrating spear has and carrying, otherwise determining that it is no-load.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510616361.3A CN105277274A (en) | 2015-09-24 | 2015-09-24 | Vibrating rod working state determination method based on noise signal analysis |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510616361.3A CN105277274A (en) | 2015-09-24 | 2015-09-24 | Vibrating rod working state determination method based on noise signal analysis |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN105277274A true CN105277274A (en) | 2016-01-27 |
Family
ID=55146686
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510616361.3A Pending CN105277274A (en) | 2015-09-24 | 2015-09-24 | Vibrating rod working state determination method based on noise signal analysis |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105277274A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110443527A (en) * | 2019-08-27 | 2019-11-12 | 水利部水利水电规划设计总院 | A kind of concrete vibrating quality evaluating method and device |
| CN110568073A (en) * | 2019-09-16 | 2019-12-13 | 四川升拓检测技术股份有限公司 | method for picking up impact signal in noise environment |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1146323A1 (en) * | 2000-04-14 | 2001-10-17 | Sagem Sa | Method and device for detecting the presence and energy of a given frequency in a sound and application to detection of knocking in combustion engines |
| CN101769787A (en) * | 2010-01-29 | 2010-07-07 | 浙江大学 | Method for monitoring fault of iron core looseness of power transformer |
| CN102449457A (en) * | 2009-04-02 | 2012-05-09 | 霍尼韦尔国际公司 | System and method for gearbox health monitoring |
| CN102721464A (en) * | 2012-06-13 | 2012-10-10 | 江苏省电力公司南京供电公司 | Method and system for detecting winding deformation faults of power transformer |
| CN102735968A (en) * | 2012-06-13 | 2012-10-17 | 江苏省电力公司南京供电公司 | GIS (Geographic Information System) fault diagnosis system and method based on vibration signal spectrum analysis |
| CN104727556A (en) * | 2015-04-02 | 2015-06-24 | 杨涛涛 | Concrete vibrating rod |
-
2015
- 2015-09-24 CN CN201510616361.3A patent/CN105277274A/en active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1146323A1 (en) * | 2000-04-14 | 2001-10-17 | Sagem Sa | Method and device for detecting the presence and energy of a given frequency in a sound and application to detection of knocking in combustion engines |
| CN102449457A (en) * | 2009-04-02 | 2012-05-09 | 霍尼韦尔国际公司 | System and method for gearbox health monitoring |
| CN101769787A (en) * | 2010-01-29 | 2010-07-07 | 浙江大学 | Method for monitoring fault of iron core looseness of power transformer |
| CN102721464A (en) * | 2012-06-13 | 2012-10-10 | 江苏省电力公司南京供电公司 | Method and system for detecting winding deformation faults of power transformer |
| CN102735968A (en) * | 2012-06-13 | 2012-10-17 | 江苏省电力公司南京供电公司 | GIS (Geographic Information System) fault diagnosis system and method based on vibration signal spectrum analysis |
| CN104727556A (en) * | 2015-04-02 | 2015-06-24 | 杨涛涛 | Concrete vibrating rod |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110443527A (en) * | 2019-08-27 | 2019-11-12 | 水利部水利水电规划设计总院 | A kind of concrete vibrating quality evaluating method and device |
| CN110443527B (en) * | 2019-08-27 | 2021-11-02 | 水利部水利水电规划设计总院 | Concrete vibration quality evaluation method and device |
| CN110568073A (en) * | 2019-09-16 | 2019-12-13 | 四川升拓检测技术股份有限公司 | method for picking up impact signal in noise environment |
| CN110568073B (en) * | 2019-09-16 | 2021-10-22 | 四川升拓检测技术股份有限公司 | Method for picking up impact signal in noise environment |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101907437B (en) | Wavelet difference algorithm-based cable fault localization method | |
| CN103760542B (en) | A kind of based on multimodal variability index CFAR object detection method | |
| CN106443203B (en) | A kind of pulse signal detection system and method | |
| CN102928514B (en) | Frequency characteristic-based nondestructive detection method of stress waves of wood | |
| CN108153989B (en) | Concrete dam compaction quality rapid prediction method considering parameter uncertainty influence | |
| CN104101751B (en) | Digital storage oscilloscope vertical resolution based on comentropy improves method | |
| CN106197650A (en) | A kind of method judging that audio signal is noiseless | |
| CN104795819B (en) | Power system state estimation system based on strong tracking set membership estimation | |
| CN108561126B (en) | Simple method for determining organic porosity of shale gas reservoir | |
| US20190376874A1 (en) | A method of estimating the number of modes for the sparse component analysis based modal identification | |
| CN110631849B (en) | Online identification method and system for wear state of friction system | |
| CN103995147B (en) | A kind of Data Post Processing System being applicable to acoustic Doppler velocimetry and application | |
| CN106845334A (en) | A kind of innovative noise extracting method based on mathematical morphology | |
| CN109239553A (en) | A kind of clustering method based on local density of partial discharge pulse | |
| CN105277274A (en) | Vibrating rod working state determination method based on noise signal analysis | |
| CN102103014A (en) | Detecting method for periodic transient component in signal | |
| CN105116301A (en) | Data auxiliary determining method based on dynamic statistics | |
| CN108507782B (en) | Method for detecting period signal crypto period under strong background noise | |
| CN112034387B (en) | Power transmission line short-circuit fault diagnosis method and device based on prediction sequence | |
| CN105891701A (en) | Delay performance test method and test device for fractional delay filter | |
| JP2018096866A (en) | Structure deterioration determination device and method thereof | |
| CN109580144A (en) | It is a kind of that the method for inspection is omitted based on the path OTPA for improving heavy phase dry analysis | |
| CN110333506B (en) | Method for extracting inhaul cable position parameters of cable force measurement radar | |
| CN110932755A (en) | Distributed low-voltage distribution network line parameter estimation method based on recursive least square method | |
| CN107389342B (en) | The instantaneous Frequency Estimation method examined based on Roberts operator and Grubbs |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB03 | Change of inventor or designer information | ||
| CB03 | Change of inventor or designer information |
Inventor after: Hua Zhengchao Inventor after: Li Dongxin Inventor after: Su Shuai Inventor before: Li Dongxin Inventor before: Su Shuai |
|
| TA01 | Transfer of patent application right | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20171121 Address after: 610000 Sichuan city of Chengdu province Pixian Pitong town street Chengguan Road No. 349 Applicant after: SINOHYDRO BUREAU 7 Co.,Ltd. Applicant after: HOHAI University Address before: Xikang Road, Gulou District of Nanjing city of Jiangsu Province, No. 1 211100 Applicant before: HOHAI University |
|
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20160127 |