Development of a Rapid Inspection Driving Cycle for Battery Electric Vehicles Based on Operational Safety

: The aim of this paper is to solve the problem for battery electric vehicles of low-precision and time-consuming inspection. A novel method of driving cycle development for battery electric vehicles’ operational safety is proposed in this paper. First, three inspection items are proposed based on relevant testing standards. The inspection calculation method of operational safety is developed based on the acceleration changing rate. Then the multi-cycle inspection method with the stable pedal mode is developed, and the Gauss ﬁltering algorithm is applied for data preprocessing. A rapid inspection driving cycle construction method based on support vector machine is proposed, and a driving cycle is built with a total time of 204 s by fusing and splicing kinematic fragments. Finally, the proposed inspection calculation method is used to validate the operational safety inspection items by tracking the established rapid inspection driving cycle based on the test bench. The results shown are those that qualiﬁed the range of acceleration changing rate for driving stability [ − 0.35, − 0.04]. The range for gliding smoothness is [0.05, 0.09]. The range for braking coordination is [ − 0.04, 0.095]. The maximum RMSE between the constructed rapid inspection segments is 9%, and the maximum RMSE between the tested driving segments is 6%. Test results meet design requirements. The thresholds for operational safety inspection items are evaluated based on the test results. We set less than 0.5 as the safety threshold for driving stability. During the experiment, gliding was less than 0.1 as the safety threshold for gliding comfort, and during braking it was less than 0.1 as the safety threshold for vehicle braking coordination.


Introduction
As one of the best approaches for environment failure and air pollution, the new energy vehicle industry has been highly valued and vigorously developed by various countries [1][2][3]. China's reasonable number of new energy vehicles reached 5 million in 2020 and will reach 17 million in 2030. Therefore, the development of new energy vehicles is regarded as a national development strategy [4][5][6].
Currently, in the study of electric vehicle driving range, battery energy management and vehicle operational safety are increasing; the driving cycle is commonly used to inspect electric vehicle driving range, battery states estimation, and energy management [7][8][9]. Several researchers proposed various applicable transient driving cycles in specific areas or on specific vehicles [10]. For example. the US Federal Test Procedure (FTP)-75 tests the fuel economy of vehicles in highway driving [11]. Japan's driving cycle is mainly used to certify vehicle emissions and fuel consumption [12]. Depending on the duration of the driving cycle in different regions, Berzi et al. [13] proposed an electric vehicle driving cycle for detecting energy consumption and efficiency evaluation based on the current traffic development in the city of Florence. Hung et al. [14] developed an optimal 2 of 16 Hong Kong urban driving cycle using vehicle performance values, vehicle velocity, and acceleration distribution to test and evaluate vehicle emission, which was easy to follow and could be substituted for other cycles. Zhao et al. [15,16] proposed a representative system construction method for urban driving cycle for battery electric vehicles in a case study in Xi'an city, which compared the driving range and energy consumption under different driving cycles, and the results showed that this driving cycle can effectively reduce energy consumption. In addition, other researchers have proposed different methods for the construction of the actual driving cycle. For example, Shi et al. [17] chose typical road driving cycles in Hefei, China, and defined 12 characteristic parameters to evaluate the constructed cycle. Shi et al. [18] designed driving cycles in Changchun, China by Markov property and laid a theoretical foundation for designing driving cycles and ECO driving (Economical and Ecological). Jing et al. [19] used linear discriminant analysis to construct a driving cycle in Tianjin, China. A new methodology was offered for building driving cycles and referenced value to related research; Wang et al. [20] analyzed driving characteristics and developed the driving cycles in 11 typical Chinese cities including Beijing and Shanghai. In addition, they proposed some important factors that lead to the significant differences in vehicle driving patterns among the cities. Peng et al. [21] gathered 18 buses' routes; constructed a cyclic condition in Fuzhou, China; and developed a 1227 s speed-time series. These driving cycles exhibited more dynamic driving characteristics and could more accurately verify the energy consumption, power battery health status, and driving range of electric vehicles. The actual driving cycle is generally divided into four different kinematic fragments, including idling, accelerating, cruising, and decelerating [22]. The driving cycle velocity varies with the time changing and can only represent transient conditions in a special area and a specifically measured single sequence of defined length [23,24]. Therefore, the electric vehicle cannot follow the driving cycle and it cannot detect the safety hazards of the electric vehicle.
In response to the problem that transient driving cycles cannot be realistically followed, most of the official recognized electric vehicle range cycle testing in China uses the New European Driving Cycle (NEDC) [25]. Vehicle Mass Analysis System (VMAS) is commonly used for testing pollutant emissions, fuel consumption, and vehicle certification of fuelpowered passenger cars [26,27]. These characteristics of two driving cycles are because of a modal cycle consisting of a series of speed-time data, which are followed in the chassis dynamometer when testing the vehicle, usually for less than 30 min. Both the area-specific transient driving cycles and the official standard driving cycles have the following disadvantages.
(1) Long cycle times. For example, the total durations of the NEDC, World Light Vehicle Test Cycle (WLTC), and China Light-Duty Vehicle Test Cycle-Passenger (CLTC-P) are 1180 s, 1800 s, and 1800 s. In addition, the amount of data and the complex combinations are difficult to find in real life. (2) The driving cycle inspection items are single and cannot complete the practical operation test, especially the dynamic test of electric vehicle operational safety.
The test of other components of electric vehicles is mainly reflected in the diagnosis of power battery, such as high and fast battery temperatures, charging, and discharging [28]. Diagnosis of these components is mainly done by sensor detection of limited characteristic parameters for faults, however, simultaneous sensor faults can be diagnosed as battery faults [29,30].
In summary, much of the current research work on electric vehicles is focused on developing a standard driving cycle, testing the energy consumption of the vehicle, or detecting hidden dangers in key electric vehicle components. However, there is very little research into the operational safety testing of electric vehicles, and there is still a gap based on the driving cycle for the operational safety testing of electric vehicles.
The intention of this paper is to solve the above problems. A novel method of developing a driving cycle is proposed based on the rapid inspection of operational safety, which could detect whether the electric vehicle operational safety system is abnormal in a short (1) We innovatively propose three inspection items according to the operational safety of national standards. Meanwhile, an inspection calculation method of operational safety is developed based on the acceleration changing rate. (2) The multi-cycle inspection method with the stable pedal mode is developed and collects sufficient driving data by OBD. Gauss filtering algorithm is applied for data preprocessing. The support vector machine is adopted to construct the rapid inspection driving cycle. (3) Validation of the developed rapid driving cycle on the test bench. The thresholds for gliding safety, driving safety, and braking safety are evaluated based on the test results.
The methodology is illustrated in the form of a flow chart as shown in Figure 1. Compared to other driving cycles, this cycle focuses more on operational safety than on driving range and energy consumption, and could quickly detect operational safety hidden dangers for electric vehicles in a short time. (1) We innovatively propose three inspection items according to the operational safety of national standards. Meanwhile, an inspection calculation method of operational safety is developed based on the acceleration changing rate. (2) The multi-cycle inspection method with the stable pedal mode is developed and collects sufficient driving data by OBD. Gauss filtering algorithm is applied for data preprocessing. The support vector machine is adopted to construct the rapid inspection driving cycle. (3) Validation of the developed rapid driving cycle on the test bench. The thresholds for gliding safety, driving safety, and braking safety are evaluated based on the test results.
The methodology is illustrated in the form of a flow chart as shown in Figure 1. Compared to other driving cycles, this cycle focuses more on operational safety than on driving range and energy consumption, and could quickly detect operational safety hidden dangers for electric vehicles in a short time.

