Experiment and simulation research of the ground-borne vibration for a high-speed train

Kai Cui1 , Dong Lv2

1Key Laboratory of High-speed Railway Engineering of the Ministry of Education, Southwest Jiaotong University, Chengdu, P. R. China

2School of Civil Engineering, Southwest Jiaotong University, Chengdu, P. R. China

1Corresponding author

Journal of Vibroengineering, Vol. 18, Issue 3, 2016, p. 1783-1797. https://doi.org/10.21595/jve.2016.16673
Received 29 November 2015; received in revised form 13 April 2016; accepted 21 April 2016; published 15 May 2016

Copyright © 2016 JVE International Ltd. 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.
Creative Commons License
Table of Contents Download PDF Acknowledgements References
Cite this article
Views 43
Reads 12
Downloads 847
WoS Core Citations 0
CrossRef Citations 0
Abstract.

In order to study the effect of the operational loads on the ground-borne vibration of the high-speed train, a train-track coupling model with considering the vertical and horizontal effects is established and applied to calculate the impact of different operational speeds on vibration acceleration. As shown in the results, the vibration acceleration is largely affected by different frequencies generated from different train speeds. By means of an indoor dynamic triaxial test, the impact of different vibration frequencies of a train on soil body is simulated. And a large number of medium- and low-frequency vibration tests are conducted according to the settings of load form, drainage requirement and vibration number of train vibration loads. The experimental results are analyzed to study the effect of different frequencies on dynamic characteristics, and a dynamic strain-time calculation formula, that takes the frequency factor into consideration, is proposed. Meanwhile, the improved formula that considered frequency is substituted into the finite element model (FEM) of the train, so as to analyze the impact of different vibration frequencies on the settlement, is applied. As shown in the results, the proposed improved formula, that considered the frequency, is good at prediction. The effect of vibration efficiency on the engineering can be reflected by a simulated high-speed train model. Based on the simulation model, a reinforcement measure is conducted for the ground-borne, and it is calculated that the settlement is obviously reduced and the service time of the train ground-borne is increased. This paper can provide a reference for a theoretical research and engineering practice.

Keywords: high-speed train, vibration frequency, indoor dynamic triaxial test, finite element model.

1. Introduction

With the economic development, the high-speed train has become one of main transport. Being stable, timely and secure, this transport means plays an important role in production and life due to its maximum operating speed which is over 200 km/h [1]. However, in the long-term vibration process, the high-speed train would impact the surrounding environment and produce uneven settlement on the soft soil ground-borne. Zarembski indicated that the long-term vibration would affect the water resistance and durability of railway buildings and structures, as well as the track regularity, ride comfort and normal operation [2]. As a result, Liang made a research on the dynamic characteristics of soft soil of the train [3].

Under the effect of vibration loads, there are many factors affecting the dynamic characteristics of soft soil [4]. In consideration of various influencing factors, Pan has proposed a method to calculate vibration loads [5]. In order to calculate vibration loads accurately, Feng has established a train-track coupling model that considers multiple factors, and compared the calculation value with the experimental value for verification [6]. In addition, the effect of different vibration factors on evaluation indexes is calculated, and a comparative analysis is conducted by Ma [7]. There are many researches regarding characteristics of the dynamic strain. Monismith C. has established an exponentially calculation formula of dynamic strain-time relationship curve, which could explain the dynamic strain-time relationship better. However, fewer factors are considered here [8]. Based on this, Li D. and Chai J. C. make improvements through considering confining pressure, vibration amplitude and other factors and conduct a verification analysis, which does not consider the frequency impact [9-11].

A wheel-track coupling model with considering the vertical and horizontal effects is presented in this paper. And by means of calculation and analysis, the impact of the train speed on the acceleration is obtained. Meanwhile, the dynamic strain-time relationship curve is improved, a dynamic strain-time relationship curve with considering the frequency factor is established, and parameter analysis is also carried out, whose calculation results can be applied to analyze the impact of different factors on characteristics of soft soil better. In addition, FEM software is used to simulate the ground-borne settlement of high-speed trains under the vibration load effect, thereby obtaining the initial dynamic stress. The settlement value is then obtained from the established dynamic strain-time formula with considering frequencies, and the comparison and verification between the experimental value and calculation value are conducted. Finally, in order to reduce sedimentation and improve the service time of soft soil, the reinforcement measure is taken for soft soil and FEM is used to calculate the settlement value. It is found in the result that the settlement can be effectively reduced and the service time of soft soil can be improved by this reinforcement measure.

