二维Euler方程的自适应非结构直角网格算法

本文用四叉树型非结构直角网格和有限体积法求解了二维Ｅｕｌｅｒ方程,在基本网格生成过程中,采用了以几何外形为基础的网格加细方法,包括物面边界网格加细和曲率加细;并提出了一个用自适应的方法在激波处及其附近局部加细基本网格过程,提高了激波分辨率,在计算中将结构网格经常采用的以中心差分为基础的Ｊａｍｅｓｏｎ的有限体积法做了推广,减少了计算工作量,节省了机时,为今后推广到三维创造条件,本文对单个翼型,多段翼     

复合材料缠绕后置处理方法研究

利用芯模表面的离散落纱点,依据空间几何和微分几何理论,针对纤维缠绕过程提出纤维缠绕出纱点轨迹的网格后置处理方法。该方法首先通过纤维缠绕轨迹的落纱点,以及缠绕过程中等悬纱长约束条件来求取出纱点,然后通过计算相邻出纱点的各坐标差得到机床运动轨迹,最后分别通过对三通管和组合回转体的缠绕成型实验,来验证该网格后置处理方法在处理非回转体和回转体缠绕出纱点轨迹的可行性和准确性。     

采用近似边界求解壁板颤振的CBS有限元算法

为了简化Euler方程求解壁板颤振问题,将近似边界条件成功地推广到CBS(characteristic-based split)有限元方法求解壁板颤振领域,把运动壁面的影响转化为Euler方程的边界条件,避免了对动网格的处理.通过对近似边界方法、动网格方法以及3阶活塞理论的数值模拟结果进行对比发现:首先,加载近似边界的CBS有限元方法操作非常简便,且模拟结果与动网格方法完全吻合;其次,Euler方程可以准确捕捉到壁板前后两端不光滑尖点所引起的压力突跳,而3阶活塞理论则无法做到;最后,采用近似边界的CBS有限元方法在求解无黏可压缩流动中细长体小变形的流-固耦合问题上具有潜力和优势.      

齿轮轮齿接触分析的分解算法

提出了齿轮轮齿接触分析算法——分解算法。传统的轮齿接触分析方法求啮合点时需要求解含5个非线性方程的方程组,求解性差;齿面接触和边缘接触的数学模型不同,需要分别进行求解,求解过程复杂。轮齿接触分析算法——分解算法,提出了瞬时共轭啮合线的概念,可有效分离传动误差,得到啮合点、瞬时接触线,求啮合点时非线性方程的个数由5个减少为2个。分解算法建立的数学模型也适用于边缘接触分析,算法简单、有效、适应性强。以一对弧齿锥齿轮为例, 对比分析了传统方法和分解算法, 结果表明: 齿面部分的印痕是一致的,传动误差幅值相差0.3″;边缘接触部分的印痕存在少许差异。      

翼身融合民机典型PRSEUS受压壁板屈曲及渐进损伤分析

由NASA和波音公司共同提出的拉挤杆缝合高效一体化结构（PRSEUS）,由于具有优异的抗压稳定性和止损/止裂等承载优势,已成为解决翼身融合布局民机非圆截面机身结构承载效率低和稳定性差等问题的主要途径。本文针对典型PRSEUS受压壁板结构,开展了线性/非线性屈曲及渐进损伤分析;提出了综合考虑蒙皮、止裂带、长桁翻边、隔框翻边等一体化缝合元件贯穿支撑构型几何关系和偏置参考面的建模方法,提高了PRSEUS受压壁板有限元模型的精度;提出了综合考虑屈曲特征值、非线性屈曲载荷等多影响因素的网格收敛性分析方法,提高了PRSEUS受压壁板屈曲分析的计算效率;提出了最小屈曲特征值、几何节点偏移以及最小屈曲特征值-几何节点偏移组合式等3种初始缺陷引入方法,提高了PRSEUS受压壁板损伤分析的计算精度;完成了基于纤维与基体损伤本构关系的典型PRSEUS受压壁板非线性屈曲损伤分析,通过与试验结果对比,给出了针对PRSEUS结构的非线性屈曲渐进损伤演化分析方法。为翼身融合布局民机PRSEUS结构的稳定性/损伤分析和设计提供了方法和技术支撑。     

风洞阻尼网周边拉力分析与计算

阻尼网是风洞稳定段内的重要设备,它可以进一步破碎蜂窝器后的气流漩涡,减小稳定段气流的湍流强度。在风洞运行时,阻尼网因受到气动载荷而变形,其变形大小与周边的拉力具有一定的关系,如果变形过大,将影响阻尼网的正常使用,起不到减小湍流强度的作用。以圆形阻尼网为例,推导出阻尼网周边的拉力计算公式及预紧力的选取方法,为风洞阻尼网的安装预紧提供了参考和借鉴。     

基于自适应网格加密的超声速可燃气热射流起爆详细反应数值模拟

采用块结构自适应网格加密开源程序AMROC,在高性能计算集群中,进行氢气、氧气、氮气详细反应机理的二维超声速热射流爆震起爆自适应网格精细数值模拟,研究在热射流的持续喷射作用下超声速可燃气热射流爆震起爆,以及形成的精细的爆震胞格结构.结果表明:超声速可燃气流中热射流类似气动斜劈,形成局部激波诱导燃烧.超声速可燃气中持续的热射流喷射会形成过驱爆震,并导致不规则爆震胞格的生成.热射流的扰动压缩作用对过驱爆震的形成以及不规则爆震胞格的产起着关键的作用.热射流的扰动压缩波以当地声速通过爆震波后的亚声速区域作用于爆震波,使爆震波处于持续过驱状态,并形成不规则的爆震胞格结构.     

高速飞行弹丸诱导的三维爆轰波数值模拟

采用三维非定常化学反应控制方程组对驻定在带攻角高速飞行弹丸上的斜爆轰波流场进行了数值模拟。对流项和化学反应项采用时间分裂法,分别用TVD格式和基元反应模型进行处理。计算结果表明三维数值模拟可以更准确的描述弹丸飞行攻角对爆轰波的形成及结构的影响,是进一步研究超驱爆轰推进技术的重要手段。     

网格简化中基于特征矩阵的二次误差测度算法

针对二次误差测度算法存在尖端特征消失、局部过度简化等缺陷,提出了基于特征矩阵的二次误差测度算法用于网格简化.通过将顶点曲率和边长引进该特征矩阵以优化误差度量,模型中各顶点便易于区分,于是具有明显几何特征区域的顶点误差度量能够被提高.这样,边折叠的顺序可以方便的得到调整,使得模型中的突出特征更多的被保留下来.仿真结果表明,本算法在保持了二次误差测度算法计算时间短、运行效率高的同时,也克服了网格分布过于均匀、无法突出模型重要特征的缺点.     

螺旋锥齿轮啮合刚度及参数振动稳定性研究

准确计算时变啮合刚度是齿轮系统动力学研究的基础.针对航空高速重载螺旋锥齿轮,基于轮齿接触分析(TCA)和轮齿加载接触分析(LTCA)通过计算瞬时接触点的轮齿变形柔度建立了时变啮合刚度数值模型;将齿轮时变啮合刚度在一个啮合周期内视为逐段线性,基于Floquet理论推导了含时变刚度参数振动系统的状态转换矩阵解析式;通过修正小轮机床调整参数设计三种接触情况,分析了算例齿轮在相同载荷工况下的接触轨迹、传动误差、重合度和时变啮合刚度;采用二自由度齿轮系统动力学模型考察工作转速范围内的周期运动不稳定区间,分析了时变啮合刚度对螺旋锥齿轮系统参数振动稳定性的影响.