Electric Vehicle Inspection Items
In recent years, abnormal acceleration and braking failure of electric vehicles have frequently led to traffic accidents. The intervention of braking energy recovery strategy during the braking process has led to mutation in speed, which can easily cause accidents for the driver. Since 2005, China has distributed a series of standards for testing the operational safety performance of motor vehicles and electric vehicles. For example: GB21670-2008, GBT 12543-2009, GBT 18385-2005, GBT 35179-2017, and GB 38900-2020, which were strict testing standards for the acceleration performance and braking performance of the vehicles, as shown in Table 1.

Electric Vehicle Inspection Items
In recent years, abnormal acceleration and braking failure of electric vehicles have frequently led to traffic accidents. The intervention of braking energy recovery strategy during the braking process has led to mutation in speed, which can easily cause accidents for the driver. Since 2005, China has distributed a series of standards for testing the operational safety performance of motor vehicles and electric vehicles. For example: GB21670-2008, GBT 12543-2009, GBT 18385-2005, GBT 35179-2017, and GB 38900-2020, which were strict testing standards for the acceleration performance and braking performance of the vehicles, as shown in Table 1.
According to the testing requirements of national standards on the safety of electric vehicles in this paper, gliding safety, driving safety, and braking safety were determined as the inspection items, respectively, which were directly related to operational safety. The main characteristics of the three inspection items are as follows: Driving stability is an important indicator to evaluate driving safety during the driving of electric vehicles. If there is a large jitter in the acceleration of the vehicle within a short period, it will be very easy to cause a traffic accident, in addition to bringing a certain amount of discomfort to the driver and passengers. Therefore, driving stability is set as the inspection item in this paper. The indicator for gliding safety during the gliding is the gliding comfort. This is because if the vehicle jitters during the gliding process, it will lead to a slow or sharp down drop in velocity, which is very likely to cause traffic accidents. Therefore, gliding comfort is set as the inspection item. Braking coordination of electromechanical braking and regenerative braking is proposed as an important indicator to evaluate braking safety during the braking of an electric vehicle. When the coordination between electromechanical braking and regenerative braking is poor in the braking process, it will affect the braking performance. Therefore, the braking coordination is set as the inspection item.

Test Site and Route Selection
The determination of the test route during the development of a conventional driving cycle requires the designation of an urban road in a certain area. However, the rapid inspection driving cycle is on a straight driving road to test in this paper. We do not need to consider road types, traffic flow, time periods, economic differences, and central business district. Therefore, the test route is chosen on the straight runway of the Chang'an University Internet of Vehicles and Intelligent Vehicles Test Field, in which the straight runway length is 1.1 km and the slope of the test route is 0 • . So, there is no effect of vehicle slope resistance. The test site and route are shown in Figure 2.

Data Collection and Test Method
We selected a battery electric vehicle as the test vehicle, which had more than 22,000 in the market and was a representative small family car among the battery electric vehicle in this paper. The technical parameters of the test vehicle are given in Table 2.

Data Collection and Test Method
We selected a battery electric vehicle as the test vehicle, which had more than 22,000 in the market and was a representative small family car among the battery electric vehicle in this paper. The technical parameters of the test vehicle are given in Table 2. According to the NEDC, full section suburban cycle maximum acceleration was ≤1.5 m/s 2 , suburban cycle maximum deceleration was ≥−1.0 m/s 2 , urban cycle acceleration was ≥0.15 m/s 2 , and urban cycle deceleration was ≤−0.15 m/s 2 [25]. So, the acceleration range of the acceleration segment was confirmed in [0.15 m/s 2 , 1.5 m/s 2 ], the deceleration range of deceleration segment was [−1.5 m/s 2 , −0.15 m/s 2 ], and the maximum speed was 65 km/h. To precisely determine the range of acceleration, a limit device was used to lock the accelerator/brake pedal opening. The driver drove the test vehicle to 65 km/h with three modes of accelerator pedal opening of 30%, 60%, and 90% respectively, then braked to 0 with three modes of brake pedal opening of 30%, 60%, and 90% respectively. Finally, the acceleration and deceleration values under different pedal modes were calculated. The results are shown in Figure 3. three modes of brake pedal opening of 30%, 60%, and 90% respectively. Finally, the acceleration and deceleration values under different pedal modes were calculated. The results are shown in Figure 3.
The results showed that the maximum acceleration for 30% accelerator pedal opening was 1.45 m/s 2 , the maximum acceleration for 60% accelerator pedal opening was 2.66 m/s 2 , and the maximum acceleration for 90% accelerator pedal opening was 3.83 m/s 2 . The maximum deceleration for 30% brake pedal opening was 1.36 m/s 2 , the maximum deceleration for 60% brake pedal opening was 2.98 m/s 2 , and the maximum deceleration for 90% brake pedal opening was 4.45 m/s 2 . Therefore, 30% accelerator/brake pedal opening satisfies the design requirements in this paper. The accelerator and brake pedal mode were set at 30% pedal opening.
Based on the set pedal opening, the test method designed in this paper was divided into two parts. The first step was designed to accelerate the vehicle to 65 km/h with the 30% accelerator pedal opening and then released the accelerator pedal gliding to the minimum steady speed by the driver in vehicle SOC between 20 and 80%. The second step was in the same vehicle state, the driver accelerated the vehicle to 65 km/h with the 30% accelerator pedal opening, then released the accelerator pedal to take the vehicle gliding to 60 km/h, and braked the vehicle to 0 km/h with the 30% brake pedal opening. The test was repeated several times until it met the test requirements.  The methods of data collection include the chase car method, on-board measurement method, and a combination of both in the urban driving cycle [31]. The chase car method means the driver randomly follows the target vehicle on a predetermined route, and if the target vehicle drives out of the test area or behaves abnormally, the test vehicle will find another target to follow [32]. The on-board measurement method installs a global position system (GPS) and on-board diagnosis (OBD) on the test vehicle to record its driving information along a predetermined route to obtain reliable data with high cost [33,34]. To develop the rapid inspection driving cycle, this paper used data from real vehicle testing at specific speeds under real driving. So, it was an inspection driving cycle for operational safety, which did not require tracking measurements on urban vehicles. Therefore, on-  The results showed that the maximum acceleration for 30% accelerator pedal opening was 1.45 m/s 2 , the maximum acceleration for 60% accelerator pedal opening was 2.66 m/s 2 , and the maximum acceleration for 90% accelerator pedal opening was 3.83 m/s 2 . The maximum deceleration for 30% brake pedal opening was 1.36 m/s 2 , the maximum deceleration for 60% brake pedal opening was 2.98 m/s 2 , and the maximum deceleration for 90% brake pedal opening was 4.45 m/s 2 . Therefore, 30% accelerator/brake pedal opening satisfies the Sustainability 2022, 14, 5079 6 of 16 design requirements in this paper. The accelerator and brake pedal mode were set at 30% pedal opening.
Based on the set pedal opening, the test method designed in this paper was divided into two parts. The first step was designed to accelerate the vehicle to 65 km/h with the 30% accelerator pedal opening and then released the accelerator pedal gliding to the minimum steady speed by the driver in vehicle SOC between 20 and 80%. The second step was in the same vehicle state, the driver accelerated the vehicle to 65 km/h with the 30% accelerator pedal opening, then released the accelerator pedal to take the vehicle gliding to 60 km/h, and braked the vehicle to 0 km/h with the 30% brake pedal opening. The test was repeated several times until it met the test requirements.
The methods of data collection include the chase car method, on-board measurement method, and a combination of both in the urban driving cycle [31]. The chase car method means the driver randomly follows the target vehicle on a predetermined route, and if the target vehicle drives out of the test area or behaves abnormally, the test vehicle will find another target to follow [32]. The on-board measurement method installs a global position system (GPS) and on-board diagnosis (OBD) on the test vehicle to record its driving information along a predetermined route to obtain reliable data with high cost [33,34]. To develop the rapid inspection driving cycle, this paper used data from real vehicle testing at specific speeds under real driving. So, it was an inspection driving cycle for operational safety, which did not require tracking measurements on urban vehicles. Therefore, onboard measurements were used to collect real-time data on driving in this paper. The GPS was used to collect the speed and time of the test vehicle to modify the OBD data, which was stored in real time on the PC. At the same time, this was undertaken to increase the validity and effectiveness of the test data and avoid spike data in the acceleration signal. Considering the OBD-collected data within the effective frequency range, the sampling frequency was set at 20 Hz in this paper. The test equipment is shown in Figure 4. In addition, the test was operated by one professional driver throughout to eliminate data jitter caused by different drivers.

