1. Introduction
Electromagnetic actuation finds application in many engineering sectors since its driving magnetic field can easily and safely propagate through many materials and reach confined spaces. These unique features make this actuation form well-suited for multiple applications broadly related to robotics, such as electromagnetic actuation for image stabilization [
1], minimally invasive microrobots for microsurgery [
2,
3,
4,
5,
6,
7,
8], targeted drug delivery inside the human body [
9,
10] or early screening and cancer treatment [
11], magnetic levitation (
e.g., vertical motion of a ball [
12]) and general motion control of ferromagnetic elements (
e.g., planar steering of a ferrofluid drop [
13]).
It is well-known that electromagnets can only pull ferromagnetic objects without the ability to push them away. Then, precisely manipulating a ferromagnetic object at a distance from the electromagnet(s) is challenging for multiple reasons. Modeling the electromagnetic force is quite uncertain because it is a nonlinear function of the actuated object position and electromagnet coil current (the magnetic flux and the resulting electromagnet’s pulling force drop quickly with the distance from its surface). Analytical expressions available to predict the electromagnetic force have limited accuracy and are not readily usable for control design due to their complexity, even when referencing idealized scenarios [
14,
15,
16,
17]. Moreover, the presence of a ferromagnetic element in the close vicinity of an electromagnet can significantly distort the magnetic flux’s distribution, an effect not normally incorporated in analytical prediction methods (refer to
Section 2.3 for more details). Fully understanding the system behavior is, therefore, hard in this context.
An additional difficulty takes over when the goal becomes achieving closed-loop position control of the actuated element. The object position must be sensed without interfering with its motion. Thus, feasible alternatives for industrial settings use laser sensors [
18] or optical tracking [
13]. Other expensive techniques exist (
e.g., magnetic resonance imaging, computed tomography, ultrasound, or infrared radiation [
2]). Still, they are more suitable for biomedical applications, such as the magnetic particle imaging used for electromagnetically navigating drug delivery [
19]. Laser sensing is costly and has some challenges, such as safety-related issues because lasers can harm the human eye. Optical tracking can be cost-effective, or at least affordable, depending on the selected camera used for visual feedback (the required frame rate, which is application-driven, plays a significant role). In any case, the latter technique is promising when cost-effectiveness is sought.
Concerning motion control in industrial applications, the steering of ferromagnetic particles (ferrofluid) was studied in 2D. The proposed form of planar manipulation relied on an array of four electromagnets actuated by feedback control based on visual feedback [
20]. The ferrofluid could follow arbitrary trajectories with reasonable accuracy, which was improved using optimal control techniques [
13,
21]. Despite the remarkable tracking abilities, this approach needs complex control design and is limited to systems with slow dynamics (the feedback loop runs at 15 Hz). Then, an approach with 3D motion capabilities in high-viscosity fluids drove a 1.46 mg magnetic cylinder via eight electromagnets [
22]. The magnetic force distribution in the workspace was not precisely known, so the controller involved a PID design whose tuning was not model-based. A physical-based control design is, in fact, very computationally demanding. A black-box approximation based on neural networks was proposed, but the system's active control was not further investigated [
23]. Moreover, the literature on the control of levitating (solid) objects is extensive. Typical techniques are complex due to the strong nonlinearities of magnetic levitating systems. They involve neural networks [
18], feedback linearization [
12,
24], model predictive control [
25], or adaptive sliding mode control [
26]. These controllers require the measurement of the actuated object’s position with a fast enough sampling time. Precise knowledge of the electromagnetic force is also crucial for control design, a rarely well-defined scenario.
Referring to motion control in the medical field, tiny dimensions are commonplace. The first example concerns a microrobot with a predefined magnetization profile [
27]. The time-varying magnetic field imposed with open-loop control around the robot generates different locomotion modes (
e.g., walking, crawling, or jumping). Then, a cylindrical magnetic millirobot can navigate in soft tissues because it is attracted to the center of a permanent magnet array [
28]; position control is achieved by x-ray imaging and a motion stage that slides and rotates the soft tissues. Three permanent magnets (only one of them is displaced by a stepper motor) enable precise control of another microrobot to perform complex trajectories and deploy cargo [
29]. This method leverages open-loop control, unlike the closed-loop technique developed for a pivot walking robot [
30]. The latter device employs visual feedback but needs a complex magnetic actuator with three nested electromagnetic Helmholtz coils.
In summary, closed-loop position control with electromagnetic actuation requires the accurate and contactless measurement of the object position, the accurate prediction of the electromagnetic force, and the implementation of a suitable control strategy to close the loop. Despite the extensive investigations carried out to date, one aspect that has not been adequately addressed is the characterization of the electromagnetic force, which most of the time has been predicted analytically or numerically with limited accuracy, restricting the actuation to small objects whose presence does not excessively distort the magnetic field. In this paper, we propose a novel and cost-effective electromagnetic set-up to perform position control via visual feedback and a model-based control law. Specifically, the set-up comprises one webcam for position tracking and one electromagnet actuating a 10 mm diameter steel ball hanging on a low-stiffness spring. Notably, the electromagnetic force has been directly measured, thereby eliminating the uncertainties associated with analytical and numerical prediction methods. The present set-up demonstrates good tracking performance with an actuated object big enough to significantly distort the electromagnetic field, even though all the components used to assemble the set-up are off-the-shelf and cost-effective.
The rest of this paper is organized as follows. After introducing our novel set-up (
Section 2.1), we validate the proposed package for optical tracking of the moving object (
Section 2.2). Then, we experimentally identify the field distribution produced by the electromagnet used in our study and the electromagnetic force it generates on the actuated object (
Section 2.3). We also characterize the friction force of the moving object in different scenarios (
Section 2.4) and propose a feedback control law to perform closed-loop position control (
Section 2.5). Finally, we report our experimental results collected under different operating conditions (
Section 3) and conclude our discussion with critical remarks and suggestions for future developments (
Section 4).
Author Contributions
Conceptualization, D.P., A.C., and A.A.; methodology, S.C., D.P., A.C., and A.A.; validation, S.C.; formal analysis, D.P. and A.C.; investigation, S.C.; data curation, S.C., D.P. and A.C.; writing—original draft preparation, D.P. and A.C.; writing—review and editing, S.C., D.P., A.C., and A.A.; visualization, S.C., D.P., and A.C.; supervision, D.P., A.C., and A.A.; funding acquisition, D.P. All authors have read and agreed to the published version of the manuscript.
Figure 1.
The novel set-up for electromagnetic actuation proposed in this paper: (a) system overview, (b) zoomed-in view of the actuator and moving object, and (c) simplified view of the moving object with its free-body diagram. .
Figure 1.
The novel set-up for electromagnetic actuation proposed in this paper: (a) system overview, (b) zoomed-in view of the actuator and moving object, and (c) simplified view of the moving object with its free-body diagram. .
Figure 2.
Block diagram showing the structure of the experimental set-up.
Figure 2.
Block diagram showing the structure of the experimental set-up.
Figure 3.
Representative results for motion tracking of the oscillating pendulum: (a) displacement time-series, (b) detailed view of displacement time-series, (c) displacement normalized histogram, (d) displacement power spectrum (PSD).
Figure 3.
Representative results for motion tracking of the oscillating pendulum: (a) displacement time-series, (b) detailed view of displacement time-series, (c) displacement normalized histogram, (d) displacement power spectrum (PSD).
Figure 4.
Electromagnet experimental characterization at 24 V supply voltage: (a) axial magnetic flux, (b) radial magnetic flux, (c) axial magnetic flux measured in the presence of the ball, (d) radial magnetic flux measured in the presence of the ball. The black circle in panels (c) and (d) represents the ball, while the coordinate system is placed on the electromagnet surface with its origin on the geometrical center of the core.
Figure 4.
Electromagnet experimental characterization at 24 V supply voltage: (a) axial magnetic flux, (b) radial magnetic flux, (c) axial magnetic flux measured in the presence of the ball, (d) radial magnetic flux measured in the presence of the ball. The black circle in panels (c) and (d) represents the ball, while the coordinate system is placed on the electromagnet surface with its origin on the geometrical center of the core.
Figure 5.
(a) Magnetic force measurements at various distances between the electromagnet surface and the ball center for different electromagnet supply voltages (4, 8, 12, 16, 20, and 24 V): the discrete data points are the measured values, while the continuous lines are the predictions of the fitting relation in Equation (7); (b) parity plot with comparison between magnetic force predictions from Equation (7) and measurements (the dashed lines are ±20% error bounds).
Figure 5.
(a) Magnetic force measurements at various distances between the electromagnet surface and the ball center for different electromagnet supply voltages (4, 8, 12, 16, 20, and 24 V): the discrete data points are the measured values, while the continuous lines are the predictions of the fitting relation in Equation (7); (b) parity plot with comparison between magnetic force predictions from Equation (7) and measurements (the dashed lines are ±20% error bounds).
Figure 6.
Measured and simulated trends of the damped free oscillations in (a) air, (b) water, and (c) glycol. The ball is in contact with the electromagnet surface at the beginning of the tests and is released with a step command that removes the applied voltage at 0 s. The plots use a different time scale for better visualization.
Figure 6.
Measured and simulated trends of the damped free oscillations in (a) air, (b) water, and (c) glycol. The ball is in contact with the electromagnet surface at the beginning of the tests and is released with a step command that removes the applied voltage at 0 s. The plots use a different time scale for better visualization.
Figure 7.
Block diagram of the proposed closed-loop position control algorithm used to generate the voltage applied to the electromagnet.
Figure 7.
Block diagram of the proposed closed-loop position control algorithm used to generate the voltage applied to the electromagnet.
Figure 8.
Closed-loop motion control of the ball surrounded by air with a constant-frequency position command (
Video S1 and control parameters
1,
2.2 V/mm,
1.8 V/mm/s, and
7 Hz): (a) ball position, (b) position error, and (c) voltage command.
Figure 8.
Closed-loop motion control of the ball surrounded by air with a constant-frequency position command (
Video S1 and control parameters
1,
2.2 V/mm,
1.8 V/mm/s, and
7 Hz): (a) ball position, (b) position error, and (c) voltage command.
Figure 9.
Closed-loop control of the ball in the air with a variable-frequency command (
Video S2 and control parameters
1,
0.05/30 1/s,
1.5 V/mm,
1.05 V/mm/s, and
15 Hz): (a) ball position, (b) position error, (c) voltage command, and (d) input frequency.
Figure 9.
Closed-loop control of the ball in the air with a variable-frequency command (
Video S2 and control parameters
1,
0.05/30 1/s,
1.5 V/mm,
1.05 V/mm/s, and
15 Hz): (a) ball position, (b) position error, (c) voltage command, and (d) input frequency.
Figure 10.
Closed-loop motion control of the ball surrounded by water with a constant-frequency command (
Video S3 and control parameters
1,
1.8 V/mm,
1.2 V/mm/s, and
25 Hz): (a) ball position, (b) position error, and (c) voltage command.
Figure 10.
Closed-loop motion control of the ball surrounded by water with a constant-frequency command (
Video S3 and control parameters
1,
1.8 V/mm,
1.2 V/mm/s, and
25 Hz): (a) ball position, (b) position error, and (c) voltage command.
Figure 11.
Closed-loop control of the ball in the water with a variable-frequency position command (
Video S4 and control parameters
1,
0.005 1/s,
1.05 V/mm,
0.7 V/mm/s, and
30 Hz): (a) ball position, (b) position error, (c) and voltage command.
Figure 11.
Closed-loop control of the ball in the water with a variable-frequency position command (
Video S4 and control parameters
1,
0.005 1/s,
1.05 V/mm,
0.7 V/mm/s, and
30 Hz): (a) ball position, (b) position error, (c) and voltage command.
Figure 12.
Closed-loop control of the ball in the glycol with a constant-frequency position command (
Video S5 and control parameters
1,
3 V/mm,
2.5 V/mm/s, and
25 Hz): (a) ball position, (b) position error, and (c) voltage command.
Figure 12.
Closed-loop control of the ball in the glycol with a constant-frequency position command (
Video S5 and control parameters
1,
3 V/mm,
2.5 V/mm/s, and
25 Hz): (a) ball position, (b) position error, and (c) voltage command.
Figure 13.
Closed-loop motion control of the ball in the glycol with a variable-frequency position command (
Video S6 and control parameters
1,
0.1/30 1/s,
0.8 V/mm,
0.7 V/mm/s, and
30 Hz): (a) ball position, (b) position error, (c) voltage command.
Figure 13.
Closed-loop motion control of the ball in the glycol with a variable-frequency position command (
Video S6 and control parameters
1,
0.1/30 1/s,
0.8 V/mm,
0.7 V/mm/s, and
30 Hz): (a) ball position, (b) position error, (c) voltage command.
Figure 14.
Comparison of the actual position error resulting from the tests reported above and performed in different conditions (air, water, and glycol): (a) constant-frequency position command, and (b) variable-frequency position command.
Figure 14.
Comparison of the actual position error resulting from the tests reported above and performed in different conditions (air, water, and glycol): (a) constant-frequency position command, and (b) variable-frequency position command.
Table 1.
Synthesis of the control settings used to perform our experimental validations.
Table 1.
Synthesis of the control settings used to perform our experimental validations.
Test |
(-) |
(-) |
(V/mm) |
(V/mm/s) |
(Hz) |
Air (constant-frequency) |
1 |
0 |
2.2 |
1.8 |
7 |
Air (variable-frequency) |
1 |
0.05/30 |
1.5 |
1.05 |
15 |
Water (constant-frequency) |
1 |
0 |
1.8 |
1.2 |
25 |
Water (variable-frequency) |
1 |
0.005 |
1.05 |
0.7 |
30 |
Glycol (constant-frequency) |
1 |
0 |
3 |
2.5 |
25 |
Glycol (variable-frequency) |
1 |
0.1/30 |
0.8 |
0.7 |
30 |
Table 2.
RMS values in mm of the actual position error resulting from the tests reported above and performed in different conditions (air, water, and glycol).
Table 2.
RMS values in mm of the actual position error resulting from the tests reported above and performed in different conditions (air, water, and glycol).
Test |
RMS value |
|
Test |
RMS value |
Air (constant-frequency) |
0.1884 |
|
Air (variable-frequency) |
0.3857 |
Water (constant-frequency) |
0.1612 |
|
Water (variable-frequency) |
0.3471 |
Glycol (constant-frequency) |
0.1491 |
|
Glycol (variable-frequency) |
0.2896 |