Van der Pol oscillator with time variable parameters

In this paper, the analytical solving procedure of the oscillator with slow time variable mass is developed. The solution is based on the Jacobi elliptic function whose properties: frequency, amplitude and modulus are obtained according to the requirements given for the amplitude of the displacement and the amplitude of the velocity of vibration and also period of vibration. The suggested procedure is applied for the solution of the time variable Van der Pol oscillator. The limit value of the initial mass of the oscillator is determined which separates the case when the limit cycle motion occurs, and the case when the amplitude of vibration tends to zero independently of the initial displacement. A numerical example is considered. The analytical solution is compared with the numerically obtained one and it is concluded that they are in good agreement.     

带三次恢复力项频率依赖于速度的非线性振子研究

带三次恢复力项频率依赖于速度（Velocity-Dependent-Frequency, VDF）的非线性振子 的周期及其性质目前没有文献讨论,且使用传统的摄动法或谐波平衡法求解这类振子一阶近似解时往往失效。特别的,其频率在有限的幅值范围内奇异。首先求得了该振子周期的积分表达式,基于积分表达式可积性条件采用谐波平衡法得到了该振子一类初始条件下的精确解;研究了该振子的周期性质,给出了由第一类完全椭圆积分表示的周期-振幅的近似解析表达式,分析了振子的方波现象及产生原因。研究表明,振子周期最终随着幅值的增大衰减到0;振子方波现象产生原因是由于系统参数 ,随着幅值的增大,方波现象更明显。最后提出使用一种Hermite插值法求解该振子的周期解,该方法将时间变量转换为新的谐振时间变量,其频率为振子频率的一半,对应的控制微分方程转变为适合于Hermite插值分析的形式,其解与数值解的对比证明了该方法的有效性。     

退化环面上的非线性jerk方程近似周期解的同伦分析方法

用同伦分析法求解退化环面上的非线性Jerk方程的近似周期和近似解析周期解。所得结果表明文中得到的一阶近似周期和一阶近似解析周期解与Gottlieb用低阶谐波平衡法求解得到的结果一样。当参数和初速度较大时,一阶近似周期与精确周期的百分比误差是4.831 8%,而二阶近似周期与精确周期的百分比误差小于0.219 9%。与数值方法给出的"精确"周期解比较,二阶近似解析周期解比一阶近似解析周期解要精确的多。因此,同伦分析法是求解非线性Jerk方程的一种非常有效的方法。     

达芬-谐波振子的改进解析逼近解

研究达芬-谐波振子的解析逼近解。所谓达芬-谐波振子是指当位移远小于1时，系统可化为三次非线性振子，而当位移远大于1时，该系统则化为线性谐波振子。通过将变形后的控制方程的线性化与谐波平衡法组合起来，我们建立了达芬-谐波振子频率及周期解的改进解板逼近。改进的解析逼近在振幅的全部取值范围内，包括振幅趋于零及无穷的情况，都有令人非常满意的精度。     

Modified higher-order wake oscillator model for vortex-induced vibration of circular cylinders

In the present study, a modified model is introduced to estimate the structural oscillation amplitude of a circular cylinder during lock-in in the vortex-induced vibration phenomenon. The modified model includes the same van der Pol term as in the classic model, while its coefficients are assumed to be variable. This modified model is utilized to bring about compliance between theoretical solution and experimental results. Then, a fifth-order aerodynamic term is added, and the coefficients are modified and optimized using a new straightforward solution method. Here, the displacement, velocity and acceleration coupling terms are used to solve the system of equations. The response of the coupled equations is assumed to be harmonic. A linear approach is adopted to simplify and derive the solutions algebraically. The results are plotted during lock-in for amplitude as a function of reduced velocity and for maximum structural oscillation amplitude versus Skop–Griffin parameter. Finally, these results are compared with those of the classic wake oscillator model. The present modified model evinces an exact compliance with experimental measurements regarding structural oscillation amplitude, lock-in range, some hysteresis and mass damping ratio attributes, and reduces sensitivity to the type of coupling term through adding the fifth-order damping term.     

