2017 Vol. 32, No. 7

Display Method:
Effect of deformation and breakup characteristic on supercooled large droplet icing process
2017, 32(7): 1537-1544. doi: 10.13224/j.cnki.jasp.2017.07.001
Abstract:
To reveal and master the effects of the droplets deformation and breakup under supercooled large droplet (SLD) conditions, a numerical solver was developed to simulate the ice accretion. The ice accretion on the NACA 0012 airfoil was predicted. For the SLD movement, the droplet deformation and breakup process was investigated using a dynamically varying drag model and Taylor analogy method. Using the above methods, the numerical simulation of ice accretion and icing effects at the deformation and breakup of SLD was performed. These numerical results were compared with the computational and experimental data in several past literatures. The results are in good agreement with experimental data, indicating that the models and methods are feasible and effective. And the comparative results show that accounting of deformation and breakup has made significant improvement to the simulation of the droplet trajectory and ice shape distribution for the SLD icing prediction.
Regularity of two-phase flow with vaporization produced by injecting water to the high temperature airflow for nozzle
2017, 32(7): 1545-1553. doi: 10.13224/j.cnki.jasp.2017.07.002
Abstract:
In order to research the gas-liquid two-phase flow problem considering vaporization phenomenon for the rocket engine, the vaporization mechanism of liquid water was studied. The vaporization model for water was established according to the saturation temperature of water at different pressures. A program to solve the vaporization phase change process for liquid water was written on the basis of vaporization model, and the vaporization phase change calculation results were added to the multi-phase flow field equations as vaporization component source term and heat source term, so the three-dimensional multiphase flow field numerical solution considering the liquid water vaporization phase change was realized. Combining with the classical example, the errors of the temperature calculating values and experimental values of the two-phase flow field for the feature points were controlled within 8.5%, verifying the reliability of the three-dimensional calculation of the program. The mass conversion rate curves of two-phase from the different distance of the rocket nozzle showed that the mass conversion rate of two-phase rose from the lowest value to the highest value and then decreased to the lowest value, the transformation mechanism of the liquid and gas phases was revealed. The method can provide direct reference for the vaporization and cooling of the water chamber in the launching stage of launch platform for launch vehicle.
Design and experiment of air/kerosene torch igniter
2017, 32(7): 1554-1560. doi: 10.13224/j.cnki.jasp.2017.07.003
Abstract:
An air/kerosene torch igniter was designed to solve the ignition problem of scramjet. The energy of igniter was 300kW, the mass flux of air and kerosene was 98.9g/s and 6.7g/s at design point respectively. Under different equivalence ratio (ER) conditions, the physical and chemical parameters of out-flow at the exit of igniter were computed by CHEMKIN 4.0. The working model of igniter was optimized as “oxygen-rich model”. Kerosene was injected from swirl injectors at one end of igniter. Air was injected through two loops, i.e. “first injector” and “second injector”. The mixture of air and kerosene was ignited by a normal spark used in the car engine. The experimental system was established. The pressure measuring and photographs showed that the igniter could work reliably within ER range of 0.3-1.3. The starting time of the igniter was about 0.9s and the length of flame was 30cm. A pulse combustion phenomenon was represented by the high frequency and low amplitude fluctuating of pressure in igniters combustor. Scramjet could be ignited by this igniter reliably in experiments.
Effect of pilot design in the TeLESSⅡ low emission combustor on emission
2017, 32(7): 1561-1568. doi: 10.13224/j.cnki.jasp.2017.07.004
Abstract:
Commercial aircraft engine needs to satisfy International Civil Aviation Organization (ICAO) emission standards, the low emission combustor research and development is the key to make the engine development successful. TeLESSⅡ low emission combustor scheme was used to research the gas emissions property of two different pilot stage designs, namely, prefilming airblast atomizer and airblast simplex atomizer, with the same main stage air flow and fuel atomizing module. The results indicate that in low condition with different pilot stage designs,the emissions and combustion efficiency changes are different. In the high condition, according to fuel staging ratio, the emissions have different change trends, and the combustion efficiency change trend is the same. Those illustrate that the combination of pilot stage and main stage design change the combustion properties.
Flow and heat transfer characteristics of interlayer liquid cooling for 3D-IC with two stacked layers
2017, 32(7): 1569-1576. doi: 10.13224/j.cnki.jasp.2017.07.005
Abstract:
With the vertical stacking of circuit layers, the power density of 3D-IC increases exponentially. Interlayer liquid cooling is a promising and scalable solution for high heat flux removal in 3D-IC (three-dimensional integration circuit). The flow and heat transfer performances of water through uniform and double-side heat flux 3D-IC of interlayer in-line micro-pin fins structure with heat transfer areas of 1cm2 and Reynolds number ranging from 150 to 900 were studied numerically. 3D-IC with in-line micro-pin fins for pitches of 200μm, diameter of 100μm and heights of 200μm was analyzed. Results show that interlayer liquid cooling of 3D-IC with in-line micro-pin fins has better heat exchange effect than that with rectangular micro-channels. For the Reynolds number of 770, the power reaches to 250W, which is equivalent to 8.3kW/cm3 volumetric heat. Compared with rectangular micro-channels,the average temperature and the maximal junction temperature of the heater surface with in-line micro-pin fins are only 46.34, 13.96K, down by 13.26, 21.34K, respectively.
Experiment of turbulent flow in rotating smooth channel using two-dimensional hotwire
2017, 32(7): 1577-1584. doi: 10.13224/j.cnki.jasp.2017.07.006
Abstract:
The averaged velocities and Reynolds stresses were measured at different streamwise positions of a rotating smooth channel with two-dimensional hotwire. The results revealed that: the relatively high local rotation parameter expanded the influence region of rotation on the averaged velocity and induced the relaminarization of the turbulent flow near the leading side. The normalized averaged velocity profiles near the trailing side were strictly arranged in the rotation number order, and presented a logarithmic distribution in the log-law region. Meanwhile, all normalized Reynolds stress components were immune to rotation near the trailing side. Owing to the relaminarization, the normalized averaged velocity profiles near the leading side cannot maintain the logarithmic distribution in the log-law region, and all Reynolds stress components decayed with the growing rotation number and radius of the streamwise position. With the u-v quadrant analysis, the immediate reason of the relaminarization was attributed to the reduction of the generation of turbulent fluctuation.
Experiment on flow boiling heat transfer of R134a in mini-channel of flat tube
2017, 32(7): 1585-1591. doi: 10.13224/j.cnki.jasp.2017.07.007
Abstract:
An experimental study about R134a refrigerant flow boiling in a horizontal multiport flat tube mini-channel with the hydraulic diameter of 0.63mm and 0.72mm was carried out. The refrigerant was heated in uniform heat flux. The experiments were performed at the mass flux in range of 82-621kg/(m2·s), heat flux in range of 9.7-64kW/m2, saturation pressure in range of 0.22-0.63MPa and vapor quality in range of 0-1. Results showed that when the vapor quality in range of 0-0.6, the local heat transfer coefficients of mini-channel were higher than those of conventional channel, nucleate boiling was the main contributing factor, the local heat transfer coefficients increased with increasing heat flux and saturation pressure, but had little dependence with mass flux. When the vapor quality was bigger than 0.6, the local transfer coefficients decreased sharply; the influences of heat flux and saturation pressure on local heat transfer coefficient were small, but closely related to mass flux. Gungor-Winterton correlation was modified based on the experimental R134a flow boiling heat transfer coefficient from present study and published literatures, with -1.71% of the mean relation deviation and 19.24% of the mean absolute relation deviation respectively, showing an obvious improvement.
Performance evaluation methods of two-stage axial swirler:Ⅰ influence of total swirling intensity
2017, 32(7): 1592-1598. doi: 10.13224/j.cnki.jasp.2017.07.008
Abstract:
Theoretical and numerical studies on the reaction flow field and combustion performance of two-stage axial swirler combustors were carried out. A performance evaluation method of multi-stage swirlers was developed. In addition, the total swirl intensity and energy utilization efficiency were put forward to evaluate the swirler performance. It was found that the total swirl intensity and flow resistance coefficient were the main influential factors to the swirler energy utilization efficiency. The interaction between the two stage swirlers was investigated numerically. And the results indicate that the inner swirl intensity has a greater impact on the width of central toroidal recirculation zone (CTRZ). The total swirl intensity plays a most important effect on combustor overall performance. When the total swirl intensity is less than 0.43, weak swirling flows will be produced; when the total swirl intensity is between 0.43 and 0.6, medium strength swirling flows will be produced and weak CTRZs are presented; when the total swirl intensity is more than 0.6, strong swirling flows will be produced and moderate CTRZs are presented at the same time; when the total swirl intensity is more than 1.03, very strong swirling flows will be generated. If the total swirl intensity is constant, a two-stage axial swirler will get better energy utilization efficiency than a single stage axial swirler. The results by the evaluation method are compared with the experimental and numerical results. Its found that the evaluation method for multi-stage axial swirler performance can accurately evaluate the performance of swirlers, providing a practical and effective method for the design of multi-stage swirlers in the future.
Measurements on laminar burning velocities of hydrogen/oxygen/diluents at elevated pressure and temperature
2017, 32(7): 1599-1604. doi: 10.13224/j.cnki.jasp.2017.07.009
Abstract:
The effects of different diluents and dilution ratios on laminar burning velocities of hydrogen/oxygen/diluents mixtures were studied at elevated pressures and temperatures with a constant volume combustion chamber and a high-speed schlieren system. Both experimental and numerical laminar burning velocities of hydrogen/oxygen/diluents mixture were obtained at equivalence ratios of 0.6-4.0, initial pressures of 0.1-0.5 MPa, and dilution ratios of 4-8. The results showed that the laminar burning velocity varied non-monotonically with the increase of the initial pressures and equivalence ratios among 1.0-3.0,and the laminar burning velocity deceased with the increase of initial pressure when the equivalence ratio was less than 1.0 or more than 3.0. Through sensitivity analysis of elemental reaction, three-body reaction(R15) was found to restrain reaction in nitrogen diluted mixture, and promote reaction in both argon and helium diluted mixtures, mainly due to different thermal physics properties of different diluents.When the equivalence ratios were within the ranges of 1.0-3.0, the non-monotonic change of laminar burning velocities with the increase of initial pressures was mainly due to the competition between chain-branching reaction R1 and chain-terminating reaction R15.
Effects of wall temperature on pseudo-shock oscillations in isolator
2017, 32(7): 1605-1612. doi: 10.13224/j.cnki.jasp.2017.07.010
Abstract:
To explore the complex flow mechanism in the scramjet isolator under different wall heat transfer, the effects of wall temperatures on the oscillation characteristics of pseudo-shock in the scramjet isolator were investigated with a two-order accuracy unsteady computational fluid dynamics code. The numerical results show that there is a close relationship between the wall pressure histories and the maximum upstream and downstream locations of the shock-train leading-edge. Moreover, it shows that the fluctuating pressure amplitudes increase and frequencies decrease with the wall temperature. The maximum value of wall pressure standard deviation occurs at the place where the shock-train leading-edge shocks/boundary layer interact. However, the maximum values of kinetic energy are mainly in the shear layer region of the separation bubble.
Numerical investigation on restart characteristics of two-dimensional twin-duct supersonic inlet at sideslip state
2017, 32(7): 1613-1622. doi: 10.13224/j.cnki.jasp.2017.07.011
Abstract:
An unstarted problem of the inlet would be caused by a high combustors pressure of the ramjet. In order to study the restart characteristics of the inlet at a small sideslip angle, an investigation on a two-dimensional, twin-duct mixed compression supersonic inlet was conducted with a numerical method. Obvious hysteresis and asymmetry was found between the inlets flow pattern in the restarted duration. Results indicate that, the windward side inlet is easier to unstart and more difficult to restart than the leeward side inlet. Moreover, by comparing the restart characteristics of the inlet at different backpressures, it is noted that when the backpressure is higher than the inlets unstarted backpressure, the inlet is restarted at the same point even if the backpressure peak changes. However, there is reversed flow in the inlet when the backpressure is too much high, and a large reversed impulse is generated. The interior drag of the inlet rises abruptly, leading to deterioration of the engine performance.
Coupled simulation model of aero-engine performance and secondary air system
2017, 32(7): 1623-1630. doi: 10.13224/j.cnki.jasp.2017.07.012
Abstract:
An approach for zero-dimensional simulation of whole engine gas path was presented. By adopting modular modeling ideology, the transient engine performance model and air-system model were integrated and solved by means of dynamic data exchange. The corresponding modular simulation program was coded by utilizing object-oriented technology. A turbofan core engine was modeled to predict and analyze the dynamic changes of the air system bleeding and returning branches as well as their influences on the overall engine performance. Results show that the dynamic interaction between the main gas path and secondary air system exerts little influence on the overall performance of the aero-engine, but this influence is strong enough to cause significant imbalanced local responses of the secondary air system. This adverse impact should not be ignored in the design of modern high performance aero engine.
Analysis of civil aircraft aerodynamic characteristics with powered effects based on CFD/CSD method
2017, 32(7): 1631-1638. doi: 10.13224/j.cnki.jasp.2017.07.013
Abstract:
Fluid structure coupling method based on three dimensional Reynolds average Navier-Stockes (RANS) equations and static equation was used for study on aerodynamic characteristics of the civil aircraft with powered effects. Firstly, numerical calculation was adopted to simulate the boundary condition of the fan-inlet and jet-exhaust and analyze the flow field of the turbine powered simulator. The numerical results agreed well with those from the wind tunnel experiment. Furthermore, the research on the performance of the aerodynamic characteristics of the aircraft with powered nacelle and with flow-through nacelle was conducted through computational fluid dynamics(CFD)/computational structure dynamics (CSD) approach. The numerical results indicate that the aerodynamic performances such as lift coefficient, drag coefficient, aerodynamic loads and pressure center location based on powered nacelles differ obviously from flow-through ones. Elastic effect can increase the differences so that it can make lift-drag ratio decrease by approximately 12.6%, lift coefficient decrease by 8.9%, and pressure center location move backwards. The numerical example showed that the aerodynamic loads close to the nacelles were influenced notably by powered effects while those far from the nacelles were affected mainly by the elastic effects. Consequently, it is necessary to consider the elastic effects when analyzing civil aircraft aerodynamic characteristics.
An analytical method for longitudinal-lateral couplings in helicopter motional modes
2017, 32(7): 1639-1647. doi: 10.13224/j.cnki.jasp.2017.07.014
Abstract:
A quantitative method to analyze helicopter longitudinal-lateral coupling characteristics was put forward, including amplitude characteristics, time characteristics, and integration characteristics of couplings among state variables of a motional mode. Relative strengths of couplings could be reflected by amplitude characteristics. Time ratios of couplings with intensifying or diminishing effects could be described by time characteristics. And divergence or convergence rates of state variables caused by couplings could be represented by integration characteristics. With the typical longitudinal oscillatory mode of a hovering Puma helicopter as an example, the couplings characteristics among state variables of the mode were analyzed through comparison between longitudinal-lateral decoupled and coupled conditions. It showed that, because of longitudinal-lateral couplings, both lateral velocity and roll rate had effects on pitch rate, and their strength was almost three times the effect of longitudinal velocity on pitch rate, and the divergence rate of pitch rate was accelerated. Longitudinal-lateral couplings had a great effect on coupling characteristics of the mode, and the dynamic stability was influenced further.
Modified verison design method of axisymmetric unmixed-flow nozzle to mixed-flow nozzle
2017, 32(7): 1648-1657. doi: 10.13224/j.cnki.jasp.2017.07.015
Abstract:
In order to achieve the increased gross thrust and different spatial distributions of infrared radiation (IR) characteristics, a modified version design method of axisymmetric unmixed-flow nozzle based on CFD numerical simulation was established to change it into mixed-flow nozzle. Results indicated that: (1) An estimating method of nozzle thrust was proposed when the core and fan flow were mixed at various degrees, the error between the estimated value and the true one was basically less than 0.005 in the parameter scope. (2) Influence on modified version design of the core and fan flow condition was analyzed, increased thrust may be achieved when the total-pressure ratio of fan flow to core flow was between 0.8 and 1.5, and the mixing degree should be larger than 0.1. (3) Increased gross thrust was achieved after the modified version design, the thrust coefficient of mixed-flow nozzles with confluent mixer increased by 0.0087-0.0126 in comparison with the unmixed-flow nozzle; when lobed mixer was adopted, the thrust coefficient increased by 0.0289. (4) The integral IR intensity of mixed-flow nozzle increased in the azimuth angles of 0°-15°, but it significantly decreased in all the other azimuth angles, the largest decline was 77.5%.
Effect of elastic restriction stiffness of blade root on flying qualities of unmanned helicopter
2017, 32(7): 1658-1665. doi: 10.13224/j.cnki.jasp.2017.07.016
Abstract:
The flying qualities of the seesaw ducted fan configurational unmanned helicopter with an elastic support on the hub was evaluated for the stability, maneuverability and inter-axis coupling based on the U.S. armys aviation design standard ADS-33E-PRF. The flight dynamic model of the sample unmanned helicopter was introduced. The stability, maneuverability and inter-axis coupling in ADS-33E-PRF were pruned applicably based on the flight control characteristics of unmanned helicopter. By taking the sample helicopter with elastically supported blade roots as an example, the characteristics in trimming and flying qualities parameters of different stiffness of elastic restriction were computed, and the flying qualities of the sample unmanned helicopter were evaluated, compared and analyzed based on the flying quality evaluation methods. The comparative and analytical results show that the stability, maneuverability and inter-axis coupling of flying qualities are all influenced by the stiffness of elastic restriction, and the stiffness coefficient about 1 is more appropriate.
Friction torque of double-row spherical roller bearings
2017, 32(7): 1666-1675. doi: 10.13224/j.cnki.jasp.2017.07.017
Abstract:
Based on the dynamic theory of rolling bearings, the dynamic differential equations and the friction torque mathematic models of double-row spherical roller bearings were established. GSTIFF (gear stiff) integer algorithm with variable steps was used to solve the equations. The influences on friction torque of double-row spherical roller bearings caused by working conditions and structural parameters were investigated. The results show that: small inner ring osculation is useful to reduce the friction torque. Increasing outer ring osculation results in reduction of the friction torque. Too large or too small pocket curvature radius is not conducive to reducing the friction torque. There is a reasonable pocket curvature radius to get the lowest friction torque. The effect of cage pocket clearance on the bearing friction torque is small. A proper bigger cage guiding clearance is beneficial to reducing the friction torque. Properly increasing the radial clearance of double-row spherical roller bearings is useful to reduce the friction torque.
General thermal analysis of helicopter gearbox based on PCA-TN
2017, 32(7): 1676-1682. doi: 10.13224/j.cnki.jasp.2017.07.018
Abstract:
Targeting the problem of current thermal analysis methods which could not meet the complexity and polytropy of structure of helicopter gearbox, a general thermal analysis method which could quickly obtain the temperature distribution of helicopter gearbox was proposed. The combination of thermal network (TN) and principal component analysis (PCA) taking typical helicopter gearboxes as sample drew a conclusion of three models of general thermal analysis unit-body. Any thermal analysis model of helicopter gearbox of different structures could be quickly and efficiently established by unit-bodies and its derivatives. As an example, the method was applied to a helicopter gearbox. Through comparison of the analytical results and the experiment results, it was found that the temperature distribution of the helicopter gearbox could be high-efficiently obtained by the method, and the maximum error between analysis and experiment was 5.07%, meeting the requirement of engineering calculation for thermal analysis of gearbox.
VMD based adaptive composite multiscale fuzzy entropy and its application to fault diagnosis of rolling bearing
2017, 32(7): 1683-1689. doi: 10.13224/j.cnki.jasp.2017.07.019
Abstract:
A rolling bearing fault diagnosis approach was proposed based on the adaptive multiscale fuzzy entropy, ILS (iterative Laplacian score) and PSO-SVM (particle swarm algorithm optimization support vector machine). In the proposed method, the variational mode decomposition was used for the decomposition and reconstruction. Then composite multiscale entropy fuzzy of the reconstructed signals was calculated. Besides, the iterative Laplacian score algorithm was used for the sensitive fault feature selection and the selected features were input to the PSO-SVM based classifier for training and recognition. Finally, the proposed method was applied to experiment data of rolling bearing. Results showed that the identifying rate of proposed method was 100%. Also, the ILS based feature selection was compared with the SFS (sequential forward selection) method; and the result indicated that the highest identifying fault rate of SFS based method was 92.86% while the identifying fault rate of the ILS based fault diagnosis method reached to 100%.
Effect of cutting parameters on average scallop height for tooth surface of spiral bevel gears
2017, 32(7): 1690-1697. doi: 10.13224/j.cnki.jasp.2017.07.020
Abstract:
Based on the spiral bevel gear finish milling theories conducted by the numerical control(NC) face milling method, the expressions of motions for the NC axes were deduced, the influencing patterns of cutting speed, rate of movement of linkage axes, and degree of single step interpolation on the average scallop heights were studied, and their coupling relations were setup. Accordingly, the cutting experiment was conducted to verify the above results. The proposed study indicates that increasing the cutting speed and the rate of movement of the linkage axes has limited effect on average scallop height, while decreasing the degree of single step interpolation is able to reduce the average scallop height dramatically, which agrees with the theoretical analysis.
Fault diagnosis of rolling bearing’ compound faults based on improved time-frequency spectrum analysis method
2017, 32(7): 1698-1703. doi: 10.13224/j.cnki.jasp.2017.07.021
Abstract:
Based on the theory of cyclostationarity and two orders cyclic statistic for the spectrum correlation (SC) or spectrum correlation density (SCD), a time-frequency analysis method, namely improved spectrum correlation (ISC), was proposed and used in fault feature extraction of rolling bearing compound fault. Results show that the proposed method can extract the modulated frequency only and has more intuitive advantage than the traditional envelope demodulation spectrum method because the latter extracts the modulated frequency and carrier frequency simultaneously. The feasibility and effectiveness of proposed method are verified through simulation and experiment.
Plain bearing friction state recognition without complete prior knowledge
2017, 32(7): 1704-1711. doi: 10.13224/j.cnki.jasp.2017.07.022
Abstract:
Given that the prior knowledge of all kinds of bearing friction degradation model cannot be attained usually, starting from the similarities of different states, a bearing friction faults state recognition algorithm without prior knowledge was proposed based on sparse representation and absolute grey relational degree of B-mode (AGRDB). First, for the defects of sparse representation without supervision, the AGRDB was involved in the sparse representation, to get normal and severe friction codes under the largest distance between classes and smallest distance within the classes. Second, sparse representation model with discriminant sex was established under the same dictionary. And current state of the bearing was identified by comparing sparse coding and reconstruction error of normal lubrication and serious friction. Finally, the results of simulation signal and diesel engine bearing experiment show that the proposed method can better identify the sliding bearing early friction state (100-216min) and serious friction state (216-384min) under the less prior knowledge. And this algorithm is suitable for plain bearing fault monitoring under less samples.
Simplification method of complex air cooling bladed disk structure’s deformation analysis model
2017, 32(7): 1712-1717. doi: 10.13224/j.cnki.jasp.2017.07.023
Abstract:
Tip deformation analysis of high pressure (HP) turbine rotor based on the HP rotor can improve the tip clearance numerical simulation precision. Due to the complex air cooling structure of HP turbine rotor blade, the number of finite element mesh is huge. The mortise structure analysis of HP turbine blade and disk belongs to nonlinear analysis, which also needs enough calculation time. A simplified method of complex air cooling blade model and mortise structure contact analysis model was presented, the calculation efficiency was improved under the premise of without affecting the calculation accuracy. Applying this method to an aero-engine turbine components deformation analysis, the results were compared with the deformation analysis results of the original complex blade model. The relative error of the largest turbine blade tip radial deformation was 0.47%, and the calculation time reduced 99%. Its proven that the simplified method and the calculation method are feasible.
Fan rotor blade flutter free design technology
2017, 32(7): 1718-1727. doi: 10.13224/j.cnki.jasp.2017.07.024
Abstract:
Flutter free design technology of a two stages fan was investigated. During fan rotor blades design, appropriate aspect ratio was selected by considering aerodynamic performance, weight and flutter stability. Energy method was used to analyze the flutter stability for both rotor blades of the fan. Blade thickness, three dimensional blade shape and tip to hub chord ratio were modified to eliminate the flutter risk. The investigation showed that the aspect ratio was not the only parameter of concern for flutter free design in conceptual design phase. The near stall points at the rotational speed when blade tip flow was transonic were more prone to flutter. Strong shock waves could enhance the coupling effect between fluid and structure, and induce complex construction of aerodynamic work distribution. Blade thickness and tip to hub chord ratio were effective parameters to improve flutter risk.
Quantitative evaluation for the structure and mechanics properties of the high bypass ratio turbofan engines
2017, 32(7): 1728-1735. doi: 10.13224/j.cnki.jasp.2017.07.025
Abstract:
The structural efficiency parameters for the rotor system, supporting system and whole engine were put forward respectively according to the structure design requirements of the high bypass ratio turbofan engine, and the relations between the structural design parameters and the mechanics property parameters were established. The evaluation results for one typical engine exhibited that the average stress coefficients for the low pressure rotor of the high bypass ratio turbofan engine were between 0.2 and 0.3, lower than those of the other types of aero gas turbo engines because of the particularly large size of the fan. There were several unavoidable critical speeds with bending modes within the operation speed range, therefore the location of the joint structures of the rotor should be designed into the low stain energy zone. During maneuvering, the interval for the clearance between the stator and the rotor was [-1.4,1.0]mm for the whole engine, and the low pressure turbine was the key position for the clearance control.
Analysis on effect of mount stiffness on whole engine coupling vibration
2017, 32(7): 1736-1746. doi: 10.13224/j.cnki.jasp.2017.07.026
Abstract:
Using an aero-engine rotor tester with casing, a finite element model of engine to simulate the actual mounting condition was built. The modal tests of the tester under installation were carried out, and the finite element model was updated and verified based on the modal test results. The first 3 order natural frequencies and modal shapes of the tester were calculated under the conditions of the free boundary, fixed boundary and flexible mounting boundary with different mount stiffness values. By defining a rotor-stator coupling factor, the effect of the mount stiffness on the coupling degree between the stator and the rotor was studied. Results show that the coupling degree between stator and rotor increases with the increasing effect of mounting condition for a certain order mode. Moreover, the effect of the mount stiffness values on the rotor-stator coupling is nonlinear. Because the coupling phenomena between stator and rotor exist in many modes of practical large scale turbofan engines, the effect of the mount stiffness values on the coupling vibration of engine cannot be ignored. Therefore, the mount conditions should be considered carefully in finite element modeling and simulating.
Modeling method and reduction of dual-rotor system with complicated structures
2017, 32(7): 1747-1753. doi: 10.13224/j.cnki.jasp.2017.07.027
Abstract:
A 3-D finite element model was established for the dual-rotor model including complicated geometries. By utilizing Craig-Bampton mode synthesis method, the model order reduction was carried out. Through comparative analyses in critical speed and eigen problems, the validity and efficiency of the reduced order model were verified. On this basis, the critical speed map and the value of critical speed under a constant ratio were calculated. Results showed that the high-dimensional solid finite model was successfully reduced to a lower one with only a few dozen degree of freedom (34 DOFs) by means of Craig-Bampton method, and the maximum error of the first two orders critical speed and the first six orders eigen value were kept less than 0.43%. Based on analysis of the framework in which the model is created by reducing the two single rotor subsystems and combining them under compatibility condition, dynamical problems such as critical speed can be easily solved.
Robust design method of rabbet joint structure in high speed assemble rotor
2017, 32(7): 1754-1761. doi: 10.13224/j.cnki.jasp.2017.07.028
Abstract:
A kind of assembly rotor connected by rabbet joint structure was introduced. Based on analyzing the stiffness loss model of rabbet joint structure, the effect of joint structure on the overall stiffness of the rotor was obtained. And then through the nonlinear finite element contact analysis, the rigidity robust design method was obtained. Results show that,by means of increasing the local stiffness of rabbet, avoid bending strain energy concentrated in the joint structure, and ensure 20% margin of the critical speed of bending from the working speed, the rigidity robust can be enhanced; by means of designing the rabbet with interference fit, improve the axial compression force and optimize the local structure of the rabbet, the contact state robust can be enhanced.
Aerodynamic stability optimization of controlled diffusion airfoil
2017, 32(7): 1762-1768. doi: 10.13224/j.cnki.jasp.2017.07.029
Abstract:
Direct method was applied with an automatically optimization tool, in which the geometrical code and the blade to blade flow solver were integrated into genetic algorithm for searching of the optimal airfoil based on a conventional controlled diffusion airfoil(CDA). The design object was considered as a compromise of losses and operating range. It showed that the attainable incidence of optimized airfoil increased from 11° to 17.5° and the negative incidence range increased 4.5°, with similar design point loss as reference CDA airfoil. Additionally, the loss curve of optimized airfoil was flatter, showing it can remain stable performance in a wider working range. It indicated that the upstream propagation of suction peak shifted from 40% chord to 20% chord, extended the diffusion area and reduced the deceleration gradient, contributing to the increase of stall margin. The increase of negative incidence range can be ascribed to two reasons: first, the optimized airfoil has larger throat area and the throat location is apart from suction peak, thus gaining larger choking range. Secondly, the pressure surface velocity distribution is flatter, even with slight acceleration until 30% chord, as a result, the early separation of boundary layer caused by the leading edge spike and the following continual diffusion can be avoided.
Design exploration of a highly loading intermediate transition duct of turbine
2017, 32(7): 1769-1781. doi: 10.13224/j.cnki.jasp.2017.07.030
Abstract:
Aerodynamic design and numerical simulation were carried out to investigate a intermediate transition duct (ITD). Results of the analyses and comparisons showed that the thickness distribution of ITD strut vane(SV) had little impact on the main flow when SV number was small,while the pressure gradient and separation of SV surface were affected mainly through the profile curvature distribution. The flow separation in the region of hub, shroud and SVs tip was controlled effectively with the combination of concave and convex curvature. The stagger angles of the SV on S1 stream surface were changed by circumferential lean to impact the effect of attack angle and depress the flow separation. The SV should be placed in the main flow with low pressure gradient as more as possible no matter part or all to decrease the pressure gradient and risk of flow separation.
Control of hub-corner stall in a compressor with stator hub gap
2017, 32(7): 1782-1792. doi: 10.13224/j.cnki.jasp.2017.07.031
Abstract:
To research the control of hub-corner stall in a compressor with stator hub gap, a 1.5-stage axial flow compressor was studied. The effects of integral hub gap and partial hub gap on the aerodynamic performance of compressor at low-operating point and design point were studied by three-dimensional numerical simulation method. The results show that the integral hub gap can restrain the hub-corner stall and improve the performance of compressor, because the energy source of the “tornado” vortex is weakened by the leakage flow, and the “tornado” vortex begins from the hub surface and ends at the suction surface, but the performance of compressor is reduced due to the leakage loss produced by stator hub gap at design point. And the partial hub gap is significantly better than the integral hub gap, with the position of partial hub gap moving backward, the performance of compressor is gradually improved at low-operating point, and the damage to the performance at design point is gradually reduced. For example, the mass flow rate of TAI2 scheme increased by 0.89kg/s and the efficiency increased by 1.25% at low-operating point, and the efficiency of design point remained unchanged. In addition, only when the partial hub gap increases to a certain size could the hub-corner stall be restrained by enough leakage flow.