Calculation of Additional Internal Forces in Post-Tensioned Prestressed Concrete Frame Columns Based on Equivalent Lateral Stiffness

: The additional internal forces in vertical members caused by prestressed tendons are typically overlooked in the design of post-tensioned prestressed concrete. A calculation method for additional internal forces in single-story multi-span prestressed concrete frame columns based on equivalent lateral stiffness is proposed in this paper. The slope-deﬂection equation for the bar element was presented using Timoshenko beam assumptions, taking into account the inﬂuence of shear and bending deformations. Subsequently, the concept of equivalent lateral stiffness and calculation equations were proposed. On this basis, the equations of the third shear and third bending moment for single-story multi-span prestressed frame columns were established. Furthermore, applying engineering examples, the method in this study was veriﬁed by ABAQUS software and previous methods. The results show that theoretical values and FEA results are in good agreement. Compared to previous methods, the method in this paper is more accurate and widely applicable. In addition, the stretching plan has a signiﬁcant path effect and time-varying effect on the interlayer distribution of the third moment. It should be considered at the building stage to check the calculation of the frame column.


Introduction
Prestressed concrete structures have long been widely used in large-scale, long-span industrial and civil buildings [1][2][3][4][5][6]. In general, prestressed tendons are arranged in horizontal members such as beams and slabs. The stretching of the prestressed tendons during post-tensioned construction results in compressive deformation of the horizontal members. Once this deformation is limited by lateral restraint, such as columns and shear walls, additional internal forces such as shears and bending moments are created in the vertical members. This internal force, in addition to the load and secondary internal force, is referred to as the "third shear" and "third bending moment" [7]. China's current code [8] does not consider the influence of this additional internal force in the design of prestressed concrete structures with lateral constraints. It is comparable to presuming that the horizontal and vertical members are in a slip connection, which is not supported by facts. Engineering problems caused by additional internal forces are not uncommon. Cracking basement walls and columns in the No.2 building of Binjiang Junyuexiangdi, Changsha, China residential community and basement shear walls in the No.2 building of Yiju Laiyincheng Changsha, China residential community are two typical examples [9]. Therefore, it is of great engineering value to research lateral restraint effects in prestressed concrete structures.
At present, most of the research on the lateral restraint effect focused on discussions of the effective preload of beams. Zhang, D et al. [10] were the first to focus on this issue and derived an equation for the effective preload of beams in the single-story prestressed frame. For multi-story prestressed frames, Jian, B [11] established equations for the effective preload of beams with different construction schemes and gave design suggestions. Subsequently, the effects of beam-column linear stiffness ratio, span number and construction scheme on the effective preload of beams were investigated by Wang, C [12], Zhang, C [13], Cao, X [14], Xia, X [15]. Furthermore, the use of lateral restraint influence coefficients to describe the effect of lateral constraints has been a popular topic [16][17][18][19][20][21]. The concept of the lateral restraint influence coefficient was first proposed by Zheng, W [16]. He established a design method for the cross-section bearing capacity of prestressed beams using the lateral restraint influence coefficient η. Subsequently, an equation for η based on the deformation compatibility method was developed and the conditions for its application were presented [17]. Moreover, Guo, Y [18] and Zhang, Q et al. [19] also obtained an equation for η based on the deformation compatibility method and the matrix displacement method. It was concluded that the deformation compatibility method was more accurate by finite element analysis (FEA). Xiong, X et al. [22], on the other hand, evaluated that using η to express the lateral constraint effect had the characteristics of clear physical meaning and easy calculation. In addition, in practical engineering the effective preload could also be obtained by establishing the overall model by structural design software [23,24].
In recent years, the influence of the lateral restraint effect on column internal force has attracted more and more attention from scholars. Zhang, D et al. [10] argued that the third bending moment was usually in the opposite direction to the bending moment caused by the external load, and that it was advantageous to disregard it for structural design. In contrast, Yan, W et al. [25] and Shen, Q [26] held the opposite view. They believed that the third bending moment should not be ignored, and it may bring potential safety hazards without consideration in design. The equation for the third bending moment applicable to single-story, single-span prestressed frames was then derived. However, the rotation angle at the column end was neglected and the calculation accuracy was poor. The existence of third bending moments was demonstrated by experimental means in an actual project [27]. In addition, it has been found that the stretching plan of prestressing tendons has a significant effect on the third bending moments in multi-story prestressed frames. However, scholars have obtained different influence laws [10,[25][26][27][28]. As for the third shear, the equation based on the force method was derived [29]. However, the calculation accuracy is not high due to the oversimplified model. Pan, L [30], on the other hand, established the equation for the third shear by introducing the rotation constraint coefficient and shear distribution coefficient. The proposal that the side column design should take into account the third shear was put forward.
Based on existing research results, it is clear that there are few methods for calculating the additional internal forces, and the accuracy of the calculation needs to be improved. Furthermore, the influence law of the stretching plan on the third bending moment is still highly controversial. Therefore, a calculation method for additional internal forces in single-story multi-span prestressed concrete frame columns is proposed in this paper. The equations for the third bending moment and third shear based on equivalent lateral stiffness were established. Applying engineering examples, the calculation methods were verified by FEA and previous methods. Finally, finite element refinement modeling was applied to analyze the influence law of the stretching plan on the interlayer distribution of the third moment. This study provides fresh ideas to broaden theoretical research on the lateral restraint effect, and the research results can provide some reference for engineering design.

Basic Assumptions
The following basic assumptions were made based on the references [25,31]. (1) The column section was still flat after being affected by shear deformation, but it no longer remained perpendicular to the axis after deformation. That is, it conformed to the assumption of Timoshenko beams. (2) The column-end rotation angle was equal to the beam-end rotation angle, noted as θ.
(3) Define the effective prestress of prestressed tendons as σ pe , the tension control stress as σ con and the total prestress loss as σ l . The equivalent load effect at the beam end is N pe = (σ con − σ l )A p , while excluding the action of the radial equivalent load and the concentrated bending moment. (4) The structural forces were in the linear elastic phase, ignoring the plastic deformation caused by shrinkage and creep of the concrete. (5) Assuming equal column heights and equal beam sections, the displacement of the column top was linearly related to the distance from the immobility point.

Equivalent Lateral Stiffness
Generally, the deformation of frame columns under horizontal load includes bending deformation and shear deformation. In frame structures, because of the large shear span ratio of columns, the cross-sectional shear deformation is small and generally negligible [32]. However, in large-span prestressed concrete frame structures, the floor heights are generally not high while column section sizes are large, thus the shear span is relatively small. Shear deformation accounts for a large proportion of the total deformation. Therefore, the lateral stiffness of the column will exhibit a large error if the effect of shear deformation is not considered. In the traditional structural mechanics displacement method, the slopedeflection equation for the bar element only takes account of bending deformation but not shear deformation [33]. Therefore, the slope-deflection equation for the bar element was proposed, where the effect of shear deformations and bending deformation were taken into account in this paper. On this basis, the concept of equivalent lateral stiffness was proposed and its calculation equation was derived.
Plane bar elements are adopted. Figure 1 shows the structural calculation sketch on the side of the deformation immobility point, and bar element EK is shown in Figure 2. For columns without axial pressure action at both ends, the rotation angle of column top is θ, displacement of column top is δ n and the chord rotation is ϕ = δ n H . EI and GA are the flexural stiffness and shear stiffness of the section. The linear stiffness of the beam is i b = EI b l , and the linear stiffness of the column is i c = EI c H . In order to simplify the calculation, θ is only the angle caused by bending deformation, excluding the angle caused by shear deformation. δ is the total displacement caused by bending deformation and shear deformation of the bar element.    Based on reference [25], the curve function yb for bending shear deformation of the bar element were established. The equ bending moment in the cross-section is: The equilibrium condition of the shearing force is:

− =
The equilibrium Equations (1) and (2) can be solved as follo where μ is the distribution in homogeneity coefficient of shear st and μ = 1.2 is taken for the rectangular cross-section. A, B and determined.
The solution of the above equation is: Equations (5) and (6) are the curve functions of bending de formation of the bar element.
According to the boundary conditions of end K, when x = H Based on reference [25], the curve function y b for bending deformation and y s for shear deformation of the bar element were established. The equilibrium condition of the bending moment in the cross-section is: The equilibrium condition of the shearing force is: The equilibrium Equations (1) and (2) can be solved as follows: where µ is the distribution in homogeneity coefficient of shear stress in the cross-section, and µ = 1.2 is taken for the rectangular cross-section. A, B and C are coefficients to be determined. According to the boundary conditions of end E, x = 0, y b (0) = 0, y s (0) = 0, y' b = θ, then A = θ, B = C = 0.
The solution of the above equation is: Equations (5) and (6) are the curve functions of bending deformation and shear deformation of the bar element.
According to the boundary conditions of end K, when x = H, y' b (H) = 0, y b (H) + y s (H) = δ n , M EK and M KE are expressed as follows: So α n can be expressed as follows: where β is the ratio of the bending moment borne by the column to the sum of the bending moment of the beam on both sides. Since the direction of the bending moment of the beam and column are opposite, β is negative. To simplify the calculation, according to engineering experience, β = − 1 3 is obtained by reference [34]. Let k = i bl +i br i c , for the same span, the correction factor of the linear stiffness ratio of beam-column is λ.
Then α n can be simplified as follows: Equation (18) is the correction factor for the linear stiffness ratio of beams and columns to correct the lateral stiffness resistance, where λ = 0.78 when the number of spans is odd, λ = 0.83 when the number of spans is even, and λ = 1 for other columns α j [34].
It should be noted that for single-story and single-span prestressed frame structures, only the rotational restraint of the beam on one side of the column needs to be considered, then i bl = 0, k = i b i c . Equation (17) can be simplified as follows: Once the correction factor of lateral stiffness is found, the equivalent lateral stiffness of the frame column can be solved by using Equation (13).

Third Shear
According to the deformation compatibility principle in structural mechanics, the column top displacement of the side column is the sum of the axial compression displacements of all span beams. δ 1 is expressed as follows: (20) According to the assumption (5), the sum of the axial compression displacement of all spans caused by stretching of prestressed steel tendons is δ bj ,as follows: The position of the immovable point is determined according to reference [4]. As shown in Figure 1, the distance from the immovable point to the side column is y 1 and let the distance from the n-th column to the immovable point be y n , then y n = y 1 − ∑ n−1 j=1 L j , as follows: The preload of the beam is reduced for each span, which is represented by N bj (j = 1, 2 . . . n), as shown in Figure 3, as follows: It should be noted that for single-story and single-span prestressed only the rotational restraint of the beam on one side of the column needs then ibl = 0, = . Equation (17) can be simplified as follows: = 0.5 + 6 + 0.67 1 + 12 2 + 6 + 0.67 1 + 12 Once the correction factor of lateral stiffness is found, the equivalen of the frame column can be solved by using Equation (13).

Third Shear
According to the deformation compatibility principle in structura column top displacement of the side column is the sum of the axial comp ments of all span beams. δ1 is expressed as follows:

= ∆
According to the assumption (5), the sum of the axial compression all spans caused by stretching of prestressed steel tendons is δbj,as follow = The position of the immovable point is determined according to shown in Figure 1, the distance from the immovable point to the side col the distance from the n-th column to the immovable point be yn, then as follows: The preload of the beam is reduced for each span, which is represen 2…n), as shown in Figure 3, as follows: Subsequently, )y n EA Equation (24) can be simplified as: Subsequently, Equation (24) can be simplified as: Based on assumption (5), it is known that the displacement of the column top is linearly related to the distance from the immobile point. As a result of , V j can be calculated by using: V j is expressed as follows: The shear of column j can be calculated by using: Equation (28) is the formula for calculating the third shear of column j for single-story and multi-span prestressed frame structures.

