非结构网格上欧拉方程的自适应逆风隐式有限元解法

本文发展了一种解二维欧拉方程的隐式逆风有限元格式。空间近似应用的是二阶精度逆风方法；时间近似上根据数值通量函数的线性化处理导出了有限元的隐格式。为减少数值振荡，引入了基于特征变量的限制因子。在经自适应处理的非结构网络上，本文对跨、超、高超声速绕流进行了数值模拟。计算实践表明发展的隐式有限元格式与显格式相比具有较高的计算效率和较强的适应性。     

非结构网格生成技术的研究

本文描述了非结构网格生成的过程。在生成方法上综合Ｄｅｌａｕｎｅｙ和Ａｄｖａｎｃｉｎｇ Ｆｒｏｎｔ两种方法之长，构造了新的网格生成方法，使非结构网格生成更方便、快捷。网格生成过程分三个部分：１．表面网格的构造和最佳逼近；２．辅助点的生成与确定；３．空间四面体网格的生成。使用上述方法，本文构造限复杂构形物体的非结构网格。     

飞机蒙皮机器人化制造曲面分片算法研究

为实现飞机曲面机器人加工的自动轨迹规划,研究了基于模型几何特征的飞机曲面分片算法。通过分析飞机构件加工的工艺参数信息,以表面加工允许的最大曲率或最大偏角为阈值,提出基于近似曲率的三角网格模型分片方法。以三角网格中任意面片为起点,沿网格邻接关系向外拓展联接,形成若干近似平面化的模型分片结构,当超过阈值时则停止拓展,从而得到一整片曲面。通过对飞机CAD模型进行分片,发现采用本文提出方法可以保证分片结果的均匀性,为后期加工轨迹自动规划规划提供保障。     

立式捏合机桨叶结构参数对混合釜流场影响的仿真分析

为揭示立式捏合机桨叶几何结构参数与混合釜流场之间的关系,以1 L两桨立式捏合机为研究对象,采用计算流体力学ANSYS Fluent 14.5软件对立式捏合机进行了三维模型仿真研究,分析了桨叶结构参数(桨桨间隙、桨叶螺旋角)的改变对混合物料挤压应力和剪切应力的影响。结果表明,增大桨桨间隙可减小桨叶对混合物料的挤压应力和剪切应力,空心桨螺旋角分别为40°、45°、50°,桨桨间隙由1.0 mm增加到3.0 mm,挤压应力、剪切应力分别平均减小了82.4%、57.23%;增大桨叶螺旋角可增大桨叶对混合物料的挤压应力和剪切应力,空心桨间隙分别为1.0、2.0 mm,桨叶螺旋角由35°增加到55°,挤压应力、剪切应力分别平均增加了92.8%、55.7%。     

无网格方法关于导数计算的程序验证

为了提高无网格方法对空间导数的数值计算能力,提出并论证了新的适用于直角坐标系、柱坐标系和球坐标系的空间导数无网格算法。基于数值解与理论解的对比分析,对数值误差和收敛性进行了后验评估。评估结果表明:该算法对函数和导数的估计均为2阶精度,支撑域尺度对数值误差的大小有一定影响,但不影响数值计算的精度等级。在离散尺度h为0.01的条件下,对所选函数及其导数的数值计算相对误差不大于0.65%。      

基于非结构动网格的跨声速非定常流场求解

对用于非结构动网格生成的弹簧近似方法进行了研究.通过采用顶点弹簧方法,分析研究了弹簧倔强系数的取值,同时通过引入防扭转、防挤压倔强系数和边界修正,对标准弹簧近似方法进行了改进.采用改进后的方法并通过耦合求解了基于ALE(Arbitrary Lagrangian-Eulerian)描述的Euler方程,模拟了做谐和振动的刚性无后掠机翼及某翼身组合体的跨声速绕流.结果表明,改进后的方法可以大大提高网格变形能力和网格质量,且计算结果与实验结果吻合良好.     

火箭助推器从芯级飞行器动态分离过程的数值模拟

利用弹簧近似和网格重构相结合的非结构动网格技术，耦合求解Euler方程及弹道方程，时间方向采用四步Runge-Kutta方法，空间方向采用改进Barth和Jespersen限制器的通量分裂方法，数值模拟火箭助推器从芯级飞行器动态分离动力学过程。首先，计算单独芯级飞行器流场，与实验数据相比，符合较好；其次，计算火箭助推器和芯级飞行器组合体流场，得到分离前状态和气动力特性；在此基础上，比较采用弹簧和火箭作为控制力的两种分离方案，研究两侧火箭助推器分离不同步、攻角、侧滑角等因素的影响。研究表明，弹簧分离初期火箭助推器位移和姿态主要取决于弹簧控制力，弹簧全部断裂后气动力的影响加快姿态发散，在给定的设计参数条件下，可以实现安全分离；火箭分离存在复杂的喷流干扰，喷流对助推器的包裹作用使得分离初期自由来流影响较小；另外，分离过程对芯级飞行器的气动干扰不容忽视。     

发动机喷管内流场对流换热系数影响因素的数值分析

基于Roe-FDS格式,采用标准κ-e湍流模型对某液氧煤油发动机喷管内部燃气流对喷管内壁的对流传热换热系数进行了数值仿真与研究,分析了网格雷诺数以及恒温壁壁温对喷管内壁对流换热系数的影响.当近壁面网格第一层法向高度为10-3 mm数量级时,仿真结果与工程估算公式的计算结果相吻合.同时,在一定范围内,对流换热系数随来流速...     

固体姿轨控针栓推力器内流场数值仿真

为研究针栓推力器内.流场的特征和规律,采用二维轴对称模型,开展了针栓运动过程的非稳态及稳态流场仿真,分析了不同开度下的流场结构、针栓所受气动力及流场对针栓几何参数的敏感性.结果表明:由于针栓的干扰作用,流场结构呈现出复杂激波系、顶端流动分离等现象;气动力系数与针栓头部锥角大小呈正相关关系,与集气室压强无关;节流面形状对气流参数分布、激波强度及位置、气动力等表征的流动模态影响很大.     

变构型航天器参数化前处理方法

针对变构型航天器结构不固定的问题,基于VC 6.0软件对Trail及Session进程文件开发的方式,研究了利用Pro/E及Patran软件对航天器进行参数化建模及网格划分的具体操作方法,通过IGES中间文件的形式,将Pro/E及Patran相连接,利用VC开发工具建立了一套变构型航天器参数化前处理平台。采用两舱构型航天器对参数化前处理方法及平台进行了验证,结果表明:该参数化建模及网格划分方法切实可行,可为变构型航天器的静力学、动力学、热、结构等分析过程的前处理方法提供参考。