A method, device and storage medium for calculating target temperature field

Abstract

The present application relates to the technical field of target temperature field calculation, and in particular to a target temperature field calculation method, device and storage medium. The application method first constructs a finite element geometric network model of the target to be calculated, and divides the target to be calculated into a number of surface elements to be calculated; then determines the material characteristics and position coordinates of each surface element to be calculated; and for each surface element to be calculated, based on the position coordinates of the surface element to be calculated, determines the first angle between the normal of the surface element to be calculated and the angle of sunlight at the time to be calculated; finally, according to the material characteristics and the first angle of each surface element to be calculated, based on the material temperature texture file pre-generated for the area where the target to be calculated is located, determines the temperature value of each surface element to be calculated, and obtains the temperature field of the target to be calculated. The calculation method provided by the present application can greatly improve the calculation efficiency of the target temperature field.

Description

Technical Field

[0001] The present invention relates to the technical field of target temperature field calculation, and in particular to a target temperature field calculation method, device and storage medium. Background Art

[0002] Theoretical modeling of target radiation characteristics is a relatively theoretical research method that comprehensively considers various effects. According to the structure of the target, various thermophysical parameters of the material, and various environmental parameters of the environment, the target is divided into micro-elements according to certain requirements, various physical processes are considered, a physical model is established, and numerical solutions are performed through mathematical methods and computers to obtain the radiation characteristics of the target. The theoretical model of target radiation characteristics can provide a fast, cost-effective means for predicting the radiation characteristics of the target, the development of the target infrared system, simulation design, etc., and has important use value.

[0003] When conducting theoretical modeling of target radiation characteristics, it is first necessary to determine the temperature field of the target. However, in the prior art, the method for determining the target temperature field is inefficient, which seriously affects the research work on the target radiation characteristics.

[0004] Therefore, a target temperature field calculation method is urgently needed to solve the above technical problems. Summary of the invention

[0005] The embodiments of the present invention provide a target temperature field calculation method, device and storage medium, which can improve the calculation efficiency of the target temperature field.

[0006] In a first aspect, an embodiment of the present invention provides a method for calculating a target temperature field, comprising:

[0007] Constructing a finite element geometric network model of the target to be calculated, and dividing the target to be calculated into a plurality of surface elements to be calculated;

[0008] Determining the material characteristics and position coordinates of each of the face elements to be calculated;

[0009] For each of the face elements to be calculated, based on the position coordinates of the face element to be calculated, determining a first angle between the normal line of the face element to be calculated and the sunlight angle at the time to be calculated;

[0010] According to the material characteristics of each of the face elements to be calculated and the first angle, based on the material temperature texture file pre-generated for the area where the target to be calculated is located, the temperature value of each of the face elements to be calculated is determined to obtain the temperature field of the target to be calculated; wherein the material temperature texture file includes the face element temperature values ​​corresponding to the second angle between the face element normal of the target to be calculated corresponding to different material characteristics and the sunlight angle; the second angle is obtained by traversing the solar azimuth and pitch angle according to a preset angle interval.

[0011] In a possible design, the material temperature texture file of the area where the target to be calculated is located is obtained in the following manner:

[0012] Determine the environmental parameters corresponding to the area where the target to be calculated is located;

[0013] According to the environmental parameters, constructing energy balance equations corresponding to the second angle between the normal of the surface element and the sunlight angle for the surface elements with different material characteristics in the area;

[0014] Solve each of the energy balance equations, calculate the surface element temperature value corresponding to the second angle of the surface element with different material characteristics, and obtain the material temperature texture file of the area where the target to be calculated is located.

[0015] In a possible design, the environmental parameters correspond to different date ranges, and the material temperature texture file of the area where the target to be calculated is located contains material temperature texture sub-files of different date ranges;

[0016] The step of determining the temperature value of each of the face elements to be calculated based on a material temperature texture file generated in advance for the area where the target to be calculated is located comprises:

[0017] In the material temperature texture file of the area where the target to be calculated is located, a material temperature texture sub-file whose date range includes the time to be calculated is obtained, and the temperature value of each face element to be calculated is determined according to the material temperature texture sub-file.

[0018] In a possible design, constructing an energy balance equation for the face elements with different material characteristics at the second angle includes:

[0019] Establishing the energy balance equation based on a one-dimensional heat transfer model;

[0020] The one-dimensional heat transfer model is:

[0021]

[0022] Wherein, T is the temperature of the preset surface element, t is the time variable, x is the space variable, and α is the heat transfer coefficient;

[0023] The energy balance equation is:

[0024] E d +E sky +E ground +E cv +G h -Q it =0

[0025] Where: Ed is the solar radiation energy absorbed by the preset surface element at the second angle (w / m 2 );E sky is the atmospheric radiation energy absorbed by the preset surface element (w / m 2 );E ground is the ground radiation energy absorbed by the preset surface element, which is divided into two parts: ground reflection and ground radiation (w / m 2 );E cv is the energy absorbed by the preset surface element through convection heat exchange with the external air (w / m 2 ), which is mainly related to air convection velocity and atmospheric temperature; G h is the heat conduction flux (w / m 2 );Q it is the preset surface element radiation emittance (w / m 2 ), which is related to the self-emissivity and temperature of the preset surface element.

