1. Introduction
The model of the diatomic chain with perfect regularity, which assumes a simple cubic lattice with two different masses and the nearest neighbor interaction only, predicts that its frequency spectrum consists of two parts: The upper optical band and the lower acoustic band [
1,
2]. The effect of the introduction of an impurity atom into a perfect crystal is to bring a local defect, which changes the mass and stiffness at that given point. The existence of a local defect generates two new eigenfrequencies, which are isolated from both the optical and acoustic bands [
1,
2]. The normal modes associated with these two new eigenfrequencies are the localized modes [
1], in which vibration concentrates around the local defect. The mechanism for the localized modes can be roughly given as follows: When a crystal vibrates with these isolated eigenfrequencies, which are the forbidden bands for most of the atoms, those atoms keep motionless except that of the local defect and few neighboring ones vibrating with much smaller amplitudes [
1]. As early as in 1877, Lord Rayleigh’s study on the effect of the parameter variations on dynamical systems showed that under right circumstances a vibration mode can split away from the band of modes and exists above the maximum frequency of the unperturbed system [
3]. In optics, the localized mode was firstly observed by Schaefer [
4] in his experiment on the absorption spectrum of potassium chloride (KCl) at the frequencies above the maximum phonon frequency, in which a sharp localized optic mode appears due to the presence of a hydrogen ion impurity. In the study of electron transport, Anderson’s model [
5] (nowadays known as the Anderson localization) shows that the electron motions becomes localized once an impurity is introduced into an “unperturbed” system, which can lead to the transition of a metallic conductor to a semiconductor or an insulator. In acoustics, the localized modes are realized by the so-called locally resonant sonic materials, which consist of two different materials: hard lead particle coated with soft silicone rubber [
6]. The transmission coefficient of the lead–rubber composite structure has the sudden decreases at two different frequencies of an incidence wave, which are referred to as the dip frequencies [
6]. At the lower dip frequency, the lead particle is similar to a mass in an oscillator, which moves as a whole as a rigid body; at the same time, the rubber displacements except those around the interface are very small. At the higher dip frequency, things are reversed: the lead particle now has almost no motion and most of the vibration is in the rubber. The motions at these two dip frequencies are referred to as the localized modes [
6]. In these localized modes, most vibrations are confined either in the rubber layer or as the oscillation of the lead particle [
6]; almost no waves are transmitted into the lead particle, i.e. the particle has no deformation. This localized mode provides another effective way of sound attenuation rather than absorption [
6]. In optomechanics, a strong optomechanical coupling results in a localized mode, in which the reflection/absorption of photons or microwave phonons is dramatically changed [
7]. Similar to in solid mechanics, the mode localization is also encountered once the periodicity/cyclicity of a structure is broken by a defect/irregularity [
8,
9].
Mode localization and veering are closely related [
10]. Some scholars think that they are the two manifestations of a same drastic phenomenon [
11]. Veering describes the phenomenon that two eigenvalues, which can physically correspond to either eigenfrequencies [
10,
11] or buckling loads [
12,
13], rapidly approach each other and then diverge abruptly without crossing. To emphasize this non-crossing property, veering is also variously called avoided crossing [
14] or anti-crossing [
15]. Coupling is the mechanism responsible for the mode localization and veering [
16]. When two eigenfrequencies are well-separated, the effect of mode coupling is usually modest, and the modal motions are essentially independent from each other [
17]. The mode coupling can be enhanced when the ratio of two eigenfrequencies is an integer, which is also called harmonic [
18]. However, most structures with small linear vibrations do not have the harmonic property; various methods, such as the step-like design [
19], are used to make a cantilever beam harmonic. Mode coupling can also be enhanced in nonlinear regime, in which large motion induced tension can effectively change the system stiffness and, therefore, tune the system eigenfrequencies to be harmonic [
17,
20]. Another alternative approach to enhance the mode coupling is to bring the two eigenfrequencies closer because veering indicates strong mode coupling [
15]. For example, two modes of a square membrane with the eigenfrequency difference of 0.07% have been used to achieve a strong mode coupling in a cavity optomechanical system [
21]. In the Kuramoto model of coupled oscillators, the effective strength of coupling is the paramount parameter of determining the synchronization [
22] and pattern formation [
23] in a dynamical system. However, extracting the effective strength of coupling directly from a real physical system is rather difficult. An advantage of studying the eigenfrequency curve veering is to provide an efficient way of evaluating/obtaining the effective strength of coupling.
The coupling between the optical and mechanical modes can result in either the damping or anti-damping force for either mode [
24]. The damping effect has been used to cool down the thermal motion of a macroscopic resonator to the quantum ground state [
25]. Quality factor is an essential factor determining the mass sensitivity in a mass resonator sensor application [
26,
27,
28]. The quality factor of a mechanical resonator can be significantly improved by the anti-damping due to the optomechanical coupling [
29]; the cavity optomechanical system is a possible candidate for the ultimate limit of mass detection [
7,
30]. One of the fast developing trends in optomechanics is to replace the optical cavity with an electrical one, which leads to electromechanics. When the materials of the mechanical end mirror are highly transparent to light, the radiation pressure exerted by photons can be weak, which is the means of the optomechanical coupling [
31]. For example, graphene is 98% transparent to light, but it is basically opaque to the electric field [
31]. Therefore, instead of an optical cavity, a microwave-regime electric cavity is introduced [
7,
31]. The vibration of the graphene resonator forms a movable capacitance together with a substrate; the electrical and mechanical energies are thus coupled [
31]. Another trend is to replace the optical or electrical cavity by a higher mode of the mechanical resonator [
32,
33]. Overhang is a simple but effective structure to realize a strong coupling between modes by bringing the two neighboring eigenfrequencies closer. For example, Okamoto et al. [
32] realized a strong mode coupling by using two modes of overhanging beams with the two closely spaced eigenfrequencies of 293.94 and 294.37 kHz, which also achieves a highly efficient energy transfer between the two modes. As shown in
Figure 1, overhang is a mechanical sledge shared by two or several (sub)structures [
34,
35,
36]. Overhang physically couples the (sub)structures and leads to the eigenfrequency veering and mode localization, which has been used as an effective mechanism to increase the sensitivity of mass sensing [
35,
36,
37,
38]. Furthermore, a veering locus also (approximately) corresponds to the maximum shift of an eigenfrequency [
37], which is a much sought-after property in a mass resonator sensor [
26,
27,
28].
Due to their simple design, overhanging structures are an important type of resonator structure [
32,
33,
34,
35,
36,
37,
38,
39,
40,
41]. Various numerical methods have been used to study the mode localization and eigenfrequency veering of the overhanging structures. However, most numerical discretization methods create or destroy the self-adjointness property of a continuous system [
42], which causes error and even unphysical (un)coupling, and therefore, leads to wrongful results. This problem becomes much more serious when two eigenfrequencies are very close to each other in the veering region [
43]. Even for a numerical discretization method which preserves the (non)self-adjointness property, such as the Galerkin method [
42], wrongful computational results can still be obtained, as demonstrated by Leissa [
44]. The finite element method fails in the computation of the three overhanging cantilevers due to the discretization asymmetry, which physically is induced by an adsorbate [
41]. The asymptotic series and perturbation method can also break down as they assume slight changes in the modes [
43]. Once a numerical discretization method is used for the computation of eigenfrequency veering, its validity is always a heated debate [
42,
44,
45]. In this study, the matrix for the eigenfrequency computation is analytically derived and the computation is carried out without any numerical discretization. A comprehensive study on the variations of the eigenfrequencies and mode shapes of two overhanging clamped–clamped beams is presented. In contrast with the deceptive simplicity of the two overhanging beams, rich variation patterns of the mode localization, eigenfrequencies and their veerings are revealed. For each eigenfrequency of a uniform beam, a new and larger eigenfrequency is generated due to the coupling effect of an overhanging structure, which resembles the above case of lattice dynamics [
1,
2]. As the overhang length increases, this newly generated eigenfrequency increases much more quickly, which causes the separation of the two eigenfrequencies. Consequently, the modes associated with these two separating eigenfrequencies undergo the so-called mode splitting [
14]. Furthermore, this new eigenfrequency also approaches the eigenfrequency of a higher mode and then diverges, which results in veering. The effects of the overhang geometry on the veering and mode localization are systematically studied. This comprehensive study should be of some help to the design of the overhanging resonator structures [
32,
33,
34,
35,
36,
37,
38,
39,
40,
41]. Recently, Kasai et al. [
46] proposed an efficient method called “virtual coupling”, which can effectively change the coupling strength of a system through a nonlinear feedback scheme. Their method has great potential to be applied to various coupled systems. Besides designing the overhanging structures, this paper also presents a rather general study on the properties of coupled systems, which can provide some help to other coupled systems, such as optomechanical [
21,
24,
29,
30,
31] and electromechanical [
33,
47,
48] ones.
2. Model Development
In the schematic diagram of
Figure 1, two beams are connected by two overhangs, which are clamped.
,
b and
h and
,
and
are the lengths, widths and thicknesses of the two overhangs, respectively. The widths and thicknesses of the two beams are
and
and
and
, respectively. The kinetic energy (
T) of the system is given as follows
where
L is the length of the whole structure;
m and
are the mass per unit length of the two overhangs; and
and
are the mass per unit length of Beams 1 and 2, respectively. For the rectangular sections,
,
,
and
(
and
and
and
are the densities of the two overhangs and two beams, respectively). The first and third terms of Equation (
1) are the kinetic energies of the two overhangs; the second term is those of the two beams. The corresponding transverse displacements of overhangs and beams are
w and
and
and
, respectively.
The potential/bending energy (
U) of the system is as follows:
where
E and
and
and
are the Young’s moduli of the overhangs and beams, respectively.
I and
and
and
are the second moments of area of the overhangs and beams, respectively. For rectangular cross sections,
,
,
and
. By applying Hamilton’s principle, i.e.,
, the following governing equations are derived:
The following sixteen boundary conditions are also derived from Hamilton’s principle:
The first four equations are the clamped boundary conditions for the two overhangs. The other twelve are at the two connecting points at
and
. Physically, these twelve equations are to ensure the continuity of displacement, slope, moment and shear force at the two connecting points [
19]. The governing equation of the overhang–beams–overhang structure is divided into four in three domains, as indicated by Equation (
3). In the domain of
, the two beams can have two different vibrations of
and
. However, the two overhangs are the parts shared by the two beams and, therefore, for each overhang, there is only one vibration, i.e.,
w for Overhang 1 and
for Overhang 2. The two overhangs couple the two beams as they share the same displacement in the overhang sections. This coupling is also reflected in the boundary conditions at the two connecting points, as given in Equation (
4). It is also noteworthy that Equations (
1)–(
4) are derived by a rather general method, which can be directly used to derive the governing equations of other type of overhanging structure with a small modification [
49].
The following dimensionless quantities are introduced [
50]
The governing equations of Equation (
3) are now non-dimensionalized as the following:
The dimensionless quantities in Equation (
6) are defined as follows:
Physically,
and
are the dimensionless mass per unit length and bending stiffness of Overhang 1 as compared with those of Beam 1.
and
are those of Overhang 2.
and
are the Beam 2 (dimensionless) deviations of mass per unit length and bending stiffness from those of Beam 1, respectively. By assuming
,
,
and
(
ia the dimensionless circular frequency) and substituting them into Equation (
6), the following solution forms are obtained
where
,
,
and
(
1–4) are the sixteen unknown constants to be determined by the sixteen boundary conditions. The dimensionless quantities of
,
,
and
are defined as the following:
The sixteen boundary conditions of Equation (
4) now become the dimensionless ones as follows:
By substituting the solution forms of Equation (
8) into the boundary conditions of (
10), an eigenvalue problem is formulated by setting the determinant of a
matrix zero, which is a nonlinear transcendental equation to be solved. Significant computation effort can be saved by utilizing the clamped boundary conditions of the two overhangs. By substituting Equation (
8) into the first four boundary conditions of Equation (
10), the following solution forms are obtained after some manipulations
Here,
,
,
and
are the four constants determined by the clamped boundary conditions of the two overhangs at
and 1, which are given in
Appendix A. By substituting Equation (
11) into the twelve boundary conditions of Equation (
10) at the connecting points of
and
, the following equation is obtained:
Here,
(
) is a
matrix and its elements are given in
Appendix A.
is a vector defined as
.
is the zero vector of
. By setting the determinant of
to zero, i.e., det(
) = 0, the eigenvalues (
s) are found numerically one by one from the lower order to the higher one by the Newton–Raphson method [
51]. For most structures, their two adjacent eigenvalues are well-separated. In contrast, the two adjacent eigenvalues of an overhanging structure are very close to each other at the veering loci. Therefore, extreme caution during the computation should be taken when solving this nonlinear transcendental equation. Otherwise, wrongful results can be easily generated during the eigenfrequency computation. After an eigenvalue is found and substituted into the above equation, the constants of
,
,
,
,
,
,
,
,
,
,
,
and
can be found by setting
. Once those constants are found and substituted into Equation (
11), the mode shape is also obtained. As given in
Appendix A,
contains both diagonal and off-diagonal terms. Mathematically, those off-diagonal terms contain the coupling information [
11,
16,
52]. Here, a noteworthy point is that though the eigenfrequency is numerically computed by setting det(
) = 0, there is no numerical discretization procedure involved.
3. Results and Discussion
In the above derivations for a general case, the material properties of the two overhangs and beams are different. For simplicity, we assume that they are made of a same material, i.e.,
and
, which is a common scenario of many micro/nano-resonators [
34,
35,
36].
In
Figure 2, the case of
,
and
is firstly examined. In this study, the two overhangs are with the same geometry except the one in Figure 9. The above parameters also indicate that the two beams are identical. In
Figure 2, the variations of the first six
s (
–6) as functions of the overhang length (
) are presented. Here,
is the square root of the
nth dimensionless eigenfrequency (
). In the vibration study, it is more common to use this
parameter rather than
[
50]. For the comparison reason, the first six eigenfrequency square roots of a uniform beam (
s) are also marked in
Figure 2 and their values are given in
Table 1. For the convenience of statement, instead of the eigenfrequency square root, we call
and
eigenfrequency hereafter. In
Figure 2, a distinct feature arises:
and
(
i = 1, 2 and 3) appear as a pair and then separate. When
, there is no overhang and the whole structure is two separate clamped–clamped beams without any coupling. Therefore,
,
and
are the first three eigenfrequencies of a uniform beam. As the overhang length increases,
and
begin to separate. When
, the two separate beams merge into one (uniform) beam. For a rectangular beam, its eigenfrequency
[
50], i.e., the beam width has no impact on its eigenfrequency, which is also the mathematical reason in the two cases of
and
,
s match
s. These two cases are also the benchmark tests to see if the eigenfrequency computation is correctly carried out. When the overhang length (
) approaches
, the beams are about to become one and the coupling is about to vanish,
approaches
and
approaches
. In summary,
and
start with
and they separate; finally,
ends up with
, and
with
. Physically,
s are the newly emerging eigenfrequencies due to the coupling of the two beams by the overhangs. Similarly, new eigenfrequencies also arise in the acoustic and optical bands of a crystal due to the introduction of an impurity [
1,
2]. In optomechanics, several important applications are associated with this phenomenon of newly emerging eigenfrequency/mode shape and their separations, which corresponds to the optomechanically induced transparency [
14].
The veering characteristics can be summarized from the results in
Figure 2: higher eigenvalues experience more veering loci. For example, there is no veering locus for
; only one for
marked as No. 1; two for
marked as Nos. 1 and 2; three for
marked as Nos. 2–4; and four for
marked as Nos. 3–6. Because only up to six eigenfrequencies are plotted, there are only two
veering loci marked as 5 and 6 shown in
Figure 2. Actually there are three more
veering loci as indicated by the arrows, which are formed by the
veering. Physically, the variation of the overhang length (
) changes the mass and stiffness distribution of the system and its impact on different mode is different [
19], which is clearly seen in the
wavy variations in
Figure 2. As the two overhangs here are identical, there is only one control parameter, i.e.,
. The eigenfrequency variation as the function of
is a curve and the corresponding veering is thus called curve veering [
42,
44]. When the number of the control parameter is two and the eigenvalue variation is a surface, the corresponding veering is thus called surface veering [
12]. Six eigenfrequency veering loci are identified and marked with Arabic numbers in
Figure 2. Locus No. 1 is the
veering; No. 2 is the
veering; Nos. 3 and 4 are the
veering; and Nos. 4 and 5 are the
veering. The smallest gap between two adjacent veering eigenvalues is called veering neck [
53] or veering width [
54]. Clearly, the veering necks of veering loci Nos. 1 and 3 are much smaller than the other three.
Figure 3 presents a closer look at these two veering loci.
Figure 3a plots the
veering around
, at which the veering width
is achieved. Similarly,
Figure 3b plots the
veering around
, at which the veering width
is achieved.
Figure 3 shows a distinctive feature of veering: The two adjacent eigenvalues rapidly converge and then abruptly diverge; however, they do not cross each other. Therefore, instead of veering, a more straightforward name of “avoided crossing” is used by physicists [
14]. It is worth mentioning that eigenvalue crossing rather than veering can occur for an overhanging structure [
49].
An important issue here is how to extract the effective strength of coupling from the the veering plots. In the modeling viewpoint, the coupling effect is incorporated in those six parameters as defined in Equation (
7) as well as the two parameters of the overhang locations (
and
). These parameters appear in both the governing equation of Equation (
6) and boundary conditions of Equation (
10), and, therefore, it is extremely difficult if not impossible to derive an analytical expression for the effective strength of coupling for this simple overhanging structure. A similar problem is also encountered in the model of internally coupled ears (ICE) [
55,
56], in which the two eardrums are modeled as two membranes connected by a cavity. In the ICE model, the coupling effect is determined by several parameters such as the tension, damping, radius, thickness and density of the membranes, and the air density, pressure and length of the cavity [
55]. The vibration of the membrane poses a movable boundary problem for the propagation of acoustic wave inside the cavity [
55,
56], which further complicates the problem and thus the evaluation on the coupling effect. A major difficulty in the development of the overhanging resonator sensor is to control/characterize the initial disorders [
40], which influence the coupling strength and result from the manufacturing error. However, those parameters in essence change the eigenfrequencies of a coupled system and the effective strength of coupling stiffness,
, can be effectively evaluated by the following equation [
57]
where
and
are the two (measured) adjacent eigenfrequencies and
.
k is the (oscillator) stiffness and in our case
k is the cantilever effective stiffness [
57], which is dependent on its mode(s) excited [
52]. According to Equation (
13), at the loci of the veering neck/width [
53,
54],
reaches its minimum and thus indicates the weakest physical coupling. However, this minimum
, i.e., the weakest physical coupling, signifies the strongest mode coupling [
15,
17]. Generally speaking, large
leads to the synchronization and small
leads to the (mode) localization [
58].
Figure 4 examines the variations of the mode shapes of
and
before and after the veering locus No. 1 around
.
Figure 4a,b plots the
mode shape variations at
and
, respectively. Because two beams can vibrate differently, there are two different deflections, which are normalized by being divided by the maximum displacement of Beam 1. At
, the peak of Beam 1 (approximately) corresponds to the valley of Beam 2, which is called that the two beams are out-of-phase [
38]. When
, a significant change occurs: Now, the peaks/valleys of Beam 1 (approximately) matches those of Beam 2, which is called that the two beams are in-phase [
38]. Similarly,
Figure 4c,d plots the
mode shape variations at
and
, respectively. Compared with those in
Figure 4a,b, things are reversed: the
mode shape experiences the transition of the in-phase mode to the out-of-phase mode. When
Figure 4a,c is examined together, physically, it means that, for a given overhanging structure with
, the lower
mode is out-of-phase and the higher
mode is in-phase. In the other configuration of overhanging structure with
, as shown in
Figure 4b,d, things are reversed again: the lower
mode is in-phase and the higher
mode is out-of-phase. This in-phase and out-of-phase pattern shift in the mode shape is often accompanied by significant amplitude variations [
40,
41].
Figure 5 examines the variations of the mode shapes of
and
before and after the veering locus No. 3 around
. Similar to the
mode shape variation before and after a veering locus in
Figure 4, the
mode shape in
Figure 5 experiences the transition of the out-of-phase mode to the in-phase mode. Meanwhile, the
mode shape in
Figure 5 experiences the transition of the in-phase mode to the out-of-phase mode, which is similar to the
mode in
Figure 4. In summary, before and after a veering locus, not only the two veering eigenvalues but also their mode shapes experience rapid changes. Because a small change of
results in large changes of eigenvalue and mode shape, eigenvalue veering is a catastrophic type phenomenon [
11]. Furthermore, the mode shape change from the in-phase mode to out-of-phase mode or vice versa is a qualitative one, which is expected to be of important use in the mass resonator sensors based on the mode shape variation [
35,
54].
Instead of just focusing the mode shape variations in a narrow area before and after a veering locus,
Figure 6 and
Figure 7 provide an overall and more comprehensive study on the mode shape variations. In
Figure 6 and
Figure 7, each row is the mode shapes of a
with different
s of
,
,
,
and
. The values of
s are given in
Table 1. For each row of a
mode shape, the mode shapes with different
s are plotted in the same scale for the comparison reason.
Figure 6 presents the mode shapes of the first three
s and
Figure 7 presents those of the last three. Two patterns can be summarized from the results in these two figures. The first is that the odd mode of
except
experiences the in-phase/out-of-phase mode transition around the first veering locus and the out-of-phase/in-phase transition around the second veering locus, and then the mode transitions of the in-phase/out-of-phase and out-of-phase/in-phase alternate. In contrast, the even mode of
reverses the transition patterns:
mode experiences the out-of-phase/in-phase mode transition around the first veering locus and the in-phase/out-of-phase transition around the second veering locus, and then the mode transitions of the out-of-phase/in-phase and in-phase/out-of-phase alternate. As discussed above,
Figure 4a,b shows the
out-of-phase/in-phase mode transition around its first and only veering locus marked as No. 1, which is around
. This
out-of-phase/in-phase mode transition around
is also captured in the
mode at
and that at
as presented in the second row of
Figure 6. However, as discussed in
Figure 2, a higher eigenvalue/mode experiences more than one veering locus. For example,
experiences two veering loci around
marked as No. 1 and around
marked as No. 2, as shown in
Figure 2. Again, the
mode in-phase/out-of-phase transition around the first veering locus of
is captured both by
Figure 4c,d and the
modes at
and
, as presented in
Figure 6. Around the second veering locus of
, the
mode out-of-phase/in-phase transition is captured in the
modes at
and
in
Figure 6.
The second one is that the wave numbers of all the
(
n = 1–6) modes except the
mode increase as
varies from 0 to
. The mechanism for this wave number increase can explained as follows: As shown in
Figure 2, the
value starts with
when
and ends up with
when
. Correspondingly, the
mode starts with the first mode of a uniform beam (the
mode) when
and ends up with the second mode of a uniform beam (the
mode) when
. Similarly, the
mode starts with the
mode and ends as the
mode, the
mode starts with the
mode and ends up with the
mode, and so on. The exception case of the
mode as shown in the first row of
Figure 6 always keep the
or
-like mode and the in-phase configuration. Again, it is explained by
Figure 2:
starts with
and ends up with the same
.
In
Figure 6 and
Figure 7, there are two benchmark characteristics of the mode shapes worth pointing out. The first one is that, at
, the deflection difference of two beams vanishes. The reason is simple: the two beams with two overhangs merge into one uniform beam and, therefore, there is only one beam deflection. At
, the structure configuration is two separate and independent beams. Therefore, the mode shape can be in two scenarios: The first one is the in-phase configurations of
, as shown in
Figure 6 and
Figure 7, in which the two-beam deflections overlap. The second one is the out-of-phase configuration, in which the two-beam deflections are in the exactly opposite phases. The second characteristic is that except the value of
where the two-beam deflections overlaps, the two-beam deflections are more or less different from other values of
. For example, in the first row of
Figure 7, there is a huge difference between the two-beam deflections of the
at
; in comparison, the difference is much smaller at
. This deflection difference of the two-beam in an eigenmode is called mode localization [
35] or confinement of vibration [
59]. Mode localization physically means that motion/energy is not equally distributed in each part of a system; the motion/energy concentration on one or some parts of a system can be harmful to cause an unexpected fatigue failure [
9].
It is worth explaining briefly how the mode localization is induced in the overhanging structure. As mentioned above, the localized modes in locally resonant materials are realized in a structure consisting of two materials with huge difference [
6]. In contrast, the two overhangs here are identical and two beams are also identical, a natural question arises: Why do the two beams vibrate differently? Many theories have been firmly established to show that the mode localization can also be caused by mistuning, i.e., irregularity of a periodic structure [
9,
11,
16,
59]. At the two connecting points of
and
, there are discontinuities of the mass and bending stiffness, which leads to a sudden change of the governing equations as indicated by Equation (
6). These discontinuities in essence are irregularities. When
and
, those irregularities are gone. When
, the two beams vibrate independently with the same shape; there is no coupling and thus no mode localization. When
, physically the two beams merge into one and mathematically the deflections of two beams overlap. Therefore, there is no mode localization, either. Overhang provides both the coupling mechanism and irregularity, which are responsible for the mode localization. The other worthy point is about the relation between the eigenfrequency veering and mode localization. Pierre presented the following statement: “Curve veering and strong mode localization are two manifestations of a single phenomenon.” [
11]. In many structures, the curve veering and (strong) mode localization indeed occur together. In this overhanging structure, the simultaneous occurrence of eigenfrequency veering and mode location is also seen. For example, in the first row of
Figure 7, there is a strong mode localization of the
mode at
, which is near the veering locus No. 2 around
. A strong mode localization of the
mode at
is also seen. However, there is no eigenfrequency veering for
at all, as shown in
Figure 2. In contrast to Pierre’s viewpoint, Natsiavas [
10] argued that the eigenfrequency veering and mode localization are not a same thing. In his analysis of a two-degree-freedom linear oscillator, Natsiavas [
10] found that the mode localization occurs for the majority of the combinations of the system parameters and in contrast, the eigenfrequency veering only occurs within a very narrow range of the system parameters. In conjunction with the eigenfrequency veering loci in
Figure 2 and mode shape variations in
Figure 6 and
Figure 7, we agree with Natsiavas on that in this particular overhanging structure, the eigenfrequency veering and mode localization are two different things; (strong) mode localization can occur without an eigenfrequency veering.
Figure 8 examines the impact of the overhang thickness on the eigenfrequency veering. Compared with the dimensions of the structure in
Figure 2, all the geometric parameters are kept the same except that the overhang thickness is doubled as
. A distinctive difference between
Figure 2 and
Figure 8 is that
now ends up with
at
, although
and
still start with
. Again, at
, the overhanging structure becomes a uniform beam and its corresponding frequency is
. As defined in Equation (
9),
. Therefore,
. Another distinctive difference is that the number of the eigenfrequency veering is significantly reduced. In
Figure 8, only three veering loci are identified.
Figure 9 examines the asymmetry of the two overhang lengths on the eigenfrequency veering. Compared with the dimensions in
Figure 8, all the geometric parameters are the same except that
is set in
Figure 9. In all other figures,
is set, which is a symmetric configuration. Now, with the asymmetric configuration, more veering loci arise in
Figure 9. In conjunction with
Figure 2,
Figure 8 and
Figure 9, it is clear that the geometry of the overhangs play the paramount role in the eigenfrequency veering, which determines the number of veering, veering loci and neck.
Figure 10 and
Figure 11 study the influences of
and
, respectively. In
Figure 10 and
Figure 11,
and
are set. In
Figure 10,
is fixed and
varies from
to
. Physically,
defined in Equation (
7) is the difference of the mass per unit length between the two beams. Larger
means more mass of the system. Therefore, all
s in
Figure 10 decrease monotonically with the increase of
. In
Figure 11,
is fixed and
varies from
to
.
is also defined in Equation (
7) and physically it is the difference of the bending stiffness between the two beams. Larger
means larger system stiffness. Therefore, all
s in
Figure 11 increase monotonically with the increase of
. From the results in
Figure 10 and
Figure 11, it is concluded that
and
can have significant impact on some eigenfrequencies of
s. The dramatic decrease and increase of
in
Figure 10 and
Figure 11 are seen. However, there is no veering behavior of abrupt converging and then diverging. The variations of
and
only change the gap distance of two (adjacent) veering eigenfrequencies, i.e.,
, and thus their veering neck. For example, the gap distance between
and
in
Figure 10 decreases with the
increase. In contrast, the gap distance between
and
increases with the
increase in
Figure 11. The setting of
in
Figure 10 and
Figure 11 is close to the veering locus of
marked as No. 1 in
Figure 8, which is the
veering. It is also noticed that, as
or
varies, the gap distance between
and
is always kept the smallest among all the gap distances.