Nonlinear Oscillation Analysis by an Orthogonal Function Method

ＮＯＮＬＩＮＥＡＲＯＳＣＩＬＩＡＴＩＯＮＡＮＡＬＹＳＩＳＢＹＡＮＯＲＴＨＯＧＯＮＡＬＦＵＮＣＴＩＯＮＭＥＴＨＯＤＳｕｎＰｉ，ｈｏｎｇ（孙丕忠）ＴａｎｇＰｉａｎｇａｎｇ（唐乾刚）ＳｕｎＳｈｉｘｉａｎ（孙世贤）（ＮａｔｉｏｎＵｎｉｖｅｒｓｉｔｙｏｆＤｅｆ...     

Nonlinear time transformation method for strong nonlinear oscillation systems

In this paper, a nonlinear time transformation method is presented for the analysis of strong nonlinear oscillation systems. This method can be used to study the limit cycle behavior of the autonomous systems and to analyze the forced vibration of a strong nonlinear system. The project partly supported by the Foundation of Zhongshan University Advanced Research Center     

离散正交法(DOM)离散点的正交分析与程序设计

对能够求解一系列线性常微分方程组边值问题的数值计算方法———离散正交法(DOM)进行了离散点的正交分析,给出了计算机实现数值计算的程序设计原理与计算流程图,指出了该方法能够克服传统计算方法由于所求函数的快速增长所引起的边界效应和局部效应的缺点,给出了得到稳定计算过程的可能性。为了推广应用,文中介绍了离散正交法的基本原理、实现方法和求解过程,讨论了采用离散正交法来求解非线性问题的处理方法。并且以承受均布载荷的环形板为例,将采用离散正交法的计算结果与经典解作了对比。     

Two-phase flow analysis based on a phase-field model using orthogonal basis bubble function finite element method

In this study, a finite element method based on a phase-field model for gas–liquid two-phase flow is proposed. MINI element based on a bubble function element stabilisation method is employed for the incompressible Navier–Stokes equations. The Cahn–Hilliard equation is employed to estimate the interface of gas and liquid. The orthogonal basis bubble function element is used to solve the Cahn–Hilliard equation. In particular, a detailed explanation for solving the Cahn–Hilliard equation based on a finite element method is given.     

Nonlinear oscillations with parametric excitation solved by homotopy analysis method

An analytical technique, namely the homotopy analysis method （HAM）, is used to solve problems of nonlinear oscillations with parametric excitation. Unlike perturbation methods, HAM is not dependent on any small physical parameters at all, and thus valid for both weakly and strongly nonlinear problems. In addition, HAM is different from all other analytic techniques in providing a simple way to adjust and control convergence region of the series solution by means of an auxiliary parameter h. In the present paper, a periodic analytic approximations for nonlinear oscillations with parametric excitation are obtained by using HAM, and the results are validated by numerical simulations.     

Nonlinear aeroelastic analysis of an airfoil-store system with a freeplay by precise integration method

The aeroelastic system of an airfoil-store configuration with a pitch freeplay is investigated using the precise integration method (PIM). According to the piecewise feature, the system is divided into three linear sub-systems. The sub-systems are separated by switching points related to the freeplay nonlinearity. The PIM is then employed to solve the sub-systems one by one. During the solution procedures, one challenge arises when determining the vibration state passing the switching points. A predictor-corrector algorithm is proposed based on the PIM to tackle this computational obstacle. Compared with exact solutions, the PIM can provide solutions to the precision in the order of magnitude of 10⁻¹². Given the same step length, the PIM results are much more accurate than those of the Runge–Kutta (RK) method. Moreover, the RK method might falsely track limit cycle oscillations (LCOs), bifurcation charts or chaotic attractors; even the step length is chosen much smaller than that for the PIM. Bifurcations and LCOs are obtained and analyzed by the PIM in detail. Interestingly, it is found that multiple LCOs and chaotic attractors can exist simultaneously. With this magnitude of precision and efficiency, the PIM could become a solution technique with excellent potential for piecewise nonlinear aeroelastic systems.     

局部正交无网格伽辽金法的研究及其在含多裂纹多孔结构中的应用

通过对无网格法中正交基函数的研究,提出局部正交无网格伽辽金法,局部正交基函数保持原正交基函数的性质,但其导数具有了通式,简洁明了,易于编程实现,计算效率高,并将其应用到求解含多裂纹多孔均匀拉伸板的应力强度因子中,计算结果与用正交无网格伽辽金法和有限元法得到的结果进行比较,证明了局部正交无网格伽辽金法的可行性和正确性。     

基于正交基无单元Galerkin法和非线性规划的安定分析方法

基于安定分析的下限定理,用正交基无单元Galerkin法建立了交交载荷作用下理想弹塑性结构安定分析的下限计算格式.在给定载荷域的载荷角点所对应的载荷作用下,采用正交基无单元Galerkin法计算相应的虚拟弹性应力场.并且利用结构在正交基无单元Galerkin法弹塑性增量分析中平衡迭代结果计算得到自平衡应力场基矢量,然后由这些基矢量的线性组合模拟自平街应力场.安定分析问题最终被归结为一系列未知变量较少的非线性数学规划子问题,通过复合形法求解.算例表明该方法的计算结果是令人满意的,并且对初始复合形顶点和用于构造自平衡应力场基矢量的载荷增量是非常不敏感的.     

