Responses of a non-linearly coupled pitch-roll ship model under harmonic excitation

Stationary responses of nonlinearly coupled pitch and roll ship modes are studied using a modified averaging method, along with two second order multiple time scale (MTS) methods for comparison. Stability of the solutions is also studied. In the case of harmonic excitation all the three methods give fairly accurate results to the original problem but the averaging method is the most efficient. Analytic solutions are obtained from the averaged equations, which can be used to predict stationary responses both for small and for large excitations. From the averaging method several qualitatively different phenomena which cannot be addressed by the first order theory have been obtained: (i) the saturation phenomenon is lost, (ii) the bifurcation points are altered and (iii) a drift term is present which, although small, appears to have a significant effect on the accuracy of the solutions.     

The Method of Averaging for Euler's Equations of Rigid Body Motion

We formulate the method of averaging for perturbations of Euler's equations of rotational motion. Euler's equations are three strongly nonlinear coupled differential equations that can be viewed as a three dimensional oscillator. The method of averaging is used to determine the long-term influence of perturbation terms on the motion by averaging about the nominal rigid body motion. The treatment is applicable to a large class of motions including precession with large nutation – it is not restricted to small motions about simple spins or nearly axi-symmetric bodies. Three examples are shown that demonstrate the accuracy of the method's predictions.     

Parametric excitation in non-linear dynamics

Consider a one-mass system with two degrees of freedom, non-linearly coupled, with parametric excitation in one direction. Assuming the internal resonance 1:2 and parametric resonance 1:2 we derive conditions for stability of the trivial solution by using both the harmonic balance method and the normal form method of averaging. If the trivial solution becomes unstable, a stable periodic solution may emerge, there are also cases where the trivial solution is stable and co-exists with a stable periodic solution; if both the trivial solution and the periodic solution(s) are unstable, we find an attracting torus with large amplitudes by a Neimark-Sacker bifurcation. The results of the harmonic balance method and averaging are compared, as well as the results on the Neimark-Sacker bifurcation obtained by the numerical software package CONTENT and by averaging. In all cases we have good agreement.     

Periodic solution of turbulent oscillating channel flows

The time-dependent turbulent Navier–Stokes equations are solved numerically by a finite element method with an algebraic eddy viscosity model (Baldwin–Lomax formulation) for oscillating turbulent channel flows. The method of averaging is used to analyse the resulting periodic motion of the fluid. Numerical results are obtained for various Strouhal numbers and relative amplitudes. A comparison is made between the numerical and published experimental results. It appears that for low relative amplitudes in a certain range of frequencies the agreement is satisfactory.     

Complex Inner Product Averaging Method for Calculating Normal Form of Ode

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

Resolving Controversies in the Application of the Method of Multiple Scales and the Generalized Method of Averaging

I compare application of the method of multiple scales with reconstitution and the generalized method of averaging for determining higher-order approximations of three single-degree-of-freedom systems and a two-degree-of-freedom system. Three implementations of the method of multiple scales are considered, namely, application of the method to the system equations expressed as second-order equations, as first-order equations, and in complex-variable form. I show that all of these methods produce the same modulation equations.I address the problem of determining higher-order approximate solutions of the Duffing equation in the case of primary resonance. I show that the conclusions of Rahman and Burton that the method of multiple scales, the generalized method of averaging, and Lie series and transforms might lead to incorrect results, in that spurious solutions occur and the obtained frequency–response curves bear little resemblance to the actual response, is the result of their using parameter values for which the neglected terms are the same order as the retained terms. I show also that spurious solutions cannot be avoided, in general, in any consistent expansion and their presence does not constitute a limitation of the methods. In particular, I show that, for the Duffing equation, the second-order frequency–response equation does not possess spurious solutions for the case of hardening nonlinearity, but possesses spurious solutions for the case of softening nonlinearity. For sufficiently small nonlinearity, the spurious solutions are far removed from the actual response. But as the strength of the nonlinearity increases, these solutions move closer to the backbone and eventually distort it. This is not a drawback of the perturbation methods but an indication of an application of the analysis for parameter values outside the range of validity of the expansion.Also, I address the problem of obtaining non-Hamiltonian modulation equations in the application of the method of multiple scales to multi-degree-of-freedom Hamiltonian systems written as second-order equations in time and how this problem can be overcome by attacking the state-space form of the governing equations. Moreover, I show that application of a variation of the method of Rahman and Burton to multi-degree-of-freedom systems leads to results that do not agree with those obtained with the generalized method of averaging.Contributed by Prof. R.A. Ibrahim.     

范式理论与平均法的等价性分析

分析了范式理论与平均法的等价性。得到的结论是：对含有一对纯虚根的二维非线性系统，使用两种方法得到的结果是等价的，并提供了两个算例来证实其结论的正确性。虽然分析是针对一类二维非线性系统，但其结论同样适合于高维非线性系统。     

具有多个极限环非线性动力系统的解析近似

应用一种新的解析方法------同伦分析法，研究了一种具有多个 极限环的Rayleigh振子问题. 与所有其他传统方法不同，该方法不依赖于小参数， 且提供了一个简便的途径以确保级数解的收敛, 因此，特别适用于强非线性问题. 将同伦分析法与平均法以及四阶的龙格库塔方法(数值解)做了比较. 结果 表明，平均法在强非线性情况失效, 四阶的龙格库塔法不能找到非稳定的极限环，而同伦分析法不仅适用于强非线性情 况，而且给出了非稳定的极限环.     

Response of quasi-integrable Hamiltonian systems with delayed feedback bang–bang control

The response of quasi-integrable Hamiltonian systems with delayed feedback bang–bang control subject to Gaussian white noise excitation is studied by using the stochastic averaging method. First, a quasi-Hamiltonian system with delayed feedback bang–bang control subjected to Gaussian white noise excitation is formulated and transformed into the Itô stochastic differential equations for quasi-integrable Hamiltonian system with feedback bang–bang control without time delay. Then the averaged Itô stochastic differential equations for the later system are derived by using the stochastic averaging method for quasi-integrable Hamiltonian systems and the stationary solution of the averaged Fokker–Plank–Kolmogorov (FPK) equation associated with the averaged Itô equations is obtained for both nonresonant and resonant cases. Finally, two examples are worked out in detail to illustrate the application and effectiveness of the proposed method and the effect of time delayed feedback bang–bang control on the response of the systems.     

Upscaling of Nonwetting Phase Residual Transport in Porous Media: A Network Approach

Based on the volume averaging method, a macroscopic model is developed for the upscaling of NAPL transport in a porous medium idealised by a network model. Under the assumption of local mass non-equilibrium, a macroscopic equation involving a dispersion tensor, additional convective terms and a linear form for the interfacial mass flux is obtained. The resolution of the two local closure problems obtained allow the determination of the local properties without adjustable parmeters. These problems are solved in a semi-analytical, semi-numerical manner on the network. The originality of this work is the association of the upscaling by volume averaging method with the network approach. The local properties, including the dispersion tensor and the mass exchange coefficient, can therefore be calculated over a large number of pore-bodies and pore-throats in a computationaly tractable manner, thus leading to more significant results. Results are presented for 3D, spatially periodic models of porous media.     

功能梯度简支矩形板的非线性动力响应

研究了功能梯度简支矩形板在横向简谐激励作用下的非线性动力响应问题。采用幂律分布规律描述功能梯度材料的等效材料参数,基于Galerkin法建立了系统广义坐标的常微分控制方程。利用平均法得到了系统的幅频响应特性,分析了功能梯度矩形薄板的非线性主共振特性。数值算例验证了平均化方法的正确性,揭示了功能梯度平板主共振响应中的多值性和跳跃现象;同时分析发现初始条件会改变功能梯度平板主共振的响应幅值。最后讨论了功能梯度材料的梯度指数对系统幅值响应的影响。     

First passage failure of quasi-partial integrable generalized Hamiltonian systems

The first passage failure of quasi-partial integrable generalized Hamiltonian systems is studied by using the stochastic averaging method. First, the stochastic averaging method for quasi-partial integrable generalized Hamiltonian systems is introduced briefly. Then, the backward Kolmogorov equation governing the conditional reliability function and the Pontryagin equation governing the conditional mean of first passage time are derived from the averaged Itô equations. The conditional reliability function, the conditional probability density and mean of the first passage time are obtained from solving these equations together with suitable initial condition and boundary conditions, respectively. Finally, one example is given to illustrate the proposed procedure in detail and the solutions are confirmed by using the results from Monte Carlo simulation of the original system.     

谐和与随机噪声联合作用下非线性系统的响应

研究了Duffing振子在谐和与随噪声联合激励下的响应和稳应性问题。用谐波平均法分析了系统在确定性谐和激励和随机激励联合作用下的响应,用随机平均法讨论了随机扰动项对系统晌应的影响。在一定条件下,系统具有两个均方响应值和跳跃现象。数值模拟表明本文提出的方法是有效的。     

Stability and chaotic motion in columns of nonlinear viscoelastic material

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｐｒｏｂｌｅｍｏｆｔｈｅｄｙｎａｍｉｃａｌｓｔａｂｉｌｉｔｙｏｆｖｉｓｃｏｅｌａｓｔｉｃｓｔｒｕｃｔｕｒｅｓｉｓｍｏｒｅｃｏｍｐｌｉｃａｔｅｄｔｈａｎｔｈａｔｏｆｅｌａｓｔｉｃｓｔｒｕｃｔｕｒｅｓ,ｓｉｎｃｅｔｈｅｓｔｒｅｓｓ＿ｓｔｒａｉｎｒｅｌａｔｉｏｎｏｆｖｉｓｃｏｅｌａｓｔｉｃｍａｔｅｒｉａｌｓｏｆｔｅｎｌｅａｄｓｔｏａｎｉｎｔｅｇｒｏ＿ｄｉｆｆｅｒｅｎｔｉａｌｅｑｕａｔｉｏｎｏｆｍｏｔｉｏｎ .Ｉｎｔｈｅｅａｒｌｙｒｅｓｅａｒｃｈｅｓｏｎｌｉｎｅａｒｖｉｓｃｏ…     

DISCRETE AND CONTINUUM MODELLING OF GRANULAR FLOW

This paper analyses three popular methods simulating granular flow at different time and length scales： discrete element method （DEM）, averaging method and viscous, elastic-plastic continuum model. The theoretical models of these methods and their applications to hopper flows are discussed. It is shown that DEM is an effective method to study the fundamentals of granular flow at a particle or microscopic scale. By use of the continuum approach, granular flow can also be described at a continuum or macroscopic scale. Macroscopic quantities such as velocity and stress can be obtained by use of such computational method as FEM. However, this approach depends on the constitutive relationship of materials and ignores the effect of microscopic structure of granular flow. The combined approach of DEM and averaging method can overcome this problem. The approach takes into account the discrete nature of granular materials and does not require any global assumption and thus allows a better understanding of the fundamental mechanisms of granular flow. However, it is difficult to adapt this approach to process modelling because of the limited number of particles which can be handled with the present computational capacity, and the difficulty in handling non-spherical particles. Further work is needed to develoo an aoorooriate aooroach to overcome these problems.     

Forced response of quadratic nonlinear oscillator: comparison of various approaches

The primary resonances of a quadratic nonlinear system under weak and strong external excitations are investigated with the emphasis on the comparison of different analytical approximate approaches. The forced vibration of snap-through mechanism is treated as a quadratic nonlinear oscillator. The Lindstedt-Poincaré method, the multiple-scale method, the averaging method, and the harmonic balance method are used to determine the amplitude-frequency response relationships of the steady-state responses. It is demonstrated that the zeroth-order harmonic components should be accounted in the application of the harmonic balance method. The analytical approximations are compared with the numerical integrations in terms of the frequency response curves and the phase portraits. Supported by the numerical results, the harmonic balance method predicts that the quadratic nonlinearity bends the frequency response curves to the left. If the excitation amplitude is a second-order small quantity of the bookkeeping parameter, the steady-state responses predicted by the second-order approximation of the LindstedtPoincaré method and the multiple-scale method agree qualitatively with the numerical results. It is demonstrated that the quadratic nonlinear system implies softening type nonlinearity for any quadratic nonlinear coefficients.     

Stochastic averaging of quasi-generalized Hamiltonian systems

A stochastic averaging method for generalized Hamiltonian systems (GHS) subject to light dampings and weak stochastic excitations is proposed. First, the GHS are briefly reviewed and classified into five classes, i.e., non-integrable GHS, completely integrable and non-resonant GHS, completely integrable and resonant GHS, partially integrable and non-resonant GHS and partially integrable and resonant GHS. Then, the averaged and FPK equations and the drift and diffusion coefficients for the five classes of quasi-GHS are derived. Finally, the stochastic averaging for a nine-dimensional quasi-partially integrable GHS is given to illustrate the application of the proposed procedure, and the results are confirmed by using those from Monte Carlo simulation.     

粘滞和粘弹性阻尼器减震结构的等效阻尼

利用积分型本构关系,建立了带支撑的一般粘滞和粘弹性阻尼器单自由度耗能结构的微分-积分混合地震响应方程;基于与随机平均分析完全相同的等效准则,推导了可直接应用反应谱的阻尼器的等效刚度和等效阻尼的解析公式;得到了带支撑广义Maxwell阻尼器和广义微分模型阻尼器的等效刚度和等效阻尼的一般结果。通过与一些典型问题的数值计算结果比较,表明了所提方法的有效性。     

凸肩叶片的非线性振动特性与运动分岔

以凸肩叶片作为研究模型, 建立了考虑凸肩摩擦力, 几何大变形与阻尼的非线性振动方程.采用Galerkin法对振动方程离散化, 应用平均法对离散后模态方程组的非线性响应进行解析分析, 得到了非线性幅频特性曲线, 与数值解比较验证了解析解, 并讨论了系统周期解的稳定性. 用非线性振动理论详细研究了平均方程组的运动分岔现象, 揭示了平均方程组周期解的变化过程及其具有的非线性动力学性质. 解析结果表明, 凸肩之间的摩擦对系统第二阶非线性振动特性影响很大. 由于凸肩之间摩擦力方向的不断改变, 系统第二阶非线性幅频特性曲线不连续, 出现两个共振频域. 随着时间的推移, 系统振动的幅值会以$T/ 4$为周期在两个频域的幅频曲线上来回跳动, 这会使叶片的振动响应大幅降低.