[0026] In a possible design, determining the temperature value of each of the facets to be calculated according to the material characteristics of each of the facets to be calculated and the first angle, based on a material temperature texture file pre-generated for the area where the target to be calculated is located, includes:

[0027] For each of the face elements to be calculated, obtaining a target material temperature texture file having the same material characteristics as the face element to be calculated from the material temperature texture file;

[0028] Selecting a second angle according to an angle difference between each second angle included in the target material temperature texture file and the first angle;

[0029] The temperature value of the bin to be calculated is determined according to the bin temperature value corresponding to the selected second angle.

[0030] In a possible design, selecting the second angle according to the angle difference between each second angle included in the target material temperature texture file and the first angle includes:

[0031] Determine whether the minimum angle difference among the angle differences is 0. If so, determine the second angle corresponding to the minimum angle difference as the selected second angle. Otherwise, determine the second angle with the minimum angle difference among the angle differences as the selected second angle.

[0032] In a possible design, determining the temperature value of the bin to be calculated according to the bin temperature value corresponding to the selected second angle includes:

[0033] If the minimum angle difference among the angle differences is 0, the temperature value of the bin corresponding to the selected second angle is determined as the temperature value of the bin to be calculated;

[0034] If the minimum angle difference among the angle differences is not 0, a linear interpolation calculation is performed according to the two selected second angles, the two surface element temperature values ​​corresponding to the second angles, and the first angle;

[0035] The bin temperature value obtained by linear interpolation calculation is determined as the temperature value of the bin to be calculated.

[0036] In a second aspect, an embodiment of the present invention further provides a target structure simulation device, comprising:

[0037] A construction module is used to construct a finite element geometric network model of the target to be calculated, and divide the target to be calculated into a plurality of surface elements to be calculated;

[0038] A first determination module is used to determine the material characteristics and position coordinates of each of the face elements to be calculated;

[0039] A second determination module is used to determine, for each of the surface elements to be calculated, a first angle between the normal line of the surface element to be calculated and the sunlight angle at the time to be calculated based on the position coordinates of the surface element to be calculated;

[0040] An acquisition module is used to determine the temperature value of each of the face elements to be calculated according to the material characteristics of each of the face elements to be calculated and the first angle, based on a material temperature texture file pre-generated for the area where the target to be calculated is located, so as to obtain the temperature field of the target to be calculated; wherein the material temperature texture file includes the face element temperature values ​​corresponding to the second angle between the face element normal of the target to be calculated corresponding to different material characteristics and the sunlight angle; the second angle is obtained by traversing the solar azimuth and pitch angle according to a preset angle interval.

[0041] In a third aspect, an embodiment of the present invention further provides a computing device, including a memory and a processor, wherein the memory stores a computer program, and when the processor executes the computer program, the method described in any embodiment of this specification is implemented.

[0042] In a fourth aspect, an embodiment of the present invention further provides a computer-readable storage medium having a computer program stored thereon, which, when executed in a computer, enables the computer to execute the method described in any embodiment of this specification.

[0043] An embodiment of the present invention provides a target temperature field calculation method, device and storage medium. The method, for a target to be calculated, first pre-generates a material temperature texture file of an area where the target to be calculated is located, the file including the surface element temperature value corresponding to the second angle between the surface element normal of the target to be calculated and the solar illumination angle corresponding to different material characteristics; then constructs a finite element geometric network model of the target to be calculated, divides the target to be calculated into a plurality of surface elements to be calculated, and determines the material characteristics and position coordinates of each of the surface elements to be calculated; then, for each of the surface elements to be calculated, based on the position coordinates of the surface element to be calculated, determines the first angle between the surface element normal to be calculated and the solar illumination angle at the time to be calculated; finally, according to the material characteristics and the first angle of each of the surface elements to be calculated, based on the pre-generated material temperature texture file, determines the temperature value of each of the surface elements to be calculated, and obtains the temperature field of the target to be calculated.

[0044] It can be seen that since the calculation method of the application has pre-generated the material temperature texture file of the area where the target to be calculated is located, there is no need to calculate the temperature value of each face element to be calculated in the target to be calculated separately. Instead, the temperature value of each face element to be calculated can be determined by querying the material temperature texture file, and then the temperature field of the target to be calculated can be determined. Therefore, the calculation method of the present application can greatly improve the calculation efficiency of the target temperature field. BRIEF DESCRIPTION OF THE DRAWINGS

[0045] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0046] Figure 1 is a flow chart of a target temperature field calculation method provided by one embodiment of the present invention;

[0047] Figure 2 is a geometric network model of a target to be calculated provided by an embodiment of the present invention;

[0048] Figure 3 It is a material temperature texture file of a region where a target to be calculated is located provided by an embodiment of the present invention;

[0049] Figure 4 It is a temperature field of a target to be calculated obtained by using the method of the present application provided by an embodiment of the present invention;

[0050] Figure 5 is a hardware architecture diagram of a computing device provided by an embodiment of the present invention;

[0051] Figure 6 It is a structural diagram of a target temperature field calculation device provided by one embodiment of the present invention. DETAILED DESCRIPTION

[0052] In order to make the purpose, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, rather than all the embodiments. Based on the embodiments in the present invention, all other embodiments obtained by ordinary technicians in this field without making creative work are within the scope of protection of the present invention.

