货运动车组装载门凹陷对隧道内会车压力波的影响

随着我国技术经济水平的提高，电子商务产业快速发展，货运动车组应运而生。采用有限体积方法和重叠网格方法，以及SST k-ω高雷诺数湍流模型，数值模拟研究时速350 km的货运动车组在隧道内交会的压力波，探究货运动车组车体结构即装载门凹陷对隧道内交会压力波的影响。实车试验结果验证了本文数值模拟方法的准确性。研究结果表明，凹陷的装载门使平直车身处的最大正压值增大约20%~30%，且有无凹陷的装载门中心最大正压值的差异百分比维持在35%~80%，最大负压值差异百分比在20%~25%。     

基于遗传算法的弹性地基加肋板肋梁无网格优化分析

基于遗传算法及一阶剪切理论,提出一种弹性地基上加肋板肋条位置优化的无网格方法.首先,通过一系列点来离散平板及肋条,并用弹簧模拟弹性地基,从而得到加肋板的无网格模型;其次,基于一阶剪切理论及移动最小二乘近似原理导出位移场,求出弹性地基加肋板总势能;再次,根据哈密顿原理导出结构的弯曲控制方程,并通过完全转换法处理边界条件;最后,引入遗传算法和改进遗传算法,以肋条的位置为设计变量、弹性地基板的中心点挠度最小值为目标函数,对肋条位置进行优化达到地基板控制点挠度最小的目的.以不同参数、载荷布置形式的弹性地基加肋板为例,与ABAQUS有限元结果及文献解进行比较.研究表明,采用所提出的无网格模型可有效求解弹性地基上加肋板弯曲问题,结果易收敛,同时基于遗传算法与改进混合遗传算法所提出的无网格优化方法均可有效优化弹性地基加肋板肋条位置,后者计算效率相对较高,只进行了三次迭代便可获得稳定的最优解,此外在优化过程中肋条位置改变时只需要重新计算位移转换矩阵,又避免了网格重构.     

《数值流形法》序

得知郑宏教授的《数值流形法》即将脱稿,激动之情难以言表.我完全同意郑教授将流形法的特点归结为“升阶”和“切割”.受郑教授的委托,我愿在此将建立流形法的历程做一简单的回忆.其实早在1985年———也就是在我完成不连续变形分析法(DDA)之前,“升阶”在我心中已完全成熟.“升阶”似乎能取得一致收敛,弥补有限元只能按能量收敛的缺憾.我没有急于发表是因为一直没能想到令我满意的解决“不连续”的办法.起初曾想过在DDA块体上引入高阶多项式,或在块体上再划网格,也就是用分片多项式来提高DDA块体的变形精度.但我对这些方案都不甚满意,因为实际计算时多项式的阶数会受到限制,引入网格又会招致有限元既有的毛病——解的精度受控于网格质量.     

基于六面体网格的向量式有限元分析及应用

基于向量式有限元基本原理，给出了八节点六面体等参实体单元的基本公式，通过投影方式将空间曲面六面体转换为投影六面体，采用参考面的逆向运动求解节点纯变形，通过单元形函的虚功方程计算节点内力；针对坐标模式和内力积分模式等关键问题提出了有效的处理方案。编制了六面体实体单元的数值计算程序，并进行工程结构算例分析。结果表明，所编制程序可有效模拟实体结构的静力、动力及大变形大位移行为分析，验证了本文理论和程序的正确性和实用性。     

三维非稳定渗流自由面边界积分项的精确数值计算

利用固定网格法分析三维非稳定渗流问题时,将要面对两项积分难题:以自由面及单元表面为边界的空间积分及以自由面为边界的曲面积分。针对常用的任意8结点6平面三维普通单元,提出采用坐标变换及等参变换技术求取空间积分项的精确数值解;至于曲面积分项,建议改用单元非饱和区部分表面作为积分边界,经过坐标变换及等参变换处理积分边界后,利用高斯数值积分可求出曲面积分项的精确数值解。通过一个普通单元及一项均质半无限边界堤坝的实例分析,表明此方法的精确性和稳定性良好。     

