Abstract

This study establishes a numerical model for beam-type steel–concrete composite specimens considering the corrosion of anchor bolts, with which studies the long-term deformation performance of the components under sustained load. The experimental results were compared with the calculated values obtained by combining the effective modulus method according to the CEB-FIP 1990 code and the ACI 209R code. Then a comparison with the calculation method of shrinkage and creep in standards (Standard Creep Method (SCM)) is made, and it shows that by supplementing the degradation of interface stiffness due to corrosion (using the defined modulus method (DMM)), the simulation results match better with the experimental results, confirming that this approach is suitable for analyzing the long-term load conditions of anchor bolt corrosion. Furthermore, based on the defined modulus method, the long-term behavior of composite beams under different loads and varying corrosion rates is studied.

1. Introduction

Steel–concrete composite bridges built in China during the early stages, which have been in service for nearly 30 years, have entered the first phase of their life cycle. They are starting to experience material degradation, corrosion, and structural damage due to accidents, all of which affect the durability of steel–concrete composite beam bridges. The issue of long-term load-carrying capacity degradation caused by bolt corrosion has gradually become apparent. Composite beams can be classified into external concrete composite beams and T-shaped composite beams according to their cross-sectional forms (Figure 1), with T-shaped composite beams being commonly used.

Wang et al. [1] conducted a long-term experiment on a steel–concrete composite beam and used finite element simulation analysis to simulate the long-term behavior of concrete using three methods: concrete material creep function, effective modulus method, and adjusted shear stiffness at steel–concrete interface. They proposed a calculation method for long-term deflection of composite beams. Ji et al. [2] derived a deflection calculation formula for steel–concrete composite beams under long-term loading by considering the interlayer slip effect between steel beams and concrete deck slabs, full-section shear deformation of steel–concrete composite beams, and shrinkage creep effect of concrete deck slabs. The derived formula was validated against experimental and finite element results, showing that the shrinkage creep effect of concrete has a significant influence on the deflection of steel–concrete composite beams. Han [3] analyzed the effects of concrete creep and shrinkage in composite beams on stress distribution and structural behavior. Through a combination of analytical calculation, numerical simulation, and experimental analysis, they derived analytical formulas for stress redistribution caused by internal constraints in statically determinate and hyperstatic composite beam sections, as well as direct calculation formulas for the time-varying stress of composite beams. They obtained a series of accurate and approximate solutions. Fan et al. [4, 5] studied the effects of concrete shrinkage creep and cracking on the long-term mechanical performance of composite beams. Their experiments showed that both long-term loading under negative bending moment and concrete shrinkage effect lead to concrete cracking and inhibit the creep deformation of concrete. Di [6] studied the influence of concrete shrinkage creep on the long-term performance of steel–concrete composite bridges by establishing a finite element analysis model. They analyzed and researched three parameters that affect the concrete shrinkage creep effect, providing reference value for the design and construction of steel–concrete composite bridges regarding the impact of concrete shrinkage creep effect. Zhang [7] used the finite element method to study the models of steel–concrete composite beams with different bolt corrosion rates. They proposed calculation formulas for the flexural bearing capacity and stiffness of corroded bolt steel–concrete composite beams. Xiong et al. [8] conducted experimental tests on eight standard connection members to study the static properties of corroded bolt connections after long-term loading and unloading, as well as the effects of corrosion and loading. They established calculation formulas for the shear bearing capacity of corroded bolt connections and the load–slip curve. Peng [9] studied the long-term performance of composite beams under the combined effects of load and bolt corrosion through long-term experiments. They proposed a time-varying formula for the long-term deflection of steel–concrete composite beams. Here are the latest research findings from scholars in 2023. Xiao [10] and others conducted experiments on single-sided specimens to study the static performance of welded nail connections under different degrees of corrosion. The results indicate that the mechanical performance of welded nail connections decreases with corrosion, leading to degradation in composite structural performance and reduced durability. Based on the experimental results, they fitted coefficient curves and proposed a calculation formula for the shear bearing capacity of corroded welded nail connections. Wang et al. [11] conducted numerical simulations by establishing ABAQUS finite element models to study the degradation law of shear performance of bolted connections under different rust distribution patterns. The results indicate that under the same rust depth conditions, as the rust height increases, the shear carrying capacity of the bolt decreases; under the same rust height conditions, as the rust depth increases, the shear carrying capacity and stiffness exhibit a linear decreasing trend; and under the same rust height and depth conditions, the shear carrying capacity and stiffness decrease most significantly due to uniform circular rust. Wang [12] studied the effects of severe corrosion on the mechanical properties of steel–concrete composite girder bridges. Through experiments, numerical analysis, and theoretical analysis, he primarily investigated the degradation law of mechanical properties after corroding the steel bars with industrial hydrochloric acid and the failure modes and shear bearing capacity of PBL shear connectors, as well as the influence of various factors on the shear bearing capacity of shear connectors, and proposed corresponding formulas and models. Wang et al. [13] investigated the influencing factors of the bending capacity of steel–concrete composite beams after corrosion and fatigue and established a calculation method based on fracture mechanics and fatigue residual strength theory. The results indicate that the corrosion rate of bolts, concrete strength, fatigue loading limit, and number of bolts are the main factors affecting the bending capacity of composite beams. You et al. [14] studied the long-term performance of continuous steel–concrete composite beams considering the effects of concrete shrinkage and creep. Using ABAQUS software, he established a finite element model and verified the reliability of the simulation results. Additionally, he conducted a parameter analysis on the influence of recycled coarse aggregates on the long-term performance of these beams. The study found that using the age-adjusted effective modulus method and considering the temperature field to account for concrete shrinkage effects can accurately simulate the long-term performance of composite beams. The incorporation of recycled coarse aggregates increases deflection to a certain extent, and different floor types should adopt different shrinkage and creep models. Zhang et al. [15] conducted a 4-year field test to investigate the temperature, shrinkage, and creep effects on steel–concrete composite girder bridges. Temperature field models (TFM) and mechanical behavior models (MBM) were verified and showed minor differences throughout the service stage. The study highlights the influence of meteorological factors on long-term behaviors of such bridges. Nan et al. [16] investigated the widespread application and economic benefits of steel–concrete continuous composite beams in buildings and bridges. Through a combination of elastic theory and numerical simulations, the study delved into the influence of slip and creep effects under long-term loads on structural performance, emphasizing the importance of considering the interaction between slip and creep. Zhao et al. [17] proposed a beam finite element model considering slip, shear lag, and time-dependent effects in steel–concrete composite box beams. The model employs a step-by-step method for time variable solutions, more accurate than the single-step algebraic method. Validation against ANSYS and existing test results on composite beam long-term performance attests to its effectiveness in instantaneous analysis. The model predicts the time-related behavior of simply supported composite beams, highlighting significant influences of concrete shrinkage and creep on structural responses. Ting et al. [18] reviewed the performance of shear connectors in steel–concrete composites exposed to hostile environments. Despite extensive use, inefficient connectors can weaken composite structures. Various connector types’ load transfer mechanisms, influences of hostile conditions like temperature and corrosion, and failure modes were evaluated, urging further research on connector deterioration. Du et al. [19] explores the shear behavior of corroded stud connectors in steel and bridge engineering. Experimental investigation reveals increased tensile deformation in corroded studs under dynamic loading, with higher strain rates enhancing shear fracture. Corrosion accelerates sensitivity to dynamic loads, yet ultimate shear stress remains unaffected.

