行星运动轨道方程的又一推导方法

 提出了推导行星运动轨道方程的一种新方法------ 矢量分析法, 较先推导比耐公式, 再求解微分方程的传统方法简单, 所得结果与传 统方法一致.     

用复矢量法研究平面机构的运动

用复数－极坐标矢量法［１］（简称复矢量法）,证明了牵连运动为平面运动时,点的加速度合成定理,给出了求解平面机构运动的新方法．     

多体系统指标2运动方程HHT方法违约校正1）

采用Cartesian绝对坐标建模方法,完整约束多体系统运动方程是指标3的微分--代数方程（differentialalgebraic equations,DAEs）,数值求解指标3的DAEs属于高指标问题,通过对位置约束方程求导,可使运动方程的指标降为2.位置约束方程求导得到的是速度约束方程.直接求解指标3的运动方程,速度约束方程得不到满足,而且高指标DAEs的数值求解存在一些问题.论文首先采用HHT（Hilber--Hughes--Taylor）直接积分方法求解降指标得到的指标2运动方程,此时速度约束方程参与离散计算,从机器精度上讲速度约束自然得到满足,而位置约束方程没有参与计算,存在“违约”.针对违约问题,采用基于Moore--Penrose广义逆理论的违约校正方法,消除位置约束方程的违约.指标2运动方程HHT方法违约校正,将HHT方法和违约校正方法很好地结合,在数值求解指标2运动方程的过程中,位置约束方程和速度约束方程都不存在违约问题,而且新方法没有引入新的未知数向量,离散得到的非线性方程组的方程数量与原指标2运动方程的方程数量相同,求解规模没有扩大.新方法的实用和有效性通过算例的数值实验得到验证,数值实验也说明新方法保持了HHT方法本身具有的数值阻尼可以控制和二阶精度的特性.最后从非线性方程组的求解规模和计算速度上与其他方法进行了比较分析,说明新方法的优势所在.     

多体系统指标2运动方程HHT方法违约校正

采用Cartesian绝对坐标建模方法,完整约束多体系统运动方程是指标3的微分-代数方程(differentialalgebraic equations,DAEs),数值求解指标3的DAEs属于高指标问题,通过对位置约束方程求导,可使运动方程的指标降为2.位置约束方程求导得到的是速度约束方程.直接求解指标3的运动方程,速度约束方程得不到满足,而且高指标DAEs的数值求解存在一些问题.论文首先采用HHT(Hilber--Hughes--Taylor)直接积分方法求解降指标得到的指标2运动方程,此时速度约束方程参与离散计算,从机器精度上讲速度约束自然得到满足,而位置约束方程没有参与计算,存在"违约".针对违约问题,采用基于Moore--Penrose广义逆理论的违约校正方法,消除位置约束方程的违约.指标2运动方程HHT方法违约校正,将HHT方法和违约校正方法很好地结合,在数值求解指标2运动方程的过程中,位置约束方程和速度约束方程都不存在违约问题,而且新方法没有引入新的未知数向量,离散得到的非线性方程组的方程数量与原指标2运动方程的方程数量相同,求解规模没有扩大.新方法的实用和有效性通过算例的数值实验得到验证,数值实验也说明新方法保持了HHT方法本身具有的数值阻尼可以控制和二阶精度的特性.最后从非线性方程组的求解规模和计算速度上与其他方法进行了比较分析,说明新方法的优势所在.     

基于参数识别视速度的双矢量粗对准方法（英文）

针对传统基于视速度双矢量粗对准中,由于传感器随机噪声的影响,存在对准精度差,收敛速度慢的缺点,提出了一种新型自适应Kalman滤波的参数识别粗对准方法。该方法通过对视速度运动进行建模,设计采用自适应Kalman滤波对模型参数进行参数识别,从而有效地消除视运动中的随机噪声,提高粗对准的精度和收敛速度。由于自适应滤波的特点,新方法不需要对传感器误差进行统计,使其在实际系统中具有更加广泛的应用价值。针对双矢量粗对准的计算特点,设计了一种矢量重构算法,从而尽可能地规避双矢量共线性问题,加快了粗对准的收敛过程。仿真与转台实验表明,与传统方法对比,新方法在相同的对准时间内具有更高的对准精度,在相同的对准精度下,具有更高的收敛速度。转台实验的最终对准精度为-0.1391°,标准差为0.012°。     

散斑照相孔径优化的研究

本文提出了三种新的散斑照相孔径方法:一维调频孔径法,交叉孔径法和旋转孔径法。给出了理论分析和实验证明。这些方法简单实用,能有效地改变衍射晕的能量分布。适当选择孔径的参数,实验者能相当任意地控制散斑图的空间频率成份。新方法兼顾了突出某一频率成份、灵敏度连续可调、获得多个方向散斑场运动信息等各方面的特性。     

等效旋转矢量系数优化的参数解析法

基于等效旋转矢量的姿态算法是减小捷联姿态计算中不可交换性误差的有效方法.针对Miller的多子样优化方法过程复杂且仅适用于优化圆锥运动的局限,提出了通过角速度矢量的高一阶近似模型并利用其各阶导数间的解析关系,实现系数优化的新方法--参数解析法.以Miller的三子样基本算法为例阐述了参数解析法的优化思路并详细推导了其优化过程,在经典圆锥运动环境下得到了与Miller法相同的优化结果.最后在规则进动和广义圆锥运动环境下对该方法的适用性进一步进行验证.该方法简单直观,无需求解误差四元数的解析式,适用范围广,为等效旋转矢量的系数优化提供了一种新思路.     

上覆岩层组合运动的力学解析特征研究

运用弹性薄板理论和复合材料力学的层合板理论,讨论了在覆岩运动中起决定作用的组合关键层的应力,为计算组合关键层的整体强度提供了途径;并且从控制开采沉陷工程实际入手,将结构塑性极限分析理论运用于组合关键层的失稳最大荷载计算,解析得出了充分采动情况下关键层塑性极限荷载的表达式,从另一个角度得到了一种判别关键层位置的新方法,为实施离层带注浆减沉提供参考。     

一种确定非线性裂纹转子解的形式的新方法

将小波变换与Ｐｏｉｎｃａｒｅ映射相结合，即用Ｐｏｉｎｃａｒｅ映射确定周期解，用谐波小波变换区分拟周期响应和混沌运动，提出了一种分析非线性裂纹转子系统解的形式随参数变化的新方法．结果表明这种方法是非常有效的，它比以前所用的计算Ｌｉａｐｕｎｏｖ指数的方法节约了计算时间，并且较易实施．     

基于群桩振动方程组的群桩基础阻抗函数求解及比较

本文首先在位移场叠加的基础上,对经典的单桩振动方程进行了修正得到群桩振动方程组;然后根据群桩振动方程组的近似求解,导出了群桩阻抗函数求解的常用方法一群桩阻抗方程方法。本文中给出的群桩阻抗方程方法在具体的表达形式上与常见的阻抗方程计算方法略有差别;而后根据振动方程组系数矩阵在复空间为正规阵的性质,提出了该四阶耦合微分方程组解耦求解的详细步骤,在直接求解微分方程组的基础上给出了群桩基础阻抗计算的一种新方法。与阻抗方程方法相比较,新方法未采用柔度因子来考虑桩与桩之间的相互影响,而是通过耦合的运动方程组对桩与桩之间的相互作用作了新的表述,并且近似地考虑了屏蔽效应。本文最后给出了上述两种计算方法的数值算例。     

聚龙一号上磁驱动铝飞片发射实验的数值分析与再设计

聚龙一号上PTS-151发次实验中,磁驱动加速370 μm厚飞片测得的最大速度为18 km/s,磁驱动加速482 μm厚飞片测得的最大速度为19 km/s。采用MDSC2程序, 对PTS-151发次实验进行了数值分析,结果表明:PTS-151发次实验中测量的最大速度的含义不同于以往文献中飞片的最大速度。以往文献中发射飞片在测试过程中自由面未被烧蚀,测试的最大速度为飞片自由面速度;PTS-151发次实验中两个飞片在测量过程中自由面被烧蚀,实验测量的最大速度为飞片被完全烧蚀前的一瞬间飞片内部最后一个固体面的速度。在飞片自由面未被烧蚀之前,370 μm厚飞片的计算最大自由面速度仅为7 km/s,482 μm厚飞片的计算最大自由面速度仅为11.8 km/s,远低于测量值。对PTS-151发次实验条件下飞片尺寸进行了再设计,飞片厚度为680 μm时最优,既能保证自由面未烧蚀,又使得飞片的速度最大,达到17.5 km/s。     

Non-planar,non-linear oscillations of a beam—I. Forced motions

Large amplitude whirling motions of a simply supported beam constrained to have a fixed length are investigated. Equations of motion taking into account bending in two planes and longitudinal deformations are developed. Using the method of harmonic balance, response curves for certain planar and non-planar steady state, forced motions are obtained. Another approximate scheme is used to study the stability of these motions. Stable regions corresponding to non-planar motions are found, thus confirming the existence of whirling motions. Numerical results are presented and discussed for several specific cases.     

两柔性梁碰撞振动类型的实验研究

基于对两弹性部件碰撞振动的理论及数值研究,设计了两柔性梁碰撞实验来研究碰撞振动类型．在实验中观察到亚谐、超谐、概周期、混沌等多种碰撞振动类型．此外,发现系统存在类似“频率锁定”性质的碰撞振动,以及同时存在性质相异的碰撞行为,如谐碰撞与混沌振动共存现象     

Stability analysis of steady rotations of a rigid body bearing two-degree-of-freedom control moment gyros with dissipation in gimbal suspension axes

To study the stability of steady rotations of a control moment gyro system with internal dissipation, we use the Barbashin-Krasovskii theorem and the relation, established in [1], between the Lyapunov function and steady motions. Taking into account the special properties of the original problem, we reduce it to a lower-dimensional problem.We give a detailed presentation of an algorithm for analyzing the stability of steady motions of a gyrostat and use this algorithm to perform a complete study for two systems consisting, respectively, of one and two gyros whose gimbal axes are parallel to the principal axis of inertia of the system. Each steady motion is identified as either asymptotically stable or unstable. We find periodic motions that exist only in the presence of dynamic symmetry and which are regular precessions. For the system with two gyros, we prove the asymptotic stability of quiescent states and prove that in the angular momentum range where these states are defined the system does not have any other stable motions.     

Free vibrations of a thin elastica by normal modes

In this work we investigate the existence, stability and bifurcation of periodic motions in an unforced conservative two degree of freedom system. The system models the nonlinear vibrations of an elastic rod which can undergo both torsional and bending modes. Using a variety of perturbation techniques in conjunction with the computer algebra system MACSYMA, we obtain approximate expressions for a diversity of periodic motions, including nonlinear normal modes, elliptic orbits and non-local modes. The latter motions, which involve both bending and torsional motions in a 2:1 ratio, correspond to behavior previously observed in experiments by Cusumano.     

高机动飞行器非指令运动及其控制的研究进展

高机动飞行器往往都是通过大攻角飞行来实现高机动科目的, 在发展高机动飞行器的过程中, 其非指令运动是伴随着大攻角飞行而常常出现的运动形态. 为此, 应在飞行器设计的早期阶段, 充分研究所设计布局的大攻角流动性态及其相应的非指令运动的形态;揭示这类运动形态的主控流动;在此基础上形成和发展流动控制新技术, 以达到抑制非指令运动的目的. 由于大攻角前体非对称涡往往与非指令运动密切相关, 为此本文首先指出前体非对称涡流动对头部微扰动十分敏感, 以致长期以来让人们误认为这类流动具有不确定性. 研究表明, 通过设置人工微扰动可使前体非对称涡流动具有可重复性, 并揭示该流动随扰动周向角变化的响应、演化规律. 通过利用大、小后掠翼两类翼身组合体的典型布局形式, 研究它们所呈现的摇滚运动形态, 揭示其摇滚运动的不同主控流动机理, 在此基础上分别发展了抑制小、大后掠翼身组合体摇滚运动的流动控制技术: 快速旋转头部扰动和适当设置扰动位使翼、身的两对非对称涡处于反相. 在抑制非指令运动的研究中, 深入理解和揭示头部微扰动对非对称涡流动的响应、演化机理是至关重要的, 应予以特别关注.     

枪械射击残留物分布密度的分形研究

为了推测发射枪械类别、弹种、枪械射击距离,利用X光摄影记录六种枪械发射子弹射击残留物的金属微粒密度与射击距离之间的关系,关系曲线表明二者是非线性的且具有分形特征;应用分形理论和线性回归原理,对射击残留物中的金属微粒分布密度进行计算和分析。结果表明,同种子弹由不同种枪发射时,其分形维数相应不同但相差较少;而不同种子弹由不同种枪发射时,其分形维数不同且相差较大。因此分形维数可以用来描述枪械发射子弹时射击残留物中的金属微粒密度分布情况。     

航空薄壁弧形结构件喷丸强化变形预测

喷丸强化是一种有效提高工件抗疲劳性能的表面处理工艺，经喷丸强化处理后的航空薄壁弧形结构件发生变形是航空制造领域的共性难题．本文基于变截面曲梁平面弯曲理论，建立了弧形件喷丸诱导应力与变形方程，将得到的构件最大变形与有限元方法计算结果进行对比，所建立理论模型能有效预测弧形件喷丸强化后的变形．讨论了喷丸压力和喷丸区域对喷丸变形的影响，结果表明：喷丸强化工艺仅作用于弧形件腹板表面时，弧形件变形以径向收缩为主；喷丸强化工艺仅作用于弧形件筋板表面时，弧形件变形以径向扩张为主；与局部喷丸相比，对弧形件所有表面喷丸强化时，弧形件弯曲变形减小，轴向伸长增大．     

一类双自由度碰振系统运动分析

基于Poincare映射方法对一类两自由度碰撞系统进行了分析。经过详细的理论演算得到单碰周期n的次谐运动的存在性判据和稳定性条件,给出计算Jacobi矩阵特征值的公式。数值模拟表明,该方法具有令人满意的结果。此外,还讨论了当不满足所提出的单碰周期n次谐运动的存在性条件时,可能会出现的运动形式。