用有限元素法验证挤压油膜短轴承解的有效性

在现代航空燃气涡轮发动机和地面高速旋转机中,要求精确地分析预计转子支统的承系振动特性。而在支承结构中带有挤压油膜轴承阻尼器的情况下,用什么方法解算油膜压力,其计算精度和速度如何,是影响整个系统特性分析预计的最主要因素之一,也是这一领域里最重要的研究课题。当前,对不带端封的挤压油膜轴承阻尼器,常采用一维短轴承近似解析法计算油膜压力,本文将通过较精确的二维有限元数值解求油膜压力分布,来验证一维短轴承近似解的有效性,并提出短轴承近似解析法所适用的运转条件。     

弹性飞机平衡的阵风外载荷计算与分析

连续紊流是适航规范中阵风载荷分析必须考虑的一种阵风模型,但其引起的阵风载荷在飞机结构设计中的应用一直存在困难.在平稳随机过程的假设下,基于von Karman连续紊流功率谱模型和线性系统连续紊流响应的功率谱方法,使用载荷累加法获得外载荷形式的阵风载荷.在此基础上,通过引入载荷响应的相关性推导并建立了一种生成多组连续紊流平衡外载荷分布的方法.垂向连续紊流情况下的算例结果表明:连续紊流外载荷分布是全机平衡的.平衡的外载荷分布可通过机动飞行载荷分析中的载荷包线方法识别载荷设计情况,从而应用于飞机结构的强度分析与设计.     

一种解析和数值相结合的机器人逆解算法

针对不存在解析逆运动学解的机器人结构,提出一种解析与数值相结合的求解算法.将传统的算法简化为寻找一个合适的目标函数,对一维关节变量进行迭代求解,其余关节变量可以用解析式求得,从而将逆运动学的多变量迭代求解问题简化为一维变量迭代求解.利用该算法对一种不存在解析逆解的5自由度机器人进行解算,耗时不足0.3ms,计算速度优于传统的算法.这种算法的实时性和准确性满足了机器人实时控制的要求,不需要进行正运动学计算,可以解决大多数关节正交结构的机器人逆运动学问题.      

直接再入大气的月地转移轨道设计

针对从月球停泊轨道出发直接再入大气的月地转移轨道设计问题,提出了一种数值求解算法。该算法由初值设计和精确解求解两部分组成。首先,根据轨道设计的相应约束,采用伪状态理论,通过简单迭代求解高精度的初值。然后,考虑精确的动力学模型,通过数值积分计算真实轨道和状态转移矩阵,并利用微分修正方法搜索精确解。该算法通过设计高精度的初值,降低了月地转移轨道的设计难度。数值仿真表明:该算法求解效率高,具有良好的鲁棒性。      

轴对称磁静平衡态Ⅰ.微变换方法

轴对称磁静平衡态被广泛用来描述不同尺度的恒星大气结构。对这类平衡态,可引入所谓磁通量函数来描述磁场,它满足一变系数、非线性椭圆型方程。迄今为止,人们只就线性情况的某些特例给出该方程的解析解;对较为复杂的情况,特别是非线性情况,还未找到有效的办法进行解析处理。本文将应用微变换方法求上述变系数、非线性方程的相似解,并给出普遍结果。有关特解的具体形式及对太阳大气的应用将在本系列的下几篇论文中另作讨论。      

UKF参数估计在三体Lambert问题中的应用

为了快速精确地求解三体Lambert问题,提出了一种新的基于无损卡尔曼滤波(UKF)参数估计的数值求解算法,该算法由初值猜测和精确解求解两部分组成.首先,基于地月系统二体模型,通过简单迭代求解三体Lambert问题的初值.然后,将三体Lambert问题对应的两点边值问题转化为参数估计问题,通过UKF滤波算法求解,可得到收敛的精确解.该算法是基于概率估计理论的,不仅避免了传统数值方法推导相关梯度矩阵的复杂性,而且降低了三体Lambert问题对初值精确度的要求,从而显著降低了三体Lambert问题求解的难度.数值仿真表明,该方法求解效率较高,具有良好的鲁棒性,与微分修正算法、二阶微分修正算法对比具有更大的收敛域.      

固体运载器飞行动力学逆问题分类

