非线性动力学方程的一种级数解

本文针对n维未知向量v的一阶微分议程v=Hv f(v,t)进行求解,其中Hv和f(v,t)分别是右端项的线性剂次部分和非线性部分,首先,将非线性部分f(v,t)在所论时刻tk处展成t-tk=τ的泰勒级数,并通过求原函数的方法直接给出每一项的积分,从而获得了待求微分方程在级数形式下的闭合争,它的具有不同精度的各次静似解可表示成τ的分段解析函数,便于研究非线性动力学行为与其物理参数的依赖关系。本文还用算例验证了各次近似解之间的数值关系,并和解析解等作了比较。     

旋转机械系统多自由度非线性动力学数值分析

首先将转子系统的动力响应问题归结为2n维未知向量v的一阶非线性动力学方程dv／dt=Hv＋f（v，t），并给出了求解这一方程的一次近似式法和三次多项式迫近法。在非稳态、非线性油膜力等作用下，以刚性Jeffcott转子与112个自由度的汽轮发电机组低压转子系统为例，用上述求解方法分析了它们的动力响应及非线性动力学特性；其问，还将计算结果与Runge—Kutta法、Newmark法的相应结果进行了比较，并深入讨论了数值稳定性问题。汽轮发电机组的算例表明对一些具有较复杂的非线性右端项，、同时规模又较大的问题，如果采用四阶Runge—Kutta法，才算几步就因数值骤然增大而失控；但若用同样步长的一次近似式，由于它是一种显式的无条件稳定算式，则计算过程迅速且结果合理可靠。     

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

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

强非线性振动系统周期解的能量迭代法

对于完全强非线性系统：x^.. g(x) f(x,x^.)x^.=0,提出求周期近似解析解以及这些解的稳定性的新方法。式中，g(x)、f(x,x^.)x^.分别是x,x、x^.的非线性函数。方法是基于能量原理，求出其一次近似解析解，然后引进牛顿迭代思想，得到周期系统数微分方程，最后根据谐波平衡原理及最小二乘法求其高次近似解，高次近似解的表达式由计算机辅助推导。计算参考文献[2]和[3]中的例题，令其中ε=1，研究该完全强非线性系统的周期解及其稳定性，本文方法与龙格-库塔数值法算得的结果对照如图1-3所示，它们表明本文方法不仅有效而且精度较高。     

对一般Hamilton体系近似解保辛条件的讨论

文献[1]给出的一般Hamilton体系近似解保辛的条件,尚需讨论。其中(3.5~6)的证明适用于v·=H(z)v的系统,要求H(z)是Hamilton矩阵,齐次方程。即使Hamilton矩阵是与位移有关的H(z,q),仍可适用。但一般Hamilton体系未必能表示为v·=H·v,例如存在有势外力的情况,此时是非齐次线性方程。故在一般情况下,近似解是否保辛的原则尚需明确。一般的Hamilton正则方程体系非线性,通常用数值积分近似求解,故时间坐标离散而成为长η的序列,离散坐标体系。分析结构力学[2]考虑了离散坐标的情况,其中证明了区段两端状态向量之间关系ζ=ζ(v)的微商S(v)是…     

非结构网格上p型多重网格法流场数值模拟

在二维、三维非结构网榕上,针对间断Galerkin方法计算量大、收敛慢的缺点将p型多重网格方法应用于该方法求解跨音速Euler方程,提高计算效率。p型多重网格方法是通过对不同阶次多项式近似解进行递归迭代求解,来达到加速收敛。文中对高阶近似(p＞0)使用显式格式,最低阶近似(p=0)采用隐式格式。NACA0012翼型和O...     

被动隔振体非线性振动的能量迭代解法

研究了由基础振动激励、弹性材料隔离的被动隔振体的强非线性动力响应。用变形的三次多项式函数表征隔振材料的非线性刚度特性,建立了被动隔振体的非线性动力学方程,得到有阻尼受迫振动Duffing方程。将求解强非线性自治系统的能量迭代方法加以改进,推广应用到强非线性非自治系统,求出周期响应的近似解析解表达式,以及幅频关系、相频关系和隔振系数的近似表达式。算例中应用本方法与Runge-Kutta方法进行了对照,结果表明求解精度较高。本文利用计算机进行了辅助推导。     

摄动法在研究弹-粘塑性波问题中的应用

采用控制金属材料宏观塑性流动的两个无量纲物理参数作为小参数,将一维弹/粘-塑性问题的解摄动展开,从而,求解非线性波动方程的问题可以转化成求解相应的齐次或非齐次电报方程的问题,用Laplace积分变换或级数展开技术首先得到零次精确解。然后,用Riemann函数方法可获得一次和高次摄动解。与非线性问题的数值解比较,在恒应力或恒速度边界条件下,一次摄动解给出了波动问题的良好近似。这就表明,摄动技术在研究一类广泛的弹/粘-塑波问题中是有效的。     

等效线性化——模态综合迭代法

本文提出一个求解具有非线性元件的转子系统次谐波响应的方法.分析和实例计算结果表明:此方法可得到近似度较高的近似解,收敛速度快,而接通常地在模态综合法中,仅考虑代数对接条件,得到的近似解近似度较低.     

一类平面Hamilton系统的大范围非线性化近似方法

在非线性动力系统的研究已经进入了占主导地位的时期,对其提出大范围的非线性化近似方法具有特别重要的意义.在本文中,我们主要对于一类典型的Hamilton系统,根据等势线有两个,或者三个交点的不同情形,给出7种不同的大范围最低次非线性化近似系统,并通过积分近似系统给出近似解(轨道).结果表明,近似椭圆周期轨道可通过线性化近似系统得到,而同(异)宿轨道则可通过2、3次非线性化近似系统得到.最后,将近似方法应用于一个具体Hamilton系统的分析.     