[0053] As mentioned above, the theoretical model of target radiation characteristics can provide a fast, cost-effective means for predicting target radiation characteristics, developing target infrared systems, and simulating designs, etc. It has important practical value. When conducting theoretical modeling of target radiation characteristics, it is first necessary to determine the temperature field of the target.

[0054] However, in the prior art, when determining the target temperature field, it is necessary to calculate the angular coefficient and occlusion relationship between the finite surface elements that make up the target, which takes a long time to calculate. Especially when the area to be calculated is large and the number of targets is large, the calculation process is more complicated and the calculation efficiency is low, which seriously affects the research work on the target radiation characteristics.

[0055] Therefore, in order to solve the above technical problems, it is possible to consider pre-generating a material temperature texture file of the area where the target to be calculated is located, and then determining the temperature field of the target to be calculated by querying the file, thereby improving the calculation speed of the target temperature field.

[0056] The specific implementation of the above concept is described below.

[0057] Please refer to Figure 1 , an embodiment of the present invention provides a method for calculating a target temperature field, the method comprising:

[0058] Step 100, constructing a finite element geometric network model of the target to be calculated, and dividing the target to be calculated into a plurality of surface elements to be calculated;

[0059] Step 102, determining the material characteristics and position coordinates of each face element to be calculated;

[0060] Step 104, for each surface element to be calculated, based on the position coordinates of the surface element to be calculated, determining a first angle between the surface element normal to be calculated and the sunlight angle at the time to be calculated;

[0061] Step 106, according to the material characteristics and the first angle of each face element to be calculated, based on the material temperature texture file generated in advance for the area where the target to be calculated is located, determine the temperature value of each face element to be calculated, and obtain the temperature field of the target to be calculated; wherein the material temperature texture file includes the face element temperature value corresponding to the second angle between the face element normal of the target to be calculated corresponding to different material characteristics and the sunlight angle; the second angle is obtained by traversing the solar azimuth and pitch angle according to a preset angle interval.

[0062] In the embodiment of the present invention, since the material temperature texture file of the area where the target to be calculated is located is generated in advance, and the file includes the surface element temperature value corresponding to the surface element normal of the target to be calculated corresponding to different material characteristics and the second angle between the sunlight angle, therefore, for each surface element to be calculated in the target to be calculated, it is no longer necessary to calculate the temperature value of each surface element separately, but only the material characteristics, position coordinates and the first angle between the surface element normal and the sunlight angle at the time to be calculated of each surface element need to be determined, and the temperature value of each surface element to be calculated can be determined by querying the material temperature texture file, thereby determining the temperature field of the target to be calculated. The calculation method provided in this embodiment can greatly improve the calculation efficiency of the target temperature field.

[0063] Described below Figure 1 How the various steps are performed.

[0064] First, with respect to step 100, a finite element geometric network model of the target to be calculated is constructed, and the target to be calculated is divided into a plurality of surface elements to be calculated.

[0065] In this embodiment, the target to be calculated can be any tangible object such as a building, a mountain, a road, a lawn, etc. The target to be calculated can be one or more. After the target to be calculated is determined, a geometric network model of the target to be calculated can be constructed according to the three-dimensional shape of the target. After the network model is completed, geometric repair can be performed as needed, such as removing overlapping parts, to make the model more simplified and accurate.

[0066] Then, the constructed network model is subjected to facet decomposition, and the target to be calculated is divided into a number of facets to be calculated. The number of facets to be calculated can be determined according to the size of the target to be calculated and the user's requirements for the calculation accuracy of the target temperature field. It can be understood that the smaller the area of ​​the facet to be calculated, the higher the calculation accuracy and the longer the calculation time. Conversely, the calculation accuracy is reduced and the calculation time is correspondingly shortened. The user can determine the size of the facet to be calculated according to the requirements for calculation accuracy and calculation time, and then determine the number of facets to be calculated according to the area of ​​the target to be calculated. This application does not make specific restrictions.

[0067] Then, for step 102, the material characteristics and position coordinates of each face element to be calculated are determined.

[0068] In this embodiment, for the target to be calculated, the material characteristics of each facet to be calculated in the target to be calculated can be determined according to the building construction standard; wherein the material characteristics include the optical and thermal parameters of the material and the material structure. Specifically, the optical and thermal parameters of the material include visible light absorptivity, full-band emissivity, infrared band absorptivity, material density, specific heat capacity, and conductivity; the material structure includes the number of material layers, the material of each layer, and the material thickness; for multi-layer materials, the material density, specific heat capacity, conductivity, and material thickness of the inner layer material also need to be determined.

[0069] The position coordinates of the face element to be calculated include the coordinate values ​​of the face element center point in the X, Y, and Z directions in a preset coordinate system.

[0070] Next, for step 104 , for each facet to be calculated, based on the position coordinates of the facet to be calculated, a first angle between the normal line of the facet to be calculated and the sunlight angle at the time to be calculated is determined.

[0071] In this embodiment, the first angle includes the angle between the normal line of the face element to be calculated and the azimuth and elevation directions of the sun. At the same calculation time, the position of the face element to be calculated is different, and the first angle with the sun's illumination angle is also different, and the corresponding face element temperature value is different. For example, if the target to be calculated is a building, at 12 noon, the sun's altitude angle is the same, but for the face element at the bottom and the face element at the top, the position coordinates are different, so the normal lines of the two face elements are also different from the first angle with the sun's illumination angle, so the corresponding face element temperature values ​​are different, which is in line with the actual law.

