Open access peer-reviewed chapter - ONLINE FIRST

Topology Optimization Methods for Flexure Hinge Type Piezoelectric Actuators

Written By

Shitong Yang, Yuelong Li, Guangda Qiao, Xiaosong Zhang and Xiaohui Lu

Submitted: February 6th, 2022 Reviewed: February 28th, 2022 Published: April 15th, 2022

DOI: 10.5772/intechopen.103983

IntechOpen
Piezoelectric Actuators Edited by Tinghai Cheng

From the Edited Volume

Piezoelectric Actuators [Working Title]

Prof. Tinghai Cheng and Dr. Jianping Li

Chapter metrics overview

33 Chapter Downloads

View Full Metrics

Abstract

Piezoelectric actuators have the obvious advantages of simple and compact structure, high precision and long stroke. However, it is difficult to satisfy the various industrial requirements. Topology optimization method can be used to find the new configurations of the compliant mechanism, and different objective function and constraint conditions can be flexibly used to determine the compliant mechanism. In the research of piezoelectric actuators, due to the advantages of compact structure, no lubrication and large displacement magnification, compliant mechanism is extremely suitable to be introduced into the design of piezoelectric actuators. In recent years, topology optimization method is frequently used to design the compliant mechanism on piezoelectric actuator, and has become a research hotspot. In this chapter, the development of topology optimization method is introduced, the design and research on the compliant mechanism of piezoelectric actuator have been summarized, and the future research direction and challenges of topology optimization design for flexure hinge type piezoelectric actuators are prospected.

Keywords

  • topology optimization
  • piezoelectric actuator
  • flexure hinge
  • compliant mechanism
  • displacement amplification

1. Introduction

Topology optimization method is firstly developed by Bendsøe and Kikuchi [1] and Bendsøe [2], which is originally used for an elastic material distribution problem in 1988. Compared with sizing and shape optimization, topology optimization is a more powerful tool, which allows new holes and connections to be generated in the design domain with a prescribed amount of material. Recently, the topology optimization method has been extended to many other fields, such as materials science [3, 4], micro-machines [5, 6], precision and ultra-precision machining [7, 8], optical focusing [9, 10], and so on. With the development of computer technology, topology optimization is developing very fast, and the calculation process of topology optimization problem has been significantly simplified.

In the research of piezoelectric actuators, piezoelectric stacks are usually used in the driving field, which requires precise but minute motion [11, 12]. Flexure hinge type piezoelectric actuators are widely applied in precision machining, in-situ mechanical measurement, biomedicine and other fields because of its compact structure, long stroke and fast response speed [13, 14, 15]. The piezoelectric actuators are usually composed of piezoelectric stack, compliant mechanism and slider [16, 17, 18, 19, 20], here the compliant mechanisms are used as transmission and amplifying mechanisms. In the field of piezoelectric actuators, the modeling and design of compliant mechanisms are key issues. For example, bridge-type [21] and parallelogram-type [22] compliant mechanisms are applied to the reach of precision engineering. Sigmund [23, 24] has designed an inverting displacement amplifier, which has been used to obtain the maximum mechanical advantage. In addition, topology optimization methods have been assembled in commercial software [25], which makes the applications of topology optimization methods much easier.

There mainly three different types of topology optimization method are used to handle the design objectives and constraints of piezoelectric actuators: density-based methods, boundary variation methods, hard-kill methods, respectively [25, 26]. (1) Density-based methods, which include Solid Isotropic Material with Penalization (SIMP) and Rational Approximation of Materials Properties (RAMP) technique. (2) Boundary variation methods (level set and phase field). (3) Hard-kill methods, typically Evolutionary Structural Optimization (ESO) method and Bidirectional Evolutionary Structural Optimization (BESO) method. Many scholars have conducted in-depth research on the above three types of topology optimization method, the design theories and methods have developed rapidly. Topology optimization has becoming one of the important methods of piezoelectric actuator design.

