用双流体—轨道模型模拟四角喷燃模型炉内三维湍流两相流动和煤粉燃烧

目前炉内两相流动和煤粉燃烧数值模拟中多半用颗粒随机轨道模型和单流体无滑移模型，这些模型都难以完整地给出三维空间内颗粒速度，浓度，湍流度分布的信息。主采用双流体－轨道模型（颗粒相连续介质－轨道模型）对一个四角喷燃模型炉内三维湍流两相流动及煤粉燃烧进行了模拟。此模型基于欧拉气相方程组、欧拉颗粒连续方程组和动量方程组以及拉氏颗粒能量方程和质量变化的方程，并使用k-ε-kp两相湍流模型，EBU－Arrhenius湍流燃烧模型，离散坐标辐射传热模型，煤粉颗粒的水分蒸发，热解挥发模型和焦炭燃烧的扩散－动力模型等。热态模拟中，为了减小为信散造成的影响，采用了扭转坐标法（将坐标扭转一定的角度使之与煤粉射流方程一致）。为了检验数值模拟，采用三维相位移普勒测速仪（PDPA）对于冷态模型炉内湍流两相流场进行了测量，得到了两相速度，湍流脉动及颗粒浓度的分布。分别对冷态模炉内两相流动和热态模型炉内三维两相流动和煤粉进行了模拟，冷态两相流动的计算与实验结果的对比表明预报的两相流场是合理的，热态模拟的结果给两相速度，气相温度、组分浓度及壁面热，显示出靠近出口处气相速度和温度分布不对称，造成一个局部高温区。     

轮盘概念设计中拓扑和形状同时优化方法

为了减小人为选择参数对优化结果的影响,针对航空发动机轮盘概念设计阶段的结构优化问题,提出了拓扑和形状同时优化(STSO)法.该方法是在变密度(SIMP)法的基础上,通过分析优化目标和约束的灵敏度,用序列二次规划优化(SQP)法进行求解.接着,以板壳结构为例,对比分析了STSO法和分步优化法的结果,说明了拓扑和形状同时优化方法的优点.最后,将同时优化方法应用于轮盘结构概念设计,对比分析了使用拓扑和形状同时优化方法与单独拓扑优化方法进行轮盘结构优化的结果,探讨了不同振型对应的频率约束对优化结果的影响.结果表明,不同振型对应的频率约束下优化结果的结构形式呈现多样性;相较于单独拓扑优化方法,STSO法收敛速度较快、结果也更精确;但是由于形状优化变量取值范围选取不当,有可能会出现网格畸变过大,而导致STSO法所得的结果无效.      

基于约束Delaunay三角剖分的筋特征识别与构建算法

为识别飞机结构件中的筋特征,提出一种基于Delaunay三角剖分的识别与构建算法。首先,根据飞机结构件腹板和平顶筋加工方式的不同以及二者几何上的相似性,引入广义腹板概念表示腹板面和平顶筋,并建立筋的表示模型;其次,利用约束Delaunay三角剖分算法,剖分广义腹板面,进行平顶筋面和腹板面的识别与区分,并提取平顶筋面的中轴线;最后,利用中轴线拆分平顶筋并将斜顶筋作为广义腹板的子特征进行识别与构建。实例结果表明,该算法是正确和有效的。      

二元机翼颤振的指令滤波反推自适应约束控制

为有效抑制二元机翼颤振现象,设计了一个指令滤波反推自适应约束控制器。考虑二元机翼气动弹性系统存在多项式结构非线性特性以及系统不确定的情况,针对二元机翼控制面偏转角度限制以及系统外部阵风扰动等问题,将前/后缘双控制面布局的二元机翼动态方程以状态空间形式描述,采用反推控制技术进行控制律整体设计,通过指令滤波环节对虚拟/实际控制信号进行幅值限制,并将滤波前后的信号差引入到自适应约束控制律设计过程中。仿真结果表明,在一定来流速度下,开环系统出现极限环颤振现象,闭环控制系统能快速达到稳定状态,二元机翼颤振现象得到有效抑制。     

基于MPSC和CPN制导方法的协同制导律

针对带有末端攻击角度约束的多导弹协同制导问题,运用模型预测扩展控制（MPSC）和协同比例制导（CPN）,设计了一种满足末端攻击角度约束的多导弹协同次优制导律。阐述了MPSC制导方法的基本理论,详细给出了控制量表达式以二次形式近似时MPSC制导律的设计过程。采用CPN对MPSC制导方法的初始控制量进行猜测,并确定协同攻击时间。仿真时考虑两枚导弹对地面静止目标进行协同攻击。仿真结果表明,两枚导弹攻击时间偏差和末端攻击角度偏差均可控制在给定范围内,即本文所设计的制导律在实现多导弹协同攻击时,还可以很好地满足末端攻击角度约束。      

附加先验轨道约束的LEO几何法实时精密定轨验证分析

低轨卫星的实时精密定轨能够极大拓展其完成复杂科学任务的能力,例如实时环境监测、机动控制和卫星自主导航等.本文根据几何法实时精密定轨模型,提出了附加LEO先验轨道约束从而改善实时定轨的精度、收敛速度和稳定性的构想.分别采用广播星历、超快速星历预报部分和实时精密星历,设计了6种实时定轨方案,并利用Swarm-A,B,C星7天的观测数据进行方案验证与分析.结果表明,使用广播星历、IGU和IGC星历的方案精度递增,附加先验轨道约束能够进一步提升精度.使用IGC星历并附加标准差为1m的先验轨道约束后,在径向、切向和法向的定轨精度分别达到6.12cm,5.55cm和4.98cm.此外,附加先验轨道约束能够显著提升收敛速度,使用IGC星历平均收敛时间约为31min,附加标准差为1m的先验轨道约束后收敛仅需约4min.      

一种适合于工程应用的准C~1连续Bézier三角曲面及其构造方法

本文比较系统地讨论了Clough-Tocher三角形分割模型在Bézier三角曲面设计中的作用。并从工程应用的角度提出了一种以Clough-Tocher分割为基础,旨在原三角形区域上构造9参数Bézier三角曲面的新方法。这种新的9参数插值不仅能够消除Clough-Tother分割产生的畸变插值区域对插值曲面品质的影响,还能够有效地减轻一般9参数三次Bézier三角曲面片之间的尖端连接。本文将从工程应用的角度把这种插值曲面称为准C~1连续的9参数Bézier三角曲面。在一些实际应用中,这种准C~1连续的插值曲面对3D离散数据的拟合效果是令人满意的。     

大攻角多分量强迫振动实验技术

北京大学低速风洞新研制的大攻角强迫振动实验设备可迫使模型绕其体轴做单自由度俯仰、偏航、或滚转振动,由六分量应变天平感受气动反作用力和力矩,测量仪器按相关滤波原理将它们分解为同相分量和正交分量。仪器具有很高的分辨度和精度,能测出组合形式的全部18个同相导数和18个正交导数。该设备具有很大的攻角范围和侧滑角范围,并具有足够大的振动频率和振幅范围,实验风速可从20米/秒至50米/秒。强迫振动风洞实验的全过程都在微机控制下自动、实时地完成。通过某教练机的强迫振动实验和不同方法的对比实验表明,该设备能较好地测量出组合形式的全部阻尼导数、交叉导数和交叉耦合导数。     

大攻角高速空间密度场显示

在大攻角空气动力学中，三维、非对称流场的密度显示是十分重要的。本文描述在超声速风洞中，应用光学纹影干涉层析术，显示大攻角钝锥的三维密度场的过程，