非光滑多体系统最大Lyapunov指数的计算方法

结合非光滑力学与非线性动力学的理论给出非光滑多体系统最大Lyapunov指数的数值计算方法.应用第一类Lagrange方法建立具有单边约束非光滑多体系统变拓扑结构动力学方程;应用线性互补理论与方法检测非光滑事件;应用混沌同步方法计算Lyapunov指数.在提高计算速度方面,提出了采用二分搜索法加快搜索速度以及采用调制参数方法加快系统达到同步状态的速度.通过算例说明本文算法的可行性与有效性.      

基于混合Gramian的柔性航天器保结构平衡降阶

对于柔性航天器动力学方程的降阶问题,传统的内平衡降阶方法破坏了原方程的动力学结构.为了解决这一问题,提出了一种牺牲严格内平衡状态,而保留原系统动力学特性的近似内平衡降阶方法.该方法利用混合Gramian描述系统的模态价值,针对柔性航天器动力学方程的特殊形式,给出了混合Gramian的闭合解析解.该方法得到的降阶模型可以保留原系统包括二阶结构,对称和正定性在内的所有动力学特性.数值算例表明该方法不但可以达到与传统一阶内平衡方法相当的降阶精度,而且降阶速度有大幅度提升.      

空间站大型伸展机构柔性多体动力学数值计算

对用ｋａｎｅ方法建立的空间站大型伸展机构柔性多体动力学数学模型,从数值计算方法和软件实现方面进行了研究。文章采用局部附着浮动坐标系描述各构件的刚体大范围伸展运动;用截断的子结构正则模态描述各构件的柔性变形;用求解刚性方程的ＴＲＥＡＮＯＲ方法进行数值积分;用混合法对系统的微分－代数动力学方程组进行求解。还对空间站伸展机构多体系统的约束方程建立、独立坐标选取和动力学方程的ｓｔｉｆｆ病态问题进行了讨论。并针对空间站伸展机构在驱动展开过程中以及展开完成后的动力学响应进行了算例计算。     

航天器空间交会过程综合变轨策略研究

文章研究了空间交会远程导引阶段综合变轨问题的变轨策略方法和实用理论模型。文章首先简要描述了综合变轨的技术条件和策略方法 ;然后重点建立了多脉冲综合变轨的轨道动力学模型 ,并给出了数值求解方法 ;最后用自行开发的仿真程序进行了远程导引多脉冲综合变轨的算例仿真。仿真结果表明 ,该模型和算法具有工程实用性     

带约束多体系统Lyapunov指数的数值计算方法

应用第1类Lagrange方程建立的带约束多体系统动力学方程为非线性微分-代数方程组.利用增广法,将其转化成了常微分方程组,并表示成矩阵形式.根据方程的特点,给出了方程的Jacobi矩阵的具体表达式,提高了计算效率.在此基础上,给出了Lyapunov指数的数值计算方法,并通过对具体的树形和非树形多体系统进行Lyapunov指数数值计算,结合相图和Poincare映射对系统的动力学特性进行了分析,表明了该方法的可行性和有效性.      

柔性机器人动力学建模的一种方法

提出一种建立具有柔性关节的多柔杆机器人动力学方程的方法.首先给出了柔性关节及柔性臂杆的简化模型,然后利用Kane方法和假设模态法推导出完整的动力学方程并给出了递推公式.在此基础上利用一个算例研究了臂杆柔性与关节柔性对机器人动力学响应的影响.结果表明:所建立的动力学模型更接近于实际的柔性机器人系统;臂杆柔性与关节柔性都起着非常重要的作用,因此在柔性机器人的动力学建模与控制中应充分考虑它们的影响.      

卫星太阳阵展开的动力学方程

研究了卫星太阳阵展开的动力学问题，将卫星太阳阵的结构视为树形多体系统，利用给出的自由形式树形多刚体系统运动微分方程建立太阳阵系统的动力学方程，应用为该系统动力学方程编制的计算程序，可对太阳阵展开过程中卫星本体姿态的改变量以及太阳帆板在展开过程中的运动参数进行数值计算。     

多梁复合结构的动力学建模及全局解析振型分析

多梁复合结构广泛用于工程领域,其动力学特性对相关装备的工作状况有着重要影响,建立准确适用的动力学模型对后续的动力学响应参数影响规律研究有决定性作用。针对工程中常见的多梁复合结构,给出了位移场的一般表达式,利用哈密顿原理得出多梁复合结构的运动控制方程和边界条件,对其全局解析振型进行理论计算,给出了计算多梁复合结构数值振型的一般过程,基于数值分析和实验分析以L型梁和Z型梁为例利用频率相对误差、振型定性比较和MAC定量分析说明了全局解析振型的求解验证方法,利用伽辽金方法和全局解析振型给出了截断原系统得到离散系统运动方程的一般形式。提出的多梁复合结构的建模方法和全局解析振型求解验证方法,对工程中此类结构动力学响应的参数影响规律研究有重要意义。     

含单边非完整约束飞机滑跑的建模与仿真方法

当飞机在地面上滑行时,其运动会受到风载荷和不对称刹车力矩的影响.为分析飞机的动力学行为,给出了一种基于非完整非光滑多体系统动力学的建模与数值仿真方法.将飞机视为由机身和前后起落架组成的多体系统,主起落架的轮子为纯滚动,飞机在地面滑跑时考虑前起落架轮子的侧向滑移.飞机在跑道上滑跑的动力学方程由Routh方程导出,用约束稳定化方法抑制约束的漂移.轮与地面间的摩擦模型为库仑摩擦模型并用集值函数描述,以用于判断前轮是否发生侧向滑移.最后通过数值仿真算例分析风载荷和不对称刹车力矩作用下飞机在跑道上滑行的动力学行为.      

考虑摩擦的周边桁架式可展天线展开动力学分析

结合一种周边桁架式可展天线的模型制作,介绍了该类天线的总体布局、展开收纳机理和设计要点;然后提出了切实有效的数学算法,有效地构造了约束方程及其Jacobi矩阵,实时地模拟了结构收纳过程的约束条件,采用广义逆矩阵的方法分析了周边可展桁架运动过程;研究了摩擦这一非线性因素在其展开过程中的作用机理,建立其展开动力学模型,实现了考虑摩擦的周边桁架式可展天线结构展开动力学分析.编制了仿真程序进行算例分析,验证了分析方法的有效性.      

基于约束线图的超导重力梯度敏感结构型综合

基于约束线图对超导重力梯度敏感结构进行型综合,分析轴向分量与交叉分量单独敏感的结构自由度与约束线图,并结合超导重力梯度测量应用提出两分量同时敏感的结构,利用自由度与约束线图确定约束类型,等效为柔性圆柱副,并选择柔性元件进行并联机构形式的布局.在型综合的基础上,通过模态仿真对比,得到两端与柔性球铰联接的刚性杆结构适于在两分量敏感结构中应用的结论.刚度特性仿真分析结果表明两分量敏感结构具有设计运动方向刚度小、寄生运动方向刚度大、运动耦合小的特点,适于在超导重力梯度测量系统中应用.     

一种变体飞行器的动力学建模与动态特性分析

对可变展长、可变后掠角的变体飞行器物理模型进行了简化,基于Kane方法,将变体运动假设为已知的可控输入,利用约束方程来表示机翼的变形运动,选取飞行器位移运动速度以及角速度在机体坐标系下的6个分量作为广义速率,建立了该变体飞行器的六自由度动力学模型.定义了附加力与附加力矩的概念用于描述变体运动对飞行器产生的动力学影响,仿真结果表明,在平稳飞行条件下,相对于变形引起的空气动力的变化,机翼变形产生的附加力和附加力矩都较小.在不同的变形速度下,对变体引起的飞行器纵向运动响应进行了仿真分析,仿真结果表明,变体过程中飞行器的高度、速度以及俯仰角等状态均会发生很大变化.      

CFE方法在导弹发射动力学建模中的应用

在传统的发射动力学微分方程基础上,添加了CFE(Constraint Force Equation)方法中的约束方程和自由度方程,建立了基于CFE方法的导弹发射动力学模型.以倾斜发射导弹系统为研究对象,定义了1个二自由度约束节点来描述导弹与发射架间的约束关系,通过建模计算得到了发射过程的约束力和约束力矩,最后利用软件ADAMS(Automatic Dynamic Analysis of Mechanical Systems)对该算例进行了计算,取得了一致的结果,从而验证了该方法在导弹发射动力学建模中应用的有效性,扩大了该方法的应用范围.     

考虑摩擦的空间抓捕过程碰撞动力学分析与控制

目前针对空间抓捕中的碰撞进行建模的方法中很少考虑摩擦因素。然而,摩擦作为碰撞过程中不可避免的现象,对碰撞动力学响应及整个空间抓捕任务有着不容忽视的影响。文章结合库伦摩擦理论和拉格朗日乘子法,建立了同时兼顾滑动和粘滞状态的摩擦模型,然后将求解到的摩擦力与碰撞力进行矢量合成,作为对碰撞模型的修正。另外,针对碰撞过程对机器人姿态产生的干扰,利用反 馈线性化设计了解耦控制器对基座姿态漂移进行稳定控制。数值仿真采用一个平面2灢DOF空间机器人,利用动量和动量矩守恒定理验证了动力学模型的正确性,同时基座姿态控制结果收敛说明了控制器的有效性。     

新型折叠翼机构设计

为了进一步缩减折叠翼飞行器（如舰载机、无人机、飞行汽车等）的展向尺寸来节省储存、运输空间,设计一种新型的折叠翼机构方案.它由内翼、中翼和外翼3部分组成,采用连杆机构作为翼面折叠机构,能进行两次折叠.简化运动模型,建立运动微分方程,分析折叠翼机构在地面（无气动载荷）情况下的运动特性;利用多刚体动力学软件进行动力学仿真分析;通过在内翼和中翼间的铰链添加扭力弹簧,改善了折叠翼的驱动性能;开展折叠翼地面折叠/展开实验,实验结果与计算机仿真基本吻合.翼面折叠/展开运动平稳,大幅缩减了折叠翼飞行器地面展向占用空间.      

空间全柔性机构位置分析的刚度矩阵法

柔性机构是一种依靠构件元素的弹性变形传输所希望运动的机构.具有集中柔度的全柔性机构是其中的一种类型.由于空间全柔性机构中存在球副,使得目前通用的伪刚体模型法受到限制,为此提出了一种扩展伪刚体模型法.并以6-RSS并联全柔性机构为例对其位置解问题进行了分析:首先利用结构分析中的位移法建立起柔性铰链的刚度模型,同时通过一系列坐标系的建立和转换,建立起机构的变形协调方程、位置闭环方程及静力平衡方程,进而求得机构的位置解.该方法充分考虑了机构中弹性构件的变形,所得结果更接近实际.