• <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
    考霸
    老熟妇仑乱视频hdxx| 亚洲无线在线观看| 免费av不卡在线播放| 亚洲av中文字字幕乱码综合| 久久天躁狠狠躁夜夜2o2o| 色在线成人网| 亚洲人成网站高清观看| 亚洲熟女毛片儿| 看免费av毛片| av黄色大香蕉| 亚洲国产精品久久男人天堂| 99久久国产精品久久久| 岛国在线免费视频观看| 偷拍熟女少妇极品色| 婷婷六月久久综合丁香| 大型黄色视频在线免费观看| 成人特级黄色片久久久久久久| 激情在线观看视频在线高清| 欧美精品啪啪一区二区三区| 亚洲av成人精品一区久久| 亚洲国产欧美一区二区综合| 久久中文看片网| 午夜免费成人在线视频| 国内少妇人妻偷人精品xxx网站 | 怎么达到女性高潮| 精品免费久久久久久久清纯| 热99在线观看视频| 18禁黄网站禁片免费观看直播| 国产精品一区二区免费欧美| 欧美日韩福利视频一区二区| 99久久无色码亚洲精品果冻| 亚洲狠狠婷婷综合久久图片| 久久久久亚洲av毛片大全| 两性午夜刺激爽爽歪歪视频在线观看| 欧美极品一区二区三区四区| 黄片大片在线免费观看| 国内毛片毛片毛片毛片毛片| 国产精品一及| 日本五十路高清| 老熟妇仑乱视频hdxx| 黄色 视频免费看| 免费一级毛片在线播放高清视频| 国产亚洲欧美98| 国模一区二区三区四区视频 | 男女那种视频在线观看| 欧美黄色淫秽网站| 久久久色成人| 给我免费播放毛片高清在线观看| 性色avwww在线观看| 这个男人来自地球电影免费观看| 最近视频中文字幕2019在线8| 又粗又爽又猛毛片免费看| netflix在线观看网站| 国产精品一区二区免费欧美| 精品久久久久久久久久免费视频| 亚洲精品乱码久久久v下载方式 | 色吧在线观看| 日本五十路高清| 手机成人av网站| 国产午夜精品论理片| 在线国产一区二区在线| 亚洲av成人一区二区三| 欧美高清成人免费视频www| 国产av麻豆久久久久久久| 精品久久久久久久久久久久久| 99热这里只有精品一区 | 久久国产精品人妻蜜桃| 婷婷丁香在线五月| 精品久久蜜臀av无| 国产视频内射| 三级毛片av免费| 天天躁日日操中文字幕| 久久精品夜夜夜夜夜久久蜜豆| 99国产极品粉嫩在线观看| www.自偷自拍.com| 久久精品人妻少妇| 757午夜福利合集在线观看| 精品一区二区三区av网在线观看| 中出人妻视频一区二区| 别揉我奶头~嗯~啊~动态视频| 亚洲黑人精品在线| 90打野战视频偷拍视频| 黄片小视频在线播放| 亚洲国产看品久久| 国产免费av片在线观看野外av| 婷婷精品国产亚洲av在线| 欧美日韩精品网址| 波多野结衣巨乳人妻| 欧美国产日韩亚洲一区| x7x7x7水蜜桃| 又爽又黄无遮挡网站| 黄色片一级片一级黄色片| 啪啪无遮挡十八禁网站| 国产高清视频在线观看网站| 日韩 欧美 亚洲 中文字幕| 久久久久久人人人人人| 黄色日韩在线| 香蕉久久夜色| 亚洲欧美日韩高清在线视频| 成熟少妇高潮喷水视频| e午夜精品久久久久久久| 美女cb高潮喷水在线观看 | 国产在线精品亚洲第一网站| 国产在线精品亚洲第一网站| 啦啦啦韩国在线观看视频| 国产久久久一区二区三区| 可以在线观看的亚洲视频| 免费观看的影片在线观看| 国产成人影院久久av| 精品欧美国产一区二区三| 成人av在线播放网站| 午夜福利视频1000在线观看| 国产亚洲av嫩草精品影院| 日本一二三区视频观看| 非洲黑人性xxxx精品又粗又长| 久久久国产成人精品二区| 成年女人看的毛片在线观看| 久久亚洲精品不卡| 国内精品一区二区在线观看| 午夜久久久久精精品| 欧美3d第一页| 男人和女人高潮做爰伦理| 精品久久久久久久毛片微露脸| 国产精品久久久人人做人人爽| 国产三级中文精品| 麻豆久久精品国产亚洲av| 欧美日韩瑟瑟在线播放| 亚洲午夜理论影院| www国产在线视频色| 毛片女人毛片| 一个人免费在线观看电影 | 18禁黄网站禁片午夜丰满| 91久久精品国产一区二区成人 | 亚洲自偷自拍图片 自拍| 国产综合懂色| 99热这里只有精品一区 | 99久久综合精品五月天人人| 97超级碰碰碰精品色视频在线观看| 国产亚洲av嫩草精品影院| 免费人成视频x8x8入口观看| 黄色视频,在线免费观看| 国产午夜精品久久久久久| 国产高清激情床上av| 国产男靠女视频免费网站| 国产精品99久久久久久久久| 亚洲中文日韩欧美视频| 欧美大码av| 国产单亲对白刺激| 熟女人妻精品中文字幕| 久久久久九九精品影院| 欧美在线一区亚洲| 日韩欧美免费精品| 亚洲av成人av| 久久中文看片网| 国产精品av视频在线免费观看| 99re在线观看精品视频| 美女大奶头视频| 国产精品影院久久| 1000部很黄的大片| 在线a可以看的网站| 搞女人的毛片| 少妇的逼水好多| 91老司机精品| 伦理电影免费视频| 蜜桃久久精品国产亚洲av| 狂野欧美激情性xxxx| 美女高潮的动态| 美女高潮的动态| 国内精品久久久久久久电影| 一个人观看的视频www高清免费观看 | 久久这里只有精品中国| 动漫黄色视频在线观看| 国产精品一区二区三区四区免费观看 | 精品国产乱子伦一区二区三区| 亚洲精品乱码久久久v下载方式 | 国产精品久久视频播放| a级毛片在线看网站| or卡值多少钱| 亚洲精品美女久久久久99蜜臀| 日韩欧美一区二区三区在线观看| 一级a爱片免费观看的视频| 国产av在哪里看| 波多野结衣高清无吗| 国产三级中文精品| www.自偷自拍.com| 日韩欧美国产一区二区入口| 色哟哟哟哟哟哟| 国产视频一区二区在线看| 亚洲国产精品合色在线| 欧美高清成人免费视频www| 亚洲av成人精品一区久久| 国产淫片久久久久久久久 | 白带黄色成豆腐渣| 在线十欧美十亚洲十日本专区| 国产精品乱码一区二三区的特点| 亚洲va日本ⅴa欧美va伊人久久| 超碰成人久久| 老熟妇乱子伦视频在线观看| h日本视频在线播放| 日韩精品中文字幕看吧| 亚洲va日本ⅴa欧美va伊人久久| 亚洲av片天天在线观看| 少妇的逼水好多| 日韩欧美在线乱码| 亚洲无线观看免费| 亚洲欧美日韩东京热| 欧美精品啪啪一区二区三区| 国产精品亚洲美女久久久| 国产精品野战在线观看| 国产激情偷乱视频一区二区| 日日摸夜夜添夜夜添小说| 可以在线观看毛片的网站| 亚洲av片天天在线观看| 欧美乱色亚洲激情| 色老头精品视频在线观看| 国产视频一区二区在线看| 夜夜看夜夜爽夜夜摸| 一区二区三区激情视频| 大型黄色视频在线免费观看| 亚洲av成人不卡在线观看播放网| 伊人久久大香线蕉亚洲五| 国产乱人视频| 国产精品久久久久久亚洲av鲁大| 免费看a级黄色片| 人人妻人人澡欧美一区二区| www日本黄色视频网| 久久天堂一区二区三区四区| 天堂动漫精品| 亚洲精品粉嫩美女一区| 国产一区二区激情短视频| 午夜福利免费观看在线| 欧美日韩乱码在线| 岛国在线免费视频观看| 中亚洲国语对白在线视频| 美女cb高潮喷水在线观看 | 香蕉久久夜色| 国产乱人伦免费视频| 又黄又粗又硬又大视频| 国产熟女xx| 99国产精品一区二区三区| 国产精品一区二区精品视频观看| 亚洲精品在线美女| 婷婷精品国产亚洲av在线| 国产精品香港三级国产av潘金莲| 国产精品av视频在线免费观看| 男人舔女人下体高潮全视频| 国产伦精品一区二区三区视频9 | 成人永久免费在线观看视频| 韩国av一区二区三区四区| 天天一区二区日本电影三级| 女同久久另类99精品国产91| 国内少妇人妻偷人精品xxx网站 | 久久久久九九精品影院| 久99久视频精品免费| 成年女人毛片免费观看观看9| 久久久久精品国产欧美久久久| 一区二区三区国产精品乱码| 免费在线观看视频国产中文字幕亚洲| 两性午夜刺激爽爽歪歪视频在线观看| 国产乱人伦免费视频| 精品国产亚洲在线| 久久久久亚洲av毛片大全| 在线观看一区二区三区| 天堂√8在线中文| 亚洲狠狠婷婷综合久久图片| 岛国在线免费视频观看| 中国美女看黄片| 麻豆国产av国片精品| 两人在一起打扑克的视频| 国产高清videossex| 国产1区2区3区精品| 国产精品国产高清国产av| 亚洲国产高清在线一区二区三| 国产极品精品免费视频能看的| 嫩草影视91久久| 亚洲国产看品久久| 亚洲精品美女久久av网站| 偷拍熟女少妇极品色| www.www免费av| 热99re8久久精品国产| 精品久久蜜臀av无| 中文在线观看免费www的网站| 亚洲在线观看片| 超碰成人久久| 老汉色∧v一级毛片| 欧美av亚洲av综合av国产av| 久久中文看片网| 久久香蕉国产精品| 久久久久久大精品| 老鸭窝网址在线观看| 久久中文字幕一级| 中国美女看黄片| 久久午夜综合久久蜜桃| 亚洲 欧美一区二区三区| 欧美激情久久久久久爽电影| 全区人妻精品视频| 亚洲真实伦在线观看| 久久99热这里只有精品18| 欧美一区二区国产精品久久精品| 色播亚洲综合网| 午夜激情欧美在线| 一本精品99久久精品77| 亚洲最大成人中文| 亚洲av成人不卡在线观看播放网| 久久精品国产亚洲av香蕉五月| 久久久久久久午夜电影| 亚洲av第一区精品v没综合| 国产精华一区二区三区| 国产精品亚洲av一区麻豆| 亚洲欧美日韩高清在线视频| 欧美不卡视频在线免费观看| 97碰自拍视频| 国产精品,欧美在线| 叶爱在线成人免费视频播放| 18禁裸乳无遮挡免费网站照片| 桃色一区二区三区在线观看| 精品久久久久久久久久免费视频| 国产精品一区二区三区四区免费观看 | 好男人电影高清在线观看| 亚洲片人在线观看| 精品国产乱码久久久久久男人| 香蕉国产在线看| 熟女少妇亚洲综合色aaa.| 啦啦啦韩国在线观看视频| 国语自产精品视频在线第100页| 欧美日本亚洲视频在线播放| 免费高清视频大片| 免费人成视频x8x8入口观看| 夜夜夜夜夜久久久久| 亚洲欧美日韩高清在线视频| 老司机深夜福利视频在线观看| 一区二区三区高清视频在线| 丰满的人妻完整版| 91av网站免费观看| 中文字幕精品亚洲无线码一区| 啦啦啦韩国在线观看视频| 国产精品亚洲一级av第二区| 欧美av亚洲av综合av国产av| 国产精品免费一区二区三区在线| 91在线精品国自产拍蜜月 | 十八禁人妻一区二区| 成人亚洲精品av一区二区| 看黄色毛片网站| 一区二区三区激情视频| 巨乳人妻的诱惑在线观看| 噜噜噜噜噜久久久久久91| 一本一本综合久久| 伦理电影免费视频| 国产精品98久久久久久宅男小说| 亚洲国产精品sss在线观看| 嫩草影院精品99| 久久九九热精品免费| 亚洲avbb在线观看| 国产成+人综合+亚洲专区| 国产精品综合久久久久久久免费| 18禁国产床啪视频网站| 99久久无色码亚洲精品果冻| 免费看日本二区| 国产亚洲精品综合一区在线观看| 精品国产超薄肉色丝袜足j| 久久精品aⅴ一区二区三区四区| 制服丝袜大香蕉在线| e午夜精品久久久久久久| 免费一级毛片在线播放高清视频| 欧美三级亚洲精品| 国产蜜桃级精品一区二区三区| 男女之事视频高清在线观看| xxxwww97欧美| 亚洲中文日韩欧美视频| 国产极品精品免费视频能看的| 国产一区二区在线观看日韩 | 黑人巨大精品欧美一区二区mp4| 久久精品夜夜夜夜夜久久蜜豆| 欧美成人性av电影在线观看| 中出人妻视频一区二区| 亚洲国产日韩欧美精品在线观看 | 可以在线观看毛片的网站| 日韩中文字幕欧美一区二区| 中文字幕人成人乱码亚洲影| 精品国产乱码久久久久久男人| 国产日本99.免费观看| 亚洲国产欧美网| 日本五十路高清| 宅男免费午夜| 国模一区二区三区四区视频 | 女人高潮潮喷娇喘18禁视频| 亚洲精品色激情综合| 男人的好看免费观看在线视频| 日韩欧美国产在线观看| 国产亚洲av高清不卡| 免费av毛片视频| 99riav亚洲国产免费| 高清在线国产一区| 三级男女做爰猛烈吃奶摸视频| 91麻豆精品激情在线观看国产| 岛国在线免费视频观看| 99久久综合精品五月天人人| 国产亚洲精品一区二区www| 亚洲国产精品sss在线观看| 国产成人精品无人区| 亚洲av电影不卡..在线观看| 岛国在线观看网站| 天天躁狠狠躁夜夜躁狠狠躁| 90打野战视频偷拍视频| 日本 欧美在线| 亚洲国产精品久久男人天堂| 国产av麻豆久久久久久久| 又紧又爽又黄一区二区| 黄片小视频在线播放| 欧美乱妇无乱码| 国产伦精品一区二区三区四那| 亚洲 国产 在线| 综合色av麻豆| e午夜精品久久久久久久| 老司机深夜福利视频在线观看| 欧美在线黄色| 99国产精品一区二区蜜桃av| 国产久久久一区二区三区| 久久香蕉精品热| 老汉色av国产亚洲站长工具| 亚洲五月婷婷丁香| 日日夜夜操网爽| 免费看光身美女| 99re在线观看精品视频| 嫩草影院入口| xxx96com| 国产精品综合久久久久久久免费| 欧美性猛交黑人性爽| 色视频www国产| 一个人看的www免费观看视频| 香蕉av资源在线| 久9热在线精品视频| 国产亚洲精品综合一区在线观看| 亚洲国产精品sss在线观看| 夜夜爽天天搞| 九色国产91popny在线| 久久久久久九九精品二区国产| 极品教师在线免费播放| 一卡2卡三卡四卡精品乱码亚洲| 麻豆成人av在线观看| 97超视频在线观看视频| 18禁黄网站禁片午夜丰满| 成人永久免费在线观看视频| 亚洲美女视频黄频| 精品欧美国产一区二区三| 国产亚洲精品av在线| 五月伊人婷婷丁香| 一边摸一边抽搐一进一小说| 成人一区二区视频在线观看| 久久久国产欧美日韩av| 欧美zozozo另类| 最近最新中文字幕大全电影3| 精品电影一区二区在线| 亚洲av第一区精品v没综合| 中文字幕最新亚洲高清| 99久久久亚洲精品蜜臀av| 久久中文字幕一级| 国产高清videossex| 男人舔女人的私密视频| 亚洲av第一区精品v没综合| 在线播放国产精品三级| 999久久久国产精品视频| 亚洲精品中文字幕一二三四区| av黄色大香蕉| 午夜福利高清视频| 亚洲欧美一区二区三区黑人| 国产精品一区二区免费欧美| 婷婷精品国产亚洲av在线| 91久久精品国产一区二区成人 | 久久中文字幕一级| www.熟女人妻精品国产| 亚洲专区中文字幕在线| 99国产精品一区二区蜜桃av| 欧美一级a爱片免费观看看| 国产精品久久久久久亚洲av鲁大| 精品久久久久久久末码| 一级毛片精品| 一级a爱片免费观看的视频| tocl精华| 国内精品一区二区在线观看| 中文字幕人妻丝袜一区二区| av片东京热男人的天堂| av福利片在线观看| 久久久久久大精品| 99久久成人亚洲精品观看| 好男人电影高清在线观看| 亚洲欧美日韩卡通动漫| 国产探花在线观看一区二区| 搡老岳熟女国产| 后天国语完整版免费观看| 又黄又粗又硬又大视频| 亚洲18禁久久av| 欧美3d第一页| 亚洲精品在线美女| 一个人免费在线观看电影 | 美女免费视频网站| 午夜精品在线福利| 成人高潮视频无遮挡免费网站| 亚洲精品美女久久av网站| 中出人妻视频一区二区| 精品一区二区三区视频在线观看免费| 欧美日韩乱码在线| 国产探花在线观看一区二区| 亚洲av免费在线观看| 怎么达到女性高潮| 哪里可以看免费的av片| 窝窝影院91人妻| 国产精品一区二区三区四区久久| 成熟少妇高潮喷水视频| 18禁观看日本| 在线观看午夜福利视频| 久久精品亚洲精品国产色婷小说| 97超视频在线观看视频| 婷婷丁香在线五月| www.自偷自拍.com| 免费观看的影片在线观看| 亚洲最大成人中文| 欧美日韩中文字幕国产精品一区二区三区| 午夜成年电影在线免费观看| 99国产极品粉嫩在线观看| 久久性视频一级片| 曰老女人黄片| 亚洲性夜色夜夜综合| 欧美一级a爱片免费观看看| 在线观看一区二区三区| 免费在线观看成人毛片| 国产99白浆流出| 老鸭窝网址在线观看| 久久精品亚洲精品国产色婷小说| 久久久水蜜桃国产精品网| 精品久久久久久久久久久久久| 午夜影院日韩av| 熟女人妻精品中文字幕| 在线观看免费视频日本深夜| а√天堂www在线а√下载| 中文字幕高清在线视频| 男女床上黄色一级片免费看| 蜜桃久久精品国产亚洲av| 久久久国产欧美日韩av| 国产美女午夜福利| 国产午夜福利久久久久久| 老鸭窝网址在线观看| 久久精品亚洲精品国产色婷小说| 19禁男女啪啪无遮挡网站| 国产综合懂色| 中文字幕人成人乱码亚洲影| aaaaa片日本免费| 婷婷精品国产亚洲av| 欧美黑人欧美精品刺激| 99热这里只有精品一区 | 小蜜桃在线观看免费完整版高清| 黄片小视频在线播放| 国产精品 欧美亚洲| 99re在线观看精品视频| 成人鲁丝片一二三区免费| 伊人久久大香线蕉亚洲五| 麻豆国产av国片精品| 精品国产超薄肉色丝袜足j| 色精品久久人妻99蜜桃| 岛国在线观看网站| 制服丝袜大香蕉在线| 久久精品国产综合久久久| 欧美一区二区国产精品久久精品| 天堂网av新在线| 亚洲,欧美精品.| 级片在线观看| 国产高清激情床上av| 午夜福利成人在线免费观看| 久久国产乱子伦精品免费另类| 成人亚洲精品av一区二区| 男人舔女人的私密视频| 中出人妻视频一区二区| 亚洲激情在线av| 可以在线观看毛片的网站| 日韩欧美免费精品| 亚洲专区国产一区二区| 欧美成人性av电影在线观看| 久久国产精品影院| 好男人在线观看高清免费视频| 大型黄色视频在线免费观看| 在线观看日韩欧美| 伊人久久大香线蕉亚洲五| 午夜免费激情av| 悠悠久久av| 国产成人一区二区三区免费视频网站| 神马国产精品三级电影在线观看| 亚洲欧美日韩高清专用| 三级国产精品欧美在线观看 | 亚洲午夜精品一区,二区,三区| 亚洲av第一区精品v没综合| 国产精品亚洲美女久久久| 嫩草影院精品99| 亚洲无线观看免费| 噜噜噜噜噜久久久久久91| netflix在线观看网站| 90打野战视频偷拍视频| 中文字幕高清在线视频| 国产一区二区三区在线臀色熟女| 国产欧美日韩精品一区二区| 视频区欧美日本亚洲| 亚洲av片天天在线观看| 国产精品精品国产色婷婷| 欧美绝顶高潮抽搐喷水| 一个人免费在线观看电影 | 亚洲av日韩精品久久久久久密| 久久久精品大字幕| 天堂√8在线中文| 日韩大尺度精品在线看网址| 国产黄a三级三级三级人| 91老司机精品| 日本精品一区二区三区蜜桃| 国产av麻豆久久久久久久| 亚洲av成人av| 看片在线看免费视频| or卡值多少钱|