非线性振动系统的多项式向量方法

对于常微分方程描述的非线性振动系统,当采用摄动方法求近似解时,先是给出满足各阶近似解的二阶常微分方程组,继而依次对每一个常微分方程进行求解,以致多自由度非线性振动系统的求解过程相当繁琐.文章针对常微分方程表示的非线性振动系统,提出了一种求解非线性振动系统近似解的多项式向量方法,该方法将二阶常微分方程组表示成一阶状态方程组,将非线性部分写成常数矩阵和多项式向量之积的形式.然后,采用直接摄动方法,获得每个幂次近似解所满足的一组状态方程,此时状态方程的非线性部分成为常数矩阵和前一幂次近似解作为元素组成的多项式向量的乘积.进一步,借助Toeplitz矩阵将多项式向量之乘法表示成矩阵形式,以解决多项式相乘带来的幂次方系数的确定问题,再根据一阶非齐次方程组的求解方法,获得状态方程组的全部近似解析解.多项式向量方法将二阶常微分描述的非线性振动求解过程转换为一阶非齐次状态方程组的求解问题,计算过程主要是矩阵和向量之间乘法运算,提高了计算效率和程序化水平.     

单桩竖向动力阻抗计算模型研究

运用土动力学和结构动力学原理,同时考虑桩周土的弱化效应和桩-土界面的相对滑移效应,利用数理方程方法分别求解单桩与桩周近场土域及远场土域的振动方程,建立了竖向荷载作用下单桩动力阻抗函数的简化计算方法,提出了一种改进的非线性动力Winkler模型,确定了模型中各物理元件的参数,进而通过对比分析验证了计算模型的合理性,从而为桩-土-上部结构耦合系统的非线性分析奠定了基础。     

表面堆载作用下群桩负摩擦研究

利用Biot固结理论和积分方程方法研究了表面有堆载的群桩负摩擦问题。根据基本解得出了群桩在圆形均布载荷作用下在时间域内的第二类Fredholm积分方程组。运用Laplace变换对上述积分方程组进行简化，求解上述积分方程组并进行相应的数值逆变换就可得出群桩在表面圆形均布载荷作用下的变形、轴力、孔压和桩侧摩阻力随时间的变化情况。     

基于叠加原理的群管桩的纵向振动

利用Winkler地基模型来模拟管桩与桩周土、管桩与桩芯土的动力相互作用,在忽略桩周土、桩芯土的径向、环向位移的情况下,得到了Winkler地基模型的刚度系数和阻尼系数。在Winkler地基模型的基础上建立了主动管桩和被动管桩的纵向振动方程,将相互作用因子的概念扩展应用到管桩-管桩纵向动力相互作用,得到了管桩-管桩纵向动力相互作用因子。在叠加原理和管桩-管桩动力相互作用因子的基础上,求解了管桩的纵向振动,得到了管桩的纵向动力阻抗。通过数值算例分析了相关参量对管桩-管桩纵向动力相互作用和群管桩纵向振动的影响。研究结果表明:管桩内外半径比、桩土剪切模量比、管桩长径比对管桩-管桩纵向动力相互作用和群管桩纵向振动有较大的影响,桩间距较大时管桩-管桩纵向动力相互作用因子和群桩动力阻抗随频率变化曲线波动越大;管桩壁厚不宜过薄。     

垂直受荷群桩在半空间饱和土中的稳态反应

研究了均质半空间饱和土中群桩在垂直稳态荷载作用下的动力反应问题。半空间饱和土采用Biot提出的三维波动原理。将桩看作是一维的弹性杆单元。利用圆形简谐载荷作用下的Biot固结方程的基本解和桩土之间的变形协调条件得到饱和土中群桩的第二类Fredholm积分方程,群桩的动力阻抗采用动力相互作用因子的方法。分析桩距等参数对群桩阻抗、轴力和孔压等的影响。本文的方法可以用于计算层状饱和土中群桩的动力反应问题。     

不可压缩平面问题的位移-压力混合重心插值配点法