Third Moment
The third moment of the column bottom section is: where η j is the height ratio of the point of changing curve, which can be obtained by reference [34]. Substituting V j into Equation (29), the third moment of the bottom section of the column j can be calculated by using: Equation (30) is the third bending moment calculation formula for the bottom section of column j for single-story and multi-span prestressed frame structures. Replace η j with (1 − η j ) to get the third bending moment of the top section of column j, noted as M t j . From Equations (28) and (30), it can be seen that the third shear and third bending moment are related to the equivalent lateral stiffness of the column, the precompression of the beam and the number of spans. The side column suffers the largest additional internal force.

Establishment of Finite Element Model
A project example was modeled and analyzed to verify the accuracy of the calculation method in this paper with ABAQUS. As shown in Figure 4a,b, a single-story and singlespan prestressed concrete frame structure with a span of 20 m and a floor height of 5 m is used for the project. Both beams and columns are made of C40 concrete. The beam cross-sectional size is b × h = 400 mm × 1300 mm, and column cross-sectional size is b × h = 600 mm × 600 mm. The prestressed beam is equipped with unbonded prestressing tendons 7φ s 15.2. The prestressing tendons are all made of high-quality steel strands with low relaxation of 1860 grade. The tensioning control stress is 0.75 f ptk , prestress loss is 20%, and the effective prestress is N pe = 1416 kN. Prestressing tendons are arranged in a double parabolic shape, tensioned at both ends with a sagittal height of 150 mm at the midspan and 150 mm at the beam-end. A total of 10 models are analyzed by varying the column section height, the effective preload value at the beam end and increasing the number of spans. Parameters of finite element models are shown in Table 1 Note: Assuming that the equivalent lateral stiffness of the frame column of YKJ1 is 1, the lateral stiffness of the frame column of the other specimens is expressed as the ratio to it.  In the modeling process, the concrete was simulated using an eight-node linear reduced-integral 3D solid unit C3D8R, and all steel bars were simulated by truss unit T3D2 [35,36]. An elastic perfectly plastic model was used for ordinary steel bars and prestressed steel strands. The concrete damage plasticity (CDP) model for concrete in ABAQUS was  Note: Assuming that the equivalent lateral stiffness of the frame column of YKJ1 is 1, the lateral stiffness of the frame column of the other specimens is expressed as the ratio to it. In the modeling process, the concrete was simulated using an eight-node linear reduced-integral 3D solid unit C3D8R, and all steel bars were simulated by truss unit T3D2 [35,36]. An elastic perfectly plastic model was used for ordinary steel bars and prestressed steel strands. The concrete damage plasticity (CDP) model for concrete in ABAQUS was used [37][38][39][40]. The constitutive was adopted according to the Chinese Standard for Design of Concrete Structures (GB 50010-2010) [8]. In defining the interaction between the steel bar and concrete, the separate modelling was used to couple the reinforcement to the concrete with degrees of freedom in the 'Embed' method, and the connection was made by means of 'Tie' constraints on the frame beams and columns. [41][42][43]. A fixed boundary condition (U1 = U2 = U3 = UR1 = UR2 = UR3 = 0) was applied to the bottom surface of the column as an initial condition. The cooling method was used to apply prestress to the truss unit T3D2 [44,45]. To improve analysis efficiency, for the multi-span prestressed frame structures YKJ8, YKJ9 and YKJ10, axisymmetric models were created by imposing symmetric boundary conditions on the split boundary of the symmetry axis [46,47]. The finite element model of YKJ1 is shown in Figure 5. Details of the material parameter values for the finite element model are shown in Table 2.
means of 'Tie' constraints on the frame beams and columns. [41][42][43]. A fixed bound condition (U1 = U2 = U3 = UR1 = UR2 = UR3 = 0) was applied to the bottom surface of column as an initial condition. The cooling method was used to apply prestress to truss unit T3D2 [44,45]. To improve analysis efficiency, for the multi-span prestres frame structures YKJ8, YKJ9 and YKJ10, axisymmetric models were created by impos symmetric boundary conditions on the split boundary of the symmetry axis [46,47]. finite element model of YKJ1 is shown in Figure 5. Details of the material parameter va for the finite element model are shown in Table 2.

