桥梁断面静风荷载的格子Boltzmann方法数值计算

通过引入反映湍流涡粘性的湍流松弛,得到了模拟高雷诺数湍流的BGK方程.在速度相空间、物理空间和时间上对BGK方程进行离散得到了三维十九速离散速度模型;结合分区计算技术,设计了格子Boltzmann并行算法;根据亚格子Smagorinsky模型,提出了直接从粒子分布函数计算湍流松弛时间的方法.用开发的并行计算程序对分体双箱截面和闭口箱梁截面的静风荷载进行了数值识别,得到的静力三分力系数和流场压力分布与风洞试验结果及CFD宏观方法计算结果吻合,并从表面压力分布入手分析了两种桥梁截面的绕流特点.     

耦合Navier-Stokes方程的可压缩边界层数值求解

为了在Reynolds-averaged Navier-Stokes(RANS)方程计算中耦合eN方法进行转捩判断,在RANS方程求解过程中耦合求解了可压缩层流边界层方程为判断转捩提供了精确、可靠的边界层信息.利用等熵关系由RANS方程求出的物面压力分布确定边界层外边界的速度分布,进一步确定出边界层外边界.边界层方程的求解使用Keller提出的二阶BOX方法.为了验证方法求解边界层方程的正确性,在低速流动状态下将计算结果和XFOIL的计算得到的边界层解进行了比较,二者吻合较好.     

有限空间内低旋流流动与燃烧特性

为了对低旋流流动特性及其对火焰传播与稳定影响进行深入理解,通过实验与数值模拟相结合的方法,对有限空间内低旋流流动与燃烧特性进行研究.结果表明:改变中心孔板小孔直径,可以形成不同旋流数的旋流喷射器;在冷态条件下,随着旋流数增加,在展向平面上逐渐形成“W”形速度矢量分布;在铅垂截面内,速度矢量分布存在一个环形高速区域、中心低速区和外侧低速区.由于燃烧气体吸热膨胀、涡量诱导和温度斜压效果,使得冷态高旋回流流动转变为了燃烧条件下低旋无回流流动;在有限空间内,低旋流燃烧没有形成脱体火焰,而是形成类似于扩散燃烧火焰;随着旋流数增大,温度分布边界和化学反应边界在展向平面逐渐缩短,而在铅垂截面逐渐扩大.      

曲率和挠率对圆截面螺旋管道粘性流动的三阶作用

利用Ｇｅｒｍａｎｏ的曲线正交坐标系下的完全Ｎ－Ｓ方程,采用摄动方法求解了圆截面螺旋管道内的粘性流动,得到完全三阶摄动解,研究了曲率和挠率对二次流动和轴向速度的一阶、二阶和三阶作用。在曲率一阶项的作用下,截面上形成方面相反的两个二次涡,轴向速度最大值移向外侧;曲率的二阶项使二次涡的中心移向外侧;在曲率和挠率的高阶项（二阶和三阶）的共同作用下,靠近下半截面的涡增加,同时,二次涡的中心和轴向速度最大值移     

旋转通道一维非定常计算模型

采用特征线法,对航空发动机二次空气系统的旋转通道结构一维非定常流动及换热的计算模型进行了研究,通过将广义的定截面、有摩擦、考虑对外热交换的旋转通道可压缩非定常流动偏微分方程组转换为常微分方程组,考虑旋转离心力,并采用Colebrook White方程定义其阻力特性,进而获得旋转通道内部沿流动方向压力、质量流量、温度等系统参数的特征线数值解,采用商业软件Fluent对其进行验证.结果表明:特征线法利用压力波传播方式对旋转通道进行求解,能够直观地描述旋转通道一维非定常流动及换热特性,计算误差不超过5%,可以作为二次空气系统旋转通道一维非定常计算研究的有效手段.      

基于Kriging模型的高超声速舵面优化设计

采用RANS方程进行了高超声速舵面的气动特性计算与平面形状优化设计研究。为了提高优化效率,采用拉丁超立方抽样试验设计方法得到样本点,并建立Kriging模型替代费时的流动数值模拟。以最大化舵面操纵力矩和最小化平面面积为目标,以阻力不增为约束条件,进行了舵面平面形状优化设计。优化结果表明,优化舵面的平面面积减少了4.78%,力矩绝对值在设计点增加了24.69%,且在较大的攻角范围均有提高。将初始舵面和优化舵面分别装配到方形截面导弹上,验证计算表明,优化舵面的操纵效率更高。     

迷宫密封的理论与实验研究

 根据相似原理。通过实验回归分析法,得出压力损失系数ζ与齿腔几何特性和流动特性的关系式。然后应用能量定律、连续方程、状态方程和Bernoulli方程与其联立,求解后可获得迷宫密封的全部特性。并可归纳成流量系数α、穿越系数β和迷宫函数ψ对全压比λ的特性曲线。应用这类曲线,可方便地计算出迷宫密封的泄漏率。     

粘性跨音速喷管流场计算

本文继用隐式近似因子分解法成功地计算无粘跨音速喷管流场之后,用同样的方法,结合任意曲线坐标系,通过非定常方程在相当长时间后的定常解,针对矩形截面喷管和轴对称喷管粘性跨音速流场,求解了可压缩层流薄层N-S方程,获得了核心流,特别是边界层中的流动参数.对于矩形截面和轴对称两种喷管进行了计算,其结果和实验数据相当一致.将本文的方法应用到粘性两相流动计算,可望较多的节省机时.     

高负荷跨音涡轮导向器环形叶栅试验研究

对一个高负荷跨音速涡轮导向器进行了环形叶栅吹风试验 ,对其根、中截面进行了平面叶栅吹风试验。本文对试验的损失特性、出气角特性、速度特性、叶型表面速度分布特性等作了对比分析。认为导向器环形叶栅试验结果与平面叶栅试验结果有可比性 ,其叶中截面相当一致 ,根、尖截面则由于导向器流动的复杂性而与平面叶栅试验结果有较大差异     

过渡式空腔流动数值计算

采用雷诺平均N-S方程数值模拟了长深比L/D=12的空腔流动,分别计算了0.85Ma和1.5Ma两种速度下的流场,得到了空腔内压力、速度、涡量等分布图,分析了空腔内部流动特性。从计算结果可以看出,空腔流动明显表现出过渡式空腔流动特性,空腔前缘和后缘都有明显的涡流区域,而且当来流速度不同时,空腔内部流动分别表现出了闭式-过渡式和开式-过渡式的空腔流动特性。     

非线性数值模拟柱面气膜密封动态特性

针对气膜密封研究中密封气膜动态特性分析难点问题,采用非线性数值模拟法对柱面气膜密封系统的动态特性进行研究[D].建立了柱面气膜密封系统的动力学模型,给出了基于有限元数值方法的瞬态气体雷诺方程（含时间项）与系统动力学方程的耦合求解流程和方法,数值仿真给出了柱面密封环的中心轨迹图、确定了一定转速下的临界无量纲质量、不同密封环无量纲质量时等量位移扰动下系统的响应图,分析了密封环无量纲质量对系统稳定性的影响.结果表明:光滑柱面气膜密封系统存在临界无量纲质量,临界时系统受扰后处于半频涡动状态,随着密封环无量纲质量的增加,受扰后系统收敛时间变长,稳定性变差.      

有限元的OpenMp并行计算技术

传统有限元串行计算技术耗时长、效率低,已远远不能满足工程实践需要。通过分析串行计算的耗时分布、程序可并行性及数据相关性,在保持串行程序不变的基础上,采用基于多线程的OpenMP技术实现单元刚度矩阵形成总体刚度矩阵的程序并行化。针对有限元计算机翼实例进行了不同线程数目下并行加速比的研究,结果表明,OpenMP技术可有效提高有限元计算效率,且存在一个最佳的线程数,使得并行加速比达到最大,此外加速比还与问题的计算量等因素有关。     

