结构非线性颤振半主动抑制

本文利用数字仿真技术，对结构非线性颤振半主动抑制－颤振驯化方案进行了探讨，仿真结果表明，对于一定的结构非线性类型和参数，利用非线性颤振的极限环特性，可使系统的颤振响应被抑制在幅值很小的稳定域内，从而达到减缓颤振的目的。     

不可压缩流中机翼颤振稳定性研究

研究两个自由度的机翼在不可压缩流作用下颤振的分支问题.运用罗司-霍维茨判据确定系统的分叉点,应用中心流形理论将四维系统降为二维系统,用直接求周期解方法对分叉点的真假中心及稳定性问题进行了分析,并研究了系统的极限环颤振.结果表明,本文研究的分叉点不是中心,而是稳定或不稳定焦点.在两个分叉点处,系统发生了超临界和亚临界Hopf分叉,产生稳定或不稳定极限环.     

大型飞机操纵面极限环颤振特性研究

操纵面间隙会引起结构非线性,容易诱发极限环颤振。根据设计需要,考虑操纵面中心间隙的影响,采用最小状态拟合技术对频域非定常气动力进行有理函数拟合,采用分段函数描述间隙引起的非线性刚度,研究操纵面在间隙作用下的极限环颤振响应的行为特点。结果表明,由于中心间隙的影响,系统会在低于线性颤振速度时就产生极限环振荡,同时,振荡幅值随飞行速度或中心间隙的增大而增大。     

大型飞机操纵面极限环颤振特性研究1)

操纵面间隙会引起结构非线性,容易诱发极限环颤振。根据设计需要,考虑操纵面中心间隙的影响,采用最小状态拟合技术对频域非定常 气动力进行有理函数拟合,采用分段函数描述间隙引起的非线性刚度,研究操纵面在间隙作用下的极限环颤振响应的行为特点。结果表明,由于 中心间隙的影响,系统会在低于线性颤振速度时就产生极限环振荡,同时,振荡幅值随飞行速度或中心间隙的增大而增大。     

LIMIT CYCLE FLUTTER CHARACTERISTICS RELATED TO CONTROL SURFACE OF LARGE AIRCRAFT1)

操纵面间隙会引起结构非线性,容易诱发极限环颤振。根据设计需要,考虑操纵面中心间隙的影响,采用最小状态拟合技术对频域非定常 气动力进行有理函数拟合,采用分段函数描述间隙引起的非线性刚度,研究操纵面在间隙作用下的极限环颤振响应的行为特点。结果表明,由于 中心间隙的影响,系统会在低于线性颤振速度时就产生极限环振荡,同时,振荡幅值随飞行速度或中心间隙的增大而增大。     

轴向流中悬臂柔性板流固耦合动力特性研究

论文从理论分析和实验研究两方面对轴向流中悬臂柔性板流固耦合动力特性作了研究.基于板的轴向不可延伸假设,采用Hamilton变分原理,推导了轴向流中悬臂板非线性运动偏微分方程,通过面元法计算悬臂板上下面压力差.采用伽辽金法离散该非线性运动偏微分方程,通过无量纲化物理参数,研究了系统无量纲颤振临界速度,给定流速下悬臂板末端响应时程曲线及板的振动响应,并与风洞模型实验结果进行了对比.研究结果表明,系统发生Hopf分叉,流体流速大于颤振临界速度时,悬臂板发生稳定的极限环颤振,悬臂板振幅、频率随流速的增加而增大.     

资金项目:国家自然科学基金资助项目

采用能计及非线性结构刚度的颤振方程为控制方程,和非定常N-S方程耦合求解,运用龙格-库塔方法在时域内求解结构响应的时间历程,从而确定颤振临界条件.计算了带结构刚度非线性的跨音速颤振特性.计算结果表明,结构刚度非线性对颤振特性有明显的影响.由于同时具有结构和气动力非线性,导致了具有复杂振荡极限环的特性.     

带外挂二元翼极限环颤振的高次线化分析

本文针对在外挂上带有初偏间隙型非线性刚度的二元翼带外挂系统的极限环颤振,进行KBM法二次渐近解等效线化分析.结果表明,在某些情况下,用二次渐近解等效线化分析极限环颤振,其结果与数值积分结果基本一致,而一次渐近解等效线化分析则得不出满意的结果.这就是说,二次解等效线化比一次解等效线化更为可靠.     

超声速飞行器壁板非线性颤振响应分析的时域法与频域法对比研究

为考查假设模态法应用于壁板非线性颤振响应时域分析的可行性,分别采用时域方法和频域方法研究了超声速飞行器壁板的非线性颤振响应,并从壁板的颤振幅值、颤振频率和颤振型态三个方面对时域和频域分析结果的一致性作了较详细的比较.首先,基于von Karman应变-位移关系和Mindlin板理论建立考虑几何非线性的壁板力学模型,应用一阶活塞理论分析壁板上单面承受的超声速准定常气动力,基于虚功原理和有限单元法推导壁板的运动微分方程.然后,用壁板的线性固有模态作为假设模态,减缩系统的自由度而得到降阶模型.采用四阶龙格-库塔法对降阶模型作时域数值积分,得到壁板的非线性颤振响应.另一方面,假设壁板的极限环颤振为简谐振荡,可对壁板的非线性刚度作等效线性化处理,进而在频域中直接在有限元(未降阶)模型的基础上分析壁板的颤振幅值、颤振频率和颤振型态.数值分析表明,当极限环颤振为简谐振荡时,时域方法和频域方法的计算结果符合一致.本文最后讨论了时域法和频域法应用在壁板非线性颤振分析中各自的优点和局限性.     

热环境下壁板非线性颤振分析

基于一阶活塞气动力理论,采用Von Karman大变形应变-位移关系建立了无限展长壁板热环境下颤振方程,采用伽辽金方法对方程进行离散处理.取温度为分叉参数,研究壁板颤振时的分叉及混沌等复杂动力学特性.结果表明:温度载荷降低了系统的颤振临界动压,改变了颤振特性.在整个分岔参数范围内,系统呈现出较为复杂的变化,包括衰减振动、极限环振动、拟周期振动和混沌型振动.当考虑材料热效应时,系统的颤振动压将进一步降低,其响应也表现出更为丰富的非线性动态力学行为.     

Parametric studies on relationships between flutter derivatives of slender bridge （Ⅰ）

Relationships between flutter derivatives of slender bridge are investigated based on our previously proposed semi-analytical flutter derivatives of flexible structure. The intrinsic relations are validated with test data of flutter derivatives of two bridges. Changes in flutter derivatives with the aerodynamic center, rotation speed, and angle variation are also studied by using a parametric method. The results show correctness of the proposed expressions of flutter derivatives given by authors in Ref. [1], and indicate that certain relations exist between these derivatives. It is also shown that semi-analytical flutter derivatives are applicable to bridges with a streamlined cross-section.     

联接刚度对机翼／外挂系统颤振边界特性的影响

本文对带外挂物二元机翼的颤振特性进行了理论和实验研究,主要分析机翼与外挂物之间俯仰联接刚度对颤振边界的影响。在大量算例的分析基础上,给出了颤振频率及颤振边界变化的一般规律及对颤振边界类型的一般判别方法,颤振模型的风洞实验验证了理论分析结果的正确性。     

高速小展弦比机翼颤振的微分求积法分析

引入微分求积法,分析高速小展弦比机翼的气动弹性问题。将小展弦比机翼等效为悬臂板,基于一阶活塞气动力理论建立机翼颤振偏微分方程,采用微分求积法将偏微分方程转化为常微分方程,根据频率重合理论对颤振问题进行求解。分析了机翼的固有频率及颤振速度,并与有限元软件计算结果进行比较,误差在2％以内,很好的验证了微分求积法求解小展弦比机翼颤振问题的有效性。分析了机翼面积、展弦比及厚度对颤振速度的影响,结果表明,小展弦比机翼的颤振速度受结构尺寸的影响较大,颤振速度随面积和展弦比的增大而减小,随机翼厚度的增大而增大。     

On bimodal flutter behavior of a flexible airfoil

The dynamic aeroelastic behavior of an elastically supported airfoil is studied in order to investigate the possibilities of increasing critical flutter speed by exploiting its chord-wise flexibility. The flexible airfoil concept is implemented using a rigid airfoil-shaped leading edge, and a flexible thin laminated composite plate conformally attached to its trailing edge. The flutter behavior is studied in terms of the number of laminate plies used in the composite plate for a given aeroelastic system configuration. The flutter behavior is predicted by using an eigenfunction expansion approach which is also used to design a laminated plate in order to attain superior flutter characteristics. Such an airfoil is characterized by two types of flutter responses, the classical airfoil flutter and the plate flutter. Analysis shows that a significant increase in the critical flutter speed can be achieved with high plunge and low pitch stiffness in the region where the aeroelastic system exhibits a bimodal flutter behavior, e.g., where the airfoil flutter and the plate flutter occur simultaneously. The predicted flutter behavior of a flexible airfoil is experimentally verified by conducting a series of systematic aeroelastic system configurations wind tunnel flutter campaigns. The experimental investigations provide, for each type of flutter, a measured flutter response, including the one with indicated bimodal behavior.     

Coupled fluid structure analysis for wing 445.6 flutter using a fast dynamic mesh technology

The flow around wing 445.6 was modelled using Navier–Stokes equations and S-A model. The wing vibration and flow mesh deformation were computed using a fast dynamic mesh technology proposed by our own group. Wing 445.6 flutter was analysed through a strong coupling between the wing vibration and flow. The reduced flutter velocity was predicted and results are in good agreement with the experimental data. It is found that the subsonic flutter is mainly induced by the flow separation and the transonic and supersonic flutter are mainly caused by the oscillating shock wave and its induced flow separation. The positive aerodynamic work increases due to the oscillating shock wave when the subsonic flow becomes transonic reducing the flutter velocity. While the positive aerodynamic work induced by the oscillating shock wave decreases when the transonic flow becomes supersonic increasing the flutter velocity. That is why the transonic dip exists.     

两种湍流模型时域颤振计算方法研究

采用时域计算分析方法进行了机翼跨音速颤振特性研究。在结构运动网格的基础上,采用格点格式有限体积方法进行空间离散和双时间全隐式方法进行时间推进求解雷诺平均N-S方程。针对流动粘性分别应用了SST湍流模型和SSG雷诺应力模型,通过对跨音速标模算例AGARD445.6机翼的计算结果与实验值的对比分析,其中应用SST湍流模型得到的颤振速度与实验值最为接近,特别是在跨音速段平均相对误差在3%以内;并且计算结果整体上反映了跨音速颤振"凹坑"物理特性,验证了方法的有效性。     

Parametric studies on relationships between flutter derivatives of slender bridge （Ⅱ）

Based on curve fitting of coefficients of three component forces of the Messina Straits Bridge, and the previously proposed semi-analytical expressions of flutter derivatives of flexible structure, the change of flutter derivatives of slender bridge cross-section with respect to its aerodynamic center, rotational speed and angle variation is studied using a parametric method. The calculated results are compared with the measured ones, and expressions of flutter derivatives of the Messina Straits Bridge are derived. The intrinsic relationships existing in flutter derivatives are validated again. It is shown that the influence of the rotational speed on flutter derivatives is not negligible. Therefore, it provides an additional semi-analytical approach for analyzing flutter derivatives of the bridge with streamlined cross-section to get its aerodynamic information.     

适用于几何非线性气动弹性分析的结构动力学降阶模型方法研究

几何非线性是壁板颤振和大展弦比机翼气动弹性等问题的一个主要特征,在进行数值仿真分析时往往需要采用商业非线性有限元求解器,存在计算量大和耦合迭代策略不易控制等问题。本文发展了一种适用于几何非线性的结构动力学降阶模型(CSD-ROM),利用广义坐标的非线性多项式表征非线性内力,采用参数识别方法获取多项式系数,并通过增加额外的线性模态来改善模型预测精度。基于此方法,分别针对壁板颤振、切尖三角翼的CFD/CSD-ROM非线性颤振问题开展了时域响应分析。计算结果表明,通过CSD-ROM计算出的壁板颤振速度为590 m/s,颤振频率为174 Hz,与有限元结果误差分别为0.8%和1.7%。马赫数0.879时切尖三角翼的颤振动压预测结果为2.25 psi,与非线性有限元相比的误差为3.8%。本文采用的非线性和线性模态基底组合方法,在保证计算精度的基础上可有效降低训练样本数量,一定程度上可替代非线性有限元开展气动弹性分析。     

Parametric studies on relationships between flutter derivatives of slender bridge (Ⅱ)

Based on curve fitting of coefficients of three component forces of the Messina Straits Bridge, and the previously proposed semi-analytical expressions of flutter deriva-tives of flexible structure, the change of flutter derivatives of slender bridge cross-section with respect to its aerodynamic center, rotational speed and angle variation is studied using a parametric method. The calculated results are compared with the measured ones, and expressions of flutter derivatives of the Messina Straits Bridge are derived. The in-trinsic relationships existing in flutter derivatives are validated again. It is shown that the influence of the rotational speed on flutter derivatives is not negligible. Therefore, it provides an additional semi-analytical approach for analyzing flutter derivatives of the bridge with streamlined cross-section to get its aerodynamic information.     

超音速气流中受热曲壁板的非线性颤振特性

基于von Karman 大变形理论及带有曲率修正的一阶活塞理论, 用Galerkin方法建立了超音速气流中受热二维曲壁板的非线性气动弹性运动方程; 采用牛顿迭代法计算得到由静气动载荷和热载荷引起的静气动弹性变形; 根据李雅谱诺夫间接法分析了壁板初始曲率与温升对颤振边界的影响; 对二维曲壁板的非线性气动弹性方程组进行数值积分求解,分析了动压参数对受热二维曲壁板分岔特性的影响, 给出了典型状态下曲壁板非线性颤振响应的时程图与相图. 分析结果表明对小初始曲率的曲壁板, 温升对其静气动弹性变形影响较大, 且随着温升的增加其颤振临界动压急剧减小; 对具有较大初始曲率的曲壁板, 温升对其静气动弹性变形的影响较弱, 且随着温升的增加颤振临界动压基本保持不变. 初始几何曲率与气动热效应使得曲壁板具有复杂的动力学特性, 不再像平壁板一样, 经过倍周期分岔进入混沌, 而会出现由静变形状态直接进入混沌运动的现象, 且在混沌运动区域中还会出现静态稳定点或谐波运动, 在大曲率情况下, 曲壁板不会产生混沌运动, 而是幅值在一定范围内的极限带振荡.