Research Article  Open Access
Dynamic Response Analysis of an Asymmetric Coupled VehicleTrack System Generated by Voided Elastic TwoBlock Sleeper
Abstract
Based on vehicletrack coupled dynamic theory, a threedimensional asymmetric vehicletrack coupling vibration model is developed to investigate the effect of voided elastic twoblock sleepers on vehicle and track system dynamic responses. For the vehicle system, one car body, two frames, and four wheel sets are assumed to be rigid, with 35 degrees of freedom (DOF). For the track system, the rails and the concrete twoblock sleepers are the main vibration components. The rails are modelled as Timoshenko beams, and the concrete twoblock sleepers are assumed to be rigid mass with vertical and lateral movement. The pads under the rails and the rubber boots under the sleepers provide greater vertical and lateral elasticity for the track. The Hertz nonlinear elastic contact theory is used to calculate the normal wheel/rail force. The wheel/rail tangent creep force is first calculated using Kalker’s linear creep theory and then modified by the ShenHedrickElkins theory. The results show that the asymmetric voided elastic twoblock sleepers have greater effects on the dynamic responses for fasteners and sleepers than on the car body and the wheel/rail forces under measured geometric irregularity and random irregularity. Two or more voided sleepers will greatly affect the vehicle running safety.
1. Introduction
In China, reinforced concrete monolithic track beds with twoblock sleeper are widely used in the urban transit railway. For vibration reduction, the concrete twoblock sleeper is usually placed into rubber boots to improve the elasticity of the track. Due to the adjusting error in the track construction, aging of the rubber pads, and the track beds being immersed in water during longterm operation, elastic twoblock sleeper becomes void and causes track stiffness irregularity, such as the sudden presence of hard spots in the rail support, as shown in Figure 1.
Voided sleepers are commonly found in ballasted tracks [1, 2], causing the sleeper to crush and the service life of the rail to be shortened. Recent studies [3, 4] have investigated the dynamic problems caused by voided long sleepers in ballasted tracks. Based on vehicletrack coupled dynamics theory, a twodimensional model was used to analyse the influence of voided sleepers or the sleeper span on the vertical displacements of the rails, the sleeper’s pressures, and the vertical wheel/rail forces. In addition, the effect of “hard spots” under the rail, induced by the supports of long sleepers being too high [5–9], on the change in the wheel/rail force and the pressure on sleepers was studied by using a threedimensional vehicletrack coupling dynamics model, and the results showed that the interaction force between the rail and the sleeper is sensitive to stiffness irregularity [10]. A threedimensional model was also used to analyse the influence of a voided sleeper on the dynamic responses of an adjacent sleeper in the ballasted track. When the sleeper is partially or completely voided, the maximum bending moment acting on the sleeper moves from the section under the rail to the middle section, from which it is easy to induce transverse cracks in the upper part of the middle section [11].
In the above research, it is obvious that the twodimensional or threedimensional vehicletrack coupling models were mainly established to analyse the influence of voided sleepers on the load increment of the adjacent sleepers in ballasted track with long sleepers. However, the influences of voided twoblock sleeper on the vehicle and track system have been hardly studied in the unballasted track of the urban transit railway. In addition, once one or more sleepers are voided, the support system under the rails will become asymmetric, which causes more dynamic problems in the vehicle and track system. To solve the problems caused by voided sleepers, based on vehicletrack coupled dynamic theory, a threedimensional asymmetric dynamic model is developed to study the coupling effect between the vehicle system and the track system with elastic twoblock sleeper monolithic beds under measured geometric irregularity and random irregularity. Through numerical simulation, the influence of a sudden change in support stiffness of the rail on the dynamic responses for the vehicle system and track system is also probed.
2. Dynamic Model
2.1. Model of VehicleTrack Coupling System
A threedimensional asymmetric vehicletrack coupling dynamic model [5] is shown in Figures 2 and 3. For the vehicle system, the car body, the frame, and the wheel sets, which are connected to each other with springs and dampers, are assumed to be rigid. The wheel set and the frame are connected by the primary suspension, while the body is supported on the frame through the secondary suspension. The primary suspension system includes the axle box springs, which have stiffness in three directions; the axle box positioning equipment, which has lateral stiffness; and the vertical primary dampers outside the axle box springs, which have lateral damping. The secondary suspension system includes the springs, which have stiffness and damping in three directions; the rubber blocks, which function as the lateral stop; the antiyaw dampers; the lateral dampers; and the vertical dampers. Therefore, the total degrees of freedom (DOF) of the vehicle model are 35, including lateral, vertical, roll, pitch, and yaw motions of the vehicle component (one car body, two frames, and four wheel sets). It should be noted that the longitudinal motions of the vehicle are ignored.
In Figures 2 and 3, , , , , and represent the lateral, vertical, roll, pitch, and yaw motions of the vehicle system, respectively; , , and are the mass of the car body, the frame, and the wheel set, respectively; and are, respectively, rolling and rotating moment of inertia of the wheel set; and are rolling and nodding moment of inertia of the frame, respectively; and are rolling and nodding moment of inertia of the car body, respectively; and are lateral and vertical stiffness of the second suspension at one side of the frame, respectively; and are lateral and vertical damping of the second suspension at one side of the frame, respectively; is the antirolling stiffness; and are lateral and vertical stiffness of each primary suspension, respectively; and and are lateral and vertical damping of each primary suspension, respectively.
The track system is modelled according to the elastic twoblock sleeper monolithic track bed. The rails and the concrete twoblock sleeper are the main vibration components. The rails are modelled as Timoshenko beams, and the concrete twoblock sleepers are assumed to be rigid mass, which have vertical and lateral movement. The pads under the rails and the rubber boots under the sleepers give the track greater vertical and lateral elasticity so that the track reduces the impact of the train loads. In Figures 2 and 3, is the mass of the rail; is the bending stiffness of the rail; is the mass of the sleeper; and are the vertical and lateral stiffness of a fastener and a pad, respectively; and are, respectively, the vertical and lateral damping of a fastener and a pad, respectively; and are the vertical and lateral stiffness of a rubber boot under the sleeper, respectively; and are the vertical and lateral damping of a rubber boot under the sleeper, respectively.
The wheel/rail interaction is the connection between the vehicle system and the track system. The rolling contact of the wheel and the rail generates the necessary conditions for a railway vehicle to run stably on the track. The rolling contact can be described by the kinematic constraints (geometric contacts) as well as the dynamic constraints (physical contacts), which are, respectively, referred to as the contact geometry relationship and contact forces between the wheel and the rail. In the present paper, the geometric problem of wheel/rail contact is solved spatially and evaluated online. The numerical method for the contact geometry calculation is discussed in detail [5].
The calculation of wheel/rail contact forces includes a normal model and a tangent model. The normal model, which characterizes the relationship law of the normal load and deformation between the wheel and the rail, is described by a Hertzian nonlinear contact theory [5]. The model by Shen et al. [12] is used as the tangent model determining the relationship between the creepages and the total creep forces of the wheel/rail. In the calculation of the tangent creep forces by Shen’s model, the creepages of the wheel and rail should be previously calculated in the transient analysis of the vehicletrack coupling dynamics. The creepage calculation is concerned with the parameters of the contact geometry of the wheel/rail, such as instant circle radii of the wheels and contact angles between the wheels and the rails.
2.2. Equations of Motion of the System
2.2.1. Equations of Motion of a Metro Vehicle
For the vehicle system, by using the coordinate system moving along the track with vehicle speed, the equations of motion of the metro vehicle can be easily derived according to D’Alembert’s principle, which can be described in the form of secondorder differential equations in time domain:where , , and are the vectors of displacement, velocity, and acceleration of the vehicle subsystem, respectively; is the mass matrix of the vehicle; and are the damping and the stiffness matrices which depend on the state of the vehicle system to describe the nonlinearities of the suspensions; and are the vectors of displacement and velocity of the track system; is the load vector acting on the vehicle system, composed of the nonlinear wheel/rail contact forces which are determined by the wheeltrack coupling model and depend on , , , and ; and is the external forces, including wheel/rail contact force, and gravitational forces and centripetal forces resulting from the vehicle running through the curve track.
2.2.2. Equations of Motion of the Track
In the track system, the rail is modelled as a Timoshenko beam. The forces acting on the rail are shown in Figure 4.
(1) The Rail. As for a rail on the right side of the track, the equation of vertical motion of the rail is
The equation of lateral motion is
The equation of vertical motion is given bywhere , , and are, respectively, the lateral, vertical, and torsion deflections of the rail, and are the slope of the deflection curve of the rail with respect to  and axes, is the rail density, is the area of the rail cross section, , is the shear modulus of the rail, is Young’s modulus of the rail material, and are, respectively, the second moment of area of the rail cross section about the  and axes, is the polar moment of inertia of the rail cross section, and , denoting, respectively, the lateral and vertical shear coefficients, which are obtained through the calculation of China’s rail cross section CN60 using the finite element method software ANSYS, and and are, respectively, the vertical and lateral wheel/rail forces between wheel and the rail, which changes with the vehicle speed , and are the force moment on the rail, is the Dirac delta function, is the coordinate of sleeper , is the coordinate of wheel , and and are the numbers of wheel sets and sleepers within the analysed rail. The subscript indicates sleeper , and wheel . and are the vertical and lateral forces between the rail and the th sleeper, respectively:
The displacements of the rail are written aswhere , , , , and are the modal displacement (coordinates) of the rail and NMR are the number of modes considered for the rail (). The frequencies of the modal shape functions of the rail are up to 1 kHz. The modal shape functions of the rail are given below:where is the calculation length of the rail.
By applying Ritz’s method, the fourthorder partial differential equation of motion of the rails then simplified to a series of secondorder ordinary differential equations in terms of the time coordinate, and , , , , and are
(2) The TwoBlock Sleeper. For concrete twoblock sleeper, the relationship between the forces is shown in Figure 5. Based on the interaction between the rail and the sleeper, the system equation is obtained.
(a) The left block
(b) The right block
The equation of the vertical motion of the left block is
The equation of the lateral motion of the left block is
The equation of the vertical motion of the right block is
The equation of the lateral motion of the right block iswhere , , , and are the reaction forces of the concrete twoblock sleeper:where and are the vertical displacements of the left and right block, respectively; and are the lateral displacement of the left and right block; and are the vertical velocities of the left and right block, respectively; and are the lateral velocities of the left and right block, respectively.
2.2.3. General Equations of Motion
All the above equations for the vehicle and track components are combined to form a set of ordinary differential equations expressed in terms of the following standard matrix form:where [], [], and [] are the generalized mass, damping, and stiffness matrices, respectively; , , and are the generalized displacement, velocity, and acceleration vectors, respectively; and is the corresponding force vector containing the wheel/rail forces ~.
2.2.4. Initial and Boundary Conditions of the Coupled VehicleTrack System
The deflections and the bending moments at the beamhinged ends are assumed to be zero. The deflection of the bed is ignored. The initial displacements and velocities of all components of the track are set to be zero. The initial displacements and the initial vertical and lateral velocities of all components of the vehicle are set to be zero, and the initial longitudinal velocity is the running speed of the vehicle, which is a constant.
2.2.5. Numerical Integration of the Coupled VehicleTrack System
Due to the strong nonlinearities involved in the equations of motion of the largescale coupled vehicletrack systems, the computation speed, accuracy, and stability of the numerical integration method are quite important when calculating the dynamic response of the coupled vehicletrack system. Zhai [13] proposed the “New Fast Numerical Integration Method” to specially solve the coupled dynamics equations of the railway vehicle and track. In this paper, Zhai’s method is used to solve the equations of motion of the coupled vehicletrack system with a time step size of s.
3. Results and Discussion
3.1. Influence of Voided Sleepers on Vehicle System
Due to the wear and damage of the rail, the nonuniformity of the sleeper span, the errors in the manufacture and installation of the fastener system, the stiffness changes in the pads, and other factors, the actual geometry of an urban rail line is usually random. In order to fully reflect the influence of voided sleepers of the elastic twoblock sleeper track on the track system and the vehicle system, the vibration responses of the track and the vehicle to voided sleepers and the typical track random irregularity are compared. The measured random irregularity, strictly controlled by construction quality and the errors of the track, is chosen as the excitation, which is the vector irregularity of the rail tread at the weld joint and the vertical profile irregularity of the rail surface, measured by the 10 m chordlength method, and the corresponding irregularity peaks of 0.3 mm and 4 mm, respectively. At present, there has been little research on the statistic features of the random irregularity of the urban transit railway, and there is no statistic random irregularity spectrum. So, the United States’ railway’s track irregularity spectrum is selected to simulate the track random irregularity of the urban transit railway in the present paper. The numerical simulation methods for various kinds of irregularities are described in [5]. In order to improve the computation efficiency, the vehicle concerned in this paper is the metroA series train, and the vehicle speed is 60 km/h in the numerical simulations. The parameters of the vehicle and track are given in [5].
The vehicletrack system is excited, respectively, by the geometric irregularity aforementioned and the stiffness irregularity caused by a single voided sleeper at the left side of the elastic twoblock sleeper track. The time history of the vertical acceleration of the mass centre of the car body is obtained by the developed model, as shown in Figure 6. It can be seen that the car body acceleration, excited by the vertical profile irregularity of the rail surface, is far greater than that excited by the voided sleeper and weld joint irregularity. The acceleration peaks of car body excited by the voided sleeper and the weld joint irregularity are equivalent, but the waveform cycle of acceleration of the car body excited by the voided sleeper is longer.
The time history of the vibration acceleration of a car body excited by different numbers of voided sleepers is shown in Figure 7. Figure 7 shows that the acceleration is sensitive to the number of voided sleepers, and the acceleration peak arises rapidly with the increase in the number of voided sleepers. By comparing Figures 6 and 7, it can be seen that even if there are three sleepers voided continuously, the vibration acceleration peak of the car body is lower than that excited by the vertical profile irregularity of the rail surface.
The maximal peaks of vibration displacement and acceleration of the car body and the frames by the stiffness irregularity caused by a single voided sleeper under the left rail, the measured geometric irregularity, and the random irregularity are, respectively, listed in Table 1. The data in the table show that the influences of the voided sleeper and the weld joint irregularity on the responses of the vehicle system are slightly different, and the influences of the voided sleeper are smaller; the vibration accelerations of the car body and the frame excited by both the vertical profile irregularity of the rail surface and the track random irregularity of the United States’ railway are relatively larger; however, the difference between the two excitations is small; the vibration accelerations of the car body and frame excited by the track random irregularity of the United States’ railway are, respectively, 8.9 times and 8.2 times that of those excited by the voided sleeper.

The direct influence of a single voided sleeper at one side of the track on the responses of the vehicle system is not obvious, compared to the measured geometric irregularity of the rail or the track random irregularity. So, it is not the factor affecting the vehicle running comfort and is not easily identified, until more sleepers on single side are voided continuously.
3.2. Influence of Voided Sleeper on Wheel/Rail Force
The wheel/rail force is an important evaluation index of dynamic interaction between the vehicle and the track. For the elastic twoblock sleeper monolithic track, a sleeper under the left rail is assumed to be voided. Deducting the mean value of the wheel/rail force, the vertical and lateral wheel/rail force time history are, respectively, shown in Figures 8 and 9. Due to the effect of the voided sleeper, both the vertical and lateral wheel/forces have peak impact, and the impact peak of the vertical wheel/rail force is far greater than that of the lateral wheel/rail force. So a single sleeper voided at one side has a greater effect on the vertical wheel/rail force than on the lateral wheel/rail force.
By comparing the left and right wheel/rail forces in Figures 8 and 9, it can be seen that although the sleeper is voided at the left of the track, the difference between the vertical and lateral wheel/rail forces is not obvious, and the amplitude of the right lateral wheel/rail force is larger than that of the left lateral wheel/rail force. At the same time, because of the poor bonding between elastic twoblock sleeper and the monolithic track bed, the influence of the voided sleepers is very strong.
The wheel/rail forces excited by the weld joint irregularity, the vertical profile irregularity, and track random irregularity of the United States’ railway are shown in Figure 10. The peaks of the vertical wheel/rail forces, from high to low, respectively, correspond to the weld joint irregularity, a voided sleeper, and sinusoidal excitation. The influence of the number of sleepers continuously voided under the left rail on the vertical wheel/rail force is shown in Figure 11. By comparing Figures 10 and 11, it can be seen that if there are three sleepers continuously voided, the peak of the vertical wheel/rail force is far greater than that excited by the weld joint irregularity. Using the developed model, the maximal vertical wheel/rail force excited by the track random irregularity of the United States’ railway is 74.7 kN, which is equivalent to that excited by the weld joint irregularity as shown in Figure 10, and is 3.8% smaller than that excited by three sleepers continuously voided. Therefore, in theory, only when there are three or more sleepers voided at one side can the change of vertical wheel/rail force be detected.
3.3. Influence of Voided Sleeper on the Adjacent Sleepers
The voided sleeper has a great effect on the fasteners and the adjacent sleepers, which can be evaluated by the pressures on the sleeper. When the ninth sleeper voided under the left rail is a voided sleeper, the maximal values of the pressures on the adjacent sleepers are shown in Figure 12. It is obvious that the pressures on the eighth, the tenth, and the eleventh sleepers increase significantly; in particular, the pressure on the tenth sleeper increases by 10 kN. The sleepers voided under the left rail have little effect on the pressure on the sleeper under the right rail.
When the ninth sleeper is voided under the two different conditions (voided at the left side only or on both sides), the pressures on the sleeper under the left rail are shown in Figure 13. It can be seen that, under the two different conditions, the influence of voided sleepers on the pressure on the sleeper is almost the same, the pressure on the tenth sleeper is the highest, and the pressures on the eighth and the eleventh sleepers are the next highest. But, under the condition of the ninth sleeper voided only under the left rail, the pressures on the eighth, the tenth, and the eleventh sleepers under the left rail are 0.3~0.4 kN higher than that under the condition of the ninth sleeper voided under both rails. Therefore, compared with the sleeper voided under both rails, the sleeper voided under one side only has a greater effect on the pressure on the sleeper.
The influence of the number of elastic twoblock sleepers continuously voided under the left rail on the pressures on the sleeper under the left rail is shown in Figure 14, which illustrates that the law of the pressure change of the adjacent sleeper caused by the sleepers continuously voided is consistent with that caused by a single voided sleeper. With the number of the sleepers continuously voided increasing, the pressure on the adjacent sleeper is higher. From Figure 14, it can be seen that the maximal pressures on a sleeper caused by a single sleeper and two and three sleepers continuously voided are found at the next adjacent sleeper, being 34 kN, 45 kN, and 61 kN, respectively. The fastener pressure is commonly designed at about 40 kN, so the pressures on the sleeper caused by a single elastic voided sleeper will theoretically satisfy the fastener pressure design requirement, but the pressures on the sleeper caused by two and more sleepers continuously voided will go beyond the design limits of the fastener pressure.
When a single elastic twoblock sleeper is voided under the left rail, the influence of different train speeds on the maximal pressure on the sleeper under the left rail is shown in Figure 15. From Figure 15, it can be seen that the relationship between the maximum pressure value and the train speed is approximately linear. The maximum pressure value increases by only about 0.14 kN with each 20 km/h increase in train speed; the sleeper pressure increases by about 0.56 kN when the train speed rises from 40 km/h to 120 km/h. So, for the urban railway with a maximum design speed of 120 km/h, it can be considered that the train speed has little effect on the maximal pressures on sleeper change caused by elastic twoblock voided sleepers.
The maximal sleeper pressures caused by different track irregularities are listed in Table 2. The vertical profile irregularity, weld joint irregularity, track random irregularity of the United States’ railway, and a single voided sleeper have linearly increasing influence on the maximal pressures on the sleeper. But the maximum pressure caused by a single voided sleeper is far greater than that caused by other irregularities.

4. Conclusions
Based on vehicletrack coupled dynamic theory, a threedimensional asymmetric vehicletrack coupling vibration model is developed to investigate the effect of voided elastic twoblock sleeper on vehicle and track system dynamic responses. The Hertz nonlinear elastic contact theory is used to calculate the normal wheel/rail force. The wheel/rail tangent creep force is first calculated using Kalker’s linear creep theory and then modified by the ShenHedrickElkins theory. Some conclusions are as follows:(1)A single voided elastic twoblock sleeper has the most obvious effect on the pressures on the sleeper, followed by the wheel/rail force, but has a very limited effect on the comfort of the car body.(2)For an urban railway with a low train speed, the pressure on a sleeper caused by it being a single voided sleeper is greater than that caused by the vertical profile irregularity, weld joint irregularity, or track random irregularity of the United States’ railway, but the pressures on the sleeper caused by only a single voided elastic sleeper still meet the design limits of the fastener pressure. However, in the case of two and more sleepers continuously voided, the pressures on the sleeper will go beyond the design limits of the fastener pressure, which may be harmful to the safe operation of the train.(3)The effects on the pressures on the sleeper, of the sleeper voided at one side and both sides, are nearly the same. But, compared to sleepers voided on both sides, the effect of a sleeper voided at one side has a slightly greater effect on the maximal pressure on the sleeper.(4)The sleeper voided at one side makes the vertical and lateral wheel/rail forces increase, with the vertical wheel/rail force increase being greater. By comparison, it can be seen that only when there are three or more sleepers voided at one side can this be detected by the change in the vertical wheel/rail force, but, at this time, both the fastener pressure and the pressure on the sleeper are beyond the design limits.(5)The voided sleepers have little effect on the vibration of the vehicle system. Even if there are three sleepers continuously voided, the vibration acceleration excitation peaks are lower than those excited by the vertical profile irregularity. Therefore, the voided sleepers are not easily identified by the vibration response of the vehicle system, and the track inspection should be enhanced during operation.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This research was supported by the National Natural Science Foundation of China (no. 51165017 and no. 51378395) and the Fundamental Research Funds for the Central Universities (Tongji University).
References
 E. L. Olsson and P. Zackrisson, LongTerm Measurement Results, Swedish National Road Administration, Borlänge, Sweden, 2002.
 S. Auguetin, G. Gudehus, G. Huber, and A. Schünemann, “Numerical model and laboratory tests on settlement of ballast track,” in System Dynamics and LongTerm Behaviour of Railway Vehicles, Track and Subgrade, vol. 6 of Lecture Notes in Applied Mechanics, pp. 317–336, Springer, Berlin, Germany, 2003. View at: Publisher Site  Google Scholar
 S. L. Grassie and S. J. Cox, “Dynamic response of railway track with unsupported sleepers,” Transport Engineering, vol. 199, no. 2, pp. 123–135, 1985. View at: Google Scholar
 M. Ishida, T. Moto, A. Kono, and Y. Jin, “Influence of loose sleeper on track dynamics and bending fatigue of rail welds,” Quarterly Report of RTRI, vol. 40, no. 2, pp. 80–85, 1999. View at: Google Scholar
 W. M. Zhai, VehicleTrack Coupled Dynamics, Science Press, Beijing, China, 3rd edition, 2007 (Chinese).
 C.H. Zou, S.H. Zhou, and B.L. Wang, “Calculation method of unsupported sleepers caused by differential subgrade settlement of ballasted track,” Journal of the China Railway Society, vol. 35, no. 1, pp. 87–92, 2013 (Chinese). View at: Publisher Site  Google Scholar
 Z. H. He, “The dynamic influences of the sleeper gap in urban rail transit,” Railway Engineering, vol. 3, pp. 86–88, 2010 (Chinese). View at: Google Scholar
 J. Xiang, H. Yang, and D. He, “Research on the responses of vertical vibration of traintrack system due to voided sleepers,” Journal of Railway Science and Engineering, vol. 4, no. 1, pp. 8–12, 2007 (Chinese). View at: Google Scholar
 A. Z. Jabbar, H. Xia, and J. J. Fan, “Effects of unsupported sleeper on dynamic responses of railway track,” Journal of Northern Jiaotong University, vol. 24, no. 1, pp. 50–55, 2000. View at: Google Scholar
 X. B. Xiao, X. S. Jin, and Z. F. Wen, “Effect of hard spot on dynamic response of railway sleeper,” Journal of Traffic and Transportation Engineering, vol. 24, no. 1, pp. 50–55, 2000 (Chinese). View at: Google Scholar
 J. Zhang, C. H. Wu, X. B. Xiao, Z. Wen, and X. Jin, “Effect of unsupported sleepers on sleeper dynamic response,” Journal of Southwest Jiaotong University, vol. 45, no. 2, pp. 203–208, 2010 (Chinese). View at: Publisher Site  Google Scholar
 Z. Y. Shen, J. K. Hedrick, and J. A. Elkins, “A comparison of alternative creep force models for rail vehicle dynamic analysis,” in Proceedings of the 8th IAVSD Symposium, pp. 591–605, Cambridge, Mass, USA, January 1984. View at: Google Scholar
 W.M. Zhai, “Two simple fast integration methods for largescale dynamic problems in engineering,” International Journal for Numerical Methods in Engineering, vol. 39, no. 24, pp. 4199–4214, 1996. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
Copyright
Copyright © 2016 Zhenxing He and Xinwen Yang. 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.