一种新的判断二维直角坐标网格角点在任意外形物面内外的算法

 根据计算机图形学中的视线遮挡原理,提出了二维平面上判断直角坐标网格角点在任意外形物面内外的一种新算法,并计算了一个 27× 15的直角坐标点阵元素与 3个等腰三角形物面的位置关系,作为该方法的一个算例。计算结果表明,该方法精度高,计算量小,通用性好,并且易于推广至三维复杂情况。     

航天器在轨全过程表面辐射热计算数值仿真研究

对在轨航天器表面辐射热计算进行了全过程数值仿真研究。航天器结构较复杂,针对不同结构进行区域分解,对几何模型进行相应的规则化,同时采用结构化网格和非结构化网格建立通用的计算网格自生成技术。仿真过程重点考虑了任意曲面的网格自动划分和任意形状交界面的数据传递,兼顾几何结构、物理过程、计算精度和计算速度。将有限元法和能束均匀分布法相结合计算角系数和辐射传递系数。将积分法和能束均匀分布法相结合计算外热流。由于在轨航天器表面多用多层隔热组件包裹,针对这部分结构采用节点网络法和控制容积法计算其表面温度,而未被包裹的结构采用有限元法计算其表面温度。对具有辐射换热关系的非连通区域温度场的有限元计算进行了分析和公式推导。最后,用Microsoft VC++6.0编程设计开发了近地轨道航天器表面辐射热计算仿真软件。     

计算网格对气动力气动热数值模拟影响的研究

以二维圆柱钝头超声速流场为研究对象,应用计算流体力学方法研究了计算网格在气动力及气动热数值模拟中的影响.计算程序中对流项、扩散项分别采用AUSM格式和中心差分格式离散,并用多步Runge-Kutta显式时间推进法求解空间离散后的控制方程.计算结果表明,壁面法向第一层网格间距、网格长宽比以及扩张度等因素都会对壁面摩阻、热流的计算产生影响,尤其是壁面法向第一层网格间距影响最为突出.边界层的计算结果与实验数据进行了仔细对比,包括壁面压力分布,摩阻系数Cf分布,斯坦顿数St分布等.高马赫数滞止气流导致壁面附近产生巨大速度梯度以及温度梯度,这些梯度的准确计算需要合适的壁面法向第一层网格间距以及合理的网格长宽比、扩张度.     

基于STL文件和有限元法的在轨航天器关键部件热计算

针对STL(stereo lithography)文件在传递复杂几何实体模型信息方面具有精度高的特点,提出了一种基于STL文件计算复杂结构角系数和外热流的方法.根据单元和节点的拓扑关系识别六面体网格边界单元并进行辐射换热计算.详细阐述了基于STL文件和有限元法进行复杂结构的角系数和外热流的计算方法.研究了导热-辐射耦合计算的有限元方法.最后利用有限元法计算航天器关键部件在轨的三维瞬态温度场.      

基于梯度下降的水滴收集率计算方法

拉格朗日方法在计算三维复杂外形的水滴撞击特性时，一般需要通过插值获得物面网格上的水滴收集率，而这一过程容易造成计算结果失真。建立了一种基于梯度下降的水滴收集率计算方法，通过确定撞击物面网格角点的水滴轨迹实现物面网格中心处的水滴收集率计算，避免插值运算。该方法的计算思路如下：根据撞击点到目标点的距离，采用梯度下降法自适应地调控水滴的初始释放位置，使撞击点落在壁面网格角点附近的给定阈值范围内；采用邻点检索算法，快速确定迎风面内所有撞击网格角点的水滴轨迹；基于流管定义，计算出网格中心处的水滴收集率。采用该方法对典型三维外形的水滴收集率进行了计算，并与相关文献结果进行对比。结果显示计算结果与文献数据吻合良好，有效克服了现有技术中因插值导致的结果失真问题，能够为飞机结冰及防除冰技术研究提供参考。     

非结构网格限制器研究

研究了6种限制器形式及其参数计算方法,其中新的自适应VenkatakrishnanM限制器能够根据当地流场特征自动调节参数从而实现不同区域不同特性.数值试验表明:非结构网格限制器对邻居集合选择方法较敏感,其中共面法比共点法精度高约15%;高马赫数时压缩型限制器比耗散型限制器计算精度高约10%~40%;参变量算法中,方向梯度法得到的流场比增量法平滑;自适应VenkatakrishnanM限制器能够模拟复杂超声速流动.      

反舰导弹三维流场数值模拟

三维流场的数值模拟是国外军事技术发达国家在飞行器设计过程中所采用的重要方法.本文建立了导弹的真实三维几何外形,然后结合实验数据,采用有限体积法,对国产某典型反舰导弹在高亚音速、不同攻角时的三维流场进行了数值模拟.计算出了导弹在不同工况下的阻力系数、升力系数和力矩系数.把计算所得结果与实验值进行了比较,用理论公式对计算结果进行了对比分析,并简单地分析了偏差的原因和解决的可行措施.文中所用的残差光顺技术和网格自适应技术能有效地提高计算的收敛速度和精度.     

基于人工黏性的二维激波小波多尺度数值计算

基于自适应小波配点法和人工黏性，构造了二维激波问题的小波数值计算格式。利用最细尺度的小波系数构造两组激波定位函数，分别用以控制两个正交方向人工黏性的大小和分布。对三类二维激波问题进行计算和验证，结果表明：自适应小波配点法利用小波阈值滤波删除变化平缓区域大量网格点，而保留变化急剧区域的网格点，比传统方法计算效率高，且分辨率越高，计算效率更高；利用最细层小波系数构建的幂函数形式的激波定位函数能准确判断激波位置和控制X、Y方向人工黏性的大小和分布，从而捕捉不同方向和不同强度的激波。     

混合网格方法在栅格翼数值模拟中的应用研究

栅格翼是采用一种蜂房结构设计的升力面和控制面,和传统舵面相比具有铰链力矩小、能在较大攻角下保持升力等特点,并得到了越来越广泛的应用。由于栅格尾翼的特殊结构,网格生成具有很大的难度。本文采用了结构、非结构混合网格的办法,在弹体和栅格翼的物面附近区域,生成适合粘性计算的大长宽比的结构网格,弹体网格和栅格翼网格之间采用非结构网格进行填充,满足非结构网格和结构网格交接面的完全对接。本文基于结构/非结构网格体系采用有限体积方法求解NS方程,对不同舵偏角下的栅格翼构形进行了数值模拟,并通过实验结果对数值方法进行了验证。     

不同喷口修形的二元收敛喷管RCS数值模拟

用等效棱边电磁流(EEC)方法和物理光学迭代(IPO)法自主开发了计算程序,用来分析不同尾缘修形喷管的雷达散射特性;在计算中考虑喷管出口边缘的绕射场对雷达散射截面(RCS)的影响;提出了一种高效的判断面元间遮挡关系的算法和喷管出口轮廓线的辨别方法;在此基础上,研究了四种不同喷口修形的二元收敛喷管的RCS特性.结果表明:该开发程序具有较好的计算精度;适当的喷口修形可降低重点姿态角下喷管的RCS.      

弹性模量对铝胶焊搭接接头应力分布的影响

运用利用弹塑性有限元方法,研究了胶粘剂弹性模量对铝合金胶焊搭接接头上工作应力分布的影响.研究结果表明:随着胶粘剂弹性模量的增加,胶焊搭接区焊点中心线上的应力分布趋于均匀,胶层-焊点边缘处的峰值应力降低,但靠近搭接区端部胶层中的各应力峰值增大;距胶层中心线0.15mm的面上与焊点中心相对应点的峰值应力随弹性模量提高而下降.为提高胶焊接头承载能力,建议采用弹性模量适中的胶粘剂.     

聚硫推进剂燃烧条件下裂纹扩展过程的研究

研究了影响聚硫推进剂试件内裂纹扩展的各种因素。对含裂纹的推进剂试件进行了大量的燃烧试验,其燃烧过程用Ｘ射线实时成像系统进行了记录;并用粘弹性有限元方法计算了试件在燃烧过程中的应力应变状态,利用Ｊ积分法对裂纹扩展的可能性进行了预估。理论分析和试验结果均表明：燃烧室增压速率、裂纹尺寸、边界条件等因素对裂纹的扩展均有较强的影响。该研究结果对制定固体火箭发动机装药失效判据有一定的参考价值。     

真空电子束钎焊温度场数值模拟

针对不锈钢毛细管板类结构件真空电子束钎焊的特点，利用ANSYS软件建立了一个三维有限元分析模型对电子束钎焊温度场进行模拟。通过将热电偶测得的工件表面温度与有限元分析结果比较，证明所建模型的合理性。利用该模型对不锈钢圆板电子束钎焊温度场进行了预测，证明通过电子束扫描加热方式可以得到钎焊所需的局部均匀温度场。     

Optimization of Preform Shapes by RSM and FEM to Improve Deformation Homogeneity in Aerospace Forgings

This article, in order to improve the deformation homogeneity in aerospace forgings, proposes an approach that combines the finite element method (FEM) and the response surface method (RSM) to optimize the preform shapes. New expressions that take into account the influences of equivalent effective strain distribution are developed to evaluate the homogeneity of deformation distribution in aerospace forgings. In order to reduce the number of design variables, the domain-division method is put forward to determine the optimal design variables. On the basis of FEM results, the RSM is used to establish an approximate model to depict the relationship between the responses (deformation homogeneity and die underfilling) and the design variables represented by geometric parameters of the preform shape. With a typical aeroengine disk as an example, the proposed method is verified by achieving an optimal combination of design variables. By comparing the preform shape obtained with the proposed method to that with the existing one, it is evidenced that the former could achieve more homogeneous deformation in forging.     

三维机翼亚临界定常流的局部线性化有限元法

本文使用了两种类型等参元,通过变分法,计算了三维机翼亚临界定常位势流中的压力分布,计算结果与试验结果符合较好。文中还使用了有限元法中的局部线性化理论,插值出元素中心点密度去处理迭代过程中元素系数阵的计算,其结果与通常的有限元计算相吻合,但计算时间却明显减少了。文中还对有限元法局部线性化理论和局部线性化有限元法做了阐述。     

燃烧条件下影响推进剂脱粘面扩展的因素

研究了固体火箭发动机推进剂与包覆层界面脱粘在燃烧条件下扩展的影响因素。对含有预制脱粘面的试件进行了大量的燃烧试验，用X射线实时成像系统记录了燃烧过程；根据组成试件的不同材料特性，分别用粘弹性和线弹性有限元方法计算了试件在燃烧过程中的应力应变状态，并利用修正后J积分法对脱粘面扩展的可能性进行了预估。理论分析和试验结果均表明：燃烧室增压速率、脱粘面尺寸、边界条件等因素对脱粘面扩展均有较强的影响。该研究结果对制定固体火箭发动机装药失效判据有一定参考价值。     

有限元法在标准直齿圆柱齿轮齿弯曲疲劳强度计算中的应用

本文利用有限元法计算了标准直齿圆柱齿轮轮齿应力，得出其分布规律以及危险截面应力最大值。与传统计算结果相比较，其结果更为准确，为轮齿弯曲疲劳强度计算提供了一种简洁、精确的计算方法。     

平箔片楔形高度对气体止推箔片轴承特性影响

针对波箔型气体止推箔片轴承,建立了箔片结构二维薄板模型,并通过有限差分法和有限元法耦合求解可压缩气体Reynolds方程和气膜厚度方程,获得了给定轴承载荷条件下轴承气膜压力分布、气膜厚度分布、平箔片变形量和功率损耗等轴承特性.通过对比楔形高度分别为25,70,100,200,300μm时的轴承特性仿真结果研究了楔形高度对轴承性能的影响.结果表明:降低楔形高度使轴承气膜压力分布更均匀,并降低了平箔片的局部集中载荷.但楔形高度存在一个最佳值,使得达到相同轴承载荷所需的最小气膜厚度最大,并且轴承具有最小功率损耗,提高了轴承的工作效率.该结果为气体止推箔片轴承的结构设计提供的理论参考.      

大涡模拟进口合成湍流方法

进口瞬时湍流脉动条件直接影响大涡模拟的精度。对比研究了更适合工程应用的合成涡和合成谱方法,并通过考虑进口速度时间相关和随机数的概率密度分布,对合成谱方法进行改进。对零压力梯度平板湍流模拟结果的对比表明:不同湍流合成方法得到的合成湍流均存在发展为真实湍流的过渡过程;考虑进口湍流脉动的时间相关有助于合成湍流更快地发展为真实湍流,需要的流向距离更短;合成谱方法采用正态分布系数优于均匀分布,改进的合成谱方法优于常规合成谱方法。