CN113215881B - Method for determining minimum safe radius of flat curve of double-lane highway in plateau area - Google Patents

Method for determining minimum safe radius of flat curve of double-lane highway in plateau area Download PDF

Info

Publication number
CN113215881B
CN113215881B CN202110390149.5A CN202110390149A CN113215881B CN 113215881 B CN113215881 B CN 113215881B CN 202110390149 A CN202110390149 A CN 202110390149A CN 113215881 B CN113215881 B CN 113215881B
Authority
CN
China
Prior art keywords
electroencephalogram
heart rate
flat curve
driver
curve
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202110390149.5A
Other languages
Chinese (zh)
Other versions
CN113215881A (en
Inventor
陈飞
余瀚林
胡飞
薄雾
张平
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southeast University
Original Assignee
Southeast University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Southeast University filed Critical Southeast University
Priority to CN202110390149.5A priority Critical patent/CN113215881B/en
Publication of CN113215881A publication Critical patent/CN113215881A/en
Application granted granted Critical
Publication of CN113215881B publication Critical patent/CN113215881B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • E—FIXED CONSTRUCTIONS
    • E01—CONSTRUCTION OF ROADS, RAILWAYS, OR BRIDGES
    • E01C—CONSTRUCTION OF, OR SURFACES FOR, ROADS, SPORTS GROUNDS, OR THE LIKE; MACHINES OR AUXILIARY TOOLS FOR CONSTRUCTION OR REPAIR
    • E01C1/00—Design or layout of roads, e.g. for noise abatement, for gas absorption
    • E—FIXED CONSTRUCTIONS
    • E01—CONSTRUCTION OF ROADS, RAILWAYS, OR BRIDGES
    • E01C—CONSTRUCTION OF, OR SURFACES FOR, ROADS, SPORTS GROUNDS, OR THE LIKE; MACHINES OR AUXILIARY TOOLS FOR CONSTRUCTION OR REPAIR
    • E01C1/00—Design or layout of roads, e.g. for noise abatement, for gas absorption
    • E01C1/002—Design or lay-out of roads, e.g. street systems, cross-sections ; Design for noise abatement, e.g. sunken road
    • E—FIXED CONSTRUCTIONS
    • E01—CONSTRUCTION OF ROADS, RAILWAYS, OR BRIDGES
    • E01C—CONSTRUCTION OF, OR SURFACES FOR, ROADS, SPORTS GROUNDS, OR THE LIKE; MACHINES OR AUXILIARY TOOLS FOR CONSTRUCTION OR REPAIR
    • E01C1/00—Design or layout of roads, e.g. for noise abatement, for gas absorption
    • E01C1/007—Design or auxiliary structures for compelling drivers to slow down or to proceed with caution, e.g. tortuous carriageway; Arrangements for discouraging high-speed or non-resident traffic

Landscapes

  • Engineering & Computer Science (AREA)
  • Architecture (AREA)
  • Civil Engineering (AREA)
  • Structural Engineering (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
  • Control Of Driving Devices And Active Controlling Of Vehicle (AREA)

Abstract

The invention discloses a method for determining the minimum safe radius of a flat curve suitable for a double-lane road in a plateau area, which comprises the following steps: (1) acquiring vehicle speed, heart rate and electroencephalogram data of a driver driving on a flat curve road section with different altitudes and radii according to an on-site driving experiment; (2) carrying out sample entropy processing on electroencephalogram and heart rate indexes of a driver, and extracting characteristic indexes of the electroencephalogram and heart rate indexes; (3) performing principal component analysis on the heart rate sample entropy and the electroencephalogram signal sample entropy to obtain a comprehensive evaluation index CISE of the psychological load of the driver; (4) constructing a regression model of a comprehensive evaluation index CISE and the transverse force acceleration of the flat curve section in different altitude intervals; (5) and determining a threshold value of the comprehensive evaluation index through a Lagrange median theorem, and finally providing a minimum safe radius threshold value of the double-lane road in different altitude intervals of the plateau area. The invention improves the driving safety of highways in plateau areas.

Description

Method for determining minimum safe radius of flat curve of double-lane highway in plateau area
Technical Field
The invention relates to the technical field of road safety, in particular to a method for determining the minimum safe radius of a flat curve of a double-lane highway in a plateau area.
Background
At present, highlands are mostly low-grade double-lane highways. Due to the influence of complex terrain conditions and engineering economic conditions, the sections of the dual-lane highways in plateau areas, such as mountains, valleys and rivals, adjacent to cliffs are more, so that the low-linear index sections of the dual-lane highways in plateau areas are high in frequency. However, the road sections with lower linear indexes are usually the road sections with easy traffic accidents, and statistics shows that the major traffic accidents in the plateau of China basically occur on the two-lane roads. Because the low-pressure and low-oxygen environment in the plateau has certain negative influence on the psychophysiological condition of a driver, the linear design safety of the low-index flat curve section of the plateau double-lane road faces higher challenges.
According to the current research situation at home and abroad, most of the previous researches focus on the influence of vehicle dynamics on the minimum radius of a horizontal curve in a plain area. The correction of the minimum radius is not sufficient from the point of view of the specific physical and psychological condition of the driver on the plateau. Due to the unique climate and geographical environment of plateaus, the driver's physiological and psychological characteristics are different from those of other regions. The adverse effect of the low-pressure and low-oxygen special environment in the plateau area on the driving safety of drivers. Therefore, the research on the minimum radius of the horizontal curve provided by the plateau area has important significance for improving the safety of the plateau road.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and provide a method for determining the minimum safe radius of a flat curve of a two-lane highway in a plateau area.
The invention adopts the following technical scheme for solving the technical problems:
the invention provides a method for determining the minimum safe radius of a flat curve of a double-lane highway in a plateau area, which comprises the following steps of:
step 1, performing a field driving experiment in a plateau area, collecting position information, a vehicle speed, a heart rate and an electroencephalogram signal of a driver driving on a flat curve section with different altitudes and flat curve radiuses, and performing smooth noise reduction processing on the vehicle speed, the heart rate and the electroencephalogram signal;
step 2, carrying out sample entropy processing on the electroencephalogram signal and the heart rate of the driver after the smoothing and noise reduction processing in the step 1 to obtain characteristic indexes of electroencephalogram and heart rate indexes, wherein the characteristic indexes of the electroencephalogram and heart rate indexes refer to electroencephalogram sample entropy and heart rate sample entropy;
step 3, carrying out principal component analysis on the electroencephalogram sample entropy and the heart rate sample entropy to obtain a comprehensive evaluation index CISE of the physiological load of the heart of the driver;
step 4, according to the position information acquired in the step 1 and the vehicle speed after smooth noise reduction processing, constructing a comprehensive evaluation index CISE of the psychological and physiological load of the driver and the lateral force acceleration a borne by the driver in the flat curve road section in different altitude intervalshThe regression model of (2);
step 5, determining the CISE and the lateral force acceleration a in different altitude intervalshThe curve radius threshold value R corresponding to the curve mutation point in the regression modelsafe，RsafeAs the minimum safe radius which meets the condition that the psychological and physiological loads of the driver are in a safe state in different altitude intervals.
As a further optimization scheme of the method for determining the minimum safe radius of the flat curve of the double-lane highway in the plateau area, the step 1 is as follows:
the division principle of the flat curve section is that the radius R of a circular curve is less than or equal to 600m, the longitudinal slope I is less than or equal to 3%, a double-lane road meeting the requirements in different altitude intervals of a plateau area is selected, an outdoor real-vehicle driving experiment is carried out, a heart physiological feedback instrument is used for collecting the heart rate and electroencephalogram signals of a driver in the experiment process, and a GPS signal machine is used for recording position information and speed data in the driving process; and carrying out smooth noise reduction on the heart rate, the electroencephalogram signals and the speed data by utilizing Matlab software, and removing outliers.
The invention relates to a method for determining the minimum safe radius of a flat curve of a double-lane highway in a plateau area, which comprises the following steps: in the step 2, the electroencephalogram signals refer to beta/theta and beta/(theta + alpha), beta is beta wave in the electroencephalogram signals, theta is theta wave in the electroencephalogram signals, alpha is alpha wave in the electroencephalogram signals, and the heart rate refers to the beating times HR of the human heart per minute; the sample entropy processing is represented by sample entropy SampEn (m, r, N), wherein m is an embedding dimension, r is a similar tolerance, SD is a standard deviation of original data, and N is a data length, so that three electroencephalogram sample entropies SaEn (beta), SaEn (beta/theta), SaEn (beta/(theta + alpha)) and heart rate sample entropies SaEn (HR) of the electroencephalogram signal are obtained respectively.
The invention relates to a method for determining the minimum safe radius of a flat curve of a double-lane highway in a plateau area, which comprises the following steps: in the step 3, the principal component analysis method is to perform principal component analysis on the electroencephalogram sample entropy and the heart rate sample entropy obtained in the step 2, select the first n principal components with the cumulative contribution rate reaching the calculation threshold to determine the weight coefficient of each sample entropy in the integrated evaluation index CISE of the cardiac physiological load, and perform weighted average normalization on the coefficients by taking the variance contribution rate of the principal components as the weight to obtain the integrated evaluation index CISE of the cardiac physiological load;
Figure GDA0003103769150000031
in the formula: a. b, c and d are all weighted coefficients after weighted average normalization.
The invention relates to a method for determining the minimum safe radius of a flat curve of a double-lane highway in a plateau area, which comprises the following steps: in step 4, the different altitude intervals refer to the altitude range (3000 m-5000 m) divided into four altitude intervals (3000m, 3500m), (3500m, 4000m), (4000m, 4500m), (4500m, 5000 m); the lateral force acceleration ahThe expression is as follows:
Figure GDA0003103769150000032
wherein V is the running speed of the smooth curve road section after the smooth noise reduction treatment, R is the radius of the smooth curve, ihWhen the altitude is 10m, g is 9.807m/s2(ii) a H is 3000-3500 m, g is 9.796m/s2(ii) a H is 3500-4000 m, g is 9.794m/s2(ii) a H is 4000-4500 m, g is 9.793m/s2(ii) a H is 4500-5000 m, g is 9.791m/s2；
Comprehensive evaluation index CISE and lateral force acceleration a borne by driverhThe regression model of (a) is abouthThe regression form of the monotonically increasing negative power function curve of (1) is as follows:
CISE＝eah f
in the formula, e and f are regression coefficients, and both e and f are greater than 0.
As a further optimization scheme of the method for determining the minimum safe radius of the flat curve of the two-lane highway in the plateau area, step 5 is to use the flat curve radius threshold value RsafeThe determination steps are as follows:
respectively calculating CISE and a in each altitude interval by utilizing Lagrange median theoremhThe curvature mutation point of the regression model of (1), and a corresponding to the curvature mutation pointhThe value is taken as the safety state of meeting the physiological load of the driver in different altitude intervalsMaximum lateral force acceleration of state ahmaxThe maximum lateral force acceleration a is calculated using the following formulahmaxCorresponding flat curve radius threshold value Rsafe：
Figure GDA0003103769150000033
Wherein V is the running speed of the smooth curve road section after the smooth noise reduction treatment, g is the gravity acceleration, ihThe road section is a flat curve section ultrahigh value.
As a further optimization scheme of the method for determining the minimum safe radius of the flat curve of the double-lane highway in the plateau area, m is 2, r is 0.2SD, and N is 256.
Compared with the prior art, the invention adopting the technical scheme has the following technical effects:
the invention solves the defect that the linear index of the curved slope combined section in the existing road geometric linear design does not consider the adverse influence of the low-pressure oxygen-poor environment of the plateau area on the characteristics of the driver on the basis of the prior art, and simultaneously, the method disclosed by the invention is simple and convenient, provides a linear index selection method for the curved slope combined section of the two-lane road, and provides reference for the safety design of the road in the plateau area.
Drawings
FIG. 1 is a flow chart of the present invention.
FIG. 2 shows the acceleration a of CISE along with the lateral force in different altitude intervalshThe trend of change of (c).
Detailed Description
The following examples further illustrate the practice of the present invention, but the present invention is not limited to the following examples.
Example 1
A typical double-lane road in a plateau region, namely a national road 318 Linzhi city and a color season mountain-drawing section, is selected to carry out an on-site driving experiment, the experimental road section is a two-way two-lane road, the lane width is 3.5m, the design speed is 40km/h, the total length of the route is about 50km, the starting point height range of the selected road section is 2979m, and the end point height is 4720 m. The experimenter selects 10 young drivers who enter a plateau area for the first time (the plateau habit period is less than or equal to 10d) as experimental samples to carry out experiments, wherein 7 males and 3 females have the average age of 25.4 years and the average driving age of 4.7 years; the double-lane road means that the lane width is 3.5m, and the whole road comprises more flat curve sections. As shown in fig. 1, the method for calculating the minimum safe radius of the flat curve section of the dual-lane highway in the plateau area mainly comprises the following steps:
step (1): collecting vehicle speed, heart rate and electroencephalogram data of a driver in a plateau area at different altitudes and different radii of a flat curve section;
the linear structure of the experimental road is mainly a flat curve section, and the linear design data of the experimental road is shown in table 1.
TABLE 1 Linear combination distribution of combined sections of curved slopes
Number of curve combined segments Flat curve radius variation range (m) Vertical curve longitudinal slope variation range (%)
49 60-600 3-6
The heart rate and electroencephalogram data of experimenters are collected by an eight-channel multi-parameter biofeedback instrument produced by Thought Technology company of Canada, the collection frequency is 256Hz, and position information and speed data in the driving process are recorded by a GPS signal machine.
Step (2): carrying out sample entropy processing on the electroencephalogram signal and the heart rate of the driver, and extracting characteristic indexes of the electroencephalogram and heart rate indexes, namely electroencephalogram and heart rate sample entropy;
the electroencephalogram signal beta wave after the filtering and noise reduction processing of the tested person, the ratio values beta/theta, beta/(theta + alpha) and the heart rate HR are processed by adopting a sample entropy algorithm, wherein the embedding dimension m is 2, the similarity tolerance r is 0.2SD, the data length N is 256 (sample data of 1 s), and the heart rate sample entropy SaEn (HR) and the electroencephalogram sample entropy SaEn (beta), SaEn (beta/theta) and SaEn (beta/(theta + alpha)) are respectively obtained, and the statistics description is shown in Table 2.
Table 2 describes statistics
Mean value Standard deviation of Analysis of N
SaEn(β) 1.132404 .202231 1030
SaEn(β/θ) 1.13638 .233624 1030
SaEn(β/(θ+α)) 1.20570 .207640 1030
SaEn(HR) .34307 .113539 1030
And (3): performing principal component analysis on the heart rate sample entropy and the electroencephalogram signal sample entropy to obtain a comprehensive evaluation index CISE of the physiological load of the heart of the driver;
and (3) carrying out principal component analysis on the heart rate and electroencephalogram sample entropy indexes obtained in the step (2), wherein the specific process is as follows:
Figure GDA0003103769150000051
in the formula: x is the number ofi(i ═ 1,2 …, p) represents the original variables for which principal component analysis is required; y isj(j ═ 1,2 …, k) is the jth principal component of the original variable; a isijReferred to as the factor load, can be found by:
Figure GDA0003103769150000052
in the formula: lambda [ alpha ]iThe ith eigenvalue of a correlation coefficient matrix of the original variable; u. ofijIs the jth component of the ith feature vector.
The number of principal components is determined by the cumulative contribution ratio:
Figure GDA0003103769150000053
in the formula: ckRepresents the cumulative contribution rate of the first k principal components, CkMore than 85% should be achieved.
The obtained analysis results of the main components are shown in tables 3 to 5.
Total variance as explained in Table 3
Figure GDA0003103769150000054
Figure GDA0003103769150000061
TABLE 4 component matrix
Figure GDA0003103769150000062
TABLE 5 component score coefficient matrix
Figure GDA0003103769150000063
From the principal component analysis results, the characteristic root of the first principal component was 1.998, which accounts for 49.948% of the total variation, and the characteristic root of the second principal component was 1.405, which accounts for 35.218% of the total variation. The characteristic root of the first two principal components is larger than 1, the cumulative contribution rate reaches 85.166%, therefore, the first principal component and the second principal component are selected to determine the weighting coefficient of each sample entropy in the comprehensive index, and the linear combination of the two principal components is as follows:
Z1＝0.355*SaEn(β)+0.426*SaEn(β/θ)+0.438*SaEn(β/(θ+α))+0.037*SaEn(HR)
Z2＝-0.168*SaEn(β)+0.069*SaEn(β/θ)-0.144*SaEn(β/(θ+α))+0.979*SaEn(HR)
taking the variance contribution rate of the principal component as the weight, and normalizing the weighted average of the coefficients of the index in each principal component linear combination to obtain an expression of the psychophysiological load comprehensive entropy index CISE as follows; CISE may also be referred to as "cardiophysiological load integrated Sample Entropy" which is collectively referred to as "complete Indicator of Sample entry":
CISE＝0.127*SaEn(β)+0.254*SaEn(β/θ)+0.229*SaEn(β/(θ+α))+0.390*SaEn(HR)
and (4): constructing a regression model of the comprehensive entropy index CISE of the cardiac physiological load and the linear combination value CA of the combined section of the bent slope in different altitude intervals;
dividing a plateau altitude range (3000 m-5000 m) into four altitude intervals of (3000m, 3500m), (3500m, 4000m), (4000m, 4500m) and (4500m, 5000 m); transverse force acceleration a of a flat curve sectionhThe expression is as follows:
Figure GDA0003103769150000071
in the formula:
v is the running speed of the flat curve road section after smooth noise reduction treatment, and the unit is km/h;
r-radius of the flat curve, in m;
ih-the flat curve section is of ultra high value, in units%;
g-acceleration of gravity, unit m/s2When H is 10m, it is 9.807m/s2(ii) a H is 3000-3500 m, and 9.796m/s is taken2(ii) a H is 3500-4000 m, and 9.794m/s is taken2(ii) a H is 4000-4500 m, and 9.793m/s is taken2(ii) a H is 4500-5000 m, and 9.791m/s is taken2。
Analysis of comprehensive sample entropy CISE and transverse force acceleration a of cardiophysiological load of driver on flat curve road sectionhThe relation of (a) is made as the acceleration a of the curve section CISE along with the transverse force in different altitude intervalshThe scatter diagram of the variation trend of (2) is shown in fig. 2.
As can be seen from FIG. 2, the regression model of the integrated entropy index CISE of the cardiac physiological load is about the lateral force acceleration ahThe regression form of the monotonically increasing power function curve of (1) is as follows:
CISE＝eCAf
in the formula:
e. f-regression coefficient (> 0).
Comprehensive sample entropy CISE and transverse force acceleration a in different altitude intervals by means of SPSShThe variation trend of the sample is regressed to obtain the comprehensive sample entropy CISE and the transverse force acceleration a in different altitude intervalshThe regression model of (2):
(1)H＝3000～3500m
CISE＝1.004ah 0.634
the significance of the model is checked, and the model judges the coefficient R20.934, significance level Sig 0.000<0.001。
(2)H＝3500～4000m
CISE＝1.177ah 0.609
The significance of the model is checked, and the model judges the coefficient R20.895, significance level Sig 0.000<0.001。
(3)H＝4000～4500m
CISE＝1.304ah 0.570
The significance of the model is checked, and the model judges the coefficient R20.944, significance level Sig 0.000<0.001。
(4)H＝4500～5000m
CISE＝1.354ah 0.516
The significance of the model is checked, and the model judges the coefficient R20.915, significance level Sig 0.000<0.001。
And (5): determining flat curve radius threshold value R corresponding to mutation points of CISE regression curve in different altitude intervalssafeThe linear design index is used as the linear design index of the flat curve road section;
flat curve radius threshold RsafeRespectively calculating mutation points in each altitude interval by using the Lagrange median theorem to serve as flat curve radius value critical points meeting the physiological load safety of a driver, and the maximum transverse force acceleration a corresponding to each mutation pointhmaxThe calculation results are shown in table 6.
TABLE 6 CISE-based CA Critical values for different altitude intervals
Altitude interval (m) Ahmax(m/s2)
3000～3500 1.03
3500～4000 1.01
4000～4500 0.98
4500～5000 0.93
In different altitude intervals, the comprehensive sample entropy CISE shows a trend of increasing rapidly and then slowly along with the increase of the transverse acceleration AhIs less than AhmaxThe driver's sensitivity to lateral forces is higher, with ahThe operation complexity of a driver on a flat curve section is gradually increased, and the increase amplitude of the physiological load index CISE of the driver is larger; when a ishGreater than ahmaxMeanwhile, the psychophysiological load of the driver is in a relatively high state, the sensitivity of the driver to the transverse force is reduced, and the increase amplitude of the psychophysiological load is small. a ishmaThe corresponding CISE is a catastrophe point of curve change and corresponds to a demarcation point of the cardiophysiological load comprehensive sample entropy CISE change of the flat curve road section of the driver in different altitude intervals, so that a catastrophe point a of the regression curve can be changedhmaxAs a critical point of the lateral force acceleration meeting the psychological and physiological load safety of the driver, ahmaxSubstituting the following formula:
Figure GDA0003103769150000081
in the formula:
v is the running speed of the flat curve road section after smooth noise reduction treatment, and the unit is km/h;
g-acceleration of gravity, unit m/s2When H is 10m, it is 9.807m/s2(ii) a H is 3000-3500 m, and 9.796m/s is taken2(ii) a H is 3500-4000 m, and 9.794m/s is taken2(ii) a H is 4000-4500 m, and 9.793m/s is taken2(ii) a H is 4500-5000 m, and 9.791m/s is taken2。
ihThe flat curve section is ultrahigh in value, in units%.
A is tohmaxSubstituting the calculated flat curve radius value into a flat curve radius threshold value R meeting the psychological and physiological load safety of the driversafeWhen the linear design of the flat curve section is carried out, the radius value R of the flat curve is preferably controlled to be larger than RsafeSo as to ensure that the psychological and physiological load of a driver is in a safe state when the driver drives on the plateau two-lane highway flat curve road section.
The foregoing detailed description of the preferred embodiments of the invention has been presented. It should be understood that numerous modifications and variations could be devised by those skilled in the art in light of the present teachings without departing from the inventive concepts. Therefore, the technical solutions available to those skilled in the art through logic analysis, reasoning and limited experiments based on the prior art according to the concept of the present invention should be within the scope of protection defined by the claims.

Claims (7)

1. A method for determining the minimum safe radius of a flat curve of a double-lane highway in a plateau area is characterized by comprising the following steps:
step 1, performing a field driving experiment in a plateau area, collecting position information, a vehicle speed, a heart rate and an electroencephalogram signal of a driver driving on a flat curve section with different altitudes and flat curve radiuses, and performing smooth noise reduction processing on the vehicle speed, the heart rate and the electroencephalogram signal;
step 2, carrying out sample entropy processing on the electroencephalogram signal and the heart rate of the driver after the smoothing and noise reduction processing in the step 1 to obtain characteristic indexes of electroencephalogram and heart rate indexes, wherein the characteristic indexes of the electroencephalogram and heart rate indexes refer to electroencephalogram sample entropy and heart rate sample entropy;
step 3, carrying out principal component analysis on the electroencephalogram sample entropy and the heart rate sample entropy to obtain a comprehensive evaluation index CISE of the physiological load of the heart of the driver;
step 4, according to the position information acquired in the step 1 and the vehicle speed after smooth noise reduction processing, constructing a comprehensive evaluation index CISE of the psychological and physiological load of the driver and the lateral force acceleration a borne by the driver in the flat curve road section in different altitude intervalshThe regression model of (2);
step 5, determining the CISE and the lateral force acceleration a in different altitude intervalshThe curve radius threshold value R corresponding to the curve mutation point in the regression modelsafe，RsafeAs the minimum safe radius which meets the condition that the psychological and physiological loads of the driver are in a safe state in different altitude intervals.
2. The method for determining the minimum safe radius of the flat curve of the two-lane highway in the plateau area according to claim 1, wherein the step 1 is as follows:
the division principle of the flat curve section is that the radius R of a circular curve is less than or equal to 600m, the longitudinal slope I is less than or equal to 3%, a double-lane road meeting the requirements in different altitude intervals of a plateau area is selected, an outdoor real-vehicle driving experiment is carried out, a heart physiological feedback instrument is used for collecting the heart rate and electroencephalogram signals of a driver in the experiment process, and a GPS signal machine is used for recording position information and speed data in the driving process; and carrying out smooth noise reduction on the heart rate, the electroencephalogram signals and the speed data by utilizing Matlab software, and removing outliers.
3. The method for determining the minimum safe radius of the flat curve of the two-lane highway in the plateau area as claimed in claim 1, wherein the method comprises the following steps: in the step 2, the electroencephalogram signals refer to beta/theta and beta/(theta + alpha), beta is beta wave in the electroencephalogram signals, theta is theta wave in the electroencephalogram signals, alpha is alpha wave in the electroencephalogram signals, and the heart rate refers to the beating times HR of the human heart per minute; the sample entropy processing is represented by sample entropy SampEn (m, r, N), wherein m is an embedding dimension, r is a similar tolerance, SD is a standard deviation of original data, and N is a data length, so that three electroencephalogram sample entropies SaEn (beta), SaEn (beta/theta), SaEn (beta/(theta + alpha)) and heart rate sample entropies SaEn (HR) of the electroencephalogram signal are obtained respectively.
4. The method for determining the minimum safe radius of the flat curve of the two-lane highway in the plateau area as claimed in claim 3, wherein the method comprises the following steps: in the step 3, the principal component analysis method is to perform principal component analysis on the electroencephalogram sample entropy and the heart rate sample entropy obtained in the step 2, select the first n principal components with the cumulative contribution rate reaching the calculation threshold to determine the weight coefficient of each sample entropy in the integrated evaluation index CISE of the cardiac physiological load, and perform weighted average normalization on the coefficients by taking the variance contribution rate of the principal components as the weight to obtain the integrated evaluation index CISE of the cardiac physiological load;
Figure FDA0003016360010000021
in the formula: a. b, c and d are all weighted coefficients after weighted average normalization.
5. The method for determining the minimum safe radius of the flat curve of the two-lane highway in the plateau area as claimed in claim 1, wherein the method comprises the following steps: in step 4, the different altitude intervals refer to the altitude range (3000 m-5000 m) divided into four altitude intervals (3000m, 3500m), (3500m, 4000m), (4000m, 4500m), (4500m, 5000 m); the lateral force acceleration ahThe expression is as follows:
Figure FDA0003016360010000022
wherein V is the running speed of the smooth curve road section after the smooth noise reduction treatment, R is the radius of the smooth curve, ihWhen the altitude is 10m, g is 9.807m/s2(ii) a H is 3000-3500 m, g is 9.796m/s2(ii) a H is 3500-4000 m, g is 9.794m/s2(ii) a H is 4000-4500 m, g is 9.793m/s2(ii) a H is 4500-5000 m, g is 9.791m/s2；
Overall evaluation index CISE andthe acceleration a of the driver's lateral forcehThe regression model of (a) is abouthThe regression form of the monotonically increasing negative power function curve of (1) is as follows:
CISE＝eah f
in the formula, e and f are regression coefficients, and both e and f are greater than 0.
6. The method for determining the minimum safe radius of the flat curve of the two-lane highway in the plateau area as claimed in claim 1, wherein the threshold R of the radius of the flat curve in the step 5 issafeThe determination steps are as follows:
respectively calculating CISE and a in each altitude interval by utilizing Lagrange median theoremhThe curvature mutation point of the regression model of (1), and a corresponding to the curvature mutation pointhThe value is taken as the maximum lateral force acceleration a which meets the condition that the psychological and physiological loads of the driver are in a safe state in different altitude intervalshmaxThe maximum lateral force acceleration a is calculated using the following formulahmaxCorresponding flat curve radius threshold value Rsafe：
Figure FDA0003016360010000023
Wherein V is the running speed of the smooth curve road section after the smooth noise reduction treatment, g is the gravity acceleration, ihThe road section is a flat curve section ultrahigh value.
7. The method for determining the minimum safe radius of the flat curve of the two-lane highway in the plateau area as claimed in claim 3, wherein m is 2, r is 0.2SD, and N is 256.
CN202110390149.5A 2021-04-12 2021-04-12 Method for determining minimum safe radius of flat curve of double-lane highway in plateau area Active CN113215881B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110390149.5A CN113215881B (en) 2021-04-12 2021-04-12 Method for determining minimum safe radius of flat curve of double-lane highway in plateau area

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110390149.5A CN113215881B (en) 2021-04-12 2021-04-12 Method for determining minimum safe radius of flat curve of double-lane highway in plateau area

Publications (2)

Publication Number Publication Date
CN113215881A CN113215881A (en) 2021-08-06
CN113215881B true CN113215881B (en) 2022-03-08

Family

ID=77087154

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110390149.5A Active CN113215881B (en) 2021-04-12 2021-04-12 Method for determining minimum safe radius of flat curve of double-lane highway in plateau area

Country Status (1)

Country Link
CN (1) CN113215881B (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101818470A (en) * 2010-04-20 2010-09-01 长安大学 Method for optimally setting expressway traffic safety facilities
CN101942790A (en) * 2010-09-10 2011-01-12 天津市市政工程设计研究院 Port area bend linear design method based on ACT-R
CN106650057A (en) * 2016-12-07 2017-05-10 东南大学 Road horizontal curve radius design scheme security assessment method based on virtual vehicle rollover and sideslip tests
CN108182310A (en) * 2017-12-25 2018-06-19 合肥工业大学 Area of heavy rainfull road radius of horizontal curve and safety speed-limit setting method

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101818470A (en) * 2010-04-20 2010-09-01 长安大学 Method for optimally setting expressway traffic safety facilities
CN101942790A (en) * 2010-09-10 2011-01-12 天津市市政工程设计研究院 Port area bend linear design method based on ACT-R
CN106650057A (en) * 2016-12-07 2017-05-10 东南大学 Road horizontal curve radius design scheme security assessment method based on virtual vehicle rollover and sideslip tests
CN108182310A (en) * 2017-12-25 2018-06-19 合肥工业大学 Area of heavy rainfull road radius of horizontal curve and safety speed-limit setting method

Also Published As

Publication number Publication date
CN113215881A (en) 2021-08-06

Similar Documents

Publication Publication Date Title
CN113197576B (en) Method for selecting linear indexes of curve slope combined section of dual-lane highway in plateau area
CN103895649B (en) A kind of driver safety driving warning method
CN109147279B (en) Driver fatigue driving monitoring and early warning method and system based on Internet of vehicles
DE102007011169A1 (en) A method and system for detecting a driver&#39;s skill in handling by the driver&#39;s steering behavior
DE102012207458A1 (en) System and method for vehicle driving style determination
CN106859674A (en) A kind of upland highway driver gives birth to mental load degree characterizing method
CN113112796B (en) Construction method and system of driving behavior characteristics and driving behavior analysis method
CN111311093B (en) Road intersection risk assessment and early warning method based on driver physiological data
CN108171361A (en) Consider the Traffic Flow Simulation Models scaling method of traffic conflict index distribution problem
CN113215881B (en) Method for determining minimum safe radius of flat curve of double-lane highway in plateau area
CN111657923A (en) Method and system for testing danger perception capability of driver
CN103955135A (en) Method for calculating locomotive optimized operation sequence based on double-layer mode curve
CN106326873A (en) Maneuvering intention method employing electromyographic signals of CACC driver&#39;s limbs for representation
CN109649399B (en) Driver distraction judgment method based on transverse driving speed
CN115500843A (en) Sleep stage staging method based on zero sample learning and contrast learning
Lv et al. Compact vehicle driver fatigue recognition technology based on EEG signal
CN112489428B (en) Road evaluation method and system based on driving behaviors
El-Shawarby et al. Evaluation of driver perception–reaction time under rainy or wet roadway conditions at onset of yellow indication
CN115192026B (en) Tunnel driving load monitoring method and terminal
Zhong et al. An expressway driving stress prediction model based on vehicle, road and environment features
CN116579702B (en) Cerebral apoplexy diagnosis ambulance operation monitoring system based on data analysis
Rosipal et al. EEG-based drivers’ drowsiness monitoring using a hierarchical Gaussian mixture model
CN115331431B (en) Pedestrian crossing psychological load test and evaluation method under signal lamp countdown condition
WO2020225970A1 (en) Disease condition determination device, disease condition determination method, program for disease condition determination device, and disease condition determination system
CN109835334B (en) Performance evaluation method and device of lane keeping system

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant