# Dynamic Behavior of Externally Prestressed Continuous Beam considering Second-Order Effect

Authors: De-Ping Fang
Journal: Mathematical Problems in Engineering
 (2020)

---

## Abstract

Considering precisely the second-order deformation of external tendon, the analytical solution of natural frequencies of 2-span externally prestressed continuous beam was obtained by the energy method. The effect of external prestress compression softening is between the zero effect of unbonded prestress compression and the effect of axial outside compression and is determined by the influence coefficient ranging within 0∼1. The influence coefficient is mainly related to the number of deviators and slightly related to tendon layout. Without deviator, the influence coefficient is 1, and the effect of external prestress compression softening is the same as the effect of axial outside compression. As the number of deviators increases, the influence coefficient gradually decreases from 1 to near 0, and the effect of external prestress compression softening is close to zero effect of unbonded prestress compression. With one or more deviators, the effect of external prestress compression softening is negligible. As the eccentricity and area of tendon increase, only the first symmetric frequency increases obviously, and other frequencies almost remain unchanged. The influence of tendon layout linear transformation on the frequency is negligible.


---

## Body

## 1. Introduction

Externally prestressed concrete structure has been widely used in the strengthening of existing concrete bridges as well as in the construction of various new structures. Lou et al. [1, 2] and Anwar et al. [3] studied the static nonlinear behavior of beam externally prestressed with steel tendon and FRP cables, and Fang [4] investigated the friction effect of externally prestressed beam. However, the dynamic behavior of externally prestressed beam was not investigated in detail. The analysis of the natural frequencies of externally prestressed simply supported or continuous beams is of vital importance in dynamic behavior. The natural frequency of a simply supported beam under axial outside compression is(1)ωi=iπl1mEIiπl2−N,i=1,2,3,…,where EI is flexural stiffness of beam, m is the mass per unit length of beam, l is beam span, i is the mode number, and N is the axial outside compression. Equation (1) reveals that N leads to compression softening (N softening effect hereinafter) and then decreases the natural frequencies. The prestress tendon applies axial inside compression to the beam. There exists a problem whether the prestress compression softening effect is the same as N softening effect. Different researchers presented different conclusions.Saiidi et al. [5] determined the natural frequencies of a prestressed concrete bridge using equation (1). Dallasta and Dezi [6] pointed out Saiidi et al.’s [5] approach to consider the prestress compression as N is incorrect, and indicated that the effect of prestress compression on the beam natural frequencies is negligible based on a linear model. Deak [7] also pointed out that prestress compression does not reduce the natural frequencies. However, the view was not supported by any analytical or mathematical proof. Jain and Goel [8] further pointed out that because the tendon becomes an integral part of the system, the prestress compression cannot be treated as N and does not affect the natural frequencies of beams. Jaiswal [9] investigated the first natural frequency of beams by finite element method (FEM). For the beams with bonded tendon, the prestress compression does not have any appreciable effect on the first natural frequency. For the beams with unbonded tendon, the first natural frequency significantly changes with the prestress compression and eccentricity of the tendon. Jaiswal’s [9] conclusions are only based on FEM results and are not supported by any analytical or mathematical proof. Kanaka [10] treated prestress compression as N and concluded that the prestress compression reduces the natural frequency of the lower modes based on a Rayleigh–Ritz formulation. Chan [11] also treated prestress compression as N and indicated that the natural frequencies of a prestressed bridge decrease as the prestress compression increases. Dallasta and Leoni [12] presented a general formulation for the vibration of concrete beams, prestressed by internal frictionless tendons using kinematic relations of small displacements. The formulation for the beams does not include the effect of the compression; however, the natural frequencies decrease as the prestress compression increases. Kerr [13] studied experimentally and analytically the dynamic response of a prestressed beam and pointed out that the magnitude of the prestress compression for a tendon that passes through the centroid of the beam cross section has zero effect on the dynamic response of the beam. The analytical model was based on a linear formulation to study the dynamic response of the prestressed beam. It was limited to straight tendons that pass through the centroid of the beam. Simsek and Kocaturk [14] indicated that the deflections of the beams increase as the prestress compression increases. The prestress compression softening decreases the beam stiffness and natural frequencies. Wang and Ren [15] considered the additional potential energy of the prestress compression and concluded that the prestress compression reduces the low transverse natural frequencies of the bridge. Jiang and Ye [16] transferred the eccentric unbonded prestress compression to N and a couple and concluded the prestress compression results in compression softening effect. Hamed and Frostig [17] derived the equations of motion for a prestressed beam using the variational principle of virtual work following Hamilton’s principle. The mathematical model is rigorous and general and is valid for any kind of boundary and continuity conditions as well as any tendon layout. It was mathematically proven that the prestress compression does not affect the natural frequencies of bonded or unbonded prestressed beams as opposed to some research conclusions.Hamed and Frostig [17] pointed out zero softening effect of prestress compression of bonded or unbonded tendon. The external tendon applies external prestress compression EP to the beam; does EP lead to softening effect? Miyamoto et al. [18] presented the dynamic behavior of prestressed simply supported beam strengthened with external tendon. EP effect is simply treated as N effect. Shi et al. [19] adopted similar method of Miyamoto et al. [18] to treat EP effect as N effect and obtained the analytical frequency expression of externally prestressed concrete continuous beam. The analytical values were compared with the measured frequencies. Simsek and Kocaturk [20] treated the eccentric EP as N and a couple, and the EP effect is equivalent to N effect. Ji et al. [21] proposed sophisticated formulas to predict the vertical bending vibration frequencies of prestressed concrete box girders with corrugated steel webs based on Hamilton’s energy variational principle. The effects of external tendon and EP were not taken into consideration.There exists the eccentricity loss in external tendon, but there is zero eccentricity loss in unbonded tendon. The loss of external tendon eccentricity decreases the flexural stiffness and consequentially decreases the natural frequencies. The loss of eccentricity depends on the layout of external tendon. Therefore, the influence of EP on the natural frequencies is also related to the layout of external tendon. It can be illustrated simply in Figure1. During the vibration of a simply supported beam, the internal unbonded tendon ACB along the beam axis bends, but its length and strain energy remain constant. The prestress compression does not influence the dynamic behavior. It conforms with Hamed and Frostig’s [17] conclusion of zero effect of prestress compression on the natural frequencies of bonded or unbonded prestressed beams. The external tendon AB anchored at two ends remains straight; therefore it shortens. The loss of tendon AB strain energy transfers to the beam, increases the strain energy and deformation of beam, and decreases the beam stiffness and natural frequency. If N is equal to EP, the difference between the loss of AB tendon strain energy and the work of N is high order negligible; EP softening effect is the same as N softening effect and influences the natural frequencies.Figure 1
Deformation of external and internal unbonded tendons.Based on energy method, adopting the vibration modes of usual simply supported beam, considering precisely the second-order effect, Fang [22] analyzed the natural frequencies of externally prestressed simply supported beam and EP softening effect. It was concluded that EP softening effect depends on the loss of external tendon eccentricity. As the number of deviators increases from zero to a large number, the layout of the external tendon is close to unbonded tendon. EP softening effect decreases from N effect to almost zero. It is consistent with the conclusion of Hamed and Frostig [17]. Fang [22] only discussed the simply supported beam. The continuous beams were widely adopted in engineering, and the analysis of continuous beam dynamic behavior is also indispensable. Using the method proposed by Fang [22], adopting the vibration modes of usual continuous beam, the natural frequencies of continuous beam were also analyzed, and similar conclusions are presented in this paper.


