Numerical simulation of a rotary engine primary compressor impacted by bird

In order to examine the process of a rotary engine primary compressor impacted by bird, a finite element model of a bird impacted on plate is developed with the explicit code PAM-CRASH. The smooth particles hydrodynamic (SPH) method is used to simulate the bird because of the SPH method showing no signs of instability and correctly modeling the breaking-up of the bird into particles. Good agreement between the simulation results and experimental results indicates that the numerical method of bird strike used in the present paper is reasonable. Then a rotary engine primary compressor impacted by three different configurations bird named straight-ended cylinder bird, quadrangular bird, hemispherical-ended bird are investigated using the numerical simulation method. It is found that the whole process of bird strike sustained about 3.5 ms and the bird is slashed by blade during the strike. The geometry configuration of bird affected the displacement and von Mises stress of some blades severely, just because the breaking bird’s mass is affected by the bird’s configuration. In the event of bird striking on the site of"up"some blades may develop plastic deformation and it is very adverse for the safety work of the engine.     

鸟体撞击结构过程的相似律研究(英文)

With dimensional analysis and similarity theory, the model similarity law of aircraft structures trader bird impact load is investigated. Numerical calculations by means of nonlinear dynamic software ANSYS/LS-DYNA are conducted on the finite element models constructed with different scaling factors. The influence of strain rate on the model similarity law is found to be dependent on the strain rate sensitivity of materials and scale factors. Specifically, materials that are not sensitive to strain rate obey the model similarity law in the bird impact process. The conclusions obtained are supposed to provide a theoretical basis for the experimental work of bird impact on aircraft structure.     

Study on Titanium Alloy TC4 Ballistic Penetration Resistance Part I:Ballistic Impact Tests

Ballistic impact test of different-scale casings is an efficient way to demonstrate the casing containment capability at the preliminary design stage of the engine. For the sake of studying the titanium alloy TC4 casing performance, the ballistic tests of flat and curved simulation casing are implemented by using two flat blades of different sizes as the projectile. The impact mechanism and failure of the target are discussed. Impact of the projectile is a highly nonlinear transient process with the large deformation of the target. On the impact, failures of the flat casing and the subscale casing are similar, concluding two parts, the global dishing and localized ductile tearing. The main localized failure mode combines plugging (shear) and petaling (shear) if the projectile perforates or penetrates, while crater (shear) if the projectile rebounds. The ballistic limit equation is verified by the test data and the results show that this empirical equation could be a practical way to estimate the critical velocity.     

Numerical evaluation of tandem rotor for highly loaded transonic fan

Transonic tandem rotor was designed for highly loaded fan at a corrected tip speed of 381m/s and another conventional rotor was designed as a baseline to evaluate the loading superiority of tandem rotor with three-dimensional (3-D) numerical simulation. The aft blade solidity and its impact on total loading level were studied in depth. The result indicates that tandem rotor has potential to achieve higher loading level and attain favorable aerodynamic performance in a wide range of loading coefficient 0.55~0.68, comparing with the conventional rotor which produced a total pressure ratio of 2.0 and loading coefficient of 0.42.      

Blade containment evaluation of civil aircraft engines

The potential hazard resulting from uncontained turbine engine rotor blade failure has always been the long-term concern of each aero engine manufacturer, and to fully contain the failed blades under critical operating conditions is also one of the most important considerations to meet the rotor integrity requirements. Usually, there are many factors involving the engine containment capability which need to be reviewed during the engine design phases, such as case thickness, rotor support structure, blade weight and shape, etc. However, the premier method to demonstrate the engine containment capability is the fan blade-off test and margin of safety (MS) analysis. Based on a concrete engine model, this paper aims to explain the key points of aero engine containment requirements in FAR Part 33, and introduces the implementation of MS analysis and fan blade-off test in the engine airworthiness certification. Through the introduction, it would be greatly helpful to the industrial community to evaluate the engine containment capability and prepare the final test demonstration in engine certification procedure.     

Applications of numerical optimization techniques to design of axial compressor blades

This paper describes the shape optimization of NASA rotor 37 and rotor and stator blades in a single-stage transonic axial compressor.Shape optimization of the blades operating at the design flow condition has been performed using the response surface method and three-dimensional Navier-Stokes analysis.Thin-layer approximation is introduced to the Navier-Stokes equations,and an explicit Runge-Kutta scheme is used to solve the governing equations.The three design variables,blade sweep,lean and skew,are introduced to optimize the three-dimensional stacking line of the blades.The objective function of the shape optimization is an adiabatic efficiency.Throughout the optimization of rotor and stator blades, optimal blade shape can be obtained.It is noted the increase of adiabatic efficiency by optimization of the blade shape with the stacking line in the single-stage transonic axial compressor is more effective in a rotor blade rather than a stator blade because of the large deformation of blade shape in the stator blade.      

GA-Based Model Predictive Control of Semi-Active Landing Gear

Semi-active landing gear can provide good performance of both landing impact and taxi situation, and has the ability for adapting to various ground conditions and operational conditions. A kind of Nonlinear Model Predictive Control algorithm （NMPC） for semi-active landing gears is developed in this paper. The NMPC algorithm uses Genetic Algorithm （GA） as the optimization technique and chooses damping performance of landing gear at touch down to be the optimization object. The valve＇s rate and magnitude limitations are also considered in the controller＇s design. A simulation model is built for the semi-active landing gear＇s damping process at touchdown. Drop tests are carried out on an experimental passive landing gear systerm to validate the parameters of the simulation model. The result of numerical simulation shows that the isolation of impact load at touchdown can be significantly improved compared to other control algorithms. The strongly nonlinear dynamics of semi-active landing gear coupled with control valve＇s rate and magnitude limitations are handled well with the proposed controller.     

Trajectory control strategy of cathodes in blisk electrochemical machining

A turbine blisk, which combines blades and a disk together, is one of the most important components of an aero engine. In the process of blisk electrochemical machining (ECM), the sheet cathode, which is usually used as a tool electrode, has a complicated structure. In addition to that, the channel between the adjacent blades is narrow and twisted, so interference is apt to happen when the sheet cathode feeds into the channel. Therefore, it is important to choose suitable trajectory control strategy. In this paper, a new trajectory control strategy of the sheet cathode is presented and corresponding simulation analysis is conducted on the basis of an actual blisk model. The simulation results demonstrate that the sheet cathode can feed into the channel by a spatial line trajectory without interference. Moreover, the verification experiments are carried out according to the simulation. The experimental results show that the cathode can move into the channel without interference. It is verified that the new trajectory control strategy is correct and can be used in the blisk ECM process successfully.     

NON-DIMENSIONAL DESIGN PARAMETERS FOR FOD TOLERANT FAN BLADES

Non-dimensional design concept for FOD tolerant fan blades is introduced based on the analyses of simplified impact models. The fan blades arc idealized as either beams or plates of elastic or rigid-plastic materials. The case of constant force impact as well as that of mass impact is analyzed. The centrifugal force effects are also considered in the beam models. The critical fracture conditions arc shown in simple npn-dimensional formulae or diagrams for each case.     

A second-order numerical method for elliptic equations with singular sources using local flter

The presence of Dirac delta function in differential equation can lead to a discontinuity,which may degrade the accuracy of related numerical methods.To improve the accuracy,a secondorder numerical method for elliptic equations with singular sources is introduced by employing a local kernel flter.In this method,the discontinuous equation is convoluted with the kernel function to obtain a more regular one.Then the original equation is replaced by this fltered equation around the singular points,to obtain discrete numerical form.The unchanged equations at the other points are discretized by using a central difference scheme.1D and 2D examples are carried out to validate the correctness and accuracy of the present method.The results show that a second-order of accuracy can be obtained in the fltering framework with an appropriate integration rule.Furthermore,the present method does not need any jump condition,and also has extremely simple form that can be easily extended to high dimensional cases and complex geometry.     

Numerical study of a 3.5-stage axial compressor at on-and off-design conditions

Numerical investigation is conducted on a 3.5-stage axial compressor,on which numerous experimental projects were carried out at the Institute during the last years and an experimental database was established.In the current study five on-and off-design operating points are simulated using a RANS solver and the results are compared with the measurement.The result shows that the compressor performance can be qualitatively predicted by the mixing-plane method.Better agreement is obtained for the on-design operating point.However,as the flow unsteadiness is insufficiently considered,the numerical method produces end-wall low-speed flow layers accumulated with the flow passing through the passage,which is in no good agreement with the experimental data.In the numerical simulation the rotor rows receive less work and this difference from the measurement increases with the rotational speed.In contrast,the stator rows increase the pressure more efficiently than the measurement.In the simulation the flow in the last stator row tends more to separate on the pressure side of the blade.For the operating points close to the surge line,the predicted separation is more intense than the experimental observation.But for the operating points close to the choke,the separation is suppressed.      

