Abstract

In the design of prestressing concrete structures, the friction characteristics between strands and channels have an important influence on the distribution of prestressing force, which can be considered comprehensively by curvature and swing friction coefficients. However, the proposed friction coefficient varies widely and may lead to an inaccurate prestress estimation. In this study, four full-scale field specimens were established to measure the friction loss of prestressing tendons with electromagnetic sensors and anchor cable dynamometers to evaluate the friction coefficient. The least square method and Bayesian quantile regression method were adopted to calculate the friction coefficient, and the results were compared with that in the specifications. Field test results showed that Bayesian quantile regression method was more effective and significant in the estimation of the friction coefficient.

1. Introduction

Prestressed concrete is used to reduce or offset the tensile stress of concrete under the action of loads, so as to improve the crack resistance and stiffness of components. Prestressed distribution plays an important role in the design and construction of prestressed concrete structures. For post-tensioned concrete structures, the friction characteristics between the steel strand and the pipe (also known as sheaths) are directly related to the elongation required by the tensioning operation and the prestress force distribution. The value of friction coefficient is usually obtained from the design specifications or the instructions provided by manufacturers. The recommended friction coefficient [14] ranges widely, which may result in the confusion and inaccuracy. The inaccurate calculation of friction loss can cause unreasonable camber and prestress reduction or even deflection of beams. Aiming to obtain the true force of the structure, many researches have been carried out to investigate by the calculation method [5] and the full-scale field tests [68]. When the friction loss data of prestress is obtained by tests, the corresponding friction coefficient can be got according to the prediction equation of prestress. In addition, the post-tensioning method cannot guarantee the accuracy of stress measurements of steel strands [9, 10]. The total prestress loss at the beam end and the stress distribution are measured by the pressure sensor at the anchorage point [11] and the electromagnetic (EM) sensor [1214], respectively.

However, it is possible for the designers and builders to determine the coefficients erroneously. Although many studies have been conducted to determine the friction coefficient by various methods, no standard calculation method has been established [1517]. The least square method [18] is the most common calculation method for estimating the friction coefficient. First and foremost, such method needs to assume that the random error term obeys the zero mean and homoscedastic distribution. Secondly, it is essentially based on the conditional mean model, mainly in the statistical analysis framework of the mean values. Nevertheless, too much pursuit of the average gives rise to the loss of much other information. Furthermore, the field test environment is less controllable than the precision laboratory, and the data are more discrete. Considering the cost factor, the sample size of the field test is often small, so the accuracy of the results estimated by the least square method is low.

In this study, a Bayesian quantile estimation method [19, 20] was proposed in order to make up for the deficiency of the least square method in the calculation of the friction coefficient. The Bayesian quantile regression method does not require mandatory distribution of test data, and the friction coefficient can still be effectively estimated even when the sample size is small. Based on the monitoring of the prestressed tension process of four full-scale test beams in the construction site, the above two calculation methods of friction coefficient were compared.

2. Friction Coefficients in Post-Tensioning System

During the tensioning process of the steel strand, the contact between the steel strand and the pipe can cause two types of friction loss: the curvature friction that resulted from the curved part of the pipe and the wobble friction due to the deformation of the pipe.

The regression analysis model [1, 21] of prestressing friction loss parameters is presented aswhere fpj is the stress in the prestressing steel at jacking; x is length of a prestressing tendon from the jacking end to any point under consideration; k is wobble friction coefficient; µ is coefficient of friction; θ is sum of the absolute values of angular change of prestressing steel path from jacking end, or from the nearest jacking end if tensioning is done equally at both ends, to the point under investigation (rad.); and e is the base of Napierian logarithms.

Taking logarithms on both sides of an equal sign equation (1), the equation can be derived as

Let

The equation can be presented as

Assuming that x and θ are deterministic variables, Y is random variables, k and µ are constants, and ε is the error term, then equation (4) can be regarded as a binary linear regression model.

Table 1 lists the friction coefficients recommended for strands in galvanized metal pipes in various provisions, which is the most widely used type of pipes. It can be found that various provisions give different ranges, some of which are ranges while others are fixed values. In general, unpractical calculation of friction loss can lead to unreasonable buckling and reduction of prestressing. Overestimating the value of friction loss may lead to excessive prestressing, which will bring about excessive buckling and shrinkage of members. This study attempted to estimate the coefficient of friction based on Bayesian quantiles and verify the current practice, shown in Table 1.