[0072] It can be seen that this embodiment determines the first angle between the normal line of the surface element to be calculated and the sunlight angle at the time to be calculated based on the position coordinates of the surface element to be calculated, which can accurately correspond to the temperature value of each surface element.

[0073] Finally, for step 106, according to the material characteristics and the first angle of each face element to be calculated, based on the material temperature texture file pre-generated for the area where the target to be calculated is located, the temperature value of each face element to be calculated is determined to obtain the temperature field of the target to be calculated; wherein the material temperature texture file includes the face element temperature values ​​corresponding to the second angle between the face element normal of the target to be calculated corresponding to different material characteristics and the sunlight angle; the second angle is obtained by traversing the solar azimuth and pitch angle according to a preset angle interval.

[0074] It should be noted that the area where the target to be calculated is located can be determined according to the type and number of the target to be calculated. For example, if the target to be calculated is the residential buildings and streets in Haidian District, Beijing, then for the entire Beijing, the longitude and latitude and environmental parameters of each area are not much different. Therefore, the area can be determined as Haidian District, Beijing, or the area can be determined as Beijing, so as to increase the versatility of the pre-generated material temperature texture file and ensure that the file can cover every material in the target to be calculated. Of course, a larger area can also be determined, but it may reduce some of the calculation accuracy.

[0075] After the area is determined, all materials that may appear in the area can be determined based on the statistical data and building standards of the relevant parts. For example, the materials that may appear include grass, land, bricks, tiles, concrete, and multi-layer materials. In this embodiment, only one face element is taken for each material with different characteristics in the area, and the face element temperature value corresponding to the second angle between the normal line of each face element and the angle of sunlight is calculated to ensure that each face element to be calculated in the target to be calculated can match a second angle that is closest to the first angle when its normal line and the angle of sunlight at the calculation time are the first angle, thereby determining the temperature value of the face element to be calculated.

[0076] In addition, when determining the second angle, the azimuth angle of the sun is usually set to 0-360° and the pitch angle range is -90°-90°, so that the sunlight angle of 24 hours a day can be covered. The preset angle interval can be 3°, 5°, 10°, 15°, etc., and 5° is preferred in this embodiment. The second angle obtained by traversal can be 5° in azimuth and 0°, 5°, 10°, 15°, 20°, etc., respectively. This can not only ensure high calculation accuracy, but also save computer resources.

[0077] In some implementations, the material temperature texture file of the area where the target to be calculated is located is obtained in the following manner:

[0078] Step A1: Determine the environmental parameters corresponding to the area where the target to be calculated is located;

[0079] Step A2: According to the environmental parameters, for the facets with different material characteristics in the region, respectively construct the energy balance equation corresponding to the second angle between the facet normal and the sunlight angle;

[0080] Step A3: Solve each energy balance equation, calculate the surface element temperature value corresponding to the second angle of the surface element with different material characteristics, and obtain the material temperature texture file of the area where the target to be calculated is located.

[0081] In step A1, the environmental parameters include the longitude and latitude of the area, the ambient temperature, humidity, wind speed and other meteorological parameters. The longitude and latitude are the longitude and latitude corresponding to the center of the area, and the meteorological parameters can be the average parameters for a certain day or a certain time period. After determining the environmental parameters and the second angle, the calculation of solar radiation energy, solar altitude angle and atmospheric radiation energy can be carried out.

[0082] With respect to step A2, in some embodiments, constructing an energy balance equation for facets with different material characteristics at the second angle includes:

[0083] Establish the energy balance equation based on the one-dimensional heat transfer model;

[0084] The one-dimensional heat transfer model is:

[0085]

[0086] Where T is the temperature of the preset surface element, t is the time variable, x is the space variable, and α is the heat transfer coefficient;

[0087] The energy balance equation is:

[0088] E d +E sky +E ground +E cv +G h -Q it =0

[0089] Where: E d is the solar radiation energy absorbed by the preset surface element at the second angle (w / m 2 );E sky is the atmospheric radiation energy absorbed by the preset surface element (w / m 2 );E ground It is the ground radiation energy absorbed by the preset surface element, which is divided into two parts: ground reflection and ground radiation (w / m 2 );E cv is the energy absorbed by the preset surface element through convection heat exchange with the external air (w / m 2 ), which is mainly related to air convection velocity and atmospheric temperature; G h is the heat conduction flux between the outer surface of the preset surface element and the inner surface of the preset surface element (w / m 2 );Q it is the preset surface element radiation emittance (w / m 2 ), which is related to the self-emissivity and temperature of the preset surface element.

[0090] In this embodiment, since the preset surface element is very small, the heat conduction of the surface element to be calculated along the vertical and horizontal directions can be ignored, and only the heat flux in the thickness direction of the surface element is considered. Therefore, the heat transfer model can be simplified into a one-dimensional model to reduce the computational complexity of the temperature field calculation.

[0091] In addition, since the surface temperature of the surface element tends to be constant in actual calculation, the heat conduction flux G of the inner and outer surfaces in the thickness direction of the preset surface element is calculated. h The inner surface temperature of the element can be set to a constant value to simplify the calculation. Of course, the inner surface temperature can also be determined by solving the energy balance equation of the inner surface of the element, and then the heat conduction flux of the inner and outer surfaces in the thickness direction of the element is calculated.

[0092] For step A3: solve the surface element temperature value of each material surface element corresponding to the second angle in turn, obtain the surface element temperature values ​​of all material surface elements corresponding to the second angle, and form a material temperature texture file of the area where the target to be calculated is located.

[0093] In some embodiments, the environmental parameters correspond to different date ranges, and the material temperature texture file of the area where the target to be calculated is located includes material temperature texture sub-files of different date ranges;

[0094] Based on the material temperature texture file generated in advance for the area where the target to be calculated is located, the temperature value of each face element to be calculated is determined, including:

[0095] In the material temperature texture file of the area where the target to be calculated is located, a material temperature texture subfile whose date range includes the time to be calculated is obtained, and the temperature value of each face element to be calculated is determined according to the material temperature texture subfile.

[0096] In this implementation, the environmental parameters are refined according to the date range. For example, the environmental parameters can be determined according to the month or season, and the material temperature texture subfiles corresponding to different date ranges are generated. For example, the material temperature texture subfiles with the date range from January to December or the first quarter to the fourth quarter are generated respectively. Then, when the target temperature field on January 15th of the region needs to be calculated, the material temperature texture subfile of January can be obtained from the material temperature texture file first, and then the temperature value of each face element to be calculated can be determined according to the material temperature texture subfile. The face element temperature value determined in this way is more accurate.

[0097] In some embodiments, according to the material characteristics and the first angle of each face element to be calculated, based on a material temperature texture file pre-generated for the area where the target to be calculated is located, determining the temperature value of each face element to be calculated includes:

[0098] Step B1: for each face element to be calculated, a target material temperature texture file having the same material characteristics as the face element to be calculated is obtained from the material temperature texture file;

[0099] Step B2: selecting a second angle according to the angle difference between each second angle included in the target material temperature texture file and the first angle;

[0100] Step B3: Determine the temperature value of the bin to be calculated according to the bin temperature value corresponding to the selected second angle.

[0101] For step B1, the material temperature texture file contains the temperature textures of various materials. Therefore, when the material of the face element to be calculated is determined, for example, the material is concrete, then only the temperature texture file with the material of concrete is obtained from the material temperature texture sub-file, which can improve the query and matching speed.

[0102] In some embodiments, step B2 includes:

[0103] Determine whether the minimum angle difference among the angle differences is 0. If so, determine the second angle corresponding to the minimum angle difference as the selected second angle. Otherwise, determine the second angle with the smallest angle difference among the angle differences as the selected second angle.

[0104] It should be noted that the angle difference includes the angle difference in azimuth and the angle difference in elevation. When the angle differences in azimuth and elevation are both minimum, it is the minimum angle difference.

[0105] In some embodiments, step B2 includes:

[0106] If the minimum angle difference among the angle differences is 0, the temperature value of the face element corresponding to the selected second angle is determined as the temperature value of the face element to be calculated; for example, if the selected second angle and the first angle are both 5° in azimuth and 20° in elevation, then the face element temperature of 25°C corresponding to the second angle can be directly assigned as the temperature of the face element to be calculated.

[0107] If the minimum angle difference among the angle differences is not 0, a linear interpolation calculation is performed based on the two selected second angles, the surface element temperature values ​​corresponding to the two second angles, and the first angle;

[0108] The surface element temperature value obtained by linear interpolation calculation is determined as the temperature value of the surface element to be calculated.

[0109] For example, the two selected second angles are 5° in azimuth, 20° in elevation and 25° in elevation, and the corresponding surface element temperatures are 25°C and 26°C respectively, while the first angle is 5° in azimuth and 23°C in elevation. Then, through sequential linear interpolation calculations, it can be determined that the temperature value of the surface element to be calculated is 25.6°C. If the second angle and the first angle are different in azimuth and elevation, two interpolation calculations are required to obtain the temperature value of the surface element to be calculated.

[0110] The following is a specific example to illustrate the temperature field calculation effect of the method of the present application.

[0111] In this embodiment, the temperature field of an office building in Beijing on January 15 needs to be calculated. The three-dimensional geometric network model of the office building is as follows: Figure 2 As shown, using the network model, the office building can be decomposed into several facets to be calculated.

[0112] In order to quickly determine the temperature field of the office building, this embodiment first generates material temperature texture files of different characteristic materials in the Beijing area in January, such as Figure 3 As shown, Figure 3 In the file, each row of texture represents the temperature value of a material surface element at different second angles. The material temperature texture file includes the temperature values ​​of all surface elements of different characteristic materials in Beijing at different second angles.

[0113] For each face element to be calculated in the office building, the first angle between the face element normal and the sunlight angle at the calculation time can be used to calculate the face element. Figure 3 The material temperature texture file in the file determines the temperature value of the face element to be calculated, and the temperature values ​​of all the face elements to be calculated are calculated in turn to form the temperature field of the office building, such as Figure 4 shown.

[0114] We know that for any object, the side facing the sun absorbs more sunlight, so its temperature will be higher than the side facing away from the sun. In the height direction, due to the high wind speed, the convection heat exchange between the surface of the object and the air is obvious, and the heat loss is relatively large. Therefore, the surface temperature of the upper layer of the object is lower than that of the lower layer. Figure 4 It can be seen that the temperature field of the target calculated by this method has a significantly higher temperature on the side facing the sun than on the side facing away from the sun, and within the same side of the target, the temperature value of the surface element gradually decreases in the vertical upward direction, which is in line with the laws of nature.

[0115] It can be seen that the calculation method of the present application can not only improve the calculation efficiency, but also the calculation results can truly reflect the surface element temperature values ​​at different positions of the target.

[0116] like Figure 5 , Figure 6 As shown, an embodiment of the present invention provides a target temperature field calculation device. The device embodiment can be implemented by software, or by hardware or a combination of software and hardware. From the hardware level, Figure 5 As shown, it is a hardware architecture diagram of a computing device where a target temperature field computing device provided by an embodiment of the present invention is located, except Figure 5 In addition to the processor, memory, network interface, and non-volatile memory shown in the figure, the computing device in which the device is located in the embodiment may also generally include other hardware, such as a forwarding chip responsible for processing messages, etc. Taking software implementation as an example, Figure 6 As shown, as a device in a logical sense, the CPU of the computing device in which it is located reads the corresponding computer program in the non-volatile memory into the memory and runs it. This embodiment provides a target temperature field calculation device, including:

[0117] A construction module 300 is used to construct a finite element geometric network model of the target to be calculated, and divide the target to be calculated into a plurality of surface elements to be calculated;

[0118] The first determination module 302 is used to determine the material characteristics and position coordinates of each face element to be calculated;

[0119] The second determination module 304 is used to determine, for each facet to be calculated, based on the position coordinates of the facet to be calculated, a first angle between the facet normal to be calculated and the sunlight angle at the time to be calculated;

[0120] The acquisition module 306 is used to determine the temperature value of each face element to be calculated according to the material characteristics and the first angle of each face element to be calculated, based on the material temperature texture file generated in advance for the area where the target to be calculated is located, and obtain the temperature field of the target to be calculated; wherein the material temperature texture file includes the face element temperature value corresponding to the second angle between the face element normal of the target to be calculated corresponding to different material characteristics and the sunlight angle; the second angle is obtained by traversing the solar azimuth and pitch angle according to a preset angle interval.

[0121] In an embodiment of the present invention, the construction module 300 can be used to execute step 100 in the above method embodiment, the first determination module 302 can be used to execute step 102 in the above method embodiment, the second determination module 304 can be used to execute step 104 in the above method embodiment, and the acquisition module 306 can be used to execute step 106 in the above method embodiment.

[0122] In one embodiment of the present invention, the material temperature texture file of the area where the target to be calculated is located is obtained in the following manner:

[0123] Determine the environmental parameters corresponding to the area where the target to be calculated is located;

[0124] According to the environmental parameters, the energy balance equation corresponding to the second angle between the normal line of the surface element and the sunlight angle is constructed for the surface elements with different material characteristics in the area;

[0125] Solve each energy balance equation, calculate the surface element temperature value of the second angle corresponding to the surface element of different material characteristics, and obtain the material temperature texture file of the area where the target to be calculated is located.

[0126] In one embodiment of the present invention, the environmental parameters correspond to different date ranges, and the material temperature texture file of the area where the target to be calculated is located includes material temperature texture sub-files of different date ranges;

[0127] Based on the material temperature texture file generated in advance for the area where the target to be calculated is located, the temperature value of each face element to be calculated is determined, including:

[0128] In the material temperature texture file of the area where the target to be calculated is located, a material temperature texture subfile whose date range includes the time to be calculated is obtained, and the temperature value of each face element to be calculated is determined according to the material temperature texture subfile.

[0129] In one embodiment of the present invention, constructing an energy balance equation for face elements with different material characteristics at the second angle includes:

[0130] Establish the energy balance equation based on the one-dimensional heat transfer model;

[0131] The one-dimensional heat transfer model is:

[0132]

[0133] Where T is the temperature of the preset surface element, t is the time variable, x is the space variable, and α is the heat transfer coefficient;

[0134] The energy balance equation is:

[0135] E d +E sky +E ground +E cv +G h -Q it =0

[0136] Where: E d is the solar radiation energy absorbed by the preset surface element at the second angle (w / m 2 );E sky is the atmospheric radiation energy absorbed by the preset surface element (w / m 2 );E ground It is the ground radiation energy absorbed by the preset surface element, which is divided into two parts: ground reflection and ground radiation (w / m 2 );E cvis the energy absorbed by the preset surface element through convection heat exchange with the external air (w / m 2 ), which is mainly related to air convection velocity and atmospheric temperature; G h is the heat conduction flux between the outer surface of the preset surface element and the inner surface of the preset surface element (w / m 2 );Q it is the preset surface element radiation emittance (w / m 2 ), which is related to the self-emissivity and temperature of the preset surface element.

[0137] In one embodiment of the present invention, according to the material characteristics and the first angle of each face element to be calculated, based on a material temperature texture file pre-generated for the area where the target to be calculated is located, determining the temperature value of each face element to be calculated includes:

[0138] For each face element to be calculated, a target material temperature texture file having the same material characteristics as the face element to be calculated is obtained from the material temperature texture file;

[0139] Selecting a second angle according to the angle difference between each second angle contained in the target material temperature texture file and the first angle;

