Rotating systems misalignment identification using fuzzy clustering method

Pająk, Michał; Muślewski, Łukasz; Landowski, Bogdan; Grządziela, Andrzej

Vibroengineering Procedia

Browse Procedia

Published: 31 August 2016

Check for updates

Rotating systems misalignment identification using fuzzy clustering method

Michał Pająk¹

Łukasz Muślewski²

Bogdan Landowski³

Andrzej Grządziela⁴

¹Department of Thermal Technology, University of Technology and Humanities, Radom, Poland

^{2, 3}Machine Maintenance Department, University of Science and Technology, Bydgoszcz, Poland

⁴Polish Naval Academy, Mechanical Electrical Faculty, Gdynia, Poland

Corresponding Author:

Michał Pająk

Cite the article Download PDF

Downloads 1246

Abstract

In the phase of operation, a technical system accomplishes goals which were the reason of the system designing and creation. The execution of operation tasks by the system causes changes in the system state. The changes in the dynamic state of rotating systems can be early identified by means of vibrations measurement. Therefore, the operational tests were performed in order to register time histories of vibrations. Next, characteristics of the collected vibration signals were recorded and the most significant of them were taken as an input for further studies. The studies were performed to formulate a method of identifying the rotating systems misalignment with high reliability. The method consists in fuzzy clustering implementation. According to the method the input signals were used to formulate the fuzzy clusters of the ability and inability states of the system. Thanks to the method, it is possible to identify the operation position of the system on base of the vibration signals analysis. In the paper the proposed method was described. The results of the method application for a gas turbine engine as an example of a rotating system are also presented. It should be emphasized that the reliability of the described operation state identification is higher than other methods implemented by the authors previously.

1. Introduction

Vibrations of propulsion units of transport means are an important research issue to be explored in terms of their operation processes. They are the main diagnostic symptoms from the point of view of damageability and reliability of these subsystems functioning. They have also effect on the comfort of transport means operators and passengers.

Vibrations and noise produced by rotating subsystems occur in practically all mechanical transport means. Depending on the transport branch and the type of object vibrations are produced by: drive shafts, axes, turbines and fans. On the basis of literature, it was found that early symptoms of rotating systems dynamic changes can be detected mainly by the system vibration measurement in particular its misalignment [1-4].

2. Characteristics of the research object

The research object of this study are eight propulsion systems of selected transport means. Each of them consists of: nominal rotational speed equal to 1550 rpm and nominal power 736 kW, hydro-kinetic clutch with nominal skid 2 % and skid control range 2-98 %, reverse reduction transmission with ratio 3.5:1 and resistance bearing equipped with three rolling bearings.

Vibrations of the main engine shafting were measured in six points of each propulsion system situated on main engine foundation, power end of the engine, a hydro-kinetic clutch, a bearing on the reverse reduction transmission inlet; a bearing on the reverse reduction transmission outlet and a resistance carrier bearing. Three identical accelerometers (B&K 4514B) were used for measurements.

They were mounted on three mutually perpendicular axes. As there was no possibility to use threaded joints the converters were fixed to the tested components of the propulsion system by means of glue. Measurement tracks were calibrated before and after measurements. All the recorded signals were synchronized by means of a four channel measurement gauge [5].

Vibrations were recorded within the range of band 3.2 kHz and for sampling frequency 8192 Hz.

As has already been mentioned the signals were collected from eight objects where vibrations were recorded separately for 4 rotational speeds (850 rpm, 1100 rpm, 1300 rpm and 1500 rpm) and three directions perpendicular to each other. The results were recorded in cvs format files. Each file included values of vibrations for a given transport means, shaft system, measurement point and rotational speed [2, 6].

The above data provided the basis for further studies of this article.

3. Analysis of experimental tests results

The values of registered vibrations were treated as discrete signals expressed in the time domain. For each signal ( $v s_{i}$ ) sixteen characteristics ( $C H (v s_{i})$ ) were calculated [6-9]:

– the signal integral – $I (v s_{i})$ ,

– the signal mean value – $M (v s_{i})$ ,

– the signal energy – $E (v s_{i})$ ,

– the signal mean power – $P (v s_{i})$ ,

