Hermite Finite Element Method for a Class of Viscoelastic Beam Vibration Problem

Beam equation can describe the deformation of beams and reflect various bending problems and it has been widely used in large engineering projects, bridge construction, aerospace and other fields. It has important engineering pra-ctice value and scientific significance for the design of numerical schemes. In this paper, a scheme for vibration equation of viscoelastic beam is developed by using the Hermite finite element. Based on an elliptic projection, the errors of semi-discrete scheme and fully discrete scheme are analyzed respectively, and the optimal L 2 -norm error estimates are obtained. Finally, a numerical example is given to verify the theoretical predictions and the validity of the scheme.


Introduction
Consider the following viscoelastic beam vibration problem: where ( ) , u x t represents the transverse displacement of the beam, also known as the deflection, EI is the bending stiffness, 0 ρ is the density, S is the cross-sectional area, 0 µ > is the damping coefficient, L is the constant, ( ) ( ) various bending problems. It has been widely used in large engineering projects, bridge construction, aerospace and other fields. Its theoretical and numerical methods have been paid much attention. In [1], the viscoelastic dynamic beam model was studied by using the finite element method and Newmark method. Pierro [2] not only studied the transverse vibration of the viscoelastic beam, but also enlightened the damping characteristics of the structure. Huang, Kouami, et al. [3] [4] presented the finite element scheme for viscoelastic sandwich beam in 2019 and 2020 respectively. Snchez [5] proposed a finite element formulation for thick arbitrarily laminated beams. In [6] the Faedo-Galerkin method was used to solve beam vibration equation with nonlinear tension, a viscoelastic damping and distributed delay term. Wu [7] studied the viscoelastic string-beam coupled system in detail. In [8], the one-dimensional finite element method and analytical method were used to study Euler-Bernoulli beams equations resting on multi-layered viscoelastic soil. A new numerical method was considered in [9] for the equations of fractional-order viscoelastic Euler-Bernoulli beams. Among the aforementioned methods, the finite element method has the characteristics of flexible mesh generation, wide application area and high calculation accuracy. The cubic Hermite element can guarantee the continuity of the first derivative of the interpolation function, and at the same time, the displacement value u and the angle value u' can be calculated [10]. However, there is little work about the Hermite finite element method for viscoelastic beam vibration problem. Therefore, in this paper, a finite element calculation scheme is established for problem (1) based on Hermite element. Using elliptic projection operator, the error analysis is carried out and the convergence order of the scheme is proved to be ( ) 4 O h τ + . Finally, a numerical example is given to demonstrate the correctness of the theoretical analysis and the validity of the scheme.
The rest of the paper is arranged as follows. In Section 2, we give the semidiscrete finite element scheme and fully discrete finite element scheme for problem (1) and conduct convergence analysis. A numerical example is given in Section 3. In section 4, we present conclusions and some future works.

Finite Element Approximation
When finite element method is used to solve this kind of equation, time variable t is usually treated as parameter. Let

Error Estimation of Semi-Discrete Finite Element Schemes
Let H I constituted by piecewise cubic Hermite type polynomials on h I . We can define the semi-discrete Galerkin finite element approximation of (4): find In order to estimate the error, we first define the elliptic projection A well known error estimate (see [11]) is the following.
where C is a constant, independent of L and h.
Proof. We can obtain easily by Lemma 2.1 that Combining (3) and (4), we get the following error equation Here the first and the second are nonnegative, setting By the triangle inequality, we deduce that Combing (8) and (13), we complete the proof.

Error Estimation of Fully Discrete Finite Element Schemes
Next, we will estimate the right end of equation (18) one by one. By using Taylor expansion, we deduce Using Schwarz inequality we can verify Then, we observe that (20), (24) and Lemma 2.1 yield      We use the finite element fully discrete scheme to solve the problem. In Table  1, by choosing 4 h τ = , we give the L 2 -error and the order of spatial convergence.
The results show that the space convergence order of the scheme is fourth order. Surface plots of the exact solution and the numerical solution are presented in Figure 1, Figure 2.

Conclusions
Based on Hermite finite element method, the vibration of viscoelastic beam is analyzed numerically. Semi-discrete and fully discrete Hermite finite element formats are given. By means of a numerical example, the L 2 -error and the convergence order between the exact solution of the original equation and the finite element solution are given. The validity of the proposed scheme in theoretical applications is verified.
Since the equations in this paper are one-dimensional in space, high-order finite element method such as finite element method for two-dimensional plate vibration problem will be the direction of our future research.

Funding
The work is supported by National Natural Science Foundation, China (11501335).