Research Article  Open Access
Improved Finite Beam Element Method to Analyze the Natural Vibration of SteelConcrete Composite Truss Beam
Abstract
Based on Hamilton’s principle, this study has developed a continuous treatment for the steelconcrete composite truss beam (SCCTB). It has also deduced the SCCTB element stiffness matrix and mass matrix, which include the effects of interface slip, shear deformation, moment of inertia, and many other influencing factors. A finite beam element method (FBEM) program for SCCTB’s natural vibration frequency has been developed and used to calculate the natural vibration frequencies of several SCCTBs with different spans and different degrees of shear connections. The FBEM’s calculation results of several SCCTBs agree well with the results obtained from ANSYS. Based on the results of this study, the following conclusions can be drawn. For the SCCTB with highorder natural vibration frequency and with short span, the effect of the shear deformation is greater. Hence, the effect of the shear deformation on the SCCTB’s highorder natural vibration frequency cannot be ignored. On the other hand, the effect of the interface slip on the SCCTB’s highorder natural vibration frequency is insignificant. However, the effect of the interface slip on the SCCTB’s loworder natural vibration frequency cannot be ignored.
1. Introduction
The SCCTB is a newtype structural member. It is developed on the basis of ordinary composite beams and uses shear connectors to connect the concrete slab to the steel truss into one entire joint work. Comparing with ordinary truss beams, SCCTB is more effective because it uses both the compressive strength of concrete and the tensile strength of steel. Hence, it is characterized by high utilization rate of indoor space, high economic efficiency, high bending stiffness, and high bearing capacity [1, 2].
Even though a beam is designed as a full composite section, due to the deformation of shear studs, it cannot respond with perfect composite action without slip. Further, unlike Euler–Bernoulli beams in which there is infinite shear stiffness, it may not be appropriate to design other beams by neglecting the shear deformation of the concrete slab and steel truss beam, for example, composite beams with a small spantodepth ratio [3–6]. Therefore, the natural vibration characteristics of SCCTBs are affected by the coupled effects of the shear deformation and the interfacial slip between the concrete slab and steel truss beam.
Hitherto, there are many studies on the shear behavior of hybrid steeltrussedconcrete beams, which are constructed by embedding the prefabricated steel trusses into a concrete core cast in situ. On the other hand, there are few studies on the mechanical properties of SCCTB. Using the software ABAQUS to carry out the finite element (FE) numerical simulations, Monaco [7] investigated the shear behavior of hybrid steeltrussedconcrete beams. Monti and Petrone [8] developed the shear capacity equations from a purposely developed mechanicsbased shear model for hybrid steeltrussedconcrete beams. Campione et al. [9] investigated a calculation method for the prediction of the shear resistance of precast composite beam; both experimental and numerical results were used to validate the developed analytical expressions.
Giltner and Kassimali [10] developed a method which can replace the trusses and the beam elements with an equivalent beam, thereby reducing the size of the computer model for analysis. Machacek and Cudejko [11, 12] investigated the distribution of longitudinal shear along the interface between the steel and the concrete of various composite truss bridges from the elastic phase to the plastic collapse. He found that the nonlinear distribution of the longitudinal shear significantly depends on the rigidity of the shear connection and the densification of the shear connectors above the truss nodes. By carrying out an experimental study and a theoretical analysis, Chan and Fong [13] concluded that the use of the effective length method in the linear analysis and design method was less convenient and less accurate than the secondorder analysis. On the basis of the literature [13], Fong et al. [14] further showed that the secondorder analysis method was not only an accurate design method, but it can also avoid the uncertain approximate value of the effective length. By comparing the results from the finite element numerical calculation method with the experimental results, Bujnak and Bouchair [15] found that the local effects of the concentrated longitudinal shear forces should be examined in SCCTB. Siekierski [16] analyzed the effects of the shrinkage of concrete slab in SCCTB and developed a set of linear equations to compute the axial forces in the flange members of the truss girder and the transverse shear forces in SCCTB.
The natural vibration characteristic analysis is the basis to investigate the dynamic characteristics of SCCTB. So, it is necessary to investigate its natural vibration characteristics. As shown by the earlier studies, the mechanical properties of SCCTB are affected by the interface slip, shear deformation, moment of inertia, and many other factors. However, there are only a few studies on the SCCTB’s natural vibration characteristics which have taken these factors into account. In this study, a continuous treatment for SCCTB has been developed. It has also deduced the SCCTB element stiffness matrix and mass matrix with the cubic Hermite polynomial shape function, which includes the effects of the interface slip, shear deformation, moment of inertia, and many other influencing factors. Hence, even with fewer degrees of freedom, the precision in the calculation is satisfactory. Based on the developed SCCTB element, this study has developed a FBEM program which can calculate the natural vibration frequency of SCCTBs commonly used in the engineering practice. Finally, it calculates the natural vibration frequencies of several SCCTBs with different spans and different degrees of shear connections. Based on the results of the analyses, some conclusions which are related to the engineering design are drawn.
2. CrossSectional Analysis of Strain and Stress of SCCTB
2.1. Displacement Model of SCCTB
In Figure 1, , , , and are the height and width of the external and internal walls of the upper chord, respectively; , , , and are the height and width of the external and internal walls of the lower chord, respectively; , , , and are the height and width of external and internal walls of the oblique web member, respectively; , , , and are the height and width of the external and internal walls of the vertical web member, respectively; , , , and are the height and width of the external and internal walls of the lower horizontal connection member, respectively; , , , and are the height and width of the external and internal walls of the oblique bracing member, respectively; and are the thicknesses of the concrete roof and the cantilever slabs, respectively; and are the widths of the concrete roof and the cantilever slabs, respectively; and , , , , , and are the lengths of the upper chord, lower chord, oblique web member, vertical web member, lower horizontal connection member, and oblique bracing member, respectively.
The strain of oblique bracing member can be expressed aswhere is the angle between lower horizontal connection member and oblique bracing member, is the axial displacement of the oblique bracing member, and is the axial displacement of the lower chord.
Further, the axial force of the oblique bracing member is given by
The longitudinal component of the axial force of the oblique bracing member is given by
Therefore, the longitudinal equivalent area of the oblique bracing member can be expressed as follows:
The longitudinal displacement of SCCTB can be expressed as follows [17]:where is the coordinate of the centroid of the steel truss beam; is the coordinate of the centroid of the concrete slab; is the crosssectional rotation angle; is the vertical deflection of the SCCTB; is the difference between the longitudinal displacement of the concrete slab centroid and that of the steel truss beam centroid; is the elastic modulus of the steel truss and is the elastic modulus of the concrete slab. Further, ; ; ; ; ; ; ; ; ; and .
The axial displacement of the oblique web member is given bywhere is the angle between the oblique and vertical web members.
2.2. Strain Model of SCCTB
According to the displacement model of SCCTB in Section 2.1, the longitudinal strain of the crosssection can be expressed aswhere are the longitudinal strains of the top slab, cantilever slab, upper chord, and lower chord, respectively; is the normal strain of the oblique web member.
Derived from (5), the longitudinal relative slip between the concrete slab and the truss iswhere ; and are the distances from the centroids of the concrete slab and steel truss to the interface, respectively.
2.3. Stress Model of SCCTB
Based on the strain models of SCCTB in Section 2.2, the stresses of SCCTB arewhere are the longitudinal stresses of the top slab, cantilever slab, upper chord, and lower chord, respectively; is the shear stress of the upper chord and the lower chord, respectively; and is the normal stress of the oblique web member.
The shearing force per unit length of the interface is given bywhere is the longitudinal shear stiffness in per unit length between the concrete slab and the steel truss beam and the unit is ; is the longitudinal shear stiffness of the shear connections, and the unit is ; is the number of shear connections in a crosswise row; is the longitudinal spacing between two studs; is the effective span of SCCTB; is the degree of shear connection; is the yield strength of steel material; and is the shear strength of the stud.
3. Finite Beam Element Method of SCCTB
3.1. Improved SCCTB Element Stiffness Matrix and Mass Matrix Including the Effects of Interface Slip and Shear Deformation
In order to satisfy the consistency requirements, for the different displacement functions to contribute to the same strain function, the degree of polynomial has to be the same after the finite element approximation. Otherwise, there are unreal geometric constraint conditions, which can significantly reduce the accuracy of the finite element method and result in shear locking [5, 18]. The highestorder derivatives of each displacement function in (8)–(11) show that the highest derivatives of and are firstorder derivatives. Hence, the corresponding shape function can only satisfy the continuity condition for . Further, (9) shows that the shape function of the deflection function must satisfy the continuity condition of . According to the principle of numerical calculation, in order to satisfy the continuity condition of , the Hermite polynomial has to be at least cubic. However, (9) shows that the contribution of the deflection function to the crosssection rotation function is a quadratic polynomial. Hence, in order to satisfy the consistency requirements, the displacement functions and have to be approximated using the quadratic polynomials. Further, the elements satisfying the consistency requirements have to have at least 10 degrees of freedom, as shown in Figure 2. Assume that the nodal displacement vector is
The Hermite polynomial shape function to satisfy the continuity condition of at the element boundary iswhere and is the element length.
Substituting (15) into (8) gives
Substituting (15) into (12)(13) gives
The kinetic energy of the SCCTB can be expressed aswhere ; and are the densities of the steel material and the concrete material, respectively; and , .
The strain energy of SCCTB can be expressed aswhere and .
According to (22) and (23), the variation forms of the strain energy and the kinetic energy can be expressed as
By substituting (18)–(20) into (24), the SCCTB element mass matrix iswhere ; ; ; ; ; ; ; ; ; ; ; ; .
By substituting (18)–(20) into (25), the SCCTB element stiffness matrix is
3.2. Solving the Natural Vibration Frequency of SCCTB
Based on the SCCTB element stiffness matrix , element mass matrix , and element displacement vector , and using the “seat by number” method, the SCCTB overall stiffness matrix , overall mass matrix , and overall degree of freedom vector can be obtained. The common boundary conditions for the SCCTB can be expressed as follows [17, 19]:
The transformational relationship between the SCCTB overall degree of freedom vector before and after applying the boundary constraint can be expressed as follows:where is the transformation matrix for the overall degree of freedom vector, which can be obtained for the given boundary conditions, and is the SCCTB overall degree of freedom vector after applying the boundary constraint.
The free vibration function for a structure with multiple degrees of freedom can be expressed as follows:
Letwhere is the structure’s natural vibration frequency.
Substituting (31) into (30) gives
Substituting (29) into (32) giveswhere
Using (33), the frequency equation of SCCTB can be obtained as follows:
The software MATLAB has been used to develop the FBEM program of the abovementioned SCCTB element. Then, the natural vibration frequency of SCCTB has been solved by including the effects of both the interface slip and shear deformation.
3.3. Simplification of SCCTB Element Stiffness Matrix
If the shear deformation is not considered, then
Using (36), the following expression can be obtained:
Using (37) and excluding the effect of the shear deformation, the element stiffness matrix and mass matrix arewhere is a transformation matrix and can be obtained using (37).
After obtaining the SCCTB element stiffness matrix and mass matrix without the effect of the shear deformation, the same method in Section 3.2 can be used to calculate the natural vibration frequency of SCCTB without the effect of the shear deformation.
4. Cases for Analyses
To verify the accuracy of the developed FBEM (Section 3), the natural frequencies of two groups of clamped supported SCCTBs (i.e., SCCTB1 and SCCTB2) have been calculated using both the finite element method and the FBEM. For each group, there are cases with five degree of shear connection () and two spans. The mechanical and geometrical parameters of the SCCTBs are as follows:
are the calculation lengths of the two groups of SCCTBs, respectively. Further, and are Poisson’s ratio of steel and concrete, respectively.
The calculations by the finite element method have been carried out using the finite element program, ANSYS. The upper chord and the lower chord have been simulated using SHELL43 shell element. The vertical web member, oblique web member, lower horizontal connection member, and oblique bracing member have been simulated using BEAM188 elements. The concrete slab has been simulated using SOLID65 solid elements. The studs have been simulated using COMBIN14 spring elements. The elastic modulus of the spring element has been calculated using (14). In the finite element models, by coupling the degrees of freedom in the vertical direction of the nodes at the same position, the interface between the concrete slab and the steel truss beam is connected in the vertical; that is, there is no vertical separation between the concrete slab and the steel truss beam. In order to simulate the clamped supported at the ends of the beams in the finite element models, there are constrains in the degrees of freedom in the vertical, transverse, and longitudinal directions.
Tables 1 and 2 show a comparison of the results of the calculated natural frequencies by the FBEM and ANSYS model. is the calculation result from the ANSYS model. is the calculation result from the FBEM including the effects of the shear deformation and the interface slip. is the calculation result without including the effect of the shear deformation. Further, is the calculation error in FBEM, and and , where and are the errors due to the effects of the shear deformation and the interface slip, respectively.


