Calculation of solutions to the Mathieu equation and of related quantities

For the Mathieu equation, we consider finding eigenvalues with a given index (on the basis of oscillation theorems for the relevant difference equations), the stability of solutions to the difference equations, correct definition and calculation of eigenvalues and Mathieu functions with noninteger numbers, correct definition and calculation of the Mathieu characteristic exponent, and the calculation of values of solutions to the Mathieu equation for large arguments. Numerical algorithms are proposed for the problems listed above.     

Cauchy problem for Mathieu’s equation at parametric resonance

Mathieu’s equation is solved by an asymptotic averaging method in the fourth approximation for the first to fourth resonance domains and in the third approximation for the zero resonance domain. The general periodic and aperiodic solutions on characteristic curves are found, and the general solution is obtained in instability domains and stability-domain areas adjacent to the characteristic curves. All the solutions are explicitly found in the form of functions of an argument without using the auxiliary parameter employed in Whittaker’s method. Simple formulas depending on two parameters of the equation are derived for the characteristic exponent in instability domains and for the frequency of slow oscillations in stability domains near the characteristic curves. The theory is developed by analyzing the resonances exhibited by Mathieu’s equation.     

Effect of an elastic core on the dynamic stability of an orthotropic cylindrical shell

A method of determining the regions of dynamic instability of an orthotropic cylindrical shell "bonded" to an elastic cylinder is proposed. An expression for the core reaction is obtained from the coupling conditions for the forces normal to the lateral surface and the radial displacements of the shell and the core at the contact surface. When the reaction is substituted in the system of equations of motion of the shell, the part corresponding to the free vibrations of the cylinder is discarded. The system of equations of motion of the shell is reduced to an equation of Mathieu type, from which transcendental equations for determining the boundaries of the regions of dynamic instability are obtained. These regions are analyzed for various modes of loss of stability and different values of the core modulus of elasticity.     

Nonlinear Mathieu equation and coupled resonance mechanism

The method of multiple-scales is used to determine a third-order solution for a cubic nonlinear Mathieu equation. The perturbation solutions are imposed on the so-called solvability conditions. Solvability conditions in the non-resonance case yield the standard Landau equation. Several types of a parametric Landau equation are derived in the neighborhood of five different resonance cases. These parametric Landau equations contain a parametric complex conjugate term or a parametric second-order complex conjugate term or a parametric complex conjugate term as well as a parametric second-order term. Necessary and sufficient conditions for stability are performed in each resonance case. Stability criteria correspond to each parametric Landau equation and are derived by linear perturbation. Stability criteria for the non-trivial steady-state response are discussed. The analysis leads to simultaneous resonance. Transition curves are performed in each case. Numerical calculations are made for some transition curves to illustrate the coupled resonance regions, where the induced stability tongues within the instability tongues are observed. The amplitude of the periodic coefficient of Mathieu equation plays a dual role in the stability criteria for nonlinear Mathieu equation.     

Output feedback vibration control of a string driven by a nonlinear actuator

In this paper, we design an observer-based output feedback controller to exponentially stabilize a system of nonlinear ordinary differential equation-wave partial differential equation-ordinary differential equation. An observer is designed to estimate the full states of the system using available boundary values of the partial differential equation. The output feedback controller is built via the combination of the ordinary differential equation backstepping which is applied to deal with the nonlinear ordinary differential equation, and the partial differential equation backstepping which is used for the wave partial differential equation-ordinary differential equation. The controller can be applied into vibration suppression of a string-payload system driven by an actuator with nonlinear characteristics. The global exponential stability of all states in the closed-loop system is proved by Lyapunov analysis. The numerical simulation illustrates the states of the actuator, string, payload and the observer errors are fast convergent to zero under the proposed output feedback controller.     

Dynamics of a congestion control model in a wireless access network

In this paper, a congestion control algorithm with heterogeneous delays in a wireless access network is considered. We regard the communication time delay as a bifurcating parameter to study the dynamical behaviors, i.e., local asymptotical stability, Hopf bifurcation and resonant codimension-two bifurcation. By analyzing the associated characteristic equation, the Hopf bifurcation occurs when the delay passes through a sequence of critical value. Furthermore, the direction and stability of the bifurcating periodic solutions are derived by applying the normal form theory and the center manifold theorem. In the meantime, the resonant codimension-two bifurcation is also found in this model. Some numerical examples are finally performed to verify the theoretical results.     

Rotating orbits of a parametrically-excited pendulum

The authors consider the dynamics of the harmonically excited parametric pendulum when it exhibits rotational orbits. Assuming no damping and small angle oscillations, this system can be simplified to the Mathieu equation in which stability is important in investigating the rotational behaviour. Analytical and numerical analysis techniques are employed to explore the dynamic responses to different parameters and initial conditions. Particularly, the parameter space, bifurcation diagram, basin of attraction and time history are used to explore the stability of the rotational orbits. A series of resonance tongues are distributed along the non-dimensionalied frequency axis in the parameter space plots. Different kinds of rotations, together with oscillations and chaos, are found to be located in regions within the resonance tongues.     

Stochastic semidiscretization method: Second moment stability analysis of linear stochastic periodic dynamical systems with delays

In this paper, an efficient numerical approach is presented, which allows the analysis of the moment dynamics, stability, and stationary behavior of linear periodic stochastic delay differential equations. The method leads to a high dimensional stochastic mapping with periodic statistical properties, from which the periodic first and second moment mappings are derived. The application of the method is demonstrated first through the analysis of the stochastic delay Mathieu equation. Then a practical case study, where the effect of spindle speed variation on the stability and the resulting surface quality of turning operations is investigated.     

具有脉冲效应的两食饵一捕食者系统分析

构建并分析了一个在固定时刻脉冲投放捕食者且具有功能性反应的两食饵一捕食者系统,应用脉冲比较定理和微分方程的分析方法,得到了食饵灭绝周期解稳定的条件和系统持续生存的条件,并数值分析了所得的理论结果.     

Finite‐region stabilization via dynamic output feedback for 2‐D Roesser models

Finite‐region stability (FRS), a generalization of finite‐time stability, has been used to analyze the transient behavior of discrete two‐dimensional (2‐D) systems. In this paper, we consider the problem of FRS for discrete 2‐D Roesser models via dynamic output feedback. First, a sufficient condition is given to design the dynamic output feedback controller with a state feedback‐observer structure, which ensures the closed‐loop system FRS. Then, this condition is reducible to a condition that is solvable by linear matrix inequalities. Finally, viable experimental results are demonstrated by an illustrative example.     

Adoption of the Numerical Method of Chebyshev Polynomials to the Stability Analysis of Delayed DEs

We propose a new way of the application of the numerical method based on Chebyshev approximation concerning the stability analysis of (periodic) delay‐differential equations. The technique is tested on the delayed Mathieu equation.     

A newton approach for long term stability studies in power systems

This paper deals with the problem of transient and long term stability of power systems. The issue of assessing both horizons of analysis is particularly focused. This is because the long term stability may be studied by a simplified algebraic model which also captures some dynamic characteristics. Such an approach is called quasi-dynamic model. The idea of analyzing the transient period and migrating to the quasi-dynamic model is addressed in this paper. The theoretical foundation is presented and some tests are carried out in order to validate the approach.     

Pure parametric excitation of a micro cantilever beam actuated by piezoelectric layers

This paper deals with the stability analysis of transverse motions of a cantilever microbeam sandwiched by two piezoelectric layers located on the lower and upper surfaces of the microbeam. Application of same DC and AC voltages to the upper and lower piezoelectric layers creates an axial force with steady and time-varying components. The eigenfunction expansion of the transverse motion equation leads to the creation of a Mathieu type parametric equation which is mostly seen in the stability analysis of the structures in the literature; using Floquet theory for single degree of freedom systems the stable and unstable regions of the problem are investigated. The effect of viscous damping and DC voltage on the stability region of the problem is also studied. The results show the stabilizing effect of the viscous damping and positive DC voltage on the behavior of the microbeam. The achieved results are finally compared with those reported in the literature.     

弹塑性杆在刚性块轴向撞击下的动力屈曲

基于能量原理,对弹塑性杆在刚性块轴向撞击下的动力屈曲问题进行了讨论．用特征线法分析了刚性块轴向撞击弹塑性直杆时应力波传播的过程．考虑了弹塑性应力波传播对屈曲的影响,建立了该问题横向扰动方程．用幂级数解法,理论上给出了该问题的级数解．分析解的性质,得到了发生屈曲时的临界条件．通过理论分析和数值计算,得到了临界速度与冲击质量、临界长度及线性强化模量间的关系．     

Asymptotic solutions and stability analysis for generalized non-homogeneous Mathieu equation

The asymptotic solutions and transition curves for the generalized form of the non-homogeneous Mathieu differential equation are investigated in this paper. This type of governing differential equation of motion arises from the dynamic behavior of a pendulum undergoing a butterfly-type end support motion. The strained parameter technique is used to obtain periodic asymptotic solutions. The transition curves for some special cases are presented and their corresponding periodic solutions with the periods of 2π and 4π are evaluated. The stability analyses of those transition curves in the ε–δ plane are carried out, analytically, using the multiple scales method. The numerical simulations for some typical points in the ε–δ plane are performed and the dynamic characteristics of the resulting phase plane trajectories are discussed.     

动载下缝端应力强度因子计算的扩展有限元法

在扩展有限元法(extended finite element methods, XFEM)的理论框架下,重点研究了动荷载作用下稳定裂纹尖端动态应力强度因子(dynamic stress intensity factors, DSIFs)的求解方法．根据XFEM的位移模式,推导了动力XFEM的支配方程,采用Newmark隐式算法进行时间积分同时,提出一种XFEM质量矩阵的集中策略,给出了求解DSIFs的相互作用积分方法,与静态问题的相互作用积分方法相比,增加了惯性项的贡献．最后,若干典型算例的计算结果表明:XFEM可以准确评价稳定裂纹尖端的DSIFs,建议的质量矩阵集中策略是有效的,为得到正确的DSIFs,惯性项的贡献不可忽略.     

An Efficient and Accurate Spectral Method for Acoustic Scattering in Elliptic Domains

An efficient and accurate method for solving the two-dimensional Helmhokz equation in domains exterior to elongated obstacles is developed in this paper. The method is based on the so called transformed field expansion (TFE) coupled with a spectral-Galerkin solver for elliptical domain using Mathieu functions. Numerical results are presented to show the accuracy and stability of the proposed method.     

非线性发展方程的H(o)lder连续惯性流形

引入非线性发展方程的H\"older连续惯性流形的概念,为原来惯性流形概念的推广和修正.惯性流形是有限维不变的Lipschiz流形,是研究发展方程解的长时间性态的合适工具,其缺点是需要谱间隙条件.提出H\"older连续惯性流形也是有限维不变的,但光滑性减弱为H\"older连续,不需要谱间隙条件.该流形与指数吸引子交集具有指数吸引性,无穷维动力系统可在H\"older连续惯性流形上约化为有限维常微分方程组.     

非线性阻尼作用下标准线性固体粘弹性Ⅲ型破裂的解析解

把非线性Rayleigh阻尼引入标准线性固体粘弹性介质的Ⅲ型破裂的控制方程中,此方程是一个偏微分积分方程;首先设法消去积分项,得到一个三阶非线性偏微分方程,然后用小参数摄动法,得出线性化的各阶渐近控制方程;把每一个具有变系数的三阶线性控制方程分解为弹性部分及剩余部份,而前者的解析解为已知,后者是一个二阶变系数线性偏微分方程;它化不成Mathieu方程,也化不成Hill方程,故采用WKBJ的方法得出其渐近的解析解。