# Sensor Placement Optimization of Vibration Test on Medium-Speed Mill

Authors: Lihua Zhu; Jun Dai; Guoliang Bai
Journal: Shock and Vibration
 (2015)

---

## Abstract

Condition assessment and decision making are important tasks of vibration test on dynamic machines, and the accuracy of dynamic response can be achieved by the sensors placed on the structure reasonably. The common methods and evaluation criteria of optimal sensor placement (OSP) were summarized. In order to test the vibration characteristic of medium-speed mill in the thermal power plants, the optimal placement of 12 candidate measuring points inX, Y, and Z directions on the mill was discussed for different targeted modal shapes, respectively. The OSP of medium-speed mill was conducted using the effective independence method (EfI) and QR decomposition algorithm. The results showed that the order of modal shapes had an important influence on the optimization results. The difference of these two methods on the sensor placement optimization became smaller with the decrease of the number of target modes. The final scheme of OSP was determined based on the optimal results and the actual test requirements. The field test results were basically consistent with the finite element analysis results, which indicated the sensor placement optimization for vibration test on the medium-speed mill was feasible.


---

## Body

## 1. Introduction

The model updating, condition identification, condition assessment, and decision making of the concerned structure require enough efficient data to reflect structural properties. It is a common way to achieve the related data by sensors. However, restricted by testing technology and cost, a practical problem is how to select a set with a minimum number of sensor locations from all possibilities [1]. Optimal sensor placement (OSP) has been one of the more important tasks in dynamic test, especially for those high-rise buildings and complicated structures.The basic principle of sensor placement optimization is that limited measuring points can get enough efficient data which can reflect the structural properties. Recently, most methods aim to achieve the best identification of structural characteristics, including the frequencies and mode shapes. Among these methods, effective independence (EfI) method is the most widely used method, which was proposed by Kammer [2]. For modal kinetic method (MKE), the target of OSP is to ensure the maximization of the measured kinetic mode energy that improves signal-to-noise ratio [3]. Besides, EfI is an iterated version of the MKE with reorthonormalized mode shapes [4]. The modal assurance criterion (MAC) method is based on QR decomposition with column pivoting, and these measuring points are chosen to minimize off-diagonal element of MAC matrix [5].So far, a multitude of researches have been conducted on sensor placement optimization for large-span bridges and high-rise buildings in civil engineering. Yi et al. conducted sensor placement optimization for Dalian International Trade Mansion [1]. Chow et al. presented optimal sensor configuration of a typical transmission tower for the purpose of structural model updating [6]. Wang et al. established the wind and structural health monitoring system implemented on the Runyang Yangtze River Bridge [7]. However, there are few researches on the sensor placement optimization for those structures which are not very high or large, while the vibration modes are very complicated. The medium-speed mill is one of the most commonly used auxiliary power items of equipment in thermal power plants. It generates significant disturbing force when it is running, so it is very important to monitor its vibration condition. The study of Zhu et al. shows that vibration modes of medium-speed mill are very complicated [8].In order to obtain the reliable dynamic response, reasonable condition assessment, and decision making of medium-speed mill, the measuring points should be placed on the locations that reflect abundant structural information of medium-speed mill. So the OSP for vibration test on medium-speed mill is extremely necessary. In this paper, a reliable finite element model (FEM) of medium-speed mill supported by spring vibration-isolated foundation is built, and two OSP methods are applied to the vibration test of medium-speed mill. Finally, considering the actual situation of vibration test, an OSP scheme is proposed in the vibration test on medium-speed mill.


## 2. The Algorithm and Evaluation Criteria of OSP

### 2.1. Effective Independence Method

The aim of the EfI method is to select measurement positions that make the mode shapes of interest as linearly independent as possible while containing sufficient information about the target modal responses in the measurements.The measured structural response can be expressed as(1)u=Φq+w=∑i=1Nqiφi+w,
where Φ is the matrix of FEM target mode shapes, q is the coefficient response vector and also is modal coordinate, and w is a sensor noise vector, assumed stationary random process with a mean value zero. There must be a deviation for corresponding real generalized coordinate q^. Assuming that this process is an effective unbiased estimation, the covariance matrix J is obtained as
(2)J=Eq-q^q-q^T=σ-2ΦTΦ-1=Q-1,
where Q is the Fisher information matrix (FIM) by assuming the measured noise is independent and has the same statistical properties. The matrix Q can be written as
(3)Q=σ-2ΦTΦ=σ-2Α.
Then, the maximization of Q is equivalent to the maximization of A, and thus A can be used to simplify FIM. Constructing matrix ED,
(4)ED=ΦΦTΦ-1ΦT,
where ED is the effective independence allocation matrix. The diagonal elements of the ED represent the contribution of the candidate points to the modal matrix linearly independent.


### 2.2. The Algorithm for OSP Based on QR Decomposition

QR decomposition is derived from the maximization of the FIM. The linearly independent row vectors are extracted from modal matrix by QR decomposition that ensures a larger norm of matrix. These degrees of freedom (DOFs) corresponding with linearly independent row vectors can be used as an ideal placement for sensors. The decomposition can be expressed as(5)A=ΦTE=QR,
where E is the unit conversion matrix, descending order according to the value of the diagonal elements of R. We choose the modal DOFs in accordance with the sequence of the elements in E.


### 2.3. The Evaluation Criteria of OSP

#### 2.3.1. Fisher Information Matrix Determinant

According to the principle of EfI method, the greater the trace or the value of FIM determinant is, the smaller the covarianceJ is and the more effective the modal coordinate is. Meanwhile, FIM also measures the abundance modal information from the response of structures, and the greater the determinant value is, the more abundant modal information it contains [5].


#### 2.3.2. Modal Assurance Criterion

The quantity of measured DOFs is far less than the whole DOFs of structural model and the responses are affected by measurement accuracy and noise; thus, the modal vector measured cannot retain the original space property. MAC is a good method to evaluate the space angle of modal vectors [9]. The formula is expressed as
(6)MACij=ϕiTϕj2ϕiTϕiϕjTϕj,
where ϕi and ϕj are the ith and jth modal vectors, respectively. The off-diagonal elements of MAC matrix should be minimized for OSP, and thus the modal matrix measured can keep better orthogonality.


## 2.1. Effective Independence Method

The aim of the EfI method is to select measurement positions that make the mode shapes of interest as linearly independent as possible while containing sufficient information about the target modal responses in the measurements.The measured structural response can be expressed as(1)u=Φq+w=∑i=1Nqiφi+w,
where Φ is the matrix of FEM target mode shapes, q is the coefficient response vector and also is modal coordinate, and w is a sensor noise vector, assumed stationary random process with a mean value zero. There must be a deviation for corresponding real generalized coordinate q^. Assuming that this process is an effective unbiased estimation, the covariance matrix J is obtained as
(2)J=Eq-q^q-q^T=σ-2ΦTΦ-1=Q-1,
where Q is the Fisher information matrix (FIM) by assuming the measured noise is independent and has the same statistical properties. The matrix Q can be written as
(3)Q=σ-2ΦTΦ=σ-2Α.
Then, the maximization of Q is equivalent to the maximization of A, and thus A can be used to simplify FIM. Constructing matrix ED,
(4)ED=ΦΦTΦ-1ΦT,
where ED is the effective independence allocation matrix. The diagonal elements of the ED represent the contribution of the candidate points to the modal matrix linearly independent.


## 2.2. The Algorithm for OSP Based on QR Decomposition

QR decomposition is derived from the maximization of the FIM. The linearly independent row vectors are extracted from modal matrix by QR decomposition that ensures a larger norm of matrix. These degrees of freedom (DOFs) corresponding with linearly independent row vectors can be used as an ideal placement for sensors. The decomposition can be expressed as(5)A=ΦTE=QR,
where E is the unit conversion matrix, descending order according to the value of the diagonal elements of R. We choose the modal DOFs in accordance with the sequence of the elements in E.


## 2.3. The Evaluation Criteria of OSP

### 2.3.1. Fisher Information Matrix Determinant

According to the principle of EfI method, the greater the trace or the value of FIM determinant is, the smaller the covarianceJ is and the more effective the modal coordinate is. Meanwhile, FIM also measures the abundance modal information from the response of structures, and the greater the determinant value is, the more abundant modal information it contains [5].


### 2.3.2. Modal Assurance Criterion

The quantity of measured DOFs is far less than the whole DOFs of structural model and the responses are affected by measurement accuracy and noise; thus, the modal vector measured cannot retain the original space property. MAC is a good method to evaluate the space angle of modal vectors [9]. The formula is expressed as
(6)MACij=ϕiTϕj2ϕiTϕiϕjTϕj,
where ϕi and ϕj are the ith and jth modal vectors, respectively. The off-diagonal elements of MAC matrix should be minimized for OSP, and thus the modal matrix measured can keep better orthogonality.


## 2.3.1. Fisher Information Matrix Determinant

According to the principle of EfI method, the greater the trace or the value of FIM determinant is, the smaller the covarianceJ is and the more effective the modal coordinate is. Meanwhile, FIM also measures the abundance modal information from the response of structures, and the greater the determinant value is, the more abundant modal information it contains [5].


## 2.3.2. Modal Assurance Criterion

The quantity of measured DOFs is far less than the whole DOFs of structural model and the responses are affected by measurement accuracy and noise; thus, the modal vector measured cannot retain the original space property. MAC is a good method to evaluate the space angle of modal vectors [9]. The formula is expressed as
(6)MACij=ϕiTϕj2ϕiTϕiϕjTϕj,
where ϕi and ϕj are the ith and jth modal vectors, respectively. The off-diagonal elements of MAC matrix should be minimized for OSP, and thus the modal matrix measured can keep better orthogonality.


## 3. OSP for Vibration Test on Medium-Speed Mill

### 3.1. Description of the Medium-Speed Mill

Medium-speed mill is an important auxiliary dynamic machine in the thermal power plants. To mitigate the vibration, the spring vibration-isolated foundation is normally used to support the mill. A spring vibration-isolated foundation for medium-speed mill is typically comprised of upper and lower bedplates and spring isolators. An elevation view of a vibration-isolated foundation for medium-speed mill is shown in Figure1. As shown, there are two buttresses over the lower bedplate, and the distance between the upper bedplate and buttresses is 330 mm. Spring isolators are placed between the two buttresses and the upper bedplate. Twelve sets of spring isolators are placed on each buttress; thus, a total of 24 sets of isolators are placed.Figure 1
Spring vibration-isolated foundation of medium-speed mill.


### 3.2. Finite Element Analysis

#### 3.2.1. Modeling

A medium-speed mill in the Yuzhou power plant (Henan Province, China) was used as the prototype for the developed finite element model. In order to provide the input data for the optimization analysis, a three-dimensional finite element coupled model of a medium-speed mill and its foundation was built using ANSYS (version 12.0) [10]. SOLID185 elements were used to simulate the foundation, the mill, and the motor; SHELL181 elements were used for the frame and chassis; BEAM188 elements were used for the pressure frame; COMBINE14 elements were used for the isolators; and MASS21 elements were used for the load.The frame, motor, and mill were rigidly connected with the upper bedplate of the foundation. The boundary condition of the lower bedplate of the foundation was simulated using spring elements, which represented the constraint effect of the soil on the foundation. For the coordinate system of the model,X represented the transverse direction of the foundation of medium-speed mill and Y represented the longitudinal direction. The mesh size of the finite element model was between 100 and 200 mm. The finite element model was shown in Figure 2, and the primary parameters of the mill and vibration isolator were listed in Table 1.Table 1
Parameters of the mill and vibration isolator.
Medium-speed millZGM113G
Total weight of the mill/t
157.40
Total weight of the foundation/t
398.74
Total weight of the mill and foundation/t
556.14
Mass ratio of the foundation to the mill
2.53
Vibration isolatorSPV2076.12.12
Vertical stiffness of the vibration isolator/kN/mm
12.10
Horizontal stiffness of the vibration isolator/kN/mm
9.50
Vertical damping factor/kNs/m
240
Horizontal damping factor/kNs/m
160
Damping ratio of the structure
0.02
Disturbing force/kN
X direction = 10
Y direction = 10
Z direction = 40Figure 2
Finite element model.


#### 3.2.2. Modal Analysis

After building a FE model of the medium-speed mill with spring vibration-isolated foundation, a modal analysis was done using the block Lanczos method. The vibration modal shapes and corresponding frequencies were shown in Figure3 and Table 2.Table 2
Vibration frequencies.
Order
Frequencyf/Hz
Vibration modal shape
1
1.828
X direction horizontal vibration
2
2.003
Y direction horizontal vibration
3
3.578
Vertical vibration
4
4.078
Rotation around horizontal axis of the mill
5
4.235
Torsion around the vertical axis
6
4.422
Horizontal oblique translational vibrationFigure 3
First three modes.
(a)
First modal shape
(b)
Second modal shape
(c)
Third modal shape


### 3.3. Selection of Candidate Measuring Points

The medium-speed mill is an axial symmetric structure. According to ANSYS analysis, the primary vibration mode shapes of medium-speed mill include horizontal longitudinal vibration, horizontal transverse vibration, and vertical vibration. And also considering the cost and technical factors, measuring points were placed along one side of the medium-speed mill. 12 measuring points were set as candidate points. Each point had DOFs ofX, Y, and Z direction; thus, there were a total of 36 DOFs. EfI and QR methods were employed to optimize the sensors placement of medium-speed mill and the first 3, 4, 5, and 6 were taken as the monitoring modal shapes, respectively. The placement of candidate measuring points was shown in Figure 4.Figure 4
Placement of candidate measuring points.36 DOFs were analyzed preliminarily using EfI method with MATLAB software. According toED, the contribution of 12 DOFs in Z direction to modal matrix linear independence was far less. There is no sufficient reason to believe that the FE model can simulate the actual structure accurately, so removing the DOFs in Z direction just based on the optimization results of 36 DOFs, the vibration characteristic of medium-speed mill may not be understood comprehensively. Motivated by the above consideration, it is determined to optimize the placement of 12 DOFs in X, Y, and Z direction, respectively.


### 3.4. Results of OSP