From Tables 1 and 2 and Figures 3–5, the following can be seen:
() By including the effects of the interface slip, shear deformation, moment of inertia, and many other influencing factors on the SCCTB, the results of the FBEM’s calculations agree with those of the ANSYS’s finite element calculations. The maximum calculation error of the first six orders of natural vibration frequency is less than 4.6%. This is an indication that the SCCTB element stiffness matrix and mass matrix developed in this study are rational and effective. According to the calculation results of ANSYS, the local vibration of SCCTB1 with small span is becoming obvious in the 6th mode, while SCCTB2 with larger span does not cause this kind of vibration in high modes. Neglecting local vibration will result in generating additional restraints on the SCCTB, which will overestimate the restrain rigidity of the SCCTB. The proposed FBEM model in this paper can be used to describe the actual behavior of SCCTBs in flexural vibrations effectively, but it does not take into account the local vibration. Therefore, the 6th natural frequency of SCCTB1 using FBEM is slightly larger than the ANSYS’ results.
() Without including the effect of the shear deformation on the SCCTB’s natural vibration frequency, the FBEM’s calculation results are greater than those of the ANSYS’s calculations. Further, SCCTB’s shear deformation effect is greater for SCCTB with the highorder natural vibration frequency. If SCCTB’s natural vibration frequency order is six, the shear deformation effect is 36.5%. Hence, the shear deformation effect on the SCCTB’s highorder natural vibration frequency cannot be ignored.
() For SCCTB with a shorter span, the shear deformation effect is greater.
() The section rotation caused by the bending moment and the shear deformation caused by the shear force both result in a bending displacement of the structural member. The shear deformation effect on SCCTB at the loworder frequency is insignificant. This is an indication that the SCCTB’s bending vibrationtype of loworder is mainly due to the bending deformation caused by the section rotation and the bending deformation caused by shear deformation is small.
() For the shear deformation effectnatural vibration frequency order curves of SCCTB under different degrees of shear connections, they overlap with each other. This is an indication that the effect of the degree of shear connections on SCCTB’s shear deformation effect is not significant. It is because the shear force of SCCTB is mainly carried by the web member of the steel truss and has little relation to the degree of shear connections.
() While the interface slip effect on the SCCTB’s highorder natural vibration frequency is insignificant, the effect on the SCCTB’s loworder natural vibration frequency can be as high as 10.1%. Hence, the effect of interface slip stiffness on the SCCTB’s loworder natural vibration frequency cannot be ignored. According to the results of the analyses, this is because SCCTB’s section bending stiffness increases with increasing interface slip stiffness, while SCCTB’s bending vibrationtype of loworder is dominated by the bending deformation caused by the section rotation.
Figures 6 and 7 show the comparison between the first six flexural mode shapes of ANSYS and FBEM models of SCCTB1 and SCCTB2. It can be seen that the first six flexural mode shapes of FBEM’s calculations agree well with those of the ANSYS’s finite element calculations when including the effects of the interface slip, shear deformation, moment of inertia, and many other influencing factors on the SCCTB, which shows that the proposed FBEM model can be used to describe the actual behavior of SCCTBs in flexural vibrations effectively and accurately.
5. Conclusions
By including the effects of the interface slip, shear deformation, moment of inertia, and many other influencing factors on SCCTB, and based on Hamilton’s principle, the SCCTB element stiffness matrix and mass matrix with cubic Hermite polynomial shape function have been developed. Based on the developed SCCTB element, this study has also developed a FBEM program which can calculate the natural vibration frequency of SCCTBs that are common in the engineering practice. The program has been used to calculate the natural vibration frequencies of several SCCTBs with different spans and different degrees of shear connections. The conclusions are as follows:(1)The FBEM’s calculation results of the natural vibration frequencies of several SCCTBs agree well with the results of ANSYS’s finite element calculations. This is an indication that the SCCTB element stiffness matrix and mass matrix developed in this study are effective. Hence, this is a basis for the further applications of the FBEM to the dynamic calculations of SCCTB.(2)The SCCTB’s bending vibrationtype of loworder is dominated by the bending deformation caused by the section rotation, and the shear deformation effect of SCCTB’s loworder natural vibration frequency is insignificant.(3)For the SCCTB with highorder natural vibration frequency, the effect of the shear deformation is greater. Hence, the shear deformation effect on the SCCTB’s highorder natural vibration frequency cannot be ignored.(4)Under different degrees of shear connections, the curves for the shear deformation effect and natural vibration frequency order of SCCTBs are overlapped with each other. This is an indication that the effect of the degree of shear connections on SCCTB’s shear deformation is not significant.(5)While the interface slip effect on the SCCTB’s highorder natural vibration frequency is insignificant, the effect on the SCCTB’s loworder natural vibration frequency is dominated by the section rotation deformation which cannot be ignored.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The research described in this paper was financially supported by the National Natural Science Foundation of China (51778630, 51408449 and 51378502), the Special Fund of Strategic Leader in Central South University of China (under Grant 2016CSU001), and the Fundamental Research Funds for the Central Universities of Central South University of China (2016zzts078).
References
 F.X. Ding, J. Liu, X.M. Liu, F.Q. Guo, and L.Z. Jiang, “Flexural stiffness of steelconcrete composite beam under positive moment,” Steel and Composite Structures, vol. 20, no. 6, pp. 1369–1389, 2016. View at: Publisher Site  Google Scholar
 J. Liu, F.X. Ding, X.M. Liu, and Z.W. Yu, “Study on flexural capacity of simply supported steelconcrete composite beam,” Steel and Composite Structures, vol. 21, no. 4, pp. 829–847, 2016. View at: Publisher Site  Google Scholar
 J. G. Nie, C. S. Cai, T. R. Zhou, and Y. Li, “Experimental and analytical study of prestressed steelconcrete composite beams considering slip effect,” Journal of Structural Engineering, vol. 133, no. 4, pp. 530–540, 2007. View at: Publisher Site  Google Scholar
 G. Ranzi and A. Zona, “A steelconcrete composite beam model with partial interaction including the shear deformability of the steel component,” Engineering Structures, vol. 29, no. 11, pp. 3026–3041, 2007. View at: Publisher Site  Google Scholar
 A. Chakrabarti, A. H. Sheikh, M. Griffith, and D. J. Oehlers, “Dynamic response of composite beams with partial shear interaction using a higherorder beam theory,” Journal of Structural Engineering, vol. 139, no. 1, pp. 47–56, 2013. View at: Publisher Site  Google Scholar
 W. Zhou and W. Yan, “Refined nonlinear finite element modelling towards ultimate bending moment calculation for concrete composite beams under negative moment,” ThinWalled Structures, vol. 116, pp. 201–211, 2017. View at: Publisher Site  Google Scholar
 A. Monaco, “Numerical prediction of the shear response of semiprefabricated steelconcrete trussed beams,” Construction and Building Materials, vol. 124, pp. 462–474, 2016. View at: Publisher Site  Google Scholar
 G. Monti and F. Petrone, “Shear resisting mechanisms and capacity equations for composite truss beams,” Journal of Structural Engineering, vol. 141, 2015, UNSP 0401505212. View at: Publisher Site  Google Scholar
 G. Campione, P. Colajanni, and A. Monaco, “Analytical evaluation of steel–concrete composite trussed beam shear capacity,” Materials and Structures/Materiaux et Constructions, vol. 49, no. 8, pp. 3159–3176, 2016. View at: Publisher Site  Google Scholar
 B. Giltner and A. Kassimali, “Equivalent beam method for trusses,” Practice Periodical on Structural Design and Construction, vol. 5, no. 2, pp. 70–77, 2000. View at: Publisher Site  Google Scholar
 J. Machacek and M. Cudejko, “Longitudinal shear in composite steel and concrete trusses,” Engineering Structures, vol. 31, no. 6, pp. 1313–1320, 2009. View at: Publisher Site  Google Scholar
 J. Machacek and M. Cudejko, “Composite steel and concrete bridge trusses,” Engineering Structures, vol. 33, no. 12, pp. 3136–3142, 2011. View at: Publisher Site  Google Scholar
 S. L. Chan and M. Fong, “Experimental and analytical investigations of steel and composite trusses,” 2011. View at: Google Scholar
 M. Fong, S. L. Chan, and B. Uy, “Advanced design for trusses of steel and concretefilled tubular sections,” Engineering Structures, vol. 33, no. 12, pp. 3162–3171, 2011. View at: Publisher Site  Google Scholar
 J. Bujnak and A. Bouchair, “Theoretical and Experimental research on SteelConcrete Composite Truss,” 2014. View at: Google Scholar
 W. Siekierski, “Analysis of concrete shrinkage along truss bridge with steelconcrete composite deck,” Steel and Composite Structures, vol. 20, no. 6, pp. 1237–1257, 2016. View at: Publisher Site  Google Scholar
 W. Zhou, L. Jiang, Z. Huang, and S. Li, “Flexural natural vibration characteristics of composite beam considering shear deformation and interface slip,” Steel and Composite Structures, vol. 20, no. 5, pp. 1023–1042, 2016. View at: Publisher Site  Google Scholar
 F.F. Sun and O. S. Bursi, “Displacementbased and twofield mixed variational formulations for composite beams with shear lag,” Journal of Engineering Mechanics, vol. 131, no. 2, pp. 199–210, 2005. View at: Publisher Site  Google Scholar
 Z. Wangbao, L. ShuJin, J. Lizhong, and Q. Shiqiang, “Vibration analysis of steelconcrete composite box beams considering shear lag and slip,” Mathematical Problems in Engineering, vol. 2015, Article ID 601757, 2015. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2017 Zhipeng Lai 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.