Mathematical Problems in Engineering

Volume 2014 (2014), Article ID 472315, 10 pages

http://dx.doi.org/10.1155/2014/472315

Research Article

Study on the Method for Collecting Vibration Signals from Mill Shell Based on Measuring the Fill Level of Ball Mill

Peng Huang, Minping Jia, and Binglin Zhong

School of Mechanical Engineering, Southeast University, Nanjing 211189, China

Received 17 April 2014; Revised 26 June 2014; Accepted 30 June 2014; Published 15 July 2014

Academic Editor: Qingsong Xu

Copyright © 2014 Peng Huang 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.

Abstract.

At present the method for measuring the fill level which used the vibration signal of mill shell shows its advantage compared with other methods. However, this method is developed late, and the technique for collecting the vibration signal from mill shell is immature. In this paper, a novel method for collecting the vibration data from mill shell is proposed. Firstly, the layout scheme of vibration sensors on mill shell is given by analyzing the axial and circumferential movement of coal powder in roller. And a special data acquisition system is developed, which can acquire vibration data from different axial and circumferential positions on mill shell. Then the sampling frequency is obtained based on impact model and hierarchical model of steel balls. At the same time, the impact region on mill shell caused by steel balls is considered as the collecting region of vibration signals. Experimental result shows that vibration signals collected by the method proposed in this paper present a high sensitivity to the changes on fill level compared with vibration data of mill bearing, which provides a reliable basis for accurate measurement of the fill level.

1. Introduction

It is well known that the ball mill pulverizing system has advantages of good adaptability for mineral aggregate of different types, and it has been widely employed in power plant, cement industry, chemical industry, metallurgy, and many other fields. However, this system needs high energy intensive. For example, in the power plant, it makes up about 20 percent of all electricity consumption of power plant [1, 2]. So making the mill operation safely and having high efficiency are important goals for enterprises. Because the ball mill system is of nonlinearity, long time delay, and time-varying, the fill level of coal powder cannot be measured accurately, which leads the mill to be operated usually under very uneconomical condition to prevent the coal blockage, and the automatic control system is also hard to accomplish steady operation of the mill. Thus the exact measurement of fill level is a key and basic factor for realizing automatic, reliable, and efficient operation of the mill system.

Amongst current techniques of measuring the fill level, the vibration method that can measure the fill level through extracting the power, energy, and vibrated root mean square (RMS) from vibration signal of the mill bearing housing has been applied widely [3–5]. However, there are many reasons and vibration sources to cause the vibration of bearing housings. In addition to the impact of steel balls in roller, there are still some other vibration sources, such as transmission system vibration, foundation vibration, vibration caused by the rotational asymmetry, and installation error. These vibration sources will affect the time and frequency domain characteristics of vibration signals of bearing housings. Furthermore, vibration transmission path also has a certain impact for vibration signals of bearing housings. So vibration signals of the bearing housings cannot accurately reflect the collision condition of steel balls in the mill and the fill level information.

Because the vibration of mill shell is caused directly by the impact of steel balls in roller, this signal can reflect directly and exactly the fill level information compared with vibration data of bearing housings. However, due to the rotating of mill shell, collecting the vibration signal from mill shell is difficult, and the study on measuring fill level which used this signal is relatively little. Among these, the work of Gugel et al. [6, 7] was noteworthy. In their work, two vibration sensors, 180 degrees apart (shown in Figure 1), were used to collect signal information from the mill shell so that instantaneous readings can be averaged together, and the averaged data was fed into a dynamic neural network and a fill level measurement from 0 to 100 percent was output. But this technique has just started, and the data acquisition method is not mature. In the process of data collection, the layout of vibration sensors on the mill shell, the determination of sampling frequency, and choosing the collection region all need to be studied in depth. And these problems will be discussed and resolved in this paper.

Figure 1: Three different sensor alignment examples for two-sensor based shell system [6]. Note: distance from sensor 1 to point of impact. distance from sensor 2 to point of impact. mill circumference/2.