Results and Discussion
The third shear and the third bending moment of side columns obtained by FEA the method of this paper are listed in Table 3. The calculation results of the method of paper are in good agreement with the finite element analysis values. The results show the relative error of the third shear is 0.36% to 7.18%, and the relative error of the th bending moment is 0.66% to 12.52%. The calculation accuracy meets engineering requ ments, which verifies the reliability of the calculation method of this paper. Except YKJ4, the calculation accuracy of the third shear obtained by the method of this pape generally higher than that of the third bending moment. The reason for this is that standard height ratio of the point of the changing curve is restricted to the situatio uniform load [34], which differs from the case of horizontal load to the column top in study. The calculated value of the method in this paper is consistently larger than the fi

. Results and Discussion
The third shear and the third bending moment of side columns obtained by FEA and the method of this paper are listed in Table 3. The calculation results of the method of this paper are in good agreement with the finite element analysis values. The results show that the relative error of the third shear is 0.36% to 7.18%, and the relative error of the third bending moment is 0.66% to 12.52%. The calculation accuracy meets engineering requirements, which verifies the reliability of the calculation method of this paper. Except for YKJ4, the calculation accuracy of the third shear obtained by the method of this paper is generally higher than that of the third bending moment. The reason for this is that the standard height ratio of the point of the changing curve is restricted to the situation of uniform load [34], which differs from the case of horizontal load to the column top in this study. The calculated value of the method in this paper is consistently larger than the finite element analysis value, which is because of the simplified rod model used in the derivation of the formula in this paper. The influence of plastic deformation caused by compression and creep of concrete is ignored.  Figure 6 shows the effect of equivalent lateral stiffness on the third shear and third bending moment. The third shear and the third bending moment rise as the equivalent lateral stiffness of the frame column increases. This is because the lateral restraint capacity of the frame column on the prestressed beam gradually enhances with the increase in lateral stiffness of the frame column. The distribution ratio of the axial prestressing force involved in the beam becomes increased, and thus the additional internal force increases. As a result, it is advantageous to improve the structure's flexibility in order to lessen the effect of additional internal forces while fulfilling the bearing capacity and normal use requirements during design.  Figure 6 shows the effect of equivalent lateral stiffness on the third sh bending moment. The third shear and the third bending moment rise as lateral stiffness of the frame column increases. This is because the lateral res of the frame column on the prestressed beam gradually enhances with the eral stiffness of the frame column. The distribution ratio of the axial pres involved in the beam becomes increased, and thus the additional internal f As a result, it is advantageous to improve the structure's flexibility in orde effect of additional internal forces while fulfilling the bearing capacity an requirements during design. According to Equation (20), the capital level displacement of the side sum of the axial compression deformation of each span caused by the stre stressed steel tendons. As shown in Figure 7, the displacement values are 0. mm, 2.281 mm and 2.699 mm for single span, double-span, three-span a frames, respectively. It can be seen that, as the number of spans increases, th According to Equation (20), the capital level displacement of the side column is the sum of the axial compression deformation of each span caused by the stretching of prestressed steel tendons. As shown in Figure 7, the displacement values are 0.971 mm, 1.746 mm, 2.281 mm and 2.699 mm for single span, double-span, three-span and four span frames, respectively. It can be seen that, as the number of spans increases, the capital level displacement increases. As a result, the axial preload force allocated to the column by the beam increases due to the third shear force, which in turn causes an increase in the third bending moment, as shown in Figure 8. In actual engineering, for continuous multispan prestressed structures, it is advisable to set post-cast zones and carry out segmental tensioning of prestressing tendons or take other measures during construction to reduce the influence of span number, especially for side columns. In addition, the values of the third shear and the third bending moment of the side columns increase with the increase in the effective prestressing force value and conform to a linear relationship, which is consistent with Equations (28) and (30).
Buildings 2022, 12, x FOR PEER REVIEW 11 of 18 displacement increases. As a result, the axial preload force allocated to the column by the beam increases due to the third shear force, which in turn causes an increase in the third bending moment, as shown in Figure 8. In actual engineering, for continuous multi-span prestressed structures, it is advisable to set post-cast zones and carry out segmental tensioning of prestressing tendons or take other measures during construction to reduce the influence of span number, especially for side columns. In addition, the values of the third shear and the third bending moment of the side columns increase with the increase in the effective prestressing force value and conform to a linear relationship, which is consistent with Equations (28) and (30).   Figure 9 shows the third shear relative errors of the YKJ1-YKJ10 side columns as determined using the method of FEA, references [25,30] and our study. It should be noted that reference [25] only gives the calculation method for single-story and single-span beam increases due to the third shear force, which in turn causes an increa bending moment, as shown in Figure 8. In actual engineering, for continuo prestressed structures, it is advisable to set post-cast zones and carry out s sioning of prestressing tendons or take other measures during construction influence of span number, especially for side columns. In addition, the valu shear and the third bending moment of the side columns increase with the effective prestressing force value and conform to a linear relationship, whic with Equations (28) and (30).   Figure 9 shows the third shear relative errors of the YKJ1-YKJ10 side c  Figure 9 shows the third shear relative errors of the YKJ1-YKJ10 side columns as determined using the method of FEA, references [25,30] and our study. It should be noted that reference [25] only gives the calculation method for single-story and single-span structures, so only the third shears of YKJ1 to YKJ7 are shown in Figure 9. Table 4 shows the comparison of the average relative error of the third shear. method of [25] are all larger. The error range is 10.44% to 83.88%, with an av 31.06% and a dispersion coefficient of 0.86. It implies that the third shear i exaggerated, with a substantial error. Especially for the YKJ4 specimen w large lateral stiffness resistance, the error is the largest. The relative error shear calculated using the method of [30] varied from 0.95% to 37.56%, wi error of 11.52% and a dispersion coefficient of 1.06. This indicates that the c curacy of this method is relatively high, but the dispersion of the calculat different specimens is large. The relative error of the third shear obtained b tion method in this paper is 0.36% to 7.18%, with an average error of 3.37% sion coefficient of 0.62. It shows that the method in this paper has a small e persion, and the highest calculation accuracy.   Figure 10 shows the relative third bending moment errors of the YKJ1umns as determined using the method of FEA, reference [25] and our study relative errors calculated by reference [25] and our methods are shown in be seen that the accuracy of calculating by this paper's method is significant that of reference [25]. The third bending moment in reference [25] is 2.45% to than FEA values, with the largest error for YKJ5. The average relative erro the dispersion coefficient is 0.81. The relative error by this paper's method 0.66% to 9.23%, with an average value of 6.41% and a dispersion coefficien clear that both the calculation error and dispersion of the method in this pa which fully demonstrates that the equation in this paper can accurately eva bending moment.  The calculation accuracy of the third shear achieved by the method proposed in this paper is much greater than that of previous methods. The third shears obtained using the method of [25] are all larger. The error range is 10.44% to 83.88%, with an average error of 31.06% and a dispersion coefficient of 0.86. It implies that the third shear is significantly exaggerated, with a substantial error. Especially for the YKJ4 specimen with relatively large lateral stiffness resistance, the error is the largest. The relative errors of the third shear calculated using the method of [30] varied from 0.95% to 37.56%, with an average error of 11.52% and a dispersion coefficient of 1.06. This indicates that the calculation accuracy of this method is relatively high, but the dispersion of the calculation results for different specimens is large. The relative error of the third shear obtained by the calculation method in this paper is 0.36% to 7.18%, with an average error of 3.37% and a dispersion coefficient of 0.62. It shows that the method in this paper has a small error, little dispersion, and the highest calculation accuracy. Figure 10 shows the relative third bending moment errors of the YKJ1-YKJ7 side columns as determined using the method of FEA, reference [25] and our study. The average relative errors calculated by reference [25] and our methods are shown in Table 5. It can be seen that the accuracy of calculating by this paper's method is significantly higher than that of reference [25]. The third bending moment in reference [25] is 2.45% to 47.62% larger than FEA values, with the largest error for YKJ5. The average relative error is 20.3% and the dispersion coefficient is 0.81. The relative error by this paper's method ranges from 0.66% to 9.23%, with an average value of 6.41% and a dispersion coefficient of 0.46. It is clear that both the calculation error and dispersion of the method in this paper are small, which fully demonstrates that the equation in this paper can accurately evaluate the third bending moment.

