# On the Critical Speed, Supercritical Bifurcation, and Stability Problems of Certain Type of High-Speed Rail Vehicle

Authors: Hao Dong; Qunsheng Wang
Journal: Shock and Vibration
 (2017)

---

## Abstract

Recently, Chinese engineers have proposed a simple way to find the vehicles’ critical speed, which is similar to the ramping method. In this article, through an example of vehicle of supercritical properties, it is proved that the new easier way is not scientifically justified and should not be used in engineering practice. In addition, the ramping way also yields inaccurate critical speed. Then one abnormal vibration phenomenon which appears on Beijing-Shanghai high-speed line is studied. The results demonstrate that it is car body hunting but not bogie hunting. Finally, through the computation and comparison of the lateral ride indices under different conditions, one stability problem about stochastic limit cycle banding is tentatively discussed.


---

## Body

## 1. Introduction

The railway vehicle stability problem has been investigated for years. Since the nonlinear bifurcation theory is utilized in solving nonlinear vehicle problems [1–3], the study on the calculation of critical speed gets a great deal of achievements [4, 5]. The right way for calculating the critical speed refers to the path following (continuation) or ramping method [6]. Besides this, the typically wrong way is also discussed. In addition, in some manufacturer and research institutes [7], especially in China with its booming high-speed vehicle demand, the third simple way based on the continuation or ramping method is proposed and generally used. That method is only calculating the stable periodical solution, similar to ramping way, and not calculating the unstable limit cycle. Specifically, first the method involves yielding the large amplitude stable limit cycle, then decreasing speed with small stepsize to find smaller amplitude stable limit cycle, and decreasing speed stepwise until converging to the trivial solution (the trivial solution here is zero solution (as shown in Figure 3), which is different from periodical solution (as shown in Figure 11)). The critical speed is the corresponding value.However, the three ways for calculating critical speed have advantages and disadvantages. The continuation method is undoubtedly rigorous but time-wasting. The other two ways are fast but existing problems, one of which may be overshooting the bifurcation point [4]. Moreover, decreasing speed is a necessary step in both the ramping and the third method, as we know; we must determine the speed which makes the periodical solution converge to the trivial one. That is another problem where we usually calculate the inaccurate critical speed especially for the supercritical Hopf bifurcation type.Regarding CRH380B, one typical supercritical Hopf bifurcation type of high-speed vehicle in China, its dynamic models are depicted in Figure1. The examples of bifurcation diagrams of CRH380B on straight lines and its special conditions (such as one fault condition) are used to give further instruction as shown in Figure 2.Figure 1
Dynamic model for CRH380B.Figure 2
Lateral bifurcation of CRH380B’s wheelset ((a) and (b): normal conditions; (c): special conditions).
(a) (b) (c)Figure 3
Slow convergence rates for four wheelsets of one vehicle.When the ramping or the third method is used, it always yields the critical speed of about 300~350 km/h, even 400 km/h. This is because, with certain speed when we solve the vehicle dynamical system, we must ensure whether the result converges to the trivial solution (the constant zero), but sometimes, the convergent rate is very slow as shown in Figure3 obtained by Simpack simulation. Obviously, the result does not converge to 0 within 20 seconds.Especially in engineering applications, the computation is huge. Therefore, some researchers propose a simple way, which depends on the small vibration amplitude [8, 9]. If the result converges so slowly for some reasons (algorithm, singularity of model, correctness or complexity of model, etc.), they set a threshold of amplitude for the limit cycle to determine the critical speed. To be exact, the threshold is about 0.5 mm or even 1 mm as shown in Figure 2(b). Hence the last step when using ramping or the third method is that if the result does not converge to the trivial solution but its amplitude is smaller than 0.5 mm, they regard the speed as the critical one. To distinguish it from the critical speed Vcr in theory, we define the above-mentioned speed which is used in engineering as VE. In most subcritical Hopf bifurcation situations, VE undoubtedly equals Vcr. But however, this is not true in the case of some supercritical bifurcation as shown in Figure 2, which leads to error estimation of Vcr for nearly one hundred kilometers per hour. So when VE is used instead of Vcr, we must pay attention to the situation Vcr = Hopf bifurcation point (linear critical speed).The main work of this paper involves the critical speed of supercritical Hopf bifurcation of types CRH380B (Prototype of ICE3) and some stability problems.


## 2. The Disadvantage ofVE

As we know, CRH3’s prototype is the ICE3 train. China’s railway authorities modify it as a higher operation speed train, CRH380B, with structure or parameters optimization. One of the most important parameters is the bogie primary radial positioning (location node of box arm) stiffnesskpx. Figure 4 shows the influence of that on VE and Vcr.Figure 4
Influence of primary location stiffness onVE and Vcr (Vcr = Vlin-cr = Vnon-cr = Hopf, where Vlin-cr is the linear critical speed and Vnon-cr is the nonlinear one).From Figure4, it can be seen that the laws of VE and Vcr are different and even opposite. Exactly, with the stiffness becoming large, VE decreases while Vcr increases. Particularly, the laws of VE are believed to be correct according to some research reports [10, 11]. As a result, considering that CRH380B is designed to obtain higher operation speed (350–380 km/h) than CRH3 based on CRH3’s structure or parameters, if VE is used instead of Vcr, it is unscientific and may be dangerous.To provide further discussion onVE and Vcr and to prove the inaccuracy of VE, we choose two cases.Case  1. kpx = 30 MN/m, Vcr = 370 km/h < VOp, VE = 510 km/h > VOp.Case  2. kpx = 115 MN/m, Vcr = VE = 386 km/h > VOp.V O p = 380 km/h is the operation speed.Now, the eigenvalues and eigenvectors of CRH380B vehicle are calculated by Simpack software.Case 1.
The results are depicted in Table1.
In Table1, the left column denotes the nature damping. This can be simply understood as corresponding eigenvalues multiplied by -C (C is a positive constant). The right side denotes the eigenfrequency. So, if a negative value of nature damping appears, similar to the positive eigenvalues, thelinearized system of CRH380B is unstable. In addition, from Figure 2, because Vcr = Vlin-cr = Vnon-cr = Hopf, thenonlinear system of CRH380B is also unstable.Table 1
Eigenvalues forkpx = 30 MN/m.
Nat. damping Frequency (Hz) 1 0.2325 1 0.7672 0.0285 1.0079 0.0606 1.8823 0.0561 1.9632 0.0404 2.0079 0.0731 2.3198 0.9899 2.5590 0.9831 3.3193 0.3636 3.4589 0.4235 3.9598 −0.0043 7.0055 −0.0086 7.0253 0.7082 9.3246 0.7076 9.3303In Table1, there are a pair of negative values of nature damping and its corresponding frequencies of about 7 Hz. The corresponding eigenvectors are shown in Figure 5.Figure 5
Eigenvectors for the unstable bogie (solved by Simpack).From Table1 and Figure 5, it can be seen that the rear bogie loses stability at this operation speed with frequency of about 7 Hz (the other negative value of nature damping corresponds to the front bogie). These imply that, in the case of kpx = 30 MN/m, the critical speed is undoubtedly less than 380 km/h and VE = 510 km/h > 380 km/h is not safe.Case 2.
The results are depicted in Table2.
There is no negative value in the left column in Table2, which means the critical speed is greater than 380 km/h.Table 2
Eigenvalues forkpx = 115 MN/m.
Nat. damping Frequency (Hz) 0.0280 1.0763 0.0547 1.9598 0.0470 2.0473 0.0347 2.0604 0.0570 2.5242 0.3682 3.4584 0.4274 3.9547 0.9621 4.4891 0.9555 4.8816 0.0075 6.6075 0.0037 6.6300 0.6168 11.571 0.6162 11.576 0.2763 11.991 0.2725 12.0062 0.4537 15.7025Combined with the analysis of Cases1 and 2, VE is unscientific and cannot be used in vehicle design.


## 3. One Possible Reason for the Difference

From Figure4, it can be seen that(1)Vcr<VE,whenkpx<115MN/mVcr>VE,whenkpx>115MN/m.Should it be doubted that some limit cycles such as those in Figure3 are not “real” by macroscopic observation? It may be caused by the algorithm accuracy.Now, set parameterskpx = 120 MN/m (designing default) and the vehicle speed as 370 km/h < Vcr = 380 km/h according to Figure 4. The solving algorithm is SODASRT [12–14], and its accuracy is, respectively, 10−5 for Case A and 10−6 for Case B. The initial condition is depicted in Figure 6, the detail of which is the lateral track irregularities in the first 3 seconds and no irregularity afterwards. The response result is depicted in Figure 7.Figure 6
Initial condition: lateral track irregularities with the first 3 seconds.Figure 7
Response for different algorithm accuracy (first wheelset on straight line).From Figure7, it is unimaginable that solving the same equation and parameters yields different results if we ignore the algorithm accuracy. Specifically, Case A converges to the small amplitude limit cycle and Case B converges to the trivial solution. However the vehicle speed is 370 km/h < Vcr. So Case B is correct. This situation is common in engineering because we usually set the accuracy as default but not higher precision in order to solve the system faster. As we know, the higher precision may make the solving process very slow. In engineering, it is usually ignored.According to the above, if we useVE instead of Vcr, the error may also be caused by algorithm accuracy. Furthermore looking at Figure 2, the interval between VE and Vcr may be one hundred kilometers per hour and the algorithm accuracy’s effect may exist. Therefore VE, which is obtained by macroscopic observation based on simple bifurcation analysis, is unreliable.There is a great deal of the research about the influence of algorithm accuracy on the response of mechanical systems, especially in aircraft [15]. It is possible to apply the same research technique to rail vehicles.


## 4. One Problem about Stability

According to the two-year test for CRH3 operating on the Beijing-Shanghai high-speed line, an abnormal vibration phenomenon of car body appears on some sections, which makes passengers feel uncomfortable and complain to the railway authorities. Thetest results for car body are depicted in Figure 8.Figure 8
Response for vehicle vibration ((a): time domain; (b): frequency domain).
(a) (b)In the beginning, this phenomenon was regarded as what was coming from bogie hunting transmission, although the bogie vibration has the same frequencies (about 1.6–1.8 Hz) and modes. However, the frequencies of about 1.6–1.8 Hz are far from hunting ones which are around 7 Hz as shown in Table1. Therefore, we make further analysis by Simpack simulation for the whole vehicle. Firstly, because this vibration phenomenon only exists on special sections, we doubt that the gauge variation was the key. The influence of gauge variation on critical speed is depicted in Figure 9. PointP corresponds to the nominal track gauge. Here in CRH380B vehicle, for its supercritical Hopf bifurcation type as shown in Figure 2(b), Vlin-cr = Vnon-cr = Hopf bifurcation point.Figure 9
Influence of gauge variation on critical speed.From Figure9, whenever the gauge is smaller than 1430 mm or greater than 1442 mm, the critical speed is smaller than 300 km/h. In particular when the gauge increases to 1444 mm, the critical speed is very low at around 150 km/h, which cannot guarantee safety. Moreover, the details of the influence of gauge on stability were studied by Zboiński and Dusza [16].Secondly, assume the gauge increases to 1450 mm and the operation speed is also 300 km/h. The results for the whole vehicle simulation are shown in Figure10.Figure 10
Time domain response on condition gauge = 1450 mm and speed = 300 km/h.Figure 11
The unstable modal vectors (solved by Simpack).
(a)
Frequency = 1.6 Hz (b)
Frequency = 1.8 HzIn Figure10, it is obvious that the simulation results of wheelset, bogie, and car body have the same resonance frequency. Specifically, the first two have nearly the same vibration wave shape for the large positioning stiffness (primary node stiffness = 120 MN/m). Besides the frequency is about 1.6 Hz, similar to what is shown in Figure 7. Although the wave shapes look like bogie hunting, the frequency is far from that. According to Figure 8, we know that this result corresponds to the instability case because Vcr < 300 km/h. Why do the results look like bogie hunting and why are the values of frequency far from the hunting ones? Then, the eigenfrequency and eigenvectors are calculated as shown in Table 3.Table 3
Eigenvalues for gauge = 1450 mm and speed = 300 km/h.
Nat. damping Frequency (Hz) 0.0482 0.8972 −0.0331 1.6035 −0.1919 1.8363 0.2934 2.0367 0.0471 2.0457 0.4043 2.0827 0.0571 2.5222 0.3675 3.4603 0.4269 3.9571 0.8586 9.9214 0.8541 10.140 0.6176 11.553 0.6169 11.559 0.2766 11.982 0.2727 12.057 0.4045 14.172From Table3, the vehicle has two negative nature damping values which correspond to positive eigenvalues and the eigenfrequencies are 1.6 and 1.8 Hz. Therefore, the simulation analysis results correspond to the test results as shown in Figure 8. Furthermore, the unstable eigen modal vectors are depicted in Figure 11.Now it is obvious that, from Figure11(a), the 1.6 Hz vibration corresponds to the car body yaw motion and the 1.8 Hz vibration corresponds to the car body rolling motion. Therefore, the abnormal vibration which makes passengers uncomfortable is car body hunting but not bogie hunting. Also this phenomenon only appears on some special sections. So according to Figure 9, in the author’s opinion, it may be caused by gauge variation on some part of the whole Beijing-Shanghai high-speed line.


## 5. One Extended and Unknown Problem

Perhaps the following problem exists in vehicle dynamic simulation. Two types of vehicles are running in different conditions as follows. Vehicle A behaves as a subcritical Hopf bifurcation, while vehicle B behaves as a supercritical one.Case 1.
Vehicle A is running below its critical speed but is excited by track irregularity (both lateral and vertical). Therefore, the time domain response is random.Case 2.
Vehicle A is running above its critical speed with initial value as shown in Figure6. Its time domain response is periodical solution.Case 3.
The same vehicle is also running at the same speed as in Case2 and also is excited by track irregularity as Case 1. Its time domain response behaves like random signal adding to periodical solution.
The results of the three cases are depicted in Figure12.
Now, the lateral ride indices for the three cases are calculated as shown in Table4. The lower value demonstrates the higher ride quality (the ride index is similar to comfort index. The former is made by Chinese standard GB5599-85 and contains separate lateral and vertical computation processes [17], while the latter refers to UIC 513 code [18]).
The results shown in Table4 seem to be reasonable because the vibration level gets worse.
However, the following cases are not the same.Table 4
Ride indexes for vehicles A and B.
Case Lateral ride index 1 1.55 2 2.33 3 2.38 4 1.53 5 1.61 6 1.56Figure 12
Time domain responses for vehicle A.Case 4.
Vehicle B is running below its critical speed but is excited by track irregularity (both lateral and vertical).Case 5.
Vehicle B is running above its critical speed with initial value as shown in Figure6.Case 6.
Vehicle B is also running at the same speed as in Case5 and also is excited by track irregularity as Case 4.
The lateral ride indices for the three cases of vehicle B are also calculated as shown in Table4.From the result, the index for Case6 is better than Case 5. We may guess that the time domain response of Case 6 is also random signal adding to periodical solution. But the results in Figure 13 demonstrate that it is not as we expect from Table 4.Figure 13
Time domain responses for vehicle B.From Figure13, it can be seen that the result of Case 6 is different from that of Case 3 (periodical solution + random signal). The curves of Case 6 seem more like random signals. Considering that the two cases operate under the same condition, only the vehicle type is different, although we know that they have different bifurcation types. Why do Cases 3 and 6 produce different results?This leads to a problem: how can we predict the behavior of a nonlinear system loss of stability when excited? Here the bogie lateral acceleration for one high line is measured as shown in Figure14.Figure 14
Field test result for the bogie lateral acceleration of CRH380B (a section on Beijing-Shanghai line with operation speed 299 km/h).It is obvious that the test result is similar to Case3 but not Case 6. According to the criterion used for online surveillance, specifically the peak value of lateral acceleration on the bogie frame, as defined in the (now invalid) version of UIC 515 [19], the limit value is seen to exceed when the value reaches 8 m/s2 for more than six consecutive cycles with a frequency of 4–8 Hz, which means the test result for the vehicle is instable. Oldrich Polach did detailed study on bogie stability assessment [20].However, although we have not measured the instability results online for the vehicle as mentioned in Case6 because it may have much higher critical speed, it can be supposed that the results may have similar wave profile as shown in Figure 13 to Case 6. From Figure 13, Case 5 demonstrates hunting, which may also be tested as instability according to the criterion of UIC 515. Despite the fact that peak value of Case 6 is obviously higher than that of Case 5, according to UIC 515, it does not satisfy the condition that is six consecutive cycles. Therefore, this means that although Case 6 has large disturbance from track irregularity and has worse ride index and while Case 5 has no excitation and has better ride index from Table 4, Case 6 may be tested as stability while Case 5 would be tested as instability. As we know, UIC 515 (1984) is now invalid, but similar standards such as TSI L 84-2008 [21] and one Chinese standard for high-speed vehicle [22] are still used. So the instability test method of measuring the peak value of lateral acceleration and its consecutive cycles on the bogie frame needs to be modified especially for the high-speed vehicle.From the analysis above, there is no doubt that the external excitation has an influence on the limit cycle, for the stochastic response adding to the periodical solution, which would make it hard to determine instability (the result in Case5 compared with that in Case 6). As we know the critical speed is independent of the external excitation. But we should investigate the vehicle stability considering the track irregularity by applying the stochastic stability theory. The latter mainly contains two aspects. One is about parametric excitation. Bigoni studied the critical speed considering that and gave each possible critical speed’s probability [23]. The other is about external excitation. This is very difficult in nonlinear stochastic theory; for example, the method of solving the stochastic differential equation (SDE) and the stability theory for SDE have been now only preliminarily studied in mathematics. So far, only Xiang and Zeng’s work [24] uses this theory to investigate the derailment problems and Zhang’s work gave the primary result that the limit cycle is not an invariant curve but becomes a banding [25, 26]. Figure 15 will give an illustration.Figure 15
Phase plane for one vehicle.In Figure15, the red dot banding is the result of a simple track irregularity simulation as White Gaussian Noise. In fact the actual track excitation may be much more difficult to simulate as stochastic process and to compute. However, evidently from this figure, according to the standards mentioned above, the black results imply instability, but the red results mean uncertainty because it may not meet the condition of more than six consecutive cycles exceeding 8 m/s2.


## 6. Conclusion

Firstly, the new fast method proposed by Chinese engineers for computingVE instead of Vcr is unscientific and inaccurate through an example of a vehicle with slow converge rate. In addition the ramping way for computing Vcr may also lead to error estimation for nearly one hundred kilometers per hour in types of some supercritical Hopf bifurcation vehicle.Besides, from the bogie lateral acceleration test and reports on the Beijing-Shanghai high-speed line, the abnormal vibration of nearly 1.8 Hz which makes passengers uncomfortable is from car body hunting, not bogie hunting, according to the eigenvalue and gauge variation analysis. In the authors’ opinion, the cause of this phenomenon may be track gauge increase.Finally, with preliminary application of the stochastic and nonlinear theory compared with pure nonlinear theory for rail vehicle system, standards such as UIC 515 and TSI L 84-2008 for assessment bogie instability by measuring the peak value of lateral acceleration need to be modified.



---

*Source: 1526373-2017-07-19.xml*