轴向运动弦线的纵向振动及其控制

综述轴向运动弦线纵向振动及其控制问题的研究进展.多种工程 系统如动力传送带、磁带、纸带、纺织纤维、带锯、空中缆车索道等均 涉及轴向运动弦线的纵向振动.对线性模型而言,除早期结果外,总结了 运动弦线的模态分析、具有复杂约束和耦合的运动弦线振动和运动弦线 参数振动的近期研究.对非线性模型而言,提出了轴向运动弦线大幅纵向 振动的运动微分方程,概述了离散化和直接近似解析分析、用黏弹性材 料模型化阻尼机制和动力传输系统的耦合振动研究的新进展.讨论了轴 向运动弦线振动主动控制的研究现状,包括能控性和能观性,控制分析的 频域方法和能量方法,振动的自适应控制和非线性振动的控制.最后指出 该研究方向今后需要研究的若干重要问题,包括运动弦线的非线性动力学 行为、黏弹性运动弦线的振动、含运动弦线的混杂系统的控制和轴向运 动弦线非线性振动的控制.     

轴向运动弦线横向振动的频域分析

轴向运动弦线是多种工程系统的模型。为明确轴向运动横向振动的频域特性，及探索频域方法的应用特点．本文用频域方法分析轴向运动弦线的横向振动。基于轴向运动弦线横向振动方程和边界条件．通过Laplace变换导出频率域中的控制方程，并将该控制方程和边界条件用状态变量表示。由状态空间中的控制方程导出特征方程，从而求出固有频率。由轴向运动弦线的矩阵函数计算得到系统的传递函数，然后用留数定理计算传递函数的Laplace逆变换．这样就可以得到时域响应。最后分析了轴向运动弦线的横向共振，若简谐外激励的频率与系统固有频率相同，系统响应将随时间无限增加。     

时变张力作用下轴向运动梁的分岔与混沌

轴向运动结构的横向参激振动一直是非线性动力学领域的研究热点之一. 目前研究较多的是轴向速度摄动的动力学模型,参数激励由速度的简谐波动产生. 但在工程应用中,存在轴向张力波动的运动结构较为广泛,而针对轴向张力摄动的模型研究较少. 本文研究了时变张力作用下轴向变速运动黏弹性梁的分岔与混沌. 考虑随着时间周期性变化的轴向张力,计入线性黏性阻尼,采用Kelvin模型的黏弹性本构关系,给出了梁横向非线性 振动的积分--偏微分控制方程. 首先应用四阶Galerkin截断方法将控制方程离散化,然后采用四阶Runge-Kutta方法计算系统的数值解,进而确定其动力学行为. 基于梁中点的横向位移和速度的数值结果,仿真了梁沿平均轴速、张力摄动幅值、张力摄动频率以及黏弹性系数变化的倍周期分岔与混 沌运动,并且通过计算系统的最大李雅普诺夫指数来识别其混沌行为. 结果表明:较小的平均轴速有助于梁的周期运动,梁在临界速度附近容易发生倍周期分岔与混沌行为. 随着张力摄动幅值的增大,梁的振动幅值的混沌区间不断增大. 较小的黏弹性系数和张力摄动频率更容易使梁发生混沌运动. 最后,给出时程图、频谱图、相图以及Poincaré 映射图来确定梁的混沌运动.      

时变张力作用下轴向运动梁的分岔与混沌

轴向运动结构的横向参激振动一直是非线性动力学领域的研究热点之一.目前研究较多的是轴向速度摄动的动力学模型,参数激励由速度的简谐波动产生.但在工程应用中,存在轴向张力波动的运动结构较为广泛,而针对轴向张力摄动的模型研究较少.本文研究了时变张力作用下轴向变速运动黏弹性梁的分岔与混沌.考虑随着时间周期性变化的轴向张力,计入线性黏性阻尼,采用Kelvin模型的黏弹性本构关系,给出了梁横向非线性振动的积分—偏微分控制方程.首先应用四阶Galerkin截断方法将控制方程离散化,然后采用四阶Runge-Kutta方法计算系统的数值解,进而确定其动力学行为.基于梁中点的横向位移和速度的数值结果,仿真了梁沿平均轴速、张力摄动幅值、张力摄动频率以及黏弹性系数变化的倍周期分岔与混沌运动,并且通过计算系统的最大李雅普诺夫指数来识别其混沌行为.结果表明:较小的平均轴速有助于梁的周期运动,梁在临界速度附近容易发生倍周期分岔与混沌行为.随着张力摄动幅值的增大,梁的振动幅值的混沌区间不断增大.较小的黏弹性系数和张力摄动频率更容易使梁发生混沌运动.最后,给出时程图、频谱图、相图以及Poincaré映射图来确定梁的混沌运动.     

轮对系统的Hopf分岔研究