## 2. Natural Frequency Equations of Externally Prestressed Beams

In the author’s view, there exist some problems in Miyamoto et al.’s [18] and Shi et al.’s [19] natural frequency formulas. The derivation procedure is rewritten briefly in this section.For free vibration, the governing partial differential equation is(2)EI∂4y∂x4+∂2Ptcy∂x2−∂2ΔMp∂x2+m∂2y∂t2=0,where y is vibration displacement; Pt and ΔPt are prestress compression and increment; Pt0 is initial prestress compression; Ptc, ΔPtc, and Ptc0 are the horizontal components of Pt, ΔPt, and Pt0; ΔMp is the bending moment due to ΔPt; e is the eccentricity; and a is the position of deviators as shown in Figure 2; the friction between the tendon and deviator is neglected. The effect of the prestress compression and the effect of external tendon are represented by the terms Ptcy and ΔMP, respectively, in equation (2)Figure 2
Externally prestressed simply supported beam.To solve equation (2), two relationships must be established, the relationship between Ptcy and y and the relationship between ΔMP and y. Miyamoto [18] neglected the term ΔPtcy in PtcyPtcy≈Ptc0y; since ymax≪e, therefore ΔPtcy≪ΔPtce (ΔPtce is included in ΔMP as shown in Figure 3). In the author’s view, ΔPtcy is negligibly small, not due to ymax≪e. But it is due to the fact that ΔPtc is proportional to y, and y can be infinitesimal; ΔPtcy becomes a high-order infinitesimal; hence, ΔPtcy can be neglected in comparison with infinitesimal Ptc0y. Acknowledge that, in some beams, where e = 0, ymax≪e is not a valid assumption. Therefore, Ptcy can be replaced by Ptc0y; then, the relationship between Ptcy and y is established. Miyamoto et al. [18] simplified ΔMP as a uniform bending moment ΔMPE as shown in Figure 3 (equivalent-area principle hereinafter):(3)ΔMP≈ΔMPE=ΔPtecosθ+asinθ.Figure 3
Bending moment diagramsΔMP,ΔMPE of simply supported beam.Miyamoto [18] assumed that ΔPt is proportional to midspan displacement y0.5l, called proportional assumption hereinafter:(4)ΔPt=ηy0.5l,(5)ΔMP≈ΔMPE=ηecosθ+asinθy0.5l,where η=24EIecosθ+asinθ/l2μ+4λcosθμ=4cosθ−3ecosθ+asinθ2; λ = I/A+EIlt/EtAtl; A is area of the beam cross section; and Et, At, and lt are the elastic modulus, cross section area, and length of tendon. Substituting equation (5) into equation (2) gives(6)EI∂4y∂x4+Ptc0∂2y∂x2−ηecosθ+asinθ∂2y0.5l∂x2+m∂2y∂t2=0,and y0.5l in equations (4)–(6) is replaced by yx in Miyamoto’s [18] equations. Due to this replacement, equation (6) can be solved for natural frequencies:(7)ωi=iπl1mEIiπl2−Ptc0+ηecosθ+asinθ,i=1,2,3,….Using a similar method, Shi et al. [19] put forward the natural frequency formula for continuous beam as shown in Figure 4. To solve equation (2), Ptcy is replaced by Ptc0y, ΔPt=ϕy0.5l by the proportional assumption. Based on the equivalent-area principle, the areas of the bending moment diagram ΔMP due to ΔPt are equal to the uniform bending moment ΔMPE, as shown in Figure 5ΔMP=ΔMPE=Hϕy0.5l. The differential equation similar to equation (6) is obtained:(8)EI∂4y∂x4+Ptc0∂2y∂x2−Hϕ∂2y0.5l∂x2+m∂2y∂t2=0.Figure 4
Externally prestressed continuous beam.Figure 5
Bending moment diagramsΔMP,ΔMPE of continuous beam.Replacingy0.5l by yx, solving equation (8) for natural frequencies of continuous beam, the expressions of H,ϕ and frequencies are detailed in the paper by Shi et al. [19].


