# Dynamic Takagi-Sugeno Model for the Control of Ultrasonic Motor

Authors: Shi Jingzhuo; Lv Lin; Zhang Yu
Journal: Journal of Control Science and Engineering
 (2011)

---

## Abstract

Model of ultrasonic motor is the foundation of the design of ultrasonic motor's speed and position controller. A two-input and one-output dynamic Takagi-Sugeno model of ultrasonic motor driving system is worked out using fuzzy reasoning modeling method in this paper. Many fuzzy reasoning modeling methods are sensitive to the initial values and easy to fall into local minimum, and have a large amount of calculation. In order to overcome these defects, equalized universe method is used in this paper to get clusters centers and obtain fuzzy clustering membership functions, and then, the unknown parameters of the conclusions of fuzzy rules are identified using least-square method. Different experimental data that are tested with different operational conditions are used to examine the validity of the fuzzy model. Comparison between experimental data and calculated data of the model indicates that the model can well describe the nonlinear characteristics among the frequency, amplitude of driving voltage and rotating speed. The proposed fuzzy model can be used to analyze the performance of ultrasonic motor driving system, and also can be used to design the speed and position controller of ultrasonic motor.


---

## Body

## 1. Introduction

The model of ultrasonic motor system is the foundation of analyzing ultrasonic motor’s performance and is also an important premise of designing motor controller and trying to improve the control performance. The energy conversion processes of ultrasonic motor have several stages [1, 2]. The first is the electromechanical energy conversion using the inverse piezoelectric effect of piezoelectric material on the stator. And then, the motion of particles on the surface of stator is transformed into rotary motion of rotor by the friction on the contact surface. Because of the nonlinear and time-varying operation of piezoelectric material and friction, the strong nonlinear characteristics and complex coupling relations of ultrasonic motor are brought out by the complicated energy conversion processes and nonlinearity of material [2]. Because of these reasons, the accurate theoretical model of ultrasonic motor is difficult to obtain.Modeling of ultrasonic motor can adopt different methods. The equivalent circuit of ultrasonic motor is a kind of model that can be used to analyze the basic characteristics of ultrasonic motors [3]. Because the electrical circuit is more familiar to us than the piezoelectric material and friction, the equivalent circuit was widely used in the early stage of the study on ultrasonic motor. But the equivalent circuit is a rude modeling method of the dynamic characteristics of ultrasonic motor. On the other hand, this type of model cannot be used to design motor’s controller directly. From the perspective of being suitable for control applications, the study gets more and more attention on the control modeling method of ultrasonic motor system. Usually, for the purpose of online realization and minimizing the amount of online calculation, the control model should be relatively simple and show the main aspects of nonlinear control of ultrasonic motor system. Due to the limitations of theoretical modeling, control modeling usually adopts identification method based on the tested data [4–11]. According to the different identification methods, the forms of model can be the transfer function, differential equation, neural network, and so forth.In recent years, fuzzy modeling method based on fuzzy reasoning is gradually arisen. The same as in the neural network model, fuzzy model is also based on experimental data, easy to show the nonlinear information. People’s fuzzy thinking process is simulated by fuzzy reasoning. So, fuzzy model is more easily mixed and makes use of people’s relevant experience knowledge. This is different from other methods. The complexity of fuzzy model is relatively low. So, another effective way is provided for the modeling of nonlinear complex system. The fuzzy modeling method is rarely used in the field of motor. The fuzzy method is mostly used to realize rotating speed and position control [12–14]. Phase difference is used as the controlled variable in [15], and a kind of fuzzy Hammerstein model is worked out. This model is made up of steady nonlinear part and dynamic linear transfer function. According to the experimental measurements, the relationship between steady rotating speed and driving frequency is realized by fuzzy steady model.This paper works out the appropriate fuzzy modeling method to obtain the dynamic model of ultrasonic motor system. Equalized universe method is used to get clusters centers and obtain fuzzy clustering membership functions. And then least-square-method is used to identify the unknown parameters of the conclusions of fuzzy rules. The two-input and one-output Takagi-Sugeno model of ultrasonic motor system is established, and the model can well show the nonlinear dynamic relationship among the amplitude of driving voltage, frequency, and rotating speed.


## 2. Experimental Data Used for Modeling