However, topology optimization of piezoelectric actuators is a complicated problem, which contains the process of determining the connectivity, shape, and location of voids inside a given design domain. Traditional researches have given some topological structure for the design of piezoelectric actuators [27, 28, 29, 30, 31, 32, 33, 34], the not solved problem is that which kind of compliant mechanisms with the flexure hinges makes the output performances best. Yang et al. [35] developed a static topology optimization method to solve the problem, and some topology optimization methods have also been extended to design the compliant mechanism of piezoelectric actuators [36, 37, 38, 39, 40]. At present, there is little research on using topology optimization method to design flexible hinges, and its theories and methods are very scarce, which has a lot of exploration space. Therefore, it is particularly important to explore the theories and methods of topology optimization for the design of flexure hinge type piezoelectric actuator.

We hope that this chapter will provide a summary of the recent advances and novel applications of topology optimization methods, provide exposure to lesser known, yet promising, techniques, and serve as a resource for those new to the field. The presentation of each method focuses on new developments and novel applications. So in Section 2, the operating principle of piezoelectric actuators is introduced, the problems are described and the topology optimization model is established, then the simulation analysis is performed. The prototype is fabricated and the experimental system is built in Section 3. Systematic experimental test is conducted to study the actual performance of the actuator in Section 4, and the conclusion and discussion of this chapter are in Section 5.

Advertisement

2. Operating principle and topology optimization method

2.1 Operation principle

The operating principle of the piezoelectric stick–slip actuator designed by the topology optimization method is shown in Figure 1. The asymmetric sawtooth wave is used as the excitation signal. A motion cycle can be divided into three phases:

Figure 1.

Operating principle: (a) initial phase. (b) Stick phase. (c) Slip phase.

Initial phase(as shown in Figure 1a): At time t0, the input voltage is zero, so the piezoelectric stack is at its original length, the actuator is at rest, and the maximum static friction force between the indenter and slider is taken as the locking force F0, point Pis the contact point between the indenter and slider.

Stick phase(as shown in Figure 1b): From time t0 to t1, the length of piezoelectric stack extends with the increases of voltage. Under the effect of flexure hinge, the driving mechanism produces oblique deformation. The vertical positive pressure Fy1 on the slider increases, and the horizontal driving force Fx1 drives the slider to move d1 distance by static friction force, resulting in a large output displacement. When it reaches time t1, the elongation of piezoelectric stack reaches the maximum.

Slip phase(as shown in Figure 1c): Contrary to the stick stage, the piezoelectric stack shrinks rapidly with the sudden drop of the input voltage, the driving mechanism returns to its original shape, and the vertical positive force Fy2 of the indenter against the slider decreases. There is relative sliding between the stator and the slider, resulting in horizontal kinetic friction force Fx2 in the opposite direction. When the action of kinetic friction force is greater than that of inertial force, a small backward displacement d2 occurs.

After one motion cycle, the slider moves a net forward distance d, and the actuator realizes a large-stroke output by repeating the above process driven by periodic asymmetric sawtooth wave.

2.2 Problem description

In order to better describe the method, the symbols used in the structural topology optimization process are expressed in Table 1.

NameSymbolNameSymbol
Displacement vectoruDesign domain volumeV0
Input displacementuinPrescribed volume fractionVf
Output displacementuoutxVolume of elementve
Parasitic displacementuoutyTotal number of elementsN
Global stiffness matrixKPenalization powerp
Element stiffness matrixkeDisplacement magnification factorλ
The element stiffness matrix for an element with unit Young’s modulusk0Evaluation indicatorτ
Global element stiffness matrixKe0Minimum relative density valueρmin
Input stiffness matrixKinObjective functiong
Output stiffness matrixKoutConstraint functiongτ, gv
Input stiffnesskinAdjoint unit load vectorLx, Ly, Lin
Output stiffnesskoutAdjoint displacement vectoru¯x, u¯y, u¯in
Input force vectorfinElement of adjoint displacement vectoru¯xe, u¯ye, u¯ine
Input forcefinElement of displacement vectorue

Table 1.

The symbols of topology optimization method.

Bendsøe and Sigmund pointed out that the ratio between output and input displacements is an important objective function for compliant mechanism [25]. The problem is how to find the optimal driving mechanism shape that can achieve the design goal in a certain design domain, and we design the driving mechanism on the basis of the theories of Ref. [25, 40].

The structural characteristics of the stick–slip actuator are considered, the design domain of the flexure hinge mechanism is revealed in Figure 2, and the dimensions of the design domain are L1 = 12.5 mm, L2 = 7 mm, L3 = 24 mm, L4 = 3.5 mm, respectively. uoutx, uouty, and uinare the output displacement at point B in xdirection, ydirection and the input displacement at point A, respectively. finis the input force, kinand koutare the stiffness of input spring and output spring, respectively.

