分布参数液体管道的分段近似状态空间模型

分别提出了分布参数层流和湍流液体管道的分段近似状态空间模型,模型考虑了稳态摩擦和频率相关摩擦损失两种情况,层流管道的分段近似状态空间模型是线性的,而湍流管道的是非线性的。包括频率相关摩擦损失的模型比只考虑稳态摩擦损失的模型更精确地描述流体管道的瞬变过程,但模型的阶次成倍地增加了。用单根管道阀门关断的水击问题仿真计算对模型进行了验证,采用２个分段的分段近似状态空间模型的计算结果已经与特征线方法的计算结果吻合较好。     

航天器推进系统管路充填过程动态特性(1)理论模型与仿真结果

研究了航天器推进系统管路充填过程动态模型与计算方法。以一维管道瞬变流理论为基础 ,采用瞬变管流湍流频率相关摩擦损失模型计算管壁摩擦流阻 ,构造了管路分支流阻计算模型 ,考虑了液流撞壁后的蒸汽泡团释出与崩溃现象。数值计算采用特征线方法 ,仿真结果表明计算模型较好地描述了充填过程中的水击和管流振荡现象     

航天器推进系统管路充填过程动态特性(Ⅰ)理论模型与仿真结果

研究了航天器推进系统管路充填过程动态模型与计算方法。以一维管道瞬变流理论为基础,采用瞬变管流湍流频率相关摩擦损失模型计算管壁摩擦流阻,构造了管路分支流阻计算模型,孝虑了液流撞壁后的蒸汽泡团释出与崩溃现象。数值计算采用特征线方法,仿真结果表明计算模型较好地描述了充填过程中的水击和管流振荡现象。     

气固两相流k-ε-kp-εp双流体模型及应用

以流体ｋ－ε双方程模型为基础，提出了一种包括颗粒湍动能输运方程和耗散率方程的ｋ－ε－ｋｐ－εｐ两相湍流模型。并对两个垂直上升管道内气固两相湍流进行了数值计算，计算结果和已有实验结果吻合得很好。     

一种新的非线性K-ε两方程湍流模型

本文基于理论分析和对标准Ｋ－ε两方程湍流模型数值模拟结果的分析,在标准Ｋ－ε两方程模型的雷诺应力表达式中加入６高阶非线性项;并在模型的涡粘性系数μ１中考虑流动各向异性的影响,将涡粘性系数中参数Ｃμ取为表征流动各向异性的参数Ａ的一个线性函数,从而得到了一个新的非线性Ｋ－ε两方程湍流模型。     

某型航空发动机环形燃烧室振动特性分析

建立了环形燃烧室薄壁多孔结构的有限元模型,模拟燃烧室在实际工作中的安装条件,分析了燃烧室在不同工作状态下的模态频率和模态振型等.利用常温自由状态条件下有限元模型计算的振动特性,与该状态下的试验结果进行对比,修正计算模型.结果表明,修正后的计算结果与试验结果基本相符,保证了环形燃烧室各工作状态下振动特性分析的有效性.研究了不同工作状态对燃烧室温度场的影响及温度场对燃烧室模态频率的影响,结果表明燃烧室的工作温度越高,其模态频率越低;另外,计算结果表明该近似对称结构系统具有非单构系统的特点,其相似振型可作为同一阶振型加以考虑.     

卫星推进系统真空充填过程数值仿真

以一维管道瞬变流理论为基础, 利用特征线有限差分法, 建立空间推进系统部件流体动力学模型和液头充填的计算模型, 瞬变管流运动方程摩擦阻力采用平均摩擦损失模型和动态摩擦损失模型, 考虑溶入气体后管流声速的变化, 对某卫星推进系统管路充填过程进行计算.与实验结果对比表明:计算模型较好地描述了发动机充填过程中的管路水击和管流振荡现象, 可为其它推进系统管路充填流动瞬变特性分析提供指导.      

基于鬃毛摩擦理论的环形橡胶减振器动态迟滞特性模型

对环形橡胶减振器的动态特性进行了理论建模。综合考虑材料的高弹特性、频率相关动态特性以及摩擦相关动态特性,并将形状因素的影响考虑在内,提出了一种基于鬃毛摩擦理论的描述减振器动态迟滞特性模型。比较以往的基于库仑摩擦理论的模型,该模型具有对摩擦速度依赖特性的描述能力,且能够预测摩擦应力在零速附近的非线性特征,使得拟合出的动态迟滞特性曲线具有很多微小的波浪特征。开展了测量减振器径向动态迟滞特性的试验,试验结果表明:环形橡胶减振器的实际动态迟滞曲线并非光滑的椭圆环,而是具有很多细微波浪的近似椭圆环,且在激振频率偏高时,波浪特征相对明显,验证了基于鬃毛摩擦理论的理论模型与实际情况的规律和趋势一致、特征相似。      

连续突变管道流动局部水头损失建模研究

管道系统中设置的阀门、限流环等十分常见,装备限流装置的管道本质就是典型的连续突变管道。突变管道流动存在十分明显的局部水头损失,因此工程设计人员在进行管道系统设计时不得不考虑因管道突变而产生的局部水头损失;目前针对单个突变管道流动局部水头损失已有较为深入的研究,并已推导出相应的计算模型,但对于连续突变管路的水头损失研究较少,工程设计中对于连续突变管路的局部水头损失计算采用基于单个突变管道水头损失的叠加计算法,由于连续突变管路流动为极其复杂的湍流,基于单个突变管道水头损失的叠加计算法而得到的连续突变管路的局部水头损失结果误差较大。鉴于此,采用理论推导与实验相结合的方法,设计四种不同管径比的管道以及23种不同进口流量工况,并以航空煤油为介质进行实验,利用实验结果进行分析并对传统的突变管道水头损失的叠加计算法进行修正,提出了连续突缩突扩管道的局部水头损失修正计算模型,其计算结果相比于传统计算方法在精度上得到了显著性提高,为工程管道系统设计中的连续突变管路局部水头损失计算提供了有力依据。     

运用低雷诺数非线性k-ε模型研究湍流分离流动

给出了一个低雷诺数非线性ｋ－ε湍流模型，其中采用了一个非线性雷诺应力表达式并引入了两个考虑湍流近壁效应的修正函数。应用该模型和标准ｋ－ε模型计算了两种不同膨胀比的后台阶分离流。与实测结果比较表明，低雷诺数非线性模型结果较标准ｋ－ε模型结果有明显改进，它对流场的各特征量包括分离—再附长度、时均速度及各雷诺应力分量均给出了较好的预测。     

旋转坐标系中可压缩湍流的基本方程和湍流模型

本文在可压缩湍流瞬时量基本方程的基础上,用质量加权平均的概念,推导了旋转坐标系中可压湍流流动的连续方程、动量方程和能量方程。进而导出了旋转坐标系中可压湍流的应力输运方程、脉动平均动能方程和湍流耗散率输运方程。这些方程中包括了由旋转引起的 Coriolis力和离心力的作用项。这些方程是在旋转坐标系中研究平均 N-S方程求解湍流流场和研究湍流模型的基础。      

标准的与非线性的k-ε模型在二维剪切流中的应用

 采用标准 k-ε和非线性 k-ε模型对二维直通道及小曲率弯道中的剪切流进行了数值计算。与实验结果对比后得到 :由于非线性 k-ε模型考虑了流体变形与应力关系中的非线性项,因而能正确地反映湍流各向异性,能正确预测二维剪切流中湍流强度及湍流动能等参数的分布     

广义奇异摄动降阶方法在控制器降阶中的应用

 首先推导了广义奇异摄动降阶和奇异摄动降阶的 Schur方法,解决了均衡变换可能出现病态的问题,并且算法具有较好的鲁棒性。然后,将其应用于控制器降阶,并同其它方法比较,认为原系统带宽附近的控制器降阶误差对降阶效果影响较大,而广义奇异摄动降阶法通过参数α的选择,能取得较好的降阶效果。可供工程设计参考。     

密集模态组间非线性相互作用的初步研究

基于平均 Lagrange函数与调制项关系的 ALFM法,研究了通过 2∶ 1内共振联系的密集模态组间复杂的非线性相互偶合作用,深入探讨了密频内共振与 2∶ 1内共振的关系,提出了密集模态组在低阶主共振条件下相互作用及整体表现的初步规律,对于具有密集模态的大型柔性航空航天结构的振动主动控制和模态分析研究具有重要的指导意义     

Identification of a 2D turbulent wind spectrum

Turbulent wind spectrum models usually involve nonrational terms. For such models the development of a Markovian time simulator relies on a rational approximation obtained by the way of an identification stage. This paper presents a general method to manage this identification stage, that provides us with a family of stable rational approximations which is proved to converge towards the true model as the dimension increases. We first give an exact but infinite dimensional state space representation of the spectrum. It is based on the use of a diffusive equation. For bi-dimensional (2D) signals it also uses decoupling by spatial Fourier transformation. Then the discretization of this exact model leads easily to stable finite dimensional approximations over a prescribed frequency range. This approach is applied to the identification of a theoretical 2D turbulent wind spectrum, and of a 1D turbulent wind spectrum estimated from in flight recorded data.     

空-空导弹攻击区的高精度快速拟合

针对一般拟合法在处理空－空导弹攻击区这一复杂的非线性函数时所存在的问题,提出多关系式拟合-插值法,降低了拟合阶次。数字仿真结果表明,该方法极大地提高了拟合的精度和速度。     

大涡PIV方法测量水平矩形槽道内湍流耗散率

湍能耗散率的准确测量对工程实际和湍流的理论研究都有着重要的意义。以往的研究多是基于单点速度测量来估算单点的湍能耗散率或者有限体积内的平均值，而不能提供整个流场范围内的耗散率的分布情况。PIV能够测量瞬时流场的速度分布情况，它更适合于用来测量整个流场内的耗散率分布。本文分析了PIV测量得到的速度与单点测量方法之间的关系，提出了大涡PIV方法，运用大涡模拟中的亚网格尺度(SCS)应力模型，对PIV测量结果进行处理得到了水平槽道内的耗散率的分布。对采用不同SCS模型的结果以及量纲分析的结果进行了比较。     

用于结构优化的一种新的两点近似函数

提出一种基于两点信息的近似函数,它在当前设计点能够保证原函数的函数值、一阶导数值和近似二阶导数值;在另一设计点保证原函数的一阶导数值。用该近似函数形成结构优化原问题的序列近似问题,再用约束变尺度法解近似问题。通过典型算例的计算,结果表明本方法具有较好的收敛性和较高的计算效率,适用于结构优化问题的求解     

基于Reynolds平均的高超声速二方程湍流模型

 在Reynolds平均的框架下推导了可压缩湍流Reynolds应力方程和湍动能方程。根据一定的假设和尺度分析简化并封闭了所推导的湍动能方程。在均匀湍流假设下,湍动能耗散率可分解成可压缩性耗散和旋度耗散,并对其中的可压缩性耗散进行了封闭;同时认为旋度耗散不受可压缩性影响,直接引用传统的Reynolds平均不可压缩湍动能耗散率模型方程。由此构造了适用于高马赫数的二方程可压缩湍流模型。应用所发展的模型计算了高超声速平板绕流,并与若干现有模型的计算、实验与半经验公式的计算结果进行了对比,验证了所发展的模型。在此基础上,通过对压缩拐角的高超声速湍流的数值模拟,对所发展的湍流模型,以及若干现有模型进行了对比,显示了不同湍流模型及可压缩性修正在计算壁面压力分布和热流分布上的特点,说明了湍流模型可压缩修正的必要性,得到了所发展模型的计算结果最接近实验结果的结论。     

一种加速的参数化模型降阶方法

参数化模型降阶(PMOR)方法离线阶段训练数据的过程比较耗时,因此提出了一种加速方法。在基于矩阵插值的PMOR方法的基础上,采用组合近似的重分析技术对基于振型向量的模型降阶(MOR)过程进行加速,利用初始刚度分解矩阵生成迭代计算基向量以对系统矩阵降阶,通过降阶矩阵生成参数化模型的振型向量,对参数空间上的采样点重复整个加速计算过程生成离线数据库。并以电磁振动台动圈为例,采用普通方法和加速方法对均布采样样本点展开仿真研究,结果表明,在保证所构建的离线模型数据库准确度的前提下,此方法能减少80%以上的MOR计算时间,且随着采样点的增多,增速越明显,可以大幅度提高参数化降阶模型离线训练效率。