The construction of non-linear normal modes for systems with internal resonance

A numerical method, based on the invariant manifold approach, is presented for constructing non-linear normal modes for systems with internal resonances. In order to parameterize the non-linear normal modes of interest, multiple pairs of system state variables involved in the internal resonance are kept as ‘seeds’ for the construction of the multi-mode invariant manifold. All the remaining degrees of freedom are then constrained to these ‘seed’, or master, variables, resulting in a system of non-linear partial differential equations that govern the constraint relationships, and these are solved numerically. The computationally-intensive solution procedure uses a combination of finite difference schemes and Galerkin-based expansion approaches. It is illustrated using two examples, both of which focus on the construction of two-mode models. The first example is based on the analysis of a simple three-degree-of-freedom example system, and is used to demonstrate the approach. An invariant manifold that captures two non-linear normal modes is constructed, resulting in a reduced order model that accurately captures the system dynamics. The methodology is then applied to a larger order system, specifically, an 18-degree-of-freedom rotating beam model that features a three-to-one internal resonance between the first two flapping modes. The accuracy of the non-linear two-mode reduced order model is verified by comparing time-domain simulations of the two DOF model and the full system equations of motion.     

N-pulse homoclinic orbits in perturbations of resonant hamiltonian systems

In this paper we develop an analytical method to detect orbits doubly asymptotic to slow manifolds in perturbations of integrable, two-degree-of-freedom resonant Hamiltonian systems. Our energy-phase method applies to both Hamiltonian and dissipative perturbations and reveals families of multi-pulse solutions which are not amenable to Melnikov-type methods. As an example, we study a two-mode approximation of the nonlinear, nonplanar oscillations of a parametrically forced inextensional beam. In this problem we find unusually complicated mechanisms for chaotic motions and verify their existence numerically.     

Lump,soliton, and interaction solutions to a generalized two-mode higher-order nonlinear evolution equation in plasma physics

This article investigates a nonlinear fifth-order partial differential equation (PDE) in two-mode waves. The equation generalizes two-mode Sawada-Kotera (tmSK), two-mode Lax (tmLax), and two-mode Caudrey–Dodd–Gibbon (tmCDG) equations. In 2017, Wazwaz [1] presented three two-mode fifth-order evolutions equations as tmSK, tmLax, and tmCDG equations for the integrable two-mode KdV equation and established solitons up to three-soliton solutions. In light of the research above, we examine a generalized two-mode evolution equation using a logarithmic transformation concerning the equation’s dispersion. It utilizes the simplified technique of the Hirota method to obtain the multiple solitons as a single soliton, two solitons, and three solitons with their interactions. Also, we construct one-lump solutions and their interaction with a soliton and depict the dynamical structures of the obtained solutions for solitons, lump, and their interactions. We show the 3D graphics with their contour plots for the obtained solutions by taking suitable values of the parameters presented in the solutions. These equations simultaneously study the propagation of two-mode waves in the identical direction with different phase velocities, dispersion parameters, and nonlinearity. These equations have applications in several real-life examples, such as gravity-affected waves or gravity-capillary waves, waves in shallow water, propagating waves in fast-mode and the slow-mode with their phase velocity in a strong and weak magnetic field, known as magneto-sound propagation in plasmas.

     

Two-mode squeezed state of nanomechanical resonators

In this paper, we introduce a flexible model for the control and measurement of NAMRs （nanomechanical resonators）. We obtain the free Hamiltonian of the dcSQUID （direct current superconducting quantum interference device） and the interaction Hamiltonian between these two NAMRs and the dc-SQUID by introducing the annihilation and creation operators under the rotating wave approximation. We can treat the mode of the dc-SQUID as a classical field. In the Heisenberg picture, the generation of two-mode squeezed states of two nanomechanical resonators is shown by their collective coordinate and momentum operators.     

The effects of fractional order on a 3-D quadratic autonomous system with four-wing attractor

In this paper, a fractional 3-dimensional (3-D) 4-wing quadratic autonomous system (Qi system) is analyzed. Time domain approximation method (Grunwald–Letnikov method) and frequency domain approximation method are used together to analyze the behavior of this fractional order chaotic system. It is found that the decreasing of the system order has great effect on the dynamics of this nonlinear system. The fractional Qi system can exhibit chaos when the total order less than 3, although the regular one always shows periodic orbits in the same range of parameters. In some fractional order, the 4 wings are decayed to a scroll using the frequency domain approximation method which is different from the result using time domain approximation method. A surprising finding is that the phase diagrams display a character of local self-similar in the 4-wing attractors of this fractional Qi system using the frequency approximation method even though the number and the characteristics of equilibria are not changed. The frequency spectrums show that there is some shrinking tendency of the bandwidth with the falling of the system states order. However, the change of fractional order has little effect on the bandwidth of frequency spectrum using the time domain approximation method. According to the bifurcation analysis, the fractional order Qi system attractors start from sink, then period bifurcation to some simple periodic orbits, and chaotic attractors, finally escape from chaotic attractor to periodic orbits with the increasing of fractional order α in the interval [0.8,1]. The simulation results revealed that the time domain approximation method is more accurate and reliable than the frequency domain approximation method.     

A monotone finite volume scheme for advection–diffusion equations on distorted meshes

A new monotone finite volume method with second‐order accuracy is presented for the steady‐state advection–diffusion equation. The method uses a nonlinear approximation for both diffusive and advective fluxes that guarantee the positivity of the numerical solution. The approximation of the diffusive flux is based on nonlinear two‐point approximation, and the approximation of the advective flux is based on the second‐order upwind method with proper slope limiter. The second‐order convergence rate for concentration and the monotonicity of the nonlinear finite volume method are verified with numerical experiments. Copyright © 2011 John Wiley & Sons, Ltd.     

Three-to-One Internal Resonances in Hinged-Clamped Beams

The nonlinear planar response of a hinged-clamped beam to a primary excitation of either its first mode or its second mode is investigated. The analysis accounts for mid-plane stretching, a static axial load and a restraining spring at one end, and modal damping. For a range of axial loads, the second natural frequency is approximately three times the first natural frequency and hence the first and second modes may interact due to a three-to-one internal resonance. The method of multiple scales is used to attack directly the governing nonlinear partial-differential equation and derive two sets of four first-order nonlinear ordinary-differential equations describing the modulation of the amplitudes and phases of the first two modes in the case of primary resonance of either the first or the second mode. Periodic motions and periodically and chaotically modulated motions of the beam are determined by investigating the equilibrium and dynamic solutions of the modulation equations. For the case of primary resonance of the first mode, only two-mode solutions are possible, whereas for the case of primary resonance of the second mode, single- and two-mode solutions are possible. The two-mode equilibrium solutions of the modulation equations may undergo a supercritical or a subcritical Hopf bifurcation, depending on the magnitude of the axial load. A shooting technique is used to calculate limit cycles of the modulation equations and Floquet theory is used to ascertain their stability. The limit cycles correspond to periodically modulated motions of the beam. The limit cycles are found to undergo cyclic-fold bifurcations and period-doubling bifurcations, leading to chaos. The chaotic attractors may undergo boundary crises, resulting in the destruction of the chaotic attractors and their basins of attraction.     

The third‐order polynomial method for two‐dimensional convection and diffusion

Using the upstream polynomial approximation a series of accurate two‐dimensional explicit numerical schemes is developed for the solution of the convection–diffusion equation. A third‐order polynomial approximation (TOP) of the convection term and a consistent second‐order approximation of the diffusion term are combined in a single‐step flux‐difference algorithm. Stability analysis confirms that the TOP‐12 scheme satisfies the CFL condition for two dimensions. Using smaller and narrower flux stencils compared to algorithms of similar accuracy, the TOP‐12 scheme is more efficient in terms of computations per single node. Numerical tests and comparison with other well‐known algorithms show a high performance of the developed schemes. Copyright © 2003 John Wiley & Sons, Ltd.     

A multi-sensing scheme based on nonlinear coupled micromachined resonators

A new multi-sensing scheme via nonlinear weakly coupled resonators is introduced in this paper, which can simultaneously detect two different physical stimuli by monitoring the dynamic response around the first two lowest modes. The system consists of a mechanically coupled bridge resonator and cantilever resonator. The eigenvalue problem is solved to identify the right geometry for the resonators to optimize their resonance frequencies based on mode localization in order to provide outstanding sensitivity. A nonlinear equivalent model is developed using the Euler–Bernoulli beam theory while accounting for the geometric and electrostatic nonlinearities. The sensor's dynamics are explored using a reduced-order model based on two-mode Galerkin discretization, which reveals the richness of the response. To demonstrate the proposed sensing scheme, the dynamic response of the weakly coupled resonator is investigated by tuning the stiffness and mass of the bridge and cantilever resonators, respectively. With its simple and scalable design, the proposed system shows great potential for intelligent multi-sensing detection in many applications.

     

Reduction of the Mathieu equation to a nonlinear equation of the first order

A second order equation with periodic coefficients is considered. It is shown that its analysis can be reduced to the study of a nonlinear equation of the first order. The second approximation is obtained for the first resonance region of the Mathieu equation. This approximation describes the behavior of solutions inside this resonance region and near it.     

Numerical solution of the incompressible Navier–Stokes equations with an upwind compact difference scheme

A new finite difference method for the discretization of the incompressible Navier–Stokes equations is presented. The scheme is constructed on a staggered‐mesh grid system. The convection terms are discretized with a fifth‐order‐accurate upwind compact difference approximation, the viscous terms are discretized with a sixth‐order symmetrical compact difference approximation, the continuity equation and the pressure gradient in the momentum equations are discretized with a fourth‐order difference approximation on a cell‐centered mesh. Time advancement uses a three‐stage Runge–Kutta method. The Poisson equation for computing the pressure is solved with preconditioning. Accuracy analysis shows that the new method has high resolving efficiency. Validation of the method by computation of Taylor's vortex array is presented. Copyright © 1999 John Wiley & Sons, Ltd.     

NONLINEAR INTERNAL RESONANCE OF FUNCTIONALLY GRADED CYLINDRICAL SHELLS USING THE HAMILTONIAN DYNAMICS

Internal resonance in nonlinear vibration of functionally graded （FG） circular cylin- drical shells in thermal environment is studied using the Hamiltonian dynamics formulation. The material properties are considered to be temperature-dependent. Based on the Karman-Donnell＇s nonlinear shell theory, the kinetic and potential energy of FG cylindrical thin shells are formu- lated. The primary target is to investigate the two-mode internal resonance, which is triggered by geometric and material parameters of shells. Following a secular perturbation procedure, the underlying dynamic characteristics of the two-mode interactions in both exact and near resonance cases are fully discussed. It is revealed that the system will undergo a bifurcation in near resonance case, which induces the dynamic response at high energy level being distinct from the motion at low energy level. The effects of temperature and volume fractions of composition on the exact resonance condition and bifurcation characteristics of FG cylindrical shells are also investigated.     

压电桁架作动器/传感器优化配置算法研究

针对自适应压电桁架结构振动控制,建立了作动器/传感器优化配置数学模型,并提出一种优化配置的新方法。为了减少结构分析次数,该方法将近似概念、对偶法和遗传算法相结合,首先采用多点近似技术建立原问题的序列近似问题,再对近似问题中的作动器/传感器位置离散变量和控制增益连续变量采用遗传算法和对偶方法分别寻优的分层优化策略。为了提高近似问题对原问题的逼近程度,本文提出一种适于离散变量结构优化的分段多点近似函数。算例表明本文方法能够以很少的结构分析次数得到最优解。     

FInite torsion of slightly compressible rubberlike circular cylinders∗

A regular perturbation procedure and the Rayleigh-Ritz method are used to study the finite torsion of cylinders made of slightly compressible neo-Hookean materials. Two cases are considered: 1. the length of the cylinder is not permitted to change during torsion and 2. the net axial force on the cylinder vanishes. The perturbation procedure fails for certain constitutive relations whereas, in principle, the Rayleigh-Ritz method has general applicability. When it works, the success of the perturbation procedure depends on prior knowledge of the problem for an incompressible material (the zeroth order nonlinear problem). The solution of problem 2. is considerably more complicated than the solution of problem 1. since the complete approximation of order n for problem 2. requires extensive work on the approximation of order (n + 1).     

A Universal Matrix Perturbation Technique for Complex Modes

ＩｎｔｒｏｄｕｃｔｉｏｎＭａｔｒｉｘｐｅｒｔｕｒｂａｔｉｏｎｍｅｔｈｏｄｓｆｏｒｔｈｅｄｙｎａｍｉｃｒｅａｎａｌｙｓｉｓｏｆｓｅｌｆ＿ａｄｊｏｉｎｔｓｙｓｔｅｍｓｈａｖｅｂｅｅｎｗｅｌｌｄｅｖｅｌｏｐｅｄ[1] .Ｈｏｗｅｖｅｒ,ｍａｎｙｓｙｓｔｅｍｓｇｉｖｅｒｉｓｅｔｏｇｅｎｅｒａｌｎｏｎ＿ｓｅｌｆ＿ａｄｊｏｉｎｔｆｏｒｍｕｌａｔｉｏｎｓ.Ｉｍｐｏｒｔａｎｔｅｘａｍｐｌｅｓａｒｅａｅｒｏｅｌａｓｔｉｃｓｔａｂｉｌｉｔｙｏｆｓｙｓｔｅｍｓ,ａｒｂｉｔｒａｒｉｌｙｄａｍｐｅｄｏｒｇｙｒｏｓｃｏｐｉｃｓｙｓ…     

Weak Non-Linear Analysis of Moderate Stefan Number Stationary Convection in Rotating Mushy Layers

The effects of rotation on a mushy layer, during the solidification of binary alloys, is considered. A near-eutectic approximation and large far-field temperature are employed in order to decouple the mushy layer from the overlying liquid melt. The current study employs a new moderate time scale for mushy layers exhibiting Stefan numbers of unit order of magnitude. The weak non-linear theory is used to evaluate the leading order amplitude. The results of the weak non-linear theory are then used to establish the nature of the bifurcation, that is whether the bifurcation is forward or inverse.     

一类高精度TVD差分格式及其应用

构造了一维非线性双曲型守恒律的一个新的高精度、高分辨率的守恒型TvD差分格式。其构造思想是：首先，将计算区间划分为若干个互不相交的小区间，再根据精度要求等分小区间，通过各细小区间上的单元平均状态变量，重构各细小区间交界面上的状态变量，并加以校正；其次，利用近似Riemann解计算细小区间交界面上的数值通量，并结合高阶Runge—Kutta TVD方法进行时间离散，得到了高精度的全离散方法。证明了该格式的TVD特性。该格式适合于使用分量形式计算而无须进行局部特征分解。通过计算几个典型的问题，验证了格式具有高精度、高分辨率且计算简单的优点。     

An approach for obtaining approximate formulas for the Rayleigh wave velocity

In this paper, we introduce an approach for finding analytical approximate formulas for the Rayleigh wave velocity for isotropic elastic solids and anisotropic elastic media as well. The approach is based on the least-square principle. To demonstrate its application, we applied it in order to obtain an explanation for Bergmann’s approximation, the earliest known approximation of the Rayleigh wave velocity for isotropic elastic solids, and used it to establish a new approximation. By employing this approach, the best approximate polynomials of the second order of the cubic power and the quartic power in the interval [0, 1] were found. By using the best approximate polynomial of the second order of the cubic power, we derived an approximate formula for the Rayleigh wave speed in isotropic elastic solids which is slightly better than the one given recently by Rahman and Michelitsch by employing Lanczos’s approximation. Also by using this second order polynomial, analytical approximate expressions for orthotropic, incompressible and compressible elastic solids were found. For incompressible case, it is shown that the approximation is comparable with Rahman and Michelitsch’s approximation, while for the compressible case, it is shown that our approximate formulas are more accurate than Mozhaev’s ones. Remarkably, by using the best approximate polynomials of the second order of the cubic power and the quartic power in the interval [0, 1], we derived an approximate formula of the Rayleigh wave velocity in incompressible monoclinic materials, where the explicit exact formulas of the Rayleigh wave velocity so far are not available.     

On the “Best” Mode Form in the Mode Approximation Technique Using the Finite Element Method

A systematic procedure for choosing the “best” mode shape is discussed by extending the mode approximation technique. The finite element method is used to calculate several valid mode forms of a given structure. Nonlinear eigenvalue problems are solved by an iterative procedure in order to obtain mode forms. Since one can choose only one mode form in the mode approximation technique due to nonlinearity of viscoplastic materials, the lower bound theorem on final time is applied to identify the best mode form. The validity of the concept is demonstrated by numerical results from two example problems of clamped beam and portal frame.     

Aspects of the transpiration model for aerofoil design

Transpiration is a technique in which extra non-physical normal flows are created on an aerofoil surface in order to form a new streamline pattern such that the surface streamlines no longer follow the aerofoil surface under inviscid flow. The transpiration model is an important technique adopted in aerofoil design either to avoid mesh regeneration when aerofoil profile co-ordinates are adjusted or to find shape corrections in inverse design methods. A first-order approximation (with respect to the normal streamline displacement) to the transpiration model is commonly adopted; it is shown that this can be a poor approximation especially in regions of high curvature. In this paper more accurate approximations are developed to address this problem and improve the accuracy.