三维机翼的静气动弹性特性数值模拟研究

发展了一种计算流体动力学（CFD）和计算结构动力学（CSD）的耦合计算方法,对三维机翼的静气动弹性进行了数值模拟研究。采用三维欧拉方程为控制方程基于直角网格计算气动力,并耦合结构静平衡方程进行静气动弹性数值模拟,设计了CFD/CSD耦合计算的数据交换的方式。以M6机翼作为算例,进行了机翼静气动弹性的数值模拟,计算结果表明所发展的三维机翼静气动弹性数值模拟方法是合理可行的,并可作为机翼结构优化设计和考虑结构弹性变形影响的气动外形优化设计的基础。     

基于CFD技术的螺旋桨气动特性研究

通过求解绝对坐标系下的雷诺平均Navier-Stokes方程数值模拟螺旋桨侧流黏性绕流,计算螺旋桨非定常气动特性.为了保证时间精度和提高计算效率,采用含牛顿型LU-SGS(lower-upper symmetric-Gauss-Seidel)子迭代的全隐式双时间法并引入多重网格技术.为了描述桨叶的相对运动,采用运动嵌套网格技术.对于轴流状态,计算结果与实验值吻合良好;对于侧流状态,与工程算法1P(once-per-revolution)载荷法进行了比较,验证了方法的有效性.应用该方法计算螺旋桨在不同来流偏角下的周期性气动特性,可有效地用于螺旋桨的设计.      

基于DES的二维和三维空腔流动特性研究

空腔流动广泛存在于航空航天工程中,对其流动特性的研究具有十分重要的工程意义。采用数值方法对比模拟了二维和三维空腔流动,控制方程采用N—S方程,空间离散采用有限体积方法,对流通量计算采用Roe格式,非定常时间离散采用双时间步长方法,湍流粘性计算采用基于SA模型的DES方法。数值计算所得的三维空腔底面压强分布与实验结果一致,所得的二维和三维空腔内的流动结构与相关文献中的分析相吻合。将二维和三维的计算结果进行比较,发现空腔底部压强分布、空腔内的压强振荡、空腔内部流线图及声压级等均有所不同,得出在空腔宽度有限时,横向流动简化了流动结构．削弱了振荡幅值．但使振荡过程复杂化．其三维效应不能忽略．     

三维复杂目标求解的多层快速多极子方法

采用多层快速多极子方法(Multilevel fast multipole algorithm,MLFMA)求解混合场积分方程(Combined field integral equation,CFIE),并选择RWG型基函数,对金属带缝锥球体、三面角反射器以及钻石体的单站RCS(Radar cross section)进行了计算,计算结果与试验吻合良好.在此基础上计算了F-22缩比模型的单站RCS,其计算量、存储量分别达到O(NlogN)量级和O(N)量级,此方法适用于带有尖点和特别细长曲面的三维复杂目标,如战斗机外形的RCS计算分析.     

用预处理方法及可压Navier-Stokes方程计算低马赫数流动

用可压方法及预处理方法求解了低马赫数紊流流动,其中空间离散格式采用了Roe格式,紊流模型采用了DES模型.以绕NACA0012的紊流流动为例,分析了该方法对提高流动速度和精度的特性,并将其用到多段翼型及绕圆柱的非定常流动中.数值模拟结果表明,该方法能有效求解低速流动及混合流场问题,且具有较高的数值精度.     

具有参数激励约简的扰动KdV方程的混沌行为

研究了具有参数激励约简的扰动KdV方程的混沌动力学行为.利用改进的Melnikov方法分析了由于同宿轨道的横截相交而产生的混沌行为.对周期外激励、周期线性参数激励和周期非线性参数激励下的扰动KdV方程的混沌行为进行了比较,发现划分混沌区与非混沌区的临界曲线是互不相同的.尤其是对非线性参数激励系统,存在"死频率".当这类系统受到该频率激励时,不论激励的振幅多大,混沌也不会发生.用时间积分法对上述系统进行了数值计算,结果与理论分析一致.     

含金属芯压电陶瓷纤维的本构方程

含金属芯压电陶瓷纤维是一种新型的压电器件。为了研究压电纤维的压电和机械性能,本文推导了其本构方程。该方程用于描述压电纤维在受到外力和电压作用时的响应。在悬臂杆式压电纤维的电极上施加电压,并在自由端施加力,计算其在热平衡状态下的内能密度,并在整根纤维上积分,得到其励参量是作用力和施加的电压,而响应参量是自由端的伸长量和电极上的电荷量,两者之间由一个2×2的压电矩阵相联整体能量。结果表明,自由端施加力作用时的广义位移为伸长量δ,电压的广义位移为电荷量Q。所建立的本构方程中的激系。     

超声速进气道与弹体一体化外形研究

通过数值求解N S方程和风洞测力试验,研究了超声速进气道的布局、剖面形状及几何参数对全弹气动力和进气道内流的影响。结果表明:细长旋成体弹身背风面和侧边不宜安放进气道;矩形进气道能够产生涡升力,升阻比明显高于半圆形进气道,而且内部流动接近二维流动,流场畸变及流动损失情况均比半圆形进气道要好;矩形进气道横截面积增加,对升力系数影响不大,但阻力系数增加明显;其高宽比增加,升力系数增大,阻力系数减小。     

直升机旋翼桨-涡干扰脉冲噪声传播特性研究

建立了一个基于Navier-Stokes方程和自由尾迹模型的高效耦合CFD方法,用于旋翼桨-涡干扰气动和噪声特性的研究。该方法将旋翼流场计算域分为两部分:在桨叶附近区域,通过求解Navier-Stokes方程来模拟流场中旋翼桨尖涡的形成;在尾迹输运区,采用自由尾迹模型表示尾迹涡的运动及其影响。噪声的计算采用基于声学类比法的FW-H方程。应用上述方法对AH-1/OLS旋翼桨-涡干扰状态进行了计算,通过对比脉冲噪声的声压时间历程,验证了方法的有效性。在此基础上,对桨-涡干扰噪声的空间传播特性进行了研究。计算结果表明:下降飞行状态的旋翼,在桨盘面会产生多处桨-涡干扰现象,桨叶片数增多,干扰也会明显增加;桨-涡干扰噪声具有较强的方向性,指向旋翼前行侧的前下方,其噪声声压级的衰减速率与距离成线性反比关系。     

后掠机翼横流驻波及其谐波研究

使用Mack方法求解三维不可压Orr-Sommerfeld方程(忽略曲率影响),得到模型在风洞最佳实验条件下空间最大放大率展向波长为4mm.以此为基础,选择大于此波长的具有代表性的展向扰动波长6mm与7mm,计算了这两个驻波及其谐波空间的放大率.研究在整个层流区主导流动的扰动波长及其发展趋势.使用丝网印刷的方式,在后掠机翼前缘添加粗糙带,人工引入6mm与7mm驻波,通过升华法验证了计算与理论分析.