Construction of Invariant Torus Using Toeplitz Jacobian Matrices/Fast Fourier Transform Approach

The invariant torus is a very important case in the study of nonlinear autonomous systems governed by ordinary differential equations (ODEs). In this paper a new numerical method is provided to approximate the multi-periodic surface formed by an invariant torus by embedding the governing ODEs onto a set of partial differential equations (PDEs). A new characteristic approach to determine the stability of resultant periodic surface is also developed. A system with two strongly coupled van der Pol oscillators is taken as an illustrative example. The result shows that the Toeplitz Jacobian Matrix/Fast Fourier Transform (TJM/FFT) approach introduced previously is accurate and efficient in this application. The application of the method to normal multi-modes of nonlinear Euler beam is given in [1].     

Exact solutions to drift-flux multiphase flow models through Lie group symmetry analysis

In the present paper, Lie group symmetry method is used to obtain some exact solutions for a hyperbolic system of partial differential equations (PDEs), which governs an isothermal no-slip drift-flux model for multiphase flow problem. Those symmetries are used for the governing system of equations to obtain infinitesimal transformations, which consequently reduces the governing system of PDEs to a system of ODEs. Further, the solutions of the system of ODEs which in turn produces some exact solutions for the PDEs are presented. Finally, the evolutionary behavior of weak discontinuity is discussed.     

Observer design and output feedback stabilization for nonlinear multivariable systems with diffusion PDE-governed sensor dynamics

This paper addresses the problems of observer design and output feedback stabilization for a class of nonlinear multivariable systems, where the nonlinear system dynamics are described by ordinary differential equations (ODEs), and the sensor dynamics are governed by diffusion partial differential equations (PDEs). Based on the Luenberger observer theory, a Luenberger-type PDE-ODE cascaded observer is derived to estimate the state variables of the system. Then, an observer-based output feedback stabilizing controller is developed. The exponential stability of both the observer error system and closed-loop control system is proven via the Lyapunov direct method. Finally, numerical examples are provided to illustrate the effectiveness of the proposed design methods.     

ARC-length method for differential equations

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

Chaotic vibrations of spherical and conical axially symmetric shells

Summary Chaotic vibrations of deterministic, geometrically nonlinear, elastic, spherical and conical axially summetric shells, subject to sign-changing transversal load using the variational principle, are analysed. The paper is motivated by an observation that variational equations of the hybrid type are suitableto solve many dynamical problems of the shells theory. It is assumed that the shell material is isotropic, and the Hook's principle holds. Intertial forces in directions tangent to mean shell surface and rotation inertia of a normal shell cross section are neglected. A transition form PDEs to ODEs (the Cauchy problem) is realized through the Ritz procedure. Next, the Cauchy problem is solved using the fourth-order Runge-Kutta method. Qualitative and quantitative analysis is carried out in the frame of both nonlinear dynamics and quantitative theory of differential equations. New scenarios from harmonic to chaotic dynamics are detected. Various vibration forms development versus control parameters (rise of arc; amplitude and frequency of the exciting force and number of vibrational modes accounted) are illustrated and discussed.     

ADAPTIVE INTERVAL WAVELET PRECISE INTEGRATION METHOD FOR PARTIAL DIFFERENTIAL EQUATIONS

IntroductionThepreciseintegrationmethod(PIM) [1],whichwasproposedforsolvingstructuraldynamicequations.Thismethodissimplerandpossesseshigherprecision .Forlinearsteadystructuraldynamicsystems,itsnumericalresultsattheintegrationpointsarealmostequaltothatoftheexactsolutioninmachineaccuracy .InthepreciseintegrationmethodforsolvingPDEs,theequationsshouldbediscretizedinthephysicalspaceforobtainingthesystemofODEsintime ,whichisoftenexecutedbythefinitedifferencemethodorthefiniteelementmethod .Inrec…     

比例边界等几何分析方法Ⅰ:波导本征问题

提出比例边界等几何方法 (scaled boundary isogeometric analysis, SBIGA), 并用以求解波导本征值问题. 在比例边界等几何坐标变换的基础上, 利用加权余量法将控制偏微分方程进行离散处理, 半弱化为关于边界控制点变量的二阶常微分方程, 即 TE 波或 TM 波波导的比例边界等几何分析的频域方程以及波导动刚度方程, 同时利用连分式求解波导动刚度矩阵. 通过引入辅助变量进一步得出波导本征方程. 该方法只需在求解域的边界上进行等几何离散, 使问题降低一维, 计算工作量大为节约, 并且由于边界的等几何离散, 使得解的精度更高, 进一步节省求解自由度. 以矩形和 L 形波导的本征问题分析为例, 通过与解析解和其他数值方法比较, 结果表明该方法具有精度高、计算工作量小的优点.     

斜拉桥拉索在轴向窄带随机激励下的振动响应

导出了拉索在考虑垂度以及索张力沿索长变化时的参激随机微分方程,进一步给出了预测拉索在窄带随机激励下响应的近似理论解------用统计矩截断法求解矩方程,获得高斯闭合解和一阶非高斯闭合解. 以南京长江二桥约330米长的A20拉索为研究对象,对以上高斯闭合解和一阶非高斯闭合解进一步进行数值求解以获得拉索的响应,并采用Monte-Carlo数值方法对求解进行验证. 分析了拉索振动的一般特征,特别分析了激励中心频率和拉索频率比为1和2时的响应随激励带宽的变化特征,得到了一些新的结论.      

Neural network as a function approximator and its application in solving differential equations

A neural network(NN) is a powerful tool for approximating bounded continuous functions in machine learning. The NN provides a framework for numerically solving ordinary differential equations(ODEs) and partial differential equations(PDEs)combined with the automatic differentiation(AD) technique. In this work, we explore the use of NN for the function approximation and propose a universal solver for ODEs and PDEs. The solver is tested for initial value problems and boundary value problems of ODEs, and the results exhibit high accuracy for not only the unknown functions but also their derivatives. The same strategy can be used to construct a PDE solver based on collocation points instead of a mesh, which is tested with the Burgers equation and the heat equation(i.e., the Laplace equation).     

Simulation of electrically driven jet using Chebyshev collocation method

The model of electrically driven jet is governed by a series of quasi 1D dimensionless partial differential equations (PDEs). Following the method of lines, the Chebyshev collocation method is employed to discretize the PDEs and obtain a system of differential-algebraic equations (DAEs). By differentiating constrains in DAEs twice, the system is transformed into a set of ordinary differential equations (ODEs) with invariants. Then the implicit differential equations solver “ddaskr” is used to solve the ODEs and post-stabilization is executed at the end of each step. Results show the distributions of radius, linear charge density, stretching ratio and also the horizontal velocity at a time point. Meanwhile, the spiral and expanding projections to X-Y plane of the jet centerline suggest the occurring of bending instability.     

Stationary random waves propagation in 3D viscoelastic stratified solid

Propagation of stationary random waves in viscoelastic stratified transverse isotropic materials is investigated. The solid was considered multi-layered and located above the bedrock, which was assumed to be much stiffer than the soil, and the power spectrum density of the stationary random excitation was given at the bedrock. The governing differential equations are derived in frequency and wave-number domains and only a set of ordinary differential equations ( ODEs) must be solved. The precise integration algorithm of two-point boundary value problem was applied to solve the ODEs. Thereafter, the recently developed pseudo-excitation method for structural random vibration is extended to the solution of the stratified solid responses.     

Reduced-order models for nonlinear vibrations of fluid-filled circular cylindrical shells: Comparison of POD and asymptotic nonlinear normal modes methods

The aim of the present paper is to compare two different methods available for reducing the complicated dynamics exhibited by large amplitude, geometrically nonlinear vibrations of a thin shell. The two methods are: the proper orthogonal decomposition (POD), and an asymptotic approximation of the nonlinear normal modes (NNMs) of the system. The structure used to perform comparisons is a water-filled, simply supported circular cylindrical shell subjected to harmonic excitation in the spectral neighbourhood of the fundamental natural frequency. A reference solution is obtained by discretizing the partial differential equations (PDEs) of motion with a Galerkin expansion containing 16 eigenmodes. The POD model is built by using responses computed with the Galerkin model; the NNM model is built by using the discretized equations of motion obtained with the Galerkin method, and taking into account also the transformation of damping terms. Both the POD and NNMs allow to reduce significantly the dimension of the original Galerkin model. The computed nonlinear responses are compared in order to verify the accuracy and the limits of these two methods. For vibration amplitudes equal to 1.5 times the shell thickness, the two methods give very close results to the original Galerkin model. By increasing the excitation and vibration amplitude, significant differences are observed and discussed. The response is investigated also for a fixed excitation frequency by using the excitation amplitude as bifurcation parameter for a wide range of variation. Bifurcation diagrams of Poincaré maps obtained from direct time integration and calculation of the maximum Lyapunov exponent have been used to characterize the system.     

环形桁架结构径向振动的等效圆环模型

在大型环形网架式可展天线中,环形桁架结构的动力学性能对于整个天线的工作状态至关重要.针对大型空间桁架结构,基于连续体等效的思想,将其动力学模型简化为简单的弹性连续体模型一直是动力学研究的热点.将环形桁架结构看作由重复的平面桁架单元构成的环形周期结构,在周期桁架单元等效梁模型的基础上,提出采用不计剪切变形和转动惯量的等效圆环模型分析环形桁架结构的径向振动,并对等效圆环模型的偏微分运动方程进行了解析求解.首先通过变量代换将描述圆环径向振动的四阶偏微分方程组降阶为一阶偏微分方程组,然后通过对降阶后的偏微分方程组进行Laplace变换将其转化为常微分方程组,并采用微分方程组的Green函数解法,获得了等效圆环模型在复频域下动力响应的解析表达式,进而得到等效圆环模型固有振动的特征方程及传递函数的表达式.最后通过数值算例对环形桁架有限元模型与等效圆环模型的固有频率和振型以及传递函数进行了对比分析,验证了等效圆环模型用于环形桁架结构径向振动分析的可行性.     

Analysis of in-plane 1:1:1 internal resonance of a double cable-stayed shallow arch model with cables' external excitations

The nonlinear dynamic behaviors of a double cable-stayed shallow arch model are investigated under the one-to-one-to-one internal resonance among the lowest modes of cables and the shallow arch and external primary resonance of cables. The in-plane governing equations of the system are obtained when the harmonic excitation is applied to cables. The excitation mechanism due to the angle-variation of cable tension during motion is newly introduced. Galerkin's method and the multi-scale method are used to obtain ordinary differential equations(ODEs) of the system and their modulation equations, respectively. Frequency-and force-response curves are used to explore dynamic behaviors of the system when harmonic excitations are symmetrically and asymmetrically applied to cables. More importantly, comparisons of frequency-response curves of the system obtained by two types of trial functions, namely, a common sine function and an exact piecewise function, of the shallow arch in Galerkin's integration are conducted.The analysis shows that the two results have a slight difference; however, they both have sufficient accuracy to solve the proposed dynamic system.     

Heat transfer characteristics of thin power-law liquid films over horizontal stretching sheet with internal heating and variable thermal coefficient

The effect of internal heating source on the film momentum and thermal transport characteristic of thin finite power-law liquids over an accelerating unsteady horizontal stretched interface is studied. Unlike most classical works in this field, a general surface temperature distribution of the liquid film and the generalized Fourier’s law for varying thermal conductivity are taken into consideration. Appropriate similarity transformations are used to convert the strongly nonlinear governing partial differential equations (PDEs) into a boundary value problem with a group of two-point ordinary differential equations (ODEs). The correspondence between the liquid film thickness and the unsteadiness parameter is derived with the BVP4C program in MATLAB. Numerical solutions to the self-similarity ODEs are obtained using the shooting technique combined with a Runge-Kutta iteration program and Newton’s scheme. The effects of the involved physical parameters on the fluid’s horizontal velocity and temperature distribution are presented and discussed.     

斜拉桥塔-索-桥面耦合参数振动模型及响应分析

针对大跨度斜拉桥拉索与桥塔、桥面的协同振动问题,考虑拉索垂度、阻尼、倾角以及重力弦向分力的影响,引入拉索高精度抛物线形,建立了桥塔-索-桥面连续非线性精细化振动模型,推导了桥塔和桥面共同激励作用下斜拉索耦合振动方程,对比分析了2种激振模式下斜拉索的参数振动特性,并编制程序研究了桥面与拉索的频率比、桥面激励幅值、索力及阻尼对结构耦合振动特性的影响规律。结果表明：桥面与拉索频率比对系统振动的影响较大,频率比为1:2和2:1时拉索均产生强烈振动,但2:1激振模式下拉索振幅更大,达到共振时间较长;随着桥面激励幅值的增大,2:1亚谐波共振模式下的拉索振幅增长速率更快;拉索振幅随索力的增大呈非线性减小趋势;斜拉索阻尼超过2%时,继续提高自身阻尼不能有效减小其振动幅值,需要通过设置附加阻尼才能更好地抑制其振动。     

行进绳索在横向流体激励下的运动

给出横向流体对行进绳索的作用力描述,建立了绳索的动力学方程。由于该方程具有零刚度特征,引入Pilipchuk变换,以自平衡状态起度量的径向振动和回转运动来描述绳的运动,获得非零刚度系统。然后,用两变量参数摄动法求得绳索关于扩张振动和回转运动的约化型,使得绳索运动可近似由一个二维动力系统来描述。最后,用数值方法讨论了行进速度和重力对绳索运动形态的影响。     

Postbuckling and mode jumping analysis of composite laminates using an asymptotically correct, geometrically non-linear theory

An analytic method is presented in this paper to study the postbuckling and mode jumping behavior of bi-axially compressed composite laminates. The governing partial differential equations (PDEs) are derived rigorously from an asymptotically correct, geometrically non-linear theory. A novel and relatively simpler solution approach is developed to solve the two coupled fourth-order PDEs, namely, the compatibility equation and the dynamic governing equation. The generalized Galerkin method is used to solve boundary value problems corresponding to antisymmetric angle-ply and cross-ply composite plates, respectively. The variety of possible modal interactions is expressed in an explicit and concise form by transforming the coupled non-linear governing equations into a system of non-linear ordinary differential equations (ODEs).

The comparison between the present method and the finite element analysis (FEA) shows a pretty good match in their numerical results in the primary postbuckling region. While the FEA may lose its convergence when solution comes close to the secondary bifurcation point, the analytic approach has the capability of exploring deeply into the post-secondary buckling realm and capture the mode jumping phenomenon for various combinations of plate configurations and in-plane boundary conditions. Free vibration along the stable primary postbuckling and the jumped equilibrium paths are also studied.     

TWO-DIMENSIONAL VORTEX-INDUCED VIBRATION OF CABLE SUSPENSIONS

Resonant responses of suspended elastic cables driven by a steady current are investigated. Phenomenological fluid force models for alternate vortex-shedding are coupled with the nonlinear partial differential equations of cable motion. Decoupled cross-flow and in-line vortex-induced vibrations (VIV) are examined first using linearized and nonlinear cable models. The linearized cable model predicts well the basic characteristics of VIV and the nonlinear cable model captures the hysteresis often observed in experiments. Next, coupled cross-flow and in-line vibrations are evaluated by considering two principal coupling mechanisms: (i) cable structural nonlinearities, and (ii) coupled fluid lift and drag. Attention is focused on a “worst-case” resonant response where the natural frequencies for cable modes in the cross-flow and in-line directions are in the same 1:2 ratio as the excitation frequencies associated with lift and drag. The inclusion of cable structural nonlinearities alone leads to coupled responses that differ qualitatively (i.e., in number and stability of periodic motions) when compared to those of the decoupled model. The inclusion of coupled fluid lift and drag produces non-planar “figure eight” motions of the cable cross-section that exhibit similar characteristics to those previously measured on spring supported cylinders.     

斜拉桥拉索-阻尼器系统非线性响应分析

考虑索的抗弯刚度、垂度及几何非线性的影响，得出了索一阻尼器系统的空间非线性振动偏微分方程，用中心差分法将微分方程在空间内离散，导出了系统的非线性振动常微分方程组。结合Newmark法及虚拟力法提出了一种用于求解非线性振动瞬态响应的杂交分析算法。并以典型的斜拉桥拉索为研究对象，给出了数值算例，并与Runge—Kutta直接积分法进行了比较，说明了杂交算法的准确性及有效性。