2. Impact factor analysis

2.1. Establishment of wheel-track coupling model

In the operation system, the wheel-track system is a complex dynamic interaction [12]. Under the vibration loads, the periodical vibration loads will cause track vibration. The vibration on the contact point between the wheel and track is also caused by the load, and it also causes the elastic wheel-track deformation in the normal direction and creeping wheel-track slip change in tangential direction of the contact point. Besides, the wheel-track coupling vibration will be affected by the vibration of wheel-track force, thus affecting the system vibration. Therefore, a dynamic model for wheel-track coupling analysis needs to be established. The wheel-track coupling vibration can be divided into the vertical direction, horizontal direction and longitudinal direction [13]. The dynamic characteristics of the linear segment are mainly studied in this paper, and the vertical and horizontal effects are primarily considered. The coupling model diagram can be seen in Fig. 1.

According to the reference [2], the wheel-track coupling dynamic model is shown as below:

(1)
x ˙ = A x + F x , υ , t ,
(2)
F x , υ , t = 0 ,       0 ,       0 ,       - P c 1 m w ,       - P c 2 m w ,       P c 1 m w ,       P c 2 m w ,       0 ,       0 ,       0 ,       0 ,       0 T ,

where x is the operating train distance, A is the coefficient matrix, υ is the train speed, and t is the operating train time. Pc1 is the force between the first wheel-set and the track, while Pc2 is the force between the second wheel-set and the track. mw is the mass of the wheel-set. As known from nonlinear Herz contact theory [4], the wheel-track force is as follows:

(3)
P c = 1 G δ y ( t ) 3 / 2 ,

where δy(t) is the compressive modulus of elasticity between the wheel and track, and G is the wheel-track contact parameter, which is related to the elastic modulus of wheel radius R, wheel and track material as well as the Poisson’s ratio. For cone-tread wheel, G=4.75R-0.149×10-8, and for worn-type wheel, G=3.86R-0.115×10-8.

The compressive modulus of elasticity between the wheel and track includes the hydrostatic pressure of the wheel, which can be directly determined by the displacement of the wheel and track at the wheel-track contact point, as shown below:

(4)
δ y t = y w i t - y r x p i , t ,

wherein, ywi(t) is the displacement of ith wheel at time t, and yr(xpi,t) is the displacement of the wheel under the ith-wheel displacement at time t. With displacement irregularity input of y0(t) in the wheel-track interface, the wheel-track force can be displayed as below:

(5)
P c i t = 1 G y w i t - y r x p i , t - y 0 t 3 2 ,
(6)
y 0 = 1 2 a 1 - c o s 2 π x l ,

wherein, a is the depth of the irregularity, while l means the length of the irregularity. The impact of track irregularity is primarily considered in the paper, and other secondary factors are thus ignored; it is regarded that the wheel and track are in consistent contact during the operation process of the train, where the track irregularity is represented by the continuously unsmooth harmonic. Two wheels of a bogie have different positions in the irregular track surface, so as the wheel-track force is also different. According to the fixed wheelbase of the bogie, the irregularities of the first and second wheel-sets are as follows:

(7)
y 01 = 1 2 a 1 - c o s 2 π υ t - l 1 l ,
(8)
y 02 = 1 2 a 1 - c o s 2 π υ t - l 1 - 2 l b l ,

wherein, l1 is the irregular distance between the first wheel-set. The spectrum of the track irregularity is obtained from experiment, which can be substituted into the above equations so as to obtain the vibration acceleration under the wheel-track coupling.

Fig. 1. Vertical and horizontal contact models between wheels and track

 Vertical and horizontal contact models between wheels and track

Fig. 2. Track irregularity experiment

 Track irregularity experiment

2.2. Impacts of different speeds

There are three dynamic performance indexes to evaluate the safety and comfort of a train, which is acceleration, limited value of wheel-rail force and rate of wheel load reduction. In the paper, the index of acceleration is focused [14].

The track irregularity can be obtained from a bogie-mounted pitch-rate gyro. Sensors were mounted on the bogie, axle-boxes, as shown in Fig. 2. The spectrum of the irregularity of two wheel-sets can be obtained under the operating condition, whose results are shown in Fig. 3.

Fig. 3. Spectrum of track irregularity of two wheel-sets

Spectrum of track irregularity of two wheel-sets

a) Spectrum of track irregularity of first wheel-set

Spectrum of track irregularity of two wheel-sets

b) Spectrum of track irregularity of second wheel-set

The mass of the same train is fixed. Under the same condition, the impact of different train speeds is actually the impact of vibration frequency on the ground-borne [15]. The spectrum of the track irregularity obtained in the experiment is imported into the mentioned wheel-track coupling model of the train, so as to obtain the impact of the train speed on the acceleration, as shown in Fig. 4.

Fig. 4(a), (b) and (c) are corresponding to the acceleration-time relationship curves at the speed of 100 km/h, 200 km/h and 400 km/h, respectively. It can be obtained that the larger the train speed is, the greater the vertical acceleration and impact will be. And as the train mass is fixed under the same condition, the impact of the train speed on the train is also the impact of the vibration frequency on the train.

3. Analysis and calculation of indoor test

3.1. Test equipment

Fig. 5 is the dynamic triaxial equipment produced by UK GDS Company, which applies 5 Hz/10 kN/38-100 mm dynamic triaxial system (DYNTTS). The maximum diameter of the sample size is 300 mm, and the maximum confining pressure is 5 MPa. Therefore, the dynamic confining pressure test can be performed. In the vibration test, the soil body can be loaded with axial vibration force and set half-sine wave and other waves.

Fig. 4. Vibration acceleration under different speed

 Vibration acceleration under different speed

a) Speed = 100 km/h

 Vibration acceleration under different speed

b) Speed = 200 km/h

 Vibration acceleration under different speed

c) Speed = 400 km/h

3.2. Test process

Sample preparation: Soil is selected from an open space [16] near the high-speed train, and the soil is divided into silt and soft soil. The static pressure method is chosen to get the soft soil by mean of a thin-wall soil sampler. The sample preparation is made from the obtained undisturbed soil with the help of a cylinder with the cutting diameter of 3.8 cm and height of 8 cm [17].

Fig. 5. Standard dynamic triaxial equipment

Standard dynamic triaxial equipment

Fig. 6. Sample preparation process

 Sample preparation process

Experiment: The prepared soil sample is placed into an experimental instrument. The counter-pressure method is applied for saturation. And the reinforcement is set as 50 kPa, 100 kPa and 200 kPa, respectively.

The calculation value of vibration load-time relationship curve is studied. Through computation, the vibration waveform changing with time can be acquired. A fitting can be further made and it can be found that the half-sine wave is most consistent with the vibration load-time relationship curve [18], as shown in Fig. 7. Fitting is conducted on the experimental curve, because the waveform of the actual vibration load is difficult simulated in the indoor test. As only sine waves are provided by the experimental instrument for loading, as a result, sine waves are thus applied in the paper for fitting. Therefore, it can be shown that in the indoor test of the high-speed train vibration load, although waveforms are not completely consistent, the minimum and maximum peak loads are basically accordant with very similar change. Among them, the fitting function is shown as Eq. (9), and parameters and standard deviations are shown as Table 1. It can be seen that the smaller the standard deviations are, the better the fitting effect will be:

(9)
y = y 0 + A s i n π x - x c w ,

wherein, A>0. Other parameters can be found in Table 1.

Fig. 7. Calculation value and fitting value of excitation force

Calculation value and fitting value of excitation force

When vibration loads are added, the amplitude is implemented normalized process. The dynamic stress ratio is firstly defined as follows:

(10)
C S R = σ d 2 σ 3 ,

where σd is the amplitude of dynamic stress. σ3 is the confining pressure of the reinforcement. CSR is set to 0.1 and 0.2, and half-sine wave is selected as the vibration waveform. In order to analyze the effects of different frequencies, the vibration frequency is set to 0.25, 0.5, 1, 2.5, and 5 Hz, respectively, with the vibration number of 3,000 times.

Table 1. Parameters of sine fitting curve

Parameters
Value
Standard deviation
x c
0.076
0.00173
x
0.079
0.00003
w
82.02
0.025
y 0
105.82
0.26

3.3. Analysis of experimental results

By the indoor dynamic triaxial test, the vibration load of the high-speed train is simulated [19] to obtain that the dynamic strain-cycles relationship curve of different frequencies and dynamic stresses is a cluster of fluctuating curve, which is studied in many papers. MONISMITH has proposed an exponential strain curve as follows:

(11)
ε d = A N b ,

where A and b are parameters of the model determined by the least square method. N is the number of vibration.

Meanwhile, Chai J. C. has made a normalized process for the stress amplitude and confining pressure as follows:

(12)
ε d = a σ d σ s m N b ,

where σd is the load amplitude, and σs is the effective confining pressure. a and b are parameters of the model determined by the least square method. N is the number of vibration.

A more classic relationship model between strain and the number of vibration is shown in Eq. (12), which has been verified by experiments since it had been proposed [11]. Such a model is applied since more comprehensive factors are considered in the model. Besides the soil, the influence factors of dynamic characteristics for soil include dynamic load amplitude, confining pressure, the number of vibration and vibration frequency. However, all the influence factors except the vibration frequency are taken into account in Eq. (12) proposed by Chai J. C. To describe the impacts of different train speeds, Eq. (12) was improved in the paper so that the impact of vibration frequency could be considered. Numerical value of the model parameters is shown in Table 2, which is determined by Eq. (12) based on the least square method.

Table 2. Parameters of CHAI J C model

Confining pressure / kPa
Frequency / Hz
A
B
m
50
0.25
0.22
0.11
1.4
0.5
0.38
0.14
1.7
1
0.37
0.13
1.5
2.5
0.43
0.21
1.5
5
0.51
0.24
1.6
150
0.25
0.47
0.19
1.8
0.5
0.46
0.21
2.3
1
0.53
0.24
2.4
2.5
0.56
0.23
2.5
5
0.59
0.29
2.3
200
0.25
0.55
0.23
3.7
0.5
0.53
0.24
3.2
1
0.64
0.27
2.8
2.5
0.67
0.30
3.2
5
0.71
0.34
3.4

In the experimental process, only the original soil structure is applied, and the GDS dynamic triaxial instrument is employed to change the vibration load amplitude, confining pressure, vibration frequency and the number of vibration. After obtaining the original experiment data, the model parameters in Eq. (12) and Table 2 are used for calculations. Finally, curves in Figs. 8-10 can be obtained.

As can be seen from Fig. 8 to Fig. 10, the strain-cycles relationship curve is divided into three parts: The first stage is the rapid and elastic growth process. The second phase is the strengthened growth process, whose curve presents a slow increase and a turning point at this time. It is the elastoplastic stage. The third stage is the steady growth segment with a stable curve. With the increase of confining pressure, the stress amplitude is increased significantly. In addition, under the same confining pressure and different vibration frequencies, the strain value becomes smaller with the increasing frequency.

Fig. 8. Stress-cycles relationship of different frequency under 50 kPa

 Stress-cycles relationship of different frequency under 50 kPa

Fig. 9. Stress-cycles relationship of different frequency under 100 kPa

 Stress-cycles relationship of different frequency under 100 kPa

Fig. 10. Stress-cycles relationship of different frequency under 200 kPa

 Stress-cycles relationship of different frequency under 200 kPa

3.4. Accumulative strain analysis with considering frequencies

The impact of frequency cannot be considered in Eq. (11). Based on CHAI J C’s formula, a strain-time model with considering frequencies is established as follows:

(13)
ε d = a 0 A f σ d σ s m N b ,

where σd/σs=0.2, Af is the functional expression of frequencies. a, b and m are selected with the following values 0.84, 0.13 and 2.0 by referring Table 3.

To establish a relationship between parameter Af and the frequency [20], the least square method is used to calculate, obtaining a relationship curve as shown in Fig. 11.

Fig. 11. Parameter Af – frequencies relationship under different pressure

 Parameter Af – frequencies relationship under different pressure

In the figure, σs/σd=0.2 is chosen for research. The curve is conducted for the calculation and analysis, thus acquiring a calculation formula of parameter Af:

(14)
A f f = m - n l n f + l ,

wherein, f is the vibration frequency; m, n and l are calculating parameters of the formula, as shown in Table 4.

Table 3. Parameters suggested by Li and Seling

Type of soil
a
b
m
CH (high liquid limit clay)
1.2
0.18
2.4
CL (low liquid limit clay)
1.1
0.16
2.0
MH (high liquid limit silt)
0.84
0.13
2.0
ML (low liquid limit clay)
0.64
0.10
1.7

Table 4. Calculation parameters

Confining pressure / kPa
m
n
l
50
39.2
8.85
–0.21
100
18.1
2.6
–0.14
200
10.6
2.1
–0.08

4. Analysis for finite element calculation

4.1. Establishment of high-speed train model

Targeting at a high-speed train project, a research is carried out in the paper. The tested road structure is the earth-filled embankment. And the geometrical dimensions of the model are as follows. The length of top ground-borne layer is 13.4 m, thickness of the ground-borne and embankment is 2.9 m, and the design slope gradient is 1:1.5, with the length of 25 m and depth of 15 m recommended according to the references [21]. Boundary conditions are as follows: the left and right displacement is constrained in X-direction, and the displacement in the bottom is fixed. The vibration load amplitude was set in the load modulus during the loading process, and the near sine load is selected. There are 1621 nodes and 2500 elements in the model. CPE4P four-node pore pressure element is applied as the element type of the model. And the plastic constitutive modules that come with the software are used as the material: namely the Mohr-Coulomb material model, which has been set in the Property modulus of ABAQUS software as shown below. The figure is not beautiful enough. In addition, considering the length of the paper, as a result, it is not included in the paper. The parameters have been chosen with the following values: cohesion c= 11 kPa, internal friction angle φ= 9°, and dilatancy angle ψ= 0°. The permeability coefficient would be k= 10-9 m/s without a dilatancy phenomenon, the elastic modulus would be 207000 MPa, and the Poisson's ratio would be 0.3 [22].

In the analysis by means of experiments and theoretical calculations, the applied load is close to sine wave as shown in Fig. 7. Therefore, the load close to the sine wave should also be applied in the finite element model for more realistic simulation of this process. The load schematic diagram as shown in Fig. 12 is also presented by setting different parameters in ABAQUS software. For clear description, the corresponding control equations are also added in the paper. These parameters can be determined by Eq. (9) and Table 1 in the paper:

(15)
F t = A 0 + n = 1 N A n c o s ω t - t 0 + B n s i n ω t - t 0 ,

wherein, A0, An, Bn are constant values; ω is the circular frequency; t is the vibration time and t0 is the initial time.

Fig. 12. Load vibration waveform added into finite element model

 Load vibration waveform added into finite element model

The acceleration vibration load is a uniform load with the amplitude of 30 kPa. In order to facilitate its comparison with the experimental part, the number of vibration is set to 10 million times.

The CPE4P four-node pore pressure element is applied in the model. According to its stability condition, meshes are set under the Mesh module, and FEM of the ground-borne is established, as shown in Fig. 13.

Fig. 13. Finite element model of the ground-borne

 Finite element model of the ground-borne

4.2. Comparative analysis of settlement value

After the initial verification, the model is submitted for the finite element calculation, and stress contour is obtained as shown in Fig. 14.

