Abstract

Comprehensive analysis on the modal characteristics of V-shaped electrothermal microactuators is presented in this paper for the first time. Considering the unique geometric characteristics of the V-shaped beam, that is, two inclined beams supporting a movable shuttle, both the lateral and longitudinal deflections are taken into account in the modal analysis. Boundary and continuity conditions are employed to obtain the frequency equation. Natural frequencies are then obtained by solving the frequency equation. Mode shapes corresponding to their natural frequencies are also calculated analytically. The theoretical modal analysis is verified with the finite element analysis using ANSYS software. Based on the model analysis, this paper further investigates the relationship between natural frequencies and the volume scaling of the V-shaped beam. Finally, comprehensive parametric studies in terms of material properties and structural dimensions are conducted to provide insights and guidance in designing the V-shaped beam electrothermal microactuators.

1. Introduction

Electrothermal microactuators have been developed for a variety of applications in MEMS (microelectromechanical systems) including RF switches [1], micropumps [2], microgrippers [3, 4], nanopositioners [5, 6], and microtesting devices [7]. Compared to other types of actuation mechanisms such as electrostatic [8], electromagnetic [9], and piezoelectric [10], electrothermal microactuators work on the thermal expansions of beam structures and have been demonstrated to be compact, stable, and large displacement and force techniques [11, 12]. Various electrothermal actuators have been demonstrated up to date in achieving in-plane [11, 13] and out-of-plane [1416] motion.

Restricted by the current bulk micromachining process, it is difficult, if not impossible, to fabricate an out-of-plane microactuator with three-dimensional features. Therefore, most out-of-plane electrothermal actuators are designed based on existing in-plane actuator structures [17]. In fact, in-plane actuators are most extensively used in MEMS devices compared to the out-of-plane actuators. The U- [18, 19], V- [20, 21], and Z-shaped [2224] beams are the most fundamental actuator types in achieving the in-plane motion. A variety of structures that have been proposed are mainly based on these three types of beams.

The U-shaped beam electrothermal microactuators are capable of producing an arc circular motion by working on the difference of the thermal expansions between the thin and the thick beams. The U-shaped beam actuators are usually used as the microgripper [18, 25] or in MEMS positioning systems [26]. Because the thin beam has to overcome the thick beam to produce a motion and force, the U-shaped beam outputs smaller displacement and force compared to the V- and Z-shaped beams [27]. In addition, the U-shaped actuators inevitably undergo an unwanted overshoot in reaching the steady state [27]. The V- and Z-shaped beam actuators, on the contrary, comprise several pairs of symmetric V-/Z-shaped beams supporting a middle shuttle and being fixed at two ends. When applying an electrical potential difference on the two ends, the beams will expand thus pushing the shuttle forward and producing motion and force. Since those pairs of beams are symmetrically configured, the actuators with V- or Z-shaped beams can generate rectilinear motion, and there is no overshoot during operations. This rectilinear motion facilitates the development of more complex devices by combining the V- or Z-shaped actuators with compliant mechanisms in achieving more flexible functions or higher performances [28, 29]. Therefore, in recent decades, these two types of beams have attracted extensive research efforts. Although the V- and Z-shaped beams are demonstrated to behave similarly in electrothermal actuation, the stiffness of the Z-shaped beam is much smaller compared to that of the V-shaped beam [30, 31]. As a result, the Z-shaped beams are often used for bidirectional actuations [32, 33].

A number of research efforts have been devoted in static modelling [34], design [35, 36], and fabrication [37] of the V-shaped beams. The V-shaped beams are demonstrated to be capable of producing large displacement and force. To the best knowledge of the authors, few research efforts have been made in dynamics of the V-shaped beams. In the our previous preliminary work [27, 38], they have established a dynamic model of the V-shaped actuators in both air and vacuum environments to describe the dynamic behavior of the V-shaped beams. However, this model takes into account only the electrothermal response, and the thermomechanical part is treated as pseudostatic with the assumption that the mechanical frequency is much higher compared to the electrothermal response. This model is problematic when the mechanical frequency decreases with increased beam length. In this case, the mechanical vibration cannot be omitted if more accurate dynamic control is required to be achieved. Research efforts have been made to partially address this problem. For instance, Burnie [39] has calculated the first-order natural frequencies of the V-shaped beam by treating the V-shaped beam as a spring. However, this oversimplification results in 7.2 percent derivations compared to simulated results. In addition, this does not allow one to compute higher-order frequencies and the related mode shapes. Stokey [40] has studied the free vibration of a pair of simply supported beams carrying a mass in the middle. This research work provides a guidance for describing the vibration of the shuttle of the V-shaped actuators. However, the beam is not inclined with an angle, and only lateral vibrations are considered based on the traditional line-shaped beam assumptions. Our preliminary work [41], a conference paper, pointed out that the V-shaped beam involves vibrations not only in the lateral direction but also in the longitudinal direction due to the existence of the inclined angle. Preliminary results showed that the errors of natural frequencies between the analytical and the simulated results are well below 2% by incorporating both the longitudinal and lateral vibrations. This demonstrates that the vibration of the V-shaped beam is essentially coupled between longitudinal and lateral directions. Furthermore, the conference paper only presented the calculation of natural frequencies. The computation of the mode shapes were only conducted numerically using finite element simulations.

A few research efforts have also reported parameter studies of the V-shaped electrothermal microactuators based on the analytical models. Zhu et al. [7, 22] and Shen and Chen [42] discussed the effect of the inclined angle and beam width on the static deflection, stiffness, and axial internal force. In our previous papers [30, 31], we provide comprehensive parameter studies on the influence of the inclined angle (), half span of beam (), beam width, beam thickness, and shuttle dimensions on the static output displacement and force. For the dynamic behavior, we also calculated the relationship between rise time and beam dimensions for both vacuum and air conditions [27]. However, up to now, no research effort has been reported on the correlations between vibration frequencies and beam dimensions.

This paper aims to provide a comprehensive modal analysis and parameter study of the V-shaped beam electrothermal microactuators. First, we establish the equations of motion in the longitudinal and lateral directions, respectively. Then, the frequency equations are derived subjected to both longitudinal and lateral constraint conditions. By solving the frequency equations, the natural frequencies and the associated mode shapes can be calculated. Analytically calculated natural frequencies and related mode shapes are compared with the results of the finite element simulations using ANSYS software to verify the presented modal analysis. With this, we extend the modal analysis of the single pair of beam to multipairs of beams. For design purposes, we also theoretically proved that the natural frequencies of the V-shaped beam are inversely proportional to the volume scaling of the beam. Furthermore, parameter analysis is carried out to reveal the correlations between the natural frequencies and both the material properties and beam dimensions. These results provide insights on the vibration characteristic analysis and design of the V-shaped beams structures. This paper presents the first insight study on the modal analysis and parameter studies of the V-shaped beam electrothermal actuators. The aim is at laying the foundation of conducting dynamic thermomechanical modelling, dynamic design, and control of the V-shaped electrothermal microactuators.

This paper is organized as follows. In Section 2, the free vibration of the V-shaped beam is formulated incorporating both longitudinal and lateral deflections. The natural frequencies and the associated mode shapes are then calculated. The analytically computed natural frequencies and the mode shapes are compared with the results of finite element simulations using ANSYS software in Section 3. The modal analysis for of a single pair of beams is further extended to a multipair case. Section 4 presents the investigation into the relationship between the natural frequencies and the volume scaling of the V-shaped beams. In Section 5, we conduct detailed parameter analysis between the natural frequencies and both the material properties and the beam dimensions of the V-shaped beams.

2. Vibration Analysis and Equations

The structure of the V-shaped beam electrothermal microactuator is shown in Figure 1, where is the inclined angle, the length of the beam is measured as , the thickness and width of the beams are denoted as and , respectively, half span of the beam is , and and represent the length and width of the shuttle, respectively. When applying a voltage at the two ends, the temperature of the beam will increase due to Joule heating. As a result, the beam will expand to push the shuttle forward.

Apparently, if the shuttle is supported with a pair of traditional line-shaped beams, that is, the inclined angle equals to zero, the deflection of the beams is only in the lateral direction [40]. In the case of the V-shaped beam, however, due to the existence of the inclined angle, besides the lateral deflections, longitudinal beam deflections are also created for this symmetric structure. In what follows, both the longitudinal and lateral vibrations will be taken into account in the modal analysis of the V-shaped beam electrothermal microactuators. In this section, we only discuss the modal analysis for the single pair of the V-shaped beam case, which can be easily extended to the multipairs of beam cases with the mass of the shuttle being evenly distributed for each single pair of beams.

2.1. Equations of Motion

In this section, the equations of motion are formulated in both lateral and longitudinal directions. Because of symmetry, half of the V-shaped beam is considered, as shown in Figure 2(a), where the axis is along the longitudinal direction and the axis represents the lateral direction. and represent the lateral and longitudinal deflections, respectively. To derive the differential equation of motion of V-shaped beams, all the forces and moments acting on the cross sections of a small element of the beam are taken into account, as shown in Figure 2(b). is the axial force, is the shear force, and is the bending moment. In the following discussions, , , , and are the density, Young’s modulus, cross-sectional area, and moment of inertia of the beam, respectively.

The equation of motion for the free lateral vibration can be expressed aswhere . The equation of motion for the free longitudinal vibration of a uniform beam iswhere .

The free vibration solution can be found using the method of separation of variables as

Substituting (3) into (1) and (4) into (2), respectively, and rearranging lead towhere and in (7) are constants that can be found from the initial conditions which will not be discussed in this work. , , , and as well as and are different constants that can be found from the boundary and continuity conditions. and relate the natural frequency as

2.2. Boundary and Continuity Conditions

Since the equations of motion for the lateral and longitudinal vibrations in (1) and (2) involve the fourth-order and second-order derivatives with respect to , respectively, six boundary or continuity conditions are required to obtain the frequency equation. The boundary and continuity conditions depict the force and deflection relationships at and , as listed in Table 1. The boundary conditions summarized in the first four rows in Table 1 are self-explained. The last two rows are the continuity conditions, which describe the relationship between longitudinal and lateral deflections or forces. Due to symmetry, the movement of the shuttle is confined in the Y direction as shown in Figure 2(a), which require that the longitudinal and lateral deflections satisfy , as shown in Figure 3(a). Using separation of variables, we have at . Since the inertia force of half of the shuttle is , where is mass of the shuttle, the force equation of motion in the direction of the movement of the shuttle at can be given as

As shown in Figure 3(b), Equation (10) can be rewritten as

Inserting (3), (4), and (7) into (11), we obtain

2.3. Frequency Equation

Applying the boundary and continuity conditions into (5) and (6) leads to a system equation as in the following equation. The constants , , and in (5) and (6) equal to zero:

For the convenience of expression, the following notations are made: , , , , , and . The elements of the coefficient matrix in (13) are , , , , , , , and . Clearly, the coefficient matrix of (13) is only with respect to the natural frequency . For nontrivial solutions of , , and , the determinant of the coefficient matrix must be zero, which gives the frequency equation as follows:

If there is no shuttle in the middle, that is, , the frequency equation becomes

The frequency equations (14) and (15) are highly coupled, which do not allow one to obtain a simple analytical form of natural frequency . The value of must then be calculated numerically. Note that the nth order of natural frequency is written as .

2.4. Mode Shapes

Once the natural frequencies are obtained, the coefficient matrix of (13) is solvable. Then, it is convenient to write and with respect to as follows:

Substituting (16) and (17) into (5) and (6) and noticing that , , and equal to zero, the mode shapes of the lateral and longitudinal vibrations are calculated, respectively, as followswhere . It is seen that given the value of the natural frequency and , the shape of the beam is determined. In contrast to the traditional line-shaped beam cases, the mode shapes of the V-shaped beam involve the coupling of the lateral and longitudinal deflections. The mode shapes can be computed as in (20) and (21) with respect to the base frame in Figure 2(a). and are lateral and longitudinal deflections calculated with (18) and (19), respectively.where .

3. Finite Element Analysis

In this section, the natural frequencies and associated mode shapes of the V-shaped beams are calculated analytically and compared with finite element analysis using ANSYS software (version, 14.0). The 3-D SOLID227 element is utilized to allow for the three-dimensional modelling and analysis of the V-shaped beam. A set of beam dimensions as well as material properties used in calculations and simulations are listed in Table 2. Note that, in the preceding section, (known as the angular frequency) is called the natural frequency. In this section, is referred to as the natural frequency in unit Hz. Since the beam is of the continuous system with the infinite number of degrees of freedom, infinite orders of natural frequencies and associated mode shapes can be obtained. In this section, only the first three orders of natural frequencies and mode shapes are computed and compared, as summarized in Table 3. The percent deviations are calculated as follows: .

It can be seen from Table 3 that the analytical results agree well with the simulation results with a maximum derivation of 1.11% for predicting the natural frequencies of the V-shaped beams with a shuttle, and a maximum derivation of −0.44% for the V-shaped beam without a shuttle. The analytically calculated mode shapes are based on the combination of the calculations of the lateral and longitudinal mode shapes, as in (20) and (21). The mode shapes calculated analytically also match the FEA simulation results very well. Clearly, the natural frequencies for the case without a shuttle are higher than the results for the case with a shuttle. This can be explained in such way that increasing the mass results in decreasing the natural frequencies for such a spring-mass system by considering the pair of beam to be a spring carrying a mass. In this comparison, the asymmetric and out-of-plane vibrations are not taken into consideration since these cases are irrelevant to the direction of actuation. For the multipairs of beam case, since each pairs of beams are identical and symmetrically configured, the mode shapes for each pairs are the same and the natural frequencies can be calculated based on the single pair of beam with the mass of the shuttle being evenly distributed for each pairs of beams.