In summary, most scholars have not considered the corrosion of connectors in the investigation of the long-term performance of composite beams. However, in practical engineering applications, corrosion of connectors is inevitable. While Zhang [7] did consider corrosion, the study only focused on the short-term flexural carrying capacity performance. Xiong et al. [8] only discussed the static performance of bolted connections after long-term load unloading due to corrosion, without addressing performance during long-term load application. Furthermore, no scholar has yet proposed an effective finite element analysis method for evaluating the long-term performance of composite beams under bolt corrosion [20]. Therefore, it is necessary to study the corrosion of connectors in composite beams. This paper investigates the long-term performance of composite beams by considering different bolt corrosion scenarios using numerical methods combined with relevant experimental data [9, 21].

2. Experiment Introduction

The long-term performance of steel–concrete composite beams was investigated by researcher Peng [9], considering different corrosion rates, presence of load, and loading age conditions. A total of 11 simply supported composite beams were designed for the experiment, divided into four categories of specimens: specimens corroded under load (SCB7∼SCB11), specimens not corroded under load (SCB2, SCB3), specimens corroded without load (SCB4, SCB5, and SCB6), and specimens not corroded without load (SCB1). Select SCB-1, SCB-2, SCB-3, SCB-4, SCB-5, SCB-7, SCB-8, and SCB-10 for numerical simulation. These specimens include both corroded and noncorroded components, with loading conditions of unloading, 7-day age loading, and 28-day age loading, in various combinations, to make the numerical simulation results more reliable. The basic parameters of the composite beams are shown in Table 1. The material parameters of the concrete test blocks from the same batch can be found in Table 2,. All experiments were conducted using C40 ordinary concrete with standard curing for 28 days, with a compressive strength not less than 40 MPa. The experiment focused primarily on the study of bolt corrosion. Using electrochemical principles, the specimens achieved the expected corrosion rate, while the maximum compressive stress σmax (at the top surface of the concrete flange at midspan) induced by long-term loading was controlled at around 0.4 fc. The long-term loading value was a concentrated load of P = 30 kN at midspan. The combined action of electrical corrosion and loading lasted for 225 days. The cross-sectional diagram and front view of the specimen can be seen in Figures 2 and 3.