一种基于直接计算高阶奇异积分的断裂力学双边界积分方程分析法

相对于有限元法,边界单元法在求解断裂问题上有着独特的优势,现有的边界单元法中主要有子区域法和双边界积分方程法.采用一种改进的双边界积分方程法求解二维、三维断裂问题的应力强度因子,对非裂纹边界采用传统的位移边界积分方程,只需对裂纹面中的一面采用面力边界积分方程,并以裂纹间断位移为未知量直接用于计算应力强度因子.采用一种高阶奇异积分的直接法计算面力边界积分方程中的超强奇异积分;对于裂纹尖端单元,提供了三种不同形式的间断位移插值函数,采用两点公式计算应力强度因子.给出了多个具体的算例,与现存的精确解或参考解对比,可得到高精度的计算结果.      

汶川地震诱发文家沟巨型滑坡 碎屑流基本特征及成因机制初步分析


2008年5月12日,汶川M80地震在四川省绵竹市清平乡文家沟内诱发一巨型滑坡。通过现场调查得知,滑坡前后缘高差455m,厚度20～30m,滑面为基岩层面,初始方量2750×107m3。滑体在运动中转化为碎屑流。滑坡-碎屑流总的水平运动距离为4022m,垂直运动距离为1443m,遗留的堆积物体积达5×107m3。滑坡距映秀—北川断裂仅36km,位于其下盘,地震烈度达XI度。滑坡导致文家沟中48人遇害,并形成一条完整的地震次生地质灾害链。初步分析表明滑坡启动速度快,滑坡向碎屑流转化过程明显、地点明确。碎屑流运动过程复杂,伴有强烈的“气垫效应”和“前缘气浪冲击效应”。作者认为,文家沟滑坡的高启动速度是长持时强烈地震动作用的结果,与山体的猛烈碰撞是导致滑体解体并转化为碎屑流的原因。     

Preface: Orbits around asteroid

One of the core issues in modern celestial mechanics is the orbital dynamics in the near-regime gravitational field of as- teroids, which provides deep insights into the mathematical nature of a class of nonlinear systems, and plays as a critical basis for in situ explorations of different science goals. Lots of efforts have been made to reveal the characteristics of orbital motion in the vicinity of asteroids, and to improve the skills of asteroid research in methodology.     

Effet de l'anisotropie élastique cristalline sur la distribution des facteurs de Schmid à la surface des polycristauxInfluence of crystalline elasticity anisotropy on Schmid factor distribution at the free surface of polycrystals

Experimental studies of the plasticity mechanisms of polycrystals are usually based on the Schmid factor distribution supposing crystalline elasticity isotropy. A numerical evaluation of the effect of crystalline elasticity anisotropy on the apparent Schmid factor distribution at the free surface of polycrystals is presented. Cubic elasticity is considered. Order II stresses (averaged on all grains with the same crystallographic orientation) as well as variations between averages computed on grains with the same crystallographic orientation but with different neighbour grains are computed. The Finite Element Method is used. Commonly studied metals presenting an increasing anisotropy degree are considered (aluminium, nickel, austenite, copper). Concerning order II stresses in strongly anisotropic metals, the apparent Schmid factor distribution is drifted towards small Schmid factor values (the maximum Schmid factor is equal to 0.43 instead of 0.5) and the slip activation order between characteristic orientations of the crystallographic standard triangle is modified. The computed square deviations of the stresses averaged on grains with the same crystallographic orientation but with different neighbour grains are a bit higher than the second order ones (inter-orientation scatter). Our numerical evaluations agree quantitatively with several observations and measures of the literature concerning stress and strain distribution in copper and austenite polycrystals submitted to low amplitude loadings. Hopefully, the given apparent Schmid factor distributions could help to better understand the observations of the plasticity mechanisms taking place at the free surface of polycrystals. To cite this article: M. Sauzay, C. R. Mecanique 334 (2006).     

Nonlinear behavior of matrix-inclusion composites under high confining pressure: application to concrete and mortar

This paper is devoted to a micromechanics-based simulation of the response of concrete to hydrostatic and oedometric compressions. Concrete is described as a composite made up of a cement matrix in which rigid inclusions are embedded. The focus is put on the role of the interface between matrix and inclusion which represent the interfacial transition zone (ITZ). A plastic behavior is considered for both the matrix and the interfaces. The effective response of the composite is derived from the modified secant method adapted to the situation of imperfect interfaces. To cite this article: T.H. Le et al., C. R. Mecanique 336 (2008).     

Dynamics and Breakup of Viscoelastic Liquids (A Review)

The phenomena of hydrodynamic breakup of liquid jets, drops, films, bridges, and filaments are reviewed for liquids with viscoelastic properties. The reasons for breakup are capillary instabilities, collisions with rigid obstacles, and other forms of dynamic action. The relationship between the properties of the liquids and the features of the breakup process is discussed.     

Flow stability analysis and excitation using pulsating jets

Classical flow stability applied to transition from laminar to turbulent flow may also describe the behavior of vorticity fluctuations created by a pulsating jet placed along a solid boundary. A numerical laminar flow experiment involving a pulsating jet placed along the surface of a duct with flow separation downstream, resulted in eliminating most part of the separated flow region. Applying the same approach to a turbulent flow, it was possible to develop a turbulent stability flow formulation and apply successfully turbulent pulsating jet flow separation control. To cite this article: D. Skamnakis, K. Papailiou, C. R. Mecanique 333 (2005).