基于Reynolds边界的滑动轴承动力学系数的计算及应用

运用有限元方法,在不需要额外求解Reynolds方程的情况下,求解了具有Reynolds边值条件的流体润滑问题,使得同时完成动力积分过程中非线性油膜力及影响Floquet乘子求解的油膜力Jacobian矩阵的计算成为可能;运用打靶法及预估-校正和打靶法相结合的延续算法考察了轴承-转子系统的非线性不平衡响应及其随轴承设计参数改变而出现的分岔现象,实现了计算量的有效减少。     

基于可观测状态的轴承-转子系统周期解计算及稳定性分析

分析了轴承-转子系统的稳定性和分岔,基于系统可观测状态信息给出1种求解系统周期解及识别周期解稳定性的方法,同时将该方法与Floquet理论相结合分析系统周期解的稳定性及失稳分岔形式,将转速作为分岔参数分析系统响应的周期、拟周期、多解共存和跳跃现象.结果表明,采用该方法计算系统周期解及稳定性时,利用系统可观测稳态和瞬态信息,即可求解出系统Jacobian矩阵而无需实时求解轴承非线性油膜力的Jacobian矩阵.与传统PNF方法相比,该方法不仅具有很高的精度而且可以节约计算量,同时可以预测追踪随控制参数变化的系统周期解及其稳定性,可用于指导轴承-转子系统的非线性动力学设计.     

两自由度机械臂λ^2=1下的Hopf分岔与混沌研究

研究了一类周期系数力学系统因周期运动失稳而产生Hopf分岔及混沌问题.首先根据拉格朗日方程给出了该力学系统的运动微分方程,并确定其周期运动的具有周期系数的扰动运动微分方程,再根据Floquet理论建立了其给定周期运动的Poincaré映射,根据该系统的特征矩阵有一对复共轭特征值从-1处穿越单位圆情况,分析该Poincaré映射不动点失稳后将发生次谐分岔、Hopf分岔、倍周期分岔,而多次倍周期分岔将导致混沌.并用数值计算加以验证.结果表明,随着分岔参数的变化,系统的周期运动可通过次谐分岔形成周期2运动,进而发生Hopf分岔形成拟周期运动,并再次经次谐分岔、倍周期分岔形成混沌运动.     

磁浮轴承-转子系统非线性动态特性分析

考虑非线性电磁力对刚性Jeffcott转子系统的影响，采用Hopf分岔理论及CPNF法对系统平衡点解和周期解进行研究，数值仿真得到系统Jacobi矩阵特征值、轴心轨迹图和Poincare映射图。转子运动呈现Hopf分岔、倍周期分岔及拟周期运动等复杂的非线性动力学特征，其结果可为磁浮轴承-转子系统设计和运行状态控制提供理论依据。     

轴承-转子系统在基础运动作用时的非线性分析

轴承-转子系统的非线性特性及其在基础运动作用时的响应,是离心机设备设计阶段必须考虑的。本文使用简化的多自由度转子模型进行模拟分析,运动方程考虑了非线性的油膜润滑轴承模型。应用自适应时间步长的Runge-Kutta-Felburg法求解微分运动方程组,将人造的正弦波加速度作为基础运动输入系统,使用Poincaré图、分岔图和瀑布图分别考察了垂直放置转子在有无基础运动作用时的动力学性质。快速傅里叶变换在频域内揭示了转动频率与基础振动频率之间的组合共振现象。计算的结果不仅给出了泵转子自身的非线性性质,也展示了泵转子在基础运动作用时的组合共振。     

线性碰振系统周期解擦边分岔的一类映射分析方法

擦边分岔是碰振机械系统的一种重要分岔行为. 以固定相位面作为Poincaré截面， 建立了线性碰振系统单碰周期$n$运动的Poincaré映射. 通过分析该映射，得到了系统 发生擦边分岔的条件和分岔方程，并以单自由度碰振系统为实例验证了分析结果的正确性. 该方法不仅可以计算线性碰振系统擦边分岔的参数值，还可以计算系统的任意周 期$n$解的分岔参数值.     

高速涡轮增压器轴承-转子系统不平衡响应及稳定性分析

