基于RBF插值技术的CFD/CSD非线性耦合分析方法研究

基于非结构混合网格的N-S方程求解器和结构柔度影响系数法,发展了一种考虑气动、结构非线性的基于RBF插值技术CFD/CSD耦合分析方法,适用于解决现代大展弦比飞机的非线性静气动弹性问题。该方法采用时间相关法(即求解非定常方程组,用长时间的渐近解趋于定常状态)求解静气弹分析时的定常流动。考虑大展弦比飞机结构变形问题为大变形小应力问题,在利用柔度系数法求解结构方程时,假设每次求解结构方程时应力与应变为线性关系,整体静气弹分析过程为非线性关系,因此每次求解结构方程时要更新柔度影响系数矩阵。在非定常N-S方程每求解一个时间步耦合一次结构有限元分析,由于结构有限元分析的时间相对于气动分析时间是很短的,所以这种方法实际上近似使用了一次求解非定常气动力的时间完成了整个静气动弹性分析的过程。对于气动网格与结构有限元网格不一致性,本文采用径向基函数(RBF)插值方法中的TPS方法进行结构弹性变形和气动载荷插值,采用虚功原理完成气动载荷数据交换。为了节省气弹分析时间,采用动网格方法对气动网格进行更新,本文基于RBF插值方法发展一种适用于混合网格(四面体、三棱柱、金字塔和六面体)变形的动网格方法,可以保证附面层网格的质量与分布从而准确模拟其流动。利用该方法对M6机翼、DLR-F6翼身组合体和某大型客机机翼进行了静气动弹性特性分析,结果验证了本文开发的非线性CFD/CSD耦合分析方法的可行性、精确性和高效性。     

基于RBF插值技术的CFD/CSD非线性耦合分析方法研究

基于非结构混合网格的N-S方程求解器和结构柔度影响系数法,发展了一种考虑气动、结构非线性的基于RBF插值技术CFD/CSD耦合分析方法,适用于解决现代大展弦比飞机的非线性静气动弹性问题。该方法采用时间相关法（即求解非定常方程组,用长时间的渐近解趋于定常状态）求解静气弹分析时的定常流动。考虑大展弦比飞机结构变形问题为大变形小应力问题,在利用柔度系数法求解结构方程时,假设每次求解结构方程时应力与应变为线性关系,整体静气弹分析过程为非线性关系,因此每次求解结构方程时要更新柔度影响系数矩阵。在非定常N-S方程每求解一个时间步耦合一次结构有限元分析,由于结构有限元分析的时间相对于气动分析时间是很短的,所以这种方法实际上近似使用了一次求解非定常气动力的时间完成了整个静气动弹性分析的过程。对于气动网格与结构有限元网格不一致性,本文采用径向基函数（RBF）插值方法中的TPS方法进行结构弹性变形和气动载荷插值,采用虚功原理完成气动载荷数据交换。为了节省气弹分析时间,采用动网格方法对气动网格进行更新,本文基于RBF插值方法发展一种适用于混合网格（四面体、三棱柱、金字塔和六面体）变形的动网格方法,可以保证附面层网格的质量与分布从而准确模拟其流动。利用该方法对M6机翼、DLR-F6翼身组合体和某大型客机机翼进行了静气动弹性特性分析,结果验证了本文开发的非线性CFD/CSD耦合分析方法的可行性、精确性和高效性。     

考虑几何非线性的复合材料机翼气动弹性分析

本文研究结构几何非线性与气动力非平面效应对大展弦比复合材料机翼的气动弹性行为的影响．将非线性有限元法与曲面涡格法结合，计算机翼静气动弹性变形；通过曲面偶极子格网法结合静气动弹性平衡位置处的结构切线刚度，建立气动弹性方程并求解得到机翼颤振速度．针对板模型机翼，分析了迎角对机翼几何非线性气动弹性特性的影响．结果表明：本文复合材料板模型机翼的颤振形式不受水平弯曲模态影响，属于经典弯扭颤振；在几何非线性的影响下，机翼扭转频率随结构变形增大而明显减小，颤振速度随迎角增大而减小．     

大展弦比无人机的静气弹问题计算及分析

普遍采用大展弦比机翼的无人机,其气动弹件的问题显得尤为突出.本文主要研究基于非结构弹性网格体系的欧拉方程CFD计算及其在大展弦比无人机静气动弹性问题中的应用.针对三维非结构运动网格技术进行了研究和开发,在此基础上,利用计算流体力学的方法,发展了一套具有一定通用性的、适用于非结构网格的Euler方程流场求解器,并综合上面的技术,进一步通过耦合结构力学方程,对大展弦比无人机的静气动弹性问题进行了,计算和分析.     

铺层参考方向选取对大展弦比机翼静气弹特性影响

以大展弦比机翼为研究对象，利用流固耦合方法对复合材料机翼铺层参考方向进行了数值模拟研究，分析了铺层参考方向轴偏角的改变对大展弦比机翼静气动弹性的影响．研究表明：铺层参考方向轴偏角的改变会对机翼气动弹性产生显著的影响．机翼的总体变形与扭转变形随着参考方向轴偏角的改变呈现周期分布；沿着机翼各个方向的挠度也会因为参考方向轴偏角的改变而产生不同的响应．     

考虑流固耦合效应的某飞行器力学性能分析

考虑流固耦合效应,研究了飞行器在服役状态下的力学特性。本文采用了一种基于CFD/CSD(Computational Fluid Dynamics/Computational Structure Dynamics)耦合的高精度通用静气弹求解器,该求解器可以同时分析飞行器结构的表面气动性能和结构的力学性能,对某一翼身组合体模型进行了静气弹分析,得出弹性飞行器和刚性飞行器的气动压力变化以及结构的应力分布。建模时流体和结构模型都使用翼身组合模型,且都考虑机身的影响。结果表明:发生静气弹变形后,飞行器的升力系数、阻力系数、力矩系数都相应下降;机翼翼根附近会产生较大的应力分布;当飞行器尤其是机翼表面受到弯矩、扭矩以及气动压力的共同作用时,翼根附近会出现复杂的应力分布状态。这说明在静气弹性分析中,机身的建模也是很有必要的。     

基于传递函数法的带多个外挂机翼颤振分析

基于带外挂机翼结构和气动特点，使用带有半解析半数值特性的传递函数方法进行处理．首先，通过机翼的运动微分方程、二元机翼非定常气动力Therdorson 模型、结合外挂挂载处的内力平衡与位移状态条件，得到了三维的带多个外挂大展弦比机翼的颤振微分方程．进而，使用传递函数方法，先将颤振计算方程整理成为状态空间方程形式，结合求解复特征值的方法，完成了带多个外挂大展弦比机翼的动气动弹性稳定性分析．对比已有文献的计算结果以及通过有限元方法进行的仿真结果，证实了文章所提计算方法的准确性和高效性，结合传递函数方法的优势，进一步将文章方法拓展到机翼的固有频率和发散速度求解．文章结尾，分析了外挂数量、外挂质量、转动惯量及位置分布等变量对带多个外挂大展弦比机翼的动气动弹性稳定性的影响．     

基于POD和代理模型的静气动弹性分析方法

静气动弹性问题考虑弹性结构与定常气动力间的相互耦合作用,对飞行器的性能和安全具有显著的影响.在现代飞行器设计阶段,计算流体力学(CFD)/计算结构力学(CSD)直接耦合方法是精确考察静气动弹性影响的重要手段.然而,基于CFD技术的气动力仿真手段在耦合过程中计算量大且耗时长,难以满足设计阶段的需求.因此,为了兼顾计算精度与效率,文章采用本征正交分解(POD)和Kriging代理模型相结合的模型降阶方法,替代CFD求解过程并耦合有限元分析(FEA)方法,建立了高效、准确的静气动弹性分析框架.相较于传统的以模态法为主的静气动弹性分析方法,该方法能够解决更为复杂的静气动弹性问题以及提供静气动弹性变形过程中的气动分布载荷.针对典型三维跨声速HIRENASD机翼模型开展的马赫数、迎角变化的算例验证表明:由建立的静气动弹性分析方法与CFD/CSD直接耦合方法计算得到机翼翼梢处的静变形量间的相对误差在5%以内;同时该方法预测静平衡位置处的气动分布载荷的误差在5%以内,静气动弹性分析的计算效率至少提升了6倍.     

静气动弹性计算方法研究

对基于结构网格的Euler方程及N-S方程求解器和基于非结构网格的Euler方程求解器,采用结构模态分析方法和柔度矩阵方法,对无人机大展弦比机翼在Ma=0.6,α=2?, 飞行高度20km的巡航状态下的静气动弹性特性进行了数值模拟. 验证了两种求解器对静气动弹性模拟的准确性. 同时,对模态分析方法和柔度矩阵方法进行了对比研究,发现柔度矩阵方法更适用于静气动弹性数值模拟. 另外,对应用物面法向偏转方法替代网格变形技术模拟静气动弹性进行了研究,计算表明物面法向偏转方法可以大大提高静气动弹性计算效率和克服机翼结构变形过大时动网格技术无法处理的不足.      

西南某大型水电站拱肩槽边坡开挖的变形响应研究

对基于结构网格的Euler方程及N-S方程求解器和基于非结构网格的Euler方程求解 器,采用结构模态分析方法和柔度矩阵方法,对无人机大展弦比机翼在Ma=0.6, α=2?, 飞行高度20km的巡航状态下的静气动弹性特性进行了数值模 拟. 验证了两种求解器对静气动弹性模拟的准确性. 同时,对模态分析方法和柔度 矩阵方法进行了对比研究,发现柔度矩阵方法更适用于静气动弹性数值模拟. 另外, 对应用物面法向偏转方法替代网格变形技术模拟静气动弹性进行了研究,计算表明 物面法向偏转方法可以大大提高静气动弹性计算效 率和克服机翼结构变形过大时动网格技术无法处理的不足.