Abstract

Safety is an important aspect of road design. In highway geometric design, road engineers use a two-stage design method due to its convenience, but the available two-dimensional design tools on the market are believed to impose limitations. Horizontal and vertical alignment parameters are considered and determined at the designer’s discretion, potentially downplaying the three-dimensional characters of spatial curves. This study focuses on the three-dimensional characteristics of highway alignments and investigates a safety evaluation method to establish the relationship between crash rate and spatial curve properties. This will necessitate an analytical investigation into the effects of higher-order properties, such as curvature and torsion, on the geometry of spatial Cartesian curves. First, these combinations of horizontal and vertical alignments were categorized into six classifications, each with its own spatial expression in mathematical form. After manipulating the curvature and torsion of the spatial curve algebraically, the correlation between geometric design variables and crash rate was ultimately established. A few cases involving geometric design data and crash facts were utilized for verification. The results revealed a considerable positive correlation between curvature or torsion variance and crashes per million vehicles kilometers, as a slight difference between curvature and torsion could also be spotted. And curvature distribution is correlated with collision frequency more closely than torsion spatial variation.

1. Introduction

Highway geometric design should not only meet the needs of traffic but also consider social, economic, environmental, and other factors. The composition of a road’s spatial alignment is a comprehensive product of an integrated part of an integrated transportation process. Since the 1930s, national road design standards, policy studies, and guidelines have been implemented to address economic, environmental, and other social issues [1, 2]. However, more and more people are affected by man-made road disasters, making road risk control crucial [3]. Among them, the safety of highway alignment design is the primary consideration, and safety is also a vital factor in road design that may have an impact on other social compositions [4]. According to the data of the World Health Organization in 2022, approximately 1.3 million people around the world have their lives terminated due to road traffic accidents every year. The concept of flexibility design believes that there are many road accidents related to the unreasonable design of highway alignment [5, 6]. This makes it necessary for designers to re-examine the importance of road design. It is likely that the economy and safety of the highway construction are in conflict, but a reasonable safety design could be the best compromise [6]. For the purpose of sustainability, research should place a greater emphasis on the consistency of dynamic route design from the perspective of improved personal and property safety. However, the traditional two-dimensional evaluation methods may not truly reflect the real situation of road alignment design.

Horizontal alignment is of paramount significance to location in highway geometric design, while vertical alignment is restricted by altitude and has great impacts on project costs, including earthwork, construction, land use, and user costs [5]. Due to the complexity of three-dimensional (3D) spatial curve design, a two-stage method is typically adopted to facilitate the combinations of horizontal and vertical alignments (spatial curve) by completing the horizontal location first and then the vertical design [7, 8]. In the two stages, quantitative reference values are assigned to design elements as control standards to meet fundamental safety and comfort requirements. However, for horizontal alignment and profile combinations, previous studies appeared to be limited to qualitative analysis.

3D highway alignment design and safety evaluation are a persistent problem for which there are no always-satisfactory implementation methods. How to comprehensively consider the coordination of horizontal and vertical alignments is challenging. Several studies have begun to contemplate and investigate combined highway alignment design [911]. To optimise highway design, Wang et al. [7] and Zhou et al. [12] created a method for constructing a 3D model of highway alignments using parameter programming. The use of visualization technology was an alternative way of creating a 3D model of highway alignments. It was believed to provide an opportunity for designers to calculate virtual perspective views and control specific parameters with visualization tools to verify the 3D alignment [13]. In addition, spatial curve fitting to highway alignments was supposed to be beneficial to the better design of the next generation of highways and future research on vehicle infrastructure integration, route guidance, and congestion management.

Although a route is a spatial curve with randomness in design, it has certain regularity consisting primarily of six types of horizontal and vertical elements combinations (including tangent segment, circular curve, and spiral curve of horizontal alignment and tangent segment and parabola curve of vertical alignment). This paper aims to build a parametric representation of highway geometric alignments based on the mathematical properties of curves in order to evaluate the interrelations between crash rate and spatial path planning.

The remaining sections of the paper are structured as follows. Literature review is presented in Section 2, discussing evaluation methods for highway geometric, limitations of 2D methods, and contributions of 3D models. In Section 3, spatial curve properties used in highway design and 3D geometric models of spatial combinations are introduced. Meanwhile, the third section provides an evaluation method for design indexes. In Section 4, the application and discussion of the new approach are illustrated through case studies, and analysis results are provided. Conclusions and significant findings are summarized in the final section.

2. Literature Review

Road design research has made significant progress in recent decades. Advanced computer simulations and new materials have enabled improvements in road quality and safety features. In the early stages of highway geometric design, it may be practical and feasible to measure and evaluate safety by comparing computations to design policy. When the design parameters are fulfilled, the road is considered to be safe and reliable. Interactive highway safety design model (IHSDM) is a typical application that provides a policy review module to compare roadway segment geometry to pertinent design policy [14]. It was discovered that the evaluation of highway geometric design consistency checks is an important issue in design and traffic operation evaluation [15]. Because only considering separated indexes that meet design standards does not guarantee the alignment fluidity and safety of the highway, the coordination research on horizontal and vertical alignments is a critical procedure that should commence with a preliminary design in order to be easily adaptable. IHSDM also supports the evaluation module for design consistency [16]. The review shows that operating speed consistency [17, 18], traffic volume and crash consistency [19, 20], vehicle stability and driver workload consistency [21], and design consistency of highway elements as well as sight distance consistency [22, 23] were considered to be the research focus.

However, there are still key issues that need to be addressed. Road accidents remain a significant global public health concern. Existing road standards and design practices do not adequately account for the complex physics. Previous studies focused more on qualitative analysis by external means, such as engineering experience, a two-dimensional (2D) analysis model, and polynomial or spline fitting [24]. In essence, highway alignments are 3D spatial curves but horizontal and vertical alignments are treated differently in 2D projections. Thus, the evaluation of geometric continuity is calculated in two planes, respectively. It cannot precisely reflect the characteristics of spatial alignment continuity in 3D Euclidean space [25]. Moreover, research studies in 2D lack explicit quantitative coordination of combined alignments [25].

3D alignment models provide a more comprehensive understanding of the driving experience. They enable designers to visualize clear-sight triangles, simulate curves and grades, identify concealed dangers, and analyze collisions from multiple angles [26]. Studies have demonstrated that 3D alignment design significantly improves safety over 2D methods [23]. For example, 3D models optimise sight distance on hills and around curves based on the eye height of the driver more effectively [27]. The occurrence of collisions can be reduced by incorporating appropriate gradients and curves on crucial sections of roadway [28]. Alternately, 3D alignment design aids in reducing the travel time of motorists and releasing the loss caused by traffic congestion [29].

3. Methodology

3.1. Basic Properties of a Spatial Curve

In highway geometric design, five line-type elements including tangent segment, transition curves (Euler spiral), circular curves in the horizontal plane (3 types), and tangent and parabola segments in the vertical plane (2 types), are used to build the complete path by permutations and combinations in the space coordinate system. Different combinations of horizontal and vertical elements will generate different parametric equations. Six types of combinations are shown in Table 1, and their space shapes are descripted [30].

The coordinate points of the spatial curve (see Figure 1) can be represented by the following parametric equation:where represents arc length (corresponding to stake number of highway).

3.1.1. Tangent

The parametric equation of a straight line in two-dimensional Cartesian space is shown in the following equation:where is a starting point for calculating in the straight line; is the distance between the calculated point and the point , m; is an angle between tangent and X-axis, .

If the line is a combination of horizontal and vertical tangent that locates in 3D Cartesian space (), the parametric equation can be written in the following equation:where is a starting point for calculating in and denotes the slope of spatial straight line.

3.1.2. Spiral Transition Curve

