Abstract

Ladder track, which has drawn increased attention in scientific communities, is an effective method for reducing vibrations from underground railways. In order to optimize the vibration reduction ability of ladder track, a new method, that is, the finite element method (FEM) coupled with ant colony optimization (ACO), has been proposed in this paper. We describe how to build the FEM model verified by the vibration tests in the Track Vibration Abatement and Control Laboratory and how to couple the FEM with ACO. The density and elasticity modulus of the sleeper pad are optimized using this method. After optimization, the vibration acceleration level of the supporting platform in the 1–200 Hz range was reduced from 102.8 dB to 94.4 dB. The optimized density of the sleeper pad is 620 kg/m3, and the optimized elasticity modulus of the sleeper pad is 6.25 × 106 N/m2.

1. Introduction

Vibrations generated by underground railways are one of the most serious engineering problems of such systems. Waves induced by the dynamic interaction between the train wheels and the rails propagate from the surrounding soils to the foundations of nearby buildings [1], resulting in structural vibrations and reradiated noise. One effective method for reducing vibrations from underground railways is to use ladder tracks.

The idea for ladder track was originally from the baulk road system and was then applied in the Leeds and Selby Railway in 1830. From the middle of the 20th century, systematic research on ladder track was started in Japan, Russia, and France. In the last 10 years, the vibration reduction performance of ladder track by the Railway Technical Research Institute of Japan has attracted great attention from Asian researchers, especially in China.

In the early 1990s, Wakui et al. [2] reported its new development in Japan. Oyado et al. [3] analyzed running performance and dynamic settlement test results of a ballasted ladder track. Initially, the ladder track was usually used as a ballasted track. Hosking and Milinazzo [4] built a simple mathematical model to study ladder track responses under steadily moving loads by employing the periodic discrete elastic support of the combined floating rails. This model was extended to include the support mass and viscous damping [5]. Hui and Ng [6] measured the vibration velocity level at the station using the ladder track. The results showed that the ladder track has better mitigation properties above 35 Hz, compared with ballast track. The first bending resonance of ladder track is 31.5 Hz under a moving train load. In the Beijing subway system, Xia et al. [79] and Inoue et al. [10] compared the dynamic response of an elevated bridge having ordinary nonballasted slab track with a bridge having ladder track. Theoretical analysis and experimental study proved that the ladder track has good vibration reduction characteristics. Jin et al. performed a modal analysis of the ladder track with different bearing forms by a numerical method [11] and a laboratory test [12]. The results showed that the first natural frequency of ladder track occurs at 33–36 Hz. To solve the problem of rail corrugation, the dynamic properties of the ladder track used in the Beijing subway were optimized by Yan et al. [13, 14].

The present contribution aims to optimize vibration reduction ability of ladder track by FEM coupled with ACO. This paper is organized as follows. Section 2 describes a vibration test of ladder track performed in the Track Vibration Abatement and Control Laboratory. In Section 3, a description is given of an FEM model that was built using the commercial software LS-DYNA and verified based on the tested ladder track. Then, the FEM model coupled with ACO is introduced in Section 4. In Section 5, the optimization of the vibration reduction performance of ladder track is explained. Conclusions are given in Section 6.

2. Vibration Testing in the Laboratory

The Track Vibration Abatement and Control Laboratory of Beijing Jiaotong University is a two-story underground structure for researching track vibration and is the only such facility in Asia. The buried depth of Lab #1 is 6 m, and the buried depth of Lab #2 is 14 m. Figure 1 shows the ground conditions around the laboratory.

One unit of ladder track was constructed in Lab #2 (Figure 2). The cross section of Lab #2 is in the shape of a horseshoe. The height is 4 m, and the width is 4 m.

The test track used 60 kg/m rail, which was 6.15 m in length. Ten WJ-2 fasteners were employed to fix each rail in place. The size of the longitudinal sleeper was 6.15 m × 0.46 m × 0.185 m. The size of the connecting beam, which was designed to connect the two longitudinal sleepers, was 0.975 m × 0.06 m × 0.06 m. There were five sleeper pads under each longitudinal sleeper. The size of the sleeper pad was 0.46 m × 0.25 m × 0.03 m. The size of the supporting platform was 0.86 m × 0.4 m × 0.28 m. Ten supporting platforms were used for the test ladder track. Figure 3 shows the plan of the test ladder track. Figure 4 illustrates the cross section of the test ladder track.

As the vibration source, an automatic falling weight machine, shown in Figure 5, was designed and was employed to provide impulse to the rail. By changing the number of mass blocks and the drop height, different impulse forces could be obtained. The material of the hammer head could also be changed to aluminum, rubber, nylon, or steel. To avoid influencing the sleeper by the mass of the setup, a scaffold was installed to support the automatic falling weight machine.

For this test, five mass blocks, 73 kg in total, were installed. The aluminum hammer head was used and the drop height was 10 cm. A force sensor was installed in the hammer head. The impulse site is the middle of one rail (see Figure 3). The sampling frequency of the force signal was 12.8 kHz. The time history of the average impulse force, with a peak value of 55 kN, is shown in Figure 6. The vibration acceleration sensor was fixed on the middle supporting platform to measure the vertical vibration, as shown in Figure 3. The sampling frequency of the acceleration signal was 1600 Hz.

3. FEM Model Building and Verification

3.1. Building the FEM Model

The commercial software LS-DYNA was employed to build the FEM model of the test ladder track. The geometry of the longitudinal sleepers, the connecting beams, and the supporting platforms are described in Section 2. Considering the calculation time of the LS-DYNA software, the cross section of the rail was simplified, as demonstrated in Figure 7.

The properties of the rails, the longitudinal sleepers, the supporting platforms, and the connecting beams are presented in Table 1.

Springs were used to simulate the fasteners and the sleeper pads. The vertical stiffness of the fastener was 60 MN/m. The vertical stiffness of the sleeper pad was 18.055 MN/m. Figure 8 shows the FEM model established by LS-DYNA. The bottom of all supporting platforms was constrained in all directions.

3.2. FEM Model Verification

The impulse force shown in Figure 6, which was obtained from laboratory testing, was used in the FEM model. Then, the vibration acceleration of the middle supporting platform was calculated. Figure 9 shows the calculated result and the test result. From Figure 9, we see that the calculated and the test results are almost the same. Therefore, the FEM model is accurate for the following calculation in Section 4.

4. FEM Coupled with ACO

4.1. Ant Colony Algorithm

With the development of computer technology, swarm intelligence optimization algorithms, including ant colony optimization (ACO) and particle swarm optimization (PSO), have attracted increased attention. ACO is one of the most successful swarm intelligence optimization algorithms. It was proposed by Colorni et al. [15] to solve the traveling salesman problem in 1991. It was named the ant system (AS), inspired by the foraging behavior of the Argentine ants [1618]. In recent years, ACO has been used to successfully solve many combinatorial optimization problems [19, 20] and is being extended to obtain the solutions of continuous problems [21]. There are many different ACOs based on the AS. In this paper, Chen’s ACO [22] was employed to solve the continuous function. Here is a function of one variable describing Chen’s ACO.

is the original objective function, and is the original design variable whose minimum (maximum) value is a (b). Let be by simple mathematical transform, in which , .

Therefore, the optimization process for function is simplified to an artificial ant that makes a selection from ten decimal numbers whenever it takes a step, except for the first place and the last place (see Figure 10).

Each artificial ant goes from the first floor (the nest) toward the last floor (the food source). is the decimal number when ant is at floor . Ant selects the decimal number of the next floor according toin which is a real random variable uniformly distributed in the internal array , is a tunable parameter controlling the influence of the pheromone, randperm is a number randomly selected from , and is the pheromone intensity laid on the path between the number to .

When an ant finishes one step, that is, an ant arrives at floor from floor , the strength of the pheromone laid on the path between the number and the number by ant should be updated as follows:where is the modified coefficient of the pheromone intensity and is the local evaporation factor of the pheromone.

After all the ants have arrived at the last place, the best solution for a given interior can be obtained. Through all interiors, the best solution of the function can be obtained.

According to the function tests, Chen’s ACO parameters are , . The interior number In is 10, and the ant number An is 10.

4.2. FEM Coupled with ACO

The calculation code for Chen’s ACO was written in Matlab. The FEM model was carried out through the DOS functions in Matlab. The detailed program flow is demonstrated in Figure 11.

5. Optimization

5.1. Optimization Objective

Vibrations induced by the interaction between the wheels and the rails are transferred from the rails to the invert and then propagate to the foundations of nearby buildings, resulting in structure vibrations and reradiated noise. Therefore, the vibration of the supporting platform reflects the structure vibration. To reduce the structure vibration, the platform vibration should first be reduced. Consequently, platform vibration acceleration was set as the optimization objective. The objective function waswhere is the vibration acceleration level of the supporting platform, is the tested number ranging from 1 Hz to 200 Hz, is the effective value of the vibration acceleration at the frequency point, and is the reference value of the vibration acceleration.

5.2. Optimization Variables

In engineering, the density and the elasticity modulus of the sleeper pad are always modified for the optimizing reduction in vibrations. Consequently, the density of the sleeper pad and the elasticity modulus of the sleeper pad were set as the optimization variables. In practical applications, the density of the sleeper pad ranges from 500 kg/m3 to 1500 kg/m3. The elasticity modulus of the sleeper pad varies from 2.6  ×  106 N/m2 to 7.8  ×  106 N/m2.

5.3. Optimization Results

The optimization results were obtained by the program presented in Section 4. The optimized density of the sleeper pad is 620 kg/m3 and the optimized elasticity modulus of the sleeper pad is 6.25  ×  106 N/m2. Before optimization, the vibration acceleration level of the middle supporting platform, ranging from 1 Hz to 200 Hz, was 120.8 dB. After optimization, the vibration acceleration level of the middle supporting platform was 94.4 dB. The vibration acceleration level decreased by 26.4 dB.

Figure 12 shows the frequency spectrum of the vibration acceleration of the middle supporting platform. Before optimization, the maximum vibration acceleration was 10.7  ×  10−4 m/s2. After optimization, the maximum vibration acceleration was 3.99  ×  10−4 m/s2. From Figure 12, the vibration acceleration was reduced across the frequency range after optimization.

6. Conclusion

To optimize the vibration reduction ability of ladder track, a new method, that is, the finite element method coupled with ant colony optimization, was proposed. We introduced how to build the FEM model and how to couple the FEM with ACO. The density of the sleeper pad and the elasticity modulus of the sleeper pad were optimized using this method. After optimization, the vibration acceleration level of the middle supporting platform was reduced from 102.8 dB to 94.4 dB. The optimized density of the sleeper pad was 620 kg/m3. The optimized elasticity modulus of the sleeper pad was 6.25  ×  106 N/m2.

Conflict of Interests

The authors declare that there is no conflict of interests with regard to this study and the publication of this paper.