## 1 Introduction

The equation of motion of a thin beam with weak-internal damping undergoing cylindrical bending can be written as

(1) |

where is the transverse displacement,

is the vector of spatial coordinates,

is the time. The aerodynamic damping term is denoted by is the nonlinear stiffness of beam, is an in-plane tensile load, and belongs to the noncylindrical domain, defined by(2) |

with is an open, limited and regular and is a - real function. All quantities are physically nondimensionalized, are fixed positive and is not necessarily positive.

In the absence of the damping term , many authors have already studied Cauchy’s problems and the mixed problems associated with the equation (1) on bounded and unlimited cylinders. The existence and uniqueness of solution was demonstrated in [9].

The Cauchy problem associated with the (1) equation in abstract terms on a Hilbert space is studied, among other authors, by [4, 18]. In addition, in the last were established results on existence, unicity and asymptotic stability of the solution.

For the one-dimensional case, omitting the term and considering the cylinder , with open and limited, we obtain the Kirchhoff equation with internal damping, namely,

(4) |

The equation (4) has been extensively studied by a wide variety of authors on -dimensional cases and generic mathematical models defined in Hilbert spaces. The existence of a local and global solution can be found in several physical-mathematical contexts, such as [14, 15].

In the article [15] a new model associated with the Kirchhoff equation (4) was introduced, called the Medeiros-Kirchhoff equation, given by

There are several areas that work with the application of the Beam equation, for example, in civil and naval engineering, as well as in the aerospace area as the reaction of rocket and satellite metal structures to various space situations, as can be observed in the articles [12, 13]. In addition to these areas, it is possible to associate Maxwell’s equations, which govern electromagnetic theory, and the Euler-Bernoulli Beam Equation, in the study presented by article [16]

The equation (1) follows the Euler-Bernoulli Beam model. In addition to this model, there are other models, extensions of the Euler-Bernoulli model, which consider rotational movements besides transverse ones, being possible to create models closer to reality. The article [11] presents a more detailed study on three of these models.

The existence and uniqueness of the solution, the asymptotic decay and the numerical simulations of the beam equation with movinf boundary were studied in [Raquel] for the one-dimensional case.

In this paper, we present a family of numerical method for the one-dimensional and two-dimensional cases, based on the finite element method and the finite difference method. In addition, we are doing numerical analysis for both the semi-discrete problem and the fully discrete problem, showing that the order of convergence is quadratic in space and time. For this we will use the Faedo-Galerkin method, using the Hermite polynomials as base functions. In addition, the nonlinear system associated with the system of ordinary differential equations is being solved by Newton’s method and for temporal discretization is applying the Neumark’s method, where we will show the convergence of the approximate solution to a family of numerical methods depending on the choice of

. Numerical simulations are presented as examples, with different types of boundary and tables showing the order of numerical convergence. In addition, we present the asymptotic decay, and the graphs of the solutions for the one-dimensional and two-dimensional cases.## 2 Analytical results

Our section will present, without demonstrating, the existence and uniqueness results of the mixed problem solution for the equation (1), considering the boundary conditions (3), given by:

Let us consider the diffeomorphism , defined by , with . Therefore, the change of variable transforms the problem(5) into the equivalent problem with cylindrical boundary

(6) |

where, for simplicity, we denote and

(7) | ||||

In the following, we present the results on the existence and uniqueness of the weak solution of problems (5) and (6) and asymptotic decay, whose demonstrations can be found in [8] for the one-dimensional case and [5] for the two-dimensional case.

###### Theorem 1.

(Existence and Uniqueness of solution in nocylindrical domain.)

###### Theorem 2.

(Existence and Uniqueness of solution in cylindrical domain)

###### Theorem 3.

(Asymptotic Decay: Two-dimensional case)

From the important theoretical results, we can now develop numerical methods to determine an approximate numerical solution and verify the consistency of the numerical methods used.

Due to the equivalence of the problems, the methods used will be for the problem (6) that is defined in a cylindrical domain.

## 3 Numerical Method

The purpose of this section is to develop a numerical method, based on the finite element method and the finite difference method to make the numerical simulations in order to obtain a approximate solution of the problem (6).

Denoting the Hilbert space , then weak formulation of the problem (6) can be written in the follows form, when we integrate by parts the second, third, fifth and sixth terms:

(8) |

where .

### 3.1 Approximated Problem

Let be a subspace generated by the first vectors of the basis of Hilbert space . Then we want to determine an approximate solution of the following approximate system, given by:

(9) | ||||

where we integrate by parts the sixth term.

Since that, , then the function can be represented by

(10) |

Taking , for some , and substituting (10) in (9), we obtain

(11) | ||||

In order to validate the method (11), we add the term source on the right side of the equation (9) and, taking , we have the following definitions for the matrices and vectors:

(12) | ||||

Considering the definitions (12) and varying , we obtain the following nonlinear system of second order ordinary differential equations:

(13) |

with

(14) | ||||

where we denote by the transpose of the matrix .

Now, to solve the system of ordinary nonlinear equations (13) in discrete time, we will use the finite difference method.

### 3.2 The Newmark’s method

In order, consider the discretize the fixed time interval , using a uniform time step of size , then , for and . Consider the Neumark’s approximations:

Substituting the approximations in the system (13) and taking the equation in the discrete time , we obtain

(15) | ||||

or equivalently

(16) |

where

#### Remark.

For the initial step, that is, for , appears in the equation the term , that does not exist. Proceeding formally, we use a quadratic order approximation for the initial data . Therefore, we consider , so we have . So, for in (16), we have

(17) |

where we assume the approximation and

considering the approximation .

### 3.3 Newton’s method

For each , consider the left side of the equation in (17). In this way, solve (17) is equivalent to finding the root of . By Newton’s method, given the initial approximation and denoting ,, we obtain a new approximate solution of , given by , where is the linear system solution . The key point of the method is the calculation of the Jacobian matrix . This is because, if it has a high computational cost, it can make the method unfeasible.

Define

(18) |

where .

For , follows

(19) |

where .

For :

(20) | ||||

and for :

(21) | ||||

where

Using (14), we obtain

where is the -th coordinate of the initial velocity vector .

From there, we apply Newton’s refinement algorithm, for , until one of the stopping criteria is satisfied: or , where is the solution obtained in the -th iteration of Newton. We also define a maximum limit for iterations in order to avoid infinite refinements.

### 3.4 Finite Element Method- Hermite Base

Let be a polygonalization family of , where , satisfying the minimum angle condition (see [7]) and indexed by the parameter , representing the maximum diameter of the elements . Given an integer , we introduce the following finite element space

spanned by the Hermite interpolation polynomials basis functions. More specifically, we let

where represents a finite-dimensional function space generated by piecewise polynomial functions of degree less or equal associated with a discretization by finite elements of size . In particular, in this paper we will use the Hermite polynomials of degree equal .

From the Lemma of Douglas-Dupont (see in [7]), it follows that, given a function , there is an interpolator such that

where denotes the seminorm over Hilbert space and .

Knowing the Hermite polynomials, then all the matrices of the system (16) can be calculated explicitly and solving the nonlinear system by Newton’s method gives the approximate solution, given by

We note that since the Hermite polynomials are orthogonal then

## 4 Numerical analysis

In this section we present error estimates for the finite element semi-discrete and fully discrete approximations of the nonlinear beam equation (8) or equivalently for the weak formulation of the nonlinear beam equation with moving ends.

Remarks on regularity: In order to get an error estimate, we require more regularity for the function and its derivatives in time:

(H2) |

### 4.1 Error estimates for the semi-discrete problem

The semi-discrete problem consists of finding such that satisfies the problem (8) restricted to , that is, to find for all such that

(22) |

For the estimates, we will need the following hypothesis about the projections of in , denoted by and defined in (27), of the initial data

(H3) |

where and are the projections of and in , respectively, and are constants independent of .

Let us assume that,

(H4) |

which requires that the speed of the end points be smaller than the characteristic speed of the equation.

Let us denote by , the error between the exact and approximate solution. Then, we have the following error estimate:

###### Theorem 4.

(Semi-discrete problem - Estimate)

where the constant is dependent on and independent of .

Comments

There are no comments yet.