Modeling Schemes and Performance Comparisons of Unbonded and Partially Bonded Tendon in Post-Tensioned Concrete Beam

: The modeling method of unbonded effects is a challenging and hot topic for the structural performance analysis of unbonded and partially bonded post-tensioned concrete beams. The main concerns accounting for the unbonded effects are the longitudinal free-slip behaviors and the vertical deformation compatibility relationship between the unbonded tendon and concrete beam. Three modeling schemes, namely, the beam–truss element model, the slipping cable element model


Introduction
The unbonded prestressing tendon has become a popular option for constructing post-tensioned concrete beam due to its significant advantages of facilitating site operation, replaceability, and preferable deformability [1,2].The main types of unbonded tendons include internal unbonded and externally prestressing tendons [3].In particular, the unbonded tendon elongation can no longer be calculated based on the plane-section assumption, but is related to the overall deformation of the concrete beam between the anchor points; this is generally called the unbonded phenomenon.The unbonded effects cause quite different structural performance, as compared to bonded members [4,5], which presents challenges for structural analysis and design [6,7].
Over the past few decades, numerous experimental studies have explored the mechanical properties of unbonded prestressed concrete beams [8][9][10][11].For the uncracked beam, the influence of the unbonded effects on the structural stiffness is small enough to ignore, because the tendon stress increment is tiny [12].Once the concrete beam cracks, the tendon elongation will significantly develop along with the structural deformation.
In general, two main issues still need to be clarified for the modeling of unbonded tendons.First, the "slipping cable element" model and the "slack spring" model are widely used schemes for modeling unbonded tendons, but the associated simulation accuracy is still not well compared and discussed.Second, the "slack spring" model is currently only used for the internal unbonded straight tendon; whether it can be popularized for curved or draped tendons, or in fully or partially unbonded conditions, should be studied.To address these two main issues, three analyses are conducted in this paper.First, the modeling schemes, including the beam-truss element, slipping cable element models, and slack spring models, are systematically reviewed and compared for applicability.Second, a detailed theoretical analysis is conducted for the "slack spring" method, and error sources and suggestions for improvement are proposed for the simulation of unbonded tendons with different profiles and bonding conditions.Further, with the suggested slack spring method, the influences of unbonded length on the flexural performance are studied by parameter analysis.Some interesting influence rules for unbonded length as to deformability, carrying capacity, and self-centering performance are concluded.

Modeling Schemes for Unbonded and Partially Bonded Tendons 2.1. Mechanical Characteristics of Post-Tensioned Concrete Beams with Different Bonding Conditions
The bonding conditions in post-tensioned concrete beams include completely bonded, unbonded, and partially bonded members, as shown in Figure 1.

Modeling Schemes and Applicability Comparisons
From the view of computational efficiency, the fiber beam element-based model is usually considered more appropriate than the solid element model for the flexural performance analysis of prestressed concrete beams [39,40].The basic modeling idea based on the nonlinear fiber beam element is shown in Figure 2.For the loading case in Figure 3a, the representative modeling schemes (semi-structural analysis model) for different bonding conditions are as depicted in Figure 3b-d.The main differences between these models for different bonding conditions are the adopted element types for tendons and the connections for concrete beams and tendons.(1) Model I is the conventional bonded prestressed concrete beams model, as shown in Figure 3b.The concrete beam is modeled by the fiber beam elements, and the tendons are modeled by the two-node truss elements.The concrete beam and tendons are meshed For a fully bonded prestressed concrete beam, the grouting in the duct makes the deformation between the tendon and the surrounding concrete compatible, as shown in Figure 1a.The axial and rotational deformation and structural deflection of the concrete beam will cause tendon elongation.The sectional deformation can satisfy the plane-section assumption, showing better bearing capacity and ductility.
For unbonded members (including internally unbonded prestressed members and externally prestressed members), the compatible deformation condition is only satisfied at the anchorage point, as shown in Figure 1b.The tendon's vertical deformation develops along with the structural deflection in the unbonded area, but there is no connection in the axial direction.Usually, the friction between tendon and concrete is slight, and the free-slip assumption is adopted.In this situation, the tendon stress is equalized between the anchorage points.The tendon elongation depends on the whole structural deformation rather than only sectional analyses like those of bonded members.The modeling method used for the unbonded effects of tendons is a hot topic, and will be further discussed in this paper.
For partially bonded members, the bonded and unbonded areas co-exist in the ducts, as shown in Figure 1c.This situation mainly corresponds to the bonding defects caused by tendon corrosion [24,38] or artificial construction seeking to improve structural performance [25].The mechanical behaviors of these members are more complicated when considering the bonded and unbonded conditions simultaneously.To the authors' knowledge, there are still very few studies on modeling methods and performance comparisons for partially bonded concrete beams, topics which will be discussed in this paper.

Modeling Schemes and Applicability Comparisons
From the view of computational efficiency, the fiber beam element-based model is usually considered more appropriate than the solid element model for the flexural performance analysis of prestressed concrete beams [39,40].The basic modeling idea based on the nonlinear fiber beam element is shown in Figure 2.For the loading case in Figure 3a, the representative modeling schemes (semi-structural analysis model) for different bonding conditions are as depicted in Figure 3b-d.The main differences between these models for different bonding conditions are the adopted element types for tendons and the connections for concrete beams and tendons.

