Next Article in Journal
On Dynamic Investigations of Cournot Duopoly Game: When Firms Want to Maximize Their Relative Profits
Previous Article in Journal
Analysis of Sound Absorption Characteristics of Acoustic Ducts with Periodic Additional Multi-Local Resonant Cavities
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Vibration Localization and Anti-Localization of Nonlinear Multi-Support Beams with Support Periodicity Defect

1
Department of Mechanics, School of Aeronautics and Astronautics, Zhejiang University, Hangzhou 310027, China
2
Hong Kong Branch of National Rail Transit Electrification and Automation Engineering Technology Research Center, Department of Civil and Environmental Engineering, The Hong Kong Polytechnic University, Hung Hom, Kowloon, Hong Kong 999077, China
*
Author to whom correspondence should be addressed.
Symmetry 2021, 13(12), 2234; https://doi.org/10.3390/sym13122234
Submission received: 25 October 2021 / Revised: 5 November 2021 / Accepted: 12 November 2021 / Published: 23 November 2021

Abstract

:
A response analysis method for nonlinear beams with spatial distribution parameters and non-periodic supports was developed. The proposed method is implemented in four steps: first, the nonlinear partial differential equation of the beams is transformed into linear partial differential equations with space-varying parameters by using a perturbation method; second, the space-varying parameters are separated into a periodic part and a non-periodic part describing the periodicity defect, and the linear partial differential equations are separated into equations for the periodic and non-periodic parts; third, the equations are converted into ordinary differential equations with multiple modes coupling by using the Galerkin method; fourth, the equations are solved by using a harmonic balance method to obtain vibration responses, which are used to discover dynamic characteristics including the amplitude–frequency relation and spatial mode. The proposed method considers multiple vibration modes in the response analysis of nonlinear non-periodic structures and accounts for mode-coupling effects resulting from structural nonlinearity and parametric non-periodicity. Thus, it can handle nonlinear non-periodic structures with a high parameter-varying wave in wide frequency vibration. In numerical studies, a nonlinear beam with non-periodic supports (resulting in non-periodic distribution parameters or periodicity defect) under harmonic excitations was explored using the proposed method, which revealed some new dynamic response characteristics of this kind of structure and the influences of non-periodic parameters. The characteristics include remarkable variation in frequency response and spatial mode, and in particular, vibration localization and anti-localization. The results have potential applications in vibration control and the support damage detection of nonlinear structures with non-periodic supports.

1. Introduction

Many engineering structures such as continuous girder viaducts and synchrotron radiation bridges can be simplified as a beam with multi-supports in dynamics analysis. The beam is usually symmetric and becomes asymmetric after a change in the stiffness of some supports, e.g., stiffness reduction due to damage. The changeability of dynamic characteristics such as frequency response and spatial mode (or mode component) is important for structural optimization design (or vibration control) and anomaly diagnosis (or damage detection). Periodic structures with spatial distribution parameters such as periodic geometric parameters, periodic physical parameters and periodic constraints have special dynamics which are different from those of common structures with uniform distribution parameters. Parametric periodicity has outstanding influence on structural dynamic responses because of the coupling of periodic time-varying (or harmonic) response and periodic distribution parameters. Numerous studies have reported characteristic frequencies [1,2,3,4], modal localization and buckling [5,6,7,8,9,10,11,12,13,14], quasi-periodic distribution parameter effects [15,16,17,18] and application [19,20,21,22,23,24,25,26,27] based on transfer matrix method, spatial (Bloch) harmonic expansion method, (Floquet–Bloch and Galerkin) double expansion method and finite element method, etc. Waves and vibration in beams with non-uniform distribution parameters have also been pursued using a fundamental solution, semi-analysis and finite element methods [28,29,30,31,32,33,34]. Perturbation and multiple scale methods were applied for weakly periodic parameter cases [35,36]. The periodic beam is an important engineering structure or structural component, and its dynamics including finite-size and boundary effects have been addressed [37,38,39,40]. Compared with periodic structures, quasi-periodic structures are likely to afford highly variable dynamic responses because period-structural dynamics have certain instability and are sensitive to periodicity perturbation (or defect). However, in the aforementioned investigations, only the dynamics of linear periodic and quasi-periodic structures and the corresponding linear dynamic characteristics were addressed.
Under certain large loads, engineering structures will produce nonlinear vibrations because of inherent geometric and physical nonlinearity. The nonlinear dynamic characteristics of periodic and non-periodic structures are a subject of research interest. Research on the wave propagation of periodic structures modeled by nonlinear single-degree-of-freedom systems has been pursued based on periodicity conditions, in which the multiple scales method, perturbation method, linearization method and harmonic balance method were adopted [41,42,43,44,45,46,47]. A nonlinear periodic structure with cycle symmetry was studied by converting it into independent nonlinear single-degree-of-freedom systems using nonlinear normal modes for frequency response analysis, and its mode localization under weak coupling was discovered by asymptotic approximation [48,49]. The nonlinear vibration of non-uniform beams and dynamic stability of sandwich beams on nonlinear elastic foundations have also been studied using the finite element and harmonic balance methods, Galerkin and iterative numerical methods, linearization and Chebyshev collocation methods, Bolotin method and Hamiltonian approach, respectively [50,51,52,53,54]. A differential quadrature method was applied to the nonlinear vibration analysis of sandwich beams where vibration modes were obtained using a numerical iterative algorithm [55,56]. Non-periodic dynamics of complex systems with fractional derivatives have been studied [57,58]. However, only a few studies have addressed the nonlinear dynamics of finite-size periodic beams. An asymptotic tolerance averaging technique was applied to the nonlinear vibration analysis of a beam with geometric nonlinearity and periodic distribution parameters, but this study was limited to low-frequency vibration and a high parameter-varying wave because of the averaging in a periodic unit [59,60]. Those analysis methods have limitations in the number of used vibration modes which can affect the accuracy and reliability of results as the nonlinearity and periodicity will induce the coupling of linear structural modes. Recently, a multimode perturbation method for the frequency response analysis of nonlinear beams with periodic distribution parameters has been proposed, and some novel dynamic response characteristics including remarkable variations of amplitude–frequency relation and spatial mode of nonlinear periodic beams were obtained [61]. However, the study only considered the dynamics of periodic beams. The dynamics of nonlinear non-periodic beams, e.g., with non-periodic supports (or support periodicity defect, that can be viewed as variation of periodicity), need to be further studied and to this end, the corresponding analysis method needs to be developed. Non-periodic or quasi-periodic beams are more realistic structures in engineering fields, and they would have dynamic responses which would greatly differ from those of periodic beams based on the results of linear periodic and quasi-periodic structures. This study of the dynamics of nonlinear non-periodic structures also has potential applications in vibration control and support damage detection [14,15,18,20,21,22,23,62,63,64,65,66].
In this paper, the vibrational amplitude–frequency characteristics of nonlinear multi-support beams with non-periodic supports (resulting in non-periodic distribution parameters or periodicity defect) under harmonic excitations was studied. A response analysis method for nonlinearly vibrational beams with spatial distribution parameters and non-periodic supports was developed by combining perturbation, periodic and non-periodic separation, Galerkin expansion and harmonic balance methods. First, a partial differential equation with space-varying parameters is elicited for a nonlinear non-periodic beam. The nonlinear equation is transformed into a set of linear partial differential equations by using a perturbation method. Second, the spatial distribution parameters are separated into periodic part and non-periodic part describing the periodicity defect. The linear partial differential equations are separated into a set of differential equations for the periodic beam with periodic part parameters and another set of differential equations for the non-periodic beam with the non-periodic part parameters. Third, the equations are converted into ordinary differential equations with coupling multiple modes by using the Galerkin method. Multiple vibration modes are considered and thus the equations are suited for the nonlinear vibration of non-periodic structures with high parameter-varying wave in a wide frequency band. Fourth, the ordinary differential equations are solved by using a harmonic balance method to obtain vibration responses of the beam with non-periodic parameters under harmonic excitations, which are used for characteristics analysis of amplitude–frequency relation and spatial mode. Finally, resulting graphs for a non-periodic beam with various support stiffness and damping are provided to illustrate effective application of the proposed method and to reveal frequency response characteristics. New response characteristics including remarkable variations in frequency response and spatial mode for the nonlinearly vibrational beam with different non-periodic supports, particularly vibration localization and vibration anti-localization, were discovered for the first time.

2. Nonlinear Vibration Equation of Beam with Non-Periodic Supports

Consider a horizontal elastic beam with spatial distribution parameters and non-periodic supports under vertical external excitation, as shown in Figure 1. Its material is isotropic, and normal stress in the z direction is small and neglected, where z denotes a vertical coordinate. Vertical displacement is assumed as unvarying with z and the cross-section is perpendicular to the central axis as the initial x, where x denotes the horizontal coordinate. Longitudinal and rotational inertias are small and neglected. For large excitation acting on a long beam, the vertical vibration in the (x,z) plane exhibits certain geometric nonlinearity, and the nonlinear longitudinal normal strain can be expressed as
ε s = z 2 w x 2 + 1 2 ( w x ) 2 + 1 2 z 2 w x 2 ( w x ) 2
where w is vertical displacement. Longitudinal normal stress is σs = s, where E is Young’s modulus. Axial force FN, shear force FS and bending moment M of the beam can be calculated using the normal stress σs. For a long beam in low and middle frequency vibration, a mechanical model based on the Bernoulli–Euler theory can be applied. Considering the vertical middle support forces, vertical inertia force and damping force, the nonlinear vibration equation of the beam is obtained as
ρ A 0 2 w t 2 + c d w t + x [ F S cos ( w x ) ] cos ( w x ) x [ F N sin ( w x ) ] cos ( w x ) + i = 1 N p ( c s i w t + k s i w ) δ ( x x s i ) = f e ( x , t )
where ρ is the mass density, A0 is the cross-sectional area, cd is the damping coefficient, ksi, csi, xsi are, respectively, the stiffness, damping and coordinate of the ith support except two end constraints, δ(·) is the Dirac-delta function, Np is the support number, fe is the external force and t is the time variable. Equation (2) can be simplified by neglecting high-order nonlinear terms and using the condition that the beam height (h) is much smaller than beam length (L). The non-dimensional nonlinear vibration equation is [61]:
m b 2 w ¯ t 2 + c b w ¯ t + 2 x ¯ 2 ( E I ¯ 2 w ¯ x ¯ 2 ) x ¯ [ E A ¯ ( w ¯ x ¯ ) 3 ] + i = 1 N p ( c s i w ¯ t + k s i w ¯ ) δ ( x ¯ x ¯ s i ) = F ( x ¯ , t )
where:
x ¯ = x / L ,   w ¯ = w / W a ,   m b = ρ A 0 L ,   c b = c d L ,   I ¯ = I 0 / L 3 ,   A ¯ = A 0 W a 2 / 2 L 3 ,   F = f e L / W a
in which I0 is the cross-sectional second moment and Wa is a small displacement amplitude under unit excitation. For the beam with uniform distribution parameters, mb, E A ¯ and E I ¯ are constants. For the case of non-periodic supports, there are supports whose coordinates x ¯ s i i / ( N p + 1 ) or the stiffness ksi and damping csi differ from the others.
The partial differential Equation (3) describes the nonlinear vibration of the beam with spatial distribution parameters and non-periodic supports. Under Wa << 1, rewrite nonlinear coefficient as E A ¯ = ε B 10 / 3 , where ε is a small parameter. Equation (3) can be expressed in the following form to facilitate a common solution:
B 1 2 w ¯ t 2 + B 2 w ¯ t + k = 1 3 B k + 2 ( 5 k ) w ¯ x ¯ ( 5 k ) + ε B 10 2 w ¯ x ¯ 2 ( w ¯ x ¯ ) 2 + ε B 11 ( w ¯ x ¯ ) 3 + i = 1 N p ( D 1 i w ¯ t + D 2 i w ¯ ) δ ( x ¯ x ¯ i ) = F ( x ¯ , t )
where coefficients Bk and Dji are determined by Equation (3), and x ¯ i = x ¯ s i . Consider two ends of the beam having simply supported constraints. The boundary conditions are:
w ¯ ( 0 , t ) = w ¯ ( 1 , t ) = 0 ,   2 w ¯ ( 0 , t ) x ¯ 2 = 2 w ¯ ( 1 , t ) x ¯ 2 = 0
Vibration responses can be obtained by solving Equation (5) with conditions (6). However, Equation (5) is a nonlinear partial differential equation with space-varying parameters. Conventional analysis methods can only consider a few modes and actually use several coupled nonlinear ordinary differential equations. The effect of limited modes on response accuracy increases with the wave number of varying parameters and supports and vibration frequency. Therefore, an improved analysis method [61] combined with periodic and non-periodic separation analysis method was developed herein. The analysis method can use multiple vibration modes and is suitable for the nonlinear vibration of non-periodic structures with a high parameter-varying wave in a wide frequency band, which is described in the following sections.

3. Perturbation and Separation Solutions to Nonlinear Vibration of Non-Periodically Supported Beam

A perturbation method was first applied to transform nonlinear partial differential Equation (5) into a set of linear partial differential equations. The dimensionless displacement w ¯ , which depends on the small parameter ε, can be expanded as
w ¯ ( x ¯ , t ) = j = 1 ε j 1 w j ( x ¯ , t )
where wj is the jth-order approximation. Substituting Equation (7) into Equation (5) (see Appendix A as Equation (A1)) and balancing each power term of the small parameter yield the following linear differential equations:
B 1 2 w 1 t 2 + B 2 2 w 1 t 2 + k = 1 3 B k + 2 ( 5 k ) w 1 x ¯ ( 5 k ) + i = 1 N p ( D 1 i w 1 t + D 2 i w 1 ) δ ( x ¯ x ¯ i ) = F
B 1 2 w 2 t 2 + B 2 w 2 t + k = 1 3 B k + 2 ( 5 k ) w 2 x ¯ ( 5 k ) + i = 1 N p ( D 1 i w 2 t + D 2 i w 2 ) δ ( x ¯ x ¯ i ) = B 10 2 w 1 x ¯ 2 ( w 1 x ¯ ) 2 B 11 ( w 1 x ¯ ) 3
B 1 2 w 3 t 2 + B 2 w 3 t + k = 1 3 B k + 2 ( 5 k ) w 3 x ¯ ( 5 k ) + i = 1 N p ( D 1 i w 3 t + D 2 i w 3 ) δ ( x ¯ x ¯ i ) = B 10 [ 2 w 2 x ¯ 2 ( w 1 x ¯ ) 2 + 2 2 w 1 x ¯ 2 w 1 x ¯ w 2 x ¯ ] 3 B 11 ( w 1 x ¯ ) 2 w 2 x ¯
The above equations can be solved to obtain an asymptotic solution. The vibration displacement is generally determined by the first several approximate wj. For the beam with spatial distribution parameters and non-periodic supports, Bk, Dji, x ¯ i and F can be regarded as periodicity defect parameters. The parameters are separated into a spatial periodic part and non-periodic part (which describes the periodicity defect) as
B k = B p k + Δ B k ,   D 1 i = D p 1 + Δ D 1 i ,   D 2 i = D p 2 + Δ D 2 i ,   y i = y p i + Δ y i ,   F = F p + Δ F
where subscript “p” denotes the periodic part and “Δ” denotes the non-periodic part. The separation yields a periodic beam with the periodic part parameters and a non-periodic beam with non-periodic part parameters (or periodicity defect). The approximate displacement wj can be separated accordingly into:
w j = w p j + Δ w j
Substituting expressions (11) and (12) into Equations (8)–(10) yields a set of differential equations for the periodic part displacement wpj of the periodic beam and another set of differential equations for the non-periodic part displacement Δwj of the non-periodic beam. The linear partial differential equations for the periodic part are:
B p 1 2 w p 1 t 2 + B p 2 2 w p 1 t 2 + k = 1 3 B p , k + 2 ( 5 k ) w p 1 x ¯ ( 5 k ) + i = 1 N p ( D p 1 w p 1 t + D p 2 w p 1 ) δ ( x ¯ x ¯ p i ) = F p
B p 1 2 w p 2 t 2 + B p 2 w p 2 t + k = 1 3 B p , k + 2 ( 5 k ) w p 2 x ¯ ( 5 k ) + i = 1 N p ( D p 1 w p 2 t + D p 2 w p 2 ) δ ( x ¯ x ¯ p i ) = B p 10 2 w p 1 x ¯ 2 ( w p 1 x ¯ ) 2 B p 11 ( w p 1 x ¯ ) 3
B p 1 2 w p 3 t 2 + B p 2 w p 3 t + k = 1 3 B p , k + 2 ( 5 k ) w p 3 x ¯ ( 5 k ) + i = 1 N p ( D p 1 w p 3 t + D p 2 w p 3 ) δ ( x ¯ x ¯ p i ) = B p 10 [ 2 w p 2 x ¯ 2 ( w p 1 x ¯ ) 2 + 2 2 w p 1 x ¯ 2 w p 1 x ¯ w p 2 x ¯ ] 3 B p 11 ( w p 1 x ¯ ) 2 w p 2 x ¯
The linear partial differential equations for the non-periodic part (or periodicity defect) are:
B 1 2 Δ w 1 t 2 + B 2 Δ w 1 t + k = 1 3 B k + 2 ( 5 k ) Δ w 1 x ¯ ( 5 k ) + i = 1 N p ( D 1 i Δ w 1 t + D 2 i Δ w 1 ) δ ( x ¯ x ¯ i ) = Δ F Δ B 1 B p 1 1 F p + ( Δ B 1 B p 1 1 B p 2 Δ B 2 ) w p 1 t + k = 1 3 ( Δ B 1 B p 1 1 B p , k + 2 Δ B k + 2 ) ( 5 k ) w p 1 x ¯ ( 5 k ) + i = 1 N p [ ( 1 + Δ B 1 B p 1 1 ) ( D p 1 w p 1 t + D p 2 w p 1 ) δ ( x ¯ x ¯ p i ) ( D 1 i w p 1 t + D 2 i w p 1 ) δ ( x ¯ x ¯ i ) ]
B 1 2 Δ w 2 t 2 + B 2 Δ w 2 t + k = 1 3 B k + 2 ( 5 k ) Δ w 2 x ¯ ( 5 k ) + i = 1 N p ( D 1 i Δ w 2 t + D 2 i Δ w 2 ) δ ( x ¯ x ¯ i ) = ( Δ B 1 B p 1 1 B p 2 Δ B 2 ) w p 2 t + k = 1 3 ( Δ B 1 B p 1 1 B p , k + 2 Δ B k + 2 ) ( 5 k ) w p 2 x ¯ ( 5 k ) + i = 1 N p [ ( 1 + Δ B 1 B p 1 1 ) ( D p 1 w p 2 t + D p 2 w p 2 ) δ ( x ¯ x ¯ p i ) ( D 1 i w p 2 t + D 2 i w p 2 ) δ ( x ¯ x ¯ i ) ] [ B 10 2 w 1 x ¯ 2 ( w 1 x ¯ ) 2 ( 1 + Δ B 1 B p 1 1 ) B p 10 2 w p 1 x ¯ 2 ( w p 1 x ¯ ) 2 ] [ B 11 ( w 1 x ¯ ) 3 ( 1 + Δ B 1 B p 1 1 ) B p 11 ( w p 1 x ¯ ) 3 ]
B 1 2 Δ w 3 t 2 + B 2 Δ w 3 t + k = 1 3 B k + 2 ( 5 k ) Δ w 3 x ¯ ( 5 k ) + i = 1 N p ( D 1 i Δ w 3 t + D 2 i Δ w 3 ) δ ( x ¯ x ¯ i ) = ( Δ B 1 B p 1 1 B p 2 Δ B 2 ) w p 3 t + k = 1 3 ( Δ B 1 B p 1 1 B p , k + 2 Δ B k + 2 ) ( 5 k ) w p 3 x ¯ ( 5 k ) + i = 1 N p [ ( 1 + Δ B 1 B p 1 1 ) ( D p 1 w p 3 t + D p 2 w p 3 ) δ ( x ¯ x ¯ p i ) ( D 1 i w p 3 t + D 2 i w p 3 ) δ ( x ¯ x ¯ i ) ] B 10 [ 2 w 2 x ¯ 2 ( w 1 x ¯ ) 2 + 2 2 w 1 x ¯ 2 w 1 x ¯ w 2 x ¯ ] + ( 1 + Δ B 1 B p 1 1 ) B p 10 [ 2 w p 2 x ¯ 2 ( w p 1 x ¯ ) 2 + 2 2 w p 1 x ¯ 2 w p 1 x ¯ w p 2 x ¯ ] 3 [ B 11 ( w 1 x ¯ ) 2 w 2 x ¯ ( 1 + Δ B 1 B p 1 1 ) B p 11 ( w p 1 x ¯ ) 2 w p 2 x ¯ ]
These equations can be solved to obtain the periodic part and non-periodic part of the asymptotic solution. The boundary conditions (6) can be expanded and separated, corresponding to displacements (7) and (12) for determining the linear vibration modes of the beam. The linear partial differential Equations (13)–(18) for the beam with partial distribution parameters and non-periodic supports can be further converted into linear ordinary differential equations by using the modes for obtaining modal displacement responses.

4. Mode Expansions of Asymptotic Solutions to Linear Partial Differential Equations for Periodic and Non-Periodic Beams

The Galerkin method was applied to convert the linear partial differential Equations (13)–(18) into a set of ordinary differential equations. Spatial distribution parameters and non-periodic supports of the beam result in the coupling of the differential equations, and multiple vibration modes need to be used in the transformation. The linear vibration modes of the beam can be determined based on boundary conditions (6). The dimensionless displacement w ¯ and then approximate displacement wj can be expanded into series using the modes. The approximate displacement expansion is:
w j ( x ¯ , t ) = n = 1 N m ϕ n ( x ¯ ) q j n ( t ) = Φ ( x ¯ ) Q j ( t )
where φn is the nth linear vibration mode, qjn is the nth modal displacement of the jth-order approximation, Φ = [φ1, φ2, …, φNm] is the mode vector, Qj = [qj1, qj2, …, qjNm]T is the modal displacement vector and Nm is mode number. For the simply supported beam with boundary conditions (6), the mode function is ϕ n ( x ¯ ) = sin n π x ¯ . The periodic part and non-periodic part of the displacement expansion (19) based on Equation (12) are:
w p j ( x ¯ , t ) = n = 1 N m ϕ n ( x ¯ ) q p j n ( t ) = Φ ( x ¯ ) Q p j ( t )
Δ w j ( x ¯ , t ) = n = 1 N m ϕ n ( x ¯ ) Δ q j n ( t ) = Φ ( x ¯ ) Δ Q j ( t )
where qpjn is the periodic part of the modal displacement qjn corresponding to the periodic beam, Δqjn is the non-periodic (or periodicity defect) part of the modal displacement qjn corresponding to the non-periodic beam, Qpj = [qpj1, qpj2, …, qpjNm]T and ΔQj = [Δqj1, Δqj2, …, ΔqjNm]T are the periodic part and non-periodic part vectors of modal displacements, respectively.
According to the Galerkin method, substituting expression (20) into Equations (13)–(15), multiplying the equations by ΦT and integrating them with respect to x ¯ yield the following linear ordinary differential equations in matrix form for the periodic part:
M p 1 2 Q p j t 2 + C p 1 Q p j t + K p 1 Q p j = F p j ( t )
where Mp1, Cp1, Kp1 are the modal mass, damping and stiffness matrices for the periodic part, respectively, and Fpj is the modal excitation vector for the periodic part. Their expressions are given in Appendix A as Equation (A2).
Similarly, substituting expressions (20) and (21) into Equations (16)–(18) yields the following linear ordinary differential equations in matrix form for the non-periodic part:
M 1 2 Δ Q j t 2 + C 1 Δ Q j t + K 1 Δ Q j = F j ( t )
where M1, C1, K1 are the modal mass, damping and stiffness matrices for the non-periodic part, respectively, and Fj is the modal excitation vector for the non-periodic part. Their expressions are given in Appendix A as Equation (A3).
The linear ordinary differential Equations (22) and (23) have coupling damping and stiffness, and can incorporate multiple vibration modes in accordance with Equations (A2) and (A3). The modal displacement Qj (Qpj with ΔQj) and then approximate displacement wj (wpj with Δwj) can be obtained by solving these equations. Solution accuracy for the high parameter-varying wave and high frequency vibration will increase with the number of used modes. The effect of support periodicity defect on the vibration response of the beam can be explored based on the displacement solution with the periodicity defect part of the non-periodic beam.

5. Harmonic Balance Solutions to Ordinary Differential Equations for Multi-Mode Coupling Vibrations of Periodic and Non-Periodic Beams

To obtain amplitude–frequency characteristics, the beam is considered subjected to harmonic excitation F. Modal excitation Fp1 in Equation (22) is harmonic due to F or Fp based on Equation (A2). Let Fp1 be:
F p 1 ( t ) = E p 11 s sin ω t
where Ep11s = [Ep11s,1, Ep11s,2, …, Ep11s,Nm]T is a constant vector and ω is the excitation frequency. According to the harmonic balance method, the stationary response of the modal displacement of the first-order approximation for the periodic beam is obtained by Equation (22) (j = 1) as
Q p 1 = R p 11 s sin ω t + R p 11 c cos ω t
where Rp11s = [Rp11s,1, Rp11s,2, …, Rp11s,Nm]T, Rp11c = [Rp11c,1, Rp11c,2, …, Rp11c,Nm]T and:
R p 11 = [ R p 11 s R p 11 c ] = Δ p 1 1 [ E p 11 s 0 ] ,   Δ p 1 = [ K p 1 ω 2 M p 1 ω C p 1 ω C p 1 K p 1 ω 2 M p 1 ]
By using Equations (A2) and (22) for j = 2, 3, …, modal excitations Fp2, Fp3, … can be determined and stationary responses of modal displacements of the high-order approximations Qp2, Qp3, … for the periodic beam can be obtained. They are:
F p j = l = 1 , 3 2 j 1 ( E p j l s sin l ω t + E p j l c cos l ω t )
Q p j = l = 1 , 3 2 j 1 ( R p j l s sin l ω t + R p j l c cos l ω t )
where Epjls = [Epjls,1, Epjls,2, …, Epjls,Nm]T and Epjlc = [Epjlc,1, Epjlc,2, …, Epjlc,Nm]T are constant vectors, Rpjls = [Rpjls,1, Rpjls,2, …, Rpjls,Nm]T, Rpjlc = [Rpjlc,1, Rpjlc,2, …, Rpjlc,Nm]T, and:
R p j l = [ R p j l s R p j l c ] = Δ p l 1 [ E p j l s E p j l c ] ,   Δ p l = [ K p 1 l 2 ω 2 M p 1 l ω C p 1 l ω C p 1 K p 1 l 2 ω 2 M p 1 ]
Similarly, for the non-periodic beam, modal excitations Fj (j = 1, 2, 3, …) can be determined and stationary responses of modal displacements of the approximations ΔQj (j = 1, 2, 3, …) for the non-periodic (or periodicity defect) part can be obtained by Equations (A3) and (23). They are:
F j = l = 1 , 3 2 j 1 ( E j l s sin l ω t + E j l c cos l ω t )
Δ Q j = l = 1 , 3 2 j 1 ( R j l s sin l ω t + R j l c cos l ω t )
where Ejls = [Ejls,1, Ejls,2, …, Ejls,Nm]T and Ejlc = [Ejlc,1, Ejlc,2, …, Ejlc,Nm]T are constant vectors, Rjls = [Rjls,1, Rjls,2, …, Rjls,Nm]T, Rjlc = [Rjlc,1, Rjlc,2, …, Rjlc,Nm]T, and:
R j l = [ R j l s R j l c ] = Δ l 1 [ E j l s E j l c ] ,   Δ l = [ K 1 l 2 ω 2 M 1 l ω C 1 l ω C 1 K 1 l 2 ω 2 M 1 ]
By using Equations (19)–(21), (28) and (31), frequency responses of the approximate displacement wj in multi-mode coupling vibration can be obtained as a function of ω and x ¯ , which is used to reveal amplitude–frequency characteristics and spatial mode. The response includes periodic part wpj (for the beam with periodic parameters and supports) and non-periodic part Δwj (for the beam with periodicity defect or non-periodic parameters and supports). Then, the dimensionless displacement response w ¯ of the nonlinearly vibrational beam with spatial distribution parameters and non-periodic supports is obtained by Equation (7) as
w ¯ ( x ¯ , t ) = j = 1 ε j 1 Φ ( x ¯ ) Q j ( t )
From Equation (12), the periodic part and non-periodic part of the displacement response are obtained as
w ¯ p ( x ¯ , t ) = j = 1 ε j 1 Φ ( x ¯ ) Q p j ( t )
Δ w ¯ ( x ¯ , t ) = j = 1 ε j 1 Φ ( x ¯ ) Δ Q j ( t )
For certain large nonlinear vibrations, principal harmonic response is dominant and the response amplitude is determined by constant vectors Rpj1s, Rpj1c, Rj1s and Rj1c (j = 1, 2, …) with the mode vector Φ. The above analysis method for the nonlinear vibration responses of non-periodic beams has no limit on the number of modes and is thus suitable for accurate numerical computation. The proposed method is particularly applicable to the nonlinear vibration of non-periodic structures with a high parameter-varying wave in wide frequency vibration. The nonlinear vibration of non-periodic structures under harmonic excitation has the coupling of multiple modes, and its dynamic stability is related to the parametrically excited stability of a multi-degree-of-freedom system with periodic time-varying parameters. The stability problem can be solved by using a direct eigenvalue analysis approach [61].

6. Results on Response Amplitude–Frequency Characteristics of Nonlinear Multi-Support Beam with Support Periodicity Defect

To illustrate the application of the proposed method and the response characteristics of non-periodic structures, consider a simply supported beam with a uniform distribution parameters and periodicity defect supports under harmonic excitation which models, e.g., a continuous girder viaduct. The nonlinear vibration of the beam is described by Equation (3) where mb, E A ¯ and E I ¯ are constants. In the case of periodic supports, the dimensionless coordinate of the ith support is x ¯ s i = i / ( N p + 1 ) (i = 1, 2, …, Np). By using dimensionless time τ = t/ω0, Equation (3) can be rewritten as
2 w ¯ τ 2 + c u w ¯ τ + 1 π 4 4 w ¯ x ¯ 4 k u 2 w ¯ x ¯ 2 ( w ¯ x ¯ ) 2 + i = 1 N p ( c u s i w ¯ τ + k u s i w ¯ ) δ ( x ¯ x ¯ s i ) = F u ( x ¯ , τ )
where:
ω 0 = π 2 E I 0 ρ A 0 L 4 ,   k u = 3 W a 2 A 0 2 π 4 I 0 ,   c u = c b m b ω 0 ,   k u s i = k s i L 3 π 4 E I 0 ,   c u s i = c s i ω 0 L 3 π 4 E I 0 ,   F u = F L 3 π 4 E I  
By comparing Equation (36) with Equation (5), there are B1 = 1, B2 = cu, B3 = 1/π4, B4 = B5 = 0, B10 = −ku/ε, B11 = 0, D1i = cusi, D2i = kusi and F is replaced by Fu. Boundary conditions are as given by Equation (6). Equation (36) is a nonlinear partial differential equation with spatial varying parameters due to non-periodic supports. The proposed method was applied to obtain the dimensionless displacement responses of the non-periodic beam in nonlinear multi-mode coupling vibration. The response amplitude is a function of dimensionless frequency ω and coordinate x ¯ , which is used to explore the characteristics of the amplitude–frequency relation and spatial mode.
Harmonic excitation on point ye is considered as F u = F 0 δ ( x ¯ y e ) sin ω τ . For two middle supports, Np = 2. Basic parameter values are ku = 0.1, ε = 0.01, kusi = 0.1, cusi = 0.02, ye = 0.4, F0 = 0.01, and the modal damping coefficient is 0.04 [14]. Results obtained by the proposed method for linear vibration are verified by exact theoretical solution and those for nonlinear vibration are verified by direct numerical simulation. The first four dimensionless natural frequencies of the linearly vibrational beam without damping and middle supports are 1, 4, 9 and 16. For dimensionless frequency ω ∈ [0, 20], the corresponding dimensionless response covers the first four resonances, and its accurate results are obtained by using the proposed method with the mode number Nm = 10 which corresponds to the natural frequency with the upper limit ω = 100.

6.1. Effects of Support Periodicity Defect on Response Amplitude–Frequency Characteristics

Results on dimensionless response amplitude–frequency relations of the nonlinearly vibrational beam with periodic supports and with non-periodic supports are given for comparison. Figure 2 shows the amplitude–frequency relations on the point x ¯ = 0.4 of the beam with periodic supports and with non-periodic supports (i.e., the left support deviation to the right Δ x ¯ s1 = 0.1/3 and left support deviation to the left Δ x ¯ s1 = −0.1/3) under the harmonic excitation. It can be seen that the first response peak increases with local span due to the support deviation. Figure 3 illustrates the relative variation of the response amplitudes of the non-periodic beam with left support deviation to the right (Δ x ¯ s1 = 0.1/3) to the beam with periodic supports under harmonic excitation. The first response peak is reduced from 0.115 to 0.110 at the dimensionless frequency ω = 1.14 with a relative reduction of 4.52%, and the maximal relative reduction is 6.06% at ω = 1.10. Therefore, the small deviation of periodic supports has a slight effect on the response amplitudes of the nonlinearly vibrational beam.
Figure 4 shows the amplitude–frequency relations on the point x ¯ = 0.4 of the beam with periodic supports and with non-periodic supports (i.e., left support stiffness kus1 = 0.6 and left support stiffness kus1 = 1.1) under the harmonic excitation. It can be seen that the first response peak decreases with the increase in the support stiffness and the first resonant frequency increases with the support stiffness. As the left support stiffness increases from 0.1 to 1.1, the first response peak is reduced from 0.115 (at ω = 1.14) to 0.068 (at ω = 1.62) and the relative reduction is 41.0%. However, high-order resonant and anti-resonant responses have relatively small effects. Figure 5 shows the amplitude–frequency relations on the point x ¯ = 0.4 of the beam with periodic supports and with non-periodic supports (i.e., left support damping cus1 = 0.12 and left support damping cus1 = 0.22) under the harmonic excitation. It can be observed that the first and second response peaks decrease with the increase in the support damping and the resonant frequencies have small variation. As the left support damping increases from 0.02 to 0.12, the first response peak is reduced from 0.115 (at ω = 1.14) to 0.055 (at ω = 1.12) and the relative reduction is 52.1%. Therefore, increasing the stiffness and damping of several supports (as support periodicity defect) can especially reduce the first and second resonant responses of the nonlinearly vibrational beam, which has potential application to nonlinear structural vibration control.
Particularly, it can be observed that the third resonant response and resonant frequency remain unchanged against the increasing support stiffness. The dynamic characteristics of the frequency response of the beam with the support periodicity defect (i.e., several support stiffness variations) will determine the outstanding variation of vibration response distribution and jump or the alteration of mode component orders. The alteration of mode component orders caused by varying support stiffness has been demonstrated for beams with periodic supports [61]. It will be further demonstrated in the following subsections for the nonlinear beam with non-periodic supports (or support periodicity defect).

6.2. Effects of Non-Periodic Support Stiffness Increase on Amplitude–Frequency Characteristics and Vibration Anti-Localization

Special dimensionless amplitude–frequency characteristics and low-order mode component jump of the nonlinearly vibrational beam with non-periodic supports (i.e., a support stiffness larger than the others) are explored by increasing a support stiffness. Figure 6 shows the amplitude–frequency relations on the point x ¯ = 0.4 of the beam with non-periodic supports under the harmonic excitation for different left support stiffnesses (i.e., left support stiffness kus1 = 10.1, kus1 = 20.1, kus1 = 30.1, kus1 = 45.1, and kus1 = 0.1 as periodic support for comparison). It is obtained that the second resonant frequency increases with the support stiffness and the increment is larger than that of the first resonant frequency. However, the third resonant frequency remains unchanged, and thus the second resonant frequency will reach the third resonant frequency as the support stiffness increases. The original third resonant response vanishes when the left support stiffness kus1 > 30. However, the original third resonant response reappears when the left support stiffness kus1 > 40, and the resonant response then becomes the second resonant response in accordance with the corresponding frequency order.
Figure 7 and Figure 8 show, in the case of the left support stiffness kus1 = 10.1, the vibration response distributions (absolute amplitude values) and the corresponding mode components (first three components relevant to φn) of the beam with non-periodic supports (i.e., the left support stiffness is larger than the others) under the harmonic excitation with a dimensionless frequency close to the second and third resonant frequencies (ω = 5.0 and 9.0), respectively. The vibration response distributions of the beam without periodic support were also given for comparison. The second and third resonant responses have dominant components of the second and third modes, respectively. The second resonant response at the support ( x ¯ = 1/3, 2/3) is large, but the third resonant response at the support ( x ¯ = 1/3, 2/3) is close to zero, as shown in Figure 7 and Figure 8. Figure 9 shows, in the case of the left support stiffness kus1 = 45.1, the vibration response distribution and the corresponding mode components of the beam with non-periodic supports under the harmonic excitation with a dimensionless frequency close to the second resonant frequency (ω = 8.9). The second resonant response has a dominant component of the third mode. Then, the original third mode component (for left support stiffness kus1 = 10.1) becomes the second mode component (for left support stiffness kus1 = 45.1), and the original second mode component becomes the third mode component. This is the mode component jump or alteration of mode component orders, which yields outstanding variation of the vibration response distribution due to the changing stiffness of a non-periodic support. The left support stiffness kus1 = 40 is called critical stiffness for the mode alteration as the support periodicity defect. The current second resonant response at the support ( x ¯ = 1/3, 2/3) is close to zero, which is called local weak coupling between adjacent spans.
By the comparison of the vibration response distributions in Figure 7 and Figure 9, the large variation of the spatial distribution of the second resonant response is obtained, particularly, at the supports as the non-periodic support stiffness increases to exceed the critical stiffness. The comparison of the amplitude–frequency relations in Figure 6 results in the large variation of the second resonant response versus frequency as the non-periodic support stiffness increases. Figure 10 shows the evolution of the vibration response distribution of the beam with increasing left support stiffness (kus1) for ω = 8.9 near the original third resonant frequency. It is observed that the third mode component changes from dominance to puniness and returns dominance in the response. The remarkable variation of frequency response and the mode component jump or alteration of mode component orders for nonlinear beams with varying non-periodic supports (as non-periodic support stiffness increases) are revealed for the first time. The vibration response distribution (for the second resonant response) changes from the non-dominant third mode component (response at the support is large) to the dominant third mode component (response at the support is close to zero); this phenomenon is called vibration anti-localization. The anti-localization has potential application, e.g., in the vibration control of nonlinearly vibrational structures with non-periodic supports.

6.3. Effects of Non-Periodic Support Stiffness Reduction on Amplitude–Frequency Characteristics and Vibration Localization

Dimensionless amplitude–frequency characteristics and the low-order mode component jump of the nonlinearly vibrational beam with non-periodic supports (there is a support with its stiffness smaller than the others) are explored here by considering changes in large support stiffness. Figure 11 shows the amplitude–frequency relations at the point x ¯ = 0.4 of the beam with periodic supports under the harmonic excitation for different large support stiffnesses (kus1 = kus2) (the support stiffness is, respectively, kus1 = 10, kus1 = 23, kus1 = 31, and kus1 = 40). It can be seen that, as the support stiffness increases, the original second resonant response vanishes when the support stiffnesses kus1 = kus2 > 24, and the original second and third resonant responses vanish when the support stiffnesses kus1 = kus2 > 31. However, the original third resonant response reappears when the support stiffnesses kus1 = kus2 > 38, but the resonant response becomes the first resonant response and the original first resonant response then becomes the second resonant response in terms of the frequency order.
The beam with periodic supports of large support stiffnesses kus1 = kus2 = 45 has local weak coupling where the first resonant response at the support is close to zero. Now consider a support stiffness reduction for the periodic beam, e.g., due to damage, which is the support periodicity defect. Figure 12 illustrates the amplitude–frequency relations at the point x ¯ = 0.4 of the beam with non-periodic supports under the harmonic excitation for different left support stiffnesses (the left support stiffness is, respectively, kus1 = 40, kus1 = 38, kus1 = 25, and kus1 = 45 in the periodic support case for comparison). It can be seen that as the support stiffness decreases, the second resonant response reappears when the left support stiffness kus1 < 40. The resonant response peak increases and the resonant frequency decreases with the left support stiffness. The original first resonant response (for the periodic beam) vanishes when the left support stiffness kus1 < 38. However, the original first resonant response reappears when the left support stiffness kus1 < 25, and the resonant response then becomes the second resonant response in terms of the frequency order.
Figure 13 shows, in the case of the periodic support stiffnesses kus1 = kus2 = 45, the vibration response distribution (absolute amplitude values) and the corresponding mode components (first three components relevant to φn) of the beam under the harmonic excitation with a dimensionless frequency close to the first resonant frequency (ω = 8.9). It is observed that the first resonant response has a dominant component of the third mode, and the resonant response at the support ( x ¯ = 1/3, 2/3) is close to zero (local weak coupling). Figure 14 shows, in the case of the non-periodic supports with left support stiffness kus1 = 25, the vibration response distribution and the corresponding mode components of the beam under the harmonic excitation with dimensionless frequency close to the first resonant frequency (ω = 8.0). It is observed that the first resonant response has dominant components of the first and second modes, and the resonant response at the left support ( x ¯ = 1/3) is large. Therefore, the original second and third mode components (with support stiffness kus1 = kus2 = 45) then become the first and second mode components (with left support stiffness kus1 = 25), and the original first mode component then becomes the third mode component. This is the mode component jump or alteration of mode component orders (support stiffness kus1 = 25 is called critical stiffness), which yields outstanding variation of the vibration response distribution due to the stiffness reduction of a non-periodic support.
By the comparison of the vibration response distributions in Figure 13 and Figure 14, the large variation in the spatial distribution of the first resonant response (particularly at the support) is obtained, and the comparison of the amplitude–frequency relations in Figure 12 results in the large variation of the first resonant response versus frequency as the non-periodic support stiffness decreases below the critical stiffness. Figure 15 shows the evolution of the vibration response distribution of the beam with the reducing left support stiffness (kus1) for ω = 8.9 near the original first resonant frequency. The observation similar to Figure 10 was obtained. The outstanding variation of the frequency response and the mode component jump or alteration of mode component orders for nonlinear beams with varying non-periodic supports (as non-periodic support stiffness decreases) are revealed for the first time. The vibration response distribution (for the first resonant response) changes from the dominant third mode component (response at the support is close to zero) to the non-dominant third mode component (response at the support is large); this is called vibration localization. The localization has a potential application, e.g., to support the damage detection of nonlinearly vibrational structures with non-periodic supports.

7. Conclusions

A response analysis method for nonlinearly vibrational beams with spatial distribution parameters and non-periodic supports (or support periodicity defect) was developed by combining perturbation, periodic and non-periodic separation, Galerkin expansion and harmonic balance methods, and applied to analyze the amplitude–frequency characteristics of non-periodic beams with multi-mode coupling vibration. The proposed method has the following main advantages: (i) it separates a non-periodic structure and its response into periodic and non-periodic (or periodicity defect) parts; (ii) it takes multiple vibration modes for the response analysis of nonlinear non-periodic structures; (iii) it considers the coupling effects of vibration modes owing to structural nonlinearity and parametric periodicity; (iv) it is applicable to nonlinear non-periodic structures with a high-parameter-varying wave in a wide frequency vibration; and (v) it is suitable for effective numerical computation. The proposed method can be extended for the dynamic response analysis of other nonlinear structures with periodic and non-periodic distribution parameters.
Nonlinear beams with non-periodic supports (resulting in non-periodic distribution parameters or periodicity defect) under harmonic excitation have been studied to demonstrate the application of the proposed method which has good convergence and accuracy with increasing expansion terms. Graphical results on the amplitude–frequency characteristics of a non-periodic beam with various support stiffness and damping demonstrate that: (i) a small deviation of periodic supports has a slight effect on the response, and the increasing stiffness and damping of several supports can remarkably reduce the first and second resonant responses which has a potential application to nonlinear structural vibration control; (ii) for the local strong coupling case, the second resonant frequency largely increases with non-periodic support stiffness, but the third resonant frequency with the response remains almost unchanged against increasing non-periodic support stiffness; (iii) for the local strong coupling case, the second resonant response and spatial distribution have remarkable variation as the non-periodic support stiffness increases to exceed the critical stiffness, and the mode component jump or alteration of mode component orders gives rise to the outstanding variation or vibration anti-localization, which has a potential application to nonlinear structural vibration control; (iv) for the local weak coupling case, the first resonant response increases and the resonant frequency largely decreases with non-periodic support stiffness when it is smaller than the critical stiffness, and there is the mode component jump or alteration of mode component orders; and (v) for the local weak coupling case, the first resonant response and spatial distribution have remarkable variation as the non-periodic support stiffness decreases below the critical stiffness, and the mode component jump yields the vibration localization, which has a potential application to support damage detection of nonlinear structures with non-periodic supports. Some salient dynamic characteristics including outstanding frequency response variation and the mode component jump for nonlinearly vibrational beams with non-periodic supports (or support periodicity defect) were discovered for the first time.

Author Contributions

Conceptualization, Z.-G.Y. and Y.-Q.N.; methodology, Z.-G.Y. and Y.-Q.N.; software, Z.-G.Y.; validation, Z.-G.Y.; writing—original draft preparation, Z.-G.Y.; writing—review and editing, Y.-Q.N. and Z.-G.Y.; project administration, Z.-G.Y. and Y.-Q.N.; funding acquisition, Z.-G.Y. and Y.-Q.N. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (grant number 12072312), the Research Grants Council of the Hong Kong Special Administrative Region (grant number R-5020-18), and the Innovation and Technology Commission of the Hong Kong Special Administrative Region (grant number K-BBY1).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

This research was supported by the National Natural Science Foundation of China (grant number 12072312), the Research Grants Council of the Hong Kong Special Administrative Region (grant number R-5020-18), and the Innovation and Technology Commission of the Hong Kong Special Administrative Region to the Hong Kong Branch of the National Rail Transit Electrification and Automation Engineering Technology Research Centre (grant number K-BBY1). This support is gratefully acknowledged.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

B 1 j = 1 ε j 1 2 w j t 2 + B 2 j = 1 ε j 1 w j t + k = 1 3 B k + 2 j = 1 ε j 1 ( 5 k ) w j x ¯ ( 5 k ) + ε B 10 ( j = 1 ε j 1 2 w j x ¯ 2 ) ( j = 1 ε j 1 w j x ¯ ) 2 + ε B 11 ( j = 1 ε j 1 w j x ¯ ) 3 + i = 1 N p ( D 1 i j = 1 ε j 1 w j t + D 2 i j = 1 ε j 1 w j ) δ ( x ¯ x ¯ i ) = F ( x ¯ , t )
M p 1 = [ M p 1 , 11 M p 1 , 12 M p 1 , 1 N m M p 1 , 21 M p 1 , 22 M p 1 , 2 N m M p 1 , N m 1 M p 1 , N m 2 M p 1 , N m N m ] C p 1 = [ C p 1 , 11 C p 1 , 12 C p 1 , 1 N m C p 1 , 21 C p 1 , 22 C p 1 , 2 N m C p 1 , N m 1 C p 1 , N m 2 C p 1 , N m N m ] K p 1 = [ K p 1 , 11 K p 1 , 12 K p 1 , 1 N m K p 1 , 21 K p 1 , 22 K p 1 , 2 N m K p 1 , N m 1 K p 1 , N m 2 K p 1 , N m N m ] F p j , m = [ F p j , 1 F p j , 2 F p j , N m ] T M p 1 , m n = 0 1 ϕ m ( y ) B p 1 ( y ) ϕ n ( y ) d y C p 1 , m n = 0 1 ϕ m ( y ) B p 2 ( y ) ϕ n ( y ) d y + i = 1 N p ϕ m ( y p i ) D p 1 ϕ n ( y p i ) K p 1 . m n = 0 1 ϕ m ( y ) k = 1 3 B p , k + 2 ( y ) ( 5 k ) ϕ n ( y ) y ( 5 k ) d y + i = 1 N p ϕ m ( y p i ) D p 2 ϕ n ( y p i ) F p 1 , m = 0 1 ϕ m ( y ) F p ( y , t ) d y F p 2 , m = 0 1 ϕ m ( y ) B p 10 ( y ) [ 2 Φ ( y ) y 2 Q p 1 ] [ Φ ( y ) y Q p 1 ] 2 d y 0 1 ϕ m ( y ) B p 11 ( y ) [ Φ ( y ) y Q p 1 ] 3 d y F p 3 , m = 0 1 ϕ m ( y ) B p 10 ( y ) { [ 2 Φ ( y ) y 2 Q p 2 ] [ Φ ( y ) y Q p 1 ] 2 + 2 [ Φ ( y ) y Q p 2 ] [ 2 Φ ( y ) y 2 Q p 1 ] [ Φ ( y ) y Q p 1 ] } d y 3 0 1 ϕ m ( y ) B p 11 ( y ) [ Φ ( y ) y Q p 2 ] [ Φ ( y ) y Q p 1 ] 2 d y y p i = x ¯ p i m ,   n = 1 ,   2 ,   ,   N m ;   j = 1 ,   2 ,   3
M 1 = [ M 1 , 11 M 1 , 12 M 1 , 1 N m M 1 , 21 M 1 , 22 M 1 , 2 N m M 1 , N m 1 M 1 , N m 2 M 1 , N m N m ] C 1 = [ C 1 , 11 C 1 , 12 C 1 , 1 N m C 1 , 21 C 1 , 22 C 1 , 2 N m C 1 , N m 1 C 1 , N m 2 C 1 , N m N m ] K 1 = [ K 1 , 11 K 1 , 12 K 1 , 1 N m K 1 , 21 K 1 , 22 K 1 , 2 N m K 1 , N m 1 K 1 , N m 2 K 1 , N m N m ] F j , m = [ F j , 1 F j , 2 F j , N m ] T M 1 , m n = 0 1 ϕ m ( y ) B 1 ( y ) ϕ n ( y ) d y C 1 , m n = 0 1 ϕ m ( y ) B 2 ( y ) ϕ n ( y ) d y + i = 1 N p ϕ m ( y i ) D 1 i ϕ n ( y i ) K 1 , m n = 0 1 ϕ m ( y ) k = 1 3 B k + 2 ( y ) ( 5 k ) ϕ n ( y ) y ( 5 k ) d y + i = 1 N p ϕ m ( y i ) D 2 i ϕ n ( y i ) F 1 , m = 0 1 ϕ m ( y ) ( Δ F Δ B 1 B p 1 1 F p ) d y + 0 1 ϕ m ( y ) ( Δ B 1 B p 1 1 B p 2 Δ B 2 ) [ Φ ( y ) Q p 1 t ] d y + 0 1 ϕ m ( y ) k = 1 3 ( Δ B 1 B p 1 1 B p , k + 2 Δ B k + 2 ) [ ( 5 k ) Φ ( y ) y ( 5 k ) Q p 1 ] d y + i = 1 N p { ϕ m ( y p i ) ( 1 + Δ B 1 B p 1 1 ) [ D p 1 Φ ( y p i ) Q p 1 t + D p 2 Φ ( y p i ) Q p 1 ] ϕ m ( y i ) [ D 1 i Φ ( y i ) Q p 1 t + D 2 i Φ ( y i ) Q p 1 ] } F 2 , m = 0 1 ϕ m ( y ) ( Δ B 1 B p 1 1 B p 2 Δ B 2 ) [ Φ ( y ) Q p 2 t ] d y + 0 1 ϕ m ( y ) k = 1 3 ( Δ B 1 B p 1 1 B p , k + 2 Δ B k + 2 ) [ ( 5 k ) Φ ( y ) y ( 5 k ) Q p 2 ] d y + i = 1 N p { ϕ m ( y p i ) ( 1 + Δ B 1 B p 1 1 ) [ D p 1 Φ ( y p i ) Q p 2 t + D p 2 Φ ( y p i ) Q p 2 ] ϕ m ( y i ) [ D 1 i Φ ( y i ) Q p 2 t + D 2 i Φ ( y i ) Q p 2 ] } 0 1 ϕ m ( y ) { B 10 ( y ) [ 2 Φ ( y ) y 2 Q 1 ] [ Φ ( y ) y Q 1 ] 2 ( 1 + Δ B 1 B p 1 1 ) B p 10 ( y ) [ 2 Φ ( y ) y 2 Q p 1 ] [ Φ ( y ) y Q p 1 ] 2 } d y 0 1 ϕ m ( y ) { B 11 ( y ) [ Φ ( y ) y Q 1 ] 3 ( 1 + Δ B 1 B p 1 1 ) B p 11 ( y ) [ Φ ( y ) y Q p 1 ] 3 } d y F 3 , m = 0 1 ϕ m ( y ) ( Δ B 1 B p 1 1 B p 2 Δ B 2 ) [ Φ ( y ) Q p 3 t ] d y + 0 1 ϕ m ( y ) k = 1 3 ( Δ B 1 B p 1 1 B p , k + 2 Δ B k + 2 ) [ ( 5 k ) Φ ( y ) y ( 5 k ) Q p 3 ] d y + i = 1 N p { ϕ m ( y p i ) ( 1 + Δ B 1 B p 1 1 ) [ D p 1 Φ ( y p i ) Q p 3 t + D p 2 Φ ( y p i ) Q p 3 ] ϕ m ( y i ) [ D 1 i Φ ( y i ) Q p 3 t + D 2 i Φ ( y i ) Q p 3 ] } 0 1 ϕ m ( y ) B 10 ( y ) { [ 2 Φ ( y ) y 2 Q 2 ] [ Φ ( y ) y Q 1 ] 2 + 2 [ Φ ( y ) y Q 2 ] [ 2 Φ ( y ) y 2 Q 1 ] [ Φ ( y ) y Q 1 ] } d y + 0 1 ϕ m ( y ) ( 1 + Δ B 1 B p 1 1 ) B p 10 ( y ) { [ 2 Φ ( y ) y 2 Q p 2 ] [ Φ ( y ) y Q p 1 ] 2 + 2 [ Φ ( y ) y Q p 2 ] [ 2 Φ ( y ) y 2 Q p 1 ] [ Φ ( y ) y Q p 1 ] } d y 3 0 1 ϕ m ( y ) { B 11 ( y ) [ Φ ( y ) y Q 2 ] [ Φ ( y ) y Q 1 ] 2 ( 1 + Δ B 1 B p 1 1 ) B p 11 ( y ) [ Φ ( y ) y Q p 2 ] [ Φ ( y ) y Q p 1 ] 2 } d y y p i = x ¯ p i ,   y i = x ¯ i m ,   n = 1 ,   2 ,   ,   N m ;   j = 1 ,   2 ,   3

References

  1. Mead, D.J. Wave propagation in continuous periodic structures: Research contributions from Southampton, 1964–1995. J. Sound Vib. 1996, 190, 495–524. [Google Scholar] [CrossRef]
  2. Hussein, M.I.; Leamy, M.J.; Ruzzene, M. Dynamics of phononic materials and structures: Historical origins, recent progress and future outlook. ASME Appl. Mech. Rev. 2014, 66, 040802. [Google Scholar] [CrossRef]
  3. Chen, R.; Wu, T. Dynamic characteristics of a periodic rib-skin structure. J. Vib. Control 2016, 22, 662–677. [Google Scholar] [CrossRef]
  4. Ying, Z.G.; Ni, Y.Q. Dynamic characteristics of infinite-length and finite-length rods with high-wave-number periodic parameters. J. Vib. Control 2018, 24, 2344–2358. [Google Scholar] [CrossRef]
  5. Pierre, C. Mode localization and eigenvalue loci veering phenomena in disordered structures. J. Sound Vib. 1988, 126, 485–502. [Google Scholar] [CrossRef] [Green Version]
  6. Nayfeh, A.H.; Hawwa, M.A. Use of mode localization in passive control of structural buckling. AIAA J. 1994, 32, 2131–2133. [Google Scholar] [CrossRef]
  7. Hunt, G.W.; Wadee, M.A. Localization and mode interaction in sandwich structures. Proc. R. Soc. Lond. A 1998, 454, 1197–1216. [Google Scholar] [CrossRef]
  8. Bendisen, O.O. Localization phenomena in structural dynamics. Chaos Solitons Fractals 2000, 11, 1621–1660. [Google Scholar] [CrossRef]
  9. Luongo, A. Mode localization in dynamics and buckling of linear imperfect continuous structures. Nonlinear Dyn. 2001, 25, 133–156. [Google Scholar] [CrossRef]
  10. Ding, L.; Zhu, H.P.; Luo, H.; Yin, T. Flexural wave propagation and localization in periodic jointed tunnels subjected to moving loads. J. Vib. Control 2016, 22, 2788–2804. [Google Scholar] [CrossRef]
  11. Cai, C.W.; Cheung, Y.K.; Chan, H.C. Mode localization phenomena in nearly periodic systems. ASME J. Appl. Mech. 1995, 62, 141–149. [Google Scholar] [CrossRef]
  12. Elishakoff, I.; Li, Y.W.; Starnes, J.H. Buckling mode localization in elastic plates due to misplacement in the stiffener location. Chaos Solitons Fractals 1995, 5, 1517–1531. [Google Scholar] [CrossRef]
  13. Xie, W.C. Buckling mode localization in rib-stiffened plates with randomly misplaced stiffeners. Comput. Struct. 1998, 67, 175–189. [Google Scholar] [CrossRef]
  14. Ying, Z.G.; Ni, Y.Q.; Kang, L. Mode localization characteristics of damaged quasiperiodically supported beam structures with local weak coupling. Struct. Control Health Monit. 2019, 26, e2351. [Google Scholar] [CrossRef]
  15. Bouzit, D.; Pierre, C. Vibration confinement phenomena in disordered, mono-coupled, multi-span beams. ASME J. Vib. Acoust. 1992, 114, 521–530. [Google Scholar] [CrossRef]
  16. Cai, G.Q.; Lin, Y.K. Statistics distribution of frequency response in disordered periodic structures. AIAA J. 1992, 30, 1400–1407. [Google Scholar] [CrossRef]
  17. Bardell, N.S.; Langley, R.S.; Dunsdon, J.M.; Klein, T. The effect of period asymmetry on wave propagation in periodic beams. J. Sound Vib. 1996, 197, 427–445. [Google Scholar] [CrossRef]
  18. Yan, Z.Z.; Zhang, C.; Wang, Y.S. Attenuation and localization of bending waves in a periodic/disordered fourfold composite beam. J. Sound Vib. 2009, 327, 109–120. [Google Scholar] [CrossRef]
  19. Ruzzene, M.; Baz, A. Control of wave propagation in periodic composite rods using shape memory inserts. ASME J. Vib. Acoust. 2000, 122, 151–159. [Google Scholar] [CrossRef]
  20. Shelley, F.J.; Clark, W.M. Active mode localization in distributed parameter systems with consideration of limited actuator placement, part 1: Theory. ASME J. Vib. Acoust. 2000, 122, 160–164. [Google Scholar] [CrossRef]
  21. Romeo, F.; Luongo, A. Vibration reduction in piecewise bi-coupled periodic structures. J. Sound Vib. 2003, 268, 601–615. [Google Scholar] [CrossRef] [Green Version]
  22. Wu, T.Y.; Wang, K.W. Periodic isolator design enhancement via vibration confinement through eigenvector assignment and piezoelectric circuitry. J. Vib. Control 2007, 13, 989–1006. [Google Scholar] [CrossRef]
  23. Song, Y.; Wen, J.; Yu, D.; Wen, X. Suppression of vibration and noise radiation in a flexible floating raft system using periodic structures. J. Vib. Control 2015, 21, 217–228. [Google Scholar] [CrossRef]
  24. Chen, J.S.; Tsai, S.M. Sandwich structures with periodic assemblies on elastic foundation under moving loads. J. Vib. Control 2016, 22, 2519–2529. [Google Scholar] [CrossRef]
  25. Lallart, M.; Yan, L.; Richard, C.; Guyomar, D. Damping of periodic bending structures featuring nonlinearly interfaced piezoelectric elements. J. Vib. Control 2016, 22, 3930–3941. [Google Scholar] [CrossRef]
  26. Wang, Y.Z.; Li, F.M.; Wang, Y.S. Influences of active control on elastic wave propagation in a weakly nonlinear phononic crystal with a monoatomic lattice chain. Int. J. Mech. Sci. 2016, 106, 357–362. [Google Scholar] [CrossRef]
  27. Harne, R.L.; Urbanek, D.C. Enhancing broadband vibration energy suppression using local buckling modes in constrained metamaterials. ASME J. Vib. Acoust. 2017, 139, 061004. [Google Scholar] [CrossRef] [Green Version]
  28. Lee, S.Y.; Ke, H.Y.; Kao, M.J. Flexural waves in a periodic beam. ASME J. Appl. Mech. 1990, 57, 779–783. [Google Scholar] [CrossRef]
  29. Lee, S.Y.; Ke, H.Y. Flexural wave propagation in an elastic beam with periodic structure. ASME J. Appl. Mech. 1992, 59, S189–S196. [Google Scholar] [CrossRef]
  30. Cheung, Y.K.; Au, F.T.K.; Zheng, D.Y.; Cheng, Y.S. Vibration of multi-span non-uniform bridges under moving vehicles and trains by using modified beam vibration functions. J. Sound Vib. 1999, 228, 611–628. [Google Scholar] [CrossRef]
  31. Au, F.T.K.; Zheng, D.Y.; Cheung, Y.K. Vibration and stability of non-uniform beams with abrupt changes of cross-section by using C1 modified beam vibration functions. Appl. Math. Model. 1999, 23, 19–34. [Google Scholar] [CrossRef]
  32. Cottrell, J.A.; Realim, A.; Bazilevs, Y.; Hughes, T.J.R. Isogeometric analysis of structural vibrations. Comput. Methods Appl. Mech. Eng. 2006, 195, 5257–5296. [Google Scholar] [CrossRef] [Green Version]
  33. Sun, Y.; Liu, S.; Rao, Z.; Li, Y.; Yang, J. Thermodynamic response of beams on Winkler foundation irradiated by moving laser pulses. Symmetry 2018, 10, 328. [Google Scholar] [CrossRef] [Green Version]
  34. Shariati, A.; Jung, D.W.; Mohammad-Sedighi, H.; Zur, K.K.; Habibi, M.; Safa, M. Stability and dynamics of viscoelastic moving Rayleigh beams with an asymmetrical distribution of material parameters. Symmetry 2020, 12, 586. [Google Scholar] [CrossRef] [Green Version]
  35. Hawwa, M.A. Reflection of flexural waves in geometrically periodic beams. J. Sound Vib. 1997, 199, 453–461. [Google Scholar] [CrossRef]
  36. McIntyre, J.S.; Rasmussen, M.L.; Bert, C.W.; Kline, R.A. Resonance in fiber-reinforced composite materials with sinusoidal stiffness properties. Wave Motion 1999, 30, 97–123. [Google Scholar] [CrossRef]
  37. Wei, J.; Petyt, M. A method of analyzing finite periodic structures, part 2: Comparison with infinite periodic structure theory. J. Sound Vib. 1997, 202, 571–583. [Google Scholar] [CrossRef]
  38. Hvatov, A.; Sorokin, S. Free vibration of finite periodic structures in pass- and stop-bands of the counterpart infinite waveguides. J. Sound Vib. 2015, 347, 200–217. [Google Scholar] [CrossRef]
  39. Junyi, L.; Balint, D.S. An inverse method to determine the dispersion curves of periodic structures based on wave superposition. J. Sound Vib. 2015, 350, 41–72. [Google Scholar] [CrossRef] [Green Version]
  40. Domadiya, P.G.; Manconi, E.; Vanali, M.; Andersen, L.V.; Ricci, A. Numerical and experimental investigation of stop-bands in finite and infinite periodic one-dimensional structures. J. Vib. Control 2016, 22, 920–931. [Google Scholar] [CrossRef]
  41. Marathe, A.; Chatterjee, A. Wave attenuation in nonlinear periodic structures using harmonic balance and multiple scales. J. Sound Vib. 2005, 289, 871–888. [Google Scholar] [CrossRef]
  42. Lazarov, B.S.; Jensen, J.S. Low-frequency band gaps in chains with attached non-linear oscillators. Int. J. Non-Linear Mech. 2007, 42, 1186–1193. [Google Scholar] [CrossRef] [Green Version]
  43. Narisetti, R.K.; Leamy, M.J.; Ruzzene, M. A perturbation approach for predicting wave propagation in one-dimensional nonlinear periodic structures. ASME J. Vib. Acoust. 2010, 132, 031001. [Google Scholar] [CrossRef]
  44. Narisetti, R.K.; Ruzzene, M.; Leamy, M.J. Study of wave propagation in strongly nonlinear periodic lattices using a harmonic balance approach. Wave Motion 2012, 49, 394–410. [Google Scholar] [CrossRef]
  45. Clementi, F.; Demeio, L.; Mazzilli, C.E.N.; Lenci, S. Nonlinear vibrations of non-uniform beams by the MTS asymptotic expansion method. Contin. Mech. Thermodyn. 2015, 27, 703–717. [Google Scholar] [CrossRef]
  46. Yousefzadeh, B.; Phani, A.S. Supratransmission in a disordered nonlinear periodic structure. J. Sound Vib. 2016, 380, 242–266. [Google Scholar] [CrossRef] [Green Version]
  47. Bitar, D.; Kacem, N.; Bouhaddi, N. Investigation of modal interactions and their effects on the nonlinear dynamics of a periodic coupled pendulums chain. Int. J. Mech. Sci. 2017, 127, 130–141. [Google Scholar] [CrossRef] [Green Version]
  48. Vakakis, A.F. Dynamics of a nonlinear periodic structure with cycle symmetry. Acta Mech. 1992, 95, 197–226. [Google Scholar] [CrossRef]
  49. Vakakis, A.F.; Cetinkaya, C. Mode localization in a class of multidegree-of-freedom nonlinear systems with cyclic symmetry. SIAM J. Appl. Math. 1993, 53, 265–282. [Google Scholar] [CrossRef] [Green Version]
  50. Freno, B.A.; Cizmas, P.G.A. A computationally efficient non-linear beam model. Int. J. Non-Linear Mech. 2011, 46, 854–869. [Google Scholar] [CrossRef]
  51. Weeger, O.; Wever, U.; Simeon, B. Isogeometric analysis of nonlinear Euler-Bernoulli beam vibrations. Nonlinear Dyn. 2013, 72, 813–835. [Google Scholar] [CrossRef] [Green Version]
  52. Jang, T.S. A general method for analyzing moderately large deflections of a non-uniform beam: An infinite Bernoulli–Euler–von Karman beam on a nonlinear elastic foundation. Acta Mech. 2014, 225, 1967–1984. [Google Scholar] [CrossRef]
  53. Masjedi, P.K.; Maheri, A. Chebyshev collocation method for the free vibration analysis of geometrically exact beams with fully intrinsic formulation. Eur. J. Mech. A/Solids 2017, 66, 329–340. [Google Scholar] [CrossRef] [Green Version]
  54. Ying, Z.G.; Ni, Y.Q. Vibrational amplitude frequency characteristics analysis of a controlled nonlinear meso-scale beam. Actuators 2021, 10, 180. [Google Scholar] [CrossRef]
  55. Pourasghar, A.; Chen, Z. Nonlinear vibration and modal analysis of FG nanocomposite sandwich beams reinforced by aggregated CNTs. Polym. Eng. Sci. 2019, 59, 1362–1370. [Google Scholar] [CrossRef]
  56. Pourasghar, A.; Chen, Z. Effect of hyperbolic heat conduction on the linear and nonlinear vibration of CNT reinforced size-dependent functionally graded microbeams. Int. J. Eng. Sci. 2019, 137, 57–72. [Google Scholar] [CrossRef]
  57. Baleanu, D.; Sajjadi, S.S.; Jajarmi, A.; Defterli, O. On a nonlinear dynamical system with both chaotic and non-chaotic behaviors: A new fractional analysis and control. Adv. Differ. Equ. 2021, 2021, 234. [Google Scholar] [CrossRef]
  58. Baleanu, D.; Sajjadi, S.S.; Asad, J.H.; Jajarmi, A.; Estiri, E. Hyperchaotic behaviors, optimal control, and synchronization of a nonautonomous cardiac conduction system. Adv. Differ. Equ. 2021, 2021, 157. [Google Scholar] [CrossRef]
  59. Domagalski, L.; Jedrysiak, J. Geometrically nonlinear vibrations of slender meso-periodic beams. the tolerance modeling approach. Compos. Struct. 2016, 136, 270–277. [Google Scholar] [CrossRef]
  60. Domagalski, L.; Swiatek, M.; Jedrysiak, J. An analytical-numerical approach to vibration analysis of periodic Timoshenko beams. Compos. Struct. 2019, 211, 490–501. [Google Scholar] [CrossRef]
  61. Ying, Z.G.; Ni, Y.Q. A multimode perturbation method for frequency response analysis of nonlinearly vibrational beams with periodic parameters. J. Vib. Control 2020, 26, 1260–1272. [Google Scholar] [CrossRef]
  62. Spencer, B.F.; Nagarajaiah, S. State of the art of structural control. ASCE J. Struct. Eng. 2003, 129, 845–856. [Google Scholar] [CrossRef]
  63. Casciati, F.; Rodellar, J.; Yildirim, U. Active and semi-active control of structures—Theory and applications: A review of recent advances. J. Intell. Mater. Syst. Struct. 2012, 23, 1181–1195. [Google Scholar] [CrossRef]
  64. Carden, E.P.; Fanning, P. Vibration based condition monitoring: A review. Struct. Health Monit. 2004, 3, 355–377. [Google Scholar] [CrossRef]
  65. Fan, W.; Qiao, P.Z. Vibration-based damage identification methods: A review and comparative study. Struct. Health Monit. 2011, 10, 83–111. [Google Scholar] [CrossRef]
  66. Yang, Y.B.; Yang, J.P. State-of-the-art review on modal identification and damage detection of bridges by moving test vehicles. Int. J. Struct. Stab. Dyn. 2017, 18, 1850025. [Google Scholar] [CrossRef]
