亚音速位势方程的并行计算

将复杂组合体亚音速气动力特性串行算法改造成并行算法，并在PAN95并行机上进行了数值试验。试验表明，并行计算不仅缩短了计时，而且可增加到18000个面元，计算结果与实验结果吻合较好。     

高阶面元法中的连接条件和相关控制点的选取方法

本文针对高阶面元法在复杂组合体中的应用,发展了相关控制点的选取方法和网片连接条件。应用该连接条件和相关点选取方法,计算复杂外形的飞机气动力,在部件连接部取得了较满意的结果。本方法实用、灵活,工程应用性较强,能适应不同的气动布局。     

天宫飞行器低轨控空气动力特性一体化建模与计算研究

对非规则板舱组合体天宫飞行器300~200 km低轨道飞行过程空气动力特性一体化计算建模,提出考虑复杂构型物面遮盖效应面元解析法与经修正的Boettcher/Legge非对称桥函数,发展基于三角形面元逼近复杂外形通用处理方法,建立适于天宫飞行器复杂物形处理与面元气动力系数计算规则;将DSMC方法与求解Boltzmann模型方程气体运动论统一算法应用于天宫飞行器简化外形,进行气动力当地化关联参数计算修正,建立针对大型复杂结构天宫飞行器低轨道飞行控制过程空气动力特性一体化快速算法与程序软件。对大尺度圆柱体外形与天宫飞行器300~200 km不同高度变轨飞行过程不同迎角与侧滑角及帆板平面与本体主轴不同夹角复杂构型气动力特性计算分析验证,表明天宫飞行器在200 km以上低轨道飞行控制过程中所受空气动力系数随飞行高度发生显著变化(8%~50%),证实长期在轨运行的大型航天器若采用统一固定的气动力系数,误差累积巨大,需要采取防护措施,低轨道飞控大气阻力仍是制约航天器定轨预报精度最关键因素。     

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

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

运输类飞机翼身组合体参数化及快速气动力计算研究

飞机方案设计阶段,迫切需要一个能够快速评估、计算不同参数组合的翼身组合体基本气动力特性的分析平台。为满足这一实际工程设计需要,本文提出运输类飞机参数化概念和具体研究方法,结合运输机气动布局设计的基本准则,应用计算机图形学进行辅助设计,利用丰富的参数曲线,如NURBS曲线、贝塞尔曲线、拉格朗日插值等方法描述复杂的翼身组合体外形,开发了快速分析计算翼身组合体气动力特性分析平台,实现了从总体方案初步设计到MGAERO软件数值模拟的无缝过渡,并结合具体算例验证了方法的可行性和工程实用价值。     

机翼-副翼-调整片的非定常气动力和颤振计算

本文应用偶极子格网法对一个机翼-副翼-调整片颤振模型进行了三元非定常气动力和颤振计算,还研究了网格数目对结果收敛性的影响。计算结果表明本方法具有较好的收敛性,和试验结果相比较表明本方法有一定的工程精度,在飞机设计中能用于操纵面-调整片构型的非定常气动力和颤振计算。     

类乘波体飞行器的气动力工程计算

为估算高超声速类乘波体飞行器的纵向气动力,以三角形面元逼近飞行器外形,根据飞行器迎/背风面情况分别采用达黑姆-巴克法、普朗特-迈耶耳法、切锥法和膨胀波方法计算无粘气动力;通过经验公式估算粘性阻力,并考虑飞行器主要部件之间的气动干扰情况,计算了翼片之间、翼身之间的气动干扰因子,得到整个飞行器的气动力。为验证该方法,以某飞行器为例进行了计算,计算结果与CFD吻合。     

增升状态下的全机位流数值模拟

研究了势流理论下的全机增升装置复杂外形气动力数值计算方法,该方法以低阶面元法为基础,具有计算量小、速度快的优点.复杂流动的粘性效应主要通过附面层粘性修正技术加以体现.采用该方法对某增升装置翼身组合体外形压力分布特性进行计算,并和风洞实验数据相比较,验证了方法的正确性和实用性.在此基础上对某型飞机真实复杂外形的压力分布特性进行了计算,为复杂外形飞行器增升装置的初步设计、方案选型与评估提供了非常理想的工程技术手段.     

民机翼身组合体参数化研究

飞机方案设计阶段,迫切需要一个能够快速评估、计算不同参数组合的翼身组合体基本气动力特性的分析平台。为满足这一实际工程设计需要,本文提出运输类飞机参数化概念和具体研究方法,开发了快速分析计算翼身组合体气动力特性分析平台,结合具体算例验证了方法的可行性和工程实用价值。     

亚音速复杂组合体表面压力分布及气动力特性数值计算方法

本文从亚音速小扰动位势方程出发,给出了计算任意复杂组合体表面压力分布和纵向气动力特性的基本解的数值计算方法。把全机表面划分成许多网格,在每个网格上布置基本解以模拟物体,用物面上没有法向气流穿透的边界条件决定基本解的未知强度,从而计算出物面压力系数,积分出总气动力特性。该方法适用于任意形状的机身及任意平面形状、任意厚度、弯度和扭转分布的机翼(尾翼),机翼(尾翼)可以带有两段后缘操纵面。为降低对电子计算机存贮量、速度和稳定时间的要求,本方法采用了机翼、机身、尾翼部件之间的循环迭代求解。     

基于修正偶极子格网法的跨音速颤振分析

针对工程上翼面的跨音速颤振计算存在的问题，阐述了修正偶极子格网法的跨音速颤振计算方法。采用较为精确的定常气动力数据，对工程上常用的偶极子格网法（DLM）的气动力影响系数矩阵（AIC）进行修正，考虑了翼型、厚度及粘性等影响非定常气动力的因素，并将改进后的非定常气动力用于某超临界机翼的跨音速域的颤振计算，结果较好地预测到了跨音速“凹坑”现象。     

直升机旋翼/机身非定常气动干扰数值分析

基于非定常面元/时间步进全展自由尾迹建立了旋翼/机身非定常气动干扰分析方法。方法中耦合了非定常面元法和时间步进自由尾迹,以准确模拟旋翼非定常气动力、旋翼尾迹及桨叶对机身的非定常干扰效应。为模拟前飞状态下具有升力的机身,将机身离散为汇/偶极子面元,并采用涡线镜像法模拟旋翼尾迹靠近机身表面产生的加速效应。通过计算前飞状态的Maryland、ROBIN(Rotor Body INteraction)旋翼/机身干扰下的非定常压力分布,并与可得到的实验值、CFD计算结果对比,验证方法的准确性。随后分析前飞速度、旋翼与机身距离对旋翼/机身非定气动干扰的影响。计算结果表明机身头部和中部非定常压力主要受桨叶的通过性影响,而机身尾梁主要受尾迹/机身干扰影响,机身非定常气动力频率为桨叶片数的倍频。随前飞速度的增加,机身非定常压力幅值增加,尾梁压力幅值先增加后减小;增加旋翼与机身距离将减小机身和尾梁非定常压力幅值。     

全弹组合体超音速颤振的一种工程计算方法

1.引言 由机翼、尾翼和弹身组合布局型式的导弹各部件间存在复杂的气动干扰,颤振分析时只考虑孤立部件就不够了,且在颤振型态上也有较复杂的组合情况。要求寻找一种有效又切实可行的工程计算方法。1964年间,北京空气动力研究所提出借用定常的干扰因子法计算非定常翼身干扰和串置翼下洗。本文在此工作基础上作了若干改进,气动力方面     

高超声速气动力及激波位置快速计算方法研究

采用快速计算方法进行高超声速气动力计算时,影响计算精度的关键问题主要在于模型面网格的划分和计算方法的选取。采用一种灵活实用的结构化面网格划分策略,使得模型的各个部件能分别选择合适的计算方法;发展一种基于近似流线的二阶激波膨胀波方法,该方法可以用于多种具有三维流场特性的部件,不仅降低对使用者的经验依赖,还能提高计算精度;配合激波位置计算方法,可以较为准确地计算模型的激波位置,保证边界层外缘参数的计算精度;粘性力计算使用基于起始面元修正的Spalding-Chi方法和参考温度方法。通过对四个典型算例的计算与分析,表明本文发展的高超声速气动力计算方法具有较高的计算精度,能够作为高超声速飞行器初步设计阶段的气动力快速分析工具。     

一种翼身组合体的气动概念设计

翼身组合体具有较高的升阻比,可进行较大范围的机动,而且还可以提高落点精度、扩大再入走廊、降低热流峰值并降低过载.本文采用模线设计方法设计横截面控制点,借鉴航天飞机气动力工程计算方法发展了一套可以预估翼身组合体飞行器纵横向气动力的工程计算方法.提出并建立了翼身组合体飞行器的优化设计模型并进行了计算,获得了带后掠下反翼的翼身组合体优化方案.对其升阻比特性、质心设计、稳定性问题、滑翔飞行特性及气动热环境进行了预测和讨论.研究表明,带后掠下反翼的翼身组合体方案可以在较小攻角时获得较大升阻比,纵横向稳定且具有较大的滑翔距离和滞空时间,是一种潜在的高超声速机动方案.     

亚超声速位流场中复杂组合体的压力分布与气动特性的数值计算

本文对F.A.伍德华得的亚超声速位流场中组合体的气动分析方法作了进一步的研究和改进。改进表面边界条件方案中前缘涡面的处理;在翼身交接区内,采用改变机身内部过渡涡位置和纵向分块数的方法,改善了超声速流中组合体上的压力分布;采用了实用的物体表面几何处理方法,扩大了计算程序的几何描述能力。本文简要描述了机翼-机身-尾翼(包括垂尾短舱等部件)复杂组合体的压力分布和气动特性的数值计算方法,介绍了多种简单组合体和复杂组合体在亚超声速位流场中的计算结果,并与实验数据作对比。理论计算和实验数据比较一致,复杂组合体的计算结果表明,计算精度基本上能够满足初步选型设计的要求,可供实际使用。     

飞行器完全非定常气动力的一种工程计算方法

本文以谐振运动非定常气动力计算理论为基础,采用偶极子格网技术进行数值计算、有理式拟合、Fourier积分和Duhamel积分等步骤,给出了飞行器任意运动和受任意大气扰动时完全非定常气动力的一种工程计算方法。该方法数学处理方便,得到的气动力表达式简单。算例结果表明:该完全非定常气动力计算方法可行,精度令人满意,具有工程实际意义。     

考虑机身影响的三角翼跨音速气动力设计

将给定翼面目标压力分布的余量修正法（一种三维机翼气动力反设计方法）推广应用于考虑机身影响的小展弦比三角翼跨音速气动力设计。通过一个翼身组合体的设计实例说明了该方法在处理这类设计问题时遇到的困难及其解决办法     

弹性飞机跨声速机动载荷计算方法

基于跨声速非定常气动力求解的重叠场源法,开展了弹性飞机跨声速机动载荷计算方法研究,为现代飞机结构强度设计提供更加可靠和精确的临界载荷计算方法。首先通过采用重叠场源法求解关于62%根弦俯仰振荡的LANN机翼在马赫数为0.822的非定常气动力并与实验结果进行对比,验证了重叠场源法对跨声速激波效应预测和非定常气动力计算的能力;其次应用频域气动力有理近似技术拟合场源法计算得到的广义气动力系数矩阵,建立了机动载荷分析的状态空间模型;然后完成了某型民用飞机俯仰机动载荷分析,研究了俯仰机动飞行情况下飞机机体状态量及飞机部件载荷响应规律。计算结果表明：考虑机体弹性变形后,机翼和平尾气动载荷响应最大值分别减小了5.1%和10.6%,升降舵气动载荷响应最大值增大了16.2%,在飞机结构强度设计中必须考虑机体弹性效应对飞机部件载荷的影响。     

直升机旋翼的非定常可压缩升力面气动力

本文得出了三元可压缩非定常连续尾迹面条件下直升机等速前飞时旋翼压力分布的解答,使压力偶极子沿尾迹面运动而不沿桨叶的真实轨迹移动以满足尾迹条件,叠加移动压力偶极子获得翼面作复杂运动时三元非定常可压缩的强奇性积分方程,得出了处理强奇性的方法使之成为带强峰值的连续积分方程,并得出了求解该方程的数值解法。在小型计算机上获得的简单算例证明了方程与解法的正确性。