Modelling and Analysis of a New Piezoelectric Dynamic Balance Regulator

In this paper, a new piezoelectric dynamic balance regulator, which can be used in motorised spindle systems, is presented. The dynamic balancing adjustment mechanism is driven by an in-plane bending vibration from an annular piezoelectric stator excited by a high-frequency sinusoidal input voltage. This device has different construction, characteristics and operating principles than a conventional balance regulator. In this work, a dynamic model of the regulator is first developed using a detailed analytical method. Thereafter, MATLAB is employed to numerically simulate the relations between the dominant parameters and the characteristics of the regulator based on thedynamic model. Finally, experimental measurements are used to certify the validity of the dynamic model. Consequently, the mathematical model presented and analysed in this paper can be used as a tool for optimising the design of a piezoelectric dynamic balance regulator during steady state operation.


Introduction
As a novel type of piezoelectric actuator, ultrasonic motors have been studied by researchers and companies all over the World for nearly 40 years. These devices use the converse piezoelectric effect of piezoceramics and convert the ultrasonic vibration of the stator into linear or rotational motion in a rotor using the friction force. They inherently posses slow speed/high torque output, high holding torque and rapid response characteristics, all of which combine for the potential to be used as a precise and accurate positioning actuator. Because standing wave-type ultrasonic motors have the advantage of accurate positioning without feedback, a new idea is proposed in this paper, an ultrasonic piezoelectric dynamic balance regulator based on the basic principles of a standing wave-type ultrasonic motor. This regulator is the key component of an in-situ dynamic balancing system for a motorised spindle to achieve high precision balancing.
There are many different types of dynamic balance regulators that have evolved over the years in the field of rotor balancing. Van de Vegte [1], who first studied in situ dynamic balance regulators for motors in 1964, proposed that a balance block, whose radial position relative to the centre of rotation was under continuous adjustment by a worm gear mechanism, could be used to achieve in situ rotor balancing. In 1978, Van de Vegte modified the dynamic balance regulator to use a mass block driven by a motor to change the angular displacement relative to the attached rotating spindle. The electrical power for the driving motor was supplied by external brush contacts [2,3]. Bishop [4] presented a similar design in 1982. Lee et al. [5] used a wireless remote-control technique to control a dynamic balance regulator in 1987. However, there were other problems with the balancing motor in addition to the difficulty in power transmission. For example, due to the significant size and inertia of the balancing motor, the dynamic balance system, which rotated with the spindle, could limit the maximum rotating speed and cause imbalance. To solve these problems, in 1998, Zeng et al. [6,7] presented a new electromagnetic dynamic balance regulator with a small axial dimension and no moving parts in the dynamic balance system. The approach of rotating a balance mass relative to the spindle was similar to the principle of a stepping motor. Up to this day, two basic types of dynamic balance devices have been used/applied in industrial fields. One type is similar to the electromagnetic dynamic balance devices described; such devices are produced by Dittel in Germany and Lord in the United States. The other typeis the sprayed-fluids dynamic balance device used primarily in grinding machines. This type of device was developed by a German company, Hofmann, in the 1960s and is still produced by some companies.
In this paper, a new type of piezoelectric dynamic balance regulator is introduced, which can overcome some of the disadvantages of existing dynamic balance regulators. As pointed out before, this regulator based on the basic principles of a standing wave-type ultrasonic motor, and combines features such as high driving torque at low rotational speed, high holding torque without an applied electric power, extremely low noise in operation, accurate positioning without feedback, simple mechanical design and rapid response. Because the regulator utilises the in-plane bending vibration mode of the piezoelectric stators to drive the rotor, a contact interface exists on the circumferential side face. This type of structure further minimises the regulator thickness and miniaturises the regulator.Meanwhile, because the regulator is driven with single-phase AC voltage, this type of drive mode simplifies the supply power and enhances the system reliability. Above all, a particularly attractive feature of the regulator is that it possesses high displacement resolution and can achieve control precision on the order of microns, thus it can achieve very high adjustment accuracy.
The topic of this paper is to perform high-efficiency dynamic modelling and to optimise the key parameters of this regulator. Analytical, numerical and experimental methods have been employed in these investigations. For simplicity, the non-linear contact mechanics are studied under the assumption that the vibrational characteristics of the resonating structure are not affected by the contact process. This analysis process simplifies the investigations considerably and quite often results in a good description of the mechanism of motion for the regulator. For an experimental illustration of the theoretical and numerical results, a sample device was fabricated and characterised based on this mechanical analysis. The experimental results reported on here confirm the validity and applicability of the regulator structure and can be used for guidance in further optimisation of the regulator.
In Section 2 of this paper, the overall structural characteristics and the operating principles of the piezoelectric dynamic balance regulator are described in detail. In Section 3, a dynamic analysis on the stator, the rotor, the friction layer and the stator-rotor contact model is performed, which permits the energy conversion characteristics of the regulator to be obtained. In Section 4, MATLAB is employed to numerically simulate the characteristics of the regulator based on the above mathematical model, and optimisation of core parameters for the regulator is conducted according to the simulation results. Afterwards, a prototype is fabricated using the optimised parameters, and its performance is tested experimentally. Comparing the experimental and simulated load torque/rotational speed relation for the regulator, the validity of the dynamic model is proven. Finally, conclusions and future research prospects are detailed in Section 5. Figure 1 shows the structure of the piezoelectric dynamic balance regulator, which can be mounted on the end of the rotor of a motorised spindle system using an interference fit or a keyway to rotate with the rotor. Except for the counterweight block, all of the components of the regulator are axially symmetric to obtain an intrinsically balanced design. There are two piezoelectric actuating devices installed in the regulator housing. Each device consists of a rotor and a piezoelectric stator with six driving teeth. The stator is fixed on the mounting surface of the regulator. The driving teeth on the stator press against the inner circumference of the rotor, which rotates in a groove on the mounting surface ( Figure 2). The stator and the rotor are pre-tensioned in the radial direction to avoid an increase in the motor thickness during operations. Additionally, a counterweight block is installed on each rotor and rotates with the rotor. When the piezoelectric stator resonates in an in-plane bending vibration mode under high-frequency sinusoidal voltage excitation, the stator will act as a friction drive to move the rotor structure to the desired position to alter the balancing vector. The combined action of the two piezoelectric actuating devices creates a resultant balancing vector that achieves dynamic balance adjustment by controlling the positions of the two counterweight blocks. The dynamic balancing principle is shown in Figure 3. The piezoelectric stator consists of an annular elastic metal body with two annular piezoceramic discs bonded to the top and bottom. Six driving teeth are distributed equally around the outer circumference of the metal body. Each bi-directionally polarised piezoceramic disc has six uniformly polarised regions, as shown in Figure 4(a). In this figure, the black arrows represent the direction of the polarisation of each region. The two identical piezoceramic discs are bonded so that the regions in the same circumferential location on the two discs have opposite polarisations. The radial centre line of each region of the piezoceramic discsis offset by a certain angle from the midline of the nearest driving tooth on the metal body, as shown in Figure 4  Under the excitation of a single-phase sinusoidal voltage at a specific frequency, a single wavelength of a standing wave is induced in each pair of neighbouring regions on each piezoceramic disc; consequently, the composite piezoelectric stator is three wavelengths long in the circumferential direction. When three identical high-frequency sinusoidal electrical signals are applied to the two piezoceramic discs, respectively, using the metal body as one electrode and the unbonded surfaces of the two piezoceramic discs as the other electrode, each disc generates its own standing wave pattern according to its polarisation pattern. The two standing waves are in phase and can excite the composite piezoelectric stator to produce the desired in-plane bending vibration modes (3,1), (3,2) and (3,3), respectively. Because the (3,1) mode can be excited at relatively low frequency and has greater level radial and circumferential displacement components than the other two modes, it is chosen as the operating vibrational mode of the regulator.

