• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Effects of axial gap and nozzle distribution on aerodynamic forces of a supersonic partial-admission turbine

    2017-12-22 06:23:06JinpengJIANGJiwenLIGuobioCAIJueWANGQiynLIN
    CHINESE JOURNAL OF AERONAUTICS 2017年6期

    Jinpeng JIANG,Jiwen LI,Guobio CAI,Jue WANG,Qiyn LIN

    aSchool of Astronautics,Beihang University,Beijing 100083,China

    bSystem Engineering Division,China Academy of Launch Vehicle Technology,Beijing 100076,China

    cBeijing Aerospace Propulsion Institute,Beijing 100076,China

    Effects of axial gap and nozzle distribution on aerodynamic forces of a supersonic partial-admission turbine

    Jinpeng JIANGa,Jiawen LIa,Guobiao CAIa,Jue WANGb,*,Qiyan LINc

    aSchool of Astronautics,Beihang University,Beijing 100083,China

    bSystem Engineering Division,China Academy of Launch Vehicle Technology,Beijing 100076,China

    cBeijing Aerospace Propulsion Institute,Beijing 100076,China

    Aerodynamic force; Axial gap; Computational fluid dynamics(CFD); Nozzle distribution; Partial admission; Turbine

    The turbine in an LH2/LOX rocket engine is designed as a two-stage supersonic partialadmission turbine.Three-dimensional steady and unsteady simulations were conducted to analyze turbine performance and aerodynamic forces on rotor blades.Different con figurations were employed to investigate the effects of the axial gap and nozzle distribution on the predicted performance and aerodynamic forces.Rotor blades experience unsteady aerodynamic forces because of the partial admission.Aerodynamic forces show periodicity in the admission region,and are close to zero after leaving the admission region.The unsteady forces in frequency domain indicate that components exist in a wide frequency region,and the admission passing frequency is dominant.Those multiples of the rotational frequency which are multiples of the nozzle number in a fulladmission turbine are notable components.Results show that the turbine efficiency decreases as the axial gap between nozzles and the 1st stage rotor(rotor 1)increases.Fluctuation of the circumferential aerodynamic force on rotor 1 blades decreases with the axial gap increasing.The turbine efficiency decreases as the circumferential spacing between nozzles increases.Fluctuations of the circumferential and axial aerodynamic forces increase as the circumferential spacing increases.As for the non-equidistant nozzle distribution,it produces similar turbine performance and amplitude frequency characteristics of forces to those of the normal con figuration,when the mean spacing is equal to that of the normal case.

    1.Introduction

    Supersonic turbines with partial admission are sometimes employed in rocket engines to avoid an extremely low blade height which would result in large secondary losses such as the tip leakage loss.In addition,partial admission con figurations are commonly applied to control the power output and get higher efficiency in the control stages of power plants and industrial steam turbines.However,additional forms of loss exist in partial admission compared with full-admission turbines,such as the mixing loss generated in the interface regions between the low-energy dead flow and the high-energy through flow,which significantly influences the efficiency.1Furthermore,circumferential non-uniformity is increased due to partial admission.The flow in a turbomachine is highly unsteady because of the interaction between adjacent rows due to the effects of wakes.2,3In a partial-admission turbine,rotor blades periodically pass through flowing regions as well as regions of no flow,inducing strong flow-exciting forces besides vane wakes,which may result in high cycle fatigue.4Rotor blade cracks once occurred in turbopump tests for some liquid rocket engines designed with a partial-admission turbine,and the failure is believed to be caused by exciting forces.It is certainly necessary to investigate the unsteady flow and aerodynamic forces in partial-admission turbines.

    Experimental investigations and theory analysis of partial admission have been conducted by several researchers.5–9Numerical simulation has also been a powerful approach widely used by researchers.Full 360°Computational Fluid Dynamics(CFD)models and 3-Dimensional(3D)transient simulations are usually required,and sometimes multistage simulation is needed,which all make numerical investigations with high computational costs.102-Dimensional(2D)1,11and quasi-3D12simulations were usually conducted in the past.With the development of CFD and the improvement of computers,3D transient numerical simulations for investigations of partial admission have been widely conducted.Effects of the axial gap between the first-stage stator and rotor and multiblocking on the performance of partial-admission turbines were numerically studied by Hushmandi and Fransson.13They illustrated that reducing the axial gap produced a better efficiency at the first stage,but the overall efficiency at the second stage was decreased.Numerical simulations were conducted for a steam turbine with different degrees of partial admission by Qu et al.14,and they pointed out that the unsteady computational result was more accurate to analyze the flow of the control stage.Newton et al.15investigated the sources of loss in a turbine with both full and partial admission numerically,and evaluated the distribution of loss with entropy production.From these previous investigations about partial admission,we can see that steam turbines has usually been studied and the turbine aerodynamic performance has mostly been focused on.However,unsteady forces on rotor blades for a rocket engine turbine have been investigated by relatively few studies,especially a turbine with partial admission.J?cker et al.16numerically studied the aerodynamic excitation mechanisms due to the unsteady stator-rotor interaction in a supersonic turbine for a rocket engine turbopump.Hudson et al.17experimentally measured and numerically simulated the surface pressure on a rocket engine turbine blade.Tokuyama et al.18conducted numerical simulation for a partial-admission turbine used for a rocket engine.The blade scaling procedure was adopted,which resulted in a simulation with one third of full passages.The unsteady aerodynamic force was analyzed.They showed that unsteady force components appeared in wide frequencies.Using a genetic algorithm,the partial admission of a supersonic turbine was optimized by Tog and Tousi.19

    The turbine in the current study is used for an expander cycle LH2/LOX rocket engine.It is a two-stage supersonic turbine with partial admission.A first-stage rotor blade cracked in a test.Besides processing defects,the excitation produced by the partial admission con figuration and supersonic nozzles was considered to result in the failure.Therefore,the unsteady flow in a full-annulus turbine was numerically investigated in this paper to study the unsteady aerodynamic forces.Furthermore,effects of the axial gap and nozzle distribution on the turbine performance and aerodynamic forces were studied,in order to search for a possible method to reduce the excitation.Con figurations with different equal circumferential spacings between nozzles and unequal circumferential spacings were analyzed to study the effect of the nozzle distribution on the turbine performance and unsteady forces on rotor blades.Con figurations with various axial gaps between nozzles and the firststage rotor were analyzed to investigate the in fluence of the axial gap on the turbine performance and aerodynamic forces.

    2.Computational procedures

    2.1.Fundamental equation

    Reynolds time-averaged Navier-Stokes equations are solved by commercial CFD software of ANSYS-CFX.The general form of governing equations for the compressible viscous unsteady flow in a Cartesian coordinate is expressed as follows:

    where ρ is density,φ is the general variable for different equations,t is time,U is velocity,Cφis the general diffusion coef ficient,and Sφis the general source term.Conservation equations for mass,momentum,and energy can be obtained by setting φ to 1,u,v,w,and h,where u,v,w are velocity components,and h is enthalpy.Turbulence models can be obtained when setting φ to turbulence kinetic energy k and turbulent frequency ω.The ideal gas equation of state is adopted for enclosure.

    The k-ω based Shear Stress Transport(SST)turbulence model was employed in this study,since it can give highly accurate predictions of the onset and the amount of flow separation under adverse pressure gradients.20–22The SST model combines the advantages of the k-ω model in the near-wall region and the k-ε model in the bulk domain,and accounts for the transport of the turbulent shear stress.The proper transport behavior can be obtained by a limit to the formulation of the eddy-viscosity as:

    where a1is equal to 5/9,S is an invariant measure of the strain rate,and F2is a blending function expressed as:

    where β′=0.09,y is the distance to the nearest wall,and ν is the kinematic viscosity.

    The governing equations were solved using the Finite Volume Method(FVM).The upwind schemes were employed to discrete advection terms,which gave results that agreed well with test results.The second-order implicit schemes were adopted for temporal discretization.

    2.2.Computational model and boundary condition

    The turbine studied in this paper is a two-stage supersonic turbine with partial admission,and consists of 4 nozzles,56 blades for the 1st stage rotor(rotor 1),28 vanes for the 2nd-stage stator(stator 2),and 58 blades for the 2nd-stage rotor(rotor 2).The rotational speed of rotors at the design point is 42431 r/min.Nozzle and blade pro files are shown in Fig.1,and the CAD(Computer Aided Design)model is shown in Fig.2.

    The 3D full-annulus model was taken in simulations.Structural grids for nozzles were generated by Gambit.Computational grids for the 1st rotor and 2nd stator and rotor passages were generated by the IGG/Autogrid software package.Structured grid systems were created in the computational domain with O-type grids near the blade surfaces and H-type grids in the other regions.A grid-dependency check was conducted by simulating the steady flow of the normal con figuration with different grid sizes.The grid size was determined to ensure that the total-to-total isentropic efficiency was stable above that size.The total number of grid nodes was approximately 4500000.

    Cases for different axial gaps include a small gap con figuration(20%of the 1st blade axial chord),a normal gap con figuration(32%of the 1st blade axial chord),and a large gap con figuration(46%of the 1st blade axial chord).Table 1 lists the axial gap sizes for the three cases,where bris the 1st blade axial chord length at midspan.Fig.3 shows the considered computational con figurations focused on the axial gap.The axial gap was changed by shifting a nozzle axially.The domain behind the nozzle was elongated or shortened correspondingly,while the rotor domain stayed the same.By doing so,the axial distances between the nozzle-rotor interface and the blade leading edge were the same for all computational con figurations.However,the distances between the nozzle outlet and the interface were different.

    Fig.4 is a sketch showing the nozzle distribution,in which nozzles are numbered from 1 to 4 along the direction of rotor rotating,and central angles α1, α2,and α3are used to present the circumferential spacings between nozzles.The circumferential distribution of turbine nozzles was changed by rotating nozzles around the axis.Four con figurations with differentuniform circumferential spacings and a case with non-uniform circumferential spacings were studied.The circumferential spacing sizes for different cases are summarized in Table 2.When the circumferential spacing is 27.5°,the nozzle outlet edges overlap.

    Table 1 Cases with different axial gaps.

    Table 2 Cases with different circumferential spacings.

    The working gas in the turbine is hydrogen-rich gas.The specific heat capacity at constant pressure is 8254 J/(kg·K),and the specific heat ratio is 1.3666.The total pressure 3.061 MPa and total temperature 834 K were specified at the inlet.The flow direction was normal to the boundary,and a 5%turbulent intensity was set.Static pressure was imposed at the outlet,which was 0.45 MPa.The mixing plane method was applied for interfaces between the stator and rotor in steady simulations,while the sliding mesh technology was used in unsteady simulations,and steady simulation results were used as initial values for unsteady simulations.No-slip and adiabatic conditions were enforced along all the solid surfaces.100 time steps were set in a nozzle-passing period which was de fined as the period of time required for one blade to pass one nozzle.The convergence criteria for the root mean square residual was set to 1×10-5.20 iterations were performed for each time step to ensure a convergence.

    3.Results and discussion

    3.1.Flow field analysis

    Steady and unsteady simulation results for the designed geometry are analyzed in this section.

    Simulation results of mass flow rate(m˙)and turbine efficiency(ηT)are compared with test results in Table 3,and the absolute deviation(Da)and relative deviation(Dr)are given.Test results were obtained from the firing tests for the LH2/LOX rocket engine.The measurement accuracy of mass flow was 0.5%.Predicted errors for mass flow rate and the turbine efficiency are both lower than 1.5%.It reveals that the simulation results are in good agreement with the test results.

    A contour plot of the Mach number at midspan is shown in Fig.5.A circumferential non-uniform distribution is produced due to the partial admission con figuration.Flow at the nozzle exit is supersonic,and the Mach number in most of the area is between 1 and 1.7.A rapid flow acceleration occurs around thenozzle exit corner due to the expansion of the working gas,and the maximum Mach number exceeds 2.The acceleration flow is hindered by the next nozzle’s exit flow which is affected in turn by the acceleration flow from the previous nozzle.The high-Mach number region is larger at the No.4 nozzle exit corner leaving the admission area since not being affected by other nozzle exit flow.

    Table 3 Comparison of turbine specific work and efficiency.

    The static pressure contours at midspan is shown in Fig.6.Low pressure occurs at the nozzle exit corner due to the expansion of the working gas.Fig.7 shows the circumferential pressure distributions at four positions between stages,and the pressure p is normalized by the inlet total pressure p0,in.The circumferential pressure downstream from rotor 2 is weak,which is close to the outlet pressure and hardly affected by upstream flow.The other pressure distributions have strong circumferential variations due to the partial admission.The pressure distributions show obvious periodicity in the admission area,corresponding to the four nozzles.Pressure is close to the outlet pressure in the non-admission region.Pressure around the region leaving the admission area seems to be lower than that in the region entering the admission area,because the supersonic gas expands and accelerates at the nozzle outlet corners.The pressure distribution between a nozzle and rotor 1 is the most complex.Many local pressure peaks are caused due to the interaction between the supersonic flow and rotor 1 blades.Strong pressure decreasing regions between a nozzle and rotor 1 occur between nozzles,because of the expansion of the supersonic gas at the nozzle exit corners.A larger area of low pressure appears behind the No.4 nozzle.As can be seen,four pressure peaks exist at positions between rotor 1 and stator 2 and between stator 2 and rotor 2,which result from the four nozzles.Pressure distributions at the admission area seem to be approximated by a sine wave.Slight fluctuations around the main wave can also be seen due to the influence of upstream wakes.

    Aerodynamic forces on a rotor blade are calculated as follows:

    where Ftand Fzare the circumferential and axial forces,respectively,Fviscorepresents the viscous force on the blade,ntand nzare the circumferential and axial unite vectors,respectively,andRspds represents the surface integral of static pressure on the rotor blade.

    The unsteady aerodynamic forces(F)on a rotor 1 blade in time domain are plotted in Fig.8.Time was normalized by the rotational period(T),and one cycle is shown.The aerodynamic forces are presented corresponding to the circumferential locations.Both circumferential and axial forces show characteristic due to partial admission.In the admission area,the unsteady forces show great fluctuations because of the nonuniform upstream flow from the nozzles,and periodicity can be seen from the figure.Force peaks occur between nozzle outlets and at the corner of the No.4 nozzle,where a large pressure difference between the pressure side and the suction side occurs.After leaving the admission area,forces on the rotor blade in both directions are close to zero.The circumferential forces are mostly positive(along the rotating direction).There are unusual fluctuations when the blade enters and leaves the admission area.Negative forces are produced around the region entering the admission area.This can be expected,because the suction side enters the high-pressure region ahead of the pressure.For a similar reason,a small positive force peak occurs when leaving the admission area.The axial forces are mostly negative(opposite to the axial flow direction).

    Fig.9 shows the results of Fast Fourier Transformation(FFT)analysis of the unsteady aerodynamic forces on the rotor 1 blade.Many frequency components exist in a wide region,which is a characteristic of a partial-admission turbine.The 1st harmonics of the rotational frequency means the admission passing frequency.The 11 multiples of the rotational frequency approximately correspond to the nozzle number in a full-admission turbine.As can be seen,the 1st harmonics is a dominating component in both directions.This can be identified in the regular circumferential distribution of the unsteady force due to the admission as shown in Fig.8.The 1st harmonics for the circumferential force is an absolutely dominating component,while it is relatively weak for the axial force in comparison with the circumferential force.The 11,10,and 12 multiples of the rotational frequency are remarkable,while the 22,21,and 23 multiples of the rotational frequency are notable.As we can see,multiples of 11 multiples components and components adjacent to them are outstanding.However,components over 30 multiples of the rotational frequency are weak.

    The unsteady aerodynamic forces on a rotor 2 blade in time domain are shown in Fig.10.Both circumferential and axial forces are also strongly affected by partial admission,just like those on a rotor 1 blade.In the admission area,the forces are mostly positive,and show remarkable periodicity.After the blade leaves the admission area,both forces are close to zero.There are unusual fluctuations for the circumferential force when the blade enters and leaves the admission area.Negative forces are produced around the region entering the admission area,and a small negative force peak occurs when leaving the admission area.

    Fig.11 shows the results of FFT analysis of the unsteady aerodynamic forces on the rotor 2 blade.Components over 15 multiples of the rotational frequency are negligible.The frequency range is much narrower compared to that of rotor 1.The 1st harmonic is dominant in both forces.The 11,10,and 12 multiples of the rotational frequency are notable,just like those on a rotor 1 blade.The 2nd harmonic is more significant than 11 multiples component compared to the force on a rotor 1 blade.

    Forces on a rotor 2 blade are much lower than forces on a rotor 1 blade,because the pressure in rotor 2 is much lower than that in rotor 1.Thus,aerodynamic forces on a rotor 1 blade are concerned in the following analysis.

    3.2.Effects of axial gap

    Steady simulations were conducted for cases with different axial gaps,and turbine performances are listed in Table 4.The differences of turbine efficiency from that of Case 0(ΔηT)were calculated,and were expressed as percentages of the normal gap case efficiency.As can be seen,the turbine efficiency decreases as the axial gap increases.This is due to the increase of wake mixing loss between the nozzle exit corners and leading edges of rotors as the axial gap increases.

    Fig.12 shows the contours of static entropy(s)on the nozzle-rotor interface plane for con figurations with different axial gaps.As the axial gap increases,the entropy between nozzles decreases.However,the main flow is affected and the entropy downstream from the nozzle outlet increases.In addition,entropy near the inner end wall in the admission region increases when the axial gap increases.As can also be seen in Fig.13 which shows the contours of total pressure(pt)on the nozzle-rotor interface plane,as the axial gap increases,the total pressure between nozzles increases,but the total pressure downstream from the nozzle outlet decreases.By comparing the contours in Figs.12 and 13,it can be seen that increasing the axial gap declined the wake strength and increased the uniformity of flow.However,increasing the axial gap increased the mixing of the wake and the main flow,and losses could be expected to increase.

    Histories of the circumferential force on a rotor 1 blade for cases with different axial gaps are plotted in Fig.14.It can be seen that as the axial gap increases,the force peak decreases and the valley increases in the admission area.The fluctuation range of the aerodynamic force apparently reduces,while the mean force barely changes.

    Results of FFT analysis of the unsteady circumferential force on a rotor 1 blade for cases with different axial gaps are compared in Fig.15.As can be seen,the amplitudes for main frequencies decrease as the axial gap increases.The amplitude for 1st harmonics slightly decreases,while those for 11,10,22,and 21 multiples of the rotational frequency significantly reduce.It indicates that high-order frequency components are more affected by the axial gap.

    Fig.16 shows histories of the axial force on a rotor 1 blade for con figurations with different axial gaps.Major force valleys show little difference among different con figurations.The con figuration with a small axial gap generates a relatively larger minor force valley in the admission area.Generally,the axial force is slightly influenced by the axial gap.

    Table 4 Turbine performances for different axial gaps.

    Results of FFT analysis of the unsteady axial force on a rotor 1 blade for cases with different axial gaps are plotted in Fig.17.When the axial gap increases,the interaction between nozzles and rotor blades weakens.Minor force valleys decrease and flatten just as shown in Fig.16,and thus highorder frequency components decrease.As can be seen in Fig.17,over the 12 multiples of the rotational frequency,amplitudes decrease when the axial gap increases.Considering that axial forces were slightly affected by the axial gap,when high-order frequency components strengthened,low-order frequency components relatively weakened.It shows that the amplitudes of main low-order frequency components(1st harmonic,11 and 10 multiples of the rotational frequency)increase as the axial gap increases.

    In general,increasing the axial gap helps to reduce the unsteady aerodynamic forces,while the turbine efficiency would be sacrificed to some extent.

    3.3.Effects of nozzle distribution

    Table 5 Turbine performances for different circumferential spacings.

    Steady simulations were conducted for con figurations with different nozzle distributions,and turbine performances are compared in Table 5.As can be seen,the turbine efficiency decreases as the circumferential spacing increases.It can be concluded that increasing the circumferential spacing tends to decrease the turbine efficiency,although the reduction is relatively small.This can result from the declining uniformity and the increasing mixing loss between nozzles as the circumferential spacing increases.The con figuration with a non-equidistant nozzle distribution has similar performance to that of the normal case.It’s probably because the former has the same average circumferential spacing as that of the latter.

    Fig.18 shows the contours of static entropy on the nozzlerotor interface plane for con figurations with different circumferential spacings.As the circumferential spacing increases,the entropy between nozzles increases and the width increases corresponding to the circumferential spacing.That indicates that both the wake strength and width increase when the circumferential spacing increases.In addition,entropy near the inner end wall in the admission region increases with the circumferential spacing increasing.

    Fig.19 shows the contours of total pressure on the nozzlerotor interface plane.As the circumferential spacing increases,the total pressure between nozzles decreases.When the circumferential spacing increases,it becomes difficult for the wake to recover to the main flow.By comparing the contours in Figs.18 and 19,it can be concluded that increasing the circumferential spacing increases the wake strength and decreases the uniformity of flow.

    Histories of the circumferential force on a rotor 1 blade for cases with different circumferential spacings are plotted in Fig.20.It can be seen that as the circumferential spacing increases,the force peak slightly increases and the valley significantly reduces in the admission area,except for the No.4 nozzle.The circumferential spacing has little influence on the force peak and valley downstream from the No.4 nozzle.As a result,the amplitude of the circumferential aerodynamic force in the admission region increases,while the mean force slightly decreases.Apparently,the decreasing uniformity between nozzle outlets increases the pressure difference,and thus results in circumferential force peaks increasing.However, flow at the corner of the No.4 nozzle is not affected by the circumferential spacing,and thus the 4th force peak changes little.As for the con figuration with a non-equidistant nozzle distribution,the force peak and valley are equal to those in cases with responding spacing.Downstream from nozzles leaving the admission area,the force peak and valley are almost the same for con figurations with different spacings.

    Fig.21 shows the results of FFT analysis of the unsteady circumferential force on a rotor 1 blade for cases with different circumferential spacings.The 1st harmonic is dominant in all cases.However,the main frequency components vary due to different circumferential spacings. Main frequency components for the con figuration with a 27.5°spacing are approximately multiples of the 13 multiples of the rotational frequency,those for the 30°spacing case are approximately multiples of 12 multiples of the rotational frequency,and those for the 35°spacing case are approximately multiples of 10 multiples of the rotational frequency.The con figuration with unequal spacings has similar amplitude-frequency characteristics to those of the normal case.The amplitude of the 1st harmonic slightly varies,while the amplitudes of main high-order frequency components apparently increase as the circumferential spacing increases.

    Histories of the axial aerodynamic force on a rotor 1 blade for con figurations with different circumferential spacings are shown in Fig.22.The force peaks in the admission area increase as the circumferential spacing rises,while the force valley value slightly decreases.Similarly,the force peak and valley downstream from the No.4 nozzle are hardly affected.It reveals that increasing the circumferential spacing could result in increasing of the fluctuation range of the axial aerodynamic force.

    Results of FFT analysis of the unsteady axial force on a rotor 1 blade for con figurations with different circumferential spacings are plotted in Fig.23.The amplitude-frequency characteristics of the axial force are influenced by the circumferential spacing in a way that the circumferential force is affected.The amplitude of the 1st harmonic decreases,while the amplitudes of main high-order frequency components increase as the circumferential spacing increases.

    Generally,increasing the circumferential spacing between nozzles could reduce the turbine efficiency and increase the fluctuation range of aerodynamic forces.It is better to use a circumferential spacing as small as possible when designing nozzles for a partial-admission turbine.

    4.Conclusions

    3D steady and unsteady flow simulations of a two-stage supersonic partial-admission turbine have been conducted for different con figurations,to investigate effects of the axial gap and circumferential spacing on the predicted aerodynamic performance and forces.Results show that:

    (1)A strong circumferential pressure distribution is created by the partial admission.Pressure distributions show obvious periodicity in the admission region,corresponding to the nozzle number.The pressure distribution between nozzles and rotor 1 is especially complex.

    (2)Rotor blades experience unsteady aerodynamic forces due to the partial admission.The circumferential force fluctuates when the blade enters and leaves the admission region.Aerodynamic forces show periodic variations due to nozzle exit flow in the admission area,and are close to zero after leaving the admission region.

    (3)Unsteady force components on a rotor 1 blade appear in a wide frequency region,and the admission passing frequency which is also the rotational frequency is dominant.Those multiples of the rotational frequency which are multiples of the nozzle number in a fulladmission turbine are notable components.

    (4)The turbine efficiency decreases as the axial gap between nozzles and rotor 1 increases.At the same time,the fluctuation range of the circumferential aerodynamic force in the admission region decreases.In frequency domain,amplitudes of main frequency components for the circumferential force on a rotor 1 blade decrease with the axial gap increasing.As for the axial force on the rotor 1 blade,the variation is relatively complex.Amplitudes of frequency components below the12multiples increase,while those over 12 multiples decrease as the axial gap increases.Generally,unsteady aerodynamic forces can be reduced by increasing the axial gap at some expense of turbine efficiency.

    (5)The turbine efficiency reduces as the circumferential spacing between nozzles increases.Fluctuation ranges of circumferential and axial aerodynamic forces increase.In frequency domain,the amplitude of the 1st harmonic decreases,while the amplitudes of highorder main frequency components increase for all forces as the circumferential spacing increases.In general,the force fluctuation in the admission region increases with the spacing rising.As for the non-equidistant nozzle distribution,it produces similar turbine performance and amplitude-frequency characteristics of forces to those of the normal con figuration,when the mean spacing is equal to that of the normal case.A circumferential spacing as small as possible is recommended for high turbine efficiency and small aerodynamic forces.

    1.He L.Computation of unsteady flow through steam turbine blade rows at partial admission.Proc Inst Mech Eng Part A:J Power Energy 1997;211(3):197–205.

    2.Zhao B,Yang C,Chen S,Qi MX,Zhou M.Unsteady flow variability driven by rotor-stator interaction at rotor exit.Chin J Aeronaut 2012;25(6):871–8.

    3.Ghenaiet A,Touil K.Characterization of component interactions in two-stage axial turbine.Chin J Aeronaut 2016;29(4):893–913.

    4.Xie YH,Gao KK,Lan JB,Xie GN.Computational fluid dynamics modeling three-dimensional unsteady turbulent flow and excitation force in partial admission air turbine.Math Probl Eng 2013;2013:251926.

    5.Yahya SM,Doyle MDC.Aerodynamic losses in partial admission turbines.Int J Mech Sci 1969;11(5):417–31.

    6.Cho CH,Cho SY,Ahn KY.A study of partial admission characteristics on a small-scale radial-in flow turbine.Proc Inst Mech Eng Part A:J Power Energy 2010;224:737–48.

    7.Beer W,Hirsch P.In fluence of partial admission on the downstream stages of a reaction turbine:An analytical approach verified by measurements.New York:ASME;2013.Report No.:GT2013-570 94661.

    8.Li JL,Guo XH,Yu C,Huang SH.Catastrophe performance analysis of steam- flow-excited vibration in the governing stage of large steam turbines with partial admission.J Energy Res Technol 2013;135(4):041601.

    9.Yahya SM.Some tests on partial admission turbine cascades.Int J Mech Sci 1969;11(10):853–66.

    10.Kalkkuhl TJ,Engelmann D,Harbecke U,Mailach R.Numerical analysis of partial admission flow in an industrial steam turbine.New York:ASME;2012.Report No.:GT2012-68482.

    11.Lampart P,Szymaniak M,Rzazdkowski R.Unsteady load of partial admission control stage rotor of a large power steam turbine.New York:ASME;2004.Report No.:GT2004-53886.

    12.Sakai N,Harada T,Imai Y.Numerical study of partial admission stages in steam turbine.JSME Int J Ser B 2006;49(2):212–7.

    13.Hushmandi NB,Fransson TH.Effects of multiblocking and axial gap distance on performance of partial admission turbines:a numerical analysis.J Turbomach 2011;133(3):031028.

    14.Qu HC,Zhang D,Xie YH,Cao SH,Wu QL.Numerical simulation of unsteady flow phenomena and exciting force in control stage of steam turbine under multiple working conditions.New 591 York:ASME;2011.Report No.:POWER2011-55261.

    15.Newton P,Copeland C,Martinez-Botas R,Seiler M.An audit of aerodynamic loss in a double entry turbine under full and partial admission.Int J Heat Fluid Flow 2012;33(1):70–80.

    16.J?cker M,Hillion FX,Fransson TH,Wa?hle′n U.Numerical unsteady flow analysis of a turbine stage with extremely large blade loads.New York:ASME;2001.Report No.:2001-GT-0260.

    17.Hudson ST,Zoladz TF,Dorney DJ.Rocket engine turbine blade surface pressure distributions:experiment and computations.J Propul Power 2003;19(3):364–73.

    18.Tokuyama Y,Funazaki K,Kato H,Shimiya N,Shimaga M,Uchiumi M.Computational analysis of unsteady flow in a partial admission supersonic turbine stage.New York:ASME;2014.604 Report No.:GT2014-26071.

    19.Tog RA,Tousi AM.Experimental and numerical investigation of design optimization of a partial admitted supersonic turbine.Propul Power Res 2013;2(1):70–83.

    20.Jia W,Liu HX.Numerical investigation of the interaction between upstream cavity purge flow and main flow in low aspect ratio turbine cascade.Chin J Aeronaut 2013;26(1):85–93.

    21.Wei ZJ,Qiao WY,Shi PJ,Chen PP,Zhao L.Tip-leakage flow loss reduction in a two-stage turbine using axisymmetric-casing contouring.Chin J Aeronaut 2014;27(5):1111–21.

    22.Zhang SW,MacManus DG,Luo JQ.Parametric study of turbine NGV blade lean and vortex design.Chin J Aeronaut 2016;29(1):104–16.

    21 October 2016;revised 19 December 2016;accepted 12 April 2017

    Available online 14 October 2017

    ?2017 Chinese Society of Aeronautics and Astronautics.Production and hosting by Elsevier Ltd.This is an open access article under the CCBY-NC-ND license(http://creativecommons.org/licenses/by-nc-nd/4.0/).

    *Corresponding author.

    E-mail address:wangjuetougao@163.com(J.WANG).

    Peer review under responsibility of Editorial Committee of CJA.

    少妇高潮的动态图| 午夜爱爱视频在线播放| 91在线精品国自产拍蜜月| 久久人人爽人人爽人人片va| 男人舔女人下体高潮全视频| 日韩av在线免费看完整版不卡| 伦理电影大哥的女人| 麻豆成人av视频| 中文字幕久久专区| 久久精品国产亚洲av天美| 一本一本综合久久| 一个人观看的视频www高清免费观看| 我要看日韩黄色一级片| 久久99热这里只有精品18| 亚洲综合色惰| 亚洲精品成人久久久久久| 小说图片视频综合网站| 欧美三级亚洲精品| 欧美日韩一区二区视频在线观看视频在线 | 我的女老师完整版在线观看| 在线天堂最新版资源| www.av在线官网国产| 午夜免费男女啪啪视频观看| 国产综合懂色| 蜜桃亚洲精品一区二区三区| 国产精品国产三级国产专区5o | 久久鲁丝午夜福利片| 天堂av国产一区二区熟女人妻| 亚洲av成人精品一区久久| 欧美成人午夜免费资源| 热99re8久久精品国产| 老司机影院成人| 精品欧美国产一区二区三| 久久午夜福利片| 亚洲国产色片| 最后的刺客免费高清国语| 日日干狠狠操夜夜爽| 91久久精品国产一区二区成人| av在线天堂中文字幕| 国产高清国产精品国产三级 | 中国国产av一级| 成年女人永久免费观看视频| 99热网站在线观看| 美女大奶头视频| 国产一区二区三区av在线| 在线观看av片永久免费下载| 国产精品一区二区在线观看99 | 中文字幕人妻熟人妻熟丝袜美| 天堂av国产一区二区熟女人妻| 亚洲精华国产精华液的使用体验| 天堂av国产一区二区熟女人妻| 国产日韩欧美在线精品| 男女那种视频在线观看| 久久精品国产鲁丝片午夜精品| 美女高潮的动态| 国产单亲对白刺激| 韩国av在线不卡| 韩国av在线不卡| 人妻少妇偷人精品九色| 熟女人妻精品中文字幕| 麻豆精品久久久久久蜜桃| 色视频www国产| 亚洲在线观看片| 免费看日本二区| 亚洲人成网站高清观看| 国产日韩欧美在线精品| 18+在线观看网站| 亚洲欧美成人精品一区二区| 免费一级毛片在线播放高清视频| 小蜜桃在线观看免费完整版高清| 联通29元200g的流量卡| 欧美一级a爱片免费观看看| 精品一区二区三区视频在线| 性插视频无遮挡在线免费观看| 久久久久精品久久久久真实原创| 亚洲成人精品中文字幕电影| 国产av在哪里看| 免费观看a级毛片全部| 又粗又爽又猛毛片免费看| 在线观看av片永久免费下载| 久久热精品热| av在线观看视频网站免费| 精品久久久久久电影网 | 国产69精品久久久久777片| 精品国内亚洲2022精品成人| 网址你懂的国产日韩在线| ponron亚洲| 亚洲真实伦在线观看| 纵有疾风起免费观看全集完整版 | 久久精品国产亚洲网站| 欧美xxxx黑人xx丫x性爽| 一夜夜www| 九九热线精品视视频播放| 国产欧美另类精品又又久久亚洲欧美| 男人舔女人下体高潮全视频| 国产av在哪里看| 国产熟女欧美一区二区| 熟妇人妻久久中文字幕3abv| 免费av观看视频| 久久久a久久爽久久v久久| 又爽又黄无遮挡网站| 99热这里只有是精品在线观看| 18禁裸乳无遮挡免费网站照片| 一级二级三级毛片免费看| 长腿黑丝高跟| 亚洲美女视频黄频| av在线观看视频网站免费| 蜜臀久久99精品久久宅男| 在线观看一区二区三区| 国产三级在线视频| 国产女主播在线喷水免费视频网站 | 国产一区二区亚洲精品在线观看| 在线观看一区二区三区| 日本与韩国留学比较| 国产乱人视频| 插逼视频在线观看| 99久久精品一区二区三区| 又爽又黄a免费视频| 国产精品永久免费网站| 国产久久久一区二区三区| 免费看日本二区| 人妻系列 视频| 最近2019中文字幕mv第一页| 国产精品嫩草影院av在线观看| 一本久久精品| 丝袜美腿在线中文| 欧美性感艳星| 99在线人妻在线中文字幕| 久久婷婷人人爽人人干人人爱| 2021天堂中文幕一二区在线观| 两个人的视频大全免费| 岛国毛片在线播放| av在线亚洲专区| 中文字幕av在线有码专区| 国产日韩欧美在线精品| 国产在线一区二区三区精 | 亚洲在线自拍视频| 亚洲久久久久久中文字幕| 少妇人妻精品综合一区二区| 国产亚洲精品久久久com| 少妇猛男粗大的猛烈进出视频 | 国产亚洲av片在线观看秒播厂 | 国产黄a三级三级三级人| 日韩欧美精品v在线| 成人毛片a级毛片在线播放| 精品不卡国产一区二区三区| 国产伦理片在线播放av一区| 亚洲av中文字字幕乱码综合| 国产高清三级在线| 九九爱精品视频在线观看| 人妻少妇偷人精品九色| 国产精品人妻久久久影院| 91aial.com中文字幕在线观看| 亚洲精品自拍成人| 国产免费一级a男人的天堂| 高清午夜精品一区二区三区| 欧美又色又爽又黄视频| 精品久久国产蜜桃| av国产久精品久网站免费入址| 人妻少妇偷人精品九色| 蜜桃亚洲精品一区二区三区| 亚洲av中文字字幕乱码综合| 亚洲成人av在线免费| 国产乱人视频| 99久久精品热视频| 久久精品综合一区二区三区| 最近2019中文字幕mv第一页| 亚洲丝袜综合中文字幕| 久久久久久久久久久免费av| 欧美日韩一区二区视频在线观看视频在线 | 3wmmmm亚洲av在线观看| 国产精品日韩av在线免费观看| 2021少妇久久久久久久久久久| 女的被弄到高潮叫床怎么办| 精品酒店卫生间| 我的女老师完整版在线观看| 国产精品嫩草影院av在线观看| 国产欧美日韩精品一区二区| 91久久精品电影网| 淫秽高清视频在线观看| 乱系列少妇在线播放| 啦啦啦韩国在线观看视频| 人妻夜夜爽99麻豆av| 国产精品三级大全| 男人和女人高潮做爰伦理| 精品久久久噜噜| 成人毛片60女人毛片免费| av在线观看视频网站免费| 久久人人爽人人爽人人片va| 99国产精品一区二区蜜桃av| 亚洲欧美一区二区三区国产| 亚洲精品乱码久久久v下载方式| 人妻制服诱惑在线中文字幕| av国产免费在线观看| 精品免费久久久久久久清纯| 欧美3d第一页| 亚洲国产日韩欧美精品在线观看| 中文字幕av成人在线电影| 长腿黑丝高跟| 亚洲最大成人手机在线| 又爽又黄a免费视频| 国产精品1区2区在线观看.| 91av网一区二区| 免费看美女性在线毛片视频| 午夜激情福利司机影院| 啦啦啦观看免费观看视频高清| 欧美三级亚洲精品| 99久久中文字幕三级久久日本| 亚洲aⅴ乱码一区二区在线播放| 亚洲精品456在线播放app| 97超碰精品成人国产| 日本免费一区二区三区高清不卡| 国产一区二区在线av高清观看| 国产精品久久久久久精品电影小说 | 亚洲av一区综合| 久久精品夜色国产| 成人综合一区亚洲| 久久精品国产亚洲av涩爱| 国产极品精品免费视频能看的| 亚洲欧美日韩无卡精品| 最近最新中文字幕免费大全7| 欧美一区二区亚洲| 亚洲综合精品二区| 99热全是精品| 老师上课跳d突然被开到最大视频| 在线播放无遮挡| 免费av不卡在线播放| 婷婷色麻豆天堂久久 | 国产一级毛片在线| 色尼玛亚洲综合影院| 国产在线男女| 国产精品久久久久久久电影| 亚洲激情五月婷婷啪啪| 男人舔女人下体高潮全视频| 国产精品一区二区三区四区久久| 舔av片在线| 成人无遮挡网站| 国产中年淑女户外野战色| 中文字幕制服av| 精品国内亚洲2022精品成人| 久久精品国产亚洲av涩爱| 黄色日韩在线| 日韩精品青青久久久久久| 久久久久精品久久久久真实原创| 长腿黑丝高跟| 欧美激情在线99| 亚洲国产精品成人久久小说| 嘟嘟电影网在线观看| 黄色配什么色好看| 日韩制服骚丝袜av| 狠狠狠狠99中文字幕| 亚洲国产色片| 午夜精品国产一区二区电影 | 久久久久国产网址| 91久久精品电影网| 亚洲欧洲国产日韩| 丰满乱子伦码专区| 欧美xxxx黑人xx丫x性爽| 建设人人有责人人尽责人人享有的 | 成人毛片a级毛片在线播放| 久久精品国产自在天天线| 亚洲av成人av| 69av精品久久久久久| 国产麻豆成人av免费视频| av线在线观看网站| 亚洲成av人片在线播放无| 精品熟女少妇av免费看| 亚洲成人中文字幕在线播放| 男女边吃奶边做爰视频| 六月丁香七月| 青春草国产在线视频| 一级爰片在线观看| 国产精品一区二区三区四区免费观看| ponron亚洲| 亚洲伊人久久精品综合 | 亚洲人与动物交配视频| 亚洲性久久影院| 国产不卡一卡二| 毛片一级片免费看久久久久| 精品人妻一区二区三区麻豆| 免费观看在线日韩| 99热网站在线观看| 1024手机看黄色片| 18+在线观看网站| 深夜a级毛片| 亚洲精品影视一区二区三区av| 99久久成人亚洲精品观看| 内射极品少妇av片p| 只有这里有精品99| 久久久精品欧美日韩精品| 99热全是精品| 一个人看视频在线观看www免费| 日本五十路高清| 亚洲国产欧美人成| 三级毛片av免费| av在线老鸭窝| 三级经典国产精品| 久久精品国产自在天天线| 18禁在线播放成人免费| 国产黄色视频一区二区在线观看 | 成人国产麻豆网| 亚洲av一区综合| 狂野欧美白嫩少妇大欣赏| 午夜精品一区二区三区免费看| 国产免费福利视频在线观看| 亚洲18禁久久av| 国产人妻一区二区三区在| a级一级毛片免费在线观看| 成人毛片60女人毛片免费| 亚洲18禁久久av| 水蜜桃什么品种好| 黄色欧美视频在线观看| 日本-黄色视频高清免费观看| 国产欧美另类精品又又久久亚洲欧美| 精品久久久久久久末码| 中文字幕免费在线视频6| 少妇熟女aⅴ在线视频| av在线观看视频网站免费| 精品久久国产蜜桃| 中文字幕精品亚洲无线码一区| 熟妇人妻久久中文字幕3abv| 精品国产三级普通话版| 岛国在线免费视频观看| 美女高潮的动态| 欧美高清成人免费视频www| ponron亚洲| 国产精品三级大全| 国产真实乱freesex| 欧美97在线视频| 夜夜看夜夜爽夜夜摸| 我的女老师完整版在线观看| 午夜老司机福利剧场| 久久国内精品自在自线图片| 国产中年淑女户外野战色| 尤物成人国产欧美一区二区三区| 99热这里只有是精品50| 人体艺术视频欧美日本| 亚洲av免费在线观看| 禁无遮挡网站| 国产 一区 欧美 日韩| 国产精品久久久久久久电影| 亚洲av.av天堂| av国产久精品久网站免费入址| 直男gayav资源| 美女脱内裤让男人舔精品视频| 亚洲美女视频黄频| 69av精品久久久久久| 国产一级毛片在线| 欧美激情国产日韩精品一区| 亚洲成人精品中文字幕电影| 成人漫画全彩无遮挡| 国产视频内射| 国产探花极品一区二区| av在线观看视频网站免费| 青春草亚洲视频在线观看| 麻豆乱淫一区二区| 少妇熟女aⅴ在线视频| 美女大奶头视频| 狠狠狠狠99中文字幕| 亚洲国产精品久久男人天堂| 人体艺术视频欧美日本| 九九久久精品国产亚洲av麻豆| 美女内射精品一级片tv| 久久人妻av系列| 日本黄色视频三级网站网址| 亚洲av免费在线观看| 国产精品久久久久久精品电影小说 | av视频在线观看入口| 高清在线视频一区二区三区 | 91av网一区二区| 日本av手机在线免费观看| 久久精品国产鲁丝片午夜精品| 亚洲成人av在线免费| 成人特级av手机在线观看| 日本三级黄在线观看| 日产精品乱码卡一卡2卡三| 国产成人精品婷婷| 国产精品熟女久久久久浪| av黄色大香蕉| 国产精品一区二区在线观看99 | 久久久精品大字幕| 国产欧美日韩精品一区二区| 色网站视频免费| 久久亚洲精品不卡| av视频在线观看入口| 久热久热在线精品观看| 熟妇人妻久久中文字幕3abv| 午夜福利在线在线| 亚洲欧洲日产国产| 1024手机看黄色片| 午夜精品一区二区三区免费看| 美女脱内裤让男人舔精品视频| 国产成人91sexporn| 欧美激情在线99| 国产精品久久久久久久久免| 欧美人与善性xxx| 久久久亚洲精品成人影院| 丰满乱子伦码专区| 秋霞在线观看毛片| 中国美白少妇内射xxxbb| 国产av在哪里看| 国产精品精品国产色婷婷| 青春草国产在线视频| 午夜免费男女啪啪视频观看| 久久欧美精品欧美久久欧美| 日韩强制内射视频| 麻豆精品久久久久久蜜桃| 国产亚洲精品av在线| 国产三级在线视频| 一个人免费在线观看电影| 国产精品蜜桃在线观看| 免费看光身美女| 精品人妻视频免费看| 搡老妇女老女人老熟妇| 亚洲欧美精品综合久久99| 九草在线视频观看| 日韩精品有码人妻一区| 天天躁日日操中文字幕| 高清在线视频一区二区三区 | 男女国产视频网站| 99久久成人亚洲精品观看| 亚洲乱码一区二区免费版| 最近手机中文字幕大全| av在线亚洲专区| 久久久久久久久久黄片| 99热全是精品| 老司机影院毛片| 观看美女的网站| 国产在线男女| 亚州av有码| 91精品国产九色| 亚洲中文字幕日韩| 69人妻影院| 久久久欧美国产精品| 长腿黑丝高跟| 高清毛片免费看| 亚洲精品日韩在线中文字幕| 亚洲av免费高清在线观看| 国产精品.久久久| 国产精品爽爽va在线观看网站| 又黄又爽又刺激的免费视频.| 精品久久久久久久久亚洲| 国产麻豆成人av免费视频| 国内精品一区二区在线观看| 国国产精品蜜臀av免费| 亚洲在线自拍视频| 亚洲欧洲国产日韩| 久热久热在线精品观看| 亚洲精品乱码久久久久久按摩| 婷婷色麻豆天堂久久 | or卡值多少钱| 九九爱精品视频在线观看| 国产精品一区二区性色av| 免费看美女性在线毛片视频| 亚洲精品乱码久久久v下载方式| av在线蜜桃| 一区二区三区四区激情视频| 亚洲在线自拍视频| 欧美日韩在线观看h| 91精品伊人久久大香线蕉| 99热精品在线国产| 精华霜和精华液先用哪个| 99久久精品国产国产毛片| 91狼人影院| 爱豆传媒免费全集在线观看| 国产乱来视频区| 三级男女做爰猛烈吃奶摸视频| 97在线视频观看| 免费黄网站久久成人精品| 在线观看av片永久免费下载| 亚洲国产精品成人综合色| 国产精品电影一区二区三区| 国产一区二区在线观看日韩| 1024手机看黄色片| 大香蕉久久网| 成人特级av手机在线观看| 亚洲国产精品国产精品| 日韩精品青青久久久久久| 好男人视频免费观看在线| 亚洲av中文字字幕乱码综合| 久久婷婷人人爽人人干人人爱| 人妻系列 视频| 国产极品精品免费视频能看的| 一夜夜www| 国产精品电影一区二区三区| 国内揄拍国产精品人妻在线| 99热6这里只有精品| 全区人妻精品视频| 一夜夜www| 日韩一区二区视频免费看| 日本wwww免费看| 国产亚洲一区二区精品| 亚洲电影在线观看av| 最近中文字幕高清免费大全6| 精品国内亚洲2022精品成人| 久久这里只有精品中国| 建设人人有责人人尽责人人享有的 | 淫秽高清视频在线观看| videos熟女内射| 亚洲精品成人久久久久久| 国产精品美女特级片免费视频播放器| av福利片在线观看| 在线观看美女被高潮喷水网站| 成人综合一区亚洲| 久久精品夜夜夜夜夜久久蜜豆| 午夜激情福利司机影院| 亚洲中文字幕日韩| 久久综合国产亚洲精品| 亚洲在线自拍视频| 热99在线观看视频| 欧美丝袜亚洲另类| av福利片在线观看| 欧美性猛交╳xxx乱大交人| 中文字幕av成人在线电影| 能在线免费看毛片的网站| 久久人人爽人人爽人人片va| 好男人视频免费观看在线| 97热精品久久久久久| 91精品国产九色| 搞女人的毛片| 久久亚洲国产成人精品v| 成人毛片60女人毛片免费| 久久久久久伊人网av| 一个人观看的视频www高清免费观看| 嫩草影院新地址| 1024手机看黄色片| 桃色一区二区三区在线观看| 国产在线一区二区三区精 | 久久亚洲国产成人精品v| 午夜福利在线在线| 嫩草影院入口| 少妇的逼水好多| 热99在线观看视频| a级毛色黄片| 两性午夜刺激爽爽歪歪视频在线观看| 国产av在哪里看| 午夜福利视频1000在线观看| 少妇熟女aⅴ在线视频| 久久久久久国产a免费观看| av黄色大香蕉| 日本色播在线视频| 久99久视频精品免费| 人体艺术视频欧美日本| 69人妻影院| 少妇的逼水好多| 中文精品一卡2卡3卡4更新| 桃色一区二区三区在线观看| 人妻夜夜爽99麻豆av| 最后的刺客免费高清国语| 黄色一级大片看看| 我要看日韩黄色一级片| 亚洲丝袜综合中文字幕| 亚洲18禁久久av| 尤物成人国产欧美一区二区三区| 成人三级黄色视频| 国产片特级美女逼逼视频| 亚洲在久久综合| 一边亲一边摸免费视频| 尤物成人国产欧美一区二区三区| www日本黄色视频网| 人妻系列 视频| 国产免费一级a男人的天堂| 日韩高清综合在线| 女人久久www免费人成看片 | 秋霞伦理黄片| 黑人高潮一二区| 久久久久久久久久久丰满| 深夜a级毛片| 久久精品国产亚洲网站| 国产女主播在线喷水免费视频网站 | 一区二区三区高清视频在线| 国产人妻一区二区三区在| 亚洲精品乱码久久久久久按摩| 精品无人区乱码1区二区| 汤姆久久久久久久影院中文字幕 | 婷婷六月久久综合丁香| 亚洲图色成人| 国产精品久久久久久精品电影| 亚洲国产精品国产精品| 亚洲av.av天堂| 中文字幕制服av| 亚洲真实伦在线观看| 一边摸一边抽搐一进一小说| 免费一级毛片在线播放高清视频| 91午夜精品亚洲一区二区三区| 好男人视频免费观看在线| 中文资源天堂在线| 免费播放大片免费观看视频在线观看 | 久久久久久国产a免费观看| 99热这里只有是精品50| 国内揄拍国产精品人妻在线| 又爽又黄无遮挡网站| 22中文网久久字幕| 亚洲av免费高清在线观看| 人体艺术视频欧美日本| 欧美成人a在线观看| 国产精品一区二区三区四区免费观看| 国产精品女同一区二区软件| 少妇的逼水好多| 久久久精品94久久精品| 天堂av国产一区二区熟女人妻| 国产精品麻豆人妻色哟哟久久 | 亚洲性久久影院| 免费搜索国产男女视频| 亚洲av中文av极速乱| 老女人水多毛片| 日本免费a在线| 日韩欧美 国产精品| 国产91av在线免费观看| 成人特级av手机在线观看| 插逼视频在线观看| 欧美成人免费av一区二区三区| 91久久精品国产一区二区三区| 九九在线视频观看精品| 欧美3d第一页|