Disclosure of Invention
The invention provides a bridge temperature measuring method, a bridge temperature measuring device, bridge temperature measuring equipment and a readable storage medium, and aims to solve the problems in the prior art.
In order to achieve the purpose, the invention is realized by the following technical scheme:
in a first aspect, the present invention provides a bridge temperature measurement method, including:
determining a target time sequence corresponding to a bridge to be detected, wherein the target time sequence is used for representing solar weather data in the environment where the bridge to be detected is located;
determining the initial temperature of the bridge to be detected and the radiation temperature of the bridge to be detected based on the target time sequence;
determining the temperature of a first surface fluid of the bridge to be detected according to a preset first condition;
determining the temperature of a second surface fluid of the bridge to be detected according to a preset second condition, wherein the second surface is a surface opposite to the first surface;
determining the temperature of the bridge to be tested based on the initial temperature, the radiation temperature, the first surface fluid temperature and the second surface fluid temperature.
Optionally, the target time sequence includes a solar radiation sequence, an atmospheric temperature sequence, a wind speed sequence, and a structural internal air temperature sequence of the bridge to be tested.
Optionally, the determining the initial temperature and the radiation temperature of the bridge to be measured based on the target time series includes:
determining the initial atmospheric temperature in the atmospheric temperature sequence as the initial temperature of the bridge to be measured;
determining the radiation temperature of the bridge to be detected according to the solar radiation sequence and the atmospheric temperature sequence;
recording the radiation temperature of the bridge to be measured as T e Said T is e The following calculation relationship is satisfied:
in the formula, T a,s Is a sequence of atmospheric temperatures, q s Is a solar radiation sequence, gamma is the radiance of the surface of the concrete structure, and h is the concreteComprehensive heat transfer coefficient of the surface of the structure, T e,s Is a radiant temperature sequence.
Optionally, the first condition satisfies the following calculation relationship:
wherein x represents the position of the outer surface of the bridge concrete member, T (x, T) represents the temperature value to be obtained of the bridge concrete member, h
1 Represents the value of the overall heat transfer coefficient of the outer surface of the bridge concrete member,
t (x, T) is shown to calculate the partial derivative of the variable x, and k represents the heat conductivity coefficient of the bridge concrete member;
the first surface fluid temperature satisfies the following calculation relationship:
f 1 (t)=T 1,v -T 1,w cos(w 1 (t-δ 1 )),t∈(t 1 ~t 2 );
in the formula (f) 1 (T) denotes the first surface fluid temperature, T 1,v Expressed as a function uniform value, T 1,w Representing the fluctuation value of the function, w 1 Representing the frequency, δ, of the function 1 Representing the phase of the function, t representing the current time, t 1 Indicates the start time, t 2 Indicating the end time.
Optionally, the second condition satisfies the following calculation relationship:
in the formula: x is the position of the inner surface of the bridge structure, h 2 The value of the comprehensive heat exchange coefficient of the inner surface of the bridge structure is obtained;
the second surface fluid temperature satisfies the following calculation relationship:
f 2 (t)=T 2,v -T 2,w cos(w 2 (t-δ 2 )),t∈(t 1 ~t 2 );
in the formula: t is 2,v For a gas temperature sample T in a bridge structure b Average air temperature of (T) 2,w Is the fluctuating air temperature in the cavity, w 2 Is the frequency of the cosine function, delta 2 Is the phase of the cosine function.
Optionally, when the temperature of the bridge to be measured is determined based on the initial temperature, the radiation temperature, the first surface fluid temperature, and the second surface fluid temperature, the following calculation relationship is satisfied:
in the formula (f)
i,n (t) a form of variation of the temperature of the fluid at the internal and external surfaces of the bridge structure, C
i,n 、
β
n 、X(β
n X) are all calculated coefficients, T
0 Alpha is the thermal diffusivity of the bridge concrete member, which is the initial temperature of the bridge structure.
In a second aspect, an embodiment of the present application further provides a bridge temperature measuring device, including:
the system comprises a first processing unit, a second processing unit and a third processing unit, wherein the first processing unit is used for determining a target time sequence corresponding to a bridge to be detected, and the second target time sequence is used for representing solar weather data in the environment where the bridge to be detected is located;
the second processing unit is used for determining the initial temperature of the bridge to be detected and the radiation temperature of the bridge to be detected based on the target time sequence;
the third processing unit is used for determining the temperature of the first surface fluid of the bridge to be detected according to a preset first condition;
the fourth processing unit is used for determining the temperature of a second surface fluid of the bridge to be detected according to a preset second condition, wherein the second surface is a surface opposite to the first surface;
and the fifth processing unit is used for determining the temperature of the bridge to be measured based on the initial temperature, the radiation temperature, the first surface fluid temperature and the second surface fluid temperature.
In a third aspect, an embodiment of the present application further provides an electronic device, which includes a processor, a memory, and a program or an instruction stored on the memory and executable on the processor, where the program or the instruction, when executed by the processor, implements the steps of the bridge temperature measurement method according to the first aspect.
In a fourth aspect, embodiments of the present application further provide a readable storage medium, on which a program or instructions are stored, and when executed by a processor, the program or instructions implement the steps of the bridge temperature measurement method according to the first aspect.
Has the advantages that:
the invention provides a bridge temperature measuring method, a device, equipment and a readable storage medium, wherein the bridge temperature measuring method comprises the steps of determining a target time sequence corresponding to a bridge to be measured, wherein the target time sequence is used for representing solar weather data in the environment where the bridge to be measured is located; determining the initial temperature of the bridge to be detected and the radiation temperature of the bridge to be detected based on the target time sequence; determining the temperature of a first surface fluid of a bridge to be measured according to a preset first condition; determining the temperature of a second surface fluid of the bridge to be detected according to a preset second condition, wherein the second surface is a surface opposite to the first surface; and determining the temperature of the bridge to be measured based on the initial temperature, the radiation temperature, the first surface fluid temperature and the second surface fluid temperature. Therefore, the temperature field distribution of different structures and different construction stages of the bridge can be rapidly calculated by combining the solar weather data in the environment where the bridge is located and the internal temperature of the bridge structure.
Detailed Description
The technical solutions of the present invention are described clearly and completely below, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
As shown in fig. 1, the present invention provides a bridge temperature measurement method, including:
step 101, determining a target time sequence corresponding to the bridge to be detected, wherein the target time sequence is used for representing solar weather data in the environment where the bridge to be detected is located.
In this embodiment, the target time series includes, but is not limited to, a solar radiation series, an atmospheric temperature series, a wind speed series, and a structural internal air temperature series of the bridge under test. The present invention is not limited to the above embodiments, and other types of parameters may be selected as the target time sequence in other possible embodiments, and only the temperature of the bridge needs to be reflected, but the embodiments of the present invention are within the scope of protection of the present invention regardless of the change of the parameters.
In a possible implementation, when acquiring the target time sequence, a small-sized weather station including a solar radiation meter, an anemometer and two temperature sensors may be installed at the bridge to be measured; the solar radiation meter is used for measuring a solar radiation sequence, the anemometer is used for measuring a wind speed sequence, and the temperature sensor is used for measuring an atmospheric temperature sequence and can also be used for measuring a structure internal air temperature sequence. The sequence refers to a period of measurement values, in other words, the sequence is a set including a plurality of measurement values, and details are not described herein. The measurement condition can be satisfied by that the small-sized meteorological station is positioned in an open field at the bridge site, is not shielded by surrounding structures, and is erected on the ground by a special meteorological support. The temperature sensor of the inside temperature of bridge structures is placed in the meteorological shutter and placed inside the bridge structures, and the temperature sensor is connected with the cables of the meteorological station. This is by way of example only and not by way of limitation.
In practice, when a target time sequence corresponding to a bridge to be measured is determined, carrying out statistical calculation on samples within a solar irradiation time range, specifically, determining an average value of atmospheric temperature and an average value of wind speed; the time series samples of the weather data are then used to determine the initial temperature and the radiant temperature of the bridge structure, for example, the atmospheric temperature at the time of rising sun may be used as the initial temperature of the bridge structure, and the radiant temperature of the bridge structure may be calculated from the time-by-time sample data, in a manner to be described in detail in the following steps.
Step 102, determining the initial temperature of the bridge to be measured and the radiation temperature of the bridge to be measured based on the target time sequence;
103, determining the temperature of a first surface fluid of the bridge to be detected according to a preset first condition;
104, determining the temperature of a second surface fluid of the bridge to be detected according to a preset second condition, wherein the second surface is a surface opposite to the first surface;
and 105, determining the temperature of the bridge to be measured based on the initial temperature, the radiation temperature, the first surface fluid temperature and the second surface fluid temperature.
The bridge temperature measuring method comprises the steps of firstly determining the initial temperature of a bridge to be measured and the radiation temperature of the bridge to be measured through a target time sequence, and then determining the first surface fluid temperature of the bridge to be measured according to a preset first condition; determining the temperature of a second surface fluid of the bridge to be detected according to a preset second condition, wherein the second surface is a surface opposite to the first surface; and determining the temperature of the bridge to be measured based on the initial temperature, the radiation temperature, the first surface fluid temperature and the second surface fluid temperature. The solar weather data in the environment where the bridge is located and the internal temperature of the bridge structure can be combined, and the temperature field distribution of different structures and different construction stages of the bridge can be rapidly calculated.
It should be noted that when determining the target time sequence corresponding to the bridge to be measured, it is necessary to satisfy the weather conditions of cloudless and breeze on a sunny day, and in this embodiment, the sunrise time, for example, 7 am, is taken as the first time, the sunset time, for example, 7 pm, is taken as the second time, and the target time sequence corresponding to the bridge to be measured within the time range of the first time and the second time is determined. Wherein, the collection frequency is once per hour, and the number of samples is m = (t) 2 -t 1 ) The intensity q of solar radiation and the atmospheric temperature T are formed a And temperature T in cavity of bridge structure b And time series samples of wind speed v.
Then, the time series samples of the atmospheric temperature are calculated statistically, and the average value T of the atmospheric temperature in one day is determined av The calculation formula satisfies the following relation:
s=1,2...m;
in the formula, T a,s Time-by-time samples for atmospheric temperature, T a =[T a,1 ,T a,2 …T a,m-1 ,T a,m ]M, m being the number of samples and s being the time series of samples.
Carrying out statistical calculation on the obtained time series samples of the wind speed, and determining the average value v of the wind speed a The calculation formula satisfies the following relation:
s=1,2...m;
in the formula, v a,s For wind speed time-by-time samples, v = [ v ] a,1 ,v a,2 …v a,m-1 ,v a,m ],s=1,2...m。
Further, the obtained solar radiation q is subjected to s Atmospheric temperature T a The time series samples are subjected to statistical calculation to determine the radiation air temperature T of the outer surface of the bridge structure e Time series samples of (a);
s=1,2...m
in the formula, T e,s For radiating a temperature sequence, T e =[T e,1 ,T e,2 …T e,m-1 ,T e,m ],s=1,2...m。
For the radiation air temperature T obtained in the above step e In the calculation formula, a value taking method of the parameter gamma is explained, wherein the gamma is the radiance of the surface of the concrete structure, the value taking range is 0.4-0.7, and the judgment is carried out according to the coating color of the surface, the cleanness of the surface and the roughness of the surface. For example, the darker the coating color of the structure surface, the larger the coefficient value; the neater the structure surface is, the smaller the coefficient value is; the rougher the structure surface, the larger the coefficient value.
Wherein the obtained radiation air temperature T c In the calculation formula of e, h is the comprehensive heat exchange coefficient of the surface of the concrete structure and the comprehensive heat exchange coefficient h of the outer surface of the concrete structure 1 Value 15.5, comprehensive heat exchange coefficient h of comprehensive heat exchange coefficient of inner surface of concrete structure 2 Value 12 is taken.
In this embodiment, the initial atmospheric temperature in the atmospheric temperature sequence is determined as the initial temperature of the bridge to be measured, for example, the atmospheric temperature T a =[T a,1 ,T a,2 …T a,m-1 ,T a,m ]Determining the initial temperature T of the bridge structure 0 Value, T 0 =T a,1 I.e. the first atmospheric temperature collected at the sunrise time of the sun is the initial temperature of the concrete element.
Then, determining a fluid temperature function formula of the outer surface of the bridge structure: the radiant air temperature T of the bridge structure obtained in the step e Fitting time series samples T for the fluid temperature of the outer surface of the bridge by using Fourier function e The second order Fourier function is as follows.
f 1 (t)=T 1,v +acos(w 1 t)+bsin(w 1 t),t∈(t 1 ~t 2 );
In the formula, T 1,v 、a、b、w 1 Fitting parameters of a Fourier function;
the above equation is converted to the form of a cosine function as follows:
f 1 (t)=T 1,v +T 1,w cos(w 1 (t-δ 1 )),t∈(t 1 ~t 2 );
in the formula, T
1,w For uniform values of the cosine function, using the formula
Calculating; delta
1 As a function of phase, using formula
Calculating; t is t
1 、t
2 Sunrise and sunset times, respectively.
It should be noted that, because the outer surface and the inner surface of the bridge structure are both in a composite heat exchange state, the heat flow boundary is nonlinear, and the equation cannot be effectively used for solving. Therefore, it is necessary to determine a convection boundary condition (a preset first condition) of the outer surface of the bridge structure and a convection boundary condition (a preset second condition) of the inner surface of the bridge structure.
In this embodiment, the functional formula f for obtaining the boundary of the outer surface (first surface) of the bridge structure is described above 1 (t), determining that the convection boundary condition of the outer surface of the bridge structure meets the following formula:
where x =0 is the outer surface position of the bridge concrete member and T (x, T) is the temperature value to be determined for the bridge concrete member, h in this embodiment 1 Taking the value of 15.5. The value of k is a thermal conductivity coefficient of the concrete member of the bridge, and the value is determined by concrete materials of actual engineering, and generally ranges from 1.2 to 1.8.
Further, the fluid temperature function of the inner surface of the bridge structure is determined as follows:
to obtain an internal temperature sequence sample T of the bridge structure b For the fluid temperature at the inner surface of the bridge, the transformation is expressed in the form of a cosine function as follows:
f 2 (t)=T 2,v -T 2,w cos(w 2 (t-δ 2 )),t∈(t 1 ~t 2 );
in the formula, T
2,v Is a sequence sample T
b Average value of (1) using
Calculating to obtain; t is a unit of
2,w Is a sequence sample T
b Is taken as the fluctuation value of->
Calculating to obtain; t is
b,max As sequence samples T
b The highest air temperature, i.e. T
b,max =max[T
b,1 ,T
b,2 …T
b,m-1 、T
b,m ];T
b,min Is a sequence sample T
b Of the lowest air temperature, i.e. T
b,min =min[T
b,1 ,T
b,2 …T
b,m-1 、T
b,m ];w
2 For frequencies based on a cosine function, use->
m is the number of samples in a day; delta
2 The phase of the middle cosine function is obtained by taking the time when the highest air temperature appears on the surface air temperature of the bridge structure, for example, when the time when the highest air temperature appears is 15
2 =15。
Determining the convection boundary conditions of the inner surface of the bridge structure as follows:
wherein x = d is the position of the inner surface of the bridge structure, T is the temperature field to be evaluated of the bridge structure, h 2 Is the comprehensive heat exchange coefficient value h of the inner surface of the bridge structure 2 Value 12 is taken. This is by way of example only and not by way of limitation, and the values may be adjusted within a certain range in other possible embodiments.
Further, according to the fluid temperature function formulas of the inner surface and the outer surface of the bridge structure, the fluid temperature function formula of the surface of the bridge structure obtained in the step is converted as follows:
when i =1, the temperature function of the fluid on the outer surface of the bridge structure is shown, and accordingly, f 1 (t)=ΔT f,1 +T f,1 (t)。
When i =2, the fluid temperature function of the inner surface of the bridge structure is indicated, and accordingly, f 2 (t)=ΔT f,2 +T f,2 (t)。
Wherein:
ΔT f,i =T i,v -T 0 ;
T f,i (t)=-T i,w cos(w i (t-δ i ));
in the formula, T i,v Is the mean value of the sequence samples, w i I =1,2 for the frequency of the functional formula. In the formula of this embodiment, when i is 1, the formula indicates a fluid temperature function formula of the outer surface of the bridge structure, i.e., i =1. When i is 2, the functional formula of the fluid temperature of the outer surface and the inner surface of the bridge structure is shown, i.e. i =2.
Establishing a relation model of a bridge structure temperature field: the relation model for predicting the bridge structure temperature field obtained according to the steps is as follows:
(4)θ(x,t)=T(x,t)-T 0 ;
in the formula, α and k are the thermal diffusivity and the thermal conductivity of concrete respectively, and θ (x, t) is the temperature of the bridge structure to be obtained, is a variable and is related to the position x inside the bridge structure and the time t. d is the position of the inner surface of the bridge structure and x is the position inside the bridge structure.
The above-described relational model is a thermal conduction partial differential equation in heat transfer science, and in this embodiment, by establishing the above-described relational model, the equation can be made solvable. The temperature of the bridge can be rapidly determined.
Specifically, a value taking method of the parameters alpha and k is explained, the specific value is determined by the thermophysical property parameter of the concrete material, and if field test conditions exist, the value can be taken in the thermophysical property test of the concrete material; if the field test condition is not present, values of α =0.0032 and k =1.8 are commonly used.
Then, based on the mathematical model of the temperature field of the bridge structure, the calculation formula of the temperature field of the bridge structure obtained by performing integral calculation on the mathematical model of the temperature field satisfies the following relation:
in the formula (f) i,n (t) is a variation of the fluid temperature at the inner and outer surfaces of the bridge structure, as follows:
i=1,2;
n=1,2,3,4,5;
wherein n represents the nth characteristic value,
for a calculation coefficient, it can be calculated by the following formula:
C i,n the calculation coefficient for one of the formulas can be calculated according to the following formula:
in the formula, h i Is the comprehensive heat exchange coefficient, beta, of the external (i = 1) or internal (i = 2) surface of the bridge structure n As a characteristic value, X (. Beta.) n x i ) As a characteristic function, N (. Beta.) n ) Is the eigenvalue.
Herein C
i,n It can be composed of several expressions:
2X(β
n x
i )、3N(β
n )
the calculation coefficient for one of the formulas can be calculated according to the following formula:
β n the calculation coefficient for one of the formulas can be calculated according to the following formula:
wherein L is the thickness of the bridge concrete member in meters/m.
X(β n x i ) The calculation coefficient for one of the formulas can be calculated according to the following formula:
X(β n ,x i )=β n cosβ n x i +H 1 sinβ n x i ;
in the formula, x i Is the external or internal surface of a bridge structure, H 1 Is the ratio of the comprehensive heat exchange coefficient and the heat conductivity coefficient of the outer surface of the bridge structure, i =1, x i =0; when i =2, x i =L。
N(β n ) The calculation coefficient for one of the formulas can be calculated according to the following formula:
the calculation formula for obtaining the bridge structure temperature field is a series superposition formula about a time variable t and a space variable x, can effectively solve the problems of distribution and prediction of the bridge structure temperature field in solar radiation, and can be conveniently applied to a high-speed railway cable-stayed bridge with large span and a complex system: and combining the actually measured meteorological data and the internal temperature of the bridge structure, and rapidly calculating the temperature field distribution of different structures and different construction stages. For example, the temperature field distribution and the temperature form of the concrete bridge deck in the construction stage, the box girder-track system in the operation stage or different members (concrete bridge pier, bridge tower and track slab) of the cable-stayed bridge at the same time are adopted.
It should be noted that the bridge temperature measurement method provided in the embodiment of the present application may be applied to, but is not limited to, cable-stayed bridge structures, which may be concrete track slabs, concrete piers, and concrete pylons.
The invention discloses a bridge temperature measuring method which has certain universality, can predict a bridge structure temperature field only by meteorological data of local regions and the internal temperature of a bridge structure at a bridge site, and simultaneously meets the prediction of the bridge structure temperature at different construction stages. The method is simple and convenient to implement and has good feasibility.
The embodiment of the present application further provides a bridge temperature measuring device, including:
the system comprises a first processing unit, a second processing unit and a third processing unit, wherein the first processing unit is used for determining a target time sequence corresponding to a bridge to be detected, and the second target time sequence is used for representing solar weather data in the environment where the bridge to be detected is located;
the second processing unit is used for determining the initial temperature of the bridge to be detected and the radiation temperature of the bridge to be detected based on the target time sequence;
the third processing unit is used for determining the temperature of the first surface fluid of the bridge to be detected according to a preset first condition;
the fourth processing unit is used for determining the temperature of a second surface fluid of the bridge to be detected according to a preset second condition, wherein the second surface is a surface opposite to the first surface;
and the fifth processing unit is used for determining the temperature of the bridge to be measured based on the initial temperature, the radiation temperature, the first surface fluid temperature and the second surface fluid temperature.
The bridge temperature measuring device can realize the steps of the bridge temperature measuring method embodiment, and can obtain the same or similar beneficial effects, and the details are not repeated here.
Embodiments of the present application further provide an electronic device, which includes a processor, a memory, and a program or an instruction stored in the memory and executable on the processor, where the program or the instruction when executed by the processor implements the steps of the bridge temperature measurement method described above.
Embodiments of the present application further provide a readable storage medium, on which a program or instructions are stored, and when executed by a processor, the program or instructions implement the steps of the bridge temperature measurement method as described above.
The foregoing detailed description of the preferred embodiments of the invention has been presented. It should be understood that numerous modifications and variations could be devised by those skilled in the art in light of the present teachings without departing from the inventive concepts. Therefore, the technical solutions available to those skilled in the art through logic analysis, reasoning and limited experiments based on the prior art according to the concept of the present invention should be within the scope of protection defined by the claims.