Research Article  Open Access
S. H. Gawande, L. G. Navale, M. R. Nandgaonkar, D. S. Butala, S. Kunamalla, "Fault Detection of Inline Reciprocating Diesel Engine: A Mass and GasTorque Approach", Advances in Acoustics and Vibration, vol. 2012, Article ID 314706, 6 pages, 2012. https://doi.org/10.1155/2012/314706
Fault Detection of Inline Reciprocating Diesel Engine: A Mass and GasTorque Approach
Abstract
Early fault detection and diagnosis for mediumspeed diesel engines are important to ensure reliable operation throughout the course of their service. This work presents an investigation of the diesel engine combustionrelated fault detection capability of crankshaft torsional vibrations. Proposed methodology state the way of early fault detection in the operating sixcylinder diesel engine. The model of six cylinders DI Diesel engine is developed appropriately. As per the earlier work by the same author the torsional vibration amplitudes are used to superimpose the mass and gas torque. Further mass and gas torque analysis is used to detect fault in the operating engine. The DFT of the measured crankshaft’s speed, under steadystate operating conditions at constant load shows significant variation of the amplitude of the lowest major harmonic order. This is valid both for uniform operating and faulty conditions and the lowest harmonic orders may be used to correlate its amplitude to the gas pressure torque and mass torque for a given engine. The amplitudes of the lowest harmonic orders (0.5, 1, and 1.5) of the gas pressure torque and mass torque are used to map the fault. A method capable to detect faulty cylinder of operating Kirloskar diesel engine of SL90 EngineSL8800TA type is developed, based on the phases of the lowest three harmonic orders.
1. Introduction
The great interest of engine manufacturing industry and research centers on the crankshaft vibration [1] and instantaneous speed analysis for the fault or misfire detection in internal combustion engines is due to the positive results obtained using different approaches. A misfire event can be considered as an impulsive excitation on the crankshaft, related to the lack of torque resulting from the missing combustion. This excitation causes a sudden engine speed decrease followed by an oscillation with characteristics amplitude and frequency associated with the torsional behavior of the engineload system.
A number of methods are devised to detect nonuniformities in the contributions of the cylinders to the total engine output by analyzing the variation of the measured crankshaft’s speed may be classified in different groups such as observerbased method, pattern recognition, statistical estimators, order domain methods, and a combination of these methods. In [2], a good overview of the earlier state of the art has been presented. The model developed in [3] has been extensively used with very good results to predict and control torsional vibrations to correlate the amplitudes of the major harmonic orders of the crankshaft’s speed variation with the average indicated mean effective pressure of the engine. Balancing the cylinderwise torque contributions of automotive and other highspeed engines is usually addressed by reconstructing the cylinderwise netindicated torques, by direct use of the angular acceleration [4]. The method suggested [5, 6], in which the measured angular speed was directly related to the nonuniform torque contribution to detect misfire. In order to determine the cylinderwise netindicated torque, the oscillating torque applied on the flywheel was reconstructed from measurements of the angular speeds of the crankshaft. The oscillating torque can be reconstructed by using a single mass engine model and one angular speed measurement [7] by assuming that the crankshaft is rigid and that the engine is sufficiently decoupled from the transmission and load,. Obtaining an estimate of the prevailing load torque [7, 8] uses the fact that the instantaneous torque from the engine is zero at topdead center (TDC) and bottomdeadcenter (BDC). In [9], misfire was simulated by disconnection of fuel supply of one cylinder of the engine which describes the measurements of vibroacoustic exhaust locomotive engine signals for misfire phenomenon simulation and some results of their nonlinear analysis. Instead of estimating the engine load torque [10] proposed an observerbased method, which uses the measured engine load torque directly to reconstruct the cylinderwise netindicated torques of a sixcylinder engine. In [11], a novel method, based on the polar coordinate system of the instantaneous angular speed (IAS) waveform, has also been introduced: to improve the discrimination features of the faults compared to the FFTbased approach of the IAS signal. In [11], two typical experimental studies were performed on 16 and 20cylinder engines, with and without faults, and diagnose the results by the proposed polar presentation method. Method for evaluating the torque nonuniformity between the various cylinders of an internal combustion engine due to statistical dispersion in manufacturing or aging in the injection system is proposed in [12], where the problem of the crankshaft torsional vibration effects is presented on the engine speed fluctuations, that make the misfire detection critical, when multiple misfires are present in the same cycle. In [13], a simple algorithm based on an energy model of the engine which requires the measurement of the instantaneous angular speed for misfire detection in reciprocating engines is proposed. By processing the engine dynamics in the angular domain, variations in the working parameters of the engine, such as external load and mean angular speed, are compensated. A dimensionless feature has been abstracted for evaluation of the combustion as well as compression process of each cylinder.
In this work, a method of detecting early fault for mediumspeed power plant diesel engines is proposed, in which as in [14, 15] torsional vibration amplitudes are used to superimpose the mass and gas torque. Further mass and gas torque analysis is used to detect fault in the operating engine.
2. Problem Definition
Problem studying in this work is to identify the early faulty cylinder by using the amplitudes of torsional vibration and lowest harmonic orders (0.5, 1, and 1.5) of the inertia masses and DFT of measured engine speed signal based on gas forces. In this paper, the authors focus and present how amplitudes of torsional vibration play an important role in performance diagnosis of internal combustion engine which will be the recent area of research for engine manufacturing industry and research centers.
3. Engine Model
Figure 1 illustrates a complete masselastic model of a fourstroke sixcylinder inline diesel engine system. There is a flywheel on the engine shaft to even the rotational speed of the engine shaft. The rectangles in Figure 1 represent masses with link rotating in relation to the shafts. For example, the mass of the piston/crank mechanism of each cylinder is illustrated by rectangles. Dampings are represented by the box/plate symbol and stiffness of shafts is represented by spring symbol. In Figure 1, are the inertia of the each cylinder and pulley, are the stiffness, and are the viscous damping’s, respectively.
4. Results and Discussion
In this work, as per [1, 6, 14, 15], detailed torsional vibration analysis was carried out and is used to calculate the torsional vibration amplitude for the given operating engine. The information in Table 1 pertains to the fourstroke, sixcylinder inline diesel engine. Tables 2, 3, and 4 show detailed torsional vibration analysis for the engine under consideration.




5. Interpretation and Calculation of Resultant Phase Vector Based on Reciprocating Masses
Here, the engine under consideration is SL90 model manufactured by Kirloskar Oil Engine Pune03. For the considered sixcylinder engine used as an example, the crankshaft has crank throws spaced 120° apart. The firing sequence selected is l53624. Considering the order number is the number of complete cycles occurring during one crankshaft revolution, and then the firstorder vectors will be spaced 120° apart in the sequence corresponding to the firing order. secondorder vectors rotate at twice crank intervals so each succeeding secondorder vector will be spaced at 2 × 120° = 240° apart in their correct sequence and so on.
5.1. Torsional Resonance Calculation by Holzer Tabulation Method
A phase vector diagram as shown in Figure 2 is drawn for torsional resonance of 1387.22 rpm for the given engine to analyze the minor and major critical harmonic orders. It is seen that the third order is the major critical orders of excitation for calculated torsional resonance. The corresponding critical speeds of engine operation would be;
Since the 3rd order excitation of 3node vibration mode falls within the operating range of 750 to 2200 rpm, hence a phase vector diagram as shown in Figure 2 is drawn for torsional resonance of 1387.33 rpm or 4162 vpm (Vibration per Minute).
Figures 2(a) and 2(b) shows the phase angle diagram with resultant vectors for the 0.5 order and l order with phase angles 60° and 120°. The cylinder number and the amplitude vectors corresponding to it are shown in Figure 2 are drawn by using data in Table 5. It is found that more than one vector has the same angular position, so these have been arithmetically added making the diagram to the right. Then, by taking cosine components and algebraically adding the total cosine (vertical) component for the resultant vector is found, and by taking the sine components the resultant sine component is similarly found. The resultant vector magnitude is then It is found that the fourthorder successive vectors are spaced at 4 × 120 = 480° = (360 + 120°) apart and, therefore, have the same pattern as for the first order, similarly for the 7th and 10th orders. Figure 2(a) shows the0.5, 2.5, 3.5, 5.5, 6.5 orders. The 1.5, 4.5, 7.5 orders as shown in Figure 2(c) all act up and down, as drawn, and have a moderately large resultant. In the case of order numbers which are multiples of half the number of cylinders, that is, 3, 6, 9, and so forth, for the considered sixcylinder engine, all the vectors act in the same direction (not shown in Figure 2), for this reason, they are called “Major Orders.”

For 0.5 Order
Phase vector angle, (sum of all vertical components/sum of all horizontal components) °.
For 1 Order
Phase vector angle, (sum of all vertical components/sum of all horizontal components) °.
For 1.5 Order
Phasevector angle (sum of all vertical components/sum of all horizontal components) °.
For given engine the resultant phase vector can be calculated as above, and this Resultant phase Vector can be represented (by arrow having magnitude equal to amplitude of torsional vibration) as in Figures 2(a), 2(b), and 2(c).
6. Interpretation and Calculation of Resultant Phase Vector Based on Gas Forces
Extensive experiments were conducted on a fourstroke six cylinder direct injection diesel engine (Kirloskar SL90SL8800TA). The engine was operated at constant speed and full loads. To simulate a early faulty cylinder, the engine was operated under normal working condition. The pressures were measured in all cylinders by piezoelectric pressure transducers. The meanindicated pressure (MIP) and the gas pressure torque (GPT) of each cylinder were calculated from the pressure traces. The GPT and the measured speed were subjected to a discrete Fourier transform (DFT) to determine the amplitudes and phases of their harmonic components. Figure 3 shows an actual pressure curve generated by the pressure transducers on all six cylinders when engine was working under normal condition. The results of the application of this technique to given DI Diesel engine are graphically shown in Figures 4 and 5.
(a) θ = 11°
(b) θ = 110°
(c) θ = 40°
It is observed that when the cylinders are uniformly contributing to the total engine torque, the first three harmonic orders (, 1, 1.5) play a significant role in the frequency spectrum of the total gaspressure torque and, consequently, appear with a very low contribution in the frequency spectrum of the crankshaft’s speed [5, 6]. If the frequency spectrum of the crankshaft’s speed corresponding to uniform cylinders operation is compared to the spectrum corresponding to a faulty cylinder, one may see that the major difference is produced by the amplitudes of the first three harmonic orders. As far as the cylinders operate uniformly, these amplitudes are maintained under a certain limit. Once a cylinder starts to reduce its contribution, the amplitudes of the first three harmonic orders start increasing. These amplitudes may be used to determine the degree by which a cylinder reduces its contribution to the total gas pressure torque. The identification of the faulty cylinder may be achieved by analyzing the phases of the lowest three harmonic orders.
Figure 4 is drawn by reconstructing the pressure traces of the six cylinders in a sequence corresponding to the firing order (153624) with engine operating under normal conditions. Figure 5 shows the lowest three harmonic orders of the measured speed, respecting the measured amplitudes and phases. It is seen that only for the expansion stroke of cylinder 5, all three harmonic curves have, simultaneously, a negative slope. In the phaseangle diagrams (Figure 5) of these orders, the resultant vector corresponding to the harmonic component of the measured speed is also represented (shown by arrow having magnitude equal to amplitude of vibration) is drawn by using data in Figure 6. One may see that, for each of the three considered orders, the vectors are pointing toward the group of cylinders that produces less work. The cylinder that is identified three times among the less productive cylinders is the faulty as shown in Figure 5.
7. Method for Detecting Fault in Operating Engine
Based on the Figures 2 and 5, the following method is developed.
(1) The phaseangle diagrams, considering the firing order of the engine are drawn for the lowest three harmonic orders placing in the top dead centre (TDC) the cylinder that fires at 0° in the considered cycle.
(2) On these phase angle diagrams, the corresponding vectors of the measured speed are represented in a system of coordinate axes.
(3) The cylinders toward which the vectors are pointing are the less contributors and receive a “” mark. If there are cylinders that receive a “” mark for all three harmonic orders they are clearly identified as less contributors to the engine total output.
(4) This method is able to identify a faulty cylinder at very early stage, that is, before starting the engine manufacturing (Table 6) and as soon as its contribution reduces from rated value to zero with respect to the contribution of the other cylinders (Table 7).


8. Conclusions
Proposed methodology state the way of early fault detection in the operating diesel engine and Tables 6 and 7 shows validation of developed methodology. Here, torsional vibration amplitudes Table 2 are used to detect early fault in the operating diesel engine. In this paper, combination of mass torque and gas torque is used to map the fault. From this work it is found that the amplitudes of the lowest harmonic orders (0.5, 1, and 1.5) of the inertia masses and DFT of measured speed based on gas forces may be effectively used to map the fault. This method is capable to detect early faulty cylinders (which is now cylinder number 5) based on the phases of the lowest three harmonic as soon as its contribution starts to drop from rated value to zero due to its inertia forces as compared to the contribution of the other cylinders. This methodology can be also used for early fault diagnosis in the operating SI engine.
References
 E. J. Nestrides, A Handbook on Torsional Vibration, B.I.C.E.R.A, Research Library, Cambridge University Press, Cambridge, UK, 1958.
 J. Williams, “An overview of misfiring cylinder engine diagnostic techniques based on crankshaft angular velocity measurements,” SAE Technical Paper 960039, 1996. View at: Google Scholar
 K. S. Chen and T. Chang, “Crankshaft torsional and damping simulationan update and correlation with test results,” SAE Paper 861226, 1986. View at: Google Scholar
 M. J. Van Nieuwstadt and I. V. Kolmanovsky, “Detecting and correcting cylinder imbalance in direct injection engines,” ASME Journal of Dynamic Systems, Measurement and Control, vol. 123, no. 3, pp. 413–424, 2001. View at: Google Scholar
 D. Taraza, N. Henein, and W. Bryzik, “Diesel engine diagnostics based on analysis of the crankshaft's speed variations,” SAE Paper 982540, 1998. View at: Google Scholar
 S. H. Gawande, L. G. Navale, M. R. Nandgaonkar D, S. Butala, and S. Kunamalla, “Detecting power imbalance in multicylinder inline diesel engine genst set,” Journal of Electronic Science and Technology, vol. 8, no. 3, pp. 273–279, 2010. View at: Google Scholar
 U. Kiencke and L. Nielson, Automotive Control Systems, for Engine, Driveline and Vehicle, Springer, Italy, 2005.
 U. Kiencke, “Engine misfire detection,” Control Engineering Practice, vol. 7, no. 2, pp. 203–208, 1999. View at: Publisher Site  Google Scholar
 P. Boguś and J. Merkisz, “Misfire detection of locomotive diesel engine by nonlinear analysis,” Mechanical Systems and Signal Processing, vol. 19, no. 4, pp. 881–899, 2005. View at: Publisher Site  Google Scholar
 Y. W. Kim, G. Rizzoni, and Y. Y. Wang, “Design of an IC engine torque estimator using unknown input observer,” Transactions of ASME Journal of Dynamic Systems Measurement and Control, vol. 121, no. 3, pp. 487–495, 1999. View at: Google Scholar
 P. Charles, J. K. Sinha, F. Gu, L. Lidstone, and A. D. Ball, “Detecting the crankshaft torsional vibration of diesel engines for combustion related diagnosis,” Journal of Sound and Vibration, vol. 321, no. 3–5, pp. 1171–1185, 2009. View at: Publisher Site  Google Scholar
 N. Cavina and F. Ponti, “Engine torque nonuniformity evaluation using instantaneous crankshaft speed signal,” ASME Journal of Engineering for Gas Turbines and Power, vol. 125, no. 4, pp. 1050–1058, 2003. View at: Publisher Site  Google Scholar
 F. V. Tinaut, A. Melgar, H. Laget, and J. L. Domínguez, “Misfire and compression fault detection through the energy model,” Mechanical Systems and Signal Processing, vol. 21, no. 3, pp. 1521–1535, 2007. View at: Publisher Site  Google Scholar
 S. H. Gawande, L. G. Navale, M. R. Nandgaonkar, and D. Butala, “Torsional frequency analysis of multicylinder inline diesel engine generator system,” in Proceedings of the International Conference on Mechanical, Industrial, and Manufacturing Technologies (MIMT '10), pp. 519–524, ASME, January 2010. View at: Google Scholar
 S. H. Gawande, L. G. Navale, M. R. Nandgaonkar, D. Butala, and S. Kunamalla, “Cylinder imbalance detection of six cylinder DI diesel engine using pressure variation,” International Journal of Engineering Science and Technology, vol. 2, no. 4, pp. 433–441, 2010. View at: Google Scholar
Copyright
Copyright © 2012 S. H. Gawande 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.