The Multi-Scale Model Method for U-Ribs Temperature-Induced Stress Analysis in Long-Span Cable-Stayed Bridges through Monitoring Data

: Temperature is one of the important factors that affect the fatigue failure of the welds in orthotropic steel desks (OSD) between U-ribs and bridge decks. In this study, a new analysis method for temperature-induced stress in U-ribs is proposed based on multi-scale ﬁnite element (FE) models and monitoring data First, the long-term temperature data of a long-span cable-stayed bridge is processed. This research reveals that a vertical temperature gradient is observed rather than a transverse temperature gradient on the long-span steel box girder bridge with tuyere components. There is a linear relationship between temperature and temperature-induced displacement, taking into account the time delay effect (approximately one hour). Then, a multi-scale FE model is established using the substructure method to condense each segment of the steel girder into a super-element, and the overall bridge temperature-induced displacement and temperature-induced stress of the local U-rib on the OSD are analyzed. The agreement between the calculated temperature-induced stresses and measured values demonstrates the effectiveness of the multi-scale modeling strategy. This approach provides a valuable reference for the evaluation and management of bridge safety. Finally, based on the multi-scale FE model, the temperature-induced strain distribution of components on the OSD is studied. This research reveals that the deﬂection of the girder continually changes with the temperature variation, and the temperature-induced strain of the girder exhibits a variation range of approximately 100 µε .


Introduction
Orthotropic steel desks (OSD) are widely used in long-span steel box girder bridges due to their lightweight and excellent bearing capacity [1]. The evaluation of OSD's fatigue performance, considering the coupling effect between vehicle loads and environmental temperature effects, has become an important research topic [2,3]. The asphalt concrete pavement on the OSD, as a temperature-sensitive material, gradually reduces its elastic modulus and the overall stiffness of the bridge deck as the temperature rises. Consequently, the reduced stiffness of the road will inevitably weaken the top plate ability to constrain out-of-plane deformation on the OSD [4]. Researchers have observed that numerous bridges, especially box girder bridges, are subject to temperature-induced stress due to the nonlinear temperature gradient of the structure caused by solar radiation, atmospheric temperature and other external factors [5][6][7][8][9][10]. However, there are limited studies on the temperature effect evaluation of U-ribs on the OSD at present. Therefore, it is necessary to analyze temperature-induced stress and evaluate the performance of the OSD on long-span steel bridges. local models. If all local models need to be concerned, the multi-scale FE model can have similar DOFs to the global fine FE model, which reduces the potential advantage.
The aforementioned issue can be solved by implementing the substructure method. A group of elements is condensed into a super-element based on the linear response in the substructure method [38]. Li et al. [39] established a substructure model considering U-ribs for the mid-span segment on Runyang Bridge. The model can be used for the global dynamic analysis and hot spot stress analysis of local welds. Zhu et al. [40] built a multi-scale FE model using the substructure method to simulate the concrete box and steel box girders in a cable-stayed bridge in ANSYS. The substructures were connected by the main DOF coupling. Comparing the other two modeling methods, the substructure method has significant advantages when applied to structures with many repetitive geometries. For such geometries, the super-element is generated and copied to different locations. This advantage enables the substructure method to effectively solve the large DOFs problem of the long-span bridge. Connected by standard segments, the girder of long-span cablestayed bridges contains a large number of repetitive geometric structures. Therefore, in this paper, a multi-scale FE model based on the substructure method is used to analyze the temperature-induced stress of U-ribs on a long-span cable-stayed bridge.
When evaluating the structural performance of the U-ribs on the Sutong Yangtze River Bridge (STB) under the influence of temperature, the length of the girder on the STB is 2088 m while the size of a U-rib on the girder is 300~400 mm, resulting in a scale difference of approximately 5000 times. If the global fine FE method is used to analyze the temperature-induced strain of the U-ribs, the model would involve a large number of DOFs, leading to significant hardware resource consumption, long computation time and even infeasibility of calculation. Therefore, the multi-scale FE model method is used to address these challenges. It is the first time a multi-scale FE modeling method based on substructure technology is used to analyze the temperature-induced stress of U-ribs on long-span bridges through monitoring data.
The framework of U-ribs state evaluation through the multi-scale FE model and monitoring data are shown in Figure 1. It can be seen from Figure 1 that the global response temperature-induced displacement and the local response temperature-induced stress can be obtained by the SHM system and multi-scale FE model. The rest of the paper is organized as follows. The long-term monitoring data of STB are analyzed in Section 2. In Section 3, the substructure method is used to establish the multi-scale FE model and the thermal field and temperature-induced effects are analyzed and discussed. Section 4 provides the concluding remarks of the paper.

STB and Its SHM System
The STB is located downstream of the Yangtze River of China and connects the cities of Nantong and Suzhou. It is a double towers cable-stayed bridge with three piers arranged on the side spans ( Figure 2a). Initially, the STB held the record for the world's largest cable-stayed bridge with a span of 1088 m, and the total length of the girder is 2088 m. The STB's steel box girder is connected to different components, including the top plate, bottom plate, diaphragms, longitudinal plates (or longitudinal trusses), top plate Uribs, bottom plate U-ribs, tuyeres (Figure 2b). The STB towers and piers are constructed using reinforced concrete structures. The bridge is equipped with a total of 272 stay cables. Figure 2c illustrates the general layout of thermometers, displacement gauges, strain gauges and cable accelerometers of the STB SHM system. Structural thermometers are installed in the mid-span section of the girder (section . Displacement gauges are placed at the end of the girder (sections 1-1 and 9-9) to measure the longitudinal deformation of expansion joints. Strain gauges are installed in three sections of the girder: the mid-span section, the 1/4 section of the main span (section 4-4), and the girder section fitted at the north tower (section 3-3).
The schematic diagram in Table 1 depicts the sensor layout and SHM data for a day of the mid-span section (section . The red and blue lines represent the test results of the temperature sensors (T-5-#) and strain sensors (S-5-#), respectively.

Structural response
Structural field analysis

SHM system
Load information

STB and Its SHM System
The STB is located downstream of the Yangtze River of China and connects the cities of Nantong and Suzhou. It is a double towers cable-stayed bridge with three piers arranged on the side spans ( Figure 2a). Initially, the STB held the record for the world's largest cable-stayed bridge with a span of 1088 m, and the total length of the girder is 2088 m. The STB's steel box girder is connected to different components, including the top plate, bottom plate, diaphragms, longitudinal plates (or longitudinal trusses), top plate U-ribs, bottom plate U-ribs, tuyeres (Figure 2b). The STB towers and piers are constructed using reinforced concrete structures. The bridge is equipped with a total of 272 stay cables. Figure 2c illustrates the general layout of thermometers, displacement gauges, strain gauges and cable accelerometers of the STB SHM system. Structural thermometers are installed in the mid-span section of the girder (section . Displacement gauges are placed at the end of the girder (sections 1-1 and 9-9) to measure the longitudinal deformation of expansion joints. Strain gauges are installed in three sections of the girder: the mid-span section, the 1/4 section of the main span (section 4-4), and the girder section fitted at the north tower (section 3-3).
The schematic diagram in Table 1 depicts the sensor layout and SHM data for a day of the mid-span section (section . The red and blue lines represent the test results of the temperature sensors (T-5-#) and strain sensors (S-5-#), respectively.

Top Plate Top Plate U-Rib Bottom Plate Bottom Plate U-Rib
In the data graph, X label is hour, blue Y label is strain (unit: microstrain), and red Y label is temperature (unit: centigrade).

Top Plate Top Plate U-Rib Bottom Plate Bottom
In the data graph, X label is hour, blue Y label is strain (unit: microstrain), and perature (unit: centigrade).

Top Plate Top Plate U-Rib Bottom Plate Bottom Plate U-Rib
In the data graph, X label is hour, blue Y label is strain (unit: microstrain), and red Y label is tem perature (unit: centigrade).

Top Plate Top Plate U-Rib Bottom Plate Bottom Plate U-Rib
In the data graph, X label is hour, blue Y label is strain (unit: microstrain), and red Y label is temperature (unit: centigrade).

Top Plate Top Plate U-Rib Bottom Plate Bottom
In the data graph, X label is hour, blue Y label is strain (unit: microstrain), and perature (unit: centigrade).

Top Plate Top Plate U-Rib Bottom Plate Bottom Plate U-Rib
In the data graph, X label is hour, blue Y label is strain (unit: microstrain), and red Y label is tem perature (unit: centigrade).

Top Plate Top Plate U-Rib Bottom Plate Bottom Plate U-Rib
In the data graph, X label is hour, blue Y label is strain (unit: microstrain), and red Y label is tem perature (unit: centigrade).

Monitoring Data and Statistical Analysis of the Temperature and Structure Response
The steel box girder of the investigated bridge is longer in the longitudinal direction compared to the lateral and vertical directions. Additionally, the longitudinal curvature is small, and the cross-section along the longitudinal direction changes uniformly. Therefore, it is assumed that the thermal field along the longitudinal direction of the bridge is consistently exposed to solar radiation and atmospheric temperature. The lateral and vertical temperature distribution of the STB steel box girder is focused. Previous studies have found that the temperatures at the top and bottom plates of the bridge remain consistent along the longitudinal direction, while the temperatures at the web, which is the vertical component connecting the top and bottom plates, show variations between the upstream and downstream of the bridge [41]. However, the analysis of the STB temperature monitoring data does not reveal any temperature gradient along the lateral direction of the top plate. The same observation holds for the bottom plate, top plate U-ribs and bottom plate U-ribs. Moreover, the web temperatures at both sides are the same as shown in Figure 3a, which is different from previous research results. The temperatures of diaphragms are consistent with those of the webs. This is because the tuyere, as a non-structural component, is directly exposed to solar radiation, while the diaphragm and the web are located inside of the tuyere and are not directly affected by solar radiation. The girder without the tuyere component generally exhibits a temperature gradient in the lateral direction. However, after the tuyere is installed, the web is not directly exposed to the sun, resulting in the disappearance of the lateral temperature gradient. Figure 3b presents the vertical temperature distribution of the mid-span cross-section. The components are arranged sequentially from high to low temperatures, namely, the top plate, top plate U-rib, diaphragm, bottom plate U-rib and bottom plate, respectively. Figure 3c illustrates the cross-section temperature gradient data, which are compared with the BS 5400 code. According to the data collected by the temperature sensor, the measured temperatures of the top plate and the top plate U-ribs comply with BS 5400, while the measured temperatures of the diaphragm and U-ribs of the bottom plate exceed the values specified by BS 5400. This raises potential concerns for the structural design of the bridge and requires more attention. Based on the above statistical analysis, the thermal field of the STB can be analyzed using a 1D vertical thermal field approach.

Monitoring Data and Statistical Analysis of the Temperature and Structure Response
The steel box girder of the investigated bridge is longer in the longitudinal direction compared to the lateral and vertical directions. Additionally, the longitudinal curvature is small, and the cross-section along the longitudinal direction changes uniformly. Therefore, it is assumed that the thermal field along the longitudinal direction of the bridge is consistently exposed to solar radiation and atmospheric temperature. The lateral and vertical temperature distribution of the STB steel box girder is focused. Previous studies have found that the temperatures at the top and bottom plates of the bridge remain consistent along the longitudinal direction, while the temperatures at the web, which is the vertical component connecting the top and bottom plates, show variations between the upstream and downstream of the bridge [41]. However, the analysis of the STB temperature monitoring data does not reveal any temperature gradient along the lateral direction of the top plate. The same observation holds for the bottom plate, top plate U-ribs and bottom plate U-ribs. Moreover, the web temperatures at both sides are the same as shown in Figure 3a, which is different from previous research results. The temperatures of diaphragms are consistent with those of the webs. This is because the tuyere, as a non-structural component, is directly exposed to solar radiation, while the diaphragm and the web are located inside of the tuyere and are not directly affected by solar radiation. The girder without the tuyere component generally exhibits a temperature gradient in the lateral direction. However, after the tuyere is installed, the web is not directly exposed to the sun, resulting in the disappearance of the lateral temperature gradient. Figure 3b presents the vertical temperature distribution of the mid-span cross-section. The components are arranged sequentially from high to low temperatures, namely, the top plate, top plate U-rib, diaphragm, bottom plate U-rib and bottom plate, respectively. Figure 3c illustrates the cross-section temperature gradient data, which are compared with the BS 5400 code. According to the data collected by the temperature sensor, the measured temperatures of the top plate and the top plate U-ribs comply with BS 5400, while the measured temperatures of the diaphragm and U-ribs of the bottom plate exceed the values specified by BS 5400. This raises potential concerns for the structural design of the bridge and requires more attention. Based on the above statistical analysis, the thermal field of the STB can be analyzed using a 1D vertical thermal field approach. The temperature-induced displacement of the girder on the STB is typically related to the average temperature [42]. Figure 4a shows a significant correlation between the average displacement and average temperature. As the girder temperature increases, the displacement also increases, and vice versa. However, this trend is not consistently observed. There is a time delay effect in which the temperature data lag behind the displacement data for a certain period. Perhaps this is because the temperature sensors are installed inside the steel box girder, and it takes a certain time for the heat to transfer to the The temperature-induced displacement of the girder on the STB is typically related to the average temperature [42]. Figure 4a shows a significant correlation between the average displacement and average temperature. As the girder temperature increases, the displacement also increases, and vice versa. However, this trend is not consistently observed. There is a time delay effect in which the temperature data lag behind the displacement data for a certain period. Perhaps this is because the temperature sensors are installed inside the steel box girder, and it takes a certain time for the heat to transfer to the thermometers. The time delay effect between the structural displacement response and the temperature load on the STB is approximately one hour. Figure 4b presents the correlation between the average temperature difference ∆T and the average displacement difference δ u . Both values are calculated as the differences from their initial values. The correlation does not follow a linear relationship. The nonlinear relationship is mainly attributed to the time delay effect. By considering the time delay effect and shifting the temperature axis in the negative direction, a linear relationship between displacement and temperature can be observed (Figure 4c). thermometers. The time delay effect between the structural displacement response and the temperature load on the STB is approximately one hour. Figure 4b presents the correlation between the average temperature difference ∆T and the average displacement difference . Both values are calculated as the differences from their initial values. The correlation does not follow a linear relationship. The nonlinear relationship is mainly attributed to the time delay effect. By considering the time delay effect and shifting the temperature axis in the negative direction, a linear relationship between displacement and temperature can be observed (Figure 4c). It is shown in Table 1 that there is a significant correlation between temperature and temperature-induced strain for a single day. Owing to the temperature of the top plates being higher than the bottom plates, the top plates (S-5-2, 4 and 6) represent compression, while the bottom plates (S-5-7 and 9) exhibit tension. The compressive strain of top plates increases continuously with increasing temperature, and vice versa. For instance, if the temperature rises by 20 °C, the compressive strain of the top plate increases to 60 µε. On the other hand, the tensile strain of bottom plates increases continuously with the rising temperature, and vice versa. When the temperature increment of the bottom plate reaches 7 °C, the tensile strain of the bottom plate increases to 80 µε. However, the temperatureinduced strain of the U-ribs is influenced by the combined temperature of the decks and U-ribs. Between 10:00 and 15:00, the top plate U-ribs are under tension. In other time periods, the top plate U-ribs are under compression. This phenomenon occurs because the temperature of the top plates is higher than the top plate U-ribs between 10:00 and 15:00. Therefore, the top plates exhibit compression and the top plate U-ribs exhibit tension. However, the bottom plate U-ribs are under tension throughout. Figure 5 shows the correlation between temperature and longitudinal temperatureinduced strain of the girder under the influence of seasonal temperature variations in one year. The temperature and longitudinal strain were obtained from the temperature sensor T-5-5 and strain gauge S-5-5 on the top plate of the mid-span section on the STB. Firstly, the temperature and longitudinal strain of the first three days of each month were selected as representative data, and the data for each hour were averaged. Then, the monthly temperature and longitudinal strain data were averaged, respectively, resulting in a total of 12 monthly average temperature and 12 monthly average strain data points. The results indicate a sinusoidal relationship between seasonal temperature and longitudinal strain within a year. The top plate temperature reached its maximum value in August and its minimum value in February. On the other hand, the top plate strain reached its maximum value in September and its minimum value in February. It is shown in Table 1 that there is a significant correlation between temperature and temperature-induced strain for a single day. Owing to the temperature of the top plates being higher than the bottom plates, the top plates (S-5-2, 4 and 6) represent compression, while the bottom plates (S-5-7 and 9) exhibit tension. The compressive strain of top plates increases continuously with increasing temperature, and vice versa. For instance, if the temperature rises by 20 • C, the compressive strain of the top plate increases to 60 µε. On the other hand, the tensile strain of bottom plates increases continuously with the rising temperature, and vice versa. When the temperature increment of the bottom plate reaches 7 • C, the tensile strain of the bottom plate increases to 80 µε. However, the temperatureinduced strain of the U-ribs is influenced by the combined temperature of the decks and U-ribs. Between 10:00 and 15:00, the top plate U-ribs are under tension. In other time periods, the top plate U-ribs are under compression. This phenomenon occurs because the temperature of the top plates is higher than the top plate U-ribs between 10:00 and 15:00. Therefore, the top plates exhibit compression and the top plate U-ribs exhibit tension. However, the bottom plate U-ribs are under tension throughout. Figure 5 shows the correlation between temperature and longitudinal temperatureinduced strain of the girder under the influence of seasonal temperature variations in one year. The temperature and longitudinal strain were obtained from the temperature sensor T-5-5 and strain gauge S-5-5 on the top plate of the mid-span section on the STB. Firstly, the temperature and longitudinal strain of the first three days of each month were selected as representative data, and the data for each hour were averaged. Then, the monthly temperature and longitudinal strain data were averaged, respectively, resulting in a total of 12 monthly average temperature and 12 monthly average strain data points. The results indicate a sinusoidal relationship between seasonal temperature and longitudinal strain within a year. The top plate temperature reached its maximum value in August and its minimum value in February. On the other hand, the top plate strain reached its maximum value in September and its minimum value in February.

STB Multi-Scale Modelling Using the Substructure Method
By providing both global and local information to assess the structural safety performance, the multi-scale FE model can not only be used to perform static and dynamic analysis at the global level but also can analyze the stress of local components such as local strain, fatigue and the fracture of cracks.
The STB multi-scale FE model was established through ANSYS. Figure 6a shows the division of the girder into segments on the STB. The segment among the stay cables is selected as a segment model considering the connections between the girder, cables, towers and piers. The diaphragm enhances the load-bearing capacity of the segment at the end of the girder and the segment at the tower. Each segment is considered a substructure. Accordingly, the segment with a longitudinal length of 6.7 m is divided into a substructure GEN3, located at the end of the girder; the segment between stay cable No.2, with a longitudinal length of 48 m, is divided into a substructure GEN2, located at the tower; the standard segment between the adjacent stay cables, with a longitudinal length of 16 m, is divided into a substructure GEN0 and the standard segment between the contiguous stay cable with a longitudinal length of 12 m, is divided into a substructure GEN01. As the longitudinal truss of the mid-span segment is designed differently from the standard segment GEN0, it was divided into a single substructure GEN1. Overall, there are 53 standard segment substructures GEN0, 11 standard segment substructures GEN01, 1 mid-span segment substructure GEN1, 1 segment substructure GEN2 under the tower, and 1 segment substructure GEN3 at the end of the girder, all positioned on the side of the symmetrical axis on the STB. Figure 6d presents the FE model of the segment substructure GEN1 on the STB, modeled according to the design drawing. The segment, including top plates, bottom plates, webs, diaphragms and U-ribs, is simulated through the SHELL63 shell element in ANSYS. The top and bottom flanges of the longitudinal truss are welded with the top and bottom plates by steel plates, and thus the SHELL63 element is used to simulate the top and bottom flanges. The diagonal braces of the longitudinal truss adopt angle steel, and therefore, the BEAM4 element is used to simulate the longitudinal braces.
In the STB multi-scale FE model, each segment substructure model (Figure 6c) should be condensed into a super-element (Figure 6b), which is the condensation of a group of elements into one element using the substructure method. The purpose of this condensation is to simplify the model and improve computational efficiency. Once the super-element matrices have been formed, they are stored in a file and can be used as normal finite elements in subsequent analyses. The master nodes (purple nodes in Figure 6c) are selected to connect segment models and the other components in the multi-scale FE model.

STB Multi-Scale Modelling Using the Substructure Method
By providing both global and local information to assess the structural safety performance, the multi-scale FE model can not only be used to perform static and dynamic analysis at the global level but also can analyze the stress of local components such as local strain, fatigue and the fracture of cracks.
The STB multi-scale FE model was established through ANSYS. Figure 6a shows the division of the girder into segments on the STB. The segment among the stay cables is selected as a segment model considering the connections between the girder, cables, towers and piers. The diaphragm enhances the load-bearing capacity of the segment at the end of the girder and the segment at the tower. Each segment is considered a substructure. Accordingly, the segment with a longitudinal length of 6.7 m is divided into a substructure GEN3, located at the end of the girder; the segment between stay cable No.2, with a longitudinal length of 48 m, is divided into a substructure GEN2, located at the tower; the standard segment between the adjacent stay cables, with a longitudinal length of 16 m, is divided into a substructure GEN0 and the standard segment between the contiguous stay cable with a longitudinal length of 12 m, is divided into a substructure GEN01. As the longitudinal truss of the mid-span segment is designed differently from the standard segment GEN0, it was divided into a single substructure GEN1. Overall, there are 53 standard segment substructures GEN0, 11 standard segment substructures GEN01, 1 mid-span segment substructure GEN1, 1 segment substructure GEN2 under the tower, and 1 segment substructure GEN3 at the end of the girder, all positioned on the side of the symmetrical axis on the STB. Figure 6d presents the FE model of the segment substructure GEN1 on the STB, modeled according to the design drawing. The segment, including top plates, bottom plates, webs, diaphragms and U-ribs, is simulated through the SHELL63 shell element in ANSYS. The top and bottom flanges of the longitudinal truss are welded with the top and bottom plates by steel plates, and thus the SHELL63 element is used to simulate the top and bottom flanges. The diagonal braces of the longitudinal truss adopt angle steel, and therefore, the BEAM4 element is used to simulate the longitudinal braces.
In the STB multi-scale FE model, each segment substructure model (Figure 6c) should be condensed into a super-element (Figure 6b), which is the condensation of a group of elements into one element using the substructure method. The purpose of this condensation is to simplify the model and improve computational efficiency. Once the super-element matrices have been formed, they are stored in a file and can be used as normal finite elements in subsequent analyses. The master nodes (purple nodes in Figure 6c) are selected to connect segment models and the other components in the multi-scale FE model. In the substructure model, all nodes (with the exception of the master nodes) are eliminated through the substructure technique. More specifically, a group of elements is condensed into one super-element, with representative nodes of the segment model interface selected as the master nodes. In the case of segment substructure model GEN2, three additional rows of master nodes are additionally selected based on the connection of the beam between the substructure and tower, as well as the connection between the substructure and stay cable No.1. The DOFs of the master nodes on the interface of each segment substructure were then coupled. The above models can be directly applied to stress field analysis. However, it i essary to investigate the thermal field analysis before analyzing the temperature-ind displacement/stress. Since the longitudinal thermal variation of the bridge is not obv only the thermal field analysis of segment substructures with different geometry fea was carried out. The thermal field was calculated using the heat-transfer analysis o segment substructures GEN0, GEN01, GEN1, GEN2 and GEN3, which had not yet condensed into super-elements. The structural elements should be replaced by th elements. The 3D thermal shell element SHELL57 replaced the elastic shell ele SHELL63. Similarly, the convection link element LINK34 was used instead of the 3D tic beam element BEAM4. The bridge deck pavement has a great influence on the r of thermal field analysis. Therefore, the SOLID70 thermal solid element was employ  Owing to the length of the girder being long, the mass block was used inside the steel box girder at the piers to endure the warping deformation of the girder at these locations. The mass block was simulated by structural mass element MASS21 at the corresponding master nodes. In addition, the secondary dead load exerts a great impact on the static and dynamic performance of the bridge structures and is non-negligible. The secondary dead load of the STB deck was also simulated via MASS21. Based on the design drawing, the secondary dead loads of the bridge deck were converted into the mass of the MASS21 elements. The girder model was built following the completion of all the above steps.
The tower and pier of the STB were modeled by the 3D beam element BEAM188, which could accurately simulate the hollow cross-section characteristics of the concrete structure. Nonlinear truss element LINK10 was employed to simulate the stay cable, which was set to tension-only. The stay cable force of the bridge design was converted into the initial strain to simulate the initial tension of the stay cables. The stay cable experiences sag due to its self-weight, and cannot be calculated using the tensile truss element. In this case, the sag effect needs to be considered by the equivalent elasticity modulus. The Ernst formula (Equation (1)), initially proposed by German scholar Ernst in 1965, is widely used for this purpose. The formula assumes that the cable weight is uniformly distributed along chord length (not arc length) and neglects the cable's ability to bear bending moments. As a result, the shape of the stay cable follows a parabolic curve.
where E l is the corrected elastic modulus of the stay cable; E e is the initial elastic modulus of the stay cable; σ is the stress of the stay cable; γ is the unit weight of the stay cable; l x is the horizontal projection length of the stay cable. As STB bearings in the piers constrain the transverse and vertical displacement of the girder, these constraints were used at the joints between the piers and the girder. With the presence of the support beam between the tower and the girder, constraints were used to connect the support beam and the girder. The foundation at the bottom of the towers and three piers were constrained by fixed constraints. The multi-scale FE model was developed by coupling the master node DOFs on the girder super-elements with the nodes on stay cables, towers, and piers. Figure 5b presents the assembly of the STB multi-scale model by various components.
The proposed model was verified using SHM data prior to its application for the temperature-induce response analysis. The results obtained from the modal analysis of the multi-scale FE model were compared with measured values reported by Zhang and Chen [43]. The first three order modes of the multi-scale FE model were identified as the first-order longitudinal mode, the first-order symmetric lateral bending mode and the first-order symmetric vertical bending mode, with corresponding frequencies of 0.0693, 0.1069, 0.1884 Hz, respectively. By comparing with the frequency in the reference [43], the errors are determined as 8.1%, 5.3%, and 3.5%, respectively. The values calculated with the proposed method are consistent with the test values, indicating the proposed modeling strategy is reliable and the model can be adopted for further analysis.
The above models can be directly applied to stress field analysis. However, it is necessary to investigate the thermal field analysis before analyzing the temperature-induced displacement/stress. Since the longitudinal thermal variation of the bridge is not obvious, only the thermal field analysis of segment substructures with different geometry features was carried out. The thermal field was calculated using the heat-transfer analysis on the segment substructures GEN0, GEN01, GEN1, GEN2 and GEN3, which had not yet been condensed into super-elements. The structural elements should be replaced by thermal elements. The 3D thermal shell element SHELL57 replaced the elastic shell element SHELL63. Similarly, the convection link element LINK34 was used instead of the 3D elastic beam element BEAM4. The bridge deck pavement has a great influence on the results of thermal field analysis. Therefore, the SOLID70 thermal solid element was employed to simulate the bridge deck pavement, which is made up of 5.5 cm thick epoxy asphalt concrete on the bridge deck. During the thermal field analysis, the secondary dead load was deleted. When the stress field model was analyzed, the temperature results calculated by the thermal field were substituted into the segment substructures GEN0, GEN01, GEN1, GEN2 and GEN3 as the element internal force, and then each segment substructure was condensed into a super-element.
The advantages of the multi-scale FE model were illustrated by comparing with the global fine FE model. As shown in Table 2, the DOFs of the global fine FE model are  11,274,550, while the DOFs of the multi-scale FE thermal field model proposed in this paper are only 471,721. Based on the ANSYS run-time stats module, the number of the nodes, elements and DOFs of the global fine thermal field FE model is about 22 times, 25 times and 24 times that of the multi-scale FE model, respectively, and the time required for thermal analysis is about 25 times longer. As shown in Table 3, the DOFs of the global fine stress field FE model are 11,274,550, while the DOFs of the multi-scale FE model proposed in this paper are only 69,272. The number of nodes, elements and DOFs of the global fine stress field FE model is approximately 113 times, 2280 times and 162 times that of the multi-scale FE model, respectively, and the time required for static self-balance analysis (under the action of gravity load and measured cable force, the calculated value of girder deformation is consistent with the measured value) and temperature-induced stress analysis is about 1630 times and 40 times longer. Therefore, if the global fine FE method is used to analyze the temperature-induced strain of U-ribs, the model has a significantly large number of DOFs, leading to the consumption of extensive hardware resources and calculation time, and in some cases, it may not even be feasible to perform the calculation. However, the multi-scale FE model analysis can save a significant amount of computing time on the premise of accuracy, which is convenient for subsequent parameter analysis and structural performance evaluation.

Thermal Field Analysis
The thermal parameters of the steel components and deck pavement asphalt concrete, including thermal conductivity, specific heat capacity and density, were determined as shown in Table 4 [44][45][46]. Generally, the atmospheric temperature one hour before sunrise is considered the initial temperature of the bridge. The temperature distribution inside and outside the STB box girder at 6:00 a.m. is relatively uniform, so it is appropriate to take the atmospheric temperature at 6:00 a.m. as the initial atmospheric temperature. The time interval for the thermal field analysis is set to one hour. Thermal field analysis was performed to determine the boundary conditions, including absolute radiation temperature T v and overall heat transfer coefficient h o . The relation between the radiation and reflection balance of the deck surface on the bridge is described as Equation (2) [44]: where q B is the structural radiation; q k is the thermal irradiation; q J is the solar direct radiation; q H is the scattered radiation; q R is the reflected radiation; q G a is the atmospheric radiation and q UR is the surface effective radiation. The equations for the above parameters can be found in reference [45]. Equation (2) is a fourth-order transcendental equation of T v and cannot be solved directly. Therefore, the following formula is used to simplify the equation [45]: where ε bl is the structural radiation coefficient; C s is the Stefan-Boltzmann constant, C s = 5.677 × 10 −4 W/m 2 K 4 ; T a is the atmospheric temperature; h o = h r + h k , h r is the coefficient of radiation heat transfer and h v is the exchange heat transfer coefficient.
The absolute radiation temperature T v can be determined as the following formulas [45]: T v,bot = T a + γ uk a bk h r +h k where T v includes T v,top , T v,bot and T v,web . T v,top is the absolute radiation temperature of the top plate; T v,bot is the absolute radiation temperature of the bottom plate; T v,web is the absolute radiation temperature of the web; a bk is the shortwave radiation absorptivity; h is solar altitude angle; ε a is the atmospheric radiation coefficient, ε a = 0.82; γ uk is the ground reflection coefficient; I J is the direct solar radiation intensity and I H is the solar scattered energy; β is included angle between web and horizontal plane; α is the azimuth of the sun, which is the included angle between the projection of the line from the observer to the sun on the ground plane and the due south direction; α w is the azimuth of the external normal of the web. The steel and asphalt shortwave radiation absorptivity were considered to be 0.685 and 0.9, and the steel and asphalt surface radiation coefficients were assumed to be 0.8 and 0.92, respectively [44]. For more detailed information on these expressions, please refer to [45]. Based on the measured meteorological data and calculation parameters of the thermal field, the thermal boundary conditions were calculated by Equation (4). Then these boundary conditions were applied to the segment substructures for transient thermal analysis. Figure 7 presents the thermal field of the substructure GEN1 calculated for one day. Through conducting statistical analysis of the errors between calculated and measured values, it was found that the maximum error of temperature for the top plate is 8.2%, and the weighted average error within a day is 3.9%. The maximum error and the weighted average error for the bottom plate are 7.8% and 3.5%, respectively. For the web, corresponding values are 7.1% and 3.5%, respectively. For the diaphragm, they are 5.6% and 4.1%, respectively. For the top plate U-rib, they are 9.9% and 6.2%, respectively. For the bottom plate U-rib, they are 5.0% and 2.1%, respectively. There is an agreement between the predicted and measured temperatures. This verifies the effectiveness of the numerical model and the heat transfer analysis. and the weighted average error within a day is 3.9%. The maximum error and the weighted average error for the bottom plate are 7.8% and 3.5%, respectively. For the web, corresponding values are 7.1% and 3.5%, respectively. For the diaphragm, they are 5.6% and 4.1%, respectively. For the top plate U-rib, they are 9.9% and 6.2%, respectively. For the bottom plate U-rib, they are 5.0% and 2.1%, respectively. There is an agreement between the predicted and measured temperatures. This verifies the effectiveness of the numerical model and the heat transfer analysis.

Temperature-Induced Structural Responses
The temperature-induced responses can be investigated by combing the temperature time history with the multi-scale FE stress field model. In this section, the focus is on the analysis of longitudinal temperature-induced displacements at the end of the girder and the longitudinal temperature-induced strains on the sensor layout cross-section. First, the results of thermal field analysis are imported into each segment substructure model (Figure 6c) and then the temperature load is treated as the internal force, each segment sub-

Temperature-Induced Structural Responses
The temperature-induced responses can be investigated by combing the temperature time history with the multi-scale FE stress field model. In this section, the focus is on the analysis of longitudinal temperature-induced displacements at the end of the girder and the longitudinal temperature-induced strains on the sensor layout cross-section. First, the results of thermal field analysis are imported into each segment substructure model ( Figure 6c) and then the temperature load is treated as the internal force, each segment substructure model is condensed into a super-element ( Figure 6b); at last, these superelements are applied to the multi-scale FE model for stress field analysis. The isotropic linear constitutive model for thermal expansion is adopted. Temperature-induced displacement can be directly calculated from the master nodes of the girder end in the multi-scale FE model. Figure 8a,b depicts the temperature-induced displacements at the end of the girder. The calculated values of the FE model are in agreement with the measured values of the D-9-1 displacement gauges. However, there is a certain deviation between the temperatureinduced displacements of the D-1-1 displacement gauges and the calculated values. It should be noted that the STB multi-scale FE model is assumed symmetry with respect to the mid-span cross-section of the girder, resulting in consistent displacements calculated by the FE model on both sides. However, this assumption does not match the actual situation. Figure 8c,d displays the correlation between the temperature and temperature-induced displacements. The statistical results indicate the accuracy of the correlation between temperature and temperature-induced displacement based on the proposed method.
ues. It should be noted that the STB multi-scale FE model is assumed symmetry with respect to the mid-span cross-section of the girder, resulting in consistent displacements calculated by the FE model on both sides. However, this assumption does not match the actual situation. Figure 8c,d displays the correlation between the temperature and temperature-induced displacements. The statistical results indicate the accuracy of the correlation between temperature and temperature-induced displacement based on the proposed method. The temperature-induced displacement can be calculated directly. However, in order to obtain the temperature-induced stress and strain, an additional step is required. At first, it is necessary to reverse the displacement of the master nodes back to the segment substructure to calculate the displacement of the slave nodes; then the element temperatureinduced stress and strain in the segment substructure are calculated based on the displacement of the slave node. Figure 9 shows that the calculated values of decks and U-ribs are consistent with the measured values. Equation (5) is used as an evaluation indicator to The temperature-induced displacement can be calculated directly. However, in order to obtain the temperature-induced stress and strain, an additional step is required. At first, it is necessary to reverse the displacement of the master nodes back to the segment substructure to calculate the displacement of the slave nodes; then the element temperature-induced stress and strain in the segment substructure are calculated based on the displacement of the slave node. Figure 9 shows that the calculated values of decks and U-ribs are consistent with the measured values. Equation (5) is used as an evaluation indicator to compare the calculated and measured values of temperature-induced strain. The evaluation indicator value ranges between 0 and 1, with a value closer to 1 indicating a better similarity between the two sets of values.
where I is the evaluation indicator, S ai is the calculated value of temperature-induced strain for i-th sampling point, S ei is the measured value of temperature-induced strain for i-th sampling point, n is the number of samples for one day. The evaluation indicator values for S-5-1 to S-5-6 in the top plate or top plate U-rib are 0.96, 0.97, 0.94, 0.95, 0.75 and 0.98, respectively. The evaluation indicator values for S-5-8 to S-5-10 in the bottom plate or bottom plate U-rib are 0.90, 0.77 and 0.93, respectively. Therefore, compared with the traditional FE model, the multi-scale FE model based on the substructure method can be adopted to analyze the temperature-induced stress of the local members (i.e., U-ribs) on long-span bridges.
The evaluation indicator values for S-5-1 to S-5-6 in the top plate or top plate U-rib are 0.96, 0.97, 0.94, 0.95, 0.75 and 0.98, respectively. The evaluation indicator values for S-5-8 to S-5-10 in the bottom plate or bottom plate U-rib are 0.90, 0.77 and 0.93, respectively. Therefore, compared with the traditional FE model, the multi-scale FE model based on the substructure method can be adopted to analyze the temperature-induced stress of the local members (i.e., U-ribs) on long-span bridges. In order to evaluate the variations in the strain distribution across the longitudinal and vertical directions of the girder cross-section for one day, strains at different cross-sections and positions were analyzed at 6:00 and 14:00, respectively. Figure 10a shows that the top plates (SA-N-1, 3, 5) are subject to compression, while the bottom plates (SA-N-2, 4, 6) suffer from tension at 14:00. Moreover, the girder exhibits an upward bending behavior. However, the inverse phenomenon is observed in Figure 10b, whereby the top and bottom deck are subject to tension and compression, respectively, at 6:00, and the girder exhibits a downward bending behavior. Similar strain trends are observed at the top and bottom plate U-ribs in Figure 10c,d, respectively. This is attributed to the greater temperature of the top girder compared to the bottom girder at 14:00, while an inverse temperature distribution along the vertical direction is present at 6:00. The maximum variation range of the temperature-induced strain at the top plate and bottom plate is approximately 100 µε over the course of one day.
bottom deck are subject to tension and compression, respectively, at 6:00, and the girder exhibits a downward bending behavior. Similar strain trends are observed at the top and bottom plate U-ribs in Figure 10c,d, respectively. This is attributed to the greater temperature of the top girder compared to the bottom girder at 14:00, while an inverse temperature distribution along the vertical direction is present at 6:00. The maximum variation range of the temperature-induced strain at the top plate and bottom plate is approximately 100 με over the course of one day.

Conclusions
The performance evaluation of OSD by considering the temperature load is very important. In previous studies, the strain of local components needed to be analyzed in the

Conclusions
The performance evaluation of OSD by considering the temperature load is very important. In previous studies, the strain of local components needed to be analyzed in the global fine finite FE model. In this paper, based on substructure technology, the temperature-induced stress of U-ribs is evaluated using a multi-scale FE model on a longspan cable-stayed bridge. The following key conclusions are obtained.
(1) The temperature-induced stress of U-ribs on the STB was analyzed based on monitoring data and the multi-scale FE method. This method can be applied to other long-span bridges to address the issue of low computational efficiency in analyzing U-ribs in the global fine FE model. (2) Analysis of monitoring data indicates that the long-span steel box bridge with the tuyere components exhibits a vertical temperature gradient rather than a transverse temperature gradient. The correlation between temperature-induced displacement and temperature demonstrates a linear relationship once the time delay effect is considered. The temperature-induced strain of the top plates and bottom plates is influenced by the temperature between them. The temperature-induced strain of U-ribs is influenced by the temperature of the decks and U-ribs. Furthermore, the seasonal temperature and longitudinal strain over time within a year exhibit a sinusoidal relationship.
(3) A multi-scale FE model, which can effectively reduce the calculation time based on the substructure method, has been established to analyze the temperature-induced stress of U-ribs on long-span bridges. The accuracy of the multi-scale FE model results for the temperature-induced stress of U-ribs has been confirmed through monitoring data. (4) By evaluating the temperature-induced strain during the highest and lowest temperatures of one day on the multi-scale FE model, it indicates that the deflection of the girder, a key index for bridge design and SHM assessment, exhibits dynamic changes in response to temperature loads. The temperature-induced strain of the top and bottom plates displays a maximum variation range of approximately 100 µε.

Recommendation
Due to the limitations of research time, the evaluation and discussion are focused on the longitudinal temperature-induced strain of specific segments on the bridge, specifically the top plate, bottom plate, and U-ribs. In future research, the temperature-induced strain of web plates, diaphragms, and diaphragms can be studied in depth, as well as more segments of the girder. Additionally, it is recommended to investigate the influence of temperature load on the fatigue behavior of U-ribs under varying thermal conditions.