DIFFERENTIAL QUADRATURE METHOD FOR VIBRATION ANALYSIS OF PRESTRESSED BEAMS

. This paper investigates natural frequencies of free transverse vibrations of prestressed beams, and the governing equations and natural frequencies of the free vibration in related literatures are discussed and corrected. The differential quadrature methods (DQ) are applied directly to the corrected governing equations to get the the values of natural frequency numerically. Under the simple supported boundary conditions, the natural frequencies of model beam are numerically studied, and the physical parameters of the beam are analyzed respectively. The numerical results show that the natural frequency values increase with the growth of concrete strength and eccentricity of prestressed steels. But with the increase of the span length of beam and values of original prestressing force, the natural frequency values decrease.


Introduction
Prestressed concrete technique can effectively improve the cracking of concrete structure, improve the stiffness and durability of the structure, and give full play to the high strength performance of materials, which has been widely used in civil engineering [1][2][3] . In recent years, many researchers have studied the static performance of prestressed concrete structure, but the investigations on the complex dynamic performance are relatively insufficient.
The differential quadrature (DQ) method is a powerful tool for dealing with dynamical problems. Firstly it was introduced for structural analysis by Bert et al. [4] . and since then it has been successfully employed for the analysis of vibration of beams [5][6] .
For the vibration characteristics of prestressed beams, the formulas of different linear vibration frequencies are derived by applying the theory of material mechanics and structural dynamics [3,7] , However, these formulas are based on simplified vibration control equations, and their accuracy is worth discussing. In this paper, the existing vibration control equations of prestressed beams are analyzed and corrected, and the differential quadrature method (DQ) is directly applied to the corrected equations. The differential equations are transformed to a generalize d eigenvalue problem, and the numerical solutions of natural frequency are obtained based on the method of DQ. Under the simply supported boundary conditions, the first four order numerical solutions of natural frequency of model beam are given, and the different effects of physical parameters of prestressed beam and prestressed steel on natural frequency are analyzed.
The present paper is organized as follows. Section 2 establishes the governing equation for the transverse free vibration of a finite elastic prestressed beam. Section 3 develops the DQ schemes to descretize the governing equation. Section 4 presents numerical examples for investigations on the effects of the different parameters of the relevant beam and prestressed steel on the natural frequencies. Section 5 ends the paper with the concluding remarks. Free transverse vibration of a prestressed beam can be described by transverse deflection of the cross section due to bending moment, depending on axial coordinate x and time t. Fig.1 shows a finite elastic prestressed concrete beam with a length l, the bending moment and the axial force can be expressed as

Equation of motion
The free vibration equation of a prestressed beam is governed by the equation [ Based on the assumption that increased prestressing forces P N  is proportional to midspan deflection y, the horizontal displacement due to the applied midspan load F can be calculated from the equation 2 The horizontal displacement caused by the unit force acting on the anchorage point can be obtained as The change of anchoring force From Eq. (7), the load F can be given in the form The deflection of the midpoint of a beam under the action of force F is Eq. (14) is the free vibration equation of prestressed beam with straight prestressed steel. It is derived by applying the theory of material mechanics and structural dynamics.
With method of separation of variables, the vibration solution to Eq. (14) can be assumed as For a simply supported beam, X(x) can be expressed as With Eq. (19), the natural circular frequencies of free vibration of a prestressed beam can be easily calculated. But it must be noticed that these formulas are based on simplified vibration control equations, and their accuracy is worth discussing.
In paper [2], Eq. (14) is rewritten as [2]   It is a complex partial differential equation with variable coefficients, which cannot be directly solved to obtain the analytical expression of natural frequency.
In paper [3] Solving the equation with method of separation of variables, the explicit formulae for natural circular frequencies can be obtained [3] . Comparing two equations, it can be found that the last two items of Eq. (21) disappear in Eq. (22). However, when x=l/2 is substituted into Eq. (21), the coefficient before y term is not zero, and this term mainly reflects the influence of eccentricity on vibration performance. So the accuracy of the analytical expression of natural frequency obtained by Eq. (22) is also worth discussing.
Considering the coefficient before y term of Eq. (21) is not zero, this paper corrected the Eq. (22)

Differential quadrature method
Eq.(26) is a complex partial differential equation, and the analytical expression of natural frequency cannot be obtained directly. In here, the method of DQ is employed to numerically solve the problem. This method requires to descretize the domain of the problem into N points.
Then the derivatives at any point are approximated by a weighted linear summation of all the functional values along the descretized domain. Consider the computational domain 0≤x≤1 of the beam. The number of sampling point is N in the x direction. Introduce N sampling point as [5]   1 0, 1, The quadrature rules for the derivatives of a function at the sampling points yield and in case of r =2, 3, ..., N-1 To overcome difficulties in the implementation of the boundary conditions, the modified weighting coefficient matrix are introduced as Substituting matrix (35) and (36) into Eq. (30), one obtain Eq. (37) leads to a generalized eigenvalue problem that can be numerically solved.

Numerical results
Consider a finite elastic prestressed concrete model beam, the parameters of the beam are listed in Table 1. Under pinned boundary conditions, from Eq. (37), the values of natural frequency can be found by numerical method. For the free motion of the model beam, the effects of parameters on the natural frequencies are respectively studied. In Table 2, the first four natural frequencies of the model beam with various concrete strength are given. For same concrete strength, with the growth of the orders, the natural frequencies increase greatly, but for the beam with higher concrete strength, the natural frequencies increase slightly. In Table 3, it shows that for same span length, with the growth of the orders, the natural frequencies increase significantly, but with the increase of span length, the natural frequencies decrease a little, and the natural frequency values of higher orders decrease more quickly. In Table 4 and Table 5, with the increase of the values of the original prestressing force, the natural frequencies decrease slightly, but with the growth of eccentricity, the natural frequencies increase a little.

Conclusions
This paper investigates natural frequencies of free transverse vibrations of a finite prestressed beams under simply supported boundary conditions, and the free vibration governing equations in related literatures are discussed and recorrected. As for the complex partial differential equation, the analytical expression of natural frequency cannot be obtained directly, the differential quadrature methods are applied to get the numerical results, and the effects of the physical parameters of the beam on the natural frequencies are respectively analyzed. The numerical results show that the natural frequency values increase with the growth of concrete strength and eccentricity of prestressed steels. But with the increase of the beam span length and values of original prestressing force, the natural frequency values decrease.