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

    Numerical and Experimental Studies of the Acoustic Scattering from an Externally Ring-Stiffened Cylindrical Shell

    2019-07-08 03:54:44ZHOUFulinFANJunWANGBinTANGWeilin
    船舶力學(xué) 2019年6期

    ZHOU Fu-lin,FAN Jun,WANG Bin,TANG Wei-lin

    (1.Collaborative Innovation Center for Advanced Ship and Deep-Sea Exploration,Shanghai Jiao Tong University,Shanghai 200240,China;2.State Key Laboratory of Ocean Engineering,Shanghai Jiao Tong University,Shanghai 200240,China)

    Abstract:Studying the interaction of sound with double shells immersed in water is essential for improving existing underwater target detection and classification algorithms.To identify the key physics,a hybrid 2-D/3-D finite element modeling technique was employed to calculate and analyze the acoustic scattering from an externally ring-stiffened cylindrical shell.Approximate analytical expressions were explicitly derived to identify Bragg and Bloch-Floquet wave scattering in frequency-angle spectra.Monostatic scattering measurements were performed for validation.In contrast to the case of a plate,there are three kinds of elastic waves in a thin cylindrical shell:the flexural wave a0,compressional wave s0and shear wave T0.The flexural wave makes an important contribution to Bloch-Floquet wave scattering.The rings indirectly influence the exterior acoustic scattered field through the dynamic response of the shell,and serve as direct acoustic scatterers that increase the target cross section,especially at oblique incidence.For the external rings at an incident angle of 45°,the geometrical scattering is enhanced due to their influence as corner reflectors.Compared with the analytical method and physical acoustic method,the hybrid 2-D/3-D FE method considers the elasticity of structure,the occlusion and multiple scattering between the rib and the shell,and the interaction between the rib and the water.It has the advantages of high computing efficiency and accuracy,which is suitable for the structural-acoustic coupling analysis of complex stiffened structure in low and middle frequency.

    Key words:acoustic scattering;cylindrical shell;external ring;finite element method

    0 Introduction

    Elastic cylindrical shells often occur as part of double-shell constructions and have been widely used in marine and aerospace research.Double shells are structurally connected to each other by common frames or rings which significantly affect the sound scattered by or radiated from the shells.Based on thin shell theory,Skelton(1992)[1-3]examined the sound scattering from two concentric cylindrical shells,in which the internal rings indirectly influence the exterior acoustic scattered field through the dynamic response of the shell.In target detection,a-coustic waves are easily transmitted through the outer shell to the interior,especially at low frequencies,and the waves are then scattered by the inner shell and the rings in the water area between double shells.In contrast to the inner shell,the outer shell can occasionally be neglected;thus,the externally stiffened shell becomes the main scatterer.Additional factors,including the coupling of water and rings,occlusion and multiple scatterings,remain to be considered[4].In this paper,the acoustic scattering from an externally ring-stiffened cylindrical shell is investigated as a simplification of double shells.

    Numerical methods are generally used to calculate the structural-acoustic response of complex structures,such as the planar element method(PEM)based on the Kirchhoff approximation[5-7]and the finite element method(FEM)[8].Notably,the Kirchhoff approximation is applicable only for geometrical scattering at high frequencies[7].For the FEM,an efficient model with perfectly matched layers(PMLs)was proposed for calculating the radiation and scattering from axisymmetric objects[8].It is called the hybrid 2-D/3-D FE model and has been applied in ocean acoustics and scattering from complex objects[9-10].

    In this study,the scattering characteristics of an externally ring-stiffened cylindrical shell are investigated.Chap.1 introduces hybrid 2-D/3-D FE modeling and presents a numerical calculation of scattering from an externally ring-stiffened cylindrical shell and a uniform shell for comparison.For the case of a finite periodically ribbed cylindrical shell,Chap.2 describes approximate analytical expressions to identify the elastic waves in the frequency-angle spectrum.Chap.3 analyzes the scattering characteristics of an externally ring-stiffened cylindrical shell.The results of a tank experiment are compared with the numerical solutions in Chap.4.Finally,conclusions are presented in Chap.5.

    1 Hybrid 2-D/3-D finite element modeling

    Considering a harmonic plane wave incident upon a thin cylindrical shell submerged in water,the incident wave is assumed to be parallel to the xOz plane with an incident angle θi(see Fig.1).The origin O is located at the center of the shell.The interior of the shell is in vacuo,and the ends of the shell are closed by flat endcaps.A set of rings is attached to the outside surface of the cylinder and seamlessly connected to the cylindrical shell.Therefore,the focus of the present research is the acoustic scattering of axisymmetric fluid-loaded structures subject to a nonsymmetric forcing field.The harmonic time dependence in eiωtwill be omitted.

    Fig.1 Schematic diagram of the externally ringstiffened cylindrical shell model

    Considering Zampolli's method[8],the model superposition technique breaks the 3-D wave equations in the fluid and linear elastodynamic equations in the solid into Fourier series in cylindrical coordinates.The 3-D problem is decomposed into 2-D Fourier modal subproblems to solve the target scattering in the near field.A generic incident acoustic field pinccan be expressed in the following form:

    where n represents the circumferential order and φ is the circumferential angle(0≤φ<2π).Consistent with the decomposition of the applied load,the solid displacement (u,v, )w and pressure p are decomposed in the same manner:

    To simulate the free field,a Berenger PML with cylindrical coordinates is implemented in the azimuthal modal decomposition FE scheme[11].Next,the partial differential equations,the PML and the vibro-acoustic coupling boundary conditions are expressed in a weak form by the variational principle[12].To enable a solution via the FE scheme,these expressions are assembled into an algebraic form that includes a stiffness matrix and mass matrix.Finally,the FE model is built up and solved in the vibro-acoustic coupling module of the software COMSOL.

    Next,the scattered pressure in the far field is calculated by using a 3-D propagation model.In accordance with the discrete form of the Helmholtz integral,the far-field pressure can be calculated as follows[10]:

    where rijis the distance between the jth point on the surface S and the far field observation point i,with the surface normal pointing outwards from the surface S;(pj, dpj/dnj)denotes the scattered pressure as well as its derivative on the surface and Gijis the Green's function in the free field.Finally,the target strength(TS)takes the following form:

    2 Approximate theoretical analysis

    Considering an infinite thin cylindrical shell insonified by acoustic plane waves,in line with the consistency condition of the wave trajectory,the vibration of the cylindrical shell should be the traveling wave in the form of eiKzin the z-direction and symmetric in the circumferential direction for the circumferential angle φ.Thus,a general solution for the shell vibration can be expressed as follows:

    where kφand K denote the circumferential and axial wavenumbers,respectively.Substituting Eq.(5)into the Donnell equations[13]for the displacements of the cylindrical shell yields the following equation containing the coefficient (u~,v~,w~)

    The characteristic matrix is written as follows[13]:

    The characteristic equation for elastic waves is given as follows:

    which is used to solve the phase velocity of elastic waves.Unlike with a thin plate,three types of elastic waves might exist in a thin cylindrical shell because Eq.(8)has three solutions.They correspond to the flexural wave a0,compressional wave s0and shear wave T0.

    At oblique incidence,the elastic waves in a cylindrical shell can be helical,and the corresponding wavenumbers are composited by the axial and circumferential wavenumbers:

    The corresponding phase velocity is cph=ω/kl.The frequency dispersion curves shown in Fig.2(a)can be obtained by solving the characteristic equation.The detailed parameters for the calculation are presented in Chap.4.

    Circumferential resonance occurs for helical waves when the azimuthal wave path equals an integer number of wavelengths,i.e.,kφ=n,where n=1,2,3….Based on the phase velocities of the elastic waves in Fig.2(a),the theoretical trajectories of the three types of circumferential waves are shown in Fig.2(b).

    Fig.2(a)Phase velocities of three types of circumferential waves in an infinite cylindrical shell:flexural wave a0(red dot),compressional wave s0(black dashed-dot)and shear wave T0(blue solid)and(b)Their theoretical trajectories in the frequency-angle spectrum.Each family of curves is identified as follows:a0(n=27~3 2 ),s0(n=1~ 4 )and T0(n=1~ 4 )

    When the infinite cylindrical shell is stiffened periodically by rings,the Bloch-Floquet waves appear in the z-direction.Consistent with the Bloch-Floquet theorem[14],the wavenumber becomes the following:

    where m=-∞,…,0,…∞ is the order of the Bloch-Floquet wave.The flexural waves can be scattered as a meridional leaky ray(in the z-direction)and then contribute to the Bloch-Floquet waves.Like a thin plate,the diagram of the flexural wave propagation is shown in Fig.3(a).The resonance condition of the Bloch-Floquet flexural wave scattering is the following:

    Fig.3(a)Diagram of the Bloch-Floquet flexural wave propagation[(1)forward and(2)backward]in an infinite stiffened cylindrical shell and(b)Its theoretical trajectories in the frequency-angle spectrum

    where the symbols(‘+'and ‘-')represent the forward and backward BFWs,respectively,as shown in Fig.3(a).The Bloch-Floquet trajectories of the flexural waves are shown in Fig.3(b).

    Notably,for the above approximate method,certain deficiencies could not be avoided:(i)the thin shell theory is not applicable to high frequencies;(ii)the detailed sizes and shapes of rings cannot be distinguished;and(iii)the radiation efficiencies of the three elastic waves are distinct but cannot be obtained here because of the neglected fluid loading.

    3 Numerical analysis

    Based on the above FE modeling technique and the approximate theoretical solutions,this chapter provides a numerical analysis of an externally ring-stiffened cylindrical shell.The externally ring-stiffened cylindrical shell used for the computations has similar parameters to those of Ref.[15].The shell has an outer radius a=3.12 cm,total length L=42 cm,thickness h1=1 mm,and ratios of h1/a=0.032 and a/L=0.149.The rings have a height H=1.3 cm,thickness h2=0.9 mm,and rib spacing d=2.1 cm.The elastic parameters are consistent with 304 stainless steel with E=191 GPa,σ=0.279 and ρ=7 570 kg/m3,and the corresponding longitudinal-wave velocity cLand transverse-wave velocity cTare 5 675 m/s and 3 141 m/s,respectively.The external fluid is water,with c0=1 470 m/s and ρ0=1 000 kg/m3.The dimensionless frequency range is 2.7<ka<32.0.

    3.1 Helical waves

    For a comparative analysis,the frequency-angle spectra of the TS for a uniform cylindrical shell and an externally ring-stiffened cylindrical shell are numerically calculated by using Eq.(4),as presented in Fig.4(a)and(b),respectively.The uniform cylindrical shell has no rings in comparison with the externally ring-stiffened cylindrical shell,but they have the same shell parameters.

    Fig.4 Numerical frequency-angle spectra of the target strength for(a)A uniform cylindrical shell and(b)An externally ring-stiffened cylindrical shell

    In the response of the uniform cylindrical shell,bright arcs are prevalent near broadside incidence.Consistent with the circumferential resonance conditions(see Eq.(1)in Ref.[16]),these arcs are identified as the circumferential coupling loci for the helical propagation of the three elastic waves,as indicated in Fig.4(a).The compressional waves and shear waves,whose behaviors at low frequencies are similar to those in a plate,can be launched within the region from a broadside incidence to cutoff angles γL=sin-1(c0/ cL)=15.0°and γT=sin-1(c0/ cT)=27.9°,where cLand cTare the compressional and shear velocities in the shell material,respectively.High backscattering levels arise from the launching of helical leaky waves,which are reflected off the shell truncation and reradiated in the backscattering direction.The coupling loci for the a0and T0waves are more visible than those for the s0wave.Further,the bright arcs for the low-order T0waves are wider and stronger while those for s0waves are covered easily.Besides,for metal shells with external fluid loading,the a0wave is divided into two waves,namely,a0+and a0-.The a0+wave is similar to the a0wave,whereas an additional a0-wave is derived from the fluid loading and is also called the Scholte-Stoneley wave.Above the critical frequency corresponding to the sound speed in ambient fluid,a0+is a flexural and a0-a fluid-borne wave,and vice versa below the frequency[14].Further details on the scattering characteristics of helical waves in a finite cylindrical shell are discussed in Refs.[16-18].

    Compared with the helical wave coupling loci in a uniform cylindrical shell,such as those in Fig.4(a),most helical waves become weak in Fig.4(b)for an externally ring-stiffened cylindrical shell because the helical propagation is blocked by rings.Near broadside incidence(0°-30°),helical waves at low frequencies(ka<15)have more influence on the exterior scattered field than those at high frequencies(ka>15).The bright arcs for the low-order T0waves are still visible.Additionally,many bright arcs appear in other regions.Considering previous studies on ribbed structures(e.g.,Refs.[3,14]),most of the waves are naturally considered to be generated via Bragg and BFW scattering,which will be discussed in the next section.

    3.2 Bragg scattering and Bloch-Floquet wave scattering

    To identify their trajectories,approximate theoretical solutions of Bragg and Bloch-Floquet scattering from an infinite ring-stiffened cylindrical shell are helpful.The same mechanisms exist for a finite cylindrical shell,with the only difference being the energy strongly scattered from the ends of the shell.According to the theoretical analysis in Chap.3,the bright arcs in Fig.4(b)can be identified as shown in Fig.5(a).

    Fig.5 Numerical frequency-angle spectra of the target strength for(a)An externally ring-stiffened cylindrical shell and(b)An internally ring-stiffened cylindrical shell

    First,Bragg scattering arises from the geometrical reflections of an incident wave on the inner and outer edges of the rings.These two parts have different arrival times,and their trajectories on the frequency-angle spectrum are consistent and can be evaluated as follows:

    where k0is the wavenumber of the incident plane wave and d is the rib spacing.A number of marked bright arcs show excellent agreement with the theoretical Bragg trajectories of the first six orders.

    Then,according to the theoretical trajectories in Fig.3(b),the Bloch-Floquet flexural waves with two low orders can be identified,as denoted as BFW ( a0)in Fig.5(a).The trajectories of flexural waves with high orders are more obvious than those with low orders.The same phenomenon occurs in the case of the internally ring-stiffened shell in Fig.5(b).In addition,according to the geometrical relationship of the ring height and ring spacing,sound waves are reflected by external rings and cannot be directly incident on the cylindrical shell when the incident angle is greater than 58°.Although the transmitted and diffracted waves can arrive on the cylindrical shell,the received echoes are weakened slightly.

    According to the theoretical analysis in Chap.2,three types of elastic waves exist in the cylindrical shell.In the frequency-angle spectrum of the ring-stiffened shell,not all of the bright arcs can be identified for various reasons.Many factors make the scattering problem more complex,such as the coupling between the rings and the shell and the blocking or partition effect of rings on the surface helical waves.The rings transmit the coupling forces and moments to the shell and contribute to the dynamic response,which indirectly influences the exterior acoustic field.Moreover,external rings become direct elastic scatterers in water.

    3.3 Effects of external rings

    To identify the contributions of the external rings,an internally ring-stiffened cylindrical shell and an externally ring-stiffened cylindrical shell with a perfectly rigid surface are numerically modeled for comparison.Other parameters,such as the ring spacing,are the same for the both shells.For the quantitative analysis,the target strength as a function of the incident angle for the three models is shown in Fig.6.

    Compared with the results for a uniform shell,the TS of the externally ring-stiffened shell is stronger as a whole and has more peaks at oblique incidence.These peaks are caused by the geometrical and elastic scattering from the rings,including the shielding effects and multiple scatterings between the rings and the shell.According to the frequency-angle spectrum in Fig.4(a),helical waves in the uniform shell can be labeled as HW in Fig.6(a)and(b),which appear in small incident angle range.

    Fig.6 Comparison of the target strength for the backscattering of sound by an externally ringstiffened shell(black solid),an internally ring-stiffened shell(red solid),and a uniform cylindrical shell(blue dot)at(a)100 kHz(ka=13.3)and(b)200 kHz(ka=26.7)

    For the two ring-stiffened shells,bright arcs of Bragg scattering from periodic rings are observed in Fig.5(a)and(b).Their peaks in Fig.6 are also labeled as Bragg orders.As expected,Bragg scattering from the external rings is generally stronger than it is from internal rings.The difference value of the TS between the two cases exceeds 10 dB at several incident angles and frequencies(e.g.,m=2~4 at 200 kHz).There are several reasons for this result.An external ring has two edges that directly reflect sound waves and contribute to strong Bragg scattering.For an internal ring,the incident waves need to be transmitted through the shell twice to encounter the ring,and thus more sound energy has to be consumed.From the aspect of mechanics,internal rings exert only forces and bending moments on the shell,thereby indirectly influencing the exterior scattered field.In addition to the indirect influence,the external rings as direct acoustic scatterers increase the target cross section,especially at oblique incidence.Moreover,a remarkable difference is observed in which the amplitude of the Bragg scattering from external rings increases progressively in the range of 0°-45°.In addition to the geometrical scattering,more peaks appear due to the elasticity of the rings,some of which are identified in Secs.4.1 and 4.2.Generally,the external rings enhance the target strength.

    In addition,the TS of the externally ringstiffened cylindrical shell at 45°is less than those at broadside and end-on but much larger than those of the internally ring-stiffened and uniform cylindrical shells.This is because the external rings at 45°serve as corner reflectors and enhance the geometrical scattering,as shown in Fig.7.Furthermore,by comparing the average TS over the frequency band(20-240 kHz)of the three models,an enhancement near 45°is also observable for the case of the externally ring-stiffened cylindrical shell.

    Fig.7 Corner reflector at an incident angle of 45°

    4 Experiment

    4.1 Experimental program and data processing

    An experiment was conducted in an acoustic water tank with dimensions of 5×5×5 m3.The experimental configuration is shown in Fig.8.A ring-stiffened cylindrical steel shell was manufactured.In consideration of the limited space of the tank and the convenience of model machining,the length of the cylindrical shell is L=21 cm,and the other parameters of the material and geometry are presented in Chap.3.The trajectories of the Bragg scattering and Bloch-Floquet wave scattering are not changed.Considering the broad bandwidth of 20-240 kHz(ka=2.7-32.0),we employed four transducers with bandwidths of 20-40 kHz,40-80 kHz,60-120 kHz,and 120-240 kHz.Tab.1 summarizes their specifications.Each transducer outputs a linear-phase finite impulse response(FIR)signal.A standard Teledyne RESON hydrophone(TC4014)was used,and its sensitivity is-186 dB±3 dB re1V/μPa over the frequency range 15 Hz to 480 kHz.The whole experimental spectra were completed by overlapping and stitching together the four slices of data.

    Fig.8 Experimental configuration

    Tab.1 Specifications of the four transducers used in the tank experiment

    To obtain the full-azimuth scattering characteristics,the model was automatically rotated in the horizontal xOz plane by drivers at a rotational speed of 60°/min.Continuous data acquisition was achieved in a time interval of 0.5 s to ensure adequate resolution of the experimental spectra and avoid measurement errors from the model's waggle to the greatest extent possible.Having obtained the direct waves and reflections from the model through appropriate time windows,we performed an FFT on the time signals for each incidence angle to obtain the experimental frequency-angle spectrum.

    4.2 Experimental analysis

    Based on previous data processing methods,the experimental frequency-angle spectrum of the externally ring-stiffened cylindrical shell is shown in Fig.9(a).For comparison,the simulated result of a model with the same size is shown in Fig.9(b).The scattering characteristics of the simulated and experimental results are generally consistent.Distinct highlights are observed in both figures that are consistent with the theoretical Bragg scattering trajectories of the first six orders in Fig.4.According to the theoretical Bloch-Floquet condition,we identi-fied the experimental spectrum,as labeled in Fig.9.Some highlights in the experimental spectrum are ambiguous because of the relatively weak echoes from rings and the low signal-tonoise ratio.

    Fig.9 Frequency-angle spectra of the externally ring-stiffened cylindrical shell(a)Experiment and(b)Simulation

    5 Conclusions

    A hybrid 2-D/3-D finite element modeling technique was employed to investigate the scattering characteristics of an externally ring-stiffened cylindrical shell.For a contrastive analysis,relevant FE models were also designed,including a nonribbed model and an internally ring-stiffened model.To identify the Bragg and Bloch-Floquet trajectories in the frequencyangle spectra,approximate analytical expressions for a periodically ring-stiffened cylindrical shell were derived.The combination of the hybrid FE model and the approximate solution revealed different aspects of the physics driving the target response.

    The effects of external rings and the structural elasticity on scattering were analyzed by comparison with a uniform elastic FE model and a perfectly rigid FE model.In addition to the indirect influence of the forces and bending moments on the shell similar to internal rings,the external rings in water act as direct acoustic objects that increase the target cross section,especially at oblique incidence.For the external rings at 45°,the geometrical scattering is greatly enhanced via their influence as corner reflectors.

    The validity of the hybrid FE models of the externally ring-stiffened cylindrical shell was confirmed by comparing the simulated frequency-angle spectrum with the monostatic scattering measured data.The Bragg and Bloch-Floquet trajectories were generally consistent,and any discrepancies might arise from the effects of the external rings.Some highlights in the experimental spectrum are ambiguous because of the relatively weak echoes from the rings and a somewhat low signal-noise-ratio during the measurement.

    Acknowledgements

    We thank Dr.Peng Zilong and Dr.Tong Yunzhe for their helpful discussions.We also appreciate the help of Xu Kai,Dr.Zhao Kaiqi and Dr.Liu Bo in performing the tank experiments.

    建设人人有责人人尽责人人享有的| 麻豆成人av在线观看| 99香蕉大伊视频| 99国产精品一区二区蜜桃av | 黄色视频不卡| 99re在线观看精品视频| 亚洲熟女毛片儿| 男男h啪啪无遮挡| 校园春色视频在线观看| 制服诱惑二区| 天天操日日干夜夜撸| 免费女性裸体啪啪无遮挡网站| 亚洲精品国产精品久久久不卡| 久久精品国产清高在天天线| 亚洲国产精品sss在线观看 | 久久久久视频综合| 在线观看午夜福利视频| 男女免费视频国产| 国产午夜精品久久久久久| 高清视频免费观看一区二区| 捣出白浆h1v1| 一本一本久久a久久精品综合妖精| 亚洲全国av大片| 国产乱人伦免费视频| 色综合婷婷激情| 成人永久免费在线观看视频| 9色porny在线观看| 久久久久国产一级毛片高清牌| 在线观看免费午夜福利视频| 亚洲一区二区三区欧美精品| 狠狠狠狠99中文字幕| 中文字幕另类日韩欧美亚洲嫩草| 亚洲精品美女久久久久99蜜臀| 18禁国产床啪视频网站| 宅男免费午夜| 天天影视国产精品| 嫩草影视91久久| 国产亚洲欧美98| 国产亚洲欧美在线一区二区| 国产淫语在线视频| 亚洲少妇的诱惑av| 精品国产超薄肉色丝袜足j| 757午夜福利合集在线观看| 看黄色毛片网站| 国产人伦9x9x在线观看| 两人在一起打扑克的视频| 欧美成狂野欧美在线观看| 亚洲av日韩在线播放| 日本黄色日本黄色录像| 在线观看免费视频网站a站| 伊人久久大香线蕉亚洲五| 欧美成人午夜精品| 搡老熟女国产l中国老女人| 午夜成年电影在线免费观看| 久久久国产成人精品二区 | 国产三级黄色录像| 热99久久久久精品小说推荐| 国产一区有黄有色的免费视频| 国产精品久久久久久精品古装| 国产视频一区二区在线看| 老司机亚洲免费影院| 免费在线观看日本一区| 亚洲欧洲精品一区二区精品久久久| e午夜精品久久久久久久| 精品无人区乱码1区二区| 99riav亚洲国产免费| 亚洲精华国产精华精| 满18在线观看网站| 国产欧美日韩精品亚洲av| 免费少妇av软件| 日韩制服丝袜自拍偷拍| av中文乱码字幕在线| 久热爱精品视频在线9| 一级黄色大片毛片| 天天躁日日躁夜夜躁夜夜| 窝窝影院91人妻| 丰满人妻熟妇乱又伦精品不卡| 国产单亲对白刺激| 操出白浆在线播放| 男人操女人黄网站| 日韩欧美免费精品| 五月开心婷婷网| 午夜福利,免费看| 精品乱码久久久久久99久播| 久久亚洲精品不卡| 丁香欧美五月| 精品人妻在线不人妻| 精品亚洲成a人片在线观看| 99国产极品粉嫩在线观看| 免费看a级黄色片| 久久国产乱子伦精品免费另类| 欧美精品高潮呻吟av久久| av天堂在线播放| 久久精品国产a三级三级三级| 纯流量卡能插随身wifi吗| 午夜视频精品福利| 午夜精品国产一区二区电影| 12—13女人毛片做爰片一| 女同久久另类99精品国产91| 国产成+人综合+亚洲专区| 老司机深夜福利视频在线观看| 欧美成人午夜精品| av线在线观看网站| 日韩制服丝袜自拍偷拍| 国产亚洲av高清不卡| 777久久人妻少妇嫩草av网站| 王馨瑶露胸无遮挡在线观看| 高清在线国产一区| 两个人看的免费小视频| 精品国内亚洲2022精品成人 | 国产精品欧美亚洲77777| 婷婷丁香在线五月| 91麻豆av在线| 亚洲第一欧美日韩一区二区三区| 男女高潮啪啪啪动态图| 50天的宝宝边吃奶边哭怎么回事| 一本大道久久a久久精品| 国产免费男女视频| 91在线观看av| av不卡在线播放| 不卡av一区二区三区| 香蕉久久夜色| av电影中文网址| 男男h啪啪无遮挡| 亚洲中文日韩欧美视频| netflix在线观看网站| 国产在线精品亚洲第一网站| 一区二区三区国产精品乱码| netflix在线观看网站| 在线十欧美十亚洲十日本专区| 岛国在线观看网站| 中文亚洲av片在线观看爽 | 国产成人一区二区三区免费视频网站| 人人妻人人澡人人看| 亚洲精品国产一区二区精华液| 成人三级做爰电影| 欧美一级毛片孕妇| 高清av免费在线| 一级片'在线观看视频| 大片电影免费在线观看免费| 成人三级做爰电影| 午夜久久久在线观看| 亚洲成国产人片在线观看| 亚洲成国产人片在线观看| 成年女人毛片免费观看观看9 | 人人妻人人添人人爽欧美一区卜| 人人妻人人添人人爽欧美一区卜| 看免费av毛片| 免费一级毛片在线播放高清视频 | 亚洲男人天堂网一区| 精品人妻熟女毛片av久久网站| 日韩三级视频一区二区三区| 欧美日韩中文字幕国产精品一区二区三区 | 欧美乱色亚洲激情| 国产亚洲精品一区二区www | 欧美日韩中文字幕国产精品一区二区三区 | 国产亚洲精品久久久久5区| 亚洲一码二码三码区别大吗| 日韩大码丰满熟妇| 中文字幕人妻熟女乱码| 18禁黄网站禁片午夜丰满| 色婷婷av一区二区三区视频| 婷婷丁香在线五月| 久久人人爽av亚洲精品天堂| 欧美精品啪啪一区二区三区| 国产精品亚洲一级av第二区| 久久久久久久国产电影| 亚洲少妇的诱惑av| 涩涩av久久男人的天堂| 在线天堂中文资源库| 久久天躁狠狠躁夜夜2o2o| 亚洲精品乱久久久久久| 9色porny在线观看| 日本wwww免费看| 亚洲成人手机| 精品久久久久久,| 电影成人av| 久久ye,这里只有精品| 国产亚洲av高清不卡| 国产欧美日韩精品亚洲av| 午夜精品久久久久久毛片777| 高潮久久久久久久久久久不卡| 丝袜在线中文字幕| 日韩三级视频一区二区三区| 国产精品秋霞免费鲁丝片| 女人爽到高潮嗷嗷叫在线视频| 日韩精品免费视频一区二区三区| 最新的欧美精品一区二区| 精品国产乱子伦一区二区三区| 国产99白浆流出| 日本黄色日本黄色录像| bbb黄色大片| 99久久99久久久精品蜜桃| 国产一区有黄有色的免费视频| 亚洲国产欧美一区二区综合| 午夜成年电影在线免费观看| 丰满饥渴人妻一区二区三| 99国产精品免费福利视频| 又大又爽又粗| 人人妻人人添人人爽欧美一区卜| 亚洲欧美激情在线| 午夜福利一区二区在线看| 久久精品亚洲av国产电影网| svipshipincom国产片| 色老头精品视频在线观看| 国产欧美日韩一区二区三| 一级片免费观看大全| 精品久久蜜臀av无| av片东京热男人的天堂| 国产欧美日韩一区二区精品| 老司机福利观看| 午夜成年电影在线免费观看| 99精品在免费线老司机午夜| 99国产精品一区二区三区| 美女视频免费永久观看网站| 欧美午夜高清在线| 午夜福利免费观看在线| 精品福利永久在线观看| 男人操女人黄网站| 美女 人体艺术 gogo| 中文字幕色久视频| 免费在线观看日本一区| 亚洲午夜理论影院| 欧美国产精品va在线观看不卡| 一级黄色大片毛片| 婷婷精品国产亚洲av在线 | 日本黄色日本黄色录像| 视频区欧美日本亚洲| 国内毛片毛片毛片毛片毛片| 成人手机av| av片东京热男人的天堂| 人人妻人人添人人爽欧美一区卜| 久久精品国产综合久久久| 欧美日韩黄片免| 国产精品久久久av美女十八| 久久久国产精品麻豆| 欧美大码av| 老司机午夜十八禁免费视频| 午夜久久久在线观看| 免费在线观看视频国产中文字幕亚洲| 身体一侧抽搐| 午夜福利,免费看| 中文欧美无线码| 最新美女视频免费是黄的| 亚洲精品一卡2卡三卡4卡5卡| 国产真人三级小视频在线观看| 一个人免费在线观看的高清视频| 国产极品粉嫩免费观看在线| 日韩欧美在线二视频 | av视频免费观看在线观看| 久久人妻福利社区极品人妻图片| 欧美乱妇无乱码| 国产片内射在线| 热99国产精品久久久久久7| 久久人人97超碰香蕉20202| 国产不卡av网站在线观看| 亚洲欧美一区二区三区黑人| 热99re8久久精品国产| 欧美大码av| 777久久人妻少妇嫩草av网站| 亚洲国产看品久久| 亚洲国产精品一区二区三区在线| 日本wwww免费看| 欧美最黄视频在线播放免费 | 国产亚洲欧美98| 欧美最黄视频在线播放免费 | 夜夜躁狠狠躁天天躁| 亚洲欧美精品综合一区二区三区| 免费黄频网站在线观看国产| 精品人妻熟女毛片av久久网站| 色精品久久人妻99蜜桃| 成人国语在线视频| av国产精品久久久久影院| 亚洲中文字幕日韩| 午夜福利乱码中文字幕| 一夜夜www| 日本vs欧美在线观看视频| 真人做人爱边吃奶动态| 黄色 视频免费看| 国产欧美日韩一区二区三| 久久久久久久久久久久大奶| 自拍欧美九色日韩亚洲蝌蚪91| 国产成+人综合+亚洲专区| 久久精品国产a三级三级三级| 99国产精品免费福利视频| 久久久久国内视频| 夫妻午夜视频| 欧美性长视频在线观看| 国产精品久久久av美女十八| 亚洲一区二区三区欧美精品| 在线播放国产精品三级| av免费在线观看网站| 国产真人三级小视频在线观看| 精品国产乱码久久久久久男人| 很黄的视频免费| 一级毛片高清免费大全| 国产精品久久久av美女十八| 国产精品自产拍在线观看55亚洲 | 又黄又爽又免费观看的视频| 亚洲美女黄片视频| 国产三级黄色录像| 欧美成狂野欧美在线观看| 欧美日韩亚洲国产一区二区在线观看 | 欧美中文综合在线视频| 老司机靠b影院| 美女高潮喷水抽搐中文字幕| 一本大道久久a久久精品| 在线视频色国产色| 亚洲七黄色美女视频| 正在播放国产对白刺激| 亚洲精品久久午夜乱码| 热re99久久国产66热| 女人高潮潮喷娇喘18禁视频| 精品久久久久久久毛片微露脸| 丝袜人妻中文字幕| 91精品三级在线观看| 一边摸一边抽搐一进一小说 | 国产视频一区二区在线看| 亚洲精品久久成人aⅴ小说| 人成视频在线观看免费观看| 亚洲欧美日韩高清在线视频| 国产成人系列免费观看| 这个男人来自地球电影免费观看| 香蕉久久夜色| www.999成人在线观看| 99久久综合精品五月天人人| 欧美亚洲日本最大视频资源| 在线看a的网站| 成人18禁高潮啪啪吃奶动态图| 国产深夜福利视频在线观看| 最近最新免费中文字幕在线| 久久人妻av系列| 男女午夜视频在线观看| 久久久国产精品麻豆| 日本撒尿小便嘘嘘汇集6| 看黄色毛片网站| 亚洲全国av大片| 叶爱在线成人免费视频播放| 精品欧美一区二区三区在线| 老司机在亚洲福利影院| 丰满的人妻完整版| 亚洲五月婷婷丁香| 人人澡人人妻人| 国产精品综合久久久久久久免费 | 99精品久久久久人妻精品| 窝窝影院91人妻| 老熟妇仑乱视频hdxx| 久久精品人人爽人人爽视色| 久久久精品国产亚洲av高清涩受| 亚洲av成人一区二区三| 免费看a级黄色片| 免费观看人在逋| 欧美一级毛片孕妇| a级毛片黄视频| 精品人妻在线不人妻| 这个男人来自地球电影免费观看| 亚洲精品美女久久久久99蜜臀| 制服诱惑二区| 啦啦啦视频在线资源免费观看| 91字幕亚洲| 丝袜美腿诱惑在线| 久久草成人影院| 国产精品一区二区精品视频观看| a级片在线免费高清观看视频| 久久午夜亚洲精品久久| 日日爽夜夜爽网站| 老熟女久久久| 亚洲av熟女| 精品国产乱子伦一区二区三区| 男人舔女人的私密视频| 国产精品一区二区在线观看99| 亚洲精品国产一区二区精华液| 9色porny在线观看| 成人av一区二区三区在线看| 国产男女超爽视频在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲成人免费电影在线观看| 国产一卡二卡三卡精品| 国产熟女午夜一区二区三区| 精品国产乱子伦一区二区三区| 久久精品国产a三级三级三级| 18禁美女被吸乳视频| 黄色丝袜av网址大全| 久久久久久人人人人人| 欧美日韩中文字幕国产精品一区二区三区 | 久久99一区二区三区| 国产高清国产精品国产三级| 超碰97精品在线观看| 欧美日本中文国产一区发布| 69精品国产乱码久久久| 一区二区三区国产精品乱码| 久久精品国产a三级三级三级| 大型av网站在线播放| 国产日韩一区二区三区精品不卡| 国产无遮挡羞羞视频在线观看| 国产成人精品无人区| 亚洲精品久久午夜乱码| 好男人电影高清在线观看| 90打野战视频偷拍视频| 久久久国产成人精品二区 | 成人永久免费在线观看视频| 国产欧美日韩一区二区三| 亚洲人成电影免费在线| 久久精品国产综合久久久| 黄色女人牲交| 国产色视频综合| 香蕉久久夜色| 午夜福利一区二区在线看| 亚洲人成电影免费在线| 亚洲自偷自拍图片 自拍| 精品国产美女av久久久久小说| 成人黄色视频免费在线看| 日韩免费高清中文字幕av| 亚洲国产看品久久| 丝袜人妻中文字幕| 久久 成人 亚洲| 天天躁狠狠躁夜夜躁狠狠躁| 夜夜爽天天搞| 国产精品一区二区免费欧美| 久久久久久久国产电影| 最近最新中文字幕大全免费视频| 9色porny在线观看| 女人高潮潮喷娇喘18禁视频| 午夜福利在线观看吧| 成人精品一区二区免费| 亚洲视频免费观看视频| 国产亚洲欧美98| 久热爱精品视频在线9| 成年人免费黄色播放视频| 久久久国产成人精品二区 | 国产99久久九九免费精品| 极品少妇高潮喷水抽搐| 久久天躁狠狠躁夜夜2o2o| 少妇猛男粗大的猛烈进出视频| 国产成人系列免费观看| 欧美日韩视频精品一区| 美国免费a级毛片| 日本精品一区二区三区蜜桃| 精品欧美一区二区三区在线| 色94色欧美一区二区| 色94色欧美一区二区| 黑人猛操日本美女一级片| 亚洲国产欧美日韩在线播放| 99国产综合亚洲精品| 欧美大码av| av中文乱码字幕在线| 视频在线观看一区二区三区| 美女福利国产在线| av视频免费观看在线观看| 美女扒开内裤让男人捅视频| 国产在线一区二区三区精| 欧美成狂野欧美在线观看| 一区二区三区激情视频| 最近最新免费中文字幕在线| 国产精品98久久久久久宅男小说| 久久99一区二区三区| 欧美精品啪啪一区二区三区| 天堂中文最新版在线下载| 精品卡一卡二卡四卡免费| 国产精品一区二区在线观看99| 国产激情欧美一区二区| 黄色女人牲交| 免费久久久久久久精品成人欧美视频| 国产一区有黄有色的免费视频| 欧美黄色淫秽网站| 国产野战对白在线观看| 两人在一起打扑克的视频| 十八禁高潮呻吟视频| 久久久精品国产亚洲av高清涩受| 精品电影一区二区在线| 欧美在线一区亚洲| 亚洲一区中文字幕在线| 亚洲九九香蕉| 成人18禁高潮啪啪吃奶动态图| 日韩 欧美 亚洲 中文字幕| 一进一出抽搐动态| 悠悠久久av| 午夜精品在线福利| 午夜两性在线视频| 久久精品91无色码中文字幕| 一本综合久久免费| 久久精品国产亚洲av香蕉五月 | 亚洲精品粉嫩美女一区| 成人18禁高潮啪啪吃奶动态图| 国产高清激情床上av| 大码成人一级视频| 久久国产精品男人的天堂亚洲| 久久精品国产99精品国产亚洲性色 | 欧美+亚洲+日韩+国产| 超色免费av| 色老头精品视频在线观看| 日韩欧美免费精品| 精品少妇久久久久久888优播| 日本欧美视频一区| 国产成人一区二区三区免费视频网站| 精品福利观看| 欧美日韩亚洲综合一区二区三区_| 久久国产精品人妻蜜桃| 欧美精品高潮呻吟av久久| 一级作爱视频免费观看| 久久香蕉精品热| 男女免费视频国产| 色播在线永久视频| 男女床上黄色一级片免费看| 免费久久久久久久精品成人欧美视频| 自线自在国产av| 青草久久国产| 欧美精品亚洲一区二区| 日韩视频一区二区在线观看| 久久人人爽av亚洲精品天堂| 国产欧美日韩一区二区三| 久久久国产成人精品二区 | 乱人伦中国视频| 国产av又大| 成人特级黄色片久久久久久久| 亚洲黑人精品在线| 黑人操中国人逼视频| 精品人妻熟女毛片av久久网站| 法律面前人人平等表现在哪些方面| 亚洲精品在线观看二区| 久久午夜亚洲精品久久| 水蜜桃什么品种好| 精品免费久久久久久久清纯 | 美女高潮喷水抽搐中文字幕| 亚洲国产中文字幕在线视频| 日本五十路高清| 在线观看日韩欧美| 一二三四在线观看免费中文在| 久久中文字幕人妻熟女| 国产99久久九九免费精品| 少妇被粗大的猛进出69影院| 中文字幕色久视频| 日韩欧美国产一区二区入口| 如日韩欧美国产精品一区二区三区| 老熟妇乱子伦视频在线观看| 免费观看人在逋| 亚洲国产欧美网| 女性生殖器流出的白浆| av国产精品久久久久影院| 国产一区二区三区综合在线观看| 最近最新中文字幕大全电影3 | 黄色成人免费大全| 精品国产乱码久久久久久男人| 看黄色毛片网站| 大香蕉久久网| 脱女人内裤的视频| 欧美激情 高清一区二区三区| 夫妻午夜视频| 久久狼人影院| 十分钟在线观看高清视频www| 999久久久国产精品视频| 精品熟女少妇八av免费久了| 久久久久久免费高清国产稀缺| 国产精华一区二区三区| 国产亚洲精品久久久久久毛片 | 黑人操中国人逼视频| 欧美激情极品国产一区二区三区| 怎么达到女性高潮| 黄片小视频在线播放| 欧美久久黑人一区二区| 亚洲男人天堂网一区| 国产精品久久电影中文字幕 | 久久久久久久久久久久大奶| 亚洲欧美日韩另类电影网站| 久久性视频一级片| 女性被躁到高潮视频| 女性生殖器流出的白浆| 啦啦啦视频在线资源免费观看| 国产一区二区三区综合在线观看| 我的亚洲天堂| 久久香蕉国产精品| 99国产精品99久久久久| 叶爱在线成人免费视频播放| 咕卡用的链子| 久久人妻福利社区极品人妻图片| 精品久久久精品久久久| 久久中文字幕人妻熟女| 午夜91福利影院| 91成年电影在线观看| 脱女人内裤的视频| 久久久久精品人妻al黑| 精品久久久久久久久久免费视频 | 中文字幕精品免费在线观看视频| 日韩视频一区二区在线观看| 亚洲av片天天在线观看| 国产精品一区二区在线观看99| 男男h啪啪无遮挡| 99国产综合亚洲精品| 好看av亚洲va欧美ⅴa在| 少妇裸体淫交视频免费看高清 | 国产免费av片在线观看野外av| 国产成人啪精品午夜网站| 久久亚洲精品不卡| 国产精品 欧美亚洲| 80岁老熟妇乱子伦牲交| 男女高潮啪啪啪动态图| 久热这里只有精品99| 国产精品 欧美亚洲| 亚洲,欧美精品.| 午夜免费鲁丝| 美女高潮到喷水免费观看| 亚洲av片天天在线观看| 中文字幕另类日韩欧美亚洲嫩草| 亚洲国产精品合色在线| 一级片免费观看大全| 亚洲熟女毛片儿| 69精品国产乱码久久久| www.自偷自拍.com| 夜夜夜夜夜久久久久| 亚洲欧美色中文字幕在线| 久久久国产一区二区| 久久精品国产清高在天天线| 丰满迷人的少妇在线观看| 亚洲精品国产精品久久久不卡| 女人爽到高潮嗷嗷叫在线视频|