Overall Structure and Actuation Principle of the Piezoelectric Dynamic Balance Regulator
In the in-plane bending vibration mode (3,1), each point on the outer circumference of the annular stator has both radial and circumferential motion components. The circumferential displacement is non-zero for all points, except those at extreme radial displacements, i.e., the crests and troughs of the stator vibration waves. Consequently, the six driving teeth are placed on the outer circumference of the annular stator so that when the teeth contact the rotor, their radial displacement produces a positive radial pressure on the rotor. Simultaneously, the circumferential displacement of the teeth produces relative motion in the rotor, and the rotor is propelled by the contact friction between the teeth and the rotor. Figure 5 shows the displacement vector of each tooth plotted with respect to its angular position, where the vibration period is T, the time is t and the outer circumference of the stator and the six driving teeth are spread along the abscissa. In the first half of the vibration period, the three teeth (1,3 and 5) in the rising state act to propel the output rotor in a counter-clockwise direction during contact. In the second half of the vibration period, the other three teeth (2, 4 and 6) act in the same driving role. Therefore, during continuous vibration, the six teeth act in two groups to alternately propel the rotor continuously in a counter-clockwise direction.

Kinematics Modelling of the Piezoelectric Dynamic Balance Regulator
Based on motion analysis in the previous section, it is obvious that the stator intermittently contacts the rotor, so the actuation direction of the regulator is unique. The analysis includes a significant number of non-linear and uncertain factors, and the whole drive process is very complicated. To simplify the model, the following assumptions are made [8,9]: (1) the rotor is a rigid body; (2) thefriction material is visco-elastic, its surface is smooth and the influence of surface roughness is neglected; (3) the Coulomb friction law is valid at the contact interface between the stator and rotor; and (4) the regulator runs in steady state (the starting and stopping processes of the regulator are not considered).

