Journal of Applied Mathematics

Volume 2014, Article ID 298092, 7 pages

http://dx.doi.org/10.1155/2014/298092

## Solvability of a Model for the Vibration of a Beam with a Damping Tip Body

Department of Mathematics and Applied Mathematics, University of Pretoria, Pretoria 0002, South Africa

Received 22 January 2014; Revised 11 April 2014; Accepted 11 April 2014; Published 15 May 2014

Academic Editor: Francisco J. Marcellán

Copyright © 2014 M. Basson et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

We consider a model for the vibration of a beam with a damping tip body that appeared in a previous article. In this paper we derive a variational form for the motion of the beam and use it to prove that the model problem has a unique solution. The proofs are based on existence results for a general linear vibration model problem, in variational form. Finite element approximation of the solution is discussed briefly.

#### 1. Introduction

In article [1], the authors model and analyze the damped vibration of a cantilever beam with an attached hollow tip body that contains a granular material. The Euler-Bernoulli theory for a beam with Kelvin-Voigt damping is used. The beam is clamped at one end and the tip body is attached to the other end. The authors state that “the problem contains more complicated boundary conditions” than problems “treated previously” (and provide references). This is due to the fact that the model is more realistic—as we explain below. It is an interesting model for more than one reason. The dynamics of the rigid body is treated in a realistic way: the fact that the center of mass of the rigid body is not at the endpoint of the beam is taken into account. The damping mechanism is also explained unlike other papers where vibration models with boundary damping are considered. There are other articles with realistic models of a beam with damping, for example [2]. However, [1] provides the only realistic model for the relevant beam-body configuration.

The existence of a unique solution for the model problem is established in [1]. To obtain the result, the problem is written as an abstract differential equation and an abstract existence result from a previous paper [3] is applied.

In this paper, we also prove the existence of a unique solution. Our approach differs from that in [1]: we write the model problem in variational form and use results from [4] where general linear vibration models in variational form are considered. The existence theorems in [4] were applied to the vibration of a complex plate beam system in [5]. It should be noted that the same variational form can be used for finite element approximations (see Section 7). Since our approach differs from that in [1], it is natural to consider differences regarding assumptions and results. There are indeed some differences but these are not substantial (see Section 6).

The mathematical model is considered in Section 2. In Section 3, we write the model problem in variational form and present the weak variational form in Section 4. Auxiliary results are proved in Section 5. The existence theorems are stated and proved in Section 6 where different methods are also compared. In Section 7, natural frequencies and modes are discussed as well as finite element approximation.

#### 2. Model Problem

The Euler-Bernoulli model for the transverse vibration of a beam is derived from the equation of motion for the beam deflection : and the relation

In these equations, denotes the density, the area of the cross section, the shear force, the bending moment, and a load (beam models are treated in [6, pages 323-324], [7, pages 337-338], and [8, pages 392–395]).

The usual constitutive equation is , where is an elastic constant (Young’s modulus) and is the area moment of inertia. Due to Kelvin-Voigt damping, it changes to where denotes the damping parameter. The partial differential equation (which we do not use) is

The constitutive equation for the moment and the relation between the moment and shear force are also used to model the interface conditions.

The left endpoint of the beam is clamped where the boundary conditions are the usual

The interface conditions at the other endpoint are determined by the interaction between the beam and the rigid body. This is explained in [1] in some detail. It is necessary to consider the equations of motion for the rigid body carefully when deriving these conditions.

The position of the center of mass of the tip body relative to the endpoint of the beam is where is the angle of the neutral plane with the horizontal (see Figure 1). Therefore, the velocity and acceleration of the center of mass are given by where denotes the length of the beam. For the linear approximation, it is assumed that the term may be neglected, , , and . Using these approximations, we have the following expressions for the vertical components of the velocity and acceleration:

In [1], the term in the expressions for the vertical component of the velocity is neglected. In our opinion, this should not be done and we motivate our point of view in the next section where we discuss the decay of energy for the system.

In the equations below, and denote damping parameters, the mass, and the moment of inertia of the rigid body. Using Newton’s second law for the motion of the center of mass, we have where is an external force that may act on the rigid body, for example, gravity. Taking moments about the center of mass, we have

Following [1], we combine (9) and (10) and find that

It is convenient to rewrite (9) and (11) as follows:

*Model Problem*. The mathematical model consists of equations of motion (1) and (2) and constitutive equation (3) for the beam, boundary conditions (5), and interface conditions (12). Initial conditions and need to be specified.

#### 3. Variational Form

Multiply (1) by an arbitrary smooth function and integrate. Using integration by parts and (2) yields

*Test Functions*. A function is a test function if , is integrable, and . The space of test functions is denoted by .

It follows that for each .

We now use the constitutive equation and the interface conditions (12) to derive the variational form of the model problem.

It is convenient to introduce the following bilinear forms:

We now have the variational form of the model problem.

*Problem PV*. Find such that, for each , and
for each , with and .

*Remark 1. *Note that the bilinear forms , , and are symmetric. The additional term in the definition of is necessary for symmetry. Problem may be used to compute finite element approximations.

*Mechanical Energy*. The mechanical energy (kinetic energy plus elastic potential energy) of the system is

Using the symmetry of and and assuming that is sufficiently smooth, we have for the homogeneous case. It is obvious that is nonnegative and not difficult to show that and are nonnegative:

As a result, . This result is to be expected from Physics. The fact that is symmetric and is nonnegative is due to the additional term.

#### 4. Weak Variational Form

Let denote the Sobolev space with weak derivatives up to order in . The inner product for is denoted by , with , the inner product for . The corresponding norms are and .

*Consider Problem PV*. We start off as usual by considering the closure of the space of test functions. Let be the closure of in ; then is a Hilbert space (being a closed subspace of a Hilbert space).

We require the so-called* trace operator* which we denote by . For , but (as is well known) can be extended to ; see, for example, [9]. Clearly is defined for .

The following product spaces are necessary for the abstract problem:

An element is written as . An obvious inner product for is and we denote the corresponding norm by .

*Definition 2 (bilinear forms). *For and in ,
and for and in ,

*Remarks. * The bilinear forms , , and are symmetric. For and in , , , and . It is essential to use product spaces since the bilinear form is not defined on .

For the weak variational form of the model problem, we need to show that the bilinear forms and are inner products for and , respectively. We use the well-known Poincare type inequalities given in the proposition below. The boundedness of is also required.

Proposition 3. *For each in ,
*

*Proof. *Using the fundamental theorem of calculus, the inequalities are easy to prove for the space of test functions . Since is the closure of with respect to the norm of , the result follows.

*Proposition 4. There exists a constant such that
*

*Proof. *It is sufficient to show that there exists a constant such that for each . From a well-known inequality,

By the definition of the moment of inertia of a rigid body, for some . Now choose and and the desired result follows.

*Corollary 5. The bilinear form is an inner product for the space .*

*Definition 6 (inertia space). *The norm is defined by . We refer to the vector space equipped with this norm as the inertia space and denote it by .

*Proposition 7. There exists a constant such that
*

*Proof. *We use Proposition 3 and the definition of the bilinear form :

*Corollary 8. The bilinear form is an inner product for .*

*Proof. *Clearly implies that and therefore and .

*Definition 9 (energy space). *The space equipped with the inner product is referred to as the energy space. The norm is defined by .

*We proceed to determine a weak variational form of the model problem. Let , , and , with , and arbitrary.*

*Problem PW*. Find such that, for each , , , , and
with and .

*Remark 10. *It is natural to think that and so forth are the correct initial conditions. This is discussed in Section 6.

*5. Auxiliary Results*

*5. Auxiliary Results*

*We need the results of this section to apply Theorems 15 and 16 in Section 6.*

*Proposition 11. Space is a dense subset of .*

*Proof. *Consider any . Since is dense in , there exists a sequence such that .

It is not difficult to construct sequences and in with the following properties:

Now, let ; then , , and .

Consequently, = and .

*Proposition 12. There exists a constant such that
*

*Proof. *We use Proposition 3.

Consider

Now apply Proposition 7.

*Proposition 13. There exists a constant such that
for each and in .*

*Proof. *We use Proposition 3.

Consider

Now use Proposition 7.

*The result above is true for . If , the bilinear form is positive definite on and this has implications for existence results.*

*Proposition 14. Consider
*

*Proof. *We have that .

*6. Existence*

*6. Existence**In this section, we apply the existence results from [4]. For convenience, we formulate the general linear vibration problem and present the relevant existence theorems. Let , , and be real Hilbert spaces with . Spaces , , and have inner products , , and and norms , , and , respectively. Consider also a bilinear form defined on .*

*Problem PG*. Find a function such that, for each , , , , and
with and .

*In Theorems 15 and 16, the following is assumed.*

*Assumptions.* (A1) is dense in and is dense in . (A2) There exists a constant such that for each . (A3) There exists a constant such that for each . (A4) The bilinear form is symmetric, nonnegative, and bounded on ; that is, for each and in .

*Theorem 15 (see [4, Theorem 1]). Suppose that assumptions (A1), (A2), (A3), and (A4) are satisfied. If(a),(b), and there exists a such that
then problem PG has a unique solution:
*

*Theorem 16 (see [4, Theorem 3]). Suppose that assumptions (A1), (A2), (A3), and (A4) are satisfied. If(a)the bilinear form is positive definite, that is, there exists a positive constant such that , for each ,(b) is locally Hölder continuous on ,(c), ,*

then problem PG has a unique solution:

If , then .

*6.1. Applying General Results*

*6.1. Applying General Results*

*Theorem 15 above is applied to the case where and Theorem 16 to the case where . Note that assumptions (A1), (A2), (A3), and (A4) are satisfied due to Propositions 11, 12, 4, and 13, respectively. In the formulation of the theorems, we denote the function by .*

*Theorem 17. Suppose that and (a) and ,(b), and there exists a such that
*

Then problem has a unique solution:

*Proof. *Clearly . The result follows from Theorem 15.

*Theorem 18. Suppose that and(a) is locally Hölder continuous on with respect to the norm of and is locally Hölder continuous on ,(b) and .*

Then problem has a unique solution:

If , then .

*Proof. *Clearly is locally Hölder continuous on with respect to the norm and the bilinear form is positive definite by Proposition 14. The result follows from Theorem 16.

*6.2. Sufficient Conditions for Existence*

*6.2. Sufficient Conditions for Existence*

*The case is trivial.*

*If , sufficient conditions on and are required to satisfy condition (b) in Theorem 17. It is obviously necessary to assume that and are in which implies that and . Suppose that and . From the definition of the bilinear form , using integration by parts, we obtain
*

*From the definition of , we have
*

*Therefore,
if and only if
for each . Since and are arbitrary, it follows that
*

*Therefore, a sufficient condition for existence is , , and
*

*The conditions for the shear force and bending moment have an interesting physical interpretation. Comparing them to (9) and (11), we see that the force and moment at the endpoint must match the force and moment due to damping.*

*6.3. Discussion*

*6.3. Discussion*

*In [1], the model problem is written in the form
equation (3.11). The existence of a unique solution then follows from the theory in [3]. It should also be possible to use [10]. Relevant abstract existence results may also be found in other publications, for example, [11].*

*Theorems 15 and 16 (existence results in [4]) are convenient for application when the model problem is in weak variational form. This is so because the assumptions are in terms of the bilinear forms , , and and it is not necessary to consider linear operators as in the publications cited above.*

*The approach in [4] is relatively new and therefore we discuss briefly how it is related to semigroup theory. Problem G is equivalent to a first order differential equation in the product space . A linear operator with domain is constructed in [4] and problem G is equivalent to an initial value problem of the form
*

*The pair if and only if and in satisfy condition () in Theorem 15. The properties of are determined by the properties of the three bilinear forms , , and . Under the assumptions in Theorem 15 is the infinitesimal generator of a semigroup of contractions and under the assumptions in Theorem 16 is the infinitesimal generator of an analytic semigroup.*

*7. Application*

*7. Application*

*7.1. Natural Frequencies*

*7.1. Natural Frequencies*

*In the second half of Section 4 in [1], the sequence of natural frequencies of the undamped system is considered. The conjecture on p1041 concerning the eigenvalues for the undamped system is indeed correct. For each , the problem , for each , has a unique solution . The mapping , defined by , is a symmetric linear operator. Since is bounded as a mapping from into and the embedding of into is compact, is compact. Considering as a bilinear form in , its eigenvalues are real and if is an eigenvalue, then is an eigenvalue of . The corresponding eigenvectors are the same. It follows that the sequence of eigenvalues tends to infinity and the sequence of eigenvectors is complete in (see, e.g., [12, Theorem 4.A, p.232]).*

*For the model problem with the damping tip body, the situation is different. To the best of our knowledge, there is no general spectral theory for systems with boundary damping. However, results are known for specific model problems. In [13], it is proved that the sequence of eigenfunctions for an Euler-Bernoulli beam with boundary damping has the Riesz basis property, but there is no attached body.*

*Galerkin Approximation for Problem *. Let be a finite dimensional subspace of .

*Problem PW*^{h}. Find such that, for each , and
with and . The functions and are suitable approximations for and in .

*Problem PW ^{h} is equivalent to a system of ordinary differential equations:
*

*The system can be used to approximate the solution of problem PW. How to construct the relevant matrices is explained in [14].*

*The quadratic eigenvalue problem
can be used to calculate approximations for the natural frequencies (see [14]).*

*Conflict of Interests*

*Conflict of Interests*

*The authors declare that there is no conflict of interests regarding the publication of this paper.*

*References*

*References*

- K. T. Andrews and M. Shillor, “Vibrations of a beam with a damping tip body,”
*Mathematical and Computer Modelling*, vol. 35, no. 9-10, pp. 1033–1042, 2002. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - J. W. Hijmissen and W. T. van Horssen, “On aspects of damping for a vertical beam with a tuned mass damper at the top,”
*Nonlinear Dynamics*, vol. 50, no. 1-2, pp. 169–190, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - K. L. Kuttler Jr., “Time-dependent implicit evolution equations,”
*Nonlinear Analysis: Theory, Methods & Applications*, vol. 10, no. 5, pp. 447–463, 1986. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - N. F. J. van Rensburg and A. J. van der Merwe, “Analysis of the solvability of linear vibration models,”
*Applicable Analysis*, vol. 81, no. 5, pp. 1143–1159, 2002. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - N. F. J. van Rensburg, L. Zietsman, and A. J. van der Merwe, “Solvability of a Reissner-Mindlin-Timoshenko plate-beam vibration model,”
*IMA Journal of Applied Mathematics*, vol. 74, no. 1, pp. 149–162, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - Y. C. Fung,
*Foundations of Solid Mechanics*, Prentice-Hall, Englewood Cliffs, NJ, USA, 1965. - D. J. Inman,
*Engineering Vibration*, Prentice-Hall, Englewood Cliffs, NJ, USA, 1994. - D. E. Newland,
*Mechanical Vibration Analysis and Computation*, Longman, Essex, UK, 1989. View at MathSciNet - J. T. Oden and J. N. Reddy,
*An Introduction to the Mathematical Theory of Finite Elements*, John Wiley & Sons, New York, NY, USA, 1976. View at MathSciNet - K. T. Andrews, K. L. Kuttler, and M. Shillor, “Second order evolution equations with dynamic boundary conditions,”
*Journal of Mathematical Analysis and Applications*, vol. 197, no. 3, pp. 781–795, 1996. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - R. E. Showalter,
*Hilbert Space Methods for Partial Differential Equations*, Pitman, London, UK, 1977. View at MathSciNet - E. Zeidler,
*Applied Functional Analysis: Applications to Mathematical Physics*, Springer, New York, NY, USA, 1995. View at MathSciNet - B.-Z. Guo and R. Yu, “The Riesz basis property of discrete operators and application to a Euler-Bernoulli beam equation with boundary linear feedback control,”
*IMA Journal of Mathematical Control and Information*, vol. 18, no. 2, pp. 241–251, 2001. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - A. Labuschagne, N. F. J. van Rensburg, and A. J. van der Merwe, “Distributed parameter models for a vertical slender structure on a resilient seating,”
*Mathematical and Computer Modelling*, vol. 41, no. 8-9, pp. 1021–1033, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus

*
*