CN101865789B - Fault detecting device of near field acoustic holography sound image mode identification and detecting method thereof - Google Patents

Abstract

The invention discloses a fault detecting device of near field acoustic holography sound image mode identification and a detecting method thereof in the field of industrial detection. The device comprises a microphone array, a scanner frame, a reference source, a clamped rib plate, an excitation source and a data acquiring system. The invention adopts sound pressure amplitude and distribution changes under mechanical normal and fault states, borrows ideas from applications of an image diagnosis technology in other fields and utilizes image processing, feature extracting and mode identifying technologies to process a near field acoustic holography image. Experiments show that the method is effective, vibration source positioning and identifying functions coordinated with the holographic imaging technology is proven, the image feather extracting and mode identifying technologies are combined with the acoustic imaging technology, thus broadening the application range of the acoustic imaging technology and bringing a set of effective acoustic fault diagnosis technologies and feasibility of being widely used in engineering thereof.

Description

Near field acoustic holography sound image pattern-recognition failure detector and detection method thereof

Technical field

What the present invention relates to is a kind of device and method of industrial detection technical field, specifically is a kind of near field acoustic holography sound image pattern-recognition failure detector and detection method thereof.

Background technology

Fault diagnosis technology based on vibration signal has limitation under some occasion, and mechanical noise is containing abundant machine state information, has the advantage of non-contact measurement, can partly substitute the means of vibration signal as fault diagnosis.The technology of utilizing noise signal to carry out fault diagnosis is referred to as the acoustics diagnose technology.Perfect acoustics diagnose technical requirement can either the effective recognition sound source position, can carry out spectrum analysis simultaneously.Traditional identification of sound source method mainly contains subjective estimate method, branch's running method, surface strength method, near field measurement method etc., but receives the restriction of physical construction and working environment easily.Acoustical holography is rapid as a kind of technological development in recent years of identification of sound source, mainly is through measuring the acoustic pressure on the two-dimensional surface (holographic facet), using restructing algorithm to come the three-dimensional sound field (comprising acoustic pressure, the sound intensity and normal direction vibration velocity) on reconstruct sound source surface.The acoustical holography method has been utilized the strength information and the phase information of sound, and visual result can be easy to noise source is positioned, quantizes, and the route of transmission of ability display noise, and acoustical holography is carried out at frequency domain simultaneously, has also inherited the characteristics of frequency-domain analysis.One Chinese patent application numbers 200910039395.5, denomination of invention is " the statistically optimal near-field acoustical holography used method and the method for operating thereof of the visual identification of air-conditioner noise sources ", discloses the application of a kind of acoustical holography in the identification of air-conditioning noise source; One Chinese patent application number 03129405.7; Denomination of invention is " method that adopts near field acoustic holography technology identification non-stationary sound source "; Further expanded the scope that near field acoustic holography is used; Through adopting the theoretical alternative traditional Fourier change technique of cyclo-stationary, select the physical quantity of spectral density function as sound field rebuilding, the near field acoustic holography technology of cyclo-stationary sound field has been proposed.

As a kind of identification of sound source algorithm, near field acoustic holography is with the obvious advantage, but the bibliographical information that directly is used for fault diagnosis seldom.Can obtain the distribution pattern of mechanical surface vibration velocity, acoustic pressure and the sound intensity behind the array acoustical signal process frequency domain near field acoustic holography, from figure, can obtain the position of vibration source under a certain frequency.

In image processing field, use very extensive based on the diagnostic techniques of image.Retrieval through to prior art is found; One Chinese patent application numbers 200710069113.7; Put down in writing a kind of " differentiating the method and apparatus of different variety green tea based on the texture analysis of multispectral image "; This technology obtains the image of different variety green tea through the 3CCD camera, calculates various green tea image texture characteristics and puts into sorter and train classification, and then unknown green tea is carried out diagnostic classification.Reference forms a kind of method for diagnosing faults based on the near field acoustic holography sound image pattern-recognition based on the diagnostic techniques of image, can expand the application of acoustical holography and promote the development that acoustics diagnose is technological.

Summary of the invention

