# An Adaptive Control Method for Ros-Drill Cellular Microinjector with Low-Resolution Encoder

Authors: Zhenyu Zhang; Nejat Olgac
Journal: Journal of Medical Engineering
 (2013)

---

## Abstract

A novel control methodology which uses a low-resolution encoder is presented for a cellular microinjection technology called the Ros-Drill (rotationally oscillating drill). It is developed primarily for ICSI (intracytoplasmic sperm injection) operations, with the objective of generating a desired oscillatory motion at the tip of a micro glass pipette. It is an inexpensive setup, which creates high-frequency (higher than 500 Hz) and small-amplitude (around 0.2 deg) rotational oscillations at the tip of an injection pipette. These rotational oscillations enable the pipette to drill into cell membranes with minimum biological damage. Such a motion control procedure presents no particular difficulty when it uses sufficiently precise motion sensors. However, size, costs, and accessibility of technology to the hardware components severely constrain the sensory capabilities. Consequently, the control mission and the trajectory tracking are adversely affected. This paper presents two contributions: (a) a dedicated novel adaptive feedback control method to achieve a satisfactory trajectory tracking capability. We demonstrate via experiments that the tracking of the harmonic rotational motion is achieved with desirable fidelity; (b) some important analytical features and related observations associated with the controlled harmonic motion which is created by the low-resolution feedback control structure.


---

## Body

## 1. Introduction and Motivation

We provide some background on the main task at hand and relevant motivation, before the control methodology is elaborated. ICSI (intracytoplasmic sperm injection) is a broadly utilized technique for artificial fertilization. This procedure is successfully performed in human oocytes as well as other species such as mouse and bovine. First, a holding pipette is used to immobilize an individual oocyte with a slight suction. Then an injection pipette (with outer diameter of about 8μm), which contains the sperm head to be injected, is forced into the cell. The piercing through the zona layer and the membrane needs to be achieved with minimal biological damage to facilitate rapid healing. A significant amount of research effort has been devoted towards developing microscopic instruments for ICSI from this perspective alone. The most popular procedure at the present is the piezo-assisted ICSI [1]. However, its piercing performance is successful only by using a small mercury droplet in the pipette tip [2, 3]. Without this addition, undesirable lateral oscillations occur at the tip and severely hamper the performance of piercing. Due to high toxicity of mercury, on the other hand, the piezo-assisted ICSI procedure is forbidden in many biological laboratories. In recent years, an improved remedial technology, called rotationally oscillating drill (Ros-Drill), is introduced [4]. This technique shows comparable results to those obtained by the piezo-assisted ICSI process, with one major difference, that Ros-Drill does not have the mercury problem [5].A schematic of the Ros-Drill assembly is shown in Figure1. The injection pipette is connected to a small-precision micromotor which is controlled to track a desired sinusoidal trajectory:
(1)θd=Adsin(2πfd·t),
where Ad is the oscillation amplitude (typically around 0.2 deg), and fd is the frequency (in the range of 500–700 Hz). These selections are based on a simple bandwidth analysis; at this frequency range of external stimuli the cell membrane is not expected to follow the pipette tip motion. The ensuing relative rotational motion between the pipette and the cell membrane creates a clean piercing action, which facilitates the rapid healing of membrane after piercing.(a) Black mice reproduced using Ros-Drill technology (white mouse is the surrogate mother) [5]. (b) Ros-Drill assembly and control system.
(a)
(b)Considerable amount of experimental effort has been invested to demonstrate the validity of Ros-Drill technology [4]. It is shown that the success rate in ICSI by Ros-Drill is comparable to that of piezo-assisted technology, provided that the pipette oscillations are maintained as close as possible to the desired trajectory in (1). Many healthy mouse offspring are produced using the Ros-Drill methodology, as shown in Figure 1(a). These biological tests are conducted by a group of experimentalists from the University of Connecticut and the University of California, Davis, USA. Two outstanding requirements are noted in these experiments towards an acceptable Ros-Drill performance: (i) rotational oscillation to track the harmonic trajectory very closely; (ii) the flexible pipette to be concentric to the rotational axis of Ros-Drill. The later condition is shown to satisfy because of the extreme bending compliance and whirling effects [2, 3]. The condition (i) is the topic of this paper. The smooth variation of the rotational motion in accordance with a harmonic function is the most natural desired trajectory. The contribution of this paper is a proper control law which can serve the objective under severe sensing constraints. Overarching restriction in this study comes from the pricing aspect. In order to make the Ros-Drill accepted by a wide range of IVF clinics, the cost of this automated device must remain within $1500 per copy; otherwise, it would be difficult to compete against the commonly-used but inefficient method adopted by trained ICSI specialists. These conditions set the tone of the key research issues.Let us take a closer look at the physics of the Ros-Drill microinjection procedure. The rotational stiffness of the pipette holder and the pipette including the extremely fine tip are assumed to be high. Consequently, the angular displacement,θ, of the injection pipette is transmitted from the shaft of the micromotor to the pipette tip without a loss. The lateral (bending) vibrations are the major concern at the tip, and the present operating scheme is intended to suppress them to ignorable levels especially when compared with those oscillations caused by the piezo-assisted ICSI [2, 3]. With these assumptions, the major objective of this study is directed to insure desirably precise harmonic motion tracking capabilities at the pipette tip [6–8] despite the very coarse sensory capabilities.For position servoing, in general, the sensors are expected to have high resolution vis-à-vis the range of the intended motion, which can yield a desirable tracking capability. For instance, in representing a harmonic trajectory, one expects to have a minimum of 10 discrete data points per cycle. However, resolution of digital encoders is limited by the number of slots on a rotating disk, through which the encoder’s light beam travels. Although advances in encoder technology have wonderfully progressed to increase the number of slots in order to improve the resolution of the encoders, the trade-off between resolution and cost is unavoidable. In such applications as Ros-Drill where the cost limitations are very stringent, sensor resolution is often compromised.We encounter considerable past research on control methods using some low-resolution sensors. Recent model-based speed observers [9] make the velocity estimate robust using an interesting disturbance observer. In another effort, Kwon et al. [10] incorporate acceleration measurement in velocity estimation and motion control. Bautista-Quintero and Pont [11] propose an H-infinity control algorithm for sensor-constrained mechatronic systems using the position sensors with relatively low resolutions. They demonstrate how this procedure allows a faithful reproduction of observed motion starting from a limited sensing ability using relatively common (and inexpensive) microcontrollers. Furthermore, the methods which deal with control applications with low-resolution sensors usually have high computational demands to compensate the sensory shortfall. This aspect quickly makes the method prohibitive in a cost-sensitive design such as Ros-Drill.We wish to familiarize the reader with the components of the first-generation Ros-Drill setup and the current improvements. The first-generation design having a 512 lines/revolution-encoder with quadrature signature is the finest selection we could find within the cost and spatial confines. Its resolution is 0.17 deg (including the quadrature signature feature). It allows a maximum of 2-step reading over the desired peak-to-peak stroke (note that the desired amplitude is 0.2 deg) [4, 12]. This sensor makes our best observation of a Ros-Drill harmonic cycle with a 2-step representation, which is a colossal handicap to perform the control. The spirit of the proposed control scheme and the focus of this paper are primarily on this crucial feature.Moving on to another component, the first trial generation of Ros-Drill employs a PLC (programmable logic controller) as its digital controller, which has a maximum sampling speed of 1 KHz. This constraint clearly limits the maximum frequency of the controlled trajectory to 500 Hz. Most recently reported ICSI tests use injection pipette oscillations up to 0.3 deg amplitude and maximum frequency of 500 Hz. These oscillations last a certain length of time which we name theduration of oscillation, D. Typically, D varies within 250–500 msec. The preliminary reports claim that Ros-Drill-assisted ICSI results in embryo survival, embryo development, birth, and weaning rates comparable to those of piezo-assisted ICSI using mercury [5]. Although these biological results are very promising, the trajectory tracking performance of the existing prototype is not satisfactory mainly because of the low resolution of position sensor and low control sampling rate of PLC. This work represents an effort to further improve this performance.In the second-generation (current) Ros-Drill microinjector prototypes, the low-sampling-frequency matter is considerably improved by replacing the PLC with a microcontroller (which brings the sampling speed from 1000 Hz to 10000 Hz) and the encoder resolution issue by selecting one-level higher-capacity sensor (1000 lines and 0.09 degree resolution after quadrature [12]). However, despite this upgrade, the position sensor with the resolution of 0.09 degrees still presents the biggest hurdle in this control system design because the required harmonic motion of 0.2 degree amplitude displays only 4-step peak-to-peak encoder recording, which is still low. In [12], a look-up-table-based adaptively tuned PID control law is used to cope with the limitations in the hardware. The most appropriate PID (proportional, integral, and derivative) control gains are a priori selected corresponding to the given operating frequency using dynamic simulations. For a given desired trajectory, θd, with Ad and fd attributes, an ad hoc 3D search routine is performed using the representative SIMULINK simulation program. This yields a set of feedback gains, which are then utilized through a look-up adaptation structure. They are regarded as nominal control gains as they are tuned based on the nominal model of Ros-Drill, and the control design is considered as the first adaptation stage towards the second-generation Ros-Drill.An important point to discuss for this first adaptation stage is that the look-up-table-based control law can handle the potential parameter uncertainties, disturbances, and the unmodeled mechanical properties of the system only via the supervisory interference from the user. These uncertainties are important and unavoidable. For instance, the external load originating from the membrane resistance torque may vary with the type of species and it may affect the control performance of the system. Still worse is the extremely stringent requirement of concentricity on the shaft of motor and pipette holder via a coupling and ensuing resistance to the servosystem. In order to achieve better updates of the control gains, we introduce a second adaptation scheme in this paper. This procedure constitutes the primary contribution of the paper. Furthermore, the low-resolution encoder is a range indicator (as discussed later) rather than a measurement device. We offer some novel observations on the actual angular motion.The text is organized as follows. In the first part of the paper, a dedicated novel adaptive feedback control method is presented to achieve a satisfactory trajectory tracking capability with poor-quality sensors. Some important definitions about encoder signal are listed in Section2. The first adaptation scheme is revisited in Section 3. In Section 4, the second-stage adaptive law is developed to update the control gains in situ, which is validated by experiments. In the second part of the paper, a stochastic perspective is presented which correlates the encoder readings with the actual angular motion. Section 5 presents an intriguing observation over the trajectory tracking process. Finally, conclusions are given in Section 6.


## 2. Some Descriptors of the Encoder Signal

The highlighted feature in this study is the unusually low-resolution sensor (encoder). Figure2 gives a description of the operating principle of an optical encoder which is our main sensor. For simplicity, we overlook the quadrature effect and depict an encoder disk as the combination of transparent and dark lines. Δ is the angular resolution of the encoder. The small circle indicates the position of light beam. The angular motion θ denotes the rotational motion of the encoder. At the start, the encoder is declared at its “zero position” and θ=0.Figure 2
Encoder schematic and functionalities.We now present a critical argument on the position detection ability of the optical encoder. It is about an encoder offset angle, which we will denote bya, 0<a<Δ. It is the angular motion needed for the first encoder pulse to register under a counterclockwise rotation (see Figure 2). “a” can be taken as a dead zone during which the encoder does not respond. Notice that this offset angle is not really an “offset position” from a baseline configuration. On the contrary, it is a measurement where we observe the first pulse. This value presents no importance at all, when the monitored motion is a few orders of magnitude larger than the resolution, Δ. However, for the particular application here the complete range of motion is composed of only a few Δ’s (e.g., 4Δ). Thus, the offset angle a plays a very critical role and we present a set of novel observations on this issue, later in the paper.Clearly,a is an unknown quantity which is random and uniformly distributed within 0<a<Δ. That is, from the starting position θ=0 the encoder does not register any reading until -a degree of rotation (in counterclockwise sense) or -a+Δ degree (in clockwise sense) is completed (see Figure 2).In Figure 3, we further depict the sensing ability of the encoder on a hypothetical oscillation θ (shown as the red line). The encoder can only register when θ reaches angular displacements of -a+Δ,-a+2Δ,…,-a+mΔ clockwise and -a,-a-Δ,…,-a-mΔ counterclockwise. Therefore, the responses to the same angular motion may vary; that is, the encoder may sense the same θ differently, depending on encoder offset angle, a.Figure 3
Oscillatory motionθ read by encoder (blue discrete lines represent the encoder signal); 0<a<Δ is the random offset angle.Of course, if the resolution of encoder is high compared to the total stroke, such an offset angle would not cause a noticeable difference in detecting the motion. We represent this behavior using a quantizer block in SIMULINK if the resolution is high. However, for low-resolution encoders (i.e., relative to the stroke), such a quantizer block cannot reveal the true response of the encoder, as depicted in Figure3. For such cases, the correct model will be introduced in Section 3.2.Some definitions are provided next, in Figure3; A1 and A2 are the upper and lower amplitudes of θ,respectively. Without loss of generality, we assume (A1>A2>0) and define the complete actual peak-to-peak stroke as A12=A1+A2. The blue line represents the rotational motion being sensed by the encoder and the encoder reading is denoted as θenc. The peak-to-peak angular stroke which is recorded by the encoder is named Aenc, which is an integer multiple of Δ (e.g., 4Δ in Figure 3). We denote the average of A12 and Aenc over a certain number of oscillations by A12¯ and Aenc¯, respectively. The bias of actual angular motion is defined as the distance of center of A12 from “zero” and it is expressed as
(2)A12-b=A1-A22.
For a certain frequency, the hypothetical actual oscillations of the pipette can be described by the bias and the amplitude of actual angular motion. Let us express θ as
(3)θ=A12-b+Asinωt,
where ω=2πf and A and f are amplitude and frequency of harmonic wave, respectively, and we have A12=2A. To summarize these definitions, we list them again as follows:(i)
resolution of position sensor =Δ;(ii)
encoder offset angle =a (deg);(iii)
upper stroke of actual angular position from rest =A1>0 (deg);(iv)
lower stroke of actual angular position from rest =A2>0 (deg);(v)
actual peak-to-peak angular stroke =A12=A1+A2 (deg);(vi)
bias of actual angular stroke =A12-b;(vii)
peak-to-peak angular stroke detected by the encoder =Aenc (deg);(viii)
average actual angular stroke over a certain number of oscillations =A12¯ (deg);(ix)
average encoder stroke over a certain number of oscillations =Aenc¯ (deg).


