Melnikov方法在输流管混沌运动研究中的应用

对基础简谐运动激励下两端固定输流管道的混沌运动进行了研究，推导出了系统的运动方程，确定了系统存在的平衡点及其稳定性，计算出了未扰系统的同宿轨道，并利用Melnikov方法得到了系统发生混沌运动时参数需满足的临界条件，同时还利用相平面图和：Poincare映射等方法对管道的混沌运动进行了数值模拟，通过比较发现，由Melnikov方法确定的临界参数值要稍小于数值模拟中首次观察到混沌运动时对应的临界值。     

基于格子Boltzmann模型的SINE-GORDON方程的数值解法

我们给出了用于求解Sine-Gordon方程的格子Boitzmann模型．通过假设与时间导数有关的平衡态分布函数，得到了具有2阶误差的Sine-Gordon方程的格子Boitzmann算法．数值例子再现了在零点附近的扭结解。     

不确定车轨耦合系统辛随机振动分析

建立了轨道不平顺作用下具有不确定参数车轨耦合系统随机振动评估方法. 车辆系统采用物理坐标下多刚体系统模型,并应用高斯随机变量模拟车体、转向架和轮对一系、二系连接系统中动力学参数具有的不确定性. 采用无穷周期结构进行弹性轨道模拟,在哈密顿状态空间下建立了典型轨道子结构的状态运动方程,通过轮轨耦合关系建立了混合 物理坐标及辛模态坐标车轨耦合系统运动方程. 应用Hermite正交多项式展开得到了耦合系统动力响应相对于不确定性参数的控制方程. 由于利用轨道周期特性建模,所获得的控制方程有效地降低了方程维度. 轮轨接触处轨道不平顺载荷模拟为完全相干多分量平稳随机过程,推广和发展虚拟激励法建立了耦合系统随机振动受不确定动力学 参数影响的量化评估方法. 通过Monte Carlo数值模拟,验证了该方法在不确定参数变异很大时也能够保持较好的精度,具有一定的工程实用性.     

Sine-Gordon方程的截断系统的同宿轨道

研究Ｓｉｎｅ Ｇｏｒｄｏｎ方程的广义渐近惯性流形上的常微分方程组，证实了在一定参数条件下存在Ｗｉｇｇｉｎｓ［１］意义下的同宿轨道．计算表明，与Ｂｉｓｈｏｐ［２］用数值计算得到的Ｓｉｎｅ Ｇｏｒｄｏｎ方程产生混沌的参数值尚有差别，考虑到同宿出现参数值往往低于混沌出现参数值，故结果在定性上正确，而且改进了文［１］中的结果．     

长型浮置板轨道隔振系统理论分析(Ⅰ)——弥散曲线及临界速度

浮置板轨道结构是近年来减振降噪效果最好的轨道结构形式.将整个轨道结构看作一个隔振黑匣子,钢轨与浮置板分别用欧拉梁来模拟,推导动力方程,用傅氏变换将其转换到频率-波数域中,求得系统的弥散方程,并给出轨道结构临界速度计算方法.定义了反映浮置板轨道结构参数的质量比例ξ_(m)、刚度比例ξ_(h)及阻尼比例ξ_(c)等三个系数.在给定实际中板上结构物理参数的条件下,深入讨论了三个比例系数对于系统弥散曲线以及临界速度的影响,进而为浮置板轨道共振频率及临界速度的初步设计提供建议.     

长型浮置板轨道隔振系统理论分析(I)——弥散曲线及临界速度

浮置板轨道结构是近年来减振降噪效果最好的轨道结构形式.将整个轨道结构看作一个隔振黑匣子,钢轨与浮置板分别用欧拉梁来模拟,推导动力方程,用傅氏变换将其转换到频率-波数域中,求得系统的弥散方程,并给出轨道结构临界速度计算方法.定义了反映浮置板轨道结构参数的质量比例ξ_(m)、刚度比例ξ_(h)及阻尼比例ξ_(c)等三个系数.在给定实际中板上结构物理参数的条件下,深入讨论了三个比例系数对于系统弥散曲线以及临界速度的影响,进而为浮置板轨道共振频率及临界速度的初步设计提供建议.     

一类含平方、立方非线性质项1：2内共振两自由度系统的全局分岔

研究了一类含平方、立方非线性项的两自由度系统全局分岔。首先应用多尺度法求解其平均方程，然后通过一系列变换得到一个近似可积的两自由度系统。应用能量-相位准则，确定了在哈密顿共振时Silnikov轨道存在的条件。通过数值计算验证了此条件。     

sine—Gordon方程的多重脉动跳跃轨道

本文利用ｓｉｎｅ－Ｇｏｒｄｏｎ方程的广义渐近惯性流形上的常微分方程，证明了在一定参数条件下存在多重脉动跳跃轨道，本工作与先前工作（１－４）证实可用ｓｉｎｅ－Ｇｏｒｏｄｎ方程的有限维形式来讨论其的动力学行为。     

有界噪声激励下单摆-谐振子系统的混沌运动

研究了具有同宿轨道和周期轨道的可积单摆-谐振子系统在弱Hamilton摄动(即弱耦合摄动)和弱非Hamilton摄动(即阻尼和有界噪声微扰)下的混沌运动．用Melnikov方程预测Hamilton系统中可能存在混沌运动的参数域，并用Poincare截面验证解析结果．用数值方法计算了有阻尼与有界噪声激励下系统的最大Lyapun0V指数和Poincare截面，结果表明有界噪声在频率上的扩散减小了引发系统产生混沌运动的效应。     

利用无网格伽辽金法分析材料稳态蠕变问题

稳态蠕变是高温环境下材料的一个重要的考虑问题，它也是材料破坏的一种重要形式。由于存在划分网格，利用传统有限元法模拟稳态蠕变有一定的不足之处。作为一种新兴的数值模拟方法，无网格法不需要划分单元，只需将求解问题离散成独立的结点，计算过程可以局部细化。利用连续介质的虚功原理，将无网格伽辽金法应用于稳态蠕变的数值模拟，推导了稳态蠕变的无网格伽辽金法控制方程。利用罚参数来实现稳态蠕变的不可压缩条件和本质边界条件，能够保证求解过程中刚度矩阵的对称正定性。通过实例的计算结果表明，无网格伽辽金法在求解稳态蠕变时具有较高的计算精度，结果与理论解结果吻合，而且前后处理较为简单。     

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

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

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


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).