Research on Manufacturing Equipment Operation State Evaluation Technology Based on Fractional Calculus

1. Introduction

Mechanical manufacturing equipment is an indispensable tool in industrial production, and the working reliability of manufacturing equipment is a prerequisite for the normal operation of an enterprise's manufacturing system. Therefore, realizing a real-time and accurate assessment of the operational status of manufacturing equipment has become an urgent problem in enterprise production. If we can achieve real-time accurate assessment of the equipment operating status, we can scientifically develop equipment maintenance programs and equipment failure points in advance to predict and prevent, ensure the reliability of the enterprise manufacturing equipment management system and the normal operation of enterprise production, and promote the development of the manufacturing industry is of great significance.

According to the operation rule and life cycle of manufacturing equipment, the operation of manufacturing equipment can be divided into three stages: normal operation, fault warning, and fault state. The evaluation of the equipment operation state at home and abroad mainly focuses on the diagnosis stage of a fault. For example, Andrei S. Maliuk et al. proposed a Gaussian mixture model-based (GMM) bearing fault band selection (GMM-WBBS) method for signal processing, which benefits reliable feature extraction using a fault frequency oriented Gaussian mixture model (GMM) window series [1]. Stator winding faults, broken rotor bar faults, rotor asymmetry, and abnormal air-gap eccentricity have been applied to the diagnosis of bearing faults in many studies [2,3,4]. Krzysztof Kecik et al. [5]. presented the problem of rolling bearing fault diagnosis based on the vibration velocity signal. To determine the optimal set of fault-relevant variables, Deng et al. [6]. proposed a fault reconstruction algorithm based on least angle regression (LARS). Referance [7,8] used a contribution plot algorithm to identify the fault-relevant variable without any prior information about the fault. Wang et al. applied the ADD-HMM algorithm to predict the occurrence of mechanical equipment faults [9]. Reng et al. implemented a deep learning waveform image recognition method for the fault diagnosis technology of bearings [10]. Yong et al. applied a Bayesian network for the fault diagnosis of weapons and equipment [11]. Bin et al. applied the variability of manic power to determine the failure time of mechanical equipment [12]. Qi et al. [13]. proposed a progressive fault diagnosis method(PFDM) for the overall diagnosis of an entire EHA system, greatly improving the safety and reliability of the double-redundancy EHA system.

It can be seen from the existing research that they have the following characteristics.

1)

Limitations of application objects: Each fault diagnosis method can only be applied to specific research objects, and cannot be extended to the diagnosis of various faults in other manufacturing equipment.

2)

Limitations of application functions: In equipment management, there are no real-time online monitoring and fault prediction functions, which affect the reliability of equipment operation and the continuity of production system work.

Therefore, there is an urgent need for a method to monitor and evaluate the running status of various online equipment in current manufacturing systems. This method should have the function of equipment operation state judgment and fault prediction in advance, and ensure the normal operation of equipment and production systems through timely troubleshooting of equipment faults. This study draws on the research results of various experts and scholars in the evaluation of the operation status of different types of equipment. To analyze the working characteristics of various types of manufacturing equipment, explore evaluation methods that can be adapted to the operating status of multiple types of manufacturing equipment, combine the advantages of IOT technology and data processing technology, and achieve the research goal of real-time and accurate evaluation of the operating status of various types of manufacturing equipment.

2. Manufacturing equipment operating condition assessment methods

2.1. fundamental principle

Manufacturing equipment is a combination of multiple motion units, each of which has its own motion law and motion parameters. Many motion units exist as long as there is a unit in a state of failure, and the equipment will be in a state of failure. Therefore, the operational status of the manufacturing equipment depends on the shortest life cycle of all motion units. As part of the manufacturing equipment, the working state of the motion unit is also divided into three stages: normal operation, fault warning, and fault state. The motion parameters under each working state have their own value intervals, which are [0,C₁], [C₁,C₂], and [C₂,C₃], respectively. If we can obtain accurate operating parameters C of the motion unit in real time, we can accurately judge the operating state of the motion unit or even the entire equipment in real time by referring to the value intervals of each working state.

2.2. Implementation steps of the method

According to the evaluation principle of the operation state of the manufacturing equipment, the evaluation of the operation state of the entire equipment can be divided into the following steps:

1)

The main motion units that affect the running state of the equipment were analyzed according to the type of manufacturing equipment.

2)

The operation law of the main motion units was analyzed, and the motion unit with the shortest life cycle was selected.

3)

The main parameters affecting the running state of the motion unit and their value ranges under different working states were analyzed.

4)

Obtain high-precision detection data for important motion parameters when the motion unit works normally in real time.

5)

The running state of the motion unit is determined by comparing the high-precision parameter values obtained in real time with the above value range.

2.3. Premise of method implementation

According to the introduction of the basic principle and implementation scheme of the evaluation method of manufacturing equipment operation status, it can be seen that the implementation of the scheme needs to solve the following three problems.

1)

The judgment technology of the shortest life cycle motion unit among many motion units;

2)

Real-time detection technology for the main parameters that affect the running state of the motion unit.

3)

High precision acquisition technology of main detection information data in complex environment.

Among the above three preconditions, the judgment technology of the shortest life cycle motion unit in the motion unit can realize the judgment function of the minimum life cycle motion unit of different manufacturing equipment through the analysis of the motion characteristics of the motion unit and the accumulation of equipment management experience. In the era of rapid development of information network technology, there is no technical problem in using the Internet of Things to obtain all kinds of detection information and data in real time; however, in the collection of detection information data, the detection value of information will inevitably be affected by factors such as equipment performance, working environment, and signal interference, which will cause detection errors of information data, resulting in the acquisition technology of high-precision information data becoming the main problem in the implementation of manufacturing equipment operation status evaluation schemes.

3. High-precision information data acquisition technology

The essence of the hardware method is to improve the detection accuracy of the information data using high-performance detection instruments. The main research results for the hardware methods are as follows: Hu et al. [14] proposed a high-precision safety valve test architecture with three testing channels and effectively solved the problems of current safety valve testing. Huijun et al. [15] developed a comprehensive sliding-separation test platform for RV reducers to realize high-precision and high-display test performance for various RV reducer parameters. Reference [16,17] proposed a load differential radiation pulse on a transient electromagnetic high-performance radiation source for pulse-scanning detection to solve the problems of urban electromagnetic interference and insufficient harmonic components emitted by the radiation sources. Reference [18,19,20] designed a hardware system based on radar and realized a real-time detection function of underground space-related information by enlarging the detection information. Jiaqi et al. [21] proposed a one-stage remote sensing image object detection model: a multi-feature information complementary detector (MFICDet), which can improve the ability of the model to recognize long-distance dependent information and establish spatial–location relationships between features.

However, in engineering applications, we found that the hardware method had the following shortcomings.

1)

The detection accuracy of information data depends on the performance of the detection equipment. With improvements in detection accuracy, the cost of the detection system is higher. Therefore, they have low-cost performance.

2)

Their essence is to reduce the signal distortion caused by energy loss and signal interference in the information transmission process by improving signal strength. However, when collecting information, the measurement error of the information data cannot be eliminated owing to the differences in equipment performance and working environment.

In recent years, most researchers have attempted to use software methods to achieve high-precision information data detection to solve the shortage of hardware methods for information data detection in complex environments. The essence of the software method is the information data-fusion algorithm. up to now, there have been many research results Common mathematical algorithms are fuzzy set theory [22], fuzzy neural networks [23], probability model [24] and particle swarm optimization algorithm [25], et al. and obtained a regrettable research review. For example, Huo et al. [26] proposed an integral infinite log-ratio algorithm (IILRA) and an integral infinity log-ratio algorithm based on signal-to-noise ratio (BSNR-IILRA) to improve the detection accuracy of the laser communication detection position in the atmosphere. Zhiyuan et al. [27] proposed a normalized-variance-detection method based on compression sensing measurements of received signals and solved the problem of fast and accurate spectrum sensing technology under the condition of a low signal-to-noise ratio. Liu et al. [28] proposed a target detection algorithm based on the improved RetinaNet, which is suitable for transmission-line defect detection and improved the intelligent detection accuracy of UAV in power systems. Cheng et al. [29] proposed a lightweight ECA-YOLOX-Tiny model by embedding an efficient channel attention (ECA) module into the lightweight YOLOX-Tiny model, which has a higher response rate for decision areas and special backgrounds, such as overlapping small target insulators, insulators obscured by tower poles, or insulators with high-similarity backgrounds. Liu Wenqiang et al. [30] introduced a point cloud segmentation and recognition method based on three-dimensional convolutional neural networks (3-D CNNs) to determine the different components of the catenary cantilever devices. Yin et al. [31] proposed a complementary symmetric geometry-free (CSGF) method that makes the detection of cycle slips more comprehensive and accurate. Lingfeng et al. [32] established a junction temperature model based on a multiple linear stepwise regression algorithm and used it to extract high-precision intersection online temperatures. However, through the analysis of various current software methods, the following deficiencies were found in the detection of information in complex environments:

1)

However, they do not improve the detection information strength and cannot solve the problems of energy loss and signal interference during information transmission. Therefore, it is difficult to apply this method in engineering practice.

2)

They did not analyze the cause of the information data detection error, the change rule of each influencing factor, or its influence on the detection value. Therefore, it is difficult to improve the detection accuracy of information data by reducing the detection error caused by various influencing factors.

Therefore, to date, an ideal high-precision detection method for information data under the joint action of multiple influencing factors has not been developed. To solve these problems, our team has been using the method of fractional calculus theory in data processing for many years [33,34,35,36,37,38,39] and found that fractional differential operators are suitable for studying nonlinear, non-causal, and non-stationary signals, and have dual functions of improving detection information and enhancing signal strength. Therefore, by fusing the differences between information and data, the information and data detection errors caused by various influencing factors can be eliminated. By improving the signal strength of the information, it can compensate for the energy loss of the signal in the transmission process and improve the anti-interference ability of the signal.Therefore, in this paper, we try to combine fractional order calculus with IoT technology to realize the function of real-time accurate evaluation of the operation status of manufacturing assembly.

4. Fractional order differentiation

4.1. Definition of fractional order differentiation

Fractional differentials, also known as fractional derivatives, extend the differential order of an integer-order differential equation to a fractional order. The concept of fractional differential emerged in 1812. In hundreds of years of development of this theory, many scholars have proposed their own definition methods and theoretical systems based on their own understanding and application fields. Therefore, a strict definition of the fractional differential is not available to date. At present, the most commonly used definitions are those of Grunwald–Letnikov (G-L), Caputo, and Riemann–Liouville (R-L).

4.1.1. R-L definition of fractional differential the fractional

According to the principle of mutual inverse operation of Cauchy's indefinite integral formula with fractional differentiation and fractional integration,, we can obtain the R-L definition of fractional differentiation:

$${}_a{}^{ }D_t^{v}f\left(t\right)=\frac{d^n}{d{t}^n}\left[{}_{a}D^{v-n}f\left(t\right)\right]=\frac{1}{\mathrm{\Gamma }\left(n-v\right)}\frac{d^n}{d{t}^n}{\displaystyle {\int }_a^t{\left(t-\tau \right)}^{n-v-1}f\left(\tau \right)d\tau }$$

(1)

In the equation above, 0 ≤ n-1 < n., Г Is the gamma function.

The main advantage of this definition is that the initial value of the Laplace transform can be obtained only by using an integer-order derivative. Its disadvantage is that it has stricter requirements for the function f (x) than other definitions, but its premise is that the integer-order derivative of the function f (x) is absolutely integrable.

4.1.2. G-L definition of fractional differential the fractional

For any real number v, where the integer part of v is denoted as [v], assuming that the function $f\left(t\right)$ has (n + 1)-order continuous derivatives within the interval [a, t], when v > 0 and n ≥ [v], the fractional v-order derivative is defined as follows:

$${}_a{}^{ }D_t^{v}f\left(t\right)=\underset{h\to 0}{\lim }{f}_h^{\left(v\right)}\left(t\right)=\underset{h\to 0}{\lim }{h}^{-v}{\displaystyle \sum _{j=0}^{\left[\frac{t-a}{h}\right]}{\left(-1\right)}^j\left(\begin{array}{c}v\\ j\end{array}\right)f\left(t-jh\right)}$$

(2)

In the above equation, 0 ≤ n-1 < n, $\left(\begin{array}{c}v\\ j\end{array}\right)$ is a binomial coefficient, and its formula is

$$\left(\begin{array}{c}v\\ j\end{array}\right)=\frac{v\left(v-1\right)\left(v-2\right)\mathrm{....}\left(v-j+1\right)}{j!}=\left[\begin{array}{c}-v\\ j\end{array}\right]$$

(3)

Based on the classical definition of an integer derivative of a continuous function, this definition extends the order of the differential from integer to fraction, which is suitable for numerical calculations.

4.1.3. Caputo definition of fractional differential the fractional

$${}_a{}^{C}D_t^{v}f\left(t\right)={}_{a}D_t^{-\left(n-v\right)}\left[D^{n}f\left(x\right)\right]$$

(4)

In the equation above, 0 ≤ n-1 < n

To simplify the calculation of the fractional differential, this definition further improves the Grunwald Letnikov definition on the basis of the basic properties of fractional calculus.

Among the three definitions, the G-L definition is widely used in engineering because of its simple calculation process and high speed. Therefore, the G-L definition of fractional partial differentials was used in this study to research the manufacturing mquipment mvaluation technology.

4.2. Properties and applications of fractional order differentiation

Detection data were obtained from the analysis and processing of the detection information. Therefore, the measurement error generated during the acquisition and transmission of detection information directly affects the accuracy of the detection data. In terms of the mathematical properties of the signal and the characteristics of the signal structure, the detection signal contains information on the fractional differential characteristics; however, this type of information is not suitable for processing by integer differential operators. For years, many mathematicians and scientists have attempted to apply fractional calculus theory in signal processing and have achieved good application effects [40,41]. It is assumed that there is a detection signal $S\left(t\right)$, where $S\left(t\right)\in L^2\left(R\right)$, and its Fourier transform is.

$$S\left(\widehat{\omega }\right)={\displaystyle \underset{R}{\int }S\left(t\right)}{e}^{-i\omega t}dt$$

(5)

Let $S^v\left(t\right)$ be the v-order differential of $S\left(t\right)$. According to the properties of the Fourier transform, we know that the v-order differential operator is equal to the multiplicative operator of the v-order differential multiplier function $\widehat{d}(\omega )={(i\omega )}^v$. Thus, we can obtain the following equation:

$$D^{v}S\left(t\right)\stackrel{FT}{\iff }{\left(\widehat{D}S\right)}^v\left(\omega \right)={\left(i\omega \right)}^{v}S\left(\widehat{\omega }\right)={\left|\omega \right|}^v{e}^{i{\theta }^v\left(\omega \right)}S\left(\widehat{\omega }\right)$$

(6)

From the perspective of signal modulation, the physical meaning of the fractional differential of a detection signal is equivalent to generalized amplitude and phase modulation. From the perspective of signal processing, the v-order fractional calculus operation of the detected signal is equivalent to establishing a linear time-invariant filtering system for the signal, and its filtering function is:

$$d^v(\widehat{\omega })={\left[\omega \right]}^v{e}^{i{\theta }_v\left(\omega \right)}$$

(7)

Through the Fourier transform defined by the G-L of fractional-order calculus, it is known that the essence of the fractional-order calculus processing of the energy signal is to filter the signal, and its filter function is $j^v\left(\widehat{\omega }\right)={\left(i\omega \right)}^v={\left|\omega \right|}^v{e}^{i{\theta }_v\left(\omega \right)}$. According to the filter function, we can draw the spectral characteristic curves of the fractional-order differential operator and fractional-order integral operator in signal processing, as shown in Figure 1. After analyzing the spectral characteristic curves, we can obtain the following characteristics of fractional-order calculus in signal processing.

From Equation (7), we know that after fractional differential processing, the signal has the following characteristics:

1)

The signal shows different levels of signal enhancement for different fractional differential operators so that the very-low-frequency components of the signal can be preserved nonlinearly.

2)

From a physics perspective, signal processing by the fractional differential operator can be understood as the generalized amplitude phase modulation of the signal. Thus, the fractional differential operator can significantly improve signal strength when processing the high-frequency part of the signal.

3)

The fractional differential operator significantly improves the high-frequency signal strength. It also improves low-frequency signals.

5. Fractional order differentiation based fusion model of equipment operation parameters

5.1. Characteristics of manufacturing equipment inspection information data

Manufacturing equipment is divided into continuous manufacturing equipment and discrete manufacturing equipment, which is mainly used in the production mode of a small variety of large quantities. In the production process, the operating state of the manufacturing equipment movement unit is a regular change, as long as the law of change in the operating state and the operating state of the evaluation criteria, the evaluation of the operating state of the continuous manufacturing equipment function can be achieved. Discrete manufacturing is the main mode of production in China's manufacturing industry, production characteristics of multi-variety small batch production, the differences in product types lead to the differences in the manufacturing process of the manufacturing equipment used, as well as the differences in the operating state of the manufacturing equipment work. The dynamic changes in the parameter values of the motion units of discrete manufacturing equipment make it impossible to apply the continuous manufacturing equipment operating state assessment scheme to the evaluation of the operating state of discrete manufacturing equipment. This paper addresses the important problem that the operating states of manufacturing equipment in two manufacturing types are very different and difficult to consider by traditional evaluation techniques, as well as a way to realize the evaluation function of the operating states of manufacturing equipment in two manufacturing modes by comparing accurate motion parameters and evaluation criteria collected in real time without considering the motion law of manufacturing equipment.

5.2. Fractional order differential based operational state evaluation model

According to the spectral characteristics of fractional-order calculus in information processing, it is known that both fractional-order differential and integral operators have the advantage of reducing the measurement error caused by external interference signals, and both can effectively fuse to remove the detection error of information data caused by a single influence factor in a one-dimensional space. Therefore, fractional-order differential or integral arithmetic can be applied to the processing of the main parameter values of the manufacturing equipment motion unit, and fractional-order differential arithmetic has the advantage of enhancing the signal strength and reliability of the detection system compared to fractional-order integral arithmetic, which has the disadvantage of weakening the information strength. Therefore, fractional-order differential or integral arithmetic can be applied to the processing of the main parameter values of the manufacturing equipment motion unit, and fractional-order differential arithmetic has the advantage of enhancing the signal strength and reliability of the detection system compared to fractional-order integral arithmetic, which has the disadvantage of weakening the information strength.

After analyzing the common fault types of manufacturing equipment, we can obtain the main fault types and their main causes, analyze the motion unit where the fault point is located, and obtain the main influencing factor x of its operation status; then, the value interval [a,b] of the influencing factor x and its corresponding detection value S_i in the interval after the experiment, fit the function formula S(x) between the influencing factor x and the detection value S Then, we can obtain the mathematical model of the manufacturing equipment operation status evaluation shown in Equation (8) according to Equation (2).

$$\begin{array}{c}{}_a{}^{ }D_x^{v}S\left(x\right)=\underset{h\to 0}{\lim }{S}_h^{\left(v\right)}\left(x\right)=\underset{h\to 0}{\lim }{h}^{-v}{\displaystyle \sum _{i=0}^{\left[\frac{b-a}{h}\right]}{\left(-1\right)}^i\left(\begin{array}{c}v\\ i\end{array}\right)S\left(x-ih\right)}\approx S\left(x\right)-vS\left(x-h\right)\\ -\frac{v\left(1-v\right)}{2!}S\left(x-2h\right)+\cdots +\frac{v\left(1-v\right)(2-v)\cdots (n-1-v)}{n!}S\left(x-nh\right)\end{array}$$

(8)

6. Application of assessment methods

6.1. Experimental platform construction

According to the analysis of manufacturing equipment failure types and their causes over the years, the common failure of manufacturing equipment from the spindle motion unit cannot meet the processing accuracy requirements. After analyzing the causes of various machine tool spindle unit failures, it is concluded that the main reason is that the fatigue damage of the spindle bearing brings a bearing clearance greater than the maximum allowable value, resulting in machine tool machining accuracy below the machining accuracy standards. Most bearings used in machine tool spindle components are rolling bearings, which have the performance characteristics of easy wear and short life cycle in the application, thus becoming the main source of spindle motion unit failure. The causes of the short bearing life cycle are fatigue damage and permanent deformation; however, fatigue damage is the main cause of bearing failure, which is the result of its accumulation due to the amount of wear during normal operation. Therefore, only by obtaining the real wear amount of spindle bearings in real time and referring to the value interval of each operating condition can the evaluation function of the operating condition of the manufacturing equipment be obtained.

The bearings used in the spindle motion unit of the equipment were the NN3046 series cylindrical roller bearings. As the equipment has been used in the production of multi-species and small-lot products, the wear of the bearings has irregular and dynamic changes, which makes the traditional method of assessing the operating status of the equipment by applying the regularity of changes in motion parameters seem The traditional method of assessing the operating condition of the equipment by applying the regularity of the motion parameters appears to be inefficient. Therefore, this study attempts to verify the feasibility and reliability of the method described in this paper for evaluating the operating conditions of the equipment.

6.2. Information Data Collection

6.2.1. Data collection methods