In consideration of the long operation for the current high-speed train, the ground-borne has been compacted densely under the long-term vibration loads [23]. Therefore, the initial vibration times in Fig. 13 are assumed to be 10 million times, and initial value of the long-term deformation is set to 0. That is to say, only vertical deformation produced from more than 10 million times of vibration loads is studied [24]. As can be seen from Fig. 16, the vertical deformation of ground-borne is occurred slowly with the increase of the vibration times, and the average vertical deformation is about 0.08 mm to 0.14 mm for every one million times. Its stress-times relationship is also shown in Fig. 15. As can be seen from Fig. 15, when the vibration times increases to one million, the stress change is close to sine curve, which is similar to the applied loads in Fig. 12. Fig. 12 is the load applied into the finite element model. As shown in the figure, the load is similar to a sine curve. When the load is just applied into the finite element model, with the increasing time of the applied load, the foundation will bear the increasing vibration times and stress, as shown in Fig. 15. However, the load applied into the foundation will decrease as time continues to go on. In the meanwhile, the foundation will bear gradually decreasing stress. In this case, the final stress curve is obtained and shown in Fig. 15. The foundation stress calculated by reference [25] is similar to the result of this paper. The above analysis proves that Fig. 15 is reliable.

Fig. 14. Results of finite element simulation for the ground-borne

Results of finite element simulation for the ground-borne

Fig. 15 is substituted into Eq. (13) established herein, and then the layering method is applied to calculate the final settlement. Finally, the calculation value is also compared with the experimental value, as shown in Fig. 16. As can be seen from the result, the experimental value is consistent with the calculated value. As a result, the improved formula is reliable.

Fig. 15. Relationship curve between stress and vibration times for the ground-borne

 Relationship curve between stress and vibration times for the ground-borne

4.3. Reinforcement measure and calculation analysis

Through the above mentioned calculation, it is found that the predicted value is closer to the experimental value, and the FEM of ground-borne can make accurate prediction and calculation. In order to control and prevent settlement, the reinforcement measurement is taken [26]. There are a lot of holes in the foundation as shown in Fig. 17. Then, the concrete is injected into the holes. The water-cement ration is 0.55 and density is about 1.62 g/cm3. In addition, 0.20 % citric acid solution is added, and the temperature is maintained at 29 C°, finding that the initial setting time of the slurry is 1.65 hours, and calculi rate of slurry solution is 96 %. The flow rate of the injected slurry reaches 12 L/min. The single injection of the slurry solution is 150 L and 225 L, respectively [27-29]. Finally, the corresponding computation results can be obtained as shown in Fig. 18.

Fig. 16. Settlement-vibration times curves of the ground-borne

 Settlement-vibration times curves of the ground-borne

Fig. 17. Schematic diagrams of reinforcement strategy

Schematic diagrams  of reinforcement strategy

Fig. 18. Results of finite element simulation for the ground-borne after reinforcement

 Results of finite element simulation  for the ground-borne after reinforcement

Fig. 19. Relationship curve between stress and vibration times for ground-borne before and after reinforcement

 Relationship curve between stress and vibration times  for ground-borne before and after reinforcement

It is obtained from the finite element calculation that the stress value of the reinforced for the ground-borne is shown as Fig. 19, reaching 22.9 kPa in the maximum value. It is about half of maximum stress value before the reinforcement. The stress-times waveforms before and after reinforcement is similar.

Similar to the above method, the finite element stress value of the reinforced ground-borne [30, 31] is substituted into Eq. (13), and layered method is applied to calculate the settlement value, as shown in Fig. 20. The settlement value is controlled obviously.

It can be obtained through the calculation and comparison that the reinforced ground-borne is reduced by half in terms of the final settlement. It thus indicates that the uneven settlement can be effectively controlled by the reinforcement measure. As a result, the impacts of long-term vibration on ground-borne can be prevented.

Fig. 20. Settlement value of ground-borne before and after reinforcement

Settlement value of ground-borne before and after reinforcement

5. Conclusions

1. The train-track coupling model considering vertical and horizontal contacts is established, and the impact of the speed on the train acceleration is then analyzed.

2. Soil is extracted on the site of the ground-borne, and the indoor dynamic triaxial test is then carried out in order to simulate the load condition of the soft soil under the vibration load effect during operation. Different frequencies and vibration loads are set to get a dynamic strain curve and conduct regular research and analysis. In addition, improvements are made based on the calculation formula of this model, the impact of vibration frequency is considered, and the parameter analysis is conducted.

3. The improved exponential strain formula with considering frequency is then used in engineering. FEM is applied to simulate ground-borne, obtaining a stress value of the ground-borne. After the substitution, the corresponding dynamic strain is calculated. The layered method is applied to calculate the settlement value during operation, which is compared with the experimental value. As shown from the result, the improved model is suitable for analyzing the settlement calculation of the ground-borne.

