Generalized parametric resonance in electrostatically actuated microelectromechanical oscillators
|
|
- Beryl Hood
- 6 years ago
- Views:
Transcription
1 Journal of Sound and Vibration 296 (26) JOURNAL OF SOUND AND VIBRATION Generalized parametric resonance in electrostaticall actuated microelectromechanical oscillators Jeffre F. Rhoads a,, Steven W. Shaw a, Kimberl L. Turner b, Jeff Moehlis b, Barr E. DeMartini b, Wenhua Zhang b a Department of Mechanical Engineering, Michigan State Universit, East Lansing, MI 48824, USA b Department of Mechanical Engineering, Universit of California - Santa Barbara, Santa Barbara, CA 9316, USA Received 8 June 25; received in revised form 28 Februar 26; accepted 8 March 26 Available online 23 Ma 26 Abstract This paper investigates the dnamic response of a class of electrostaticall driven microelectromechanical (MEM) oscillators. The particular sstems of interest are those which feature parametric ecitation that arises from forces produced b fluctuating voltages applied across comb drives. These sstems are known to ehibit a wide range of behaviors, some of which have escaped eplanation or prediction. In this paper we eamine a general governing equation of motion for these sstems and use it to provide a complete description of the dnamic response and its dependence on the sstem parameters. The defining feature of this equation is that both the linear and cubic terms feature parametric ecitation which, in comparison to the case of purel linear parametric ecitation (e.g. the Mathieu equation), significantl complicates the sstem s dnamics. One consequence is that an effective nonlinearit for the overall sstem cannot be defined. Instead, the sstem features separate effective nonlinearities for each branch of its nontrivial response. As such, it can ehibit not onl hardening and softening nonlinearities, but also mied nonlinearities, wherein the response branches in the sstem s frequenc response bend toward or awa from one another near resonance. This paper includes some brief background information on the equation of motion under consideration, an outline of the analtical techniques used to reach the aforementioned results, stabilit results for the responses in question, a numerical eample, eplored using simulation, of a MEM oscillator which features this nonlinear behavior, and preliminar eperimental results, taken from an actual MEM device, which show evidence of the analticall predicted behavior. Practical issues pertaining to the design of parametricall ecited MEM devices are also considered. r 26 Elsevier Ltd. All rights reserved. 1. Introduction The recent analsis of parametricall ecited microelectromechanical (MEM) oscillators, b the authors and others, has led to the consideration of a certain class of nonlinear equations of motion [1 3]. The defining feature of these equations is that parametric ecitation acts on both the linear and cubic terms. As such, the dnamics of these oscillators are quite comple in comparison to a tpical Mathieu sstem, which features Corresponding author. Tel.: ; fa: address: rhoadsje@msu.edu (J.F. Rhoads) X/$ - see front matter r 26 Elsevier Ltd. All rights reserved. doi:1.116/j.jsv
2 798 ARTICLE IN PRESS J.F. Rhoads et al. / Journal of Sound and Vibration 296 (26) onl linear parametric ecitation. One interesting result of this complication is that a single effective nonlinearit for the overall sstem cannot be defined. Rather, the sstem features separate effective nonlinearities for each branch of its nontrivial response. Accordingl, the sstem can ehibit not onl the tpical hardening and softening nonlinearities, but also mied nonlinearities wherein the response branches in the sstem s frequenc response bend toward or awa from one another near resonance. In fact, such behavior has been observed in eperiments, as described subsequentl. The paper begins in Section 2 with a brief description of the tpes of MEM devices of interest, and a derivation of their governing equation of motion. General analtical results for the undamped sstem, including information on the eistence, stabilit, and bifurcations of the stead-state solutions, are given in Sections 3 and 4. These results include detailed descriptions of the resulting frequenc response curves and the associated phase planes. To facilitate understanding of the aforementioned results, a representative numerical eample is eplored in Section 5. In Section 6, the effects of light damping on the sstem are considered. Section 7 eplores the given results in the contet of a specific MEM oscillator, provides the pertinent results of preliminar eperimental investigations, and discusses practical design issues associated with parametricall ecited MEM devices. The paper concludes in Section 8 with a summar of the present work and some comments on possible directions of future work in this area. It should be noted that the bulk of this paper describes a detailed analsis of the governing equation of motion; this ma be of interest to the nonlinear vibrations communit, but of less direct interest to the MEMS communit. Those readers not interested in the details of the analsis can focus on Section 2, in which the MEMS model is emploed, and Section 7, which presents the general response features of parametricall ecited MEMS, including an approach for evaluating and predicting the response features of these devices. 2. Parametricall ecited MEM oscillators The impetus of the present stud was the analsis of single-degree-of-freedom (sdof) parametricall ecited MEM oscillators, such as the one shown in Fig. 1. Such oscillators consist of a shuttle mass, namel the oscillator s backbone, B, connected to a substrate via four folded beam springs, S, and ecited b a pair of non-interdigitated electrostatic comb drives, N, which are powered b a voltage source. These devices are Fig. 1. A representative parametricall ecited MEM oscillator. The oscillator s backbone is labeled B, S designates the mechanical beam springs, and N designates the non-interdigitated comb drives, which are used to provide the parametric ecitation.
3 J.F. Rhoads et al. / Journal of Sound and Vibration 296 (26) being proposed for use in a number of sensing and filtering applications [1,3,4]. Previous works, including those b the authors, have verified that these parametricall ecited oscillators ehibit motion consistent with a rather simple-looking equation of motion, presented in Eq. (1) [1,3,5,6]. This is based on the following development. To begin, the oscillator s motion is assumed to be described b the movement of the shuttle mass in one direction. The forces acting on the sstem include those from the springs, the electrostatic forces from the comb drives, and a force that represents dissipative effects, which arise primaril from aerodnamic effects. These dissipative effects are generall small, especiall when these devices are operated in a vacuum, and are tpicall modeled b a simple viscous damping force, although more complicated models have been proposed (see, for eample, Refs. [7 1]). The resulting equation of motion is given b m þ c _ þ F r ðþþf es ð; tþ ¼, (1) where F r ðþ represents the elastic restoring force provided b the springs and F es ð; tþ represents the electrostatic restoring and driving forces provided b the non-interdigitated electrostatic comb drives. As previous works have shown, the aforementioned forces can be accuratel represented b cubic functions of displacement [1,3]. The restoring force from the mechanical springs is modeled b F r ðþ ¼k 1 þ k 3 3, (2) where k 1 and k 3 are the effective linear and nonlinear spring coefficients from the beams. The electrostatic forces arising from the comb drives can accuratel be modeled for small displacements as F es ð; tþ ¼ðr 1A þ r 3A 3 ÞV 2 ðtþ, (3) where r 1A and r 3A are electrostatic coefficients that depend on the phsical dimensions of the electrostatic comb drives, and VðtÞ is the voltage applied across the drives. In order to provide harmonic ecitation to the device, an AC voltage of the form pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi VðtÞ ¼V A 1 þ cos ot (4) is applied across the comb drives [1,3,6]. Substituting these forces back into Eq. (1) results in an equation of motion for the oscillator s backbone of the form [1,6]: m þ c _ þ k 1 þ k 3 3 þðr 1A þ r 3A 3 ÞV 2 Að1 þ cos otþ ¼. (5) Time and displacement in Eq. (5) are rescaled according to t ¼ o t, (6) where o is the purel elastic natural frequenc defined as rffiffiffiffiffi k 1 o ¼ (7) m and z ¼, (8) where is a characteristic length of the sstem (e.g. the length of the backbone). This scaling, assuming that the nondimensional damping, mechanical nonlinearities, and electrostatic forces are small (which is entirel consistent with the operation of such oscillators, especiall near resonance), ields a nondimensional equation of motion of the form z þ 2ezz þ zð1 þ el 1 þ el 1 cos OtÞþez 3 ðw þ l 3 þ l 3 cos OtÞ ¼, (9) where the new derivative operator and nondimensional parameters are defined as stated in Table 1 [1,6] and e represents a small parameter introduced for analtical purposes. The equation of motion considered in this work is a generalization of Eq. (9), wherein the effective stiffness and ecitation coefficients (the l s) are allowed to differ, in both the linear and nonlinear terms. This arises
4 8 ARTICLE IN PRESS J.F. Rhoads et al. / Journal of Sound and Vibration 296 (26) Table 1 Nondimensional parameter definitions [1,6] Definition ðþ ¼ dðþ dt ez ¼ c 2mo el 1 ¼ r 1AV 2 A k 1 O ¼ o o ew ¼ k 3 2 k 1 el 3 ¼ 2 r 3AV 2 A k 1 Nondimensional parameter Scaled time derivative Scaled damping ratio Linear electrostatic ecitation amplitude Nondimensional ecitation frequenc Nonlinear mechanical stiffness coefficient Nonlinear electrostatic ecitation coefficient naturall in the analsis of parametricall ecited MEM oscillators with two sets of non-interdigitated comb drives, one of which provides an AC voltage, as described above, while the other provides a DC bias to the sstem. This bias does not offset the sstem equilibrium in the usual sense, but allows one to independentl var the effective electrostatic stiffness and ecitation coefficients, resulting in Eq. (1). This arrangement of comb drives allows for more fleibilit in tuning the response of such devices, so that designers can achieve certain desirable response features [1,5,6]. More generall, this governing equation of motion is representative of a number of parametricall ecited dnamical sstems, as highlighted in Refs. [11 13]. While works such as [1,3,5,6] eamine the nonlinear dnamics of these parametricall ecited MEM oscillators, their nonlinear response has not been full eplored, and a number of apparent anomalies have been observed in eperimental investigations of these devices [2,3,14]. The analsis in the following sections provides a description for the entire range of possible responses for these sstems (under the small damping and small response amplitude assumptions). This is followed b a specific eample using the parameters from a device similar to that shown in Fig. 1. One of the more interesting features of the equation of motion presented in Eq. (9) is that several of the coefficients of the equation, including the effective linear and nonlinear stiffnesses, depend on the amplitude of the AC voltage, V A. This results in the fact that the qualitative nature of the sstem s response generall depends on the amplitude of the ecitation, and ma differ for different input amplitudes. The analsis presented here describes a complete picture of this behavior, and provides designers of these devices with useful predictive tools. In Section 7, an eample of this process and preliminar eperimental results that verif the utilit of the approach are given. 3. The equation of motion and perturbation analsis The equation of motion of interest here is of the form z þ 2ezz þ zð1 þ en 1 þ el 1 cos OtÞþez 3 ðg 3 þ l 3 cos OtÞ ¼, (1) where e is a scaling parameter introduced for analsis and prime designates the derivative with respect to t [1,5,6,15,16]. Note that this is a special case of generalized parametric ecitation in that there is not a phase difference in the ecitation seen b the linear and nonlinear terms. However, in man sstems the ecitation arises from a single source, resulting in zero phase difference (see, for eample, Refs. [11 13,17,18]). This is the case for the present stud, and this equation proves sufficient for capturing the dnamic response of the MEM devices under consideration.
5 J.F. Rhoads et al. / Journal of Sound and Vibration 296 (26) The method of averaging is emploed to analze Eq. (1). To aid this, a standard coordinate transformation is introduced to convert the equation to amplitude and phase coordinates, namel, zðtþ ¼aðtÞ cos Ot 2 þ cðtþ, (11) z ðtþ ¼ aðtþ O Ot sin 2 2 þ cðtþ. (12) Note that this conversion requires that a certain constraint equation involving a and c be satisfied, just as in the method of variation of parameters; see, for eample, Ref. [19]. In order to capture the dnamics near the principal parametric resonance, an ecitation frequenc detuning parameter s is introduced, defined b, O ¼ 2 þ es, (13) where s measures the closeness of the ecitation frequenc to the principal parametric resonance condition. Averaging the resulting equations, now in terms of amplitude and phase coordinates, over the period 4p=O in the t-domain produces the sstem s averaged equations [1,5,6]. These are given b a ¼ 1 8 ae½ 8z þð2l 1 þ a 2 l 3 Þ sin 2cŠþOðe 2 Þ, (14) c ¼ 1 8 e½3a2 g 3 þ 4n 1 4s þ 2ðl 1 þ a 2 l 3 Þ cos 2cŠþOðe 2 Þ. (15) Note that the presence of parametric ecitation in the cubic term, captured b the parameter l 3, introduces additional coupling in these equations that is not present in the case of purel linear parametric ecitation. Consequentl, the nontrivial stead-state solutions of the sstem take on a more complicated form and the sstem has much more interesting response characteristics. 4. Stead-state responses, their stabilit, and their effective nonlinearities Since the characteristic form of the sstem nonlinearit (i.e., hardening or softening) is a critical feature of the response for current purposes, and this is largel unaffected b damping, zero damping ðz ¼ Þ is assumed to simplif the analsis. (Note that while closed-form solutions are obtainable for nonzero damping, the are inherentl complicated, and thus not well suited for detailed analsis or eas interpretation.) We will consider the effects of small damping subsequentl. Stead-state responses are captured b setting ða ; c Þ¼ð; Þ. Eamining the resulting equations reveals that the sstem has a trivial solution and three distinct sets of nontrivial solution branches. The first two nontrivial sets appear in pairs with amplitudes and associated phases given b sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi 4s þ 2l 1 4n 1 a 1 ¼, (16) 3g 3 2l 3 c 1 ¼ p 2, (17) and sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi a 2 ¼ 4s 2l 1 4n 1, 3g 3 þ 2l 3 (18) c 2 ¼. (19) Note that each pair represents the same phsical response; the appear this wa due to the phase relation of this subharmonic response relative to the ecitation. Likewise, all responses with magnitude p phase shifts represent the same phsical response, and thus are not included. For each of these solutions, the signs of the terms under the square root determine the frequenc range over which these branches are real valued, and thus phsicall meaningful. The numerators of these terms are dictated b linear frequenc terms, and the
6 82 ARTICLE IN PRESS J.F. Rhoads et al. / Journal of Sound and Vibration 296 (26) denominators are set b nonlinear terms. Of direct interest here is the role plaed b the nonlinearities, which differ for each of these frequenc-dependent branches. To describe these effects, effective nonlinearities Z 1 and Z 2 are defined that dictate the behavior of these branches. These nonlinearities are taken from the denominator of the epressions above and are given b [1,5] Z 1 ¼ 3g 3 2l 3 (2) and Z 2 ¼ 3g 3 þ 2l 3. (21) The case Z 1 o and Z 2 o results in the usual softening nonlinearit. Likewise, the usual hardening nonlinearit occurs for Z 1 4 and Z 2 4. However, due to the nature of this equation of motion, two mied cases also eist, namel, Z 1 4 and Z 2 o, and Z 1 o and Z 2 4, which correspond to the two branches delineated in Eqs. (16) and (18) bending toward or awa from one another, as determined b l 1 and n 1. This result is summarized in Fig. 2, which delineates the various response regions within the g 3 l 3 parameter space. Note that a transition between response tpes occurs when one or both of the effective nonlinearities equals zero, in which case the bifurcation associated with crossing the boundar of the parametric instabilit (given b an Arnold tongue) becomes degenerate. Accordingl, phenomena near these degeneracies cannot be captured without the inclusion of high-order nonlinear terms, which are not considered here. The solutions for the third set of nontrivial responses have constant amplitude (for zero damping) and are given, with their phase, b sffiffiffiffiffiffiffiffiffiffiffi 2l 1 a 3 ¼, (22) l 3 c 3 ¼ 1 2 arccos 3g 3l 1 þ 2l 3 n 1 2l 3 s. (23) l 1 l 3 Note that there will be four such solutions, onl two of which will be phsicall distinct (as evident in the net section), corresponding to each branch of the respective inverse cosine function. Accordingl, these solutions will appear in either the upper or lower half-plane of the parameter space delineated in Fig. 2, as dictated b the sign of l 1. As such, the sstem undergoes a bifurcation across the l 3 ¼ ais in Fig. 2, which corresponds to the creation or annihilation of the constant amplitude solution. It should be noted however, that this creation (or annihilation) is a smooth process, as the solution appears from positive or negative infinit due to the presence of l 3 in the solution s denominator. As will be shown shortl, the presence of these additional solution branches significantl complicates the sstem s dnamics within certain parameter regimes, especiall λ η2 < η2> III IV IIb Va Vb γ 3 IIa I VI η2 < η2> η 1 < η 1 < η 1 > η 1 > Fig. 2. The g 3 l 3 parameter space [1]. Those regions designated I and VI ehibit hardening (or quasi-hardening) nonlinear characteristics, those designated III and IV ehibit softening (or quasi-softening) nonlinear characteristics, and those designated IIa, IIb, Va, and Vb ehibit mied nonlinear characteristics.
7 J.F. Rhoads et al. / Journal of Sound and Vibration 296 (26) when viewed from a frequenc response perspective. Before turning to the generation of frequenc response curves, the stabilit of the various response branches is first determined. Though the local stabilit of stead-state responses is generall addressed b considering the local linear behavior of the averaged equations near the stead-state responses in the ða; cþ polar coordinate space, singularities at the origin require that a coordinate change be introduced. As such, the following coordinate change is evoked, which converts the polar coordinates ða; cþ into Cartesian coordinates ð; Þ: ðtþ ¼aðtÞ cos½cðtþš, (24) ðtþ ¼aðtÞ sin½cðtþš. (25) This coordinate change results in averaged equations of the form ¼ 1 8 e½ð2l 1 4n 1 þ 4sÞ 3g 3 2 þð 3g 3 þ 2l 3 Þ 3 8zŠþOðe 2 Þ, (26) ¼ 1 8 e½ð2l 1 þ 4n 1 4sÞ þ 3g 3 2 þð3g 3 þ 2l 3 Þ 3 8zŠþOðe 2 Þ. (27) Using these equations, with z ¼, local stabilit results are inferred b considering the response of the sstem when linearized about the stead states. Letting the state vector be " XðtÞ ¼ ðtþ # (28) ðtþ and the stead-state values be given b X ¼ " #, (29) the local linearized equation of motion for the sstem can be written as Y ðtþ ¼Jj X YðtÞ, (3) where YðtÞ ¼XðtÞ X (31) and J is the Jacobian matri of the averaged equations, Eqs. (26) and (27), evaluated at stead state, X. The stabilit of the stead-state response is governed b the eigenvalues of the corresponding Jacobian. To simplif the analsis, we will consider the stabilit in terms of the trace and determinant of the 2 2 Jacobian [2]. In particular, the eigenvalues can be epressed in terms of the trace, T, and the determinant, D, as pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi L 1;2 ¼ 1 2 ðt T 2 4DÞ. (32) For the undamped case under consideration here, the trace of the Jacobian for each of the stead-state responses is zero, thus onl two generic equilibrium tpes are possible for nontrivial eigenvalues: saddles (unstable) and centers (marginall stable). Which of these equilibrium tpes eists depends solel on the sign of the determinant. In particular, when D4 the equilibrium is a center, and when Do the equilibrium is a saddle. The remaining case, D ¼, corresponds to the degenerate case of two identicall zero eigenvalues, and as such is used here onl to identif where stabilit changes occur. The trivial solution features a determinant of the form D ¼ 1 16 e2 ½l 2 1 4ðn 1 sþ 2 Š. (33) B finding the roots of D ¼, the critical frequenc values (detuning values) at which stabilit changes occur can be shown to be s 1 ¼ n 1 l 1 2 (34)
8 84 ARTICLE IN PRESS J.F. Rhoads et al. / Journal of Sound and Vibration 296 (26) and s 2 ¼ n 1 þ l 1 2, (35) respectivel. The stabilit in each frequenc and detuning region is then determined b evaluating the determinant at a single point within the domain or b calculating the derivative of the determinant at the critical frequenc value. This shows that the trivial solution appears as a center for so minðs 1 ; s 2 Þ and s4 maðs 1 ; s 2 Þ, and as a saddle for minðs 1 ; s 2 Þoso maðs 1 ; s 2 Þ. These results are, of course, unaffected b the nonlinear ecitation term and are closel related to the results obtained for a standard Mathieu equation. The distinct difference here, however, is that wedge of instabilit which appears in the l 1 O parameter space is generall distorted awa from the tpical configuration seen in Refs. [21,22], for eample, due to the presence of the frequenc-shifting n 1 term (this is considered etensivel in Ref. [1]). Using similar techniques, the constant amplitude solutions of Eq. (22) can be shown to be saddles for all frequenc values where the solutions eist, as determined b considering parameter values where their phases are real-valued. It is found that these solutions eist between the detuning values and s ¼ 3l 1g 3 l 1 l 3 þ 2l 3 n 1 2l 3 (36) s ¼ 3l 1g 3 þ l 1 l 3 þ 2l 3 n 1. (37) 2l 3 Note that each of these values ma serve as an upper or lower bound, depending on the sign of l 1. The remaining two branches are slightl more complicated in that their stabilit is heavil dependent on which region of the g 3 l 3 parameter space (see Fig. 2) the sstem falls within. Regardless, the stabilit for each of the branches can be found b evoking the fact that the response branch presented in Eq. (16) has a determinant given b D ¼ e2 ðl 1 2n 1 þ 2sÞð3l 1 g 3 l 1 l 3 2l 3 n 1 þ 2l 3 sþ, (38) 12g 3 8l 3 which features critical frequenc values of s ¼ n 1 l 1 2, (39) corresponding to the bifurcation at which this branch emanates from the origin, and s ¼ 3l 1g 3 þ l 1 l 3 þ 2l 3 n 1, (4) 2l 3 corresponding to the bifurcation which connects this branch to the constant amplitude branches. The response branch presented in Eq. (18) has a determinant given b D ¼ e2 ðl 1 þ 2n 1 2sÞð3l 1 g 3 þ l 1 l 3 2l 3 n 1 þ 2l 3 sþ, (41) 12g 3 þ 8l 3 which features critical frequenc values of s ¼ n 1 þ l 1 2, (42) corresponding to the bifurcation at which this branch emanates from the origin, and s ¼ 3l 1g 3 l 1 l 3 þ 2l 3 n 1, (43) 2l 3 corresponding to the bifurcation which connects this branch to the constant amplitude branches. Combining these results with the stead-state solutions derived above ields a complete picture of the oscillator s possible frequenc responses, as described in detail in Section 5.
9 J.F. Rhoads et al. / Journal of Sound and Vibration 296 (26) A B C 5 4 a Fig. 3. Representative frequenc response plot showing response amplitude versus frequenc detuning for Region I (g 3 ¼ :5, l 3 ¼ :5) in Fig. 2. Note the response s hardening nonlinear characteristics. Also note that here, as well as in Figs. 4 14, solid lines indicate stable response branches and dashed lines unstable response branches. 5. Sstem frequenc responses To facilitate understanding of the analsis presented in Section 4, consider a representative sdof dnamical sstem modeled b Eq. (1). In particular, consider an undamped version (z ¼ ) with n 1 ¼ l 1 ¼ 1. (The results can easil be etended for other values of n 1 and l 1, both positive and negative.) B appling the results of the aforementioned section to this tpical case, frequenc response plots, such as those shown in Figs. 3 14, are easil generated for each of the regions depicted in Fig. 2. A brief eplanation of each of these plots follows. To begin, consider the amplitude versus frequenc response curves depicted in Figs. 3 and 6, corresponding to Regions I and III in Fig. 2, respectivel. As a quick eamination reveals, these regions appear consistent with normal softening and hardening nonlinear behavior. That is, the nontrivial solutions branch off in two distinct pitchfork bifurcations, one subcritical and the other supercritical, and all solutions remain globall bounded. (Note that the term globall bounded is used here to designate the fact that all initial conditions result in bounded responses.) Though not shown, the phase response will be equall mundane in that it remains constant across the branches, as dictated b Eqs. (17) and (19). The frequenc responses shown in Figs. 4 and 5, corresponding to the topologicall equivalent regions labeled IIa and IIb in Fig. 2, are slightl more interesting. 1 Here the nontrivial solutions also branch off, but in each instance a subcritical bifurcation occurs. The net result is that the two response curves actuall bend awa from one another and some solutions are globall unbounded (specificall, all nonzero initial conditions lead to unbounded motions over the frequenc range where the trivial solution is unstable, and where the trivial solution is stable, initial conditions with sufficientl large energies lead to unbounded responses). Of course, here, and elsewhere, these globall unbounded solutions will be bounded b higher-order nonlinear effects not included here. The phase responses of the solutions within these regions, much like those considered above, remain invariant in frequenc and thus are omitted. 1 Note that these regions are given different designations onl due to the fact that the nontrivial response branches to the left and right of the resonance have different amplitudes, which are reversed in these two regions. This is not phsicall important, since both branches are unstable. However, in other regions, this distinction will have consequences in terms of the observed frequenc response behavior.
10 86 ARTICLE IN PRESS J.F. Rhoads et al. / Journal of Sound and Vibration 296 (26) D E F a Fig. 4. Representative frequenc response plot showing response amplitude versus frequenc detuning for Region IIa (g 3 ¼ :5, l 3 ¼ :1) in Fig. 2. Note the response s mied nonlinear characteristics G H I a Fig. 5. Representative frequenc response plot showing response amplitude versus frequenc detuning for Region IIb (g 3 ¼ :5, l 3 ¼ :1) in Fig. 2. Note the response s mied nonlinear characteristics. In the lower half-plane of the parameter space depicted in Fig. 2, the eistence of the additional nontrivial, constant amplitude, solutions complicates matters. First, consider Figs. 7 and 13, corresponding to Regions IV and VI in Fig. 2, respectivel. Here the near zero amplitude responses resemble those of Figs. 3 and 6, respectivel. However, as the frequenc is swept further awa from zero detuning (s ¼ ), the sstem undergoes two additional local bifurcations, corresponding to the creation and annihilation of the constant amplitude response, which connects the other two nontrivial branches. Though this still results in a globall
11 J.F. Rhoads et al. / Journal of Sound and Vibration 296 (26) J K L 5 4 a Fig. 6. Representative frequenc response plot showing response amplitude versus frequenc detuning for Region III (g 3 ¼ :5, l 3 ¼ :5) in Fig. 2. Note the response s softening nonlinear characteristics C M J K L 3 25 a Fig. 7. Representative frequenc response plot showing response amplitude versus frequenc detuning for Region IV (g 3 ¼ :5, l 3 ¼ :5) in Fig. 2. Note the response s quasi-softening nonlinear characteristics. bounded sstem with quasi-softening or hardening characteristics, it is distinct from the usual hardening and softening response curves of the upper half-plane, due to the fact that the stabilit in the two nonconstant branches switches via the constant amplitude branch. This switching can be further reconciled b considering the phase response plots shown in Figs. 8 and 14, which show the transition from one phase to another as the constant amplitude branch is traversed. The responses shown in Figs. 9 and 11, corresponding to the regions labeled Va and Vb in Fig. 2, are also quite interesting. Again, the response features four distinct local bifurcations, two of the pitchfork variet
12 88 ARTICLE IN PRESS J.F. Rhoads et al. / Journal of Sound and Vibration 296 (26) C M J K L ψ (rad) Fig. 8. Representative frequenc response plot showing phase angle versus frequenc detuning for Region IV (g 3 ¼ :5, l 3 ¼ :5) in Fig F N Q D 5 O P 4 a Fig. 9. Representative frequenc response plot showing response amplitude versus frequenc detuning for Region Va (g 3 ¼ :5, l 3 ¼ :1) in Fig. 2. Note the response s mied nonlinear characteristics. where the nontrivial solutions branch off from zero, and two corresponding to the creation and annihilation of the constant amplitude solution, again connecting the other two nontrivial branches. The net result is a globall unbounded sstem (that is, sufficientl large initial energies will lead to unbounded responses) with nontrivial response branches which bend towards one another and are both initiall stable. Both of these branches lose stabilit via bifurcations where the connect with the constant amplitude branch. It should also be noted that a global bifurcation, in particular, a saddle connection, occurs within these two regimes as well.
13 J.F. Rhoads et al. / Journal of Sound and Vibration 296 (26) F N Q D O P ψ (rad) Fig. 1. Representative frequenc response plot showing phase angle versus frequenc detuning for Region Va (g 3 ¼ :5, l 3 ¼ :1) in Fig I R U G 5 S T 4 a Fig. 11. Representative frequenc response plot showing response amplitude versus frequenc detuning for Region Vb (g 3 ¼ :5, l 3 ¼ :1) in Fig. 2. Note the response s mied nonlinear characteristics. This can be verified geometricall b considering the phase portraits presented in Fig. 15 below, which clearl show the connection of the manifolds originating from the saddle at the origin with the manifolds originating from the saddles corresponding to the constant amplitude solutions. This can also be confirmed analticall b performing a simple invariant manifold calculation. In particular, it is observed that when the saddle connections eist, the are linear manifolds, and the attendant and dnamics of the sstem at the specified
14 81 ARTICLE IN PRESS J.F. Rhoads et al. / Journal of Sound and Vibration 296 (26) I R U G S T ψ (rad) Fig. 12. Representative frequenc response plot showing phase angle versus frequenc detuning for Region Vb (g 3 ¼ :5, l 3 ¼ :1) in Fig A B C V J a Fig. 13. Representative frequenc response plot showing response amplitude versus frequenc detuning for Region VI (g 3 ¼ :5, l 3 ¼ :5) in Fig. 2. Note the response s quasi-hardening nonlinear characteristics. bifurcation point can be related according to ¼ a (44) and ¼ a, (45)
15 J.F. Rhoads et al. / Journal of Sound and Vibration 296 (26) A B C V J ψ (rad) Fig. 14. Representative frequenc response plot showing phase angle versus frequenc detuning for Region VI (g 3 ¼ :5, l 3 ¼ :5) in Fig (a) (b) Fig. 15. Representative phase planes corresponding to the frequenc (detuning) boundaries (a) OP and (b) ST where the saddle connection global bifurcations occur. where a is a to be determined constant that represents the slopes of the saddle-connection invariant manifolds. Substituting these relationships into the undamped versions of Eqs. (26) and (27) and equating like powers of (or ) ields both a constraint on the magnitude of a and a frequenc, or more accuratel, a detuning condition, at which the bifurcation occurs. In particular, it can be shown that the global bifurcation occurs for s ¼ 3g 3l 1 þ 4l 3 n 1, (46) 4l 3 which interestingl enough also corresponds to the frequenc value at which the two stable nontrivial branches intersect in the specified amplitude versus frequenc response plots (for Regions Va and Vb). It should be noted that the phase dnamics within these regions follow the general trends outlined above for Regions IV and VI.
16 812 ARTICLE IN PRESS J.F. Rhoads et al. / Journal of Sound and Vibration 296 (26) (a) (b) (c) (d) (e) (f) Fig. 16. Representative phase planes corresponding to (a) Region A, (b) Region B, (c) Region C, (d) Region D, (e) Region E, and (f) Region F in Figs To provide additional insight into the bifurcations described above, representative phase planes corresponding to each of the frequenc response regimes delineated in Figs are presented. These plots, corresponding to the conservative sstem eplored in this section, are shown in
17 J.F. Rhoads et al. / Journal of Sound and Vibration 296 (26) (a) (b) (c) (d) (e) (f) Fig. 17. Representative phase planes corresponding to (a) Region G, (b) Region H, (c) Region I, (d) Region J, (e) Region K, and (f) Region L in Figs Figs Note that the capital letters designate topologicall distinct phase portraits, and that these are associated with the similarl labeled frequenc ranges in the frequenc response plots shown in Figs
18 814 ARTICLE IN PRESS J.F. Rhoads et al. / Journal of Sound and Vibration 296 (26) (a) (b) (c) (d) (e) (f) Fig. 18. Representative phase planes corresponding to (a) Region M, (b) Region N, (c) Region O, (d) Region P, (e) Region Q, and (f) Region R in Figs While the results presented above are for a special case of the equation of motion in question, it should be noted, as previousl mentioned, that the results are easil etendible for arbitrar n 1 and l 1. In particular, other values of these parameters tpicall result in two notable differences. First, selecting different values for
19 J.F. Rhoads et al. / Journal of Sound and Vibration 296 (26) (a) (b) (c) (d) Fig. 19. Representative phase planes corresponding to (a) Region S, (b) Region T, (c) Region U, and (d) Region V in Figs n 1 and l 1 results in a distortion of the Arnold tongue associated with the stabilit of the solution s trivial solution. This is easil reconciled through consideration of Eqs. (34) and (35). In addition, specifing the value of l 1 to be negative results in a smmetric reflection of the g 3 l 3 parameter space about the l 3 ¼ ais. This can be reconciled b noting that the eistence of a 3 depends on the sign of the ratio l 1 =l 3, and thus for l 1 o the constant amplitude solutions will appear in the upper half-plane of the parameter space instead of the lower half-plane. It should also be noted that, in addition to the phenomena described above, the sstem ehibits a number of additional interesting features. First, note that in the absence of damping the phase portraits show a distinct smmetr. The breaking of this smmetr is considered in the following section, where the effects of small sstem damping on the sstem response are considered. Also, and related to the smmetr, over certain frequenc (detuning) ranges, an invariant circle eists, given b 2 þ 2 ¼ 2l 1. (47) l 3 Note that the radius of this circle corresponds to the magnitude of the constant amplitude branches and, therefore, these solutions lie on this invariant circle. Furthermore, these solutions are connected to one another via invariant manifolds consisting of circular arcs. These analtical results, and the graphical representations epressed in terms of frequenc response curves and phase portraits, give a complete picture of the response of the sstem for zero damping. Small damping will destro the smmetr that permits such a thorough analsis, but the resulting responses can be inferred b
20 816 J.F. Rhoads et al. / Journal of Sound and Vibration 296 (26) straightforward etrapolations of the undamped results. Steps in this direction, including results from numerical analsis of the damped equations of motion, are considered in the net section. 6. The effects of damping Thus far, our analsis of Eq. (1) has considered the case of zero damping (z ¼ ). To gain a better understanding of the full oscillator dnamics, we now present numerical bifurcation results for the sstem with nonzero damping and with the detuning (s) treated as a bifurcation parameter. These results are obtained using AUTO [23], which requires that we rewrite the non-autonomous Eq. (1) as a sstem of autonomous first-order equations. Rewriting this as first-order equations is accomplished b defining v ¼ z and considering equations for z and v. We make this sstem autonomous b augmenting Eq. (1) with equations for a nonlinear oscillator with a stable periodic orbit with one variable equal to cos Ot, and replacing the appropriate terms of Eq. (1) with this variable. Most of our attention will be on Region IV from Fig. 2, as this displas the most interesting bifurcation behavior. We will also present tpical results for Regions IIb, III, and Va; Regions I, IIa, Vb, and VI are related to the presented regions b smmetr. As a tpical point in Region IV, we take l 1 ¼ n 1 ¼ 1 and g 3 ¼ l 3 ¼ :5, as in Fig. 7. Throughout this section, we take ¼ :1. The no-motion state (z ¼ z ¼ ) eists for all values of z and s, and undergoes bifurcations giving rise to branches of periodic orbits which ma be followed with AUTO. These periodic orbits displa response at half the driving frequenc O, and have the smmetr propert ðz; z Þðt þ T=2Þ ¼ð z; z ÞðtÞ, (48) where T ¼ 2p=O is the period of the response. Alread for small damping we find the interesting result that the constant amplitude branch shown in Fig. 7 splits such that the branches with the distinct phases now have slightl different amplitudes; see Fig. 2 for results with z ¼ :1. Here the stable periodic orbit branch which bifurcates from the trivial solution at s 1:5 undergoes a saddle-node bifurcation at s 1 before returning to the trivial state in the bifurcation at s :5. The branch containing the stable periodic orbit for s ¼ 1 undergoes a saddle-node bifurcation at s, and is disconnected from the no-motion state. Periodic orbits on this disconnected branch also displa the smmetr given in Eq. (48) z ma Fig. 2. Bifurcation diagram showing the maimum value of the displacement (z ma ) for z ¼ :1 and other parameters as described in the tet. In this and other bifurcation diagrams, diamonds, squares, and triangles indicate crossing the boundar of an Arnold tongue, period-doubling, and smmetr-breaking pitchfork bifurcations, respectivel. Furthermore, solid and dashed lines represent stable and unstable solutions, respectivel.
21 J.F. Rhoads et al. / Journal of Sound and Vibration 296 (26) z ma Fig. 21. Bifurcation diagram for z ¼ :1 and other parameters corresponding to Region IV, as described in the tet. For larger damping, we find additional bifurcations on the disconnected branch; see Fig. 21 for z ¼ :1 (a ver large value for a MEM resonator). A periodic orbit with the smmetr of Eq. (48) is stable from :57oso 3:739; a tpical orbit is shown in Fig. 22(a). At the upper limit it undergoes a saddle-node bifurcation, while at the lower limit it undergoes a smmetr-breaking pitchfork bifurcation. The smmetrbreaking bifurcation leads to two smmetr-related asmmetric periodic orbits, which lack the smmetr given in Eq. (48); see Fig. 22(b) for an eample of one of these (the other is obtained b letting ðz; z Þ!ð z; z Þ). Each of these smmetr-related asmmetric periodic orbits undergoes a period-doubling bifurcation at s ¼ 3:884; one of the period-doubled orbits is shown in Fig. 22(c). This is the beginning of a period-doubling cascade to chaos as s decreases, with chaotic attractor as shown in Fig. 22(d). This is one of two smmetr-related chaotic attractors (the other is obtained b letting ðz; z Þ!ð z; z Þ). We note that the basin of attraction for this chaotic attractor is quite small, making it somewhat difficult to locate (and likel unobservable in phsical eperiments). We also note that from [24], it is necessar for a smmetric periodic orbit such as that shown in Fig. 22(a) to first undergo a smmetr-breaking bifurcation before a perioddoubling cascade can occur. Fig. 23 shows results from a two-parameter (s and z) bifurcation stud for the Region IV parameters used above. The no-motion state is unstable inside the dashed curve. The width of this instabilit region decreases as the damping z increases, consistent with the effect of damping for the related Mathieu equation [21]. The dot-dashed lines show bifurcation sets for the saddle-node bifurcations of the smmetric periodic orbits. The crossing of these at ðs; zþ ð :3; :5Þ does not correspond to a codimension-two bifurcation; the two saddlenode bifurcations are just (independentl) occurring at the same parameter values. One dash-dotted curve (corresponding to a saddle-node bifurcation of the periodic orbit branch which bifurcates from the no-motion branch) terminates on the dashed curve at ðs; zþ ð:796; :222Þ. Here the periodic orbit branch emerging from the no-motion branch bifurcates verticall. For smaller z values, there is an unstable periodic orbit branch bifurcating to smaller s values (which then turns around in a saddle-node bifurcation), while for larger z values there is a stable branch which bifurcates to larger s values. We now consider the point in Region Va with l 1 ¼ n 1 ¼ 1 and g 3 ¼ :5; l 3 ¼ :1, as in Fig. 9. For small z we similarl find that the constant amplitude solution branch shown in Fig. 9 splits, giving a periodic orbit branch which is disconnected from the no-motion state; see Fig. 24 for z ¼ :1. This disconnected branch persists for larger z; see Fig. 25 for z ¼ :1, for which we do not find an further bifurcations. Note that we obtained unstable periodic orbits on the disconnected branch using a shooting method, and these were then followed with AUTO.
22 818 J.F. Rhoads et al. / Journal of Sound and Vibration 296 (26) z z -3-4 z (a) z τ z z z (b) z τ z z z -3-4 (c) z τ z z (d) z z Fig. 22. Stable solutions for parameters as in Fig. 21 and with forcing cos Ot: (a) smmetric periodic orbit for s ¼ 1:5; (b) asmmetric periodic orbit for s ¼ 3:85; (c) period doubled orbit for s ¼ 3:895; (d) chaotic solution for s ¼ 3:9235, with right panel zoomed in. The timeseries in (a), (b), and (c) show one period.
23 J.F. Rhoads et al. / Journal of Sound and Vibration 296 (26) ζ Fig. 23. Bifurcation sets for the parameters corresponding to Region IV. The no-motion state is unstable inside the dashed curve, while the dot-dashed lines show bifurcation sets for the saddle-node bifurcations of the smmetric periodic orbits. z ma Fig. 24. Bifurcation diagram for z ¼ :1 and other parameters corresponding to Region Va, as described in the tet. A similar period-doubling cascade to chaos, as found above for Region IV, is present for z ¼ :1 for the point in Region III with l 1 ¼ n 1 ¼ 1 and g 3 ¼ :5; l 3 ¼ :5, as in Fig. 6; see Fig. 26. In this region, there are no constant amplitude solutions for z ¼, as there were for Region IV. We close this section b noting that no interesting bifurcation behavior is found for z ¼ :1 for the point in Region IIb with l 1 ¼ n 1 ¼ 1andg 3 ¼ :5; l 3 ¼ :1, as in Fig. 5.
24 82 ARTICLE IN PRESS J.F. Rhoads et al. / Journal of Sound and Vibration 296 (26) z ma Fig. 25. Bifurcation diagram for z ¼ :1 and other parameters corresponding to Region Va, as described in the tet z ma Fig. 26. Bifurcation diagram for z ¼ :1 and other parameters corresponding to Region III, as described in the tet. 7. Application to MEM resonators The results developed in the previous sections are directl applicable to the equations of motion for parametricall ecited MEM oscillators of the tpe developed in Section 2. However, in this case, the coefficients in the equation of motion that are related to the electrostatic forces depend on the amplitude of the AC voltage, V A. Accordingl, the qualitative nature of both the linear and nonlinear behavior of the oscillator depend on V A. The effect of varing V A on the linear behavior of the sstem is easil noted b considering the undamped MEM oscillator s linear stabilit chart, as presented in Fig. 27. As clearl shown, as V A is increased, the width (in terms of detuning) of the linear instabilit zone increases. As noted in Section 4, this is largel consistent
25 J.F. Rhoads et al. / Journal of Sound and Vibration 296 (26) Stable Unstable Stable 25 V A (V) f (Hz) Fig. 27. Linear stabilit chart in the AC ecitation voltage amplitude (V A ) versus ecitation frequenc (f) parameter space for an undamped MEM oscillator with sstem parameters as specified in Table 2. Eperimentall derived instabilit points are also shown. Note that in each case the noted stabilit corresponds to the trivial solution. with the behavior of a tpical Mathieu sstem, but with two notable eceptions. First, due to the nonlinear relationship between l 1 and V A the boundaries of the instabilit feature a distinct curvature. In addition, due the presence of the frequenc shift caused b the electrostatic stiffness term (l 1 ), the instabilit region is distorted from its nominall smmetric configuration. Each of these differences is eamined etensivel in Ref. [1]. From a practical point of view, it is also important to note that in the presence of damping, the bottom of the instabilit wedge in Fig. 27 pulls up from zero voltage and becomes rounded off, as shown eperimentall in both Fig. 27 and Ref. [15]. The effect of V A on the nonlinear behavior of the sstem can be realized b noting that as the ecitation amplitude is varied, the sstem can transition from one response regime to another (in Fig. 2). This can be confirmed b eamining the form of the nonlinearit in the g 3 l 3 parameter space. Taking the definitions of l 3 and g 3 from Table 1 it is seen that at zero voltage the sstem starts at ðg 3 ; l 3 Þ¼ðw; Þ (the point corresponding to the purel mechanical hardening nonlinearit), which lies on the boundar between Regions I and VI, and moves along a straight line in the g 3 l 3 parameter space given b l 3 ¼ g 3 w (49) as V A is increased [1,5]. To validate this, consider the oscillator design described b the parameters in Table 2, which were obtained, as subsequentl described, from a device similar to that shown in Fig. 1. AsFig. 28 shows, the nonlinearit for this device will move through regions VI, Vb, Va, and IV as the amplitude of the AC voltage is varied. Thus, the sstem transitions from a quasi-hardening nonlinearit, to a mied nonlinearit, and finall to a quasi-softening nonlinearit as the oscillator s input voltage increases. Figs show response curves for the eperimental parameter values for three different values of V A.Note that the response indeed transitions from hardening, to mied, to softening, as predicted. With the aforementioned behavior in mind, a designer can account for these transitions and estimate whether or not the will occur in a given device s operating voltage range. The following relationships, derived b setting Z 1 and Z 2 equal to zero, provide epressions for the transition voltages, that is, the critical voltage values corresponding to qualitative changes in the sstem s nonlinearit. For a given oscillator, sffiffiffiffiffiffiffiffiffiffiffi 3k 3 V A;C1 ¼ (5) 5r 3A
520 Chapter 9. Nonlinear Differential Equations and Stability. dt =
5 Chapter 9. Nonlinear Differential Equations and Stabilit dt L dθ. g cos θ cos α Wh was the negative square root chosen in the last equation? (b) If T is the natural period of oscillation, derive the
More informationVibrational Power Flow Considerations Arising From Multi-Dimensional Isolators. Abstract
Vibrational Power Flow Considerations Arising From Multi-Dimensional Isolators Rajendra Singh and Seungbo Kim The Ohio State Universit Columbus, OH 4321-117, USA Abstract Much of the vibration isolation
More informationCONTINUOUS SPATIAL DATA ANALYSIS
CONTINUOUS SPATIAL DATA ANALSIS 1. Overview of Spatial Stochastic Processes The ke difference between continuous spatial data and point patterns is that there is now assumed to be a meaningful value, s
More informationNew approach to study the van der Pol equation for large damping
Electronic Journal of Qualitative Theor of Differential Equations 2018, No. 8, 1 10; https://doi.org/10.1422/ejqtde.2018.1.8 www.math.u-szeged.hu/ejqtde/ New approach to stud the van der Pol equation for
More information10 Back to planar nonlinear systems
10 Back to planar nonlinear sstems 10.1 Near the equilibria Recall that I started talking about the Lotka Volterra model as a motivation to stud sstems of two first order autonomous equations of the form
More information* τσ σκ. Supporting Text. A. Stability Analysis of System 2
Supporting Tet A. Stabilit Analsis of Sstem In this Appendi, we stud the stabilit of the equilibria of sstem. If we redefine the sstem as, T when -, then there are at most three equilibria: E,, E κ -,,
More informationStrain Transformation and Rosette Gage Theory
Strain Transformation and Rosette Gage Theor It is often desired to measure the full state of strain on the surface of a part, that is to measure not onl the two etensional strains, and, but also the shear
More informationUNFOLDING THE CUSP-CUSP BIFURCATION OF PLANAR ENDOMORPHISMS
UNFOLDING THE CUSP-CUSP BIFURCATION OF PLANAR ENDOMORPHISMS BERND KRAUSKOPF, HINKE M OSINGA, AND BRUCE B PECKHAM Abstract In man applications of practical interest, for eample, in control theor, economics,
More informationA route to computational chaos revisited: noninvertibility and the breakup of an invariant circle
A route to computational chaos revisited: noninvertibilit and the breakup of an invariant circle Christos E. Frouzakis Combustion Research Laborator, Paul Scherrer Institute CH-5232, Villigen, Switzerland
More informationFinding Limits Graphically and Numerically. An Introduction to Limits
60_00.qd //0 :05 PM Page 8 8 CHAPTER Limits and Their Properties Section. Finding Limits Graphicall and Numericall Estimate a it using a numerical or graphical approach. Learn different was that a it can
More informationFinding Limits Graphically and Numerically. An Introduction to Limits
8 CHAPTER Limits and Their Properties Section Finding Limits Graphicall and Numericall Estimate a it using a numerical or graphical approach Learn different was that a it can fail to eist Stud and use
More informationNONLINEAR DYNAMICS AND CHAOS. Numerical integration. Stability analysis
LECTURE 3: FLOWS NONLINEAR DYNAMICS AND CHAOS Patrick E McSharr Sstems Analsis, Modelling & Prediction Group www.eng.o.ac.uk/samp patrick@mcsharr.net Tel: +44 83 74 Numerical integration Stabilit analsis
More informationAdditional Topics in Differential Equations
6 Additional Topics in Differential Equations 6. Eact First-Order Equations 6. Second-Order Homogeneous Linear Equations 6.3 Second-Order Nonhomogeneous Linear Equations 6.4 Series Solutions of Differential
More informationC. Non-linear Difference and Differential Equations: Linearization and Phase Diagram Technique
C. Non-linear Difference and Differential Equations: Linearization and Phase Diaram Technique So far we have discussed methods of solvin linear difference and differential equations. Let us now discuss
More information10.2 The Unit Circle: Cosine and Sine
0. The Unit Circle: Cosine and Sine 77 0. The Unit Circle: Cosine and Sine In Section 0.., we introduced circular motion and derived a formula which describes the linear velocit of an object moving on
More informationA new chaotic attractor from general Lorenz system family and its electronic experimental implementation
Turk J Elec Eng & Comp Sci, Vol.18, No.2, 2010, c TÜBİTAK doi:10.3906/elk-0906-67 A new chaotic attractor from general Loren sstem famil and its electronic eperimental implementation İhsan PEHLİVAN, Yılma
More informationAnalytic Geometry in Three Dimensions
Analtic Geometr in Three Dimensions. The Three-Dimensional Coordinate Sstem. Vectors in Space. The Cross Product of Two Vectors. Lines and Planes in Space The three-dimensional coordinate sstem is used
More informationChapter 4 Analytic Trigonometry
Analtic Trigonometr Chapter Analtic Trigonometr Inverse Trigonometric Functions The trigonometric functions act as an operator on the variable (angle, resulting in an output value Suppose this process
More informationSection 1.5 Formal definitions of limits
Section.5 Formal definitions of limits (3/908) Overview: The definitions of the various tpes of limits in previous sections involve phrases such as arbitraril close, sufficientl close, arbitraril large,
More informationDynamics of a mass-spring-pendulum system with vastly different frequencies
Dynamics of a mass-spring-pendulum system with vastly different frequencies Hiba Sheheitli, hs497@cornell.edu Richard H. Rand, rhr2@cornell.edu Cornell University, Ithaca, NY, USA Abstract. We investigate
More informationAdditional Topics in Differential Equations
0537_cop6.qd 0/8/08 :6 PM Page 3 6 Additional Topics in Differential Equations In Chapter 6, ou studied differential equations. In this chapter, ou will learn additional techniques for solving differential
More informationStudy on Nonlinear Vibration and Crack Fault of Rotor-bearing-seal Coupling System
Sensors & Transducers 04 b IFSA Publishing S. L. http://www.sensorsportal.com Stud on Nonlinear Vibration and Crack Fault of Rotor-bearing-seal Coupling Sstem Yuegang LUO Songhe ZHANG Bin WU Wanlei WANG
More informationMethods of Solving Ordinary Differential Equations (Online)
7in 0in Felder c0_online.te V3 - Januar, 05 0:5 A.M. Page CHAPTER 0 Methods of Solving Ordinar Differential Equations (Online) 0.3 Phase Portraits Just as a slope field (Section.4) gives us a wa to visualize
More informationExperimental Calibration and Head Loss Prediction of Tuned Liquid Column Damper
Tamkang Journal of Science and Engineering, Vol. 8, No 4, pp. 319-35 (005) 319 Experimental Calibration and Head Loss Prediction of Tuned Liquid Column Damper Jong-Cheng Wu Department of Civil Engineering,
More informationLATERAL BUCKLING ANALYSIS OF ANGLED FRAMES WITH THIN-WALLED I-BEAMS
Journal of arine Science and J.-D. Technolog, Yau: ateral Vol. Buckling 17, No. Analsis 1, pp. 9-33 of Angled (009) Frames with Thin-Walled I-Beams 9 ATERA BUCKING ANAYSIS OF ANGED FRAES WITH THIN-WAED
More informationEstimation Of Linearised Fluid Film Coefficients In A Rotor Bearing System Subjected To Random Excitation
Estimation Of Linearised Fluid Film Coefficients In A Rotor Bearing Sstem Subjected To Random Ecitation Arshad. Khan and Ahmad A. Khan Department of Mechanical Engineering Z.. College of Engineering &
More informationSecond-Order Linear Differential Equations C 2
C8 APPENDIX C Additional Topics in Differential Equations APPENDIX C. Second-Order Homogeneous Linear Equations Second-Order Linear Differential Equations Higher-Order Linear Differential Equations Application
More informationRECENTLY, microelectromechanical oscillators exploiting
310 JOURNAL OF MICROELECTROMECHANICAL SYSTEMS, VOL. 16, NO. 2, APRIL 2007 Linear and Nonlinear Tuning of Parametrically Excited MEMS Oscillators Barry E. DeMartini, Student Member, IEEE, Jeffrey F. Rhoads,
More informationFunctions of Several Variables
Chapter 1 Functions of Several Variables 1.1 Introduction A real valued function of n variables is a function f : R, where the domain is a subset of R n. So: for each ( 1,,..., n ) in, the value of f is
More informationAE/ME 339. K. M. Isaac Professor of Aerospace Engineering. December 21, 2001 topic13_grid_generation 1
AE/ME 339 Professor of Aerospace Engineering December 21, 2001 topic13_grid_generation 1 The basic idea behind grid generation is the creation of the transformation laws between the phsical space and the
More informationNonlinear Systems Examples Sheet: Solutions
Nonlinear Sstems Eamples Sheet: Solutions Mark Cannon, Michaelmas Term 7 Equilibrium points. (a). Solving ẋ =sin 4 3 =for gives =as an equilibrium point. This is the onl equilibrium because there is onl
More informationTwinkling Phenomena In Snap-Through Oscillators
Twinkling Phenomena In Snap-Through Oscillators Brian F. Feeny Department of Mechanical Engineering Michigan State University East Lansing, MI 4884 6 Alejandro R. Diaz Department of Mechanical Engineering
More information10.5 Graphs of the Trigonometric Functions
790 Foundations of Trigonometr 0.5 Graphs of the Trigonometric Functions In this section, we return to our discussion of the circular (trigonometric functions as functions of real numbers and pick up where
More informationResearch Article Chaotic Attractor Generation via a Simple Linear Time-Varying System
Discrete Dnamics in Nature and Societ Volume, Article ID 836, 8 pages doi:.//836 Research Article Chaotic Attractor Generation via a Simple Linear Time-Varing Sstem Baiu Ou and Desheng Liu Department of
More informationAPPENDIX D Rotation and the General Second-Degree Equation
APPENDIX D Rotation and the General Second-Degree Equation Rotation of Aes Invariants Under Rotation After rotation of the - and -aes counterclockwise through an angle, the rotated aes are denoted as the
More information8.1 Exponents and Roots
Section 8. Eponents and Roots 75 8. Eponents and Roots Before defining the net famil of functions, the eponential functions, we will need to discuss eponent notation in detail. As we shall see, eponents
More informationSMALL bodies of the solar system, such as Phobos, cannot
INSTITUTO SUPERIOR TÉCNICO, LISBON, PORTUGAL, JUNE 013 1 Controlled Approach Strategies on Small Celestial Bodies Using Approimate Analtical Solutions of the Elliptical Three-Bod Problem: Application to
More informationy R T However, the calculations are easier, when carried out using the polar set of co-ordinates ϕ,r. The relations between the co-ordinates are:
Curved beams. Introduction Curved beams also called arches were invented about ears ago. he purpose was to form such a structure that would transfer loads, mainl the dead weight, to the ground b the elements
More informationSymmetry Arguments and the Role They Play in Using Gauss Law
Smmetr Arguments and the Role The la in Using Gauss Law K. M. Westerberg (9/2005) Smmetr plas a ver important role in science in general, and phsics in particular. Arguments based on smmetr can often simplif
More informationControl Schemes to Reduce Risk of Extinction in the Lotka-Volterra Predator-Prey Model
Journal of Applied Mathematics and Phsics, 2014, 2, 644-652 Published Online June 2014 in SciRes. http://www.scirp.org/journal/jamp http://d.doi.org/10.4236/jamp.2014.27071 Control Schemes to Reduce Risk
More information= F ( x; µ) (1) where x is a 2-dimensional vector, µ is a parameter, and F :
1 Bifurcations Richard Bertram Department of Mathematics and Programs in Neuroscience and Molecular Biophysics Florida State University Tallahassee, Florida 32306 A bifurcation is a qualitative change
More informationAdditive resonances of a controlled van der Pol-Duffing oscillator
Additive resonances of a controlled van der Pol-Duffing oscillator This paper has been published in Journal of Sound and Vibration vol. 5 issue - 8 pp.-. J.C. Ji N. Zhang Faculty of Engineering University
More informationSection 8.5 Parametric Equations
504 Chapter 8 Section 8.5 Parametric Equations Man shapes, even ones as simple as circles, cannot be represented as an equation where is a function of. Consider, for eample, the path a moon follows as
More informationMathematics 309 Conic sections and their applicationsn. Chapter 2. Quadric figures. ai,j x i x j + b i x i + c =0. 1. Coordinate changes
Mathematics 309 Conic sections and their applicationsn Chapter 2. Quadric figures In this chapter want to outline quickl how to decide what figure associated in 2D and 3D to quadratic equations look like.
More information5.6 RATIOnAl FUnCTIOnS. Using Arrow notation. learning ObjeCTIveS
CHAPTER PolNomiAl ANd rational functions learning ObjeCTIveS In this section, ou will: Use arrow notation. Solve applied problems involving rational functions. Find the domains of rational functions. Identif
More informationThe Hopf Bifurcation Theorem: Abstract. Introduction. Transversality condition; the eigenvalues cross the imginary axis with non-zero speed
Supercritical and Subcritical Hopf Bifurcations in Non Linear Maps Tarini Kumar Dutta, Department of Mathematics, Gauhati Universit Pramila Kumari Prajapati, Department of Mathematics, Gauhati Universit
More informationSurvey of Wave Types and Characteristics
Seminar: Vibrations and Structure-Borne Sound in Civil Engineering Theor and Applications Surve of Wave Tpes and Characteristics Xiuu Gao April 1 st, 2006 Abstract Mechanical waves are waves which propagate
More informationSolving Variable-Coefficient Fourth-Order ODEs with Polynomial Nonlinearity by Symmetric Homotopy Method
Applied and Computational Mathematics 218; 7(2): 58-7 http://www.sciencepublishinggroup.com/j/acm doi: 1.11648/j.acm.21872.14 ISSN: 2328-565 (Print); ISSN: 2328-5613 (Online) Solving Variable-Coefficient
More informationA simple feedback control for a chaotic oscillator with limited power supply
Journal of Sound and Vibration 299 (2007) 664 671 Short Communication A simple feedback control for a chaotic oscillator with limited power supply JOURNAL OF SOUND AND VIBRATION S.L.T. de Souza a, I.L.
More informationx y plane is the plane in which the stresses act, yy xy xy Figure 3.5.1: non-zero stress components acting in the x y plane
3.5 Plane Stress This section is concerned with a special two-dimensional state of stress called plane stress. It is important for two reasons: () it arises in real components (particularl in thin components
More informationEP225 Note No. 4 Wave Motion
EP5 Note No. 4 Wave Motion 4. Sinusoidal Waves, Wave Number Waves propagate in space in contrast to oscillations which are con ned in limited regions. In describing wave motion, spatial coordinates enter
More information14.1 Systems of Linear Equations in Two Variables
86 Chapter 1 Sstems of Equations and Matrices 1.1 Sstems of Linear Equations in Two Variables Use the method of substitution to solve sstems of equations in two variables. Use the method of elimination
More informationStability Analysis of a Geometrically Imperfect Structure using a Random Field Model
Stabilit Analsis of a Geometricall Imperfect Structure using a Random Field Model JAN VALEŠ, ZDENĚK KALA Department of Structural Mechanics Brno Universit of Technolog, Facult of Civil Engineering Veveří
More informationDynamics of Few Co-rotating Vortices in Bose-Einstein Condensates
Dnamics of Few Co-rotating Vortices in Bose-Einstein Condensates R. Navarro 1, R. Carretero-González 1, P.J. Torres, P.G. Kevrekidis 3, D.J. Frantzeskakis 4, M.W. Ra 5, E. Altuntaş 5, and D.S. Hall 5 1
More information3.7 InveRSe FUnCTIOnS
CHAPTER functions learning ObjeCTIveS In this section, ou will: Verif inverse functions. Determine the domain and range of an inverse function, and restrict the domain of a function to make it one-to-one.
More informationOn Range and Reflecting Functions About the Line y = mx
On Range and Reflecting Functions About the Line = m Scott J. Beslin Brian K. Heck Jerem J. Becnel Dept.of Mathematics and Dept. of Mathematics and Dept. of Mathematics and Computer Science Computer Science
More informationStructure and stability of two-dimensional Bose- Einstein condensates under both harmonic and lattice confinement
Universit of Massachusetts Amherst ScholarWorks@UMass Amherst Mathematics and Statistics Department Facult Publication Series Mathematics and Statistics 28 Structure and stabilit of two-dimensional Bose-
More informationIdentifying second degree equations
Chapter 7 Identifing second degree equations 71 The eigenvalue method In this section we appl eigenvalue methods to determine the geometrical nature of the second degree equation a 2 + 2h + b 2 + 2g +
More informationAnd similarly in the other directions, so the overall result is expressed compactly as,
SQEP Tutorial Session 5: T7S0 Relates to Knowledge & Skills.5,.8 Last Update: //3 Force on an element of area; Definition of principal stresses and strains; Definition of Tresca and Mises equivalent stresses;
More information1.2 Relations. 20 Relations and Functions
0 Relations and Functions. Relations From one point of view, all of Precalculus can be thought of as studing sets of points in the plane. With the Cartesian Plane now fresh in our memor we can discuss
More informationResearch Article A New Four-Scroll Chaotic Attractor Consisted of Two-Scroll Transient Chaotic and Two-Scroll Ultimate Chaotic
Mathematical Problems in Engineering Volume, Article ID 88, pages doi:.//88 Research Article A New Four-Scroll Chaotic Attractor Consisted of Two-Scroll Transient Chaotic and Two-Scroll Ultimate Chaotic
More informationExercise solutions: concepts from chapter 5
1) Stud the oöids depicted in Figure 1a and 1b. a) Assume that the thin sections of Figure 1 lie in a principal plane of the deformation. Measure and record the lengths and orientations of the principal
More informationLocal Phase Portrait of Nonlinear Systems Near Equilibria
Local Phase Portrait of Nonlinear Sstems Near Equilibria [1] Consider 1 = 6 1 1 3 1, = 3 1. ( ) (a) Find all equilibrium solutions of the sstem ( ). (b) For each equilibrium point, give the linear approimating
More information3 First order nonlinear equations
Separable equations 3 First order nonlinear equations d f dx = where f(x,) is not (in general) a linear function of. ( x, ) The equation is autonomous if there is no explicit dependence on the independent
More informationParabolas. Example. y = ax 2 + bx + c where a = 1, b = 0, c = 0. = x 2 + 6x [expanding] \ y = x 2 + 6x + 11 and so is of the form
Parabolas NCEA evel Achievement Standard 9157 (Mathematics and Statistics.) Appl graphical methods in solving problems Methods include: graphs at curriculum level 7, their features and their equations
More informationSection 1.2: Relations, from College Algebra: Corrected Edition by Carl Stitz, Ph.D. and Jeff Zeager, Ph.D. is available under a Creative Commons
Section.: Relations, from College Algebra: Corrected Edition b Carl Stitz, Ph.D. and Jeff Zeager, Ph.D. is available under a Creative Commons Attribution-NonCommercial-ShareAlike.0 license. 0, Carl Stitz.
More informationHÉNON HEILES HAMILTONIAN CHAOS IN 2 D
ABSTRACT HÉNON HEILES HAMILTONIAN CHAOS IN D MODELING CHAOS & COMPLEXITY 008 YOUVAL DAR PHYSICS DAR@PHYSICS.UCDAVIS.EDU Chaos in two degrees of freedom, demonstrated b using the Hénon Heiles Hamiltonian
More informationMath Notes on sections 7.8,9.1, and 9.3. Derivation of a solution in the repeated roots case: 3 4 A = 1 1. x =e t : + e t w 2.
Math 7 Notes on sections 7.8,9., and 9.3. Derivation of a solution in the repeated roots case We consider the eample = A where 3 4 A = The onl eigenvalue is = ; and there is onl one linearl independent
More informationPolynomial and Rational Functions
Polnomial and Rational Functions Figure -mm film, once the standard for capturing photographic images, has been made largel obsolete b digital photograph. (credit film : modification of work b Horia Varlan;
More informationSolutions to Two Interesting Problems
Solutions to Two Interesting Problems Save the Lemming On each square of an n n chessboard is an arrow pointing to one of its eight neighbors (or off the board, if it s an edge square). However, arrows
More informationLESSON #42 - INVERSES OF FUNCTIONS AND FUNCTION NOTATION PART 2 COMMON CORE ALGEBRA II
LESSON #4 - INVERSES OF FUNCTIONS AND FUNCTION NOTATION PART COMMON CORE ALGEBRA II You will recall from unit 1 that in order to find the inverse of a function, ou must switch and and solve for. Also,
More informationAdditional Material On Recursive Sequences
Penn State Altoona MATH 141 Additional Material On Recursive Sequences 1. Graphical Analsis Cobweb Diagrams Consider a generic recursive sequence { an+1 = f(a n ), n = 1,, 3,..., = Given initial value.
More informationDynamics of quantum vortices in a toroidal trap
PHYSICAL REVIEW A 79, 4362 29 Dnamics of quantum vortices in a toroidal trap Peter Mason and Natalia G. Berloff Department of Applied Mathematics and Theoretical Phsics, Universit of Cambridge, Wilberforce
More informationInterspecific Segregation and Phase Transition in a Lattice Ecosystem with Intraspecific Competition
Interspecific Segregation and Phase Transition in a Lattice Ecosstem with Intraspecific Competition K. Tainaka a, M. Kushida a, Y. Ito a and J. Yoshimura a,b,c a Department of Sstems Engineering, Shizuoka
More informationMAE 323: Chapter 4. Plane Stress and Plane Strain. The Stress Equilibrium Equation
The Stress Equilibrium Equation As we mentioned in Chapter 2, using the Galerkin formulation and a choice of shape functions, we can derive a discretized form of most differential equations. In Structural
More informationMULTIPLE BIFURCATIONS OF A CYLINDRICAL DYNAMICAL SYSTEM
Journal of Theoretical and Applied Mechanics, Sofia, 2016, vol 46, No 1, pp 33 52 MULTIPLE BIFURCATIONS OF A CYLINDRICAL DYNAMICAL SYSTEM Ning Han, Qingjie Cao Centre for Nonlinear Dnamics Research, Harbin
More information12.4. Curvature. Introduction. Prerequisites. Learning Outcomes
Curvature 12.4 Introduction Curvature is a measure of how sharpl a curve is turning. At a particular point along the curve a tangent line can be drawn; this tangent line making an angle ψ with the positive
More informationKINEMATIC RELATIONS IN DEFORMATION OF SOLIDS
Chapter 8 KINEMATIC RELATIONS IN DEFORMATION OF SOLIDS Figure 8.1: 195 196 CHAPTER 8. KINEMATIC RELATIONS IN DEFORMATION OF SOLIDS 8.1 Motivation In Chapter 3, the conservation of linear momentum for a
More informationNormalization and Zero-Point Energy The amplitude A 2 in Eq can be found from the normalizing equation, 1106 CHAPTER 39 MORE ABOUT MATTER WAVES
116 CHAPTER 39 MORE ABOUT MATTER WAVES Fig. 39-4 A dot plot of the radial probabilit densit P(r) for the hdrogen atom in a quantum state with a relativel large principal quantum number namel, n 45 and
More informationNUMERICAL MODELING OF THE DYNAMIC RESPONSE OF A MULTI- BILINEAR-SPRING SUPPORT SYSTEM
Universit of Kentuck UKnowledge Universit of Kentuck Master's Theses Graduate School NUMERICAL MODELING OF THE DYNAMIC RESPONSE OF A MULTI- BILINEAR-SPRING SUPPORT SYSTEM Tre D. Gilliam Universit of Kentuck,
More informationP 0 (x 0, y 0 ), P 1 (x 1, y 1 ), P 2 (x 2, y 2 ), P 3 (x 3, y 3 ) The goal is to determine a third degree polynomial of the form,
Bezier Curves While working for the Renault automobile compan in France, an engineer b the name of P. Bezier developed a sstem for designing car bodies based partl on some fairl straightforward mathematics.
More informationDifference Resonances in a controlled van der Pol-Duffing oscillator involving time. delay
Difference Resonances in a controlled van der Pol-Duffing oscillator involving time delay This paper was published in the journal Chaos, Solitions & Fractals, vol.4, no., pp.975-98, Oct 9 J.C. Ji, N. Zhang,
More information1.5. Analyzing Graphs of Functions. The Graph of a Function. What you should learn. Why you should learn it. 54 Chapter 1 Functions and Their Graphs
0_005.qd /7/05 8: AM Page 5 5 Chapter Functions and Their Graphs.5 Analzing Graphs of Functions What ou should learn Use the Vertical Line Test for functions. Find the zeros of functions. Determine intervals
More informationJones & Bartlett Learning, LLC, an Ascend Learning Company. NOT FOR SALE OR DISTRIBUTION
stefanel/shutterstock, Inc. CHAPTER 2 First-Order Differential Equations We begin our stu of differential equations with first-order equations. In this chapter we illustrate the three different was differential
More informationProperties of Limits
33460_003qd //04 :3 PM Page 59 SECTION 3 Evaluating Limits Analticall 59 Section 3 Evaluating Limits Analticall Evaluate a it using properties of its Develop and use a strateg for finding its Evaluate
More informationBroken symmetry in modified Lorenz model
Int. J. Dnamical Sstems and Differential Equations, Vol., No., 1 Broken smmetr in modified Lorenz model Ilham Djellit*, Brahim Kilani Department of Mathematics, Universit Badji Mokhtar, Laborator of Mathematics,
More information17.3. Parametric Curves. Introduction. Prerequisites. Learning Outcomes
Parametric Curves 17.3 Introduction In this section we eamine et another wa of defining curves - the parametric description. We shall see that this is, in some was, far more useful than either the Cartesian
More informationMAT 1275: Introduction to Mathematical Analysis. Graphs and Simplest Equations for Basic Trigonometric Functions. y=sin( x) Function
MAT 275: Introduction to Mathematical Analsis Dr. A. Rozenblum Graphs and Simplest Equations for Basic Trigonometric Functions We consider here three basic functions: sine, cosine and tangent. For them,
More informationCHAPTER 2: Partial Derivatives. 2.2 Increments and Differential
CHAPTER : Partial Derivatives.1 Definition of a Partial Derivative. Increments and Differential.3 Chain Rules.4 Local Etrema.5 Absolute Etrema 1 Chapter : Partial Derivatives.1 Definition of a Partial
More informationReview of Prerequisite Skills, p. 350 C( 2, 0, 1) B( 3, 2, 0) y A(0, 1, 0) D(0, 2, 3) j! k! 2k! Section 7.1, pp
. 5. a. a a b a a b. Case If and are collinear, then b is also collinear with both and. But is perpendicular to and c c c b 9 b c, so a a b b is perpendicular to. Case If b and c b c are not collinear,
More informationRegular Physics - Notes Ch. 1
Regular Phsics - Notes Ch. 1 What is Phsics? the stud of matter and energ and their relationships; the stud of the basic phsical laws of nature which are often stated in simple mathematical equations.
More informationDynamical aspects of multi-round horseshoe-shaped homoclinic orbits in the RTBP
Dnamical aspects of multi-round horseshoe-shaped homoclinic orbits in the RTBP E. Barrabés, J.M. Mondelo, M. Ollé December, 28 Abstract We consider the planar Restricted Three-Bod problem and the collinear
More informationApproximate step response of a nonlinear hydraulic mount using a simplified linear model
Journal of Sound and Vibration 99 (007) 656 663 Short Communication JOURNAL OF SOUND AND VIBRATION Approximate step response of a nonlinear hydraulic mount using a simplified linear model Song He, Rajendra
More informationClearly the passage of an eigenvalue through to the positive real half plane leads to a qualitative change in the phase portrait, i.e.
Bifurcations We have already seen how the loss of stiffness in a linear oscillator leads to instability. In a practical situation the stiffness may not degrade in a linear fashion, and instability may
More information7.7 LOTKA-VOLTERRA M ODELS
77 LOTKA-VOLTERRA M ODELS sstems, such as the undamped pendulum, enabled us to draw the phase plane for this sstem and view it globall In that case we were able to understand the motions for all initial
More informationOne Dimensional Dynamical Systems
16 CHAPTER 2 One Dimensional Dynamical Systems We begin by analyzing some dynamical systems with one-dimensional phase spaces, and in particular their bifurcations. All equations in this Chapter are scalar
More informationc 1999 Society for Industrial and Applied Mathematics
SIAM J. SCI. COMPUT. Vol., No. 6, pp. 978 994 c 999 Societ for Industrial and Applied Mathematics A STUDY OF MONITOR FUNCTIONS FOR TWO-DIMENSIONAL ADAPTIVE MESH GENERATION WEIMING CAO, WEIZHANG HUANG,
More information4 Strain true strain engineering strain plane strain strain transformation formulae
4 Strain The concept of strain is introduced in this Chapter. The approimation to the true strain of the engineering strain is discussed. The practical case of two dimensional plane strain is discussed,
More informationFragment Orbitals for Transition Metal Complexes
Fragment Orbitals for Transition Metal Complees Introduction so far we have concentrated on the MO diagrams for main group elements, and the first 4 lectures will link into our Main Group chemistr course
More informationWavelet-Based Linearization Method in Nonlinear Systems
Wavelet-Based Linearization Method in Nonlinear Sstems Xiaomeng Ma A Thesis In The Department of Mechanical and Industrial Engineering Presented in Partial Fulfillment of the Requirements for the Degree
More information