复杂组合体亚音速气动力计算改进的面元法

发展一种低阶面元法计算复杂组合体的亚音速气动力。本方法采用两种形式的速度势的内部第一类边界条件，通过对速度势的数值微分得到速度，进而计算气动力参数。本方法克服了早期低阶面元法对面元剖分敏感和高阶面元法的运算量大的缺点，兼有高阶面元法的精度高和早期低阶面元法的运算量小的优点。用本文给出的方法，计算了多个模型（包括型号）的气动特性，计算结果和实验结果有很好的一致性，表明本方法灵活、稳定、准确和有效。     

面元法求解的一种新尝试

如何提高面元法的求解效率是面元法应用研究的一个重要课题。本文利用面元法积分核的变量分离技术及网格组合的网片串技术,对飞机气动力计算方面应用颇广的求解技术做了一定的改进,提出了一种面元法的快速算法,从在DEC7620机上运行的一个实例来看,其结果是令人满意的。     

多段翼型分析和反设计位流计算

本文采用第二代低阶面元法给出了多段翼型分析和反设计位流计算方法。按照曲率变化情况划分板块可很好地提高计算效率。     

适用于整体求解吊舱推进器的定常面元法

为了方便开展吊舱推进器设计和混合式CRP(Contra-Rotating propeller)的面元法分析,提出了一种将吊舱和螺旋桨作为整体进行定常面元法求解的可行办法,对相关的面元法基本积分公式进行了分析,说明了定常面元法整体求解的结果与实际结果之间的关系。从吊舱桨尾流区的实际流场状况和计算效果出发,分析了整体求解吊舱推进器的面元法尾涡模型,并对此种尾涡模型的修正予以详细说明。在进行相关的理论分析后,编制了整体求解吊舱推进器的定常面元法程序,对吊舱推进器进行计算分析,并与试验结果和RANS(Reynolds-averaged Navier-Stokes)方法计算结果进行了对比,最后对本文所提出的定常整体面元法求解方法进行了误差分析。计算结果表明,整体面元法求解吊舱推进器的误差与试验值相比在5%以内;在计算耗时方面,将吊舱和螺旋桨分开的迭代求解面元法是整体求解面元法的1.7倍。     

基于最小F-范数修正面元法的高速直升机平尾静气弹计算

与常规构型直升机相比,高速直升机平尾面临更加明显的静气弹效应,工程面元法是静气弹数值计算领域广泛使用的气动建模方法。为了改进面元法计算精度方面的缺陷,同时保持其形式简单且计算高效的优势,提出了一种最小F-范数修正面元法。利用以原始刚性模型和预设扭转模型的CFD数据为基础获得的修正系数矩阵,可使平尾弹性变形状态下气动载荷的面元法与CFD计算结果保持较高一致性。基于所提修正面元法,结合结构有限元模型开展了高速直升机平尾的静气弹数值计算,证明了所提方法在静气弹分析中的适用性和有效性。     

机翼／机身低亚声速最大升力系数预测方法

本文研究了大展弦比机翼／机身低亚声速最大升力系数的半经验预测方法。该方法采用便宜可靠的势流面元法／分离泡／边界层粘性耦合技术与最大升力状态下流动现象的观察经验（压力偏差律）相结合。其中势流面元法采用第二代低阶面元法，每个面元上常值源汇和常值仍极子奇点分布；机翼后缘分离流采用简单有效的虚拟固体边界包围分离区模拟方法。它包括二维和三维间的信息交换，以及粘性／势流间的耦合迭代。ＧＸ－Ⅱ机翼／机身算例表明该方法能够计算机翼／机身详细的亚声速气动载荷和表面压力分布。计算的Ｃｌ－α曲线结果与实验数据有较好的一致性。     

高阶面元法计算多段翼型之位流

本文给出了一种求解多段翼型位流的高阶面元法。面元和奇点的分布分别采用了二阶和三阶的形式。由于高阶面元法较为精确地反映了绕流体的几何形状和奇点表面分布,本文给出了比低阶面元法更为精确的解,并且消除了低阶面元法计算中经常出现的翼型后缘处解的波动。     

基于面向对象的面元法

本文介绍了八十年代出现的面向对象的程序设计方法，讨论了用这种方法开发空气动力数值计算中的面元法的优点，并详细介绍了采用面向对象的程序设计开发面元法的设计思想。     

CLT桨敞水性能的势流预报方法

为了提高面元法对尾流收缩叶梢负载（CLT）桨敞水性能的预报精度,通过修正常规面元法中的螺旋桨计算模型、尾涡模型和网格,改进预报CLT桨敞水性能的面元法计算理论。以CLT标准桨（P1727）为研究对象,从CFD方法模拟结果中,提取不同进速下的尾涡螺距和尾涡半径,构建双尾涡模型;提出混合网格划分方法,提高计算精度和效率。将面元法预报结果与国际拖曳水池会议（ITTC）实验数据对比,结果表明,在设计进速附近,推力系数和转矩系数的相对误差在5%以内,效率的相对误差在2%以内,满足计算精度。     

基于修正面元法的机翼焦点位置变化研究

针对静气动弹性问题中机翼焦点位置求解精度与效率的矛盾,提出采用分段斜率修正面元法来计算气动力,该方法以刚体机翼在不同迎角下的CFD数据作为修正基础,对弹性变形后的机翼下洗采用分段修正,克服了传统面元法无法考虑非线性因素的缺点。通过M6机翼气动焦点位置的计算表明:相对于CFD方法,该方法保持了较高的准确度并且计算效率较高;相对于传统面元法,该方法能够准确地计算焦点位置和弹性机翼焦点位置变化量,适用于工程设计。     

机翼平面形状参数摄动面元法

在内部Ｄｉｒｉｃｈｉｅｔ边界条件的基础上,藉助于对奇异积分的解析微分,可用解析法计算偶强密度对机翼平面形状乡数的敏感性偏导致（建立一般敏感性方程所需）。机翼平面形状参数摄动后其偶强分布即可藉线性外推得出。从而,压力分布、升力和俯仰力矩系数即可快速确定。振动面元法与相应的低阶面无法计算结果表明,该方法有良好的准确度且使计算工作量大大减少（摄动外推所需ＣＰＵ时间比相应的低阶面元法少两个量级）。     

不同截面机身的RCS及气动特性研究

对三个“板块”多边形截面机身及圆截面机身模型进行了RCS隐身特性和低速气动特性的测试与研究,发现多边形截面机身不但具有良好的隐身性能,特别是机身侧面,无论是水平极化或是垂直极化,其RCS值均比相同截面面积的圆截面机身降低一个数量级,而且其气动特性不比圆截面机身差,其升力特性与最大升阻比均比圆截面机身好。最后还给出了考虑多边形侧边缘夹角的横向绕流粘性效应的经验修正公式,可供多边形截面机身的气动特性估算参考。     

计算飞航导弹RCS的混合方法

本文提出了一种计算复杂形状物体RCS的算法。这个方法在应用物理光学原理求解散射场时,把散射体表面分割成许多小曲面元,用平面板块柔代替曲面元,求出板块元的散射场矢量之和作为物体散射场的近似值。由于求板块元散射场的积分运算转化成了代数运算,使这种方法具有运算速度快,应用灵活的特点。经与实际测量结果对比,计算结果是较为满意的。     

基于偏导数的全局灵敏度指标的高效求解方法

在全局灵敏度分析领域,基于偏导数的全局灵敏度指标由于其优良的特性得到了广泛的关注。针对目前求解基于偏导数的全局灵敏度指标计算效率低的问题,提出了一种高效的求解方法。该方法利用乘法降维近似地将响应函数展开为连乘积的形式,从而将求解偏导数全局灵敏度的高维积分问题转化为一维积分的连乘积,在保证计算精度的前提下大大降低了求解基于偏导数的全局灵敏度指标的计算量。在求解特定点的偏导数时,采用了复数步长方法,提高了计算精度。最后,通过算例验证了所提方法的准确性和高效性。     

水轰五飞机大型整体壁板的设计与制造

 本文介绍水轰五飞机大型整体壁板的结构设计和经济效果及其数控加工,喷丸成形和毛坯扁挤压制造工艺特点。     

复合材料加筋壁板压缩稳定性工程算法验证研究

复合材料加筋壁板的主要失效模式是丧失稳定性,为确保其结构安全同时有效降低结构重量,有必要 对其稳定性校核时的边界条件和蒙皮受载边宽度选取进行系统的研究。通过对两种边界条件、三种蒙皮受载 边宽度选取进行研究,提出考虑长桁影响的蒙皮受载边宽度选取方法,进行加筋壁板稳定性计算,将计算结果 与试验结果进行验证。结果表明:边界条件选取四边简支,蒙皮受载边宽度选取考虑长桁影响的工程稳定性计 算方法准确度最高,计算结果与试验结果最小偏差0.29%,最大偏差7.16%,且全部低于试验值;该方法用于 结构强度稳定性校核是安全可靠的。     

Dirichlet边界条件低阶面无法的改进

对基于内部Ｄｉｒｉｃｈｌｅｔ边界条件的低阶面元法的速度分布提出一个改进算法。由极限过程结合高阶奇异积分的解析计算直接计算物面速度分布取代传统的偶强梯度法,以避免偶极强度边缘外推的困难。改进算法改善了ＶＳＡＥＲＯ分析中诱导阻力系数估值过低的问题。对后掠翼、前掠翼和翼身组合体等算例的计算表明,计算结果的准确度与高阶面元法相近。此外,直接积分法还可正确地计算外部流场速度分布。这是偶强梯度法不能完成的。     

耦合高效配平策略的旋翼气动特性分析方法

建立了一个耦合高效配平策略和计算流体力学（CFD）技术的旋翼气动特性分析方法.采用CFD求解器进行雅克比矩阵和旋翼流场的计算,建立了一个包含直接配平法的旋翼气动特性分析平台;引入差量法思想对配平策略进行简化,发展了一种差量配平法,解决了直接配平法中雅克比矩阵计算效率低下的问题,显著缩短了计算时间.旋翼流场求解中控制方程采用Navier-Stokes方程,空间离散格式选取2阶迎风Roe格式,时间推进采用隐式lower-upper symmetric Gauss-Seidel (LU-SGS)格式,湍流模型为Baldwin-Lomax (B-L)模型.分别采用建立的基于直接配平方法和差量配平方法的气动特性分析平台对AH-1G和Helishape 7A直升机的旋翼进行了配平计算和研究,对比分析了两种配平方法的计算效率和精度及背景网格的网格量对配平结果的影响.结果表明:与直接配平法相比,耦合差量配平法的气动特性分析方法计算圈数更少,可节省至少60%的计算时间;使用疏化的背景网格虽然不利于捕捉流场细节,但是可以在兼顾配平结果准确性的同时,降低配平计算代价,提高旋翼配平和流场计算的效率.      

飞行器非圆因子和翼身干扰因子修正系数的数值解

提出了飞行器任意形状截面表观面积的数值计算模型和方法,从而求解非圆因子(即非圆剖面与当量圆剖面局部法向力系数之比)和翼身干扰因子的修正系数。计算结果与解析解进行了比较,可看出该计算方法具有精度高和广泛通用性的特点。它对复杂飞行器大迎角(含小迎角)气动特性工程估算有现实的使用价值。     

超低轨航天器气动特性快速预测的试验粒子Monte Carlo方法

超低轨(LEO)卫星气动特性的快速准确计算是对其进行轨道预测和控制的关键输入条件。基于真空技术领域中计算管道分子流率的试验粒子Monte Carlo(TPMC)方法,结合自由分子流理论,发展了一套快速准确预测低轨卫星气动特性的TPMC方法,给出了其模拟步骤及主要关键技术点,并采用该方法模拟了带电池翼超低轨卫星的气动特性和航天器典型构件之间的多次反射效应。结果表明:TPMC方法在计算超低轨航天器气动力、力矩时具有较高的可靠性和对工程复杂外形的适用性;该方法能够准确模拟自由分子流理论无法求解的多次反射问题,给出正确的气动力系数;该方法的计算速度比直接模拟Monte Carlo(DSMC)快3~4个量级,存储量要求也比后者低1~2个量级,是超低轨航天器气动特性快速准确预测的一个理想方法。