Data Processing
The main function of data processing is to confirm velocity time series as well as to remove noise and duplicate data from the original data. This process had two steps. The first step was to interpolate the original data with a polynomial interpolation algorithm for 10 ms. The sampling frequency of OBD to collect the data was sampled at 20 Hz and after a single data collection test, the test results showed that the time interval between the collected data was 10 ms. So, the original data was interpolated at 10 ms. Due to the limitations of the OBD collection data, the collected velocity and time samples will have stopped repetitive data and therefore need to be interpolated to obtain a new speed-time relationship. The second step was to smooth the original data. The original data collected by the sensor often show burrs or sudden changes in the data, due to the driving environmental disturbances and buildings, and some disturbances in the driving environment can also affect the test data. A Gaussian smoothing algorithm was introduced to filter and denoise the data. Gaussian smoothing is a linear smoothing filter, which has a softer smoothing effect and better edge retention than mean filtering. The original and denoising data for the driving-gliding state speed is shown in Figure 5. According to the data results, the minimum speed of the braking energy recovery is 8 km/h. The driving-braking state speed is shown in Figure 6.  The motion data was divided into three different kinematic segments, including driving segments, gliding segments, and braking segments after data processing. To increase the accuracy of testing the operational safety of the electric vehicle, the three different kinematic segments were combined with the idling and cruising kinematic segments in the construction of the rapid inspection driving cycle.
To effectively inspect the safety of the driving, gliding, and braking according to the determined operational safety testing items, the velocity mutation in speed during the sliding time window was purposed as an evaluation index. Acceleration can be used as an important feature to evaluate velocity mutation, but as acceleration only reflects the velocity variation during the whole process, it cannot characterize its operational stability. Therefore, the first derivative of the acceleration in the sliding time window for each kin-  The motion data was divided into three different kinematic segments, including driving segments, gliding segments, and braking segments after data processing. To increase the accuracy of testing the operational safety of the electric vehicle, the three different kinematic segments were combined with the idling and cruising kinematic segments in the construction of the rapid inspection driving cycle.
To effectively inspect the safety of the driving, gliding, and braking according to the determined operational safety testing items, the velocity mutation in speed during the sliding time window was purposed as an evaluation index. Acceleration can be used as an important feature to evaluate velocity mutation, but as acceleration only reflects the velocity variation during the whole process, it cannot characterize its operational stability. Therefore, the first derivative of the acceleration in the sliding time window for each kin- The motion data was divided into three different kinematic segments, including driving segments, gliding segments, and braking segments after data processing. To increase the accuracy of testing the operational safety of the electric vehicle, the three different kinematic segments were combined with the idling and cruising kinematic segments in the construction of the rapid inspection driving cycle.
To effectively inspect the safety of the driving, gliding, and braking according to the determined operational safety testing items, the velocity mutation in speed during the sliding time window was purposed as an evaluation index. Acceleration can be used as an important feature to evaluate velocity mutation, but as acceleration only reflects the velocity variation during the whole process, it cannot characterize its operational stability. Therefore, the first derivative of the acceleration in the sliding time window for each kinematic segment, i.e., the acceleration changing rate, is used as an indicator to evaluate its safety in this paper.
After testing the experimental data, t was used as the sliding time window, and the ∆t was determined as the velocity sampling interval according to the data sampling frequency. The acceleration calculation formula was established for all velocity segments, which can be written as: Based on the acceleration results obtained, we selected t as the sliding time window, and determined the ∆t as the acceleration sampling interval. Finally, the acceleration changing rate was established for all kinematic segments, as shown in the following equation: where: i = (1, 2, 3, · · · , n − ∆t).