3. Numerical Model Establishment and Theoretical Background

A numerical modeling analysis is conducted based on the corrosion of anchor bolts in commonly used composite structures. The feasibility of finite element analysis is verified by comparing the results with experimental data. On this basis, more working conditions are extended to study and analyze the long-term behavior of steel–concrete composite beams with corroded anchor bolts.

3.1. Finite Element Model

In this study, a simply supported constraint is adopted for the composite beam model, and DIANA software is used for numerical simulation. At one end of the flange of the I-beam at the support location, all translational degrees of freedom are constrained (Tl = T2 = T3 = 0), while at the other end, the lateral and vertical translational degrees of freedom are constrained (T1 = T3 = 0). The self-weight of the structure is applied through a gravity field, and the long-term load is applied as a surface load on the top surface of the concrete slab (Figure 4). The finite element model is a 1:1 scale model. In the figure, some parts of the concrete slab are hidden to facilitate a better observation of the steel and bolts. Eight-node hexahedral elements (HX24L) are used to simulate the steel and concrete in the specimen model. For the concrete and steel plates, hexahedral elements are used to ensure good adaptability and convergence of the elements. Steel reinforcement beam elements are used for the anchor bolts, which are embedded in the concrete to withstand shear forces. Through simulations with various grid sizes, taking into account the analysis accuracy and computational efficiency of the model, it is finally concluded that a grid size of 20 mm achieves better accuracy and higher computational speed. The convergence criteria include force and displacement convergence, indicating convergence of the iterative calculation results as long as either force or displacement converges during the iterative process. The convergence tolerance is set to 0.01, with a maximum of 50 iterations allowed. The mesh size of the steel and concrete is set to 20 mm. The total number of solid elements is 11,742, and the total number of steel reinforcement elements is 61.

3.2. Interface Element and Parameter Selection

The simulation of the interface behavior in composite beams includes two parts: the interface simulation between concrete and steel sections and the interface simulation between bolts and concrete. The interface simulation between concrete and steel sections is conducted using surface-to-surface contact elements with Coulomb friction properties. In the simulated experiments, a layer of epoxy resin is applied to the surface of the steel, resulting in minimal friction. Therefore, it can be assumed that the steel–concrete interface only provides vertical support. Hence, a large value is chosen for the normal stiffness modulus of the steel–concrete interface, while small values are selected for the two shear directions to simulate the behavior. The normal stiffness modulus of the interface between the steel and concrete is set to 2,000 N/mm3, while the tangential stiffness modulus is set to 0.001 N/mm3. The friction angle is set to 0.38 rad [22], and the cohesion is set to 0.001 N/mm2. The interface element opening model selected is the gap model.

Two approaches can be used to model the interface of bolts: (1) The first one is using embedded beam elements for the bolts, as shown in Figure 5. The bolts are set as beam elements and penetrate the interface elements between the steel section and concrete, being embedded in the steel flange [23]. The bolts are modeled as cylindrical geometries, and a cohesive-sliding model is applied for the interface. (2) The second one is using solid elements for the bolts, where the root of the bolt is fixed to the steel beam and other parts of the bolt are embedded in the concrete slab [24]. After modeling and analysis using both approaches, it was found that the results were consistent in terms of the long-term behavior of the composite beams (e.g., deflection–time curve). The only difference lies in the local stress distribution of the bolts. However, using discrete reinforcement elements to simulate the bolts proved to be more convenient, efficient, and significantly faster in terms of computational speed compared to using solid elements. Therefore, the later numerical analyses in this study adopted the embedded beam elements to simulate the bolts.

Embedded beam element settings: In the geometric properties setting of the anchor bolt, select the reinforcement type as circular beam bond-slip, with a diameter of 13 mm. In the data properties setting, select Bondslip Reinforcement/Foundation Pile as BEAM. Choose the interface continuity in the reinforcement properties interface, and in the analysis options under the model assessment, select to evaluate the reinforcement of the interface element. The anchor bolts do not need to be merged with the concrete nodes.

3.3. Simulation Method for Long-Term Behavior of Steel–Concrete Composite Beams

This study focuses on the long-term behavior of composite beams under normal service conditions, but localized regions in the concrete may reach the plastic stage. Therefore, the mechanical behavior of concrete is simulated using the total strain crack model [25]. The tensile curve is Hordijk’s tensile softening curve, the compressive curve is Maekawa’s cracking concrete curve, the tensile fracture energy is 0.1465 N/mm, and the Poisson’s ratio is 0.15. Other material properties can be found in Table 2 based on experimental data. The constitutive relationship is shown in Figure 6. The constitutive model selected for reinforcing bars, structural steel, and bolts is the von Mises plastic yield model, which incorporates plastic hardening. The stress–strain curve is represented by a trilinear form, as shown in Equations (1) and (2), with material properties detailed in Table 1. The stress–strain curve is illustrated in Figure 7.

In the above equation, and represent the yield strength and yield strain, respectively, while denotes the end strain of the yield plateau, with a value of 0.025. and represent the ultimate strength and ultimate strain, respectively, with set to 0.15. represents the elastic modulus of the reinforcing bars, while represents the linear slope of the reinforcing bars in the strain-hardening stage, calculated using Equation (2).

The simulation method for the long-term behavior of composite beams is as follows:(1)Calculation Method of Shrinkage and Creep in Standards (referred to as Standard Creep Method, SCM). The standard method considers two types of shrinkage in concrete [5, 26]. One type is inherent autogenous shrinkage, which is caused by the volume of cement hydration products being smaller than the volume of cement and water involved in the hydration reaction. The autogenous shrinkage curve is shown in Figure 8. The other type considered is drying shrinkage, which occurs when the distance between gel particles in the hydrated cement paste is less than the thickness of 10 water molecules, which leads to a tensile force generated by adsorbed water molecules to balance the molecular attraction between the particles, resulting in volume expansion [27]. Hence, the loss of this adsorbed water in drying will instead cause volume shrinkage. The shrinkage curve is shown in Figure 9 [28].