Modeling Schemes and Applicability Comparisons
From the view of computational efficiency, the fiber beam element-based model is usually considered more appropriate than the solid element model for the flexural performance analysis of prestressed concrete beams [39,40].The basic modeling idea based on the nonlinear fiber beam element is shown in Figure 2.For the loading case in Figure 3a, the representative modeling schemes (semi-structural analysis model) for different bonding conditions are as depicted in Figure 3b-d.The main differences between these models for different bonding conditions are the adopted element types for tendons and the connections for concrete beams and tendons.(1) Model I is the conventional bonded prestressed concrete beams model, as shown in Figure 3b.The concrete beam is modeled by the fiber beam elements, and the tendons are modeled by the two-node truss elements.The concrete beam and tendons are meshed to ensure they have the same longitudinal coordinates, and the rigid beams are built between them to model the connections.With the rigid beam, the tendon elongation will only be controlled by the deformation of the connected beam element, and the stress of the adjacent tendon elements will be independent.Therefore, the unbonded effects cannot be considered in Model I.
(2) In Model II, the new element type named the "slipping cable element" is used to model the fully unbonded tendons.For unbonded members, the critical problem is to model the constant stress mechanical characteristics within the whole unbonded tendon length.This element type takes the tendon between anchorage points as an entire element, one which consists of multiple segments and nodes, as shown in Figure 3c.The element elongation will be calculated during each iterative step to ensure the same stress within the tendon.The finite element formulas have been reported in references [33][34][35].This method has been verified to be accurate for modeling unbonded tendons.However, this new element type has not been developed in the general finite element software, and the self-compiled program is needed [41].Moreover, the downside is that this method can only be used for fully unbonded members, not partially bonded members, due to the coexisting bonded and unbonded conditions beyond the element assumption.(1) Model I is the conventional bonded prestressed concrete beams model, as shown in Figure 3b.The concrete beam is modeled by the fiber beam elements, and the tendons are modeled by the two-node truss elements.The concrete beam and tendons are meshed to ensure they have the same longitudinal coordinates, and the rigid beams are built between them to model the connections.With the rigid beam, the tendon elongation will only be controlled by the deformation of the connected beam element, and the stress of the adjacent tendon elements will be independent.Therefore, the unbonded effects cannot be considered in Model I.
(2) In Model II, the new element type named the "slipping cable element" is used to model the fully unbonded tendons.For unbonded members, the critical problem is to model the constant stress mechanical characteristics within the whole unbonded tendon length.This element type takes the tendon between anchorage points as an entire element, one which consists of multiple segments and nodes, as shown in Figure 3c.The element elongation will be calculated during each iterative step to ensure the same stress within the tendon.The finite element formulas have been reported in references [33][34][35].This method has been verified to be accurate for modeling unbonded tendons.However, this new element type has not been developed in the general finite element software, and the self-compiled program is needed [41].Moreover, the downside is that this method can only be used for fully unbonded members, not partially bonded members, due to the co-existing bonded and unbonded conditions beyond the element assumption.
(3) Model III is an improvement of Model I.The tendons are also modeled by the traditional truss element.However, the additional auxiliary nodes were built at the location of the tendon centroid to model the surrounding concrete.The rigid beam elements were set between concrete beam nodes and auxiliary nodes to transform the structural deformation to the surrounding concrete.The connections between the tendon nodes and auxiliary nodes were coupled for vertical degree-of-freedom (DOF) (k v equals a considerable value) and set to zero-stiffness spring (k t = 0) in the longitudinal direction (called "slack spring method" in this paper), as shown in Figure 3d.In this way, the stresses in the adjacent truss elements can be realized as approximately equal.Further, if we set k t to be a considerable value, the model can also be applied to the analysis of bonded members.Chen et al. [37] applied this modeling method for the internal unbonded tendons using the existing element types in OpenSees software [42].This is a worthwhile method, because it can be realized by the built-in element types in most finite element software.However, there are still two issues that need to be clarified.First, Chen et al. [37] only analyzed the straight unbonded prestressed members.It is necessary to discuss the applicability to more complicated tendon profiles, such as curvilinear internal unbonded and draped external tendons.Then, it needs to be determined whether the "slack spring method" can be popularized to analyze the behaviors of partially bonded members by adjusting the different values of k t .This paper mainly contributes to researching and validating these two problems and further discusses the effects of different bonding conditions on structural performance.auxiliary nodes were coupled for vertical degree-of-freedom (DOF) (kv equals a considerable value) and set to zero-stiffness spring (kt = 0) in the longitudinal direction (called "slack spring method" in this paper), as shown in Figure 3d.In this way, the stresses in the adjacent truss elements can be realized as approximately equal.Further, if we set kt to be a considerable value, the model can also be applied to the analysis of bonded members.Chen et al. [37] applied this modeling method for the internal unbonded tendons using the existing element types in OpenSees software [42].This is a worthwhile method, because it can be realized by the built-in element types in most finite element software.However, there are still two issues that need to be clarified.First, Chen et al. [37] only analyzed the straight unbonded prestressed members.It is necessary to discuss the applicability to more complicated tendon profiles, such as curvilinear internal unbonded and draped external tendons.Then, it needs to be determined whether the "slack spring method" can be popularized to analyze the behaviors of partially bonded members by adjusting the different values of kt.This paper mainly contributes to researching and validating these two problems and further discusses the effects of different bonding conditions on structural performance.For Model III, taking a tendon node as an analytic target, the force diagram can be drawn as in Figure 4.The static equilibrium equations built in x and y directions are shown in Equation (1).The combination of two equations yields Equation (2).If unbonded effects are satisfied, the condition T 1 = T 2 ̸ = 0 must be satisfied.Equation (2) can be transformed into Equation (3).This is simple but valuable.The prerequisite for the unbonded effect to be satisfied by Model III is that the direction of the zero-length element should be on the bisector line of the angle between the truss elements.For the straight unbonded tendon, the direction obtained by Equation (3) is vertical.However, for the tendons with curved or draped profiles, the suggested direction of the zero-length element should be determined according to the slope of the tendon profile rather than the descriptions adopted by Chen et al. [37].To verify this conclusion, we build two models for Model III, named Model III-1 (zero length element along with the vertical direction) and Model III-2 (zero length element along with the angular bisector), respectively.
or draped profiles, the suggested direction of the zero-length element should be determined according to the slope of the tendon profile rather than the descriptions adopted by Chen et al. [37].To verify this conclusion, we build two models for Model III, named Model III-1 (zero length element along with the vertical direction) and Model III-2 (zero length element along with the angular bisector), respectively.
0, sin sin sin 0 0, cos cos cos 0 ( )  However, we noticed that the tendon's dip angle and structural deformation will constantly change under loading.So, even though we define the zero-length element by the direction of Equation ( 3), the angular deviation δ will inevitably appear following the beam deflection.Defining the coefficient β = 1 − (T1/T2) describes the degree of non-uniform distribution between two truss elements.Substituting θ = (α1 + α2)/2 + δ into Equation (2), Equation ( 4) can be obtained.With the increase of δ caused by structural deformation, the nonuniformity relative to the adjacent truss element will increase.The suggested element direction in Equation (3) can decrease the initial deviation.In addition, the angular deviation δ is mainly caused by structural deflections, which is small compared with (α1 + α2)/2, for purposes of practical engineering.The effects of δ will be analyzed in the following section.
Table 1 shows the summary and comparisons between different modeling schemes.It can be seen that Model III is a promising method for members with different bonding conditions.The model validation and detailed discussion are presented in the following section.However, we noticed that the tendon's dip angle and structural deformation will constantly change under loading.So, even though we define the zero-length element by the direction of Equation ( 3), the angular deviation δ will inevitably appear following the beam deflection.Defining the coefficient β = 1 − (T 1 /T 2 ) describes the degree of non-uniform distribution between two truss elements.Substituting θ = (α 1 + α 2 )/2 + δ into Equation (2), Equation ( 4) can be obtained.With the increase of δ caused by structural deformation, the nonuniformity relative to the adjacent truss element will increase.The suggested element direction in Equation (3) can decrease the initial deviation.In addition, the angular deviation δ is mainly caused by structural deflections, which is small compared with (α 1 + α 2 )/2, for purposes of practical engineering.The effects of δ will be analyzed in the following section.
Table 1 shows the summary and comparisons between different modeling schemes.It can be seen that Model III is a promising method for members with different bonding conditions.The model validation and detailed discussion are presented in the following section.

