处于非均匀温度场中含裂纹变厚度旋转盘的应力强度因子广义J积分解法

本文首先从直法线假设出发导出了处于非均匀温度场中变厚度各向同性钣元素本构方程并证明该方程与等厚度正交各向异性钣元素的本构方程相同。然后本文从上述本构方程出发,引用自由能与熵函数以及固体力学所有基本方程导出了与积分路径无关的广义J积分J_8,该广义J积分可用于处在非均匀温度场中的含裂纹变厚度旋转盘。最后本文应用这个广义J积分与有限元素法计算了涡轮盘应力强度.因子K_1并得到了令人满意的结果。     

圆轴承三维热流体动力润滑的研究

联立求解了广义雷诺方程、三维能量方程、三维固体热传导方程及载荷平衡方程 ,并考虑粘度随温度及压力、密度随粘度及压力的变化 ,在油膜与轴瓦界面使用热流连续性边界条件 ,得到了圆轴承油膜及轴瓦的三维温度场。结果表明最高温度并不出现在油膜与轴瓦交界面 ,不同圆周位置沿膜厚方向的最高温度出现位置不同。证实了考虑三维温度得到的轴承的润滑性能与考虑二维温度得到的轴承润滑性能有很大的差异。采用Newton-Raphson法求解广义雷诺方程和载荷平衡方程 ,采用快速扫描法求解由有限差分离散的能量方程和热传导方程 ,在计算中发现所使用的算法收敛速度快 ,大大节约了计算时间。     

压电材料修正后的变分原理及智能板的精确解

建立了压电材料修正后的H-R混合变分原理，并从该变分原理直接推导了压电材料的状态向量方程，同时给出了与广义位移和广义应力相关的边界条件表达式。为了验证本文的理论，给出了简支压电材料矩形板静力和自由振动问题的无穷级数精确解形式。并对两个具有代表性的实例进行了数值分析，部分结果与已有的文献进行了比较，并给出了分广义位移和广义应力的分布图，它们直观地展示了精确解的优点。本文的变分公式将用采构造新颖的适应工程领域复杂压电材料问题的半解析法。     

液体火箭发动机稳态特性仿真软件的改进

对作者先前编写的液体火箭发动机稳态特性仿真软件进行了改进, 使仿真软件通过调用用户编写的动态链接库, 添加新的方程和求解算法, 提高了仿真软件的通用性和可扩充性.同时, 仿真软件允许用户直接调整模型中参数的数量和种类, 以及需要求解的方程数量, 使仿真软件能够在不更改模型输入文件的条件下, 进行故障分析、系统参数调整等工作.      

广义矩阵法及其在流体力学中的应用

本文提出了由作者发展的、建立在非笛卡儿张量分析基础上的广义矩阵法,并根据选择的任意坐标系和基,在非正交的一般曲线坐标系中,导出三维空间的流体力学的一系列方程。作为示例,在任意的、不均匀螺旋距的螺旋线坐标系中,用广义矩阵法导出了包括粘性耗散函数表达式在内的、完整的可压缩粘性流动方程。还可导出在两种流体的分界面上,带有粘性及表面张力的应力边界条件。用该法导出的方程,文献已证实是正确的。     

Fourier伪谱方法在不可压缩平板边界层研究中的应用

以三维不可压缩平板边界层为研究对象,扰动形式N-S方程为控制方程,从空间模式的角度,直接数值模拟了三维不稳定T-S波传播的过程.时间离散采用三阶精度混合显隐分裂格式,空间离散则结合Fourier伪谱方法及高精度紧致有限差分逼近,法向采用非等间距网格坐标变换,出口边界条件采用嵌边函数法,程序采用MPI(Message passing interface)并行方法编写.实例验证,该方法计算结果与流动稳定性分析的结果一致.      

基于气动力分段表达的颤振求解方法与模态跟踪技术的应用

在频域颤振方程的求解中,PQI法采用分段气动力的表达形式将颤振方程进行分段直接求解,避免了迭代。由于广义气动力分段插值函数PQI的构建依赖折合频率的选取,将构造的分段插值函数嵌入颤振方程中求解存在模态窜支问题,对此发展了一种基于特征值摄动理论的预测跟踪法来确定各阶模态对应的特征值随参数变化的前后对应关系。文中介绍了PQI颤振分析的基本方法,并采用P-K法进行对比分析来验证算法的计算精度,最后将发展方法应用到某飞机全机轴对称状态颤振分析中来验证算法在实际工程颤振分析中的有效性,数值结果表明该方法计算精度高,模态跟踪效果良好。     

