CN113063524A - Method for constructing temperature field of high-temperature abnormal area of horizontal warehouse by using limited temperature measuring points - Google Patents

Method for constructing temperature field of high-temperature abnormal area of horizontal warehouse by using limited temperature measuring points Download PDF

Info

Publication number
CN113063524A
CN113063524A CN202110291498.1A CN202110291498A CN113063524A CN 113063524 A CN113063524 A CN 113063524A CN 202110291498 A CN202110291498 A CN 202110291498A CN 113063524 A CN113063524 A CN 113063524A
Authority
CN
China
Prior art keywords
temperature
sensor
ijk
heating source
grain
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
Application number
CN202110291498.1A
Other languages
Chinese (zh)
Other versions
CN113063524B (en
Inventor
刘晓鹏
张永林
刘羽祚
余宇骁
王立宗
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Jingzhou Yuzhong Food Machinery Co ltd
Wuhan Polytechnic University
Original Assignee
Jingzhou Yuzhong Food Machinery Co ltd
Wuhan Polytechnic University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Jingzhou Yuzhong Food Machinery Co ltd, Wuhan Polytechnic University filed Critical Jingzhou Yuzhong Food Machinery Co ltd
Priority to CN202110291498.1A priority Critical patent/CN113063524B/en
Publication of CN113063524A publication Critical patent/CN113063524A/en
Application granted granted Critical
Publication of CN113063524B publication Critical patent/CN113063524B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K13/00—Thermometers specially adapted for specific purposes
    • G01K13/10—Thermometers specially adapted for specific purposes for measuring temperature within piled or stacked materials
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K1/00—Details of thermometers not specially adapted for particular types of thermometer
    • G01K1/02—Means for indicating or recording specially adapted for thermometers
    • G01K1/022—Means for indicating or recording specially adapted for thermometers for recording
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K1/00—Details of thermometers not specially adapted for particular types of thermometer
    • G01K1/02—Means for indicating or recording specially adapted for thermometers
    • G01K1/026—Means for indicating or recording specially adapted for thermometers arrangements for monitoring a plurality of temperatures, e.g. by multiplexing
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K1/00—Details of thermometers not specially adapted for particular types of thermometer
    • G01K1/14—Supports; Fastening devices; Arrangements for mounting thermometers in particular locations
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08—Thermal analysis or thermal optimisation

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Arrangements For Transmission Of Measured Signals (AREA)

Abstract

The method for constructing the temperature field of the high-temperature abnormal area of the horizontal warehouse by using the limited temperature measuring points comprises the following steps: s1 setting a temperature measurement network in the horizontal warehouse; s2, determining a heating source point in the horizontal warehouse; s3, combining a temperature measurement network, and calculating according to the heating source point to obtain the attenuation rate; s4 calculating the layer height as z0Average temperature of grain in plane
Figure DDA0002982844670000011
S5 obtains a temperature field map. This method of utilizing limited temperature measurement point to found abnormal regional temperature field of flat storehouse high temperature, it sets up the sensor to each storehouse of grain storehouse equidistant in the flat storehouse, utilizes the limited temperature measurement point data in the storehouse and change, calculates and obtains the temperature field picture, can ventilate or turn over the storehouse to abnormal region to the temperature field picture and handle, realizes rapid cooling, guarantees the grain safety of storing up.

Description

Method for constructing temperature field of high-temperature abnormal area of horizontal warehouse by using limited temperature measuring points
[ technical field ] A method for producing a semiconductor device
The invention belongs to the technical field of grain storage, and particularly relates to a method for constructing a temperature field of a high-temperature abnormal area of a horizontal warehouse by using a limited temperature measuring point.
[ background of the invention ]
China is a large country for grain production, and grain storage is the fundamental guarantee for national safety and social stability. The grain storage in China mainly adopts a horizontal warehouse type, and measures such as moisture control, killing of grain insects, ventilation or turning over the warehouse for cooling need to be taken to maintain the grain in the warehouse in a proper storage environment. However, the grain seeds are used as organisms, and part of the seeds can generate heat due to biological activities of the grains, so that the temperature of a local grain stack is increased, the quality of stored grains is deteriorated, and the loss of the stored grains is caused. For a normal-temperature grain storage bin, the normal grain storage temperature is required to be within 28 ℃ (within 15 ℃ of a low-temperature bin), if the grain storage temperature is 28-30 ℃, a ventilation cooling measure needs to be taken for the heating grain stack, and if the grain storage temperature is higher than 30 ℃, the heating grain stack needs to be manually or mechanically turned over. Some thermometry networks are also presented nowadays, but the thermometry networks cannot directly acquire data to construct a high-temperature region temperature field.
[ summary of the invention ]
In order to solve the above problems, the present invention provides a method for constructing a temperature field in an abnormal high temperature area of a flat-type storehouse using a limited temperature measurement point, in which sensors are provided at equal intervals in each grain storehouse in the flat-type storehouse, a temperature field map is calculated and obtained using the limited temperature measurement point data and changes thereof in the storehouse, and the abnormal area can be ventilated or turned over with respect to the temperature field map, thereby realizing rapid cooling and ensuring grain storage safety.
The invention is realized by the following technical scheme, and provides a method for constructing a temperature field of a high-temperature abnormal area of a horizontal warehouse by using a limited temperature measuring point, which comprises the following steps:
s1 setting a temperature measurement network in the horizontal warehouse;
s2, determining a heating source point in the horizontal warehouse;
s3, combining a temperature measurement network, and calculating according to the heating source point to obtain the attenuation rate;
s4 calculating the layer height as z0Average temperature of grain in plane
Figure BDA0002982844650000023
S5 obtains a temperature field map.
Specifically, the step S1 is implemented according to the following method:
each grain storehouse with the temperature measurement network is located in the flat storehouse, the temperature measurement network includes four at least layers of sensors that set up from top to bottom, and every layer of sensor comprises a plurality of rows, a plurality of rows of sensors, and wherein the sensor interval between ranks is dx＝dyM, and the distance between two adjacent layers of sensors is dzThe unit is m, each sensor corresponds to a number, the numbers of the rows, the columns and the layers are i, j and k respectively, the physical positions of the sensors are expressed by a Cartesian coordinate system, and the sensors are temperature measurement sensors.
Specifically, the step S2 is implemented according to the following method:
sensor a in step S1ijkReporting abnormal state and highest temperature, wherein the temperature value is KijkDetermining the coordinates of the point location of the heat source as O (a, b, c) and the temperature value as K by measuring the highest temperature within the radius R of the neighborhood of the sensor0The calculation formula of the sensor neighborhood radius R is as follows:
Figure BDA0002982844650000021
specifically, the step S3 is implemented according to the following method:
s31 selecting the highest temperature sensor a according to the step S2ijkAnd the adjacent row, column and layer are numbered as ai-1jk、ai+1jk、aij-1k、aij+1k、aijk-1、aijk+1Calculating the temperature values reported by 7 sensors in total;
s32 linear attenuation of grain pile temperature along x, y and z directions with attenuation rate of rhox、ρy、ρzThe unit is ℃/m, and the temperature and the attenuation rate meet the following formula:
Figure BDA0002982844650000022
s33, establishing coordinate system Ox ' y ' z ' by the heating source point, overlooking the plane direction of Ox ' y ' to obtain the heating source plane, if the heating source O is positioned at the sensor a of the highest temperature sensorijkUpper region, then sensor aijkAnd a sensor ai+1jkIs positioned at the same side of the heating source, and a sensor ai+1jkBit phase ratio sensor ai-1jkThe distance from the heating source is farther, so the temperature of the heating source and the temperature of the heating source meet the following conditions: ki+1jk<Ki-1jk(ii) a If the heating source O is positioned at the sensor a of the highest temperature sensorijkLower region, sensor aijkAnd a sensor ai+1jkIs positioned at the same side of the heating source, and a sensor ai-1jkBit phase ratio sensor ai+1jkThe distance from the heating source is farther, so the temperature of the heating source and the temperature of the heating source meet the following conditions: ki-1jk<Ki+1jk(ii) a The sensor a is known from the geometrical relationship of the circumference division of the coordinate systemijkAdjacent sensor a at the same side of heating sourcei+1jkOr ai-1jkDistance d ofxIs approximately equally divided into line segments dxIs projected to the heat generation source plane through an isotherm to form a line segment d'xFrom the geometric relationship: dxAnd d'xApproximately equal, indicating sensor aijkAnd a sensor ai+1jkOr ai-1jkThe temperature between the two approximately satisfies the linear attenuation change rule in the x or x' direction, and is approximately equal to the attenuation rate rho of the heating sourcexEquality, the following equation is obtained:
Figure BDA0002982844650000031
in the formula (3), ρx0Is aijkAnd ai+1jk(or a)i-1jk) The decay rate therebetween;
s34 likewise, the sensor a can be obtainedijkAnd a sensor ai+1jkOr ai-1jkThe temperature between the two temperature sensors also satisfies a linear attenuation change rule in the y or y ', z or z' direction, and the attenuation rate along the y and z directions is as follows:
Figure BDA0002982844650000032
in the formula (4), ρy0Is aijkAnd aij+1kOr aij-1kAttenuation ratio between, pz0Is aijkAnd aijk+1Or aijk-1The decay rate therebetween.
Specifically, the step S4 is implemented according to the following method:
in the horizontal warehouse, under the condition of normal grain storage, a temperature gradient (figure 7) exists along the height direction (namely the z direction), and the normal temperature mean value of the 4-layer temperature measurement network is respectively calculated to be K1、K2、K3、K4The layer heights are respectively z1、z2、z3、z4From the above data, a function of temperature variation with height can be fitted as:
Figure BDA0002982844650000033
in equation (5), m ≠ n, according to function KzF (z) the layer height z can be calculated0The average temperature of the grain in the plane is
Figure BDA0002982844650000041
Specifically, the step S5 is implemented according to the following method:
s51 has a pair layer height of z0The grain in the plane is influenced by the heat source more and more as the grain is closer to the heat source O, the temperature of the grain is higher, and the heat source O respectively takes rhox、ρy、ρzThe speed of the heat source is linearly attenuated along the directions of x, y and z axes, so that the influence range of the heat source O on the plane is limited, and when the temperature of the heat source is attenuated to the average temperature of the grain on the layer
Figure BDA0002982844650000042
When the heat source is equal, the heat source does not conduct any more, and the influence area is an ellipse according to the formula (2), and the function expression of the ellipse area is as follows:
Figure BDA0002982844650000043
in the formula, LxThe length dimension of the grain pile along the x direction, m;
Lyis the length dimension of the grain pile along the y direction, m.
S52 setting layer height as z0Substituting the coordinates of the points in the plane into a formula (6), and if the calculation result H (x, y) is less than or equal to 0, calculating the temperature according to a formula (2); if the calculation result is (x, y)>0, calculating the temperature of the point according to the formula (5), and taking the average temperature of the layer as
Figure BDA0002982844650000044
S53 fitting the data obtained by S52 with Origin software to obtain the layer height z ═ z0Cloud plot of planar temperature field (FIG. 8), when z0The value range is [0, 6.2 ]]And then, temperature field cloud pictures of grains in different layer high planes of the horizontal warehouse can be respectively obtained.
The invention provides a method for constructing a temperature field of a high-temperature abnormal area of a horizontal warehouse by using a limited temperature measuring point, which is characterized in that a starting heat source point is determined by setting a sensor in the horizontal warehouse, temperature field cloud pictures of grains in different layer high planes of the horizontal warehouse are finally obtained by using Origin software fitting according to data obtained by calculation, and the abnormal area can be purposefully ventilated or turned over according to the obtained temperature field cloud pictures, so that the rapid cooling is realized, the grain storage safety is ensured, and the grain storage loss is avoided.
[ description of the drawings ]
FIG. 1 is a view showing a state of a flat house storehouse and its stored grain in a method of constructing a temperature field of an abnormal high temperature area of the flat house storehouse according to the invention using a limited temperature measuring point;
FIG. 2 is a front view showing a temperature measuring network arrangement of the flat house storehouse according to the method of the present invention for constructing a temperature field of an abnormal high temperature area of the flat house storehouse using a limited temperature measuring point;
FIG. 3 is a plan view showing a temperature measuring network arrangement of the flat house storehouse according to the method of the present invention for constructing a temperature field of an abnormal high temperature area of the flat house storehouse using a limited temperature measuring point;
FIG. 4 is a diagram showing a distribution of the positions of the sensors in the high-temperature area of the heat source in the method for constructing the temperature field in the high-temperature abnormal area of the horizontal warehouse by using the limited temperature measuring points according to the present invention;
FIG. 5 is a gradient diagram of the temperature field along the x direction in the method for constructing the temperature field in the high-temperature abnormal area of the horizontal warehouse by using the limited temperature measuring points according to the present invention;
FIG. 6 is a diagram showing a distribution of sensor positions near a heat source in a method for constructing a temperature field in a high-temperature abnormal area of a horizontal warehouse by using a limited temperature measuring point according to the present invention;
FIG. 7 is a diagram of the temperature gradient in the direction of the upper layer height in the method for constructing the temperature field in the high-temperature abnormal region of the horizontal warehouse by using the limited temperature measuring points according to the invention;
FIG. 8 is a schematic diagram of the calculation of the temperature field of the high temperature region in the method for constructing the temperature field of the high temperature abnormal region of the horizontal warehouse by using the limited temperature measuring points according to the present invention;
fig. 9 is a graph of temperature fields of high-temperature abnormal regions of different layers calculated according to an embodiment of the method for constructing a temperature field of a high-temperature abnormal region of a horizontal warehouse by using a limited temperature measuring point.
[ detailed description ] embodiments
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is further described in detail with reference to the following embodiments.
The invention provides a method for constructing a temperature field in a high-temperature abnormal area of a flat storehouse by using a limited temperature measuring point, referring to fig. 1-3, a temperature measuring network of the flat storehouse or the grain storehouse shown in fig. 1-3 is composed of 192 sensors in total of 6 rows, 8 columns and 4 layers, each temperature measuring sensor corresponds to a definite number (i, j and k are respectively a row number, a column number and a layer number), and the physical position of the temperature measuring network is expressed by a Cartesian coordinate system. In the temperature measurement network, the distance between the sensors in the rows and the columns is dx＝d y4m, the distance between two adjacent layers of sensors is dz＝1.8m。
The specific implementation steps are as follows:
(1) if sensor a in thermometric network432Reporting abnormal conditions and the highest temperature (the temperature value is 32 ℃), determining the coordinates of the point position of the heat source to be O (14.4, 8, 4.5) by measuring the highest temperature of the sensor with the neighborhood radius of R being within 3m, and recording the temperature value K at the position0The temperature was 34.8 ℃.
(2) Referring to FIG. 4, the highest temperature sensor a is selected432And the adjacent row, column and layer are numbered as a332、a532、a422、a442、a431、a433Temperature sensors of (1), the temperatures of which are respectively indicated as K332＝26.7℃、K532＝29.1℃、K422＝30.5℃、K442＝25.8℃、K431＝30.9℃、K43328.8 ℃, calculated according to formula (2) and formula (3):
Figure BDA0002982844650000061
(3) under the condition of normal storage of grains in the horizontal warehouse, a temperature gradient exists along the height direction (z direction), and the normal temperature mean value of the 4-layer temperature measurement network is calculated to be K according to the formula (2)1、K2、K3、K4The layer heights are respectively 0.3m, 2.1m, 3.9m and 5.7m, and the data can be fitted according to the data and the formula (5)The temperature as a function of height is:
Kz＝0.003z3-0.064z2+1.695z-18.297；
(4) this example is intended to calculate z1＝2.1m、z2＝3.3m、z3＝4.5m、z4The average normal temperature of each of the four layers is calculated from the formula (5) as a temperature field having a height of 5.7m layers
Figure BDA0002982844650000062
Figure BDA0002982844650000063
Will z1Substituting the coordinate data of each 2.1m layer into a function H (x, y) of formula (6), and calculating the temperature according to formula (2) if the result satisfies that H (x, y) is less than or equal to 0; if calculating H (x, y)>0, the point temperature is taken as the average temperature of the layer
Figure BDA0002982844650000064
Similarly, z can be calculated2＝3.3m、z3＝4.5m、z4A 5.7m layer high temperature field; the temperature field of the high temperature abnormal region with the layer height is fitted by Origin software as shown in FIG. 9.
According to the temperature field obtained by the calculation, the height of the measurement layer in the grain depot is z1＝2.1m、z2＝3.3m、z3＝4.5m、z4Temperature measurement is carried out in a high-temperature area of 5.7m, actual temperature measurement is carried out at 50 points sampled in the high-temperature area, and theoretical errors between the measured temperature and a fitted temperature field diagram are within 3 percent.