## 3. Some Problems in Proportional Assumption and Equivalent-Area Principle

In the derivation of Miyamoto et al.’s [18] and Shi et al.’s [19] analytical expressions of the natural frequencies, in the author’s view, there are 4 problems in the proportional assumption and equivalent-area principle.(1)
In proportional assumption, the increment of tendon tensionΔPt is proportional to the vibration displacement at midspan of beam y0.5l; the proportional coefficient is a constant in equation (4). But for different vibration mode, the same y0.5l results in different deformation and ΔPt of tendon, the proportional assumption is not valid.(2)
In order to solve equations (6) and (8) conveniently, y0.5l is replaced by yx; i.e., ΔPt is not proportional to y0.5l but is proportional to yx. This replacement is wrong because ΔPt is constant along the whole tendon without friction, but the value yx varies along x.(3)
Also for solving convenience, the equivalent-area principle simplifiesΔMP into uniform ΔMPE, which results in negligible error in simply supported beam for small a in Figure 2, but it would result in obvious error in continuous beam. In Figure 5, both midspan bending moment above axis and bending moment below axis near the middle support contribute to stiffness of continuous beam in spite of plus-minus sign. The calculation of ΔMPE is the area of ΔMP above axis minus the area of ΔMP below axis; two bending moments originally strengthen together and now wrongly counteract mutually. The continuous beam deflects upward under uniform ΔMPE, which does not meet zero displacement at middle support, so obvious error is expected in the equivalent-area principle.(4)
Ptc0y in equations (6) and (8) implies EP effect is treated as N effect without considering the influence of deviators, which is not consistent with Fang [22] conclusion.Therefore, it is imperative to develop more accurate analytical solution for continuous beam frequencies considering the second-order effects without the above problems.


## 4. Energy Method to Analyze the Natural Frequencies of Continuous Beam

In the present study, the energy method considering the second-order effects is adopted to solve the problems that existed in Shi et al.’s [19] solution, and to find the natural frequencies of the continuous beam shown in Figure 4. The assumptions in the energy method are as follows:(1)
Theith vibration mode yi=AiXixsinωit; Ai: amplitude of ith mode; antisymmetric ith mode: Xix=sinkix,kil=iπ,i=1,2,3… symmetric ith mode: Xix=sinkix−fisinhkix,fi=sinkil/sinhkil,kil=i+0.25π,i=1,2,3,… both symmetric and antisymmetric modes are the same as modes of usual continuous beam.(2)
The beam is straight after applying prestress and before vibrating.(3)
The axial deformation of the beam due to axial compression and the mass of tendon and deviators are neglected.Wheny˙i=max or yi=0, the strain energy of the beam is(9)Ub0=12EI∫02lPt0mt2dx,where Pt0mt is the bending moment due to Pt0 and mt is the bending moment due to the unit tension of external tendon; mt=ΔMP/ΔPt in Figure 5.The kinetic energy of the beam is(10)Kb0=12∫02lmy˙i2dx=Ai2ωi2rb0m.For antisymmetricith mode, rb0=l/2; for symmetric ith, mode,rb0=l/2−1/4ki+fi2sinh2kil−2kil/4ki−fisinkilcoshkil−coskilsinhkil/ki. The strain energy of external tendon is(11)Ut0=Pt02lt2EtAt.Whenyi=max or y˙i=0, the bending moment of beam is Mb1=−EIyi″−Pt0mt, and the strain energy of the beam is(12)Ub1=12EI∫02lMb12dx=EIAi2rb1+∫02lyi″Pt0mtdx+Ub0.For antisymmetricith mode, rb1=ki4l/2; for symmetric ith mode, rb1=ki3kil/2−1/4+fi2sinh2kil−2kil/4+fisinkilcoshkil−coskilsinhkil.As for calculating the deformation of the external tendon, the deformation of line OA in Figure6 is investigated firstly. The deformations of line OA in x and y axes are δx and δy; considering the second differentiation, the deformation of line OA is(13)OA¯−OA=sinθδy+cosθδx+sinθδx22lA+cosθδy22lA−sinθcosθδxδylA,where θ and lA are angle and length of line OA. It is indispensable to consider the second differentiation in analyzing EP softening. There are n contact points between the beam and external tendon, the first and the last contact points are two anchorages at two ends of the beam, and the remaining others are at the deviators. In the jth segment of external tendon as shown in Figure 7, the coordinates of the jth and (j+1)th contact points are xj,ej and xj+1,ej+1 before vibrating. Considering the horizontal displacement of the point at beam axis due to the axis bending, the coordinates of the contact points x¯j,y¯j and x¯j+1,y¯j+1 when y˙i=0 are(14)x¯j=xj−12∫0xjyi′2dx−ejyi′xj,x¯j+1=xj+1−12∫0xj+1yi′2dx−ej+1yi′xj+1,y¯j=ej+yixj,y¯j+1=ej+1+yixj+1.Figure 6
Deformation of line OA.Figure 7
Deformation of external tendon segment.The deformations of thejth segment in x and y axes δxj and δyj are(15)δxj=x¯j+1−x¯j−xj+1−xj=bjAi2+cjAi,δyj=y¯j+1−y¯j−ej+1−ej=djAi,bj=fxj−fxj+1,cj=ejXi′xj−ej+1Xi′xj+1,dj=Xixj+1−Xixj,fx=12∫0xXi′2dx.For antisymmetricith mode, fx=0.125ki2kix+sin2kix; for symmetric ith, modefx=0.125ki2kix+sin2kix+fi2sinh2kix+2kix−4fisinhkixcoskix+coshkixsinkix. Substituting δxj and δyj into equation (13), and neglecting the terms involving high-order infinitesimal Ai3 and Ai4 (Ai can be infinitesimal), give the deformation of the jth segment of external tendon:(16)Δlj=αjAi+βjAi2αj=djsinθj+cjcosθj, βj=bjcosθj+cjsinθj2/2lj+djcosθj2/2lj−cjdjsinθjcosθj/lj; θj and lj are the length and angle of the jth segment of external tendon, respectively. The whole deformation of the external tendon is(17)Δlt=Aiψi+Ai2ζi,where ψi=∑j=1n−1αj, ζi=∑j=1n−1βjAiψi is the first-order deformation due to the ith vibration mode; i.e., Aiψi=−∫02lyi″mtdx. For antisymmetric ith mode, Aiψi=0. Ai2ζi is the second-order deformation, which involves not only the horizontal projection variation of contact point due to the axis bending, i.e., bj term, but also cj2,dj2,cjdj terms in equation (16). If considering only bj term, it results in incorrect conclusion: EP effect is the same as N effect. Precisely considering cj2,dj2,cjdj terms leads to correct conclusion: EP effect is between N effect and zero effect of unbonded prestress compression. The tension of external tendon is(18)Pt1=Pt0+EtAtAiψi+Ai2ζilt.The strain energy of external tendon isUt1=Pt12lt/2EtAt; neglecting the terms involving high-order infinitesimal Ai3 and Ai4 gives(19)Ut1=Ut0+Pt0Aiψi+Ai2ζi+EtAtAiψi22lt.Substituting each term in equations (9)–(12) and (19) into the equation of energy method, Ub0+Kb0+Ut0=Ub1+Ut1, and deducting the term ∫02lyi″Pt0mtdx in Ub1 and the term Pt0Aiψi=−Pt0∫02lyi″mtdx in Ut1 mutually give(20)Ai2ωi2rb0m=EIAi2rb1+Pt0Ai2ζi+EtAtAiψi22lt.The natural frequency isωi=1/rb0mEIrb1+Pt0ζi+EtAtψi2/2lt. For antisymmetric ith mode, ψi=0; ωi can be simplified as(21)ωi=iπl1mEIiπl2−CpiPt0,i=1,2,3….The influence coefficientCpi=−2lζi/iπ20<Cpi≤1; it reflects EP softening effect. For symmetric ith mode, there is no concise analytical frequency solution of usual continuous beam under N. Cpi is defined as follows: when Pt0=0, yi=max under N; the beam axial horizontal length shortens 2Ai2fl due to axis bending; the equation of energy method: Kb0=Ut1+Ub1−2Ai2flN, obtaining ωi=1/rb0mEIrb1−N2fl+EtAtψi2/2lt; when Pt0=N, the influence coefficient Cpi=−ζi/2flωi for the symmetric mode is(22)ωi=1rb0mEIrb1−CpiPt02fl+EtAtψi22lt,i=1,2,3….In comparison of equations (21) and (22) with equation (1) under N, EP or Pt0 does affect the natural frequencies of beam, but its effect depends on the coefficient Cpi which is related to the layout of external tendon. When Cpi=1, EP effect is the same as N effect, similar to the straight external tendon AB in Figure 1 when Cpi=0, zero EP softening effect, similar to the internal unbonded tendon ACB in Figure 1.


## 5. Numerical Examples

For the externally prestressed continuous concrete beam shown in Figure4, l = 5 m, with width 0.4 m and height 0.15 m of rectangular cross section; E = 32.5 GPa; the mass per unit length of the beam m = 0.8 t/m; external tendon Et = 200 GPa; the initial prestress σt0=1000MPa; l1 = l2 = l3 = l4 = 1.25 m; and the deviators are distributed uniformly spaced along the span. 8 beams were studied with 2 layouts of tendon (A, B), 2 two eccentricities (large eccentricity L, small eccentricity S), and two areas of the tendon (1 tendon At1 = 137 mm2, 2 tendons At2 = 274 mm2), which are listed in Table 1. For all layouts of tendon, e1 = 0 (see Figure 4); all deviators are located at parabola. A-layout can be linearly transformed into B-layout. Large eccentricity L is twice as large as small eccentricity S.Table 1
Layouts and eccentricities of external tendons.
Beam codee2 (mm)e3 (mm)e4 (mm)e5 (mm)AL406500281250AS203250141125BL34437594500BS17218847250Table2 shows the relationship between influence coefficient Cpi and the number of contact points n for straight tendon, the straight tendon is along the axis of beam, and the contact points are uniformly distributed. Table 3 shows the relationship between influence coefficient Cpi and the number of contact points n for parabolic tendon; the contact points are also uniformly distributed and located at the same parabola of AL tendon. The results of other parabolic tendon layouts (AS, BL, BS) are almost the same as the results of AL tendon. It is found that the influence coefficient Cpi is mainly related to the contact points n and slightly related to tendon layout. Table 2 indicates that, for the beam with only 3 contact points, i.e., without deviator along span, Cp = 1, EP effect is the same as N effect, because the loss of eccentricity is maximum, similar to the straight external tendon AB in Figure 1. As the number of contact points increases, the loss of eccentricityreduces, and the external tendon is gradually close to the unbonded tendon. The influence coefficient decreases from 1 to almost 0, and EP softening effect also decreases to almost 0. It conforms with the conclusion pointed out by Hamed and Frostig [17], i.e., zero effect of prestress compression of bonded or unbonded tendon on the natural frequencies. As the number of vibration modes increases, the distance between points of inflection (y'′=0) shortens, and the number of contact points between two inflection points decreases, which is equivalent to the reduction of contact points in the first antisymmetric mode; hence, the influence coefficient increases.Table 2
Relationship between influence coefficientCpi and the number of contact points n for straight tendon.
nAntisymmetric modeSymmetric modeCp1Cp2Cp3Cp1Cp2Cp3311111150.18910.910.2750.970.92870.0880.31610.1610.3350.93990.050.1890.3850.0980.2410.372110.0330.1250.2630.0650.1690.296130.0230.0880.1890.0460.1220.225170.0130.050.110.0260.0720.137210.0080.0330.0720.0170.0470.091Table 3
Relationship between influence coefficientCpi and the number of contact points n for AL tendon.
nAntisymmetric modeSymmetric modeCp1Cp2Cp3Cp1Cp2Cp350.1970.8980.9670.3040.9010.94470.0360.3350.960.150.3620.92190.0250.1270.4280.0780.2140.412110.0120.0450.2080.0370.1350.271130.010.020.1130.0130.0810.195170.0060.0110.0180.010.0190.093210.0030.0060.010.0050.010.037M①, M②, and M③ represent the methods in this paper, Shi et al. [19], and FEM, respectively. Generally speaking, the validity of theoretical results would be tested by experimental results. However, in the field and laboratory experiments, Saiidi et al. [5] found that the natural frequencies increase as the prestress compression increases, which is not consistent with the theoretical conclusion. For this disparity, they opined that the prestress compression causes closure of microcracks in the concrete, which increases the flexural stiffness and natural frequencies. In the experiment, the frictionless slip does not occur; even friction slip hardly occurs due to tight compression between deviator and external tendon under small vibration amplitude, so the experimental results may not validate the theoretical results. For this reason, the theoretical results are validated by FEM results in the present study.There exists the frictionless slip between the external tendon and deviator; the tension in external tendon is constant along whole length and is decided by whole beam deformation, which results in an additional difficult in FEM. Fang [22] proposed a simple method to solve this difficulty by adding a contact link element at angle bisector between external tendon and deviator as shown in Figure 8. The angle bisector position ensures the uniform increment of tension of external tendon, so the contact link element can simulate the frictionless slip. The length of contact link element is 1 mm, small enough to maintain the tendon layout unchanged; its axial stiffness is EtAt. The external tendon is subjected to only axial tension and is simulated by the link element. The deviator is simulated by the rigid beam element. The lumped mass method is adopted at two ends of the beam element. The beam is discretized into 16 beam elements with the same length 0.625 m, each lumped mass 0.5 t. The equation to calculate the natural frequencies is(23)ωi2δMφi=φi,where φi is the ith natural mode, M is the diagonal lumped mass matrix, and δ is the vertical flexibility matrix. The influence coefficients of Cpi and Pt0 in equations (21) and (22) are taken into consideration; hence, the stiffness matrix of beam element is(24)K=EIl3126l−126l6l4l2−6l2l2−12−6l12−6l6l2l2−6l4l2−CpiPt030l363l−363l3l4l2−3l−l2−36−3l36−3l3l−l2−3l4l2.CpiPt0 reduces the stiffness of beam element, lowers the frequency, and reflects EP softening. δ is computed by standard procedures of structural analysis and involves various parameters, such as EI, l, Cpi in beam stiffness matrix EtAt,lt and the layout of external tendon, the frictionless slip simulated by contact link element, etc. ω1∼ω14 are determined by MATLAB software for 14 lumped masses. The lumped mass method provides high accuracy in calculating low-order natural frequencies, so the high accuracy of ω1∼ω4 for 8 beams AS1∼BL2 listed in Table 4 is expected. Hence, the results of M③ in Table 4 can be used to validate the accuracy of the results of M① and M②.Figure 8
Contact link element.Table 4
The 1st∼4th natural frequency calculation results of externally prestressed continuous beam.
Beam codeM① (rad/s)M② (rad/s)M③ (rad/s)ω1ω2ω3ω4ω1ω2ω3ω4ω1ω2ω3ω4AS126.6544.31106.54134.8825.5840.87105.67134.1826.6444.26106.5134.8AS226.6146.77106.31134.6424.3740104.51133.9326.646.6106.28134.57AL126.6651.07106.6134.9525.8941.11105.97134.4426.6750.92106.56134.88AL226.6358.97106.43134.824.7840.29104.9133.5226.6458.49106.4134.74BS126.6344.28106.55134.8725.5840.87105.66134.1826.6444.22106.51134.79BS226.5846.7106.34134.6224.3740104.51133.1926.5746.55106.32134.55BL126.6250.71106.58134.8825.8941.1105.96134.4426.6150.56106.55134.8BL226.5658.34106.4134.6524.7740.29104.89133.5226.5758.19106.37134.57Note: M①, M②, and M③ represent the methods of this paper, Shi et al. [19], and FEM.The first and second frequencies of antisymmetric modes areω1, ω3; the first and second frequencies of symmetric modes are ω2, ω4; in Table 4, ω1∼ω4 values of M① are very close to the values of M③, indicating the accuracy of M① and M③. In M② method, the equivalent-area principle as shown in Figure 5 wrongly assumes the mutually strengthening effect of bending moment above axis and bending moment below axis on beam stiffness as mutually counteracting effect, which significantly underestimates the effect of external tendon. The equivalent-area principle approximately neglects the effect of external tendon for all frequencies. ωi of M② is slightly less than ωi of beam 0 without tendon in Table 5, the compression softening effect of Pt0 on frequencies is small, and the effect of external tendon is significantly underestimated. The error in M② method is obvious. In M① method, for symmetric mode ψ1≫ψi≈0,i≥2 in equation (17). The first-order tendon deformation is caused almost only by the first symmetric mode, very little by other symmetric modes; therefore, the stiffness of tendon EtAt mainly increases the frequency of the first symmetric mode (ω2) and slightly influences the frequency of the other symmetric modes (ω4, ω6...) in equation (22). The stiffness of tendon EtAt does not influence the frequencies of the antisymmetric modes (ω1, ω3...) in equation (21). Comparing ω2 values of M① among 8 beams, it shows that as the eccentricity and tendon area increase, ω2 values increase, but ω2 values are slightly related to the tendon A or B layout, so the influence of tendon layout linear transformation on the frequency is negligible. In M① method, the effect of external tendon on ω1, ω3 is zero, and the effect of external tendon on ω4 is negligible. In M② method, the effect of external tendon on all ωi is negligible due to the equivalent-area principle. Comparing results of M① and M②, it is found that ω1, ω3, and ω4 of M① are slightly greater than the results of M②, because M② method wrongly treats EP effect as N effect in equation (8); i.e., Cpi=1; ω2 of M① is obviously greater than ω2 of M②, because M① considers the effect of external tendon ψ1 in equation (21); M② method still neglects the effect of external tendon due to the equivalent-area principle.Table 5
ω2, ω4 of beams without prestress in M①.
Beam0AS1AS2AL1AL2BS1BS2BL1BL2ω2 (rad/s)41.744.3946.9251.1359.0744.3646.8650.7958.48ω4 (rad/s)135.11135.13135.16135.19135.26135.13135.15135.14135.18Note: beam 0 denotes beam without external tendon;ω1 = 26.69 rad/s, ω3 = 106.76 rad/s.Due to zero effect of external tendon onω1, ω3 and no prestress compression Pt0=0, ω1, ω3 for every beam are the same in Table 5. The frequencies of 8 beams AS1∼BL2 without prestress in Table 5 are quite close to M① frequencies in Table 4 with the prestress, indicating EP effect is negligible. In equation. (21), Cp1 is much less than 1 for one or more deviators (n≥5 in Tables 2 and 3) along beam span; usually Pt0≪EIπ/l2 the critical load of stability, EP effect, or Cp1Pt0 on ω1 is slight. For high-order frequencies ωi,i≥2, Cpi values increase, but the high-order critical loads of stability EIiπ/l2 increase more rapidly, and EP effect on high-order frequency is also slight. Equation (22) also presents the same conclusion. In a word, EP effect on all order frequencies is negligible for the beam with one or more deviators. Comparing ω1∼ω4 in M① Table 4 and in beam 0 Table 5, it is found that the effect of external tendon on ω2 is obvious, and the effect of external tendon on other frequencies is slight and negligible.


## 6. Conclusions

The conclusion mathematically rigorously proven by Hamed and Frostig [17] was introduced, i.e., zero effect of both bonded and unbonded prestress compression on the natural frequencies of beams as opposed to some research works. It is pointed out that there does exist EP effect on the natural frequencies of beams unlike the zero effect of both bonded and unbonded prestress compression. Miyamoto et al.’s [18] and Shi et al.’s [19] analytical solutions for externally prestressed simply supported and continuous beams were introduced and discussed. The energy method considering the second-order effects was adopted to obtain more accurate solution. The continuous beam frequency results of the energy method in this study, Shi et al.’s [19] analytical solution, and FEM were compared. The results obtained from this study are summarized below:There exist 4 problems in Miyamoto et al.’s [18] and Shi et al.’s [19] analytical solutions: 2 problems in the proportional assumption that the increment of tendon tension is proportional to the vibration displacement at midspan of beam; 1 problem in the equivalent-area principle that ΔMP is simplified into uniform ΔMPE; 1 problem in the treatment of EP as N.The influence coefficient of EP, i.e., EP softening effect, mainly depends on the number of contact points and slightly depends on tendon layout. Without deviator, the influence coefficient is equal to 1, and EP effect is the same asN effect. As the number of contact points increases from minimum value 3 to a large number in continuous beam, the deformation of external tendon is gradually close to unbonded tendon, the loss of external tendon eccentricity reduces, and the influence coefficient decreases from 1 to almost 0, which is consistent with the conclusion pointed out by Hamed [17]. The loss of external tendon eccentricity results in EP softening effect.For the beam without deviator, EP effect is the same asN effect, and it may not be negligible. For the beam with one or more deviators, EP effect is much less than N effect and is negligible. As one or more deviators always exist along the continuous beam span, therefore EP effect on frequency is negligible in continuous beam.The effect of tendon layout linear transformation on the frequency and influence coefficient is negligible.The first-order tendon deformation is not caused by all antisymmetric modes but is caused mainly by the first symmetric mode, and slightly by other symmetric modes. Therefore, the influence of external tendon on the first symmetric frequency (ω2) is obvious; the influence on other frequencies is negligible. As the eccentricity and area of tendon increase, the tendon improves the beam stiffness, and the first symmetric frequency (ω2) increases, but other frequencies almost remain unchanged.



---

*Source: 4397262-2020-12-28.xml*