针对轮对系统中的非线性动力学问题, 本文基于Hopf分岔代数判据得到考虑陀螺效应的轮对系统Hopf分岔点解析表达式, 即轮对系统蛇形失稳的线性临界速度解析表达式. 基于分岔理论得到轮对系统的第一、第二Lyapunov系数表达式, 并结合打靶法分别得到不同纵向刚度下, 考虑陀螺效应与不考虑陀螺效应的轮对系统分岔图. 通过对比有无陀螺效应的轮对系统分岔图发现, 在同一纵向刚度下, 考虑陀螺效应的轮对系统线性临界速度和非线性临界速度均大于不考虑陀螺效应的轮对系统, 即陀螺效应可以提高轮对系统的运动稳定性. 基于Bautin分岔理论, 以纵向刚度和纵向速度作为参数, 分别得到考虑陀螺效应和不考虑陀螺效应的轮对系统, 从亚临界Hopf分岔到超临界Hopf分岔, 再从超临界Hopf分岔到亚临界Hopf分岔的迁移机理拓扑图. 通过对比有、无陀螺效应的轮对系统Bautin分岔拓扑图发现, 陀螺效应将改变轮对系统的退化Hopf分岔点, 但对于轮对系统Bautin分岔拓扑图的影响不大.      

轴向移动局部浸液单向板的1:3内共振分析

考虑单向板的轴向速度、轴向张力、流固耦合作用以及阻尼等因素, 基于由 von Kármán薄板大挠度方程得到的轴向移动局部浸液单向板的非线性振动方程, 研究了外激励作用下单向板在1:3内共振情况时的非线性振动特性. 首先利用Galerkin法对非线性振动方程离散化, 然后分别应用数值法和近似解析法对离散后模态方程组进行求解, 获得了系统内共振情况下复杂的幅频特性曲线, 并讨论了周期解的稳定性. 最后研究了1:3内共振系统平均方程组的运动分岔现象.     

宏观交通流模型的余维2分岔分析

交通流特性是混合交通流建模的一个重要因素. 交通流模型中的分岔现象是导致复杂交通现象的因素之一. 交通流的分岔, 涉及复杂的动力学特征且研究较少. 因此, 提出了一个最优速度模型来研究驾驶员记忆对驾驶行为的影响. 基于带有记忆的最优速度连续交通流模型, 利用非线性动力学, 分析和预测了复杂交通现象. 推导了鞍结 (LP) 分岔存在条件, 并通过数值计算得到了余维1 Hopf (H) 分岔、LP分岔和同宿轨 (HC) 分岔以及余维2广义Hopf (GH) 分岔、尖点 (CP) 分岔和Bogdanov-Takens (BT) 分岔等多种分岔结构. 根据双参数分岔区域的特点, 研究了记忆参数对单参数分岔结构的影响, 分析了不同分岔结构对交通流的影响, 并用相平面描述了平衡点附近轨迹的变化特征. 选择Hopf分岔和鞍结分岔作为密度演化的起点, 描述了均匀流、稳定和不稳定的拥挤流以及走走停停现象. 结果表明, 驾驶员记忆对交通流的稳定性有重要影响; 动力学行为很好地解释了交通拥堵现象; 考虑余维2分岔的影响, 能更好地理解交通拥堵产生的根源, 并为制定有效抑制拥堵的方法提供一定的理论依据.      

轴向运动弦线横向振动的线性反馈控制

研究轴向运动弦线和作动器组成的耦合系统的横向振动控制。此系统被作动器分成未控和受控两部分,通过作用在作动器上的控制力对受控部分的横向振动进行控制。采用能量方法获得反馈控制规律,得到两种控制力,并用半群理论证实受控弦线横向振动的渐近稳定性和指数稳定性。在初始扰动和激励力作用下,通过数值仿真证实控制规律的有效性。     

轮对非线性动力学模型稳定性和分岔特性研究

为了探究轮对系统的横向失稳问题,考虑了陀螺效应和一系悬挂阻尼的影响作用,建立非线性轮轨接触关系的轮对动力学模型,研究轮对系统的蛇行稳定性、Hopf分岔特性及迁移转化机理.通过稳定性判据获得了轮对系统失稳临界速度.采用中心流形定理和规范型方法对轮对动力学模型进行化简,得到与轮对系统分岔特性相同的一维复变量方程,理论推导求得轮对系统的第一Lyapunov系数的表达式,根据其符号即可判断轮对系统的Hopf分岔类型.讨论了不同参数对轮对系统Hopf分岔临界速度的影响,探究了轮对系统的超临界、亚临界Hopf分岔域在二维参数空间的分布规律.利用数值模拟得到轮对系统的3种典型Hopf分岔图,验证了轮对系统超临界、亚临界Hopf分岔域分布规律的正确性.结果表明,轮对系统的临界速度随着等效锥度的增大而减小,随着一系悬挂的纵向刚度和纵向阻尼的增大而增大,随着纵向蠕滑系数的增大呈先增大后减小.系统参数变化会引起轮对系统Hopf分岔类型发生改变,即亚临界与超临界Hopf分岔相互迁移转化.轮对系统Hopf分岔域在二维参数空间的分布规律对于轮对系统参数匹配和优化设计具有一定的指导意义.     

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

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

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

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

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

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

Preface: Current research progress on mechanics of high speed rail

The construction and operation of high-speed rail（HSR）grid within the past two decades in China,in terms of the scale,may be possibly comparable to any national-wide construction in the history of China,even the Great Wall.By counting railways with commercial train service at the speed of200 km/h,China has the world’s longest HSR network with over 19 369.8 km of track in service today and this number     

ACTA MECHANICA SINICA 2014 ANNUAL SUBJECT INDEX

正~~     

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