直升机用新型发动机燃烧性能的数值研究

针对液态辛烷/空气两相混合物燃烧波在一维封闭空间内传播的爆燃向爆震转变(DDT)问题,运用化学反应流动基本理论、颗粒轨道模型和两步反应模型,在拉格朗日质量坐标系中建立了对该现象的数学表达,并开发了数值模拟程序CTPD,其中气体动力学方程采用显式Sγα,β差分格式,化学反应项使用二阶精度的Adams方法.通过数值模拟研究了液态燃料爆震波的基本特性,以及点火能量及初始温度对于爆震波结构与发展的影响作用.研究发现,点火能量及初始温度的提高显著促进了爆震的转变,爆震强度随着点火能量的提高而增大,随着初始温度的提高而减小.研究中计算结果与实验波形符合良好,说明本文为包括DDT现象在内的脉冲爆震发动机工作过程提供了一种可行的计算方法.     

基于内禀模态高频积分能量法的煤油发动机爆震因子计算方法

提出了基于内禀模态高频积分能量法的爆震因子计算方法.该方法对电控二冲程煤油发动机缸内燃烧压力信号进行经验模态分解,自适应得到若干内禀模态函数分量与残余函数.对内禀模态函数分量采用快速傅里叶变换得到爆震高频信号分量,选取残余函数峰值对应的曲轴转角为爆震窗口起始曲轴转角,对爆震窗口宽度分析,得到合理爆震窗口宽度持续曲轴转角为30°.利用内禀模态高频积分能量法对无爆震、轻微爆震、中度爆震与强烈爆震工作循环的缸内燃烧压力信号进行爆震因子计算,得到4个工作循环的爆震因子分别为0.6642,1.8191,3.0275,5.3717,可表征爆震的强弱.研究表明基于内禀模态高频积分能量法的爆震因子计算方法简便、快速有效.      

相位延迟边界条件在叶轮机械颤振分析中的应用

基于带相位延迟的周期边界条件,建立了某跨声速转子的双通道高效气动阻尼计算模型.数值计算了该转子的气动性能、颤振边界和叶片模态,和实验数据吻合较好.通过传统的多通道能量法以及双通道方法计算了叶片在一弯模态,不同叶片间相位角条件下的气动阻尼,获得了基本一致的计算结果,而双通道方法相比于传统的多通道能量法计算效率提升约7.7倍,内存需求约为后者的0.45倍.不同叶片振幅对气动阻尼结果的影响研究表明,对于较小的叶片振幅,流动非线性对气动阻尼计算结果仍然有显著的影响.不同工况的计算结果表明:叶片间相位角对转子叶片的气动阻尼有显著的影响,对于该转子最小的气动阻尼均在叶片间相位角为-42.4°时得到;同时,在近颤振状态,不同叶片间相位角对应的气动阻尼均小于近设计状态.     

基于机翼-帆尾的高纬度跨年驻留太阳能飞机总体参数设计方法

拓展太阳能飞机在较高纬度地区的跨年驻留性能有助于促成太阳能飞机的广泛实用化。建立了适用于任意高度、任意纬度、任意指向的光伏组件面功率模型,并考虑了光伏组件的温度效应,通过能量仿真得出：在方位角跟踪方式下,滚偏角为90°的主动式光伏组件的日均面功率最优。然后在布局与能源综合设计思想指导下,建立了一套基于机翼-帆尾的太阳能飞机总体参数设计方法,其组成模块包括各部件质量方程、气动效率方程、用于构建气动布局参数与全机光伏组件面功率特性之间映射关系的Kriging代理模型,以及参与总体参数匹配优化设计的量子粒子群优化（QPSO）算法及其多目标评价函数。面向高纬度与跨年驻留的设计指标,开展了机翼-帆尾太阳能飞机的方案实例设计,其中驻留纬度与高度指标分别为45°N和18 km。详细分析了此方案在23.5°N~55°N纬度域内的可持续高度包线。研究结果表明：与传统布局形式相比,机翼-帆尾布局形式大幅提升了高纬度地区冬季附近的光伏组件面功率,有效地减小了翼展尺度、机翼面积并提升了巡航速度,具有良好的应用优势。方案设计实例也验证了基于机翼-帆尾的太阳能飞机总体参数设计方法的可行性。     

血管壁内应力分布研究

应力一直被认为是影响动脉血管平滑肌细胞组织性质变化的重要因素.利用零应力状态和张开角的概念,推导受内压的血管壁的内力分布的解析解,并对具体的实例进行计算分析.计算结果表明,张开角对于沿血管壁厚度方向变化的圆周方向应力分布具有非常强烈的影响.      

氢频标贮存泡涂敷材料特性及工艺研究

氢原子频标中贮存泡内壁的成膜状态直接决定跃迁信号的产生及指标,贮存泡的涂敷和烧结工艺成为氢原子频标研制中的一项关键技术。通过分析聚四氟乙烯和聚全氟乙丙稀树脂的热处理特性以及氢原子频标对成膜的要求,分别针对蓝宝石氢原子频标和被动型氢原子频标贮存泡的特点,给出了相应的涂敷和烧结工艺。实际氢原子频标系统的测试表明,该涂膜工艺满足氢原子频标的使用要求。     

硝酸酯增塑高能推进剂爆炸性能研究

对我国硝酸酯增塑高能推进剂的爆炸性能、爆炸当量等进行了试验和评估。简述了共试验装置、方案、试样、试验过程及结果。结果表明：该推进剂危险性分类与美国国防部划入１．１级危险品的高能推进剂相似；ＴＮＴ当量１０ｍ之内为１．０７～１．２６，１０～５５ｍ之间约为０．８９～１．２１，也与国外相关试验值相似；另外观测到１０ｋｇ该推进剂爆炸时，距３５ｍ以外建筑物无大影响，５５ｍ外无影响。     

叶间相位角对叶片颤振的影响

发展了一种基于能量法并计及叶间相位角(IBPA)影响的数值计算方法.建立了全环振荡叶栅模型,将每一个扇区分为叶片表面附近的可动域及其外围的固定域,采用有限元形函数将计算结构动力学(CSD)网格点位移的传递到计算流体动力学(CFD)的耦合面网格点上,求解了全环振荡叶栅由k-ε湍流模型封闭的Reynolds平均Navier-Stokes(RANS)方程.以NASA 67转子叶片为例,在设计转速下,基于能量法分析了叶间相位角对叶片颤振的影响.计算结果表明:叶间相位角对叶片颤振有显著影响.      

复合材料层合板低速冲击损伤方法研究

在试验的基础上,分析了冲击能量与凹坑深度之间的关系,拟合出了冲击能量与凹坑深度曲线方程。曲线方程表明,凹坑深度的变化是与冲击能量的变化过程相适应的,从而在凹坑深度已知的情况下可以计算得到层合板所受的外载,有此外载从便对层合板进行相应的模拟。利用AN-SYS有限元程序对复合材料层合板横向低速冲击进行了模拟,模拟采用瞬态分析方法。对冲击后的试件进行了C扫描,在此基础上对损伤的分布形式及大小做了详细的分析。计算结果和试验结果表明,方法是可行的,拟合曲线是正确的,且适合于层合板冲击后的损伤评估。     

结构动力学方程的显式与隐式数值计算

从理论上详细推导了动力学方程的显式隐式时间积分数值计算方法,针对经典中心差分格式在求解阻尼结构计算效率降低的问题,研究提出了一种新的显式积分格式,理论证明该数值积分格式能够处理各种复杂阻尼状态而不降低计算效率。对于不同的时间步长通过设定调节加速度参数,能得到比经典显式计算格式和NEWMARK逐步时间积分更高的精度,最后通过数值算例对三种计算格式的理论分析进行了验证。