1. Introduction
Railway transportation is regarded as the most sustainable means of moving passengers and cargo in modern society. The advent of high-speed railways (HSRs) has transformed the landscape of railway travel in recent decades, with freight logistics increasingly integrating rail as a core component. This shift is driven by environmental concerns and economic efficiencies. However, railway systems, despite their numerous advantages, generate ground-borne vibrations that can have significant impacts on their surrounding environment. These vibrations originate from the wheel–rail contact, and propagate through track components, subgrade, and soils, affecting surrounding nearby structures [
1].
The sources of vibration effects include quasi-static rolling stock moving loads, deformations induced by wheel and/or rail unevenness, punctual track defects, axles passing by, and other track components. To predict vibrations accurately, it is crucial to define the transmission medium between the source and receiver. In addition, a comprehensive analysis of wave propagation and potential attenuation involves examining not only the time domain but also the frequency domain. Geometrical attenuation, influenced by the distance between the source and measurement point, is unrelated to the frequency domain. On the other hand, material typology attenuation, influenced by the properties of traversed materials, is directly associated with the frequency response [
2]. Nevertheless, predicting and assessing these ground-borne vibrations is paramount, especially when planning new railway lines or expansions, as these vibrations can profoundly impact human well-being, particularly within the frequency range of 15 to 200 Hz [
3].
Two predominant families of prediction models exist, distinguished by their underlying calculation methods: analytical and numerical. These models diverge in terms of the mechanisms they incorporate to simulate excitation (e.g., quasi-static axle loads, parametric excitation, transient excitation due to rail joints and wheel imperfections, and excitation stemming from wheel and rail irregularities) and their approach to modelling dynamic structure–soil interactions [
4]. In this work, a range of relevant studies were analyzed in order to build the base of the state-of-the-art of latest-generation numerical prediction models for ground-borne vibrations.
Initial findings and studies of prediction models assumed homogeneity in the track direction for the ground and structures. By combining Boundary Element Method (BEM) and Finite Element Method (FEM) models to reduce computational efforts, some researchers sought to further enhance the analysis [
5]. Other researchers have introduced an all-encompassing 3D Boundary Element Method (BEM) model designed to assess soil motion and the effects of high-speed train (HST) passage on both surface and underground structures, facilitating complete coupling between the soil and adjacent structures. Nevertheless, it is important to highlight that this method heavily depends on environmental factors, leading to an extensive simplification of track elements. This oversimplification may lead to a lack of detailed information concerning events occurring within the track [
6].
These kinds of vibrations are intricately linked to various operational factors, such as the condition of the infrastructure, the train speed, and the use of adjacent areas. Primarily, ground-borne vibrations stem from the metallic contact between the train wheels and rail. Critical to numerical predictions in this domain is the selection of an appropriate method to represent wheel and rail unevenness. Various methods have been explored in relevant literature to incorporate them as an additional source of realistic vibrations in prediction models [
7,
8]. Without these considerations, the resulting signals may lack fidelity.
After conducting this analysis, the decision was made to adopt the standardized representation [
9] from various options presented in the literature for addressing such irregularities. This approach was deemed the optimal choice for ensuring reproducibility. In addition to developing the track model, careful attention must be given to defining the representation of the rolling stock. It is essential not only to create a representative model but also to ensure its stability in calculations and realism by incorporating genuine characteristics of the rolling stock. Having identified references that provide detailed descriptions of a 3D rigid–flexible multibody system [
10], the modeling process can benefit from a thorough understanding and incorporation of the nuances involved.
Load simulations are generally modeled as a birth-and-death system that includes loads that are punctual and time-dependent [
11]. Models must be designed based on the calculation time step to accommodate loads at existing nodes. Loads are distributed as triangular pulses between three nodes. As loads are modeled and vibrations propagate through the tracks and into the ground, they inevitably affect the foundations of nearby structures. Typically, the intensity of ground vibrations diminishes as one moves further from the railway track [
12].
There are many types of vibration-prediction models and load simulations, depending on the degree of vehicle modeling and the desired level of representation of the track and environment. In [
13], a schematic and detailed summary of these types is presented, with the references cited therein analyzed, identifying those of interest for the detailed modeling of rolling stock and the 3D representation of the track and environment using FEM. It is precisely in this latter case where the fewest reference works are found. Other researchers [
14] employ computational contact techniques at the core of their investigation to represent wheel–rail interactions, including unevenness. The main difference is the consideration of simplifications by assuming linear calculations, which assumes the invariance of the system matrices. Additionally, another approach introduces an explicit modelling method [
15], with the introduction of point forces through the use of springs.
All these numerical models were developed and published to serve as the foundation for the present discussion and to highlight the potential implications of these advancements in the field of railway track engineering for effectively controlling vibration emissions. The primary objective of this paper is to introduce and validate a numerical Finite Element (FE) model based on a 3D spatial framework. This model has the versatility to accommodate a wide range of variations along the railway track, thereby enabling a comprehensive understanding of the vibration dynamics. It is important to note that the model will focus on the analysis of vibrations in the immediate proximity of the source, the rail. The goal is to precisely characterize the vibrational signal across the entire track system, allowing for the identification and potential redesign of specific segments with attenuation objectives as needed.
The proposed model incorporates computational contact tools to simulate the transmission of waves throughout the entire train-track system seamlessly, effectively eliminating discontinuities in the application of loads [
2,
11]. Working with a 3D model and a small time step of 1024 calculations per elapsed second (emulating field measurements with piezoelectric accelerometers), as well as the use of computational contact theories, allows for the continuous introduction of solicitations and the acquisition of a sufficiently dense FFT spectrum up to 512 Hz. This advancement not only enhances time-domain accuracy but also improves the frequency-domain data that can be predicted.
Furthermore, the main novelty lies in simulating a continuous model from the rolling stock to the track. This enables the representation of complex dynamic phenomena, such as those arising from the connections between the carriages of the rolling stock in their complete longitudinal composition and their kinematic motion linkages, which could not be simulated using decoupled models like load-overlay models. By employing this innovative approach, this research aims to contribute to the field of railway engineering by offering a robust framework for assessing and mitigating ground-borne vibrations. This framework is particularly valuable during the infrastructure design and planning stages, as it empowers decision makers to proactively address the potential impact of railway-induced vibrations on sustainability and the quality of life for nearby communities.
2. Experimental Data
The genesis of ground-borne vibrational prediction models stemmed from the imperative need for a foundational dataset, serving as a reference benchmark for model development. This underscores the necessity for any subsequent model to have an initial reference point or a reliable comparison framework for calibration and result validation. At this juncture, a fundamental reference for the calibration and validation of such vibrational models can be found in [
4]. The study presents measurements from an actual railway track, using accelerometers to monitor a high-speed rail (HSR) line with a double-track configuration, considering the rail, sleeper, ballast, and terrain. In this model, the scenario of a train traveling at 300 km/h was specifically chosen as the most representative of a high-speed rail (HSR) line.
In the particular case presented in [
4], the Brussels–Paris HSR line was monitored as a conventional ballast track, as illustrated in
Figure 1. The layered structure of the track and subgrade system consists of the following components:
Subgrade: The foundation layer of the railway system, representing the ground beneath the track.
Subballast: Positioned above the subgrade, this layer acts as an intermediary component between the subgrade and ballast.
Ballast: The layer in which the track is embedded, providing essential support and vibration dampening.
Sleepers: Components strategically placed along the track, supporting the rails and distributing loads.
Under-Rail Pad: A crucial element for vibration isolation, positioned between the rail and sleeper.
Rail: The track element that interacts directly with the rolling stock.
Each of these components is defined by a series of inherent characteristics, specific properties that describe the behavior of each material and its response to various conditions, ultimately influencing the transmission of ground-borne vibrations generated by the passage of rolling stock over them. Subsequently, the paper outlines the key characteristics of each component used in [
4]. In the proposed study, this assumption is crucial for developing a finite element model that accurately simulates the actual behavior of the materials and components related to the high-speed rail (HSR) line studied in the previously cited research.
The track features continuously welded UIC 60 rails, each with a mass-per-unit length of 60 kg/m and a moment of inertia (I) of 3038 × 10−6 m4. These rails are affixed using a Pandrol E2039 rail fixing system onto precast, prestressed concrete monoblock sleepers measuring 2.5 m (m) in length, 0.285 m in width, and 0.205 m in height (under the rail), with a mass of 300 kg.
Flexible rail pads with a thickness (t) of 0.01 m and a static stiffness of approximately 100 MN/m, suitable for loads ranging from 15 to 90 kN, are positioned under the rail. The track is supported by a 30 cm layer of ballast, followed by a layer of limestone or porphyry (subballast), and finally, a limestone supporting layer (subgrade).
In addition, the measurements obtained from [
4] were conducted during the passage of a Thalys–Alstom high-speed train (Spanish Renfe S100), which in this study serves as the model load.
Figure 2 displays the arrangement of the Thalys high-speed train (HST), comprising two locomotives and eight carriages, with a total train length of 200.18 m. Each locomotive is supported by two bogies and is equipped with four axles. In proximity to the locomotives, adjacent carriages share a bogie, while the remaining six carriages share bogies with neighboring ones. This brings the total number of bogies to 13, resulting in a train with a total of 26 axles.
Moreover,
Table 1 provides a summary of key parameters, including the axle mass (
Mu), bogie mass (
Mb), car mass (
Mc), primary and secondary stiffnesses (
k1 and
k2), and primary and secondary damping coefficients (
C1 and
C2) obtained from [
16].
This initial stage of our investigation drew upon data obtained from the literature, encompassing the current state of numerical models and actual measurements from the track. The subsequent stages of the process involved the development, calibration, and validation of the model against on-site measurements using the comparison method. This method primarily entails the graphical comparison of the obtained results with the measured ones for assessment. This calibration process is essential to ensure the accuracy and reliability of our model in replicating real-world ground-borne vibrations.
In the subsequent sections, the developed numerical model is revealed in detail, and the implications of this research are analyzed within the context of railway engineering.
3. Materials and Methods: Numerical Model
The proposed model is a Finite Element Method (FEM) model of a layered approach (presented in
Figure 1) to represent the subgrade and track system accurately in a 3D environment, as shown in
Figure 3. The FEM model generated in this work represents and incorporates all of the components described in
Section 2, defining these components with the actual characteristics and properties of the railway infrastructure and superstructure to be modelled.
To perform the dynamic analysis of the railway system, a general multibody analysis approach is employed, in which two important points need to be considered:
The equations of motion (described in
Section 3.3) govern the translational and rotational motions of the rolling stock and track elements. These equations consider the forces generated by the train’s weight distribution, wheel–rail interactions, and the response of the track components to external loads.
The discrete equations of the system (both flexible and rigid parts), corresponding to the conservation of linear momentum, are obtained following the standard displacement approach in FEM. These equations are completed by adding the constraint equations associated with the joints that connect different bodies.
The Newmark method (described in
Section 3.4) is the time-integration method used because of its simplicity and proven capacity to capture transient dynamic responses accurately.
By incorporating these mathematical formulations and computational techniques, the model can provide a comprehensive representation of the complex dynamics involved in ground-borne vibrations induced by railway operations. In the subsequent subsections, a detailed explanation of all the subsystems that compose the entire model and their fundamental equations is described.
3.1. Railway System and Terrain
The spatial framework of our model is meticulously designed to simulate distinct layers that capture the behavior of the railway track and its surroundings as shown in
Figure 3. The structural components of the railway system and the adjoining subgrade undergo meticulous representation through the utilization of solid continuum elements, specifically employing eight-node hexahedron elements. Given the assumed constitutive behavior of linearity, isotropy, and compressibility, a pure displacement formulation is employed to accurately capture the dynamic response of the system. Contact interactions between different parts are modelled using surface-to-surface contact by target and contact elements. Geometric contact detection is performed over target elements, and the role of the contact elements is to provide the points of contact.
The contact algorithm used is the augmented Lagrangian method [
17]. Frictionless contact with sliding and separation is considered between the sleepers and ballast, bonded contact between the under-rail pads and sleepers and rail, and separation without sliding contact between the rolling stock wheels and rail.
Wheel to rail: To simulate the interaction between the rolling stock and rail, non-glued computational contacts are programmed. These contacts allow for the dynamic interaction between the train and track, accounting for factors such as wheel–rail interaction forces and rail irregularities.
Rail to sleeper: The connection between the rail and sleeper is modelled using a “glued contact”. This modelling technique replicates the effect of the fastening system, ensuring that the rail and sleeper act as a cohesive unit in response to dynamic forces.
Infinite media: The ground is modelled as an infinite media. To prevent wave reflections that could distort the results, this infinite medium is designed to be wave-absorbent. The simulation employs both linear and viscous damping, following the approach outlined in [
18], where damping elements introduce the shear and normal components of ground forces along its boundary. These calculations are based on ground data inputs, specifically the velocities of P-waves and S-waves, which are dependent on the mechanical properties of the soil.
The properties of the track and soil are detailed in
Table 2.
Figure 4 illustrates a side view of the contacts between wheel and rail and rail and sleeper through the under-rail pad.
3.2. Rolling Stock
The rolling stock is modelled as a rigid multibody system. Each rigid part is a compound of a pilot node, which drives the motions of the rest of nodes of the rigid body. Bodies are connected to each other through joints, imposing kinematic constraints between the two nodes they connect. Additionally, spring-damper elements are introduced between bodies to model the suspension system, as shown in
Figure 5.
The approach used to define rigid bodies involves a master node that controls the motion of all other nodes within the body, as shown in
Figure 5. The rolling stock model is represented using rigid elements connected with masses and springs, reflecting the physical characteristics of the train. Two types of joints are employed: translational and spherical. The translational joints impose five degrees of freedom (DOFs) between the nodes—two for translations and three for rotations—while the rotational joints impose three relative rotations between the nodes. These constraints are implemented using the Lagrange multiplier method [
19].
A key feature of the model is its ability to incorporate rail irregularities, which significantly impact ground-borne vibrations. Rail irregularities include various imperfections and deviations along the rail profile. It is important to note that the model simplifies the wheel’s behavior by excluding rotational movement, focusing instead on translational motion.
Additionally, the method proposed for modeling the continuous introduction of loads and wheel–rail irregularities in the track model involves constructing three types of nodes: one for introducing wheel–rail irregularities and two for simulating the wheel–rail contact. These nodes are as follows (see
Figure 6):
Physical node: Represents a characteristic point of the wheel.
Joint node: Defines the translational joint element where wheel–rail irregularities are applied, in conjunction with the physical node of the wheel.
Master node: Defines the element through which contact between the rolling stock and the track model is introduced, in conjunction with the physical node of the wheel.
Figure 6 shows the previously described nodes.
In summary, irregularities and interaction of contact are modelled as follows:
Modelling of wheel and rail irregularities: Achieved by incorporating a translational joint that enforces a relative displacement between the wheel and rail nodes, introducing the corresponding irregularity at each instance. This is the element between the joint node and the physical node.
Modelling wheel–rail contact interaction: Accomplished using computational node-segment contact. This contact allows all forces or vibrations generated by the rolling stock to be transmitted to the track system, represented by the element between the physical node and the master node.
A wheel–rail roughness spectrum, which combines data from both the wheel and rail, can be used to describe the statistically consistent excitation of a specific vehicle and track. This allows for the prediction of the resulting power spectral density (PSD) of the response using random process theory [
20].
The method chosen for reproducing wheel and rail irregularities, as detailed by [
9], is implemented at this stage.
Figure 7 shows the deflection values versus space for these surfaces. These elements are integrated into the overall mechanical system described, serving as an additional factor to consider.
The interaction between the train wheels and the rail is simulated using point-to-surface contact, incorporating a frictionless interface with capabilities for both sliding and separation. The contact algorithm used in this simulation is based on the augmented Lagrangian method [
17].
A key aspect of this model is the use of computational contacts to accurately represent the connections within the railway system. Unlike the previously common approach of segmenting loads, which results in discontinuous representations, this method introduces a significant improvement by converting these discontinuous load representations into continuous ones. This marks a notable shift in methodology.
3.3. Governing Equations
To express the dynamic behavior of the railway system within the model, fundamental equations of motion and dynamic analysis techniques are employed. The system of discrete equations for a multibody system subjected to holonomic constraints is as follows [
19]:
is the mass matrix, which is a function of the displacements because of the presence of rigid bodies.
is a vector, which include the effects of external forces (gravity) and internal forces from dampers.
and
are the constraint matrix and vector corresponding to the constraint equations introduced by the joints and the contact algorithm, and
is the Lagrange multipliers vector. Linearization of the system of equations (Equations (3) and (4)) gives [
19] the following:
where:
3.4. Loading Description
The simulation is divided into two steps. In the first step (initial stress introduction), inertia effects are not considered. This phase focuses on introducing initial stress due to the self-weight of both the terrain and rolling stock. The rolling stock is assumed to be on a rigid platform, so the terrain’s stress state is solely due to gravity. This step involves a nonlinear static analysis, where geometric nonlinearity accounts for the initial deformation of the rolling stock caused by its self-weight, ensuring consistency with equilibrium stresses. In the second step (transient analysis), inertia effects are included, and a full nonlinear transient analysis is performed. This step captures the dynamic behavior of the system, considering the effects of inertia on the rolling stock and terrain.
Both steps are solved as part of the same transient analysis, with the first step focusing on static effects and neglecting dynamic effects. In this initial step, displacement and velocity are set to zero to establish initial conditions. During the second step, dynamic effects are introduced, and velocity is applied to the master node of the traction vehicle.
The time-integration method used in the second step is the Newmark method [
21]. This method assumes that:
where the super index
indicates a value at the end of a time increment, and the superindex
t indicates a value at the beginning of the time increment. Time integration parameters of the Newmark method,
and
, are expressed in terms of the amplitude decay factor
, which is introduced to filter the influence of high frequencies:
In the present simulations, due to the high-speed vibrations, the value of used for this model is 0.005.
3.5. Boundary Conditions
To ensure realistic simulations, our model adheres to specific boundary conditions:
Additionally, the model incorporates a critical feature proposed by [
18]. This involves placing several springs around the platform at the vertical planes of the model’s beginning and end, as well as beneath the lower horizontal plane, to ensure the effective absorption of wave reflections.
By incorporating these specific boundary conditions and absorption mechanisms, the developed model achieves a robust representation of railway-induced ground vibrations. This setup enables the simulation of realistic scenarios while effectively mitigating boundary effects that might otherwise compromise the accuracy of the results. The subsequent section presents the findings from these numerical simulations.
4. Results and Discussion
In this section, the results obtained from the numerical model are presented, which are a valuable tool for quantifying its environmental impact and elucidating the complex dynamics of ground-borne vibrations in railway systems.
The results detail the effects observed when a train passes a specified reference point.
Figure 8 illustrates the time-domain accelerations recorded along the
y-axis in the rail at a specific location—the precise center of the entire model.
It is crucial to highlight that, in the displayed graph, a clear correspondence is evident between the accelerogram peaks (passing of the axles) and the central areas among them, emphasizing the geometry and load distribution of the rolling stock. This observation signifies the visible representation of the energy generated as the rolling stock passes.
The central regions between the peaks illustrate wheel and rail irregularities, as well as other movements caused by the train’s motion. The smaller intermediate peaks, characterized by less energetic vibrations, serve as a representation of the enrichment introduced by the continuous modelling of the entire system. This encompasses the implementation of a continuous system and continuous loads. As previously discussed in the literature review, the information presented here must be complementarily showcased in both the time and frequency domains to provide a comprehensive understanding of the situation.
These nuanced results, complementing those obtained in the time domain, are illustrated in
Figure 9. It is noteworthy that the entire spectrum is fully captured, capturing the most representative frequencies—peaks, attributed to the reduction of the time step and the continuous way of introducing the loads, with calculations fixed at a frequency of 1024 Hz—which are identical to the referenced measured vibration spectra. This translates to an average of over 4000 full-model calculations for every four seconds of simulation (1024 calculations per second), covering the entire passage of the train.
Regarding the representation of results in the frequency domain, a comprehensive frequency representation is observed across the entire spectrum to have a complex image of the predicted phenomena. The model’s groundbreaking advancements become evident in its representation of results within the frequency domain. It offers a comprehensive view across the entire spectrum, providing a nuanced understanding of predicted phenomena. These insights significantly enhance our grasp of vibrations stemming from railway operations, particularly when integrated into three-dimensional numerical models capable of accommodating various track system configurations.
The wealth of information provided by this model is indispensable for supplementing data in the time domain, especially when considering the design of specialized track components or attenuation measures to mitigate vibrational pollution. Notably, the ability to identify prominent or energetic vibrational frequencies facilitates an informed material selection for optimal effectiveness.
In essence, this model represents a significant leap forward, offering unparalleled insights that can revolutionize our approach to mitigating vibrational impacts in railway environments.
5. Experimental Validation
Thanks to published and discussed data [
4], this model underwent validation through a comparison with a scenario similar to the conditions mentioned in the experimental data section. The validation of the model’s results was carried out through direct comparison. This method of assessing results was implemented due to its widespread use in this type of predictive models by highly relevant authors in the field [
8,
23,
24,
25,
26,
27,
28,
29].
The objective of these predictions is not to acquire precise data in absolute terms relative to the initial data series, but rather to conform to the characteristics of the power spectral density (PSD) produced by the phenomenon of the rolling stock traversing the track. This phenomenon is inherently difficult to precisely adjust, as its most influential effect is the wheel–rail irregularity, which comprises a unique and variable signature at each spatial position of both the wheel and the rail, as irregularities are introduced statistically [
20]. On-site measured data from the referenced study are depicted in
Figure 10 and
Figure 11, representing the time and frequency domains, respectively.
These presented results conclusively demonstrate that the main graph peaks share the same order of magnitude, and the passage of the axles is consistently depicted, particularly in the time domain (
Figure 8 vs.
Figure 10). This signifies that the simulated values align closely with the experimental on-site measurements, establishing the presented model as validated. Shifting to the frequency domain, the entire spectrum is encompassed by the assumed standardized wheel–rail unevenness, which is generic. This successful coverage affirms the correctness and effectiveness of the employed method in this domain.
Additionally,
Figure 12 is presented below to provide a comparison of the data in both the time and frequency domains, offering a more robust validation of the model’s performance.
6. Conclusions
These deductions are obtained from rigorous simulations and meticulous data analyses obtained from the analyzed bibliography. These findings shed light on the complex interplay between the rolling stock, track components, and surrounding environment, ultimately enhancing our understanding of the factors contributing to ground-borne vibrations. Through these results, we embark on a journey to untangle the intricacies of railway-induced vibrations, paving the way for informed decision making and effective mitigation strategies.
The main conclusions of this article could be grouped by the following statements:
Continuity in load representation has yielded substantial benefits, particularly in the faithful depiction of frequency characteristics, as shown in the Results section.
Advancements in frequency domain accuracy: A pivotal outcome of our research is the substantial enhancement of numerical models, particularly in the frequency domain.
The accurate representation of measured data signifies a significant advancement in railway engineering, facilitating more dependable assessments and predictions of ground-borne vibrations caused by train operations.
This approach facilitates a thorough depiction of phenomena from the source through the entire track system, enabling the effective simulation of various elements’ influence so as to attenuate these vibrations.
The model’s proficiency in capturing frequency-specific data offers a profound understanding of the vibrational characteristics linked to railway-induced ground vibrations.
This model is capable of quantifying environmental effects, which is essential to allow the precise planning and assessment of future railway infrastructures, their vibrational emissions, and their overall sustainability.
A possible limitation of this model mainly stems from its dependency on rail irregularities, which have been identified as the most influential factor in the frequency domain.
In summary, transitioning from discontinuous to continuous load simulations, our model enhances the precision of frequency representations, providing a more comprehensive and accurate assessment of ground-borne vibrations induced by railway operations. The advancements made in the presented numerical model and the seamless integration between time and frequency domains stand as a testament to the strides taken toward more accurate, sustainable, and resilient railway infrastructures.
7. Further Developments
These findings hold immense promise for the future planning and design of railway systems, paving the way for more efficient and environmentally conscious transportation networks. Some potential areas for further development could include:
Utilization of advanced models in railway superstructure design: The utilization of these advanced numerical models holds promise in designing future railway superstructures, leading to more robust and vibration-resilient railway infrastructures.
Improved sustainability assessment and design metrics: This enhanced index holds the potential to serve as a cornerstone for evaluating sustainability aspects in railway construction, enabling better-informed decisions and more sustainable infrastructure designs.
A possibility to study the influence of the quality of track maintenance and rolling stock wheels, with a focus on different scenarios using actual measured rail irregularity data, preferably obtained directly from track measurements, to further enhance the model’s accuracy and applicability.