基于近似边界条件的三维机翼气动弹性数值模拟研究

物面边界条件采用一阶近似边界条件,计算网格采用静止的笛卡尔网格;通过求解Euler方程,得到三维机翼的非定常气动力,耦合机翼运动方程进行气动弹性的计算;流场和结构之间的数据通过无限平板样条法(IPS)进行交换.使用文中方法计算了三维气动弹性标准模型AGARD445.6机翼的颤振边界,计算结果与实验结果吻合,表明该方法可以准确高效地求解三维机翼的气动弹性问题.     

机翼-机身-尾翼结构的跨音速颤振分析研究

利用无限插值方法(TFI)生成三维多块贴体运动网格，以Navier-Stokes方程为控制方程，求解机翼-机身-尾翼结构的跨音速非定常气动力，并与颤振方程耦合迭代计算，求解飞行器广义位移响应的时间历程，根据广义位移的时间历程的衰减、等幅和发散振荡等情况确定飞行器跨音速颤振临界条件。为了提高计算效率，研制了以多台微机组成的分布式计算系统。开发了颤振分析并行软件。经算例验证，计算结果与实验结果和理论分析相吻合。     

一种流体-结构耦合计算问题的网格数据交换方法

气动／结构耦合数值模拟是研究非线性气动弹性的基础。数据交换和插值是非线性气动弹性仿真问题的关键。目前的插值方法不能满足非线性气动弹性问题。本提出了一种有限元四节点(FEFN)插值方法。该方法是一种局部插值方法，并不依赖于结构模型带来的整体信息。以圆柱体为具体算例，插值结果与有限平板插值方法(IPS)进行了算例对比，表明FEFN方法更能代表计算物体的表面，且计算简单、计算量小、误差小，是一种适合计算流体力学(CFD)／计算结构动力学(CSD)耦合仿真的界面数据交换工具。     

基于欧拉方程的一种机翼气动弹性计算方法

利用一种双时间方法求解三维非定常欧拉方程，采用无限插值理论生成O－H型代数网格，考虑了机翼变形时的网格生成问题，通过与气动力方程的联立求解，在时间域内用二阶龙－库塔方法求解机翼弹性运动方程。计算结果表明，本计算方法具有较高的计算效率，所计算的颤振临界速度与风洞实验一致。     

间隙非线性对二元翼段颤振特性影响

为了研究间隙非线性对颤振特性的影响,本文采用基于非定常雷诺平均方程的非定常气动力求解方法,耦合结构运动方程建立了时域气动弹性分析系统,并运用该系统计算三自由度无间隙二元翼段构型的颤振速度。采用描述函数法对间隙问题进行处理,得到了间隙非线性所带来的极限环振荡现象,并分析亚跨音速阶段间隙大小对颤振特性的影响。通过研究预加载对颤振特性的影响,得出预加载能够减弱间隙非线性影响,有效提高系统颤振速度。     

高超声速机翼瞬态气动加热下热颤振分析

研究了一种气动力(热)/结构耦合的高超声速机翼热颤振的时域数值分析方法。研究目的是分析结构内部温度梯度对结构固有特性的影响机制和结构发生颤振的规律。采用非定常计算流体力学耦合结构传热的算法,获得结构的瞬态温度场;通过有限元计算得到瞬态温度场不同时刻下的热结构模态;利用结构模态叠加法建立结构动力学模型,结合Volterra级数建立的非定常气动力模型进行气动/结构耦合计算获得颤振动压。对马赫数为5、高度13 km的一小展弦比机翼进行了颤振分析。验证了该方法的可实现性。研究表明,随着温度梯度的增大结构固有频率减小,颤振动压最小值发生在结构主频率差值最小处。     

基于流固耦合方法的大展弦比机翼非线性颤振特性分析