汽车涡轮增压器广泛采用浮环轴承支承的小型轻质转子系统,以实现100 000~300 000 r/min的工作转速,提高发动机功率和动力性能,并降低燃油消耗和排放. 在此超高速工况下,动压油膜的强非线性作用和转子固有的不平衡效应使该系统呈现出复杂的动力学现象,其中油膜涡动、振荡、跳跃、倍周期分岔和混沌等非线性动力学行为对增压器的健康运转意义重大,因而备受关注. 本文作者从摩擦学动力学耦合的角度出发,基于流体动压轴承润滑理论和有限差分法计算非稳态油膜压力,结合达朗贝尔原理和传递矩阵法建立了转子离散化动力学方程,提出了一种由双油膜浮环支承的涡轮增压器转子系统动力学模型,并从转子轨迹、轴承偏心率、频谱响应、庞加莱映射和分岔特性等方面比较分析,描述了该非线性轴承-转子系统的不平衡效应及油膜失稳特征. 结果表明:转子一般在相对低速下作稳定的单周期不平衡振动,在高转速下其被油膜失稳引起的次同步涡动所抑制,但不平衡量的增加可阻碍转子以拟周期运动通向混沌运动的路径;适当不平衡补偿下,由于内、外油膜间交互的非线性刚度和阻尼作用,在油膜失稳区间之间的中高速区会出现适合增压器健康运转的稳定区间.      

可倾瓦径向滑动轴承完整动特性系数的分析模型

基于可倾瓦径向滑动轴承瓦块的扰动特性,提出了计算轴承完整动力系数的数学解析模型。在由单块瓦和轴颈构成的子系统上建立局部动坐标参考系,全局广义位移向量可以通过简练的步骤转换为局部动坐标系下轴颈的位移向量,利用求解固定瓦轴承动力特性的方法求得的局部动坐标系下的油膜力又可以精确地转换为全局坐标系下的表达形式,全局坐标系下的油膜力向量关于广义位移和广义速度的Jocabian矩阵形成轴承的完整动力特性系数;利用Newton-Raphson方法同时求解瓦块和轴颈在给定的静态载荷下的平衡位置,其中所需用到的系数矩阵恰好为轴承油膜力关于广义位移的Jocabian矩阵的负值,在得到平衡位置的同时可以获得轴承完整的刚度和阻尼矩阵。数值算例证明了此方法的有效性。     

三自由度齿轮系统的双参数分岔与全局特性研究

以三自由度齿轮系统为研究对象,通过构造参数平面内不同运动类型的边界线算法,得到了系统在参数平面内的分岔曲线。为了判断分岔曲线的分岔类型,构造了三自由度齿轮系统Poincaré映射的Jacobi矩阵及Floquet乘子算法。结合系统的分岔图、最大Lyapunov指数图(TLE)、相图、Poincaré映射图和Floquet理论,讨论了双参数平面上系统的分岔特性以及参数平面内系统动力学特性的演变,并利用胞映射法对系统随啮合频率变化下的全局动力学特性进行了研究。结果表明:系统在参数平面k-ξ33内存在倍化分岔曲线、鞍结分岔曲线、Hopf分岔曲线等;阻尼系数越大,综合误差越小,系统运动越稳定;鞍结分岔对系统的全局稳定性影响较大,而Hopf分岔对系统的全局稳定性影响较小。研究结果可为齿轮系统设计和参数选择提供理论依据,研究方法也适用于其它非线性系统的双参数分岔分析。     

复杂非线性转子—轴承系统动力特性数值分析

研究非线性高维复杂转子－轴承系统的动力特性。针对系统的局部非线性特征,给出了一种降阶及配套动力积分方法。降阶系统仍保持局部非线性特征,非线性响应数值积分所需的迭代只需在局部非线性的维数上执行。对于油膜力无封闭解的实际轴承,采用变分不等方程有限元法求解Reynolds边值问题,使得油膜力及其Jacobian矩阵的计算变得非常简单明了且与具有协调一致的精度。应用上述方法计算分析了一双跨、椭圆轴承－转子系统的不平衡响应,数值结果展现了系统丰富复杂的非线性现象。     

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

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

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

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