Numerical Testing of Switch Point Dynamics—A Curved Beam with a Variable Cross-Section

The article presents mathematical considerations on the dynamics of the springing switch point being an element of the railway junction. Due to the structure of the switch point, mathematical analysis was divided into two stages: The first stage refers to the analysis of the dynamics of the switch point as a beam of variable rectilinear stiffness to which three forces (coming from three closures of switch drives) placed in the initial section of the switch point are applied. The next stage of the analysis concerns an identical beam, but curved, with a variable cross-section. In both cases, the beam is subjected to a vertical force resulting from forces from the rail vehicle. The calculations refer to a switch point of 23 m length and a curvature radius R = 1200 m. The first stage of the switch point analysis refers to the movement of a rail vehicle on a straight track, and the second stage concerns the rail vehicle movement on a reverse path. This article also provides an analysis of mode vibrations of a curved beam with a variable cross-section, and variable inertia and stiffness moments (further in the article the changes will be referred to as beam parameter changes). It is assumed that the beam is loaded with vertical forces (coming) from a rail vehicle. The solution was found by applying the Ritz method, which served to present the fourth-order partial equations as ordinary differential ones. The numerical research whose results are given aimed to define how the changes in beam parameters and vertical load affect mode vibrations of the beam.


Introduction
Free vibrations of continuous mechanical systems, which include a switch point-a rail turnout element-do not solely depend on a load model (conservative or non-conservative one) or a type of support structures, but also on the change in both the cross-section, inertia moment and stiffness along its length. An additional quantity to consider in mathematical analysis is a switch point curve.
The dynamics of a nonlinear beam on a linear elastic foundation was tested by Akour [1], who investigated a simply supported beam with harmonic distributed load, on a one-parameter elastic foundation with damping included. The impact of damping coefficient, beam nonlinearity, and vibration frequency on beam stability was tested. An issue of nonlinear vibrations for Euler-Bernoulli beams was presented in [2]. Two methods were used to determine the frequency of vibrations: the variation iteration method (VIM) and the parameterized perturbation method (PPM). Both methods allow obtaining results without an introduction of a small amplitude parameter to the main equation that describes transverse vibrations of the beam.
The perturbation method was also used in paper [3]. They focused on a beam with a concentrated load perpendicular to its longitudinal axis, when the beam is located on a rigid one-parameter elastic The perturbation method was also used in paper [3]. They focused on a beam with a concentrated load perpendicular to its longitudinal axis, when the beam is located on a rigid one-parameter elastic foundation. Apart from the analysis of transverse vibrations, the authors provided the ways of analyzing the areas where a loaded beam breaks away from the rigid foundation used.
Transverse vibration analysis was also conducted in works [4,5], where analytical solutions to vibrations of a beam supported on springs (rotational or translational) on both ends, placed on two-parameter elastic foundation, were presented. The paper provides an explicit conclusion that the frequency of a system mode vibrations increases with the rising stiffness of the foundation and support springs.
An issue of transverse vibrations of two-segment beams was presented, among others, in [6]. The systems analyzed were divided into three basic groups depending on the shape of the cross-sectional area. The values of three first frequencies of mode vibrations at various mounting conditions were indicated. An identical analysis of changes in mode vibration frequency was carried out in [7], where the systems of three and more segments were taken into consideration. In the case of a system built of any finite number of segments, with an added discrete element in the form of concentrated mass or translation spring, the properties of Green's function [8] were used to determine the changes in mode values. Publications [9] dealt with systems of beams of a linearly variable cross-section, but only one of the main dimensions of the cross section was subject to change.
This paper presents theoretical and numerical tests on the stability and free vibrations of a curved beam with a variable cross section at vertical load. A curved beam of a variable cross-section, variable inertia moment, and variable stiffness along its length is examined. Motion equations and boundary conditions of analyzed systems are formulated on the grounds of total mechanical energy. The paper also includes the results of simulation tests on the course of changes in mode vibrations over the whole switch point structure: dimensionless load parameter-dimensionless mode frequency parameter. None of the publications mentioned take into consideration a beam with four simultaneously occurring variable quantities (variable cross-section, variable inertia moment, variable stiffness, and beam curve).

