1. Introduction
Modern research and manufacturing systems, environment scanning systems, and micro/nano satellites rely on high-accuracy optical and laser systems. These systems create a demand for high-accuracy motors and actuators [
1,
2]. Moreover, additional requirements, such as small mass, limited footprint area or mounting space, as well as flexible driving and mounting options, are applied for such motors and actuators [
3,
4]. Usually, electromagnetic motors and actuators are used for optical and laser systems. However, electromagnetic drives have drawbacks, such as electromagnetic interference with electronic devices that affect their application option close to sensitive electronics. Also, the scalability of electromagnetic motors and actuators is limited because of the manufacturing and technological aspects. In most cases, electromagnetic motors and actuators cannot drive the payload directly; therefore, a transmission system is needed which increases the mass of the whole electro-mechanical system [
5,
6,
7].
To address these drawbacks and limitations effectively, alternative sources of motion must be introduced and utilized.
Piezoelectric drives emerge as the most promising alternative to replace electromagnetic motion actuators and motors in high-precision motion systems [
8]. Piezoelectric motors and actuators have good scalability and operate without magnetic fields, enabling them to achieve nano or micro-level motion resolution without the need for transmission or additional support systems. Moreover, piezoelectric drives possess a self-locking capability without any power consumption and can meet strict requirements of footprints, mass, and mounting space [
9,
10,
11,
12,
13]. Furthermore, piezoelectric actuators and motors can be designed as multipurpose drives that provide several degrees of freedom and generate rotary and linear motion of the payload using a single actuator [
14,
15]. However, it must be noted that multi-degree-of-freedom piezoelectric actuators have main drawbacks, such as a more complex design compared to single-degree-of-freedom actuators [
16,
17].
Piezoelectric actuators can be classified based on their operational principles into the following categories: traveling wave, inchworm, and inertial actuators [
18]. Although traveling wave and inchworm-type actuators have certain advantages, inertial actuators are the most commonly employed method for achieving linear load motion. This preference is attributed to its simplicity of excitation and its ability to attain high motion resolutions, rapid motion speeds, and small response times [
19,
20]. While there are notable advantages of already developed inertial linear motion actuators, there continues to be a demand for piezoelectric inertial linear motion actuators that are characterized by simple design, small size and volume, flexible driving, and mounting.
Lu et al. introduced an inertial linear piezoelectric actuator based on an asymmetric flexure hinge with a triangular displacement amplification mechanism driven by a multilayer piezoelectric actuator [
21]. The operation of the actuator is based on low-frequency longitudinal displacements of the multilayer piezoelectric actuator, which are inducted by the sawtooth electric signal. The longitudinal displacements of the actuator generate bending vibrations of the whole stator because of the asymmetry of the hinges and are amplified by a mechanical amplification mechanism. The actuator operates based on the slick-stick operation principle when a sawtooth waveform excitation signal is applied. Motion is transferred to the slider that is preloaded against the actuator. The authors carried out numerical and experimental investigations and showed that the maximum speed of the slider is 20.17 mm/s under the driving voltage of 100 V
p-p while the driving frequency was 610 Hz with a preload force of 2.5 N.
Yang et al. reported on the piezoelectric inertial linear motion actuator featuring an arc-shaped driving tip that is driven by a multilayer piezoelectric actuator via flexible hinges [
22]. Longitudinal displacements of the multilayer actuator are excited using a sawtooth electric signal and are transferred to an arc-shaped driving tip; therefore, inertial stick–slip operation is obtained. The authors conducted numerical and experimental investigations and found that the actuator is able to provide 2754 μm/s under a driving voltage of 120 V while the preload force was 5N.
Ko et al. introduced a piezoelectric linear motion actuator based on a bimorph piezoelectric disc that operates based on the inertial principle [
23]. The actuator consists of a bimorph piezoelectric disc clamped in a special housing. A guidance rail with a slider is placed in the center of the disc. The actuator is excited by a sawtooth signal at the frequency equal to or close to the first out-of-plane bending mode of the bimorph disc. Displacements of the slider along the rail are obtained because of the different vibration velocities of the disc when it bends in different directions. It was shown that the actuator can provide linear speed up to 15 mm/s output force of 140 mN.
This paper represents numerical and experimental investigations of a novel inertial piezoelectric linear actuator that is based on a square cross-section-shaped rod stator. The design and operating principle of the actuator is described in
Section 2. The results of numerical modeling and simulation are presented in
Section 3. Prototype and experimental measurements are provided in
Section 4. Finally,
Section 5 concludes this work.
2. Design and Operation Principle of the Actuator
The proposed linear piezoelectric actuator consists of a square cross-section-shaped rod made of beryllium bronze, eight piezoceramic plates, a T-shaped clamping element, a guidance rail, and a slider (
Figure 1 and
Figure 2). The aforementioned rod contains an asymmetric hollowed cut and T-shaped clamping element in the central zone of the rod. The cut is made with an offset from the longitudinal axis of symmetry of the rod, resulting in the formation of a slender wall. Piezoceramic plates are glued on the four surfaces of the rod while the polarization direction is pointed away from the surface. The cylindrical guidance rail is located in the center of the slender wall, perpendicular to the surface, and aligned with the transversal symmetry axis of the rod. Rail is made from fiber carbon and is used as a guidance for the slider. Finally, a T-shaped clamping element is clamped on the PCB using bolts (
Figure 1). The design of the actuator and its stator are shown in
Figure 1,
Figure 2 and
Figure 3, respectively.
The geometrical parameters of the actuator are given in
Table 1. Values of geometrical parameters were obtained during the optimization study. A detailed description of the design optimization of the stator is given in
Section 3.
The operation of the actuator is based on the inertial stick—slip principle. The first longitudinal vibration mode of the square cross-section rod is excited using two sawtooth waveform electric signals with a phase difference of π. The vibration mode of the first longitudinal mode is induced and generates transversal deformations in the central zone of the rod. Asymmetrical longitudinal vibrations are employed to obtain bending vibrations of the slender wall and longitudinal motion of the guidance rail. In order to obtain the linear motion of the slider, the asymmetrical motion of the rail is used.
In order to amplify the amplitude of the transversal vibrations in the nodal zone of the square cross-section rod, the hollowed cut is used. The position of the cut is selected in such a way that it forms the thin wall that produces bending vibrations when longitudinal vibrations of the rod are excited. Moreover, the asymmetrical cut and T-shaped clamping ensure asymmetrical stiffness of the square rod, and as a result, bending in-plane deformations of square cross-section rods are induced, which increases displacement amplitudes as well as creates the possibility of control linear motion.
Two sawtooth electric signals with a phase difference of π and frequency close to the resonance frequency of the first longitudinal deformation mode of the square cross-section-shaped rod are used for actuator driving. Eight piezo ceramic plates are divided into three groups and ensure excitation of the first longitudinal vibration mode and bending in-plane deformations of the rod. The groups of the piezoceramic plates and the actuator excitation schematic, are shown in
Figure 4. The following excitation regime is used, i.e., the first group of piezo ceramic plates is used to excite the first longitudinal vibrations (
Figure 2 and
Figure 4) when a sawtooth electric signal is applied. The other two groups of electrodes are used to induce bending vibrations, control motion direction of the slider, and increase in-plane bending deformations of the rod (
Figure 3 and
Figure 4). These groups are driven by two sawtooth signals with the phase difference of π. Therefore, by such excitation schematics, longitudinal vibration mode is excited and ensures output force and linear motion speed while, in plane, bending deformations, which appear due asymmetrical design of the stator, are enlarged by two groups of piezo ceramic plates. The motion direction of the slider is changed via swapping of signals between groups of the electrodes used for inducing bending deformations.
3. Numerical Investigation of the Actuator
Numerical investigations were performed in order to confirm the operation principle and indicate the mechanical and electromechanical characteristics of the actuator. Also, an optimization study was performed to obtain optimal geometrical and mechanical characteristics of the actuator. The Finite Element Model (FEM) was built using the COMSOL Multiphysics software with strict respect to the actuator geometry, which is represented in
Figure 2. Moreover, the boundary conditions were established as follows: the T-shaped clamping was rigidly fixed. Electrical boundary conditions were set with respect to the schematics shown in
Figure 4. Finally, the material properties were included in the model. The C17200 beryllium bronze was used for the rod and T-shaped clamping element, the properties of NCE81 (CTS Corp, USA) hard piezo ceramic were applied to piezoceramic plates, and the carbon fiber characteristics for the guide guidance rail (
Table 2).
The first step of the numerical investigation was dedicated to finding the optimal geometric parameters of the actuator. For this purpose, a frequency domain study was made in the frequency range of 35 kHz to 75 kHz with a step of 5 Hz, while the excitation signal was set to 100 V
p-p. Optimization of the actuator was performed with the goal of maximizing displacement amplitudes of the guidance rail in the
y direction by changing the geometrical parameters of the asymmetrical cut and the length of the square cross-section-shaped rod. The amplification ratio between the output and input displacement amplitudes
was used as optimization criteria. The description of the optimization problem is as follows:
subjected to
here,
ux is the longitudinal displacement of the rod;
uy is the longitudinal displacement amplitude of the guidance rail in the
y direction;
is the length of the hollowed cut excluding the radius of endings;
R is the radius of hollowed cut endings;
W is the total length of the squared-shaped rod;
is the lowest value of asymmetrical cut length;
is the highest value of asymmetrical cut length;
is the lowest value of the radius of asymmetrical cut endings; and
is the highest value of the radius of asymmetrical cut endings. In order to carry out the linear optimization study, the values of
and
were set to 3 mm and 7 mm with an increment step of 0.5 mm.
and
were set to 0.5 mm and 1.75 mm with an increment step of 0.25 mm was used. Equation (4) shows that the total length of the rod (
W) depends on the length of piezo ceramic plates, radius
R, and length
It means that the total length of the rod will be changed at every step of calculations. The results of the study are shown in
Figure 5.
Analyzing the results, it can be found that the highest longitudinal displacement amplitude of the guide rail (
uy) was obtained, while
Lcut was equal to 3.5 mm and
R was equal to 1 mm and reached 76.4 μm. On the other hand, as
Figure 5b shows, the average longitudinal displacement of rod (
ux) amplitudes is around 23.5 μm. In addition, it can be found that the
ux values have an almost even distribution among all possible combinations of geometrical parameters. Therefore, it can be stated that optimal geometrical parameters of the actuator were obtained. However, considering that the actuator design foresees amplification of displacements of the first longitudinal vibration mode, the ratio
was calculated in the full range of
Lcut and
R parameters. The results are given in
Figure 6.
The best ratio
value of 3.45 was obtained when
Lcut was 3.5 mm, and
R was 1 mm. It means that the vibration amplitude of the guidance rail was amplified 3.45 times compared to the amplitude of the first longitudinal vibrations of the square-shaped rod. More detailed results of the actuator length (
W) calculations at the resonance frequencies of the first longitudinal vibration mode are shown in
Figure 7.
It can be seen that the resonance frequency of 59.12 kHz was obtained when the optimal values of
Lcut and
R were calculated. This frequency can be called the operating frequency of the actuator. In addition,
Figure 7b shows the lengths (
W) of the rod while different geometric parameters were applied. As it can be found, it has linear behavior while its value was directly related to
Lcut and
R and calculated according to Equation (4). Therefore, considering the plot, the length (
W) of the square-shaped rod is 29.5 mm. Finally, to fully represent the results of the calculations, a summary of the results was made (
Table 3).
The geometric parameters given in
Table 3 will be used for further numerical investigations that will focus on the confirmation and indication of the mechanical and electromechanical characteristics of the actuator.
The next step of the numerical investigation was dedicated to modal analysis of the actuator. In order to carry out this study, the geometric parameters given in
Table 3 were applied to the numerical model of the actuator. The goal of the analysis was to indicate the modal shape of the actuator and to confirm the presence of bending deformations, which are inducted by the asymmetric design. The results of the modal analysis are given in
Figure 8. As can be seen, the modal shape of the actuator confirms its operation in the first longitudinal vibration mode of the square-shaped rod. In addition, it can be found that despite longitudinal vibrations of the rod, it also has bending deformations, which confirms the assumption that the asymmetrical design of the actuator, the usage of cut and T-shaped clamping, enables this type of vibration. Finally, the modal shape shows that the thin wall operates as a mechanical displacement amplifier and allows to increase in displacement amplitudes generated by longitudinal vibrations.
Impedance and phase frequency characteristics of the actuator were analyzed as well. For this purpose, a frequency domain study for a range from 59.7 to 60 kHz with the step of 5Hz was set, while mechanical and electrical boundary conditions were set as in previous studies. The results of the calculations are shown in
Figure 9.
As can be found in
Figure 9, the resonance frequency of the actuator was 59.82 kHz. A slight mismatch between resonant values obtained by modal analysis and harmonic response studies was caused by minor differences in the mesh as well as in discreet steps used in the frequency domain study. The electromechanical effective coupling
keff of the actuator was calculated (Equation (5)) and reached a value of 0.045.
here,
keff is electromechanical effective coupling;
fa is antiresonance frequency; and
fr is resonance frequency.
The time domain study was performed to indicate the operation sequence of the actuator and to confirm the operation principle. The study lasted for one period (T) of operation frequency. Mechanical boundary conditions were set as in the previous case, while the excitation signal was set to a sawtooth-shaped signal. The operation sequence of the actuator is shown in
Figure 10 and confirms the inertial operation principle of the actuator. It can be found that longitudinal displacements of the square rod are directly related to excitation signal waveform and lead to asymmetrical displacements of the rod and cylindrical guidance rail. Also, it can be noticed that displacements of the rod, despite longitudinal deformations, contain bending in-plane displacements that are inducted because of the asymmetrical design of the actuator. Finally, it can be found that the thin wall is bent by two opposite forces that are created by longitudinal deformations of the rod.
The numerical study was established with the objective of indicating the displacements and amplification characteristics of the actuator at different amplitudes of the excitation signal. To carry out this study, the frequency range was set from 59.4 to 60 kHz with a step of 2 Hz while the amplitude of the excitation signal changes from 40 V
p-p to 200 V
p-p with an increment step of 20 V
p-p. The results of the simulation are given in
Figure 11.
It can be seen that the displacement amplitudes of the cylindrical guidance rail in the y direction depend on the amplitude of the excitation signal. The lowest displacement amplitude of the rail was obtained at 40 Vp-p and reached 30.5 μm or 762.5 nm/Vp-p. The highest rail displacement amplitude was obtained at 200 Vp-p and reached 152.9 μm or 764.5 nm/Vp-p. On the other hand, the displacement amplitudes of a square cross-section-shaped rod in the x direction when the voltage of 40 Vp-p and 200 Vp-p is applied are 9.2 μm and 46.12 μm or 230 nm/Vp-p and 231 nm/Vp-p, respectively. Moreover, it can be found that the amplification ratio of the square-shaped rod displacement is approximately 3.31 times at the whole range of excitation signal amplitudes. These results show that the actuator can provide stable and well-predictable output displacements and a wide range of excitation signal amplitudes. Stable displacement characteristics ensure the possibility of using flexible and simple excitation systems for actuator control and driving.
4. Experimental Investigation of the Actuator
Experimental investigations were carried out to validate the results of numerical investigations. The actuator prototype was made with strict respect to the geometrical parameters and material properties indicated in the section related to numerical investigations (
Figure 12).
Firstly, impedance and phase characteristics were measured in the frequency domain. Hewlett-Packard 4294A impedance analyzer (Hewlett-Packard, CA, USA) was used. The results of the measurements are given in
Figure 13.
Results showed that the resonant frequency of the first longitudinal vibration mode is obtained at 59.45 kHz. The difference between the calculated and measured values is less than 1%, while the impedance values have notable differences. In addition, the electromechanical effective coupling coefficient keff was found and equaled 0.061 (Equation (5)). Differences in keff values occur due to differences in clamping and the glue layer that were neglected in the numerical model and minor mismatches in material characteristics, manufacturing, and assembly errors. However, the results of experimental and numerical investigations are in good agreement, and further experimental investigations can be performed.
Linear motion speeds of the slider were measured when different loads and excitation signal amplitudes were applied. For this purpose, an experimental setup was assembled (
Figure 14).
The experimental setup included a computer, function generator WW5064 (Tabor Electronics, Nesher, Israel), a power amplifier PD200 × 4 (Piezo Drive, Shortland, Australia), an oscilloscope DL2000 (Yokogawa, Tokyo, Japan), a displacement sensor ILD 2300 (Micro-Epsilon, Ortenburg, Germany) and a custom-made signal junction box for control of motion direction. The results of the speed measurements are shown in
Figure 15 while
Supplementary Material Video S1 represents operation of the actuator.
As can be found in
Figure 15, the highest linear motion speed of 45.6 mm/s was obtained at a voltage of 200 V
p-p while the load was 3.8 g. The motion speed per voltage ratio was 0.228 mm/s/V
p-p. On the other hand, at the same load conditions and excitation signal amplitude of 40 V
p-p, the lowest linear speed reached 9.83 mm/s or 0.245 mm/s/V
p-p. The lowest linear motion speeds were obtained when a load of 2 g was applied, and the voltage was 200 V
p-p. The linear speed of the slider reached 41.3 mm/s or 0.206 mm/s/V
p-p, while the voltage of 40 V
p-p speed reached 7.35 mm/s or 0.184 mm/s/V
p-p. Therefore, it can be stated that an actuator can provide stable and well-predictable linear motion speed at a wide range of voltage and load values.
The next step of investigations was dedicated to measuring output force at different excitation voltages and loads. The results of the measurements are shown in
Figure 16.
The highest output force of 115.4 mN was obtained at the voltage of 200 Vp-p while the load of 3.8 g was applied. The force per voltage ratio reached 577 µN/Vp-p. The lowest output force was obtained under the same load condition while the excitation signal was set at 40 Vp-p and reached 31.6 mN or 790 µN/Vp-p. In addition, it can be found that the lowest output force characteristics were obtained at a load of 2 g. The force values of 91.3 mN or 456.5 µN/Vp-p were obtained at 200 Vp-p. When the voltage was set to 40 Vp-p, the force values reached 27.9 mN or 697.5 µN/Vp-p. Despite minor fluctuations in output force values, it can be concluded that the actuator can provide stable output forces at a wide range of excitation signal amplitudes and different loads.