Next Article in Journal
Optimizing Biodegradable Starch-Based Composite Films Formulation for Wound-Dressing Applications
Next Article in Special Issue
Magnetically Driven Manipulation of Nonmagnetic Liquid Marbles: Billiards with Liquid Marbles
Previous Article in Journal
Precise Control of Glioma Cell Apoptosis Induced by Micro-Plasma-Activated Water (μ-PAW)
Previous Article in Special Issue
Nonsingular Integral Terminal Sliding Mode Control for Resonant Frequency Tracking of Electromagnetic Acoustic Transducers (EMATs) Based on Fixed-Time Strategy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

2D Magnetic Manipulation of a Micro-Robot in Glycerin Using Six Pairs of Magnetic Coils

1
College of Internet of Things Engineering, Jiangnan University, Wuxi 214000, China
2
Shenzhen International Graduate School, Tsinghua University, Shenzhen 518055, China
3
Jiangsu Key Laboratory of Advanced Food Manufacturing Equipment and Technology, Wuxi 214000, China
*
Author to whom correspondence should be addressed.
Micromachines 2022, 13(12), 2144; https://0-doi-org.brum.beds.ac.uk/10.3390/mi13122144
Submission received: 22 October 2022 / Revised: 20 November 2022 / Accepted: 30 November 2022 / Published: 4 December 2022
(This article belongs to the Special Issue Micro- and Nano-Systems for Manipulation, Actuation and Sensing)

Abstract

:
This paper demonstrates the control system of a single magnetic micro-robot driven by combined coils. The combined coils consist of three pairs of Helmholtz coils and three pairs of Maxwell coils. The rotating magnetic field, gradient magnetic field, and combined magnetic field model of the combined coils were analyzed. To make the output magnetic field quickly converge to the reference point without steady-state error, the discrete-time optimal controller was designed based on the auto disturbance rejection technology. We have designed a closed-loop controller based on a position servo. The control system includes the position control and direction control of the micro-robot. To address problems with slow sampling frequency in visual feedback and inability to feed real-time position back to the control system, a Kalman filter algorithm was used to predict the position of the micro-robot in two-dimensional space. Simulations and experiments were carried out based on the proposed structure of combined coils and control scheme. The experimental results demonstrated the uniformity and excellent dynamic performance of the generated magnetic field.

1. Introduction

In recent years, the broad attention on micro-robots has resulted in their rapid growth. Because of their small volume and flexibility in narrow and complex environments [1,2], unrestricted and controllable micro-robots are widely used in various fields. In industry, they can be used for micro-manipulation, transportation, classification, and micro-assembly [3] of micro-objects. The high liquidity and repeatability of micro-robots make them effective tools for interacting with biological cells, which can be combined with cell operations. For example, micro-robots can be operated at the cellular or subcellular level, allowing effective in vitro interactions for moving and classifying cells [4,5,6]. The integration of microfluidic chips based on micro-electromechanical system (MEMS) with micro-robot technology is an innovation in biomedicine [7,8,9]. High-precision motion control and high-propulsion pumping machines enable micro-robots to carry multiple objects and transport them to a desired location in microfluidic chips [10]. Applications of micro-robots in vivo are broad and include minimally invasive surgery, targeted drug delivery, brachytherapy, hyperthermia, the removal of obstacles by mechanical means, and serving as fundamental structures in the body. Micro-robots can be manipulated remotely in the human body at the micro/nanoscale, including the circulatory system, urinary system, and central nervous system, facilitating applications of capsule endoscopy [11,12], ophthalmic surgery [13], and cardiovascular interventional surgery [14,15]. Therefore, the development of micro-robots will undoubtedly contribute to multi-field developments.
Micro-robots are small and cannot be powered by a built-in power supply, so they can only be driven externally. A number of driving methods have been developed, such as the use of dielectrophoretic force generated by electric fields to drive micro-scale dielectric particles [16], piezoelectric driving [17], thermal driving [18], electro-osmotic force driven [19] by biological bacteria [20,21], and micro-motors [22] driven by chemical fuel. All of these driving methods have certain limitations, especially when applied to biomedicine. Specifically, piezoelectric driving requires high voltage, the thermal driving method can easily result in excess temperatures, and the biological bacterial driving method requires monitoring to maintain low cytotoxicity. External magnetic fields [23] are also viable power sources for driving micro-robots. Magnetic fields have proven to be harmless to organisms [23,24] and have the ability to penetrate deep tissues. Therefore, the magnetic force and moment generated by magnetic fields can be applied to micro-robots without being affected by biological fluids [25].
There are usually two methods to realize magnetic field driving of micro-robots. The first method is generating the required magnetic field with a permanent magnet. For instance, the in vivo capsule endoscope can be driven with an in vitro handheld permanent magnet [11]. Nguyen and others developed the Stereotaxis Niobe system (Niobe®, Stereotaxis, Inc., St. Louis, MO, USA), which adopts the magnetic field driving method to drive a permanent magnet for catheter operation [26]. However, the permanent magnet driving method normally requires an additional manipulator to control the movement and rotation of the permanent magnet; compared with the second driving method based on electromagnetic coils, this method cannot generate any real-time controlled magnetic field. Therefore, many studies have selected magnetic systems driven by electromagnetic coils [27,28,29,30,31,32,33,34,35]. In [29], an enhanced electromagnetic driving device with a large workspace was developed, which could realize the trajectory tracking of magnetic particles with a proportion integration differentiation (PID) controller in two-dimensional and three-dimensional space. In [36], a position estimator was designed according to the dynamic model of the micro-robot, and trajectory tracking in the two-dimensional space of the micro-robot was realized with a self-made electromagnetic operation system.
However, the two electromagnetic driving devices mentioned above cannot produce constant magnetic flux density and magnetic field gradient. They can only be estimated approximately when analyzing the magnetic force model such that the precise dynamic model of the micro-robot cannot be obtained, which poses a higher challenge to the precise motion control of the micro-robot. As Helmholtz coils can produce a magnetic field with constant flux density, when energizing a sinusoidal current, they can produce a rotating magnetic field, whereas Maxwell coils can produce a constant magnetic field gradient. In context, micro-robots can therefore be rotated by the rotating magnetic field of a Helmholtz coil, and driven by the magnetic field gradient generated by a Maxwell coil. In [37], an intravascular electromagnetic drive system was proposed that consisted of three pairs of Helmholtz coils and one pair of Maxwell coils. The rotating magnetic field generated by the Helmholtz coils was used for the steering and rotating motion of the micro-robot, and the gradient magnetic field generated by the Maxwell coils was used for the forward motion of the micro-robot. The feasibility of the scheme was subsequently verified via experiment. In our study, three pairs of Helmholtz coils and three pairs of Maxwell coils were used to generate uniform magnetic field and uniform gradient magnetic field in space, respectively.
Some researchers have studied the motion control methods of micro-robots. In one such study, Helmholtz electromagnetic coils were manufactured to guide the micro-robot to a desired position [38]. Based on analysis of the dynamic model of the micro-robot, a controller based on neuro-fuzzy network was designed. Finally, the performance of the control system was verified by experiments. An optimal path planner based on particle swarm optimization algorithm and a robust model predictive trajectory tracking controller were designed. Fluorescent imaging was used for visual feedback, and the automatic control of the magnetic micro-robot was realized [39]. To prove the controllability and observability of the system, a micro-robot motion control method combining high-gain observer and adaptive backstepping control was proposed in [40]. While simulations and experiments verify the effectiveness and robustness of this method, some important issues have been ignored in these studies. To achieve high-precision driving in the micro-robot, it is necessary to accurately control the input current of the coil [41]. However, the inductance of the coils causes a large time delay in the system, such that when it is necessary to generate a high-frequency magnetic field, the tracking performance of the micro-robot will become unsatisfactory. As the working time increases, the coil temperature rises, and the linear relationship between the input current and the output magnetic field is destroyed, causing some interference in the system. Another problem worth studying is that industrial cameras are widely used to directly obtain the position information of micro-robots. Due to the influence of the sampling frequency in the visual feedback, the real-time position of a micro-robot cannot be immediately fed back to the control system, so the precise path tracking of the micro-robot cannot be realized.
In this paper, we propose a full-state feedback controller based on auto disturbance rejection technology to solve the asynchronous problem between the magnetic field generated by the combined coils and the input current. A circuit board based on a power inverter was used to control the voltage applied to the coil, and simulation analysis was performed in MATLAB software. At the same time, the step response, sinusoidal response, and impulse response were simulated for separate tracking. Compared with traditional PID control and open-loop control, the designed controller achieved two important functionalities: (1) the output magnetic field and magnetic field gradient could converge to the reference point; (2) optimal transient response was achieved. The control system includes a direction control module and a position control module. The direction control module was configured as a closed loop of the direction variable, and the position control module was configured as three closed loops of position, velocity, and acceleration variables, respectively. As the sampling frequency in visual feedback was slow and the real-time position could not be fed back to the control system immediately, a Kalman filter algorithm was used to predict the position of the micro-robot.
The rest of the paper is organized as follows. In Section 2, we introduce the structure of the electromagnetic coil and the dynamic model of the electromagnetic coil and the micro-robot. In Section 3, we describe the designs of the driving circuit and auto disturbance rejection controller as well as the driving model of the electromagnetic coils. In Section 4, we propose the dual closed-loop control, with which a Kalman filter algorithm is used for position prediction. Section 5 presents the simulation results of the Maxwell and Helmholtz coils, illustrates the motion control experiments for driving the micro-robot along a square and S-shape and plots for the trajectory tracking error, and provides an evaluation of the performance of the Kalman filter. The experimental results demonstrate that the proposed Kalman filter algorithm can predict the position of the control scheme effectively and improve the dynamic performance of the system. Finally, the conclusion of this work is given in Section 6.

2. Electromagnetic Coil Structure and Model

2.1. Electromagnetic Coil Structure

The coil structure that was used in the magnetic driving system is shown in Figure 1. The combined coils included three pairs of Helmholtz coils and three pairs of Maxwell coils, which were mutually perpendicular along the x-axis, y-axis, and z-axis. The distance between the Helmholtz coils was equal to the radius of the coils. Currents of equal value and direction that were injected into the two coils could produce constant magnetic flux density in the central area. If the currents of the two coils were in opposite directions, a rotating magnetic field was generated that could align the micro-robot in a desired direction, as shown in Figure 2a. The distance between the Maxwell coils was 3 times of the radius, and the two coils were supplied with currents of equal value and opposite directions to generate a constant magnetic field gradient, which was employed to drive the micro-robot, as shown in Figure 2b. Therefore, a uniform magnetic torque and magnetic force could be generated in the central area of the combined coils by combining Helmholtz coils and Maxwell coils [42,43].

2.2. Electromagnetic Coil Model

The micro-robot could be controlled by magnetic torque and magnetic force in the magnetic field. Torque and magnetic force are closely related to the magnetic flux density and magnetic field gradient. The force and the magnetic torque can be expressed as follows [44,45]:
F = V M B
T = V M × B
where F is the magnetic force, T is the magnetic torque, V is the volume of the micro-robot, B is the magnetic flux density (T), and M is the magnetization of the micro-robot (A/m). The torque of the micro-robot is determined by the magnetic field generated by the Helmholtz coil, and the magnetic force is determined by the gradient magnetic field generated by the Maxwell coil. Therefore, it is necessary to analyze these two kinds of magnetic fields to control the steering and position of the micro-robot.
In addition, the Helmholtz coil magnetic field can produce a rotating magnetic field. The normal vector of the rotating magnetic field is defined as n = (nx, ny, nz)T, and the magnetic field can be expressed as
B H = B H X B H Y B H Z = B 0 cos ( 2 π f t ) u n y n x 0 + B 0 sin ( 2 π f t ) v n x n z n y n z n x 2 n y 2
where BHX, BHY, and BHZ are respectively expressed as the magnetic flux density generated by the x-axis, y-axis, and z-axis Helmholtz coils. f is the frequency of the rotating magnetic field, B0 is the modulus of the magnetic flux density, u = n y 2 + n x 2 , v = ( n x n z ) 2 + ( n y n z ) 2 + ( n x 2 + n y 2 ) 2 [46,47]. The magnetic flux density B0 produced by the uniaxial Helmholtz coil in a uniform region is proportional to the input current [48]:
B 0 = ( 4 5 ) 3 / 2 μ 0 N a 2 I H
where N and a respectively represent the number of turns per coil and the radius of the Helmholtz coils, μ 0 is the vacuum permeability, and IH is the current energized by the Helmholtz coils.
When nx = ny = 0, the micro-robot will rotate in two-dimensional space, meeting the research requirements of this paper. In this scenario, only the x-axis and y-axis Helmholtz coils are generating magnetic fields:
B H = B H X B H Y B H Z = B 0 cos ( 2 π f t ) B 0 sin ( 2 π f t ) 0
The magnetic flux density generated by the x-axis Maxwell coil in the region of the central uniform gradient magnetic field is [49]
B M X = k I M X 0.5 I M Y 0.5 I M Z x y z
where k = 16 3 ( 3 7 ) 5 / 2 μ 0 N a 2 , N, and a represent the number of turns per coil and the radius of the coil, μ 0 is the vacuum permeability, and IMX, IMY, and IMZ are respectively expressed as the currents injected by the x-axis, y-axis, and z-axis Maxwell coils.
Considering the magnetic field in other directions, the synthetic magnetic field of the Maxwell coil is
B M = B M X + B M Y + B M Z = 1 0.5 0.5 0.5 1 0.5 0.5 0.5 1 k M X I M X k M Y I M Y k M Z I M Z x y z

2.3. Dynamic Model of the Micro-Robot