Structure Turnouts Railways
The turnout is one of the most important elements of the railway network: it is the element where railway disasters occur most often. There are two elements in the turnout which determine proper functioning of the railway turnout. These are a switch point and a frog [10,11] shown in Figure 1. The basic components and geometrical elements of the regular turnout are: beginning of turnout 1 (in pre-switch contact point); end of turnout 2 (in contact point behind the frog); turnout geometric center, 0, which is the intersection of the stock rail axis with the closure rail axis; radius R of the closure rail curve; turnout angle α (angle between the axes of the stock rail and the closure rail), slant of the turnout-tangent α-expressed as a vulgar fraction with one in the numerator; switch points, 4; rheostats, 5; frog, 6; wing rails, 7; guard rails, 8; connecting rails, 9; AOB triangle. The basic components and geometrical elements of the regular turnout are: beginning of turnout 1 (in pre-switch contact point); end of turnout 2 (in contact point behind the frog); turnout geometric center, 0, which is the intersection of the stock rail axis with the closure rail axis; radius R of the closure rail curve; turnout angle α (angle between the axes of the stock rail and the closure rail), slant of the turnout-tangent α-expressed as a vulgar fraction with one in the numerator; switch points, 4; rheostats, 5; frog, 6; wing rails, 7; guard rails, 8; connecting rails, 9; AOB triangle.
Turnouts for high-speed rail, V > 160 km/h, have curvature radii greater than 1000 m. In these turnouts, switch points have the following features; they are elastic elements with a permanently fixed Materials 2020, 13, 701 3 of 18 radius of curvature-from the point of view of mechanics that means that it is a curved beam with a variable cross-section, variable moment of inertia, and variable stiffness. In turnouts with smaller curvature radii, these are usually simple elements connected to the rail with an articulated element. The examination of such a beam (for R ≥ 1200 m) requires mathematical modeling to take into account all the elements that occur in the real object. In these systems, there are a minimum of 3 adjusting closures, in which the holding force for each should be greater than 7.5 kN. A mathematical model of a straight beam with a variable cross-section can be found in [12]. The results of these considerations were used to describe the switch point with smaller radii [13][14][15].
The switch point allows rail vehicle movement to be diverted into a straight or closure rail. The mathematical modeling process of the switch point needs to include forces coming from subsequent wheel sets; for example, for a rail vehicle composed of two bogies during the passage of the vehicle through the switch point, there are four forces in the resulting sequence: between sets in the bogie and between sets in two further bogies. These elements are considered in the modeling of the switch point as a beam subjected to a variable load [16][17][18].

Parameters Characterizing the Switch Point
Considering the switch point as a trapezoidal curved beam (Figure 2), lying on a continuous elastic substrate, and subjected to a vertical load, which passes through its neutral axis, let F mark the shear force in the cross-section of the beam at a distance x from the origin of the coordinates, and M is the moment of external forces relative to this cross-section. By F we understand the sum of all vertical forces acting on the beam from the side of the cross-section on which the origin of the coordinates is located. In [19][20][21], one can find mathematical models describing the dynamics of a beam with a variable cross-section, but they lack the impulse of forces acting on the surface in the vertical direction (wheel pressure on the switch point) and horizontal (the force holding a switch drive) [22,23]. Turnouts for high-speed rail, V > 160 km/h, have curvature radii greater than 1000 m. In these turnouts, switch points have the following features; they are elastic elements with a permanently fixed radius of curvature-from the point of view of mechanics that means that it is a curved beam with a variable cross-section, variable moment of inertia, and variable stiffness. In turnouts with smaller curvature radii, these are usually simple elements connected to the rail with an articulated element. The examination of such a beam (for R ≥ 1200 m) requires mathematical modeling to take into account all the elements that occur in the real object. In these systems, there are a minimum of 3 adjusting closures, in which the holding force for each should be greater than 7.5 kN. A mathematical model of a straight beam with a variable cross-section can be found in [12]. The results of these considerations were used to describe the switch point with smaller radii [13][14][15].
The switch point allows rail vehicle movement to be diverted into a straight or closure rail. The mathematical modeling process of the switch point needs to include forces coming from subsequent wheel sets; for example, for a rail vehicle composed of two bogies during the passage of the vehicle through the switch point, there are four forces in the resulting sequence: between sets in the bogie and between sets in two further bogies. These elements are considered in the modeling of the switch point as a beam subjected to a variable load [16][17][18].