Motion of the Stator
The displacement behaviour for the annular stator can be derived approximately using the analytical method for the in-plane vibration of the thin elastic plate [10]. First, a cylindrical coordinate system is established, with the annular stator axis as the z-axis and the middle annular surface of the stator as the r-θ plane. The resultant radial and circumferential displacement functions, u r and u θ , of the in-plane bending vibration mode of the annular stator structure are derived as follows [11]: Young's modulus, G is the shear modulus, ρ is the material density, ω is the resonant angular frequency, n is the number of the wave or nodal diameter of the in-plane vibration, A, B, C and D are coefficients determined by the boundary conditions of the annular plate structure and J n and Y n are Bessel functions of the first and second type of the n-th order, respectively. The annular stator is considered to be a laminated annular plate with clamped boundary conditions at the inner edge and free boundary conditions at the outer edge. The boundary conditions are adopted to derive the resonant frequency of the in-plane vibration mode of the piezoelectric stator and are as follows: ( ) where a and bare the outer and inner radii of the annular stator structure, respectively, and T rr and T rθ are the radial normal stress and tangential shearing stress, respectively. Substituting the displacement function in Equation (1) into the boundary conditions in Equation (2), the resonant frequency f r for the annular stator is: where α n,m is the natural frequency constant of the annular plate, the subscript m represents the number of the nodal circle and n represents the number of the wave or nodal diameter. The ANSYS finite element software is used to analyse the corresponding modal shapes of the in-plane bending vibration of the piezoelectric stator, as shown in Figure 6. The structural dimensions of the FEM model are shown in Table 1, and the material parameters are shown in Table 2.   Figure 6, the whole circumference is outspread as a sine curve, and the displacement of the outer circumference of the stator can be given by: where u rmax and u θmax are the maximum radial and circumferential displacements, respectively.
To simplify the analysis, an orthogonal coordinate system is set up, as shown in Figure 7. A local spread diagram expanded in the circumferential direction of the stator is drawn in the coordinate system, where λ is the wavelength of the standing wave. Two points are employed to perform the dynamic analysis and are located at the root and top of the stator tooth, respectively (the two points are called the root point and the top point in this paper, respectively). The variables A 1 and A 0 are the positions of the root point with and without deformation, respectively. The variables B 1 and B 0 are the positions of the top particle with and without deformation, respectively. The variable α is the swing angle of the tooth, and h T is the length of the tooth. The deformed positionA 1 (x-u θ , y+u r ) is displaced from the initial position A(x, y) by (u θ ,u r ) when the piezoceramic stator is excited by a sinusoidal electric signal. The general expressions for the vertical deformation w r and the horizontal deformation w θ of the top pointare given as Equations (5) and (6), respectively: For small values of α, the following approximations hold true: Accordingly, Equations (5) and (6) can be simplified as follows: Additionally, the swing angle α can be given as: According to the above analysis, the vertical deflection of outer circumferential surface of the stator can be expressed as: For sinnu θ ≈ nu θ , and cosnu θ ≈ 1, Equation (12) According to Equation (9), the tangential velocity of the top point can be obtained from:

Motion of Particles on Surface of Friction Layer
In the piezoelectric dynamic balance regulator, the high-frequency intermittent friction contact between the stator and the rotor is essential to generate the driving force. Tribological processes occurring in high-frequency friction contacts determine the torque-speed characteristics, lifetime and long-term behaviour of the regulator. The appropriate choice of materials for the stator-rotor interface is important in the design of the regulator. Until now, intermittent-contact type ultrasonic motors have successfully used hard contact materials such as Al 2 O 3 . However, in this paper, a special hard composite material coating is proposed for the regulator, which is softer than the base material and has a high wear resistance and stable mechanical properties with respect to temperature and environmental changes. The composite material employs an epoxy resin as base material, MoS 2 and Al 2 O 3 as padding materials. The curing temperature is possessed at 80 °C.
If the contact mechanics are studied assuming that the vibration characteristics of the stator are not affected by the contact processes, deformation will be produced on the friction layer when the rotor is pressed against the stator. The intermittent contact condition is assumed during steady state, and thus, there is not a contact gap between the interface at the driving teeth and the friction layer. Therefore, the displacement in the y-direction at a point on the friction layer is equal to the y-direction displacement of the top point. As mentioned above, the rotor is assumed to be a rigid body. Additionally, the friction layer solidifies directly on the inner circumferential surface of the rotor. Therefore, the angular velocity of the rotor is equal to the angular velocity of the friction layer: where ω R and ω Cθ are the angular velocities of the rotor and the friction layer, respectively.

Motion of the Rotor
Because the friction layer is softer than the base material of the rotor, the contact deformation is limited to the friction layer region. Thus, the rotor can be seen as a rigid body rotating on a fixed axis, and the equation describing the rotational inertia of the rotor can be given as: where m i is the quality of the point of the rotor, r i is the rotational radius of the point m i and J R is the sum of the inertia moments of the rotor.
The momentum equation for a rigid rotor rotating on a fixed axis is given as: where the right side of the equation expresses the momentum sum of all the external forces.
The dynamic differential equation describing the rotor used in this research is written as: where M D is the driving torque from the stator and M L is the load torque.

Stator-Rotor Contact Model Analysis
The mechanics of the stator-rotor contact for the regulator are complicated due to the many parameters that must be taken into account. For the case of the rotor being regarded as a rigid body and the motion of the stator is assumed independent of the contact conditions, a visco-elastic foundation model can be employed to describe the stator-rotor contact characteristics. In the proposed regulator, the normal contact resultant force F N between the stator and rotor consists of two parts; one is the pre-load force F C0 between the stator and the rotor, and the other is generated by means of in-plane vibration waves from the piezoceramic stator. In this paper, we restrict our attention to the contact problem between the stator and the rotor. For simplicity, we assume that the in-plane vibration waves are generated by a suitable distributed force F CN acting on the stator [12]. The normal contact resultant force is therefore given as: where the units of F C0 and F CN are N/m. To further simplify the analysis, the essential assumption that the distributed contact force locally depends only on the local deformation at the contact points is used. The Kelvin-Voigt constitutive equation is therefore used to describe the stress-strain relationship of the friction material and is given by: where E C and are the Young's modulus and the viscosity coefficient of the friction material, respectively. The variables σ(t) and ε(t) are the stress and the strain of the friction layer, respectively. Additionally, ε(t) can be given as: where h C is the thickness of the friction layer. For a contact area with a width of b C , the normal contact force F CN (t) is obtained as: Substituting Equation (23) into (25), the following expression is obtained: Because the friction forces depend mainly on factors such as the coefficient of static/dynamic friction, and considering the forces and stresses at the interface between the stator and the rotor, the normal contact pressure and the relative motion of interfaces, two different types of "friction" should be distinguished [8,13].
The maximum tangential resisting force at the contact interface before sliding begins is called the maximum static friction force f S and in most cases f S is proportional to the normal force F N : where μ S is the coefficient of static friction. Once relative motion between the stator and the rotor begins, a certain tangential resistance force f d , called the dynamic friction force, arises at the contact interface. In most cases, f d is also proportional to the normal force F N : where μ d is the coefficient of dynamic friction. At the contact interface, while the linear velocity of the points on top of the driving teeth are greater than the corresponding points on the surface of the friction layer, the rotor is driven by the stator, and it can be thought that the friction force creates positive work, as shown in Figure 8(a). Conversely, the rotor is slowed by the stator and the friction force can be thought to create negative work, as shown in Figure 8(b). Thus, the tangential forces on rotor can be given as: where ν ST and ν R are the linear velocities of the top point in the x-direction and the corresponding point on surface of the friction layer, respectively. Because the piezoelectric stator has six drive teeth and makes use of only three teeth to drive the rotor, according to the Coulomb friction law, the driving torque M D can be therefore calculated as: where r p is the inner radius of rotor.

