Abstract
In the biomedical field, a wireless microrobot in a pipe which can move smoothly in water or other aqueous mediums has been urgently demanded. In this paper, several methods of designing a novel microrobot with symmetrical motion characteristics have been discussed and a new kind of wireless microrobot has been developed. According to the modelling analysis, we considered two kinds of common cases occurring in vertical motion, which required gravity compensation. Based on two groups of simulations and experiments on forward-backward motion, upward-downward motion and inclined plane motion, the results and dynamic error evaluation indicated that the wireless microrobot with symmetrical structure could realize similar kinematic characteristics in the horizontal motion. The gravity compensation played an important role in the design process, and the performance of the vertical motion had been improved by gravity compensation. With this method, we made the wireless microrobot realize symmetrical motion characteristics, and simplified the control strategies. Finally, a control panel for our system was designed, which could control the current motion states more intuitively and far more easily through the buttons. The developed wireless microrobot would be very useful in the industrial application and microsurgery application.
1. Introduction
Endovascular intervention has become more and more popular in recent times. The study of microrobots in a pipe is an important research branch of MEMS (Micro-Electro-Mechanical System). Studies on the locomotion methods of microrobots which move within organic tubes in the human body can be used in practical medical fields, such as diagnosis and treatment of body organs. The traditional endoscopy uses a long, thin and flexible tube which is required for connecting the endoscope with the workstation outside the patient's body [1]. However, the traditional endoscope has two problems, in that they cause pain to the patients and are limited to the view of the small intestine. Therefore, some techniques are required that thoroughly improve the capabilities of the endoscopy field. To minimize the suffering of patients, a capsule endoscope has been developed [2]. With the aid of peristalsis, the capsule moves passively through the GI tract. However, the capsule cannot control its moving direction and moving speed itself, due to the loss of large quantities of data. So the capsule endoscope needs an autonomous moving function.
A number of mechanisms for an actuator have been proposed, which are used for providing locomotion in the pipe, such as piezoelectric elements [3, 4], air cylinders [5], electrorheological fluids [6], shape memory alloys [7, 8], electromagnetic motors [9–14] and globular magnetic actuators capable of locomotion in a pipe by a combination of mechanical vibration and electromagnetic force [15, 16].
Many studies on the microrobot and microswimmer also had been developed. In 2013, Jian Feng and Sung Kwon Cho proposed a kind of two-dimensionally steering microswimmer propelled by oscillating bubbles [17]. In 2012, Yan-Hom Li and He-Ching Lin proposed a novel steering technology of magnetic micro-swimmers under an oscillating field [18]. In 2011, J. Scogna and J. Olkowski proposed a kind of microswimmer which could be controlled in one dimension for drug delivery and the size was only 7–8 μm [19]. In 2011, Li Zhang and Tristan Petit proposed rotating ferromagnetic nanowires (NWs) which were capable of self-propulsion near a solid boundary using a uniform rotating magnetic field [20].
In our previous study, wireless microrobots controlled by an external magnetic field were safe and reliable, and could be carried deep within the tissues of living organisms in the body fluids. Thus, a new kind of magnetic field model to easily drive the microrobot is required [21–22]. With the different kinds of external magnetic field driven microrobots, many locomotion forms are proposed as shown in Figure 1. For example, the fish-like motion [23], paddle motion [24], propeller-driven motion [25] and spiral motion [26–27]. In many kinds of motions, we identified the kind of spiral motion which not only could obtain the maximum driving force, but also had the highest efficiency in a very small space [28].

Several kinds of wireless in-pipe microrobots
However, the control algorithm is high level, non-linear and very complex [29], which is mainly due to the asymmetry of kinematic performance. It is also difficult to realize bidirectional motion and upward motion. So in this paper, we focus on developing a kind of wireless microrobot with symmetrical motion characteristics and simplify the control strategies.
2. Rotation Magnetic Field Model
2.1 Design of the Rotation Magnetic Field
The rotation magnetic field model was shown in Figure 2, which contained two pairs of coils. When we inputted the currents with the phase difference of 90° in each coil, the rotation magnetic field was generated.