The OSP of 12 candidate measuring points was conducted using EfI method and QR decomposition algorithm for different target mode shapes inX, Y, and Z direction, respectively, and the results of OSP were presented in Table 3. It was noted that those measuring points in brackets were suboptimal ones. The suboptimal measuring points were not included in the optimization results in this paper, while the referred optimization results meant that the number of measuring points was consistent with the number of target mode shapes.Table 3
Optimization results of measuring points.
Order
Direction
X (EfI)
X (QR)
Y (EfI)
Y (QR)
Z (EfI)
Z (QR)
6
1
1
3
1
1
1
4
4
4
4
2
2
7
8
5
5
4
6
8
9
9
9
5
7
9
10
10
10
10
10
11
11
11
11
11
11
5
1
1
3
1
1
1
4
4
4
4
2
2
8
9
5
5
6
6
9
10
10
10
10
10
11
11
11
11
11
11
(7)
(2)
(7)
(2)
(4)
(8)
4
1
1
3
1
1
1
4
4
4
4
4
6
10
10
7
7
10
10
11
11
10
10
11
11
(9)
(8)
(1)
(8)
(6)
(8)
(6)
(2)
(5)
(2)
(9)
(2)
3
1
1
3
1
1
1
6
6
4
4
10
6
11
11
10
10
11
10
(5)
(7)
(1)
(7)
(6)
(7)
(3)
(8)
(2)
(8)
(9)
(8)
(4)
(2)
(5)
(2)
(4)
(2)In Table3, it was found that two optimization results did not agree with each other exactly. The optimization results of the first 3 and 4 mode shapes in X direction were consistent exactly, but they were different for one and two measuring points in the first 5 and 6 mode shapes, respectively. The optimization results in X direction had the trend that the fewer the number of target mode shapes was, the smaller the difference in two results was. Although the same trend did not appear in Y and Z directions, the optimization results were still consistent basically. Therefore, EfI and QR decomposition algorithms in OSP for medium-speed mill were feasible.FIM determinant behavior was recorded in terms of percentage of its initial value against the decrease in candidate measuring points [11], and the results were shown in Figure 5. It was noted that the FIM determinant in lower orders of target mode shapes had larger value in case of the same number of measuring points, which meant that these measuring points were more efficient in getting concerned information of the medium-speed mill. On the other hand, the value of FIM determinant increased more significantly when adding a measuring point in case of the same number of measuring points, which meant that the optimization results in larger orders were more sensitive to the number of measuring points. So the target modal played a significant role in the final scheme of the OSP.Figure 5
Fisher information matrix determinant.
(a)
X direction
(b)
Y direction
(c)
Z directionThe average values of off-diagonal elements of MAC matrix were shown in Figure6. The average values of off-diagonal elements of MAC matrix were greater even if when 12 candidate measuring points were all placed. It was explained that the vibration modes of medium-speed mill were very complicated and even 12 measuring points cannot capture its dynamic characteristics expectedly. But the average value of off-diagonal elements of MAC increased a little in Y direction using these two optimal methods to decrease the number of measuring points, while the average values in X and Z directions decreased, in particular in Z direction. In other words, the expected optimal results could be achieved. Meanwhile, it could be confirmed that EfI method was much better due to smaller average values of MAC.Figure 6
Average values of MAC matrix off-diagonal elements.
(a)
X direction
(b)
Y direction
(c)
Z directionIt was easy to find that the selection of suboptimal measuring points using two methods was significantly different. The values of FIM determinant after adding the suboptimal points to the optimization results and the final results were shown in Table4. The value of FIM using EfI method was greater than that of QR decomposition method. The effect of optimization using EfI method was better, and the corresponding measuring points had a stronger ability to get the information of the medium-speed mill.Table 4
Results of Fisher information matrix determinant.
Direction
Order
X (EfI)
X (QR)
Y (EfI)
Y (QR)
Z (EfI)
Z (QR)
5
3.53
×
10
-
17
1.03
×
10
-
1
7
7.17
×
10
-
23
5.71
×
10
-
2
3
2.61
×
10
-
20
1.24
×
10
-
19
4
2.49
×
10
-
11
1.19
×
10
-
11
4.28
×
10
-
20
3.14
×
10
-
20
8.60
×
10
-
15
5.50
×
10
-
15
3
1.27
×
10
-
8
1.17
×
10
-
8
1.51
×
10
-
13
1.09
×
10
-
13
3.41
×
10
-
9
1.35
×
10
-
9


## 3.1. Description of the Medium-Speed Mill

Medium-speed mill is an important auxiliary dynamic machine in the thermal power plants. To mitigate the vibration, the spring vibration-isolated foundation is normally used to support the mill. A spring vibration-isolated foundation for medium-speed mill is typically comprised of upper and lower bedplates and spring isolators. An elevation view of a vibration-isolated foundation for medium-speed mill is shown in Figure1. As shown, there are two buttresses over the lower bedplate, and the distance between the upper bedplate and buttresses is 330 mm. Spring isolators are placed between the two buttresses and the upper bedplate. Twelve sets of spring isolators are placed on each buttress; thus, a total of 24 sets of isolators are placed.Figure 1
Spring vibration-isolated foundation of medium-speed mill.


## 3.2. Finite Element Analysis

### 3.2.1. Modeling

A medium-speed mill in the Yuzhou power plant (Henan Province, China) was used as the prototype for the developed finite element model. In order to provide the input data for the optimization analysis, a three-dimensional finite element coupled model of a medium-speed mill and its foundation was built using ANSYS (version 12.0) [10]. SOLID185 elements were used to simulate the foundation, the mill, and the motor; SHELL181 elements were used for the frame and chassis; BEAM188 elements were used for the pressure frame; COMBINE14 elements were used for the isolators; and MASS21 elements were used for the load.The frame, motor, and mill were rigidly connected with the upper bedplate of the foundation. The boundary condition of the lower bedplate of the foundation was simulated using spring elements, which represented the constraint effect of the soil on the foundation. For the coordinate system of the model,X represented the transverse direction of the foundation of medium-speed mill and Y represented the longitudinal direction. The mesh size of the finite element model was between 100 and 200 mm. The finite element model was shown in Figure 2, and the primary parameters of the mill and vibration isolator were listed in Table 1.Table 1
Parameters of the mill and vibration isolator.
Medium-speed millZGM113G
Total weight of the mill/t
157.40
Total weight of the foundation/t
398.74
Total weight of the mill and foundation/t
556.14
Mass ratio of the foundation to the mill
2.53
Vibration isolatorSPV2076.12.12
Vertical stiffness of the vibration isolator/kN/mm
12.10
Horizontal stiffness of the vibration isolator/kN/mm
9.50
Vertical damping factor/kNs/m
240
Horizontal damping factor/kNs/m
160
Damping ratio of the structure
0.02
Disturbing force/kN
X direction = 10
Y direction = 10
Z direction = 40Figure 2
Finite element model.


### 3.2.2. Modal Analysis

After building a FE model of the medium-speed mill with spring vibration-isolated foundation, a modal analysis was done using the block Lanczos method. The vibration modal shapes and corresponding frequencies were shown in Figure3 and Table 2.Table 2
Vibration frequencies.
Order
Frequencyf/Hz
Vibration modal shape
1
1.828
X direction horizontal vibration
2
2.003
Y direction horizontal vibration
3
3.578
Vertical vibration
4
4.078
Rotation around horizontal axis of the mill
5
4.235
Torsion around the vertical axis
6
4.422
Horizontal oblique translational vibrationFigure 3
First three modes.
(a)
First modal shape
(b)
Second modal shape
(c)
Third modal shape


## 3.2.1. Modeling

A medium-speed mill in the Yuzhou power plant (Henan Province, China) was used as the prototype for the developed finite element model. In order to provide the input data for the optimization analysis, a three-dimensional finite element coupled model of a medium-speed mill and its foundation was built using ANSYS (version 12.0) [10]. SOLID185 elements were used to simulate the foundation, the mill, and the motor; SHELL181 elements were used for the frame and chassis; BEAM188 elements were used for the pressure frame; COMBINE14 elements were used for the isolators; and MASS21 elements were used for the load.The frame, motor, and mill were rigidly connected with the upper bedplate of the foundation. The boundary condition of the lower bedplate of the foundation was simulated using spring elements, which represented the constraint effect of the soil on the foundation. For the coordinate system of the model,X represented the transverse direction of the foundation of medium-speed mill and Y represented the longitudinal direction. The mesh size of the finite element model was between 100 and 200 mm. The finite element model was shown in Figure 2, and the primary parameters of the mill and vibration isolator were listed in Table 1.Table 1
Parameters of the mill and vibration isolator.
Medium-speed millZGM113G
Total weight of the mill/t
157.40
Total weight of the foundation/t
398.74
Total weight of the mill and foundation/t
556.14
Mass ratio of the foundation to the mill
2.53
Vibration isolatorSPV2076.12.12
Vertical stiffness of the vibration isolator/kN/mm
12.10
Horizontal stiffness of the vibration isolator/kN/mm
9.50
Vertical damping factor/kNs/m
240
Horizontal damping factor/kNs/m
160
Damping ratio of the structure
0.02
Disturbing force/kN
X direction = 10
Y direction = 10
Z direction = 40Figure 2
Finite element model.


## 3.2.2. Modal Analysis

After building a FE model of the medium-speed mill with spring vibration-isolated foundation, a modal analysis was done using the block Lanczos method. The vibration modal shapes and corresponding frequencies were shown in Figure3 and Table 2.Table 2
Vibration frequencies.
Order
Frequencyf/Hz
Vibration modal shape
1
1.828
X direction horizontal vibration
2
2.003
Y direction horizontal vibration
3
3.578
Vertical vibration
4
4.078
Rotation around horizontal axis of the mill
5
4.235
Torsion around the vertical axis
6
4.422
Horizontal oblique translational vibrationFigure 3
First three modes.
(a)
First modal shape
(b)
Second modal shape
(c)
Third modal shape


## 3.3. Selection of Candidate Measuring Points

The medium-speed mill is an axial symmetric structure. According to ANSYS analysis, the primary vibration mode shapes of medium-speed mill include horizontal longitudinal vibration, horizontal transverse vibration, and vertical vibration. And also considering the cost and technical factors, measuring points were placed along one side of the medium-speed mill. 12 measuring points were set as candidate points. Each point had DOFs ofX, Y, and Z direction; thus, there were a total of 36 DOFs. EfI and QR methods were employed to optimize the sensors placement of medium-speed mill and the first 3, 4, 5, and 6 were taken as the monitoring modal shapes, respectively. The placement of candidate measuring points was shown in Figure 4.Figure 4
Placement of candidate measuring points.36 DOFs were analyzed preliminarily using EfI method with MATLAB software. According toED, the contribution of 12 DOFs in Z direction to modal matrix linear independence was far less. There is no sufficient reason to believe that the FE model can simulate the actual structure accurately, so removing the DOFs in Z direction just based on the optimization results of 36 DOFs, the vibration characteristic of medium-speed mill may not be understood comprehensively. Motivated by the above consideration, it is determined to optimize the placement of 12 DOFs in X, Y, and Z direction, respectively.