Nonlinear vibration analysis by an extended averaged equation approach

The paper presents an extended averaged equation approach to the investigation of nonlinear vibration problems. The proposed method is applied to some free and self-excited oscillators, the Duffing's forced oscillators including main resonance, subharmonic resonance and super harmonic resonance. The results in analyzing the vibration systems with arbitrary non-linearity show advantages of the method.     

Transonic limit cycle oscillation behavior of an aeroelastic airfoil with free-play

The limit cycle oscillation (LCO) behaviors of an aeroelastic airfoil with free-play for different Mach numbers are studied. Euler equations are adopted to obtain the unsteady aerodynamic forces. Aerodynamic and structural describing functions are employed to deal with aerodynamic and structural nonlinearities, respectively. Then the flutter speed and flutter frequency are obtained by V-g method. The LCO solutions for the aeroelastic airfoil obtained by using dynamically linear aerodynamics agree well with those obtained directly by using nonlinear aerodynamics. Subsequently, the dynamically linear aerodynamics is assumed, and results show that the LCOs behave variously in different Mach number ranges. A subcritical bifurcation, consisting of both stable and unstable branches, is firstly observed in subsonic and high subsonic regime. Then in a narrow Mach number range, the unstable LCOs with small amplitudes turn to be stable ones dominated by the single degree of freedom flutter. Meanwhile, these LCOs can persist down to very low flutter speeds. When the Mach number is increased further, the stable branch turns back to be unstable. To address the reason of the stability variation for different Mach numbers at small amplitude LCOs, we find that the Mach number freeze phenomenon provides a physics-based explanation and the phase reversal of the aerodynamic forces will trigger the single degree of freedom flutter in the narrow Mach number range between the low and high Mach numbers of the chimney region. The high Mach number can be predicted by the freeze Mach number, and the low one can be estimated by the Mach number at which the aerodynamic center of the airfoil lies near its elastic axis. Influence of angle of attack and viscous effects on the LCO behavior is also discussed.     

粒子法中压力振荡问题的研究

传统移动粒子半隐式法MPS(Moving Particle Semi-implicit Method)中一直存在压力振荡问题,针对此问题对MPS方法进行改进。改进的MPS方法,采用一种新型抑制压力振荡的压力泊松方程离散格式;在核函数的选择方面,采用能够增加计算稳定性的二次样条核函数;并且针对MPS方法中粒子插值不完整问题,对粒子插值不完整性进行了修正。应用改进的MPS方法对溃坝问题进行数值模拟验证。结果表明,应用改进的MPS方法能够得到更为光滑的压力场空间分布。对模拟过程中的检测点压力进行采集,并且与实验值进行对比分析,发现改进的MPS方法能够有效地抑制模拟过程中的压力振荡,而且与实验值接近。同时应用改进的MPS方法对静水问题进行验证模拟,发现改进的MPS方法能够有效地抑制模拟过程中的压力振荡,而且监测点的压力与理论解接近。改进的MPS方法对今后应用MPS方法模拟实际工程问题,并且获得准确稳定的压力值有着重要的意义。     

Nonlinear quasi-conforming finite element method

The nonlinear quasi-conforming FEM is presented based on the basic concept of the quasiconforming finite element. First, the incremental principle of stationary potential energy is discussed. Then, the formulation process of the nonlinear quasi-conforming FEM is given. Lastly, two computational examples of shells are given.     

NONLINEAR ANALYSIS OF EQUATORIAL EASTERN PACIFIC AIR-SEA COUPLING OSCILLATION AND A LIMIT- CYCLE THEORY FOR ENSO CYCLE

The troposphere and ocean mixed layer were considered as two components of a dynamic system operated by solar radiation as the constant source of energy, where upon an air-sea coupling self-exited coupling oscillation model was based with the aid of a locally averaged thermodynamic climate model, resulting mathematically in a closed self-governed dynamic system, a so-called El Nino-Southern Oscillation(ENSO) system. With the limit cycle solution of the system. It is shown that the essential physics of the coupled system can be described by the ENSO system. Compared with the observations, the theoretical limit cycle orbit matches the observed phase loop qualitatively. The ENSO system provides a useful theoretical framework for study of interannual variation of the tropical climate system.     

自激振荡脉冲射流喷嘴装置系统频率特性理论研究

根据相似系统原理和流体网络理论建立了自激振荡脉冲射流喷嘴装置的等效网络模型,利用系统传递函数推导了系统频率特性方程并进行了数值计算。结果表明:喷嘴装置的固有频率主要由喷嘴形状、结构参数、入口流速、射流中压力扰动波波速决定;自激振荡腔腔径、自激振荡腔腔长、上喷嘴直径、下喷嘴直径都对系统频率特性影响很大。提出了相应的自激振荡脉冲射流喷嘴设计准则,即喷嘴装置在最佳阻尼比下产生谐波共振。     

板弯曲问题的群论方法

利用群论方法提出周期区域的分片正交多项式连续函数,在周期区域内利用正交分片多项式逼近位移函数可以大大地降低计算量。非周期区域的板可以通过附加边界载荷将其扩展成周期区域的板。利用此方法可以大大降低板弯曲问题的计算量。     

Improved L-P method for solving strongly nonlinear problems

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｓｔｒｏｎｇｌｙｎｏｎｌｉｎｅａｒｏｓｃｉｌｌａｔｉｏｎｐｒｏｂｌｅｍ¨ｕ ω20 ｕ =εｆ(ｕ , ｕ ,¨ｕ) , (1 )ｗｈｅｒｅｔｈｅｓｙｍｂｏｌ“·”ｄｅｎｏｔｅｓｔｈｅｄｉｆｆｅｒｅｎｔｉａｔｉｏｎｗｉｔｈｒｅｓｐｅｃｔｔｏｔｈｅｉｎｄｅｐｅｎｄｅｎｔｖａｒｉａｂｌｅｔａｎｄω0ａｎｄεａｒｅａｌｌｔｈｅａｒｂｉｔｒａｒｙｃｏｎｓｔａｎｔｓ,ｈａｖｅｉｎｔｈｅｒｅｃｅｎｔｙｅａｒｓｓｅｖｅｒａｌｉｍｐｒｏｖｅｄＬ＿Ｐｍｅｔｈｏｄｓｗｈｅｎεｉｓｎｏｔｓｍａｌｌ (ｓｔｒｏ…     

Nonlinear numerical simulation method for galloping of iced conductor

Based on the principle of virtual work, an updated Lagrangian finite element formulation for the geometrical large deformation analysis of galloping of the iced conductor in an overhead transmission line is developed. In numerical simulation, a three-node isoparametric cable element with three translational and one torsional degrees-of-freedom at each node is used to discretize the transmission line. The nonlinear dynamic system equation is solved with the Newmark time integration method and the Newton-Raphson nonlinear iteration. Numerical examples demonstrate the efficiency of the presented method and the developed finite element program. A new possible galloping mode, which may reflect the saturation phenomenon of a nonlinear dynamic system, is discovered under the condition that the lowest order of vertical natural frequency of the transmission line is approximately two times of the horizontal one.     

A New Method for Approximate Analytical Solutions to Nonlinear Oscillations of Nonnatural Systems

This paper deals with nonlinear oscillations of a conservative,nonnatural, single-degree-of-freedom system with odd nonlinearity. Bycombining the linearization of the governing equation with the method ofharmonic balance, we establish approximate analytical solutions for thenonlinear oscillations of the system. Unlike the classical harmonicbalance method, the linearization is performed prior to proceeding withharmonic balancing thus resulting in linear algebraic equations insteadof nonlinear algebraic equations. Hence, we are able to establish theapproximate analytical formulas for the exact period and periodicsolution. These approximate solutions are valid for small as well aslarge amplitudes of oscillation. Two examples are presented toillustrate that the proposed formulas can give excellent approximateresults.     

NONLINEAR THREE-DIMENSIONAL ANALYSIS OF COMPOSITELAMINATED PLATES

ＮＯＮＬＩＮＥＡＲＴＨＲＥＥ－ＤＩＭＥＮＳＩＯＮＡＬＡＮＡＬＹＳＩＳＯＦＣＯＭＰＯＳＩＴＥＬＡＭＩＮＡＴＥＤＰＬＡＴＥＳ￥(江晓禹，张相周，陈百屏)ＪｉａｎｇＸｉａｏｙｕ；（ＳｏｕｔｈｗｅｓｔｅｒｎＪｉａｏｔｏｎｇＵｎｉｖｅｒｓｉｔｙ，Ｃｈｅｎｇｄｕ６...     

Pattern analysis of viscous oscillation in a shear flow induced by an exothermic reaction

An investigation is made of oscillatory phenomenon induced by an exothermic reaction. This oscillatory phenomenon occurs in a very thin mixing layer between two miscible and reacting fluids. A qualitative model based on the interaction among reaction, molecular momentum transfer, molecular heat transfer, molecular mass transfer and forced convective transfer was proposed by Kuroda and Ogawa [1994. Nonlinear waves in a shear flow with a diffusive exothermic reaction and its qualitative reasoning. Chem. Eng. Sci. 49(16), 2699-2708]. It is experimentally shown in the present study that the oscillatory patterns change in the flow direction. In the upstream area where oscillatory patterns are nearly straight stripes, effects of viscosity on those stripes are investigated. In the downstream area where stripes are wavy and disordered, fractal analysis is introduced to investigate the relationship between the transition of oscillatory flow patterns and process factors, i.e. the viscosity ratio, the entire viscosity, flow rate and the flow rate ratio. Fractal analysis is also applied to temperature oscillation, and it is confirmed that the characteristic patterns of oscillation become obscure as the entire viscosity increases. The entire viscosity is an important factor for controlling both oscillatory patterns of flow and temperature in this reactive flow system.