Figure 2.

The design domain of driving mechanism.

The magnification factor λis defined as the objective function, which is considered the relationship of output displacement uoutxand input displacement uin, while the parasitic displacement uoutyis small, here the evaluation indicator τis defined as a constraint condition to restrict the parasitic displacement uouty.

The displacement magnification factor λand evaluation indicator τare given as follows:

λ=uoutxuinE1
τ=uoutyuoutxE2

2.3 Mathematical formulation

A basic engineering goal of piezoelectric stick–slip actuators is to maximize the stroke amplification in structures and mechanisms, a serious of different objective functions can be formulated. The basic topology optimization model of piezoelectric actuator can be given as:

findρ=ρ1ρ2ρNTRnmingρs.t.gjρ0j=12mgvρ=ΩρdVVfV0ρi=0or1,i=1,2,,n.E3

where ρis the vector of design variables, gρis the objective function, gjρis the constraint functions, gvρis the volume constraint function respectively.

The optimization models described in Eq. (3) are implemented using the MATLAB programming language. Usually, the design aim of the piezoelectric actuators is to maximize mechanical advantage, so the objective functions can be chosen as magnification factor and mechanical efficiency [23, 24], which can be used for static response. Although the compliant hinge is only a part of the compliant mechanism, and its motion and configuration are relatively simple, we can regard the compliant hinge as a simple compliant mechanism, so the topology optimization idea of the compliant mechanism is also applicable to the design of the compliant hinge. Sigmund [23, 24] developed an inverting displacement amplifier, and in [25] mechanical efficiency (ME) is defined as an objective function.

The objective function of the piezoelectric stick–slip actuator topology optimization problem is to obtain largest magnification factor λ, and the evaluation indicator τis used to restrict the parasitic displacement uouty. So the topology optimization model is given as follows:

findρ=ρ1ρ2ρNTRNmingρ=uoutxuins.t.Ku=fingτρ=uoutyuoutx2τgvρ=e=1NveρeVfV00<Vf<10<ρminρe1,e=1,2,,N.E4

where Kis the global stiffness matrix, uis the global displacement vector and finis input force vector, respectively. ρis the design variables vector, ρmin(usually ρmin=0.001) is introduced to avoid singularity case. veis the volume of element e, V0and Vfare the design domain volume and the prescribed volume fraction, Nis the number of elements used to discretize the design domain, respectively.

The element stiffness matrix kecan be written as

ke=ρepk0E5

where prepresents the penalization power (typically p=3), and k0indicates the element stiffness matrix for an element with unit Young’s modulus.

The global stiffness matrix Kcan be expressed as

K=e=1NρepKe0+Kin+KoutE6

where Ke0is the element stiffness matrix for an element with unit Young’s modulus in global sense, Kinis the stiffness matrix of the input spring kinat the global level, Koutis the stiffness matrix of the output spring koutat the global level, respectively.

2.4 Sensitivity analysis

In this part, how to determinate the output displacement uoutx, parasitic displacement uoutyand input displacement uinand their derivatives of uoutxρe, uoutyρeand uinρe, is a key problem. In addition, the relationships of uoutx, uoutyand uincan be considered as.

uoutx=LxTu,uouty=LyTu,uin=LinTuE7

where Lx, Lyand Linare the vectors for which the inner product with uproduces the relevant output displacement uoutx, uoutyand the relevant input displacement uin(Lx, Lyand Linare interpreted as a (unit) load vector), respectively.

The derivatives of uoutx, uoutyand uincan be expressed as

uoutxρe=LxTuρeE8
uoutyρe=LyTuρeE9
uinρe=LinTuρeE10

Differentiating the static equation Ku=fin, we have

uρe=K1finρeKρeuE11

Due to finis a permanent load vector, we obtain

finρe=0E12

Differentiating the global stiffness matrix, we have

Kρe=pρep1Ke0=pρep1k0E13

Next, we need to solve the adjoint equilibrium equations.

Ku¯x=Lx,Ku¯y=Ly,Ku¯in=LinE14

where u¯x, u¯yand u¯inare the adjoint displacement vectors, which are determined by adjoint (unit) load vector Lx, Lyand Lin, respectively.

Combined with Eq. (8), Eqs. (11)(14), uoutxρecan be written as

uoutxρe=LxTuρe=LxTK1finρeKρeu=K1LxTfinρeKρeu=u¯xTKρeu=pρep1u¯xTKe0u=pρep1u¯xeTk0ueE15

Then we have

uoutyρe=pρep1u¯yeTk0ueE16
uinρe=pρep1u¯ineTk0ue=pρep1ueTk0uefinE17

where u¯xe, u¯ye, u¯ineand ueare the elements of displacement vectors of u¯x, u¯y, u¯inand u, respectively.

The sensitivity of objective function gis found as

gρe=uoutxρeuinuinρeuoutxuin2E18

and the sensitivity of constraint functions gτand gvcan be expressed as

gτρe=2uoutyuoutxuoutyρeuoutxuoutxρeuoutyuoutx2E19
gvρe=veE20

To avoid the checkerboards patterns and mesh dependencies phenomena, some restriction on the design must be imposed. Here a filtering technique is used to modify the sensitivity of gρeas follows:

gρê=1ρek=1NkHkek=1NkHkeρkgρkE21

where Nkis the set of elements e, distekimplies the center-to-center distance of element eand element k, which is smaller than the filter radius rmin, and the convolution operator (weight factor) Hkein Eq. (21) is defined as

Hke=max0rmindistekE22

From Eq. (7) and Eqs. (15)(17), the solutions of Eqs. (18), (19) are obtained.

Figure 3 shows the flow chart of topology optimization procedures. Initially, the design domain, boundary conditions and material are defined, the design domain is discretized into finite element meshes and solve the static problem. Then the sensitivities of the objective function and constraint conditions are obtained, design variables are updated by MMA algorithm. Here the MMA algorithm is employed by Svanberg [41, 42] as the optimizer. Finally, check the convergence of the results, if it converges, the iteration will stop; if not, go to step 3.

Figure 3.

Flow chart of topology optimization procedures.

2.5 Numerical results and structural design

Figure 2 shows that the design domain is divided into 64 × 48 square four node elements, the side length of the element is 1, and the structural dimensions and material units are dimensionless. Based on the practical engineering experiments, the parameters in topology optimization model are chosen as follows: the volume fraction Vfis set to 0.2, τis set to 104, and the input stiffness kinand output stiffness koutin Eq. (6) are chosen as 10 and 0.001, respectively. and the input force is a unit force. For more details, we refer readers to Refs. [25, 40].

Figure 4 shows that the evaluation indicator τdecrease significantly with the increase of iteration number. Therefore, the first constraint restriction in Eq. (4) markedly eliminates the effect of parasitic displacement uoutyon the output displacement uoutx. As shown in Figure 5, the displacement magnification factor λraises with the increases of iteration number. With the continuous updating of design variables, the related topology optimization results are gradually converged.

Figure 4.

The relationship between the iteration number and evaluation indicator τ.

Figure 5.

The relationship between the iteration number and displacement magnification factor λ.

Figure 6 shows the topology optimization result, which is similar to the traditional four-bar mechanism, and the positions of the four flexure hinges and the tilt angles of the beams will be the main consideration in the mechanism design. Due to the gray unit in the topology image and the inaccurate contour extraction, the details of the hinges and boundaries of the structure are modified. Figure 7 showns the main parameters of the driving mechanism, the overall dimension of the flexure hinge mechanism is 32 mm × 32 mm × 9 mm, and an indenter with radius of 2.5 mm is selected to ensure that the driving mechanism can provide sufficient contact. The thicknesses of the four special-shaped flexure hinges are 0.3 mm, 0.3 mm, 0.4 mm and 0.4 mm, respectively.

Figure 6.

The result of the topology optimization.

Figure 7.

Structural parameters of the driving mechanism. (unit: mm).

2.6 Finite element analysis

