1. Introduction
The V-shaped electrothermal microactuators (VEM) were firstly proposed by Michael J. Sinclair in 2000 [
1]. Since then, they have been studied, analyzed, and widely applied in micro-electro-mechanical systems (MEMS) technology. Major improvements and applications can be mentioned as micro-gripper [
2], linear micro-motors [
3,
4], rotary micro-motors [
5,
6], etc. To facilitate the analysis and calculation, the authors have built models and simulated by using the finite element method [
6,
7,
8]. The mechanical parameters such as temperature, stress, and force were surveyed relatively well. However, the modeling and design of controllers for the actuator is still a challenging problem that many researchers would like to study and develop.
Belonging to actuators in general, or VEM in particular, force and displacement are two important control variables. For thermal actuators, thermal forces depend on the displacement, and therefore, designers aim to control the displacement of the system [
9].
The most appropriately conventional control methods for VEM are open loop, such as the publications presented in [
10,
11,
12]. They are simple in their layout and hence very economical, stable due to their simplicity, and easier to construct. All of these methods require a sufficiently accurate mathematical model of VEM, including its parameters. Unfortunately, this requirement cannot always be satisfied and therefore, in principle, they are not applicable if there are some uncertainties contained in the models, as well as not being able to track the control to a desired reference, and if disturbances occur in the systems.
To overcome these disadvantages, a few conventional closed loop controllers have been investigated. A classical PID regulator is firstly considered due to its simplicity and capability of providing an acceptable performance [
13,
14]. Nevertheless, they must try to set limits to represent the plant utilizing transfer functions or state differential equations. This can lead to simulation results that are very different from those in the real object.
Hence, a linear matrix inequality-based robust PID control, a combination of a robust PID with a feedforward controller and an intelligent pre-filter, was designed in [
15]. Using this method, though, a sufficiently precise mathematical model is required. Although numerous studies have been completed in the past decades [
16,
17], obtaining an accurate model is still a difficult task. In this situation, the machine learning control method (MLC) seems to be a potential option. In MLC, a mathematical model of controlled devices is not always required [
18,
19].
Iterative learning control (ILC) is a special control concept of the intelligent MLC paradigm. The theory of iterative cybernetics was first introduced in 1978 by a Japanese author, Uchiyama [
20]. From single ideas, after about 15 years of research, ILC has become a potential research topic with many positive application results. ILC has been successfully applied to industrial robots [
21,
22], computer numerical control (CNC) machine tools [
23], wafer stage motion systems [
24], injection-molding machines [
25], chain conveyor systems [
26], rapid thermal processing [
27], etc.
In the ILC concept, the controller acts based on the observation of system inputs and outputs to force the changing of the system behavior from trial to trial. Thus, the word “iterative” is used there. After each trial, the refinements of the system output tracking error will be made via the observation of the system reaction in the past, as well as on its past performance, until the desired performance is reached. Hence, it is called “learning” [
28].
This paper aims to figure out whether the introduced ILC algorithm can be applied to output tracking the VEM with a given desired tracking performance and if it can promisingly overcome the disadvantages of these previous methods. The answer to it will be obtained by carrying out a few computational simulations.
In this work, the structure and working principle of the V-shaped electrothermal microactuators are introduced. Then, the model of the V-shaped beam system using the Simscape tool is presented in
Section 3. The iterative learning control model design utilizing for V-shaped thermal actuator is applied and simulation results are discussed in
Section 4. Finally, some directions and remarks are given in the conclusion.
2. Configuration and Working Principle of the V-Shaped Electrothermal Microactuator
Although the application of the ILC controller does not require any precise mathematical model of its controlled plan, a necessary so-called priori-information about its dynamic is compulsory, such as whether the controlled plan is stable or not [
21,
22,
28,
29,
30,
31,
32]. Therefore, the verifying of VEM stability before applying the ILC approach is a prerequisite. For this purpose, especially in a virtually real environment, the paper will use a Simscape tool, since this tool allows us to create models of physical component systems based on physical connections that directly integrate with block diagrams and other modeling paradigms [
33].
Around VEM, the Simscape tool is used to build a model for it based on the separation of thin beams into numerous electro–thermal elements linked together. Therefore, it is necessary to describe the physical characterization of the object.
The typical structure of the V-shaped actuator is shown in
Figure 1a. It consists of
n beam pairs (2) arrayed as a “V” letter with a slope angle of α in the
x-direction. They are connected to two fixed electrodes (called bond pads-1). The shuttle (3) is at the center and suspended by the V-shaped beam system.
The power supply for the system is in the form of a square or half cycle of sine waves. In the driving period (i.e., first half of the cycle), the voltage is applied to the fixed electrodes (1), the current transmits through the beam system (2), and then it produces heat. The beams will be heated and expand along the beams length , hence generating a displacement of the shuttle (3) in the y-direction.
In the returning period, the voltage goes to zero and the temperature of the beams decreases gradually. They will shrink and pull the shuttle (3) until it returns to the initial position. The geometric parameters of
Figure 1a are given in
Table 1.
Displacement of the beam tip, B, (called
) can be shown as a scheme in
Figure 2.
The vertical displacement (i.e., in the
y-direction) of the single beam can be calculated as:
The displacement of the actuator system
equals the vertical displacement of a single beam plus a half of the expansion of the shuttle
(assuming that the shuttle expands on both tips B). We have:
3. Simscape Model of V-Shaped Beam System
As mentioned before, Simscape is a multi-field simulation tool (electricity, magnetism, thermal, mechanical, hydraulic, pneumatic, other physical fields, etc.) in Simulink. It helps users to build object models relying on the built-in physical models. Hence, this tool allows us to create a multi-field transformation model of physical signals on the same interface as designing and real systems, without using the mathematical model of the controlled objects [
33].
Moreover, Simscape helps to develop control systems and check the input–output operation of the system. Based on the working principle, it can be illustrated by the physical process of converting energy in the V-shaped beam in a series of electric–thermal-mechanics. In this case, we use a thermal resistor to describe the electric–thermal conversion process, conductive/convective heat transfer elements to describe the heat transfer process in the beam (Ks1, Ks2…), and the process of heat convection to the air (Ka1, Ka2…), respectively. Hence, with the Simscape tool, it is possible to analyze VEM system performance effectively without having to use any strictly mathematical model of it.
By using the Simscape tool for modeling VEM in a virtually real environment, a single beam is separated into tiny elements along the beam length, with the assumption that each is uniform and their temperature distributions are the same. Thus, we can model each differential element of the V-shaped beam as shown in
Figure 3.
Accuracy of the model strongly depends on the number of segments and length
. However, the large segment-number will cause a longer simulation time. In this work, corresponding to the fixed length,
µm, of a single beam, each beam will be approximately separated into 15 elements, with the length of one segment
µm, as illustrated in
Figure 4.
The geometric parameters of the VEM and the material properties of silicon are given in
Table 1 and
Table 2.
The model of this actuator has been simulated with a square signal (amplitude:
V). A temperature response of the system is exhibited in
Figure 5. It represents relatively from this simulation result that VEM is a bounded input-bounded output stable object (BIBO), which met completely with the conclusions given in [
8,
9,
11,
13], as well as with our previous work presented in [
9]. Hence, the ILC approach now can be applied confidently to tracking control VEM.
Furthermore, simulation results show that the output response matches closely to the physical calculation. Here, the highest temperature (Tmax) is concentrated in the shuttle and is about 425 °C, and the average temperature (Tave.) at the center of the single beam is 256 °C. A simulation by Simscape is also performed with different input signals, such as DC or half-of-sine voltage at various values of driving voltage. The obtained results match well with analysis and simulation done by ANSYS.
4. Iterative Learning Controller Design for VEM, Simulations and Discussion
The block diagram and the working principles of iterative learning controllers are shown in
Figure 6 and
Figure 7 below.
The main purpose of the ILC scheme is to determine a learning function (also called the ILC algorithm [
31,
32]):
for controlled systems, denoted by
, which work in a repetitive manner during a finite time interval:
, so that the following condition of output tracking performance:
with any constant
and any suitably defined norm, will be satisfied after a finite trial number
, where the subscript
indicates the trial number by learning, and
is the tracking error belonging to the
trial.
is the desired reference for the controlled system
. Since the VEM is a controlled system that is not only stable (
Figure 5), but also operates in a repetitive mode, the ILC algorithm (3) can be applied therefore directly to output tracking that controls it.
The operation of the ILC algorithm (3) above could be shortly explained as follows [
28,
32]:
First, during the trial, the whole input , with 0 ≤ t ≤ T is sent to the controlled plant for producing its complete response .
Afterward, at the moment when the current trial is finished, the output tracking error (5) is calculated based on a particularly chosen learning algorithm to determining the input
for the next
trial, which should produce a smaller tracking error than the previous input
. In other words, a suitably chosen learning function (3) is that, which should generate continuously:
for
For a particular circumstance, that the input–output operator
:
as well as the chosen learning function (3) of an additional structure:
are linear, then the requirement (6) will be satisfied, if the following sufficient condition:
where
denotes the identity operator.
Moreover, it is remarkably here for the requirement (6) that in all trials, the same initial values have to be assigned to the controlled system .
- 3.
Finally, the learning process will be continued repetitively for , until the condition (4) is satisfied.
A meaningful version of ILC algorithms (3), which can be also applied during the whole time without having to check the required terminative condition (4) is the following uninterrupted one:
where
and
if
. In contrast to the common context of ILC algorithm (3), this uninterrupted learning algorithm allows the learning process to be faster and be continued forever without any degradation.
Now, it is time for coming back to the design of the ILC algorithm as seen in Equations (3) or (11) for VEM. Since Uchiyama proposed his concept of iterative control for the first time in [
20] based on the ILC learning algorithm in Equation (3), the ILC principle was afterwards developed very rapidly with a huge number of ILC algorithms having been published. Many of them can be found in [
21,
22,
29,
30,
31,
32,
33]. A good summary of ILC development is presented in [
34,
35,
36,
37], in which a few of the following practicably powerful linear ILC algorithms of the alternative structure of Equation (11) are given (see [
21,
22] for more examples):
or [
29,
33]:
and
which are often known as the Derivative (D-type) and Proportional Derivative (PD-type) algorithms.
One of the most general linear ILC algorithms (11) is the one in [
21] which introduced the PID-type algorithm as follows:
Of course, all learning parameters
, and
of the ILC algorithms above have to be determined accordingly to the required condition (4).
The quality of the controller depends on the learning formula, controller coefficients, and number of learning iterations. In this paper, PD-type (14) and PID-type algorithm (15) will be applied to design the iterative learning controller for output tracking control of the VEM.
In this subject, the discrete time version of PD-type (14) and of PID-type (15) will be used to design ILC controllers for continuous-time VEM. It means that the following learning function:
with
and
will be applied. These learning parameters are chosen based on a sufficiently convergent condition for the learning process presented in [
28,
29,
35,
37], under a theoretical assumption that the VEM mathematical model given in [
9] is exact.
The simulation results for the displacement of the top of the beam by using the discrete time PD-type learning procedure (16) are given in
Figure 8 and
Figure 9 for both cases without and with an additional output disturbance. It is exhibited extraordinary here in the obtained results, especially in
Figure 8a (without external disturbances) and
Figure 9a (with an additional output disturbance), which illustrate the dependence of maximal tracking error on applied learning step number, that the tracking error does not decrease monotonously if the learning number is increased, as proven in [
28,
29,
30,
31,
32,
33,
34,
35,
36,
37] by using a mathematical model of VEM. A reason for it could be that an uncertain difference between the mathematical models given in [
5,
8,
9] and the here used Simscape model of VEM has appeared. It can be deduced that these mathematical models may not be precise enough for using it to design a conventional open-loop controller for VEM.
Next,
Figure 10 and
Figure 11 illustrate the simulation results by using the following discrete time PID learning algorithm:
with the chosen learning parameters
and
for the two cases without and with output disturbances, respectively. These learning parameters are also chosen correspondingly to the required convergent condition presented in [
28,
29,
35,
37] and the mathematical model of VEM given in [
9]. Likely by using the PD-type, as exhibited in
Figure 10a and
Figure 11a for both cases without and with an additional disturbance in system output, respectively, the tracking error does not decrease monotonously by increasing the learning steps as usual. It implies, therefore here also, that the mathematical model presented in [
9] may not be precise enough for using it to design a conventional open-loop controller for VEM.
Figure 8c,
Figure 9c,
Figure 10c, and
Figure 11c illustrate the output tracking error. These simulation results confirm that the designed ILC controllers (16) and (17) have met the tracking performance as expected.
The obtained simulation results in
Figure 8b,
Figure 9b,
Figure 10b, and
Figure 11b show that the chosen learning function and their coefficients have produced the output that converges gradually to the desired reference with a minimal tracking error after a finite number of learning steps. Concretely, for a given trapezoidal reference, the output of the system is almost identical to the reference after 100 iterations with the minimal tracking norm = 0.1994 (Euclidean norm) by using the PD learning function (16) and after 62 iterations to reach the minimal tracking norm = 0.1811 with the PID learning function (17). When an external disturbance occurs additionally in the system, then an amplitude of about 10% compared to the amplitude of the set value is added. The output is relatively close to the desired result after 85 iterations with a minimal tracking norm of 0.2296 by using the PD learning function (16), and after 54 iterations to the minimal norm of 0.2138 with the PID learning function (17).
Both learning algorithms (16) and (17) produce a good disturbance resistance. If the disturbance is predetermined, the results will converge faster.
Figure 8,
Figure 9,
Figure 10 and
Figure 11 show that the PID learning algorithm (17) produces a relatively faster convergent rate than the PD learning algorithm (16). Additionally, if we reduce the coefficients
, and
, the learning process may be longer, but the error will be smaller.
After all, the simulation also shows that the Simscape model reflects the physical characterization of VEM more precisely than its mathematical model presented in [
9]. The control behavior obtained by simulating with Simscape model of VEM becomes, therefore, more accurate.
5. Conclusions
The paper proposes two types of discrete time iterative learning controllers to control the displacement of the continuous-time V-shaped electrothermal microactuator (VEM): the PD-type and the PID-type. To verify the output tracking performance of these so called “sampled data control systems” in a virtually real environment, the Simscape tool has been used for modeling the thermal transmission/emission in a V-shaped beam via separating a single beam into segments along the beam length. Hence, every segment could then be considered to be identical in both electrical and thermal properties. By using this model, the dynamic system behavior will be simulated nearly as in a real actuator system, which is to be controlled. Therefore, we can be confident that the introduced controller will act fast in the same way as in the real environment.
The obtained simulation results for both cases, without and with an external disturbance, showed that the designed controller is robust with disturbances, completely suitable for heat transfer in a thin beam of VEM, and met the required system performance as expected. Moreover, it helped the system to be more stable, faster in calculation, and produce a smaller tracking error.
Finally, it can be seen here that the introduced control method based on a modified PD-type and PID-type iterative learning could also be applied easily for controlling different MEMS devices such as hydraulic, pneumatic, or other physical models.