Claims (6)

1. The method for constructing the temperature field of the high-temperature abnormal area of the horizontal warehouse by using the limited temperature measuring points is characterized by comprising the following steps of:
s1 setting a temperature measurement network in the horizontal warehouse;
s2, determining a heating source point in the horizontal warehouse;
s3, combining a temperature measurement network, and calculating according to the heating source point to obtain the attenuation rate;
s4 calculating the layer height as z0Average temperature of grain in plane
Figure FDA0002982844640000011
S5 obtains a temperature field map.
2. The method for constructing the temperature field of the high-temperature abnormal area of the horizontal warehouse according to the claim 1, wherein the step S1 is implemented according to the following method:
each grain storehouse with the temperature measurement network is located in the flat storehouse, the temperature measurement network includes four at least layers of sensors that set up from top to bottom, and every layer of sensor comprises a plurality of rows, a plurality of rows of sensors, and wherein the sensor interval between ranks is dx＝dyM, and the distance between two adjacent layers of sensors is dzThe unit is m, each sensor corresponds to a number, the numbers of the rows, the columns and the layers are i, j and k respectively, the physical positions of the sensors are expressed by a Cartesian coordinate system, and the sensors are temperature measurement sensors.
3. The method for constructing the temperature field of the high-temperature abnormal area of the horizontal warehouse according to the claim 2, wherein the step S2 is implemented according to the following method:
sensor a in step S1ijkReporting abnormal state and highest temperature, wherein the temperature value is KijkDetermining the coordinates of the point location of the heat source as O (a, b, c) and the temperature value as K by measuring the highest temperature within the radius R of the neighborhood of the sensor0The calculation formula of the sensor neighborhood radius R is as follows:
Figure FDA0002982844640000012
4. the method for constructing the temperature field of the high-temperature abnormal area of the horizontal warehouse according to the claim 3, wherein the step S3 is implemented according to the following method:
s31 selecting the highest temperature sensor a according to the step S2ijkAnd the adjacent row, column and layer are numbered as ai-1jk、ai+1jk、aij-1k、aij+1k、aijk-1、aijk+1Calculating the temperature values reported by 7 sensors in total;
s32 linear attenuation of grain pile temperature along x, y and z directions with attenuation rate of rhox、ρy、ρzThe unit is ℃/m, and the temperature and the attenuation rate meet the following formula:
Figure FDA0002982844640000021
s33, establishing coordinate system Ox ' y ' z ' by the heating source point, overlooking the plane direction of Ox ' y ' to obtain the heating source plane, if the heating source O is positioned at the sensor a of the highest temperature sensorijkUpper region, then sensor aijkAnd a sensor ai+1jkIs positioned at the same side of the heating source, and a sensor ai+1jkBit phase ratio sensor ai-1jkThe distance from the heating source is farther, so the temperature of the heating source and the temperature of the heating source meet the following conditions: ki+1jk<Ki-1jk(ii) a If the heating source O is positioned at the sensor a of the highest temperature sensorijkLower region, sensor aijkAnd a sensor ai+1jkIs positioned at the same side of the heating source, and a sensor ai-1jkBit phase ratio sensor ai+1jkThe distance from the heating source is farther, so the temperature of the heating source and the temperature of the heating source meet the following conditions: ki-1jk<Ki+1jk(ii) a The sensor a is known from the geometrical relationship of the circumference division of the coordinate systemijkAdjacent sensor a at the same side of heating sourcei+1jkOr ai-1jkDistance d ofxIs approximately equally divided into line segments dxIs projected to the heat generation source plane through an isotherm to form a line segment d'xFrom the geometric relationship: dxAnd d'xApproximately equal, indicating sensor aijkAnd a sensor ai+1jkOr ai-1jkThe temperature between the two approximately satisfies the linear attenuation change rule in the x or x' direction, and is approximately equal to the attenuation rate rho of the heating sourcexEquality, the following equation is obtained:
Figure FDA0002982844640000022
in the formula (3), ρx0Is aijkAnd ai+1jkOr ai-1jkThe decay rate therebetween;
s34 likewise, the sensor a can be obtainedijkAnd a sensor ai+1jkOr ai-1jkThe temperature between the two temperature sensors also satisfies a linear attenuation change rule in the y or y ', z or z' direction, and the attenuation rate along the y and z directions is as follows:
Figure FDA0002982844640000031
in the formula (4), ρy0Is aijkAnd aij+1kOr aij-1kAttenuation ratio between, pz0Is aijkAnd aijk+1Or aijk-1The decay rate therebetween.
5. The method for constructing the temperature field of the high-temperature abnormal area of the horizontal warehouse according to the claim 4, wherein the step S4 is implemented according to the following method:
in the horizontal warehouse, under the condition of normal grain storage, a temperature gradient exists along the height direction, and the normal temperature mean value of the 4-layer temperature measurement network is respectively calculated to be K1、K2、K3、K4The layer heights are respectively z1、z2、z3、z4From the above data, a function of temperature variation with height can be fitted as:
Figure FDA0002982844640000032
in equation (5), m ≠ n, according to function KzF (z) the layer height z can be calculated0The average temperature of the grain in the plane is
Figure FDA0002982844640000033
6. The method for constructing the temperature field of the high-temperature abnormal area of the horizontal warehouse according to the claim 5, wherein the step S5 is implemented according to the following method:
s51 has a pair layer height of z0The grain in the plane is influenced by the heat source more and more as the grain is closer to the heat source O, the temperature of the grain is higher, and the heat source O respectively takes rhox、ρy、ρzThe speed of the heat source is linearly attenuated along the directions of x, y and z axes, so that the influence range of the heat source O on the plane is limited, and when the temperature of the heat source is attenuated to the average temperature of the grain on the layer
Figure FDA0002982844640000034
When the heat source is equal, the heat source does not conduct any more, and the influence area is an ellipse according to the formula (2), and the function expression of the ellipse area is as follows:
Figure FDA0002982844640000041
in the formula (6), LxIs the length dimension of the grain pile along the x direction and has the unit of m and LyThe length dimension of the grain pile along the y direction is m;
s52 setting layer height as z0Substituting the coordinates of the points in the plane into a formula (6), and if the calculation result H (x, y) is less than or equal to 0, calculating the temperature according to a formula (2); if the calculation result is (x, y)>0, calculating the temperature of the point according to the formula (5), and taking the average temperature of the layer as
Figure FDA0002982844640000042
S53 fitting the data obtained by S52 with Origin software to obtain the layer height z0Temperature field of plane, when z0The value range is [0, 6.2 ]]And then, temperature field cloud pictures of grains in different layer high planes of the horizontal warehouse can be respectively obtained.
CN202110291498.1A 2021-03-18 2021-03-18 Method for constructing temperature field of high-temperature abnormal region of horizontal warehouse by using limited temperature measuring points Active CN113063524B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110291498.1A CN113063524B (en) 2021-03-18 2021-03-18 Method for constructing temperature field of high-temperature abnormal region of horizontal warehouse by using limited temperature measuring points

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110291498.1A CN113063524B (en) 2021-03-18 2021-03-18 Method for constructing temperature field of high-temperature abnormal region of horizontal warehouse by using limited temperature measuring points

Publications (2)

Publication Number Publication Date
CN113063524A true CN113063524A (en) 2021-07-02
CN113063524B CN113063524B (en) 2023-10-10

Family

ID=76562007

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110291498.1A Active CN113063524B (en) 2021-03-18 2021-03-18 Method for constructing temperature field of high-temperature abnormal region of horizontal warehouse by using limited temperature measuring points

Country Status (1)

Country Link
CN (1) CN113063524B (en)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6278083B1 (en) * 2016-09-07 2018-02-14 マツダ株式会社 Thermal conductivity measurement method and apparatus
CN108805501A (en) * 2018-06-07 2018-11-13 吉林大学 A kind of granary storage state monitoring and managing method based on correlation statistics
CN108896197A (en) * 2018-05-22 2018-11-27 安徽瑞鑫自动化仪表有限公司 A kind of intelligent thermometric regulator control system under hot environment
CN109060873A (en) * 2018-03-16 2018-12-21 河南工业大学 A kind of grain matter heat transmitting and condense parameter detection device and method
CN110309610A (en) * 2019-07-08 2019-10-08 重庆科技学院 A kind of provision store reconstructing temperature field method based on improved adaptive GA-IAGA
CN110631639A (en) * 2019-10-24 2019-12-31 富德康(北京)科技股份有限公司 Method for judging grain safety according to grain temperature and humidity
CN110704512A (en) * 2019-10-22 2020-01-17 吉林大学 Granary ventilation time interval judgment method based on historical grain situation data
CN112067162A (en) * 2020-09-04 2020-12-11 李灵芝 Granary temperature distribution mode detection method based on spatial interpolation and BIM

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6278083B1 (en) * 2016-09-07 2018-02-14 マツダ株式会社 Thermal conductivity measurement method and apparatus
CN109060873A (en) * 2018-03-16 2018-12-21 河南工业大学 A kind of grain matter heat transmitting and condense parameter detection device and method
CN108896197A (en) * 2018-05-22 2018-11-27 安徽瑞鑫自动化仪表有限公司 A kind of intelligent thermometric regulator control system under hot environment
CN108805501A (en) * 2018-06-07 2018-11-13 吉林大学 A kind of granary storage state monitoring and managing method based on correlation statistics
CN110309610A (en) * 2019-07-08 2019-10-08 重庆科技学院 A kind of provision store reconstructing temperature field method based on improved adaptive GA-IAGA
CN110704512A (en) * 2019-10-22 2020-01-17 吉林大学 Granary ventilation time interval judgment method based on historical grain situation data
CN110631639A (en) * 2019-10-24 2019-12-31 富德康(北京)科技股份有限公司 Method for judging grain safety according to grain temperature and humidity
CN112067162A (en) * 2020-09-04 2020-12-11 李灵芝 Granary temperature distribution mode detection method based on spatial interpolation and BIM

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
王志山,范明,毛广卿,陈正权,邢九菊: "绘制粮食仓库温度分布图,快速确定局部操作作业点", 粮食与饲料工业 *
程龙;曹阳;李光涛;张来林;: "圆球导热法测定粮食导热系数研究", 中国粮油学报 *
闫艳霞;曹玲芝;: "粮仓谷物局部发热温度场数学模型研究与应用", 粮食与食品工业 *
陈桂香;刘超赛;蒋敏敏;陈家豪;王海涛;张宏伟;: "储料竖向压力对粮仓中小麦粮堆湿热传递的影响", 农业工程学报 *