Parameters Characterizing the Switch Point
Considering the switch point as a trapezoidal curved beam ( Figure 2), lying on a continuous elastic substrate, and subjected to a vertical load, which passes through its neutral axis, let F mark the shear force in the cross-section of the beam at a distance x from the origin of the coordinates, and M is the moment of external forces relative to this cross-section. By F we understand the sum of all vertical forces acting on the beam from the side of the cross-section on which the origin of the coordinates is located. In [19][20][21], one can find mathematical models describing the dynamics of a beam with a variable cross-section, but they lack the impulse of forces acting on the surface in the vertical direction (wheel pressure on the switch point) and horizontal (the force holding a switch drive) [22,23]. The initial assumptions were made to simplify the process of applying differential equations to the description of the movement of a curved beam with a variable cross-section, which is an element of a beam associated with two degrees of freedom. Thus, one rotary movement and one translational movement take place at both ends of the beam. The action of F force causes deformations (displacements) in the direction x of the narrower part of the switch point. In Figure 2, the beginning of the switch point is marked by A, the end by B. In turn, the section depth of the beam with a variable cross-section is indicated by ga and gb. The length of the beam is marked by L [24,25]. The initial assumptions were made to simplify the process of applying differential equations to the description of the movement of a curved beam with a variable cross-section, which is an element of a beam associated with two degrees of freedom. Thus, one rotary movement and one translational movement take place at both ends of the beam. The action of F force causes deformations (displacements) in the direction x of the narrower part of the switch point. In Figure 2, the beginning of the switch point is marked by A, the end by B. In turn, the section depth of the beam with a variable cross-section is indicated by g a and g b . The length of the beam is marked by L [24,25].

Switch Point Geometry as a Beam with a Variable Cross-Section
As for the structure of the beam geometry, the following assumptions have been made.

•
The y and z axes are the main central axes of inertia of the cross-section.

•
The cross-section and longitudinal section has variable quantities and its height h changes linearly along the beam (Figure 3).

•
The initial load in the coordinate system has zero speed, which changes later along the beam.

•
The beam stiffness of the constrained torsion is equal to zero.

•
The bending stiffness with respect to the y axis is slightly greater than the stiffness with respect to the z axis.

Switch Point Geometry as a Beam with a Variable Cross-Section
As for the structure of the beam geometry, the following assumptions have been made.

•
The y and z axes are the main central axes of inertia of the cross-section.

•
The cross-section and longitudinal section has variable quantities and its height h changes linearly along the beam (Figure 3).

•
The initial load in the coordinate system has zero speed, which changes later along the beam.

•
The beam stiffness of the constrained torsion is equal to zero.

•
The bending stiffness with respect to the y axis is slightly greater than the stiffness with respect to the z axis.  According to the cross-section of the beam shown in Figure 4a, the dependence on the cross-section area [26,27] is determined as follows, The moment of inertia I in the x-axis is defined by In mathematical considerations, the geometrical parameters of the cross-section of the beam marked with the red envelope (Figure 4a), were used. The profiles change all the values describing the structure along the x-axis, which directly affects the change of the cross-sectional area of the beam. This, in turn, changes the moments of inertia occurring at individual sections defined in the xyz axes in Figure 4b.

Switch Point Geometry as a Beam with a Variable Cross-Section
As for the structure of the beam geometry, the following assumptions have been made.

•
The y and z axes are the main central axes of inertia of the cross-section.

•
The cross-section and longitudinal section has variable quantities and its height h changes linearly along the beam (Figure 3).

•
The initial load in the coordinate system has zero speed, which changes later along the beam.

•
The beam stiffness of the constrained torsion is equal to zero.

•
The bending stiffness with respect to the y axis is slightly greater than the stiffness with respect to the z axis.  According to the cross-section of the beam shown in Figure 4a, the dependence on the cross-section area [26,27] is determined as follows, The moment of inertia I in the x-axis is defined by According to the cross-section of the beam shown in Figure 4a, the dependence on the cross-section area [26,27] is determined as follows, where, The moment of inertia I in the x-axis is defined by The moment of inertia I in the z-axis is described by the following equation.
The moment of inertia I in the y-axis was defined by the dependence

Beam with a Variable Cross-Section, Switch Point, Model of Assumption, and Real Object
For the switch point modeling process, the function characterizing the change in the quantities of moments of inertia along the length was determined. The results are shown in Figure 5.
The moment of inertia I in the z-axis is described by the following equation.
The moment of inertia I in the y-axis was defined by the dependence

Beam with a Variable Cross-Section, Switch Point, Model of Assumption, and Real Object
For the switch point modeling process, the function characterizing the change in the quantities of moments of inertia along the length was determined. The results are shown in Figure 5. These results were obtained based on measurements in the real object. The switch point shown in Figure 6a for a turnout with a radius of 1200 m will be shown as part of a curved beam. The assumption was that the switch point would be a curved beam of length L with a variable cross-section whose main axis runs along the central point of its cross-section [28,29]. Where the radius of curvature R defining the axis of the beam is large in relation to its section width h 5 , it can be assumed that under working stresses of a cross-section defined in this way, the beams produce identical effects as in the case of a straight beam ( → ∞) . Therefore, in further These results were obtained based on measurements in the real object. The switch point shown in Figure 6a for a turnout with a radius of 1200 m will be shown as part of a curved beam. The assumption was that the switch point would be a curved beam of length L with a variable cross-section whose main axis runs along the central point of its cross-section [28,29]. Where the radius of curvature R defining the axis of the beam is large in relation to its section width h R h > 5 , it can be assumed that under working stresses of a cross-section defined in this way, the beams produce identical effects as in the case of a straight beam (R → ∞). Therefore, in further considerations, it has been adopted that the curved beam meets the assumption R h < 5 .
Materials 2020, 13, 701 6 of 17 It was assumed that the radius would be larger than 1200 m and a mathematical description would be presented.
The switch point will be considered as a beam of a defined curvature and length s. The cross-section is a variable parameter in relation to its length in the planes x1 and x3 cut perpendicularly by two planes relative to the main axis. In the discussed model, F(L) indicated the forces acting on the beam (switch point), whereas all the moments of inertia appearing in the considered beam are marked by m(L). Both the F(L) and the m(L) parameters are determined per unit length of the main axis of the beam [30][31][32]. The phenomena occurring on such a defined beam can be divided into three components: peripheral ( circular force and circular moment of inertia ), radial (radial force and moment of radial inertia ), and normal (normal forces and moments of inertia). Figure 6b shows the model of a switch point defined by a curved beam with applied components. Forces and moments (normal, peripheral, and radial) acting on individual segments of the beam on the left side were marked by N, Q2, and , i , respectively. For the right side, the individual components are defined as follows; + Δ , + Δ , and + Δ , and the moments are defined by values + Δ , + Δ , and + Δ . The values in which the symbol Δ is displayed are considered variable. Parameter Δ determines the force vector in the switch point. As you can see in Figure 6b, these values will always be positive [33][34][35].
To obtain the state of balance in the analyzed model, the following equations must be satisfied, and It was assumed that the radius would be larger than 1200 m and a mathematical description would be presented.
The switch point will be considered as a beam of a defined curvature and length s. The cross-section is a variable parameter in relation to its length in the planes x 1 and x 3 cut perpendicularly by two planes relative to the main axis. In the discussed model, F(L) indicated the forces acting on the beam (switch point), whereas all the moments of inertia appearing in the considered beam are marked by m(L). Both the F(L) and the m(L) parameters are determined per unit length of the main axis of the beam [30][31][32].
The phenomena occurring on such a defined beam can be divided into three components: peripheral (p θ circular force and circular moment of inertia M θ ), radial (radial force p ζ and moment of radial inertia M ζ ), and normal (normal forces and moments of inertia). Figure 6b shows the model of a switch point defined by a curved beam with applied components. Forces and moments (normal, peripheral, and radial) acting on individual segments of the beam on the left side were marked by N, Q 2 , Q ζ and M θ , M 2 i M ζ , respectively. For the right side, the individual components are defined as follows;N + ∆N, Q 2 + ∆Q 2 , and Q ζ + ∆Q ζ , and the moments are defined by values M θ + ∆M θ , M 2 + ∆M 2 , and M ζ + ∆M ζ . The values in which the symbol ∆ is displayed are considered variable. Parameter ∆s determines the force vector in the switch point. As you can see in Figure 6b, these values will always be positive [33][34][35].
To obtain the state of balance in the analyzed model, the following equations must be satisfied, Assuming that and where R(L) is the radius of the switch point from the axis of beam curvature. Dividing expressions (6) and (7) by ∆s at the limit ∆s → 0 and θ → 0 allowed obtaining dependencies so that the moments of inertia of the switch point M ζ , M z were determined. Therefore, the following equations are true [36], and Considering the mutual relations between Equations (10) and (11), we obtained [37,38] dM θ Expressions (11) and (12) determine the state of balance for a flat switch point. For a straight beam → ∞ , s → x 1 , in turn the indexes ζ → 3 and thus θ → 1 . From the equations presented above, the moment of inertia M 2 described by the equation was determined [39][40][41]: where What is more, from the same relationships, you can get expressions for moments M θ and M ζ described as follows [42][43][44], According to the analysis carried out in the case of a switch point, individual moments of inertia are given directly from curved places and in the case of a straight rail at its ends. It should be added that the switch point is exposed to the same phenomena as the straight rail. Therefore, we decided to carry out an analysis determining the bending effect of the beam. Differential equations that describe the dynamics of the switch point movement (beam with a variable cross-section) are solved to obtain the required stiffness of the surface and beam deformations, and can also be used to estimate the mode vibrations described by the mathematical equations of the object [45][46][47].
To determine motion equations for a curved beam with a variable cross-section, first, the model of the tested object should be correctly defined. For this purpose, a curved beam with a variable cross-section and with the vertical force from the wheelset operating on its upper surface is considered. This situation is illustrated in Figure 7, where the value of u means deflections of the beam in three directions: longitudinal, vertical, and transverse. The radius determining the beam curvature is defined by R, and its length by L. Among x and y axes, which coincide with the main axis of the system defining the cross-section of the tested element, the x-axis is tangent to the switch point arc. Dislocations (deflections) u x , u y , and u z of the center of gravity should be determined for each section of the beam cross-section. The sizes θ x , θ y , and θ z describe rotations around the three axes [48][49][50]. According to the analysis carried out in the case of a switch point, individual moments of inertia are given directly from curved places and in the case of a straight rail at its ends. It should be added that the switch point is exposed to the same phenomena as the straight rail. Therefore, we decided to carry out an analysis determining the bending effect of the beam. Differential equations that describe the dynamics of the switch point movement (beam with a variable cross-section) are solved to obtain the required stiffness of the surface and beam deformations, and can also be used to estimate the mode vibrations described by the mathematical equations of the object [45][46][47].
To determine motion equations for a curved beam with a variable cross-section, first, the model of the tested object should be correctly defined. For this purpose, a curved beam with a variable cross-section and with the vertical force from the wheelset operating on its upper surface is considered. This situation is illustrated in Figure 7, where the value of u means deflections of the beam in three directions: longitudinal, vertical, and transverse. The radius determining the beam curvature is defined by R, and its length by L. Among x and y axes, which coincide with the main axis of the system defining the cross-section of the tested element, the x-axis is tangent to the switch point arc. Dislocations (deflections) ux, uy, and uz of the center of gravity should be determined for each section of the beam cross-section. The sizes θx, θy, and θz describe rotations around the three axes [48][49][50]. In further considerations, the accepted principle is that all deformations occurring in the analyzed model are small, so in further mathematical analysis a linear theory was used. It was also assumed that the beam has very little deformation resistance. Thus, differential equations defining the movement of a curved beam with a variable cross-section are as follows.
The movement specifying an axial shift: The movement specifying radial dislocation: In further considerations, the accepted principle is that all deformations occurring in the analyzed model are small, so in further mathematical analysis a linear theory was used. It was also assumed that the beam has very little deformation resistance. Thus, differential equations defining the movement of a curved beam with a variable cross-section are as follows.
The movement specifying an axial shift: The movement specifying radial dislocation: The movement defining a vertical displacement of a curved beam with a variable cross-section: The motion defining the torsional movement: where the entire description of the movement dynamics of the curved beam with a variable cross-section is shown for the longitudinal axis x, the quantities E and G denote, respectively, the moduli of elasticity of the beam material and its stiffness. A: cross-sectional area, and I y and I z are, respectively, moments of inertia determined around the y-axis and z. J: constant torsional. From the deductions, note that the differential equations determining the movement for the plane of displacements u x and u z are independent of the displacements occurring outside the plane u y and θ x . In addition, the differential equations determined for the axial displacement u x and for the radial shift u z are reciprocally incorporated, and the situation looks the same for the equation determining the vertical displacements of u y and torsional movement. Next, it is assumed that both ends of the beam are supported, which translates into the lack of dynamics of the movement of the flexible displacement and the rotational movement. Thus, their first derivatives are not "zero". In general, this can be triggered by the action of a rail vehicle traveling on a beam. In mathematical considerations, the vehicle movement may be replaced by a vertical moving force as shown in Figure 6. The movement in the horizontal direction along a curved beam with a variable cross-section on the upper surface of the tested element is considered as a perpendicular force F moving at v velocity and acting on the upper surface of the beam. Load of the object moving on the beam is expressed by F 0 = −mg, where g is gravitational acceleration. The equations of vertical vibrations of a curved beam with a variable cross-section are determined based on the dependence (16)(17)(18)(19). However, in this case, to better reflect the actual conditions in modeling the phenomena of beam motion dynamics, nonlinearity is taken into account.
where m stands for the mass per unit length, ρ is the density of the curved beam, and δ(t) is the Dirac impulse. It can be seen that the right side of the first Equation (16) can be expressed as the sum of a series of sinusoidal functions that satisfy the boundary conditions: where q yi means the generalized coordinate of the tested beam section for the i-th note. The expression of the steering angle is marked by θ x and must be defined by the Equation (16). Subsequently, by replacing Equation (17) with the last Equation (18) and taking into account the boundary conditions regarding the angle of twist of the beam, we obtained where the parameter ; then, as a result of the mutual relationship, we where q 0i means i-th coordinate for the steering angle θ x . Further considerations refer to defining the deflection of the beam and determining the turning angle: where by q y1 and q 01 , respectively, the first generalized coordinates for displacement u y and steering angle θ x have been determined. The nominal model used to determine the vibrations of a curved beam with a variable cross-section is shown in Figure 8. where means i-th coordinate for the steering angle . Further considerations refer to defining the deflection of the beam and determining the turning angle: ( , ) = sin , ( , ) = sin (24) where by and , respectively, the first generalized coordinates for displacement and steering angle have been determined. The nominal model used to determine the vibrations of a curved beam with a variable cross-section is shown in Figure 8. To solve the differential equations presented in (16), the Galerkin method should be used. Thus, we multiply both sides of the first equation by variable and the second equation by variable.
Then, we apply approximation to the first form for quantity and in Equation (3.18). In the next step, we perform the operation of combining the two differential equations, each with respect to x in the range from 0 to L. The result of this operation is the mathematical relationships: where = − 1 + = − 1 1 + Because variables and are arbitrary, Equation (25) is reduced to the following, To solve the differential equations presented in (16), the Galerkin method should be used. Thus, we multiply both sides of the first equation by δu y variable and the second equation by δθ x variable. Then, we apply approximation to the first form for quantity u y and θ x in Equation (3.18). In the next step, we perform the operation of combining the two differential equations, each with respect to x in the range from 0 to L. The result of this operation is the mathematical relationships: where Because variables δq y1 and δq 01 are arbitrary, Equation (25) is reduced to the following, The general solution to the Equation (27) which is a homogeneous one is q y1 = q y1h + q y1p q 01 = q 01h + q 01p (28) where indices h and p denote the solution of Equation (27), respectively. Thus, the solutions can be presented as follows, where ω v1 is the basic vibration frequency determined in the vertical direction of a curved beam with a variable cross-section, and the parameters h 1 , h 2 , k 1 , k 2 are, respectively, coefficients that should be determined on the basis of the initial boundary conditions. In the next step, on the basis of the expression (25), which is substituted into Equation (27), the F 0 value was determined: Then, by performing an operation of substituting the determinant equal to zero into Equation (24), the fundamental frequency ω v1 can be determined from the relationship To obtain individual higher harmonic frequencies, it can be assumed that where p y1 and p 01 are the value of the amplitude of vibrations occurring in a given object. Next, by carrying out the substitution (32) to Equation (27) the following was obtained, From Equation (33), the solution for p y1 was determined: where S v1 specifies the velocity parameter for vertical vibrations occurring on a curved beam, which is defined as the ratio of the frequency of generated vibrations πv L transmitted by the load to the fundamental frequency ω v1 .
In a further stage of mathematical considerations it was assumed that the initial displacement and the speed of the load are zero if the object does not move. Such conditions should be taken into account when solving the differential equations defined in (28) or performing summation of solutions of equations (27) and (29). The coefficients in the solution defined by h 1 = −p 1 S v1 and h 2 = 0 result from the solution of the coefficient q y1 from According to the assumption that F 0 = −m v g and Equation (34), the above equation can be rewritten as where the amplitude function ψ v1 (t) is determined by As a result, vertical displacements of the curved beam with a variable cross-section are expressed by the equation For the middle section x = L 2 of a curved beam with a variable cross-section, the vertical displacements are: The presented mathematical analysis, which aimed to determine the equations of the movements for a curved beam with a variable cross-section refers only to obtaining a solution that takes into account the adopted boundary conditions. In practice, the load caused by a rail vehicle passing through the switch point is so short that it induces a temporary process. The results regarding higher vibration frequencies occurring on the tested element depend on the location of the load on the switch point and the speed of the rail vehicle. In the further part of the article, results will be presented in a form of higher harmonic vibrations for a 1200-m turnout.

Results of the Simulation Study
The numerical solutions were obtained with the finite element calculation program-ANSYS, which allows modeling the behavior of structural systems, while taking into account the diverse but complementary material properties of beam objects in the area of construction and mutual cooperation of these construction materials. Computer simulations of this type contribute to the reduction of costs associated with conducting experimental design tests. They do not replace these tests but allow more effective planning of experiments. They enable the analysis of complex issues of beam behavior during loading in the entire strain range from the linear-elastic range, through the elastic deformation phases in the areas of tensile stress causing initial and advanced scratch states, and in areas of compressive stress causing switch point deflections that lead to the destruction of the railway turnout construction elements.
Based on the presented nominal models, mathematical models were determined and calculations were made for a beam having a turnout switch point with a radius of 1200 m, and the forces that arise when the rail vehicle travels through the turnout on stock and closure rail were taken into account. The results of simulation tests were performed in the ANSYS software.
For computer simulations, a curved beam with a variable cross-section in the shape of a rail with the profile: 49 E1 (S49), 60E1 (UIC60) of steel grade 900A was used. Technical data of the UIC60 profile can be found in publication [51]. In the numerical model under study, rounding and chamfering of the switch point were omitted to reduce the calculation time. On the curved beam model with a variable cross-section, a tetra-type finite element grid of 24082 elements was applied. The vibration analysis was obtained by combining strength and frequency analysis. The results of calculations from the strength analysis were implemented as input data for frequency analysis [52]. To induce mechanical loads, the model was fixed on the left side and the supports were added, one fixed and one movable. The parameters used in the simulation are shown in Table 1. The calculations were carried out in one bending case: • Bending in the direction of (x axis). The beam is curved in the plane of action of the concentrated force with the value F 0 = 12 kN. The force is applied at the middle connector of the top edge of the free end of the beam and works in the direction of the x axis (negative). Theoretical value of the maximum displacement in this direction according to the work is-0.2623 cm.
Simulation research was also carried out based on work [52]. Own forms were determined: first, second, third, and fourth for x equal to the length of the switch point, for x = 0 the switch point is free (Figure 9). The simulation results were determined for different nominal forces affecting the switch points. Figures 10 and 11 show the designated own forms for a curved beam with variable stiffness (change of section area, stiffness and moments of inertia) characterizing the properties of the switch point.
The purpose of the calculations was to determine the vibrations caused by the action of a 12 kN load. Figure 10 shows the waveforms of individual harmonics (from the first to the sixth one) obtained from computer simulations. However, the distribution of beam deflections is shown in Figure 11. The obtained results of vibrations and deflections indicate the influence of the load on the formation of strains and displacements in the analyzed case. The numerical simulation performed with the ANSYS program allowed us to determine displacements at the load acting on the upper surface of the curved beam with a variable cross-section. The maximum deflection of the beam, i.e., elastic displacement in the x-axis direction, occurs at the end of the tested element and amounts to 1.1398 mm, whereas the maximum longitudinal displacements of the beam with a value of 3.2194 mm appear at the end of it.
Simulation research was also carried out based on work [52]. Own forms were determined: first, second, third, and fourth for x equal to the length of the switch point, for x = 0 the switch point is free (Figure 9). The simulation results were determined for different nominal forces affecting the switch points. Figures 10-11 show the designated own forms for a curved beam with variable stiffness (change of section area, stiffness and moments of inertia) characterizing the properties of the switch point.   Simulation research was also carried out based on work [52]. Own forms were determined: first, second, third, and fourth for x equal to the length of the switch point, for x = 0 the switch point is free (Figure 9). The simulation results were determined for different nominal forces affecting the switch points. Figures 10-11 show the designated own forms for a curved beam with variable stiffness (change of section area, stiffness and moments of inertia) characterizing the properties of the switch point.    The purpose of the calculations was to determine the vibrations caused by the action of a 12 kN load. Figure 10 shows the waveforms of individual harmonics (from the first to the sixth one) obtained from computer simulations. However, the distribution of beam deflections is shown in Figure 11. The obtained results of vibrations and deflections indicate the influence of the load on the formation of strains and displacements in the analyzed case. The numerical simulation performed with the ANSYS program allowed us to determine displacements at the load acting on the upper surface of the curved beam with a variable cross-section. The maximum deflection of the beam, i.e., elastic displacement in the x-axis direction, occurs at the end of the tested element and amounts to 1.1398 mm, whereas the maximum longitudinal displacements of the beam with a value of 3.2194 mm appear at the end of it.