Accuracy Analysis
Reference [25] assumed consolidation at the column end and did no effect of the rotation angle at the column end. At the same time, the solu stiffness resistance only considered the influence of bending deformation o which led to its large calculated value. Therefore, the lateral restraint capa umn was seriously overestimated, resulting in large errors in the calculated ditional internal forces. The column end restraint and shear distribution we by adopting the rotational constraint coefficient and the shear distribution reference [30]. However, the lateral stiffness did not take into account the deformation of the column section. Poor calculation accuracy for structures eral stiffness indicated that it is not universally applicable. In contrast, the e eral stiffness proposed in this paper could accurately assess the lateral restr of columns, by considering both shear and bending deformation effects. could reasonably determine the distribution ratio of the frame column in axial preload of the beam and ensure the accuracy of calculations of the addi  Table 5. Average relative errors and dispersion coefficient of the third moment.

Statistical Quantities Method of [25] Method of This Paper
Average relative error 20.30% 6.41% Dispersion coefficient 0.81 0.46

Accuracy Analysis
Reference [25] assumed consolidation at the column end and did not consider the effect of the rotation angle at the column end. At the same time, the solution of lateral stiffness resistance only considered the influence of bending deformation of the column, which led to its large calculated value. Therefore, the lateral restraint capacity of the column was seriously overestimated, resulting in large errors in the calculated values of additional internal forces. The column end restraint and shear distribution were considered by adopting the rotational constraint coefficient and the shear distribution coefficient in reference [30]. However, the lateral stiffness did not take into account the effect of shear deformation of the column section. Poor calculation accuracy for structures with high lateral stiffness indicated that it is not universally applicable. In contrast, the equivalent lateral stiffness proposed in this paper could accurately assess the lateral restraint resistance of columns, by considering both shear and bending deformation effects. This stiffness could reasonably determine the distribution ratio of the frame column involved in the axial preload of the beam and ensure the accuracy of calculations of the additional internal force value. Therefore, compared with the existing calculation methods, the method in this paper has higher accuracy. Equivalent lateral stiffness is a key factor in determining the accuracy of additional internal force calculations.

Influence of Tensioning Order during the Construction Phase
For multi-story prestressed frames, the current design method employs the "integral casting, one-time tensioning" structural model in overall structural analysis [48]. It ignores the effect of the stretching plan of prestressed tendons on structural internal forces, resulting in the difference between the design value and the actual value, particularly for the bending moment [49]. However, it is very difficult to establish a theoretical formulation considering the dynamic changes in the construction process. Therefore, in this paper, finite element simulations of the construction process with different stretching plans for multi-story prestressed frames were carried out. We hope to find out the influence law of stretching plans on the interlayer distribution of the third bending moment.
A finite element model of a single-span three-story prestressed concrete frame was built using the fundamental characteristics of the YKJ2 specimen. The impacts of the several-layer and layer-by-layer stretching plan on the third bending moment were investigated. Figures 11 and 12 show X-direction displacement cloud images of the several-layer stretching plan and layer-by-layer stretching plan. The distribution of the third bending moment of the side columns between layers for the several-layer method and the layer-bylayer method are shown in Figures 13 and 14. It should be noted that the third moment value of the column after tensioning of the third layer of prestressed tendons in Figure 13 is also the final value calculated by the current design method. A finite element model of a single-span three-story prestressed concrete frame was built using the fundamental characteristics of the YKJ2 specimen. The impacts of the several-layer and layer-by-layer stretching plan on the third bending moment were investigated. Figures 11 and 12 show X-direction displacement cloud images of the several-layer stretching plan and layer-by-layer stretching plan. The distribution of the third bending moment of the side columns between layers for the several-layer method and the layerby-layer method are shown in Figures 13 and 14. It should be noted that the third moment value of the column after tensioning of the third layer of prestressed tendons in Figure 13 is also the final value calculated by the current design method.  A finite element model of a single-span three-story prestressed concrete frame was built using the fundamental characteristics of the YKJ2 specimen. The impacts of the several-layer and layer-by-layer stretching plan on the third bending moment were investigated. Figures 11 and 12 show X-direction displacement cloud images of the several-layer stretching plan and layer-by-layer stretching plan. The distribution of the third bending moment of the side columns between layers for the several-layer method and the layerby-layer method are shown in Figures 13 and 14. It should be noted that the third moment value of the column after tensioning of the third layer of prestressed tendons in Figure 13 is also the final value calculated by the current design method.  The tensioning scheme has a path effect on the interlaminar distributi bending moment between layers. Compared to the total bending moment u loads in reference [25], the distribution of the final value of the third bendi the layer-by-layer method is staggered with it, with odd layers mutually aba and even layers superimposing unfavorably. The layer-by-layer method is o direction of the external load moment, and each layer of the column mome tually reduced. For the same structure, different construction schemes corr ferent structural response histories and final responses of the third bending Tensioning scheme has a time-varying effect on the distribution of the The tensioning scheme has a path effect on the interlaminar distributio bending moment between layers. Compared to the total bending moment u loads in reference [25], the distribution of the final value of the third bendin the layer-by-layer method is staggered with it, with odd layers mutually abat and even layers superimposing unfavorably. The layer-by-layer method is o direction of the external load moment, and each layer of the column mome tually reduced. For the same structure, different construction schemes corre ferent structural response histories and final responses of the third bending Tensioning scheme has a time-varying effect on the distribution of the moment between layers [50]. The bottom and top of the column at the botto The tensioning scheme has a path effect on the interlaminar distribution of the third bending moment between layers. Compared to the total bending moment under external loads in reference [25], the distribution of the final value of the third bending moment in the layer-by-layer method is staggered with it, with odd layers mutually abating favorably and even layers superimposing unfavorably. The layer-by-layer method is opposite to the direction of the external load moment, and each layer of the column moment can be mutually reduced. For the same structure, different construction schemes correspond to different structural response histories and final responses of the third bending moment.
Tensioning scheme has a time-varying effect on the distribution of the third bending moment between layers [50]. The bottom and top of the column at the bottom layer of the construction of the several-layer method are the control section. The third bending moment at the top of the bottom column of the first layer beam after prestressing tendon tensioning is the largest, which is 1.5 times the final value and taken as the design control value. The control sections of layer-by-layer method construction are the bottom and top layer column bottom, the middle layer column bottom, and the column top. After stretching the prestressing tendons of the third layer beam, the third bending moment at the bottom section of the top layer column is the largest, which is 3.8 times the final value calculated by the current design method and is the design control value. Therefore, the tensioning construction scheme is different, the time and size of the design control value is different, and the control section position of the third bending moment is different.
It is advisable to use the third bending moment design control value of each story column for crack width and bearing capacity verification during the construction stage. Generally, prestressed concrete frame structures have few layers. The vertical loads often play a controlling role in the service phase, and the design is usually on the safe side without considering the effect of the third bending moment. However, during the construction stage, the prestressed beam formwork and supports are not removed, and the vertical load effect is very small. At this time, if the third bending moment value is large, it may lead to shear and bending failure at the column end. Therefore, it is recommended to check the calculation for the construction stage.