Also Published As

Publication number Publication date
CN113063524B (en) 2023-10-10

Similar Documents

Publication Publication Date Title
CN111427100B (en) Typhoon center positioning method and device and typhoon path generation method
US11317562B2 (en) Generating a yield map for an agricultural field using classification and regression methods
KR20140007298A (en) Apparatus and method for detecting temperature distribution
WO2014098698A1 (en) Sensor calibration method, computer program and computer readable medium
CN108764802A (en) A kind of ABC strategy process suitable for the supervision of grain storage number
Oke The temperature profile near the ground on calm clear nights
CN107687899A (en) A kind of infrared measurement of temperature method and system
CN108362311A (en) A kind of distance measuring sensor segmented compensated distance method
CN108981960A (en) A kind of AD strategy process suitable for grain storage supervision
CN107133286A (en) A kind of three-dimensional map generalization of computer room temperature parameter distribution and analysis method and system
CN110810027B (en) Grain temperature-based circulation temperature control method in granary
CN107424947A (en) The temperature testing method and system of hull cell process equipment
CN113063524A (en) Method for constructing temperature field of high-temperature abnormal area of horizontal warehouse by using limited temperature measuring points
CN111637974A (en) Non-contact temperature measuring device and temperature measuring method
CN112362693A (en) Method and system for calculating evapotranspiration amount based on soil heat flux estimation
CN108764749A (en) A kind of 3I strategy process for preferably depositing evaluation suitable for grain storage
CN113008413A (en) Method for accurately reversely calculating target point of heating grain by using existing temperature measurement network
CN111579086A (en) Remote infrared temperature measurement precision correction method based on distance compensation
CN105259209B (en) A kind of measuring system and computational methods of LNG evaporation rates
CN108549068B (en) Three-dimensional scanning data processing method and data processing system
CN110631721B (en) Granary heat insulation judgment method based on grain condition big data
CN113190633A (en) Dual-distance-based climate space partitioning method
CN111947631A (en) Triangular coordinate transfer measurement method applied to inside of tunnel
CN111551259A (en) Temperature measurement correction method based on cooperation of UWB and infrared sensor
CN106204525B (en) The detection unified approach of tree crown Light interception-rate

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant