多层手征介质圆波导的传输特性

阐述多层手征介质圆波导传输特性的递推特征函数解 ,导出了 N层手征介质圆波导传播常数的特征方程 ,提供了若干例子的数值结果     

温度载荷下装药内外径比和长径比对结构完整性影响的规律研究

在固体火箭发动机装药设计中通常将药柱的外径与内径之比定义为m数(m=R/r)为研究装药设计参数m数和药柱长径比在温度载荷下对装药结构完整性的影响, 采用基于Total Lagrangian方法的热粘弹性大变形增量本构关系, 通过选取9种m数以及8种药柱长径比, 共72个计算模型, 对装药结构完整性进行了计算, 结果表明温度载荷下药柱最大Von Mises等效应力、应变值受药柱m数以及长径比共同影响.并得出当长径比大于3时, 药柱最大等效应力、应变值主要取决于m数的结论.      

二维翼型绕流的CE/SE方法

采用时空守恒元和解元CE/SE(space-time Conservation Element and Solution Element method)法,求解二维Euler方程,开展了翼型绕流的无粘数值模拟研究.用非敏感克朗数计算格式消除克朗数过小引起的数值耗散对解的污染,结合当地时间步长法,解决网格不均匀引起的当地克朗数变化跨度大的问题.对NACA0012翼型的无激波流场进行了二维数值模拟,并与AGARD算例做了对比.结果表明:CE/SE方法的计算结果与AGARD结果吻合得很好,为该数值计算方法对翼型绕流数值模拟的进一步应用奠定了基础.      

斜切反喷管性能分析

固体火箭发动机前端斜切反喷管，其结构简单、作用时间短、气动型面具有尖点，并在超音速区有台阶，喷管内存在一系列激波，并伴有流动分离现象。本文从雷诺平均的非定常Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程出发，结合采用Ｂｏｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型，利用时间相关法及ＭａｃＯｃｏｒｍａｒｋ两步显格式求解，模拟了斜切反喷管流场。计算得到的壁面压强分布与风洞吹风实验测得的压强分布相一致。该方法可应用于斜切反喷管的     

二维折角管道内超声速流动的多重网格计算

本文应用Euler模型和隐式TVD格式计算了二维折角管道内的超声速流动问题,为了有效地模拟激波和膨胀波系的相互干扰及壁面的反射,本文试将多重网格法和TVD格式结合起来进行数值求解。考虑到目前绝大部分文章采用V-循环多重网格法,我们对多重网格法中的几种多层次循环形式(V-循环、W-循环)进行了比较性的数值计算,旨在从中发现合适的多层次循环求解的形式,以期充分发挥多重网格法在Euler方程求解过程中的效益。同时,本文也是对多重网格法和TVD格式结合求解的一次初步尝试。     

基于CFD/CSD耦合方法的旋翼气动弹性载荷计算分析

为提高直升机前飞状态下旋翼非定常气动弹性载荷的预估精度,在旋翼气动弹性综合分析方法中引入旋翼CFD模块,建立了一套基于CFD/CSD松耦合分析的计算方法和程序。为高效解决流固耦合方法中由于桨叶挥舞、扭转等弹性变形带来的旋翼贴体网格变形问题,采用基于代数变换方法的网格变形技术,桨叶运动变形量和旋翼气动力信息通过流固交接面传递。旋翼流场分析方法中,主控方程采用耦合S-A湍流模型的Navier-Stokes方程,围绕旋翼流场的网格采用结构嵌套网格方法生成,无黏通量计算采用Roe格式,时间推进采用双时间法。旋翼结构分析中,考虑旋翼配平,基于Hamilton变分原理和20自由度Timoshenko梁模型求解弹性旋翼非线性运动方程。分别对CSD和CFD方法进行验证,在此基础上,计算了SA349/2旋翼桨叶在前飞状态下的非定常气动力、挥舞弯矩、摆振弯矩和扭转力矩,并与飞行测试数据进行了对比。计算表明:CFD/CSD耦合方法可以显著提高旋翼非定常气动弹性载荷的分析精度,精确捕捉桨叶表面压强峰值、激波位置等,表明本文发展的旋翼CFD/CSD耦合方法可以有效地运用到旋翼气动弹性载荷的预测分析中。     

多面体网格在CFD中的应用

介绍了多面体网格技术在计算流体力学(CFD)中的应用.考虑到粘性计算的需要,将混合网格技术引入多面体网格中,发展了一套多面体网格生成及求解技术,并用于M6机翼、DLR-F6复杂外形的气动特性模拟中.计算结果表明,多面体网格技术具有准确度较高、效率较高及适用性广等优点.     

基于全速势方程的超临界翼型设计

基于全速势方程,完成了超临界翼型设计。首先分析了最优化设计和反设计两种常用的翼型设计方法的特点。综合考虑两种方法的利弊,选用最优化设计法。结果显示,设计理论和实施方法完全适合现代翼型设计,设计结果可以满足设计指标要求,诸如巡航状态下的升力系数、阻力系数、升阻比、攻角、翼型最大厚度及厚度位置、巡航效率、压力分布形态、激波强度等等。     

稀疏非结构网格上的亚声速流高精度数值模拟

采用高阶间断有限元法在非结构网格上数值求解二维亚声速Euler方程。数值结果表明,尽管采用的非结构网格非常稀疏,但通过采用真实曲线物面边界和高阶的基函数仍可得到高精度的数值解。另外,对于超低速情况,方程无需经过任何特殊处理就可以得到收敛的数值解。由于采用牛顿一般最小余量法(Newton-GMRES)时通常需要较好的初始值,本文设计了一种阶谱循环过程来提高数值求解的鲁棒性。     

轮履复合式军用地面探测车运动学建模及分析

为更好的对军用探测车运动过程进行控制,提出了轮履复合式探测车运动学模型的构建方法.通过对轮履复合式移动系统结构进行分析,首先建立了探测车三维运动学模型,并推出了其运动方程;其次建立了探测车车体速度与车轮速度之间的雅可比矩阵,利用最小二乘法对车体速度进行求解,提供了获得探测车越障过程中位置和方位的方法.最后对探测车行进过程进行了运动仿真,数值仿真分析验证了运动学模型的正确性及轮履复合式移动系统的运动特性.