Mathematical Problems in Engineering

Volume 2014 (2014), Article ID 679817, 13 pages

http://dx.doi.org/10.1155/2014/679817

Research Article

Analytical Method of Modelling the Geometric System of Communication Route

Wladyslaw Koc

Faculty of Civil and Environmental Engineering, Gdansk University of Technology, Narutowicza 11/12, 80-233 Gdansk, Poland

Received 10 October 2013; Accepted 14 April 2014; Published 11 May 2014

Academic Editor: Youqing Wang

Copyright © 2014 Wladyslaw Koc. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The paper presents a new analytical approach to modelling the curvature of a communication route by making use of differential equations. The method makes it possible to identify both linear and nonlinear curvature. It enables us to join curves of the same or opposite signs of curvature. Solutions of problems for linear change of curvature and selected variants of nonlinear curvature in polynomial and trigonometric form were analyzed. A comparison of determined horizontal transition curves was made and examples of negotiating these curves into a geometric system were given.

1. Introduction

A rapid progress in computer calculation technique algorithms based on current knowledge and inspired by further prospecting is also taken into account in the field of geometrical shaping of communication routes. In the sphere of shaping the communication routes, vehicular roads and railway lines, the development of satellite measuring technique GNSS [1, 2] also plays a crucial role. A particular significance here will be attached to the question of vehicular dynamics. Therefore the modelling of curvature will obtain a key role.

This subject also includes the design of elements of the route of diversified curvature [3]. The elaboration of a method for a precise determination of the coordinate route points seems to be the most appropriate area of action. To perform this activity use is made of tools necessary to provide analytical solutions, that is, the most advantageous for practical application. The junction of elements of a vehicular route or railway line of diversified curvature should ensure a continuous change of the unbalanced side acceleration, in an advantageous way for interaction dynamics within the road-vehicle system. This is strictly connected with the proper shaping of the curvature.

The most developed investigation branch concerning this problem has been for years the study of transition curves connecting straight lengths of route with circular arcs. The problem of transition curves is related to vehicular roads and railway routes. However, it is easy to note that there is a distinct disproportion in taking interest in this problem. A search for new solutions is going on this sphere on vehicular roads (e.g., [4–12]). As regards the railway lines, the situation is entirely different. Publications on new transition curves are relatively rare and in the majority were published long ago (e.g., [13–17]).

The application of transition curves is aimed at ensuring a continuous change of the unbalanced side acceleration between intervals of motor road (or railway route) of diversified curvature in an advantageous way for the interaction dynamics within the road-vehicle system. Such a requirement relates to all sorts of transition curves. Under this situation it might appear that there is one particular algorithm to create them, common to all the analyzed curves. In fact most of the solutions that have been known so far appear independently and bear various names (sometimes coming from the name of their author). A general knowledge of the determination of transition curves equations would make it possible to compare various forms of curves with each other and to make an assessment of their practical application.

Twenty years ago this issue had already been sufficiently explained with respect mainly to railway routes [16]. Attention was then concentrated on working out a technique of identification of unbalanced accelerations occurring on various types of transition curves. It was based on a comparative analysis of some selected transition curves with the use of a dynamic model. In the method acceleration was a factor that excited the transverse vibrations of the carriage [18]. The main conclusion resulting from the considerations was to prove the relation existing between the response of the system and the class of the excitation function. The dynamic interactions were smaller (i.e., more advantageous) if the class of the function was higher. It turned out that evidently the largest acceleration values were noted on the cubical parabola (class of function C⁰). With respect to Bloss curve and a cosine curve (class C¹) they are significantly smaller, whereas on sine curve (class C²) they are the smallest.

After years the mentioned identification method of unbalanced accelerations [16] provided a source of inspiration for the elaboration of a new technique for modelling the geometric system of the communication route based on joining two points of the route with diversified curvature.

2. Analytical Method of Modelling the Curvature

A measure of the route bending is the ratio of the angle which determines the direction of the vehicle’s longitudinal axis after covering a certain arc (Figure 1). The curvature of curve at point is called the boundary which is aimed at by the relation of acute angle between tangents to curve at points and , to the length of arc when point tends along curve to point

Figure 1: Schematic diagram for the explanation of the notion of curvature.

If the operation procedure by the use of rectangular coordinates , is to be continued, it is necessary to take into consideration the sign of the curvature. As illustrated in Figure 1 angle ; the curvature is assumed here that it has a positive value which, as can be seen, accompanies the curves with downward convexities. A negative value of curvature corresponds to angle and appears on curves of an upward convexity.

Making a generalization of the identification method applied to unbalanced accelerations occurring on various types of transition curves [16], it is possible to search for curvature function among the solutions of the differential equation

with conditions for the transition curve at the outset (for ) and at the end (for )

where and indicate the curvature magnitudes on both ends of the curve (taking note of an adequate sign). The order of the differential equation (2) is , and the obtained function is of class in the range , where .

The application of the method under consideration makes it possible to combine various geometric elements, for example, a straight with a circular arc, and also, after introducing an adequate sign of curvature, two circular arcs of a consistent run or converse arcs.

On determining the curvature function a fundamental task is to find the coordinates that correspond to a curve in the rectangular coordinate system , . A solution for such a system is required by the satellite measuring technique GNSS [19–21], which becomes an aid in the determination of the coordinate points in a uniform, 3D-system of coordinates WGS 84 (the world geodetic system 1984). Next the measured ellipsoidal coordinates (GPS) are transformed to Gauss-Kruger conformal coordinates [22]. An assumption is made that the beginning of system , is at the terminal point of the input curve (with curvature ), while the axis of abscissae is tangent to this curve at this point.

The equation of the sought after connection can be written in parametric form

Parameter is the position of a given point along the curve length. Function is determined on the basis of the formula

The presented method of modelling the curvature has a universal character. It can be applied to both vehicular roads and railway routes. In the case of railway routes there arises an additional possibility, for modelling in a similar way, the superelevation ramp, understood here as a determined difference of heights of the rail. The principle of operation is similar to curvature modelling. The only differences are that has to be replaced by the value of superelevation at the beginning of the superelevation ramp and has to be replaced by the value of superelevation at the end of the superelevation ramp. The ramp length, of course, corresponds to the length of the transition curve. Using the same procedure it is possible to find the unbalanced acceleration .

3. Transition Curve of Linear Curvature

3.1. Finding Out the Curvature Equation

Linear change of curvature along a defined length is obtained by assuming two principal conditions

and a differential equation

After the determination of constants the solution of the differential problem (6), (7) is as follows:

Function is obtained from formula (5). In the case under consideration

3.2. The Determination of the Transition Curve Length (for Railway Lines)

With regard to rail tracks the transition curve must satisfy two kinematic conditions

where is the maximum value of acceleration increment on transition curve in m/s³, is permissible value of acceleration increment in m/s³, is the maximum speed of lifting the wheel on superelevation ramp in mm/s, and is permissible value of speed of lifting the wheel in mm/s.

The unbalanced acceleration occurring on the transition curve results from speed of the train, magnitude of curvature, and ordinates of the superelevation ramp; its course is similar to and . Analogously with (8)

where and are accelerations at ends of the transition curve

( and are ordinates of the superelevation ramp at its ends).

On the assumption that the speed of the train = const., the formula for acceleration increment is as follows:

As can be seen the acceleration increment is here a constant value. However, from condition (10) it follows that

The superelevation ramp equation in the case under consideration is similar to (8) and (12)

The speed of lifting the wheel on the superelevation ramp (assuming that the travelling speed of the train = const.) is determined by the equation

Thus, also is here a constant value, but from condition (11) it follows that

The assumed length of the transition curve must fulfill the condition

3.3. Coordinates of Transition Curve Connecting Uniform Curvatures

The determination of and by using (4) will need the expansion of integrands into Maclaurin series. After completing the whole procedure the following parametric equations are obtained:

Equations (20) for describe the curve in the form of clothoid used to connect a straight with a circular arc.

3.4. Coordinates of Transition Curve Connecting Inverse Curvatures

At inverse arcs (i.e., with diversified curvature signs) there arises the problem related to the determined coordinates and , due to the course of function . For function is a monotonic function, whereas for in the diagram of function there appears an extremum at point , where

Values and follow from the relations:

Figure 2 illustrates a scheme of angle for the connection of inverse circular arcs of radii  m and  m. The length of the transition curve  m was determined for speed value  km/h using the procedure as described at Section 3.2 (the assumed superelevation values were  mm and  mm).

Figure 2: Diagram of angle defined by (9) for circular arcs connection of curvatures  rad/m and  rad/m by means of transition curve of length  m.

In this situation the parametric equations of transition curve (20) are effective for ; for functions and should be expanded to Taylor series. After integration of the equations take the following form:

On completing the parametric equations (20) and (23) along with (24) a correct solution of the task for the case of diversified curvature signs is obtained (Figure 3). From a practical point of view the determination of the magnitude of the tangent slope angle at the end of the transition curve is very important. It amounts to

Figure 3: Diagram of horizontal ordinates defined by (20), (23), and (24) for the joint of circular arcs of curvatures  rad/m and  rad/m with transition curve of length  m (using different horizontal and vertical scales).

A knowledge of makes it possible to add the transition curve to the other one having curvature at its initial point (satisfying the tangency condition of both the curves).

4. Transition Curve of Polynomial Curvature

From the point of view of vehicles’ dynamics, the transition curve of linear change of curvature is not the most advantageous solution [23]. A definitely better one is nonlinear solution whose characteristics are interesting enough to be compared with the Bezier curves preferable lately in literature [4, 7].

The differential equation (2) allows to appoint an unlimited number of curves with nonlinear curvature change. In the case of geometric layout of communication routes it will be appropriate to consider the transition curves of polynomial and trigonometric curvature. In fact such forms of curves are used in the standard ways to connect a straight line with a circular arc.

4.1. Formulation of the Curvature Equation

To formulate the curvature equation in the polynomial form, let us assume the following boundary conditions:

and the differential equation

After solving the differential problem (26), (27) the following formula for curvature is obtained:

The form of function determined on the basis of (5) and (28) is as follows:

On rail routes the length of the transition curve is determined like at Section 3.2 for solution of linear curvature. The following conditions should be fulfilled here:

The adopted length of the transition curve must satisfy condition (19).

4.2. Coordinates of Transition Curve Connecting Uniform Curvatures

On expanding functions and into Maclaurin series and after integration, the following parametric equations are acquired:

Equations (31) for describe the so-called Bloss curve which can be used to join a straight to a circular arc.

4.3. Coordinates of Transition Curve Connecting Inverse Curvatures

When the arcs are inverse (i.e., the curvature signs are different, Figure 4) there occurs a problem relating to the determined coordinates and which results from function . For function is a monotonic function, while for in the diagram of function there appears extremum at point , where

Figure 4: Diagram of curvature described by (28) for the connection of circular arcs of curvatures  rad/m and  rad/m with the use of transition curve of length  m.

Parametric equations of transition curve (31) are valid for . The determination of value makes it necessary to solve (32) which is cubic. The sought-after value , satisfying conditions of the problem is defined by the formula

where angle follows from the relation

Under such circumstances in order to determine parametric equations of the transition curve for , functions and should be expanded into Taylor series; after integration the following parametric equations of the transition curve are obtained:

where

Making use of parametric equations (31) and also (35) one can obtain a correct solution of the problem with diversified signs of curvature.

Value of the tangent slope at the end of the transition curve makes it possible to connect to it a curve of curvature satisfying simultaneously the tangency condition of both the curves. The above value is the same as in the case of linear curvature. It is defined by (25).

5. Transition Curve of Curvature in Trigonometric Form

5.1. Formulation of Curvature Equation

As before boundary conditions (26) and another differential equation are assumed.

As soon as the differential problem (26), (37) has been solved, a formula for curvature, is obtained

The form of function determined by the use of (5) and (38) is as follows:

The length of the transition curve on railway routes follows from conditions:

The assumed length of the curve must satisfy the condition (19).

5.2. Ordinates of Transition Curve Connecting Uniform Curvatures

After the expansion of functions and into Maclaurin series to be next integrated the following parametric equations are obtained:

Equations (41) for describe the curve known as cosinusoid which can be used to connect a straight with a circular arc.

5.3. Transition Curve Coordinates Connecting Inverse Curvatures

Near inverse arcs in the diagram of function there appears the extremum at point , where

The transition curve parametric equations (41) are effective for . To determine value one must solve (42). The equation can be written in the form

It is easy to prove that for it is necessary to satisfy the condition , thus . The formula for length is as follows:

For functions and should be expanded into Taylor series. The integration provides the following equations:

where