引入人工压力变量,将弹性本构方程以应力、应变和压力表达,建立求解不可压缩平面弹性问题的位移-压力方程和不可压缩条件方程的耦合偏微分方程组。利用张量积型重心Lagrange插值近似二元函数,得到计算插值节点处偏导数的偏微分矩阵。采用配点法离散不可压缩弹性控制方程,利用偏微分矩阵直接离散弹性力学控制方程为矩阵形式方程组。利用插值公式离散位移和应力边界条件,将离散边界条件与离散控制方程组合为新的方程组,得到求解弹性问题的过约束线性代数方程组;利用最小二乘法求解线性方程组,得到弹性力学问题位移数值解。数值算例验证了所提方法的数值计算精度为10-14~10-10。     

电场作用下石英晶体板非线性高频振动分析

建立了强电场作用下考虑材料和几何非线性的石英晶体板厚度剪切振动和弯曲振动的非线性方程组。基于厚度剪切振动近似理论的线性位移解,利用伽辽金法将非线性偏微分方程组转化为关于时间变量的常微分方程组。接着利用逐次近似法获得了石英晶体板厚度剪切振动和弯曲振动耦合的频率响应关系,并绘制了不同电压下厚度剪切振动的频率响应曲线图。数值计算结果表明,电压对石英晶体板厚度剪切振动的频率有着显著影响,其引起的频率漂移值超过了常见压电谐振器的允许值。所建立的方程可以用于石英晶体谐振器的非线性有限元分析和偏场效应求解。     

非线性动力学常微分方程组高精度数值积分方法

建立了一种求解非线性动力学常微分方程组初值问题的新方法．若非线性函数一阶导数存在，则给出解的积分方程表达式，计算得到按规定误差要求的高精度数值解．引入一般自治或非自治非线性系统的首次近似Jacobi矩阵，不作任何假设重构等价的非线性常微分方程组，简捷而有广泛的适应性，不改变方程的本质，但其主项构成线性化方程组，其它项则代表非线性函数高阶余项而不涉及Taylor级数展开计算，给出该方程组初值问题的Duhamel卷积分解析表达式，在时间步长内进行数值积分选代求解，在指定误差内快速收敛，逐步递推获得非线性常微分方程的瞬态响应和全时域高精度数值解．积分解连续满足微分方程组而不是在离散的步长端点上满足代数方程组，打破了传统用增量法在离散点上建立的代数方程组迭代求解，从而使传统Euler型逐步积分法的各种差分格式算法改变成真正的积分格式算法．数值计算中给出指数矩阵递增展开式，变矩阵乘法为乘积系数的加法，避免了大量矩阵自乘而大大提高计算效率．算法验证为无条件稳定，则保证对线性常微分方程而言，计算中舍入误差的传播不会扩散，不出现计算机字长有限而引起舍入误差导致计算不确定性问题．基于以上理论和数值方法，计算了线性非线性算例并进行了分析，验证了本方法简捷而有广泛的适应性，可以有足够的精确性．     

抗滑桩与滑坡相互作用分析的常微分方程组边值法

将抗滑桩与滑坡的相互作用抽象为一阶常微分方程组在特定边界条件约束下的定解问题,通过龙格-库塔(Runge-Kutta)差分法求解该方程组可得桩身内力与变形及滑床抗力,该方法不同于常规的基于桩体挠曲四阶微分方程的级数解法、差分法及有限元法,可以一次性解出桩身内力及变形值,无需进行二次换算;可对全桩进行整体分析,无需像传统方法那样以滑动面为界将桩身分成受荷段和锚固段分别计算;该方法一改传统的求解高阶微分方程为解低阶微分方程组,符合计算力学的优化思想,并且可以方便地考虑桩身的剪切变形,为抗滑桩的计算分析提供了一种切实可行的新思路,可作为传统抗滑桩内力分析方法的有效补充。本文还编写了基于该法的全桩内力计算和图形处理程序。工程算例表明,该方法与传统方法的计算结果能很好吻合,且程序运行效率更高。     

层状土中考虑频域内轴向力分担的群桩水平动力阻抗

以动力文克尔梁模型为基本理论，在充分考虑了土的分层以及桩顶轴向力参与作用后导出了一种新的计算层状介质中水平动力相互作用因子的方法，进一步扩展了群桩-土相互作用简化法的应用范围，可以计算任意构型、桩数、土层属性下的群桩系统阻抗反应，其计算工作量较小，是其它数值方法难以比拟的。数值计算结果表明；桩顶轴向力的引入不能从本质上改变桩基振动形式，但考虑竖向振动附加影响的水平激振作用会产生相当的动力P-△效应，致使水平单、群桩阻抗效应降低，位移反应增大，同时桩间动力相互作用影响也较单-水平力作用时有明显规律性的增强。     

Accurate analytical perturbation approach for large amplitude vibration of functionally graded beams

The present work derives the accurate analytical solutions for large amplitude vibration of thin functionally graded beams. In accordance with the Euler–Bernoulli beam theory and the von Kármán type geometric non-linearity, the second-order ordinary differential equation having odd and even non-linearities can be formulated through Hamilton's principle and Galerkin's procedure. This ordinary differential equation governs the non-linear vibration of functionally graded beams with different boundary constraints. Building on the original non-linear equation, two new non-linear equations with odd non-linearity are to be constructed. Employing a generalised Senator–Bapat perturbation technique as an ingenious tool, two newly formulated non-linear equations can be solved analytically. By selecting the appropriate piecewise approximate solutions from such two new non-linear equations, the analytical approximate solutions of the original non-linear problem are established. The present solutions are directly compared to the exact solutions and the available results in the open literature. Besides, some examples are selected to confirm the accuracy and correctness of the current approach. The effects of boundary conditions and vibration amplitudes on the non-linear frequencies are also discussed.     

Response of a non-linear system with strong damping to multifrequency excitations

Summary In this paper, a new asymptotic calculation method is presented for a class of nonlinear nonautonomous vibration systems with multiple external periodic interferences. Simple calculation formulae of an asymptotic solution for resonance and off-resonance vibration are derived. The paper is concerned with a vibration system representing a class of nonlinear oscillators. Consequently, the calculation of a class of nonlinear oscillators is routinized. Two different Duffing equations are verified which shows that the results are completely in accordance with the solutions of references [3, 4, 6]. The derivation of solutions of Duffing equations becomes easier and simpler. In addition, some errors in reference [6] are pointed out.     

Vibration and stability of an aircraft tail under simultaneous primary-combined and internal resonance

This paper adds a negative velocity feedback to the dynamical system of twin-tail aircraft to suppress the vibration. The system is represented by two coupled second-order nonlinear differential equations having both quadratic and cubic nonlinearities. The system describes the vibration of an aircraft tail subjected to both multi-harmonic and multi-tuned excitations. The method of multiple time scale perturbation is adopted to solve the nonlinear differential equations and obtain approximate solutions up to the third order approximations. The stability of the proposed analytic solution near the simultaneous primary, combined and internal resonance is studied and its conditions are determined. The effect of different parameters on the steady state response of the vibrating system is studied and discussed by using frequency response equations. Some different resonance cases are investigated numerically     

A new matrix perturbation method for analytical solution of the complex modal eigenvalue problem of viscously damped linear vibration systems

A new matrix perturbation analysis method is presented for efficient approximatesolution of the complex modal quadratic generalized eigenvalue problem of viscouslydamped linear vibration systems.First,the damping matrix is decomposed into the sum of aproportional-and a nonproportional-damping parts,and the solutions of the real modaleigenproblem with the proportional dampings are determined,which are a set of initialapproximate solutions of the complex modal eigenproblem.Second,by taking thenonproportional-damping part as a small modification to the proportional one and using thematrix perturbation analysis method,a set of approximate solutions of the complex modaleigenvalue problem can be obtained analytically.The result is quite simple.The new methodis applicable to the systems with viscous dampings-which do not deviate far away from theproportional-damping case.It is particularly important that the solution technique be alsoeffective to the systems with heavy,but not over,dampings.The solution formul     

用薄层法分析层状地基中各种基础的阻抗函数

薄层法是分析和模拟弹性波在层状介质中传播的一种半解析半数值方法.采用薄层单元和傍轴边界分别模拟层状地基和弹性半空间.利用薄层法位移基本解和容积法推导了层状地基中基础-地基动力相互作用方程及块式基础、桩基础和承台群桩基础阻抗函数的统一计算公式.通过计算半无限弹性地基中桩基础、块式基础和二层地基上基础的阻抗函数验证了方法的适用性.进而计算了某实际层状地基中承台群桩基础的阻抗函数,并与试验结果进行对比,两者吻合较好.本文方法可用于分析弹性层状地基中各种基础的阻抗函数.     

