Abstract
Deformation of partially composite beams under distributed loading and free vibrations of partially composite beams under various boundary conditions are examined in this paper. The weak-form quadrature element method, which is characterized by direct evaluation of the integrals involved in the variational description of a problem, is used. One quadrature element is normally sufficient for a partially composite beam regardless of the magnitude of the shear connection stiffness. The number of integration points in a quadrature element is adjustable in accordance with convergence requirement. Results are compared with those of various finite element formulations. It is shown that the weak form quadrature element solution for partially composite beams is free of slip locking, and high computational accuracy is achieved with smaller number of degrees of freedom. Besides, it is found that longitudinal inertia of motion cannot be simply neglected in assessment of dynamic behavior of partially composite beams.
1. Introduction
Composite structures of different materials with interlayer slip have found a wide range of applications in structural engineering, such as in building, bridge and shelter construction. Common examples of such structures in civil engineering are steel-concrete composite beams (or decks) in buildings and bridges, wood-concrete floor or wall systems, and so forth. Shear connectors are usually used to connect different portions of these composite structures such as headed studs for steel-concrete composite beams. Since the shear connector is not completely rigid in practice, an interlayer slip, also known as partial interaction, often develops. In some cases, the interlayer slip significantly affects the mechanical behavior of composite systems. Thus, the structural model must account for the interlayer slip between the components, and this is very important in analysis and design of such partially composite structural systems.
The problem of partially composite structures has been studied for several decades. One of the first analytical theories for partial interaction of composite beams was developed by Newmark et al. [1] who focussed their attention on steel-concrete composite beams with flexible shear connection based on Euler-Bernoulli beam hypothesis. Newmark’s model was later extended by many researchers. Recently, analytical solutions for static and dynamic responses of partially composite beams (PCBs) were developed by Girhammar et al. [2, 3]. Analytical solutions for static and dynamic responses of partially composite beams based on Timoshenko beam theory were obtained by Xu and Wu [4]. However, these analytical solutions are restricted to beams with single span and simple loads. Nonlinear analysis and linear elastic analysis of the engineering problems under more general circumstances have to resort to numerical tools [5–18] among which the finite element method (FEM) appears to be widely used [5–7]. It is well known that low-order finite element models of composite beams experience the so-called “slip locking” phenomenon for large values of the shear connection stiffness [8–10]. Possible strategies in finite element analysis to reduce or even to eliminate slip locking are lowering the degree of interpolation functions for the slip or introducing elements with larger numbers of degrees of freedom. These strategies proceed unexceptionally at the expense of increasing the total degrees of freedom of the finite element model to achieve an acceptable accuracy. Furthermore, the demand of more degrees of freedom is exacerbated when nonlinear analysis of partially composite beams is performed.
Another issue that may be worth mentioning is the evaluation of dynamic behavior of partially composite beams. In the available information about dynamic behavior of partially composite beams, the longitudinal inertia of motion has been neglected a priori. In the opinion of the authors, this is inherited from dynamic analysis of uniform beams to acquire mathematical convenience. Otherwise, analytical solutions of dynamic equations of partially composite beams would be rather difficult to develop due to the coupling of flexural and axial motions.
In this paper, the weak-form quadrature element method is applied to static and free vibrational analyses of partially composite beams. The weak-form quadrature element method (QEM) [11–14] has been shown to be rather successful in improving computational efficiency. So far, the QEM has been applied to various structural problems [11–14]. Representative cases with various levels of shear connection stiffness are studied in this paper and the results are compared with those of finite element formulations [10]. The rapid convergence, high accuracy, and slip-locking-free characteristics of the weak-form quadrature element method are highlighted. In addition, free vibration analysis of partially composite beams is conducted with the longitudinal inertia of motion included. It is shown that the longitudinal inertia of motion cannot be simply neglected since low-order frequencies may be significantly affected.
2. Kinematic Model and Quadrature Element Formulation
In this section, quadrature element formulation will be derived based on Newmark kinematical model.
2.1. Newmark Kinematical Model
For simplicity, partially composite beams are modeled as two sub-Euler-Bernoulli beams linked by a deformable connection which allows for relative axial displacement between the two layers without transverse separation [1]. The assumption of small displacements, rotations, and strains is made for the two layers. In this paper, planar flexural and free vibrational analyses of partially composite beams are considered.
A typical composite beam with a symmetric cross-section is shown in Figure 1, where the top subcomponent is formed by a concrete slab while the bottom subcomponent is formed by a steel joist. Without losing generality, the formulation is derived for a beam segment of length and about an arbitrary reference axis located at distances and from the centroids of two subcomponents, respectively.
The displacement field of the partially composite beam consists of transverse displacement , axial displacements and at the centroids of the top and bottom sub-components. For simplicity, these displacements will be abbreviated as , , and . The expressions for the rotation and the curvature at a cross-section of the beam, denoted as , , respectively, can be obtained by differentiating the deflection with respect to coordinate .
Under the assumptions of Bernoulli-Euler beams, only axial strains exist, that is, while the slip at the interface is expressed by the following relation: where is the distance between the two subcomponent centroids (Figure 1). The parameters describing the strains can be collected in a vector: whose corresponding internal forces can be expressed in a vector as with the axial forces in two subcomponents , and the summation of the bending moments of two subcomponents . Then, the linear generalized constitutive relation between strains and internal forces can be written in the form: The elasticity matrix is where , , and are the cross-section area, the second moment of area about the neutral axis and the elastic modulus of the material of the top subcomponent; symbols with superscript designate relevant properties of the bottom subcomponent. The interface shear force can be written as where is the shear connection stiffness.
2.2. Quadrature Element Formulation
In quadrature element analysis, the following normalized dimensionless coordinate is introduced in a quadrature-beam element: The strain energy of the element is then written as where with The displacement vector for the interlayer slip is defined as with The interlayer elasticity matrix is given by It is noted that quantities with an overbar are designated as variables in the normalized coordinate system.
An efficient numerical scheme is chosen first to evaluate the strain energy. Thus, (2.9) can be rewritten as where denotes the weighting coefficients in the numerical integration scheme. Usually, Lobatto rule [19] is employed in a quadrature element whose sampling point distribution in the normalized domain is where is the zero of the first-order derivative of the -order Legendre polynomial .
Denote the slopes at the two ends of the beam as and and note that The nodal displacement vector of the quadrature element in the physical domain is written as (see Figure 2) which is related to its counterpart in the normalized dimensionless coordinate system through where With differential quadrature analogs [20], the normalized strain at a node (integration point) in (2.10) can be related to the normalized nodal displacement vector as Similarly, the normalized displacement vector for the interlayer slip can also be expressed in terms of the normalized nodal displacement vector as The essence of the conventional differential quadrature analog is that the derivative of a function at a node is expressed as weighted linear summation of function values at all nodes in the domain of interest [20].
Symbolically, the mth-order derivative of a function with respect to a variable at a node is approximated by where are the weighting coefficients for mth-order derivatives determined by using Lagrange interpolation functions as basis functions [21] and is function value at node . Details of differential quadrature analogs can be referred to [22, 23].
For Bernoulli-Euler beams, an element end condition necessitates the consideration of the slope as well as the deflection of the beam. Thus, a generalized differential quadrature analog is introduced as where are weighting coefficients for mth order derivatives in the generalized differential quadrature analog and determined by using Hermitian interpolation functions as basis functions [13]. All derivatives involved in (2.11) and (2.13) are approximated using the differential quadrature analogs given in (2.23) and (2.24), that is, Substitution of (2.25) into (2.21) and (2.22) yields the expressions of and , that is, where is the Kronecker delta, implying Then, (2.15) is rewritten as where
The kinetic energy of the partially composite beam consists of three portions: longitudinal kinetic energy, flexural energy, and rotational energy which are given, respectively, by where the single dot stands for the first-order partial derivative with respect to time variable. Introduction of the same numerical integration and the normalized nodal displacement vector turns the kinetic energy of the beam into where with The energy expression of external forces (see Figure 3) applied on the beam element is The same numerical integration scheme was introduced into (2.34), where is the element force vector given by with Hamilton principle dictates the following condition of the element: where and are two instants of time. Substitution of (2.28), (2.31), and (2.34) into (2.37) yields where the double dot represents the second-order differentiation with respect to time variable. Introducing the coordinate transformation (see (2.19)), the element equilibrium equations in physical domain become where Applying Hamilton principle to the entire beam yields the global dynamic equilibrium equations of the system where is the global nodal displacement vector; , , and are the global stiffness matrix, the global mass matrix, and the global force vector, respectively.
There are three typical boundary conditions that are usually considered in analysis of composite beams: clamped end (C), simply supported end (S), and free end (F). For a clamped end, the following conditions are enforced: For a simply supported end, it is required that The displacement variables or internal forces at an arbitrary position can be obtained using Langrange interpolation and differential quadrature analogs after solving for the nodal displacements in (2.41).
3. Numerical Results
To illustrate the effectiveness of the QEM in analysis of partially composite beams, examples are studied in this section.
3.1. Flexural Analysis
In static analysis, neglecting the inertial terms, (2.41) is reduced to For comparison, a typical partially composite beam reported in [9] is re-examined firstly. The partially composite beam of 20 meters in length () is subjected to a uniformly distributed load of 35 kN/m. The two subcomponents of the beam are a rectangular reinforced concrete slab (2300 mm × 230 mm) with 1% reinforcement and steel joist with an I-shaped cross-section. The steel joist has both flanges of 500 mm × 40 mm and a web of 1120 mm × 16 mm. The elastic modulus for the steel reinforcement and joist is 210,000 MPa while the elastic modulus for the concrete is 34,200 MPa. The densities of the two materials are 2400 Kg/m3 and 7850 Kg/m3, respectively. The beam with various levels of shear connection is considered to evaluate the proposed quadrature element formulation for varying from 1 to 40 where is the dimensionless shear stiffness coefficient first identified by Girhammar and Pan [2] and defined as The chosen values of in the paper are as typical as those used in many other references [8–10].
For brevity, two capital letters which stand for the boundary conditions at the two ends of the beam are used to designate the beam. Thus, four typical boundary conditions of the beam: SS, FC, SC, and CC are considered in the present investigation.
The relative errors of the maximum deflection and the maximum slip in the beam are defined as respectively. The exact values of the two displacements in (3.3) under four boundary conditions are computed using the exact solution given in [2].
The computed results of the maximum deflection for various shear connection stiffnesses are listed and compared with finite element solutions in Tables 1, 2, 3, and 4. Acronyms are used to refer to the results for various formulations and discretizations. For instance, 1Q8.26DOF is used to designate one quadrature element with 8 nodes and 26 total degrees of freedom; 4F8.20DOF stands for four finite elements, each with eight degrees of freedom, and 20 degrees of freedom in total. The nodes in an element are equally spaced for all the finite elements as proposed by Dall’Asta and Zona [10] and the beam is divided into elements of the same size.
From Tables 1–4, it is seen that both the quadrature element and finite element formulations are able to provide a good estimate of the maximum deflection under various levels of connection stiffness. It is found that quadrature element formulation offers high computational efficiency than finite element formulation does, especially for large shear connection stiffness, say, . For instance, the error of one four-node quadrature element with 14 DOF in total is less than that of eight 8DOF finite elements with 36DOF in total. It is observed that given the same number of degrees of freedom, the relative error of the maximum deflection is in the following ascending sequence for the four boundary conditions: SS, FC, SC and CC in example one. This is believed to be ascribed to the high-order feature due to the enhancement of constraints at the two ends of the composite beam.
The computed results of the maximum slip for various shear connection stiffnesses of the two examples are listed and compared with the finite element solutions in Tables 5, 6, 7, and 8. For computation of the slip, low-order finite element solutions have been afflicted with large error for partially composite beams with large shear connection stiffness regardless of the end conditions. This phenomenon has been known as “slip locking.” In order to overcome the problem, finite elements with more degrees of freedom were proposed [10] by introducing additional internal nodes into low-order finite elements. From Tables 5–8, it is observed that the accuracy of slip for beams with high shear connection stiffness using high-order finite elements indeed improves. In contrast, quadrature element solutions of the slip are less affected by the sliplocking especially for elements with large .
It is seen that the quadrature element solution of the relative error of the maximum slip fluctuates slightly with the increase of shear connection stiffness. Given the same number of DOF, large shear connection stiffness in general results in large maximum slip error and the largest relative error of the maximum slip occurs in the SC beam. For partially composite beams with large shear connection stiffness, say, , the number of integration points in a quadrature element has to exceed twelve to keep the relative error of the maximum slip below 1.0%.
It is noted that the change of shear connection stiffness has significant influence on the variation of slip. The variation of slip along the beam length for a partially composite SC beam is given in Figure 4 for and . In general, quadrature element solutions and higher-order finite element formulations provide acceptable slip predictions for different shear connection levels. Nonetheless, the quadrature element solutions exhibit higher accuracy than those of finite element solutions given the same total number of degrees of freedom. It is found that low-order finite element solutions suffer oscillations especially for beams with large shear connection stiffness which were also reported and explained in detail in [10].
(a) |
(b) |
To further validate the capacity of the QEM, a continuous two-span composite beam of 40 meters in length () with the same geometric and material data as those of example one is investigated. The external loads are two concentrated load of kN applied at the two midspans, respectively. The computed results of the maximum deflection and the maximum slip for various shear connection stiffnesses are listed and compared with finite element solutions in Tables 9 and 10. It is seen that similar observations as in example one are acquired and the high computational accuracy and efficiency of the QEM are verified.
3.2. Free Vibration Analysis
For free vibration analysis of beams, it is assumed that Substitution of (3.4) and into (2.41) yields which is a generalized eigenvalue problem. According to [24], a dimensionless frequency of the partially composite beam can be defined as For simplicity, a partially composite beam with the same geometric and material data as in Subsection 3.1 is adopted for free vibration analysis. First, an SS (both ends simply supported) beam is considered to validate the convergence of the quadrature element solution. The first and the third dimensionless natural frequencies for two levels of connection stiffness are given in Figures 5 and 6.
(a) |
(b) |
(a) |
(b) |
It is seen that the quadrature element solution converges more rapidly than 8DOF and 16DOF finite element solutions and they agree with each other with the increase of total degrees of freedom. The 16DOF finite element can be used to acquire results with accuracy comparable with that of quadrature element solutions. However, a number of finite elements is needed in contrast to only one element in quadrature element analysis of partially composite beams. The dashed lines in Figures 5 and 6 are the analytical solutions given in [3] neglecting longitudinal inertias of motion ( in (2.31)). Obviously, omission of the longitudinal inertias of motion of the composite beam overestimates natural frequencies. Although the fundamental frequency of the beam is overestimated by less than 1%, the third-order frequency is overestimated by 35% and 44% for small shear connection stiffness case () and large shear connection stiffness case (), respectively.
To further investigate the dynamic behavior of partially composite beams, the effects of rotary inertia ( in (2.31)) and longitudinal inertia ( in (2.31)) on the frequencies and mode shapes of the beam are assessed. A quadrature element with 15 nodes (47 total degrees of freedom) is adopted to ensure the accuracy of solution. Although an allowance for rotary inertia results in reduction of frequencies of the beam, the amount of reduction is small especially for slender beams like the one in question. It is observed that the reduction of the fifth frequency of the simply supported composite beam, the largest of the first five frequencies, is less than 1.5%, indicating the negligible effect of the rotary inertia for low-order frequencies of the beam. The longitudinal inertia, on the other hand, has to be included in dynamic analysis since it can lead to significant change of vibration mode shapes as well as frequency characteristics of the beam.
As shown in Figure 7, ignoring longitudinal inertia leads to overestimation of frequencies; the higher the frequency order, the larger frequency overestimation which increases with shear connection stiffness and gradually levels off. It is noted that the first two frequencies for SS and SC beams and the first three frequencies for FC and CC beams are virtually not affected by the longitudinal inertia of the beam.
(a) SS
(b) FC
(c) SC
(d) CC
In Figures 8, 9, 10, and 11, the first five longitudinal and flexural modes of the beam under four different boundary conditions are displayed for small shear connection stiffness and large shear connection stiffness , respectively. It is seen that in most cases longitudinal and flexural modes are coupled despite occasional predominance of either longitudinal or flexural mode for a certain-order frequency of partially composite beams. In addition, the longitudinal mode and the flexural mode for a certain-order frequency are more likely coupled with the increase of the shear connection stiffness or beam length. Since the longitudinal mode may appear at low-order frequencies of the composite beam (see Figure 10), it is concluded that longitudinal inertia of motion cannot be simply neglected in dynamic analysis of partially composite beams. The effect of cross-sectional dimensions of the joist on the first five frequencies and mode shapes are shown in Figure 12. It can be seen that longitudinal modes may appear in low-order frequencies when the flexural rigidity of the joist is relatively large.
(a) Longitudinal mode shape |
(b) Longitudinal mode shape |
(c) Flexural mode shape |
(d) Flexural mode shape |
(a) Longitudinal mode shape |
(b) Longitudinal mode shape |
(c) Flexural mode shape |
(d) Flexural mode shape |
(a) Longitudinal mode shape |
(b) Longitudinal mode shape |
(c) Flexural mode shape |
(d) Flexural mode shape |
(a) Longitudinal mode shape |
(b) Longitudinal mode shape |
(c) Flexural mode shape |
(d) Flexural mode shape |
(a) Longitudinal mode shape |
(b) Longitudinal mode shape |
(c) Flexural mode shape |
(d) Flexural mode shape |
4. Concluding Remarks
Static and free vibrational analyses of partially composite beams have been studied using the weak-form quadrature element method. It has been shown that the weak-form quadrature element method works very well for partially composite beams with any level of connection stiffness. One quadrature element is usually sufficient for a partially composite beam under distributed load and the number of nodes is adjustable in accordance with the convergence requirement. However, division is clearly needed for quadrature element analysis of a multi-span beam as well as a single-span beam under complex loadings like stepwise distributed loads. It is found that longitudinal inertia of motion cannot be simply neglected in dynamic analysis of partially composite beams since low-order frequencies as well as the mode shapes may be significantly affected. Although the QEM applications presented in this paper are limited to static and free vibration analysis, the extension of the present work to further dynamic analysis, geometric, and material nonlinear analysis of partially composite beams is straightforward and will be reported in an upcoming paper.
Acknowledgment
This study was undertaken under the support of the National Natural Science Foundation of China (no. 51178247).