为了降低问题分析及求解难度,对固体运载器出现的飞行动力学逆问题进行了分类研究.通过归纳逆问题解的特点,将其分为3类,并给出了每类问题的统一提法.进一步提出了逆问题的解析求解和优化求解2种基本解法.解析解法通常求解较快,但精度不高.优化求解可以得到较为精确的解,但通常耗时.针对每类动力学逆问题选取有代表性的典型工程问题,即主发动机需用推力设计、合外力辨识和误差分离问题,分别进行了分析并得到了求解结果.典型工程问题的分析和计算结果表明了分类的合理性.     

一种地月转移轨道中途修正的制导算法

针对地月转移轨道中途修正问题,提出了一种求解修正速度增量的制导算法。该算法由初值设计和精确解求解两部分组成。首先,利用伪状态理论,通过简单迭代设计中途修正的初值,并通过Vinti预报方法修正了地球扁率的影响。然后,在求解精确解时,提出了一种基于伪状态理论的状态转移矩阵解析算法。该算法通过设计高精度的初值,降低了求解地月转移轨道中途修正问题的难度,而且避免了传统数值方法计算状态转移矩阵的复杂性。数值仿真结果表明,该算法可有效求解中途修正问题。     

弹性欧拉梁大变形分析的椭圆积分统一形式

一端固支、一端受集中载荷的欧拉梁受载问题是一种基础的力学模型,具有重要的理论研究意义。针对传统的线性求解方法在大变形分析中不适用及无法计算中心定向受压杆在载荷系数超过临界值后屈曲变形的问题,提出一种非线性精确解来进行受集中载荷梁的大变形计算方法。通过椭圆积分形式来推导受集中载荷梁的变形表达,考虑在固支梁自由端加载任意角度下的定向载荷及随动载荷,给出形式统一的梁大变形方程,求解一定载荷因子系数及载荷方位角组合下的变形结果;同时利用此形式对定向受压杆的平衡分支解问题进行了分析。所提方法计算结果准确,可以应用于弹性欧拉梁受定向及随动载荷的大变形分析。      

Von Karman模型三维大气紊流仿真理论与方法

基于三维紊流场相关函数矩阵,用蒙特卡罗法仿真生成Von Karman模型三维空间大气紊流场,实现了对复杂的Von Karman模型的精确推导和应用,仿真流场的相关特性在理论上可以任意逼近Von Karman理论模型.通过提高白噪声序列的白化程度及优化蒙特卡罗算法的方法提高了仿真精度和效率.数值仿真结果证明:所生成的空间大气紊流场符合模型特性,具有较好的统计特性和仿真精度,满足实时飞行仿真对紊流场数据的要求.     

基于积谱的雷达目标高分辨率距离像识别

针对飞机目标的识别问题,将一种用于语音识别中的积谱特征引入雷达目标一维距离像识别领域.积谱定义为功率谱与群延迟的乘积,该特征能够充分利用目标一维距离像的幅度谱与相位谱信息.在目标分类阶段,选择基于弹性传播(RPROP)算法的多层前馈神经网络作为分类器.利用4种飞机模型的重点散射源二维分布测试数据和频率步进法得到目标的一维纵向距离像,对距离像积谱的分类性能进行了测试,结果表明基于积谱的特征对于一维距离像具有较高的识别率,并具有较好的抗噪性能.      

大气紊流数字仿真的改进方法

赵震炎等提出了一种数字摸拟Dryden大气紊流的方法。该方法使用方便,并具有一定的理论基础。但是在某些步长情况下它的相关函数和Dryden模型的理论值吻合得还不够令人满意。本文从理论上分析了产生误差的原因,并且提出了改进的方法。通过相关性检验可以看到改进后的仿真结果非常符合Dryden谱的大气紊流模型。     

正交各向异性矩形薄板稳定问题的直接解法

用分离变量法求解了正交异性矩形薄板的稳定问题,得到了3种边界条件下的解析解,结果与相关文献完全相同,求解方法比传统的半逆法要简单.分析了力边界条件给实施分离变量法带来的困难.给出了不产生泊松效应并能消除板刚体位移的方法.将叠层材料板等效为正交各向异性板,用解析方法求解了面内单向压曲问题,结果与有限元分析结果吻合很好,证明了分离变量方法的正确性和实用性.     

铁木辛柯梁固有振动频率的边界元解法

根据傅里叶变换推导具有两个广义位移的铁木辛柯梁固有振动的基本解.利用加权残量方法,从控制微分方程出发建立边界积分方程,进而根据边界条件得到频率方程,采用代数特征值方法和影响系数方法求解频率,并分析了两种方法的特点.以杆为例证明了对于一维均匀结构,对不同的边界条件利用边界元方法(BEM,Boundary Element Method)都可以得到精确频率.将铁木辛柯梁的BEM结果与有限元结果和精确解进行了比较.     

飞机突防飞行非线性解耦控制系统的混合仿真研究

飞机低空突防作地形跟踪／回避飞行时，会出现纵向和横向运动的非线性交叉耦合，利用非疆性逆系统理论和现代最优控制理论可以设计一个飞行控制规律来解开这些耦合，主要讨论了这种线性解耦控制系统的混合仿真，给了用于仿真的数学模型以及仿真方法，对系统有无大气紊流扰动的解耦控制、系统对指令输入的跟踪、飞行品质等方面进行了探讨。结果表明，该控制规律在有无大气紊流影响时都能实现解耦，准确地跟踪期望的三维突防航迹，并具     

基于改进多目标遗传算法的舰尾紊流模拟方法

为提高舰尾紊流自由大气紊流分量仿真的可信度,提出了一种紊流数值模拟的新方法.首先,使用带有修正系数的Euler前向差分格式表示各个方向的紊流序列,同时结合智能算法的思想,把紊流相关性检验中的均方差误差和相关函数误差作为待优化目标函数,将修正系数的选择看成一个多目标优化问题,并采用改进的多目标遗传算法进行求解.最后,通过仿真算例验证了本文方法的正确性与合理性,计算结果表明该方法可以根据不同的采样步长灵活地生成所需紊流.尤其在小步长情况下,亦可得到很好符合理论值的紊流序列,可以满足虚拟飞行实验的要求.      

有热辐射边界条件的热传导问题

本文研究在有热辐射的非线性边界条件下,在圆柱体内的温度分布问题。本问题可以化为非线性积分方程问题,并最终地利用非线性代数方程组求解,从而得到了温度场的级数形式的解析解。     

微重力环境下蒸发液层热毛细对流的数值模拟

提出了一种新模型来研究由单一物质构成的液层在其纯蒸气中的蒸发.液层置于微重力环境中并且受到水平方向温度梯度的作用,液层的热毛细对流和蒸发耦合在一起,使得气液界面的传热传质规律更加复杂.用理论分析的方法求解了不考虑热毛细效应的纯蒸发模型,得出温度场分布和界面质量流量的解析表达式.对于热毛细对流和蒸发耦合情况,采用有限差分的投影算法同时求解Navier-Stokes方程和能量方程,得到了不同蒸发Biot数和Marangoni数下流场和温度场的稳态数值解.论述了蒸发Biot数和Marangoni数对界面传热传质的影响,提出并解释了蒸发和热毛细对流耦合的三种模式.      

ESF大尺度湍谱分布特性

本文从考虑磁剪切效应的电子连续性方程出发,直接进行两点重正化,建立并解析求解了电子密度涨落关联函数<δn(1)δn(2)>的演化方程.得到的谱<(δn)2>k在谱形,强度变化和最小方差拟合的谱指数等方面均与赤道扩展F大尺度湍谱分布实验结果符合较好.统一解释了通常分为Medium,Intermediate两区的实验结果,同时说明了F峰上、下的谱形相似和垂直与平行于地面的两个方向上的谱形相似问题.     

卫星定位和精度因子的改进方法

在卫星导航定位系统中,在精度因子计算和采用最小二乘法进行定位求解时,传统上采用测量矩阵直接求逆方法来进行.为了克服矩阵求逆带来的计算量大和数值稳定性差的不足,利用测量矩阵的对称正定性,提出了一种基于矩阵 UTDU 分解的定位解算和精度因子计算方法.改进方法具有严格的数学理论基础,保证了方法的正确性和有效性.数值分析结果表明,相对直接求逆的传统方法而言,在定位解算时,该方法能降低约 60%的运算量,而在精度因子计算中,约能降低36%的运算量.且改进方法能大大降低求解矩阵的条件数,提高了求解的数值稳定性.