非线性模态的分类和新的求解方法

引入不可分偶数维不变流形的概念来定义非线性模态．在此基础上，揭示出了一种新的模态——耦合非线性模态，并对实际系统中各种可能的模态进行了分类．这种分类可能是新的构筑非线性模态理论的框架．用此方法构造非线性模态，得到的模态振子具有范式的形式，形式最简、却能反映原系统在平衡点附近的主要动力学行为，且易于得到非线性频率及非线性稳定性等方面的信息．不仅适用于分析一般的多自由度系统，还可用于分析奇数维系统；不仅可构造内共振系统的非耦合模态，还可用于构造内共振耦合模态．从掌握的资料看，以前的方法还不能解决上述所有问题     

连结子结构与自由界面模态综合法在非线性动力分析中的应用

本文提出了一种由线性连接元和非线性连接元组成的连接子结构，并将这种连接子结构用于自由界面的模态综合技术。利用非线性振动理论的渐近方法，求得经模态综合法降维后系统方程的近似解析解。从而将具有连接子结构的自由界面的模态综合技术推广应用到具有局部非线性的复杂结构系统的动力分析，为利用非线性振动理论的渐近方法及动力系统理论进一步研究高维非线性动力学系统的振动特性、分岔及混沌行为创造了一种新的途径。算例表明，该方法具有足够的精度。     

关于一类非线性阻尼模型控制问题的一个注记

利用模态方法进行柔性体控制，就必须作模态截断。在线性粘性阻尼模型情形，阻尼正比于频率，高阶模态会很快衰减。但对于非线性阻尼模型，情况就复杂了。结果表明，高阶模态响应不总能很快衰减，并依赖于初始条件、因为对于非线性阻尼模型的控制，如何进行模态截断，是一个需要谨慎处理的问题。     

基于假设模态加权的全局模态提取方法及其应用

假设模态法在单一梁、杆、索、板等柔性结构动力学建模中有广泛应用，但在处理组合结构振动问题时，常常因无法反映各部件之间的耦合作用使其应用受限。通过假设模态建立组合结构的近似动力学模型，利用近似模型求得系统的固有频率和相应的特征向量，据此可以有效地获得系统的全局模态。本文以跨中带有多个弹性支撑的简支梁为例，通过假设模态加权来提取系统的全局模态，从而建立系统的动力学模型。对系统进行固有特性分析的结果表明，通过假设模态加权可以方便地获得系统的全局模态；对系统动态响应分析的结果表明，采用本文提出的全局模态建立的非线性动力学模型可以有效地反映系统的非线性动力学特性。     

作大范围运动矩形板的动力学建模理论研究

研究了初应力法的作大范围运动矩形板的建模理论。根据连续介质理论,考虑应变-位移中的非线性项,用一致质量有限元法对柔性板进行离散,基于Jourdain速度变分原理导出定轴转动下大范围运动为自由的柔性板刚-柔耦合动力学方程。从其刚柔耦合动力学方程出发,考虑在大范围运动已知情况下的结构动力学方程。通过引入准静态概念,把其结构动力学方程转化为准静态方程。对纵向和横向变形节点坐标进行坐标分离,解出与纵向变形相关的准静态方程,得到准静态时的纵向应力表达式,从而获得附加刚度项;并对此非惯性系下作大范围运动柔性板的结构动力学方程进行数值仿真,验证了采用初应力法柔性板的动力学建模方法来计算经历大范围运动的不规则柔性板的动力学响应是可行的,体现了初应力法对柔性板建模的优越性。     

具有局部非线性动力系统周期解及稳定性方法

对于具有局部非线性的多自由度动力系统，提出一种分析周期解的稳定性及其分岔的方法该方法基于模态综合技术，将线性自由度转换到模态空间中，并对其进行缩减，而非线性自由度仍保留在物理空间中在分析缩减后系统的动力特性时，基于Ｎｅｗｍａｒｋ法的预估－校正－局部迭代的求解方法，与Ｐｏｉｎｃａｒé映射法相结合，推导出一种确定周期解，并使用Ｆｌｏｑｕｅｔ乘子判定其稳定性及分岔的方法     

有间隙折叠舵面的振动实验与非线性建模研究

针对折叠舵面内、外舵铰接处存在的间隙对地面振动响应的影响及间隙处的非线性建模方法展开研究.消除间隙,利用锤击法对线性折叠舵面进行模态实验,得到了前五阶模态参数;打开间隙,进行振动台扫频基础激励,实验结果表明间隙的存在会使结构的动力学响应产生非线性现象,如正反向扫描差异、跳跃、多谐波及频率漂移.非线性的影响主要体现在一阶弯曲模态上,激励量级的增大和间隙的减小均会使基频增大,且逐渐趋向于无间隙的结果,但对第二阶扭转模态的影响与第一阶相比较小.建立了折叠舵面的有限元模型.提出了一种适用于具有集中非线性的折叠机构的模型缩减方法,并对舵面进行了模态缩减.根据Hertz接触理论,用具有线性和3/2次刚度组合形式的非线性扭转弹簧来模拟铰接处的间隙和接触.通过比较锤击实验与数值计算得到的前四阶频率和振型对模型的线性部分进行验证.通过Bathe两子步隐式复合算法计算基础激励下非线性结构的动力学响应,得到的传递函数可以模拟实验中出现的频率变化特征,验证了连接处非线性建模方法的合理性.     

有间隙折叠舵面的振动实验与非线性建模研究

针对折叠舵面内、外舵铰接处存在的间隙对地面振动响应的影响及间隙处的非线性建模方法展开研究.消除间隙,利用锤击法对线性折叠舵面进行模态实验,得到了前五阶模态参数;打开间隙,进行振动台扫频基础激励,实验结果表明间隙的存在会使结构的动力学响应产生非线性现象,如正反向扫描差异、跳跃、多谐波及频率漂移.非线性的影响主要体现在一阶弯曲模态上,激励量级的增大和间隙的减小均会使基频增大,且逐渐趋向于无间隙的结果,但对第二阶扭转模态的影响与第一阶相比较小.建立了折叠舵面的有限元模型. 提出了一种适用于具有集中非线性的折叠机构的模型缩减方法,并对舵面进行了模态缩减.根据Hertz接触理论,用具有线性和3/2次刚度组合形式的非线性扭转弹簧来模拟铰接处的间隙和接触.通过比较锤击实验与数值计算得到的前四阶频率和振型对模型的线性部分进行验证.通过Bathe两子步隐式复合算法计算基础激励下非线性结构的动力学响应,得到的传递函数可以模拟实验中出现的频率变化特征,验证了连接处非线性建模方法的合理性.      

挠性航天器太阳翼全局模态动力学建模与实验研究

现代柔性航天器通常安装有大型太阳翼为其在轨运行提供所需动力. 航天器入轨后太阳翼展开并锁定成为铰链连接多板结构, 此类结构质量轻、跨度大、刚度低的特点使其低频振动和非线性振动问题越来越凸显. 分析和处理此类结构出现的复杂振动问题的关键在于建立系统精确的非线性动力学模型. 为此, 本文提出铰链连接多板结构解析全局模态的提取方法, 获取太阳翼的固有频率和解析函数表征的全局模态. 提出可变刚度的扭转弹簧等效模型, 考虑铰链非线性刚度及摩擦力矩等因素, 通过全局模态离散得到系统的低维高精度非线性动力学模型, 研究了太阳翼在周期激励作用下的非线性特性. 开展太阳翼地面振动实验研究, 采用锤击法获取系统模态, 利用振动台施加正弦扫频激励, 将物理实验结果与理论结果进行对比, 从而验证全局模态动力学建模方法的合理性与准确性. 结果表明, 铰链刚度等结构参数对系统固有特性的影响较大, 铰链的存在会使太阳翼的动态响应出现跳跃等非线性现象. 全局模态动力学建模方法能很好地解决多板结构在非经典边界下解析全局模态求解的困难, 系统全局模态反映的是系统各个部件弹性振动的真实模态, 所建立的动力学模型具有低维高精度的特点, 对于复杂组合结构非线性动力学建模具有重要的参考价值.      

空间网格结构动力分析的非线性模态方法

空间网格结构因自由度数多且无简化的力学模型,非线性动力分析通常要耗费大量时间．传统的非线性模态方法用于求解多高层结构的局部非线性问题已获得良好的效果,但对系统非线性问题的应用尚缺少研究．对比分析多高层结构和空间网格结构动力性能差异,指出网格结构动力非线性分析存在的问题．以主振型理论和切线刚度分离法为基础,将非线性模态方法用于几何非线性效应显著的空间网格结构动力分析．通过对运动方程的非线性恢复力进行拆分,形成线性表达形式,然后解耦到主振型所在的广义坐标系,以达到缩减自由度数量的目的．并通过实例验证非线性模态方法的高效性与适用性．     

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

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

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

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