盘─鼓结构体─壳过渡元的“虚假应力”现象及改进方法

分析了体—壳过渡元中产生“虚假应力”现象的根本原因，给出了两种改进方法。方案一是在采用体元应力应变模型的基础上提出了一种单元内位移、应变及应力的直接修改方法；方案二是将采用壳体应变应力状态的三维体—壳过渡元推广到轴对称，给出了具体公式。通过几个算例将上述两种方法与改进前方法进行了比较，结果表明本文提出的两种方案均大大改善了轴对称问题体壳过渡元的“虚假应力”现象，相比之下，前者公式简单，更易于实现。     

子母弹抛壳过程非定常流动的数值模拟

利用弹簧近似和网格重构相结合的非结构动网格技术,耦合求解Euler方程及弹道方程,模拟飞行在46千米高空、0°攻角、马赫数6.5、壳片分离初速分别为6m/s,12m/s、壳片分离角速度分别为0°/s,50°/s,150°/s的子母弹抛壳动态过程。研究表明,抛壳初期,壳片与母弹周围形成复杂流动,在壳片内表面存在激波的生成、发展、减弱等非定常过程,导致壳片受到气动力随着分离距离增加表现为很强的非线性特征;壳片分离初始角速度对壳片姿态的影响不是线性的,尽管分离初始角速度增加到150°/s,由于转动惯量小、俯仰力矩大,壳片实际飞行姿态达不到设计要求,壳片攻角较小,气动阻力小,导致壳片伴随母弹较长时间飞行;增大分离初速是实现快速分离的有效措施。     

不同厚度航行器高速入水冲击载荷及壳体变形特性

航行器高速入水过程中,壳体受冲击载荷作用将产生大变形,而大变形又反过来影响航行体的受力状态。为了了解航行体结构受入水冲击的动态响应特性,基于结构化任意拉格朗日-欧拉（structured arbitrary Lagrange-Euler, S-ALE）算法研究了航行器撞水阶段壳体所受冲击载荷与变形的关联性,获取了壳体变形模式,并分析了壳体厚度对变形的影响规律。结果表明：壳体内凹大变形将增大航行器所受法向载荷,载荷峰值的脉宽虽为“毫秒级”,但对壳体变形具有较大影响。入水速度较低时,壳体变形模式以弹性变形为主,变形区域在上凹和下凸间不断转换,增大壳厚,结构主要通过减小脉宽来抑制变形;入水速度较高时,壳体变形模式以塑性变形为主,依次出现内凹区、拉伸区和卷曲压缩区,随形变量增大,拉伸区不断扩张,而卷曲压缩区不断后移,其中内凹变形受斜入水影响向y+方向偏斜,而卷曲压缩区在壳厚较小时出现S形卷曲对,增大壳厚,结构主要通过减小峰值来抑制变形。航行器斜入水内凹变形非对称性的内在机理为弹性应变与塑性应变的叠加效应,当冲击能量不变时,增大壳厚提升了内凹区壳体对冲击能量的吸收能力,减小了拉伸和卷曲压缩形变量...     

采用壳体应变—应力状态的体—壳过渡元

由于初始热应变表达式不当及不能完全实现泊松效应,以往建立的体-壳过渡元会导致不真实的大应力。本文采用壳体应变-应力状态作为过渡元中的应变应力状态,给出~种新的15节点体-壳过渡元列式,从而解决了以上两个方面的问题。算例表明了这种过渡元的有效性与实用性。本文方法既适用于其它类型的过渡元,对于其它处理体-壳组合结构的方法也有重要参考价值。     

高超声速风洞多体干扰与分离试验技术

在FL-31风洞中进行了某高超声速飞行器的多体干扰与分离试验技术研究,成功建立了多体干扰与分离试验技术。试验模型是某典型构型的可重复使用航天飞行器,由助推器以及再入体两部分组成。利用风洞上下投放机构实现两模型间的相对运动,采用两台天平对模型的气动力进行测量,同时利用纹影仪记录模型分离过程中的激波干扰情况。结果表明：试验系统设计合理,能准确模拟物体间分离过程,并能精确测量多体干扰的气动力特性,激波干扰清晰可见。     

高超声速多体干扰与分离试验

进行了高超声速飞行器多体系统分离过程中存在的气动力干扰试验研究.试验模型是某构型的可重复使用航天飞行器,由助推器以及再入体两部分组成.研究在FL-31风洞中进行,试验马赫数为Ma=6.97.试验结果表明:分离过程中助推器和再入体之间存在复杂的激波干扰现象,多体系统分离过程中的气动干扰本质上是激波干扰引发的.      

