Measuring the speed of rotation or reciprocation of a mechanical component using one or more cameras
11373317 · 2022-06-28
Assignee
Inventors
- Jeffrey R. Hay (Prospect, KY, US)
- Kenneth Ralph Piety (Knoxville, TN, US)
- Mark William Slemp (Tellico Plains, TN, US)
- Joseph A. Vrba (Clinton, TN)
Cpc classification
G06T7/246
PHYSICS
International classification
Abstract
Present embodiments pertain to systems, apparatuses, and methods for analyzing and reporting rotational or reciprocating movements in mechanical structures, machines, and machine components, including measuring the speed of rotation or reciprocation of a component on the structure, through the use of an acquired video recording representing a plurality of cycles of motion, by measuring intensity values of a subset of pixels contained in a region of interest within the video recording in a plurality of frames of the video recording, thereby determining distance per time period of rotational or reciprocal motion of a machine or machine component, and the application of numerical algorithm to the repeating patterns in the intensity waveform enables the determination of the average speed value or with the use of a gear wheel or graduated tape, instantaneous values for speed can be derived and the torsional vibration characteristics of the component determined.
Claims
1. A method for calculating a speed of a machine or machine component which is rotating or reciprocating using a video recording acquired of a scene that includes the machine or machine component, comprising: storing the video recording in a memory operatively connected to a processor that executes computer-readable program instructions, wherein the video recording comprises a plurality of video frames; configuring the processor to receive user input through a graphical user interface (GUI) to identify a region of interest (ROI) in the scene which includes at least a portion of the rotating or reciprocating machine or machine component; configuring the processor to extract a waveform representing motion of the machine or machine component, wherein the waveform is based on a plurality of cycles as the machine or machine component rotates or reciprocates; and calculating and displaying the speed of the machine or machine component as an average of speed values for at least one subset of pixels contained in the ROI, wherein each of the pixels in the at least one subset of pixels represents an intensity of light reflected from the machine or machine component exhibiting the periodicity of the rotational or reciprocal motion of the machine or machine component.
2. The method of claim 1, further comprising configuring the processor to receive an approximated speed of the machine or machine component that is input by the user.
3. The method of claim 1, wherein the processor is not provided an estimate of the speed of the machine or machine component and the processor applies an autocorrelation algorithm to one or more possible speed intervals until it locates a speed interval containing the true speed of the machine or machine component.
4. The method of claim 1, wherein the processor is not provided an estimate of the speed of the machine or machine component and the processor applies an FFT algorithm to one or more possible speed intervals until it locates a speed interval containing the true speed of the machine or machine component.
5. The method of claim 1, further comprising configuring the processor to allow the user to accept or reject the speed of the machine or machine component.
6. The method of claim 1, further comprising configuring the processor to automatically adjust at least one of a field of view of the video recording, camera exposure, a frame rate of the camera, and a duration of the video recording.
7. The method of claim 1, wherein calculating the average speed values for the at least one subset of pixels comprises: analyzing a periodicity of the intensity of light reflected in a plurality of pixels associated with the rotating or reciprocating motion of the machine component by determining the spacing of a plurality of the largest peaks in an autocorrelation function of the waveform; determining a coefficient of variance of the separation of the peaks in the autocorrelation function; and retaining the speed value obtained from one or more retained pixels that are less than a predetermined threshold for coefficient of variance and using said retained speed values to calculate the average speed values for the at least one subset of pixels.
8. The method of claim 7, wherein the number of the largest peaks is at least 6.
9. The method of claim 7, wherein the coefficient of variance is less than 5.0.
10. The method of claim 1, wherein calculating the average speed values for the at least one subset of pixels comprises: analyzing a periodicity of the intensity of light reflected in a plurality of pixels associated with the rotating or reciprocating motion of the machine component by determining the spacing of a plurality of the largest peaks in the FFT function of the waveform; locating a set of harmonically related peaks in the set of largest peaks in the spectrum; determining the fundamental frequency of this harmonic family and retaining this as the speed value obtained from one or more retained pixels; and using said retained speed values to calculate the average speed values for the at least one subset of pixels.
Description
BRIEF DESCRIPTION OF DRAWINGS
(1) The drawings, schematics, figures, and descriptions contained in this application are to be understood as illustrative of steps, structures, features and aspects of the present embodiments. Accordingly, the scope of embodiments is not limited to features, dimensions, scales, and arrangements shown in the figures.
(2)
(3)
(4)
(5)
(6)
(7)
(8)
(9)
(10)
(11)
(12)
(13)
(14)
(15)
(16)
(17)
(18)
(19)
(20)
MULTIPLE EMBODIMENTS AND ALTERNATIVES
(21) In some aspects of the present disclosure, in accordance with multiple embodiments and alternatives, a user in the operation of present embodiments measures the speed of rotation or reciprocation of an element on a mechanical structure, a machine or machine component without contacting the structure using the same camera that will be employed to investigate the dynamic motion of the mechanical structure.
(22) In
(23) Machines often present a complex picture on a video recording that represents the operating speed of the machine itself, as well as potentially different operating speeds of components of the machine as they undergo periodic movements (rotational or reciprocating). By way of simple and non-limiting illustration, a machine might include a motor with components captured on the video rotating at first frequency, operatively connected to a rotating shaft that drives one or more gears through a series of rotations at a second frequency. Therefore, it is extremely important to know the exact frequency of rotation or reciprocation associated with the current operation of a machine. This allows an analyst to identify other frequencies present in the vibration measured on the mechanical structure. When diagnosing a specific fault condition, it is important to know which frequencies on a mechanical structure are sub-synchronous, synchronous, or non-synchronous with respect to the operation speed. A synchronous peak which occurs at 12 times the running speed might be associated with a 12-tooth gear; however, a nonsynchronous peak occurring at 12.08 times running speed may be associated with a defect in an anti-friction bearing. The ability to measure the rate of rotation or reciprocation using the same non-contacting camera that will make the vibration measurements on the structure is very cost effective and efficient for the analyst.
(24) In one or more exemplary embodiments, the user follows the typical steps necessary to obtain good recordings of a machine or portion of machine having a rotating or reciprocating component:
(25) 1. Position the camera to acquire the perspective of the equipment of interest, containing at least a portion of the rotating or reciprocating component,
(26) 2. Focus the camera, and
(27) 3. Adjust the aperture; this may require the addition of external light or shielding the field of view in the presence of bright light conditions to achieve acceptable lighting conditions for recording.
Example—Measuring Speed of Asset
(28) After the camera is set up as desired or needed (
Example—Automatically Adjusting Exposure of Camera Before Measuring Asset Speed
(29) The system software in some embodiments can automatically adjust the exposure to improve the asset speed measurement by modifying the brightness and gain settings on the camera. One embodiment providing an automated exposure adjustment method is outlined in the flowchart provided in
(30) 1. Gain is initialized to 0 (
(31) 2. A frame is acquired and the percentage of pixels above half intensity (i.e. 2048 for 12-bit data) is computed.
(32) 3. If less than 10% of pixels are above half intensity (
(33) a. If Brightness is <95%, then Brightness is increased by 20% (limited to 100%) and then repeat from “2” above (
(34) b. Otherwise, if Gain is <95%, then Gain is increased by 10% (limited to 100%) and then repeat from “2” above (
(35) 4. If more than 10% of pixels are above half intensity, then exposure is accepted (
(36) Automated adjustment of the exposure minimizes user interaction and the time required to perform the speed measurement. Any other algorithms which result in an acceptable exposure could be employed in alternate embodiments and would remain in the scope of the embodiments described herein. If the system cannot successfully make the speed measurement, then the user will need to select another location on the shaft or manually adjust the exposure.
Example—Autocorrelation Algorithm Facilitating Measurement of Asset Speed
(37) In some embodiments, after the exposure has been automatically adjusted, then a recording is acquired which will provide waveforms with approximately 12 revolutions of the shaft (128-1024 samples per waveform). For a set of N pixels (camera pixels or virtual pixels) (
(38) 1. Intensity waveforms are formed, and “DC” is removed from them (
(39) 2. The autocorrelation waveform is computed (
(40) 3. The largest peaks in the autocorrelation waveform are located (
(41) 4. The peaks are sorted in time order and statistics of their spacing computed. If the coefficient of variance of the peaks is >−5.0, the pixel is discarded (
(42) 5. The average peak spacing yields the speed from that pixel (
(43) 6. The asset speed is average of the retained pixel speeds ((
(44) The resulting asset speed is displayed to the user who then accepts or rejects it. A flowchart of this algorithm is shown in
Example—Peak Location Algorithm Based on Frequency Spectrum
(45) In an alternate embodiment, the practice of which is described below in exemplary fashion, the asset speed is determined by locating the peaks in FFT frequency spectrum of the waveforms with approximately 12 revolutions of the shaft (128-1024 samples per waveform). For a set of N pixels (camera pixels or virtual pixels) (
(46) 1. Intensity waveforms are formed, and “DC” is removed from them (
(47) 2. The FFT frequency spectrum of the intensity waveform is computed (
(48) 3. The largest 5 peaks in the FFT spectrum are located accurately using FFT windowing parameters (
(49) 4. The peaks are tested to find harmonically related peaks and the fundamental frequency of the family is calculated (
(50) 5. The mean value of the retained fundamental harmonic frequency values is formed and any pixels whose value fundamental frequency differs by more than 1 Hz is discarded from the set,
(51) 6. The asset speed is average of the fundamental harmonic frequencies for the retained pixels (
(52) The resulting asset speed is displayed to the user who then accepts or rejects it. An exemplary flowchart of this algorithm is shown in
(53) This process is improved by using located peak values rather than the nominal peak value which is identified by the frequency line in the spectrum with the highest amplitude. In a spectrum calculated with 1 Hz resolution, the nominal peak value might be 22 Hz because the frequency line at 22 Hz has the highest amplitude of 2. The amplitude value at 21 Hz might be 0.05 and the line at 23 Hz might be 1.90 indicating that the true peak frequency lies between 22 Hz and 23 Hz. It is well known in signal processing art, that the true frequency value can be estimated more accurately by applying formulas that consider the windowing function used when calculating the FFT frequency spectrum. In the case above, the true value would be about halfway between the two lines giving a located peak frequency of 22.4 Hz. The FFT could be constructed using any number of windows such as the Uniform, Hanning, Hamming, Blackman-Harris, Kaiser-Bessel, or others. More accurate frequency estimates of the peak location can be calculated using the parameters that are characteristic of the respective windows. An improved location of the peak frequency can also be accomplished by applying any number of well-known fitting algorithms to the center line in the peak and the 2 lines on either side. The generic fitting algorithms are generally not as accurate as using the algorithm that takes into account the FFT windowing functions. Demonstrating the broad nature of the descriptions herein, alternative embodiments could use any of the methods discussed or those obtaining equivalent improvements in locating the peak frequency values. Nominal peak frequency values can also be used but would not provide as reliable results in some cases. When calculating the ratio of the frequency of two peaks to determine if they are harmonically related (integer multiples of the fundamental frequency), located peak values will result in a closer match to integer values and thus correctly identify harmonic family members.
(54) The intensity waveform from one of the pixels in the region of interest located on the rotating or reciprocating component is shown in
(55) In an alternate embodiment, the algorithms described above could be applied to the regular video recording which are captured to measure the orthogonal vibration in the field of view. In the situation where the rotating or reciprocating element is in the field of view and the data acquisition parameters are sufficient to enable resolving the rate of rotational or reciprocation, the rate of repetition could be measured for a ROI identified by the user.
(56) In other aspects of the present disclosure, in accordance with multiple embodiments and alternatives, a user in the operation of present embodiments measures the instantaneous speed of rotation of an element on a mechanical structure, a machine or machine component without contacting the structure using the same camera that will be employed to investigate the dynamic motion of the mechanical structure. Measurement of the instantaneous speed at various points during each revolution of the shaft allows the user to investigate the angular or torsional vibration characteristics of the mechanical system. This measurement requires access to a gear wheel which is integral or attached to the shaft of interest or the application of a patterned or graduated tape to a visible position on the shaft.
(57) When using a video camera to make torsional measurements, the tape applied could have graduated lines or a repeating pattern, such as shown in
(58) In a use case shown in
(59) 1. On a selected frame, the user optionally identifies the edge of component by drawing a reference line (212 in
(60) 2. On the selected frame, the user identifies the region where the torsional motion will be measured by drawing the reference rectangle (identified as 211 in
(61) 3. The system software can optionally match the references to the side of the disk and the closest lines of the selected division to achieve a finer resolution in the position of user selected references.
(62) 4. Count number of whole divisions (N) and measure fraction division (F), visually or from video,
(63) 5. Determine circumference (C) and diameter (D) of component being monitored,
(64) 6. Set the sampling rate to be the larger of:
(65) a. 1.5*N.Math.F*RPS or
(66) b. 2.5 times the number of orders of interest,
(67) 7. Collect the number of samples or frames equal to 3 times the reciprocal of the lowest sub order (1/M) of interest times the sampling rate: Total Samples=3*SR*M,
(68) 8. Determine change in arc for the division passing through the reference rectangle for each frame (sample) recorded:
(69) a. Use the fractional movement of the lines across this reference zone if processing a full tape division,
(70) b. If a partial division (tape overlap) enters the zone, then move to preceding full tape division to determine distance moved in the interval or the arc of the motion,
(71) 9. Movement of the wheel due to vibration can be removed by measuring the motion at the reference line (212 in
(72) Further, in accordance with the embodiments herein,
(73) It will be appreciated that this embodiment describes the use of targets or markings affixed to the shaft. In some instance measurements may be possible on pre-existing markings on the shaft. These marking may be from but not limited to normal wear on the shaft, damage to the shaft, scuff or scratches or marking made on the shaft by an individual.
(74) In some instances, the displacements of a marking on the shaft will be made directly from the apparent distance traveled by the mark as seen by the camera. It will be appreciated that a shaft measurement can be corrected for the curvature of the shaft such that apparent motion on the horizon of the shaft can be compared to the motion at the nearest point of the shaft to the camera. When a mark is moving near the horizon of the shaft a given angular displacement will appear to travel a shorter distance as viewed from the camera as compared to the same angular displacement moving at the nearest point of the shaft to the camera. By accounting for this curvature, the angular displacement can be normalized and compared across the entire shaft. The radius of the shaft may be entered into the software to make this correction or the camera may make the measurement of the shaft radius from the image itself if the shaft is visible and the scale of the image is known.
(75) Several exemplary claims are set forth herein but are not intended to place boundaries on the full range of embodiments and alternatives described and provided for herein, nor are these intended to waive or otherwise circumscribe any potential claims that could be pursued in a later application claiming the benefit of the teachings and disclosures herein. A claim expressed in this filing as a method also could represents a system that performs such a method or other methods, and a system claim if recited also could represent a method of operation executed by said system.
(76) It will be understood that the embodiments described herein are not limited in their application to the details of the teachings and descriptions set forth, or as illustrated in the accompanying figures. Rather, it will be understood that the present embodiments and alternatives, as described and claimed herein, are capable of being practiced or carried out in various ways. Also, it is to be understood that words and phrases used herein are for the purpose of description and should not be regarded as limiting. The use herein of such words and phrases as “including,” “such as,” “comprising,” “e.g.,” “containing,” or “having” and variations of those words is meant to encompass the items listed thereafter, and equivalents of those, as well as additional items.
(77) Accordingly, the foregoing descriptions of embodiments and alternatives are meant to illustrate, rather than to serve as limits on the scope of what has been disclosed herein. The descriptions herein are not meant to limit the understanding of the embodiments to the precise forms disclosed. It will be understood by those having ordinary skill in the art that modifications and variations of these embodiments are reasonably possible in light of the above teachings and descriptions.