Abstract

Gears are the most important parts of rotating machinery and power transmission devices. When gears are engaged in meshing transmission, vibration will occur due to factors such as gear machining errors, meshing rigidity, and meshing impact. The traditional FxLMS algorithm, as a common active vibration algorithm, has been widely studied and applied in gear transmission system active vibration control in recent years. However, it is difficult to achieve good performance in convergence speed and convergence precision at the same time. This paper proposes a variable-step-size multichannel FxLMS algorithm based on the sampling function, which accelerates the convergence speed in the initial stage of iteration, improves the convergence accuracy in the steady-state adaptive stage, and makes the modified algorithm more robust. Simulations verify the effectiveness of the algorithm. An experimental platform for active vibration control of the secondary gear transmission system is built. A piezoelectric actuator is installed on an additional gear shaft to form an active structure and equipped with a signal acquisition system and a control system; the proposed variable-step-size multichannel FxLMS algorithm is experimentally verified. The experimental results show that the proposed multichannel variable-step-size FxLMS algorithm has more accurate convergence accuracy than the traditional FxLMS algorithm, and the convergence accuracy can be increased up to 123%.

1. Introduction

As a general-purpose part, gear is an important transmission device in various mechanical equipment systems and plays a very key role in equipment in machinery, transportation, chemical, aviation, aerospace, shipbuilding, and other industries. In actual applications, the gears inevitably have manufacturing and installation errors, and the gear teeth will deform under the action of the load. These errors and deformations destroy the meshing relationship of the gear transmission and make the position of the gear me\shing deviate from its theoretical position so that the instantaneous transmission ratio changes, causing the collision and impact between the teeth and forming the gear meshing. Internal excitation results in vibration and noise [14].

The vibration of the gear system will not only produce noise and lead to the instability of the transmission system, but also accelerate the fatigue damage of the transmission system and cause its failure and serious consequences [5]. For example, one of the main excitation sources of the helicopter body is the vibration from the transmission system [6]. Excessive vibration levels will cause discomfort to the crew, cause fatigue, affect work efficiency, and reduce the precision and reliability of precision instruments and weapon aiming systems, or even fail to work normally. Besides, vibration and noise will cause fatigue damage and damage to some components, which directly affects the flight reliability of the helicopter [79]. The problem of vibration and noise of a helicopter is more prominent than that of other aircraft. If part of the vibration and noise can be eliminated during the operation of the helicopter, it has important practical significance to ensure the flight safety of the helicopter and improve the ride comfort [10, 11].

For the study of the vibration of the gear transmission system, in 1967, British scholar H. Optiz published an article on the dynamic performance of spur gear and the helical gear. According to the mechanism of gear vibration and noise, it is concluded that the vibration of gearbox and noise is a function of gear transmission error and precision, and some valuable analysis curves are given. Since the mid-1970s, the research objects of many scholars are mainly flexible rotors under high speed, such as rotors on high-speed high-precision machine tools and centrifugal units, which are usually under supercritical speed. How to restrain the rotor vibration, prevent the overall instability, and ensure the safe operation of the rotor has become a key problem. In the aspect of active control actuator, electromagnetic bearing and electromagnetic damper are studied more. As early as 1980, Haberman and Brunet [12] studied the stiffness characteristics of electromagnetic bearing and proposed the “bathtub” shaped complex stiffness curve. Since then, most researchers have analyzed the supporting characteristics of the electromagnetic bearing system based on this. Zhang and Zhang [13] reviewed in detail the research progress of unbalanced control technology of electromagnetic bearing. Maslen et al. [14] made a comprehensive review and summary of the existing unbalanced vibration control methods and the latest research results in the field of electromagnetic bearing. Huang and Lin [15] used a fuzzy control method to conduct theoretical and experimental research on unbalanced vibration control of conical electromagnetic bearing. Jiang and Zmood [16] designed a controller based on H∞ for self-excitation and external interference. Electromagnetic dampers have the advantages of noncontact, no mechanical friction and wear, no lubrication, long service life, and controllable stiffness and damping. Kasarda et al. [17] applied the electromagnetic damper to the active control of subsynchronous resonance of high-speed single-disk rotor system, which reduced the resonance amplitude by 93%. Nonami et al. applied the electromagnetic damper to the vibration control of the flexible rotor shaft. The control signal of the damper was provided by the analog controller, which effectively reduced the resonance peak of the flexible rotor shaft. Elbuken et al. [18] applied the electromagnetic damper to the vibration control of the rotating suspension and effectively reduced the vibration of the rotating suspension by exerting control force on the rotating suspension through the electromagnetic damper. Using the idea of finite element model and considering the influence of meshing error, Guan et al. [19, 20] established the coupled vibration model of the gear transmission system with bending and torsion and deduced the transmission function of the vibration from the gear excitation to the gearbox body, which laid a foundation for the vibration control of the rear gear transmission system. According to the criteria of modal control force and energy consumption, Guan evaluated four different ways of actuator installation in gearbox vibration active control system. The results show that the actuator is connected by extra bearings on the shaft in the best way and is more convenient to install. In the research of gear vibration, many scholars [2123] also adopt this actuator installation method.

Montague et al. [24] as early as the active control method used in gear vibration control in 1994, they adopt the method of feedforward control needed for the inhibition of gear vibration control was deduced, by manual adjustment of phase shift amplifier in the experiment to adjust the voltage signal of driving voltage actuator, experimental results show that the meshing frequency is 4500 Hz, with 75% of vibration attenuation. Subsequently, Sutton et al. [7] clamped the three magnetostrictive actuators on the support and controlled the vibration of the gear of the helicopter transmission in the frequency range of 250–1250 Hz to minimize the vibration kinetic energy of the supporting rod. Rebbechi et al. [25] adopt the feedforward controller to adjust the phase and amplitude of the control force to realize the vibration control of the gearbox box body. The experiment shows that the vibration at the meshing fundamental frequency, double frequency, and triple frequency of the gearbox and the radiation noise at the meshing fundamental frequency is greatly reduced. Chen and Brennan’s [26] simple gear mesh model has set up the control needed to control the size of the torsional vibration being deduced, based on the nonlinear system design of the adaptive controller to control the harmonic frequency. The experiment showed that stiffness modulation factor is proportional to the load on the gear,and the torsional vibration damping is about 7 dB at meshing frequencies of 150–350 Hz. Guan et al. [27] proposed a direct hybrid adaptive control method based on the Lyapunov stability theorem to simultaneously adjust feedback gain and feedforward gain to control gear vibration, and the results showed that this method was not sensitive to gear meshing frequency error. Li et al. [28, 29] installed the actuator on the gear shaft with additional support bearings, combined with the phase-locked loop for frequency estimation of reference signals, and used the filtered LMS algorithm to control the gear vibration. After that, many scholars have improved the filtering LMS algorithm and applied it to active vibration control and achieved good control results [3039]. Wu and Lin [40] used a digital processor (DSP) to combine three different control algorithms to form a hybrid controller to control the vibration of the gearbox, but it was too complex to be applied in practice. Li et al. [29], based on frequency estimation technology, presented an improved filtering LMS algorithm that is used to control the vibration of the gearbox. To avoid out-of-band overadjustment, an integrated adaptive linear enhancer is also used.

Most of the above researches on vibration control of gearbox adopt a single gear transmission system, so single-input and single-output control methods are adopted. In this paper, a two-stage gearbox is taken as the controlled object, and there is two gear meshing excitation in the system. To achieve global vibration suppression of the system, multiple sensors are needed to capture the dynamic vibration of the structure and multiple actuators to suppress the energy generated by excitation [4143]. Therefore, this paper will adopt a multichannel control method [44, 45] to control the vibration of the two-stage gear transmission system actively. As a common model, the traditional FxLMS algorithm has been widely studied and applied in vibration control in recent years, but it is difficult to achieve good performance in convergence speed and convergence accuracy. To solve the problems existing in traditional FxLMS algorithms, a multichannel FxLMS algorithm with variable step size is proposed in this paper. The algorithm accelerates the convergence speed in the initial stage of convergence, improves the convergence accuracy in the steady-state adaptive stage, and makes the system have certain robustness [46].

Traditional vibration control methods such as modal control, pole assignment, and optimal control are based on the precise model of controlled structure, but the gear transmission system is a vibration system with serious uncertainty and difficult to accurately model. It is difficult to design a steady feedback system to achieve the desired control performance by using traditional control methods. However, the adaptive control method is less dependent on prior knowledge and does not depend on specific models, so it is an ideal choice to apply the adaptive control method to the active control of gear vibration.

2. Establishment of the Two-Stage Transmission Gearbox Model

UG software can get different specifications of parts by modifying the parameters. It only needs to input some basic parameters of gears, such as tooth number, modulus, pressure angle, and tooth surface width. UG software can automatically generate the required gear model. This method can also be used for gear shaft and bearing, which greatly improves design efficiency. The main parameters of the gears in the gearbox are calculated by the formula derivation [47]. The basic data of gears are shown in Table 1.

The shaft is designed according to the size of the gear. The main radial size of the shaft is shown in Figure 1. Among them, (a) is the high-speed shaft and low-speed shaft; (b) is the intermediate shaft. The size of the box is determined according to the size of the shaft and the size of the gear. The internal size of the box is 282 mm × 135 mm × 270 mm. Use UG to generate a three-dimensional model according to the designed size, and the established three-dimensional model is shown in Figure 2. The motor is an adjustable speed motor, which can adjust the speed of the gearbox in a large range. The load is a magnetic powder brake; on the one hand, it can adjust the load, allowing a more comprehensive analysis of the dynamic characteristics of gear load. On the other hand, it could make gear mesh closer by applying a load. Due to gear wear and insufficient precision, there will be gaps between the gear components. Under low or light load conditions, the gear meshing may swing back and forth, which will have a certain impact on the acquired signal. This effect can be reduced by applying a load.

The details of the active control structure inside the gearbox are shown in Figure 3. The active control structure is composed of an actuator bracket, a piezoelectric actuator, an antiload bar, and a support bearing. The bottom of the bracket is bolted to the wall of the box, and one end of the actuator is clamped with the top sleeve. The other end of the actuator is connected with the antiload bar, and the antiload bar is fixed on the support bearing. The antiload bar is a slender shaft with a certain degree of flexibility for overload protection. Besides, to improve the control accuracy, the direction of the control force output of the piezoelectric ceramic actuator is parallel to the direction of the meshing line of the gear.

The 3D model of the two-stage gearbox was imported into ADAMS, and the virtual prototype model was constructed by adding constraints. The functional expression of contact force iswhere K is the stiffness coefficient, q0 is the reference collision distance of two contact objects, q is the actual collision distance of two contact objects, e is the index in the rigidity, C is the damping rate, and d is the distance that the damping rate reaches the maximum.

The gear is made of no. 45 steel, and the bearing is made of type 6003 deep groove ball bearing. For the contact collision between the gears, the contact collision parameters between the bearing ball and the inner and outer rings are set as shown in Table 2.

Add the constraints of each component, set the speed on the input shaft, and apply load on the output shaft, and the virtual prototype of the two-stage gear transmission system can be obtained, as shown in Figure 4.

To verify the usability of the virtual prototype and provide a theoretical basis for vibration active control joint simulation, dynamic simulation in ADAMS, first set the multistage gearbox virtual prototype input speed to 2540 rpm and the load to 1 N·m. At this time, the meshing fundamental frequency of low-speed shaft gear is 500 Hz, and that of high-speed shaft gear is 804 Hz. The simulation time is set as 0.5 s, and the simulation step length is 0.0001 for simulation verification. Figures 5(a)5(d) show the vibration acceleration of the low-speed axis and the high-speed axis in the time domain and frequency domain, respectively.

3. Description of the Improved FxLMS Algorithm

3.1. Traditional FxLMS Algorithm

The online identification method of secondary channels is an additional random noise method first proposed by Eriksson and Allie [48], as shown in Figure 6. P(n) is the transfer function from the primary source to the control point. S(n) is a transfer function of the secondary path from the actuator to the control point and is a very complex electro-acoustic coupling system, which is mainly composed of D/A conversion, loudspeaker, secondary acoustic channel, and A/D conversion. X(n) is the reference signal. E(n) refers to the residual vibration response signal detected by the error sensor at the control point.

Random noise (n) is used as the input of the secondary channel adaptive online identification filter S'(n). When the identification filter converges, S'(n) can uniquely converge to S(n). In the adaptive control process, S'(n) can be used as an estimated transfer function of the secondary channel S(n). However, this method will add random noise to the input of the actuator y(n) and cause additional vibration of the control point through the secondary channel, and the vibration cannot be eliminated in the whole identification process of the secondary channel.

The calculation process of the algorithm is as follows:

According to the LMS algorithm, the weight iteration of online identification of secondary channels can be expressed aswhere is the step size of filter S’(n), andwhere M is the length of filter S’(n).

At this time, the weight update of the FxLMS algorithm is expressed by the following formula:where μ is the step size of the FxLMS algorithm and N is the length of the weight W(n) of the FxLMS algorithm.

The traditional FxLMS algorithm updates iteratively with a fixed step size. When the step size factor is larger, the convergence rate is faster in the early stage, but in the steady-state adaptive stage, the step size is too large, and the oscillation occurs, resulting in a large error. On the contrary, when the step size factor is small, the steady-state error is small, but the convergence rate in the early stage will decrease greatly. Moreover, the traditional FxLMS algorithm mainly studies the single-channel problem. With the increasing degree of mechanical complexity, it is difficult to meet the control requirements of a single channel.

3.2. Improved Multichannel Variable-Step-Size FxLMS Algorithm

To solve the problem of the traditional fixed step size algorithm, many kinds of research have been made on the variable-step-size LMS algorithm. In this paper, three classical variable-step-size LMS are compared. Based on the variable-step-size LMS algorithm based on the sampling function, a certain change is made to the form of the variable-step-size LMS algorithm. In the case of little impact on the convergence rate, the steady-state error is smaller, and it has a better anti-interference ability.

The mechanism of variable-step-size LMS algorithm is a function model between step size factor and error. Three typical function models are as follows.

Step factor μ(n) is the Sigmoid function of error e(n), as shown in the following equation:

Step factor μ(n) is the tongue function of error e(n), as shown in the following equation:

Step factor μ(n) is the sampling function of error e(n), as shown in the following equation:where and take the same parameters as 4 and 8, respectively, to obtain the nonlinear curve of step factor and error signal of the above algorithm.

From Figure 7, we can see that compared with the sigmoid function and tongue function, the LMS algorithm based on the sampling function has a larger step size in the early stage, so its convergence rate is relatively fast. In the steady-state adaptive stage, the steady-state error is small because of the smaller dμ/de(n) value. On this basis, the function is further improved.where the values of and are still 4 and 8; the curves are compared in Figure 8.

It can be seen from Figure 8 that the improved LMS algorithm based on sampling function does not reduce the step size significantly in the early stage but significantly reduces the step size in the steady-state adaptive stage, which will lead to a decrease in the convergence rate, but significantly improves the convergence accuracy. Moreover, the autocorrelation value of the error signal is used to adjust the step factor in the functional relation between the step factor and the error so that the step factor has the anti-interference ability.

The iterative formula of the improved LMS adaptive algorithm based on the variable step size of the sampling function is as follows:

From (10), we can get

According to the principle of adaptive filtering, the above equation can be rewritten aswhere is the optimal weight coefficient of the filter, and the mean value of is 0, which is unrelated to the input signal.

Substitute (12) into equation (10) to getwhere is the weight coefficient deviation.

So, we get

According to the fact that (n) is a Gaussian white noise signal, its correlation is poor, that is,  = 0, and is not related to the input signal X(n), the following can be obtained:

From equation (15), we can see that in the function of step factor and error signal obtained by the improved LMS adaptive algorithm based on variable step size of sampling function, the influence of noise signal is eliminated, so the algorithm is not affected by noise signal and has a strong anti-interference ability.

The improved LMS algorithm is applied to the multichannel active vibration control system. The improved FxLMS algorithm is the multichannel variable step size FxLMS algorithm. The schematic diagram of the active control structure is shown in Figure 9, where x(n) represents the reference signal, S represents the secondary channel between all actuators and sensors, W represents the channel variable step FxLMS algorithm controller in the multichannel vibration control system, and V represents the white noise output to the secondary channel. Suppose the system has J active controllers, K sensors, the reference signal is written as x(n), and the FxLMS controller W is written aswhere is the weight vector of the jth controller, which is to write the ownership value coefficient in the FxLMS algorithm as a weight vector.

The multichannel active vibration control system has J × K secondary channels; secondary channel S can be written aswhere represents the secondary channel between the jth active actuator and the kth sensor. represents the filter coefficient identified using the secondary channel between the jth active actuator and the kth sensor.

4. The Cosimulation

Combined with the multichannel FxLMS algorithm structure in Figure 10 and the active control structure of the two-stage gearbox in Figure 2, the active vibration control system of the gearbox is designed as shown in Figure 10 below. The vibration caused by excitation and the gear pair is the control target. and , respectively, represent two secondary channels. To eliminate the influence of secondary channels, corresponding models and are introduced into the control system to filter the input signals of the controller. In adaptive active vibration control, the selection of reference signals has a great influence on the control effect. To obtain the reference signals with high correlation with the vibration of the gear system, the frequency estimation of the vibration signals collected in real time is carried out through the second-order IIR adaptive notch filter [49] to obtain the meshing frequency of the gear pair and then synthesize the reference signals online.

Through simulation, the effectiveness of the improved FxLMS algorithm in the application of two-gearbox vibration active control is verified. In Simulink to build vibration active control system, its core is the improved FxLMS adaptive filter, due to the multiple-input multiple-output filter, the Level-2 function improved FxLMS algorithm, the adaptive trap of reference signal synthesizer using the same Level-2 function to carry on the design, use of FIR filter structure simulation of the secondary channel, connect the controller, the reference signal synthesizer, and the secondary channel model, and get the vibration active control system. The virtual prototype model of the gear transmission system in ADAMS was imported into Simulink as a mechanical module, which was connected with the control system to obtain the active vibration control cosimulation model of gear transmission integrated with the mechanical system and the control system, as shown in Figure 11.

The input speed of the gearbox is 2540 rpm. The meshing frequency of the high-speed shaft gear pair is 800 Hz and the meshing frequency of the low-speed shaft gear pair is 500 Hz. The load is 1 N·m. Set the simulation parameters as follows: the adaptive filter order as the controller is set to 80, the step size factor of the algorithm adopts the variable-step-length factor based on the sampling function, and the sampling frequency is selected as 10 kHz. Joint simulation is carried out, and the simulation results are shown in Figure 12.

Figures 12(a) and 12(c) show that compared with the traditional FxLMS algorithm, the improved FxLMS algorithm has a faster convergence rate and a smaller steady-state error. In Figure 12(a), the FxLMS control system reaches stability after about 0.15 s, while the improved FxLMS control system reaches stability after about 0.07 s. In Figure 12(c), the FxLMS control system reaches stability after about 0.12 s, while the improved FxLMS control system reaches stability after about 0.08 s. Moreover, it can be seen from Figures 12(a) and 12(c) that the improved FxLMS control system has a smaller vibration acceleration after stabilization is achieved.

Figures 12(b) and 12(d) show that the suppression of vibration by the active vibration control system is mainly at the 2nd and 3rd harmonics. In Figure 12(b), the vibration reduction of the FxLMS control system at the 2nd harmonic 1000 Hz and the 3rd harmonic 1500 Hz is 10 dB and 14 dB, respectively. The vibration reduction of the improved FxLMS control system at 1000 Hz and 1500 Hz is 23 dB and 31 dB, respectively. In Figure 12(d), the vibration reduction of the FxLMS control system at 1608 Hz and 2412 Hz of the second- and third-order harmonics is 17 dB and 28 dB, respectively. The vibration reduction of the improved FxLMS control system at 1608 Hz and 2412 Hz is 24 dB and 37 dB, respectively. It is further proved that compared with the traditional FxLMS algorithm, the improved FxLMS algorithm has a smaller steady-state error.

5. Design and Verification of Experiments

The experimental platform of the active vibration control system of the gearbox system can be divided into three parts, namely, the gear transmission system, the vibration measurement processing system, and the active control system. The established platform is shown in Figure 13.

The gear transmission system consists of a DC motor, motor governor, coupling, gear transmission part, and magnetic powder brake. The gear transmission part is composed of two pairs of meshing gears, the modulus is 2, the pressure angle is 20°, the number of teeth of high-speed shaft meshing gear is 19 and 37, respectively, and the number of teeth of low-speed shaft meshing gear is 23 and 35. Adjustable speed motor speed range is of 0–1800 rpm. The two couplings are, respectively, connected with the motor and the high-speed shaft, the magnetic powder brake, and the low-speed shaft. Magnetic powder brake can provide 0–15 N·m torque.

The vibration measurement processing systemis used to detect the vibration signal of the gearbox and the error signalinput to the controller. The acceleration sensor is used to collect the signalof the casing, and the collected signal is displayed in real time with theanalysis software DASP-V11 and is input to the controller as a calculatederror signal to calculate and output the control signal.

The active control system includes a controller, a piezoelectric actuator, a piezoelectric actuator power supply, and an actuator clamp. For the control system, this paper selects the Speedgoat real-time simulation system as the control system of active vibration control and establishes a semiphysical simulation control platform. The Speedgoat real-time simulation system provides a MATLAB/Simulink driver module that converts the algorithm control system established on Simulink into the required C code and transmits the generated object code to the processor, which calculates the voltage signal of the piezoelectric actuator and controls the piezoelectric actuator.

To verify the effectiveness of the proposed multichannel active vibration control structure for the vibration control of the gear transmission system, the proposed active control algorithm was programmed, compiled, and downloaded to the real-time target machine as the controller of the active control system on the MATLAB/Simulink platform, and the multichannel variable-step-size FxLMS algorithm was adopted in this paper. To ensure the effectiveness of the control system, the sampling time of the controller is set to 0.0001 s. The motor speed was adjusted to 1200 rpm, the load was adjusted to 1 NM, and the active vibration control was started. The vibration signals of the gearbox body were collected through the acceleration sensor, as shown in Figure 14. The vibration signals before and after the control were Fourier transformed to obtain the vibration frequency domain control diagram, which was compared with the traditional FxLMS algorithm, as shown in Figure 15. The meshing fundamental frequencies of the low-speed shaft and high-speed shaft are 236 Hz and 380 Hz, respectively.

It can be seen from the time-domain diagram of vibration control in Figure 14 that the vibration attenuation begins at 8 s when the control begins, and the convergence is completed at the 12th second, and the vibration attenuation reaches 70% before the control. As can be seen from the local enlargement of Figure 15, the control effect is obvious at the basic frequency and multiple frequencies of gear meshing. The specific meshing frequency data is shown in Table 3 below. It can be seen from Table 3 that the decay of the meshing frequency of the low-speed shaft is about 3-4 dB, and that of the high-speed shaft is 3–5.5 dB. According to the data in Table 3, it can be calculated that the convergence accuracy has increased by 23%–82%.

When the speed is adjusted to 1530 rpm, the meshing fundamental frequency of the low-speed shaft and high-speed shaft is 301 Hz and 485 Hz, respectively. The vibration generated by the gear transmission system is controlled, and the control frequency domain curve is shown in Figure 16. It can be seen from the figure that the attenuation at low-speed shaft meshing frequency is 4.5-5 dB, and the vibration attenuation at high-speed shaft meshing frequency is about 8-9 dB. The specific meshing frequency data is shown in Table 4. According to the data in Table 4, it can be calculated that the convergence accuracy has increased by 36%–123%.

Based on the active vibration control results of the gear under two rotating speeds, we can know that the proposed multichannel variable-step-size FxLMS algorithm has an effective suppression effect on the vibration at the high harmonic in the gear transmission system, and the vibration attenuation range is 3–9 dB. It can be seen that the attenuation of the meshing frequency of the high-speed shaft is higher than that of the low-speed shaft. Compared with the traditional FxLMS algorithm, the proposed multichannel variable-step-size FxLMS algorithm has a better control effect. The control accuracy can be increased up to 123%.

6. Conclusion

In this paper, a novel active vibration control platform based on a piezoelectric actuator is constructed, which is composed of hardware and software such as a two-stage transmission device, control system, signal acquisition, and processing system. The multichannel variable-step-size FxLMS algorithm is used as the controller to control the vibration of the gearbox. Experimental results show that the multichannel variable-step-size FxLMS algorithm is effective for vibration control of the gear system, and its control effect is better than the traditional FxLMS control, which can reduce the vibration by 3–9 dB at the target. The control accuracy is improved by 23%–123%.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by China National Natural Science Funds (51405169) and the Fundamental Research Funds for the Central Universities (no. ZQN-704).