伞翼无人机精确建模与控制

前缘切口以及后缘下偏是影响伞衣气动力计算的关键因素。为实现伞翼无人机（UAV）的精确控制，从提高翼伞系统动力学模型的精度入手，在升力线理论的基础上，基于计算流体动力学方法，综合考虑前缘切口以及后缘下偏的影响，计算了不同切口尺寸模型的升力、阻力系数。利用最小二乘法辨识了升力、阻力系数与迎角、切口尺寸以及下偏量的关系，实现了翼伞气动力的精确计算，改进了伞翼无人机的六自由度动力学模型。对改进的动力学模型进行轨迹跟踪控制的仿真，通过与空投试验数据的对比，验证了改进翼伞系统动力学模型方法的准确性，对于伞翼无人机的仿真和控制器设计具有重要意义。     

关于Fleury算法的一点注记

为了使用Ｆｌｅｕｒｙ的算法，在每一步都必须去判断图Ｇ－ｅ的连通性［１］。本文将给出一个十分简单的判断图的连通性的线性算法。为了证明它的正确性，本文将证明以下三个条件是等价的：（１）图Ｇ是连通的；（２）Ｍ的任意ｎ－１行都是线性无关的（这里Ｍ为图Ｇ的关联矩阵）；（３）在Ｍ中存在ｎ－１个线性无关的行。     

高频不稳定燃烧的声学数值仿真

高频不稳定燃烧一直是液体火箭发动机研制过程中所要面临的重大难题之一.采用具有低色散低耗散特点的计算气动声学方法,对自燃推进剂变轨发动机(OME)的高频不稳定燃烧进行时域下的数值仿真.由Crocco的压力时滞模型对燃烧热释放和声波之间的耦合进行模拟,并对不同的时滞模型参数对稳定性结果的影响进行了分析,给出发动机的稳定性极限图,确定一阶切向及一阶径向振型为主的不稳定振型,与地面试车实验捕捉到的不稳定振型相一致.结果表明:采用计算气动声学方法对带有Crocco压力时滞模型的声波扰动方程进行时域下的数值求解,可以对发动机的高频不稳定燃烧进行成功地预测.      

CFD数学模型的线性化方法及其应用

计算流体力学(CFD)方法不仅仅起到数值模拟的作用,它本身是一个复杂的非线性系统。在流动稳定性分析、气动弹性分析、优化设计以及流动控制等领域,从系统的角度出发,对CFD数学模型线性化后,可以对模型的系统矩阵进行定量分析获得更多的系统特性。但是CFD数学模型往往非常复杂且阶数很高,因此其线性化系统矩阵的获得比较困难。鉴于此,采用人工编程和自动微分相结合,构造有限体积法并行CFD模型的线性化系统矩阵。其中自动微分只被用来得到每个界面通量的局部雅可比矩阵,而采用人工编程方法来实现并行环境下的稀疏雅可比矩阵的组装。线性化系统的并行求解采用了块雅可比预处理的广义最小残量法,每个并行进程内部则采用零填充不完全LU分解预处理。为了验证这种线性化方法,上述方法被用于:① NACA 0012翼型的非定常绕流线性系统构造与求解;② NACA 0012翼型稳态流动的伴随方程构造与求解;③ AGARD wing 445.6机翼颤振问题降阶建模。上述三个算例的结果与CFD模拟的吻合一致。     

高温燃气进入舵机舱过程仿真与流动机理分析

采用计算流体力学(CFD)方法研究了火箭发动机工作拖尾段高温发动机燃气进入舵机舱的物理现象。结合导弹实际飞行弹道参数变化特点和超声速流场扰动不向前传递的空气动力学理论,提出了简化而不失真的非定常流场仿真方案,显著缩短了仿真周期;复现了某型导弹实际飞行时舵机舱先被"抽气"再进高温燃气的动态过程,并分析了高温发动机燃气进入舵机舱的流动机理,即在发动机工作段,导弹底端面压强低于舵机舱内压强,舵机舱被"抽气",在拖尾段随着燃烧室总压降低,喷口附近的马赫盘向导弹底端面移动,使导弹底端面压强增大且高于舵机舱内压强,高温燃气进入舵机舱烧毁电路致使导弹折断;明确了某型导弹折断故障产生的诱因,提出了改进措施和检测方法,并得到了大量飞行靶试的验证,解决了舵机舱热防护结构可靠性问题。     

旋翼VPM/CFD耦合模型的建立及其在小前进比气动分析中的应用

为了提高旋翼计算流体力学（CFD）流场计算效率及克服其在尾迹捕捉上的不足,将旋翼黏性涡粒子方法（VPM）与CFD分析相结合,建立了一个新的旋翼VPM/CFD耦合气动分析模型。在该模型中,采用VPM分析以实现对旋翼尾迹中黏性涡的高效捕捉而不引入数值耗散,而采用CFD分析用于精确地模拟旋翼桨叶近体区域内复杂的流动现象,同时也为VPM分析提供一个较高精度的涡源模型。至于两者信息交换,则使用集中涡源法将CFD信息传递至VPM分析,而VPM计算得到的尾迹信息则通过边界修正施加至CFD域的远场边界上,从而可以鲁棒地实现CFD域与VPM域的耦合计算。在此基础上,对“Helishape 7A旋翼”小前进比前飞桨-涡干扰（BVI）状态进行了较为深入的气动分析,计算结果表明：与全CFD计算比较,建立的VPM/CFD耦合分析模型可以有效地避免旋翼尾迹区桨尖涡的数值耗散,从而更加可靠地捕捉桨-涡干扰状态下的桨叶非定常气动载荷脉动,同时对于本文算例,计算效率可以提高30%以上。     

防范秘密攻击的安全计算的博弈论实现

在博弈论中,惩罚博弈模拟了参与者试图欺骗但又不想被抓住,即安全计算中秘密攻击者的情形。针对密码学的计算博弈模型,本文对Halpern与Rafael提出的能否在计算具有成本的惩罚博弈与具有一定威慑度的防范秘密攻击的安全计算之间建立联系的问题给出肯定的回答,提出威慑度为1/2的防范秘密攻击的安全是计算博弈中错误可忽略的调解人的通用实现。     

火星大气环境模拟装置设计及仿真分析研究

对火星表面大气环境特性进行了研究,通过选取合适的计算方法并结合FLUENT流体有限元计算软件对火星表面稀薄气体内部环流进行了模拟仿真分析,提出了以动量源模拟风扇段内流的仿真方法,并进行了可行性讨论。进一步实现了针对圆柱形模拟装置多工况下的内部气体流场稳态和非稳态计算仿真,并对计算结果进行了分析讨论,为火星大气环境模拟装置的设计提供了技术支持和参考。     

严格随机正则(3,s)-SAT模型及其相变现象

研究变元和文字出现次数受限制的规则3-SAT问题,提出了一种严格随机正则（3,s）-SAT问题,并给出了该问题的实例产生模型--SRR模型。结合一阶矩方法和生成函数展开项系数的渐近近似技术,证明了严格随机正则（3,s）-SAT问题相变点的上界,即当变元规模N较大且变元出现次数s>11时,严格随机正则（3,s）-SAT实例是高概率不可满足的。实验结果表明：由SRR模型所生成的随机实例中,当N>60且s>11时,所有的（3,s）-SAT实例均是不可满足的,而当N>150且s<11时,所有的（3,s）-SAT实例均是可满足的,即严格随机正则（3,s）-SAT实例的相变点位于s=11处,且在s=11处（子句变元比为11/3）的严格随机正则（3,s）-SAT实例,比在相变点（子句变元比）4.267处同规模的均匀随机3-SAT实例更难求解,因此,SRR模型可以很方便地在s=11处构造难解的随机3-SAT实例。     

密切曲面锥导乘波体的设计与理论分析

将密切锥导乘波体技术应用于一般曲面锥流场,以期获得升阻比、容积率更高的密切曲面锥导乘波体。首先,设定不同乘波体后缘激波型线,在两种代表性曲面锥流场内,生成了四种乘波体构型,利用数值方法验证乘波体设计方法可行性;然后构造了三类典型单一控制变量的曲面锥流场,对比分析了对应密切曲面锥导乘波体性能变化规律。研究表明：（1）密切曲面锥导乘波体流面压力分布、后缘激波型线与理论设计吻合,关键位置压力与理论值相对偏差约1%,说明利用密切技术在曲面锥流场中“截取”乘波体的方法是可行的。（2）曲面锥流场控制参数不同,可获得升阻系数、容积率、压缩量增加,升阻比降低的乘波体,也可获得升阻比、容积率增加,升阻系数降低的乘波体。