Acessibilidade / Reportar erro

Analysis of the effect of the secondary moment on curved beams of full cross section

Análise do efeito do hiperestático de protensão em vigas curvas de seções transversais maciças

abstract

This paper analysis the effect of the secondary moment on curved beams using the equivalent nodal loads method. A case study was carried out applying the equivalent nodal loads method to a two-span prestressed curved beam, using a curved finite element as the structural model, analyzing the effect of the secondary moment for each case and comparing it with its equivalent in straight beam. It was found that the stiffness parameters EI and GJ influence the secondary moment. The results demonstrates that the beam opening angle reduces the effect of secondary moment, and that the greater the angle, the greater is the reduction that occurs in its secondary moment compared to its equivalent straight beam.

Keywords:
prestressed concrete; secondary moment; continuous beam; curved beam; curved finite element


resumo

Neste artigo estuda-se o efeito do hiperestático de protensão em vigas curvas, por meio do método das cargas nodais equivalentes. É realizado um estudo de caso, aplicando-se o método das cargas nodais equivalentes em uma viga curva protendida de dois vãos, utilizando-se como modelo estrutural um elemento finito curvo, analisando-se o efeito do hiperestático de protensão para cada caso e comparando-o com o seu equivalente em viga reta. Constatou-se analiticamente que os parâmetros de rigidez EI e GJ influenciam na obtenção do hiperestático de protensão. Os resultados demonstram que o ângulo de abertura da viga reduz o efeito do hiperestático de protensão e que, quanto maior esse ângulo, maior é a redução que ocorre em seu hiperestático de protensão em comparação ao seu equivalente em viga reta.

Palavras-chave:
concreto protendido; hiperestático de protensão; viga contínua; viga curva; elemento finito curvo

1 INTRODUCTION

Prestressed concrete has increasingly become the solution for the design of certain structures. Its use is due to the great structural advantages obtained, managing to overcome large spans, reduce the dead weight of the structures, and lower the cost of foundations.

One of the most frequent uses in prestressed concrete is in slabs and bridge beams, as the geometry of bridges often requires curved sections that make it difficult to design, considering that curved elements have peculiar characteristics that complicate structural analysis. However, the study of prestressed concrete in curved beams is still little known in scientific and academic environment, especially in the case of statically indeterminate structures. Among the tests performed on curved prestressed elements, those of Zhang et al. [11 N. Zhang, C. C. Fu, and H. Che, "Experiment and numerical modeling of prestressed concrete curved slab with spatial unbonded tendons," Eng. Struct., vol. 33, no. 3, pp. 747–756, Mar 2011.], Amorn et al. [22 W. Amorn, C. Y. Tuan, and M. K. Trados, "Curved, precast, pretensioned concrete I-girders bridges," Civ. Eng. Fac. Publ., vol. 11, no. 11, pp. 48–66, Dec 2008.] and Choi et al. [33 J. Choi, C. R. Woods, T. D. Hrynyk, and O. Bayrak, "Behavior of curved post-tensioned concrete structures without through-thickness reinforcement," ACI Struct. J., vol. 114, no. 114, pp. 983–994, Aug 2017.] stand out.

The calculation of prestressing in statically indeterminate structures must take into account the effect of secondary moments. Cohn and Frostig [44 M. Z. Cohn and Y. Frostig, "Inelastic behavior of continuous prestressed concrete beams," J. Struct. Eng., vol. 109, no. 10, pp. 2292–2309, Oct 1983.] found that the secondary moments is of relevant importance in the ultimate limit state, which makes it necessary to obtain its magnitude more accurately.

Some methods proposed for the determination of the secondary moment are mentioned: Kong [55 J. Kong, "A practical method for the calculation of secondary prestress moments," Hong Kong Inst. Eng., vol. 11, no. 4, pp. 49–52, Dec 2004.] developed a method for calculating the secondary moments based on the influence line; Raju [66 K. Raju, Prestressed Concrete, 4th ed. India: McGraw-Hill, 2013.] presents the usual way of calculating the secondary moment, using the equivalent distributed load method and the flexibility method; Cunha [77 T. Cunha, “Torsional moment analysis in prestressed curved beam,” M.S. thesis, Civ. Eng. Postgrad. Program, Fluminense Fed. Univ., Rio de Janeiro, RJ, Brazil, 2018.] presented a method that is based on the concept of equivalent nodal loads.

In this article, the effect of the secondary moment on statically indeterminate curved beams is studied, comparing the result with its straight beam equivalent, in addition to demonstrating the influence that the curvature of the structure has in determining the equivalent nodal loads.

2 EQUIVALENT NODAL LOADS

For the calculation of the secondary moment in curved beams, the concept of equivalent nodal load developed by Cunha [77 T. Cunha, “Torsional moment analysis in prestressed curved beam,” M.S. thesis, Civ. Eng. Postgrad. Program, Fluminense Fed. Univ., Rio de Janeiro, RJ, Brazil, 2018.] is used. The determination of equivalent nodal loads in prestressed structures is carried out similarly on to obtaining values of equivalent nodal load in the matrix analysis of structures. Soriano [88 H. L. Soriano, Análise de Estruturas – Formulação Matricial e Implementação Computacional. Rio de Janeiro: Ciência Moderna Ltda, 2005.], Martha [99 L. F. Martha, Análise de Estruturas: Conceitos e Métodos Básicos, 1a ed. Rio de Janeiro: Elsevier, 2010.] and Przemieniecki [1010 J. S. Przemieniecki, Theory of Matrix Structural Analysis. New York: McGraw-Hill, 2012.] demonstrate how load values are obtained for external loads and charts with their respective values.

A curved prestressed element with fixed ends is considered and analyzed as a beam element, thus the flexibility method is applied to obtain the load values. Figure 1 illustrates the prestressed bar element whose layout has a generic configuration, in order to make the method applicable to any situation. The load value to be obtained is the equivalent bending moment, as torsion due to prestressing loading is disregarded.

Figure 1
Curved fixed ends bar subjected to prestressing.

According to Cunha [77 T. Cunha, “Torsional moment analysis in prestressed curved beam,” M.S. thesis, Civ. Eng. Postgrad. Program, Fluminense Fed. Univ., Rio de Janeiro, RJ, Brazil, 2018.], the vertical reaction that appears in the fixed support should be disregarded from the calculations, because only the effect of the bending moment due to prestressing generates the secondary moment.

The deduction of equivalent nodal loads is obtained through the flexibility method (Figure 2). Considering that it is a curved element, when applying the unitary bending moment, the torsion moment occurs. The diagram of primary bending moment due to prestressed is given in a generic configuration, with the goal of making the method usable for any problem. In this study the warping of the section is disregarded.

Figure 2
Application of the flexibility method for obtaining equivalent nodal loads.

Figure 3 shows a differential curved bar element used as a reference for deducting the differential equations of equilibrium, where the elements of higher order are neglected. The deduction of all differential equations of equilibrium for curved beam is presented in detail in Cunha [77 T. Cunha, “Torsional moment analysis in prestressed curved beam,” M.S. thesis, Civ. Eng. Postgrad. Program, Fluminense Fed. Univ., Rio de Janeiro, RJ, Brazil, 2018.] and Montánari [1111 I. Montánari, “Protensão em vigas curvas,” Ph.D. dissertation, Fed. Univ. São Carlos, São Carlos, SP, Brazil, 1980.].

Figure 3
Differential curved bar element.

The curved beam equilibrium equation that relates the bending moment with the torsion moment can be obtained from the sum of the moments in x direction (normal axis to the section of the element) and from geometric relationship ds=rdθ, obtaining:

dT ds + M y r = t (1)

Considering that only the effect of prestressing on the structure is studied, the external torsion loading is null, so t=0, thus:

dT ds + M y r = 0 (2)

The bending moment is unitary and constant along the structure (Figure 2). Replacing its value and integrating it into Equation 2, the torsion moment acting on the structure is obtained due to virtual unitary load.

T 1 = L r (3)

Equations 5 and 6 depict the flexibility coefficients of the flexibility method

Δ 11 = 0 l T 1 2 GJ ds + 0 l M 1 2 EI ds (4)
Δ 11 = L 3 r 2 GJ + L EI (5)
Δ 10 = 0 L m m ¯ EI = 0 L Pz s ds EI (6)

that gathered in the compatibility equation

Δ 10 + x 1 Δ 11 = 0 (7)
L 3 r 2 GJ + L EI x 1 + 0 l Pz s 1 ds EI = 0 (8)