The present invention is directed to the above-mentioned deficiency that prior art exists; A kind of near field acoustic holography sound image pattern-recognition failure detector and detection method thereof are provided; Acoustic pressure distribution or vibration velocity distribution based on machinery a certain frequency under normal and malfunction can change; Thereby in acoustic image, produce corresponding the variation; In noise source position, location, use the Feature Extraction Technology extraction sensitive features of Flame Image Process to characterize the acoustic image under normal and the failure condition through acoustic image, by the SVMs in the pattern-recognition these proper vectors are trained classification and then carried out fault diagnosis and on-line monitoring then.

The present invention realizes through following technical scheme:

The present invention relates to a kind of failure detector of near field acoustic holography sound image pattern-recognition; Comprise: scanning support, microphone array, reference source and data acquisition system (DAS); Wherein: reference source is close to machinery, is arranged in the possible noise source position of machinery, and scanning support is positioned at sound source one side; Microphone array is fixed on the scanning support according to survey sound source sound field position by different way; Can carry out side and top scanning survey to sound source, reference source microphone and scanning array microphone are connected with data acquisition system (DAS) simultaneously, storage of collected to time-domain signal and calculate testing result.

The present invention relates to the detection method of the failure detector of above-mentioned near field acoustic holography sound image pattern-recognition, may further comprise the steps:

The first step, the machinery with a plurality of noise sources is carried out normal divide and setting with the fail operation state, the situation ratio of vibration source number and position generation significant change is easier to from hologram, judge during for fault;

Second step, according to sound source face size (lx, ly) confirm holographic measurement face size (Lx, Ly), and estimate sensitive frequency F scope confirm holographic facet and source face apart from the SI between Z and the holographic measurement face microphone (Δ x, Δ y);

Described definite holographic measurement face size is meant: Lx >=1.2*lx&Ly >=1.2*ly

Be raising holographic reconstruction precision, Δ x≤λ/6, Δ y≤λ/6, Z≤λ/6 and Z>=max (Δ x, Δ y), wherein λ is the wavelength of being confirmed by F.

The 3rd step, the sound pressure signal P (t) in measurement mechanical under the consistent condition of ground unrest under the normal and malfunction; And establishment sample storehouse S; Carry out spectrum analysis and find sensitive frequency f (f≤F), each normal sound field sample, fault sound field sample are adopted near field acoustic holography.

The sound pressure signal of any point p can be expressed as the form of Helmholtz integral equation in the described sound field:

$p\left({\stackrel{\→}{r}}_p\right)=\underset{S}{\∫\∫}[p\left({\stackrel{\→}{r}}_Q\right)\frac{\∂G({\stackrel{\→}{r}}_p,{\stackrel{\→}{r}}_Q)}{\∂n}+\mathrm{j\ρck}{u}_n\left({\stackrel{\→}{r}}_Q\right)G({\stackrel{\→}{r}}_p,{\stackrel{\→}{r}}_Q)]\mathrm{dS}\left({\stackrel{\→}{r}}_Q\right)---\left(1\right)$

Wherein: S representes the sound source surface, and Q representes the border, U _nRepresentation is to vibration velocity, and G representes Green function.According to the suitable holographic algorithm of sound source structure choice: the acoustic pressure on the acoustical holography utilization measurement face of quadrature conformal structure is the convolution of source face acoustic pressure and Green function; Fast two-dimensional fourier transformation is used for the Helmholtz equation, has realized rebuilding the distribution of acoustic pressure, vibration velocity and the sound intensity on the face of source by the acoustic pressure of holographic measurement face; The near field acoustic holography of arbitrary shape structure has BEM-based NAH and the equivalent source method based on numerical solution Helmholtz integral equation.

Calculate frequency sound source surface vibration velocity distribution matrix A when being F through near field acoustic holography, A is a complex matrix.Choose appropriate gray shade level n, matrix A is visualized as image B, and image B is the vibration velocity distribution pattern, position that therefrom can the main vibration source of identification.

The 4th step, the near field acoustic holography result is carried out feature extraction, structural attitude vector sample storehouse is convenient to sorter and is trained identification;

Described feature extraction is meant: adopt in the following dual mode any one:

A) directly the reconstruct matrix A is carried out svd and obtains singular value features:

$A=U{\mathrm{\ΛV}}^T=\underset{i=1}{\overset{k}{\mathrm{\Σ}}}{\mathrm{\α}}_i{u}_i{v}_i^T=\underset{i=1}{\overset{k}{\mathrm{\Σ}}}{\mathrm{\α}}_i{A}_i---\left(2\right)$