[0140] The temperature value of the bin to be calculated is determined according to the bin temperature value corresponding to the selected second angle.

[0141] In one embodiment of the present invention, selecting the second angle according to the angle difference between each second angle included in the target material temperature texture file and the first angle includes:

[0142] Determine whether the minimum angle difference among the angle differences is 0. If so, determine the second angle corresponding to the minimum angle difference as the selected second angle. Otherwise, determine the second angle with the smallest angle difference among the angle differences as the selected second angle.

[0143] In one embodiment of the present invention, determining the temperature value of the bin to be calculated according to the bin temperature value corresponding to the selected second angle includes:

[0144] If the minimum angle difference among the angle differences is 0, the temperature value of the bin corresponding to the selected second angle is determined as the temperature value of the bin to be calculated;

[0145] If the minimum angle difference among the angle differences is not 0, a linear interpolation calculation is performed based on the two selected second angles, the surface element temperature values ​​corresponding to the two second angles, and the first angle;

[0146] It is to be understood that the structure illustrated in the embodiment of the present invention does not constitute a specific limitation on a target temperature field calculation device. In other embodiments of the present invention, a target temperature field calculation device may include more or fewer components than shown in the figure, or combine some components, or split some components, or arrange the components differently. The components shown in the figure may be implemented in hardware, software, or a combination of software and hardware.

[0147] The information interaction, execution process and other contents between the modules in the above-mentioned device are based on the same concept as the embodiment of the method of the present invention. For the specific contents, please refer to the description in the embodiment of the method of the present invention, and no further description is given here.

[0148] An embodiment of the present invention further provides a computing device, including a memory and a processor, wherein the memory stores a computer program, and when the processor executes the computer program, a target temperature field calculation method in any embodiment of the present invention is implemented.

[0149] An embodiment of the present invention further provides a computer-readable storage medium, on which a computer program is stored. When the computer program is executed by a processor, the processor executes a target temperature field calculation method in any embodiment of the present invention.

[0150] Specifically, a system or device equipped with a storage medium can be provided, on which software program code that implements the functions of any of the above-mentioned embodiments is stored, and a computer (or CPU or MPU) of the system or device can be enabled to read and execute the program code stored in the storage medium.

[0151] In this case, the program code itself read from the storage medium can realize the function of any one of the above-mentioned embodiments, and thus the program code and the storage medium storing the program code constitute a part of the present invention.

[0152] The storage medium embodiments for providing the program code include a floppy disk, a hard disk, a magneto-optical disk, an optical disk (such as CD-ROM, CD-R, CD-RW, DVD-ROM, DVD-RAM, DVD-RW, DVD+RW), a magnetic tape, a non-volatile memory card, and a ROM. Alternatively, the program code can be downloaded from a server computer by a communication network.

[0153] In addition, it should be clear that the functions of any of the above embodiments can be implemented not only by executing the program code read by the computer, but also by enabling an operating system operating on the computer to complete part or all of the actual operations based on instructions from the program code.

[0154] In addition, it can be understood that the program code read from the storage medium is written to a memory provided in an expansion board inserted into the computer or to a memory provided in an expansion module connected to the computer, and then based on the instructions of the program code, a CPU installed on the expansion board or expansion module is enabled to perform part or all of the actual operations, thereby realizing the functions of any of the above-mentioned embodiments.

[0155] It should be noted that, in this article, relational terms such as first and second are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Moreover, the terms "include", "comprise" or any other variants thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or device. In the absence of further restrictions, the elements defined by the statement "comprise a ..." do not exclude the presence of other identical factors in the process, method, article or device including the elements.

[0156] A person of ordinary skill in the art can understand that all or part of the steps of implementing the above method embodiments can be completed by hardware related to program instructions, and the aforementioned program can be stored in a computer-readable storage medium. When the program is executed, it executes the steps of the above method embodiments; and the aforementioned storage medium includes: ROM, RAM, magnetic disk or optical disk, etc., various media that can store program codes.

[0157] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technical features therein. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for calculating a target temperature field, characterized in that: include: Constructing a finite element geometric network model of the target to be calculated, and dividing the target to be calculated into a plurality of surface elements to be calculated; Determining the material characteristics and position coordinates of each of the face elements to be calculated; For each of the face elements to be calculated, based on the position coordinates of the face element to be calculated, determining a first angle between the normal line of the face element to be calculated and the sunlight angle at the time to be calculated; According to the material characteristics of each of the face elements to be calculated and the first angle, based on the material temperature texture file pre-generated for the area where the target to be calculated is located, the temperature value of each face element to be calculated is determined to obtain the temperature field of the target to be calculated; wherein the material temperature texture file includes the face element temperature values ​​corresponding to the second angle between the face element normal and the solar illumination angle corresponding to different material characteristics of the target to be calculated; the second angle is obtained by traversing the solar azimuth and pitch angle according to a preset angle interval; The material temperature texture file of the area where the target to be calculated is located is obtained in the following way: Determine the environmental parameters corresponding to the area where the target to be calculated is located; According to the environmental parameters, constructing energy balance equations corresponding to the second angle between the normal of the surface element and the sunlight angle for the surface elements with different material characteristics in the area; Solving each of the energy balance equations, calculating the surface element temperature values ​​of the surface elements with different material characteristics corresponding to the second angle, and obtaining a material temperature texture file of the area where the target to be calculated is located; The determining the temperature value of each of the face elements to be calculated according to the material characteristics of each of the face elements to be calculated and the first angle, based on a material temperature texture file pre-generated for the area where the target to be calculated is located, comprises: For each of the face elements to be calculated, obtaining a target material temperature texture file having the same material characteristics as the face element to be calculated from the material temperature texture file; Selecting a second angle according to an angle difference between each second angle included in the target material temperature texture file and the first angle; The temperature value of the bin to be calculated is determined according to the bin temperature value corresponding to the selected second angle.

