CN103940517A - Method for obtaining three-dimensional temperature field in metal structure - Google Patents
Method for obtaining three-dimensional temperature field in metal structure Download PDFInfo
- Publication number
- CN103940517A CN103940517A CN201410163865.XA CN201410163865A CN103940517A CN 103940517 A CN103940517 A CN 103940517A CN 201410163865 A CN201410163865 A CN 201410163865A CN 103940517 A CN103940517 A CN 103940517A
- Authority
- CN
- China
- Prior art keywords
- metal construction
- temperature
- finite element
- user
- element model
- 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
Landscapes
- Radiation Pyrometers (AREA)
Abstract
The invention discloses a method for obtaining a three-dimensional temperature field in a metal structure. At first, the mapping relation between geometrical nodes of the front surface and the rear surface of a three-dimensional finite element model of the metal structure and nodes of a two-dimensional diagram of a thermal imager of the metal structure is determined by means of a calibration method, and a temperature distribution diagram of the surface of the metal structure can be obtained through the mapping relation; secondly, user assigned sections in the metal structure are positioned through geometrical information of the three-dimensional finite element model of the metal structure, then temperature values of the nodes on the user assigned sections in the metal structure are calculated according to a linear computation method, and therefore temperature distribution diagrams of all the user assigned sections are obtained; at last, the temperature distribution diagrams are combined together, so that the three-dimensional temperature field in the metal structure is obtained.
Description
Technical field
The present invention relates to infrared temperature measuring technology, particularly a kind of method of obtaining metal construction interior three-dimensional Temperature Distribution.
Background technology
Infrared thermography technology is widely used in all types of industries such as machinery, space flight and aviation, building materials, metallurgy, chemical industry, oil.Temperature survey is one of important content of modern industrial control system, and its performance quality directly affects the quality of controlled device.By the accurate measurement control heating-up temperature of temperature, make it to reach rapidly, exactly setting value for improving production quality, the very important realistic meaning of enhancing productivity.
In temperature controlled processes, temperature is also one of wherein main controlled parameter.The Temperature Distribution of typical components inside has directly embodied the stability of this equipment work, can help staff to observe intuitively, carries out temperature control according to actual conditions.The factor analysis such as physical dimension, material of Temperature Distribution and typical components, thus it also can reflect the design of typical components and manufacture rationally whether.From temperature controlled processes, the automatic control technology being based upon on the Temperature Distribution analytical technology of typical components will provide effective guarantee for security, high efficiency and the low contaminative of equipment operation.
But, due to the restriction of measurement means, the Temperature Distribution situation that existing surveying instrument can only measurement component surface, and utilize existing measurement means directly to know about the Temperature Distribution situation of components interior.In order to obtain the Temperature Distribution situation of components interior, feasible scheme is all to utilize the technology of D temperature profiles three-dimensional temperature field at present.Such as, there are many research workers to realize the reconstruction of three-dimensional temperature field with the contrary method for solving of Image Processing Structure radiation delivery, and in the coal-burning boiler three-dimensional temperature field in power station is visual, are applied.2007, DRESOR method based on Monte Carlo technology is calculated the radiation delivery process in stove, according to the linear relationship between three-dimensional temperature field in the flame temperature image of having set up and stove, adopt improved Tikhonov regularization method from flame temperature image, to rebuild the three-dimensional temperature field in stove.In recent years, people adopted acoustic method to propose some three-dimensional temperature field reconstruction techniques.These methods roughly can be divided into two classes: a class is the three-dimensional temperature field method for reconstructing based on two-dimentional Typical Planar reconstruction and three-dimensional interpolation.Another kind of is Direct Three-dimensional reconstruction of temperature field method.First kind method, first get several typical aspects in three-dimensional detected space, then arrange as far as possible uniformly sound wave transceiver in these typical aspect surroundings, form through many in plane equally distributed sound wave flight paths, calculate or measure the sound wave flight time in these typical aspects, utilize D temperature profiles algorithm to rebuild the temperature field of these typical aspects, the method of recycling three-dimensional interpolation just can draw the temperature field of relevant position in all the other aspects, thereby obtains the reconstructed results of whole three-dimensional temperature field.And Direct Three-dimensional reconstruction of temperature field method, first arrange sound wave transceiver in detected space surrounding, form through many equally distributed sound wave flight paths in same aspect and the different aspects in this tested region, and calculate the sound wave flight time, utilize Direct Three-dimensional reconstruction of temperature field algorithm to reconstruct three-dimensional temperature field.
At present these methods just rebuild for the three-dimensional temperature field of transparent gas or flame the feasible scheme providing, for material not the method for reconstructing of the three-dimensional temperature field of transparent solid structure also do not find pertinent literature.
Summary of the invention
The object of this invention is to provide a kind of method of the D temperature profiles interior three-dimensional temperature field by metal construction outside surface.
For reaching above object, the present invention takes following technical scheme to be achieved:
Obtain the method that metal construction interior three-dimensional temperature field distributes, it is characterized in that, comprise the steps:
The first step, first, sets up the three-dimensional finite element model of metal construction, and determines two outside surfaces that parallel with this model geometric coordinate system XY axle, i.e. front surface and rear surface, and wherein Z-direction is pointed to front surface by rear surface; Then adjust the camera site of thermal imaging system, the XY axle of the XY axle of pixel coordinate system of the thermal imaging system X-Y scheme of the forward and backward surface temperature distribution of guarantee metal construction and the geometric coordinate of three-dimensional finite element model system in the same way; Finally shoot respectively former and later two surface temperature distribution figure of metal construction, form the X-Y scheme of two these metal structure surface Temperature Distribution;
Second step, utilizes the X-Y scheme of three-dimensional finite element model and two metal structure surface Temperature Distribution of metal construction, sets up the mapping relations between the two, and concrete steps are as follows:
(1) set up the mapping relations between the geometric coordinate of metal construction three-dimensional finite element model and the thermal imaging system two dimensional image vegetarian refreshments coordinate of metal construction front surface Temperature Distribution, shown in (1), because X-Y scheme does not have Z axis coordinate, so do not consider in the time setting up mapping relations, formula (2) in like manner;
Wherein, X
1q, Y
1qthe transverse and longitudinal coordinate of a pixel on thermal imaging system X-Y scheme, unit: pixel; X
2q, Y
2qbe the transverse and longitudinal coordinate of the corresponding node of front surface in the three-dimensional finite element model of metal construction, unit is construction standard unit; X
0, Y
0the initial point of three-dimensional finite element model of the metal construction transverse and longitudinal coordinate in thermal imaging system X-Y scheme, unit: pixel; L
1for the three-dimensional finite element model of the metal construction length in thermal imaging system X-Y scheme, unit: pixel; L
2for the structure length of the three-dimensional finite element model of metal construction, unit is the physical unit that structure is used;
(2) set up the relation between the geometric coordinate of metal construction three-dimensional finite element model and the X-Y scheme pixel of temperature field, metal construction rear surface distribution, shown in (2):
Wherein, X
1h, Y
1hthe transverse and longitudinal coordinate of a node on thermal imaging system X-Y scheme, unit: pixel; X
2h, Y
2hbe the transverse and longitudinal coordinate of the corresponding node in rear surface in the three-dimensional finite element model of metal construction, unit is construction standard unit;
(3) utilize formula (1) and formula (2) that the pixel of the thermal imaging system X-Y scheme of the system point of metal construction front surface or rear surface and metal construction front and rear surfaces Temperature Distribution is mapped, and then, the temperature value of node the metal construction front surface obtaining from thermal imaging system or rear surface is given on the counter structure point of metal construction rear surface or rear surface; Finally, obtain the temperature profile of metal structure surface;
The 3rd step, the temperature profile based on metal structure surface is set up metal construction interior three-dimensional temperature field, and concrete steps are as follows:
1) utilize the geometric coordinate of metal construction three-dimensional finite element model to be, metal construction inner section is positioned, determine user's specified cross-section of metal construction inside; Concrete localization method is: first, obtain the degree of depth of user's specified cross-section, and the third direction maximum coordinates value of node in metal construction three-dimensional finite element model, deduct again the actual grade of user's input by this third direction maximum coordinates value, just can obtain the third direction coordinate of user's specified cross-section, navigate on user's specified cross-section;
2) temperature value of each node on calculating user specified cross-section, concrete grammar is: for any point on user's specified cross-section, find the corresponding node that transverse and longitudinal coordinate is identical with the transverse and longitudinal coordinate of this point according to its transverse and longitudinal coordinate in the front and rear surfaces structured data of the three-dimensional finite element model of metal construction; Then, from metal construction data, obtain their third direction coordinate; Utilize again the method for second step according to these two corresponding node coordinates, obtain their temperature value, utilize formula (3), calculate the temperature value of this node:
T
i=T
iq+(T
iq-T
ih)|Z-Z
iq|/|Z
iq-Z
ih| (3)
Wherein, T
ifor the temperature of any point on user's specified cross-section, T
iq, T
ihthe temperature that is respectively the corresponding node of temperature and rear surface of the corresponding node of front surface of the three-dimensional finite element model of metal construction, Z is user's specified cross-section third direction coordinate, Z
iq, Z
ihbe respectively the third direction coordinate of corresponding node in the front and rear surfaces of three-dimensional finite element model of metal construction; Finally, according to same method, calculate the temperature value of all nodes on user's specified cross-section;
3) utilize formula 1), find corresponding pixel points in X-Y scheme by all nodes on user's specified cross-section, and compose step 2 to corresponding pixel points) temperature value that obtains; Utilize the method for neighbor point interpolation, obtain the temperature value of other pixel on user's specified cross-section, draw the two-dimension temperature distribution plan of user's specified cross-section;
4) the two-dimension temperature distribution plan of the inner all user's specified cross-sections of metal construction is combined, just obtain the three-dimensional temperature field of metal construction inside.
Compared with prior art, the invention has the beneficial effects as follows the two-dimensional temperature field distributed data that utilizes the forward and backward surface of measurable metal construction, according to metal construction three-dimensional structure information, calculate its inside configuration user and specify temperature data and the user of any node on arbitrary section to specify the temperature profile on arbitrary section again.
Embodiment
A kind of method of obtaining the distribution of metal construction interior three-dimensional temperature field:
The first step, first, sets up the three-dimensional finite element model of metal construction, and determines that two outside surfaces that parallel with the XY axle of its geometric coordinate system are front surface and rear surface, and wherein Z-direction is pointed to front surface by rear surface.Then adjust the camera site of thermal imaging system, to ensure the XY axle of thermal imaging system two dimensional image vegetarian refreshments coordinate system of the forward and backward surface temperature distribution of metal construction and the XY axle of the geometric coordinate of metal construction three-dimensional finite element model system in the same way.Finally shoot respectively former and later two surface temperature distribution figure of metal construction, form the X-Y scheme of two these metal structure surface Temperature Distribution.
Second step, utilizes the X-Y scheme of three-dimensional finite element model and two these metal structure surface Temperature Distribution of metal construction, sets up the mapping relations between the two.Concrete steps are as follows:
First, set up the mapping relations between the geometric coordinate of metal construction three-dimensional finite element model and the thermal imaging system two dimensional image vegetarian refreshments coordinate of metal construction front surface Temperature Distribution, shown in (1).Because X-Y scheme does not have Z axis coordinate, so do not consider in the time setting up mapping relations.Formula (2) in like manner.
Wherein, X
1q, Y
1qthe transverse and longitudinal coordinate (unit: pixel) of a pixel on thermal imaging system X-Y scheme, X
2q, Y
2qthe transverse and longitudinal coordinate (unit is construction standard unit) of the corresponding node of front surface in the three-dimensional finite element model of metal construction, X
0, Y
0the initial point of three-dimensional finite element model of the metal construction transverse and longitudinal coordinate (unit: pixel) in thermal imaging system X-Y scheme, L
1for the three-dimensional finite element model of the metal construction length (unit: pixel) in thermal imaging system X-Y scheme, L
2for the structure length (unit is the physical unit that structure is used) of the three-dimensional finite element model of metal construction.
Then, set up the relation between the geometric coordinate of metal construction three-dimensional finite element model and the X-Y scheme pixel of temperature field, metal construction rear surface distribution, shown in (2).
Wherein, X
1h, Y
1hthe transverse and longitudinal coordinate (unit: pixel) of a node on thermal imaging system X-Y scheme, X
2h, Y
2hthe transverse and longitudinal coordinate (unit is construction standard unit) of the corresponding node in rear surface in the three-dimensional finite element model of metal construction, X
0, Y
0the initial point of three-dimensional finite element model of the metal construction transverse and longitudinal coordinate (unit: pixel) in thermal imaging system 2-dimentional photo, L
1for the three-dimensional finite element model of the metal construction length (unit: pixel) in thermal imaging system X-Y scheme, L
2for the structure length (unit is construction standard unit) of the three-dimensional finite element model of metal construction.
Utilize formula (1) and (2) that the pixel of the thermal imaging system X-Y scheme of the system point of metal construction front surface or rear surface and metal construction front and rear surfaces Temperature Distribution is mapped.And then, the temperature value of pixel the metal construction front surface obtaining from thermal imaging system or rear surface can be given on the counter structure point of metal construction rear surface or rear surface.Finally, obtain the temperature profile of metal structure surface, in order to facilitate us that the above-mentioned method of setting up mapping relations between the system point of metal construction front surface or rear surface and the pixel of the thermal imaging system X-Y scheme of metal construction front and rear surfaces Temperature Distribution is called to standardization.
The 3rd step, the Temperature Distribution based on metal structure surface is set up metal construction interior three-dimensional temperature field, and concrete steps are as follows:
1) utilize the geological information of metal construction three-dimensional finite element model, its inside configuration cross section is positioned.Concrete localization method is: first, obtain the degree of depth of user's specified cross-section, and the third direction maximum coordinates value of node in metal construction three-dimensional finite element model, deduct again the actual grade of user's input by this third direction maximum coordinates value, the third direction coordinate that just can obtain user's specified cross-section, has navigated on user's specified cross-section.Just can determine user's specified cross-section of metal construction inside based on the method.
2) when having determined after user's specified cross-section of metal construction inside, then calculate the temperature value of each node on user's specified cross-section, and then obtain the Temperature Distribution on user's specified cross-section.Specific practice is: for any point on user's specified cross-section, find the corresponding node that transverse and longitudinal coordinate is identical with the transverse and longitudinal coordinate of this point according to its transverse and longitudinal coordinate in the front and rear surfaces structured data of the three-dimensional finite element model of metal construction.Then, from metal construction data, obtain their third direction coordinate.Utilize again the method for second step according to these two corresponding node coordinates, obtain their temperature value, utilize formula (3), calculate the temperature value of this node.
T
i=T
iq+(T
iq-T
ih)|Z-Z
iq|/|Z
iq-Z
ih| (3)
Wherein, T
ifor the temperature of any point on user's specified cross-section, T
iq, T
ihthe temperature that is respectively the corresponding node of temperature and rear surface of the corresponding node of front surface of the three-dimensional finite element model of metal construction, Z is user's specified cross-section third direction coordinate, Z
iq, Z
ihbe respectively the third direction coordinate of corresponding node in the front and rear surfaces of three-dimensional finite element model of metal construction.Finally, according to same method, calculate the temperature value of all nodes on user's specified cross-section.
3) utilize formula (1) just can find corresponding pixel points in two-dimension picture by all nodes on user's specified cross-section, and compose temperature value to corresponding pixel points.Utilize the method for neighbor point interpolation, just can obtain the temperature value of other pixel on user's specified cross-section.And can draw the two-dimension temperature distribution plan of user's specified cross-section.
Adopt said method, just obtained the Temperature Distribution of the inner all user's specified cross-sections of metal construction, be exactly the distribution of three-dimensional temperature of metal construction, and then can also provide the two-dimension temperature distribution plan of any user's specified cross-section.The two-dimension temperature distribution plan of the inner all user's specified cross-sections of metal construction is combined, just obtain the three-dimensional temperature field of metal construction inside.
Illustrate said method of the present invention with a metal typical structure-rolling bearing below.
The first step, first, the geological information of the three-dimensional finite element model based on rolling bearing is demarcated camera or plane, finds the mapping relations of node in the node of rolling bearing front surface and thermal imaging system plane photo.Taking a node on rolling bearing front surface as example, the horizontal ordinate of this node is X
2=-4.09493 millimeters, ordinate is Y
2=-4.429293 millimeters, the horizontal ordinate X of three-dimensional finite element model initial point in thermal imaging system plane photo
0=158 pixels, ordinate is Y
0=111 pixels, the length L of three-dimensional finite element model in thermal imaging system plane photo
1=212 pixels, the structure length L of three-dimensional finite element model
2=90 millimeters.Utilize respectively two formulas in formula (1) can calculate the horizontal ordinate X of corresponding node on thermal imaging system plane photo
1=148 pixels, ordinate Y
1=100 pixels.
Then, utilize horizontal ordinate and the ordinate (148 of corresponding node on thermal imaging system plane photo, 100), by the upper node (148 of two-dimension temperature distribution of typical structure bearing upper surface, 100) temperature value 12oC gives on the system point (4.09493 ,-4.429293) of rolling bearing front surface.Utilize same method, can obtain rolling bearing front surface or rear surface other temperature value a little.Finally, obtain rolling bearing surface three dimension temperature profile.
Second step, based on the result of the first step, calculates the three-dimensional temperature field data of rolling bearing inside.Concrete steps are following two small steps: the first small step, three-dimensional finite element model geological information based on rolling bearing, rolling bearing inner section is positioned, concrete localization method is: first, 5 millimeters of the degree of depth of acquisition user specified cross-section, then, obtain 39.77539 millimeters of maximum third direction coordinate figures according to the coordinate of all surface node of rolling bearing, deduct the degree of depth of user's specified cross-section with this maximum ordinate, can obtain the third direction coordinate Z=34.77539 millimeter of user's specified cross-section, locate user's specified cross-section.
The second small step, when having determined after user's specified cross-section, calculates the two-dimension temperature data of all nodes on user's specified cross-section below.First, for any point on cross section, for example node (.313396,-3.16591), in the forward and backward surface structure data of the three-dimensional finite element model of rolling bearing, finding transverse and longitudinal coordinate according to its coordinate is all (.313396,-3.16591) corresponding node obtains their third direction coordinate Z from the three-dimensional structure data of rolling bearing
io=30.8138 and Z
ib=-38.7618.Then, utilize the method for the first step according to these two corresponding node coordinates, the temperature value that obtains them is respectively T
io=49.6 DEG C and T
ib=44.7 DEG C.Recycling formula (3), the temperature value that calculates this node (.313396 ,-3.16591) is 53.56159 DEG C.Finally, according to same method, calculate the temperature value of all nodes on user's specified cross-section.
The 3rd small step, utilizes formula (1) just can find corresponding pixel points by all nodes on user's specified cross-section, and composes temperature value to corresponding pixel points.For example, by node (.313396 ,-3.16591), utilize two formulas in formula (1) can calculate the horizontal ordinate X1=157 pixel of corresponding pixel points in two-dimension picture, ordinate Y1=103 pixel.Finally utilize the method for neighbor point interpolation, just can obtain the temperature value of other pixel on user's specified cross-section.And can draw the two-dimension temperature distribution plan of user's specified cross-section.
Adopting said method, just obtain finishing the Temperature Distribution of the inner all user's specified cross-sections of structure, is exactly the distribution of three-dimensional temperature of rolling bearing.Due to length, in table 1, provide the degree of depth and be that on user's specified cross-section of 5 millimeters, ordinate is from 200 pixels of the 100th pixel to the, horizontal ordinate is from the temperature value of 114 pixel regions of the 105th pixel to the.And then can also provide the two-dimension temperature distribution plan of any user's specified cross-section.
Table 1 (unit: pixel)
Transverse and longitudinal coordinate | 105 | 106 | 107 | 108 | 109 | 110 | 111 | 112 | 113 | 114 |
100 | 63.155 | 66.53999 | 69.395 | 72.165 | 74.05 | 75.92001 | 77.20499 | 78.00501 | 78.79 | 79.275 |
101 | 76.915 | 78.785 | 79.955 | 81.225 | 82.01001 | 82.795 | 83.38001 | 83.86501 | 84.07999 | 84.565 |
102 | 80.635 | 80.635 | 79.955 | 81.225 | 82.01001 | 82.795 | 83.38001 | 83.86501 | 84.07999 | 84.565 |
103 | 81.4 | 80.3 | 79.3 | 78.20001 | 77.6 | 77.185 | 76.885 | 76.685 | 76.685 | 76.9 |
104 | 68.085 | 65.91499 | 64.43 | 63.045 | 62.245 | 61.645 | 61.16 | 61.06 | 61.06 | 61.275 |
105 | 56.345 | 54.075 | 52.605 | 51.22 | 50.52 | 49.835 | 49.35 | 49.25 | 49.25 | 49.365 |
106 | 45.78 | 44.295 | 43.21 | 42.125 | 41.625 | 41.04 | 40.74 | 40.64 | 40.64 | 40.84 |
107 | 36.77 | 35.68501 | 34.9 | 34.115 | 33.715 | 33.415 | 33.03 | 32.93 | 32.93 | 33.03 |
108 | 31.145 | 30.345 | 29.66 | 29.06 | 28.76 | 28.375 | 28.09 | 27.99 | 27.99 | 27.99 |
109 | 25.79 | 25.09 | 24.32 | 23.92 | 23.535 | 23.335 | 23.05 | 22.95 | 22.95 | 22.95 |
110 | 22.45 | 21.95 | 21.465 | 20.98 | 20.78 | 20.495 | 20.395 | 20.295 | 20.295 | 20.295 |
111 | 19.31 | 18.91 | 18.525 | 18.325 | 18.125 | 17.84 | 17.74 | 17.825 | 17.825 | 17.74 |
112 | 17.34 | 17.04 | 16.755 | 16.555 | 16.455 | 16.17 | 16.255 | 16.07 | 16.07 | 16.07 |
113 | 15.685 | 15.485 | 15.47 | 15.37 | 15.27 | 15.17 | 14.985 | 14.985 | 14.985 | 14.985 |
114 | 14.6 | 14.4 | 14.485 | 14.57 | 14.385 | 14.285 | 14.1 | 14.1 | 14.1 | 14.1 |
115 | 14 | 13.9 | 13.9 | 13.985 | 13.8 | 13.7 | 13.7 | 13.785 | 13.785 | 13.6 |
116 | 13.5 | 13.13 | 13.13 | 13.215 | 13.215 | 13.3 | 13.3 | 13.385 | 13.2 | 13.2 |
117 | 13.115 | 12.93 | 12.93 | 13.015 | 13.1 | 13.1 | 13.185 | 13 | 13 | 12.915 |
118 | 12.915 | 12.73 | 12.73 | 12.815 | 12.815 | 12.985 | 12.985 | 12.715 | 12.715 | 12.715 |
119 | 12.63 | 12.63 | 12.53 | 12.53 | 12.615 | 12.615 | 12.615 | 12.615 | 12.615 | 12.615 |
120 | 12.43 | 12.43 | 12.43 | 12.43 | 12.43 | 12.43 | 12.43 | 12.43 | 12.43 | 12.43 |
121 | 12.33 | 12.33 | 12.33 | 12.33 | 12.33 | 12.245 | 12.245 | 12.33 | 12.33 | 12.33 |
122 | 12.415 | 12.415 | 12.415 | 12.23 | 12.23 | 12.415 | 12.415 | 12.415 | 12.415 | 12.415 |
123 | 12.415 | 12.415 | 12.23 | 12.23 | 12.23 | 12.23 | 12.23 | 12.415 | 12.23 | 12.23 |
124 | 12.23 | 12.23 | 12.23 | 12.23 | 12.23 | 12.23 | 12.23 | 12.23 | 12.23 | 12.13 |
125 | 12.13 | 12.045 | 12.045 | 12.045 | 12.045 | 12.045 | 12.13 | 12.13 | 12.13 | 12.13 |
126 | 12.13 | 12.045 | 12.045 | 12.045 | 12.045 | 12.13 | 12.13 | 12.13 | 12.13 | 11.945 |
127 | 12.045 | 12.045 | 12.045 | 12.045 | 12.045 | 12.13 | 12.13 | 12.13 | 12.13 | 11.945 |
128 | 12.045 | 12.045 | 12.045 | 12.045 | 12.13 | 12.13 | 11.945 | 12.13 | 12.13 | 12.03 |
129 | 12.13 | 12.13 | 12.13 | 12.13 | 12.13 | 12.13 | 12.13 | 11.945 | 12.03 | 12.03 |
130 | 12.13 | 12.13 | 11.945 | 11.945 | 11.945 | 12.13 | 12.03 | 12.03 | 12.03 | 11.845 |
131 | 11.945 | 12.03 | 12.03 | 12.03 | 12.03 | 12.03 | 12.03 | 12.03 | 11.845 | 11.845 |
132 | 11.845 | 11.845 | 12.03 | 12.03 | 12.115 | 11.845 | 11.845 | 11.845 | 11.845 | 11.845 |
133 | 11.845 | 11.845 | 11.845 | 11.845 | 11.93 | 11.845 | 11.845 | 11.845 | 11.93 | 11.93 |
134 | 11.845 | 11.845 | 11.845 | 11.845 | 11.93 | 11.845 | 12.03 | 12.115 | 11.93 | 11.93 |
135 | 12.03 | 12.03 | 11.845 | 11.845 | 11.93 | 11.93 | 12.115 | 12.115 | 12.115 | 11.93 |
136 | 12.03 | 12.03 | 11.845 | 11.93 | 11.93 | 11.93 | 11.93 | 12.115 | 12.115 | 11.93 |
137 | 12.03 | 11.845 | 11.845 | 11.745 | 11.745 | 11.93 | 11.93 | 11.93 | 11.93 | 11.845 |
138 | 11.845 | 11.845 | 11.845 | 11.745 | 11.745 | 11.745 | 11.745 | 11.93 | 11.845 | 11.845 |
139 | 11.845 | 11.845 | 11.93 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.66 |
140 | 11.845 | 11.93 | 11.93 | 11.93 | 11.93 | 11.93 | 11.745 | 11.745 | 11.745 | 11.745 |
141 | 11.93 | 11.93 | 11.93 | 11.93 | 11.93 | 11.93 | 11.93 | 11.745 | 11.745 | 11.745 |
142 | 11.93 | 11.93 | 11.93 | 11.93 | 11.845 | 11.845 | 11.845 | 11.745 | 11.745 | 11.745 |
143 | 11.93 | 11.845 | 11.845 | 11.845 | 11.845 | 11.845 | 11.845 | 11.745 | 11.745 | 11.745 |
144 | 11.845 | 11.845 | 11.845 | 11.845 | 11.845 | 11.845 | 11.66 | 11.745 | 11.745 | 11.745 |
145 | 11.93 | 11.745 | 11.745 | 11.745 | 11.745 | 11.66 | 11.66 | 11.745 | 11.745 | 11.745 |
146 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 |
147 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 |
148 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 |
149 | 11.93 | 11.93 | 11.93 | 11.93 | 11.93 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 |
150 | 11.93 | 11.93 | 11.93 | 11.845 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 |
151 | 11.745 | 11.93 | 11.93 | 11.66 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 |
152 | 11.66 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.93 |
153 | 11.66 | 11.66 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.93 |
154 | 11.66 | 11.66 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 |
155 | 11.845 | 11.845 | 11.745 | 11.93 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 |
156 | 11.745 | 11.93 | 11.93 | 11.93 | 11.93 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 |
157 | 11.745 | 11.93 | 11.93 | 11.93 | 11.93 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 |
158 | 11.745 | 11.93 | 11.93 | 11.93 | 11.745 | 11.66 | 11.66 | 11.66 | 11.745 | 11.745 |
159 | 11.93 | 11.93 | 11.93 | 12.015 | 11.83 | 11.745 | 11.66 | 11.66 | 11.745 | 11.745 |
160 | 11.845 | 11.845 | 11.93 | 11.83 | 11.83 | 11.745 | 11.745 | 11.66 | 11.745 | 11.745 |
161 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 |
162 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 |
163 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.83 | 11.83 | 11.745 |
164 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.66 | 11.93 | 11.93 | 12.015 | 11.93 |
165 | 11.93 | 11.745 | 11.745 | 11.745 | 11.745 | 11.845 | 11.845 | 11.93 | 11.93 | 11.93 |
166 | 11.845 | 11.93 | 11.745 | 11.745 | 11.745 | 11.745 | 11.845 | 11.845 | 11.845 | 11.93 |
167 | 11.845 | 11.93 | 11.745 | 11.745 | 11.745 | 11.66 | 11.845 | 11.845 | 11.845 | 11.845 |
168 | 11.845 | 11.93 | 11.745 | 11.745 | 11.66 | 11.66 | 11.66 | 11.845 | 11.845 | 11.845 |
169 | 11.93 | 11.93 | 11.745 | 11.66 | 11.66 | 11.66 | 11.66 | 11.845 | 11.845 | 11.845 |
170 | 11.93 | 11.745 | 11.745 | 11.745 | 11.66 | 11.66 | 11.66 | 11.845 | 11.845 | 11.845 |
171 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.66 | 11.66 | 11.845 | 11.845 | 11.845 |
172 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.66 | 11.845 | 11.845 | 11.93 | 11.745 |
173 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.845 | 11.845 | 11.93 | 11.745 | 11.745 |
174 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.845 | 11.845 | 11.745 | 11.745 | 11.745 |
175 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.66 | 11.66 | 11.745 | 11.745 | 11.745 |
176 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.66 | 11.66 | 11.745 | 11.745 | 11.745 |
177 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.66 | 11.66 | 11.745 | 11.745 | 11.745 |
178 | 11.66 | 11.745 | 11.745 | 11.745 | 11.745 | 11.845 | 11.845 | 11.745 | 11.745 | 11.745 |
179 | 11.845 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 |
180 | 11.93 | 11.93 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 |
181 | 11.93 | 11.93 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.93 | 11.93 |
182 | 11.93 | 11.93 | 11.745 | 11.745 | 11.745 | 11.745 | 11.745 | 11.93 | 11.93 | 11.845 |
183 | 11.93 | 11.93 | 11.93 | 11.745 | 11.745 | 11.745 | 11.93 | 11.845 | 11.845 | 11.845 |
184 | 11.93 | 11.93 | 11.93 | 11.93 | 11.93 | 11.93 | 11.93 | 11.845 | 11.845 | 11.845 |
185 | 11.745 | 11.93 | 11.93 | 11.93 | 11.93 | 11.93 | 11.93 | 11.845 | 11.845 | 11.845 |
186 | 11.845 | 11.93 | 11.93 | 11.93 | 11.93 | 11.93 | 11.93 | 11.845 | 11.845 | 11.845 |
187 | 11.845 | 11.93 | 11.93 | 11.93 | 11.93 | 11.93 | 11.93 | 11.845 | 11.845 | 11.845 |
188 | 11.845 | 11.845 | 11.93 | 11.93 | 11.93 | 11.93 | 11.93 | 11.845 | 11.845 | 11.845 |
189 | 11.845 | 11.845 | 11.845 | 11.845 | 11.93 | 11.93 | 11.66 | 11.66 | 11.66 | 11.66 |
190 | 11.845 | 11.845 | 11.845 | 11.845 | 11.845 | 11.845 | 11.66 | 11.66 | 11.66 | 11.66 |
191 | 11.845 | 11.845 | 11.845 | 11.845 | 11.845 | 11.845 | 11.845 | 11.66 | 11.66 | 11.66 |
192 | 11.845 | 11.845 | 11.845 | 11.845 | 11.845 | 11.845 | 11.845 | 11.66 | 11.66 | 11.66 |
193 | 11.845 | 11.845 | 11.845 | 11.845 | 11.845 | 11.845 | 11.845 | 11.66 | 11.66 | 11.745 |
194 | 11.845 | 11.845 | 11.845 | 11.845 | 11.845 | 11.845 | 11.845 | 11.66 | 11.66 | 11.745 |
195 | 11.845 | 11.845 | 11.845 | 11.845 | 11.845 | 11.845 | 11.845 | 11.66 | 11.66 | 11.845 |
196 | 11.76 | 11.76 | 11.76 | 11.845 | 11.845 | 11.845 | 11.845 | 11.845 | 11.845 | 11.845 |
197 | 11.86 | 11.86 | 11.945 | 11.945 | 12.03 | 11.93 | 11.93 | 11.845 | 11.845 | 12.03 |
198 | 12.045 | 12.045 | 12.045 | 11.945 | 12.03 | 12.03 | 12.03 | 12.03 | 12.03 | 12.03 |
199 | 12.045 | 12.045 | 12.045 | 12.045 | 11.945 | 12.03 | 12.03 | 12.03 | 12.03 | 12.03 |
200 | 12.045 | 12.045 | 12.045 | 12.045 | 12.13 | 12.13 | 12.13 | 12.13 | 12.215 | 12.215 |
Claims (1)
1. obtain the method that metal construction interior three-dimensional temperature field distributes, it is characterized in that, comprise the steps:
The first step, first, sets up the three-dimensional finite element model of metal construction, and determines two outside surfaces that parallel with this model geometric coordinate system XY axle: front surface and rear surface, and wherein Z-direction is pointed to front surface by rear surface; Then adjust the camera site of thermal imaging system, make the XY axle of pixel coordinate system of the forward and backward surface temperature distribution X-Y scheme of metal construction and the XY axle of three-dimensional finite element model geometric coordinate system in the same way; Finally shoot respectively this former and later two surface temperature distribution figure, form the X-Y scheme of two these metal structure surface Temperature Distribution;
Second step, utilizes the X-Y scheme of three-dimensional finite element model and two metal structure surface Temperature Distribution of metal construction, sets up the mapping relations between the two, and concrete steps are as follows:
(1) set up the mapping relations between the geometric coordinate of metal construction three-dimensional finite element model and the thermal imaging system two dimensional image vegetarian refreshments coordinate of metal construction front surface Temperature Distribution, shown in (1), because X-Y scheme does not have Z axis coordinate, so do not consider in the time setting up mapping relations, formula (2) in like manner;
Wherein, X
1q, Y
1qthe transverse and longitudinal coordinate of a pixel on thermal imaging system X-Y scheme, unit: pixel; X
2q, Y
2qbe the transverse and longitudinal coordinate of the corresponding node of front surface in the three-dimensional finite element model of metal construction, unit is construction standard unit; X
0, Y
0the initial point of three-dimensional finite element model of the metal construction transverse and longitudinal coordinate in thermal imaging system X-Y scheme, unit: pixel; L
1for the three-dimensional finite element model of the metal construction length in thermal imaging system X-Y scheme, unit: pixel; L
2for the structure length of the three-dimensional finite element model of metal construction, unit is the physical unit that structure is used;
(2) set up the relation between the geometric coordinate of metal construction three-dimensional finite element model and the X-Y scheme pixel of temperature field, metal construction rear surface distribution, shown in (2):
Wherein, X
1h, Y
1hthe transverse and longitudinal coordinate of a node on thermal imaging system X-Y scheme, unit: pixel; X
2h, Y
2hbe the transverse and longitudinal coordinate of the corresponding node in rear surface in the three-dimensional finite element model of metal construction, unit is construction standard unit;
(3) utilize formula (1) and formula (2) that the pixel of the thermal imaging system X-Y scheme of the system point of metal construction front surface or rear surface and metal construction front and rear surfaces Temperature Distribution is mapped, and then, the temperature value of node the metal construction front surface obtaining from thermal imaging system or rear surface is given on the counter structure point of metal construction front surface or rear surface; Finally, obtain the temperature profile of metal structure surface;
The 3rd step, the temperature profile based on metal structure surface is set up metal construction interior three-dimensional temperature field, and concrete steps are as follows:
1) utilize the geometric coordinate of metal construction three-dimensional finite element model to be, metal construction inner section is positioned, determine user's specified cross-section of metal construction inside; Concrete localization method is: first, obtain the degree of depth of user's specified cross-section, and the third direction maximum coordinates value of node in metal construction three-dimensional finite element model, deduct again the actual grade of user's input by this third direction maximum coordinates value, just can obtain the third direction coordinate of user's specified cross-section, navigate on user's specified cross-section;
2) temperature value of each node on calculating user specified cross-section, concrete grammar is: for any point on user's specified cross-section, find the corresponding node that transverse and longitudinal coordinate is identical with the transverse and longitudinal coordinate of this point according to its transverse and longitudinal coordinate in the front and rear surfaces structured data of the three-dimensional finite element model of metal construction; Then, from metal construction data, obtain their third direction coordinate; Utilize again the method for second step according to these two corresponding node coordinates, obtain their temperature value, utilize formula (3), calculate the temperature value of this node:
T
i=T
iq+(T
iq-T
ih)|Z-Z
iq|/|Z
iq-Z
ih| (3)
Wherein, T
ifor the temperature of any point on user's specified cross-section, T
iq, T
ihthe temperature that is respectively the corresponding node of temperature and rear surface of the corresponding node of front surface of the three-dimensional finite element model of metal construction, Z is user's specified cross-section third direction coordinate, Z
iq, Z
ihbe respectively the third direction coordinate of corresponding node in the front and rear surfaces of three-dimensional finite element model of metal construction; Finally, according to same method, calculate the temperature value of all nodes on user's specified cross-section;
3) utilize formula (1), find corresponding pixel points in X-Y scheme by all nodes on user's specified cross-section, and compose step 2 to corresponding pixel points) temperature value that obtains; Utilize the method for neighbor point interpolation, obtain the temperature value of other pixel on user's specified cross-section, draw the two-dimension temperature distribution plan of user's specified cross-section;
4) the two-dimension temperature distribution plan of the inner all user's specified cross-sections of metal construction is combined, just obtain the three-dimensional temperature field of metal construction inside.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410163865.XA CN103940517B (en) | 2014-04-22 | 2014-04-22 | A kind of method obtaining metal structure interior three-dimensional temperature field |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410163865.XA CN103940517B (en) | 2014-04-22 | 2014-04-22 | A kind of method obtaining metal structure interior three-dimensional temperature field |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103940517A true CN103940517A (en) | 2014-07-23 |
CN103940517B CN103940517B (en) | 2017-01-04 |
Family
ID=51188283
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410163865.XA Active CN103940517B (en) | 2014-04-22 | 2014-04-22 | A kind of method obtaining metal structure interior three-dimensional temperature field |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103940517B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107845100A (en) * | 2017-10-23 | 2018-03-27 | 浙江工业大学 | A kind of energy-storage system lithium battery heap 3 D stereo reconstructing temperature field method |
CN113984233A (en) * | 2021-10-18 | 2022-01-28 | 厦门大学 | Three-dimensional furnace wall temperature field visualization method and system |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080082288A1 (en) * | 2006-09-25 | 2008-04-03 | Raad Peter E | Thermography measurement system for conducting thermal characterization of integrated circuits |
US7365330B1 (en) * | 2006-09-18 | 2008-04-29 | Uchicago Argonne, Llc | Method for thermal tomography of thermal effusivity from pulsed thermal imaging |
CN101457267A (en) * | 2009-01-05 | 2009-06-17 | 中南大学 | Intelligent extraction method of blast burden temperature field isotherm |
US20100016711A1 (en) * | 2008-07-21 | 2010-01-21 | University Of South Carolina | Membrane-Deformation Mapping Technique |
-
2014
- 2014-04-22 CN CN201410163865.XA patent/CN103940517B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7365330B1 (en) * | 2006-09-18 | 2008-04-29 | Uchicago Argonne, Llc | Method for thermal tomography of thermal effusivity from pulsed thermal imaging |
US20080082288A1 (en) * | 2006-09-25 | 2008-04-03 | Raad Peter E | Thermography measurement system for conducting thermal characterization of integrated circuits |
US20100016711A1 (en) * | 2008-07-21 | 2010-01-21 | University Of South Carolina | Membrane-Deformation Mapping Technique |
CN101457267A (en) * | 2009-01-05 | 2009-06-17 | 中南大学 | Intelligent extraction method of blast burden temperature field isotherm |
Non-Patent Citations (2)
Title |
---|
ZHIFENG HUANG ET AL.: "Measurement on flame temperature and its 3D distribution in a 660 MWe arch-fired coal combustion furnace by visible image processing and verification by using an infrared pyrometer", 《MEASUREMENT SCIENCE AND TECHNOLOGY》 * |
马国军 等: "生物组织传热红外热成像的三维温度场重建", 《中国组织工程研究与临床康复》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107845100A (en) * | 2017-10-23 | 2018-03-27 | 浙江工业大学 | A kind of energy-storage system lithium battery heap 3 D stereo reconstructing temperature field method |
CN107845100B (en) * | 2017-10-23 | 2020-05-05 | 浙江工业大学 | Energy storage system lithium battery stack three-dimensional temperature field reconstruction method |
CN113984233A (en) * | 2021-10-18 | 2022-01-28 | 厦门大学 | Three-dimensional furnace wall temperature field visualization method and system |
Also Published As
Publication number | Publication date |
---|---|
CN103940517B (en) | 2017-01-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wei et al. | A novel 1D target-based calibration method with unknown orientation for structured light vision sensor | |
Tang et al. | Large deformation measurement scheme for 3D digital image correlation method | |
Tang et al. | Vision‐based three‐dimensional reconstruction and monitoring of large‐scale steel tubular structures | |
CN105021124A (en) | Planar component three-dimensional position and normal vector calculation method based on depth map | |
CN104990501B (en) | A kind of systematic parameter calibration method of there-dimensional laser scanning device | |
Zhang et al. | Novel calibration method for a multi-sensor visual measurement system based on structured light | |
Zhou et al. | A novel way of understanding for calibrating stereo vision sensor constructed by a single camera and mirrors | |
Zhou et al. | A measurement system based on internal cooperation of cameras in binocular vision | |
Zhang et al. | A novel 3D optical method for measuring and evaluating springback in sheet metal forming process | |
Wang et al. | An improved measurement model of binocular vision using geometrical approximation | |
Zhengchun et al. | Point cloud uncertainty analysis for laser radar measurement system based on error ellipsoid model | |
Jiang et al. | Calibration and uncertainty analysis of a combined tracking-based vision measurement system using Monte Carlo simulation | |
Wei et al. | Calibration approach for structured-light-stripe vision sensor based on the invariance of double cross-ratio | |
Zeng et al. | An error similarity-based error-compensation method for measurement in the nonuniform temperature field | |
Dong et al. | The calibration method of the circle-structured light measurement system for inner surfaces considering systematic errors | |
CN103940517A (en) | Method for obtaining three-dimensional temperature field in metal structure | |
Moghadam et al. | Heatwave: the next generation of thermography devices | |
CN103913114B (en) | Method for position association between binocular vision systems | |
García-Moreno et al. | Error propagation and uncertainty analysis between 3D laser scanner and camera | |
Li et al. | A performance evaluation method to compare the multi-view point cloud data registration based on ICP algorithm and reference marker | |
Franceschini et al. | Multivariate control charts for monitoring internal camera parameters in digital photogrammetry for LSDM (Large-Scale Dimensional Metrology) applications | |
Xue et al. | Bias estimation and correction for triangle-based surface area calculations | |
Fu et al. | A flexible approach to light pen calibration for a monocular-vision-based coordinate measuring system | |
Jia et al. | 3D temperature distribution model based on thermal infrared image | |
González-Jorge et al. | Metrological comparison of terrestrial laser scanning systems Riegl LMS Z390i and Trimble GX |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |