栅格翼与机身发动机组合体的气动特性计算

本文给出了栅格翼和机身、机身和发动机以及机身的台阶与其前后柱体等复杂气动力干扰的计算方法，考虑栅格翼与机身气动干扰时，假定栅格翼处在机身的向流动的上洗流场中，用翼片理论计算有机身影响时栅格翼上载荷分布，考虑机身和发动机气动干扰时，有杉一个机身和上个发动机组成的十字型五体组合体模型，用保角变换算出干扰因子，考虑机身台阶与其前后柱体的气动干扰时，采用锥柱体和截锥结合实验修正的方法，用本方法计算了栅格翼     

栅格翼和机身组合体的气动特性计算与分析

给出计算栅格翼和机身组合体亚超音速气动力干扰理论方法.亚音速采用涡格法计算栅格翼气动力,用一个二维偶极子和圆柱映像涡系模拟机身的上洗流场和栅格翼迁移到机身上气动载荷.超音速假定栅格翼处在机身横向流动的上洗流场中,用翼片理论计算机身对栅格翼气动力干扰的修正因子.用本文方法计算了几个栅格翼和机身组合体的气动特性的到了比较满意的结果.      

气动力计算的积分技术讨论

讨论了"尾迹面积分法"的气动力计算方法及其改进技术.基于欧拉方程求解流场,通过对飞行器以外的一个控制体采用动量定理而得,是传统远场积分法的改进.在飞行器所在流场的下游一垂直于来流速度的截面上进行积分计算,得到升力、诱导阻力、激波阻力和总阻力.这种方法的优点是有利于外形复杂物体的气动力积分计算,并可将总阻力按产生的物理机理进行分解,以使设计师对飞行器的气动特点有更为明确的了解.详细讨论了影响这种方法计算精度和效率的多种因素及解决途径.各种数值模拟结果证明了该方法和改进技术的正确性和实用性.      

基于导重法的自重载荷下悬臂梁结构拓扑优化

针对自重载荷悬臂梁结构拓扑优化末端区域材料分布不收敛的问题，提出了一种拓扑优化方法。根据虚功相等原理，利用四节点矩形单元的形函数、单元体积密度和质量的关系建立了载荷等效方法。根据优化模型的Kuhn-Tucker条件推导出了导重准则，以及目标函数的灵敏度公式和考虑自重载荷拓扑优化的迭代算法。针对自重载荷作用下悬臂梁结构拓扑优化存在的末端区域材料分布模糊问题，研究了变密度法和非结构质量相结合的求解策略，揭示了典型因素对拓扑结构的影响规律。结果表明:所提方法能够解决自重载荷下悬臂梁末端区域材料分布模糊问题。      

非结构三角形网格的一个ENO型有限体积方法

基于双曲型守恒律方程,对非结构三角形网格给出了一种ENO(Essentially Nonoscillatory Scheme)型有限体积格式,方法的主要思想是先对每一个三角形单元构造一个加权的二次插值多项式,而在计算交界面的流通量时采用了两点高斯积分公式以保证格式的整体精度,时间离散采用三阶TVD Runge-Kutta方法.最后给出了该格式收敛的数值阶,并对前台阶问题进行了计算.      

非结构三角形网格的一个ENO型有限体积方法

基于双曲型守恒律方程,对非结构三角形网格给出了一种ＥＮＯ（ＥｓｅｎｔｉａｌｙＮｏｎｏｓｃｉｌａｔｏｒｙＳｃｈｅｍｅ）型有限体积格式,方法的主要思想是先对每一个三角形单元构造一个加权的二次插值多项式,而在计算交界面的流通量时采用了两点高斯积分公式以保证格式的整体精度,时间离散采用三阶ＴＶＤＲｕｎｇｅＫｕｔａ方法．最后给出了该格式收敛的数值阶,并对前台阶问题进行了计算     

最优气动力辅助空间拦截和交会

航天器从高轨道向位于低地球圆轨道的靶目标实施拦截或进行交会 ,并用间接求解的直接伴随方法 ,即用D型拉格朗日方法求解其优化问题。拦截条件是指在最终时刻两者的位置相同 ,而交会条件是指不仅位置相同且最终的速度矢量亦相等。文章讨论了最小时间与最小脱靶量 ,或最小时间与最小燃料消耗的复合性能指标 ,以及推力协同机动的优化解。推力协同下 ,对于特定情况 ,如在热流约束边界部分推力弧是奇异的 ,气动力控制是正常的。这种复杂性和奇异问题的求解在拦截和交会问题中都是值得慎重处理的。     

谱载荷下三维构件疲劳分析的损伤力学方法

建立了工程中常见的三维构件的一个损伤力学守恒积分,并利用此积分的守恒性与小范围损伤的条件,推导出计及损伤耦合效应时集中的应力与应变所应满足的方程.根据以损伤驱动力表示的损伤演化方程,推导了在谱载荷作用下,疲劳裂纹萌生寿命预估的解析表达式,建立了以等效应力表示的非分离变量型损伤演化方程,并在短周期加载条件下得到积分形式的疲劳寿命闭合解.对某飞机起落架为代表的三维构件在谱载荷作用下的疲劳寿命进行了预报.      

基于压敏漆的多羽流气动力效应试验研究

多股羽流相互作用会形成复杂流场形态的干扰羽流，为了对干扰羽流开展气动力效应试验研究，采用压敏漆（PSP）表面压力光学测量技术，对以常温空气为工质的单喷嘴羽流和双喷嘴干扰羽流撞击平板模型的气动力进行了高分辨率的全场测量，分析了喷嘴和平板间不同入射距离和入射角对羽流气动力作用强度和范围的影响。试验结果表明，以高透氧聚合物为基层的压敏漆在羽流气动力测量中具有快速响应的特性，能分辨Pa级的微小压力变化。相对传统离散测压孔方法，压敏漆温度敏感度低，能捕捉大梯度的压力变化，准确描述羽流气动力效应。通过对比，发现干扰羽流对气动力具有明显增强作用，且增强作用随着入射距离的降低而减弱，干扰羽流的气动力不能直接用单股羽流的气动力进行线性叠加；羽流相互作用增强羽流返流，在航天器设计中需要考虑多股羽流相互作用后的羽流防护问题。      

不确定气动载荷计算的区间分析方法

研究微小型飞行器在大气飞行时受到连续阵风干扰后气动升力的预测问题,考虑von Karman谱参数的不确定性,将其用区间向量定量化,基于区间扩张理论和Taylor级数展开,并结合非定常气动力理论,提出了一种考虑飞行环境不确定性的非定常气动载荷计算的区间分析方法.该方法只需要知道不确定参数所在范围的界限,而不需要其他概率信息,为解决不确定气动载荷计算这种复杂问题提供了一个途径.通过数值算例,将区间分析方法与概率方法的结果进行了比较,显示了区间分析方法是可行并有效的.     

捆绑火箭气动载荷分布不确定性分析

针对带捆绑火箭气动载荷分布受飞行状态及本身外形参数变化影响存在波动的现象,提出了依据多项式混沌理论对捆绑火箭气动载荷分布特征进行全局灵敏度分析及不确定性量化的方法,并以两助推构型火箭为例对所提方法进行验证。首先,提出了捆绑火箭气动载荷分布不确定性分析的方法,并给出仿真分析流程。其次,以两助推构型火箭为例对所提方法进行验证,建立火箭气动外形参数化模型,验证气动特性分析结果。最后,对该模型开展影响因素灵敏度分析及载荷分布不确定性分析,得到了不同因素的影响程度,以及气动轴力和法向力的不确定性分布形式,分析了流场流动情况及气动载荷波动的主要原因。分析结果为捆绑火箭气动载荷波动控制提供了一定参考,通过定量描述气动载荷分布不确定性,可以有效降低安全系数冗余,为开展精确结构设计提供依据。      

螺旋桨/大柔性机翼静气动弹性快速分析方法

旋转的螺旋桨滑流掠过机翼将使机翼的气动特性发生改变,在高空超长航时无人机的设计中有必要对大柔性机翼气动弹性问题的螺旋桨滑流影响进行分析.运用Prandtl修正的动量叶素理论分析螺旋桨滑流及面内载荷;采用兰金涡核模拟滑流对机翼的诱导速度;采用三维升力线方法计算机翼定常气动力,利用曲面样条插值方法解决结构/气动耦合问题,并结合非线性有限元静力学计算方法,建立了螺旋桨滑流及面内载荷作用下大柔性机翼静气动弹性问题的快速迭代求解方法.以某大展弦比螺旋桨机翼为例,采用文中所建立方法对其静气动弹性特性进行计算研究.结果表明,旋转的滑流改变了机翼绕流当地攻角,从而影响了机翼气动力和变形分布,且在小前进比时影响更大.所建立的分析方法简便高效,在初步设计阶段有较好的应用前景翼绕流当地攻角,从而影响了机翼气动力和变形分布,且在小前进比时影响更大.所建立的分析方法简便高效,在初步设计阶段有较好的应用前景.      

超高速飞行器平尾大迎角气动弹性特性研究

临近空间超高速飞行器在飞行过程中受到外部干扰作用时会出现大迎角飞行姿态,此时需大角度偏转全动平尾进行配平,带来平尾大迎角下的气动弹性问题。采用计算流体力学/计算固体力学/计算热力学（CFD/CSD/CTD）耦合方法分析了一种超高速飞行器全动平尾的气动弹性特性,重点研究了大迎角下平尾的气动响应及结构变形特点。结果表明:各迎角时的气动力曲线均出现波动,随时间变化逐渐衰减至平衡位置。迎角越大,初始振幅越大,气动力系数减小的比例越大,但随时间衰减得越快。平尾存在弯曲/扭转耦合现象,结构变形导致表面压力分布发生变化,使得整体压力减小、升力系数降低,迎角越大现象越明显。平尾最大应力在迎角30°时达1.2 GPa,已达到所用镍合金材料的屈服强度极限。应在结构设计时在翼轴与平尾接触部位附近加强,或在控制方案设计时限制全动平尾的工作角度。结构发生轴向与法向变形,轴向变形主要由气动热引起,法向变形由气动力和气动热共同引起。      

一种高超声速滑翔再入在线轨迹规划算法

为了提高滑翔再入飞行器响应动态任务的能力并提高其制导系统的鲁棒性，建立了高超声速滑翔再入轨迹规划问题的一种在线求解算法。该算法以一列凸子问题来逼近原问题。通过选择以弧长为自变量并引入对数速度代替速度作为状态，使得动力学方程的非线性大为减弱，同时使动压和热流约束完全成为线性约束。通过使用类似于混合整数规划中割平面的思想处理禁飞区约束，尽可能避免了不必要的计算。通过直接以气动系数和大气密度的乘积作为控制量，构造伪线性控制模型，进一步减弱非线性。非凸约束被适当松弛，以保证子问题的可行性。为了避免过度松弛，利用离线给定的高度和速度上下边界估计对应的参数，以加速收敛。以X-33再入任务为例验证了所提算法的有效性。该算法可以以简单的常值函数为初值并在数次迭代后收敛。      

基于试验气动力的纵向机动飞行载荷分析

提出了一种基于非线性风洞试验数据(简称试验气动力)同步进行纵向机动飞行过程和飞行载荷的计算分析方法.通过曲面样条插值将试验气动力的升力系数和俯仰力矩系数引入纵向机动飞行过程的计算,同时基于试验气动力的压力分布进行全机飞行载荷的静气动弹性修正.弹性载荷的计算采用线性气动力影响系数矩阵处理.采用偶极子格网法提供气动力(简称理论气动力).将基于这两种气动力的计算结果进行了对比.结果表明:基于试验气动力的机动飞行过程更为合理;基于理论气动力的载荷随机动飞行过程的变化趋势有可能与基于试验气动力的载荷变化趋势不一致.      

基于准则的大展弦比飞翼气动设计

从设计实际出发,为切实提高气动性能,开展了大展弦比飞翼无人机(UAV)的气动设计及分析研究.在设计分析过程中,依据飞翼无人机的特征,提出了气动设计准则;基于设计准则,采用更新设计的策略,结合变可信度数值模拟、代理模型优化方法构建了优化设计框架;针对飞翼无人机开展了参数化表达、无限插值网格自动生成以及多轮更新优化,得到了优化推荐构型;应用γ-Re_θt转捩模型方法对优化构型的气动性能进行了细致地验证分析.研究结果表明:通过气动设计,飞翼无人机设计构型很好地契合了设计准则,其巡航升阻比相比最初的原始构型提高了14%,γ-Re_θt转捩模型能较细致地分析大展弦比飞翼的流动特征.      

双尾撑布局飞机飞行载荷分析

对双尾撑布局飞机进行飞行载荷和气动导数分析,重点分析尾撑纵向弯曲刚度和飞行动压变化,对纵向气动力导数、各操纵面效率以及载荷分布的影响.尾撑纵向弯曲刚度变化会引起纵向静稳定性、平尾升降舵的操纵效率及反效速度、平尾的载荷分布的显著改变.尾段刚度尤其是尾撑的纵向弯曲刚度往往会构成双尾撑布局飞机气动弹性设计的薄弱环节,在设计时应该严格控制尾撑的刚度指标.