CN110595811B - Method for constructing health state characteristic diagram of mechanical equipment - Google Patents

Method for constructing health state characteristic diagram of mechanical equipment Download PDF

Info

Publication number
CN110595811B
CN110595811B CN201910860506.2A CN201910860506A CN110595811B CN 110595811 B CN110595811 B CN 110595811B CN 201910860506 A CN201910860506 A CN 201910860506A CN 110595811 B CN110595811 B CN 110595811B
Authority
CN
China
Prior art keywords
frequency
time
spectrum
fmap
characteristic diagram
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.)
Active
Application number
CN201910860506.2A
Other languages
Chinese (zh)
Other versions
CN110595811A (en
Inventor
郭远晶
杨友东
钱少明
林勇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhijiang College of ZJUT
Original Assignee
Zhijiang College of ZJUT
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Zhijiang College of ZJUT filed Critical Zhijiang College of ZJUT
Priority to CN201910860506.2A priority Critical patent/CN110595811B/en
Publication of CN110595811A publication Critical patent/CN110595811A/en
Application granted granted Critical
Publication of CN110595811B publication Critical patent/CN110595811B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01HMEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H17/00Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M99/00Subject matter not provided for in other groups of this subclass
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/044Recurrent networks, e.g. Hopfield networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/045Combinations of networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods

Abstract

The method for constructing the health state characteristic diagram of the mechanical equipment comprises the steps of firstly carrying out time-frequency transformation on vibration signals collected on the mechanical equipment, carrying out compression, normalization and graying on the obtained time-frequency spectrum, obtaining a time-frequency characteristic gray value, and filling the time-frequency characteristic gray value into the health state characteristic diagram; then, Fourier transform is carried out on the vibration signal, frequency band equal division and graying processing are carried out on the frequency spectrum, a frequency domain characteristic gray value is obtained, and the frequency domain characteristic gray value is filled into a health state characteristic diagram; and finally, calculating a time domain dimensionless index of the vibration signal, acquiring a time domain characteristic gray value, and filling the time domain characteristic gray value into the health state characteristic diagram. The health state characteristic diagram obtained by the invention is small in size, can fully reflect the health state information of mechanical equipment, and improves the fault diagnosis accuracy.

Description

Method for constructing health state characteristic diagram of mechanical equipment
Technical Field
The invention relates to a method for constructing a health state characteristic diagram of mechanical equipment.
Background
With the continuous development of technologies such as computers, sensors, communication, artificial intelligence and the like, intelligent health state monitoring and fault diagnosis technologies of mechanical equipment are also rapidly developed. The intelligent fault diagnosis technology overcomes the contradiction between the massive health state monitoring data of the mechanical equipment operation and the relative scarcity of fault diagnosis professionals, can automatically process and analyze the health state data of the mechanical equipment, and can give out an accurate fault diagnosis result. In recent years, deep learning has been rapidly developed in academic and industrial fields by virtue of its strong pattern recognition capability. The fault diagnosis of mechanical equipment is also essentially pattern recognition, so that a plurality of scholars and experts at home and abroad research and utilize a deep learning technology to realize intelligent fault diagnosis of the equipment.
The basic network models commonly used for deep learning include a Deep Belief Network (DBN), a Convolutional Neural Network (CNN), a stacked automatic coding machine (SAE), a Recurrent Neural Network (RNN), a generative countermeasure network (GAN), and the like. Among them, CNN is very suitable for processing and learning of mass data as a supervised deep learning model, and has been successfully implemented in commercial applications at present. Since the input of the CNN is a two-dimensional graph, and the two-dimensional graph is usually an N × N square graph for the sake of processing simplicity, one of the key points in applying the CNN to intelligent fault diagnosis of mechanical equipment is how to obtain a two-dimensional graph that sufficiently reflects health status information from equipment health status data (usually vibration data).
Wuchuanzhi, Qujianling, Populi Hongbai and the like adopt one-dimensional time domain vibration signals as input of CNN, and are respectively used for fault diagnosis of a gear box, a rolling bearing and a reciprocating compressor. Luyang king et al utilize vibration signal frequency domain data as CNN input for gearbox health monitoring. The Liheng and Wanglihua utilize short-time Fourier transform time frequency spectrum as CNN input, and are respectively used for fault diagnosis of bearings and motors. Chenrenxiang and the like utilize discrete wavelet transformation to construct a two-dimensional time-frequency matrix which is used as CNN input for fault diagnosis of the rolling bearing. And the wavelet time-frequency graph is input as CNN by Yuan-Jian and the like and is used for fault diagnosis of the rolling bearing. And constructing a time-frequency diagram by using a Fourier transform method for the Liu (registered public network) set and the like as the input of the CNN, wherein the time-frequency diagram is used for diagnosing the fault of the rolling bearing.
At present, the input of a fault diagnosis CNN model of mechanical equipment usually only considers one of three analysis domains, namely a vibration signal time domain, a vibration signal frequency domain or a vibration signal time-frequency domain, and characteristic information in each analysis domain cannot be comprehensively considered, so that the health state of the mechanical equipment is difficult to be fully reflected, and the accuracy of the CNN fault diagnosis model for judging the health state of the mechanical equipment is influenced.
Disclosure of Invention
The invention provides a method for constructing a health state characteristic diagram of mechanical equipment, aiming at overcoming the defects of a method for acquiring the health state characteristic diagram of the mechanical equipment in the prior art.
The invention discloses a method for constructing a health state characteristic diagram of mechanical equipment, which comprises the following steps of:
step one, presetting the dimension of a health state characteristic diagram FMap as NXN;
step two, collecting a vibration signal x (L) with the length of L on mechanical equipment, wherein L is 1,2, …, L;
step three, acquiring a time-frequency characteristic gray value of the vibration signal x (l) and filling the time-frequency characteristic gray value into a health state characteristic diagram FMap, wherein the method specifically comprises the following steps:
(31) performing time-frequency transform processing on x (L) to obtain a time spectrum TFP (t, f), wherein t is 1,2, …, L, f is 1,2, …, L/2;
(32) the time axis N of the time spectrum TFP (t, f) is equally divided, the time point number of each equal division is delta t equal to L/N, the frequency point number of each equal division is delta f for the frequency axis N/21For (L/2)/(N/2) L/N, so that for time spectrum TFP (t, f) the block is divided, defining the position (t)d,fd) The time-frequency block of (d) is TFP ((t)d-1)·△t+i,(fd-1)·△f1+j1) Wherein t isd=1,2,…,N,fd=1,2,…,N/2,i=1,2,…,△t,j1=1,2,…,△f1
(33) The element values in each time-frequency block are summed to obtain a compressed time frequency spectrum TFPC(td,fd) Wherein the coordinate (t)d,fd) Has an element value of
Figure GDA0002829571930000021
(34) Calculating a compressed time spectrum TFPC(td,fd) Maximum value of (P)cpThe compressed time spectrum TFPC(td,fd) Each element in (1) divided by PcpTo obtain a normalized spectrum NTFPC(td,fd);
(35) Calculating normalized spectra NTFPC(td,fd) Maximum value of (P)ncpAnd a minimum value VncpNormalized spectrum NTFPC(td,fd) Conversion to grey spectrum GTFPC(td,fd) Coordinate (t)d,fd) At a gray value of
Figure GDA0002829571930000022
Wherein [. cndot. ] represents a rounding symbol;
(36) will grey scale spectrum GTFPC(td,fd) Filling the state of health characteristic diagram FMap from line N to line N/2+ 1;
step four, acquiring frequency domain characteristic gray values of the vibration signals x (l) and filling the frequency domain characteristic gray values into the health state characteristic diagram FMap, wherein the method specifically comprises the following steps:
(41) performing fourier transform processing on x (L) to obtain a frequency spectrum x (f), wherein f is 1,2, …, L/2;
(42) equally dividing the frequency axis N of the frequency spectrum X (f) to obtain N frequency bands, wherein the frequency point number of each frequency band is delta f2(L/2)/N, and the kth frequency band is defined as X ((k-1) · Δ f)2+j2) Where k is 1,2, …, N, j2=1,2,…,△f2
(43) Calculating the energy ratio of each frequency band in the whole frequency spectrum X (f), wherein the energy ratio of the k-th frequency band is
Figure GDA0002829571930000031
(44) Calculating an energy ratio sequence Er(k) Maximum value of (P)ErAnd a minimum value VErComparing the energy to the sequence Er(k) Conversion to Gray spectra GEr(k),GEr(k) The gray value of is calculated by
Figure GDA0002829571930000032
Wherein [. cndot. ] represents a rounding symbol;
(45) gray scale GEr(k) Copying and filling into the lines 6 to N/2 of the health state characteristic diagram FMap, wherein the specific method is that the h line FMap (h, k) ═ GE of FMapr(k) Wherein h is 6,7, …, N/2;
step five, acquiring a time domain characteristic gray value of the vibration signal x (l) and filling the time domain characteristic gray value into the health state characteristic diagram FMap, wherein the method specifically comprises the following steps:
(51) 5 time domain dimensionless indexes of x (l) are calculated, and are respectively peak value indexes CFPulse index IFKurtosis index KFMargin index MFSum waveform index WFThereby obtaining a dimensionless index set Fnd={CF,IF,KF,MF,WF};
(52) F is to bendConverted into a time-domain index gray value set GFnd={GCF,GIF,GKF,GMF,GWFThe concrete method is that
Figure GDA0002829571930000033
Wherein iF=1,2,…,5;
(53) Gray value GF of time domain indexnd(iF) Copy and fill to (6-i) th of the health status feature map FMapF) In the middle of the line, i.e., FMap (6-i)F) All values of a line are GFnd(iF) Thus, a complete health status characteristic diagram FMap is obtained.
The invention has the following positive effects: in the health state characteristic diagram, time-frequency characteristic gray values depict time-frequency domain detail information of vibration signals, the frequency domain characteristic gray values represent frequency domain information of the vibration signals, and the time-domain characteristic gray values represent time-domain overall information of the vibration signals.
Drawings
FIG. 1 is a flow chart of an embodiment of the present invention;
FIG. 2 is a time domain waveform of vibration data of a rolling bearing according to an embodiment of the present invention;
fig. 3 is a health state characteristic diagram constructed by vibration data of the rolling bearing in the embodiment of the invention.
Detailed Description
The invention is further illustrated by the following figures and examples.
Referring to fig. 1, the method for constructing the health state characteristic diagram of the mechanical equipment comprises the following steps:
step one, presetting the dimension of a health state characteristic diagram FMap as NXN;
step two, collecting a vibration signal x (L) with the length of L on mechanical equipment, wherein L is 1,2, …, L;
step three, acquiring a time-frequency characteristic gray value of the vibration signal x (l) and filling the time-frequency characteristic gray value into a health state characteristic diagram FMap, wherein the method specifically comprises the following steps:
(31) performing time-frequency transform processing on x (L) to obtain a time spectrum TFP (t, f), wherein t is 1,2, …, L, f is 1,2, …, L/2;
(32) the time axis N of the time spectrum TFP (t, f) is equally divided, the time point number of each equal division is delta t equal to L/N, the frequency point number of each equal division is delta f for the frequency axis N/21For (L/2)/(N/2) L/N, so that for time spectrum TFP (t, f) the block is divided, defining the position (t)d,fd) The time-frequency block of (d) is TFP ((t)d-1)·△t+i,(fd-1)·△f1+j1) Wherein t isd=1,2,…,N,fd=1,2,…,N/2,i=1,2,…,△t,j1=1,2,…,△f1
(33) The element values in each time-frequency block are summed to obtain a compressed time frequency spectrum TFPC(td,fd) Wherein the coordinate (t)d,fd) Has an element value of
Figure GDA0002829571930000041
(34) Calculating a compressed time spectrum TFPC(td,fd) Maximum value of (P)cpThe compressed time spectrum TFPC(td,fd) Each element in (1) divided by PcpTo obtain a normalized spectrum NTFPC(td,fd);
(35) Calculating normalized spectra NTFPC(td,fd) Maximum value of (P)ncpAnd a minimum value VncpNormalized spectrum NTFPC(td,fd) Conversion to grey spectrum GTFPC(td,fd) Coordinate (t)d,fd) At a gray value of
Figure GDA0002829571930000051
Wherein [. cndot. ] represents a rounding symbol;
(36) will grey scale spectrum GTFPC(td,fd) Filling the state of health characteristic diagram FMap from line N to line N/2+ 1;
step four, acquiring frequency domain characteristic gray values of the vibration signals x (l) and filling the frequency domain characteristic gray values into the health state characteristic diagram FMap, wherein the method specifically comprises the following steps:
(41) performing fourier transform processing on x (L) to obtain a frequency spectrum x (f), wherein f is 1,2, …, L/2;
(42) equally dividing the frequency axis N of the frequency spectrum X (f) to obtain N frequency bands, wherein the frequency point number of each frequency band is delta f2(L/2)/N, and the kth frequency band is defined as X ((k-1) · Δ f)2+j2) Where k is 1,2, …, N, j2=1,2,…,△f2
(43) Calculating the energy ratio of each frequency band in the whole frequency spectrum X (f), wherein the energy ratio of the k-th frequency band is
Figure GDA0002829571930000052
(44) Calculating an energy ratio sequence Er(k) Maximum value of (P)ErAnd a minimum value VErComparing the energy to the sequence Er(k) Conversion to Gray spectra GEr(k),GEr(k) The gray value of is calculated by
Figure GDA0002829571930000053
Wherein [. cndot. ] represents a rounding symbol;
(45) gray scale GEr(k) Copying and filling into the lines 6 to N/2 of the health state characteristic diagram FMap, wherein the specific method is that the h line FMap (h, k) ═ GE of FMapr(k) Wherein h is 6,7, …, N/2;
step five, acquiring a time domain characteristic gray value of the vibration signal x (l) and filling the time domain characteristic gray value into the health state characteristic diagram FMap, wherein the method specifically comprises the following steps:
(51) 5 time domain dimensionless indexes of x (l) are calculated, and are respectively peak value indexes CFPulse index IFKurtosis index KFMargin index MFSum waveform index WFThereby obtaining a dimensionless index set Fnd={CF,IF,KF,MF,WF};
(52) F is to bendConverted into a time-domain index gray value set GFnd={GCF,GIF,GKF,GMF,GWFThe concrete method is that
Figure GDA0002829571930000061
Wherein iF=1,2,…,5;
(53) Gray value GF of time domain indexnd(iF) Copy and fill to (6-i) th of the health status feature map FMapF) In the middle of the line, i.e., FMap (6-i)F) All values of a line are GFnd(iF) Thus, a complete health status characteristic diagram FMap is obtained.
The invention is applied to the construction of the health state characteristic diagram of the SKF 6205-2RS type rolling bearing. The bearing vibration data is sourced from a bearing data center website of the university of western storage, and drive end bearing fault data is selected. The inner ring of the test bearing was machined with a dimple having a fault diameter of 0.18 mm. The load is set to be 0.75kW, the rotating speed is 1772r/min, the vibration data sampling frequency is 12k, and the steps of constructing the health state characteristic diagram of the rolling bearing by using the method are as follows.
Step one, the dimension of the preset health state characteristic diagram FMap is 32 multiplied by 32.
And step two, acquiring a vibration signal x (L) with the length L of 2048 on a mechanical device, wherein L is 1,2, … and 2048.
Step three, acquiring a time-frequency characteristic gray value of the vibration signal x (l) and filling the time-frequency characteristic gray value into a health state characteristic diagram FMap, wherein the method specifically comprises the following steps:
(31) performing short-time fourier transform on x (l) to obtain a time-frequency spectrum TFP (t, f), wherein t is 1,2, …,2048, f is 1,2, …, 1024;
(32) the time axis 32 of the time-frequency spectrum TFP (t, f) is equally divided, the number of time points per equal division is 64, the frequency axis 16 is equally divided, the number of frequency points per equal division is 64, and therefore the position (t, f) is defined for the TFP (t, f) blocksd,fd) The time-frequency block of (d) is TFP ((t)d-1)×64+i,(fd-1)×64+j1) Wherein t isd=1,2,…,32,fd=1,2,…,16,i=1,2,…,64,j1=1,2,…,64;
(33) The element values in each time-frequency block are summed to obtain a compressed time-frequency spectrum TFPC(td,ffd) Wherein the coordinate (t)d,fd) Has an element value of
Figure GDA0002829571930000071
(34) Calculating a compressed time spectrum TFPC(td,fd) Maximum value Pcp0.0092, the compressed time spectrum TFPC(td,fd) Each element in (1) divided by PcpTo obtain normalized time frequency spectrum NTFPC(td,fd);
(35) Calculating normalized spectra NTFPC(td,fd) Maximum value Pncp1 and minimum value VncpNormalized spectrum NTFP is normalized to 0.0021C(td,fd) Converting into time-frequency gray spectrum GTFPC(td,fd) Coordinate (t)d,fd) At a gray value of
Figure GDA0002829571930000072
Wherein [. cndot. ] represents a rounding symbol;
(36) time-frequency gray scale spectrum GTFPC(td,fd) Wherein t isd=1,2,…,32,f d1,2, …,16, filled into lines 32 through 17 of the state of health feature map FMap;
step four, acquiring frequency domain characteristic gray values of the vibration signals x (l) and filling the frequency domain characteristic gray values into the health state characteristic diagram FMap, wherein the method specifically comprises the following steps:
(41) fourier transforming x (l) to obtain a spectrum x (f), wherein f is 1,2, …, 1024;
(42) dividing the frequency axis 32 of the frequency spectrum X (f) equally to obtain 32 frequency bands, wherein the frequency point number of each frequency band is delta f2The kth band is X ((k-1) × 32+ j) × 322) Where k is 1,2, …,32, j2=1,2,…,32;
(43) Calculating the energy ratio of each frequency band in the whole frequency spectrum X (f), wherein the energy ratio of the k-th frequency band is
Figure GDA0002829571930000073
(44) Calculating an energy ratio sequence Er(k) Maximum value of (P)Er0.0471 and minimum value VEr0.0005, energy is compared to sequence Er(k) Conversion to frequency domain gray spectrum GEr(k),GEr(k) Has a gray value of
Figure GDA0002829571930000074
(45) GE frequency domain gray scale spectrumr(k) Copying and filling the data into lines 6 to 16 of the health state characteristic diagram FMap by lines, wherein the specific method is that the h line FMap (h, k) ═ GE of the health state characteristic diagram FMapr(k) Wherein h is 6,7, …, 16;
step five, acquiring a time domain characteristic gray value of the vibration signal x (l) and filling the time domain characteristic gray value into the health state characteristic diagram FMap, wherein the method specifically comprises the following steps:
(51) calculating 5 time domain dimensionless indexes of x (l), wherein the peak index CF4.8716 pulse index IF250.8961 kurtosis index KF5.4728, margin index MF8.6028 and waveform index WF51.5022, obtaining a dimensionless index set Fnd={4.8716,250.8961,5.4728,8.6028,51.5022};
(52) F is to bendConverted into a time-domain index gray value set GFnd={GCF,GIF,GKF,GMF,GWFThe concrete method is that
Figure GDA0002829571930000081
Wherein iF=1,2,…,5;
(53) Gray value GF of time domain indexnd(iF) Copy and fill to (6-i) th of the health status feature map FMapF) In the row, i.e., (6-i) of FMapF) All values of a line are GFnd(iF) Thus, a complete health status characteristic diagram FMap is obtained.
The embodiments described in this specification are merely illustrative of implementations of the inventive concept and the scope of the present invention should not be considered limited to the specific forms set forth in the embodiments but rather by the equivalents thereof as may occur to those skilled in the art upon consideration of the present inventive concept.