2. The method according to claim 1, characterized in that The environmental parameters correspond to different date ranges, and the material temperature texture file of the area where the target to be calculated is located contains material temperature texture sub-files of different date ranges; The step of determining the temperature value of each of the face elements to be calculated based on a material temperature texture file generated in advance for the area where the target to be calculated is located comprises: In the material temperature texture file of the area where the target to be calculated is located, a material temperature texture sub-file whose date range includes the time to be calculated is obtained, and the temperature value of each face element to be calculated is determined according to the material temperature texture sub-file.

3. The method according to claim 1, characterized in that Constructing the energy balance equation of the surface elements with different material characteristics at the second angle includes: Establishing the energy balance equation based on a one-dimensional heat transfer model; The one-dimensional heat transfer model is: In the formula, T is the temperature of the surface element, t is the time variable, x is the space variable, and α is the heat transfer coefficient; The energy balance equation is: E d +E sky +E ground +E cv +G h -Q it =0 Where: E d is the solar radiation energy absorbed by the surface element at the second angle; E sky is the atmospheric radiation energy absorbed by the surface element; E ground It is the ground radiation energy absorbed by the surface element, which is divided into two parts: ground reflection and ground radiation; E cv It is the energy absorbed by the surface element through convection heat exchange with the external air, which is mainly related to the air convection speed and atmospheric temperature; G h is the heat conduction flux between the outer surface of the element and the inner surface of the element; Q it It is the surface element radiation emittance, which is related to the surface element's own emissivity and temperature.

4. The method according to claim 1, characterized in that The selecting the second angle according to the angle difference between each second angle included in the target material temperature texture file and the first angle comprises: Determine whether the minimum angle difference among the angle differences is 0. If so, determine the second angle corresponding to the minimum angle difference as the selected second angle. Otherwise, determine the second angle with the minimum angle difference among the angle differences as the selected second angle.

5. The method according to claim 4, characterized in that The step of determining the temperature value of the bin to be calculated according to the bin temperature value corresponding to the selected second angle includes: If the minimum angle difference among the angle differences is 0, the temperature value of the bin corresponding to the selected second angle is determined as the temperature value of the bin to be calculated; If the minimum angle difference among the angle differences is not 0, a linear interpolation calculation is performed according to the two selected second angles, the two surface element temperature values ​​corresponding to the second angles, and the first angle; The bin temperature value obtained by linear interpolation calculation is determined as the temperature value of the bin to be calculated.

6. A target temperature field calculation device, characterized in that: include: A construction module is used to construct a finite element geometric network model of the target to be calculated, and divide the target to be calculated into a plurality of surface elements to be calculated; A first determination module is used to determine the material characteristics and position coordinates of each of the face elements to be calculated; A second determination module is used to determine, for each of the surface elements to be calculated, a first angle between the normal line of the surface element to be calculated and the sunlight angle at the time to be calculated based on the position coordinates of the surface element to be calculated; An acquisition module is used to determine the temperature value of each of the facets to be calculated according to the material characteristics of each facet to be calculated and the first angle, based on a material temperature texture file generated in advance for the area where the target to be calculated is located, and obtain the temperature field of the target to be calculated; wherein the material temperature texture file includes the facet temperature value corresponding to the second angle between the facet normal of the target to be calculated corresponding to different material characteristics and the solar illumination angle; the second angle is obtained by traversing the solar azimuth and pitch angle according to a preset angle interval; The material temperature texture file of the area where the target to be calculated is located is obtained in the following way: Determine the environmental parameters corresponding to the area where the target to be calculated is located; According to the environmental parameters, constructing energy balance equations corresponding to the second angle between the normal of the surface element and the sunlight angle for the surface elements with different material characteristics in the area; Solving each of the energy balance equations, calculating the surface element temperature values ​​of the surface elements with different material characteristics corresponding to the second angle, and obtaining a material temperature texture file of the area where the target to be calculated is located; The determining the temperature value of each of the face elements to be calculated according to the material characteristics of each of the face elements to be calculated and the first angle, based on a material temperature texture file pre-generated for the area where the target to be calculated is located, comprises: For each of the face elements to be calculated, obtaining a target material temperature texture file having the same material characteristics as the face element to be calculated from the material temperature texture file; Selecting a second angle according to an angle difference between each second angle included in the target material temperature texture file and the first angle; The temperature value of the bin to be calculated is determined according to the bin temperature value corresponding to the selected second angle.

7. A computing device, comprising a memory and a processor, wherein the memory stores a computer program, and when the processor executes the computer program, the method according to any one of claims 1 to 5 is implemented.

8. A computer-readable storage medium having a computer program stored thereon, which, when executed in a computer, causes the computer to execute the method according to any one of claims 1 to 5.

Images