Driving Cycle Construction
The construction of a rapid inspection driving cycle is crucial for testing the operational safety of electric vehicles, and the nonlinear classification of speed-time directly affects the test result. The support vector machine (SVM) algorithm is an effective method for solving nonlinear pattern recognition and is mainly used to solve classification and regression problems. The support vector machine was referenced to solve the problem of regression curves for the rapid driving cycle kinematic segments in this paper. The core idea of the regression problem requires sample data to construct bounded training sets, which can be represented as (X i ,Y i ) (i = 1, 2, 3, . . . , n). Where n is the number of samples, X is the input parameter, and Y is the output parameter. Regression aims to find the optimal function f (x), which can be written as follows: where , is the dot product; ω is the parameter vector; and b is the bias vector. The dot product of two numbered sets that satisfy the Mercer condition can be approximated by a kernel function to reduce the complexity of the calculation in reproducing the kernel Hilbert space.
For the centralization of experimental datasets, the Gaussian radial basis kernel function has good generalization ability and good local data fitting ability [35]. So, the Gaussian radial basis kernel function was chosen as the kernel function of the support vector machine, which is as follows: where σ is the dot product. We set ω as a parameter vector, and ξ i and ξ i * as slack variables. The optimization problems will recast as follows: where C is the penalty coefficient and C needs to be greater than 0; ε is an insensitive coefficient to control the width of the regression function for the region of insensitivity of the data sample [36].
To take the validity of the regression fitting data, we set up the root mean square error (RMSE) to indicate the curve dispersion, which is shown as follows: where yi is the regression value, and f (xi ) is the true value. A support vector machine algorithm was used to regress the original driving cycle curve. The results showed that the RMSE was less than 10%, therefore, the candidate driving cycles were acceptable.
To further establish a completely electric vehicle operation safety rapid inspection driving cycle and improve the accuracy of the inspection, three kinematic segments of driving, gliding, and braking segment were set. The kinematic segments are shown in Figure 7. The total duration of the driving kinematic segment was 11.3 s, the gliding kinematic segment was 14.8 s, and the braking kinematic segment was 12.2 s. The different kinematic segments were spliced together through the idling and cruising segments, so the driver could have sufficient preparation time to switch and follow each segment of the driving cycle. According to the time of NEDC, the time of idling and cruising segments was set at 10 s. According to the minimum braking energy recovery speed of 10 km/h for the inspection items, the speed of the cruising segments was set at 10 km/h in this paper.
of the data sample [36].
To take the validity of the regression fitting data, we set up the root mean square error (RMSE) to indicate the curve dispersion, which is shown as follows: A support vector machine algorithm was used to regress the original driving cycle curve. The results showed that the RMSE was less than 10%, therefore, the candidate driving cycles were acceptable.
To further establish a completely electric vehicle operation safety rapid inspection driving cycle and improve the accuracy of the inspection, three kinematic segments of driving, gliding, and braking segment were set. The kinematic segments are shown in Figure 7. The total duration of the driving kinematic segment was 11.3 s, the gliding kinematic segment was 14.8 s, and the braking kinematic segment was 12.2 s. The different kinematic segments were spliced together through the idling and cruising segments, so the driver could have sufficient preparation time to switch and follow each segment of the driving cycle. According to the time of NEDC, the time of idling and cruising segments was set at 10 s. According to the minimum braking energy recovery speed of 10 km/h for the inspection items, the speed of the cruising segments was set at 10 km/h in this paper.
The driving kinematic segments were divided into 0-65 km/h and 10-65 km/h, where the 10-65 km/h segment was the vehicle driving to 65 km/h and then for gliding inspection. The gliding kinematic segments were divided into 65-60 km/h and 60-10 km/h, where 65-60 km/h took the vehicle speed decline to the inspection speed and also provided preparation time for the driver to the next stage of the driving operation. The 10-65 km/h driving segment and 65-60 km/h gliding segment are shown in Figure 8. The idling was for the driver to adjust the driving state, so the waiting segment was set 5 s before the inspection started. The rapid inspection driving cycle of the electric vehicle was constructed by splice and reconstruction between the kinematic segments. The developed rapid inspection driving cycle is shown in Figure 9.  The driving kinematic segments were divided into 0-65 km/h and 10-65 km/h, where the 10-65 km/h segment was the vehicle driving to 65 km/h and then for gliding inspection. The gliding kinematic segments were divided into 65-60 km/h and 60-10 km/h, where 65-60 km/h took the vehicle speed decline to the inspection speed and also provided preparation time for the driver to the next stage of the driving operation. The 10-65 km/h driving segment and 65-60 km/h gliding segment are shown in Figure 8. The idling was for the driver to adjust the driving state, so the waiting segment was set 5 s before the inspection started. The rapid inspection driving cycle of the electric vehicle was constructed by splice and reconstruction between the kinematic segments. The developed rapid inspection driving cycle is shown in Figure 9.  As can be seen from Figure 9, the total duration of the rapid inspection driving cycle is 204 s. There were three 0-65 km/h driving segments for inspecting the safety of the driving stability of the electric vehicle, three 65-10 km/h gliding segments for the inspection of the safety of the gliding comfort of the electric vehicle, three 60-0 km/h braking segments for the inspection of the safety of the braking coordination of the electric vehicle, and three 65-60 km/h gliding segments and three idling segments. The gliding segment was only set for this test vehicle, which cannot restrict the gliding time for different vehicles. Compared to the NEDC, this model avoided the problems of a high number of parked vehicles, long duration, and speed instability. This driving cycle provides a solution for accurate and rapid inspection of the operational safety of electric vehicles.
Based on the analysis of the testing data, the sliding time window was 0.5 s. Since the sampling frequency was 20 Hz, the speed sampling interval t Δ was 50. The range of acceleration changing rate was [−0.35 −0.04] for the driving segment, the range of acceleration changing rate for the gliding segment was [0.05 0.09], and the range of acceleration changing rate for the braking segment was [−0.04 0.095]. Results are shown in Figure 10.  As can be seen from Figure 9, the total duration of the rapid inspection driving cycle is 204 s. There were three 0-65 km/h driving segments for inspecting the safety of the driving stability of the electric vehicle, three 65-10 km/h gliding segments for the inspection of the safety of the gliding comfort of the electric vehicle, three 60-0 km/h braking segments for the inspection of the safety of the braking coordination of the electric vehicle, and three 65-60 km/h gliding segments and three idling segments. The gliding segment was only set for this test vehicle, which cannot restrict the gliding time for different vehicles. Compared to the NEDC, this model avoided the problems of a high number of parked vehicles, long duration, and speed instability. This driving cycle provides a solution for accurate and rapid inspection of the operational safety of electric vehicles.
Based on the analysis of the testing data, the sliding time window was 0.5 s. Since the sampling frequency was 20 Hz, the speed sampling interval t Δ was 50. The range of acceleration changing rate was [−0.35 −0.04] for the driving segment, the range of acceleration changing rate for the gliding segment was [0.05 0.09], and the range of acceleration changing rate for the braking segment was [−0.04 0.095]. Results are shown in Figure 10. As can be seen from Figure 9, the total duration of the rapid inspection driving cycle is 204 s. There were three 0-65 km/h driving segments for inspecting the safety of the driving stability of the electric vehicle, three 65-10 km/h gliding segments for the inspection of the safety of the gliding comfort of the electric vehicle, three 60-0 km/h braking segments for the inspection of the safety of the braking coordination of the electric vehicle, and three 65-60 km/h gliding segments and three idling segments. The gliding segment was only set for this test vehicle, which cannot restrict the gliding time for different vehicles. Compared to the NEDC, this model avoided the problems of a high number of parked vehicles, long duration, and speed instability. This driving cycle provides a solution for accurate and rapid inspection of the operational safety of electric vehicles.
Based on the analysis of the testing data, the sliding time window was 0.5 s. Since the sampling frequency was 20 Hz, the speed sampling interval ∆t was 50. The range of acceleration changing rate was [−0.35, −0.04] for the driving segment, the range of acceleration changing rate for the gliding segment was [0.05, 0.09], and the range of acceleration changing rate for the braking segment was [−0.04, 0.095]. Results are shown in Figure 10. As shown in Figure 10, the driving segment has the largest range of accelera changing rate. As the vehicle accelerates from 0 to 65 km/h, the velocity differ becomes progressively smaller under the sliding time window , which causes the re to be a negative value. The trend of gradual flattening of the increasing velocity w stable pedal opening is verified, and the qualified range of driving stability for the dri segment is [−0.35 −0.04]. The second interval range of acceleration changing rate braking segment. When the vehicle enters the braking state from gliding, electromechanical brake starts to intervene to completely stop. The velocity differenc the sliding time window range becomes larger, so the acceleration changing fluctuates more, while there is no sudden change. The qualified range of bra coordination for the braking segment is [−0.04 0.095]. The gliding segment has the sma range of acceleration changing rate, which is mainly due to the intervention of regenerative braking during the gliding process. The external dominant factors of velo reduction are reduced, and the vehicle can find stability in gliding. The qualified rang the changing rate for gliding smoothness for the gliding segment is [0.05 0.09].

Driving Cycle Verification
To verify the representativeness and typicality of the developed driving cycle, d ing stability, gliding comfort, and braking coordination were detected for inspec based on the established operational safety test calculation method. The test vehicle used to follow the rapid driving cycle on a testing bench in Chang'an university Auto tive Performance Test Laboratory, which is shown in Figure 11. The vehicle inertia loading force were matched before following the rapid driving cycle. Therefore, the be could simulate the road test. The maximum permissible error range for the rapid dri cycle velocity and following velocity is [−2 km/h 2 km/h]. If the following velocity exc the maximum permissible error range, the following test data is considered ineffic Simultaneously, the driver selected the same professional person to operate the ele vehicle and obtained a large number of efficiency data after following the test. As shown in Figure 10, the driving segment has the largest range of acceleration changing rate. As the vehicle accelerates from 0 to 65 km/h, the velocity difference becomes progressively smaller under the sliding time window, which causes the result to be a negative value. The trend of gradual flattening of the increasing velocity with a stable pedal opening is verified, and the qualified range of driving stability for the driving segment is [−0.35, −0.04]. The second interval range of acceleration changing rate is a braking segment. When the vehicle enters the braking state from gliding, the electromechanical brake starts to intervene to completely stop. The velocity difference in the sliding time window range becomes larger, so the acceleration changing rate fluctuates more, while there is no sudden change. The qualified range of braking coordination for the braking segment is [−0.04, 0.095]. The gliding segment has the smallest range of acceleration changing rate, which is mainly due to the intervention of the regenerative braking during the gliding process. The external dominant factors of velocity reduction are reduced, and the vehicle can find stability in gliding. The qualified range of the changing rate for gliding smoothness for the gliding segment is [0.05, 0.09].

