卫星大型旋转载荷动平衡控制应用研究

卫星携带的转动体的动不平衡对星体产生非期望的干扰力矩,进而影响卫星姿态,为此在卫星高品质运动控制中必须对转动体的动平衡特性进行控制。本文在建立卫星多挠性体姿态动力学模型的同时,分析了动平衡检测在转动体残余干扰力矩分析中的应用,提出了大型旋转载荷动平衡控制方法和流程,最后通过实际工程应用情况,验证了该动平衡控制方法和流程的实用性和有效性。     

对转双转子系统临界转速求解方法的改进与对比

为解决对转双转子航空发动机临界转速求解的难题,满足转子动力学设计需求,结合商业有限元软件,给出直接法和完全法两种临界转速求解方法,并结合算例详细论述了方法的理论基础和求解过程,最终以谐响应分析结果加以验证.研究结果表明:基于商业有限元软件特征值分析的完全法不需要任何振型信息,就可对正反进动曲线进行区分,极大简化了对转双...     

航空发动机高压转子试验模型的动力学相似设计

为了获得能够替代航空发动机转子原型进行动力学试验的模型，提出了发动机转子系统试验模型的动力学相似设计方法。采用相似理论推导了转子系统缩比模型的动力学相似准则，继而建立缩比模型与原型的动力学相似关系。通过动力学优化方法对转子缩比模型的动力学等效模型进行修正，目标函数由固有频率、临界转速和振型组成。以双转子发动机高压转子系统的动力学相似设计为例，通过有限元仿真对设计方法的有效性进行了验证。结果表明，在设计转速范围内，动力学相似模型与原型的前两阶临界转速相似误差分别为3.92%和0.86%，对应振型的相关性均大于0.98，且模型与原型的模态应变能分布基本一致。通过该方法获得的动力学相似模型能够较好地预测原型的动力学特性，并有效降低模型试验的成本和难度。     

Nonlinear interval analysis of rotor response with joints under uncertainties

The dynamic influence of joints in aero-engine rotor systems is investigated in this paper. Firstly, the tangential stiffness and loss factor are obtained from an isolated lap joint setup with dynamic excitation experiments. Also, the influence of the normal contact pressure and the excitation level are examined, which revel the uncertainty in joints. Then, the updated Thin Layer Elements (TLEs) method with fitted parameters based on the experiments is established to simulate the dynamic properties of joints on the interface. The response of the rotor subjected to unbalance excitation is calculated, and the results illustrate the effectiveness of the proposed method. Meanwhile, using the Chebyshev inclusion function and a direct iteration algorithm, a nonlinear interval analysis method is established to consider the uncertainty of parameters in joints. The accuracy is proved by comparison with results obtained using the Monte-Carlo method. Combined with the updated TLEs, the nonlinear Chebyshev method is successfully applied on a finite model of a rotor. The study shows that substantial attention should be paid to the dynamical design for the joint in rotor systems, the dynamic properties of joints under complex loading and the corresponding interval analysis method need to be intensively studied.     

Effects of tooth bending damage on the leakage performance and rotordynamic coefficients of labyrinth seals

Tooth bending damage resulting from an intense impact by the rotor sometimes occurs in the transient operation. To investigate the influence of after-damage clearance and tooth bending length on the leakage performance and rotordynamic coefficients of labyrinth seals, three tooth bending damages were taken into consideration, including the unbent tooth damage (abbreviated as Unbent), the partial tooth bending damage (abbreviated as Pbent) and the complete tooth bending damage (abbreviated as Cbent). The transient CFD solution was utilized to calculate the leakage flow rates and rotordynamic coefficients of labyrinth seals with clearances of 0.3, 0.4, 0.5, 0.6 mm for three tooth bending damages. The obtained result shows that the Unbent tooth damage leaks least while the Pbent tooth bending damage leaks most, and an increase of 6.1% for Cbent tooth bending damage and an increase of 19.4% for Pbent tooth bending damage are discovered at the tooth clearance of 0.6 mm in comparison with the Unbent tooth damage. Compared to the Unbent tooth damage, the effective damping for Pbent tooth bending damage and Cbent tooth bending damage is lower and drops by 9.7%–33.6% and 8.5%–22.6% respectively at the tooth clearance of 0.6 mm, suggesting that Pbent tooth bending damage or Cbent tooth bending damage tends to weaken the seal stability when compared to the Unbent tooth damage.     

Nonlinear characteristic investigation of magnetorheological damper-rotor system with local nonlinearity

Magnetorheological (MR) dampers show superior performance in reducing rotor vibration, but their high nonlinearity will cause nonsynchronous response, resulting in fatigue and instability of rotors. Herein, we are devoted to the investigation of the nonlinear characteristics of MR damper mounted on a flexible rotor. First, Reynolds equations with bilinear constitutive equations of MR fluid are employed to derive nonlinear oil film forces. Then, the Finite Element (FE) model of rotor system is developed, where the local nonlinear support forces produced by MR damper and its coupling effects with the rotor are considered. A hybrid numerical method is proposed to solve the nonlinear FE motion equations of the MR damper-rotor system. To validate the proposed model, a rotor test bench with two dual-coil MR dampers is constructed, upon which experimental studies on the dynamic characteristics of MR damper-rotor system are carried out. The effects of different system parameters, including rotational speed, excitation current and amount of unbalance, on nonlinear dynamic behaviors of MR damper-rotor system are evaluated. The results show that the system may appear chaos, jumping, and other complex nonlinear phenomena, and the level of the nonlinearity can be effectively alleviated by applying suitable excitation current and oil supply pressure.     

Integrated optimization analyses of aerodynamic/stealth characteristics of helicopter rotor based on surrogate model

Based on computational fluid dynamics(CFD)method,electromagnetic high-frequency method and surrogate model optimization techniques,an integration design method about aerodynamic/stealth has been established for helicopter rotor.The developed integration design method is composed of three modules:integrated grids generation(the moving-embedded grids for CFD solver and the blade grids for radar cross section(RCS)solver are generated by solving Poisson equations and folding approach),aerodynamic/stealth solver(the aerodynamic characteristics are simulated by CFD method based upon Navier–Stokes equations and Spalart–Allmaras(S–A)turbulence model),and the stealth characteristics are calculated by using a panel edge method combining the method of physical optics(PO),equivalent currents(MEC)and quasi-stationary(MQS),and integrated optimization analysis(based upon the surrogate model optimization technique with full factorial design(FFD)and radial basis function(RBF)),an integrated optimization analyses on aerodynamic/stealth characteristics of rotor are conducted.Firstly,the scattering characteristics of the rotor with different blade-tip swept and twist angles have been carried out,then time–frequency domain grayscale with strong scattering regions of rotor have been given.Meanwhile,the effects of swept-tip and twist angles on the aerodynamic characteristic of rotor have been performed.Furthermore,by choosing suitable object function and constraint condition,the compromised design about swept and twist combinations of rotor with high aerodynamic performances and low scattering characteristics has been given at last.     

Aerodynamic shape optimization for alleviating dynamic stall characteristics of helicopter rotor airfoil

In order to alleviate the dynamic stall effects in helicopter rotor, the sequential quadratic programming(SQP) method is employed to optimize the characteristics of airfoil under dynamic stall conditions based on the SC1095 airfoil. The geometry of airfoil is parameterized by the class-shape-transformation(CST) method, and the C-topology body-fitted mesh is then automatically generated around the airfoil by solving the Poisson equations. Based on the grid generation technology, the unsteady Reynolds-averaged Navier-Stokes(RANS) equations are chosen as the governing equations for predicting airfoil flow field and the highly-efficient implicit scheme of lower–upper symmetric Gauss–Seidel(LU-SGS) is adopted for temporal discretization. To capture the dynamic stall phenomenon of the rotor more accurately, the Spalart–Allmaras turbulence model is employed to close the RANS equations. The optimized airfoil with a larger leading edge radius and camber is obtained. The leading edge vortex and trailing edge separation of the optimized airfoil under unsteady conditions are obviously weakened, and the dynamic stall characteristics of optimized airfoil at different Mach numbers, reduced frequencies and angles of attack are also obviously improved compared with the baseline SC1095 airfoil. It is demonstrated that the optimized method is effective and the optimized airfoil is suitable as the helicopter rotor airfoil.     

航空发动机双转子系统临界转速求解方法

针对航空发动机转子动力学设计需求,将有限元法与振型筛选法相结合,提出适于工程应用的临界转速求解方法,结合算例详细论述了完全法和缩减法的求解过程。研究结果表明:两种求解方法的结果完全一致,而缩减法求解过程迭代次数少,更为方便;某双转子系统,顺转时存在5阶临界转速,而反转时在相同的转速范围内将增加n=3,n=4和n=5的反进动临界转速,临界转速的密集性导致反转转子的动力学设计更为复杂。     

Rotor airfoil aerodynamic design method and wind tunnel test verification

Well-designed airfoil is very important for high-performance rotor. This paper developed an efficient multi-objective and multi-constraint optimization design system for rotor airfoils based on RANS analysis, and verified the performance of the optimized airfoil. Using CRA09-A as the baseline rotor airfoil, the CRA09-B optimized rotor airfoil was designed successfully. Combined with the foundation of high-precision rotor airfoil stationary test technology, the CRA09-B and CRA09-A rotor airfoils were tested in the S3MA high-speed wind tunnel of ONERA. In order to correct the aerodynamic data, a single parameter linear wall pressure method is used to consider the tunnel effects. The results indicate that multi-objective and multi-constraint optimization design method developed in this study is reliable, and that CRA09-B optimized airfoil provides better stationary performance than CRA09-A airfoil in terms of maximum lift coefficient and lift-to-drag ratio.