arxiv: v2 [math.ds] 9 Jun 2016

Transcription

1 Stability Tests for Second Order Linear and Nonlinear Delayed Models arxiv: v [math.ds] 9 Jun 16 Leonid Berezansky Dept. of Math, Ben-Gurion University of the Negev, Beer-Sheva 8415, Israel Elena Braverman Dept. of Math & Stats, University of Calgary, 5 University Dr. NW, Calgary, AB, Canada TN1N4 and Lev Idels Dept. of Math, Vancouver Island University, 9 Fifth St. Nanaimo, BC, Canada V9S5J5 Abstract For the nonlinear second order Lienard-type equations with time-varying delays ẍt)+ f k t,xt),ẋg k t)))+ l s k t,xh k t))) =, global asymptotic stability conditions are obtained. The results are based on the new sufficient stability conditions for relevant linear equations and are applied to derive explicit stability conditions for the nonlinear Kaldor-Kalecki business cycle model. We also explore multistability of the sunflower non-autonomous equation and its modifications. AMS Subject Classification: 34K, 9D5, 34K45, 34K1, 34K5 Keywords: Second order delay differential equations; global asymptotic stability; boundedness of solutions; Lienard-type nonautonomous linear and nonlinear delay differential equations; sunflower equation; the Kaldor-Kalecki business cycle model; variable delays 1 Introduction The second order delay differential equation ẍ+ft,xt),ẋt τ))+gt,xt),xt τ)) = 1.1) has a more than 65-year history of study, and was used to examine aftereffects in mechanics, physics, biology, medicine and economics see, for example, [18]). Recently, these models have been used to mimic regenerative vibrations in a milling process, a balancing motion and chatter vibrations. For example, a one degree of freedom milling equation ẍt)+aẋt)+bxt) = α[xt) xt τt))] 1.) Corresponding author. maelena@math.ucalgary.ca. Fax 43) Phone 43)

2 was introduced in [36]. The milling model with several delays ẍt)+aẋt)+bxt)+ p α k [xt) xt τ k )] k = was recently studied, mostly numerically, in [13, 19, ]. The following milling models with variable parameters were derived and examined in [18, 8, 9, 36, 37, 38]: ẍt)+aẋt)+bt)xt) = ct)xt τt)), 1.3) ẍt)+aẋt)+bxt)+ p α k t)[xt) xt τ k t))] k =. In economics, the well-known Kaldor-Kalecki business cycle model expressed as the delayed system of two nonlinear equations [15], in some cases can be reduced to the second order equation see, for example, [3]) ẍt)+[α βp xt))]ẋt)+γ[pxt)) ηxt)]+δpxt τ)) =. 1.4) Here px) is a frequently used in mathematical economics sigmoid function [15], e.g. px) = A A 1+e bx, and all coefficients are nonnegative constants. Different techniques were applied to study second-order delay equations in [6, 7, 11, 16, 1, 3, 6] and [31] [35]. Characteristic quasipolynomials were broadly used for local stability analysis of autonomous models, see, for example, [18]). The fixed point technique for second order differential and functional equations was pioneered by T. A. Burton [8, 9]. In the paper [1] explicit and easilyverifiable tests were obtained for the autonomous model ẍt) = p 1 ẋt)+p ẋt τ)+q 1 xt)+q xt τ). 1.5) Theorem 1.1. [1] Assume that at least one of the following conditions holds: 1) p 1 p >,q 1 >,q > or ) p 1 >,p >,q 1 >,q <. Then equation 1.5) is unstable. Theorem 1.. [1] Assume p 1 = p =, q > and denote B = τ q 1, D = τ q. Equation 1.5) is asymptotically stable if and only if q 1 < and there exists k N such that kπ < B < k +1)π, D < min { k) π B,k +1) π +B }. Example 1.3. The second-order delay equation ẍt) = 49xt)+7xt 1) 1.6) is asymptotically stable by Theorem 1.. Based on the algorithmic tests presented in [1], the equation ẍt) =.6ẋt)+.3ẋt 1) xt)+xt 1) 1.7) is asymptotically stable. It is interesting to note that equations 1.6) and 1.7) without delays are unstable. This illustrates a very interesting feature of second-order delay differential equations, i.e. delays may improve asymptotic properties of a given equation, whereas delays in first-order linear equations have mostly destabilizing effects or do not change stability of the model.

3 Several stability tests for non-autonomous linear models with variable delays ẍt) + at)ẋgt)) + bt)xht)) =, 1.8) ẍt)+at)ẋt)+bt)xt)+a 1 t)ẋgt))+b 1 t)xht)) =, 1.9) were obtained in our recent paper [3], under the assumptions: a, a 1, b and b 1 are Lebesgue measurable and essentially bounded functions on [, ); at) a >, bt) b >, t ht) τ, t gt) δ, a t) 4bt), t gt) as)ds < 1/e. Below is the norm in the space L [t, ). Theorem 1.4. [3, Theorem 5.1] If for some t δ a a b b a )+τ + b b a < 1, then equation 1.8) is exponentially stable. Theorem 1.5. [3, Theorem 5.3] Suppose for some t a 1 < 1, a a 1 b then equation 1.9) is exponentially stable. b a + b 1a 1 a b a b < 1, In the present paper, a specially designed substitution transforms linear second order equations into a system, with a further application of the M-matrix method. This and the linearization techniques are used to devise new global stability tests for nonlinear non-autonomous models. These results are explicit, easily verifiable and can be applied to a general class of second order non-autonomous equations. Some of the theorems of the present paper complement our earlier results [, 3], as well as the tests obtained in recent papers [1, 11, 16]. The paper is organized as follows. Section contains stability results for linear second order non-autonomous equations with several delays. To illustrate efficiency of the results obtained each stability test is accompanied by numerical examples. In Section 3 the tests for linear models are applied to nonlinear Lienard-type equations of the second order. Applications incorporate a global stability test for the non-autonomous business cycle model. Section 4 includes the study of bounds and multistability properties for the sunflower model and its generalizations. In particular, sufficient conditions for convergence to one of an infinite number of equilibrium points are presented, and existence of unbounded linearly growing solutions is illustrated. Final remarks are presented in Section 5. Stability tests for linear Lienard equations The technique in this section involves parlaying a second order equation into two first order equations. Consider a linear equation of the second order ẍt)+ a k t)ẋh k t))+ t b k t) ẋs)ds + g k t) c k t)xr k t)) =..1) 3

4 Together with equation.1), for any t we consider the initial condition xt) = ϕt), ẋt) = ψt), t t..) Henceforth, we assume that the following assumptions are satisfied: a1) a i,b i,c i,i = 1,...,m are Lebesgue measurable and essentially bounded on [, ); a) h i,g i,r i are Lebesgue measurable functions, h i t) t,g i t) t,r i t) t, lim h it) =, lim g i t) =,lim t r i t) =, i,j = 1,...,m; t t a3) ϕ and ψ are Borel measurable bounded functions. Definition.1. A function x : R R with locally absolutely continuous on [t, ) derivative ẋ is called a solution of problem.1),.) if it satisfies equation.1) for almost every t [t, ) and equalities.) for t t. We quote a useful lemma that will play a major role in the proofs. Lemma.. [5] Consider the system x i t) = a i t)x i t)+ l ij b k ij t)x jh k ij t)), i = 1,...,m,.3) j=1 where a i t) α i >, b k ijt) L k ij, t h k ijt) σ k ij. If the matrix B = b ij ) m i,j=1, with b ii = 1 lii an M-matrix, then system.3) is exponentially stable. L k ii )/ l ij α i, b ij = L k ij / α i, i j, is We recall that a matrix B = b ij ) m i,j=1 is a nonsingular) M-matrix if b ij, i j and one of the following equivalent conditions holds: either there exists a positive inverse matrix B 1 > or all the principal minors of the matrix B are positive. Further proofs will also require the following lemma. Lemma.3. Consider the system l ij ) t x i t) = a i t)x i t)+ c k ij t)x jgij k t))+dk ij t) x j s)ds, i = 1,...,m,.4) h k ij t) j=1 where a i t) α i >, d k ijt) L k ij, c k ijt) Cij, k t h k ijt) σij, k t gijt) k τ. If the matrix B = b ij ) m i,j=1, with b l ii ) l ij ) ii = 1 L k ii σk ii +Ck ii /α i, b ij = L k ij σk ij +Ck ij /α i, i j, is an M-matrix, then system.4) is exponentially stable. Proof. Let xt) be a solution of.4). Since x j t) are continuous then for any i,j,k and t there exists p k ij t) hk ij t),t) such that x jp k ijt))t h k ijt)) = x j s)ds. h k ij t) Thus x j are solutions of system.3) with b k ij t)x jh k ij t)) being replaced by ck ij x jgij kt)) + d k ij t)t hk ij t))x jp k ij t)). Wehave ck ij t) Ck ij, dk ij t)t hk ij t)) Lk ij σk ij, i j. Theapplication of Lemma. validates the proof. 4 t

5 Note that a different proof of Lemma.3 involves application of the Halanay-type inequalities see, for example, []). To examine the equation we assume ẍt)+at)ẋt)+bt)xt)+ c k t)xh k t)) =.5) < a at) A, < b bt) B, c k t) C k, t h k t) τ. Theorem.4. Suppose at least one of the following conditions holds: 1) B a m 4, C k < b a A a), ) b a A a ), C k < a B. Then equation.5) is exponentially stable. Proof. Substituting ẋ = a x+y,ẍ = a into equation.5), we arrive at ẋ+ẏ Condition 1) yields a the matrix ẋ = a x+y [ a ẏ = at) a ) ] bt) xt) at) a ) bt) a a a 4 B, a c k t)xh k t)) at) a ) yt). 1 ) a A a ) b+ C k 1 is an M-matrix. By Lemma. equation.5) is exponentially stable. If condition ) holds then bt) a at) a ) b a A a ), bt) a at) a ) B a a a ) = B a /4, and the matrix 1 a B a m ) a 4 + C k 1 is an M-matrix. By Lemma. equation.5) is exponentially stable..6) at) a ) bt) a A a ) b. Hence Remark.5. Application of the classical substitution ẋ = y is not useful in our stability investigation, since for the system obtained after this substitution, the matrix B in Lemma. is not an M-matrix. For equation.5) with constant a and b, Ct) C and m = 1 ẍt)+aẋt)+bt)xt)+ct)xht)) =,.7) 5

6 we compare two substitutions ẋt) = λxt)+yt), λ >,.8) and ẋt) = a xt) + yt). By Theorem.4 equation.7) is exponentially stable, if at least one of the following conditions holds: 1) b a 4 ) b > a 4,C < b,,c < a b. Whereas application of.8) by the same token yields a slight improvement: equation.7) is exponentially stable, if at least one of the following conditions holds: 1), ) or 3) b a 4 [,C < a b and the following two intervals have a nonempty intersection a a 4b, a+ a 4b ] [ a a b+c),a+ ] a b+c). Implementation of.8) for equation.5) with nonconstant coefficients at) and bt) will produce a more complicated condition 3). Trading-off these options, we prefer the substitution ẋt) = a xt)+yt). The following numerical examples illustrate the application of Theorem.4. Example.6. Consider the delay equation ẍt)+aẋt)+bxt)+cxt h sin t ) =.9) a) a = 3,b = 1.1,c =.8,h =. Condition 1) of Theorem.4 holds, condition ) does not hold. Equation.9) is asymptotically stable. b) a =,b = 1.1,c =.8,h =. Condition ) of Theorem.4 holds, condition 1) does not hold. Equation.9) is asymptotically stable. c) a =.1,b = 1.5,c = 1.45,h =. Conditions of Theorem.4 do not hold, and equation.9) is unstable. Let us note that in c) the coefficient of the non-delay term exceeds the one of the delayed term: ct) < bt). This is in contrast to the result for the equation with the second derivative omitted aẋt)+bt)xt)+ct)xht)) =, which is exponentially stable if a >, bt) b >, ct) < bt), t τ ht) t for τ >. Consider the equation ẍt)+at)ẋt)+bt)xt)+ c k t)ẋh k t)) =,.1) where < a at) A, < b bt) B, c k t) C k, t h k t) τ. 6

7 Theorem.7. Suppose that at least one of the following conditions holds: 1) B a m 4, C k < b aa a), a ) b a A a ), C k < a B. a Then equation.1) is exponentially stable. Proof. The substitution ẋ = a x+y,ẍ = a ẋ+ẏ into equation.1) yields ẋ = a x+y [ a ẏ = at) a If condition 1) holds, we have a at) a a at) a ) bt) a A a ) ] bt) xt)+ a c k t)yh k t)) at) a ) bt) a 4 ) b. Hence the matrix c k t)xh k t)) ) yt). B, [ 1 ] a a A a ) b+ a C k 1 a a is an M-matrix. By Lemma. equation.1) is exponentially stable. If the inequalities in ) hold then bt) a at) a ) b a A a ), bt) a at) a ) B a a a ) = B a /4. Thus the matrix 1 [ ] a B a a 4 + a C k 1 C k a is an M-matrix. By Lemma. equation.1) is exponentially stable. C k.11) Example.8. Consider the equation ẍt)+aẋt)+bxt)+cẋt h sin t ) =.1) To illustrate Theorem.7, we examined: a) a =.1, b = 1, c =.4, h =. Condition 1) of Theorem.7 holds, condition ) does not hold. Equation.1) is asymptotically stable. b) a = 4, b = 5, c =.7, h =. Condition ) of Theorem.7 holds, condition 1) does not hold. Equation.1) is asymptotically stable. c) a = 1, b = 1.5, c =.8, h =. Conditions of the Theorem.7 do not hold, and, as can be seen from numerical simulations, equation.1) is unstable. Hence, in general, the conditions at) a >,bt) b >,m = 1, ct) < at) are not sufficient for stability of equation.1). 7

