arxiv: v1 [math.fa] 3 Jan 2019

Transcription

1 DAMPED AND DIVERGENCE EXACT SOLUTIONS FOR THE DUFFING EQUATION USING LEAF FUNCTIONS AND HYPERBOLIC LEAF FUNCTIONS A PREPRINT arxiv:9.66v [mah.fa] Jan 9 Kazunori Shinohara Deparmen of Mechanical Sysems Engineering Daido Universiy - Takiharu-cho, Minami-ku, Nagoya , Japan shinohara@6.alumni.u-okyo.ac.jp January 4, 9 ABSTRACT According o he wave power rule, he second derivaive of a funcion () wih respec o he variable is equal o negaive n imes he funcion () raised o he power of n. Solving he ordinary differenial equaions numerically resuls in waves appearing in he figures. The ordinary differenial equaion is very simple; however, waves, including he regular ampliude and period, are drawn in he figure. In his sudy, he funcion for obaining he wave is called he leaf funcion. Based on he leaf funcion, he eac soluions for he undamped and unforced Duffing equaions are presened. In he ordinary differenial equaion, in he posiive region of he variable (), he second derivaive d ()/d becomes negaive. Therefore, in he case ha he curves vary wih he ime under he condiion () >, he gradien d()/d consanly decreases as ime increases. Tha is, he angenial vecor on he curve of he graph (wih he abscissa and he ordinae ()) changes from he upper righ direcion o he lower righ direcion as ime increases. On he oher hand, in he negaive region of he variable (), he second derivaive d ()/d becomes posiive. The gradien d()/d consanly increases as ime decreases. Tha is, he angen vecor on he curve changes from he lower righ direcion o he upper righ direcion as ime increases. Since he behavior occurring in he posiive region of he variable () and he behavior occurring in he negaive region of he variable () alernaely occur in regular inervals, waves appear by hese ineracions. In his paper, I presen seven ypes of damped and divergence eac soluions by combining rigonomeric funcions, hyperbolic funcions, hyperbolic leaf funcions, leaf funcions, and eponenial funcions. In each ype, I show he derivaion mehod and numerical eamples, as well as describe he feaures of he waveform. Keywords Leaf funcions Hyperbolic leaf funcions Lemniscae funcions Jacobi ellipic funcions Ordinary differenial equaions Duffing equaion Nonlinear equaions Inroducion. Hyperbolic Leaf Funcions In our previous repor, he eac soluions for he Duffing equaion under free vibraion were presened using leaf funcions []. Leaf funcions comprise wo ypes of funcions wih iniial condiions. One of hem saisfies he following ordinary differenial equaion and iniial condiions. d d = n()n () - Takiharu-cho, Minami-ku, Nagoya , Japan

2 A PREPRINT - JANUARY 4, 9 () = () d() = () d The variable n represens an ineger. In he paper, he variable was named as he basis. A funcion ha saisfies he above equaions is defined as follows. () = sleaf n () (4) Anoher funcion saisfies he following ordinary differenial equaions and iniial condiions. d d = n()n (5) () = (6) d() = (7) d The variable n represens an ineger. A funcion ha saisfies he above equaions is defined as follows. () = cleaf n () (8) For he basis n =, he funcion sleaf () and cleaf () represen he rigonomeric funcions sin() and cos(), respecively. For he basis n =, he funcions sleaf () and cleaf () represen he lemniscae funcions sl() and cl(), respecively. Based on resuls obained from solving he ordinary differenial equaion numerically, Eqs. ()-() or Eqs. (5)-(7) yield waves wih respec o an arbirary n. The ordinary differenial equaion is very simple. However, waves including he regular ampliude and period are drawn in he figure. On he righ side of Eq. () and Eq. (5), if he minus changes o plus, he waves disappear wih respec o any arbirary n. The variable () increases monoonically as he variable increases. These ordinary differenial equaions are as follows: d d = n()n (9) () = () d() = () d In his sudy, a funcion ha saisfies he above equaions is defined as sleafh n () []. The basis n represens an ineger. () = sleafh n () () For he basis n =, he funcion sleafh () represens he hyperbolic funcion sinh(). In case of he basis n =, he funcion sleafh () represens he Hyperbolic Lemniscae Funcion slh() or sinhlemn() [, 4]. We furher discuss combinaions of he following differenial equaions and iniial condiions. d d = n()n () () = (4) d() = (5) d The basis n represens an ineger. A funcion ha saisfies he above equaions is defined as follows [5]. () = cleafh n () (6) For he basis n =, he funcion cleafh () represens he hyperbolic funcion cosh(). For he basis n =, he leaf funcion ha saisfies Eq. () - (5) is cleafh (). In he lieraure, he corresponding funcions wih respec o he funcion cleafh () canno be found [4]. For n =, Ramanujan suggesed an inverse funcion using a series []. Apar from his ype of soluion, eac soluions have no been presened. By applying an imaginary number o he phase of a rigonomeric funcion, he relaion beween he rigonomeric funcion and he hyperbolic funcion can be obained. Similar analogies eis beween leaf funcions and hyperbolic leaf funcions. Using he imaginary number for he phase of his leaf funcion, he relaion beween he leaf funcion and he hyperbolic leaf funcion can be derived []. Firs, we describe he relaion beween he leaf funcion sleaf n () and he hyperbolic leaf funcion sleafh n (). When he basis n(= m ) is an odd number, he following relaion is obained. sleaf m (i ) = i sleafh m ()(m =,,, ) (7)

3 A PREPRINT - JANUARY 4, 9 If he basis n(= m) is an even number, he following relaion is obained. sleaf m (i ) = i sleaf m ()(m =,,, ) (8) sleafh m (i ) = i sleafh m ()(m =,,, ) (9) Wih respec o an arbirary basis n, he relaion beween he leaf funcion cleaf n () and he hyperbolic funcion cleafh n () is given as follows: cleafh n (i ) = cleaf n ()(m =,,, ) () Using he above relaions, we can derive he relaion equaion beween he hyperbolic leaf funcion sleafh n () and cleafh n () using he relaion equaion beween he leaf funcions sleaf n () and cleaf n () a he basis n =, and. We can also derive he addiion heorem for he hyperbolic leaf funcions sleafh n () and cleafh n () using he addiion heorem beween sleaf n () and cleaf n () [, 6]. These relaion equaions and addiion heorems are boh heoreically and numerically consisen.. Comparison of Legacy funcions wih boh leaf funcions and hyperbolic leaf funcions The leaf funcions and hyperbolic leaf funcions based on he basis n = are as follows: sleaf () = sin() () cleaf () = cos() () sleafh () = sinh() () cleafh () = cosh() (4) The lemniscae funcions are presened by Gauss [7, 8]. The relaion equaions beween his funcion and he leaf funcion are as follows: sleaf () = sl() (5) cleaf () = cl() (6) sleafh () = slh() (7) The definiion of he funcion slh() of Eq.(7) can be confirmed by references [9, ]. For he funcion corresponding o he hyperbolic leaf funcion cleafh (), no clear descripion can be found in he lieraure. In he case ha he basis of a leaf funcion or a hyperbolic leaf funcion becomes n = or higher, i canno be represened by he lemniscae funcion. We can obain he second derivaives of he Jacobi ellipic funcion, he lemniscae funcion, he leaf funcion, or he hyperbolic leaf funcion. These funcions are defined as (). These second derivaives of he funcion () can be described by he sum of erms obained by raising he original funcion, (). d () d = c + c () + c () + c () (c c : consans) (8) By deciding he values of he coefficiens c c, one of eiher funcions in he lemniscae funcion, he leaf funcion and he hyperbolic funcion is deermined. The maimum value of he eponen in he above equaion is. Hisorically, discussions have no been had in he lieraure wih respec o eponens above 4. Using definiions based on he leaf funcion or he hyperbolic leaf funcion, we can discuss differenial equaions composed of higher order eponens such as he following d () d = c + c () + c () + c () + c 4 () 4 + c 5 () 5 (9) + + c n () n + c n () n (c c n : consans). The Duffing equaion For he basis n =, boh leaf funcions and hyperbolic funcions are applied o he damped or divergence Duffing equaion. The eac soluions consis of seven ypes of damped and divergence eac soluions. In each case, he derivaion mehod and numerical eamples are shown and he feaures of he waveform described. The Duffing equaion consiss of he second derivaive of he unknown funcion, he firs derivaive of he unknown funcion, he unknown funcion, he cube of he unknown funcion and a rigonomeric funcion. The equaion is presened as follows: d () d + δ d() + α() + β() = F cos(ω) () d The variable () in he above equaion represens an unknown funcion. The unknown funcion depends on ime,. The sign d()/d and he sign d ()/d represen he firs and second derivaives of (),respecively. The coefficiens δ, α, β, F and ω do no depend on ime,. If he above equaion is regarded as a mahemaical model of mechanical vibraion, he firs, second, hird and fourh erms on he lef side of Eq. () represen he inerial, damping, rigidiy, and nonlinear siffness erms, respecively. The erm on he righ side of Eq. () represens he eernal force erm.

4 A PREPRINT - JANUARY 4, 9.4 Solving he Duffing Equaion The mehod for finding he soluion of he damped Duffing equaion is roughly divided ino wo ypes: numerical soluions and eac soluions. In numerical soluions, a simple echnique based on Taylor epansion is applied o deermine an approimae soluion for a nonlinear Duffing oscillaor wih damping effec under differen iniial condiions []. The nonlinear problem of he Duffing equaion has been considered using he homoopy analysis echnique. Compared o he perurbaion mehods [], homoopy reamen does no require any small parameers [].The dynamic behaviors and effecs of random parameers are invesigaed. The Chebyshev orhogonal polynomial approimaion mehod is applied o reduce he random parameer [4]. The harmonic oscillaions of freedom of a Duffing oscillaor wih large damping are invesigaed using a simple poin collocaion mehod [5].The improved consrained opimizaion harmonic balance mehod is employed o solve he Duffing oscillaor. Analyical gradiens of he objec funcion are formulaed. The sensiiviy informaion of he Fourier coefficiens can also be obained [6]. Conversely, eac soluions of he Duffing equaion using he Jacobi ellipic funcion are discussed. The eac soluion for he cubic-quinic Duffing oscillaor based on he use of Jacobi ellipic funcions is presened [7]. Furhermore, he eac soluions for he damped Duffing equaion are presened by eending he parameers of he eac soluion [8]. The auhors in Ref. [9] do no assume any epression for he soluion bu eacly solved he nonlinear differenial equaion, unlike he Elias-Zuniga procedure. An analyical soluion for he damped Duffing equaion was derived by seing he ime-dependen modulus in he Jacobi ellipic funcion []..5 Originaliy and Purpose In our previous paper, eac soluions for undamped and unforced Duffing equaions were presened by Ref. [] using he leaf funcion. Conversely, he divergence phenomena (or damped vibraion) can be obained in he Duffing equaion. The purpose of his paper is o presen more eac soluions for he Duffing equaion by combining hyperbolic leaf funcions, leaf funcions, and eponenial funcions. These combining funcions make i possible o produce divergence phenomena or damped vibraion for he Duffing equaion. To represen he eac soluion ha decay wih ime, he eponenial funcion is placed in he phase of he leaf funcion. In his ype of soluion, he descripion of he eac soluions could no be found in lieraure. In conras, in he eac soluions ha diverge wih ime, he original funcions in he eac soluions also need o diverge wih ime. However, here are no eamples defined as soluions for he Duffing equaion using he lemniscae funcion slh(). In his paper, o represen he eac soluion for he Duffing equaion, he inegral funcion of he leaf funcion or he hyperbolic leaf funcion is placed in he phase of he rigonomeric funcions or hyperbolic funcions. The descripion of he eac soluions could also no be found in he lieraure. In his paper, seven ypes of he eac soluion are presened, ou of which five are divergence soluions and wo are damped soluions wihou eernal forces. In an eac soluion using he Jacobi ellipic funcion, muliple parameers of he Jacobi ellipic funcion cumulaively influence he periods and ampliudes. Therefore, he parameer based on a period is hard o separae from he parameer based on an ampliude. To deermine hese parameers, we need o solve a sih order equaion [8]. The eac soluions presened in his paper have he feaure ha he coefficiens of he Duffing equaion can be represened in boh ampliude and phase by wo parameers using a leaf funcion, a hyperbolic leaf funcion, and an eponenial funcion. Numerical daa for he Hyperbolic leaf funcion In he hyperbolic leaf funcion described in he previous secion, curves and numerical daa were described. The eac soluions are applied o he hyperbolic leaf funcions or leaf funcions based on he basis, n =. The graph of he hyperbolic leaf funcion based on he basis n = is shown in Fig.. The verical and horizonal aes represen he variables () and, respecively. For he basis n =, he numerical daa are summarized in Table. The following equaions are obained from he iniial condiions in Eq. () and Eq. (4). () = sleafh () = () () = cleafh () = () The hyperbolic leaf funcion sleafh () is an odd funcion whereas he hyperbolic leaf funcion cleafh () is an even funcion. The following relaions are obained. sleafh ( ) = sleafh () () cleafh ( ) = cleafh () (4) The hyperbolic leaf funcions sleafh () and cleafh () have limis. Given ha each limi is ζ and η, hey are obained by he following equaions [, 5]. ζ = du =.8547 (5) + u 4 4

5 A PREPRINT - JANUARY 4, 9 Table : Numerical daa for sleafh (), cleafh (), sleafh (u)du and cleafh (u)du sleafh () cleafh () sleafh (u)du cleafh (u)du (Noe) Numerical daa under he inequaliy < are obained using Eqs. () and (4). η = du =. (6) u4 From hese ereme values, he following relaions are obained by he limis ζ and η. lim sleafh () = + (7) ζ lim sleafh () = (8) ζ lim cleafh () = + (9) ±η Eac soluion for he Duffing equaion We consider he case where δ = and F = in he Duffing equaion of (). Three ypes of divergence eac soluions were presened in he above equaion. We define hose ordinary differenial equaions and iniial condiions as ype (VII I) o ype (XIV). (See Ref. [] for ypes (I) o (VII)). The inegral funcion of he leaf funcions and he hyperbolic leaf funcions are defined as follows: SL () = SLH () = CLH () = sleaf (u)du (4) sleafh (u)du (4) cleafh (u)du (4) The variables (), A,, ω, u, φ and B represen he displacemen, ampliude, ime, angular frequency, parameer, iniial phase and consan, respecively. 5

6 A PREPRINT - JANUARY 4, sleafh u du sleafh. cleafh u du cleafh Figure : Curves of hyperbolic leaf funcions and heir inegral funcions. Type (VIII) (See Appendi (VIII )in deail) Eac Soluion: Ordinary Differenial Equaion: Iniial Condiion: Iniial Velociy: d() d () = A cosh(clh (ω + φ)) (4) d () d + ω () 4( ω A ) () = (44) () = A cosh(clh (φ)) (45) = A sinh(clh (φ)) ω cleafh (φ) (46). Type (IX) (See Appendi (IX) in deail) Eac Soluion: Ordinary Differenial Equaion: Iniial Condiion: Iniial Velociy: d() d () = A sinh(clh (ω + φ)) (47) d () d ω () 4( ω A ) () = (48) () = A sinh(clh (φ)) (49) = A cosh(clh (φ)) ω cleafh (φ) (5) 6

7 A PREPRINT - JANUARY 4, 9. Type (X) (See Appendi (X) in deail) Eac Soluion: () = A cos(sl (ω + φ)) A sin(sl (ω + φ)) + A cosh(clh (ω + φ)) (5) Ordinary Differenial Equaion: Iniial Condiion: Iniial Velociy: d () d + ω ( + )() ( ω A ) () = (5) () = { A cos(sl (φ) + π } 4 ) + cosh(clh (φ)) d() = A sin(sl (φ)) ω sleaf (φ) A cos(sl (φ)) ω sleaf (φ) d + A sinh(clh (φ)) ω cleafh (φ) (5) (54).4 Type (XI) (See Appendi (XI) in deail) Ne, consider he case where F = in he Duffing equaion(eq. ()). In he eac soluion ha saisfies he Duffing equaion, negaive damping (divergence) based on ype (XI) and ype (XII) are presened. The ordinary differenial equaions and iniial condiions are as follows. Eac Soluion: Ordinary Differenial Equaion: () = A e ω sleafh (B e ω + φ) (55) Iniial Condiion: Iniial Velociy: d() d d () d ω d() d + ω () ( Bω A ) () = (56) () = A sleafh (B + φ) (57) = Aω sleafh (B + φ) + ABω + (sleafh (B + φ)) 4 (58).5 Type (XII) (See Appendi (XII) in deail) Eac Soluion: Ordinary Differenial Equaion: Iniial Condiion: Iniial Velociy: d() d d () d () = A e ω cleafh (B e ω + φ) (59) ω d() d.6 Type (XIII) (See Appendi (XIII) in deail) + ω () ( Bω A ) () = (6) () = A cleafh (B + φ) (6) = Aω cleafh (B + φ) + ABω (cleafh (B + φ)) 4 (6) In he eac soluion ha saisfies Eq. () under he condiion F =, damping vibraions based on ype (XIII) and ype (XIV) are presened. The ordinary differenial equaions and iniial condiions are as follows. Eac Soluion: () = Ae ω sleaf (B e ω + φ) (6) 7

8 A PREPRINT - JANUARY 4, 9 Ordinary Differenial Equaion: Iniial Condiion: Iniial Velociy: d() d d () d ω d() d + ω () + ( Bω A ) () = (64) () = A sleaf (B + φ) (65) = Aω sleaf (B + φ) + ABω (sleaf (B + φ)) 4 (66).7 Type (XIV) (See Appendi (XIV) in deail) Eac Soluion: Ordinary Differenial Equaion: Iniial Condiion: Iniial Velociy: 4 Analysis of resuls d() d d () d 4. Divergence soluion of ype(viii) () = Ae ω cleaf (B e ω + φ) (67) ω d() d + ω () + ( Bω A ) () = (68) () = A cleaf (B + φ) (69) = Aω cleaf (B + φ) ABω (cleaf (B + φ)) 4 (7) Le us consider a case where A =, ω = and φ = in Eq. (4), he inegral funcion, he hyperbolic funcion cosh() and he hyperbolic leaf funcion cleafh () are shown in Fig.. The verical and horizonal aes represen he variables () and, respecively. As shown in Eq. (6) and (9), he hyperbolic leaf funcion cleafh () has limis [5]. Therefore, he hyperbolic leaf funcion cleafh () increases sharply near he limi = ±η. The inegral of he hyperbolic leaf funcion cleafh () is obained as follows [5]: CLH () = cleafh (u)du = ln( (cleafh ()) + + (cleafh ()) ) ln( ) (7) The funcion cleafh () is included in he inegral funcion, so ha he inegral funcion also increases sharply near he limi = ±η. In he eac soluion of ype (VIII) including he inegral funcion l cleafh ()d, he (VIII) eac soluion increases monoonically near he limi = ±η. In he case ha he curves vary wih he ampliude A under he condiions ω = and φ =, he curves obained by he (VIII) soluion are shown in Fig.. The limis = ±η do no vary even if he parameer A varies. For he inequaliy A >, he (VIII) soluion increases monoonically wih ime,. For he inequaliy A <, he curve of he (VIII) soluion decreases monoonically wih. The iniial posiion () a = (Eq. (45)) also varies according o he parameer A. In case ha he curves vary wih he phase ω under he condiions A = and φ =, he curves obained by he (VIII) soluion are shown in Fig.4. The limi is = ± η ω. The limi varies according o he phase ω. As he phase ω increases, he absolue value of he limi becomes smaller and he possible domain of he variable becomes narrower. As he phase ω decreases, he absolue value of he limi becomes larger and he possible domain of he variable become wider. 4. Divergence soluions of he ype(ix) Le us consider a case where A =, ω = and φ = in Eq. (47), he inegral funcion, he hyperbolic funcion sinh(), and he hyperbolic leaf funcion cleafh () are shown in Fig. 5. The verical and horizonal aes represen he variables () and, respecively. The hyperbolic leaf funcion cleafh () has limis [5]. Therefore, he eac soluion 8

9 A PREPRINT - JANUARY 4, 9 cosh cosh cleafh u du cleafh u du () cleafh Figure : Curves obained by he (VIII) soluion, he hyperbolic funcion cosh(), he hyperbolic leaf funcion cleafh () and he funcion cleafh ()d of ype (IX) increases monoonically near he limi = η, and decreases monoonically near he limi = η. In case ha curves vary wih he ampliude A under he condiions ω = and φ =, he curves obained by he (IX) soluion are shown as in Fig.6. The limis of = ±η do no vary even if he parameer A varies. For he inequaliy A >, he (IX) soluion increases monoonically wih ime,. For he inequaliy A <, he (IX) soluion decreases monoonically wih. The curve passes hrough () =. In case ha he curves vary wih he ampliude ω under he condiions A = and φ =, he curves obained by he (IX) soluion are shown as in Fig.7. The limi is = ± η ω. The limi varies according o he phase ω. As he phase ω increases, he absolue value of he limi decreases and he possible domain of he variable becomes narrower. As he phase ω decreases, he absolue value of he limi increases and he possible domain of he variable becomes wider. 4. Divergence soluion of he ype (X) Le us consider a case where A =, ω = and φ = in Eq. (5). The (X) eac soluion, he hyperbolic funcion cosh() and he hyperbolic leaf funcion cleafh () are shown in Fig. 8. The verical and horizonal aes represen he variables () and, respecively. As shown in Fig. 8, hese funcions are even funcions ha are symmerical abou he () ais. The limis eis in he (X) eac soluion. The (X) eac soluion is ransformed as follows: () = cos(sl ()) sin(sl ()) (7) Using he Eq. (III.4) in Ref. [], squaring boh sides of he above equaion yields he following equaion. ( ()) = (cos(sl ()) + (sin(sl ()) sin(sl ()) cos(sl ()) = sin( SL ()) = (sleaf ()) (7) The variable () is obained as follows: Ne, he following equaion is ransformed: () = ± (sleaf ()) (74) () = cosh(clh ()) (75) Using he Eq. (H5) in Ref. [5], squaring boh sides of he above equaion gives he following equaion. ( ()) = (cosh(clh ())) = + cosh(clh ()) = + (cleafh ()) (76) 9

10 A PREPRINT - JANUARY 4, 9 () A A 5 4 A Acosh cleafh u du A 4 5 A A Figure : Curves obained by he (VIII) eac soluion as a funcion of variaions in he ampliude A (A = ±, ±, ±) (Se ω = and φ = ) As shown in Fig. 8, he inequaliy () > is obvious. The following equaion is obained. Therefore, he (X) eac soluion is obained as follows: () = + (cleafh ()) (77) () = () + () = ± (sleaf ()) + + (cleafh ()) (78) Because he (X) eac soluion conains he funcion cleafh (), is eac soluion also has limis according o he limis of cleafh (). The domain of (X) is as follows. η < < η (79) The sign η is a consan described by Eq. (6). The sign η has he following relaion given by Eq. (4) in Ref. [5]. η = π ( 4 d = d) (8) 4 In he case ha he curves vary wih he ampliude A under he condiions ω = and φ =, he curves obained by he (X) soluion are shown in Fig. 9. The limi = ±η does no vary even if he parameer A varies. Conversely, he iniial condiion (Eq. (5)) varies wih he ampliude, A. In case ha he curves vary wih he ampliude ω under he condiions A = and φ =, he curves obained by he (X) soluion are shown in Fig.. The limi is = ± η ω. The limi varies wih he phase, ω. As he phase ω increases, he absolue value of he limi decreases and he possible domain of he variable becomes narrower. As he phase ω decreases, he absolue value of he limi increases and he possible domain of he variable becomes larger. 4.4 Divergence soluion of he ype (XI) Le us consider a case where A =, ω =, φ = and B = in Eq. (55), he (XI) eac soluion, he hyperbolic funcion sinh() and he hyperbolic leaf funcion sleafh () are shown in Fig.. Based on Eq. (57), he curve passes hrough () =. For he inequaliy <, he variable () converges asympoically o zero. For >, he variable () increases monoonically. The limi eiss for he hyperbolic leaf funcion sleafh (e ) (See Ref. []). The possible range of e is as follows: ζ < e < ζ (8)

11 A PREPRINT - JANUARY 4, 9 cosh cleafh u du () Figure 4: Curves obained by he (VIII) eac soluion as a funcion of variaions in he angular frequency ω (ω = ±, ±, ±) (Se A = and φ = ) The following equaion is obained from he above equaion. ζ + < e < ζ + (8) The consan ζ is obained by he Eq. (5) (See Ref. []). In Eq. (8), he inequaliy ζ + =.8547 < e always holds for any arbirary ime,. We can ake he logarihm of boh sides in he inequaliy e < ζ + in Eq. (8). Solving for he variable yields he following equaion. < ln(ζ + ) (= ln( ) =.487 ) (8) The above inequaliy represens he domain of he variable in Eq. (8). In case ha curves vary wih he ampliude A under he condiions ω =, φ = and B =, he curves obained by he (XI) soluion are shown in Fig.. The limis do no vary even if he parameer A varies. = ln(ζ + ) = ln( ) =.487 (84) For he inequaliy A >, he eac soluion () of (XI) diverges o plus infiniy. On he oher hand, for he inequaliy A <, he eac soluion () of (XI) diverges o minus infiniy. In case ha he curves vary wih he phase ω under he condiions A =, φ = and B =, he curves obained by he (XI) soluion are shown in Fig.. Le us consider he limi of he (XI) eac soluion when he parameer ω varies. The limi of he variable is obained by he following equaion. e ω = ±ζ = ±.8547 (85) When he righ side of Eq. (85) is.8547, solving he above equaion for he variable yields he following equaion. Le us consider a case where he righ side of Eq. (85) is = ω ln(ζ + ) (86) e ω =.8547 (87) The variable ha saisfies he above equaion does no eis. Therefore, he condiion for he eisence of he limi is ζ = When ω > in Eq. (86), he limi of he variable is generaed in he domain where >. When ω < in Eq. (86), he limi of he variable is produced in he domain <. In case ha he curves vary wih he

12 A PREPRINT - JANUARY 4, 9 Table : Numerical daa for he (VIII) soluion under he condiions A =, φ = and ω = () () d d phase A under he condiions ω =, φ = and B =, he curves obained by he (XI) soluion are shown in Fig.4. The possible range of e in he above equaion is obained as follows: ζ < e < ζ (88) The above equaion is ransformed as an inequaliy of he variable. The following equaion is obained. ln( ζ ) < < ln( + ζ ) (89) Because he consan ζ is defined by he Eq. (5), he following equaion is obained by a numerical value..67 < <.655 (9) Therefore, one or wo limis depend on he variables B and φ. For A, he eac soluion () of he ype (XI) passes hrough zero and increases monoonously. For A, he eac soluion () of he ype (XI) passes hrough zero and decreases monoonously. The limis =.67 and =.6559 do no vary even if he parameer A varies. In case ha he curves vary wih he ampliude ω under he condiions A =, φ = and B =, he curves obained by he (XI) soluion are shown in Fig.5. The possible range of he variable in he above equaion is as follows: ω ln( ζ ) < < ω ln( + ζ ) (9) As he absolue value of ω increases in he above equaion, he limi of he variable approaches zero, and he domain ha he variable can ake is narrowed. As he absolue value of ω decreases, he limi of he variable goes away from. Therefore, he domain of he variable widens. 4.5 Negaive Damping Soluion of he ype (XII) Using he hyperbolic leaf funcion, he negaive damping soluion of he ype (XII) is presened in his secion. In he (X II) eac soluion, he damped erm for he Duffing equaion d()/d does no work in he direcion of he suppressing

13 A PREPRINT - JANUARY 4, 9 cleaf sinh u du sinh cleafh u du () cleafh Figure 5: Curves obained by he (IX) eac soluion, he hyperbolic funcion sinh(), he leaf funcion cleafh () and he funcion cleafh (u)du vibraion. Le us consider when A =, ω =, φ = and B = in Eq. (59), he eac soluion of he ype (XII) is shown in Fig. 6. For he inequaliy <, he variable () converges asympoically o zero. The curve passes hrough () =. In he domain of he inequaliy >, he variable () diverges o infiniy. The hyperbolic leaf funcion cleafh (e ) has a limi. Therefore, he range ha he variable e can ake is as follows. The following equaion is obained from he above equaion. η < e < η (9) η + < e < η + (9) The consans η is obained by he Eq. (6) (See Ref. [5]). In Eq. (9), he inequaliy η + =. < e always holds for arbirary ime. In he inequaliy e < η + in Eq. (9), we can ake he logarihm of boh sides. Solving for he variable yields he following equaion. < ln(η + ) = ln(. ) =.8769 (94) The above inequaliy represens he domain of he variable. In case ha he curves vary wih he ampliude A under he condiions ω =, φ = and B =, he curves obained by he (XII) soluion are shown in Fig. 7. The limis do no vary even if he parameer A varies. For he inequaliy A >, he eac soluion () for (XII) diverges o posiive infiniy. In case ha he curves vary wih he phase ω under he condiions A =, φ = and B =, he curves obained by he (XII) soluion are shown in Fig.8. Le us consider he limi of he eac soluion for (XII) when he parameer ω varies. The limi value for he variable is obained using he following equaion. e ω = ±η = ±. (95) When he righ side of Eq. (95) is., solving he above equaion for he variable yields he following equaion. = ω ln(η + ) =.6 ln(. + ) = (96) ω ω When ω > in Eq. (95), he limi of he variable eiss in he domain >. When ω <, he limi of he variable eiss in he domain <. As he absolue value of ω increases in he above equaion, he limi of he variable approaches zero. The domain of possible values of becomes narrow. As he absolue value of ω decreases, he limi of

14 A PREPRINT - JANUARY 4, 9 A A A () Asinh cleafh u du A A A Figure 6: Curves obained by he (IX) eac soluion as a funcion of variaions in he ampliude A (A = ±, ±, ±) (Se ω = and φ = ) he variable goes away from. Therefore, he domain of he variable widens. When he righ side of equaion (95) is -., solving he above equaion for he variable yields he following equaion. = ω ln( η + ) = ln(. ) ln(. + ) = (97) ω ω Because he logarihm ake he negaive value (. ) in parenhesis, he variable ha saisfies he above equaion does no eis. In case ha he curves vary wih he phase A under he condiions ω =, φ = and B =, he curves obained by he (XII) soluion are shown in Fig. 9. The possible range of e in he above equaion is as follows: η < e < η (98) The above equaion is ransformed as an inequaliy of he variable o yield he following equaion. ln( η ) < < ln( + η ) (99) Because he consan η is defined by Eq. (6), he following equaion is obained < <.549 () One or wo limis depend on he variables B and φ. The eac soluion () for (XII) passes hrough () = A. In a siuaion where A, he eac soluion diverges o plus infiniy near boh =.6569 and =.549. In a case where A, he eac soluion diverges o minus infiniy near boh =.6569 and =.549. These limis do no vary even if he parameer A varies. In case ha he curves vary wih he ampliude ω under he condiions A =, φ = and B =, he curves obained by he (XII) soluion are shown in Fig.. The possible range of he variable in he above equaion is as follows: ω ln( η ) < < ω ln( + η ) () As he absolue value of ω increases in he above equaion, he domain of possible values for becomes narrower and narrower as shown in he Eq. (). As he absolue value of ω decreases, he limi of he variable goes away from. Therefore, he domain of he variable becomes wider and wider. 4.6 Damping Soluion of he ype (XIII) Le us consider Eq. (6) ha saisfies he Duffing equaion. Physically, he eac soluion of he ype (XIII) corresponds o he nonlinear spring model of he damping sysem. 4

15 A PREPRINT - JANUARY 4, () sinh cleafh u du Figure 7: Curves obained by he (IX) eac soluion as a funcion of variaions in he ampliude ω (ω = ±, ±, ±) (Se A = and φ = ) 4.6. Period When A =, ω =, φ = and B = π /, he curves obained by he (XIII) soluion are shown in Fig.. The horizonal ais represens ime whereas he verical ais represens displacemen, (). As shown in Fig., he period of he wave varies wih ime. In he Fig., he T, T,, T 6 represens he firs period, he second period,, and he sih period, respecively. Firs, le us consider he firs period T. Subsiuing - ino he variable, he value of he ype (XIII) can be obained as follows: ( ) = e sleaf ( π e ) = sleaf ( π ) = () On he oher hand, subsiuing ln4 ino he variable, he value of he ype (XIII) can be obained as follows: (ln4) = e ln4 sleaf ( π eln4 ) = 4 sleaf (π ) = () The firs period T is as follows: T = ln4 = (4) The domain of he firs period T is as follows: < ln4 (5) Ne, le us consider he second period T. Subsiuing ln8 ino he variable, he value of he ype (XIII) can be obained as follows: (ln8) = e ln8 sleaf ( π eln8 ) = 8 sleaf (4π ) = (6) The second period T is as follows: T = ln8 ln4 (7) The domain of he second period T is as follows: ln4 ln8 (8) In his way, he m h period T m can be generalized as follows: T m = ln(4m) ln(4m 4) = ln m (m =,,, ) (9) m The domain of he m h period T m is as follows: ln(4m 4) ln(4m)(m =,,, ) () 5

16 A PREPRINT - JANUARY 4, 9 Table : Numerical daa for he (IX) soluion under he condiions A =, φ = and ω = () () d d Ampliude A ime wih respec o he conve upward of wave is obained by using he gradien of he ype (XIII). Using he condiion d()/d = in he Eq. (8), we can obain he equaion as follows: B e ω (sleaf (B e ω )) 4 + (sleaf (B e ω )) B e ω = () Unforunaely, we can no obain he eac soluion of he variable by he above equaion. Using he numerical analysis, he variable ha saisfy wih he Eq. () is obained. Under he condiions B = π /, A =, ω = and φ =, he ampliude () wih respec o he ime are summarized in he Table 7 and Table 8. The curve () = ±e is added in Fig.. As shown in Fig., he curve e inersecs he curve sleaf ( π e ) a he condiion sleaf ( π e ) = ±. In he case where he condiion sleaf ( π e ) =, he ime is obained as follows: The following equaion is obained by he above equaion. π e = 4k π (k =,, ) () k = ln(4k )(k =,,, ) () Subsiuing k ino he Eq. (6) for he soluion of ype (XIII), he following equaion is obained. ( k ) = ep( k ) sleaf ( π ep( k)) = ep(ln(4k )) sleaf ( π ep(ln(4k ))) = (4k )sleaf ( π (4k )) = 4k (k =,,, ) (4) 6

17 A PREPRINT - JANUARY 4, 9 cleafh cos sleaf u du sin sleaf u du cosh cleafh u du 5 4 () u cosh cleafh du cleafh cosh cos sleaf u du sin sleaf u du sleaf Figure 8: Curves obained by he (X) eac soluion, he hyperbolic funcion cosh(), he leaf funcion cleafh (), he funcion + (cleafh ()) and he funcion (sleafh ()) In he case where he condiion sleaf ( π e ) =, he ime is obained as follows: The following equaion is obained by he above equaion. Subsiuing k ino he Eq. (6), he following equaion is obained. π e = 4k π (k =,, ) (5) k = ln(4k )(k =,,, ) (6) ( k ) = ep( k ) sleaf ( π ep( k)) = ep(ln(4k )) sleaf ( π ep(ln(4k ))) = (4k )sleaf ( π (4k )) = (4k )(k =,,, ) (7) These daa boh he ime (by he Eq. () and he Eq. (6) ) and he ampliude () (by he Eq. (4) and he Eq. (7)) are summarized in Table 9 and Table. Since he condiion sleaf ( π e ) = ± is no saisfied wih he Eq. (), he values obained by he Eq. () and he Eq. (6) are no he eremal value in he ype (XIII). However, comparing he numerical daa in he able 7(or he able 8) wih he numerical daa in he able 9 (or he able ), he approimae ampliude (he able 9 and he able ) are almos agreed wih he eac ampliude by numerical analysis (he able 7 and he able 8). As an approimae value of he ampliude, ampliude can be easily obained from he Eq. (4) and he Eq. (7). Fig. and Fig. show he variaion of waves as he parameer A vary. As he parameer A increases, he heigh of he wave increases wih keeping he wave period. Le us consider he parameer ω under he condiions B = π /, φ = and A =. The curves for he (XIII) soluion are shown in he Fig. 4. In he negaive domain of he variable, he waves have no period. Conversely, in he posiive domain of he variable, he waves show periodiciy. As he magniude of ω increases, a period becomes shorer. As he magniude of ω decreases, a period becomes longer. As shown in Fig 5, for ω <, he wave has no periodiciy in he posiive domain. The eac soluion for (XIII) decreases monoonically as ime, increases. In he negaive domain of he variable, he wave shows periodiciy. 7

18 A PREPRINT - JANUARY 4, () Acos sleaf u du Asin sleaf u du Acosh cleafh u du A A A A A A Figure 9: Curves obained by he (X) eac soluion as a funcion of variaions in he ampliude A (A = ±, ±, ±) (Se ω = and φ = ) 4.7 Damping Soluion of he ype (XIV) Le us consider Eq. (6) ha saisfies he Duffing equaion. Physically, he eac soluion of he ype (XIV) corresponds o he nonlinear spring model of he damping sysem Period When A =, ω =, φ = and B = π /, he curves obained by he (XIV) are shown in Fig.6. The horizonal ais represens ime whereas he verical ais represens displacemen, (). As shown in Fig. 6, he period of he wave varies wih ime. In he Fig.6, he T, T,, T 6 represens he firs period, he second period,, and he sih period, respecively. Firs, le us consider he firs period T. Subsiuing - ino he variable, he value of he ype (XIV) can be obained as follows: ( ) = e cleaf ( π e ) = cleaf ( π ) = (8) On he oher hand, subsiuing ln4 ino he variable, he value of he ype (XIV) can be obained as follows: (ln4) = e ln4 cleaf ( π eln4 ) = 4 cleaf (π ) = 4 (9) The firs period T is as follows: T = ln4 = () The domain of he firs period T is as follows: < ln4 () Ne, le us consider he second period T. Subsiuing ln8 ino he variable, he value of he ype (XIV) can be obained as follows: (ln8) = e ln8 cleaf ( π eln8 ) = 8 cleaf (4π ) = 8 () The second period T is as follows: T = ln8 ln4 () The domain of he second period T is as follows: ln4 ln8 (4) 8

19 A PREPRINT - JANUARY 4, cos sleaf u du sin sleaf u du cosh cleafh u du () Figure : Curves obained by he (X) eac soluion as a funcion of variaions in he parameer ω (ω = ±, ±, ±) (Se A = and φ = ) In his way, he m h period T m can be generalized as follows: T m = ln(4m) ln(4m 4) = ln m (m =,,, ) (5) m The domain of he m h period T m is as follows: 4.7. Ampliude ln(4m 4) ln(4m)(m =,,, ) (6) A ime wih respec o he conve upward of wave is obained by using he gradien of ype (XIV). Using he condiion d()/d = in he Eq. (9), we can obain he equaion as follows: B e ω (cleaf (B e ω )) 4 + (cleaf (B e ω )) B e ω = (7) We can no obain he eac soluion of he variable by he above equaion. Using he numerical analysis, he variable ha saisfy wih he Eq. (7) is obained. The value and he ampliude () are summarized in he Table and Table. The curve () = ±e is added in Fig. 6. As shown in Fig. 6, he curve e inersecs he curve cleaf ( π e ) a he condiion cleaf ( π e ) = ±. In he case where he condiion cleaf ( π e ) =, he ime is obained as follows: π e = kπ (k =,, ) (8) The following equaion is obained by he above equaion. k = ln(4k)(k =,,, ) (9) Subsiuing k ino he Eq. (67) for he soluion of ype (XIV), he following equaion is obained. ( k ) = ep( k ) cleaf ( π ep( k)) = ep(ln(4k)) cleaf ( π ep(ln(4k))) = 4k cleaf ( π 4k) = 4k(k =,,, ) 9 ()

20 A PREPRINT - JANUARY 4, 9 Table 4: Numerical daa for he (X) soluion under he condiions A =, φ = and ω = () () d d In he case where he condiion cleaf ( π e ) =, he ime is obained as follows: π e = (k )π (k =,, ) () The following equaion is obained by he above equaion. Subsiuing k ino he Eq. (67), he following equaion is obained. k = ln(4k )(k =,,, ) () ( k ) = ep( k ) cleaf ( π ep( k)) = ep(ln(4k )) cleaf ( π ep(ln(4k ))) = (4k )cleaf ( π (4k )) = (4k )(k =,,, ) These daa boh he ime (by he Eq. (9) and he Eq. () ) and he ampliude () (by he Eq. () and he Eq. ()) are summarized in Table 4 and Table 5. Since he condiion cleaf ( π e ) = ± is no saisfied wih he Eq. (7), he values obained by he Eq. (9) and he Eq. ( ) are no he eremal value in he ype (XIV). However, comparing he numerical daa in he able 4(or he able 5) wih he numerical daa in he able (or he able ), he approimae ampliude (he able 4 and he able 5) are almos agreed wih he eac ampliude by numerical analysis (he able and he able ). As an approimae value of he ampliude, ampliude can be easily obained from he Eq. () and he Eq. (). Fig. 7 and Fig. 8 show he variaion of waves as he parameer A vary. As he parameer A increases, he heigh of he wave increases wih keeping he wave period. Le us consider he parameer ω under he condiions B = π /, φ = and A =. The curves for he (XIV) soluion are shown in he Fig. 9. In he negaive domain of he variable, he waves have no period. Conversely, in he posiive domain of he variable, he ()

21 A PREPRINT - JANUARY 4, 9 Table 5: Numerical daa for he (XI) soluion under he condiions A =, B =, φ = and ω= () () d d Table 6: Numerical daa for he (XII) soluion under he condiions A =, B =, φ = and ω = () () d d

22 A PREPRINT - JANUARY 4, 9 Table 7: Eac ampliude () a conve upward for (XIII) Number in Fig. Time Eac ampliude () (U) (U) (U) (4U) (5U) (6U) Table 8: Eac ampliude () a conve downward for (XIII) Number in Fig. Time Eac ampliude () (D) (D) (D) (4D) (5D) (6D) Table 9: Approimae ampliude () a conve upward for (XIII) Number in Fig. Time by Eq. () Approimae ampliude () by Eq. (4) (U). (=ln()). (U) (=ln(5)) 5. (U) (=ln(9)) 9. (4U) (=ln()). (5U).844 (=ln(7)) 7. (6U) (=ln()). Table : Approimae ampliude () a conve downward for (XIII) Number in Fig. Time by Eq. (6) Approimae ampliude () by Eq. (7) (D) (=ln()) -. (D) (=ln(7)) -7. (D) (=ln()) -. (4D).785 (=ln(5)) -5. (5D) (=ln(9)) -9. (6D) (=ln()) -.

23 A PREPRINT - JANUARY 4, 9 Table : Numerical daa for he (XIII) soluion under he condiions A =, B =, φ = and ω = () () d d Table : Eac ampliude () a conve upward for (XIV) Number in Fig. 6 Time Eac ampliude () (U) (U) (U) (4U) (5U) (6U) Table : Eac ampliude () a conve downward for he (XIV) Number in Fig. 6 Time Eac ampliude () (D) (D) (D) (4D) (5D) (6D) Table 4: Approimae ampliude () a conve upward for (XIV) Number in Fig. 6 Time by Eq. (9) Approimae ampliude () for (XIV) by Eq. () (U) (=ln(4)) 4. (U) (=ln(8)) 8. (U) (=ln()). (4U) (=ln(6)) 6. (5U) (=ln()). (6U).7858 (=ln(6)) 4.

24 A PREPRINT - JANUARY 4, 9 e sleafh e () sinh() sleafh Figure : Curves obained by he (XI) eac soluion, he hyperbolic funcion sinh() and he leaf funcion sleafh () Table 5: Approimae ampliude () a conve downward for (XIV) Number in Fig. 6 Time by Eq. () Approimae ampliude () for (XIV) by Eq. () (D) (=ln()) -. (D) (=ln(6)) -6. (D).5859 (=ln()) -. (4D).6957 (=ln(4)) -4. (5D) (=ln(8)) -8. (6D).9445 (=ln()) -. waves show periodiciy. As he magniude of ω increases, a period becomes shorer. As he magniude of ω decreases, a period becomes longer. As shown in Fig, for ω <, he wave has no periodiciy in he posiive domain. The eac soluion for (XIV) decreases monoonically as ime, increases. In he negaive domain of he variable, he wave shows periodiciy. 5 Conclusion In our previous paper, eac soluions ha saisfy undamped and unforced Duffing equaions were presened using he leaf funcion. The waveform obained from he eac soluion ehibis periodiciy. The waveform has a differen waveform from he sine and cosine waveforms. In his paper, eac soluions for divergence and damped Duffing soluions were presened using leaf funcions, hyperbolic leaf funcions, and eponenial funcions. When he soluion diverges wih ime, he original funcion also needs o diverge wih ime. However, here is no eample in he lieraure showing he soluion of he Duffing equaion using he lemniscae funcion slh(). The funcion slh() diverges as ime, increases. To represen he eac soluion for divergence in he Duffing equaion, he inegral funcion of he hyperbolic leaf funcion is implemened in he phase of he rigonomeric funcion. The hyperbolic leaf funcion has a limi. Depending on he limi of he hyperbolic leaf funcion, he eac soluion for he divergen ype also has a limi. One or wo limis occur depending on he iniial condiions. In he eac soluion for he damped Duffing equaion, he eponenial funcion was implemened in he phase of he leaf funcion o describe he eac soluion. Over ime, one period of he wave changes wih ime, depending on he condiion of he parameer. Finally, he periodic sae of he wave ransiions o a sae where here is no periodiciy. 4

25 A PREPRINT - JANUARY 4, 9 Table 6: Numerical daa for he (XIII) soluion under he condiions A =, B =, φ = and ω = () () d d Table 7: Numerical daa for he (XIV) soluion under he condiions A =, B =, φ = and ω = () () d d

26 A PREPRINT - JANUARY 4, 9 () A= A= A= Ae sleafh e ln A=- A=- A=- Figure : Curves obained by he (XI) eac soluion as a funcion of variaions in he parameer A (A = ±, ±, ±) (Se B =, ω = and φ = ) Appendi (VIII) The eac soluion of he ype (VIII) saisfies he cubic Duffing equaion. The soluion diverges wih ime,. The firs derivaive of he Eq. (4) is obained as follows: d() d The second derivaive of he above equaion is obained as follows. = A sinh(clh (ω + φ)) ω cleafh (ω + φ) (4) d () d = A cosh(clh (ω + φ)) ω (cleafh (ω + φ)) +A sinh(clh (ω + φ)) ω (cleafh (ω + φ)) 4 The following equaion is derived from Eq. (6) in Ref. [5]. (5) (cleafh (ω + φ)) = cosh( CLH (ω + φ)) (6) Using he above equaion, he equaion can be ransformed as follows. (cleafh (ω + φ)) 4 = (cosh( CLH (ω + φ)) = sinh( CLH (ω + φ)) Subsiuing he above wo equaions ino equaion (5) yields he following equaion. d () d = Aω cosh(clh (ω + φ)) cosh( CLH (ω + φ)) + Aω sinh(clh (ω + φ)) sinh( CLH (ω + φ)) = Aω cosh(clh (ω + φ) + CLH (ω + φ)) = Aω cosh( CLH (ω + φ)) = Aω cosh(clh (ω + φ)) + 4( ω A ){ A cosh(clh (ω + φ)) } (7) (8) 6

27 A PREPRINT - JANUARY 4, 9 ω=-. ω=-. ω=-. () ω=. ω=. ω= ln.9 ln 5.. ln e sleafh e ln ln.. ln Figure : Curves obained by he (XI) eac soluion as a funcion of variaions in he parameer ω (ω = ±., ±., ±.)(Se B =, A = and φ = ) Subsiuing equaion (4) ino he above equaion leads o he following equaion. d () d = ω () + 4( ω A ) (()) (9) The coefficiens α and β in Eq. () are as follows. α = ω (4) β = 4( ω A ) (4) Appendi (IX) The eac soluion of he ype (IX) saisfies he cubic Duffing equaion. The soluion diverges wih ime,. The firs derivaive of he above equaion is obained as follows. d() = A cosh(clh (ω + φ)) ω cleafh (ω + φ) (4) d The second derivaive of he above equaion is obained as follows. d () d = A sinh(clh (ω + φ)) ω (cleafh (ω + φ)) +A cosh(clh (ω + φ)) ω (cleafh (ω + φ)) 4 Subsiuing (6) and (7) ino he above equaion gives he following equaion. d () d = Aω sinh(clh (ω + φ)) cosh( CLH (ω + φ)) + Aω cosh(clh (ω + φ)) sinh( CLH (ω + φ)) = Aω sinh(clh (ω + φ) + CLH (ω + φ)) = Aω sinh( CLH (ω + φ)) = Aω sinh(clh (ω + φ)) + 4( ω A ){ A sinh(clh (ω + φ)) } (4) (44) 7