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

    Experimental study of a controlled variable double-baffle distortion generator engine test rig

    2018-09-27 07:08:14AiguoXIAXudongHUANGWeiTUOMingZHOU
    CHINESE JOURNAL OF AERONAUTICS 2018年9期

    Aiguo XIA,Xudong HUANG,*,Wei TUO,Ming ZHOU

    aSchool of Aerospace Engineering,Tsinghua University,Beijing 100084,China

    bBeijing Aeronautical Technology Research Center,Beijing 100076,China

    KEYWORDS Distortion generator;Engine stability;Inlet flow;Inlet total-pressure distortion;Turbofan engines

    Abstract In order to explore the total-pressure distortion test assessment method for a turbofan engine,a Controlled Variable Double-Baffle Distortion Generator(CVDBDG)with a horizontal symmetry moving form was developed,which can adjust the steady-state and time–variant distortion separately in real time.The inlet total-pressure distortion test was conducted on an afterburner turbofan engine.The distortion parameters of CVDBDG and the instability characteristics of the engine were measured.The experimental data were modeled and analyzed by using back propagation artificial neural networks,and the work envelope of CVDBDG was obtained.Based on the analysis of the data on the engine’s instability,the properties of CVDBDG used for the stability assessment were preliminarily evaluated.The results show that CVDBDG can simulate both steady-state and time–variant distortions simultaneously in a range determined by three envelopes.Under the condition of symmetric double baffles,a critical depth of insertion exists,beyond which the symmetric baffles will generate an asymmetric flow field.In the case of double baffles,compared to a single baffle,the engine exhibited different instability characteristics.Based on CVDBDG,it is expected that more efficient engine stability and durability assessment methods can be developed.?2018 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/).

    1.Introduction

    Inlet total-pressure distortion assessment is an important part of an aero-engine test.Inlet total-pressure distortion can cause stalling or surges of the fan/compressor,induces synchronous vibration response of the blades of the fan and compressor blades,and change the temperature distribution of the combustor and turbine,which can lead to higher thermal stress on the turbine blades.1–3Since the end of the 1960s,after the F-111 aircraft engine surge failure,the problem of inlet distortion has been the subject of continuous research and attention in various countries,and extensive research has been conducted on a large number of ground-based and aerial testing technologies.4–7In recent years,more and more mobility and stealth performance have been designed into advanced fighters.In addition,large mobile flight and the use of S-bend and other shapes of inlet stealth designs have made the problem of air flow distortion more prominent.8–10

    Due to the adverse effects of inlet total-pressure distortion on engine performance,structural safety,and reliability,the assessment of inlet total-pressure distortion has become one of the most important aspects of aero-engine assessment.11Compared with the experiments of flight tests,propulsion wind tunnel tests,and free jet altitude tests,the use of the inlet total pressure distortion simulator in direct-connect altitude test facilities or ground test beds has better safety and economic performance.12–17

    Inlet total-pressure distortion assessment includes steadystate distortion and the time–variant distortion test.18The steady-state distortion simulators that have been developed mainly include a distortion screen19and air jet distortion generators.20,21Time–variant distortion simulators mainly include a controllable plug and ramp random frequency generators,discrete frequency generators,22and a splitting airfoil transient total-pressure distortion generator.23The distortion simulators,which can produce steady-state distortion and time–variant distortion simultaneously,mainly include a simulation board and a baffle distortion generator.24Currently,the distortion screen,simulation board,and baffle distortion generator are the most extensively used types in the research and engineering application of inlet total-pressure distortion.The distortion screen is simple to use,which only requires limited steady-state measurement,and can be designed for a selected distortion map.The Society of Automotive Engineers(SAE)has published a series of standards on inlet total-pressure distortion,18,25–27the main content of which is based on the distortion screen.The main disadvantage of the distortion screen is that it can only simulate steady-state distortion and cannot do real-time simulation.Thus,the simulation of several distortion maps is time-consuming and laborious.Both the simulation board and the baffle distortion generator can simulate both steady-state and time–variant distortion of total pressure,but the simulation board is a fixed shape that only applies to specific assessment situations.The movable baffle distortion generator utilizes a movable mechanism to enable real-time state conversion,and the simulated distortion state has a good similarity to the actual distortion characteristics.It is considered to be the most advanced inlet total-pressure distortion simulator developed in the last 100 years.28The movable baffle distortion generator was first used by Russia;however,there are very few public literature references on the distortion characteristics of the baffle distortion generator and the instability of the engine’s performance with a baffle distortion generator.In China,there also have been some studies on baffle distortion generators,including tests of engine parts,29,30testing in a ground-level test facility,1,31,32testing in direct-connect altitude test facilities,33flight tests,34and other related research work.Most of these experimental studies were based on one movable baffle,in a baffle position and engine operating condition,which can only produce a specific corresponding steady-state distortion and time–variant distortion,and it can only produce one low-pressure zone,which limits the baffle distortion generator for simulating actual complex and variable flight conditions. In practical applications, the phenomenon often occurs that the surge caused by the generator is inconsistent with the actual features.This is a problem that cannot be solved by a baffle distortion generator,which only assist in additional assessment of the special simulation board.In this case,adding another movable baffle is a natural idea.However,there are several questions that must be considered,e.g.,what the distortion characteristics of the two baffles would be,whether it can be used for stability evaluation,and what the instability characteristics of the engine would be.In the public literature,there are few studies on multi-baffle distortion generators.In Ref.35,four kinds of baffle structure forms,including a single flat plate,two symmetrical flat plates,and a fan plate,and two flat plates with a certain angle were tested on a two-stage,low-speed,axial- flow compressor test rig.It was found that the complex total-pressure distortion descriptor of single flat plate baffle was the largest by the same blockage coefficient,and using it as the standard baffle of distortion assessment is conservative and strict.The literature drew this conclusion without considering the performance of compressor stability for different baffles,so whether this conclusion conforms to reality must be analyzed further and studied in an engine test.Ref.24gave a kind of work envelope of the double-baffle distortion generator when the flow function was 0.5,but it did not explain the meaning of the boundary line mean or how the line influenced the engine’s stability.

    Based on the above introduction,a baffle distortion generator,compared with other inlet total-pressure distortion simulators,can simultaneously simulate steady-state and time–variant total-pressure distortion in real time.It is easy to use,and the simulated distortion is very close to the actual distortion characteristics.Thus,the basic data and theoretical guidance required for engineering applications can be acquired by conducting the development and testing of a two-baffle distortion generator to obtain the distortion characteristics and its effect on the stability of the engine.It can also provide a large number of real and credible distortion test cases available for the CFD simulation of flow field distortion.Based on the CVDBDG and the engine test system developed in this work,a turbofan engine was used as a test object in the ground test bed to test the distortion characteristics of the double baffles at different combinations,engine operating responses,and surge characteristics.In this paper,we mainly analyze the distorted flow field characteristics and engine stability with double baffles.

    2.Test equipment and engine

    2.1.Engine test bed

    The test was conducted on the indoor engine test bed at Beijing Aeronautical Technology Research Center.The test bed was a U-shaped structure,with an intake cross-section of 12 m×12 m.The maximum thrust test capacity was 250 kN,and the engine was suspended from the roof in the center of the intake air flow.

    2.2.Test engine

    The test engine was a two-spool after burning turbofan engine.The fan had an adjustable inlet guide vane,which is controlled according to the corrected rotational speed of the low-pressure spool.The inlet and first two stages of the high-pressure compressor had adjustable guide vanes,and the regulation law was determined according to the corrected rotational speed of the high-pressure spool.The nozzle was an adjustable Laval nozzle,and the regulation law of the nozzle of the throttling state was determined according to the corrected rotational speed of the high-pressure spool.The main parameters that were measured to determine the engine’s performance included the speed of the high pressure spool,the speed of the low pressure spool,thrust,fuel consumption,temperature and total pressure after the low pressure turbine,vibration of the inlet casing,vibration of the intermediate casing,total pressure after the fan,and total pressure and static pressure after the high pressure compressor.

    2.3.Controlled variable double-baffle distortion generator and flow field measurement

    Fig.1 shows the layout of the distortion test system and the engine.The distortion test system was in front of the engine,and it included the bell mouth,air meter,CVDBDG,Aerodynamic Interface Plane(AIP)measurement cross-section and air flow duct.

    The designed CVDBDG was installed at 3D in front of the engine,where D is the diameter of engine’s inlet,which was 905 mm.The distortion generator had two horizontal baffles.Both baffles were rectangular plates that were driven by two servo cylinders moving horizontally in a sealed box.Both Baffle 1 and Baffle 2 could travel to the relative insertion depth,with the ratio of length L1or L2to the inlet diameter D being 0.5.Fig.2 shows the location layout of the two baffles.

    Compared with a traditional electric motor and an ordinary hydraulic cylinder drive,the servo cylinder drive improved the control of the speed of the baffles and the accuracy of their positioning.The maximum adjusting speed of the servo cylinder was 100 mm/s,and the control error can reach 0.05%for different control speeds.The distortion generator was supported by a fixed platform that was fixed on the ground.There was a clearance of 2 mm between the rear of the distortion generator and the inlet section in the front of the engine,and this clearance was filled and sealed with flexible materials to isolate the distortion generator and the engine.This ensured that the air flow channel was smooth and the engine’s net thrust could be obtained during the distorted air flow.

    The inlet section between the engine and the distortion generator was provided with an engine inlet measuring section,which we referred to as the AIP,and was marked as Crosssection 2.The distance between the measuring section and the engine was 270 mm.There were a total of six pressure test probes and six wall static pressure measurement points uniformly distributed along the circumferential direction.The pressure test probes and the measurement points for the wall static pressure were interlaced,and the space was 12.5°.Each total-pressure test probe included six steady-state total pressure measuring points that were located in the centroids of equal areas,and one pulsating total-pressure measuring point,arranged at 0.9 times the radius of the flow channel.The average of the 36-point steady-state total-pressure measurements on the six probes was recorded as Pt2fav.The measured average of the 6-point wall static pressures was recorded as Ps2fav.

    Fig.2 Location layout of two baffles(viewed from engine).

    The air meter(named as Cross-section 1)and bell mouth sections were in the front of the distortion generator.The distance between the air meter measuring section and the distortion generator was 2D.Four total-pressure probes and four wall static pressure measurement points were distributed uniformly along the circumferential direction of the air meter with spaces of 45°.Each total-pressure probe included six steadystate total-pressure measuring points that were located in the centroids of equal areas.In total,there were 24 total pressure measurement points located in Cross-section 1,with the average value recorded as Pt1fav.The average of the 4-point wall static pressure was recorded as Ps1fav.The engine inlet air flow rate was calculated by the total pressure and static pressure of Cross-section 1.The bell mouth section had a standard,double twist line inlet,and the distance from the tangent point of the bell mouth to Cross-section 1 was 0.3D.The bell section,the air flow meter,and the distortion generator were connected rigidly and supported on the mounting platform.Four atmospheric temperature measurement points were distributed evenly outside the inlet port,giving the test parameters marked as Tt0and the average temperature marked as

    Among the tested parameters,the six-point pulsating totalpressure parameter of the AIP section and the vibration parameter of the inlet casing were collected and processed by a Gen3i high-speed acquisition system,and the sampling frequency was 25000 Hz.The other parameters were collected and processed by a PXI(PCI eXtensions for Instrumentation)acquisition system configured for the test bed,and it had a sampling frequency of 50 Hz.

    Fig.1 Diagram of engine and location of distortion test system.

    3.Definition of parameters and calculation method

    3.1.Distortion parameter

    There were two main characteristic parameters of the double baffles,i.e.,the relative insertion depth H1/H2of the baffle,defined as the ratio of length L1/L2to the inlet diameter D,and the relative insertion area S of the baffle,defined as the ratio of the insertion occlusion area to the inlet cross-section area.To facilitate understanding,the following S was expressed as a percentage of the relative insertion area and H1/H2was expressed as the percentage of the relative insertion depth.

    The main inlet total-pressure distortion evaluation indexes mainly included the AIP cross-section area-averaged totalpressure recovery coefficient σav,the steady-state circumferential distortion descriptorthe maximum ring steady-state radial distortion descriptorthe area-averaged turbulence εav,and the complex total-pressure distortion descriptor W.

    The area-averaged total-pressure recovery coefficient of AIP was calculated by

    The steady-state circumferential distortion descriptorwas calculated by

    where σ0is the averaged total-pressure recovery coefficient of AIP low-pressure sector.

    Area-averaged turbulence εav,was calculated by

    where Ndis the measurement points of the pulsation total pressure,and ε(i)is the turbulence of the ith pulsation totalpressure measurement point,calculated according to the following equation:

    where(ΔPt2)RMSis the root mean square of the pulsation total pressure,and Pt2is the time-averaged value of pulsation total pressure,with the averaging period of two seconds.

    The complex total-pressure distortion descriptor,W,which represents the total characteristic of the total-pressure distortion,is composed of a steady-state circumferential distortion descriptor and the area-averaged turbulence.

    The maximum ring steady-state radial distortion descriptoris the maximum steady-state radial distortion descriptor of the six total-pressure test rings of the AIP section.The ring steady-state radial distortion descriptorcan be calculated by

    where(Pt2av)iis the average value of six steady-state total pressure in the ith ring.

    3.2.Engine parameters

    The inlet air flow rate Wais calculated by the total pressure,the static pressure of Cross-section 1,and the ambient air temperature,based on the following equation:

    where KWis the boundary layer correction factor for the airflow meter,and K is an aerodynamic constant(for the conditions used in this paper,K is a constant with the value of 0.0404.);A is the area of Cross-section 1;λ1is the velocity coefficient of the air flow meter;q(λ1)is the flow function of the airflow meter.

    Corrected air flow rate is calculated by

    The air flow rate corrected to Cross-section 1 is identified as Wac1,and when it is corrected to AIP,it is identified as Wac2.In order to facilitate a relative comparison,the air flow rate parameters are represented by the ratio of the corrected air flow rate to the design air flow rate,and the corresponding parameters were represented with a line over the original name,e.g.,Wac2and

    The total-pressure ratio of the fan was calculated by

    where Pt21is the total pressure downstream from the fan.

    In this paper,stability index,KS,is used to reflect the situation of stability margin.It is defined as the ratio of fan totalpressure ratio on the stable boundary of the engine and the relative corrected air flow rate of the AIP section of the engine:

    4.Test method

    The main purpose of the experiment was to analyze the flow field characteristics and the engine’s response to the two baffles at different locations,so this paper takes two relative insertion depths of the baffle as the names of the test cases,i.e.,Case H1/H2.In order to facilitate the expression,the percentage symbol of H1/H2are removed.The comparison case of the benchmark is Case 0/0.In the experiment,we tested the quasisteady-state response of the engine from minimum to maximum flow conditions under different combinations of baffles.The performance of the 90%corrected high-pressure spool speed engine was tested by changing the positions of the baffles until the surge of the engine was tested without changing the engine’s throttle.

    5.Test results and analysis

    5.1.Matching work of engine and air inlet with CVDBDG

    The engine and the air inlet with CVDBDG are analyzed as a system.The rotating engine sucks air in from the air inlet with CVDBDG;at the same time,the change in the position of the adjustable baffle will also cause a change in the inlet flow field of the engine,resulting in a change of the state of the engine.In this paper,we select three parameters to jointly analyze the combining work of the engine and the distortion generator,i.e.,the corrected low-pressure spool rotational speed Nlc,the corrected air flow rateof Cross-section 1 at the entrance to the air inlet,and the area-averaged total-pressure recovery coefficient of AIP σav.Fig.3 shows the variation of Wac1and σavfor 18 different combinations of baffle positions in the test in the case thatchanges from 50%to 80%(analysis interval of 5%).The baffle position,H1/H2,and the corresponding percentage of relative insertion area,S,in the 18 test cases are included in the graphical representation.In order to facilitate comparison,the same color and logo are used to reflect the state of the single-baffle and the double-baffle with similar relative insertion areas.

    Fig.3 shows that both Wac1and σavtend to decrease gradually with the insertion of the baffle while the state of the engine is kept unchanged(i.e.,Nlcremains unchanged).As the position of the baffle remains unchanged,Wac1gradually increases,while σavgradually decreases as the condition of the engine improves.In the case of a low air flow rate and/or a low insertion area,the rate of the reduction of σavis close to a gentle linear decrease.Later,as the flow/insertion depth increases,there is a transition zone of the rate at which σavdeclines,and after going beyond this zone,σavdeclines faster.

    5.2.Results and analysis of distorted flow field

    Fig.4 shows the total-pressure recovery coefficient(σ)maps of the AIP section at 18 different baffle positions with an engine corrected high-pressure spool rotational speed Nhcof 88%.The relative insertion depth of the baffle for the corresponding location of the map is listed in the left and right positions of each map.According to the AIP section test,the 36-point total-pressure values are interpolated by the inverse-distance interpolation method.The maximum radius of these maps corresponds to the total-pressure measurement point in the maximum radial position.

    In this paper,the distorted flow field affected by the baffles is divided into three regions for analysis.The high-pressure region,where the total-pressure recovery coefficient is higher than the average total-pressure recovery coefficient,is either less affected or not affected by the distortion.The low pressure zone is where total-pressure recovery coefficient is lower than the average total-pressure recovery coefficient.In addition,the transition area is the junction area between the high-and low-pressure areas.Fig.5 shows the difference between the total-pressure recovery coefficient and the average total-pressure recovery coefficient of the AIP section(Δσ).All of the areas below zero are low-pressure areas.The range and strength of the low-pressure area can be seen more intuitively on the difference map.

    The map shows that,as the depth of a single-baffle increases,the total-pressure recovery coefficient of the affected area is decreasing.The low-pressure region extends gradually from the edge to the center of the circle,the range of the low-pressure region has its largest value near the edge,and the range decreases gradually from the edge to the center.The hump-type low-pressure region forms gradually on the side of the baffle.The variable trend of the high-pressure region is opposite that of the low-pressure area.With the insertion of the baffle,the range that is not influenced by distortion in the high-pressure region is reduced gradually.The totalpressure recovery coefficient in the high-pressure region also decreases gradually,and the saddle-shaped high-pressure concentration zone is formed gradually on the non-baffle side.There is an obvious transition area between the high-and lowpressure regions.With the gradual insertion of the baffle,the transition zone is gradually stabilized near the longitudinal center line of the section,forming a strip that has a totalpressure recovery coefficient that is basically equivalent to that of the average total pressure of the surface.

    Fig.4 Total-pressure recovery coefficient maps of AIP section.

    Fig.5 Total-pressure recovery coefficient difference maps.

    When double baffles are inserted simultaneously,the AIP section’s total-pressure recovery coefficient atlas also has three types of flow field characteristics,i.e.,(A)an asymmetric flow field under the asymmetric baffles,(B)a symmetric flow field under symmetric baffles,and(C)an asymmetric flow field under the symmetric baffles.When the baffles are inserted asymmetrically,two low-pressure areas are formed at the positions of the two baffle occlusions,and the low-pressure area has a hump shape.The high-pressure area is sandwiched between two low-pressure areas.When the baffles on both sides are inserted to a lower depth(in this case,the insertion depth is less than 15%),the high-pressure area is a symmetrical saddle around the longitudinal centerline of the cross section.When either side is inserted deeper than a certain position(in this case,the insertion depth is greater than 20%),the high-pressure area deviates the longitudinal centerline of the cross section and is offset to the side with lower insertion depth;the saddle-shaped high-pressure zone is also separated into two parts that are essentially symmetrical up and down.In the case of a symmetric baffle,there is a critical baffle depth.The flow field is basically symmetrical when the insertion depth of the bilateral baffle is below the critical depth(in this case,the insertion depth is below 15%),beyond which the flow field of the high-pressure district begins to favor one side,and it is separated into two parts with up and down symmetry.Similar to the single baffle,the double baffles also create a strip region in the longitudinal center of the cross-section where the total-pressure recovery coefficient is roughly equivalent to the coefficient of the area-average total-pressure recovery of the section.The mechanism of the asymmetric flow field under symmetric double baffles requires further study.

    5.3.Results of distorted parameters

    In this part,we tested the distortion parameters under different baffle positions from the engine idle state to the maximum air flow rate state.Because some baffle locations will surge before they reach the maximum engine air flow rate,only part of the baffle states can work on the maximum air flow state of the engine.

    Figs.6–9 show the variation of,W andwithfor 18 cases.In order to facilitate comparison,the same color and logo are used to reflect the single-baffle and double-baffle states with similar relative insertion areas.The fourth-order polynomial fitting method is used to fit W,and εavof the test points with different baffle positions.The fitted curves well reflect the trends of the variations of the parameters and their tendencies.Compared to Δσ0and W,and εavhave some jump fluctuations in several test points,which also reflects that these two parameters have fluctuations in the test.In general,the figures show that the four distortion indexes increase gradually with.As the depth of the baffle increases,the rate of variation also increases.

    Fig.6 compares the variation ofwithin the singlebaffle and double-baffle states.The figure shows that,with the same baffle position,gradually increases asincreases;with the same,increases as the relative insertion area increases.From the comparison of the states that the relative insertion areas are approximate in the figure,it is apparent that,for the same relative insertion area and,of the veneer status is higher than the status of the doublebaffle.For example,the relative insertion area of veneer status Case 25/0 is 19.6%,and thevalue for the correspondingexceeds both status Case 10/20(S=19.4%)and status Case 15/20(S=23.6%).This can be explained based on the effect of the baffle on σav,as explained above.As the relative insertion area increases,the σavnonlinearity decreases rapidly,i.e.,the total pressure decreases rapidly in a nonlinear fashion.In this way,for the same relative insertion area,the state of the veneer will produce a lower local low-pressure state,resulting in an increase in the degree of circumferential distortion.

    Fig.7 shows that,under the same baffle position,εavalso increases gradually with Wac1;with the sameWac1,εavincreases gradually as the relative insertion area increases.For the same relative insertion area andWac1,the εavvalue corresponding the single baffle is similar to or slightly higher than the status of the double-baffle,while the degree of excess is not as much as that of.Since W is the sum of and εav,W of a singlebaffle will be higher than that of a double-baffle under the same Wac1and relative insertion area.Fig.8 also shows this characteristic.

    The characteristic features ofare different from the previous three distortion parameters.Fig.9 shows thatshows a gradually increasing trend as the relative insertion area and Wac1increase for either a single-or a double-baffle.The variation ofis small when the relative insertion area is between 10%and 30%,and the variation ofis severe when the relative insertion area exceeds 30%.In a state in which the relative insertion area is less than about 10%,for the same Wac1,in the double-baffle state is larger than what it is in the single-baffle state.The two are close to each other between 10%and 30%,andof the single-baffle state is larger than that of the double-baffle state when the relative insertion area is greater than about 30%.

    Selecting some states with similar relative insertion areas as the object of the research,we studied the weight of the circumferential distortion descriptor and turbulence in the complex total-pressure distortion descriptor with the variation of the air flow rate.Fig.10 shows that the weight of the steadystate circumferential distortion descriptor in the complex total-pressure distortion descriptor increases gradually as the air flow rate increases.For the same relative insertion area,the single-baffle has a larger steady-state weight coefficient than the double-baffle.In the double-baffle states,the double-baffle state with a larger one-side baffle depth has a trend of having a larger steady-state weight coefficient.From another perspective,the weight of turbulence in the complex total-pressure distortion descriptor produced by the doublebaffle is higher than that produced by the single-baffle.

    5.4.Modeling and analysis of distortion test data

    In order to research the relationship and rules between the distortion parameters and the location of the baffle,we performed some data modeling and analysis using the experimental data.

    Due to the excellent data fitting ability of Back-Propagation Artificial Neural Networks(BPANN),we selected BPANN with a 3-layer structure to model the distortion test data, consisting of three input nodes, 10 hidden layer nodes, and three output nodes. The input parameters were Wac1,H1andH2,and the out put parameters were εav,and W.The input and output parameters were normalized from 0 to 1.We selected 2341 testing points from the test data,and we obtained an additional 2341 testing points by exchanging the data of H1and H2,so the total modeling data consisted of4682points.70%of the data were selected as training data,15%were the verification data,and the remaining 15% were the test data.Using the

    Fig.6 Variation of Δ

    Fig.7 Variation of εavwith Wac1.

    σ0with Wac1.Levenberg-Marquardt training algorithm,the network error performance obtained after training is shown in Fig.11.The abscissa is the difference between the model output data and the experimental data, and the ordinate is the number of samples for the corresponding error. The figure shows that the distortion network model based on BPANN has high enough precision for both the training data and the test data.The model can be used for the analysis of the test data.

    Fig.8 Variation of W with

    Wac1.

    Fig.9 Variation of Δ

    σrmwith Wac1.

    Fig.10 Weight of steady-state circumferential distortion descriptor in complex total-pressure distortion descriptor varying with air flow rate.

    Fig.11 Distortion network model based on BPANN error histogram.

    The BPANN model presented above was used to conduct a comparative analysis of the single and double baffles.We selected the state in whichwas 0.9 and H1and H2were changing in the range of 0–20%to study the working characteristics of double baffles.Fig.12 shows the variation of εavandwith different baffle combinations.The pink scatter marks in the figure are the distortion parameter combination range formed by the two baffles working in the range of 0–20%relative insertion depths.Such an area is an enclosed area surrounded by three border lines.The analysis shows that the lower line marked with a solid red line is a work envelope with a baffle position of 0 and another baffle changing from 0 to 20%.This line reflects the ability of a single-baffle,which is defined as the single-baffle working line in this article.The right blue boundary line is the one-sided baffle relative insertion depth of 20%,and the other side is from 0 to 20%;in this paper,we define this line as the single largest baffle working line.The left boundary line of cyan is the working line of a symmetric baffle that changes simultaneously from 0 to 20%,which is defined as the symmetrical baffle working line.

    Fig.12 Example of variation of εavand with different combinations of baffles(Wac1=0.9).

    In addition,the difference between the single-baffle and double-baffle was analyzed from the view of the relative insertion area.The black°-mark in the lower right corner is a single-baffle status with one side baffle at 0 and the other side at 20%,and the relative insertion area of the state is 14.2%.Based on the relative insertion area of 14.2%,we can look for some double-baffle states for which the total insertion area is 14.2%.The black solid line in Fig.12 shows the characteristic line that consists of various double-baffle states with a relative insertion area of 14.2%.It can be seen that the equal relative insertion area states formed by the double baffles are a working line that extends from the working point of the single baffle to the left boundary.In this process, εavanddecline gradually,anddeclines faster.This means that a single-baffle state will correspond to a symmetrical doublebaffle state,with minimum εavand.Similarly,we can analyze the relationship between single-baffle and double-baffle in the case of the other insert areas.The equal relative insertion area working lines of Cases 0/18 and 0/16 are also shown in Fig.12,and it is apparent that they are similar to Case 0/20.The difference is that,as the insertion area decreases,the degree of reduction of the turbulence decreases.

    Further analysis was performed on the influence of the status of the engine on the work envelope.Fig.13 shows the change of the working envelope of Wac1from 0.5 to 0.9 when the relative insertion depth of the baffle is kept at 20%or less.The figure shows that,with the same baffle relative insertion depth,as the engine state increases,the double-baffle simulating envelope expands gradually.

    The analysis above shows that,after adding another baffle,the simulation capability of the baffle distortion generator is extended from only one line of a single baffle to a working area surrounded by three working lines.In the example of Fig.12,the maximumsimulation capability of the double-baffle increases by 198%compared with that of the single-baffle;εavhas an increase of 120%,and the simulation range of the distortion is increased significantly.

    5.5.Results and analysis of engine instability test

    In the test,we also conducted a comparative study of engine instability with different baffle forms.We chose the condition of Nhcaround 90%in this study,and we conducted engine surge test with different forms,i.e.,single-baffle,symmetrical double-baffle,and asymmetric double-baffle.Table 1 lists the main parameters of the engine and the distortion at the last operating point before the engine became unstable under the four conditions.Based on the data in the table,we conclude the following:

    Fig.13 Change of working envelope of Wac1from 0.5 to 0.9(for a relative insertion depth not more than 20%).

    (1)The operating conditions of the engines in the table are close to each other.The difference between the maximum and minimum values of Nlcis 2.7%,and the difference between the maximum and minimum corresponding critical complex total-pressure distortion descriptor(the value at the last operating point before the engine becomes unstable)is 11%.A complex totalpressure distortion descriptor is used to evaluate the stability of the engine with a single baffle.The experimental results show that the critical complex total-pressure distortion descriptor also has a certain degree of engineering accuracy if the double-baffle is used for the evaluation of the stability of the engine.

    (2)By comparing the two cases of the single-baffle state,i.e.,Case 50/0 and Case 43.5/0,with the two doublebaffle cases,i.e.,Case 31.5/31.5 and Case 42.9/10,it is apparent that,irrespective of whether the states are single-baffle or double-baffle states,when the engine is unstable,the operating conditions and the corresponding critical complex total-pressure distortion descriptor are inversely proportional to each other.As the engine’s operating speed increases,its critical complex totalpressure distortion descriptor decreases.

    (3)By comparing single-baffle Case 43.5/0 with two doublebaffle cases,i.e.,Case 31.5/31.5 and Case 42.9/10,it is apparent that the critical complex total-pressure distortion descriptor is comparatively higher when the operating speed is slightly higher.As a result,it can be inferred that,for the same engine operating condition,the critical complex total-pressure distortion descriptor in the double-baffle state will be smaller than that in the single-baffle state when the engine loses stability.

    (4)No obvious correlation is found between the obtainedand stability.

    (5)According to the stability index KS,Case 31.5/31.5 is about 3%lower than the other states.Although there is no further support for the test data,based on this example,we can infer that,if double-baffle is used for the stability assessment for the same engine operating condition,the estimated engine stability margin will tend to decrease gradually as the baffle symmetry increases.

    The following is a further analysis of the characteristics of air flow when the engine is unstable.According to wellknown researchers,the surge process can be divided into classic surge and deep surge,depending on whether there is back flow during the surge.In order to investigate back flow at the inlet to the engine,we set the ratio of the area-averaged total pressure and static pressure on the AIP as a basis for judging the back flow.A ratio of less than one indicates that back flow has occurred.Fig.14 shows the changing state of Pt2fav/Ps2favin the four cases in Table 1 during the unstable operation of the engine.Fig.14 shows that all of the cases except Case 42.9/10 have back flow,which should have occurred in deep surge states.Fig.15 shows the changing process of total pressure on AIP(take the outermost measure point of No.2 steady-state total-pressure probe as a representativeexample,labeled as Pt2_62),the total pressure after the fan Pt21,and the total pressure after the high-pressure compressor(labeled as Pt3)during engine instability.The four samples in Fig.15 show three different characteristics.Cases 50/0 and 43.5/0 show that Pt21and Pt3continue to decrease until the surge disappears(note:All of the experiments have an automatic anti-surge system,and surges will automatically exit just after it occurs).In Case 31.5/31.5,Pt3continues to decline until the surge disappears,while Pt21has a process of decline,i.e.,it declines,then rises,and then declines again.Case 42.9/10 shows a synchronous oscillation of Pt2_62,Pt21and Pt3.In addition,from Case 50/0 in Fig.15,it is apparent that Pt21declines before Pt3,and its decline is accompanied by the increase of Pt2_62;thus,it can be concluded that the surge occurs first at the fan.In Cases 43.5/0 and 42.9/10,the total pressure on AIP increases earlier than the decline of Pt21begins.We also infer that the surge occurs first at the fan.In Case 31.5/31.5,Pt2_62increases,and Pt21and Pt3decrease almost simultaneously.It is difficult to distinguish whether the surge occurs at the fan or the high-pressure compressor from the existing data that are sample at 50 Hz.Further assessment of this issue will require the analysis of dynamic pressure data.According to Fig.14 and Fig.15,it can be assumed that Cases 50/0 and 43.5/0 indicate a deep surge that occurs first at the fan before it occurs at the high-pressure compressor.Case 42.9/10 shows classic surge characteristics with rotational stall.Case 31.5/31.5 shows a deep surge,but it differs from the performance of Cases 50/0 and 43.5/0.

    Table 1 Main parameters of engine and distortion at the last operating point before engine becoming unstable.

    Fig.14 Change of Pt2fav/Ps2favduring unstable state of engine.

    6.Conclusions and outlook

    In this paper,using a CVDBDG and an after burning turbofan engine,we conducted research on the inlet total-pressure distortion on an indoor engine test bed,and we measured the distortion simulation characteristics of the CVDBDG and the instability characteristics of the engine.The method of back propagation artificial neural networks was used to model the distortion test data and to analyze the working envelope of CVDBDG.Through the analysis of the engine instability data,we obtained guidance concerning the use of CVDBDG to evaluate the stability of an engine.Through the above experimental and modeling analyses,we can draw the following conclusions and prospects for future engineering applications and research of the CVDBDG:

    (1)CVDBDG can simulate steady-state and time–variant total-pressure distortion in real time.The process of changing two baffles in the same maximum insertion depth range can form a simulated work area of steady state and time–variant total-pressure variations that consist of a single-baffle working line,a symmetrical double-baffle working line,and a single maximum baffle working line.

    Fig.15 Change of Pt2_62,Pt21and Pt3during unstable state of engine.

    (2)Compared with the movable single-baffle distortion generator,both the capability and range of the distortion simulation of steady-state and time–varianttotal pressure of CVDBDG are extended significantly.When another baffle is added,a simulated working point of the movable single-baffle distortion generator is expanded to a simulated working line while the relative insertion area is kept unchanged.When the relative insertion depth is kept constant,an analog working line of a movable single-baffle distortion generator is expanded into a simulated working area.

    (3)When the engine becomes unstable during different combination styles of the double-baffle simulation conditions,the characteristics are different from the single-baffle.Based on this feature,double-baffle is expected to be used in evaluating some kinds of engine surges to verify whether the engine has some potential destabilization hazards through different kinds of baffle combinations.This is a direction worth further studying.

    (4)Under the same engine operating conditions,the critical complex total-pressure distortion descriptor that led to instability of the engine decreases with the increase of the symmetry of the double-baffle.In addition,the estimated surge margin also shows a downward trend.From this perspective,if the CVDBDG is used to assess stability,the evaluation results tend to become more stringent as the symmetry increases.

    (5)Under the condition of a symmetric double-baffle,there exists a critical depth of insertion beyond which the symmetric baffle will generate an asymmetric flow field.Further study of the mechanism of the critical insertion depth can achieve more accurate simulation of analog control using the double-baffle approach.

    (6)CVDBDG can produce more complex and diverse flow fields,which also means that excitation air flows with a wide range of frequencies can be generated in this way.Perhaps it can be a new way to assess the durability of engines.

    (7)The CVDBDG greatly extends its capabilities compared with a single movable baffle distortion generator,and it is expected to become a new and more effective means of assessing the stability and durability of engines.

    Acknowledgements

    The engine test was supported by the Beijing Aeronautical Technology Research Center.The authors express their appreciation for everyone who was involved in and contributed to the experiment.

    黄片大片在线免费观看| 久久婷婷人人爽人人干人人爱 | 热re99久久国产66热| 热re99久久国产66热| 久久香蕉激情| 欧美成人一区二区免费高清观看 | 亚洲一卡2卡3卡4卡5卡精品中文| 叶爱在线成人免费视频播放| АⅤ资源中文在线天堂| 国产高清激情床上av| 色播在线永久视频| 色综合欧美亚洲国产小说| 日韩有码中文字幕| 色播亚洲综合网| 亚洲中文av在线| 国产精品免费一区二区三区在线| 午夜精品国产一区二区电影| 亚洲电影在线观看av| 久久精品成人免费网站| 性色av乱码一区二区三区2| 夜夜看夜夜爽夜夜摸| 日韩国内少妇激情av| 日韩大尺度精品在线看网址 | 亚洲一码二码三码区别大吗| 两个人视频免费观看高清| 在线视频色国产色| 欧美亚洲日本最大视频资源| 欧美国产精品va在线观看不卡| 中文字幕色久视频| 成年人黄色毛片网站| 最近最新中文字幕大全免费视频| 丰满的人妻完整版| 亚洲国产毛片av蜜桃av| 韩国精品一区二区三区| 国产1区2区3区精品| 亚洲一区二区三区不卡视频| 黄色成人免费大全| 亚洲一区二区三区不卡视频| 久久久久久免费高清国产稀缺| 嫩草影视91久久| 91成人精品电影| 欧美日本亚洲视频在线播放| 日韩成人在线观看一区二区三区| 免费在线观看影片大全网站| 精品少妇一区二区三区视频日本电影| 成人18禁在线播放| 97超级碰碰碰精品色视频在线观看| 国产国语露脸激情在线看| 日本a在线网址| 美女高潮到喷水免费观看| 欧美激情高清一区二区三区| 每晚都被弄得嗷嗷叫到高潮| 777久久人妻少妇嫩草av网站| 亚洲专区字幕在线| 久久久久久久久中文| 一级毛片精品| 首页视频小说图片口味搜索| 亚洲五月婷婷丁香| 免费看十八禁软件| 啦啦啦免费观看视频1| 在线视频色国产色| 日韩精品免费视频一区二区三区| 欧美+亚洲+日韩+国产| ponron亚洲| 可以在线观看的亚洲视频| 亚洲国产欧美一区二区综合| 黄色女人牲交| 精品久久久久久,| 我的亚洲天堂| 精品欧美一区二区三区在线| 久久精品aⅴ一区二区三区四区| av在线播放免费不卡| 久久国产精品男人的天堂亚洲| 18禁观看日本| 免费一级毛片在线播放高清视频 | 一区二区三区高清视频在线| 侵犯人妻中文字幕一二三四区| 欧美乱色亚洲激情| 麻豆国产av国片精品| 国产熟女午夜一区二区三区| 亚洲午夜理论影院| 99re在线观看精品视频| 给我免费播放毛片高清在线观看| 男男h啪啪无遮挡| 国产av在哪里看| 天天添夜夜摸| 欧美成狂野欧美在线观看| 亚洲男人的天堂狠狠| 日韩欧美国产在线观看| 国产成人精品久久二区二区免费| 日韩精品中文字幕看吧| 激情视频va一区二区三区| ponron亚洲| av视频免费观看在线观看| 一进一出好大好爽视频| 精品国产乱码久久久久久男人| 女人被狂操c到高潮| 成人永久免费在线观看视频| 日韩成人在线观看一区二区三区| 久久精品成人免费网站| 一区二区三区国产精品乱码| 国产真人三级小视频在线观看| 国产av在哪里看| 人人妻人人澡欧美一区二区 | 熟妇人妻久久中文字幕3abv| 欧美色欧美亚洲另类二区 | 国产激情久久老熟女| 亚洲avbb在线观看| 又黄又爽又免费观看的视频| 日本 av在线| 欧美日韩亚洲综合一区二区三区_| 69精品国产乱码久久久| 成人欧美大片| 国产精品久久久久久亚洲av鲁大| 亚洲精品国产精品久久久不卡| www.www免费av| 精品欧美一区二区三区在线| 久久精品亚洲精品国产色婷小说| 宅男免费午夜| 欧美色欧美亚洲另类二区 | 满18在线观看网站| 老司机午夜福利在线观看视频| 亚洲色图综合在线观看| 精品国产亚洲在线| 中文字幕人妻丝袜一区二区| 一边摸一边抽搐一进一小说| 久久精品国产99精品国产亚洲性色 | 黑人操中国人逼视频| 老熟妇仑乱视频hdxx| 激情视频va一区二区三区| 亚洲国产看品久久| 色播亚洲综合网| 亚洲免费av在线视频| 黄频高清免费视频| 久久久精品欧美日韩精品| 精品无人区乱码1区二区| 在线播放国产精品三级| 国产免费男女视频| 精品国产国语对白av| 欧美一区二区精品小视频在线| 亚洲五月色婷婷综合| 欧美黑人精品巨大| 一个人免费在线观看的高清视频| 亚洲黑人精品在线| www.精华液| 三级毛片av免费| 日日夜夜操网爽| 精品国产亚洲在线| 久久久久久国产a免费观看| 中亚洲国语对白在线视频| 成人国语在线视频| 亚洲国产欧美网| 亚洲国产欧美日韩在线播放| 亚洲精品在线观看二区| 天天添夜夜摸| 女人被狂操c到高潮| 真人做人爱边吃奶动态| 欧美+亚洲+日韩+国产| 午夜a级毛片| 国产精品久久久久久人妻精品电影| www国产在线视频色| 精品一区二区三区av网在线观看| 亚洲午夜精品一区,二区,三区| 国产蜜桃级精品一区二区三区| 日本精品一区二区三区蜜桃| 男人操女人黄网站| 啦啦啦韩国在线观看视频| 一区二区三区激情视频| 在线观看一区二区三区| 夜夜爽天天搞| 欧美中文综合在线视频| 欧美日韩瑟瑟在线播放| 国产免费av片在线观看野外av| 人成视频在线观看免费观看| 欧美大码av| 女同久久另类99精品国产91| 欧美午夜高清在线| 亚洲aⅴ乱码一区二区在线播放 | 欧美大码av| 欧美日韩精品网址| 欧美日本亚洲视频在线播放| 亚洲欧美一区二区三区黑人| 国产麻豆成人av免费视频| 国产激情欧美一区二区| 日本黄色视频三级网站网址| 久久久国产精品麻豆| 免费在线观看黄色视频的| 91老司机精品| 日韩欧美国产一区二区入口| 一a级毛片在线观看| 精品免费久久久久久久清纯| 别揉我奶头~嗯~啊~动态视频| 黄片播放在线免费| 亚洲 欧美 日韩 在线 免费| 久久香蕉国产精品| 国产精品亚洲av一区麻豆| 国产成+人综合+亚洲专区| 色在线成人网| cao死你这个sao货| 人人妻人人澡人人看| 如日韩欧美国产精品一区二区三区| 黑人欧美特级aaaaaa片| 亚洲国产精品合色在线| 自拍欧美九色日韩亚洲蝌蚪91| 夜夜爽天天搞| 麻豆国产av国片精品| 中文字幕最新亚洲高清| 人妻丰满熟妇av一区二区三区| 欧美中文日本在线观看视频| 欧美成人免费av一区二区三区| 久久午夜亚洲精品久久| 给我免费播放毛片高清在线观看| e午夜精品久久久久久久| 免费女性裸体啪啪无遮挡网站| 亚洲色图 男人天堂 中文字幕| 91国产中文字幕| 在线av久久热| 亚洲精华国产精华精| 极品人妻少妇av视频| 国产极品粉嫩免费观看在线| 国内精品久久久久久久电影| 99久久综合精品五月天人人| 很黄的视频免费| 很黄的视频免费| 搞女人的毛片| 国产精品久久久av美女十八| 精品久久久久久久毛片微露脸| 婷婷丁香在线五月| 日韩欧美国产在线观看| 欧美日韩黄片免| 1024香蕉在线观看| 久久久久国内视频| 50天的宝宝边吃奶边哭怎么回事| 国产三级在线视频| 日本a在线网址| 黄色丝袜av网址大全| 黄色成人免费大全| 久久精品人人爽人人爽视色| 日本撒尿小便嘘嘘汇集6| 91麻豆av在线| 亚洲在线自拍视频| 国产三级在线视频| 色播亚洲综合网| 国产精品一区二区在线不卡| 在线av久久热| 乱人伦中国视频| 又大又爽又粗| 欧美精品啪啪一区二区三区| 欧美日本中文国产一区发布| 亚洲欧美激情综合另类| 中文字幕高清在线视频| 国产精品99久久99久久久不卡| 国产野战对白在线观看| 桃红色精品国产亚洲av| 可以在线观看毛片的网站| 亚洲av成人一区二区三| 欧美日本亚洲视频在线播放| 少妇裸体淫交视频免费看高清 | 精品国内亚洲2022精品成人| 亚洲三区欧美一区| 99国产精品免费福利视频| 别揉我奶头~嗯~啊~动态视频| 女人爽到高潮嗷嗷叫在线视频| 国产一区在线观看成人免费| av福利片在线| 国产精品国产高清国产av| 国产欧美日韩一区二区三区在线| 搡老岳熟女国产| 成人手机av| 美国免费a级毛片| 亚洲一区二区三区色噜噜| 少妇裸体淫交视频免费看高清 | 日韩一卡2卡3卡4卡2021年| 亚洲国产精品合色在线| 久久香蕉国产精品| 三级毛片av免费| 老司机深夜福利视频在线观看| 两个人免费观看高清视频| 成年人黄色毛片网站| 欧美乱码精品一区二区三区| 亚洲第一青青草原| 亚洲成a人片在线一区二区| 精品一区二区三区四区五区乱码| 香蕉丝袜av| 中文字幕最新亚洲高清| 午夜影院日韩av| www.www免费av| 在线国产一区二区在线| 精品国产美女av久久久久小说| 男女下面插进去视频免费观看| 精品久久久久久成人av| 12—13女人毛片做爰片一| 俄罗斯特黄特色一大片| 国产亚洲欧美精品永久| 久久婷婷成人综合色麻豆| 亚洲精品美女久久av网站| 非洲黑人性xxxx精品又粗又长| 757午夜福利合集在线观看| 村上凉子中文字幕在线| 香蕉国产在线看| 成人欧美大片| 国产精品久久久久久精品电影 | 亚洲欧美精品综合一区二区三区| 狠狠狠狠99中文字幕| 国产1区2区3区精品| 人人妻人人爽人人添夜夜欢视频| 午夜久久久久精精品| 在线观看免费视频日本深夜| 国产精品 国内视频| 一级毛片高清免费大全| 中文字幕人妻丝袜一区二区| 国产成人av激情在线播放| 久久天堂一区二区三区四区| 日本三级黄在线观看| 久久人人精品亚洲av| 精品熟女少妇八av免费久了| 亚洲精品中文字幕在线视频| 国产欧美日韩一区二区精品| 国产99白浆流出| 99riav亚洲国产免费| 欧美国产精品va在线观看不卡| 成人免费观看视频高清| 国产成人免费无遮挡视频| 国产免费av片在线观看野外av| 成人三级做爰电影| 啦啦啦免费观看视频1| 国内精品久久久久久久电影| 黑人巨大精品欧美一区二区mp4| 色播亚洲综合网| 在线十欧美十亚洲十日本专区| 俄罗斯特黄特色一大片| 国产精品亚洲美女久久久| videosex国产| 99国产精品99久久久久| 久久久久久久午夜电影| 亚洲国产欧美网| 视频区欧美日本亚洲| 日韩大尺度精品在线看网址 | 91九色精品人成在线观看| 亚洲精品国产精品久久久不卡| 亚洲精品国产色婷婷电影| 老熟妇仑乱视频hdxx| 亚洲精品中文字幕一二三四区| 狂野欧美激情性xxxx| 中文字幕高清在线视频| 成人特级黄色片久久久久久久| 婷婷丁香在线五月| 亚洲avbb在线观看| 久久精品91蜜桃| 男女做爰动态图高潮gif福利片 | 法律面前人人平等表现在哪些方面| 久久香蕉精品热| 亚洲精品国产精品久久久不卡| 日韩欧美三级三区| 国产麻豆成人av免费视频| 亚洲午夜精品一区,二区,三区| 男人舔女人下体高潮全视频| 99香蕉大伊视频| 亚洲va日本ⅴa欧美va伊人久久| 黄片小视频在线播放| 久久久精品国产亚洲av高清涩受| 国产精品亚洲美女久久久| 欧美黄色片欧美黄色片| 国产区一区二久久| 女性被躁到高潮视频| 午夜福利,免费看| 亚洲第一欧美日韩一区二区三区| 亚洲人成网站在线播放欧美日韩| 亚洲中文av在线| 人妻久久中文字幕网| 国产成人精品久久二区二区免费| 亚洲自拍偷在线| 欧美乱妇无乱码| 亚洲熟女毛片儿| a级毛片在线看网站| 亚洲午夜精品一区,二区,三区| 69av精品久久久久久| 久久精品亚洲精品国产色婷小说| 国产精品自产拍在线观看55亚洲| 两性午夜刺激爽爽歪歪视频在线观看 | 久久青草综合色| 国产精品一区二区三区四区久久 | 无限看片的www在线观看| 国产精品精品国产色婷婷| 精品高清国产在线一区| 久久久精品欧美日韩精品| 咕卡用的链子| 最近最新中文字幕大全免费视频| 亚洲欧美精品综合一区二区三区| 大香蕉久久成人网| 亚洲电影在线观看av| 男人的好看免费观看在线视频 | 手机成人av网站| 电影成人av| 精品不卡国产一区二区三区| 国产av又大| 中文亚洲av片在线观看爽| a在线观看视频网站| 在线观看www视频免费| 中国美女看黄片| 国产精品免费视频内射| 欧美乱妇无乱码| 国产高清视频在线播放一区| 一进一出抽搐动态| 女人高潮潮喷娇喘18禁视频| 成年女人毛片免费观看观看9| 精品国产一区二区三区四区第35| 日本三级黄在线观看| 国产欧美日韩一区二区三| 九色国产91popny在线| 人人妻人人爽人人添夜夜欢视频| 窝窝影院91人妻| 久久久久久国产a免费观看| 亚洲国产欧美网| 免费看美女性在线毛片视频| 精品福利观看| 久久午夜亚洲精品久久| 神马国产精品三级电影在线观看 | 国产国语露脸激情在线看| 神马国产精品三级电影在线观看 | 国产精品秋霞免费鲁丝片| 乱人伦中国视频| 少妇的丰满在线观看| 美女免费视频网站| 午夜久久久久精精品| 99国产精品免费福利视频| 久久青草综合色| 国产一区二区三区在线臀色熟女| 欧美 亚洲 国产 日韩一| 免费搜索国产男女视频| 亚洲成人免费电影在线观看| 国产极品粉嫩免费观看在线| 久久国产乱子伦精品免费另类| 欧美日韩亚洲国产一区二区在线观看| www.999成人在线观看| 中文字幕色久视频| 天天躁夜夜躁狠狠躁躁| 国产色视频综合| 日本在线视频免费播放| 色老头精品视频在线观看| 亚洲精品美女久久av网站| √禁漫天堂资源中文www| 亚洲欧美精品综合一区二区三区| 露出奶头的视频| 成人欧美大片| 长腿黑丝高跟| 欧美丝袜亚洲另类 | 色播在线永久视频| 在线免费观看的www视频| 午夜福利成人在线免费观看| 涩涩av久久男人的天堂| 69精品国产乱码久久久| 亚洲国产欧美日韩在线播放| 男人的好看免费观看在线视频 | 国产熟女午夜一区二区三区| 国产精品一区二区精品视频观看| 欧美日韩亚洲综合一区二区三区_| 亚洲av片天天在线观看| 亚洲久久久国产精品| 久久草成人影院| 国产精品久久久久久亚洲av鲁大| 午夜精品久久久久久毛片777| 欧美精品啪啪一区二区三区| 99国产精品免费福利视频| 亚洲天堂国产精品一区在线| 日本免费一区二区三区高清不卡 | 亚洲精品国产区一区二| 亚洲精华国产精华精| 欧美日韩精品网址| 满18在线观看网站| 丁香六月欧美| 美女高潮喷水抽搐中文字幕| 国产片内射在线| 久久婷婷人人爽人人干人人爱 | 一a级毛片在线观看| 久久久久国产一级毛片高清牌| 男人舔女人的私密视频| 99riav亚洲国产免费| 国产亚洲av高清不卡| 中国美女看黄片| 亚洲欧美日韩另类电影网站| a级毛片在线看网站| 超碰成人久久| 日日干狠狠操夜夜爽| 欧美乱码精品一区二区三区| 99re在线观看精品视频| 欧美老熟妇乱子伦牲交| 欧美丝袜亚洲另类 | 欧美日韩精品网址| 亚洲一码二码三码区别大吗| 欧美日韩瑟瑟在线播放| 大型av网站在线播放| 欧美 亚洲 国产 日韩一| 搡老妇女老女人老熟妇| 精品久久久精品久久久| 成在线人永久免费视频| 久久国产亚洲av麻豆专区| 国产亚洲欧美在线一区二区| 久久国产亚洲av麻豆专区| 女生性感内裤真人,穿戴方法视频| 国产精品一区二区免费欧美| 久久精品成人免费网站| 国产免费男女视频| 18禁美女被吸乳视频| 日日摸夜夜添夜夜添小说| 在线天堂中文资源库| 日日夜夜操网爽| 国产一区二区三区综合在线观看| 中文字幕久久专区| 亚洲精品粉嫩美女一区| 精品国内亚洲2022精品成人| 亚洲精品国产精品久久久不卡| 国产1区2区3区精品| 久久久久久久午夜电影| 国产亚洲精品综合一区在线观看 | 在线观看免费视频日本深夜| 这个男人来自地球电影免费观看| 国产视频一区二区在线看| 日韩有码中文字幕| avwww免费| 久久影院123| 天天添夜夜摸| 亚洲欧美激情在线| 国产亚洲欧美98| 日韩欧美一区二区三区在线观看| 午夜日韩欧美国产| 这个男人来自地球电影免费观看| 婷婷精品国产亚洲av在线| 国产单亲对白刺激| 精品人妻在线不人妻| 精品第一国产精品| 国产一区二区三区综合在线观看| 免费一级毛片在线播放高清视频 | 天堂动漫精品| 国产精品 欧美亚洲| 国产免费av片在线观看野外av| 夜夜躁狠狠躁天天躁| 欧美日韩中文字幕国产精品一区二区三区 | 国产国语露脸激情在线看| 嫩草影视91久久| 三级毛片av免费| 久久欧美精品欧美久久欧美| 亚洲一卡2卡3卡4卡5卡精品中文| 精品电影一区二区在线| 欧美色视频一区免费| 91字幕亚洲| 欧美黄色淫秽网站| 久久狼人影院| 亚洲人成伊人成综合网2020| 男女午夜视频在线观看| 亚洲avbb在线观看| 日韩欧美三级三区| 国产精品一区二区免费欧美| 亚洲无线在线观看| 中文字幕精品免费在线观看视频| 国产精品 国内视频| 精品一区二区三区视频在线观看免费| 禁无遮挡网站| 国产精品综合久久久久久久免费 | 欧美性长视频在线观看| 亚洲国产精品成人综合色| 夜夜看夜夜爽夜夜摸| 午夜精品久久久久久毛片777| 91字幕亚洲| 精品欧美一区二区三区在线| av视频免费观看在线观看| 国产精品电影一区二区三区| 黄色视频不卡| 国产成人精品久久二区二区免费| 黄频高清免费视频| 一边摸一边做爽爽视频免费| 黄色片一级片一级黄色片| 老熟妇仑乱视频hdxx| 男人的好看免费观看在线视频 | 精品久久久久久成人av| 国产一卡二卡三卡精品| 国产日韩一区二区三区精品不卡| 国产一区二区三区综合在线观看| 99久久99久久久精品蜜桃| 亚洲精品中文字幕一二三四区| 侵犯人妻中文字幕一二三四区| 久久久久久久久久久久大奶| 亚洲电影在线观看av| 日本一区二区免费在线视频| 久久国产乱子伦精品免费另类| 国产激情久久老熟女| 日韩欧美三级三区| 真人做人爱边吃奶动态| 国产高清有码在线观看视频 | 在线观看午夜福利视频| 国产精品免费视频内射| 男人操女人黄网站| 午夜免费成人在线视频| 9热在线视频观看99| 成人永久免费在线观看视频| 久久中文字幕一级| 亚洲精品一卡2卡三卡4卡5卡| 如日韩欧美国产精品一区二区三区| 黑人巨大精品欧美一区二区蜜桃| 国产97色在线日韩免费| 久久婷婷人人爽人人干人人爱 | 国产三级在线视频| 亚洲精品国产一区二区精华液| 精品乱码久久久久久99久播| 日日摸夜夜添夜夜添小说| 日本a在线网址| 中文亚洲av片在线观看爽| 精品国产一区二区久久| 一a级毛片在线观看| 亚洲欧美日韩高清在线视频| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲专区字幕在线| 免费看十八禁软件| 欧美一区二区精品小视频在线| 久久久久国产一级毛片高清牌|