整体加筋壳体化铣网格结构参数的概率统计评定

本文应用概率统计理论,对不同的铝合金壁板的化铣网格结构参数做出了正态性分布检验,并根据理论计算和统计图表的综合分析,得出化铣网格结构参数符合正态分布的结论。归纳出在常规的化铣工艺条件和一定概率下的网格参数的偏差值。可作为优化设计时制定化铣网格参数偏差的参考依据。针对化铣工艺特点,阐明了化铣网格参数值和网格重量的关系。     

飞机外挂物投放轨迹捕获的数值模拟研究

本文结合飞机外挂物投放时的运动规律,通过利用复杂外形网格生成技术、六自由度动网格技术和对三维N—S方程的求解实现对外挂物投放时的运动轨迹的数值模拟,分析了不同时间步长对外挂物运动轨迹的影响。     

不同矢量偏角下二元喷管内流特性的数值研究

应用多重网格MacCormack显式格式对不同矢量角下二元喷管的内部流动进行了模拟分析。计算时, 在细网格上采用了雷诺平均N-S方程与k-ε湍流模型方程;在粗网格上则采用Euler方程。数值计算给出了较清晰的流谱, 壁面压力分布的计算也与实验结果相吻合。通过计算结果分析了矢量喷管几何参数对流场的影响, 并对矢量喷管的推力特性进行了数值研究。      

结构自动离散的一种新方法

本文采用等参元素中的几何插值方法,以少量的输入数据,自动生成全部元素的节点坐标;采用“棋盘”网格,使元素和节点的编号简单易行;同时,采用“波前法”解方程组,使编号的程序简化。 本方法具有输入数据少,改变网格容易,边界节点坐标足够准确等优点,已用于轴对称体二维网格及透平叶片的三维网格,结果满意。     

二元矢量喷管内流特性的数值模拟

应用多重网格ＭａｃＣｏｒｍａｃｋ显示格式对不同矢量角下二元喷管的流动进行了模拟分析。计算时,在细网络上采用了雷诺平均Ｎ－Ｓ方程和Ｋ－ε湍流模型方程;在粗网格上则采用Ｅｕｌｅｒ方程。数值计算给出了较清晰的流谱。壁面压力分布的计算值也与实验结果吻合。通过计算结果分析了喷管几何参数对流场的影响。所得结果显示出多重网格算法高效可靠,可用于开展参数研究     

航天器高超声速气动力的数值计算

本文把航天器的外形用数以百计的平面四边形元素来近似。利用高超声速绕流时面积元素之间的干扰可忽略的假设,分别计算出各面积元素的气动力,积分后得出航天器的气动力,算例的结果是令人满意的。本文也描述了航天器表面网格的生成技术和计算机绘制航天器在不同视角下外形图的原理。本文方法可供航天器初步气动设计时用。     

高超声速数值激波失稳的网格依赖性

高超声速强激波的稳定捕捉仍极具挑战性。目前工程计算中普遍应用的数值格式在模拟多维强激波时通常会遭遇明显的激波不稳定现象,且数值格式的激波稳定性对计算网格表现出严重的依赖性。基于矩阵稳定性分析法,对比了具有不同耗散性质的数值格式稳定捕捉激波的能力,分析了空间二阶精度格式的激波稳定性及限制器对激波稳定性的影响。在此基础上,重点探究了计算网格对激波稳定性的影响规律,研究了网格长宽比和畸变角度对激波稳定性的影响。结果显示,在激波附近通过增大网格长宽比或改变网格畸变角度可有效改善激波捕捉的稳定性;相比于增大网格长宽比,改变网格畸变角度提升激波捕捉稳定性的效果更加明显。在此基础上,结合数值耗散分析对高马赫数下数值激波失稳现象的网格依赖性机制进行了探讨。     

基于网格参数化变形的单级涡轮多学科可靠性设计优化

基于自由变形技术发展了一种可以协调实现流场、结构分析网格参数化变形的方法,并实现单级涡轮的多学科可靠性设计优化。针对流场分析网格展向积叠的特征,将三维网格变形转化为不同叶高二维截面的网格变形;根据二维流场分析网格的拓扑结构,设计二维控制面,基于叶型设计参数建立叶型控制点、其他控制点伴随叶型控制点移动,实现了网格的参数化变形,避免了网格畸变、提高了网格变形的质量;不同叶高的二维控制面组成控制体,实现三维流场分析网格的参数化变形。采用相同的控制体进行结构分析网格的变形,实现了耦合界面网格协调变形。基于Kriging代理模型进行单级涡轮的多学科可靠性优化,在满足可靠性约束的情况下效率提高4.97%、寿命提高40.86%,证明了方法的有效性。      

交互式棱柱网格生成中翘曲现象形成机制及消除算法

在基于径向基函数插值的交互式棱柱网格生成中,当物面网格单元空间法矢变化范围较大或变化剧烈时,棱柱网格生成会发生翘曲问题,降低棱柱网格质量。针对该问题,首先采用DLR-F6机翼外形算例对翘曲现象进行了介绍和初步分析,并结合平板外形棱柱网格生成,通过不同波长推进位移扰动模拟法矢变化对推进位移的影响,研究棱柱网格翘曲现象的形成机制。不同位移扰动下平板棱柱网格生成结果证实边界棱线网格点推进位移不光滑会通过空间插值放大,引起棱柱网格发生翘曲,且网格翘曲的表现形式与推进位移的分布有关。从而间接表明棱柱网格生成中翘曲现象是由法矢变化带来推进位移分布不光滑引起的。然后从翘曲现象的形成机制出发,基于拉普拉斯方法发展了推进位移光顺技术,并通过DLR-F6机翼算例进行了验证。测试结果表明,光顺技术能够完全消除棱柱网格的翘曲,获得光顺的高质量棱柱网格。最后通过DLR-F6翼身组合体外形棱柱网格生成,展示了添加光顺技术的棱柱网格生成算法的应用潜力。     

基于工程应用的有限元网格划分研究

首先介绍和分析了划分有限元网格的基本原则,然后基于MSC.Patran有限元软件,结合3个不同的工程实际问题进行网格划分分析,总结了划分有限元网格的一般过程,最后对划分有限元网格进行了小结。     

非结构化网格生成中背景网格对源的响应

背景网格对源的响应方法对非结构网格的生成具有重要意义,如果背景网格参数能恰当地反映出源的变化,则所生成的非结构网格就能满足网格设计的要求。将背景网格的三角形单元分为7个子域,根据源在单元中的位置对网格进行局部加密或变形,使每一个源点都与一个背景网格结点相对应,再通过求解热传导方程,使背景网格参数分布光滑,基本反映网格设计的要求,从而提高了非结构化网格的生成质量。     

混合动网格方法研究

基于Delaunay图映射的动网格方法是一种快速、鲁棒的新型方法,其具有无需迭代、效率高、适用于任意拓扑结构的网格等特点。在分析影响Delaunay图映射方法变形能力主要因素的基础上,对Delaunay图映射方法进行了改进,通过在Delaunay图中增加辅助点,并结合弹簧原理动网格方法,消除了Delaunay图在大变形时出现交叉的问题。通过三维翼身组合体大变形算例验证了发展的混合动网格方法的可靠性和处理气弹问题时的计算效率。     

无蒙皮复合材料网格结构设计与分析

无蒙皮复合材料网格结构是复合材料结构中承载效率最高的结构形式,但其设计成型较其他复合材料网格结构更为困难。本文以工程应用为目的,考虑工艺成型问题,定性与定量分析相结合。按网格形式选择、典型结构工程计算、有限元优化计算、确定工程方案及试验验证五个层次分级优化。最终达到结构形式、设计计算、工艺成型等各方面综合最优的效果。生产出的无蒙皮复合材料网格结构满足工程应用要求,同时找到了工艺成型薄弱环节,为制造最优结构创造了条件。     

电子回旋共振推力器C/C复合材料栅极的力学性能

为了研究电子回旋共振推力器在搭载火箭升空过程中,过载和振动环境对栅极带来的影响,采用有限元分析软件建立了10 cm C/C复合材料栅极的有限元分析模型,计算分析了不同C/C复合材料栅极的力学性能、频率响应和振动模态。结果表明:20g过载下,栅极的最大变形量为6.12×10-4mm;在5~800 Hz激振下,栅极频率响应的最大位移为2.3×10-6mm;过载对栅极的影响要大于振动对它的影响,但是这两个值都在安全门限内;综合分析栅极频响和模态分析的计算结果,栅极仅在其第一阶模态处于外界激振发生共振。     

结构分析程序中几个问题的处理方法

JY851程序是一个复合材料、金属混合结构应力分析及强度校核程序。从1986年投入使用以来,已推广到几个厂、所使用,计算了几个型号的机翼、垂尾、方向舵等复合材料结构。现就编写元素刚度矩阵中所遇到的几个问题的处理方法作一简单的介绍。 1.复合材料板弹性模量阵的计算 复合材料板弹性模量阵A、B、D的计算中,关键是求偏轴模量阵[ō],层压板中第     

有限元网格模型装配技术研究

针对复杂装配体有限元网格模型手工装配难度大的问题,提出了CAE环境下基于工程语义的带装配约束关系的有限元网格装配方法。对有限元网格模型进行自动分组,通过建立CAD装配本体和CAE装配本体之间的映射关系,将CAD的装配约束类型转换为几何位置关系,将这些关系转换为网格模型的装配定位变换矩阵,在对PATRAN二次开发的基础上,通过实例验证了方法的有效性。     

直接统计模拟位置元方法中的分子表面反射确定论判据

在直接模拟Ｍｏｎｔｅ－Ｃａｒｌｏ（ＤＳＭＣ）方法的位置元方法中准确判断分子是否与某一表面元相撞,直接关系着气动力热的计算精度,是方法需要解决的关键问题。本文发展的分子表面反射判据,使这个问题得到了有效的解决。此外,还讨论了分子空间位置坐标的记录方法、分子碰撞对的极限距离、网格自适应等实际模拟技术。与圆球自由分子流准确解和过渡领域贴体网格ＤＳＭＣ模拟结果的比较,证明了上述判据和技术的有效性。     

基于非结构重叠网格的二维N—S方程求解与应用研究

本文发展了非结构重叠网格技术.通过构造任意单元类型网格重叠的网格间边界确定算法、宿主单元搜索算法以及网格间插值算法,把非结构重叠网格推广应用于求解NS方程,并结合湍流模型模拟二维湍流流场.N-S方程的求解采用有限体积法,通量计算采用Osher格式,湍流模型采用SA模型.为验证本文方法的正确性和有效性,模拟了RAE2822翼型和某三段翼型绕流湍流流场,结果表明采用非结构重叠网格模拟复杂外形湍流流场不但网格生成效率高,而且容易推广应用于具有相对运动的非定常问题,具有较高的工程应用价值.     

飞机主承力结构细节分析

工程上往往会遇到如飞机、轮船、汽车等大型结构,对这些大型结构进行有限元分析有时须同时考虑整体以及局部。对于大型结构采用较密的网格,这样会耗费大量时间、资源等,有时会导致计算不收敛;对于局部结构须划分较细的网格才能得到局部构件的细节应力。本文采用子结构法（超单元法）对机翼整体模型中的主承力构件梁进行网格细化,其中梁筋条与肋以及梁缘条与机翼蒙皮采用钉元（弹簧元）模拟多钉连接。在危险载荷工况下对梁进行强度校核和钉载分析,结果表明子结构法能很好地反映大型结构的整体以及细节应力分布,钉元法能很好地模拟机翼多钉连接。     

一种CFD/CSD耦合计算方法

针对柔性大展弦比机翼气动弹性分析和主动弹性机翼(AAW)设计发展了一种计算流体动力学(CFD)和计算结构动力学(CSD)的耦合计算方法。其主要思想是采用在同一物理时间弱耦合求解CFD/CSD技术。气动力采用非定常N-S方程的双时间有限体积求解技术,结构响应则采用有限元数值求解技术。CFD和CSD耦合计算的边界信息(气动力和网格)由所设计的界面程序传输。网格信息传输采用守恒体积转换(CVT)方法将CSD计算结构响应位移插值到CFD网格点上。变形已有的CFD网格技术用以确定CFD的变形网格。以位移或载荷的迭代误差为判断耦合计算的收敛标准。最后得到了机翼在Ma=0.8395,α=5.06°时CFD/CSD耦合计算的收敛值。针对计算结果分析了机翼受静气动弹性过程中结构响应和气动特性随时间变化的效应。初步研究结果表明:这种弱耦合方法求解非线性气动弹性问题是可行的。