求解非线性方程及方程组的模拟退火算法

探讨了用模拟退火算法求解非线性方程及方程组的方法,该方法把方程组求解问题转化为函数优化问题.计算中不需要使用目标函数的导数信息,在Matlab环境下实现了该算法,数值实验结果表明了该算法的有效性.     

一类非线性系统的混沌控制

 描述航天器、陀螺和气浮台等刚体姿态运动的欧拉动力学方程,是一个具有广泛代表意义的三阶非线性方程。当该方程中的参数取不同值时,可得到著名的Lorenz系统、Rssler系统、Newton-Leipnik系统、Chen系统及Lü系统。在不同的外力矩作用下,该动力学系统会呈现出相当复杂的动力学行为。从该系统中,发现了一大类新的混沌吸引子。本文分析了这一类混沌吸引子具有的共同特征,并采用基于输出反馈的PI型控制器将一种新的混沌运动稳定于指定平衡点。仿真结果表明,该控制器能够有效地抑制混沌,能将系统稳定于任意指定的不稳定平衡点。     

液氢贮箱停放过程中的力热分析

为了寻找到运载火箭长时间停放过程中液氢贮箱的最经济液位,用于制定合理的发射流程以及紧急处置方法,采用计算流体力学(CFD)技术,对某型运载火箭停放期间液氢贮箱的力、热情况进行了仿真计算和分析。计算选用了VOF(Volume-of-fluid)两相流模型以及Lee相变模型,为了提高Lee模型在不同压力情况下对相变过程的模拟精度,采用安托因方程修正了该模型。修正后的模型首先由试验数据校验了其精确性,随后开展的液氢贮箱停放过程仿真结果表明：贮箱的竖直方向与径向均存在温度分层的现象,液相内会形成大的漩涡,该漩涡会使得冷热流体不断进行热交换,并导致贮箱内部的液氢出现气化。贮箱停放期间蒸发率最大值超过2 m 3/h,发生在停放4 h左右;而贮箱液位充填至37 m 3以上或17 m 3以下时蒸发率较低,最小值接近1 m 3/h。     

An improved numerical method for constructing Halo/Lissajous orbits in a full solar system model

An improved numerical method that can construct Halo/Lissajous orbits in the vicinity of collinear libration points in a full solar system model is investigated. A full solar system gravitational model in the geocentric rotating coordinate system with a clear presentation of the angular velocity relative to the inertial coordinate system is proposed. An alternative way to determine patch points in the multiple shooting method is provided based on a dynamical analysis with Poincare′sections. By employing the new patch points and sequential quadratic programming, Halo orbits for L_1, L_2, and L_3 points as well as Lissajous orbits for L_1 and L_2 points in the EarthMoon system are generated with the proposed full solar system gravitational model to verify the effectiveness of the proposed method.     

单推力航天器交会对接轨迹规划及跟踪控制

针对单推力航天器交会对接问题,提出一种轨迹规划及跟踪算法。首先,考虑到追踪航天器只沿本体X轴安装推力器,且推力方向固定,为了实现从起始位置转移至期望位置并满足姿态要求,基于三维螺旋线设计两阶段转移轨迹,根据初末位置以及末端速度方向要求,求解螺旋线参数。该螺旋线可以保证在初末速度方向固定情况下,曲率积分最小。其次,为了降低轨迹跟踪难度并减小初始时刻的位置跟踪控制力,需要将转移轨迹初始速度与追踪星X轴重合。传统螺旋线无法满足该约束条件。本文对传统螺旋线进行改进,提出一种旋转螺旋线轨迹设计方法。通过引入姿态旋转矩阵,将螺旋线在三维空间旋转,在不改变曲线形状的前提下满足初末位置及速度方向要求。然后,为了跟踪转移轨迹以及跟踪期望推力方向,提出基于CLF（Control Lyapunov Function）的滑模控制策略,当追踪星X轴与期望推力方向夹角较大时,采用CLF,保证最优性;当姿态误差收敛至滑模面附近时,切换为滑模控制,以提升系统鲁棒性。最后,通过仿真验证旋转螺旋线相比于传统螺旋线的优势。     

多芯片组件（MCM）用金锡焊带材料电阻率及热导率温度特性研究

对牌号为Au₈₀Sn₂₀的金锡焊带材料在208~423 K的电阻率及热导率与温度的函数关系进行了研究,并对其在多芯片组件（MCM）中的传热效果进行了评估。分别对材料在208~423 K中5个温度点的电阻率及4个温度点的热导率进行了测试,基于理论模型建立电阻率/热导率随温度变化的函数关系,最终采用模拟热扩散数值方法评估材料在高温下的传热能力。结果表明,采用修正函数模型后,金锡焊带材料在208~423 K下热导率与电阻率的关系符合测试结果,随芯片表面温度的边界条件从208 K升高至423 K,采用变温热导率模型得到的热流密度模拟计算结果相比理想化恒定热导率模型的差异性逐渐升高至5.5%。综上,金锡焊带材料热导率与电阻率的关系符合Wiedemann-Franz法则修正后的Smith-Palmer方程,在该材料传热设计时应考虑其热导率温度效应。     

某机第4级压气机转子叶片榫头折断故障分析

阐明了某机第4级压气机转子叶片在榫头部位断裂的故障特点。故障分析表明：经长期使用的叶片在修理厂大修时，为满足叶片摆动量要求在榫头底面涂尼龙胶，会使叶片自振频率有较大幅度下降；为排除3级叶片折断故障而提高的慢车转速又缩小了4级叶片的共振裕度，结果造成少数4级叶片落入慢车共振。共振发生时，不规则的尼龙胶和装配过程不规范引起的非正常接触容易导致微动磨损的加剧，从而大大降低了疲劳寿命，这是引起榫头断裂故障的主要原因。     

涡喷发动机加速控制系统数学模型及数字仿真

 本文介绍了航空涡轮喷气发动机加速控制系统数学模型建立、数字仿真方法。由于加速过程是一个大偏离动态过程,因而发动机与加速控制器模型均用非线性方程描述,对非线性方程组求解和系统迭代方法本文进行了讨论。算例中分析了在不同飞行条件下发动机的加速规律。     

叶轮机械三维不可压Euler方程的直接解法

本文用近似因式分解交替方向隐式格式直接求解叶轮机械中原参数三维不可压Euler方程组。在连续方程中引入“拟压缩性”之后,将原来类型不确定的控制方程变成双曲型方程,给定初边值之后用时间推进法求得稳定解,对一个单级压气机转子内部流场进行了计算,结果与实验值比较符合。     

战斗机弹射救生系统数值模拟研究

对弹射救生系统这种几何外形复杂、不规则非流线型钝体,以Jameson有限体积法为基础发展了一套基于高质量混合网格的Euler方程解算器,在Euler方程解算器基础上,采用S-A湍流模型发展了一套同样基于混合网格的N-S方程解算器,初步建立了适合弹射救生系统这类几何外形复杂物体的数值模拟技术平台,并采用多种加速收敛措施和合适的耗散项模型使得所发展的解算器在软件实现上鲁棒性很强。通过选取合适的计算起始迎角,对弹射救生系统成功地进行了大迎角和大侧滑角绕流计算,获得了与风洞试验比较吻合的计算结果,并且可以清楚地观察流场中马赫数分布、弹射救生系统各部位压力分布等流场细节,对弹射救生系统空气动力特性分析具有重要的指导意义。该计算平台可作为风洞试验的合理补充、加强和完善。