Mathematical Problems in Engineering
Volume 2016 (2016), Article ID 2465025, 12 pages
http://dx.doi.org/10.1155/2016/2465025

Research Article

A Trigonometric Analytical Solution of Simply Supported Horizontally Curved Composite I-Beam considering Tangential Slips

Qin Xu-xi, Liu Han-bing, Wu Chun-li, and Gu Zheng-wei

College of Traffic, Jilin University, Changchun 130025, China

Received 27 March 2016; Revised 31 July 2016; Accepted 14 August 2016

Academic Editor: Leonid Shaikhet

Copyright © 2016 Qin Xu-xi et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

This paper presents an analytical solution of the simply supported horizontally composite curved I-beam by trigonometric series considering the effect of partial interaction in the tangential direction. Governing equations and boundary conditions are obtained by using the Vlasov curved beam theory and the principle of minimum potential energy. The deflection functions and the Lagrange multiplier functions are expressed as trigonometric series to satisfy the governing equations and the simply supported constraints at both ends. The numerical results of deflections and forces which are obtained by this method are compared with both FEM results and experimental results, and the inaccuracy between the analytical solutions in this paper and the FEM results is small and reasonable.

1. Introduction

A composite beam can be defined as “partially composite beam” when the number of shear connectors is less than the required number for fully composite design; therefore, the interface shear force is limited by the strength of shear connectors. In contrast to fully composite beam, slip between layers can be significant and results in a decrease in the elastic stiffness of partially composite beam. Composite beam exhibiting partial shear interaction will stand larger deflection than the beam exhibiting full shear interaction. By assuming complete shear interaction, the calculation of deflection for partial shear interaction beam is maybe underestimated. Because serviceability issues often govern the structural design of composite section, the accurate calculation of deflection is critical. The partial interaction is applied not only in steel-concrete composite beam, but also in other types of composite beams, such as layered wooden beams, wood-concrete floor systems, and other multilayered laminated composite structures [1, 2].

Earlier studies on behavior of the partially composite beam are mostly focused on the straight composite beam. The first paper dealing with the analysis of composite beam with partial interaction has been completed by Newmark et al. [3]. After that, Goodman and Popkov [4, 5] have conducted the analytical and numerical research on the relative slip between layers and found that the relative slip between layers has a significant effect on the overall characteristics of the composite beam with the reduction of shear connectors’ stiffness. Girhammar et al. [6, 7] have applied a partial shear interaction theory for composite beam subjected to static loads and dynamic loads. Wang [8] has developed a method to calculate deflection of partially composite beam based on the stiffness of the shear connectors. Dall’Asta [9] has developed a three-dimensional theory for composite beam with partial shear interaction dealing with combination of bending in the symmetry plane, torsion, and transverse bending in the plane parallel to the shear connector interface. Nie and Cai [10] have studied the effects of shear slip on the deflection of steel-concrete composite beam. Ranzi and Bradford [11] have presented analytical solutions for time dependent behavior of partially composite beam. Liu et al. [12] have found out the solution of shearing slip for steel-concrete composite beam under the concentrated load. Campi and Monetto [13] have presented a new formulation to analyze two-layer linearly elastic Timoshenko beam with interlayer slip. In the aspect of numerical simulations studies, different kinds of numerical and finite element formulations for the analysis of composite beam with interlayer slip have been suggested [1, 14–19].

Many significant researches have been accomplished in regard to the behavior of straight partially composite beam. In the same time, various scholars have done researches on the curved beam theory. One of the earliest works dealing with the stability behaviors of curved beam has been put forward by Vlasov [20]. After that, scholars covered different extensions and enhancements to the Vlasov model [21–26]. These researches are helpful to the development of the composite beam theory. However, very little literature has focused on the curved partially composite beam such as horizontally composite curved steel I-beam bridge. Thevendran et al. [27] and Shanmugam et al. [28] have conducted experiments on the steel-concrete composite curved beam to investigate the ultimate load behavior. In their study, the finite element software ABAQUS was used to analyze the behavior of test specimens. Full composite action between steel beam and the concrete slab was assumed. The results of deformations, stress distributions, and ultimate strengths obtained by finite element analysis were found to be in good agreement with the experimental results. After that, Topkaya et al. [29] have conducted experimental and numerical studies to establish the behavior of composite curved bridge during construction. In their study, two FEM models were established by different kinds of software to predict the behavior of the curved steel trapezoidal box-girder, and the authors have drawn a conclusion that the reasonable finite element model is able to accurately capture girder behavior during construction. Giussani and Mola [30] have developed an analytical equation for elastic composite beam curved in-plane with the long term behavior. But the partial interaction between the steel girder and concrete slab was not considered in the study. Erkmen and Bradford [31] have solved the equation of the composite curved beam considering the two-layer partial interaction by providing a highly efficient 3D beam finite element. The results demonstrate that the developed formulation is accurate and effective in capturing the behavior of composite beams curved in-plane.

