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

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.


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][15][16][17][18][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][22][23][24][25][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 inplane 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.

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 Igirder 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 [N/m 2 ].
(3) The uplift between the slab and I-girder is neglected.The radius of curvature is constant along the beam.
The deflections in tangential direction, vertical direction, and radial direction (-direction, -direction, and direction) are indicated as , , and V, 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 (, ,  1 , V 1 ) and (, ,  2 , V 2 ), 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  0 and torsion moment  0 .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 By eliminating  1 and  2 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

Equilibrium at the Interface
At the interface, considering the deflections described in Figure 3, shear force  2 described in Figure 2 can be written as where  =  1 +  2 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 Considering deflection described in Figure 4, equilibrium at the interface in the radial direction can be written as 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 V 2 and  2 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 V 2 and  2 by inserting these two derivative equations into the second formulas of (1), using (2), one can obtain the following: where

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  1 and  2 . where is the torsional strain energy due to the internal torsion moments  1 and  2 . where is the axial strain energy due to the internal axial forces  1 and  2 .
is the bending strain energy due to the internal bending moments  1 and  2 .
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), , , V 1 ,  1 ,   , and   are all independent variables.The variation of Π * is where Γ  ( = 1∼6) and   ( = 1∼10) are calculated from the variation of Π * , which are given in Appendix A. Each term in (21) must be identically zero.It makes Γ 5 = 0, Γ 6 = 0,  5 = 0,  6 = 0,  7 = 0,  8 = 0, and  9 = 0; that is, the undermined Lagrange multipliers   and   are Equation ( 22) notes that the Lagrange multipliers and the deflections ( 1 and V 1 ) 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 Γ 1 = 0, Γ 2 = 0, Γ 3 = 0, and Γ 4 = 0, 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.

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  ( = 1, . . ., ) and  is span length.The applied distributed load  0 and  0 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 (0, ), 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   ( = 1, . . ., 5) are the relevant coefficients which are given in Appendix B.

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  = 150 KN vertical load at the midspan section.
The material properties and dimensions of the curved I-beam (SP4 beam) are shown in Table 2.
The finite element software ANSYS is used to model the described structure (Figure 5).In the FEM model, both the

Type
The boundary conditions R Free end

Space beam element Rigid link
Interface node and shear connector 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.The vertical deflection , the tangential deflection , the radial deflection V, 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  = 10 10 N/m 2 .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.
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  1 ( 1 =   ), the axial force  1 ( 1 =   ), the bending moment   , the shear force   , and the torsion moment   , respectively.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 10 10 N/m 2 .Moreover, the   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 10 7 to 10 9 N/m 2 .The reduction of shear connectors' stiffness between layers has a significant effect on the vertical deflection and torsional angle at midspan section.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.

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 Consider

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

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

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

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

Figure 10 :
Figure 10: Bending moment   of the beam.

Figure 11 :
Figure 11: Shear force   of the beam.

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

Table 1 :
Common support conditions.

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