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

    Green Aircraft Taxiing Strategy Based on Multi-scenario Joint Optimization

    2016-09-14 01:16:41WanLiliHuMinghuaTianYongYeBojia

    Wan Lili,Hu Minghua,Tian Yong,Ye Bojia

    College of Civil Aviation,Nanjing University of Aeronautics and Astronautics,Nanjing 211106,P.R.China(Received 5June 2015;revised 18August 2015;accepted 28August 2015)

    ?

    Green Aircraft Taxiing Strategy Based on Multi-scenario Joint Optimization

    Wan Lili,Hu Minghua*,Tian Yong,Ye Bojia

    College of Civil Aviation,Nanjing University of Aeronautics and Astronautics,Nanjing 211106,P.R.China
    (Received 5June 2015;revised 18August 2015;accepted 28August 2015)

    The issue of green aircraft taxiing under various taxi scenarios is studied to improve the efficiency of aircraft surface operations and reduce environmental pollution around the airport from aircraft emissions.A green aircraft taxi programming model based on multi-scenario joint optimization is built according to airport surface network topology modeling by analyzing the characteristics of aircraft operations under three different taxiing scenarios:all-engine taxi,single-engine taxi,and electronic taxi.A genetic algorithm is also used in the model to minimize fuel consumption and pollutant emissions.The Shanghai Pudong International Airport is selected as a typical example to conduct a verification analysis.Compared with actual operational data,the amount of aircraft fuel consumption and gas emissions after optimization are reduced significantly through applying the model.Under an electronic taxiing scenario,fuel consumption can be lowered by 45.3%,and hydrocarbon(HC)and carbon dioxide(CO)emissions are decreased by 80%.The results show that a green aircraft taxiing strategy that integrates taxiway optimization and electronic taxiing can effectively improve the efficiency of airport operations and reduce aircraft pollution levels in an airport′s peripheral environment.

    electronic taxi;fuel consumption;emission;genetic algorithm

    0 Introduction

    In the wake of rapid development in the air transportation industry,it is increasingly difficult to meet the ever growing demands of air traffic with the relatively limited supply of airport capacity.Capacity flow imbalances cause frequent airport surface congestion and flight delays,which cause a substantial increase in aircraft surface taxiing time.In addition,because the fuel efficiency of aircrafts is the lowest during the surface taxi phase,the gas emissions index is the subsequently highest,and the emissions directly affect the airport′s peripheral environment[1].Issues such as aircraft fuel consumption,emissions,and noise have received widespread concern within the aviation industry as well as from the public.Under a backdrop of ever deepening concepts in global en-vironmental protection,energy savings,and emissions reduction,methods of using highly effective scientific surface taxi strategies to maximize conflict avoidance in aircraft operations and reduce aircraft surface taxiing time have important implications for improving the efficiency of aircraft operations,lowering airline operating costs,and reducing pollution in the airport′s peripheral environment.

    Studies in China on airport surface taxi optimization have primarily concentrated on optimizing the taxiway.Under the premise of ensuring the safety of aircraft operation,taxiway optimization takes airport network topology modeling as backdrop,and utilizes mathematical programming models,multi-agent technology,Petri net theory,etc.to search for the optimal path for a departing scheduled flight that includes the shortesttaxiway or the least amount of taxi time[2-3]. However,aircraft all make use of main engine propulsion when taxiing,thus,the aircraft is in a state of low combustion efficiency and high emissions over the entire distance that it idles,and such operations have a relatively large impact on the airport′s peripheral environment.Current studies in China on surface taxi optimization have not yet fully considered the needs of green aviation development and rarely taken aircraft fuel consumption and emissions as optimization objectives.However,studies conducted in other countries related to environmentally friendly aircraft taxiing have set up multi-scenario taxi systems for all-engine taxi,single-engine taxi,and electronic taxi,etc.to maximize the reduction of work time during engine idling[4-7].The electronic taxi system,which uses an auxiliary power unit(APU),drives the electric motors that are installed on the landing gears,thus allowing the aircraft to taxi without being towed by a tug or propelled by the main engine.This method also has low fuel consumption,low emissions,and low noise levels,etc.Currently,companies such as Honeywell and WheelTug are developing this type of electronic taxi system,and certain products are already in an identification phase[8,9].

    Aircraft fuel consumption and emissions model constructed under different taxi scenarios is presented,and it integrates the taxiway optimization method,which is applied to the Pudong International Airport for verification and the comparative analysis of environmental pollution under various taxi scenarios.The work presented in this paper will enable air traffic control authorities to set and provide methods of supporting highly effective,and environmentally friendly green aircraft taxiing strategies.

    1 Multi-scenario Fuel Consumption and Emissions Models

    The aircraft taxi process includes taxi-out and taxi-in processes.Taxi-out is the ground movement process by which an aircraft moves from the gate to the head of the runway prior to takeoff.Taxi-in is the ground movement process by which an aircraft moves to a gate after landing and leaving the runway.When taxiing,an aircraft can be powered by the engine or APU.According to selected power,the scenarios may include all-engine taxi,single-engine taxi,and electronic taxi[10].

    1.1 Scenario one:All-engine taxi

    All-engine taxi refers to a process in which all of the engines are turned on and remain in an idling state during aircraft taxiing,and it is the most commonly adopted method in aircraft operations at present.In past studies,fuel calculations were simplified by assigning an average value for the engine fuel flow rate.However,to accurately calculate fuel consumption and emissions during taxiing,fluctuations in the engine state and fuel flow rate as the aircraft moves must be considered.Studies have clearly shown that stopping the engine to wait decreases the fuel flow rate and subsequent acceleration significantly increases the fuel flow rate compared with taxiing at an even speed[11].The fuel consumption FAiof aircraft i during all-engine taxi can be expressed as follows

    where Fijis the fuel consumption amount of aircraft i under engine state j,Tijthe taxi time for aircraft i under engine state j,Nithe number of engines on aircraft i,and fijthe fuel flow rate of a single engine on aircraft i under state j.

    The hydrocarbon(HC),carbon dioxide(CO),and mono-nitrogen oxides(NOx)[12]emitted by the aircraft will directly impact the local air quality index.After exposure to the sun′s ultraviolet rays,HC can generate toxic photochemical smog that severely irritates the eyes,nose,and throat,whereas CO is a colorless and odorless poisonous gas,and NOxcan easily form acid rain and damage the ozone layer.Therefore,the aforementioned three gaseous pollutants are usedto build an aircraft emissions model in this paper.

    The amount of gaseous pollutant emissions from an aircraft is directly proportional to the fuel consumption amount and emissions indices.However,each gaseous pollutant has a different emissions index.The emissions amount EAikof gaseous pollutant type k from aircraft i during all-engine taxi can be expressed as follows

    where Iijkis the emissions index of gaseous pollutant type k when aircraft iis under engine state j.

    1.2 Scenario two:Single-engine taxi

    Single-engine taxi is a taxiing method in which some of the engines are turned off during aircraft taxi while the rest are turned on and put in an idling state.During single-engine aircraft taxi,only the engines that are turned on will consume fuel and produce emissions.If single-engine taxi is used during aircraft taxi-out,the engines that are off will need to be started up prior to entering the runway for takeoff.The engine startup time(ESUT)typically takes 2—5min and is related to the model and engine model number of the aircraft and time of the engine has been shut off.If the engine has been shut off for six hours or more,ESUT will increase to 10—15min.During single-engine aircraft taxi-in,the engines that are not used for taxiing after landing must also be cooled down.The engine cool-down time(ECDT)is similar to ESUT.Because the engines continue to consume fuel and produce emissions during the start-up and cool-down processes,these two time intervals must be considered when calculating single-engine taxi fuel consumption and emissions.The single-engine taxi fuel consumption FSiof aircraft i can be expressed as follows

    where Tiis the taxi time for aircraft i,NSithe number of engines turned on by aircraft i during taxiing,TEthe time required for engine start-up or cool-down,and fidlethe fuel flow rate of a sin-gle engine on aircraft i during engine start-up or cool-down.

    The emissions amount ESikof gaseous pollutant type k from aircraft i during single-engine taxi can be expressed as follows

    1.3 Scenario three:Electronic taxi

    Electronic taxi is a process in which the engines are turned off during aircraft taxiing and only APU is used to drive the electric motors installed on the landing gears for taxiing.Thus,an electronic taxi system is used for taxiing.

    APU is an aircraft auxiliary power unit typically located at the tail end of the aircraft.Although its small turbine engine core consumes fuel during power output,the fuel consumption and emissions are far lower than that of the main engines.In addition,there are notable advantages in its use,such as low noise and extended engine life,etc.The fuel consumption FEiof aircraft i when using the electronic taxi system is as follows

    where fAPUiis the fuel flow rate of APU on aircraft i.

    The emissions amount EEikof gaseous pollutant type k from aircraft i during electronic taxi can be expressed as follows

    where IEikis the emissions index of gaseous pollutant type k when aircraft i is using the electronic taxi system to taxi.

    Similar to single-engine taxi,when an aircraft uses the electronic taxi system to taxi,all engines that are off prior to entering the runway must be started for takeoff and all engines must be cooled off after landing.If these two segments are executed independently,the taxi time is increased and may even lead to extra delays.To shorten taxi time,APU-driven taxiing can be integrated with the engine start-up and cool-down segments.Taking a twin-engine aircraft as an example,the process of using APU-driven electronic taxiing is indicated in Fig.1[13].

    Fig.1 Process diagram of use of APU-driven electronic taxiing

    During the taxi-out process,aircraft engines can provide power to the remainder of ESUT prior to takeoff.Otherwise,the aircraft is driven by APU during taxi time.During the taxi-in process,the aircraft changes to APU-driven taxiing only after ECDT following landing.The fuel consumption of aircraft i when using the electronic taxi system is as follows

    where fidleiis the fuel flow rate of aircraft iin an idling state.

    2 Taxiway Optimization

    2.1 Taxiway optimization model

    To obtain the optimal green taxi strategy,the taxi method must be integrated with route optimization.When building a taxiway optimization model,the topological network diagram G=(R,E)can be used to indicate the airport surface taxi system,where Rrepresents the set of taxiway nodes r and Erepresents the edge set of taxiway intersections[14].The constructed taxiway optimization model is presented as follows

    where m∈Mis the set of runways and p∈Pthe set of gates.

    Eq.(8)takes the minimum total fuel consumption of all aircraft as an objective function. Each aircraft can select Eqs.(1),(3),(7)to calculate the fuel consumption Ficorresponding to the different taxi scenarios.Eq.(9)is the set of taxiway nodes Uiof aircraft i,where uih∈R,(uih,ui,h+1)∈E.Eq.(10)is the total taxi time Tiof aircraft i,where Tiuihis the taxi time of aircraft i on taxi segment(uih,ui,h+1).In Eq.(11),|uihui,h+1|is the length of taxi segment(uih,ui,h+1).In Eq.(12),vihis the average taxiing velocity of aircraft i on taxi segment(uih,ui,h+1),and vihis the velocity of aircraft i when passing taxiway node uih.In Eqs.(13),(14),λimand gipare the decision variables for runway and gate. Eqs.(15),(16)are the unique constraints in runway and gate assignment.In addition,taxiing aircraft must satisfy relevant airport operating rules and safety intervals.

    2.2 Optimization algorithm

    When an aircraft encounters another aircraft during the taxi process,situations such as deadlock conflict,intersection conflict,and one aircraft overtaking another may occur.When optimizing aircraft taxiway,all types of conflicts must be considered and solutions must be provided so that conflicts between aircrafts do not occur and each aircraft can continue to taxi along the relief path to a predetermined position.In order toimprove the overall efficiency of the operation of the taxiway and meet other special requirements,the speed adjustment strategy is adopted according to the flight priority.It can ensure that the aircraft which priority is higher or who first reach the conflict position preferentially pass.The taxiway optimization model is a non-deterministic polynomial(NP)problem and difficult to find the global optimal solution by algorithm precise in acceptable time.Genetic algorithms are classic intelligent optimization methods that simulate biological survival of the fittest evolutionary processes,and they are frequently applied in the field of route optimization because of their high efficiency in search for the global optimal solution.The specific steps of the genetic algorithm in this paper are listed below.

    Step 1 Obtain flight plan data,including flight arrival and departure times,information of gates and runways.

    Step 2 Build a topological network diagram of the airport surface taxiway system and establish relationship matrices of taxiway node sets,taxiway edge sets,and vertex edges.

    Step 3 Generate the initial parent population according to the flight information of gate and runway and randomly assign taxiways.

    Step 4 Perform the genetic manipulation,such as selection,crossover,and mutation,etc.,on the inital parent population according to the vertex edge relationship matrix,then form a progeny population.

    Step 5 Calculate the time at which a flight passes each taxiway node according to the flight arrival or departure time and perform conflict detection and resolution.

    Step 6 Determine whether the current evolutionary process satisfies the termination condition.If so,end the algorithm;Otherwise,return to Step 5and continue the iteration.

    3 Case Study

    The flight plan data from July 25,2013at Shanghai Pudong International Airport is selected for a case study.Shanghai Pudong International Airport is a famous hub airport with intensive air traffic flow and complex surface taxiing that can effectively support the verification analysis of the method presented in this paper.On the day in question,Pudong International Airport used three runways and 193gates and ran a total of 1 089flights,of which 545are arrivals and 544 are departures.A topological structure is constructed(See Fig.2)according to the structure of the Pudong International Airport surface taxi system for calculating fuel consumption and emissions while integrating each taxi scenario with the taxiway optimization model and comparing simulated data with actual operating data.

    Fig.2 Surface topological map of Shanghai Pudong International Airport

    3.1 Fuel consumption

    The fuel consumption by an aircraft can be calculated according to the fuel flow rate of all aircraft models under various states,which is provided by the International Civil Aviation Organization(ICAO)Engine Emissions Databank[15].The total amount of fuel consumption by all aircraft that taxied on the surface at Pudong International Airport can be minimized through simulated optimizations.Fig.3illustrates genetic algorithm optimization process of Scenario 1.

    The population optimization curve within 100 generations shows that the total amount of fuel consumption rapidly decreases from 237 071.8kg to 172 764.9kg,whereas between 100to 400 generations,the total amount of fuel consumption changes less rapidly,falling to 156 187.5kg,andafter 400generations,the changes are small and gradually converge after 450generations.In the optimal solution corresponding to the maximum evolution generation,the total fuel consumption is 154 087.9kg.

    Fig.3 Trends of fuel consumption in population evolution process

    Similarly,optimization can be performed on Scenarios 2and 3.A comparison of the fuel consumption prior to taxiway optimization and after optimization under different taxi scenarios is shown in Fig.4.The total fuel consumption prior to taxiway optimization is calculated based on the surface surveillance data.

    Fig.4 Comparison of fuel consumption pre-optimization and post-optimization under different taxi scenarios

    Fig.4shows that using the same all-engine taxi scenario,taxiway optimization saves 13.5% in fuel compared with the amount pre-optimization.Using single-engine scenario and electronic taxi senario,the fuel consumption can be further reduced by 34.9%and 45.3%.

    In addition,the fuel savings of different aircraft models during the optimization process are also different.Fig.5represents the average fuel consumption of different aircraft models at Pudong International Airport.Because the fuel flow rate in the main engines of heavy aircraft is grea-ter during all-engine and single-engine taxi,the fuel savings after using the electronic taxi system is greater compared with that of mid-sized aircraft,such as aircraft models A340,A380,B747,MD11,etc.Because the fuel flow rate in the main engines of aircraft models E90and CRJ is smaller,the amount of fuel consumption will present a greater relative increase during electronic taxi.

    Fig.5 Average fuel consumed by each aircraft model under different scenarios

    3.2 Taxi time

    After optimizing the aircraft surface taxiway,fuel consumption is reduced and flight taxi time is optimized.Fig.6shows the taxi time distribution of each aircraft model after optimization under Scenario 1,where pre-optimal taxi time(PET)represents the average pre-optimization taxi time and post-optimal taxi time(POT)represents the average taxi time after optimization.The average aircraft taxi time after optimization is within a 10-to 20-minute range and represents a reduction of approximately 7%compared with pre-optimization amounts.

    Fig.6 Comparison of taxi time for scenario one before and after optimization

    Based on ESUT before takeoff and ECDT after landing in Scenarios 2and 3,the taxi time after optimization is slightly longer than that of Scenario 1.However,compared with pre-optimization time,there is still a reduction in average taxi time for all other aircraft models,as shown in Fig.7,with the exception of A300and A310 aircraft models,which have a small optimization opportunity because of the small number of flights.Thus,Increased taxi time is presented.

    Fig.7 Average taxi time of each aircraft model in different taxi scenarios

    3.3 Gaseous pollutant emissions

    The amount of gaseous pollutant emissions from each aircraft model operating on the surface of Pudong International Airport can be calculated according to the HC,CO,and NOxgaseous pollutant emission indices for each aircraft model provided by the ICAO Engine Emissions Databank.The calculation results are shown in Fig.8.

    Fig.8 Amount of pollutant emissions at Shanghai Pudong International Airport by different taxi methods

    The amounts of HC and CO gaseous pollutant emissions clearly decrease after surface optimization,with approximately 13%and 35%in Scenarios 1and 2,respectively,and steep decrease of 80%in Scenario 3.The average airport surface taxi time of A320and B737series aircraft is 10—15min,which provides sufficient time for engine start-up and taxiing and reduces the a-mount of gaseous pollutant emissions by approximately 30%.Because a three-engine or four-engine aircraft is driven by more than one engine,the reduction in gaseous pollutant emissions is relatively smaller in proportion,with a decrease of 10%—15%,as shown in Figs.9,10.

    The reduction in emissions proportion of NOxgaseous pollutant emissions in Scenarios 1 and 2after optimization is similar to that of HC and CO(See Fig.11).However,the higher NOxemissions index during APU operation causes an increase in NOxgaseous pollutant emissions in Scenario 3by approximately 6%compared with the amount pre-optimization and approximately 22.64%and 62.29%compared with the amounts in Scenarios 1and 2,respectively.

    Fig.9 Amount of HC emissions of each aircraft model in different taxi scenarios

    Fig.10 Amount of CO emissions of each aircraft model in different taxi scenarios

    Fig.11 Amount of NOxemissions of each aircraft model in different taxi scenarios

    3.4 Electronic taxi system with fuel addition

    Although an electronic taxi system can sig-nificantly reduce fuel consumption and emissions during the surface operations phase,the drive system that must be installed on the landing gears will increase the basic operating empty weight(OEW)and the extra fuel consumption for the aircraft in other flight phases.Assuming an electronic taxi system makes OEW average increase in 400kg weight,the amount of extra fuel consumption of each aircraft model caused by the increase in OEW must be calculated according to the aircraft performance indices provided by the BADA3.11database.The results are shown in Fig.12.

    Fig.12 Changes in aircraft fuel consumption when using an electronic taxi system

    The amount of fuel saved when using an electronic taxi system is approximately 45%of the taxi fuel,and the increased fuel consumption in the cruise phase is approximately 16%of the taxi fuel.Thus,the proportion of fuel saved when using APU to taxi is far greater than the increased fuel consumption in the cruise phase,and the aircraft′s total fuel consumption decreases. Although the increment of fuel consumption caused by increase of OEW is larger,the fuel consumption saved in taxi phase for a heavy aircraft is relatively more.Therefore,the use of an electronic taxi system in a heavy aircraft can conserve more fuel over the entire flight period.

    4 Conclusions

    A green aircraft taxi strategy based on multiscenario joint optimization is proposed by analyzing aircraft fuel consumption and exhaust emissions under different taxi scenarios,such as allengine taxi,single-engine taxi,and electronic taxi,and by integrating these scenarios with the taxiway optimization model to improve the efficiency of flight operations,reduce delays,and decrease negative impacts on the environment.The example of Shanghai Pudong International Airport shows that the amount of aircraft fuel consumption and gas emissions after taxiway optimization can be lowered by about 13%,combined with environmentally friendly technology,such as single-engine taxi and electronic taxi,the performance of airport surface operations can be further improved,F(xiàn)or example,electronic taxi can reduce fuel consumption by 45.3%,reduce HC and CO emissions by 80%,and increase NOxemissions by 6%.

    References:

    [1] PENNER J,LISTER D,GRIGGS D,et al.Aviation and the global atmosphere:A special report of the intergovernmental panel on climate change[M].Cambridge UK:Cambridge University Press,1999:29-165.

    [2] ZHU Xinping,TANG Xinmin,HAN Songchen. Aircraft initial taxiing route planning based on petri net and genetic algorithm[J].Journal of Southwest Jiaotong University,2013,48(3):565-573.(in Chinese)

    [3] TANG Yong,HU Minghua,HUANG Rongshun,et al.Aircraft taxi routes planning based on free time windows and multi-agent for ASMGCS[J].Acta Aeronautica et Astronautica Sinica,2015,36(5):1627-1638.(in Chinese)

    [4] NHUT J P,REN Liling.Continuous descent approach:Design and flight test for louisville international airport[J].Journal of Aircraft,2008,41(5):1054-1066.

    [5] TODD J,CHRISTOPHER C,MICHAEL K,et al. Investigating the impact of off-nominal events on high-density″Green″arrivals[C]∥30th Digital Avionics Systems Conference.Seattle:IEEE,2011:2D5-1-2D5-16.

    [6] DEONANDAN I,BALAKRISHNAN H.Evaluation of strategies for reducing taxi-out emissions at airports[C]∥Proceedings of the AIAA Aviation Technology,Integration,and Operations(ATIO)Conference.Texas:AIAA,2010:13-15.

    [7] SIMAIAKIS I,BALAKRISHNAN H,KHADILKAR H,et al.Demonstration of reduced airport congestion through pushback rate control[J].Transportation Research Part A,2014,66:251-267.

    [8] WIJNTERP C,ROLING P C,de WILDE W.Electric taxi systems:An operations and value estimation[C]∥14th AIAA Aviation Technology,Integration,and Operations Conference.Atlanta:AIAA,2014:1-16.

    [9] RE F,de CASTRO R.Energetically optimal path following for electric aircraft taxi systems based on convex optimization[C]∥ESARS 2014International Electric Vehicle Conference.Italy:IEEE,2014:16-18.

    [10]GUO Rui,ZHANG Yu,WANG Qing.Comparison of emerging ground propulsion systems for electrified aircraft taxi operations[J].Transportation Research Part C,2014,44:98-109.

    [11]TASOS N,GAUTAM G,MATTHEW K,et al. Detailed estimation of fuel consumption and emissions during aircraft taxi operations at Dallas/Fort Worth International Airport[J].Transportation Research Part D,2011,16(4):302-308.

    [12]SUN Jianzhong,ZUO Hongfu,LIU Pengpeng,et al.Aircraft engine pollutant emissions estimation methods[J].Journal of Traffic and Transportation Engineering,2012,12(2):53-61.(in Chinese)

    [13]ROLING P,SILLEKENS P,CURRAN R.The effects of electric taxi systems on airport surface congestion[C]∥Proceedings of the AIAA Aviation Technology,Integration,and Operations(ATIO)Conference.Texas:AIAA,2015:1-10.

    [14]TIAN Yong,F(xiàn)U Jianjun,WANG Yanjun.Research on airport ground capacity evaluation[J].Journal of Nanjing University of Aeronautics &Astronautics,2006,38(5):619-622.(in Chinese)

    [15]Eurocontol Experimental Center.Model accuracy summary report for the base of aircraft data(BADA):Revision 3.11,EEC Technical/Scientific Report No.12/04/10-48[R].Bretigny:EEC,2012.

    Ms.Wan Lili is currently a lecturer in College of Civil Aviation at Nanjing University of Aeronautics and Astronautics(NUAA).She received her M.S.degree at NUAA.Her research interest is green air traffic management.

    Mr.Hu Minghua is currently aprofessor and doctoral supervisor at NUAA.His research interests are air traffic flow management and airspace management.

    Dr.Tian Yong is currently an associate professor at NUAA.He received his first degree and Ph.D.degree at NUAA.His research interests are air traffic planning and performance analysis.

    Dr.Ye Bojiais currently a lecturer at NUAA.He received his first degree and Ph.D.degree at NUAA.His research interests are collaborative air traffic management and system modeling and simulation.

    (Executive Editor:Xu Chengting)

    V351.11 Document code:A Article ID:1005-1120(2016)02-0215-09

    *Corresponding author,E-mail address:minghuahu@nuaa.edu.cn.

    How to cite this article:Wan Lili,Hu Minghua,Tian Yong,et al.Green aircraft taxiing strategy based on multi-scenario joint optimization[J].Trans.Nanjing Univ.Aero.Astro.,2016,33(2):215-223.

    http://dx.doi.org/10.16356/j.1005-1120.2016.02.215

    中文亚洲av片在线观看爽| 男女之事视频高清在线观看| 国内精品美女久久久久久| xxx96com| 女警被强在线播放| 国产伦人伦偷精品视频| 老熟妇仑乱视频hdxx| 免费看十八禁软件| 黄片小视频在线播放| 国内少妇人妻偷人精品xxx网站 | 国产av麻豆久久久久久久| 久久久色成人| 九九热线精品视视频播放| 91av网一区二区| 精品国产超薄肉色丝袜足j| 国产高清三级在线| av女优亚洲男人天堂 | 99久久综合精品五月天人人| cao死你这个sao货| 亚洲av美国av| av天堂在线播放| 麻豆成人av在线观看| 亚洲在线观看片| 午夜精品久久久久久毛片777| 中文字幕人妻丝袜一区二区| 2021天堂中文幕一二区在线观| 不卡av一区二区三区| 级片在线观看| 天天躁日日操中文字幕| 美女高潮喷水抽搐中文字幕| 日韩精品青青久久久久久| 91在线观看av| 亚洲欧美精品综合久久99| 男人的好看免费观看在线视频| 久久久久久国产a免费观看| 午夜a级毛片| 1000部很黄的大片| 亚洲专区字幕在线| 国产成人啪精品午夜网站| 国产亚洲精品久久久com| 国产精品精品国产色婷婷| 国内揄拍国产精品人妻在线| 久久精品91蜜桃| 麻豆一二三区av精品| 精品久久蜜臀av无| 麻豆国产av国片精品| 久久中文字幕一级| 午夜福利18| 99久国产av精品| 十八禁人妻一区二区| 亚洲狠狠婷婷综合久久图片| 国产蜜桃级精品一区二区三区| 亚洲中文字幕一区二区三区有码在线看 | tocl精华| 99久久无色码亚洲精品果冻| 亚洲精品粉嫩美女一区| 亚洲乱码一区二区免费版| 国产精品1区2区在线观看.| 看片在线看免费视频| 久久伊人香网站| 国产成人一区二区三区免费视频网站| 亚洲九九香蕉| 久久天堂一区二区三区四区| 女人高潮潮喷娇喘18禁视频| 精品国产亚洲在线| 欧美不卡视频在线免费观看| 亚洲 欧美一区二区三区| 久久草成人影院| tocl精华| 午夜福利在线观看吧| 黄色视频,在线免费观看| 国产伦精品一区二区三区视频9 | 久久久久国产一级毛片高清牌| 一个人观看的视频www高清免费观看 | 亚洲国产精品sss在线观看| 精品国产超薄肉色丝袜足j| 91在线观看av| 亚洲av成人精品一区久久| 精品久久久久久成人av| 午夜激情欧美在线| 国产精品 国内视频| 国产一区二区在线av高清观看| 啦啦啦免费观看视频1| www.自偷自拍.com| 成人高潮视频无遮挡免费网站| 日韩大尺度精品在线看网址| 国产视频内射| 国产亚洲精品久久久com| 亚洲aⅴ乱码一区二区在线播放| 变态另类丝袜制服| 国产91精品成人一区二区三区| 国产高潮美女av| 亚洲欧美日韩东京热| 一级毛片高清免费大全| 黄色丝袜av网址大全| 国产综合懂色| 免费搜索国产男女视频| 欧美性猛交黑人性爽| 亚洲av美国av| 国产蜜桃级精品一区二区三区| 国内少妇人妻偷人精品xxx网站 | 色综合站精品国产| 男女做爰动态图高潮gif福利片| ponron亚洲| 久久精品91蜜桃| 成年版毛片免费区| 国产亚洲精品av在线| 99热这里只有是精品50| 色尼玛亚洲综合影院| www国产在线视频色| 亚洲乱码一区二区免费版| 亚洲av五月六月丁香网| 日韩欧美在线乱码| 成人欧美大片| 脱女人内裤的视频| 高清在线国产一区| 欧美日韩亚洲国产一区二区在线观看| 麻豆国产av国片精品| www.www免费av| 18禁美女被吸乳视频| 天天添夜夜摸| 视频区欧美日本亚洲| 亚洲av成人av| 岛国视频午夜一区免费看| 国产精品一及| 日本 av在线| 97人妻精品一区二区三区麻豆| 大型黄色视频在线免费观看| 九色国产91popny在线| 亚洲第一欧美日韩一区二区三区| 1024手机看黄色片| 亚洲成av人片免费观看| 桃红色精品国产亚洲av| 国产精品久久久久久人妻精品电影| 亚洲av日韩精品久久久久久密| 法律面前人人平等表现在哪些方面| 欧美日韩福利视频一区二区| 国模一区二区三区四区视频 | 国产高清videossex| 欧美中文综合在线视频| 老司机福利观看| 99久久综合精品五月天人人| 日本撒尿小便嘘嘘汇集6| 非洲黑人性xxxx精品又粗又长| 国产淫片久久久久久久久 | 亚洲人成伊人成综合网2020| 精品免费久久久久久久清纯| 听说在线观看完整版免费高清| 亚洲av成人av| 五月玫瑰六月丁香| 国产91精品成人一区二区三区| 午夜福利欧美成人| 日韩成人在线观看一区二区三区| 久久久久久人人人人人| 午夜视频精品福利| 欧美日韩黄片免| 小说图片视频综合网站| 九九在线视频观看精品| 一个人看的www免费观看视频| 国产精品亚洲美女久久久| 国产免费男女视频| 国产伦一二天堂av在线观看| 看黄色毛片网站| 在线国产一区二区在线| 人人妻人人看人人澡| 一级毛片高清免费大全| 欧美日本视频| 9191精品国产免费久久| 亚洲中文日韩欧美视频| 日韩欧美在线乱码| 国产熟女xx| 91在线观看av| ponron亚洲| 日韩免费av在线播放| 一区二区三区高清视频在线| 欧美一区二区精品小视频在线| 精品国产超薄肉色丝袜足j| 亚洲中文字幕一区二区三区有码在线看 | 亚洲精品一卡2卡三卡4卡5卡| 少妇熟女aⅴ在线视频| 不卡av一区二区三区| 99国产精品一区二区三区| 亚洲国产精品合色在线| 国产毛片a区久久久久| 曰老女人黄片| 中文亚洲av片在线观看爽| 国产亚洲av高清不卡| 欧美三级亚洲精品| 一进一出抽搐gif免费好疼| 国产视频一区二区在线看| 夜夜看夜夜爽夜夜摸| av在线蜜桃| 成人国产一区最新在线观看| 国产野战对白在线观看| 无遮挡黄片免费观看| 日韩欧美精品v在线| 搡老岳熟女国产| 国产一区在线观看成人免费| 日本 av在线| 国产精品一及| 亚洲精品国产精品久久久不卡| 欧美黄色片欧美黄色片| 国产高清视频在线观看网站| 99在线人妻在线中文字幕| 一卡2卡三卡四卡精品乱码亚洲| 首页视频小说图片口味搜索| 久久久久国产一级毛片高清牌| 欧美丝袜亚洲另类 | а√天堂www在线а√下载| 国产黄色小视频在线观看| 最近最新免费中文字幕在线| 一本久久中文字幕| 黄色片一级片一级黄色片| 很黄的视频免费| 久久久国产成人精品二区| 一卡2卡三卡四卡精品乱码亚洲| 色综合欧美亚洲国产小说| 99久久成人亚洲精品观看| 真实男女啪啪啪动态图| 成人一区二区视频在线观看| 小蜜桃在线观看免费完整版高清| 成人精品一区二区免费| 亚洲av成人一区二区三| 欧美激情久久久久久爽电影| 黄色成人免费大全| 国产精品综合久久久久久久免费| 日本黄大片高清| 日本五十路高清| 18禁黄网站禁片免费观看直播| 久久亚洲精品不卡| 亚洲av免费在线观看| 天天添夜夜摸| 老鸭窝网址在线观看| 后天国语完整版免费观看| 免费在线观看亚洲国产| 色av中文字幕| 国产高清三级在线| 亚洲国产欧美网| 国语自产精品视频在线第100页| av天堂中文字幕网| 日韩免费av在线播放| 欧美av亚洲av综合av国产av| 日本 av在线| 老汉色av国产亚洲站长工具| 极品教师在线免费播放| 国产成人精品无人区| av女优亚洲男人天堂 | 国产亚洲精品一区二区www| 婷婷亚洲欧美| 老熟妇仑乱视频hdxx| 丁香欧美五月| 88av欧美| 精品国产乱码久久久久久男人| 中文字幕熟女人妻在线| 岛国在线免费视频观看| 国产淫片久久久久久久久 | www.精华液| 亚洲欧美日韩高清在线视频| 熟妇人妻久久中文字幕3abv| 国产视频一区二区在线看| 国内少妇人妻偷人精品xxx网站 | 免费看美女性在线毛片视频| 一个人观看的视频www高清免费观看 | 久久久国产精品麻豆| 搡老妇女老女人老熟妇| 午夜精品在线福利| 久久精品国产99精品国产亚洲性色| 性欧美人与动物交配| 午夜精品一区二区三区免费看| bbb黄色大片| 国产成人av教育| 国产av一区在线观看免费| 看黄色毛片网站| 亚洲人成网站在线播放欧美日韩| 国产99白浆流出| 少妇人妻一区二区三区视频| 美女高潮喷水抽搐中文字幕| 黄频高清免费视频| 日本一二三区视频观看| 成年免费大片在线观看| 国产高清三级在线| 久久精品亚洲精品国产色婷小说| 一个人看的www免费观看视频| 亚洲一区二区三区色噜噜| 欧美av亚洲av综合av国产av| 久久中文看片网| 91在线精品国自产拍蜜月 | 国产一区二区在线av高清观看| 国产精品女同一区二区软件 | 一a级毛片在线观看| 亚洲国产欧洲综合997久久,| 免费搜索国产男女视频| 午夜激情欧美在线| 中文字幕最新亚洲高清| 超碰成人久久| 老熟妇仑乱视频hdxx| 国产伦精品一区二区三区视频9 | 国产毛片a区久久久久| 亚洲专区中文字幕在线| 日韩欧美一区二区三区在线观看| 成年女人永久免费观看视频| 久久精品91无色码中文字幕| 色视频www国产| 脱女人内裤的视频| 国产黄片美女视频| 久久人人精品亚洲av| 国产欧美日韩精品亚洲av| 嫩草影院入口| 国产精品亚洲av一区麻豆| 怎么达到女性高潮| 国产精品一区二区精品视频观看| 欧美一区二区精品小视频在线| 久久精品夜夜夜夜夜久久蜜豆| 久久久久久国产a免费观看| 性色avwww在线观看| 国产v大片淫在线免费观看| 99国产精品99久久久久| 一二三四在线观看免费中文在| 国产高潮美女av| www.999成人在线观看| 听说在线观看完整版免费高清| 嫩草影院入口| 国产成人欧美在线观看| 真人做人爱边吃奶动态| h日本视频在线播放| 日韩成人在线观看一区二区三区| 亚洲av免费在线观看| 中亚洲国语对白在线视频| www国产在线视频色| 日韩欧美在线乱码| 成人av一区二区三区在线看| 99久久精品热视频| 成年女人永久免费观看视频| a级毛片a级免费在线| 亚洲精华国产精华精| a级毛片在线看网站| 国产v大片淫在线免费观看| 精品久久久久久久人妻蜜臀av| 国产精品一区二区精品视频观看| 欧美又色又爽又黄视频| 男女床上黄色一级片免费看| 国产精品久久视频播放| 亚洲av电影不卡..在线观看| 一本综合久久免费| 国产精品久久久久久亚洲av鲁大| 国产91精品成人一区二区三区| 日韩欧美免费精品| 美女 人体艺术 gogo| 日日夜夜操网爽| 成人特级黄色片久久久久久久| 91av网一区二区| 我要搜黄色片| 亚洲成a人片在线一区二区| 波多野结衣巨乳人妻| 成人av一区二区三区在线看| 国产精品,欧美在线| 欧美xxxx黑人xx丫x性爽| www.精华液| 国产精品国产高清国产av| 成人三级黄色视频| 舔av片在线| aaaaa片日本免费| 真人做人爱边吃奶动态| 人妻夜夜爽99麻豆av| 久久久精品欧美日韩精品| 精品一区二区三区四区五区乱码| 国产探花在线观看一区二区| 亚洲国产欧美网| 99视频精品全部免费 在线 | 久久精品91无色码中文字幕| 亚洲avbb在线观看| 香蕉久久夜色| 色在线成人网| 国产精品av久久久久免费| 欧美日韩一级在线毛片| 亚洲av成人一区二区三| 午夜福利高清视频| 亚洲精品色激情综合| 毛片女人毛片| 中文在线观看免费www的网站| 在线观看免费视频日本深夜| 国产不卡一卡二| 69av精品久久久久久| 久久性视频一级片| 中文字幕最新亚洲高清| 悠悠久久av| 岛国视频午夜一区免费看| 亚洲国产色片| 欧美最黄视频在线播放免费| 国产午夜福利久久久久久| 久久香蕉国产精品| 日本一二三区视频观看| 亚洲第一电影网av| 一二三四在线观看免费中文在| av中文乱码字幕在线| 精品免费久久久久久久清纯| 国产成人欧美在线观看| 亚洲国产日韩欧美精品在线观看 | 久久精品91蜜桃| 美女高潮喷水抽搐中文字幕| 一边摸一边抽搐一进一小说| svipshipincom国产片| 亚洲午夜精品一区,二区,三区| 国产不卡一卡二| 观看美女的网站| 亚洲第一电影网av| 国产亚洲av高清不卡| 一进一出好大好爽视频| 欧美色视频一区免费| 人妻丰满熟妇av一区二区三区| 12—13女人毛片做爰片一| 国产黄a三级三级三级人| 国产淫片久久久久久久久 | 精品福利观看| 美女高潮喷水抽搐中文字幕| 一本久久中文字幕| 精品国产乱子伦一区二区三区| 国产一级毛片七仙女欲春2| 亚洲av免费在线观看| 日韩有码中文字幕| 可以在线观看毛片的网站| 大型黄色视频在线免费观看| 搡老妇女老女人老熟妇| 午夜福利成人在线免费观看| 国产99白浆流出| 国产欧美日韩精品亚洲av| 97超级碰碰碰精品色视频在线观看| 久久久久久久久中文| 18禁黄网站禁片免费观看直播| 亚洲aⅴ乱码一区二区在线播放| 国产单亲对白刺激| 国产毛片a区久久久久| 在线免费观看的www视频| 久久午夜亚洲精品久久| 搡老熟女国产l中国老女人| 久久人人精品亚洲av| 性色avwww在线观看| 久久久久久久久免费视频了| ponron亚洲| 亚洲av第一区精品v没综合| 成熟少妇高潮喷水视频| 三级国产精品欧美在线观看 | 亚洲国产精品久久男人天堂| 变态另类成人亚洲欧美熟女| 麻豆国产97在线/欧美| 亚洲成a人片在线一区二区| 老司机深夜福利视频在线观看| 久久这里只有精品中国| 免费大片18禁| 亚洲第一电影网av| 亚洲 国产 在线| 国产一区二区三区在线臀色熟女| 人妻久久中文字幕网| 国产精品98久久久久久宅男小说| 亚洲欧美精品综合一区二区三区| 亚洲国产精品sss在线观看| 青草久久国产| 国产欧美日韩精品一区二区| 国内揄拍国产精品人妻在线| 青草久久国产| 特级一级黄色大片| 色综合亚洲欧美另类图片| 国产欧美日韩一区二区精品| 日韩欧美在线乱码| 嫩草影视91久久| 性色avwww在线观看| 国产精品电影一区二区三区| 男人舔女人下体高潮全视频| 在线观看午夜福利视频| 亚洲成av人片在线播放无| 国产一区在线观看成人免费| 国产亚洲精品久久久com| 亚洲无线观看免费| 久久天堂一区二区三区四区| 在线观看日韩欧美| 老司机在亚洲福利影院| 欧美黑人欧美精品刺激| 亚洲欧美精品综合久久99| 麻豆成人av在线观看| 香蕉丝袜av| 国产av一区在线观看免费| 免费看a级黄色片| 嫩草影院入口| 国产av麻豆久久久久久久| 精品国产三级普通话版| 午夜精品在线福利| 久久久久国产精品人妻aⅴ院| 少妇丰满av| 在线免费观看不下载黄p国产 | svipshipincom国产片| 国产三级中文精品| 免费人成视频x8x8入口观看| 成人欧美大片| 婷婷六月久久综合丁香| 99精品久久久久人妻精品| 免费av不卡在线播放| 床上黄色一级片| 级片在线观看| 最好的美女福利视频网| 亚洲五月婷婷丁香| 国产亚洲精品一区二区www| 长腿黑丝高跟| 国产v大片淫在线免费观看| 999久久久国产精品视频| 久久亚洲真实| 日本五十路高清| 国产成人影院久久av| 一本综合久久免费| 亚洲av成人av| 久久伊人香网站| 午夜影院日韩av| 最近在线观看免费完整版| 欧美日韩瑟瑟在线播放| 90打野战视频偷拍视频| svipshipincom国产片| 精品久久久久久,| 青草久久国产| 亚洲精品一卡2卡三卡4卡5卡| 岛国视频午夜一区免费看| 亚洲avbb在线观看| 久久久久久九九精品二区国产| 99久久无色码亚洲精品果冻| 日韩大尺度精品在线看网址| 国产免费av片在线观看野外av| 欧美日韩福利视频一区二区| 亚洲午夜精品一区,二区,三区| 成年女人永久免费观看视频| 成人亚洲精品av一区二区| 午夜影院日韩av| 午夜精品在线福利| 九九在线视频观看精品| 99国产综合亚洲精品| 国产成人精品久久二区二区免费| 啦啦啦免费观看视频1| 欧美色欧美亚洲另类二区| 全区人妻精品视频| 久久天躁狠狠躁夜夜2o2o| 亚洲av成人一区二区三| 19禁男女啪啪无遮挡网站| 欧美日韩亚洲国产一区二区在线观看| 麻豆久久精品国产亚洲av| 亚洲精品久久国产高清桃花| 老司机深夜福利视频在线观看| 国产精品一区二区三区四区久久| 国产成人欧美在线观看| 亚洲 国产 在线| 嫁个100分男人电影在线观看| 免费在线观看视频国产中文字幕亚洲| 宅男免费午夜| 国产精品永久免费网站| 日韩成人在线观看一区二区三区| 午夜精品在线福利| 黄色日韩在线| 母亲3免费完整高清在线观看| 美女被艹到高潮喷水动态| АⅤ资源中文在线天堂| 日本在线视频免费播放| 一区福利在线观看| 午夜激情欧美在线| 男人舔女人的私密视频| 日韩高清综合在线| 一个人看的www免费观看视频| 国产亚洲欧美在线一区二区| 99热精品在线国产| 亚洲激情在线av| 中文在线观看免费www的网站| 亚洲精品一区av在线观看| 久久久久性生活片| 精品国产亚洲在线| 香蕉国产在线看| 婷婷精品国产亚洲av| 中亚洲国语对白在线视频| 亚洲自偷自拍图片 自拍| 制服丝袜大香蕉在线| 99国产精品99久久久久| 成年女人永久免费观看视频| 日本一本二区三区精品| 深夜精品福利| 在线播放国产精品三级| 久久久成人免费电影| 人妻久久中文字幕网| 精华霜和精华液先用哪个| 夜夜夜夜夜久久久久| 最好的美女福利视频网| 在线看三级毛片| 伊人久久大香线蕉亚洲五| bbb黄色大片| 日本免费一区二区三区高清不卡| 免费搜索国产男女视频| 曰老女人黄片| 国产成人精品无人区| h日本视频在线播放| 亚洲欧美一区二区三区黑人| 欧美日韩瑟瑟在线播放| 香蕉丝袜av| 久久久久国产精品人妻aⅴ院| or卡值多少钱| 国产精品亚洲美女久久久| 天天添夜夜摸| 一本综合久久免费| 香蕉丝袜av| 亚洲无线观看免费| 欧美日韩乱码在线| 亚洲无线观看免费| 狠狠狠狠99中文字幕| 日韩国内少妇激情av| 亚洲国产精品成人综合色| 国内精品一区二区在线观看| 天天躁日日操中文字幕| 国产毛片a区久久久久| 午夜福利在线观看吧| 在线观看66精品国产| 一本久久中文字幕| 日本黄色视频三级网站网址| 给我免费播放毛片高清在线观看|