Generally, the Euler spiral is used in the design of spiral transition curves [31]. The curvature radius varies from infinite at the tangent end of the spiral to the radius of the circular arc at the end that adjoins the circular arc (). The radius of curvature at any point on an Euler spiral varies inversely with the distance measured along the spiral.

Defining as the length of a spiral curve beginning with a straight segment, as the corresponding radius, and as the deflection angle of this segment in radians, Euler’s spiral is formulated as follows:where is the spiral parameter; is the segment length of spiral curve; and is a parameter of the spiral equation that denotes the azimuth angle of the located point along the tangent direction.

Referring to the fundamental theorem of calculus, the geometrical relationship of Euler spiral parameters can be expressed as follows:

The parametric equation of the spiral curve is shown in equation (6) through definite integral to in the interval .

3.1.3. Circular Curve

The most common type of horizontal curves used to connect tangent segments of highways is the circular curve. The parametric equation of circular curve is written in the following equation:where denotes the center of a circle; represents the radius of circular curve; and is a parameter of circle equation that denotes the azimuth angle.

3.1.4. Parabola

The vertical curve is the transition section between the deflection points connecting adjacent slopes. Building vertical curves can ensure the smooth connection of geometric alignments and adequate visibility for driving safety. Generally, a quadratic parabola is used to complete the transition.

The equation of a parabola curve can be expressed as , where represents the curve length of a point from the basic point on the parabola. The tangent line of a parabola can be written as . Supposing the starting point of the parabola curve is , and denote the slope of tangents before and after the parabola, respectively. represents the total length of the parabola. Thus, two parameters of the equation can be solved to gain and . The equation of parabola curve is shown in the following equation:

3.2. 3D Highway Geometric Model

According to the selection and explanation of vertical alignments and shape descriptions in Section 3.1, six different combinations are defined in this study (Table 2), where TT represents a spatial combination alignment consisting of straight line in the horizontal projection plane and straight line in the vertical projection plan (Figure 2), TP denotes a spatial combination alignment consisting of straight line in the horizontal projection plane and parabola in the vertical projection plane, ST means a spatial combination alignment consisting of Euler spiral in the horizontal projection plane and straight line in the vertical projection plane, SP means a spatial combination alignment consisting of Euler spiral in the horizontal projection plane and parabola in the vertical projection plane, CT signifies a spatial combination alignment consisting of circular curve in the horizontal projection plane and straight line in the vertical projection plane, and CP signifies a spatial combination alignment consisting of circular curve in the horizontal projection plane and parabola in the vertical projection plane.

Figure 2 shows 6 type combinations in 3D, as already stated in Table 2. Based on the 6 types of line element combinations, any point coordinate in the space along the highway center line can be determined by algebraic representation. For example, the equations of a horizontal spiral or circular curve with a vertical parabola (type SP or CP) can be expressed as formula (9). Besides, the different manifestations of line segments on different planes (xy and zy planes) are defined into four types (T, S, C, and P) to describe the changes in the curvature of a line segment. For example, two S curves were connected by a C curve, which could be considered as a part of the circle. Then, a curve segment could translate smoothly into a straight segment like the T segment shown.

3.3. Design Indexes
3.3.1. Curvature

In a spatial parametric curve, the computational formula of curvature can be written as shown in equation (10) [32]. The terms in the equation represent that the computing result is a spatial curvature. If the terms were removed, the equation would be reduced to planar parametric curves. In a planar curve, the computing value denotes the local reciprocal of curvature and generates an indication of the bending in the curve. Thus, a zero value represents a straight line.where and .

Algebraic operation shows that, for a parametric curve, the curvature value in the equation will always be positive as opposed to being a signed value as in the planar case. It is believed that a spatial curve can technically bend in an infinite number of directions, as opposed to just two directions in a planar curve [33]. As the curve is traversed, the tangent vector can either move in the clockwise or counterclockwise directions, and these two directions respectively correspond to positive and negative curvature. However, in a spatial curve, the tangent vector can rotate in any direction. Thus, curvature is defined as the magnitude of the bending without a direction, and an inflection point in a spatial curve cannot be described as a point where the sign of curvature changes.