## 3. Review of Earlier Work: The First-Stage Adaptation Scheme

This section presents a brief review of the earlier work [12], which establishes the departure points of the present effort.


### 3.1. Control Objective and the Sensitivity Analysis

In Ros-Drill application, the desired angular trajectory requires high control sampling frequency,fs, in the microcontroller (such as 10 KHz) in order to perform a meaningful tracking. This imposes some further constraints on the limited computational capabilities of selected microcontroller. Keeping these restrictions in mind, a proportional integral and derivative (PID) control logic is adopted here.The Ros-Drill microinjector can be considered as a simple rotational mass, which is attached to a torque-generating DC servomotor (Figure1). This yields a transfer function of Ros-Drill system and the corresponding frequency response creates the magnification factor as a function of ω, Kp, Ki, and Kd:
(4)M=|G(s=jω)|=M(fd,Kp,Ki,Kd),
where ω=2π·fd and Kp, Ki, and Kd are the control gains. The objective of the control is to achieve a flat response over a given range of operating frequencies. That is, M(fd,Kp,Ki,Kd)≈1. We wish to emphasize that the phase angle between the desired harmonic input and the resulting output is not considered as a part of the performance description. Another important point to stress is that since we have limited access to the actual rotationally oscillating motion via the coarse measurements of the encoder, this flat response characteristic can be achieved, at best, by enforcing the peak-to-peak strokes of the encoder readings to be equal to Aenc=floor(2Ad/Δ)·Δ deg.It is shown in [12] that M is much more sensitive to the Kd variations rather than Kp and Ki. We compare the effects of variation of these three control gains around their nominal values (ΔK/K). This results in some percentage variation of magnification factor (ΔM/M). Among the three gains (Kd, Kp, and Ki), the most effective one is Kd, from this perspective. Thus, the tuning of the magnification factor can be achieved more directly by utilizing only Kd variations. This tuning methodology needs to be further adapted to the actual sample data system when using a low-resolution encoders. The logic steps of this tuning methodology, or the adaptive update laws of Kd, form the main contributions of this study, and they will be detailed in the following sections.


### 3.2. SIMULINK Model

In this section, we present a realistic dynamic model of the system using MATLAB-SIMULINK platform. Figure4 shows this model for the sample data PID control system. The encoder block is created to reflect the limitations of the position sensor as closely as possible to the reality. Two pieces of information are used from the quad signature encoder outputs: one is the angular step counter, and the other is the velocity of the angular motion. These data divide the angular movement space into four sections:
(5)(i)θ>0,θ˙>0(ii)θ>0,θ˙<0(iii)θ<0,θ˙<0(iv)θ<0,θ˙>0.
Notice that only one of these four cases will be active at any sampling period; thus, the summation of the four outputs will actually declare the one-step angular accumulation. The angular measurements, Ad, in zero-order-hold (Z-O-H) mode is created with sampling time Ts just to simulate the sample data procedures of the microcontroller, see Figures 5(a) and 5(b). The integration step size for the simulation routine is chosen to be one-tenth of the sampling period of microcontroller (10 μs), which is sufficient to represent the transient behavior within a sampling period with a desirable fidelity.Figure 4
SIMULINK model of the spatially and temporally discrete control system.(a) Encoder reading scheme. (b) Model of the spatially discreteθ.
(a)
(b)The amplifier gain,K, in Figure 4 is chosen by the user based on experimental knowledge such that the saturation in the D/A converter is prevented. The section enclosed by dotted lines in Figure 4 is the dynamic model for Ros-Drill microinjector itself. Using this SIMULINK tool, we first determine some starting values of the control gains Kp, Ki, and Kd using a pole-placement-based method as described in [12]. Then, we systematically vary Kd. The objective in this tuning procedure is to obtain M=1 as close as possible for the particular interval of frequencies. For a desired motion of θd=Adsin(2πfd·t), M=1 implies a peak-to-peak stroke of 2Ad. This is the continuous rotational angle and the corresponding peak-to-peak encoder recording should be Aenc=floor(2Ad/Δ)·Δ.We continue the systematic variations ofKd until the frequency response amplitude condition, M=1,is achieved for the given operating frequency. We then repeat the same operation for other fd values on a list of potential (i.e., biologically required) operating frequencies. This set of Kd gains will form a look-up table which can be used by the real-time control program on the microcontroller. Once the user identifies a preferred harmonic frequency, the program selects a Kd feedback gain from a given list and adaptively sets the new control logic. This completes the first adaptation scheme, which follows the steps below:(a)
declare the desired amplitude and frequency via GUI (Ad, fd);(b)
select the corresponding control gain (Kd) using the look-up table;(c)
observe the encoder registrations of peak-to-peak strokes,Aenc;(d)
evaluate the average ofAenc over a predetermined duration D, Aenc¯;(e)
signal the operator the direction of the errore=floor(2Ad/Δ)·Δ-Aenc¯ when it is outside a tolerance range, |e|> tolerance. Notice that Ad is typically not an integer multiple of Δ. Therefore, we deploy 2Ad/Δ to create a comparable basis with Aenc¯;(f)
manually adjust the gain,Kd, as supervisory fine tuning and repeat the steps (c)–(f) until no violation of tolerance in (e) remains.This procedure is shown to perform the task so that the peak-to-peak strokes of Ros-Drill are guaranteed to stay desirably close toAenc=floor(2Ad/Δ)·Δ [12].


### 3.3. Experimental Setup Used for All the Tests

One common setting was utilized during the different phases of progress. Figure6 shows that experimental setup, the mechanical device along with the controller box. A FAULHABER series 2342 micromotor is the actuator driving the pipette holder. An optical encoder with 1000 lines is utilized as position sensor, enabling the angular resolution at 0.09 deg (with quad-signature characteristics). It is important to clarify that this is the best option for an encoder in terms of size, resolution, and costs demanded by the application. A control box contains the components that handle the logical operations between the sensor (encoder) and the actuator (micromotor). Its main unit (the CPU and I/O device) is a Silicon Lab’s C8051F121 μ controller. This box also contains the necessary peripheral circuits such as the converter of encoder’s quadrature-signatures into the rotational pulse counter and direction determination of the rotation, optoisolators, and power output chip which feeds directly into the DC micromotor.Experimental setup: (a) microinjector. (b) Control box.
(a)
(b)We wish to give an idea to the reader about the microcontroller’s time management structure. For the aimed Ros-Drill application, 10 KHz SISO (single input-single output) control sampling rate is selected, and the C8051F121 micro-controller can easily accommodate this speed. Out of this 100μs total loop time, 71.8 μs is used for sensing, control logic evaluation, and D/A conversion; the remaining 28.9 μs is the idle time which can be devoted to other applications (such as filtering) to be added later.In order to crossvalidate the capabilities of the tracking control, we use an independent angular motion monitoring tool. It is a dSPACE 1104 DSP card which simply performs encoder decoding duties. The motion validations provided in the entire text are from this sensing channel instead of the unit which handled the same task in the loop.


## 3.1. Control Objective and the Sensitivity Analysis

In Ros-Drill application, the desired angular trajectory requires high control sampling frequency,fs, in the microcontroller (such as 10 KHz) in order to perform a meaningful tracking. This imposes some further constraints on the limited computational capabilities of selected microcontroller. Keeping these restrictions in mind, a proportional integral and derivative (PID) control logic is adopted here.The Ros-Drill microinjector can be considered as a simple rotational mass, which is attached to a torque-generating DC servomotor (Figure1). This yields a transfer function of Ros-Drill system and the corresponding frequency response creates the magnification factor as a function of ω, Kp, Ki, and Kd:
(4)M=|G(s=jω)|=M(fd,Kp,Ki,Kd),
where ω=2π·fd and Kp, Ki, and Kd are the control gains. The objective of the control is to achieve a flat response over a given range of operating frequencies. That is, M(fd,Kp,Ki,Kd)≈1. We wish to emphasize that the phase angle between the desired harmonic input and the resulting output is not considered as a part of the performance description. Another important point to stress is that since we have limited access to the actual rotationally oscillating motion via the coarse measurements of the encoder, this flat response characteristic can be achieved, at best, by enforcing the peak-to-peak strokes of the encoder readings to be equal to Aenc=floor(2Ad/Δ)·Δ deg.It is shown in [12] that M is much more sensitive to the Kd variations rather than Kp and Ki. We compare the effects of variation of these three control gains around their nominal values (ΔK/K). This results in some percentage variation of magnification factor (ΔM/M). Among the three gains (Kd, Kp, and Ki), the most effective one is Kd, from this perspective. Thus, the tuning of the magnification factor can be achieved more directly by utilizing only Kd variations. This tuning methodology needs to be further adapted to the actual sample data system when using a low-resolution encoders. The logic steps of this tuning methodology, or the adaptive update laws of Kd, form the main contributions of this study, and they will be detailed in the following sections.


## 3.2. SIMULINK Model

In this section, we present a realistic dynamic model of the system using MATLAB-SIMULINK platform. Figure4 shows this model for the sample data PID control system. The encoder block is created to reflect the limitations of the position sensor as closely as possible to the reality. Two pieces of information are used from the quad signature encoder outputs: one is the angular step counter, and the other is the velocity of the angular motion. These data divide the angular movement space into four sections:
(5)(i)θ>0,θ˙>0(ii)θ>0,θ˙<0(iii)θ<0,θ˙<0(iv)θ<0,θ˙>0.
Notice that only one of these four cases will be active at any sampling period; thus, the summation of the four outputs will actually declare the one-step angular accumulation. The angular measurements, Ad, in zero-order-hold (Z-O-H) mode is created with sampling time Ts just to simulate the sample data procedures of the microcontroller, see Figures 5(a) and 5(b). The integration step size for the simulation routine is chosen to be one-tenth of the sampling period of microcontroller (10 μs), which is sufficient to represent the transient behavior within a sampling period with a desirable fidelity.Figure 4
SIMULINK model of the spatially and temporally discrete control system.(a) Encoder reading scheme. (b) Model of the spatially discreteθ.
(a)
(b)The amplifier gain,K, in Figure 4 is chosen by the user based on experimental knowledge such that the saturation in the D/A converter is prevented. The section enclosed by dotted lines in Figure 4 is the dynamic model for Ros-Drill microinjector itself. Using this SIMULINK tool, we first determine some starting values of the control gains Kp, Ki, and Kd using a pole-placement-based method as described in [12]. Then, we systematically vary Kd. The objective in this tuning procedure is to obtain M=1 as close as possible for the particular interval of frequencies. For a desired motion of θd=Adsin(2πfd·t), M=1 implies a peak-to-peak stroke of 2Ad. This is the continuous rotational angle and the corresponding peak-to-peak encoder recording should be Aenc=floor(2Ad/Δ)·Δ.We continue the systematic variations ofKd until the frequency response amplitude condition, M=1,is achieved for the given operating frequency. We then repeat the same operation for other fd values on a list of potential (i.e., biologically required) operating frequencies. This set of Kd gains will form a look-up table which can be used by the real-time control program on the microcontroller. Once the user identifies a preferred harmonic frequency, the program selects a Kd feedback gain from a given list and adaptively sets the new control logic. This completes the first adaptation scheme, which follows the steps below:(a)
declare the desired amplitude and frequency via GUI (Ad, fd);(b)
select the corresponding control gain (Kd) using the look-up table;(c)
observe the encoder registrations of peak-to-peak strokes,Aenc;(d)
evaluate the average ofAenc over a predetermined duration D, Aenc¯;(e)
signal the operator the direction of the errore=floor(2Ad/Δ)·Δ-Aenc¯ when it is outside a tolerance range, |e|> tolerance. Notice that Ad is typically not an integer multiple of Δ. Therefore, we deploy 2Ad/Δ to create a comparable basis with Aenc¯;(f)
manually adjust the gain,Kd, as supervisory fine tuning and repeat the steps (c)–(f) until no violation of tolerance in (e) remains.This procedure is shown to perform the task so that the peak-to-peak strokes of Ros-Drill are guaranteed to stay desirably close toAenc=floor(2Ad/Δ)·Δ [12].


## 3.3. Experimental Setup Used for All the Tests

One common setting was utilized during the different phases of progress. Figure6 shows that experimental setup, the mechanical device along with the controller box. A FAULHABER series 2342 micromotor is the actuator driving the pipette holder. An optical encoder with 1000 lines is utilized as position sensor, enabling the angular resolution at 0.09 deg (with quad-signature characteristics). It is important to clarify that this is the best option for an encoder in terms of size, resolution, and costs demanded by the application. A control box contains the components that handle the logical operations between the sensor (encoder) and the actuator (micromotor). Its main unit (the CPU and I/O device) is a Silicon Lab’s C8051F121 μ controller. This box also contains the necessary peripheral circuits such as the converter of encoder’s quadrature-signatures into the rotational pulse counter and direction determination of the rotation, optoisolators, and power output chip which feeds directly into the DC micromotor.Experimental setup: (a) microinjector. (b) Control box.
(a)
(b)We wish to give an idea to the reader about the microcontroller’s time management structure. For the aimed Ros-Drill application, 10 KHz SISO (single input-single output) control sampling rate is selected, and the C8051F121 micro-controller can easily accommodate this speed. Out of this 100μs total loop time, 71.8 μs is used for sensing, control logic evaluation, and D/A conversion; the remaining 28.9 μs is the idle time which can be devoted to other applications (such as filtering) to be added later.In order to crossvalidate the capabilities of the tracking control, we use an independent angular motion monitoring tool. It is a dSPACE 1104 DSP card which simply performs encoder decoding duties. The motion validations provided in the entire text are from this sensing channel instead of the unit which handled the same task in the loop.


## 4. New Procedure: Second-Stage Adaptation Scheme

We present the main novelty of this paper in this section. The look-up-table-based control law offered in [12], as we summarized above, can handle the potential parameter uncertainties, disturbances, and the unmodeled mechanical properties of the system only via the supervisory interference from the user, see step (f) above. For a typical application such as Ros-Drill, these uncertainties originate mainly from friction-related sources, which are impossible to model accurately. We include them in our SIMULINK model (Figure 4) as some arbitrary combinations of viscous and coulomb friction effects, just to get a feel.The control objective is again to achieve a sustained peak-to-peak stroke ofAenc=floor(2Ad/Δ)·Δ as the encoder registers it (as an integer multiple of Δ). Remember that this is the only real-time measuring ability we are given. The average deviation of the Aenc (over a predetermined number of cycles) from this desired Aenc=floor(2Ad/Δ)·Δ is used as the error for tuning of Kd. One must pay attention, however, to another angle of these encoder readings. During the experiments, it is commonly observed that Aenc fluctuates typically between two successive steps (say 3Δ and 4Δ). This feature occurs due to the random offset angle a as mentioned above, but also stems from the transient regimes where the oscillations are not yet set into periodically repeated format. To discriminate the two causes of the same effect is an impossible task. Throughout this study, we consider the motion to have reached a steady regime for simplicity, therefore disqualifying the latter cause. The effects of the transient behavior are left for a future document.In this paper, we confine the execution of supervisory gain adjustments: (a) through a windowing and averaging procedure; (b) with sufficient update-and-wait period. In short, before we adaptively selectKd we allow sufficient length of oscillations to be recorded. Therefore, Kd can change only after a certain number of oscillations. In the ICSI experiment, typically entire duration of oscillations, D, lasts less than 500 ms. The tuning of Kd is expected to be completed within a small fraction of D and in the very early stages of the period D. After extensive tests using SIMULINK with additional disturbances imposed, it is observed that the angular position of the servo system (θ) reaches the steady state no later than 5 oscillations for all the operating frequencies (from 400 to 1000 Hz). Therefore, we use the average of peak-to-peak strokes within 15 cycles of oscillations, which we denote by  Aenc¯, to assess the performance of current control gains and to update them. Obviously, if  Aenc¯  is equal to Aenc=floor(2Ad/Δ)·Δ that is considered to be satisfactory for the objectives of the control.In summary, this adaptation scheme follows the steps listed below:(a,b,c)
are identical to the first adaptation scheme above;(d)
evaluate the average ofAenc over 15 cycles, Aenc¯;(e)
utilize an update law forKd adjustments if the error, e, is outside a tolerance range, that is, |e|> tolerance. If the error is within the tolerance, no Kd update is needed.(f)
repeat step (d).These steps will again guarantee the execution of the desiredAenc using an adaptive control gain update law. The details of this adaptation are given below. For the sake of simplicity, from this point onwards we will take Ad=0.2 degrees and the corresponding Aenc=4Δ=0.36 degrees.


### 4.1. The Adaptive Update Law