Computational Procedure
All of the analysis models are built in OpenSees V2.5.0 software [42] because of its open-source advantages and powerful nonlinear analysis capability.In particular, the slipping cable element in Model II is built employing the program pre-developed by the authors in OpenSees v2.5.0 [32,33].The other models are all built using the built-in element type in OpenSees v2.5.0 software.In the numerical model, the concrete beam element size is 50 mm, and the tendon element is meshed with the same element amounts and nodal coordinate in the longitudinal direction as the concrete beam.
The materials used in the analysis mode include concrete, steel reinforcement, and prestressing tendon.The concrete is simulated by the "Concrete 02" model [43].The "Rein-forcingSteel" model is employed to simulate the material nonlinearity of steel reinforcement according to the suggestion of reference [44].The "Steel 02" model is adopted to simulate steel tendons [33].
The numerical models are loaded from two analysis steps.First, the prestressing force is applied to the tendon by way of initial stress according to the control tension stress.Then, externally applied forces are loaded on the model by displacement-controlled loading up to failure.

Externally Prestressed Concrete Beam
The specimen named SSS1, in the experiment conducted by Au et al. [45], is taken as an example in order to verify the feasibility of the mentioned modeling schemes for the unbonded behavior simulation of the external tendon.The specimen was a simply supported concrete beam with a span length of 4.5 m.The structural design details and loading schema are shown in Figure 5.The material characteristics of concrete, rebar, and tendon are listed in Table 2.The external tendons were designed with a harped profile, the effective tendon depth d p of which was 112 mm at the beam end and 168.8 mm at the midspan.The external tendon was post-tensioned at 28 days of age, then monotonically loaded up to failure on the same day.Under four-point loading, the failure is caused by the tensile reinforcement yield and the top flange's concrete crushing at the midspan cross-section.Three different models mentioned in Section 2 are employed to predict the whole-process structural behavior.

Computational Procedure
All of the analysis models are built in OpenSees V2.5.0 software [42] because of its open-source advantages and powerful nonlinear analysis capability.In particular, the slipping cable element in Model II is built employing the program pre-developed by the authors in OpenSees v2.5.0 [32,33].The other models are all built using the built-in element type in OpenSees v2.5.0 software.In the numerical model, the concrete beam element size is 50 mm, and the tendon element is meshed with the same element amounts and nodal coordinate in the longitudinal direction as the concrete beam.
The materials used in the analysis mode include concrete, steel reinforcement, and prestressing tendon.The concrete is simulated by the "Concrete 02" model [43].The "Re-inforcingSteel" model is employed to simulate the material nonlinearity of steel reinforcement according to the suggestion of reference [44].The "Steel 02" model is adopted to simulate steel tendons [33].
The numerical models are loaded from two analysis steps.First, the prestressing force is applied to the tendon by way of initial stress according to the control tension stress.Then, externally applied forces are loaded on the model by displacement-controlled loading up to failure.