三维无网格Galerkin结构分析方法及其应用研究

针对无网格Galerkin法在三维复杂几何形状的结构分析中存在的刚度矩阵的稀疏存储实现难、六面体背景网格适应性差等问题,本文采用逐节点对法组装刚度矩阵,利用CSR格式存储刚度矩阵,提出了一种基于四面体背景积分的改进的三维无网格Galerkin法。通过采用罚函数法施加对称约束和给定位移值约束,并推导出了施加这两种位移约束的统一格式。利用所提算法完成了三维悬臂梁的计算,所得结果与其理论解相吻合;完成了轴流式风机轮毂的结构分析,得到的位移与应力分布结果与其有限元解相吻合。这表明本文所提方法能满足工程应用中的计算要求,并适用于具有复杂形状的几何模型分析。     

等径角挤压后Ti5553钛合金的冲蚀磨损机理演变

本文中系统研究了经等径角挤压(ECAE)处理后的Ti5553钛合金在海砂环境中的冲蚀磨损机理演变历程和失效原因.发现Ti5553钛合金在高角度冲蚀工况下,其冲蚀磨损机理从初期的微切削转变为冲击挤压变形,并保持稳定;合金显微组织中弥散的α相以及ECAE工艺影响了冲蚀机理演变过程及冲蚀磨损程度;提高合金的强度、韧塑性对于抵抗冲蚀磨损起到了促进作用.     

HyperFLOW软件非结构网格亚跨声速湍流模拟的验证与确认

计算流体力学(computational fluid dynamics,CFD)数值模拟在航空航天等领域发挥越来越重要的作用,然而CFD数值模拟结果的可信度仍然需要通过不断地验证与确认来提高.本文给出了从制造解精度测试、简单到复杂外形湍流模拟网格收敛性研究等三个方面开展CFD软件验证与确认的方法,并对自主研发的CFD软件平台HyperFLOW在非结构网格上模拟亚跨声速湍流问题的能力进行了验证与确认.首先通过基于Euler方程和标量扩散方程的制造解精度测试,分别验证了HyperFLOW在非结构网格上对Euler方程和黏性项的求解精度,结果表明其能够在任意非结构网格上达到设计的二阶精度. 其次,通过NASATurbulence Modeling Resource中的湍流平板、二维翼型近尾迹流动、二维Bump等几个典型的亚声速湍流算例的网格收敛性研究,量化考察了数值结果的观测精度阶和网格收敛性指数,并与国外知名CFD解算器CFL3D,FUN3D的计算结果进行了对比,验证了HyperFLOW对简单湍流问题的模拟能力,且具有良好的网格收敛性和计算精度(阶). 最后,通过NASA CommonResearchModel标模定升力系数的网格收敛性研究和升阻极曲线预测,验证了软件在复杂外形亚跨声速湍流流动数值模拟中也具有良好的可信度.      

旋转射弹高速倾斜入水多相流场与弹道数值模拟

基于VOF多相流模型和有限体积法求解水、汽、气多相流动的RANS方程，结合重叠网格技术和six DOF算法对某一型号舰载射弹倾斜入水过程进行数值模拟研究。首先基于该方法研究了射弹旋转效应对射弹运动特性及流体动力特性的影响，然后对不同入水角下倾斜入水过程进行分析，得到不同倾角下旋转射弹入水空泡形态发展规律、弹体运动特征及流体动力特性变化规律。研究结果表明:射弹的旋转有利于弹体在初始对称面内的弹道稳定性，但会降低弹体侧向稳定性，使射弹受到的阻力系数、俯仰力矩系数变小；入水角越小，形成的空泡越不对称，由射弹运动状态的改变引起的空泡形态变化越明显，在超空泡航行阶段，弹体运动较稳定，不同角度下流体动力系数差别很小，当弹体下表面刺破空泡壁沾湿时，弹体运动状态发生较大变化，流体动力系数迅速增大，此时入水角度过小，弹体容易失稳；弹体的沾湿对空泡形态、弹体运动稳定性和流体动力特性有着重要的影响。