Wherein: k is the order of matrix A, α _iBe i singular value of matrix A, make α ₁>=a ₂>=...>=α _n, m singular value constitutive characteristic vector is used for the training classification before getting.

B) be to utilize image processing techniques to extract the texture statistics characteristic of image B:

The known image gray level is n, the second order joint probability density function P through asking for image (d θ) obtains gray level co-occurrence matrixes for i, j:

M(d，θ)＝[P(i，j，d，θ)] (3)

Wherein: function P has described image distance on the θ direction has the appearance of gray scale i and j respectively for a pair of pixel of d probability; D=1; θ gets 0 °, 45 °, 90 ° and 135 ° respectively; Obtain the gray level co-occurrence matrixes of four direction; And ask for 12 kinds of textural characteristics coefficients based on gray level co-occurrence matrixes: angle second moment value, contrast value, correlation, entropy, variance yields, contrary gap value and mean value and variance yields and entropy, difference mean value, difference variance and difference entropy, and be used for pattern-recognition.

The 5th the step, obtain proper vector sample storehouse after; Put into SVMs and train classification; Contrast the discrimination under the various conditions, obtain best classifier parameters: punishment parameters C and kernel function coefficient gamma, and feature extraction mode; The array signal that collects is in real time carried out feature extraction and adopts the sorter after optimizing to discern, realize fault diagnosis.

The present invention utilize machinery normal with malfunction under the variation of sound pressure amplitude and distribution, reference diagnostic imaging technology is in the application of other field, employing Flame Image Process, feature extraction and mode identification technology are handled the near field acoustic holography image.Experiment shows the validity of method; Simultaneous verification cooperate vibration source location, the recognition function of holographic imaging technology; Combine image treatment features extraction and mode identification technology with the acoustic imaging technology; Widen acoustic imaging The Application of Technology scope, form the feasibility of effective Acoustic Based Diagnosis technology of a cover and widespread use in engineering.

Description of drawings

Fig. 1 is a structural representation of the present invention.

Fig. 2 is exciting point position and floor synoptic diagram.

Fig. 3 is a schematic flow sheet of the present invention.

Fig. 4 is that the embodiment test sound is as comparison diagram.

Embodiment

Elaborate in the face of embodiments of the invention down, present embodiment provided detailed embodiment and concrete operating process, but protection scope of the present invention is not limited to following embodiment being to implement under the prerequisite with technical scheme of the present invention.

As shown in Figure 1, present embodiment comprises: scanning support 1, microphone array 2, prop up floor 3, reference source 4, exciting source 5, data acquisition system (DAS) 6 and computer 7 admittedly, wherein: microphone array is arranged in the side of floor sound source, apart from sound source surface Z; Because floor exciting source excited frequency is identical, three vibration sources are relevant, adopt a reference source; Position through scanning support control microphone array forms the holographic measurement face; Gather the sound field information of floor sound source through data acquisition system (DAS) in each scanning step.

As shown in Figure 2, described exciting point position and floor: embodiment carries out 3 excitings to floor, and each exciting point distance is greater than the resolution of near-field holography; Floor props up on experimental stand admittedly.

As shown in Figure 3, said apparatus carries out work in the following manner:

Experiment is carried out in semianechoic room, adopts multi-point exciting to prop up the overriding noise source that floor is simulated machinery admittedly, and is as shown in Figure 2.The thickness of floor and rib all is 4mm; Floor area 640mm * 440mm; Rib 1 area 360mm * 40mm; Rib 2 area 560mm * 60m.The floor center is a true origin, exciting point 1 (180mm, 87mm), exciting point 2 (20mm, 87mm), exciting point 3 (120mm ,-87mm).

The first step, the machinery with a plurality of noise sources is carried out normal dividing and setting with the fail operation state:

Connect through force transducer between vibrator and the floor; Simulate the normal and malfunction of each vibration source through measuring and adjust the exciting force size; As: vibration source 1 exciting force is in normal operating condition when between 6N～8.8N, changing, be malfunction greater than the situation of 8.8N; Vibration source 2 exciting forces are in normal operating condition when between 4N～5.8N, changing, be malfunction greater than the situation of 5.8N; Vibration source 3 exciting forces are in normal operating condition when between 5N～7.4N, changing, be malfunction greater than the situation of 7.4N.Set 30 in normal condition sample altogether, 35 of various failure operation state sample.

