Introduction

According to the World Health Organization (WHO), one million people are killed each year, three million people are seriously injured for life and thirty million are injured in road accidents. Road accident is expected to be the third significant adding to the global burden of disease and injury in 2020. Furthermore, road accidents in Malaysia are one of the formidable problems when addressing this country's road traffic. Malaysia is also one of the Asia's developing countries suffering from high road accidents which are around 30,000 accidents recorded and for every year 1% of its fatal accidents. Vehicle crashes in Malaysia have risen over the last 30 years with a mean incidence of about 9.7% once a year, consistent with data from the Royal Malaysia Police (PDRM). The vast amount of traffic on the road will lead to collisions involving various types of vehicles [1], the most probable collision between motorcyclists and vehicles. A motorcyclist is 17 times more at risk than a passenger driver; therefore, the driver will be more at risk of injury or death because the chance of injury is greater than that of a passenger car. In addition, the probability of injury during an accident is estimated to be 12 times greater than that of a passenger vehicle [15]. Research shows that three major types of accidents are passenger vehicle collisions, other vehicle crashes and single-motorcycle injuries in Malaysia. Several cases of cars and motorcycles accidents, most of the injured were passengers, not drivers or motorcyclists where 97% of motorcyclists were killed or injured compared to 50.5% of the drivers. Moreover, 95.4% of motorcyclists were involved in high-speed motorcycles (single vehicle accidents), and  0.9% are likely to be injured and collieded with cars [5].

Some reports confirmed that most road accidents occur because of driver fault and mechanical failures such as breakdown and split of the tire. Yet, due to the road conditions, the high number of road accidents occur [1]. The road accident factor continues to increase in Malaysia due to the lack of geometric road design accuracy where drivers make mistakes on lane geometric features. There is a lack of continuity when driving along the geometric features of the road, which will increase the chance of accident [3, 6].

85th percentile of operating speed model

Several models were developed over the past 50 years to estimate the operating speed of vehicles based on traffic flow and geometric characteristics [2]. Several studies have been done by using several types of tools to measure the 85th percentile of speed at the two-lane rural roads along the horizontal curves and tangents, which were based on speed limits and types of vehicles. Most of the models produced are supported spot speed data obtained at the actual locations; therefore, the speed calculation used may be a radar gun that can provide less accuracy because the angle of reading is different from zero and drivers can drive (reduce the speed) inconsistently.

The study was carried out in some studies, which calculated the operating velocity profile from the circular curve radius, the angles of deflection and the relation of curve length. Previous research works have considered only passenger vehicles operating under free-flow conditions, while many previous studies specified curve radius because of the impact parameter on the estimation of operating speeds [8]. In addition, some studies have established operating speed prediction models by supporting continuous operating speed information collected through VBox.

Design consistency model

The consistency of design is calculated by estimating the variance of speed along successive geometric elements and assessing the deviation of speed from the standard velocity in a geometric feature section of the lane. This could cause drivers to make errors in speed that result in an increased risk of collision [7]. A road with the best quality level will influence the driving force in a safe and careful way to drive, while high speed can be affected by low consistency along a specific path between different drivers and accidents will end [13]. The features of cars, drivers and traffic affect the geometric designs. The volume and width of lanes, the existence and width of shoulders and highway medians are elements of geometric design that affect traffic operations. Therefore, the horizontal and vertical alignment of the movement of vehicles' highways and conflict streams could lead to accidents [14]. The main criteria used to evaluate design consistency are referred to as the 85th percentile vehicle speed evaluation of operating speed, and are obtained through the utilization of predictive models of operating speed. Researchers have expanded in defining the use of continuous speed profiles (area of normalized bounded by altitude). The area is used to assess the variance of vehicle speed along the road sections which are used to determine the profile design consistency for the whole road sections [4, 7, 9]. Methods of consistency assessment are developed by estimating the operating speed profile utilized in other countries such as the United States, Switzerland and Germany. In addition to the previous consistency approach [7], the following new approaches take into account the consistency of road design. The proposed model represents the parameters as the bounded area between the profiles that are based on the average speed along a two-lane highway segment. It can be seen while the consistency of the design increases, the crash rates will show a significant decrease. The equations correspond with a function of Ra (1)