Even though some previous researches of these evaluation methodologies have been accomplished, there is still lack of the fundamental understanding of the system-level behavior on the overall performance of composite curved beam with partial shear interaction. Besides, although exact solutions can be obtained by the 3D finite element model, it is very complex and time consuming. Therefore, the FEM analysis procedures are not suitable for the initial design. This study aims to provide an analytical theory of the horizontally composite curved beam considering the partial interaction in tangential direction. The beam is assumed to be statically determinate with a constant radius of curvature along the longitudinal axis. Governing equations and boundary conditions are obtained by using both the Vlasov curved beam theory and the energy variation principle. The undetermined vertical deflection, torsional deflection, and Lagrange multipliers will be approximated by Fourier series to solve the governing equations of the partial interaction composite beam theory in the procedures. The numerical results of deflections and forces obtained by using proposed theory are presented and compared with FEM results and experimental results. Comparison results show that the calculation can be easily and accurately handled which is a great advantage of this method.

2. Basic Assumptions and Conditions

This study is on a horizontally composite curved I-beam where the model is shown as in Figure 1, and the noteworthy features of this research are shown as follows:

(1)The slab and I-girder are linear-elastic different materials, and the cross-sections made of both materials are rigid in their plane. The effects of shear deformation, warping deformation, and distortion deformation are neglected in this research. Structural analysis of the beam is based on the Vlasov curved beam theory (for each part of the beam).(2)The interlayer connectors between the slab and I-girder are flexible, and they are continuous in tangential direction and rigid in radial direction. The load-slip behavior (per unit length) of connectors in tangential direction is described in a linear-elastic range with a constant slip modulus .(3)The uplift between the slab and I-girder is neglected. The radius of curvature is constant along the beam.

Figure 1: Model of a horizontally composite curved I-beam.

The deflections in tangential direction, vertical direction, and radial direction (-direction, -direction, and -direction) are indicated as , , and , respectively. The I-girder has the same torsional deflection and vertical deflection as the slab has. In this paper, subscript “1” refers to the cross-section of the slab and subscript “2” refers to the cross-section of I-girder, such that the deflections of the slab and I-girder can be represented as and , respectively. The geometric relationships between strains and deflections may be written as

where is the axial strain, is the curvature in -direction, is the curvature in -direction, and is the curvature in -direction. Where , .

The constitutive equations are as follows:

where , , , and are internal axial force in -direction, internal bending moment in -direction, internal bending moment in -direction, and internal torsion moment in -direction, respectively.

The microunit’s force diagram of the horizontally composite curved I-beam is shown in Figure 2, which is under the external distribution force and torsion moment . The balance equations of the slab are

The balance equations of the I-girder are

where is the distributed bending moment caused by shear force , is distributed torque caused by shear force , is the shear force in -direction, and is the shear force in -direction. For the shear forces , , and , there are

Figure 2: A microunit’s force diagram of the horizontally composite curved I-beam.

By eliminating and from the first three equations of (3) and (4), we can get

When there are no axial force and bending moment in -direction of the I-beam caused by external forces, we can get

3. Equilibrium at the Interface

At the interface, considering the deflections described in Figure 3, shear force described in Figure 2 can be written as

where and is the amount of the deflection at the interface in the tangential direction. So, the equilibrium at the interface in the tangential direction can be written as

Figure 3: Deflections of differential elements for the beam in - plane.

Considering deflection described in Figure 4, equilibrium at the interface in the radial direction can be written as

Figure 4: Deflection of the composite curved I-beam in - plane.

In order to solve the problem under consideration, it is convenient to rearrange (9) and (10) as detailed in what follows. Firstly, taking the derivative of (9) one time with respect to component, then eliminating and by inserting the derivative equation and (10) into the first formula of (1), and at last using the first and second formulas of (2), one can obtain the following:

Analogously, taking the derivative of differential equation (9) one time and taking the derivative of (10) two times with respect to component firstly and then eliminating and by inserting these two derivative equations into the second formulas of (1), using (2), one can obtain the following:

where and .

4. Problem Formulation

The governing differential equations of the beam will be derived by the principle of minimum potential energy. The potential energy of the beam takes the following form:

is the bending strain energy due to the internal bending moments and .

where .

