
Citation: Disha Mohan Bangalore, Ingrid Tessmer. Unique insight into protein-DNA interactions from single molecule atomic force microscopy[J]. AIMS Biophysics, 2018, 5(3): 194-216. doi: 10.3934/biophy.2018.3.194
[1] | Vincent Renault, Michèle Thieullen, Emmanuel Trélat . Optimal control of infinite-dimensional piecewise deterministic Markov processes and application to the control of neuronal dynamics via Optogenetics. Networks and Heterogeneous Media, 2017, 12(3): 417-459. doi: 10.3934/nhm.2017019 |
[2] | Simone Göttlich, Stephan Martin, Thorsten Sickenberger . Time-continuous production networks with random breakdowns. Networks and Heterogeneous Media, 2011, 6(4): 695-714. doi: 10.3934/nhm.2011.6.695 |
[3] | Yilun Shang . Group pinning consensus under fixed and randomly switching topologies with acyclic partition. Networks and Heterogeneous Media, 2014, 9(3): 553-573. doi: 10.3934/nhm.2014.9.553 |
[4] | Ciro D'Apice, Olha P. Kupenko, Rosanna Manzo . On boundary optimal control problem for an arterial system: First-order optimality conditions. Networks and Heterogeneous Media, 2018, 13(4): 585-607. doi: 10.3934/nhm.2018027 |
[5] | Xianmin Geng, Shengli Zhou, Jiashan Tang, Cong Yang . A sufficient condition for classified networks to possess complex network features. Networks and Heterogeneous Media, 2012, 7(1): 59-69. doi: 10.3934/nhm.2012.7.59 |
[6] | Charles Bordenave, David R. McDonald, Alexandre Proutière . A particle system in interaction with a rapidly varying environment: Mean field limits and applications. Networks and Heterogeneous Media, 2010, 5(1): 31-62. doi: 10.3934/nhm.2010.5.31 |
[7] | Michael Herty, Lorenzo Pareschi, Sonja Steffensen . Mean--field control and Riccati equations. Networks and Heterogeneous Media, 2015, 10(3): 699-715. doi: 10.3934/nhm.2015.10.699 |
[8] | Guy Barles, Emmanuel Chasseigne . (Almost) Everything you always wanted to know about deterministic control problems in stratified domains. Networks and Heterogeneous Media, 2015, 10(4): 809-836. doi: 10.3934/nhm.2015.10.809 |
[9] | Jésus Ildefonso Díaz, Tommaso Mingazzini, Ángel Manuel Ramos . On the optimal control for a semilinear equation with cost depending on the free boundary. Networks and Heterogeneous Media, 2012, 7(4): 605-615. doi: 10.3934/nhm.2012.7.605 |
[10] | Karoline Disser, Matthias Liero . On gradient structures for Markov chains and the passage to Wasserstein gradient flows. Networks and Heterogeneous Media, 2015, 10(2): 233-253. doi: 10.3934/nhm.2015.10.233 |
The PHEV has the characteristics of traditional hybrid electric vehicle and pure electric vehicle, which is an important development direction of new energy vehicles [1]. In order to achieve better energy-saving effect of the PHEV, energy management strategy (EMS) is needed to determine the working mode of various components of the power transmission system and the energy distribution of different power sources to adapt to diverse driving environments [2].
In terms of EMS goals for PHEV, some studies focus on battery aging and its thermal dynamics [3,4], and most studies is devoted to reducing fuel consumption. The way to achieve EMS can be divided into two types: rule-based and optimization-based [5,6]. The rule-based strategy is mainly based on the engineering experience and the characteristics of each component in the power system to design the working mode of the system and the energy distribution of different power sources. The rules are versatile and practical, but the potential for energy saving is limited [6,7,8,9]. Optimization-based strategy can be divided into instantaneous optimization, local optimization and global optimization.
The equivalent consumption minimization strategy (ECMS) is the representative of the instantaneous optimization strategy [7,10,11]. ECMS transforms the electric energy consumption into the equivalent fuel consumption by introducing the EF, so as to optimize the instantaneous equivalent fuel consumption [7,12,13]. This strategy has strong real-time performance, but poor robustness and limited optimization effect.
The EMS of local optimization is mainly based on model predictive control (MPC). This strategy calculates the optimal control parameters in the prediction domain, with small amount of calculation, real-time control and significant optimization effect [8,9].
The global optimization strategy is represented by Reinforcement Learning (RL) and DP. Learning-based EMS has a strong self-learning ability and it is model-free which relies on a mass of real-world sample for training [14,15]. DP can obtain the multi-stage optimization decision and the theoretical global fuel economy optimization, if the global driving cycle is known beforehand. However, the disadvantage is that the calculation amount is large, and the future driving cycle is difficult to obtain in advance, which limits the use of DP in real time [16,17].
In order to apply approximately optimization decision by DP for real-time control, future driving cycles need to be known, which can be obtained by two ways: one is to build a future speed prediction model based on the historical driving cycle data by neural network, Markov or Kalman filter [9,18,19]. And the other is to obtain the traffic information by high-tech [11,16,17], such as Navigation System, Global Positioning System (GPS), Geographic Information System (GIS), Internet of Vehicles and Intelligent Transportation System (ITS). Reference [16] combines GPS and ITS to obtain driving information, and adjusts the EF in A-ECMS accordingly.
The PHEV controller can plan the optimization decision by DP offline to adjust the control parameters in real-time [11]. In reference [17], an EMS of hierarchical control is proposed. The upper layer obtains real-time traffic data with GPS and ITS, and makes long-term planning of SOC optimal consumption trajectory, while the lower layer uses MPC to realize short-term time-domain speed prediction. Some research build A-ECMS based on DP combining initial EF optimized by genetic algorithm [20], or combining driving cycle recognition [21].
When using DP to optimize EMS, it is necessary to recognize driving cycle for more targeted controlling of various driving cycles. In recent years, many researchers [22,23,24] have made great efforts to study the driving cycles of typical cities or specific vehicles. K-means cluster analysis [25], fuzzy c-means (FCM) [26], and neural network [27] are common methods to construct typical driving cycle by historical driving data. Driving cycles directly affect fuel consumption, so how to optimize the energy management of PHEV in combination with driving cycles will be a research trend [17]. Some research uses existing cycles to optimize the EMS, but few combine personalized travel characteristics for a specific driver.
Based on the above studies, it can be seen that current PHEV energy management strategy has achieved multiple research results, a variety of optimization methods have been applied to PHEV energy management. How to make the strategy to get better adaptability to driving cycles and enough real-time optimization ability need to be further resolved.
As discussed above, how to combine ITS to globally optimize PHEV energy management is the research focus. However, long-term traffic flow information from ITS can only reflect the average speed and cannot be adjusted according to the actual speed, lacking robustness. Few studies combine ITS and short-term speed prediction for EMS. Therefore, the main purpose of this study is to improve the performance of A-ECMS by tuning EF instantaneously based on driving cycle prediction and personalized travel characteristics. Two contributions have been realized in this paper, which is: An A-ECMS was proposed, with an EF jointly correction algorithm combining long-term target driving cycle and short-term vehicle speed prediction. Moreover, the DP algorithm is applied to the offline energy optimization process of A-ECMS based on typical driving cycles constructed according to personalized travel characteristics. The improved A-ECMS formulated in this paper can optimize EF in real time based on mileage, SOC, long-term driving cycle and real-time vehicle speed. Meanwhile, the online driving cycle recognition can be practically applied, and the adaptability of A-ECMS to actual driving can be enhanced.
This paper is organized as follows. In Section 2, an A-ECMS with driving cycle prediction and personalized travel characteristics merging DP algorithm is described. In Section 3, the methodology of building typical driving cycles for a specific driver and planning SOC reference trajectory algorithm is introduced. Typical driving cycle construction and simulation verification of A-ECMS are in Section 4. And conclusions are subsequently discussed in Section 5.
In this paper, the P2 configuration PHEV is taken as the research object. This configuration refers to the coaxial connection of the engine and the motor, while P2 means that the motor is placed at the front end of the transmission, between the engine and the transmission, and using the clutch between the engine and the motor to control the power coupling and disconnection. The structure diagram is shown in Figure 1 and the main parameters are shown in Table 1.
Component | Parameter value |
Engine | Peak power: 64kW; Maximum torque :153Nm |
Motor | Peak power: 32Kw; Maximum torque 160Nm |
Battery capacity | 27.9Ah |
Rated voltage of battery | 320V |
Final drive ratio | 3.944 |
Transmission speed ratio | 3.538/2.045/1.333/1.028/0.82 |
Vehicle | Mass: $ m=1568 $ kg; Vehicle equivalent cross section: $ A=2.13{m}^{2} $; Aerodynamic drag coefficient: $ {C}_{d}=0.29 $; Tire radius: $ r=312mm $ |
The characteristics of the engine and the motor are shown in Figure 2 and Figure 3 respectively, and the SOC calculation is realized by using the ampere hour integration method [20].
DP, a global optimization algorithm, can calculate the optimal solution of PHEV energy management problem, so as to obtain the optimal fuel economy, but the premise is to know the global driving cycle [12].
In this paper, SOC of battery is selected as the system state variable, motor torque $ {T}_{m} $ and transmission gear $ gear\left(k\right) $ are selected as the system control variables, then the state and control quantity of the system can be given as
$\left\{ {x(k)=SOC(k)u(k)=[Tm(k),gear(k)] } \right.$
|
(1) |
According to the dynamic relationship between power assemblies, the engine torque $ {T}_{e} $ obtained from the motor torque $ {T}_{m} $ can be given as follows
${T_e} = {T_r}/({i_0} \cdot {i_{gear(k)}}) - {T_m}$ | (2) |
where $ {T}_{r} $ is the required torque of the whole vehicle, $ {i}_{0} $ is the main reduction ratio, and $ {i}_{gear\left(k\right)} $ is the transmission ratio of the current gear.
Since the selected state variable is the SOC of battery, the state equation under the dynamic programming algorithm can be given as follows
$SOC(k + 1) = SOC(k) - \frac{{{U_{oc}}(k) - \sqrt {U_{oc}^2(k) - 4 \cdot {R_0}(SOC(k)) \cdot {P_{bat}}(k)} }}{{2 \cdot {R_0}(SOC(k)) \cdot {Q_{bat}}}}$ | (3) |
where $ {U}_{oc} $ is the open circuit voltage, $ {R}_{0} $ is the internal resistance of the battery, $ {Q}_{bat} $ and $ {P}_{bat} $ are the battery capacity and the battery power, respectively.
Since the system takes fuel consumption as the objective function, the objective function under the DP algorithm can be given as
${J^*} = \min J = \min \left\{ {\sum\limits_{k = 1}^N {L(x(k), u(k))} } \right\} = \min \left\{ {\sum\limits_{k = 1}^N {fue{l_k} + \alpha {{(SOC(k) - SO{C_f})}^2}} } \right\}$ | (4) |
where $ SO{C}_{f} $ is the SOC value at the termination time, $ fue{l}_{k} $ is the fuel consumption of the whole vehicle at the current stage. And $ \alpha \left(SOC\right(k)-SO{C}_{f}{)}^{2} $ is the penalty function of SOC in DP, $ \alpha $ is the penalty coefficient, which is to correct the SOC value of battery and the target SOC value at the termination time when DP calculation ends.
At the same time, when DP algorithm is applied to solve the PHEV energy management problem, it will be limited by the hybrid system itself, so some constraints in the iterative calculation process are shown as follows
$\left\{ {ne_min⩽ne(k)⩽ne_maxTe_min⩽Te(k)⩽Te_maxnm_min⩽nm(k)⩽nm_maxTm_min⩽Tm(k)⩽Tm_maxSOCmin⩽SOC(k)⩽SOCmax } \right.$
|
(5) |
where $ {n}_{e}\left(k\right) $ is the engine speed, $ {n}_{m}\left(k\right) $ is the motor speed, $ {n}_{e\_\mathit{min}\mathrm{ }} $ and $ {n}_{e\_\mathit{max}\mathrm{ }} $ are the engine minimum and maximum speed respectively; $ {T}_{e\_\mathit{min}\mathrm{ }} $ and $ {T}_{e\_\mathit{max}\mathrm{ }} $ are the engine minimum and maximum torque respectively; $ {n}_{m\_\mathit{min}\mathrm{ }} $ and $ {n}_{m\_\mathit{max}\mathrm{ }} $ are the motor minimum and maximum speed respectively; $ {T}_{m\_\mathit{min}\mathrm{ }} $ and $ {T}_{m\_\mathit{max}\mathrm{ }} $ are the motor minimum and maximum torque respectively; $ SO{C}_{\mathit{min}} $ and $ SO{C}_{max} $ are the upper and lower limits of battery SOC respectively.
In order to study the optimal SOC trajectory under different driving mileage of PHEV based on DP, this paper selects the different number NEDC driving cycles as test cycles. In the strategy, the initial value of the battery SOC is set to 0.9 and the target value is 0.3. Respectively, simulating under 1 to 8 NEDC driving cycles and result is shown in Figure 4.
The driving range is a key factor affecting the EMS for PHEV. As can be seen from Figure 4, the trajectory of battery SOC will vary with the distance, and SOC presents an approximately linear decline curve. When the trip ends, the SOC value drops to near the target value. But DP algorithm requires a large amount of calculation, so it cannot be applied to real-time control. However, SOC linear descent curve optimized by DP can provide a basis for reference SOC trajectory planning in EMS.
Taking multiple NEDC cycles as an example above, it is found that the optimal SOC reference trajectory is an approximate straight line related to mileage. However, it is difficult to use the SOC trajectory under NEDC or a certain existing driving cycle as a reference for actual driving condition. Therefore, if cluster analysis is performed according to the historical driving data of a specific driver, $ n $ typical driving cycles can be obtained, which can better reflect the personalized travel characteristics of the driver. Then studying the optimal SOC consumption under each typical driving cycle by DP.
As is shown in Figure 5, according to the characteristic parameters, the driving condition is segmented and corresponding classified into different typical driving cycles. The SOC reference trajectory can be linearly composed of the optimal SOC consumption of each section. The specific allocation method and analysis are described in Section 3.2.
SOC reference trajectory described above determines the target driving cycle in advance which requires the use of ITS to generate. This paper extracts the traffic flow based on ITS, generates the target driving cycle of the travel road as the traffic information of the road section, and obtains the average traffic flow speed of the front road section in real time [28].
The actual navigation route is segmented every 200m, and the average traffic speed of each section can be calculated as follows
${V_{seg(i)}}(t) = \frac{{\sum\nolimits_{k = 1}^N {V_k^{seg(i)}(t)} }}{N}$ | (6) |
where $ {V}_{seg}\left(i\right) $ is the average speed of the section marked $ i $ at time $ t $, $ {V}_{k}^{seg\left(i\right)} $ is the average speed of the $ k $ th vehicle passing through the section at time $ t $, and $ N $ is the total number of vehicles recorded in the section at that time.
Based on the above method, the route to the destination and the traffic flow speed information can be obtained by the navigation system and ITS, and then the possible vehicle speed information in the next trip can be estimated by fitting according to the traffic flow speed of the sub sections. For a certain section of route, the curve of the average speed and the mileage of each section is shown in Figure 6.
However, the road network only obtains the relationship between mileage and vehicle speed, it also needs to be converted to the time domain, as shown in Eqs. (7).
${T_{seg(i)}} = \frac{{{L_{seg(i)}}}}{{{V_{seg(i)}}}}$ | (7) |
where $ {L}_{seg\left(i\right)} $ and $ {T}_{seg\left(i\right)} $ are the length and the time of passing through section $ i $.
After conversion to time domain, the speed-time curve is a step curve, which does not conform to the continuity characteristics of the vehicle speed, so it cannot be regarded as the real driving cycle. Therefore, this paper uses the LOWESS filter function to smooth the step speed signal, so as to convert it into a driving cycle curve that conforms to the normal driving speed rule, as shown in Figure 7.
In the energy management optimization strategy of PHEV, this is the most likely driving cycle of the vehicle. The speed curve can be taken as the target driving cycle for long-term prediction.
The long-term target driving cycle based on traffic flow information is the global predicted driving cycle of the travel route, but there may be deviation in the real-time process, and the short-term vehicle speed prediction is more consistent with the real-time driving cycle. Therefore, in the prediction of PHEV driving cycle, the global prediction and short-term prediction can be combined.
The NAR neural network has advantages of solving time variant and nonlinear problems, making it suitable for vehicle speed prediction [29]. Therefore in this paper, NAR neural network with feedback is built to predict the speed, the structure diagram is shown in Figure 8. Taking the historical driving cycle data as the input sample, through trial-and-error method and repeated tests, the number of neurons in the hidden layer is determined to be 20, and the delay order of the output feedback is determined to be 35 [29,30]. The neural network is trained and the training results are obtained. Figure 9 is the training process, and the ordinate is the mean square error (MSE). It can be seen from the figure that when the number of training reaches 16 times, the mean square error of the verification set is the smallest and converges to 0.0019219, so the training ends.
Figure 10 is an error autocorrelation graph of NAR neural network, which describes the correlation of prediction errors in time. Theoretically, the autocorrelation function should have only one maximum value and at Lag = 0, indicating that the prediction error is completely uncorrelated. If there is a significant correlation in the prediction error, the training effect of the network is not good [31]. It can be seen from Figure 10 that the NAR network in this paper has the largest correlation at Lag = 0, and the correlation of other points is within the confidence interval, indicating that the training results and prediction performance of the network are good at this time.
It can be seen from the Figure 11 that most of the error values of single-step prediction are in the interval of [-0.2, 0.2], and the error is very small, indicating that the single-step prediction results are accurate and the selection of neurons number is reasonable.
Then the prediction effect of the neural network is studied as follows. Root mean square error (RMSE) is used to evaluate the accuracy of NAR neural network and follow the rule of the higher the RMSE value, the lower the prediction accuracy of the model. In order to evaluate the prediction effects on speed, the results in different prediction duration are shown in Figure 12,and the corresponding RMSE values are shown in Table 2.
Prediction durations | 5s | 10s | 15s | 20s |
RMSE | 2.6456 | 5.5578 | 7.8314 | 10.0521 |
It can be seen from Figure 12 and Table 2 that when the prediction time domain is 5 seconds, the overall prediction effect of the test cycle is good, and its RMSE value is also small. With the increase of prediction time domain from 5 to 20 seconds, RMSE increased significantly, while the prediction accuracy decreased gradually. On the whole, vehicle speed prediction time domain is selected for 5 seconds.
The essence of ECMS is to convert the electric energy consumption by the battery into the equivalent fuel consumption through the equivalent factor in each instant. The energy distribution rule with the minimum sum of the instantaneous fuel consumption of the engine and the equivalent fuel consumption of the motor as the optimal solution is conducive to reducing the energy consumption and improving the economy [13].
The instantaneous equivalent fuel consumption of PHEV at a certain time $ t $ can be given as follows
${\dot m_{eq}}(t) = {\dot m_{fuel}}(t) + {\dot m_{ele}}(t)$ | (8) |
where $ {\dot{m}}_{eq} $ is the total equivalent fuel consumption; $ {\dot{m}}_{fuel} $ is the fuel consumption of the engine; $ {\dot{m}}_{ele} $ is the equivalent fuel consumption after the electric energy conversion of the power battery, and its calculation formula is shown as follows
${\dot m_{ele}}({P_m}(t)) = k \cdot s \cdot \frac{{{P_m}(t)}}{{{\eta _{dis}}{Q_{lhv}}}} + (1 - k) \cdot s \cdot \frac{{{\eta _{char}}{P_m}(t)}}{{{Q_{lhv}}}}$ | (9) |
where $ k = 0.5\times \left(1+sign\left({P}_{m}\left(t\right)\right)\right) $, indicating whether the working state of the motor is discharging or charging; $ {P}_{m}\left(t\right) $ is the power of the motor at time $ t $, $ {Q}_{lhv} $ is the calorific value constant of the gasoline mass, $ {\eta }_{dis} $ and $ {\eta }_{char} $ are the discharge and charge efficiency of the battery, and $ s $ is oil-electricity conversion equivalent factor (EF).
According to the energy management optimization problem of PHEV, the objective function can be obtained as follows
$\min J(t) = \int_{{t_1}}^{{t_n}} {{{\dot m}_f}} (u(t))dt$ | (10) |
where $ J\left(t\right) $ is the fuel consumption and $ u\left(t\right) $ is the control variable at time $ t $.
For the ECMS strategy, the EF is fixed. However, the driving cycle information and the SOC of the battery also affect the EF. If a fixed EF is adopted globally, the optimization effect may be poor [32,33]. Based on ECMS, an adaptive control strategy based on real-time correction of EF by SOC reference trajectory and vehicle speed is proposed in this section. An adaptive equivalent consumption minimization strategy (A-ECMS) is established by introducing the SOC reference trajectory and short-term driving speed prediction. The formula for changing the EF can be given as follows
$s(t) = {s_{ref\_i}} + {K_p} \cdot (SO{C_{ref}} - SOC(t))$ | (11) |
where $ {s}_{ref\_i} $ which is the initial equivalent factor of various typical driving cycles under the current vehicle speed, $ {K}_{p} $ is the scale factor, $ SO{C}_{ref} $ is the reference SOC value, and $ SOC\left(t\right) $ is the real SOC value at time $ t $.
SOC reference trajectory can modify the EF in real time through the proportional control method. Constant part of the adaptive equivalent factor $ {s}_{ref\_i} $, or initial EF, is also very important in the method. Therefore, this paper considers the combination of short-term vehicle speed prediction to optimize initial EF of four typical driving cycles under different SOC values. In the application, the constant part of the adaptive EF is modified according to the type of driving cycle, and compared with the EF calculated based on SOC reference trajectory. If there is too much difference between the predicted EF and the current EF, then the constant part of the current EF is modified.
In summary, A-ECMS based on driving cycle prediction and personalized travel characteristics is proposed in this paper and its process is shown in Figure 13. First of all, the driver's historical driving data is clustered and analyzed offline to extract personalized travel characteristics and build the typical driving cycles. Then real-time traffic flow information is obtained through ITS to plan the optimal SOC reference trajectory, and combining the short-term vehicle speed prediction to modify the EF, so as to control the vehicle to achieve the purpose of improving fuel economy.
The development and application of EMS should be able to select appropriate control parameters according to the actual driving cycles of the vehicle, so as to distribute the energy of each power source more reasonably, achieve the best control effect and save energy as much as possible [28,33]. Meanwhile, driving style of each driver and actual driving conditions are diverse. This section will explain the methodology of building typical driving cycles by analyze personalized travel characteristics, and planning the SOC trajectory by driving cycle recognition.
In order to make the PHEV energy management strategy more targeted, the historical driving cycle data by a specific driver should be analyzed firstly, and building the typical driving cycles for the driver which is the basis of SOC reference trajectory analysis. The analysis process is shown in Figure 14.
The data of the driver's travel is collected by GPS equipment, and the sampling frequency is 1Hz. And preprocessing is needed to remove the excessively long idle speed section and abnormal section in the data. In the process of driving on urban roads, the driving cycles of vehicles are usually affected by different types of road conditions and traffic environment. There are many idle, acceleration, deceleration and other states [34]. Therefore, after preprocessing the historical data, the next important step is to divide the kinematic segments, also known as short segments. As shown in Figure 15, the kinematic segment is defined as the operating range from one idle to the next [35,36].
By analyzing the driving characteristic parameters of each kinematic segment, the travel characteristics of the vehicle are obtained. There are 13 characteristic parameters representing the characteristics of driving cycles to be selected, as shown in Table 3.
Numbers | Driving cycle characteristic parameter | Unit |
1 | Average speed:$ {v}_{ave} $ | km/h |
2 | Average driving speed:$ {u}_{ave} $ | km/h |
3 | Standard deviation of vehicle speed:$ {v}_{std} $ | km/h |
4 | Maximum speed:$ {v}_{\mathit{max} } $ | km/h |
5 | Average acceleration:$ {a}_{ave} $ | m/s2 |
6 | Average deceleration:$ {a}_{dec} $ | m/s2 |
7 | Maximum acceleration:$ {a}_{\mathit{max} } $ | m/s2 |
8 | Maximum deceleration:$ {a}_{\mathit{min} } $ | m/s2 |
9 | Standard deviation of acceleration:$ {a}_{std} $ | m/s2 |
10 | Idle time ratio:$ P\left({t}_{idle}\right) $ | % |
11 | Acceleration time ratio:$ P\left({t}_{acc}\right) $ | % |
12 | Deceleration time ratio:$ P\left({t}_{dec}\right) $ | % |
13 | Uniform time ratio:$ P\left({t}_{uni}\right) $ | % |
If the 13 characteristic parameters are directly used for clustering analysis, the calculation will be large and affect the clustering results, and there will be information redundancy between the various feature parameters. Therefore, it is necessary to reduce the dimension of the original driving cycle matrix. In this paper, principal component analysis (PCA) method is used for extracting the key ingredients and reducing the dimensionality.
Then by analyzing the principal component of the kinematic segments, they can be clustered into several representative driving cycles for optimization of EMS. In this paper, the fuzzy c-means algorithm (FCM) is used for clustering analysis. If samples $ X = [{x}_{1}, {x}_{2}, ..., {x}_{n}] $ are given, and $ c $ represents the number of clusters, $ [{\mathrm{\alpha }}_{1}, {\mathrm{\alpha }}_{2}, ..., {\mathrm{\alpha }}_{C}] $ represents cluster centers of various categories. Then the objective function can be expressed by Eqs. (12).
${J_b}(U, \alpha ) = \sum\limits_{i = 1}^n {\sum\limits_{k = 1}^c {{{({\delta _{ik}})}^b}{{({d_{ik}})}^2}} } $ | (12) |
where $ b $ is weighting parameters. $ {\delta }_{ik} $ and $ {d}_{ik} $ are membership and Euclidean distance between $ i $ th sample and $ k $ th category. And Euclidean distance can be calculated by Eqs. (2). $ m $ is the dimensions of the sample.
${d_{ik}} = {\left[ {\sum\limits_{j = 1}^m {{{({x_{ij}} - {\alpha _{kj}})}^2}} } \right]^{1/2}}$ | (13) |
Membership between $ i $ th sample and $ k $ th category can be calculated by Eqs. (3).
${\delta _{ik}} = {{1\;} \mathord{\left/ {\vphantom {{1\;} {\sum\limits_{j = 1}^c {{{(\frac{{{d_{jk}}}}{{{d_{ik}}}})}^{\frac{2}{{b - 1}}}}} }}} \right. } {\sum\limits_{j = 1}^c {{{(\frac{{{d_{jk}}}}{{{d_{ik}}}})}^{\frac{2}{{b - 1}}}}} }}$ | (14) |
However, if the initial value of FCM is unreasonable, the algorithm may converge to the local minimum point, which will affect the clustering results [25,37]. To solve this problem, this paper applies simulated annealing algorithm (SA) and genetic algorithm (GA) to FCM, making clustering more efficient and accurate [25]. The process of FCM based on genetic simulated annealing algorithm (GA-SA) is shown in Figure 16.
The type of vehicle driving cycle has a great influence on the energy management strategy of PHEV. Therefore, the algorithm of driving cycle recognition is designed based on the results of above clustering analysis, in order to identify the current type of driving cycle in the process of vehicle driving, and provide help for the later SOC reference trajectory planning. The driving cycle recognition process is shown in Figure 17.
The main idea of driving cycle recognition based on cluster analysis is to extract the characteristic parameters of the driving cycle during the actual driving process, to calculate the principal component and the Euclidean proximity of the driving cycle to the clustering center of various typical driving cycles, and to classify the driving cycle as the cluster of a certain category if it is closest to the clustering center of this category [34]. Thus, the classification of the current driving cycle can be recognized.
Let the cluster center of typical driving cycles be $ {Y}_{n}(n = \mathrm{1, 2}, \dots, n) $ and the current actual driving cycle of the vehicle to be identified is $ X $. Then the Euclidean proximity between the current actual driving cycle and the typical driving cycle in the historical driving cycle database can be given as follows
$\sigma ({Y_n}, X) = 1 - \frac{1}{{\sqrt m }}{(\sum\limits_{k = 1}^m {({Y_n}(k) - X(k)} )^2}{)^{1/2}}$ | (15) |
where $ m $ is the number of characteristic parameters of the driving cycle.
According to the Euclidean proximity formula, the approach degree between $ {Y}_{n} $ and $ X $ can be calculated. If there is
$\sigma ({Y_n}, X) = \max \left\{ {\sigma ({Y_1}, X), \sigma ({Y_2}, X), \sigma ({Y_3}, X), \sigma ({Y_4}, X)} \right\}$ | (16) |
It can be considered that the current driving cycle $ X $ and $ {Y}_{n} $ have the highest similarity, that is, the current driving cycle belongs to this kind of driving cycle.
After classifying the driving cycle fragments into typical driving cycles, DP algorithm is needed to solve the optimal SOC consumption law of various typical driving cycles in this paper. By setting different initial value and end value of SOC in the DP algorithm, study the relationship between the SOC consumption per kilometer and the high efficiency utilization rate of engine under typical driving cycles, and the comparison between the SOC consumption per kilometer and the average fuel consumption rate. The calculation of high efficiency area utilization rate is given as follows [21].
${\eta _{eng\_efficient}} = \frac{{b{e_{\min }}}}{N} \cdot \sum\limits_{i = 1}^N {\frac{1}{{be(i)}}} \cdot 100\% $ | (17) |
where $ N $ represents the length of time when the engine is working, and $ {be}_{min} $ represents the minimum fuel consumption rate of the engine.
If the whole journey is within the pure electric mileage of PHEV, the pure electric mode is used for the journey without SOC trajectory planning. But if it exceeds the pure electric mileage that can be driven under the current battery state, in order to improve fuel economy, reasonable SOC planning is required. The optimal SOC consumption of each typical driving cycle can be determined according to the relationship between SOC, fuel consumption, and high efficiency utilization rate of engine, so as to formulate SOC trajectory planning rule reasonably. This rule needs to traffic flow information in advance which can be obtained by navigation system and ITS, and then allocate power to each section based on the type of typical driving cycle and the current SOC state. The process of SOC trajectory planning is shown in Figure 18.
The driving data of a specific driver is collected by GPS and the sampling frequency is 1Hz. The data of vehicle speed and time information is extracted, which is shown in the Figure 19.
There are some abnormal acceleration and deceleration data points caused by GPS equipment or vehicle, abnormal points with discontinuous time, excessive idle data points. These data points need to be eliminated and the historical driving cycle data after preprocessing is shown in Figure 20.
1468 segments of driving cycle are extracted from historical data. 13 driving cycle characteristic parameters, such as average speed, average driving speed and standard deviation of speed, are respectively represented by $ {x}_{1}, {x}_{2}, \dots, {x}_{13} $, and the index of driving cycle segment is represented by $ i = \mathrm{1, 2}, \dots, 1468 $, while the characteristic parameter $ {x}_{1}, {x}_{2}, \dots, {x}_{13} $ of the $ i $ th driving cycle segment is recorded as $ [{x}_{i1}, {x}_{i2}, ..., {x}_{i13}] $. Then there is the full information matrix $ X = ({x}_{ij}{)}_{1468\times 13} $ of the actual driving cycle. Then using PCA regroups the original parameters into a new set of unrelated integrated parameters to reduce dimension and the result is shown in Figure 21.
As a result, the accumulative contribution rate of the first four principal components reached 91.18%. Generally speaking, when it is higher than 85%, it can be considered that the number of principal components selected can retain most of the original data information, and achieve the purpose of dimensionality reduction [34]. Therefore, the first four principal components are selected in this paper. The dimension of historical driving cycle data has been reduced from 13 dimensions to 4 dimensions, which greatly reduces the calculation amount. The principal component score matrix of driving cycle segments obtained by principal component analysis can be used for cluster analysis, as shown in Table 4.
Index | F1 | F2 | F3 | F4 |
1 | -1.1058 | -1.3408 | -0.2677 | 1.9947 |
2 | 0.8687 | -1.7607 | -0.2811 | 0.8793 |
3 | -0.1886 | -1.6428 | 0.5574 | 0.7264 |
4 | -0.7973 | -1.8836 | 2.0468 | 0.0025 |
⋮ | ⋮ | ⋮ | ⋮ | ⋮ |
1467 | 3.6596 | -0.3058 | 0.8410 | -1.4611 |
1468 | 3.0708 | -0.3280 | -0.5293 | -0.1305 |
Four principal components of each kinematic segment have obtained above which can be used as the feature vector of FCM samples. And four cluster centers are selected. Iterative process of clustering is shown in Figure 22. It can be seen the convergence speed and objective function value of GASA-FCM are both better than FCM.
In order to make the clustering results intuitive and clear, the two-dimensional graph composed of principal component 1 and 2 and the three-dimensional graph composed of principal component 1, 2 and 3 are described to show the classification situation, as shown in Figure 23. The black fork in the graph is cluster center.
Then it is necessary to extract typical segments of four types of driving cycles. According to the principle that 10 groups of short segments of driving cycles closest to the clustering center are selected as candidate characteristic driving cycles, four types of typical driving cycles are selected, as shown in Figure 24.
Type-1 driving cycle can be defined as medium speed driving cycles. The vehicle speed and the idle time are both moderate, which is generally a relatively smooth urban driving cycle. Type-2 driving cycle can be defined as high-speed driving cycle. The average speed is the highest of the four driving cycle types. Meanwhile, the idle time ratio is the lower of the four driving cycle types, and the maximum speed is also higher, indicating that the road condition of this kind of driving cycle is smooth. Generally, it belongs to the smooth urban expressway section or suburban section. Type-3 driving cycle can be defined as urban medium and low speed driving cycle. It belongs to the urban smooth road condition, which is reflected in the low average driving speed, low idle ratio, and little difference in the proportion of acceleration state, deceleration state and uniform state. Type-4 of driving cycle can also be defined as urban medium and low speed driving cycle. However, the road condition is congested and the start and stop are frequent, which is reflected in the low average driving speed. The idle state accounts for the highest proportion among the four states of acceleration, deceleration, uniform speed and idle.
The comparison between the SOC consumption per kilometer and the high efficiency utilization rate of engine under four typical driving cycles, and the comparison between the SOC consumption per kilometer and the average fuel consumption rate are in Figure 25 and Figure 26.
The main purpose of calculating the relationship among SOC consumption per kilometer, engine efficient area utilization rate and average fuel consumption rate under various driving cycles is to provide reasonable SOC distribution value for each driving cycle during SOC trajectory planning. At the same time, calculate the SOC consumption value per kilometer of vehicle under pure electric driving under the condition of sufficient power, as shown in Table 5.
Driving cycle type | 1 | 2 | 3 | 4 |
SOC distribution value per km(%/km) | 1.3229 | 1.3643 | 1.2180 | 1.0602 |
It can be seen that Type-1 typical driving cycles, when the SOC consumption reaches 1.1759% per kilometer, the utilization ratio of the engine efficient area is still over 90%, and the average fuel consumption rate is at the critical point of sudden increase, so this point can be regarded as the critical point of the reasonable SOC distribution range. At the same time, the highest point of high efficiency utilization rate of engine is selected as the optimal SOC distribution value under this driving cycle, which is also at the lowest point of fuel consumption, so $ \Delta SO{C}_{opt\_cycle\_1} = \text{0.7349}\mathrm{\%}/km $.
For Type-2 typical driving cycle, the average speed of this driving cycle is higher, and the start and stop are less. It can be seen the utilization rate of the high-efficiency area of the engine is not high as a whole. In Figure 23, the utilization rate of the high-efficiency area of the engine suddenly rises because it is about to reach the critical value of the pure electric mode, and the engine has little work at this time. Moreover, the change trend of the average fuel consumption rate is not obvious.
For Type-3 typical driving cycle, the critical value shall be selected according to the comprehensive consideration so that the engine efficiency can be as high as possible and the average fuel consumption should not be too high, so the optimal SOC value of this type of driving cycle is $ \Delta SO{C}_{opt\_cycle\_3}=\text{0.8526}\mathrm{\%}/km $.
For Type-4 typical driving cycle, this driving cycle belongs to the congestion driving cycle, with many starts and stops, and the speed is not high, which is not conducive to the engine work. It can be seen engine working efficiency of this kind of driving cycle is not high and the fuel consumption is increased. Therefore, try to operate in pure electric mode under this driving cycle. In conclusion, the reasonable distribution range of SOC under various typical driving cycles is shown in Table 6.
Driving cycle type | 1 | 2 | 3 | 4 |
SOC max distribution value per km(%/km) | 1.1759 | 1.2993 | 0.9744 | 1.0602 |
Optimal SOC max distribution value per km(%/km) | 0.7349 | 0.3898 | 0.8526 | 1.0602 |
SOC min distribution value per km(%/km) | 0 | 0 | 0 | 1.0602 |
Generally speaking, when the speed is low and the engine starts and stops frequently, the motor should be allowed to work first, because in this case, it is not conducive to the engine working in the efficient area. Therefore, according to the priority order of typical driving cycle 4-3-1-2, the battery power is allocated. This SOC reference trajectory planning rule needs to obtain the road condition information of the trip in advance, and based on the SOC state at this time, then the power allocation of the sub sections can be carried out. The SOC trajectory planning algorithm is shown in Table 7.
1: Identify the road section belonging to Type-4 typical driving cycle; |
2: Allocate the SOC of the roads of Type-4 typical driving cycle:$ \mathrm{\Delta }SO{C}_{cycle\_4}=\mathrm{\Delta }SO{C}_{EV\_cycle\_4}\times {S}_{cycle\_4} $; |
3: Identify the road section belonging to Type 3 of typical driving cycles; |
4:Allocate the SOC of the roads of Type-3 typical driving cycle:$ \mathrm{\Delta }SO{C}_{cycle\_3}=\mathrm{\Delta }SO{C}_{opt\_cycle\_3}\times {S}_{cycle\_3} $; |
5:Identify the road section belonging to Type-1 of typical driving cycles; |
6:Allocate the SOC of the roads of Type-1 typical driving cycle:$ \mathrm{\Delta }SO{C}_{cycle\_1}=\mathrm{\Delta }SO{C}_{opt\_cycle\_1}\times {S}_{cycle\_1} $; |
7:Identify the road section belonging to Type 2 of typical driving cycle; |
8:Allocate the SOC of the roads of Type-2 typical driving cycle:$ \mathrm{\Delta }SO{C}_{cycle\_2}=\mathrm{\Delta }SO{C}_{total}-(\mathrm{\Delta }SO{C}_{cycle\_1}+\mathrm{\Delta }SO{C}_{cycle\_3}+\mathrm{\Delta }SO{C}_{cycle\_4}) $; |
9:if $ 0\le \mathrm{\Delta }SO{C}_{cycle\_2}\le \mathrm{\Delta }SO{C}_{\mathit{max}\_cycle\_2}\cdot {S}_{cycle\_2} $ Indicate that the allocation is reasonable, SOC trajectory planning ending; |
10:else $ \mathrm{\Delta }SO{C}_{cycle\_2}>\mathrm{\Delta }SO{C}_{\mathit{max}\_cycle\_2}\cdot {S}_{cycle\_2} $ Allocate the SOC of the roads of Type-2 typical driving cycle:$ \mathrm{\Delta }SO{C}_{cycle\_2}=\mathrm{\Delta }SO{C}_{\mathit{max}\_cycle\_2}\times {S}_{cycle\_2} $ |
11:Reallocate the SOC of the roads of Type-1 typical driving cycle:$ \mathrm{\Delta }SO{C}_{cycle\_1}=\mathrm{\Delta }SO{C}_{total}-(\mathrm{\Delta }SO{C}_{cycle\_2}+\mathrm{\Delta }SO{C}_{cycle\_3}+\mathrm{\Delta }SO{C}_{cycle\_4}) $ |
12:SOC trajectory planning ending. |
The relevant symbolic meanings and calculation methods are shown in Table 8.
Symbol | Meaning |
$ \mathrm{\Delta }SO{C}_{total} $ | SOC consumption value of the whole driving cycle |
$ \mathrm{\Delta }SO{C}_{cycle\_i}, i=\mathrm{1, 2}, \mathrm{3, 4} $ | SOC distribution value of type i driving cycle |
$ {S}_{cycle\_i}, i=\mathrm{1, 2}, \mathrm{3, 4} $ | Mileage of each type of driving cycle |
$ \mathrm{\Delta }SO{C}_{EV\_cycle\_i}, i=\mathrm{1, 2}, \mathrm{3, 4} $ | SOC distribution value in EV mode of each type of driving cycle |
$ \mathrm{\Delta }SO{C}_{opt\_cycle\_i}, i=\mathrm{1, 2}, \mathrm{3, 4} $ | Optimal SOC distribution value in EV mode of each type of driving cycle |
$ \mathrm{\Delta }SO{C}_{\mathit{max}\_cycle\_i}, i=\mathrm{1, 2}, \mathrm{3, 4} $ | Maximum SOC distribution value of each type of driving cycle |
In order to verify whether the above SOC trajectory planning algorithm is reasonable and effective, this paper collects a section of data through ITS which only represents the information on a certain route. and uses the method in Section 2.2 to convert to the target driving cycle for simulation and verification. And compares the planned SOC trajectory with the SOC trajectory obtained by DP.
As shown in Figure 27, the driving cycle is divided into 9 sections. Then, it is necessary to calculate the mileage of each driving cycle, and according to the above analysis results of SOC consumption of four typical driving cycle segments, the reasonable SOC distribution value of each typical driving cycle segment can be obtained. The SOC of each driving cycle can be allocated according to the mileage and driving cycle category of each driving cycle, and the detailed information of each driving cycle can be obtained, as shown in Table 9.
Segment numbers | 1 | 2 | 3 | 4 | 5 |
Period (s) | 1-363 | 364-1355 | 1356-1727 | 1728-2492 | 2493-2966 |
Distance (km) | 1.7469 | 11.2342 | 1.3181 | 10.5460 | 2.2507 |
Type of driving cycle | 1 | 2 | 3 | 2 | 1 |
$ \mathrm{\Delta }SOC $ (%) | 2.3110 | 9.8721 | 1.4449 | 9.2674 | 2.9775 |
Segment numbers | 6 | 7 | 8 | 9 | |
Period (s) | 2967-3299 | 3300-4095 | 4096-4824 | 4825-5120 | |
Distance (km) | 0.3181 | 11.8117 | 2.7869 | 0.3678 | |
Type of driving cycle | 4 | 2 | 3 | 4 | |
$ \mathrm{\Delta }SOC $ (%) | 0.3212 | 10.3796 | 3.0550 | 0.3713 |
As shown in Figure 28, the reference SOC trajectory is close to the SOC trajectory curve obtained based on DP algorithm, which shows that the SOC reference trajectory obtained by the SOC trajectory planning algorithm is approximately the optimal SOC consumption trajectory.
The driving cycle in Figure 27 is taken as the target driving cycle to carry out simulation comparison of different control strategies, mainly including A-ECMS based on vehicle speed prediction, A-ECMS without vehicle speed prediction, DP and CD-CS. The SOC reference trajectory is shown in the Figure 28 and the simulation results are shown in the Figure 29-32 and Table 10.
Control strategy | DP | A-ECMS based on speed prediction | A-ECMS without speed prediction | CD-CS |
The initial value of SOC | 0.80 | |||
The final value of SOC | 0.40 | 0.4017 | 0.4035 | 0.3956 |
Equivalent fuel consumption (L/100km) | 2.2543 | 2.3601 | 2.4127 | 2.5841 |
Savings of fuel economy | 12.8% | 8.7% | 6.6% | - |
Figure 29 shows vehicle speed prediction in 5s by NAR, the predicted speed and the real speed can be well matched. Figure 30 shows the SOC decline curve under different control strategies. The A-ECMS takes the SOC trajectory obtained from the SOC planning rules mentioned in this paper as the SOC reference trajectory, and modifies the EF in real time. At the same time, combined with the short-term prediction of vehicle speed, when the speed changes greatly in the future, the EF is slightly different from the optimal EF, then the constant part of the adaptive EF is adjusted to the appropriate value. For example, in Figure 31, it can be seen that the EF increases greatly between about 270-500s, which is due to the high-speed driving cycle. The figure also shows that the fluctuation law of the EF is similar to the speed change curve of the simulation driving cycle.
It can also be seen from Figure 30 that the SOC descent trajectory of A-ECMS strategy based on speed prediction is closer to that of DP and can follow the SOC reference trajectory well, while the SOC trajectory of A-ECMS strategy without speed prediction is quite different. It can also be seen from Table 10 that the A-ECMS based on speed prediction and the DP based strategy can better control the SOC terminal value near the target value. For the A-ECMS without speed prediction, the slight difference between the SOC terminal value and the target value is also acceptable. The final value of SOC of CD-CS is more different. It can be seen from the comparison of engine operating points based on the A-ECMS based on speed prediction and CD-CS strategy in Figure 31 that both rotate speeds are concentrated at 1500-3000 rpm, but it is obvious that A-ECMS based on speed prediction engine operating points are more concentrated in the efficient area, which is conducive to the fuel economy of the whole vehicle.
Finally, results in Table 10 show that the fuel saving performance of A-ECMS based on vehicle speed prediction is 8.7% higher than that of CD-CS and same as this is 6.6% in A-ECMS without prediction, which indicates that when the future trip mileage is higher than the pure electric driving mileage under SOC state, the power consumption can be used more reasonably to achieve better fuel saving effect, which is helpful to improve the fuel economy of PHEV.
In order to solve the energy management problems for PHEV, such as the lack of self-adaptability of driving cycle and the inability of real-time online control, the research on the driving cycle prediction and adaptive control strategy of PHEV is carried out.
First of all, this paper proposes an A-ECMS for PHEV energy management, combined with long-term driving cycle recognition by ITS and short-term vehicle speed prediction by NAR neural network to real-time modify EF. Then through the analysis of historical driving cycles by PCA and FCM clustering analysis based on genetic simulated annealing algorithm, this paper obtain four typical driving cycles of the specific driver. Meanwhile, the relationship between SOC consumption, fuel consumption and engine efficiency under typical driving cycle by DP is studied which contributes to establish optimal SOC reference trajectory. Finally, the effect of A-ECMS based on speed prediction proposed in this paper is verified by simulation. Simulation results show that compared with CD-CS and A-ECMS without speed prediction, the fuel consumption of A-ECMS based on speed prediction is reduced by 8.7% and 2.1%. SOC trajectory under A-ECMS basically coincides with SOC trajectory under DP which shows that the algorithm in this paper can apply the optimized results of DP to real-time control and take full advantage of energy saving potential for PHEV.
In this paper, driving data is obtained by the specific driver driving on a substantially fixed route. Then in the further study, it will expand the database and combine Internet of Vehicles and big data analysis to study the versatility of the algorithm and use it for real vehicle control.
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by the National Key R & D Program of China (2018YFB0106200) and Natural Science Foundation of China (51775229).
[1] |
Binnig G, Quate CF, Gerber C (1986) Atomic Force Microscope. Phys Rev Lett 56: 930–933. doi: 10.1103/PhysRevLett.56.930
![]() |
[2] |
Tessmer I, Kaur P, Lin JG, et al. (2013) Investigating bioconjugation by atomic force microscopy. J Nanobiotechnol 11: 1–17. doi: 10.1186/1477-3155-11-1
![]() |
[3] |
Werten PJL, Remigy HW, de Groot BL, et al. (2002) Progress in the analysis of membrane protein structure and function. FEBS Lett 529: 65–72. doi: 10.1016/S0014-5793(02)03290-8
![]() |
[4] |
Muller DJ, Sapra KT, Scheuring S, et al. (2006) Single-molecule studies of membrane proteins. Curr Opin Struct Biol 16: 489–495. doi: 10.1016/j.sbi.2006.06.001
![]() |
[5] |
Gorle S, Pan YG, Sun ZQ, et al. (2017) Computational Model and Dynamics of Monomeric Full-Length APOBEC3G. ACS Cent Sci 3: 1180–1188. doi: 10.1021/acscentsci.7b00346
![]() |
[6] | Sander B, Tria G, Shkumatov AV, et al. (2013) Structural characterization of gephyrin by AFM and SAXS reveals a mixture of compact and extended states. Acta Crystallogr 69: 2050–2060. |
[7] |
Ishino S, Yamagami T, Kitamura M, et al. (2014) Multiple interactions of the intrinsically disordered region between the helicase and nuclease domains of the archaeal Hef protein. J Biol Chem 289: 21627–21639. doi: 10.1074/jbc.M114.554998
![]() |
[8] |
Shinozaki Y, Sumitomo K, Tsuda M, et al. (2009) Direct Observation of ATP-Induced Conformational Changes in Single P2X(4) Receptors. PLos Biol 7: e1000103. doi: 10.1371/journal.pbio.1000103
![]() |
[9] |
Lemaire PA, Tessmer I, Craig R, et al. (2006) Unactivated PKR exists in an open conformation capable of binding nucleotides. Biochemistry 45: 9074–9084. doi: 10.1021/bi060567d
![]() |
[10] |
Kinoshita E, van Rossum-Fikkert S, Sanchez H, et al. (2015) Human RAD50 makes a functional DNA-binding complex. Biochimie 113: 47–53. doi: 10.1016/j.biochi.2015.03.017
![]() |
[11] |
Bonazza K, Rottensteiner H, Seyfried BK, et al. (2014) Visualization of a protein-protein interaction at a single-molecule level by atomic force microscopy. Anal Bioanal Chem 406: 1411–1421. doi: 10.1007/s00216-013-7563-0
![]() |
[12] |
Shlyakhtenko LS, Gall AA, Filonov A, et al. (2003) Silatrane-based surface chemistry for immobilization of DNA, protein-DNA complexes and other biological materials. Ultramicroscopy 97: 279–287. doi: 10.1016/S0304-3991(03)00053-6
![]() |
[13] |
Scheuring S, Muller DJ, Ringler P, et al. (1999) Imaging streptavidin 2D crystals on biotinylated lipid monolayers at high resolution with the atomic force microscope. J Microsc 193: 28–35. doi: 10.1046/j.1365-2818.1999.00434.x
![]() |
[14] |
Whited AM, Park PS (2014) Atomic force microscopy: A multifaceted tool to study membrane proteins and their interactions with ligands. Biochim Biophys Acta 1838: 56–68. doi: 10.1016/j.bbamem.2013.04.011
![]() |
[15] |
Dubrovin EV, Schachtele M, Klinov DV, et al. (2017) Time-Lapse Single-Biomolecule Atomic Force Microscopy Investigation on Modified Graphite in Solution. Langmuir 33: 10027–10034. doi: 10.1021/acs.langmuir.7b02220
![]() |
[16] |
Oliveira Brett AM, Chiorcea Paquim AM (2005) DNA imaged on a HOPG electrode surface by AFM with controlled potential. Bioelectrochemistry 66: 117–124. doi: 10.1016/j.bioelechem.2004.05.009
![]() |
[17] |
Garcia R, Perez R (2002) Dynamic atomic force microscopy methods. Surf Sci Rep 47: 197–301. doi: 10.1016/S0167-5729(02)00077-8
![]() |
[18] |
Dufrene YF, Ando T, Garcia R, et al. (2017) Imaging modes of atomic force microscopy for application in molecular and cell biology. Nat Nanotechnol 12: 295–307. doi: 10.1038/nnano.2017.45
![]() |
[19] |
Hansma HG, Sinsheimer RL, Groppe J, et al. (1993) Recent Advances in Atomic-Force Microscopy of DNA. Scanning 15: 296–299. doi: 10.1002/sca.4950150509
![]() |
[20] |
Balamurugan S, Obubuafo A, Soper SA, et al. (2008) Surface immobilization methods for aptamer diagnostic applications. Anal Bioanal Chem 390: 1009–1021. doi: 10.1007/s00216-007-1587-2
![]() |
[21] |
Knecht S, Ricklin D, Eberle AN, et al. (2009) Oligohis-tags: Mechanisms of binding to Ni2+-NTA surfaces. J Mol Recognit 22: 270–279. doi: 10.1002/jmr.941
![]() |
[22] |
Ritzefeld M, Walhorn V, Anselmetti D, et al. (2013) Analysis of DNA interactions using single-molecule force spectroscopy. Amino Acids 44: 1457–1475. doi: 10.1007/s00726-013-1474-4
![]() |
[23] |
Rief M, Clausen-Schaumann H, Gaub HE (1999) Sequence-dependent mechanics of single DNA molecules. Nat Struct Biol 6: 346–349. doi: 10.1038/7582
![]() |
[24] |
Carrion-Vazquez M, Oberhauser AF, Fowler SB, et al. (1999) Mechanical and chemical unfolding of a single protein: A comparison. Proc Natl Acad Sci U. S. A 96: 3694–3699. doi: 10.1073/pnas.96.7.3694
![]() |
[25] |
Woodside MT, Block SM (2014) Reconstructing Folding Energy Landscapes by Single-Molecule Force Spectroscopy. Annu Rev Biophys 43: 19–39. doi: 10.1146/annurev-biophys-051013-022754
![]() |
[26] |
Hughes ML, Dougan L (2016) The physics of pulling polyproteins: A review of single molecule force spectroscopy using the AFM to study protein unfolding. Rep Prog Phys Phys Soc 79: 076601. doi: 10.1088/0034-4885/79/7/076601
![]() |
[27] |
Fisher TE, Marszalek PE, Fernandez JM (2000) Stretching single molecules into novel conformations using the atomic force microscope. Nat Struct Biol 7: 719–724. doi: 10.1038/78936
![]() |
[28] |
Beckwitt EC, Kong M, Van Houten B (2018) Studying protein-DNA interactions using atomic force microscopy. Semin Cell Dev Biol 73: 220–230. doi: 10.1016/j.semcdb.2017.06.028
![]() |
[29] |
Kasas S, Dietler G (2018) DNA-protein interactions explored by atomic force microscopy. Semin Cell Dev Biol 73: 231–239. doi: 10.1016/j.semcdb.2017.07.015
![]() |
[30] |
Lyubchenko YL, Shlyakhtenko LS (2016) Imaging of DNA and Protein-DNA Complexes with Atomic Force Microscopy. Crit Rev Eukaryot Gene Expr 26: 63–96. doi: 10.1615/CritRevEukaryotGeneExpr.v26.i1.70
![]() |
[31] |
Halford SE, Marko JF (2004) How do site-specific DNA-binding proteins find their targets? Nucleic Acids Res 32: 3040–3052. doi: 10.1093/nar/gkh624
![]() |
[32] |
Schneider R, Grosschedl R (2007) Dynamics and interplay of nuclear architecture, genome organization, and gene expression. Genes Dev 21: 3027–3043. doi: 10.1101/gad.1604607
![]() |
[33] |
Lambert SA, Jolma A, Campitelli LF, et al. (2018) The Human Transcription Factors. Cell 172: 650–665. doi: 10.1016/j.cell.2018.01.029
![]() |
[34] |
Sancar A, Lindsey-Boltz LA, Unsal-Kacmaz K, et al. (2004) Molecular mechanisms of mammalian DNA repair and the DNA damage checkpoints. Annu Rev Biochem 73: 39–85. doi: 10.1146/annurev.biochem.73.011303.073723
![]() |
[35] |
Kapanidis AN, Strick T (2009) Biology, one molecule at a time. Trends Biochem Sci 34: 234–243. doi: 10.1016/j.tibs.2009.01.008
![]() |
[36] |
Larson MH, Landick R, Block SM (2011) Single-Molecule Studies of RNA Polymerase: One Singular Sensation, Every Little Step It Takes. Mol Cell 41: 249–262. doi: 10.1016/j.molcel.2011.01.008
![]() |
[37] |
Ghodke H, Wang H, Hsieh CL, et al. (2014) Single-molecule analysis reveals human UV-damaged DNA-binding protein (UV-DDB) dimerizes on DNA via multiple kinetic intermediates. P Natl Acad Sci USA 111: E1862–E1871. doi: 10.1073/pnas.1323856111
![]() |
[38] |
Buechner CN, Maiti A, Drohat AC, et al. (2015) Lesion search and recognition by thymine DNA glycosylase revealed by single molecule imaging. Nucleic Acids Res 43: 2716–2729. doi: 10.1093/nar/gkv139
![]() |
[39] |
Bewley CA, Gronenborn AM, Clore GM (1998) Minor groove-binding architectural proteins: Structure, function, and DNA recognition. Annu Rev Biophys Biomol Struct 27: 105–131. doi: 10.1146/annurev.biophys.27.1.105
![]() |
[40] |
Perez-Rueda E, Hernandez-Guerrero R, Martinez-Nunez MA, et al. (2018) Abundance, diversity and domain architecture variability in prokaryotic DNA-binding transcription factors. PLos One 13: e0195332. doi: 10.1371/journal.pone.0195332
![]() |
[41] |
Richards FM, Kundrot CE (1988) Identification of Structural Motifs from Protein Coordinate Data-Secondary Structure and 1st-Level Supersecondary Structure. Proteins 3: 71–84. doi: 10.1002/prot.340030202
![]() |
[42] |
Shanahan HP, Garcia MA, Jones S, et al. (2004) Identifying DNA-binding proteins using structural motifs and the electrostatic potential. Nucleic Acids Res 32: 4732–4741. doi: 10.1093/nar/gkh803
![]() |
[43] |
Bustin M, Reeves R (1996) High-mobility-group chromosomal proteins: Architectural components that facilitate chromatin function. Prog Nucleic Acid Res Mol Biol 54: 35–100. doi: 10.1016/S0079-6603(08)60360-8
![]() |
[44] |
Smith NC, Matthews JM (2016) Mechanisms of DNA-binding specificity and functional gene regulation by transcription factors. Curr Opin Struct Biol 38: 68–74. doi: 10.1016/j.sbi.2016.05.006
![]() |
[45] |
Saravanan M, Vasu K, Nagaraja V (2008) Evolution of sequence specificity in a restriction endonuclease by a point mutation. Proc Natl Acad Sci U. S. A 105: 10344–10347. doi: 10.1073/pnas.0804974105
![]() |
[46] |
Ferredamare AR, Prendergast GC, Ziff EB, et al. (1993) Recognition by Max of Its Cognate DNA through a Dimeric B/Hlh/Z Domain. Nature 363: 38–45. doi: 10.1038/363038a0
![]() |
[47] |
Morgunova E, Yin Y, Jolma A, et al. (2015) Structural insights into the DNA-binding specificity of E2F family transcription factors. Nat Commun 6: 10050. doi: 10.1038/ncomms10050
![]() |
[48] |
Li J, Rodriguez JP, Niu F, et al. (2016) Structural basis for DNA recognition by STAT6. P Natl Acad Sci USA 113: 13015–13020. doi: 10.1073/pnas.1611228113
![]() |
[49] |
Rudolph MJ, Gergen JP (2001) DNA-binding by Ig-fold proteins. Nat Struct Biol 8: 384–386. doi: 10.1038/87531
![]() |
[50] |
Doherty AJ, Serpell LC, Ponting CP (1996) The helix-hairpin-helix DNA-binding motif: A structural basis for non-sequence-specific recognition of DNA. Nucleic Acids Res 24: 2488–2497. doi: 10.1093/nar/24.13.2488
![]() |
[51] |
Burgers PMJ, Kunkel TA (2017) Eukaryotic DNA Replication Fork. Annu Rev Biochem 86: 417–438. doi: 10.1146/annurev-biochem-061516-044709
![]() |
[52] |
Raghunathan S, Kozlov AG, Lohman TM, et al. (2000) Structure of the DNA binding domain of E-coli SSB bound to ssDNA. Nat Struct Biol 7: 648–652. doi: 10.1038/77943
![]() |
[53] |
Theis K, Chen PJ, Skorvaga M, et al. (1999) Crystal structure of UvrB, a DNA helicase adapted for nucleotide excision repair. EMBO J 18: 6899–6907. doi: 10.1093/emboj/18.24.6899
![]() |
[54] |
Waters TR, Eryilmaz J, Geddes S, et al. (2006) Damage detection by the UvrABC pathway: Crystal structure of UvrB bound to fluorescein-adducted DNA. FEBS Lett 580: 6423–6427. doi: 10.1016/j.febslet.2006.10.051
![]() |
[55] |
Chai N, Li WX, Wang J, et al. (2015) Structural basis for the Smad5 MH1 domain to recognize different DNA sequences. Nucleic Acids Res 43: 9051–9064. doi: 10.1093/nar/gkv848
![]() |
[56] |
Maiti A, Morgan MT, Pozharski E, et al. (2008) Crystal structure of human thymine DNA glycosylase bound to DNA elucidates sequence-specific mismatch recognition. P Natl Acad Sci USA 105: 8890–8895. doi: 10.1073/pnas.0711061105
![]() |
[57] |
Bruner SD, Norman DPG, Verdine GL (2000) Structural basis for recognition and repair of the endogenous mutagen 8-oxoguanine in DNA. Nature 403: 859–866. doi: 10.1038/35002510
![]() |
[58] |
Ando T, Kodera N, Takai E, et al. (2001) A high-speed atomic force microscope for studying biological macromolecules. P Natl Acad Sci USA 98: 12468–12472. doi: 10.1073/pnas.211400898
![]() |
[59] |
Ando T (2018) High-speed atomic force microscopy and its future prospects. Biophys Rev 10: 285–292. doi: 10.1007/s12551-017-0356-5
![]() |
[60] |
Sawicka M, Aramayo R, Ayala R, et al. (2017) Single-Particle Electron Microscopy Analysis of DNA Repair Complexes. Methods Enzymol 592: 159–186. doi: 10.1016/bs.mie.2017.03.010
![]() |
[61] |
Sun JC, Yuan ZN, Bai L, et al. (2017) Cryo-EM of dynamic protein complexes in eukaryotic DNA replication. Protein Sci 26: 40–51. doi: 10.1002/pro.3033
![]() |
[62] |
Tessmer I, Moore T, Lloyd RG, et al. (2005) AFM studies on the role of the protein RdgC in bacterial DNA recombination. J Mol Biol 350: 254–262. doi: 10.1016/j.jmb.2005.04.043
![]() |
[63] |
Lohr D, Bash R, Wang H, et al. (2007) Using atomic force microscopy to study chromatin structure and nucleosome remodeling. Methods 41: 333–341. doi: 10.1016/j.ymeth.2006.08.016
![]() |
[64] |
Hizume K, Kominam H, Kobayash K, et al. (2017) Flexible DNA Path in the MCM Double Hexamer Loaded on DNA. Biochemistry 56: 2435–2445. doi: 10.1021/acs.biochem.6b00922
![]() |
[65] |
Noort JV, Heijden TVD, Dutta CF, et al. (2004) Initiation of translocation by Type I restriction-modification enzymes is associated with a short DNA extrusion. Nucleic Acids Res 32: 6540–6547. doi: 10.1093/nar/gkh999
![]() |
[66] |
Maurer S, Fritz J, Muskhelishvili G, et al. (2006) RNA polymerase and an activator form discrete subcomplexes in a transcription initiation complex. EMBO J 25: 3784–3790. doi: 10.1038/sj.emboj.7601261
![]() |
[67] |
Verhoeven EEA, Wyman C, Moolenaar GF, et al. (2002) The presence of two UvrB subunits in the UvrAB complex ensures damage detection in both DNA strands. EMBO J 21: 4196–4205. doi: 10.1093/emboj/cdf396
![]() |
[68] |
Cellai S, Mangiarotti L, Vannini N, et al. (2007) Upstream promoter sequences and alpha CTD mediate stable DNA wrapping within the RNA polymerase-promoter open complex. EMBO Rep 8: 271–278. doi: 10.1038/sj.embor.7400888
![]() |
[69] | Umemura K, Okada T, Kuroda R (2005) Cooperativity and intermediate structures of single-stranded DNA binding-assisted RecA-single-stranded DNA complex formation studied by atomic force microscopy. Scanning 27: 35–43. |
[70] |
Hamon L, Pastre D, Dupaigne P, et al. (2007) High-resolution AFM imaging of single-stranded DNA-binding (SSB) protein-DNA complexes. Nucleic Acids Res 35: e58. doi: 10.1093/nar/gkm147
![]() |
[71] |
Li BS, Goh MC (2010) Direct visualization of the formation and structure of RecA/dsDNA complexes. Micron 41: 227–231. doi: 10.1016/j.micron.2009.10.011
![]() |
[72] |
Li BS, Wei B, Goh MC (2012) Direct visualization of the formation of RecA/dsDNA complexes at the single-molecule level. Micron 43: 1073–1075. doi: 10.1016/j.micron.2012.04.016
![]() |
[73] |
Tessmer I, Melikishvili M, Fried MG (2012) Cooperative cluster formation, DNA bending and base-flipping by O6-alkylguanine-DNA alkyltransferase. Nucleic Acids Res 40: 8296–8308. doi: 10.1093/nar/gks574
![]() |
[74] |
Miyagi A, Ando T, Lyubchenko YL (2011) Dynamics of Nucleosomes Assessed with Time-Lapse High-Speed Atomic Force Microscopy. Biochemistry 50: 7901–7908. doi: 10.1021/bi200946z
![]() |
[75] |
Sanchez H, Suzuki Y, Yokokawa M, et al. (2011) Protein-DNA interactions in high speed AFM: Single molecule diffusion analysis of human RAD54. Integr Biol 3: 1127–1134. doi: 10.1039/c1ib00039j
![]() |
[76] |
Endo M, Sugiyama H (2014) Single-Molecule Imaging of Dynamic Motions of Biomolecules in DNA Origami Nanostructures Using High-Speed Atomic Force Microscopy. Acc Chem Res 47: 1645–1653. doi: 10.1021/ar400299m
![]() |
[77] |
Lee AJ, Endo M, Hobbs JK, et al. (2018) Direct Single-Molecule Observation of Mode and Geometry of RecA-Mediated Homology Search. Acs Nano 12: 272–278. doi: 10.1021/acsnano.7b06208
![]() |
[78] |
Suzuki Y, Shin M, Yoshida A, et al. (2012) Fast microscopical dissection of action scenes played by Escherichia coli RNA polymerase. FEBS Lett 586: 3187–3192. doi: 10.1016/j.febslet.2012.06.033
![]() |
[79] |
Buechner CN, Tessmer I (2013) DNA substrate preparation for atomic force microscopy studies of protein-DNA interactions. J Mol Recognit 26: 605–617. doi: 10.1002/jmr.2311
![]() |
[80] |
Yang Y, Sass LE, Du C, et al. (2005) Determination of protein-DNA binding constants and specificities from statistical analyses of single molecules: MutS-DNA interactions. Nucleic Acids Res 33: 4322–4334. doi: 10.1093/nar/gki708
![]() |
[81] |
Sukhanova MV, Abrakhi S, Joshi V, et al. (2016) Single molecule detection of PARP1 and PARP2 interaction with DNA strand breaks and their poly(ADP-ribosyl)ation using high-resolution AFM imaging. Nucleic Acids Res 44: e60. doi: 10.1093/nar/gkv1476
![]() |
[82] |
Buechner CN, Heil K, Michels G, et al. (2014) Strand-specific Recognition of DNA Damages by XPD Provides Insights into Nucleotide Excision Repair Substrate Versatility. J Biol Chem 289: 3613–3624. doi: 10.1074/jbc.M113.523001
![]() |
[83] |
Wirth N, Gross J, Roth HM, et al. (2016) Conservation and Divergence in Nucleotide Excision Repair Lesion Recognition. J Biol Chem 291: 18932–18946. doi: 10.1074/jbc.M116.739425
![]() |
[84] |
Doniselli N, Rodriguez-Aliaga P, Amidani D, et al. (2015) New insights into the regulatory mechanisms of ppGpp and DksA on Escherichia coli RNA polymerase-promoter complex. Nucleic Acids Res 43: 5249–5262. doi: 10.1093/nar/gkv391
![]() |
[85] |
Nettikadan S, Tokumasu F, Takeyasu K (1996) Quantitative analysis of the transcription factor AP2 binding to DNA by atomic force microscopy. Biochem Biophys Res Commun 226: 645–649. doi: 10.1006/bbrc.1996.1409
![]() |
[86] |
Timofeeva OA, Chasovskikh S, Lonskaya I, et al. (2012) Mechanisms of Unphosphorylated STAT3 Transcription Factor Binding to DNA. J Biol Chem 287: 14192–14200. doi: 10.1074/jbc.M111.323899
![]() |
[87] |
Zhang P, Xia JH, Zhu J, et al. (2018) High-throughput screening of prostate cancer risk loci by single nucleotide polymorphisms sequencing. Nat Commun 9: 2022. doi: 10.1038/s41467-018-04451-x
![]() |
[88] |
Huang Q, Whitington T, Gao P, et al. (2014) A prostate cancer susceptibility allele at 6q22 increases RFX6 expression by modulating HOXB13 chromatin binding. Nat Genet 46: 126–135. doi: 10.1038/ng.2862
![]() |
[89] |
Gao P, Xia JH, Sipeky C, et al. (2018) Biology and Clinical Implications of the 19q13 Aggressive Prostate Cancer Susceptibility Locus. Cell 174: 576–589. doi: 10.1016/j.cell.2018.06.003
![]() |
[90] |
Crampton N, Bonass WA, Kirkham J, et al. (2006) Collision events between RNA polymerases in convergent transcription studied by atomic force microscopy. Nucleic Acids Res 34: 5416–5425. doi: 10.1093/nar/gkl668
![]() |
[91] |
Countryman P, Fan Y, Gorthi A, et al. (2018) Cohesin SA2 is a sequence-independent DNA-binding protein that recognizes DNA replication and repair intermediates. J Biol Chem 293: 1054–1069. doi: 10.1074/jbc.M117.806406
![]() |
[92] |
Schneider SW, Larmer J, Henderson RM, et al. (1998) Molecular weights of individual proteins correlate with molecular volumes measured by atomic force microscopy. Pflug Arch Eur J Phy 435: 362–367. doi: 10.1007/s004240050524
![]() |
[93] |
Ratcliff GC, Erie DA (2001) A novel single-molecule study to determine protein-protein association constants. J Am Chem Soc 123: 5632–5635. doi: 10.1021/ja005750n
![]() |
[94] |
Wang H, Dellavecchia MJ, Skorvaga M, et al. (2006) UvrB domain 4, an autoinhibitory gate for regulation of DNA binding and ATPase activity. J Biol Chem 281: 15227–15237. doi: 10.1074/jbc.M601476200
![]() |
[95] |
Roth HM, Romer J, Grundler V, et al. (2012) XPB helicase regulates DNA incision by the Thermoplasma acidophilum endonuclease Bax1. DNA Repair 11: 286–293. doi: 10.1016/j.dnarep.2011.12.002
![]() |
[96] |
Fuentes-Perez ME, Dillingham MS, Moreno-Herrero F (2013) AFM volumetric methods for the characterization of proteins and nucleic acids. Methods 60: 113–121. doi: 10.1016/j.ymeth.2013.02.005
![]() |
[97] | Amidani D, Tramonti A, Canosa AV, et al. (2016) Study of DNA binding and bending by Bacillus subtilis GabR, a PLP-dependent transcription factor. Biochim Biophys Acta Gen Subj 1861: 3474–3489. |
[98] |
Rivetti C, Guthold M, Bustamante C (1996) Scanning force microscopy of DNA deposited onto mica: Equilibration versus kinetic trapping studied by statistical polymer chain analysis. J Mol Biol 264: 919–932. doi: 10.1006/jmbi.1996.0687
![]() |
[99] |
Cassina V, Manghi M, Salerno D, et al. (2016) Effects of cytosine methylation on DNA morphology: An atomic force microscopy study. Biochim Biophys Acta 1860: 1–7. doi: 10.1016/j.bbagen.2015.10.006
![]() |
[100] |
Scipioni A, Anselmi C, Zuccheri G, et al. (2002) Sequence-dependent DNA curvature and flexibility from scanning force microscopy images. Biophys J 83: 2408–2418. doi: 10.1016/S0006-3495(02)75254-5
![]() |
[101] |
Moukhtar J, Faivre-Moskalenko C, Milani P, et al. (2010) Effect of genomic long-range correlations on DNA persistence length: from theory to single molecule experiments. J Phys Chem B 114: 5125–5143. doi: 10.1021/jp911031y
![]() |
[102] |
Jager MD, Noort JV, Gent DCV, et al. (2001) Human Rad50/Mre11 is a flexible complex that can tether DNA ends. Mol Cell 8: 1129–1135. doi: 10.1016/S1097-2765(01)00381-1
![]() |
[103] |
Tessmer I, Yang Y, Zhai J, et al. (2008) Mechanism of MutS searching for DNA mismatches and signaling repair. J Biol Chem 283: 36646–36654. doi: 10.1074/jbc.M805712200
![]() |
[104] |
Bosch D, Campillo M, Pardo L (2003) Binding of proteins to the minor groove of DNA: What are the structural and energetic determinants for kinking a basepair step? J Comput Chem 24: 682–691. doi: 10.1002/jcc.10200
![]() |
[105] |
Kong MW, Liu LL, Chen XJ, et al. (2016) Single-Molecule Imaging Reveals that Rad4 Employs a Dynamic DNA Damage Recognition Process. Mol Cell 64: 376–387. doi: 10.1016/j.molcel.2016.09.005
![]() |
[106] |
Wang H, Yang Y, Schofield MJ, et al. (2003) DNA bending and unbending by MutS govern mismatch recognition and specificity. Proc Natl Acad Sci U. S. A 100: 14822–14827. doi: 10.1073/pnas.2433654100
![]() |
[107] |
Chen L, Haushalter KA, Lieber CM, et al. (2002) Direct visualization of a DNA glycosylase searching for damage. Chem Biol 9: 345–350. doi: 10.1016/S1074-5521(02)00120-5
![]() |
[108] |
Lamers MH, Perrakis A, Enzlin JH, et al. (2000) The crystal structure of DNA mismatch repair protein MutS binding to a G center dot T mismatch. Nature 407: 711–717. doi: 10.1038/35037523
![]() |
[109] |
Koroleva ON, Dubrovin EV, Yaminsky IV, et al. (2016) Effect of DNA bending on transcriptional interference in the systems of closely spaced convergent promoters. Biochim Biophys Acta 1860: 2086–2096. doi: 10.1016/j.bbagen.2016.06.026
![]() |
[110] |
Fronczek DN, Quammen C, Wang H, et al. (2011) High accuracy FIONA-AFM hybrid imaging. Ultramicroscopy 111: 350–355. doi: 10.1016/j.ultramic.2011.01.020
![]() |
[111] |
Sanchez H, Kertokalio A, van Rossum-Fikkert S, et al. (2013) Combined optical and topographic imaging reveals different arrangements of human RAD54 with presynaptic and postsynaptic RAD51-DNA filaments. P Natl Acad Sci USA 110: 11385–11390. doi: 10.1073/pnas.1306467110
![]() |
[112] |
Frederickx W, Rocha S, Fujita Y, et al. (2018) Orthogonal Probing of Single-Molecule Heterogeneity by Correlative Fluorescence and Force Microscopy. Acs Nano 12: 168–177. doi: 10.1021/acsnano.7b05405
![]() |
[113] |
Schmucker SW, Kumar N, Abelson JR, et al. (2012) Field-directed sputter sharpening for tailored probe materials and atomic-scale lithography. Nat Commun 3: 935. doi: 10.1038/ncomms1907
![]() |
[114] |
Pfreundschuh M, Alsteens D, Hilbert M, et al. (2014) Localizing Chemical Groups while Imaging Single Native Proteins by High-Resolution Atomic Force Microscopy. Nano Lett 14: 2957–2964. doi: 10.1021/nl5012905
![]() |
[115] |
Monig H, Hermoso DR, Arado OD, et al. (2016) Submolecular Imaging by Noncontact Atomic Force Microscopy with an Oxygen Atom Rigidly Connected to a Metallic Probe. Acs Nano 10: 1201–1209. doi: 10.1021/acsnano.5b06513
![]() |
[116] |
Senapati S, Lindsay S (2016) Recent Progress in Molecular Recognition Imaging Using Atomic Force Microscopy. Acc Chem Res 49: 503–510. doi: 10.1021/acs.accounts.5b00533
![]() |
[117] | Khan Z, Leung C, Tahir BA, et al. (2010) Digitally tunable, wide-band amplitude, phase, and frequency detection for atomic-resolution scanning force microscopy. Rev Sci Instrum 81: 197. |
[118] |
Calo A, Eleta-Lopez A, Stoliar P, et al. (2016) Multifrequency Force Microscopy of Helical Protein Assembly on a Virus. Sci Rep 6: 21899. doi: 10.1038/srep21899
![]() |
[119] |
Wu D, Kaur P, Li ZM, et al. (2016) Visualizing the Path of DNA through Proteins Using DREEM Imaging. Mol Cell 61: 315–323. doi: 10.1016/j.molcel.2015.12.012
![]() |
1. | Nathalie Krell, Émeline Schmisser, Nonparametric estimation of jump rates for a specific class of piecewise deterministic Markov processes, 2021, 27, 1350-7265, 10.3150/20-BEJ1312 | |
2. | Alessandro Calvia, Stochastic filtering and optimal control of pure jump Markov processes with noise-free partial observation, 2020, 26, 1292-8119, 25, 10.1051/cocv/2019020 | |
3. | Elena Bandini, Michèle Thieullen, Optimal Control of Infinite-Dimensional Piecewise Deterministic Markov Processes: A BSDE Approach. Application to the Control of an Excitable Cell Membrane, 2021, 84, 0095-4616, 1549, 10.1007/s00245-020-09687-y | |
4. | Benoîte de Saporta, Aymar Thierry d’Argenlieu, Régis Sabbadin, Alice Cleynen, Inés P. Mariño, A Monte-Carlo planning strategy for medical follow-up optimization: Illustration on multiple myeloma data, 2024, 19, 1932-6203, e0315661, 10.1371/journal.pone.0315661 |
Component | Parameter value |
Engine | Peak power: 64kW; Maximum torque :153Nm |
Motor | Peak power: 32Kw; Maximum torque 160Nm |
Battery capacity | 27.9Ah |
Rated voltage of battery | 320V |
Final drive ratio | 3.944 |
Transmission speed ratio | 3.538/2.045/1.333/1.028/0.82 |
Vehicle | Mass: $ m=1568 $ kg; Vehicle equivalent cross section: $ A=2.13{m}^{2} $; Aerodynamic drag coefficient: $ {C}_{d}=0.29 $; Tire radius: $ r=312mm $ |
Prediction durations | 5s | 10s | 15s | 20s |
RMSE | 2.6456 | 5.5578 | 7.8314 | 10.0521 |
Numbers | Driving cycle characteristic parameter | Unit |
1 | Average speed:$ {v}_{ave} $ | km/h |
2 | Average driving speed:$ {u}_{ave} $ | km/h |
3 | Standard deviation of vehicle speed:$ {v}_{std} $ | km/h |
4 | Maximum speed:$ {v}_{\mathit{max} } $ | km/h |
5 | Average acceleration:$ {a}_{ave} $ | m/s2 |
6 | Average deceleration:$ {a}_{dec} $ | m/s2 |
7 | Maximum acceleration:$ {a}_{\mathit{max} } $ | m/s2 |
8 | Maximum deceleration:$ {a}_{\mathit{min} } $ | m/s2 |
9 | Standard deviation of acceleration:$ {a}_{std} $ | m/s2 |
10 | Idle time ratio:$ P\left({t}_{idle}\right) $ | % |
11 | Acceleration time ratio:$ P\left({t}_{acc}\right) $ | % |
12 | Deceleration time ratio:$ P\left({t}_{dec}\right) $ | % |
13 | Uniform time ratio:$ P\left({t}_{uni}\right) $ | % |
Index | F1 | F2 | F3 | F4 |
1 | -1.1058 | -1.3408 | -0.2677 | 1.9947 |
2 | 0.8687 | -1.7607 | -0.2811 | 0.8793 |
3 | -0.1886 | -1.6428 | 0.5574 | 0.7264 |
4 | -0.7973 | -1.8836 | 2.0468 | 0.0025 |
⋮ | ⋮ | ⋮ | ⋮ | ⋮ |
1467 | 3.6596 | -0.3058 | 0.8410 | -1.4611 |
1468 | 3.0708 | -0.3280 | -0.5293 | -0.1305 |
Driving cycle type | 1 | 2 | 3 | 4 |
SOC distribution value per km(%/km) | 1.3229 | 1.3643 | 1.2180 | 1.0602 |
Driving cycle type | 1 | 2 | 3 | 4 |
SOC max distribution value per km(%/km) | 1.1759 | 1.2993 | 0.9744 | 1.0602 |
Optimal SOC max distribution value per km(%/km) | 0.7349 | 0.3898 | 0.8526 | 1.0602 |
SOC min distribution value per km(%/km) | 0 | 0 | 0 | 1.0602 |
1: Identify the road section belonging to Type-4 typical driving cycle; |
2: Allocate the SOC of the roads of Type-4 typical driving cycle:$ \mathrm{\Delta }SO{C}_{cycle\_4}=\mathrm{\Delta }SO{C}_{EV\_cycle\_4}\times {S}_{cycle\_4} $; |
3: Identify the road section belonging to Type 3 of typical driving cycles; |
4:Allocate the SOC of the roads of Type-3 typical driving cycle:$ \mathrm{\Delta }SO{C}_{cycle\_3}=\mathrm{\Delta }SO{C}_{opt\_cycle\_3}\times {S}_{cycle\_3} $; |
5:Identify the road section belonging to Type-1 of typical driving cycles; |
6:Allocate the SOC of the roads of Type-1 typical driving cycle:$ \mathrm{\Delta }SO{C}_{cycle\_1}=\mathrm{\Delta }SO{C}_{opt\_cycle\_1}\times {S}_{cycle\_1} $; |
7:Identify the road section belonging to Type 2 of typical driving cycle; |
8:Allocate the SOC of the roads of Type-2 typical driving cycle:$ \mathrm{\Delta }SO{C}_{cycle\_2}=\mathrm{\Delta }SO{C}_{total}-(\mathrm{\Delta }SO{C}_{cycle\_1}+\mathrm{\Delta }SO{C}_{cycle\_3}+\mathrm{\Delta }SO{C}_{cycle\_4}) $; |
9:if $ 0\le \mathrm{\Delta }SO{C}_{cycle\_2}\le \mathrm{\Delta }SO{C}_{\mathit{max}\_cycle\_2}\cdot {S}_{cycle\_2} $ Indicate that the allocation is reasonable, SOC trajectory planning ending; |
10:else $ \mathrm{\Delta }SO{C}_{cycle\_2}>\mathrm{\Delta }SO{C}_{\mathit{max}\_cycle\_2}\cdot {S}_{cycle\_2} $ Allocate the SOC of the roads of Type-2 typical driving cycle:$ \mathrm{\Delta }SO{C}_{cycle\_2}=\mathrm{\Delta }SO{C}_{\mathit{max}\_cycle\_2}\times {S}_{cycle\_2} $ |
11:Reallocate the SOC of the roads of Type-1 typical driving cycle:$ \mathrm{\Delta }SO{C}_{cycle\_1}=\mathrm{\Delta }SO{C}_{total}-(\mathrm{\Delta }SO{C}_{cycle\_2}+\mathrm{\Delta }SO{C}_{cycle\_3}+\mathrm{\Delta }SO{C}_{cycle\_4}) $ |
12:SOC trajectory planning ending. |
Symbol | Meaning |
$ \mathrm{\Delta }SO{C}_{total} $ | SOC consumption value of the whole driving cycle |
$ \mathrm{\Delta }SO{C}_{cycle\_i}, i=\mathrm{1, 2}, \mathrm{3, 4} $ | SOC distribution value of type i driving cycle |
$ {S}_{cycle\_i}, i=\mathrm{1, 2}, \mathrm{3, 4} $ | Mileage of each type of driving cycle |
$ \mathrm{\Delta }SO{C}_{EV\_cycle\_i}, i=\mathrm{1, 2}, \mathrm{3, 4} $ | SOC distribution value in EV mode of each type of driving cycle |
$ \mathrm{\Delta }SO{C}_{opt\_cycle\_i}, i=\mathrm{1, 2}, \mathrm{3, 4} $ | Optimal SOC distribution value in EV mode of each type of driving cycle |
$ \mathrm{\Delta }SO{C}_{\mathit{max}\_cycle\_i}, i=\mathrm{1, 2}, \mathrm{3, 4} $ | Maximum SOC distribution value of each type of driving cycle |
Segment numbers | 1 | 2 | 3 | 4 | 5 |
Period (s) | 1-363 | 364-1355 | 1356-1727 | 1728-2492 | 2493-2966 |
Distance (km) | 1.7469 | 11.2342 | 1.3181 | 10.5460 | 2.2507 |
Type of driving cycle | 1 | 2 | 3 | 2 | 1 |
$ \mathrm{\Delta }SOC $ (%) | 2.3110 | 9.8721 | 1.4449 | 9.2674 | 2.9775 |
Segment numbers | 6 | 7 | 8 | 9 | |
Period (s) | 2967-3299 | 3300-4095 | 4096-4824 | 4825-5120 | |
Distance (km) | 0.3181 | 11.8117 | 2.7869 | 0.3678 | |
Type of driving cycle | 4 | 2 | 3 | 4 | |
$ \mathrm{\Delta }SOC $ (%) | 0.3212 | 10.3796 | 3.0550 | 0.3713 |
Control strategy | DP | A-ECMS based on speed prediction | A-ECMS without speed prediction | CD-CS |
The initial value of SOC | 0.80 | |||
The final value of SOC | 0.40 | 0.4017 | 0.4035 | 0.3956 |
Equivalent fuel consumption (L/100km) | 2.2543 | 2.3601 | 2.4127 | 2.5841 |
Savings of fuel economy | 12.8% | 8.7% | 6.6% | - |
Component | Parameter value |
Engine | Peak power: 64kW; Maximum torque :153Nm |
Motor | Peak power: 32Kw; Maximum torque 160Nm |
Battery capacity | 27.9Ah |
Rated voltage of battery | 320V |
Final drive ratio | 3.944 |
Transmission speed ratio | 3.538/2.045/1.333/1.028/0.82 |
Vehicle | Mass: $ m=1568 $ kg; Vehicle equivalent cross section: $ A=2.13{m}^{2} $; Aerodynamic drag coefficient: $ {C}_{d}=0.29 $; Tire radius: $ r=312mm $ |
Prediction durations | 5s | 10s | 15s | 20s |
RMSE | 2.6456 | 5.5578 | 7.8314 | 10.0521 |
Numbers | Driving cycle characteristic parameter | Unit |
1 | Average speed:$ {v}_{ave} $ | km/h |
2 | Average driving speed:$ {u}_{ave} $ | km/h |
3 | Standard deviation of vehicle speed:$ {v}_{std} $ | km/h |
4 | Maximum speed:$ {v}_{\mathit{max} } $ | km/h |
5 | Average acceleration:$ {a}_{ave} $ | m/s2 |
6 | Average deceleration:$ {a}_{dec} $ | m/s2 |
7 | Maximum acceleration:$ {a}_{\mathit{max} } $ | m/s2 |
8 | Maximum deceleration:$ {a}_{\mathit{min} } $ | m/s2 |
9 | Standard deviation of acceleration:$ {a}_{std} $ | m/s2 |
10 | Idle time ratio:$ P\left({t}_{idle}\right) $ | % |
11 | Acceleration time ratio:$ P\left({t}_{acc}\right) $ | % |
12 | Deceleration time ratio:$ P\left({t}_{dec}\right) $ | % |
13 | Uniform time ratio:$ P\left({t}_{uni}\right) $ | % |
Index | F1 | F2 | F3 | F4 |
1 | -1.1058 | -1.3408 | -0.2677 | 1.9947 |
2 | 0.8687 | -1.7607 | -0.2811 | 0.8793 |
3 | -0.1886 | -1.6428 | 0.5574 | 0.7264 |
4 | -0.7973 | -1.8836 | 2.0468 | 0.0025 |
⋮ | ⋮ | ⋮ | ⋮ | ⋮ |
1467 | 3.6596 | -0.3058 | 0.8410 | -1.4611 |
1468 | 3.0708 | -0.3280 | -0.5293 | -0.1305 |
Driving cycle type | 1 | 2 | 3 | 4 |
SOC distribution value per km(%/km) | 1.3229 | 1.3643 | 1.2180 | 1.0602 |
Driving cycle type | 1 | 2 | 3 | 4 |
SOC max distribution value per km(%/km) | 1.1759 | 1.2993 | 0.9744 | 1.0602 |
Optimal SOC max distribution value per km(%/km) | 0.7349 | 0.3898 | 0.8526 | 1.0602 |
SOC min distribution value per km(%/km) | 0 | 0 | 0 | 1.0602 |
1: Identify the road section belonging to Type-4 typical driving cycle; |
2: Allocate the SOC of the roads of Type-4 typical driving cycle:$ \mathrm{\Delta }SO{C}_{cycle\_4}=\mathrm{\Delta }SO{C}_{EV\_cycle\_4}\times {S}_{cycle\_4} $; |
3: Identify the road section belonging to Type 3 of typical driving cycles; |
4:Allocate the SOC of the roads of Type-3 typical driving cycle:$ \mathrm{\Delta }SO{C}_{cycle\_3}=\mathrm{\Delta }SO{C}_{opt\_cycle\_3}\times {S}_{cycle\_3} $; |
5:Identify the road section belonging to Type-1 of typical driving cycles; |
6:Allocate the SOC of the roads of Type-1 typical driving cycle:$ \mathrm{\Delta }SO{C}_{cycle\_1}=\mathrm{\Delta }SO{C}_{opt\_cycle\_1}\times {S}_{cycle\_1} $; |
7:Identify the road section belonging to Type 2 of typical driving cycle; |
8:Allocate the SOC of the roads of Type-2 typical driving cycle:$ \mathrm{\Delta }SO{C}_{cycle\_2}=\mathrm{\Delta }SO{C}_{total}-(\mathrm{\Delta }SO{C}_{cycle\_1}+\mathrm{\Delta }SO{C}_{cycle\_3}+\mathrm{\Delta }SO{C}_{cycle\_4}) $; |
9:if $ 0\le \mathrm{\Delta }SO{C}_{cycle\_2}\le \mathrm{\Delta }SO{C}_{\mathit{max}\_cycle\_2}\cdot {S}_{cycle\_2} $ Indicate that the allocation is reasonable, SOC trajectory planning ending; |
10:else $ \mathrm{\Delta }SO{C}_{cycle\_2}>\mathrm{\Delta }SO{C}_{\mathit{max}\_cycle\_2}\cdot {S}_{cycle\_2} $ Allocate the SOC of the roads of Type-2 typical driving cycle:$ \mathrm{\Delta }SO{C}_{cycle\_2}=\mathrm{\Delta }SO{C}_{\mathit{max}\_cycle\_2}\times {S}_{cycle\_2} $ |
11:Reallocate the SOC of the roads of Type-1 typical driving cycle:$ \mathrm{\Delta }SO{C}_{cycle\_1}=\mathrm{\Delta }SO{C}_{total}-(\mathrm{\Delta }SO{C}_{cycle\_2}+\mathrm{\Delta }SO{C}_{cycle\_3}+\mathrm{\Delta }SO{C}_{cycle\_4}) $ |
12:SOC trajectory planning ending. |
Symbol | Meaning |
$ \mathrm{\Delta }SO{C}_{total} $ | SOC consumption value of the whole driving cycle |
$ \mathrm{\Delta }SO{C}_{cycle\_i}, i=\mathrm{1, 2}, \mathrm{3, 4} $ | SOC distribution value of type i driving cycle |
$ {S}_{cycle\_i}, i=\mathrm{1, 2}, \mathrm{3, 4} $ | Mileage of each type of driving cycle |
$ \mathrm{\Delta }SO{C}_{EV\_cycle\_i}, i=\mathrm{1, 2}, \mathrm{3, 4} $ | SOC distribution value in EV mode of each type of driving cycle |
$ \mathrm{\Delta }SO{C}_{opt\_cycle\_i}, i=\mathrm{1, 2}, \mathrm{3, 4} $ | Optimal SOC distribution value in EV mode of each type of driving cycle |
$ \mathrm{\Delta }SO{C}_{\mathit{max}\_cycle\_i}, i=\mathrm{1, 2}, \mathrm{3, 4} $ | Maximum SOC distribution value of each type of driving cycle |
Segment numbers | 1 | 2 | 3 | 4 | 5 |
Period (s) | 1-363 | 364-1355 | 1356-1727 | 1728-2492 | 2493-2966 |
Distance (km) | 1.7469 | 11.2342 | 1.3181 | 10.5460 | 2.2507 |
Type of driving cycle | 1 | 2 | 3 | 2 | 1 |
$ \mathrm{\Delta }SOC $ (%) | 2.3110 | 9.8721 | 1.4449 | 9.2674 | 2.9775 |
Segment numbers | 6 | 7 | 8 | 9 | |
Period (s) | 2967-3299 | 3300-4095 | 4096-4824 | 4825-5120 | |
Distance (km) | 0.3181 | 11.8117 | 2.7869 | 0.3678 | |
Type of driving cycle | 4 | 2 | 3 | 4 | |
$ \mathrm{\Delta }SOC $ (%) | 0.3212 | 10.3796 | 3.0550 | 0.3713 |
Control strategy | DP | A-ECMS based on speed prediction | A-ECMS without speed prediction | CD-CS |
The initial value of SOC | 0.80 | |||
The final value of SOC | 0.40 | 0.4017 | 0.4035 | 0.3956 |
Equivalent fuel consumption (L/100km) | 2.2543 | 2.3601 | 2.4127 | 2.5841 |
Savings of fuel economy | 12.8% | 8.7% | 6.6% | - |