In expressions with double signs and , the upper one corresponds to the case , , while the lower one refers to , .

The value of tangent slope at the end of the transition curve makes it possible to connect to it a curve of curvature satisfying the tangency condition of both the curves. The value is the same as in the previous cases. It is determined by (25).

6. Comparison of Horizontal Ordinates of Determined Transition Curves

6.1. The Shaping of Horizontal Ordinates of Comparable Curves

If the analysis is carried out for railway routes, then the necessity of meeting kinematic conditions (10) and (11) will require a diversification of the lengths of the compared with each other curves (using the same speed of trains). Moreover, it is necessary to individually consider cases of connecting homogeneous curvatures () and inverse arcs ().

It is assumed that at joints of consistent arcs with radii  m and  m an attempt is made to attain a speed  km/h. Hence it is necessary to employ on arcs superelevation of  mm and  mm, respectively. The retention of the same kinematic parameters can provide diversified lengths of the transition curve:

(i)for a curve of linear curvature  m,(ii)for a curve of polynomial curvature  m,(iii)for a curve of trigonometric curvature  m.

The above length relations result directly from the equations given at Sections 3.2, 4.1, and 5.1.

Figure 5 gives diagrams of horizontal ordinates for compared transition curves of positive and negative curvature values. As seen the curvature ordinates of linear curvature evidently deviate from the ordinates of other curves. In turn, curves of polynomial and trigonometric curvature are in principle consistent. A difference can only be noted in the final area. And so, the class of the function describing the curvature plays here a significant part: with respect to a curve of linear curvature it is a function of class C⁰, while regarding the other two curves it is a function of class C¹.

Figure 5: Diagrams of horizontal ordinates of curves joining consistent arcs for positive and negative values of curvature (using different horizontal and vertical scales): curve of linear curvature: line in violet colour, curve of polynomial curvature in green colour, and curve of trigonometric curvature in red colour.

At the joint of inverse arcs of radii  m and  m an attempt is made to reach the speed of  km/h. Hence it follows that it is necessary to apply superelevation along arcs attaining  mm and  mm, respectively. By the use the same kinematic parameters it is possible to obtain the following lengths of the transition curve:

(i)for a curve of linear curvature  m,(ii)for a curve of polynomial curvature  m,(iii)for a curve of trigonometric curvature  m.

Figure 6 presents diagrams of horizontal ordinates of compared transition curves for positive and negative values of curvature.

Figure 6: Diagrams of horizontal ordinates of curves connecting inverse arcs for positive and negative values of curvature (using different horizontal and vertical scales): the curve of linear curvature is in violet color, the curve of polynomial curvature is in green color, and the curve of trigonometric curvature is in red color.

As can be seen, like in the case of consistent arcs, the curve ordinates of linear curvature deviate significantly from the ordinates of the remaining curves, but the shape of these curves also differs one from the other. The class of function describing the curvature is still important. However, in the range of the same class various geometric solutions are possible.

6.2. Negotiation of Transition Curves in a Geometric System

The transition curves do not appear alone in a geometric system but must be connected with some arcs in the neighborhood. If the input arc is circular, then its equation is as follows:

(i)for

(ii)for

 where .

However, if the output arc is also a circular arc, then it is expressed in the following equation:

(i)for

(ii)for

 where .

In (49) and (50) denotes the abscissae of the transition curve end, is the ordinate of the transition curve end, while is the value of tangent at the end of the transition curve. Values for numerical examples under consideration are given in Table 1.

Table 1: Comparison of values , , and for numerical examples analyzed.

Curve
[m]
[rad]

[m]
[m]

Case  1:  rad/m and  rad/m

Curve C⁰
40
0,045238
0,045269
39,988
0,825

Curve C¹ (p)
60
0,067857
0,067961
59,963
1,821

Curve C¹ (t)
62,832
0,071060
0,071180
62,796
1,940

Case  2:  rad/m and  rad/m

Curve C⁰
60
−0,017857
−0,017859
59,999
0,143

Curve C¹ (p)
90
−0,026786
−0,026792
89,991
0,627

Curve C¹ (t)
94,248
−0,028050
−0,028057
94,238
0,714

Case  3:  rad/m and  rad/m

Curve C⁰
60
0,01786
0,017859
59,999
−0,143

