Consider the following viscoelastic beam vibration problem:
where represents the transverse displacement of the beam, also known as the deflection, EI is the bending stiffness, is the density, S is the cross-sectional area, is the damping coefficient, L is the constant, and are smooth functions which are given known.
Beam is the most common component of mechanical equipment and engineering structure. Beam equation can describe the deformation of beams and reflect 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 , the viscoelastic dynamic beam model was studied by using the finite element method and Newmark method. Pierro  not only studied the transverse vibration of the viscoelastic beam, but also enlightened the damping characteristics of the structure. Huang, Kouami, et al.   presented the finite element scheme for viscoelastic sandwich beam in 2019 and 2020 respectively. Snchez  proposed a finite element formulation for thick arbitrarily laminated beams. In  the Faedo-Galerkin method was used to solve beam vibration equation with nonlinear tension, a viscoelastic damping and distributed delay term. Wu  studied the viscoelastic string-beam coupled system in detail. In , 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  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 . 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 . 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 semi-discrete 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.
2. Finite Element Approximation
When finite element method is used to solve this kind of equation, time variable t is usually treated as parameter. Let , introducing the Sobolev space . For any and for fixed t, multiply both sides of equation (1) by and integrate, we have
Integrate by parts with respect to and pay attention to the boundary
conditions and the value of . Based on (2), we can obtain the weak formulation of (1): find , satisfying
2.1. Error Estimation of Semi-Discrete Finite Element Schemes
Let be an uniform partition of the interval ,
, , . Let be a finite-dimensional subspace of
constituted by piecewise cubic Hermite type polynomials on . We can define the semi-discrete Galerkin finite element approximation of (4): find satisfying
In order to estimate the error, we first define the elliptic projection such that
A well known error estimate (see  ) is the following.
Lemma 2.1. denotes the degree of finite element space , for any , we have
Now, we prove the estimate for the error between the solutions of the semidiscrete and continuous problems. Denote and , the first important conclusion of this paper is as follows.
Theorem 2.1. Let u and be the solutions of (1) and (4), respectively. . Then
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
belongs to . Thus, choosing in (9), we conclude that
Using ε-Cauchy inequality gives
Integrate the above inequality over
Here the first and the second are nonnegative, setting , yields
It follows as above that
By the triangle inequality, we deduce that
Combing (8) and (13), we complete the proof.
2.2. Error Estimation of Fully Discrete Finite Element Schemes
Let be an uniform partition of the interval ,
, , . Replacing the first and the second time de-
rivative in equation (4) by backward difference quotient and central difference quotient respectively, full-discrete finite element form of problem (1) can be define: find such that
Denote and in analogy with definitions before. We obtain theorem 2.2 which estimates the error between the solutions of the fully discrete and continuous problems.
Theorem 2.2. Let u and be the solutions of (1) and (15), respectively. . Then
where C is a constant, independent of T and .
Proof. Minusing (3) and (15), we get
Setting in (17) and using ε-Cauchy inequality, we observe that
Next, we will estimate the right end of equation (18) one by one. By using Taylor expansion, we deduce
Therefore, from (19), it is immediate to realize that
Using Schwarz inequality we can verify
Then, we observe that (20), (24) and Lemma 2.1 yield
Combining (21) and (24) and Lemma 2.1, we have
Analogously, combining (22), (24) and Lemma 2.1, we deduce
Now we combine (23), (25), (26) and (27) to have
Multiplying equation (28) by and taking the sum of n, we derive
On the other hand, by Lemma 2.1, we have
Using triangle inequality and the above estimates gives the proof.
3. Numerical Experiment
In order to prove the validity of the theoretical analysis, we give the problem for which the exact solution is given. We apply the full-discrete scheme proposed above to solve this problem.
Example. Let the variables and T in problem (1) be equal to 1. Consider the following initial-boundary value problem:
The exact solution is .
Table 1. Errors and convergence order of Example.
Figure 1. Exact solution with .
Figure 2. Solution of fully discrete finite element schemes with .
We use the finite element fully discrete scheme to solve the problem. In Table 1, by choosing , we give the L2-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.
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 L2-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.
The work is supported by National Natural Science Foundation, China (11501335).
The authors would like to thank editor and referees for their valuable advice for the improvement of this article.
 Putrevu, H., Subramanian, H. and Mulay, S.S. (2021) On the Viscoelastic Dynamic Beam Modelling. International Journal of Advances in Engineering Sciences and Applied Mathematics, 13, 18-32. https://doi.org/10.1007/s12572-021-00288-8
 Huang, Z.C., Wang, X.G., Wu, N.X., Chu, F.L. and Luo, J. (2019) A Finite Element Model for the Vibration Analysis of Sandwich Beam with Frequency-Dependent Viscoelastic Material Core. Materials, 12, No. 20.
 Kouami, K., Foudil, M., Erasmo, C. and Mostafa, D.E. (2021) A Finite Element Approach for the Static and Vibration Analyses of Functionally Graded Material Viscoelastic Sandwich Beams with Nonlinear Material Behavior. Composite Structures, 274, Article ID: 114315.
 Snchez, E.D., Nallim, L.G., Bellomo, F.J. and Oller, S.H. (2019) Generalized Viscoelastic Model for Laminated Beams Using Hierarchical Finite Elements. Composite Structures, 235, Article ID: 111794.
 Lekdim, B. and Khemmoudj, A. (2021) Existence and Energy Decay of Solution to a Nonlinear Viscoelastic Two-Dimensional Beam with a Delay. Multidimensional Systems and Signal Processing, 32, 915-931.
 Wu, Q.L. and Qi, G.Y. (2020) Viscoelastic String-Beam Coupled Vibro-Impact System: Modeling and Dynamic Analysis. European Journal of Mechanics-A/Solids, 82, Article ID: 104012.
 Elhuni, H. and Basu, D. (2019) Dynamic Soil Structure Interaction Model for Beams on Viscoelastic Foundations Subjected to Oscillatory and Moving Loads. Computers and Geotechnics, 115, Article ID: 103157.
 Yu, C., Zhang, J., Chen, Y., Feng, Y.J. and Yang, A.M. (2019) A Numerical Method for Solving Fractional-Order Viscoelastic Euler-Bernoulli Beams. Chaos, Solitons & Fractals, 128, 275-279.