Rotation magnetic field model
We assumed that the currents flowing through the coil M1 and coil M2 had the same amplitude, but the phase difference was 90°, which was shown in Figure 2(a). Then the strengths of a magnetic field generated by the coil M1 and M2 were given by:
In order to conveniently analyse the magnetic field, we assumed B10=B20=B0. The rotation vector directions of B1 and B2 were perpendicular to each other, so the absolute value and the direction of the magnetic field after the superposition of the two magnetic field vectors were:
From the Equation 3 and Equation 4, we could find the value of the synthesis magnetic field strength was a constant, and it was equal to the maximum value B1 (or B2). In addition to that, the direction of B changed along with the periodic time changing. So the synthesis of vector B was a rotation vector, and the angular velocity of rotation was ω, which was a rotation magnetic field.
2.2 Parameter of the Rotation Magnetic Field
In our study, we needed to establish a rotation magnetic field to drive the developed microrobot. Based on the magnetic field theory discussed, we had designed a new external superposition magnetic field in horizontal and vertical directions. The experimental model of the rotation magnetic field was shown in Figure 2(c).
We set the two pairs of coil M1 and coil M2 at 90° to each other. In order to generate the rotation magnetic field, we made the solenoid by using enamelled wire. In particular, we inserted the magnetic cores to enhance the magnetic flux density. It was a simple design to obtain the rotation magnetic field, and the specification of each coil was shown in Table 1.
Specification of each coil
Following the Biot-Savart Law, the relationship between the magnetic flux density H and current i could be measured by a gauss meter, which was shown in Figure 3, and Equation (5) of the curve could be fitted as well:

Relationship between the magnetic density and current
3. Design of the Microrobot Structure
3.1 The Functional Analysis of the Microrobot
When considering practical application, the microrobot developed should realize the following basic movements:
Moving forward and backward in the horizontal direction
Moving upward and downward in the vertical direction
Stopping and running at any position we need
The microrobot was mainly driven and controlled by an external magnetic field, and the properties of the magnetic field were generated by the current input coil. In order to simplify and optimize the control strategies, the microrobot should have symmetrical motion characteristics.
The symmetrical motion characteristics can be obtained by the most intuitive method, which is the symmetry of the mechanical structure and functions. In addition, as the spiral motion had the highest efficiency in a very small space [21], we should fully consider the symmetrical mechanical structure and the spiral structure.
3.2 The Structure of the Microrobot
The spiral type of microrobot had been designed in our previous study [30], which contained four pieces of permanent magnets and the rotation mechanism was separate with the main body of the microrobot. Based on the analysis aforementioned, we had proposed a new kind of spiral type of microrobot, which only contained one piece of permanent magnet, so it had a more compact structure [31].
Figure 4 s the structure of the prototype of the wireless microrobot. It consists of two main parts. In order to ensure the balance of the microrobot, we adopted a completely symmetrical structure, which also could achieve the same kinematic characteristics. Regardless of the posture and position of the microrobot in the pipe, the centre of mass could always be fixed in the centre of the microrobot (the position of the permanent magnet) as shown in Figure 5.

The structure of the proposed microrobot

The centre of mass distribution
3.3 The Dynamic Analysis of the Microrobot
With two input currents flowing into the coils, a rotation magnetic field can be generated; the kinematic characteristic of the microrobot are shown in Figure 6. The permanent magnet would rotate following the rotation of the magnetic field. With its spiral structure, the microrobot could thus move forward. By changing the order of the two-phase input currents, we can obtain an opposite direction of rotation of the magnetic field immediately. Then the microrobot would obtain an opposite direction of the driving force and move towards the opposite direction. The moving speed of the microrobot could be controlled by adjusting the frequency of the input currents.