Conclusions
To determine the additional internal forces of columns caused by prestressed tendons more accurately, a calculation method based on equivalent lateral stiffness was proposed in the paper. Shear and bending deformation were considered in the equivalent lateral stiffness of the column. The deformation compatibility and force equilibrium conditions were established according to the equivalent lateral stiffness, and equations for the third shear and third bending moment were proposed based on these conditions. The proposed method was verified by FEA and compared with previous methods to demonstrate its superior performance. In addition, the effect of the stretching plan on the interlayer distribution of the third bending moment was analyzed by a single-span three-story finite element model. The following conclusions can be drawn:

1.
The slope-deflection equation for the bar element was derived considering the effect of shear deformations and bending deformation. On this basis, the equation for the equivalent lateral stiffness of frame columns was derived. The equation could effectively evaluate the lateral restraint capacity of the column and had good applicability.

2.
The third shear and third bending moment values calculated by the proposed method were in good agreement with FEA results. The average relative error of the third shear is 3.37%, and that of the third moment is 6.41%. Compared to the existing methods, the calculation method in this paper had higher calculation accuracy.

3.
The proposed method can be directly applied to solve the additional internal forces in columns for post-tensioned single-story multi-span prestressed concrete frames with prestressed tendons stretched at beam both ends. The accuracy of the calculation meets engineering requirements. Considering the effect of equivalent lateral stiffness, span number and preload value, the method has more reliability and applicability. 4.
The stretching plan has a significant path effect and time-varying effect on the interlayer distribution of the third moment. It is advisable to use the third bending moment design control value of each story column for crack width and bearing capacity verification during the construction stage.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.

Conflicts of Interest:
The authors declare no conflict of interest.