高阶间断有限元法的并行计算研究

根据间断有限元法的数据结构特点,基于METIS网格分区技术,设计并行计算策略,在非结构网格上实现了并行高阶间断有限元法。控制方程的数值通量项使用Local Lax-Friedrichs(LLF)格式计算。设计了并行的牛顿-块高斯赛德尔法(Newton-Block GS)来加速收敛,提高迭代效率。并行性能分析表明,所设计的并行算法能够得到较好的加速比和并行效率,有效地节省计算时间,合理分配内存。这使得采用高阶间断有限元法计算更为复杂的问题成为可能。     

常开洞工业厂房风致内压试验研究

通过一系列不同条件下的单个常开洞单跨厂房刚性模型风洞测压试验,分析了风向角、开洞位置、开洞形状以及开洞大小对内压平均值与脉动值的影响,并与基于稳态理论与非定常理论的内压预测值进行了比较。研究发现,单个常开洞厂房内压的空间分布是均匀的,几乎处处相等,仅在开洞附近测点的风压系数值稍有差别;内压的脉动特性主要来源于来流湍流作用、旋涡脱落与尾流干扰以及赫姆霍兹共振等三个方面的影响;内压的均值与均方根值受开洞条件以及风向角等因素影响较大;稳态理论计算得到的内压系数均值与试验所得非常接近,但是均方根值的预测偏保守。非定常理论推导得到的内压响应方程在风速方向正对开洞时可以给出良好的预测,但由于受旋涡脱落、尾流干扰等影响,斜风向时运用该方法预测内压是不合适的。     

应用多块对接结构网格方法的直升机涵道尾桨气动特性分析

基于湍流Navier-Stokes方程,建立了一套直升机涵道尾桨流场及气动特性的CFD分析方法.在该方法中,针对直升机涵道尾桨流场的特点以及动量源方法对网格系统的要求,提出了一种多块对接网格生成方法,分块贴体网格采用求解Poisson方程的方法生成.通过结合涵道尾桨的运动方式、几何特征及气动特征,建立一个包含动量源项的N-S方程的涵道尾桨流场计算方法和迭代流程;为较好地捕捉涵道壁附近存在的旋涡及分离现象,采用了一方程S-A湍流模型.通过对ONERA M6机翼的绕流特性的计算分析,以及对涵道尾桨在悬停状态下的拉力、诱导速度、载荷分布等的计算分析,验证了该CFD方法的有效性.在此基础上,采用该方法进行直升机悬停、侧飞和前飞状态下涵道尾桨流场与气动特性的数值分析,得到了关于涵道尾桨流场和气动特性的一些有意义的结论.     

新型超高强Al-Zn-Mg-Cu铝合金热压缩变形的流变应力行为

采用热力模拟试验方法进行热压缩变形试验,研究了一种新型Al-Zn-Mg-Cu合金在变形温度为300～450℃,应变速率为0.001～1s-1,压缩变形程度为50%条件下的流变应力行为.结果表明,变形温度和应变速率对合金流变应力的大小影响显著,流变应力随变形温度的升高而降低,随应变速率的提高而增大.采用统计回归方法建立了...     

改进联邦滤波器在小卫星姿态确定中的应用

采用微小卫星姿态运动学、动力学方程,以磁强计／太阳敏感器为姿态敏感器进行姿态确定,研究了一种基于融合反馈的改进联邦滤波算法,提出了一种基于观测新息修改的信息分配原则。给出了由四元数描述卫星姿态的误差状态方程,并对观测矢量进行推导,获得了基于角速度的观测方程,丰富了观测模型。仿真结果表明：基于角速度观测方程的引入,增强了观测信息的利用率,进而提高了姿态确定的精度,降低了计算量,保证了可靠性。     

悬停状态直升机桨叶扭转分布的优化数值计算

建立了一套基于高精度计算流体力学(CFD)技术和代理模型优化算法的旋翼气动外形设计方法。在该方法中,旋翼流场气动性能的计算采用了基于Navier-Stokes/Euler方程的CFD方法,并根据流场特点、精度和效率的要求采用Baldwin-Lomax(B-L)湍流模型,通量计算采用Roe-MUSCL格式进行。为提高网格生成质量和便于流场控制方程的求解,将流场分成两个区域,即围绕旋翼的黏性区和无黏的背景网格区。其中,桨叶网格使用了基于二维翼型网格的参数化方法生成,数值计算结果表明该方法有效地提高了网格生成质量及效率。在参考旋翼流场及桨叶细节流动分析的基础上给出设计变量及范围,有效减小了优化问题的规模;为满足优化和机理分析的需要,将基于置换遗传算法优化的拉丁超立方(PermGA LHS)方法和径向基函数(RBF)的代理模型优化方法引入到桨叶外形的优化设计中。首先以Helishape 7A旋翼为算例,检验了数值模拟方法的准确性。然后,应用所建立的优化方法针对旋翼负扭转分布进行了优化计算,结果表明优化后的旋翼悬停气动性能比优化前有了明显提高。     

固体推进剂/衬层界面脱粘裂纹的三区域界面层模型

将固体推进剂/衬层粘结界面层按照剪切模量的变化分成三部分,由此建立了其界面脱粘裂纹的三区域界面层模型。应用Fourier积分变换和传递矩阵方法,推导并数值求解Cauchy奇异积分方程,得到了反平面剪切状态下裂纹问题的半解析解,并讨论了界面层模量和厚度对断裂参量的影响。结果表明,界面层模量随贮存时间的延长而降低时,应力强度因子逐渐减小,但裂纹张开位移变化不大;界面层厚度对应力强度因子的影响则不明显。     

点隐式龙格-库塔方法的应用研究

为了提高求解器的效率,在显式龙格-库塔时间推进的欧拉方程求解器之上,发展了点隐式龙格-库塔时间推进格式。给出了其推导过程和非结构网格下中心格式和迎风格式（Roe和Van Leer格式）预处理算子的构造方法。NACA0012翼型和RAE2822翼型的跨音速无粘流动模拟表明：与显式龙格-库塔方法相比,方法能提高求解效率且内存需求相当,具有一定的工程应用价值。     

结合响应面法与网格法的低速翼型优化设计

本文首先了介绍响应面法和网格法的基本原理,翼型绕流场分析使用有限体积法空间离散和RUNGE-KUTTA时间迭代格式进行N-S方程计算,由翼型表面压力和摩擦力数值积分得到升力和阻力.然后,以NACA0012和GAw-1为算例,结合响应面法和网格法进行低速翼型的减阻优化设计,算例表明:在应用响应面法进行低速翼型优化时,采用...     

机动飞行器中等攻角高超声速无粘绕流数值模拟

本文采用ＮＮＤ差分格式，利用推进－迭代方法，数值模拟了机动飞行器的中等攻角高超声速无粘绕流流场，发展了一种适用于处理中等攻下背风面流场的技术，大量计算表明，结果准确、可靠。气动力力系数和壁面压力分布与实验数据一致，流场波系结构模拟正确。