Driving Cycle Verification
To verify the representativeness and typicality of the developed driving cycle, driving stability, gliding comfort, and braking coordination were detected for inspection based on the established operational safety test calculation method. The test vehicle was used to follow the rapid driving cycle on a testing bench in Chang'an university Automotive Performance Test Laboratory, which is shown in Figure 11. The vehicle inertia and loading force were matched before following the rapid driving cycle. Therefore, the bench could simulate the road test. The maximum permissible error range for the rapid driving cycle velocity and following velocity is [−2 km/h, 2 km/h]. If the following velocity exceeds the maximum permissible error range, the following test data is considered inefficient. Simultaneously, the driver selected the same professional person to operate the electric vehicle and obtained a large number of efficiency data after following the test. Based on the analysis of the following test data, six driving segments were sele for calculation and obtained the range of the acceleration changing rate. The re showed the maximum driving acceleration changing rate was −0.49, and the minim driving acceleration changing rate was −0.03. The maximum RMSE in the accelera changing rate between the constructed rapid inspection driving segments and the te driving segments was 9%. The maximum RMSE in the acceleration changing rate betw the tested driving segments was 6%. Then, the range of acceleration changing rate o six gliding segments was calculated. The result showed the maximum gliding accelera changing rate was 0.13, and the minimum gliding acceleration changing rate was The maximum RMSE in the acceleration changing rate between the constructed rapi spection gliding segments and the tested gliding segments was 2.7%. The maxim RMSE in the acceleration changing rate between the tested gliding segments was 4% errors results were shown in Table 3. Finally, the range of the acceleration changing of six braking segments was calculated. The result showed a maximum braking acce tion changing rate of 0.1 and a minimum braking acceleration changing rate of 0.02. maximum RMSE in the acceleration changing rate between the constructed rapid ins tion braking segments and the tested braking segments was 6.9%. The maximum R in the acceleration changing rate between the tested braking segments was 3.4%. The results were valid and errors are in the calculated permissible range, which is show Figure 12.  Based on the analysis of the following test data, six driving segments were selected for calculation and obtained the range of the acceleration changing rate. The results showed the maximum driving acceleration changing rate was −0.49, and the minimum driving acceleration changing rate was −0.03. The maximum RMSE in the acceleration changing rate between the constructed rapid inspection driving segments and the tested driving segments was 9%. The maximum RMSE in the acceleration changing rate between the tested driving segments was 6%. Then, the range of acceleration changing rate of the six gliding segments was calculated. The result showed the maximum gliding acceleration changing rate was 0.13, and the minimum gliding acceleration changing rate was 0.03. The maximum RMSE in the acceleration changing rate between the constructed rapid inspection gliding segments and the tested gliding segments was 2.7%. The maximum RMSE in the acceleration changing rate between the tested gliding segments was 4%. The errors results were shown in Table 3. Finally, the range of the acceleration changing rate of six braking segments was calculated. The result showed a maximum braking acceleration changing rate of 0.1 and a minimum braking acceleration changing rate of 0.02. The maximum RMSE in the acceleration changing rate between the constructed rapid inspection braking segments and the tested braking segments was 6.9%. The maximum RMSE in the acceleration changing rate between the tested braking segments was 3.4%. The test results were valid and errors are in the calculated permissible range, which is shown in Figure 12. on the sudden change of speed and the operability of the driver, the rapid condition designed in this paper has a low speed and a small pedal opening, which is conducive to the following test. On the other hand, we chose a professional driver to operate the vehicle, and after several simulation tests, it followed the rapid driving cycle well, and the speed error reached the minimum. The testing results showed that the proposed rapid inspection method was very effective in verifying the operational safety of the electric vehicle. The changing rate of each acceleration was within a small range of variation, which could characterize the safety of the driving, gliding, and braking movement. The maximum changing rate of acceleration during driving at less than 0.5 is the safety threshold for driving stability. The maximum changing rate of acceleration during gliding is less than 0.1 as the safety threshold for gliding comfort. The maximum changing rate of acceleration during braking is less than 0.1 as the safety threshold for vehicle braking coordination.

