CN108229080A - The optimization method in unknown-model gas turbine blades digital radial subregion transillumination direction - Google Patents
The optimization method in unknown-model gas turbine blades digital radial subregion transillumination direction Download PDFInfo
- Publication number
- CN108229080A CN108229080A CN201810254256.3A CN201810254256A CN108229080A CN 108229080 A CN108229080 A CN 108229080A CN 201810254256 A CN201810254256 A CN 201810254256A CN 108229080 A CN108229080 A CN 108229080A
- Authority
- CN
- China
- Prior art keywords
- transillumination
- subregion
- blade
- defect
- image
- 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
- 238000000034 method Methods 0.000 title claims abstract description 29
- 238000005457 optimization Methods 0.000 title claims abstract description 26
- 230000007547 defect Effects 0.000 claims abstract description 71
- 230000008859 change Effects 0.000 claims abstract description 18
- 230000004044 response Effects 0.000 claims abstract description 17
- 238000012545 processing Methods 0.000 claims abstract description 4
- 239000000463 material Substances 0.000 claims description 7
- 238000009738 saturating Methods 0.000 claims description 5
- 238000013139 quantization Methods 0.000 claims description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims 1
- 238000004451 qualitative analysis Methods 0.000 abstract description 2
- 238000004445 quantitative analysis Methods 0.000 abstract description 2
- 238000001514 detection method Methods 0.000 description 9
- 238000010586 diagram Methods 0.000 description 9
- 238000005516 engineering process Methods 0.000 description 7
- 230000008901 benefit Effects 0.000 description 3
- PXHVJJICTQNCMI-UHFFFAOYSA-N Nickel Chemical compound [Ni] PXHVJJICTQNCMI-UHFFFAOYSA-N 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000002591 computed tomography Methods 0.000 description 2
- 230000007812 deficiency Effects 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000009659 non-destructive testing Methods 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 208000037656 Respiratory Sounds Diseases 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000005266 casting Methods 0.000 description 1
- 238000002485 combustion reaction Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 229910052759 nickel Inorganic materials 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 239000000825 pharmaceutical preparation Substances 0.000 description 1
- 229940127557 pharmaceutical product Drugs 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000002601 radiography Methods 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 229910000601 superalloy Inorganic materials 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
The invention discloses a kind of optimization methods in unknown-model gas turbine blades digital radial subregion transillumination direction, transillumination is carried out to flat panel detector, between the linear response area of gray scale for determining flat panel detector, to blade, each subregion carries out transillumination from different perspectives, the relation curve for drawing each subregion gray scale minimum of blade and angle change determines the optimal transillumination direction of each subregion, transillumination image after each subregion is converted splices to obtain the transillumination image of blade, each subregion overlapping region image is handled, across subregion defect overlapping volume is calculated using approximate data, the sum of volume the defects of being obtained from each subregion, volume the defects of subtracting overlapping region, obtain blade defect true volume size, realize the optimization processing in transillumination direction.The present invention improves the detector efficiency to tiny defect, avoids the design of complicated furniture, reduces testing cost, and versatility is high, and the positioning, qualitative and quantitative analysis for defect are laid a good foundation.
Description
Technical field
The invention belongs to industrial x-ray technical field of nondestructive testing, and in particular to a kind of unknown-model gas turbine blades number
The optimization method in word ray subregion transillumination direction.
Background technology
Gas turbine is a kind of using the gas continuously flowed as working medium impeller high speed rotation, so as to by the energy of fuel
It is changed into the rotary vane type dynamic power machine of useful work, is widely used in every field.Gas turbine blades generally use is accurate
The method manufacture of casting, and need to interact with the gas medium of high temperature and pressure high speed during the work time, therefore it is being made
It makes and is on active service the stage, be likely to generate shrinkage cavity and porosity, crackle, the various forms the defects of such as be mingled with.These defects will be tight
Ghost image rings working performance, service life and the security reliability of operation of gas turbine complete machine.Therefore, gas turbine leaf is studied
The detection technique of piece defect, timely and accurately finds defect, has to the normal work of entire gas turbine group very important
Meaning.
Since gas turbine blades belong to complex free curved surface class part, and making material is usually nickel base superalloy material
Material, therefore generally use ray carries out non-destructive testing it.
Industry CT (Computed Tomography) technology be capable of clear and intuitive accurate display defect size and
Shape, but equipment and testing expense are excessively high, a crop leaf measuring cost about 300,000 or so, and detection efficiency is low,.Industry
Although photograph detection has many advantages, such as imaging resolution height, high sensitivity, intuitive and reliable, need to consume in detection process
A large amount of film and treatment of pharmaceutical products etc., cost is higher, and due to needing to handle film, testing efficiency is relatively low.
DR (Digital Radiography) technology refers to directly is digitized x-ray imaging under the control of the computer
The X-ray information of workpiece is converted into digital signal, and by computer reconstruction figure by technology using digital radial flat panel detector
Picture and the technology for carrying out a series of post processing of image.DR technologies detect in resolution ratio already close to photograph, and transillumination cost
Low, detection efficiency is high, can digitize storage, therefore is detected using DR technologies.Since blade dimensions are larger and thickness distribution
Unevenness, so needing to carry out blade subregion, traditional transillumination method carries out transillumination to all subregions from same direction, then right
Transillumination image is spliced.Transillumination direction is not clearly required, it is random larger, and from same direction transillumination, Wu Fada
To the optimization in each subregion transillumination direction.Based on the above problem, the invention discloses a kind of gas turbine leaves of unknown-model
Piece digital radial subregion transillumination direction optimization method, selects the optimal transillumination direction of each subregion.
Invention content
In view of the above-mentioned deficiencies in the prior art, the technical problem to be solved by the present invention is that provide a kind of unknown-model
The optimization method in gas turbine blades digital radial subregion transillumination direction.
The present invention uses following technical scheme:
The optimization method in unknown-model gas turbine blades digital radial subregion transillumination direction carries out flat panel detector saturating
According between the linear response area of gray scale for determining flat panel detector, carrying out subregion to blade using galvanized wire, each subregion is never to blade
Transillumination is carried out with angle, the relation curve for drawing each subregion gray scale minimum of blade and angle change determines that each subregion is optimal
Transillumination direction, the transillumination image after each subregion is converted splices to obtain the transillumination image of blade, to each subregion overlapping region
Image is handled, and defect volume is determined according to defect gray scale and defect thickness relationship curve, using approximate data calculate across
Subregion defect overlapping volume, the sum of volume the defects of being obtained from each subregion, volume the defects of subtracting overlapping region to get to
Blade defect true volume size realizes the optimization processing in transillumination direction.
Specifically, changing tube voltage from minimum to peak, transillumination is carried out to flat panel detector, according under each voltage
The average gray value of the transillumination image obtained draws out average gray value and the relation curve of voltage squared, according to average ash
The range of linearity of angle value and voltage squared relation curve determines the tonal range corresponding to the linear response section of flat panel detector.
Further, according to the output gray level value of flat panel detector and the linear determining tablet of transmitted intensity of input
The output gray level value of detector and square linear, the transmitted intensity I at X-ray tube focal length F of tube voltageF
It calculates as follows:
IF=α ZiV2/F2
Wherein, α is proportionality coefficient, and Z is the atomic number of target, and i is tube current, and V is tube voltage.
Specifically, by measure on blade the datum level of integral shroud and leaf root part and flat panel detector in the horizontal plane institute into
Angle, the transillumination direction of blade is quantified, finds out transillumination image of the tonal range in the linear range of flat panel detector,
Draw each subregion gray scale minimum of blade and the relation curve of angle change.
Further, it determines the best transillumination angular range of each blade subregion, changes blade placing angle, best saturating
Increase by 0.5 ° in the range of irradiation angle every time from minimum to peak, change blade shroud part datum level and exist with flat panel detector
Angle carries out transillumination, the quantization in complete blade pair transillumination direction to blade in horizontal plane.
Further, on the basis of the highest region of intra vane thickness of transilluminated area, when determining that region scanning thickness is minimum
Transillumination angle for optimal transillumination direction, transillumination image is higher for transmitted intensity, the bigger negative-appearing image of gray value.
Further, the best transillumination angle of the 1st subregion of blade is 30 °~60 °.
Specifically, on the basis of the transillumination direction of the 1st subregion of blade to blade carry out transillumination, obtain blade the 2nd, 3,4,5,
Width of 6 subregions under the 1st subregion transillumination direction, then by blade the 2nd, 3,4,5,6 subregions under respective optimal transillumination direction
The transillumination image width adjusting that transillumination is obtained is to the width under the 1st subregion transillumination direction, then on the basis of the 1st subregion,
Blade is spliced.
Specifically, making one piece of voussoir with blade same material, carry out transillumination and obtain defect gray scale and defect thickness relationship
Curve, each subregion overlapping region is approximately triangle, is converted by the gray value of the transillumination image area limit of two subregions
The length on triangle both sides is obtained, angle is the difference of the transillumination angle of two subregions, and then obtains the area of triangle, is based on
The area of triangle is multiplied by pixel width by Finite Element Principle, obtains the pixel position overlapping volume, subregion is then obtained
The overlapping area of each pixel position of defect in boundary line carries out the defects of summation obtains obtained overlapping region total volume.
Further, the 1st subregion and volume v the defects of the 2nd subregion overlapping region12It calculates as follows:
Wherein, a be the flat panel detector pixel length of side, θ1For the best transillumination angle of the 1st subregion, θ2It is most preferably saturating for the 2nd subregion
Irradiation angle, h1kFor in the 1st sectional image at line of demarcation pixel the corresponding thickness of gray value, h2kFor in the 2nd sectional image
The corresponding thickness of the gray value of pixel at line of demarcation, i are the number of pixels of defect on the 1st, 2 subregion lines of demarcation.
Compared with prior art, the present invention at least has the advantages that:
The optimization method in unknown-model gas turbine blades digital radial subregion transillumination direction disclosed by the invention, passes through combustion
The optimization in each subregion transillumination direction of blade of gas turbine, improves the detector efficiency to tiny defect, by measuring integral shroud and leaf
The angle of root datum level and flat panel detector obtains blade transillumination direction, avoids the design of complex proprietary fixture, method letter
It is single easy and at low cost, image mosaic is carried out to each subregion overlapping region and eliminates overlapping volume, is blade defect into one
Step positioning, qualitatively and quantitatively processing provides premise.
Further, transillumination image gray scale is digital radial transillumination in the linear response tonal range of flat panel detector
Most basic premise, transillumination image of the gray scale in the linear response range of flat panel detector is only effective transillumination image, therefore needs
The tonal range corresponding to the linear response section of flat panel detector is obtained in advance.
Further, transillumination image gray scale definition graph within the linear response range of flat panel detector seems effective transillumination figure
Picture draws a subregion gray scale minimum and the relation curve of angle change, the angle corresponding to curve peak on this basis
It is worth for the best transillumination angle of subregion.Because the gray value of transillumination image is inversely proportional with thickness of workpiece, therefore gray scale minimum is located at work
Part thickness, the most thin transillumination angle of angle i.e. workpiece thickness corresponding to curve peak.Workpiece is thinner, to tiny
The detector efficiency of defect is higher, therefore selects the angle as the best transillumination angle of subregion.
Further, traditional arrangement is to carry out transillumination from same transillumination direction to each subregion of entire blade,
The scanning thickness for leading to each subregion is not minimum, can not realize the optimization of tiny defect detector efficiency, therefore to each point
Area selects different transillumination directions respectively, improves necessary to the detector efficiency of tiny defect.And transillumination direction is carried out excellent
Change the standard quantified firstly the need of determining transillumination direction.
Further, scanning thickness is thicker, lower to the detector efficiency of tiny defect, the tiny defect in the most thick region of blade
Detector efficiency is minimum, and it is maximum to reach tiny defect detector efficiency minimum for blade when the region scanning thickness is minimum, therefore with this
For optimal transillumination direction.
Further, transillumination image is spliced, must ensure the spliced size constancy of entire blade first, thus with
Under the transillumination direction of first subregion on the basis of each subregion transillumination image width of gained, change remaining each zoned width and arrive this
Value, it is ensured that spliced blade dimensions are constant.
Further, when carrying out volume calculating to trans-regional defect, if only defect volume in each subregion folded
Add, segmental defect volume can be caused to be repeated calculating, therefore take nearly algorithm herein, calculate defect overlapping volume, will repeat to count
The defects of calculation, volume cut, and reduced the error that volume calculates.
In conclusion the present invention improves the detector efficiency to tiny defect, the design of complicated furniture is avoided, is reduced
Testing cost, and versatility is high, the positioning, qualitative and quantitative analysis for defect are laid a good foundation.
Below by drawings and examples, technical scheme of the present invention is described in further detail.
Description of the drawings
Fig. 1 is gas turbine blades digital radial detection system connection diagram;
Fig. 2 is blade angle acquisition methods schematic diagram;
Fig. 3 is blade subregion schematic diagram;
Fig. 4 is image lap handling principle schematic diagram;
Fig. 5 schematic diagrames between the linear response area of gray scale of flat panel detector;
Fig. 6 is the Digit Display Angle Ruler schematic diagram with air level;
Fig. 7 is gray scale minimum and the relation curve schematic diagram of angle change;
Fig. 8 is defect thickness and gray-scale relation curve synoptic diagram.
Specific embodiment
It is directed in practical engineering application and digital radial detection is carried out to large-scale gas turbine, subregion first is carried out to blade,
Then from same direction, to blade, each subregion carries out transillumination, and transillumination direction is not up to optimal deficiency, the present invention provides
A kind of optimization method in unknown-model gas turbine blades digital radial subregion transillumination direction, since scanning thickness is bigger, transillumination
Sensitivity is lower, lower to the Detection capability of tiny defect, therefore the transillumination side that the criterion optimized is each subregion scanning thickness minimum
To for optimal transillumination direction, using the original partition method of blade and the transmitting illumination parameter of each subregion, obtained transillumination image
It is higher for transmitted intensity, the bigger negative-appearing image of gray value.
Referring to Fig. 1, turntable side sets radiographic source, opposite side setting detector, turntable, radiographic source and detector difference
It is connect with control system, blade is arranged on turntable, forms gas turbine blades digital radial detection system.
A kind of optimization method in unknown-model gas turbine blades digital radial subregion transillumination direction of the present invention, specific steps
It is as follows:
S1, it determines between the linear response area of flat panel detector gray scale
It keeps other parameter constant, changes tube voltage from minimum to peak, transillumination is carried out to flat panel detector.
Transmitted intensity I at X-ray tube focal length FFFor
IF=α ZiV2/F2 (1)
Wherein, α is proportionality coefficient, (1.1~1.4) × 10-6, Z is the atomic number of target, and i is tube current (mA), and V is
Tube voltage (kV);
The average gray value of the transillumination image obtained under each voltage is calculated, draws out average gray value and voltage
Square variation relation curve.
By above-mentioned formula (1) it is found that in the case where other parameter is constant, transmitted intensity and tube voltage square into just
Than, and the linear response of flat panel detector is linear for the output gray level value and the transmitted intensity of input of flat panel detector,
Therefore in the linear response range of flat panel detector, the output gray level value of flat panel detector also should be square linear with tube voltage
Relationship.
Using the average gray value of drafting and the relation curve of voltage squared, the range of linearity therein is found out, determines tablet
Tonal range corresponding to the linear response section of detector.
S2, each optimal transillumination angle of subregion is obtained
To blade, each subregion carries out transillumination from different perspectives, finds out tonal range in the linear range of flat panel detector
Transillumination image draws each subregion gray scale minimum of blade and the relation curve of angle change.
Transillumination direction Optimality Criteria is that the transillumination direction of scanning thickness minimum is optimal transillumination direction, since blade is complicated
Free form surface class part, the thickness of each point is different in transilluminated area, therefore using the highest region of part thickness in transilluminated area as base
Standard, it is optimal transillumination direction to find out transillumination angle when making the region scanning thickness minimum.
The transillumination direction of blade is corresponded with placement angle of the blade on turntable, thus by measure on blade integral shroud and
The datum level of leaf root part and flat panel detector angulation in the horizontal plane, quantify the transillumination direction of blade, such as
Shown in Fig. 2.
S3, subregion is carried out to blade using galvanized wire, the atomic number of lead is larger, big to the attenuation rate of ray, can be in transillumination
Highlighted feature is formed on image, in this, as blade subregion and the mark of image mosaic.Blade subregion schematic diagram is as shown in Figure 3.
Due to carrying out transillumination to same subregion from different directions, obtained transillumination image width is inconsistent, so needing
First on the basis of the transillumination direction of the 1st subregion, transillumination is carried out to blade, it is saturating in the 1st subregion to obtain 2,3,4,5,6 subregion of blade
It is then that 2,3,4,5, the 6 subregions transillumination image that transillumination is obtained under respective optimal transillumination direction is wide according to the width under direction
Degree is adjusted to the width under the first subregion transillumination direction, and then on the basis of the 1st subregion, blade is spliced.
S4, across subregion defect overlapping region volume calculate
One piece of voussoir with blade same material is made first, is carried out transillumination, is obtained between defect gray scale and defect thickness
Image deflects gradation conversion is thickness, using approximate data, calculates across subregion defect overlapping volume size by relation curve.
When the image adjacent to left and right splices, there are laps for spliced image, are lacked calculating across subregion
It needs to handle image lap when falling into volume, principle is as shown in Figure 4.
Overlapping region can approximation regard a triangle as, determined between gray scale and thickness when calculating defect volume
Correspondence then by the gray value of the transillumination image area limit of two subregions, can convert and obtain the length on triangle both sides
Degree, angle is the difference of the transillumination angle of two subregions.
It, can be in the hope of the area of triangle, based on finite element in the case that known to two edge lengths of triangle and its angle
The area of triangle is multiplied by pixel width by principle, then can obtain the pixel position overlapping volume, subregion is similarly obtained
The overlapping area of each pixel position of defect in boundary line, the defects of summing, then can obtain obtained overlapping region
Total volume will calculate the defects of the sum of the defects of obtaining volume subtracts overlapping region volume, you can obtain defect from each subregion
True volume size.
Purpose, technical scheme and advantage to make the embodiment of the present invention are clearer, below in conjunction with the embodiment of the present invention
In attached drawing, the technical solution in the embodiment of the present invention is clearly and completely described, it is clear that described embodiment is
Part of the embodiment of the present invention, instead of all the embodiments.The present invention being described and shown in usually here in attached drawing is real
Applying the component of example can be arranged and be designed by a variety of different configurations.Therefore, the present invention to providing in the accompanying drawings below
The detailed description of embodiment be not intended to limit the range of claimed invention, but be merely representative of the selected of the present invention
Embodiment.Based on the embodiments of the present invention, those of ordinary skill in the art are obtained without creative efforts
The every other embodiment obtained, shall fall within the protection scope of the present invention.
1st, between the linear response area of gray scale for determining flat panel detector
Using the Focal distance of 1000mm, the electric current of 0.6mA is constant, and voltage is transformed to from the minimum voltage 5kV of X-ray machine X
Ceiling voltage 220kV, increases 1kV every time, and transillumination is carried out to flat panel detector.Obtain transillumination image f5~f220Totally 216 width, meter
Calculate the average gray value of each transillumination image, draw transillumination image average gray value and shining voltage square relation curve
As shown in figure 5, by analyzing the range of linearity in curve, obtain between the linear response area of gray scale of flat panel detector as [A, B].
2nd, the best transillumination direction of each subregion is obtained
Blade is divide into upper part and lower part, and top half is divided into 1,2,3 three subregion, and lower half portion is divided into 4,5,6 three points
Area.Blade shroud and leaf root part datum level and flat panel detector are measured using the Digit Display Angle Ruler with air level in the horizontal plane
Angle, the basic structure of Digit Display Angle Ruler is as shown in fig. 6, reachable ± 0.1 ° of the angle precision measured by bevel protractor.
Here by taking the 1st subregion as an example, the acquisition methods in the best transillumination direction of each subregion of blade are told about.
The transillumination angle of 1st subregion is by measuring integral shroud part datum level and flat panel detector angle in the horizontal plane
It obtains.By visual observations, the best transillumination angle of the subregion can be substantially determined in the range of 30 °~60 °, revolving-turret changes
Become blade placing angle, making blade shroud part datum level, angle increases to from 30 ° in the horizontal plane with flat panel detector
60 °, increase by 0.5 ° every time, transillumination is carried out to the 1st subregion of blade, obtains 61 transillumination images of subregion 1.It produces each
It is as shown in Figure 7 according to the relation curve of gradation of image minimum and angle change.
Angle, θ in curve corresponding to peak is the best transillumination angle of the 1st subregion, corresponding to best transillumination angle
Transillumination direction is optimal transillumination direction.It can similarly obtain the optimal transillumination direction of 2,3,4,5,6 subregion of blade.
3rd, each subregion splicing
Under the optimal transillumination direction of the 1st subregion, transillumination is carried out to 1,2,3,4,5,6 subregion of blade, obtain at this time 2,3,
4th, 5,6 subregion transillumination image width are respectively wid12、wid13、wid14、wid15、wid16。
Under the optimal transillumination direction of the 2nd subregion, transillumination is carried out to blade subregion 2, obtains the transillumination figure of blade subregion 2
Picture, it is wid to adjust transillumination image width using the imresize functions in matlab12。
Similarly the transillumination image width that transillumination is obtained under respective optimal transillumination direction of subregion 3,4,5,6 is adjusted respectively
For wid13、wid14、wid15、wid16。
Then on the basis of the 1st subregion transillumination image, the transillumination image after 2,3,4,5,6 subregions are converted respectively is spliced to
On the transillumination image of 1st subregion, the transillumination image of entire blade is obtained.
4th, across subregion defect overlapping region volume calculates
One piece of voussoir with blade same material is made first, is carried out transillumination, is obtained between defect gray scale and defect thickness
Relation curve, as shown in Figure 8.
By taking subregion 1 and 2 overlapping region of subregion as an example, best transillumination angle is respectively θ1And θ2, flat panel detector pixel
The length of side is a, and defect shares i pixel on the 1st, 2 subregion lines of demarcation, in the 1st sectional image at line of demarcation pixel gray value
For gray11, gray12... ..., gray1i, corresponding thickness is respectively h11, h12... ..., h1i;In 2nd sectional image at line of demarcation
The gray value of pixel is gray21, gray22... ..., gray2i, corresponding thickness is respectively h21, h22... ..., h2i.Then the 1st subregion
With the 2nd subregion lap defect volume v12For
The present invention changes the size of tube voltage, makes the grey scale curve of flat panel detector by keeping other parameter constant,
The linear response section of flat panel detector is obtained, the transillumination image for ensureing to obtain later is within the range of linearity, for follow-up step
It is rapid that theoretical premise is provided.When determining the best transillumination direction of each subregion, visual observations are first passed through, obtain the best transillumination of blade subregion
The approximate range in direction, substantially reduces workload.Then by revolving-turret, change the transillumination angle of blade, to blade into
Row transillumination, and each transillumination image gray scale minimum and the relation curve of angle change are produced, obtain each point by curve
The best transillumination direction in area, method are simple and practicable.On the basis of the first subregion transillumination image, it is big to adjust each subregion transillumination image
It is small, then blade is spliced, ensure that spliced leaf image size is constant.By making voussoir, transillumination figure is obtained
As gray scale and the correspondence of thickness, it is thickness by image deflects gradation conversion, using approximate data, calculates across subregion defect
Overlapping volume size can greatly reduce the error of defect volume calculating.
More than content is merely illustrative of the invention's technical idea, it is impossible to protection scope of the present invention is limited with this, it is every to press
According to technological thought proposed by the present invention, any change done on the basis of technical solution each falls within claims of the present invention
Protection domain within.
Claims (10)
1. the optimization method in unknown-model gas turbine blades digital radial subregion transillumination direction, which is characterized in that visited to tablet
It surveys device and carries out transillumination, between the linear response area of gray scale for determining flat panel detector, subregion is carried out to blade using galvanized wire, it is each to blade
A subregion carries out transillumination from different perspectives, and the relation curve for drawing each subregion gray scale minimum of blade and angle change determines respectively
A optimal transillumination direction of subregion, the transillumination image after each subregion is converted splices to obtain the transillumination image of blade, to each point
Area overlapping area image is handled, and is determined defect volume according to defect gray scale and defect thickness relationship curve, is calculated using approximation
Method calculates across subregion defect overlapping volume, the sum of volume the defects of being obtained from each subregion, body the defects of subtracting overlapping region
Product realizes the optimization processing in transillumination direction to get to blade defect true volume size.
2. a kind of optimization side in unknown-model gas turbine blades digital radial subregion transillumination direction according to claim 1
Method, which is characterized in that change tube voltage from minimum to peak, transillumination is carried out to flat panel detector, according under each voltage
The average gray value of the transillumination image obtained draws out average gray value and the relation curve of voltage squared, according to average ash
The range of linearity of angle value and voltage squared relation curve determines the tonal range corresponding to the linear response section of flat panel detector.
3. a kind of optimization side in unknown-model gas turbine blades digital radial subregion transillumination direction according to claim 2
Method, which is characterized in that visited according to the output gray level value of flat panel detector and the linear determining tablet of transmitted intensity of input
Survey the output gray level value of device and square linear, the transmitted intensity I at X-ray tube focal length F of tube voltageFMeter
It calculates as follows:
IF=α ZiV2/F2
Wherein, α is proportionality coefficient, and Z is the atomic number of target, and i is tube current, and V is tube voltage.
4. a kind of optimization side in unknown-model gas turbine blades digital radial subregion transillumination direction according to claim 1
Method, which is characterized in that by measure on blade the datum level of integral shroud and leaf root part and flat panel detector in the horizontal plane institute into
Angle, the transillumination direction of blade is quantified, finds out transillumination image of the tonal range in the linear range of flat panel detector,
Draw each subregion gray scale minimum of blade and the relation curve of angle change.
5. a kind of optimization side in unknown-model gas turbine blades digital radial subregion transillumination direction according to claim 4
Method, which is characterized in that determine the best transillumination angular range of each blade subregion, change blade placing angle, in best transillumination
Increase by 0.5 ° in angular range every time from minimum to peak, change blade shroud part datum level with flat panel detector in water
Angle carries out transillumination, the quantization in complete blade pair transillumination direction to blade in plane.
6. a kind of optimization side in unknown-model gas turbine blades digital radial subregion transillumination direction according to claim 5
Method, which is characterized in that on the basis of the highest region of intra vane thickness of transilluminated area, determine saturating when region scanning thickness is minimum
Irradiation angle is optimal transillumination direction, and transillumination image is higher for transmitted intensity, the bigger negative-appearing image of gray value.
7. a kind of optimization side in unknown-model gas turbine blades digital radial subregion transillumination direction according to claim 5
Method, which is characterized in that the best transillumination angle of the 1st subregion of blade is 30 °~60 °.
8. a kind of optimization side in unknown-model gas turbine blades digital radial subregion transillumination direction according to claim 1
Method, which is characterized in that transillumination is carried out to blade on the basis of the transillumination direction of the 1st subregion of blade, obtains blade the 2nd, 3,4,5,6
Width of the subregion under the 1st subregion transillumination direction, then by blade the 2nd, 3,4,5,6 subregions under respective optimal transillumination direction thoroughly
It is right then on the basis of the 1st subregion according to the transillumination image width adjusting obtained to the width under the 1st subregion transillumination direction
Blade is spliced.
9. a kind of optimization side in unknown-model gas turbine blades digital radial subregion transillumination direction according to claim 1
Method, which is characterized in that make one piece of voussoir with blade same material, carry out transillumination and obtain defect gray scale and defect thickness relationship song
Line, each subregion overlapping region is approximately triangle, is converted by the gray value of the transillumination image area limit of two subregions
To the length on triangle both sides, angle is the difference of the transillumination angle of two subregions, and then obtains the area of triangle, based on having
The area of triangle is multiplied by pixel width by the first principle of limit, obtains the pixel position overlapping volume, subregion circle is then obtained
The overlapping area of each pixel position of defect on line carries out the defects of summation obtains obtained overlapping region total volume.
10. a kind of optimization in unknown-model gas turbine blades digital radial subregion transillumination direction according to claim 9
Method, which is characterized in that volume v the defects of the 1st subregion and the 2nd subregion overlapping region12It calculates as follows:
Wherein, a be the flat panel detector pixel length of side, θ1For the best transillumination angle of the 1st subregion, θ2For the best transillumination angle of the 2nd subregion
Degree, h1kFor in the 1st sectional image at line of demarcation pixel the corresponding thickness of gray value, h2kTo demarcate in the 2nd sectional image
The corresponding thickness of the gray value of pixel at line, i are the number of pixels of defect on the 1st, 2 subregion lines of demarcation.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810254256.3A CN108229080B (en) | 2018-03-26 | 2018-03-26 | Optimization method for model unknown gas turbine blade digital ray subarea transillumination direction |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810254256.3A CN108229080B (en) | 2018-03-26 | 2018-03-26 | Optimization method for model unknown gas turbine blade digital ray subarea transillumination direction |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108229080A true CN108229080A (en) | 2018-06-29 |
CN108229080B CN108229080B (en) | 2020-06-19 |
Family
ID=62660031
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810254256.3A Expired - Fee Related CN108229080B (en) | 2018-03-26 | 2018-03-26 | Optimization method for model unknown gas turbine blade digital ray subarea transillumination direction |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108229080B (en) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108956655A (en) * | 2018-07-06 | 2018-12-07 | 西安交通大学 | A kind of gas turbine high temperature leaf digital ray detection partition method of unknown-model |
CN110705097A (en) * | 2019-09-29 | 2020-01-17 | 中国航发北京航空材料研究院 | Method for removing duplicate of nondestructive testing data of rotating part of aircraft engine |
CN111915569A (en) * | 2020-07-09 | 2020-11-10 | 西安交通大学 | Method, equipment and medium for screening digital radiographic image areas of free-form surface type parts |
CN112446123A (en) * | 2019-08-28 | 2021-03-05 | 电子科技大学 | Measuring head pose planning method for blisk three-coordinate measuring machine |
CN113390902A (en) * | 2020-02-26 | 2021-09-14 | 丰田自动车株式会社 | Method and apparatus for inspecting membrane electrode assembly |
CN113670958A (en) * | 2021-09-02 | 2021-11-19 | 西安交通大学 | Gas turbine blade defect identification method based on X-ray attenuation coefficient difference |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002373328A (en) * | 2001-05-02 | 2002-12-26 | Agilent Technol Inc | Method and device for equalizing luminance of image picked up under point light source illumination |
CN104730091A (en) * | 2015-02-10 | 2015-06-24 | 西安交通大学 | Gas turbine blade defects extraction and analysis method based on region segmenting detection |
CN106645223A (en) * | 2016-12-13 | 2017-05-10 | 山西平阳重工机械有限责任公司 | Nondestructive testing method of carbon fiber composite multi-blade integral rotor |
CN106845011A (en) * | 2017-02-16 | 2017-06-13 | 西安交通大学 | A kind of large-scale gas turbine leaf digital ray subregion detection method |
CN106932416A (en) * | 2017-02-16 | 2017-07-07 | 西安交通大学 | Gas turbine blades internal flaw three-dimensional parameter extracting method based on digital radial |
-
2018
- 2018-03-26 CN CN201810254256.3A patent/CN108229080B/en not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002373328A (en) * | 2001-05-02 | 2002-12-26 | Agilent Technol Inc | Method and device for equalizing luminance of image picked up under point light source illumination |
CN104730091A (en) * | 2015-02-10 | 2015-06-24 | 西安交通大学 | Gas turbine blade defects extraction and analysis method based on region segmenting detection |
CN106645223A (en) * | 2016-12-13 | 2017-05-10 | 山西平阳重工机械有限责任公司 | Nondestructive testing method of carbon fiber composite multi-blade integral rotor |
CN106845011A (en) * | 2017-02-16 | 2017-06-13 | 西安交通大学 | A kind of large-scale gas turbine leaf digital ray subregion detection method |
CN106932416A (en) * | 2017-02-16 | 2017-07-07 | 西安交通大学 | Gas turbine blades internal flaw three-dimensional parameter extracting method based on digital radial |
Non-Patent Citations (1)
Title |
---|
李兵 等: "复合式叶片型面测量系统的误差分析与补偿", 《中国工程科学》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108956655A (en) * | 2018-07-06 | 2018-12-07 | 西安交通大学 | A kind of gas turbine high temperature leaf digital ray detection partition method of unknown-model |
CN108956655B (en) * | 2018-07-06 | 2019-12-20 | 西安交通大学 | Digital ray detection partitioning method for high-temperature blades of gas turbine with unknown model |
CN112446123A (en) * | 2019-08-28 | 2021-03-05 | 电子科技大学 | Measuring head pose planning method for blisk three-coordinate measuring machine |
CN112446123B (en) * | 2019-08-28 | 2022-12-30 | 电子科技大学 | Measuring head pose planning method for blisk three-coordinate measuring machine |
CN110705097A (en) * | 2019-09-29 | 2020-01-17 | 中国航发北京航空材料研究院 | Method for removing duplicate of nondestructive testing data of rotating part of aircraft engine |
CN113390902A (en) * | 2020-02-26 | 2021-09-14 | 丰田自动车株式会社 | Method and apparatus for inspecting membrane electrode assembly |
CN111915569A (en) * | 2020-07-09 | 2020-11-10 | 西安交通大学 | Method, equipment and medium for screening digital radiographic image areas of free-form surface type parts |
CN111915569B (en) * | 2020-07-09 | 2022-04-22 | 西安交通大学 | Method, equipment and medium for screening digital radiographic image areas of free-form surface type parts |
CN113670958A (en) * | 2021-09-02 | 2021-11-19 | 西安交通大学 | Gas turbine blade defect identification method based on X-ray attenuation coefficient difference |
Also Published As
Publication number | Publication date |
---|---|
CN108229080B (en) | 2020-06-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108229080A (en) | The optimization method in unknown-model gas turbine blades digital radial subregion transillumination direction | |
CN106845011B (en) | Large-scale gas turbine blade digital ray partition detection method | |
Sun et al. | A calibration method for misaligned scanner geometry in cone-beam computed tomography | |
US7015473B2 (en) | Method and apparatus for internal feature reconstruction | |
CN102711613A (en) | A computed tomography imaging process and system | |
KR102073803B1 (en) | Three-dimensional surface roughness evaluating device, three-dimensional surface roughness evaluating method, three-dimensional surface roughness data acquiring device, and three-dimensional surface roughness data acquiring method | |
CN102032875B (en) | Image-processing-based cable sheath thickness measuring method | |
CN106932416B (en) | Gas turbine blades internal flaw three-dimensional parameter extracting method based on digital radial | |
CN1797027A (en) | Method and system of three-dimensional computed tomography | |
CN108388713A (en) | Gas turbine blades search coverage adaptivenon-uniform sampling and exposure parameter optimization method | |
CN105806858A (en) | CT detection method and CT device | |
CN108956655B (en) | Digital ray detection partitioning method for high-temperature blades of gas turbine with unknown model | |
CN105510362A (en) | Micro-CT-based rice tillering characteristic nondestructive measurement device and measurement method thereof | |
CN104936526A (en) | X-ray CT device, and image reconfiguration method | |
US4907152A (en) | Method of improving CT resolution | |
CN104596449A (en) | Accurate interstice measuring method based on CT (computed tomography) images | |
CN108362716A (en) | A kind of historical relic material detection determination method and detection device based on Medical CT | |
CN104515473A (en) | Online diameter detection method of varnished wires | |
JPS6021440A (en) | Method for measuring distribution of local void rate | |
JP2010204060A (en) | X-ray inspection device, and inspection method of the same | |
CN110400253A (en) | A method of transmitting chromatography weight matrix is determined based on bilinear interpolation principle | |
CN115791830A (en) | Steel plate detection system, steel plate detection method and electronic equipment | |
CN109357664B (en) | Method for monitoring three-dimensional structure of disc type mirror surface in real time | |
Du et al. | Data reduction for time-of-flight small-angle neutron scattering with virtual neutrons | |
CN100545642C (en) | A kind of radiographic inspection compensation method |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200619 |
|
CF01 | Termination of patent right due to non-payment of annual fee |