Exponential stability of Euler-Bernoulli beam under boundary controls in rotation and angular velocity

Alemdar Hasanov
2023-07-04
Abstract:This paper addresses the analysis of a boundary feedback system involving a non-homogeneous Euler-Bernoulli beam governed by the equation $m(x)u_{tt}+\mu(x)u_{t}$$+\left(r(x)u_{xx}\right)_{xx}=0$, subject to the initial $u(x,0)=u_0(x)$, $u_t(x,0)=v_0(x)$ and boundary conditions $u(0,t)=0$, $\left (-r(x)u_{xx}(x,t)\right )_{x=0}=-k^{-}_r u_{x}(0,t)-k^{-}_a u_{xt}(0,t)$, $u(\ell,t)=0$, $\left (-r(x)u_{xx}(x,t)\right )_{x=\ell}=-k^{+}_r u_{x}(\ell,t)-k^{+}_a u_{xt}(\ell,t)$, with boundary control at both ends resulting from the rotation and angular velocity. The approach proposed in this study relies on the utilization of regular weak solutions, energy identity, and a physically motivated Lyapunov function. By imposing natural assumptions concerning physical parameters and other inputs, which ensure the existence of a regular weak solution, we successfully derive a uniform exponential decay estimate for the system's energy. The decay rate constant featured in this estimate is solely dependent on the physical and geometric properties of the beam. These properties encompass crucial parameters such as the viscous external damping coefficient $\mu(x)$, as well as the boundary springs $k^{-}_r,k^+_r $ and dampers $k^{-}_a,k^+_a$. To illustrate the practical effectiveness of our theoretical findings, numerical examples are provided. These examples serve to demonstrate the applicability and relevance of our derived results in real-world scenarios.
Optimization and Control,Mathematical Physics,Analysis of PDEs
What problem does this paper attempt to address?
This paper aims to solve the problem of exponential stability analysis of non - uniform Euler - Bernoulli beams under rotational and angular velocity boundary control. Specifically, the research focuses on how to ensure the exponential stability of the system energy described by the equation \(m(x)u_{tt}+\mu(x)u_t+(r(x)u_{xx})_{xx} = 0\) under minimal conditions. Here, \(m(x)\) represents the mass density of the beam multiplied by the cross - sectional area, \(\mu(x)\) is the viscous damping coefficient, and \(r(x)\) is the flexural rigidity. ### Research Background and Problem In many practical applications, such as submarine pipelines and long bridges, these structures can be regarded as elastic beams with rotational and angular velocity boundary control at both ends. However, existing research usually assumes relatively ideal support conditions for pipelines or bridges, for example, assuming that the pipeline is located on a rigid seabed without penetration. This assumption is not always valid in practice. Therefore, a more realistic model is needed to consider the influence of rotational springs and angular velocity dampers at both ends of the beam. ### Research Method The author uses the theory of regular weak solutions, energy identities, and Lyapunov functions to derive the exponential decay estimate of the system energy. By naturally assuming physical parameters and other inputs, the existence of regular weak solutions is ensured, and the energy decay rate of the system is successfully obtained. This decay rate only depends on the physical and geometric parameters of the beam, including the viscous external damping coefficient \(\mu(x)\) and the stiffness \(k_r^-\), \(k_r^+\), \(k_a^-\) and \(k_a^+\) of the boundary springs and dampers. ### Main Results 1. **Energy Identity**: \[ E(t)+\int_0^t\int_0^\ell\mu(x)u_t^2(x,\tau)dx d\tau=E(0)-k_a^-\int_0^t u_x^2(0,\tau)d\tau - k_a^+\int_0^t u_x^2(\ell,\tau)d\tau \] where \(E(t)\) is the total energy of the system, defined as: \[ E(t)=\frac{1}{2}\int_0^\ell\left[m(x)u_t^2(x,t)+r(x)u_{xx}^2(x,t)\right]dx+\frac{1}{2}k_r^-u_x^2(0,t)+\frac{1}{2}k_r^+u_x^2(\ell,t) \] 2. **Exponential Stability Estimate**: \[ E(t)\leq M e^{-\sigma t}E(0) \] where \(M\) and \(\sigma\) are respectively: \[ M = \frac{1+\beta_1\lambda}{1-\beta_0\lambda},\quad\sigma=\frac{2\lambda}{1+\beta_1\lambda} \] and \(\lambda\) satisfies \(0 < \lambda<\min\left(\frac{1}{\beta_0},\frac{\gamma m_0}{2m_1}\right)\). ### Numerical Examples To verify the validity of the theoretical results, the paper provides numerical examples. These examples show the behavior of system energy decay under different damping parameters. In particular, when the stiffness \(k_a^-\) and \(k_a^+\) of the rotational dampers are not zero, the change in the energy decay rate \(\sigma\). ### Conclusion This research proposes a method based on regular weak solutions, energy identities, and Lyapunov functions for analyzing the exponential stability of Euler - Bernoulli beams under rotational and angular velocity boundary control. By naturally assuming physical parameters and other inputs, the existence of regular weak solutions is ensured, and successfully...