Results in

x 1 = r 2 GJ L L 2 EI + r 2 GJ 0 L Pz s ds (9)

In equation 9 using the following notation we have

μ = r 2 GJ L 2 EI + r 2 GJ (10)

where m = bending moment due to prestressing; m¯=bendingmomentduetounitload; P=prestressingforceapplied; L=spanlength; r=beamradiusofcurvature; EI=spanflexuralstiffness; GJ=spantorsionalrigidity; Δ10=displacementatpoint1duetotheloadappliedatpoint0; Δ11=displacementatpoint1duetounitloadatpoint1; zs=heightoftheprestressingcable.

The value μ shown in Equation 10 is called by Cunha [77 T. Cunha, “Torsional moment analysis in prestressed curved beam,” M.S. thesis, Civ. Eng. Postgrad. Program, Fluminense Fed. Univ., Rio de Janeiro, RJ, Brazil, 2018.] as a reduction factor due to curvature. Thus, the equivalent nodal load is given by

M 1 = μ 0 L Pz s ds L (11)

Equation 11 is valid only when the prestressing force is centered, that is, when there is no resulting eccentricity in the section that generates torsion. If there is an eccentricity, the effect of torsion due to prestressing should be considered in determining the load values.

Figure 4 shows the equivalent nodal loads due to prestressing for fixed ends situation.

Figure 4
Prestressing equivalent nodal load.

The orientation of the nodal load is a function of the 0LPzsds portion, which, if positive, the orientation is the same as shown in Figure 4. If on the contrary, the orientation goes in the opposite direction.

The 0Lzsds portion is the area between the centroid of the section and the prestressing profile along its entire length. The strategy for using this method consists of dividing the beam into sections, calculating the area between the prestressing profile and the centroid in each section and multiplying it by the average prestressing force of the section, considering prestressing losses.

It is observed that the equivalent nodal load is presented as a parameter to assess the influence that the secondary moment has on the structure. As this is the calculation of the area between the centroid of the section and the profile, if the sum of the areas is null, the secondary moment does not occur, Raju [66 K. Raju, Prestressed Concrete, 4th ed. India: McGraw-Hill, 2013.] defines this characteristic as a concordant prestressing path.

3 STIFFNESS MATRIX OF CURVED FINITE ELEMENT

The stiffness matrix of a curved element of constant section and curvature to be used was deduced by Palaninathan and Chandrasekharan [1212 R. Palaninathan and P. S. Chandrasekharan, "Curved beam element stiffness matrix formulation," Comput. Struc., vol. 21, no. 4, pp. 663–669, 1985.], based on Lee’s studies [1313 H.-P. Lee, "Generalized stiffness matrix of a curved-beam element," AIAA J., vol. 7, no. 10, pp. 2043–2045, Oct 1969.]. The original element has six degrees of freedom per node (Figure 5).

Figure 5
Curved beam element, Palaninathan and Chandrasekharan [1212 R. Palaninathan and P. S. Chandrasekharan, "Curved beam element stiffness matrix formulation," Comput. Struc., vol. 21, no. 4, pp. 663–669, 1985.].

However, the number of degrees of freedom was reduced to three, being the vertical force, the bending moment and the torsion moment, transforming it into a curved grid finite element (Figure 6).

Figure 6
Curved grid finite element.

The reduction in the number of degrees of freedom aims at simplifying the structural calculation, since the inversion of the 12 × 12 stiffness matrix would result in a large computational work, in addition to the other variables involved.

The internal forces and the moments acting at point P are expressed in values of the forces from node i as follows

V' = V 2 (12)
T' = T 1 cos M 3 sen V 2 r 1 cos (13)
M' = T 1 sen + M 3 cos + V 2 rsen (14)

The strain energy in the element, disregarding the section warping effect, is given by

U = r 2 0 T' 2 GJ + M' 2 EI d (15)

Using Castigliano's theorem we obtain

U 1 = U T' = r GJ 0 T' T' T' d + r EI 0 M' M' T' d (16)
U 2 = U V' = r GJ 0 T' T' V' d + r EI 0 M' M' V' d (17)
U 3 = U M' = r GJ 0 T' T' M' d + r EI 0 M' M' M' d (18)

