This section describes the proposed MS, including all the theoretical models used. As previously explained, the MS proposed here was developed to allow the simulation of the layered elements of a railway infrastructure (e.g., ballast layer and subgrade layer), and a moving railway vehicle, in the context of an MBS model. More specifically, the main goal was to be able to obtain realistic results from the vehicle–track dynamic interaction, with a focus on the vertical resonance frequency of the underlying soil layer of the substructure (i.e., the first subgrade layer).
2.1. Underlying Premise and Main Concept
Before describing the main concept behind the proposed MS, this subsection presents an important premise permitting the underlying logic that was used in its design, in terms of physical and mechanical principles. The premise in question is: “The values of the natural frequencies of the soil layers beneath a railway track do not directly depend on the presence of a railway vehicle on said track”. This means that the observable values of the natural frequencies of the soil layers do not directly change due to the presence of a vehicle above those layers, and do not change based on the modal parameters of the vehicle. This statement does not disregard the expected smaller changes in the values of the resonance frequencies from the soil layers due to the dynamic interaction between the two mechanical systems in question (i.e., the railway infrastructure and vehicle). Instead, it just states that there is no direct change in these natural frequencies caused by the presence of a railway vehicle, since the characteristic frequencies of a linear mechanical system are independent from the properties of the excitation source used to assess them [
18]. This concept is applicable here as an approximation of this complex multi-element system, for the purpose of modal analysis, in the same way that it is acceptable in the context of analyzing other common engineering structures [
20]. This premise has been clearly presented here first, due to its relevance to the topic under study, since this is something that is not commonly stated in the bibliography that deals with the natural frequencies of soil layers in the railway context [
4,
19,
21].
In an MBS model, using only rigid-type body elements, there is an inherent detachment between the mechanical properties of a body element and its geometric dimensions. While overall mass can be established using a body’s geometry and its density, parameters related with other mechanical properties cannot be directly linked to body elements (e.g., Young’s modulus). These other mechanical properties are usually simulated through other element types, such as springs and dampers. In the present context, this means that there is a disconnection between the geometry/mass of a subgrade layer body element and its equivalent stiffness characteristics. Furthermore, this means that the natural frequency values of the simulated subgrade layer are dependent on the combination of the imposed mass of the subgrade body element and the stiffness of the associated spring element, which have no intrinsic relationship. This is in opposition to other modeling software (e.g., FEM model), where the mechanical properties of a subgrade layer can be intrinsically established during model setup. Thus, in an MBS model, the subgrade’s mechanical properties are not fully linked to enable an adequate simulation of the subgrade’s natural frequencies. Another aspect related to this issue is the difficulty in obtaining realistic values for the modal properties of a soil layer, as opposed to a common steel beam, which hinders the option of simply imposing these values on an MBS model.
The proposed MS overcomes these issues by providing a set of theoretical models that result in a fully determined system for the subgrade layer element’s geometry to be used in an MBS model. More specifically, by selecting a cylindrical shape for the body elements of the railway’s layered elements (e.g., a subgrade layer and a ballast layer) and some input parameters related with their mechanical properties, this MS provides an equation that gives the diameter of these cylindrical elements as a function of subgrade layer height. As previously explained, the proposed MS represents a potential innovation in the specific field of soil layer modeling in the railway context, which is denoted by both the combined model and the geometry shape choice, since the authors did not find any similar solution in the consulted bibliography. As will be presented in
Section 2.2, soil layer height is a major parameter in the process of estimating soil layer natural frequencies using the established theoretical model [
12,
19], hence its importance in this combined model.
2.2. Proposed Modeling Solution
The MS proposed here is a combination of four distinct theoretical models, which together create a deterministic solution for the layered elements’ geometry in an MBS model.
Table 1 presents a list of the main geometrical and mechanical parameters that relate to these theoretical models. The control parameter of this combined model, to verify the realistic scenario conditions, is the natural frequency of the first vertical vibration mode of the subgrade layer. The four theoretical submodels (SM) are described below.
The first submodel (SM-Freq) is the theoretical equation of the natural frequencies of a surface-level soil layer from a natural soil deposit [
12,
19]. For each vibration wave type (e.g., longitudinal waves or shear waves), the general form of this equation can estimate the frequency values for the
n vibration modes of a soil layer. The kinematics of each vibration mode depend on the type of vibration wave created by the excitation source and its position relative to the soil layer. The specific form of this equation used here (Equation (1)) gives the frequency value of the first vibration mode (
f1), due to a vertical excitation source, based on the soil’s characteristic longitudinal wave velocity (
cp) and the soil layer’s height (
h). This corresponds to the first vertical vibration mode that is dynamically excited by the presence of an excitation source above the soil layer, such as a railway vehicle, which is the most relevant situation in the present context [
12].
The second submodel (SM-Wave) is the equation for the characteristic longitudinal wave velocity of a soil layer (
cp) [
12,
19]. This equation (Equation (2)) provides a relationship between relevant mechanical properties of the soil to be simulated, namely, its shear modulus (
G), its density (
ρ), its Young’s modulus (
E) and its Poisson’s ratio (
ν). Equations (3) and (4) present the formula for the Lamé parameters (
λ,
μ), which are relevant for Equation (2). SM-Freq and SM-Wave represent a deterministic pair of equations for the natural frequency of a soil layer, for a given set of the relevant input parameters. This theoretical formulation is well known in the geotechnical field [
21]. However, it is not enough to define the geometric shape of the railway’s layered elements in an MBS model, since it only establishes layer height. Therefore, there was a need for additional models to be added to the overall formulation.
The third sub-model (SM-Lumped) is a lumped-elements type model that can provide an adequate estimate of the subgrade’s resonance frequency of interest obtainable from an MBS simulation [
17]. This theoretical model was developed in the context of the new TCMM that was previously mentioned. It is a two Degree-of-Freedom (2-DoF) model, containing an assembly of two spring-mass subsystems in series, which represents the vehicle, per axle, and the railway infrastructure beneath it. The lower mass element represents the unsprung mass of one vehicle axle and the modal mass of the infrastructure section that is being dynamically excited by the passing vehicle (
m1), and the second represents the vehicle’s sprung mass that is supported by one axle (
m2). These body elements are linked by a spring that represents the vehicle’s suspension system (
k2). The remaining element is a spring that characterizes the modal stiffness of the railway infrastructure (
k1). This last element aggregates all stiffness contributions from the entire railway infrastructure, including from the subgrade layer. The outputs of this model are two natural frequencies (
ω1 and
ω2) [
20]. Frequency
ω2 is mainly influenced by the vehicle’s modal characteristics and corresponds to the vertical oscillatory motion of the vehicle’s sprung mass, which has a frequency range between 1 and 5 Hz for common passenger carriages [
10]. On the other hand, frequency
ω1 is influenced by the modal parameters of the railway infrastructure elements. This is the frequency that represents an adequate estimate of the subgrade’s first vertical resonance frequency obtainable from an MBS model (
fsim) in the presented conditions [
17]. Morais et al. (2021) provides a more in-depth description of this TCMM and this concept [
16].
The formula for frequency
ω1 is the equation associated with SM-Lumped (Equation (5)). In this equation,
k2 represents the vertical stiffness of the primary suspension of the vehicle, even for a vehicle with two suspension levels, since it is the spring element that separates the sprung and unsprung mass elements of the vehicle. Equation (5) also shows how
ω1 is calculated with the relevant parameters from
Table 1.
Having integrated the layered elements’ mass (included in m1) in the overall formulation, the conditions to establish their geometry were still not met. The system of equations was still undetermined due to the vertical stiffness term of the subgrade layer (ki.su).
The fourth sub-model (SM-Stiffness) is an equation that gives an estimate of the vertical stiffness that a platform or foundation would perceive when placed on a single soil layer (
kf). This sub-model (Equation (6)) defines the soil layer as a homogeneous half-space and is solely based on the in-plane dimensions of the platform and the mechanical properties of the soil layer [
22]. This equation will be used to provide an estimate for the subgrade layer vertical stiffness (
ki.su), a concept that has already been used in another railway-related context to study vibration propagation from a railway track [
12]. In the railway context, the platform’s dimensions used are the length (2
a) and width (2
b) of a sleeper. The soil’s mechanical properties relevant for this equation are its shear modulus (
G) and its Poisson’s ratio (ν).
With this additional equation, the combined model now represents a deterministic solution for the geometry of the railway’s layered elements to be used in an MBS model. However, before describing in detail how this combined model is implemented, the following list presents some relevant aspects that are related to the required modeling conditions:
This formulation works with a vehicle control volume of one axle and the respective portion of vehicle sprung mass (e.g., for a vehicle with four axles, the sprung mass to be considered here would be ¼ of the total mass);
Only a single sleeper is being considered as part of each set of infrastructure elements;
A cylindrical shape is being used for all the railway layered elements to be modeled in the MBS software, and their respective heights are imposed by the user;
The equivalent vertical stiffness and density of the railway layered elements are also selected by the user. From a practical point of view, estimates of these vertical stiffness values can be obtained based on experimental data, or obtained from theoretical track deformation values and the respective imposed testing load.
Finally, this combined model is implemented by pairing SM-Freq with SM-Lumped, via imposing equality between the frequency values obtainable from these two submodels (
f1 and
ω1). This results in an equation that can be solved in respect to the diameter of the layered elements of the railway infrastructure (
θ). Its final form is shown in Equation (7), which makes use of the relevant parameters from
Table 1. This equation represents the main output of this combined model.
Of note, while the modeling scenario considered here involves a railway vehicle placed above a railway infrastructure, this formulation should still work in other scenarios with just some adjustments. For example, to represent only a soil layer and no vehicle would only require an adjustment to the lumped elements model used in SM-Lumped by considering only a 1-DoF model (
m1 and
k1). Naturally, this MS had to be tested to verify the validity of the concept. Some results regarding this topic are shown in
Section 3.
2.3. Simpack® Numerical Model
This subsection presents the Simpack
® (version 2021x) model (
Figure 1) that was developed based on the proposed MS to assess its validity and perform a sensitivity analysis on the parameters that influence Equation (7). This will be checked by comparing the subgrade’s theoretical resonance frequency value given by the combined model with the equivalent resonance frequency obtained from the Simpack
® model (
fsim). The sensitivity analysis results for some of these parameters will also be presented here, along with a discussion on the overall validity of this concept. The core elements of this MBS model were a basic railway passenger carriage, with two bogies, and four sets of elements to represent the railway infrastructure presence below each vehicle axle. The geometry and dimensions of the layered components of these railway infrastructure element sets were defined based on the proposed MS by selecting subgrade layer heights from a predefined testing range.
The generic vehicle model included all the elements relevant for these tests, namely, springs and dampers elements for the two levels of suspension systems, and mass elements to represent both the sprung and unsprung mass components of the vehicle. A passenger carriage was chosen as the model for these tests since it represents one of the more common railway vehicle types and is the intended subject for application of the proposed TCMM.
Table 2 presents the nominal values of the main modal parameters of the vehicle model. These parameter values are in accordance with standard values used on similar numerical models [
10,
23].
The railway infrastructure was described by four equal sets of body elements, each including two rail elements, one sleeper, and two cylinders placed on top of each other. These cylinders, which were set according to the proposed MS, represented a ballast layer (grey cylinder) and a subgrade layer (brown cylinder). This part of the model also contained spring-damper elements to simulate the equivalent vertical stiffness behavior from each component. During a standard simulation, these body element sets follow their corresponding vehicle axle as the vehicle moves along the cartographic profile of the track, thus enabling the application of varying track stiffness profiles on the track at each infrastructure level (i.e., sleeper, ballast and subgrade levels). The sleeper spacing effect is simulated by a relatively small harmonic oscillation of the vertical stiffness at that level. However, these operational details were not relevant for the sensitivity analysis since those results were obtained in a standstill model simulation. But they are essential for the validation efforts of the proposed TCMM using numerical data, as shown in Morais et al. (2023).
Table 3 presents the nominal values used for the modal parameters of the infrastructure elements sets. The implemented modal stiffness and modal damping values are within the respective ranges used in the consulted bibliography, as are the density values used to calculate the mass for each layered element of the infrastructure [
23].
The main objective of this model was to obtain the resonance frequency value of the first vertical vibration mode of the simulated subgrade layer (fsim), which should be related with the natural frequency value of a real subgrade layer with similar characteristics as the ones being simulated. This value was obtained by performing modal analyses on the presented numerical model, via a tool in Simpack® called “Eigenvalues Analysis”. This tool calculates each resonance frequency present on the model, based on its initial conditions, and provides a representation of the kinematics of each vibration mode.
As presented in Morais et al. (2023), the subgrade vibration mode of interest from a MBS model with this overall configuration corresponds to the full-body oscillation of all infrastructure elements, coupled with the unsprung mass of the vehicle, over the adjoining spring elements. This vibration mode has the highest sensitivity to the subgrade’s modal properties. Thus, it is the best candidate to represent the equivalent vibration mode of the first resonance frequency of the subgrade layer. Further details on this connection are available in Morais et al. (2023), including a more thorough explanation of why this choice was made and its adequacy [
17].
Based on the nominal values selected for the model’s modal parameters (
Table 2 and
Table 3), an initial successful assessment of the validity of the proposed MS was performed by comparing the subgrade’s theoretical natural frequency value given by the combined model with the equivalent resonance frequency obtained by Simpack
®. The obtained frequency was sufficiently close to the theoretical model’s estimate, with the error being under 5%. Then, to further verify this topic, a sensitivity analysis was performed to check the direct or indirect influence that several parameters could have on the quality of this comparison; that is, to check if the obtainable error margin was still adequate in a broader range of the relevant model parameters, thus assessing the MS’s robustness. The parameters that could directly influence this comparison are those that enter into the combined model (Equations (1)–(7)). The indirect influencing parameters correspond to those that are not specifically mentioned in the combined model but could potentially influence the results.
Table 4 presents the values of the parameters that were used to obtain the sensitivity analysis results that will be shown in
Section 3. The only indirectly influencing parameter shown in this table is the car body mass (
mv.c). These specific ranges were selected with the intent of testing the limits of the MS from a mathematical standpoint, with less focus on trying to replicate a realistic situation. But each parameter range was still centered on the respective nominal values presented in
Table 2 and
Table 3. Each of these parameters was changed separately in each step of the sensitivity analysis, with the remaining parameters kept in the nominal values presented in
Table 2 and
Table 3.
For each of the parameter values presented in
Table 4, theoretical and numerically calculated frequency values were determined for multiple subgrade layer heights, namely, for [
2,
3,
4,
5,
6,
10] meters. This range, which includes the nominal subgrade layer height mentioned in
Table 3, was chosen to further test the capabilities of the proposed MS by simulating an even wider range of scenarios.