Externally Prestressed Concrete Beam
The specimen named SSS1, in the experiment conducted by Au et al. [45], is taken as an example in order to verify the feasibility of the mentioned modeling schemes for the unbonded behavior simulation of the external tendon.The specimen was a simply supported concrete beam with a span length of 4.5 m.The structural design details and loading schema are shown in Figure 5.The material characteristics of concrete, rebar, and tendon are listed in Table 2.The external tendons were designed with a harped profile, the effective tendon depth dp of which was 112 mm at the beam end and 168.8 mm at the midspan.The external tendon was post-tensioned at 28 days of age, then monotonically loaded up to failure on the same day.Under four-point loading, the failure is caused by the tensile reinforcement yield and the top flange's concrete crushing at the midspan cross-section.Three different models mentioned in Section 2 are employed to predict the whole-process structural behavior.Figure 6 compares the predictions and the test results of load-deflection and deflectiontendon stress increment curves.The comparisons clearly show the noteworthy influences of the slip effect of the external tendon at the deviators on flexural performance.For Model I, the slip effects at the deviators are not considered.The tendon elongation is independent of the different tendon segments, which is subtly different from the actual behavior of the external tendon.The tendon elongation of Model I mainly occurred at the middle tendon segment, "Seg-2", causing the higher tendon stress increment and structural flexural capacity.Model I overestimates the structural stiffness and flexural capacity after the mild rebar yielding, without considering the slip effects.For Model II, the external tendon is taken as a whole element, and the slip effects are naturally considered during the deriving of the finite element formula.On the basis of Model I, the longitudinal relaxation spring determined by Equation ( 3) is introduced at the deviators in Model III to consider the slip effects.The tension redistribution phenomenon among the tendon segments will occur in Model III due to the zero "slack spring" stiffness in the longitudinal direction.In equilibrium, the tendon stress is approximately equal anywhere within the tendon length, satisfying the assumption of unbonded effects.Model II and Model III analysis results are very close and agree well with the tests for flexural capacity and tendon stress increments.
Model I, the slip effects at the deviators are not considered.The tendon elongation is independent of the different tendon segments, which is subtly different from the actual behavior of the external tendon.The tendon elongation of Model I mainly occurred at the middle tendon segment, "Seg-2", causing the higher tendon stress increment and structural flexural capacity.Model I overestimates the structural stiffness and flexural capacity after the mild rebar yielding, without considering the slip effects.For Model II, the external tendon is taken as a whole element, and the slip effects are naturally considered during the deriving of the finite element formula.On the basis of Model I, the longitudinal relaxation spring determined by Equation ( 3) is introduced at the deviators in Model III to consider the slip effects.The tension redistribution phenomenon among the tendon segments will occur in Model III due to the zero "slack spring" stiffness in the longitudinal direction.In equilibrium, the tendon stress is approximately equal anywhere within the tendon length, satisfying the assumption of unbonded effects.Model II and Model III analysis results are very close and agree well with the tests for flexural capacity and tendon stress increments.Further, let us increase the midspan dp to 1/10 span length (dp = 450 mm), and the dip angle θ of the draped tendon to 35.3°.For this situation in which the included angle between the tendon segments is large, the orientation of the relaxation spring at the deviators in Model III has non-negligible effects on the flexural behaviors.Three numerical models are built and analyzed for comparisons, including Model II with slipping cable element, Model III with longitudinal relaxation spring along the beam length (named Model III-1), and Model II with tangential relaxation spring (named Model III-2).The analysis results are plotted in Figure 7 for comparison.It can be seen that the orientation of the relaxation spring will influence the results of the analysis, and merits attention.Model III-1 slightly underestimates the flexural capacity, by 3%.The reason for the deviation in Model II-1 is that the relaxation spring orientation is in the longitudinal direction, in which the horizontal component of Seg-1 equalizes the tendon force of Seg-2.So the tendon force of the middle segment of Seg-1 is smaller than that of Seg-2, causing a Further, let us increase the midspan d p to 1/10 span length (d p = 450 mm), and the dip angle θ of the draped tendon to 35.3 • .For this situation in which the included angle between the tendon segments is large, the orientation of the relaxation spring at the deviators in Model III has non-negligible effects on the flexural behaviors.Three numerical models are built and analyzed for comparisons, including Model II with slipping cable element, Model III with longitudinal relaxation spring along the beam length (named Model III-1), and Model II with tangential relaxation spring (named Model III-2).The analysis results are plotted in Figure 7 for comparison.It can be seen that the orientation of the relaxation spring will influence the results of the analysis, and merits attention.Model III-1 slightly underestimates the flexural capacity, by 3%.The reason for the deviation in Model II-1 is that the relaxation spring orientation is in the longitudinal direction, in which the horizontal component of Seg-1 equalizes the tendon force of Seg-2.So the tendon force of the middle segment of Seg-1 is smaller than that of Seg-2, causing a prediction of lower flexural capacity.With the suggested tangential relaxation orientation, which is perpendicular to the angle bisector of tendon segments (Model III-2), the analysis results are close to the slipping cable element model (Model II).It should be noted that the method exemplified by Model III-2 can be easily realized in the general finite element software without additional program development, like that of Model II.
prediction of lower flexural capacity.With the suggested tangential relaxation orientation, which is perpendicular to the angle bisector of tendon segments (Model III-2), the analysis results are close to the slipping cable element model (Model II).It should be noted that the method exemplified by Model III-2 can be easily realized in the general finite element software without additional program development, like that of Model II.

