基于超网格的重叠网格守恒插值方法

保证重叠网格边界数据插值的守恒性是计算流体力学面临的一大难题。基于超网格方法与格心格式有限体积法,发展了一种新型混合重叠网格守恒插值方法。详细研究了用网格切割技术在重叠网格边界构造局部超网格,用网格求交算法合理地扩大贡献单元模板,建立了一种适用于任意网格类型的隐式并行重叠网格守恒插值方法,以超网格为媒介可实现重叠网格二阶精度的守恒插值。数值结果表明,本文方法对二阶分布的流场变量具有严格的守恒性,相比三线性插值方法和逆向距离权插值方法,本文方法减小了数值误差,提高了重叠网格边界的插值精度,加快了计算收敛速度,改善了重叠区域网格尺度相差较大时流场的光滑性和连续性。     

复合材料网格结构的研究现状

介绍了复合材料网格结构的网格分类、特点以及国内外复合材料网格结构的发展应用概况.介绍了复合材料网格结构的设计方法,以及各种成型工艺的优缺点.叙述了复合材料网格结构型模的分类和各种型模的优缺点,并指出了修正方法.同时,介绍了复合材料网格结构性能评估和测试方法.     

基于非结构背景网格的动网格方法研究

借鉴弹簧动网格方法和Delaunay背景图方法,提出一种基于非结构背景网格的动网格生成方法。先用少量的物面控制点生成相对较粗的非结构背景网格,再用弹簧法将动边界的变化量作用到该背景网格,然后利用计算网格与背景网格的映射关系插值生成对应的动网格,最后对物面网格点及其附近受影响的密网格点进行局部网格修正,获得最终计算所需的动网格。所提方法与直接弹簧法比较,由于弹簧法作用的背景网格远粗于实际的计算网格,提高了计算效率,此外,由于网格的运动受控于受弹簧法作用的背景网格,边界变形或运动能较为均匀地传递到内部网格上,使网格单元过渡光滑,提高了动网格的质量。给出若干算例,展示了方法的效率和可行性。     

基于Delaunay背景网格插值和局部网格重构的变形体动态混合网格生成技术

建立了一种基于Delaunay背景网格插值方法和局部网格重构方法相结合的变形体动态混合网格生成方法.首先利用作者发展的定态混合网格生成技术生成变形体初始网格,即变形体附近利用结构化的四边形网格,外场采用自适应Cartesian网格,中间由三角形网格过渡.当物体运动或变形时, 首先利用Liu和Qin提出的基于Delaunay背景网格的插值方法对变形体附近的网格进行变形;当物体运动或变形位移很大导致局部网格质量急剧下降、甚至网格相交时,则在局部重新生成网格.利用该方法生成了多种变形体外形的动态混合网格,如单个鱼体、双鱼串列巡游,昆虫单翼、双翼扑动,多段翼型变形等,并在鱼体巡游的非定常数值模拟中得到了应用.     

基于径向基函数与混合背景网格的动态网格变形方法

发展了一种基于径向基函数(RBF)插值和混合背景网格映射的动态网格变形方法。混合背景网格由运动物面附近的各向异性网格和偏离运动物面的各向同性网格组成。物面变形或运动后,包含在各向异性背景网格中的所有物面网格点的位移由指定的物面运动规律得到,而各向同性背景网格点的位移通过径向基函数插值获得。然后变形后的计算网格利用计算网格与混合背景网格之间固定的映射关系代数插值得到。各向同性网格距离物面一定的距离,消除了物面保形对径向基函数插值的要求,从而可以减少径向基函数基点的数量,以提高插值的效率和可靠性。最后,通过二维和三维算例对该方法的变形网格质量、参数影响、网格拓扑依赖性、复杂外形变形能力及网格变形效率进行了测试。测试结果表明,该方法能够高效生成大变形下高质量的变形网格。     

基于网格缩小的非结构网格梯度重构及制造解精度测试与验证

