A multi-GPU parallel computing method for 3D random vibration of train-track-soil dynamic interaction

Abstract

In this paper, an efficient computation method based on a multi-GPU parallel algorithm is proposed to overcome the low efficiency in random calculation of the train-track-soil coupled system (TTSCS). Firstly, for the large time consumption caused by solving multiple independent equations of TTSCS at different frequency points in serially random vibration analysis, the multi-GPU parallel algorithm is proposed and programmed based on the OpenMP-CUDA algorithm. The tasks of solving multiple linear equations for random vibration analysis of the TTSCS are distributed to different GPUs for parallel execution. On each GPU, the large sparse linear equations of TTSCS are solved by the CUDA-based parallel preconditioned conjugate gradient (PCG) method, and the large sparse matrix is stored in the compressed sparse row (CSR) format to save memory space. Then, the parallel computing program is implemented on the MATLAB-CUDA hybrid platform. Finally, numerical examples show that the efficiency of solving large sparse linear equations based on the multi-GPU parallel algorithm implemented on a 4-GPU node and the GPU-accelerated PCG algorithm implemented on a personal computer with a single GPU is 22.59 times and 3.75 times that of the multi-point synchronization algorithm (MPSA), respectively.

Keywords

random vibration / parallel computing / multi-GPU / three-dimensional finite element method / train-track-soil couple model

References

[1] OuakkaS, VerlindenO, KouroussisG. Railway ground vibration and mitigation measures: Benchmarking of best practices [J]. Railway Engineering Science, 2022, 30(1): 1-22

[2] ZhaiW-m, HeZ-x, SongX-lin. Prediction of high-speed train induced ground vibration based on train-track-ground system model [J]. Earthquake Engineering and Engineering Vibration, 2010, 9(4): 545-554

[3] PowrieW, YangL, ClaytonC R. Stress changes in the ground below ballasted railway track during train passage [J]. Proceedings of the Institution of Mechanical Engineers, Part F: Journal of Rail and Rapid Transit, 2007, 221(2): 247-262

[4] HallL. Simulations and analyses of train-induced ground vibrations in finite element models [J]. Soil Dynamics and Earthquake Engineering, 2003, 23(5): 403-413

[5] MandhaniyaP, ShahuJ T, ChandraS. An assessment of dynamic impact factors for ballasted track using finite element method and multivariate regression [J]. Journal of Vibration Engineering & Technologies, 2022, 10(7): 2609-2623

[6] MandhaniyaP, ShahuJ T, ChandraS. A parametric study of embedded slab track system for highspeed applications on cohesive subgrade [J]. Transportation Infrastructure Geotechnology, 2023, 10(1): 1-16

[7] AlvesC P, CalçadaR, SilvaC A. Track-ground vibrations induced by railway traffic: In-situ measurements and validation of a 2.5D FEM-BEM model [J]. Soil Dynamics and Earthquake Engineering, 2012, 32(1): 111-128

[8] AndersenL, NielsenS R K. Boundary element analysis of the steady-state response of an elastic half-space to a moving force on its surface [J]. Engineering Analysis with Boundary Elements, 2003, 27(1): 23-38

[9] KouroussisG, VerlindenO, ContiC. Free field vibrations caused by high-speed lines: Measurement and time domain simulation [J]. Soil Dynamics and Earthquake Engineering, 2011, 31(4): 692-707

[10] GalvínP, RomeroA, DomínguezJ. Fully three-dimensional analysis of high-speed train-track-soil-structure dynamic interaction [J]. Journal of Sound and Vibration, 2010, 329(24): 5147-5163

[11] MANDHANIYA P, SHAHU J T, CHANDRA S. Analysis of dynamic response of ballasted rail track under a moving load to determine the critical speed of motion [J]. Journal of Vibration Engineering & Technologies, 2022: 1 - 17. DOI: https://doi.org/10.1007/s42417-022-00741-3

[12] KouroussisG, ParysL V, ContiC, et al. . Using three-dimensional finite element analysis in time domain to model railway-induced ground vibrations [J]. Advances in Engineering Software, 2014, 70: 63-76

[13] ConnollyD P, KouroussisG, GiannopoulosA, et al. . Assessment of railway vibrations using an efficient scoping model [J]. Soil Dynamics and Earthquake Engineering, 2014, 58: 37-47

[14] XuQ-y, OuX, AuF T K, et al. . Effects of track irregularities on environmental vibration caused by underground railway [J]. European Journal of Mechanics-A, 2016, 59: 280-293

[15] SiL T, ZhaoY, ZhangY H, et al. . Random vibration of an elastic half-space subjected to a moving stochastic load [J]. Computers & Structures, 2016, 168: 92-105

[16] ZengZ-p, LiuF-s, LouP, et al. . Formulation of three-dimensional equations of motion for train-slab trackbridge interaction system and its application to random vibration analysis [J]. Applied Mathematical Modelling, 2016, 40(11–12): 5891-5929

[17] WangL-d, ZhuZ-h, BaiY, et al. . A fast random method for three-dimensional analysis of train-track-soil dynamic interaction [J]. Soil Dynamics and Earthquake Engineering, 2018, 115: 252-262

[18] WestphalE, SinghS P, HuangC C, et al. . Multiparticle collision dynamics: GPU accelerated particle-based mesoscale hydrodynamic simulations [J]. Computer Physics Communications, 2014, 185(2): 495-503

[19] ChenX, JieY-x, YuY-zhen. GPU-accelerated iterative solutions for finite element analysis of soil-structure interaction problems [J]. Computational Geosciences, 2013, 17(4): 723-738

[20] LinS-z, XieZ-qiang. A Jacobi_PCG solver for sparse linear systems on multi-GPU cluster [J]. The Journal of Supercomputing, 2017, 73(1): 433-454

[21] ZhangFan. Research on parallel computing performance visualization based on MPI [C]. 2010 2nd International Conference on Advanced Computer Control, 2010, Shenyang, China, IEEE: 323327

[22] ZhangY-j, YuanJ-b, LuX-w, et al. . Multi-GPU parallel computing and task scheduling under virtualization [J]. International Journal of Hybrid Information Technology, 2015, 8(7): 253-266

[23] HouZ-l, HuangD-nian. Multi-GPU parallel algorithm design and analysis for improved inversion of probability tomography with gravity gradiometry data [J]. Journal of Applied Geophysics, 2017, 144: 18-27

[24] LuF, LinJ H, KennedyD, et al. . An algorithm to study non-stationary random vibrations of vehicle-bridge systems [J]. Computers & Structures, 2009, 87(3–4): 177-185

[25] KrishnasamyE, SourouriM, CaiXing. Multi-GPU implementations of parallel 3D sweeping algorithms with application to geological folding [J]. Procedia Computer Science, 2015, 51: 1494-1503

[26] KoricS, GuptaA. Sparse matrix factorization in the implicit finite element method on petascale architecture [J]. Computer Methods in Applied Mechanics and Engineering, 2016, 302: 281-292

[27] WangS, YanX-k, ZhangY, et al. . Research on EBE-FEM realized by CUDA applying to electromagnetic field analysis [C]. 2018 IEEE Student Conference on Electric Machines and Systems, 2019, Huzhou, China, IEEE: 14

[28] VuT, WangY-h, XiaJun. Optimizing photoacoustic image reconstruction using cross-platform parallel computation [J]. Visual Computing for Industry, Biomedicine, and Art, 2018, 112