Developing the integrals and applying Castigliano’s theorem we have the following system of equations

u 1 u 2 u 3 = a 11 a 12 a 13 a 12 a 22 a 23 a 13 a 23 a 33 V' T' M' (19)

By developing the integral of the strain energy and using the following equation to simplify the values obtained, results

λ = EI GJ (20)
Δ AA = 1 EI r 2 2 λCr + Br r 2 2 B λV r 2 2 λS D r 2 2 B λV r 2 λN + B r 2 λD D r 2 2 λS D r 2 λD D r 2 λB + N (21)

Where

N = θ + sen 2 θ 2 (22)
B = θ sen 2 θ θ (23)
C = 3 θ + sen 2 θ θ 4 sen θ (24)
S = 3 4 cos θ + cos 2 θ 4 (25)
V = 2 sen θ θ sen 2 θ 2 (26)
D = cos 2 θ 2 1 2 (27)

The stiffness matrix of the element is given by the joining of matrices KAA,KAB,KBA e KBB, thus

k = K AA K AB K BA K BB (28)

The KAA matrix is obtained by inverting the flexibility matrix ΔAA

K AA = Δ AA 1 (29)

The matrix T relates the forces of the final node j to the forces of the initial node i

T = 1 0 0 r 1 cos θ cos θ sen θ rsen θ sen θ cos θ (30)

The KBA matrix is obtained by

K BA = TK AA (31)

In a similar way, the flexibility matrix is obtained

Δ BB = 1 EI r 2 2 λCr + Br r 2 2 B λV r 2 2 λS D r 2 2 B λV r 2 λN + B r 2 λD D r 2 2 λS D r 2 λD D r 2 λB + N (32)

The KBB matrix is obtained by inverting the flexibility matrix

K BB = Δ BB 1 (33)

being that

K AB = TK BB (34)

With this data, the stiffness matrix of the element can be assembled.

4 NUMERICAL EXAMPLE

The effect of the secondary moment on a continuous beam with two spans of constant curvature is studied and the results are compared with straight continuous beam of length equivalent to the radius and angle adopted in each case. The support conditions for both beams are of restriction to vertical and horizontal reaction. The study is carried out on a beam with a radius of 40 m, varying its angle from 10 ° to 10 ° increments, up to 90 °.

Mathcad software was used to calculate the example, basing its implementation routine on studies carried out by Vaz [1414 L. E. Vaz, Método dos Elementos Finitos em Análise de Estruturas. Rio de Janeiro: Elsevier, 2011.] and Soriano [88 H. L. Soriano, Análise de Estruturas – Formulação Matricial e Implementação Computacional. Rio de Janeiro: Ciência Moderna Ltda, 2005.], using as a basis the matrix analysis of structures and the finite element method.

In order to simplify the study of the effect of the secondary moment on curved beams, it is settled that, in all cases, the lengths of the spans are equal, thus, L1=L2.

The prestressing path is straight with an applied load of 1000 kN. In order to facilitate the study, a constant prestressing loss of 20% was also admitted along its entire path and its eccentricity e=40cm. The portion corresponding to the prestressing of the equivalent nodal load (Equation 11) is the following value.

0 L 1 Pz s ds L 1 = 0 L 2 Pz s ds L 2 = 800 x 0.4 = 320 kN m (35)

Equation 35 is constant in all cases, since there is no change in the eccentricity or in load path profile. However, the equivalent nodal load (Equation 11) changes in all cases as a function of the reduction factor due to the curvature μ.

The modified curved finite element of Palaninathan and Chandrasekharan [1010 J. S. Przemieniecki, Theory of Matrix Structural Analysis. New York: McGraw-Hill, 2012.] was used for the calculation of the curved beam and for the straight beam the bar element.

Two span continuous curved beam

Figure 7 shows the section data.

Figure 7
Prestressed cross section, measured in cm.

The ratio between flexural stiffness and torsional stiffness is equal to 1 for all cases presented.

The boundary conditions adopted in the structure are shown in Figure 8, which are the restriction to vertical displacement in all nodes of the structure.

Figure 8
Two-span curved continuous beam.

Table 1 presents the results obtained in each case and Figure 9 shows, in a generic way, the effect of the secondary moment in the curved beam.

Table 1
Two-span curved continuous beam subjected to centered prestressingr=40m.
Figure 9
Effect of the secondary moment in curved beam.

5 ANALYSIS OF RESULTS

It is assessed that the secondary moment in a curved structure has its value reduced as function of the beam opening angle. For small opening angles, between 10 ° and 30 °, the reduction that occurs in its secondary moment is insignificant, reaching a maximum of 7.2% when compared to its straight equivalent beam. Garrett and Cochrane [1515 R. J. Garrett and R. A. Cochrane, "The analysis of prestressed beams curved in plan with torsional restraint at the supports," Struct. Eng., vol. 48, no. 3, pp. 128–132, Mar 1970.] arrived at similar analytical results, concluding that a curved beam with an opening angle of less than 30° can be calculated as a straight beam with good precision of the results. Tests carried out by Amorn et al. [22 W. Amorn, C. Y. Tuan, and M. K. Trados, "Curved, precast, pretensioned concrete I-girders bridges," Civ. Eng. Fac. Publ., vol. 11, no. 11, pp. 48–66, Dec 2008.] and compared by a straight beam model in Cunha [77 T. Cunha, “Torsional moment analysis in prestressed curved beam,” M.S. thesis, Civ. Eng. Postgrad. Program, Fluminense Fed. Univ., Rio de Janeiro, RJ, Brazil, 2018.] prove that curved beams with a small opening angle, including prestressed ones, can be approximated by a straight beam with a negligible difference in result, which corroborates the results obtained.

The radius of the structure does not influence the obtaining of factor μ, as this depends directly on the opening angle of the structure and the stiffness parameters EI and GJ. The reduction factor μ has a greater influence on the structure when its opening angle increases, because the larger its opening, the greater the interaction between the bending moment and the torsion moment, as can be seen in Equation 14.

The secondary moment has its value reduced by down to approximately 42% (90° angle), showing that in curved structures, with a high opening angle, the secondary moment has significantly less influence than in rectilinear structures.

It should be noted that the results obtained took into account that the flexural stiffness has a numerical value equal to the torsional stiffness, it can be deduced that if λ>1, that is, the flexural stiffness is superior to the torsional stiffness, a reduction of the factor μ occurs and if λ<1 result in its increase. It should be highlighted that the stiffness parameters, in addition to being important in determining the μ factor, are also relevant in determining the coefficients of the curved finite element stiffness matrix.

It is noteworthy that this study did not take into account the distortion due to shear, as well as the effect of prestressing in thin-walled sections, whether open or closed. The results obtained are valid only for full sections, where shear has little influence on the distortion of the section.

6 CONCLUSIONS

In curved beams, the reduction factor due to curvature μ is the main cause of the decrease in secondary moment. It has been proven that beams with a high opening angle have a greater reduction in their secondary moment. It is worth noting that the relationship between flexural stiffness and torsional stiffness also contribute to the variation of the μ factor intensity, it was found that if the flexural stiffness is greater than the torsional stiffness it results in a reduction of the μ factor and if otherwise, it increases. Therefore, it is inferred that the main variables for the calculation of the secondary moment, in a curved beam, are the opening angle of the structure and its stiffness parameters EI and GJ.

By the method of equivalent nodal load it is verified that, in prestressed beams whose tracing is straight, the secondary moment occurs accentuated, since the area between the centroid and the prestressing cable is larger than in beams whose tracing is variable.

The results conveyed that in curved beams with a small opening angle (between 10 ° and 30 °) the difference that occurs between the secondary moment of the curved beam and its equivalent of straight beam is negligible in the centered prestressing, reaching less than 7.2%, even with reduction factor applied. There is a significant decrease in secondary moment, mainly in beams whose opening angle is high.

The curved element used for the beam calculation is more complicated to be computationally implemented. While the stiffness matrix of the structure stays smaller than if discretization was used in bar elements, the number of variables involved ends up significantly impairing its computational performance.

The proximity of the results obtained, when comparing the result of straight beam with that of curved beam, allows us to consider it reasonable to discretize the structure instead of using curved elements for small opening angles, being recommended to use curved element for higher opening angles.

Ultimately, it is worth mentioning that the reduction factor must be used in curved beams whose approach has been modeled from straight bar elements. The fact that the modeling occurs with bar elements does not exempt the application of the reducing factor, which is a factor that seeks to consider the interaction between the torsion moment and the bending moment.

  • Financial support: This study was financed partially by Coordenação de Aperfeiçoamento de Pessoal de Nível Superior - Brasil (CAPES) Finance Code 001.
  • How to cite: T. C. Silva and E. S. Sánchez Filho, “Analysis of the effect of the secondary moment on curved beams of full cross section,” Rev. IBRACON Estrut. Mater., vol. 14, no. 2, e14207, 2021, https://doi.org/10.1590/S1983-41952021000200007

CITATIONS

  • 1
    N. Zhang, C. C. Fu, and H. Che, "Experiment and numerical modeling of prestressed concrete curved slab with spatial unbonded tendons," Eng. Struct., vol. 33, no. 3, pp. 747–756, Mar 2011.
  • 2
    W. Amorn, C. Y. Tuan, and M. K. Trados, "Curved, precast, pretensioned concrete I-girders bridges," Civ. Eng. Fac. Publ., vol. 11, no. 11, pp. 48–66, Dec 2008.
  • 3
    J. Choi, C. R. Woods, T. D. Hrynyk, and O. Bayrak, "Behavior of curved post-tensioned concrete structures without through-thickness reinforcement," ACI Struct. J., vol. 114, no. 114, pp. 983–994, Aug 2017.
  • 4
    M. Z. Cohn and Y. Frostig, "Inelastic behavior of continuous prestressed concrete beams," J. Struct. Eng., vol. 109, no. 10, pp. 2292–2309, Oct 1983.
  • 5
    J. Kong, "A practical method for the calculation of secondary prestress moments," Hong Kong Inst. Eng., vol. 11, no. 4, pp. 49–52, Dec 2004.
  • 6
    K. Raju, Prestressed Concrete, 4th ed. India: McGraw-Hill, 2013.
  • 7
    T. Cunha, “Torsional moment analysis in prestressed curved beam,” M.S. thesis, Civ. Eng. Postgrad. Program, Fluminense Fed. Univ., Rio de Janeiro, RJ, Brazil, 2018.
  • 8
    H. L. Soriano, Análise de Estruturas – Formulação Matricial e Implementação Computacional. Rio de Janeiro: Ciência Moderna Ltda, 2005.
  • 9
    L. F. Martha, Análise de Estruturas: Conceitos e Métodos Básicos, 1a ed. Rio de Janeiro: Elsevier, 2010.
  • 10
    J. S. Przemieniecki, Theory of Matrix Structural Analysis. New York: McGraw-Hill, 2012.
  • 11
    I. Montánari, “Protensão em vigas curvas,” Ph.D. dissertation, Fed. Univ. São Carlos, São Carlos, SP, Brazil, 1980.
  • 12
    R. Palaninathan and P. S. Chandrasekharan, "Curved beam element stiffness matrix formulation," Comput. Struc., vol. 21, no. 4, pp. 663–669, 1985.
  • 13
    H.-P. Lee, "Generalized stiffness matrix of a curved-beam element," AIAA J., vol. 7, no. 10, pp. 2043–2045, Oct 1969.
  • 14
    L. E. Vaz, Método dos Elementos Finitos em Análise de Estruturas. Rio de Janeiro: Elsevier, 2011.
  • 15
    R. J. Garrett and R. A. Cochrane, "The analysis of prestressed beams curved in plan with torsional restraint at the supports," Struct. Eng., vol. 48, no. 3, pp. 128–132, Mar 1970.

Edited by

Editors: Leandro Mouta Trautwein, José Luiz Antunes de Oliveira e Sousa, Guilherme Aris Parsekian.

Publication Dates

  • Publication in this collection
    25 Jan 2021
  • Date of issue
    2021

History

  • Received
    10 July 2018
  • Accepted
    01 July 2020
IBRACON - Instituto Brasileiro do Concreto Instituto Brasileiro do Concreto (IBRACON), Av. Queiroz Filho, nº 1700 sala 407/408 Torre D, Villa Lobos Office Park, CEP 05319-000, São Paulo, SP - Brasil, Tel. (55 11) 3735-0202, Fax: (55 11) 3733-2190 - São Paulo - SP - Brazil
E-mail: arlene@ibracon.org.br