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

    Analysis of Typical Earth-Mars Launch Trajectory on Mission Profiles and Launch Window Extension

    2018-03-21 08:11:22GengGuangyouWangJueSongQiangZhangZhiguoWangJianming
    關(guān)鍵詞:北京航空航天大學(xué)宇航力學(xué)

    Geng Guang-you , Wang Jue, Song Qiang, Zhang Zhi-guo, Wang Jian-ming

    ?

    Analysis of Typical Earth-Mars Launch Trajectory on Mission Profiles and Launch Window Extension

    Geng Guang-you1,2, Wang Jue1,3, Song Qiang2, Zhang Zhi-guo2, Wang Jian-ming2

    (1. School of Astronautics, Beihang University, Beijing, 100191; 2. Beijing Institute of Astronautical Systems Engineering, Beijing, 100076; 3.China Academy of Launch Vehicle Technology, Beijing, 100076)

    Taking into account characteristics and range safety of different launch vehicles implemented in recent years, there are obvious differences in launch trajectories of recent typical Mars exploration missions. Through the principle of flight mechanics on direct Earth-Mars transfer trajectory, characteristics of recent Mars exploration departure trajectories are analyzed and summarized, mainly on mission profiles and how to extend the launch window.

    Earth-Mars launch trajectory; Transfer trajectory; Launch window; Hyperbolic escape trajectory

    0 Introduction

    The average distance between Mars and the Sun is 1.524Au[1]. The sidereal period of revolution of Mars is 687 days. The period of conjunction with the Earth is 780 days. The included angle between Mars’ orbit and the ecliptic is about 1.85o. The eccentricity of Mars’ orbit is about 0.093 (the orbital eccentricity of the Earth is about 0.017). Different Earth-Mars relative positions caused by different launch years, lead to the complicated design of Earth-Mars transfer orbit, which mainly includes the heliocentric elliptical arc and the hyperbolic transfer trajectories at both ends (Earth or Mars). The transfer orbits characterized by non-coplanar, small eccentricity and multi-perturbation have been studied for a long time as a nonlinear or typical multi-extremum problem. However, the success rate of Mars exploration is only near 50% so far, including flyby, orbiting and landing (about 22successes).

    In addition, the latitudes of most launch complex are different apparently; their scopes of launch azimuth and range safety also constrain the trajectory launched. So in recent years, the departure hyperbolic trajectories toward Mars are very different, here we’ll analyze and summarize why the characteristics of recent Mars exploration launch trajectories are so different.

    1 Earth-Mars Transfer Flight Mechanics Model

    Ignoring the influence of other perturbation forces and the flight phase inside the gravitational sphere of influence (SOI), the Earth-Mars transfer orbit is a typical Lambert's problem. The mathematic expression of the two-point transfer orbit is as follows[1~4]:

    Where μ is the gravitational constant of the centric celestial body; a is the semi-major axis of the transfer orbit; r1, r2 and c are distances between the departure point, the arrival point and the centric celestial body, respectively; the angle between r1 and r2 is θ, Fig.1 is illustration.

    From the north pole of Ecliptic, see Fig.2.

    Fig.2 Earth-Mars Transfer Orbit Illustration

    In general, the exact mechanics model of the Earth-Mars transfer orbit involves the gravity (the whole process) of the Sun, the Moon and eight planets,the dominant2terms of the Earth and Mars and the higher order non-spherical terms, the latter are not so important and can be omitted, as well as the impact of light pressure terms; without loss of generality, here the mechanics model is given to consider the Earth and Mars2coefficients:

    Based on the above analysis, it can be seen that for the Earth-Mars transfer orbit of the certain year, the DLA of the hyperbolic trajectory at the gravity boundary of SOI affects the launch vehicle's azimuth (viz. orbital inclination); the velocity change ΔVPE+ΔVPM affects the departure and arrival energy, and right ascension angle (RLA) affects the launch time. Although launch trajectories of Mars exploration are different due to different years, the characteristics of the departure and arrival energy and the DLA are similar. Therefore, in this case, for the Earth-Mars transfer orbit with a transfer time about 5~10 months launched in 2018, DLA and velocity change ΔVPE+ΔVPM analysis results are presented as Fig.3 and Fig.4. To be convenient, here we assume ΔVPEmeans the velocity increased from the departing 200 km altitude parking circle orbit to the escape hyperbolic trajectory of the Earth, andΔVPM means the velocity decreased from the reentry hyperbolic trajectory to the 500 km altitude circle orbit of Mars.

    Fig.4 DLA &RLA Varies with Launch Date andTransfer Time in 2018

    Then we find the optimal window of the total velocity change is from 2458240 JD (corresponding to May 1, 2018) to 2458270 JD, i.e. the optimal launch period lasts about 30 days, and with the transfer time of about 200 days, that means it’s a type I transfer.

    Fig.4 indicates that the DLA and RLA are determined by the departure and arrival date, thus the launch vehicle’s trajectory inclination (viz. azimuth demand) and launch time are determined as well, if we suppose the departure date and transfer days of the trajectory.

    Through the above, we could get the departure hyperbolic trajectory toward Mars, the below shows the typical launch trajectories for Mars exploration.

    2 Typical Launch Trajectory for Mars Exploration

    At present, most Martian probes’ orbits are polar or with large obliquity. Typical trajectory profiles launched and patterns on how to extend launch window are analyzed below.

    2.1 ExoMars / Schiaparelli Mission Launched by ESA/Russia[5]

    This is a Mars exploration projected by ESA and Russian Space Agency, with ExoMars as the Trace Gas Orbiter (TGO) and Schiaparelli as the Mars Lander. The total mass is 4332 kg. Finally, ExoMars Orbiter successfully circled around Mars, but Schiaparelli failed during landing, hitting the surface of Mars directly.

    The mission is launched instantaneously by the storable propellant rocket, Proton.The launch period is from 7thto 27thJanuary or from 14thto 25thMarch 2016. Breeze M is a good choice to perform deep space exploration because it works very efficiently due to its specific impulse up to 325 s, compact structure design and annular throwaway propellant tank. Located at the latitude of 45.6o, Baikonur launch site cannot resort the Earth's rotation speed much, but with the high DLA ability attained, thereby making optimal launch window possible.

    This mission is Proton’s first successful deep space launch in recent 20 years. Taking 10 hours and 44 minutes from the liftoff to the separation of the probe, the rocket's first three stage flight profiles are similar to common GTO launches. After that, Breeze M ignited four times, transferring from the three parking orbits to the hyperbolic escape trajectory toward Mars, the three parking orbits were 51.55°×175 km circle, 51.58°×292 km×5272 km, 51.55°×693 km×21079 km. After releasing the probe, Breeze M entered into the abandoned orbit to avoid the interference.

    The actual launching date is 14thMarch 2016 and the launch azimuth is 61.3o. The corresponding trajectory inclination is 51.5o and the probe arrived at Mars about 7 months later. On 16thOctober, Schiaparelli is separated from the orbiter at 5833.3 m/s speed and entered the atmosphere of Mars. Due to failure of the parachute braking device, the lander crashed. On 19thOctober 15:24 UTC, TGO orbiter entered the elliptical orbit, and later, finally circled around Mars with the orbit altitude of 400 km on March 2018.

    2.2 Mars Exploration Missions Launched by Atlas V of U.S.[6]

    In recent years, Atlas V series launched several Mars probes. Atlas V541 launched MAVEN in 2013, Curiosity in 2011, and reconnaissance orbiter in 2005. Curiosity even made a successful landing with parachutes. The other two were launched by Atlas V401 configuration (both probes entered large elliptical orbits of Mars). Atlas V401 is a two stage rocket using liquid oxygen–kerosene at first stage, and the second stage, Centaur with LO2/LH2. Atlas V541 added four solid boosters.

    As a flagship launch vehicle of United States, it has abundant launch capability. The launch period is generally about 20 days, with the daily launch window varies from 30 minutes to 2 hours. According to the available information, to ensure the launch window, the coast time varies dynamically according to the change of the liftoff time. As the super-cryogenic propellant upper stage could ignite more than 3 times within 6 hours, it can be speculated that a fixed azimuth is selected according to the trajectory inclination required by deep space exploration, depending on the liftoff time, a dog-leg maneuver with the loss of launch capability to some extent would be taken. The data show that 94o or 104o azimuths are used for different missions respectively, with corresponding orbital inclination of 28.7o or 35.5o (the latter with a primary yaw program firstly). Escape declination (DLA), right ascension angle (RLA), and launch energy C3 are taken as target parameters for Earth escape hyperbolic trajectory design.

    NASA's MAVEN mission launched in 2013 is mainly used to study the evolution of the Martian atmosphere and water presence. The launch period lasted 20 days, from 18thNovember to 7thDecember, 2013, see Tab.1. While the daily window lasted up to 2 hours; on special situations, the launch period can be postponed until 23rdDecember (up to a maximum of 36 days ). Main launch period restricted the DLA within 28° (so there is no extra loss of launch capability caused by increasing the azimuth). The Earth's hyperbolic escape trajectory was achieved by the fixed azimuth and dynamical adjustment of the coasting time. Centaur upper stage would optimize the escape trajectory phase, if the liftoff time changed.

    Tab.1 Earth Escape C3, RLA, DLA, Arrival Date of Mars

    Launch Date18th Nov.Launch Period Open27th Nov.LP Medium7th Dec.LP Day 20 C3/(km2·s-2)12.0710.309.40 DLA/(o)13.3818.8426.93 RLA/(o)198.26200.58200.88 Arrival Date2014-09-242014-09-242014-09-24

    On 18thNovember 2013 18:28 UTC, at the beginning of the 2 hours launch window, Atlas V 401 adopted a 94o azimuth (corresponding to 28.6° orbital inclination), lifted a 2454 kg Mars rover MAVEN into the escape trajectory from Cape Canaveral Air Force Base. During the actual launch sequence, coasting time between Centaur upper stage’s two ignitions was 1656 s (with the rolling thermal control program during the coasting phase), with 570 s and 329 s for two burns respectively. The perigee of transfer trajectory was 192 km. The entire flight took 308 days, with 4 orbital corrections. MAVEN arrived at Mars on 02:24 UTC, 22thSeptember 2014 (the close of the arrival window was next 28thSeptember, corresponding to the launch window closed on 7thDecember). The 1321 N engine (main engine 6×200N) consumed more than half of the fuel in more than 36 minutes, generating a velocity of 1233 m/s and entering a large elliptical orbit with periapsis of 380 km, orbital inclination of 75o and period of 35 hours. Later it maneuvered into a 150 km×6275 km, 4.5 hours period elliptical orbit. One-way speed-of-light time from Mars to Earth is 4~20 minutes which depends on the Earth's relative position.

    2.3 Mars Exploration Missions Launched by Delta II of U.S.[7]

    Delta II series launched a lot of Mars probes, including Mars Global Surveyor and Mars Pathfinder respectively in 1996, Mars Climate Orbiter in 1998, Mars Polar Lander in 1999, Mars Odyssey in 2001, MER-A Spirit and MER-B Opportunity respectively in 2003, Dawn and Phoenix Mars Lander respectively in 2007.

    Among missions launched in recent years, Delta 7925H launched Dawn in 2007 and MER-B Opportunity in 2003, while 7925 launched the remaining missions. Delta II 7925H used nine solid rocket boosters and first stage of LO2/Kerosene, second stage of N2O4/Hydrazine50, and the solid upper stage. 7925H employed longer solid boosters and larger fairings than 7925. Due to its limited launch capability, Delta II series have a typical launch window of about 19 minutes to 44 minutes per day for a launch period of 20 days. According to the available information, Delta II series adopt two fixed azimuth corresponding to the two lift off time to achieve the launch; for example, Phoenix launched in 2007, MER-A Spirit and MER-B Opportunity launched in 2003, all clearly adopted the method.

    The data show that 93o and 99o azimuths are used for different missions, respectively, see Tab.2. The corresponding orbital angles are 28.6o and 29.8o, respectively. The following is a description of the Phoenix Mars mission launched in 2007 by Delta II.

    Phoenix probe weighs 680 kg. The main task is to explore water ice and early life on the surface of Mars. Its launch period is 3rdAugust 2007.

    On 3rdAugust launch window, the first launch opportunity was 5:35:18 EDT (US Eastern Time); the second launch opportunity is 6:11:24 EDT (Eastern US time), due to the relative position of celestial bodies, et al; the next window is generally about 10 minutes earlier than the previous day. For orbital parameters of 3rdAugust, at separation of the probe, the altitude is 2255 km and the velocity is 11 021.5 m/s according to the 93o azimuth.

    Tab.2 Planned Launch Sequence of the Phoenix Timing

    EventFlight Time/(min: sec) 93°Azimuth99° Azimuth Liftoff00:00.000:00.0 Stage I-II Separation04:31.304:31.3 First Cutoff – Stage II (SECO-1)09:20.509:21.0 First Restart – Stage II73:47.271:51.3 Second Cutoff – Stage II (SECO-2)76:02.374:06.3 Target Interface Point (TIP)87:42.885:46.9

    2.4 MOM Probe Launched by Indian PSLV-XL[8]

    Indian PSLV is a four and a half stage-rocket (boosters, core 1ststage and 3rdstage using solid propulsion systems, core 2ndstage and 4thstage using liquid propulsion systems). PSLV injected the probe into the sub-GTO orbit firstly, then multiple perigee burns were performed to compensate the liftoff time offsets, which was almost the same idea used in United States’ first lunar orbit mission or Chang'e I, China's first lunar exploration in 2007. It was originally planned to launch on 28thOctober, 2013, but was delayed for about 1 week. The whole flight profile is as follows.

    On 5thNovember, 2013, at 14:38 (09:08 UT), Indian PSLV-XL launch vehicle lifted off. It sent the Mars probe to an elliptical orbit with orbital inclination of 19.27o (about 104o launch azimuth), perigee of 250 km, and apogee of 23 500 km. Then after igniting 6 times, the probe is sent into a hyperbolic transfer orbit. The fourth maneuver only raised the apogee to 78 276 km (originally planned 100 000 km) due to the power solenoid conflict, and two days later an additional maneuver raised the apogee to 118 642 km. About 300 days later, the probe arrived at Mars on 24thSeptember, 2014. Its 440 N main engine (and 8×22N reaction control system) ignited for 1454 seconds and consumed 249.5 kg propellant, decelerated 1098.7 m/s into an elliptical orbit with periapsis of 423 km and apoapsis of 80 000 km to Mars.

    2.5 Rosetta Probe Launched by Ariane V[9]

    The purpose of Rosetta probe is to detect the comet 67P/Churyumov-Gerasimenko. The first obstacle to be overcome is the limited launch window. In order to meet up with Comet Wirtanen, Rosetta must be launched within a period of about 20 days, starting on 26thFebruary 2004 (postponed from the January 2003). On six of those days within the period, there must be a launch window of just 20~30 minutes. The remaining days account for a roll back of the launch vehicle for replenishing of the cryogenic fuel. If the opportunity is missed, this mission will have to be postponed while another target is selected.

    After flying by the Earth 3 times and Mars once, the probe entered into the orbit heading to the asteroid and comet, which is Ariane's first deep space exploration launched in recent years.

    The Ariane V launched at 07:17 UTC on 2ndMarch 2004 from Kourou, which is near the equator, most likely used a “fixed azimuth + dog-leg maneuver + coasting time adjustment” method to achieve the launch window. The flight lasted 1 hour and 45 minutes from liftoff to the separation of probe which weighs 3187 kg. When solid boosters stopped working, the onboard computer began to optimize the target of the parking orbit and the escape hyperbolic trajectory.

    The hyperbolic escape trajectory parameters are as follows:∞= 3.545 km/s, DLA= 2o. As the launch is delayed, the actual parking orbit and target orbital parameters are adjusted; in this situation, the first stage entered into a 200 km×4000 km parking orbit, and about 2 hours later, the upper stage entered into the hyperbolic trajectory toward Mars.

    3 Characteristics of Typical Launch Trajectory Toward Mars in Recent Years

    Relevant information for launch vehicles of U.S., ESA, etc shows that, the key points of the Earth escape trajectory towards Mars targets are: RLA, DLA, and escape C3 parameters. Based on the above analysis of typical Mars exploration trajectory launched in recent years, and coupled with the mechanics of flight, the following conclusions are drawn:

    a) Since Russia’s launch complex is located at high latitude, and Proton (including Breeze M) with storable propellants, among its approximate 20 or 11 days launch period, Proton launch instantaneously. Multi-perigee burns are performed by Breeze M to achieve Mars exploration hyperbolic escape orbit in 10 more hours.

    b) For early Delta II series, among its approximate 20 days launch period, each launch window is generally about 19 minutes to 44 minutes, and the data clearly show that Delta II set two fixed azimuth corresponding to two liftoff moments, respectively, composing a launch window to ensure the liftoff (and possibly using a similar launch mode of Atlas V when there is sufficient launch capability), and injects the probe directly to the hyperbolic transfer orbit toward Mars.

    c) Atlas V, as a flagship of U.S. launch vehicle now, for the sufficient launch capability and guidance control ability, 20 days launch period is generally selected, and daily launch window varies from 30 minutes to 2 hours Maven even reached a daily window of 2 hours. According to the orbital inclination required by Mars exploration, Atlas V can accommodate dynamical change of the liftoff time within the launch window via a fixed azimuth, adjustment of dog-leg maneuver and coasting time. The data show that 94o or 104o azimuth are used for different missions with corresponding orbital inclination of 28.7o or 35.5o, respectively (the latter includes a primary yaw program). Atlas V sent the probe directly into the hyperbolic transfer orbit toward Mars.

    d) Due to PSLV’s limited launch capability, the probe is sent into a sub-GTO orbit firstly, then carries out multi-perigee burns lasting for several days, finally entering Mars exploration hyperbolic escape trajectory. This can be catalogued into a modified version of the Russian launch mode.

    e) Six days during the 20 days launch period are selected by ESA, and a deep-space escape hyperbolic trajectory with a launch window of 20~30 minutes per day is designed. It is likely that the method of ‘fixed azimuth + dog-leg maneuver + coast time adjustment’ is adopted to ensure the launch window.

    Overall, except for different relative positions between the Earth and Mars, typical mission profiles and launch window extensions of Earth-Mars launch trajectories depend closely on each rocket’s launch capability and guidance ability of navigation systems onboard.

    [1] Liu L, Hou X Y. Orbit mechanics of deep space explorations[ M]. Beijing: PHEI Press, 2015.

    [2] Izzo D. Revisiting lambert’s problem[J]. Celestial Mechanics and Dynamical Astronomy, 2015, 121: 1-15.

    [3] Battin R H. An introduction to the mathematics and methods of astrodynamics[M]. New York: AIAA, 1999.

    [4] Hu W D. Fundamental spacecraft dynamics and control[M]. Wiley: Singapore, 2015.

    [5] Graham W. Proton-M successfully launches first ExoMars spacecraf[EB/OL]. (2017-07-03 )[2017-07-03].https://www.nasaspaceflight. com/2016/03/proton-m-first-exomars-spacecraft/.

    [6] Spaceflight101. MAVEN - Atlas V 401 - Ascent Info[EB/OL]. (2017-07-03) [2017-07-03].http://www.spaceflight101.net/atlas-v-401-maven-ascent-info.html.

    [7] Delta Launch Vehicle Programs.PHOENIX.[EB/OL]. (2017-08-11)< [2017-08-11].www.ulalaunch.com.

    [8] Nampoothiri M V. PSLV-C25: the vehicle that launched the Indian Mars Orbiter[J]. Current Science, 2015, 109.

    [9] Arianespace. Ariane flight 158[EB/OL]. (2017-07-02) [2017-07-02]. http://www.arianespace.com/wp-content/uploads/2017/06/04_feb_26-en.pdf.

    世界典型火星探測發(fā)射軌道及窗口拓展分析

    耿光有1,2,王 玨1,3,宋 強(qiáng)2,張志國2,王建明2

    (1. 北京航空航天大學(xué)宇航學(xué)院,北京,100191;2. 北京宇航系統(tǒng)工程研究所,北京,100076;3.中國運(yùn)載火箭技術(shù)研究院,北京,100076)

    由于運(yùn)載火箭自身特點(diǎn)及飛行航、落區(qū)約束的差異,致使世界典型運(yùn)載火箭直接火星探測任務(wù)的發(fā)射軌道存在明顯不同。從直接火星轉(zhuǎn)移的飛行力學(xué)原理的角度進(jìn)行了分析,重點(diǎn)關(guān)注了典型火星探測任務(wù)的軌道剖面設(shè)計(jì)及窗口拓展方法,完成了對當(dāng)前火星探測任務(wù)典型出發(fā)軌道特點(diǎn)的分析總結(jié)。

    地球-火星發(fā)射軌道;轉(zhuǎn)移軌道;發(fā)射窗口;雙曲逃逸軌道

    V41

    A

    2018-01-03

    10.7654/j.issn.1004-7182.20180104

    1004-7182(2018)01-0018-07

    耿光有(1972-),男,博士,研究員,主要研究方向?yàn)檫\(yùn)載火箭飛行力學(xué)及總體設(shè)計(jì)研究

    Geng Guangyou (1972-), male, Ph.D, research professor, mainly focus on flight mechanics and systems design of launch vehicle

    猜你喜歡
    北京航空航天大學(xué)宇航力學(xué)
    力學(xué)
    《北京航空航天大學(xué)學(xué)報(bào)》征稿簡則
    《北京航空航天大學(xué)學(xué)報(bào)》征稿簡則
    《北京航空航天大學(xué)學(xué)報(bào)》征稿簡則
    《北京航空航天大學(xué)學(xué)報(bào)》征稿簡則
    弟子規(guī)·余力學(xué)文(十)
    快樂語文(2021年11期)2021-07-20 07:41:32
    弟子規(guī)·余力學(xué)文(四)
    快樂語文(2020年30期)2021-01-14 01:05:28
    力學(xué) 等
    我的宇航夢
    我的宇航夢
    国产精品久久久久久亚洲av鲁大| 成人永久免费在线观看视频| 久久精品国产亚洲av天美| 三级毛片av免费| 久99久视频精品免费| 男插女下体视频免费在线播放| 国产亚洲av嫩草精品影院| 免费搜索国产男女视频| 国产黄色小视频在线观看| 国产精品精品国产色婷婷| 欧美性猛交╳xxx乱大交人| 亚洲最大成人av| 噜噜噜噜噜久久久久久91| 精品久久国产蜜桃| 男人舔女人下体高潮全视频| 97碰自拍视频| 日韩人妻高清精品专区| 亚洲国产高清在线一区二区三| .国产精品久久| 九九在线视频观看精品| 欧美+日韩+精品| 日韩在线高清观看一区二区三区| 99久久精品一区二区三区| 可以在线观看毛片的网站| 日本一本二区三区精品| 久久久a久久爽久久v久久| 一本精品99久久精品77| 如何舔出高潮| 高清毛片免费看| 夜夜夜夜夜久久久久| 99久国产av精品国产电影| 一级毛片我不卡| 最后的刺客免费高清国语| 久久天躁狠狠躁夜夜2o2o| 日韩av在线大香蕉| 一区二区三区四区激情视频 | 精品久久国产蜜桃| 激情 狠狠 欧美| 日本a在线网址| 亚洲aⅴ乱码一区二区在线播放| 国产精品嫩草影院av在线观看| 日本免费一区二区三区高清不卡| 欧美区成人在线视频| 国产精品亚洲一级av第二区| 1000部很黄的大片| 永久网站在线| 精品熟女少妇av免费看| 欧美三级亚洲精品| 国产极品精品免费视频能看的| 久久热精品热| 最好的美女福利视频网| eeuss影院久久| 最新中文字幕久久久久| 久久久成人免费电影| 久久久成人免费电影| 黑人高潮一二区| 嫩草影院精品99| 亚洲av熟女| 久久久久九九精品影院| 国产真实乱freesex| 国产亚洲欧美98| 欧美极品一区二区三区四区| 亚洲国产精品久久男人天堂| 日本一本二区三区精品| 中文资源天堂在线| 白带黄色成豆腐渣| 两个人的视频大全免费| 亚洲精品日韩在线中文字幕 | 天天躁夜夜躁狠狠久久av| 国产精品无大码| 成人国产麻豆网| 日本精品一区二区三区蜜桃| 国产精品三级大全| 校园人妻丝袜中文字幕| 午夜福利在线观看吧| 亚洲国产精品成人久久小说 | 在线播放无遮挡| 色综合色国产| 日本-黄色视频高清免费观看| 自拍偷自拍亚洲精品老妇| 精品人妻视频免费看| 五月伊人婷婷丁香| av免费在线看不卡| 国产一区二区在线av高清观看| 亚洲av一区综合| 美女内射精品一级片tv| 给我免费播放毛片高清在线观看| 1024手机看黄色片| 国产精品伦人一区二区| 精品久久久久久久久久久久久| 91在线观看av| 国产高清不卡午夜福利| 老熟妇仑乱视频hdxx| 99久久精品热视频| 俄罗斯特黄特色一大片| 精品久久久久久久末码| 美女 人体艺术 gogo| 国内久久婷婷六月综合欲色啪| 人妻丰满熟妇av一区二区三区| 在线播放无遮挡| 熟女电影av网| 国产高清视频在线播放一区| 欧美又色又爽又黄视频| 久久国产乱子免费精品| 一级毛片电影观看 | 啦啦啦观看免费观看视频高清| 国产精品亚洲美女久久久| 小说图片视频综合网站| 日日啪夜夜撸| 日韩制服骚丝袜av| 欧美性感艳星| 国产精品一区www在线观看| 精品99又大又爽又粗少妇毛片| 97超碰精品成人国产| 精品福利观看| 国产av在哪里看| 日韩制服骚丝袜av| 嫩草影院精品99| 国产精品久久视频播放| 精品一区二区三区视频在线| 女的被弄到高潮叫床怎么办| 成人欧美大片| 看免费成人av毛片| av在线蜜桃| 少妇丰满av| 亚洲无线在线观看| 亚洲精品影视一区二区三区av| 国产片特级美女逼逼视频| 欧美中文日本在线观看视频| 国产单亲对白刺激| 精品久久久噜噜| 不卡视频在线观看欧美| 我的老师免费观看完整版| 亚洲成人精品中文字幕电影| 国产高清视频在线观看网站| 两个人视频免费观看高清| 能在线免费观看的黄片| 国产精品永久免费网站| 亚洲无线观看免费| 精品不卡国产一区二区三区| 搡老岳熟女国产| 长腿黑丝高跟| 一夜夜www| 亚洲精品国产av成人精品 | 国产高清视频在线播放一区| 欧美日韩精品成人综合77777| 99久久中文字幕三级久久日本| 欧美性感艳星| 好男人在线观看高清免费视频| 成人二区视频| 99热精品在线国产| 欧美最新免费一区二区三区| 波多野结衣高清无吗| 久久中文看片网| 俄罗斯特黄特色一大片| 国内精品宾馆在线| 中出人妻视频一区二区| 淫妇啪啪啪对白视频| av福利片在线观看| av天堂在线播放| 99久久中文字幕三级久久日本| 搡老妇女老女人老熟妇| 国产真实乱freesex| 久久精品国产鲁丝片午夜精品| 色5月婷婷丁香| 三级国产精品欧美在线观看| 久久这里只有精品中国| 我要看日韩黄色一级片| 插阴视频在线观看视频| 久久精品国产99精品国产亚洲性色| 国产av麻豆久久久久久久| 在线a可以看的网站| 精品人妻视频免费看| 十八禁网站免费在线| 日本a在线网址| 好男人在线观看高清免费视频| 精品午夜福利视频在线观看一区| 国产在线精品亚洲第一网站| 日韩国内少妇激情av| 波多野结衣高清无吗| 久久久精品94久久精品| 国产片特级美女逼逼视频| 在线播放国产精品三级| 亚洲欧美日韩东京热| 精品人妻视频免费看| 精品久久国产蜜桃| 性色avwww在线观看| 午夜精品国产一区二区电影 | 亚洲av免费高清在线观看| 亚洲av电影不卡..在线观看| 免费av观看视频| 91久久精品电影网| 国产成人91sexporn| 精品久久久久久久末码| 成年av动漫网址| 欧美日本视频| 天天躁夜夜躁狠狠久久av| 国产精品国产三级国产av玫瑰| 国产精品女同一区二区软件| 禁无遮挡网站| 在线播放国产精品三级| 天天一区二区日本电影三级| 大型黄色视频在线免费观看| 久久精品国产亚洲av香蕉五月| 日韩欧美免费精品| 美女 人体艺术 gogo| 久久久久精品国产欧美久久久| 在线观看免费视频日本深夜| 在线免费观看的www视频| 亚洲精品影视一区二区三区av| 99热这里只有精品一区| 人人妻人人看人人澡| 神马国产精品三级电影在线观看| 亚洲自偷自拍三级| 免费人成视频x8x8入口观看| 亚洲av电影不卡..在线观看| 亚洲色图av天堂| 国产三级在线视频| 丰满乱子伦码专区| 国产片特级美女逼逼视频| 亚洲图色成人| 免费av不卡在线播放| 久久婷婷人人爽人人干人人爱| 一夜夜www| 男人狂女人下面高潮的视频| 久久久久久久久中文| 身体一侧抽搐| 国产av不卡久久| 在线免费观看不下载黄p国产| 中文字幕av在线有码专区| 成人av一区二区三区在线看| 有码 亚洲区| 男人和女人高潮做爰伦理| 老师上课跳d突然被开到最大视频| 免费av观看视频| 日日啪夜夜撸| or卡值多少钱| 成年女人毛片免费观看观看9| 欧美日韩国产亚洲二区| 国产一区二区在线观看日韩| 亚洲七黄色美女视频| 在线播放国产精品三级| 国产精品1区2区在线观看.| 少妇丰满av| 男插女下体视频免费在线播放| 啦啦啦观看免费观看视频高清| 精品人妻偷拍中文字幕| 亚洲综合色惰| 中国美女看黄片| 精品久久久久久久久久免费视频| 身体一侧抽搐| 成人特级av手机在线观看| 亚洲精品456在线播放app| 欧美激情国产日韩精品一区| 欧美bdsm另类| videossex国产| 久久久久精品国产欧美久久久| 亚洲人成网站在线播放欧美日韩| 床上黄色一级片| 成人鲁丝片一二三区免费| 免费av不卡在线播放| 深夜精品福利| 俺也久久电影网| 亚洲三级黄色毛片| 22中文网久久字幕| 狂野欧美激情性xxxx在线观看| 99热这里只有是精品在线观看| 精品少妇黑人巨大在线播放 | 久久精品国产99精品国产亚洲性色| 国产精品一区二区三区四区免费观看 | 成人毛片a级毛片在线播放| 欧美精品国产亚洲| 国内精品宾馆在线| 亚洲欧美日韩高清在线视频| 我的老师免费观看完整版| 日本精品一区二区三区蜜桃| 亚洲国产精品成人久久小说 | 啦啦啦韩国在线观看视频| 久久精品人妻少妇| 国产精品爽爽va在线观看网站| 日韩 亚洲 欧美在线| 欧美+日韩+精品| 村上凉子中文字幕在线| 免费大片18禁| 亚洲久久久久久中文字幕| 精品不卡国产一区二区三区| 香蕉av资源在线| 久久国内精品自在自线图片| 淫妇啪啪啪对白视频| 免费在线观看影片大全网站| 免费av毛片视频| 亚洲va在线va天堂va国产| 国产高清三级在线| 悠悠久久av| 亚洲无线在线观看| 亚洲欧美精品综合久久99| 欧美激情久久久久久爽电影| 国产精品免费一区二区三区在线| 插阴视频在线观看视频| 久久久久精品国产欧美久久久| 99视频精品全部免费 在线| 国产欧美日韩精品亚洲av| 99热这里只有是精品在线观看| 国产毛片a区久久久久| 搡女人真爽免费视频火全软件 | 免费观看的影片在线观看| 久久久久久久亚洲中文字幕| 免费观看在线日韩| 亚洲av美国av| 日韩 亚洲 欧美在线| 亚洲,欧美,日韩| 久久亚洲精品不卡| 成人无遮挡网站| 午夜精品一区二区三区免费看| 国产亚洲91精品色在线| 久久久久九九精品影院| 十八禁国产超污无遮挡网站| 欧美一区二区精品小视频在线| 色综合站精品国产| 久久99热6这里只有精品| 精品少妇黑人巨大在线播放 | 99久国产av精品| 国产老妇女一区| 长腿黑丝高跟| 在线免费观看不下载黄p国产| 不卡视频在线观看欧美| 高清日韩中文字幕在线| 免费观看的影片在线观看| 看片在线看免费视频| 国产一区二区三区av在线 | 日韩高清综合在线| 国产探花在线观看一区二区| 亚洲精品一区av在线观看| 嫩草影院入口| 国产欧美日韩精品一区二区| 村上凉子中文字幕在线| 99riav亚洲国产免费| 国产精品不卡视频一区二区| 久久亚洲精品不卡| 亚洲最大成人手机在线| 波多野结衣巨乳人妻| 成人三级黄色视频| 插阴视频在线观看视频| 日日摸夜夜添夜夜添小说| 亚洲欧美日韩无卡精品| 色播亚洲综合网| 天堂网av新在线| 精品人妻视频免费看| 毛片女人毛片| 免费一级毛片在线播放高清视频| 黄片wwwwww| 国产欧美日韩一区二区精品| 欧美+日韩+精品| www.色视频.com| 桃色一区二区三区在线观看| 伊人久久精品亚洲午夜| 久久久久久伊人网av| 人妻少妇偷人精品九色| 黄色配什么色好看| 久久人人爽人人爽人人片va| a级毛色黄片| 全区人妻精品视频| 我要看日韩黄色一级片| 直男gayav资源| 不卡一级毛片| 综合色av麻豆| 亚洲美女搞黄在线观看 | 91久久精品国产一区二区成人| 人妻丰满熟妇av一区二区三区| 亚洲精品日韩av片在线观看| 亚洲精品影视一区二区三区av| 在线播放无遮挡| 国产精品无大码| 国产亚洲精品久久久com| 国产精品99久久久久久久久| 高清日韩中文字幕在线| 亚洲成a人片在线一区二区| 国产成人精品久久久久久| 日本a在线网址| 99久久精品一区二区三区| 97人妻精品一区二区三区麻豆| 97热精品久久久久久| 中文字幕久久专区| 日产精品乱码卡一卡2卡三| 国产成人a∨麻豆精品| 搡老熟女国产l中国老女人| 色综合亚洲欧美另类图片| a级毛色黄片| 亚洲人成网站在线播放欧美日韩| 九九热线精品视视频播放| 欧美3d第一页| 免费观看在线日韩| 午夜亚洲福利在线播放| 日韩大尺度精品在线看网址| www日本黄色视频网| 在线免费十八禁| 国产熟女欧美一区二区| 嫩草影院精品99| 国产 一区 欧美 日韩| 欧美日韩精品成人综合77777| 中文字幕精品亚洲无线码一区| av.在线天堂| 亚洲国产精品成人综合色| a级一级毛片免费在线观看| 男女做爰动态图高潮gif福利片| 国产中年淑女户外野战色| 亚洲专区国产一区二区| 九色成人免费人妻av| 国产av一区在线观看免费| 国产真实伦视频高清在线观看| 欧美色视频一区免费| 黄色配什么色好看| 久久久久久国产a免费观看| 人妻久久中文字幕网| 91久久精品电影网| 国产精品久久视频播放| 午夜久久久久精精品| 久久久久国产精品人妻aⅴ院| 欧美日韩在线观看h| 看非洲黑人一级黄片| 美女黄网站色视频| 少妇的逼好多水| 观看美女的网站| 给我免费播放毛片高清在线观看| 老司机影院成人| 最新在线观看一区二区三区| 老熟妇乱子伦视频在线观看| 国产亚洲91精品色在线| 成年免费大片在线观看| 高清午夜精品一区二区三区 | 亚洲熟妇中文字幕五十中出| 久久草成人影院| av天堂中文字幕网| 舔av片在线| 精品人妻一区二区三区麻豆 | 久久久精品欧美日韩精品| 久久精品国产亚洲av天美| 亚洲美女黄片视频| 亚洲国产色片| 亚洲av二区三区四区| 国产av不卡久久| 亚洲av不卡在线观看| 亚洲激情五月婷婷啪啪| 日本一二三区视频观看| 亚洲国产精品合色在线| 欧美一区二区亚洲| 啦啦啦韩国在线观看视频| 久久久成人免费电影| 久久久国产成人精品二区| 国产三级中文精品| 日韩欧美精品免费久久| 欧美激情国产日韩精品一区| 久久精品影院6| 99热这里只有精品一区| 99国产极品粉嫩在线观看| 久久久国产成人免费| 精品不卡国产一区二区三区| 91在线精品国自产拍蜜月| 亚洲欧美成人精品一区二区| 国产高清激情床上av| 香蕉av资源在线| 国产久久久一区二区三区| 国产一区二区三区av在线 | 在线免费十八禁| .国产精品久久| 综合色丁香网| 国产综合懂色| 在线观看一区二区三区| 欧美色欧美亚洲另类二区| 麻豆久久精品国产亚洲av| 高清毛片免费观看视频网站| 少妇熟女aⅴ在线视频| 欧美性感艳星| 国产在线男女| 哪里可以看免费的av片| 精品久久久久久成人av| 亚洲成人中文字幕在线播放| 能在线免费观看的黄片| 简卡轻食公司| 欧美激情国产日韩精品一区| 久久久精品94久久精品| 啦啦啦观看免费观看视频高清| 99久久精品热视频| 欧美极品一区二区三区四区| 91av网一区二区| aaaaa片日本免费| 日韩人妻高清精品专区| 男人和女人高潮做爰伦理| 国产精品亚洲一级av第二区| 欧美zozozo另类| 91狼人影院| 中国国产av一级| 国产精品野战在线观看| 亚洲18禁久久av| 国产大屁股一区二区在线视频| 亚洲欧美成人综合另类久久久 | 九九久久精品国产亚洲av麻豆| 亚洲美女搞黄在线观看 | 欧美xxxx性猛交bbbb| 色哟哟哟哟哟哟| 久久午夜亚洲精品久久| 你懂的网址亚洲精品在线观看 | 国产免费一级a男人的天堂| 午夜老司机福利剧场| 伦理电影大哥的女人| 欧美成人a在线观看| 国产精品1区2区在线观看.| 中国国产av一级| 国产高清三级在线| 日日撸夜夜添| 天堂影院成人在线观看| 亚洲美女视频黄频| 午夜福利成人在线免费观看| 搡老岳熟女国产| 免费看av在线观看网站| 欧美日韩在线观看h| 一级av片app| 无遮挡黄片免费观看| 久久精品国产99精品国产亚洲性色| 深夜a级毛片| 啦啦啦啦在线视频资源| 午夜影院日韩av| 成人鲁丝片一二三区免费| 久久人人精品亚洲av| 免费看av在线观看网站| www日本黄色视频网| 一级av片app| 亚洲成人中文字幕在线播放| 国产单亲对白刺激| 午夜影院日韩av| 久久精品国产99精品国产亚洲性色| 午夜影院日韩av| 成人二区视频| 亚洲精品日韩av片在线观看| 久久精品国产99精品国产亚洲性色| 深夜a级毛片| 蜜桃久久精品国产亚洲av| 男女下面进入的视频免费午夜| 内地一区二区视频在线| 日韩精品有码人妻一区| 中文字幕免费在线视频6| 久久精品国产自在天天线| 国产91av在线免费观看| 一级a爱片免费观看的视频| 极品教师在线视频| 美女大奶头视频| 亚洲久久久久久中文字幕| 久久亚洲精品不卡| 又爽又黄a免费视频| 免费高清视频大片| 久久久久久久久中文| 欧美一区二区亚洲| 欧美+日韩+精品| 欧美成人a在线观看| 成人特级黄色片久久久久久久| 久久久久国内视频| 久久精品91蜜桃| 午夜福利视频1000在线观看| 我的老师免费观看完整版| 午夜激情欧美在线| 国产av在哪里看| 老司机午夜福利在线观看视频| 亚洲成人精品中文字幕电影| 啦啦啦啦在线视频资源| 一级a爱片免费观看的视频| 国产91av在线免费观看| 国产精品av视频在线免费观看| 中文在线观看免费www的网站| 国产私拍福利视频在线观看| 国产精品一区二区免费欧美| 亚洲国产欧洲综合997久久,| 国产探花在线观看一区二区| 少妇的逼好多水| 波多野结衣巨乳人妻| av在线天堂中文字幕| 国产免费男女视频| 日日摸夜夜添夜夜爱| 亚洲成人av在线免费| 给我免费播放毛片高清在线观看| 人妻夜夜爽99麻豆av| 久99久视频精品免费| 亚洲精品在线观看二区| 在线看三级毛片| 亚洲成人精品中文字幕电影| 美女高潮的动态| 成人国产麻豆网| 最近最新中文字幕大全电影3| 国产一区亚洲一区在线观看| 99精品在免费线老司机午夜| 国产一区二区在线av高清观看| 看黄色毛片网站| 一夜夜www| 婷婷精品国产亚洲av| 最近在线观看免费完整版| 中文字幕av成人在线电影| 成人二区视频| 热99在线观看视频| 精品人妻一区二区三区麻豆 | 欧美成人a在线观看| 国产亚洲精品久久久久久毛片| 一本一本综合久久| 亚洲精品亚洲一区二区| 欧美丝袜亚洲另类| 亚洲精品在线观看二区| 色av中文字幕| 色吧在线观看| 成人精品一区二区免费| 欧美一区二区亚洲| 黄色一级大片看看| 性色avwww在线观看| 日韩欧美 国产精品| 日本撒尿小便嘘嘘汇集6| 亚洲成人中文字幕在线播放| 直男gayav资源| 精品久久久久久久久av|