$$R_{{\text{a}}} = \frac{{\left( {\sum a_{i} } \right)}}{L}$$
(1)

\(R_{{\text{a}}}\) is the relative area in m/s measuring consistency, \(\sum a_{i}\) is the total of i bounded area between profiles of speed and operating speed (m2/s), L is the segment length (m) in total

$$\sigma \, = \left\{ {\frac{{( V_{j} - V_{{{\text{avg}}}} }}{n}} \right\}^{0.5}$$
(2)

where σ is the standard deviation of vehicle operating speeds (km/h), \(V_{j }\) is the operating speed of the jth geometric element in tangent/curve (km/h), \(V_{{{\text{avg}}}}\) is the length average weighted of operating speed at a segment (km/h), n is the geometric elements along a road segment (km/h)

$$V_{{{\text{avg}}}} = \frac{{\sum V_{i} L_{i} }}{{\sum L_{i} }}$$
(3)

where \(V_{{{\text{avg}}}}\)   \(V_{i}\) is the operating speed (km/h), and \(L_{i}\) is the length of the segment (km).

The model of consistency (C) and the model integrated design consistency are as follows

$$C = 2.808.\exp^{{\left( { - 0.278.R_{{\text{a}}} .\sigma } \right)}}$$
(4)

C is the basic consistency of a road segment, \(R_{{\text{a}}}\) is the normalized bounded area based on the average speed profile of cars–motorcycles and the average operating speed (m/s), σ is the standard deviation of car speeds (m/s).

$${\text{IC}} = \left[ {2.808*\exp^{{\left( { - 0.278*R_{{\text{a}}} *\sigma } \right)}} } \right]* \exp^{{\left( { - 0.01* A_{{{\text{CT}}}} } \right)}}$$
(5)

where IC is the integrated consistency of road sections, σ is the standard deviation of car speeds and motorcycles, and ACT is the normalized bounded area between the speed profile of cars and motorcycles (m/s).

Data collection

The objective of this research is determining and deriving the performance level of road safety based on the development of linear speed profiles of cars–motorcycles. Based on vehicle size, location, vehicle data velocity using the Global Positioning System (GPS) and road geometric characteristics, the data collection was calculated. The study was conducted on the Johor Federal Road (F0050) from 21 to 25 km at Batu Pahat Community Nurses College in Taman Manis, Parit Raja, Johor as shown in Fig. 1.

Fig. 1
figure 1

Location 21–25 km (F0050)

During the daytime, the primary speed data collection was carried out under free flow and off peak, along the selected segment with the GPS DG-200 for vehicle tracking movements. GPS was activated for moving vehicles: cars and motorcycles during which the two vehicles were used (two cars, two motorcycles). A hundred (100) test drivers and riders were instructed to use the same standards for vehicles as drivers and riders driving along the sections selected. In order to avoid the effect of vehicle speed shifts and free flow velocity reduction, speed data from the signalized intersection were omitted at approximately 300 m space distance [16]. The GPS system has been designed to identify the vehicle locations and speeds at intervals of time in 1 s. The data collected required two surveyors, drivers and GPS trackers. For each segment/road section, the average vehicle travel speeds along the traveled space are determined by GPS [10,11,12]. In addition, a replacement method was used to generate the continuous speed profile information. Therefore, the index of integrated road design consistency (IC) can be determined based on the speed profiles of cars–motorcycles.

There are secondary sets that provide data from the incident report, traffic volumes and traffic characteristics. Data collection of accident data from the police (PDRM) was obtained, which involves the type of accident, location and number of vehicles involved. Average annual daily traffic (AADT) collected from agencies such as the Highway Design Unit (HPU) is included in the collection of traffic data. Road characteristics include features of the road and direction, such as curve length, tangent length, lane width, shoulder width and entry point. In order to determine the accuracy of the road design model, information was obtained from authorities such as Public Ministry Works (JKR) or by visiting the locations [7].

Result and discussions

Model of tangent and curve 85th percentile based on operating speed

Operating speed models are developed based on the current method of analyzing the vehicles' speed profiles. The GPS model of DG-200 was used to collect the speed-time data which are used to plot the profiles. Furthermore, vehicles' speed and road features were obtained (on site) to produce the 85th percentile of speed. The multiple linear regression approach is used to use the statistical software (Minitab) to construct the 85th percentile speed prediction model. The tangent and curve sections were also considered to estimate the operating speed model for motorcycles and cars. The equation is developed using different variables: length of tangent, curve length, shoulder width, lane width and point of access. The subsequent model shows that the amount of vehicles/volumes of traffic (AADT) did not contribute to the V85 speed. The statistical findings show typical linear functions.

Motorcycles-85th percentile models:

$$V_{85t} = \, 72.0 \, + \, 6.99{\text{ SW }} + \, 11.7 + \, 0.276{\text{ AP }}{-} \, 2.61{\text{ LW}}\quad R^{2} = \, 0.98$$
(6)
$$V_{85c} = \, 20.6 \, + \, 4.16{\text{ SW }} + \, 102 + \, 0.82{\text{ AP}}\quad R^{2} = 0.79$$
(7)

Cars-85th percentile models:

$$V_{85t} = \, 88.9 \, + \, 27.4{\text{ SW }} - \, 1.26 + \, 1.07{\text{ AP }} - \, 16.9{\text{ LW}}\quad R^{2} = \, 0.99$$
(8)
$$V_{85c} = \, - \, 31.1 \, + \, 4.1{\text{ SW }}202 + \, 1.38{\text{ AP}}\quad R^{2} = \, 0.83$$
(9)

\(V_{85t}\) is the tangent speed (operating speed), \(V_{85c}\) is the curve speed (operating speed), \(L_{{\text{t}}}\) is the tangent section (length), \(L_{{\text{c}}}\) = curve sections (length), SW = shoulder width, LW = lane width, and AP = road access point.

Validation of linear profile models

Validation of the models developed is performed to validate the models' prediction and accuracy. The data of speed observed and the operating speed model were plotted  which results in well compared between data and the model. Equations (6) and (8) show the negative (−) effects of the LW to the vehicles' speed behaviors. This could be explained as the typical road geometric features F0050 road with high activities along the edge of a major road. Vehicles tend to slow down and reduce the speed while driving at the wider major road with roadside activities such as shops and illegal parking. However, the study has found that AP would contribute for the higher speed (+). This is true that vehicles did not reduce the speed while driving at the major road F0050 even though several access points/minor roads are found along the side of the major road. Therefore, Eqs. (69) show a high number of R2 which means of high and homogenous data and correlation among the variables; V85, L, SW, LW and AP. The figures of the linear line fit to show well compared of observed and predicted linear speed model, which are shown in Figs. 2 and 3. Comparisons of observed velocity versus predicted motorcycle speed at 21–25 km are shown in Fig. 2, and comparisons of observed velocity versus predicted car speed at 21–25 km are shown in Fig. 3.

Fig. 2
figure 2

Comparison observed speed versus predicted speed of motorcycles (21–25 km)

Fig. 3
figure 3

Comparison observed speed vs. predicted speed of cars (21–25 km)

Development design consistency

The consistency of the design is verified by using the models of continuous speed profile operations. ACAD software was used to develop a road design consistency model. Plots of the model operating speed profile (85th percentile) along the length are called as speed profiles and are used in this study. The graph line represents the early-developed models of the operating speed along the jth geometric elements that represented tangent and curve of roads. The motorcycle speed profile of F0050 21–25 km is shown in Fig. 4, and the car speed profile of 21–25 km is shown in Fig. 5. Using the 85th percentile of the operating speed model, the continuous operating speed profile models were plotted. The linear lines pattern shows the motorcycles–cars speed along the segment of roads that depend on the road environmental factors: traffic signal location, access point and intersection along segments. Furthermore, architectural factors were also considered effects on motorcycle and vehicle speeds: shoulder width, lane width, and median and access point.

Fig. 4
figure 4

Speed profile F0050 of motorcycles (21–25 km)

Fig. 5
figure 5

Speed profile F0050 of cars (21–25 km)

The following model is developed based on the speed profiles establishment which relies on two main parameters: speed and road geometric features. The design consistency of every road sections is based on Eqs. (1), (2) and (3). The equations are assessed for every road segment per km. Furthermore, the bounded area between the profiles, the vehicles' average speed and the standard deviation of speeds along a two-lane highway segment was defined. A consistency model was created based on Eq. (4). Table 1 displays motorcycles profile performances at 21–25 km per km while cars at 21–25 km per km are shown in Table 2.

Table 1 Design consistency motorcycles (21–25 km) per km
Table 2 Design consistency cars (21–25 km) per km

Integrated design consistency model

An integrated-consistency model is required in developing the function of car and motorcycle speed profiles. The model is included with the effect of motorcycles speed profiles on design performances. Figure 6 shows an F0050 speed profile at 21–25 km between motorcycles and cars. The normalized bounded area between the motorcycles and cars speed profiles is calculated by means of the speed profiles. Table 3 demonstrates an integrated feature consistency at 21–25 km per km between motorcycles and cars. Each kilometer road section shows value-integrated consistency in design and their threshold per km. The accuracy of integrated design consistency is evaluated using Eq. (5).

Fig. 6
figure 6

Speed profile F0050 of motorcycles and cars (21–25 km)

Table 3 Integrated design consistency motorcycle and cars (21–25 km) per km

The speed consistency model was developed for every kilometer (21–25 km). The consistency model along 21 km of motorcycles design consistency shows an acceptable design and good performances for cars. Motorcycles–cars integrated speed consistency shows good performances along the road sections surrounding areas with mosque, restaurants. Intersections and access points with no median/divider were found along the sections which result in conflict movements and accidents on the road. The bounded area of speed consistency motorcycles–cars gives an acceptable performance at 22 km and the integrated consistency cars–motorcycles perform suitable profiles. The study found suitable consistency along the road Ayer Hitam–Batu Pahat even though the road does not have a median that can cause vehicle conflicts, which contribute to the accidents.

Results on IC bounded profiles at 23 and 24 km for cars–motorcycles show poor design consistency while the integrated design consistency also shows poor design for both types of vehicles. This can be concluded based on areas surrounding that the road geometric design considered are not designed properly which is not suitable for traffic flow streams; the roads are located close to several road activities such as shop and illegal parking at road shoulders. Car drivers and motorcyclists commonly drive at the high speed along the federal roads, which can be derived with good performances of integrated speed consistency for cars–motorcycles. However, performances decrease as cars drive through 25 km with poor consistency same as the integrated design consistency cars–motorcycles. Vehicles tend to drive with higher speed along the section of 25 km, and cars–motorcycles travel with comfort due to straight road section which mostly results in a high accident rate.

Conclusions

The following study has conducted the analysis of speed design consistency at federal road F0050 in Johor. The following study of the integrated design consistency between the motorcycles and cars along  km 21 up to km 25 indicates that  poor road design consistency would likely results in more accident occurs. The study found that there are several access points along the section along a road ( km21–km25) and mostly vehicles coming from minor roads cross the main road without hesitation, which creates the potential of more conflict traffic streams. The conflicts would results in number of accidents on road traffic movements with particularly motorcycles and cars. The un-availability of road dividers (two-way directions) along the roads could be the worst conditions, in which drivers tend to cross the road from a different direction without any control/no separations.

The standard of road design would impress the driving force and behaviors of drivers and motorcyclist to drive along the road safely. Therefore, road geometric design performance level for suitable vehicles-motorcycle movements should be designed with results in sophisticated integrated design consistency between cars–motorcycles where good, fair/reasonable and poor design consistencies are specified as the road performance level.