Geometrically, curvature is a measure of how quickly the unit tangent is moving with respect to distance along the curve direction. Physically, it represents “bending” in the space of a curve. The reciprocal of curvature is called the radius of curvature denoted by

March [32] argued that the curvature is the magnitude of the change in direction of the tangent along the curve. When the curvature is zero, the tangent vector will not change, and the curve will continue in a straight line. On the other hand, a high value of curvature will result in the rapid direction changing of the tangent vector. An infinite curvature, as in the case of a cusp, represents the discontinuity in the tangent vector along the curve.

3.3.2. Torsion

Torsion is another property of spatial curves that is used to measure the tendency of a curve to twist out of the plane. For a curve in parametric, the computational formula is shown in equation (12) [32], where . Obviously, a curve in the planar will have a zero value torsion. Because the coordinate value of the terms is zero, thus there will be a zero from first-order derivative to three-order derivative in the direction. One thing to notice is that, unlike curvature, the torsion of a curve in space has a signed value (see Figure 3).

Physically, torsion is a measurement of the changing rate of the osculating plane ( plane) relative to the governing parameter (corresponding to stake number location). The osculating plane is defined as the plane spanned by the curve tangent and normal vectors, namely, and . Therefore, a constant zero torsion denotes a planar curve lying in the osculating plane.

The parametric representation of each spatial curve segment is built in Table 3, and their properties are calculated based on curvature and torsion. This representation provides valuable geometric meaning to the spatial shape of the curve, as the physical meaning of curvature and torsion is well understood.

3.4. Evaluation Model

In order to determine the interrelations of highway geometric alignments and crash rates, the selected highway is divided into several segments in terms of spatial curve combinations. Related characteristic functions revealing the influence of curvature/torsion differences on crashes per million vehicles-kilometers are used to evaluate the quality of highway geometric design.

3.4.1. Local Evaluation

The curvature/torsion difference of adjacent segments is shown in the following equation:where denotes the curvature/torsion difference of adjacent combination segments and ; represents the average value of curvature/torsion in section ; denotes crash rate of the section ; and means that there is a certain correlation between and .

Crashes per million vehicles-kilometers are calculated to represent the crash rate. The results show each combination segment has a corresponding value of the crash rate ( to , where denotes curvature/torsion difference of adjacent section and ). Thus, the rationality of alignment coordination in design can be determined by the contrast analysis.

3.4.2. Global Evaluation

For the whole highway, the arithmetic mean values of spatial curvature/torsion are calculated as follows:

Then, the standard deviation of spatial curvature/torsion is shown in the following equation:

To evaluate the dispersion of a probability distribution of each spatial curve property (curvature/torsion), the coefficient of variation is introduced, as shown in the following equation:

The coordination of horizontal and vertical alignments, through the above calculation, can be described and adjusted to cut down the crash rate.

4. Results and Discussion

The design documents and nearly five years of crash data from 21 roads were collected for study analysis, including 6 expressways and 15 ordinary highways in China. Each road was orderly divided into many spatial segments in terms of 6 types of combination rules (see Figure 4). The relationship between crash rate and spatial curvature and torsion was analyzed through statistics and regression analysis. The applications of local and global evaluation in the form of case studies are shown as follows, respectively.

The total length of the selected roads is 1072.07 kilometers. All roads were divided into 1345 space curve segments based on the abovementioned 6 types of combination rules. According to the analysis of the crash data, some accidents are not relevant to the road alignment changing or the correlation is weak and are excluded from the statistics, such as drunk driving or misoperation. In addition, some nonstatistical accidents are also removed to ensure the correlation between accidents and space alignments. The detailed analysis results are described as follows. The effect analysis of crash rate and curvature difference is shown in Figure 5. And the effect analysis of crash rate and torsion difference is shown in Figure 6.