Internal Unbonded Prestressed Concrete Beam
A 4.8 m + 4.8 m continuous test beam with internal unbonded tendons, named YLA2 [46], is selected to verify and determine the applicability of the above-mentioned methods for modeling the internal unbonded tendons.The details of the structural dimensions, cross-section, and loading pattern are shown in Figure 8, and the related material parameters are listed in Table 3.The diameter of the steel reinforcement used was 12 mm, and its elastic modulus and yield strengths were 361 MPa and 200 GPa, respectively.The prestressing tendon's yield strength, ultimate strength, and elastic modulus are 1680 MPa, 1941 MPa, and 197 GPa, respectively.The tendon profile was designed using the piecewise quadratic equations [46].The relation between the beam deformation and the tendon elongation was more complicated than that of the straight tendon profile.The internal unbonded tendon was prestressed at 28 days of age, and then third-point loading was applied at each span up to failure.With the increasing load, the top mild rebar As2 at the middle support cross-section was first-yielded, and then flexural failure occurred at the midspan cross-section almost simultaneously.The tendon was still elastic at the ultimate, which differs from those of bonded members.The analysis results for load-deflection and load-tendon stress increment as determined by different modeling schemes, and their comparisons with the test results, can be seen in Figure 9.The method of Model I equivalently takes the specimen as bonded members, overestimating the flexural capacity and the tendon stress.The results of Model II and Model III-2, with the unbonded effects considered, are matched well with the tests for both flexural capacity and tendon stress increment.The slight deviations between predictions and tests are mainly caused by the possible friction and material randomness.

Internal Unbonded Prestressed Concrete Beam
A 4.8 m + 4.8 m continuous test beam with internal unbonded tendons, named YLA2 [46], is selected to verify and determine the applicability of the above-mentioned methods for modeling the internal unbonded tendons.The details of the structural dimensions, cross-section, and loading pattern are shown in Figure 8, and the related material parameters are listed in Table 3.The diameter of the steel reinforcement used was 12 mm, and its elastic modulus and yield strengths were 361 MPa and 200 GPa, respectively.The prestressing tendon's yield strength, ultimate strength, and elastic modulus are 1680 MPa, 1941 MPa, and 197 GPa, respectively.The tendon profile was designed using the piecewise quadratic equations [46].The relation between the beam deformation and the tendon elongation was more complicated than that of the straight tendon profile.The internal unbonded tendon was prestressed at 28 days of age, and then third-point loading was applied at each span up to failure.With the increasing load, the top mild rebar A s2 at the middle support cross-section was first-yielded, and then flexural failure occurred at the midspan cross-section almost simultaneously.The tendon was still elastic at the ultimate, which differs from those of bonded members.The analysis results for load-deflection and load-tendon stress increment as determined by different modeling schemes, and their comparisons with the test results, can be seen in Figure 9.The method of Model I equivalently takes the specimen as bonded members, overestimating the flexural capacity and the tendon stress.The results of Model II and Model III-2, with the unbonded effects considered, are matched well with the tests for both flexural capacity and tendon stress increment.The slight deviations between predictions and tests are mainly caused by the possible friction and material randomness.The effects of the orientation of the relaxation springs in Model III on the tendon stress distributions are further discussed here.For comparison purposes, we collect the tendon stress results from different models under the same midspan deflection of 40 mm.The discussed models include Model II with slipping cable element, Model III-1 with relaxation springs defined along the beam length direction, and Model III-2 with relaxation springs defined perpendicular to the angle bisector of tendon segments.The comparisons of the results are shown in Figure 10.The tendon stress in Model II is uniformly distributed within the whole element, which fully satisfies the frictionless and free-slip hypothesis.For Model III, there are some differences in the tendon stress along the beam length.The standard deviation of the tendon stress by the Model III-1 method is 5.54 MPa, while for Model III-2 with the suggested relaxation springs orientation, the standard deviation is reduced to 1.25 MPa, which is closer to the goal hypothesis.The relaxation springs method with the suggested orientation (Model III-2) can be adopted to simulate the behavior of the internally prestressed concrete beams, and the analysis results are close to those of the slipping cable element model (Model II).The effects of the orientation of the relaxation springs in Model III on the tendon stress distributions are further discussed here.For comparison purposes, we collect the tendon stress results from different models under the same midspan deflection of 40 mm.The discussed models include Model II with slipping cable element, Model III-1 with relaxation springs defined along the beam length direction, and Model III-2 with relaxation springs defined perpendicular to the angle bisector of tendon segments.The comparisons of the results are shown in Figure 10.The tendon stress in Model II is uniformly distributed within the whole element, which fully satisfies the frictionless and free-slip hypothesis.For Model III, there are some differences in the tendon stress along the beam length.The standard deviation of the tendon stress by the Model III-1 method is 5.54 MPa, while for Model III-2 with the suggested relaxation springs orientation, the standard deviation is reduced to 1.25 MPa, which is closer to the goal hypothesis.The relaxation springs method with the suggested orientation (Model III-2) can be adopted to simulate the behavior of the internally prestressed concrete beams, and the analysis results are close to those of the slipping cable element model (Model II).

Partially Bonded Prestressed Concrete Beam
Losanno et al. [24] conducted a series of specimen tests to discuss the effects of grouting conditions on the flexural behaviors of prestressed concrete beams.Three specimens with different bonded states are collected here to verify the modeling method for partially

Partially Bonded Prestressed Concrete Beam
Losanno et al. [24] conducted a series of specimen tests to discuss the effects of grouting conditions on the flexural behaviors of prestressed concrete beams.Three specimens with different bonded states are collected here to verify the modeling method for partially bonded members, namely, B-HP (full grouting for bonded members), U-HP (no grouting for unbonded members), and PB-M-HP (no grouting at midspan for partially bonded members), respectively.The specimens were simply supported concrete beams with span lengths of 6.6 m.The design details of the specimens are shown in Figure 11, and the material properties are listed in Table 4. Tendon prestressing was imposed after 28 days of concrete curing.The actuator with 3000 kN loading capacity applied two-point loads at the midspan, spaced 850 mm apart, and all specimens were monotonically loaded to failure.The load cells were arranged at the beam ends to monitor the tendon stress variations, and the LVDT, with maximum strokes of 500 mm, was placed at the midspan to test the deflection change with increasing loads.From the above analysis cases, we can conclude that Model III-2 can be used to model the behaviors of prestressed concrete beams with different bonded conditions by setting different connection stiffness values for the relaxation springs.For bonded members, the tangential and normal stiffness of the springs can be set to a significant value (such as 1 × 10 18 N/mm, as in this paper).For unbonded members, the tangential stiffness can be set to zero for all springs.The zero-stiffness spring can be set in the unbonded area for partially unbonded members.Model I can only be used to model the bonded members, and Model II can only be used to model the unbonded members.Therefore, in this section, only Model III-2 is employed to analyze the behaviors of three specimens.
Figure 12 shows the analysis results and the comparisons with the tests.With the degradation of the bonding state (from B-HP, PB-M-HP, to U-HP), the structural flexural capacity and ultimate displacement are all decreased.The modeling method of Model III-2 can successfully capture the effects of bonding conditions on these responses.
Buildings 2024, 14, x FOR PEER REVIEW 13 of 19 unbonded members.Model I can only be used to model the bonded members, and Model II can only be used to model the unbonded members.Therefore, in this section, only Model III-2 is employed to analyze the behaviors of three specimens.
Figure 12 shows the analysis results and the comparisons with the tests.With the degradation of the bonding state (from B-HP, PB-M-HP, to U-HP), the structural flexural capacity and ultimate displacement are all decreased.The modeling method of Model III-2 can successfully capture the effects of bonding conditions on these responses.As reported in reference [24], all of the specimens displayed flexural failure at the midspan cross-section, but the failure modes had some differences, due to the different bonding conditions.A significant opening crack was observed at the midspan cross-section, dominating the failure of unbonded specimen U-HP.The concrete beam was almost As reported in reference [24], all of the specimens displayed flexural failure at the midspan cross-section, but the failure modes had some differences, due to the different bonding conditions.A significant opening crack was observed at the midspan cross-section, dominating the failure of unbonded specimen U-HP.The concrete beam was almost elastic in the other regions.For bonded specimen B-HP, the cracks were distributed in a broader range, and the crack spacing was smaller than in U-HP.The failure mode was between U-HP and B-HP for partially bonded specimen PB-M-HP, the major cracks of which appeared in the unbonded region and were slight in the other areas.
Based on the results of the analysis, the ultimate sectional curvature distribution at the peak loads is collected and compared with the experimentally observed failure modes, as shown in Figure 13.The locations of the more significant sectional curvature represent the failure cross-section with large plastic deformation.It can be clearly seen that the plastic deformation is concentrated at the midspan for unbonded specimen U-HP, and the other parts are almost elastic.The plastic deformations are distributed in a broader range for bonded specimen B-HP.Interestingly, for partially bonded specimen PB-M-HP, the sectional curvature distribution in the bonded area is close to the bonded members and close to unbonded members in the unbonded area.The results of the analyses match well with the observed failure mode from the tests, validating the effectiveness of the modeling method of Model III-2.
Buildings 2024, 14, x FOR PEER REVIEW 14 of 19 with the observed failure mode from the tests, validating the effectiveness of the modeling method of Model III-2.
Based on the validation of these prestressed concrete beams with different bonding conditions, we can conclude that Model III (slack spring method) is an effective model for predicting structural deformation, carrying capacity, and failure modes.Moreover, we suggest setting the orientation of the zero-length element on the angular bisector line between the adjacent truss elements (Model III-2).The Model III-2 can improve the prediction accuracy of tendon stress distributions, compared with Model III-1.

Discussion
Based on the validated model, the Model III-2 method is further used to discuss the effects of the degree of unbonded conditions on the flexural performance.The specimen PB-M-HP by Losanno et al. [24] is employed here as the prototype beam for numerical analysis.The structural dimensions and the design details can be seen in Figure 11.The length of the unbonded area αL is taken as the only variable to be discussed.The unbonded area length αL is set from 0 to 6600 mm, with a step size of 600 mm (α ranges from 0 to 1).The control stress for prestressing is 800 MPa.The focused performance indicia include deformability (maximum deflection), carrying capacity, and self-centering behaviors.
With the specified loading pattern in Figure 11, all models are loaded up to failure, and the load-deflection curves are shown in Figure 14a.To compare the models' self-centering behaviors, all numerical beams are firstly loaded to the same midspan deflection Based on the validation of these prestressed concrete beams with different bonding conditions, we can conclude that Model III (slack spring method) is an effective model for predicting structural deformation, carrying capacity, and failure modes.Moreover, we suggest setting the orientation of the zero-length element on the angular bisector line between the adjacent truss elements (Model III-2).The Model III-2 can improve the prediction accuracy of tendon stress distributions, compared with Model III-1.

Discussion
Based on the validated model, the Model III-2 method is further used to discuss the effects of the degree of unbonded conditions on the flexural performance.The specimen PB-M-HP by Losanno et al. [24] is employed here as the prototype beam for numerical analysis.The structural dimensions and the design details can be seen in Figure 11.The length of the unbonded area αL is taken as the only variable to be discussed.The unbonded area length αL is set from 0 to 6600 mm, with a step size of 600 mm (α ranges from 0 to 1).The control stress for prestressing is 800 MPa.The focused performance indicia include deformability (maximum deflection), carrying capacity, and self-centering behaviors.
With the specified loading pattern in Figure 11, all models are loaded up to failure, and the load-deflection curves are shown in Figure 14a.To compare the models' self-centering behaviors, all numerical beams are firstly loaded to the same midspan deflection (100 mm) and then unloaded to zero force.The residual deformation res 100 is obtained for comparisons, as shown in Figure 14b.The main indices of the results from these models, including maximum midspan deflection, peak load, and residual deformation res 100 , are listed in Table 5.The curves of these indices, along with the changing of the unbonded area length parameter α, are plotted in Figure 15.In Table 5 and Figure 15, the definition of the RES index follows the calculation of Chen et al. [37], and RES = 1 − res 100 /100.Based on the comparisons, we find some interesting conclusions: (1) Before the cracking of the concrete beam, the effects of the unbonded effects can be ignored.This result is consistent with the research conclusion of Bonopera et al. [12].Their test showed that the specimen's stiffness with a parabolic unbonded tendon and a bonded one is similar.However, the situation changes after the concrete beam cracks.The tangent slope of the load-deflection curve decreases, along with an increase of the unbonded area length.During this stage, the stiffness contributions of the tensile area are mainly from steel reinforcements and tendons.Along with the increase of unbonded area Based on the comparisons, we find some interesting conclusions: (1) Before the cracking of the concrete beam, the effects of the unbonded effects can be ignored.This result is consistent with the research conclusion of Bonopera et al. [12].Their test showed that the specimen's stiffness with a parabolic unbonded tendon and a bonded one is similar.However, the situation changes after the concrete beam cracks.The tangent slope of the load-deflection curve decreases, along with an increase of the unbonded area length.During this stage, the stiffness contributions of the tensile area are mainly from steel reinforcements and tendons.Along with the increase of unbonded area length, the tensile force increments of the tendon will decrease, causing a decrease of the tangent slope.
(2) As expected, the flexural capacity shows a decreasing trend with the increase of unbonded length, and the peak load of the unbonded member is 23.1% lower than that of the bonded member.Interestingly, we find that the variation of maximum deflection with unbonded length is not monotonic.The maximum deflection decreases for cases of α ranging from 0 to 0.4, and then increases for cases of α ranging from 0.4 to 1.0.The reason is revealed by the maximum deflection prediction equation of Peng et al. [16], as shown in Figure 16 and Equation (5).Values ε cu and c u denote the concrete ultimate compressive strain and neutral axis height, respectively.L p is the plastic hinge length of the idealized curvature distribution curve.From Equation (5), we find that the maximum deflection is mainly affected by L p and c u .Based on the previous analysis in Figure 13a, we know that L p and c u both decrease with the increase of unbonded length, causing the above-mentioned non-monotonic phenomenon.For cases with small unbonded lengths, the decrease of L p dominates the reduction of maximum deflection.With a large unbonded length, the ultimate tendon stress decrease causes the reduction of c u and increase of maximum deflection.
Buildings 2024, 14, x FOR PEER REVIEW 16 of 19 length, the tensile force increments of the tendon will decrease, causing a decrease of the tangent slope.
(2) As expected, the flexural capacity shows a decreasing trend with the increase of unbonded length, and the peak load of the unbonded member is 23.1% lower than that of the bonded member.Interestingly, we find that the variation of maximum deflection with unbonded length is not monotonic.The maximum deflection decreases for cases of α ranging from 0 to 0.4, and then increases for cases of α ranging from 0.4 to 1.0.The reason is revealed by the maximum deflection prediction equation of Peng et al. [16], as shown in Figure 16 and Equation (5).Values εcu and cu denote the concrete ultimate compressive strain and neutral axis height, respectively.Lp is the plastic hinge length of the idealized curvature distribution curve.From Equation (5), we find that the maximum deflection is mainly affected by Lp and cu.Based on the previous analysis in Figure 13a, we know that Lp and cu both decrease with the increase of unbonded length, causing the above-mentioned non-monotonic phenomenon.For cases with small unbonded lengths, the decrease of Lp dominates the reduction of maximum deflection.With a large unbonded length, the ultimate tendon stress decrease causes the reduction of cu and increase of maximum deflection.
(3) The residual deformation is mainly caused by the plastic of concrete, rebar, and (3) The residual deformation is mainly caused by the plastic of concrete, rebar, and prestressing tendon.Under the same structural deflection, the unbonded tendon stress is smaller than the bonded ones and even still elastic without plastic strain.Therefore, the self-centering performance of unbonded members is better than that of bonded members, as shown in Figure 15c.
In summary, the unbonded length significantly influences partially bonded concrete beams' deformability, carrying capacity, and self-centering behaviors.Moreover, the influence rule is not simply monotonic for the deformability.Currently, study on the partially prestressed concrete beam is still limited, due to its more complicated behaviors and construction scheme.However, perhaps we can set an appropriate unbonded length according to the balance of deformability, flexural capacity, and self-centering properties to improve the overall structural performance, a possibility which should be further studied.

