超静超稳卫星碰振动力学建模

针对分离式卫星载荷模块(PM)受到扰动时可能与服务模块(SM)发生碰撞的问题，综合音圈电机反电动势(back-EMF)和柔性线缆动力学的效应，基于牛顿欧拉法建立了分离式卫星(DFP)载荷模块动力学模型。基于Hertz接触理论，推导了分离式卫星碰撞过程中连续接触力模型，并分析了碰撞过程中产生的接触力对载荷模块指向精度和指向稳定度的影响。数值仿真结果表明，碰撞使得载荷模块指向精度和指向稳定度下降5个数量级，碰撞后载荷模块可再次恢复到超静超稳工作状态，恢复时间超过1400 s。本文建立的碰撞模型对研究分离式卫星碰撞规避和碰撞控制具有重要意义。     

某飞行器推力矢量试验测力装置研制

为满足某飞行器推力矢量试验的测试要求,采用两台六分量应变天平和一个空气桥组成测力装置来实现飞机全机气动力和喷管推进特性同时分别测量。基于有限元软件,对天平的应变及空气桥对天平的干扰进行了分析。结果表明:空气桥对天平力分量的干扰值优于5%,对力矩(滚转力矩除外)分量的干扰值优于15%,达到了设计指标。通过校准获得了单独天平、天平带空气桥(充压、不充压)状态下的校准公式,校准结果表明:两台天平各分量的综合加载误差均优于0.3%,天平带空气桥(充压、不充压)状态下各分量的综合加载误差优于0.5%,空气桥对天平的干扰量值与有限元分析结果一致。理论分析与实测结果证明:研制的天平及空气桥达到了预定目标,它的测量精准度高,满足推力矢量风洞试验需求。     

遗传算法在机翼气动优化设计中的应用

把最优化方法与机翼的气动力求解相结合,进行跨音速机翼的气动优化设计。采用最优化方法为遗传算法,机翼的气动力由三维欧拉方程的数值解来提供。与基准机翼相比较,优化设计的机翼其气动性能有较大幅度的提高,表明遗传算法在机翼优化设计中的可行性和有效性     

翼-身组合体绕流的Euler方程数值模拟

用代数方法生成翼-身组合体H-O型网格,并用有限体积法研制出翼-身组合体绕流三维Euler方程计算程序。该方法的特点是改进了机身与机翼表面网格点分布,机翼后缘有后掠时也能保证后缘与网格线一致。程序除能按常规提供横流截面展向压强分布外,还能提供弦向压强分布。对NASA TND-712的翼-身组合体模型的计算结果与实验符合很好。     

高精度差分求解气动方程的几个问题

本文探讨了发展高精度格式的必要性，研究了高精度格式与网格以及熵增条件的关系，并发展一种半离散化的空间为三阶精度的格式。从此半离散化的格式出发，可建立多步显式格式及隐式格式。模型问题的计算表明，该三阶精度的格式具有好的精度，且激波附近基本上没有波动。     

TA17钛合金板状连接件的疲劳裂纹扩展行为研究

TA17钛合金是航空航天等工程中重要的结构材料,其疲劳裂纹扩展性能直接影响整体结构的安全性和完整性。本文基于Cr2Ni2MoV钢材紧凑拉伸试样斜裂纹疲劳裂纹扩展的试验结果,验证了ABAQUS软件XFEM模块用于分析平面应力状态下疲劳裂纹扩展的有效性,然后采用XFEM模块分析了TA17钛合金及其连接件的疲劳裂纹扩展性能。结果表明降低最大荷载、减小加载比和采用长圆孔型均可有效提高TA17钛合金连接件的疲劳寿命。采用XFEM模块可以预测连接件疲劳寿命的S-N曲线,预测TA17钛合金结构载荷增大33.3%,疲劳寿命减小为原来的1/4,可为工程疲劳设计提供参考依据。     

求解多维欧拉方程的二阶旋转输运格式

提出了一个基于旋转近似Riemann求解器的二阶精度迎风型有限体积方法.不同于"网格相关"(Gridaligned)的有限体积方法或者维数分裂的有限差分方法,本格式在求解Riemann问题时不依赖网格的方向,而是沿具有一定的物理意义的两个方向(称为迎风方向).我们发现当迎风方向取为控制体界面两侧速度差矢量方向(及与之正交的方向)时,该格式能够完全消除基于Riemann求解器的通量差分裂格式存在的激波不稳定或者所谓"红斑"(carbuncle)现象.为了提高格式的时间和空间精度,我们通过在控制体界面处求解线化的广义Riemann问题的方法体现输运过程对通量计算的影响.这种方案,使得我们有可能以此为基础,构造真正多维的有限体积型迎风格式.     

AUSM+-up格式在无网格算法中的推广

将AUSM -up格式推广运用到了无网格算法中.在点云离散的基础上,运用曲面拟合和AUSM -up格式求得数值通量;采用四步龙格-库塔方法进行时间推进,并引入当地时间步长和残值光顺等加速收敛措施.为了验证本文方法的计算精度和鲁棒性,对NACA0012翼型跨声速流动、某三段翼型低速绕流、平行错位NACA0012双翼超声速、高超声速流动进行了数值模拟.     

多重网格法求解二元跨音速欧拉流

 近年来,采用时间相关欧拉方程或N-S方程求解跨音速流得到了迅速发展。为了获得二元或三元复杂形体正确的表面气动参数和足够高分辨率的流场物理特性,必须采用十分细密的网格,因而导致计算时间成倍增长。为此,除继续发展各种差分格式,探讨各种自适应网格技术外,运用多重网格技术加速收敛过程得到了国外普遍关注并取得了一定的成功。国内也已有采用多重网格法解跨音速欧拉流的文章,然而均局限于两重网格的计算。为探索多重网格法和焓阻尼技术对加速收敛的作用,我们以NACA0012翼型的跨音速流动为例,采用MacCormack两步显式格式进行了数值实验。     

矩的显式积分算法研究与应用

矩的求解通常被用于求解有限元、体积、惯性矩等问题中.基于矩的叠加性,首先给出了在三维空间中计算域的离散方式,并推导了矩的显式积分公式,随后将其推广到n维空间中,该表达式易于在计算机上实现;设计了矩的并行计算算法,并通过Fortran和Python混编的方式,实现了矩的并行计算;对多重精度下的样例数据给出了一个算例,实现了零阶矩和二阶矩的计算,并和串行算法、逐次降维算法作出比较,进行了效率分析和误差分析.结果显示,矩的显式积分并行计算算法易于程序实现,并且在效率上高于串行算法,能够很容易推广到高维空间,该算法具有高度可并行性,误差主要来自计算域离散.