Advanced Control and Optimization with Applications to Complex Automotive SystemsView this Special Issue
Research Article | Open Access
Hongwen He, Henglu Tang, Ximing Wang, "Global Optimal Energy Management Strategy Research for a Plug-In Series-Parallel Hybrid Electric Bus by Using Dynamic Programming", Mathematical Problems in Engineering, vol. 2013, Article ID 708261, 11 pages, 2013. https://doi.org/10.1155/2013/708261
Global Optimal Energy Management Strategy Research for a Plug-In Series-Parallel Hybrid Electric Bus by Using Dynamic Programming
Energy management strategy influences the power performance and fuel economy of plug-in hybrid electric vehicles greatly. To explore the fuel-saving potential of a plug-in hybrid electric bus (PHEB), this paper searched the global optimal energy management strategy using dynamic programming (DP) algorithm. Firstly, the simplified backward model of the PHEB was built which is necessary for DP algorithm. Then the torque and speed of engine and the torque of motor were selected as the control variables, and the battery state of charge (SOC) was selected as the state variables. The DP solution procedure was listed, and the way was presented to find all possible control variables at every state of each stage in detail. Finally, the appropriate SOC increment is determined after quantizing the state variables, and then the optimal control of long driving distance of a specific driving cycle is replaced with the optimal control of one driving cycle, which reduces the computational time significantly and keeps the precision at the same time. The simulation results show that the fuel economy of the PEHB with the optimal energy management strategy is improved by 53.7% compared with that of the conventional bus, which can be a benchmark for the assessment of other control strategies.
In recent years, the problems of energy shortage and environmental pollution have greatly promoted the development of electric vehicles (EVs). Among the EVs, the pure electric vehicles (PEVs) run with zero emissions and renewable electricity, but their disadvantages, such as the short operation range, high battery price, and long battery charging time, have limited the user’s acceptability. The hybrid electric vehicles (HEVs) have longer operation range and higher performance than PEVs, but the electricity that keeps the battery state of charge (SOC) in a narrow window is still from the onboard fossil fuel [1, 2]. While the plug-in hybrid electric vehicle (PHEVs), with larger battery capacity, can run a long pure electric mileage and make full use of the cheap power from gird, hence it is more competitive than EVs and charge sustainable HEVs .
The energy management strategy is one of the key factors that influence the fuel economy and power performance of the PHEVs. In the PHEVs, in order to make full use of the electricity energy stored in batteries, it is preferred that the battery energy drops to its minimum when the vehicle arrives at the destination. Therefore the energy management strategy becomes more complicated than that of the HEVs. Similar to HEVs, the energy management strategies in PHEVs can be usually classified into two categories: rule-based control strategies and optimization-based control strategies . The main idea of rule-based control strategies is to make each component work in efficient area individually [5–7]. The reference  put forward a PHEV rule-based control strategy after considering the all-electric range and charge depletion range operations. The reference  proposed a PEHV rule-based energy management strategy by using the ADVISOR. The rule-based control strategies are simple and easy to implement, but they cannot safeguard the systematic optimization and cannot fully exploit the advantages of PHEVs.
The optimization-based control strategies include global optimization and real-time optimization. The real-time optimization realizes a local optimum step by step real timely and loses the potential to get a global optimum. The adaptive control is a good example of real-time optimization, and the H∞ control theory is powerful in adaptive control [8–11]. The global optimization finds an optimal solution for the whole process, which is suitable for energy management issues of the PHEB with regular driving cycle. For example, dynamic programming (DP) algorithm, which is effective to solve the constrained and nonlinear optimization problems, is selected to realize a global optimization of energy management for HEVs . The reference  studied an optimal energy management of a parallel HEV with the known driving cycle using DP algorithm. The reference  built the driving cycle model using traffic information with the help of intelligent transportation systems and utilized the DP algorithm to study the global energy management optimization of a parallel plug-in hybrid electric sport utility vehicle (SUV). The reference  presented a way on how to implement the DP algorithm in the optimization of HEVs and carried out a global optimization with Toyota Prius as an example. The reference  determined an optimal energy management law for a two-clutch single-shaft parallel HEV by using optimization software, named KOALA. By comparisons, the DP algorithm has been proved to be powerful and effective in the global optimization of control strategies in HEVs. In this paper, a global optimization of the energy management strategy for a plug-in series-parallel hybrid bus (PHEB) is explored, and the battery energy state control is specially considered and discussed. The PHEB model was built in Section 2 and the global optimization problem with DP algorithm was put forward in Section 3. The DP numerical computation method was discussed and put forward in Section 4 and the simulation results were given in Section 5. Section 6 gives the main conclusions.
2. Plug-In Series-Parallel Hybrid Electric Bus Modeling
2.1. Plug-In Series-Parallel Hybrid Electric Bus Configuration
Figure 1 shows a schematic view of the PHEB powertrain, which includes a diesel engine, an integrated starter generator (ISG) motor, and a main drive motor. The ISG is connected to engine through a torsion damper. There is an on-off mode clutch between the ISG motor and the main drive motor. The PHEB works in parallel mode or engine-only mode when the mode clutch is in “ON” condition, and the diesel engine, ISG motor, and main drive motor drive the wheel mechanically. While the PHEB works in its series mode or all-electric mode when the mode clutch is in “OFF” condition, the main drive motor drives the wheels directly, and the diesel engine drives the ISG motor to generate electricity or not according to the battery state of charge. The main specific parameters of the PHEB are listed in Table 1.
2.2. Plug-In Series-Parallel Hybrid Electric Bus Simulation Models
There are two modeling methods in PHEV simulations. One is forward modeling, which is more accurate but with heavier computational burden, and is always used to test the vehicle dynamic performance and drivability. The other one is the backward modeling, which is calculated with fixed time steps ignoring the dynamics of the powertrain components and usually is used to evaluate the vehicle fuel economy . Since the global optimization is based on the fixed driving cycle and the DP problem is solved backward from the terminal of the driving cycle, we built the facing-backward simulation models as follows.
The diesel engine is modeled as a 3-dimension look-up table, where the inputs are the engine torque and speed and the output is the fuel consumption rate, as shown in Figure 2. The fuel consumption efficiency map is based on the experimental data.
The ISG motor is modeled as a 3-dimension look-up table based on the experimental efficiency map as shown in Figure 3. Due to the limit of battery power, the output torque of ISG is described as follows: where is the required ISG torque, is the ISG speed, SOC is the battery state of charge (SOC), and are the maximum ISG torque when driving and generating, respectively, and and are the torque limits due to battery current limits in the charging and discharging modes, which are functions of ISG speed and torque.
The main drive motor is modeled as a 3-dimension look-up table based on the experimental efficiency map as shown in Figure 4. Due to the limit of battery power, the output torque of the main drive motor is as follows: where is the required torque of the main drive motor, is the speed of the main drive motor, and are the maximum torque when driving and regenerative braking, respectively, and and are the torque limits due to battery current limits when charging and discharging, respectively, which are the functions of the main drive motor speed and torque.
The static equivalent circuit battery model described in  is used. The model inputs are the speed and torque of ISG and main drive motor, the model output is the battery SOC, which is calculated by (3) where is the internal resistance, is the open-circuit voltage, and are the function of , is the maximum battery capacity, is the terminal resistance, and are the efficiencies of the main drive motor and ISG accordingly, and denotes the calculation step in discretization way.
The mode clutch works in three conditions: disengaged, where ; engaged, where ; and half-engaged, where . The transition duration is very small, less than the time interval of the sample point of the driving cycle, so only two working conditions are considered in the dynamic optimization. If , the torque of engine and ISG motor can be delivered to final drive completely. If , only the main drive motor drives the vehicle.
When the PHEB works in series hybrid mode, taking the efficiency of engine and ISG into consideration, the minimum fuel consumption curve can be found as shown in Figure 5. For each required , only the points on the minimum fuel consumption line are considered.
The time interval is set to 1 second. Assuming that the torques of the engine, main drive motor, and ISG remain constant in one time interval, we can calculate the vehicle dynamics as where is the total required torque as the input of the final drive, is the mechanical powertrain efficiency, is the final drive ratio, is the vehicle speed, is the dynamic tire radius, is the effective mass of the vehicle, and and are aerodynamic drag force and rolling resistance force, respectively.
According to (4), we can get the as the driving cycle is known in advance. For the PHEB powertrain, the relationship between and the powertrain components can be expressed as where is the hydraulic brake torque and is the engine torque.
3. Global Optimal Energy Management Strategy Modeling
For PHEB, DP aims to find the control of each stage to minimize the cost function over the whole driving cycles. The control variables and state variables are determined before DP problem is formulated. The state variables, including vehicle speed and SOC, reflect the operating state of the system. As the driving cycle and the vehicle speed in every stage are known, SOC is chosen as the state variable. There are many control variables in the PHEB such as engine torque , engine speed , motor torque , ISG torque , and hydraulic brake torque , but only three of them are independent. Here the , , and are chosen as the independent control variables.
In the discrete-time format, the PHEB system can be expressed as where and are state vector and control vector, respectively.
For PHEB, the price of electricity is very low compared with that of diesel when used to drive the same distance , so we focus our research on minimizing the fuel consumption. The cost function is built as follows: where is the duration of the driving cycle, and is the instantaneous cost; fuel denotes the diesel consumption.
If the SOC drops below the lower limit, the battery will not supply electricity to the main drive motor. To avoid the condition that the engine cannot supply the required torque, another cost function should be considered besides (7) then the cost function rewrites as follows: where is a positive weighting factor.
4. A Numerical Computation for the DP Problem
Based on Bellman’s Principle of Optimization, the global optimization problem can be solved by dealing with a sequence of subproblems of optimization backward from the terminal of the driving cycle [17, 18]. Then the DP problem can be described by the recursive equation (10)-(11). The subproblem for step is For step (), the subproblem is where is the optimal cost-to-go function at state from stage to the end of the driving cycle and is the state in stage after the control is applied to state at stage according to (6).
4.1. Solving the DP Problem Backward
The recursive equation (10)-(11) is solved backward, and quantization and interpolation are needed to solve the equation. The continuous state SOC is discretized into finite grids first, and the number of discretized state is where is the increment of the discretized SOC and and are upper and lower constrains of SOC.
Then find all possible control solutions at every state of each stage. The function at every grid points of SOC is evaluated, and is evaluated by interpolation if the calculated value of admissible in (3) does not fall exactly on grid points. The way of interpolation is shown in . The procedure of solving the DP problem backward is shown in Figure 6, where the required speed is the same as driving cycle and required torque is determined by inversely solving vehicle dynamic model as shown in (4).
4.1.1. To Find All Possible Control Solutions
It is very important to find all possible control solutions in the procedure of solving the DP problems backward. The possible control solutions are the possible combination of the discrete torque of the components in each state of the driving cycle which meets the torque need of the vehicle. The number of the control solutions influences the accuracy of the optimization greatly, and the way to search for all the control solutions influences the computational burden significantly. To get a compromise, here we find the possible working modes of the PHEB first, and then we find all possible control solutions in every mode; finally we get all the control solutions at every grid point of SOC.
The PHEB works in many modes, such as engine-only mode, battery electric (EV) mode, engine-ISG parallel mode, engine-motor parallel mode, and series mode. The ISG is used as a starter and generator and will not drive the bus directly. Only when the torque required by the vehicle exceeds the maximum torque that can be provided by the main drive motor and engine together, the ISG will provide the remaining torque.
According to (6), the state variables in the may exceed the range of SOC, as shown in Figure 7, where the state at stage exceeds the range of SOC with the control variables and . To avoid this situation, the control variables should be limited. We divide SOC into three areas: , , and . The initial SOC and terminal SOC are usually in the area .(A)When the SOC is higher than high limit , the motor drives the bus without regenerative braking. Only when the torque required by the vehicle exceeds the maximum torque that can be provided by drive motor, ISG will supply positive torque to drive the bus. If more driving power is required, the engine comes to work to supply the remaining torque.(B)When the SOC drops below the low limit , the battery will not supply electric energy any more. According to the required torque and required speed areas of final drive shown in Figure 8, the PHEB works in different modes as listed in Table 2. The blue line in Figure 8 represents the maximum output torque of the motor with the power supplied by engine, ISG, when PHEB works in series mode. (C)If , the torque that can be supplied by the powertrain components is shown in Figure 9. The possible working modes are shown in Table 3.
If the PHEB works in series mode, discretize the minimum fuel consumption curve into finite points, and the engine/ISG works on these points. If the PHEB works in engine-motor parallel mode, such as the PHEV which works on area 3 in Figure 9, the flow chart to find all possible controls is shown in Figure 10(a). If the PHEB works in the same mode on areas 4, 5, and 6 in Figure 9, the initial condition of is set to be . If the PHEB works in engine-ISG parallel mode, the way to find the possible controls is shown in Figure 10(b).
(a) For engine-motor parallel mode
(b) For engine-ISG parallel mode
4.2. To Find the Optimal Control Path Forward
The optimal controls at every state point of every stage are obtained by solving the DP problem backward; if the initial SOC is specified, the optimal control path will be found forward. The interpolation is also needed to find the optimal control path as shown in Figure 11. If the optimal control at stage is , the optimal control at stage is got through interpolation between the controls and , which are the optimal controls at state grid points and , respectively, at stage .
5. Simulation Results
For the PHEB, it is reasonable to make full use of the battery energy. Considering the health and the efficiency of the battery, the low level, the high level, and the initial SOCs were selected to be 0.3, 0.8, and 0.6, respectively. The driving cycle is the Chinese typical urban drive cycle (CTUDC) as shown in Figure 12. The total distance of one CTUDC is 5.897 km, and the duration of one CTUDC is 1314 seconds.
The battery capacity of PHEB is much higher than that of HEVs, and the PHEB drives with the mode of one day one charge, so the PHEB would drive for many consecutive driving cycles. To show the energy distribution between the engine, ISG, and main drive motor, the model is simulated with the input of 15 consecutive CTUDC cycles. The increment of discretized SOC, , is selected to be 0.001, the increment of engine torque is selected to be 5 Nm, the weighting factor is selected to be 100, and the PHEB weight of simulation is set to be the gross weight.
Figure 13 shows the simulation result of SOC under the DP optimal control for 15 consecutive CTUDC cycles. The SOC decreased to 0.314 when the bus reached the destination. Because of the regenerative braking at the end of the cycle, the terminal SOC is a little higher than the low level, but it is still very close to the low level, so the bus can make full use of the battery energy with the optimal control.
The SOC decreases evenly from the initial SOC to the low level of the SOC. When the optimal control was applied to the 15 consecutive CTUDC cycles, the SOC reduction would be 0.02 for one cycle on average. To relieve the heavy computational burden, the optimal control for one CTUDC cycle can be used as the optimal control for 15 consecutive CTUDC driving cycles through restricting the initial SOC and terminal SOC. The initial SOC and desired terminal SOC are selected to be 0.5 and 0.48, respectively. To ensure that the SOC at final time is the desired value, an additional terminal constrain on SOC needs to be imposed and the cost function would be where is positive weighting factor and is the desired SOC at the end of driving cycle.
The increment of SOC, , is selected to be 0.001, the increment of engine torque is selected to be 5 Nm, and the weighting factors and are selected to be 100 and 1 × 107, respectively. The difference between simulation result of velocity and desired vehicle velocity is very small as shown in Figure 14. Figure 15 shows the simulation result of SOC with the optimal control for one driving cycle, and the terminal SOC turned out to be 0.4804, which is very close to the desired value.
The simulation results for 15 consecutive driving cycles and one driving cycle are shown in Table 4. The fuel consumption per 100 km of 15 driving cycles is 19.72 L, and the fuel consumption per 100 km of one driving cycle with restricted terminal SOC is 19.90 L. The fuel consumption per 100 km increased by 0.91%, while the computation time decreased by 94.7%. Considering that the internal resistance of the battery is a function of SOC, the optimal control of one driving cycle with restricted terminal SOC is applied to the simulation for 15 consecutive driving cycles, and the simulation result is shown in Table 5.
With the optimal control for one driving cycle, the fuel consumption per 100 km increased by 0.91% and the electric energy consumption increased by 3.2%, but the terminal SOC is still higher than the low level. The computation time of solving DP problem to find the optimal control decreased significantly, so it is feasible to find the optimal controls for consecutive driving cycles by solving the DP problem for one driving cycle with restricted initial SOC and terminal SOC.
The state increment in (12) also influences the accuracy of the optimization. If the is smaller, the quantized search area will be larger, hence the computational burden will be heavier. To study the tradeoff between accuracy of the optimization and computation time, is selected to be 0.0005, 0.001, and 0.005, respectively. The SOC simulation results are shown in Figure 16. The terminal SOC dropped to 0.4802, 0.4804, and 0.4806, respectively, which are very close to the desired value. When is selected to be 0.001, the curve of SOC is very similar to the curve when is selected to be 0.0005.
The fuel consumption and computation time results are summarized in Table 6. Compared with the results when is 0.0005, the fuel consumption increased by 0.61% and the computation time decreased by 53.9% when is 0.001, while the fuel consumption increased by 7.03% and the computation time decreased by 91.0% when is 0.005. Considering the tradeoff between fuel consumption and computation time, it is feasible to set to be 0.001. And the simulation results in this case are shown in Figure 17.
The output torque of the PHEB components is shown in Figure 17. It can be seen that the ISG seldom works as a generator to charge the battery, and most of the negative power is from regenerative braking. If the ISG works as a generator, there are engine efficiency losses, ISG efficiency losses, main drive motor efficiency losses, and battery efficiency losses, hence the system is ineffective, so in most cases the optimal control strategy based on DP avoids the situation when the ISG works as a generator. At the end of the cycle, to force the terminal SOC to be the desired value, the ISG supplied negative power to charge the battery.
The working modes of the PHEB is shown in Figure 18, where means that the bus stops and no powertrain components is working; means regenerative braking mode, means EV mode, means engine-only mode, means parallel mode, and means series mode. It can be seen that the system does not work in series mode in case of the low efficiency.
When the PHEB works in full load condition and the load rate of engine is high, the engine can work in high efficiency area without need of the load regulation of ISG, so the ISG seldom supplies negative torque. The working points of engine and main drive motor are shown in Figure 19. It shows that the engine works in high efficiency areas in most cases. The main drive motor drives the bus alone when the vehicle speed is low. The engine works if the vehicle required speed is high and the main drive motor supplies the remaining torque to maintain the engine working in high efficiency. The ISG seldom works under the optimal control based on DP, but it does not mean that the ISG is useless, because the energy control strategy of the PHEB is based on rules in reality, the driving distance is not a fixed value, and the ISG is needed to charge the battery to maintain the level. There is no standard PHEB fuel consumption for rule-based energy management control strategy, and the PHEB may consume more fuel than traditional bus if the control rule parameters were not properly designed. So we made a comparison of the fuel consumption between the PHEB with optimized control strategy and the prototype conventional diesel bus. The results show that the experimented fuel consumption of the prototype conventional diesel bus is 43 L/100 km in CTUDC driving cycle, and the fuel consumption of the optimized PHEB is 19.90 L/100 km, with additional electricity consumption of 11.61 kWh/100 km, the fuel consumption decreased by 53.7%. The optimal control based on DP can improve fuel economy significantly. It should be noted that the fuel consumption results given by the optimal control based on DP are maximum potential gains, and they cannot be reached in a real vehicle, because the entire driving cycle is known in advance, and neither comfort constraints nor highly dynamic phenomena are taken into account .
It is very complicated to determine the energy management strategy for a series-parallel PHEB, and the dynamic programming is a powerful tool to get global optimization results. The backward simulation model of the series-parallel PHEB was built. Then, to explore the potential of fuel economy, the dynamic programming algorithm is utilized to realize an optimal control on a known-in-advance driving cycle. The procedure of DP for the series-parallel powertrain topology is introduced in detail. An appropriate method is proposed to improve the computational efficiency which can reduce the computation burden greatly and keep the precision of DP.
The simulation results show that with the global optimal control, the battery SOC can reach its lower limit at the end of the cycle, which means that the bus can make full of the battery energy. Meanwhile, the ISG seldom works in generation mode under given cycle and SOC interval, which avoids the inefficient situation. It is proved that the optimal control based on DP can reduce the fuel consumption greatly.
The drawback of optimal control based on DP is that the driving cycle should be known in advance, and the computational burden is still very heavy, so it is difficult to be applied in a real vehicle. In the further study, a near-optimal control law will be extracted according to the global optimization results.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
This work was supported by the National High Technology Research and Development Program of China (2011AA11A228, 2012AA111603, and 2011AA11A290) in part, the International Cooperation Research Program of Chinese Ministry of Science and Technology (2011DFB70020) in part, and the Program for New Century Excellent Talents in University (NCET-11-0785) in part. The authors would also like to thank the reviewers for their corrections and helpful suggestions.
- M. Ehsani, Y. Gao, and A. Emadi, Modern Electric, Hybrid Electric, and Fuel Cell Vehicles: Fundamentals, Theory, and Design, CRC Press, New York, NY, USA, 2009.
- Y. Gao and M. Ehsani, “Design and control methodology of plug-in hybrid electric vehicles,” IEEE Transactions on Industrial Electronics, vol. 57, no. 2, pp. 633–640, 2010.
- C. H. Stephan and J. Sullivan, “Environmental and energy implications of plug-in hybrid-electric vehicles,” Environmental Science and Technology, vol. 42, no. 4, pp. 1185–1190, 2008.
- S. G. Wirasingha and A. Emadi, “Classification and review of control strategies for plug-in hybrid electric vehicles,” IEEE Transactions on Vehicular Technology, vol. 60, no. 1, pp. 111–122, 2011.
- H. Banvait, S. Anwar, and Y. Chen, “A rule-based energy management strategy for plugin hybrid electric vehicle (PHEV),” in Proceedings of the American Control Conference (ACC ’09), pp. 3938–3943, St. Louis, Mo, USA, June 2009.
- N. J. Schouten, M. A. Salman, and N. A. Kheir, “Fuzzy logic control for parallel hybrid vehicles,” IEEE Transactions on Control Systems Technology, vol. 10, no. 3, pp. 460–468, 2002.
- P. Pisu and G. Rizzoni, “A comparative study of supervisory control strategies for hybrid electric vehicles,” IEEE Transactions on Control Systems Technology, vol. 15, no. 3, pp. 506–518, 2007.
- H. Zhang, Y. Shi, and M. Liu, “step tracking control for networked discrete-time nonlinear systems with integral and predictive actions,” IEEE Transactions on Industrial Informatics, vol. 9, no. 1, pp. 337–345, 2013.
- H. Zhang, Y. Shi, and A. S. Mehr, “Robust PID control for multivariable networked control systems with disturbance/noise attenuation,” International Journal of Robust and Nonlinear Control, vol. 22, no. 2, pp. 183–204, 2012.
- H. Zhang, Y. Shi, and A. S. Mehr, “Robust weighted filtering for networked systems with intermittent measurements of multiple sensors,” International Journal of Adaptive Control and Signal Processing, vol. 25, no. 4, pp. 313–330, 2011.
- H. Zhang, Y. Shi, and A. S. Mehr, “Robust static output feedback control and remote PID design for networked motor systems,” IEEE Transactions on Industrial Electronics, vol. 58, no. 12, pp. 5396–5405, 2011.
- C.-C. Lin, H. Peng, J. W. Grizzle, and J.-M. Kang, “Power management strategy for a parallel hybrid electric truck,” IEEE Transactions on Control Systems Technology, vol. 11, no. 6, pp. 839–849, 2003.
- Q. Gong, Y. Li, and Z.-R. Peng, “Trip-based optimal power management of plug-in hybrid electric vehicles,” IEEE Transactions on Vehicular Technology, vol. 57, no. 6, pp. 3393–3401, 2008.
- W. Wang and S. M. Lukic, “Dynamic programming technique in hybrid electric vehicle optimization,” in Proceedings of the IEEE International Electric Vehicle Conference, pp. 4–8, Greenville, SC, USA, March 2012.
- J. Scordia, M. Desbois-Renaudin, R. Trigui, B. Jeanneret, F. Badin, and C. Plasse, “Global optimisation of energy management laws in hybrid vehicles using dynamic programming,” International Journal of Vehicle Design, vol. 39, no. 4, pp. 349–367, 2005.
- W. Zheng, Reseach on power train parameter matching method and control strategy for hybrid electric vehicle [Ph.D. thesis], Harbin Institute of Technology, Harbin, China, 2010.
- Y. Zou, S. J. Hou, and E. L. Han, “Dynamic programming-based energy management strategy optimization for hybrid electric commercial vehicle,” Automotive Engineering, vol. 34, no. 8, pp. 663–668, 2012.
- R. Bellman, Dynamic Programming, Princeton University Press, Princeton, NJ, USA, 1957.
Copyright © 2013 Hongwen He et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.