Curve C1 (p)
90
0,026786
0,026792
89,991
−0,627

Curve C¹ (t)
94,248
0,028050
0,028057
94,238
−0,714

Case  4:  rad/m and  rad/m

Curve C⁰
40
−0,045238
−0,045269
39,988
−0,825

Curve C¹ (p)
60
−0,067857
−0,067961
59,969
−1,821

Curve C¹ (t)
62,832
−0,071060
−0,071180
62,976
−1,940

Figures 7 and 8 illustrate the geometric systems made up of two circular arcs of consistent curvature, connected with various types of transition curves. As seen the application of a given type of transition curve leads to an adequate position of the output circular arc. The position of this arc by the use of a linear transition curve clearly deviates from such cases where advantage is taken of curves with polynomial and trigonometrical curvature. With regard to the latter a mutual position of the circular arc varies insignificantly.

Figure 7: Examples of geometric systems consisting of two circular arcs with consistent curvature (brown line) connected with different transition curves, for positive values of curvature (using different horizontal and vertical scales): curve of linear curvature in violet line, curve of polynomial curvature in green, and curve of trigonometric curvature in red.

Figure 8: Examples of geometric systems consisting of two circular arcs with consistent curvature (brown line) connected with different transition curves, for negative values of curvature (using different horizontal and vertical scales): curve of linear curvature in violet line, curve of polynomial curvature in green, and curve of trigonometric curvature in red.

Figures 9 and 10 illustrate geometric systems made up of two inverse circular arcs connected with various types of transition curves.

Figure 9: Examples of geometric systems consisting of two inverse circular arcs (brown line), connected with different types of transition curves, for positive values of input curvature (using different horizontal and vertical scales): curve of linear curvature in violet line, curve of polynomial curvature in green, and curve of trigonometric curvature in red.

Figure 10: Examples of geometric systems consisting of two inverse circular arcs (brown line), connected with different types of transition curves, for negative values of input curvature (using different horizontal and vertical scales): curve of linear curvature in violet line, curve of polynomial curvature in green, and curve of trigonometric curvature in red.

Here the layout of the output circular arc is much more diversified. The position of the arc obtained by the application of the transition curve of linear curvature diverts by several dozen meters from situations where curves of polynomial and trigonometric curvature have been used. The difference for curves of polynomial and trigonometric curvature attains, in the numerical example analyzed, an order of several meters.

The presented examples are taken, in the very nature of things, at random. They reveal wide possibilities for the presented analytical method in modelling geometric systems. Its application makes it possible to generate diversified layout of the communication routes adapted to the field conditions and other circumstances. Simultaneously the method ensures correctness of the obtained solution in view of the kinematic parameters magnitudes in force.

7. Conclusions

The presented universal method of modelling geometrical system of a communication route is based on the determination of an adequate curvature equation by using differential equations. The technique enables us to design connections of various geometrical elements of the route indicating diversification in its curvature. It is characterized by an analytical way of recording certain functions, and moreover it is possible to specify in advance the class of these functions.

The sought-after connection is recorded in the form of parametric equations and , where parameter is the position of a given point along the length of the curve. The determination of and needs the expansion of integrand functions into Maclaurin series and with regard to inverse curvatures into Taylor series.

The determination of coordinates of an adequate transition curve in rectangular coordinate system , makes it possible for an easy transfer of it to Gauss-Kruger conformal coordinates and its subsequent layout in field making use of the satellite measuring technique GNSS.

The method can be applied to identify various types of transition curves. The paper gives examples of solving the problem involving linear change of curvature and the curvature in polynomial and trigonometric form.

The performed comparative analysis reveals a wide range of possibilities which are offered by the presented design method. Its application enables us to generate diversified layout of the route. It simultaneously ensures a complete control of correctness of the obtained solutions in view of satisfying the geometric and kinematic conditions.

The presented method of modelling curvature has a universal character. It can be used for vehicular roads and railway lines, as well. In the case of railway lines there arises an additional possibility of modelling, in a similar way, the superelevation ramp understood here as a determined difference in the height of the rail courses.

Conflict of Interests

The author declares that he has no conflict of interests regarding the publication of this paper.

References

W. Koc and C. Specht, “Application of the Polish active GNSS geodetic network for surveying and design of the railroad,” in Proceedings of the 1st International Conference on Road and Rail Infrastructure (CETRA '10), pp. 757–762, May 2010.

P. Lipar, M. Lakner, T. Maher, and M. Žura, “Estimation of road centerline curvature from raw GPS data,” Baltic Journal of Road and Bridge Engineering, vol. 6, no. 3, pp. 163–168, 2011.

L. Tasci and N. Kuloglu, “Investigation of a new transition curve,” The Baltic Journal of Road and Bridge Engineering, vol. 6, no. 1, pp. 23–29, 2011.

H.-H. Cai and G.-J. Wang, “A new method in highway route design: joining circular arcs by a single C-Bézier curve with shape parameter,” Journal of Zhejiang University: Science A, vol. 10, no. 4, pp. 562–569, 2009.

S. Dimulyo, Z. Habib, and M. Sakai, “Fair cubic transition between two circles with one circle inside or tangent to the other,” Numerical Algorithms, vol. 51, no. 4, pp. 461–476, 2009.

G. Farin, “Class A Bézier curves,” Computer Aided Geometric Design, vol. 23, no. 7, pp. 573–581, 2006.

Z. Habib and M. Sakai, “G² pythagorean hodograph quintic transition between two circles with shape control,” Computer Aided Geometric Design, vol. 24, no. 5, pp. 252–266, 2007.

G. Harary and A. Tal, “The natural 3D spiral,” Computer Graphics Forum, vol. 30, no. 2, pp. 237–246, 2011.

K. T. Miura, “A general equation of aesthetic curves and its self-affinity,” Computer-Aided Design and Applications, vol. 3, no. 1–4, pp. 457–464, 2006.

S. H. Yahaya, M. S. Salleh, and J. M. Ali, “Spur gear design with an S-shaped transition curve application using mathematica and CAD tools,” in Proceedings of the International Conference on Computer Technology and Development (ICCTD '09), pp. 273–276, Dubai, United Arab Emirates, November 2009.

N. Yoshida and T. Saito, “Quasi-aesthetic curves in rational cubic Bézier forms,” Computer-Aided Design and Applications, vol. 4, no. 1–6, pp. 477–486, 2007.

R. Ziatdinov, N. Yoshida, and T.-W. Kim, “Analytic parametric equations of log-aesthetic curves in terms of incomplete gamma functions,” Computer Aided Geometric Design, vol. 29, no. 2, pp. 129–140, 2012.

R. J. Grabowski, “Smooth curvilinear transitions along vehicular roads and railway lines,” Science Papers of AGH University of Science and Technology, no. 82, Cracow, Poland, 1984 (Polish)

L. T. Klauder, “Track transition curve geometry based on gegenbauer polynomials,” in Proceedings of the 6th International Conference on Railway Engineering, May 2003.

W. Koc, Elements of Design Theory on Track Systems, Gdansk University of Technology, Gdansk, Poland, 2004, (Polish).

E. Mieloszyk and W. Koc, “General dynamic method for determining transition curve equations,” Rail International—Schienen der Welt, vol. 10, pp. 32–40, 1991.

B. Kufver, “Realigning railway in track renewals—linear versus S-shaped superelevation ramps,” in Proceedings of the 2nd International Conference on Railway Engineering, May 1999.

W. Koc and E. Mieloszyk, “Comparative analysis of selected transition curves by the use of dynamic model,” Archives of Civil Engineering, vol. 33, no. 2, pp. 239–261, 1987 (Polish).

S. Gleason and D. Gebre-Egziabher, GNSS Applications and Methods, Artech House Publishers, Norwood, NJ, USA, 2009.

B. Hoffmann-Wellenhof, H. Lichtenegger, and E. Wasle, GNSS Global Navigation Satellite Systems—GPS, GLONASS, Galileo & More, Springer, Wien, Austria, 2008.

P. Misra and P. Enge, Global Positioning System: Signals, Measurements, and Performance, Ganga-Jamuna Press, Lincoln, Mass, USA, 2012.

C. Specht, “Accuracy and coverage of the modernized Polish Maritime differential GPS system,” Advances in Space Research, vol. 47, no. 2, pp. 221–228, 2011.

W. Koc and K. Palikowska, “Dynamic properties evaluation of the selected methods of joining route segments with different curvature,” Technika Transportu Szynowego, vol. 9, pp. 1785–1807, 2012 (Polish).