针对大展弦比机翼的柔性大、变形大的特点,基于非定常涡格法求解机翼的非定常气动力,考虑了大展弦比机翼的几何非线性效应,提出了计算大展弦比机翼非线性颤振分析的新方法。以某平板机翼为例:计算了机翼静气动弹性变形、振动特性和颤振特性随着攻角增大的变化规律;比较了颤振结果线性解与非线性解的差别。相关的分析结果表明:大展弦比机翼颤振分析需同时考虑几何非线性效应和气动网格变形。     

猴子岩水库色玉滑坡涌浪灾害链CFD–DEM耦合数值模拟

滑坡涌浪是一种常见的地质灾害现象,由于滑坡体与水体之间存在复杂的流固耦合作用,使得传统的单一介质模型无法进行准确求解。为此,介绍一种基于计算流体力学方法（CFD）与离散单元法（DEM）的流固耦合模型CFD–DEM,采用计算流体力学方法（CFD）求解水体流动,采用离散单元法（DEM）模拟散粒体滑坡运动,充分利用不同计算模型的优势,对滑坡及涌浪演进过程进行数值模拟分析。首先,利用该耦合模型对Robbe–Saule开展的颗粒堆积体坍塌–涌浪试验进行了相同工况下的数值计算,从颗粒坍塌运动过程、涌浪高度演化过程等方面进行对比,结果表明模拟结果与试验结果吻合很好,验证了CFD–DEM流固耦合模型的有效性。然后,将该方法应用于四川省猴子岩水库色玉滑坡–涌浪灾害的演进过程分析,重现了该事件滑坡失稳运动、涌浪产生及传播、涌浪爬升、涌浪回流的全过程,计算结果显示：计算得到的电站进水口处涌浪高度与实测数据较为接近;色玉滑坡从失稳运动至静止堆积的持续时间约为20 s,颗粒平均速度最大达到16.12 m/s;滑坡引起的涌浪约在滑坡失稳10 s后传播到对岸,之后开始沿坡面向上爬升,最大爬升高度达到27.32 m。...     

船舶螺旋桨流固耦合稳态求解算法

为了研究复合材料变形特性对螺旋桨水动力性能的影响,运用计算流体力学(CFD)理论,结合雷诺时均纳维斯托克斯(PANS)方程和3种湍流模型,计算了螺旋桨的敞水性能并将其和试验结果进行了对比,两者吻合较好,验证了流体计算的准确性.推导了变形后的螺旋桨纵倾、侧斜、叶切面螺距以及拱度等几何参数,并将CFD软件计算和有限元求解耦合起来,即首先利用流体计算软件计算桨叶压力载荷,再将流体载荷通过插值技术传递给有限元,计算出螺旋桨的变形,然后再将变形后的桨叶几何输入流体计算软件,反复迭代,直至结果收敛的过程,实现了螺旋桨稳态性能求解的流固耦合算法.该算法可以计算螺旋桨在不同转速和进速下的应力分布和变形,能够求解螺旋桨变形后的压力分布和水动力性能,适用于复合材料螺旋桨的设计.     

几何大变形太阳能无人机非线性气动弹性稳定性研究

大柔性太阳能无人机在气动载荷的作用下产生较大的弯曲变形,机翼结构的刚度、质量分布等特性亦发生较大改变,线性理论无法满足这类飞机气动弹性稳定性分析的精度要求。基于Co-rotational(CR)理论,推导了结构变形后的切线刚度矩阵和质量矩阵,建立了大柔性机翼结构动力学模型;采用建立在局部气流坐标系下的片条非定常气动力模型,建立了考虑几何非线性效应的大柔性无人机气动弹性运动方程。引入准模态假设,采用P-k法研究了几何大变形对类"太阳神"布局太阳能无人机的气动弹性稳定性的影响。研究结果表明:随着弯曲变形的增加,非线性颤振速度可降低10%以上,非线性颤振频率可下降8%;合理的增加扭转刚度、前移弹性轴、前移剖面质心等,均可以有效改善几何大变形引起的不利影响。研究工作对大柔性飞机的气动弹性设计具有一定的参考意义。