In the standard method, the concrete creep model adopts the Kelvin chain model [29]. A Kelvin element consists of an elastic component and a viscous component connected in parallel. Its physical significance can be described as a series of interconnected elastic components (springs) and viscous components (dampers). This is shown in Figure 10. En represents the elastic modulus of the nth Kelvin element, and ηn represents the viscous damping coefficient of the nth Kelvin element. The Kelvin chain is formed by connecting multiple Kelvin elements in series. The Kelvin chain is based on the aging theory of concrete and is often used to analyze time-dependent problems in reinforced concrete structures. The elastic modulus of concrete in the Kelvin chain is one of its important parameters, which is closely related to the age of the element, temperature, and humidity.(2)The defined modulus method (DMM) is used to account for the variation in normal stiffness modulus at the interface between corroded anchor bolts and concrete. After long-term corrosion, gaps exist at the interface between the anchor bolt and concrete, resulting in a decrease in interface stiffness. This leads to additional deflection in the composite beam. When linked with the degree of corrosion, the normal stiffness modulus then can be referred to for the study of the time-dependent behavior of corroded anchor bolts in composite beams.

4. The Numerical Comparative Study on the Long-Term Performance of Composite Beams Considering the Influence of Corroded Force on Interface Stiffness

Based on the experimental data from the reference [2], a study is conducted using different methods to simulate the long-term behavior of composite beams.

In the analysis using the SCM, many researchers adopted a multi-Kelvin element model to fit the experimental data. Due to its ability to capture both elastic and creep deformation which has a stable limit, the three-parameter model [25] is suitable for components with initial elastic deformation that undergo long-term sustained loading. Therefore, this study selects the three-parameter Kelvin model, which consists of an elastic component and a Kelvin element connected in series. The three parameters are denoted as E0, E1, τ1, and τ represents the time damping, where τ = η/E.

In the research utilizing the DMM to supplement the simulation of the corrosion effect on anchor bolts, the normal stiffness modulus is defined as the ratio of normal stress to normal relative displacement. After anchor bolt corrosion, the diameter of the anchor bolt decreases, and therefore, the relative displacement is considered with respect to the reduction in diameter (Δd). As steel reinforcement corrosion can cause rust expansion stress [20], the curve of rust expansion stress with respect to Δd is shown in Figure 11. The peak rust expansion force is calculated using the following formula: P = ft (0.94 c/d + 0.91) [20], where c is the thickness of the protective layer and d is the diameter of the anchor bolt. The normal stress after corrosion is determined by the initial normal stress and the rust expansion stress, thus determining the corresponding normal stiffness modulus in the corroded state. The deflection cloud diagram of the composite beam under long-term load simulated using the DMM is shown in Figure 12.

The comparative results are shown in Figure 13. SCB-1, SCB-4, and SCB-5 represent unloaded specimens, while the rest are loaded specimens. SCB-1, SCB-2, and SCB-3 are specimens without corrosion. When only concrete shrinkage and creep occur but without corrosion, the simulation results of both methods align closely and are represented by the same line. Table 3 presents the comparison between the long-term deflection (in millimeters) calculated by the two methods and the experimental results for a 225-day loading period. Currently, in China, when simplifying calculations for long-term deformations of composite structures, it is common practice to reduce the long-term elastic modulus of concrete. According to the “Code for Design of Steel Structures,” the long-term elastic modulus of concrete is specified to be half of the initial elastic modulus, while in the “Code for Design of Highway Bridges and Culverts Steel Structures,” it is taken as 0.4 times the initial elastic modulus. This paper will adopt the recommended reduction values from the “Code for Design of Steel Structures” to calculate the development of concrete strain and compare it with experimental measurements, simulation values, and the calculated values obtained from the formula proposed by researcher Peng [9], as shown in Table 3 for comparison.

From Figure 13, it can be observed that the DMM yields good agreement with the experimental results for all loading conditions, as it takes into account the influence of anchor bolt corrosion on stiffness and long-term deflection. On the other hand, when simulating using the SCM, the results deviate significantly from the experimental data under the corroded condition and only align well in the absence of corrosion, though not as well as the DMM method.

The reason for this difference is that the deflection of the composite beam is influenced by both concrete shrinkage and creep, and the effect of anchor bolt corrosion on long-term deflection needs to be considered. In the unloaded state, although the anchor bolts experience corrosion, the creep caused by their self-weight is relatively small compared to the influence of concrete shrinkage. Therefore, the error of the SCM method is relatively small. However, when considering loaded conditions, the failure to account for the effects of corrosion leads to larger deviations.

Currently, the calculation method for the long-term deflection of steel–concrete composite beams is based on the principle of effective modulus method. This method involves reducing the elastic modulus of concrete after long-term loading, thereby increasing the ratio of elastic modulus between steel and concrete. The calculation is performed using the formula for elastic deflection [1]. However, this method can only be applied to noncorroded components, and it cannot be used to calculate deflection for corroded components. The calculation of effective modulus values is related to the creep coefficient. The creep models in CEB-FIP 1990 and ACI 209R perform well. Therefore, in this study, the deflection of noncorroded components (SCB-1, SCB-2, and SCB-3) is calculated using the effective modulus method combined with the creep models from these two standards. The results are as shown in Table 3.The calculation results indicate that there is a significant discrepancy between the calculated values from CEB-FIP 1990 and the experimental values, suggesting an overly conservative approach. On the other hand, the calculated values from ACI 209R are close to the experimental values. Therefore, for noncorroded components, the ACI 209R code can be adopted for calculation.

In both DMM and SCM methods, the results from the DMM are closest to the measured values, with very small relative errors. Therefore, using DMM can provide better simulation. The SCM is more suitable for conditions where the anchor bolts are not corroded under long-term loading. The proposed DMM in this study is applicable to the study of concrete shrinkage and creep under both loaded and unloaded conditions. Moreover, it takes into account the corrosion of anchor bolts and its impact on the long-term behavior of composite beams. Therefore, the subsequent numerical studies in this paper will be conducted using the DMM.

5. The Impact of Load Strength, Bolt Corrosion Rate, and Bolt Diameter on the Long-Term Behavior of Composite Beams

5.1. Influence of Anchor Bolt Corrosion on the Long-Term Behavior of Composite Beams under Different Loading Strengths

Considering the limited number of specimens and the limited variation in load levels and anchor bolt corrosion on the long-term behavior of composite beams under sustained loading conditions, to further improve the analysis, a numerically validated method is used to supplement the experimental study and investigate the long-term deformation capacity of steel–concrete composite beams under more levels of anchor bolt corrosion and loading strengths. Through numerical analysis, it was found that when the corrosion rate exceeds 15%, the deflection of the composite beam is almost identical to that of a composite beam without bolts, indicating bolt failure. This finding is consistent with the conclusion reached by Zhang [7] in the study of the flexural load-carrying capacity of corroded bolted composite beams, where it was suggested that bolt failure occurs when the corrosion rate exceeds 13%. Therefore, the parametric analysis in this section will focus on studying the deflection for corrosion rates ranging from 5% to 15%, which is also the range studied by most scholars.

Figure 14 depicts the variation of midspan deflection with time under different loading strengths. The deflection of the composite beam increases with increasing duration under different levels of anchor bolt corrosion and tends to stabilize at various loading strengths. The larger the load, the more significant the variation of deflection with the rust rate. This means that the range between 5% rust rate and 15% rust rate, as shown in the graph, is wider. It indicates that under the coupled effect of load and rust rate, the deflection of composite beams will increase significantly. As shown in Table 4, at 225 days, under a load level of 20 kN, the increase is 105.42%, while under a load level of 50 kN, the increase reaches 131.41%. However, as the loading strength increases, the increase in deflection with time becomes more significant within the same range of corrosion rate limits. The relationship between deflection and corrosion rate for the same load but different corrosion rates is shown in Table 4.

Figure 15 illustrates the variation of deflection of a loaded specimen at a 7-day age (corrosion rate: 7.5%) under different loading durations. From the results, it can be observed that the deflections almost follow a linear relationship with load changes for the same duration. Moreover, as the loading duration increases, the increases in the deflections become smaller and tend to reach a stable value. It can be observed that from 100 days onward, the magnitude of deflection change decreases significantly. This is because shrinkage and creep mainly occur within the first 100 days, while the variation in deflection after 100 days is primarily due to the slow growth of creep and the effects of bolt corrosion.

