1. Introduction
Vibration machines (VMs) and technologies play a crucial role across various industries, including manufacturing and agriculture. It is challenging to envision conducting tasks such as material separation, vibratory transportation, grinding, rolling, mixing, compaction, pile driving, and more without the aid of VMs and associated devices [
1]. Their utilization not only yields significant technical and economic benefits but also contributes to enhancing working conditions.
Among the primary sources for inducing oscillations in vibration machines, unbalanced actuators stand out as the most commonly employed [
2,
3,
4]. These actuators find application in various technologies such as vibratory conveyors, lifters, screens, and rammers, where circular, elliptical, or directional vibrations are necessary. Circular vibrations are typically produced by an unbalanced rotor, or two unbalanced rotors, revolving in phase in the same direction, or anti-phase—in the opposite directions. When there exists a phase discrepancy between the driving forces of two-shaft vibration actuators, elliptical oscillations manifest themselves.
In the VM applications for tasks like screening, crushing, and bulk material transportation, enhancing productivity often hinges on maintaining a stable synchronous rotation mode of unbalanced rotors (vibration actuators). The foundations of VM synchronization theory were established in the 1960s by I.I. Blekhman [
3,
5,
6]. This theory laid the groundwork for synchronization mechanisms in traditional rigid or flexible transmissions, gradually showing the advantages of vibration equipment powered by two or more actuators. Such systems enable the vibration mashine to produce various trajectories through a synchronization phenomenon, see [
4]. The multiple synchronization modes are based on the self-synchronization phenomenon, wherein vibration actuators rotate at average angular velocities that are multiples of each other. This mode leads to more efficient vibrating behavior, especially in challenging tasks such as transporting dusty, sticky, or wet granular media.
Studying the flow field variation in twin-roll strip casting processes induced by oscillating rollers was carried out in [
7]. This research is focused on achieving a homogeneous distribution of strip impurities, reducing grain size, and thereby enhancing strip quality through the adoption of the proposed vibrating casting.
Electro-hydraulic vibration equipment finds extensive application in various vibration environment simulation tests. Liu et al. [
8] addressed the problems related to bandwidth, waveform distortion control, stability, offset control, and complex waveform generation in high-frequency vibration conditions.
Shi et al. [
9] concentrated on the vibration synchronization of mechanical structures propelled by two vibroactuators. They introduced a synchronization mechanism combining space and plane patterns through the extension of vibrator installation to arbitrary directions in spatial axes. Altshul et al. [
10] and Panovko et al. [
11] examined variations in the arrangement of technological loads on the machine’s working body. It was demonstrated how shifts in the load center from the machine’s constructive axis of symmetry lead to changes in resonant frequencies and mutual phasing of unbalances, and the feasibility of controlling resonant VMs by adjusting vibro-actuator power frequencies to compensate for uncontrolled load shifts was discussed.
In the domain of dual-mass vibration systems, Liu et al. [
12] investigated vibratory synchronization phenomena driven by two actuators. They employed averaging methods to deduce synchronization and stability criteria in synchronous states, quantitatively analyzing parameters such as actuator rotational speeds, phase differences, and responses in sub-resonant and super-resonant states, accompanied by an engineering example. Using asymptotic and average methods, Zhang et al. [
13] derived theoretical conditions for implementing multiple synchronization and stability, with a particular focus on analyzing synchronization for four actuators. Leniowska and Sierże [
14] studied the controlled vibration damping of a round plate utilizing a controller with phase shift adjustment in the feedback loop. A novel controller was proposed, combining regulator structures with positive position feedback and strain rate feedback. Experimental studies involve measuring plate vibrations using a laser vibrometer, with a control signal applied to the plate via a Macro Fiber Composite (MFC) disk affixed to its center. Their findings demonstrate the feasibility of employing this solution to mitigate plate vibration effectively.
The phenomenon of transporting solid and granular bodies over oscillating rough surfaces has found numerous technical applications. Vibrational movement arises from an average directed motion of material particles relative to horizontally oscillating uniformly rough surfaces, stemming from the asymmetry of surface oscillation shapes, as expressed in the inequality of time intervals between consecutive extrema of surface acceleration oscillations, see [
15,
16,
17]. Theoretical investigations by Blekhman et al. [
17] focused on the oscillatory transportation of solid particles over flat surfaces under non-translational vibrations. This work, with experimental validation conducted on a vibration platform, revealed that the existing theory, with additional parameters introduced, can effectively set up the speed of vibrational transportation. Advancements towards intelligent technological systems involve employing mechatronics principles, particularly computer-controlled vibration setups with feedback mechanisms [
17].
The introduction of feedback control in manipulating the vibration field has been addressed by Fradkov et al. [
18], who were initially focused on simple synchronization cases. However, it became clear that the appearance of multiple synchronous modes introduces asymmetry, enhancing vibration movement efficacy by creating complex trajectories for rotors and carrier bodies. Nevertheless, the stability of multiple-synchronous movements remains a challenge, necessitating further exploration into controlling vibration using feedback control for vibroactuators.
This paper studies the vibrational system controlled by the method proposed by Fradkov et al. [
18] from the perspective of creating the various vibration fields for the operational platform, driven by the actuators with unbalanced rotors. Experimental validation conducted at the Multi-resonance Mechatronic Laboratory Setup confirms the approach’s efficacy, demonstrating its practical applicability in controlling the type of vibration fields of the working platform, crucial for bulk material transportation processes.
The subsequent sections are structured as follows:
Section 2 briefly outlines the Mechatronic Setup SV-2M utilized for experiments.
Section 3 presents the bidirectional control law for multiple synchronization of the unbalanced rotors. The experimental results along with the data processing algorithm are presented in
Section 4. Concluding remarks and future work intentions are outlined in
Section 5.
3. Bidirectional Control Law for Multiple Synchronization of Unbalanced Rotors
In [
18], the problem of controlled synchronization for unbalanced rotors’ phase shift was considered. The control law presented in [
18] for the multiple synchronization of unbalanced rotors has the following form:
where
,
are the motors’ reference velocities;
is the normalization coefficient;
,
are the errors of motor velocities;
,
,
are the integrals or tracking errors
,
,
(respectively), which are used in the PI-controllers; PI-controllers for the velocities of the left and right motors are described by (5), (6), respectively, where
,
are for the integral, and
,
are for the proportional controller gains; in (7),
and
are the phase angles of the rotors, and
denotes the normalized phase shift between the rotors; the phase shift error is denoted by
, where
is the reference phase shift angle between the rotors;
and
in (10) denote the control signals applied to the left and right drive systems produced by the corresponding PI-controllers.
In [
18], the linearized system stability analysis for model (
2) and control law (
3)–(
10) was performed. The results were confirmed by the computer simulations and the numerous real-world experiments on the setup SV-2M.
4. Experimental Results
In the present paper, the experimental study results of [
18] are extended to studying the corresponding vibration platform deflections and summarized in the form of the spectral densities of selected points of platform movements and the “vibration fields” of the platform.
The main aim of this study is to demonstrate the vibration fields by the example of moving the platform points #1–#4 (see
Figure 3).
4.1. Data Processing Algorithm
The following processing algorithm is implemented to obtain the type of vibration fields based on experimental data.
Platform position sensors data are acquired with a sampling interval and saved to a file for further offline post-processing. Given the range of platform oscillation frequencies, an interval of s was chosen, corresponding to a Nyquist frequency rad/s. At the end of the experiment, the data array was saved to the hard disk drive (HDD) for subsequent processing.
Data acquisition from the sensors is performed using Simulink Desktop Real-Time toolbox, which includes Simulink I/O driver blocks enabling closed-loop control of physical systems from a desktop computer, allowing connections to sensors, actuators, and other devices. Data Type Conversion blocks are used to convert the measurements to MATLAB double type data for further processing. Simulink blocks “To Workspace” are used to save the data to the workspace. After the experiment, using the MATLAB save routine, the data obtained are saved on the HDD in the form of a MATLAB mat-file.
Since the measurement characteristics of optical sensors
s1,
s6, given by
,
, are close to the inversely proportional mappings presented in
Figure 2, the displacements
, at times
,
of the corresponding platform points
s1s6 are restored based on the sensor readings. Before the experiments, optical sensors position had been calibrated to obtain their real measurement characteristics and the results were tabulated in 20 points inside the working range (see
Figure 2). Sensor measurements sometimes include random spikes, which are erroneous readings outside the normal range. Therefore, these anomalies were removed before further processing. For this purpose, each sensor’s entire array of measurements was subjected to statistical analysis using partitions the values of
into
M bins, and returning the count in each bin.
bins was set for data processing. Next, values that occur in the measurement array less than the specified threshold value
are selected (
was set) and are replaced by the arithmetic mean of the measurement array. To recover displacement values
from sensor signals
, an spline-interpolation procedure was used on the every data processing step
.
For the obtained data arrays, spectrograms and the filtered centered platform displacements are calculated and displayed to the user by the following sequence of operations.
- 1.
Find time-average values of the displacement processes of the sensor positions.
- 2.
Apply Discrete Fourier Transform (DFT) to the resulting centered processes .
- 3.
Calculate the squared modulus, normalized by the number of points.
- 4.
Convert the resulting complex-valued array of images into the array of spectral densities .
- 5.
Select frequencies , , where bursts were observed.
- 6.
Clear sensor measurements using the DFT concerning .
- 7.
Find the set of frequencies in the vicinities of such as: , where is a chosen bandwidth relative to the peak frequencies .
- 8.
Form new complex-valued arrays by zeroing the values of except those that correspond to the frequencies from , which values are copied from , and their spectral densities were found.
- 9.
Find as the complex-valued inverse Fourier transforms to .
- 10.
Find the real part of .
- 11.
Apply the Inverse Discrete Fourier Transform (IDFT) to .
- 12.
Calculate the filtered centered platform displacements as , where stands for the Euclidean norm of the complex-valued vector argument.
As an example of the described procedure application, consider the case of filtering data of
,
, obtained for
rad/s,
rad/s,
, taking
rad/s,
rad/s, and
rad/s. The resulting time histories for
,
along with
are depicted in
Figure 4. The graphs obtained for different values of
show the influence of this parameter on the result of measurement filtering. The choice of frequency filter parameters is always a compromise between noise alleviation and possible signal distortion in the area of essential frequencies. The decision is assigned to the researcher, making it based on the study of processes for specific conditions.
After
are found, the platform points of interest #1–#4 (see
Figure 3) displacements
,
,
on the vertical plane
are calculated as follows:
The experimental results for various revolving speeds of rotors and phase shifts are presented below. For the legibility of the graphs, the trajectories are displayed only in the final section of the experiment with a duration of , which is about two characteristic periods of oscillation.
4.2. Case of Identical Rotation Velocities
At first, consider the case when both rotors revolve with the same desired frequency. Namely, let
rad/s be taken (i.e., in (
3)–(
10) parameter
be set). The desired phase shift between the rotors
is set to zero. The spectral densities for centered displacements
,
of sensors #2, #3, are pictured in
Figure 5. As is seen from the plot, the “outstanding” peak on frequency close to
rad/s appears. It corresponds to the low-frequency resonant mode of the platform. The calculated value of the resonance frequency is quite close to the experimental one, which is confirmed by the Formula (
1) for given in [
19] values
N/m,
kg,
kg,
rad/s. The appearance of a small oscillation peak in the region of the frequency 1 rad/s can be treated as a manifestation of oscillations arising in the process of the phase shift control.
The motion of the various points of the platform on the vertical plane
(where
stands for the horizon axis,
denotes the vertical one), are plotted in
Figure 6. The time interval for plotting the curves is taken as
s. This interval corresponds to six periods of motion with a frequency of 80 rad/s.
The upper plots series correspond zero phase shift (
) between the rotors. Respectively, the lower plot is related to phase shift
. It is seen that the vibration field, having, in general, a similar shape to that in the case of
, changes its inclination. This property can be used for control of the vibrational transportation speed and direction. Of course, for practical application, it is essential what kind of movement the processed stuff acquires. This dependence can be quite complex and determined not only by the platform movement but also by the physical characteristics of the substance. Some preliminary analysis is contained in the works [
15,
16,
17], but it is not exhaustive and does not cover the whole variety of vibration fields that can be provided by the method proposed here. The studies conducted show its potential and a necessity to be continued for real machines and situations.
4.3. Case of Different Rotation Velocities
Additionally, similar experiments were performed for the case of different rotation frequencies of rotors. Namely,
rad/s and
rad/s were taken (i.e., in (
3)–(
10) parameter
was set). The desired phase shift between the rotors
is set to zero. The corresponding spectral densities for centered displacements
,
of sensors #2, #3 are pictured in
Figure 7. As is seen from the plot, the spectrograms have three peaks: on the frequencies of 45 rad/s, 90 rad/s, corresponding the rotation velocities of the drives, and the peak on 20 rad/s appears related to the natural oscillating of the platform.
The corresponding vibration fields of the platform sections #1–#4 on the vertical plane
are plotted in
Figure 8. The time interval for plotting the curves is taken as
s. This interval corresponds to six periods of motion with a frequency 90 rad/s. The platform motion is for the considered multiple synchronization case, which is more complex compared with the single-frequency one due to interference both from the slow oscillations with the natural frequency of the platform (20 rad/s) and two rotation speeds of the drives with different values (90 and 45 rad/s).
The upper plots in
Figure 8 correspond to zero phase shift (
) between the rotors, while the lower plot is related to phase shift
. It is seen that the vibration field is more complex than in the single-frequency case and also varies with the phase shift, set between the rotors. This gives additional flexibility in control of the platform vibration fields and can be used in developing smart vibrational technologies.
As for the coincident frequencies considered in the previous subsection, and for different frequencies, the results for which are presented here, it is clear that by changing the phase shift it is possible to significantly tune the vibration fields. The desired type of field depends on specific conditions and needs to be studied from a technological point of view. However, according to the obtained graphs, it is clear that the conclusion made based on simulations in [
19], that in the case of the double synchronous mode the ordinate of the lower point of the trajectory is lower than in the case of the single synchronous mode, is confirmed experimentally, which improves the efficiency of unloading and prevents congestion formation (see plots
in
Figure 8, right column).
5. Conclusions
The paper demonstrates the possibility of using controlled multiple synchronization of unbalanced rotors driven by induction motors to form vibration fields in the working platform of a vibration mechatronic setup. The experimental results obtained at the MMLS SV-2M of the IPME RAS are presented, demonstrating the efficacy of the proposed approach. The method proposed here obviously leads to simpler algorithms than in [
19] which simplifies its implementation but makes the system rougher, without the possibility of fine-tuning.
Expanding the capabilities of controlled synchronization can be achieved by employing unidirectional (master–slave synchronization). Since the control method using a PI controller may be insufficient, especially in theoretical research, the alternative control methods including adaptive control, variable-structure control on sliding-mode, speed gradient techniques, Neural Networks, nonlinear correction methods, Field-Oriented Control (FOC) or Linear Quadratic Regulator (LQR), could be employed in future investigations, representing promising areas for further research.