高效计算网格壁面距离的KD树方法

计算流体动力学的一些领域中经常要用到流场点到壁面的最小距离,为了提高壁面距离计算的精度和效率,将KD树应用于非结构网格的壁面距离计算,计算了若干三维外形的壁面距离,结果表明:基于KD树的网格壁面距离计算方法鲁棒性好,计算效率和结果精度高,适应复杂外形的能力强;算法的通用性好,可以拓展应用到多种网格类型。     

基于新型非结构有限体积解耦算法的激波诱导燃烧数值模拟

把基于结构网格有限差分方法建立的化学非平衡流动模拟的新型解耦算法推广到基于非结构网格的有限体积法,从而可用于模拟复杂几何构型下的化学反应流动.对H2/ Air预混气体中激波诱导振荡燃烧的Lehr实验进行了数值模拟,计算得到的振荡频率与实验结果符合很好,表明计算方法具有时间和空间二阶精度.通过对不同几何外形发射体的计算发...     

大功率船用减速器斜齿圆柱齿轮啮合仿真模型

利用大型有限元分析软件MSC.Marc的直接约束法,建立了某大功率船用减速器斜齿圆柱齿轮副多齿啮合的三维有限元非线性接触分析模型。基于该模型,对齿轮副进行了一个啮合周期内的准静态啮合仿真分析。计算表明,本文得到的结果非常符合实际啮合规律,也验证了模型的正确性。由于模型的特点,有望进一步应用于齿轮的动态啮合分析。     

复杂平面区域的三角网格生成算法

为生成复杂平面区域的有限元网格,提出了基于网格细化的三角网格生成算法.该算法首先采用耳尖移除法对区域边界做三角划分,得到粗略的初始网格.提出Delaunay优化平分方法,根据网格密度细化初始网格,该网格细化方法结合最长边平分技术与Delaunay边交换技术,可有效提高内点生成与单元细分的质量.实验表明,基于Delaunay 优化平分的三角网格生成算法可对任意平面域进行网格剖分,生成符合有限元计算要求的高质量三角网格.     

全隐LU-SGS算法在高超声速热化学非平衡流刚性问题中的应用

在高超声速热化学非平衡流动计算中,当地气体能量松弛时间、化学反应特征时间与流动时间推进步长量级差异过大会带来严重数值刚性问题,且在高雷诺数条件下,壁面、拐角等强干扰区网格加密使得该问题加剧,导致初始最大CFL数极小,收敛速度缓慢.原始LU-SGS算法仅考虑化学反应源项和对流项的隐式处理,通过推导黏性项Jacobian矩...     

面向复杂几何模型的多级并行四面体网格生成算法

高性能计算机的快速发展为数值模拟提供了必需的硬件环境,数值模拟领域对网格的需求已高达数十亿到数百亿量级,而网格生成作为数值模拟的重要一环,其发展则相对滞后,很难满足并行数值模拟求解器对大规模网格的需求.本文面向复杂几何模型提出一种多级并行四面体网格生成算法,该算法首先基于模型的几何特征建立网格的尺寸场,并基于尺寸场和几...     

某型运载火箭地震响应分析及减震设计

竖立状态的运载火箭一旦遭遇地震,极易造成箭体倾倒等灾难性后果。为了确保运载火箭在地震中的安全,针对某型运载火箭,采用时程分析法对其进行地震响应分析,得到了火箭的弯矩、位移响应峰值沿箭体纵向的分布规律。分析了火箭关键位置处弯矩响应的频域特性,设计了火箭的减震措施,并分析了火箭与发射台的连接刚度、隔震支座刚度对火箭地震响应的影响规律。结果表明:火箭弯矩响应的最大值发生在靠近火箭尾段的位置,位移响应的最大值发生在火箭头部位置;改变火箭与发射台连接刚度以及在发射台底部加装隔震支座的方式可以有效地减小火箭的地震响应。     

认知无线Mesh网络中资源节约的最小信道冲突值组播路由算法

认知无线Mesh网络中满足服务质量约束的多目标优化组播路由问题比单目标优化组播问题更加复杂,为了快速求解认知无线Mesh网络中满足服务质量约束的、以最小化资源消耗与最小化信道冲突值为目标的多目标优化组播路由问题,提出一种基于粒子群优化算法的问题求解框架,包括问题描述、粒子编码与粒子初始化、适应度函数、粒子飞行、粒子变异、粒子消环。粒子由表示节点之间连接关系的带权邻接矩阵表示,重新定义了用于粒子飞行的3种运算规则以及粒子飞行运算和粒子变异运算。仿真结果表明提出的算法能达到预定目标,资源消耗较低且能获得较低的信道冲突值。     

渐开线环形齿球齿轮传动时变刚度系数

为研究渐开线环形齿球齿轮传动过程中齿轮副间的接触力,基于Hertz弹性接触理论定义了齿面接触点的刚度系数;利用微分几何相关知识,分析了球齿轮的齿面特征,得到了实时啮合点主曲率半径计算公式;研究了重合度系数对啮合刚度的影响,结合球齿轮机构传动特点,推导出任意位置反向啮合点位置计算公式;最后通过一个实例说明球齿轮机构接触点刚度系数具有典型的时变特征.     

激光熔覆搭接过程残余应力应变的有限元计算

激光熔覆过程中产生的残余应力与应变对熔覆层的裂纹开裂倾向有重要影响。利用有限元法对激光多道搭接过程熔覆层的残余应力和应变进行计算,计算过程中主要考虑瞬时温度变化引起的热应力、σ-ε曲线、有限元网格、力学边界条件4个方面。热应力主要在温度场计算结果的基础上求得;σ-ε曲线采用的是线性强化材料的弹塑性曲线;有限元网格主要采用分区划分的形式进行。通过计算,得到了搭接后的熔覆层上关键点的残余应力和应变分别为1 400 MPa和2.2×10^-2。结果表明：搭接熔覆层交界处的点产生焊趾裂纹的倾向最大,第1道熔覆顶点处产生横向裂纹的倾向次之。