特殊形状橡胶密封圈的性能分析

采用非线性有限元法对特殊截面形状的橡胶密封圈进行了分析.通过轴对称有限元分析模型, 分析了其在矩形沟槽内装配压缩过程中的变形和应力变化, 以及应力的分布情况, 同时还对承受单侧压力的密封情况进行了研究.对于密封圈的抗扭转特性以及在特殊压力负载情况下的密封特性也进行了模拟并且和普通O形圈进行了对比.结果表明, 所采用的方法能够预测特殊形状的密封圈在压缩中的变形和应力等特征参数, 加深了对其密封性能的了解, 并对同类密封结构设计有一定的指导意义.      

篦齿封严装置流动诱导力的计算

本文以文献[1、2]的模型和方法为基础,对篦齿封严装置的流动诱导力作了一些研究导出有些不同的线化公式,算例结果也有所不同。      

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

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

航空发动机主泵金属封严装置漏油故障分析

针对航空发动机主燃油泵金属封严装置漏油问题,分析了导致漏油产生的各种因素;研究了金属封严圈组合研喜及检测、铅铟润滑扩散镀昙工艺控制和主燃油泵装配工艺3个方面对封严的影响,提出了针对性的改进措施。进行了相应的试验．验证了改进措施的有效性。     

动静压混合式气体密封追随性及主动调控振动特性数值分析

动静压混合式气体密封(DHHGS)运行中,静环对动环轴向窜动和角向摆动的跟踪响应(追随性)能够降低外界干扰对密封稳定性造成的影响。将DHHGS简化为弹簧-阻尼-质量系统,基于摄动法求解了DHHGS的动态非线性Reynolds方程,得到了表征密封动态特性的气膜刚度和阻尼系数。研究了两种DHHGS (泵入式和泵出式)在3个方向简谐激励作用下的追随性,并得到了阻封气压力对追随性的影响规律,分析了主动调控时静环的轴向自振稳定性,并给出了轴向自振的临界失稳判据。研究结果表明:即使激励振幅足够大,DHHGS仍具有良好的追随性;阻封气压力增大,密封的追随性增强;主动调控时,静环自由振动是非往复的衰减运动,密封仍能稳定运转。     

基于电液伺服控制燃油泵的航空发动机控制与仿真

为探索轻质化燃油系统结构,基于电调燃油变量泵的航空发动机转速控制系统,构建了柱塞泵斜盘位置电液伺服控制系统,油泵出口燃油直接输入电液伺服阀;建立了电液伺服阀线性化模型。通过数字仿真,研究了电液伺服阀工作特性,并得到了其适应性模型;在航空发动机特性半物理试验系统上,对斜盘位置电液伺服控制系统实物进行了验证试验,并与航空发动机模型一起构成了发动机转速闭环控制系统。结果表明：变输入压力的燃油电液伺服位置控制系统有效可行,变量泵工作稳定可靠,电液伺服阀模型能够准确反映实际工作状况;基于变参数PI控制算法的转速闭环控制初步取得成效。     

Analysis of crack opening stresses for center- and edge-crack tension specimens

Accurate determination of crack opening stress is of central importance to fatigue crack growth analysis and life prediction based on the crack-closure model. This paper studies the crack opening behavior for center- and edge-crack tension specimens. It is found that the crack opening stress is affected by the crack tip element. By taking the crack tip element into account, a modified crack opening stress equation is given for the center-crack tension specimen. Crack surface displace- ment equations for an edge crack in a semi-infinite plate under remote uniform tension and partially distributed pressure are derived by using the weight function method. Based on these displacements, a crack opening stress equation for an edge crack in a semi-infinite plate under uniform tension has been developed. The study shows that the crack opening stress is geometry-dependent, and the weight function method provides an effective and reliable tool to deal with such geometry depen- dence.