旋成面叶栅气动反命题的新解法(Ⅱ)

本文利用映象平面ξ-η的特点, 将旋成面叶栅气动正命题的微分-积分解法与叶型修正的“喷气模型”结合, 建立了一种适用于旋成面叶栅气动反命题的新解法。叶片表面给定的气动参数分布作为反命题的边界条件, 在求解过程中直接得到满足。通过对一例实际透平叶型的计算验证了本文解法的可靠性和实用性。      

一种新的适用于叶栅绕流计算的网格生成法及其应用

本文对平面叶栅气动计算网格生成方法作了简要评述。在此基础上,根据叶栅槽道流动的特点,提出一种以等比级数为基础,用几何方法构作H网格的数值方法。用这种网格进行流场计算的结果表明,所生成的网格在叶片前后缘区有很高的分辨能力,能够算出速度分布的细微结构。本文的计算结果与亚音速及跨音速叶栅速度图解符合得很好。     

基于小生境遗传算法和RANS方程的平面叶栅气动优化设计

发展了一种将小生境遗传算法与RANS方程数值求解相结合的平面叶栅全局优化设计方法.该方法针对遗传算法的搜索原理和平面叶栅的气动外形特点, 提出了与之相适应的Bezier曲线参数化表达, 构造了用于优化设计的适应度函数.优化设计变量是平面叶栅的参数化Bezier曲线特征多边形控制点坐标, 目标函数为极大化叶栅的升阻比, 约束条件综合考虑了叶栅的进出口几何安装角及叶片强度和性能等要求, 约束条件的处理采用松紧罚函数法.优化得到的叶栅升阻比比初始叶栅提高了8.3%.      

应用嵌入式H型网格技术计算孤立翼型绕流场

提出一种嵌入式H型网格技术,用于生成孤立翼型和叶栅的计算网格。与传统的H型网格相比,采用嵌入式H型网格时,翼型头部及尾部附近的网格畸变很容易得到控制,有益于提高计算的准确性。      

航空发动机叶片颤振机理探究

航空燃气涡轮发动机中,由颤振引发的叶片疲劳失效是导致发动机机械故障的主要原因之一。为了研究叶片颤振这一非定常气动现象,在能量法的基础上采用弱耦合的计算方法,模拟了发动机风扇叶片在不同转速和不同出口反压条件下的振动,得到了叶片的颤振边界。其中流场与叶片建模用到三维线性插值法转化振动位移,流场动网格建模用到了多层动网格技术。计算表明,叶片表面气动功的分布与叶栅内的激波分布有很大的关系,激波在叶栅通道中的振荡是引起叶片颤振的原因之一。激波在叶栅中前后振荡,使得叶片表面上单位面积做正功的区域发生变化。做正功的区域扩大时,叶片表面气动功绝对值减小,从而气动阻尼比减小,更靠近颤振边界。     

涡轮叶片气动设计软件BladeDesign

涡轮叶片气动设计软件BladeDesign将涡轮气动设计中迭代最频繁的叶栅几何设计、s1流面计算、叶片积叠三个环节集成起来,极大地提高了涡轮部件设计的效率。叶栅型线设计采用目前流行的Bezier曲线,叶栅造型方法充分考虑了工程实际。集成的S1流面计算网格划分采用ANSYS ICEM CFD 11．0,分析采用ANSYS...     

涡轮平面叶栅非轴对称端壁优化设计

开发了一套造型灵活直观、网格生成速度快的涡轮平面叶栅非轴对称端壁优化设计工具,该工具的核心技术是非均匀有理B样条(NURBS)曲面造型和网格变形.在此基础上以商业软件Isight为优化驱动器,以CFX为求解器,搭建了非轴对称端壁优化设计流程.以Pack B涡轮平面叶栅为例,对其进行了非轴对称端壁优化设计.优化后涡轮平面叶栅总压损失系数减小了12.96%.结果表明:涡轮平面叶栅端部的静压分布改变削弱了涡轮平面叶栅通道中马蹄涡、通道涡的强度,提高了涡轮平面叶栅的气动性能.      

一种平面叶栅反设计方法的改进与验证

介绍了一种平面叶栅反设计方法,假设因压力作用有一流体微元沿槽道中心流线运动,设计参数包括:叶栅进/出口气流角、截面半径、叶型载荷分布和厚度分布等.为完善原设计方法,提出了新的理论,用于获得必要的密度分布曲线;改进了前、后倒圆造型算法,提升了叶栅气动性能;为验证改进后的设计方法,重新设计了德国Aachen涡轮静叶叶栅.      

任意旋成面叶栅的跨声速绕流计算

本文中的方法是把时间推进有限体积法和流线迭代法结合起来,用以解决任意旋成面叶栅跨声速绕流的计算问题。在旋成流面上,网格是由一族流线及一族固定的节距线所组成的。流线的位置预先并不知道,在计算过程中需要逐步加以修正。采用流线作为一族网格线可使网格线的分布更为合理,并使差分格式简化。同时,为了加速计算过程的收敛,在本计算中采用逐步加密网格技术。     

基压对跨声速涡轮叶栅尾迹损失影响的研究

跨声速涡轮尾迹损失是叶栅损失的主要部分,它大约占总损失的三分之一.跨声速尾迹气流是十分复杂的.必须要了解其基本的气流流动模型.目前计算损失的方法大多是根据经验的公式,但这个方法限制了计算损失的准确性.国外的一些研究表明:基压和损失,基压和反压都存在着一定的关系.本文就是利用超、跨声速平面叶栅风洞在近二十年中所做的叶栅试验数据,进行分类整理,寻找出基压对反压和基压对损失的简便的经验公式,为叶型设计的气动计算提供叶栅损失系数和叶片表面马赫数分布的预估.     

平面热应力问题的组合杂交有限元方法

利用组合变分原理可以增强杂交元方法解的稳定性,建立了平面热应力问题的基于区域分解的组合变分原理.与弹性力学问题相比,右端项改变;并且提出组合杂交元来离散新的变分原理、形成刚度矩阵;进而引入能量协调条件,不仅简化了变分原理及其相应的刚度矩阵,而且减少有限元解的误差.数值结果表明,能量协调的4节点组合杂交元数值性能最佳:实现粗网格高精度,对网格畸变不敏感且可以克服闭锁现象.     

一种求解全位势方程的跨音松弛法

 本文发展了一种适用于求解二维跨音速叶栅流场的松弛方法。采用流线和平行于y-轴的直线族组成的非正交网格系统,在物理平面中求解全位势方程,沿流线构造相关差分格式。数值结果表明,本方法具有收敛速度快、差分格式稳定的优点,松弛150步以后,便可获得基本满意的结果。计算结果和物理分析一致,和试验数据相近。     

前后缘型线同时可控的乘波体设计

提出了一种前后缘型线同时可控的乘波体设计方法,在马赫数可控的外锥形曲面基准流场中,结合流线追踪技术和混合函数,实现了椭圆前缘转椭圆后缘的乘波体设计,并在设计点(Ma=6.0)和接力点(Ma=4.0)对其进行数值仿真研究。前后缘同时可控的乘波体在型面剧烈过渡处产生了较弱的激波,出口两侧存在高温高压区,后部对称面附近的激波形状由圆弧变为平直线且出口处流场基本均匀,非常有利于与进气道匹配设计。另外,该乘波体具有较高的容积率和预压缩效率,附面层修正后的容积率为0.24,设计点时乘波特性较好,接力点时前部完全乘波,具有较高的升阻比,有黏条件下设计点和接力点的升阻比分别为2.54和2.41。此外,与给定前缘的乘波体相比,其升力、阻力、俯仰力矩和出口增压比都有明显增加,但是升阻比和出口总压恢复系数有所降低,在设计点无黏升阻比由3.56降为3.00。以上研究表明,本文的设计方法可行且更加灵活,拓宽了乘波体的选择范围。     

基于欧拉方程的二维振荡机翼非定常气动设计反命题方法

建立了基于欧拉方程组的二维振荡机翼非定常气动设计反命题方法的数学模型.通过一系列变换,将物理时空中的求解域转换成映象坐标系中的规范区域;并导出了映象坐标下的欧拉方程,结合非定常反命题的边界条件,便可用有限差分方法求解.本方法引入映象坐标系解决反命题边界形状的不确定性,并能利用欧拉方程现有的各种差分格式于反命题求解.     

一种求解涡轮平面叶栅的跨音松驰法

一、引　　言跨音松驰法是七十年代初发展的一种新方法。文献 [1 ]是一篇关于跨音松驰法的奠基性文章 ,文章提出了求解小扰动势方程的混合差分格式。文献 [2 ]发展了求解全位势方程的旋转混合差分格式。近年来 ,求解跨音全位势方程的方法有了更全面的发展 [4、5]。作者于 80年发展了一种求解大弯度二维叶栅的跨音松驰法[3] 。本文将该法推广到求解涡轮平面叶栅跨音带激波流场 ,直接求解全位势方程 ,进一步提高边界区差分格式的精度 ,并注意初场的给法。本法具有节约机时和内存的优点。数值算例表明 ,对于头部不大的涡轮平面叶栅 ,应用本方法…     

翼型—扰流片的分离气动特性计算

本文介绍了用涡面元法模拟带扰流片的翼型低速无粘分离绕流。在翼型和扰流片的面元上分布线性变化涡。在扰流片后的上下分离流线的面元上分布等强度的涡。上分离流线始自扰流片的梢部，下分离流线自翼型的后缘引出。分离所泡由两离散涡结尾。气泡内总压为常值，它与涡强大小一同求解。分离气泡的形状在迭代求解过程中确定。压强分布和升力系数的计算值与现存文献的数值结果和实验数据是一致的。     

翼型前缘变形对动态失速效应影响的数值计算

翼型或机翼的动态失速效应所引起的低头力矩和正气动阻尼限制了飞行器气动性能的提高,甚至可能诱导发生不稳定运动。应用于小尺寸机翼的前缘动态变形（DDLE）技术,通过实时改变前缘形状,能够改善翼型前缘区域的速度梯度,进而抑制动态失速效应。采用转捩剪切应力输运（SST）黏性模型结合分区混合动态网格技术,研究了这种前缘变形对机翼俯仰运动所引起的非定常流动的影响,得到通过小幅度前缘变形抑制和延迟动态失速的方法,从而提高翼型的气动性能。翼型NAC A0012的数值模拟结果与动态失速风洞试验结果比较表明：所使用的数值计算方法能够较为准确地模拟翼型在动态失速过程中升力系数与俯仰力矩系数的变化情况,可用于研究前缘变形对翼型俯仰运动所引起的非定常流动的影响。前缘动态变形翼型俯仰运动过程的非定常流场的数值模拟表明：在大迎角下不同幅度的前缘下垂运动能够抑制流动分离的发生,从而抑制动态失速,但在大迎角下小幅度高频率的前缘下垂变形能更高效地抑制动态失速;前缘变形幅度以及变形沿中弧线的分布对升力系数和俯仰力矩系数的影响并不明显。     

双锥Bump压缩面设计及气动特性

采用反设计方法研究双锥Bump压缩面设计技术,该方法实质是融合运动间断边界在运动网格条件下对轴对称欧拉方程组求解.为避免激波捕获法对激波型面位置求解误差,应用计算域分块的方法预估第2道激波,在此基础上采用流线追踪法生成Bump压缩面;结合平面机身设计一双锥Bump实例,运用计算流体力学仿真手段对其进行黏性数值模拟.研究结果表明:①该模型流场结构仍保持较强的附面层扫掠能力;②在来流马赫数为2.0条件下,相比传统的正圆锥乘波体Bump压缩面设计,新型设计方法可使Bump外压缩系统总压恢复系数提高0.04左右,为Bump进气道性能提高奠定基础.      

Topology optimization of compliant adaptive wing leading edge with composite materials

An approach for designing the compliant adaptive wing leading edge with composite material is proposed based on the topology optimization. Firstly, an equivalent constitutive relationship of laminated glass fiber reinforced epoxy composite plates has been built based on the symmetric laminated plate theory. Then, an optimization objective function of compliant adaptive wing leading edge was used to minimize the least square error（LSE） between deformed curve and desired aerodynamics shape. After that, the topology structures of wing leading edge of different glass fiber ply-orientations were obtained by using the solid isotropic material with penalization（SIMP） model and sensitivity filtering technique. The desired aerodynamics shape of compliant adaptive wing leading edge was obtained based on the proposed approach. The topology structures of wing leading edge depend on the glass fiber ply-orientation. Finally, the corresponding morphing experiment of compliant wing leading edge with composite materials was implemented, which verified the morphing capability of topology structure and illustrated the feasibility for designing compliant wing leading edge. The present paper lays the basis of ply-orientation optimization for compliant adaptive wing leading edge in unmanned aerial vehicle（UAV） field.     

翼身融合飞机参数化几何模型

本文研究翼身融合飞机的参数化几何模型生成的方法。翼身融合飞机各部件均视作翼面类部件,运用形状函数和分类函数变换方法描述飞机各个剖面翼型。用三次施密特曲线描述飞机沿展向的厚度变化与后缘曲线形状,用分段曲线描述复杂的前缘形状。翼梢小翼划分为过渡段与主段两部分分别描述,并通过变换矩阵计算主翼段关键参数坐标。应用CATIA二次开发技术实现参数化几何模型的自动生成。