稀薄气体动力学的非线性耦合本构方程理论及验证

作为一种新兴的气体动力学本构方程理论体系,非线性耦合本构方程(NCCR)理论的创新之处在于黏性应力和热传导中抛弃了广义牛顿定律和傅里叶热传导定律,而是考虑熵条件从Boltzmann方程直接推导出了黏性应力和热传导非线性耦合输运方程即NCCR模型。NCCR模型在连续区域与广义牛顿定律和傅里叶热传导定律一致,但是在稀薄区域其非线性关系逐渐增强,即NCCR模型大大扩展了应力-应变和热传导-温度梯度的本构关系,为稀薄气体流动模拟提供了新的途径。为解决NCCR模型强非线性难题,发展了混合模态间断伽辽金求解NCCR和流动守恒方程的数值算法,成功避免了NCCR边界条件高阶量赋值的难题。并对典型亚声速、超声速NACA0012翼型绕流、高超声速圆柱绕流、极高马赫数圆柱绕流、微尺度激波-涡干涉、连续稀薄渐变算例、方腔流动进行了数值计算和验证。结果表明,在稀薄区域,NCCR模型准确捕捉到了流场信息,吻合于蒙特卡罗直接模拟(DSMC)或实验结果,包括:压力分布、速度分布、温度分布、壁面热流等。对圆柱绕流的进一步研究发现NCCR在低努森数下与Navier-Stokes方程结果相同,随着努森数升高两者差距逐渐扩大且在高努森数下NCCR吻合于DSMC和实验结果,从侧面证明了基于NCCR理论用同一套方程解决连续稀薄耦合流动的可能性。     

高效欧拉方程解法求解轴对称引射喷管三股流问题

推导了三维流动的广义欧拉方程黎曼增量梯度形式的欧拉方程及其轴对称形式表达，针对引射喷管三股流问题的特点在运动网格下作了H变换。在此基础上将高效欧拉方程解法应用于求解轴对称引射喷管三股流问题，为解决次流通道狭窄带来的问题，采用了弱耦合法处理滑流边界。     

压电材料壳的广义变分原理和状态向量方程

将三维弹性材料的H-R变分原理引入到具有机-电耦合效应的三维压电弹性材料圆柱壳问题中,建立了对应于正交各向异性压电材料圆柱壳的Hamilton型广义变分原理。通过变分运算得到了柱坐标系下的状态向量方程,该状态向量方程也是Hamilton正则方程。最后指出纯正交各向异性弹性材料壳Hamilton型广义变分原理是变分公式的特例。变分原理的建立将有利于压电材料圆柱壳静力学和动力学问题的有限元方法或半解析法的推导。     

有限元法求解QPNS方程模拟高速旋成体气动力和热

从ＦＮＳ方程中略去流向粘性导数项导出熵变量形式ＱＰＮＳ方程（拟抛物化Ｎ－Ｓ方程）其离散解能自动满足熵增不等式。利用ＧＬＳ（伽辽金／最小二乘法）有限元法构造出ＱＰＮＳ方程的弱解式,所形成的整体方程组用ＧＭＲＥＳ（广义最小余数法）求解。采用空间推进求解ＱＰＮＳ方程,能明显降低计算机内存和运算时间。算例给出典型旋成体的物面压强、剪应力、热流量分布以及流场中密度、温度分布。     

旋翼气动载荷的简化求法

目前,在工程应用上迫切要求发展一种简化的旋翼载荷估算法。 本文从旋翼的广义经典涡流理论出发,通过环量对诱导速度的贡献的简化处理,推导出旋翼诱导速度的二阶分布。然后,根据叶素理论,建立起封闭形式的环量方程,计及桨叶挥舞条件,得出了估算旋翼载荷的简化公式;特别是在挥舞系数的表达式中,考虑了诱导速度的变化,但直接与桨叶参数和飞行参数有联系。 最后,算例表明,本文这些公式适用于直升机气动分析和初步设计。     

用广义乘子法求解航天器最优平面再入轨迹