Claims (1)

1. The method for constructing the health state characteristic diagram of the mechanical equipment comprises the following steps:
step one, presetting the dimension of a health state characteristic diagram FMap as NXN;
step two, collecting a vibration signal x (L) with the length of L on mechanical equipment, wherein L is 1,2, …, L;
step three, acquiring a time-frequency characteristic gray value of the vibration signal x (l) and filling the time-frequency characteristic gray value into a health state characteristic diagram FMap, wherein the method specifically comprises the following steps:
(31) performing time-frequency transform processing on x (L) to obtain a time spectrum TFP (t, f), wherein t is 1,2, …, L, f is 1,2, …, L/2;
(32) the time axis N of the time spectrum TFP (t, f) is equally divided, the time point number of each equal division is delta t equal to L/N, the frequency point number of each equal division is delta f for the frequency axis N/21For (L/2)/(N/2) L/N, so that for time spectrum TFP (t, f) the block is divided, defining the position (t)d,fd) The time-frequency block of (d) is TFP ((t)d-1)·△t+i,(fd-1)·△f1+j1) Wherein t isd=1,2,…,N,fd=1,2,…,N/2,i=1,2,…,△t,j1=1,2,…,△f1
(33) The element values in each time-frequency block are summed to obtain a compressed time frequency spectrum TFPC(td,fd) Wherein the coordinate (t)d,fd) Has an element value of
Figure FDA0002829571920000011
(34) Calculating a compressed time spectrum TFPC(td,fd) Maximum value of (P)cpThe compressed time spectrum TFPC(td,fd) Each element in (1) divided by PcpTo obtain a normalized spectrum NTFPC(td,fd);
(35) Calculating normalized spectra NTFPC(td,fd) Maximum value of (P)ncpAnd a minimum value VncpNormalized spectrum NTFPC(td,fd) Conversion to grey spectrum GTFPC(td,fd) Coordinate (t)d,fd) At a gray value of
Figure FDA0002829571920000012
Wherein [. cndot. ] represents a rounding symbol;
(36) will gray scaleSpectrum GTFPC(td,fd) Filling the state of health characteristic diagram FMap from line N to line N/2+ 1;
step four, acquiring frequency domain characteristic gray values of the vibration signals x (l) and filling the frequency domain characteristic gray values into the health state characteristic diagram FMap, wherein the method specifically comprises the following steps:
(41) performing fourier transform processing on x (L) to obtain a frequency spectrum x (f), wherein f is 1,2, …, L/2;
(42) equally dividing the frequency axis N of the frequency spectrum X (f) to obtain N frequency bands, wherein the frequency point number of each frequency band is delta f2(L/2)/N, and the kth frequency band is defined as X ((k-1) · Δ f)2+j2) Where k is 1,2, …, N, j2=1,2,…,△f2
(43) Calculating the energy ratio of each frequency band in the whole frequency spectrum X (f), wherein the energy ratio of the k-th frequency band is
Figure FDA0002829571920000021
(44) Calculating an energy ratio sequence Er(k) Maximum value of (P)ErAnd a minimum value VErComparing the energy to the sequence Er(k) Conversion to Gray spectra GEr(k),GEr(k) The gray value of is calculated by
Figure FDA0002829571920000022
Wherein [. cndot. ] represents a rounding symbol;
(45) gray scale GEr(k) Copying and filling into the lines 6 to N/2 of the health state characteristic diagram FMap, wherein the specific method is that the h line FMap (h, k) ═ GE of FMapr(k) Wherein h is 6,7, …, N/2;
step five, acquiring a time domain characteristic gray value of the vibration signal x (l) and filling the time domain characteristic gray value into the health state characteristic diagram FMap, wherein the method specifically comprises the following steps:
(51) 5 time domain dimensionless indexes of x (l) are calculated, and are respectively peak value indexes CFPulse index IFKurtosisDegree index KFMargin index MFSum waveform index WFThereby obtaining a dimensionless index set Fnd={CF,IF,KF,MF,WF};
(52) F is to bendConverted into a time-domain index gray value set GFnd={GCF,GIF,GKF,GMF,GWFThe concrete method is that
Figure FDA0002829571920000023
Wherein iF=1,2,…,5;
(53) Gray value GF of time domain indexnd(iF) Copy and fill to (6-i) th of the health status feature map FMapF) In the middle of the line, i.e., FMap (6-i)F) All values of a line are GFnd(iF) Thus, a complete health status characteristic diagram FMap is obtained.
CN201910860506.2A 2019-09-11 2019-09-11 Method for constructing health state characteristic diagram of mechanical equipment Active CN110595811B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910860506.2A CN110595811B (en) 2019-09-11 2019-09-11 Method for constructing health state characteristic diagram of mechanical equipment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910860506.2A CN110595811B (en) 2019-09-11 2019-09-11 Method for constructing health state characteristic diagram of mechanical equipment