Let us defineKd(n) as the nth update of the derivative gain and e(n) as ensuing amplitude error which is obtained over 15 oscillatory cycles (an informed selection based on observation of experimental data), after Kd(n) is applied. That is,
(6)e(n)=4Δ-Aenc¯.
For the next 15-cycle period, we use a new feedback gain with the following update law:(7a)Kd(n)=Kd(n-1)+ΔKd(n),
where
(7b)ΔKd(n)=C·e(n-1),and C is an update constant. This update process is performed in the following sequences:(i)
start withKd(0) which is taken from a look-up table as explained in Section 3.(ii)
Determine the resultinge(0) at the end of the following 15 cycles.(iii)
Forn=1, use ΔKd(1)=C·e(0) and evaluate Kd(1) from (7a).(iv)
Again after 15 cycles of using this control gain, determinee(1).(v)
(A)
If|e(1)|≤|e(0)|, evaluate ΔKd(2)=C·e(1), Kd(2)=Kd(1)+ΔKd(2), assign Kd-ref=Kd(1) and eref=e(1), and go to (vi).(B)
If|e(1)|>|e(0)|, evaluate ΔKd(2)=0.5C·e(0), Kd(2)=Kd(0)+ΔKd(2) and determine e(2). If |e(2)|≤|e(0)|, assign Kd-ref=Kd(2) and eref=e(2), go to (vi). If |e(2)|>|e(0)|, repeat (vB) but this time with ΔKd(3)=0.5*0.5C·e(0) and continue until the absolute error falls below e(0). Assign the current Kd(n) value to Kd-ref, current e(n) value to erefand move to (vi).(vi)
After 15 cycles, determine the newe(n) and go to (vA). Use comparison of |e(n)|≤|eref|.Some nuances on this gain adaptation process are discussed later in this section over an example data set along with the role ofKd-ref and eref.The following portion is devoted to the selection of the update constantC. Uncertain friction term is modeled as some combination (in SIMULINK model of Figure 4) of viscous and coulomb frictions:
(8)Tf(t)=Tc·sgn(θ˙)+B·θ˙,
where Tc is the coulomb friction torque and B is the viscous friction coefficient. Note that different friction scenarios refer to the model in Figure 4 with the same system parameters but with combination of different Tc and B. Scenarios with various combinations of frictions are artificially created by varying Tc and B components in (8). The selection of C is performed offline via the following steps:(a)
estimate the bounds ofTc and B based on experimental studies.(b)
DeployKd(0) in SIMULINK and find the corresponding e(0) over a certain number of cycles (say 15).(c)
DetermineΔKd(1)=C·e(0), using systematically increasing C values (starting from zero), so that the resulting error e(1) becomes desirably small (e.g., 75% as we used in our tests) compared with e(0) under most adverse friction conditions. At the same time, we note that Kd should be upper-bounded as it multiplies the encoder-based angular speed and tends to saturate the actuator input (U in Figure 4). Note that the determination of C becomes critical from the concern of minimizing the number of gain updates before the error converges to zero. Using the proposed method, determination of constant C is done offline.The following portion is devoted to the convergence ofe to zero and the determination of Kd over some experimental data. We give an interpretation of the earlier assigned variable Kd-ref. “Reference Kd, Kd-ref” is the updated derivative gain which leads to smaller amplitude error |eref| than that which evolves under the previous Kd-ref. And the corresponding new error is denoted by eref. To ensure that e converges to zero, the update law of Kd always moves the gain in the direction of smaller |eref|. The update of Kd-ref is based on eref (i.e., ΔKd=C·eref). If the resultant |e|<|eref|, this update results in a new Kd-ref and eref. Otherwise, we can continue reducing this ΔKd by 50% until a new Kd-ref as described in the updating sequence earlier.This update process ofKd is shown over an example experimental case study. In this case, experiment is done on our Ros-Drill setup for the desired amplitude of 0.2 degrees and frequency of 500 Hz. The update of Kd is performed every 15 oscillations and throughout D (i.e., 500 ms). Figure 7 shows the experimental result. Following the logic above, C=0.0256 is selected for the present operating conditions. The microcontroller is programmed to monitor the average peak-to-peak value of θenc over 15 oscillations (i.e., Aenc¯). Table 1 shows the update process of Kd for this experiment. In the first round of tuning (marked as k=0 in Table 1 and Figure 7), Kd(0)=0.0036 is selected from the look-up table for operating frequency of 500 Hz. However, e(0)=1.9Δ is large (see the inset in Figure 7), meaning that the peak-to-peak strokes are much smaller than 4Δ. In the second round of tuning (k=1 in Table 1 and Figure 7), C·e(0) is used to update Kd and Kd(1)=Kd-ref+C·e(0)=0.008. Because Kd(1) results in |e(1)|>|eref|, ΔKd(1) is bisected for the third round of tuning (k=2 in Table 1 and Figure 7), Kd(2)=Kd-ref+ΔKd(1)/2. Because |e(2)|<|eref|, Kd(2) and e(2) are taken as new Kd-ref and eref, respectively. Tuning updating process terminates here. The controlled oscillation of pipette lasts for a predetermined duration and then the pipette is returned to starting “zero position” in order to prevent the wrap-around effect (on the tubing attachments). In Figure 7, at 0.98 sec, oscillation of pipette ends and begins to return to “zero position.” Figure 8 shows the flowchart of this monitoring and on-line tuning mechanism.Table 1
Update process ofKd for Figure 7.
k
K
d
(
k
)
e
(
k
)
e
ref
K
d
-ref
Δ
K
(
k
+
1
)
0
0.0036
1.9Δ
1.9Δ
0.0036
C
·
e
(
0
)
1
0.0080
−2.66Δ
1.9Δ
0.0036
Δ
K
d
(
1
)
/
2
2
0.0058
0
0
0.0058
0Figure 7
Experimental results with second adaptation scheme (k=0,1,2 refer to the stages in Table 1).Figure 8
Flowchart of monitoring and on-line tuning mechanism.The first part of the paper, which is related to the new, second-stage adaptation scheme, ends here. In what follows, we will demonstrate the influence mechanism for the actual angular motion, if the encoder reading of the stroke isAenc=4Δ just to provide a better insight to the reader.


## 4.1. The Adaptive Update Law