Before the experimental verification, it is necessary to ensure that the flexure hinge driving mechanism can achieve the desired displacement amplification effect, thus the finite element analysis is carried out. AL7075 is selected as the material of the driving mechanism, the Young’s modulus, density and Poisson’s ratio of this material are 7.17 × 104 MPa, 2810 kg/m3 and 0.33, respectively. Take the contact point Pat the top of the driving mechanism as the reference, and the two holes are defined as fixed constraints. The elongation of the piezoelectric stack is set to 10 μm to simulate the deformation of piezoelectric elements. The static simulation deformation is displayed in Figure 8a, the total displacement from point Pto P′is 52.415 μm, the displacement in the positive direction of y-axis direction is 10.006 μm, and the displacement along the negative direction of the x-axis direction is 51.195 μm. Then the analysis of equivalent stress is shown in Figure 8b, the maximum equivalent stress of the driving mechanism is 128.68 MPa near the point Q, which is less than the allowable stress of AL7075 to ensure the safe operation of the actuator. The static simulation results verify the feasibility and reliability of the flexure hinge driving mechanism. Further, the first-order mode of the mechanism is shown in Figure 8c. To ensure the operation stability of the actuator, the driving frequency should be lower than this value as much as possible.

Figure 8.

Simulation analysis by FEA: (a) the static simulation deformation. (b) Analysis of equivalent stress. (c) First-order mode of the mechanism.

Advertisement

3. Prototype configuration and experimental system construction

The prototype configuration of the developed actuator is shown in Figure 9. It mainly includes a slider, a stator, a preload mechanism and a base. The stator contains a piezoelectric stack, a driving mechanism, a preload screw and a shim block. The locking force is the maximum static friction force between the indenter and the slider, it can be adjusted by the preload mechanism. One end of the linear-motion guide is a slider for precise long-stroke motion, and the other end is fixed on the base. The base is used for fixing and supporting.

Figure 9.

Prototype configuration of the developed actuator.

As shown in Figure 10, an experimental system is constructed to study the motion characteristics of the actuator. The experimental system mainly includes the prototype, a waveform generator, a power amplifier, a laser displacement sensor, a PC, a pulley, and a vibration-isolation platform. The asymmetric sawtooth wave is sent out by the waveform generator (WF 1974, Negative Feedback Corporation), then adjusted by the power amplifier (HSA 4051, Negative Feedback Corporation) as the excitation signal for the piezoelectric stack (AE0505D16). The motion characteristics are detected by the laser displacement sensor (LK-H020, Keyence Corporation), and saved in the PC. The pulley and wire are used to connect the slider and the weight, and the weight is used to calibrate the locking force and apply the load.

Figure 10.

The experimental system and the prototype.

Advertisement

4. Characteristic experiments

The driving voltage is set to 100 Vp-p, and the frequency characteristics of the actuator under different locking forces are tested and shown in Figure 11. Under the locking force of 1 N, the maximum velocity is 15.25 mm/s at the frequency of 650 Hz; when the locking forces are 2 N and 3 N, the maximum velocities can be obtained at the frequency values of 700 Hz and 690 Hz, and the velocities are 12.48 mm/s and 9.67 mm/s, respectively.

Figure 11.

Frequency characteristics under different locking forces.

Thus, there are corresponding different optimal frequencies under the different locking forces conditions. As displayed in Figure 12, the voltage characteristics under different locking forces are explored at the corresponding optimal frequency. It can be seen that the minimum starting voltage of the actuator increases with the increases of the locking force, and the minimum starting voltage of the actuator is 21.5 Vp-p under the locking force of 1 N. In addition, the elongation of the piezoelectric stack increases as the voltage raises, the single step displacement of the actuator increases, thereby increasing the velocity. The motion resolution is also an important parameter of the actuator, which reflects the precise positioning characteristics of the actuator. It can be seen from Figure 13, the single step motion resolution of the actuator reaches 96 nm under the locking force of 1 N.

Figure 12.

Voltage characteristics under different locking forces.

Figure 13.

The motion resolution under the locking force of 1 N.

Maintaining the voltage is 100 Vp-p, the displacement characteristics under different locking forces are plotted, as shown in Figure 14. It is obvious that the velocity is the fastest at the locking force of 1 N, the friction resistance is small at this time, so the backward motion is minimum. The load characteristics of the actuator is emerged in Figure 15, as the locking force gradually increases, the maximum load of the actuator increases significantly. Within a certain adjustment range, the greater locking force can increase the friction driving force, which improves the load capacity of the actuator. It can be seen that the velocity decreases almost linearly with the load increases, and the maximum load mass of the actuator exceeds 330 g under 3 N locking force.

Figure 14.

Displacement characteristics under different locking forces.

Figure 15.

Load characteristics under different locking forces.

The efficiency ηis usually introduced to evaluate the output capacity of the actuator, which can be calculated by

η=PoutPin=F×vPin=mg×vPin×100%E23

where Poutrepresents the output power, Pinrepresents the input power, Fis the gravity load, and vis the output velocity, respectively. The input power can be detected by the power analyzer, and according to the load characteristics of actuator (in Figure 15), the output power Poutcan be calculated.

Figure 16 shows that the efficiency of the actuator increases first and then decreases with the increases of load. It can be seen that when the locking force is 3 N, the efficiency of the actuator reaches the highest value 0.70% under the load of 240 g.

Figure 16.

Efficiencies under different locking forces.

Advertisement

5. Conclusions

A topology optimization work of piezoelectric actuators is illustrated to improve their output performances. As shown in this chapter, topology optimization is a powerful tool for the design of piezoelectric actuators. The aim of this work is to serve as an introduction for those who are interested in the piezoelectric actuator optimization research field, and provide a reference work for scholars. It has to be pointed out that there remain a series of unsolved problems in the field of designing piezoelectric actuators using topology optimization. In addition to the potential future work mentioned in each section, the following aspects should be considered to further promote the development of this research field: (1) The static topology optimization methods provide an efficient way to find the optimal design of piezoelectric actuators. However, some manufacturing-oriented static topology optimization methods for the design of piezoelectric actuators need to receive more attention. (2) Topology optimization of piezoelectric actuators considering dynamic problems is another important topic. When the loads changed rapidly, such as “stick–slip” type loads, the dynamic response problems of piezoelectric actuators should be concerned. (3) More topology optimization codes should be released, especially for the design of piezoelectric actuators. (4) Multi-materials topology optimization problem can be applied to piezoelectric actuators, the related theory should be considered by many outstanding scholars in the further.

Advertisement

Acknowledgments

This work was financially supported by the Science and Technology Development Plan of Jilin Province (Nos. 20200201057JC and 20190201108JC), and the Technology Research Planning Project of Education Department of Jilin Province (No. JJKH20220690KJ).

The authors would like to thank K. Svanberg for providing the Matlab code of MMA optimizer.

Advertisement

Conflict of interest

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