Publications (2)

Publication Number Publication Date
CN110595811A CN110595811A (en) 2019-12-20
CN110595811B true CN110595811B (en) 2021-04-09

Family

ID=68858904

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910860506.2A Active CN110595811B (en) 2019-09-11 2019-09-11 Method for constructing health state characteristic diagram of mechanical equipment

Country Status (1)

Country Link
CN (1) CN110595811B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113822257B (en) * 2021-11-24 2022-06-07 航天智控(北京)监测技术有限公司 Abnormal point detection method based on combination of dimensionless features and virtual samples

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102222324B (en) * 2011-06-17 2012-08-29 电子科技大学 Symmetry property-based method for detecting salient regions of images
US10762894B2 (en) * 2015-03-27 2020-09-01 Google Llc Convolutional neural networks
CN105654037B (en) * 2015-12-21 2019-05-21 浙江大学 A kind of electromyography signal gesture identification method based on deep learning and characteristic image
US10334357B2 (en) * 2017-09-29 2019-06-25 Apple Inc. Machine learning based sound field analysis
CN108564005B (en) * 2018-03-26 2022-03-15 电子科技大学 Human body falling identification method based on convolutional neural network
CN108764601B (en) * 2018-04-03 2021-07-16 哈尔滨工业大学 Structural health monitoring abnormal data diagnosis method based on computer vision and deep learning technology
CN109934206A (en) * 2019-04-08 2019-06-25 中国矿业大学(北京) A kind of rotary machinery fault diagnosis method under non-stationary operating condition
CN109994203B (en) * 2019-04-15 2023-03-24 江南大学 Epilepsia detection method based on EEG signal depth multi-view feature learning
CN110197127B (en) * 2019-05-06 2022-10-18 安徽继远软件有限公司 Wireless signal detection and electromagnetic interference classification system and method based on deep learning
CN110108431B (en) * 2019-05-22 2021-07-16 西安因联信息科技有限公司 Mechanical equipment fault diagnosis method based on machine learning classification algorithm

Also Published As

Publication number Publication date
CN110595811A (en) 2019-12-20

Similar Documents

Publication Publication Date Title
Wang et al. A deep learning method for bearing fault diagnosis based on time-frequency image
CN111160167A (en) Spindle fault classification and identification method based on S-transform deep convolutional neural network
CN106124212B (en) Fault Diagnosis of Roller Bearings based on sparse coding device and support vector machines
CN107702922B (en) Rolling bearing fault diagnosis method based on LCD and stacking automatic encoder
CN112633245B (en) Planetary gear box fault diagnosis method based on deep reinforcement learning model
CN112418013A (en) Complex working condition bearing fault diagnosis method based on meta-learning under small sample
CN112232244A (en) Fault diagnosis method for rolling bearing
CN110702411A (en) Residual error network rolling bearing fault diagnosis method based on time-frequency analysis
Liang et al. Intelligent fault diagnosis of rolling element bearing based on convolutional neural network and frequency spectrograms
CN107024271A (en) Mechanical oscillation signal compression reconfiguration method and system
CN113780412A (en) Fault diagnosis model training method and system and fault diagnosis model training method and system
Wang et al. An end-to-end fault diagnostics method based on convolutional neural network for rotating machinery with multiple case studies
CN107966287B (en) Weak fault feature extraction method for self-adaptive electromechanical equipment
CN110595811B (en) Method for constructing health state characteristic diagram of mechanical equipment
CN113642414A (en) Method for predicting residual service life of rolling bearing based on Transformer model
CN111610024A (en) Gear crack fault diagnosis method based on vibration signals
Zhang et al. Compound Fault diagnosis for gearbox based using of euclidean matrix sample entropy and one-dimensional convolutional neural network
CN116465623B (en) Gearbox service life prediction method based on sparse converter
CN113408578A (en) Harmonic reducer fault diagnosis method and system based on generative countermeasure network
CN115808312A (en) Rolling bearing fault diagnosis method based on deformable volume and Transformer
Zhong et al. Parallel learning attention-guided CNN for signal denoising and mechanical fault diagnosis
CN114462446A (en) Rolling bearing fault diagnosis method based on vibration signal and electronic equipment
Liu et al. Sinc-based multiplication-convolution network for equipment intelligent edge diagnosis under small samples
CN117708574B (en) CNN variable-speed rolling bearing fault diagnosis method embedded with physical information
Song et al. Bearing fault diagnosis method using a spatio-temporal neural network based on feature transfer learning

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant