非线性动力系统的频率近似法

给出非线性动力系统周期振动的频率近似法，本法将描述动力系统的非线性微分方程，化为以相角为自变量，振动频率为未知函数的积分方程，将弹性恢复力表示为线性及非线性两部分，从而得到积分方程的近似解，即频率的近似表达式。     

无网格局部Petrov-Galerkin方法在弹塑性断裂力学问题中的应用

采用无网格局部Petroy-Galerkin方法来分析弹塑性断裂力学问题.这种无网格方法采用移动最小二乘法(MLS)来构造近似试函数和采用Heaviside函数作为加权残值法中的权函数,由于近似函数不满足KroneckerDelta条件,因此采用直接插值法来施加本质边界条件.如果不考虑体力,所形成的整体刚度矩阵只包含局部边界积分,而不包含局部域积分和奇异积分.采用增量Newton-Raphson迭代法来求解弹塑性增量形式的局部Petrov-Galerkin方程.数值算例结果表明,该文方法对于弹塑性断裂力学问题的求解是可行的和有效的,并且所得到的结果具有较好的精度.     

非齐次动力方程Duhamel项的精细积分

提出了不需要矩阵求逆运算的求解Duhamel积分项的精细积分方法.通过将精细积分法的关键思想--加法定理和增量存储--直接应用于Duhamel积分响应矩阵的求解,可给出当非齐次项分别为多项式、正弦/余弦以及指数函数等基本形式时Duhamel积分在计算机上的精确解.特别的,该算法不依赖于系统矩阵(或相关矩阵)的形态.当系统矩阵奇异或接近奇异时,其优越性更为显著.算例验证了该算法的有效性.     

多层地基条带基础动力刚度矩阵的精细积分算法

提出应用精细积分算法计算多层地基的动力刚度问题. 精细积分是计算层状介质中波传播的高效而精确的数值方法. 利用傅里叶积分变换将层状地基的波动方程转换为频率-波数域内的两点边值问题的常微分方程组, 运用精细积分方法求解格林函数, 最后再将得到的频率-波数域内地基表面的动力刚度矩阵转换到频率-空间域内, 进而得到刚性条带基础频率域的动力柔度或刚度矩阵. 所建议的精细积分算法, 可以避免一般传递矩阵计算中的指数溢出问题, 对各种情况有广泛的适应性, 计算稳定, 在高频段可以保障收敛性, 并能达到较高的计算精度.      

基于对偶变量变分原理的完整约束系统保辛算法

基于对偶变量变分原理,选择积分区间两端位移为独立变量,构造了求解完整约束哈密顿动力系统的高阶保辛算法。首先,利用拉格朗日多项式对作用量中的位移、动量及拉格朗日乘子进行近似;然后,对作用量中不包含约束的积分项采用Gauss积分近似,对作用量中包含约束的积分项采用Lobatto积分近似,从而得到近似作用量;最后,在此近似作用量的基础上,利用对偶变量变分原理,将求解完整约束哈密顿动力系统问题转化为一组非线性方程组的求解。算法具有保辛性和高阶收敛性,能够在位移的插值点处高精度地满足完整约束。算法的收敛阶数及数值性质通过数值算例验证。     

考虑频率依赖性的涂层复合结构固有特性求解

粘弹性阻尼材料的力学特性参数会随着频率的变化而改变,即具有频率依赖性,因而传统的动力学建模及分析方法不能满足实际涂层结构优化设计的需要。在简要介绍粘弹性阻尼材料频率依赖性的基础上,本文提出用特征向量增值法来求解涂层复合结构的固有特性,并详细推导了特征向量增值法的求解原理。由此,提出了特征向量增值法的计算流程,包括计算无阻尼系统的固有特性,用Fox and Kapoor或者Nelson方法计算复特征向量增量;用Rayleigh熵法求解复特征值。最后,以涂敷粘弹性阻尼材料的钛基薄板为例,求解了该复合结构的固有特性,并与经典的模态应变能法进行了比较,证明了所提方法的正确性。     

A PRECISE INTEGRATION APPROACH FOR THE DYNAMIC-STIFFNESS MATRIX OF STRIP FOOTINGS ON A LAYERED MEDIUM 1）

提出应用精细积分算法计算多层地基的动力刚度问题.精细积分是计算层状介质中波传播的高效而精确的数值方法.利用傅里叶积分变换将层状地基的波动方程转换为频率-波数域内的两点边值问题的常微分方程组,运用精细积分方法求解格林函数,最后再将得到的频率-波数域内地基表面的动力刚度矩阵转换到频率-空间域内,进而得到刚性条带基础频率域的动力柔度或刚度矩阵.所建议的精细积分算法,可以避免一般传递矩阵计算中的指数溢出问题,对各种情况有广泛的适应性,计算稳定,在高频段可以保障收敛性,并能达到较高的计算精度.     

外部激励Stuart-Landau模型的数值研究

将非线性常微分方程组周期解的求解看作一个边值问题 ,运用Newton迭代构造求解这组方程的数值方法。利用上述方法求得了激励Stuart Landau方程的周期解 ,研究了圆柱振动对圆柱后Karman涡街的抑制现象 ,和振动的频率锁定现象 ,证明了激励Stuart Landau方程描写钝体尾迹动力系统的有效性     

VISAR测速中的信号丢失及丢失条纹数的确定

论述了VISAR测速中信号频率与被测速度增量的关系和光电倍增管、数字示波器所能响应的最高速度增量.分析了信号丢失的原因,给出了丢失条纹数的确定方法.最后对VISAR应用中如何正确选择条纹常数提出了建议.     

同震断层三维裂纹扩展波动时域超奇异积分法

应用波动时域超奇异积分法将P波、S波和磁电热弹多场耦合作用下同震断层任意形状三维裂纹扩展问题转化为求解以广义位移间断率为未知函数的超奇异积分方程组问题;定义了广义应力强度因子,得到裂纹前沿广义奇异应力增量解析表达式;应用波动时域有限部积分概念及体积力法,为超奇异积分方程组建立了数值求解方法,编制了FORTRAN程序,以三维矩形裂纹扩展问题为例,通过典型算例,研究了广义应力强度因子随裂纹位置变化规律;分析了同震断层裂纹扩展中力、磁、电场辐射规律.      

A new numerical technique for the analysis of parametrically excited nonlinear systems

A new computational scheme using Chebyshev polynomials is proposed for the numerical solution of parametrically excited nonlinear systems. The state vector and the periodic coefficients are expanded in Chebyshev polynomials and an integral equation suitable for a Picard-type iteration is formulated. A Chebyshev collocation is applied to the integral with the nonlinearities reducing the problem to the solution of a set of linear algebraic equations in each iteration. The method is equally applicable for nonlinear systems which are represented in state-space form or by a set of second-order differential equations. The proposed technique is found to duplicate the periodic, multi-periodic and chaotic solutions of a parametrically excited system obtained previously using the conventional numerical integration schemes with comparable CPU times. The technique does not require the inversion of the mass matrix in the case of multi degree-of-freedom systems. The present method is also shown to offer significant computational conveniences over the conventional numerical integration routines when used in a scheme for the direct determination of periodic solutions. Of course, the technique is also applicable to non-parametrically excited nonlinear systems as well.     

1/3 subharmonic solution of elliptical sandwich plates

IntroductionInrecentyears,withtheessentialadvantageoflightweightandhighrigidity ,sandwichplatesandshellshavebeenusedasanimportantpatternofstructuralelementsinaeronautical,astronauticalandnavalengineering .However,nonlinearproblemsforsandwichplatesandshellsareonlyinvestigatedbyafewbecauseofthedifficultiesofnonlinearmathematicalproblems.LiuRen_huaiandXuJia_chu[1,2 ]andothershavemadesomeinvestigationsinthisfield .Bifurcationofnonlinearvibrationforsandwichplateshasnotyetbeeninvestigated .Inthisp…     

Nonlinear vibration of a two-mass system with nonlinear stiffnesses

An analytical approach is developed for nonlinear free vibration of a conservative, two-degree-of-freedom mass–spring system having linear and nonlinear stiffnesses. The main contribution of the proposed approach is twofold. First, it introduces the transformation of two nonlinear differential equations of a two-mass system using suitable intermediate variables into a single nonlinear differential equation and, more significantly, the treatment a nonlinear differential system by linearization coupled with Newton’s method and harmonic balance method. New and accurate higher-order analytical approximate solutions for the nonlinear system are established. After solving the nonlinear differential equation, the displacement of two-mass system can be obtained directly from the governing linear second-order differential equation. Unlike the common perturbation method, this higher-order Newton–harmonic balance (NHB) method is valid for weak as well as strong nonlinear oscillation systems. On the other hand, the new approach yields simple approximate analytical expressions valid for small as well as large amplitudes of oscillation unlike the classical harmonic balance method which results in complicated algebraic equations requiring further numerical analysis. In short, this new approach yields extended scope of applicability, simplicity, flexibility in application, and avoidance of complicated numerical integration as compared to the previous approaches such as the perturbation and the classical harmonic balance methods. Two examples of nonlinear two-degree-of-freedom mass–spring system are analyzed and verified with published result, exact solutions and numerical integration data.     

Analytical approximation to large-amplitude oscillation of a non-linear conservative system

This paper deals with non-linear oscillation of a conservative system having inertia and static non-linearities. By combining the linearization of the governing equation with the method of harmonic balance, we establish analytical approximate solutions for the non-linear oscillations of the system. Unlike the classical harmonic balance method, linearization is performed prior to proceeding with harmonic balancing, thus resulting in a set of linear algebraic equations instead of one of non-linear algebraic equations. Hence, we are able to establish analytical approximate formulas for the exact frequency and periodic solution. These analytical approximate formulas show excellent agreement with the exact solutions, and are valid for small as well as large amplitudes of oscillation.     

Accurate Higher-Order Analytical Approximate Solutions to Large-Amplitude Oscillating Systems with a General Non-Rational Restoring Force

A new approximate analytical approach for accurate higher-order nonlinear solutions of oscillations with large amplitude is presented in this paper. The oscillatory system is subjected to a non-rational restoring force. This approach is built upon linearization of the governing dynamic equation associated with the method of harmonic balance. Unlike the classical harmonic balance method, simple linear algebraic equations instead of nonlinear algebraic equations are obtained upon linearization prior to harmonic balancing. This approach also explores large parameter regions beyond the classical perturbation methods which in principle are confined to problems with small parameters. It has significant contribution as there exist many nonlinear problems without small parameters. Through some examples in this paper, we establish the general approximate analytical formulas for the exact period and periodic solution which are valid for small as well as large amplitudes of oscillation.     

含外激励van der Pol-Mathieu方程的非线性动力学特性分析

本文针对含有自激励, 参数激励和外激励等三种激励联合作用下van der Pol-Mathieu方程的周期响应和准周期运动进行分析, 发现其准周期运动的频谱中含有均匀边频带这一新的特性. 首先, 采用传统的增量谐波平衡法(IHB法)分析了van der Pol-Mathieu方程的周期响应, 得到了其非线性频率响应曲线; 再利用Floquet理论对周期解进行稳定性分析, 得到了两种类型的分岔及它们的位置. 然后, 基于van der Pol-Mathieu方程准周期运动的频谱中边频带相邻频率之间是等距的且含有两个不可约的基频的特性(其中一个基频是已知的, 另一个基频事先是未知的), 推导了相应的两时间尺度IHB法, 精确计算出van der Pol-Mathieu方程的准周期运动的另一个未知基频和所有的频率成份及其对应的幅值, 尤其在临界点附近处的准周期运动响应. 得到的准周期运动结果和利用四阶龙格-库塔(RK)数值法得到的结果高度吻合. 最后, 研究发现了含外激励van der Pol-Mathieu方程在不同激励频率时的一些丰富而有趣的非线性动力学现象.      

Construction of approximate analytical solutions to strongly nonlinear damped oscillators

An analytical approximate method for strongly nonlinear damped oscillators is proposed. By introducing phase and amplitude of oscillation as well as a bookkeeping parameter, we rewrite the governing equation into a partial differential equation with solution being a periodic function of the phase. Based on combination of the Newton’s method with the harmonic balance method, the partial differential equation is transformed into a set of linear ordinary differential equations in terms of harmonic coefficients, which can further be converted into systems of linear algebraic equations by using the bookkeeping parameter expansion. Only a few iterations can provide very accurate approximate analytical solutions even if the nonlinearity and damping are significant. The method can be applied to general oscillators with odd nonlinearities as well as even ones even without linear restoring force. Three examples are presented to illustrate the usefulness and effectiveness of the proposed method.     

New explicit and exact traveling wave solutions for two nonlinear evolution equations

In this paper, with the aid of computer symbolic computation system such as Maple, an algebraic method is firstly applied to two nonlinear evolution equations, namely, nonlinear Schrodinger equation and Pochhammer–Chree (PC) equation. As a consequence, some new types of exact traveling wave solutions are obtained, which include bell and kink profile solitary wave solutions, triangular periodic wave solutions, and singular solutions. The method is straightforward and concise, and it can also be applied to other nonlinear evolution equations in mathematical physics.     

求解非线性振动问题的一种新方法

首先把描述非线性振动的微分方程归结为一非线性积分微分方程，然后把此积分微分方程的求解转化为一无穷阶的非线性代数方程组的求解。从理论上讲，可得到满足任何精度要求的周期解。本文用此方法对Ｄｕｆｉｎｇ系统进行了分析。     

An enriched multiple scales method for harmonically forced nonlinear systems

This article explores enrichment to the method of Multiple Scales, in some cases extending its applicability to periodic solutions of harmonically forced, strongly nonlinear systems. The enrichment follows from an introduced homotopy parameter in the system governing equation, which transitions it from linear to nonlinear behavior as the value varies from zero to one. This same parameter serves as a perturbation quantity in both the asymptotic expansion and the multiple time scales assumed solution form. Two prototypical nonlinear systems are explored. The first considered is a classical forced Duffing oscillator for which periodic solutions near primary resonance are analyzed, and their stability is assessed, as the strengths of the cubic term, the forcing, and a system scaling factor are increased. The second is a classical forced van der Pol oscillator for which quasiperiodic and subharmonic solutions are analyzed. For both systems, comparisons are made between solutions generated using (a) the enriched Multiple Scales approach, (b) the conventional Multiple Scales approach, and (c) numerical simulations. For the Duffing system, important qualitative and quantitative differences are noted between solutions predicted by the enriched and conventional Multiple Scales. For the van der Pol system, increased solution flexibility is noted with the enriched Multiple Scales approach, including the ability to seek subharmonic (and superharmonic) solutions not necessarily close to the linear natural frequency. In both nonlinear systems, comparisons to numerical simulations show strong agreement with results from the enriched technique, and for the Duffing case in particular, even when the system is strongly nonlinear.