基于多点约束的大展弦比机翼静气动弹性计算

基于雷诺平均N-S方程和多块结构网格技术计算大展弦比机翼的气动力.机翼采用梁模型结构,利用有限元方法计算结构变形.建立基于多点约束(MPC,Multi-Point Constrain)的气动、结构数据的双向传递方法.采用MPC数据插值方法快捷地实现自主研发的计算流体力学(CFD,Computational Fluid Dynamics)程序和NASTRAN计算结构力学(CSD,Computational Structural Dynamics)软件的耦合计算,发展出可以考虑结构非线性的静气动弹性CFD/CSD耦合计算方法,该方法比只考虑单方向变形的柔度法更精确.开展了某大展弦比机翼的静气动弹性数值计算分析,结果显示只考虑机翼纵向变形而忽略展向变形对计算结果会有一定影响,而该算例几何非线性影响较小.     

水平弯曲刚度对大展弦比机翼颤振的影响

高空长航时飞机普遍采用小后掠角大展弦比机翼,飞行时机翼具有较大变形．几何非线性加剧机翼结构水平弯曲与扭转的刚度耦合,水平弯曲模态在颤振分析中的作用已不可忽略,对机翼气动弹性特性影响显著．将结构在给定来流攻角的静平衡位置附近线化,以求解结构的固有振动特性,并应用考虑翼面变形的片条理论计算非定常气动力,用 p-k 法计算颤振速度．以金属梁式机翼为对象的颤振计算结果表明,水平一弯模态参与耦合的机翼颤振速度低于线性颤振速度．增大水平弯曲刚度有助于这类颤振速度的提高,而扭转频率的影响也要加以考虑．     

锥形後掠机翼的应力分析

本文参■计算锥形机翼基本应力的方法,■根据导师专家在无锥度后掠机翼中计算二次应力的方法,结合起来,试图对锥形后掠机翼进行应力分析。若有错误之处,敬希指正。 本文是在专家和王德荣教授二位导师指导下完成的,初稿又经王德荣教授审阅指正,使本文有许多改进,特在此向二位导师表示衷心的感谢。     

大展弦比柔性机翼的气动弹性分析

高空长航时飞机普遍具有大展弦比机翼,其气动弹性问题尤为突出.考虑机翼结构的几何非线性对结构刚度的影响和机翼平面变形对非定常气动力的影响,采用准模态的思想进行大变形下的振动工程分析,并利用片条理论计算变形机翼的非定常气动力,然后进行颤振计算.以某型飞机机翼为例的计算表明,随变形幅度增大,机翼的颤振速度有明显的下降.      

用三力矩方程计算小展弦比空心机翼

维廉曾用变形法计算小展弦比空心机翼,该法将节点挠度作为求知数,设节点较多,解联立方程系便较困难。本文目的有二:1.用三力矩方程计算由翼梁及翼肋传来的节点载荷,以简化节点载荷的计算;2.用初参数法来计算所得方程系,最后联立方程数目可减少一半以上,因而使计算大为简化。用一个具体机翼来说明计算步骤。 作者在写稿时,曾承苏联专家费阿发诺夫同志指教;脱稿后又承叶逢培和诸德超两同志指正,在此一并向他们表示衷心的感谢。     

充气式机翼的颤振特性分析

充气式机翼的结构刚度由内充气压决定,其颤振特性需要建立静、动力学耦合的分析方法.机翼结构刚度和固有振动特性需要在静力分析基础上计算,进一步计算非定常气动力,从而采用传统的颤振计算方法分析其颤振特性.针对某一充气式机翼采用膜单元建立了有限元模型.在不同内充压条件下,对充气机翼进行了静力分析得到其结构刚度;然后对机翼进行模态计算和颤振分析.研究表明:各阶模态的频率随内充气压的升高而升高;除典型的弯扭模态外,充气机翼的弦向弯曲模态频率较低;充气机翼的颤振形式除常规的弯扭模态耦合外,弦向弯曲模态同样会发生颤振;机翼的临界颤振速度随内充压的变化近似呈分段线性变化;临界颤振模态及耦合分支在一定气压范围内保持不变.      

基于Euler方程的三维机翼厚度与扭角优化设计

使用Euler方程作为流场解算器,结合Powell算法,讨论了在确定的机翼平面形状和翼型的条件下,以最大升阻比为目标的三维机翼截面翼型最大厚度与扭角的优化设计.设计中以机翼沿展向0.1,0.2,0.3,0.8,1.0倍半展长的翼剖面的厚度作为优化控制量,再选取翼梢剖面的扭角增量为第6个控制量,对机翼作了数值优化设计计算,得到了在亚音速时具有相对较大升阻比的机翼其厚度及扭角的优化分布.针对Lockheed-AFOSR Wing B的亚临界和超临界算例结果表明,厚度的非线性分布和负的扭角会改善机翼流场的流动状态,使机翼的升阻比得到提高,优化设计方法是可行的.      

考虑几何非线性的气动弹性模型缩比方法