2. Layout of Vibration Sensors

Because the working situation in mill roller cannot be observed directly, the distribution of coal powder cannot be accurately obtained. At present it is generally thought that the inside of drum is divided by a belt of demarcation [8]. Above the belt of demarcation, there is a flowing space of coal powder airflow. Under the belt of demarcation, there is a motion space of mixed materials, and the mixed materials are composed of coal cinder, coal powder, and steel balls. The motion of mixed materials can be described by the axial and circumferential movements. For the axial motion, the amount of coal at the inlet of mill will be increased firstly along with the increase of the coal inputting, and then the amount of coal at the outlet of mill begins to rise when coal is sent to the outlet of the mill by airflow. From the perspective of longitudinal section of mill, the dividing line is declining from the inlet to the outlet, as shown in Figure 2(a) [9]. The circumferential motion of mixed materials is shown in Figure 2(b) because of the influence of rotation of mill.

(a) The longitudinal section of mill

(b) The cross section of mill

Figure 2: The schematic diagram of coal powder distribution in roller [9].

It is clearly shown in Figure 2(a) that the distribution of coal powder is nonuniform in roller. So information from only a vibration sensor installed on the mill shell cannot describe the amount of coal completely. In this study, according to the distribution characteristic of coal powder, four acceleration sensors are installed at different positions that locate at the same axis on the mill shell, as shown in Figure 3 (sensors 1, 2, 3, and 4). In the conventional vibration method for measuring fill level, two acceleration sensors (sensors 5 and 6, as shown in Figure 3) are installed at forward and back bearing housing, separately.

Figure 3: Schematic diagram of the fixed position of acceleration sensors.

3. Data Acquisition System of the Shell Signal

Sensors that are installed on the mill shell will rotate together with the rotation of mill, which is not the same with the traditional vibration method that sensors can be directly connected with data acquisition computer by cable. So a special data acquisition system was developed in this study. This system is also installed on the mill shell, which was connected with sensors by sensor cables, and can transfer the signals to computer by wireless way. This system consists of host computer and slave computer two-level system. The function of the slave system is data acquisition and transferring the data to computer. The function of host system is receiving, saving, displaying, and analyzing signals. Figure 4 shows the overall structure of this data acquisition system.

Figure 4: Structure of data acquisition system.

Because sensors mounted on the mill shell will rotate together with the mill and the vibration strength differentiates significantly at different circumferential positions on mill shell, sensors need to be controlled to collect signals at the designated circumferential position. This difficult problem is solved by photoelectric switch in this research, and Figure 5 shows its working principle.

Figure 5: The work functional diagram of photoelectric switch.

In Figure 5, a trigger device of photoelectric switch is placed below the mill. When the photoelectric switch passes through this device, a switch signal will be generated, and this trigger position is considered as the initial phase of sensors. And the collecting position on mill shell can be confirmed by the delay time transferred from host system. Figure 6 shows a real photo of the slave system.

Figure 6: Installation image of the slave system mounted on the mill shell.

4. Confirming the Sampling Frequency of the Vibration Signal

In the process of signal acquisition, if the sampling frequency is set too low, the frequency spectrum of signal will produce overlapping phenomenon, which leads to the distortion of signal. So the appropriate sampling frequency is the key to guarantee the accuracy of data acquisition. According to the Nyquist sampling theory, when sampling frequency () is 2 times higher than the highest frequency () of signal, the sampled signal can accurately reflect the information of original signal [10]. However for the impact vibration signal, the highest frequency in the spectrum may be infinite, and the determination of sampling frequency is difficult.

If the high frequency part of signal can be cut, the highest frequency of signal is considered as limited value. At the same time, the loss of signal information is little, and the sampling frequency is not too high in accordance with the sampling theory [11]. For the vibration signal of mill shell, the main scope of frequency response caused by the impact of steel ball needs to be obtained, and the maximum frequency of this scope is considered as the highest frequency of signal to obtain the sampling frequency. In this research, this frequency band is gained by the theoretical derivation, which is combined with the actual vibration signal of mill shell to determine the sampling frequency.

