• <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

    国产v大片淫在线免费观看| 一区福利在线观看| 日韩欧美精品v在线| 男女下面进入的视频免费午夜| 精品国产美女av久久久久小说| 美女黄网站色视频| 十八禁网站免费在线| 亚洲无线在线观看| 免费大片18禁| 又大又爽又粗| 在线观看日韩欧美| 亚洲精品456在线播放app | 一本精品99久久精品77| 国内精品久久久久精免费| 91在线观看av| 精品福利观看| 97碰自拍视频| 一进一出好大好爽视频| 国产成人福利小说| 欧美日韩乱码在线| www.www免费av| 一进一出抽搐动态| 校园春色视频在线观看| 国产成人av激情在线播放| 好看av亚洲va欧美ⅴa在| 久久国产精品影院| 999久久久国产精品视频| www日本黄色视频网| 国产精品免费一区二区三区在线| 国产精品爽爽va在线观看网站| 色播亚洲综合网| 日日夜夜操网爽| 狂野欧美激情性xxxx| 国产高清三级在线| 国产精品99久久久久久久久| 黑人巨大精品欧美一区二区mp4| 好男人电影高清在线观看| 亚洲国产看品久久| 老汉色av国产亚洲站长工具| 精品免费久久久久久久清纯| 亚洲欧美日韩高清专用| 国产单亲对白刺激| 男女那种视频在线观看| 搞女人的毛片| 999精品在线视频| 国产又色又爽无遮挡免费看| 亚洲精品美女久久av网站| 成人午夜高清在线视频| 久久久色成人| 成人三级做爰电影| 成人欧美大片| 国产毛片a区久久久久| 亚洲 国产 在线| 国产主播在线观看一区二区| 一个人免费在线观看的高清视频| 中国美女看黄片| 三级毛片av免费| 男女那种视频在线观看| 国产毛片a区久久久久| 又粗又爽又猛毛片免费看| 国产精品av视频在线免费观看| 老汉色av国产亚洲站长工具| 亚洲九九香蕉| 免费人成视频x8x8入口观看| 久久午夜亚洲精品久久| 免费在线观看成人毛片| 制服人妻中文乱码| 亚洲国产日韩欧美精品在线观看 | 成人国产综合亚洲| 亚洲成人免费电影在线观看| 99国产极品粉嫩在线观看| 久久婷婷人人爽人人干人人爱| 麻豆国产97在线/欧美| 欧美在线一区亚洲| 久久久色成人| 亚洲av电影不卡..在线观看| 看免费av毛片| 老汉色∧v一级毛片| www.999成人在线观看| 亚洲aⅴ乱码一区二区在线播放| 中亚洲国语对白在线视频| 99久国产av精品| 国产成人av激情在线播放| 国产精品久久电影中文字幕| 很黄的视频免费| av黄色大香蕉| 欧美中文综合在线视频| 国产淫片久久久久久久久 | 黄色成人免费大全| 日韩人妻高清精品专区| 国产精品一及| 精品国产乱码久久久久久男人| 麻豆国产av国片精品| 国产人伦9x9x在线观看| 亚洲最大成人中文| 午夜两性在线视频| 美女大奶头视频| 国产欧美日韩精品亚洲av| 香蕉丝袜av| 欧美3d第一页| 脱女人内裤的视频| 日本成人三级电影网站| 亚洲欧美日韩无卡精品| 亚洲av美国av| 亚洲av熟女| 久久中文字幕一级| 女生性感内裤真人,穿戴方法视频| 最近最新免费中文字幕在线| 桃色一区二区三区在线观看| 看免费av毛片| 日日夜夜操网爽| 免费在线观看日本一区| 午夜两性在线视频| 国产精品,欧美在线| 国产精品国产高清国产av| 国产伦精品一区二区三区四那| 99国产综合亚洲精品| 变态另类成人亚洲欧美熟女| 18美女黄网站色大片免费观看| 91av网一区二区| 国产极品精品免费视频能看的| 国产又色又爽无遮挡免费看| a级毛片a级免费在线| 桃红色精品国产亚洲av| 99久久综合精品五月天人人| 亚洲av中文字字幕乱码综合| 国产精品亚洲一级av第二区| 亚洲国产看品久久| 黄色视频,在线免费观看| 黄片大片在线免费观看| 欧美另类亚洲清纯唯美| 亚洲成av人片在线播放无| 国产99白浆流出| 国产精品香港三级国产av潘金莲| 亚洲狠狠婷婷综合久久图片| 久久久久久久精品吃奶| 欧美激情在线99| cao死你这个sao货| av女优亚洲男人天堂 | 久久草成人影院| 亚洲人成电影免费在线| 黄色片一级片一级黄色片| 日本一本二区三区精品| 日日夜夜操网爽| 日韩中文字幕欧美一区二区| 精品一区二区三区视频在线 | 中文字幕最新亚洲高清| 国内精品一区二区在线观看| 少妇的逼水好多| 午夜视频精品福利| 精品国产三级普通话版| 亚洲欧美日韩东京热| 亚洲成人精品中文字幕电影| 欧美最黄视频在线播放免费| 亚洲熟女毛片儿| 国产精品亚洲一级av第二区| 男女做爰动态图高潮gif福利片| 久久久国产精品麻豆| 人人妻人人看人人澡| 国内精品美女久久久久久| 婷婷亚洲欧美| 亚洲真实伦在线观看| 人人妻人人澡欧美一区二区| 九色成人免费人妻av| 亚洲aⅴ乱码一区二区在线播放| 别揉我奶头~嗯~啊~动态视频| 亚洲成a人片在线一区二区| 亚洲va日本ⅴa欧美va伊人久久| 国产精品影院久久| 国产精品爽爽va在线观看网站| 黄色片一级片一级黄色片| 叶爱在线成人免费视频播放| 国产亚洲欧美98| 午夜日韩欧美国产| 免费看光身美女| 99精品在免费线老司机午夜| 国产熟女xx| www.999成人在线观看| 日本熟妇午夜| 欧美成人免费av一区二区三区| 精品一区二区三区av网在线观看| 动漫黄色视频在线观看| 在线免费观看的www视频| 十八禁人妻一区二区| 黄色丝袜av网址大全| 搞女人的毛片| 村上凉子中文字幕在线| 精品国产乱码久久久久久男人| 成人鲁丝片一二三区免费| 久久九九热精品免费| 又爽又黄无遮挡网站| 可以在线观看的亚洲视频| 成人一区二区视频在线观看| 99国产精品一区二区三区| 免费看a级黄色片| av国产免费在线观看| 亚洲九九香蕉| 每晚都被弄得嗷嗷叫到高潮| 国产真人三级小视频在线观看| 在线a可以看的网站| 国产私拍福利视频在线观看| 久久天堂一区二区三区四区| 久久精品91无色码中文字幕| 国产成人精品久久二区二区91| 国产视频内射| 亚洲av中文字字幕乱码综合| 国产精品一区二区免费欧美| 精品不卡国产一区二区三区| 国产精品99久久99久久久不卡| 色在线成人网| 一进一出好大好爽视频| 精品乱码久久久久久99久播| 久久久色成人| 这个男人来自地球电影免费观看| 丰满人妻熟妇乱又伦精品不卡| 久久草成人影院| 精品久久久久久久人妻蜜臀av| 精品久久久久久成人av| 丰满人妻熟妇乱又伦精品不卡| 人人妻人人看人人澡| 精品福利观看| 欧美日韩黄片免| 亚洲熟妇熟女久久| 操出白浆在线播放| 首页视频小说图片口味搜索| 高潮久久久久久久久久久不卡| 免费无遮挡裸体视频| 人人妻,人人澡人人爽秒播| 亚洲一区二区三区色噜噜| 色综合站精品国产| 日本一二三区视频观看| 久久九九热精品免费| 丁香六月欧美| 三级国产精品欧美在线观看 | 国产成人精品无人区| 亚洲七黄色美女视频| 久9热在线精品视频| 日韩有码中文字幕| 国产成人精品久久二区二区免费| 1000部很黄的大片| 亚洲精品美女久久久久99蜜臀| 九九热线精品视视频播放| 十八禁人妻一区二区| 亚洲av成人精品一区久久| 国内精品美女久久久久久| 国产一区在线观看成人免费| 国模一区二区三区四区视频 | 国产熟女xx| 日本熟妇午夜| 久久欧美精品欧美久久欧美| 久久99热这里只有精品18| 人人妻人人澡欧美一区二区| 一个人免费在线观看的高清视频| 我要搜黄色片| 久久中文字幕一级| 色av中文字幕| 中文字幕最新亚洲高清| 国产伦一二天堂av在线观看| 国产精品久久久久久亚洲av鲁大| 三级男女做爰猛烈吃奶摸视频| 中文字幕人妻丝袜一区二区| 噜噜噜噜噜久久久久久91| 在线观看免费午夜福利视频| 亚洲男人的天堂狠狠| 波多野结衣巨乳人妻| 大型黄色视频在线免费观看| 亚洲成人中文字幕在线播放| cao死你这个sao货| 欧美成人一区二区免费高清观看 | 18禁黄网站禁片免费观看直播| 精华霜和精华液先用哪个| 99riav亚洲国产免费| 日韩大尺度精品在线看网址| 欧美高清成人免费视频www| 人人妻人人澡欧美一区二区| 久久天堂一区二区三区四区| 嫩草影院入口| 成人一区二区视频在线观看| 国产精品一区二区三区四区免费观看 | 国产亚洲欧美在线一区二区| 日本黄色片子视频| 免费在线观看视频国产中文字幕亚洲| 欧美激情久久久久久爽电影| 人人妻,人人澡人人爽秒播| 琪琪午夜伦伦电影理论片6080| 日韩av在线大香蕉| 最近最新中文字幕大全电影3| 日日夜夜操网爽| 在线观看66精品国产| 免费看光身美女| 亚洲美女黄片视频| 999久久久国产精品视频| 黄频高清免费视频| 99热6这里只有精品| 国产伦在线观看视频一区| 丁香欧美五月| 久久久成人免费电影| 麻豆国产av国片精品| 亚洲中文日韩欧美视频| 国产伦在线观看视频一区| 草草在线视频免费看| 欧美日韩黄片免| 舔av片在线| 欧美日韩亚洲国产一区二区在线观看| 老司机深夜福利视频在线观看| 在线观看一区二区三区| 久久这里只有精品中国| 69av精品久久久久久| 亚洲五月婷婷丁香| 可以在线观看毛片的网站| 精华霜和精华液先用哪个| 亚洲欧美日韩卡通动漫| 成人18禁在线播放| 曰老女人黄片| 久久香蕉精品热| 亚洲狠狠婷婷综合久久图片| 99久久综合精品五月天人人| 欧美三级亚洲精品| 亚洲精品国产精品久久久不卡| 视频区欧美日本亚洲| 精品一区二区三区视频在线 | 青草久久国产| 天堂影院成人在线观看| av福利片在线观看| 亚洲av成人不卡在线观看播放网| 日本a在线网址| 亚洲精品乱码久久久v下载方式 | 成年女人看的毛片在线观看| 99久久精品热视频| tocl精华| 欧美绝顶高潮抽搐喷水| 成人高潮视频无遮挡免费网站| 男女床上黄色一级片免费看| www.www免费av| 狂野欧美激情性xxxx| 日本黄色视频三级网站网址| 国产精品自产拍在线观看55亚洲| 亚洲国产中文字幕在线视频| 中文亚洲av片在线观看爽| 91麻豆精品激情在线观看国产| 十八禁人妻一区二区| 国产探花在线观看一区二区| 99热这里只有精品一区 | 久久婷婷人人爽人人干人人爱| 99国产精品一区二区蜜桃av| 九色国产91popny在线| 国产精品电影一区二区三区| 99视频精品全部免费 在线 | av女优亚洲男人天堂 | 淫秽高清视频在线观看| 亚洲国产精品久久男人天堂| 亚洲一区二区三区色噜噜| 国内少妇人妻偷人精品xxx网站 | 老熟妇乱子伦视频在线观看| 看片在线看免费视频| 亚洲成人中文字幕在线播放| 黄频高清免费视频| 国产精品久久电影中文字幕| 日本成人三级电影网站| 国产极品精品免费视频能看的| 国产免费av片在线观看野外av| 久9热在线精品视频| 国产乱人伦免费视频| 国产黄a三级三级三级人| 欧美黄色淫秽网站| 免费在线观看日本一区| 一a级毛片在线观看| 在线观看日韩欧美| 日本撒尿小便嘘嘘汇集6| 欧美日本视频| 国产精品久久久av美女十八| 99riav亚洲国产免费| 亚洲专区字幕在线| 99在线人妻在线中文字幕| 亚洲精华国产精华精| 一本综合久久免费| 免费电影在线观看免费观看| 亚洲狠狠婷婷综合久久图片| bbb黄色大片| 天堂av国产一区二区熟女人妻| 在线永久观看黄色视频| 亚洲一区高清亚洲精品| 又黄又粗又硬又大视频| 此物有八面人人有两片| 嫁个100分男人电影在线观看| 真实男女啪啪啪动态图| 此物有八面人人有两片| 国产亚洲av嫩草精品影院| 日韩精品中文字幕看吧| 午夜视频精品福利| 高清毛片免费观看视频网站| 亚洲欧美日韩高清在线视频| 小说图片视频综合网站| 久久久久久久午夜电影| 男人的好看免费观看在线视频| 国产成人精品久久二区二区免费| 99re在线观看精品视频| 欧美色视频一区免费| 精品久久久久久久久久免费视频| 精品一区二区三区视频在线观看免费| 91av网站免费观看| tocl精华| 免费无遮挡裸体视频| 亚洲国产看品久久| 非洲黑人性xxxx精品又粗又长| 白带黄色成豆腐渣| 美女黄网站色视频| 亚洲人成电影免费在线| 亚洲va日本ⅴa欧美va伊人久久| 19禁男女啪啪无遮挡网站| 免费在线观看亚洲国产| 成人高潮视频无遮挡免费网站| 久久精品91蜜桃| 亚洲av电影在线进入| 国产v大片淫在线免费观看| 国产伦一二天堂av在线观看| 亚洲av日韩精品久久久久久密| 国产 一区 欧美 日韩| 欧美一区二区精品小视频在线| av视频在线观看入口| 亚洲av中文字字幕乱码综合| 熟女电影av网| a级毛片a级免费在线| 国内精品久久久久精免费| 黄色 视频免费看| 亚洲五月天丁香| 久久热在线av| 久久精品夜夜夜夜夜久久蜜豆| 日韩高清综合在线| 国产精华一区二区三区| 亚洲国产精品合色在线| 久久久久亚洲av毛片大全| 真人一进一出gif抽搐免费| 午夜成年电影在线免费观看| av国产免费在线观看| 精品免费久久久久久久清纯| 亚洲国产欧美网| 久久天堂一区二区三区四区| 欧美3d第一页| 高潮久久久久久久久久久不卡| 亚洲中文字幕一区二区三区有码在线看 | x7x7x7水蜜桃| 亚洲av中文字字幕乱码综合| 色在线成人网| 岛国在线免费视频观看| 香蕉久久夜色| 国产欧美日韩精品亚洲av| 中文字幕最新亚洲高清| 色噜噜av男人的天堂激情| bbb黄色大片| 香蕉久久夜色| 亚洲中文av在线| 最近最新中文字幕大全电影3| 亚洲成人久久性| tocl精华| 九九在线视频观看精品| 丁香欧美五月| 国产97色在线日韩免费| 精品久久久久久久人妻蜜臀av| 国产av不卡久久| 香蕉丝袜av| 91九色精品人成在线观看| 精品一区二区三区视频在线观看免费| 久久久久国内视频| 国产成+人综合+亚洲专区| 一二三四社区在线视频社区8| 制服人妻中文乱码| 免费av毛片视频| 亚洲男人的天堂狠狠| 9191精品国产免费久久| 变态另类丝袜制服| 免费在线观看日本一区| 两个人的视频大全免费| 99在线视频只有这里精品首页| www日本在线高清视频| 免费看美女性在线毛片视频| 久久伊人香网站| 成人三级黄色视频| 久久久国产成人精品二区| 成人三级黄色视频| 欧美成人性av电影在线观看| 搡老熟女国产l中国老女人| 国产精品av久久久久免费| 给我免费播放毛片高清在线观看| 午夜福利在线观看免费完整高清在 | 黄片小视频在线播放| 日韩欧美在线二视频| 久久精品aⅴ一区二区三区四区| 久久久久国产一级毛片高清牌| 亚洲欧美精品综合一区二区三区| 天天躁日日操中文字幕| 亚洲国产欧美网| 久久久久久国产a免费观看| 麻豆一二三区av精品| 国产精品一区二区三区四区久久| 亚洲aⅴ乱码一区二区在线播放| 亚洲国产日韩欧美精品在线观看 | 99久久无色码亚洲精品果冻| 又黄又粗又硬又大视频| 国产免费av片在线观看野外av| 国产高清videossex| 欧美午夜高清在线| 日本黄色片子视频| 91av网一区二区| 色综合欧美亚洲国产小说| www.自偷自拍.com| 午夜福利免费观看在线| 怎么达到女性高潮| 国产精品香港三级国产av潘金莲| 欧美成人一区二区免费高清观看 | 制服丝袜大香蕉在线| 中文字幕av在线有码专区| 国产高潮美女av| 国产成人av教育| www.自偷自拍.com| 色精品久久人妻99蜜桃| 久久香蕉精品热| 桃红色精品国产亚洲av| 两性午夜刺激爽爽歪歪视频在线观看| 国产成人精品无人区| 成年版毛片免费区| or卡值多少钱| 99国产精品一区二区三区| av在线天堂中文字幕| 午夜免费激情av| 三级国产精品欧美在线观看 | 亚洲成av人片在线播放无| 长腿黑丝高跟| 国产精品电影一区二区三区| 天天躁狠狠躁夜夜躁狠狠躁| www国产在线视频色| 免费av不卡在线播放| 免费搜索国产男女视频| 麻豆成人av在线观看| 久久午夜亚洲精品久久| 免费观看精品视频网站| 国产一区在线观看成人免费| 欧美日韩亚洲国产一区二区在线观看| 男人舔女人的私密视频| 久久久色成人| 国产亚洲精品av在线| 日本a在线网址| 亚洲国产日韩欧美精品在线观看 | 一区福利在线观看| 午夜成年电影在线免费观看| 男人舔奶头视频| 国产成人欧美在线观看| 欧美日本亚洲视频在线播放| 国产亚洲精品久久久com| 欧美一区二区精品小视频在线| 很黄的视频免费| 在线视频色国产色| 午夜福利成人在线免费观看| 18禁国产床啪视频网站| 成人一区二区视频在线观看| or卡值多少钱| 老司机午夜福利在线观看视频| 99久久精品热视频| 夜夜躁狠狠躁天天躁| 久久精品影院6| 国产探花在线观看一区二区| 婷婷六月久久综合丁香| xxx96com| 波多野结衣高清无吗| 久久草成人影院| 亚洲成人精品中文字幕电影| 国产久久久一区二区三区| 一区二区三区激情视频| 在线观看免费午夜福利视频| 人人妻,人人澡人人爽秒播| 日本黄色片子视频| 日本黄大片高清| 免费在线观看亚洲国产| 亚洲 国产 在线| 欧美日韩福利视频一区二区| 久久久久亚洲av毛片大全| 国产视频内射| 国产又色又爽无遮挡免费看| 老司机午夜福利在线观看视频| 亚洲欧美日韩无卡精品| 欧美极品一区二区三区四区| 啦啦啦免费观看视频1| 午夜久久久久精精品| 久久精品国产清高在天天线| 日日摸夜夜添夜夜添小说| 国产一区二区在线观看日韩 | 免费看十八禁软件| 精品电影一区二区在线| 亚洲成人久久爱视频| 国产精品久久久av美女十八| 麻豆av在线久日| 午夜视频精品福利| 91麻豆av在线| 久久久久久久久久黄片| 最新美女视频免费是黄的| 国产精品爽爽va在线观看网站| 亚洲专区字幕在线| 这个男人来自地球电影免费观看| 亚洲国产精品合色在线| 国产精品久久久久久人妻精品电影| 身体一侧抽搐| 色视频www国产| 久久天堂一区二区三区四区| 真人做人爱边吃奶动态| 美女高潮喷水抽搐中文字幕| 国产熟女xx| 国产精品精品国产色婷婷| 成年版毛片免费区| 淫秽高清视频在线观看| 老司机在亚洲福利影院| 久久久国产成人免费| 悠悠久久av| av女优亚洲男人天堂 | 国产淫片久久久久久久久 | 国产私拍福利视频在线观看|