## 3.4. Results of OSP

The OSP of 12 candidate measuring points was conducted using EfI method and QR decomposition algorithm for different target mode shapes inX, Y, and Z direction, respectively, and the results of OSP were presented in Table 3. It was noted that those measuring points in brackets were suboptimal ones. The suboptimal measuring points were not included in the optimization results in this paper, while the referred optimization results meant that the number of measuring points was consistent with the number of target mode shapes.Table 3
Optimization results of measuring points.
Order
Direction
X (EfI)
X (QR)
Y (EfI)
Y (QR)
Z (EfI)
Z (QR)
6
1
1
3
1
1
1
4
4
4
4
2
2
7
8
5
5
4
6
8
9
9
9
5
7
9
10
10
10
10
10
11
11
11
11
11
11
5
1
1
3
1
1
1
4
4
4
4
2
2
8
9
5
5
6
6
9
10
10
10
10
10
11
11
11
11
11
11
(7)
(2)
(7)
(2)
(4)
(8)
4
1
1
3
1
1
1
4
4
4
4
4
6
10
10
7
7
10
10
11
11
10
10
11
11
(9)
(8)
(1)
(8)
(6)
(8)
(6)
(2)
(5)
(2)
(9)
(2)
3
1
1
3
1
1
1
6
6
4
4
10
6
11
11
10
10
11
10
(5)
(7)
(1)
(7)
(6)
(7)
(3)
(8)
(2)
(8)
(9)
(8)
(4)
(2)
(5)
(2)
(4)
(2)In Table3, it was found that two optimization results did not agree with each other exactly. The optimization results of the first 3 and 4 mode shapes in X direction were consistent exactly, but they were different for one and two measuring points in the first 5 and 6 mode shapes, respectively. The optimization results in X direction had the trend that the fewer the number of target mode shapes was, the smaller the difference in two results was. Although the same trend did not appear in Y and Z directions, the optimization results were still consistent basically. Therefore, EfI and QR decomposition algorithms in OSP for medium-speed mill were feasible.FIM determinant behavior was recorded in terms of percentage of its initial value against the decrease in candidate measuring points [11], and the results were shown in Figure 5. It was noted that the FIM determinant in lower orders of target mode shapes had larger value in case of the same number of measuring points, which meant that these measuring points were more efficient in getting concerned information of the medium-speed mill. On the other hand, the value of FIM determinant increased more significantly when adding a measuring point in case of the same number of measuring points, which meant that the optimization results in larger orders were more sensitive to the number of measuring points. So the target modal played a significant role in the final scheme of the OSP.Figure 5
Fisher information matrix determinant.
(a)
X direction
(b)
Y direction
(c)
Z directionThe average values of off-diagonal elements of MAC matrix were shown in Figure6. The average values of off-diagonal elements of MAC matrix were greater even if when 12 candidate measuring points were all placed. It was explained that the vibration modes of medium-speed mill were very complicated and even 12 measuring points cannot capture its dynamic characteristics expectedly. But the average value of off-diagonal elements of MAC increased a little in Y direction using these two optimal methods to decrease the number of measuring points, while the average values in X and Z directions decreased, in particular in Z direction. In other words, the expected optimal results could be achieved. Meanwhile, it could be confirmed that EfI method was much better due to smaller average values of MAC.Figure 6
Average values of MAC matrix off-diagonal elements.
(a)
X direction
(b)
Y direction
(c)
Z directionIt was easy to find that the selection of suboptimal measuring points using two methods was significantly different. The values of FIM determinant after adding the suboptimal points to the optimization results and the final results were shown in Table4. The value of FIM using EfI method was greater than that of QR decomposition method. The effect of optimization using EfI method was better, and the corresponding measuring points had a stronger ability to get the information of the medium-speed mill.Table 4
Results of Fisher information matrix determinant.
Direction
Order
X (EfI)
X (QR)
Y (EfI)
Y (QR)
Z (EfI)
Z (QR)
5
3.53
×
10
-
17
1.03
×
10
-
1
7
7.17
×
10
-
23
5.71
×
10
-
2
3
2.61
×
10
-
20
1.24
×
10
-
19
4
2.49
×
10
-
11
1.19
×
10
-
11
4.28
×
10
-
20
3.14
×
10
-
20
8.60
×
10
-
15
5.50
×
10
-
15
3
1.27
×
10
-
8
1.17
×
10
-
8
1.51
×
10
-
13
1.09
×
10
-
13
3.41
×
10
-
9
1.35
×
10
-
9


## 4. The Experimental Results of OSP

### 4.1. Final Scheme of OSP

Based on the analysis of Section3.2.2, the X-translational vibration was the primary vibration mode of medium-speed mill, so the optimal results of X direction in the first 3 mode shapes using EfI method were employed. Considering the stability of data and security of sensors, the measuring point 1 was changed to point 2, and they had the same testing condition. In addition, the data collected by the measuring point 2 was an important indicator to decide the working state of medium-speed mill. In order to reduce the effects of environmental and operational variability and noise contamination, the measuring points 2 and 6 with more modal kinetic energy were selected to improve the signal-to-noise ratio [3, 12–14]. The analysis of OSP showed that the contribution of the actual measuring point 12 in X, Y, and Z directions is minimal, and it should be removed in theory. Considering that one of the aims of vibration test is to understand the isolation effect of the spring vibration-isolated foundation, the measuring point 12 was retained. From Table 3, the contribution of the measuring point 11 on the vibration of medium-speed mill in three directions was great, and the vertical vibration of medium-speed mill also should be tested. Therefore, the measuring points 2, 6, 11, and 12 were chosen as the final scheme of OSP.In order to place the sensors in site easily, the measuring points were renumbered, which were shown in Figure7. Among these test points, the HTP1 sensor was used to measure the horizontal vibration of the lower bedplate of the foundation or ground; the HTP2 sensor was attached to the upper bedplate; the HTP3 sensor was attached to the middle of the mill; and the HTP4 sensor was attached to the top of the mill. The VTP1 sensor was used to measure the vertical vibration of the lower bedplate of the foundation or ground; and the VTP2 sensor was attached to the upper bedplate of the foundation.Figure 7
Placement of sensors (HTTP: horizontal test point, VTP: vertical test point).


### 4.2. Field Tests

In order to verify the feasibility of OSP results, the medium-speed mill in the Yuzhou thermal power plant was tested in site. In the field test, 891-II type vibration measurement instruments including six sensors, amplifiers with a frequency response ranging from 0.5 Hz to 100 Hz, an INV306 intelligent signal acquisition and processing analysis system, and a DASP data processing program were used. The sensors were placed according to the final scheme of OSP. The sensors used in the experiment were manufactured by the Institute of Engineering Mechanics, China Earthquake Administration. The sensors have four gears, which are used to test the large displacement, small displacement, velocity, and acceleration. Actually, they are displacement sensors. In this test, the gear of small displacement was used. The elevation view of the tested medium-speed mill was shown in Figure8, and the sensor in close-up was shown in Figure 9.Figure 8
The elevation view of the tested medium-speed mill.Figure 9
The sensor shown in close-up.The vibration test on the medium-speed mill was conducted in the plant, so the variability of physical property with thermal effect can be ignored [13]. The primary noise contamination was sourced from other nearby dynamic machines. The vibration data was recorded when the medium-speed mill did not run, which was used in the next filtering [12].Field tests were conducted during the working steady stage of the medium-speed mill, and the vibration displacements were tested by the described sensors. The collection time was 60 seconds per sample. In this condition, the vibrations of the foundation and mill body were assumed to be the ergodic and stationary random processes; the collected samples can thus describe the vibration characteristics of the mill. The noise was eliminated using filter techniques with DASP software and the typical vibration curves for the mill obtained from the test were shown in Figure10. It should be noted that only vibration displacement ranging from 0 to 10 second were shown in Figure 10.Figure 10
Horizontal vibration displacement inX direction.


### 4.3. Comparison of Vibration Displacement

Harmonic response analysis at the service load was performed by ANSYS. The analytical results of the vibration amplitude of each measuring point were compared with the test results, which was presented in Table5. The test results in Table 5 were the maximum vibration displacements of these measuring points. The calculation results generally agreed with the field test results in X direction, which indicated that the FE model could simulate the dynamic characteristics of the medium-speed mill. It was also noted that the simulation in Z and Y directions was less effective than that in X direction. From the location of measuring points, the test point 4 was more satisfactory, and the test points 2 and 3 were followed. Abundant dynamic characteristics of medium-speed mill were achieved based on the final sensors placement scheme, and the finite element analysis results were basically consistent with the field test results, which indicated the sensor placement optimization for vibration test on the medium-speed mill was effective and economical.Table 5
Comparison of vibration displacement between the finite element analysis results and the field test results.
Measuring points
2(X)/mm
3(X)/mm
4(X)/mm
2(Y)/mm
3(Y)/mm
4(Y)/mm
2(Z)/mm
Test results
0.062
0.132
0.421
0.103
0.226
0.246
0.106
Calculation results
0.051
0.101
0.385
0.080
0.169
0.202
0.082
Error (%)
17.5
23.5
8.6
22.4
25.1
18.3
22.6


## 4.1. Final Scheme of OSP

Based on the analysis of Section3.2.2, the X-translational vibration was the primary vibration mode of medium-speed mill, so the optimal results of X direction in the first 3 mode shapes using EfI method were employed. Considering the stability of data and security of sensors, the measuring point 1 was changed to point 2, and they had the same testing condition. In addition, the data collected by the measuring point 2 was an important indicator to decide the working state of medium-speed mill. In order to reduce the effects of environmental and operational variability and noise contamination, the measuring points 2 and 6 with more modal kinetic energy were selected to improve the signal-to-noise ratio [3, 12–14]. The analysis of OSP showed that the contribution of the actual measuring point 12 in X, Y, and Z directions is minimal, and it should be removed in theory. Considering that one of the aims of vibration test is to understand the isolation effect of the spring vibration-isolated foundation, the measuring point 12 was retained. From Table 3, the contribution of the measuring point 11 on the vibration of medium-speed mill in three directions was great, and the vertical vibration of medium-speed mill also should be tested. Therefore, the measuring points 2, 6, 11, and 12 were chosen as the final scheme of OSP.In order to place the sensors in site easily, the measuring points were renumbered, which were shown in Figure7. Among these test points, the HTP1 sensor was used to measure the horizontal vibration of the lower bedplate of the foundation or ground; the HTP2 sensor was attached to the upper bedplate; the HTP3 sensor was attached to the middle of the mill; and the HTP4 sensor was attached to the top of the mill. The VTP1 sensor was used to measure the vertical vibration of the lower bedplate of the foundation or ground; and the VTP2 sensor was attached to the upper bedplate of the foundation.Figure 7
Placement of sensors (HTTP: horizontal test point, VTP: vertical test point).


## 4.2. Field Tests

In order to verify the feasibility of OSP results, the medium-speed mill in the Yuzhou thermal power plant was tested in site. In the field test, 891-II type vibration measurement instruments including six sensors, amplifiers with a frequency response ranging from 0.5 Hz to 100 Hz, an INV306 intelligent signal acquisition and processing analysis system, and a DASP data processing program were used. The sensors were placed according to the final scheme of OSP. The sensors used in the experiment were manufactured by the Institute of Engineering Mechanics, China Earthquake Administration. The sensors have four gears, which are used to test the large displacement, small displacement, velocity, and acceleration. Actually, they are displacement sensors. In this test, the gear of small displacement was used. The elevation view of the tested medium-speed mill was shown in Figure8, and the sensor in close-up was shown in Figure 9.Figure 8
The elevation view of the tested medium-speed mill.Figure 9
The sensor shown in close-up.The vibration test on the medium-speed mill was conducted in the plant, so the variability of physical property with thermal effect can be ignored [13]. The primary noise contamination was sourced from other nearby dynamic machines. The vibration data was recorded when the medium-speed mill did not run, which was used in the next filtering [12].Field tests were conducted during the working steady stage of the medium-speed mill, and the vibration displacements were tested by the described sensors. The collection time was 60 seconds per sample. In this condition, the vibrations of the foundation and mill body were assumed to be the ergodic and stationary random processes; the collected samples can thus describe the vibration characteristics of the mill. The noise was eliminated using filter techniques with DASP software and the typical vibration curves for the mill obtained from the test were shown in Figure10. It should be noted that only vibration displacement ranging from 0 to 10 second were shown in Figure 10.Figure 10
Horizontal vibration displacement inX direction.


## 4.3. Comparison of Vibration Displacement

Harmonic response analysis at the service load was performed by ANSYS. The analytical results of the vibration amplitude of each measuring point were compared with the test results, which was presented in Table5. The test results in Table 5 were the maximum vibration displacements of these measuring points. The calculation results generally agreed with the field test results in X direction, which indicated that the FE model could simulate the dynamic characteristics of the medium-speed mill. It was also noted that the simulation in Z and Y directions was less effective than that in X direction. From the location of measuring points, the test point 4 was more satisfactory, and the test points 2 and 3 were followed. Abundant dynamic characteristics of medium-speed mill were achieved based on the final sensors placement scheme, and the finite element analysis results were basically consistent with the field test results, which indicated the sensor placement optimization for vibration test on the medium-speed mill was effective and economical.Table 5
Comparison of vibration displacement between the finite element analysis results and the field test results.
Measuring points
2(X)/mm
3(X)/mm
4(X)/mm
2(Y)/mm
3(Y)/mm
4(Y)/mm
2(Z)/mm
Test results
0.062
0.132
0.421
0.103
0.226
0.246
0.106
Calculation results
0.051
0.101
0.385
0.080
0.169
0.202
0.082
Error (%)
17.5
23.5
8.6
22.4
25.1
18.3
22.6


## 5. Conclusion

In this paper, the common methods and evaluation criteria for OSP had been summarized. The OSP of vibration test on medium-speed mill supported by spring vibration-isolated foundation was conducted using EfI and QR decomposition methods. It was found that these two methods for vibration test on medium-speed mill were all feasible, although this was different for some measuring points in different directions, and the difference became smaller with the decrease of the number of target mode shapes. For the evaluation criteria of FIM determinant, the effect was better in the case of small number of target mode shapes when the measuring points were certain; on the other hand, the effects in more target mode shapes were much more sensitive to the number of measuring points. For MAC, EfI method was more efficient than QR decomposition method. In addition, the OSP based on EfI method can achieve more efficient suboptimal measuring points.The final scheme of OSP was based on the result that the first three mode shapes inX direction were set as the targeted modal shapes. Six sensors were placed to achieve the dynamic response of the medium-speed mill, and the effect of the environmental and operational variability and noise contamination was reduced by filtering technique. The field test results were basically consistent with the calculation results by ANSYS, and the FEM of the medium-speed mill with spring vibration-isolated foundation could simulate the vibration characteristic well. The sensor placement optimization for vibration test on the medium-speed mill was effective and economical.



---

*Source: 690196-2015-05-05.xml*