网格收敛性研究是验证与确认中的重要手段,本文针对一般非结构网格难以自相似加密的问题,实现了基于网格缩小的精度测试方法。在典型的各向同性和各向异性拉伸网格上,利用网格缩小精度测试方法分别考察了梯度重构精度以及制造解流动模拟精度,将网格缩小精度测试结果与传统的网格加密精度测试结果以及理论分析结果进行对比,验证了网格缩小精度测试方法与网格加密精度测试方法及理论分析的一致性。最后将网格缩小精度测试方法应用到各向异性拉伸网格粘性制造解精度测试中,得到了预期中的结果,并初步考察了梯度重构方法、网格类型对数值模拟精度的影响,显示出网格缩小精度测试方法在一般非结构网格精度测试上的优势,具有一定应用前景。     

基于非结构动网格的非定常流数值模拟方法

给出了基于非结构动网格的非定常流动问题的数值模拟方法,重点说明两种动态非结构网格方法的原理以及在非定常流模拟中的算法.一种是基于非结构网格的重叠网格方法,该方法结合了非结构网格方法和重叠网格方法的优点,以物面距为准则参数对网格节点进行分类,进而建立重叠网格的网格问插值关系,通过采用动态非结构重叠网格算法,可以处理多体间...     

基于径向基函数改进的Delaunay图映射动网格方法

为了提高大变形下动网格生成的效率和质量,基于径向基函数插值技术对原始Delaunay图映射动网格方法进行了改进。首先通过带多项式的径向基函数插值方法计算出背景网格远场边界点的位移;然后将背景网格位移插值到计算网格;最后利用衰减函数将计算网格远场位移衰减为零。基于矩形旋转网格变形实例比较了改进方法与原始Delaunay图映射动网格方法之间的差别,并研究了径向基基点数目和衰减函数类型对变形网格质量的影响。矩形旋转网格变形实例说明该方法能够精确恢复出Delaunay背景网格的旋转特性。NACA 0012翼型、NLR 7301两段翼和M6机翼网格变形算例进一步证明,通过添加控制点,该方法能够不重构背景网格实现大变形下高质量动网格的生成。     

基于嵌套网格的CHN-T1标模数值模拟

为响应第一届航空CFD可信度研讨会(AeCW-1),对CHN-T1标模翼身-尾翼组合体构型采用两种RANS求解器进行了基于结构化嵌套网格的网格收敛性研究和抖振特性计算,并与对接网格结果进行了对比分析。参考AeCW-1网格生成指南自主研发了一族以约3.3倍规模增长的稀-中等-密-特密嵌套网格系列。相比于分区点对接网格,融合了贴体网格和笛卡尔网格优势的嵌套网格从根本上降低了网格生成难度,单块网格拓扑形式更为合理且正交性较好,同时能够有效地平衡近场和远场的网格量。对于跨音速运输机构型,同等规模下嵌套网格的物面网格密度较对接网格更大,能够对激波和分离等复杂流动进行更理想的预测。对比计算过程中,自研求解器OFS3D表现出了较高的嵌套网格计算效率和可靠性,而NASA CFL3D求解器的嵌套网格计算能力有待进一步挖掘和验证。计算还发现,不同湍流模型得到的表面压力分布、分离区大小以及最大升力系数存在差别,而风洞模型支撑和机翼静气弹变形对飞行器力矩特性预测影响较大。     

一种基于约束框架的棱柱网格生成方法

提出了一种基于约束框架的棱柱网格生成方法,即利用物面网格节点与背景结构框架网格之间的单一映射关系线性插值生成棱柱网格。背景结构框架网格基于TFI方法生成。物面网格节点与背景框架网格之间的映射关系通过一种线性投影方法建立。然后通过V型槽、圆球、M6机翼和F6翼身组合体四种不同外形对该棱柱网格生成方法进行了测试,测试结果表明该方法能够避免棱柱网格生成中的网格交叉现象,保证棱柱网格的成功生成。该方法在棱柱网格质量控制、棱柱网格局部修改等方面具有较好的效率和优势,但背景框架网格的生成依然是该方法的一个瓶颈,需要进一步的研究。     

用刚粘塑性有限元法分析叶片精锻时的金属流动和变形组织

本文依据现行生产的某航空发动机压气机二级叶片的有关数据,采用刚粘塑性有限元法对叶身截面的成形过程进行了模拟,得到了纯铅和TC4钛合金两种材料在不同工艺条件下的金属流动和应变分布的规律,讨论了摩擦条件和毛边桥尺寸对金属流动的影响比较了锻件内不同区域的应变量所对应的组织。计算网格与实验网格比较基本吻合。     

斜齿轮齿面柔度矩阵与修形的有限元计算

提出了一套斜齿轮修形用的有限元计算程序。该程序自动生成轮齿有限元网格和齿面接触线, 计算得到接触线上结点的柔度矩阵;然后考虑支承结构与齿轮的综合变形, 由数学规划获得最佳齿面修形。计算与实验结果对比证明该程序操作简便、计算精度高, 具有工程应用价值。      

一种层压复合材料组合界面单元及有限元模型

 为了高效模拟计算复合材料层合板的层间损伤,提出了一种由刚性元和零厚度的内聚力单元组合而成的新型组合界面单元,并提出了以弯曲板单元作为子层单元,与上述界面单元结合构成的层合板三维有限元模型,通过界面单元刚度的连续衰减模拟分层损伤累积过程。有限厚度的八结点组合界面单元每个结点有5个自由度,能考虑子层结点平动和转动对层间损伤的作用。通过对双悬臂梁(DCB)和末端缺口弯曲(ENF)试验的计算模拟表明,新型组合界面单元和有限元模型能很好地模拟复合材料层合板的分层损伤过程,对不同的网格划分有很好的适应性。     

基于PCL的导弹吊挂强度分析的参数化方法

基于Patran/Nastran有限元分析软件的后台文件,采用PCL（Patran Comm and Language）对Patran/Nastran软件进行了二次开发,建立了用于导弹吊挂强度分析的参数化、自动化方法。提出的方法实现了吊挂结构强度分析过程中的参数化建模-自定义网格划分-材料边界选择-强度分析-结果查看等整个过程的自动化。方法的应用可以方便开展吊挂结构随意改变尺寸、材料参数、边界等的强度分析。并且方法对于任意结构在有限元分析中尺寸、网格、材料、边界等参数化建模及有限元分析均可适用。     

轴对称正挤压的刚塑性有限元模拟

 采用刚塑性有限元法,模拟具有锥底凹模和平底凹模的轴对称正挤压过程。利用自编程序对正挤压时的挤压力,应力场和应变场分布以及挤压件前端部的“双峰”现象等进行了研究。此外,还对轴对称正挤压时可能出现的中心缩孔和表面裂纹缺陷进行了模拟研究,得到了预测上述缺陷的曲线或数据。 在计算过程中,对刚塑性有限元计算技巧也作了一些改进,例如提出了处理八节点四边形单元奇异点和质点跟踪技术的新方法等。 试验结果表明,理论计算结果与试验结果相当一致。     

用有限元素法分析转子-轴承系统动力学

 本文采用有限元素法分析转子-轴承系统的动力学问题。给出了模拟转子-轴承系统的动力分析程序（FEPRS程序）,该系统由刚性盘、弹性分布参量有限转子元素及离散轴承组成。用该程序能计算任何复杂转子系统在各种进动状态下的进动频率,并由计算机自动绘制相应的振型曲线。     

复合材料层合板、壳结构的承载能力

本文以有限单元法研究了复合材料层合板、壳结构的承载能力,所采用的单元是考虑耦合效应的高精度的离散性克希霍夫（Kirchhoff）假定层合板和层合壳三角形单元,非线性迭代采用变步长的载荷增量法以及等刚度和变刚度相结合的交替迭代法,编制了计算机程序DDJ-FZJF,并以多例证明该程序和算法是有效的。     

有限元法在高速碰撞模拟中的应用

在试验数据分析的基础上,应用LS-DYNA三维动力有限元程序对高速弹丸撞击铝板的三维数值模拟。通过计算结果与实验结果的比较,证明采用合理的材料参数和物理模型有限元法能够直观地分析弹靶的详细相互作用过程,且两者在宏观物理特征上基本吻合。     

飞机外挂物离机过渡过程的有限元算法.

研究飞机外挂物在按给定轨迹下沉旋转的离机过渡过程中，流场非结构网格的调整和再生成方法以及过渡过程欧拉方程有限元数值解法。中指出外挂物运动时流场网格调整或再生成仅需在围绕在外挂物的一小区内进行以减少计算量，而且固定网格的欧拉方程的有限元数值解法也可适用于运动网格情况，只要将固定网格的守恒通量项Fk用运动网格的Ek来代替，并建立了Fk和Ek之间的关系。末给出了带外挂物飞机的实例计算结果。