随着飞机性能和需求的提高,大展弦比高柔性机翼逐渐成为新型飞机的主要结构形式。这类机翼具有高升阻比、大变形和重量轻等特性,几何非线性效应明显。然而机翼的大展弦比高柔性会带来更大的机翼变形,而机翼大变形则会引起相关的非线性气动弹性行为。为了评估这些非线性气动弹性行为并同时降低设计风险和成本,一般要使用缩比模型进行风洞试验以研究和确认真实飞机的气动弹性特性。基于此,首先使用了传统线性缩比方法来进行缩比,通过刚度质量耦合匹配模态响应法与刚度质量解耦匹配模态响应法这2种线性缩比方法,不断优化缩比结构的设计参数来满足目标缩比值。同时,提出一种动力学有限元模型的非线性静响应-模态协同优化方法,该方法是基于等效静态载荷法的几何非线性气动弹性模型缩比方法,通过2个不同的优化子程序分别匹配全尺寸飞机的非线性静响应和模态振型。结果表明,相比于传统线性缩比模型,考虑几何非线性的缩比模型能够更好地再现全尺寸飞机的非线性气动弹性行为。      

跨音压气机叶栅流场的附面层修正

本文发展了二维压气机叶栅亚音、跨音流中粘性——无粘相互作用程序。无粘计算在弦向大扰动、弦法向小扰动的假设条件下用跨音松弛法求解;粘性附面层用Nash方法计算。考虑了叶栅尾迹的影响。最后,连同通道区和尾迹区一起,求得了包括冲波影响的叶栅的总损失系数。计算结果表明,在跨音情况下,粘性对冲波强度、位置以及全流场的速度分布都有不可忽视的影响。     

基于变弯度后缘的机翼阵风响应减缓数值研究

针对带有变弯度后缘的机翼建立了阵风响应分析的数学模型，并开展了阵风响应减缓的仿真研究。采用计算流体力学（CFD）方法计算变弯度后缘给定动态偏转运动下的广义非定常气动力，并基于状态观测器法辨识CFD数据建立后缘动态偏转下的广义气动力模型，采用面元法计算模态运动、阵风引起的广义气动力，利用广义预测控制（GPC）方法设计阵风减缓控制律，在此基础上对变弯度后缘与传统铰链舵面动态气动力特性进行对比。仿真结果表明：基于变弯度后缘的GPC方法能够有效减缓由阵风引起的机翼翼尖加速度响应，翼尖加速度减缓效率为44.33%；相比传统铰链舵面，变弯度后缘偏转时弦向剖面上下表面压力分布更连续，相同舵偏下对机翼动态气动力影响更大，阵风响应减缓效率也更高，采用变弯度后缘进行阵风减缓具有更为广阔的应用前景。     

利用硬钣弯曲理论简化三角形机翼用变形法的解法

在应用变形法解空心三角形机翼时,通常把这机翼分解为翼梁、翼肋以及盒段,再分别利用这些构件的影响系数写出构件节点位移为外力的函数,而后解联立方程得出该构件外力为节点位移的函数,即所谓求逆矩阵。最后根据节点的平衡条件得出一组变形法的正则方程,如文献所述。由於在计算过程中需要多次求逆矩阵,所以相当费事,并且这种方法无法致虑到蒙皮的波桑系数,对於较厚蒙皮的机翼就不适宜。本文根据硬钣弯曲理论,用有限差分方程代替偏微分方程,从而直接得出机翼节点外载荷为节点位移的函数。这样就避免了多次求逆矩阵的计算,并且,也把蒙皮的波桑系数致虑进去,因而适用厚蒙皮的情况。本文给出了各种在计算小展弦比机翼时可能遇到的边界条件。     

基于升力线理论的机翼几何扭转设计方法

针对飞机设计对于减小诱导阻力的需求,结合升力线理论中三维机翼诱导阻力的相关研究,提出机翼的几何扭转设计方法.此设计方法以总升力、机翼平面形状、翼型气动特性为输入参数,以机翼的几何扭转为输出.在低速情况下,利用此方法得到的机翼具有椭圆形升力分布,从而具有最小的诱导阻力.采用远场法计算诱导阻力,应用参数化方法和CATIA实现机翼的自动生成,应用网格生成工具ICEM的脚本语言实现气动网格自动生成,应用MATLAB的符号运算功能编制程序.该设计方法适用于低速有粘流中平直翼的几何扭转设计,算例的CFD(Computational Fluid Dynamics)计算结果表明:利用该方法设计的机翼展向升力分布为椭圆,并能有效降低机翼的诱导阻力.     

基于升力线理论的机翼几何扭转设计方法

