飞机主承力结构细节分析

工程上往往会遇到如飞机、轮船、汽车等大型结构,对这些大型结构进行有限元分析有时须同时考虑整体以及局部。对于大型结构采用较密的网格,这样会耗费大量时间、资源等,有时会导致计算不收敛;对于局部结构须划分较细的网格才能得到局部构件的细节应力。本文采用子结构法（超单元法）对机翼整体模型中的主承力构件梁进行网格细化,其中梁筋条与肋以及梁缘条与机翼蒙皮采用钉元（弹簧元）模拟多钉连接。在危险载荷工况下对梁进行强度校核和钉载分析,结果表明子结构法能很好地反映大型结构的整体以及细节应力分布,钉元法能很好地模拟机翼多钉连接。     

小型昆虫“合拢-打开”运动机制的增升机理数值研究

通过动态混合网格技术和基于虚拟压缩方法的不可压缩流非定常计算方法,对"合拢-打开"(clap-fling)机制下小型昆虫翼的周期性运动进行了数值模拟,并与单个翼拍动时的升力系数、流场结构进行了对比分析,同时对比研究了雷诺数对力学特性、流场结构的影响,探讨了"clap-fling"机制能够增强升力的流动机理.对比发现,"clap-fling"机制增升的关键在于合拢-打开过程中的三个升力峰值以及之后平动阶段的"不失速机制";且当雷诺数较小时,"clap-fling"机制增升效果最为明显.     

二维间断有限元混合网格混合算法研究

针对二维流场进行结构／非结构混合网格划分,在靠近物面处生成结构化网格,在结构化网格以外和远场之间以非结构化网格填充。使用间断有限元方法（DGM）离散Euler方程。在混合网格混合算法中,结构化网格区域使用结构DGM,非结构化网格区域使用非结构DGM,保证结构／非结构对接面处的通量守恒,并加入moment限制器消除激波前后伪震荡。对NACA0012和RAE2822翼型跨音速无粘流动的数值模拟结果说明了混合算法的有效性以及在捕捉间断方面的能力。     

三维非结构贴体直角网格的N-S方程数值模拟

基于八叉树数据结构,采用直接分层光顺的方法,解决了非结构直角网格在物面不正交贴体的难点,发展了任意复杂外形的三维非结构贴体直角网格,并且可以根据流场特征进行自适应处理。在此基础上采用最小二乘法对N-S方程进行了MUSCL有限体积法的非结构网格的粘性求解,提出了带约束边界条件的最小二乘法的处理办法。计算表明,三维非结构贴体直角网格的粘性求解,可以适应各种高超声速的复杂外形,计算结果与试验吻合较好。     

Euler方程的自适应多重非结构网格计算

介绍一种基于四面体网格的多重网格技术,用于求解三维Euler方程定常解,具有十分明显的加速收敛效果。求解Euler方程的基本算法采用有限体积中心差分格点格式。多重网格技术采用一组几何上相互独立的粗细不同的网格。流动变量、余量和修正量在相邻粗细网格间的相互传递通过线性插值实现。采用自适应技术对网格进行加密,可以获得较优的网格,并作为一套细网格,进行多重网格计算。数值实验表明本文自适应多重网格法十分有效。     

立式捏合机混合釜内推进剂药浆混合过程数值仿真

使用POLYFLOW软件对固体推进剂药浆在捏合机混合釜内的流动过程进行了数值仿真,采用混合釜壁面相对桨叶反方向旋转的方法实现公转,桨叶只需要绕自身轴线白转即可实现行星运动;运用网格叠加技术并结合相应的网格加密技术,完成了复杂边界条件下的网格划分。实现了混合过程巾药浆流动的可视化,量化了捏合机混合过程中釜内各点的速度、剪...     

轴向入射激波反射聚焦的实验和数值模拟

为研究两级脉冲爆震发动机凹腔入口宽度对第2级凹腔中激波聚焦过程的影响,设计轴向入射激波聚焦的实验系统,对不同入口宽度下激波反射聚焦的过程进行了实验和数值模拟。采用纹影系统拍摄凹腔中的流场结构,并测量了凹腔顶点处的动态压力。采用保持强稳定性（SSP）的Runge Kutta格式、weighed essential non oscillation（WENO)格式和块结构网格自适应加密（SAMR）算法对激波反射、衍射和聚焦的过程进行了数值模拟,数值模拟结果和实验吻合较好。通过分析比较相同马赫数激波通过不同宽度凹腔的反射聚焦过程,发现较小宽度凹腔入口下激波反射和聚焦较强,顶点的压力峰值较高,激波聚焦诱导的射流较强。      

一种基于径向基函数和峰值选择法的高效网格变形技术

基于径向基函数的网格变形方法因其具有诸多优点,而被广泛应用于气动外形优化设计等领域。对于大规模网格或复杂构型,该方法所需计算量是难以承受的。为了提高网格变形效率,可以通过减少建立插值模型所需支撑点数目来实现。为此,提出一种高效的选点算法--峰值选择法。该算法在选点过程的每个迭代步中对边界节点处的误差进行分析,从物面节点中选取多个峰值点来更新支撑点集,减少迭代步数,提高选点效率。在该算法的基础上,实现了网格的高效变形。三段翼型的网格变形算例证明：该方法可以在保证网格质量的同时实现复杂网格的变形。以DLR-F6复杂模型（约1 000万网格）的刚性运动和弹性大变形为算例对该方法的变形效率和变形后网格质量做了进一步评估：当相对误差设置为5.0×10^-7时,在保证变形后网格质量的前提下,该方法变形效率最快比传统贪婪算法提高了13倍,其中在选点效率方面最快提高了31倍。     

子母弹不同舱段分离流场特性及运动特性研究

为研究子弹药在不同装配舱段下抛撒分离过程的干扰气动特性,以2层弹舱轴向排布子母弹系统为模型,基于有限体积法,结合非结构动网格技术,通过耦合求解任意拉格朗日欧拉(ALE)描述下的三维非定常可压缩N-S流体控制方程及六自由度刚体运动方程,对时序抛撒方式下子弹在不同舱段分离的三维非定常流场进行了数值模拟。得到了不同舱段下子母弹分离流场的干扰特性及子弹药气动参数变化曲线,分析了子弹分离过程的运动特性,揭示了子弹与母弹激波在不同分离阶段的相互作用过程。为进一步研究子母弹分离干扰流场机理提供了依据。     

考虑动力影响的民用飞机静气动弹性分析方法

采用了一种基于多块网格的N-S方程和结构柔度影响系数法,考虑气动、结构非线性的基于RBF插值和RBFDelaunay动网格变形技术的静气动弹性分析方法对喷流对弹性机翼的气动力影响进行了研究。利用DLR F6翼身组合体构型对静气动弹性方法进行验证,保证了计算的可信性。采用该方法对比分析了某民用飞机无喷流/有喷流构型的静气动弹性特性,表明发动机喷流会给机翼带来一个正的扭转效应,抵消一部分机翼后掠效应的影响,使机翼前后缘挠度均会有所增大,弹性变形引起的多数剖面的附加扭转角有所减小。研究表明:喷流影响会使刚性机翼表面的压力分布发生变化,升力系数有所损失;考虑喷流的机翼静气动弹性变形是一个耦合效应,发动机喷流区主要受喷流影响,外翼段主要受弹性变形影响。数值模拟结果表明:无喷流影响时机翼的弹性变形使升力系数下降约16%,升阻比下降8.4%,考虑喷流影响时,升力系数下降达到18%,升阻比下降36%。因此,对于大展弦比机翼,考虑喷流影响的静气动弹性分析十分必要。