Second the step, estimate sensitive frequency scope F, be lower than the first natural frequency 473Hz that props up floor admittedly.By sensitive frequency F confirm holographic facet and source face apart from the SI between Z and the holographic measurement face microphone (Δ x, Δ y), according to sound source face size confirm holographic measurement face size (Lx, Ly);

As shown in Figure 1: because sound source is regular, embodiment adopts the FFT-based NAH based on the quadrature conformal structure; Owing to the restriction of acquisition instrument equipment, the method that embodiment adopts linear array scanning to measure is gathered acoustical signal simultaneously; Arranging the microphone linear array apart from sound source face Z=0.1m place, linear array is made up of 9 microphones, and adjacent microphone is at a distance of Δ x=Δ y=0.1m; The sound source face is 0.64m * 0.44mm, scanning holography face 0.9m * 0.8m; Sound field is concerned with, and near vibration source 2, places microphone source as a reference; The acoustic field signal of each 10 passages of synchronous acquisition, and the volume coordinate position of record microphone.

The 3rd step, measurement mechanical under the consistent condition of ground unrest normal with malfunction under sound pressure signal P (t), and set up sample storehouse S, set sensitive frequency f=300Hz, to each normal sound field sample, fault sound field sample employing near field acoustic holography.

The acoustic pressure of any point p can be expressed as the form of Helmholtz integral equation in the sound field:

$p\left({\stackrel{\→}{r}}_p\right)=\underset{S}{\∫\∫}[p\left({\stackrel{\→}{r}}_Q\right)\frac{\∂G({\stackrel{\→}{r}}_p,{\stackrel{\→}{r}}_Q)}{\∂n}+\mathrm{j\ρck}{u}_n\left({\stackrel{\→}{r}}_Q\right)G({\stackrel{\→}{r}}_p,{\stackrel{\→}{r}}_Q)]\mathrm{dS}\left({\stackrel{\→}{r}}_Q\right)---\left(1\right)$

Wherein: S representes the sound source surface, and Q representes the border, U _nRepresentation is to vibration velocity, and G representes Green function.According to the suitable holographic algorithm of sound source structure choice: the acoustic pressure on the acoustical holography utilization measurement face of quadrature conformal structure is the convolution of source face acoustic pressure and Green function; Fast two-dimensional fourier transformation is used for the Helmholtz equation, has realized rebuilding the distribution of acoustic pressure, vibration velocity and the sound intensity on the face of source by the acoustic pressure of holographic measurement face; The near field acoustic holography of arbitrary shape structure has BEM-based NAH and the equivalent source method based on numerical solution Helmholtz integral equation.

Calculate frequency sound source surface vibration velocity distribution matrix A when being f through near field acoustic holography, A is a complex matrix.Choose appropriate gray shade level n, matrix A is visualized as image B, and image B is the vibration velocity distribution pattern, position that therefrom can the main vibration source of identification.

65 sound field samples are carried out FFT-based NAH calculate source image restructuring matrix and acoustic image, Fig. 4 a is the acoustic image under normal and the malfunction.Wherein can be referring to patent 03129405.7 based on the near field acoustic holography theory of algorithm of scanning method, patent name is " method that adopts near field acoustic holography technology identification non-stationary sound source ", steadily the cycle frequency of sound field is 0.The vibration velocity restructuring matrix is the complex matrix that comprises phase information, and gray level n get respectively 8,16,32,64 and 128 o'clock the visual acoustic image that changes into of matrix.Observe the vibration source position, with actual exciting position versus, because floor single order Mode Shape and influence of measurement error, vibration source and actual exciting position do not overlap, but can accurately locate the general location of main vibration source from the vibration velocity distributed image.

The 4th goes on foot, restructuring matrix is directly extracted singular value features, calculates the textural characteristics based on gray level co-occurrence matrixes to confirming the acoustic image after the gray level, forms proper vector sample storehouse.

The 5th the step, obtain proper vector sample storehouse after; Put into SVMs and train classification; Contrast the discrimination under the various conditions, obtain best classifier parameters: punishment parameters C and kernel function coefficient gamma, and feature extraction mode; The array signal that collects is in real time carried out feature extraction and adopts the sorter after optimizing to discern, realize fault diagnosis.

Use Libsvm to sample training classification, wherein as shown in table 1 based on the classification results of textural characteristics, as shown in table 2 based on the classification results of singular value features:

Table 1 is based on discrimination contrast (%) under each gray level of statistical nature

Can see that by table 1 discrimination is the highest when gray level is 32,45 ° of directions, reach 90.7%, punish parameters C=2048, kernel function coefficient gamma=0.0078125 this moment.

Table 2 is based on the discrimination contrast (%) of singular value

By table 2, the singular value number is that 8 o'clock recognition effects preferably reach 92.3%, punishment parameters C=32768, kernel function coefficient gamma=0.0078125.

Under the present embodiment condition, it is best to choose the recognition effect that preceding 8 big singular value constitutive characteristic vectors train classification to obtain, and has reached 92.3%, the limitations affect of sample size recognition effect, but surpass the validity that 90% discrimination has been explained method.

Claims (1)

1. the detection method of the failure detector of a near field acoustic holography sound image pattern-recognition, this device comprises: scanning support, microphone array, reference source and data acquisition system (DAS), wherein: reference source is close to machinery; Be arranged in the possible noise source position of machinery; Scanning support is positioned at noise source one side, and microphone array is fixed on the scanning support according to survey noise source sound field position by different way, can carry out side and top scanning survey to noise source; The reference source microphone is connected with data acquisition system (DAS) with the scanning array microphone simultaneously; Storage of collected to time-domain signal and calculate testing result, it is characterized in that, said method comprising the steps of:

The first step, the machinery with a plurality of noise sources is carried out normal divide and setting with the fail operation state, the situation ratio of vibration source number and position generation significant change is easier to from hologram, judge during for fault;

Second step, estimate the scope of sensitive frequency F, confirm holographic facet and source face apart from SI Δ x and Δ y between Z and the holographic measurement face microphone, according to big or small definite holographic measurement face size Lx of noise source face and Ly;

The 3rd step, measurement mechanical under the consistent condition of ground unrest normal with malfunction under sound pressure signal P (t), and set up sample storehouse S, carry out spectrum analysis and find sensitive frequency f, to each normal sound field sample, fault sound field sample employing near field acoustic holography;

The 4th step, the near field acoustic holography result is carried out feature extraction, structural attitude vector sample storehouse is convenient to sorter and is trained identification;

The 5th the step, obtain proper vector sample storehouse after; Put into SVMs and train classification; Contrast the discrimination under the various conditions, obtain best classifier parameters: punishment parameters C and kernel function coefficient gamma, and feature extraction mode; The array signal that collects is in real time carried out feature extraction and adopts the sorter after optimizing to discern, realize fault diagnosis;

Described definite holographic measurement face size is meant: Lx>=1.2*lx&Ly>=1.2*ly, wherein: Δ x≤λ/6, Δ y≤λ/6; Z≤λ/6 and Z>=max (Δ x; Δ y), λ is the wavelength of being confirmed by F, and wherein: lx and ly are noise source face size; Δ x and Δ y are the SI between the holographic measurement face microphone, and F is for estimating sensitive frequency;

Described sound pressure signal is:

$p\left({\stackrel{\→}{r}}_p\right)=\∫\underset{S}{\∫}[p\left({\stackrel{\→}{r}}_Q\right)\frac{\∂G({\stackrel{\→}{r}}_p,{\stackrel{\→}{r}}_Q)}{\∂n}+\mathrm{j\ρck}{u}_n\left({\stackrel{\→}{r}}_Q\right)G({\stackrel{\→}{r}}_p,{\stackrel{\→}{r}}_Q)]\mathrm{dS}\left({\stackrel{\→}{r}}_Q\right),$

Wherein: S representes the noise source surface, and Q representes the border, u _nRepresentation is to vibration velocity; G representes Green function, and n is the gradation of image number, and c is the punishment parameter; K is the order of reconstruct matrix A; According to the suitable holographic algorithm of noise source structure choice: the acoustic pressure on the acoustical holography utilization measurement face of quadrature conformal structure is used for the Helmholtz equation for the convolution of source face acoustic pressure and Green function with fast two-dimensional fourier transformation, has realized being rebuild by the acoustic pressure of holographic measurement face the distribution of acoustic pressure, vibration velocity and the sound intensity on the face of source; The near field acoustic holography of arbitrary shape structure has BEM-based NAH and the equivalent source method based on numerical solution Helmholtz integral equation.

Images