基于混合笛卡尔网格的翼型绕流数值模拟

在对重叠网格技术研究的基础上,一套以笛卡尔网格为主体、贴体的结构网格为特殊固壁边界的混合笛卡尔网格(Hybrid Cartesian Grid,HCG)生成方法得到了发展,同时基于流场特征的笛卡尔网格自适应技术和有限体积方法的雷诺平均Navier-Stokes(RANS)的数值求解方法也得到了开发。通过对经典跨音速RAE2822翼型、双NACA0012翼型和高升力两段NLR-7301翼型进行绕流数值模拟,获得了较好的结果,验证了方法的可行性与精确性,为研究二维复杂外形流动问题提供了一种解决方法。     

预混气爆轰马赫反射实验研究

为研究斜爆轰马赫反射的实验现象以及不稳定性对马赫反射的影响,在矩形爆轰管道中,对3种预混气的马赫反射现象进行了实验研究,取得烟膜记录。烟膜结果显示,三波点高度在楔角附近与无反应三激波理论类似,随后逐渐偏向反应三激波理论,三波点实际高度由反应实际放热量决定,三种气体偏转始末距离与结束偏转距离之比平均值分别为0.42,0.49,0.31,与三者活化能大小关系相符合。胞格结构随初始压力增大而渐趋规律,胞格尺寸也逐渐减小,由CJ区域到过驱区域胞格结构存在明显的转变。横波是影响斜爆轰稳定的重要因素,压力增大可抑制分子的无规则运动,使横波出现的位置更加集中,胞格更为规律。预混气爆轰不稳定度大的气体横波出现的时空位置更不易确定,胞格结构更不规律,胞格尺寸较小。     

浅谈子超格

目的研究超格的子超格.方法利用格的等价定义及超代数运算,引入了超格和子超格的概念.结果给出了超格的理想的概念,并且研究了它的一些相关性质.结论超格和子超格概念的引入,对于超格理想的研究是有实际意义的,并且完善了超格的商结构.     

旋翼诱速的横向分布

 直升机前飞时,旋翼平面处诱速的横向分布,究竟前行桨叶半边大还是后行半边大,迄今无定论。本文根据经典涡流理论,全面分析环量Γ_o及Γ_(1s)对诱速横向系数υ_(1s)的贡献,然后利用桨叶挥舞条件,通过Γ_o及Γ_(1s),不同分布的组合情况进行估算,发现存在着一个临界半径位置。在临界半径以外υ_(1s)为负,即前行桨叶的诱速小于后行半边的;在临界半径以内υ_(1s)为正,前行桨叶的诱速大于后行桨叶的。上述结果得到实验证实。     

气相爆震波马赫反射转变机理的实验研究

在矩形爆震管道中,选用具有规则胞格结构的高浓度氩气稀释稳态气体C2H2-2.5O2-8.17Ar及不规则胞格结构的非稳态气体C2H2-5N2O,利用纹影和烟膜实验对CJ爆震状态下的爆震波发生马赫反射的转变机理进行了研究。在楔角为30°,5k Pa初始条件下,C2H2-2.5O2-8.17Ar马赫反射转变时马赫杆临界高度约为1.25cm,过驱胞格结构出现的最小距离略小于胞格长度的1/3;C2H2-5N2O马赫反射转变时马赫杆临界高度约为1.5cm,由于其胞格结构不规则,故难以衡量过驱胞格结构出现的位置。实验结果表明:爆震波马赫反射三波点运动轨迹特征为在初始阶段遵循无反应冲击波理论,随后马赫反射发生转变使得运动轨迹逐渐平行于反应冲击波理论,从而验证了爆震波开始形成前导激波马赫反射,由于CJ区域内三波点的介入引起扰动,最终在马赫反射区形成过驱爆震的过程。     

二维多段翼型网格生成与数值模拟研究

为提高多段翼型的网格生成效率和数值模拟精度,发展了一套自适应混合笛卡尔网格(AHCG)生成方法和基于有限体积方法的雷诺数平均Navier-Stokes (RANS)的数值求解技术.混合笛卡尔网格由围绕物体几何外形的贴体结构网格和填充流场其他区域的笛卡尔网格构成,两套网格之间的信息传递由“贡献单元”提供,且“贡献单元”由基于ADT(Alternating digital tree)技术的搜寻方法获得.为更准确地捕捉流场信息,采用了基于流场特征的网格自适应技术.数值模拟结果显示,AHCG方法能够准确且高效地模拟高升力多段翼型绕流问题.     

基于CFD动网格的内装式空射分离研究

研究了使用飞机发射运载火箭时的机箭分离过程。提出了一种可以解决内装式空射分离问题的整体解决方案,主要基于CAD三维实体建模方法、多体动力学方法和结合六自由度运动的计算流体力学(CFD)动网格方法。以CFD动网格计算为核心,进行了空射分离过程关键问题的研究,得到了火箭在载机气动干扰环境下的运动特性,为轨道卫星和弹道导弹的快速、廉价发射打下了研究基础。     

爆轰胞格尺寸的统计分析

爆轰胞格尺寸作为可燃系统的本征值,可用来确定爆轰传播时的临界尺寸,比如直管传播时的临界管径、管道突扩时的临界尺寸等。然而,该尺寸的测量具有较强的主观性。为了减少人工测量带来的不确定性,文章采用两种统计方法:概率密度函数法和自相关函数法,对数值模拟得到的不同规则程度的爆轰胞格进行了统计分析。爆轰胞格的不规则程度采用单步反应的活化能来体现。当活化能较低时(Ea=15),胞格较规则,两种统计方法得到的胞格尺寸较一致。随着活化能的增加(Ea=20),胞格开始不规则,两种统计方法也开始出现偏差。当活化能进一步增加(Ea=25、27),胞格极不规则,开始出现横波的合并,两种统计方法得到的结果出现较大的偏差。这是由于概率密度法处理爆轰波阵面三波点规则时,没有考虑迹线上能量的分布,而自相关函数法根据三波点上的涡量值对胞格尺寸进行统计分析,因此,该方法得到的爆轰胞格尺寸比概率函数法更能真实体现可燃系统的爆轰特性。     

基于三阶龙格库塔法的铣削稳定性半解析法预测

针对使用再生颤振理论建立的铣削动力学模型,提出了一种基于三阶龙格库塔法用于预测铣削稳定性的半解析方法.首先,以状态空间方程的形式表示动力学微分方程;其次,利用三阶龙格库塔法推导出传递矩阵;最后利用Floquet理论判断特定切削状态下的稳定性,进而获得铣削稳定性叶瓣图.通过与半离散法的仿真结果进行对比发现,基于三阶龙格库塔法的铣削稳定性求解方法具有更高的预测精度和计算效率.     

基于CFD动网格的内装式空射分离研究

研究了使用飞机发射运载火箭时的机箭分离过程.提出了一种可以解决内装式空射分离问题的整体解决方案,主要基于CAD三维实体建模方法、多体动力学方法和结合六自由度运动的计算流体力学(CFD)动网格方法.以CFD动网格计算为核心,进行了空射分离过程关键问题的研究,得到了火箭在载机气动干扰环境下的运动特性,为轨道卫星和弹道导弹的快速、廉价发射打下了研究基础.     

基于二维紧凑型四阶格子模型的激波管流动模拟

为了提高格子的稳定性，使用Hermite展开方法，构建了新的二维四阶紧凑型格子模型，即D2Q37A。比较了D2Q37A和与Philippi给出的紧凑型格子模型(D2Q37B)的稳定性。在相同的碰撞频率下，与D2Q37B相比，D2Q37A可以模拟初始密度比更高的一维激波管流动。这表明D2Q37A与现有格子模型相比，具有更好的稳定性。详细给出了适用于高阶格子模型的边界条件实现方式。此边界条件实现方式保留了体现LBM（lattice Boltmann method）粒子特性的迁移 碰撞机制。用以上给出的格子模型和边界条件处理方式模拟激波管流动，得到的模拟结果和解析解吻合得很好。这表明所给出的边界处理方式是可行的。此边界格式同样可以用于其他类型的流动和边界。      

点阵单元类型和排布角度对传热特性的影响

基于瞬态液晶测试技术和数值仿真,在相同孔隙率条件下,研究了桁架型点阵结构的单元类型和排布角度对矩形通道壁面对流传热特性的影响.结果表明:点阵单元类型对壁面传热性能有决定性的影响,其中交叉排布的体心立方(body centered cubic-0,BCC-0)和Kagome-0点阵归因于出色的传热强化效果,雷诺数在5 7...     

模拟MKDV方程的格子Boltzmann方法

目前,格子Boltzmann方法已被广泛应用于模拟各种非线性物理方程。文中用D1Q2模型给出MKDV方程的带修正项的BGK型格子Boltzmann方法。通过对演化方程的泰勒展开并应用多尺度技术恢复了宏观方程。数值模拟表明文中所建立的模型是有效的。     

一类边传递图

令G=（α,β｜α^n=β^2=1,α^β=α^r）,r＆^2≡1（modn），是图Г的一个自同构群。目的是研究关于G一边传递图的性质．运用置换群和代数图论的相关理论，获得了这类图的完全分类，它们是一些互不相交的圈和完全二部图的并。     

关于雨流计数法的讨论及修正

在疲劳载荷的统计处理中,雨流计数法一直被认为是比较合理的循环计数方法,在国内外得到了广泛应用。本文结合对应力-应变响应的研究指出了它的不足,提出了修正雨流计数法。此外,对于现有的两种利用有效系数模拟应力-应变响应的方法也进行了讨论,并提出了修正有效系数法。     

压力对有序LixTiS2电子性质的影响

采用基于密度泛函理论的局域密度近似（LDA）,计算了压力对LixTiS2（x=0.0．25,0．5,0．75,1）的晶格结构和电子结构的影响。研究发现,对于同一“含量的LixTiS2,随着压力的增大,LixTiS2的晶格参数有规律的减小,插层化合物中Li的含量越高,材料越难被压缩;在压强等于1GPa左右时各体系的总能量达到最小值。各体系的费米能级随着压力的增大而降低;各原子的电子能态密度产生了漂移,压力还导致了S原子和．Ti原子p-d杂化减弱。“原子和S原子的共价相互作用得到了增强     

碳纤维增强点阵夹芯结构的屈曲强度

考虑到点阵芯层的结构形式及杆单元以拉压为主的变形模式,基于连续介质等效理论,提出了一种全新的位移场假设,并以此编写了计算点阵夹芯复合材料屈曲强度的有限元程序.通过与不同试样的压缩试验的比较,对碳纤维增强点阵夹芯结构的屈曲行为进行了分析,并且验证了程序的有效性.同时,采用程序,讨论了点阵参数(包括杆单元长度、半径及倾斜角)对点阵夹芯复合材料结构屈曲强度的影响规律.结果表明:芯层体积分数相同的情况下,金字塔型和四面体型点阵夹芯结构的屈曲强度相差不大,后者略高.所得结论对点阵夹芯复合材料结构设计具有一定的指导意义.      

直升机六力素识别

本文介绍载荷识别技术。提出了“应变传递矩阵法”载荷识别的基础理论、不同性质载荷的识别方法和公式。叙述了直升机飞行载荷-六力素识别的三个主要阶段:传递矩阵的确定（地面标定）、结构响应实测（空中飞行试验）、载荷计算（六力素计算）。文中还就随机载荷识别进行了讨论。     

WAM-FM遥测系统

 根据正交分割理论,现在流行的频率分割制（FDM）和时间分割制（TDM）都属于正交分割的范畴。满足正交条件的信号不仅仅只有频分和时分两种形式。1923年美国数学家沃尔什（WALSH）提出的沃尔什函数就是一种正交完备的函数系。本文在理论上提出了一种以沃尔什函数为基础的序率分割制（SDM）,并以实验证明了新方案是可行的。文中详细地介绍了这种SDM系统的设计思想、方框图、工作原理以及具体实施中的技术问题。实验室研制成的样机是沃尔什调幅——调频制的,称为WAM-FM遥测系统。为了阐明新系统的优越性及其潜力,把它与现行的频分遥测系统作了简单的比较。说明新系统不仅具有一定的学术价值,而且具有一定的实用价值。     

增材制造钛合金微桁架夹芯板低速冲击响应

增材制造技术能够制造复杂点阵结构。相比于传统的加工工艺,可以一次成型,克服了低速冲击下传统工艺芯层与面层在连接点处易发生脱粘的问题。利用低速落锤试验装置对增材制造面心立方（FCC）夹芯板和体心立方（BCC）夹芯板进行了低速冲击试验,获得了两种微桁架点阵夹芯板的破坏模式和冲击响应曲线。低速冲击下,微桁架夹芯板上面层在冲击部位产生局部凹坑,并出现裂纹,其余部位没有大变形。试验结果表明在相同能量冲击下,BCC夹芯板的凹坑深度要小于FCC夹芯板,BCC夹芯板的抗冲击性能要优于FCC夹芯板;建立有限元模型,较好地表征了低速冲击过程中微桁架结构的损伤。发现在低速冲击过程中,对于两种微桁架点阵夹芯板,冲击能量主要由上面层和芯层吸收;冲击能量改变,夹芯板各部分吸能百分比变化较小。BCC夹芯板和FCC夹芯板结构稳定,整体性好;低速冲击下,FCC夹芯板最先发生破坏的部位是上面层与芯层连接处;而BCC夹芯板最先发生破坏的部位是中间竖直桁架。