从纵向平面无量纲航天器质点运动方程出发，引入了能量参数对运动方程进行推导，使得运动方程和优化问题易于处理。通过将优化变量转化为分段连续的线性函数和Runge—Kutta数值计算方法，将轨迹优化问题转化为非线性规划问题，应用广义乘子法对其进行数值分析，给出了不同航程要求下．总加热量最小，满足终端要求以及过载、动压、热流约束的航天器最优轨迹，并分析了其特点。     

喷管型面对固体火箭发动机性能的影响

用时间相关法数值计算跨音速喷管流动时,其控制方程分别采用守恒型和非守恒型方程.通过对JPL喷管进行计算发现:在相同的精度下,非守恒型方程的收敛速度比守恒型方程的收敛速度快一倍多.采用修正阻尼法(DASHCO法)进行数值计算,其收敛速度均提高六倍以上.最后,采用直接优化法研究了喷管收敛段型面及喉部型面对发动机性能的影响.     

用有限元法解稳态蠕变问题的新途径

本文给出直接求解稳态蠕变应力解的有限元方法 ,针对不同材料特性 ,提出了初应变率法与广义切线刚度法。并写出计算公式与步骤 ,研制了分析程序。有算例说明本方法的优越性。     

大展弦比弹性飞机着陆特性影响因素分析

安装点广义位移、结构阻尼系数和机体弹性振动频率对着陆冲击能量耗散效率有显著影响。以大展 弦比全球鹰无人机为例,基于拉格朗日方法建立刚弹耦合的弹性飞机动力学方程,对各模态质量、广义安装位 移、频率和结构阻尼系数等参数下的着陆特性进行仿真,分析各模态质量、广义安装位移、频率和结构阻尼系数 参数对缓冲器载荷和能量耗散时间的影响。结果表明:模态质量越小,安装点广义位移越大,缓冲器载荷越小; 安装点广义位移越大,结构阻尼系数越大,弹性振动频率越小,着陆冲击耗散速率越快。     

中航总科研院所领导视察金航工程

元月15日,中国航空工业总公司的科研院所领导约50人,在总公司张彦仲副总经理带领下视察了设在中国航空信息中心的金航工程网络运行中心,总公司毛德华总师、石川副总师和科技局何怿晋局长,崔德刚副局长和华俊副局长一起参加了视察.金航网络运行人员向各位领导演示了北京、上海、沈阳、成都、西安和闫良六大航空基地主干网联通情况,演示得到了各位领导的好评.     

GMRES算法在悬停旋翼数值模拟中的应用

发展了一种基于广义极小残差(GMRES)算法的悬停旋翼数值模拟方法,并对GMRES算法中矩阵与向量乘积的两种计算方法进行了分析和讨论。应用该方法在旋转坐标系中采用非结构混合网格和格点格式有限体积法对以绝对速度为变量的欧拉方程进行了直接求解,其中对流项的离散应用了基于Roe的Riemann近似解的迎风格式。对Caradonna-Tung旋翼跨声速悬停流场进行了数值模拟,计算结果与相关实验数据吻合较好,并与LU-SGS方法进行了对比,表明GMRES算法可以有效地加速流场的收敛,提高计算效率。     

参数不确定空间机械臂系统的增广自适应控制

讨论了载体位置与姿态均不受控制的漂浮基两杆空间机械臂系统的控制问题,为此对系统的运动学、动力学作了分析。结果表明 :结合系统动量守恒及动量矩守恒关系得到的系统广义Jacobi关系、以及系统的动力学方程将为系统惯性参数的非线性函数。证明了借助于增广变量法可以将增广广义Jacobi矩阵以及系统动力学方程表示为一组适当选择的惯性参数的线性函数。在此基础上,给出了系统参数未知时,空间机械臂末端抓手跟踪惯性空间期望轨迹的增广自适应控制方案。通过仿真运算,证实了方法的有效性。     

先进战斗机控制分配方法研究综述

综述了先进战斗机控制分配方法的研究和应用状况。文中首先介绍了4种主要的控制分配方法———广义逆法、串接链法、直接分配法和数学规划法,给出各种方法在飞行控制系统分配应用中的限制,然后简单介绍控制分配方法在现实飞行控制系统中的应用情况,最后对控制分配中的几个关键问题进行了分析总结。