Let us defineKd(n) as the nth update of the derivative gain and e(n) as ensuing amplitude error which is obtained over 15 oscillatory cycles (an informed selection based on observation of experimental data), after Kd(n) is applied. That is,
(6)e(n)=4Δ-Aenc¯.
For the next 15-cycle period, we use a new feedback gain with the following update law:(7a)Kd(n)=Kd(n-1)+ΔKd(n),
where
(7b)ΔKd(n)=C·e(n-1),and C is an update constant. This update process is performed in the following sequences:(i)
start withKd(0) which is taken from a look-up table as explained in Section 3.(ii)
Determine the resultinge(0) at the end of the following 15 cycles.(iii)
Forn=1, use ΔKd(1)=C·e(0) and evaluate Kd(1) from (7a).(iv)
Again after 15 cycles of using this control gain, determinee(1).(v)
(A)
If|e(1)|≤|e(0)|, evaluate ΔKd(2)=C·e(1), Kd(2)=Kd(1)+ΔKd(2), assign Kd-ref=Kd(1) and eref=e(1), and go to (vi).(B)
If|e(1)|>|e(0)|, evaluate ΔKd(2)=0.5C·e(0), Kd(2)=Kd(0)+ΔKd(2) and determine e(2). If |e(2)|≤|e(0)|, assign Kd-ref=Kd(2) and eref=e(2), go to (vi). If |e(2)|>|e(0)|, repeat (vB) but this time with ΔKd(3)=0.5*0.5C·e(0) and continue until the absolute error falls below e(0). Assign the current Kd(n) value to Kd-ref, current e(n) value to erefand move to (vi).(vi)
After 15 cycles, determine the newe(n) and go to (vA). Use comparison of |e(n)|≤|eref|.Some nuances on this gain adaptation process are discussed later in this section over an example data set along with the role ofKd-ref and eref.The following portion is devoted to the selection of the update constantC. Uncertain friction term is modeled as some combination (in SIMULINK model of Figure 4) of viscous and coulomb frictions:
(8)Tf(t)=Tc·sgn(θ˙)+B·θ˙,
where Tc is the coulomb friction torque and B is the viscous friction coefficient. Note that different friction scenarios refer to the model in Figure 4 with the same system parameters but with combination of different Tc and B. Scenarios with various combinations of frictions are artificially created by varying Tc and B components in (8). The selection of C is performed offline via the following steps:(a)
estimate the bounds ofTc and B based on experimental studies.(b)
DeployKd(0) in SIMULINK and find the corresponding e(0) over a certain number of cycles (say 15).(c)
DetermineΔKd(1)=C·e(0), using systematically increasing C values (starting from zero), so that the resulting error e(1) becomes desirably small (e.g., 75% as we used in our tests) compared with e(0) under most adverse friction conditions. At the same time, we note that Kd should be upper-bounded as it multiplies the encoder-based angular speed and tends to saturate the actuator input (U in Figure 4). Note that the determination of C becomes critical from the concern of minimizing the number of gain updates before the error converges to zero. Using the proposed method, determination of constant C is done offline.The following portion is devoted to the convergence ofe to zero and the determination of Kd over some experimental data. We give an interpretation of the earlier assigned variable Kd-ref. “Reference Kd, Kd-ref” is the updated derivative gain which leads to smaller amplitude error |eref| than that which evolves under the previous Kd-ref. And the corresponding new error is denoted by eref. To ensure that e converges to zero, the update law of Kd always moves the gain in the direction of smaller |eref|. The update of Kd-ref is based on eref (i.e., ΔKd=C·eref). If the resultant |e|<|eref|, this update results in a new Kd-ref and eref. Otherwise, we can continue reducing this ΔKd by 50% until a new Kd-ref as described in the updating sequence earlier.This update process ofKd is shown over an example experimental case study. In this case, experiment is done on our Ros-Drill setup for the desired amplitude of 0.2 degrees and frequency of 500 Hz. The update of Kd is performed every 15 oscillations and throughout D (i.e., 500 ms). Figure 7 shows the experimental result. Following the logic above, C=0.0256 is selected for the present operating conditions. The microcontroller is programmed to monitor the average peak-to-peak value of θenc over 15 oscillations (i.e., Aenc¯). Table 1 shows the update process of Kd for this experiment. In the first round of tuning (marked as k=0 in Table 1 and Figure 7), Kd(0)=0.0036 is selected from the look-up table for operating frequency of 500 Hz. However, e(0)=1.9Δ is large (see the inset in Figure 7), meaning that the peak-to-peak strokes are much smaller than 4Δ. In the second round of tuning (k=1 in Table 1 and Figure 7), C·e(0) is used to update Kd and Kd(1)=Kd-ref+C·e(0)=0.008. Because Kd(1) results in |e(1)|>|eref|, ΔKd(1) is bisected for the third round of tuning (k=2 in Table 1 and Figure 7), Kd(2)=Kd-ref+ΔKd(1)/2. Because |e(2)|<|eref|, Kd(2) and e(2) are taken as new Kd-ref and eref, respectively. Tuning updating process terminates here. The controlled oscillation of pipette lasts for a predetermined duration and then the pipette is returned to starting “zero position” in order to prevent the wrap-around effect (on the tubing attachments). In Figure 7, at 0.98 sec, oscillation of pipette ends and begins to return to “zero position.” Figure 8 shows the flowchart of this monitoring and on-line tuning mechanism.Table 1
Update process ofKd for Figure 7.
k
K
d
(
k
)
e
(
k
)
e
ref
K
d
-ref
Δ
K
(
k
+
1
)
0
0.0036
1.9Δ
1.9Δ
0.0036
C
·
e
(
0
)
1
0.0080
−2.66Δ
1.9Δ
0.0036
Δ
K
d
(
1
)
/
2
2
0.0058
0
0
0.0058
0Figure 7
Experimental results with second adaptation scheme (k=0,1,2 refer to the stages in Table 1).Figure 8
Flowchart of monitoring and on-line tuning mechanism.The first part of the paper, which is related to the new, second-stage adaptation scheme, ends here. In what follows, we will demonstrate the influence mechanism for the actual angular motion, if the encoder reading of the stroke isAenc=4Δ just to provide a better insight to the reader.


## 5. Comprehension of Actual Angular Motion from Coarse Encoder Readings