– 1st moment of the signal – $M_{1}^{'} (v s_{i})$ ,

– 2nd moment of the signal – $M_{2}^{'} (v s_{i})$ ,

– 1st central moment of the signal – $C_{1} (v s_{i})$ ,

– 2nd central moment of the signal – $C_{2} (v s_{i})$ ,

– normalized 1st moment – $N_{1} (v s_{i})$ ,

– normalized 2nd moment – $N_{2} (v s_{i})$ ,

– normalized 1st central moment – $N C_{1} (v s_{i})$ ,

– normalized 2nd central moment – $N C_{2} (v s_{i})$ ,

– abscissa of signal square gravity centre – $G (v s_{i})$ ,

– the signal square variance – $V (v s_{i})$ ,

– the signal equivalent diameter – $E D (v s_{i})$ ,

– the distance from reference signal – $δ (v s_{i})$ .

Next, the mean value and standard deviation of each characteristic were calculated separately for rotating systems remaining in ability state and in inability state in division by measurement point, vibration axis and rotation speed (PAS).

According to the theory presented in [4] characteristics were analysed from the concentration and unequivocally point of view using Eqs. (1) and (2):

1

C H_{C O N} \Leftrightarrow σ_{C H (P A S)} \leq \bar{C H} (P A S) \cdot 0.15,

where $C H_{C O N}$ – concentrated characteristic; $σ_{C H (P A S)}$ – standard deviation of characteristic $C H$ ; $\bar{C H} (P A S)$ – the mean value of the characteristic $C H$ calculated over $P A S$ group of signals.

2

C H_{U E} \Leftrightarrow \{C H \in C H_{C O N} \land σ_{C H (P A S_{I A})} + σ_{C H (P A S_{A})} \leq |\bar{C H} (P A S_{I A}) - \bar{C H} (P A S_{A})|\},

where: $C H_{U E}$ – unequivocal characteristic; $C H_{C O N}$ – concentrated characteristic; $σ_{C H (P A S)}$ – standard deviation of characteristic $C H$ calculated over $P A S$ group of signals; $\bar{C H} (P A S)$ – the mean value of the characteristic $C H$ calculated over $P A S$ group of signals; $P A S_{I A}$ – $P A S$ group of signals registered on rotating systems remaining in inability state; $P A S_{A}$ – $P A S$ group of signals registered on rotating systems remaining in ability state.

On this basis the characteristics which fulfil the concentration and unequivocal conditions to relatively high degree, were chosen for further analysis. These were the signal energy, the signal mean power, abscissa of the signal square gravity centre, the signal square variance and the signal distance from the reference one.

The characteristics formulated five-dimensional space in which fuzzy clusters of ability and inability states can be identified. In the space each signal was expressed as a point and their coordinates were the values of the characteristics. The identification procedure was carried out according to the FCM (fuzzy c-means) method [10].

Analysing the signals, it was stated that ranges of the dimensions of the formulated space were completely different. Therefore, prior to the clusters identification the minimum and maximum values for each dimension was calculated and on this basis the values of each signal coordinates were normalized according to Eq. (3):

3

C H (v s_{i}) = \frac{C H (v s_{i})}{\max (C H (V S)) - \min (C H (V S))},

where: $V S$ – the set of vibration signals, $m a x (C H (V S))$ – maximum value of the characteristic $C H$ of $V S$ ; $m i n (C H (V S))$ – minimum value of the characteristic $C H$ of $V S$ ; $v s_{i}$ – vibration signal no. $i$ , $C H (v s_{i})$ – value of the characteristic $C H$ of the vibration signal $v s_{i}$ .

Thanks to the coordinates normalisation each dimension of the space has the same influence on the clusterisation process.

In the first step of identification the signals were divided into two groups. One group consisted of 17792 signals recorded on rotating systems remaining in ability state and 2854 signals recorded on rotating systems remaining in inability state. This group, called learning set, was used in centre of clusters identification process. The second group consisted of 5670 signals recorded on rotating systems remaining in ability state and 1386 signals recorded on rotating systems remaining in inability state. This group, called testing set, was used for testing the misalignment identification method reliability. Next, the FCM method was applied with the learning set.

The method requires three starting parameters to be set. The number of clusters ( $c n o$ ), coefficient which define the fuzziness of the cluster membership function ( $q$ ) and $ε$ – procedure stop condition. At the beginning, the initial random values of the membership matrix are generated. The matrix defines the membership degree of each signal to each cluster. Later, the clusters centre points were recalculated according to for Eq. (4):

4

c_{j, k} = \frac{\sum_{i = 1}^{I} μ_{i, j}^{q} \cdot C H_{k} (v s_{i})}{\sum_{i = 1}^{I} μ_{i, j}^{q}},

where: $c_{j, k}$ – the dimension no. $k$ of the cluster no. $j$ , $I$ – the amount of the signals $v s$ , $μ_{i, j}$ – the membership degree of the signal no. $i$ to the cluster no. $j$ , $q$ – the coefficient which define the fuzziness of the cluster membership function, $C H_{k} (v s_{i})$ – value of the characteristic $C H$ no. $k$ of the vibration signal $v s_{i}$ , $v s_{i}$ – vibration signal no. $i$ .

Then the membership matrix values are recalculated according to Eq. (5):

5

μ_{i, j s}^{} = \frac{1}{\sum_{j = 1}^{c n o} {[\frac{\sqrt{\sum_{k = 1}^{c h n o} (c_{j s, k} - C H_{k} (v s_{i}))^{2}}}{\sqrt{\sum_{k = 1}^{c h n o} (c_{j, k} - C H_{k} (v s_{i}))^{2}}}]}^{\frac{2}{q - 1}}},

where: $μ_{i, j s}$ – the membership degree of the signal no. $i$ to the cluster no. $j s$ , $c n o$ – amount of clusters, $c h n o$ – amount of dimensions, $q$ – the coefficient which define the fuzziness of the cluster membership function, $c_{j, k}$ – the dimension no. $k$ of the cluster no. $j$ , $C H_{k} (v s_{i})$ – value of the characteristic $C H$ no. $k$ of the vibration signal $v s_{i}$ , $v s_{i}$ – vibration signal no. $i$ .

The calculations are performed until the stop condition, described by the Eq. (6), is met:

6

\sum_{j = 1}^{c n o} \sum_{i = 1}^{v s n o} |μ_{i, j}^{} (t) - μ_{i, j}^{} (t - 1)| < ε,

where: $μ_{i, j}$ – the membership degree of the signal no. $i$ to the cluster no. $j$ , $c n o$ – amount of clusters, $v s n o$ – amount of signals $v s$ , $t$ – iteration number, $ε$ – procedure stop condition.

During the research the clusters identification was performed several times using different values of the cluster number and the membership function fuzziness coefficient in order to account for the signals noisiness and the signals arrangement in the space. The parameters of the identification process were accepted only if the cardinality of the weakest cluster found was not lower than 30 % of the strongest cluster cardinality Eq. (7):

7

\forall c_{j} \in C_{i d e n t} : c a r d (c_{j}) \geq 0.3 \cdot m a x (c a r d (C_{i d e n t})),

c a r d (c_{j}) = \sum_{i = 1}^{I} μ_{i, j},

where: $c_{j}$ – the cluster no. $j$ , $C_{i d e n t}$ – the set of the identified clusters, $c a r d (c_{j})$ – cardinality of the cluster no. $j$ , $m a x (c a r d (C_{i d e n t}))$ – maximum value of the cardinality among the clusters from set $C_{i d e n t}$ , $I$ – the amount of the signals $v s$ , $μ_{i, j}$ – the membership degree of the signal no. $i$ to the cluster no. $j$ .

Table 1The parameters of the clusters identification process

	Number of clusters $c l n o$	The coefficient which define the fuzziness of the cluster membership function $q$	Procedure stop condition $ε$
Signals collected on rotating systems remaining in ability state	2	3	0.1
Signals collected on rotating systems remaining in inability state	2	2	0.1

Table 2The cardinality and coordinates of the identified clusters

	Cardinality	Coordinates
Signals collected on rotating systems remaining in ability state	8084,79	0,068; 0,067; 0,479; 0,451; 0,934
	9707,2	0,025; 0,024; 0,481; 0,457; 0,956
Signals collected on rotating systems remaining in inability state	1026,88	0,176; 0,175; 0,484; 0,459; 0,979
	1827,12	0,078; 0,077; 0,482; 0,458; 0,983

As a result, the value of identification process parameters where chosen (Table 1). Using them, the clusters for learning set were identified. The cardinality and the coordinates of the clusters are presented in Table 2.

For each subgroup, according to the FCM method, one cluster was identified. Next, the obtained clusters were analysed. For each of them, the distance from the clusters identified for the learning set of signals collected on rotating systems remaining in ability state and the distance from the clusters identified for the learning set of signals collected on rotating systems remaining in inability state were calculated according to Eq. (8) and (9):

8

d_{a s} (c t s t_{l s}) = \sum_{j = 1}^{c a n o} \sqrt{\sum_{k = 1}^{c h n o} (c a_{j, k} - c t s t_{l s, k})^{2}},

where: $c t s t_{l s}$ – the cluster identified for the test set no. $l s$ , $d a s$ – distance from clusters identified for signals collected on rotating systems remaining in ability state, $c a n o$ – amount of clusters identified for signals collected on rotating systems remaining in ability state, $c h n o$ – amount of dimensions, $c a_{j, k}$ – the dimension no. $k$ of the ability state cluster no. j, $c t s t_{j, k}$ – the dimension no. $k$ of the test cluster no. $j$ :

9

d_{i n s} (c t s t_{l s}) = \sum_{j = 1}^{c i n n o} \sqrt{\sum_{k = 1}^{c h n o} (c i n_{j, k} - c t s t_{l s, k})^{2}},

where: $c t s t_{l s}$ – the cluster identified for testing set no. $l s$ , $d i n s$ – distance from clusters identified for signals collected on rotating systems remaining in inability state, $c i n n o$ – amount of clusters identified for signals collected on rotating systems remaining in inability state, $c h n o$ – amount of dimensions, $c i n_{j, k}$ – the dimension no. $k$ of the inability state cluster no. $j$ , $c t s t_{j, k}$ – the dimension no. $k$ of the test cluster no. $j$ .

If the distance between specified test cluster and the ability state clusters was bigger than the distance between the cluster and the inability state clusters, then the signals of the analysed test cluster were interpreted as signals of the rotating system remaining in inability state Eq. (10). Otherwise, they were interpreted as signals of the rotating system remaining in ability state:

10

d_{a s} (c t s t_{l s}) > d_{i n s} (c t s t_{l s}),

where: $c t s t_{l s}$ – the cluster identified for testing set no. $l s$ , $d_{i n s}$ – distance from clusters identified for signals collected on rotating systems remaining in inability state, $d_{a s}$ – distance from clusters identified for signals collected on rotating systems remaining in ability state.

In Table 3 the results of the testing set analysis are presented.

Table 3The results of the testing set analysis

	Number of signals tested	Number of signals identified as signals collected on rotating systems remaining in ability state	Number of signals identified as signals collected on rotating systems remaining in inability state
Signals collected on rotating systems remaining in ability state	5670	4536 (80 %)	1134 (20 %)
Signals collected on rotating systems remaining in inability state	1386	126 (9 %)	1260 (91 %)

Presented results indicates that the efficiency of the proposed method is higher in case of the signals collected on the rotating systems remaining in the inability state than in case of signals collected on rotating systems remaining in ability state. Comparing the efficiency of the method with the earlier works of the authors [4] in can be noticed that the efficiency of the inability state identification is slightly lower but still relatively high (higher than 90 %) and the efficiency of the ability state identification is much higher. The overall efficiency of the rotating systems misalignment identification is higher – over 82 %.

4. Conclusions

The conclusions of the carried out research can be expressed in the following list:

– the signal energy, the signal mean power, abscissa of the signal square gravity centre, the signal square variance and the signal distance from the reference one can be used to formulate the method of the rotating systems misalignment identification,

– the efficiency of the presented identification method is – 80 % for signals collected on rotating systems remaining in the ability state and 91 % for signals collected on rotating systems remaining in the inability state,

– the total efficiency of the proposed method is higher than of other methods used in previous works of the authors.

It should be emphasized that the presented method can a universal tool to be used for analysis of any rotating system vibrations.

References

Girtler J. The semi-Markov model of the process of appearance of sea-going ship propulsion system ability and inability states in application to determining the reliability of these systems. Polish Maritime Research, Vol. 4, Issue 80, 2013, p. 18-24.

Search CrossRef
Grządziela Muślewski A. Ł. High quality simulation of the effects of underwater detonation impact. Journal of Vibroengineering, Vol. 15, Issue 1, 2013, p. 106-113.

Search CrossRef
Gurr C., Rulfs H. Influence of transient operating conditions on propeller shaft bearings. Journal of Marine Engineering and Technology, Vol. 12, 2008, p. 3-7.

Search CrossRef
Grządziela A., Musiał J., Muślewski Ł., Pająk M. A method for identification of non-coaxiality in engine shaft lines of a selected type of naval ships. Polish Maritime Research, Vol. 1, Issues 85-2015, 22, p. 65-71.

Search CrossRef
Bielawski P. Diagnostics of marine propeller shafts. Journal of Polish CIMAC, Vol. 6, Issue 2, 2011, p. 31-40.

Search CrossRef
Muślewski Ł., Pająk M., Grządziela A., Musiał J. Analysis of vibration time histories in the time domain for propulsion systems of minesweepers. Journal of Vibroengineering, Vol. 17, Issue 3, 2015, p. 1309-1316.

Search CrossRef
Izydorczyk J., Pionka G., Tyma G. Theory of Signals. Introduction. HELION, Gliwice, 2006.

Search CrossRef
Stefański T. Teoria Sterowania. Tom II. Układy Dyskretne, Nieliniowe, Procesy Stochastyczne Oraz Optymalizacja Statyczna I Dynamiczna. Kielce University of Technology, Kielce, 2001.

Search CrossRef
Szabatin J. Basics of Signal Theory. WKŁ, Warszawa, 2007.

Search CrossRef
Piegat A. Fuzzy Modelling and Control. EXIT, Warszawa, 1999.

Search CrossRef

About this article

Received

01 July 2016

Accepted

04 July 2016

Published

31 August 2016

SUBJECTS

Fault diagnosis based on vibration signal analysis

Keywords

rotating system

misalignment identification

vibration analysis

fuzzy clustering

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.