Condition monitoring is expected to assist in the prevention of machine failures and enhance the reliability with lower maintenance cost. With advancement in sensor and computer technologies, it is now possible to acquire and process a large amount of machine response data in order to extract the characteristic features which provide the indication about health condition of the system. In real applications, the machine has infinite node positions but sensors can be only placed at a finite number of locations. Hence selection of optimal node positions is a challenging task and needs to be addressed. The objective of sound sensor placement optimization is to obtain a sensor layout that gives as much information of the dynamic system as possible for condition monitoring. This paper proposes a methodology for placing sound sensor on a fixed-axis gearbox to obtain high-quality information regarding the dynamic characteristics of machine. The gearbox is operated under no load and load with varying levels of gear faults and sound sensor placement positions. Loudness, loudness level, and sharpness of sound has been considered as response parameters. Mathematical relations and models between input variables and response parameters are developed. Sound sensor placement is optimized for maximum Loudness, loudness level, and sharpness values. The models are experimentally validated and tested. Results indicate that an overall accuracy of 92.2% is achieved and the approach has significant utility in industrial environment where system complexity makes the choice of sensor placement vital for condition monitoring.
optimal sensor placement; response surface methodology (RSM); fault diagnosis; sound analysis; fixed-axis gearbox
Condition monitoring is expected to assist in the prevention of machine failures and enhance the reliability with lower maintenance cost. In real applications, the machine has infinite node positions but sensors can be placed at a finite number of locations. Hence selection of optimal node positions is a challenging task and needs to be addressed. The objective of the current study is to obtain a sensor layout that gives as much information of the dynamic system as possible for condition monitoring. This paper proposes a methodology for placing sound sensor on a fixed-axis gearbox to obtain high-quality information regarding the dynamic characteristics of machine. Mathematical relations and models between input variables and response parameters are developed. Results indicate that an overall accuracy of 92.2% is achieved and the approach has significant utility in industrial environment where system complexity makes the choice of sensor placement vital for condition monitoring.
Intelligent fault diagnosis system is aimed to predict the current status of the machines, which facilitate timely preventive steps to increase the reliability with lower maintenance cost. With advancement in technology, an increase in system complexity leads to increase in probability of system failure. Over more than three decades, many condition monitoring techniques have been developed which are based on measurement of dynamic responses like acoustics, vibration, eddy current, thermal fields, and radiography, acquired using various sensors which are mounted near the vicinity of the machines. Stability and accuracy of the fault diagnosis system are entirely dependent on the quality of information captured by sensors.
Gearbox plays most essential roles in the modern machinery. Most machine breakdowns, related to gearbox, are as a result of improper operating conditions and loading, hence leads to failure of the whole mechanism (Li, Yan, Yuan, Peng, & Li, [
Significant research has recently been conducted in speech recognition, biomedical, and ocean study applications considering sound as a parameter (Bourlard & Morgan, [
In a typical industrial environment, sound signals undergo many reflections before being absorbed, hence considering the presence of nearby reflecting surfaces is mandatory for sound-based condition monitoring of a machine. The sound conditions around machines varies with direction and locations, thus sound-based condition monitoring is generally influenced by noise generated from the surrounding environment which further lead to unreliable data. Sound attenuation could lead to inaccurate readings in lower frequency range of signals.
In order to improve the quality of sensor data, high signal-to-noise ratio needs to be maintained, which can be achieved by either appropriate selection of sensor sensitivity or by optimizing the placement of sensors. Generally, more the sensors are placed on the machine, more the information could be obtained. However, the number of sensors is strictly constrained by the high installation and maintenance cost. This issue is related to optimal sensor placement on a machine under investigation and needs to be addressed. Various efforts for optimizing the sensor placement have been reported in the literature such as effective independence based on fishers information matrix (Stephan, [
However, sound sensors optimal placement study for fixed axis gearbox is still relatively new. The use of Response Surface Methodology (RSM), a statistical tool, for optimizing machining parameters had showed significant results for better surface finish and tool life (Bhogal, Sindhu, Dhami, & Pabla, [
This paper addressed RSM-based sound sensor placement optimization methodology. The main objective under the problem posed is to identify the maximum possible sound sensor node position away from the gearbox without loosing characteristic information. This paper is organized as follows. Section 2 describes sensor placement as an optimization problem, in which experimental setup and data acquisition, experimental design using RSM, model generation, and analysis of the generated mathematical model are represented. Section presents influence of input variables on response parameters, generation of optimal solution, and model validation. Finally, concluding remarks are stated in section followed by research trends.
The basic idea of sensor placement optimization is to identify the optimal region in the vicinity of a machine that gives as much information of the dynamic system as possible. To realize this objective, the problem can be formulated as an optimization problem with optimizing multiple response parameters related to the system. The response parameters are expressed by a generalized equation, as a function of different input variables, and its mathematical representation is given in Equation (
subject to constraints
where Z = response parameters; = input variables, namely fault, load, and location of sensor, respectively, ‘’ = error of the response measured, lb and ub = lower and upper bounds of respective response parameters, and = sets of real numbers.
Figure 1(a) shows the experimental setup for condition monitoring of two stage gearbox. The setup has been designed to simulate real working conditions of a gearbox. The gearbox was driven by a 3 HP, 3-phase induction motor. The rotational speed of the setup could be varied from 100 to 3,600 rpm using variable speed drive. For external loading, magnetic particle braking system was provided with loading capacity from 0.126 to 24.85 N-m. The gearbox consists of single stage fixed axis shafts mounted with spur gears with speed reduction of 3.44:1. Table 1 summarizes the specifications of the gearbox.
For acoustic signal acquisition, three microphones with frequency range from 5 Hz to 20 kHz were used. The faults were induced on the pinion mounted on the input shaft. In this paper, depth wise damage has been simulated on the spur gear tooth by grinding the tooth in steps based on percent removal of teeth depth. A total of three conditions of the gear have been investigated: healthy gear and gear with two stages of depth wise tooth removal i.e. 30 and 50% across the tooth width as shown in Figure 1(b). For all simulated gear conditions, sound signals were acquired. The sampling frequency of the data acquisition system was 12.8 kHz and 30,000 data points were collected for each case at 2,100 rpm and each experiment is repeated five times in order to obtain the average value.
Gearbox specifications
Specification Values Length 27.5 cm Width 19 cm Height 26.5 cm Number of teeth on driver 29 Number of teeth on driven 100 Spur gear pressure angle (ap)
Gear fault (F), load (L), and location of sound sensor (x) was considered as input variables for the experimental setup to investigate the effect on response parameters viz. sound loudness , loudness level , and sharpness and obtain the optimized sensor location. The positions of sound sensors were varied vertically and horizontally at different sensor node positions as shown in Figure 1(c) and (d), respectively. The notations of all running condition is listed in Table 2. Three levels of gear faults were considered i.e. H, CTL30, and CTL50. Experiments were designed for two levels of load i.e. no load and 10.5 N-m load. RSM was utilized to design the experiments for different sensors considering three input variables with different levels as shown in Table 3. The experimental design matrix for the experimentation is shown in Table 4. The results were further analyzed using analysis of variance (ANOVA) and optimized for maximum , and simultaneously. The Design Expert© 8.0 software package was utilized for analysis and optimization of input variables for response parameters.
Fault description
Fault description Notations Healthy H Chipped tooth 30% CTL30 Chipped tooth 50% CTL50
Levels of input variables selected for experimental design
Sensor Level Lowest Center Highest Coding 0 1 Sensor 1 Fault H CTL30 CTL50 Location (position) 15 cm 30 cm 45 cm Load (N-m) 0 - 10.5 Sensor 2 Fault H CTL30 CTL50 Location (position) 50 cm 100 cm 150 cm Load (N-m) 0 - 10.5 Sensor 3 Fault H CTL30 CTL50 Location (position) 50 cm 100 cm 150 cm Load (N-m) 0 - 10.5
This section discusses the results obtained from the experiments conducted for the optimization of sensor location using RSM-based methodology.
Experimental design matrix
Runs Coded parameters A B C 1 2 0 3 1 4 0 5 0 0 6 1 0 7 1 8 0 1 9 1 1 10 1 11 0 0 1 12 1 1 13 0 1 14 0 0 1 15 1 0 1 16 1 1 17 0 1 1 18 1 1 1
RSM-based full factorial 18 experiments for each sensors were conducted. Table 5 shows the regression analysis of the responses of each sensors with standard deviation, adjusted , and predicted values. The adeq precision (press) measures the signal-to-noise ratio. A ratio greater than four is desirable. Obtained press values are greater than required values hence confirms an adequate signal. These results indicate that the model is significant and provides good predictions of the outcomes. Further analysis using ANOVA for 2F1 (two-factor interaction) multi response model was performed to investigate the significance of F-probability value for each sensor model as shown in Table 6.
Accuracy of the RSM for each sensor
Sensor no. Responses Standard deviation Adjusted Predicted Press Sensor 1 0.9847 0.9764 0.9568 34.627 0.9842 0.9755 0.9547 34.068 0.013 0.9374 0.924 0.8917 25.02 Sensor 2 0.9411 0.909 0.8637 18.35 0.9367 0.9022 0.8504 17.889 0.011 0.9656 0.9468 0.9127 22.589 Sensor 3 0.9016 0.8479 0.7891 13.565 0.8967 0.8403 0.7782 13.319 0.011 0.9606 0.9391 0.8988 19.548
ANOVA of 2FI model of multi response parameters
Sensor no. Source SS df MS F-value Sensor 1 Model 0.00568 0.035 6 6 3 0.000947 0.012 118.26 114 69.85 0.018 1 1 1 0.018 532.15 513.09 106.32 1 1 1 21.73 21.39 3.93 0.0007 0.0007 0.0675 0.017 1 1 1 0.017 140.3 135.97 99.3 - 1 1 - - 2.9 1.68 - 0.1168 0.2212 - - 1 1 - - 1.93 0.47 - 0.1919 0.5056 - - 1 1 - - 10.54 11.4 - 0.0078 0.0062 - Sensor 2 Model 0.000625 0.035 6 6 6 0.000104 29.29 27.13 51.46 (x) 1 1 1 82.07 75.37 45.58 (F) 1 1 1 2.89 4.55 9.19 0.117 0.0564 0.0114 (L) 0.025 1 1 1 0.025 62.81 57.2 221.11 1 1 1 0.27 0.15 0.012 0.6154 0.7051 0.9162 1 1 1 3.44 1.11 1.05 0.0907 0.3.147 0.3274 1 1 1 24.23 24.41 31.79 0.0005 0.0004 0.0002 Sensor 3 Model 0.000644 0.034 6 6 6 0.000107 16.8 15.91 44.66 (x) 1 1 1 47.57 46.37 38.25 (F) 1 1 1 0.63 1.2 .121 0.4436 0.296 0.2943 (L) 0.026 1 1 1 0.026 38.18 35.23 203.74 1 1 1 0.37 0.32 1.42 0.5571 0.5834 0.2586 1 1 1 2.14 0.64 0.95 0.1719 0.4396 0.3519 1 1 1 11.91 11.71 22.41 0.0054 0.0057 0.0006
The model shows an acceptable value of F-probability i.e. less than 0.5 for all the sensor models, thus suggesting all the three models for different sensors to be significant. Consequently, the response equations for all the three responses for each sensor were developed as represented in Equations (
Sensor 1:
Sensor 2:
Sensor 3:
The multi-response parameters were analyzed for simultaneous maximization to obtain relevant data with no loss of characteristic information, utilizing the error values obtained in actual value and the predicted value from the empirical model. Table 7 shows a maximum prediction error of 9.14% in case of , 2.53% for and 8.57% being the maximum in case of predicting value. Thus, it has been found that the model is efficient in prediction of , followed by and least for , with overall confidence level above 95%. Moreover, the confirmation of data generated from residual curves for sensor 1 is shown in Figure 2, which represents the interaction between actual and predicted data. The proximity of all the data points to the inclined line indicates the validity of model and proves its adequacy. Predicted versus actual curve of responses for sensor 2 and sensor 3 shows the similar results (see Appendix 1).
Experimental design matrix for individual sensor with actual, predicted and percent error for all responses
Sensor no. Runs Coded parameters Loudness, (sone) Loudness level, (phone) Sharpness, (acum) A B C Actual Predicted % Error Actual Predicted % Error Actual Predicted % Error Sensor 1 1 96.64 97.73 105.95 106.30 4.42 4.43 2 0 66.81 67.82 100.62 100.92 3.92 3.79 3.56 3 1 47.99 45.72 4.72 95.85 96.05 3.22 3.30 4 0 105.75 103.62 2.01 107.25 106.23 0.95 5.01 4.58 8.57 5 0 0 71.86 70.55 1.82 101.67 101.18 0.49 4.11 3.89 5.19 6 1 0 51.48 50.92 1.09 96.86 96.58 0.29 3.33 3.39 7 1 94.27 95.19 105.59 106.15 4.20 4.34 8 0 1 68.98 67.42 2.27 101.08 101.43 3.95 4.01 9 1 1 52.44 50.14 4.39 97.13 97.12 0.01 3.51 3.47 0.98 10 1 114.86 109.98 4.25 108.44 108.85 5.80 6.05 11 0 1 89.45 84.32 5.73 104.83 103.50 1.26 5.08 4.91 3.30 12 1 1 56.69 55.92 1.35 98.25 98.66 3.95 4.13 13 0 1 123.90 113.54 8.36 109.53 110.40 6.20 6.34 14 0 0 1 91.83 92.64 105.21 105.26 5.34 5.10 4.55 15 1 0 1 65.93 64.02 2.90 100.43 100.58 4.05 4.26 16 1 1 158.03 162.33 113.04 111.99 0.93 6.91 6.65 3.82 17 0 1 1 98.83 99.51 106.27 107.07 5.13 5.30 18 1 1 1 78.36 79.10 102.92 102.57 0.34 4.79 4.40 8.11 Sensor 2 1 44.54 45.69 94.77 94.92 2.99 3.02 2 0 32.53 35.50 90.24 91.65 2.73 2.82 3 1 29.18 29.02 0.54 88.67 88.60 0.08 2.66 2.64 0.69 4 0 47.13 43.69 7.28 95.58 93.85 1.82 3.15 2.95 6.31 5 0 0 34.77 33.42 3.89 91.20 90.79 0.45 2.78 2.76 0.97 6 1 0 29.51 27.88 5.52 88.83 87.93 1.02 2.59 2.59 0.13 7 1 37.82 38.34 92.41 92.79 2.79 2.88 8 0 1 29.46 31.57 88.81 89.95 2.65 2.70 9 1 1 26.59 26.83 87.33 87.27 0.07 2.56 2.54 1.00 10 1 46.15 45.01 2.47 95.28 95.04 0.25 3.36 3.36 0.05 11 0 1 37.51 37.76 92.29 92.42 3.23 3.17 1.65 12 1 1 34.06 32.52 4.54 90.90 89.94 1.06 3.07 3.00 2.29 13 0 1 54.35 52.75 2.94 97.64 97.29 0.36 3.65 3.69 14 0 0 1 39.81 42.71 93.15 94.70 3.30 3.46 15 1 0 1 34.90 37.31 91.25 92.25 3.17 3.27 16 1 1 64.16 63.71 0.70 100.04 99.64 0.39 4.17 4.08 2.26 17 0 1 1 51.84 51.88 96.96 97.10 3.85 3.81 1.00 18 1 1 1 47.22 43.76 7.33 95.61 94.68 0.98 3.64 3.58 1.74 Sensor 3 1 47.16 48.69 95.59 95.87 3.14 3.12 0.85 2 0 33.09 34.79 90.48 92.20 2.76 2.87 3 1 29.49 29.57 88.82 88.81 0.01 2.65 2.65 0.04 4 0 51.82 47.44 8.44 96.95 94.50 2.53 3.13 2.96 5.19 5 0 0 34.52 34.34 0.51 91.09 91.23 2.74 2.77 6 1 0 31.76 29.39 7.45 89.89 88.17 1.91 2.66 2.60 2.06 7 1 39.38 39.22 0.41 92.99 93.18 2.79 2.82 8 0 1 29.00 31.20 88.58 90.27 2.58 2.68 9 1 1 26.95 27.31 87.52 87.54 2.61 2.56 2.14 10 1 52.67 49.15 6.68 97.19 96.44 0.77 3.56 3.55 0.46 11 0 1 37.04 40.23 92.11 93.43 3.23 3.29 12 1 1 35.69 34.06 4.58 91.57 90.60 1.07 3.19 3.07 3.60 13 0 1 61.09 56.21 7.98 99.33 98.25 1.08 3.78 3.73 1.31 14 0 0 1 40.68 42.96 93.46 95.44 3.34 3.51 15 1 0 1 39.31 38.88 1.11 92.97 92.78 0.20 3.28 3.31 16 1 1 66.33 65.66 1.02 100.52 100.13 0.38 4.07 3.94 3.10 17 0 1 1 51.13 53.60 96.76 97.54 3.70 3.76 18 1 1 1 48.13 45.29 5.91 95.89 95.08 0.84 3.64 3.60 1.28
The mathematical models developed for each sensor using RSM were analyzed and were found significant for all response parameters. The effects of various input parameters on individual response parameters were studied to analyze the level of influence and are listed in Table 8. The results have been represented graphically and comparison charts were generated accordingly. Finally, the desired conditions for optimization were simultaneously evaluated for predicting the optimal results.
The input variables, that is, location, fault, and load, were analyzed to study their individual effect on , , and . The results were represented in the form of multidimensional curves representing varied slopes, thus indicating variation in level of influence on individual response variables.
The loudness of sound signal is a perceptual measure of the effect of the energy content of sound, which is directly connected to the change in physical system (Arnold, [
The increase in load is the next most influential parameter in increasing value followed by fault. The behavior of is seen to be almost constant with an increase in fault when there was no load on the system as shown in Figure 3(a). However, an increasing trend in value is observed with an increase in fault under loading condition. shows a decreasing trend with an increase in location. Table 8 lists the values of at different points shown in Figure 3(a) and (b). An increase of 20.5 and 51.03% in value of for point A and B, respectively, was observed which demonstrate the effectiveness of load for fault identification. For sensor 2 and sensor 3, similar trending behavior of was observed (see Appendix 2 and 3). However, the maximum amplitude of for sensor 2 and sensor 3 decreased as compared to sensor 1. The reason behind this is the near field position of sensor 1 as compared to far-field position of sensor 2 and sensor 3.
Effect of input parameters on response variables for sensor 1
Response parameter Point Status Value (sone) A minimum L and F, maximum x 48.67 B minimum L, maximum F and x 97.89 A’ minimum F, maximum L and x 58.65 B’ maximum F, L and x 147.85 (phone) C minimum L and F, maximum x 96.05 D minimum L, maximum F and x 106.15 C’ minimum F, maximum L and x 98.66 D’ maximum F, L and x 111.99 (acum) E minimum L and F, maximum x 3.30 F minimum L, maximum F and x 4.73 E’ minimum F, maximum L and x 4.132 F’ maximum F, L and x 6.64
Loudness level of sound signals was measured to depict the sound amplitude in a particular frequency band of interest, which shows change in gear-mesh frequency amplitude. Hence any change in gear-mesh frequency and its harmonics directly influence the loudness level. Figure 4(a) and (b) shows the effect of location and fault for maximum and minimum load on for sensor 1. Figure 4(c) shows that sensor location exhibits maximum influence on followed by load and fault. Table 8 list the values of at different points shown in Figure 4(a)-(b). An increase of 2.71 and 5.50% in value of for point C and D, respectively, was observed which indicate enhancement of fault frequency identification under load using loudness level. However, a linear increasing trend of value with increasing fault was observed for both minimum and maximum loads. Similar trend of was observed for sensor 2 and senor 3 (see Appendix C2 and 3).
The sharpness of sound signal reflects the ratio of high frequency level to overall level. Sounds that exhibit higher energy in high frequencies tend to be sharper. It is a crucial parameter which is physically detected by human ears. In general, for any change in physical system properties, there is a noticeable change in sharpness of sound. From perturbation curve shown in Figure 5(c), the sensor location is the most influential factor with negative slope for as sensor location transition from near-field position to free-field position thereby dissipating higher energy component to the surroundings. However, load tends to increase in the presence of fault as faulty teeth mesh with healthy teeth with more impact thereby generating high-energy frequency component as shown in Figure 5(b). Table 8 list the values of at different points shown in Figure 5(a)-(b). An increase of 25.21 and 40.38% in value of for point E and F, respectively, was observed indicating load effectiveness with increasing fault. Similar trend of was observed for sensor 2 and senor 3 (see Appendix 2).
The effect of all the input factors on individual responses was modeled. However, to obtain the maximum information about the condition of gearbox, all these response parameters needed to be maximized based on input variables. Also, it is necessary to find the maximum possible sensor node position away from the gearbox that can be achieved without compromising the information loss. Considering that an optimal solution was developed for maximum , maximum , and maximum . For this, varied conditions of input variables were modeled and tested. The models which generated an overall highest desirability for all the responses simultaneously were with maximum sensor location, in-range fault, and in-range load. The optimal values thus obtained for sensor 1 are 24-cm vertical location, CTL50 fault, and 10.5 N-m load. For sensor 2 and sensor 3, the values are 112.5-cm horizontal location, CTL50 fault, and 10.5 N-m load.
The set of input parameters for individual sensors were computed for the optimal model generation and were experimentally validated. For the individual sensors, three repetitions of experiments were performed and the average value for each response parameters was considered. The results were compared with those predicted by the RSM model. The graphical representation of optimal solutions obtained experimentally, compared with predicted responses is shown in Figure 6. The results showed a maximum of 7.8% error in case of R1. So, an overall accuracy of 92.2% was achieved.
Sensor placement is of key importance in condition monitoring of gearbox. In the present work, RSM-based sound sensor location optimization scheme is presented. The optimization accounts for maximum sensor location, both in horizontal and vertical directions, to maximize the loudness, loudness level, and sharpness considering load on the gearbox and fault in gear as input variables. The approach is designed such that maximum dynamic information of the system can be retained for the fault diagnosis of gearbox. The use of statistical modeling technique for predicting multi-process factors to find the optimal sensor location is an effective technique. All response parameters viz. loudness, loudness level, and sharpness are influenced mostly by sensor location, which indicates critical information loss as distance increases. Load on the system enhances the fault information as there is an increase in response parameters at maximum load conditions for all considered ranges of faults. The optimal model for maximum loudness, maximum loudness level, and maximum sharpness was envisaged, and out of various possible models, the one with maximum sensor location with in-range load and in-range fault was selected. The responses proposed by the optimal model were at 24-cm vertical sensor location and 112.5-cm horizontal location with maximum fault and maximum load. Finally, the authenticity of model was confirmed experimentally and a maximum error of 7.8% was obtained, which suggested an acceptable liaison with the predicted model. Thus, the utilization of RSM for finding optimal sensor location is recommended.
Future research may focus on testing of proposed methodology for various types of faults (e.g. root crack) to demonstrate the generalization of the sensor placement optimization using RSM in the reliable classification of various faults. Authors also suggest to evaluate the robustness of the methodology using field data collected in real-world applications where other rotating machines are working near the vicinity of the testing machine. The proposed method can be employed in non stationary operation like robot condition monitoring.
Source: Author.
PHOTO (COLOR): Figure 1. Experimental setup with horizontal, vertical sound sensor node positions and different gear conditions.
PHOTO (COLOR): Figure 2. Predicted vs. actual curve for all three responses for sensor 1 (a) Loudness (b) Loudness level (c) Sharpness.
PHOTO (COLOR): Figure 3. Effect of input parameters on for sensor 1.
PHOTO (COLOR): Figure 4. Effect of input parameters on for sensor 1.
PHOTO (COLOR): Figure 5. Effect of input parameters on for sensor 1.
PHOTO (COLOR): Figure 6. Comparison chart of actual versus predicted responses for optimal conditions for (a) Sensor 1 (b) Sensor 2 (c) Sensor 3.
PHOTO (COLOR): Figure A1. Sensor 2 predicted versus actual curve for all three responses.
PHOTO (COLOR): Figure A2. Sensor 3 predicted versus actual curve for all three responses.
PHOTO (COLOR): Figure B1. Effect of input parameters on for sensor 2.
PHOTO (COLOR): Figure B2. Effect of input parameters on for sensor 2.
PHOTO (COLOR): Figure B3. Effect of input parameters on for sensor 2.
PHOTO (COLOR): Figure B4. Effect of input parameters on for sensor 3.
PHOTO (COLOR): Figure B5. Effect of input parameters on for sensor 3.
PHOTO (COLOR): Figure B6. Effect of input parameters on for sensor 3.
By Vanraj; S.S. Dhami and B.S. Pabla
Edited By Ding