Investigation of vertical vibrations of rigid disk on saturated poroelastic stratum under high-frequency excitation

The forced high-frequency vertical vibrations of a rigid circular disk on a saturated poroelastic stratum supported by a rigid rock foundation are studied. An exact formulation is developed for the dual integral equations governing the mixed boundary-value problem of the vertical vibration. By finding an approximate solution to the exact governing dual integral equations, an evaluation of the amplitude response of the rigid disk is presented for the case of a frictionless contact both between the disk and the stratum and between the stratum and the foundation. The influences of parameters regarding the disk and the elastic stratum on the amplitude response of the disk are discussed. It is concluded that the amplitude response of the disk is dependent on the exciting frequency, the depth of the stratum, the Poison’s ratio of the solid skeleton, the porosity of the medium and the radius and mass of the disk.     

非线性奇异振子的解析逼近解

利用修正的牛顿-谐波平衡法建立了非线性奇异振子的解析逼近周期和周期解．通过改写控制方程和选取简单、合适的校正项对牛顿-谐波平衡法进行了修正.构造的两个解析逼近周期和周期解不仅在振幅和参数全部取值范围内有效且能快速地收敛到精确解；两个逼近周期与精确周期的百分比误差分别低于0.92%和0.09%，后者比已有结果精度高。     

Optimization-oriented exponential-polynomial-closure approach for analyzing nonlinear stochastic oscillators

A novel method named optimization-oriented exponential-polynomial-closure approach is proposed in this article. The main idea of this attempt is to extend the original exponential-polynomial-closure solution procedure methodologically by minimizing the resulted residual error square of the governing equation, which is achieved after an exponential polynomial is adopted as the approximate solution. The objective function for computing the parameters in the approximate solutions of nonlinear random oscillators is then formulated. The probabilistic solutions of the oscillators obtained by the presented approach are verified by the exact solutions in some special cases or by Monte Carlo simulation. Numerical examples indicate that the solutions attained by the presented approach match with the exact or Monte Carlo simulation solutions. The advantage of the presented solution procedure is that it can provide a much more accurate solution than the Equivalent Linearization approach and it is much more efficient than Monte Carlo simulation as demonstrated by the numerical examples.     

传导薄板的非线性磁弹性振动问题

研究了磁场环境中传导薄板的非线性磁弹性振动问题。由虚功原理,给出了磁场中薄板的磁弹性耦合运动方程。并根据薄板薄壳的磁弹性基本假设及麦克斯威尔方程,得到了方程中电磁力及力矩的表达式。作为具体问题,采用多尺度法求出了横向磁场中条形板非线性振动的近似解析解。通过算例,分析了磁场环境对振动周期和幅值的影响。从而证实了,通过改变磁场因素,可达到控制该磁场环境中传导薄板振动的目的。所得的结论,对工程实际将有较大的应用参考价值。     

Approximate analytic solutions of adiabatic capillary tube

Approximate analytic solution of capillary tube is valuable for theoretical analysis and engineering calculation. In this work, two kinds of approximate analytic solutions of adiabatic capillary tube have been developed. One is the explicit function of capillary tube length. Another is the explicit function of refrigerant mass flow rate. In these solutions, the choked flow condition is taken into account without iterative calculations. The approximate predictions are found to agree reasonably well with experimental data in open literatures.     

应用优化的同伦分析法求解非线性Jerk方程

应用优化的同伦分析法计算了具有三次非线性项的三阶微分方程(Jerk)的近似周期和近似解析周期解。文中给出一个算例说明由优化的同伦分析法可以容易得到精确的二阶近似周期解。当初速度 比较大时,一阶近似周期与精确周期的百分比误差为-0.415%,而二阶近似周期与精确周期的百分比误差为-0.0298%。与数值方法给出的“精确”周期解比较,一阶近似解析周期解和二阶近似周期解的精度很高。这个说明同伦分析法对求解非线性Jerk方程非常有效     

圆柱壳体的非线性磁弹性振动问题

导出了纵向磁场中圆柱薄壳体的非线性轴对称振动基本方程。利用多尺度法和摄动法求出了两端简支壳体振动的近似解析解,并对振动幅值和频率进行了讨论。     

Forced vibrations of a body supported by viscohyperelastic shear mountings

The damped, finite-amplitude forced vibration of a rigid body supported symmetrically by simple shear springs and by a smooth inclined bearing surface is studied. The spring material is characterized as a compressible or incompressible, homogeneous and isotropic viscohyperelastic material for which the shear response function in a simple shear deformation is a quadratic function of the amount of shear. The trivial case of constant shear response is included. The equation for the damped motion of the load is a nonlinear, ordinary differential equation of the forced Duffing type with a constant static shift term due to gravity, and for which an exact solution is unknown. An approximate solution is obtained by the method of harmonic balance. Results for the motion of the load relate the system design parameters to the amplitude-frequency response and to the amplitude-driving force intensity response of the system. Regions of stable motion are identified in terms of the amplitude of the motion, driving-force intensity, driving frequency, and system design parameters. Geometrical characterizations of the motion are related schematically to certain cross-sections through the full three-dimensional solution surfaces for the amplitude and for the phase of the motion. A simple diagram maps the loci of all bifurcation points against the static shear deflection, which serves as the system design parameter for the inclined motion. An infinitesimal stability analysis shows that the bifurcation points of the inclined motion fall on the stability boundaries of the numerical solution of a three-parameter Hill equation. The solution provides information that illustrates how the system design parameters affect the motion of the load and how these may be chosen to control the amplitude of the oscillations and the stability of the system. The results are valid for all compressible or incompressible, homogeneous and isotropic, viscohyperelastic materials in the aforementioned class.     

随机参数结构最优控制的闭环响应分析

在模态坐标下对结构降阶进行最优控制,用近似离散化的方法得到结构闭环响应的近似解来代替精确解。考虑结构中物理参数和几何参数的随机性,将这些随机参数用随机因子来表示。在此基础上,利用求解随机变量函数矩的方法,导出了在最优控制的作用下,随机结构位移闭环响应的均值和方差的计算表达式。通过算例考察了结构各个参数的随机性对闭环响应的影响,经与Monte Carlo数值模拟法结果比较,验证了文中理论分析和计算方法的正确性。     

非线性流滞阻尼器耗能结构随机地震响应和首超时间分析

对非线性流滞阻尼器耗能结构在Kanai-Tajimi谱地震激励下的随机响应及其随机失效时间和动力可靠性进行了系统研究。首先建立了结构的非线性运动方程;然后,基于随机平均法,将结构响应幅值近似为一维markov扩散过程,获得了扩散过程漂移系数和扩散系数的解析表达式;其次,利用扩散过程与FPK方程的对应关系,获得了幅值平稳概率密度函数和幅值任意阶矩的解析表达式;再次,利用幅值与结构位移和速度的相互转化关系,获得了结构位移与速度的平稳联合概率密度函数和位移、速度方差以及位移期望穿越率的解析表达式;最后,利用扩散过程的后向Kolmogrov方程,基于首超失效模型,建立了结构动力可靠性函数方程和结构随机失效时间统计矩方程,并利用一维扩散过程的边界分类性质,将统计矩方程的奇异定性边界条件转化为等价的定量边界条件,进而获得了失效时间任意阶统计矩的解析解,并利用此矩,对结构动力可靠性和失效时间概率分布函数进行了近似分析,给出了算例,从而建立了结构非线性随机地震响应及其随机失效时间和动力可靠性的分析方法。     

A variational approach to the problem of solidification of a metal semi-infinite thermally nonlinear body

Summary Based upon the Gauss differential variational principle and the method of optimal linearization, analytical approximate solutions are developed for the problem of solidification of a metal semi-infinite body whose coefficient of thermal conductivity depends linearly on the temperature. For the case of constant thermal properties, the approximate solution obtained by means of the Gauss variational principle agrees remarkable well with the exact solution obtained by Neumann.With 1 Figure     

Non-linear vibration of frames by the incremental dynamic stiffness method

The dynamic stiffness method is extended to large amplitude free and forced vibrations of frames. When the steady state vibration is concerned, the time variable is replaced by the frequency parameter in the Fourier series sense and the governing partial differential equations are replaced by a set of ordinary differential equations in the spatial variables alone. The frequency-dependent shape functons are generated approximately for the spatial discretization. These shape functions are the exact solutions of a beam element subjected to mono-frequency excitation and constant axial force to minimize the spatial discretization errors. The system of ordinary differential equations is replaced by a system of non-linear algebraic equations with the Fourier coefficients of the nodal displacements as unknowns. The Fourier nodal coefficients are solved by the Newtonian algorithm in an incremental manner. When an approximate solution is available, an improved solution is obtained by solving a system of linear equations with the Fourier nodal increments as unknowns. The method is very suitable for parametric studies. When the excitation frequency is taken as a parameter, the free vibration response of various resonances can be obtained without actually computing the linear natural modes. For regular points along the response curves, the accuracy of the gradient matrix (Jacobian or tangential stiffness matrix) is secondary (cf. the modified Newtonian method). However, at the critical positions such as the turning points at resonances and the branching points at bifurcations, the gradient matrix becomes important. The minimum number of harmonic terms required is governed by the conditions of completeness and balanceability for predicting physically realistic response curves. The evaluations of the newly introduced mixed geometric matrices and their derivatives are given explicitly for the computation of the gradient matrix.     

A Semi-analytic Approach for Coherent Laser Radar System Efficiency

Abstract

A nearest Gaussian approximation (NGA) is proposed to approximate any shape for a single mode laser beam by a Gaussian shape. The application considered is a determination of the system efficiency in heterodyne coherent laser radar (HCLR). For an actual beam its NGA is defined by three parameters: the waist spot size and location, and an amplitude coefficient. These parameters are computed by a maximization of the norm of the scalar product written for the actual and Gaussian beams. In the case of the truncated Gaussian beam, particularly relevant to HCLR, the waist location can be analytically calculated, and only two parameters remain unknown: the waist spot size and amplitude coefficient. Using numerical applications, it is shown the NGA is in good agreement with Fresnel integral solution. The NGA combines a good accuracy and capability of analytical solutions. It can treat a variation in system efficiency owing to a misalignment angle between the transmitter and local oscillator.     

广义分数阶van der Pol-Duffing振子的动力学响应与隔振效果研究

用平均法研究了含分数阶导数项的van der Pol-Duffing振子的动力学行为和力传递率.得到了主共振时振子的一阶解析解、定常解幅频曲线和相频曲线的解析表达式,进一步通过与数值解作对比,验证了解析解的正确性,分析了不同参数对幅频曲线和力传递率的影响.结果 表明:解析解与数值解吻合良好;在无量纲情况下,共振区分数阶...     

非线性尘埃等离子体孤立子波变分迭代解

在全球气候变暖的极端反常的情形下，大气尘埃的扩散现象会带来巨大的灾害．本文研究了大气尘埃等离子体扩散的一类广义非线性孤立子波模型．首先对非扰动情形下利用待定系数法得到孤立子波解的解析表示式．其次用广义变分迭代的方法求出对应的变分乘子并构造变分迭代式，依次求出孤子波的各次迭代解．然后用行波变换得到广义非线性尘埃等离子体扰动模型的孤立子波的各次近似解．最后，由得到解的近似函数序列据变分理论知，在自变量的一定区域内此序列为一致收敛的．因此便证明了迭代解的极限函数是尘埃等离子体低频振动非线性方程的精确解．本文得到的近似解是尘埃等离子体的低频振动孤立子波的近似解析解，据它可用解析运算来求出相关量的物理性态，如孤立子波的波峰值．可以根据本文理论采取相应措施，避免出现电荷超高密度的聚集而导致放电击穿现象等．