4. Volume Scaling

For MEMS machining process, electrothermal microactuators can be fabricated with the beam length ranging from hundreds to thousands of microns and even larger. One of the problems that has to be addressed at the early stage during the design process is to determine the volume scale of the V-shaped electrothermal microactuators. In general, in order to achieve larger displacement and force, it is desirable if we can just simply scale up the volume. However, this will decrease the electrothermal response, and it may not be suited to fast actuation applications. On the contrary, scaling up the volume will result in decreasing of the natural frequencies and the impact of the mechanical vibrations can increase.

In this section, we investigate the relationship between natural frequencies and volume scale. First, we introduce the scale factor , and the beam dimensions can be written with respect to the reference dimensions as , , , , , , and , where , , , , , and are called the reference beam dimensions. The cross-sectional area, moment of inertia, and mass of the shuttle of V-shaped beam corresponding to the reference beam dimensions are calculated as , , and , respectively. For a particular V-shaped beam with a factor , we have , , and .

Note that the inclined angle remains to be the same with the beam scaling up or down as . With this, (8) and (9) can be rewritten as

Writewhere and . We have

Denoting , , , , , and , we have , , , , , and . Clearly, since and are the functions of , , , , , , and are the functions of . The last row in (13) is expressed as

Considering (24) and (25), , , and can be rewritten aswhere , , and . Introducing (29), (30), and (31) into (28) leads to

Considering (32), (13) can be written with respect to the reference beam dimensions, and all the elements in the coefficient matrix are the functions of . Therefore, the determinant of the coefficient matrix must be a function of , which corresponds to the reference beam dimensions. That is, the frequency equations (14) and (15) are both functions of :

Comparing (14) and (15) with (33) and (34), we can conclude that is the natural frequency corresponding to the reference beam dimensions, that is, . Hence, we havewhere the subscript indicate the orders of the natural frequencies. Equation (35) is important for predicting the natural frequencies of the V-shaped beam actuator. If the natural frequencies of one particular V-shaped actuator are given, the natural frequencies for the scaled actuator with the scale factor can also be obtained. Clearly, the natural frequencies are inversely proportional to the volume scale for the same dimension proportions and same material properties.

This relationship is demonstrated in Figure 4, (35). Figures 4(a) and 4(b) are based upon the V-shaped beam dimensions listed in Table 2, which correspond to the scale factor . The material properties used in Figures 4(a) and 4(b) are silicon (listed in Table 2) and SU-8, respectively. Young’s modulus and density of SU-8 are and . It is seen that the calculated natural frequencies are indeed inversely proportional to the volume scale . This suggests that when the V-shaped beam is scaled up, the impact of the mechanical vibration on the dynamics of the V-shaped beam increases. Compared to silicon, the natural frequencies of the V-shaped beam made of SU-8 is nearly 10 times smaller. This indicates that the material properties also have a large impact on the vibrations of the V-shaped beam microactuators. In what follows, we will discuss in detail the influence of the material and dimensions of the structure on the natural frequencies.

5. Parameter Analysis

In this section, we investigate the influence of material properties and beam dimensions on the natural frequencies in order to provide insight and guidance on the further design of the V-shaped beam electrothermal microactuators.

5.1. Material Properties

It is seen from the frequency equation that only Young’s modulus and density of the actuator determine the natural frequencies if beam dimensions are fixed. may be called the characteristics material property for vibration. For a reference characteristics material property , relates to with a scale factor as follows:

Introducing (36) into (8) and (9) leads toand comparing (39) and (40) with (8) and (9), we notice that and are the functions of for characteristics material properties , whereas and are the functions of for the reference characteristics material properties . Obviously, we have

Denoting , , , , , and , we have , , , , , and . The mass of the shuttle is calculated as , in which is the volume of the shuttle. For the reference density , mass of the shuttle is calculated as . Clearly,

Applying (43) into (28), , , and can be expressed aswhere , , and . Introducing (44), (45), and (46) into (28) leads to