针对飞机设计对于减小诱导阻力的需求,结合升力线理论中三维机翼诱导阻力的相关研究,提出机翼的几何扭转设计方法.此设计方法以总升力、机翼平面形状、翼型气动特性为输入参数,以机翼的几何扭转为输出.在低速情况下,利用此方法得到的机翼具有椭圆形升力分布,从而具有最小的诱导阻力.采用远场法计算诱导阻力,应用参数化方法和CAT-IA实现机翼的自动生成,应用网格生成工具ICEM的脚本语言实现气动网格自动生成,应用MATLAB的符号运算功能编制程序.该设计方法适用于低速有粘流中平直翼的几何扭转设计,算例的CFD(Computational Fluid Dynamics)计算结果表明:利用该方法设计的机翼展向升力分布为椭圆,并能有效降低机翼的诱导阻力.     

双尾撑布局弹性飞机配平诱导阻力分析与优化

基于能量方法推导了弹性飞机静气弹配平方程,提出了用定常涡格法结合Trefftz平面理论与阻力计算经验公式来求解弹性飞机定直平飞配平状态诱导阻力与配平阻力的分析思路,建立了通过机翼几何扭转角优化配置减小全机诱导阻力的分析方法.以某双尾撑无人机为例进行研究,分析发现尾撑纵向弯曲刚度对该布局飞机的全机诱导阻力、升降舵配平舵偏型阻及全机阻力等参数有重要影响.优化分析结果表明:通过机翼扭转角的合理优化,可使全机展向升力分布形式更接近于椭圆分布,进而有效减小全机诱导阻力,且对其他阻力分量也有一定的减缓效果.     

无限多自由度连续梁的自由频率方程

本文系利用三转角方程组的行列式解法,渐近解法和力矩分配法,系统地分析力矩分配和剖面旋转等动作,进而说明剖面抗弯刚度和它的修正系数的意义。根据三转角方程组的特点,照高斯消去法,进行系数矩阵行的变换,把它变为三角形矩阵,导出自由频率的一般方程。并可据此用计算机来计算自由频率。     

柔性飞机大变形状态飞行载荷分析

飞行载荷分析中综合考虑了结构大变形的几何非线性效应以及曲面气动力效应.飞机结构由相互连接的可以表征任意变形的几何非线性欧拉梁表示,升力面由展向以及弦向分布的涡环网格表示.计算时通过曲面网格动态跟踪结构的变形不断修改气动力影响系数矩阵,反复迭代求解气动力和结构变形直至结构变形收敛.给出了风洞试验机翼模型以及无约束定常平飞模型算例.计算结果表明:在小变形阶段该方法与线性计算方法的结果基本一致,非线性效应不明显;大变形阶段由于曲面气动力的非线性效应,计算结果与线性方法的有显著差别.分析表明该方法在大变形阶段的计算结果比线性结果更为合理,数值计算时间短,适合于工程快速分析.     

复合材料机翼翼梁与蒙皮间应力分析

本文采用三维有限元模式，用“ＡＢＡＱＵＳ”分析程序，进行了复合材料机翼翼梁和蒙皮结构的应力分析，对所获得的缘条和蒙皮之间六个应力分量的分布曲线，与采用二维有限元模式的结果做了对比分析，得出了一些有益的结论。     

鱼骨柔性翼段线性/非线性静气动弹性对比分析

鱼骨柔性翼作为一种性能优越的主动变弯度机翼机构形式，具有弦向抗弯刚度低、翼型厚度方向刚度高的特点，其在进行大弯度主动变形时结构存在较强的几何非线性且气动弹性效应显著。针对传统线性静气动弹性分析方法并不适用于鱼骨柔性翼段的气动弹性分析问题，以Bristol大学的公开鱼骨柔性翼段模型为研究对象，采用曲面涡格法（VLM）和非线性有限元耦合的非线性静气动弹性方法，以及传统气动弹性分析中常用的平面涡格法和线性有限元耦合的线性静气动弹性方法，分别对鱼骨柔性翼段进行大变形下的静气动弹性分析，并进行结果对比。对比验证了所用曲面涡格法与XFOIL软件气动计算结果。算例结果表明:鱼骨柔性翼段大变形下气动弹性效应显著，相比传统线性静气动弹性分析方法，非线性静气动弹性分析方法得到的鱼骨柔性翼段在大变形状态下升力系数最多减少8.28%，力矩系数最多减少6.86%，且能准确快速得到真实变形结果，更具有实际工程应用价值。      

快速计算虚拟物体之间精确接触位置的算法

在触觉、力觉交互的虚拟现实系统中,快速计算虚拟物体之间精确接触位置是触觉、力觉交互碰撞检测系统的前提.在这样的系统中虚拟物体大多是多面体模型且是三角形模型,采用层次包围盒法判断物体是否相交,如果相交,采用本文提出的二级检测算法快速计算出物体的接触位置.该算法首先判断那些相交三角形是异面还是同面,如果是异面,求出一个三角形T₂与另外一个三角形T₁所在平面的交线,然后再判断该交线与T₁之间的位置关系,进而求出有效的相交线段;如果是同面,需要判断两个三角形之间相互位置关系,然后有针对性地求出相关线段之间的交点,从而求出三角形的相交区域.最后用实例验证了该算法的有效性和鲁棒性.