The block diagram of experimental system is shown in Figure1. An optoelectric encoder, E, is used to measure the rotating speed of ultrasonic motor. And the voltage amplitude controllers of phases A and B are designed to make the voltage amplitude of the two phases equal to the given value by adjusting the duty ratio of PWM signals. In the figure, Nref and Uref are the given value of rotating speed and amplitude of driving voltage, respectively.Figure 1
Structure of the experimental system for speed control.Supply voltage of the circuit is DC12V. Phase-shift PWM method is used in the circuit to control the H-bridge driver as shown in Figure2. In Figure 2, Q1~Q8 are power MOSFETs, T1 and T2 are transformers, and L1 and L2 are series compensating inductances. Figure 3 is the equivalent circuit of ultrasonic motor’s stator. Cd is the capacitance due to the piezoelectric element's dielectric properties called the “tank capacitance.” Rm,Cm, and Lm are known as the motional impedance, and they are combined in the form (Rm+1/jωCm+jωLm). Therefore, from the viewpoint of the driving circuit, the ultrasonic motor is a kind of capacitive load. To compensate this capacitive load, series compensating inductances L1 and L2 are used here to increase the driving performance.Figure 2
Structure of the driving and control circuit.Figure 3
Equivalent circuit of ultrasonic motor’s stator.PWM1~PWM8 are the PWM control signals forQ1~Q8. The traveling ultrasonic motor used here is the commercial used ultrasonic motor, Shinsei USR60. The principal specification of USR60 motor is shown in Table 1. In the control circuit, DSP is used to implement the control strategies shown in Figure 1, which is type DSP56F801. The CPLD is used to produce the phase-shift PWM control signals PWM1 ~ PWM8 for H-bridge.Table 1
Specification of USR60 motor.
ItemValueDriving voltage130 VrmsRated torque0.5 NmRated output power5.0 WRated rotating speed100 r/minMax. torque1.0 NmTemperature range−10 ~ + 55°CWeight260 gThe controllable variables that can be used for rotating speed control are the amplitude of driving voltage, frequency, and phase difference of the two-phase voltages. The theoretical analyses and related experimental studies show that if the phase difference of two-phase voltages is set to ±90°, better operation performance can be maintained and the energy conversion efficiency of ultrasonic motor is higher. In this paper, the phase difference of two-phase voltages is set to ±90°. The remainder controllable variables, amplitude and frequency of driving voltage can be adjusted to realize the aim of speed control. Therefore, the dynamic data, which contain the amplitude of driving voltage, frequency, and speed, should be measured in the experiments. A dynamic fuzzy model, which uses the amplitude of driving voltage and frequency as input and rotating speed as output, is established.Experiment process is designed as follows. The step response of rotating speed is measured by settingNref. The amplitude of driving voltage is changed by adjusting Uref. With the change of voltage amplitude, the frequency will be changed by speed closed-loop controller in order to maintain the rotating speed at Nref. Simultaneously, the dynamic change processes of Uref, frequency and rotating speed are recorded. A group of data is shown in Figure 4. The adjustment method of Uref and value of Nref are changed, respectively. The above processes are repeated. Many groups of data can be obtained for fuzzy modeling.Figure 4
Tested data of speed control.


## 3. Fuzzy Modeling Method of Ultrasonic Motor

The basic structure of fuzzy model is shown in Figure5. In fuzzy models, T-S model which can approximate to any nonlinear system is a kind of typical dynamic model. The form of T-S model is as follows. The premise section is fuzzy, and the conclusion section is an exact value. So, the defuzzification is not necessary in T-S model(1)Ri:If(x1isAi1),…,(xkisAik),…,(xrisAir),thenyi=ai0+ai1x1+⋯+aikxk+⋯+airxr.
Here,Ri represents the ith fuzzy rule, xk is the kth input variable, Aik is one of the fuzzy subsets of xk, yi is the output of the ith rule Ri, and aik is the unknown parameter of the rule’s conclusion, k=0,…,r.Figure 5
Basic structure of fuzzy model.According to the tested input and output data, structure identification and parameter identification are required in establishing the T-S model of ultrasonic motor system. Determining the appropriate model structure is the basis of parameter identification. The first step of structure identification is to confirm the input variables and the input space. This paper adopts equalized universe method to divide the fuzzy input space. First, set a clustering number which is the number of membership functions and rules. Secondly, divide the input space. Equalized universe method is used to determine the characters of Gaussian membership function, such as value of the center and width. The structure identification of T-S model is completed. The unknown parameters of conclusion section of T-S model can be determined using the least square method. After the model is established, the validation calculation should be carried on. If the error of model output is too large, the number of rules can be increased to improve the precision.The steps of this fuzzy modeling method can be summarized as follows.(1)
Determine the input variables, the numberp of input variables, and the number c of clusters.(2)
According to the analysis of tested data, determine the corresponding domain space[M-,M+] of driving voltage amplitude, frequency, and rotating speed.(3)
Determine the centersVi¯ and widths hi of Gaussian membership functions using equalized universe method
(2)Vi¯=M-+M+-M-c+1i,(i=1,2,…,c),hi=M+-M-c+1.(4)
Calculate the distancedik between the each data point xk and clustering  center Vi¯(3)dik=‖xk-vi‖,(k=1,2,…,n).(5)
Calculate the membership degree of modeling data(4)μik=1∑j=1c(dik/djk)2/(m-1).
Here, m is a constant, usually set m=2.(6)
The unknown parameters of conclusion section of T-S model are obtained by the least-square method.(7)
If the precision of model is not satisfactory after validation calculation, change the numberc of clusters and/or the number p of input variables, and go to step (2).


## 4. Dynamic Fuzzy Model for Speed Control of Ultrasonic Motor

Dynamic fuzzy model reflects the dynamic characteristics. The characteristics of ultrasonic motor system for speed control are related with time. Output variable of the model is the rotating speed valuen(k) at the current time. Input variables contain the value of driving voltage amplitude (peak-peak value) u(k-m) and frequency f(k-m) at the current and previous time, where m=0,1,…, and also may contain the value of rotating speed n(k-l) at the previous time, where l=1,2,….The 9 groups of data which are obtained through the above experiments are used for the fuzzy modeling. Among them, the 7 groups are used to model, and the other 2 groups are used as validation data. By studying the tested data, the domain space ofu, f, and n can be set as [60,300], [41.9572,45.5482] and [0,120], respectively. To simplify the modeling process, the quantification factors of variables are set to 1, and the quantified variables are defined asU, F, and N, respectively.In accordance with the above modeling steps, the process of T-S dynamic fuzzy model is an optimization process. According to the experience and tested data analysis, the value of rotating speedn(k) is related with the value of driving voltage amplitude (peak-peak value) u(k-m), frequency f(k-m) at the current and previous time, and rotating speed n(k-l) at the previous time, where m=0,1,…, and l=1,2,…. In order to find the proper fuzzy model, u(k-1),u(k),f(k-1),f(k),n(k-1) are firstly selected as input variables. The output variable is n(k). The number p of input variables is 5. The number c of rules is 5. The centers and widths of Gaussian membership functions are obtained by equalized universe method. The membership functions of U are shown in Figure 6.Figure 6
Membership functions ofU.After parameter identification, the five fuzzy rules are shown in Table2. In Table 2, the membership functions of input variables and output variables are Ai, Bi, and Ci, in the rule Ri, respectively. The premise section is fuzzy. The conclusion section is an exact value. The execution time to estimate the six parameters of fuzzy rule R1 is 0.2498 s. This time is appropriate for real-time control.Table 2
Fuzzy rules.
RiR1R2R3R4R5Premise partu(k-1)A1A2A3A4A5u(k)A1A2A3A4A5f(k-1)B1B2B3B4B5f(k)B1B2B3B4B5n(k-1)C1C2C3C4C5Conclusion parta0654.5386.7187.0−40.247.896a1−0.0298−0.1323−0.03210.0062−0.0855a20.05350.15610.0522−0.02740.1099a338.9429.295.8569.296−0.8417a4−54.12−38.35−10.23−8.2570.5271a50.78830.91240.93081.0110.9904The comparisons of model output data and tested data are shown in Figures7 and 8. Figure 7 is step response curves when the given value ofNref is 30 r/min. The model output and tested data have the same change trend, but the model output is significantly higher than the tested value. The maximum error is 3.76r/min, and appears at the initial stage of step response. In Figure 8, Nref is 90 r/min and the model error is less than Figure 7. In the rising stage of step response, the model output is slightly lower than the tested data. In steady stage, the output data is in good agreement with tested data. The maximum error is 12.92 r/min as shown in Figure 8.Figure 7
Comparison between model output and tested data (Nref=30 r/min).Figure 8
Comparison between model output and tested data (Nref=90 r/min).The precision of the model is not high enough as shown in Figures7 and 8. In order to improve the accuracy, the number p of input variables and the number c of clusters can be changed. After several attempts, compare the model output data with the tested data in different value of p and c, and choose the best structure of fuzzy model. The ultimately decision is that the input variables are u(k-3),u(k-2),u(k-1),u(k),f(k-3),f(k-2),f(k-1),f(k),n(k-2),  and n(k-1). The output variable is n(k). The number p of input variables is 10. The numberc of rules is 10. The membership functions of U are obtained by programming as shown in Figure 9. The fuzzy rules are shown in Table 3. The execution time to estimate the eleven parameters of fuzzy rule R1 is 0.2549 s. This time is appropriate for real-time control.Table 3
Final fuzzy rules.
RiR1R2R3R4R5R6R7R8R9R10Premise partu(k-3)A1A2A3A4A5A6A7A8A9A10u(k-2)A1A2A3A4A5A6A7A8A9A10u(k-1)A1A2A3A4A5A6A7A8A9A10u(k)A1A2A3A4A5A6A7A8A9A10f(k-3)B1B2B3B4B5B6B7B8B9B10f(k-2)B1B2B3B4B5B6B7B8B9B10f(k-1)B1B2B3B4B5B6B7B8B9B10f(k)B1B2B3B4B5B6B7B8B9B10n(k-2)C1C2C3C4C5C6C7C8C9C10n(k-1)C1C2C3C4C5C6C7C8C9C10Conclusion parta0180.91249145.2255.8169.2130.113.530.7859.0438.81a10.0937−0.041480.20.12590.43780.36070.14640.27380.86691.252a2−0.517−0.009459−0.7707−0.273−1.328−1.252−0.5266−1.185−2.115−5.656a30.66630.0075950.71370.13531.2661.4160.63141.5861.6177.234a4−0.23660.08016−0.14120.0179−0.38−0.5232−0.247−0.6777−0.3582−2.806a532.8937.2155.949.5688.8918.670.844422.93134.822.62a6−7.497−108.4−342.1−96.72−87.4344.22−21.45−6.138−163.2−9.541a7−20.24173.1275.553.49−70.97−129.351.14−44.23−64.11−42.66a8−9.343−130.9−92.69−12.2565.6363.41−30.8626.7491.0928.54a90.1530.4280.4712−0.5714−0.13220.03390.83060.38260.1745−0.0896a100.79010.15860.4861.491.0790.92050.16080.60790.8071.076Figure 9
Membership functions ofU.The model output data are given in Figures10 and  11 with no load. The first point is not painted in Figure 11. Compared with Figure 7, the model error is improved obviously. When load torque is 0.2 Nm, Figure 12 shows the speed response as the speed reference is 90 r/min. There is small error in the rising stage. Figures 13 and 14 show the speed response when load torque is 0.2 Nm and Urefis 300 V. The speed response of Figure 13(b) has fluctuations in steady state. Slope response is shown in Figures 15 and 16. The load torque is 0.2 Nm in Figure 16. The speed responses are in good agreement with the given value in Figures 15(b) and 16(b).Step response comparison between model output and tested data (Nref=30 r/min, no load).
(a)
The curve of drive voltage and frequency(b)
Comparison between model output and tested dataStep response comparison between model output and tested data (Nref=90 r/min, no load).
(a)
The curve of drive voltage and frequency(b)
Comparison between model output and tested dataStep response comparison between model output and tested data (Nref=90 r/min, load torque = 0.2 Nm).
(a)
The curve of drive voltage and frequency(b)
Comparison between model output and tested dataStep response comparison between model output and tested data (Nref=30 r/min, Uref=300 V, load torque = 0.2 Nm).
(a)
The curve of drive voltage and frequency(b)
Comparison between model output and tested dataStep response comparison between model output and tested data (Nref=90 r/min, Uref=300 V, load torque = 0.2 Nm).
(a)
The curve of drive voltage and frequency(b)
Comparison between model output and tested dataSlope response comparison between model output and tested data (Uref=300 V, no load).
(a)
The curve of drive voltage and frequency(b)
Comparison between model output and tested dataSlope response comparison between model output and tested data (Uref=300 V, load torque = 0.2 Nm).
(a)
The curve of drive voltage and frequency(b)
Comparison between model output and tested data


## 5. Conclusions

Based on the experimental data, the two-input and one-output dynamic model for speed control of ultrasonic motor is worked out using dynamic fuzzy modeling method. The model structure and parameters are identified using equalized universe method and the least-square method, respectively. Using the proposed modeling methods, a few fuzzy rules are needed. The structure of model is simpler. The amount of online calculation is smaller. Comparison between model output and tested data shows that the model can well simulate the nonlinear relationship among the amplitude of driving voltage, frequency, and rotating speed.This work shows that the fuzzy modeling method is simple. Besides the tested data, membership functions and fuzzy rules can be obtained without other prior knowledge. Fuzzy model has a strong ability of approximating to the nonlinear characteristics  and is suitable for ultrasonic motor which is a strong nonlinear system.



---

*Source: 219278-2011-09-22.xml*