完整系统的多刚体系统离散时间传递矩阵法研究

基于多刚体系统离散时间传递矩阵法，采用提高计算精度的方法，研究具有大运动、非线性特征的完整系统在平面、空间中的动力学响应。提出了对部分变量重新赋值的违约修正方法，计算机仿真表明了其有效性。多刚体系统离散时间传递矩阵法不须进行违约修正，体现了该方法建模灵活性较强、程式化程度较高的优点。     

控制一类磁性刚体航天器的混沌姿态运动

本文研究一类磁性航天器的混沌姿态运动及其控制,建立了在近地球赤道面圆轨道上运动受万有引力矩、磁力矩作用磁性刚体航天器姿态运动的动力学方程。采用时间历程、Poincare截面、Lyapunov指数和功率谱对系统的动力学行为进行数值识别,结果表明随着磁场参数的增大系统动力学行为由准周期环面破裂而出现混沌。利用输入－输出反馈精确线性化的方法将航天器的混沌姿态控制运动控制为姿态静止和按给定的周期规律运动,数值结果表明该控制方法的有效性。     

多刚体系统分离策略及释放动力学研究

紧密连接的多刚体系统可在脱离运载航天器后在轨自主分离,无需多次利用航天器发射装置或在航天器中安装多个发射装置进行分离释放,从而有效提高运载航天器空间利用率, 简化分离释放操作和降低碰撞风险.本文针对多刚体系统的在轨分离释放问题, 研究在轨分离策略及释放过程动力学.首先, 考虑刚体相对运动及姿态变化,基于虚功原理及自然坐标方法建立单个刚体的动力学模型.考虑多刚体系统在轨分离释放阶段的轨道运动和连接约束变化,计入分离时刚体间的相互作用,利用拉格朗日乘子法获得含连接约束的非线性动力学模型. 考虑到实际工程应用,在多刚体系统分离释放阶段,通过安装在刚体间每个接触表面4个角上的弹射装置实现自主分离. 其次,为保证分离过程中刚体之间无碰撞发生, 规划了多刚体系统的分离时序,并基于不同弹射方向及分离顺序设计了两种分离释放方案. 最后,通过算例研究分析了在轨分离释放过程中刚体的非线性动力学行为,验证了分离释放方案的有效性.      

刚体平面运动学问题的复数解法

刚体平面运动学问题的复数解法金明德（南京气象学院基础科学系，南京２１００４４）１引言在理论力学教材中，讨论刚体作平面运动的问题，总是将刚体的平面运动看成随基点的平动和绕基点转动的合成。因此，求刚体上任一点速度的方法也采用基点法（或由基点法派生的速度瞬...     

角速度张量及其应用

本文用旋转张量表示绕定点转动的刚体的角速度。给出计算刚体绕体轴系列、固定轴系列或体轴与固定轴交叉出现的一系列转动时其合成运动的角速度的统一公式。     

多刚体系统动力学的正则方程及数值解法

本文研究了多刚体系统动力学的正则方程,并由此论证了多刚体系统动力学方程的运动稳定性态.在比较不同数值积分方法的基础上,得出结论:基于正则方程的数值积分方法能较好地保证计算精度.     

汽车悬架闭环多体系统运动学分析

本文将多刚体系统动力学的R-W方法应用于具有闭合链系统的汽车悬架运动分析.首先通过假想切除铰链将闭环系统简化为树系统,然后引入约束方程使减缩系统的运动符合原系统的实际情况,得到运动学封闭解.文中给出了计算机程序框图及汽车两种典型悬架的分析计算结果.工程应用的实践表明,用此程序进行汽车悬架运动分析,所得结果精确,计算迅速,可方便地进行悬架设计方案的比较和参数的选择.     

树形带球铰的多刚体系统动力学方程的建立方法

本文提出了一种建立树形带球铰多刚体系统动力学方程的新方法,运用约束系统动力学研究成果,提出了广义刚体的概念以代替多刚体系统的子系统,并借助于广义刚体的不断扩充,求得了多刚体系统动力学方程建立的递推方法,该方法简单直观,几何概念清楚,并允许多刚体系统的扩充,且便于进行计算机符号处理。     

完全笛卡尔坐标描述的多体系统动力学

用完全笛卡尔坐标描述多体系统的运动学和动力学在提高计算效率方面有突出优点．导出用完全笛卡尔坐标表示的刚体及多体系统的动量和动量矩的解析式，给出与之对应的广义惯量矩阵概念，建立无力矩状态下用完全笛卡尔坐标描述的多体系统动力学的一阶微分方程组，用于多体航天器的姿态运动分析     

燃料消耗下充液航天器等效动力学建模与分析

在轨航天器贮腔内的液体可能表现出多种不同的运动模式, 主要包括液体相对于贮腔的整体性刚体运动、自由液面横向晃动、液体起旋后逐步发生明显的旋转晃动及液体自旋运动; 复合三自由度刚体摆晃动模型能够较为全面地描述这些液体运动模式, 同时为研究起旋阶段的液体晃动动力学问题提供了有效手段. 本文对非线性液体晃动刚体摆复合模型作进一步发展, 考虑模型等效参数随贮腔充液比的变化, 提出了变参数的刚体摆复合模型, 该模型适用于研究燃料消耗下非线性晃动类充液航天器大范围运动耦合动力学问题. 采用刚体摆复合模型对球形贮腔内的液体晃动进行等效后, 基于混合坐标意义下的拉格朗日方程推导了一类充液航天器轨道-姿态-晃动全耦合的动力学方程组, 并展开了充液航天器大角度三轴稳定姿态机动和零冲量轨道机动仿真以及航天器耦合动力学响应特性分析. 研究表明: 液体相对于贮腔的运动会造成航天器主刚体位置发生偏移, 当航天器在执行零冲量机动时, 燃料消耗会造成航天器的轨道平动速度无法收敛到零; 贮腔偏心布放时, 航天器在执行轨道机动过程中贮腔内液体易发生剧烈而且形式复杂的晃动行为, 进而可能造成航天器刚体运动的不稳定.      

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

正~~     

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

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

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

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

人造地震动合成中的位移误差分析

三角函数级数法是合成人工地震动常用算法之一,但是通过对加速度积分求取位移时,却存在与零线漂移相类似的位移漂移现象。