着水初始条件对水陆两栖飞机着水性能的影响

针对着水初始条件对水陆两栖飞机静水面降落性能影响的问题,使用非定常雷诺平均Navier-Stokes方程和标准k-ω两方程模型求解非定常流场,耦合流体体积模型捕捉自由液面,结合改进的流场域网格划分策略,并借助任意拉格朗日欧拉方法处理机体、网格和水气交界面三者关系。基于上述方法,探讨初始姿态角、下降速度和前飞速度对着水性能的影响,其中驾驶舱与重心位置加速度、气动力与水动力、俯仰角与吃水深度等参数变化被考虑在内,逐步得到最优着水初始条件。结果表明：增大初始姿态角有利于减缓飞机各位置除初次触水过载峰值外过载峰值的变化。当飞机俯仰力矩出现极值时,驾驶舱垂向过载显著大于其他方向和位置过载。减小初始下降速度,各加速度值也随之减小,值得注意的是重心垂向加速度值与初始垂向下降速度平方值存在线性变化关系。着水前飞速度的选取应使得气动升力略大于重力,飞机才能获得更佳的着水性能。     

三维散乱点的c~1曲面插值及应用

本文介绍了一种对三维散乱点插值的二次Bernstein-Bezier c~1曲面构造方法。首先,不规则分布的3D数据点{(x_i,y_i,z_i),i=1,2,3,…,N}被投影到X-Y平面上,并按本文提出的能够处理任意区域内不规则分布点的三角化算法,自动形成平面三角形插值网络。然后按照所形成的三角形网络和网络结点处的函数值,分别估计出每一数据点上的一阶导数值。最后本文给出了用二次多项式表示的三角形网络上的Bernstein-Bezier c~1曲面插值公式,并指出了这一曲面插值模型在某些应用领域的广阔前景。     

基于顶点分类的曲面三角网格模型自适应光顺研究

逆向工程中由测量数据得到的三角网格模型往往含有大量的噪声、扰动及不规则三角片，需对其进行光顺处理，以满足后续处理的使用要求。本文首先提出了一种新的网格模型顶点法矢计算公式，该公式采用网格模型中三角片面积与顶角角度综合加权，可同时反映三角片面积与顶角角度对顶点法矢的影响。在此基础上，对网格模型顶点进行了分类处理，提出了能够反映网格顶点特征性质的顶点势概念。势为1的点为特征点，势为0的点为普通点。最后，提出了一种新的自适应曲面三角网格模型光顺方法，综合了普通拉普拉斯光顺法与平均曲率法的优点。在该方法中，顶点调整方向为拉普拉斯光顺矢量在被调整顶点切平面上的分量与该点法矢的加权合成；顶点调整幅度根据顶点类别的不同而不同，特征点的调整幅度小，普通点的调整幅度大，从而可保护原有特征。实例表明，与现有方法相比，该方法在有效去除噪声，匀化三角片的同时，充分保留了初始网格模型上的特征。     

非结构网格的高超声速粘性MHD流场数值模拟

应用数值方法在非结构网格上对磁场干扰下的二维高超声速钝头体粘性绕流进行了数值模拟。控制方程为N-S方程耦合Maxwell方程的粘性MHD方程组,空间离散采用有限体积方法,对流项用AUSM格式计算,粘性项用中心格式求解,时间推进用显式5步Runge-Kutta格式,引用双曲型散度清除技术加强▽.B=0的条件。计算模型为二维钝头体,在高超声速来流条件下,分别对有、无均匀磁场干扰下的流动进行了数值模拟。计算结果表明,在均匀磁场干扰下,激波脱体距离显著增加,物体表面压力急剧下降。对比表明文中发展的计算方法可以准确地进行二维粘性MHD流场的数值模拟。     

基于点簇合并的多边形表面模型简化算法研究

 在计算机图形的动态实时显示技术中, 将高度密集的多边形表面模型用尽可能少的多边形逼近成为近年研究的焦点, 也出现了一批简化算法。针对存在的问题, 文中以累进网格为基础, 并将点簇合并引起的几何逼近作为研究重点, 将多种度量方法通过权因子统一到基于二次误差的模型逼近度量方法中, 并提出了改进的基于二次误差的多边形表面简化算法。实例表明该算法可快速、高质的实现逼近。该算法也适用于实时环境下的多边形网格的动态多分辨率表示。     

非结构网格高超声速化学非平衡MHD流场数值模拟

 在非结构网格上对考虑化学非平衡效应的二维高超声速磁流体绕钝头体流动进行了数值模拟。控制方程由二维理想磁流体动力学（MHD）方程和组元连续方程两部分组成,化学动力学模型为5组元17反应模型。MHD方程空间离散采用AUSM格式,时间推进采用显式5步龙格-库塔格式,并通过弱耦合的方式与化学反应控制方程结合在一起。计算模型为二维钝头体,外加磁场为偶极子场,磁场源位于钝头体内部。在高超声速来流条件下,对有、无磁场干扰,是否考虑化学反应下的4种工况进行了数值计算,得到了满意的结果,并与有限的参考文献进行了对比。结果表明本文发展的方法能准确地模拟考虑化学非平衡效应的高超声速MHD流场。     

叶栅流场中的气动阻尼研究

从气动阻尼产生的机理出发,建立合理的气动阻尼模型;利用动网格技术,结合等效模态气动阻尼比的计算式,以Rotor37跨声叶栅为例,用弱耦合方法数值计算得到一弯、一扭与弯扭耦合模态的气动阻尼比,并考虑流场压力水平及叶片振幅对气动阻尼的影响.计算结果表明,在保证进出口压比不变的情况下,气动阻尼与压力成线性关系;在小振幅下,气动阻尼基本与叶片振幅无关.      

多重网格法在混合网格中的应用

设计了一种高效的混合网格聚合方法,在此基础上应用多重网格技术来加速Navier-Stokes方程的收敛过程.方程求解采用Jamson的中心格式有限体积法.对NACA0012翼型和ONERA M6机翼的粘性绕流进行了数值模拟,算例结果与试验结果符合良好,加速收敛效果显著.     

任意单元间断Galerkin有限元计算方法研究

基于龙格库塔间断Galerkin(RKDG)有限元法的构造思想,通过局部坐标变换,发展了非正交单元DG有限元计算方法;借鉴非结构网格有限体积隐式计算方法,发展了适应于DG有限元方法的隐式计算方法;借鉴一维和二维(三角形单元)DG有限元限制器构造方法,提出了非正交三棱柱单元限制器方法;利用上述方法数值模拟了球头和双椭球的高超声速粘性绕流,得到了清晰的流场结构,并得到了较好的压力和热流分布,表明该方法在复杂高超声速流动的数值模拟方面具有广阔的应用前景。     

双参数弹性地基上四边自由矩形薄板精确解

将弹性地基视为Vlazov双参数模型,利用二维有限域积分变换的方法推导出了Vlazov双参数弹性地基上四边自由矩形薄板在任意荷载作用下挠度和内力的精确解。在求解过程中,不需要预先人为选取挠度函数,而是直接从弹性薄板的基本方程出发。通过集中荷载计算实例验证了该方法的正确性。该方法计算简便,适用于不同边界的薄板问题求解。