#### Abstract

The measurement of hydraulic cylinder displacement has been addressed from different fields. The detection principle of magnetic grating is able to realize the high integration and accuracy. In this paper, a signal response quality evaluation algorithm for devising and optimizing a high-accuracy displacement measuring system is proposed. On the basic of signal response quality evaluation method, structure variables are optimized to enhance the working performance. By defining the parameters, an optimum structure cylinder prototype is made and tested to provide better estimates. Experimental results on working characteristic are presented to verify the effectiveness of the optimized structure. The efficiency of the proposed signal response quality evaluation function is therefore demonstrated through the working performance.

#### 1. Introduction

In the past decades, the development of electro-hydraulic actuators (EHA) allows flight-control systems to balance requirements of high power, lightweight, safety, fast response, and continuity of service [1]. Actuation system has introduced to meet the control surface loading and response demands. The Airbus A380 employs EHA in its flight surface control system as a state-of-the-art example for exclusion of complex hydraulic lines from controlling systems [2]. In other words, the architecture of EHA is sufficiently robust to failures so as to ensure flight safety and minimize the size of the control system [3]. Within the design procedure, considering system integration, miniaturization, and measuring precision, the displacement measuring methodology of hydraulic cylinders is one of the most fundamental and promising research settings in the field of EHA.

Studies emerged have paved a way for research on cylinder displacement sensing methods of EHA. In the industrial practice, two kinds of position sensing devices are used: one is linear variable inductor displacement transducer (LVDT), and the other is magnetostrictive displacement sensor [4, 5]. The resolution and accuracy of such a sensor are mainly determined by measuring principle and working conditions. Since the distribution of actuation system requires high integration under a severe working environment, both measurements are difficult for installation and maintenance due to their large volumes and complex structures. In such cases, it is advisable to introduce a more appropriate method for cylinder stroke measurement.

Previously, Murakami and Kitsunai applied the idea of equipping the hydraulic cylinder with a stroke sensor based on a magnetoresistor sensor to detect positions [6]. Their results indicated that a high magnetic field and a large field length are required for sensing devices working. Compared with the approaches mentioned above, magnetic-grating-like stroke-sensing cylinder, which is also known as magnetic-grating-like hydraulic cylinder integrated displacement sensor, provides a more effective alternative. This sensing device is to use simple working elements that will minimize the size of EHA by minimizing the dimension of the cylinder [7]. With magnetic-grating-like stroke-sensing cylinder, a combination of cylinder and sensor can be designed to handle long working distance measurement as well as matching with the structure of EHA. On the other hand, Yang et al. proposed a stroke-sensing cylinder based on Hall sensors associated with the automation of heavy construction machinery, which reduced the effect of noise with specific electric circuits [8]. The magnetic-grating-like stroke-sensing cylinder investigated in this research is composed of a piston rod, permanent magnet, and a Hall sensor, as shown in Figure 1.

The permanent magnet generates magnetic fields while the sensor detects the magnetic field intensity and thereby measures the displacement. The piston rod is made of ferromagnetic materials with the repeating grooves on the body. The intensity variation is sampled as a function of the relative position between the permanent magnet and the piston rod by positioning the receiving magnetic force across a range of lateral offsets. When the piston rod is moving, the permanent magnet forms constant magnetic scales on the piston rod in combination with the repeating grooves. The magnetic induction line is thus periodically modulated. Supposing the piston rod is moving with a fixed speed, the induction line along with the sensor is simulated, as shown in Figure 2.

**(a) At the position of 0 mm**

**(b) At the position of 0.4 mm**

**(c) At the position of 0.8 mm**

The sensor can effectively detect the magnetic field depending on the magnetic scale, in which the displacement is measured in a magnetic-grating sensing way. The Hall sensor is employed to sense signals on the cylinder’s surface with magnetic scale [9]. The Hall elements measure the variation of the magnetic flux density which is produced by the slots grooved on the surface of the cylinder. The sensing outcome of Hall sensor is likely to be harmonic signal within the working period (Figure 3). The signal glitch is because of the discrete block given through modeling and simulation, which stands for the dimensional precision of working elements.

The problem of magnetic-grating-like stroke-sensing cylinder employment arises in precision measurement, since current laboratory research is still not able to address this issue. Specifically, with the added sliding mode control, Yang et al. did not give the emphasis on accuracy improvement. Most of the previous research aimed at setting up a preliminary system for sensing the stroke of the cylinder. To the best of our knowledge, specific structure design and optimization have not yet been established and no quantitative information is given.

In this paper, we will consider the reason for the low accuracy of magnetic-grating-like hydraulic cylinder integrated displacement sensor. A set of structural optimization method is proposed for addressing the design of the sensor. The research contributions of this paper can be summarized as follows: (a)Response quantification criteria setup: the study develops a function which combined the signal subdivision principle and sensor output used for evaluating the measurement accuracy.(b)Structural parameter optimization: define the structural variables by referring to the response quantification function and compare the influence of each component on optimization results.(c)Performance testing design: to prove the technical efficacy of proposed optimization framework, the prototype is produced and experimental platform is designed. Through the characteristic testing, the measurement accuracy of magnetic-grating-like stroke-sensing cylinder can be obtained.

This paper will introduce signal subdivision-related knowledge and present specific response quality assessment method in Section 2, illustrate parameter optimization process and the optimal structure in Section 3, show the performance testing results achieved of experimental prototype and analyze the measurement accuracy in Section 4, and present future research expectation in Section 5.

#### 2. Response Quality Evaluation Algorithm

##### 2.1. Signal Subdivision Principle

The fundamental principle of stroke measurement is the evaluation of the relative displacement. When the piston rod moves with the distance of one slot, one complete period of harmonic signal will therefore be generated. The signal, which is sensed by the Hall sensor, reflects the variation of magnet field intensity. In this way, the measurement accuracy means the length of one-slot range. According to previous study, the responses of the grating sensor are sine and cosine signals [10]. Hence, the outcome can be normalized by the transformed circumferential angle from 0° to 360° in all cases. Suppose that the minimal length value to be measured is defined by unorganized angles. To start the accurate measurement process, a subdivision 0° to 360° should be generated from every single angle to attain a higher accuracy. Although the magnetic-grating-like stroke-sensing cylinder is inherently less sensitive to changes in displacement, the angle-segmented approach of the harmonic signal offers the advantage that only the response of the sensor is required to demodulate displacement from measured piston rod. In precision measurement, the angle of harmonic signal can be segmented into a far smaller angle change, based on the basic theory of Moiré Fringe signal subdivision technology [11, 12]. At every given angle, the increment to the next is a controllable fixed amount, which can be regulated by the subdivision method. An analysis of the working performance evaluation criteria is presented in this section. If the position of cylinder within one slot is a standard harmonic signal, the testing outcome can be converted to the corresponding segmentation angle. The smaller the angle is divided, the higher precision can be acquired.

##### 2.2. Tangent Method for Signal Subdivision

Lot of research has been conducted to address signal subdivision for performance criteria. Tangent-cotangent subdivision technique is a state-of-the-art method for grating signal segmentation according to Moiré Fringe electric subdivision technique [13, 14]. The fundamental principle of tangent-cotangent subdivision is to investigate the tangent and cotangent value of the two orthogonal signals of different angles, when the output is reflected from a harmonic signal, that is,

For simplicity, the expression is written in terms of angle, rather than radian.

In this way, the displacement can be calculated with the angle got by arctangent algorithm, which eliminates the low precision due to amplitude measurement. Considering the demand of high resolution, the angles will be segmented into even smaller division values. Aiming at establishing the response quantification criteria, the theoretical output of sensor can be calculated in advance. Based on the variation of structure, the response to system input varies. Thus, a function for evaluating the signal response should be set up for system optimization. The response is where is the angle of optimized modeling structure and is the theoretical output angle corresponding to the displacement of piston rod. Accordingly, the maximum error between and stands for the quality of response. Due to the nonlinearity of the response the resolution of (2) will inevitably bring large error to the result. Distinct unit increments of tangent and cotangent value range within the whole cycle (Figure 4).

In line with Figure 2, the change of tangent/cotangent is defined as a nonlinear function according to the angle increments from 0° to 360°, which cannot be resolved by (2). However, the angle transformation period of objective function is from 0° to 90°, where the output can be mathematically revised for calculation. Further, two parts of the curve, 0–45° and 45°–90°, are symmetric to each other along the *y*-axis. Thus, the tangent/cotangent increment of 45°–90° is defined as
where ranges from 45° to 90°.

During the parameter modification process, the variable is revised according to parameter distribution and it determines the quality of response within its angle range domain. of the whole range is taking the cyclic variation of angle into account. where mod represents taking the remainder of target values. The combined objective is defined by

We refer readers to [13–16] for more detailed information on Moiré Fringe electric subdivision and its application.

In order to verify the correctness of the adopted response quality evaluation method, we compute the error of predefined samples with its numerical change in combination with sin/cosine signal. In this test, white noise is added into the sinusoidal signal and cosine signal (Figure 5(a)). The shape of outcome angle error is defined by response evaluation function . Thus, the white noise error can be extracted with the evaluation function. The error response of different signals is given in Figures 5(b) and 5(c). As shown in Figure 4, the real deviation and calculating error can be distinguished only at the points of 45°, 90°, 135°, 225°, and 315°. After the modification, the identification of white noise error is substantially improved in Figure 5(c).

**(a)**

**(b)**

**(c)**

For the purpose of structure optimization, the algorithm is implemented as the following steps:

*Step 1. *Take one working period of a specific sensor as the research sample.

*Step 2. *On the basic of signal subdivision principle, set the sensing response as harmonic signal.

*Step 3. *Divide the period into even parts which correspond to angle data and compute the theoretical tangent/cotangent value based on (1).

*Step 4. *Take (2) to conduct the response evaluation on the sensor.

*Step 5. *Modify the sensing response with (6) to find a minimized difference between the theoretical and optimal structure.

##### 2.3. GA

Intelligent algorithms are employed for strategy facilitating and cost reducing [17, 18]. Genetic algorithm, which is a class of optimization heuristics based on natural selection and genetic evolution, can be applied to problems where the variables to be optimized (genes) can be encoded to form a string (chromosome). Each string represents a trial solution of the problem. The GA operators exchange information between the strings to evolve new and better solutions [19]. Because the structure of magnetic-grating-like stroke-sensing cylinder determines the working performance, this paper takes GA as the optimization model for variable design.

A multiparameter optimization problem using genetic algorithm is formulated by defining the system figure and the objective function of response quality, respectively. The modal parameters (i.e., natural modes) of the real structure are numerically obtained using its finite element model, from which the reduced numerical modes correspond to the target structure. The original related parameters are discussed in Section 3, where the characteristics of working elements are also described. The algorithm we use is implemented as follows: (1)Initialization: consider a population of individuals, and each individual is described by one chromosome with genes, where each gene has alleles. Therefore, it has possible individuals [20]. Actual parameter information is contained within the components. The initial parameter values are randomly assigned.(2)Fitness calculation: the optimum configuration is sought by minimizing the value of signal response quality. The fitness function of each individual in current population can be calculated according to (6). The fitness of total population is evaluated based on the outcome from the proposed function.(3)Selection: the selection mechanism allows a small portion of chromosomes from the population for further evaluation [21]. Selection is accomplished by holding a fitness function among random competitors. The winner of the fitness function is the individual with the highest fitness, and the winner is then inserted into the mating pool [22]. In other words, chromosomes with higher fitness values are more likely selected to produce the next generation.(4)Crossover: offspring is generated by two randomly selected parents exchanging genetic information with each other. The crossover step aims at leading to minor differences between parents and children, strengthening the exploitative power of reproduction [23].(5)Mutation: mutation provides an opportunity to search new areas of the solution space. By randomly altering the alleles of genes, the GAs can effectively avoid trap situations and maintain sufficient variance in the population [24, 25], whereas the probability of mutation is around 1%, since a high probability of mutation reduces the GA to a random search function.(6)Replace the initial population with new individuals of best chromosome. An outcome value closer to the target value via repetitive iterative process can be generated. Consequently, after a large number of iterations, the best chromosome in the population is translated as the selected solution.

#### 3. Parameter Optimization

A sensor’s performance characteristics are governed by its specific architecture. The parameters with large uncertainty in the design process can affect the measurement accuracy significantly. Aiming at investigating the structure design for increasing working performance, optimization techniques are used to determine the design parameters. According to paragraph II, the key of measuring is the response of the sensor output. That is, in other words, given a parameterized representation of a sensor, we are capable of predicting its measurement accuracy.

