基于不可逆热力学的宏细观粘弹性损伤本构模型

为建立复合固体推进剂的损伤本构模型,基于不可逆热力学和叠加原理,通过引入宏观损伤效应张量,推导出一个通过有效应力表征损伤的蠕变型损伤本构方程。假设材料为初始各向同性,进一步引入细观标量损伤效应函数表征材料对称性的改变,进而得到由细观损伤效应函数表征宏观损伤效应张量的一般表达式。通过选取合适的细观损伤效应函数,文中建立的本构方程可以用于表征材料的正交各向异性损伤、横观各向同性损伤和各向同性损伤。最后,基于Schapery粘弹性微裂纹扩展模型,选取相对微裂纹密度为损伤内变量,建立了一个由损伤热力学对偶力表征的幂律型损伤演化方程。数值结果表明,建立的模型能够较好地反映材料损伤的率相关性和温度依赖性,具有良好的预测精度。     

旋转弹体马格努斯效应数值模拟方法研究

旋转的弹体通常会存在马格努斯效应,严重影响飞行稳定性和弹道轨迹。提出基于旋转壁面法的多对称面弹体旋转效应数值模拟方法,无需借助动网格和多参考系模型,仅需要将旋转速度附加壁面,即可使用定常方法模拟马格努斯效应;同时,使用SOCBT标准弹体模型对该方法进行验证分析。结果表明:计算值与试验值吻合较好;由于旋转引起的附面层堆积以及大攻角分离效应、激波的干扰,计算时附面层的网格节点要足够密,并且应该针对不同的攻角选择精度与效率兼顾的计算模型。     

局部催化特性差异对气动热环境影响的计算分析

高温气体非平衡效应及其壁面催化效应对高超声速飞行器气动热环境造成显著影响,是当前高超声速飞行器气动热环境预测和热防护设计的关键问题之一。考虑高温空气离解与电离等化学反应、气体分子热力学激发、流动中的非平衡效应和壁面催化效应,通过数值求解三维热化学非平衡Navier-Stokes方程和壁面处质量、能量平衡关系,完善了高温气体热化学非平衡流场有限催化气动热环境数值计算方法和计算程序,采用典型算例进行了考核验证。在此基础上,开展了不同条件下高超声速飞行器热化学非平衡流场气动热环境数值模拟,分析局部催化特性差异对气动热环境的影响。研究表明：所建立的高超声速飞行器热化学非平衡流场有限催化气动热环境数值计算方法及程序,其数值模拟结果与飞行试验、文献符合;局部催化特性差异会导致热流跳变,其热流跳变量与催化特性差异量、材料分布方式等有关;催化特性差异较大时,局部区域热流可能远远高于飞行器全表面完全催化的热流结果,此时将飞行器在全表面完全催化（FCW）和完全非催化（NCW）条件下的数值模拟结果作为实际飞行过程中表面热流的上、下限这一简化处理方式,是不可取的。     

剪切稀化非牛顿射流撞击液膜破碎直接数值模拟

液态射流撞击是液体火箭推进系统中广泛采用的一种燃料雾化方法,其破碎特征直接影响燃料最终的掺混及燃烧效率。采用直接数值模拟（DNS）工具,研究了低雷诺数（Re_l=41）和中等韦伯数（We_l=163）条件下剪切稀化非牛顿射流撞击液膜破碎的问题,着重分析了对角液膜的三维结构、破碎特征和非牛顿特性等。研究结果表明：在所研究的射流参数下,该非牛顿撞击液膜破碎属于Open Rim类别,破碎过程具有三维特性并伴随液丝与边缘的融合、液丝向液滴的转变等时域流动特征。液体的总表面积随时间不断增长,但单位表面积随液膜破碎的发生而下降,液膜扩张半角随时间逐渐增加并趋于恒定值43°,而后部液膜的长度几乎不随时间发生变化。此外,撞击液膜表现出明显的剪切稀化特性,液体内部最低黏性系数仅为零剪切黏性系数的1/5。     

多电极偶接对金属大气腐蚀影响的试验与仿真

搭建薄液膜厚度测量与控制装置,采用微距参比电极后置法组建三电极体系,测量了厚度为100 μm的3.5wt% NaCl液膜下2A12铝合金、TC18钛合金和30CrMnSiNi2A高强钢的极化曲线和耦合体系的电偶电流;建立了基于薄壳电流分布的Comsol腐蚀仿真模型,得到了100 μm液膜下电极电位分布和电流密度分布图,通过对电极表面进行局部电流面积分计算得到了其电偶电流,对不同电极进行了面积的参数化扫描,讨论分析了不同面积比对于电偶电流的影响。结果表明：利用仿真模型得到的电偶电流值与试验值吻合较好,在2A12/30CrMnSiNi2A/TC18多电极体系中2A12充当阳极,TC18/30CrMnSiNi2A充当阴极,且阴阳极极性不随面积比的变化而发生转变,电偶电流主要集中在30CrMnSiNi2A与2A12之间,电偶电流与TC18面积呈线性关系,其对数与2A12和30CrMnSiNi2A面积的对数之间呈线性关系。利用建立的仿真模型对模拟搭接件大气腐蚀进行了仿真研究,仿真发现垫圈、2A12与30CrMnSiNi2A偶接接缝处及2A12与TC18偶接铆钉处电流密度分布更为集中,易发生腐蚀,与暴晒试验结果相符,说明该仿真模型对于多电极偶接的腐蚀分布具有很好的预测作用。     

内吹式襟翼环量控制翼型升力响应特性

传统尖尾缘翼型通过控制迎角,综合利用襟翼、缝翼来改变升力,升力对迎角变化的时间响应历程可以用Wagner函数来描述,而内吹式襟翼（IBF）主要通过控制分离来拓展最大升力,并在一定范围内通过调节射流强度改变驻点位置和环量来对升力进行有效控制,其升力随吹气动量变化的时间响应尺度是否与传统尖尾缘翼型相同还不是很清楚。本文主要研究内吹式襟翼升力响应过程,并将其与传统尖后缘翼型升力响应特性进行对比。首先通过某襟翼偏角为30°的双圆弧环量控制翼型对数值方法进行验证,再对某最大厚度为18%弦长的亚声速翼型内吹式襟翼定常吹气控制下的流场进行非定常数值模拟,并分析了其中的瞬态特征。结果表明内吹式襟翼环量控制翼型对激励响应的时间依赖特征与Wagner函数有很好的相互关系,并可以用该函数来描述。     

针式铟场发射电推力器的研制及实验特性 *

为了研制针式铟场发射电推力器实验样机,采用微细电化学腐蚀法制备了场发射电推力器的核心部件发射针,并对发射针进行了表面粗化以促进推进剂的输运,研究了真空环境中发射针的浸润方法,实现了推力器的点火,并测得了不同针尖尺寸的发射针的I-V特性曲线和推力值。实验结果表明,发射针针尖半径决定了点火起始电压的大小,针尖半径越小,点火起始电压越低;推力器样机有连续的推力输出,最大推力可达30μN。     

基于矩化GO法的某型电动飞机扰流板系统可靠性

扰流板冗余系统提升了新能源电动飞机的可靠性。采用GO法建立电动飞机扰流板系统可靠性模型,设计SIMULINK仿真界面,针对操作符11给出一种新的运算模块。针对扰流板冗余系统设计一种矩化操作符,给出运算流程,有效减少传统GO法在冗余系统建模工作量,利用故障树分析法对结果进行验证。结果表明,本文提出的方法切实可行,为新能源电动飞机扰流板系统可靠性分析提供了新的思路。     

前缘钝化和化学非平衡效应对斜爆震发动机进气道特性的影响

为了探究前缘钝化、化学非平衡效应对斜爆震发动机进气道的工作性能及出口温度边界层分布的影响,采用热完全气体、化学非平衡气体两种模型对不同顶板、唇口前缘钝化半径下斜爆震发动机进气道进行数值模拟,结果表明：相比基准进气道,钝化后进气道上壁面温度边界层厚度较大,下壁面温度边界层厚度较小;在化学非平衡气体模型下,顶板前缘钝化半径(R1)>2mm时进气道的顶板附近和分离区内离解反应较为明显,唇口前缘钝化半径(R2)>2mm时进气道的唇罩、唇口板附近离解反应较为明显;当钝化半径≥4mm时,两种气体模型下进气道出口总压恢复系数和静温的相对变化量绝对值大于0.5%,有必要考虑化学非平衡气体效应对进气道出口性能的影响。     

侧板边界层对正激波/边界层干扰特性影响的数值模拟研究

为了研究外压式进气道处于临界工况时结尾正激波同时与压缩面和侧板上边界层的相互干扰,专门设计了三组具有不同前伸长度侧板的简化构型作为研究对象,利用数值模拟手段评估了侧板边界层厚度对正激波/边界层干扰特性的影响。仿真结果表明,当没有侧板边界层参与干扰时流动呈现较好的准二维特性,但当侧板边界层参与干扰后将形成较强的角区干扰结构,该角区三维干扰结构与对称面上的准二维干扰结构存在此消彼长的关系。此外,波系的空间形态也将由“准二维λ波”结构变为“双λ波+角区压缩波”结构,波系形态的改变则进一步导致壁面回流区分布以及摩擦力线拓扑结构变得更加复杂。