参-强激励联合作用下输流管的分岔和混沌行为研究

研究输送脉动流的两端固定输流管道在其基础简谐运动激励下的分岔和混沌行为,考虑管道变形的几何非线性和管道材料的非线性因素,推导了系统的非线性运动方程,并应用Galerkin方法对其进行了离散化处理。通过采用数值模拟方法,对系统的运动响应进行仿真,重点探讨了流体平均流速、流速脉动振幅以及基础简谐运动激励振幅对系统动态特性的影响。结果表明,系统在不同的参数下会发生围绕不同平衡点的周期和混沌等运动,并在系统中发现了两条通向混沌运动的途径:倍周期分岔和阵发混沌运动。     

磁场中轴向运动条形板强非线性超谐共振及混沌运动

针对磁场环境中周期外载作用下轴向运动导电条形板的非线性振动及混沌运动问题进行研究。应用改进多尺度法对横向磁场中条形板的强非线性振动问题进行求解,得到超谐波共振下系统的分岔响应方程。根据奇异性理论对非线性动力学系统的普适开折进行分析,求得含两个开折参数的转迁集及对应区域的拓扑结构分岔图。通过数值算例,分别得到以磁感应强度、轴向拉力、激励力幅值和激励频率为分岔控制参数的分岔图和最大李雅普诺夫指数图,以及反映不同运动行为区域的动力学响应图形,讨论分岔参数对系统呈现的倍周期和混沌运动的影响。结果表明,可通过相应参数的改变实现对系统复杂动力学行为的控制。     

压电复合材料层合梁的分岔、混沌动力学与控制

研究了简支压电复合材料层合梁在轴向、横向载荷共同作用下的非线性动力学、分岔和混沌动力学响应. 基于vonKarman理论和Reddy高阶剪切变形理论,推导出了压电复合层合梁的动力学方程. 利用Galerkin法离散偏微分方程,得到两个自由度非线性控制方程,并且利用多尺度法得到了平均方程. 基于平均方程,研究了压电层合梁系统的动态分岔,分析了系统各种参数对倍周期分岔的影响及变化规律. 结果表明,压电复合材料层合梁周期运动的稳定性和混沌运动对外激励的变化非常敏感,通过控制压电激励,可以控制压电复合材料层合梁的振动,保持系统的稳定性,即控制系统产生倍周期分岔解,从而阻止系统通过倍周期分岔进入混沌运动,并给出了控制分岔图.      

黏弹性传动带1:3内共振时的周期和混沌运动

研究了参数激励作用下黏弹性传动带在1:3内共振时的周期解分岔和混沌动力学. 同时考虑传动带的线性外阻尼因素和材料内阻尼因素. 首先建立了具有线性外阻尼情况下的黏弹性传动带平面运动时的非线性动力学方程, 黏弹性材料的本构关系用Kelvin模型描述. 然后考虑黏弹性传动带的横向振动问题, 利用多尺度法和Galerkin离散法得到黏弹性传动带系统在1:3内共振时的平均方程. 最后利用数值模拟方法研究了黏弹性传动带系统的周期振动和混沌动力学, 得到了系统在不同参数下的混沌运动. 数值模拟结果说明黏弹性传动带系统存在周期分岔, 概周期运动及混沌运动.     

不平衡量对非线性多转子系统动力特性的影响

用近代非线性动力学理论分析了弹性支承有间隙和摩擦的非线性刚性多转子系统的复杂运动．建立了支座有间隙和有摩擦的弹性支承的力学模型．导出了这类多转子系统的运动微分方程组．用数值方法得到系统在某些参数区域内的轴心轨迹图，Poincare映射图和分岔图等．以转子不平衡量为控制参数讨论了进出混沌区的不同路径和系统各种形式的拟周期，倍周期和混沌运动．分析结果为定性地改善转子系统的稳定运行状态提供了理论依据．     

叶片振动影响下双盘转子-轴承系统的稳定性与分岔

研究叶片与转子-轴承系统的耦合非线性振动,建立了一个带叶片的双盘转子-轴承系统的非线性动力学模型,其中包含一个弹性转轴、两个滑动轴承、两个刚性圆盘和两组弹性叶片.为了分析叶片的惯性影响,将其简化为单摆模型.采用4阶Runge-Kutta法进行了数值模拟,并利用分岔图、三维谱图、轴心轨迹和Poincaré映射图等方法分析了系统的非线性动力学特性.研究发现,随着转速的变化,系统响应演化出了倍周期运动、概周期运动、混沌运动和倍周期分岔等典型的非线性动力学行为.在与忽略了叶片振动的转子系统对比后发现,叶片振动使转子发生混沌运动的转速区域增大.在某些参数条件下,采用不同的叶片刚度,叶片振动可能引起转子系统产生混沌运动.     

悬臂输流管道的运动分岔现象和混沌运动

研究受约束悬臂输流管道的稳定性和运动分岔问题.在静态与动态失稳区域边界上的一个交叉点附近,用理论分析的方法详细研究了该系统可能发生的复杂运动和运动分岔现象.在动态失稳区域内发现了输流管道的概周期运动和由于概周期运动环面破裂而导致混沌的现象.理论分析结果与数值模拟结果相吻合.     

悬臂输流管道在基础激励与脉动内流联合作用下的参激振动

首先建立了悬臂输流管道在基础激励与脉动内流联合作用下的运动方程;然后基于Galerkin法研究了该系统的非线性动力学行为,分析了系统运动状态随激励频率和相位差的变化,以及混沌百分比随频率比(基础激励频率与脉动频率之比)和相位差的变化。结果表明,无论以激励频率还是以相位差为分岔参数,系统都具有多种形式的动态响应,包括周期运动、概周期运动和混沌运动,但进入和脱离混沌的途径不同。相位差和频率比对系统的混沌百分比有重要影响:相位差为π/2时系统混沌百分比最大;频率比为1时系统混沌百分比最小,频率比较小或较大时系统混沌百分比与只有基础激励时接近。     

非自治时滞反馈控制系统的周期解分岔和混沌

研究时滞反馈控制对具有周期外激励非线性系统复杂性的影响机理,研究对应的线性平衡态失稳的临界边界,将时滞非线性控制方程化为泛函微分方程,给出由Hopf分岔产生的周期解的解析形式．通过分析周期解的稳定性得到周期解的失稳区域,使用数值分析观察到时滞在该区域可以导致系统出现倍周期运动、锁相运动、概周期运动和混沌运动以及两条通向混沌的道路：倍周期分岔和环面破裂．其结果表明,时滞在控制系统中可以作为控制和产生系统的复杂运动的控制“开关”．     

外激励作用下亚音速二维壁板分岔及响应研究

研究了亚音速流中二维壁板在外激励作用下的分岔和响应问题. 采用Galerkin方法将非线性运动控制方程离散为常微分方程组. 采用Runge-Kutta数值方法进行了数值计算,研究了壁板系统非单周期区在参数空间的分布情况. 结果表明: 在参数空间中, 非单周期区和单周期区会交替出现; 在不同的单周期区内, 系统运动轨线也在有规律的变化; 系统由单周期运动进入混沌运动是经过一系列周期倍化分岔产生的.      

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

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

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

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