Mass and performance optimization of an airplane wing leading edge structure against bird strike using Taguchi-based grey relational analysis

Collisions between birds and aircraft are one of the most dangerous threats to flight safety. In this study, smoothed particles hydrodynamics (SPH) method is used for simulating the bird strike to an airplane wing leading edge structure. In order to verify the model, first, experiment of bird strike to a flat aluminum plate is simulated, and then bird impact on an airplane wing lead-ing edge structure is investigated. After that, considering dimensions of wing internal structural components like ribs, skin and spar as design variables, we try to minimize structural mass and wing skin deformation simultaneously. To do this, bird strike simulations to 18 different wing structures are made based on Taguchi’s L18 factorial design of experiment. Then grey relational analysis is used to minimize structural mass and wing skin deformation due to the bird strike. The analysis of variance (ANOVA) is also applied and it is concluded that the most significant parameter for the performance of wing structure against impact is the skin thickness. Finally, a validation simu-lation is conducted under the optimal condition to show the improvement of performance of the wing structure.     

Study on ballistic penetration resistance of titanium alloy TC4, Part II: Numerical analysis

Enhancing containment capability and reducing weight are always great concerns in the design of casings. Ballistic tests can help to mitigate a catastrophic event after a blade out, yet taking time and costing money. A wise way is to hunt for a validated numerical simulation technology, through which the material dynamic behavior over the strain rate range in the ballistic tests should be represented and reasonable failure strain should be defined. The simulation results show that the validation of the numerical simulation technology based on the test data can accurately estimate the absorption energy, describe the physical process and failure mode during the penetration, as well as the failure mechanism. It is found that energy dissipation of projectiles is in manner of compression stage, energy conversion stage, and interactive scrap stage. An effect indicator is proposed, where the factors of critical velocity including impact orientation and mass of projectiles and thickness of casings are considered. The critical velocity presents a linear relation with the effect indicator, which implies the critical velocity obtained by the flat casing could underestimate the capability of the real casing.     

Effect of blade tip winglet on the performance of a highly loaded transonic compressor rotor

The tip leakage flow has an important influence on the performance of transonic com-pressor. Blade tip winglet has been proved to be an effective method to control the tip leakage flow in compressor, while the physical mechanisms of blade tip winglet have been poorly understood. A numerical study for a highly loaded transonic compressor rotor has been conducted to understand the effect of varying the location of blade tip winglet on the performance of the rotor. Two kinds of tip winglet were designed and investigated. The effects of blade tip winglet on the compressor over-all performance, stability and tip flow structure were presented and discussed. It is found that the interaction of the tip winglet with the flow in the tip region is different when the winglet is located at suction-side or pressure-side of the blade tip. Results indicate that the suction-side winglet (SW) is ineffective to improve the performance of compressor rotor. In addition, a significant stall range extension equivalent to 33.74% with a very small penalty in efficiency can be obtained by the pressure-side winglet (PW). An attempt has been made to explain the fundamental mechanisms of blade tip winglet in detail.     

Dynamic Response Analysis of Bird Strike on Aircraft Windshield Based on Damage-modified Nonlinear Viscoelastic Constitutive Relation

Damage-modified nonlinear viscoelastic constitutive equation and failure criterion are introduced and the three-dimensional incremental forms are deduced based on the updated Lagrangian approach. A simple tensile test model and a split Hopkinson pressure bar model are built to verify the accuracy of the subroutine implemented within the non-linear finite element program LS-DYNA. A numerical model of bird strike on windshield is established to study the responses of windshield under three different bird velocities at three sites. The bird is represented by a cylinder with a hemisphere at each end and the contact-impact coupling algorithm is used in this study. It is found that the implemented subroutine can properly describe the mechanical behavior of polymethyl methacrylate under low and high strain rates and large deformation, and can be used validly.     

Cathode design investigation based on iterative correction of predicted profile errors in electrochemical machining of compressor blades

