计算跨音速三维机翼非定常载荷的核函数法

 本文给出了在跨音速流动中,计算三维机翼非定常载荷分布的核函数法和计算结果。此方法是基于线化升力面理论,利用在亚、超音速情况下,联系载荷分布与下洗的积分方程的核函数,采用局部线化法来处理跨音速的非定常绕流问题。在计算中引进了线分布的冲波偶极子,并嵌入适当的正冲波边界条件。 我们从统一的超、亚音速核函数表达式出发,对亚音速区数值积分作了新的简便的处理,对于亚、超音速区相互的诱导下洗给出了自己的数值积分过程。我们算了三个例题,和文献上的数值符合较好。     

亚音速升力面理论的一种新解法

本文提供了亚音速升力面理论的一种新解法,可求解亚音速流中有限翼展机翼的载荷分布。计算中采用传统的弦向级数和展向分条。证明了Mangler二重奇异积分可交换积分次序,从而使积分可解析地求出。对弦向积分的被积函数作多项式分段逼近,可使核函数的表达式成为封闭形式的有限和,导致性态很好的系数矩阵。计算结果与实验及其它方法的结果作了比较,并给出了计算效率。     

有地面干扰的动导数计算方法

提出了一种有地面干扰情况的飞行器动导数计算方法,其基本思路是从非定常升力面气动力理论出发,利用飞行器的低频运动的气动力,通过振动频率的级数展开,选用适当的气动力部分进行积分,得到所需的气动力导数,为阿了面的影响,将传统的镜像法结合到非定常用功力理论理论方法中,计算结果表明,地面干扰随距地面的高度变化而迅速变化,而且某些动导数的符号都发生了变化,说明地面干扰城飞行动力学的分析是应该加以考虑的。     

计算Cauchy和Mangler积分主值的新方法

对于Cauchy积分主值曾提出过一些数值积分法,主要是抽去奇异性方法和内插型或Gauss型方法。对于前者,要求积分核在积分区间上有连续的一阶导数;对于后者,实际上只能求出节点(某类函数零点)处的积分主值,且当积分核不是多项式时     

准静弹性飞机动态特性的计算研究

本文在小扰动和准静弹性假设下对弹性飞机纵向动态特性进行了计算和分析。对弹性变形对动态特性的影响进行了初步讨论。本文采用格林函数法求得飞机的气动力,用正则模态法计入弹性变形影响并引入准静弹性假设,得到了研究纵向动态特性所需的全部刚体飞机气动导数和弹性飞机气动导数。利用这些导数,对准静弹性飞机动态特性进行计算研究。并将计算结果与刚化模型的同一飞机进行比较和分析。     

基于共轭方程法的跨音速机翼气动力优化设计

 设计状态的机翼气动力特性是设计人员最为关心的指标, 应用控制理论设计方法进行了有升力约束情形下跨音速机翼阻力优化设计研究, 根据给定的目标函数推导了相应的共轭方程和边界条件, 研究了共轭方程的数值求解方法, 以及计算目标函数对设计变量的敏感性导数时所涉及的度量矩阵变分求解问题, 研究了流场计算、共轭方程数值求解、敏感性导数求解和拟牛顿优化算法这几个主要方面的有效结合问题, 发展出了一种跨音速机翼气动力优化设计方法, 进行了跨音速机翼气动力优化设计研究验证, 优化后机翼气动力特性有一定程度的改善, 阻力系数能减少20%左右, 而升力系数有所增大, 说明所发展的设计方法是成功的, 该设计方法在跨音速及复杂外形气动设计方面比以往设计方法具有更好的适用性和优越性。     

平尾上反角和静升力对T尾颤振的影响

T尾颤振分析在计算非定常气动力时要计及平尾上反角和静升力等参数的影响。基于片条理论计算由静升力引起的附加非定常气动力，将其引入到用MSC／NASTRAN常规T尾颤振计算得到的广义气动力当中，再用PK法求解T尾颤振，从而建立T尾的颤振工程计算方法。     

等速上仰翼型绕流结构的数值模拟

通过数值求解涡量-流函数形式的N-S方程研究等速上仰翼型的分离流动结构。用ADI方法解涡量方程;用Poisson方程直接法解流函数方程。在计算得到的流场中除可看到实验中观察到的多涡结构外,还发现一些实验中未观察到的现象。将分离流动结构的计算结果与翼型气动力变化相联系,并结合涡动力学理论,定性地解释了上仰翼型产生高升力的原因。     

一种基于共轭方程法求解黏性反问题的简化方法

 对于给定压力分布的黏性气动反问题,考虑到壁面微小扰动造成的压力变化主要由势流作用引起,因此可以简化用以获得目标函数对设计变量敏感性导数的共轭方程。将黏性流场插值到粗网格中作为彻体力模型方程的解,则其相应的共轭方程将以简单的源项取代原方程中复杂的黏性项。由于在粗网格中求解,网格数减少,同时收敛速度加快,简化的共轭方程计算时间可以减少到黏性方程的十分之一。典型的算例结果表明,对于附着的边界层流动简化方法计算得到的敏感性导数具有较高的精度,能够有效完成反问题设计且减少总的计算耗时。     

螺旋桨性能计算的升力线与升力面方法

基于改进升力线方法，编制了螺桨升力线方法计算程序，并以此为基础，在桨尖外一层虚网格和虚控制点，得升力线理论中导出的附着涡和自由涡的积分公式，发展了一种螺旋动性能计算的升力面方法，该方法具有简洁、方便，实用等特点，用本文发展的升力线和升力面方法对多个螺旋桨算例进行了计算，计算结果与实验结果符合良好。     

一种求解核函数方程的新方法—辅助核法

本文提出了求解核函数方程的辅助核法。此法克服了偶极子栅格法对网格分布的限制,具有更大的适用范围。数值计算表明,此法对求解绕翼型的跨声速激波流动是很有利的。     

绕三维物体亚音速振荡位流及其在飞行器动导数计算上的应用

 本文介绍处理不同外形三维物体亚音速振荡绕流一种统一的方法。本方法的主要特点是采用有限元素法直接解由格林定理导出的物面速势积分-微分方程以求得物面的速势分布,然后再用有限差分法对速势进行微分求物面的压力分布。 由于本方法理论上比较严格,适用于复杂外形物体绕流的计算,所得结果又比较准确,因此近几年来在国外得到越来越广泛的应用。在本文中,采用了与有关文献相同的基本方程,但在气动影响系数的计算上略有不同,本文并将这一方法应用于飞行器动导数的计算,所得的结果与实验结果符合。     

渐变结构系统模糊可靠性灵敏度分析的矩方法(英文)

For a degradable structural system with fuzzy failure region, a moment method based on fuzzy reliability sensitivity algorithm is presented. According to the value assignment of performance function, the integral region for calculating the fuzzy failure probability is first split into a series of subregions in which the membership function values of the performance function within the fuzzy failure region can be approximated by a set of constants. The fuzzy failure probability is then transformed into a sum of products of the random failure probabilities and the approximate constants of the membership function in the subregions. Furthermore, the fuzzy reliability sensitivity analysis is transformed into a series of random reliability sensitivity analysis, and the random reliability sensitivity can be obtained by the constructed moment method. The primary advantages of the presented method include higher efficiency for implicit performance function with low and medium dimensionality and wide applicability to multiple failure modes and nonnormal basic random variables. The limitation is that the required computation effort grows exponentially with the increase of dimensionality of the basic random vari- able; hence, it is not suitable for high dimensionality problem. Compared with the available methods, the presented one is pretty competitive in the case that the dimensionality is lower than 10. The presented examples are used to verify the advantages and indicate the limitations.     

基于升力面自由尾迹的直升机旋翼悬停性能参数影响研究

 最大悬停效率(FM_max)作为衡量旋翼悬停性能的常用指标,反映了旋翼能达到的最大悬停效率,但不能反映旋翼在一定桨叶载荷范围内保持高悬停效率的能力,本文给出了旋翼悬停保持能力的定义.为更准确地反映桨叶涡量分布,建立了基于升力面理论的桨叶气动模型;考虑有弯度翼型的影响,将涡量布置于翼型中弧线,随后基于自由尾迹模型、耦合刚性桨叶挥舞运动方程、翼型动态失速模型以及二阶精度时间步进格式建立了升力面自由尾迹方法.通过计算模型旋翼在不同桨尖马赫数下的悬停效率,并与试验数据对比,验证了方法的准确性.相比于升力线自由尾迹方法,建立的升力面自由尾迹分析方法能显著提高旋翼悬停效率计算精度.最后分析旋翼关键设计参数对悬停性能的影响,得到设计参数影响旋翼悬停保持能力的新规律.     

稳定性导数计算方法研究

稳定性导数的准确计算对于飞行器的操稳特性具有重要意义.应用计算流体力学(CFD)方法中的非结构化动网格技术建立能够模拟飞行器做周期性俯仰运动的强迫振荡法,以国际动导数标模Finner导弹为验证算例,获得不同马赫数下俯仰力矩系数的迟滞环曲线,进而计算Finner导弹在不同马赫数下的静稳定性导数和动稳定性导数.结果表明:本文计算得到的俯仰静稳定性导数与试验数据非常接近;在亚音速和超音速范围内,动稳定性导数计算结果与试验值很接近,但在跨音速范围内,本文计算结果与试验曲线的规律性一致,但误差较大.     

高超声速滚转阻尼导数数值模拟

采用非定常Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程描述物体简谐振动流场,并在Ｅｔｋｉｎ理论下给出绕定轴转动时滚转阻尼导数的计算公式,定常流场的计算采用空间二阶精度的交替方向隐式分解的ＮＮＤ格式,非定常流场的计算采用时、空二阶精度的Ｒｕｎｇｅ－Ｋｕｔｔａ多步格式,采用代数和方法生成物体静、动网格。最后给出高超声速下钝体外形滚转阻尼导数的计算结果,以及滚转力矩系数随瞬时振幅的变化曲线。     

桥联重积分方程一种新的数值方法

纤维增强复合材料断裂力学中的桥联问题常常可归结为一个二重积分方程。常见的求解方法有矩阵迭代法和积分迭代法，但这两种方法都产复杂。本文给出一种新的分段插值方法，对此积分方程进行求解。数值计算结果表明，该方法效果令人满意。     

模糊随机可靠性灵敏度分析方法

 依据模糊随机可靠性灵敏度的基本概念,提出了模糊随机可靠性灵敏度分析的通用数字模拟法,推导了模糊随机可靠性灵敏度数字模拟解的方差和变异系数的计算公式。当模糊变量的隶属函数为正态型时,可以将模糊随机可靠性灵敏度转化为随机可靠性灵敏度,并利用复合函数求导法则求解模糊随机失效概率对正态型隶属函数分布参数的灵敏度。对于工程中常用的对称三角型隶属函数,提出了两种近似等价正态化变换基础上的模糊随机可靠性灵敏度分析方法,这两种将对称三角型隶属函数近似等价为正态型的方法分别是“3σ规则”法和“最大最小”法。算例结果表明：由于“最大最小”法比“3σ规则”法所得的等价正态隶属函数能在形状上更好地近似对称三角型隶属函数,因而“最大最小”变换法更适用于对称三角型隶属函数情况下模糊随机可靠性灵敏度的近似计算。