JOURNAL OF MICROELECTROMECHANICAL SYSTEMS, VOL. 16, NO. 4, AUGUST
|
|
- Juliet Waters
- 5 years ago
- Views:
Transcription
1 JOURNAL OF MICROELECTROMECHANICAL SYSTEMS, VOL. 16, NO. 4, AUGUST Influence of Boundary Conditions on the Dynamic Characteristics of Squeeze Films in MEMS Devices Ashok Kumar Pandey, Rudra Pratap, Member, IEEE, and Fook Siong Chau Abstract Micromechanical structures that have squeeze-film damping as the dominant energy dissipation mechanism are of interest in this paper. For such structures with narrow air gap, the Reynolds equation is used for calculating squeeze-film damping, which is generally solved with trivial pressure boundary conditions on the side walls. This procedure, however, fails to give satisfactory results for structures under two important conditions: 1) for an air gap thickness comparable to the lateral dimensions of the microstructure and 2) for nontrivial pressure boundary conditions such as fully open boundaries on an extended substrate or partially blocked boundaries that provide side clearance to the fluid flow. Several formulas exist to account for simple boundary conditions. In practice, however, there are many micromechanical structures such as torsional microelectromechanical system (MEMS) structures that have nontrivial boundary conditions arising from partially blocked boundaries. Such boundaries usually have clearance parameters that can vary due to fabrication. These parameters, however, can also be used as design parameters if we understand their role on the dynamics of the structure. We take a MEMS torsion mirror as an example device that has large air gap and partially blocked boundaries due to static frames. We actuate the device and experimentally determine the quality factor Q from the response measurements. Next, we model the same structure in ANSYS and carry out computational fluid dynamics analysis to evaluate the stiffness constant K, the damping constant D, and the quality factor Q due to the squeeze film. We compare the computational results with experimental results and show that without taking care of the partially blocked boundaries properly in the computational model, we get unacceptably large errors. [ ] Index Terms Computational fluid dynamics (CFD), experimental quality factor, microelectromechanical system (MEMS), partially blocked boundaries, squeeze-film damping, torsional motion. I. INTRODUCTION EVALUATION of fluid damping and spring characteristics is critical in the design and analysis of dynamic microelectromechanical system (MEMS) devices where the motion Manuscript received September 26, 2006; revised May 1, This work was supported in part by the Department of Science and Technology (DST), New Delhi, India, and in part by the National University of Singapore, Singapore. Subject Editor N. Aluru. A. K. Pandey and R. Pratap are with the Department of Mechanical Engineering, Indian Institute of Science, Bangalore , India ( ashok@mecheng.iisc.ernet.in; pratap@mecheng.iisc.ernet.in). F. S. Chau is with the Department of Mechanical Engineering, National University of Singapore, Singapore ( mpecfs@nus.edu.sg). Color versions of one or more of the figures in this paper are available online at Digital Object Identifier /JMEMS of the structure is in a fluid medium such as air or some other gas. For MEMS devices such as a torsional mirror or an accelerometer, in which the main sensing/actuating element executes transverse motion with respect to a fixed substrate, fluid damping is mainly due to the squeeze film and drag on the lateral and top surfaces [1], [2]. The squeeze-film damping in MEMS devices is generally calculated by solving the conventional Reynolds equation. The analytical results for the squeeze-film damping and spring forces in a system with a plate executing small amplitude harmonic motion are given in [3], and the corresponding torques in the case of torsional motion are derived in [4]. These damping and spring forces in complex structures are also calculated by solving Reynolds equation numerically for rigid [5] and flexible structures [6] executing transverse motion. Since the computation methods are intense and time consuming, Hung et al. [7], [8] have proposed a low-order model for fast dynamical simulation by extracting basis functions from a few finite-element method (FEM) simulation results. In all these studies, the Reynolds equation is solved with trivial pressure boundary condition (P =0)on the free edges and trivial flow conditions (V =0)on the fixed edges. Although the trivial pressure boundary condition hardly represents the actual condition around a moving structure, it is helpful in predicting the nature of squeeze-film effects for structures with large length to air gap ratio. For comparable length to air gap ratio, the effect of surrounding pressure on the squeeze-film damping can be significant. Thus, a nontrivial boundary condition should be used while calculating the effect of the squeeze film. There are a few studies ([9] [12]) that present various ways to consider nontrivial boundary conditions on squeeze-film damping under various conditions. To apply a nontrivial pressure boundary condition on the plate edges, Vemuri et al. [9] have performed numerical simulation on a 3-D domain by extending the control volume beyond the edges of the oscillating plate. The extension to the control volume is taken such that it is greater than the distance from the plate edges where the pressure difference becomes zero. Subsequently, the extended length and width of the plate are found by comparing the numerical solution with the analytical damping and spring constants obtained using the trivial boundary condition for the effective plate with elongated length and width. The required extension to the plate dimensions is a function of the characteristic dimension of the squeeze film, i.e., the characteristic length to air gap ratio, and the frequency of operation. Veijola et al. [10] have proposed two /$ IEEE
2 894 JOURNAL OF MICROELECTROMECHANICAL SYSTEMS, VOL. 16, NO. 4, AUGUST 2007 Fig. 1. (a) Side view and (b) top view of the microfabricated cantilever beam and torsional structure are shown with the partially blocked boundaries. Fig. 2. Schematic diagram of the (a) double-gimballed torsion mirror and (b) its sectional view along A A. To evaluate the effect of different factors on fluid damping in the structure shown in (a) and (b), an appropriate boundary condition is used. The extended boundaries (shown dotted) are taken in four different ways to capture different effects: (c) squeeze-film effect; (d) squeeze-film effect, end effect, and drag effect only on the lateral side of the vibrating structure; (e) squeeze-film effect, end effect, and drag effect on the side and the top of the oscillating structure; and (f) in addition to the effects captured in (e), the effect of a static frame separated by a small distance from the boundaries of the oscillating structure is also considered. models to consider the effect of open boundaries on squeezefilm damping: 1) the border flow channel model where imaginary flow channels are used at the edges of the plate; and 2) the surface extension model where the edges are extended to accommodate the effect of open borders. They have calculated the appropriate extension to the original dimensions of the plate by comparing the numerical squeeze-film resistance obtained using trivial pressure boundary conditions on the extended
3 PANDEY et al.: DYNAMIC CHARACTERISTICS OF SQUEEZE FILMS IN MEMS DEVICES 895 volume with the analytical resistance value for the extended plate. They have found that in the continuum flow regime, the validity of Reynolds equation can be extended to systems with characteristic structural dimensions comparable to the air gap thickness if the structural dimensions are increased by 1.3 times the air gap thickness. The model is also modified to include the rarefaction effect in the slip flow regimes. In torsional motion, the maximum relative error under this condition is found to be less than 10% for the length-to-air-gap ratio greater than 10. Since the dimensional extension is obtained under the assumption of incompressible flow, the model is reported to be valid below 500 khz. The above studies are mainly done for fully open boundaries, as shown in case 2 of Fig. 1(b). In order to model the effect of partially blocked boundaries, where the fluid flows through a narrow opening on the side of the plate, as shown in case 1 of Fig. 1(b), Darling et al. [12] have proposed a compact model under arbitrary venting conditions by solving acoustic equations. They have obtained an analytical solution for the squeeze-film damping and spring forces by solving the conventional Reynolds equation for trivial as well as nontrivial boundary conditions. In order to consider the effect of side walls, Hao et al. [11] have obtained the squeeze-film damping and spring torques for torsional motion by taking trivial pressure boundary conditions on the side walls rather than on the edges of the oscillating plate. However, the zero-pressure boundary condition on the side walls is not a good idea when the side walls are placed closer to the oscillating plate. The effect of nontrivial boundary conditions in partially blocked boundaries requires thorough investigation in order to find the distance from the edges of the oscillating structure where the pressure difference becomes zero. The models presented in the last two references [1], [12] are restricted by the upper frequency limit below which inertial effects are negligible. Also, these models do not apply to cases with large air gaps (where the Reynolds equation is not applicable due to the violation of its basic assumption of small air gap to length ratio). These studies are also not applicable to structures with complex boundary conditions such as the partially blocked boundaries that exist in the MEMS device under consideration [which is shown in Fig. 2(a)]. Partially blocked boundaries have clearance parameters that can have small enough values to be significantly affected by fabrication processes. The gap between a side wall (a fixed boundary) and a vibrating structure is one such parameter. The geometric parameters that characterize the fluid domain between the fixed boundaries and the moving structure are of interest here. These parameters can also be used as effective design parameters if we understand their role on the dynamics of the structure. In this paper, we take an example structure with partially blocked boundaries and carry out computational fluid dynamics (CFD) simulations using several models for boundary conditions to evaluate the dynamic characteristics of the structure. In particular, we compute the quality factor with each of the boundary models and compare the results with experimentally obtained values. We also compare the numerical and experimental results with the available analytical formulas and discuss the validity of analytical formulas under different conditions. II. THEORY A. Governing Equations In a system where a structure vibrates near a fixed substrate in a fluid medium, the fluid behavior can be modeled using the Navier Stokes equation [14], [15]. For a low vibration frequency, assumptions of incompressibility and laminar flow hold good. However, as the excitation frequency increases, compressibility and inertial effects become significant. Under isothermal flow conditions and constant dynamic viscosity, we get the following set of equations for fluid flow [14], which are solved using ANSYS: ρ +.(ρ u) =0 (1) t ρ D u Dt = ρ g p + µ 2 u (2) ρ p (3) where u is the fluid velocity, p is the fluid pressure, ρ is the fluid density, g is the acceleration due to gravity, and µ is the dynamic viscosity. In problems related to MEMS devices, the effect of body forces such as ρ g can be neglected without any substantial error. The time rate of change of density can be expressed in terms of the time rate of change of pressure as ρ/ t =( ρ/ p)( p/ t), where ρ/ p = ρ/κ, and κ is the bulk modulus of the fluid. In this paper, the fluid considered is air. For air, at temperature T = 298 K and pressure p a = Pa, µ N-s/m 2, ρ a =1.2 kg/m 3, and the isothermal bulk modulus κ T = p a = Pa. We must add here that at sufficiently high frequencies or velocities, adiabatic conditions may have to be used in place of isothermal conditions. B. Boundary Conditions Based on the type of opening on the side of the vibrating structure, the boundary conditions can be mainly of the following three kinds. 1) Fully closed boundaries: In this case, the fluid remains inside the cavity. For pure transverse motion, the fluid behaves like a spring, while for torsional motion, there is also some loss due to the lateral fluid motion in the air gap. The magnitude of this effect depends on the frequency of excitation [12]. Under this condition, the fluid velocity along the boundaries is zero. 2) Fully open boundaries: In this case, the fluid can flow freely from the boundaries, and the pressure variation is taken as zero along the boundaries, i.e., p =0.However, because of the sudden flow expansion from the closed cavity to the open surroundings, the pressure on the boundary is different from the ambient pressure. This difference in pressure is called the open border effect or the end effect [17]. This effect is negligible for a system with small air gap h 0 to length L (or width W ) ratio. Under this condition, trivial pressure boundary conditions can be used to calculate damping. For incompressible
4 896 JOURNAL OF MICROELECTROMECHANICAL SYSTEMS, VOL. 16, NO. 4, AUGUST 2007 Fig. 3. Picture of double-gimballed MEMS torsion mirror. flow, the damping constant under this condition is given by [4] C φ = 192WL5 µ π 6 h 3 0 m n even odd 1 (mn) 2 (m 2 + n 2 ( L W ) 2 ) (4) where µ is the fluid viscosity. For comparable air gap and structural length, the open border effect can be modeled by applying either a nontrivial boundary condition on the original boundaries or a trivial boundary condition on the extended boundaries of the oscillating structure. For low-frequency torsional motion about an axis bisecting the plate length, Veijola et al. [10] have proposed the extension of the original length l of the oscillating plate by 1.3h 0, i.e., l eff = l +1.3h 0, where h 0 is the initial air gap thickness. Consequently, for an air gap to length ratio of 1 in a structure with fully open boundaries, (4) can be used to calculate damping by replacing L and W with L eff = L +1.3h 0 and W eff = W +1.3h 0.After calculating the damping force in the two different cases, the value of the moment of inertia I and the angular velocity ω of the original structure are used to calculate the quality factor from Q anal = Iω. (5) C φ 3) Partially open boundaries: Unlike the previous cases, the boundaries are neither completely open nor completely closed. In this case, αp + β p =0is used on the boundaries, where α and β are constants whose values depend on the closed and open portions of the boundaries. C. Solution Methodology To study the effect of boundaries on the dynamic characteristics of the squeeze film, we choose a MEMS torsion mirror (shown in Fig. 3) as a test structure whose dimensions are given in Table I. Our methodology involves experimental and numerical investigations that include different boundary effects. In this section, we outline the procedure we have followed for these studies. 1) Experimental Procedure: In order to evaluate the quality factor of the torsional microstructure under ambient conditions, we build a test setup shown schematically in Fig. 4. The dynamics of the torsional structure involves rotational vibration of the inner plate about the inner torsional bars. The experimental studies in this case require a function generator, an oscilloscope, two dc power sources, a LASER source, a LASER pointer, a position sensing device (PSD), a platform, and three adjustable stages to mount the PSD, LASER pointer, and test structure, respectively. The procedure of determining the quality factor experimentally involves the following steps. 1) An input voltage V = V input sin 2πft, where f is the input frequency, and t is the time, and a bias voltage V bias are applied to the electrodes (E 1 and E 4 ) and the inner plate [see Fig. 2(a)], respectively. Here, V bias =20 V, and V input =0.5 V. The static angular deflection due to the applied bias voltage is about rad, which is negligibly small. The application of a sinusoidal voltage results in angular oscillation φ of the inner plate about the inner torsional bars. 2) A LASER source and a LASER pointer are used to focus the LASER beam at a point on the top of the vibrating plate. After getting reflected from the plate, the LASER beam falls on the PSD, as shown in Fig. 5(a). 3) As the plate vibrates, the voltage reading of the PSD is taken from an oscilloscope. At steady state, the voltage reading is taken, which in turn gives the angular displacement of the oscillating plate from the relation φ =(V p /d)(in rad). Here, V p is the PSD output voltage (in the present case, the relation between voltage and displacement is set as 1 V 1 mm), and d (in millimeters) is the center distance between the PSD and the LASER spot on the vibrating plate. 4) By repeating steps 1 3 for different frequencies and by reading the corresponding maximum angular displacements, we get the frequency response curve as shown in Fig. 5. 5) Using the half-width method [20], we measure the quality factor Q exp from the frequency response curve. The expression for quality factor is given by Q exp 1 2ξ = f d f 2 f 1 (6) where f 1 and f 2 are the half power frequencies. By following the above experimental procedure, we get the response curve for the angular oscillation of the inner plate, as shown in Fig. 5. The experiment is repeated ten times, and the average of all data points is taken. From the response curve, we get the quality factor Q exp 9.19 and the damping factor ξ =0.054 from f 1 = Hz and f 2 = 558 Hz. The damped natural frequency comes out to be f d = Hz. The maximum angular displacement at resonant frequency is found to be rad. 2) Numerical Procedure: The CFD procedure to calculate the dynamic characteristics of the squeeze film under different boundary conditions involves three major steps: 1) modeling and discretizing the fluid domain; 2) solving for the pressure
5 PANDEY et al.: DYNAMIC CHARACTERISTICS OF SQUEEZE FILMS IN MEMS DEVICES 897 TABLE I DIMENSIONS OF MEMS TORSION MIRROR SHOWN IN FIG.2 Fig. 4. Outline sketch of the experimental setup. distribution with the given load and boundary conditions using (1) (3); and 3) computing the desired forces and other quantities from the computed pressure distribution and velocity field. An outline of the different computational substeps is given in Fig. 11. A few important aspects and expressions related to modeling and estimation of the spring force, damping force, and quality factor are discussed below. 1) Issues related to modeling: First, a 3-D fluid domain around the test structure is created. The outer boundaries, which are shown dotted in Fig. 2(f), are extended from the side and the top surface of the oscillating structure. The extension factor depends on various factors such as the air gap thickness to the length ratio [10] and the Fig. 5. Frequency response (in terms of angular displacement φ) ofmems torsion mirror vibrating about the y-axis for bias voltage V bias =20Vand V input =0.5V. nature of the boundary condition such as fully or partially blocked boundaries. The reasons for extending the boundaries in the open and partially blocked conditions are mainly twofold: first, to accommodate the effect of sudden expansion and contraction through the side and top openings; and second, to account for the drag forces on the side walls and the top surface of the oscillating structure (see Fig. 2(a) and (b) for the openings around the inner plate). The main criteria for selecting the extension of surfaces is that it should be greater than the distance where the solution converges for pressure distribution
6 898 JOURNAL OF MICROELECTROMECHANICAL SYSTEMS, VOL. 16, NO. 4, AUGUST 2007 Fig. 6. (a) Comparison between the numerical and analytical damping constants for different air-gap thicknesses. (b) Percentage error in analytical result with respect to numerical result with different air-gap thickness. around the structure. For example, the boundaries in Fig. 2(e), which are similar to the condition around the experimental test structure, are extended by l p /2 on the left and right sides of the oscillating plate, and 3h 0 /4 from the upper surface of the vibrating plate on the top to capture the aforementioned effects, where l p is the length of the plate, and h 0 is the air-gap thickness. 2) Calculation of the quality factor: After getting the pressure distribution p on the walls of the oscillating plate and integrating the moment of the total back force about the axis of oscillation, we estimate the total back torque on the oscillating structure. Then, after calculating the phase difference γ between the total back torque and the angular velocity, we calculate the spring torque, the damping torque, and the quality factor based on the total fluid damping. If the angular displacement of the oscillating structure is φ(t) =Φsinωt, where ω =2πf =2π/τ is the angular frequency of oscillation, τ is the time period of oscillation, and γ is the phase difference as defined above, then the angular velocity and the total torque can be written as φ(t) =Φω cos ωt and T = T 0 cos(ωt + γ), respectively. Under the assumptions of viscous damping and linear spring, the spring and the damping torques are proportional to the angular displacement and the angular velocity, respectively. Therefore, the expressions for the damping and spring torques are given by T d (t) =T d0 cos ωt T s (t) =T s0 sin ωt (7) where T d0 = T max cos γ, and T s0 = T max sin γ. Thus, the expressions for the damping constant D and the spring constant K are D = T d0 Φω Now, the quality factor is defined as [20] 1 2 IΦ2 ω 2 Q num =2π τ 0 T d(t) K = T s0 Φ. (8) φ(t)dt = IΦω2 T d0 = Iω D (9) where I is the moment of inertia about the axis of rotation. For the present structure, as shown in Fig. 2(a) and (b), I = ρl p b p t p (lp 2 + t 2 p)/12 without perforations, and ρl p b p t p (lp 2 + t 2 p)/12 j=n j=1 (nρl hl h t p (lh 2 + t2 p)/12 k=+(n+1)/2 k= (n+1)/2 (ρl2 h t pk 2 p 2 h )) with n n equally spaced perforations, where n is an odd number and equals 13 in the present case. The I value corresponding to n =13is about kg-m 2. The net torsional stiffness of the structure is from as k φ = ω 2 I =(2πf d ) 2 I = N-m, where f d = Hz is the resonant frequency obtained from Fig. 5. III. RESULTS AND DISCUSSION First, we discuss the validation of the numerical results by comparing them with the experimental results for the MEMS torsion mirror shown in Fig. 2. We discuss various effects on the dynamic performance of the test structure associated with
7 PANDEY et al.: DYNAMIC CHARACTERISTICS OF SQUEEZE FILMS IN MEMS DEVICES 899 Fig. 8. Variation of percentage error in calculating T max with the number of elements N e (here, T ref = N-m is the converged value of T max corresponding to elements). Fig. 7. (a) Finite elements in the domain bounded by the bottom surface and the section A A. (b) Sectional view of the fluid domain passing through the middle of the oscillating plate that is bounded by the top surface, side surfaces, and bottom surface. different boundary conditions. Before going into the details of numerical modeling, we characterize the nature of the air flow around the given structure using popular characteristic numbers. 1) Modified Reynolds number Re = ρh 2 ω/µ: It is a measure of the inertial effects. Using the characteristic flow length h = h 0 =80µm, the dynamic viscosity of air µ = N-s/m 2, the density of air ρ =1.2 kg/m 3, and the angular frequency at resonance ω =2πf, where f = Hz, we get Re = 1.4. Here, the Reynolds number is greater than 1 because of the large air gap thickness. Under this condition, Reynolds equation, which is used to model the squeeze-film flow, does not hold good. To incorporate the effect of the large air gap thickness, we solve the 3-D Navier Stokes equation numerically. 2) Squeeze number σ =12µl 2 pω/h 2 0p a :Forl p = 400 µm, p a = Pa, and the same values of h 0, ω, and µ as mentioned earlier, we get σ = Since σ 1, the compressibility effect can be neglected. 3) Knudsen number Kn = λ/h: For ambient pressure and temperature, the mean free path λ of the air molecules is about 64 nm. For the flow through the air-gap thickness, the characteristic flow length h = h 0, and Kn = ; whereas for the flow through the side opening, h = h w =19µm, and Kn = Since for both flows Kn < 0.01, the fluid flow around the structure is in the continuum flow regime [21]. Therefore, the no-slip boundary conditions at the boundaries of the oscillating structure can be used to solve the Navier Stokes equation in the fluid domain around the test structure. 4) Perforation ratio P r = A perf /A plate :Itistheratioof the perforated area to the plate area. For the perforation area A perf = N p 3 µm 3 µm, where N p =13 13 is the number of perforations, and the plate area A plate = 400 µm 400 µm, we get P r = Moreover, the ratio of the perforation size to the air-gap thickness P h = l h /h 0 =0.037 is also very small. Since the perforation size is very small compared to the pitch of the perforation distribution and the air gap thickness, we neglect the effect of perforations in the modeling. A. Validation of Numerical Procedure To validate the numerical procedure, we compare the numerically computed damping constant with the analytical value obtained from (4) for L = 400 µm, W = 400 µm, and h 0 = 10 to 80 µm. Air is used as the surrounding medium. The density and viscosity of air are taken to be 1.2 kg/m 3 and N-s/m 2, respectively. While the excitation frequency is Hz, the amplitude of angular displacement is taken as rad, which corresponds to 1% of the maximum allowable angular displacement when the tip of the electrode touches the substrate. The numerical damping constant is obtained by solving the Navier Stokes equation in 3-D domain with the trivial boundary condition on the plate edges as explained in the previous section. It is found that the analytical results compare very well with the numerical result for h 0 /L 0.1 with less than 10% error as shown in Fig. 6. Beyond 40-µmgap thickness, the effect of large air gap increases due to the nonuniform pressure variation across the gap thickness. At 80 µm, the analytical formula gives an error of about 33%. Thus, we validate our numerical procedure that matches well with the analytical result for small air gap to length ratios. Subsequently, we do further comparisons with the experimental results under different boundary conditions in the next section. Another important aspect in numerical modeling is the selection of the oscillation amplitude, which should be comparable to the experimental values. The experimental value of the oscillation amplitude is rad at Hz (see Fig. 5), which is about 8.75% of the maximum allowable angular amplitude ( 0.4 rad when the tip of the electrode touches the substrate). In numerical simulations, we vary the oscillation amplitude from 1% (0.004 rad) to 10% (0.04 rad) of the maximum allowable angular displacement and compute the corresponding damping constant. In this range, we find the damping constant to be independent of the amplitude as expected (for small amplitude motion) with a variation as small as 0.1%. Based on these results, we select rad as the amplitude of angular
8 900 JOURNAL OF MICROELECTROMECHANICAL SYSTEMS, VOL. 16, NO. 4, AUGUST 2007 Fig. 9. (a) Pressure distribution on the top surface of the fluid domain of Fig. 7(b) (here, the pressure variation is found to be zero). (b) Pressure distribution on the bottom surface. (c) Pressure distribution on the oscillating plate and the inner frame in the domain bounded by the bottom surface and section A A shown in Fig. 7(b) (black marks around the structure represent the edges of FLUID142 elements cut by the section A A ). displacement for all simulations in Section III-B, which is 1% of the maximum allowable angular displacement (0.4 rad). B. Comparison of Numerical, Experimental, and Analytical Results By following the numerical procedure outlined earlier, we calculate the quality factor of the MEMS torsion mirror. The mirror is subjected to an angular displacement φ =Φsin2πft, where f = Hz is the damped natural frequency at atmospheric pressure as found from the experimental result showninfig.5,andφ=0.004 rad is the amplitude of angular displacement, which is taken as 0.01 times the maximum possible angular displacement (when the plate touches the bottom electrode at the far end). Since the bottom electrode thickness (0.3 µm) and the static angular displacement ( rad) due to the bias voltage are very small, we neglect their effects on the nominal gap thickness h 0. From the numerical results, we find that for the present structure, the perforations in the plate can be neglected without any significant loss in accuracy in calculating the total reactive torque and the quality factor. The percentage error introduced by neglecting perforations is 0.087% and 0.79% in the evaluation of total torque and quality factor, respectively. This is largely because the perforation size (3 µm) is very small compared to the pitch ( 32.5 µm) and air gap thickness (80 µm) [22], [23]. Therefore, in the subsequent discussions, we neglect the perforations in the modeling. 1) Comparison of Numerical and Experimental Results: To validate the numerical procedure, we first create a 3-D fluid domain around the torsional mirror of the given dimensions (see Table I). The lateral boundaries of the fluid domain are fixed at l p /2 from the longer external edge and l p /10 from the shorter external edge of the inner frame, and the top is fixed at 3h 0 /4 from the top of the movable plate (see Fig. 7). The extension is obtained by gradually moving the boundary outwards until the obtained solution converges for the pressure distribution around the structure. The fluid domain is meshed with tetrahedral FLUID 142 elements of ANSYS. The zero-pressure boundary condition is applied on the extended surfaces; zero velocity condition is applied on the bottom surface and the surfaces around the static inner frame; and the sinusoidal angular velocity is applied to the surfaces binding the moving structure. Transient response calculations are carried out for six cycles because six cycles are found to be sufficient to attain steady state. On the time axis, 30 points are used per cycle. At each time step, local iterations are done to achieve convergence of order five [i.e., (P i+1 P i )/P i 10 5, where i denotes the local iteration step] in the calculation of pressure variation. The under-relaxation factors for iterating velocity and pressure are taken as 0.5 to ensure a fast and smooth convergence [24]. After calculating the total reactive torque and subsequently breaking it into its damping and spring components, we calculate the quality factor from (9). To ensure the accuracy of the numerical method, we calculate the maximum torque T max around the oscillating structure for different finite-element discretizations, which are used to mesh the fluid domain around the structure as shown in Fig. 7. Convergence of the finite-element calculations is clearly seen by plotting the percentage error against the number of elements (Fig. 8), where the error is computed with respect to the converged value of T max.ifn e is the number of elements, then the percentage error seems to reduce as 1/N 2 e (see inset in Fig. 8). The numerical error is below 1% if the number of elements is over Fig. 9 shows the pressure distribution on the top surface, the bottom surface, and the section A A [see Fig. 7(b)] along the bottom surface of the oscillating structure in the fluid domain, which is meshed with elements. The converged values of the maximum torque T max and the damping torque T d0 are N-m and N-m (= 100 µn-nm), respectively, which correspond to the angular oscillation of φ = sin( t). Finally, we get Q num =9.27, which is very close to the experimentally determined quality factor Q exp =9.19, givingan error of about 0.87%. The effect of squeeze-film flow between the vibrating plate and the fixed plate, and the effect of air flow through the side walls and the top surface of the vibrating plate, are compared and analyzed for different boundary conditions, as described in Fig. 10. In case (a), the zero-pressure boundary condition is taken on the edges of the vibrating structure where only squeeze-film damping is considered. In case (b), the zeropressure boundary condition is applied on the extended boundary so as to include the open border effect and the effect of the side walls on fluid flow (although the effect of the side wall is negligible in the present case). In case (c), the effect of the side and top surfaces on fluid flow is taken by applying the zero-pressure boundary condition on the surfaces extended from the side and top surfaces of the vibrating plate. Finally,
9 PANDEY et al.: DYNAMIC CHARACTERISTICS OF SQUEEZE FILMS IN MEMS DEVICES 901 Fig. 10. Patterns of average pressure variation and quality factor of a 3-D structure executing torsional motion under open and partially blocked boundary conditions. Different boundary conditions in (a) (d) are used to capture different effects as mentioned in Fig. 1(c) (f). Here, trivial pressure boundary condition is applied on the dotted lines, zero velocity condition is applied on the lines around the blank rectangle, and sinusoidal velocity is applied on the lines around the grey-colored rectangle. Here, the movement of mesh is scaled-up for better view. TABLE II PERCENTAGE ERROR (Q num Q exp/q exp 100) OBTAINED UNDER DIFFERENT BOUNDARY CONDITIONS WITH RESPECT TO EXPERIMENTAL VALUE OF QUALITY FACTOR Q exp =9.19 FOR THE STRUCTURE SHOWN IN FIG. 1 TABLE III ERROR OBTAINED BY USING CONVENTIONAL ANALYTICAL SOLUTION WITH OR WITHOUT OPEN BORDER EFFECT FOR THE GIVEN STRUCTURE.HERE, error = Q anal Q exp/q exp 100 we imitate the actual boundary conditions that exist around the MEMS device (shown in Fig. 3) in case (d) by putting a static frame near the boundaries of the vibrating plate. In this case, the free flow of the fluid from the boundaries of the vibrating plate gets partially obstructed by the presence of the nearby static structure. In Table II, the quality factor calculated numerically for the aforementioned cases is compared with the experimental value of the quality factor Q exp =9.19 for the device shown in Fig. 3. It is found that the condition in case 1, based on which most of the compact formulas are derived, gives an error of about 203%, while the case 3 that considers the effect of the side and top surfaces (i.e., fully open boundary condition) of the vibrating structure reduces the error to about 11%. However, when the actual condition is considered, the error reduces even further and gives a value of about 1% only. Although the boundary condition considered in case (d) is specific to the present MEMS device, it represents a more generic case of partially blocked boundaries. It is clear from the results that such boundaries need to be modeled carefully for reasonably accurate results. 2) Comparison of Analytical Models With Experimental Results: In Table III, the experimental result is compared with the analytical solution obtained by solving the conventional Reynolds equation [15] with zero-pressure boundary condition on the boundaries of the oscillating plate [4]. The relevant dimensions are taken from Table I. In case 1 of Table III, the zero-pressure boundary condition is applied on the original boundary of the oscillating inner plate, as shown in case (a) of Fig. 10, thereby ignoring the open border effect. The analytical quality factor is calculated using (4) and (5). It is found that the analytical result gives an error of about 356% when compared with the experimental result. The error in this case is greater than that in case 1 of Table II. This is so because in addition to the error shown
10 902 JOURNAL OF MICROELECTROMECHANICAL SYSTEMS, VOL. 16, NO. 4, AUGUST 2007 Fig. 11. Algorithm to calculate the quality factor numerically. in Table II, the error due to the violation of one of the most basic assumptions very small air-gap-to-length ratio used in deriving the Reynolds equation also exists in this case, i.e., the 2-D flow assumption of Reynolds equation is violated. In case 2, the zero-pressure boundary condition is applied on the boundary of the surface extended by (1.3/2)h 0 [10] to account for the fully open boundary effect. The analytical quality factor is then calculated using (4) and (5). Here, the effective dimensions are calculated by L eff = L +1.3h 0 and W eff = W +1.3h 0, whereas the moment of inertia I = kg-m 2 is calculated using the original dimensions of the structure and ω = rad/s from the experimental result. The analytically computed result differs from the experimentally measured value by 14%, which is probably due to the fact that the factor 1.3h 0 does not account for the partially blocked boundaries such as the case in the given structure. Since the extension factor given by Veijola [10] accounts for fully open boundary conditions, Q anal matches well with the corresponding Q num in case 3 of Table II with an error of about 3%. It also validates our computational result. Finally, from the comparative studies of the experimental, numerical, and analytical results, we find that the extension factor 1.3h 0 is indeed useful to capture border effects of fully open boundaries. To the best of our knowledge, this is the first experimental verification of this factor.
11 PANDEY et al.: DYNAMIC CHARACTERISTICS OF SQUEEZE FILMS IN MEMS DEVICES 903 IV. CONCLUSION CFD is suitable for studying the various effects of different boundary conditions such as fully open, fully closed, and partially closed boundaries in oscillating MEMS devices. We also use it to study a problem where the air gap to plate length ratio is large, which makes the conventional Reynolds equation inapplicable. We take a fairly complex torsional MEMS structure with large air gap and complicated boundary conditions, and simulate the effect of squeeze-film damping by progressively refining the model for the boundary conditions. We compare the simulation results with experimentally obtained results (captured by the Q-factor of the device) and show how the numerically computed results get progressively close to the experimental value with refinements in the boundary condition modeling. Starting from the most common and naive model of the trivial pressure boundary condition on the fluid domain enclosed by the oscillating structure, which gives an error of 203% compared to the experimental value, we get to the most realistic model of the partially blocked boundary that gives results within 1% of the experimental value. We then compare the numerical as well as the experimental result with relevant analytical models and show that Veijola s analytical model with an appropriate correction for open border effects comes closest with an error of 14%. APPENDIX The numerical procedure used to calculate the quality factor due to fluid damping is given in Fig. 11. REFERENCES [1] H. Hosaka, K. Itao, and S. Kuroda, Damping characteristics of beam-shaped micro-oscillators, Sens. Actuators A, Phys., vol. 49, no. 9, pp , [2] J. Yang, T. Ono, and M. Esashi, Energy dissipation in submicrometer thick single-crystal silicon cantilevers, J. Microelectromech. Syst., vol. 11, no. 6, pp , Dec [3] J. J. Blech, On isothermal squeeze films, J. Lubr. Technol., vol. 105, pp , [4] F. Pan, J. Kubby, E. Peeters, A. T. Tran, and S. Mukherjee, Squeeze film damping effect on the dynamic response of a MEMS torsion mirror, J. Micromech. Microeng., vol. 8, no. 3, pp , [5] Y.-J. Yang and S. D. Senturia, Numerical simulation of compressible squeezed-film damping, in Proc. Solid-State Sensor and Actuator Workshop, Jun. 1996, pp [6] Y.-J. Yang, M. A. Gretillat, and S. D. Senturia, Effect of air damping on the dynamics of nonuniform deformations of microstructures, in Proc. Transducers, Int. Conf. Solid-State Sens. Actuators, 1997, vol. 2, pp [7] E. S. Hung, Y.-J. Yang, and S. D. Senturia, Low-order models for fast dynamical simulation of MEMS microstructures, in Proc. Transducers, Int. Conf. Solid-State Sens. Actuators, 1997, vol. 2, pp [8] E. S. Hung and S. D. Senturia, Generating efficient dynamical models for microelectromechanical systems from a few finite-element simulation runs, J. Microelectromech. Syst., vol. 8, no. 3, pp , Sep [9] S. Vemuri, G. K. Fedder, and T. Mukherjee, Low-order squeeze film model for simulation of MEMS devices, in Proc. 3rd Int. Conf. Model. and Simul. Microsyst., 2000, pp [10] T. Veijola, A. Pursula, and P. Ráback, Extending the validity of squeezedfilm damper models with elongations of surface dimensions, J. Micromech. Microeng., vol. 15, no. 9, pp , Sep [11] Z. Hao, R. Clark, J. Hammer, M. Whitley, and B. Wingfield, Modeling air-damping effect in a bulk micromachined 2D tilt mirror, Sens. Actuators A, Phys., vol. 102, no. 1, pp , Dec [12] R. B. Darling, C. Hivick, and J. Xu, Compact analytical modelling of squeeze film damping with arbitrary venting conditions using a Green s function approach, Sens. Actuators A, Phys., vol. 70, no. 1/2, pp , [13] ANSYS 9, Finite Element Solver for Multiphysics Problems. [Online]. Available: [14] G. K. Batchelor, An Introduction to Fluid Dynamics. Cambridge, U.K.: Cambridge Univ. Press, [15] O. Pinkus and B. Sternlicht, Theory of Hydrodynamic Lubrication. New York: McGraw-Hill, [16] W.-L. Li, Analytical modelling of ultra-thin gas squeeze film, Nanotechnology, vol. 10, no. 4, pp , [17] T. Veijola, K. Ruokonen, and I. Tittonen, Compact squeezed-film damping model including the open border effects, in Proc. MSM, 2001, pp [18] L. E. Kinsler, A. R. Frey, A. B. Coppens, and J. V. Sanders, Fundamentals of Acoustics. New York: Wiley, [19] T. Veijola, Acoustic impedance elements modeling oscillating gas flow in micro channels, in Proc. MSM, 2001, pp [20] S. S. Rao, Mechanical Vibration. New York: Wesley, [21] G. A. Bird, Molecular Gas Dynamics and the Direct Simulation of Gas Flows. Oxford, U.K.: Oxford Univ. Press, [22] M. Bao, H. Yang, Y. Sun, and P. J. French, Modified Reynolds equation and analytical analysis of squeeze film air damping of perforated structures, J. Micromech. Microeng., vol. 13, no. 6, pp , Nov [23] S. S. Mohite, K. Haneesh, V. R. Sonti, and R. Pratap, Analytical solutions for the stiffness and damping coefficients of squeeze films in MEMS devices with perforated back plates, J. Micromech. Microeng., vol. 15, no. 11, pp , Nov [24] S. V. Patankar, Numerical Heat Transfer and Fluid Flow. New York: Hemisphere, Ashok Kumar Pandey received the B.E. degree in mechanical engineering from Bhilai Institute of Technology, Durg, India, in 2001, and the M.S. degree from the Indian Institute of Science, Bangalore, India, in He is currently working toward the Ph.D. degree in MEMS at the Indian Institute of Science. His current research focus is on the investigation of different energy dissipation mechanisms in micro and nano devices. Rudra Pratap (M 07) received the B.Tech. degree from the Indian Institute of Technology (IIT), Kharagpur, India, in 1985, the M.S. degree in mechanics from the University of Arizona, Tucson, in 1987, and the Ph.D. degree in theoretical and applied mechanics from Cornell University, Ithaca, NY, in From 1993 to 1996, he was with the Sibley School of Mechanical and Aerospace Engineering, Cornell University. In 1996, he was with the Indian Institute of Science, Bangalore, India. From 1997 to 2000, he was an Adjunct Assistant Professor at Cornell University. He is currently an Associate Professor in the Department of Mechanical Engineering, Indian Institute of Science. His research interests are MEMS design, computational mechanics, nonlinear dynamics, structural vibration, and vibroacoustics. Dr. Pratap is a member of the Institute of Smart Structure and System and is on the Editorial Board of the Journal of Computers, Materials and Continua (CMC). Fook Siong Chau received the Ph.D. degree from Nottingham University, Nottingham, U.K. After working as a Research Scientist and Project Engineer in a research laboratory, in 1981, he joined the National University of Singapore, Singapore, where he is currently an Associate Professor in the Department of Mechanical Engineering. He also heads the Applied Mechanics Academic Group. His main research interests are in the modeling, simulation, and characterization of microsystems, particularly MOEMS, and the development and application of optical techniques for nondestructive evaluation of components. Dr. Chau is a member of the Institution of Engineers, Australia, and the Institution of Engineers, Singapore.
ME 237: Mechanics of Microsystems : Lecture. Modeling Squeeze Film Effects in MEMS
ME 237: Mechanics of Microsystems : Lecture Squeeze Film Effects in MEMS Anish Roychowdhury Adviser : Prof Rudra Pratap Department of Mechanical Engineering and Centre for Nano Science and Engineering
More informationAnalytical solution of the modified Reynolds equation for squeeze film damping in perforated MEMS structures
Sensors and Actuators A 135 (2007) 839 848 Analytical solution of the modified Reynolds equation for squeeze film damping in perforated MEMS structures Ashok Kumar Pandey a, Rudra Pratap a,, Fook Siong
More informationChapter 8. Model of the Accelerometer. 8.1 The static model 8.2 The dynamic model 8.3 Sensor System simulation
Chapter 8. Model of the Accelerometer 8.1 The static model 8.2 The dynamic model 8.3 Sensor System simulation 8.2.1 Basic equations 8.2.2 Resonant frequency 8.2.3 Squeeze-film damping 8.2 The dynamic model
More informationDesign and characterization of in-plane MEMS yaw rate sensor
Sādhanā Vol. 34, Part 4, August 2009, pp. 633 642. Printed in India Design and characterization of in-plane MEMS yaw rate sensor K P VENKATESH, NISHAD PATIL, ASHOK KUMAR PANDEY and RUDRA PRATAP CranesSci
More informationCHAPTER 4 DESIGN AND ANALYSIS OF CANTILEVER BEAM ELECTROSTATIC ACTUATORS
61 CHAPTER 4 DESIGN AND ANALYSIS OF CANTILEVER BEAM ELECTROSTATIC ACTUATORS 4.1 INTRODUCTION The analysis of cantilever beams of small dimensions taking into the effect of fringing fields is studied and
More informationAerodynamic behavior of the bridge of a capacitive RF MEMS switch
Aerodynamic behavior of the bridge of a capacitive RF MEMS switch Dragos ISVORANU 1,*, Stefan Sorohan 2, Gabriela CIUPRINA 3 * Corresponding author: Tel.: ++40 (0)213250704; Fax: ++40 (0) 213181007; Email:
More informationEffect of flexural modes on squeeze film damping in MEMS cantilever resonators
IOP PUBLISHING JOURNAL OF MICROMECHANICS AND MICROENGINEERING J. Micromech. Microeng. 17 (7) 475 484 doi:1.188/96-1317/17/1/13 Effect of flexural modes on squeeze film damping in MEMS cantilever resonators
More informationEFFICIENT compact models for the squeeze-film analysis
JOURNAL OF MICROELECTROMECHANICAL SYSTEMS, VOL. 17, NO. 3, JUNE 2008 709 A Compact Squeeze-Film Model Including Inertia, Compressibility, and Rarefaction Effects for Perforated 3-D MEMS Structures Suhas
More informationIMPROVED SENSITIVITY OF RESONANT MASS SENSOR BASED ON MICRO TILTING PLATE AND MICRO CANTILEVER
Michigan Technological University Digital Commons @ Michigan Tech Dissertations, Master's Theses and Master's Reports - Open Dissertations, Master's Theses and Master's Reports 2015 IMPROVED SENSITIVITY
More informationDesign of a MEMS Capacitive Comb-drive Accelerometer
Design of a MEMS Capacitive Comb-drive Accelerometer Tolga Kaya* 1, Behrouz Shiari 2, Kevin Petsch 1 and David Yates 2 1 Central Michigan University, 2 University of Michigan * kaya2t@cmich.edu Abstract:
More information874. The squeeze film effect on micro-electromechanical resonators
874. The squeeze film effect on micro-electromechanical resonators Shih-Chieh Sun 1, Chi-Wei Chung, Chao-Ming Hsu 3, Jao-Hwa Kuang 4 1,, 4 Department of Mechanical and Electromechanical Engineering National
More informationCoding isothermal Reynolds equation in Matlab and comparing with analytical solution for the OOOO boundary condition
Mechanics of Microsystems : Micro/nano Mechanics (NE 211) Course Project Coding isothermal Reynolds equation in Matlab and comparing with analytical solution for the OOOO boundary condition More S. K.
More informationViscous Damping Effect on the CMUT Device in Air
Journal of the Korean Physical Society, Vol. 58, No. 4, April 2011, pp. 747 755 Viscous Damping Effect on the CMUT Device in Air Seung-Mok Lee Micromachined Sensing Laboratory, Ingen MSL Inc., Ayumino
More informationDamping Ratio Analysis of a Silicon Capacitive Micromechanical Accelerometer
Wireless Sensor Network, 017, 9, 178-188 http://www.scirp.org/journal/wsn ISSN Online: 1945-3086 ISSN Print: 1945-3078 Damping Ratio Analysis of a Silicon Capacitive Micromechanical Accelerometer Yuming
More informationEE C245 / ME C218 INTRODUCTION TO MEMS DESIGN FALL 2009 PROBLEM SET #7. Due (at 7 p.m.): Thursday, Dec. 10, 2009, in the EE C245 HW box in 240 Cory.
Issued: Thursday, Nov. 24, 2009 PROBLEM SET #7 Due (at 7 p.m.): Thursday, Dec. 10, 2009, in the EE C245 HW box in 240 Cory. 1. Gyroscopes are inertial sensors that measure rotation rate, which is an extremely
More informationMicrostructure cantilever beam for current measurement
264 South African Journal of Science 105 July/August 2009 Research Articles Microstructure cantilever beam for current measurement HAB Mustafa and MTE Khan* Most microelectromechanical systems (MEMS) sensors
More informationFinite Element Analysis of Piezoelectric Cantilever
Finite Element Analysis of Piezoelectric Cantilever Nitin N More Department of Mechanical Engineering K.L.E S College of Engineering and Technology, Belgaum, Karnataka, India. Abstract- Energy (or power)
More informationDesign and Analysis of dual Axis MEMS Capacitive Accelerometer
International Journal of Electronics Engineering Research. ISSN 0975-6450 Volume 9, Number 5 (2017) pp. 779-790 Research India Publications http://www.ripublication.com Design and Analysis of dual Axis
More informationModeling and simulation of multiport RF switch
Journal of Physics: Conference Series Modeling and simulation of multiport RF switch To cite this article: J Vijay et al 006 J. Phys.: Conf. Ser. 4 715 View the article online for updates and enhancements.
More informationVORTEX LEVITATION. Toshiharu Kagawa 1 and Xin Li 2
VORTEX LEVITATION Toshiharu Kagawa 1 and Xin Li ABSTRACT In this paper, a new pneumatic levitation method, called vortex levitation, is introduced. Vortex levitation can achieve non-contact handling by
More informationIDENTIFICATION OF TEST STRUCTURES FOR REDUCED ORDER MODELING OF THE SQUEEZE FILM DAMPING IN MEMS. Aurelio Somà, Giorgio De Pasquale
Stresa, Italy, 25-27 April 2007 IDENTIFIATION OF TEST STRUTURES FOR REDUED ORDER MODELING OF THE SQUEEZE FILM DAMPING IN MEMS Aurelio Somà, Giorgio De Pasquale Laboratory of Microsystems, Department of
More informationEFFICIENT MULTI-PHYSICS MODELING OF THE DYNAMIC RESPONSE OF RF-MEMS SWITCHES
EFFICIENT MULTI-PHYSICS MODELING OF THE DYNAMIC RESPONSE OF RF-MEMS SWITCHES Jeroen Bielen 1, Jiri Stulemeijer 1 1 EPCOS Netherlands Deepak Ganjoo 2, Dale Ostergaard 2, Stephen Scampoli 2 2 Ansys Inc.
More informationNatural vibration frequency of classic MEMS structures
Natural vibration frequency of classic MEMS structures Zacarias E. Fabrim PGCIMAT, UFRGS, Porto Alegre, RS, Brazil Wang Chong, Manoel Martín Pérez Reimbold DeTec, UNIJUI, Ijuí, RS, Brazil Abstract This
More informationModeling and Design of MEMS Accelerometer to detect vibrations on chest wall
Modeling and Design of MEMS Accelerometer to detect vibrations on chest wall P. Georgia Chris Selwyna 1, J.Samson Isaac 2 1 M.Tech Biomedical Instrumentation, Department of EIE, Karunya University, Coimbatore
More informationEE C245 / ME C218 INTRODUCTION TO MEMS DESIGN FALL 2011 C. Nguyen PROBLEM SET #7. Table 1: Gyroscope Modeling Parameters
Issued: Wednesday, Nov. 23, 2011. PROBLEM SET #7 Due (at 7 p.m.): Thursday, Dec. 8, 2011, in the EE C245 HW box in 240 Cory. 1. Gyroscopes are inertial sensors that measure rotation rate, which is an extremely
More information10 Measurement of Acceleration, Vibration and Shock Transducers
Chapter 10: Acceleration, Vibration and Shock Measurement Dr. Lufti Al-Sharif (Revision 1.0, 25/5/2008) 1. Introduction This chapter examines the measurement of acceleration, vibration and shock. It starts
More informationEE C245 - ME C218 Introduction to MEMS Design Fall Today s Lecture
EE C45 - ME C8 Introduction to MEMS Design Fall 3 Roger Howe and Thara Srinivasan Lecture 9 Energy Methods II Today s Lecture Mechanical structures under driven harmonic motion develop analytical techniques
More informationModel for Gas Damping in Air Gaps of RF MEMS Resonators
Model for Gas Damping in Air Gaps of RF MEMS Resonators T. Veijola, A. Lehtovuori To cite this version: T. Veijola, A. Lehtovuori. Model for Gas Damping in Air Gaps of RF MEMS Resonators. DTIP 2007, Apr
More informationESTIMATION OF ACOUSTIC LOSSES IN THE QUALITY FACTOR OF A MICROMECHANICAL 2D RESONATOR
ESTIMATION OF ACOUSTIC LOSSES IN THE QUALITY FACTOR OF A MICROMECHANICAL D RESONATOR Santhosh Doreswamy Vishwakarma, and Rudra Pratap Dept. of Mechanical Engineering, Indian Institute of Science, Bangalore,
More informationDESIGN AND FABRICATION OF THE MICRO- ACCELEROMETER USING PIEZOELECTRIC THIN FILMS
DESIGN AND FABRICATION OF THE MICRO- ACCELEROMETER USING PIEZOELECTRIC THIN FILMS JYH-CHENG YU and FU-HSIN LAI Department of Mechanical Engineering National Taiwan University of Science and Technology
More informationNONLINEAR CHARACTERISTICS OF THE PILE-SOIL SYSTEM UNDER VERTICAL VIBRATION
IGC 2009, Guntur, INDIA NONLINEAR CHARACTERISTICS OF THE PILE-SOIL SYSTEM UNDER VERTICAL VIBRATION B. Manna Lecturer, Civil Engineering Department, National Institute of Technology, Rourkela 769008, India.
More informationEE C245 - ME C218. Fall 2003
EE C45 - E C8 Introduction to ES Design Fall Roger Howe and Thara Srinivasan ecture 9 Energy ethods II Today s ecture echanical structures under driven harmonic motion develop analytical techniques for
More informationDesign of Electrostatic Actuators for MOEMS Applications
Design of Electrostatic Actuators for MOEMS Applications Dooyoung Hah 1,, Hiroshi Toshiyoshi 1,3, and Ming C. Wu 1 1 Department of Electrical Engineering, University of California, Los Angeles Box 951594,
More informationEE C245 ME C218 Introduction to MEMS Design Fall 2007
EE C245 ME C218 Introduction to MEMS Design Fall 2007 Prof. Clark T.-C. Nguyen Dept. of Electrical Engineering & Computer Sciences University of California at Berkeley Berkeley, CA 94720 Lecture 17: Energy
More informationDesign and Simulation of Comb Drive Capacitive Accelerometer by Using MEMS Intellisuite Design Tool
Design and Simulation of Comb Drive Capacitive Accelerometer by Using MEMS Intellisuite Design Tool Gireesh K C 1, Harisha M 2, Karthick Raj M 3, Karthikkumar M 4, Thenmoli M 5 UG Students, Department
More informationTAILORING THE GEOMETRY OF MICRON SCALE RESONATORS TO OVERCOME VISCOUS DAMPING MARGARITA VILLA
TAILORING THE GEOMETRY OF MICRON SCALE RESONATORS TO OVERCOME VISCOUS DAMPING MARGARITA VILLA Thesis submitted to the Faculty of the Virginia Polytechnic Institute and State University in Partial Fulfillment
More informationSTRUCTURAL CONTROL USING MODIFIED TUNED LIQUID DAMPERS
STRUCTURAL CONTROL USING MODIFIED TUNED LIQUID DAMPERS A. Samanta 1 and P. Banerji 2 1 Research Scholar, Department of Civil Engineering, Indian Institute of Technology Bombay, Mumbai, India, 2 Professor,
More informationAnalysis of Hydrodynamic Journal Bearing Using CFD and FSI Technique
Analysis of Hydrodynamic Journal Bearing Using CFD and FSI Technique Priyanka Tiwari M.E. Student of Government Engineering College Jabalpur, M.P.-India Veerendra Kumar Principal of Government Engineering
More informationDSMC-Based Shear-Stress/Velocity-Slip Boundary Condition for Navier-Stokes Couette-Flow Simulations
DSMC-Based Shear-Stress/Velocity-Slip Boundary Condition for Navier-Stokes Couette-Flow Simulations J. R. Torczynski and M. A. Gallis Engineering Sciences Center, Sandia National Laboratories, P. O. Box
More informationDesign and Analysis of Various Microcantilever Shapes for MEMS Based Sensing
ScieTech 014 Journal of Physics: Conference Series 495 (014) 01045 doi:10.1088/174-6596/495/1/01045 Design and Analysis of Various Microcantilever Shapes for MEMS Based Sensing H. F. Hawari, Y. Wahab,
More informationE05 Resonator Design
POLITECNICO DI MILANO MSC COURSE - MEMS AND MICROSENSORS - 2018/2019 E05 Resonator Design Giorgio Mussi 16/10/2018 In this class we will learn how an in-plane MEMS resonator handles process variabilities,
More informationAnalytical Design of Micro Electro Mechanical Systems (MEMS) based Piezoelectric Accelerometer for high g acceleration
Analytical Design of Micro Electro Mechanical Systems (MEMS) based Piezoelectric Accelerometer for high g acceleration Arti Arora 1, Himanshu Monga 2, Anil Arora 3 Baddi University of Emerging Science
More informationContactless Excitation of MEMS Resonant Sensors by Electromagnetic Driving
Excerpt from the Proceedings of the COMSOL Conference 2009 Milan Contactless Excitation of MEMS Resonant Sensors by Electromagnetic Driving M. Baù *, V. Ferrari, D. Marioli Department of Electronics for
More informationPressure Dependent Quality Factor of Micron Scale Energy Harvesters A.T. Brimmo, M.I. Hassan, A. N. Chatterjee
Pressure Dependent Quality Factor of Micron Scale Energy Harvesters A.T. Brimmo, M.I. Hassan, A. N. Chatterjee Version 1 1 Overview Problem statement Approach Experimentation Computational Model Validation
More informationThe student will experimentally determine the parameters to represent the behavior of a damped oscillatory system of one degree of freedom.
Practice 3 NAME STUDENT ID LAB GROUP PROFESSOR INSTRUCTOR Vibrations of systems of one degree of freedom with damping QUIZ 10% PARTICIPATION & PRESENTATION 5% INVESTIGATION 10% DESIGN PROBLEM 15% CALCULATIONS
More informationTHEORETICAL AND NUMERICAL STUDIES OF THE AIR DAMPING OF MICRO-RESONATORS IN THE NON-CONTINUUM REGIME. Sarne M. Hutcherson
THEORETICAL AND NUMERICAL STUDIES OF THE AIR DAMPING OF MICRO-RESONATORS IN THE NON-CONTINUUM REGIME A Thesis Presented to The Academic Faculty by Sarne M. Hutcherson In Partial Fulfillment of the Requirements
More informationA novel fluid-structure interaction model for lubricating gaps of piston machines
Fluid Structure Interaction V 13 A novel fluid-structure interaction model for lubricating gaps of piston machines M. Pelosi & M. Ivantysynova Department of Agricultural and Biological Engineering and
More informationLIQUID FILM THICKNESS OF OSCILLATING FLOW IN A MICRO TUBE
Proceedings of the ASME/JSME 2011 8th Thermal Engineering Joint Conference AJTEC2011 March 13-17, 2011, Honolulu, Hawaii, USA AJTEC2011-44190 LIQUID FILM THICKNESS OF OSCILLATING FLOW IN A MICRO TUBE Youngbae
More informationMEMS PARALLEL PLATE ACTUATORS: PULL-IN, PULL-OUT AND OTHER TRANSITIONS
MEMS PARALLEL PLATE ACTUATORS: PULL-IN, PULL-OUT AND OTHER TRANSITIONS Subrahmanyam Gorthi, Atanu Mohanty and Anindya Chatterjee* Supercomputer Education and Research Centre, Indian Institute of Science,
More informationExperimental Modal Analysis of a Flat Plate Subjected To Vibration
American Journal of Engineering Research (AJER) 2016 American Journal of Engineering Research (AJER) e-issn: 2320-0847 p-issn : 2320-0936 Volume-5, Issue-6, pp-30-37 www.ajer.org Research Paper Open Access
More informationSound Pressure Generated by a Bubble
Sound Pressure Generated by a Bubble Adrian Secord Dept. of Computer Science University of British Columbia ajsecord@cs.ubc.ca October 22, 2001 This report summarises the analytical expression for the
More informationLecture 20. Measuring Pressure and Temperature (Chapter 9) Measuring Pressure Measuring Temperature MECH 373. Instrumentation and Measurements
MECH 373 Instrumentation and Measurements Lecture 20 Measuring Pressure and Temperature (Chapter 9) Measuring Pressure Measuring Temperature 1 Measuring Acceleration and Vibration Accelerometers using
More informationCHAPTER 5 FIXED GUIDED BEAM ANALYSIS
77 CHAPTER 5 FIXED GUIDED BEAM ANALYSIS 5.1 INTRODUCTION Fixed guided clamped and cantilever beams have been designed and analyzed using ANSYS and their performance were calculated. Maximum deflection
More informationInstitute for Electron Microscopy and Nanoanalysis Graz Centre for Electron Microscopy
Institute for Electron Microscopy and Nanoanalysis Graz Centre for Electron Microscopy Micromechanics Ass.Prof. Priv.-Doz. DI Dr. Harald Plank a,b a Institute of Electron Microscopy and Nanoanalysis, Graz
More informationTheory and Practice of Rotor Dynamics Prof. Rajiv Tiwari Department of Mechanical Engineering Indian Institute of Technology Guwahati
Theory and Practice of Rotor Dynamics Prof. Rajiv Tiwari Department of Mechanical Engineering Indian Institute of Technology Guwahati Module - 7 Instability in rotor systems Lecture - 4 Steam Whirl and
More informationThe Behaviour of Simple Non-Linear Tuned Mass Dampers
ctbuh.org/papers Title: Authors: Subject: Keyword: The Behaviour of Simple Non-Linear Tuned Mass Dampers Barry J. Vickery, University of Western Ontario Jon K. Galsworthy, RWDI Rafik Gerges, HSA & Associates
More informationSimulation of a Micro-Scale Out-of-plane Compliant Mechanism
Simulation of a Micro-Scale Out-of-plane Compliant Mechanism by Arpys Arevalo PhD Candidate in Electrical Engineering Computer, Electrical and Mathematical Sciences and Engineering (CEMSE) King Abdullah
More informationEffect of Liquid Viscosity on Sloshing in A Rectangular Tank
International Journal of Research in Engineering and Science (IJRES) ISSN (Online): 2320-9364, ISSN (Print): 2320-9356 Volume 5 Issue 8 ǁ August. 2017 ǁ PP. 32-39 Effect of Liquid Viscosity on Sloshing
More informationEvaluation of a surface acoustic wave motor with a multi-contact-point slider
Smart Mater. Struct. 7 (1998) 305 311. Printed in the UK PII: S0964-1726(98)91230-7 Evaluation of a surface acoustic wave motor with a multi-contact-point slider Minoru Kuribayashi Kurosawa, Makoto Chiba
More informationNONLINEAR ANALYSIS OF PULL IN VOLTAGE IN MICRO- CANTILEVER BEAM
International Workshop SMART MATERIALS, STRUCTURES & NDT in AEROSPACE Conference NDT in Canada 011-4 November 011, Montreal, Quebec, Canada NONLINEAR ANALYSIS OF PULL IN VOLTAGE IN MICRO- CANTILEVER BEAM
More informationIntroduction to Continuous Systems. Continuous Systems. Strings, Torsional Rods and Beams.
Outline of Continuous Systems. Introduction to Continuous Systems. Continuous Systems. Strings, Torsional Rods and Beams. Vibrations of Flexible Strings. Torsional Vibration of Rods. Bernoulli-Euler Beams.
More informationImplementing a Partitioned Algorithm for Fluid-Structure Interaction of Flexible Flapping Wings within Overture
10 th Symposimum on Overset Composite Grids and Solution Technology, NASA Ames Research Center Moffett Field, California, USA 1 Implementing a Partitioned Algorithm for Fluid-Structure Interaction of Flexible
More informationU-sequence in electrostatic microelectromechanical systems (MEMS)
Proc. R. Soc. A (26) 462, 3435 3464 doi:1.198/rspa.26.1733 Published online 31 May 26 U-sequence in electrostatic microelectromechanical systems (MEMS) BY SUDIPTO K. DE ANDN. R. ALURU* Department of Mechanical
More informationDynamics of structures
Dynamics of structures 2.Vibrations: single degree of freedom system Arnaud Deraemaeker (aderaema@ulb.ac.be) 1 Outline of the chapter *One degree of freedom systems in real life Hypothesis Examples *Response
More informationAn optimum design of a double pendulum in autoparametric resonance for energy harvesting applications
An optimum design of a double pendulum in autoparametric resonance for energy harvesting applications Taizoon Chunawala 1, Maryam Ghandchi-Tehrani 2, Jize Yan 2 1 Birla Institute of Technology and Science-Pilani,
More informationSound Propagation through Media. Nachiketa Tiwari Indian Institute of Technology Kanpur
Sound Propagation through Media Nachiketa Tiwari Indian Institute of Technology Kanpur LECTURE-13 WAVE PROPAGATION IN SOLIDS Longitudinal Vibrations In Thin Plates Unlike 3-D solids, thin plates have surfaces
More informationMagneto-Mechanical Modeling and Simulation of MEMS Sensors Based on Electroactive Polymers
Magneto-Mechanical Modeling and Simulation of MEMS Sensors Based on Electroactive Polymers F.J.O. RODRIGUES, L.M. GONÇALVES, J.H. CORREIA, P.M. MENDES University of Minho, Dept. Industrial Electronics,
More information1859. Forced transverse vibration analysis of a Rayleigh double-beam system with a Pasternak middle layer subjected to compressive axial load
1859. Forced transverse vibration analysis of a Rayleigh double-beam system with a Pasternak middle layer subjected to compressive axial load Nader Mohammadi 1, Mehrdad Nasirshoaibi 2 Department of Mechanical
More informationMODELING AND SIMULATION OF HYDRAULIC ACTUATOR WITH VISCOUS FRICTION
MODELING AND SIMULATION OF HYDRAULIC ACTUATOR WITH VISCOUS FRICTION Jitendra Yadav 1, Dr. Geeta Agnihotri 1 Assistant professor, Mechanical Engineering Department, University of petroleum and energy studies,
More informationMEMS Tuning-Fork Gyroscope Mid-Term Report Amanda Bristow Travis Barton Stephen Nary
MEMS Tuning-Fork Gyroscope Mid-Term Report Amanda Bristow Travis Barton Stephen Nary Abstract MEMS based gyroscopes have gained in popularity for use as rotation rate sensors in commercial products like
More informationSensors & Transducers 2016 by IFSA Publishing, S. L.
Sensors & Transducers, Vol. 96, Issue, January 206, pp. 52-56 Sensors & Transducers 206 by IFSA Publishing, S. L. http://www.sensorsportal.com Collapse Mode Characteristics of Parallel Plate Ultrasonic
More informationStep 1: Mathematical Modeling
083 Mechanical Vibrations Lesson Vibration Analysis Procedure The analysis of a vibrating system usually involves four steps: mathematical modeling derivation of the governing uations solution of the uations
More informationFinite Element Static, Vibration and Impact-Contact Analysis of Micromechanical Systems
Finite Element Static, Vibration and Impact-Contact Analysis of Micromechanical Systems Alexey I. Borovkov Eugeny V. Pereyaslavets Igor A. Artamonov Computational Mechanics Laboratory, St.Petersburg State
More informationPROJECT 1 DYNAMICS OF MACHINES 41514
PROJECT DYNAMICS OF MACHINES 454 Theoretical and Experimental Modal Analysis and Validation of Mathematical Models in Multibody Dynamics Ilmar Ferreira Santos, Professor Dr.-Ing., Dr.Techn., Livre-Docente
More informationFluid Dynamics for Ocean and Environmental Engineering Homework #2 Viscous Flow
OCEN 678-600 Fluid Dynamics for Ocean and Environmental Engineering Homework #2 Viscous Flow Date distributed : 9.18.2005 Date due : 9.29.2005 at 5:00 pm Return your solution either in class or in my mail
More informationAnswers to questions in each section should be tied together and handed in separately.
EGT0 ENGINEERING TRIPOS PART IA Wednesday 4 June 014 9 to 1 Paper 1 MECHANICAL ENGINEERING Answer all questions. The approximate number of marks allocated to each part of a question is indicated in the
More informationThe basic principle to be used in mechanical systems to derive a mathematical model is Newton s law,
Chapter. DYNAMIC MODELING Understanding the nature of the process to be controlled is a central issue for a control engineer. Thus the engineer must construct a model of the process with whatever information
More informationDesign and Analysis of a Triple Axis Thermal Accelerometer
Design and Analysis of a Triple Axis Thermal Accelerometer DINH Xuan Thien a and OGAMI Yoshifumi b Ritsumeikan University, Nojihigashi, Kusatsu, Shiga 525 8577 Japan a thien@cfd.ritsumei.ac.jp, b y_ogami@cfd.ritsumei.ac.jp
More information1820. Selection of torsional vibration damper based on the results of simulation
8. Selection of torsional vibration damper based on the results of simulation Tomasz Matyja, Bogusław Łazarz Silesian University of Technology, Faculty of Transport, Gliwice, Poland Corresponding author
More informationOptimal Design of Capacitive Micro Cantilever Beam Accelerometer
Optimal Design of Capacitive Micro Cantilever Beam Accelerometer Othman Sidek Collaborative Microelectronic Design Excellence Centre Universiti Sains Malaysia, Engineering Campus, 14300 Nibong Tebal Sri
More informationChapter 2 Nonlinear Modeling of Squeeze-Film Phenomena
Chapter 2 Nonlinear Modeling of Squeeze-Film Phenomena in Microbeam MEMS Reza N. Jazar Abstract Oscillating microplates attached to microbeams is the main part of many microresonators and micro-electro-mechanical
More informationAir Damping of Oscillating MEMS Structures: Modeling and Comparison with Experiment
Air Damping of Oscillating MEMS Structures: Modeling and Comparison with Experiment Sergey Gorelick *1, James R. Dekker 1, Mikko Leivo 1, and Uula Kantojärvi 1 1 VTT Technical Research Centre of Finland,
More informationResearch Article. Slip flow and heat transfer through a rarefied nitrogen gas between two coaxial cylinders
Available online wwwjocprcom Journal of Chemical and Pharmaceutical Research, 216, 8(8):495-51 Research Article ISSN : 975-7384 CODEN(USA) : JCPRC5 Slip flow and heat transfer through a rarefied nitrogen
More informationChapter a. Spring constant, k : The change in the force per unit length change of the spring. b. Coefficient of subgrade reaction, k:
Principles of Soil Dynamics 3rd Edition Das SOLUTIONS MANUAL Full clear download (no formatting errors) at: https://testbankreal.com/download/principles-soil-dynamics-3rd-editiondas-solutions-manual/ Chapter
More informationEE C245 ME C218 Introduction to MEMS Design Fall 2012
EE C245 ME C218 Introduction to MEMS Design Fall 2012 Prof. Clark T.-C. Nguyen Dept. of Electrical Engineering & Computer Sciences University of California at Berkeley Berkeley, CA 94720 Lecture EE C245:
More informationJOURNAL OF MICROELECTROMECHANICAL SYSTEMS, VOL. 13, NO. 5, OCTOBER Sudipto K. De and N. R. Aluru, Member, IEEE, Associate Member, ASME
JOURNAL OF MICROELECTROMECHANICAL SYSTEMS, VOL. 13, NO. 5, OCTOBER 2004 737 Full-Lagrangian Schemes for Dynamic Analysis of Electrostatic MEMS Sudipto K. De N. R. Aluru, Member, IEEE, Associate Member,
More informationIntroduction to structural dynamics
Introduction to structural dynamics p n m n u n p n-1 p 3... m n-1 m 3... u n-1 u 3 k 1 c 1 u 1 u 2 k 2 m p 1 1 c 2 m2 p 2 k n c n m n u n p n m 2 p 2 u 2 m 1 p 1 u 1 Static vs dynamic analysis Static
More informationEE 5344 Introduction to MEMS CHAPTER 6 Mechanical Sensors. 1. Position Displacement x, θ 2. Velocity, speed Kinematic
I. Mechanical Measurands: 1. Classification of main types: EE 5344 Introduction MEMS CHAPTER 6 Mechanical Sensors 1. Position Displacement x, θ. Velocity, speed Kinematic dx dθ v =, = ω 3. Acceleration
More informationAbstract. 1 Introduction
In R. A. Adey et al., eds., Simulation and Design of Microsystems and Microstructures (Proceedings of the 1st International Conference on Simulation and Design of Microsystems and Microstructures), Computational
More informationA new cantilever beam-rigid-body MEMS gyroscope: mathematical model and linear dynamics
Proceedings of the International Conference on Mechanical Engineering and Mechatronics Toronto, Ontario, Canada, August 8-10 2013 Paper No. XXX (The number assigned by the OpenConf System) A new cantilever
More informationWORK SHEET FOR MEP311
EXPERIMENT II-1A STUDY OF PRESSURE DISTRIBUTIONS IN LUBRICATING OIL FILMS USING MICHELL TILTING PAD APPARATUS OBJECTIVE To study generation of pressure profile along and across the thick fluid film (converging,
More informationPRESSURE AND VELOCITY AMPLITUDES OF THE INCOMPRESSIBLE FLUID IN CONCENTRIC ANNULAR PASSAGE WITH OSCILLATORY BOUNDARY: TURBULENT FLOW
Journal of Engineering Science and Technology Vol. 9, No. 2 (2014) 220-232 School of Engineering, Taylor s University PRESSURE AND VELOCITY AMPLITUDES OF THE INCOMPRESSIBLE FLUID IN CONCENTRIC ANNULAR
More informationA MICROBEARING GAS FLOW WITH DIFFERENT WALLS TEMPERATURES
Mili}ev, S. S., et al.: A Microbearing Gas Flow with Different Walls Temperatures THERMAL SCIENCE, Year 01, Vol. 16, No. 1, pp. 119-13 119 A MICROBEARING GAS FLOW WITH DIFFERENT WALLS TEMPERATURES by Snežana
More informationWEEKS 8-9 Dynamics of Machinery
WEEKS 8-9 Dynamics of Machinery References Theory of Machines and Mechanisms, J.J.Uicker, G.R.Pennock ve J.E. Shigley, 2011 Mechanical Vibrations, Singiresu S. Rao, 2010 Mechanical Vibrations: Theory and
More informationEE C245 ME C218 Introduction to MEMS Design
EE C245 ME C218 Introduction to MEMS Design Fall 2007 Prof. Clark T.-C. Nguyen Dept. of Electrical Engineering & Computer Sciences University of California at Berkeley Berkeley, CA 94720 Lecture 23: Electrical
More informationMicroelectromechanical Systems (MEMs) Applications Fluids
ROCHESTER INSTITUTE OF TEHNOLOGY MICROELECTRONIC ENGINEERING Microelectromechanical Systems (MEMs) Applications Fluids Dr. Lynn Fuller Webpage: http://people.rit.edu/lffeee 82 Lomb Memorial Drive Rochester,
More informationDESIGN AND OPTIMIZATION OF BULK MICROMACHINED ACCELEROMETER FOR SPACE APPLICATIONS
INTERNATIONAL JOURNAL ON SMART SENSING AND INTELLIGENT SYSTEMS, VOL. 1, NO. 4, DECEMBER 008 DESIGN AND OPTIMIZATION OF BULK MICROMACHINED ACCELEROMETER FOR SPACE APPLICATIONS Thampi Paul 1, Jaspreet Singh,
More informationA Vertical Electrostatic Actuator with Extended Digital Range via Tailored Topology
A Vertical Electrostatic Actuator with Extended Digital Range via Tailored Topology Yanhang Zhang and Martin L. Dunn Department of Mechanical Engineering University of Colorado at Boulder Boulder, CO 80309
More informationEE C245 - ME C218 Introduction to MEMS Design Fall Today s Lecture
EE C45 - ME C18 Introduction to MEMS Design Fall 003 Roger Howe and Thara Srinivasan Lecture 11 Electrostatic Actuators II Today s Lecture Linear (vs. displacement) electrostatic actuation: vary overlap
More informationEE C245 ME C218 Introduction to MEMS Design
EE C245 ME C218 Introduction to MEMS Design Fall 2007 Prof. Clark T.-C. Nguyen Dept. of Electrical Engineering & Computer Sciences University of California at Berkeley Berkeley, CA 94720 Lecture 16: Energy
More information