CN103267567A - Measuring device and method for vibration of flexible cantilever on basis of machine vision - Google Patents
Measuring device and method for vibration of flexible cantilever on basis of machine vision Download PDFInfo
- Publication number
- CN103267567A CN103267567A CN2013102248219A CN201310224821A CN103267567A CN 103267567 A CN103267567 A CN 103267567A CN 2013102248219 A CN2013102248219 A CN 2013102248219A CN 201310224821 A CN201310224821 A CN 201310224821A CN 103267567 A CN103267567 A CN 103267567A
- Authority
- CN
- China
- Prior art keywords
- flexible cantilever
- led luminotron
- image
- cantilever beam
- camera
- 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.)
- Granted
Links
Images
Landscapes
- Length Measuring Devices By Optical Means (AREA)
Abstract
The invention discloses a measuring device and method for the vibration of a flexible cantilever on the basis of machine vision. The measuring device comprises a fixed bracket, the flexible cantilever, a camera bracket, a CCD (Charge Coupled Device) camera, a lens, a plurality of LED (Light-Emitting Diode) light-emitting tubes and a PC (Personal Computer), wherein one end of the flexible cantilever is fixed on the fixed bracket; all the LED light-emitting tubes are arranged on the upper surface of the flexible cantilever in sequence; the top of the fixed bracket is provided with the camera bracket; the CCD camera is fixed on the camera bracket and is provided with the lens; and an output port of the CCD camera is connected with the PC. The CCD camera measures each frame image of the vibration of the LED light-emitting tubes, and conveys the images to the PC; and the PC processes a detected image sequence, extracts the position of the mass center of the light spots of the LED light-emitting tubes and acquires the vibration displacement of all the LED light-emitting tubes and parameters for reflecting the low-frequency vibration of the flexible cantilever structure. The measuring device and method disclosed by the invention has the advantages that the non-contact effect is achieved, the measuring range is wide and the vibration characteristic of a measured object is not changed and the like, and can be widely applied.
Description
Technical field
The present invention relates to the measuring method of large-scale flexible semi-girder vibration, particularly a kind of measurement mechanism and method of the flexible cantilever beam vibration based on machine vision.
Background technology
Solar energy sailboard is a kind of collection solar device, is generally used for the energy supply of satellite, spaceship.The modal damping of this class large-scale flexible structure is little, in case be subjected to the effect of certain exciting force, it significantly vibrates last very long.This not only can influence the work of space structure, also will make structure produce too early fatigure failure, influences the serviceable life of structure, or causes the damage of instrument in the structure.Therefore need to measure and suppress the vibration of this type of spacecraft flexible structure, particularly at the structural vibration of the flexible cantilever beam of analog solar windsurfing.
In the prior art, the vibration survey of the flexible cantilever beam of analog solar windsurfing mainly contains piezoelectric patches, acceierometer sensor, angular rate gyroscope sensor, photoelectrical position sensor (Position Sensitive Detector) and the fiber-optic grating sensor methods such as (Fiber Grating Sensor) of adopting: when piezoelectric patches is measured because multiple links such as signal amplifications, filtering can cause signal delay and phase place hysteresis, measuring speed slowly, be limited in scope; Acceierometer sensor and angular rate gyroscope sensor are to noise-sensitive, exist sluggish and temperature such as floats at its precision of problems affect, and can only measure on the object certain any displacement, want to obtain full detail, must on semi-girder, arrange several accelerometers or angular rate gyroscope respectively; The measurement range of photoelectrical position sensor is less, complex structure, operation easier is big, computation process is loaded down with trivial details and cost is expensive; The topmost problem of fiber-optic grating sensor is the demodulation of transducing signal, because fiber grating is relatively more fragile, is very easy to destroy in abominable working environment, thereby could uses packaging technology and safeguard measure complex structure after need encapsulating it.
Summary of the invention
Measurement mechanism and the method based on the flexible cantilever beam vibration of machine vision that the purpose of this invention is to provide a kind of simple in structure, stable performance, efficiently and accurately, realize to the flexible cantilever beam vibration displacement noncontact, accurately measure in real time.
The technical solution that realizes the object of the invention is, a kind of measurement mechanism of the flexible cantilever beam vibration based on machine vision, comprise fixed support, flexible cantilever beam, this measurement mechanism also comprises camera support, CCD camera, camera lens, a plurality of LED luminotron and PC; Described flexible cantilever beam one end is fixed on the fixed support, and the other end is free end, and each LED luminotron is along the upper surface that is successively set on flexible cantilever beam from the fixed support end to free-ended direction; The fixed support top is provided with camera support, and the CCD camera is fixed on the camera support, and the CCD camera configuration has camera lens, and the output port of CCD camera is connected with PC.
Measuring method based on the flexible cantilever beam vibration of machine vision is characterized in that step is as follows:
The 1st step: according to the parameter of CCD camera, each LED luminotron is tied to the horizontal transformation factor K of world coordinate system from image coordinate
x, vertical transformation factor K
yDemarcate;
The 2nd step: the CCD camera is gathered the vibrational image sequence of flexible cantilever beam, and exports the vibrational image sequence to PC by data line;
The 3rd step: according to the position of each LED luminotron in first two field picture, entire image is divided into the subregion corresponding with each LED luminotron, and determines the region of interest ROI of each subregion;
The 4th step: adopt the OTSU thresholding method to extract the hot spot of the interior corresponding LED luminotron of each region of interest ROI of flexible cantilever beam, and determine the centroid position of each hot spot;
The 5th step: according to the centroid position of each hot spot, obtain horizontal shift Δ x and the vertical displacement delta y of corresponding LED luminotron in the image;
The 6th step: according to horizontal shift Δ x, the vertical displacement delta y of each LED luminotron in the image and the horizontal transformation factor K of the 1st step demarcation
x, vertical transformation factor K
y, obtain the real standard displacement x of each LED luminotron, actual perpendicular displacement y;
The 7th step: according to the acquisition order of the 2nd step vibrational image sequence, draw out the position curve of each LED luminotron, Fast Fourier Transform (FFT) is carried out in the displacement of each LED luminotron, obtain single order model frequency and the second-order modal frequency of flexible cantilever beam vibration.
The present invention compared with the prior art, have following remarkable advantage: (1) utilizes Machine Vision Detection can be easy to, obtain intuitively the low-frequency vibration information of flexible cantilever beam, the measurement noise is little, and adopts the simple, convenient enforcement of light source detection of following the tracks of the LED luminotron; (2) utilize vision-based detection that the advantage of multimetering is arranged, can measure multi-point displacement and deformation, real-time and accuracy height; (3) by picture portion, by each LED luminotron of row labels, choose suitable interesting image regions ROI rapidly, reduced image deal with data amount, greatly promoted the speed of image characteristics extraction; (4) under the situation of the vibration characteristics that does not change measured object, follow the tracks of light source by machine vision, detect the multiple spot dynamic displacement of flexible structure, have the advantages such as vibration characteristics that noncontact, measurement range are wide, do not change measured object.
Description of drawings
Fig. 1 is the structural representation of measurement mechanism that the present invention is based on the flexible cantilever beam vibration of machine vision.
Fig. 2 is the vertical view of flexible overarm arm upper surface.
Fig. 3 is CCD camera imaging principle schematic.
Fig. 4 is the level of vibration displacement curve figure of the 3rd LED luminotron in the embodiment of the invention.
Fig. 5 is the rumble spectrum figure of the 3rd LED luminotron in the embodiment of the invention.
Embodiment
The present invention is further illustrated below in conjunction with the drawings and specific embodiments.The present invention obtains the vibration displacement of flexible cantilever beam and the devices and methods therefor that carries out the low frequency modal vibration analysis for the method that adopts machine vision by the trace labelling light source.
As shown in Figure 1, the measurement mechanism based on the flexible cantilever beam vibration of machine vision comprises fixed support 1, flexible cantilever beam 5, camera support 2, CCD camera 3, camera lens 4, a plurality of LED luminotron 6 and PC 7; Described flexible cantilever beam 5 one ends are fixed on the fixed support 1, and the other end is free end, and each LED luminotron 6 is along the upper surface that is successively set on flexible cantilever beam 5 from the fixed support end to free-ended direction; Fixed support 1 top is provided with camera support 2, and CCD camera 3 is fixed on the camera support 2, and CCD camera 3 disposes camera lens 4, and the output port of CCD camera 3 is connected with PC 7.Described each LED luminotron equidistantly is fixed on the center line of flexible cantilever beam 5 upper surfaces, and all LED luminotrons 6 are all in the visual field of CCD camera 3.The quantity of LED luminotron 6 is determined according to the size of flexible cantilever beam 5 and the mode number of needs consideration and control.
Fig. 2 is the upper surface vertical view of flexible cantilever beam 5, comprises that flexible cantilever beam upper surface, flexible cantilever beam upper surface center line, CCD camera optical axis, optical axis are at the projection line of property semi-girder upper surface.The position of CCD camera 3 relative flexibility semi-girders 5 stiff ends and attitude are determined, the CCD plane of CCD camera 3 and the upper surface of flexible cantilever beam 5 13 angled α, α determines according to the visual angle of camera, the height that camera is placed, the length of flexible cantilever beam, the α scope is 0 °~90 °, and all LED luminotrons 6 are in the field range of CCD camera 3 all the time when making flexible cantilever beam 5 vibrations; As shown in Figure 2,15 one-tenth 90 °-α of flexible cantilever beam upper surface 13 and CCD camera optical axis angle when static, and CCD camera optical axis and last plane center line intersect, CCD camera optical axis at the projection line on plane on the flexible cantilever beam on the center line of last plane.
Based on the measuring method of the flexible cantilever beam vibration of machine vision, step is as follows:
The 1st step: according to the parameter of CCD camera 3, each LED luminotron 6 is tied to the horizontal transformation factor K of world coordinate system from image coordinate
x, vertical transformation factor K
yDemarcate, concrete steps are as follows:
(1.1) with the camera model approximate expression actual imaging of no perspective distortion, as shown in Figure 3, set up coordinate system: with the mid point O of flexible cantilever beam 5 upper surface fixed edges
wBe initial point, set up world coordinate system O
w-X
wY
wZ
wPhotocentre O with CCD camera 3
CBe initial point, setting up camera coordinates is O
c-X
cY
cZ
cBe initial point O with first pixel of the image upper left corner, set up image coordinate system O-UV; Central point O with image
dBe initial point, setting up the image planes physical coordinates is O
d-X
dY
d, O wherein
dCoordinate be (c
x, c
y); The optical axis of CCD camera 3 and X
cO
cY
c, UOV, X
dO
dY
dThe plane is all vertical, the angled α of upper surface of the CCD imaging plane of CCD camera 3 and flexible cantilever beam 5, the optical axis of CCD camera 3 and X
WO
WY
WThe angled 90 °-α in plane, wherein X
WAxle, X
CAxle, X
dAxle is parallel, and CCD camera 3 positions and attitude are fixed, camera coordinates system and world coordinate system only exist around X-axis rotate, translation relation;
(1.2) world coordinate system O
w-X
wY
wZ
wBe O to camera coordinates
c-X
cY
cZ
cTransformation relation be:
(x wherein
c, y
c, z
c) for camera coordinates be O
c-X
cY
cZ
cCoordinate, (x
w, y
w, z
w) be world coordinate system O
w-X
wY
wZ
wCoordinate,
In the formula, R is 3 * 3 rotational transform matrix, and θ is the anglec of rotation around X-axis; T is 3 * 1 translation vector, and b is shift value;
(1.3) camera coordinates is O
c-X
cY
cZ
cBe O to the image planes physical coordinates
d-X
dY
dCoordinate transform be:
In the formula, (x
d, y
d) for the image planes physical coordinates be O
d-X
dY
dCoordinate, f is that effective focal length is the distance of the plane of delineation and optical center of lens;
(1.4) the image planes physical coordinates is O
d-X
dY
dTo being transformed to of image coordinate system O-UV:
In the formula, (u v) is the coordinate of image coordinate system O-UV, d
xBe in the image pixel along the physical size of x axle, d
yBe in the image pixel along the physical size of y axle;
(1.5) composite type (1)~(3) are by world coordinate system O
w-X
wY
wZ
wHomogeneous coordinates to image coordinate system O-UV conversion are expressed as:
In the formula, S is scale-up factor, and formula (4) is transformed to:
(5)
Because LED luminotron 6 is only on the plane
Internal vibration, so following formula (5) is deformed into:
Get image coordinate system O-UV to world coordinate system O by formula (6)
w-X
wY
wZ
wThe horizontal transformation factor K
xWith vertical transformation factor K
y, be shown below respectively:
In the formula, d
x, d
y, f, c
x, c
yBe CCD camera inner parameter, θ, z
cBe external parameter, inner parameter is relevant with the camera inner structure, and external parameter obtains by demarcation.
The 2nd step: CCD camera 3 is gathered the vibrational image sequence of flexible cantilever beams 5, and exports the vibrational image sequence to PC 7 by data line; By the PC program to the vibrational image sequence identification of being correlated with handle the displacement of obtaining LED luminotron 6 places on the flexible cantilever beam 5.
The 3rd step: according to the position of each LED luminotron 6 in first two field picture, entire image is divided into the subregion corresponding with each LED luminotron 6, and determines the region of interest ROI of each subregion; According to the maximum perpendicular displacement of LED luminotron 6 each subregion is chosen the region of interest ROI of image, reduce the data volume of handling, the LED luminotron mainly is horizontal displacement x, and y is less for the perpendicular displacement amount.
The 4th step: adopt the OTSU thresholding method to extract the hot spot of the interior corresponding LED luminotron 6 of each region of interest ROI of flexible cantilever beam 5, and determine the centroid position of each hot spot; The OTSU thresholding method is a kind of method of determining threshold value automatically that makes the inter-class variance maximum, and its basic thought is that image pixel is divided into background and target two classes, calculates the inter-class variance maximal value by search, obtains optimal threshold, and the specific implementation step is as follows:
(4.1) number of pixel in interior each gray level of statistics region of interest ROI;
(4.2) number of pixels that calculates each gray level accounts for ratio and the total gray average of this region of interest ROI of this region of interest ROI;
(4.3) the traversal gray level is calculated and is determined maximum inter-class variance σ, and the gray-scale value of corresponding this maximum between-cluster variance is threshold value T:
σ
2=p
A(ω
A-ω
0)
2+p
B(ω
B-ω
0)
2 (8)
In the formula, p
A, ω
ABe respectively probability and gray average that background classes A occurs, p
B, ω
BBe respectively probability and gray average that target class B occurs, ω
0Be the total gray average of this region of interest ROI;
(4.4) according to threshold value T with image binaryzation, obtain the hot spot of corresponding LED luminotron 6 in this region of interest ROI.
The 5th step: according to the centroid position of each hot spot, obtain horizontal shift Δ x and the vertical displacement delta y of corresponding LED luminotron 6 in the image; Barycenter horizontal level M
x, upright position M
yComputing formula be:
In the formula, (j is to be positioned at that (height represents the pixels tall of region of interest ROI for j, i) gray values of pixel points, and width represents the pixel wide of region of interest ROI i) to f.
Determine horizontal shift Δ x, the vertical displacement delta y of flexible cantilever beam 5 upper end LED luminotrons 6 in the image:
Δx=M
x-x
0,Δy=M
y-y
0 (11)
Wherein, x
0The horizontal coordinate of representing flexible cantilever beam 5 upper end LED luminotrons 6 when static, y
0The vertical coordinate of representing flexible cantilever beam 5 upper end LED luminotrons 6 when static.Horizontal shift Δ x, vertical displacement delta y represent on level after 6 imagings of LED luminotron, the vertical direction the shared pixel size of displacement in image respectively.
The 6th step: according to horizontal shift Δ x, the vertical displacement delta y of each LED luminotron 6 in the image and the horizontal transformation factor K of the 1st step demarcation
x, vertical transformation factor K
y, obtain the real standard displacement x of each LED luminotron 6, actual perpendicular displacement y:
x=K
xΔx,y=K
yΔy (12)
The 7th step: according to the acquisition order of the 2nd step vibrational image sequence, draw out the position curve of each LED luminotron 6, Fast Fourier Transform (FFT) is carried out in the displacement of each LED luminotron 6, obtain the single order mode frequencies omega of flexible cantilever beam 5 vibrations
1With the second-order modal frequencies omega
2
Embodiment 1
Flexible cantilever girder construction with the analog solar windsurfing is example, and measurement mechanism and the method for flexible cantilever beam vibration that the present invention is based on machine vision is specific as follows.
In conjunction with Fig. 1, the present invention is based on the measurement mechanism of the flexible cantilever beam vibration of machine vision, comprise fixed support 1, flexible cantilever beam 5, camera support 2, CCD camera 3, camera lens 4,6 LED luminotrons 6 and PC 7; Described flexible cantilever beam 5 one ends are fixed on the fixed support 1, and the other end is free end, and the dimensioning modest ability of flexible cantilever beam 5 * wide * height is 1200mm * 80mm * 2mm; LED luminotron 6 is along be respectively a LED luminotron 6-1, the 2nd LED luminotron 6-2, the 3rd LED luminotron 6-3, the 4th LED luminotron 6-4, the 5th LED luminotron 6-5 and the 6th LED luminotron 6-6 from the fixed support end to free-ended direction, be successively set on the upper surface of flexible cantilever beam 5, equidistantly paste 6 red LED luminotrons along the every 200mm of the center line of flexible cantilever beam upper surface.Fixed support 1 top is provided with camera support 2, and CCD camera 3 is fixed on the camera support 2, and the model of CCD camera 3 is Cognex In-Sight5400, sensor type is 1/3 inch CCD, resolution is 640 * 480, dark 256 grey levels in position, per second 60 two field pictures; CCD camera 3 disposes camera lens 4, and the C mouth mirror head 4 of 8mm focal length is installed on the CCD camera; Under this configuration, the visual angle of camera is 38.5 °.The output port of CCD camera 3 is connected with PC 7 by data line.Utilize CCD camera 3 to detect the low-frequency vibration image of flexible cantilever beam 5 as sensor.The CCD camera that adopts in the present embodiment, d
x=d
ySo pixel is identical with the physical size of y axle along x in the image.
Fig. 2 is the upper surface vertical view of flexible cantilever beam 5, comprises that flexible cantilever beam upper surface 13, flexible cantilever beam upper surface center line 14, CCD camera optical axis 15, optical axis are at the projection line 16 of property semi-girder upper surface.The position of CCD camera 3 relative flexibility semi-girders 5 stiff ends and attitude are determined, all LED luminotrons 6 are in the field range of CCD camera 3 all the time when making flexible cantilever beam 5 vibrations, the CCD plane of CCD camera 3 and the upper surface of flexible cantilever beam 5 13 angled α, α=60 °, as shown in Figure 2, flexible cantilever beam upper surface 13 and CCD camera optical axis ° angle, 15 one-tenth 90 °-α=30 when static, and CCD camera optical axis and last plane center line intersect, CCD camera optical axis at the projection line on plane on the flexible cantilever beam on the center line of last plane.
Based on the measuring method of the flexible cantilever beam vibration of machine vision, step is as follows:
The 1st step: according to the parameter of CCD camera 3,6 LED luminotrons 6 are tied to the horizontal transformation factor K of world coordinate system from image coordinate
x, vertical transformation factor K
yDemarcate, because in the actual imaging, the K of diverse location place on the flexible cantilever beam 5
x, K
yValue changes, so 66 of LED luminotrons are located all will demarcate, obtains and need demarcate 6 LED luminotron places, obtains (K
X1, K
Y1), (K
X2, K
Y2), (K
X3, K
Y3), (K
X4, K
Y4), (K
X5, K
Y5), (K
X6, K
Y6), because d
x=d
y, then formula (7) changes into:
In the image on the flexible cantilever beam plane midline a bit at image coordinate system O-UV and world coordinate system O
w-X
wY
wZ
wIn initial position be respectively (u
1, v
1) and (x
W1, y
W1, 0), the position in a certain moment is (u during vibration
2, v
2) and (x
W2, y
W2, 0), this pass between the displacement under image coordinate system and the world coordinate system is:
The 2nd step: CCD camera 3 is gathered the vibrational image sequence of flexible cantilever beams 5, and exports the vibrational image sequence to PC 12 by data line; By the PC program to the vibrational image sequence identification of being correlated with handle the displacement of obtaining LED luminotron 6 places on the flexible cantilever beam 5.
The 3rd step: according to the position of each LED luminotron 6 in first two field picture, entire image is divided into 6 sub regions corresponding with each LED luminotron 6, and determines the region of interest ROI of each subregion; According to the maximum perpendicular displacement of LED luminotron 6 each subregion is chosen the region of interest ROI of image, reduce the data volume of handling, the LED luminotron mainly is horizontal displacement x, and y is less for the perpendicular displacement amount.The ROI that the whole horizontal directions of each LED luminotron sub-region processes, vertical 50 pixel datas are formed, namely the ROI size is 640 * 50.
The 4th step: adopt the OTSU thresholding method to extract the hot spot of the interior corresponding LED luminotron 6 of each region of interest ROI of flexible cantilever beam 5, and determine the centroid position of each hot spot.
The 5th step: according to the centroid position of each hot spot, obtain horizontal shift Δ x and the vertical displacement delta y of corresponding LED luminotron 6 in the image, acquisition order according to the 2nd step vibrational image sequence, the horizontal shift curve of the 3rd LED luminotron 6-3 as shown in Figure 4, horizontal ordinate is that the time, (unit: second), ordinate was displacement (unit: pixel).
The 6th step: according to horizontal shift Δ x, the vertical displacement delta y of each LED luminotron 6 in the image and the horizontal transformation factor K of the 1st step demarcation
x, vertical transformation factor K
y, obtain the real standard displacement x of each LED luminotron 6, actual perpendicular displacement y.
The 7th step: according to the acquisition order of the 2nd step vibrational image sequence, draw out the position curve of each LED luminotron 6, Fast Fourier Transform (FFT) is carried out in the displacement of each LED luminotron 6, obtain the single order mode frequencies omega of flexible cantilever beam 5 vibrations
1With the second-order modal frequencies omega
2The horizontal shift of the 3rd LED luminotron 6-3 is carried out curve that Fast Fourier Transform (FFT) obtains as shown in Figure 5, and horizontal ordinate is frequency (unit: hertz), and ordinate is amplitude.
In sum, the constructed machine vision vibration measurement device of method that the present invention proposes can be followed the tracks of the method for light source by machine vision under the situation of the vibration characteristics that does not change measured object, detect the multiple spot dynamic displacement of flexible structure, and obtain the low-frequency vibration modal information.This method has the plurality of advantages such as vibration characteristics that noncontact, measurement range are wide, do not change measured object.
Claims (6)
1. measurement mechanism based on the flexible cantilever beam vibration of machine vision, comprise fixed support (1), flexible cantilever beam (5), it is characterized in that this measurement mechanism also comprises camera support (2), CCD camera (3), camera lens (4), a plurality of LED luminotron (6) and PC (7); Described flexible cantilever beam (5) one ends are fixed on the fixed support (1), and the other end is free end, and each LED luminotron (6) is along the upper surface that is successively set on flexible cantilever beam (5) from the fixed support end to free-ended direction; Fixed support (1) top is provided with camera support (2), and CCD camera (3) is fixed on the camera support (2), and CCD camera (3) disposes camera lens (4), and the output port of CCD camera (3) is connected with PC (7).
2. the measurement mechanism of the flexible cantilever beam vibration based on machine vision according to claim 1, it is characterized in that: the quantity of described LED luminotron (6) is 6, be respectively a LED luminotron (6-1), the 2nd LED luminotron (6-2), the 3rd LED luminotron (6-3), the 4th LED luminotron (6-4), the 5th LED luminotron (6-5) and the 6th LED luminotron (6-6), and each LED luminotron is fixed on equidistantly on the center line of flexible cantilever beam (5) upper surface.
3. the measurement mechanism of the flexible cantilever beam vibration based on machine vision according to claim 1, it is characterized in that: described CCD camera (3) setting angle meets the following conditions: each LED luminotron (6) is all in the visual field of CCD camera (3).
4. measuring method based on the flexible cantilever beam vibration of machine vision is characterized in that step is as follows:
The 1st step: according to the parameter of CCD camera (3), each LED luminotron (6) is tied to the horizontal transformation factor K of world coordinate system from image coordinate
x, vertical transformation factor K
yDemarcate;
The 2nd step: CCD camera (3) is gathered the vibrational image sequence of flexible cantilever beam (5), and exports the vibrational image sequence to PC (12) by data line;
The 3rd step: according to the position of each LED luminotron (6) in first two field picture, entire image is divided into the subregion corresponding with each LED luminotron (6), and determines the region of interest ROI of each subregion;
The 4th step: adopt the OTSU thresholding method to extract the hot spot of the interior corresponding LED luminotron (6) of each region of interest ROI of flexible cantilever beam (5), and determine the centroid position of each hot spot;
The 5th step: according to the centroid position of each hot spot, obtain horizontal shift Δ x and the vertical displacement delta y of corresponding LED luminotron (6) in the image;
The 6th step: according to horizontal shift Δ x, the vertical displacement delta y of each LED luminotron (6) in the image and the horizontal transformation factor K of the 1st step demarcation
x, vertical transformation factor K
y, obtain the real standard displacement x of each LED luminotron (6), actual perpendicular displacement y;
The 7th step: according to the acquisition order of the 2nd step vibrational image sequence, draw out the position curve of each LED luminotron (6), Fast Fourier Transform (FFT) is carried out in displacement to each LED luminotron (6), obtains single order model frequency and the second-order modal frequency of flexible cantilever beam (5) vibration.
5. the measuring method of the flexible cantilever beam vibration based on machine vision according to claim 4 is characterized in that, the described horizontal transformation factor K that each LED luminotron (6) is tied to world coordinate system from image coordinate of the 1st step
x, vertical transformation factor K
yDemarcate, concrete steps are as follows:
(1.1) set up coordinate system: with the mid point O of flexible cantilever beam (5) upper surface fixed edge
wBe initial point, set up world coordinate system O
w-X
wY
wZ
wPhotocentre O with CCD camera (3)
CBe initial point, setting up camera coordinates is O
c-X
cY
cZ
cBe initial point O with first pixel of the image upper left corner, set up image coordinate system O-UV; Central point O with image
dBe initial point, setting up the image planes physical coordinates is O
d-X
dY
d, O wherein
dCoordinate be (c
x, c
y); Optical axis and the X of CCD camera (3)
cO
cY
c, UOV, X
dO
dY
dThe plane is all vertical, the angled α of upper surface of the CCD imaging plane of CCD camera (3) and flexible cantilever beam (5), 0 °<α<90 °, optical axis and the X of CCD camera (3)
WO
WY
WThe angled 90 °-α in plane, wherein X
WAxle, X
CAxle, X
dAxle is parallel;
(1.2) world coordinate system O
w-X
wY
wZ
wBe O to camera coordinates
c-X
cY
cZ
cTransformation relation be:
(x wherein
c, y
c, z
c) for camera coordinates be O
c-X
cY
cZ
cCoordinate, (x
w, y
w, z
w) be world coordinate system O
w-X
wY
wZ
wCoordinate,
In the formula, R is 3 * 3 rotational transform matrix, and θ is the anglec of rotation around X-axis; T is 3 * 1 translation vector, and b is shift value;
(1.3) camera coordinates is O
c-X
cY
cZ
cBe O to the image planes physical coordinates
d-X
dY
dCoordinate transform be:
In the formula, (x
d, y
d) for the image planes physical coordinates be O
d-X
dY
dCoordinate, f is that effective focal length is the distance of the plane of delineation and optical center of lens;
(1.4) the image planes physical coordinates is O
d-X
dY
dTo being transformed to of image coordinate system O-UV:
In the formula, (u v) is the coordinate of image coordinate system O-UV, d
xBe in the image pixel along the physical size of x axle, d
yBe in the image pixel along the physical size of y axle;
(1.5) composite type (1)~(3) are by world coordinate system O
w-X
wY
wZ
wHomogeneous coordinates to image coordinate system O-UV conversion are expressed as:
In the formula, S is scale-up factor, and formula (4) is transformed to:
Because LED luminotron (6) is only on the plane
Internal vibration, so following formula (5) is deformed into:
Get image coordinate system O-UV to world coordinate system O by formula (6)
w-X
wY
wZ
wThe horizontal transformation factor K
xWith vertical transformation factor K
y, be shown below respectively:
6. the measuring method of the flexible cantilever beam vibration based on machine vision according to claim 4 is characterized in that, described OTSU thresholding method of the 4th step, and step is as follows:
(4.1) number of pixel in interior each gray level of statistics region of interest ROI;
(4.2) number of pixels that calculates each gray level accounts for ratio and the total gray average of this region of interest ROI of this region of interest ROI;
(4.3) the traversal gray level is calculated and is determined maximum inter-class variance, and the gray-scale value of corresponding this maximum between-cluster variance is threshold value;
(4.4) according to threshold value with image binaryzation, obtain the hot spot of corresponding LED luminotron (6) in this region of interest ROI.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310224821.9A CN103267567B (en) | 2013-06-05 | 2013-06-05 | Based on measurement mechanism and the method for the flexible cantilever beam vibration of machine vision |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310224821.9A CN103267567B (en) | 2013-06-05 | 2013-06-05 | Based on measurement mechanism and the method for the flexible cantilever beam vibration of machine vision |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103267567A true CN103267567A (en) | 2013-08-28 |
CN103267567B CN103267567B (en) | 2015-09-30 |
Family
ID=49011207
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310224821.9A Active CN103267567B (en) | 2013-06-05 | 2013-06-05 | Based on measurement mechanism and the method for the flexible cantilever beam vibration of machine vision |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103267567B (en) |
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103954221A (en) * | 2014-05-08 | 2014-07-30 | 哈尔滨工业大学 | Binocular photogrammetry method of large flexible structure vibration displacement |
CN104048744A (en) * | 2014-07-08 | 2014-09-17 | 安徽常春藤光电智能科技有限公司 | Non-contact real-time online vibration measurement method based on images |
CN104089697A (en) * | 2014-07-08 | 2014-10-08 | 安徽常春藤光电智能科技有限公司 | Real-time online visual vibration measurement method based on thread pool concurrent technology |
CN105157938A (en) * | 2015-05-04 | 2015-12-16 | 南京理工大学 | Device and method for adjusting modal frequency of flexible cantilever |
CN105953906A (en) * | 2016-04-25 | 2016-09-21 | 上海航天控制技术研究所 | Ultra large flexible structure perfect information measuring system and method thereof |
CN106124034A (en) * | 2016-09-07 | 2016-11-16 | 湖南科技大学 | Thin-wall part operation mode based on machine vision test device and method of testing |
CN106153331A (en) * | 2016-09-18 | 2016-11-23 | 安徽理工大学 | The measurement apparatus vibrated based on the flexible doublejointed mechanical arm of machine vision and method |
CN106644034A (en) * | 2016-12-08 | 2017-05-10 | 南京理工大学 | Non-contact high-speed rail road foundation vibration measurement system |
CN106885622A (en) * | 2017-02-07 | 2017-06-23 | 上海理工大学 | A kind of big visual field multiple spot three-dimensional vibrating measuring method |
CN107271026A (en) * | 2017-07-07 | 2017-10-20 | 河南科技大学 | A kind of method of steel wire rope oscillation crosswise measurement |
CN107817044A (en) * | 2017-11-30 | 2018-03-20 | 安徽理工大学 | The measurement apparatus and method of sheet material vibration based on machine vision |
CN108593087A (en) * | 2018-03-29 | 2018-09-28 | 湖南科技大学 | A kind of thin-wall part operational modal parameter determines method and system |
CN108731788A (en) * | 2018-05-22 | 2018-11-02 | 河海大学常州校区 | A kind of working at height arm low-frequency vibration vision inspection apparatus and method |
CN108760022A (en) * | 2018-06-20 | 2018-11-06 | 航天海鹰(镇江)特种材料有限公司 | A kind of digitized measurement of wing class product vibration frequency and the method for Dynamic Model |
CN108952682A (en) * | 2018-08-15 | 2018-12-07 | 中国地质大学(武汉) | A kind of wellbore annulus pressure sensor based on machine vision |
CN109727267A (en) * | 2018-12-29 | 2019-05-07 | 中国计量科学研究院 | A kind of virtual sinusoidal straight-line oscillation measurement method of standard |
CN109883533A (en) * | 2019-01-21 | 2019-06-14 | 哈尔滨工业大学(深圳) | Low frequency vibration measurement method based on machine vision |
CN110849314A (en) * | 2018-12-26 | 2020-02-28 | 中国航空工业集团公司北京长城计量测试技术研究所 | Dynamic strain calibration method with continuously adjustable frequency |
CN111366094A (en) * | 2020-03-05 | 2020-07-03 | 山东大学 | Method for detecting and calibrating large deformation space pose of flexible body by fiber bragg grating |
CN112444365A (en) * | 2020-11-30 | 2021-03-05 | 哈尔滨工业大学 | Satellite solar wing substrate unfolding low-frequency modal testing method based on force hammer swing method and laser Doppler method |
CN114518161A (en) * | 2022-01-11 | 2022-05-20 | 武汉理工大学 | Measuring method for transverse vibration of end part of cylindrical simply supported beam based on computer vision |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101387546A (en) * | 2008-10-29 | 2009-03-18 | 华南理工大学 | Space sailboard bend and torsion low frequency mode detecting method and apparatus based on vision |
CN101650179A (en) * | 2009-09-14 | 2010-02-17 | 中南大学 | Detection method and system of offset of contact net |
CN102322938A (en) * | 2011-08-22 | 2012-01-18 | 华南理工大学 | Visual measuring device for vibration of piezoelectric flexible board structure and control method |
-
2013
- 2013-06-05 CN CN201310224821.9A patent/CN103267567B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101387546A (en) * | 2008-10-29 | 2009-03-18 | 华南理工大学 | Space sailboard bend and torsion low frequency mode detecting method and apparatus based on vision |
CN101650179A (en) * | 2009-09-14 | 2010-02-17 | 中南大学 | Detection method and system of offset of contact net |
CN102322938A (en) * | 2011-08-22 | 2012-01-18 | 华南理工大学 | Visual measuring device for vibration of piezoelectric flexible board structure and control method |
Cited By (35)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103954221B (en) * | 2014-05-08 | 2016-08-17 | 哈尔滨工业大学 | The binocular photogrammetric survey method of large-size pliable structure vibration displacement |
CN103954221A (en) * | 2014-05-08 | 2014-07-30 | 哈尔滨工业大学 | Binocular photogrammetry method of large flexible structure vibration displacement |
CN104048744A (en) * | 2014-07-08 | 2014-09-17 | 安徽常春藤光电智能科技有限公司 | Non-contact real-time online vibration measurement method based on images |
CN104089697A (en) * | 2014-07-08 | 2014-10-08 | 安徽常春藤光电智能科技有限公司 | Real-time online visual vibration measurement method based on thread pool concurrent technology |
CN104089697B (en) * | 2014-07-08 | 2017-02-15 | 安徽常春藤光电智能科技有限公司 | Real-time online visual vibration measurement method based on thread pool concurrent technology |
CN105157938A (en) * | 2015-05-04 | 2015-12-16 | 南京理工大学 | Device and method for adjusting modal frequency of flexible cantilever |
CN105157938B (en) * | 2015-05-04 | 2018-12-14 | 南京理工大学 | The device and method for adjusting flexible cantilever beam modal frequency |
CN105953906B (en) * | 2016-04-25 | 2018-08-10 | 上海航天控制技术研究所 | A kind of ultra-large type flexible structure full information measuring system and method |
CN105953906A (en) * | 2016-04-25 | 2016-09-21 | 上海航天控制技术研究所 | Ultra large flexible structure perfect information measuring system and method thereof |
CN106124034A (en) * | 2016-09-07 | 2016-11-16 | 湖南科技大学 | Thin-wall part operation mode based on machine vision test device and method of testing |
CN106153331A (en) * | 2016-09-18 | 2016-11-23 | 安徽理工大学 | The measurement apparatus vibrated based on the flexible doublejointed mechanical arm of machine vision and method |
CN106644034A (en) * | 2016-12-08 | 2017-05-10 | 南京理工大学 | Non-contact high-speed rail road foundation vibration measurement system |
CN106885622A (en) * | 2017-02-07 | 2017-06-23 | 上海理工大学 | A kind of big visual field multiple spot three-dimensional vibrating measuring method |
CN106885622B (en) * | 2017-02-07 | 2019-08-30 | 上海理工大学 | A kind of big visual field multiple spot three-dimensional vibrating measurement method |
CN107271026A (en) * | 2017-07-07 | 2017-10-20 | 河南科技大学 | A kind of method of steel wire rope oscillation crosswise measurement |
CN107271026B (en) * | 2017-07-07 | 2020-01-03 | 河南科技大学 | Method for measuring transverse vibration of steel wire rope |
CN107817044A (en) * | 2017-11-30 | 2018-03-20 | 安徽理工大学 | The measurement apparatus and method of sheet material vibration based on machine vision |
CN107817044B (en) * | 2017-11-30 | 2023-08-22 | 安徽理工大学 | Device and method for measuring plate vibration based on machine vision |
CN108593087A (en) * | 2018-03-29 | 2018-09-28 | 湖南科技大学 | A kind of thin-wall part operational modal parameter determines method and system |
CN108731788A (en) * | 2018-05-22 | 2018-11-02 | 河海大学常州校区 | A kind of working at height arm low-frequency vibration vision inspection apparatus and method |
CN108760022A (en) * | 2018-06-20 | 2018-11-06 | 航天海鹰(镇江)特种材料有限公司 | A kind of digitized measurement of wing class product vibration frequency and the method for Dynamic Model |
CN108952682B (en) * | 2018-08-15 | 2020-08-07 | 中国地质大学(武汉) | Shaft annulus pressure sensor based on machine vision |
CN108952682A (en) * | 2018-08-15 | 2018-12-07 | 中国地质大学(武汉) | A kind of wellbore annulus pressure sensor based on machine vision |
CN110849314B (en) * | 2018-12-26 | 2021-06-22 | 中国航空工业集团公司北京长城计量测试技术研究所 | Dynamic strain calibration method with continuously adjustable frequency |
CN110849314A (en) * | 2018-12-26 | 2020-02-28 | 中国航空工业集团公司北京长城计量测试技术研究所 | Dynamic strain calibration method with continuously adjustable frequency |
CN109727267B (en) * | 2018-12-29 | 2020-02-11 | 中国计量科学研究院 | Standard virtual sine linear vibration measurement method |
CN109727267A (en) * | 2018-12-29 | 2019-05-07 | 中国计量科学研究院 | A kind of virtual sinusoidal straight-line oscillation measurement method of standard |
CN109883533B (en) * | 2019-01-21 | 2021-05-25 | 哈尔滨工业大学(深圳) | Low-frequency vibration measurement method based on machine vision |
CN109883533A (en) * | 2019-01-21 | 2019-06-14 | 哈尔滨工业大学(深圳) | Low frequency vibration measurement method based on machine vision |
CN111366094A (en) * | 2020-03-05 | 2020-07-03 | 山东大学 | Method for detecting and calibrating large deformation space pose of flexible body by fiber bragg grating |
CN111366094B (en) * | 2020-03-05 | 2021-02-12 | 山东大学 | Method for detecting and calibrating large deformation space pose of flexible body by fiber bragg grating |
CN112444365A (en) * | 2020-11-30 | 2021-03-05 | 哈尔滨工业大学 | Satellite solar wing substrate unfolding low-frequency modal testing method based on force hammer swing method and laser Doppler method |
CN112444365B (en) * | 2020-11-30 | 2023-08-29 | 哈尔滨工业大学 | Satellite solar wing substrate unfolding low-frequency mode testing method |
CN114518161A (en) * | 2022-01-11 | 2022-05-20 | 武汉理工大学 | Measuring method for transverse vibration of end part of cylindrical simply supported beam based on computer vision |
CN114518161B (en) * | 2022-01-11 | 2023-02-24 | 武汉理工大学 | Measuring method for transverse vibration of end part of cylindrical simply supported beam based on computer vision |
Also Published As
Publication number | Publication date |
---|---|
CN103267567B (en) | 2015-09-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103267567B (en) | Based on measurement mechanism and the method for the flexible cantilever beam vibration of machine vision | |
Sony et al. | A literature review of next‐generation smart sensing technology in structural health monitoring | |
CN101387546B (en) | Space sailboard bend and torsion low frequency mode detecting method and apparatus based on vision | |
CN104729665A (en) | Flexible cantilever slab vibration detection device and method based on binocular vision | |
Bartilson et al. | Target-less computer vision for traffic signal structure vibration studies | |
CN102589523A (en) | Method and equipments for remotely monitoring displacement of building | |
CN102289306A (en) | Attitude sensing equipment and positioning method thereof as well as method and device for controlling mouse pointer | |
CN102506727B (en) | Method and system for measuring vibration displacement of arm support and concrete pump truck comprising same | |
CN203178143U (en) | Automatic bending rigidity test device | |
CN107153382A (en) | Flexible hinged plate vibration control apparatus and method based on Binocular vision photogrammetry | |
CN103727876A (en) | Strip width and center measurement system and method based on parallel laser rays | |
CN102798456A (en) | Method, device and system for measuring working amplitude of engineering mechanical arm support system | |
CN106153331B (en) | The measuring device and method of flexible doublejointed mechanical arm vibration based on machine vision | |
CN106483330A (en) | One kind is based on reflective silk thread attitude angle visual identity two-D wind speed wind direction method of testing | |
CN102788572A (en) | Method, device and system for measuring attitude of engineering machinery lifting hook | |
CN103487011A (en) | Method for detecting attitude angle of data glove | |
CN106291567A (en) | A kind of laser radar, laser radar system and data processing method | |
Bai et al. | UAV based accurate displacement monitoring through automatic filtering out its camera's translations and rotations | |
CN104913775A (en) | Method for measuring height of transmission line of unmanned aerial vehicle and method and device for positioning unmanned aerial vehicle | |
Luo et al. | A novel static deformation measurement and visualization method for wind turbine blades using home-made LiDAR and processing program | |
CN102306054B (en) | Attitude sensing equipment and positioning method and device thereof and method and device for controlling mouse pointer | |
CN103090863B (en) | Method for measuring posture and height of dynamic platform | |
CN206892604U (en) | Flexible hinged plate vibration control apparatus based on Binocular vision photogrammetry | |
Wang et al. | Real-time vibration visualization using GPU-based high-speed vision | |
Lim et al. | Visual measurement of pile movements for the foundation work using a high-speed line-scan camera |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |