Tensioning-Phase Box Girder Deformation Prediction Model Based on Ant Colony Algorithm and Residual Correction

,


Introduction
It has become a worldwide problem that midspan lag and box girder cracks are common in prestressed concrete (PSC) continuous box girder bridges [1,2].The method of strengthening with a stay cable system (SCS) can fairly well solve this kind of problems, ensures the safety of the bridge structure, and extends the service life of the bridge; this method gradually gains the attention of the engineering community [3].However, affected by the natural environment and anthropic factors, the deformation mechanism of existing bridge is still unclear [4], which brings great difficulty to the deformation prediction of main girder strengthened by such system in tension phase.As a result, how to predict the deformation of the main girder becomes one of the key technical problems in the cable-stayed system reinforcement project.
At present, the deformation calculation methods of bridge structure mainly include theoretical analysis method, model test method, and numerical simulation method, intelligent prediction method.Among them, the theoretical analysis often adopts the idealized mechanical model, resulting in the relatively big difference between the theoretical calculation results and the measured values [5].Although the model test can directly reflect the deformation of the main girder, it has a high cost and some limitations in the deformation prediction of existing bridges.Numerical simulation method can well simulate the variation rules of the main girder deformation; however the selection of the calculation parameters and the constitutive model has a relatively big impact on the calculation results; therefore this method needs high techniques and does not have universal applicability [6][7][8].Therefore, based on the monitoring data, using intelligent forecasting model to predict the deformation of the bridge is a more practical approach; this method not only saves costs, but also improves the accuracy and speed of the calculation to a great extent.
Unequal interval grey model (1, 1) (UIGM (1, 1)) is a grey dynamic prediction model consisting of a single-variable, first-order differential equation; it is a popular model of grey system theory and has no special requirements and limitations on experimental data, so it has a very wide application field [9,10].The essence of establishing this model is making a single accumulated generating operation of raw data so that the generated data could have some regularity, and to obtain the fitting curve by establishing the differential equation model so that the system can be predicted [11].However, UIGM (1,1) can solve the change problem of single measuring point position only but not analyze the change rules of multiple positions associated with each other at the same time, so that the prediction model has the defects of low prediction efficiency, low prediction accuracy, and unstable overall performance.In practical application, numerical simulation method and UIGM (1,1) output predicted values lower than the measured values, failing to achieve high prediction accuracy.
Ant colony algorithm (ACA), a new heuristic bionic algorithm put forward by Dorigo M. et al., has the capacity of solving the problems of "complexity", "nonlinearity" and "uncertainty", "high cost" and obtains good results in longterm settlement prediction of tunnel [12].However, when solving the large-scale prediction problem, the model also has a low prediction accuracy and its residual sequence is usually alternatively positive and negative and has regularity which is not so regular.This phenomenon means the original data column often implies random components which make original sequence present the wave characteristic.In order to improve the prediction accuracy of the model, it is necessary to set up reasonable residual correction model to continue fixing it, so that the prediction accuracy could be improved.
According to the analyses above, this paper proposes a prediction method based on the combination of ant colony algorithm and residual combination correction model (ACAC model) to predict the deformation of main girder.First, ant colony algorithm is used to solve the problems of "complexity", "nonlinearity" and "uncertainty", "high cost", excavate the inherent rules of sampling data, and realize the quantitative prediction of the change in the next phase.Secondly, the periodic sequence function and sine function generated by harmonic variation are used to fit the residual sequence, which is used to modify ant colony algorithm model prediction value and compensate the deflect of low precision of the single prediction model, so that the model could achieve a higher prediction accuracy.The effect of this model was verified by the measured data of the main girder in the tensioning phase.The example analysis shows that the prediction accuracy of this model is better than that of UIGM (1, 1) model and ant colony algorithm model; this model has a certain theoretical and engineering application value.

The Basic Principle and Mathematical Model of Ant Colony
Algorithm.Ant colony algorithm [13] is inspired by ants feeding behaviors; scientists found that ants have no vision, and they cannot grasp the relevant geographic information in the region, but the ant colony can find the optimal path from the nest to the food source.Study finds that ants release a kind of pheromones on the path of travel when feeding, through the pheromones, the ant colony transmit messages with each other to find the optimal path.The ants are able to perceive the presence and intensity of such substances in their foraging process and guide their movement in the direction; they tend to move in the direction toward high intensity of such substances.In addition, pheromone is a kind of volatile substance that evaporates over time.Therefore, more of the pheromones in the path which is away from the nest has been evaporated, then the ants tend to walk on shorter path, and the collective foraging behavior of ant colony shows a positive feedback phenomenon.The shorter a path is, the more ants will pass it, then the greater the intensity of the pheromones left behind is, the greater the probability that the latter will choose the path is.The ant individuals work together through this kind of pheromones and finally choose the shortest path.
Combining the Traveling Salesman Problem (TSP), this paper introduces ant colony algorithm; its mathematical model can be described as follows: to set the number of cities in the plane to be ,  g () is the number of ants in city g at time t; that is, the total number of ant colonies is Suppose   ℎ () is the information function; it represents the residual pheromones concentration of ant k on the city line (g, h) at time t.Generally, it is assumed that the pheromones concentration  gh (0) on each path at the initial moment is equal.
Taboo list tabu k (k=1,2,. ..,m) is used to record the city where the ant k is passing; the collection tabu k is constantly updated with the search process of the ant.During the traversal, the ant calculates state transition probability based on the number of pheromones and the elicitation information in each route [14].  ℎ () represents the state transition probability of ant k transferring from city g to city h at time t, and its expression is as follows: In the formula,  is the relative importance of the trajectory,  ∈ [0, 5], and the larger  is, the more likely the ant is to choose the path that other ants have already passed; it indicates that the ants are more cooperative;  is the relative importance of visibility,  ∈ [0, 5]; it reflects the importance of the heuristic factor to the path selection of ant in the search process, the larger  is, the closer the state transition probability is to the greedy rule; allowed  = {1, 2, . . ., }tabu k represents the collection of candidate cities that ant k has not yet traveled.
The heuristic function  gh (t) represents the degree of expectation that the ant moves from city g to city h, and its expression is In the formula,  ℎ represents the distance between city g and city h.The smaller  gh is, the bigger  gh (t) is, the bigger   ℎ () is.
When the ants are traveling in the city, the pheromones are constantly updated according to formulae (4) and (5).
In the formula,  is the persistence coefficient,  ∈ (0, 1).The greater the value of  is, the greater the chance that the path that has been traversed will be chosen again is, and it will affect the globality and randomness of the algorithm.If  is too small, it will reduce the convergence rate of the algorithm.Δ  ℎ represents the increment of pheromones left by the ant k in the path (g, h) during this cycle.
Three models were developed by Colorni A. et Ant-Quantity model: , if the th ant moves from city  to city ℎ between  and  + 1 0, otherwise Ant-Density model: where Q is a constant, denoting the strength of pheromones; when its value is too large, it is easy to make the algorithm converge to the local minimum value, and if it is too small, it will slow the global convergence rate [16]; L k is the path length of the ant k.
The Ant-Cycle model uses the overall information; that is, the ant updates the pheromone on each path after a complete tour.The latter two models rely on local information; that is, the ant updates the pheromone in every step.The previous research has shown that the Ant-Cycle model outperforms the other two in dealing with the TSP [17,18], because the pheromone increment in this model only depends on the path length L k of ant k in this cycle and has nothing to do with the path length between g and h.In addition, the ants going through the shorter path will arrive at the destination earlier and release pheromone on the way back from the original path.With a relatively high pheromone concentration, the shorter path will attract even more ants.

Construction of Deformation Prediction Model Based on
Ant Colony Algorithm.According to the deformation curve of the main girder in three continuous tensioning phases, the deformation difference between the last two phases of a measuring point is closely related to (1) the deformation difference between the first two phases of this measuring point and (2) the result of the deformation difference between the first two phases of another measuring point minus that of yet another measuring point.This relationship is consistent with the positive feedback mechanism of the ACA; that is, the present behaviour can strengthen the future behaviour.Hence, the deformation difference between the first two phases of the main girder was analogized as the present behaviour and taken as the input parameters of the ACA prediction model.On this basis, it is possible to predict the deformation difference between the last two phases in an accurate manner.
Establishing reasonable heuristic function and information function is the key of constructing the deformation prediction model of the reinforced girder.Among them, the information function reflects the influence degree of the whole ant colony's cumulative experience on ants' path selection, and the inspiration function reflects the degree of influence of the ants' individual experience on selecting path.In order to apply the basic idea of ant colony algorithm into the deformation of the girder strengthened by a stay cable system, it is necessary to establish a similar colony path selection mechanism according to the actual deformation of the early tension phase of the main girder.In this paper, the monitoring points of the bridge are considered as the cities where the ants travel in TSP; in the course of selecting the path, the ant tends to move towards the measuring point of the larger deformation.As a result, the cumulative deformation of a site is regarded as an information function which can reflect the whole ant colony's cumulative experience.However, the difference of cumulative deformation of a certain measurement point in each period can be regarded as an inspiration function which reflects the experience of an individual ant.After calculating the transition probability based on the information function and the heuristic function of the path, ants move to the next point and release a certain concentration of pheromones on the path that they have passed in order to increase the cumulative experience of the whole ant colony.Thus, the positive feedback and autocatalytic mechanism of path selection of ant colony algorithm is accomplished.
Based on this, the initial information function  gh (0) and the heuristic function  gh (t) are established as follows: In the formula, Δ h and Δ g , respectively, represent the difference between the deformation of measure point h and measure point g in the two adjacent tension phases.ΔW min and ΔW sum , respectively, represent the minimum value and the sum of the measure point deformation difference between two adjacent tension phases in the search space.From formula (9), we can see that the initial information function is positively correlated with the deformation difference of h measuring point.In formula (10),  gh (t) reflects the desired value that the ant moves from point g to point h; that is, the larger the deformation difference of h point compared with that of g point, the easier the ant transfers h point.
Unlike the TSP, the ant movement of the ACA depends on the deformation difference of measuring point between tensioning phases rather than the path length.The pheromone update method is similar to that in Ant-Quantity model; i.e., the ants update the pheromone in each step.Thus, the concentration of pheromone Δ  ℎ (,  + 1) released by ant k on path (, ℎ) between time t and t+1 can be calculated by It can be seen from formula (11) that Δ  ℎ (,  + 1) depends on the cumulative deformation difference between tensioning phases   and  +1 when the ant arrives at measuring point ℎ.In other words, the value of Δ  ℎ is positively correlated with that of Δ ℎ .Hence, the greater the cumulative deformation difference when ant  moves from point  to point ℎ, the higher the pheromone concentration on path (, ℎ), and the easier it is for this path to attract ants.As the pheromone evaporates on each path over time, the pheromone concentration on path (, ℎ) at time  + 1 can be defined as Obviously, the value of  ℎ () is positively correlated with the likelihood for ant  to move from point  to point ℎ and the times that it arrives at point ℎ (i.e., the future deformation at point ℎ).
By using the parallel distributed search capability of ant colony algorithm, several ants are placed in the girder deformation area.These ants are ordered to search and update pheromones based on the aforementioned path selection mechanism and pheromones updating mechanism and also take self-organizing evolution of the algorithm.
At the end of the ant search, the number of times the ants arrive at each measuring point between tensioning phases  −2 and  −1 can be expressed as The total number of times the ants arrive at all measuring points is denoted as  sum Then, the predicted deformation difference of each measuring point between tensioning phases  −1 and   can be calculated as [19] where Δ sum  −1 is the total deformation difference of all measuring points between tensioning phases  −2 and Then, obtain the predicted cumulative deformation in tensioning phase   by adding up Δ Ŵ  (  ) and   −1 (  ): where   −1 (  ) is the measured deformation values of measuring point  in tensioning phase  −1 .

Parameters Selection Method of Ant Colony Algorithm.
The performance of ant colony algorithm prediction model is affected by parameters , , , , and ; the global search ability and convergence speed of ant colony algorithm should be considered to find the parameters' optimal values.At present, the selection methods of these parameters are not yet determined, and there is still no general method to determine the optimal parameter combination.The optimal combination of the parameters in this paper is determined based on the "three-step" method summarized by [12]; it is calculated out after many times' trial.The specific steps are as follows.
Step 1. Determine the ants' total number m.A larger number of ants can improve the global search ability and the stability of algorithm, but at the same time, it can make the pheromones concentration change at path evener and decrease the positive feedback effect of the algorithm.
According to the experience, the performance of the algorithm is more balanced when the ratio of city size N to the number of ants m is around 3/2, and the number of ants is determined based on it.
Step 2. Adjust the parameters with a larger value range in turn; they are , , and .Set  ∈ [1, 10000] and the step length is 50; set ,  ∈ [0, 5] and the step length is 0.5.
Step 3. Adjust the parameter  which has a smaller value range.Set  ∈ (0, 1), and the iterative step length is 0.1.
Ant colony algorithm prediction model can realize the prediction of the main girder deformation in the third equal interval tensioning phase with only two sets of deformation measured values at adjacent equal interval phases as input data.Because there might be large difference between the predicted deformation values of ant colony algorithm and the measured values, it is necessary to take a step further to use the residual correction model to modify the model and make up for the defect of low accuracy and poor adaptability of single prediction model.

Residual Combination Correction Model
Generally, the residual sequence of the main girder deformation is usually alternatively positive and negative and has periodicity which is not so regular.This phenomenon means the original data column has some random components and periodic components.In order to further improve the prediction accuracy of the model, generated by harmonic variation, the periodic sequence function and sine function are used to modify the model, and the former can reflect most of the periodic components of the residual sequence waveform, while the latter can reflect the remaining small fraction of random fluctuations in the waveform and be supplement to the former [20].The concrete implementation steps of the residual composite correction model are as follows.
Step 1. Build residual sequence in  k tension phase.Assume the longitudinal sequence of the bridge floor measurement point positions is So, residual sequence of  k tensioning phase is as follows: In the formula,    (  ) and Ŵ  (  ), respectively, represent the measured value and predicted value of the measurement point i in the F k tension phase.
Step 2. Fit    (  ) by sinusoidal function Ê1   (  ) generated by the average characteristics of the residual sequence waveform.
In the formula,  = (1/)∑  =1 |   (  )|; T is the average of the residual sequence    (  ) wave time span; its value follows the principle of the minimum quadratic sum 1 min of the fitting values residual for sine function; that is, Step 3. Fit    (  ) by periodic sequence function Ê2   (  ) generated by harmonic change. where Step 4. Introduce the weight coefficient (0 ≤  ≤ 1) and pair ( 18) and ( 19) to combine the weights, and the residual composite correction value Ê  (  ) is obtained.
In the formula, the value of  and  follows the principle of the minimum quadratic sum 2 min of the fitting values residual by the periodic function generated by harmonic transform and sine function; that is, Step 5. Based on the nearby principle of data, the most recent data can reflect the nature of main girder deformation.F k+1 phase prediction value is revised by the fitting value of the residual sequence at F k phase [21].Like that, add up the ant colony algorithm model prediction value Ŵ +1 (  ) and the correction value Ê  (  ); the prediction value of the final ant colony algorithm correction model is obtained as follows:

The Specific Implementation Steps of the Ant Colony Algorithm Correction Model
Step 1. Initialize the parameters , , , and .
Step 2. Determine the initial position of ants.First of all, according to the number of measuring points N, determine the scale of ant colony (assume the total number of ants is m).Then, calculate the deformation difference between each measuring point in the two adjacent tensioning phases.
According to the order of values, assign m ants to the measuring points from large to small and determine the initial locations of ants.
Step 3. Calculate the initial information function  gh (0) (formula ( 9)), the heuristic function  gh (t) (formula (10)), and the state transition probability   ℎ () (formula (2)), determine the measuring points that the ants will attain in next station, update the pheromones continuously according to rules of formula (12), and record the number of times that the ants arrive at each measuring point by the taboo table Tabu.
Step 4. After the termination of search, the predicted deformation values of measuring points in next phase can be calculated by formulae (14) and (15).If the measuring points in the predicted area are not many, it is advisable to set all m ants to be polymerized to a certain measuring point as an end condition.If the measuring points are many, it is best to set most ants to be polymerized to a measuring point as an end search condition [19].
Step 5.According to a certain step length, parameter T is set, and the sine function (formula (18)) is used to correct the residual sequence    , and  best corresponding to the minimum value of the cumulative residual error is output.
Step 6. Select T= best and set parameters ,  according to a certain step length.The residual combination correction model is used to modify the residual sequence, and the correction sequence Ê  (formula (20)),  best and  best corresponding to the minimum value of the cumulative residual error are output.
Step 7. Add up the ant colony algorithm model prediction value Ŵ +1 and correction value Ê  ; the final correction model prediction value Ŵ1  +1 can be obtained by formula (21).
The flowchart of deformation forecasting process of SCSenhanced main girder in tensioning phase by ACA correction model is shown in Figure 1.In order to grasp the deformation of the main bridge during the tensioning phases, 46 vertical displacement measuring points were set on the upstream and downstream of the whole bridge.This paper selects the upstream monitoring data to study, as shown in Table 1.
There is a certain deviation between the predicted and measured deformation values in the 1,785kN tensioning phase.This is attributable to the following two reasons.
First, the ants are not sensitive to the measuring points with small deformation in the early phase; thus, some measuring points are not visited by the ants, leading to a large residual value.Second, it can be seen from the calculation results of the residual sequence of 1,785kN tensioning phase that this residual sequence is usually alternatively positive and negative and has periodicity which is not so regular.This phenomenon means the original data column has some random components and periodic components.
In order to further improve the prediction accuracy of the ant colony algorithm prediction model, this paper takes fitting calculation on residual sequence of 1,575kN tensioning phase by residual combination correction model and takes it as the correction sequence of the predicted value of the 1,785kN tensioning phase.
Firstly, based on the ant colony algorithm prediction model, the deformation prediction value Ŵ1,575 of 1,575kN tensioning phase is obtained as follows: Ŵ1,575 =(-  Secondly, the residual sequence is modified by sine function, and the parameter  best is calculated.Set  ∈ [1, 10000] and the step length is 1; take iterative computation.When T=220s, 1 min =35.40mm 2 , the calculation result is best.The first 300 iteration results are shown in Figure 3.
Then, the periodic sequence function and sine function generated by harmonic variation are used to modify the residual sequence in combination.Set  ∈ [1, 10000] and the step length is 1; take iterative computation.When T=220s, =15 and =0.65,2 min =2.61mm 2 , the calculation result is best; the first 300 iteration results are shown in Figure 4.

Analysis Results of UIGM (1,1) and Its Two Residual
Correction Models.This paper regards the cable force data in tensioning phase as the original time sequences, regards the vertical deformation data obtained from the bridge floor survey points as the original date sequences, and builds the prediction model of main girder deformation strengthened by a stay cable system based on UIGM (1, 1) [22].This model takes the deformation data of the first four tension phases of every point as input values and predicts the deformation of the 1,785kN tensioning phase.
In addition, the modified model provided in [23,24] was adopted to correct the predicted value of UIGM (1, 1).The one-time residual correction UIGM (1, 1) (OC -UIGM (1, 1)) is introduced in [23].During prediction, this model first establishes the UIGM (1, 1) model of the residual value and then compensates the original predicted value with the sum of the predicted value of the residual UIGM (1, 1) and the original predicted value.
Under various factors, the main girder deformation and its trend cannot be determined completely, which is consistent with the highly random nature of Markov residual prediction.Through Markov residual prediction, it is possible to improve the random features of the sequence and optimize the predicted results.Therefore, the UIGM (1, 1) predicted value is modified by the Markov chain in [24].
(1) Predicted Results of OC -UIGM (1,1).Based on the predicted results of UIGM (1, 1), the deformation residuals of the first four tensioning phases were calculated.Taking the tensioning force   as the time sequence and the residual sequence    (  ) is used as the original sequence, UIGM (1, 1) was reestablished to obtain the predicted residual value Ê +1 (  ) of tensioning phase  +1 .Finally, the predicted residual values Ê +1 (  ) and Ŵ +1 (  ) were added to obtain the final predicted value Ŵ2  +1 (  ).
Here, measuring point C3 is cited as an example to describe the prediction process.According to the predicted results of UIGM (1, 1), the residual sequence of the first four tensioning phases was identified as    (  )=(0.00,-0.40, 0.46, 3.11).Then, the UIGM (1, 1) was reestablished with   =(630, 1,155, 1,365, 1,575) as the time sequence of the grey model and the residual sequence (0.00, -0.40, 0.46, 3.11) as the original data sequence.However, the modelling cannot be implemented by the common method, because the sign of the residual sequence is either positive or negative.In light of this, the translation transformation method in [25] was adopted as follows.
Step 1.The residual sequence was transformed into a nonnegative one by adding a translation quantity H.
Step 2. The new residual was predicted again.
Step 3. The predicted value was subtracted by the translation quantity H.
Thus, the modified predicted deformation value is as follows: Ŵ2 1,785 = Ŵ1,785 + Ê1,785 = 48.85+ 1.79 = 50.64 (2) Predicted Results of Markov Chain Residual Correction UIGM (1, 1) (MC -UIGM (1,1)).Markov chain prediction mainly forecasts the possible state of system variables at a certain time in future according to the current state and the change trend.According to a rational self-made standard, the residual sequence was divided into  states, namely,  1 ,  2 , ⋅ ⋅ ⋅ ,   .Thus, the transition probability of  steps can be expressed as being the number of data in the residual sequence transferred from state   to state   through the  steps and ∑  =1  ()  being the total number of original data in the residual sequence in state   .Then, the state transition matrix can be defined as where ∑  =1  ()  = 1.Generally, the state intervals are classified based on the mean and standard deviations of the data sequence, and the mean   of the two endpoint values of the interval is the centre of the i-th state interval.Here, the residual sequence    (  ) is divided into n states S, and r transfer probability vectors are used.Transition probability vectors (the rows of a matrix) are needed to consider the transition processes in r steps.The probability of residual state at the target moment is the sum of the r vectors, and the vector element is denoted as is the state weight and   is the centre of the state interval.Therefore, the predicted value of UIGM (1, 1), modified by Markov, can be expressed as Ŵ3  +1 (  ) = Ŵ +1 (  ) + M +1 (  ).
The residual sequence was divided into three state intervals []-1.9,]-0.6], []-0.6,]+0.6], and []+0.6,]+1.9] according to the mean ] and standard deviation  of the residual sequence of this point.The three state intervals correspond to states  1 ,  2 , and  3 and cover the entire residual range.For measuring point C3, the 3 state intervals were [-0.72, -0.03], [-0.03, 1.62], and [1.62, 3.40], and the corresponding states are shown in Table 2.The residual value of the 1,785kN tensioning phase was predicted by the state transitions of 630kN∼1,575kN tensioning phases, with a total of 3 steps of transitions.The state transition matrix can be expressed as Then, the transition probabilities are listed in Table 3.The state centres were calculated as  1 = −0.38, 2 = 0.79, and  3 = 2.51; the state weights were calculated as  1 =  1 / ∑  =1   = 0/1 = 0,  2 =  2 / ∑  =1   = 1/1 = 1, and  3 =  3 / ∑  =1   = 0/1 = 0.Then, the Markov residual of the 1,785kN tensioning phase can be derived as Thus, the final predicted value of the deformation in the said phase is For other measuring points, the deformation in the 1,785kN tensioning phase can be obtained similarly with the OC model and MC model.

Analysis Results of FEM.
The previous studies have attributed the deformation of PSC box girder bridges to the following factors: concrete shrinkage, creep, prestress loss, and the mismatch between actual and theoretical concrete amounts [26,27].Excessive deformation of such bridges may exacerbate the cracking and reduce the shear stiffness of box girder web, which, in turn, aggravates the deformation of the main girder.
The deformation features of the SCS-enhanced main girder in the tensioning phase were discussed using the finiteelement software.A finite-element model was established  There are many diagonal cracks, whose width falls within the allowable range.
The crack length is between 1/3 and 1/2 of the sectional size. 4 There are numerous diagonal cracks, whose width exceeds the allowable range.
The crack length is greater than 1/2 of the section size, and the mean crack interval is less than 30cm.

5
There are numerous through diagonal cracks, whose width far exceeds the allowable range.
The crack width is greater than 1.0mm and the mean crack interval is less than 20 cm.Elastic connection was adopted between the stay cable and the main girder.The finite-element model is presented in Figure 5.
The deformation increment caused by the diagonal cracks of the web was simulated by the shear stiffness degradation model [28] and the crack grading method [29] through the following steps.
(1) The main bridge was divided into several units.Then, the grade   of the diagonal cracks in each unit was determined according to the diagonal cracks of the web measured in 2013 [30] and the diagonal crack grading standard (Table 4).
(2) The limit  u on shear stiffness degradation factor (SSDF) was calculated by the following equation: where  V is the stirrup rate of the web;  is the diagonal crack angle.
(3)The SSDFs  were determined for all crack grades.When there is no diagonal crack, the elastic shear stiffness of the web can be determined as   =    V =    V /2(1 + ), with   ,   ,  V , and  ( = 0.20) being the shear modulus, elastic modulus, shear area, and Poisson's ratio of concrete.The SSDF  describes the degradation of shear stiffness in general cracking state ( =   ).The value of  can be calculated by (4) The modified Poisson's ratios were determined for all crack scales: (5) The shear stiffness degradation was simulated by modifying Poisson's ratio of each unit in the finite-element model, laying the basis for deformation analysis of the main girder.
The calculated SSDF and modified Poisson's ratio of our example are listed in Table 5.
According to the FEM results, the deformation of the measuring points in the 1,785kN tensioning phase amounts to ant colony algorithm correction model, the calculation results show that C = 0.060, p =1.00, and it shows that the prediction accuracy for main girder deformation of this paper's model is in excellent range and superior to that of the other two prediction models.
To sum up, the proposed model can achieve fast computation and high prediction accuracy using a few measured data through easy programming and simple operation.The case study reveals that the model is an excellent tool for predicting SCS-enhanced main girder deformation in the tensioning phase, shedding new light on main girder deformation control and prediction.
In the same way, the deformation values of the next tensioning phase (1,995kN) are extrapolated.Table 8 contains the prediction values obtained by FEM, UIGM (1,1), OC-UIGM (1,1), MC-UIGM (1,1), ACA model, and ACAC model.The results show that the forecast values of first five models are generally low.The C3 deformation values obtained by ACAC model are largest when the cable tension is drawn to 95% design tension value, and its value is 55.27mm.

Conclusion
It is difficult to accurately predict the girder deformation of PSC continuous box girder bridges enhanced by the SCS, owing to the complex deformation mechanism.The existing prediction models, lacking the modelling based on measured data, face low accuracy and poor adaptability in application.To overcome these defects, this paper combines the ant colony algorithm and the combined residual correction model into a new prediction model, which not only tackles the complexity, nonlinearity, uncertainty, and high cost, but also mines out the inherent deformation rules of main girder under the tensioning force.The new model was compared with finite element model, UIGM (1, 1), one-time residual correction UIGM (1,1), Markov chain residual correction UIGM (1,1), and ant colony algorithm model.The main conclusions are as follows.
(1) Compared with the unequal interval GM (1,1), onetime residual correction UIGM (1,1), and Markov chain residual correction UIGM (1,1), ant colony algorithm prediction model takes into account the time-space effect of deformation development of multiple points, and it can better dig out the inherent rules of the original data and make the prediction results more accurate, and its mean absolute error is 5.90%.
(2) Using residual combination correction model to modify the predicted result of the ant colony algorithm can well simulate the trend components and random fluctuation compositions of the main girder deformation, and it has high prediction accuracy and the mean absolute error of perdition values is 3.39%, the posterior error ratio is 0.060, and the accuracy level reaches level 1.The prediction precision is improved greatly; the paper's model has a certain theoretical and engineering application value.
(3) The model needs less measured data, the principle of the model is simple, and it is easy for programming, but the theory of optimal parameter selection needs to be further improved in the follow-up work so as to improve its forecasting efficiency and expand the range of adaptation.Unlike the finite element model, the model must be modelled based on three continuous and of equal interval measured data.

Figure 1 :
Figure 1: Deformation prediction flowchart of SCS-enhanced main girder in tensioning phase by ACA correction model.

Figure 2 :
Figure 2: Strengthening of Dongming Huanghe River Highway Bridge by application of a stay cable system.

5. 1 .
Project Profile.The main bridge of Dongming Huanghe River Highway Bridge belongs to PSC continuous rigid frame-continuous beam composite structure system; it has nine holes and one link, and the cross-diameter combination is (75 + 7 × 120 + 75) m with a total length of 990 m.The main girder and the bridge pier are poured with C50 and C25 concrete, respectively.To curb the excessive deformation and web cracking, the main girder was enhanced in 2003 by increasing the thickness of the web and adding external tendons.Nevertheless, new cracks and deformation were observed in some of the spans of the main bridge after the enhancement.In order to effectively improve the stress condition of the main bridge and raise the height of the midspan section, it was proposed to reinforce this bridge by the cable-stayed system in 2014, and this bridge was the first PSC continuous box girder bridge reinforced by this project in China[3].The specific connection construction of the bridge tower, the main beam, and the cable were as follows: the top part of the bridge tower was set as a fixed end of the cable by installing the steel anchor boxes and the sliding cable saddles.The trimmer beams, which were set laterally beneath the main girder, were attached to the main girder by based plates on the bottom of the box girder.The stay cable force was transmitted to the main beam through the cooperation of trimmer beams and based plates, as shown in Figure2.From the middle towers (61# and 62#) to the side towers (58# and 65#), from the long cables to the short cables, cables were tensioned symmetrically by 8 lifting jacks by level during the tension construction.After tensioning construction, the actual tensioning cable values of short and long cables are 90% of the design cable values of the bridge, namely, 1,890kN and 2,430kN.On the basis of the tensioning values of short cables (630kN, 1,155kN, 1,365kN, 1,575kN, and 1,785kN), predicted results are analyzed and studied in this article.

Figure 5 :
Figure 5: Finite element model of main girder strengthening by a SCS.

Table 1 :
Vertical deformation of the measuring points under different tensioning phases.

Table 2 :
State of residual sequence.

Table 6 :
Comparison of predicted and measured values.

Table 7 :
Comparison of evaluation results under different prediction models.

Table 8 :
Comparison of predicted values in 1,995kN tensioning phase.