In this investigation,some different approaches are implemented for analyzing a generalized forced damped complex Duffing oscillator,including the hybrid homotopy perturbation method(H-HPM),which is sometimes called t...In this investigation,some different approaches are implemented for analyzing a generalized forced damped complex Duffing oscillator,including the hybrid homotopy perturbation method(H-HPM),which is sometimes called the Krylov-Bogoliubov-Mitropolsky(KBM)method and the multiple scales method(MSM).All mentioned methods are applied to obtain some accurate and stable approximations to the proposed problem without decoupling the original problem.All obtained approximations are discussed graphically using different numerical values to the relevant parameters.Moreover,all obtained approximate solutions are compared with the 4thorder Runge-Kutta(RK4)numerical approximation.The maximum residual distance error(MRDE)is also estimated,in order to verify the high accuracy of the obtained analytic approximations.展开更多
We propose a novel computational framework that is capable of employing different time integration algorithms and different space discretized methods such as the Finite Element Method,particle methods,and other spatia...We propose a novel computational framework that is capable of employing different time integration algorithms and different space discretized methods such as the Finite Element Method,particle methods,and other spatial methods on a single body sub-dividedintomultiple subdomains.This is in conjunctionwithimplementing thewell known Generalized Single Step Single Solve(GS4)family of algorithms which encompass the entire scope of Linear Multistep algorithms that have been developed over the past 50 years or so and are second order accurate into the Differential Algebraic Equation framework.In the current state of technology,the coupling of altogether different time integration algorithms has been limited to the same family of algorithms such as theNewmarkmethods and the coupling of different algorithms usually has resulted in reduced accuracy in one or more variables including the Lagrange multiplier.However,the robustness and versatility of the GS4 with its ability to accurately account for the numerical shifts in various time schemes it encompasses,overcomes such barriers and allows a wide variety of arbitrary implicit-implicit,implicit-explicit,and explicit-explicit pairing of the various time schemes while maintaining the second order accuracy in time for not only all primary variables such as displacement,velocity and acceleration but also the Lagrange multipliers used for coupling the subdomains.By selecting an appropriate spatialmethod and time scheme on the area with localized phenomena contrary to utilizing a single process on the entire body,the proposed work has the potential to better capture the physics of a given simulation.The method is validated by solving 2D problems for the linear second order systems with various combination of spatial methods and time schemes with great flexibility.The accuracy and efficacy of the present work have not yet been seen in the current field,and it has shown significant promise in its capabilities and effectiveness for general linear dynamics through numerical examples.展开更多
The main objective of this paper is to investigate the influence of inertia of nonlinear springs on the dispersion behavior of discrete monoatomic chains with lumped and distributed masses.The developed model can repr...The main objective of this paper is to investigate the influence of inertia of nonlinear springs on the dispersion behavior of discrete monoatomic chains with lumped and distributed masses.The developed model can represent the wave propagation problem in a non-homogeneous material consisting of heavy inclusions embedded in a matrix.The inclusions are idealized by lumped masses,and the matrix between adjacent inclusions is modeled by a nonlinear spring with distributed masses.Additionally,the model is capable of depicting the wave propagation in bi-material bars,wherein the first material is represented by a rigid particle and the second one is represented by a nonlinear spring with distributed masses.The discrete model of the nonlinear monoatomic chain with lumped and distributed masses is first considered,and a closed-form expression of the dispersion relation is obtained by the second-order Lindstedt-Poincare method(LPM).Next,a continuum model for the nonlinear monoatomic chain is derived directly from its discrete lattice model by a suitable continualization technique.The subsequent use of the second-order method of multiple scales(MMS)facilitates the derivation of the corresponding nonlinear dispersion relation in a closed form.The novelties of the present study consist of(i)considering the inertia of nonlinear springs on the dispersion behavior of the discrete mass-spring chains;(ii)developing the second-order LPM for the wave propagation in the discrete chains;and(iii)deriving a continuum model for the nonlinear monoatomic chains with lumped and distributed masses.Finally,a parametric study is conducted to examine the effects of the design parameters and the distributed spring mass on the nonlinear dispersion relations and phase velocities obtained from both the discrete and continuum models.These parameters include the ratio of the spring mass to the lumped mass,the nonlinear stiffness coefficient of the spring,and the wave amplitude.展开更多
A perturbation algorithm Multiple Scales Modified Lindstedt–Poincare(MSMLP),combination of method of Multiple Scales and modified Lindstedt–Poincare is proposed for the solution of Quintic Duffing equation which com...A perturbation algorithm Multiple Scales Modified Lindstedt–Poincare(MSMLP),combination of method of Multiple Scales and modified Lindstedt–Poincare is proposed for the solution of Quintic Duffing equation which combines the advantages of both the methods.Solution obtained by the MSMLP method is compared with the Multiple Scales method and accurate closed form approximate solution of the Quintic Duffing equation.The proposed method produces better results for a wide range of amplitude values of oscillations and strong nonlinearities.Numerical simulation has been performed in MATHEMATICA 7.0.展开更多
In this article, the nonlinear dynamic responses of sandwich functionally graded(FG) porous cylindrical shell embedded in elastic media are investigated. The shell studied here consists of three layers, of which the o...In this article, the nonlinear dynamic responses of sandwich functionally graded(FG) porous cylindrical shell embedded in elastic media are investigated. The shell studied here consists of three layers, of which the outer and inner skins are made of solid metal, while the core is FG porous metal foam. Partial differential equations are derived by utilizing the improved Donnell's nonlinear shell theory and Hamilton's principle. Afterwards, the Galerkin method is used to transform the governing equations into nonlinear ordinary differential equations, and an approximate analytical solution is obtained by using the multiple scales method. The effects of various system parameters,specifically, the radial load, core thickness, foam type, foam coefficient, structure damping,and Winkler-Pasternak foundation parameters on nonlinear internal resonance of the sandwich FG porous thin shells are evaluated.展开更多
In this paper,we used time delay feedback to minimize the vibrations of a hybrid Rayleigh–van der Pol–Duffing oscillator.This system is a one-degree-offreedom containing the cubic and fifth nonlinear terms and an ex...In this paper,we used time delay feedback to minimize the vibrations of a hybrid Rayleigh–van der Pol–Duffing oscillator.This system is a one-degree-offreedom containing the cubic and fifth nonlinear terms and an external force.We applied the multiple scales method to get the solution from first approximation.Graphically and numerically,we studied the system before and after adding time delay feedback at the primary resonance case(ffi!).We used MATLAB program to simulate the efficacy of different parameters and the time delay on the main system.展开更多
In this paper a modified L-P method and multiple scale method are used to solve sub-harmonic resonance solutions of strong and nonlinear resonance of general Van der Pol equation with parametric and external excitatio...In this paper a modified L-P method and multiple scale method are used to solve sub-harmonic resonance solutions of strong and nonlinear resonance of general Van der Pol equation with parametric and external excitations by parametric transformation. Bifurcation response equation and transition sets of sub-harmonic resonance with strong nonlinearity of general Van der Pol equation with parametric and external excitation are worked out.Besides, transition sets and bifurcation graphs are drawn to help to analysis the problems theoretically. Conclusions show that the transition sets of general and nonlinear Van der Pol equation with parametric and external excitations are more complex than those of general and nonlinear Van der Pol equation only with parametric excitation, which is helpful for the qualitative and quantitative reference for engineering and science applications.展开更多
An electromagnetic nonlinear energy harvester(NEH)based on a rotating system is proposed,of which the host system rotates at a constant speed and vibrates harmonically in the vertical direction.This kind of device exh...An electromagnetic nonlinear energy harvester(NEH)based on a rotating system is proposed,of which the host system rotates at a constant speed and vibrates harmonically in the vertical direction.This kind of device exhibits several resonant phenomena due to the combinations of the rotating and the vibration frequencies of the host system as well as the cubic nonlinearity of the NEH.The governing equation of motion for the NEH is derived,and the dynamic responses and output power are investigated with the multiple scale method under the 1:1 primary and 2:1 superharmonic resonant conditions.The effects of system parameters including the nondimensional external frequency,the rotating speed,and the nonlinear stiffness on the responses of free vibration for the system are studied.The results of the primary resonance show that the responses exhibit not only the resonant characteristics but also the nonlinear dynamic characteristics such as the saddle-node(SN)bifurcation.The coexistence of multiple solutions and the varying trends of responses are verified with the direct numerical simulation.Moreover,the effects of system parameters on the average output power are investigated.The results of the analyses on the two resonant conditions indicate that the large power can be harvested in two resonant frequency bands.The effect of resonance on the output power is dominant for the 2:1 superharmonic resonance.Moreover,the results also show that introducing the nonlinearity can increase the value of the output power in large frequency bands and induce the occurence of new frequency bands to harvest the large power.The efficiency of the harvested power could be improved by the combined effects of the resonance as well as the nonlinearity of the NEH device.Suitable parameter conditions could help optimize the power harvesting in design.展开更多
The nonlinear interactions of a microarch resonator with 3:1 internal resonance are studied.The microarch is subjected to a combination of direct current(DC)and alternating current(AC)electric voltages.Thin piezoelect...The nonlinear interactions of a microarch resonator with 3:1 internal resonance are studied.The microarch is subjected to a combination of direct current(DC)and alternating current(AC)electric voltages.Thin piezoelectric layers are thoroughly bonded on the top and bottom surfaces of the microarch.The piezoelectric actuation is not only used to modulate the stiffness and resonance frequency of the resonator but also to provide the suitable linear frequency ratio for the activation of the internal resonance.The size effect is incorporated by using the so-called modified strain gradient theory.The system is highly nonlinear due to the co-existence of the initial curvature,the mid-plane stretching resulting from clamped anchors,and the electrostatic excitation.The eigenvalue problem is solved to conduct a frequency analysis and identify the possible regions for activating the internal resonance.The effects of the piezoelectric actuation,the electric excitation,and the small-scale effect are investigated on the internal resonance.Exclusive nonlinear phenomena such as Hopf bifurcation and hysteresis are identified in the microarch response.It is shown that by applying appropriate piezoelectric actuation,one is able to activate microarch internal resonance regardless of the initial rise level of the microarch.It is also disclosed that among all the parameters,AC electric voltage has the greatest effect on the energy exchange between the interacting modes.The results can be used to design resonators and internal resonance based micro-electro-mechanical system(MEMS)energy harvesters.展开更多
An externally excited Duffing oscillator under feedback control is discussed and analyzed under the worst resonance case.Multiple time scales method is applied for this system to find analytic solution with the existe...An externally excited Duffing oscillator under feedback control is discussed and analyzed under the worst resonance case.Multiple time scales method is applied for this system to find analytic solution with the existence and nonexistence of the time delay on control loop.An appropriate stability analysis is also performed and appropriate choices for the feedback gains and the time delay are found in order to reduce the amplitude peak.Different response curves are involved to show and compare controller effects.In addition,analytic solutions are compared with numerical approximation solutions using Rung-Kutta method of fourth order.展开更多
Micro-electro-mechanical system(MEMS)gyroscopes are an important sort of inertial sensor for identifying parameters of spinning structures,such as the spinning speed and angular deviation,based on the Coriolis effect....Micro-electro-mechanical system(MEMS)gyroscopes are an important sort of inertial sensor for identifying parameters of spinning structures,such as the spinning speed and angular deviation,based on the Coriolis effect.In this paper,the nonlinear mechanism of MEMS vibratory ring gyroscopes is analyzed by applying a fully coupled nonlinear model,in which the gyroscopic coupling and geometrically and structurally nonlinear couplings are all taken into account.The coupled differential equations governing the drive and sense motions are established via the Lagrangian equations.Numerical simulation is conducted,and the key nonlinear components and energy transfer behaviors between the drive and sense modes are elucidated.It is revealed that the cubic rigidity nonlinearity is another significant factor leading to the coupling between the drive and sense modes other than the gyroscopic coupling.Perturbation analysis is also carried out by using the method of multiple scales.The nonlinear frequency-amplitude responses of the drive and sense vibrations are obtained,and comprehensive parametric studies are performed.The significant effects of system damping,excitation amplitude,drive amplitude and spinning speed on the responses are discussed,which will facilitate to improve the nonlinear performance and sensitivity of the gyroscope.展开更多
The nonlinear transverse vibrations of ordered and disordered twodimensional(2D)two-span composite laminated plates are studied.Based on the von Karman’s large deformation theory,the equations of motion of each-span ...The nonlinear transverse vibrations of ordered and disordered twodimensional(2D)two-span composite laminated plates are studied.Based on the von Karman’s large deformation theory,the equations of motion of each-span composite laminated plate are formulated using Hamilton’s principle,and the partial differential equations are discretized into nonlinear ordinary ones through the Galerkin’s method.The primary resonance and 1/3 sub-harmonic resonance are investigated by using the method of multiple scales.The amplitude-frequency relations of the steady-state responses and their stability analyses in each kind of resonance are carried out.The effects of the disorder ratio and ply angle on the two different resonances are analyzed.From the numerical results,it can be concluded that disorder in the length of the two-span 2D composite laminated plate will cause the nonlinear vibration localization phenomenon,and with the increase of the disorder ratio,the vibration localization phenomenon will become more obvious.Moreover,the amplitude-frequency curves for both primary resonance and 1/3 sub-harmonic resonance obtained by the present analytical method are compared with those by the numerical integration,and satisfactory precision can be obtained for engineering applications and the results certify the correctness of the present approximately analytical solutions.展开更多
基金the Deputyship for Research&Innovation,Ministry of Education in Saudi Arabia for funding this research work through the project number RI-44-0143
文摘In this investigation,some different approaches are implemented for analyzing a generalized forced damped complex Duffing oscillator,including the hybrid homotopy perturbation method(H-HPM),which is sometimes called the Krylov-Bogoliubov-Mitropolsky(KBM)method and the multiple scales method(MSM).All mentioned methods are applied to obtain some accurate and stable approximations to the proposed problem without decoupling the original problem.All obtained approximations are discussed graphically using different numerical values to the relevant parameters.Moreover,all obtained approximate solutions are compared with the 4thorder Runge-Kutta(RK4)numerical approximation.The maximum residual distance error(MRDE)is also estimated,in order to verify the high accuracy of the obtained analytic approximations.
文摘We propose a novel computational framework that is capable of employing different time integration algorithms and different space discretized methods such as the Finite Element Method,particle methods,and other spatial methods on a single body sub-dividedintomultiple subdomains.This is in conjunctionwithimplementing thewell known Generalized Single Step Single Solve(GS4)family of algorithms which encompass the entire scope of Linear Multistep algorithms that have been developed over the past 50 years or so and are second order accurate into the Differential Algebraic Equation framework.In the current state of technology,the coupling of altogether different time integration algorithms has been limited to the same family of algorithms such as theNewmarkmethods and the coupling of different algorithms usually has resulted in reduced accuracy in one or more variables including the Lagrange multiplier.However,the robustness and versatility of the GS4 with its ability to accurately account for the numerical shifts in various time schemes it encompasses,overcomes such barriers and allows a wide variety of arbitrary implicit-implicit,implicit-explicit,and explicit-explicit pairing of the various time schemes while maintaining the second order accuracy in time for not only all primary variables such as displacement,velocity and acceleration but also the Lagrange multipliers used for coupling the subdomains.By selecting an appropriate spatialmethod and time scheme on the area with localized phenomena contrary to utilizing a single process on the entire body,the proposed work has the potential to better capture the physics of a given simulation.The method is validated by solving 2D problems for the linear second order systems with various combination of spatial methods and time schemes with great flexibility.The accuracy and efficacy of the present work have not yet been seen in the current field,and it has shown significant promise in its capabilities and effectiveness for general linear dynamics through numerical examples.
基金the support of Texas A&M University at Qatar for the 2022 Sixth Cycle Seed Grant Project。
文摘The main objective of this paper is to investigate the influence of inertia of nonlinear springs on the dispersion behavior of discrete monoatomic chains with lumped and distributed masses.The developed model can represent the wave propagation problem in a non-homogeneous material consisting of heavy inclusions embedded in a matrix.The inclusions are idealized by lumped masses,and the matrix between adjacent inclusions is modeled by a nonlinear spring with distributed masses.Additionally,the model is capable of depicting the wave propagation in bi-material bars,wherein the first material is represented by a rigid particle and the second one is represented by a nonlinear spring with distributed masses.The discrete model of the nonlinear monoatomic chain with lumped and distributed masses is first considered,and a closed-form expression of the dispersion relation is obtained by the second-order Lindstedt-Poincare method(LPM).Next,a continuum model for the nonlinear monoatomic chain is derived directly from its discrete lattice model by a suitable continualization technique.The subsequent use of the second-order method of multiple scales(MMS)facilitates the derivation of the corresponding nonlinear dispersion relation in a closed form.The novelties of the present study consist of(i)considering the inertia of nonlinear springs on the dispersion behavior of the discrete mass-spring chains;(ii)developing the second-order LPM for the wave propagation in the discrete chains;and(iii)deriving a continuum model for the nonlinear monoatomic chains with lumped and distributed masses.Finally,a parametric study is conducted to examine the effects of the design parameters and the distributed spring mass on the nonlinear dispersion relations and phase velocities obtained from both the discrete and continuum models.These parameters include the ratio of the spring mass to the lumped mass,the nonlinear stiffness coefficient of the spring,and the wave amplitude.
基金support provided by the University Grant Commission,New Delhi,Government of India,under research grant no.37-515/2009(SR).
文摘A perturbation algorithm Multiple Scales Modified Lindstedt–Poincare(MSMLP),combination of method of Multiple Scales and modified Lindstedt–Poincare is proposed for the solution of Quintic Duffing equation which combines the advantages of both the methods.Solution obtained by the MSMLP method is compared with the Multiple Scales method and accurate closed form approximate solution of the Quintic Duffing equation.The proposed method produces better results for a wide range of amplitude values of oscillations and strong nonlinearities.Numerical simulation has been performed in MATHEMATICA 7.0.
基金Project supported by the National Natural Science Foundation of China (No. 11972204)。
文摘In this article, the nonlinear dynamic responses of sandwich functionally graded(FG) porous cylindrical shell embedded in elastic media are investigated. The shell studied here consists of three layers, of which the outer and inner skins are made of solid metal, while the core is FG porous metal foam. Partial differential equations are derived by utilizing the improved Donnell's nonlinear shell theory and Hamilton's principle. Afterwards, the Galerkin method is used to transform the governing equations into nonlinear ordinary differential equations, and an approximate analytical solution is obtained by using the multiple scales method. The effects of various system parameters,specifically, the radial load, core thickness, foam type, foam coefficient, structure damping,and Winkler-Pasternak foundation parameters on nonlinear internal resonance of the sandwich FG porous thin shells are evaluated.
文摘In this paper,we used time delay feedback to minimize the vibrations of a hybrid Rayleigh–van der Pol–Duffing oscillator.This system is a one-degree-offreedom containing the cubic and fifth nonlinear terms and an external force.We applied the multiple scales method to get the solution from first approximation.Graphically and numerically,we studied the system before and after adding time delay feedback at the primary resonance case(ffi!).We used MATLAB program to simulate the efficacy of different parameters and the time delay on the main system.
基金Supported by the National Natural Science Foundation of China(11201118)
文摘In this paper a modified L-P method and multiple scale method are used to solve sub-harmonic resonance solutions of strong and nonlinear resonance of general Van der Pol equation with parametric and external excitations by parametric transformation. Bifurcation response equation and transition sets of sub-harmonic resonance with strong nonlinearity of general Van der Pol equation with parametric and external excitation are worked out.Besides, transition sets and bifurcation graphs are drawn to help to analysis the problems theoretically. Conclusions show that the transition sets of general and nonlinear Van der Pol equation with parametric and external excitations are more complex than those of general and nonlinear Van der Pol equation only with parametric excitation, which is helpful for the qualitative and quantitative reference for engineering and science applications.
基金the National Natural Science Foundation of China(Nos.11702111,11572205,and 11732014)the Natural Science Foundation of Shandong Province of China(Nos.ZR2017QA005 and ZR2017BA031)。
文摘An electromagnetic nonlinear energy harvester(NEH)based on a rotating system is proposed,of which the host system rotates at a constant speed and vibrates harmonically in the vertical direction.This kind of device exhibits several resonant phenomena due to the combinations of the rotating and the vibration frequencies of the host system as well as the cubic nonlinearity of the NEH.The governing equation of motion for the NEH is derived,and the dynamic responses and output power are investigated with the multiple scale method under the 1:1 primary and 2:1 superharmonic resonant conditions.The effects of system parameters including the nondimensional external frequency,the rotating speed,and the nonlinear stiffness on the responses of free vibration for the system are studied.The results of the primary resonance show that the responses exhibit not only the resonant characteristics but also the nonlinear dynamic characteristics such as the saddle-node(SN)bifurcation.The coexistence of multiple solutions and the varying trends of responses are verified with the direct numerical simulation.Moreover,the effects of system parameters on the average output power are investigated.The results of the analyses on the two resonant conditions indicate that the large power can be harvested in two resonant frequency bands.The effect of resonance on the output power is dominant for the 2:1 superharmonic resonance.Moreover,the results also show that introducing the nonlinearity can increase the value of the output power in large frequency bands and induce the occurence of new frequency bands to harvest the large power.The efficiency of the harvested power could be improved by the combined effects of the resonance as well as the nonlinearity of the NEH device.Suitable parameter conditions could help optimize the power harvesting in design.
文摘The nonlinear interactions of a microarch resonator with 3:1 internal resonance are studied.The microarch is subjected to a combination of direct current(DC)and alternating current(AC)electric voltages.Thin piezoelectric layers are thoroughly bonded on the top and bottom surfaces of the microarch.The piezoelectric actuation is not only used to modulate the stiffness and resonance frequency of the resonator but also to provide the suitable linear frequency ratio for the activation of the internal resonance.The size effect is incorporated by using the so-called modified strain gradient theory.The system is highly nonlinear due to the co-existence of the initial curvature,the mid-plane stretching resulting from clamped anchors,and the electrostatic excitation.The eigenvalue problem is solved to conduct a frequency analysis and identify the possible regions for activating the internal resonance.The effects of the piezoelectric actuation,the electric excitation,and the small-scale effect are investigated on the internal resonance.Exclusive nonlinear phenomena such as Hopf bifurcation and hysteresis are identified in the microarch response.It is shown that by applying appropriate piezoelectric actuation,one is able to activate microarch internal resonance regardless of the initial rise level of the microarch.It is also disclosed that among all the parameters,AC electric voltage has the greatest effect on the energy exchange between the interacting modes.The results can be used to design resonators and internal resonance based micro-electro-mechanical system(MEMS)energy harvesters.
文摘An externally excited Duffing oscillator under feedback control is discussed and analyzed under the worst resonance case.Multiple time scales method is applied for this system to find analytic solution with the existence and nonexistence of the time delay on control loop.An appropriate stability analysis is also performed and appropriate choices for the feedback gains and the time delay are found in order to reduce the amplitude peak.Different response curves are involved to show and compare controller effects.In addition,analytic solutions are compared with numerical approximation solutions using Rung-Kutta method of fourth order.
基金This study was supported by the National Natural Science Foundation of China(Grant Nos.1167218911772009 and 12072311)and High-Level Talents Program of Yangzhou University.
文摘Micro-electro-mechanical system(MEMS)gyroscopes are an important sort of inertial sensor for identifying parameters of spinning structures,such as the spinning speed and angular deviation,based on the Coriolis effect.In this paper,the nonlinear mechanism of MEMS vibratory ring gyroscopes is analyzed by applying a fully coupled nonlinear model,in which the gyroscopic coupling and geometrically and structurally nonlinear couplings are all taken into account.The coupled differential equations governing the drive and sense motions are established via the Lagrangian equations.Numerical simulation is conducted,and the key nonlinear components and energy transfer behaviors between the drive and sense modes are elucidated.It is revealed that the cubic rigidity nonlinearity is another significant factor leading to the coupling between the drive and sense modes other than the gyroscopic coupling.Perturbation analysis is also carried out by using the method of multiple scales.The nonlinear frequency-amplitude responses of the drive and sense vibrations are obtained,and comprehensive parametric studies are performed.The significant effects of system damping,excitation amplitude,drive amplitude and spinning speed on the responses are discussed,which will facilitate to improve the nonlinear performance and sensitivity of the gyroscope.
基金This research is supported by the National Natural Science Foundation of China(Nos.11572007 and 11172084).
文摘The nonlinear transverse vibrations of ordered and disordered twodimensional(2D)two-span composite laminated plates are studied.Based on the von Karman’s large deformation theory,the equations of motion of each-span composite laminated plate are formulated using Hamilton’s principle,and the partial differential equations are discretized into nonlinear ordinary ones through the Galerkin’s method.The primary resonance and 1/3 sub-harmonic resonance are investigated by using the method of multiple scales.The amplitude-frequency relations of the steady-state responses and their stability analyses in each kind of resonance are carried out.The effects of the disorder ratio and ply angle on the two different resonances are analyzed.From the numerical results,it can be concluded that disorder in the length of the two-span 2D composite laminated plate will cause the nonlinear vibration localization phenomenon,and with the increase of the disorder ratio,the vibration localization phenomenon will become more obvious.Moreover,the amplitude-frequency curves for both primary resonance and 1/3 sub-harmonic resonance obtained by the present analytical method are compared with those by the numerical integration,and satisfactory precision can be obtained for engineering applications and the results certify the correctness of the present approximately analytical solutions.