4.1. Model of Steel Ball Impact

In the operational process of mill, the mill shell is impacted continuously by steel balls, and the impact force can be seen as a series of instantaneous force [12]. So in this study, the situation that the single steel ball impacts the mill shell needs to be analyzed simply, and the characteristics of frequency spectrum for the shell vibration will be obtained. Figure 7 shows the radial vibration of shell caused by impact of a single steel ball.

Figure 7: The model of shell vibration caused by a single steel ball.

The action spot on the shell caused by impact of a single ball is , and the instantaneous impact force at time can be described by

where is the impact force of steel ball, is unit impulse function whose variable is angle, and is unit impulse function whose variable is time. can be obtained by the theorem of impulse:

where is the mass of steel ball and is velocity variation of steel ball when this ball runs time.

According to the vibration theory, the radial displacement of the shell meets the Poisson equation [13], which is expressed by

where is the propagation velocity of vibration in the mill shell, which can be calculated as follows:

where    is the density of mill shell (7.8 × 10³ kg/m³) and is Young’s modulus.

In the cylindrical coordinate system (see Figure 7), the following equation is obtained by (1) combined with (3):

At both sides of (5) for integral of , and combined with (2), can be determined by

Taking into account the boundary condition, the whole problem can be represented as

The radial displacement of the shell is obtained by segregation variable method as follows:

where is mode of vibration and is the distance between the steel ball and the axis of mill shell.

The vibration velocity of mill shell is gained by the derivation of (8):

The vibration acceleration of mill shell can be denoted as

It can be seen from (10) that the vibration frequency of mill shell caused by impact of single steel ball in the frequency spectrum is distributed. And its value connects with the propagation velocity of vibration in the mill shell , the distance between the steel ball and the axis of mill shell , and mode of vibration . The discrete frequency caused by the impact of steel ball can be represented as

where is the angular frequency of impact vibration, which can be obtained by (9) and (10).

Practical experience shows that the first 15 modes of vibration can express the vibration condition completely for the structure of high freedom degree [14]. So in (11), can take the first 15 modes of vibration. However, the different steel balls may move at different steel ball layers in roller, and value will be verified with different layers. The hierarchical model of steel balls needs to be built to calculate the value.

4.2. Hierarchical Model of Steel Balls

In roller, steel balls are layered ideally at any cross section, and Figure 8 shows the layered condition of steel balls. In this study, the motion conditions of other layer balls are considered similarly with those of the first layer. The ball at the first layer moves the circular motion along with the rotation of mill firstly and then moves the parabolic motion when this ball is promoted to a certain height. In Figure 9, the radius of circular motion for each layer is denoted by . In the experiment, the radius of mill shell is 1.75 m, and the diameter of steel ball is 0.05 m. Therefore, the radii of circular motion from the first layer to th layer are 1.75 m, 1.70 m,  m, respectively.

Figure 8: Schematic diagram of layered condition of steel balls.

Figure 9: Schematic diagram of motion trajectory of each layer balls in roller.

For steel ball at the first layer, the following equations can be obtained by its motion trajectory of circular and parabolic motions:

where is the rotation speed of mill and is the acceleration of gravity, 9.8 m/s².

Because the motion trajectories of other layers balls are similar to those of the first layer, the following equations can be gained:

The runtime of parabolic motion for the ball of th layer can be calculated as follows:

The time interval that two adjacent balls are thrown from the starting point of parabolic motion (in Figure 9) at the th layer can be represented as

where is the angular velocity of mill. And at the th layer the number of balls on the parabolic segment is calculated from

where means to round numbers. And at the th layer the number of balls on the circular segment is calculated from

where is the arc length of the circular motion, which is obtained by

So the number of balls at the th layer is expressed by

The number of balls layers can be solved by the following equation:

where is the total mass of steel balls, is the mass of single steel ball, and is the lateral length of mill shell and is equal to 6 m.

In the experiment, the total mass of steel balls is 38 tones, and steel balls are distributed evenly along the 120 axis sections of mill shell. Because the motion conditions of balls for every section are almost identical, the hierarchical condition of balls needs to be analyzed for any section simply [15]. Table 1 shows the hierarchical condition of balls at a certain axis section of mill shell.

Table 1: The hierarchical condition of steel balls.

The serial number of balls layer
(10⁻² s)

1
1.55
68
90
158

2
1.60
65
85
150

3
1.65
62
80
142

4
1.75
59
75
134

The impact condition of mill shell caused by steel balls of different layers can be simplified as the impact condition that steel balls impact the shells of different radii. According to (11), the scope of vibration frequency of mill shell caused by impact of balls at th layer is ~ when takes the first 15 modes of vibration. Table 2 shows the impact discrete frequency of each layer ball.

Table 2: Vibration frequency of mill shell caused by impact of each layer balls.

The serial number of balls layer
Frequency at first mode of vibration (Hz)
Frequency at 15th mode of vibration (Hz)
Scope of vibration frequency of mill shell (Hz)

1
460
6900
460~6900

2
474
7110
474~7110

3
488
7320
488~7320

4
504
7560
504~7560

In Table 2, the impact vibration frequency of mill shell caused by a single ball at different layers is distributed in the scope 460~7560 Hz. In practical work process of mill, mill shell is impacted by many layers balls, and every layer has a great quantity of balls. So the impact vibration frequency of mill shell can be regarded as continuous, and its scope is 460~7560 Hz.

4.3. Confirming the Sampling Frequency

According to the above information, the theoretical highest frequency () of vibration signal is 7560 Hz, and the sampling frequency () should be 2 times higher than based on Nyquist sampling theory. However in practical engineering application, the sampling frequency () should be 5~10 times of to describe the signal accurately [16]. In order to choose the appropriate sampling frequency, in experiment the sampling frequencies () were set at 20 KHz, 40 KHz, 50 KHz, and 60 KHz to collect vibration signals, respectively. And the sampling point was set at 4096 points. In experiment, the steel balls load was 38 tones, and coal feed amount was 55 t/h. Figure 10 shows the frequency spectrums of vibration signals at different sampling frequency conditions.

(a) = 20000 Hz

(b) = 40000 Hz

(c) = 50000 Hz

(d) = 60000 Hz

Figure 10: Frequency spectrums of vibration signals at four types of sampling frequencies.

It can be seen from Figure 10(a) that the amplitude is still considerable at the biggest frequency (10000 Hz) when sampling frequency (20 KHz) is twice as much as the theoretical highest frequency (). This shows that this sampling frequency (20 KHz) cannot describe the original vibration information accurately. In Figures 10(b), 10(c), and 10(d), the amplitude of high frequency part is little when the sampling frequency exceeds 5 times, and the signal can describe the original impact information of steel balls entirely. When the sampling demand is met, the sampling frequency should be set at a low value to enhance frequency resolution. Therefore in this study, the sampling frequency is set at 40 KHz.

5. Confirming the Region of Signal Acquisition

As sensors rotate with the rotation of mill, the collecting region of sensors on mill shell needs to be determined. Due to a close relationship between the vibration of mill shell and fill level, in this research the impact region of the first layer steel balls on mill shell is chosen as the region of data acquisition. In the author’s previous study, the charge motion of mill has been simulated by DEM (Discrete Element Method) [17, 18]. And the impact region of the first layer steel balls is shown in Figure 11. In order to reduce the acquisition group and time of vibration signals, the region of data acquisition is in the range of 310°~340° on mill shell, as shown in Figure 11. In other words, only when the sensors rotate in this region will it be in a work state.

Figure 11: The region of data acquisition on mill shell.

According to sampling frequency and sampling point , the time span of collecting a set of data can be obtained by

When a set of data is collected, the angular of mill shell running is calculated by

According to (22), , , and , value is 10.80°. In other words, this vibration data contains vibration information of about 10 degrees arc length of mill shell. In order to obtain the vibration information of collected region (between 310° and 340° on mill shell), sensors need to work at angular positions 310°, 320°, and 330° on mill shell, respectively.

6. Results and Discussion

6.1. Calculating the Sensitivity Index of Reflecting Fill Level Change

The vibration signals of mill shell can reflect information of fill level directly and accurately compared with vibration data of bearing housings, which is demonstrated by the vibrated mechanism of mill shell and bearing housing that is introduced in the above content. However this conclusion is only the result of theoretical analysis, and it needs to be verified by analysis of actual data. In this research, the superior signal for measuring fill level is confirmed by comparing the sensitivity that vibration signals of mill shell and bearing housing reflect the change of fill level. The sensitivity index of reflecting fill level change is calculated by the RMS of these two kinds of vibration signals.

During the working condition of ball mill, the amount of fill level can be controlled by the coal feeder. When the coal feeder is set higher, the fill level is higher, and vice versa. So the different coal feed amount conditions can denote different fill level conditions. In experiment, the vibration signals were collected under eight stable coal feed amount conditions, that is, 0 t/h, 20 t/h, 30 t/h, 35 t/h, 40 t/h, 45 t/h, 50 t/h, and 55 t/h. The sampling frequency of bearing housing signals was set at 19.2 kHz, and the sampling point was set at 2048 points. In the following, the experimental steps and process of collecting data are illustrated.

(1)At normal work state of mill, the coal feeder was stopped, and during this operation only coal powder was output from the outlet and no coal was input from the inlet by coal feeder. This process lasted 10 minutes, and then there was no coal load (0 t/h) in the mill.(2)Under the stable condition of 0 t/h, the shell and the bearing vibration signals were collected simultaneously, and thirty pieces of data were collected at the forward and back bearing housing, separately.(3)During the process of collecting bearings signals, four sensors installed on mill shell collected thirty pieces of data at angular positions 310°, 320°, and 330°, separately.(4)With the same method, the shell and bearing vibration signals were obtained under the 20 t/h, 30 t/h, 35 t/h, 40 t/h, 45 t/h, 50 t/h, and 55 t/h conditions, respectively.

The RMS of forward and back bearing housing signals at every coal feed amount condition is the average value of RMS, which is obtained by

where is the number of sampling point and is equal to 2048; is the discrete vibration value in time domain of th data.

In the process of calculating RMS of mill shell signals, the RMS is expressed by the average value of RMS of three angular position (310°, 320°, and 330°) signals, which is expressed by

where is the number of sampling point and is equal to 4096; is the discrete vibration value in time domain of th data collected at angular position 310°; is the discrete vibration value in time domain of th data collected at angular position 320°; is the discrete vibration value in time domain of th data collected at angular position 330°.

Because signals amplitudes have very big differences between mill shell and bearing housing, these two kinds of signals are carried out by the normalization method. This normalization method is expressed by

where is the data that are obtained by the normalization processing, is the RMS, is the maximum RMS value in every condition at every vibration station, and is the minimum RMS value in every condition at every vibration station.

The sensitivity index of reflecting fill level change is obtained by RMS variety of the adjacent coal feeder amount conditions, which is calculated by

where is the sensitivity index of reflecting fill level change, is obtained by the normalization processing for the RMS of the previous condition, and is obtained by the normalization processing for the RMS of the latter condition.

6.2. Comparison between the Vibration Signals of Mill Shell and Bearing Housings

Table 3 shows the sensitivity index of reflecting fill level change at every vibration station.

Table 3: The sensitivity index of reflecting the fill level change for different sensors.

Adjacent conditions
Forward bearing (sensor 5)
Back bearing (sensor 6)
Sensor 1
Sensor 2
Sensor 3
Sensor 4

0 t/h, 20 t/h
0.17535
0.38583
0.17277
0.075801
0.06955
0.034087

20 t/h, 30 t/h
0.54767
0.10914
0.17
0.2019
0.18071
0.15551

30 t/h, 35 t/h
0.082084
0.12782
0.10083
0.12498
0.12462
0.12549

35 t/h, 40 t/h
0.52685
0.09317
0.11314
0.11106
0.17316
0.18405

40 t/h, 45 t/h
0.13752
0.061331
0.30897
0.32467
0.25355
0.28187

45 t/h, 50 t/h
0.1194
0.67477
0.3774
0.4382
0.40363
0.43218

50 t/h, 55 t/h
0.096971
0.15895
0.52832
0.48956
0.54763
0.53186

Average
0.24084
0.23014
0.25306
0.25231
0.25041
0.24929

It is clearly seen from Table 3 that the sensitivity index of mill shell signals is higher than that of bearing housings which indicates that the vibration signals of mill shell are sensitive to detect the small changes of fill level and can also measure the great changes of fill level. So the vibration signals of mill shell are the superior signals for measuring fill level. At the same time it can be seen from Table 3 that the sensitivity index of forward bearing housing is higher than that of back bearing housing, which has a good agreement with existing research conclusion [19].

6.3. The Correlation between the Vibration of Mill Shell and Fill Level

Figure 12 shows the variation in the RMS amplitude of vibration signals of mill shell with coal feed amount. And RMS amplitudes are calculated by (24). In all sensors it is found that RMS amplitude decreases as coal feed amount (fill level) increases. The higher the coal load in the mill is, the more impact energy of steel balls is absorbed by coal powder, and the lower impact energy reaches the mill shell to cause relatively small vibration, and vice versa. This conclusion has a good agreement with the result of theoretical analysis. It shows that vibration strength of mill shell can reflect accurately the information of fill level.

Figure 12: Variation in the RMS amplitude of vibration signals of mill shell with coal feed amount at various sensors.

At the same time in Figure 12, when coal feed amount increases from 0 t/h (0% fill level), at the same coal feed amount condition, measure point (such as sensor 1, as shown in Figure 3) is more close to the inlet of mill; its RMS amplitude is lower. Measure point (such as sensor 4, as shown in Figure 3) is more close to the outlet of mill, and its RMS amplitude is higher. It means that amount of coal powder near the inlet of mill is more, and amount of coal powder near the outlet of mill is little. This result is consistent with the distributional conclusion of coal powder in Section 2 (as shown is Figure 2(a)).

7. Conclusions

(1)The result that the nonuniform distribution of coal powder in roller is obtained by analyzing the axial and circumferential motion of coal powder, and the scheme that multiple acceleration sensors need to be installed on mill shell to acquire the information of fill level is confirmed.(2)According to the characteristic of shell signals, a special data acquisition system was developed in this study to collect the vibration signals at different circumferential positions the mill shell.(3)Based on impact model and hierarchical model of steel balls, the vibration frequency band of mill shell caused by impact of steel balls is obtained, which is combined with sampling theorem and actual data of mill shell to confirm the sampling frequency.(4)Due to a close relationship existing in the vibration of mill shell and fill level, in this research the impact region of steel balls on mill shell was chosen as the region of data acquisition.(5)The vibration signals of mill shell are compared with those of bearing housing according to the sensitivity of change of fill level, and this comparison reveals that the vibration signals of mill shell present higher sensitivity to changes on fill level. It indicates that the method of data acquisition proposed in this paper is valid.(6)The relation between the vibration of mill shell and fill level shows that vibration strength of mill shell could well reflect the change of fill level, which provides a solid basis for the accurate measurement of fill level.

Because fill level of mill is a very important monitoring parameter and control variable for the optimal control of mill, optimal, automatic, and online control of ball mill system is studied based on fill level in author’s research, which will be discussed in the next work.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors would like to thanks Professor Fei-yun Xu and Professor Jian-zhong Hu for the helpfulness of experiments. This work was funded under a Grant from the National Natural Science Foundation of China (no. 51005047).

References

V. K. Gupta and S. Sharma, “Analysis of ball mill grinding operation using mill power specific kinetic parameters,” Advanced Powder Technology, vol. 78–99, pp. 33–36, 2014.

M. H. Wang, R. Y. Yang, and A. B. Yu, “DEM investigation of energy distribution and particle breakage in tumbling ball mills,” Powder Technology, vol. 223, pp. 83–91, 2012.

G. Si, H. Cao, Y. Zhang, and L. Jia, “Experimental investigation of load behaviour of an industrial scale tumbling mill using noise and vibration signature techniques,” Minerals Engineering, vol. 22, no. 15, pp. 1289–1298, 2009.

J. Tang, L. Zhao, J. Zhou, H. Yue, and T. Chai, “Experimental analysis of wet mill load based on vibration signals of laboratory-scale ball mill shell,” Minerals Engineering, vol. 23, no. 9, pp. 720–730, 2010.

Z.-G. Su, P.-H. Wang, and X.-J. Yu, “Immune genetic algorithm-based adaptive evidential model for estimating unmeasured parameter: estimating levels of coal powder filling in ball mill,” Expert Systems with Applications, vol. 37, no. 7, pp. 5246–5258, 2010.

K. Gugel and R. M. Moon, “Automated mill control using vibration signal processing,” in Proceedings of the IEEE Cement Industry Technical Conference, pp. 17–25, Charleston, SC, USA, April 2007.

K. Gugel, G. Palacios, J. Ramirez, and M. Parra, “Improving ball mill control with modern tools based on digital signal processing (DSP) technology,” in Proceedings of the IEEE-IAS/PCA 2003 Cement Industry Technical Conference, pp. 311–318, Dallas, Tex, USA, May 2003.

W. Heng and J. Minping, “A fuzzy control method for ball mill system based on fill level soft sensor,” in Proceedings of the Chinese Control and Decision Conference (CCDC '09), pp. 5888–5891, Guilin, China, June 2009.

Z.-G. Su, X.-J. Yu, and Z.-Z. Lv, “Judgment of ball-mills maximal coal pulverizing capacity based on gray incidence analysis of shaft vibration,” Themal Power Generationr, no. 11, pp. 29–33, 2006.

A. V. Oppenheim and A. S. Willsky, Signals and Systems, Prentice Hall, New York, NY, USA, 2000.

C.-T. Chen, Digital Signal Processing Spectral Computation and Filter Design, Oxford University Press, Oxford, UK, 2001.

J. Chen, The Technique of Noise Control for the Ball Mill, China Power Press, Beijing, China, 2002.

B.-L. Zhong and R. Huang, Introduction to Machine Fault Diagnosis, China Machine Press, Beijing, China, 2007.

S.-Y. Zhou, L. Yi, and J.-S. Li, “Analysis on fatigue of ball mill based on finite element method,” Mining Processing Equipment, vol. 41, no. 1, pp. 60–63, 2013.

L.-T. Chen, Y.-J. Xue, and Y. Liu, “Modal analysis of oversize ball mill tube based on ANSYS,” Advanced Material Research, vol. 889-890, pp. 62–65, 2014.

X.-Y. J, Q.-K. Sun, and L.-J. Yin, “Application of equipment condition monitoring and fault diagnosis,” Machine Tool & Hydraulics, vol. 39, no. 22, pp. 129–131, 2011.

P. Huang, M.-P. Jia, and B.-L. Zhong, “New method to measure the fill level of the ball mill I—Theoretical analysis and DEM simulation,” Chinese Journal of Mechanical Engineering, vol. 23, no. 4, pp. 460–467, 2010.

P. Huang, M.-P. Jia, and B.-L. Zhong, “Analysis of the relation between charge motion and fill level of the ball mill using the DEM simulation,” in Proceedings of the International Conference on Information Science and Engineering, pp. 1647–1650, Yangzhou, China, September 2011.

H. Wang, M. Jia, and Z. Chen, “Fill-level measuring and control system based on multi-sensor information fusion for ball mill,” Electric Power Automation Equipment, vol. 32, no. 9, pp. 138–141, 2012.