Energy Conversion Characteristics of the Regulator
A common feature of all piezoelectric actuators is their two-stage energy conversion process. In the first stage, the piezoelectric elements convert electric power energy to mechanical vibration power by the piezoelectric ceramics employed in the stator to induce standing waves in the stator at frequencies in the ultrasonic range. In the second stage, the wave energy in the stator is transferred to the rotor by means of the friction contact force between them. The energy conversion efficiency of the second stage is discussed in detail in this paper. Because the intermittent contact between the stator and the rotor is cyclical, it is possible to choose a contact cycle to analyse the energy conversion characteristics [13]. The power transferred to the rotor is: The frictional losses at the interfaces are: where L C is the length of the contact area between a stator tooth and the rotor. The energy dissipation from damping in the friction layer is: where C S is the equivalent damping coefficient of the friction material. Finally, the conversion efficiency of the interface can be given by:

Effects of Some Parameters on Regulator Characteristics
In the proposed research on the piezoelectric dynamic balance regulator, the qualitative relationships between the structural parameters and the energy conversion characteristics must be determined. Therefore, MATLAB is employed to numerically simulate characteristics of the regulator based on the above mathematical model. The model parameters are found in Table 3. When some parameters are analysed, they do not use the values from Table 1, but vary between certain limits. In the following analysis, P out is the output power and (P loss + Pdiss) as a whole is the loss power.

Effects of Dynamic Friction Coefficient on Regulator Characteristics
The coefficient of dynamic friction is one of the most important properties of the friction material, and it is necessary to analyse its influence on the output characteristics of the regulator. By using the related parameters of the regulator depicted in Table 3 and Equations (21), (29), (30), (31), (32), (33) and (34), the relationship between the coefficient of dynamic friction and the regulator characteristics can be calculated and shown in Figure 9. With an increase in the coefficient of dynamic friction, the loss power exhibits nearly linear growth, and the conversion efficiency curve remains approximately unchanged for coefficients of dynamic friction varying from 0.25 to 0.35. Therefore, a friction material possessing a coefficient off riction in the above range is the best choice for this analysis.

Effects of Young's Modulus of Friction Material on Regulator Characteristics
Young's modulus is also an important property of the friction material, and it is important to research its effects on the regulator characteristics. By using the related parameters of the regulator depicted in Table 3

Efficiency
As a result, plots of the output power, loss power and conversion efficiency versus the Young's modulus are shown in Figure 10. It can be seen that an increase in the elastic modulus causes the power output and the loss power to increase at an accelerated rate, whereas the conversional efficiency falls gradually. As a result, if a relatively value for high conversion efficiency is demanded, the Young's modulus of the friction material should be chosen to be between 1 × 10 9 and 2 × 10 9 N/m 2 .

Effects of Thickness of Friction Layer on the Regulator Characteristics
By using the related parameters of the regulator depicted in Table 3 and Equations (21), (26), (29), (30), (31), (32), (33) and (34), the relationship between the Young's modulus and the regulator characteristics can be calculated. As a result, the output power, loss power and conversion efficiency for the regulator plotted versus the thickness of the friction layer (varying between 0.2 mm and 1 mm) are shown in Figure 11. It can be seen that as the thickness increases, the output power and loss power show similar downward tendencies. The conversion efficiency exhibits significant higher values when the thickness of the friction layer varies between 0.5 mm to 0.7 mm. On the basis of these analyses, it is confirmed that the friction layer thickness should be controlled to obtain higher conversion efficiencies. Therefore, 0.6 mm is chosen to be the thickness of the friction layer in this research. Figure 11. Diagram of the effects of the thickness of the friction layer on the regulator characteristics.

Effects of Pre-Load Force on the Regulator Characteristics
By using the related parameters of the regulator depicted in Table 3 and Equations (21), (22), (29), (30), (31), (32), (33) and (34), the relationship between the pre-load force and the regulator characteristicsfor a friction layer thickness of 0.6 mm can be calculated and shown in Figure 12.
The plots indicate that the output power and loss power increase with the pre-load force in a certain range. Simultaneously, the conversion efficiency decreases gradually. High pre-load forces not only

Effects of Load Torque on Regulator Characteristics
By using the related parameters of the regulator depicted in Table 3 and Equations (21), (31), (32), (33) and (34), the relationship between the load torque and the regulator characteristics for a pre-load force of 20 N and a friction layer thickness of 0.6 mm can be calculated and shown in Figure 13.