References

  1. 1. Bendsøe MP, Kikuchi N. Generating optimal topologies in structural design using a homogenization method. Computer Methods in Applied Mechanics and Engineering. 1988;71:197-224. DOI: 10.1016/0045-7825(88)90086-2
  2. 2. Bendsøe MP. Optimal shape design as a material distribution problem. Structural optimization. 1989;1:193-202. DOI: 10.1007/BF01650949
  3. 3. Bhowmick S, Hintsala E, Stauffer D, Syed Asif SA. In-situ SEM and TEM nanomechanical study of wear and failure mechanisms. Microscopy and Microanalysis. 2018;24:1934-1935. DOI: 10.1017/S1431927618010152
  4. 4. Breguet JM, Driesen W, Kaegi F, Cimprich T. Applications of Piezo-Actuated Micro-Robots in Micro-Biology and Material Science. Harbin, China: 2007 International Conference on Mechatronics and Automation; 5–8 August 2007, IEEE; 2007. pp. 57-62
  5. 5. Wang L, Chen W, Liu J, Deng J, Liu Y. A review of recent studies on non-resonant piezoelectric actuators. Mechanical Systems and Signal Processing. 2019;133:106254. DOI: 10.1016/j.ymssp.2019.106254
  6. 6. Wang S, Rong W, Wang L, Pei Z, Sun L. Design, analysis and experimental performance of a novel stick-slip type piezoelectric rotary actuator based on variable force couple driving. Smart Materials and Structures. 2017;26:055005. DOI: 10.1088/1361-665X/aa64c3
  7. 7. Liu Y, Wang L, Gu Z, Quan Q, Deng J. Development of a two-dimensional linear piezoelectric stepping platform using longitudinal-bending hybrid actuators. IEEE Transactions on Industrial Electronics. 2019;66:3030-3040
  8. 8. Tian Y, Zhang D, Shirinzadeh B. Dynamic modelling of a flexure-based mechanism for ultra-precision grinding operation. Precision Engineering. 2001;35:554-565
  9. 9. Cherepanov V, Coenen P, Voigtländer B. A nano-positioner for scanning probe microscopy: The KoalaDrive. Review of Scientific Instruments. 2012;83:023703. DOI: 10.1063/1.3681444
  10. 10. Kratochvil BE, Dong L, Nelson BJ. Real time rigid-body visual tracking in a scanning electron microscope. The International Journal of Robotics Research. 2009;28:498-511. DOI: 10.1177/0278364908099849
  11. 11. Cheng T, He M, Li H, Lu X, Zhao H, Gao H. A novel trapezoid-type stick–slip piezoelectric linear actuator using right circular flexure hinge mechanism. IEEE Transactions on Industrial Electronics. 2017;64:5545-5552. DOI: 10.1109/TIE.2017.2677318
  12. 12. Lu X, Gao Q, Li Y, Yu Y, Zhang X, Qiao G, et al. A linear piezoelectric stick-slip actuator via triangular displacement amplification mechanism. IEEE Access. 2020;8:6515-6522. DOI: 10.1109/ACCESS.2019.2963680
  13. 13. Oubellil R, Voda A, Boudaoud M, Régnier S. Mixed stepping/scanning mode control of stick-slip SEM-integrated nano-robotic systems. Sensors and Actuators A: Physical. 2019;285:258-268. DOI: 10.1016/j.sna.2018.08.042
  14. 14. Gao X, Zhang S, Deng J, Liu Y. Development of a small two-dimensional robotic spherical joint using a bonded-type piezoelectric actuator. IEEE Transactions on Industrial Electronics. 2021;68:724-733. DOI: 10.1109/TIE.2019.2959475
  15. 15. Le T, Danh LT, Jeon JU. A proposal of a piezo rotary positioning device: Design, modeling and experiments. Smart Materials and Structures. 2019;28:115032. DOI: 10.1088/1361-665X/ab4736
  16. 16. Canfield S, Frecker M. Topology optimization of compliant mechanical amplifiers for piezoelectric actuators. Structural and Multidisciplinary Optimization. 2000;20:269-279. DOI: 10.1007/s001580050157
  17. 17. Claeyssen F, Letty RL, Barillot F, Sosnicki O. Amplified piezoelectric actuators: Static and dynamic applications. Ferroelectrics. 2007;351:3-14. DOI: 10.1080/00150190701351865
  18. 18. Guo Z, Tian Y, Zhang D, Wang T, Wu M. A novel stick-slip based linear actuator using bi-directional motion of micropositioner. Mechanical Systems and Signal Processing. 2019;128:37-49. DOI: 10.1016/j.ymssp.2019.03.025
  19. 19. Huang H, Zhao H, Fan Z, Zhang H, Ma Z, Yang Z. Analysis and experiments of a novel and compact 3-DOF precision positioning platform. Journal of Mechanical Science and Technology. 2013;27:3347-3356. DOI: 10.1007/s12206-013-0856-6
  20. 20. Hunstig M, Hemsel T, Sextro W. Stick–slip and slip–slip operation of piezoelectric inertia drives. Part I: Ideal excitation. Sensors and Actuators A: Physical. 2013;200:90-100. DOI: 10.1016/j.sna.2012.11.012
  21. 21. Li J, Zhou X, Zhao H, Shao M, Hou P, Xu X. Design and experimental performances of a piezoelectric linear actuator by means of lateral motion. Smart Materials and Structures. 2015;24:065007. DOI: 10.1088/0964-1726/24/6/065007
  22. 22. Li J, Zhou X, Zhao H, Shao M, Li N, Zhang S, et al. Development of a novel parasitic-type piezoelectric actuator. IEEE/ASME Transactions on Mechatronics. 2017;22:541-550. DOI: 10.1109/TMECH.2016.2604242
  23. 23. Sigmund O. On the design of compliant mechanisms using topology optimization. Mechanics of Structures and Machines. 1997;25:493-524. DOI: 10.1080/08905459708945415
  24. 24. Sigmund O. A 99 line topology optimization code written in Matlab. Structural and Multidisciplinary Optimization. 2001;21:120-127. DOI: 10.1007/s001580050176
  25. 25. Bendsøe MP, Sigmund O. Topology Optimization: Theory: Methods and Applications. Berlin: Springer; 2004
  26. 26. Da D. Topology Optimization Design of Heterogeneous Materials and Structures. Hoboken: Wiley; 2019
  27. 27. Li Y, Li H, Cheng T, Lu X, Zhao H, Chen P. Note: Lever-type bidirectional stick-slip piezoelectric actuator with flexible hinge. Review of Scientific Instruments. 2018;89:086101. DOI: 10.1063/1.5038640
  28. 28. Li J, Huang H, Morita T. Stepping piezoelectric actuators with large working stroke for nano-positioning systems: A review. Sensors and Actuators A: Physical. 2019;292:39-51. DOI: 10.1016/j.sna.2019.04.006
  29. 29. Hunstig M. Piezoelectric inertia motors: A critical review of history, concepts, design, applications, and perspectives. Actuators. 2017;6:7. DOI: 10.3390/act6010007
  30. 30. Iqbal S, Malik A. A review on MEMS based micro displacement amplification mechanisms. Sensors and Actuators A: Physical. 2019;300:111666. DOI: 10.1016/j.sna.2019.111666
  31. 31. Kang D, Lee MG, Gweon D. Development of compact high precision linear piezoelectric stepping positioner with nanometer accuracy and large travel range. Review of Scientific Instruments. 2007;78:075112. DOI: 10.1063/1.2756627
  32. 32. Kim JH, Kim SH, Kwak YK. Development of a piezoelectric actuator using a three-dimensional bridge-type hinge mechanism. Review of Scientific Instruments. 2003;74:2918. DOI: 10.1063/1.1569411
  33. 33. Lau GK, Du HJ, Guo NQ, Lim MK. Systematic design of displacement-amplifying mechanisms for piezoelectric stacked actuators using topology optimization. Journal of Intelligent Material Systems and Structures. 2000;11:685-695. DOI: 10.1106/17V9-DR0Q-6BMU-RKUY
  34. 34. Lee HJ, Kim HC, Kim HY, Gweon DG. Optimal design and experiment of a three-axis out-of-plane nano positioning stage using a new compact bridge-type displacement amplifier. Review of Scientific Instruments. 2013;84:115103. DOI: 10.1063/1.4827087
  35. 35. Yang S, Li Y, Xia X, Ning P, Ruan W, Zheng R, et al. A topology optimization method and experimental verification of piezoelectric stick-slip actuator with flexure hinge mechanism. Archive of Applied Mechanics. 2022;92:271-285. DOI: 10.1007/s00419-021-02055-4
  36. 36. Ling M, Cao J, Zeng M. Enhanced mathematical modeling of the displacement amplification ratio for piezoelectric compliant mechanisms. Smart Materials and Structures. 2016;25:075022. DOI: 10.1088/0964-1726/25/7/075022
  37. 37. Liu K, Tovar A. An efficient 3D topology optimization code written in Matlab. Structural and Multidisciplinary Optimization. 2014;50:1175-1196. DOI: 10.1007/s00158-014-1107-x
  38. 38. Schlinquer T, Mohand-Ousaid A, .Rakotondrabe M. Displacement Amplifier Mechanism for Piezoelectric Actuators Design Using SIMP Topology Optimization Approach. In: 2018 IEEE International Conference on Robotics and Automation (ICRA); 21–25 May 2018; Brisbane, QLD, Australia: IEEE; 2018. p. 4305–4311
  39. 39. Yang S, Xia X, Liu X, Qiao G, Zhang X, Lu X. Improving velocity of stick-slip piezoelectric actuators with optimized flexure hinges based on SIMP method. IEEE Access. 2020;8:213122-213129. DOI: 10.1109/ACCESS.2020.3039857
  40. 40. Zhang X, Zhu B. Topology Optimization of Compliant Mechanisms. Singapore: Springer; 2018
  41. 41. Svanberg K. The method of moving asymptotes—A new method for structural optimization. International Journal for Numerical Methods in Engineering. 1987;24:359-373. DOI: 10.1002/nme.1620240207
  42. 42. Svanberg K. MMA and GCMMA—Two Methods for Nonlinear Optimization [Internet]. 2007. Available from:https://people.kth.se/∼krille/mmagcmma.pdf. [Accessed: May 2, 2022]

Written By

Shitong Yang, Yuelong Li, Guangda Qiao, Xiaosong Zhang and Xiaohui Lu

Submitted: February 6th, 2022 Reviewed: February 28th, 2022 Published: April 15th, 2022