Conclusions
This paper systematically presents three modeling schemes for analyzing post-tensioned concrete beams with different bonding conditions: the conventional beam-truss element model (Model I), the slipping cable element (Model II), and the "slack spring method" (Model III).The slack spring model is recommended as the generally applicable model for analyzing post-tensioned concrete beams with different bonding conditions.The direction and stiffness of the slack spring in Model III are proposed to improve the accuracy of the analysis.The applicability of these models is analyzed and compared by verifications with the experimental beams, and the effects of unbonded length on flexural performance are parametrically analyzed.The conclusions are as follows: (1) Model I is adequate for the analysis of bonded prestressed members.Model II can capture the flexural behaviors of unbonded members well, but it cannot be used to analyze partially bonded members.Moreover, the application of Model II goes beyond the function of most finite element software, and a self-complied procedure is needed.(2) Model III is applicable for analyzing prestressed concrete beams with different bonding conditions (bonded, unbonded, and partially bonded members).The prediction accuracy can be improved by setting the orientation of the zero-length element, along with the angular bisector between the adjacent truss elements (Model III-2).( 3) With increased unbonded length, the flexural capacity decreases, but the self-centering performance is improved.The effects of unbonded length on structural deformability are not monotonic, showing a trend of declining and then rising along with the increase in unbonded length.(4) An experimental study on the flexural performance of post-tensioned concrete beams with different bonding conditions needs to be conducted in the future to verify the modeling schemes thoroughly.More accurate information on the practicability of the proposed FE models can thus be provided.The influences of unbonded length on flexural capacity, deformability, and self-centering performance will be further experimentally verified.

Figure 2 .
Figure 2. Truss-fiber beam element model for the prestressed concrete beam.

Figure 1 .
Figure 1.Prestressed concrete beams with different bonding conditions.

Figure 1 .
Figure 1.Prestressed concrete beams with different bonding conditions.

Figure 2 .
Figure 2. Truss-fiber beam element model for the prestressed concrete beam.

Figure 2 .
Figure 2. Truss-fiber beam element model for the prestressed concrete beam.

Figure 3 .
Figure 3. Modeling schemes for prestressed concrete beams with different bonding conditions.Figure 3. Modeling schemes for prestressed concrete beams with different bonding conditions.

Figure 3 .
Figure 3. Modeling schemes for prestressed concrete beams with different bonding conditions.Figure 3. Modeling schemes for prestressed concrete beams with different bonding conditions.
Connection of tendon node (b) Force decomposition at tendon node

Figure 4 .
Figure 4. Connection and force decomposition of tendon node for Model III.

Figure 4 .
Figure 4. Connection and force decomposition of tendon node for Model III.
for tendon force monitoring LVDT for midspan deflection monitoring 100kN load cell for tendon force monitoring

Figure 5 .
Figure 5. Dimension details and loading schema of specimen SSS1.Figure 5. Dimension details and loading schema of specimen SSS1.

Figure 5 .
Figure 5. Dimension details and loading schema of specimen SSS1.Figure 5. Dimension details and loading schema of specimen SSS1.

Figure 6 .
Figure 6.Comparisons using the experimental results of specimen SSS1.

Figure 6 .
Figure 6.Comparisons using the experimental results of specimen SSS1.

2 Figure 7 .
Figure 7. Effects of the relaxation spring orientation on the flexural behaviors.

Figure 7 .
Figure 7. Effects of the relaxation spring orientation on the flexural behaviors.

Figure 8 .
Figure 8. Details of the specimens YLA2 and YLC2.

Figure 8 .
Figure 8. Details of the specimens YLA2 and YLC2.

Figure 9 .
Figure 9. Analysis results for comparisons between different modeling schemes.

Figure 9 .
Figure 9. Analysis results for comparisons between different modeling schemes.

Figure 10 .
Figure 10.Tendon stress distribution along the beam length under the same deflection.

Figure 10 .
Figure 10.Tendon stress distribution along the beam length under the same deflection.

425Figure 11 .
Figure 11.Design details of the specimens with different grouting conditions.

Figure 12 .
Figure 12. Analysis results of comparisons for cases with different bonded conditions.

Figure 12 .
Figure 12. Analysis results of comparisons for cases with different bonded conditions.

4 − 4 − 5 Figure 13 .
Figure 13.Failure modes comparison based on the results of the analysis.

Figure 13 .
Figure 13.Failure modes comparison based on the results of the analysis.

Figure 14 .
Figure 14.Load-deflection analysis results of numerical models with different bonding conditions.

Figure 15 .
Figure 15.Curves of maximum deflection, peak load, and RES index, along with the parameter α.

Figure 15 .
Figure 15.Curves of maximum deflection, peak load, and RES index, along with the parameter α.

Figure 16 .
Figure 16.Failure mode and idealized curvature distribution of the prestressed concrete beam.

Figure 16 .
Figure 16.Failure mode and idealized curvature distribution of the prestressed concrete beam.

Table 1 .
Comparisons between the modeling schemes for different bonding conditions.

Table 3 .
Parameters for the continuous test beams with internal unbonded tendons.

Table 3 .
Parameters for the continuous test beams with internal unbonded tendons.

Table 3 .
Parameters for the continuous test beams with internal unbonded tendons.

Table 4 .
Parameters for partially bonded prestressed concrete beam.

Table 4 .
Parameters for partially bonded prestressed concrete beam.

Table 5 .
Analysis results of the numerical models. α

Table 5 .
Analysis results of the numerical models.