Efficiency
It can be seen that the output power and conversion efficiency show similar upward trends with an increase in the load torque. However, the loss power decreases gradually. Therefore, an increase in the load torque will not increase the loss power, which is produced by the friction between the stator and the rotor.

Comparison between Simulated and Test Results
Through the analysis of the dynamic model of the regulator by means of MATLAB, the simulation results allow one to determine the key parameters of the regulator. Using the parameters shown in Table 2, a prototype was fabricated, and its performance was tested experimentally. To verify the analytical and numerical modal analyses of the regulator, an experimental test bench is established primarily to measure the operating frequency of the in-plane bending vibration of the piezoelectric stator, as shown in Figure 14. The test bench includes a frequency response analyzer FRA5097 to generate the excitation signal, a power amplifier HSA4051 to amplify the signal to excite the stator, and both instruments are products of the NF Company. The output voltage and current of the stator are measured simultaneously by the analyzer to obtain the frequency response of the impedance. The impedance-frequency characteristic of the piezoelectric stator measured in the experiment is shown in Figure 15.  Figure 15 shows the impedance characteristics of the regulator. It is can be seen that the actual resonant frequency of the in-plane bending vibration mode (3,1) of the stator, i.e., the operating frequency of the regulator, is approximately 67 kHz. This value is slightly higher than the computational results. After that, an experimental test bench is established to measure the load torque/rotational speed relation of the regulator, as shown in Figure 16. The test bench includes a magnetic powder brake ZX-0.3YN-24 to generate the the load torque, a torque sensor (with built-in speed sensor) JN338-A to measure the load torque and rotational speed in real time, and the two equipments are produced by the MITSUBISHI Company and BEIJING SANJING Company respectively. Figure 17 shows the experimental and simulated load torque/rotational speed relation for the regulator. The round dots are the test results, and the solid line is the simulated results which is calculated by using the related parameters of the regulator depicted in Table 3.

Test result Simulated result
It can be seen that the simulation accuracy is approximately 10% lower than the measured experimental rotational speed of the regulator. The simulated result of the stall torque is approximately 3% lower than the measured experimental stall torque. Several reasons can be given for this error between experimental and simulated values. The proposed dynamic model is an approximation of the actual regulator and a parametric model in which design and optimisation processes are possible. To achieve a balance between computational efficiency and accuracy, some of the necessary assumptions may cause modelling error. For example, a possible reason for the modelling error could arise from the fact that the rotor model does not include the rotor flexion because it was assumed to behave as a rigid body. Moreover, the in-plane bending vibration waves are assumed to be generated by a distributed force acting on the stator. This approximation also causes some degree of modelling error. However, the results from the simulations and experiments generally exhibit the same trend, and thus, the validity of the dynamic model is proven. To summarise, the mathematical model presented and analysed in this paper can be used as a tool for optimising the design of a piezoelectric dynamic balance regulator in steady state operation.

Conclusions
In this paper, an ultrasonic piezoelectric dynamic balance regulator using the basic principles of a standing wave-type ultrasonic motor is proposed. Its overall structure and the actuation principle are introduced. The FEM and analytical analysis based modelling studies of the regulator are investigated in detail. The proposed model is simple and suitable for the study on friction behavior at the contact interface of regulator, and can simulate and analyze the effects of thickness and Young's modulus of friction layer, dynamic friction coefficient of contact interface, pre-load and load torque, etc. on the performances of regulator, and can be used to optimize design of the regulator.Based on the above model, MATLAB software is employed to numerically simulate characteristics of the regulator. Through the analysis of the dynamic model of the regulator by means of MATLAB, the simulation results allow one to determine the key parameters of the regulator. Using these parameters, a self-developed regulator prototype is fabricated, and its performance is tested experimentally. The experimental results are in good agreement with the simulation results, which demonstrates the effectiveness of the proposed method. Thus it can be seen that the numerical modelling approach to this complex problem has provided valuable insight into the mechanics and dynamics of the regulator. The model presented here lays the foundation for a general framework capable of the existing ultrasonic piezoelectric dynamic balance regulator and serving as a useful design tool for optimizing the key parameters of future regulators. The research on the model of the new piezoelectric dynamic balance regulator is still in its initial stages. For future applications of the regulator, systematic research and studies will be necessary. For example, the effects of stator teeth and the starting and stopping of the motor are not considered. These effects will be further researched in the future.