Considering (47), all the elements in the coefficient matrix in (13) are the functions of , and therefore, the determinant of the coefficient matrix, that is, the frequency equation, is the function of . The frequency (14) and (15) are expressed as

Comparing (48) and (49) with (14) and (15), we can conclude that represents the natural frequencies corresponding to the reference characteristics material properties, that is, . Hence, we havewhere the subscript indicate the orders of the natural frequencies. Equation (50) demonstrates that the natural frequencies are inversely proportional to the scale factor . Combining (35) and (50) leads towhere represents the natural frequencies associated to the reference beam dimensions and characteristics material properties. Though and are both called the scale factor, they are different. refers to the scaling of the volume of the actuator, and refers to scaling the value of the characteristics material property .

Table 4 lists six materials that are commonly used for fabrication of electrothermal microactuators. The characteristics material property of the polysilicon is the smallest and was considered as the reference material. The materials relate the polysilicon with the scale factor from 1.00 to 8.70 according to the relative value of the characteristics material properties .

Figure 5 shows the natural frequencies with respect to the scale factor , in which the solid line represents the natural frequencies derived from (50) with the natural frequencies of polysilicon as the reference value. The asterisks represent the analytically calculated natural frequencies of particular material properties (as listed in Table 4) based upon the modal analysis in Section 2. The beam dimensions are chosen as in Table 2. Clearly, the asterisks fall exactly on the solid lines as predicted by (50). For the same beam dimensions, the natural frequencies of V-shaped beam made of polysilicon are the largest, meaning that, in this case, the dominant dynamics is determined by the electrothermal part and the influence of the mechanical vibration is small and negligible. The V-shaped beam made of SU-8 is most likely to undergo mechanical vibrations since the natural frequencies are the smallest.

5.2. Structural Parameters

In this section, the effect of varying the structural parameters, that is, , , and , of the V-shaped beam on its natural frequencies is studied. The thickness of the V-shaped beam is not studied since it is irrelevant to the natural frequencies as demonstrated in (14) and (15). In the following discussions, we consider a single pair of beam with a shuttle for illustration, and similar conclusions can be extended to the cases of multipairs of beam and beam without a shuttle.

The effect of changing on the first-order natural frequencies is shown in Figure 6. The natural frequency increases to the maximum value followed by a drop as increases. In the meantime, the mode shape also changes with increasing , as shown in Figure 7. Obviously, the motion of the shuttle contributes to the vibration mode for less than around 60 . As increases, the vibration mode of the structure tends not to involve the motion of the shuttle. For the first case, as shown in Figure 7(a), the V-shaped beam can be thought of as a spring-mass system in which the mass of the shuttle represents the equivalent mass of the system and the stiffness of that system is the stiffness in the direction of the motion of the shuttle. Apparently, increasing leads to larger longitudinal deflections thus increasing the stiffness of the system. For a spring-mass system, increasing the stiffness will lead to increase in the natural frequencies. For the latter case, as shown in Figure 7(b), the V-shaped beam can be thought of as a fixed-fixed line-shaped beam for the half of the structure. With an increase of , the stiffness for this fixed-fixed beam decreases resulting in decreasing of the natural frequencies. Clearly, Figures 6 and 7 demonstrate that varying will not only influence the natural frequency but also influence the mode shape.

The effect of varying on the first-order natural frequency is illustrated in Figure 8. The green and the red curves represent the two classes of vibration modes, that is, the spring-mass and the fixed-fixed cases, respectively. For both vibration modes, as shown in Figure 8, the first-order natural frequency decreases with increasing . Obviously, increasing results in the decrease of the stiffness of the structure. Figure 9 shows the first-order natural frequencies with respect to beam width . The red and green curves represent the two vibration modes. Clearly, for both the modes, increasing the beam width results in stiffening the beam thus increasing the natural frequencies.

6. Conclusions

A comprehensive analytical modal analysis and a parameter study of V-shaped beam electrothermal microactuators have been conducted for the first time in this work. Both the longitudinal and lateral vibrations have been included in deriving the natural frequencies and the corresponding mode shapes of the V-shaped beam. Analytically obtained natural frequencies and mode shapes were compared to the results of finite element simulations using ANSYS to verify the analytical modal analysis. In addition, this paper has analytically proved that the natural frequencies are inversely proportional to the volume scale and the scale factor of characteristic material properties. This provides the guidance on how to determine the volume scale and to choose the materials of the V-shaped beam electrothermal microactuators. Finally, the parameter analysis was conducted to provide insights on the design and optimization of the V-shaped electrothermal microactuators.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (no. 51575006).