is the torsional strain energy due to the internal torsion moments and .

where .

is the axial strain energy due to the internal axial forces and .

is the bending strain energy due to the internal bending moments and .

is the strain energy due to the connector deflections.

is the potential energy due to the external loading.

where , , and are the bending moment, shear force, and torsion moment of the beam, respectively. Using the Lagrange multiplier method, the principle of minimum potential energy in augmented form now may be written as

where the parameters and are the Lagrange multipliers and and are the equilibrium conditions at the interface corresponding to (11) and (12), respectively. In (20), , , , , , and are all independent variables. The variation of is

where and are calculated from the variation of , which are given in Appendix A. Each term in (21) must be identically zero. It makes , , , , , , and ; that is, the undermined Lagrange multipliers and are

Equation (22) notes that the Lagrange multipliers and the deflections ( and ) are related. So, there will be four rather than six independent variables in this problem. We can take , , , and as independent variables. The rest of governing equations are , , , and , which can be rearranged and rendered in terms of matrix form, as shown below.

We can get the pertaining boundary conditions:

It is evident that (23) is tenth-order ordinary differential equation. We have ten boundary conditions for the beam, that is, five boundary conditions for each end. For the boundary conditions (see (24)) to hold, the following alternatives at each end shown as in Table 1 are conducted.

Table 1: Common support conditions.

Type
The boundary conditions

5. Analytical Solution by Trigonometric Series

In this section, a closed form solution of the simply supported horizontally composite curved I-beam is obtained. The deflection functions and the Lagrange multiplier functions are expressed as undetermined coefficients and known trigonometric series to satisfy the governing equations and the rigid torsion constraints at each end. The deflection and the Lagrange multiplier are assumed to be the following form:

where , , , and are unknown Fourier coefficients to be determined for each and is span length. The applied distributed load and are expanded in single trigonometric series as

where and are the Fourier coefficients. Multiplying both sides of (26) by each trigonometric series and then integrating them over , finally, the trigonometric coefficients and can be readily determined with the orthogonality properties. Substituting (25)-(26) into (23) and solving (23), the unknown parameters , , , and can be determined. The unknown Fourier coefficients are given as follows:

where , , and are the relevant coefficients which are given in Appendix B.

6. Numerical Example

A steel-concrete composite curved I-beam calculated by Thevendran et al. [27] is taken as a numerical example in order to show the accuracy and reliability of the present trigonometric solution of the beam. The beam is simply supported at both ends as shown in Figure 1. And the beam is subjected to vertical load at the midspan section. The material properties and dimensions of the curved I-beam (SP4 beam) are shown in Table 2.

Table 2: Material properties and dimensions of the curved I-beam.

Young’s modulus
(GPa)
Poisson’s ratio
Density
(kg/m)
Width
(mm)
Thickness
(mm)
Central angle
Radius of curvature
(m)

Steel girder
206
0.3
7850
Flange: 332
Web: 172
Flange: 13
Web: 8

24

Concrete slab
26
0.27
2400
1500
100

The finite element software ANSYS is used to model the described structure (Figure 5). In the FEM model, both the concrete slab and the steel I-girder are modeled as space beam elements and the connection between the slab and I-girder is simulated by multiple-point constraints (MPC) which are modeled by two rigid links connected through interface nodes between the slab and I-girder. The spring elements are used in the tangential direction to allow for the possibility of movement. Coupling degrees of freedom are used to prevent the uplifting issue.

Figure 5: Finite element model diagram of the composite curved I-beam.

The vertical deflection , the tangential deflection , the radial deflection , and the torsional deflection are restrained to satisfy the simply supported constraints at both ends. Although the BEAM4 element in ANSYS cannot consider the warping effect, the element meets the theory of this paper quite well. So, 120 BEAM4 elements are selected as the spatial beam elements. The element of MPC 184 with the number of 122 is modeled as the multiple-point constraints element. In addition, the element of COMBIN 14 with the number of 61 is used when tangential slip is considered.

For the full interaction case, the slip parameters are taken as . Figure 6 shows the vertical deflection based on the FEM, the solution calculated by this paper, and Thevendran et al. experimental results, respectively. For this paper solution, the number of Fourier coefficients (in (25)) can be chosen as large as we can to achieve a required accuracy. As an illustration of convergence of the results, we show the results from this paper by taking the number as 1, 5, and 10 terms, respectively. It can be seen in Figure 6 that the vertical deflections results solved from this paper and from Thevendran et al. experiment are in reasonable agreement (the difference is within 12%), and the solved results from this paper and from FEM are in good agreement.

Figure 6: Vertical deflection of the beam.

Considering only vertical deflections are available from the experiments of Thevendran et al., the other results of the solution calculated by this paper and by the FEM are compared in Figures 7–12. The results are the torsional angle , the bending moment (), the axial force (), the bending moment , the shear force , and the torsion moment , respectively.

Figure 7: Torsional angle of the beam.

Figure 8: Bending moment of the beam.

Figure 9: Axial force of the beam.

Figure 10: Bending moment of the beam.

Figure 11: Shear force of the beam.

Figure 12: Torsion moment of the beam.

It can be seen from Figures 7–12 that the results solved by this paper can be quickly convergent, and one can get very good results by taking more than 5 terms. When the stiffness of the shear connectors is changed, the vertical deflection and torsional angle at midspan section are shown in Figures 13 and 14, respectively. In Figure 13, it can be seen that the vertical deflections converge to the minimum when the slip parameter is greater than . Moreover, the changes with slip parameter indicate the shear connecter’s stiffness in experimental model of Thevendran et al. with the range from to . As it can be seen from Figures 13 and 14, the performance of the curved composite I-beam becomes significantly rigid when the shear connecter’s stiffness is increased from to . The reduction of shear connectors’ stiffness between layers has a significant effect on the vertical deflection and torsional angle at midspan section.

Figure 13: Vertical deflection at midspan section when the stiffness of the shear connecter is changed.

Figure 14: Torsional angle at midspan when the stiffness of the shear connecter is changed.

The inaccuracy of the vertical deflection and torsional angle between analytical solutions in this paper and the FEM results in Figures 6–14 is small and reasonable. Thus, the model created in this paper can be applied sufficiently for practical purposes.

7. Conclusions

In this paper, an analytical solution has been developed and presented for the simply supported horizontally composite curved I-beam. It is primarily used to solve a static problem of a two-layered composite curved beam with flexible shear connection in the tangential direction. The trigonometric series are adopted in solution expression in terms of span coordinate. Both governing equations and boundary conditions are obtained by using the Vlasov curved beam theory and the principle of minimum potential energy. The accuracy of calculation results based on the theory presented in this paper is verified by comparing with the FEM and the Thevendran et al. experimental results. Therefore, the analytical expression and model established in this paper possess certain actual engineering significance.

Appendix

A. Component of and

Consider

B. Component of , , and

Consider

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

The research described in this paper was financially supported by the National Natural Science Foundation of China (51378236) and project supported by the transportation science and technology plan of Jilin province (2012-4-1-4).

References

A. Chakrabarti, A. H. Sheikh, M. Griffith, and D. J. Oehlers, “Analysis of composite beams with longitudinal and transverse partial interactions using higher order beam theory,” International Journal of Mechanical Sciences, vol. 59, no. 1, pp. 115–125, 2012.

E. Oñate, A. Eijo, and S. Oller, “Simple and accurate two-noded beam element for composite laminated beams using a refined zigzag theory,” Computer Methods in Applied Mechanics and Engineering, vol. 213-216, pp. 362–382, 2012.

N. M. Newmark, C. P. Seiss, and I. M. Veist, “Tests and analysis of composite beams with incomplete interactions,” Proceedings of the Society for Experimental Stress Analysis, vol. 4, no. 1, pp. 75–92, 1951.

J. R. Goodman and E. P. Popkov, “Layered beam systems with interlayer slip,” Journal of Structural Division, vol. 24, no. 11, pp. 2535–2547, 1968.

J. R. Goodman and E. P. Popov, “Layered wood systems with inter-layer slip,” Wood Science, vol. 1, no. 3, pp. 148–158, 1969.

U. A. Girhammar and V. K. A. Gopu, “Composite beam-columns with interlayer slip—exact analysis,” Journal of Structural Engineering, vol. 199, no. 4, pp. 1265–1282, 1993.

U. A. Girhammar and D. Pan, “Dynamic analysis of composite members with interlayer slip,” International Journal of Solids and Structures, vol. 30, no. 6, pp. 797–823, 1993.

Y. C. Wang, “Deflection of steel-concrete composite beams with partial shear interaction,” Journal of Structural Engineering, vol. 124, no. 10, pp. 1159–1165, 1998.

A. Dall'Asta, “Composite beams with weak shear connection,” International Journal of Solids and Structures, vol. 38, no. 32-33, pp. 5605–5624, 2001.

J. Nie and C. S. Cai, “Steel–concrete composite beams considering shear slip effects,” Journal of Structural Engineering, vol. 129, no. 4, pp. 495–506, 2003.