Electrochemical machining (ECM) is an effective and economical manufacturing method for machining hard-to-cut metal materials that are often used in the aerospace field. Cathode design is very complicated in ECM and is a core problem influencing machining accuracy, especially for complex profiles such as compressor blades in aero engines. A new cathode design method based on iterative correction of predicted profile errors in blade ECM is proposed in this paper. A math-ematical model is first built according to the ECM shaping law, and a simulation is then carried out using ANSYS software. A dynamic forming process is obtained and machining gap distributions at different stages are analyzed. Additionally, the simulation deviation between the prediction profile and model is improved by the new method through correcting the initial cathode profile. Further-more, validation experiments are conducted using cathodes designed before and after the simulation correction. Machining accuracy for the optimal cathode is improved markedly compared with that for the initial cathode. The experimental results illustrate the suitability of the new method and that it can also be applied to other complex engine components such as diffusers.     

Panel/full-span free-wake coupled method for unsteady aerodynamics of helicopter rotor blade

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.     

Study on the Behavior of Solar Array Deployment with Root Hinge Drive Assembly

In this paper, a method of using a root hinge drive assembly (RHDA) to control the solar array deployment is provided and a multi-DOF mechanism dynamic model of the system is established. In this way, the root hinge torque can be calculated iteratively. Then taking the predicted torque as a reference, a RHDA is designed for a large multiple-stage packaging and deployable solar array system. The control effect of the drive assembly is validated by ground tests. The test results indicate that the solar arrays can be deployed smoothly, and the deployment velocities are restricted by the drive assembly as expected. During the tests, the RHDA output speed and output torque are obtained. In order to examine the impact force when the yoke is lock-up with a hard stop, dynamics simulations are performed according to the actual behavior. The simulation result indicates that the designed RHDA reduces the impact force significantly and improves the lock-up reliability effectively.     

Modeling of a space flexible probe–cone docking system based on the Kane method

Recent developments in micro- and nano-satellites have attracted the interest of the research community worldwide. Many colleges and corporations have launched their satellites in space. Meanwhile, the space flexible probe–cone docking system for micro- and nano-satellites has become an attractive topic. In this paper, a dynamic model of a space flexible probe–cone docking system, in which the flexible beam technology is applied, is built based on the Kane method. The curves of impact force versus time are obtained by the Lagrange model, the Kane model, and the experimental method. The Lagrange model was presented in the reference and verified by both finite element simulation and experiment. The results of the three methods show good agreements on the condition that the beam flexibility and the initial relative velocity change. It is worth mentioning that the introduction of vectorial mechanics and analytical mechanics in the Kane method leads to a large reduction of differential operations and makes the modeling process much easier than that of the Lagrange method. Moreover, the influences of the beam flexibility and the initial relative velocity are discussed. It is concluded that the initial relative velocity of space docking operation should be controlled to a certain value in order to protect the docking system. a 2014 Production and hosting by Elsevier Ltd. on behalf of CSAA BUAA.     

A fully automated framework for helicopter rotor blades design and analysis including aerodynamics,structure, and manufacturing

This study describes an integrated framework in which basic aerospace engineering aspects (performance, aerodynamics, and structure) and practical aspects (configuration visualiza-tion and manufacturing) are coupled and considered in one fully automated design optimization of rotor blades. A number of codes are developed to robustly perform estimation of helicopter config-uration from sizing, performance analysis, trim analysis, to rotor blades configuration representa-tion. These codes are then integrated with a two-dimensional airfoil analysis tool to fully design rotor blades configuration including rotor planform and airfoil shape for optimal aerodynamics in both hover and forward flights. A modular structure design methodology is developed for real-istic composite rotor blades with a sophisticated cross-sectional geometry. A D-spar cross-sectional structure is chosen as a baseline. The framework is able to analyze all realistic inner configurations including thicknesses of D-spar, skin, web, number and ply angles of layers of each composite part, and materials. A number of codes and commercial software (ANSYS, Gridgen, VABS, PreVABS, etc.) are implemented to automate the structural analysis from aerodynamic data processing to sec-tional properties and stress analysis. An integrated model for manufacturing cost estimation of composite rotor blades developed at the Aerodynamic Analysis and Design Laboratory (AADL), Aerospace Information Engineering Department, Konkuk University is integrated into the framework to provide a rapid and dynamic feedback to configuration design. The integration of three modules has constructed a framework where the size of a helicopter, aerodynamic performance analysis, structure analysis, and manufacturing cost estimation could be quickly investigated. All aspects of a rotor blade including planform, airfoil shape, and inner structure are considered in a multidisciplinary design optimization without an exception of critical configuration.