Mechanics & Industry
Volume 19, Number 3, 2018
|Number of page(s)||11|
|Published online||04 October 2018|
Real-time simulation test-bed for an industrial gas turbine engine’s controller
Systems Simulation and Control Laboratory, School of Mechanical Engineering, Iran University of Science and Technology,
Accepted: 18 May 2018
A hardware-in-the-loop (HIL) test for a control unit of an industrial gas turbine engine is performed to evaluate the designed controller. Although the dynamic performance of the studied gas turbine is strictly related to the variable inlet guide vain (VIGV) position, one of the main challenges is to develop an engine model considering VIGV variations. The model should also be capable of real time simulation. Accordingly, the gas turbine is numerically modeled using bond graph concepts. To demonstrate the operational reliability of the engine’s control strategy, the control algorithm is implemented on an industrial hardware as an embedded system. This is then put into a HIL test along with the engine model. The actual component (controller) and the virtual engine model are the hardware and software parts of the HIL test, respectively. In this experiment, the interaction between the real part and the rest of the system is compared with that of the completely numerical model in which the controller is a simulated software-based model as is the engine itself. Finally, the results indicate that the physical constraints of the engine are successfully satisfied through the implementation of control algorithms on the utilized hardware.
Key words: Bond graph / industrial gas turbine engine / electronic control system / hardware-in-the-loop simulation
© AFM, EDP Sciences 2018
A key factor in designing complex systems and structures is the ability to test every subsystem in each design stage to ensure a convergent design strategy. In this regard, various theoretical and experimental procedures have been considered. The hardware-in-the-loop (HIL) simulation has proven to be one of most efficient method for testing complicated and costly systems. HIL simulation has been used extensively in a variety of fields for real-time testing and development of interconnected physical components of a system replaced virtually by computer models. Interaction between hardware and software during the test is accomplished via electric signals transferred by data acquisition cards. Consequently, HIL simulation is widely used in numerous applied fields and industries. Using this approach, it is possible to test the performance of real mechanical parts of a system along with software-based simulated models of other parts in real time [1,2].
One of the requirements of HIL test is designing and constructing a test-system so as to experiment with different parts of a complex system in accordance with the defined rules. Furthermore, to make HIL more cost-effective, the main system in which we are interested to design should be numerically modeled. These hardware and software models are then put into a loop together for experiment.
Hanselman  benefited from HIL simulation in the control development of electronic control units (ECUs) used in engines, vehicles and other components. Cao et al.  verified the validity of their control scheme based on adaptive network-based fuzzy inference engine using HIL test. Gans et al.  presented another HIL simulation to control unmanned vehicles. In that study, a real camera captured pictures of a virtual 3D environment which would later be used in the control system. Aerospace is another field where HIL has progressively been utilized as Maclay  enumerated several examples. Canadian Space Agency successfully applied HIL in meticulously tuning controllers used in the International Space Station. A non-gravity environment for the controllers was simulated using HIL. The simulation results were acceptable compared to practical tests . In the gas turbine industry, HIL applications are innumerable, specifically for designing, testing and performance verification of the gas turbine engine’s ECU or fuel control unit (FCU). HIL simulation studies have been reported for rapid prototyping of ECU of turbofan engines [8,9] and turbojet engines [10,11]. An HIL simulation is reported by Montazeri-Gh et al.  for ECU performance verification of a turbo-shaft engine. Another HIL simulation is presented by Montazeri-Gh et al.  for testing a fuel control unit of a jet engine.
Despite numerous reports on HIL simulation of a gas turbine engine’s ECU or FCU, there has not been any publication regarding HIL simulation of a two-shaft industrial gas turbine engine modeled using bond graph methodology. As an example of bond graph power in modeling gas turbine engines, Novinzadeh et al. used this method to simulate an ideal turbocharger . Montazeri et al.  showed how bond graph approach can be utilized for modeling the cold start phase of a microjet engine.
Moreover, Krikelis and Papadakis  modeled a simple cycle of the single-shaft gas turbine using bond graph model. By linearization of the model around an operation point, they designed a PI controller for it. In addition, they utilized several parameters including the pressure, temperature and torque as the effort variables, as well as the mass flow rate and engine speed (rpm) as the flow variables.
Sanei et al.  considered the effects of kinetic energy and momentum (in the convergent-divergent nozzles with supersonic fluid flows) using the pseudo-bond graph approach.
Uddin and Gravdahl  developed the bond graph model of a radial compressor system, and complemented it with a control system. They also developed certain methods to prevent the surge in the compressor. Montazeri and Miran-F  presented the bond graph approach application for modeling the industrial gas turbine engine. The bond graph model developed in that study was applicable to the real time implementation. A more comprehensive study using bond graph was later accompanied by the modeling and simulation of the propulsion system of a two-shaft gas turbine including a plate-type clutch  as well as another related study on JetQuard aerial robot . This was similar to other gas turbine nonlinear models [22–24] involving thermodynamics equations as well as compressor and turbine performance maps. Therefore, the bond graph model seems appropriate for controlling purposes during an HIL simulation as is the case in present study.
In this article, the designed control system of a two-shaft industrial gas turbine engine, along with the engine’s ECU (as a physical model), is tested in an HIL simulation. The utilized ECU is an electronic hardware called PC/104, which is an embedded system on which the control algorithm has been implemented via C++ programming language. The engine model is actually a bond graph model simulated on a personal computer using 20SIM and Matlab/Simulink. The interface between hardware and software subsystems is realized by an I/O data acquisition card.
To examine the performance of ECU to ensure the accuracy of its operation, the results are compared with the simulation results of the same system in which the ECU is numerically simulated. Such simulation is called software-in-the-loop (SIL) simulation.
The schematic of the system as well as the interaction of different parts in HIL and SIL test is shown in Figure 1a–c. It is a feedback control system composed of two main subsystems including the plant and the controller. Only the plant is simulated numerically, while the ECU is a hardware within the HIL simulation. The plant, a two-shaft industrial engine (SGT600)  with two outputs and two inputs, is a gas turbine engine model. The engine specifications at the design point are given in Table 1 . The plant’s inputs are the required fuel flow rate to the combustion chamber of the gas turbine engine and the required variable inlet guide vain (VIGV) position. The engine outputs are the angular velocities of the power gas turbine engine shaft (NPT) and the gas generator shaft (NGG). The gas turbine engine has two bleed valves which are considered to be closed during the engine modeling process at a specific state in this article. This state is demonstrated in Figure 2, indicating that their effect is negligible. The plant controller, known as the ECU, calculates the appropriate fuel flow rate and the required VIGV at every moment based on its inputs, and accordingly generates control signals towards the plant. VIGV position is a function of NGG and thus, the effects of VIGV variations are considered in the engine’s performance maps. The ECU has three inputs, two of which are the angular velocities of the engine that are fed back to the ECU, and a reference angular velocity which is set by the user. RRV is a preferred angular velocity with which the power turbine shaft rotates. The outputs are the fuel flow rate and VIGV, as stated previously.
Schematics of (a) complete system, (b) HIL test, and (c) SIL test
The engine components  (compressor, combustor and turbine) are simulated according to the bond graph theory. Every component is assumed to be an energy field interacting with other components, i.e. mass, energy and work are permanently exchanged in all components. The gas turbine engine is considered to be a thermo-fluid system, and the effort and flow variables are selected as introduced by Karnopp [19,26]. According to this variable selection, the pseudo-bond graph is advantageous over the true bond graph. The effort and flow variables of the compressor and turbine are torque and engine speed. Mass flow rate and pressure are the flow and effort variables in one pseudo-bond while energy flow and temperature are the variables for another pseudo-bond. In order to obtain a gas turbine model using the bond graph theory, it is necessary to take the following criteria into consideration:
inclusion of the control model in engine model to observe the physical limitations of the engine;
dynamic operation and steady state performance predictability;
considering the change in the composition of working medium.
Generally, two different approaches are employed to develop the model of compressor or turbine. One is based on the turbo-machinery fundamental equations. Such approach requires not only a thorough understanding of the machine, but also considering some simplifying assumption that may lead to limited accuracy.
The other approach requires taking the state performance maps into account. These maps contain valuable performance information of the machine over a large operational range. Also, no mass and energy accumulation within the compressor and turbine models are considered (quasi-steady assumption). As a result, the performance map for steady state condition is valid even during transient operation.
The isentropic efficiency and corrected mass flow rate of the compressor vary in accordance with three compressor parameters: inlet guide vanes, pressure ratio and corrected rotor speed [19,22]. These are described as (1) (2) (3) (4) where ΓC, πC, NC,cor, ηis,C, MC and are, respectively, the corrected mass flow rate, pressure ratio, corrected rotor speed, isentropic efficiency, torque and air mass flow rate of the compressor. Moreover, NGG is the gas generator speed, Pout is the outlet pressure, (Pin, hin) are the inlet pressure and enthalpies, his,out is the outlet isentropic enthalpy, (Tout, Tin) are the outlet and inlet temperatures and γ is the specific heat ratio. Finally, E being the internal energy, (, ) describe the energy flow in and out of the compressor.
By definition, the dimensionless pressure and temperature are defined as and , respectively, where Pref and Tref represent the standard pressure and temperature (ISA).
For the combustion chamber, the following assumptions are considered: volume of the chamber is constant and physical and chemical properties of fuel and air mixture are the same throughout the chamber [19,22]. Based on the conservation laws of mass (Eq. (5)) and energy (Eq. (6)), the combustion chamber equations are expressed as (5) (6) where is the heat value of fuel, mf is the fuel mass, mair is the air mass and ηcc is the combustor efficiency. By differentiation the above equations, one can obtain (7)where Cv is the specific heat at constant volume. Setting equation (6) equal to equation (7), and then using the ideal gas law, the temperature and pressure of chamber exhaust are described as (8) (9)
In a similar manner, the isentropic efficiency and corrected mass flow rate of the turbine are dependent on two turbine parameters: expansion ratio and corrected rotor speed [19,22], written in the form (11) (12) (13) (14) (15) where g represents a function, i.e. the corrected mass flow rate (ΓT) and isentropic efficiency (ηis,T) are both functions of expansion ratio and corrected rotor speed.
The plenum is taken as an isentropic passage in which energy and flow speed are not significant and thus are neglected. The governing equations to obtain the plenum pressure in addition to the temperature variation caused by the mass accumulation can be written as (16) where Vp is the plenum volume. Moreover, cp and ρ are the constant pressure heat capacity, respectively. δ can be estimated by the specific heat ratio [19,22].
The gas generator shaft acceleration is due to the difference between the turbine output shaft power and the input power to the compressor. It should be noted that the variations in the load applied to the power turbine shaft causes changes in the power turbine speed. This leads to the acceleration of the connector shaft between the power turbine and generator (load). Accordingly, the gas generator shaft along with the power turbine shaft can be described using mathematical relations as in (17) (18) (19) where ηmech,GG and ηmech,PT represent the mechanical efficiency of, respectively, the gas generator and power turbine. In addition, (Mfric,GG, MGG.Turb) and (Mfric,PT, MPower.Turb) are the friction and turbine torque of, respectively, the gas generator and power turbine. Furthermore, ML signify the consumed torque as a result of applied load on the power turbine shaft.
Figure 3 shows the component models. The compressor, combustion chamber and turbine are modeled via modulated energy fields (MR, MC and MR). IGV and fuel-to-air-ratio (f) signals are transferred to the compressor and turbine by informative bonds. When the bond graphs of all sub-models are coupled, a complete gas turbine engine dynamic model would be constructed, as shown in Figure 4.
Complete bond graph model of the engine.
ECU is the master mind of the system. It computes the necessary amount of required fuel as well as the appropriate IGV position for the engine to provide a satisfactorily operation and ensure a safe performance. The fuel control algorithm is based on the Min-Max control strategy. The IGV control algorithm is a function of the angular velocity of the gas generator shaft. Upon designing the control algorithm in 20SIM  and Matlab/Simulink, its precise and non-destructive performance is tested by the computer. As in every gas turbine engine design and construction process, ECU eventually needs to be implemented on a hardware. In this study, a microprocessor called PC/104 has been used. In addition, the Min-Max control strategy for fuel controlling along with the IGV control function are implemented using C++ language. The models of plant  and ECU [12,28] as well as PC/104 specifications used for this study were discussed in detail in previous studies.
The Min-Max controller is composed of five transient control loops and a single steady state control loop, as depicted in Figure 5. In each loop, the required fuel to fulfill the needs of that loop is calculated. Next, a Min-Max algorithm is used for fuel selection throughout the engine operation. This algorithm is written as 
The description of each parameter is presented in Table 2.
In order to assess the performance of this controller, an input load is applied to the engine, as demonstrated in Figure 6. As shown in Figures 7–9, NGG (NGG), NPT (NPT) and their rate of change have been kept within allowable limits. Also, NPT approximately remains at its desired set value.
Load input to the gas turbine model.
Variation of the gas generator turbine shaft acceleration.
Variation of the gas generator shaft speed.
Variation of the power turbine shaft speed.
In order to evaluate the performance accuracy of the ECU and its proper implementation, a real time HIL simulation test-bed is prepared as displayed in Figure 10.
The gas turbine engine model, created using 20SIM and Simulink/Matlab, is loaded on a PC labeled “1” in the figure. The control algorithm is embedded on VDX6354 PC/104 (labeled “2”) as the ECU, via a C++ code.
The VDX-6354 family of controllers is designed as a plug-in replacement to support legacy software and help extend the existing product life cycle without heavy re-engineering. VDX-6354 is suitable for a broad range of data-acquisition tasks, industrial automation, process control, automotive controller, AVL, intelligent vehicle management device, medical device, human machine interface, robotics, machinery control, in addition to applications that require small footprint, low-power and low-cost hardware with open industry standard such as PC/104. Transmission of signals from software-based engine model to the hardware-based ECU model and vice versa is rendered by a data acquisition card, Advantech PCL-812PG I/O (labeled “3”). PC/104 VGA output is connected to a monitor (labeled “4”) to enable the user to edit the C++ code of the control algorithm.
HIL simulation of the test-bed.
To conduct an HIL simulation, an input load is designated by the user. In this study, the load is selected to be a combination of some ramp and step inputs, as indicated in Figure 11. The load starts from its minimum value and remains constant for about 20 s, then increases to its maximum load value in two consecutive steps after 60 s and finally returns to its initial value. This should be considered as the worst-case scenario of load since it involves acceleration along with deceleration for a very short period of time with a relatively steep trend. Should the assumed control strategies keep the system operating in the desired range, they will most likely do for other load function.
The HIL simulation results are compared with those of SIL simulation when subjected to the same load.
Figure 12 shows the result of NGG in SIL and HIL simulations. It shows that NPT closely follows the load trend, and a rise or fall in the load will cause the same effect in NGG plot. In addition, NGG changes based on the change in load with virtually no observed lag for both case of SIL and HIL simulations. This indicates that the control strategies on NGG are satisfactory and the implementation is properly carried out for HIL test.
Normalized NPT from HIL and SIL simulations are illustrated in Figure 13. NPT tends to remain at its desired set value (RRV). However, the ECU designates a new fuel demand for the engine when a sudden change in the load occurs. As a result, some oscillations occur in NPT until it reaches its defined value. Fuel mass flow rate for HIL simulation and SIL tests are demonstrated in Figure 14. Upon load deviation, the fuel flow to the engine is decreased and vice versa, as indicated by the simulation results.
Sampling time plays a key role in HIL simulations for the physical part of the test system, i.e. PC/104 microprocessor. Choosing an appropriate sampling time for every HIL simulation is an essential part of the test. Figure 15 shows the effect of different sampling times (53 ms, 75 ms and 102 ms) for the same HIL simulation described in the previous section.
As can be observed, the higher the sampling time, the more oscillatory the test results are. Nevertheless, relative accuracy is reached for the sampling time of about 53 ms. To put it more accurately, sampling time greatly affects the results stability. By choosing a small sampling time, the simulation time for a sampled data may not match the require time and addition of new sampled data may end in results divergence. A similar discussion can be made on long simulation time where the simulation may be left with no data input for a period, which may lead to instability. There are recommendations for choosing sampling time, but the best one is obtained through trial and error, as carried out in this study.
Applied load to the HIL simulation.
HIL and SIL simulation results for NGG signal.
HIL and SIL simulation results for NPT signal.
HIL and SIL simulation results for fuel flow signal.
Fuel flow signal for three HIL simulations with different sampling times.
In this article, an industrial gas turbine engine with two shafts and variable IGVs was studied. The VIGV position and the engine inlet fuel were considered as the controlling parameters of the engine. Since the performance characteristic of the engine is highly dependent on VIGVs, a modeling procedure was chosen to take this effect into account (bond graph modeling). The ECU of this gas turbine engine provides the required fuel flow as well as VIGV operational position. To evaluate the accuracy of the designed control system, the control strategy was implemented on an electronic hardware and tested in real time via some comprehensive HIL simulations. Every physical constraint of the engine is satisfied by fuel and VIGV position, indicating the successful implementation of the control algorithm on the PC/104 hardware. Changes in bleed valves can be the subject and future studies, while HIL test can be carried out at lower speed and in the start-up phase of gas turbine. A small rounding of the load ramp corner is also necessary to obtain better results. Finally, a more advanced control system such as model predictive control can be used to perform HIL test.
- H.K. Fathy, Review of hardware-in-the-loop simulation and its prospect in the automotive area, Society of Photo-optical Instrumentation Engineers, Proc. SPIE 6228 (2006) [Google Scholar]
- R. Isermann, J. Schaffnit, S. Sinsel, Hardware-in-the-loop simulation for the design and testing of engine-control systems, Control Eng. Pract. 7 (1999) 643–653 [CrossRef] [Google Scholar]
- H. Hanselman, Hardware-in-the-loop simulation testing and its integration into a CACSD toolset, in: The IEEE International Symposium on Computer-Aided Control System Design, Dearborn, Michigan, USA, 1996, pp. 152–156 [EDP Sciences] [Google Scholar]
- Y. Cao, W. Teng, H. Zhang, Hardware-in-the-loop simulation for engine idle speed control based on Adaptive Neural Fuzzy Inference Engine, IEEE, 2008 [Google Scholar]
- N.R. Gans, W.E. Dixon, R. Lind, A. Kurdila, A hardware in the loop simulation platform for vision-based control of unmanned air vehicles, Mechatronics 19 (2009) 1043–1056 [CrossRef] [Google Scholar]
- D. Maclay, Si mulation gets into the loop, IEE Review (1997) 109–112 [CrossRef] [Google Scholar]
- J.C. Piedboeuf, M. Doyon, R. L’Archeveque, E. Martin, Simulation environments for space robot design and verification, in: Advanced Space Technologies for Robotics and Automation, Noordwijk, The Netherlands, 2000 [Google Scholar]
- H. Wang, Y. Guo, J. Lu, Design and validation of aeroengine control system with non-fully recovering LQG/LTR method, in: Second International Conference on Innovative Computing, Information and Control, ICICIC0, 2007, Art. No. 4428111, 2008 [Google Scholar]
- L. Jun, G. Ying-Qing, W. Hai-Quan, Rapid prototyping real-time simulation platform for digital electronic engine control, in: Second International Symposium on Systems and Control in Aerospace and Astronautics, ISSCAA, 2008, Art. No. 4776230, 2008 [Google Scholar]
- M. Montazeri-Gh, M. Nasiri, S. Jafari, Real-time multi-rate HIL simulation platform for evaluation of a jet engine fuel controller, Simul. Model. Pract. Theory 19 (2011) 996–1006 [CrossRef] [Google Scholar]
- A. Watanabe, S.M. Ölçmen, R.P. Leland, K.W. Whitaker, L.C. Trevino, C. Nott, Soft computing applications on a SR-30 turbojet engine, Fuzzy Sets Syst. 157 (2006) 3007–3024 [CrossRef] [Google Scholar]
- M. Montazeri-Gh, S. Abyaneh, S. kajemnejad, Hardware-in-the-loop simulation of two-shaft gas turbine engine’s electronic control unit, Proc. Inst. Mech. Eng. Part I: J. Syst. Control Eng. 230 (2016) 512–521 [CrossRef] [Google Scholar]
- M. Montazeri-Gh, M. Nasiri, M. Rajabi, M. Jamshidfard, Actuator-based hardware-in-the-loop testing of a jet engine fuel control unit in flight conditions, Elsevier, Simul. Model. Pract. Theory 21 (2012) 6577 [Google Scholar]
- A.S. Movaghar, A. Novinzadeh, Ideal Turbo Charger Modeling and Simulation Using bond graph Approach. in: ASME 2011 Turbo Expo: Turbine Technical Conference and Exposition, American Society of Mechanical Engineers, 2011, pp. 871–879 [Google Scholar]
- M. Montazeri-Gh, S.A. Miran-F, Application of bond graph method in microjet engine cold start modeling to investigate the idea of injecting compressed air, in: Applied Mechanics and Materials, Trans Tech Publications, Vol. 799, 2015 [Google Scholar]
- N.J. Krikelis, F. Papadakis. Gas turbine modelling using pseudo-bond graphs, Int. J. Syst. Sci. 19 (1988) 537–550 [CrossRef] [Google Scholar]
- A. Sanei, A.B. Novinzadeh, M. Habibi. Addition of momentum and kinetic energy effects in supersonic compressible flow using pseudo bond graph approach. Math. Comput. Model. Dyn. Syst. 20 (2014) 491–503 [CrossRef] [MathSciNet] [Google Scholar]
- N. Uddin, J.T. Gravdahl, Bond graph modeling of centrifugal compression systems. Simulation 91 (2015) 998–1013 [CrossRef] [Google Scholar]
- M. Montazeri-Gh, S.A. Miran-F, Application of bond graph approach in dynamic modelling of industrial gas turbine, Mech. Ind. 18 (2017) 410 [CrossRef] [Google Scholar]
- M. Montazeri-Gh, S.A. Miran-F, Modeling and simulation of a two-shaft gas turbine propulsion system containing a frictional plate–type clutch, Proc. Inst. Mech. Eng. Part M: J. Eng. Maritime Environ. (2018), DOI: 10.1177/1475090218765378 [Google Scholar]
- S.A. Miran-F, M. Montazeri-Gh, Modeling and simulation of Jet Quad aerial robot, knowledge-based engineering and innovation (KBEI), in: 2017 IEEE 4th International Conference on. IEEE, 2017 [Google Scholar]
- S. Camporeale, B. Fortunato, M. Mastrovito, A modular code for real time dynamic simulation of gas turbines in simulink, ASME J. Eng. Gas Turbines Power 128 (2006) 506–517 [CrossRef] [Google Scholar]
- R. Chacartegui, D. Sánchez, A. Muonoz, T. Sánchez, Real time simulation of medium size gas turbines, Energy Convers. Manag. 52 (2011) 713–724 [CrossRef] [Google Scholar]
- N.U. Rahman, J.F. Whidborne. Real-time transient three spool turbofan engine simulation: a hybrid approach. J. Eng. Gas Turbines Power 131 (2009) 051602 [CrossRef] [Google Scholar]
- “SGT-600 Industrial Gas Turbine,” Siemens Industrial Turbomachinery, Duisburg, Germany, 2005, http://www.energy.siemens.com/ru/pool/hq/power-generation/gas-turbines/SGT-600/downloads/SGT-600_GT_PowerGen_EN.pdf [Google Scholar]
- D.C. Karnopp, D.L. Margolis, R.C. Rosenberg, System dynamics: modeling, simulation, and control of mechatronic systems, John Wiley and Sons, Hoboken, NJ, 2012 [CrossRef] [Google Scholar]
- J.F. Broenink, 20-sim software for hierarchical bond graph/block-diagram models. Simul. Pract. Theory 7 (1999) 481–492 [CrossRef] [Google Scholar]
- M. Montazeri-Gh, S. Abyaneh, Real-time simulation of a turbo-shaft engine's electronic control unit, Mech. Ind. 18 (2017) 503 [CrossRef] [Google Scholar]
- https://www.slideshare.net/AliRafiei2/gas-turbine-training-power-point-sample. Last check: 2017-08-19 [Google Scholar]
Cite this article as: M. Montazeri-Gh, S.A. Miran Fashandi, S. Abyaneh, Real-time simulation test-bed for an industrial gas turbine engine’s controller, Mechanics & Industry 19, 311 (2018)
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.