The best solution for measuring the parameter value of the object is to place the corresponding sensor at the corresponding position of the object; however, in the detection of rolling bearing wear of precision machine tools, the existence of disassembly of the main moving parts of the machine tool will lead to a decrease in accuracy. The detection value of the clearance between the outer ring of the rolling bearing and the spindle is affected by roller interference and other factors, resulting in the above method in the detection of spindle bearing wear of precision machine tools being difficult to implement. Because the front bearing of the machining center is the component with the greatest force and wear in the main motion system, based on the analysis of the structure of the main drive system, the experiment used a displacement sensor arranged on the inside of the front bearing end cover of the spindle near the spindle and indirectly realized the detection function of the front bearing wear by detecting the value of the external runout of the spindle during operation. Because it is difficult to predict the magnitude and direction of the force on the spindle during machining and the speed is too fast to bring about the loss of experimental data owing to the lack of sensitivity of the testing instrument, the experiment was carried out under idling conditions with a spindle speed of 750r/min. When the spindle is idle, the clearance caused by bearing wear is reflected in the upper part of the bearing bore under the action of self-weight; therefore, the effective detection of bearing wear can be realized by arranging displacement sensors only in the upper part of the spindle. Based on a comprehensive consideration of the installation space and detection accuracy, five L3002-12.7 LVDT displacement sensors were placed in the upper part of the bearing end cover, and the configuration position of the testing instrument is shown in Figure 2.

To ensure real-time detection data, the experiment used IoT technology to transmit the information data collected by each sensor in real time. IoT technology has two methods for the transmission of detection information: wireless networks and wired networks. Wireless networks have the advantage of real-time collection of mobile equipment production information, but also have a shortage of information detection errors; wired networks are suitable for real-time collection of production information of fixed equipment in close proximity, and have the double advantage of real-time and accuracy in the collection of fixed equipment production information. Combining the advantages and disadvantages of the two transmission methods and the working characteristics of the equipment in this case, the experiment applied the wired network to transmit the real-time collected bearing clearance information to the information processing center for the expert system to analyze and judge.

6.2.2. Experimental data collection

In a limited area, the variability of similar testing information data mainly originates from the variability in the performance of the testing equipment. Although the standard deviation is an important indicator reflecting the performance of the testing equipment, the performance of each sensor changes dynamically owing to the influence of the service life and the working environment. To obtain the true standard deviation of the sensors during the test, the equipment spindle at 750r/min idle speed, the application of each sensor in the C point position to collect test data, the spindle every two revolutions sampled a total of six times, and the relevant test data and its standard deviation values are listed in Table 1.

After obtaining the standard deviation of each sensor, the spindle continues to rotate, and the five sensors are sampled six times at the detection points A, B, C, D and E according to the above sampling frequency, and the average of the six detection values is taken as the detection value of the sensors at each point, and the detection values of each sensor at different points are shown in Table 2.

6.3. Analysis and processing of detection data

Although the data shown in Table 2 can be intuitively seen, sampling point D is the location of the largest amount of wear, because the differences in the performance of the testing equipment lead to a large difference between the detection data values, resulting in a large error between the detection data of different sensors at the same sampling point, such as sampling point A, there are 2# and 4# sensor measurement errors of more than 5%; therefore, the data in the Table 1 are not easy to use in the evaluation of the wear amount. Therefore, the evaluation of the spindle bearing wear must be based on the high accuracy of the measured values at each sampling point to lay the foundation for the comparative method in the identification of the maximum wear point and to avoid misjudgment in other cases.

6.3.1. Selection of the influence factor of the detection value

According to the above analysis, information detection data in the collection process will be affected by a variety of factors such as equipment performance, working environment, and signal interference; however, whether it is mobile equipment or fixed equipment, its working interval is within a limited range, and the working environment and signal interference factors on the same detection value are basically the same, so the variability between the detection data mainly comes from the variability between the performance of the detection equipment. Considering that the standard deviation is the best parameter for measuring the performance of the testing equipment, the standard deviation Si of the sensor can be used as the influence factor x of the bearing clearance detection value. The influence of the influencing factor on the detection value of each testing point can be explored according to the standard deviation of each sensor and its measured value at each testing point, as shown in Table 2.

6.3.2. The functional relationship between the detection value F_i and the impact factor x_i

Among various fitting algorithms, the least squares method has the advantage of not requiring a priori data information in the data processing process. It is widely used to fit the function polynomial of the measurement data, and can obtain the ideal data fusion accuracy, which is suitable to be applied to fit the equation of the function F(x) between the detected value F_i and each influence factor x_i (standard deviation S_i) at each point in the experiment. Assuming that the expression of the function F(x) is:

$$F(x)=a_0+a_{1}x+a_2{x}^2+\cdots +a_n{x}^m$$

(9)

According to Equation (9), we can determine the parameters in the equation to fit the required functional relationship. According to the data values shown in Table 2, the Polyfit function in MATLAB software can be applied to fit the values of each sampling point shown in Equation (10), and the functional relationship between each influence factor x_i $F_a\left(x\right)、F_b\left(x\right)、F_c\left(x\right)$、$F_d\left(x\right)$, and $F_e\left(x\right)$

$$\left\{\begin{array}{l}{F}_a(x)=-3.6416x+0.1269\\ F_b(x)=1.0405x+0.1562\\ F_c(x)=-3.0058x+0.1874\\ F_d(x)=1.0116x+0.2027\\ F_e(x)=5.0x+0.1265\end{array}\right.$$

(10)

6.3.3. Selection of fractional order v and step size h values

As shown by Equation (8), establishing a detection data fusion model based on fractional differentiation under the definition of G-L requires two problems to be solved: determination of the fractional order v and selection of step h.

1)

Selection of order v:

After analyzing the amplitude-frequency characteristics of the fractional-order differential operator, it can be seen that when the differential order v ∈ [0,1], the signal intensity in the high-frequency stage increases with an increase in the fractional order. However, with an increase in the fractional order, the difference between the differential operators of different orders in the enhancement value of the signal shows a decreasing trend. Therefore, from the viewpoint of saving space and facilitating calculations, the experiment explores the application effect of the fractional-order differential operator in the processing of machining center spindle bearing wear detection data when the fractional order v is taken as 0.5, which is the middle value of [0,1].

2)

Selection of step h:

In the processing of information data, the signal frequency depends on the value interval [a, b] of the influencing factor and its step value h. According to the amplitude-frequency characteristics of the fractional-order differential operator, it is known that the step length and the value interval of the influence factor are associated with the frequency of the signal; however, in the high-frequency region of the signal, the difference between the signal enhancement effect of the differential operator in the same order is negligible. The smaller the step size h, the higher the fusion accuracy of the data, but it also causes the shortcoming of decreasing the fusion efficiency owing to the increase in the computation step n. Therefore, the fusion speed is considered in this case. Therefore, in this case, considering both the fusion speed and fusion accuracy, the fusion step h = 0.0001 is taken according to the value range [0.0037, 0.0048] of the influence factor x (sensor standard deviation) shown in Table 2, and the number of steps n = 11 needs to be calculated.

6.3.4. Information Data Processing Techniques Based on Fractional Order Differential Operators

According to the fusion model of manufacturing equipment operation parameters under the G-L definition shown in Equation (3), the fusion processing model of the machining center spindle bearing wear data can be expressed as

$$\begin{array}{c}{}_a{}^{ }D_x^v{F}_k\left(x\right)=\underset{h\to 0}{\lim }{h}^{-v}{\displaystyle \sum _{i=0}^{\left[\frac{b-a}{h}\right]}{\left(-1\right)}^i\left(\begin{array}{c}v\\ i\end{array}\right)F_k\left(x-ih\right)}\approx F_k\left(x\right)-v{F}_k\left(x-h\right)\\ -\frac{v\left(1-v\right)}{2!}{F}_k\left(x-2h\right)+\cdots +\frac{v\left(1-v\right)(2-v)\cdots (n-1-v)}{n!}{F}_k\left(x-nh\right)\end{array}$$

(11)

where: v = 0.5, n = 11, h = 0.0001,k = a, b, c, d, e.

In Equation (11), x is the standard deviation $S_i$of each sensor, and $F_k$ is its detection value at each position. Now, the parameters (v, h, n) are substituted into Equation (6), and combined with Equation (5), the detection data $F_{ij}$ of each sensor at different positions shown in Table 3 can be obtained using MATLAB software.

From the data shown in Table 3, it can be seen that the fused data values show two characteristics compared to the data before fusion processing: the measured values of different sensors at each detection point are evenly distributed around the mean value and the variability between the data has been significantly reduced; the fused values increase significantly, and the amplification factor K is as high as 17.05 compared to that before fusion. It was verified that the fractional-order differential operator has amplitude-frequency characteristics that enhance the information intensity and reduce the variability among the information data. To compare the accuracy of the data before and after the fusion of the machining center spindle bearing wear, the data shown in Table 3 were divided by the amplification factor K to obtain the final value of the data shown in Table 2 after the fusion of the 0.5 order differential operator, and the results are shown in Table 4.

6.4. Analysis of operating condition assessment results

6.4.1. Evaluation criteria for bearing wear

Machinery manufacturing equipment in the important parts of the life cycle is equipped with normal operation [0, T₁], failure warning [T₁, T₂], and failure state [T₂, T₃] in three time periods, corresponding to the main parameter values for [0, C₁], [C₁, C₂], and [C₂, C₃], and the state for standard parts and important parts of the main parameters have been developed corresponding to the value of the standard (parameter value C ≤ C₃). However, C₁ and C₂ parameter values have not been issued by the state corresponding to the value of the standard; each enterprise can only be determined according to the performance of the equipment and production characteristics. In this case, considering that the performance of the machining center is slightly higher than that of similar machine tools, and it is mainly used for processing small- and low-precision parts, with reference to the national standard of C₃ ≤ 0.3 mm for bearing clearance of diameter more than 100 mm, the parameters C₁ = 0.280 mm and C₂ = 0.295 mm are set in the experiment to judge the operation of the machining center under the existing conditions.

6.4.2. Evaluation of operation status

Combining the data shown in Table 2 and Table 4, it can be seen that the standard deviation between the fused data is significantly reduced compared to that before fusion. As shown at Detection Point D, the standard deviation of the fused detection data is 0.0013, and the precision is close to 5 times that before fusion, which is significantly lower than the 0.0055 of the least squares method and 0.0036 of the particle swarm optimization [25]. The data distribution diagram of detection point D before and after fusion is shown in Figure 3. It can be observed that the data are randomly distributed near its average value after being processed by the 0.5-order differential operator, and the discreteness is significantly reduced, greatly improving the accuracy of the detection data of the spindle bearing wear of the machining center. According to the data shown in Table 4, the test data shown at test point D were obviously higher than the measured values of the other test points. According to the evaluation method of the fault points mentioned above, test point D is the maximum point of the wear amount of the spindle bearing of the machining center in this case, and the wear amount is δ _D = 0.215 mm. With reference to the above set operation state evaluation standard, the wear amount δ _D is within [0, 0.28], indicating that the equipment is in normal operation.

6.5. Application analysis of experimental results

According to the preconditions for the implementation of the evaluation method of the manufacturing equipment operation state described above, this case, through the analysis of the operation characteristics of the current typical manufacturing equipment - machining center, we conclude that the wear value of the main shaft bearing is an important basis for evaluating the operation state of the machining center and realizes the judgment function of the shortest life cycle motion unit among many motion units. In view of the fact that the position of the equipment in the case is static during operation, the detection information transmission mode based on LAN is adopted, which avoids the lack of information distortion caused by various interference factors in the transmission process and realizes the real-time detection function of the manufacturing equipment movement unit parameters. Through research on fractional calculus theory, the application of the fractional differential operator has the dual advantages of enhancing signal strength and improving the accuracy of detection data, and realizes the high-precision acquisition function of detection information data. Therefore, the entire experimental process meets the three prerequisites for the evaluation of the manufacturing equipment operation status, and through the application of the method described in this paper, the research goal of real-time and accurate evaluation of the manufacturing equipment operation status is achieved.

Figure 4. Distribution of D-point detection data before and after 0.5 order differential fusion.

Figure 4. Distribution of D-point detection data before and after 0.5 order differential fusion.

7. Conclusions

Through an analysis of the operation characteristics of the manufacturing equipment, it is concluded that the operation status of the manufacturing equipment depends on the operation characteristics of the weakest motion unit in the equipment. It is proposed that the application of real-time acquisition of accurate values of important motion parameters in the motion unit can realize the concept of real-time evaluation of the operational status of manufacturing equipment by comparing it with its evaluation criteria. Through the research of fractional calculus theory, this paper proposes to combine the advantages of IoT technology and fractional differential algorithm, use IoT technology to realize the real-time collection function of detection information data, apply a fractional differential operator to significantly reduce the difference between detection data, and solve the technical problem of real-time acquisition of high-precision detection values of the main motion parameters of the manufacturing equipment movement unit, which lays a foundation for the implementation of real-time and accurate evaluation schemes for manufacturing equipment operation status. Through application in a real-time evaluation experiment of the machining center running state, the effectiveness and feasibility of the method described in this paper in engineering practice are verified.