Symmetric Properties of Eigenvalues and Eigenfunctions of Uniform Beams

In this paper, the models of Euler–Bernoulli beams on the Winkler foundations are considered. The novelty of the research is in consideration of the models with an arbitrary variable coefficient of foundation. Qualitative results that influence the symmetry of the coefficient of foundation on the spectral properties of the corresponding problems are obtained, for which specific variable coefficients of foundation are tested using numerical calculations. Three types of fixing at the ends are studied: clamped-clamped, hinged-hinged and free-free. The conditions of the stiffness and types of beam fixing have been found for the set of eigenvalues of boundary value problems on a full segment and can be represented as two groups of the eigenvalues of certain problems on a half segment. Such qualitative spectral properties of a mechanical system can contribute to the creation of various algorithms for nondestructive testing, which are widely used in technical acoustics.


Introduction
In this paper, the eigenvalues of a uniform Euler-Bernoulli beam lying on a Winkler's type foundation at three types of fixings at the ends are studied: clamped-clamped, hinged-hinged and free-free [1]. The spectral properties of a mechanical system under any type of oscillations resting on elastic foundations have a wide range of applications in different science and engineering areas. These properties are used for modeling railway tracks, bridges, pavements, runways, pipelines, piles, etc. (see [1][2][3][4]).
Let k(x), x ∈ (0; 1) be a real-valued summable function, and λ is an eigenvalue. A force of foundation reaction is equal: where k(x) is the variable coefficient of foundation (kN m −1 ) and w(x, t) is the transverse displacement (m). The properties of the eigenvalues and eigenfunctions of boundary value problems for the Euler-Bernoulli beam with a constant coefficient of foundation were studied well. The comprehensive review on different theoretical elastic and viscoelastic foundation models in oscillatory systems can be found in [4]. The foundation of a mechanical system is generally described by a variable coefficient in numerous applied mathematical models. If the foundation coefficient is variable, then the question of studying the qualitative behaviors of the eigenvalues and eigenfunctions arises. Analytical solutions are not always obtained in an explicit form with a variable coefficient, as in the case of a constant coefficient.
Attempts to estimate the coefficients of foundation of a mechanical system were carried out in [5][6][7]. The approximate solutions of free vibrations in the integral form and Taylor series expansion for a uniform beam on a linear variable coefficient of foundation were considered in [2] (p. 109) and [8] (p. 26). Clastornik et al. [9] studied the approximation solution for finite beams with a polynomial coefficient k(x) = ∑ ∞ n=1 k n x n . The case of the variable coefficient of foundation where it is given as the random stiffness following form k(x) = k 0 + k ε (x), where k 0 = const, k ε (x) is the longitudinal variations of stiffness, was considered in [10][11][12].
Natural frequencies of the Euler-Bernoulli beam for the linear k(x) = k 0 (1 − αx), 0 ≤ α ≤ 1 and the nonlinear k(x) = k 0 (1 − βx 2 ), 0 ≤ β ≤ 1 functions were investigated in [13,14]. In [13], the general solution of the vibrations of uniform beams on variable Winkler elastic foundations were derived. In [14], the method of differential transformation was applied. In [3], the explicit form of the solution of the fourth order differential equation for k(x) = 1 (c 0 +c 1 x) 4 with hinged-hinged of fixing and at constant external load was obtained. In [15], the explicit solution static bending response of a finite uniform free-free Euler-Bernoulli beam resting on a Winkler elastic foundation in terms of special functions was derived. The boundary value problems for ordinary high-order differential operators in which the spectrum is absent or the spectrum is a countable set were considered in the case with constant coefficients in [16] and with variable coefficients in [17]. The odd order differential operators for the similar questions were investigated in [18,19]. The eigenvalues of the two-point boundary value problem of the ordinary differential operator for a triple differentiation with asymmetric weight were studied in [19]. In [20,21], these ones were studied for the second-order differential operators. The class of Volterra three-point problems for the Sturm-Liouville operator on the closed interval were described in [21]. Here, the Volterra was related to the double symmetry of the potential. The research methods used in this work develop the methods of works from [19,21].
It is important to develop the numerical methods together with the qualitative methods to study the spectral properties of the differential operators and allow one to test the correctness of quality results. Recently, based on the numerical methods for estimation of the eigenvalues, the variational iteration methods [22], the Haar wavelet method [23], the differential quadrature method and the Boubaker polynomials expansion scheme [24], the polynomial expansion and the integral technique methods [25], the Chebyshev spectral Method [26] was developed.
Basically, the effective methods for the numerical calculation of the eigenvalues of the fourth-order differential equation with definite variable coefficients were developed in [8][9][10][11][12][13][14][22][23][24][25][26]. However, they did not study the qualitative properties of the eigenvalues. In particular, these works did not consider the influence of the coefficients of the equation to the spectrum factorization. It is difficult to reveal such spectral properties of the system using numerical calculations. Here, it is necessary to use qualitative properties (symmetry).
In science, questions often arise to generalize well-known classical results for clearer understanding of the system properties and more detailed analysis. For instance, the first four alternating symmetric and asymmetric eigenfunctions with respect to the middle of uniform clamped-clamped Euler-Bernoulli beam without an elastic foundation (Winkler's coefficient is identically equal to zero) were defined in [27] (p. 270). An assumption was formulated that this tendency was preserved for all other eigenfunctions. A research problem arises, then: which nonzero coefficients of the Winkler foundation of this property are preserved? The novelty of this work is defining a condition for a nonzero variable coefficient of foundation under the symmetric properties of the eigenfunctions of the beam are preserved (see Theorem 1). Factorization of the set of eigenvalues was carried out. Similar results for other fixing types of beams were obtained (see Theorems 2 and 3).
The goal of this study is to find conditions for a variable coefficient of foundation and fixing types of a beam under which it is possible to establish the qualitative behavior of the eigenvalues and eigenfunctions. According to the authors' assumption, an important factor to achieve the goal is the symmetric properties of the variable coefficient of foundation and corresponding fixing types of beam.
The work is organized as follows. After the introduction, a mathematical model of free transverse vibration of the Euler-Bernoulli beam lying on a Winkler's type foundation are formulated. Eigenvalues and eigenfunctions with three types of fixings at the ends are investigated: clamped-clamped, hinged-hinged and free-free. In the third section, the examples are given to illustrate the results. Finally, main conclusions are outlined in closing.

The Problem of Transverse Vibrations of a Beam Lying on a Winkler's Type Foundation
Consider a beam of unit length endowed with a right-handed system of rectangular coordinates (x, y, z) with the origin at the center of gravity of the left cross-section of the beam; the x-axis is taken along the axis of the beam, while the yand z-axes are taken along the principal axis of cross section, the positive direction of the y-axis being oriented vertically downward (see Figure 1). The equation of free transverse vibrations of a beam of length = 1 has the following form: where ρ is the density of material, (kg m −3 ); A is the cross-sectional area, (m 2 ); E is the elastic modulus of material, (N m −2 ); J is the moment of inertia of the cross-sectional area of relative to around the z axis, (m 4 ). The problem of transverse vibrations of a beam of unit length after replacement w(x, t) = v(λ, x)sin(ωt) reduces to the following spectral problem: where v(λ, x) are the eigenfunctions of the transverse static deflection of the beam; λ = ρAω 2 are the eigenvalues, ω are the frequency parameters, (Hz); I p are the finite intervals. In this study, the bending stiffness, EJ, is assumed to be a constant; such a beam is called uniform. We present the necessary boundary conditions on finite intervals and their physical interpretations for convenience of presentation in Table 1 (see [27][28][29]). The main question of this work is to study the characteristics and the connections of the set of eigenvalues and eigenfunctions of the fourth order differential equation with a variable coefficient in Equation (1) with two-point boundary conditions in Equations (2), (5) and (8).

Clamped-Clamped Euler-Bernoulli Beam
First, consider the behavior of the eigenvalues of the clamped-clamped Euler-Bernoulli beam. Let σ(A 1 ), σ(B 1 ), σ(C 1 ) be a set of eigenvalues of problems A 1 − λI, B 1 − λI, C 1 − λI generated by Equation (1) on finite intervals by boundary conditions in Equations (2)-(4), respectively; the data are shown in Table 1. Table 1. The types of boundary conditions on the finite intervals and their interpretations.

Intervals
The Physical Interpretation The Types of the Boundary Conditions

Number of Conditions
Hinged-hinged at the points Main results. The first result of the work is Theorem 1. Let k(x) be a symmetric function with respect to the point x = 1 2 : The following statements are true: If λ ∈ σ(B 1 ) or λ ∈ σ(C 1 ), then the eigenfunctions of problems A − λI corresponding to the eigenvalues λ are symmetric or asymmetric with respect to the middle of the beam at the point x = 1 2 on the interval (0, 1), respectively.
An essential role is played by the following assertion.
Proof of Lemma 1. Consider the Cauchy problem with the zero initial conditions at the point . The initial conditions acquire the form From Equation (13), it follows that the first and third derivatives of any symmetric (α = 1) solution w + (λ, x) of Equation (1), as well as any asymmetric (α = −1) solution w − (λ, x) together with its second derivative, are zero at the point x = 1 2 . In particular, these properties hold for the fundamental system of solutions {v i (λ, x)}, (i = 1, 4), generated by the following conditions at the point Thus, when the condition of symmetry of Equation (11) is satisfied, there are two symmetric Note that the existence of a fundamental system of solutions at an arbitrary point x 0 ∈ [0, 1] for Equation (1) was proved in [30].
The proof is complete.
Proof of Theorem 1. We show correctness of the relation ) be a fundamental system of solutions generated by the conditions of Equation (14). Then, we get the characteristic determinant as [30] (pp. 1-27) Taking the relation in Equation (12) into account, the characteristic determinant ∆ A 1 (λ) can be rewritten in the following form: The characteristic determinants of the problems B 1 − λI and C 1 − λI are obtained by direct substitution and by taking Equation (14) into account in the following form: is the solution of Equation (1) and satisfies the condition in Equation (3): This function is the eigenfunction of problem B 1 − λI on the interval 1 2 , 1 . Moreover, its characteristic determinant is given by the formula The functionsṽ + (λ , x),ṽ 1 (λ , x) andṽ 3 (λ , x) will continue symmetrically on the interval 0, 1 2 if the condition of symmetry, Equation (11), is satisfied and if the Lemma 1 is taken into account. These are denoted as v + (λ , x), v 1 (λ , x) and v 3 (λ , x), respectively. Furthermore, we consider these functions on the entire interval (0, 1). The conditions at the point x = 1 are satisfied automatically for v + (λ , x).
Verification is required only for the conditions at the point Substituting v + (λ , x) into the condition in Equation (16), the characteristic determinant is obtained as: Taking the symmetry of the solutions v 1 (λ , x), v 3 (λ , x) into account with respect to x = 1 2 , this expression can be reduced to the form Since λ is a zero of ∆ B 1 (λ), from representation in Equation (15), it follows that it is a zero of ∆ A 1 (λ). Thus, the execution of the condition in Equation (16) is proved. Therefore, the continued function v + (λ , x) is the eigenfunction of the operator A 1 − λ I on the interval (0, 1). The second part of Statement 2 will be proved similarly.
The proof is complete.

Hinged-Hinged Euler-Bernoulli Beam
Consider the behavior of the eigenvalues of the hinged-hinged Euler-Bernoulli beam (see Figure 2). Let σ(A 2 ), σ(B 2 ), σ(C 2 ) be a set of eigenvalues of problems A 2 − λI, B 2 − λI, C 2 − λI generated by Equation (1) on finite intervals by boundary conditions in Equations (5)-(7), respectively; the data are shown in Table 1. Theorem 2. Let k(x) be a symmetric function with respect to the point x = 1 2 , i.e., the condition in Equation (11) holds. The following statements are true: If λ ∈ σ(B 2 ) or λ ∈ σ(C 2 ), then the eigenfunctions of problems A 2 − λI corresponding to the eigenvalues λ are symmetric or asymmetric with respect to the middle of the beam at the point x = 1 2 on the interval (0, 1), respectively.

Proof of Theorem 2. We show correctness of the relation ∆
Let {v i (λ, x)}, (i = 1, 4) be a fundamental system of solutions generated by the conditions in Equation (7). Then, we get the characteristic determinant as [30] (pp. : Taking the relation in Equation (12) into account, the characteristic determinant ∆ A 2 (λ) can be rewritten in the following form: The characteristic determinants of the problems B 2 − λI and C 2 − λI are obtained by direct substitution and by taking Equation (14) in the following form into account: is the solution of Equation (1) and satisfies the condition in Equation (6): This function is the eigenfunction of problem B 2 − λI on the interval 1 2 , 1 . Moreover, its characteristic determinant is given by the formula The functionsṽ + (λ , x),ṽ 1 (λ , x) andṽ 3 (λ , x) will continue symmetrically on the interval 0, 1 2 if the condition of symmetry in Equation (11) is satisfied and if the Lemma 1 is take into account. These are denoted as v + (λ , x), v 1 (λ , x) and v 3 (λ , x), respectively. Furthermore, we consider these functions on the entire interval (0, 1). The conditions at the point x = 1 are satisfied automatically for v + (λ , x). Verification is required only for the conditions at the point Substituting v + (λ , x) into the condition in Equation (18), the characteristic determinant is obtained as: Taking the symmetry of the solutions v 1 (λ , x), v 3 (λ , x) into account with respect to x = 1 2 , this expression can be reduced to the form Since λ is a zero of ∆ B 2 (λ), from representation in Equation (17), it follows that it is a zero of ∆ A 2 (λ). Thus, the execution of the condition in Equation (18) is proved. Therefore, the continued function v + (λ , x) is the eigenfunction of the operator A 2 − λ I on the interval (0, 1). The second part of Statement 2 will be proved similarly.
The proof is complete.

Free-Free Euler-Bernoulli Beam
Consider the behavior of the eigenvalues of the free-free Euler-Bernoulli beam (see Figure 3). Let σ(A 3 ), σ(B 3 ), σ(C 3 ) be a set of eigenvalues of problems A 3 − λI, B 3 − λI, C 3 − λI generated by Equation (1) on finite intervals by the boundary conditions in Equations (8)-(10), respectively; the data are shown in Table 1. Theorem 3. Let k(x) be a symmetric function with respect to the point x = 1 2 , i.e., the condition in Equation (11) holds. The following statements are true: If λ ∈ σ(B 3 ) or λ ∈ σ(C 3 ), then the eigenfunctions of problems A 3 − λI corresponding to the eigenvalues λ are symmetric or asymmetric with respect to the middle of the beam at the point x = 1 2 on the interval (0, 1) respectively.

Proof of Theorem 3. We show correctness of the relation ∆
Let {v i (λ, x)}, (i = 1, 4) be a fundamental system of solutions generated by the conditions in Equation (14). Then, we get the characteristic determinant as [30] (pp. : Taking the relation in Equation (12) into account, the characteristic determinant ∆ A 3 (λ) can be rewritten in the following form: The characteristic determinants of the problems B 3 − λI and C 3 − λI are obtained by direct substitution and by taking Equation (14) in the following form into account: , the proof of Statement 1 of Theorem 3 follows. (1) has the following form:

Proof of the Statement 2. The general solution of Equation
for the interval 1 2 , 1 . Here, {ṽ i (λ, x)}, (i = 1, 4) are the fundamental system of solutions of Equation (1) given by the conditions in Equation (14) on the interval 1 2 , 1 . Let λ ∈ σ B 3 (λ) be the eigenvalues; then, a linear combination of the fundamental solutions is the solution of Equation (1) and satisfies the condition in Equation (9): This function is the eigenfunction of problem B 3 − λI on the interval 1 2 , 1 . Moreover, its characteristic determinant is given by the formula The functionsṽ + (λ , x),ṽ 1 (λ , x) andṽ 3 (λ , x) will continue symmetrically on the interval 0, 1 2 if the condition of symmetry in Equation (11) is satisfied and if Lemma 1 is taken into account. These are denoted as v + (λ , x), v 1 (λ , x) and v 3 (λ , x), respectively. Furthermore, we consider these functions on the entire interval (0, 1). The conditions at the point x = 1 are satisfied automatically for v + (λ , x). Verification is required only for the conditions at the point Substituting v + (λ , x) into the condition in Equation (20), the characteristic determinant is obtained as: Taking the symmetry of the solutions v 1 (λ , x), v 3 (λ , x) into account with respect to x = 1 2 , this expression can be reduced to the form Since λ is a zero of ∆ B 3 (λ), from the representation in Equation (19), it follows that it is a zero of ∆ A 3 (λ). Thus, the execution of the condition in Equation (20) is proved. Therefore, the continued function v + (λ , x) is the eigenfunction of the operator A 3 − λ I on the interval (0, 1). The second part of Statement 2 will be proved similarly.
The proof is complete.

Examples
In this section, we calculate approximately the five eigenvalues of boundary value problems The numerical results of first five eigenvalues for Example 1 are shown in Table 2. Example 2. Let k(x) = 4x(1 − x) and EJ = 1.
The numerical results of first five eigenvalues for Example 2 are shown in Table 3. The calculations which represent in Examples 1 and 2 confirm the validity of Statement 1 of the Theorem 1 on the factorization of the set of eigenvalues. The physical interpretation of Statement 2 was given in [27] (p. 270) for k(x) ≡ 0. Example 3. Let k(x) = 4x(1 − x) and EJ = 1. The function k(x) satisfies the symmetry condition in Equation (11).
The violation of the regularity of factorization of eigenvalues in Example 4 is due to the asymmetry of the function k(x). The above calculations in Examples 3 and 4 confirm the validity of Statement 1 of Theorem 2. The numerical results of the first five eigenvalues for Examples 3 and 4 are shown in Tables 4 and 5, respectively. Example 5. Let k(x) = 10 − 40x(1 − x) and EJ = 1. The function k(x) satisfies the symmetry condition in Equation (11).
The numerical results of the first five eigenvalues for Example 5 are shown in Table 6. The calculations represented in Example 5 confirm the validity of Statement 1 of the Theorem 3. The eigenfunctions that confirm Statement 2 were given in [22] (Example 3.2) for k(x) ≡ 0. According to the authors, the qualitative behavior of the spectrum is revealed on the basis of the qualitative theory of differential equations. In this work, numerical methods give the possibility to test hypotheses and are used as an additional tool. We used a numerical method for calculating the eigenvalues at variable coefficient of k(x) with the polynomial expansion and integral techniques [25].
The degree of the polynomial was chosen N = 25, which ensures high accuracy of calculations. Numerical calculations were carried out in the Maple computer mathematics system.
The new results can be useful for solving the inverse coefficient problems using known natural frequencies (see [31]). In the future, the practical interest will be to investigate the behavior of a beam (thermodynamic response) with a variable foundation coefficient when the beam is irradiated by moving laser pulses. This issue was investigated in [32] for the constant coefficient. The results can be generalized to investigate the dynamic design of thin films on compliant substrates [33].

Conclusions
In this paper, the problems for determining the eigenvalues of the uniform Euler-Bernoulli beam lying on the Winkler's type foundation at three types of fixings at the ends were solved: clamped-clamped, hinged-hinged and free-free. The types of boundary value problems for a beam with a symmetric coefficient of foundation are described, and the set of eigenvalues on a full segment can be represented as two groups of the eigenvalues of certain problems on a half segment. If λ belongs to the first group, the eigenfunctions of the problems on the complete segment corresponding to the eigenvalues λ are symmetric with respect to the middle of the beam of the point x = 1 2 if the second ones are asymmetric. Numerical calculations are carried out to check the validity of Statement 1 of Theorems 1-3. The proposed method can be useful for solving inverse problems of restoring the elasticity coefficient of foundation from natural frequencies and developing algorithms for nondestructive testing devices.