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

根据间断有限元法的数据结构特点,基于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%条件下的流变应力行为.结果表明,变形温度和应变速率对合金流变应力的大小影响显著,流变应力随变形温度的升高而降低,随应变速率的提高而增大.采用统计回归方法建立了...     

考虑喷流效应的机载导弹发射及气动干扰数值模拟

研究比较了机载导弹发射时是否考虑导弹尾喷流情况下导弹运动特性及机-弹气动干扰影响.通过动态结构嵌套网格技术实现了导弹运动的模拟,结合刚体六自由度运动方程求解三维非定常Euler方程得到流场信息.导弹推力由发动机燃烧室喷流计算得到;当不考虑喷流效应时,推力通过直接在导弹尾部给定非定常作用力实现.应用此方法模拟了基于类"全...     

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

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

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

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

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

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

跨声速风扇叶片变形对气动性能的影响

为了研究跨声速风扇叶片变形对气动性能的影响,发展了流固耦合方法计算风扇叶片在气动力作用下变形的数值程序,耦合求解了三维可压缩雷诺平均的Navier-Stokes (N-S)方程和描述叶片振动的结构动力学方程,并模拟了跨声速风扇Rotor 67转子叶片在气动载荷下的变形,重点分析了由于叶片变形带来的气动性能的变化.