Conclusions
The operational safety of the electric vehicle is the core in full vehicle safety, and accurate inspection of its safety hazard is a premise for ensuring the operational safety of electric vehicles. However, there is no testing of relevant items according to specific follow-on testing cycles in traditional electric vehicle operation safety inspection. To resolve this issue, an in-depth analysis of operational safety inspection technology at home and abroad, a rapid inspection method was proposed for electric vehicle operational safety. According to the inspection results, the conclusion is as follows: In order to better reflect the results of the following test, this paper is built from both the construction of the rapid driving cycle and the following test. On the one hand, based on the sudden change of speed and the operability of the driver, the rapid condition designed in this paper has a low speed and a small pedal opening, which is conducive to the following test. On the other hand, we chose a professional driver to operate the vehicle, and after several simulation tests, it followed the rapid driving cycle well, and the speed error reached the minimum.
The testing results showed that the proposed rapid inspection method was very effective in verifying the operational safety of the electric vehicle. The changing rate of each acceleration was within a small range of variation, which could characterize the safety of the driving, gliding, and braking movement. The maximum changing rate of acceleration during driving at less than 0.5 is the safety threshold for driving stability. The maximum changing rate of acceleration during gliding is less than 0.1 as the safety threshold for gliding comfort. The maximum changing rate of acceleration during braking is less than 0.1 as the safety threshold for vehicle braking coordination.

Conclusions
The operational safety of the electric vehicle is the core in full vehicle safety, and accurate inspection of its safety hazard is a premise for ensuring the operational safety of electric vehicles. However, there is no testing of relevant items according to specific follow-on testing cycles in traditional electric vehicle operation safety inspection. To resolve this issue, an in-depth analysis of operational safety inspection technology at home and abroad, a rapid inspection method was proposed for electric vehicle operational safety. According to the inspection results, the conclusion is as follows: (1) According to the existing electric vehicles inspection methods and means, combined with the relevant standards for electric vehicles and the accident causes, the electric vehicle safety inspection items were identified, which included gliding comfort, driving stability, and braking coordination. Then, an acceleration changing rate was purposed as an items inspection method based on the identified inspection items. Using the basic theory of the driving cycle, the accelerator/brake pedal opening, acceleration, and velocity were selected as test parameters. Then, a multi-cycle operational safety test method was proposed with stable pedal mode and set several accelerator/braking cycle tests were set to obtain the operational safety test data. (2) Advanced filtering algorithms were used to remove burrs and duplicate data from the collected data. A support vector machine was used for regression cycles and fused with a spliced rapid inspection driving cycle with a total time of 204 s for electric vehicles. Then the range of three kinematic segments stability was determined. Finally, this rapid inspection driving cycle was verified on the testing bench. The verification results showed that the maximum driving acceleration changing rate error was 9%. The maximum tested driving acceleration changing rate error was 6%, and the maximum verification gliding acceleration changing rate error was 2.7%. The maximum tested gliding acceleration changing rate error was 4%, and the maximum verification braking acceleration changing rate error was 6.9%. The maximum tested braking acceleration changing rate error was 3.4%. All three inspection items had a good calculation effect. The safety thresholds for driving acceleration changing rate were less than 0.5, the gliding acceleration changing rate was less than 0.1, and the braking acceleration changing rate was less than 0.1. (3) A rapid inspection driving cycle was established for the operational safety of electric vehicles, which provided some technical support for the "annual inspection" of electric vehicles. However, there are also certain limitations. For example, there are many models of electric vehicles with widely varying braking strategies. A single driver and changing different drivers may increase the calculation error. Therefore, further analysis can be performed for different braking energy recovery strategies, different vehicles and multi-cycle following correction. Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.

Conflicts of Interest:
The authors declare no conflict of interest.