梁-桩-土竖向耦合振动特性分析

针对桩与连续梁组合这一基础形式,研究了梁—桩—土竖向耦合振动特性. 首先,将桩假设为一维黏弹性杆件,采用平面应变模型模拟桩侧成层土对桩的动力作用. 同时,桩顶对梁的支撑简化为竖向点支撑. 然后,分别根据一维杆件纵向振动理论和Timoshenko梁理论给出桩和梁的竖向振动控制方程并求解,进而借助离散傅里叶逆变换获得在梁上瞬时激振下梁和土层以上桩段的时域响应半解析解. 通过与有限元模拟结果对比,验证了解的合理性. 在此基础上,讨论了梁几何参数和桩身缺陷对梁—桩—土竖向耦合振动特性的影响.      

Recursive geometric integrators for wave propagation in a functionally graded multilayered elastic medium

The differential equations governing transfer and stiffness matrices and acoustic impedance for a functionally graded generally anisotropic magneto-electro-elastic medium have been obtained. It is shown that the transfer matrix satisfies a linear 1st order matrix differential equation, while the stiffness matrix satisfies a nonlinear Riccati equation. For a thin nonhomogeneous layer, approximate solutions with different levels of accuracy have been formulated in the form of a transfer matrix using a geometrical integration in the form of a Magnus expansion. This integration method preserves qualitative features of the exact solution of the differential equation, in particular energy conservation. The wave propagation solution for a thick layer or a multilayered structure of inhomogeneous layers is obtained recursively from the thin layer solutions. Since the transfer matrix solution becomes computationally unstable with increase of frequency or layer thickness, we reformulate the solution in the form of a stable stiffness-matrix solution which is obtained from the relation of the stiffness matrices to the transfer matrices. Using an efficient recursive algorithm, the stiffness matrices of the thin nonhomogeneous layer are combined to obtain the total stiffness matrix for an arbitrary functionally graded multilayered system. It is shown that the round-off error for the stiffness-matrix recursive algorithm is higher than that for the transfer matrices. To optimize the recursive procedure, a computationally stable hybrid method is proposed which first starts the recursive computation with the transfer matrices and then, as the thickness increases, transits to the stiffness matrix recursive algorithm. Numerical results show this solution to be stable and efficient. As an application example, we calculate the surface wave velocity dispersion for a functionally graded coating on a semispace.     

内共振条件下轴向运动梁磁弹性超谐共振研究

以载流导线激发的磁场中轴向运动梁为研究对象,同时考虑外激励力作用,推导出梁的磁弹性非线性振动方程．通过位移函数的设定和伽辽金积分法,将非线性振动方程离散为常微分方程组．采用多尺度法得到系统的近似解析解．应用Matlab 和Mathematica 软件求解幅频响应方程,并对稳态解进行稳定性判定．通过具体算例得到前两阶假设模态的响应幅值随不同参数的变化规律．结果发现：系统在内共振条件下发生超谐波共振时,二阶假设模态幅值明显小于一阶;随着外激励的增大,多值解区域范围明显缩小;随着电流强度增加,振动幅值减小,表明载流导线能够起到控制共振的作用．     

Nonlinear resonance of a subsatellite on a short constant tether

The nonlinear resonant behavior of a subsatellite on a short constant tether during station-keeping phase is investigated in this paper. The nonlinear dynamic equations of in-plane motion of the system are derived based on Kane’s method first. Then an approach of multiple scales expressed in matrix form is employed in solving the simplified nonlinear system of cubic nonlinearity near its local equilibrium position. Analysis shows that there exists a three-to-one resonance in such a nonlinear system with two degrees of freedom. Afterward, the approximate solution up to third order determined analytically by the Weierstrass elliptic function is obtained and the comparison between the approximate and numerical solutions presented as well. The results show that the approximate solution is coincide well with the numerical solution of original system. The nonlinear resonance of the subsatellite on short tether exhibits coexistent quasiperiodic motions or a quasiperiodic oscillation near local equilibrium position.