纤维增强复合材料细观温度场计算

对高阶理论的原始算法进行了改进,采用界面的平均温度替代了假设温度函数的系数做为求解的未知量,并利用了子胞的热传导方程,建立了热流与平均温度间的关系,进一步求解出温度场。改进后使求解的方程数目减少了大约60%,求解时间大大缩短;并且消除了亚子胞的概念。结果表明,该方法计算的复合材料细观温度场结果与理论解和有限元解具有较好的一致性。      

基于任意四边形子胞的多尺度温度场计算方法(英文)

基于功能梯度材料的高阶理论,提出了基于任意四边形子胞的多尺度温度场计算方法。以二维问题为研究对象,首先采用任意四边形子胞离散求解域,并且采用二次多项式来近似子胞内的温度场。通过采用子胞边界温度代替假设温度函数的系数减少未知量数目,应用Fourier热传导定律,建立单个子胞边界平均热流密度与边界平均温度的关系。将此关系代入边界条件以及子胞间的平均热流平衡方程和平均温度连续条件,消去子胞边界热流密度。然后以子胞边界平均温度作为未知量对方程组进行求解。从而求解出整个物体的温度场分布。本文提出的方法突破了原高阶理论中子胞必须是规则长方形的限制,提高了高阶理论的使用范围。通过与有限元的计算结果进行比较,证明了本文方法计算的有效性。     

改进的三维高精度通用单胞模型

基于复合材料细观力学周期性假设,对高精度通用单胞模型进行了改进,并推广到三维。模型中用界面的平均量代替细观位移函数中解系数,并利用细观单元力学方程的分析与求解,建立了细观平均量与复合材料宏观平均量间的联系.改进的高精度通用单胞模型的求解方程数目减少了大约55%,求解时间缩短,并且消除了亚子胞的概念.算例分析表明该模型的计算结果与细观有限元结果具有较好的一致性.      

基于离散伴随方程的三维雷达散射截面几何敏感度计算

针对有限差分法计算雷达散射截面（RCS）梯度效率低，采用高精度雷达散射截面评估时计算代价高的问题，提出了一种基于麦克斯韦积分方程离散伴随方程的RCS梯度高效计算方法。基于伴随方程的梯度计算可以通过一次雷达散射截面求解、一次伴随方程求解获得RCS关于所有设计变量的梯度。其中麦克斯韦积分方程离散伴随方程的形式与原方程基本一致，可以采用矩量法（MOM）及多层快速多极子算法（MLFMA）求解。伴随方程求解计算量与直接雷达散射截面评估基本一致，存储量在直接雷达散射截面评估的基础上增加不明显。通过双椎体模型、导弹模型对基于矩量法、多层快速多极子算法的伴随梯度进行验证，证明了基于伴随方法的RCS梯度计算可以实现复杂外形中RCS梯度的高效、高精度求解，为基于梯度的高精度气动/隐身一体化优化提供了基础。     

柱面气相散心爆轰波胞格演化规律研究

本文采用二阶精度NND格式,应用改进的二阶段化学反应模型,通过求解二维Euler方程对柱面气相散心爆轰波胞格演化过程进行了数值模拟.计算结果表明在传播过程中,空间尺度的扩张导致了散心爆轰波后气流的自然膨胀,使得多波结构的爆轰阵面呈现出显著的胞格自组织特性.根据计算结果与理论分析,本文归纳了五种胞格演化模式,分别命名为内凹波阵面会聚、波阵面扭结、褶皱波面失稳、胞格自合并和三波点滑移,并定义了各种模式的物理特征,分析了其相关的演化机制和规律.     

航空发动机故障诊断主因子选取相似系数法

针对目前航空发动机气路故障诊断主因子模型中计算量大、故障方程多重共线性的问题,提出了主因子选取的相似系数法。通过计算故障模式和测量参数偏差之间的相似系数来选取故障方程的主因子,避免了对全部可能的主因子组合的计算。讨论了利用方法求解时的方程共线性问题,并对比了直接利用主因子模型求解的计算结果,结果表明,提高了故障诊断结果的可判断性,使故障诊断更加迅速可靠。     

耦合Navier-Stokes方程的可压缩边界层数值求解

为了在Reynolds-averaged Navier-Stokes(RANS)方程计算中耦合eN方法进行转捩判断,在RANS方程求解过程中耦合求解了可压缩层流边界层方程为判断转捩提供了精确、可靠的边界层信息.利用等熵关系由RANS方程求出的物面压力分布确定边界层外边界的速度分布,进一步确定出边界层外边界.边界层方程的求解使用Keller提出的二阶BOX方法.为了验证方法求解边界层方程的正确性,在低速流动状态下将计算结果和XFOIL的计算得到的边界层解进行了比较,二者吻合较好.     

弱界面黏结通用单胞模型理论及应用

将通用单胞模型与弱界面分离模型结合.首先推导了增量型通用单胞模型计算公式, 然后采用界面位移公式对平均位移连续条件进行改造, 建立子胞平均应变增量和宏观应变增量的函数关系.应用子胞增量本构方程得到子胞平均应力增量表达式, 采用均匀化方法获得宏观应力增量表达式.最后采用渐进适应界面模型模拟了长纤维复合材料的横向拉伸过程.算例分析表明, 将界面分离模型和通用单胞模型结合能较好地模拟复合材料在拉伸载荷下的宏观响应.      

一种强耦合Spalart-Allmaras湍流模型的RANS方程的高效数值计算方法

在工程实际中,一方程湍流模型或两方程湍流模型的求解通常和雷诺平均Navier-Stockes (RANS)方程的求解是解耦的,也称之为松耦合求解.在松耦合求解过程中,RANS方程和湍流模型方程通常采用不同的数值方法异步求解.这种求解方式很容易产生因两者计算精度不一致而引起的额外数值耗散.为了消除这种耗散,将RANS方程与Spalart-Allmaras模型方程耦合成一个系统方程——强耦合RANS方程,并发展了一种用于求解该系统方程的高效强耦合算法,其中对流项离散采用了Roe格式,时间项的离散采用了隐式LU-SGS(Lower-Upper Symmetric Gauss-Seidel)格式,为了提高计算效率,采用了三层V循环多重网格方法.通过翼型/机翼和振荡翼型/机翼等算例验证了本文发展的强耦合算法不仅具有较好的收敛性,而且计算精度明显优于松耦合算法,特别对于阻力的预测,强耦合算法更加准确.     

考虑传质传热效应的翼型积冰计算

本文采用计算流体力学的方法对NACA0012翼型前缘涉及传质传热效应的glaze ice积冰进行了预测.在流场计算中,用有限体积法对二维定常不可压粘流的时均N-S方程进行离散求解;采用四阶龙戈-库塔(RK)法求解水滴运动方程;通过求解在翼型表面控制体积内积冰过程中所遵循的质量和能量守恒方程,并且假定冰沿着与当地翼型表面法向一致的方向增长,来预测积冰的形状.文中同时对结有glaze ice冰翼型的气动特性进行了计算与分析.     

弹流润滑螺旋锥齿轮热摩擦行为分析

针对直升机主减速器中的高速重载螺旋锥齿轮,建立了点接触热弹流润滑分析数学模型,包括Reynolds方程、能量方程和膜厚方程等,采用数值方法求解弹流润滑状态下的齿面摩擦因数.模型中考虑了润滑油黏度和密度随压力、温度的变化,并通过轮齿承载接触分析(LTCA)获得齿面真实载荷和卷吸速度、相对滑动速度等运动学变量.基于真实齿面点建立了螺旋锥齿轮单齿模型,考虑滑动摩擦生热和不同表面上的热边界条件,通过有限元稳态热分析和瞬态热分析得到了轮齿本体温度场和接触点瞬时闪温,并与现有文献和算例齿轮台架试验结果进行对比.      

考虑辐射效应的多孔材料传热有限元算法

应用有限元法求解一维半透明多孔介质辐射对流传热问题。通过对辐射传递方程与边界条件中积分项的离散 ,推导出有限元求解列式。对于多孔介质传热控制方程和气体传热控制方程的求解 ,采用时间积分技术对两相温度统一进行求解。而对于边界方程处理 ,则通过传热方程与边界方程的循环求解 ,求出任意时刻温度场的瞬态解。给出数值算例 ,得出了放置在通风口、同时受到伴随辐射作用的半透明多孔介质传热的瞬态解 ,讨论了部分参数对瞬态温度场以及换热效果的影响      

航空燃气轮机挤压油膜阻尼器油膜温度场的理论分析

油膜温度是挤压油膜阻尼器的一项重要工作参数。本文运用有限元法对挤压油膜阻尼器油膜温度场进行分析,即根据油膜边界条件,联立求解适用于阻尼器的雷诺方程、能量方程、滑油粘温特性方程得到油膜温度分布。文中进一步分析航空发动机阻尼器的工作环境,提出启动阶段,滑油温升由于阻尼减振发热所致;稳态阶段主要由于油膜环向滑油的强迫对流换热,并提出用试凑法来求解稳态滑油工作温度。     

直升机动力舱流场及温度场的模拟

应用非结构网格,建立了直升机发动机动力舱的有限元模型,并根据舱内实际冷却气体的流通情况,采用SIMPLE算法,求解了流动和传热的控制方程,模拟了动力舱内流场与温度场的分布,研究了流场、压力场及温度场之间的依变关系。在此基础上针对现有设计存在的问题,提出了动力舱冷却系统的改进方案,并对不同方案的改进效果进行了分析比较。      

基于有限单元热网络法的航空发动机主轴轴承热分析

航空发动机主轴轴承温度场的精准预测对保证发动机润滑系统及整机稳定运转至关重要.针对传统热网络法精度低的问题,提出一种基于有限单元思想的有限单元热网络法.将轴承结构通过网格划分为有限个单元体,每个单元体设置一个温度节点代表该单元体的集总参数,节点间构建热阻关系并形成热网络,通过求解以温度为未知量的大型稀疏矩阵线性方程组来...     

倾转旋翼航空器建模方法研究

结合倾转旋翼航空器的整体气动结构设计,建立了该型航空器的数学模型。在考虑旋翼滑流的情况下,给出各气动面的气动力和力矩方程,并采用叶素法,建立了旋翼的气动力方程。该数学模型可为倾转旋翼航空器的控制系统设计与仿真提供参考。     

燃气轮机叶轮非线性瞬态温度场的边界元分析

首先采用克希荷夫变换将包含变物性及混合和辐射边界条件的热扩散方程转换成线性方程，进一步由加权余量法和离散化将方程转变成边界元方程组。用所开发的二维线性—非线性瞬态温度场边界元程序分析计算了燃气轮机叶轮起动工况下的非稳态温度场，与有限元数值计算结果比较表明，边界元法应用于这类非线性工程问题的分析是可靠、初始精度高，数值稳定性好及计算输入数据简单等优点。     

基于非协调单元的液体火箭发动机推力室热结构分析

用巴兹公式和发动机的热力计算数据,得到了推力室内的燃气壁面对流传热系数和燃气温度的分布,并用有限单元法计算了推力室的温度场.在温度场计算基础上,采用Wilson非协调单元,对推力室的应力场进行了计算,其中位移函数采用二阶形式,将计算精度提高到二阶.结果表明:相比常规单元,Wilson非协调单元在计算结构比较复杂,温度梯度比较大的结构时,精度比较高,结果合理.      

有限元传热分析中燃气辐射处理方法及其在飞机结构热分析中的应用

 研究了有限元传热分析中有吸收-发射性介质( 或称弱吸收性介质, 例如燃气) 参与辐射传热时的处理方法。通过对有、无吸收-发射性介质的封闭包壳的辐射能量交换方程的分析, 提出了用表面边界单元替代其辐射效应的准则; 研究了有吸收-发射性介质参与辐射时表面间辐射角系数计算方法及编制了相应的计算程序。最后, 对某型飞机机尾罩的温度场作了计算并对结果进行了分析。