用时间相关法计算喷管跨声速流场

本文用时间相关法完成了定常、无粘、跨声速喷管流动的计算。内部点参数用MacCormack差分格式计算,壁边界点参数用简化了的特征边界条件计算,轴上点参数用反射原理计算。计算表明:计算是收敛的,计算结果与实验数据符合良好。     

低裂解度正癸烷物性快速计算方法

为研究正癸烷轻微裂解时物性计算方法,基于广义对应态法则建立起正癸烷裂解过程中的密度、黏度、导热系数和定压比热容的三种计算方法,包括直接计算、物性库插值计算、物性子库加权计算,在验证直接计算方法的准确度基础上,对比三种方法的计算精度、计算内存和计算时间来对其综合评估,并探究裂解度对高温物性的影响。计算温度变化范围为300~1020K,压力变化范围为1~15MPa,裂解度变化范围为0~0.25。结果表明:不同物性计算方法均能定性预测热物性的特殊变化趋势;以物性直接计算结果为基准,物性子库加权计算的误差与裂解度大小成正比;在不考虑计算内存的情况下,物性库插值计算方法能无损加速计算,低裂解度下快速工程计算可选择物性子库加权方法;裂解度对高温物性的影响显著且不可忽略。     

飞机主馈线单相交流电阻解析计算方法

以大型飞机变频交流(AC)电源系统为研究背景,针对主馈线敷设方式缺乏有效的交流电阻解析计算方法的现状,研究导体间具有通用性的邻近因子解析计算方法,能够理论计算宽变频范围内两根导体在不同材料、不同半径、不同中心距离下的交流电阻值,该计算方法具有求解过程简单、物理意义明确、适用范围广、精度高等优点。通过对比航空电缆交流电阻的解析计算、有限元仿真计算和实验计算结果,验证了该方法的有效性和正确性。该计算方法解决了传统使用有限元仿真工具计算交流电阻时过程繁琐、耗时长的问题,简化了计算步骤,具有重要的工程应用价值,并为多导体主馈线敷设方式交流电阻解析计算方法的推导奠定基础。     

具有颤振约束的结构优化计算

为了提高新产品设计水平,我们编制了具有颤振约束的结构优化计算程序。主要内容包括:柔度系数插值、自振计算、非定常气动力计算、颤振方程求解和优化准则。在计算方法上,本程序作了如下改进:自振计算采用L_R型方法;对非定常气动力计算中与型态有关的部分采用预先存储的办法;颤振计算采用零阻尼插值等。从而大大地提高了本程序的计算速度。在709计算机上,利用本程序计算一个颤振点,仅需4分钟。我们用本程序对某产品进行了试算,结果表明,本程序是可靠的。     

动态尾涡间隔计算工具的开发与实现

基于建立的尾涡消散模型、运动模型、遭遇模型,以Delphi 7.0为平台,采用面向对象的设计思想,开发了动态尾涡间隔计算工具,包括前机诱导尾涡流场计算、后机遭遇响应过程分析、安全间隔计算等功能模块。根据输入的飞机参数、气象数据、机场数据,实现对尾涡间隔的实时、高效计算。基于开发的计算工具,以B777-200和B737-800为前后机进行了验证计算。     

基于多核并行计算的超燃冲压发动机一维模型实时性研究

在发动机控制系统设计中,为了缩短设计周期、降低研发成本,需要建立面向控制的、较为精确的、实时性高的超燃冲压发动机性能计算模型,以保证模型精度、提高计算速度为研究目标,基于多核高性能计算仿真平台,开展了面向控制的超燃冲压发动机一维模型实时性优化工作。运用简化计算流程、改进C语言程序、开拓缓存区等方法有效提高了一维模型计算速度。创新性地尝试了计算流体力学并行化方法,对隔离段和燃烧室一维模型进行结构分解。计算网格平衡分配至多个中央处理器,并借助核间数据通讯实现多核并行计算。与串行模型计算结果对比,七核并行计算模型性能参数偏差不超过0.1%,全工况仿真时间小于30ms,计算耗时较优化前缩短了75%以上。实时性优化后的多核并行模型计算精度高、速度快、收敛性好,可以作为超燃冲压发动机控制系统设计和半实物仿真验证平台。     

计算气动光学研究进展

计算气动光学是采用数值计算的方法研究高速空气动力流场对光波传输和光学成像影响的一门交叉学科,在红外成像制导设计等工程应用中发挥着重要的作用。总结了计算气动光学研究的三种基本计算方法,包括基于CFD简化方法的光学统计估算、基于RANS流场的计算和基于LES/DNS瞬态流场的计算,对这三种计算方法的基本理论、计算过程进行了介绍,综述了新近的一些计算气动光学研究进展,最后指出了未来的研究重点和方向。     

基于CFD数值仿真技术的飞行器动导数计算

为了研究表征飞行器动态稳定特性的动导数,结合计算流体力学方法,基于滑移网格技术,采用小幅度强迫振动法和差分法两种非定常方法,推导了俯仰组合动导数以及滚转阻尼导数的计算表达式,并结合多参考系模型采用准定常方法计算滚转阻尼导数。建立了基于计算流体力学技术的飞行器动导数计算方法。对国际动导数标模Finner导弹进行了计算验证。两种非定常方法对俯仰组合动导数的计算误差分别为0.65%、6.13%,对滚转阻尼导数的计算误差分别为2.5%、0.06%。求解滚转阻尼的准定常方法对于滚转阻尼导数的计算误差为2.67%。三种方法的计算结果与风洞试验结果吻合的很好。非定常方法计算精度高,准定常方法可以快捷迅速得到结果。本文的数值方法不仅适用于超声速,同样可以用于亚、跨声速时的动导数计算。     

高性能计算在飞机设计中的应用

介绍了高性能计算在国内外的发展现状,分析了高性能计算在飞机设计中的重要应用方向与发展需求,包括:气动力预测、噪音分析、滑流研究、雷达散射截面计算、飞行器多学科设计优化、空投空降和多体分离研究等。根据当前需求,建议发展百万亿次每秒的计算能力,并支持国有自主大型高性能计算应用软件开发工作。     

考虑涡轮传热性能的气动设计耦合计算方法研究

为了在涡轮叶片气动型线设计过程中同时提高气动效率并保证传热性能,提出了一种基于一维管道网络法和三维CFD的耦合计算方法,分别采用管网/三维耦合计算方法和全三维耦合计算方法对MARK-II冷却叶片多个工况进行计算,两种数值计算方法计算结果与实验数据交叉对比,以验证本文计算方法可行性。计算结果表明,两种数值计算得到的叶片型面压力、温度、换热系数和实验值都比较吻合,但管网/三维耦合计算得到的壁面温度相比全三维耦合计算结果整体略微偏低,最大偏差不超过3.89%。基于管网/三维耦合计算方法对某航空发动机涡轮第二级动叶叶片型线优化设计,气动效率提高0.34%,壁面平均温度几乎没有变化。     

斜削梁协调元素族及其在梁固有振动分析上的应用

本文给出了位移函数为三次、五次、七次幂多项式时斜削梁元素刚度与质量矩阵的显表达式以及计算任意奇次幂高阶斜削梁元素族刚度与质量矩阵显表达式的一般方法。然后应用高阶协调元素计算了六种线性斜削梁在悬臂支持、二端简支、二端固支等情况下的头五阶固有频率及振型。结果表明,在所有情况下,与低阶元素相比,高阶元素均显出有较好的收敛性。但若施加以力的边界条件,则在使用少量元素和梁的斜削度较大的情况下,高阶元素的结果反而不如低阶元素的结果好。     

正交各向异性材料的三维弹塑性分析

本文根据Hill屈服准则以张量形式推导出正交各向异性材料的弹塑性刚度,并进而导出在三维和二维情况下对应的弹塑性刚度矩阵的显表达式。所得到的公式在形式上与各向同性材料的弹塑性刚度矩阵十分类似。所以。根据本文的公式,用有限元求解时,就能很方便地将现有的各向同性体弹塑性分析的实用程序改换为正交各向异性体弹塑性分析的实用程序。作为例子,选择了一个非线性有限元通用程序,用本文公式在其三维元素的材料库中建立了新模块,并就三维悬臂梁进行了计算。     

旋翼桨叶动力分析及内力计算中有限元法的改进

提高用有限元法进行桨叶动力分析、特别是动内力计算的精度是本文工作的主要目的。本文建立了变剖面旋转梁协调单元族;提出利用单元动刚度矩阵直接由节点位移计算桨叶内力（简称为计算内力的“动刚度法”）。作为例子,用所导出的不同精度的有限单元和建议的内力计算方法计算了具有不连续结构特性的变剖面桨叶弯曲振动固有频率、振型和模态内力。计算结果表明,工作可以达到预期的目的。     

基于接触刚度分布特性的叶片振动响应分析

根据接触刚度在锯齿形叶冠结合面上非均匀分布的特性,提出一种基于定义和有限元计算相结合的接触刚度计算方法。在此基础上,将微-宏滑动摩擦模型作为叶冠结合面处的摩擦模型,推导微滑动和宏滑动状态下摩擦力表达式。利用谐波平衡法将非线性摩擦力转化为等效刚度和等效阻尼进行振动特性分析。针对叶冠结合面相对位移幅值动态变化的特点,提出了一种迭代求解的振动响应分析方法。与文献提供的带摩擦阻尼结构悬臂梁振动响应实验数据相比,共振峰附近的振幅误差为1.99 mm,相对误差为3.9%。振幅的最大误差为5.53 mm,出现在远离共振峰的位置,验证了响应分析方法的可行性。将该方法应用于带冠叶片上,结果表明当激振力频率为812.3 Hz时,振幅为0.56 mm。     

复合材料臂杆刚度性能研究

为获得大型复合材料臂杆刚度特性用于结构动力学仿真,采用有限元建模和试验测试两种方法,分别得到试验件一阶模态频率与弯曲及扭转刚度。分析结果与试验结果并不吻合,对产品刚度降低原因进行了分析,不良成型工艺和纤维打断对复合材料产品的刚度特性有着较大的影响。     

基于MSC.Nastran优化梁元偏置的分析

在有限元结构设计分析中,对梁的受力形式以及在空间放置方式等的理解直接影响整个分析计算结果。文中将加筋板离散为板元和梁元,用大型通用有限元商业软件MSC.Nastran的SOL200优化了梁元考虑偏心和不偏心的4种工况,其结果满足位移和应力约束。对几种工况中优化后的加强筋的不同尺寸,表明存在巨大差异,且符合真实物理模型的有限元模型,即带梁元偏心的情况质量达到了最轻。将优化结果参数用于原模型,用MSC.M arc进行了模态分析,并比较了前三阶模态。最后,得出了结论,如设计中不考虑梁元偏心,会给梁元尺寸的优化结果带来极大的误差,并对刚度的影响也不可忽略。考虑梁元偏置的加筋板不仅具有最小质量,而且具有广泛的实用价值。     

复合材料闭剖面加筋薄壁梁刚度系数计算分析

复合材料的可设计性为通过弹性剪裁来获得想要的变形模式带来了优势,其结构的耦合特性如拉扭耦合和弯扭耦合可以利用复合材料的各向异性来实现。以薄壁壳结构力学理论为基础,采用多种理论方法研究分析正交各向异性壳体本构关系,给出适合复合材料层压板矩形闭剖面薄壁梁截面的刚度系数解析表达;采用逐阶近似方法,完整设计复杂闭剖面的刚度系数数值算法;针对两种典型的复合材料铺层的矩形闭剖面梁截面布局构型,进行刚度特性的计算分析与讨论,获得铺层角度对复合材料薄壁梁弯曲、扭转及弯扭耦合刚度特性的定量结果。对分析设计闭剖面薄壁结构刚度条件具有应用价值。     

拉杆组合转子的刚度修正及动力学建模

为研究拉杆组合转子的动力学特性并比较其与整体式转子的差异,基于集总参数法,提出了一种拉杆组合转子动力学建模的方法.该方法通过对整体式转子梁模型抗弯刚度的修正,可以考虑接触、预紧力和拉杆等因素的影响.首先分别分析了这3个因素对等效梁静态抗弯刚度的影响,给出了相应的修正系数.研究结果表明:接触界面和拉杆对轴段抗弯刚度的影响很大,而轴向预紧力的影响较小.然后,针对一模化实验拉杆组合转子,计算了各轴段的修正系数并建立了动力学模型.其临界转速计算值和实验结果吻合,误差仅0.8%,表明了方法的正确性.研究可为其他组合式结构的动力学分析和设计提供参考和借鉴.      

转子-定子系统碰摩可靠性灵敏度研究

针对仅考虑质量偏心的Jeffcott转子-定子系统碰摩的动力学方程,进行了旋转机械转子-定子系统的碰摩可靠性灵敏度分析,应用Kronecker代数,矩阵微分理论,四阶矩技术,矩阵摄动理论,概率统计方法,可靠性理论和灵敏度技术系统地发展了转子-定子系统碰摩的可靠性灵敏度分析方法,研究了转轴刚度和阻尼、定子刚度和阻尼、偏心距和定子径向刚度对可靠度灵敏度的影响,并给出了数值分析解.     

充气机翼的褶皱与失效行为研究

在经典工程梁理论的基础上,结合张力薄膜的应力状态分析,提出充气机翼褶皱失稳的判据。计入表面薄膜褶皱引起的刚度退化效应,将机翼等效处理为一个变截面刚度的梁,建立了充气悬臂机翼的等效梁模型,并采用微分求积法进行充气机翼弯曲变形分析。计算结果与充气机翼的静力弯曲试验结果相吻合,验证了充气机翼弯曲变形分析方法的有效性。应用片条理论引入气动力模型,并与所建立的等效梁模型相耦合,建立充气机翼的静气动弹性耦合模型,并用迭代算法进行求解。考虑起皱和失稳两种判据,并计算获取了试验机翼的起皱动压和皱褶失稳动压形式,计算结果与风洞试验结果一致。根据所建立的充气机翼静气动弹性分析方法,可以预测充气机翼表面褶皱区的扩展和弯曲变形,进而绘制充气机翼的静气弹许用包线,为充气机翼的设计提供必要的安全边界参考依据。