The above control scheme is designed to perform an adaptively updating logic to assure a givenAenc peak-to-peak stroke. Let us take Aenc=4Δ for the sake of simplified arguments, without loss of generality. We query which operational conditions are satisfied if Aenc=4Δ reading is guaranteed, next. This quest results in some interesting observations which are stated in this section. Note that the derivations pertaining to Aenc=4Δ in what follows can be generalized with ease to cases when desired control performance is Aenc=2mΔ, for m=1,2,….To clarify some definitions about encoder reading, we list them as follows:(i)
a duration of oscillations =D;(ii)
absolute deviation ofA12 from its nearest odd integer multiple of Δ=ε.A very important observation is stated next. One can see from Figure3 that for the encoder to register 4-step stroke (i.e., Aenc=4Δ), the following are the necessary and sufficient conditions:(9a)-a+2Δ<A1<-a+3Δ,(9b)a+Δ<A2<a+2Δ.By summing (9a) and (9b), we obtain
(10)3Δ<A12<5Δ.
Here, we define ε as the absolute deviation of A12 from the nearest odd integer multiple of Δ and 0<ε<Δ. For 3Δ<A12<5Δ, we have(11a)either3Δ<A12<4Δ,ε=A12-3Δ(11b)or4Δ≤A12<5Δ,ε=5Δ-A12.For the simplicity of expressions, let us denote
(12)A12-b¯=Δ-2a2.
Note that the offset angle is an uniformly distributed random variable between zero and Δ. Hence, A12-b¯ is also an uniformly distributed random variable but between -Δ/2 and Δ/2.Proposition 1.
For any harmonic motion with(2m-1)Δ<A12<(2m+1)Δ, peak-to-peak encoder reading will always be Aenc=2mΔ if and only if the bias of stroke (A12-b) satisfies
(13)-ε2<A12-b-A12-b¯<ε2.
Furthermore if A12-b=A12-b¯, for such a harmonic motion the encoder reading will always be 2mΔ.Proof.
Without loss of generality, let us takem=2 and prove the sufficiency of (13) for Aenc=4Δ. By substituting (2) and (12) into (13), we obtain(14a)Δ-2a-ε<A1-A2<Δ-2a+ε.
(a) First consider the interval 3Δ<A12<4Δ. Note from (11a):
(14b)A12=A1+A2=3Δ+ε.Using (14a) and (14b) once for A1 and again for A2, we obtain(15a)-a+2Δ<A1<-a+2Δ+ε,(15b)a+Δ<A2<a+Δ+ε.Keeping in mind that 0<ε<Δ, one can see that satisfying (15a) and (15b) automatically satisfies (9a) and (9b), which are the necessary and sufficient conditions for Aenc=4Δ measurement (as mentioned at the beginning of Section 5).
(b) One should next consider the interval4Δ<A12<5Δ. Note from (11b):
(16)A12=A1+A2=5Δ-ε.
Using (14a) and (16) once for A1 and again for A2, we obtain(17a)-a+3Δ-ε<A1<-a+3Δ,(17b)a+2Δ-ε<A2<a+2Δ.With the condition 0<ε<Δ, (17a) and (17b) guarantee the fulfillment of (9a) and (9b), respectively. Again they are the necessary and sufficient conditions for Aenc=4Δ measurement. This completes the proof of sufficiency of Proposition 1.Next we handle the necessity clause of the proposition. We take into account thatAenc=4Δ and, therefore, ((9a) and (9b)) conditions hold. We wish to show that this necessitates (13). Let us again consider 3Δ<A12<4Δ first. Under this assumption, the right inequalities of ((9a) and (9b)) are satisfied automatically. The left inequality of (9a) is -a+2Δ<A1 and it can be expressed as(18a)-(A1+A2-3Δ)<A1-A2-(Δ-2a).
Similarly, the left inequality of (9b) is a+Δ<A2 and it can be rewritten as
(18b)A1-A2-(Δ-2a)<A1+A2-3Δ.Using the definitions in (2) and (12) and the constraint (14b), the combined inequalities of ((18a) and (18b)) render exactly the conditions given in (13).Let us now focus on the interval4Δ<A12<5Δ. We also claim that this condition, this time, together with the right inequalities of ((9a) and (9b)), will satisfy the left inequalities of ((9a) and (9b)). Following the similar procedure above, with the right inequality of (9a), A1<-a+3Δ, we arrive at(19a)A1-A2-(Δ-2a)<5Δ-A1-A2.
The right inequality of (9b), A2<a+2Δ, creates
(19b)-(5Δ-A1-A2)<A1-A2-(Δ-2a).Once again, using the definitions in (2) and (12) and the constraint (16), the combined inequalities of ((19a) and (19b)) render exactly the conditions given in (13). This completes the necessity clause of the proposition.Furthermore, ifA12-b=A12-b¯, (13) automatically holds for any A12 within 3Δ<A12<5Δ; therefore, the causality follows. This proof for m=2 can be extended for cases when Aenc=2mΔ, for m=1,2,….QEDThe important implications of Proposition1 can be summarized in the following logical sequences:(i)
A12-b¯ is a random number (due to the term “a”). However, “a” remains unchanged during the oscillatory period D; thus, A12-b¯ remains constant in the same period.(ii)
Bounds expressed by (13) suggest that the bias, A12-b, of each cycle during D stays within a bounded distance from A12-b¯  (i.e., ε/2) if and only if we maintain the encoder reading of Aenc=2mΔ.(iii)
In essence, if a control mechanism which is described earlier assures theAenc=2mΔ, this also brings a guarantee for the boundedness of |A12-b-A12-b¯|<ε/2, where ε is defined as (11a) and (11b). This implies that the bias A12-b has an attractive bonding to A12-b¯. It can fluctuate around this value but cannot run away from it.


## 6. Conclusions

A novel control system with low-resolution encoder for the desired harmonic trajectory is studied on a cellular microinjection technology called the Ros-Drill (rotationally oscillating drill). In the first part of paper, a novel adaptive control logic is developed to facilitate the tracking of the harmonic rotational motion under uncertainties, especially frictions. We demonstrate via dynamics simulations first, followed by experiments that the tracking of the harmonic rotational motion is achieved with desirable fidelity. In the second part, a stochastic analysis connecting the actual motions and their low-resolution sensory recordings is presented. It is observed that when the control structure guarantees a fixed peak-to-peak stroke, the bias of actual angular motion is bounded in a range and it is attracted to a certain predefined value.



---

*Source: 418068-2013-03-06.xml*