It can be seen from Equation (5) that the micro-robot can be rotated in any direction in two-dimensional space by injecting dynamic current into the Helmholtz coil. However, to realize the movement of the micro-robot in the plane, the gradient magnetic field drive is required. The Maxwell coils can provide magnetic force to drive the micro-robot, and when the magnetization direction of the micro-robot is consistent with the magnetic field gradient, it has the maximum magnetic force. Therefore, in our study, the uniform magnetic field was used to change the direction of the magnetic torque of the micro-robot and the direction of the desired magnetic force.
The micro-robot, when in glycerol solution, will be affected by magnetic force, gravity, viscous resistance, friction and other forces. For a micro-robot moving in two-dimensional space, the force in the vertical direction does not need to be considered. Since a cylindrical micro-robot was used, the contact area between its surface and the bottom of the container was very small. The friction in this case is close to zero but cannot be ignored, so it was regarded as a disturbance in the system. Therefore, the dynamic equation for the magnetic micro-robot that we used is
F m + F d + F f = m   p ¨ .
Here, Fm is the magnetic force applied to the micro-robot by the external magnetic field generated by the Maxwell coil. Fd represents the resistance of the micro-robot moving in glycerol solution. Ff is the friction between the micro-robot and the bottom of the container. The moving medium of the micro-robot in the experiment was glycerol, which is a solution of low Reynolds number. The viscous resistance of a micro-robot placed in a low-Reynolds-number solution is [50]
F d = 6 π η R v ,
where η represents the dynamic viscosity of the glycerol solution, R is the radius of the micro-robot, and v represents the speed of the micro-robot.
It is assumed that the coil produces a magnetic field Bk(p) at point p such that for a current injected into the coil ik, the kth coil produces a magnetic field B ¯ k ( p ) at point p when a current of 1A is passed. Therefore, the magnetic field generated by the combined coils at point p can be expressed by the following formula:
B ( p ) = k = 1 6 B k ( p ) = k = 1 6 B ¯ k ( p ) i k = [ B 1 ( p ) B 6 ( p ) ] i 1 i 6 = β ( p ) I
where k = 1, 2,…, 6 respectively represent the x-axis, y-axis and z-axis of Helmholtz coils and the x-axis, y-axis, and z-axis of Maxwell coils. I denotes i 1 i 6 T . The gradient components of the generated magnetic field are
B ( p ) x = B ¯ 1 ( p ) x B ¯ 6 ( p ) x i 1 i 6 = β x ( p ) I B ( p ) y = B ¯ 1 ( p ) y B ¯ 6 ( p ) y i 1 i 6 = β y ( p ) I B ( p ) z = B ¯ 1 ( p ) z B ¯ 6 ( p ) z i 1 i 6 = β z ( p ) I
As the direction of the internal magnetic torque of the magnetic micro-robot will be consistent with the direction of the external magnetic field, and as the magnetic field generated by the Maxwell coil in the central area can be safely ignored, according to (1) and (2), the magnetic torque and magnetic force of the micro-robot in the magnetic field can be calculated as
B F = β ( p ) m T β X ( p ) m T β Y ( p ) m T β Z ( p ) i 1 i 6 = Λ B , F ( m ,   p ) I
Given the desired path of the micro-robot, the required magnetic field Bdes and magnetic force Fdes can be determined according to the desired motion direction and position. The input current of each coil can be obtained by solving Equation (12):
I = Λ 1 B , F ( m , p ) B d e s F d e s

3. Optimal Control of Magnetic Field Drive

3.1. Design of Driving Circuit

The coil drive module was designed based on the power inverter. The control program is executed by the microprocessor, at which point the control input calculated in each loop is converted by the DSP into a pulse width modulation (PWM) signal through the power inverter (switching frequency is 20 kHz) on a full bridge circuit control board as shown in Figure 3.

3.2. Design of Auto Disturbance Rejection Controller

The design of this component is based on auto disturbance rejection control technology to solve the time delay problem caused by the self-induced coil currents and realize the rapid response of the coil-generated magnetic fields to the input current. Figure 4 shows the structure diagram of the auto disturbance rejection controller, which is composed of a tracking differentiator (TD), an extended state observer (ESO), and a state error feedback control law (NLSEF).
For current control, the discrete equation for the transition process of the tracking differentiator is expressed as
x 1 ( k + 1 ) = x 1 ( k ) + T x 2 ( k ) x 2 ( k + 1 ) = x 2 ( k ) + T f s t [ x 1 ( k ) x ( k ) , x 2 ( k ) , r , h ]
In Formula (14), x represents the input current as the input signal of the auto disturbance rejection controller, x1(k) is the tracking of the input signal x(k), x2(k) is the differential of x1(k), and T is the period of discrete control. r and h are adjustable parameters in the tracking differentiator TD. The larger that r is, the faster the tracking speed becomes, but this may exacerbate overshooting. h is the filter parameter: the noise filtering effect improves with the increase of h. Thus, we adjusted the parameters through the output waveform of the tracking differentiator. The optimal control function of the system is
f s t x 1 k x k , x 2 k , r , h = r a / d , a d , d = r h r sgn a , a > d , d = r h
a = x 2 k + a 0 d 2 sgn y k , y k > d 0 , d 0 = d h x 2 k + y k h , y k d 0 , d 0 = d h
y k = x 1 k x k + h x 2 k
a 0 = d 2 + 8 r | y ( k ) |
The second-order extended state observer is
e ( k ) = z 1 ( k ) y ( k ) z 1 ( k + 1 ) = z 1 ( k ) + T [ z 2 ( k ) β 01 f a l ( e ( k ) , α 1 , δ 1 ) + b u ( k ) ] z 2 ( k + 1 ) = z 2 ( k ) T β 02 f a l ( e ( k ) , α 1 , δ 1 )
In the above Formula (19), y(k) is the output current based on the active disturbance rejection controller. The estimation of output current is represented by z1(k), and z2(k) estimates the total disturbance of the system. α1, δ1, β01, and β02 are the four adjustable parameters of the expanded state observer; we usually set α1∈(0,1). To make the system disturbance observable by the extended state observer in time, the fal function is added to the extended state observer to achieve rapid convergence. Note that it is necessary to make α1 as small as possible to achieve faster tracking. δ1 is the filtering parameter of the ESO; in our case, usually set as δ1 = 0.01. The dynamic performance of the system is largely affected by β01 and β02. The estimation of the system state variables is mainly controlled by β01, while the estimation of the system disturbance is controlled by β02. Larger β01 and β02 are required to obtain a faster estimation convergence rate, but overlarge β01 and β02 will generate high-frequency noise and produce oscillation divergence. Usually, we set β01 = 1/T and β02 = 1/(5T2). The state variables of the system can be accurately observed by z1(k), and the disturbance of the system can also be accurately estimated by z2(k) after the parameters are set completely. The nonlinear combination function fal is
f a l ( e ( k ) , α 1 , δ 1 ) = | e ( k ) | α 1 sgn ( e ( k ) ) , | e ( k ) | > δ 1 e ( k ) δ 1 α 1 1 , | e ( k ) | δ 1
According to the output of the tracking differentiator and the ESO, the first-order nonlinear state feedback control law can be constructed as follows:
e 1 ( k ) = x 1 ( k ) z 1 ( k ) u 0 ( k ) = β 1 f a l ( e 1 ( k ) , α 01 , δ 01 ) u ( k ) = u 0 ( k ) z 2 ( k ) / b
where u(k) is the output control quantity of the auto disturbance rejection controller and α01, δ01, and β1 are the three adjustable parameters in the nonlinear feedback control law. Generally, α01∈(0,1) and δ01 = 0.01; their actual definition is similar to the parameters for α1 and δ1 in the ESO. The response speed of the system can be enhanced by increasing β1, but it will cause oscillation, so it is necessary to adjust β1 according to the overall control performance of the system.

4. Closed-Loop Control System

In our study, we generally adopted a dual closed-loop control structure. The specific control system diagram is shown in Figure 5.
In Figure 5, θdes and Cdes are the expected movement direction and position of the micro-robot on the basis of the planned path. The Kalman filter algorithm is used to predict the motion states of the micro-robot to address real-time delay in feeding back position due to the slow sampling frequency of visual feedback. The difference between the filtered real-time feedback motion direction θreal and the feedback real-time position Creal is calculated to acquire the momentary motion direction error θerr and position error Creal. The real-time output magnetic field and magnetic field gradient are applied to the micro-robot via direction and position control, respectively. The control system employs a dual closed-loop structure, which is composed of a position tracking loop with the real-time position and a direction control loop with the movement direction. The micro-robot position closed loop is composed of a position controller, a speed controller, and an acceleration controller. The output of the previous controller is the input of the latter. The final output current is injected into the Helmholtz coil and Maxwell coil through the designed ADRC to quickly generate the magnetic field for steering and movement, driving the micro-robot to follow the planned path. The dual closed-loop structure was adopted for the control system because, in the high-precision real-time control process of the micro-robot, the position error and the direction error cannot be completely eliminated at the same time in most cases. In addition, the position controller, speed controller, and acceleration controller all adopt traditional PID control. In our study, the working environment of the micro-robot was a glycerin solution with low Reynolds number; the inertia of the micro-robot when moving in it could basically be ignored, so there was no need to use the differential link to realize advance correction. Therefore, the direction controller adopted PI control.
Figure 6 is a schematic diagram of the path tracking of the micro-robot. Assuming that the position of the micro-robot is not on the planned path PnPn+1, it is located at C, but the movement direction is consistent with the expected path direction. Moving in the direction of the path, the direction error θerr will inevitably increase. Similarly, when the micro-robot is in the desired position, the direction of movement is not consistent with the desired direction, such that the position error Creal will inevitably increase. Consequently, use of a dual closed-loop structure can achieve high-precision path tracking.

Location Prediction

In the above-mentioned closed-loop control process, the position of the micro-robot is fed back using vision-based method during the movement. The actual position is compared with the expected position to obtain the next control signal. Therefore, the accuracy of the feedback position is directly related to the precise movement of the micro-robot. In actual visual feedback, since the frequency of the visual feedback is less than the control frequency, acqusition of position data at the control moment is unreliable between the two picture frames. We assumed a kth control signal in the two adjacent frames of i − 1 and i such that, in the actual movement process, the micro-robot will have passed the moment in the i − 1 picture to the kth control signal. At this moment, data that can be visually fed back are still the data in the i − 1th picture, such that the position information feedback of the micro-robot is not accurate. Therefore, predicting the position of the micro-robot at the time of the control signal between the i − 1 and i frame pictures is highly necessary to improve control accuracy.
To address this problem, we employed a Kalman filter to predict the position of the micro-robot between the two picture frames, which mainly include the prediction stage and the update stage. In the prediction stage, the momentary state is estimated by the state of the micro-robot at the previous moment, and then the state information of the micro-robot is forwarded to the update stage for updating the predicted state information to obtain more accurate position information. The Kalman filter does not predict all information at each control time; rather, the feedback information is used when the control time overlaps with the image feedback time.
Details of the adopted Kalman filter algorithm and the derived state transition matrix of the micro-robot are provided in the S1 section of the supplementary material.

5. Simulation and Experiments

5.1. Simulation

5.1.1. Magnetic Field of Combined Coils

To improve the accuracy of rotation and movement of the micro-robot, the uniformity of the magnetic field generated by the electromagnetic coils is very important. In this section, we discuss our finite-element analysis of the constant magnetic field, gradient magnetic field, and rotating magnetic field generated by the combined coil. The analysis was conducted using the multi-physics simulation software COMSOL according to the structure model of the combined coils.
The specific parameters of the combined coils are shown in Table 1.
In Table 1, HX represents the x-axis Helmholtz coil, and MX represents the x-axis Maxwell coils.
Figure 7 shows the simulation results of the magnetic field generated by the x-axis Helmholtz coil in the plane when a 4 A current was applied. Figure 7a shows the modulus field diagram of the magnetic flux density. To show the direction and intensity of the magnetic flux density, the vector field diagram of the magnetic flux density was simulated, as shown in Figure 7b. As shown in Figure 7c, a constant magnetic flux density of approximately 1 mT could be generated in the center of the workspace. Figure 8 shows the gradient distribution of the magnetic field generated by the x-axis Helmholtz coil on the central plane of space. Figure 8 shows that the magnetic field gradient is approximately 0. The simulation results demonstrate that Helmholtz coils can produce uniform magnetic flux density on the central plane. Therefore, Helmholtz coils could regulate the direction of the micro-robot by rotating it, but it could not apply propulsive force.
Given that Helmholtz coils can produce a uniform magnetic field, we injected a constant 4A current into three pairs of Helmholtz coils to analyze the resulting magnetic field. As shown in Figure 9, the superposed magnetic field generated by three pairs of coils in the central area is also uniform.
To verify the effectiveness of the rotating magnetic field, the z-axis is selected as the rotation axis, and the magnetic flux density is obtained by Equation (5):
B H = B H X B H Y B H Z = B 0 cos ( 2 π f t ) B 0 sin ( 2 π f t ) 0
From the simulation results of the uniaxial coil, the linear relationship between the magnetic flux density generated by each pair of Helmholtz coils and the input current can be approximated as follows:
B H X = 0.25 I H X B H Y = 0.25 I H Y B H Z = 0.25 I H Z
Combining (22) and (23) and setting B0 = 1 mT and f = 1 Hz, the input current of coil can be obtained by
I H X = 4 cos ( 2 π t ) I H Y = 4 sin ( 2 π t ) I H Z = 0
The simulation results of the rotating magnetic field rotating around the z-axis are shown in Figure S1. Since the period of the rotating magnetic field is 1 s, we recorded the magnetic flux density distribution in the x–y plane every 0.125 s.
Figure 10 and Figure 11 present the simulation results of the x-axis Maxwell coil, which shows that the gradient of the magnetic flux density was uniform and that a magnetic field gradient of 0.05 T/m was generated when 4A current was injected. Therefore, it was possible to obtain different magnetic forces for driving the micro-robot by changing the supply current in the Maxwell coil. Note, therefore, that reversal of the direction of the injected current between the two Maxwell coils will cause the direction of the magnetic field gradient to be reversed.
Figure 12 shows the magnetic flux density distributions generated by combining the x-axis Helmholtz coil and the x-axis Maxwell coil. The current for the Helmholtz coil and Maxwell coil was set at 4 A. The results show that a gradient magnetic field along the x direction could be generated.
To obtain the force in any direction on the plane of the micro-robot, the combined coils need to generate gradient magnetic fields in different directions. Figure 13 shows that the combined coils can generate gradient magnetic fields in different directions on the x–y plane.

5.1.2. Combined Coil Drive Based on ADRC

In Section 3, we described a combined-coil current-control strategy based on ADRC to improve the fast response performance of the magnetic field. This section describes the verification of the effectiveness of this method via simulation. A simulation platform of the electromagnetic coil system was built using the MATLAB/SIMULINK environment. Since any discrete time process consists of a step response, and as any continuous signal can be approximated by a finite number of discrete steps, a discrete system model with calibrated noise was used to obtain the step response, sinusoidal response, and impulse response, respectively, to evaluate the dynamic response performance of the system. The open loop control, PI control, and auto disturbance rejection control of the combined-coils control system were respectively simulated and then compared. The control frequency of the combined coils was set to 20 kHz.
The simulation results of step response, sinusoidal response, and impulse response of the Helmholtz coils and Maxwell coils are shown in Figure 14. The Helmholtz coils and Maxwell coils were used to generate a uniform magnetic field and gradient magnetic field, so the ordinates are the magnetic flux density value and the magnetic flux density gradient, respectively.
The magnetic flux density given in Figure 14a is 1 mT, and the step value was set as t = 0 s. Upon comparing the three control methods of the x-axis Helmholtz coil, under open-loop control, although the coil response had no overshoot, the time to reach stability was longer than PI control and the auto disturbance rejection control methods, 400 ms in duration. When the steady state was reached, the response was slower and there was a larger steady state error. Under PI control, the coil response performance improved, but maximum overshoot was 16% and there was a steady-state error. In addition, both the PI control and open-loop control exhibited oscillations, which are harmful to the system. Finally, under the auto disturbance rejection control method proposed in this paper, the response performance was further accelerated with only 20 ms required, and there was basically no overshoot, no steady-state error, and small oscillation. The setting magnetic flux density gradient of the Maxwell coil shown in Figure 14b is 0.05 T/m. The open-loop response performance was the worst, having the largest steady state error, and the steady state time reached up to about 500 ms. PI control improved the response time to 300 ms and the maximum overshoot to 18%. Thus, the response time and steady-state performance of the coil under auto-disturbance control were greatly improved. The performance of the sinusoidal response and impulse response were also improved.

5.1.3. Closed-Loop Control System

The results from our simulations show that the auto disturbance control strategy can make the magnetic field generated by the energized coil respond quickly enough to solve the time-delay problem of the combined-coil drive. This section discusses the simulation of the closed-loop control system described in Section 4 to verify the effectiveness of the proposed closed-loop control method and the effectiveness of the micro-robot path-tracking when the Kalman filter is implemented. In the simulation setting, the micro-robot tracks the step signal and the sinusoidal signal with and without the Kalman filter algorithm. The step signal and the sinusoidal signal were respectively set as
y 1 = 0.001 m y 2 = 0.001 sin 2 π t m
The response of the micro-robot tracking the step signal is shown in Figure 15a. According to the simulation results, the micro-robot basically has no overshoot or steady-state error under the Kalman filter algorithm while tracking the step signal. Compared with the result without Kalman filter algorithm, this approach shows better dynamic performance and steady-state performance. The response of the micro-robot tracking the sinusoidal signal is shown in Figure 15b. According to the simulation results, the micro-robot has obvious tracking error without the Kalman filter algorithm when tracking the sinusoidal signal.

5.2. Experiments

The micro-robot magnetic driving force generation and control system experimental platform was mainly composed of an electromagnetic coils module, a coils driving module, and a visual feedback module. The experimental environment is shown in Figure 16. The picture on the left is a partial view of the working space. The combined-coil parameters of the magnetic drive system in Figure 16 are the same as those in Table 1 for the simulation part. A cylindrical rubidium iron boron magnet with a diameter of 1 mm and a height of 1 mm was selected for the driving object. The micro-robot had a circular two-dimensional workspace with a diameter of 60 mm, and glycerin liquid with a kinematic viscosity of 350 cs was used for the moving medium of the micro-robot.

5.2.1. Combined Coils Drive Based on ADRC

In Section 5.1.2, the effectiveness of the combined-coils current-control strategy based on ADRC was verified through simulation. In this section, we describe our experiments for the further verification of the performance of this method. The dynamic performance of the coil was determined based on its movement trajectory and time during the actual micro-robot movement. In the experiments, the initial movement direction of the micro-robot s along the x-axis direction. When returning the micro-robot to the origin, application of force of equivalent magnitude along both the y-axis and x-axis should result in a y = x trajectory. Similar to the simulation in Section 3, open-loop control, PI control, and the proposed auto-disturbance control methods were used to realize the movement of the micro-robot, respectively. The experimental results are shown in Figure 17.
The motion trajectory of the micro-robot under three control methods is presented in Figure 17. At the beginning of the linear motion stage, the motion of the micro-robot under the auto-disturbance control is more stable, and there is basically no time delay for the change of response. Compared with open-loop control and PI control methods, the adjustment distance of the micro-robot under auto-disturbance control method was the shortest. There was a sudden change of adjustment distance under PI control due to overshoot, and the open-loop control method generated the longest adjustment distance. Based on whole movement process observation, the micro-robot movement under auto-disturbance control was the smoothest and basically had no vibration.

5.2.2. Closed-Loop Motion Control

To demonstrate the performance of the closed-loop control system, the micro-robot was driven to move along the square path and S-shape path, respectively. The performance of the Kalman filter was evaluated by observing the trajectory and error of the micro-robot.
First, we drove the micro-robot to move along a square path with a side length of 10 mm. Figure 18 shows the stages of the movement process of the micro-robot.
Figure 19a shows the trajectory of the micro-robot moving along the square path. The important role that the Kalman filter played in the control system is evident in that the motion of the micro-robot was significantly smoother and the error was smaller after the Kalman filter was added as the position predictor.
Figure 19b and c respectively show the errors in the x-axis direction and the y-axis direction. It is obvious that the errors in the two directions decreased significantly under the Kalman filter. When the Kalman filter was applied in the control system, the x-axis direction error was reduced by about 57.3%, and the y-axis direction error was reduced by about 23.5%.
Figure 20 shows the process of the micro-robot moving along the S-shape path. Figure 21a shows the trajectory of the micro-robot moving along the S-shape path. Use of the Kalman filter in the control system is apparent in these figures, as the movement of the micro-robot is more smooth and lacks excessive jitter. Using the filter, the overall dynamic performance of system was improved.
Figure 21b and c show the errors in the x-axis direction and the y-axis direction, respectively. Use of the Kalman filter in the control system resulted in the errors in the x-axis direction being reduced by about 53.1% and the errors in the y-axis direction being reduced by about 22.5%.

6. Conclusions

This paper demonstrates the use of combined coils to drive a magnetic micro-robot on a two-dimensional plane, for which we proposed a time-optimal control strategy based on auto disturbance rejection control technology. Compared with traditional PID controllers, the auto disturbance rejection control had better dynamic performance, and its effectiveness was verified by simulation analysis. The Kalman filter algorithm was used to predict the position of the micro-robot. Analysis of the experimental results shows that use of the Kalman filter algorithm substantially reduced the tracking errors of the micro-robot. Future studies will focus on the motion control of micro-robots in three-dimensional space and attempt experiments in more complex environments or in the context of living organisms.

Supplementary Materials

The following supporting information can be downloaded at: www.mdpi.com/article/10.3390/mi13122144/s1.

Author Contributions

Q.F., J.L. and J.Q. designed and performed the experiments, analyzed the data, and wrote the paper. J.J. contributed to the polish of the manuscript and provided insightful comments. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported in part by Six Talent Peaks Project in Jiangsu Province (grant numbers: GDZB-138) and Research Collaboration Project in Jiangsu Province (grant numbers: 2018-178) and National First-class Discipline Program of Food Science and Technology (grant numbers: JUFSTR20180302), and in part by the National Natural Science Foundation of China under Grant 52103283 and the Open Fund (Grant No. FM-202101) of Jiangsu Key Laboratory of Advanced Food Manufacturing Equipment & Technology.

Data Availability Statement

The datasets used or reported in the current work are available from the corresponding author on reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wei, Y.-K.; Li, Y.-J.; Zhang, Y.; Luo, X.-T.; Li, C.-J. Corrosion resistant nickel coating with strong adhesion on AZ31B magnesium alloy prepared by an in-situ shot-peening-assisted cold spray. Corros. Sci. 2018, 138, 105–115. [Google Scholar] [CrossRef]
  2. Nakadate, R.; Iwasa, T.; Onogi, S.; Arata, J.; Oguri, S.; Okamoto, Y.; Akahoshi, T.; Eto, M.; Hashizume, M. Surgical Robot for Intraluminal Access: An Ex Vivo Feasibility Study. Cyborg Bionic Syst. 2020, 2020, 1–9. [Google Scholar] [CrossRef]
  3. Bolopion, A.; Xie, H.; Haliyo, D.S.; Regnier, S. Haptic Teleoperation for 3-D Microassembly of Spherical Objects. IEEE/ASME Trans. Mechatron. 2010, 17, 116–127. [Google Scholar] [CrossRef]
  4. Kim, S.; Qiu, F.; Kim, S.; Ghanbari, A.; Moon, C.; Zhang, L.; Nelson, B.J.; Choi, H. Fabrication and characterization of magnetic microrobots for three dimensional cell culture and targeted transportation. Adv. Mater. 2013, 25, 5863–5868. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Sitti, M.; Ceylan, H.; Hu, W.; Giltinan, J.; Turan, M.; Yim, S.; Diller, E. Biomedical Applications of Untethered Mobile Milli/Microrobots. Proc. IEEE 2015, 103, 205–224. [Google Scholar] [CrossRef]
  6. Fong, B. Self-Cognizant Bionic Liquid Sensor for Pathogen Diagnosis. Cyborg Bionic Syst. 2021, 2021, 1–8. [Google Scholar] [CrossRef]
  7. Suter, M.; Zhang, L.; Siringil, E.C.; Peters, C.; Luehmann, T.; Ergeneman, O.; Peyer, K.E.; Nelson, B.J.; Hierold, C. Superparamagnetic microrobots: Fabrication by two-photon polymerization and biocompatibility. Biomed. Microdevices 2013, 15, 997–1003. [Google Scholar] [CrossRef]
  8. Kye, H.G.; Park, B.S.; Lee, J.M.; Song, M.G.; Song, H.G.; Ahrberg, C.D.; Chung, B.G. Dual-neodymium magnet-based microfluidic separation device. Sci. Rep. 2019, 9, 1–10. [Google Scholar] [CrossRef] [Green Version]
  9. Garg, A.A.; Jones, T.H.; Moss, S.M.; Mishra, S.; Kaul, K.; Ahirwar, D.K.; Ferree, J.; Kumar, P.; Subramaniam, D.; Ganju, R.K.; et al. Electromagnetic fields alter the motility of metastatic breast cancer cells. Commun. Biol. 2019, 2, 1–16. [Google Scholar] [CrossRef] [Green Version]
  10. Huang, J.-P.; Ge, X.-H.; Xu, J.-H.; Luo, G.-S. Controlled formation and coalescence of paramagnetic ionic liquid droplets under magnetic field in coaxial microfluidic devices. Chem. Eng. Sci. 2016, 152, 293–300. [Google Scholar] [CrossRef]
  11. Mahoney, A.W.; Abbott, J.J. Five-degree-of-freedom manipulation of an untethered magnetic device in fluid using a single permanent magnet with application in stomach capsule endoscopy. Int. J. Robot. Res. 2016, 35, 129–147. [Google Scholar] [CrossRef]
  12. Ohuchida, K. Robotic Surgery in Gastrointestinal Surgery. Cyborg Bionic Syst. 2020, 2020, 1–7. [Google Scholar] [CrossRef]
  13. Kratochvil, B.E.; Kummer, M.P.; Abbott, J.J.; Borer, R.; Ergeneman, O.; Nelson, B.J. OctoMag: An electromagnetic system for 5-DOF wireless micromanipulation. IEEE Trans. Robot. 2010, 26, 1006–1017. [Google Scholar]
  14. Heunis, C.; Sikorski, J.; Misra, S. Flexible Instruments for Endovascular Interventions: Improved Magnetic Steering, Actuation, and Image-Guided Surgical Instruments. IEEE Robot. Autom. Mag. 2018, 25, 71–82. [Google Scholar] [CrossRef] [Green Version]
  15. Wang, Y.; Li, W.; Togo, S.; Yokoi, H.; Jiang, y. Survey on Main Drive Methods Used in Humanoid Robotic Upper Limbs. Cyborg Bionic Syst. 2021, 2021, 1–12. [Google Scholar] [CrossRef]
  16. Kharboutly, M.; Gauthier, M.; Chaillet, N. Modeling the trajectory of a microparticle in a dielectrophoresis device. J. Appl. Phys. 2009, 106, 114312. [Google Scholar] [CrossRef] [Green Version]
  17. Kosa, G.; Shoham, M.; Zaaroor, M. Propulsion Method for Swimming Microrobots. IEEE Trans. Robot. 2007, 23, 137–150. [Google Scholar] [CrossRef]
  18. Liew, L.-A.; Bright, V.; Dunn, M.; Daily, J.; Raj, R. Development of SiCN Ceramic Thermal Actuators; IEEE: New York, NY, USA, 2003. [Google Scholar]
  19. Hwang, G.; Braive, R.; Couraud, L.; Cavanna, A.; Abdelkarim, O.; Robert-Philip, I.; Beveratos, A.; Sagnes, I.; Haliyo, S.; Régnier, S. Electro-osmotic propulsion of helical nanobelt swimmers. Int. J. Robot. Res. 2011, 30, 806–819. [Google Scholar] [CrossRef]
  20. Wang, R.; Chow, Y.T.; Chen, S.; Ma, D.; Luo, T.; Tan, Y.; Sun, D. Magnetic Force-Driven in Situ Selective Intracellular Delivery. Sci. Rep. 2018, 8, 14205. [Google Scholar] [CrossRef] [Green Version]
  21. Martel, S.; Felfoul, O.; Mathieu, J.-B.; Chanu, A.; Tamaz, S.; Mohammadi, M.; Mankiewicz, M.; Tabatabaei, N. MRI-based Medical Nanorobotic Platform for the Control of Magnetic Nanoparticles and Flagellated Bacteria for Target Interventions in Human Capillaries. Int. J. Robot. Res. 2009, 28, 1169–1182. [Google Scholar] [CrossRef]
  22. Gao, W.; Pei, A.; Wang, J. Water-Driven Micromotors. ACS Nano 2012, 6, 8432–8438. [Google Scholar] [CrossRef] [PubMed]
  23. Ongaro, F.; Pane, S.; Scheggi, S.; Misra, S. Design of an Electromagnetic Setup for Independent Three-Dimensional Control of Pairs of Identical and Nonidentical Microrobots. IEEE Trans. Robot. 2018, 35, 174–183. [Google Scholar] [CrossRef] [Green Version]
  24. Xu, T.; Gao, W.; Xu, L.-P.; Zhang, X.; Wang, S. Fuel-Free Synthetic Micro-/Nanomachines. Adv. Mater. 2016, 29, 1603250. [Google Scholar] [CrossRef] [PubMed]
  25. Fahrni, F.; Prins, M.W.J.; van Ijzendoorn, L.J. Micro-fluidic actuation using magnetic artificial cilia. Lab a Chip 2009, 9, 3413–3421. [Google Scholar] [CrossRef] [PubMed]
  26. Carpi, F.; Pappone, C. Stereotaxis Niobe® magnetic navigation system for endocardial catheter ablation and gastrointestinal capsule endoscopy. Expert Rev. Med Devices 2009, 6, 487–498. [Google Scholar] [CrossRef]
  27. Wang, X.; Luo, M.; Wu, H.; Zhang, Z.; Liu, J.; Xu, Z.; Johnson, W.; Sun, Y. A Three-Dimensional Magnetic Tweezer System for Intraembryonic Navigation and Measurement. IEEE Trans. Robot. 2017, 34, 240–247. [Google Scholar] [CrossRef]
  28. Oulmas, A.; Andreff, N.; Regnier, S. Closed-Loop 3D Path Following of Scaled-Up Helical Microswimmers; IEEE: New York, NY, USA, 2016. [Google Scholar]
  29. Niu, F.; Li, J.; Ma, W.; Yang, J.; Sun, D. Development of an Enhanced Electromagnetic Actuation System with Enlarged Workspace. IEEE/ASME Trans. Mechatronics 2017, 22, 2265–2276. [Google Scholar] [CrossRef]
  30. Jing, W.; Chowdhury, S.; Cappelleri, D. Magnetic mobile microrobots for mechanobiology and automated biomanipulation. In Microbiorobotics; Elsevier: Amsterdam, The Netherlands, 2017; pp. 197–219. [Google Scholar]
  31. Mohammadi, A.; Samsonas, D.; Leong, F.; Tan, Y.; Thiruchelvam, D.; Valdastri, P.; Oetomo, D. Modeling and control of local electromagnetic actuation for robotic-assisted surgical devices. IEEE/ASME Trans. Mechatron. 2017, 22, 2449–2460. [Google Scholar] [CrossRef]
  32. Salerno, M.; Firouzeh, A.; Paik, J. A low profile electromagnetic actuator design and model for an origami parallel platform. J. Mech. Robot. 2017, 9, 041005. [Google Scholar] [CrossRef]
  33. Yu, J.; Zhao, Y.; Chen, G.; Gu, Y.; Wang, C.; Huang, S. Realizing Controllable Physical Interaction Based on an Electromagnetic Variable Stiffness Joint. J. Mech. Robot. 2019, 11, 054501. [Google Scholar] [CrossRef]
  34. Santamato, G.; Chiaradia, D.; Solazzi, M.; Frisoli, A. A Lightweight Robotic Device Based on a Micro-Macro Actuation Concept for the Inspection of Railway Pantograph. J. Mech. Robot. 2020, 12, 061002. [Google Scholar] [CrossRef]
  35. Hu, Y. Self-Assembly of DNA Molecules: Towards DNA Nanorobots for Biomedical Applications. Cyborg Bionic Syst. 2021, 2021, 1–3. [Google Scholar] [CrossRef]
  36. Dong, D.; Lam, W.S.; Sun, D. Electromagnetic Actuation of Microrobots in a Simulated Vascular Structure with a Position Estimator Based Motion Controller. IEEE Robot. Autom. Lett. 2020, 5, 6255–6261. [Google Scholar] [CrossRef]
  37. Jeong, S.; Choi, H.; Cha, K.; Li, J.; Park, J.-O.; Park, S. Enhanced locomotive and drilling microrobot using precessional and gradient magnetic field. Sensors Actuators A Phys. 2011, 171, 429–435. [Google Scholar] [CrossRef]
  38. Salehi, M.; Pishkenari, H.N.; Zohoor, H. Position control of a wheel-based miniature magnetic robot using neuro-fuzzy network. Robotica 2022, 40, 3895–3910. [Google Scholar] [CrossRef]
  39. Yang, L.; Zhang, Y.; Wang, Q.; Chan, K.-F.; Zhang, L. Automated Control of Magnetic Spore-Based Microrobot Using Fluorescence Imaging for Targeted Delivery with Cellular Resolution. IEEE Trans. Autom. Sci. Eng. 2019, 17, 490–501. [Google Scholar] [CrossRef]
  40. Arcese, L.; Fruchard, M.; Ferreira, A. Adaptive Controller and Observer for a Magnetic Microrobot. IEEE Trans. Robot. 2013, 29, 1060–1067. [Google Scholar] [CrossRef] [Green Version]
  41. Jang, B.; Nam, J.; Lee, W.; Jang, G. A Crawling Magnetic Robot Actuated and Steered via Oscillatory Rotating External Magnetic Fields in Tubular Environments. IEEE/ASME Trans. Mechatronics 2017, 22, 1465–1472. [Google Scholar] [CrossRef]
  42. Gao, Y.; Huang, X.; Mann, I.S.; Su, H.-J. A Novel Variable Stiffness Compliant Robotic Gripper Based on Layer Jamming. J. Mech. Robot. 2020, 12, 051013. [Google Scholar] [CrossRef]
  43. Han, J. From PID to Active Disturbance Rejection Control. IEEE Trans. Ind. Electron. 2009, 56, 900–906. [Google Scholar] [CrossRef]
  44. Nelson, B.J.; Kaliakatsos, I.K.; Abbott, J.J. Microrobots for Minimally Invasive Medicine. Annu. Rev. Biomed. Eng. 2010, 12, 55–85. [Google Scholar] [CrossRef] [PubMed]
  45. Diller, E.; Sitti, M. Micro-scale mobile robotics. Found. Trends Robot. 2013, 2, 143–259. [Google Scholar] [CrossRef]
  46. Yu, C.; Kim, J.; Choi, H.; Choi, J.; Jeong, S.; Cha, K.; Park, J.-O.; Park, S. Novel electromagnetic actuation system for three-dimensional locomotion and drilling of intravascular microrobot. Sensors Actuators A Phys. 2010, 161, 297–304. [Google Scholar] [CrossRef]
  47. Cheang, U.K.; Kim, M.J. Self-assembly of robotic micro- and nanoswimmers using magnetic nanoparticles. J. Nanoparticle Res. 2015, 17, 1–11. [Google Scholar] [CrossRef]
  48. Song, S.; Song, S.; Meng, M.Q.H. Electromagnetic Actuation System Using Stationary Six-Pair Coils for Three-Dimensional Wireless Locomotive Microrobot; IEEE: New York, NY, USA, 2017. [Google Scholar]
  49. Li, D.; Niu, F.; Li, J.; Li, X.; Sun, D. Gradient-Enhanced Electromagnetic Actuation System with a New Core Shape Design for Microrobot Manipulation. IEEE Trans. Ind. Electron. 2019, 67, 4700–4710. [Google Scholar] [CrossRef]
  50. Pankhurst, Q.A.; Connolly, J.; Jones, S.K.; Dobson, J. Applications of magnetic nanoparticles in biomedicine. J. Phys. D Appl. Phys. 2003, 36, R167–R181. [Google Scholar] [CrossRef]
