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

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

THE DYNAMIC BUCKLING OF AN ORTHOTROPIC CYLINDRICAL THIN SHELL UNDER TORSION VARYING AS A POWER FUNCTION OF TIME

The subject of this investigation is to study the buckling of orthotropic cylindricalthin shells under torsion,which is a power function of time.The dynamic stability and compati-bility equations are obtained first.These equations are subsequently reduced to a time dependentdifferential equation with variable coefficient by using Galerkin's method.Finally,the critical dy-namic and static loading,the corresponding wave numbers,the dynamic factors,critical time andcritical impulse are found analytically by applying the Ritz type variational method.Using thoseresults,the effects of the variations of the power of time in the torsion load expression,of theloading parameter,the ratio of the Young's moduli and the ratio of the radius to thickness onthe critical parameters are studied numerically.It is observed that these factors have appreciableeffects on the critical parameters of the problem in the heading.     

Non-linear normal modes,invariance, and modal dynamics approximations of non-linear systems

Non-linear systems are here tackled in a manner directly inherited from linear ones, that is, by using proper normal modes of motion. These are defined in terms of invariant manifolds in the system's phase space, on which the uncoupled system dynamics can be studied. Two different methodologies which were previously developed to derive the non-linear normal modes of continuous systems — one based on a purely continuous approach, and one based on a discretized approach to which the theory developed for discrete systems can be applied-are simultaneously applied to the same study case-an Euler-Bernoulli beam constrained by a non-linear spring-and compared as regards accuracy and reliability. Numerical simulations of pure non-linear modal motions are performed using these approaches, and compared to simulations of equations obtained by a classical projection onto the linear modes. The invariance properties of the non-linear normal modes are demonstrated, and it is also found that, for a pure non-linear modal motion, the invariant manifold approach achieves the same accuracy as that obtained using several linear normal modes, but with significantly reduced computational cost. This is mainly due to the possibility of obtaining high-order accuracy in the dynamics by solving only one non-linear ordinary differential equation.     

Space-time perturbation modes for non-linear dynamic analysis of beams

In this work an attempt is made at bridging the powerful perturbation methods of analytical dynamics to the versatile finite element techniques which can readily handle arbitrarily complex structures. The proposed analysis methodology has two distinguishing features. First, a space-time finite element formulation is used, and hence the concept of modes is here naturally extended to that of space-time modes, where the time dependency is implied in the assumed modes. As a result, the partial differential equations of motion are directly reduced to purely algebraic non-linear simultaneous equations. Second, perturbation modes, rather than the usual vibration mode shapes are used and shown to be an appropriate basis for non-linear dynamic analysis. These modes bring information about the non-linearities of the system through the higher order derivatives of the strain and kinetic energies. The procedure is illustrated on non-linear beam problems and the results are compared with those of a full finite element model, i.e., when all the degrees of freedom are considered, as well as with analytical results, when available.     

复合结构剪切板的屈曲与后屈曲研究

基于Ｒｅｉｓｓｎｅｒ剪切板理论，本文将现有理论模型推进了一步，导出了控制复合材料和构造矩形板平衡的广义Ｋａｒｍａｎ大挠度方程。此外，运用摄动技术，构造了一个单向压缩简支剪切矩形板的后屈曲渐近解，适用于面板和夹心材料都具正交各向异性的夹层板以及计及横向剪切变形对称铺设的层合板。文中给出了典型算例，并将当前结果与其它理论或数值解作了比较，显示出相当好的一致性。     

充满水的金属圆柱壳经受轴向撞击时屈曲过程的实验研究

本文报道了五种规格充满水的圆柱壳经受大质量低速度轴向撞击作用时屈曲过程的实验结果。研究表明由于液体的存在，圆柱壳的屈曲模态呈轴对称型的环带波纹，整个撞击过程可以明显地分为振荡的动力加载、平稳发展的后屈曲和迅速下降的弹性卸载三个阶段。并讨论了试件壁厚和高度等几何参数对其屈曲性能的影响。     

水平圆管中受压扭作用管桩屈曲后的解析解

首先建立了考虑自重、受水平圆管约束的管柱在压担组合作用下的屈曲方程：此方程为四阶强非线性常微分方程，难以直接精确求解；通过管柱的正弦屈由分析，找到了方程中的小参数；利用小参数摄动法求解了强非线性屈曲方程，得到了管柱螺旋屈曲后的实际屈曲构形解析解，由此可进一步得到管柱发生螺旋屈曲的临界载荷；所得解析结果与数值计算结果及本文结果退化后与文献中的相关结果比较都有良好的一致性．     

A sharp condition for existence of an inertial manifold

It is shown that a perturbation argument that guarantees persistence of inertial (invariant and exponentially attracting) manifolds for linear perturbations of linear evolution equations applies also when the perturbation is nonlinear. This gives a simple but sharp condition for existence of inertial manifolds for semi-linear parabolic as well as for some nonlinear hyperbolic equations. Fourier transform of the explicitly given equation for the tracking solution together with the Plancherel's theorem for Banach valued functions are used.     

瞬态热传导问题的一阶对称SPH方法模拟

为提高传统光滑粒子动力学(smoothed particle hydrodynamics, SPH)方法模拟瞬态热传导问题的精度和稳定性,本文提出了一种一阶对称光滑粒子动力学(first order symmetric SPH, FO-SSPH)方法.该方法将具有二阶热传导方程分解成两个一阶偏微分方程,然后基于梯度离散和Taylor级数展开思想,对一阶核梯度形式进行修正,并将得到的局部矩阵对称化.数值结果表明:与传统SPH方法相比,FO-SSPH方法精度高、数值稳定性好; 该方法能较准确地直接施加混合边值 关键词： 瞬态热传导 光滑粒子动力学 非线性     

Axial instability of double-nanobeam-systems

This Letter considers the axial instability of double-nanobeam-systems. Eringen's nonlocal elasticity is utilized for modelling the double-nanobeam-systems. The nonlocal theory accounts for the small-scale effects arising at the nanoscale. The small-scale effects substantially influence the instability (or buckling) of double-nanobeam-systems. Results reveal that the small-scale effects are higher with increasing values of nonlocal parameter for the case of in-phase (synchronous) buckling modes than the out-of-phase (asynchronous) buckling modes. The increase of the stiffness of the coupling elastic medium in double-nanobeam-system reduces the small-scale effects during the out-of-phase (asynchronous) buckling modes. Analysis of the scale effects in higher buckling loads of double-nanobeam-system with synchronous and asynchronous modes is also discussed in this Letter. The theoretical development presented herein may serve as a reference for nonlocal theories as applied to the instability analysis of complex-nanobeam-system such as complex carbon nanotube system.