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

    Research on active arc-ignition technology as a possible residual-energy-release strategy in electromagnetic rail launch

    2021-08-05 08:29:58XiyuanCAI蔡喜元SaiTAN譚賽JunyongLU魯軍勇XiaoZHANG張曉andYunGUO郭赟
    Plasma Science and Technology 2021年8期
    關(guān)鍵詞:張曉

    Xiyuan CAI (蔡喜元), Sai TAN (譚賽), Junyong LU (魯軍勇),Xiao ZHANG (張曉) and Yun GUO (郭赟)

    National Key Laboratory of Science and Technology on Vessel Integrated Power System,Naval University of Engineering, Wuhan 430033, People’s Republic of China

    Abstract In order to solve the problem of the original residual energy release strategy being unsuitable for high-energy and fast-firing electromagnetic rail launch, this work has explored the applicability of active arc-ignition technology(AAT).The results obtained from the comparison of launching experiments show that AAT has no influence on the acceleration of the armature and is capable of quickly releasing the residual energy.Based on the theory of magnetohydrodynamics, this work has also made numerical simulation of the muzzle arc, analyzed the influence of AAT on the muzzle arc flow field, electromagnetic (EM) field and temperature field, and evaluated the performance of AAT according to the projectile initial disturbance, the EM impact on guidance devices and the rail ablation.The results show that AAT is now one of the most practicable strategies for residual energy release.

    Keywords: muzzle arc, initial disturbance, guidance device, arc ablation, backflow

    1.Introduction

    During the electromagnetic rail launch (EMRL), the energy remaining in the launch system is more than a quarter or even a half of the exit kinetic energy after the armature leaves the bore[1], and it needs to be released.However, the energy level of EMRL is continuously increased to 32 MJ, even 64 MJ, and the speed of continuous firing is up to six shots per minute[2,3].How does one release the residual energy but not affect the launch safety? Firing accuracy and service life of the launcher is one of the key points in the design of high energy level and fast firing rate electromagnetic (EM) rail launchers.

    The main criteria for judging whether a residual energy release strategy is good or not are as follows: (1) the release speed of residual energy; (2) the impact of residual energy release on the acceleration of the projectile in the chamber;(3)the impact on the initial disturbance of the projectile; (4) the electromagnetic impact on the guidance devices; and(5) the ablation of the rails.

    At present, one of the strategies for releasing residual energy is the arc suppression technology [4–7].When the projectile is accelerated in bore, the arc suppression impedance has a shunt effect on the rail current [7], thus reducing the driving force of the projectile.Because of the high energy Joule heat, the arc suppression device must meet the high requirements for temperature resistance and recovery ability.Another strategy is the passive arc-ignition technology [8].This method does not affect the acceleration of the projectile due to the fact that the arc can quickly release a large amount of residual energy and the air medium has strong capacity of recovery, but it has a seriously ablative effect on the rails,thereby affecting the service life of the launcher.The above strategies are not suitable for high-energy and fast continuous-firing EM rail launchers.

    Figure 1.Schematic diagram of the launcher and the AID.

    Active arc-ignition technology (AAT) is a new one for residual energy release.The arc-ignition device (AID) connected in parallel with the rails at the muzzle end releases the residual energy in the form of arc.The metal vapor generated by friction and wear of the aluminum armature in the bore enters the arcing gap and reduces the breakdown strength.Shortly after the armature leaves the bore, the muzzle overvoltage breaks through the gap between the arcing contacts,so that the arc between the armature and the rails will rapidly transfer to the AID.AAT is not used for suppressing the arc—a medium for rapid energy release,but for introducing it into the low-cost and replaceable AID.

    The influence of AAT on the rate of residual energy release and the acceleration of the projectile remains to be examined.What is worth noting is that the influence of AAT on the initial disturbance of the projectile mainly depends on the distribution of the pressure and the velocity of the muzzle arc, the electromagnetic impact caused by AAT on guidance device is determined by the distribution of the EM field, and the ablation caused by AAT to the trajectory is related to the temperature and airflow velocity.Therefore,it is necessary to study the influence of AAT on the launching process as well as the flow field, EM field and temperature field of the muzzle arc.

    In EMRL experiment, the observation of muzzle arc tends to be affected by the radiation of thermal plasma,especially the characteristics of arc flow field are difficult to obtain.Weimeret al[9] used an EM rail launcher with an emission energy level of 12 MJ and a peak current of 1.5 MA as an object of experiment.It is found that a stream of highintensity light, smoke, flame front and luminescent particles flows from the muzzle back to the chamber,with the backflow temperature in the range of 2400–4000 K, and the alkali metals, AlO and CuO existing in the gas.

    The increase of the computing capability by computer leads to the continuous progress in the numerical simulation of arc plasma.Using the magnetohydrodynamics (MHD)model of arc, it is possible to give a good description of the characteristics of the distribution of arc temperature field,EM field and flow field [10].Gaoet al[11] conducted a simulation study on the muzzle arc, but the muzzle arc root in their model only exists between the rails, and the current is only about 200 A, which is much lower than the magnitude of current 100 kA of the EMRL when AAT is adopted.This case is quite different from the muzzle arc produced in the use of AAT.The muzzle arc behavior needs to be further studied.

    This paper has described the current characteristics of the EM rail launcher in experiment when AAT is adopted and not, investigated the influence of AAT on the release rate of residual energy and the acceleration of the projectile and simulated the muzzle arc.According to the characteristics of muzzle arc, an MHD model has been established, and the finite volume method is adopted for solution.The model can be used to analyze the flow field, EM field and temperature field of the muzzle arc, and to evaluate the influence of AAT on the initial disturbance of the projectile, the EM impact of the guidance device and the rail ablation.

    2.Experiment

    An EM rail launcher with the caliber of 30 mm×30 mm and an AID,as shown in figure 1,has been designed.The AID is mainly composed of two pairs of arcing contacts and two copper connection boards, and the arcing contacts are connected with the rails via the connection board.After the armature leaves the bore, the current will be transferred from the rails to the arcing contacts, and the arc between the armature and the rails will be switched to the arcing gap.Two groups of comparative launching experiments are carried out respectively with and without the AID.Rogowski coils are used to measure the currents of the breech and the arcing contacts.

    Figure 2 shows the currents of the breech, the AID and the armature during the launching.It can be seen that the breech current is the same when the AID is used or not.The discharge characteristics of the power supply are not affected.Before the armature leaves the bore,there is no current in the device.Attout, the total current in the device rapidly rises from zero to the magnitude of the breech current.Aftertout,both the currents are basically the same.This shows that(1)in the acceleration process of the armature,the AID which is in a dormant state, does not affect the driving force of the armature; (2) at the moment of the armature leaving the bore, the device is quickly activated and receives all the current in the rail to ensure the safe disconnection between the armature and the rails;and(3)after the armature comes out of the bore,the device continues to work to finish releasing the residual energy in a short time.

    Figure 2.Current curves in experiment.

    3.Calculation model

    3.1.Hypotheses

    As the physical and chemical processes in the muzzle arc of an EM rail launcher are too complex for simulation calculation, it is necessary to simplify the physical process by introducing some reasonable hypotheses.The specific hypotheses are as follows:

    ? The arc plasma is considered to be in the local thermodynamic equilibrium (LTE) state.

    ? Ignoring the interaction between the arc and the electrode as well as the electrode sheath,considering the air area at the muzzle as the calculation domain.

    3.2.Governing equations

    MHD basic equations include mass conservation,momentum conservation,energy conservation,electric field and magnetic field equations.All these equations are unified in the steady form of the Patankar equation:

    where Φ is the variable to be solved,ΓΦis the diffusion term,andSΦis the source term.The variables and parameters of the MHD equations are listed in table 1.Among them,P,Hand V are respectively the plasma pressure, specific enthalpy and velocity vector;u,vandware the velocity components inx,yandzdirections, respectively;Ax,AyandAzare the magnetic vector potential components;Jx,JyandJzare the current density components;Bx,By,Bzare the magnetic field strength components; ρ,Cp, λ, μ and σ are respectively the plasma density, specific heat at constant pressure, thermal conductivity and viscosity;kbis Boltzmann constant,eis the electron charge, μ0is the permeability;qradis the radiation source term, which is related to the radiation model.

    The net emission coefficient εNis used to describe the radiation heat loss of plasma per unit volume, which is regarded as the source term of energy conservation equation,namely

    The realizablek-ε model is used to describe the turbulent flow state of muzzle arc.The related parameters and expressions can be obtained from [12] and [13].

    The temperature and the pressure near the wall are the highest and the gradient is large.Therefore,the enhanced wall function method is used for treatment.

    3.3.Geometry and mesh

    The geometry and calculation domain is shown in figure 3.According to the symmetry, a quarter of the air area is taken for the calculation model, in which the in-bore length is 1000 mm,and the radius and length of the air area outside the bore are both 1000 mm.The bolts are neglected.The cuboid is the arcing contact and L-type geometry is half of the armature.According to the features of the geometrical model,the computational domain is divided by the multi-block structural mesh.In order to increase the convergence of the calculation model and the calculation accuracy, the core area is super-refined.According to the requirements of the selected turbulence model and the adjacent-wall method on they+value, repeated trial calculations as well as grid adjustment are conducted and finally they+ value is about 1.They+value is within five,near the insulator wall in the breech end.Considering that the Reynolds number is small and far away from the muzzle end in this place,they+value is considered to be acceptable.The grids with the quality of 0.85–0.90 and 0.90–0.95 account for 0.032% and 6.379% respectively, and the quality of the other grids is above 0.95.The orthogonality of the mesh is greater than 36 degrees,and finally the number of grids is up to 10 million.

    3.4.Plasma properties

    The aluminum vapor tends to affect the physical properties of plasma medium in the AID.In this paper, the proportion of aluminum and air in the muzzle arc medium are assumed to be 10% and 90%, respectively.Murphy and others [14–19]offered the data,such as thermodynamic properties,transport coefficients and net emission coefficients of Al-air plasma,in the temperature range of 300–50000 K and in the pressure range of 0.01–20 MPa.Plasma properties are related to temperature and pressure.The binary three-point interpolation is adopted to obtain the physical parameters of plasma under the condition of specific temperature and pressure.

    Table 1.MHD basic equations.

    Figure 3.Calculation domain.

    3.5.Boundary conditions

    Simonyan’s and Kats’studies[20]show that,under the action of plasma, the maximum temperature on the surface of the electrode is close to the boiling point of electrode material.Based on that, the surface temperature of the rail, armature and arcing contact is set as the corresponding material boiling point (2500 K).In experiment, the insulation material was found to be hardly ablated, so it could be treated as the adiabatic wall [21].The outlet pressure is 101.325 kPa and the temperature is 300 K.

    After the armature leaves the bore,the arcing contacts on the single side each share half of the rail current.The current is distributed on the interface between the arcing contact and the arcing gap.The current is distributed parabolically.The potential of the armature surface is set to be zero.The potential and magnetic vector components of the outlet are zero,and the magnetic field on the wall is calculated by Biot–Savart law.

    4.Results and analysis

    The equations mentioned above are solved by the adapted commercial code Fluent, which is based on the finite volume method.For convenience, the forward direction of the armature is defined as thexaxis, and the launcher end is the plane ofx=0.The horizontal direction is set asyaxis, and the symmetry of left and right is the plane ofy=0.The vertical direction is the forward direction of thezaxis.

    Figure 4.Airflow rate distribution of the muzzle area.

    4.1.Flow field distribution

    Figure 4 shows the flow rate distribution at the muzzle area with and without AAT.When AAT is not adopted, the airflow outside the chamber rushes to the end of the rails and dashes upon the wall,thereby forming a complex turbulence.A part of airflow is emitted in the forward direction of the armature, and another part flows back to the inner chamber.The high-velocity region of the muzzle arc plasma is near the end of the rails, and the velocity is up to 8400 m s?1.When AAT is adopted, the airflow above the arc contacts will flow to the arcing gap.It will emit after it is offset in the arcing gap and the airflow moving to the plane ofy=0 will flow out along the vertical center of the launcher.The high velocity region of the muzzle arc plasma is in the arcing gap, and its magnitude will decrease to 5700 m s?1.The velocity of the airflow near the armature and rail end will drop down to 300 m s?1.

    Figure 5 shows the distribution of pressure in the muzzle area with and without AAT.When AAT is not used, the maximum pressure is 282 kPa and the pressure difference is 214 kPa,which appear near the muzzle end of the rails.When AAT is adopted,the maximum pressure and pressure gradient in the muzzle area will be shown in the arcing gap.The maximum pressure is 245 kPa and the pressure difference is 161 kPa.Furthermore, the pressure at the armature and rail end is uniform, being very close to the atmosphere.

    Airflow is the main reason that the initial disturbance of the projectile is affected by the muzzle arc.The use of AAT can lead to a decrease in the airflow velocity, turbulence,pressure and pressure difference in the muzzle area,especially around the armature, which will be beneficial to the stability of the projectile and reduce its initial disturbance, thereby improving the firing accuracy.

    4.2.EM field distribution

    Figures 6 and 7 show the distributions of EM field with and without AAT.When AAT is not applied,the EM field of the muzzle arc is mainly concentrated near the armature and rail end,with the peak current density up to 300 MA m?2and the magnitude of the magnetic field up to 0.67 T.When AAT is used,the muzzle arc is transferred to the AID.The maximum current density of the arcing gap is 230 MA m?2and the peak magnetic field is 0.41 T.The current near the armature is greatly reduced to about 0.2 MA m?2, and the magnetic field drops down to about 0.03 T.

    Figure 5.Pressure distribution of the muzzle area.

    Figure 6.Current distribution of the muzzle area.

    Figure 7.Magnetic distribution of the muzzle area.

    The AID has the ability to shunt the current, transfer the muzzle arc.When AAT is not adopted, the current flows directly from one rail to the other, and the current is more concentrated.But if it is applied, each of the two pairs of arcing contacts will share only a half of the residual current.As a result, the arc current on the single side will be greatly reduced.In addition, the arcing contacts are distributed on both sides of the rail,so the muzzle arc is enlarged.The AID expands the muzzle arc from two aspects of current and space,and reduces the EM field in the muzzle area,especially near the armature.Therefore, the application of AAT can reduce the EM impact on the internal guidance device of the projectile to prevent it from losing its vitality.

    4.3.Temperature distribution

    Figure 8 shows the distribution of temperature in the muzzle area with and without AAT.When AAT is not used, the muzzle arc concentrates between the armature and the rails,with less expansion in the direction perpendicular to the armature movement.In the direction of the armature movement, the muzzle arc covers the armature when expanding forward and forms a backflow in the bore when expanding backwards, with its maximum temperature being 28000 K.When AAT is applied, the muzzle arc is concentrated in the arcing gap.Its maximum temperature is 25000 K, while the temperature of the backflow at the end of the rails and in the chamber is greatly reduced to 6000 K or so.

    The application of AAT can alleviate the ablation of the rails in the muzzle end.When subjected to the high-temperature arc plasma, the rail surface will undergo evaporation erosion and droplet splashing, which lead to a loss of materials.Because of the arc transfer of the AID,the arc column is shifted from between the rails to the arcing gap.Besides, the shunting effect of the AID can reduce the current and the temperature of single-sided arc, so the heat flux to the rail surface will be greatly reduced.According to the analysis in section 4.1, the adoption of AAT makes it possible to reduce the airflow velocity around the rails,and to weaken the impact of the thermal flow on the rail surface as well.Thus, the alleviation of the rail ablation will help prolong the service life of the device.

    5.Conclusions

    Because the original strategy of residual energy release is not suitable for high-energy and rapid EMRL, this paper proposed a new approach named AAT and investigated its performance.Through experiments and numerical calculations,AAT is proved to be applicable for dealing with the residual energy release, the acceleration of the armature, the initial disturbance of the projectile, the EM impact on the guidance devices, and the ablation of rail materials.

    Figure 8.Temperature distribution of the muzzle area.

    AAT has the following characteristics: when the armature accelerates in the bore, AID is in a dormant state; when the armature leaves the bore, the device quickly transfers the current, without affecting the power discharge and the armature acceleration.What is more, the residual energy can be released shortly after the armature leaves the bore.

    AAT has the ability to transfer the muzzle arc.It reduces the airflow velocity and pressure in the area around the rail end and armature, distributes the pressure more uniformly,weakens the EM field,and decreases the temperature greatly.As a result,it is possible to decrease the initial disturbance of the projectile,weaken the EM impact on the guidance device inside the projectile and alleviate the rail ablation, thereby improving the performance and firing accuracy of the launcher and prolonging the service life of the rail.In short,AAT is now one of the most practicable strategies of releasing the residual energy for high-energy, and rapid EMRL.

    Acknowledgments

    This work was supported in part by National Natural Science Foundation of China (Nos.51522706, 51877214 and 51607187) and in part by the National Basic Research Program of China (973 Program) (No.613262).

    猜你喜歡
    張曉
    Photothermal-chemical synthesis of P–S–H ternary hydride at high pressures
    Synthesis of ternary compound in H–S–Se system at high pressures?
    Investigation of electronic,elastic,and optical properties of topological electride Ca3Pb via first-principles calculations*
    Dielectric breakdown properties of Al-air mixtures
    “你可以逼我代孕,那為什么我不能借精生子?”
    Quaternary antiferromagnetic Ba2BiFeS5 with isolated FeS4 tetrahedra
    煮餃子
    春雨
    THE BEACH PARADOX
    考霸
    亚洲专区中文字幕在线| 中文字幕最新亚洲高清| 免费看十八禁软件| 久久亚洲真实| 大香蕉久久成人网| 国产精品美女特级片免费视频播放器 | 国产欧美亚洲国产| 一级片免费观看大全| 国产欧美日韩一区二区精品| 欧美大码av| 日日摸夜夜添夜夜添小说| 99精品在免费线老司机午夜| 丰满人妻熟妇乱又伦精品不卡| 国产激情久久老熟女| 首页视频小说图片口味搜索| 亚洲精品中文字幕一二三四区 | 一本久久精品| 黄色视频在线播放观看不卡| 激情视频va一区二区三区| 日韩视频在线欧美| 91麻豆精品激情在线观看国产 | 欧美日韩视频精品一区| 这个男人来自地球电影免费观看| 国产精品免费大片| 岛国在线观看网站| 久久久精品免费免费高清| 日韩欧美免费精品| 亚洲国产av影院在线观看| av国产精品久久久久影院| 国产亚洲欧美精品永久| 国产精品 欧美亚洲| 99国产精品一区二区蜜桃av | 亚洲一卡2卡3卡4卡5卡精品中文| 久久 成人 亚洲| 婷婷成人精品国产| 99精品欧美一区二区三区四区| 一级片免费观看大全| 欧美一级毛片孕妇| 亚洲精品一卡2卡三卡4卡5卡| 国产精品麻豆人妻色哟哟久久| 午夜福利免费观看在线| 久久久久久亚洲精品国产蜜桃av| 老司机亚洲免费影院| 国产精品亚洲一级av第二区| 亚洲午夜精品一区,二区,三区| 亚洲成人手机| 91大片在线观看| 国产高清videossex| 久久国产精品人妻蜜桃| 国产无遮挡羞羞视频在线观看| 91麻豆精品激情在线观看国产 | 欧美中文综合在线视频| 亚洲国产av影院在线观看| bbb黄色大片| 99香蕉大伊视频| 超色免费av| 大型av网站在线播放| 午夜福利在线免费观看网站| 午夜福利,免费看| 欧美在线黄色| 99国产极品粉嫩在线观看| 久久久久久免费高清国产稀缺| 在线观看66精品国产| 欧美激情久久久久久爽电影 | 国产片内射在线| 黄色a级毛片大全视频| 国产成人免费观看mmmm| 天天操日日干夜夜撸| 黄色视频不卡| 捣出白浆h1v1| 宅男免费午夜| 蜜桃在线观看..| 日韩欧美一区视频在线观看| a级毛片黄视频| 久久香蕉激情| 亚洲人成电影免费在线| 国产精品偷伦视频观看了| 免费人妻精品一区二区三区视频| 久久九九热精品免费| 丰满迷人的少妇在线观看| 黑人巨大精品欧美一区二区mp4| 午夜福利一区二区在线看| 丝袜人妻中文字幕| 91精品国产国语对白视频| 狠狠精品人妻久久久久久综合| 国产一区二区三区综合在线观看| 国产极品粉嫩免费观看在线| 一区二区三区国产精品乱码| 黑人巨大精品欧美一区二区mp4| 宅男免费午夜| 一区在线观看完整版| 黄色怎么调成土黄色| 老司机靠b影院| 成年人免费黄色播放视频| 91国产中文字幕| 免费久久久久久久精品成人欧美视频| 中文字幕av电影在线播放| 亚洲精品中文字幕在线视频| 日韩免费av在线播放| 国产免费现黄频在线看| 女人久久www免费人成看片| 色94色欧美一区二区| 90打野战视频偷拍视频| 国产福利在线免费观看视频| 如日韩欧美国产精品一区二区三区| 久久精品亚洲av国产电影网| 欧美日韩成人在线一区二区| av线在线观看网站| 涩涩av久久男人的天堂| 国产欧美日韩精品亚洲av| 51午夜福利影视在线观看| 国产av又大| 黑人猛操日本美女一级片| 亚洲中文日韩欧美视频| 成人特级黄色片久久久久久久 | 悠悠久久av| 精品人妻1区二区| 少妇被粗大的猛进出69影院| 亚洲精品国产精品久久久不卡| 国产精品av久久久久免费| 久久免费观看电影| 日本撒尿小便嘘嘘汇集6| 国产亚洲午夜精品一区二区久久| 看免费av毛片| 成年人午夜在线观看视频| 一区二区三区激情视频| 巨乳人妻的诱惑在线观看| 日韩熟女老妇一区二区性免费视频| 成人黄色视频免费在线看| 免费一级毛片在线播放高清视频 | 一本—道久久a久久精品蜜桃钙片| 久久精品亚洲精品国产色婷小说| 久久中文看片网| 午夜精品久久久久久毛片777| 少妇猛男粗大的猛烈进出视频| 精品亚洲乱码少妇综合久久| 久久99一区二区三区| 久久久久精品国产欧美久久久| 动漫黄色视频在线观看| 欧美日韩亚洲综合一区二区三区_| 亚洲综合色网址| av免费在线观看网站| 深夜精品福利| 久久久久国产一级毛片高清牌| 在线观看66精品国产| 国产在线视频一区二区| 大码成人一级视频| 久久婷婷成人综合色麻豆| 久久久国产精品麻豆| 国产单亲对白刺激| 美女视频免费永久观看网站| 最近最新中文字幕大全免费视频| 黑人欧美特级aaaaaa片| 女人高潮潮喷娇喘18禁视频| 一区二区三区精品91| 中文字幕av电影在线播放| tube8黄色片| 最近最新免费中文字幕在线| 色综合婷婷激情| 国产真人三级小视频在线观看| 桃红色精品国产亚洲av| 欧美 亚洲 国产 日韩一| 午夜福利一区二区在线看| 男人操女人黄网站| 岛国在线观看网站| 19禁男女啪啪无遮挡网站| 国产成人欧美在线观看 | avwww免费| 青草久久国产| 无限看片的www在线观看| 久久人人97超碰香蕉20202| 免费女性裸体啪啪无遮挡网站| 一级毛片电影观看| 午夜精品国产一区二区电影| av网站免费在线观看视频| 久久久久久久久免费视频了| 少妇精品久久久久久久| 亚洲av国产av综合av卡| 亚洲国产欧美网| 久久亚洲精品不卡| 亚洲人成伊人成综合网2020| 一级黄色大片毛片| 宅男免费午夜| 日韩一卡2卡3卡4卡2021年| 色老头精品视频在线观看| 亚洲久久久国产精品| 成人国产av品久久久| 久久人妻福利社区极品人妻图片| 黄色a级毛片大全视频| 男女午夜视频在线观看| 两性夫妻黄色片| 王馨瑶露胸无遮挡在线观看| 国产精品二区激情视频| 亚洲欧洲日产国产| 好男人电影高清在线观看| 黄色视频,在线免费观看| 在线观看免费日韩欧美大片| 99久久国产精品久久久| 国产亚洲午夜精品一区二区久久| 黑人操中国人逼视频| 亚洲全国av大片| 欧美激情 高清一区二区三区| 在线看a的网站| 久久精品亚洲av国产电影网| 在线播放国产精品三级| 老司机亚洲免费影院| 巨乳人妻的诱惑在线观看| 交换朋友夫妻互换小说| 色播在线永久视频| 99re6热这里在线精品视频| 色婷婷av一区二区三区视频| 不卡一级毛片| 一级片免费观看大全| 国产1区2区3区精品| 99re在线观看精品视频| 免费观看a级毛片全部| 免费日韩欧美在线观看| www.999成人在线观看| 老司机午夜福利在线观看视频 | 免费观看人在逋| 丰满人妻熟妇乱又伦精品不卡| 97人妻天天添夜夜摸| 老司机靠b影院| 国产人伦9x9x在线观看| 精品国产国语对白av| 丝瓜视频免费看黄片| 一级,二级,三级黄色视频| 色播在线永久视频| 亚洲熟妇熟女久久| 高清在线国产一区| 天天影视国产精品| 69精品国产乱码久久久| 亚洲成人免费电影在线观看| 亚洲一码二码三码区别大吗| 亚洲欧美一区二区三区久久| 亚洲成av片中文字幕在线观看| av又黄又爽大尺度在线免费看| 每晚都被弄得嗷嗷叫到高潮| 亚洲一区二区三区欧美精品| 免费日韩欧美在线观看| 99re6热这里在线精品视频| 欧美精品一区二区大全| 免费在线观看黄色视频的| 欧美另类亚洲清纯唯美| 啦啦啦中文免费视频观看日本| 亚洲欧美精品综合一区二区三区| 91成人精品电影| 美女高潮喷水抽搐中文字幕| 久久久久久久精品吃奶| 老司机午夜福利在线观看视频 | 在线观看66精品国产| 十八禁网站免费在线| 少妇的丰满在线观看| 亚洲精品在线美女| 久久免费观看电影| 岛国在线观看网站| 少妇裸体淫交视频免费看高清 | 最新的欧美精品一区二区| 91九色精品人成在线观看| 高清视频免费观看一区二区| av一本久久久久| 亚洲精品在线观看二区| 免费人妻精品一区二区三区视频| 亚洲精品一卡2卡三卡4卡5卡| 成在线人永久免费视频| 国产精品电影一区二区三区 | 1024视频免费在线观看| 欧美精品啪啪一区二区三区| 日韩视频一区二区在线观看| 久久 成人 亚洲| 在线永久观看黄色视频| 久久久久视频综合| 精品亚洲成国产av| 老司机午夜十八禁免费视频| 成人av一区二区三区在线看| 久久精品国产亚洲av香蕉五月 | 亚洲精品国产色婷婷电影| 老熟女久久久| 免费女性裸体啪啪无遮挡网站| 老鸭窝网址在线观看| 女性生殖器流出的白浆| 国产熟女午夜一区二区三区| 看免费av毛片| 少妇 在线观看| 亚洲国产看品久久| 免费女性裸体啪啪无遮挡网站| 久久香蕉激情| 免费av中文字幕在线| 国产亚洲精品一区二区www | 看免费av毛片| 国产精品1区2区在线观看. | 黄色怎么调成土黄色| 水蜜桃什么品种好| 丁香六月天网| 一夜夜www| 国产精品亚洲av一区麻豆| 国产成+人综合+亚洲专区| 乱人伦中国视频| 国产一区二区三区综合在线观看| 丝袜人妻中文字幕| 法律面前人人平等表现在哪些方面| 一本综合久久免费| 老熟妇仑乱视频hdxx| 久久久欧美国产精品| 最近最新中文字幕大全免费视频| e午夜精品久久久久久久| 99国产精品99久久久久| 999久久久精品免费观看国产| 成年人免费黄色播放视频| 十八禁人妻一区二区| 正在播放国产对白刺激| 久久久久久免费高清国产稀缺| 啦啦啦 在线观看视频| 无人区码免费观看不卡 | 精品亚洲乱码少妇综合久久| 亚洲成人国产一区在线观看| 亚洲av欧美aⅴ国产| 国产男女内射视频| 水蜜桃什么品种好| 午夜福利,免费看| 黄片大片在线免费观看| 天堂8中文在线网| 国产深夜福利视频在线观看| 国产欧美日韩一区二区精品| 国产成人欧美| 久久香蕉激情| 国产欧美日韩一区二区三区在线| 在线观看免费午夜福利视频| 欧美日韩亚洲高清精品| 波多野结衣一区麻豆| 午夜两性在线视频| 视频在线观看一区二区三区| 国产精品99久久99久久久不卡| 久久精品国产亚洲av高清一级| 丰满少妇做爰视频| 叶爱在线成人免费视频播放| 欧美黄色淫秽网站| 又黄又粗又硬又大视频| 国产在线观看jvid| 51午夜福利影视在线观看| 国产精品 欧美亚洲| 日韩一区二区三区影片| 亚洲精品久久成人aⅴ小说| 性少妇av在线| av天堂在线播放| 中文字幕人妻熟女乱码| 狠狠精品人妻久久久久久综合| 成人黄色视频免费在线看| 精品国产亚洲在线| 一区二区三区精品91| 欧美成人午夜精品| 国产精品成人在线| 亚洲国产欧美日韩在线播放| 成人黄色视频免费在线看| 精品国产超薄肉色丝袜足j| 在线看a的网站| 50天的宝宝边吃奶边哭怎么回事| 少妇 在线观看| 美女福利国产在线| 午夜精品国产一区二区电影| 蜜桃在线观看..| 日韩欧美一区视频在线观看| 日韩人妻精品一区2区三区| 制服人妻中文乱码| 三上悠亚av全集在线观看| 女人久久www免费人成看片| 国产精品久久电影中文字幕 | 老司机午夜十八禁免费视频| 欧美激情久久久久久爽电影 | 不卡av一区二区三区| 亚洲性夜色夜夜综合| 一区二区三区国产精品乱码| 国产欧美亚洲国产| 高清av免费在线| 老熟女久久久| 久久性视频一级片| 老司机午夜十八禁免费视频| 一本—道久久a久久精品蜜桃钙片| 国产亚洲av高清不卡| 亚洲av日韩精品久久久久久密| 中文字幕色久视频| 悠悠久久av| 美女视频免费永久观看网站| 大香蕉久久网| 久久久国产一区二区| 午夜视频精品福利| 天天操日日干夜夜撸| 久久99热这里只频精品6学生| 国产又色又爽无遮挡免费看| 国产97色在线日韩免费| 免费少妇av软件| 人妻久久中文字幕网| 日本欧美视频一区| 色尼玛亚洲综合影院| 日韩欧美一区视频在线观看| 最近最新中文字幕大全电影3 | 妹子高潮喷水视频| 18禁美女被吸乳视频| 亚洲色图综合在线观看| 大型黄色视频在线免费观看| 亚洲五月天丁香| 99久久99久久久精品蜜桃| 国产蜜桃级精品一区二区三区| 可以在线观看的亚洲视频| 亚洲精品一卡2卡三卡4卡5卡| 制服丝袜大香蕉在线| 丁香六月欧美| 精品国产美女av久久久久小说| 久久久国产成人免费| 黄色日韩在线| 国产精品自产拍在线观看55亚洲| 国产日本99.免费观看| 中文亚洲av片在线观看爽| 亚洲国产欧美一区二区综合| 又紧又爽又黄一区二区| 男人和女人高潮做爰伦理| 国产精品综合久久久久久久免费| aaaaa片日本免费| 女人高潮潮喷娇喘18禁视频| 亚洲欧美日韩高清专用| 国产视频内射| 国内精品久久久久精免费| 成人三级黄色视频| 香蕉国产在线看| 午夜日韩欧美国产| 在线观看一区二区三区| 给我免费播放毛片高清在线观看| 成人欧美大片| 久久久国产欧美日韩av| 国模一区二区三区四区视频 | 九九在线视频观看精品| 亚洲精品美女久久av网站| 18禁黄网站禁片午夜丰满| 国产免费男女视频| 色综合亚洲欧美另类图片| 别揉我奶头~嗯~啊~动态视频| 在线国产一区二区在线| 桃色一区二区三区在线观看| 国产精品一区二区免费欧美| 国产爱豆传媒在线观看| 免费观看人在逋| 天堂影院成人在线观看| 97碰自拍视频| 亚洲aⅴ乱码一区二区在线播放| 看免费av毛片| 91av网一区二区| 成年人黄色毛片网站| 国产伦在线观看视频一区| 色吧在线观看| 亚洲一区高清亚洲精品| 国产一区二区三区视频了| 欧美成人性av电影在线观看| 国产av一区在线观看免费| 午夜视频精品福利| 欧美激情在线99| 中文字幕高清在线视频| 91av网一区二区| 亚洲欧美激情综合另类| 夜夜看夜夜爽夜夜摸| 最近最新中文字幕大全免费视频| 日韩欧美在线乱码| 中文字幕人妻丝袜一区二区| 亚洲一区高清亚洲精品| 999久久久精品免费观看国产| 黄频高清免费视频| 超碰成人久久| 俺也久久电影网| 999久久久国产精品视频| xxxwww97欧美| 国产主播在线观看一区二区| 午夜精品一区二区三区免费看| 麻豆国产97在线/欧美| 亚洲一区二区三区不卡视频| 精品久久久久久久末码| 国产精品综合久久久久久久免费| 狂野欧美激情性xxxx| 色综合亚洲欧美另类图片| 99热精品在线国产| 男女之事视频高清在线观看| 中文亚洲av片在线观看爽| 国产成人aa在线观看| 国产精品久久久av美女十八| 欧美大码av| 欧美一区二区国产精品久久精品| 少妇熟女aⅴ在线视频| 国产成人影院久久av| 亚洲 国产 在线| 91麻豆av在线| 日本精品一区二区三区蜜桃| 国产亚洲精品av在线| 母亲3免费完整高清在线观看| 国产一区二区三区视频了| 99久久精品国产亚洲精品| 久久精品aⅴ一区二区三区四区| 精品一区二区三区av网在线观看| 俄罗斯特黄特色一大片| 精品不卡国产一区二区三区| 日韩有码中文字幕| 啦啦啦免费观看视频1| xxx96com| 国产精品av久久久久免费| av福利片在线观看| 亚洲av成人av| 成人性生交大片免费视频hd| 亚洲中文日韩欧美视频| 搞女人的毛片| 19禁男女啪啪无遮挡网站| 国产成人啪精品午夜网站| 日韩 欧美 亚洲 中文字幕| 久久久久性生活片| 午夜福利在线在线| 国产成人av教育| 精品电影一区二区在线| 欧美日韩福利视频一区二区| 一区二区三区高清视频在线| 桃红色精品国产亚洲av| 国产精品久久久久久亚洲av鲁大| 日本在线视频免费播放| 亚洲av日韩精品久久久久久密| 国产精品一区二区三区四区免费观看 | 成人三级黄色视频| 99热这里只有是精品50| 香蕉丝袜av| 夜夜看夜夜爽夜夜摸| xxxwww97欧美| 久久中文看片网| 在线观看舔阴道视频| 午夜福利在线观看吧| www.自偷自拍.com| 国内久久婷婷六月综合欲色啪| 亚洲精品粉嫩美女一区| 欧美不卡视频在线免费观看| 精品一区二区三区视频在线 | 日本a在线网址| 国产一级毛片七仙女欲春2| 成年免费大片在线观看| 欧美午夜高清在线| 中文资源天堂在线| 日本一二三区视频观看| 日韩成人在线观看一区二区三区| 91老司机精品| 一夜夜www| 身体一侧抽搐| 中出人妻视频一区二区| 国产成人精品久久二区二区免费| 亚洲aⅴ乱码一区二区在线播放| 12—13女人毛片做爰片一| 欧美成人免费av一区二区三区| АⅤ资源中文在线天堂| 欧美不卡视频在线免费观看| 成人性生交大片免费视频hd| av欧美777| 精品久久久久久久人妻蜜臀av| 69av精品久久久久久| 一区二区三区高清视频在线| 国产精品久久电影中文字幕| 久久久国产精品麻豆| 在线a可以看的网站| 午夜激情福利司机影院| 国产精品久久久人人做人人爽| 99久久99久久久精品蜜桃| a在线观看视频网站| 日本五十路高清| 九色国产91popny在线| 国产伦人伦偷精品视频| 老汉色∧v一级毛片| 成熟少妇高潮喷水视频| 亚洲18禁久久av| 757午夜福利合集在线观看| tocl精华| 美女午夜性视频免费| 精品国产乱子伦一区二区三区| 国产高清有码在线观看视频| 国内少妇人妻偷人精品xxx网站 | 十八禁网站免费在线| 久久久久国产精品人妻aⅴ院| 可以在线观看的亚洲视频| 99久久精品热视频| 好看av亚洲va欧美ⅴa在| 男人和女人高潮做爰伦理| 色噜噜av男人的天堂激情| 可以在线观看毛片的网站| 最近最新免费中文字幕在线| 欧美又色又爽又黄视频| 黑人欧美特级aaaaaa片| 亚洲九九香蕉| 午夜福利在线观看免费完整高清在 | 午夜a级毛片| 国产成人精品久久二区二区免费| 高清毛片免费观看视频网站| 丁香欧美五月| 看片在线看免费视频| 国产成人影院久久av| 亚洲成人中文字幕在线播放| 99久久精品一区二区三区| 久久久久国内视频| 丁香欧美五月| 少妇的逼水好多| 久久性视频一级片| 一a级毛片在线观看| 看黄色毛片网站| 青草久久国产| 精品国产三级普通话版| 欧美乱码精品一区二区三区| 男人和女人高潮做爰伦理| av欧美777| 欧美乱码精品一区二区三区| 亚洲欧美日韩高清专用| 九九久久精品国产亚洲av麻豆 | 女生性感内裤真人,穿戴方法视频| 变态另类成人亚洲欧美熟女| 最新在线观看一区二区三区| 床上黄色一级片| 九九在线视频观看精品| 成人三级做爰电影|