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

    Anechoic Effect of Coating Covered on Target Submerged in Typical Shallow-water Waveguides

    2012-09-22 07:15:52CHENYanTANGWeilinFANWeiFANJun
    船舶力學(xué) 2012年3期

    CHEN Yan,TANG Wei-lin,FAN Wei,FAN Jun

    (School of Naval Architecture,Ocean and Civil Engineering,Shanghai Jiao Tong University,Shanghai 200240,China)

    1 Introduction

    Anechoic coating is one of the most important stealth measures of modern submarines.The reduction of the Target Strength(TS)in free space was usually used to indicate the anechoic effect of coating in previous estimation and prediction research.This can lead two problems as follows.First,what is the relation between the reduction of the TS in free space and the real echo reduction in various shallow water waveguide?Second,what the difference in the echo reduction of target with varied anechoic coatings in shallow-water channel is?The problems mentioned above are caused by different formation mechanisms of echo from objects in free space and in shallow-water environment.The ray acoustic method shows the echo formation mechanism clearly.In free space,the echo is the combination of the direct in-

    cident wave and the direct reflection wave,as shown in Fig.1(a).But in shallow water environment,the multi-path effect has a great impact on the propagation.The first incident wave is reflected by boundaries and turns into infinite incident waves,and each incident wave is then reflected by the target and the boundaries and turns into infinite reflection waves.That means the scattering in shallow-water waveguide involves the double multi-path effects.Fig.1(b)and(c)give the picture of rays reflected by surface and bottom for only one time.In fact,there are infinite incident and reflection rays.The ray angles relative to reflection areas of target are different from each other.

    Fig.1 Rays in free space and in shallow-water waveguide

    In 1987,Ingenito[1]derived the scattered field expression in terms of the normal modes and plane-wave scattering functions in free space,which was a groundwork for wave acoustic method in shallow-water waveguide.Makris promoted Ingenito’s method and put forward the spectral approach in the article about the detection of a submerged object insonified by surface noise in an ocean waveguide in 1994.The advantage of this method is that the waveguide Green’s function is decomposed into plane wave spectrum,including both discrete spectrum(normal mode)and continuous spectrums,so as to enable to calculate the scattered field much closer to the object.After that,the team of Makris did further studies of spectral approach.They developed a unified model for reverberation and submerged target scattering in a stratified medium,and discussed the validity of the sonar equation for object scattering in a shallow water waveguide[3-6].In recent years,we developed a simple and practical ray acoustic method for sonar engineering application[7-8].Since the sound propagation is calculated by the ray-tracing technique and the scattering is expressed through Kirchhoff approximation,the method can be used to deal with the scattering of complex targets in shallow-water waveguide with arbitrary velocity profile in principle.

    The main characteristic of coating covered on the target surface is the reflection coefficient,which changes with frequency and incident angle.In free space,the echo is the direct incident-reflected wave,so the coating only works on a reflection area illuminated by the direct incident wave.But in shallow water waveguide,the echo is the summation of reflection waves from various combinations of rays and each combination has the varied incident and reflection angles,which causes the different anechoic effect of coating.Especially in the situation of the negative gradient velocity profile as shown in Fig.1(c),when the horizontal range is larger than a skip distance,there is no direct echo from target.That will lead to larger difference of the action mechanism of coating between in free space and in shallow-water environment.

    A numerical procedure for calculating echoes from objects with anechoic coating submerged in shallow-water waveguide based on ray-tracing technique and physical acoustic method is developed.The coating has no effect on propagation of the waveguide but changes the target scattering function,the Kirchhoff approximation is extended to treat the scattering of objects with anechoic coating.Taking the Benchmark submarine model and two supposed coatings with different reflection coefficient,we discussed the echo characteristics of target with rigid and coated surface in two typical shallow water environment:weak positive gradient profile(W.P.G-profile)and negative gradient profile(N.G-profile).The analysis of the influence of velocity profile on the anechoic effect of coating is given and the difference of the anechoic effect in free space and shallow-water waveguide is obtained.

    2 Numerical approach

    The references[7-8]gave the numerical approach of the scattering of target submerged in shallow-water waveguide.The scattered field is the summation of all the monostatic and bistatic scattered fields.

    where(φs)l=mis the monostatic scattered field with the same incident ray and reflection ray,(φs)l≠mis the bistatic scattered field with different incident ray and reflection ray.It is only under the condition of waveguide environment,the bistatic scattering can come back to the source and contribute to the echo.Using the ray-tracing technique,the monostatic echo is expressed as

    where ΔΦz(mì)lis the path difference from arrival pointr,()z of incident ray on target surface to reference pointαzlis the included angle between incident ray and the normalat the arrival point of surface.Bl,lis the bright region for lth incident-reflection rays.The bistatic scattered field is:

    where Sl,mis the bistatic scattering function.If the target is rigid,then we get

    where Bl,mis the common bright region for both lth incident ray and mth reflection ray,and Bl,m=Bm,l,Sl,m=Sm,land (φs)l,m= (φs)m,lcan be easily obtained.

    The ray-tracing technique is used to calculate the transmission factors of incident and reflection rays and the Kirchhoff approximation is applied to obtain the monostatic or bistatic scattering function from each combination of rays.The echo in shallow water waveguide is the summation of all the monostatic and bistatic scattering.The summation process is convergent because that only limited rays work on the echo in the far field.

    In this paper,we promoted the numerical procedure in reference[8]so as to use it for target with coating.The real effect of coating in two typical ocean environments is discussed.When calculating the object with coating in shallow water environment,the transmission characteristic of waveguide is changeless but the target scattering function changes relatively.The reflection coefficient V()α is used to present the anechoic behavior of coating,where α is the incident angle relative to the surface normal.According to the physical acoustic method for nonrigid surface[9],the scattering function is related to the reflection coefficient at the incident angle.The monostatic scattering function is modified as

    The bistatic scattering function is modified as

    It is necessary to indicate that

    When the incident ray and reflection ray swap,the bistatic scattering function changes,that means (φs)l,m≠ (φs)m,l.Even then the bistatic scattering can be composed as

    When dealing with each component of the scattered field,the transmission and scattering can be decoupled.However,this does not mean that the total scattered field can be decoupled,since the summation of each scattered component can not define the total transmission factor and total scattering function.In traditional sonar equation,the waveguide transmission loss and target scattering are decoupled artificially and the echo level is defined as follows:

    where SL is the source level;TL is the transmission loss in the shallow-water waveguide,which usually indicates the transmission loss form the source to the center of the target;TS is the Target Strength in free space.But the echo level can not be represented as in formula(10)in the shallow water environment because the total transmission factor and total scattering function can not be decoupled.According to formula(2),the Target Strength in free space is expressed as

    There has no bistatic scattering in the Eq.(11).In order to evaluate the error due to traditional sonar equation in shallow-water waveguide,a physical quantity-Equivalent Target Strength ETS is defined as[7]:

    In this paper,the TS in free space,the ETS in waveguide and the TS reduction and ETS reduction after covering various anechoic coatings are analyzed.The ETS reduction can be used to evaluate the anechoic effect of coating in shallow-water waveguide because in the same shallow-water environment the echo reduction equals the ETS reduction.We discussed the influence of velocity profile on the anechoic effect of coating and the difference of the anechoic effect in free space and in shallow-water waveguide.

    3 Results and analysis

    The reflection coefficient of two supposed coatings as a function of incident angle at the frequency of 5 kHz is shown in Fig.2.It can be seen that reflection coefficients of the two coatings are equal with each other at the zero incident angle when the ray is vertical incidence.But the increasing rate of the reflection coefficient of the coating E1 is larger than the coating E2 when incident angle getting larger.The anechoic behavior of the coating E1 is worse than the coating E2.

    The shallow-water waveguide is range-independent with 80 m deep as illustrated in Fig.3.The transducer is in the depth of 10 m.The target is in the depth of 40 m.The frequency is 5 kHz and the absorption loss of sea water is 0.3862 dB/km.The sediment layer is clay with constant density of 1420 kg/m3and sound speed of 1510 m/s.The attenuation of the sediment is 7.5×10-5dB/m·Hz.Considering the sea surface is release.

    Fig.2 Reflection coefficient of coatings E1 and E2

    Fig.3 Positions of transducer and target in shallow-water waveguide

    In the winter and summer shallow-water waveguide,the typical velocity profile and the corresponding eigen rays from transducer to the target with the range of 1 km are shown in Fig.4.

    Fig.4 Typical velocity profile and eigen rays in winter and summer shallow-water waveguide

    The isothermal layer in winter leads to W.P.G-profile,which makes the rays slowly bending upward.There exists thermocline in summer waveguide,which leads to N.G-profile and makes the eigen rays so bending downward that there is no direct ray when the distance over the skip distance of limiting ray.

    3.1 The echo reduction of cylinder with coating in shallow water waveguide

    The references[7-8]pointed out that when a rigid cylinder is placed at cross direction with its axis of symmetry perpendicular to the range and depth coordinates,the ETS is less influenced by shallow water waveguide;when the cylinder is set at longitudinal direction with its axis of symmetry parallel to the depth coordinates,the shallow water environment has great effect on the ETS.The cylinder set at longitudinal direction is a typical target which cannot use the traditional sonar equation to calculate its echo level in shallow water waveguide.Therefore,we analyze the echo characteristics from the cylinder with coating firstly.The cylinder is 5 m long with radius of 1m.Its position used in the following example is shown in Fig.3.The range varies from 500 m to 2000 m.At the frequency of 5 kHz,the ratio of the cylinder length to wave length is about 16.7.

    Fig.5gives the anechoic effects of the two coatings when the cylinder is set at cross direction.Both the TS reduction(△TS)in free space and the ETS reduction(△ETS)in shallow water waveguide caused by the two coatings are close with each other.The anechoic effects caused by coating E2 is larger than that by coating E1 only 0.5 dB in free space and 1 dB in shallow water waveguide.It means that for a cylinder with anechoic coating set at cross direction the shallow-water environments has little effect on the anechoic effects and the advantage of the coating E2 is enlarged slightly in waveguide.

    Fig.5 The anechoic effects of the two coatings when cylinder setting at cross direction

    Fig.6gives the anechoic effects of the two coatings when the cylinder is set at longitudinal direction.The TS reductions caused by coatings E1 and E2 are almost the same in free space compared with Fig.5.But in shallow-water environment,the waveguide make the ETS reduction of the coating E1 worse,especially in N.G-profile,there is a difference of 1.5 dB between the coating E1 and the coating E2.

    The main reason that brings out the difference between coatings E1 and E2 are:①The echo of a cylinder set at longitudinal direction is sensitive to the incident angle,it just so happened that the multi-path effect leads to larger incident and reflection angle,as shown in Fig.1(b)and(c);②The reflection coefficient of the coating E1 increases rapidly in small angle range(0°to 30°approximately),which makes the anechoic effect of the coating E1 behave worse at larger incident angle;③In the situation of N.G-profile in shallow water waveguide,there is no direct eigen ray.All the incident-reflection angles of eigen rays are larger than zero,so that the advantage of coating E2 is more obvious.

    Fig.6 The anechoic effects of the two coatings when cylinder setting at longitudinal direction

    In this paper,we focus on the difference of the anechoic effect of coatings between in free space and in shallow-water waveguide.After covering a coating,the absolute value of ETS of cylinder still differs greatly between set at cross direction and longitudinal direction,just like the result in reference[8],not repeat any more.

    3.2 The echo characteristics of the Benchmark submarine in shallow-water waveguide

    The Benchmark submarine model is the standard model for checking echo calculation in some countries of NATO,as shown in Fig.7.The size of the model is listed in reference[12].Fig.3 gives the position of the transducer and model and their range is 1000 m.

    The echo levels EL of the Benchmark submarine with various coating as a function of azimuth angle are shown in Fig.8.No matter whether the surface is rigid or covered with anechoic coating,the EL in shallow-water waveguide with N.G-profile is always smaller than that in W.P.G-profile.This is because that N.G-profile makes the rays downward and enlarge the transmission loss.The EL in waveguide with N.G-profile is 9~10 dB lower than the EL in W.P.G-profile for rigid submarine,but for submarine with coating,the difference of EL is only about 4~5 dB,especially near the abeam direction.

    Fig.7 Benchmark submarine model

    Fig.8 The EL of Benchmark submarine with various surface in shallow-water waveguide with N.G-profile or W.P.G-profile(5 kHz)

    Next,we analyze the ETS of the Benchmark submarine with various coatings in shallowwater environment with N.G-profile or W.P.G-profile.When the target is rigid,the ETS changes a little in shallow-water waveguide,which is consistent with the result in reference[8].However when the target is covered with anechoic coating,the ETS in waveguide with N.G-profile is larger than that in waveguide with W.P.G-profile.According to the definition EL=ETS-2TL when the source level is 0 dB.The ETS is about 5 dB higher in waveguide with N.G-profile than with W.P.G-profile,as illustrated in Fig.9.At the same time the 2TL is 9 dB larger because the 2TL is 112.8 dB in waveguide with N.G-profile and 103.8 dB in waveguide with W.P.G-profile if the frequency is 5 kHz and the range is 1000 m.As a result the EL in waveguide with N.G-profile is 4dB lower than the EL in waveguide with W.P.G-profile.

    Fig.9 The ETS of Benchmark submarine with various surfaces in shallow-water waveguide with N.G-profile or W.P.G-profile(5 kHz)

    3.3 The echo reduction of the Benchmark submarine caused by coatings in shallow-water waveguide

    In this section,we will analyze the difference of echo reduction caused by coatings in various environments.The TS of the Benchmark submarine model with rigid and coated surface as a function of azimuth angle in free space is shown in Fig.10.The ETS of the Benchmark submarine with rigid and coated surface as function of azimuth angle in shallow-water waveguide with the two velocity profiles is illustrated in Fig.11 and Fig.12,respectively.

    Fig.10 The TS of Benchmark submarine with rigid and coated surface in free space

    It can be seen from Fig.10 to Fig.12 that the coatings E1 and E2 cause the echo reduction obviously both in free space and shallow-water waveguide.Since the echo of the subma-rine as a function of azimuth angle is complex,to analyze the anechoic effect directly is hard.To make a further quantitative analysis,we calculate the TS reduction in free space and the ETS reduction in waveguide with W.P.G-profile or N.G-profile,as shown in Fig.13.

    Fig.11 The ETS of Benchmark submarine with rigid and coated surface in shallow-water waveguide with W.P.G-profile

    Fig.12 The ETS of Benchmark submarine with rigid and coated surface in shallow-water waveguide with N.G-profile

    Fig.13 The TS reduction and the ETS reduction of Benchmark submarine with anechoic coating E1 or E2

    Fig.13(a)illustrates the TS reduction of the submarine with anechoic coatings E1 and E2 in free space and Figs.13(b)and 13(c)show the ETS reduction of submarine with anechoic coatings E1 and E2 in waveguide with W.P.G-profile and N.G-profile,respectively.We find out that:①The anechoic effect of the coating E2 is a little better than that of the coating E1 in free space,but the difference is small.② The anechoic effect of the coating E1 is about 2~3 dB worse than that of the coating E2 in shallow-water waveguide.The results are due to the difference of echo formation mechanisms in free space and in waveguide:

    (1)In free space,the echo is caused by the highlight on the target surface,because there is only the direct incident-reflection ray.In the examples of this paper,the incident angle of the direct ray is about 1.7°.At this angle the reflection coefficient of the coatings E1 and E2 is nearly equal,this makes the difference of the anechoic effect of coatings E1 and E2 very small.

    (2)In shallow-water waveguide,the echo comes from not only the direct incident-reflection ray but also the combination of various incident and reflection rays with larger angles caused by multi-path effect.For the later the reflection coefficient of large incident angle played a greater role.The reflection coefficient of the coating E1 increases rapidly in small angle,but the reflection coefficient of coating E2 changes very little in the angle range of 0°~30°approximately.That is why the anechoic effect of the coating E1 is not as good as the one of the coating E2 in shallow-water waveguide.

    4 Conclusions

    The Benchmark submarine model and two anechoic coatings with different reflection coefficients are used in this paper.We calculate the echo characteristics of rigid and coated target in various shallow water environments,especially for two typical velocity profiles in winter and summer.The influence of velocity profiles on the anechoic effect of coatings and also the difference of the anechoic effect in free space and in shallow-water waveguide are analyzed.The following results can be obtained:

    (1)The difference of reflection coefficient between the coatings E1 and E2 affects little on the TS reduction in free space.It is because that there is only direct mirror reflection ray in free space but no incident-reflection rays with larger angle.At small angles the difference of the reflection coefficient between the coatings E1 and E2 is very little.

    (2)The difference of reflection coefficient between the coatings E1 and E2 makes larger effect on the ETS reduction in shallow-water waveguide.The anechoic effect of the coating E2 is 2~3dB larger than the coating E1.It is because that the echo in waveguide contains contributions from various ray combinations with larger incident-reflection angles and the reflection coefficient of the coating E2 is smaller at these larger incident angles.In order to obtain a larger anechoic effect in shallow-water environment,the reflection coefficient of the coating should be maintained in small value not only at zero incident angle but also in a range of angle from 0 degree to 30 degrees approximately.

    (3)The anechoic effect of a coating in free space can not entirely reflect the real effect in shallow-water waveguide.We should pay a special attention to this result in both experimental study and theoretical evaluation.

    [1]Ingenito F.Scattering from an object in a stratified medium[J].J Acoust.Soc.Am.,1987,82(6):2051-2059.

    [2]Makris N C,Ingenito F,Kuperman W A.Detection of a submerged object insonified by surface noise in an ocean waveguide[J].J Acoust.Soc.Am.,1994,96(3):1703-1724.

    [3]Makris N C.A spectral approach to 3-D object scattering in layered media applied to scattering from submerged spheres[J].J Acoust.Soc.Am.,1998,104(4):2105-2113.Errata,J Acoust.Soc.Am.,1999,106(1):518.

    [4]Makris N C,Latilal P.A unified model for reverberation and submerged object scattering in a stratified ocean waveguide[J].J Acoust.Soc.Am.,2001,109(3):909-941.

    [5]Purniama Ratilal,Nicholas C Makris.Extinction theorem for object scattering in a stratified Medium[J].J Acoust.Soc.Am.,2001,110(6):2924-2945.

    [6]Latilal P,Lai Yisan,Makris N C.Validity of the sonar equation and Babinet’s principle for scattering in a stratified medium[J].J Acoust.Soc.Am.,2002,112(5):1797-1816.

    [7]Chen Yan,Tang Weilin,Fan Jun.The geometrical acoustic method for calculating the echo of targets submerged in a shallow water waveguide[J].ACTA ACOUSTICA,2010,35(3):335-342.

    [8]Chen Yan,Tang Weilin,Fan Jun,Fan Wei.Echo calculation of targets submerged in a shallow water waveguide based on ray-tracing technique[J].ACTA ACOUSTICA,in publish.

    [9]Tang Weilin.Calculation of acoustic scattering of a nonrigid surface using physical acoustic method[J].ACTA ACOUSTICA,1993,18(1):45-53.

    [10]Michael Jones R,Riley J P,Georges T M.HARPO-A versatile three-dimensional hamiltonian ray-tracing program for acoustic waves in an ocean with irregular bottom[R].U.S.DEPARTMENT OF COMMERCE,National Oceanic and Atmospheric Administration,Environmental Research Laboratories,1986.

    [11]Fan Jun,Tang Weilin.The planar element method for computing target strength(TS)of sonar[C].Acoustical Society of China Youth Conference,[CYCA’99],1999.

    [12]Nell CW,Gilroy L E.An improved BASIS model for the BeTSSi submarine[R].DRDC AtlanticTR 2003-199,2003.

    一边摸一边做爽爽视频免费| 午夜视频精品福利| 村上凉子中文字幕在线| 黑丝袜美女国产一区| 亚洲精品久久国产高清桃花| 最新在线观看一区二区三区| 亚洲国产日韩欧美精品在线观看 | 亚洲电影在线观看av| 精华霜和精华液先用哪个| 嫩草影院精品99| 亚洲精品久久国产高清桃花| 国产精品一区二区免费欧美| 777久久人妻少妇嫩草av网站| 婷婷精品国产亚洲av在线| 久久狼人影院| 国产三级黄色录像| 久久亚洲真实| 麻豆成人午夜福利视频| 校园春色视频在线观看| 脱女人内裤的视频| 成人永久免费在线观看视频| 国产精品亚洲美女久久久| 在线十欧美十亚洲十日本专区| 老司机在亚洲福利影院| 日本黄色视频三级网站网址| 人人妻人人澡欧美一区二区| 国产不卡一卡二| 美女高潮到喷水免费观看| 国产精品自产拍在线观看55亚洲| 极品教师在线免费播放| 狠狠狠狠99中文字幕| 精品无人区乱码1区二区| 啦啦啦韩国在线观看视频| 精品国产超薄肉色丝袜足j| 久久天堂一区二区三区四区| 这个男人来自地球电影免费观看| 国产成人av激情在线播放| 国产亚洲精品久久久久久毛片| 国产主播在线观看一区二区| 母亲3免费完整高清在线观看| 十分钟在线观看高清视频www| 黄色丝袜av网址大全| 亚洲av中文字字幕乱码综合 | 欧美性猛交╳xxx乱大交人| 男人操女人黄网站| 人人妻人人看人人澡| 91字幕亚洲| 亚洲中文av在线| 亚洲人成电影免费在线| 黑人操中国人逼视频| 男人舔奶头视频| 色综合站精品国产| 首页视频小说图片口味搜索| 国产亚洲精品一区二区www| 精品人妻1区二区| 亚洲av电影在线进入| 午夜激情福利司机影院| 亚洲专区字幕在线| 亚洲美女黄片视频| 亚洲国产精品成人综合色| 成年女人毛片免费观看观看9| 国产免费男女视频| 一区二区日韩欧美中文字幕| 高清毛片免费观看视频网站| 日日爽夜夜爽网站| 男女午夜视频在线观看| 亚洲国产精品成人综合色| 午夜久久久在线观看| 欧美日韩精品网址| 久久久久久久久久黄片| 青草久久国产| 高潮久久久久久久久久久不卡| 色哟哟哟哟哟哟| 男女床上黄色一级片免费看| 免费在线观看亚洲国产| 岛国在线观看网站| 精品乱码久久久久久99久播| 午夜精品久久久久久毛片777| 亚洲 欧美一区二区三区| 国产一区二区三区视频了| 午夜免费观看网址| 麻豆国产av国片精品| 日韩大尺度精品在线看网址| 日韩欧美国产一区二区入口| 日本a在线网址| 狂野欧美激情性xxxx| 天天一区二区日本电影三级| 国产一卡二卡三卡精品| 91国产中文字幕| 精品一区二区三区av网在线观看| 老汉色∧v一级毛片| 久久久水蜜桃国产精品网| 精品熟女少妇八av免费久了| 无限看片的www在线观看| 搡老岳熟女国产| 亚洲国产精品成人综合色| 精品乱码久久久久久99久播| 午夜激情福利司机影院| 激情在线观看视频在线高清| 在线免费观看的www视频| 91国产中文字幕| 窝窝影院91人妻| 男人舔奶头视频| 免费观看人在逋| 午夜福利成人在线免费观看| 久久久久久免费高清国产稀缺| 波多野结衣高清作品| 别揉我奶头~嗯~啊~动态视频| 国产精品野战在线观看| 人人妻人人澡欧美一区二区| 久久精品成人免费网站| 怎么达到女性高潮| 亚洲欧美日韩无卡精品| 日韩大码丰满熟妇| 精品久久久久久久人妻蜜臀av| 999精品在线视频| 成人亚洲精品一区在线观看| 久久草成人影院| 精品福利观看| 在线天堂中文资源库| 久久这里只有精品19| 欧美乱码精品一区二区三区| 两个人免费观看高清视频| 桃红色精品国产亚洲av| 精品高清国产在线一区| 午夜影院日韩av| 国产成人精品无人区| 久久伊人香网站| 亚洲 欧美 日韩 在线 免费| 午夜福利免费观看在线| www.熟女人妻精品国产| 国产精品综合久久久久久久免费| 国产成人精品无人区| 国产aⅴ精品一区二区三区波| 悠悠久久av| 嫩草影院精品99| 成人亚洲精品av一区二区| 精品少妇一区二区三区视频日本电影| 国产一区二区在线av高清观看| 免费在线观看成人毛片| 麻豆成人av在线观看| 亚洲成av片中文字幕在线观看| 国产成人av教育| 欧美绝顶高潮抽搐喷水| 亚洲,欧美精品.| 黄色片一级片一级黄色片| 搞女人的毛片| 天堂√8在线中文| 国产亚洲欧美98| 亚洲色图av天堂| 波多野结衣高清作品| 日韩一卡2卡3卡4卡2021年| 免费在线观看视频国产中文字幕亚洲| 国产亚洲精品第一综合不卡| 岛国视频午夜一区免费看| 一本精品99久久精品77| 又紧又爽又黄一区二区| 91成人精品电影| 国产不卡一卡二| 99国产精品一区二区蜜桃av| 日本三级黄在线观看| 正在播放国产对白刺激| 亚洲专区字幕在线| 久久久久亚洲av毛片大全| 亚洲无线在线观看| 欧美zozozo另类| 国产av又大| 国内精品久久久久久久电影| 亚洲国产欧美网| 国产精品综合久久久久久久免费| 国产精品一区二区三区四区久久 | 一夜夜www| 亚洲色图 男人天堂 中文字幕| 亚洲中文日韩欧美视频| 一a级毛片在线观看| 国产精华一区二区三区| 国产伦在线观看视频一区| 久久婷婷人人爽人人干人人爱| 别揉我奶头~嗯~啊~动态视频| 国产免费男女视频| 免费高清视频大片| 国产97色在线日韩免费| 国产精品亚洲一级av第二区| 少妇裸体淫交视频免费看高清 | 欧美日韩乱码在线| 非洲黑人性xxxx精品又粗又长| 高清在线国产一区| av免费在线观看网站| 一区二区三区高清视频在线| 国产av一区二区精品久久| 成人午夜高清在线视频 | 午夜影院日韩av| 色综合站精品国产| 国产成人av教育| 日韩国内少妇激情av| 一区福利在线观看| 亚洲七黄色美女视频| 黄色片一级片一级黄色片| 一本精品99久久精品77| 日韩国内少妇激情av| 日本五十路高清| 国产精品九九99| 1024手机看黄色片| 成年女人毛片免费观看观看9| 在线观看舔阴道视频| 久久久久亚洲av毛片大全| a级毛片a级免费在线| 久久这里只有精品19| av电影中文网址| 亚洲国产精品sss在线观看| 免费女性裸体啪啪无遮挡网站| 首页视频小说图片口味搜索| 欧美绝顶高潮抽搐喷水| 亚洲国产精品999在线| 两个人视频免费观看高清| 99久久综合精品五月天人人| 又黄又爽又免费观看的视频| 99re在线观看精品视频| 国产精品久久久久久人妻精品电影| 熟女少妇亚洲综合色aaa.| 岛国在线观看网站| 亚洲成av片中文字幕在线观看| 亚洲成国产人片在线观看| 一区二区三区激情视频| 叶爱在线成人免费视频播放| 国产99久久九九免费精品| 在线观看一区二区三区| 亚洲成a人片在线一区二区| 亚洲国产中文字幕在线视频| 亚洲人成网站在线播放欧美日韩| 欧美激情高清一区二区三区| 看片在线看免费视频| 一区福利在线观看| 性欧美人与动物交配| 日本黄色视频三级网站网址| 一级片免费观看大全| 久久国产亚洲av麻豆专区| 国产精品久久久人人做人人爽| 老司机午夜十八禁免费视频| 欧美一区二区精品小视频在线| 欧美日韩乱码在线| 国产成人精品久久二区二区免费| 亚洲av美国av| 欧美中文综合在线视频| 人人妻人人澡欧美一区二区| 操出白浆在线播放| 国产99白浆流出| 久久久久久久精品吃奶| 久久久久免费精品人妻一区二区 | 国产精品亚洲av一区麻豆| 日本黄色视频三级网站网址| 亚洲精华国产精华精| 精品卡一卡二卡四卡免费| 国产一卡二卡三卡精品| 一卡2卡三卡四卡精品乱码亚洲| 一区二区三区高清视频在线| 久久久久国内视频| 不卡av一区二区三区| 欧美激情极品国产一区二区三区| 亚洲在线自拍视频| 久久天堂一区二区三区四区| 欧美中文综合在线视频| 巨乳人妻的诱惑在线观看| 欧美日韩福利视频一区二区| 777久久人妻少妇嫩草av网站| 亚洲性夜色夜夜综合| 91大片在线观看| 国产一区在线观看成人免费| 最新在线观看一区二区三区| 男人舔奶头视频| 两性夫妻黄色片| 国语自产精品视频在线第100页| 久久中文字幕一级| 精品卡一卡二卡四卡免费| 国产精品二区激情视频| 国产久久久一区二区三区| 久久人人精品亚洲av| 精品熟女少妇八av免费久了| 久久婷婷成人综合色麻豆| 久久精品91无色码中文字幕| 国产黄片美女视频| 免费人成视频x8x8入口观看| 欧美绝顶高潮抽搐喷水| av在线天堂中文字幕| 人妻丰满熟妇av一区二区三区| 一区福利在线观看| 欧美另类亚洲清纯唯美| 成人三级黄色视频| 国产区一区二久久| 婷婷精品国产亚洲av在线| 桃红色精品国产亚洲av| 精品高清国产在线一区| 国产不卡一卡二| 国产亚洲精品久久久久5区| 少妇被粗大的猛进出69影院| 成年人黄色毛片网站| 婷婷丁香在线五月| 日韩免费av在线播放| 日本免费a在线| 国产精品亚洲av一区麻豆| 精品不卡国产一区二区三区| 日韩三级视频一区二区三区| 两人在一起打扑克的视频| 少妇熟女aⅴ在线视频| 亚洲男人天堂网一区| 久久久久国产一级毛片高清牌| 日韩免费av在线播放| 久久久久久久精品吃奶| 欧美最黄视频在线播放免费| 国产一区在线观看成人免费| 欧美日韩一级在线毛片| 欧美乱色亚洲激情| 色综合亚洲欧美另类图片| 亚洲色图 男人天堂 中文字幕| 亚洲色图av天堂| 日本a在线网址| 国产熟女xx| videosex国产| 一区二区三区国产精品乱码| 久久久久久国产a免费观看| 夜夜躁狠狠躁天天躁| 国产91精品成人一区二区三区| 国产成人啪精品午夜网站| 国产在线观看jvid| 可以免费在线观看a视频的电影网站| av天堂在线播放| 可以在线观看的亚洲视频| 欧美亚洲日本最大视频资源| 精品欧美国产一区二区三| 麻豆成人av在线观看| a级毛片在线看网站| 757午夜福利合集在线观看| 黄频高清免费视频| 岛国在线观看网站| 99久久无色码亚洲精品果冻| 老司机在亚洲福利影院| 久久久久久久精品吃奶| 国产黄片美女视频| 一边摸一边抽搐一进一小说| 久久香蕉精品热| 1024视频免费在线观看| 黑人巨大精品欧美一区二区mp4| 亚洲黑人精品在线| 少妇 在线观看| cao死你这个sao货| 亚洲专区国产一区二区| 高清毛片免费观看视频网站| 精品国产亚洲在线| 国产熟女午夜一区二区三区| 午夜福利高清视频| 女警被强在线播放| 日韩欧美国产一区二区入口| 成年免费大片在线观看| 国产在线观看jvid| 亚洲国产中文字幕在线视频| 免费av毛片视频| 色尼玛亚洲综合影院| 亚洲精品久久成人aⅴ小说| 美女扒开内裤让男人捅视频| 久久99热这里只有精品18| 精品久久久久久成人av| 国产免费av片在线观看野外av| 久久伊人香网站| 久久99热这里只有精品18| 丰满人妻熟妇乱又伦精品不卡| 国产av一区二区精品久久| 精品午夜福利视频在线观看一区| 欧美日韩黄片免| 亚洲精品粉嫩美女一区| 欧美又色又爽又黄视频| av在线天堂中文字幕| 在线av久久热| 国产爱豆传媒在线观看 | 国产成人av教育| 97碰自拍视频| 亚洲成a人片在线一区二区| 国产三级在线视频| 久久久久久亚洲精品国产蜜桃av| 久久国产乱子伦精品免费另类| 国产区一区二久久| 国产一区在线观看成人免费| 欧美国产日韩亚洲一区| 久久中文看片网| 久久国产精品人妻蜜桃| 免费高清视频大片| 色播在线永久视频| 制服丝袜大香蕉在线| 欧美绝顶高潮抽搐喷水| 老鸭窝网址在线观看| 亚洲精品粉嫩美女一区| 一本一本综合久久| 伊人久久大香线蕉亚洲五| 三级毛片av免费| 久久久久亚洲av毛片大全| 两人在一起打扑克的视频| 亚洲欧美日韩无卡精品| 美女扒开内裤让男人捅视频| 18禁国产床啪视频网站| 国产精品电影一区二区三区| 免费高清在线观看日韩| 欧美性长视频在线观看| 1024视频免费在线观看| 国产亚洲欧美精品永久| 国产亚洲精品第一综合不卡| 国产成人欧美在线观看| 婷婷六月久久综合丁香| 精品一区二区三区av网在线观看| 琪琪午夜伦伦电影理论片6080| 人人妻人人看人人澡| 国内毛片毛片毛片毛片毛片| 女性被躁到高潮视频| 色播亚洲综合网| 精品少妇一区二区三区视频日本电影| 欧美色视频一区免费| 亚洲全国av大片| 桃色一区二区三区在线观看| 精品国产一区二区三区四区第35| 精品久久久久久久久久久久久 | 老司机午夜十八禁免费视频| 精华霜和精华液先用哪个| 一卡2卡三卡四卡精品乱码亚洲| 亚洲精品久久成人aⅴ小说| 日韩欧美国产在线观看| 哪里可以看免费的av片| 亚洲av成人一区二区三| а√天堂www在线а√下载| 日韩国内少妇激情av| 老鸭窝网址在线观看| 亚洲成人免费电影在线观看| 1024手机看黄色片| 女人高潮潮喷娇喘18禁视频| 三级毛片av免费| 精品久久久久久,| 色av中文字幕| 亚洲男人天堂网一区| 丁香六月欧美| 搡老妇女老女人老熟妇| 又紧又爽又黄一区二区| 亚洲人成电影免费在线| 在线视频色国产色| 国产久久久一区二区三区| 久久天堂一区二区三区四区| 99热这里只有精品一区 | 精品少妇一区二区三区视频日本电影| 男人操女人黄网站| 午夜免费成人在线视频| 欧美另类亚洲清纯唯美| 美女高潮喷水抽搐中文字幕| 精品卡一卡二卡四卡免费| 亚洲欧美一区二区三区黑人| 国语自产精品视频在线第100页| 一区二区三区精品91| 18禁观看日本| 少妇的丰满在线观看| 日本 av在线| 黑人欧美特级aaaaaa片| 色综合站精品国产| 国产精品国产高清国产av| 十八禁人妻一区二区| 日本免费一区二区三区高清不卡| 成年女人毛片免费观看观看9| 黄色女人牲交| 91老司机精品| 欧美性猛交╳xxx乱大交人| 又黄又爽又免费观看的视频| 亚洲成av人片免费观看| 女人高潮潮喷娇喘18禁视频| 久久伊人香网站| АⅤ资源中文在线天堂| 国产成+人综合+亚洲专区| 欧美黑人巨大hd| 老汉色∧v一级毛片| xxxwww97欧美| 一区二区日韩欧美中文字幕| 日日夜夜操网爽| 免费女性裸体啪啪无遮挡网站| 少妇 在线观看| 欧美日韩一级在线毛片| www日本黄色视频网| 精品第一国产精品| 听说在线观看完整版免费高清| 欧美精品啪啪一区二区三区| 久久九九热精品免费| 成人国产一区最新在线观看| 99热这里只有精品一区 | 天天躁夜夜躁狠狠躁躁| 成人三级做爰电影| 十八禁网站免费在线| av电影中文网址| 免费看十八禁软件| 好男人电影高清在线观看| 精品日产1卡2卡| 免费电影在线观看免费观看| 人人妻人人看人人澡| 免费搜索国产男女视频| 99久久久亚洲精品蜜臀av| 久久精品国产99精品国产亚洲性色| 日韩欧美一区二区三区在线观看| 丁香六月欧美| 日韩av在线大香蕉| 亚洲中文字幕一区二区三区有码在线看 | 国产av在哪里看| 午夜久久久在线观看| 久久久水蜜桃国产精品网| 超碰成人久久| 1024手机看黄色片| 国产又爽黄色视频| 免费观看精品视频网站| 国产又爽黄色视频| 在线av久久热| 一级黄色大片毛片| 一本一本综合久久| 每晚都被弄得嗷嗷叫到高潮| 在线观看免费午夜福利视频| 波多野结衣av一区二区av| 欧美黄色淫秽网站| 久久精品国产清高在天天线| 亚洲在线自拍视频| videosex国产| 中文字幕精品免费在线观看视频| 18禁裸乳无遮挡免费网站照片 | 精品久久蜜臀av无| 男女做爰动态图高潮gif福利片| 午夜久久久久精精品| 国产成年人精品一区二区| 久久久久久久久中文| 男女之事视频高清在线观看| 黄色成人免费大全| 欧美日韩一级在线毛片| 美女午夜性视频免费| 淫妇啪啪啪对白视频| 叶爱在线成人免费视频播放| 一级a爱视频在线免费观看| 欧美日韩福利视频一区二区| 搡老岳熟女国产| 久久久久久久久免费视频了| 免费看美女性在线毛片视频| 精品少妇一区二区三区视频日本电影| 一进一出抽搐动态| 国产真人三级小视频在线观看| 精品久久久久久久久久免费视频| a级毛片a级免费在线| 欧美 亚洲 国产 日韩一| 欧美成人一区二区免费高清观看 | 一级黄色大片毛片| 亚洲欧美激情综合另类| 久久精品国产亚洲av香蕉五月| 日日夜夜操网爽| 亚洲一码二码三码区别大吗| 久久久久久免费高清国产稀缺| 国产成人一区二区三区免费视频网站| 国产主播在线观看一区二区| 国产1区2区3区精品| 两个人视频免费观看高清| 国产一级毛片七仙女欲春2 | 香蕉av资源在线| 午夜日韩欧美国产| 在线观看午夜福利视频| 午夜免费成人在线视频| 中文字幕最新亚洲高清| 国产精品精品国产色婷婷| 国产亚洲av高清不卡| 欧美日韩福利视频一区二区| 日本 av在线| 99国产精品99久久久久| 亚洲五月天丁香| 此物有八面人人有两片| 在线播放国产精品三级| 一个人免费在线观看的高清视频| 精品福利观看| 久久久国产成人免费| 一a级毛片在线观看| 午夜久久久在线观看| 国产亚洲欧美在线一区二区| 777久久人妻少妇嫩草av网站| 男人操女人黄网站| 不卡一级毛片| 精品午夜福利视频在线观看一区| 国产免费av片在线观看野外av| 欧美中文综合在线视频| 久久久久久国产a免费观看| 亚洲专区国产一区二区| 国内精品久久久久久久电影| 亚洲激情在线av| 香蕉久久夜色| 精品免费久久久久久久清纯| 亚洲成人久久性| 亚洲男人的天堂狠狠| 欧美性猛交黑人性爽| 亚洲专区国产一区二区| 国产亚洲欧美98| 99精品欧美一区二区三区四区| 99国产精品一区二区三区| 欧美人与性动交α欧美精品济南到| 久久狼人影院| 日本熟妇午夜| 欧美色视频一区免费| 精品国产一区二区三区四区第35| 免费一级毛片在线播放高清视频| 中亚洲国语对白在线视频| 国产av在哪里看| 黑人欧美特级aaaaaa片| 亚洲欧美一区二区三区黑人| 少妇粗大呻吟视频| 国产高清videossex| 香蕉丝袜av| 悠悠久久av| 久久伊人香网站| 99国产综合亚洲精品| 欧美亚洲日本最大视频资源| 亚洲精品国产一区二区精华液| 国产精品爽爽va在线观看网站 | 亚洲一码二码三码区别大吗| 此物有八面人人有两片| 国产在线精品亚洲第一网站|