4. In order to control and prevent settlement, the reinforcement measure is conducted. And the verified FEM is applied to calculate the reinforced stress and settlement values and compare them with the original results. As shown from the result, the performance of the ground-borne is significantly improved after the reinforcement.

Acknowledgements

This project is supported by the National Natural Science Foundation of China (Grant No. 41572245).

References

  1. Wu Chengjie The influence of regional subsidence in Suzhou-Wuxi-Changzhou and Shanghai areas on Beijing-Shanghai high-speed railway and its prevention countermeasures. Journal of Railway Engineering Society, 2007, p. 9-12, (in Chinese). [Search CrossRef]
  2. Zarembski A. N. Dynamic loading of the track structure. RT&S, Vol. 85, Issue 10, 1989, p. 11-13. [Search CrossRef]
  3. Liang B., Zhu D., Cai Y. Dynamic analysis of the vehicle sub-grade model of a vertical coupled system. The Journal of Sound and Vibration, Vol. 245, 2001, p. 79-92. [Publisher]
  4. Wei W., Yang X. L., Zhou B., et al. Combined energy minimization for image reconstruction from few views. Mathematical Problems in Engineering, Vol. 16, Issue 7, 2012, p. 2213-2223. [Publisher]
  5. Li Jian-qiang, He Suiqiang, Ming Zhong An intelligent wireless sensor networks system with multiple servers communication. International Journal of Distributed Sensor Networks, Vol. 7, 2015, p. 1-9. [Publisher]
  6. Feng Jun-he, Yan Wei-ming Numerical simulation for train stochastic vibration loads. Journal of Vibration and Shock, Vol. 27, Issue 2, 2008, p. 49-52, (in Chinese). [Search CrossRef]
  7. Ma Liheng, Liang Qinghuai, Gu Aijun, Jiang Hui Experimental study and numerical analysis on vibrations of sub-grades of Shanghai-Nanjing intercity high-speed railway. Journal of the China Railway Society, Vol. 36, Issue 1, 2014, p. 88-93, (in Chinese). [Search CrossRef]
  8. Chen J., Lin Q., Shen L. L. An immune-inspired evolution strategy for constrained optimization problems. International Journal on Artificial Intelligence Tools, Vol. 20, Issue 3, 2011, p. 549-561. [Publisher]
  9. Li D., Selig E. T. Cumulative plastic deformation for fine-grained sub grade soils. Journal of Geotechnical Engineering, Vol. 122, Issue 12, 1996, p. 1006-1013. [Publisher]
  10. Wei W., Yang X. L., Shen P. Y., et al. Holes detection in anisotropic sensornets: topological methods. International Journal of Distributed Sensor Networks, Vol. 21, Issue 9, 2012, p. 3216-3229. [Publisher]
  11. Chai J. C., Miura N. Traffic-load-induced permanent deformation of road on soft subsoil. Journal of Geotechnical and Geoenvironmental Engineering, Vol. 128, Issue 10, 2002, p. 907-916. [Publisher]
  12. Zhu Zexuan, Xiao Jun, Li Jian Qiang, Wang Fangxiao, Zhang Qingfu Global path planning of wheeled robots using multi-objective memetic algorithms. Integrated Computer-Aided Engineering, Vol. 22, 2015, p. 387-404. [Publisher]
  13. Heckle M., Hauck G., Wettschureck R. Structure-borne sound and vibration from rail traffic. Journal of Sound and Vibration, Vol. 193, Issue 1, 1996, p. 175-184. [Publisher]
  14. Yan Q., Yu F. R., Gong Q., Li J. Software-Defined Networking (SDN) and Distributed Denial of Service (DDoS) attacks in cloud computing environments: a survey, some research issues, and challenges. IEEE Communications Survey and Tutorials, Vol. 18, Issue 1, 2016, p. 602-622. [Publisher]
  15. Jinbao Yao, Xia He, Chen, et al. The impact of running trains on nearby buildings vibration test research and numerical analysis. China Railway Science, Vol. 30, Issue 5, 2009, p. 129-134. [Search CrossRef]
  16. Li Guohe, Jing Zhidong, Xu Zailiang A discussion of the correlation between land subsidence and groundwater level variation along the Beijing-Shanghai high speed railway. Hydrogeology and Engineering Geology, Vol. 35, Issue 6, 2008, p. 90-94, (in Chinese). [Search CrossRef]
  17. Wei W., Qiang Y., Zhang J. A bijection between lattice-valued filters and lattice-valued congruences in residuated lattices. Mathematical Problems in Engineering, Vol. 36, Issue 8, 2013, p. 4218-4229. [Publisher]
  18. Huang Bo, Ding Hao, Chen Yun-min cumulative deformation behavior of softclay in cyclic undrained tests. Chinese Journal of Geotechnical Engineering, Vol. 27, Issue 2, 2008, p. 331-338, (in Chinese). [Search CrossRef]
  19. Tang Yi-qun, Zhao Hua, Wang Yuan-dong, Li Ren-jie Characteristics of strain accumulation of reinforced soft clay around tunnel under subway vibration loading. Journal of Tongji University (Natural Sciences), Vol. 7, 2011, p. 972-977, (in Chinese). [Search CrossRef]
  20. Liu Ruimin, Ye Yincan, Chen Xiao Ling Study on the dynamic behavior of silt in shallow layer of Hang Zhou Bay. Journal of Marine Science, Vol. 31, Issue 3, 2013, p. 49-54, (in Chinese). [Search CrossRef]
  21. Dai Renping Study on cumulative settlement of cross-river tunnel of Hangzhou metro line 1. Modern Urban Transit, Vol. 3, 2014, p. 57-60, (in Chinese). [Search CrossRef]
  22. Lin Q., Chen J. A novel micro-population immune multiobjective optimization algorithm. Computers and Operations Research, Vol. 40, Issue 6, 2013, p. 1590-1601. [Publisher]
  23. Wang Jun, Cai Yuanqiang, Xu Changjie, et al. Study on strain degradation model in saturated soft clay under cyclic loading. Chinese Journal of Rock Mechanics and Engineering, Vol. 26, Issue 8, 2007, p. 1713-1719, (in Chinese). [Search CrossRef]
  24. Liu Tian-jun, Mo Hai-hong Strain rate of saturated soft clay under long term cyclic loading. Journal of South China University of Technology (Natural Science Edition), Vol. 36, Issue 10, 2008, p. 37-42, (in Chinese). [Search CrossRef]
  25. Wang Tingting Study on the Dynamic Characteristics of Soft Soil and the Long-Term Settlement Induced by the Dynamic Load. Thesis, Southeast University, Nanjing, 2014, (in Chinese). [Search CrossRef]
  26. Xu Yiqing, Tang Yiqun Experimental study on dynamic elastic modulus of reinforcing soft clay around subway tunnel under vibration load. Engineering Mechanics, Vol. 7, 2012, p. 250-255, (in Chinese). [Search CrossRef]
  27. Satoh T., Fushimi M., Tatsumi Y. Inversion of strain-dependent nonlinear characteristics of soils using weak and strong motions observed by borehole sites in Japan. Bulletin of the Seismological Society of America, Vol. 91, Issue 2, 2001, p. 365-380. [Publisher]
  28. Wong K. W., Lin Q., Chen J. Error detection in arithmetic coding with artificial markers. Computers and Mathematics with Applications, Vol. 62, Issue 1, 2011, p. 359-366. [Publisher]
  29. Du Z. H., Zhu Y. Y., Liu W. X. Combining quantum-behaved PSO and K2 algorithm for enhance gene network construction. Current Bionfinormatics, Vol. 8, Issue 1, 2013, p. 133-137. [Search CrossRef]
  30. Chen J. Y., Lin Q. Z., Hu Q. B. Application of novel clonal algorithm in multiobjective optimization. International Journal of Information Technology and Decision Making, Vol. 9, Issue 2, 2010, p. 239-266. [Publisher]
  31. Lin Q., Wong K. W., Chen J. An enhanced variable-length arithmetic coding and encryption scheme using chaotic maps. Journal of Systems and Software, Vol. 86, Issue 5, 2013, p. 1384-1389. [Publisher]