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

    Hypersonic impact flash characteristics of a long-rod projectile collision with a thin plate target

    2021-03-23 13:58:12YijiangXueQingmingZhangDanyangLiuRenrongLongYangyuLuTianfeiRenLiangfeiGong
    Defence Technology 2021年2期

    Yi-jiang Xue,Qing-ming Zhang,Dan-yang Liu,Ren-rong Long,Yang-yu Lu,Tian-fei Ren,Liang-fei Gong

    State Key Laboratory of Explosion Science and Technology,Beijing Institute of Technology,Beijing 100081,China

    Keywords:Hypersonic impact flash Radiation intensity Impact flash mechanism Environmental luminescence

    ABSTRACT Impact flash occurs when objects collide at supersonic speeds and can be used for real-time damage assessment when weapons rely on kinetic energy to destroy targets.However,the mechanism of impact flash remains unclear.A series of impact flash experiments of flat-head long-rod projectiles impacting thin target plates were performed with a two-stage light gas gun.The impact flash spectra for 6061 aluminum at 1.3-3.2 km/s collision speeds were recorded with a high-speed camera,a photoelectric sensor,and a time-resolved spectrometer.The intensity of the impact flash exhibited a pulse characteristic with time.The intensity(I)increased with impact velocity(V0)according to I∝Vn0,where n=4.41 for V0>2 km/s.However,for V0<2 km/s,n=2.21,and the intense flash duration is an order of magnitude less than that of higher V0.When V0>2 km/s,a continuous spectrum(thermal radiation background)was observed and increased in intensity with V0.However,for V0<2 km/s,only atomic line spectra were detected.There was no aluminum spectral lines for V0<2 km/s,which indicated that it had not been vaporized.The initial intense flash was emission from excited and ionized ambient gases near the impact surface,and had little relationship with shock temperature rise,indicating a new mechanism of impact flash.?2020 The Authors.Production and hosting by Elsevier B.V.on behalf of China Ordnance Society.This is an open access article under the CC BY-NC-ND license(http://creativecommons.org/licenses/by-nc-nd/4.0/).

    1.Introduction

    The study of impact flashes originated from the Apollo program[1,2],and is still an important topic in astrophysics[3-10].The characteristics and spectrum of the impact flash can be used to determine the material composition and the surface structure of celestial bodies[11-13].It also has important applications in realtime damage assessment for kinetic energy weapons,such as armor-piercing projectiles[14]and missile defense systems[15],where velocities exceed 1 km/s.Impact flash is the most effective means of performing real-time damage assessment of the target and its authenticity,which can then provide a basis for decision making concerning second strikes.

    Many hypervelocity impact flash experiments using different materials have been performed[16-26].The impact flash radiation intensity has a pulse characteristic with time,and appears almost simultaneous with the projectile/target contact.Several or tens of microseconds later(determined by the impact velocity and target material),the radiation intensity reaches a peak value;then it decays to 20%(or less)of the peak intensity,after which it decays for milliseconds at a low intensity.The impact flash radiation intensity I and the impact velocity V0exhibit a power function relationship given by I∝Vn0.For V0in the range 3-7 km/s,n=8 for a copper ball and copper target collision[17],n=4.0-4.5 for an aluminum ball and aluminum target collision[20],and n=3.8-5 for iron particle collisions with a gold target[18].These results are very important references for the damage assessment of kinetic-energy-based weapons.However,the speed of the weapons is usually lower(V0<3 km/s)than that of celestial bodies.There are no systematic studies of impact flashes in this speed range,and almost no experimental data is available for reference.

    Plasma is the main source of the impact flash when V0>10 km/s[27,28].For V0=5-7 km/s,atomic aluminum lines at 396.15 nm and 394.40 nm were observed[17,26],and well as the copper atomic line at 510.55 nm[23],indicating gaseous aluminum or copper.Because the impact velocities of these tests were very high,it was inferred that the shock temperature vaporizes the material.However,the mechanism of the impact flash is unknown when the collision speed is reduced and the shock temperature rise is insufficient to vaporize the material.

    Here,a two-stage light gas gun was used to accelerate projectiles into targets.A high-speed camera,a photoelectric sensor,and a time-resolved spectrometer equipped with an intensified charge-coupled device(ICCD)were used to record impact flashes and their spectra for V0in the range 1.3-3.2 km/s.Time and frequency domain diagnostic equipment was used to multidimensionally analyze the characteristics of the hypersonic impact flash.The spectral analysis was used to identify the emitting material and its state,which could be used to infer the mechanism of the impact flash.

    2.In situ experiments

    2.1.Projectile and target

    Previous studies generally used spherical projectiles and thick targets[26].Here,flat-head,long-rod projectiles and thin targets were examined.According to Ref.[29],the original design of the long-rod projectile and the sabot structure are depicted in Fig.1(a),using sub-caliber launch technology.After multiple launches,the projectile cannot meet the experimental requirements,as shown in Fig.1(b).If the end of the long-rod is the surface that is mainly subjected to the force,it is easy to amplify the error caused by machining during the driving process.Thus,the flight posture is difficult to control after ejection of the long-rod projectile,and a very significant attack angle occurs.A long-rod projectile and sabot structure used after many improvements are shown in Fig.1(c).The force-bearing surface of the long-rod projectile and the sabot is parallel to the launch tube,and the shape centroid of the sabot is also coincident with the mass centroid of the long-rod projectile.The optimized experimental result is shown in Fig.1(d),where the long-rod projectile maintains the flight posture after leaving the tube muzzle.

    The long-rod projectile used here is shown in Fig.1(e).It has a diameter of 6 mm and a 50-mm length.The sabot has a diameter of 14.5 mm and a 20-mm length.The thickness of the target plate is 4 mm.Both the projectile and target are 6061 aluminum.The target plate is perpendicular to the projectile trajectory and is placed in a vacuum chamber where the ambient pressure is 100 Pa.The range of V0is 1.3-3.2 km/s.The long-rod projectile recovered by steelplate retardation is shown in Fig.1(f),the right image is the collision surface and the left is the tail of the projectile.

    Sabot detachment could easily affect the flight posture of the long-rod projectile,which in turn affects the repeatability of the impact flash experiments.The intense flash at the moment of collision is the main concern here,and there is a delay time for the sabot/target collision.Therefore,we did not take off the sabot,and successfully closed the electronic shutter of the ICCD spectrometer before the sabot hit the target plate,which avoided its affection on the spectra.The specific experimental results in Section 3.1 proved the feasibility of this scheme.

    Fig.1.(a)Original design of a long-rod projectile and the sabot structure.(b)Flight trajectory of the original design photographed with a high-speed camera.(c)Optimized design of the long-rod projectile and sabot.(d)Flight trajectory of the optimized design photographed with a high-speed camera.(e)6061 aluminum long-rod projectile and sabot used here.(f)The collision surface is at right,and the tail of the projectile is at left.

    2.2.Experimental system

    As shown in Fig.2,the experimental arrangement consists of three subsystems:launching,diagnostics,and triggering.The launching system consists of four parts:a two-stage light gas gun,a flying chamber,a target chamber,and a three-channel electromagnetic velocity measurement device(EVMD)at the exit of the launching tube to measure V0.The EVMD generates current when the projectile passes through.V0can be calculated from the time intervals of the currents in the three channels recorded by a digital storage oscilloscope,and the 50-mm distance between each channel.Because the ambient pressure of the flying chamber and the target chamber is only 100 Pa,speed attenuation of the projectile can be ignored and the EVMD data is directly used as the velocity V0at which the projectile collides with the plate target.

    The diagnostic system includes a high-speed camera(Photron Corp.,Japan),a silicon-based photoelectric sensor(Thorlabs Corp.,USA),and an ICCD spectrometer(Andor Corp.,UK).The photoelectric sensor has a 350-1050 nm range and a response time less than 10 ns.The impact flash irradiance is converted by the sensor into a voltage signal that is recorded over time with a digital storage oscilloscope.The ICCD spectrometer has a 200-1100 nm range,with a resolution less than 0.02 nm.Wavelength and intensity calibrations were performed separately before each experiment.The fiber optic probe of the ICCD spectrometer and the photoelectric sensor are symmetrically arranged on both sides of the impact point,and respectively record the spectral signal of impact flash and its intensity over time.The high-speed camera is outside the target chamber,and the impact flash is recorded through a window of the chamber for a side/front view of the collision,as depicted in Fig.2.The diagnostic equipment is synchronized by the triggering device,and the high-speed camera,photoelectric sensor,and ICCD spectrometer can be simultaneously trigged to record signals when the projectile is hitting the target.

    2.3.Triggering system

    A 10-μm-thick insulated copper foil is attached to the surface of the target plate.The foil and the metal target plate are connected to a resistance-capacitance(RC)pulse circuit that turns on instantly when the projectile hits the plate.A current pulse is thus formed and triggers the digital delay generator(DDG)to output multiple TTL-level signals,which simultaneously start the high-speed camera,the photoelectric sensor,and the ICCD spectrometer.The response time of the whole process is less than 100 ns,and the delay time between the DDG and the ICCD spectrometer is zero during experiments.For V0=3.0 km/s,signals with and without the copper foil output by the photoelectric sensor are plotted in Fig.3.The copper foil has no significant effect on the impact flash.

    Fig.3.Effect of copper foil on the impact flash.

    3.Results

    3.1.Impact flash process

    The impact and associated flash were recorded with the highspeed camera and the photoelectric sensor.Photographs of the impact flashes for V0=1.5 km/s and V0=2.5 km/s are shown in Fig.4.As shown in Fig.2,the camera is has a side/front view of the plate target.There is no background light source;therefore,the camera exposure relies entirely on the flash intensity.At time t=0,the projectile is starting to impact the target.For V0=1.5 km/s,a circular emission area appeared at the moment of collision(Fig.4a t=0).Because the interval between each camera frame was too long,changes in the emitting region were not observed.In the second photograph at 11.5μs,the impact flash has decayed,while in the third photograph at 23μs,material is ejected from the collision point,with almost no light emission.

    For V0=2.5 km/s(Fig.4(b)),the circular illuminating area was about twice the size of that observed for V0=1.5 km/s.The entire field of view in the second photograph at 16μs was overexposed because of the intense flash.The“overexposure”lasted less than 32μs,because the third photograph started to reveal some detail.In the fourth photograph at 96μs,the flash intensity was much weaker,and“dust plume”-like weak luminescence appeared.More high-speed camera results at different impact velocities are listed in Table 1.The results indicated that higher V0produced longer“overexposures”and flash durations.

    Impact flash signals vs.time recorded from the photoelectric sensor for various V0are shown in Fig.5.The sensor output voltage increases with the intensity of the impact flash.Because the photoelectric sensor was not calibrated,only relative values were obtained.The signal changes with time for different V0,and all have pulsed characteristics.The entire impact flash process can be separated into a pulse stage and an attenuation stage,as shown in Fig.5(b).More sensor results for different V0are shown in Table 1.They indicate that as V0increases,the peak radiation intensity increases along with the duration of the intense pulse flash.However,the time at which the peak intensity occurs(peak time)lags.When compared with the spherical aluminum projectile colliding with the plate target[1],the long-rod aluminum projectile colliding with the thin plate has a slightly longer pulse width,but still has the same order of magnitude.

    To analyze the time when the sabot hits the plate target at different V0,a model using the experimental dimensions was established for smoothed particle hydrodynamics(SPH)numerical simulations of the impact process.The calculated results for different V0are shown in Fig.6.The material model of 6061 aluminum was reported previously[30,31].During the collision,the plate target is greatly deformed by the projectile and many plate fragments are ejected.Some of the fragments move in a direction opposite to the impact velocity and eventually collide with the sabot.The time when the sabot and the debris began to interact was 7.4μs,4.8μs,and 3.8μs after impact at V0=1.5 km/s,2.4 km/s,and 3.0 km/s,respectively.The results in Fig.5 indicate that the flash signals were hardly affected by the sabot during the first half of the peak time.Therefore,in order to avoid or minimize sabot effects on the spectra,the electronic shutter width of the ICCD spectrometer was set to half of the peak time of the sensor signal for the corresponding impact velocity,as shown in Fig.5(c).

    3.2.Spectral characteristics

    Based on the relative intensity calibration,the spectrum of the V0=1.43 km/s impact flash captured by the ICCD spectrometer is plotted in Fig.7.The entire spectrum over the range 300-800 nm consisted of atomic lines;no bands or continuum spectra were observed.In addition to the two strong spectral lines,the average amplitude of the other spectral lines was 1600(a.u.).

    Fig.4.High-speed camera photographs of impact flashes taken from a side/front view of the plate target.The exposure time of each photograph was 1.5μs(a)Impact flash for V0=1.5 km/s at 100 Pa ambient pressure.(b)Impact flash for V0=2.4 km/s at 100 Pa ambient pressure.

    Table 1High-speed camera and photoelectric sensor results.

    The spectrum for V0=2.07 km/s is shown in Fig.8.It still had a large number of spectral lines,but there was also a continuous thermal radiation background over the range 400-800 nm.The average amplitude of the spectral lines was around 2000(a.u.),and that of thermal radiation was about 5%of that.For V0=3.12 km/s,the spectrum shown in Fig.9 had an average amplitude for the spectral lines of about 6000(au),and a thermal background amplitude that was 30%of that.

    The average irradiance measured by the optical fiber probe of the ICCD spectrometer was obtained by integrating the spectrum over 300-800 nm and dividing the result by the exposure time.Because the distance L0between the probe and the collision point was fixed and the light entering the area S of the optical fiber was much smaller than L02,the solid angle of each experiment was approximately the same.Thus,the irradiance can be considered equivalent to the relative intensity of the impact flash,and was plotted in Fig.10 for various V0.The black dots are for V0<2 km/s,and the red are for V0>2 km/s.In both regions,V0and the intensity I of the impact flash can be fit with a power function,but the relationship is not continuous.For V0<2 km/s,the fit is I=AV2.210,whereas,for V0>2 km/s,the fit is I=BV4.410,where A and B are constants.At 2 km/s,there is a transition region(blue dot)which cannot be fit by the two relations.

    In the fits,the larger exponent indicated a steeper increase in flash intensity with V0.The fit for V0>2 km/s was very close to the results for an aluminum ball projectile collision with an aluminum plate for V0<3-7 km/s[20].Even though projectiles had different shapes,the relationship between I and V0had the same fit with a similar exponent.Hence,the exponent likely characterizes aluminum projectiles,independent of shape.This finding is helpful for the study of inverse problems.

    3.3.Spectral analysis

    When the Solis spectrometer software(Andor Corp)was used to analyze the spectra for V0=1.43 km/s,none of the spectral lines could be attributed to aluminum,but rather to environmental gases.

    Fig.5.Photoelectric sensor signals vs.time for various impact velocities.

    Fig.6.SPH simulation results.(a)At t=7.4μs for V0=1.5 km/s.(b)At t=4.8μs for V0=2.4 km/s.(c)At t=3.8μs for V0=3.0 km/s.The simulation platform is ANSYS 14.0.

    Fig.11 shows the analysis results for nitrogen and oxygen.It can be seen that lots of spectral lines can be identified by neutral nitrogen or oxygen atoms and their ion states.It indicates the excitation and ionization of the ambient gas may be the main source of the impact flash.The other lines of the spectrum plot have the possibility to be identified by some of the inert gas elements.However,since the amount of inert gas is very little,the decision on the excitation of inert gas element should be cautious until further results are obtained.

    Because 6061 aluminum is an alloy with possible impurities and other elements,impact flash tests of pure aluminum(>99.9999%)were performed under the same conditions.The spectral analysis is shown in Fig.12,where again no aluminum lines were identified.Also tested was 45 steel at V0=1.30 km/s;the spectral analysis is shown in Fig.13 which still consists of line emissions by gaseous elements in the environment.These results indicate that ambient gas emission is likely to be common in hypersonic impact flashes of metallic materials.

    Fig.7.Spectrum for impact flash at V0=1.43 km/s.The ICCD gate width was 2μs.

    Fig.8.Spectrum for impact flash at V0=2.07 km/s.The ICCD gate width was 2μs.

    Fig.10.Radiation intensity I vs.impact velocity V0.

    Fig.11.Spectral analysis of impact flash for 6061 aluminum at V0=1.43 km/s.Nitrogen and oxygen spectral lines were identified.The ICCD gate width was 2μs.

    Fig.9.Spectrum for impact flash at V0=3.12 km/s.The ICCD gate width was 4μs.

    Fig.12.Spectral analysis of impact flash for pure aluminum(99.9999%)at V0=1.40 km/s.Nitrogen and oxygen spectral lines were identified.The ICCD gate width was 2μs.

    The spectral analysis of the impact flash for V0<2 km/s was similar.When V0>2 km/s,the continuous spectrum representing high-temperature condensed matter increases in intensity with V0.Thus,it is likely to be the result of a shock temperature increase of the collision material.However,visible light emission of multiphase materials is very complex and is the result of photon,electron,and ion interactions.

    Fig.13.Spectral analysis of impact flash for 45 steel V0=1.30 km/s.Nitrogen and oxygen spectral lines were identified.The ICCD gate width was 2μs.

    4.Discussion

    4.1.Impact flash emission source and its characteristics

    Impact flash is an interdisciplinary problem of collision dynamics and light emission,and accurate measurements of impact flashes are critical.The experimental results here revealed impact flashes that are more complex than previously thought.We first define the light source,and then summarize its characteristics and suggest further research.

    We define the illuminating area formed by the collision as a light source,as shown in Fig.4(a and b)at t=0.In the high-speed photographs,the emitting area is continuously expanding;thus,it is a dynamic light source,as shown in Fig.14,the ds is the area element of photoelectric sensor(observation point),and the L0is the original distance between light source and observation point.The luminous zone expands faster at higher V0.At t=16μs in Fig.4(b),the light source exceeded the field of view of the camera(about 80×80 mm)and overexposed the image.

    Fig.14.Dynamic light source at different time.(a)Original time;(b)t=tn.

    The dynamic feature of the light source has two main effects on the impact flash test.First,from a fixed observation point,the solid angle changes with time,which greatly increases the difficulty of performing absolute measurements.Second,from different observation points,the observed impact flash is constant.In Table 1,the“overexposure”state recorded by the camera was significantly longer than the pulse width recorded by the photoelectric sensor.This was most likely because the observation point of the highspeed camera was much farther from the collision point than was the sensor.

    The spectral characteristics of the various V0suggest that the impact flash was visible light emission from multiphase materials,which is the second feature of the impact flash light source.Gaseous atoms emit sharp line spectra,gaseous molecules emit band spectra,and condensed matter(liquid sand solids)emit continuous spectra[32].For V0<2 km/s,the spectra mainly consisted of atomic lines.For V0>2 km/s,in addition to a large number of line spectra,a continuous thermal radiation background appeared,indicating a mixture of gaseous and condensed states.This makes it very difficult to determine the temperature of the condensed phase from the impact flash.Previous results suggested that the temperature of the impact flash determined with a multichannel pyrometer was significantly higher than the shock temperature rise of the collision materials[22].

    4.2.The mechanism of impact flash

    When V0<2 km/s,the atomic spectral lines were not aluminum,which is significantly in contrast to previous results.In the spherical aluminum projectile impact flash at V0=5-6 km/s,there were two characteristic atomic aluminum lines at 396.15 nm and 394.40 nm[26],derived from two excited states.It was therefore generally believed that the collision material exhibits gasification and luminescence[17,21].The energy transport process is depicted in Fig.15(a).From shock wave compression to the impact flash,the energy carriers are all based on aluminum.

    Fig.15.Energy transport processes and carrier.(a)Previous understanding;(b)New mechanism.

    However,our spectral analysis results show that excited and ionized ambient gases are the main sources of the intense flash at the moment of collision at V0<2 km/s.The energy carriers of this new impact flash mechanism are shown in Fig.15(b),where the energy transport process is unknown.It is not clear whether the new mechanism persists or dominates when V0>2 km/s.The first excited states of nitrogen and oxygen atoms are at 10.88 eV and 9.48 eV,respectively;whereas,that of aluminum atoms is at 3.13 eV.Thus,nitrogen and oxygen atoms are more difficult to excite than aluminum atoms.Hence,the energy that causes the ambient gas to emit light is not from the shock temperature rise.The new impact flash mechanism occurs when the shock residual temperature is not sufficient to vaporize aluminum.

    4.3.Spectra from V0<2 km/s

    When V0<2 km/s,no continuous spectra were observed.According to black-body thermal physics,condensed matter at a temperature greater than 0 K will radiate a spectrum of electromagnetic energy that correlates with the temperature[33].At low temperatures,the light may not be visible to the human eye,nor detected with optical devices such as ICCDs.Here,we roughly estimate the temperature threshold of condensed matter with the ICCD spectrometer.

    The sensitivity of the photosensitive ICCD element is about 1×10-6Lx(photometry unit),that corresponds to Ee=4.6×10-9w/m2/nm(Irradiance of the photosensitive element)monochromatic emission at 555 nm[34].The original distance L0(Fig.14(a))between the optical fiber probe and the collision point is 100 mm.The relationship between expanding velocity(Ve)of emitting area and V0is Ve=αV0,and theαis about 4-4.7 with different V0according to Fig.4.When V0=2 km/s,theαis 4.4.Then Rn=Vetn+R0=20.6 mm at t=2μs is calculated,where the R0is the radius of the projectile.

    As luminous zone expanding(Fig.14(b)),the relationship between the radiance of light source(Me)and Eeof ds is:

    whereλis the emission wavelength,ηis the radiation efficiency,c1=3.742×10-16w?m2is the first radiation constant,and c2=1.439×10-2m K is the second radiation constant.

    The calculation results for the radiance and various temperatures of the condensed matter source are plotted in Fig.16,where the vertical axis is Me,the horizontal axis is the wavelength,theis 1.12×10-7w/m2/nm and theηis 0.8.Emission at 555 nm and radiance equal to or greater thancan only be detected by the ICCD spectrometer for condensed matter temperatures above 822 K.The calculation results have two meanings.First,when V0<2 km/s,the thermal radiation could not be detected by ICCD spectrometer so the relationship between radiation intensity and impact velocity is not continuous.Second,if this assumption is true,then the calculated impact flash intensity for V0<2 km/s may need to be corrected,depending on the application requirements.But this does not affect the conclusion regarding the flash mechanism.Also,it can be inferred that the impact flash for V0<2 km/s has little relationship with shock temperature rise,because the aluminium material could not emit line spectrum below 822 K.

    Fig.16.Results of radiance and temperature calculations of condensed matter emission.

    5.Conclusions

    Here,a series of impact flash experiments for a flat-head longrod projectile impacting a thin target plate were performed.A highspeed camera,a photoelectric sensor,and an ICCD spectrometer were used to characterize the flash and mechanism for V0=1.3-3.2 km/s.There were two characteristics of the impact flash light source.First,it was dynamic,where its geometric size continuously increased with time.Second,the light source was a mixture of gaseous and condensed materials.These two features make it challenging to accurately identify and evaluate the flash source with photometry.Conclusions were drawn as follows.

    The impact flash intensity I of the long-rod projectile exhibited a pulse characteristic with time.The duration increased with V0.I was a power function of V0,i.e.I∝.For V0>2 km/s,n=4.41,which is close to previous results for spherical aluminum projectiles at V0=3-7 km/s.However,for V0<2 km/s,n=2.21,indicating less of an increase in flash intensity with V0.The intense flash durations were also different in that the value for V0<2 km/s was an order of magnitude less than that for V0>2 km/s.

    When V0>2 km/s,continuous spectra were observed and increased in intensity with V0.However,when V0<2 km/s,the intense flash spectra mainly consisted of atomic lines.Spectral analysis indicated that the lines derived from the excitation and ionization of ambient gases near the impact surface.In contrast to shock-induced vaporization luminescence of the collision materials,the finding here indicates a new mechanism for impact flashes.Although the energy transport process is uncertain,it has provided a new direction for future studies.Because the new mechanism was also observed in impact flash experiments with pure aluminum and 45 steel,ambient gas emission may be a common feature of hypersonic impact flashes of colliding metallic materials.

    Declaration of competing interest

    The authors declare no conflict of interests.

    Acknowledgments

    The authors would like to thank Mrs Li Chen,Mr Siyuan Ren,Mr Cheng Shang,Mr.Weijiang Hao,Mr.Yan Liu,Mr.Yaoxin Wang and others for their help in conducting the two-stage light-gas gun tests.This work was supported by the National Key R&D Program of China(Grant No.2016YFC0801204)and the National Program on Key Basic Research Project of China(973 Program,Grant No.613312).

    国产成人午夜福利电影在线观看| 久久精品夜色国产| 国产精品久久久久久精品电影小说| a 毛片基地| 久久人人爽人人爽人人片va| 2018国产大陆天天弄谢| 亚洲美女视频黄频| 极品人妻少妇av视频| 99香蕉大伊视频| 伊人亚洲综合成人网| 精品熟女少妇av免费看| 欧美日韩一区二区视频在线观看视频在线| 纯流量卡能插随身wifi吗| 一区二区三区乱码不卡18| 午夜日本视频在线| 纵有疾风起免费观看全集完整版| 国产免费福利视频在线观看| 亚洲精品aⅴ在线观看| 亚洲av日韩在线播放| 精品熟女少妇av免费看| 国产永久视频网站| 最黄视频免费看| 一本色道久久久久久精品综合| 国产日韩一区二区三区精品不卡| 高清不卡的av网站| 在线观看免费高清a一片| 又黄又粗又硬又大视频| 日韩 亚洲 欧美在线| 性色av一级| 久久精品国产亚洲av天美| 国产一区有黄有色的免费视频| 捣出白浆h1v1| 日韩精品免费视频一区二区三区 | 久久久久精品人妻al黑| 十八禁网站网址无遮挡| 精品酒店卫生间| 亚洲精品久久成人aⅴ小说| 亚洲精品国产av成人精品| 成人国语在线视频| 欧美日韩一区二区视频在线观看视频在线| 丝袜喷水一区| 一二三四在线观看免费中文在 | 97在线视频观看| 永久免费av网站大全| 男女无遮挡免费网站观看| 丝瓜视频免费看黄片| 国产色爽女视频免费观看| 久久久精品免费免费高清| www.av在线官网国产| 国产女主播在线喷水免费视频网站| 日本色播在线视频| 日韩中文字幕视频在线看片| 综合色丁香网| 成年人午夜在线观看视频| 激情视频va一区二区三区| 国产色爽女视频免费观看| 成年人免费黄色播放视频| 欧美精品人与动牲交sv欧美| 精品一区二区三区四区五区乱码 | 天美传媒精品一区二区| 欧美xxxx性猛交bbbb| 午夜精品国产一区二区电影| 亚洲精品第二区| 午夜福利影视在线免费观看| 久久久久国产精品人妻一区二区| 日本91视频免费播放| 国产伦理片在线播放av一区| 精品酒店卫生间| 亚洲av欧美aⅴ国产| 成人国产麻豆网| 免费av中文字幕在线| 精品卡一卡二卡四卡免费| 中文字幕人妻熟女乱码| 亚洲av男天堂| 精品第一国产精品| 亚洲精品国产av成人精品| 亚洲av综合色区一区| 欧美97在线视频| 国产在线免费精品| 国产精品偷伦视频观看了| 精品久久蜜臀av无| 观看美女的网站| 一级毛片我不卡| 曰老女人黄片| 国产又色又爽无遮挡免| 少妇被粗大猛烈的视频| 久久韩国三级中文字幕| 国产高清三级在线| 国产精品国产三级专区第一集| 制服人妻中文乱码| 国产爽快片一区二区三区| 亚洲成色77777| 少妇人妻久久综合中文| 亚洲国产精品成人久久小说| 国产精品国产av在线观看| 免费高清在线观看视频在线观看| 天美传媒精品一区二区| 丁香六月天网| 亚洲精品自拍成人| 亚洲 欧美一区二区三区| 久久久久久人人人人人| 欧美xxⅹ黑人| 亚洲国产欧美在线一区| 高清不卡的av网站| 国产精品嫩草影院av在线观看| 国产精品国产三级国产专区5o| 亚洲五月色婷婷综合| 在线亚洲精品国产二区图片欧美| 国产在线免费精品| 草草在线视频免费看| 好男人视频免费观看在线| 国产深夜福利视频在线观看| 在线观看免费视频网站a站| 国产精品一区二区在线观看99| 久久99蜜桃精品久久| 精品一区二区三区四区五区乱码 | 亚洲精品久久久久久婷婷小说| 18禁裸乳无遮挡动漫免费视频| 在线观看免费高清a一片| a级毛色黄片| 亚洲第一av免费看| 国产在线免费精品| 亚洲精品,欧美精品| 国产国拍精品亚洲av在线观看| 热re99久久国产66热| 久久精品久久久久久久性| 日韩电影二区| 亚洲第一区二区三区不卡| 精品一区二区免费观看| 日韩免费高清中文字幕av| 巨乳人妻的诱惑在线观看| 精品人妻偷拍中文字幕| 91精品国产国语对白视频| videos熟女内射| 国产免费一级a男人的天堂| 多毛熟女@视频| 欧美3d第一页| 国产男女超爽视频在线观看| 久久人人爽人人爽人人片va| 国产无遮挡羞羞视频在线观看| 性高湖久久久久久久久免费观看| 爱豆传媒免费全集在线观看| 99热6这里只有精品| 免费看光身美女| 国产成人91sexporn| 在线 av 中文字幕| 欧美人与善性xxx| 亚洲国产av新网站| 午夜激情av网站| 欧美激情国产日韩精品一区| 少妇人妻久久综合中文| 天堂俺去俺来也www色官网| 国产精品.久久久| 久久久久视频综合| 少妇人妻精品综合一区二区| 久久久久国产网址| 亚洲成国产人片在线观看| 一二三四在线观看免费中文在 | 国产亚洲最大av| 999精品在线视频| 成人亚洲精品一区在线观看| av又黄又爽大尺度在线免费看| 一边摸一边做爽爽视频免费| 亚洲成色77777| 高清av免费在线| 欧美3d第一页| 天天操日日干夜夜撸| 黑人欧美特级aaaaaa片| 一个人免费看片子| 哪个播放器可以免费观看大片| 最后的刺客免费高清国语| 一区二区日韩欧美中文字幕 | 韩国av在线不卡| 精品酒店卫生间| 亚洲综合色网址| 国产69精品久久久久777片| 亚洲经典国产精华液单| 九色成人免费人妻av| 99九九在线精品视频| 三上悠亚av全集在线观看| 亚洲美女搞黄在线观看| 少妇的逼水好多| 国产黄频视频在线观看| 国产国拍精品亚洲av在线观看| 毛片一级片免费看久久久久| 一边摸一边做爽爽视频免费| av黄色大香蕉| av有码第一页| 99热国产这里只有精品6| 久久人人爽av亚洲精品天堂| 一二三四在线观看免费中文在 | 国产极品天堂在线| 一级毛片电影观看| 水蜜桃什么品种好| 午夜免费观看性视频| 亚洲一码二码三码区别大吗| 一本色道久久久久久精品综合| 最近手机中文字幕大全| 国产亚洲av片在线观看秒播厂| 香蕉丝袜av| 精品久久蜜臀av无| 色5月婷婷丁香| 日本与韩国留学比较| 国产免费一区二区三区四区乱码| 国产不卡av网站在线观看| 99热国产这里只有精品6| 国产片特级美女逼逼视频| 精品人妻熟女毛片av久久网站| 国产精品国产三级国产av玫瑰| a级毛色黄片| 99视频精品全部免费 在线| av天堂久久9| videos熟女内射| 一区二区三区精品91| 午夜福利影视在线免费观看| 亚洲在久久综合| 黑人高潮一二区| 妹子高潮喷水视频| 久久免费观看电影| 亚洲av电影在线观看一区二区三区| 精品国产国语对白av| 一区在线观看完整版| 91在线精品国自产拍蜜月| 啦啦啦视频在线资源免费观看| 亚洲av在线观看美女高潮| 纵有疾风起免费观看全集完整版| av免费观看日本| 久久韩国三级中文字幕| 国产乱来视频区| 99国产精品免费福利视频| 不卡视频在线观看欧美| 久久久亚洲精品成人影院| 久久人人爽人人爽人人片va| 三级国产精品片| 捣出白浆h1v1| av免费在线看不卡| 亚洲国产精品成人久久小说| 亚洲三级黄色毛片| 日本免费在线观看一区| 男女国产视频网站| 超碰97精品在线观看| 久久久久精品人妻al黑| 美女国产视频在线观看| 人人澡人人妻人| 免费黄网站久久成人精品| 秋霞伦理黄片| 99国产综合亚洲精品| 美女主播在线视频| 夫妻午夜视频| 青春草国产在线视频| 久久久国产欧美日韩av| 欧美日韩一区二区视频在线观看视频在线| 这个男人来自地球电影免费观看 | 高清不卡的av网站| 我要看黄色一级片免费的| 麻豆乱淫一区二区| www.av在线官网国产| 亚洲av中文av极速乱| 国产av国产精品国产| 97精品久久久久久久久久精品| 男人添女人高潮全过程视频| a级毛片在线看网站| 高清av免费在线| 免费日韩欧美在线观看| 国产精品人妻久久久久久| 国产高清不卡午夜福利| 精品一区二区三区四区五区乱码 | 99热网站在线观看| 男女国产视频网站| 欧美日韩亚洲高清精品| 男人添女人高潮全过程视频| 岛国毛片在线播放| 精品国产一区二区久久| 亚洲国产最新在线播放| 妹子高潮喷水视频| 人体艺术视频欧美日本| 在线观看国产h片| 又黄又粗又硬又大视频| kizo精华| 精品国产一区二区三区久久久樱花| 一本色道久久久久久精品综合| 美女xxoo啪啪120秒动态图| 亚洲国产成人一精品久久久| 亚洲,欧美精品.| 国产成人午夜福利电影在线观看| av播播在线观看一区| 成人18禁高潮啪啪吃奶动态图| 精品少妇内射三级| 成年动漫av网址| 大香蕉久久成人网| 免费久久久久久久精品成人欧美视频 | 丝袜人妻中文字幕| 国产av国产精品国产| 午夜老司机福利剧场| 国产高清三级在线| 老司机影院毛片| 日韩不卡一区二区三区视频在线| av电影中文网址| 精品一区二区三卡| 久久国产精品大桥未久av| 大码成人一级视频| 久久韩国三级中文字幕| 三上悠亚av全集在线观看| 日本黄色日本黄色录像| 久久久久精品久久久久真实原创| 久久午夜福利片| 秋霞伦理黄片| 日韩av在线免费看完整版不卡| 黄色毛片三级朝国网站| 亚洲精品一区蜜桃| a级毛片在线看网站| 久久久国产一区二区| 在线观看www视频免费| videos熟女内射| 亚洲国产精品一区二区三区在线| 日韩人妻精品一区2区三区| 精品第一国产精品| 亚洲精品中文字幕在线视频| 久久久国产精品麻豆| 国产精品.久久久| 黄网站色视频无遮挡免费观看| 久久久a久久爽久久v久久| 免费女性裸体啪啪无遮挡网站| 夫妻午夜视频| 国产欧美日韩一区二区三区在线| 日本-黄色视频高清免费观看| 伦理电影免费视频| 午夜精品国产一区二区电影| 日韩视频在线欧美| 飞空精品影院首页| 91精品三级在线观看| 两性夫妻黄色片 | 人人澡人人妻人| 三级国产精品片| 黄色配什么色好看| 日韩中文字幕视频在线看片| 国产精品嫩草影院av在线观看| 久久国产精品男人的天堂亚洲 | 日日爽夜夜爽网站| 极品少妇高潮喷水抽搐| 一级,二级,三级黄色视频| 国产成人免费无遮挡视频| 精品一区二区三卡| 观看av在线不卡| 国内精品宾馆在线| 人妻一区二区av| 美女国产视频在线观看| 99热全是精品| 精品国产一区二区三区久久久樱花| 激情五月婷婷亚洲| 免费黄网站久久成人精品| av网站免费在线观看视频| 久久97久久精品| 精品视频人人做人人爽| 亚洲天堂av无毛| 中文乱码字字幕精品一区二区三区| 大香蕉久久网| 精品一区在线观看国产| 国产午夜精品一二区理论片| 欧美精品av麻豆av| 亚洲,欧美,日韩| 国产日韩一区二区三区精品不卡| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 激情五月婷婷亚洲| 尾随美女入室| 久久精品久久精品一区二区三区| www.av在线官网国产| 亚洲精品av麻豆狂野| 国产在线视频一区二区| 美女福利国产在线| 色94色欧美一区二区| 美女内射精品一级片tv| 在线天堂最新版资源| 美国免费a级毛片| 色网站视频免费| 天堂8中文在线网| 亚洲性久久影院| 日韩成人av中文字幕在线观看| 尾随美女入室| 高清不卡的av网站| 最近中文字幕2019免费版| 国语对白做爰xxxⅹ性视频网站| 久久国产精品大桥未久av| av国产久精品久网站免费入址| 90打野战视频偷拍视频| 中文字幕精品免费在线观看视频 | 亚洲精品第二区| 汤姆久久久久久久影院中文字幕| 亚洲av综合色区一区| 91aial.com中文字幕在线观看| 婷婷成人精品国产| 精品国产一区二区久久| tube8黄色片| 五月玫瑰六月丁香| 少妇被粗大的猛进出69影院 | 久久狼人影院| 亚洲欧洲精品一区二区精品久久久 | 丝袜人妻中文字幕| 午夜福利影视在线免费观看| 成人亚洲精品一区在线观看| 丰满迷人的少妇在线观看| 国产色婷婷99| 亚洲av.av天堂| 亚洲av在线观看美女高潮| 亚洲熟女精品中文字幕| 午夜激情av网站| 永久网站在线| 99久久综合免费| 色94色欧美一区二区| 卡戴珊不雅视频在线播放| 亚洲性久久影院| 欧美国产精品一级二级三级| 亚洲第一av免费看| 99久久中文字幕三级久久日本| 久久精品国产亚洲av天美| 国产一区二区三区av在线| 99热6这里只有精品| 成人毛片60女人毛片免费| 欧美国产精品va在线观看不卡| 插逼视频在线观看| 亚洲婷婷狠狠爱综合网| 国产毛片在线视频| 丰满迷人的少妇在线观看| 国产av一区二区精品久久| 亚洲精品色激情综合| 亚洲少妇的诱惑av| 26uuu在线亚洲综合色| 在线观看www视频免费| 2022亚洲国产成人精品| 97人妻天天添夜夜摸| 国产成人精品久久久久久| 精品少妇久久久久久888优播| 亚洲欧美成人精品一区二区| 国国产精品蜜臀av免费| 丁香六月天网| 国产麻豆69| 女人被躁到高潮嗷嗷叫费观| 亚洲精品日韩在线中文字幕| 成年人午夜在线观看视频| 免费观看在线日韩| 亚洲成人手机| 伦理电影大哥的女人| 国产成人av激情在线播放| 欧美变态另类bdsm刘玥| 色哟哟·www| 免费人成在线观看视频色| 色吧在线观看| 久久精品国产鲁丝片午夜精品| 国产精品秋霞免费鲁丝片| 午夜精品国产一区二区电影| 哪个播放器可以免费观看大片| 亚洲精品中文字幕在线视频| 欧美精品一区二区免费开放| 久久久久精品性色| 最近最新中文字幕大全免费视频 | 高清在线视频一区二区三区| 亚洲天堂av无毛| 欧美成人精品欧美一级黄| 婷婷色综合大香蕉| 美女视频免费永久观看网站| 亚洲美女黄色视频免费看| 菩萨蛮人人尽说江南好唐韦庄| 性色avwww在线观看| 大香蕉久久网| 黄色怎么调成土黄色| 国产又色又爽无遮挡免| 日本猛色少妇xxxxx猛交久久| 99久久人妻综合| 美女主播在线视频| 国产成人欧美| 蜜臀久久99精品久久宅男| 黄色视频在线播放观看不卡| 在线看a的网站| 人妻人人澡人人爽人人| 免费av不卡在线播放| 久久国产亚洲av麻豆专区| 久久鲁丝午夜福利片| 黄色怎么调成土黄色| 91精品三级在线观看| 18禁在线无遮挡免费观看视频| 男女免费视频国产| 国产日韩欧美视频二区| 国产精品秋霞免费鲁丝片| 在线观看人妻少妇| 一本色道久久久久久精品综合| 亚洲国产毛片av蜜桃av| 成人国产麻豆网| 精品国产乱码久久久久久小说| 一边摸一边做爽爽视频免费| 男女啪啪激烈高潮av片| 午夜激情av网站| 亚洲一码二码三码区别大吗| 午夜福利影视在线免费观看| 久久精品人人爽人人爽视色| 精品少妇内射三级| 人人妻人人澡人人爽人人夜夜| 国产成人免费无遮挡视频| 国产欧美亚洲国产| 国精品久久久久久国模美| 日韩,欧美,国产一区二区三区| 校园人妻丝袜中文字幕| 免费看光身美女| 啦啦啦中文免费视频观看日本| 精品少妇黑人巨大在线播放| 免费大片18禁| 欧美+日韩+精品| 久久人人爽av亚洲精品天堂| 国产精品不卡视频一区二区| 制服人妻中文乱码| 人人妻人人爽人人添夜夜欢视频| 黄色怎么调成土黄色| 一级毛片黄色毛片免费观看视频| 三级国产精品片| 美女国产视频在线观看| 亚洲欧美成人精品一区二区| 欧美日韩国产mv在线观看视频| 亚洲一区二区三区欧美精品| 亚洲伊人色综图| 精品一区在线观看国产| 老司机影院毛片| 亚洲五月色婷婷综合| 国产精品秋霞免费鲁丝片| 精品人妻在线不人妻| 亚洲欧美日韩卡通动漫| 亚洲av欧美aⅴ国产| 国产欧美日韩一区二区三区在线| 亚洲精品乱久久久久久| 一级爰片在线观看| 精品久久国产蜜桃| 亚洲人成网站在线观看播放| 国产成人av激情在线播放| 欧美日本中文国产一区发布| √禁漫天堂资源中文www| 精品卡一卡二卡四卡免费| 久久鲁丝午夜福利片| 欧美日韩视频高清一区二区三区二| 99国产精品免费福利视频| xxxhd国产人妻xxx| 色94色欧美一区二区| 亚洲av免费高清在线观看| 中文天堂在线官网| 天堂俺去俺来也www色官网| 波多野结衣一区麻豆| 老熟女久久久| 国产男人的电影天堂91| 国产精品人妻久久久久久| 中文字幕av电影在线播放| av有码第一页| 精品第一国产精品| 国产毛片在线视频| 日日爽夜夜爽网站| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 性色av一级| 人人妻人人爽人人添夜夜欢视频| 国产男女内射视频| 国产高清三级在线| 五月玫瑰六月丁香| 久久精品aⅴ一区二区三区四区 | 免费观看a级毛片全部| 久久久久久久国产电影| av免费在线看不卡| 成年女人在线观看亚洲视频| 亚洲欧洲精品一区二区精品久久久 | 亚洲人成网站在线观看播放| 国产精品一区二区在线不卡| 免费大片黄手机在线观看| 一级毛片我不卡| 高清黄色对白视频在线免费看| 麻豆精品久久久久久蜜桃| 99热国产这里只有精品6| 久久精品国产鲁丝片午夜精品| 亚洲一级一片aⅴ在线观看| 男人添女人高潮全过程视频| av一本久久久久| 欧美人与善性xxx| 美女内射精品一级片tv| 亚洲欧美精品自产自拍| 狂野欧美激情性bbbbbb| 18在线观看网站| 三上悠亚av全集在线观看| av又黄又爽大尺度在线免费看| 18禁动态无遮挡网站| 最新中文字幕久久久久| 新久久久久国产一级毛片| 日韩中文字幕视频在线看片| tube8黄色片| 久久综合国产亚洲精品| 日韩熟女老妇一区二区性免费视频| 亚洲婷婷狠狠爱综合网| 亚洲av日韩在线播放| 亚洲欧美中文字幕日韩二区| 男女下面插进去视频免费观看 | 国产精品 国内视频| 久久99热6这里只有精品| 美女主播在线视频| 国产日韩欧美亚洲二区| 激情五月婷婷亚洲| 90打野战视频偷拍视频| av女优亚洲男人天堂| 日韩av在线免费看完整版不卡| 18禁动态无遮挡网站| 免费播放大片免费观看视频在线观看| 午夜福利视频在线观看免费| 久久热在线av| 日本免费在线观看一区| 人妻少妇偷人精品九色| 成人黄色视频免费在线看| 国产亚洲精品第一综合不卡 | 欧美变态另类bdsm刘玥| 五月玫瑰六月丁香| 国产精品久久久久成人av| 日韩人妻精品一区2区三区| 欧美丝袜亚洲另类| 国产一区二区在线观看av| 欧美97在线视频|