Gas turbine rotors are complex dynamic systems with high-dimensional,discrete,and multi-source nonlinear coupling characteristics.Significant amounts of resources and time are spent during the process of solving dynam...Gas turbine rotors are complex dynamic systems with high-dimensional,discrete,and multi-source nonlinear coupling characteristics.Significant amounts of resources and time are spent during the process of solving dynamic characteristics.Therefore,it is necessary to design a lowdimensional model that can well reflect the dynamic characteristics of high-dimensional system.To build such a low-dimensional model,this study developed a dimensionality reduction method considering global order energy distribution by modifying the proper orthogonal decomposition theory.First,sensitivity analysis of key dimensionality reduction parameters to the energy distribution was conducted.Then a high-dimensional rotor-bearing system considering the nonlinear stiffness and oil film force was reduced,and the accuracy and the reusability of the low-dimensional model under different operating conditions were examined.Finally,the response results of a multi-disk rotor-bearing test bench were reduced using the proposed method,and spectrum results were then compared experimentally.Numerical and experimental results demonstrate that,during the dimensionality reduction process,the solution period of dynamic response results has the most significant influence on the accuracy of energy preservation.The transient signal in the transformation matrix mainly affects the high-order energy distribution of the rotor system.The larger the proportion of steady-state signals is,the closer the energy tends to accumulate towards lower orders.The low-dimensional rotor model accurately reflects the frequency response characteristics of the original high-dimensional system with an accuracy of up to 98%.The proposed dimensionality reduction method exhibits significant application potential in the dynamic analysis of highdimensional systems coupled with strong nonlinearities under variable operating conditions.展开更多
The protection effectiveness of traditional Lightning Strike Protection(LSP)for composite rotor blade of helicopter can be diminished due to the explosion risk in overlapping attachment under lightning strike,so a new...The protection effectiveness of traditional Lightning Strike Protection(LSP)for composite rotor blade of helicopter can be diminished due to the explosion risk in overlapping attachment under lightning strike,so a new protection method based on Air Breakdown and insulating adhesive layer(AB-LSP method)was designed to avoid it.In this study,a numerical method was developed to simulate the electrical breakdown,and verified by experiment results.Based on this method,a Finite Element Model(FEM)was established to investigate the effect of two factors(breakdown strength and initial ablation temperature of adhesive layer)on the LSP effectiveness.The results show that the breakdown strength impacts more to the ablation damage in composite than that of high-temperature resistance.Then,another FEM was established to predict the ablation damage by lightning strike in the AB-LSP method protected composite rotor blade.The mechanisms and potential key parameters(magnitude of lightning current,discharge channel location,adhesive layer thickness,and air gap width)that could affect the protection effectiveness were analyzed.The introduction of air breakdown changes the current conduction path and reduces explosion risk.After rational design,this method can offer effective lightning protection for composite helicopter rotor blade and other composite structures.展开更多
The Contra-Rotating Open Rotor(CROR)design confronts significant noise challenges despite being one of the possible options for future green aeroengines.To efficiently estimate the noise emitted from a CROR,a three-di...The Contra-Rotating Open Rotor(CROR)design confronts significant noise challenges despite being one of the possible options for future green aeroengines.To efficiently estimate the noise emitted from a CROR,a three-dimensional unsteady prediction model based on the meshless method is presented.The unsteady wake flow and the aerodynamic load fluctuations on the blade are solved through the viscous vortex particle method,the blade element momentum theory and vortex lattice method.Then,the acoustic field is obtained through the Farassat’s formulation 1A.Validation of this method is conducted on a CROR,and a mesh-based method,e.g.,Nonlinear Harmonic(NLH)method,is also employed for comparison.It is found that the presented method is three times faster than NLH method while maintaining a comparable precision.A thorough parametric analysis is also carried out to illustrate the effects of rotational speed,rotor-rotor spacing and rear rotor diameter on the noise level.The rotor speed is found to be the most influencing factor,and by optimizing the speed difference between the front and rear rotors,a notable noise reduction can be expected.The current findings not only contribute to a deeper comprehension of the CROR’s aeroacoustic properties but also offer an effective tool for engineering applications.展开更多
This paper analyzes the nonlinear dynamic characteristics and stability of Aero-Engine Dual-Rotor(AEDR)systems under high-frequency excitation,based on the Adaptive Harmonic Balance with the Asymptotic Harmonic Select...This paper analyzes the nonlinear dynamic characteristics and stability of Aero-Engine Dual-Rotor(AEDR)systems under high-frequency excitation,based on the Adaptive Harmonic Balance with the Asymptotic Harmonic Selection(AHB-AHS)method.A finite element dynamic equation for the AEDR system is introduced,considering complex nonlinearities of the intershaft bearing,unbalanced excitations,and high-frequency excitation.A solving strategy combining the AHB-AHS method and improved arclength continuation method is proposed to solve highdimensional dynamic equations containing complex nonlinearities and to track periodic solutions with parameter variations.The Floquet theory is used to analyze the types of bifurcation points in the system and the stability of periodic motions.The results indicate that high-frequency excitation can couple high-order and low-order modes,especially when the system undergoes superharmonic resonance.High-frequency excitation leads to more combination frequency harmonics,among which N_(f)ω_(1)-2ω_(2)dominates.Furthermore,changing the parameters(amplitude and frequency)of high-frequency excitation widens or shifts the unstable regions of the system.These findings contribute to understanding the mechanism of high-frequency excitation on aero-engines and demonstrate that the proposed AHB-AHS method is a powerful tool for analyzing highdimensional complex nonlinear dynamic systems under multi-frequency excitation.展开更多
Rotor blade is one of the most significant components of helicopters. But due to its highspeed rotation characteristics, it is difficult to collect the vibration signals during the flight stage.Moreover, sensors are h...Rotor blade is one of the most significant components of helicopters. But due to its highspeed rotation characteristics, it is difficult to collect the vibration signals during the flight stage.Moreover, sensors are highly susceptible to damage resulting in the failure of the measurement.In order to make signal predictions for the damaged sensors, an operational modal analysis(OMA) together with the virtual sensing(VS) technology is proposed in this paper. This paper discusses two situations, i.e., mode shapes measured by all sensors(both normal and damaged) can be obtained using OMA, and mode shapes measured by some sensors(only including normal) can be obtained using OMA. For the second situation, it is necessary to use finite element(FE) analysis to supplement the missing mode shapes of damaged sensor. In order to improve the correlation between the FE model and the real structure, the FE mode shapes are corrected using the local correspondence(LC) principle and mode shapes measured by some sensors(only including normal).Then, based on the VS technology, the vibration signals of the damaged sensors during the flight stage can be accurately predicted using the identified mode shapes(obtained based on OMA and FE analysis) and the normal sensors signals. Given the high degrees of freedom(DOFs) in the FE mode shapes, this approach can also be used to predict vibration data at locations without sensors. The effectiveness and robustness of the proposed method is verified through finite element simulation, experiment as well as the actual flight test. The present work can be further used in the fault diagnosis and damage identification for rotor blade of helicopters.展开更多
In this paper, an identification method to estimate the unbalances is introduced, which is based on the boundary element method (BEM). By using the vibration response measured at some points on the flexible rotor the ...In this paper, an identification method to estimate the unbalances is introduced, which is based on the boundary element method (BEM). By using the vibration response measured at some points on the flexible rotor the unbalances can be identified conveniently. Therefore, the rotor can be balanced without test runs.展开更多
Recently,the finite element method(FEM) has been commonly applied in the engineering analysis of rotor dynamics.Gyroscopic moments,rotary inertia,transverse shear deformation and gravity can be included in computation...Recently,the finite element method(FEM) has been commonly applied in the engineering analysis of rotor dynamics.Gyroscopic moments,rotary inertia,transverse shear deformation and gravity can be included in computational models of rotor-bearing systems.In this paper,a finite element model and its solution method are presented for the calculation of the dynamics of dual rotor systems.A typical structure with two rotor shafts is discussed and the procedure for obtaining the coupling motion equations of the subsystems is illustrated.A computer program is developed to solve critical speeds and to simulate the transient motion.The influence of gyroscopic moments on co-rotation and counter-rotation is analyzed,and the effect of the speed ratio on critical speed is studied.The dynamic characteristics under different conditions of increasing speed during start-up are demonstrated by comparison with transient nodal displacements.The presented model provides a complete foundation for further investigation of the dynamics of dual rotor systems.展开更多
A full-span free-wake method is coupled with an unsteady panel method to accurately predict the unsteady aerodynamics of helicopter rotor blades in hover and forward flight. The unsteady potential-based panel method i...A full-span free-wake method is coupled with an unsteady panel method to accurately predict the unsteady aerodynamics of helicopter rotor blades in hover and forward flight. The unsteady potential-based panel method is used to consider aerodynamics of finite thickness multi-bladed rotors, and the full-span free-wake method is applied to simulating dynamics of rotor wake. These methods are tightly coupled through trailing-edge Kutta condition and by converting doublet-wake panels to full-span vortex filaments. A velocity-field integration technique is also adopted to overcome singularity problem during the interaction between the rotor wake and blades. Helicopter rotors including Caradonna–Tung, UH-60A, and AH-1G rotors, are simulated in hover and forward flight to validate the accuracy of this approach. The predicted aerodynamic loads of rotor blades agree well with available measured data and computational fluid dynamics （CFD） results, and the unsteady dynamics of rotor wake is also well simulated. Compared to CFD, the present method obtains accurate results more efficiently and is suitable to rotorcraft aeroelastic analysis.展开更多
A hybrid Euler/full potential/Lagrangian wake method,based on single-blade simulation,for predicting unsteady aerodynamic flow around helicopter rotors in hover and forward flight has been developed.In this method,an ...A hybrid Euler/full potential/Lagrangian wake method,based on single-blade simulation,for predicting unsteady aerodynamic flow around helicopter rotors in hover and forward flight has been developed.In this method,an Euler solver is used to model the near wake evolution and transonic flow phenomena in the vicinity of the blade,and a full potential equation（FPE） is used to model the isentropic potential flow region far away from the rotor,while the wake effects of other blades and the far wake are incorporated into the flow solution as an induced inflow distribution using a Lagrangian based wake analysis.To further reduce the execution time,the computational fluid dynamics（CFD） solution and rotor wake analysis（including induced velocity up-date） are conducted parallelly,and a load balancing strategy is employed to account for the information exchange between two solvers.By the developed method,several hover and forward-flight cases on Caradonna-Tung and Helishape 7A rotors are per-formed.Good agreements of the loadings on blade surface with available measured data demonstrate the validation of the method.Also,the CPU time required for different computation runs is compared in the paper,and the results show that the pre-sent hybrid method is superior to conventional CFD method in time cost,and will be more efficient with the number of blades increasing.展开更多
Rotor airfoil design is investigated in this paper. There are many difficulties for this highdimensional multi-objective problem when traditional multi-objective optimization methods are used. Therefore, a multi-layer...Rotor airfoil design is investigated in this paper. There are many difficulties for this highdimensional multi-objective problem when traditional multi-objective optimization methods are used. Therefore, a multi-layer hierarchical constraint method is proposed by coupling principal component analysis(PCA) dimensionality reduction and e-constraint method to translate the original high-dimensional problem into a bi-objective problem. This paper selects the main design objectives by conducting PCA to the preliminary solution of original problem with consideration of the priority of design objectives. According to the e-constraint method, the design model is established by treating the two top-ranking design goals as objective and others as variable constraints. A series of bi-objective Pareto curves will be obtained by changing the variable constraints, and the favorable solution can be obtained by analyzing Pareto curve spectrum. This method is applied to the rotor airfoil design and makes great improvement in aerodynamic performance. It is shown that the method is convenient and efficient, beyond which, it facilitates decision-making of the highdimensional multi-objective engineering problem.展开更多
Current research in broken rotor bar （BRB） fault detection in induction motors is primarily focused on a high-frequency resolution analysis of the stator current. Compared with a discrete Fourier transformation, the...Current research in broken rotor bar （BRB） fault detection in induction motors is primarily focused on a high-frequency resolution analysis of the stator current. Compared with a discrete Fourier transformation, the parametric spectrum estimation technique has a higher frequency accuracy and resolution. However, the existing detection methods based on parametric spectrum estima- tion cannot realize online detection, owing to the large computational cost. To improve the efficiency of BRB fault detection, a new detection method based on the min-norm algorithm and least square estimation is proposed in this paper. First, the stator current is filtered using a band-pass filter and divided into short overlapped data windows. The min-norm algorithm is then applied to determine the fre- quencies of the fundamental and fault characteristic com- ponents with each overlapped data window. Next, based on the frequency values obtained, a model of the fault current signal is constructed. Subsequently, a linear least squares problem solved through singular value decomposition is designed to estimate the amplitudes and phases of the related components. Finally, the proposed method is applied to a simulated current and an actual motor, the results of which indicate that, not only parametric spectrum estimation technique.展开更多
Rotor system supported by nonlinear bearing such as squeeze film damper(SFD)is widely used in practice owing to its wide range of damping capacity and simplicity in structure.In this paper,an improved and effective In...Rotor system supported by nonlinear bearing such as squeeze film damper(SFD)is widely used in practice owing to its wide range of damping capacity and simplicity in structure.In this paper,an improved and effective Incremental transfer matrix method(ITMM)is first presented by combining ITMM and fast Fourier transform(FFT).Afterwards this method is applied to calculate the dynamic characteristics of a Jeffcott rotor system with SFD.The convergence dificulties incurred caused by strong nonlinearities of SFD has been dealt by adopting a control factor.It is found that for the more general boundary problems where the boundary conditions are not at input and output ends of a chain system,the supplementary equation is necessarily added.Additionally,the Floquet theory is used to analyze the stability and bifurcation type of the obtained periodic solution.The semi-analytical results,including the periodic solutions of the system,the bifurcation points and their types,are in good agreement with the numerical method.Furthermore,the involution mechanism of the quasi-periodic and chaotic motions near the first-order translational mode and the second order bending mode of this system is also clarified by this method with the aid of Floquet theory.展开更多
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 verifi...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 S3 MA 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.展开更多
In the present paper, a series of hierarchical warping functions is developed to analyze the static and dynamic problems of thin walled composite laminated helicopter rotors composed of several layers with single clos...In the present paper, a series of hierarchical warping functions is developed to analyze the static and dynamic problems of thin walled composite laminated helicopter rotors composed of several layers with single closed cell. This method is the development and extension of the traditional constrained warping theory of thin walled metallic beams, which had been proved very successful since 1940s. The warping distribution along the perimeter of each layer is expanded into a series of successively corrective warping functions with the traditional warping function caused by free torsion or free beading as the first term, and is assumed to be piecewise linear along the thickness direction of layers. The governing equations are derived based upon the variational principle of minimum potential energy for static analysis and Rayleigh Quotient for free vibration analysis. Then the hierarchical finite element method. is introduced to form a,. numerical algorithm. Both static and natural vibration problems of sample box beams axe analyzed with the present method to show the main mechanical behavior of the thin walled composite laminated helicopter rotor.展开更多
In order to develop a general calculating rotor’s torsional stiffness based on stiffness influence coefficient for different rotor assembling, the calculation method of the torsional stiffness influence coefficient o...In order to develop a general calculating rotor’s torsional stiffness based on stiffness influence coefficient for different rotor assembling, the calculation method of the torsional stiffness influence coefficient of equal thickness disc is researched in this paper at first. Then the torsional stiffness influence coefficient λ of equal thickness disc is fit to a binary curved face and a calculation equation is obtained based on a large quantity of calculating data, which lays the foundation for research on a general calculating method of rotor torsional stiffness. Thirdly a simplified calculation method for equivalent stiffness diameter of stepped equal thickness disc and cone disc in the steam turbine generators is suggested. Finally a general calculating program for calculating rotor’s torsional vibration features is developed, and the torsional vibration features of a verity of steam turbine rotors are calculated for verification. The calculating results show that stiffness influence coefficient λ of equal-thickness disc depends on parameters of B and H, as well as the stiffness influence coefficient λ;and discs with complex structure can be simplified to equal-thickness discs with little error by using the method suggested in this paper;error can be controlled within 1% when equivalent diameter of stiffness is calculated by this method.展开更多
A Large balancing weight may be obtained by calculation, when the ordinary influence coefficient method is used. In this paper, a new residual value balance method is recommended. The fundamental idea of the new metho...A Large balancing weight may be obtained by calculation, when the ordinary influence coefficient method is used. In this paper, a new residual value balance method is recommended. The fundamental idea of the new method is on the principle of getting the minimized weight on condition that the vibration values at various points are less than the admitted ones. Experimental results show that this method is effective for balancing the insensitive rotor and rotors with special dynamic characteristics.展开更多
An implicit higher ? order discontinuous Galerkin(DG) spatial discretization for the compressible Euler equations in a rotating frame of reference is presented and applied to a rotor in hover using hexahedral grids. I...An implicit higher ? order discontinuous Galerkin(DG) spatial discretization for the compressible Euler equations in a rotating frame of reference is presented and applied to a rotor in hover using hexahedral grids. Instead of auxiliary methods like grid adaptation,higher ? order simulations(fourth ? and fifth ? order accuracy) are adopted.Rigorous numerical experiments are carefully designed,conducted and analyzed. The results show generally excellent consistence with references and vigorously demonstrate the higher?order DG method's better performance in loading distribution computations and tip vortex capturing, with much fewer degrees of freedom(DoF). Detailed investigations on the outer boundary conditions for hovering rotors are presented as well. A simple but effective speed smooth procedure is developed specially for the DG method. Further results reveal that the rarely used pressure restriction for outlet speed has a considerable advantage over the extensively adopted vertical speed restriction.展开更多
In order to overcome the efficiency problem of the conventional gradient-based optimal design method,a highly-efficient viscous adjoint-based RANS equations method is applied to the aerodynamic optimal design of hover...In order to overcome the efficiency problem of the conventional gradient-based optimal design method,a highly-efficient viscous adjoint-based RANS equations method is applied to the aerodynamic optimal design of hovering rotor airfoil.The C-shaped body-fitted mesh is firstly automatically generated around the airfoil by solving the Poisson equations,and the Navier-Stokes(N-S)equations combined with Spalart-Allmaras(S-A)one-equation turbulence model are used as the governing equations to acquire the reliable flowfield variables.Then,according to multi-constrained characteristics of the optimization of high lift/drag ratio for hovering rotor airfoil,its corresponding adjoint equations,boundary conditions and gradient expressions are newly derived.On these bases,two representative rotor airfoils,NACA0012 airfoil and SC1095 airfoil,are selected as numerical examples to optimize their synthesized aerodynamic characteristics about lift/drag ratio in hover,and better aerodynamic performance of optimal airfoils are obtained compared with the baseline.Furthermore,the new designed rotor with the optimized rotor airfoil has better hover aerodynamic characteristics compared with the baseline rotor.In contrast to the baseline airfoils optimized by the finite difference method,it is demonstrated that the adjoint optimal algorithm itself is practical and highly-efficient for the aerodynamic optimization of hover rotor airfoil.展开更多
This paper implemented cooling configuration design on certain gas turbine HP rotor using parameterized method.It is convenient for complicated gas turbine blade modeling using parameters and also benefit for the geom...This paper implemented cooling configuration design on certain gas turbine HP rotor using parameterized method.It is convenient for complicated gas turbine blade modeling using parameters and also benefit for the geometry modify in later period.Parameterized modeling is the foundation of air cooling turbine blade design method engineering application.Mesh quality can be awarded when generated complicated cooling configuration blade grids,and also the increase of calculation error can arise by many mesh blocks.Film cooling and serpentine passage can effectively enhance the cooling effectiveness and protect blade.展开更多
For accurate aeroelastic analysis,the unsteady rotor flowfield is solved by computational fluid dynamics(CFD)module based on RANS/Euler equations and moving-embedded grid system,while computational structural dynamics...For accurate aeroelastic analysis,the unsteady rotor flowfield is solved by computational fluid dynamics(CFD)module based on RANS/Euler equations and moving-embedded grid system,while computational structural dynamics(CSD)module is introduced to handle blade flexibility.In CFD module,dual time-stepping algorithm is employed in temporal discretization,Jameson two-order central difference(JST)scheme is adopted in spatial discretization and B-L turbulent model is used to illustrate the viscous effect.The CSD module is developed based on Hamilton′s variational principles and moderate deflection beam theory.Grid deformation is implemented using algebraic method through coordinate transformations to achieve deflections with high quality and efficiency.A CFD/CSD loose coupling strategy is developed to transfer information between rotor flowfield and blade structure.The CFD and the CSD modules are verified seperately.Then the CFD/CSD loose coupling is adopted in airloads prediction of UH-60A rotor under high speed forward flight condition.The calculated results agree well with test data.Finally,effects of torsional stiffness properties on airloads of rotors with different tip swept angles(from 10° forward to 30° backward)are investigated.The results are evaluated through pressure distribution and airloads variation,and some meaningful conclusions are drawn the moderated shock wave strength and pressure gradient caused by varied tip swept angle and structural properties.展开更多
基金supported by the China Postdoctoral Science Foundation(No.2024M764171)the Postdoctoral Research Start-up Funds,China(No.AUGA5710027424)+1 种基金the National Natural Science Foundation of China(No.U2341237)the Development and construction funds for the School of Mechatronics Engineering of HIT,China(No.CBQQ8880103624)。
文摘Gas turbine rotors are complex dynamic systems with high-dimensional,discrete,and multi-source nonlinear coupling characteristics.Significant amounts of resources and time are spent during the process of solving dynamic characteristics.Therefore,it is necessary to design a lowdimensional model that can well reflect the dynamic characteristics of high-dimensional system.To build such a low-dimensional model,this study developed a dimensionality reduction method considering global order energy distribution by modifying the proper orthogonal decomposition theory.First,sensitivity analysis of key dimensionality reduction parameters to the energy distribution was conducted.Then a high-dimensional rotor-bearing system considering the nonlinear stiffness and oil film force was reduced,and the accuracy and the reusability of the low-dimensional model under different operating conditions were examined.Finally,the response results of a multi-disk rotor-bearing test bench were reduced using the proposed method,and spectrum results were then compared experimentally.Numerical and experimental results demonstrate that,during the dimensionality reduction process,the solution period of dynamic response results has the most significant influence on the accuracy of energy preservation.The transient signal in the transformation matrix mainly affects the high-order energy distribution of the rotor system.The larger the proportion of steady-state signals is,the closer the energy tends to accumulate towards lower orders.The low-dimensional rotor model accurately reflects the frequency response characteristics of the original high-dimensional system with an accuracy of up to 98%.The proposed dimensionality reduction method exhibits significant application potential in the dynamic analysis of highdimensional systems coupled with strong nonlinearities under variable operating conditions.
文摘The protection effectiveness of traditional Lightning Strike Protection(LSP)for composite rotor blade of helicopter can be diminished due to the explosion risk in overlapping attachment under lightning strike,so a new protection method based on Air Breakdown and insulating adhesive layer(AB-LSP method)was designed to avoid it.In this study,a numerical method was developed to simulate the electrical breakdown,and verified by experiment results.Based on this method,a Finite Element Model(FEM)was established to investigate the effect of two factors(breakdown strength and initial ablation temperature of adhesive layer)on the LSP effectiveness.The results show that the breakdown strength impacts more to the ablation damage in composite than that of high-temperature resistance.Then,another FEM was established to predict the ablation damage by lightning strike in the AB-LSP method protected composite rotor blade.The mechanisms and potential key parameters(magnitude of lightning current,discharge channel location,adhesive layer thickness,and air gap width)that could affect the protection effectiveness were analyzed.The introduction of air breakdown changes the current conduction path and reduces explosion risk.After rational design,this method can offer effective lightning protection for composite helicopter rotor blade and other composite structures.
基金the financial support from the National Natural Science Foundation of China(Nos.52276045 and 52206062)the Fundamental Research Funds for the Central Universities,China(Nos.3122019171,3122021087 and 3122022QD06).
文摘The Contra-Rotating Open Rotor(CROR)design confronts significant noise challenges despite being one of the possible options for future green aeroengines.To efficiently estimate the noise emitted from a CROR,a three-dimensional unsteady prediction model based on the meshless method is presented.The unsteady wake flow and the aerodynamic load fluctuations on the blade are solved through the viscous vortex particle method,the blade element momentum theory and vortex lattice method.Then,the acoustic field is obtained through the Farassat’s formulation 1A.Validation of this method is conducted on a CROR,and a mesh-based method,e.g.,Nonlinear Harmonic(NLH)method,is also employed for comparison.It is found that the presented method is three times faster than NLH method while maintaining a comparable precision.A thorough parametric analysis is also carried out to illustrate the effects of rotational speed,rotor-rotor spacing and rear rotor diameter on the noise level.The rotor speed is found to be the most influencing factor,and by optimizing the speed difference between the front and rear rotors,a notable noise reduction can be expected.The current findings not only contribute to a deeper comprehension of the CROR’s aeroacoustic properties but also offer an effective tool for engineering applications.
基金the financial support from the National Key R&D Program of China(No.2023YFE0125900)National Natural Science Foundation of China(Nos.12372008 and 12102234)+1 种基金Natural Science Foundation of Heilongjiang Province,China(No.YQ2022A008)Taif University,Saudi Arabia,for supporting this work through Project number(TU-DSPP-2024-73).
文摘This paper analyzes the nonlinear dynamic characteristics and stability of Aero-Engine Dual-Rotor(AEDR)systems under high-frequency excitation,based on the Adaptive Harmonic Balance with the Asymptotic Harmonic Selection(AHB-AHS)method.A finite element dynamic equation for the AEDR system is introduced,considering complex nonlinearities of the intershaft bearing,unbalanced excitations,and high-frequency excitation.A solving strategy combining the AHB-AHS method and improved arclength continuation method is proposed to solve highdimensional dynamic equations containing complex nonlinearities and to track periodic solutions with parameter variations.The Floquet theory is used to analyze the types of bifurcation points in the system and the stability of periodic motions.The results indicate that high-frequency excitation can couple high-order and low-order modes,especially when the system undergoes superharmonic resonance.High-frequency excitation leads to more combination frequency harmonics,among which N_(f)ω_(1)-2ω_(2)dominates.Furthermore,changing the parameters(amplitude and frequency)of high-frequency excitation widens or shifts the unstable regions of the system.These findings contribute to understanding the mechanism of high-frequency excitation on aero-engines and demonstrate that the proposed AHB-AHS method is a powerful tool for analyzing highdimensional complex nonlinear dynamic systems under multi-frequency excitation.
基金supported by grants from the High-Level Oversea Talent Introduction Plan,Chinathe Special Fund for Basic Scientific Research in Central Universities of China-Doctoral Research and Innovation Fund Project,China(No.3072023CFJ0206).
文摘Rotor blade is one of the most significant components of helicopters. But due to its highspeed rotation characteristics, it is difficult to collect the vibration signals during the flight stage.Moreover, sensors are highly susceptible to damage resulting in the failure of the measurement.In order to make signal predictions for the damaged sensors, an operational modal analysis(OMA) together with the virtual sensing(VS) technology is proposed in this paper. This paper discusses two situations, i.e., mode shapes measured by all sensors(both normal and damaged) can be obtained using OMA, and mode shapes measured by some sensors(only including normal) can be obtained using OMA. For the second situation, it is necessary to use finite element(FE) analysis to supplement the missing mode shapes of damaged sensor. In order to improve the correlation between the FE model and the real structure, the FE mode shapes are corrected using the local correspondence(LC) principle and mode shapes measured by some sensors(only including normal).Then, based on the VS technology, the vibration signals of the damaged sensors during the flight stage can be accurately predicted using the identified mode shapes(obtained based on OMA and FE analysis) and the normal sensors signals. Given the high degrees of freedom(DOFs) in the FE mode shapes, this approach can also be used to predict vibration data at locations without sensors. The effectiveness and robustness of the proposed method is verified through finite element simulation, experiment as well as the actual flight test. The present work can be further used in the fault diagnosis and damage identification for rotor blade of helicopters.
文摘In this paper, an identification method to estimate the unbalances is introduced, which is based on the boundary element method (BEM). By using the vibration response measured at some points on the flexible rotor the unbalances can be identified conveniently. Therefore, the rotor can be balanced without test runs.
文摘Recently,the finite element method(FEM) has been commonly applied in the engineering analysis of rotor dynamics.Gyroscopic moments,rotary inertia,transverse shear deformation and gravity can be included in computational models of rotor-bearing systems.In this paper,a finite element model and its solution method are presented for the calculation of the dynamics of dual rotor systems.A typical structure with two rotor shafts is discussed and the procedure for obtaining the coupling motion equations of the subsystems is illustrated.A computer program is developed to solve critical speeds and to simulate the transient motion.The influence of gyroscopic moments on co-rotation and counter-rotation is analyzed,and the effect of the speed ratio on critical speed is studied.The dynamic characteristics under different conditions of increasing speed during start-up are demonstrated by comparison with transient nodal displacements.The presented model provides a complete foundation for further investigation of the dynamics of dual rotor systems.
文摘A full-span free-wake method is coupled with an unsteady panel method to accurately predict the unsteady aerodynamics of helicopter rotor blades in hover and forward flight. The unsteady potential-based panel method is used to consider aerodynamics of finite thickness multi-bladed rotors, and the full-span free-wake method is applied to simulating dynamics of rotor wake. These methods are tightly coupled through trailing-edge Kutta condition and by converting doublet-wake panels to full-span vortex filaments. A velocity-field integration technique is also adopted to overcome singularity problem during the interaction between the rotor wake and blades. Helicopter rotors including Caradonna–Tung, UH-60A, and AH-1G rotors, are simulated in hover and forward flight to validate the accuracy of this approach. The predicted aerodynamic loads of rotor blades agree well with available measured data and computational fluid dynamics （CFD） results, and the unsteady dynamics of rotor wake is also well simulated. Compared to CFD, the present method obtains accurate results more efficiently and is suitable to rotorcraft aeroelastic analysis.
基金National Natural Science Foundation of China (10872094)
文摘A hybrid Euler/full potential/Lagrangian wake method,based on single-blade simulation,for predicting unsteady aerodynamic flow around helicopter rotors in hover and forward flight has been developed.In this method,an Euler solver is used to model the near wake evolution and transonic flow phenomena in the vicinity of the blade,and a full potential equation（FPE） is used to model the isentropic potential flow region far away from the rotor,while the wake effects of other blades and the far wake are incorporated into the flow solution as an induced inflow distribution using a Lagrangian based wake analysis.To further reduce the execution time,the computational fluid dynamics（CFD） solution and rotor wake analysis（including induced velocity up-date） are conducted parallelly,and a load balancing strategy is employed to account for the information exchange between two solvers.By the developed method,several hover and forward-flight cases on Caradonna-Tung and Helishape 7A rotors are per-formed.Good agreements of the loadings on blade surface with available measured data demonstrate the validation of the method.Also,the CPU time required for different computation runs is compared in the paper,and the results show that the pre-sent hybrid method is superior to conventional CFD method in time cost,and will be more efficient with the number of blades increasing.
基金supported by the National Natural Science Foundation of China (No. 11402288 and 11372254)the National Basic Research Program of China (No. 2014CB744804)
文摘Rotor airfoil design is investigated in this paper. There are many difficulties for this highdimensional multi-objective problem when traditional multi-objective optimization methods are used. Therefore, a multi-layer hierarchical constraint method is proposed by coupling principal component analysis(PCA) dimensionality reduction and e-constraint method to translate the original high-dimensional problem into a bi-objective problem. This paper selects the main design objectives by conducting PCA to the preliminary solution of original problem with consideration of the priority of design objectives. According to the e-constraint method, the design model is established by treating the two top-ranking design goals as objective and others as variable constraints. A series of bi-objective Pareto curves will be obtained by changing the variable constraints, and the favorable solution can be obtained by analyzing Pareto curve spectrum. This method is applied to the rotor airfoil design and makes great improvement in aerodynamic performance. It is shown that the method is convenient and efficient, beyond which, it facilitates decision-making of the highdimensional multi-objective engineering problem.
基金Supported by National Natural Science Foundation of China(Grant No.51607180)
文摘Current research in broken rotor bar （BRB） fault detection in induction motors is primarily focused on a high-frequency resolution analysis of the stator current. Compared with a discrete Fourier transformation, the parametric spectrum estimation technique has a higher frequency accuracy and resolution. However, the existing detection methods based on parametric spectrum estima- tion cannot realize online detection, owing to the large computational cost. To improve the efficiency of BRB fault detection, a new detection method based on the min-norm algorithm and least square estimation is proposed in this paper. First, the stator current is filtered using a band-pass filter and divided into short overlapped data windows. The min-norm algorithm is then applied to determine the fre- quencies of the fundamental and fault characteristic com- ponents with each overlapped data window. Next, based on the frequency values obtained, a model of the fault current signal is constructed. Subsequently, a linear least squares problem solved through singular value decomposition is designed to estimate the amplitudes and phases of the related components. Finally, the proposed method is applied to a simulated current and an actual motor, the results of which indicate that, not only parametric spectrum estimation technique.
文摘Rotor system supported by nonlinear bearing such as squeeze film damper(SFD)is widely used in practice owing to its wide range of damping capacity and simplicity in structure.In this paper,an improved and effective Incremental transfer matrix method(ITMM)is first presented by combining ITMM and fast Fourier transform(FFT).Afterwards this method is applied to calculate the dynamic characteristics of a Jeffcott rotor system with SFD.The convergence dificulties incurred caused by strong nonlinearities of SFD has been dealt by adopting a control factor.It is found that for the more general boundary problems where the boundary conditions are not at input and output ends of a chain system,the supplementary equation is necessarily added.Additionally,the Floquet theory is used to analyze the stability and bifurcation type of the obtained periodic solution.The semi-analytical results,including the periodic solutions of the system,the bifurcation points and their types,are in good agreement with the numerical method.Furthermore,the involution mechanism of the quasi-periodic and chaotic motions near the first-order translational mode and the second order bending mode of this system is also clarified by this method with the aid of Floquet theory.
基金supported by the National Natural Science Foundation of China(No.11902335)。
文摘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 S3 MA 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.
基金The project supported by the National Natural Science Foundation of China (19932030)
文摘In the present paper, a series of hierarchical warping functions is developed to analyze the static and dynamic problems of thin walled composite laminated helicopter rotors composed of several layers with single closed cell. This method is the development and extension of the traditional constrained warping theory of thin walled metallic beams, which had been proved very successful since 1940s. The warping distribution along the perimeter of each layer is expanded into a series of successively corrective warping functions with the traditional warping function caused by free torsion or free beading as the first term, and is assumed to be piecewise linear along the thickness direction of layers. The governing equations are derived based upon the variational principle of minimum potential energy for static analysis and Rayleigh Quotient for free vibration analysis. Then the hierarchical finite element method. is introduced to form a,. numerical algorithm. Both static and natural vibration problems of sample box beams axe analyzed with the present method to show the main mechanical behavior of the thin walled composite laminated helicopter rotor.
文摘In order to develop a general calculating rotor’s torsional stiffness based on stiffness influence coefficient for different rotor assembling, the calculation method of the torsional stiffness influence coefficient of equal thickness disc is researched in this paper at first. Then the torsional stiffness influence coefficient λ of equal thickness disc is fit to a binary curved face and a calculation equation is obtained based on a large quantity of calculating data, which lays the foundation for research on a general calculating method of rotor torsional stiffness. Thirdly a simplified calculation method for equivalent stiffness diameter of stepped equal thickness disc and cone disc in the steam turbine generators is suggested. Finally a general calculating program for calculating rotor’s torsional vibration features is developed, and the torsional vibration features of a verity of steam turbine rotors are calculated for verification. The calculating results show that stiffness influence coefficient λ of equal-thickness disc depends on parameters of B and H, as well as the stiffness influence coefficient λ;and discs with complex structure can be simplified to equal-thickness discs with little error by using the method suggested in this paper;error can be controlled within 1% when equivalent diameter of stiffness is calculated by this method.
文摘A Large balancing weight may be obtained by calculation, when the ordinary influence coefficient method is used. In this paper, a new residual value balance method is recommended. The fundamental idea of the new method is on the principle of getting the minimized weight on condition that the vibration values at various points are less than the admitted ones. Experimental results show that this method is effective for balancing the insensitive rotor and rotors with special dynamic characteristics.
基金co-supported by the National High Technology Research and Development Program of China(No.2015AA015303)the National Natural Science Foundation of China(No.11272152)+1 种基金the Aeronautical Science Foundation of China(No.20152752033)the Open Project of Key Laboratory of Aerodynamic Noise Control
文摘An implicit higher ? order discontinuous Galerkin(DG) spatial discretization for the compressible Euler equations in a rotating frame of reference is presented and applied to a rotor in hover using hexahedral grids. Instead of auxiliary methods like grid adaptation,higher ? order simulations(fourth ? and fifth ? order accuracy) are adopted.Rigorous numerical experiments are carefully designed,conducted and analyzed. The results show generally excellent consistence with references and vigorously demonstrate the higher?order DG method's better performance in loading distribution computations and tip vortex capturing, with much fewer degrees of freedom(DoF). Detailed investigations on the outer boundary conditions for hovering rotors are presented as well. A simple but effective speed smooth procedure is developed specially for the DG method. Further results reveal that the rarely used pressure restriction for outlet speed has a considerable advantage over the extensively adopted vertical speed restriction.
文摘In order to overcome the efficiency problem of the conventional gradient-based optimal design method,a highly-efficient viscous adjoint-based RANS equations method is applied to the aerodynamic optimal design of hovering rotor airfoil.The C-shaped body-fitted mesh is firstly automatically generated around the airfoil by solving the Poisson equations,and the Navier-Stokes(N-S)equations combined with Spalart-Allmaras(S-A)one-equation turbulence model are used as the governing equations to acquire the reliable flowfield variables.Then,according to multi-constrained characteristics of the optimization of high lift/drag ratio for hovering rotor airfoil,its corresponding adjoint equations,boundary conditions and gradient expressions are newly derived.On these bases,two representative rotor airfoils,NACA0012 airfoil and SC1095 airfoil,are selected as numerical examples to optimize their synthesized aerodynamic characteristics about lift/drag ratio in hover,and better aerodynamic performance of optimal airfoils are obtained compared with the baseline.Furthermore,the new designed rotor with the optimized rotor airfoil has better hover aerodynamic characteristics compared with the baseline rotor.In contrast to the baseline airfoils optimized by the finite difference method,it is demonstrated that the adjoint optimal algorithm itself is practical and highly-efficient for the aerodynamic optimization of hover rotor airfoil.
基金Sponsored by the National Natural Science Foundation of China(Grant No. 50476028)
文摘This paper implemented cooling configuration design on certain gas turbine HP rotor using parameterized method.It is convenient for complicated gas turbine blade modeling using parameters and also benefit for the geometry modify in later period.Parameterized modeling is the foundation of air cooling turbine blade design method engineering application.Mesh quality can be awarded when generated complicated cooling configuration blade grids,and also the increase of calculation error can arise by many mesh blocks.Film cooling and serpentine passage can effectively enhance the cooling effectiveness and protect blade.
文摘For accurate aeroelastic analysis,the unsteady rotor flowfield is solved by computational fluid dynamics(CFD)module based on RANS/Euler equations and moving-embedded grid system,while computational structural dynamics(CSD)module is introduced to handle blade flexibility.In CFD module,dual time-stepping algorithm is employed in temporal discretization,Jameson two-order central difference(JST)scheme is adopted in spatial discretization and B-L turbulent model is used to illustrate the viscous effect.The CSD module is developed based on Hamilton′s variational principles and moderate deflection beam theory.Grid deformation is implemented using algebraic method through coordinate transformations to achieve deflections with high quality and efficiency.A CFD/CSD loose coupling strategy is developed to transfer information between rotor flowfield and blade structure.The CFD and the CSD modules are verified seperately.Then the CFD/CSD loose coupling is adopted in airloads prediction of UH-60A rotor under high speed forward flight condition.The calculated results agree well with test data.Finally,effects of torsional stiffness properties on airloads of rotors with different tip swept angles(from 10° forward to 30° backward)are investigated.The results are evaluated through pressure distribution and airloads variation,and some meaningful conclusions are drawn the moderated shock wave strength and pressure gradient caused by varied tip swept angle and structural properties.