8 Consider the equation ẍt)+at)ẋt)+ b k t)xh k t)) =,.13) where < a at) A, < b k b k t) B k, t h k t) τ k. Theorem.9. Suppose at least one of the following conditions holds: 1) B k a 4, a A a) < m b k a B k τ k, ) b k a A a ), B k 1+aτ k ) < a. Then equation.13) is exponentially stable. Proof. With the substitution ẋ = a x+y,ẍ = a ẋ+ẏ into equation.13), we arrive at ẋ = a [ x+y a ẏ = + at) a If condition 1) holds, we have a a at) a ) b k t) a ) t b k t) h k t) ] b k t) xt) [ a xs)+ys) ] ds at) a ) yt). at) a ) b k t) a 4 m B k, A a ) b k. Hence the off-diagonal entries of the matrix 1 [ ] a a A a ) b k + a B k τ k 1 B k τ k a a.14) are non-positive, and the inequalities in 1) yield that it is an M-matrix. By Lemma.3 equation.13) is exponentially stable. Assumption ) implies b k t) a at) a ) b k a A a ), b k t) a at) a ) B k a a a ) = B k a /4, therefore the matrix 1 [ a m ] B k a a 4 + a B k τ k 1 B k τ k a is an M-matrix. By Lemma.3 equation.13) is exponentially stable. 8

9 Corollary.1. Suppose at) a >,b k t) b k >, and at least one of the following conditions holds: 1) b k a 4, m b k 1 aτ k ) >, ) b k a m 4, b k 1+aτ k ) < a. Then equation.13) is exponentially stable. Example.11. Consider the equation ẍt)+aẋt)+bxt h sin t ) =..15) To illustrate Theorem.9, we consider numerical examples: a) a =, b =.9, h =.4. Condition 1) of Theorem.9 holds, condition ) does not hold. Equation.15) is asymptotically stable. b) a =, b = 1.1, h =.4. Condition ) of Theorem.9 holds, condition 1) does not hold. Equation.15) is asymptotically stable. c) a = 1, b = 1.1, h =.5. Conditions of Theorem.9 do not hold. Equation.14) is unstable which can be confirmed numerically. Consider the equation ẍt)+at)ẋt)+bt)xt) = where < a at) A, < b b k t) B, c k t) C k, t h k t) τ k. Theorem.1. Suppose at least one of the following conditions holds: 1) B a m 4, C k τ k < b aa a), a ) b a A a ), C k τ k < a B. a Then equation.16) is exponentially stable. Proof. After rewriting equation.16) in the form c k t)[xt) ixh k t))],.16) ẍt)+at)ẋt)+bt)xt) = we apply the same argument as in the proof of Theorem.7. t c k t) ẋs)ds, h k t) Theorem.4 gives delay-independent stability conditions for equation.5). The following statement contains delay-dependent stability conditions for this equation. Theorem.13. Assume that < a at) A, < b bt)+ c k t) B, c k t) C k,t h k t) τ k and at least one of the conditions of Theorem.1 holds. Then equation.5) is exponentially stable. 9

10 Proof. Rewrite equation.5) in the form ) ẍt)+at)ẋt)+ bt)+ c k t) xt) = t c k t) ẋs)ds. h k t)) The end of the proof is a straightforward imitation of the proof of Theorem.7. 3 Stability tests for nonlinear Lienard equations In this section we examine several nonlinear delay differential equations of the second order which have the following general form ẍt)+ f k t,xp k t)),ẋg k t)))+ with the following initial function l s k t,xh k t))) =, 3.1) xt) = ϕt),ẋt) = ψt),t t, t 3.) where f k t,u 1,u ),k = 1,...,m, s k t,u), are Caratheodory functions which are measurable in t and continuous in all the other arguments, condition a) holds for delay functions p k,g k,h k ; ϕ and ψ are Borel measurable bounded functions. The definition of the solution of the initial value problem 3.1)-3.) is the same as for problem.1),.). We will assume that the initial value problem has a unique global solution on [t, ) for all nonlinear equations considered in this section. Definition 3.1. Suppose the number K is an equilibrium of equation 3.1). We will say that K is an attractor of this equation if for any solution x of the problem 3.1), 3.) we have lim t xt) = K. Theorem 3.. Consider the equation where ẍt)+ft,xt),ẋt))+st,xt))+ s k t,xt),xh k t))) =, 3.3) ft,v,) =,st,) =,s k t,v,) =, < a ft,v,u) A, u < b st,u) B, s k t,v,u) u u C k,u, t h k t) τ. If at least one of the following conditions holds: 1) B a m 4, C k < b a A a ), ) b a A a ), C k < a B, then zero is a global attractor for all solutions of problem 3.3), 3.). 1

11 Proof. Suppose x is a fixed solution of problem 3.3), 3.). Rewrite equation 3.3) in the form ẍt)+at)ẋt)+bt)xt)+ c k t)xh k t)) =, where at) = c k t) = { ft,xt),ẋt)) ẋt), ẋt), a, ẋt) =, bt) = { st,xt)) xt), xt), b, xt) =, { sk t,xt),xh k t))) xh k t)), xh k t)),, xh k t)) =. Hence the function x is a solution of the linear equation ÿt)+at)ẏt)+bt)yt)+ c k t)yh k t)) =, 3.4) which is exponentially stable by Theorem.4. Thus for any solution y of equation 3.4) we have lim yt) =. Since x is a solution of 3.4), we have lim xt) =. t t The previous proof is readily adapted to the proof of the following theorems. Theorem 3.3. Consider the equation where ẍt)+ft,xt),ẋt))+st,xt))+ s k t,xt),ẋh k t))) =, 3.5) ft,v,) =, st,) =, s k t,v,) =, < a ft,v,u) A, u < b st,u) B, s k t,v,u) u u C k,u, t h k t) τ. Suppose at least one of the following conditions holds: 1) B a m 4, C k < b a A a ), a ) b a A a ), C k < a B. a Then zero is a global attractor for all solutions of problem 3.5),3.). Theorem 3.4. Consider the equation ẍt)+ft,xt),ẋt))+ s k t,xh k t)),ẋt)) =, 3.6) where ft,v,) =,s k t,,u) =, < a ft,v,u) A, u < b k s kt,v,u) B k,u, t h k t) τ. v 11

12 Suppose at least one of the following conditions holds: 1) B k a 4, a A a m ) < b k a B k τ k, ) b k a A a ), B k 1+a τ k ) < a. Then zero is a global attractor for all solutions of problem 3.6),3.). Theorem 3.5. Consider the equation where ẍt)+ft,xt),ẋt))+st,xt)) = c k t)xt) xh k t))), 3.7) ft,v,) =,st,) =, < a ft,v,u) A, u < b st,u) B, c k t) C k, u, t h k t) τ k. u Suppose at least one of the following conditions holds: 1) B a m 4, C k τ k < b a A a ), a ) b a A a ), C k τ k < a B a. Then zero is a global attractor for all solutions of problem 3.7),3.). Example 3.6. To illustrate Part ) of Theorem 3.4, consider the equation ẍt) sin xt))ẋt) cos xt))xt.19sin t) =. 3.8) We have m = 1, a = 1.8, A =, b = 1, B = 1., τ =.19; therefore, all conditions of the theorem hold, hence zero is a global attractor for all solutions of equation 3.8). Motivated by model 1.4), consider a generalized Kaldor-Kalecki model ẍt)+ [ αt) βt)p xt)) ] ẋt)+st,xt)) = pxt)) pxht))), 3.9) where α, β are locally essentially bounded functions, s is a Caratheodory function, p is a locally absolutely continuous nondecreasing function, < α αt) α 1, < β βt) β 1, p t) C, α β 1 C >, < b st,u) B, t ht) τ. u Denote a = α β 1 C. Theorem 3.7. Suppose at least one of the following conditions holds: 1) B a 4, Cτ < b a α 1 a ), a ) b a α 1 a ), Cτ < a B. a Then zero is a global attractor for all solutions of problem 3.9),3.). 1

13 Proof. Suppose x is a fixed solution of problem 3.9),3.). There exists a function ξt) such that pxt)) phxt)) = p ξt))xt) xht))). Denote αt) βt)p xt)) = at),p ξt)) = ct). Hence x is a solution of the following equation ÿt)+at)ẏt)+st,yt)) = ct)yt) yht))). 3.1) Since p x) then < α β 1 C at) α 1. Equation 3.1) has a form 3.7) with ft,xt),ẋt)) = at)ẋt),m = 1. All conditions of Theorem 3.5 hold, hence for any solution of 3.1) we have lim t yt) =. Then also lim t xt) =. 4 Sunflower model and its modifications The sunflower equation was introduced in 1967 by Israelson and Johnson in [17] as a model for the geotropic circumnutations of Helianthus annuus and studied in [1, 4, 7]. Historically, it was derived from the following first order delay equation u+ b t τ τ ea1 t/τ) e as/τ sin us)ds =. 4.1) Taking the derivative of 4.1) we arrive at the sunflower equation ẍ+ a τẋ+ b τ sin xt τ) =, 4.) for which evidently the results of the previous section are not applicable. Remark 4.1. It is interesting to note that a non-delayed version of 4.) ẍ+aẋ+bsin xt) =, 4.3) has a long history, see, for example, [5]). However, many important questions for delayed model 4.) are still left unanswered. Consider a generalization of model 4.1) with the initial conditions du dt +b ht) Kt,s)sin us) ds =, 4.4) ut) = ϕt), t, 4.5) under the following assumptions: b1) ht) t τ for some τ > ; b) K, ) is Lebesgue measurable, Kt,s), there exists a > such that Kt,s) 1 { τ exp a } ht) τ t s τ) and Kt,s) ds dt = ; b3) ϕ : [,] R is a continuous bounded function. 13

14 Theorem 4.. Suppose that b1)-b3) hold, b > and the characteristic equation λ τ aλ+be λτ = 4.6) has a positive root λ >. Then any solution of 4.4)-4.5) with the initial conditions satisfying either ϕt) πk,πk +π), k N, or ϕt) πk π,πk), k N, together with ϕt) πk ϕ)e λ t, t <, tends to πk as t. Moreover, for ϕt) πk,πk + π) the solution is monotone decreasing, while for ϕt) πk π,πk) it is monotone increasing. Proof. First assume that ϕt),π), t, u is a solution of 4.4). Let us prove that i) ut) is positive and non-increasing function; ii) ut) satisfies the inequality iii) ut) tends to zero as t. ut) u)e λ t, t ; 4.7) Denote ut) = ϕt) for t as well, then by the assumptions of the theorem, ut) u)e λ t, t <. We start verifying ii) by induction. First, we prove that ut) u)e λ t for t [,τ], and then proceed to any segment [nτ,n+1)τ]. In the inequalities below, we use the estimates of K in b), the fact that sin u u for u > and ut) u)e λ t for t < to evaluate the derivative of u on [,τ]: du dt = b ht) ϕ) b τ t τ t τ Kt,s)sin us)) ds b exp { a } τ t s τ) e λs ds = ϕ) b τ exp { a τ t τ) } t τ exp{ a τ λ b = ϕ) a λ τ e λ t τ) = ϕ)λ e λt, Kt, s)ϕs) ds since a λ τ = b λ e λτ by 4.6). Since u t) u)λ e λt, for any s, t satisfying t s τ, the solution on [t,s] is not below the curve us) = ut)e λs t) on [,τ], and us) ut)e λs t). Taking s = τ, we obtain uτ) ut)e λτ t), or ut) uτ)e λt τ), t [,τ]. Hence uτ)e λ t τ) is an upper bound of ut) on [,τ], it is also a bound on,τ] since u) uτ)e λ τ and for t, Thus ut) uτ)e λ t τ) is valid on,τ]. ) } s ut) u)e λ t uτ)e λ t τ). 4.8) ds 14

15 Consider further the initial problem with a shifted initial point t = τ instead of t =, we get the same estimate as in 4.8) for any t,nτ] by induction. Hence, ut) unτ)e λ t nτ) u)e λ t >, [nτ,n+1)τ], and the induction step proves 4.7) and justifies ii). Thus the solution is positive for any t. From the form of the equation and nonnegativity of K and u,π) follow that the solution is also non-increasing, which justifies i). Since u is non-increasing for t and positive there is lim t ut) = d. Assuming d > we obtain from ht) Kt,s) ds dt = in b) that lim t ut) =, which is a contradiction, thus iii) is also valid. A similar process proves the case ϕt) π,). If ϕt) πk π,πk), we apply the same argument to u πk. Note that sharp conditions when all solutions of characteristic equation 4.6) have positive real parts can be found in [7, Lemma 3.1, p. 47]. Corollary 4.3. Let τ < a 4b e a/ 4.9) and ϕt) πk ϕ)e λ t, t <, then any solution of 4.4)-4.5) with the initial conditions satisfying ϕt) πk,πk +π), k N, is monotone decreasing and tends to πk as t. Any solution with ϕt) πk π,πk), k N tends to πk as t. Proof. Let fλ) = τλ aλ + be λτ, then f) = b >. Inequality 4.9) implies fa/τ)) = a /4τ) + be a/ <, so equation 4.6) has a positive solution. We invokef Theorem 4. to conclude the proof. The following example illustrates that conditions b1)-b3) do not guarantee boundedness of the solutions of equation 4.4) with the generalized kernel. Example 4.4. Let a = 1 ) 4 3 ln, b =, τ = π, π Kt,s) = { 1, t [k 1)π,k +1)π], s [k 3)π,k )π], 4, t [k 1)π,k +1)π], s [k 3)π,k )π]. Then obviously Kt,s) = for s > t π = t τ, and also for t s > 4π. The exponential estimate has the form Kt,s) 1 π e 1 3π ln4/π)t s π) = 1 ) 4 t s π)/3π), π π but as t s π 3π whenever Kt,s), the right-hand side is not less than 1 ) 4 1 = 1 π π 4, thus Kt,s) has an exponential estimate as in b). Further, ut) = t is an unbounded solution of 15

16 4.4). In fact, let ut) = t, t [ π,π]. Then for t [π,3π] we have du dt = 1 sin t) dt = 1, π 4 so ut) = t on [ π,3π]. Due to the periodicity of the sine function and K, we have du 1. Thus dt the solution is a linear function ut) = t and it is unbounded. In the following theorem we will prove that for a non-autonomous case the solution of the sunflower equation is bounded by a linear function. Consider the non-autonomous sunflower equation ẍt) + at)ẋt) + bt)sin xht)) =. 4.1) Theorem 4.5. Suppose at) a >, bt) b. For any solution xt) of equation 4.1) we have the estimates xt) xt ) + ẋ) + b ) t, ẋt) ẋ) + b. a a Proof. Denote ẋ = y,ft) = bt)sin xht)), where ft) b. Then ẏt) + at)yt) + ft) =, hence yt) = y)+ t t e s aτ)dτ fs)ds. Then ẋt) ẋ) + xt) = x)+ t t e a t s) fs)ds ẋ) + b a, ẋs)ds, xt) xt ) + ẋ) + b ) t. a Local stability conditions for equation 4.1) one can find in the following theorem. The following lemma is a corollary of [1, Theorem 8.3] Lemma 4.6. Suppose that ordinary differential equation ẍt)+at)ẋt)+bt)xt), at),bt), has a positive fundamental function, then the equation ẍt)+at)ẋt)+b 1 t)xt), where b 1 t) bt) also has a positive fundamental function. Theorem 4.7. Suppose < a at) A, < b bt) B, t ht) τ and at least one of the following conditions hold: 1) B a 4, a A a) < b abτ, ) b a A a ), B1+aτ) < a. Then any equilibrium xt) = kπ,k =,... of equation 4.1) is locally asymptotically stable. Any equilibrium xt) = k +1)π,k =,... is not asymptotically stable. 16

17 Proof. For the equilibrium xt) = kπ, the linearization of equation 4.1) has the form ÿt)+at)ẏt)+bt)yht)) =, which is exponentially stable by Theorem.9. It is well known see, for example, [4]) that exponential stability of a linearized equation implies asymptotic stability of the nonlinear equation, in our case equation 4.1). For the equilibrium xt) = k +1)π, the linearized equation for 4.1) has the form Consider now the ordinary differential equation ÿt) + at)ẏt) bt)yht)) =. 4.11) zt) + at)żt) =. 4.1) The fundamental function of equation 4.1) the solution of initial value problem with z) =,z ) = 1) has the form zt) = t e s aτ)dτ ds, which is a positive function for t > with a nonnegative derivative. By Lemma 4.6, for the fundamental function yt) of equation 4.11) we have yt) >, y t) for t >. Hence yt) does not tend to zero, and thus equation 4.11) is not asymptotically stable. 5 Concluding Remarks The technique of reduction of a high-order linear differential equation to a system by the substitution x k) = y k+1 is quite common. However, this substitution does not depend on the parameters of the original equation, and therefore does not offer new insight from a qualitative analysis point of view. Instead, we proposed a substitution which exploits the parameters of the original model. By using that approach, a broad class of the second order non-autonomous linear equations with delays was examined and explicit easily-verifiable sufficient stability conditions were obtained. There is a natural extension of this approach to stability analysis of high-order models. For the nonlinear second order non-autonomous equations with delays we applied the linearization technique and the results obtained for linear models. Our stability tests are applicable to some milling models, e.g. models 1.) and 1.3), and to a non-autonomous Kaldor Kalecki business cycle model. Several numerical examples illustrate the application of the stability tests. We suggest that a similar technique can be developed for higher order linear delay equations, with or without non-delay terms. For a non-autonomous version of a classical sunflower model, we verified that the derivative is bounded and thus the solution has a linear bound. Example 4.4 illustrates the existence of an unbounded linearly growing solution for the generalized sunflower equation. We also obtained sufficient conditions under which a solution tends to one of the infinite number of the equilibrium points. Solution of the following problems will complement the results of the present paper: 17

18 1. In all stability conditions obtained, we used lower and upper bounds of the coefficients and the delays. It is interesting to obtain stability conditions in an integral form, for instance, in the assumptions of Theorem.9 replace the term aτ k by, generally, a smaller term t h k t) as) ds.. Apply the technique used in the paper to examine delay differential equations of higher order. 3. Is it possible to generalize Theorem 4. to the case when the initial function ϕt) πk π,πk +π) and characteristic equation 4.6) has a solution with a positive real part? 4. Establish necessary stability conditions for the equations considered in this paper. 5. For the sunflower equation and its modifications establish set of conditions to guarantee boundedness of all solutions. 6 Acknowledgments The authors also would like to extend their appreciation to the anonymous referees for helpful suggestions which have greatly improved this paper. L. Berezansky was partially supported by Israely Ministry of Absorption, E. Braverman was partially supported by the NSERC grant RGPIN/ , L. Idels was partially supported by a grant from VIU. The authors are grateful to the anonymous reviewer whose comments significantly contributed to the presentation of the paper. References [1] R. P. Agarwal, L. Berezansky, E. Braverman, A. Domoshnitsky, Nonoscillation Theory of Functional Differential Equations with Applications. Springer, New York, 1. [] L. Berezansky, J. Diblik, Z. Smarda, Positive solutions of second-order delay differential equations with a damping term, Comput. Math. Appl. 6 1) [3] L. Berezansky, E. Braverman, A. Domoshnitsky, Stability of the second order delay differential equations with a damping term, Differ. Equ. Dyn. Syst. 16 8) [4] L. Berezansky, E. Braverman, On stability of some linear and nonlinear delay differential equations. J. Math. Anal. Appl ) [5] L. Berezansky, E. Braverman and L. Idels, New global exponential stability criteria for nonlinear delay differential systems with applications to BAM neural networks, Appl. Math. Comput ), [6] T. A. Burton, Second order boundedness criteria, Ann. Mat. Pura Appl. 4) ) [7] T. A. Burton, B. Zhang, Boundedness, periodicity, and convergence of solutions in a retarded Lienard equation, Ann. Mat. Pura Appl. 4) )

19 [8] T. A. Burton, Fixed points, stability, and exact linearization, Nonlinear Anal. 61 5) [9] T. A. Burton, Stability by Fixed Point Theory for Functional Differential Equations, Dover Publications, Inc., Mineola, NY, 6. [1] B. Cahlon, D. Schmidt, Stability criteria for certain second-order delay differential equations with mixed coefficients, J. Comput. Appl. Math. 17 4) [11] B. Cahlon, Unconditional stability for certain delay differential equations, Dynam. Syst. Appl ), no. 4, [1] A. Casal, A. Somolinos, Forced oscillations for the sunflower equation, entrainment, Nonlinear Anal. Theor. Meth. Appl., 6 198) [13] Z. Dombovari, D. Barton, R. Wilson, G. Stepan, On the global dynamics of chatter in the orthogonal cutting model, Int. J. Non-Linear Mech ) [14] J. Hale, S. M. Verduyn Lunel, Introduction to Functional-Differential Equations, Applied Mathematical Sciences, 99. Springer-Verlag, New York, [15] G. Gabish, H. Lorenz, Business Cycle Theory: A Survey of Methods and Concepts, Lecture Notes in Economics and Mathematics, Vol. 83, Springer-Verlag, Germany, [16] I. Gyori, F. Hartung, Asymptotically exponential solutions in nonlinear integral and differential equations, J. Differ. Equ. 49 1), [17] D. Israelson, A. Johnson, A theory of circumnutatins of Helianthus annus, Physiol. Plants 1967), [18] V. Kolmanovskii, A. Myshkis, Applied theory of functional-differential equations, Mathematics and its Applications Soviet Series), 85. Kluwer Academic Publishers Group, Dordrecht, 199. [19] P. Kim, J. Jung, S. Lee, J. Seok, Stability and bifurcation analyses of chatter vibrations in a nonlinear cylindrical traverse grinding process, J. Sound Vib ) [] L. Liu, T. Kalmár-Nagy, High-dimensional harmonic balance analysis for second-order delaydifferential equations, J. Vib. Control 16 1) [1] B. Liu, L. Huang, Boundedness of solutions for a class of Lienard equations with a deviating argument, Appl. Math. Lett. 1 8) [] E.Liz, S. Trofimchuk, Existence and stability of almost periodic solutions for quasilinear delay systems and the Halanay inequality, J. Math. Anal. Appl. 48 ) [3] Long, Wei; Zhang, Hong-Xia Boundedness of solutions to a retarded Lienard equation, Electron. J. Qual. Theory Differ. Equ. 1, no. 4, 9 pp. [4] M. Lizana, Global analysis of the sunflower equation with small delay, Nonlinear Anal. Theor. Meth. Appl )

20 [5] J. Mawhin, Seventy-five years of global analysis around the forced pendulum equation, in Proc. of the Equadiff 9 Conf. on Diff. Equat. and Their Appl. Brno, Czech Republic, 5-9 August 1997; 1998) [6] D. Pi, Study the stability of solutions of functional differential equations via fixed points, Nonlinear Anal ) [7] A. Somolinos, Periodic solution of the sunflower equation ẍ+a/r)ẋ+b/r)sin xt r) =, Quart. Appl. Math ) [8] G. Stepan, Delay effects in the human sensory system during balancing, Philos. Trans. R. Soc. Lond. Ser. A Math. Phys. Eng. Sci ) [9] R. Szalai, G. Stepan, Period doubling bifurcation and center manifold reduction in a timeperiodic and time-delayed model of machining, J. Vib. Control 16 1) [3] M. Szydlowski, A. Krawiec, The stability problem in the Kaldor Kalecki business cycle model, Chaos Solitons Fractals 5 5) [31] C. Tunq, Stability to vector Lienard equation with constant deviating argument, Nonlinear Dynam ) [3] B. Zhang, Necessary and sufficient conditions for boundedness and oscillation in the retarded Lienard equation, J. Math. Anal. Appl. 1996) [33] B. Zhang, Boundedness and stability of solutions of the retarded Lienard equation with negative damping, Nonlinear Anal. 1993) [34] B. Zhang, On the retarded Lienard equation, Proc. Amer. Math. Soc ) [35] B. Zhang, Fixed points and stability in differential equations with variable delays, Nonlinear Anal. 63 5) e33 e4. [36] Q. Xie, Q. Zhang, Stability predictions of milling with variable spindle speed using an improved semi-discretization method, Math. Comput. Simulat. 85 1) [37] M. Wan, W. Zhang, J. Dang, Y. Yang, A unified stability prediction method for milling process with multiple delays, Int. J. Machine Tools Manufacture 5 1) [38] J. Yan, Existence of oscillatory solutions of forced second order delay differential equations, Appl. Math. Lett. 4 11)