Figure 1. Combined coil structure.
Figure 1. Combined coil structure.
Micromachines 13 02144 g001
Figure 2. (a) Helmholtz coil driving principles. (b) Maxwell coil driving principles.
Figure 2. (a) Helmholtz coil driving principles. (b) Maxwell coil driving principles.
Micromachines 13 02144 g002
Figure 3. Schematic diagram of coils drive board.
Figure 3. Schematic diagram of coils drive board.
Micromachines 13 02144 g003
Figure 4. Block diagram of auto disturbance rejection controller.
Figure 4. Block diagram of auto disturbance rejection controller.
Micromachines 13 02144 g004
Figure 5. Schematic diagram of micro-robot planar path-tracking control.
Figure 5. Schematic diagram of micro-robot planar path-tracking control.
Micromachines 13 02144 g005
Figure 6. Schematic diagram of micro-robot path tracking.
Figure 6. Schematic diagram of micro-robot path tracking.
Micromachines 13 02144 g006
Figure 7. x-axis Helmholtz coil simulation. (a) The modulus of magnetic flux density. (b) The vector of magnetic flux density. (c) Magnetic induction value distribution in the working space.
Figure 7. x-axis Helmholtz coil simulation. (a) The modulus of magnetic flux density. (b) The vector of magnetic flux density. (c) Magnetic induction value distribution in the working space.
Micromachines 13 02144 g007
Figure 8. Magnetic field gradient of x-axis Helmholtz coil.
Figure 8. Magnetic field gradient of x-axis Helmholtz coil.
Micromachines 13 02144 g008
Figure 9. Three-axis Helmholtz coils magnetic flux density distribution.
Figure 9. Three-axis Helmholtz coils magnetic flux density distribution.
Micromachines 13 02144 g009
Figure 10. x-axis Maxwell coil simulation. (a) The modulus of magnetic flux density. (b) The vector of magnetic flux density. (c) Magnetic induction value distribution in the working space.
Figure 10. x-axis Maxwell coil simulation. (a) The modulus of magnetic flux density. (b) The vector of magnetic flux density. (c) Magnetic induction value distribution in the working space.
Micromachines 13 02144 g010
Figure 11. Magnetic field gradient of x-axis Maxwell coil.
Figure 11. Magnetic field gradient of x-axis Maxwell coil.
Micromachines 13 02144 g011
Figure 12. Single-axis combined coils simulation. (a) The modulus of magnetic flux density. (b) The vector of magnetic flux density. (c) Magnetic induction value distribution in the working space.
Figure 12. Single-axis combined coils simulation. (a) The modulus of magnetic flux density. (b) The vector of magnetic flux density. (c) Magnetic induction value distribution in the working space.
Micromachines 13 02144 g012
Figure 13. (a) Helmholtz x, y, z axes under 4 A current. Maxwell x, y, z axes under −4 A current. (b) Helmholtz x, z axes under −4 A current, Helmholtz y-axis under 4 A current. Maxwell x, z axes under 4 A current, Maxwell y-axis under −4 A current. (c) Helmholtz x, y, z axes under −4 A current, Maxwell x, y, z axes under 4 A current. (d) Helmholtz x, z axes under 4 A current, Helmholtz y-axis under −4 A current. Maxwell x, z axes under −4 A current, Maxwell y-axis under 4 A current.
Figure 13. (a) Helmholtz x, y, z axes under 4 A current. Maxwell x, y, z axes under −4 A current. (b) Helmholtz x, z axes under −4 A current, Helmholtz y-axis under 4 A current. Maxwell x, z axes under 4 A current, Maxwell y-axis under −4 A current. (c) Helmholtz x, y, z axes under −4 A current, Maxwell x, y, z axes under 4 A current. (d) Helmholtz x, z axes under 4 A current, Helmholtz y-axis under −4 A current. Maxwell x, z axes under −4 A current, Maxwell y-axis under 4 A current.
Micromachines 13 02144 g013
Figure 14. (a) Step-response simulation of Helmholtz coil. (b) Step-response simulation of Maxwell coil. (c) Impulse-response simulation of Helmholtz coil. (d) Impulse-response simulation of Maxwell coil. (e) Sinusoidal-response simulation of Helmholtz coil. (f) Sinusoidal-response simulation of Maxwell coil.
Figure 14. (a) Step-response simulation of Helmholtz coil. (b) Step-response simulation of Maxwell coil. (c) Impulse-response simulation of Helmholtz coil. (d) Impulse-response simulation of Maxwell coil. (e) Sinusoidal-response simulation of Helmholtz coil. (f) Sinusoidal-response simulation of Maxwell coil.
Micromachines 13 02144 g014
Figure 15. Simulation response of micro-robot tracking two signals. (a) Micro-robot tracking step signal. (b) Micro-robot tracking sinusoidal signal.
Figure 15. Simulation response of micro-robot tracking two signals. (a) Micro-robot tracking step signal. (b) Micro-robot tracking sinusoidal signal.
Micromachines 13 02144 g015
Figure 16. Experimental environment settings.
Figure 16. Experimental environment settings.
Micromachines 13 02144 g016
Figure 17. Steering tracking.
Figure 17. Steering tracking.
Micromachines 13 02144 g017
Figure 18. Square movement process.
Figure 18. Square movement process.
Micromachines 13 02144 g018
Figure 19. (a) Square trajectory. (b) Motion errors in x-axis. (c) Motion errors in y-axis.
Figure 19. (a) Square trajectory. (b) Motion errors in x-axis. (c) Motion errors in y-axis.
Micromachines 13 02144 g019
Figure 20. ‘S’ movement process.
Figure 20. ‘S’ movement process.
Micromachines 13 02144 g020
Figure 21. (a) ‘S’ motion trajectory. (b) Motion errors in x-axis. (c) Motion errors in y-axis.
Figure 21. (a) ‘S’ motion trajectory. (b) Motion errors in x-axis. (c) Motion errors in y-axis.
Micromachines 13 02144 g021
Table 1. Parameters of combined coils.
Table 1. Parameters of combined coils.
ParameterHXMXHYMYHZMZ
resistance (Ω)2.311.51.33.00.81.2
Inductance (mH)4.397207.21.51414.820.12080.511
Coil average distance (mm)1702941001735087
Equivalent radius (mm)1701701001005050
Coil outer diameter (mm)400400230230120120
Coil inner diameter (mm)2482481511516666
Number of turns52449311561439
Single copper wire diameter (mm)1.41.41.41.41.41.4
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Fan, Q.; Lu, J.; Jia, J.; Qu, J. 2D Magnetic Manipulation of a Micro-Robot in Glycerin Using Six Pairs of Magnetic Coils. Micromachines 2022, 13, 2144. https://0-doi-org.brum.beds.ac.uk/10.3390/mi13122144

AMA Style

Fan Q, Lu J, Jia J, Qu J. 2D Magnetic Manipulation of a Micro-Robot in Glycerin Using Six Pairs of Magnetic Coils. Micromachines. 2022; 13(12):2144. https://0-doi-org.brum.beds.ac.uk/10.3390/mi13122144

Chicago/Turabian Style

Fan, Qigao, Jiawei Lu, Jie Jia, and Juntian Qu. 2022. "2D Magnetic Manipulation of a Micro-Robot in Glycerin Using Six Pairs of Magnetic Coils" Micromachines 13, no. 12: 2144. https://0-doi-org.brum.beds.ac.uk/10.3390/mi13122144

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop