机翼颤振的区间有限元分析(英文)

The influences of uncertainties in structural parameters on the flutter speed of wing are studied. On the basis of the deterministic flutter analysis model of wing, the uncertainties in structural parameters are considered and described by interval numbers. By virtue of first-order Taylor series expansion, the lower and upper bound curves of the transient decay rate coefficient versus wind velocity are given. So the interval estimation of the flutter critical wind speed of wing can be obtained, which is more reasonable than the point esti- mation obtained by the deterministic flutter analysis and provides the basis for the further non-probabilistic interval reliability analysis of wing flutter. The flow chart for interval fmite element model of flutter analysis of wing is given. The proposed interval finite element model and the stochastic finite element model for wing flutter analysis are compared by the examples of a three degrees of freedom airfoil and fuselage and a 15° sweptback wing, and the results have shown the effectiveness and feasibility of the presented model. The prominent advantage of the proposed interval finite element model is that only the bounds of uncertain parameters are required, and the probabilistic distribution densities or other statistical characteristics are not needed.     

Investigation of rotor control system loads

This paper concentrates on the aeroelasticity analysis of rotor blade and rotor control systems. A new multi-body dynamics model is established to predict both rotor pitch link loads and swashplate servo loads. Two helicopter rotors of UH-60A and SA349/2, both operating in two critical flight conditions, high-speed flight and high-thrust flight, are studied. The analysis shows good agreements with the flight test data and the calculation results using CAMRAD II. The mechanisms of rotor control loads are then analyzed in details based on the present predictions and the flight test data. In high-speed conditions, the pitch link loads are dominated by the integral of blade pitching moments, which are generated by cyclic pitch control. In high-thrust conditions, the positive pitching loads in the advancing side are caused by high collective pitch angle, and dynamic stall in the retreating side excites high-frequency responses. The swashplate servo loads are predominated by the rotor pitch link loads, and the inertia of the swashplate has significant effects on high-frequency harmonics of the servo loads.     

Aerodynamic Design and Analysis of a Low-reaction Axial Compressor Stage

There is introduced a new low-reaction, highly-loaded axial compressor design concept which is coupled with boundary layer suction method. The characteristic features of the concept are made clear through its comparison with the MIT boundary layer suction compressor. Also are pointed out the potential applications of this concept as well as its key technological problems. Based on this concept, a single-stage, low-reaction and low-speed axial compressor is constructed in association with analysis and computation of boundary layer suction on vanes with the aid of a three-dimensional numerical approach. The results attest to the effectiveness of this way to control separation in blade cascades by the boundary layer suction and the feasibility of this proposed design concept.     

Multi Objective Robust Active Vibration Control for Flexure Jointed Struts of Stewart Platforms via H∞ and μ Synthesis

Active vibration control is needed for future space telescopes, space laser communication and other precision sensitive payloads which require ultra-quiet environments. A Stewart platform based hybrid isolator with 6 hybrid struts is the effective system for active/passive vibration isolation over 5-250 Hz band. Using an identification transfer matrix of the Stewart platform, the coupling analysis of six channels is provided. A dynamics model is derived, and the rigid mode is removed to keep the signal of pointing control. Multi objective robust H∞ and μ synthesis strategies, based on singular values and structured singular values respectively, are presented, which simultaneously satisfy the low frequency pointing and high frequency disturbance rejection requirements and take account of the model uncertainty, parametric uncertainty and sensor noise. Then, by performing robust stability test, it is shown that the two controllers are robust to the uncertainties, the robust stability margin of H, controller is less than that of μ controller, but the order of μ controller is higher than that of H, controller, so the balanced controller reduction is provided. Additionally, the μ controller is compared with a PI controller. The time domain simulation of the μ controller indicates that the two robust control strategies are effective for keeping the pointing command and isolating the harmonic and stochastic disturbances.     

CALCULATION OF HELICOPTER ROTOR FLAPPING ANGLES AND COMPARISON WITH MEASURED DATA

Helicopter rotor flapping angles from hover to low-speed forward flight are calculated and compared with the measured data in this paper. The analytical method is based on a second order lifting-line/full-span free wake model as well as a fully coupled rotor trim model. It is shown that, in order to accurately predict the lateral flapping angle at low advance ratio, it is necessary to use free wake analysis to account for the highly non-uniform inflow induced by the distorted wake geometryat rotor disc plane.     

具有电磁约束阻尼层梁的振动主动控制研究计算(英文)

This paper investigates vibration control of beam through electro-magnetic constrained layer damping （EMCLD） which consists of electromagnet layer, permanent magnet layer and viscoelastic damping layer. When the coil of the electromagnet is electrified with proper control strategy, the electromagnet can exert magnetic force opposite to the direction of structural deformation so that the structural vibration is attenuated. A mathematical model is developed based on the equivalent current method to calculate the electromagnetic control force produced by EMCLD. The governing equations of the system are obtained using Hamilton＇s Principle and then reduced with the assumed-mode method. A simulation on vibration control of a cantilever beam is conducted under the velocity proportional feedback to demonstrate the energy dissipation capability of EMCLD, and the beam system with the same parameter is experimented. The results of experiment and simulation are compared and the results show that the EMCLD is an effective means for suppressing modal vibration. The results also indicate that the beam system has better control performance for larger control current. The EMCLD method presented in this paper provides an applicable and efficient tool for the vibration control of structures.     

Aeroelastic two-level optimization for preliminary design of wing structures considering robust constraints

An aeroelastic two-level optimization methodology for preliminary design of wing struc- tures is presented, in which the parameters for structural layout and sizes are taken as design vari- ables in the first-level optimization, and robust constraints in conjunction with conventional aeroelastic constraints are considered in the second-level optimization. A low-order panel method is used for aerodynamic analysis in the first-level optimization, and a high-order panel method is employed in the second-level optimization. It is concluded that the design of the abovementioned structural parameters of a wing can be improved using the present method with high efficiency. An improvement is seen in aeroelastic performance of the wing obtained with the present method when compared to the initial wing. Since these optimized structures are obtained after consideration of aerodynamic and structural uncertainties, they are well suited to encounter these uncertainties when they occur in reality.     

Numerical Investigation of Inlet Distortion on an Axial Flow Compressor Rotor with Circumferential Groove Casing Treatment

On the base of an assumed steady inlet circumferential total pressure distortion, three-dimensional time-dependent numerical simulations are conducted on an axial flow subsonic compressor rotor. The performances and flow fields of a compressor rotor, either casing treated or untreated, are investigated in detail either with or without inlet pressure distortion. Results show that the circumferential groove casing treatment can expand the operating range of the compressor rotor either with or without inlet pressure distortion at the expense of a drop in peak isentropic efficiency. The casing treatment is capable of weakening or even removing the tip leakage vortex effectively either with or without inlet distortion. In clean inlet circumstances, the enhancement and forward movement of tip leakage vortex cause the untreated compressor rotor to stall. By contrast, with circumferential groove casing, the serious flow separation on the suction surface leads to aerodynamic stalling eventually. In the presence of inlet pressure distortion, the blade loading changes from passage to passage as the distorted inflow sector is traversed. Similar to the clean inlet circumstances, with a smooth wall casing, the enhancement and forward movement of tip leakage vortex are still the main factors which lead to the compressor rotor stalling eventually. When the rotor works trader near stall conditions, the blockage resulting from the tip leakage vortex in all the passages is very serious. Especially in several passages, flow-spillage is observed. Compared to the clean inlet circumstances, circumferential groove casing treatment can also eliminate the low energy zone in the outer end wall region effectively.     

跨声速弯掠动叶压气机非定常流场的数值研究(英文)

The unsteady 3D flow fields in a single-stage transonic compressor under designed conditions are simulated numerically to investigate the effects of the curved rotors on the stage performance and the aerodynamic interaction between the blade rows. The results show that, compared to the compressor with unurved rotors, the compressor under scrutiny acquires remarkable increases in efficiency with significantly reduced amplitudes of the time-dependent fluctuation. The amplitude of the pressure fluctuation around the stator leading edge decreases at both endwalls, but increases at the mid-span in the curved rotors. The pressure fluctuation near the stator leading edge, therefore, becomes more uniform in the radial direction of this compressor. Except for the leading edge area, the pressure fluctuatinn amplitude declines remarkably in the tip region of stator surface downstream of the curved rotor, but hardly changes in the middle and at the hub.     

基于数字虚拟飞行的民用飞机纵向地面操稳特性评估

针对民用飞机设计方案纵向地面操稳特性的评估需求,面向适航标准的要求,提出了一种基于数字虚拟飞行的评估方法。基于适航条例要求提出了纵向地面操稳特性的量化判定准则,建立了飞机的地面运动模型和驾驶员操纵模型,以实现起降等特定地面运行任务的数字虚拟飞行,最终依据数字虚拟飞行结果和判定准则对飞机设计方案的地面操稳特性做出评估。应用此方法研究了某大型运输类飞机的纵向地面操稳特性。数字虚拟飞行结果表明:前翻倾向的严重情况发生在起降过程的高速滑行段,主轮刹车引起的机身前翻倾向是显著的,起落架纵向定位参数设计以及飞行进近参数选择均会对飞机的纵向地面操稳特性产生影响。     

直升机急拉杆机动飞行仿真建模与验证

 针对直升机大机动飞行仿真,建立了一个非线性的飞行动力学模型,考虑了翼型非定常/动态失速、机动飞行引起的动态尾迹畸变、桨叶弹性变形效应和发动机动态特性。采用基于有限元分析的挥舞-摆振-扭转耦合的弹性桨叶模型,并利用一种新的数值方法将旋翼/机体耦合运动方程表示为显式形式,整个飞行动力学模型表示为状态空间格式。以UH-60A直升机在高速飞行条件下的急拉杆机动飞行为例进行仿真计算,并与飞行试验数据进行对比验证。分析表明,仿真结果与试验结果吻合,高速飞行条件下机体抬头过程中前行桨叶非定常气动载荷的计算误差是引起旋翼和机体运动仿真误差的主要原因。     

直升机非线性运动方程及其数值分析

基于牛顿法建立了单旋翼带尾桨直升机的非线性运动方程，提出了一种求解直升机非线性运动方程的数值方法。建立方程时，用欧拉角描述直升机在空中的姿态，旋翼采用有挥舞铰且在铰上带有弹性约束的模型，桨盘处的入流采用线性模型，通过对一算例直升机的操纵响应进行数值仿真，将非线性与小扰动线性化的直升机运动方程进行了对比分析。结果表明：从操纵输入到随后的2s左右，非线性和小扰动线性化的直升机运动方程的计算结果非常接近，但在此之后，非线性运动方程的计算结果更能反映直升机的运动规律。     

Experimental and numerical research on the aerodynamics of unsteady moving aircraft

For the experimental determination of the dynamic wind tunnel data, a new combined motion test capability was developed at the German–Dutch Wind Tunnels DNW for their 3 m Low Speed Wind Tunnel NWB in Braunschweig, Germany, using a unique six degree-of-freedom test rig called ‘Model Positioning Mechanism’ (MPM) as an improved successor to the older systems. With that cutting-edge device, several transport aircraft configurations including a blended wing body configuration were tested in different modes of oscillatory motions roll, pitch and yaw as well as delta-wing geometries like X-31 equipped with remote controlled rudders and flaps to be able to simulate realistic flight maneuvers, e.g., a Dutch Roll. This paper describes the motivation behind these tests and the test setup and in addition gives a short introduction into time accurate maneuver-testing capabilities incorporating models with remote controlled control surfaces. Furthermore, the adaptation of numerical methods for the prediction of dynamic derivatives is described and some examples with the DLR-F12 configuration will be given. The calculations are based on RANS-solution using the finite volume parallel solution algorithm with an unstructured discretization concept (DLR TAU-code).     

基于全局灵敏度分析的侧向气动导数不确定性对侧向飞行载荷的影响

考虑飞行载荷计算中使用的气动导数存在不确定性,利用基于方差的全局灵敏度分析(GSA)方法,结合偏航机动的动力学模型,分析了侧向气动导数不确定性对侧向飞行载荷的影响。以某型飞机为例,运用该方法得到侧向气动导数的全局灵敏度排序。结果表明:侧滑角、阻尼贡献的垂尾侧向载荷及垂尾侧向总载荷受全机侧力系数对侧滑角导数的影响最大,受航向静稳定导数及方向舵操纵效能的影响次之;方向舵偏度贡献的垂尾侧向载荷只受全机侧力系数对方向舵偏度导数的影响;无尾飞机侧向载荷主要受航向静稳定导数、方向舵操纵效能及无尾飞机侧力系数对侧滑角导数的影响;偏航阻尼导数基本不影响各侧向飞行载荷。同时也验证了方法的有效性,对提高飞行载荷的计算精度有一定的指导意义。     

垂直/短距起降飞机的控制仿真技术

针对垂直/短距起降飞机的特点,通过在某型飞机上虚拟加装升力风扇系统,建立了垂直/短距起降飞机的动力学模型,提出了垂直起降阶段垂向、纵向、横向和航向的控制方式并完成了控制律的设计与验证。研究结果表明,所建立的垂直起降动力学模型能够描述垂直起降飞机的动力学特点,提出的控制方式和设计的控制律能够有效地控制飞机实现垂直升降、侧飞、偏航、前飞和后飞等运动,可用于垂直起降飞机的飞行品质、起降程序的设计和验证等相关方面的研究。     

共轴式直升机飞行动力学仿真数学模型研究

 在讨论共轴双旋翼气动干扰模型的基础上,建立了一种共轴式直升机飞行动力学仿真数学模型。以某型共轴式直升机为样例机,在定直水平飞行条件下进行配平计算,并与国外计算结果作了对比分析,两者基本一致。给出了共轴式直升机在不同操纵输入下的动态响应数值模拟结果。仿真计算表明,所建模型能够刻画共轴式直升机运动的基本特征。     

控制器非线性对飞机纵向飞行品质的影响分析

针对飞机控制系统执行机构具有非线性这一特点,选用某歼击机作为研究对象,分析各类非线性环节对飞机纵向短周期特性的影响。首先分析了死区、间隙和速率饱和等非线性环节的基本特性,建立了相应的包含非线性环节的飞行动力学模型,然后,基于Chalk准则计算不同非线性参数下的飞行品质指标,进行飞行品质评价。通过对不同飞控系统构型的对比研究,探讨了执行机构非线性对飞行品质的影响。研究结果表明,各非线性环节都将对系统的快速性产生不利影响,使系统有效上升时间增加,相位延迟增加。其中,间隙还会造成系统阻尼减小,可导致系统不稳并出现极限环振荡。     

飞机拦阻钩碰撞动力学和拦阻钩纵向阻尼器性能

 在考虑舰载机降落平台纵摇和横摇的基础上,建立了飞机拦阻钩六自由度碰撞反弹模型,得到了拦阻钩接触道面后反弹的动力学性能。分析了航母纵摇和横摇下拦阻钩碰撞反弹成因,并分别考虑了纵摇角和横摇角对拦阻钩反弹角速度及机身与道面给予拦阻钩碰撞冲量的影响;研究了拦阻钩碰撞后的反弹位移,在考虑拦阻索能顺利上钩的前提下,分析了拦阻钩纵向阻尼器的缓冲阻尼特性。结果表明：因航母横摇,碰撞后拦阻钩出现了左右的反转角速度及碰撞冲量;拦阻钩反弹后在自身重力作用下不能使拦阻索顺利上钩,在加入纵向阻尼器情况下,拦阻钩第1次反弹高度及回落时间均满足拦阻索上钩的条件。     

鸭式旋翼/机翼飞机悬停状态飞行动力学特性

针对鸭式旋翼/机翼(Canard Rotor/Wing,CRW)飞机独特的气动布局,常规的分析方法及经验公式很难准确地对CRW飞机进行飞行动力学研究,通过飞行辨识对CRW飞机悬停状态特性进行了研究。首先,设计了飞行试验并获得了高质量的飞行数据,基于频率响应对CRW飞机的状态空间模型进行了简化。然后,在频域内对飞机的动力学参数进行了拟合优化,获得了CRW飞机悬停状态的动力学模型,并用飞行数据对所建模型进行了验证。最后,用辨识所得参数与常规直升机悬停状态时的参数进行了对比。结果显示悬停时CRW飞机的操纵导数和阻尼导数均比常规直升机小,经分析,操纵导数的减小主要是CRW飞机独特的旋翼设计所致,阻尼导数减小的原因主要是旋翼气动影响以及鸭翼、平尾、垂尾的结构影响。动力学特性分析结果为CRW飞机旋翼模式总体设计的进一步优化提供了指引和参考,所建立的模型可用于控制系统设计。