Stability of thermoelastic Timoshenko beam with suspenders and time-varying feedback

This paper considers a one-dimensional thermoelastic Timoshenko beam system with suspenders, general weak internal damping with time varying coefficient, and time-varying delay terms. Under suitable conditions on the nonlinear terms, we prove a general stability result for the beam model, where exponential and polynomial decay are special cases. We also gave some examples to illustrate our theoretical finding.


Introduction
In this paper, we consider a thermoelastic Timoshenko beam with suspension cables, weak internal damping, and a time-varying delay damping of the form ⎧ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎩ ρu tt (x, t)αu xx (x, t)λ(ϕu)(x, t) + γ 1 a(t)g 1 (u t (x, t)) + γ 2 a(t)g 2 (u t (x, tτ (t))) = 0, ρ 1 ϕ tt (x, t)k(ϕ x + ψ) x (x, t) + λ(ϕu)(x, t) + γ 3 ϕ t (x, t) = 0, ρ 2 ψ tt (x, t)bψ xx (x, t) + k(ϕ x + ψ)(x, t)mθ x (x, t) = 0, ρ 3 θ t (x, t)βθ xx (x, t)mψ xt (x, t) = 0, (1.1) where (x, t) ∈ (0, 1) × (0, ∞), ϕ represents the transverse displacement (in vertical direction) of the beam cross section, ψ is the angle of rotation of a cross-section. The vertical displacement of the vibrating spring (the cable) is represented by the function u, θ depicts the thermal moment of the beam, λ > 0 is the common stiffness of the string, and α > 0 is the elastic modulus of the string (holding the cable to the deck). The positive constants ρ, ρ 1 , ρ 2 are the density of the mass material of the cable, the mass density, and the moment of mass inertia of the beam, respectively. Also, b, k, β, m represent the rigidity coefficient of the cross-section, the shear modulus of elasticity, the thermal diffusivity, and the coupling coefficient which depends on the material properties, respectively. The function τ (t) > 0 is the time-varying delay, γ 1 and γ 2 are real positive damping constants, g 1 and g 2 are the damping functions, and a(t) is a nonlinear weight function. We supplement (1.1) with the boundary conditions (1.2) and the initial data ψ(x, 0) = ψ 0 (x), θ (x, 0) = θ 0 (x), in (0, 1), ψ t (x, 0) = ψ 1 (x), in (0, 1), u t (x, tτ (0)) = f 0 (x, tτ (0)), in (0, 1) × (0, τ (0)). (1. 3) The stability of the above thermoelastic Timoshenko system with suspension cables would be our penultimate focus in this work. The Timoshenko beam model is arguably very popular and most used when the vibration of a beam exhibits significant transverse shear strain. A model to describe this phenomenon was introduced by Timoshenko [35] in 1921, see also [15,18]. The nonlinear vibration of suspension bridges have captured the attention of different researchers and a number of research articles were written on the topic. The somewhat unpredictable large oscillations of suspension bridges have been modeled in diverse ways, one may see [1,14,25]. In any attempt to adequately describe the complicated dynamics of a suspension bridge, a robust model would be one with a considerable amount of degrees of freedom. Without prejudice, some simplified models have been considered, but do not account for a number of realistic behavior of suspension bridges, e.g., torsional oscillations. Of an advantage is the fact that rigorous mathematical analysis is easily carried out with such simpler models. A typical simplified model is the onedimensional vibrating beam model, where torsional motion is neglected by ignoring sectional dimensions of the beam when they are negligible compared to length of the beam. The emergence of string-beam systems which model a nonlinear coupling of a beam and main cable (the string) were born out of the pioneering works of Lazer, McKenna, and Walter [23,25,26] (see also [7] and its references). Though initially modeled through the classic Euler-Bernoulli beam theory, the Timoshenko beam theory is also proven to perform better in predicting a beam response to vibrations than a model based on the classical Euler-Bernoulli beam theory. Indeed, the Timoshenko beam theory takes into account both rotary inertia and shear deformation effects, these are often neglected when applying Euler-Bernoulli beam theory.
In the Timoshenko beam with suspenders which is modeled by (1.1), the suspenders are cables which are elastic in nature and are attached to the beam with elastic springs. The temperature dissipation here is assumed to be governed by the Fourier law of heat conduction. For a(t) ≡ 1, g 1 (s) ≡ s and γ 2 ≡ 0, g 2 ≡ 0 in system (1.1), Bochichio et al. [6] proved a well-posedness and an exponential stability result. A number of works have been done on different thermoelastic Timoshenko models without suspenders (see [10,12,16,17,28] and references in them). Time delays occur in systems modeling many phenomena in areas such as biosciences, medicine, physics, robotics, economics, chemical, thermal, and structural engineering. These phenomena depend on both present and some past history of occurrences, see [8,9,13,21,34] and the examples therein. In the case of constant delay and constant weight, the delay term usually accounts for the past history of strain, only up to some finite time τ (t) ≡ τ .
A step further involves results in the literature about constant weights (γ 1 a(t) ≡ γ 1 , γ 2 a(t) ≡ γ 2 constants) and time-varying delay τ (t). Works presenting the exponential stability result for wave equation with boundary or internal time-varying delay appeared in Nicaise et al. [32,33]. Enyi and Mukiawa in [11] presented a general decay result for a viscoelastic plate equation under the condition |γ 2 | < |γ 1 | √ (1d). Furthermore, in [4,24], the authors presented some existence and stability results for wave equation with internal time-varying delay and time-varying weights; and for suspension bridge models, see Mukiawa [3,27,29,30].
Motivated by the works in [3,6,29], in the current paper, we are concerned with the stability result for the thermoelastic Timoshenko system with suspension cables, timevarying internal feedback, and time-varying weight given in (1.1)-(1.3). The result in [6] is a particular case of our result in this paper.
We arrange this paper in the following manner. In Sect. 2, we state the needed assumptions. In Sect. 3, we present the proof of some technical and needed lemmas for our main result. In the last Sect. 4, we present and prove our main stability result. Throughout this paper, c and c i , i = 1, 2, . . . , are generic positive constants, which are not necessarily the same from line to line.

Functional settings and assumptions
In this section, we state some needed assumptions on the damping coefficients, nonlinear functions, and the time-varying delay. As in [5,32,33], we assume the following conditions: (A 1 ) Function a : [0, +∞) → (0, +∞) is a nonincreasing C 1 -function such that there exists a positive constant C satisfying +∞ 0 a(t) dt = +∞.
Similarly, as in Nicaise and Pignotti [31], we introduce the following change of variable: It follows that subjected to the boundary conditions and initial data We introduce the following spaces: For completeness, we state without proof the existence and uniqueness result for problem (1.1)-(1.3). The result can be established using the Faedo-Galerkin approximation method, see [5] or standard nonlinear semigroup method, see [19,20].
be given such that

Technical lemmas
In this section, we prove some important lemmas which will be essential in establishing the main result. Letμ be a positive constant satisfying and set The energy functional of system (2.14)-(2.16) is defined by Proof Multiplying (2.14) 1 by u t , (2.14) 2 by ϕ t , (2.14) 3 by ψ t , and (2.14) 4 by θ , integrating the outcome over (0, 1), and applying integration by parts and the boundary conditions, we get Adding (3.4)-(3.7), we arrive at Now, multiplying equation (2.14) 5 by μ(t)g 2 (z(x, σ , t)) and integrating over (0, 1) × (0, 1), we obtain On account of (2.6), we can write Therefore, (3.9) becomes (3.11) It follows that Recalling the definition of the energy functional (3.2), and adding (3.8) and (3.12), we obtain (3.13) On the account of (A 1 ) and (2.5), we get (3.14) Now, we consider the convex conjugate of G defined by Therefore, on account of (2.5) and (3.17), we have

Lemma 3.4 The functional
satisfies, along the solution of system (2.14)-(2.16), the estimate (3.26) Using the last equation in (2.14), we can express the last term on the right hand-side of (3.26) as

Stability result
In this section, we are concerned with the main stability result, and is stated as follows. where ds and χ 0 (t) = tχ (r 0 t).

Examples
We end this section by giving some examples to illustrate the obtained result. Let be a strictly increasing function such that g 0 (0) = 0 and, for some positive constants c 1 , c 2 and r, the function g 1 satisfies g 0 |s| ≤ g 1 (s) ≤ g -1 0 |s| , ∀|s| ≤ r, c 1 s 2 ≤ sg 1 (s) ≤ c 2 s 2 , ∀|s| ≥ r. It follows that χ is a C 2 -strictly convex function on (0, r] when g 0 is nonlinear and therefore satisfies condition (A 2 ). Now, we give some examples of g 0 such that g 1 satisfies (5.1) near 0.

Conclusion
In this work, we obtained some general decay results for a thermoelastic Timoshenko beam system with suspenders, general weak internal damping, time-varying coefficient, and time-varying delay terms. The damping structure in system (2.14)-(2.16) is sufficient enough to stabilize the system without any additional conditions on the coefficient parameters as it is the case with many Timoshenko beam systems in the literature. The result of the present paper generalizes the one established in Bochichio et al. [6] and allows a large class of functions that satisfy condition (A 2 ). We also gave some examples to illustrate our theoretical finding.