用于阻抗变换层的Li2O-Al2O3-SiO2玻璃陶瓷制备及其微波介电特性

 采用热压和凝浇注无压烧结的方法制备了Li₂O-Al₂O3-SiO₂ 玻璃陶瓷阻抗变换层,研究表明烧结的工艺方法和参数对该玻璃陶瓷的致密度和复介电常数有很大的影响。在较高烧结温度和压力下,该玻璃陶瓷的复介电常数保持不变。降低烧结温度和压力时,能够得到一系列复介电常数在较小范围内变化的Li₂O-Al₂O₃-SiO₂ 玻璃陶瓷,但复介电常数不随频率发生变化。将该玻璃陶瓷作为吸波材料的阻抗变换层后,大幅度提高了材料的吸波性能,而且其制备工艺对吸波性能也有一定的影响,可以通过改变该玻璃陶瓷的烧结条件来满足不同吸波材料设计对阻抗匹配层的多种要求。     

机翼亚音速非定常气动力计算

分析机翼非定常气动力特性既可用实验方法[1～3],也可用数值计算方法[4]。现采用非定常涡格时间历程法对机翼沉浮、俯仰、偏航和滚转振荡时的气动力进行计算。1　基本原理与计算公式设图1中t=0时机翼作定常运动;t=Δt时,机翼攻角突变1°,则机翼上产生一涡环ΔΓ1,并以V0向下游运动;t=2Δt时,机翼上又产生一涡环ΔΓ2;随着时间的推进,不断产生以V0向下游运动的涡环,从而模拟机翼的非定常绕流。当涡环以V0运动时(图2),其速度势为φP0=V0Γ4πlimy→0y∫b/2-b/2dz∫x1+V0tx11rdx上式对x、y、z分别求导得到非定常涡环的诱导速度Vx=V0Γb4πx…     

液体燃料火箭发动机性能可靠性的统计评定

 一批液体燃料火箭发动机的性能指标X是。N（μ,σ）,σ²=σ_A²+σ_E²,σ_A²是发动机性能的台间方差,σ_E²是台内方差。设计要求是“μ∈[μ_L,μ_U]、σ≤B”。本文提供了根据发动机的试验数据,能否以置信度γ,判定批产品性能满足设计要求的方法。     

航空发动机高压转子“可容模态”设计及实验验证

为降低航空发动机高压转子的振动，提高高压转子对复杂工况的可容度，本文建立带双阻尼器的高压转子动力学模型，开展模态分析和响应分析，研究高压转子关键结构参数对“可容模态”的影响规律。构造了可容度评价函数η，建立了“可容模态”设计方法。设计了高压转子模拟实验器，通过模态校核实验验证了计算模型和计算方法的准确性，通过“可容模态”实验验证了设计方法的可行性。研究发现，利用当量临界转速ωˉ估计转子系统前两阶临界转速的方法可行，高压转子实验器的临界转速满足ωcr1≤ωˉ，ωˉ≤ωcr2≤2ωˉ；高压转子实验器在两阶“可容模态”下连续“共振”时长达到7125s，实验器转子关键部位的振动位移峰值不超过79.42μm，振动速度有效值不超过2.33mm/s，最大波动不超过10.59%。结果表明，所建立的发动机高压转子“可容模态”设计方法是行之有效的。     

预冷组合循环发动机吸气式模态建模与性能分析

为了对预冷组合循环发动机开展性能分析，以协同吸气式火箭发动机（SABRE 4）为研究对象，采用部件法建立了发动机稳态模型，计算获得了SABRE 4发动机在吸气式模态下沿飞行弹道的性能参数变化规律。然后对发动机的高度和速度特性进行研究，得到了发动机的飞行包线。计算结果表明，在吸气式飞行弹道内，核心机推力和比冲的变化分别为488~680kN和34786~46954m/s。SABRE 4发动机具备推力大和比冲高的性能优势。在预冷器工作过程中，随着飞行马赫数增大，预冷器换热量不断增大，进入预燃室的氢流量减小，预燃室总温降低，HX3的吸热量减小。与其他压气机和涡轮相比，空气压气机和氦涡轮的工作参数变化较大。SABRE 4发动机通过对来流空气进行预冷，可实现在大空域和宽速域内工作。由于空气压气机的喘振和堵塞边界限制，发动机的高度和速度特性分别存在飞行高度和飞行马赫数的限制。     

梯度热障涂层的设计

 采用电子束物理气相沉积方法 (EB PVD)制备了梯度热障涂层,其结构设计为NiCoCrAlY粘结层 /Al₂O₃ YSZ过渡层 /YSZ陶瓷层。YSZ陶瓷层的结构为柱状晶结构,Al₂O₃ YSZ梯度过渡层为梯度微孔结构。采用有限元方法对梯度热障涂层进行热应力分析,优化了Al₂O₃ YSZ梯度过渡层的组成。计算结果表明,梯度涂层的内应力显著降低,而且界面及其附近应力和应变变化较平缓。过渡层厚度的增加有利于降低涂层内应力和缓和涂层界面处的应力集中。     

δ函数与C_nH_m-空气系统通用燃气热力性质表

 本文提出了一个新的热力性质——δ函数。 δ=θ₂-θ₁式中θ₁与θ₂为任意选定的两种基准燃料的θ函数。 对于C_nH_m-空气系统,δ函数只与温度有关,而与燃料氢碳比以及燃气当量比无关。借助于δ函数,编制了适用于任意燃料氢碳比（包括碳与氢）的通用燃气热力性质表。利用此表可以方便而准确地计算C_nH_m-空气系统在任意当量比下的燃气热力性质。文中给出了该表的样表。     

离心力场对内燃波转子火焰传播特性的影响

为了揭示离心力对内燃波转子火焰传播特性的影响，建立了内燃波转子燃烧过程的简化模型，数值模拟了不同离心力场强度下波转子通道内的燃烧过程。研究结果表明，离心力只在一定范围内对火焰传播起到加速作用，且这种影响与预混气当量油气比Φ有关，当Φ=1.0时，作用范围为［0～250g］，在Φ=0.6，0.8时，作用范围为［0～700g］；离心力场对内燃波转子火焰传播特性的影响机制是导致火焰锋面变形，进而引起火焰传播速度变化。     

一种平面埋入式进气道气动特性的试验

 针对一种平面埋入式进气道开展了高速吹风试验研究,获得了沿程静压分布、出口总压恢复图谱、基本气动性能和压力脉动特性。结果表明：沿程静压分布曲线显示,气流绕前唇口流动是先膨胀加速后减速扩压,进入内通道后上壁面静压均高于下壁面;巡航状态Ma₀=0.7,α=2°,β=0°时,进气道总压恢复系数σ=0.951,综合畸变指数W=3.55%,具有较高的性能;当Ma₀=0.6~0.8,α₀=-4°~6°,β=0°~4°范围内,σ在0.912~0.964之间,综合畸变指数在2.68%~7.43%之间,表明该平面埋入式进气道能够在较宽广的飞行包线内以较高的性能安全工作;脉动压力分析表明,出口总压脉动频谱均呈现出白噪声特征,无明显窄带信号出现,这对进气道/发动机匹配工作是有利的。     

Cr对Nb-16Si-22Ti-2Al-2Hf合金显微组织与高低温力学性能的影响

 以Nb-16Si-22Ti-2Al-2Hf成分为基础,分别添加了2%和17%(原子分数,后同)的Cr元素替代Nb(分别称2Cr和17Cr合金),研究了Cr含量对Nb-16Si-22Ti-2Al-2Hf合金相组成、显微组织形貌、室温断裂韧性和高温强度的影响,分析了高低温失效机制。结果表明,铸态和1 375 ℃×100 h热处理后2Cr合金由Nbss和Nb₅Si₃两相组成;当Cr含量为17%时,出现了具有C15结构的Laves Cr₂Nb相,合金由Nbss、Nb₅Si₃和Cr₂Nb三相组成,热处理后在Nb₅Si₃中还析出了球状Cr₂Nb相。随着Cr含量由2%提高到17%,热处理合金的室温断裂韧性K_Q由14.32 MPa·m^1/2下降到10.30 MPa·m^1/2。合金强度与Cr含量的关系受温度影响,随Cr含量提高,室温和1 150 ℃时合金的硬度或强度增高,而1 250 ℃和1 350 ℃时合金强度降低。如1 150 ℃时2Cr和17Cr合金的屈服强度σ_0.2分别为349 MPa和387 MPa;1 350 ℃时分别为306 MPa和74 MPa。     

跨音速非定常流动的Euler方程隐式解

在随时间变化的贴本坐标系中给出求解非定常Euler方程的连续通量分裂法。在此基础上建立了可用于跨音速非定常流动的Euler方程隐式求解法。采用特征向量变换,可在保证原方程组离散化精度的条件下使计算大为简化。针对振动翼绕流特点建立了固连于物体的动坐标与固定坐标间的关系。数值计算在动坐标中进行,既简化网格生成又保证在物面上满足边晃条件。对NACA64A-10冀型绕1/4弦点做简谐俯仰振动的非定常气动力进行了计算,给出了与实验结果基本相符的计算结果。此外,还给出翼型做沉浮及同时进行沉浮与俯仰二自由度振动的非定常气动力的计算结果。     

跨声速机翼操纵面非定常气动力Euler方程计算

采用代数法与椭圆型方程优化相结合的结构生成方法为带操纵面F5机翼生成面对接块结构网格;弹性变形技术生成动态网格;有限体积法、双时间推进法求解非定常Euler方程.给出一种适于面对接块结构网格计算的守恒量插值方法.数值计算了F5机翼操纵面作简谐运动时的非定常气动力,计算值与试验值吻合较好.     

二维非定常欧拉方程的有限元解法

本文通过求解时间相关的欧拉方程，对振荡翼型的非定常二维无粘可压缩绕流作了数值模拟。求解方程时采用了有限体积积分的方法，并用ＶａｎＬｅｅｒ逆风通量分裂格式来计算控制体通量。隐式格式的采用，大大提高了计算效率；同时，通过插值使得空间离散达到了二阶精度。另外，本文采用基于非结构三角形网格的动态网格来处理刚性物体的移动，从而顺利地解决了复杂几何外形物体的非定常绕流的问题。文章最后给出了ＮＡＣＡ００１２翼型绕四分之一弦线作简谐俯仰振动的计算结果。一个周期内的瞬时压强分布以及升力和力矩系数都与实验结果吻合良好，     

一种适合迭代求解的反馈力浸入边界法

对Goldstein提出的反馈力浸入边界法进行了新的思考，改进了其对力源项的计算，拓展了该浸入边界法的使用范围。传统的反馈力浸入边界法在进行力源项的计算时，含有对速度误差的时间积分项，只能用于含时间项的Navier-Stokes （N-S）方程的求解，且在显式时间推进时有严格的时间步长限制。本文改进的方法则直接通过迭代过程中的速度误差求和来计算力源项，避免了时间相关的参数，使其不仅能适合非定常隐式时间推进，还能与不含时间项的定常N-S方程求解方法结合。为了验证改进方法的可靠性，对二维静止圆柱绕流、静止流体中的振荡圆柱、运动椭圆翼以及三维静止圆球的流场进行了计算，计算结果均与文献结果符合较好，表明本文改进的方法是有效的。得出的结论为：可以直接基于迭代次数进行反馈力源项的计算，改进的反馈力浸入边界法不仅可与非定常N-S方程结合，进行隐式求解，还可以与定常N-S方程结合用于定常流动的模拟，可将改进的方法运用到更多的流动问题当中。     

高超声速滚转阻尼导数数值模拟

采用非定常Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程描述物体简谐振动流场,并在Ｅｔｋｉｎ理论下给出绕定轴转动时滚转阻尼导数的计算公式,定常流场的计算采用空间二阶精度的交替方向隐式分解的ＮＮＤ格式,非定常流场的计算采用时、空二阶精度的Ｒｕｎｇｅ－Ｋｕｔｔａ多步格式,采用代数和方法生成物体静、动网格。最后给出高超声速下钝体外形滚转阻尼导数的计算结果,以及滚转力矩系数随瞬时振幅的变化曲线。     

基于CFD/CSD耦合方法的旋翼气动弹性载荷计算分析

为提高直升机前飞状态下旋翼非定常气动弹性载荷的预估精度,在旋翼气动弹性综合分析方法中引入旋翼CFD模块,建立了一套基于CFD/CSD松耦合分析的计算方法和程序。为高效解决流固耦合方法中由于桨叶挥舞、扭转等弹性变形带来的旋翼贴体网格变形问题,采用基于代数变换方法的网格变形技术,桨叶运动变形量和旋翼气动力信息通过流固交接面传递。旋翼流场分析方法中,主控方程采用耦合S-A湍流模型的Navier-Stokes方程,围绕旋翼流场的网格采用结构嵌套网格方法生成,无黏通量计算采用Roe格式,时间推进采用双时间法。旋翼结构分析中,考虑旋翼配平,基于Hamilton变分原理和20自由度Timoshenko梁模型求解弹性旋翼非线性运动方程。分别对CSD和CFD方法进行验证,在此基础上,计算了SA349/2旋翼桨叶在前飞状态下的非定常气动力、挥舞弯矩、摆振弯矩和扭转力矩,并与飞行测试数据进行了对比。计算表明:CFD/CSD耦合方法可以显著提高旋翼非定常气动弹性载荷的分析精度,精确捕捉桨叶表面压强峰值、激波位置等,表明本文发展的旋翼CFD/CSD耦合方法可以有效地运用到旋翼气动弹性载荷的预测分析中。     

超声速平面钝体流动的快速解法

本文介绍了一种平面流动的快速欧拉方程解法。该方法将原参数非定常欧拉方程组重新组合成以广义称曼变量表示的欧拉方程组，再用二点二步迎风格式离散解。针对钝体流动，本文先建立了动网格下的方程，构造了动网格下的算法。提出了一种简单的激波处理方法。计算结果表明，该方法速度快，稳定性好，对初场不敏感。     

欧拉方程多重网格方法及三维跨声速大迎角流动数值计算

本文给出一种粗、细网格时间步长可调整匹配至最优组合的多重网格方法以获得在最细网格上最大的加速收敛效益，二重网格中节省机时可达３５％。文中的数值方法可用于飞行器跨声速大迎角无粘有旋的涡干扰流场的模拟研究。     

基于动态非结构重叠网格法的直升机前飞非定常流场数值模拟研究

 发展了非结构重叠网格技术，提出适合于直升机旋翼运动的动态非结构重叠网格算法，该算法使得在采用非结构重叠网格模拟直升机流场时实现程序自动处理网格，不需任何人工干预。流动控制方程采用非定常可压缩Euler方程，空间离散采用Jameson有限体积法，时间推进采用双时间步长法，为加速收敛，在双时间步长内循环迭代中采用预处理技术。为验证本发展的预处理解算器以及动态非结构重叠网格的算法的可信性和有效性，开展一系列的数值实验，模拟了不同速度范围的NACA0012翼型绕流,Caradonna旋翼悬停和前飞流场以及Robin直升机机身-旋翼组合的前飞非定常流场，结果表明采用本文算法得到的结果是可信的，网格处理是自动高效的，可以应用于直升机全机流场模拟，为直升机设计提供可信指导数据。     

Trajectory and attitude deviations for internal store separation due to unsteady and quasi-steady test method

The deviations of trajectory and attitude angle for internal store separation are evaluated by two wind tunnel test methods. One is the Freedrop Test(FDT), which is known as unsteady and time-dependent method of scaled model. The other is the Captive Trajectory System(CTS) test,which is usually regarded as a quasi-steady and time-averaged test technology. The result shows that there is a streamwise adverse pressure gradient on the cavity resulting in a nose-up pitching moment coefficient(>0) ...