The statistical results showed there is a positive correlation between curvature/torsion difference and crashes per million vehicles-kilometers (Table 4), especially for curvature difference (crash rate is more sensitive to curvature change of adjacent segments). With the increase in the curvature difference, the trend of crashes is ascending notably. And with the increase in curvature and torsion differences, not only the crash rate increases but also the influence of curvature and torsion on the crash rate increases. The results indicate that controlling curvature/torsion differences can reduce crash rates.

From the perspective of 3D space, the crash rate shows a rising trend as curvature and torsion difference increases simultaneously. It is also verified that the changes in physical properties of spatial alignments, curvature, and torsion are significantly positively correlated with the crash rate. However, a crash is a random event in general, especially when displayed in the beginning space within the curvature difference range from 0 to and torsion difference range from 0 to . As shown in Figure 7, those special points can be observed that fail to meet the integral tendency. The analysis shows that although these accidents are affected by the factors of road alignment space changes, other accident factors may play a leading role, such as the driver ’s fatigue driving in a small radius curve section or speeding through the curve.

After highway path construction, the spatial alignment of any road and its alignment combination have been determined. In spatial curve building, each road has a fixed combination of curvature and torsion. The average, standard deviation, and coefficient of variation of curvature/torsion, as well as crash rates can be calculated according to the collected data. The difference in curvature/torsion of each road will correspond to different crashes. Table 5 and Figure 8 reveal the relationship of the crash rate and spatial curve under different road conditions.

Obviously, computed results show crashes have a significant fluctuation following the changes in curvature and torsion. On the whole, there is a positive correlation between crash rate and mathematical statistical information on curvature/torsion (Table 6). The crash rate has a stronger correlation with the AV of curvature than the AV of torsion. The SD of torsion is more important to the accident rate than curvature. The correlation between CV of torsion and crash rate is weaker than that of CV of curvature. Among the parameters related to curvature, the crash rate has the strongest correlation with the AV of curvature, and the CV of curvature has the weakest correlation. Among the relevant parameters of torsion, the crash rate has the strongest correlation with the SD of torsion, and the correlation with CV is the weakest. And the nonlinear relationship between road parameters and crash rate is further estimated by the regression model (Figure 8 and Table 7).

Interrelations between the crash rate and curvature are shown in Figures 8 1(a)–1(c). For curvature, the crash rate increases with the increase in AV, SD, and CV. However, with the increase in curvature, the influence of curvature on the crash rate is somewhat reduced. What can also be noticed is that curvature change has a significant influence on crashes within a certain range, such as average in the intervals [0 to 0.002] and standard deviation in [0 to 0.006], seen from Figures 8 1(b) and 1(c). However, the number of accidents may reduce as the dangerous degrees of the road continue to increase (generally, a greater discrete degree of curvature change indicates more dangerous and difficult road conditions), the reason for which may be that poor road conditions usually require drivers with a high level of caution and a lower operating speed.

Figures 8 2(a)–2(c) are compared with crash rate with the torsion change, which mainly influences the longitudinal gradient in the vertical plane. A larger torsion angle will bring a higher risk of collision, as shown in Figure 8 2(a). Frequent and sharp changes of slope are detrimental to driving safety, with greater torsion bringing a surge in collision risk as shown in Figures 8 2(b) and 2(c). In other words, the higher the discrete degree of torsion change is, the more crashes there will be. However, for the AV of torsion, the greater the AV of torsion, the lower the impact of the AV of torsion on the crash rate. For the SD and CV of torsion, with the increase of SD and CV, not only the accident rate will increase but also the effect of SD and CV on the accident rate will become stronger.

5. Conclusions

The spatial configuration should not be determined separately, as the location of a highway is comprised of horizontal alignment and profile. Curvature and torsion, two invariants of a curve in space, can define the spatial characteristics of highway alignments. This study develops an evaluation method for determining highway geometric design and safety operations based on spatial curve properties. The detailed methods and findings of the research are summarized as follows.(1)The spatial curve path of a road is segmented into six types based on spatial combination regularities. The three-dimensional location of each point on the curve is written by a parametric representation. Thus, three-dimensional numerical models of highway alignments are built.(2)Based on a parametric representation of a spatial curve, two geometric invariants of curvature and torsion are calculated to determine the spatial shape and trends. The dispersion of curvature/torsion change between adjacent segments on a low-grade highway is greater than on a motorway, according to the results of a calculation.(3)Using global and local approximations, the interrelationships between collisions and curvature/torsion distribution are determined. The results of an analysis show that the regularity of curvature/torsion change has significant impacts on crash frequency. By comparison, curvature distribution has more obvious correlation to crashes than torsion spatial variation. The segments with a higher dispersion of curvature/torsion difference correspond to a higher crash rate.(4)The proposed method can be used to identify the spatial structure of curvature/torsion for in-plan or in-service highways, providing a basis for the determination of crash frequency and optimizing highway geometric design to improve safety operation levels. The evaluation model building will allow for improved accessibility for future research to explore the 3D space design of route.

The safety of road alignment design is a primary consideration in this study. The physical properties of spatial alignment, curvature, and torsion are used to characterize the quality of road alignment design, beginning with the actual situation of road alignment space change. This study establishes a safety evaluation method for analysing the interrelation between crash rate and spatial curve properties. It also offers a new perspective for evaluating road alignment design from a three-dimensional state. Then, a genuine three-dimensional evaluation of the road’s safety is realized. However, the design of roads does not account for the equilibrium between all sustainability factors. The social, economic, and environmental aspects of road design should also be examined in depth, despite the fact that safety is a significant factor that directly affects people and indirectly affects other social factors, such as economic loss and social concern. Thus, although safety is a top priority, economic, environmental, and social factors need to be considered more thoroughly in future research to promote urban sustainability.

Notations

2D:Two-dimensional
3D:Three-dimensional
&:Superposition of two line elements
(x, y, z):Any coordinate point of spatial curve
(x0, y0, z0):A starting point for calculating
l:Length, m
i:Slope of spatial straight line, %
α:Angle, degree
R:Corresponding radius, m
β:Deflection angle of this segment in radians, degree
A:Spiral parameter
(m, n):Center of a circle
θ:Parameter of circle equation that denotes the azimuth angle
a, b:Coefficient of a parabolic equation
κ:Curvature, m−1
ρ:Radius of curvature, m
τ:Torsion, m−1
:The curve tangent vector
:The curve normal vector
γ:Changing rate of grade
S:Length of vertical curve, m
i 1, i2:Grades of the former and latter tangents, %
:The curvature/torsion difference of adjacent combination segments, m−1
:The average values of curvature/torsion in section i or i − 1, m−1
:The crash rate of the section, i
:A certain correlation
:Mean value of spatial curvature/torsion, m−1
σ:Standard deviation of spatial curvature/torsion, m−1
CV:Coefficient of variation
AV:Average
SD:Standard deviation
CR:Crashes per million vehicle-kilometers
TT:A spatial combination alignment consisting of straight line in the horizontal projection plane and straight line in the vertical projection plane
TP:A spatial combination alignment consisting of straight line in the horizontal projection plane and parabola in the vertical projection plane
ST:A spatial combination alignment consisting of Euler spiral in the horizontal projection plane and straight line in the vertical projection plane
SP:A spatial combination alignment consisting of Euler spiral in the horizontal projection plane and parabola in the vertical projection plane
CT:A spatial combination alignment consisting of circular curve in the horizontal projection plane and straight line in the vertical projection plane
CP:A spatial combination alignment consisting of circular curve in the horizontal projection plane and parabola in the vertical projection plane.

Data Availability

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research was financially supported by the General Project of Philosophy and Social Science Research in Colleges and Universities in Jiangsu Province (Grant no. 2022SJYB2224) and the National Natural Science Foundation of China (Grant no. 72071097).