Conclusions
The numerical analysis carried out in this work made it possible to determine the displacements and vibrations of higher harmonics in the model of a curved beam with variable cross-section subjected to load. The presented method can be used for testing the dynamics of a rail vehicle moving on the switch point with variable cross-section, variable stiffness, and variable inertia moment for greater radii of curvature. The results obtained in the conducted simulations enable the comparison of the switch point mode frequency with own vibration frequency of a rail vehicle mass, which may lead to the creation of dynamic couplings between the turnout and the rail vehicle (Medelsztam theory).
The presented results of the numerical analysis were obtained in the ANSYS 16 Workbench module.
Solving continuous beams with analytical methods is to a large extent complicated and time-consuming, whereas ANSYS enables quick numerical calculations using FEM (Finite Element Method).
The conducted numerical tests allowed concluding that greater system stiffness and the stiffness of applied Winkler foundation result in a higher frequency of system mode vibration. According to Figure 9 diagrams, the frequency of transverse vibrations depends mainly on the applied support. Also, a generated vertical force and a very location of a segment of a curved beam with variable parameters in relation to supports significantly affect vibration frequency. Figure 11. The fourth own form for a curved beam with a variable cross-section (turnout switch point).

Conclusions
The numerical analysis carried out in this work made it possible to determine the displacements and vibrations of higher harmonics in the model of a curved beam with variable cross-section subjected to load. The presented method can be used for testing the dynamics of a rail vehicle moving on the switch point with variable cross-section, variable stiffness, and variable inertia moment for greater radii of curvature. The results obtained in the conducted simulations enable the comparison of the switch point mode frequency with own vibration frequency of a rail vehicle mass, which may lead to the creation of dynamic couplings between the turnout and the rail vehicle (Medelsztam theory).
The presented results of the numerical analysis were obtained in the ANSYS 16 Workbench module. Solving continuous beams with analytical methods is to a large extent complicated and time-consuming, whereas ANSYS enables quick numerical calculations using FEM (Finite Element Method).
The conducted numerical tests allowed concluding that greater system stiffness and the stiffness of applied Winkler foundation result in a higher frequency of system mode vibration. According to Figure 9 diagrams, the frequency of transverse vibrations depends mainly on the applied support. Also, a generated vertical force and a very location of a segment of a curved beam with variable parameters in relation to supports significantly affect vibration frequency.
In certain cases, the equations derived may be used for the verification of computer calculations and, in the authors' opinion, have some educational value. They may also be applied in calculations of other structures whose shape corresponds to a linear or quadratic relation between the inertia moment and the length of a beam.
The approximation of variable parameter beam shape affects the mode vibration frequency value depending on parameters that describe a mechanical system-a switch point and values of a load acting on it. A change of cross-section and beam curve may result either in an increase or a decrease of mode vibration frequency values. The maximum load value being equal to zero value of basic mode vibration frequency depends mainly on the value of parameters describing a beam's shape, geometrical parameters and its load structure. Maximum forces obtained from the kinetic criterion of stability loss are reflected in the results of numerical analyses of the system, which were conducted on the grounds of the static criterion of stability loss.
In further research work, the authors will focus on testing the impact of vertical force on the possibility of modification of transverse vibration frequencies when other support techniques are used and a perpendicular force from a switch drive is considered.
Author Contributions: All authors contributed equally to this article. All authors (J.K. and R.K.) have read and agreed to the published version of the manuscript. All authors have read and agreed to the published version of the manuscript.