椭球形舞台壳体风荷载风洞实验研究

通过风洞模拟试验研究了椭球形舞台壳体，及其三种开启状态时风荷载特性。结果表明就整体风荷载而言椭球体的阻力系数不是很大，但由于气流的分离和再附，椭球体的局部会出现很大的负压，系数可达－2．0以上，另外，周围建筑对椭球形舞台壳体的风荷载影响很大。     

采用壳体应变应力状态的正交异性体─壳过渡元

采用壳体应变-应力状态的过渡元推广于正交各向异性材料, 以解决壳体法线与一个材料主轴重合的正交异性体-壳组合结构的应力分析问题。算例表明:采用本文建立的过渡元计算正交异性体-壳组合结构是有效的, 精度能满足工程要求, 并能节省计算时间。      

固体火箭发动机壳体及喷管模态分析

发动机模态参数是检验发动机有限元动力模型、进行振动响应分析的基本数据.本文选用锥壳单元,介绍了轴对称旋转壳模态分析方法.计算了发动机壳体和喷管的固有频率和振型,并与实验结果进行了对比分析.     

发动机机匣、支板系统结构动力分析

发展了在环向只能处理连续变量的常规有限条法。分析在环向由若干个离散分布的支板和支架支承的发动机机匣回转壳。各离散支板和支架回转壳刚度的贡献, 是通过支板和壳体连接处位移相等的变换, 将各支板的刚度阵叠加到回转壳的总体刚度阵上实现的。同时证明了位移在环向的Fourier级数展开式中, 只有n=0, n=1项对轴系的支承刚度有贡献。支板及支板附近的壳体的局部应力, 是由局部子结构空间板壳有限元分析得到, 该子结构的边界位移由整体的回转壳分析给出。算例表明, 本文提出的方法能有效地计算轴系的支承动刚度和壳体的局部应力。      

高超声速升力体气动力气动热数值计算

升力体气动布局，由于其良好的高超声速气动特性、有效的内部空间和有益的防热特性，越来越受到气动工作者的重视，本文应用张涵信院士的NND格式思想，用数值模拟的方法对类似于美国X－33的升力体外形进行了气动力、气动热研究。特别是在粘性项的离散和表面热流率的计算中应用了散度和梯度的积分定义，避免了数值计算的奇异性，保证了通量守恒，提高了计算效率。结果表明气动热计算结果与美国同类飞行器结果接近，驻点热流率与经典理论结果一致。     

NND格式在多维理想磁流体方程组中的应用

采用修正的四步Runge-Kutta方法求解三维一般曲线坐标系下的理想磁流体方程组,为克服数值振荡,加特征型NND格式进行后处理.特征型NND格式推广到求解三维磁流体(MHD)问题需要知道雅可比通量的左右特征矩阵,在具体计算时需要克服矩阵的奇性.本文用三维程序采用推广的特征NND格式计算了一维MHD激波管和二维(MHD)喷管流动,计算结果表明,特征NND格式保持了TVD格式高精度的优点,又具有计算简单的特点,在包括强弱间断等复杂波系的定常和非定常MHD流场数值模拟中是成功的.     

矩形/非结构混合网格技术及在二维/三维复杂无粘流场数值模拟中的应用

建立了一套模拟复杂无粘流场的矩形／非结构混合网格技术,其中非结构网格仅限于物面附近,发挥非结构网格的几何灵活性,以少量的网格模拟复杂外形;同时在以外的区域采用矩形结构网格,发挥矩形网格计算简单快速的优势,有效地克服全非结构网格计算方法需要较大内存量和较长ＣＰＵ时间的不足。混合网格系统由修正的四分树（２Ｄ）／八分树（３Ｄ）方法生成。本文将ＮＮＤ有限差分与ＮＮＤ有限体积格式有机地融合于混合网格计算,消除了全矩形网格模拟曲面边界的台阶效应,同时保证了网格间的通量守恒。该算法被成功地运用于ＮＡＣＡ００１２翼型的跨、低超声速绕流,二维弯管内运动激波的非定常流,以及三维捆绑火箭的超声速绕流的数值模拟。计算结果表明本方法在模拟复杂无粘流场方面的灵活性和高效性。     

再人体复杂外形无粘流场数值计算研究

本文在文[1,2,3]的基础上研究了非定常Euler方程的推进迭代方法,并将无波动、无自由参数的耗散差分格式(下称NND格式)发展为隐式、迎风格式,用以计算有差曲控制翼的再入复杂飞行器的育攻角、有侧滑角的超声速无粘绕流流场。其结果是准确的。本文的方法具有以下特点: 1.可以在同一计算程序中实现亚跨超流场的计算:在超声速区域,采用推进技术;在亚跨声速区,则自动选择推进迭代求解技术。 2.NND格式不仅具有公式简单,无自由参数的特点,而且过激波时满足熵条件,不会产生非物理解,也不发生振荡。 3.隐式差分格式的无条件稳定和Gauss-Seidel迭代的快速收敛,计算省时。 本文的方法曾用于烧蚀头部凹陷外形和航天飞机简化外形的计算,能自动选择迭代和推进区域。可以预言,如计算机条件许可,本方法可以用于大攻角的无粘流场计算。     

非结构网格上新型的NND有限元格式

张涵信发展的ＮＮＤ差分格式是由中心差分格式,二阶迎风格式和一阶迎风格式混合组成的杂交型格式。众所周知,和中心差分格式相对应的是Ｇａｌｅｒｋｉｎ有限元格式。通过对中心型有限元格式加修正项的方法本文成功地构造出二阶迎风型有限元格式和一阶迎风型有限元格式。     

超音速来流中侧向喷流干扰流场的数值模拟

通过求解三维可压缩N S方程组,对尖锥/圆柱组合体模型在超音速绕流中的侧向喷流干扰流场进行了模拟,数值方法采用有限体积法,差分格式采用NND格式。分析了干扰流场的结构,研究了干扰效应对气动性能的影响,得到来流Ma∞=3 3和4 5时力放大因子和干扰力矩系数随攻角变化的规律,计算结果与实验一致。     

悬停旋翼下洗流对导弹的气动特性和初始弹道的干扰研究

用广义尾流模型计算了悬停旋翼下洗流场,在此基础上用工程方法计算了有旋翼尾流干扰和无旋翼尾流干扰两种条件下的某空空导弹的气动特性.将气动力计算和初始弹道计算相耦合,计算了导弹的初始弹道,给出了气动参数沿初始弹道的变化曲线,并对计算结果进行了分析.数值计算结果表明旋翼下洗流对导弹的气动特性和初始弹道有显著的影响,在尾迹区内导弹的气动参数变化较大,尾迹区外气动参数变化比较平缓;弹道参数与无下洗流干扰条件下相比也有很大不同.本文揭示了旋翼尾流对导弹的气动干扰规律和导弹的初始运动规律,研究成果对改进机载导弹的发射技术和弹道控制有重要的意义.     

亚跨声速风洞机弹干扰试验技术

介绍了一种利用小尺寸风洞进行亚跨声速挂弹/载机干扰试验的试验技术,该试验技术是在小尺寸风洞内进行挂弹/载机局部干扰试验,可以测量挂弹投放过程中作用于挂弹及其折叠尾翼上的气动力和舵效率。为验证该试验技术,选取了某空射型导弹进行了挂机干扰试验验证,试验前对局部干扰流场和完整干扰流场进行了数值计算,计算结果表明局部干扰流场与完整干扰流场流态基本一致。试验结果表明,小风洞对局部流场的模拟可以取得合理的试验结果,该试验技术具有可行性。     

旋翼尾流对空空导弹的气动干扰研究

直升机发射导弹时,旋翼尾流对导弹产生很强的气动干扰,研究旋翼尾流对导弹的气动干扰对于改进机载导弹发射技术具有重要意义。结合动量理论和涡流理论,推导了桨叶环量方程,用广义尾流模型计算了旋翼下洗流场。并在此基础上用工程方法计算了有旋翼尾流干扰和无旋翼尾流干扰两种条件下的某空空导弹气动特性,揭示了旋翼尾流对导弹的气动干扰规律,为直升机载导弹的发射控制提供了参考依据。     

跨声速轴流压气机特性全三维粘性流动分析

应用Steger-Ｗarming通量分裂技术,将守恒型方程中的流通向量分裂成两部分,应用 NND格式对无粘项进行数值离散,粘性项采用中心差分格式,采用LU-SGS隐式推进迭代建立了跨声速轴流压气机性能及内流场三维N-S方程的高分辨率和高效率的数值方法。数值模拟给出NASA ROTOR37和NASA ROTOR67性能及内流场流动结构。并进一步基于混合平面法建立了多级轴流压气机性能计算模型,详细分析了某一级跨声速轴流压气机性能和内流特性,计算与实验结果比较表明了本模型与方法是可行的和可靠的。