Table 5 presents the deflection (mm) and its rate of change under different loads at a corrosion rate of 7.5%. It reflects the relationship between anchor bolt corrosion rate and deflection, which shows a significant linear growth. This indicates that the corrosion rate has a notable impact on the long-term deflection of composite beams. As the number of loading days increases, the rate of change also increases. This is because the corrosion rate varies with the number of loading days. At the beginning of loading, the corrosion rate is 0%. However, as the number of loading days reaches 225 days, the corrosion rate gradually increases to 7.5%. The coupling effect of corrosion and load becomes more severe, leading to an increase in the rate of change.

Figure 16 illustrates the variation of deflection with loads for different corrosion rates. As shown in the figure, the deflection increases approximately linearly with increasing corrosion rate under long-term sustained loading. Additionally, at high load levels, the change in deflection with corrosion rate is more significant compared to the low load conditions, indicating that corrosion rate has a larger impact on deflection under high load conditions but a smaller impact under low load conditions.

5.2. Influence of Different Bolt Diameters on the Long-Term Deflection of Composite Beams under Varying Corrosion Rates

As bolt diameter serves as a prominent parameter in the research conducted by numerous scholars on composite beams, this section will investigate the influence of bolt diameter on the long-term deflection of composite beams under bolt corrosion conditions. Currently, bolt specifications include φ8, φ10, φ13, φ16, φ19, φ22, φ25, and φ28, totaling eight types. The research results for these eight bolt specifications are as follows.

From Figure 17, it can be observed that under different corrosion rates, the trend of long-term deflection with bolt diameter remains relatively consistent. However, as the bolt diameter increases, the long-term deflection actually increases. This is because under normal operating conditions, for the same corrosion rate, a larger bolt diameter results in a greater reduction in effective cross-sectional area, thereby leading to an increase in long-term deflection. This does not imply that larger bolt diameters are inherently undesirable. In practical engineering, under the same corrosion environment, it takes more time for a larger diameter bolt to reach the same corrosion rate as a smaller diameter bolt. It is important to determine the bolt diameter based on actual engineering conditions and to implement effective corrosion prevention measures.

6. Conclusions and Future Perspectives

In this study, the defined creep method (DCM) and DMM were used to simulate eight tested specimens from previous research, analyzing the long-term behavior of steel–concrete composite beams with anchor bolts having different corrosion rates and loading ages. The results indicate that the DMM can effectively simulate the long-term deflection of composite beams. The DCM is suitable for simulating the long-term deflection of composite beams under loading conditions, but it fails to model the deflection under self-weight conditions, showing a deviation of over 60% from the experimental values. By using the DMM, the long-term deflection of composite beams under various corrosion rates and loading levels was studied. For every 10 kN increase in load, the deflection increased by 10%–35% compared to the previous load level. When the increment of corrosion rate was the same, the higher the load, the greater the influence of corrosion rate on deflection.

The corrosion rate of the bolts simulated in this study ranges from 5% to 15%, targeting simply supported steel–concrete composite beams under low corrosion rates. Through numerical analysis, it was found that when the corrosion rate exceeds 15%, the bolts begin to enter a failure state. This study did not delve further into this aspect. It is hoped that researchers will conduct experimental studies on the failure of bolts under high corrosion rates. Additionally, further analysis and validation are needed for the long-term performance of continuous beams with multiple spans. The long-term deflection calculation method for composite beams considering anchor bolt corrosion can be proposed based on the DMM. This study only focused on the long-term deflection under normal service conditions, and further research is needed to investigate the ultimate limit state of composite beams with anchor bolts subjected to long-term corrosion.

Data Availability

The data that support the findings of this study are available in [9].

Conflicts of Interest

The authors declare that they have no conflicts of interest.