Figure 1. Beam with spatial distribution parameters and non-periodic middle supports under external excitation.
Figure 1. Beam with spatial distribution parameters and non-periodic middle supports under external excitation.
Symmetry 13 02234 g001
Figure 2. Response amplitudes W of the beam with periodic supports and with non-periodic supports under excitation frequency ω (black dotted line: beam with periodic supports; blue dashed line: non-periodic beam with left support deviation to the right Δ x ¯ s1 = 0.1/3; and red solid line: non-periodic beam with left support deviation to the left Δ x ¯ s1 = −0.1/3).
Figure 2. Response amplitudes W of the beam with periodic supports and with non-periodic supports under excitation frequency ω (black dotted line: beam with periodic supports; blue dashed line: non-periodic beam with left support deviation to the right Δ x ¯ s1 = 0.1/3; and red solid line: non-periodic beam with left support deviation to the left Δ x ¯ s1 = −0.1/3).
Symmetry 13 02234 g002
Figure 3. Relative variation of the response amplitudes W of non-periodic beam with left support deviation to the right (Δ x ¯ s1 = 0.1/3) to the beam with periodic supports under excitation frequency ω.
Figure 3. Relative variation of the response amplitudes W of non-periodic beam with left support deviation to the right (Δ x ¯ s1 = 0.1/3) to the beam with periodic supports under excitation frequency ω.
Symmetry 13 02234 g003
Figure 4. Response amplitudes W of the beam with periodic supports and with non-periodic supports (different left support stiffnesses kus1) under excitation frequency ω (black dotted line: beam with periodic supports; green dashed line: non-periodic beam with kus1 = 0.6; and blue dash-dotted line: non-periodic beam with kus1 = 1.1).
Figure 4. Response amplitudes W of the beam with periodic supports and with non-periodic supports (different left support stiffnesses kus1) under excitation frequency ω (black dotted line: beam with periodic supports; green dashed line: non-periodic beam with kus1 = 0.6; and blue dash-dotted line: non-periodic beam with kus1 = 1.1).
Symmetry 13 02234 g004
Figure 5. Response amplitudes W of the beam with periodic supports and with non-periodic supports (different left support damping cus1) under excitation frequency ω (black dotted line: beam with periodic supports; green dashed line: non-periodic beam with cus1 = 0.12; and blue dash-dotted line: non-periodic beam with cus1 = 0.22).
Figure 5. Response amplitudes W of the beam with periodic supports and with non-periodic supports (different left support damping cus1) under excitation frequency ω (black dotted line: beam with periodic supports; green dashed line: non-periodic beam with cus1 = 0.12; and blue dash-dotted line: non-periodic beam with cus1 = 0.22).
Symmetry 13 02234 g005
Figure 6. Response amplitudes W of the beam with non-periodic supports under excitation frequency ω for different left support stiffnesses kus1 (black dotted line: kus1 = 10.1; green dashed line: kus1 = 20.1; blue dash-dotted line: kus1 = 30.1; red solid line: kus1 = 45.1; and grey dash-dot-dashed line: kus1 = 0.1 as periodic support for comparison).
Figure 6. Response amplitudes W of the beam with non-periodic supports under excitation frequency ω for different left support stiffnesses kus1 (black dotted line: kus1 = 10.1; green dashed line: kus1 = 20.1; blue dash-dotted line: kus1 = 30.1; red solid line: kus1 = 45.1; and grey dash-dot-dashed line: kus1 = 0.1 as periodic support for comparison).
Symmetry 13 02234 g006
Figure 7. Vibration response distribution W(y,ω) (ω = 5.0) and mode components qfi(y) of the beam with non-periodic supports for left support stiffness kus1 = 10.1 (black solid line: response amplitude W; red dotted line: first mode component qf1; blue dashed line: second mode component qf2; green dash-dotted line: third mode component qf3; and grey dotted line: response amplitude without periodic support) (y = x ¯ ).
Figure 7. Vibration response distribution W(y,ω) (ω = 5.0) and mode components qfi(y) of the beam with non-periodic supports for left support stiffness kus1 = 10.1 (black solid line: response amplitude W; red dotted line: first mode component qf1; blue dashed line: second mode component qf2; green dash-dotted line: third mode component qf3; and grey dotted line: response amplitude without periodic support) (y = x ¯ ).
Symmetry 13 02234 g007
Figure 8. Vibration response distribution W(y,ω) (ω = 9.0) and mode components qfi(y) of the beam with non-periodic supports for left support stiffness kus1 = 10.1 (black solid line: response amplitude W; red dotted line: first mode component qf1; blue dashed line: second mode component qf2; green dash-dotted line: third mode component qf3; and grey dotted line: response amplitude without periodic support) (y = x ¯ ).
Figure 8. Vibration response distribution W(y,ω) (ω = 9.0) and mode components qfi(y) of the beam with non-periodic supports for left support stiffness kus1 = 10.1 (black solid line: response amplitude W; red dotted line: first mode component qf1; blue dashed line: second mode component qf2; green dash-dotted line: third mode component qf3; and grey dotted line: response amplitude without periodic support) (y = x ¯ ).
Symmetry 13 02234 g008
Figure 9. Vibration response distribution W(y,ω) (ω = 8.9) and mode components qfi(y) of the beam with non-periodic supports for left support stiffness kus1 = 45.1 (black solid line: response amplitude W; red dotted line: first mode component qf1; blue dashed line: second mode component qf2; and green dash-dotted line: third mode component qf3) (y = x ¯ ).
Figure 9. Vibration response distribution W(y,ω) (ω = 8.9) and mode components qfi(y) of the beam with non-periodic supports for left support stiffness kus1 = 45.1 (black solid line: response amplitude W; red dotted line: first mode component qf1; blue dashed line: second mode component qf2; and green dash-dotted line: third mode component qf3) (y = x ¯ ).
Symmetry 13 02234 g009
Figure 10. Vibration response distribution W(y,ω) (ω = 8.9) of the beam with non-periodic supports for varying left support stiffness kus1 (y = x ¯ ).
Figure 10. Vibration response distribution W(y,ω) (ω = 8.9) of the beam with non-periodic supports for varying left support stiffness kus1 (y = x ¯ ).
Symmetry 13 02234 g010
Figure 11. Response amplitudes W of the beam with periodic supports under excitation frequency ω for different support stiffnesses kus1 = kus2 (black dotted line: kus1 = 10; green dashed line: kus1 = 23; blue dash-dotted line: kus1 = 31; and red solid line: kus1 = 40).
Figure 11. Response amplitudes W of the beam with periodic supports under excitation frequency ω for different support stiffnesses kus1 = kus2 (black dotted line: kus1 = 10; green dashed line: kus1 = 23; blue dash-dotted line: kus1 = 31; and red solid line: kus1 = 40).
Symmetry 13 02234 g011
Figure 12. Response amplitudes W of the beam with non-periodic supports under excitation frequency ω for different left support stiffnesses kus1 (green dashed line: kus1 = 40; blue dash-dotted line: kus1 = 38; red solid line: kus1 = 25; and black dotted line: kus1 = 45 as periodic support for comparison).
Figure 12. Response amplitudes W of the beam with non-periodic supports under excitation frequency ω for different left support stiffnesses kus1 (green dashed line: kus1 = 40; blue dash-dotted line: kus1 = 38; red solid line: kus1 = 25; and black dotted line: kus1 = 45 as periodic support for comparison).
Symmetry 13 02234 g012
Figure 13. Vibration response distribution W(y,ω) (ω = 8.9) and mode components qfi(y) of the beam with periodic supports for support stiffness kus1 = kus2 = 45 (black solid line: response amplitude W; red dotted line: first mode component qf1; blue dashed line: second mode component qf2; and green dash-dotted line: third mode component qf3).
Figure 13. Vibration response distribution W(y,ω) (ω = 8.9) and mode components qfi(y) of the beam with periodic supports for support stiffness kus1 = kus2 = 45 (black solid line: response amplitude W; red dotted line: first mode component qf1; blue dashed line: second mode component qf2; and green dash-dotted line: third mode component qf3).
Symmetry 13 02234 g013
Figure 14. Vibration response distribution W(y,ω) (ω = 8.0) and mode components qfi(y) of the beam with non-periodic supports for left support stiffness kus1 = 25 (black solid line: response amplitude W; red dotted line: first mode component qf1; blue dashed line: second mode component qf2; and green dash-dotted line: third mode component qf3).
Figure 14. Vibration response distribution W(y,ω) (ω = 8.0) and mode components qfi(y) of the beam with non-periodic supports for left support stiffness kus1 = 25 (black solid line: response amplitude W; red dotted line: first mode component qf1; blue dashed line: second mode component qf2; and green dash-dotted line: third mode component qf3).
Symmetry 13 02234 g014
Figure 15. Vibration response distribution W(y,ω) (ω = 8.9) of the beam with non-periodic supports for varying left support stiffness kus1 (y = x ¯ ).
Figure 15. Vibration response distribution W(y,ω) (ω = 8.9) of the beam with non-periodic supports for varying left support stiffness kus1 (y = x ¯ ).
Symmetry 13 02234 g015
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ying, Z.-G.; Ni, Y.-Q. Vibration Localization and Anti-Localization of Nonlinear Multi-Support Beams with Support Periodicity Defect. Symmetry 2021, 13, 2234. https://doi.org/10.3390/sym13122234

AMA Style

Ying Z-G, Ni Y-Q. Vibration Localization and Anti-Localization of Nonlinear Multi-Support Beams with Support Periodicity Defect. Symmetry. 2021; 13(12):2234. https://doi.org/10.3390/sym13122234

Chicago/Turabian Style

Ying, Zu-Guang, and Yi-Qing Ni. 2021. "Vibration Localization and Anti-Localization of Nonlinear Multi-Support Beams with Support Periodicity Defect" Symmetry 13, no. 12: 2234. https://doi.org/10.3390/sym13122234

APA Style

Ying, Z. -G., & Ni, Y. -Q. (2021). Vibration Localization and Anti-Localization of Nonlinear Multi-Support Beams with Support Periodicity Defect. Symmetry, 13(12), 2234. https://doi.org/10.3390/sym13122234

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop