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

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

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

为实现飞机曲面机器人加工的自动轨迹规划,研究了基于模型几何特征的飞机曲面分片算法。通过分析飞机构件加工的工艺参数信息,以表面加工允许的最大曲率或最大偏角为阈值,提出基于近似曲率的三角网格模型分片方法。以三角网格中任意面片为起点,沿网格邻接关系向外拓展联接,形成若干近似平面化的模型分片结构,当超过阈值时则停止拓展,从而得到一整片曲面。通过对飞机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%。      

基于节点改善策略的伪谱轨迹优化

基于Legendre Gauss Lobatto（LGL）伪谱法(PSM)求解最优控制问题的原理,研究了有限推力轨迹快速优化设计以及计算节点和效率改善策略。对远程变轨动力学模型进行了无量纲化处理,并设计了初值生成和串行优化求解策略,用来提高多变量多约束大规模稀疏非线性规划问题的收敛性。结合极小值原理推导了问题的最优性必要条件,基于协态映射原理设计了数值解的最优性验证方法。利用二阶微分方程技术消去部分状态量,建立了节点和计算效率改进模型,并讨论了该策略的相关数值处理技术。仿真结果表明,本文的算法可以快速地提供应用于实际飞行任务的最优解,节点改善策略在保证计算精度的同时,有效降低了计算收敛时间。     

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

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

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

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

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

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

铝合金盘型件等温辗压过程中晶粒尺寸演化的数值模拟与实验验证

 为了揭示盘型件双面辗压成形过程中微观组织的演化特征,应用数值模拟方法研究了这一成形过程中所发生的动态再结晶和晶粒长大规律.所研究的盘型件材料为铝合金6061,辗压成形温度为350~500 ℃.晶粒尺寸演化规律的数值模拟结果与实验结果基本符合.与整体锻造不同,双面辗压成形没有变形死区,并且越靠近盘型件表面晶粒细化效果越好;辗压头通过对盘型件厚度方向压缩和表面环向剪切两种方式驱动盘型件变形.计算结果表明,即使压下量零增量辗压也会在沟槽附近产生显著变形和晶粒细化,这说明盘型件变形和组织演化主要源于辗压头对盘型件表面的环向剪切;多遍次辗压过程中,第1遍次辗压对晶粒细化的贡献最大,后续辗压的贡献较小;提高成形温度可能会出现动态晶粒长大区,在这种情况下,盘型件上的点将交替发生动态再结晶和动态晶粒长大,最终晶粒尺寸会有一定程度增加.     

壁面函数在激波诱导分离流动中的应用

以数值模拟激波-附面层干扰引起的流动分离问题为研究背景,发展了基于有限体积方法的雷诺平均Navier—Stokes（RANS）方程的流场数值模拟方法。利用壁面函数模型得到壁面剪切应力,通过修正壁面粘性通量,构造了一种新的湍流边界处理方法,并将其耦合到RANS方程和SSTk-ω湍流模型的数值求解中;同时,针对激波诱导引起的附面层流动分离问题,提出一种附面层网格加密技术,能够自适应加密分离区内附面层网格,使得在流动分离区域也能够使用壁面函数模型。数值算例表明,壁面函数模型能够降低数值模拟结果对网格的依赖性;同时也验证了壁面函数耦合附面层网格自适应方法,在处理激波诱导引起的附面层流动分离问题时的有效性和准确性。