G. Ranzi and M. A. Bradford, “Analytical solutions for the time-dependent behaviour of composite beams with partial interaction,” International Journal of Solids and Structures, vol. 43, no. 13, pp. 3770–3793, 2006.

H.-B. Liu, W.-H. Liu, and Y.-L. Zhang, “Calculation analysis of shearing slip for steel-concrete composite beam under concentrated load,” Applied Mathematics and Mechanics: English Edition, vol. 26, no. 6, pp. 735–740, 2005.

F. Campi and I. Monetto, “Analytical solutions of two-layer beams with interlayer slip and bi-linear interface law,” International Journal of Solids and Structures, vol. 50, no. 5, pp. 687–698, 2013.

A. Ayoub, “A two-field mixed variational principle for partially connected composite beams,” Finite Elements in Analysis and Design, vol. 37, no. 11, pp. 929–959, 2001.

G. Ranzi, M. A. Bradford, and B. Uy, “A direct stiffness analysis of a composite beam with partial interaction,” International Journal for Numerical Methods in Engineering, vol. 61, no. 5, pp. 657–672, 2004.

B. Čas, M. Saje, and I. Planinc, “Non-linear finite element analysis of composite planar frames with an interlayer slip,” Computers and Structures, vol. 82, no. 23–26, pp. 1901–1912, 2004.

F. Gara, G. Ranzi, and G. Leoni, “Displacement-based formulations for composite beams with longitudinal slip and vertical uplift,” International Journal for Numerical Methods in Engineering, vol. 65, no. 8, pp. 1197–1220, 2006.

W. Chung and E. D. Sotelino, “Three-dimensional finite element modeling of composite girder bridges,” Engineering Structures, vol. 28, no. 1, pp. 63–71, 2006.

Y. Majdi, C.-T. T. Hsu, and M. Zarei, “Finite element analysis of new composite floors having cold-formed steel and concrete slab,” Engineering Structures, vol. 77, no. 15, pp. 65–83, 2014.

V. Z. Vlasov, Thin Walled Elastic Beams, Israel Program for Scientific Translation, Jerusalem, Israel, 1961.

M.-Y. Kim, S.-B. Kim, and N.-I. Kim, “Spatial stability of shear deformable curved beams with non-symmetric thin-walled sections. I: stability formulation and closed-form solutions,” Computers and Structures, vol. 83, no. 31-32, pp. 2525–2541, 2005.

M. H. Kim, S. B. Kim, and N. I. Kim, “Spatial stability of shear deformable curved beams with non-symmetric thin walled sections, II: F.E. solutions and parametric study,” Computers & Structures, vol. 83, no. 31-32, pp. 2542–2558, 2005.

A. M. Yu, X. G. Yang, and G. H. Nie, “Generalized coordinate for warping of naturally curved and twisted beams with general cross-sectional shapes,” International Journal of Solids and Structures, vol. 43, no. 10, pp. 2853–2867, 2006.

F. N. Gimena, P. Gonzaga, and L. Gimena, “3D-curved beam element with varying cross-sectional area under generalized loads,” Engineering Structures, vol. 30, no. 2, pp. 404–411, 2008.

F. N. Gimena, P. Gonzaga, and L. Gimena, “Stiffness and transfer matrices of a non-naturally curved 3D-beam element,” Engineering Structures, vol. 30, no. 6, pp. 1770–1781, 2008.

L. Gimena, F. N. Gimena, and P. Gonzaga, “Structural analysis of a curved beam element defined in global coordinates,” Engineering Structures, vol. 30, no. 11, pp. 3355–3364, 2008.

V. Thevendran, N. E. Shanmugam, S. Chen, and J. Y. R. Liew, “Experimental study on steel-concrete composite beams curved in plan,” Engineering Structures, vol. 22, no. 8, pp. 877–889, 2000.

N. E. Shanmugam, M. Mahendrakumar, and V. Thevendran, “Ultimate load behaviour of horizontally curved plate girders,” Journal of Constructional Steel Research, vol. 59, no. 4, pp. 509–529, 2003.

C. Topkaya, E. B. Williamson, and K. H. Frank, “Behavior of curved steel trapezoidal box-girders during construction,” Engineering Structures, vol. 26, no. 6, pp. 721–733, 2004.

F. Giussani and F. Mola, “Service-stage analysis of curved composite steel-concrete bridge beams,” Journal of Structural Engineering, vol. 132, no. 12, pp. 1928–1939, 2006.

R. E. Erkmen and M. A. Bradford, “Nonlinear elastic analysis of composite beams curved in-plan,” Engineering Structures, vol. 31, no. 7, pp. 1613–1624, 2009.