Abstract
Using the monitoring temperature data from a large-span steel-concrete composite girder bridge, the spatial distribution laws of the transversal and vertical temperature differences in the steel-concrete composite girder are analyzed; then, the probabilistic statistical characteristics of the daily extreme values of temperature differences are analyzed using probability statistical method, and then cycle-life simulation is carried out to obtain the yearly maximum or minimum values in the whole bridge service life; furthermore, the standard values of positive and negative temperature differences are calculated, and finally the correlation between the standard values and the monitoring extreme values are explored using correlation analysis method. The results show that: (1) there are significant transversal and vertical temperature differences in the steel-concrete composite girder, and the maximum positive temperature difference in the transversal direction of this girder reaches 7.40°C, and the minimum negative temperature difference in the vertical direction reaches −6.80°C; (2) the generalized extreme value distribution (GEVD) has a smaller fitting error 0.0016 compared to the extreme value distribution (EVD) and the normal distribution (ND), so GEVD is selected as the optimal probability density function; (3) the maximum and minimum values of standard values for transversal temperature differences are 10.74°C and −6.38°C, respectively; the maximum and minimum values of standard values for vertical temperature differences are 4.81°C and −7.77°C, respectively; (4) there is a good linear correlation between the measured values and the standard values, with a linear correlation coefficient of 0.9044; (5) there is a big difference between current bridge codes and the calculation results by comparison, hence the actual standard values of temperature difference for steel-concrete composite girder bridges is relatively complex, and the calculation results can provide an important reference for the current bridge codes.
Keywords
Introduction
For large-span bridge structures in the whole service life, they are affected by environmental factors such as solar radiation, daily temperature differences and climate change, which may cause severe problems such as cracking and fatigue damage to bridge structures (Han et al., 2022; Huang et al., 2018, 2023; Jing et al., 2024; Li et al., 2023; Wang et al., 2021; Yang et al., 2024; Zhou et al., 2024). Therefore, the temperature field of bridge structures has a significant impact on structure damage, making it an indispensable factor in the load design of bridge structures (Keles et al., 2024; Li et al., 2023; Li et al., 2023, 2024; Shan et al., 2023). For example, Huang et al. (2018) pointed out that the changes in structural dynamic parameters resulting from damage may be masked by the change induced from variations in environmental effects, where the effect of temperature changes is deemed most important, so in structural health monitoring (SHM) of civil infrastructure systems it is of great importance to eliminate the effects of environmental changes, especially those of temperature, from the damage detection process. Accurate calculation of temperature action is of great importance for safety assessment of bridge structures.
There are many large-span bridge types, and among these bridge types the steel-concrete composite girder bridges are widely constructed in recent years because of their outstanding advantages in material utilization, economic benefit, construction convenience, durability, and seismic performance. Currently, the temperature action model of steel-concrete composite girders has been specified by the bridge design codes home and abroad (Ministry of Transport of the People’s Republic of China, 2020; American Association of State Highway and Transportation Officials, 2017; Brussels: European Committee of Standardization, 2003). For example, China and USA codes specify that the temperature action of the steel components in steel-concrete composite girders is uniform, without consideration of temperature gradients (Ministry of Transport of the People’s Republic of China, 2020; American Association of State Highway and Transportation Officials, 2017). Hence, current research on the temperature field of steel-concrete composite girders commonly treat the temperature field of the steel components as a uniform temperature field, without considering temperature gradients (Keles et al., 2024; Nikola et al., 2023; Wang et al., 2019, 2023; Zhu et al., 2021).
Some new research results based on the latest monitoring data show that there are temperature gradients in the steel components of steel-concrete composite girders (Huang et al., 2024a; Meng et al., 2023; Song et al., 2023; Wang et al., 2024; Xia et al., 2020; Zhou et al., 2020; Zhu et al., 2023, 2024). For example, Meng et al. (2023) analyzed a continuous steel-concrete composite railway girder bridge and found that the temperature gradient effect of the steel components has a more adverse effect on the local stress compared with the uniform temperature effect; Song et al. (2023) analyzed the temperature field data of a high-pier long-span steel-concrete composite girder bridge and found that there are certain temperature gradients in the steel components in addition to the daily and yearly temperature changes; Wang et al. (2024) analyzed the monitoring temperature field data of a steel-concrete composite girder bridge and found that the steel components has noticeable transversal and vertical temperature differences, and the transversal positive temperature difference reaches 14.75°C and the vertical negative temperature difference reaches −14.53°C; Zhu et al. (2024) analyzed the monitoring temperature field data of a railway suspension bridge with steel-concrete composite girder and found that the temperature field of the steel components have obvious non-uniform characteristics due to the influence of sun shadowing.
However, current research results mainly focuses on the monitoring period, which cannot fully reflect the most unfavorable conditions throughout the entire service life, so the analysis results cannot be used as standard values for thermal action design. According to current bridge codes (Ministry of Transport of the People’s Republic of China, 2020; American Association of State Highway and Transportation Officials, 2017; Brussels: European Committee of Standardization, 2003), the standard value is the most unfavorable temperature difference value in the whole bridge service life, and yearly maximum or minimum values in the whole bridge service life should be used for calculation of standard values, but the monitoring data is difficult to cover the whole bridge service life.
Therefore, this research puts forward a calculation method of standard values of temperature differences through probabilistic statistics and cycle-life simulation analysis using the monitoring data of a large-span steel-concrete composite girder bridge. First, the spatial distribution patterns of the transversal and vertical temperature differences in steel-concrete composite girder bridges are revealed; second, the probability statistical characteristics of the daily positive and negative temperature differences are studied and their cycle-life simulation in 100 years is carried out; finally, the standard values of positive and negative temperature differences are calculated and then compared with current bridge codes home and abroad. The research results have important reference significance for temperature action design.
Bridge health monitoring system
To study the spatial distribution characteristics of the temperature field of large-span steel-concrete composite girder bridges, the monitoring data of the temperature field at Hongwan Waterway Main Channel Bridge were collected using the bridge health monitoring system to analyze the temporal variation characteristics of the temperature field of the steel-concrete composite girder bridge.
The bridge is a five-span steel-concrete composite cable-stayed girder bridge in a northwest-southeast direction, which is located in Xiangzhou District, Zhuhai City. It starts from the connection line of the Hong Kong-Zhuhai-Macao Bridge and ends at the intersection of the Jiangmen-Zhuhai Expressway and the Zhuhai-Hege Expressway. The steel components are composed of rectangular and double-side I-shaped cross-sections, and the steel components is connected to the overlying concrete bridge deck through shear keys, forming a steel-concrete composite girder structure, as shown in Figure 1. To monitor the temporal variation characteristics of the steel-concrete composite girder, five temperature sensors are arranged at the top of the steel-concrete composite girder, denoted as S1 to S5; in addition, three temperature sensors are arranged at the bottom of steel-concrete composite girder, denoted as S6 to S8. The temperature sensors have a sampling frequency of 1 Hz and continuously collect data for 1 year. Temperature sensor arrangement.
The long-term monitoring law of temperature field of steel-concrete composite girder
Long-term monitoring rules of temperature values
The temperature values collected by the temperature sensor at position S
i
is represented by T
Time history of monitoring temperatures of TS4. Daily trend of monitoring temperatures of TS4.

Long-term monitoring rules of temperature difference values
Extreme values of temperature differences of hongwan waterway bridge (°C).
With regard to the maximum transversal temperature difference TS3,S4 and the maximum vertical temperature difference TS1,S6, their time histories from September 1st, 2023 to August 31st, 2024 are further plotted in Figure 4(a) and (b). It can be seen that the time histories show a uniform random characteristic without obvious seasonal variation features; the transversal temperature difference at TS3,S4 is significantly higher in positive values than in negative values, and the vertical temperature difference at TS1,S6 is significantly higher in negative values than in positive values. This further confirms the existence of significant transversal and vertical temperature differences between measurement points, suggesting that temperature differences should be considered in the design of steel-concrete composite girder bridges. Time histories of temperature differences (a) TS3,S4 and (b) TS1,S6.
Probabilistic statistical characteristics and cycle-life simulation analysis
The monitoring period of measured temperature values is relatively short, which cannot fully reflect the most unfavorable conditions throughout the entire service life, so the analysis of probabilistic statistical characteristics and cycle-life simulation is necessarily carried out to evaluate the temperature values in the whole bridge service life.
Probabilistic statistical characteristics analysis
According to the monitoring rules of temperature difference values above, the time histories show a uniform random characteristic, which is suitable to be used for probabilistic statistical characteristics analysis. The key point is to determine the probability density function of temperature difference values. The temperature difference values are divided into positive temperature differences and negative temperature differences. The daily maximum and minimum temperature differences are selected for analysis, corresponding to positive and negative temperature differences respectively. Several possible probability density functions are discussed and compared to determine the best probability density function to fit the probability density value with least errors.
The possible probability density functions are Normal Distribution Function (ND), Generalized Extreme Value Distribution Function (GEVD) and Extreme Value Distribution Function (EVD), and their equations are written as follows:
With regard to the parameter values of the three probability density functions, the maximum likelihood estimation method is used to obtain the optimal values. For convenience, the MATLAB Probability Density Fitting Toolbox can be used for quick solutions. In order to compare the fitting effect of the three probability density functions to determine the best probability density function with least fitting errors, the fitting errors error are calculated by the measured cumulative probability values vmeasured minus the fitted cumulative probability values vfitted as follows:
According to the method above, the probability density values of daily maximum and minimum temperature differences (denoted by TSi,Sj,maximum and TSi,Sj,minimum respectively) including transversal temperature differences TS1,S2, TS2,S3, TS3,S4, TS4,S5, TS6,S7, TS7,S8 and vertical temperature differences TS1,S6, TS3,S7, TS5,S8 are fitted by the three probability density functions ND, GEVD and EVD, as shown in Figure 5. The fitted probability density functions of daily extreme temperature differences (a) TS1,S2,maximum, (b) TS1,S2,minimum, (c) TS2,S3,maximum, (d) TS2,S3,minimum, (e) TS2,S4,maximum, (f) TS3,S4,minimum, (g) TS4,S5,maximum, (h) TS4,S5,minimum, (i) TS6,S7,maximum, (j) TS6,S7,minimum, (k) TS7,S8,maximum, (l) TS7,S8,minimum, (m) TS1,S6,maximum, (n) TS1,S6,minimum, (o) TS3,S7,maximum, (p) TS3,S7,minimum, (q) TS5,S8,maximum, and (r) TS5,S8,minimum.
The fitting errors of three probability density functions.
Cycle-life simulation analysis
One important part of temperature field analysis is to calculate the standard values of temperature differences for thermal effect analysis. According to current bridge codes (Ministry of Transport of the People’s Republic of China, 2020; American Association of State Highway and Transportation Officials, 2017; Brussels: European Committee of Standardization, 2003), the standard value is the most unfavorable temperature difference value in the whole bridge service life with the exceedance probability 2%. What should be mentioned is that during calculation of standard values the yearly maximum or minimum values in the whole bridge service life need to be known in advance, but the monitoring period for the temperature data is actually not long enough (e.g. the monitoring period of this bridge is only 1 year), so cycle-life simulation analysis is necessary to obtain the the yearly maximum or minimum values in the whole bridge service life.
In this research, a uniform random sampling method is employed to simulate the yearly maximum or minimum values in 100 years for this bridge by virtual of cumulative distribution functions of GEVD above. In detail, the daily maximum or minimum values in 100 years are firstly simulated using the cumulative distribution functions of GEVD above, and then the yearly maximum or minimum values are selected from the daily maximum or minimum values.
Hence, the key step is how to simulated the daily maximum or minimum values in 100 years. First, the number of simulated daily maximum or minimum values in 100 years is calculated, which is 36,500; then, uniform random sampling in the interval (0,1) is carried out 36,500 times to obtain the sampling values; third, each sampling value is treated as the cumulative probability value p
i
, so the corresponding daily maximum or minimum values TS,i are calculated as follows: Cycle-life simulation results (a) TS1,S2,maximum, (b) TS1,S2,minimum, (c)TS2,S3,maximum, (d)TS2,S3,minimum, (e)TS3,S4,maximum, (f)TS3,S4,minimum, (g) TS4,S5,maximum, (h) TS4,S5,minimum, (i) TS6,S7,maximum, (j) TS6,S7,minimum, (k) TS7,S8,maximum, (l) TS7,S8,minimum, (m) TS1,S6,maximum, (n) TS1,S6,minimum, (o)TS3,S7,maximum, (p) TS3,S7,minimum,(q) TS5,S8,maximum, and (r)TS5,S8,minimum.

Standard values of temperature differences of steel-concrete composite girders
Calculation method of standard values
According to current bridge codes (Ministry of Transport of the People’s Republic of China, 2020; American Association of State Highway and Transportation Officials, 2017; Brussels: European Committee of Standardization, 2003), the yearly maximum and minimum temperature difference values Tmax and Tmin should not exceed the positive and negative standard values Sp and Sn in the whole bridge service life, respectively:
According to the cycle-life simulation results above, the number of Tmax and Tmin is 100, so these Tmax and Tmin should satisfy the following inequalities:
Consider that the values of Tmax and Tmin in each year are random and independent which follow the same probability density function, so equations (10a) and (10b) are further written as follows:
With regard to fp (Tmax) or fp (Tmin), 100 yearly maximum or minimum values are firstly selected from the 100-year cycle-life simulated temperature difference data; then, the probability density values are calculated using through probabilistic and statistical analysis; finally, GEVD is used to fit the probability density values to determine the parameter values as shown in Figure 7. What should be mentioned is that GEVD is selected as the probability density function of Tmax or Tmin; the reason is that if the best probability density function of daily maximum or minimum values is GEVD as analyzed above, the best probability density function of yearly maximum or minimum values is also GEVD according to the probabilistic and statistical theory. The fitted probability density functions of yearly extreme temperature differences (a) TS1,S2,maximum, (b) TS1,S2,minimum, (c) TS2,S3,maximum,(d) TS2,S3,minimum, (e) TS3,S4,maximum, (f) TS3,S4,minimum, (g) TS4,S5,maximum, (h) TS4,S5,minimum, (i) TS6,S7,maximum, (j) TS6,S7,minimum, (k)TS7,S8,maximum, (l)TS7,S8,minimum, (m) TS1,S6,maximum, (n)TS1,S6,minimum, (o) TS3,S7,maximum(p) TS3,S7,minimum, (q) TS5,S8,maximum, and (r) TS5,S8,minimum.
Calculation results
Standard values TN and TP of temperature differences (°C).

The spatial distribution of the standard values of temperature differences.
It can be seen that the maximum and minimum values of standard values for transversal temperature differences are 10.74°C and –6.38°C, respectively; the maximum and minimum values of standard values for vertical temperature differences are 4.81°C and –7.77°C, respectively. Currently, the China and USA bridge codes treat the standard values as 0°C, and the European bridge code treats the standard values as ±15°C, which are different from the calculation results by comparison (Ministry of Transport of the People’s Republic of China, 2020; American Association of State Highway and Transportation Officials, 2017; Brussels: European Committee of Standardization, 2003). Hence, the actual standard values of temperature difference for steel-concrete composite girder bridges is relatively complex, and the calculation results can provide an important reference for the current bridge codes.
Furthermore, there is a one-to-one correspondence between the monitoring values and the standard values for temperature differences. For example, the maximum value of monitoring positive temperature difference for TS1,S2 is 4.91°C and its corresponding standard value is 5.42°C. Hence, the correlation scatter plot between the monitoring values and the standard values of the temperature differences for the Hongwan Waterway Bridge is shown in Figure 9. It can be seen that there is a good linear correlation between the measured values and the standard values, with a linear correlation coefficient of 0.9044. A linear function is used to fit the linear correlation which is expressed as S = 1.18 M + 0.39, where M denotes the monitoring value and S denotes the standard value. The correlation between measured extreme values and standard values.
Conclusion
Using the monitoring data of a large-span steel-concrete composite girder bridge, this research puts forward a calculation method of standard values of temperature differences in the whole service life is put forward by combination of probabilistic statistical characteristics analysis and life-cycle simulation analysis. The main conclusions are drawn as follows: (1) The monitoring results show that there are significant temperature differences between measuring points, where the maximum positive transversal temperature difference between measuring points S3 and S4 can reach 7.40°C, and the maximum negative vertical temperature difference between measuring points S1 and S6 can reach –6.80°C. This is different from current bridge codes which only consider uniform temperature action on the steel part of the steel-concrete composite bridge; (2) ND, GEVD and EVD are discussed and compared to determine the best probability density function to fit the probability density value with least errors. GEVD has the smallest average fitting errors 0.0016 compared with the other probability density functions, so in this study GEVD is selected as the optimal probability density functions of positive or negative temperature differences. Furthermore, cycle-life simulation method is put forward to obtain the yearly maximum or minimum values in the whole bridge service life using GEVD; (3) The maximum and minimum values of standard values for transversal temperature differences are 10.74°C and –6.38°C, respectively; the maximum and minimum values of standard values for vertical temperature differences are 4.81°C and –7.77°C, respectively. There is a good linear correlation between the measured values and the standard values, with a linear correlation coefficient of 0.9044, so a linear function is used to fit the linear correlation which is expressed as S = 1.18 M + 0.39; (4) Currently, the China and USA bridge codes treat the standard values as 0°C, and the European bridge code treats the standard values as ±15°C, which are different from the calculation results by comparison. Hence, the actual standard values of temperature difference for steel-concrete composite girder bridges is relatively complex, and the calculation results can provide an important reference for the current bridge codes.
Future work
This research put forwards a calculation method of standard values of temperature differences in the whole service life, which can provide a reference for thermal action design to evaluate the temperature effect of large-span steel-concrete composite girder bridges. What should be mentioned is that the monitoring results of Ts1,s6, Ts3,s7, and Ts5,s8 reveals the big temperature differences, but the number of temperature sensors may be not enough to reflect the whole vertical temperature gradients (Huang et al., 2024b). Hence, in the future more temperature sensors will be installed to obtain a more detailed temperature gradient.
Footnotes
Funding
China Railway Group Limited Science and Technology Research and Development Program (2022 - Key - 44).
Declaration of conflicting interests
The author declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Data Availability Statements
The data that support the findings of this study are available from the corresponding author upon reasonable request.
