CN112304493B - CCD camera-based optical pressure-sensitive paint amplitude-frequency characteristic detection method - Google Patents
CCD camera-based optical pressure-sensitive paint amplitude-frequency characteristic detection method Download PDFInfo
- Publication number
- CN112304493B CN112304493B CN202011186794.7A CN202011186794A CN112304493B CN 112304493 B CN112304493 B CN 112304493B CN 202011186794 A CN202011186794 A CN 202011186794A CN 112304493 B CN112304493 B CN 112304493B
- Authority
- CN
- China
- Prior art keywords
- ccd camera
- pressure
- frequency
- standing wave
- amplitude
- 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
Links
- 230000003287 optical effect Effects 0.000 title claims abstract description 33
- 239000003973 paint Substances 0.000 title claims description 16
- 238000001514 detection method Methods 0.000 title description 9
- 239000011248 coating agent Substances 0.000 claims abstract description 38
- 238000000576 coating method Methods 0.000 claims abstract description 38
- 238000000034 method Methods 0.000 claims abstract description 33
- 238000002073 fluorescence micrograph Methods 0.000 claims description 60
- 230000008569 process Effects 0.000 claims description 11
- 230000008859 change Effects 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 5
- 238000004519 manufacturing process Methods 0.000 claims description 3
- 230000003247 decreasing effect Effects 0.000 claims description 2
- 238000009826 distribution Methods 0.000 abstract description 4
- 239000000523 sample Substances 0.000 description 32
- 238000005259 measurement Methods 0.000 description 6
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 5
- 230000008901 benefit Effects 0.000 description 5
- 229910052760 oxygen Inorganic materials 0.000 description 5
- 239000001301 oxygen Substances 0.000 description 5
- 230000005540 biological transmission Effects 0.000 description 3
- 230000005281 excited state Effects 0.000 description 3
- 238000013500 data storage Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000005283 ground state Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000000171 quenching effect Effects 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 238000009530 blood pressure measurement Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000004020 luminiscence type Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- QVRVXSZKCXFBTE-UHFFFAOYSA-N n-[4-(6,7-dimethoxy-3,4-dihydro-1h-isoquinolin-2-yl)butyl]-2-(2-fluoroethoxy)-5-methylbenzamide Chemical compound C1C=2C=C(OC)C(OC)=CC=2CCN1CCCCNC(=O)C1=CC(C)=CC=C1OCCF QVRVXSZKCXFBTE-UHFFFAOYSA-N 0.000 description 1
- 230000010349 pulsation Effects 0.000 description 1
- 238000010791 quenching Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L25/00—Testing or calibrating of apparatus for measuring force, torque, work, mechanical power, or mechanical efficiency
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/28—Investigating the spectrum
- G01J3/2803—Investigating the spectrum using photoelectric array detector
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/28—Investigating the spectrum
- G01J3/2823—Imaging spectrometer
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/28—Investigating the spectrum
- G01J3/44—Raman spectrometry; Scattering spectrometry ; Fluorescence spectrometry
- G01J3/4406—Fluorescence spectrometry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L1/00—Measuring force or stress, in general
- G01L1/24—Measuring force or stress, in general by measuring variations of optical properties of material when it is stressed, e.g. by photoelastic stress analysis using infrared, visible light, ultraviolet
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/62—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
- G01N21/63—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
- G01N21/64—Fluorescence; Phosphorescence
- G01N21/645—Specially adapted constructive features of fluorimeters
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
Abstract
The invention provides a CCD camera-based method for detecting the amplitude-frequency characteristic of an optical pressure-sensitive coating, which comprises the following steps: the signal generator is connected with an external trigger interface of the laser light source and controls the working state of the laser light source; meanwhile, the signal generator is connected with an external trigger interface of the CCD camera to control the working state of the CCD camera, and when the frequency of the sinusoidal pressure standing wave is f1While obtaining the amplitude Am(f1) (ii) a Keeping the amplitude of the sinusoidal pressure standing wave stable and unchanged, and controlling the frequency of the sinusoidal pressure standing wave from f1Gradually increasing to obtain the cut-off frequency and the limiting frequency of the PSP coating. According to the method, the time sequence of the CCD camera is controlled, and the low-frame-rate CCD camera collects the fluorescent image sequence of the pressure sensitive coating under the high-frequency pulsating pressure, so that the dynamic pressure distribution of the universe is obtained.
Description
Technical Field
The invention belongs to the technical field of pressure-sensitive paint amplitude-frequency characteristic detection, and particularly relates to an optical pressure-sensitive paint amplitude-frequency characteristic detection method based on a CCD camera.
Background
The pressure is one of three parameters (pressure, temperature and flow) of thermal engineering for automatic control, and plays an important role in measurement and control. There are many ways of measuring pressure, but each requires calibration of the measuring device in order to obtain measurement data, such as the functional relationship between electrical signals, optical signals, etc. and pressure, and the characteristics of the measuring system, such as sensitivity, etc.
The Pressure Sensitive Paint (PSP) technique based on computer vision and image processing technique is an important breakthrough of non-contact flow display technique. Compared with the traditional dot matrix measurement technology in China at present, the optical pressure-sensitive measurement technology can make up for the damage of pressure probe holes, pressure sensors and the like to a model and the interference to a flow field and the complexity of a data transmission mode of a traditional method, greatly improves the measurement range, has the advantages of no contact, continuous measurement, relatively low test cost, time saving and the like, and is favored by vast experiment workers. The basic principle of the optical pressure-sensitive pressure measurement technology is as follows: the pressure sensitive paint is uniformly covered on the surface of the tested model, and consists of photosensitive molecules and an oxygen permeable substrate. When excited by light of a specific wavelength, the photosensitive molecules in the paint transition from an originally stable ground state to an excited state of a high energy level. The photosensitive molecules in the unstable excited state are collided by oxygen molecules diffused from the measured surface, the energy of the excited state is lost, and the photosensitive molecules are inactivated and returned to the ground state, and the process does not generate radiant light, so that the luminous intensity is reduced, and an oxygen quenching phenomenon is formed. The greater the concentration of oxygen molecules, i.e.: the higher the pressure in the atmosphere, the stronger the quenching effect of oxygen, and the darker the coating will emit under a certain light. Therefore, under the irradiation of light, the luminous intensity of the pressure sensitive coating can reflect the pressure value on the surface of the measured model. And taking an image picture of the surface of the measured model under the light irradiation, and analyzing the image picture to obtain the pressure distribution of the surface of the measured model. Optical pressure sensitive manometry requires pre-calibration to obtain paint properties.
In the prior art, the detection technology of the amplitude-frequency characteristic of the optical pressure-sensitive paint mainly has the following problems:
in the detection process of the amplitude-frequency characteristic of the optical pressure-sensitive coating, in order to obtain the dynamic pressure distribution of the universe, the acquisition of images must be completed through a camera, and due to the fact that the frame rate of a CCD camera is low, when a pressure-sensitive coating fluorescence image sequence under high-frequency pulsating pressure is acquired, a large system error can be introduced, and therefore the detection accuracy of the amplitude-frequency characteristic of the optical pressure-sensitive coating is reduced.
Disclosure of Invention
Aiming at the defects in the prior art, the invention provides a CCD camera-based method for detecting the amplitude-frequency characteristic of an optical pressure-sensitive coating, which can effectively solve the problems.
The technical scheme adopted by the invention is as follows:
the invention provides a CCD camera-based method for detecting the amplitude-frequency characteristic of an optical pressure-sensitive coating, which comprises the following steps:
the CCD camera (5) and the laser light source (3) are aligned to the PSP sample wafer (2) through an optical window (7) of the acoustic standing wave tube (1), and an optical filter (6) is fixedly arranged in front of a lens of the CCD camera (5); the 2 nd channel of the signal generator (8) is connected with the CCD camera (5); the 3 rd channel of the signal generator (8) is connected with the laser light source (3);
fixing a dynamic pressure sensor (4) on the cross section where the PSP sample wafer (2) is located;
the output end of the CCD camera (5) and the output end of the dynamic pressure sensor (4) are both connected to a computer (10);
the camera and the light source are aligned to the PSP sample through the optical window, and the dynamic pressure sensor is fixed on the cross section of the PSP sample and used for measuring the real-time dynamic pressure and providing a trigger signal;
the continuous output period of the signal generator (8) is T1Frequency of f1Is transmitted to a power amplifier (9), the power amplifier (9) drives a sound source (11) to emit a sinusoidal signal with a period T1Frequency of f1The sinusoidal sound wave of (2);
the sine sound wave acts through the acoustic standing wave tube (1), so that the surface of the PSP sample wafer (2) is subjected to stable period T1Frequency of f1A sinusoidal pressure standing wave of (a);
step 3.1, for a period T1Frequency of f1The sinusoidal pressure standing wave of (1) equally dividing a sinusoidal pressure standing wave into n phases, which are sequentially expressed as phasesPhase positionPhase positionThe time intervals between adjacent phases are all
Presetting the exposure time t of each time of the CCD camerae(ii) a Wherein, te≤T1/(n-1);
Presetting the shooting time interval T of the CCD camera according to the frame rate of the CCD cameraC=m×T1+T1V (n-1); wherein m represents the number of camera shooting cycle intervals; shooting time interval T of cameraCGreater than the reciprocal of the maximum frame rate of the camera;
step 3.2, where t is tkAt the moment, the sine pressure standing wave is in the phase of the 1 st periodThe state of (1); the signal generator (8) controls the laser light source (3) to be turned on and the CCD camera (5) to be turned on simultaneously;
at the moment, after the laser light source (3) is turned on, the long-open mode is entered, namely: after the laser light source (3) is turned on, the stable power is kept to continuously emit light to irradiate the surface of the PSP sample wafer (2), so that the PSP coating is excited;
after the CCD camera (5) is opened, the 1 st exposure is carried out, and the exposure time is te(ii) a When the exposure time is reached, i.e.: when t is equal to tk+teWhen the CCD camera (5) is started, the CCD camera is closed; at this time, the CCD camera (5) outputs and phasesCorresponding fluorescence image Q1;
Step 3.3, start from the 1 st exposure start time of the camera, i.e. from t to tkThe time interval T of the CCD camera is counted from the beginningCAnd then, namely: t is tk+m×T1+T1When the sinusoidal pressure standing wave is in the m +1 th cycle phase (n-1)At this time, the CCD camera (5) is turned on, and the 2 nd exposure is performed for the exposure time te(ii) a When the exposure time is reached, i.e.: when t is equal to tk+m×T1+T1/(n-1)+teWhen the CCD camera (5) is started, the CCD camera is closed; at this time, the CCD camera (5) outputs and phasesCorresponding fluorescence image Q2;
Step 3.4, starting from the 2 nd exposure starting time of the camera, passing through the shooting time interval T of the CCD cameraCThen, the sinusoidal pressure standing wave is phase in the 2m +1 th periodIn the state (2), the CCD camera (5) is turned on to perform the 3 rd exposure, and the output and phase of the CCD camera (5) are adjustedCorresponding fluorescence image Q3;
And so on, in the process that the laser light source (3) keeps stable power and continuously emits light to irradiate the surface of the PSP sample wafer (2), the shooting time interval T of the CCD camera is set everyCThe CCD camera (5) is exposed once and outputs a fluorescence image corresponding to the current phase of the sine pressure standing wave; until the output of the CCD camera (5) and the last phaseCorresponding fluorescence image Qn;
To this end, when the frequency of the sinusoidal pressure standing wave is f1Then, a fluorescence image sequence corresponding to a complete cycle is obtained, which is respectively: and phaseCorresponding fluorescence image Q1And phaseCorresponding fluorescence image Q2…, and phaseCorresponding fluorescence image Qn;
step 4.1, to fluorescenceImage QjPerforming image processing to obtain sum phaseCorresponding light intensity Ij;
Step 4.2, obtaining a fluorescence image QjIn the acquisition process, the exposure time of the CCD camera (5) is searched, and the change curve between the real pressure value and the time on the surface of the PSP sample wafer (2) is searched according to the exposure time of the CCD camera (5) to obtain a fluorescence image QjCorresponding true pressure value Pj;
And 4.3, when the sinusoidal pressure standing wave is not applied to the PSP sample (2), namely: the PSP sample (2) is brought to atmospheric pressure, and the image of the PSP sample (2) is captured, so that a reference pressure P is obtainedrefAnd light intensity at a reference pressure Iref;
since j is 1, 2.., n, when j is 1, one equation for a and B is obtained; when j is 2, an equation for a and B is obtained; by analogy, when j is equal to n, one equation for a and B is obtained; thus, a total of n equations for A and B are obtained; solving n equations about A and B by adopting a least square method to obtain final values of A and B;
a and B are constants, and the values of A and B and reference pressure PrefAnd light intensity at a reference pressure IrefSubstituting the calibration equation (1) to obtain a calibration equation (2):
wherein, in the calibration equation (2), P'jIs equal to the light intensity IjA corresponding calibration pressure value;
continuously increasing the frequency of the sine pressure standing wave if the current frequency fzAmplitude A ofm(fz) When the value is reduced to 0, thenFront frequency fzIs the limiting frequency of the PSP coating; therefore, the amplitude-frequency characteristic of the PSP coating is detected.
Preferably, step 4.1 is specifically:
fluorescence image QjW pixel points are provided; each pixel point corresponds to a light intensity value, therefore, w light intensity values are obtained in total, and the average value of the w light intensity values is the phase positionCorresponding light intensity Ij。
The CCD camera-based method for detecting the amplitude-frequency characteristic of the optical pressure-sensitive paint provided by the invention has the following advantages:
according to the method, the time sequence of the CCD camera is controlled, and the low-frame-rate CCD camera collects the fluorescent image sequence of the pressure sensitive coating under the high-frequency pulsating pressure, so that the dynamic pressure distribution of the universe is obtained.
Drawings
FIG. 1 is a schematic structural diagram of a CCD camera-based optical pressure-sensitive paint amplitude-frequency characteristic detection system provided by the invention;
fig. 2 is a schematic control timing diagram of the CCD camera according to the present invention.
Detailed Description
In order to make the technical problems, technical solutions and advantageous effects solved by the present invention more clearly apparent, the present invention is further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
The invention acquires fluorescent images emitted by the optical pressure-sensitive coating under the sinusoidal pressure fluctuation with the same amplitude and different frequencies by a CCD camera with a low frame rate based on a phase locking method to acquire the light intensity pulsation, and obtains the amplitude-frequency characteristic of the coating according to the measured amplitude after the light intensity of the coating under different frequencies is converted into pressure.
The invention provides a CCD camera-based method for detecting the amplitude-frequency characteristic of an optical pressure-sensitive coating, which comprises the following steps:
the CCD camera 5 and the laser light source 3 are aligned to the PSP sample 2 through an optical window 7 of the acoustic standing wave tube 1, and an optical filter 6 is fixedly installed in front of a lens of the CCD camera 5; the 2 nd channel of the signal generator 8 is connected with the CCD camera 5; the 3 rd channel of the signal generator 8 is connected with the laser light source 3;
fixing a dynamic pressure sensor 4 on the cross section of the PSP sample 2;
the output end of the CCD camera 5 and the output end of the dynamic pressure sensor 4 are both connected to a computer 10;
the camera and the light source are aligned to the PSP sample through the optical window, and the dynamic pressure sensor is fixed on the cross section of the PSP sample and used for measuring the real-time dynamic pressure and providing a trigger signal;
the continuous output period of the signal generator 8 is T1Frequency of f1Is transmitted to the power amplifier 9, the power amplifier 9 drives the sound source 11 to emit sound with a period T1Frequency of f1The sinusoidal sound wave of (2);
the sine sound wave acts through the acoustic standing wave tube 1, so that the stable period of the surface of the PSP sample wafer 2 is T1Frequency of f1A sinusoidal pressure standing wave of (a);
step 3.1, for a period T1Frequency of f1The sinusoidal pressure standing wave of (1) equally dividing a sinusoidal pressure standing wave into n phases, which are sequentially expressed as phasesPhase positionPhase positionThe time intervals between adjacent phases are all
Presetting the exposure time t of each time of the CCD camerae(ii) a Wherein, te≤T1/(n-1);
Presetting the shooting time interval T of the CCD camera according to the frame rate of the CCD cameraC=m×T1+T1V (n-1); wherein m represents the number of camera shooting cycle intervals; shooting time interval T of cameraCGreater than the reciprocal of the maximum frame rate of the camera;
step 3.2, where t is tkAt the moment, the sine pressure standing wave is in the phase of the 1 st periodThe state of (1); the signal generator 8 controls the laser light source 3 and the CCD camera 5 to be turned on simultaneously;
at this time, after the laser light source 3 is turned on, the long on mode is entered, that is: after the laser light source 3 is turned on, the stable power is kept to continuously emit light to irradiate the surface of the PSP sample 2, so that the PSP coating is excited;
after the CCD camera 5 is opened, the 1 st exposure is carried out for the exposure timeIs te(ii) a When the exposure time is reached, i.e.: when t is equal to tk+teWhen so, the CCD camera 5 is turned off; at this time, the CCD camera 5 outputs and phasesCorresponding fluorescence image Q1;
Step 3.3, start from the 1 st exposure start time of the camera, i.e. from t to tkThe time interval T of the CCD camera is counted from the beginningCAnd then, namely: t is tk+m×T1+T1When the sinusoidal pressure standing wave is in the m +1 th cycle phase (n-1)At this time, the CCD camera 5 is turned on, and the 2 nd exposure is performed for the exposure time te(ii) a When the exposure time is reached, i.e.: when t is equal to tk+m×T1+T1/(n-1)+teWhen so, the CCD camera 5 is turned off; at this time, the CCD camera 5 outputs and phasesCorresponding fluorescence image Q2;
Step 3.4, starting from the 2 nd exposure starting time of the camera, passing through the shooting time interval T of the CCD cameraCThen, the sinusoidal pressure standing wave is phase in the 2m +1 th periodIn the state (2), the CCD camera 5 is turned on to perform the 3 rd exposure, and the output and phase of the CCD camera 5 are adjustedCorresponding fluorescence image Q3;
And so on, in the process that the laser light source 3 continuously emits light to irradiate the surface of the PSP sample 2 with stable power, the shooting time interval T of the CCD camera is setCThe CCD camera 5 is exposed once and outputs a fluorescence image corresponding to the current phase of the sine pressure standing wave; until the output of the CCD camera 5 reaches the maximumThe latter phaseCorresponding fluorescence image Qn;
To this end, when the frequency of the sinusoidal pressure standing wave is f1Then, a fluorescence image sequence corresponding to a complete cycle is obtained, which is respectively: and phaseCorresponding fluorescence image Q1And phaseCorresponding fluorescence image Q2…, and phaseCorresponding fluorescence image Qn;
For ease of understanding step 3.1-step 3.4, an example is set forth below in conjunction with FIG. 2:
assuming a period T of sinusoidal pressure standing waves1Is 1s, wherein s is unit second; dividing a sinusoidal pressure standing wave into n-5 phases on average, wherein the n-5 phases are respectively as follows: phase positionPhase positionPhase positionPhase positionPhase positionThe time interval between adjacent phases isOf course, in practical applications, a sinusoidal pressure standing wave may be divided into 10 phases, 20 phases, etc. on average according to practical requirements, which is not limited by the present invention.
Presetting the exposure time t of each time of the CCD camerae=0.1s;
Presetting shooting time interval T of CCD cameraC=m×T1+T1(n-1), assuming that m is 3, then TC=3.25s;
First, when t is equal to tkAt the moment, the sinusoidal pressure standing wave is in the phase of the 1 st cycleThe 1 st period at this time is a relative period, and at this time, the signal generator 8 controls the laser light source 3 to be turned on and the CCD camera 5 to be turned on at the same time;
at this time, after the laser light source 3 is turned on, the long on mode is entered, that is: after the laser light source 3 is turned on, the stable power is kept to continuously emit light to irradiate the surface of the PSP sample 2, so that the PSP coating is excited;
after the CCD camera 5 is opened, the 1 st exposure is carried out for te0.1 s; when the exposure time is reached, i.e.: when t is equal to tkWhen +0.1s, the CCD camera 5 is closed; at this time, the CCD camera 5 outputs and phasesCorresponding fluorescence image Q1;
Then, through TCAfter 3.25s, the sinusoidal pressure standing wave is phase position in the 4 th cycleIn the state (2), the CCD camera 5 is turned on, the 2 nd exposure is performed, and the output and phase of the CCD camera 5 are determinedCorresponding fluorescence image Q2;
And so on every TCThe CCD camera 5 performs one exposure for 3.25s, fromSequentially output and phaseCorresponding fluorescence image Q3And phaseCorresponding fluorescence image Q4And phaseCorresponding fluorescence image Q5。
It follows that, if the conventional method is adopted, one period T is used1When the sinusoidal pressure standing wave of 1s is divided into n-5 phases on average, the time interval between adjacent phases isThus, when a sinusoidal pressure standing wave is reached, ideally, the camera needs to take 5 fluorescence images at a frame rate of 0.25s, namely: the camera completes one exposure, output and phase within 0-0.25 sCorresponding fluorescence image Q1(ii) a Then, completing one exposure within 0.25 s-0.5 s, outputting and phaseCorresponding fluorescence image Q2(ii) a Repeating the steps until one exposure is finished within 1-1.25 s, and outputting and phase positionsCorresponding fluorescence image Q5. In general, the frame rate of the camera is difficult to reach 0.25s, and therefore, a large system error is introduced.
In the present invention, the camera completes one exposure, output and phase within 0-0.25 sCorresponding fluorescence image Q1(ii) a Then, through TCThe camera completes one exposure within 3.25 s-3.5 s at a time interval of 3.25s, and outputs and phasesCorresponding fluorescence image Q2(ii) a And so on every TCThe camera performs one exposure at a time interval of 3.25s, and outputs a fluorescence image corresponding to the corresponding phase. Thereby completing the acquisition of the fluorescence image corresponding to one period phase. In addition, in practical applications, the specific value of the time interval can be adjusted according to the actual condition of the camera frame rate. Therefore, the invention collects the fluorescent image sequence of the pressure sensitive coating under the high-frequency pulsating pressure through the CCD camera with the low frame rate, avoids the error caused by the low frame rate of the camera, and improves the detection accuracy of the amplitude-frequency characteristic of the optical pressure sensitive coating. In addition, the image acquisition mode of the invention can avoid data storage/transmission congestion and ensure the smoothness of data storage and transmission.
step 4.1, to the fluorescence image QjPerforming image processing to obtain sum phaseCorresponding light intensity Ij;
The step 4.1 is specifically as follows:
fluorescence image QjW pixel points are provided; each pixel point corresponds to a light intensity value, therefore, w light intensity values are obtained in total, and the average value of the w light intensity values is the phase positionCorresponding light intensity Ij。
Step 4.2, obtaining a fluorescence image QjIn the acquisition process, the exposure time of the CCD camera 5 is obtained, and then a change curve between the real pressure value and the time on the surface of the PSP sample 2 is searched according to the exposure time of the CCD camera 5Obtaining a fluorescence image QjCorresponding true pressure value Pj;
And 4.3, when the sinusoidal pressure standing wave is not applied to the PSP sample 2, namely: the PSP sample 2 is brought to atmospheric pressure, and the image of the PSP sample 2 is taken, thereby obtaining a reference pressure PrefAnd light intensity at a reference pressure Iref;
since j is 1, 2.., n, when j is 1, one equation for a and B is obtained; when j is 2, an equation for a and B is obtained; by analogy, when j is equal to n, one equation for a and B is obtained; thus, a total of n equations for A and B are obtained; solving n equations about A and B by adopting a least square method to obtain final values of A and B;
a and B are constants, and the values of A and B and reference pressure PrefAnd light intensity at a reference pressure IrefSubstituting the calibration equation (1) to obtain a calibration equation (2):
wherein, in the calibration equation (2), P'jIs equal to the light intensity IjA corresponding calibration pressure value;
it should be noted that the main function of the calibration equation (1) is to obtain the values of a and B, so the pressure of the calibration equation (1) is the true pressure value P acquired by the dynamic pressure sensor 4j. After obtaining the values of a and B, further obtaining a calibration equation (2), wherein the main function of the calibration equation (2) is: calculating a calibration pressure value of the light intensity corresponding to each fluorescence image, since the calculated calibration pressure value is characteristic of the reaction coatingThe pressure value is not used, but the real pressure value collected by the dynamic pressure sensor 4.
continuously increasing the frequency of the sine pressure standing wave if the current frequency fzAmplitude A ofm(fz) Decreasing to 0, the current frequency fzIs the limiting frequency of the PSP coating; therefore, the amplitude-frequency characteristic of the PSP coating is detected.
Therefore, in the invention, the sinusoidal pressure standing wave with stable frequency and amplitude is obtained through the acoustic standing wave tube 1; then, according to the frequency of the sinusoidal pressure standing wave, the exposure time and the frame rate of the CCD camera are controlled, and the next phase of the sinusoidal pressure standing wave corresponding to the image shot each time is ensured, so that the obtained image sequence can just form a period; in the process, the excitation light source keeps stable and continuous luminescence, and the CCD camera starts image sequence collection of different phases.
The invention provides a CCD camera-based method for detecting the amplitude-frequency characteristic of an optical pressure-sensitive coating, which has the following advantages:
1. collecting a pressure sensitive coating fluorescence image sequence under high-frequency pulsating pressure by a low frame rate CCD camera;
2. keeping the amplitude of the sine pressure standing wave unchanged, and obtaining fluorescent image sequences under different conditions by changing the frequency of the sine pressure standing wave, thereby finally obtaining the amplitude-frequency characteristic of the coating, namely: the cut-off frequency and the limiting frequency of the PSP coating corresponding to the amplitude of the sinusoidal pressure standing wave.
3. The invention is a method for detecting the amplitude-frequency characteristic of the optical pressure-sensitive paint, which has the advantages of simple structure, low processing cost, strong usability and strong anti-interference capability, breaks through the prior method that PMT acquired by a single pixel is adopted to obtain the dynamic fluorescence intensity emitted by the optical pressure-sensitive paint, and realizes the advantage that the low frame rate CCD camera obtains the amplitude-frequency characteristic of the paint based on the PSP paint cut-off frequency and the limit frequency under higher frequency pressure obtained by the CCD camera.
It will be understood by those skilled in the art that all or part of the processes of the methods of the embodiments described above may be implemented by hardware associated with computer program instructions, and the above programs may be stored in a computer readable storage medium, and when executed, may include the processes of the embodiments of the methods described above. The storage medium may be a magnetic disk, an optical disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), or the like.
The foregoing is only a preferred embodiment of the present invention, and it should be noted that, for those skilled in the art, various modifications and improvements can be made without departing from the principle of the present invention, and such modifications and improvements should also be considered within the scope of the present invention.
Claims (2)
1. A CCD camera-based method for detecting the amplitude-frequency characteristic of an optical pressure-sensitive coating is characterized by comprising the following steps:
step 1, fixing a PSP sample wafer (2) at the position of a bottom cover of an acoustic standing wave tube (1), wherein after the acoustic standing wave tube (1) is screwed, the plane of the bottom cover is parallel to the plane of a sound source (11); a 1 st channel of the signal generator (8) is connected with a sound source (11) through a power amplifier (9), and a sound production end of the sound source (11) is connected with the acoustic standing wave tube (1);
the CCD camera (5) and the laser light source (3) are aligned to the PSP sample wafer (2) through an optical window (7) of the acoustic standing wave tube (1), and an optical filter (6) is fixedly arranged in front of a lens of the CCD camera (5); the 2 nd channel of the signal generator (8) is connected with the CCD camera (5); the 3 rd channel of the signal generator (8) is connected with the laser light source (3);
fixing a dynamic pressure sensor (4) on the cross section where the PSP sample wafer (2) is located;
the output end of the CCD camera (5) and the output end of the dynamic pressure sensor (4) are both connected to a computer (10);
the camera and the light source are aligned to the PSP sample through the optical window, and the dynamic pressure sensor is fixed on the cross section of the PSP sample and used for measuring the real-time dynamic pressure and providing a trigger signal;
step 2, starting from t being equal to 0, starting a signal generator (8) and a dynamic pressure sensor (4) at the same time; the dynamic pressure sensor (4) measures the real pressure value on the surface of the PSP sample wafer (2) in real time, and transmits the real pressure value on the surface of the PSP sample wafer (2) to the computer (10), so that the computer (10) obtains a change curve between the real pressure value on the surface of the PSP sample wafer (2) and time;
said signal sendingThe generator (8) has a continuous output period of T1Frequency of f1Is transmitted to a power amplifier (9), the power amplifier (9) drives a sound source (11) to emit a sinusoidal signal with a period T1Frequency of f1The sinusoidal sound wave of (2);
the sine sound wave acts through the acoustic standing wave tube (1), so that the surface of the PSP sample wafer (2) is subjected to stable period T1Frequency of f1A sinusoidal pressure standing wave of (a);
step 3, the signal generator (8) is connected with an external trigger interface of the laser light source (3) to control the working state of the laser light source (3); meanwhile, the signal generator (8) is connected with an external trigger interface of the CCD camera (5) to control the working state of the CCD camera (5); the specific control mode is as follows:
step 3.1, for a period T1Frequency of f1The sinusoidal pressure standing wave of (1) equally dividing a sinusoidal pressure standing wave into n phases, which are sequentially expressed as phasesPhase position…, phaseThe time intervals between adjacent phases are all
Presetting the exposure time t of each time of the CCD camerae(ii) a Wherein, te≤T1/(n-1);
Presetting the shooting time interval T of the CCD camera according to the frame rate of the CCD cameraC=m×T1+T1V (n-1); wherein m represents the number of camera shooting cycle intervals; shooting time interval T of cameraCGreater than the reciprocal of the maximum frame rate of the camera;
step 3.2, where t is tkTime of day, sineThe pressure standing wave is phase in the 1 st periodThe state of (1); the signal generator (8) controls the laser light source (3) to be turned on and the CCD camera (5) to be turned on simultaneously;
at the moment, after the laser light source (3) is turned on, the long-open mode is entered, namely: after the laser light source (3) is turned on, the stable power is kept to continuously emit light to irradiate the surface of the PSP sample wafer (2), so that the PSP coating is excited;
after the CCD camera (5) is opened, the 1 st exposure is carried out, and the exposure time is te(ii) a When the exposure time is reached, i.e.: when t is equal to tk+teWhen the CCD camera (5) is started, the CCD camera is closed; at this time, the CCD camera (5) outputs and phasesCorresponding fluorescence image Q1;
Step 3.3, start from the 1 st exposure start time of the camera, i.e. from t to tkThe time interval T of the CCD camera is counted from the beginningCAnd then, namely: t is tk+m×T1+T1When the sinusoidal pressure standing wave is in the m +1 th cycle phase (n-1)At this time, the CCD camera (5) is turned on, and the 2 nd exposure is performed for the exposure time te(ii) a When the exposure time is reached, i.e.: when t is equal to tk+m×T1+T1/(n-1)+teWhen the CCD camera (5) is started, the CCD camera is closed; at this time, the CCD camera (5) outputs and phasesCorresponding fluorescence image Q2;
Step 3.4, starting from the 2 nd exposure starting time of the camera, passing through the shooting time interval T of the CCD cameraCThen, the sinusoidal pressure standing wave is phase in the 2m +1 th periodIn the state (2), the CCD camera (5) is turned on to perform the 3 rd exposure, and the output and phase of the CCD camera (5) are adjustedCorresponding fluorescence image Q3;
And so on, in the process that the laser light source (3) keeps stable power and continuously emits light to irradiate the surface of the PSP sample wafer (2), the shooting time interval T of the CCD camera is set everyCThe CCD camera (5) is exposed once and outputs a fluorescence image corresponding to the current phase of the sine pressure standing wave; until the output of the CCD camera (5) and the last phaseCorresponding fluorescence image Qn;
To this end, when the frequency of the sinusoidal pressure standing wave is f1Then, a fluorescence image sequence corresponding to a complete cycle is obtained, which is respectively: and phaseCorresponding fluorescence image Q1And phaseCorresponding fluorescence image Q2…, and phaseCorresponding fluorescence image Qn;
Step 4, for each fluorescence image Qj1, 2., n, all treated in the following way:
step 4.1, to the fluorescence image QjPerforming image processing to obtain sum phaseCorresponding light intensity Ij;
Step 4.2, obtainingTaking a fluorescent image QjIn the acquisition process, the exposure time of the CCD camera (5) is searched, and the change curve between the real pressure value and the time on the surface of the PSP sample wafer (2) is searched according to the exposure time of the CCD camera (5) to obtain a fluorescence image QjCorresponding true pressure value Pj;
And 4.3, when the sinusoidal pressure standing wave is not applied to the PSP sample (2), namely: the PSP sample (2) is brought to atmospheric pressure, and the image of the PSP sample (2) is captured, so that a reference pressure P is obtainedrefAnd light intensity at a reference pressure Iref;
Step 5, the light intensity I is measuredjTrue pressure value PjReference pressure PrefAnd light intensity at a reference pressure IrefSubstituting the following calibration equation (1):
since j is 1, 2.., n, when j is 1, one equation for a and B is obtained; when j is 2, an equation for a and B is obtained; by analogy, when j is equal to n, one equation for a and B is obtained; thus, a total of n equations for A and B are obtained; solving n equations about A and B by adopting a least square method to obtain final values of A and B;
a and B are constants, and the values of A and B and reference pressure PrefAnd light intensity at a reference pressure IrefSubstituting the calibration equation (1) to obtain a calibration equation (2):
wherein, in the calibration equation (2), PjIs equal to the light intensity IjA corresponding calibration pressure value;
step 6, for each fluorescence image QjJ 1,2, n, each corresponding to a phaseValue of (D) and light intensity IjA value of (d); the light intensity IjSubstituting into the calibration equation (2) after steady state calibration to obtain the corresponding calibration pressure value Pj'; thus, a phase is obtainedAnd a calibrated pressure value Pj' a corresponding relationship value; since j is 1, 2.. times.n, a total of n sets of phases are obtainedAnd a calibrated pressure value Pj' a corresponding relationship value; therefore, n discrete points are drawn in a coordinate system with the abscissa as the phase and the ordinate as the calibration pressure value; fitting the n discrete points to form a pressure phase curve; analyzing the pressure phase curve to obtain an amplitude Am(f1) The meaning is as follows: when the frequency of the sinusoidal pressure standing wave is f1While obtaining the amplitude Am(f1);
Step 7, keeping the amplitude of the sinusoidal pressure standing wave stable and unchanged, and controlling the frequency of the sinusoidal pressure standing wave from f1Increase to f2Obtaining the corresponding amplitude A by adopting the mode of the step 2 to the step 6m(f2) (ii) a Judging the amplitude Am(f2) Whether or not to reduce to amplitude Am(f1) And if so, the amplitude Am(f2) The corresponding frequency is the cut-off frequency of the PSP coating; if not, the frequency of the sinusoidal pressure standing wave is further increased until the current frequency f is madexAmplitude A ofm(fx) Reduced to amplitude Am(f1) Up to half the time, at this time, the current frequency fxIs the cut-off frequency of the PSP coating;
continuously increasing the frequency of the sine pressure standing wave if the current frequency fzAmplitude A ofm(fz) Decreasing to 0, the current frequency fzIs the limiting frequency of the PSP coating; therefore, the amplitude-frequency characteristic of the PSP coating is detected.
2. The CCD camera-based method for detecting the amplitude-frequency characteristic of the optical pressure-sensitive paint as claimed in claim 1, wherein the step 4.1 is specifically as follows:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011186794.7A CN112304493B (en) | 2020-10-29 | 2020-10-29 | CCD camera-based optical pressure-sensitive paint amplitude-frequency characteristic detection method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011186794.7A CN112304493B (en) | 2020-10-29 | 2020-10-29 | CCD camera-based optical pressure-sensitive paint amplitude-frequency characteristic detection method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112304493A CN112304493A (en) | 2021-02-02 |
CN112304493B true CN112304493B (en) | 2022-04-15 |
Family
ID=74332292
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011186794.7A Active CN112304493B (en) | 2020-10-29 | 2020-10-29 | CCD camera-based optical pressure-sensitive paint amplitude-frequency characteristic detection method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112304493B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117990269A (en) * | 2024-04-07 | 2024-05-07 | 中国空气动力研究与发展中心低速空气动力研究所 | Dynamic calibration parameter acquisition method |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5355392A (en) * | 1991-12-13 | 1994-10-11 | Sharp Kabushiki Kaisha | Digital data detector for reducing errors due to frequency variations |
JP2006340446A (en) * | 2005-05-31 | 2006-12-14 | Hitachi Via Mechanics Ltd | Digital servo controller and laser processor |
CN102564694A (en) * | 2010-12-27 | 2012-07-11 | 中国航空工业第一集团公司沈阳空气动力研究所 | Device for measuring response time of pressure-sensitive paint on pressure change |
CN105205028A (en) * | 2015-09-15 | 2015-12-30 | 珠海全志科技股份有限公司 | Method and device for searching sampling point of data |
CN105940612A (en) * | 2014-02-03 | 2016-09-14 | 奥斯兰姆奥普托半导体有限责任公司 | Coding method for data compression of power spectra of an optoelectronic component and decoding method |
CN106054139A (en) * | 2016-06-01 | 2016-10-26 | 中国科学院电子学研究所 | Data acquisition method and device |
CN106411798A (en) * | 2016-08-31 | 2017-02-15 | 白昀 | Free space asynchronous image reception method of multiple optical signals |
CN106872100A (en) * | 2017-01-12 | 2017-06-20 | 西北工业大学 | Standing wave cast optical pressure sensitive coating Dynamic pressure calibration cabin |
CN208270147U (en) * | 2018-05-15 | 2018-12-21 | 中车青岛四方机车车辆股份有限公司 | Pressure sensitive coating calibrating installation and calibration system |
CN111076852A (en) * | 2020-01-03 | 2020-04-28 | 西北工业大学 | Rotor blade surface dynamic pressure optical measurement device and method based on photoelectric sensor |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB0302489D0 (en) * | 2003-02-04 | 2003-03-05 | Bae Systems Plc | Improvements relating to pressure sensitive paint |
JP5810667B2 (en) * | 2011-06-23 | 2015-11-11 | セイコーエプソン株式会社 | Optical device and detection apparatus |
CN104316262B (en) * | 2014-10-08 | 2017-05-03 | 西北工业大学 | Dual-purpose cabin for dynamic pressure calibration of optical pressure sensitive coatings |
CN111174958B (en) * | 2020-01-03 | 2021-09-07 | 西北工业大学 | Rotor blade surface universe dynamic pressure measuring device and method based on phase locking method |
-
2020
- 2020-10-29 CN CN202011186794.7A patent/CN112304493B/en active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5355392A (en) * | 1991-12-13 | 1994-10-11 | Sharp Kabushiki Kaisha | Digital data detector for reducing errors due to frequency variations |
JP2006340446A (en) * | 2005-05-31 | 2006-12-14 | Hitachi Via Mechanics Ltd | Digital servo controller and laser processor |
CN102564694A (en) * | 2010-12-27 | 2012-07-11 | 中国航空工业第一集团公司沈阳空气动力研究所 | Device for measuring response time of pressure-sensitive paint on pressure change |
CN105940612A (en) * | 2014-02-03 | 2016-09-14 | 奥斯兰姆奥普托半导体有限责任公司 | Coding method for data compression of power spectra of an optoelectronic component and decoding method |
CN105205028A (en) * | 2015-09-15 | 2015-12-30 | 珠海全志科技股份有限公司 | Method and device for searching sampling point of data |
CN106054139A (en) * | 2016-06-01 | 2016-10-26 | 中国科学院电子学研究所 | Data acquisition method and device |
CN106411798A (en) * | 2016-08-31 | 2017-02-15 | 白昀 | Free space asynchronous image reception method of multiple optical signals |
CN106872100A (en) * | 2017-01-12 | 2017-06-20 | 西北工业大学 | Standing wave cast optical pressure sensitive coating Dynamic pressure calibration cabin |
CN208270147U (en) * | 2018-05-15 | 2018-12-21 | 中车青岛四方机车车辆股份有限公司 | Pressure sensitive coating calibrating installation and calibration system |
CN111076852A (en) * | 2020-01-03 | 2020-04-28 | 西北工业大学 | Rotor blade surface dynamic pressure optical measurement device and method based on photoelectric sensor |
Non-Patent Citations (3)
Title |
---|
Characterization of Frequency Response of Pressure-Sensitive Paints;Tamao Sugimoto 等;《AIAA JOURNAL》;20170430;第55卷(第4期);第1460-1464页 * |
Experimental investigation on static/dynamic characteristics of a fast-response pressure sensitive paint;LI,RY 等;《CHINESE JOURNAL OF AERONAUTICS》;20180630;第31卷(第6期);第1198-1205页 * |
正弦波型高频动态压力光学校准系统及其应用;高丽敏 等;《航空学报》;20201025;第41卷(第10期);第(123667-1)-(123667-10)页 * |
Also Published As
Publication number | Publication date |
---|---|
CN112304493A (en) | 2021-02-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106501227B (en) | Measurement method based on pressure sensitive coating probe molecule fluorescence lifetime | |
CN102353519B (en) | Resolving power evaluation method for three-generation dim light image intensifier | |
TW486568B (en) | Method and apparatus for quantitative nondestructive evaluation of metal airfoils using high resolution transient thermography | |
US20120050537A1 (en) | Thermal imaging method and appratus for evaluating coatings | |
CN112304492B (en) | Pressure-sensitive paint pressure calibration method based on CCD camera long exposure light source stroboscopic | |
KR20020035106A (en) | Synthetic reference thermal imaging method | |
CN112304493B (en) | CCD camera-based optical pressure-sensitive paint amplitude-frequency characteristic detection method | |
CN104280670A (en) | Corona detection method based on solar-blind ultraviolet imager | |
CN104634261A (en) | Line laser source based medium plate shape inspection system and method | |
CN103335717B (en) | A kind of thermal infrared imager high precision temperature resistance drift temp measuring method based on becoming integral mode | |
CN112378576B (en) | Optical pressure sensitive coating pressure calibrating device based on CCD camera | |
WO2018232752A1 (en) | Method, apparatus and device for detecting substance | |
JP2005528985A (en) | Improvements related to pressure sensitive paints | |
CN102607805A (en) | Device and method for testing divergence angle and uniformity of illumination laser by using CCD (Charge Coupled Device) | |
CN109323851B (en) | Terahertz focal plane response rate and response rate nonuniformity test system and method | |
US7805251B2 (en) | Two slope reference for synthetic thermal time of flight imaging | |
CN112304491B (en) | Pressure-sensitive paint pressure calibrating device based on CCD camera long exposure light source stroboscopic | |
CN117268578A (en) | Device and method for non-contact temperature detection based on phosphorescent material | |
CN111366967A (en) | System and method for measuring irradiation dose rate of strong radiation field | |
CN114646663A (en) | High-efficiency measurement system and method for thermal radiation characteristics of materials with different thicknesses of high-temperature infrared hood | |
Scharf et al. | Four-band fiber-optic radiometry for determining the “true” temperature of gray bodies | |
CN113281002A (en) | Flow field measuring system and measuring method of transient flow field under flow fluctuation working condition | |
TWI494549B (en) | A luminance inspecting method for backlight modules based on multiple kernel support vector regression and apparatus thereof | |
CN111504497A (en) | Temperature measurement method based on fluorescent optical fiber | |
CN102519585A (en) | Apparatus for detecting strong-laser long distance surface intensity |
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 | ||
CB03 | Change of inventor or designer information | ||
CB03 | Change of inventor or designer information |
Inventor after: Gao Limin Inventor after: Ge Ning Inventor after: Jiang Heng Inventor after: Yang Guanhua Inventor after: Cai Ming Inventor before: Ge Ning Inventor before: Gao Limin Inventor before: Jiang Heng Inventor before: Yang Guanhua Inventor before: Cai Ming |