Kinematic characteristics of the proposed microrobot
To propel the wireless microrobot, the propulsive force and torque was provided by rotation of the magnetic field. We could change the magnetic force and torque by following Equation 6 and 7 to overcome fluid resistance in the pipe:
where μ0 is the permeability of vacuum, V represents the volume of the body, M is the average magnetization, ∇ represents a gradient operator, H is the magnetic flux density, B is the magnetic induction intensity.
Due to the symmetrical structure, the propulsive force along the axis could be described by a symmetric propulsion matrix as shown in Equation 8 [32]:
where, f is the non-fluidic applied force, τ is the non-fluidic applied torque, n is the number of thread units, R and θ are as shown in Figure 8, k1 and k2 are the constants, which are the viscous drag coefficients for the wireless capsule microrobot along the axis, υ is forward velocity and ω is angular speed.

The modelling of the microrobot in the horizontal direction

The modelling in vertical direction
When the θ was equal to 0 or π, the first and last term in Equation 8 would be zero.
3.4 Modelling and Analysis of the Microrobot
3.4.1 The Modelling in Horizontal Direction
According to the Right-hand Rule, we established the coordinate system for the microrobot movement, which is shown in Figure 7.
In the right-hand coordinate system, Va was the axial velocity, Uc was the circumferential velocity. In order to simplify the force analysis of the microrobot, we assumed the direction of the screw thread was the x-axis, the thickness direction of the screw thread was the y-axis, the direction which was vertical to the microrobot body was the z-axis, and θ was the spiral angle. Uc and Va could be composed of ux (x-component) and ωz (z-component).
Based on Hydromechanical Lubrication Theory and Newton Viscous Law, similar force analysis has been discussed [33], [34]. According to the microrobot we designed in the pipe, the main force would be evaluated as shown in Table 2 in the horizontal direction, which consisted of fx, fz, Faxial, Fcircum fl-resist and fr-fricion.
Explanation of the main force
According to Newton Viscous Law, we could obtain the fx and fz as following:
where n is the number of thread unit, μ is the viscosity coefficient of the liquid, A is the contact area, du/dy is the velocity gradient.
Then we could compute the axial traction force, circumferential friction force and resistance from liquid as following:
where k is the damping coefficient of the liquid, ρ is the density of the liquid, S is the front of the contact area of the microrobot, η is the coefficient of rolling friction, G is the weight of the microrobot, Vg is the volume of the microrobot, ρg is the density of material used to design the microrobot. The mg can be obtained from an electronic balance easily, then G and Fu can be calculated.
In order to simplify the equations above, we assumed β=a/a+b, γ=h/c, so we could obtain the equations as:
Therefore, the total force of the microrobot in the horizontal direction could be expressed as:
3.4.2 The Modelling in Vertical Direction
In the vertical motion, we also established the coordinate system for the microrobot's movement as shown in Figure 8. The significant difference was the microrobot was also influenced by its gravity and the buoyancy from the liquid, and the G and Fu should be considered separately.
If we assumed ΔF=Fu-G, the total force of the microrobot in the vertical direction was:
As the effect of the F*r-friction is very small, for the sake of calculation, we assumed F*r-friction = Fr-friction. Then the Fvertical was only related to Fhorizontal and ΔF.
4. Simulation
Due to the symmetrical spiral structure, the microrobot had similar kinematic characteristics in the horizontal direction motion, so we would focus on the study of symmetry motion in the vertical direction.
As a result of the limitation of the external magnetic density, the driving force is also limited. Therefore, there would be two common cases occurring in the vertical direction and needing gravity compensation. One was the Fvertical, a negative value, and the other was Fhorizontal, a positive value which was very small. We designed two groups of experiments to illustrate the two cases, and evaluated the motion effects with different design parameters. Firstly, we simulated the performance in the horizontal direction. Secondly, we focused on the simulation in the vertical direction.
4.1 Simulation on Group I
In the first group, we simulated the performance effect of the spiral depth h to Faxial, Fcircum, f l-resist, f r-friction, Fhorizontal and Fvertical. Referring to the structure parameters of the microrobot we developed, as shown in Table 3; it could be assumed β=0.5, μ=1, u=0.005, v=0.05, n=4, k=1, ρ=1, R=0.005, θ=π/4. Considering the size of the microrobot, we fixed the value of γ, ranging from 0 to 10. The simulation results in the horizontal direction are shown in Figure 9 and Figure 10.
Specifications of the microrobot in first group

The curves of Faxial, Fcircum, f l-resist, f r-friction, Fhorizontal when γ changes

The curves of Fvertical when γ changes
From Figure 9, we can conclude that when h/c>=1 and faxial>0, the traction force of the capsule was positive, and when h/c=1, fcircum had the minimum value. We could change the value c and h by changing the volume of the microrobot in the different applications. Usually we set a suitable value c, and changed the value of h. The changing of the volume could be approximately neglected because of the small scale of the microrobot.
From Equation 22 and Figure 10, we can find that Fvertical was just related to Fhorizontal and ΔF, but the Fvertical was a negative value. So the microrobot could not realize upward motion, and needed gravity compensation.
4.2 Simulation on Group II
In the second group, we simulated the performance effect of the spiral angle θ to Faxial, Fcircum, f l-resist, f r-friction, Fhorizontal and Fvertical. Referring to the structure parameters of the microrobot designed in Table 5, we assumed β=0.5, γ=1, μ=1, u=0.005, v=0.05, n=4, k=1, ρ=1, R=0.003. Considering the size of the microrobot, we fixed the value of θ, ranging from 0° to 90°. The simulation results are demonstrated in Figure 11 and Figure 12.
Specifications of the materials
Specifications of the microrobots in the second group

The curves of Faxial, Fcircum, f l-resist, f r-friction, Fhorizontal when θ changes

The curves of Fvertical when θ changes
The simulation results provided the tendency of force applied on the microrobot, which indicated that when θ was from 10° to 80°, the Fvertical was a positive value, but it was a very small value. So the microrobot could achieve upward motion without gravity compensation in theory. When θ was 45°, the microrobot could obtain the maximum total force, so the spiral angle of the microrobot would be designed with 45°.
5. Experimentation
5.1 System Architecture
The experimental system architecture is shown in Figure 13. In the system, the function generator outputted sinusoidal control signals which had a phase difference of 90°. Then by using the amplifying unit, the control signals were able to drive the magnetic field model to generate a rotation of the magnetic field, which controlled the status of the microrobot. During the experiments, we could change the amplitude, frequency and phase of the control signals by adjusting the function generator.

Experimental system architecture
The amplifying unit was an extremely significant part in the system and the circuit is shown in Figure 14. With the aid of the amplifying unit, the weak voltage signal was converted to a strong current signal without any frequency and phase changing. The relationship of conversion could be obtained as following:

Circuit design of the amplifying unit
If we choose the R1=R2=1Ω, R3=R4=10Ω, Rref=10Ω, the output current would follow with the input signal:
5.2 Experiments on Group I
5.2.1 Forward and Backward Motion
In the detection process, the microrobot should complete forward-backward and upward-downward motion in the pipe. Based on the theoretical analysis of the proposed microrobot in our previous study [31], we had proved that it can move forward and backward in a pipe.
By using the measurement system shown in Figure 15, the following characteristics of the moving speed were measured. We carried out the experiments by changing the frequency of input currents from 0 Hz to 20 Hz, and the amplitude of input currents were fixed at 0.7 A.

Measurement system
In the first group, the prototype of the developed microrobot (Microrobot A) is shown in Figure 16, and the design parameters are shown in Table 3.

Prototype of the microrobot
We inserted a piece of axial permanent magnets in the centre of the microrobot body, for which the size was Φ5×2 mm, as shown in Figure 17.

Prototype of the permanent magnet
Based on the measurement system in Figure 15(a), the microrobot could perform a rotation motion. By changing the frequency and phase of input currents, the moving speed and direction could be adjusted. The experimental results of average moving speed of the microrobot that can be obtained are shown in Figure 18. We had also obtained the maximum moving speed, 22.7 mm/s around 12Hz.

Moving speed of forward-backward motion
5.2.2 Upward and Downward Motion
According to our preceding analysis and the simulation results shown in Figure 11, the Fvertical was a negative value, so the microrobot could not realize upward motion in theory. We should use a counterweight for the microrobot to make ΔF=0, which could allow the microrobot to suspend in water.
The material of the microrobot body was made of EVA (Ethylene Vinyl Acetate), therefore, we selected another material which had a smaller density as gravity compensation. The specifications of materials are shown in Figure 19 and Table 4.

The gravity compensation of the microrobot
In addition, the EPE (Expandable Polyethylene) was a kind of soft material, which would be installed at the head and end of the microrobot to buffer the microrobot from impacting with the pipeline. In this way, the safety in the detection process could be improved. The microrobot with gravity compensation (Microrobot B) is shown in Figure 16.
Based on the measurement shown system in Figure 16(b), the microrobot could achieve upward-downward motion in a pipe. The experimental results of average moving speed between the microrobot before and after the gravity compensation are shown in Figure 20.

Moving speed of upward-downward motion
The experimental results indicated that the microrobot with gravity compensation could realize upward motion easily, and the dynamic characteristics had been clearly improved, which had similar kinematic characteristics in the vertical motion. We had also obtained the maximum moving speed, 6.3 mm/s around 11Hz.
5.2.3 Inclined Plane Motion
In order to evaluate the kinematic characteristics of the microrobot further, we also conducted a set of experiments on the motion with an inclined plane of 45°. The experimental results of average moving speed of the microrobot could be obtained in Figure 21. This indicated that the microrobot with gravity compensation also had similar characteristics in inclined plane motion, and the maximum moving speed was 13.1mm/s around 12 Hz.

Moving speed of inclined plane at 45°
5.3 Experiments on Group II
5.3.1 Forward and Backward Motion
We reduced the size of the microrobot in order to establish the second case of gravity compensation. The prototype of this microrobot (Microrobot C) is shown in Figure 22, and the design parameters are shown in Table 5.

Prototype of the microrobots C and D
Similar to group I, by using the measurement system shown in Figure 15(a), the average moving speed of the microrobot was obtained, as shown in Figure 23. We had also obtained the maximum moving speed, 36.5mm/s around 14Hz.

Moving speed of forward-backward motion
5.3.2 Upward and Downward Motion
According to the simulation results in Figure 13, we could conclude that when θ was from 10° to 80°, the Fvertical was a positive value, albeit a very small one. The experimental results on upward motion shown in Figure 24 indicate that the microrobot could barely realize the upward motion with a very low speed. In order to simplify the control strategy, we should make the kinematic characteristics as symmetrical as possible, so we also need a counterweight for the microrobot. The microrobot with gravity compensation (Microrobot D) is shown in Figure 22.

Moving speed of upward-downward motion without gravity compensation
Based on the measurement system shown in Figure 15(b), the experimental results of average moving speed of the microrobot could be obtained as shown in Figure 25, and the maximum moving speed was 4.7mm/s around 16Hz.

Moving speed of upward-downward motion
5.3.3 Inclined Plane Motion
The experiment of inclined plane motion had also been taken. The experimental results of average moving speed of the microrobot are shown in Figure 26, and the maximum moving speed was 14.8mm/s around 15 Hz. We found that the microrobot with gravity compensation maintained similar kinematic characteristics in the inclined plane motion.

Moving speed of inclined plane at 45°
5.4 Discussion
We had simulated and experimented, based on the two cases, the effect on the motion performance in the horizontal direction, vertical direction (with or without gravity compensation) and inclined plane motion with gravity compensation. The results are shown in Figure 18, Figures 20–21 and Figures 23–26.
In order to gain a more intuitive evaluated result on the affect of gravity compensation in the vertical motion, we calculated the dynamic error as the following [35]:
The dynamic error analysis and comparison is shown in Figure 27. From Figure 18 and Figure 23, we could conclude that the kind of microrobot we designed realized similar kinematic characteristics in the horizontal direction motion. Furthermore, the ΔF's limited effect on the horizontal motion and the symmetrical structure ensured the centre of mass would always stay at the centre of the microrobot. Viewing the experimental results shown in Figures 20–21 and Figures 25–26, it can be seen that the gravity compensation mechanism would improve the kinematic characteristics in the vertical direction, and the microrobots with symmetrical structure also had similar kinematic characteristics in the vertical motion and inclined plane motion.

Dynamic error
In addition, we found that the moving speed of the microrobot would decline to zero in the region of high level frequencies in each experimental result. This was because when the applied magnetic field rotated faster than the microrobot's step-out frequency (the frequency requiring the entire available magnetic torque to maintain synchronous rotation), the velocity of the microrobot dramatically declined. So the wireless microrobots could not rotate following the rotation magnetic field at a high level frequency and the velocity would fall to zero [36].
So this new kind of microrobot with symmetrical structure and gravity compensation had the ideal dynamic characteristics for the forward-backward, upward-downward motion and inclined plane motion. The dynamic error of the microrobot with gravity compensation was far less than that of the microrobot without gravity compensation. These advantages would make the control strategies symmetrical and much simpler.
6. Design of control panel for the system
We used a function generator to obtain sine wave signals at the experimental stage, but the device was not portable and compact enough. Due to the symmetrical kinematic characteristics of the microrobot, the control strategies could become symmetrical and much simpler.
Based on the experimental results, the control panel we designed is shown in Figure 28. Through selection of the buttons, we could change the required status of the wireless microrobot. The MCU would output the corresponding signals to the amplifying unit, and then display the current status. The control panel could realize stopping, running and changing the direction of motion, which could be adjusted through the button control. In addition, in the running status, we could choose low-speed motion, mid-speed motion and high-speed motion according to the actual operation needs. The moving speed could be referenced by the experimental results.

Control panel (Button control: 1. Start/Stop, 2. Motion direction control, 3. Low-speed motion, 4. Mid-speed motion, 5. High-speed motion)
With this control panel, our system would be portable and compact, and we could select the microrobot motion states more intuitively and easily by using the buttons.
7. Conclusion and future work
In this paper, a method to design a microrobot with symmetrical motion characteristics have been discussed and designed. Through modelling analysis, we have considered two kinds of common cases occurring in the vertical motion, which would need gravity compensation. Based on two groups of experiments, the results and dynamic error evaluation have been obtained. This indicated that the wireless microrobot with symmetrical structure would realize similar kinematic characteristics in the horizontal motion, and the gravity compensation would play an important role in the design process. In addition, the kinematic performance of the motion had been improved by a gravity compensation mechanism, which provided the wireless microrobot with symmetrical motion characteristics and simplified the control strategies. At last, a control panel for our system was designed, which would make the current motion states more intuitive and much easier to operate.
The results of two group experiments have been summarized in Table 6. The results indicated that the symmetrical kinematic characteristics of the wireless microrobot in a pipe can be designed following the method we proposed. In addition, we could obtain the ideal dynamic characteristics by the gravity compensation mechanism, especially in the vertical direction motion. These advantages made the control strategies symmetrical and much simpler.
Experimental results
In the future, we will focus on establishing a more DOFs experimental system in order to realize the potential for greater microrobot DOFs motion in the pipe, and make it apply to more practical related fields.
Footnotes
8. Acknowledgments
This research is supported by the Key Research Program of the Natural Science Foundation of Tianjin (13JCZDJC26200) and the General Research Program of the Natural Science Foundation of Tianjin (13JCYBJC38600) and The Project-sponsored by SRF for ROCS, SEM.