To assist the optimization process, there are seven variables involved in this system (Table 1), which are diameter of permanent magnet (*S*_{MD}), height of permanent magnet (*S*_{HM}), distance from bottom of permanent magnet to centre of sensor (*S*_{DM}), distance from surface of piston rod to centre of sensor (*S*_{DC}), slot shape of piston rod (*f*(*y _{t}*)), diameter of piston rod (

*S*

_{CYD}), and groove width of piston rod (

*S*

_{PT}). Based on the working principle, a novel magnetic-grating-like stroke-sensing cylinder is designed and its structural parameters are illustrated in Figure 6. The magnetic-grating-like stroke-sensing cylinder is for detecting the displacement of piston rod under working circumstance. Careful design of the structural variables can substantially improve measurement accuracy. Hence, the optimal structure design issue is formulated as an optimization problem. All the following parameters are transformed from the modeling and simulation methods into the function value of .

We notice that the diameter of piston rod in a suspended state, which avoids the direct contact between the magnet and the Hall sensor [26, 27], achieves the purpose of design without optimization. So is set as defined value of 50 mm. Due to the definition of measurement accuracy, the groove width of piston rod represents one working period for segmentation. Hence, the smaller the groove width is, the higher accuracy can be obtained. Since the highest cutting quality is 1 mm, is given a balance value of 1.25 mm for machining.

The conceptual framework for optimizing a candidate design is shown in Figure 7. The optimal parameter value is searched using the signal response evaluation function, and the performance testing is deployed for model verification. The details of the signal response evaluation-based processes are described in what follows.

##### 3.1. Effect of Distance from Surface of Piston Rod to Centre of Sensor

To validate every analytical result further, other parameters are assumed fixed during the optimizing process. The sensor exists in the position between the permanent magnet and the piston rod, and the distance is named as . Typical parameters are quantified by selecting the proper value of ordinary used sensing systems, where , , , , and is the shape of equally spaced rectangular, which is further optimized for showing that how these values are suitable for design purposes. The variation of is typically taken from 0 to 5 mm with the step of 0.1 mm. Figure 8 shows the response curve of the Hall sensor with different values of by simulation.

It can be seen that the signal intensity is affected by the distance between the Hall sensor and the piston rod. As seen here, a smaller distance has a stronger sensor response while the value over 0.7 mm changes only slightly the strength of sensor response. However, results in the severe signal distortion. To evaluate the quality of signal response, Figure 9 shows the derivation results depending on the variation of with Equation (6) applied. The parameter value from 0.1 mm to 0.6 mm is shown to be in the desired position.

##### 3.2. Effect of Diameter of Permanent Magnet

The variation of permanent magnet shape has a direct effect on the magnetic induction, which is involved with the diameter and the height . Thereby, both of them can be analyzed by this method. We take according to the aforementioned outcome while other parameters are set as before. The signal response of variable parameter is shown in Figure 10. At the beginning, the sensor response is small when is small. The sensing signal response increases with the increasing of diameter and reaches a maximum value when the parameter is in the range 16–20 mm, but the maximum value is not influenced very much.

The signal response quality is found to be improved and stabilized as diameter increases to over 10 mm as shown in Figure 11. Considering the sensing signal strength outcome, the optimal value is taken.

##### 3.3. Effect of Height of Permanent Magnet

For the height of permanent magnet , a similar analysis result is yielded as shown in Figure 12. In the range of 0–30 mm, the output of Hall sensor increases while the difference in the response strength between 30 mm and 40 mm is not distinguished.

Figure 13 depicts the relation between signal response quality and the height of permanent magnet. The quality of sensing response fluctuates significantly in the range of 0–5 mm while slightly in 5 mm–35 mm. Particularly, it decreases slightly when . In line with actual requirement of high accuracy, the optimal set at the value of 35 mm can be suggested.

##### 3.4. Effect of Distance from Bottom of Permanent Magnet to Surface of Piston Rod

The mechanical position of permanent magnet plays an important role in changing the sensing signal as well. The physical distance between the permanent magnet and the piston rod is taken as typically designed parameter. For the typically selected parameters optimized above, the effect of on performance, with variation step of 0.5 mm, can be observed in Figure 14.

As seen here, the decreasing distance causes the increasing of signal response. However, the computing results of parameter values illustrate that the characteristic of piston rod can be clearly described at all values. In the practical design, different structures correspond to different locations of permanent magnet. The objective function curve of signal response quality is shown in Figure 15. When the distance increases (≥5.5 mm), the function value changes much greater than that of short distance, which indicates that a small value of is more appropriate for practical use. The optimal is given as 5 mm.

##### 3.5. Effect of Slot Shape of Piston Rod

The shape of piston rod slot is optimized by setting up the coordinate system of in Figure 16. The slot is described in plane *X*. Axis *Y* is the axle centre of piston rod, and axis *Z* is the start of periodic variation in slot shape. For a specific architecture, however, it is possible to define the best slot which is defined with the optimum “set” of the geometrical parameters [28]. Therefore, nine anchor points () from to are picked up for analysis. For simplicity, the anchor points are evenly distributed between the dot and . Considering the machining process, the values from to decrease monotonically. The optimization problem is formulated based on the design of these nine parameters. Constraints in (7) are then formulated to reject infeasible slot geometries.
where and is the variable quantity of *Z*-direction. Initially, the relationship of can also be written as
whereas the variation of these parameters is not constrained but depends on the optimization outcome. Since the values of , , and are defined earlier, the sensor response determined by optimal parameter set can be expressed as
where is the response and is the specific working moment of Hall element. Aiming at acquiring a harmonic signal, the response evaluation function is introduced as the objective function.

In the relationship, stands for the tangent/cotangent of response incorporated which can be derived as where is the signal outcome period. Thus, the shape optimization of the slot is transferred into detecting the minimum value of .

A GA is employed to search the optimum solution of . A flowchart that describes the GA optimization procedure that has been adapted in this paper is shown in Figure 17. The input to GA is , and the output is the geometry of the slot.

The process starts with randomly generated population. Each chromosome represents a possible set of parameters associated with a fitness value. The total population of each generation is evaluated based on slot geometry where GA carries out a fitness-based selection. In our GA, the size of the population is chosen to be 40 on the basic of modeling and modification. Processes are subject to recombination to form a successor population. The single-point crossover is used, and the probability of mutation is set as 0.17%. During recombination, chromosomes with the higher fitness values are more likely to be selected from the population, which passed on to the successor population. This is an iterative process, through which various generations are evolved until some stopping criterion is specified, and an optimal solution is reached based on an increase. The GA search process for the quality evaluation function is shown in Figure 18. At the beginning of evolution process (about less than 40 generations), the differences between every individual are relatively large, which results in individuals with low fitness being eliminated and ones with higher fitness being saved. After certain generations (about larger than 100 generations), the variation in target value is small, finally approaching to a constant value of 0.04.

Accordingly, a piston rod with optimized variables is obtained. Configuration uses the following parameter values: Figure 19 shows the photograph of the manufactured piston rod as well as its actual signal response with optimal pattern. Some signal fluctuation appears randomly within the sine signal response period due to the machining error, which can be eliminated by improving the accuracy of machining technology.

**(a)**

**(b)**

#### 4. Experiment and Analysis

##### 4.1. Experimental System

The experiments are conducted to verify the working performance of magnetic-grating-like stroke-sensing cylinder. The hardware used to perform the experiments is shown in Figure 20. Together, the linear slide and the fixed bearing house the piston rod and the permanent magnet, respectively, on the supporting base. A high-accuracy grating is placed parallel with the piston rod for precisely detecting the displacement. The sensor output is recorded with the comparison of the grating. During the test, the permanent magnetic is fixed while the piston rod is in reciprocating linear motion. As explained in Section 1, the displacement is then calculated from sensing signal of Hall sensor. Experiments are carried out in an air-conditioned room maintained almost at constant temperature, humidity, and pressure. According to the JJG 644-2003 “Verification Regulation of Vibration Displacement Transducer,” much work is devoted to assess the static and dynamic performance.

##### 4.2. Static Testing Configuration

###### 4.2.1. Stroke Test

Stroke tests contain short-term test and long-term test. Both are basically designed for testing the error between theoretical stroke and experimental stroke. The working distance of short-term test is one grating, which is one thread lead of the piston rod. The output of Hall sensor is detected every 5°, and the subdivision points is 96. Figure 21 shows the results of short-term test, where is the sampled data point and is the actual one. The maximum drift of short-term test is .

Similarly, the long-term test takes 100 continuous grating as the working distance. The maximum error of long-term is . The distribution of all the recorded points is presented in Figure 22, which shows the range of output error in forward and backward stroke. The numerical difference is 0.045 mm in maximum.

###### 4.2.2. Stability Test

The magnetic flux density is addressed relying on highly stabilized system output. In the power supply driving mode, the sensor is of different position of the stroke. The pass/fail of stability is based on the output of the sensor. The voltage is kept at 24.0 V for 1 hour by use of highly stabilized power supply, and the displacement is measured every 2 hours. The corresponding displacements , obtained after sensing steps, are given in Table 2.

Arithmetic average difference can be calculated during relevant periods. The value of magnetic flux density stability error is 0.0195 mm in the power supply driving mode.

###### 4.2.3. Reproducibility Test

The reproducibility test is designed to identify the repetitive feature of sensor output within the same stroke. The calibration is conducted three times while the piston rod repeats in a defined path. The recorded values, for each measurement point, are shown in Table 3 and indicated, at each point, the maximum difference within all test results characterizing the reproducibility error which can be written as where and represent the measurement data set of these calibration tests.

###### 4.2.4. Static Error

The absolute static error expresses the difference between the real and ideal characteristics, which has the form as follows. where , , , , and are errors caused by short-term variations, long-term variations, sensing instability, data agreement uncertainty, and variance of machining accuracy, respectively. The actual length of piston rod is , and the parallelism error is according to machining technology. Hence, the measurement error in relationship with mechanical accuracy is

The error determined by machining deviation is , which is however an insignificant contributing factor. Thus, the static error of magnetic-grating-like stroke-sensing cylinder is .

##### 4.3. Dynamic Testing Configuration

###### 4.3.1. Dynamic Error

In this research, the dynamic test deals with the working accuracy of different measuring modes. The setup for evaluating the dynamic error is as follows: a low-speed (0.1 m/s) and a high-speed (0.3 m/s) measurement stroke are carried out three times with the sampling rate of 100 Ksps.

Dynamic error is evaluated in accordance with the definition in [29]. The target position is where identifies the particular position among other selected target positions along the stroke. is the measured position reached by the functional point on the *j*th approach to the *i*th target position. Then, the positioning deviation is
and the mean positioning deviation is

The standard uncertainty of the positioning deviations obtained by a series of approaches at a position is

Taking a coverage factor , the positioning repeatability at a position is

We take the maximum value of the positioning repeatability at any position as the sensor repeatability:

The difference between the algebraic maximum and minimum of the positioning deviations at any position is the systematic positioning error , which can be written as

Now, the coverage factor , the sensor positioning accuracy derived from the combination of systematic positioning errors and the estimator for positioning repeatability, is

The standard deviation of measured values is

In the working displacement of 0–90 mm, recordings corresponding to speed and position variation are shown in Figure 23.

**(a)**

**(b)**

**(c)**

**(d)**

The position deviation of proposed sensing stroke is within the range of , and the error is found to be stable and repeatable. Based on the equations above, the positioning accuracy of different working distances is deduced (Figure 24).

As shown in Figure 24, different position tracking results in the variation of positioning accuracy. The working errors at any point involved with stroke sensor are and . It can also be obtained that the sensor positioning accuracy is and standard deviation of measured values is .

#### 5. Conclusion

In this paper, the devising and deploying of magnetic-grating-like stroke-sensing cylinder are specifically studied, according to its working mechanism and distinctive structure. Motivated by the significance of developing a high-accuracy displacement sensing device, capable of operating in different working distances and optimizing its performance, we set out to define the signal response quality of working parameters. The signal response quality evaluation algorithm is constructed to characterize sensor response in developing a stroke-sensing cylinder, which leads to the evaluation of structure variables mathematically. The optimal displacement sensor combines high response and state-of-the-art signal segmentation method.

The prototype is processed and tested. The proposed testing workbench, together with the magnetic-grating-like stroke-sensing cylinder, is set up in laboratory. Based on the static and dynamic testing, experimental outputs are recorded and demonstrate that it has a stable, repeatable performance and high measurement accuracy for a wide range of stroke. The testing results characterize the magnetic-grating-like stroke-sensing cylinder, which can replace current measurement devices and fulfill the demand of high-accuracy displacement detection.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.

#### Acknowledgments

This research is a general project supported by Shanxi Province Science Foundation for Youths (2015021123) and Shanxi Provincial Key Laboratory Open Fund (XJZZ201605).