3. Methodologies to Derive Friction Coefficients

Using EM sensors to obtain the prestress distribution of the tendon, the friction coefficient can be derived as follows.

To start with, the differences between the anchor point and other measurement points are substituted into equation (4), which are overconditioned simultaneous equations because k and µ are unknown; a number of equations can be formed. Therefore, the unique coefficient usually cannot be obtained in this way, and a statistical approach is required instead. In an attempt to solve this problem, two equations were chosen among many equations as shown in equation (5), and then k and µ were calculated. Choosing another set of two equations may cause k and µ to be slightly different. Then, the average coefficients of different simultaneous equations are used as the friction coefficients:

The least square method can be adopted. By modifying equation (4), the sum of squares can be formed as displayed in equation (6), including each possible combination of two arbitrary points. Then, the minimum two coefficients of this equation can be regarded as k and µ:

However, the application of least square method analysis has some limitations. First of all, it is necessary to assume that the random error term is independent and obeys the distribution of zero mean and same variance, which is often difficult to meet the requirements of field test data. Besides, it is essentially based on the conditional mean model, mainly in the framework of the mean statistical analysis, ignoring the information of the distribution tail. The result of being too mean is to lose a lot of other information. Hence, Bayesian quantile regression method [19, 20] can be adopted to estimate the parameters k and µ.

3.1. Quantile Regression

Since quantile regression was proposed [23], its theoretical research and practical application have been developed rapidly [2427]. The calculation model is based on a combined conditional quantile function model. It can characterize the influence of the quantiles of the dependent variable, which covers for the shortcomings of only focusing on the mean value.

For a given xi, condition P quantile of yi iswhere is a coefficient vector that only depends on P. The P-regression quantile can be defined as an arbitrary solution of the quantile regression minimization lower equation problem, which is presented aswhere

Equation (9) is defined as loss function, as known as “kernel function”.

3.2. Asymmetric Laplace Distribution (ALD)

Research results illustrate that the traditional symmetric Laplace laws are a subset of ALD [28]. The application of ALD in quantile regression has been supplemented and improved in theory [20, 2830].

If a random variable Yσ,µ is considered to obey an asymmetric Laplacian distribution, it should satisfy the following conditions; the position parameter , the scale parameter , and the characteristic function has the following form:and its distribution is followed as Yσ,µALD(σ, µ).

Maximum likelihood estimation is one of the most popular methods in statistical inference. For equation (4), assuming that the error of the model obeys normal distribution and the mean value is 0, the standard deviation σ is known, i.e., ε∼N(0, σ2), the likelihood function based on parameter β of sample is as follows:

Maximization of likelihood function can be achieved by least square estimation. If the probability density function of the error term ε of the model is assumed to bewhere is the loss function. Based on this, considering the asymmetric Laplace's standard density function, the form is as follows:

Minimizing the loss function equation (8) is equivalent to maximizing the likelihood function of the asymmetric Laplacian density of the independent joint distribution. If the random variable µ conforms to an asymmetric Laplace distribution, its probability density satisfies the formula as follows:

3.3. Bayesian Quantile Regression

At present, a very common phenomenon is that Bayesian statistical inference is applied to quantile regression [19, 20, 2931]. Scholars [3235] have developed it into a systematic theory and method of statistical inference and applied it in various fields [36, 37]. Bayesian analysis method introduces the prior information of relevant parameters to improve the effect of statistical inference. And the posterior distribution of the parameters obtained is compared with the parameters estimated by the traditional quantile regression regarding parameters as random variables.

Assuming that the model error obeys an asymmetric Laplace distribution, the basic idea of Bayesian quantile regression is that the posterior density of the parameters is finally obtained by (1) transforming the problem of minimizing the loss function into a likelihood function that maximizes the asymmetric Laplace density of the independent joint distribution of the error terms and (2) assuming the prior distribution of the parameters and adopting the appropriate sampling algorithm.

If conditional quantiles are to be computed, no matter what the initial distribution of data is, two assumptions need to be satisfied in the generalized linear model:(1)is asymmetric Laplace density.(2)For any , the connection function is . For a given observation value , the posterior distribution can be presented aswhere is a prior distribution of β, is a likelihood function, which can be written by the Bayesian formula and a position parameter :

