DETC THE EFFECT OF HIGH ORDER NON-LINEARITIES ON SUB-HARMONIC EXCITATION WITH PARALLEL PLATE CAPACITIVE ACTUATORS

Size: px
Start display at page:

Download "DETC THE EFFECT OF HIGH ORDER NON-LINEARITIES ON SUB-HARMONIC EXCITATION WITH PARALLEL PLATE CAPACITIVE ACTUATORS"

Transcription

1 Proceeings of the ASME 27 International Design Engineering Technical Conferences & Computers an Information in Engineering Conference IDETC/CIE 27 September 4-7, 27, Las Vegas, Nevaa, USA DETC THE EFFECT OF HIGH ORDER NON-LINEARITIES ON SUB-HARMONIC EXCITATION WITH PARALLEL PLATE CAPACITIVE ACTUATORS Alexaner A. Trusov MicroSystems Laboratory Dept. of Mechanical an Aerospace Engineering University of California, Irvine Irvine, CA Anrei M. Shkel MicroSystems Laboratory Dept. of Mechanical an Aerospace Engineering University of California, Irvine Irvine, CA ABSTRACT Electrostatic actuation of motion is commonly use in resonant MEMS. Drive signal fee-through is unesirable as it masks the etection signal. This paper reports analysis an emonstration of a resonant motion excitation scheme, which uses a combination of a DC bias with a sinusoial AC voltage, frequency of which is twice of the mechanical resonant frequency. This configuration is experimentally emonstrate to excite a gyroscope with a parallel plate rive capacitor into high amplitue perioic vibrations. The fee-through has a frequency higher than the main motional harmonic an thus can be eliminate by simple low pass filtering. Full nonlinear ynamics along with several higher-orer approximations are consiere. Analysis of the effect of the approximation orer on the frequency response accuracy shows that the complete nonlinear equation shoul be use for moeling of high amplitue actuation. Two istinct types of frequency responses are examine as functions of riving AC an DC voltages an amping. INTRODUCTION Several important classes of MEMS evices, such as resonators [], gyroscopes [2], an chemical sensors [3], rely on resonance phenomenon in their operation. In these evices, resonant motion nees to be actuate, etecte, an controlle. Capacitive phenomena are commonly use for transuction in vi- Aress all corresponence to this author. bratory MEMS evices ue to the ease of fabrication, low sensitivity to temperature changes, an other practical avantages [4], [5]. Conventional electrostatic actuation schemes use a combination of DC an AC voltages to excite mechanical motion at the resonant frequency of the structure. This is typically achieve by matching the frequency of the rive AC voltage to the resonant frequency of the evice. In this case, the amplitue of motion is ictate by the amplitues of the riving voltages an the quality factor of the mechanical resonator. For this type of actuation, the electrical rive fee-through from rive to sense electroes cause by parasitic capacitance has the same frequency content as the mechanical motion. Thus, the motional signal is maske by the fee-through, an intricate etection schemes, such as electromechanical amplitue moulation [6], are require to extract an measure the mechanical response. These factors increase the overall system complexity an limit the performance. Nonlinearity of parallel plate capacitors presents esign opportunities for introucing new electrostatic actuation schemes. Issues of large amplitue parallel plate actuation using a conventional combination of DC an AC voltage were extensively stuie, for example in [7 9]. Parametric excitation phenomenon in MEMS was reporte in []. In its simplest form, parametric excitation using electrostatic parallel plate force is escribe by a classical linear time variant LTV) Mathieu equation. Dynamics of a clampeclampe beam uner parametric excitation was stuie in [] Copyright c 27 by ASME

2 using small amplitue assumption. The effect of an aitional thir orer term was stuie in [2]. An interesting parametrically enhance harmonic actuation scheme using two separate sets of capacitors an AC voltages was propose in [3]. Most of the relevant literature stuies parallel plate parametric effects using approximate equations of motion. The focus of this paper is actuation of high amplitue vibrations in capacitive MEMS using a combination of a DC bias with a single AC sinusoial voltage, frequency of which is twice of the mechanical resonant frequency. In Section we introuce a general electro-mechanical moel of a MEMS resonator with parallel plate actuation an erive the equations of motion. The phenomenon of high amplitue parametric excitation is emonstrate experimentally in Section 2 using a prototype of a capacitive MEMS gyroscope. In Section 3 we efine two ifferent norms of solutions use for frequency response analysis for the nonlinear system. Section 4 presents a comparative numerical stuy of system ynamics using the complete nonlinear moel an its several finite orer approximations. Section 5 summarize the obtaine results. PARAMETRIC EXCITATION Figure shows a general schematic of a capacitive microresonator, a basic element of various micro-sensors. The electromechanical iagram inclues the mechanical resonator an the electrostatic rive an sense electroes. Mass m of the resonator is suspene by a spring k an is constraine to move only along the horizontal x-axis. The variable sense capacitance is efine as C s x), an the rive capacitance as C x), where x is the isplacement. Typically in MEMS evices, the rive an sense terminals are not completely insulate, but are electrically couple by stray parasitic capacitors an resistors [6,4,5]. In this paper we assume, without loss of generality, that the parasitic circuit consists of a single lumpe capacitor C p. An AC riving voltage V t) = v sinω t) is applie to the rive electroe, an a DC bias V c is applie to the mass voltage values are reference with respect to a common groun). The sense capacitor C s is forme between the mobile mass an the fixe sense electroe. The sense electroe is connecte to the inverting input of an operational amplifier which is configure as a trans-impeance amplifier, [6]. Due to motion the sense capacitance C s x) changes, causing a flow of motional current I s = C sv s ) t, where V s is the sensing voltage across the sense capacitor. The total pick-up current It) = I s t)+i p t) consists of both the motional an the parasitic fee-through currents an is converte to the final output voltage V t) with trans-impeance gain R. Parasitic current is inuce by the rive voltage V an therefore has the same frequency ω. The total sensing voltage is the DC bias V c. Accoring to the laws of electrostatics, the total pick-up current that flows through the fee-back resistor of the trans-impeance amplifier is It) = t [V t)c p + V c C s t)]. Drive ~ Vt) Amplitue, B Vc Figure. Drive Electroe C Parasitic Circuit Mechanical resonance x Moveable Mass Vc Cp Sense Electroe Cs Ist) Sense Itt)=Is+Ip R - V=-RI Feethrough Ipt) Pick-up Voltage motional + parasitic Frequency + sub harmonic resonance no resonance, fee through Transimpeance Amp. Schematic of a capacitive MEMS resonator with parallel plate sub-harmonic excitation. Since the parasitic capacitance C p an the voltage V c are constant, It) = V c C s t) t The motion of the resonator is governe by ẍ + ω n Q ẋ + ω2 nx = 2m V t) +C p. ) t C x) V c +V t)) 2, 2) x k where ω n = m is the unampe natural frequency an Q is the quality factor. Let us enote the initial gap between plates in the parallel plate riving capacitor at rest by g, the overlap length of iniviual parallel plate pairs by L, height of the plates i.e. structural layer thickness) by y, an permittivity of the meia by ε. The total overlap area in the rive capacitor is given by A = NLy, where N is a number of parallel plate pairs in the rive capacitor. Then, the variable rive capacitance is C x) = εa g x = εa g x/g. 3) 2 Copyright c 27 by ASME

3 We introuce the nominal rive capacitance C n = εa/g an gap-normalize imensionless isplacement χ = x/g <. From Eqn. 3), using Taylor series expansion with respect to χ, the rive capacitance is C = C n χ = C n n= χ n. 4) The electrostatic force prouce by the rive capacitor is proportional to the graient of the capacitance C x = ) Cn = C n x χ g χ) 2 = C n g n= n + )χ n. 5) Finally, from Equations 2) an 5), normalize isplacement χ of the resonator is governe by χt) + ω n Q χt) + ω2 nχt) = C n 2mg 2 V c +V ) 2 χ) 2, 6) or equivalently, by n= χt) + ω n Q χt) + ω2 nχt) = C n 2mg 2 { Vc v2 ) + 2V cv sinω t) 2 v2 cos2ω t) n + )χ n 7) }. This equation has a nonlinear an time epenent right han sie. The nonlinearity is ue to the parallel plate capacitance graient, Fig. 2; the voltage-square term causes the time epenency of the excitation. Let us consier several special cases. A linear time invariant LTI) approximate moel can be obtaine by neglecting the linear time-epenent an the nonlinear terms on the right han sie. The resulting moel is ientical to the case of classical lateral comb actuation. For this moel, the resonant motion can be achieve in the following two cases: riving AC voltage v t) is applie at the resonance frequency or at its half. In the first case, ω ω n an the forcing term C n V mg 2 c v sinω t) provies resonant excitation. In the secon case, ω 2 ω n an C n v 2 4mg 2 cos2ω t) provies resonant excitation. The rawback of ω ω n is that motion an parasitics occur at the same frequency. The main isavantage of ω ω n /2 is the low actuation gain, cause by the 2 v2 coefficient. The next level of approximation can be obtaine in the following way. First, we assume V c V an neglect terms proportional to v 2. Secon, we linearize the capacitance graient Normalize capacitive graient C χ) /χ, log Figure 2. Parallel plate rive complete an approximate moels complete / χ) 2 orer 5 approx. orer 3 approx. orer 2 approx. linear approx. constant approx. 2 Gap normalize isplacement χ=x/g, log orer approximations. Nonlinear capacitive graient of parallel plates an its finite C n g + 2χ). Finally, Eqn. 7) is simpli- using Eqn. 5) as C x fie to which is equivalent to χt) + ω n Q χt) + ω2 nχt) = C n 2mg 2 { V 2 c + 2V c v sinω t) } + 2χ), 8) χt) + ω n Q χt) + ω 2 n V ) cc n mg 2 [V c + 2v sinω t)] χt) = = V cc n 2mg 2 [V c + 2v sinω t)]. 9) This equation obtaine by neglecting v 2 an linearization with respect to χ) is very similar to the classical Mathieu equation [], which escribes zero-input parametric instability, cause by harmonic moulation of the effective stiffness. However, Eqn. 9) is not homogenous an is simultaneously force with a static force an a sinusoial force at ω frequency. This equation, common in MEMS literature on parallel plate actuation [], qualitatively explains two phenomena inherent to parallel-plate actuation: active resonant frequency tuning an parametric resonance. The effective natural frequency for a parallel-plate actuate motion is given by ω 2 n V cc n ) [V mg 2 c + 2v sinω t)]. The stiffness reuces ) quaratically with applie DC bias as k ω 2 n C n V 2 mg 2 c. This effect is known as frequency tuning an the corresponing mathematical moel is known as negative 3 Copyright c 27 by ASME

4 CDIP-4 Boning pas 89 Frequency Tuning EFAB DMG) Biase Resonant Frequency, Hz experimantal parabolic fit Unbiase resonant frequency is Hz Figure 3. EFAB ie Package prototype of istribute mass gyroscope DMG) fabricate in EFAB process. electrostatic spring. In Section 2 frequency tuning of a MEMS gyroscope prototype is experimentally characterize. Classical parametric resonance is a property of a linear timevarying ynamic system to evelop unstable motion from nontrivial initial conitions uner harmonic moulation of the effective stiffness. Parametric resonance for both non-ampe an ampe systems governe by Mathieu equation is a well stuie problem in both cases motion with infinitely growing amplitue occurs when the stiffness moulation frequency satisfies the conition ω = 2 k ω n, k =,2,...,. The case of k =, ω = 2ω n, is of the greatest practical interest, because it is the only k for which the frequency of the fee-through is above the frequency of motion. Feasibility of this scheme of excitation is experimentally stuie in Section 2. Numerical analysis of ifferent orer approximations of Eqn. 7) is presente in Section 4. 2 DEMONSTRATION OF SUB-HARMONIC EXCITA- TION The simplifie LTV approximation given by Eqn. 9) allows to preict phenomenon of parametric resonance in parallel plate actuate MEMS. Specifically, we can expect the system to evelop significant amplitues of motion, when riving AC voltage is applie at ouble of the resonant frequency taking into account frequency own-tuning ue to the DC voltage). The motion occurs at the tune resonant frequency, which is half of the riving AC voltage frequency. In this section we experimentally characterize this nonlinear actuation metho. Figure 4. Characterization of parallel plate frequency tuning in the evice DC Bias Voltage, V 2. Test-Be A prototype of a Distribute Mass Gyroscope DMG) [5] was use for experimental emonstration of parametric actuation. The evice was fabricate using EFAB TM [7], a commercial process from Microfabrica. The core of the EFAB process is a sequential electroeposition of multiple structural layers. Square ies 5 5) mm were forme on.5 mm thick alumina Al 2 O 3 ) substrate. The first µm thick Ni-Co structural layer was eposite using electroplating an patterne. Then, blanket Cu layer was electroforme an planarize to the thickness of the structural layer. The cycle of these two steps was repeate 2 times, after which the sacrificial Cu layer was chemically etche away. As a result, a 2 µm high 3D metal structures on a nonconuctive substrate were forme. For the fabricate gyroscope, the minimal line-with was 2 µm an minimal spacing was 3 µm. Figure 3 shows a package an wire-bone evice. The rive moe of the gyroscope was use in this paper to stuy high amplitue actuation base on nonlinear parametric resonance phenomenon. 2.2 Structural Characterization In orer to characterize frequency tuning an the unbiase natural frequency of the evice, resonance was visually etecte at ifferent values of V c. Figure 4 shows the collecte ata along with the curve fitting results. Structural properties of the evice were also characterize in air an in vacuum using electrostatic metho of parasitics-free half-frequency sweep [5]. Figure 5 shows the experimental ata an the curve fitting results. The resonant frequency was 878 Hz for a DC bias of 6 V; the quality factor was 275 in air an 86 in 2 mtorr vacuum. 4 Copyright c 27 by ASME

5 8 Structural Characterization of EFAB DMG in Air fixe plate mobile plate Amplitue Response, Linear Scale fn =878.4 Hz Q =274.7 experimental analytical fit 25 m 2 m Amplitue Response, Linear Scale Drive frequency, Hz a) In air Structural Characterization of EFAB DMG in Vacuum fn =878.8 Hz Q =859.3 experimental analytical fit Drive Frequency, Hz Figure 5. b) In 2 mtorr vacuum Structural characterization of the evice using electrostatic actuation an capacitive etection. 2.3 High Amplitue Parametric Resonance Nonlinear parametric excitation of motion at resonant frequency by a twice higher frequency AC voltage is experimentally characterize in this section. This actuation metho provies high amplitue of motion an robustly separates motion from the parasitics in frequency omain. The experimental setup consiste of a package gyroscope assemble on a breaboar together with a trans-impeance amplifier. The teste moe of vibration was characterize utilizing parallel plate actuation an sensing, as shown in Fig.. The ca- Figure 6. 3 m mobile plate at rest Superimpose optical photographs of parallel plate evice at rest an at resonant motion parametrically excite by a twice higher frequency AC voltage. The blurre area shows amplitue of motion 8% of the gap). pacitive gap between 2 µm wie parallel plates was 3 µm. The DC bias voltage was applie to the boy of the resonator. The riving AC voltage was impose onto the anchore parallel plate structure. Sensing of the motion was one using anchore parallel plate structure, which was virtually groune by means of trans-impeance amplifier with a MΩ feeback resistor. The sense capacitor between fixe an moveable parallel plates ha a fixe DC voltage ifference across it. Velocity epenent motional current flows through the feeback resistor of the transimpeance amplifier. Also, capacitive coupling between rive an sense electroe contributes to the output current. The total pick-up current was converte to voltage an fe into a signal analyzer HP 35665A for spectral analysis an ata capturing. During the experiments, the ω frequency of the riving AC voltage V t) was manually swept in the vicinity of 2ω n. The amount of motion was etecte visually using a microscope, Fig. 6. At the same time, the Fast Fourier Transform FFT) of the pick-up signal was continuously performe, upate, isplaye an store by the ynamic signal analyzer. We compare results of the FFT at resonance an off-resonance in our experimental setup this means that motional amplitues were smaller than optically resolvable.25µm). From this comparison, shown in Fig. 7, we etermine the frequency content of the fee-through signals an istinct features of the pick-up signal in the presence of parametrically excite motion. 5 Copyright c 27 by ASME

6 Amplitue, BVrms Excitation of resonant motion by a superharmonic AC voltage motion 2n motion+ st rive at resonance no resonance 4th motion+ 2n rive 6th motion+ 3r motion 3n rive 5th motion 3 FREQUENCY RESPONSE NORMS Small amplitue parametrically excite motion can be moele by classical Mathieu equation []. The focus of this paper, however, is the steay-state large amplitue vibrations excite using parallel plate capacitors. Parametric excitation of large amplitue vibrations was experimentally emonstrate in Section 2. In this section we moel an analyze the ynamics using the complete Eqn. 7) along with several orers of approximations. These moels of the ynamics can be uniformly escribe by 8 Figure Frequency, Hz Experimental pick-up signal PSD for parametrically excite resonator. The riving AC voltage is provie at ouble of the resonant frequency. The strong frequency separation between first main) motional component an parasitics is achieve. We registere the large amplitue of motion when the riving AC voltage was applie at frequency close to twice of the DC tune resonant frequency i.e. ω 2ω n, f 2 85 =.7 khz). The amplitue of motion was several times higher than in the case of conventional LTI excitation i.e. ω ω n ) using the same combination of DC an AC voltages. Due to the capacitive parasitic coupling an non-ieality of the AC voltage generator, the fee-through of the rive occurs at all multiple harmonics of the riving AC frequency i.e. at n.7 =.7,3.4,5.,... khz). In other wors, the rive voltage leakage is present at all even orer multiples of the resonant frequency. However, the motion mainly occurs at the resonant frequency. Thus, the main vibrational moe is sub-harmonic with respect to the riving AC voltage. In aition to the main moe at the resonant frequency, the motion also occurs at the resonant frequency an all of its multiples i.e. at n.85 =.85,.7,2.55,3.4,4.25,5.,... khz). To conclue, nonlinear parametric excitation of resonant motion by a higher frequency voltage was emonstrate using a parallel plate actuate evice. This non-lti actuation approach experimentally exhibits several important avantages over the conventional LTI riving schemes. The avantages inclue actuation of large amplitues of motion 8% of the gap) using relatively small voltages, an robust frequency separation between the motion an the rive AC voltage fee-through. Parametric excitation of perioic resonant motion by a higher frequency AC voltage can potentially be use to provie resonant evices with wier banwith an higher gain while simplifying the etection of motion. N n= χt) + ω n Q χt) + ω2 nχt) = C n 2mg 2 { Vc v2 ) + 2V cv sinω t) 2 v2 cos2ω t) n + )χ n ) }, where N efines the orer of the Taylor sum approximation for the nonlinear capacitive graient on the right han sie. Orer N = correspons to the often use LTV moel, while infinite summation correspons to the complete nonlinear moel of C χ)/ χ = / χ) 2. Equation ) was solve numerically in the complete nonlinear form along with N =,2,3,5 orer approximations using integration routine oe45 in MATLAB. Base on the sensitivity analysis, the time increment was set to 2 of the resonant vibrations perio to ensure precision of the numerical solutions. The total solution time interval was s. The output of the coe was a multi-imensional array of time omain solutions, corresponing to ifferent orers N an values of the system parameters. Frequency responses were calculate from the steay-state response parts of the time omain solutions by a separate MAT- LAB coe. It is important to explicitly efine the notion of the frequency response, since for the non-lti system escribe by Eqn. ) the motion occurs at an infinite number of frequencies. We efine two alternative norms of the χt) solutions for a given rive AC frequency ω : χt) max an χt) ω. The total amplitue of motion norm χt) max is efine as the maximum span of χt) on a certain time interval. This norm etermines if the solution is boune an physically feasible. This approach to frequency response analysis is typically use to stuy instability conitions of parametrically excite micro-resonators, for example []. Operation of many resonant MEMS sensors, however, is base on a certain moe of vibration. For instance, in vibratory gyroscopes the input angular rate stimulus is amplituemoulate an subsequently etecte at the frequency of the resonant vibrational moe. If parametric phenomenon is use to actuate such evices, the amplitue of only the main resonant motion harmonic efines the scale factor an other performance 6 Copyright c 27 by ASME

7 characteristics of the evice. Therefore, for micro-evices base on sinusoial motion an Fourier signal processing, a ifferent efinition of norm shoul be use to stuy frequency responses. This norm is the amplitue of the sinusoial component of the motion at the resonant frequency. In other wors, χt) ω is the norm of the orthogonal projection of χt) onto the Fourier basis function e i ω t 2, i.e. χ ω = R t t χt)e i ω 2 t) t R t t e iωt). ) t To summarize, two ifferent types of frequency response analyses were numerically stuie for a parametrically excite micro-evice. It was conclue, that rive frequency ω epenency of χt) max escribes the instability regions, while χ ω efines performance characteristics of sinusoial motion base resonant micro-evices. 4 MODELING OF NON-LINEARITIES In this section we iscuss the results of ynamics moeling. The moel of the resonator is base on the evice, which was characterize in Section 2; the evice has the following nominal parameters: f n = 85 Hz, Q = 27, effective m = 9.3e-7 kg, C n =.5 pf, V c = 2 V, v = 5 V pk. 4. Effect of Approximation Orer Figure 8a) shows χ max frequency responses of the resonator, obtaine using N =, 2, 3 orer approximations of Eqn. ) an a complete nonlinear moel. Mathieu equation type LTV moel shows classical unboune parametric resonance, which is ifferent from experimentally observe high amplitue vibration; the frequency omain instability region given by this moel is also ifferent from the actual. Secon orer moel shown in Fig. 8a) is also quantitatively ifferent from the complete moel, however the preicte resonant region is very close to the actual. Thir orer moel is close to the actual, however the amplitue of the total χ max response is overestimate. In general, finite orer N approximations of Eqn. ) result in overestimate amplitues of motion; amount of the overestimation reuces with higher orer approximation orer N. Figure 8b) shows a comparison between χ max an χt) ω frequency responses for a complete moel an a thir orer approximation. The frequency omain resonant regions coincie, an the profiles of response curves are similar. Most of the motion is prouce by a harmonic at χt) ω frequency. The harmonic frequency responses o not exhibit a character- Normalize amplitue of motion, χ = x /g Normalize amplitue of motion, χ = x /g orer 3 moel Parametric rive ifferent moels orer 2 moel LTV moel complete moel Normalize frequency of V ac, f Figure 8. a) Maximum amplitue of motion χt) max. Parametric rive effective frequency responses Nonlin3 total Nonlin3 resonant Complete total Complete resonant Normalize frequency of V ac, f b) Comparison of χt) max an χt) ω frequency responses. Frequency responses of parametrically excite vibrations efine by two ifferent norms of solutions. istic peak, which is observe in the total amplitue frequency response. 4.2 Effect of AC Drive Voltage Amplitue v Figure 9 shows χt) ω frequency responses of the resonator for ifferent AC riving voltages v = 2,3,4,5,8,2 V base on the complete moel. Depening on the amplitue of the riving AC voltage, there are three regimes of the ynamics. For small AC amplitues in this case, approximately below V ), there is no parametric resonance. For a range of voltages in this case, approximately between an 6 V ), there is a well pronounce high amplitue pass-ban with a sharp transi- 7 Copyright c 27 by ASME

8 Normalize amplitue of motion, χ = x /g Parametric rive effect of Vac 8V 5V 4V 2V Normalize frequency of V ac, f Figure 9. Resonant harmonic frequency response χt) ω for parametric excitation. Effect of the AC riving voltage amplitue. 3V 2V Normalize amplitue of motion, χ = x /g V Parametric rive effect of Vc V c 3 V 25 V 2 V 5 V V Normalize frequency of V ac, f Figure. Resonant harmonic frequency response χt) ω for parametric excitation. Effect of the DC bias. tion on the left sie an a smooth slope on the right sie. Finally, for voltages above a certain bifurcation threshol, the pass-ban breaks into two ifferent regions. The left part of the pass-ban shows instability, while the right part allows for stable high amplitue motion. This behavior is specific to the complete nonlinear moel of parallel plate parametric excitation. 4.3 Effect of DC Bias Voltage V c Figure shows the effect of the DC bias on the χt) ω frequency response. As expecte from the negative electrostatic spring effect, the pass-ban shifts to the left with increasing V c. Aitionally, the pass-ban becomes wier an shorter. This is ue to the reuction of the effective parallel plate gap, cause by the static isplacement of the mass in response to the DC force. 4.4 Effect of Quality Factor Q In classical Mathieu equation, amping is known to efine parametric instability regions without saturating the amplitue of motion []. Figure shows the effect of Q-factor, or, equivalently the amping, on the χt) ω frequency response of a complete nonlinear parallel plate moel. Frequency response has a low sensitivity to the quality factor in the wie range of values 2 4. This relaxes packaging requirements for parametrically actuate evices []. Quality factors in the range of 2 3 are typical of vibratory gyroscopes in air an appear favorable for the parametric excitation. Normalize amplitue of motion, χ = x /g Parametric rive effect of Q Q=256 Q=24 Q=248 Q= Normalize frequency of V ac, f Figure. Resonant harmonic frequency response χt) ω for parametric excitation. Effect of the quality factor Q. 5 CONCLUSIONS In this paper we stuie high amplitue actuation of capacitive vibratory evices using parametric resonance. The actuation scheme uses a simple combination of a DC bias with a sinusoial AC voltage, frequency of which is twice of the resonator s resonant frequency. Experiments emonstrate excitation of high amplitue motion in the parallel plate rive moe of a gyroscope. Since motion mostly occurs at half of the riving AC voltage frequency, fee-through signals can be efficiently eliminate by simple low-pass filtering of the pick-up signal. This removes the nee for intricate etection methos an simplifies the overall Q 8 Copyright c 27 by ASME

9 system esign. In its linear approximation, the phenomenon is represente by a non-homogenous LTV Mathieu equation. This approach can be use to stuy small amplitue motion an instability conitions. However, it is not aequate for stuy of large amplitue responses; higher orer terms of the parallel plate capacitive graient nee to be taken into account. In this paper we presente a etaile analysis of two types of frequency responses of parametrically riven resonators using the complete nonlinear moel. Maximum isplacement χt) max frequency response etermines if the solution is boune an physically feasible. A ifferent type of frequency response, efine as the amplitue of the main sinusoial motion at half of the AC voltage frequency χt) ω, prescribes the scale factor of such resonant evices as vibratory gyroscopes. Finite orer approximations of the complete nonlinear electro-mechanical ynamics result in significant overestimation of the actual steay-state response amplitue. Also, lower orer approximations escribe the resonant region of the frequency response imprecisely. The full nonlinear moel shoul be use in orer to precisely moel the ynamics. The escribe actuation approach is insensitive to amping an can be use to improve performance of vibratory gyroscopes an other resonant MEMS by proviing high amplitues of motion using small voltages an simplifying the etection by isolating parasitics at higher frequencies. ACKNOWLEDGMENT This research was partially supporte by the National Science Founation Grant CMS We woul like to acknowlege Cenk Acar for esign of the teste evice, Microfabrica for the fabrication of the prototypes, an Aam R. Schofiel an Professor Faryar Jabbari for iscussions. REFERENCES [] Stemme, G., 99. Resonant silicon sensors. IOP Journal of Micromechanics an Microengineering, vol., pp [2] Yazi, N., Ayazi, F., an Najafi, K., 998. Micromachine inertial sensors. Proceeings of the IEEE, vol. 86 issue 8. [3] Turner, K. K., an Zhang, W. Design an analysis of a ynamic mem chemical sensor. Proceeings of the American Control Conference 2, Arlington, VA, USA, June 25-27, 2. [4] Baxter, L. K., 996. Capacitive Sensors: Design an Applications. Wiley-IEEE Press. [5] Kovacs, G. T., 998. Micromachine Transucers Sourcebook. McGraw-Hill. [6] Cao, J., an Nguyen, C.-C. Drive amplitue epenence of micromechanical resonator series motional resistance. Digest of Technical Papers, -th International Conference on Soli-State Sensors an Actuators, Senai, Japan, June 7 -, 999. [7] Lu, M. S.-C., an Feer, G. K. Close-loop control of a parallel-plate microactuator beyon the pull-in limit. Soli-State Sensor, Actuator an Microsystems Workshop, Hilton Hea Islan, South Carolina, USA, June 2-6, 22. [8] Seeger, J., an Boser, B. Parallel-plate riven oscillations an resonant pull-in. Soli-State Sensor, Actuator an Microsystems Workshop, Hilton Hea Islan, South Carolina, USA, June 2-6, 22. [9] Fargas-Marques, A., an Shkel, A. On electrostatic actuation beyon snapping conition. IEEE Sensors Conference, Irvine, California, USA, October 3 - November [] Turner, K. L., Miller, S. A., Hartwell, P. G., MacDonal, N. C., Strogatz, S. H., an Aams, S. G., 998. Five parametric resonances in a micromechanical system. Nature, vol. 396, November. [] Nayfeh, A. H., an Younis, M. I., 25. Dynamics of mems reonators uner superharmonic an subharmonic excitations. IOP Journal of Micromechanics an Microengineering, vol. 5, pp [2] Zhang, W., Baskaran, R., an Turner, K. L., 22. Effect of cubic nonlinearity on auto-parametrically amplifie resonant mems mass sensor. Sensors an Actuators A, vol. 3578, pp.. [3] Gallacher, B. J., Buress, J. S., an Harish, K. M., 26. A control scheme for a mems electrostatic resonant gyroscope excite using combine parametric excitation an harmonic forcing. IOP Journal of Micromechanics an Microengineering, vol. 6, pp [4] Cagaser, B., Jog, A., Last, M., Leibowitz, B. S., Zhou, L., Shelton, E., Pister, K. S., an Boser, B. E. Capacitive sense feeback control for mems beam steering mirrors. Soli- State Sensor, Actuator an Microsystems Workshop, Hilton Hea Islan, South Carolina, USA, June 6 -, 24. [5] Trusov, A., Acar, C., an Shkel, A. M., 26. Comparative analysis of istribute mass micromachine gyroscopes fabricate in scs-soi an efab. Proceeings of SPIE Volume: 674, Smart Structures an Materials 26: Sensors an Smart Structures Technologies for Civil, Mechanical, an Aerospace Systems. [6] Horowitz, P., an Hill, W., 989. The Art of Electronics. Cambrige University Press. [7] Cohen, A., Frois, U., Tseng, F.-G., Zhang, G., Mansfel, F., an Will, P., 999. Efab: low-cost automate electrochemical batch fabrication of arbitrary 3- microstructures. SPIE Conference on Micromachining an Microfabrication Process Technology V. 9 Copyright c 27 by ASME

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21 Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting

More information

Structure Design and Fabrication of Silicon Resonant Micro-accelerometer Based on Electrostatic Rigidity

Structure Design and Fabrication of Silicon Resonant Micro-accelerometer Based on Electrostatic Rigidity Proceeings of the Worl Congress on Engineering 9 Vol I WCE 9, July 1 -, 9, Lonon, U.K. Structure Design an Fabrication of Silicon Resonant Micro-accelerometer Base on Electrostatic Rigiity ZHANG Feng-tian,

More information

A Simple Model for the Calculation of Plasma Impedance in Atmospheric Radio Frequency Discharges

A Simple Model for the Calculation of Plasma Impedance in Atmospheric Radio Frequency Discharges Plasma Science an Technology, Vol.16, No.1, Oct. 214 A Simple Moel for the Calculation of Plasma Impeance in Atmospheric Raio Frequency Discharges GE Lei ( ) an ZHANG Yuantao ( ) Shanong Provincial Key

More information

Optimum design of tuned mass damper systems for seismic structures

Optimum design of tuned mass damper systems for seismic structures Earthquake Resistant Engineering Structures VII 175 Optimum esign of tune mass amper systems for seismic structures I. Abulsalam, M. Al-Janabi & M. G. Al-Taweel Department of Civil Engineering, Faculty

More information

ELEC3114 Control Systems 1

ELEC3114 Control Systems 1 ELEC34 Control Systems Linear Systems - Moelling - Some Issues Session 2, 2007 Introuction Linear systems may be represente in a number of ifferent ways. Figure shows the relationship between various representations.

More information

ECE 422 Power System Operations & Planning 7 Transient Stability

ECE 422 Power System Operations & Planning 7 Transient Stability ECE 4 Power System Operations & Planning 7 Transient Stability Spring 5 Instructor: Kai Sun References Saaat s Chapter.5 ~. EPRI Tutorial s Chapter 7 Kunur s Chapter 3 Transient Stability The ability of

More information

4. CONTROL OF ZERO-SEQUENCE CURRENT IN PARALLEL THREE-PHASE CURRENT-BIDIRECTIONAL CONVERTERS

4. CONTROL OF ZERO-SEQUENCE CURRENT IN PARALLEL THREE-PHASE CURRENT-BIDIRECTIONAL CONVERTERS 4. CONRO OF ZERO-SEQUENCE CURREN IN PARAE HREE-PHASE CURREN-BIDIRECIONA CONVERERS 4. A NOVE ZERO-SEQUENCE CURREN CONRO 4.. Zero-Sequence Dynamics he parallel boost rectifier moel in Figure.4 an the parallel

More information

Introduction Basic principles Finite element formulation Nonlinear algorithms Validation & examples Oofelie::MEMS, driven by SAMCEF Field Perspectives

Introduction Basic principles Finite element formulation Nonlinear algorithms Validation & examples Oofelie::MEMS, driven by SAMCEF Field Perspectives Non linear behavior of electrostatically actuate micro-structures Dr. Ir. Stéphane Paquay, Open Engineering SA Dr. Ir. Véronique Rochus, ULg (LTAS-VIS) Dr. Ir. Stefanie Gutschmit, ULg (LTAS-VIS) Outline

More information

Design A Robust Power System Stabilizer on SMIB Using Lyapunov Theory

Design A Robust Power System Stabilizer on SMIB Using Lyapunov Theory Design A Robust Power System Stabilizer on SMIB Using Lyapunov Theory Yin Li, Stuent Member, IEEE, Lingling Fan, Senior Member, IEEE Abstract This paper proposes a robust power system stabilizer (PSS)

More information

A Novel Decoupled Iterative Method for Deep-Submicron MOSFET RF Circuit Simulation

A Novel Decoupled Iterative Method for Deep-Submicron MOSFET RF Circuit Simulation A Novel ecouple Iterative Metho for eep-submicron MOSFET RF Circuit Simulation CHUAN-SHENG WANG an YIMING LI epartment of Mathematics, National Tsing Hua University, National Nano evice Laboratories, an

More information

Harmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method

Harmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method 1 Harmonic Moelling of Thyristor Briges using a Simplifie Time Domain Metho P. W. Lehn, Senior Member IEEE, an G. Ebner Abstract The paper presents time omain methos for harmonic analysis of a 6-pulse

More information

AN INTRODUCTION TO AIRCRAFT WING FLUTTER Revision A

AN INTRODUCTION TO AIRCRAFT WING FLUTTER Revision A AN INTRODUCTION TO AIRCRAFT WIN FLUTTER Revision A By Tom Irvine Email: tomirvine@aol.com January 8, 000 Introuction Certain aircraft wings have experience violent oscillations uring high spee flight.

More information

Experimental Robustness Study of a Second-Order Sliding Mode Controller

Experimental Robustness Study of a Second-Order Sliding Mode Controller Experimental Robustness Stuy of a Secon-Orer Sliing Moe Controller Anré Blom, Bram e Jager Einhoven University of Technology Department of Mechanical Engineering P.O. Box 513, 5600 MB Einhoven, The Netherlans

More information

Chapter 6. Electromagnetic Oscillations and Alternating Current

Chapter 6. Electromagnetic Oscillations and Alternating Current hapter 6 Electromagnetic Oscillations an Alternating urrent hapter 6: Electromagnetic Oscillations an Alternating urrent (hapter 31, 3 in textbook) 6.1. Oscillations 6.. The Electrical Mechanical Analogy

More information

On Using Unstable Electrohydraulic Valves for Control

On Using Unstable Electrohydraulic Valves for Control Kailash Krishnaswamy Perry Y. Li Department of Mechanical Engineering, University of Minnesota, 111 Church St. SE, Minneapolis, MN 55455 e-mail: kk,pli @me.umn.eu On Using Unstable Electrohyraulic Valves

More information

Situation awareness of power system based on static voltage security region

Situation awareness of power system based on static voltage security region The 6th International Conference on Renewable Power Generation (RPG) 19 20 October 2017 Situation awareness of power system base on static voltage security region Fei Xiao, Zi-Qing Jiang, Qian Ai, Ran

More information

An inductance lookup table application for analysis of reluctance stepper motor model

An inductance lookup table application for analysis of reluctance stepper motor model ARCHIVES OF ELECTRICAL ENGINEERING VOL. 60(), pp. 5- (0) DOI 0.478/ v07-0-000-y An inuctance lookup table application for analysis of reluctance stepper motor moel JAKUB BERNAT, JAKUB KOŁOTA, SŁAWOMIR

More information

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013 Survey Sampling Kosuke Imai Department of Politics, Princeton University February 19, 2013 Survey sampling is one of the most commonly use ata collection methos for social scientists. We begin by escribing

More information

PRACTICE 4. CHARGING AND DISCHARGING A CAPACITOR

PRACTICE 4. CHARGING AND DISCHARGING A CAPACITOR PRACTICE 4. CHARGING AND DISCHARGING A CAPACITOR. THE PARALLEL-PLATE CAPACITOR. The Parallel plate capacitor is a evice mae up by two conuctor parallel plates with total influence between them (the surface

More information

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control 19 Eigenvalues, Eigenvectors, Orinary Differential Equations, an Control This section introuces eigenvalues an eigenvectors of a matrix, an iscusses the role of the eigenvalues in etermining the behavior

More information

EE 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.

EE 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 information

CHARACTERISTICS OF A DYNAMIC PRESSURE GENERATOR BASED ON LOUDSPEAKERS. Jože Kutin *, Ivan Bajsić

CHARACTERISTICS OF A DYNAMIC PRESSURE GENERATOR BASED ON LOUDSPEAKERS. Jože Kutin *, Ivan Bajsić Sensors an Actuators A: Physical 168 (211) 149-154 oi: 1.116/j.sna.211..7 211 Elsevier B.V. CHARACTERISTICS OF A DYNAMIC PRESSURE GENERATOR BASED ON LOUDSPEAKERS Jože Kutin *, Ivan Bajsić Laboratory of

More information

Chapter 4. Electrostatics of Macroscopic Media

Chapter 4. Electrostatics of Macroscopic Media Chapter 4. Electrostatics of Macroscopic Meia 4.1 Multipole Expansion Approximate potentials at large istances 3 x' x' (x') x x' x x Fig 4.1 We consier the potential in the far-fiel region (see Fig. 4.1

More information

Designing Information Devices and Systems I Spring 2018 Lecture Notes Note 16

Designing Information Devices and Systems I Spring 2018 Lecture Notes Note 16 EECS 16A Designing Information Devices an Systems I Spring 218 Lecture Notes Note 16 16.1 Touchscreen Revisite We ve seen how a resistive touchscreen works by using the concept of voltage iviers. Essentially,

More information

Simple Electromagnetic Motor Model for Torsional Analysis of Variable Speed Drives with an Induction Motor

Simple Electromagnetic Motor Model for Torsional Analysis of Variable Speed Drives with an Induction Motor DOI: 10.24352/UB.OVGU-2017-110 TECHNISCHE MECHANIK, 37, 2-5, (2017), 347-357 submitte: June 15, 2017 Simple Electromagnetic Motor Moel for Torsional Analysis of Variable Spee Drives with an Inuction Motor

More information

LQG FLUTTER CONTROL OF WIND TUNNEL MODEL USING PIEZO-CERAMIC ACTUATOR

LQG FLUTTER CONTROL OF WIND TUNNEL MODEL USING PIEZO-CERAMIC ACTUATOR 5 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES LQG FLUTTER CONTROL OF WIND TUNNEL MODEL USING PIEZO-CERAMIC ACTUATOR Tatsunori Kaneko* an Yasuto Asano* * Department of Mechanical Engineering,

More information

Chapter 11: Feedback and PID Control Theory

Chapter 11: Feedback and PID Control Theory Chapter 11: Feeback an D Control Theory Chapter 11: Feeback an D Control Theory. ntrouction Feeback is a mechanism for regulating a physical system so that it maintains a certain state. Feeback works by

More information

inflow outflow Part I. Regular tasks for MAE598/494 Task 1

inflow outflow Part I. Regular tasks for MAE598/494 Task 1 MAE 494/598, Fall 2016 Project #1 (Regular tasks = 20 points) Har copy of report is ue at the start of class on the ue ate. The rules on collaboration will be release separately. Please always follow the

More information

NONLINEAR ANALYSIS OF PULL IN VOLTAGE IN MICRO- CANTILEVER BEAM

NONLINEAR 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 information

Optimization of Geometries by Energy Minimization

Optimization of Geometries by Energy Minimization Optimization of Geometries by Energy Minimization by Tracy P. Hamilton Department of Chemistry University of Alabama at Birmingham Birmingham, AL 3594-140 hamilton@uab.eu Copyright Tracy P. Hamilton, 1997.

More information

Chapter 31: RLC Circuits. PHY2049: Chapter 31 1

Chapter 31: RLC Circuits. PHY2049: Chapter 31 1 Chapter 31: RLC Circuits PHY049: Chapter 31 1 LC Oscillations Conservation of energy Topics Dampe oscillations in RLC circuits Energy loss AC current RMS quantities Force oscillations Resistance, reactance,

More information

Switching Time Optimization in Discretized Hybrid Dynamical Systems

Switching Time Optimization in Discretized Hybrid Dynamical Systems Switching Time Optimization in Discretize Hybri Dynamical Systems Kathrin Flaßkamp, To Murphey, an Sina Ober-Blöbaum Abstract Switching time optimization (STO) arises in systems that have a finite set

More information

Static Analysis of Electrostatically Actuated Micro Cantilever Beam

Static Analysis of Electrostatically Actuated Micro Cantilever Beam Available online at www.scienceirect.com Proceia Engineering 51 ( 13 ) 776 78 Chemical, Civil an Mechanical Engineering Tracks of 3r Nirma University International Conference on Engineering (NUiCONE 1)

More information

Nested Saturation with Guaranteed Real Poles 1

Nested Saturation with Guaranteed Real Poles 1 Neste Saturation with Guarantee Real Poles Eric N Johnson 2 an Suresh K Kannan 3 School of Aerospace Engineering Georgia Institute of Technology, Atlanta, GA 3332 Abstract The global stabilization of asymptotically

More information

THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE

THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE Journal of Soun an Vibration (1996) 191(3), 397 414 THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE E. M. WEINSTEIN Galaxy Scientific Corporation, 2500 English Creek

More information

Design of an Industrial Distributed Controller Near Spatial Domain Boundaries

Design of an Industrial Distributed Controller Near Spatial Domain Boundaries Design of an Inustrial Distribute Controller Near Spatial Domain Bounaries Stevo Mijanovic, Greg E. Stewart, Guy A. Dumont, an Michael S. Davies ABSTRACT This paper proposes moifications to an inustrial

More information

induced _ electric _ field = = = = σ V

induced _ electric _ field = = = = σ V The Figure shows that the open-circuit voltage V (an hence the fiel strength ) is proportional to compressive stress T up to a maximum of 5 kv, which occurs at a stress of 50 MPa. Inuce electrical fiel

More information

Time-Optimal Motion Control of Piezoelectric Actuator: STM Application

Time-Optimal Motion Control of Piezoelectric Actuator: STM Application Time-Optimal Motion Control of Piezoelectric Actuator: STM Application Yongai Xu, Peter H. Mecl Abstract This paper exaes the problem of time-optimal motion control in the context of Scanning Tunneling

More information

05 The Continuum Limit and the Wave Equation

05 The Continuum Limit and the Wave Equation Utah State University DigitalCommons@USU Founations of Wave Phenomena Physics, Department of 1-1-2004 05 The Continuum Limit an the Wave Equation Charles G. Torre Department of Physics, Utah State University,

More information

EE 370L Controls Laboratory. Laboratory Exercise #7 Root Locus. Department of Electrical and Computer Engineering University of Nevada, at Las Vegas

EE 370L Controls Laboratory. Laboratory Exercise #7 Root Locus. Department of Electrical and Computer Engineering University of Nevada, at Las Vegas EE 370L Controls Laboratory Laboratory Exercise #7 Root Locus Department of Electrical an Computer Engineering University of Nevaa, at Las Vegas 1. Learning Objectives To emonstrate the concept of error

More information

Placement and tuning of resonance dampers on footbridges

Placement and tuning of resonance dampers on footbridges Downloae from orbit.tu.k on: Jan 17, 19 Placement an tuning of resonance ampers on footbriges Krenk, Steen; Brønen, Aners; Kristensen, Aners Publishe in: Footbrige 5 Publication ate: 5 Document Version

More information

This module is part of the. Memobust Handbook. on Methodology of Modern Business Statistics

This module is part of the. Memobust Handbook. on Methodology of Modern Business Statistics This moule is part of the Memobust Hanbook on Methoology of Moern Business Statistics 26 March 2014 Metho: Balance Sampling for Multi-Way Stratification Contents General section... 3 1. Summary... 3 2.

More information

PARALLEL-PLATE CAPACITATOR

PARALLEL-PLATE CAPACITATOR Physics Department Electric an Magnetism Laboratory PARALLEL-PLATE CAPACITATOR 1. Goal. The goal of this practice is the stuy of the electric fiel an electric potential insie a parallelplate capacitor.

More information

Math 342 Partial Differential Equations «Viktor Grigoryan

Math 342 Partial Differential Equations «Viktor Grigoryan Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite

More information

Chapter 11: Feedback and PID Control Theory

Chapter 11: Feedback and PID Control Theory Chapter 11: Feeback an D Control Theory Chapter 11: Feeback an D Control Theory. ntrouction Feeback is a mechanism for regulating a physical system so that it maintains a certain state. Feeback works by

More information

Two Dimensional Numerical Simulator for Modeling NDC Region in SNDC Devices

Two Dimensional Numerical Simulator for Modeling NDC Region in SNDC Devices Journal of Physics: Conference Series PAPER OPEN ACCESS Two Dimensional Numerical Simulator for Moeling NDC Region in SNDC Devices To cite this article: Dheeraj Kumar Sinha et al 2016 J. Phys.: Conf. Ser.

More information

Chapter 9 Method of Weighted Residuals

Chapter 9 Method of Weighted Residuals Chapter 9 Metho of Weighte Resiuals 9- Introuction Metho of Weighte Resiuals (MWR) is an approimate technique for solving bounary value problems. It utilizes a trial functions satisfying the prescribe

More information

Semiclassical analysis of long-wavelength multiphoton processes: The Rydberg atom

Semiclassical analysis of long-wavelength multiphoton processes: The Rydberg atom PHYSICAL REVIEW A 69, 063409 (2004) Semiclassical analysis of long-wavelength multiphoton processes: The Ryberg atom Luz V. Vela-Arevalo* an Ronal F. Fox Center for Nonlinear Sciences an School of Physics,

More information

PCCP PAPER. 1 Introduction. A. Nenning,* A. K. Opitz, T. M. Huber and J. Fleig. View Article Online View Journal View Issue

PCCP PAPER. 1 Introduction. A. Nenning,* A. K. Opitz, T. M. Huber and J. Fleig. View Article Online View Journal View Issue PAPER View Article Online View Journal View Issue Cite this: Phys. Chem. Chem. Phys., 2014, 16, 22321 Receive 4th June 2014, Accepte 3r September 2014 DOI: 10.1039/c4cp02467b www.rsc.org/pccp 1 Introuction

More information

Role of parameters in the stochastic dynamics of a stick-slip oscillator

Role of parameters in the stochastic dynamics of a stick-slip oscillator Proceeing Series of the Brazilian Society of Applie an Computational Mathematics, v. 6, n. 1, 218. Trabalho apresentao no XXXVII CNMAC, S.J. os Campos - SP, 217. Proceeing Series of the Brazilian Society

More information

Nonlinear Dielectric Response of Periodic Composite Materials

Nonlinear Dielectric Response of Periodic Composite Materials onlinear Dielectric Response of Perioic Composite aterials A.G. KOLPAKOV 3, Bl.95, 9 th ovember str., ovosibirsk, 639 Russia the corresponing author e-mail: agk@neic.nsk.su, algk@ngs.ru A. K.TAGATSEV Ceramics

More information

Homework 7 Due 18 November at 6:00 pm

Homework 7 Due 18 November at 6:00 pm Homework 7 Due 18 November at 6:00 pm 1. Maxwell s Equations Quasi-statics o a An air core, N turn, cylinrical solenoi of length an raius a, carries a current I Io cos t. a. Using Ampere s Law, etermine

More information

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation

More information

RETROGRADE WAVES IN THE COCHLEA

RETROGRADE WAVES IN THE COCHLEA August 7, 28 18:2 WSPC - Proceeings Trim Size: 9.75in x 6.5in retro wave 1 RETROGRADE WAVES IN THE COCHLEA S. T. NEELY Boys Town National Research Hospital, Omaha, Nebraska 68131, USA E-mail: neely@boystown.org

More information

3-D FEM Modeling of fiber/matrix interface debonding in UD composites including surface effects

3-D FEM Modeling of fiber/matrix interface debonding in UD composites including surface effects IOP Conference Series: Materials Science an Engineering 3-D FEM Moeling of fiber/matrix interface eboning in UD composites incluing surface effects To cite this article: A Pupurs an J Varna 2012 IOP Conf.

More information

Lecture 6: Control of Three-Phase Inverters

Lecture 6: Control of Three-Phase Inverters Yoash Levron The Anrew an Erna Viterbi Faculty of Electrical Engineering, Technion Israel Institute of Technology, Haifa 323, Israel yoashl@ee.technion.ac.il Juri Belikov Department of Computer Systems,

More information

Optimal LQR Control of Structures using Linear Modal Model

Optimal LQR Control of Structures using Linear Modal Model Optimal LQR Control of Structures using Linear Moal Moel I. Halperin,2, G. Agranovich an Y. Ribakov 2 Department of Electrical an Electronics Engineering 2 Department of Civil Engineering Faculty of Engineering,

More information

6 General properties of an autonomous system of two first order ODE

6 General properties of an autonomous system of two first order ODE 6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x

More information

Modeling time-varying storage components in PSpice

Modeling time-varying storage components in PSpice Moeling time-varying storage components in PSpice Dalibor Biolek, Zenek Kolka, Viera Biolkova Dept. of EE, FMT, University of Defence Brno, Czech Republic Dept. of Microelectronics/Raioelectronics, FEEC,

More information

Capacitive Sensor Interfaces

Capacitive Sensor Interfaces Capacitive Sensor Interfaces Bernhard E. Boser Berkeley Sensor & Actuator Center Dept. of Electrical Engineering and Computer Sciences University of California, Berkeley Capacitive Sensor Interfaces 1996

More information

The Principle of Least Action

The Principle of Least Action Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of

More information

A Comparison between a Conventional Power System Stabilizer (PSS) and Novel PSS Based on Feedback Linearization Technique

A Comparison between a Conventional Power System Stabilizer (PSS) and Novel PSS Based on Feedback Linearization Technique J. Basic. Appl. Sci. Res., ()9-99,, TextRoa Publication ISSN 9-434 Journal of Basic an Applie Scientific Research www.textroa.com A Comparison between a Conventional Power System Stabilizer (PSS) an Novel

More information

Laplacian Cooperative Attitude Control of Multiple Rigid Bodies

Laplacian Cooperative Attitude Control of Multiple Rigid Bodies Laplacian Cooperative Attitue Control of Multiple Rigi Boies Dimos V. Dimarogonas, Panagiotis Tsiotras an Kostas J. Kyriakopoulos Abstract Motivate by the fact that linear controllers can stabilize the

More information

A Review of Feedforward Control Approaches in Nanopositioning for High-Speed SPM

A Review of Feedforward Control Approaches in Nanopositioning for High-Speed SPM Garrett M. Clayton Department of Mechanical Engineering, Villanova University, Villanova, PA 1985 e-mail: garrett.clayton@villanova.eu Szuchi Tien Department of Mechanical Engineering, National Cheng Kung

More information

The canonical controllers and regular interconnection

The canonical controllers and regular interconnection Systems & Control Letters ( www.elsevier.com/locate/sysconle The canonical controllers an regular interconnection A.A. Julius a,, J.C. Willems b, M.N. Belur c, H.L. Trentelman a Department of Applie Mathematics,

More information

Essential Considerations for Buckling Analysis

Essential Considerations for Buckling Analysis Worlwie Aerospace Conference an Technology Showcase, Toulouse, France, Sept. 24-26, 2001 Essential Consierations for Buckling Analysis 2001-120 Sang H. Lee MSC.Software Corporation, 2 MacArthur Place,

More information

Advanced Partial Differential Equations with Applications

Advanced Partial Differential Equations with Applications MIT OpenCourseWare http://ocw.mit.eu 18.306 Avance Partial Differential Equations with Applications Fall 2009 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.eu/terms.

More information

Damage detection of shear building structure based on FRF response variation

Damage detection of shear building structure based on FRF response variation , pp.18-5 http://x.oi.org/10.1457/astl.013.3.05 Damage etection of shear builing structure base on FRF response variation Hee-Chang Eun 1,*, Su-Yong Park 1, Rae-Jung im 1 1 Dept. of Architectural Engineering,

More information

VIRTUAL STRUCTURE BASED SPACECRAFT FORMATION CONTROL WITH FORMATION FEEDBACK

VIRTUAL STRUCTURE BASED SPACECRAFT FORMATION CONTROL WITH FORMATION FEEDBACK AIAA Guiance, Navigation, an Control Conference an Exhibit 5-8 August, Monterey, California AIAA -9 VIRTUAL STRUCTURE BASED SPACECRAT ORMATION CONTROL WITH ORMATION EEDBACK Wei Ren Ranal W. Bear Department

More information

Chapter 11: Feedback and PID Control Theory

Chapter 11: Feedback and PID Control Theory Chapter 11: Feeback an D Control Theory Chapter 11: Feeback an D Control Theory. ntrouction Feeback is a mechanism for regulating a physical system so that it maintains a certain state. Feeback works by

More information

Predictive control of synchronous generator: a multiciterial optimization approach

Predictive control of synchronous generator: a multiciterial optimization approach Preictive control of synchronous generator: a multiciterial optimization approach Marián Mrosko, Eva Miklovičová, Ján Murgaš Abstract The paper eals with the preictive control esign for nonlinear systems.

More information

H 2 GUARANTEED COST CONTROL DESIGN AND IMPLEMENTATION FOR DUAL-STAGE HARD DISK DRIVE TRACK-FOLLOWING SERVOS. Richard Conway

H 2 GUARANTEED COST CONTROL DESIGN AND IMPLEMENTATION FOR DUAL-STAGE HARD DISK DRIVE TRACK-FOLLOWING SERVOS. Richard Conway Proceeings of the ASME 211 Dynamic Systems an Control Conference DSCC211 October 31 - November 2, 211, Arlington, VA, USA DSCC211-695 H 2 GUARANTEED COST CONTROL DESIGN AND IMPLEMENTATION FOR DUAL-STAGE

More information

Four Degrees-of-Freedom Micromachined Gyroscope

Four Degrees-of-Freedom Micromachined Gyroscope Microsystems Laboratory Technical Report Four Degrees-of-Freedom Micromachined Gyroscope Cenk Acar 23 October 2001 Technical Report No: MSL-01003 cfl2001 Cenk Acar Contents Contents List of Figures Abstract

More information

Chaos, Solitons and Fractals Nonlinear Science, and Nonequilibrium and Complex Phenomena

Chaos, Solitons and Fractals Nonlinear Science, and Nonequilibrium and Complex Phenomena Chaos, Solitons an Fractals (7 64 73 Contents lists available at ScienceDirect Chaos, Solitons an Fractals onlinear Science, an onequilibrium an Complex Phenomena journal homepage: www.elsevier.com/locate/chaos

More information

State-Space Model for a Multi-Machine System

State-Space Model for a Multi-Machine System State-Space Moel for a Multi-Machine System These notes parallel section.4 in the text. We are ealing with classically moele machines (IEEE Type.), constant impeance loas, an a network reuce to its internal

More information

Dynamics of a Microgrid Supplied by Solid Oxide Fuel Cells 1

Dynamics of a Microgrid Supplied by Solid Oxide Fuel Cells 1 Bulk Power System Dynamics an Control - VII, August 19-4, 007, Charleston, South Carolina, USA 1 Dynamics of a Microgri Supplie by Soli Oxie Fuel Cells 1 Eric M. Fleming Ian A. Hiskens Department of Electrical

More information

APPPHYS 217 Thursday 8 April 2010

APPPHYS 217 Thursday 8 April 2010 APPPHYS 7 Thursay 8 April A&M example 6: The ouble integrator Consier the motion of a point particle in D with the applie force as a control input This is simply Newton s equation F ma with F u : t q q

More information

Quantum optics of a Bose-Einstein condensate coupled to a quantized light field

Quantum optics of a Bose-Einstein condensate coupled to a quantized light field PHYSICAL REVIEW A VOLUME 60, NUMBER 2 AUGUST 1999 Quantum optics of a Bose-Einstein conensate couple to a quantize light fiel M. G. Moore, O. Zobay, an P. Meystre Optical Sciences Center an Department

More information

water adding dye partial mixing homogenization time

water adding dye partial mixing homogenization time iffusion iffusion is a process of mass transport that involves the movement of one atomic species into another. It occurs by ranom atomic jumps from one position to another an takes place in the gaseous,

More information

ensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y

ensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y Ph195a lecture notes, 1/3/01 Density operators for spin- 1 ensembles So far in our iscussion of spin- 1 systems, we have restricte our attention to the case of pure states an Hamiltonian evolution. Toay

More information

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments 2 Conference on Information Sciences an Systems, The Johns Hopkins University, March 2, 2 Time-of-Arrival Estimation in Non-Line-Of-Sight Environments Sinan Gezici, Hisashi Kobayashi an H. Vincent Poor

More information

NOISE ANALYSIS IN PARAMETRIC RESONANCE BASED MASS SENSING. Wenhua Zhang and Kimberly L. Turner

NOISE ANALYSIS IN PARAMETRIC RESONANCE BASED MASS SENSING. Wenhua Zhang and Kimberly L. Turner Proceedings of IMECE04 004 ASME International Mechanical Engineering Congress and Exposition November 1-0, 004, Anaheim, California USA IMECE004-6141 NOISE ANALYSIS IN PARAMETRIC RESONANCE BASED MASS SENSING

More information

Experimental analysis of spring hardening and softening nonlinearities in. microelectromechanical oscillators. Sarah Johnson

Experimental analysis of spring hardening and softening nonlinearities in. microelectromechanical oscillators. Sarah Johnson Experimental analysis of spring hardening and softening nonlinearities in microelectromechanical oscillators. Sarah Johnson Department of Physics, University of Florida Mentored by Dr. Yoonseok Lee Abstract

More information

SYNCHRONOUS SEQUENTIAL CIRCUITS

SYNCHRONOUS SEQUENTIAL CIRCUITS CHAPTER SYNCHRONOUS SEUENTIAL CIRCUITS Registers an counters, two very common synchronous sequential circuits, are introuce in this chapter. Register is a igital circuit for storing information. Contents

More information

IERCU. Institute of Economic Research, Chuo University 50th Anniversary Special Issues. Discussion Paper No.210

IERCU. Institute of Economic Research, Chuo University 50th Anniversary Special Issues. Discussion Paper No.210 IERCU Institute of Economic Research, Chuo University 50th Anniversary Special Issues Discussion Paper No.210 Discrete an Continuous Dynamics in Nonlinear Monopolies Akio Matsumoto Chuo University Ferenc

More information

Effects of anharmonicity of current-phase relation in Josephson junctions

Effects of anharmonicity of current-phase relation in Josephson junctions Low Temperature Physics/Fizika Nizkikh Temperatur, 05, v. 4, No. 4, pp. 35 337 Effects of anharmonicity of current-phase relation in Josephson junctions (Review Article) I.N. Askerzae Computer Engineering

More information

Adaptive Gain-Scheduled H Control of Linear Parameter-Varying Systems with Time-Delayed Elements

Adaptive Gain-Scheduled H Control of Linear Parameter-Varying Systems with Time-Delayed Elements Aaptive Gain-Scheule H Control of Linear Parameter-Varying Systems with ime-delaye Elements Yoshihiko Miyasato he Institute of Statistical Mathematics 4-6-7 Minami-Azabu, Minato-ku, okyo 6-8569, Japan

More information

Sensors & Transducers 2015 by IFSA Publishing, S. L.

Sensors & Transducers 2015 by IFSA Publishing, S. L. Sensors & Transucers, Vol. 184, Issue 1, January 15, pp. 53-59 Sensors & Transucers 15 by IFSA Publishing, S. L. http://www.sensorsportal.com Non-invasive an Locally Resolve Measurement of Soun Velocity

More information

DAMAGE DETECTIONS IN NONLINEAR VIBRATING THERMALLY LOADED STRUCTURES 1

DAMAGE DETECTIONS IN NONLINEAR VIBRATING THERMALLY LOADED STRUCTURES 1 11 th National Congress on Theoretical an Applie Mechanics, 2-5 Sept. 2009, Borovets, Bulgaria DAMAGE DETECTIONS IN NONLINEAR VIBRATING THERMALLY LOADED STRUCTURES 1 E. MANOACH Institute of Mechanics,

More information

Simulation of Angle Beam Ultrasonic Testing with a Personal Computer

Simulation of Angle Beam Ultrasonic Testing with a Personal Computer Key Engineering Materials Online: 4-8-5 I: 66-9795, Vols. 7-73, pp 38-33 oi:.48/www.scientific.net/kem.7-73.38 4 rans ech ublications, witzerlan Citation & Copyright (to be inserte by the publisher imulation

More information

NONLINEAR BEHAVIOR OF A PARAMETRIC RESONANCE-BASED MASS SENSOR

NONLINEAR BEHAVIOR OF A PARAMETRIC RESONANCE-BASED MASS SENSOR Proceedings of IMECE00 ASME International Mechanical Engineering Congress & Exposition November 17-, 00, New Orleans, Louisiana IMECE00-3361 NONLINEAR BEHAVIOR OF A PARAMETRIC RESONANCE-BASED MASS SENSOR

More information

Last lecture. Today s menu. Capacitive sensing elements. Capacitive sensing elements (cont d...) Examples. General principle

Last lecture. Today s menu. Capacitive sensing elements. Capacitive sensing elements (cont d...) Examples. General principle Last lecture esistive sensing elements: Displacement sensors (potentiometers). Temperature sensors. Strain gauges. Deflection briges. Toay s menu Capacitive sensing elements. Inuctive sensing elements.

More information

ELECTRONIC SYSTEMS. Real operational amplifiers. Electronic Systems - C2 28/04/ DDC Storey 1

ELECTRONIC SYSTEMS. Real operational amplifiers. Electronic Systems - C2 28/04/ DDC Storey 1 Politecnico i Torino ICT school Lesson C2 Real Op mp parameters an moel ELECTRONIC SYSTEMS C OPERTIONL MPLIFIERS C.2 Real Op mp an examples» Real Op mp parameters an moel» Static an ynamic parameters»

More information

LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION

LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION The Annals of Statistics 1997, Vol. 25, No. 6, 2313 2327 LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION By Eva Riccomagno, 1 Rainer Schwabe 2 an Henry P. Wynn 1 University of Warwick, Technische

More information

Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing

Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing Course Project for CDS 05 - Geometric Mechanics John M. Carson III California Institute of Technology June

More information

Designing Information Devices and Systems II Spring 2019 A. Sahai, J. Roychowdhury, K. Pister Midterm 1: Practice

Designing Information Devices and Systems II Spring 2019 A. Sahai, J. Roychowdhury, K. Pister Midterm 1: Practice EES 16B Designing Information Devices an Systems II Spring 019 A. Sahai, J. Roychowhury, K. Pister Miterm 1: Practice 1. Speaker System Your job is to construct a speaker system that operates in the range

More information

Adaptive Control of the Boost DC-AC Converter

Adaptive Control of the Boost DC-AC Converter Aaptive Control of the Boost DC-AC Converter Carolina Albea-Sanchez, Carlos Canuas De Wit, Francisco Gorillo Alvarez To cite this version: Carolina Albea-Sanchez, Carlos Canuas De Wit, Francisco Gorillo

More information

Damage identification based on incomplete modal data and constrained nonlinear multivariable function

Damage identification based on incomplete modal data and constrained nonlinear multivariable function Journal of Physics: Conference Series PAPER OPEN ACCESS Damage ientification base on incomplete moal ata an constraine nonlinear multivariable function To cite this article: S S Kourehli 215 J. Phys.:

More information

Quasi optimal feedforward control of a very low frequency high-voltage test system

Quasi optimal feedforward control of a very low frequency high-voltage test system Preprints of the 9th Worl Congress The International Feeration of Automatic Control Quasi optimal feeforwar control of a very low frequency high-voltage test system W. Kemmetmüller S. Eberharter A. Kugi

More information

Vibration Analysis of Railway Tracks Forced by Distributed Moving Loads

Vibration Analysis of Railway Tracks Forced by Distributed Moving Loads IJR International Journal of Railway Vol. 6, No. 4 / December 13, pp. 155-159 The Korean Society for Railway Vibration Analysis of Railway Tracks Force by Distribute Moving Loas Sinyeob Lee*, Dongkyu Kim*,

More information