Abstract
Due to the sheltering effect of the upper structures and the significant influence of temperature changes in the cables, the steel truss cable-stayed bridge with two-layer decks exhibits complex spatial thermal effects under solar radiation and environmental temperature. In order to study the temperature distribution characteristics and thermal deformation during the construction stage, the thermal effect of a steel truss cable-stayed bridge with two-layer decks under construction is analyzed through numerical simulation. Firstly, a temperature field analysis algorithm, which considers geo-meteorological factors and the sheltering relationship between members, was used to simulate the structure's transient temperature field and thermal deformation. Then, the temperature distribution characteristics and thermal deformation under the influence of seasons were discussed. The results indicated that the web members and the lower deck showed a significant temperature difference when shaded by the upper structure, with a maximum temperature difference of 7.6°C; The transverse temperature difference of the lower chord was the main factor controlling the transverse deformation of the steel truss girder; During a day in summer, the maximum vertical, longitudinal, and transverse deformations of the main members at the closure during the maximum single cantilever stage were 7.2 cm, 6.4 cm, and 5.9 cm, respectively.
Keywords
Introduction
Steel conducts heat quickly and is more sensitive to changes in environmental temperature. Under the influence of solar radiation and environmental temperature changes, the steel bridge will undergo significant temperature changes. (Chen, 2008; Liu et al., 2016; Wang et al., 2017). The structural temperature variations and inhomogeneities may significantly alter the static and dynamic performance of the bridge, which increases the bridge design and construction difficulty (Roeder, 2003; Xia et al., 2012). When the bridge is completed and opened to traffic, the negative impact of temperature changes on the bridge is even greater than the vehicle load (Catbas et al., 2008; Priestley, 1976; Tong et al., 2001; Xia et al., 2017).
Currently, the research methods for thermal effects of bridge structures mainly include numerical simulation and experimental testing (Li et al., 2023; Zhou and Yi, 2013). Compared to numerical simulation, the data obtained from experimental testing is more authentic and reliable. Therefore, some scholars use this method to study the bridge's temperature distribution and temperature response. Based on the measured temperature data of a constructed composite box-girder bridge for about 20 months, Chang and Im (2000) defined major thermal loading parameters that characterize the temperature profile, and detailed the seasonal behavior of these parameters. Roberts-Wollman et al. (2002) monitored temperature on a concrete box girder bridge for two and a half years, and presented equations to predict positive temperature differentials. Xu et al. (2010) analyzed the temperature characteristics of the Tsing Ma Bridge using years of monitoring data. Ding et al. (2012) proposed a critical temperature difference models in cross-section for thermal stress calculation through relevant analysis of long-term monitoring data of flat steel box girder bridge. As for the steel truss cable-stayed bridge with two-layer decks, Wang et al. (2021) conducted corresponding research. The long-term monitoring data of the Zhengzhou Yellow River Bridge revealed a significant temperature difference between the steel trusses. However, there are still some problems in studying the thermal effects of bridges through actual measurements. Firstly, there are few cross-sectional measurement points arranged, and the obtained temperature data is difficult to directly apply to the design calculation of bridges; Secondly, it is hard to obtain the temperature effect of the bridge structure. Therefore, using precise numerical simulation methods is still an important supplement to grasp the spatial temperature field and temperature effects of complex bridges.
Many scholars have conducted the finite element simulation of the bridges with simple structures in the form of calculating thermal boundary conditions. Emanual and Hulsey (1978) calculated the temperature variation over time of a highway composite beam bridge. Chen (2008) conducted a long-term finite element simulation of the temperature field of composite beams for 45 years and obtained sufficient temperature action samples for extreme value analysis. Zhou et al. (2016b) determined the initial thermal conditions of the model through pre-analysis, and then carried out transient heat transfer analysis to investigate the vertical and transverse temperature differences of the steel box girder. For rapid calculation of temperature loads, Fan et al. (2022) deduced and established vertical discrete and dimensionality-reduced model of a steel–concrete composite bridge. At the same time, outdoor sunlight experiments were conducted to study the temperature field. In addition, for the first time, Shan et al. (2023) integrated heat transfer analysis and field monitoring data, and established a 3D finite element model of a navigation channel bridge of the Hong Kong‒Zhuhai‒Macao Bridge, comprehensively and accurately studying the global 3D temperature distribution of the long-span cable-stayed bridge. Although many studies have been conducted, there is still a lack of research on the numerical simulation of long-span steel truss bridge. The reason for this is the difficulty of dynamically identifying shaded areas due to bridge components and terrain obstruction (Liu et al., 2019). For this purpose, Yin et al. (2014) compiled an ANSYS subroutine using ray tracing method to accurately simulate the main beam's covering effect on the concrete arch. bridge's arch. box. Zhu and Meng (2017) proposed a three-dimensional (3D) sunlight-sheltering algorithm based on a sun ray tracing method and applied the substructure method to thermo-mechanical coupling analysis, accurately predicting the temperature field and temperature effects of cable-stayed bridges. Zhu et al. (2023) calculated the sunlight shadow coefficient between different members using the bounding box algorithm and the elimination algorithm to simulate the non-uniform temperature field caused by occlusion between components.
In summary, for the temperature field of simple structure bridge, its numerical simulation technology has been more complete. However, for the steel truss bridge with two-layer decks, a long-span bridge with complex structures, the solution of the exact temperature field is still lacking. The problem is that the dynamic identification of shadow areas caused by bridge members is difficult. On the one hand, in solving the non-uniform temperature field caused by sunlight, it is difficult to identify the shadow areas caused by bridge members dynamically. This article calculates the thermal boundary conditions of a steel truss cable-stayed bridge with two-layer decks based on geo-meteorological factors, and employs a three-dimensional ray occlusion algorithm to identify its shadow area dynamically. During this period, the mutual sheltering effect of the upper deck on the web member, the upper deck and the web member on the lower deck is considered. Finally, the simulation of the transient temperature is achieved. On the other hand, when solving the thermal deformation of a steel truss cable-stayed bridge with two-layer decks, only the influence of uniform temperature field is generally considered. Therefore, based on the non-uniform temperature field simulated above and considering the temperature changes of the cables, this article establishes a thermo-mechanical coupling model to analyze the thermal deformation at the closure during the maximum single cantilever construction stage.
Thermal analysis method for bridge structures
Thermal environment of steel truss
The temperature distribution of the bridge structure under sunlight is basically determined by both internal heat transfer and external heat exchange. The internal heat conduction is controlled by the Fourier heat transfer equation. As for the conditions for the solution of the three-dimensional transient heat transfer partial differential equation, instantaneous temperature and boundary conditions are included (Lienhard and Lienhard, 2003). As shown in Figure 1, the heat exchange between the bridge and the external natural environment under sunlight includes three forms: solar radiation, convective heat transfer with the atmospheric environment, and radiation heat transfer on the structural surface. Heat exchange diagram with external environment: (a) The spatial relationship between the bridge and sun; (b) Heat exchange between the main girder and the external environment.
Heat-Transfer Analysis
Heat-Transfer theory
The internal heat transfer of the bridge structure follows Fourier's law, and the three-dimensional transient heat transfer partial differential equation is presented in equation (1):
Thermal Boundary Conditions
The thermal boundary conditions are generally divided into three types (Lienhard and Lienhard, 2003). Based on the known ambient temperature
During the process of sunlight passing through the atmosphere and reaching the earth's surface, it is continuously affected by various atmospheric components, such as absorption, reflection, and scattering. This is called solar radiation, which consists of direct, diffuse, and reflected solar radiation.
Direct solar radiation
Attenuated by the atmosphere, the solar radiation reaching the structure's surface is direct solar radiation
If the surface is inclined, the direct solar radiation can be expressed as follows:
Diffuse solar radiation
Scattered by air, aerosol molecules, and other particles, the solar radiation reaching the structure's surface is diffuse solar radiation
If the surface is inclined, the diffuse solar radiation can be expressed as follows:
Reflected solar radiation
The portion of direct and scattered solar radiation, which reaches the ground (water) surface and is reflected on the bridge's surface, is called reflected solar radiation
If the surface is inclined, the reflected solar radiation can be expressed as follows:
Overall, the total solar radiation heat flux density that the surface of the bridge structure can absorb is shown as follows:
Differing from solar radiation, which is short-wave radiation, the essence of radiation heat transfer is electromagnetic radiation generated by heat, belonging to long-wave radiation. It absorbs radiation from the atmosphere and surface and releases radiation toward the surrounding environment. It includes atmospheric radiation, surface environmental radiation, and diffuse surface radiation.
Atmospheric radiation
The atmosphere radiates energy outward by its temperature, and the part returned to the ground (water) surface is called atmospheric reverse radiation
If the surface is inclined, the atmospheric radiation can be expressed as follows:
Surface environmental radiation
The radiation caused by the exchange of energy between the ground (water) surface and the bridge surface is called surface environmental radiation
If the surface is inclined, the atmospheric radiation can be expressed as follows:
Diffuse surface radiation
The radiation emitted by the bridge surface in the form of electromagnetic waves is called diffuse surface radiation
Thus, the radiation heat flux density
The heat flow of Convective heat transfer is calculated by Newton's heat transfer equation:
Three-dimensional ray sheltering theory
Due to the motion of the sun, some surfaces of the bridge will be in shadow, and the thermal boundary conditions of that surface will also be affected. As shown in Figure 2, the conditions for generating shadows contain the following three types. Firstly, if the connection between the bridge surface and the sun passes through its own structure at any time, then the surface is in a permanent sheltering state. Secondly, if there is some time when the solar incident angle is greater than 90°, then the surface is in a self-sheltering state. Thirdly, if the bridge does not belong to the above two sheltering states, but is unable to be directly radiated by the sun due to the obstruction of other surfaces, then the surface is in a mutual sheltering state. The specific algorithm implementation process is given in Figure 3 (Meng, 2019). Schematic diagram of three types of sheltering. Flowchart of the 3D sheltering algorithm.

Initial temperature filed
When simulating the temperature field, the selection of the initial temperature field must be considered. Commonly, it is assumed that the initial temperature of the bridge is uniform. The ambient temperature at 3:00 before sunrise is employed as the initial value for the numerical simulation of the bridge, and this temperature is regarded as the lowest value. Then, through several days of cyclic calculation, errors caused by the initial value selection are eliminated. (Kim et al., 2015; Moorty and Roeder, 1992). This study calculated the temperature changes of a steel truss cable-stayed bridge with two-layer decks and found that the temperature showed periodic changes starting from the third day. The variation in air temperature is expressed as follows:
Numerical analysis and validation
Numerical analysis process
The analysis process of the thermal effect of the bridge mainly includes two parts, the solution of transient temperature field and thermal deformation. The transient temperature field is solved using the DFLUX and FILM subroutine modules in the finite element software ABAQUS. The subroutine simulates the process of solar motion based on the geo-meteorological information of the area where the bridge is located. It also simulates its radiation intensity in view of thermal boundary conditions and 3D sheltering algorithm. For thermal deformation analysis, a thermo-mechanical coupling model is established according to the temperature field obtained from the previous calculation.
Validation
Thermophysical Parameters of Experimental Modal.
The temperature distribution of the section where the temperature measuring points are located is shown in Figure 4. Here, the shadow appeared at the lower deck because of the obstruction of the upper deck and web members. The selected temperature measuring points included points 1 and 2 on the lower chord, and measuring points 3 and 4 on the lower deck. Take the average value of points 1 and 2 as the actual measuring result of the lower chord, and the average value of points 3 and 4 as the actual measuring result of the lower deck. The measuring result and the numerical result are plotted in Figure 5. Through comparison, the trend of the two changes was the same, and the data was in good agreement. The maximum error was 1.8°C, which meets the accuracy requirements of actual engineering. Temperature distribution and shadow region. Measured and simulated temperatures of the bridge for July 22,2022. (a) Lower chord;(b) Lower deck.

Case
Description of experimental bridge
The thermal effect analysis was conducted on an ongoing steel truss cable-stayed bridge with two-layer decks, spanning the Yellow River in Jinan, China. The bridge has a 1028m-long (90 + 180+488 + 180 + 90) steel truss girder, a 37.5m-wide upper deck and a 15m-wide lower deck. It is located at 116.1° E, 39.2° N, and has an azimuth angle of 8.8° northwest. Figure 6 shows the basic information of the experimental bridge location. The bridge crosses the Yellow River, and the surrounding terrain is flat. No mountains or high-rise buildings are blocking the sunlight. Experimental bridge.
The main girder of this bridge is composed by two-layer decks. As shown in Figure 7, the main girder consists of the upper chord (M2 and M3), lower chord (M5 and M6), side girder (M1 and M4), web members, inclined strut, upper deck system, and lower deck system. For the convenience of subsequent discussions, Figure 7 also represents the sun’s exposure to the steel truss girder at sunrise and sunset. Cross-section of the main girder.
Finite element simulation
The finite element software ABAQUS was employed to study the transient temperature field and thermal deformation. The flow chart is presented in Figure 8. The primary analysis step was to use DFLUX and FILM subroutine interfaces for the secondary development of finite element software ABAQUS, during which a three-dimensional occlusion algorithm was used to identify shadow areas and accurately simulate the heat exchange between the structure surface and the external environment. Flow chart for temperature field simulation.
Figure 9 shows the established finite element model of a steel truss cable-stayed bridge with two-layer decks. In the heat transfer analysis, Shell (DS4) element was utilized for the steel truss girder, and Solid (DC3C8) was utilized for the bridge tower. When calculating thermal deformation, the elements used for steel, concrete, and cables were Shell (S4R), Solid (C3D8R), and Truss (T3D2), respectively. Tables 2 and 3 list the thermophysical and mechanical parameters. Finite element model. Thermophysical Parameters. Mechanical Parameters.
Several issues needed to be noted in the heat transfer analysis of the steel truss cable-stayed bridge with two-layer decks. Firstly, steel has good thermal conductivity, and the plates are relatively thin, so the temperature difference in the thickness direction of the plates was generally neglected. Secondly, the upper and lower chords have a closed box-shaped cross-section, and the internal air velocity is approximately 0, so the convective heat transfer between the internal surface and the air was not considered.
In order to study the typical seasonal temperature behavior of the bridge, four representative weather conditions at the bridge's location were selected as the research objects. The selected dates were March 21 (spring equinox), June 22 (summer solstice), September 22 (autumn equinox), and December 23 (winter solstice), 2022. According to the information provided by the China Meteorological Data Network, the temperatures on that day were 9°C to 15°C, 22°C to 35°C, 19°C to 26°C, −1°C to 8°C, and the average wind speeds were 4.1 m/s, 3 m/s, 1.9 m/s, and 2.5 m/s, respectively.
Temperature distribution characteristics of the steel truss girder
This section analyzed the vertical temperature gradient, transverse temperature gradient, and temperature difference of the lower chord under the influence of different seasons.
Figure 10 illustrates the simulated transient temperature distribution of steel truss girder with two-layer decks on the summer and winter solstices. It can be found that there are apparent shadow areas at the web members and the lower deck. On the winter solstice day, the shadow generated by the web member at the lower deck was longer, because the solar altitude was lower than the summer solstice. Transient temperature distribution of steel truss girder at 16:00. (a) Summer solstice (b) Winter solstice.
Vertical temperature gradient of the girder
The vertical temperature gradient of the steel truss girder on the summer solstice is shown in Figure 11(a). Every moment, there was a significant temperature difference at the junction of the side girder and the web members, with a maximum temperature difference of 15.4°C at 9:30. This condition was attributed to the direct shining on the east side of the side girder, while the end of the web member was entirely in shadow. From 8:30 to 10:30, the upper deck covered the upper part of the web members, resulting in a significant temperature difference. The maximum temperature difference occurred at 8:30, at 8.1 m from the upper deck. Meanwhile, the temperature of the web member decreased from 41.8°C to 34.2°C, and the sheltering effect caused the maximum transient temperature difference of the same member to reach 7.6°C. Thus, the sheltering effect of the upper deck had a significant impact on the structural temperature field and cannot be ignored. Vertical temperature gradient of the steel truss girder: (a) 1 day on June 22, 2022; (b) largest temperature gradient.
The most adverse vertical temperature gradient of the steel truss girder in four seasons is plotted in Figure 11(b). The vertical temperature difference on the west side was more prominent than on the east, because of the higher air temperature when the sun shines directly on the west side girder. In summer, the maximum temperature difference on the west side was 17.3°C. In winter, the maximum temperature difference on the east side was 8.8°C.
Transverse gradient of the girder
The transverse temperature gradient of the upper deck was uniform under the sunlight, but the lower deck was different. The reason was that the upper deck and web members surrounded the lower deck. For convenient representation, the middle position of the lower deck was positioned at coordinate 0, with the farthest coordinates on the east and west sides being −7m and 7m, as given in Figure 12. Transverse temperature gradient of the steel truss girder.
From 7:30 to 8:30, there was a significant temperature difference of 2.8°C at the transverse coordinates of −5.5 m and 1.5 m. The lower deck within the transverse coordinates of 1.5 m to 7m and −7.5 m to 5.5 m was obstructed by the upper deck and web members, whereas the range of −5.5 m to 1.5 m was not obstructed, resulting in a significant temperature difference at the junction. At other times, the temperature distribution was even because the lower deck was entirely in shadow.
Transverse difference of the girder
The research reported in the existed reference had revealed that for similar steel truss cable-stayed bridges with two-layer decks, the transverse temperature difference of the lower chords significantly impacted the horizontal rotation angle (Wang and Ding, 2019). Therefore, it was necessary to study the variation of the transverse temperature difference of the lower chords. It should be said that the transverse temperature difference in this simulation refers to the difference between the lower chords on the east and west side.
As is demonstrated in Figure 13 (a), the variation pattern of the transverse temperature difference in different seasons is similar. Before sunrise, the temperature difference was almost 0. The temperature difference gradually increased as the sun rose and shined directly on the eastern lower chord. At noon, both sides of the lower chord were exposed to the same solar radiation, and the temperature difference fell to 0. The sun continued to move to the west side, and the temperature difference grew again. Transverse temperature difference of the lower chord: (a) different seasons; (b) different orientations.
The influence of bridge direction (azimuth) on the transverse temperature difference of the lower chord was studied based on the steel truss girder. As shown in Figure 13 (b), four typical directions were selected, namely −45°, 0°, 45°, and 90°, all in summer. When the azimuth was 0°, the temperature difference was the largest, 13.1°C. As the azimuth of the bridge increased to 90°, the maximum temperature difference declined to 1.76°C. This can be attributed to the inability of direct solar radiation to reach the lower chord on both sides. As a consequence, the direction of the bridge has a considerable effect on the maximum temperature difference.
Thermal deformation of steel truss girder during construction stage
The maximum cantilever stage is one of the most adverse construction stages before the bridge closure. After being affected by sunlight, the plane bending deformation of the steel truss beam is large. This section studied the vertical, transverse, and longitudinal thermal deformations of the upper chord (M2 and M3), side girder (M1 and M4), and lower chord (M5 and M6) during this stage through thermo-mechanical coupling. The specific location is given in Figure 14. Location of M1–M6: (a) longitudinal view of maximum single cantilever stage; (b) cross-section 1-1.
For cables, the cross-section is so small that the cable temperature is uniform. Zhu and Meng (2017) found that compared to the vertical temperature gradient inside the beam, the effect of temperature changes inside the steel cable on deflection was more significant. Therefore, this section considered that the temperature inside the cable varied sinusoidal, and there was a lag time between its extreme and the air temperature (Cao et al., 2011).
Vertical deformation of the girder
On the summer solstice and winter solstice, the real-time vertical deformation of members M1–M6 is shown in the Figure 15. As the sun rose, the temperature of the cable and the vertical temperature difference of the steel truss girder continued to increase. Under the combined effect, each member's downward vertical deformation varied greatly. At midday, the change in vertical deformation of the members was small, because of the decreasing vertical temperature difference (Figure 11). On the summer solstice, the maximum vertical deformation occurs at the upper chord M2, which is 7.2 cm. On the winter solstice, the maximum vertical deformation of the upper chord M2 was 5.1 cm. Vertical deformation of members: (a) summer solstice; (b) winter solstice.
Longitudinal deformation of the girder
The longitudinal deformation of the girder in the real-time temperature field is presented in Figure 16. In the morning, the longitudinal deformation of the east side members M1–M2 was larger; In the afternoon, the value of the west side members M3–M4 was larger. Moreover, the deformation of members M5–M6 was generally smaller than that of members M1–M4. This can be attributed to the sheltering effect of the upper structure on the lower structure, resulting in low temperature of M5–M6. In Figure 16(a), at 15:00 on the summer solstice, the deformation of the side girder M4 reached the maximum value of 6.4 cm among all members, while that of the lower chord M5 was 4.7 cm. The difference in deformation between members M4 and M5 was 1.7 cm. In Figure 16(b), during a day in winter, the maximum longitudinal deformation of the side girder M4 was 4.2 cm, and the maximum difference between the members was 1.3 cm. Longitudinal deformation of members: (a) summer solstice; (b) winter solstice.
Transverse deformation of the girder
The variation of the transverse deformation of the girder on the summer solstice and winter solstice is given in Figure 17. The overall deformation trend of members M1–M6 was the same, with little change at night and intense change at daytime. On the summer solstice, the maximum deformation during the day reached 5.9 cm; On the winter solstice, the maximum deformation was 3.9 cm. By comparing the transverse temperature difference of the lower chord calculated in Figure 13 (a), it was found that the variation between the transverse deformation and the transverse temperature difference was consistent, and they almost reached their maximum and minimum values simultaneously. Transverse deformation of members: (a) summer solstice; (b) winter solstice.
The Results of Correlation Analysis Between the Transverse Deformation of the Girder at Cross-section 1-1 and the Transverse Temperature Difference.
Conclusions
By considering geo-meteorological factors and the sheltering relationships, a temperature field analysis algorithm was developed, accurately simulating the temperature field and thermal deformation of the bridge during construction. The main conclusions of this study can be summarized as follows: (1) By comparing with experimental data, the accuracy of the method for analyzing the sunlight temperature field of the steel truss cable-stayed bridge with two-layer decks was verified. (2) Due to the obstruction of the upper structure, apparent shadow areas appeared at the web members and lower deck, causing a maximum temperature difference of 7.6°C in the bridge structure. It can be found that sheltering effect of the upper structure has a significant impact on the temperature field, which cannot be ignored. (3) Study on the transverse temperature difference of the lower deck shows that the maximum transverse temperature difference of this north-south oriented bridge can reach 13.1°C, which is much larger than that of the same type of east-west oriented bridge, which is 1.8°C; and the comparison of transverse temperature difference and transverse deformation of different seasons reveals that the transverse temperature difference of the lower chord is the main factor controlling the transverse deformation of the steel truss girder during the construction stage. (4) Sunlight has a large impact on the deformation of steel truss girder during construction. On the summer solstice, the maximum vertical, longitudinal, and transverse deformations of the main members are 7.2 cm, 6.4 cm, and 5.9 cm, respectively.
Footnotes
Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by the National Natural Science Foundation of China (Grant No. 52078333).