Any prior distribution can be applied to equation (16), and Bayesian quantile regression transforms the solution of the original model into the solution of the likelihood function of equation (16).

3.4. MCMC Method

The advent of Markov chain Monte Carlo (MCMC) methods [3840] largely solves the problem of computing the high dimensional integrals, which is necessary to characterize the posterior distributions for most models. To begin with, a Markov chain is constructed by sampling. The limiting distribution of the Markov chain is the actual distribution f(θ) of a parameter θ. With the convergent Markov chain, parameters of the quantile model are then estimated. The Gibbs sampling, which is one of the most popular MCMC algorithms in the Bayesian analysis [41, 42], is adopted herein to facilitate the calculation.

The parameter estimation process is considered by the principle of Gibbs sampling algorithm, and the estimated parameters are k and µ. Let X represent the available data set and M represent the accepted model, supposing that the following two conditional distributions are known as and , where X is obtained from full-scale field tests.

Let µ0 be any initial value of µ, and the process of Gibbs sampling is as follows:Step 1: fromStep 2: from

This completes a Gibbs iteration, and the parameters change to and . Using the new parameter as the initial value, repeat the previous random sampling iteration, thus completing another Gibbs iteration process. After repeating for m times, a series of random samples are obtained: .

By repeating n large enough iterations and discarding the previous m random samples, Gibbs samples will be obtained: .

4. Full-Scale Test for Friction Coefficients

4.1. Test Design

In order to study the friction characteristics of prestressing loss, four prefabricated box girders with span of 20 m, 25 m, and 30 m (two box girders) were selected as test beams. The strength design grade of structural concrete was C50. The 15.2-1860 standard prestressing strands were adopted as steel strands, whose tension control stress under anchor and outside anchor was 1350 MPa and1437 MPa, respectively, and the elastic modulus of prestressing tendon was Ep = 1.95∗105 MPa. The prestressing force was symmetrically distributed along the vertical central line of the cross section of the beam.

Four beams were chosen as samples for the full-scale field tests to monitor tendons. Two web bundles (N1 and N3) and one floor bundle (N4 for 20 m and N5 for 25 m and 30 m) were deployed with sensors (anchor cable dynamometer and EM sensors) in each test beam. The cross sections of prestressed reinforcement distribution in the field test beam are described in Figure 1.

4.2. Layout of Measuring Points

The anchor cable dynamometers were embedded in the anchorage end of each test prestressing tendon so as to monitor the true stress under anchorage, while the electromagnetic sensors were embedded in different positions of each test prestressing tendon. Figure 2 outlines the layout of measuring points, and the parameters (a to h) of diverse test beams are shown in Table 2.

4.3. Anchor Cable Dynamometer

The anchor cable dynamometer is a kind of piercing multistring pressure sensor, which can be used for cable force measurement and construction control during cable tension. Moreover, cable dynamometer may be applied for Jack force calibration when tensioning prestressing steel bar. Apart from that, it is well-adapted to long-term monitoring and automatic monitoring. Cable dynamometer is an intelligent temperature type, measuring results of which are calibrated according to the measured temperature. Six-string intelligent anchor cable dynamometer is used in this test, as shown in Figure 3.

4.4. Electromagnetic (EM) Sensor

The electromagnetic sensor is mainly composed of primary and secondary coils, where a temperature sensor is built-in for temperature compensation. Figure 4 indicates the structure sketch of the electromagnetic sensor. When the primary coil passes into the pulse current, the secondary coil generates instantaneous current, thus obtaining the induced voltage. The induced voltage magnitude depends on the permeability of the specimen, which may exert an influence on the stress of the specimen. Therefore, the relationship between the induced voltage and the stress can be measured. Electromagnetic sensor was applied to PSC girders, and the sensor and reader are shown in Figure 5.

5. Results and Discussion

The experiments obtained values as parameters from equation (2) are shown in Table 3, where the value of tension control stress is fpj, the value of friction loss stress is ΔfpF, the value of distance is x, and the value of sag ratio is θ. Twenty sets of data of the first three beams (beam = 20 m, 25 m, 30 m) were regressed by least square method and Bayesian quantile, respectively, and the parameters were tested by five sets of data of the fourth beam.

Figure 6 only reveals the sample traces of the estimated parameters at the quantile of 0.5 because they can represent the rest of figure. A standard normal distribution was selected as its prior distribution for the estimated parameters [33], and the Gibbs sampling algorithm was adopted. To ensure that the chain was stable and the results were more accurate, the number of iterations was set to 12000.

In order to check the convergence of the MCMC chain, the posterior distribution of the edges can be visually displayed by drawing a histogram of the sampling parameters (µ and k). The histograms of the posterior distribution of regression parameters are overlaid with the density of the normal approximation of the posterior distribution as shown in Figure 7. Two thousand draws were excluded as burn-in draws, aiming to make sure only draws from the stationary posterior distribution are plotted.

The Bayesian quantile estimates and least squares estimates of the estimated parameters at different quantiles ( = 0.05, 0.25, 0.5, 0.75, 0.95) are listed in Table 4. And the simulations were performed on the programming package R software.

Figure 8 gives information of the quantile plots for the k and μ variables, plotted according to data in Table 4. There is a negative impact when μ is below the quantile of 0.25 whereas a positive impact when μ is above the quantile of 0.25, and this effect becomes more pronounced for higher quantiles. And the estimated value obtained by the OLS method is larger than the Bayesian median. The expected value of μ for the overall quantile is dramatically lower than the mean value. The variable k responds insignificantly to the expected median, but the quantile plot shows that it has a negative effect on lower quantiles and a positive effect on higher quantiles, indicating that k should not be ignored. And both the estimated value calculated by the OLS method and Bayesian median are steadily close to zero, which means that, for the average of the predicted values, the expected value of k in the lower quantile is greatly lower than the average.

Figure 9 demonstrates the range of friction coefficients in the four codes and the friction coefficients obtained by regression. The results show that the friction coefficients obtained in the actual PSC girder bridge are different from those determined in the tests, and the full-scale tests contain many uncertainties. Specifically, all of the regulations of curvature friction coefficient values in the specifications are range values. The curvature friction coefficient values obtained by the least square method are larger than those given by the four specifications, while the µ values obtained by Bayesian quantile regression approach the lower limit in JTG 3362-2018, in the range given by ACI 318-08, ACI 318-14, and AASHTO specifications. On the other hand, the wobble friction coefficient values obtained by both methods are in the middle of ACI 318-08 [13]. To be more exact, the wobble friction coefficient values obtained by Bayesian regression are larger than those obtained by least squares method, closer to the middle, but still larger than those specified in the other three specifications. Generally speaking, the values given by ACI 318-08 are relatively conservative and have a wider range for the friction coefficient. The friction coefficients obtained by different calculation methods are quite different, so it is difficult to give a single fixed value.

The last five sets of data in Table 3 were used to test the sample mean hypothesis of the regression parameters obtained by the two methods. The test statistics of the two methods were compared with each other, and the results are shown in Table 5.

From the data in the table, under the significance level of α = 0.05, there was no obvious difference between Bayesian quantile regression predicted value and test value, whilst the least squares predicted value and test value presented an enormous difference. Therefore, the Bayesian quantile regression parameters are more reasonable than the least squares regression parameters.

6. Conclusion

Four prefabricated prestressed box girder bridges with spans of 20, 25, and 30 meters were adopted for full-scale field tests, and the friction coefficient was investigated, which was very rare in previous relevant studies. The least square method and Bayesian quantile regression method were used to analyze the actual friction coefficient data sets of real bridges. Field test results indicated that there was a significant difference between the prediction value and the test value of the least square method. Nevertheless, the predicted value and the test value of Bayesian quantile regression showed no obvious difference. The value of curvature friction coefficient obtained by least squares regression was 0.308/rad, larger than the specified value, so it was not adopted. However, the value of curvature friction coefficient obtained by Bayesian quantile regression was 0.197/rad, which was similar to the median value recommended by ACI 318-08 and could be adopted. The wobble friction coefficient values obtained by regression of the two methods were 0.004/m and 0.005/m, respectively. Besides, the regression values were within the scope of ACI 318-08 and were larger than those recommended in other specifications. As a consequence, in the practical engineering, when the number of samples is small, the Bayesian quantile algorithm proposed in this paper should be used to obtain the friction coefficient for the purpose of higher accuracy.

Data Availability

The data supporting the conclusions of the study can be obtained in the manuscript.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The authors would like to acknowledge the financial support from Chinese Universities Scientific Fund no. 300102218516 and Guangdong Provincial Transportation and Transportation Office Technology Funding Project (2016-02-016).