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

    Numerical Estimation of Crack Depth in Curved Surfaces of Polymer Bonded Explosives Using Rayleigh Surface Wave

    2020-08-10 06:43:42JIANGChangYANGZhanfengLIWeibinZHANGWeibinTIANYong
    含能材料 2020年8期

    JIANG Chang,YANG Zhan?feng,LI Wei?bin,ZHANG Wei?bin,TIAN Yong

    (1. Institute of Chemical Materials,China Academy of Engineering Physics,Mianyang 621999,China;2. School of Aerospace Engineering,Xiamen University,Xiamen 361005,China)

    Abstract:Cracks in polymer bonded explosives(PBXs)may occur during manufacturing,transportation and storage stages. The surface crack of explosives has an important influence on their mechanical properties and detonation performance. Quantitative detection of the depth of the surface crack is of great significance to the process of surface crack removal. The numerical simula?tion of surface crack depth estimation was investigated by using Rayleigh surface waves which propagate along the cylindrical and spherical surfaces of PBXs. The propagation patterns of surface waves were studied on these two types of surfaces. The reflec?tion and transmission of surface waves at different frequencies in PBXs with crack of varied depth were simulated through the fi?nite element method(FEM).The FEM results indicate that the depth of crack in curved surfaces can be deduced according to the amplitude variations of reflected and transmitted waves.

    Key words:Rayleigh wave;curved surface;polymer bonded explosive;crack depth;simulation

    1 Introduction

    Polymer bonded explosives(PBX),as an im?portant type of national defense and military prod?uct,possesses the advantages of high energy densi?ty,great mechanical properties and superior securi?ty. However,cracks often occur in the process of forming,machining and post?processing,leading to performance degradations of the explosives[1-5].Ultra?sonic inspection is an important nondestructive tech?nique for assessment of many types of macro and mi?cro defects in contrast to the industrial computed to?mography[6-8]. Given that cracks usually initiate from the surface,ultrasonic surface wave seems to be rel?atively more feasible and appropriate for the evalua?tion of crack depth in PBXs providing more sensitivi?ty than others with most of the wave energy con?fined in the vicinity of the free surface within a depth of the specimen under test[9].

    Researches have been reported on ultrasonic crack detection in plates,but references to crack depth determination with respect to curved surfaces of PBXs are still limited for the moment[10]. In a re?cent investigation,C. He[11]presented a quantitative characterization method for estimating surface crack depth employing electromagnetic acoustic transduc?ers. Other studies have also demonstrated consider?able results on the interaction of Rayleigh waves and cracks[12]. However,the above studies investigate cracks only on the planar surface,while distinctive propagation patterns of surface waves are supposed to turn up concerning curved surfaces. To date,few works on crack detection on curved surfaces using surface waves can be found.

    This paper aims to provide solutions of ultrason?ic method for quantitatively estimating crack depth on cylindrical and spherical surfaces of PBXs. An in?version method for crack depth estimation involving a benchmark of the transmission and reflection curves is proposed. By means of the finite element method(FEM),wave motion of surface wave on cy?lindrical and spherical PBXs are studied,along with the reflection and transmission waves resulted from the interaction of Rayleigh waves and cracks of dif?ferent depths. Results are analyzed considering the theoretical predictions. Comparison of the simula?tion results with certain references is made to ensure the credit of this numerical study.

    2 Theory

    2.1 Rayleigh surface wave

    Rayleigh waves are a type of surface wave which propagate near the surface of solids and de?crease exponentially in amplitude as distance from the surface increases. Many researches have studied the wave velocity of Rayleigh surface wave regard?ing different characteristics of materials[13]. The ma?jor axis of the elliptic orbits is vertical to the surface,and the in?plane motion of a particle is retrograde at shallow depths and prograde at greater depths. The earliest approximate expression was given bycR=(0.87 + 1.12ν)( 1+ν)·c2,wherec2is the shear wave velocity.

    2.2 Dispersion in curved surfaces

    The Rayleigh wave dispersion is usually observed when the properties of the materials alter with depth such as propagations in the earth and multi?layered structures. On the other hand,the influence of cur?vature radii on the wave number and wave velocity has been theoretically analyzed. The Rayleigh wave velocity in curved surface is given by[13]

    wherecRdenotes the Rayleigh wave velocity along a plane surface,ωis the angular frequency,ραandρβare the curvature radii parallel and perpendicular to the propagation direction,respectively.A,BandCare expressed as follows

    wherec1andc2are the longitudinal and transverse wave velocities,respectively,andλ,μare Lamé constants. It is noteworthy that the phase velocity of Rayleigh wave along a curved surface is dependent on the wave frequency and curvature radii besides the elastic constants and density of the material. It can be inferred that suitable frequencies of surface wave are supposed to be applied for different mea?surement scales.

    3 Model development

    The following simulation was implemented based on the COMSOL Multiphysics V5.3a. The sol?id mechanics module was employed for the Ray?leigh surface wave studies.

    mother took her to the kitchen. She filled three pots with water and placed each on a high fire. Soon the pots came to a boil. In the first, she placed carrots, in the second she placed eggs and in the last she placed ground coffee beans. She let them sit and boil without saying a word. In about twenty minutes, she turned off the burners. She fished the carrots out and placed them in a bowl. She pulled the eggs out and placed them in a bowl. Then she ladled the coffee out and placed it in a bowl.

    3.1 Geometry & material

    The FEM simulation was conducted for existing products of PBXs of specific sizes. The computation?al model as regards cylindrical surface was estab?lished by a two?dimensional geometry. The left side of the geometry was reserved for study of the surface wave propagation along cylindrical surface,and a crack perpendicular to the surface was set on the right side(45° direction). The width of the crack is constantly set to 0.1 mm.

    According to the specific material properties provided by the Institute of Chemical Materials(Chi?na Academy of Engineering Physics),the approxi?mately isotropic homogeneous material is defined as follows:the density is 1860 kg·m-3,the Young's modulus is 10 GPa and the Poisson's ratio is 0.32.Note that certain PBXs are made of a high?volume fraction of explosive particles and a very low?vol?ume fraction of binder,especially for the HMX(e.g.PBX 9501)type of PBX that comprises a 95% explo?sive particle and only 5% binder[14]. It has been proved that the explosive particles are linear elastic at or below room temperature[15]. In addition,the monoclinic crystals (explosive particles) are ran?domly oriented in PBX[15],so the isotropic elastic properties can be approximately used for this study.

    For the spherical surface,a two?dimensional ax?isymmetric model was built with the same material properties and similar geometry(a crack set on the 45° direction). Fig.1 shows the schematics of the above configurations. It should be noted that the crack drawn in the figure has been enlarged for clari?ty. The space around the crack is filled with air.

    3.2 Structural mechanics & mesh

    The ultrasonic point source was determined by prescribed displacement. A sequence of domain point probes was set on the surface of the geometric model as shown in Fig.1 to receive the ultrasonic sig?nals with respect to the normal and tangential dis?placement components. Given that surface waves are constrained to a certain depth,two meshes of different densities are deployed to the whole do?main. The mesh size of the outer surface layer isλ·10-1whereλis the wavelength,and a predefined extremely fine mesh is set up for the rest. The surface layer is twice as thick as the crack. In order to re?duce the wave interference and facilitate the analysis of the signal,a low reflection boundary was deter?mined for the inner surface of the models.

    Since the Rayleigh surface waves are excited by point sources determined on the outer faces of the models by prescribed displacements(Hanning win?dowed sin signals of multi?cycles,the maximum am?plitudes are 1× 10-7m)towards the normal direc?tion of the surfaces,ultrasonic waves propagating along nearly half space direction are introduced into the models. Except for the expected Rayleigh surface waves which propagate along the outer surfaces,the other inconsequential waves can be diminished to a great extent by the absorption of the low reflection boundary conditions. Note that this excitation meth?od is in certain degree an imitation of laser ultrason?ic excitation,which has been used by a few re?searchers with reasonable results[16-18]. In addition,the ultrasonic near?field effect,which is not perti?nent to the scope of this study,is diminutive enough to be neglected according to this particular mechan?ic load.

    3.3 Solver configurations

    The numerical simulation involves ultrasonic in?put signals of 0.25,0.5,1.0,1.5 MHz and 2.0 MHz.The time steps of the transient solver are set as 0.1/f,wherefdenotes the frequencies of the excited Ray?leigh waves. The above configurations can ensure moderate accuracy of solutions to the specified problems.

    Fig.1 Geometric models and physical setups

    4 Results & discussion

    4.1 In cylindrical surface

    According to the model developments,an input signal at 1 MHz was brought into the model from the above of the geometry. The prescribed displace?ment vertically downward causes surface waves to be generated on the top of the geometry and propa?gated towards left and right. Firstly,signals picked up by probes on the left side are presented in Fig.2a.As one can see,the normal component of the parti?cle displacement in the cylindrical surface decreases firstly to a certain degree,then rises symmetrically to almost the same values it has reached before. The minimum of the normal displacement is presumed to occur at the 45° position on the left side. In the meanwhile,the trends of decrease and increase are most likely to follow a linear pattern. The tangential displacements seem to be identical as shown in Fig.2b. The above phenomenon can be attributed to the curve of the Rayleigh surface wave guide. For the cylindrical surface,the radii of curvature of the outer surface isRfor direction parallel to the figure,and infinite for direction perpendicular to the figure.The wave velocity of Rayleigh surface wave at 1 MHz propagating in cylindrical surface of PBX is calculated as 1302.5 m·s-1.

    Fig.2 Rayleigh surface waves(1 MHz)received on the cy?lindrical surface at different positions

    Fig.3 Snapshots of Rayleigh surface waves(1 MHz)on the cylindrical surface at different times

    The wavelength of the input ultrasonic signal at the frequency of 1 MHz is about 1.3 mm,and the crack depth ranging from 0.1 mm to 1.1 mm was de?termined to achieve the above numerical simula?tions. The ultrasonic surface wave propagated along the curved cylindrical surface and interacted with the surface crack,resulting in the reflected and trans?mitted waves which were received at the fixed posi?tions(110° and 160°)for the investigation of assess?ment of crack depth. Fig.4 shows the maximum am?plitudes of the reflected waves and transmission waves with respect to different crack depth. As seen,the maximum amplitudes of the reflected waves in?crease with the crack depth within certain range(d≤0.4 mm). In a similar way,the maximum ampli?tudes of the transmission waves decrease with the crack depth withind≤0.4 mm. Therefore,the depth of a crack can be evaluated by measuring the maxi?mum amplitudes of the reflected and transmission waves and concurrently referring to a pre?measured baseline of the relation curves. Note that the maxi?mum amplitude of the reflected waves decreases when the crack depth increases above 0.6 mm,this can be due to the energy scattering of multi wave packets generated from the interaction with the crack.

    Fig.4 Relation between the maximum amplitudes of reflect?ed and transmission waves and the crack depth(input signal at 1 MHz)

    4.2 In spherical surface

    4.2.1 Signal analysis

    Regarding the spherical surface,similar FEM simulations are achieved adopting the 2D axisym?metric model presented in Fig. 1b,which can be used to enhance the calculation efficiency of the simulation(i.e. by reducing a 3D problem to 2D axi?symmetric one). This setup is believed to be reason?able with regard to serval published works[19-20]and the COMSOL Users Guide. First,an ultrasonic sig?nal was introduced into the model in the absence of the crack for study of the wave propagation. As shown in Fig.5,both the normal and tangential com?ponents of the particle displacement are altering while the wave propagates along the surface. Espe?cially,the normal displacement of the wave in the spherical surface decays and rises over an approxi?mate exponential trend,while the tangential dis?placement varies over an approximate linear trend.The minimums of the displacement components are both observed at the 0° position. The curvature of the spherical surface mainly accounts for the above amplitude changes. In addition to the excitation and reception of a signal at the frequency of 0.25 MHz as shown in Fig.5,signal at other frequencies(0.5,1,1.5 MHz,and 2 MHz)were also introduced into the model for comparison. An interesting fact is that the amplitude of the Rayleigh surface wave in?creased back to larger values after propagating through the equator of the model (0° position)when the frequency of the input signal is 0.25 MHz and 0.5 MHz,while the amplitude of the surface wave increased much less after that when the fre?quency is above 1 MHz. Especially,the amplitude of the surface wave at the frequency of 2 MHz in?creased a few while propagated to the -80° position compared to that at the 0° position. Besides,the number of cycles of the time domain signal observed from the 0° to -80° positions increased. This phe?nomenon,which is not presented in figure here,can be attributed to the geometric dispersion of high?frequency ultrasonic signals,embodying the law of conservation of energy.

    Fig.5 Rayleigh surface waves(0.25 MHz)received on the spherical surface at different positions

    4.2.2 Motion pattern

    Fig.6 Analysis of geometrical dispersion impact on the am?plitude of Rayleigh surface wave

    In order to understand the motion of particles in more detail,on the other hand,the amplitude ratios of the received surface waves in the normal and tan?gential directions at different positions are calculat?ed. As shown in Fig.7a,the amplitude ratio(normal component to tangential component) of the Ray?leigh surface wave decreased from the wave source to 45° position,and increased to 1.6 at the equator of the model(0° position). It is interesting to find that the trends of the four different input signals are consis?tent in terms of the amplitude ratio. In addition,the amplitude ratio reached 1.0 when the surface wave propagated to the 45° position,which suggests that the orbit of the surface particle there is a standard cir?cle.More specifically and vividly,F(xiàn)ig.7b demonstrat?ed the elliptical orbits of particle motion at different positions. It is noted that the range of motion of the surface particle reduces since the Rayleigh wave dis?perses from the excitation point source.

    Fig.7 Motional pattern of the surface waves

    The skin effect of the Rayleigh surface wave in the spherical surface is also verified by setting a se?ries of domain point probes at the 45° position un?der the surface of the model with different depths. As presented in Fig.8,the normal amplitudes of the ul?trasonic signals at different frequencies decreases along the depth beneath the surface of the model. Es?pecially,the higher the frequency of the input sig?nal,the smaller the area Rayleigh surface wave can influence. As one can see,the amplitude of the re?ceived signal detected 1 mm beneath the surface is rather smaller for input signal at the frequency of 1.5 MHz than others. This fact provides guidance for the frequency selection in the estimation of surface crack depth in the subsequent simulations and future experiments. In addition,the result is in accordance with theory relating to surface wave that the primary wave field is within a wavelength long beneath the surface[13].

    Fig.8 Skin effect of the Rayleigh surface waves at different frequencies

    4.2.3 Crack depth inversion

    Based on the above analyses and the goal to as?sess the crack depth in PBXs,ultrasonic signals at frequencies of 0.25,0.5,1,1.5 MHz and 2 MHz are brought into the model shown in Fig.1b to inves?tigate the wave reflections and transmissions,respec?tively. The surface wave detected at the 80°,70°and 60° are extracted for the analysis of the wave re?flections,and signals received at the 30°,20° and 10° are processed for study of the wave transmis?sions. As shown in Fig.9,the normal amplitudes of the reflected and transmitted waves at different fre?quencies vary with the crack depth. Note that the ab?scissa is labeled asd/λsince the FEM results suggest the wave reflections share almost the same pattern except for slight discrepancies as the frequency of the input signal varies,as well as the wave transmis?sions. As seen from Fig.9a,amplitudes of signals re?ceived near the ultrasonic source(80° for instance)are larger than that away from the source(60° for in?stance),which suggests that one can set the receiv?ing probes befittingly away from the equator of the future specimen for higher signal?to?noise ratio. Be?sides,the relation of approximately monotonic de?creasing variation between the amplitude of the transmitted wave and the crack depth shown in Fig.9b seems to provide a more feasible and practi?cal way to estimate the crack depth in the spherical surface than the reflected waves. The amplitudes of the reflected and transmitted waves remain almost the same(call it a saturation)when the crack depth reaches to a certain degree,and the local irregulari?ties shown in Fig.9a may be attributed to the interfer?ences of the back scattered waves. As for the tangen?tial components of the reflected waves,F(xiàn)ig.10 shows a similar relationship to Fig. 9a. From a practical point of view,however,the normal displacement is most likely to be measured in experiments other than the tangential one,even though the tangential component possesses a less variation of amplitude while propagating along the curved surface.

    Fig.9 Relation between the crack depth and the normal am?plitudes of the waves received at different positions. d,λ de?note the crack depth and wavelength of the excited surface wave,respectively

    Fig.10 Relation between the crack depth and the tangential amplitudes of the reflected waves at different positions. d,λ denote the crack depth and wavelength of the excited surface wave,respectively

    For sufficient accuracy of FEM results ranging from 0 to about 0.3d/λ(ddenotes the crack depth),more simulations have been conducted as shown in Fig.11. The frequency of the input signal was deter?mined as 0.25 MHz(wavelength 12 mm). As one can see,due to the relative monotonic variation trend(d/λ <0.23,F(xiàn)ig. 11b). The relation between the crack depth and the normal amplitude of the transmitted wave can be utilized for crack depth esti?mation. In addition,the reflected waves can also be referred to for reliability considerations. These wave reflection and transmission curves indicate that one can measure a series of benchmarks as references for depth estimation of unknown cracks in the curved surface,and different frequencies can be used for characterization at different length scales. Thus,the inversion method for the crack depth estimation in curved surfaces of this PBX can be described as fol?lows:

    Step 1:For a PBX with unknown depth of crack,introduce a relatively very low?frequency ultrasonic Rayleigh surface wave into the material (e. g.0.1 MHz for a crack of which depth is below 1mm to ensured/λ <0.04). In other words,the detection starts from using a long wavelength surface wave to observe the reflected and transmitted waves.

    Step 2:By frequency tuning(increase the fre?quency of the input signal),the wavelength decreas?es,and the amplitude of transmitted wave drops,along with the increase of amplitude of the reflected wave.

    Step 3:When the amplitude of the reflected wave reaches its first extremum(e.g. a local maxi?mum as shown in Fig.11a),continue the frequency increasing with attention.

    Fig.11 Relation between the crack depth and the normal am?plitudes of the waves received at different positions. The fre?quency of the input signal is 0.25 MHz(wavelength 12 mm)

    Step 4:When the amplitude of the reflected wave reaches its first maximum as shown in Fig.9a,the crack depth can be estimated as 0.23λ. Noted that this can be confirmed by repeated implementa?tion of a large range of frequency,as well as the ref?erence of amplitude of the transmitted wave.

    It is important to note that the above simulations have not taken the material attenuation into consid?eration for investigation of surface acoustic waves in?fluenced by geometry and material separately. How?ever,the material attenuation could be determined as constant when the propagation path and proper?ties of the surface wave remain consistent for experi?mental measurements. Besides,the results are in cer?tain way similar to a published work which was con?ducted in a cracked plate by electromagnetic acous?tic transducer[15]. Thus,the proposed crack estima?tion method can potentially be promising with credit.

    5 Conclusions

    (1)The surface wave propagations in cylindri?cal and spherical surfaces of PBXs are investigated.The general disperse patterns of Rayleigh surface wave motion along these curved surfaces are pre?sented and analyzed.

    (2)The wave reflections and transmissions with a single crack by introducing ultrasonic signals at dif?ferent frequencies are recorded and processed for study of the influences on surface waves by the in?duced rectangle cracks with various depths.

    (3) The numerical simulations indicated that the proposed crack depth estimation method involv?ing transmission waves and reflected waves is poten?tially viable tools to quantitatively estimate the depth of cracks in the curved surface of PBXs.

    (4)Future related experiments are of great sig?nificance to achieve the crack location and depth measurement in curved surfaces with the method and results provided in this paper.

    av女优亚洲男人天堂| 一a级毛片在线观看| 在线观看免费视频日本深夜| 搞女人的毛片| 午夜视频国产福利| 欧美潮喷喷水| 蜜桃久久精品国产亚洲av| 最后的刺客免费高清国语| 99九九线精品视频在线观看视频| 九九热线精品视视频播放| 久久精品夜夜夜夜夜久久蜜豆| 国产三级中文精品| 精品人妻偷拍中文字幕| 中文亚洲av片在线观看爽| 国产精品一区二区性色av| 亚洲成人中文字幕在线播放| 听说在线观看完整版免费高清| 丝袜美腿在线中文| 最新在线观看一区二区三区| 波多野结衣巨乳人妻| 黄色日韩在线| av国产免费在线观看| 亚洲精品粉嫩美女一区| 免费观看的影片在线观看| 久久精品国产99精品国产亚洲性色| 国产在线精品亚洲第一网站| 国产精品1区2区在线观看.| 国产三级中文精品| 国产av不卡久久| 色播亚洲综合网| 在线免费观看不下载黄p国产| 国内久久婷婷六月综合欲色啪| 亚洲欧美成人精品一区二区| 天天躁夜夜躁狠狠久久av| 男女做爰动态图高潮gif福利片| a级毛片a级免费在线| 人妻少妇偷人精品九色| 国产男人的电影天堂91| 国产极品精品免费视频能看的| 久久久久久久久中文| 亚洲欧美中文字幕日韩二区| 日韩欧美国产在线观看| 国产中年淑女户外野战色| 99riav亚洲国产免费| 日本免费a在线| 我要看日韩黄色一级片| 一级毛片久久久久久久久女| 深爱激情五月婷婷| 亚洲av免费在线观看| 国产乱人视频| 国产伦一二天堂av在线观看| 久久99热6这里只有精品| 最近视频中文字幕2019在线8| 久久人人精品亚洲av| 成人高潮视频无遮挡免费网站| 99久国产av精品| 六月丁香七月| 国产高清视频在线观看网站| 亚洲乱码一区二区免费版| 欧美激情久久久久久爽电影| 啦啦啦啦在线视频资源| 国产精品嫩草影院av在线观看| 99热6这里只有精品| 色av中文字幕| 欧美国产日韩亚洲一区| 内地一区二区视频在线| 久久久久久久久中文| 波野结衣二区三区在线| 在线观看免费视频日本深夜| 精品久久久噜噜| 69人妻影院| 一级毛片电影观看 | 亚洲精品日韩av片在线观看| 国产精品无大码| 日韩中字成人| h日本视频在线播放| 日韩在线高清观看一区二区三区| 国产麻豆成人av免费视频| 91久久精品国产一区二区三区| 久久久久久久久久黄片| 国产精品日韩av在线免费观看| 日本撒尿小便嘘嘘汇集6| 一级av片app| 成年av动漫网址| 男人狂女人下面高潮的视频| 国产亚洲欧美98| 久久精品夜色国产| 热99re8久久精品国产| 日韩 亚洲 欧美在线| 日本撒尿小便嘘嘘汇集6| 国产精品福利在线免费观看| 美女大奶头视频| 日韩一本色道免费dvd| 亚洲人成网站在线观看播放| 一夜夜www| 国产精品日韩av在线免费观看| 哪里可以看免费的av片| 国产国拍精品亚洲av在线观看| 在线免费观看的www视频| 91麻豆精品激情在线观看国产| 午夜久久久久精精品| 99国产精品一区二区蜜桃av| 国产熟女欧美一区二区| 少妇高潮的动态图| 亚洲中文日韩欧美视频| 国产成人freesex在线 | 高清毛片免费看| 别揉我奶头~嗯~啊~动态视频| 日韩av在线大香蕉| 久久久久久国产a免费观看| 亚洲成av人片在线播放无| 国产单亲对白刺激| 欧美性感艳星| 中文亚洲av片在线观看爽| 国产免费男女视频| 免费人成在线观看视频色| 日韩欧美精品免费久久| 久久亚洲精品不卡| 99视频精品全部免费 在线| 久久人妻av系列| 国产淫片久久久久久久久| 俄罗斯特黄特色一大片| 日本-黄色视频高清免费观看| 深夜a级毛片| 在线播放无遮挡| 日本黄色片子视频| 国产精品嫩草影院av在线观看| 欧美zozozo另类| 一夜夜www| 久99久视频精品免费| 亚洲图色成人| 国产精品1区2区在线观看.| 看黄色毛片网站| 成年版毛片免费区| 亚洲人成网站在线播放欧美日韩| 日日啪夜夜撸| 男女做爰动态图高潮gif福利片| 久久久久久大精品| 久久精品影院6| 搡老妇女老女人老熟妇| 噜噜噜噜噜久久久久久91| 亚洲国产精品成人综合色| 亚洲高清免费不卡视频| 久久久久久久午夜电影| 激情 狠狠 欧美| 午夜福利视频1000在线观看| 久久久国产成人精品二区| 成人美女网站在线观看视频| 欧美+日韩+精品| 日本黄大片高清| 色噜噜av男人的天堂激情| 久久久久久久午夜电影| 变态另类成人亚洲欧美熟女| 国产高清不卡午夜福利| 国产男靠女视频免费网站| 国产高清视频在线观看网站| 亚洲五月天丁香| 亚洲真实伦在线观看| 亚洲精品国产av成人精品 | 免费无遮挡裸体视频| 神马国产精品三级电影在线观看| 国产男靠女视频免费网站| 亚洲真实伦在线观看| 人妻久久中文字幕网| 欧美高清性xxxxhd video| 九九爱精品视频在线观看| 精品久久久久久久久久免费视频| 久久精品国产鲁丝片午夜精品| 免费黄网站久久成人精品| 可以在线观看毛片的网站| 亚洲经典国产精华液单| 人人妻,人人澡人人爽秒播| 五月玫瑰六月丁香| 午夜a级毛片| 国产久久久一区二区三区| 九九爱精品视频在线观看| 午夜视频国产福利| 精品不卡国产一区二区三区| 亚洲成人精品中文字幕电影| 国产一区二区亚洲精品在线观看| 高清毛片免费观看视频网站| 可以在线观看的亚洲视频| 精品99又大又爽又粗少妇毛片| 99久久精品一区二区三区| 波野结衣二区三区在线| 午夜精品一区二区三区免费看| 国产午夜精品论理片| 日韩欧美 国产精品| 黄片wwwwww| 床上黄色一级片| 91久久精品国产一区二区成人| 在线免费观看不下载黄p国产| 日本一本二区三区精品| 国产精品乱码一区二三区的特点| 久久精品国产亚洲av涩爱 | 中文在线观看免费www的网站| 欧美日韩国产亚洲二区| 亚洲av五月六月丁香网| 国产爱豆传媒在线观看| 亚洲美女视频黄频| 中国国产av一级| 69人妻影院| 麻豆乱淫一区二区| 久久精品国产鲁丝片午夜精品| 日韩成人伦理影院| 夜夜看夜夜爽夜夜摸| 深夜a级毛片| 精品久久久久久久久av| 亚洲国产精品成人综合色| 亚洲人与动物交配视频| 91精品国产九色| 少妇的逼水好多| 伊人久久精品亚洲午夜| 免费av不卡在线播放| 国产爱豆传媒在线观看| 男人舔女人下体高潮全视频| 女生性感内裤真人,穿戴方法视频| а√天堂www在线а√下载| 国产成人a∨麻豆精品| 日本爱情动作片www.在线观看 | 国产免费男女视频| .国产精品久久| 久久久久久九九精品二区国产| 蜜臀久久99精品久久宅男| av国产免费在线观看| 亚洲国产欧洲综合997久久,| 亚洲欧美日韩卡通动漫| 99久久精品国产国产毛片| 欧美xxxx性猛交bbbb| 日本撒尿小便嘘嘘汇集6| 中文字幕免费在线视频6| 午夜精品国产一区二区电影 | 久久久久久久久久成人| 天天一区二区日本电影三级| 久久久久国产网址| 在线观看av片永久免费下载| 国产高清激情床上av| 亚洲av免费在线观看| 美女内射精品一级片tv| 国产中年淑女户外野战色| 桃色一区二区三区在线观看| 成年免费大片在线观看| 91精品国产九色| 中文在线观看免费www的网站| 国内精品久久久久精免费| 成人美女网站在线观看视频| 免费一级毛片在线播放高清视频| 免费av观看视频| 日日摸夜夜添夜夜添小说| 色5月婷婷丁香| 伦理电影大哥的女人| 精品一区二区三区视频在线观看免费| 男人舔奶头视频| 男人舔女人下体高潮全视频| 免费高清视频大片| 我要看日韩黄色一级片| videossex国产| 国产精品美女特级片免费视频播放器| 少妇熟女欧美另类| 啦啦啦啦在线视频资源| 国产91av在线免费观看| 久久人妻av系列| 啦啦啦韩国在线观看视频| 天堂av国产一区二区熟女人妻| 久久久久国产网址| 大型黄色视频在线免费观看| 村上凉子中文字幕在线| 国内揄拍国产精品人妻在线| 亚洲精品粉嫩美女一区| av视频在线观看入口| 老司机影院成人| 中文字幕精品亚洲无线码一区| 99热这里只有是精品50| 赤兔流量卡办理| 天天一区二区日本电影三级| 亚洲一区二区三区色噜噜| 搡老岳熟女国产| 亚州av有码| 免费观看人在逋| 久久午夜亚洲精品久久| 久久精品国产亚洲网站| 久久精品国产鲁丝片午夜精品| 午夜日韩欧美国产| 国产精品亚洲一级av第二区| 12—13女人毛片做爰片一| 村上凉子中文字幕在线| 国产av不卡久久| or卡值多少钱| 一级毛片我不卡| 能在线免费观看的黄片| 亚洲天堂国产精品一区在线| 久久久久久久久中文| 中文亚洲av片在线观看爽| 亚洲成人中文字幕在线播放| 色综合站精品国产| 日韩成人av中文字幕在线观看 | 亚洲精华国产精华液的使用体验 | 精品少妇黑人巨大在线播放 | 久久6这里有精品| 亚洲欧美中文字幕日韩二区| 国产在线精品亚洲第一网站| 日韩精品中文字幕看吧| 亚洲精品一区av在线观看| 99久久中文字幕三级久久日本| 男女之事视频高清在线观看| 国产淫片久久久久久久久| 中文在线观看免费www的网站| 小说图片视频综合网站| 国产伦精品一区二区三区视频9| 91在线精品国自产拍蜜月| 麻豆av噜噜一区二区三区| 天堂av国产一区二区熟女人妻| 日本免费一区二区三区高清不卡| 天堂√8在线中文| 国产成人a∨麻豆精品| 蜜桃久久精品国产亚洲av| 熟妇人妻久久中文字幕3abv| 村上凉子中文字幕在线| 噜噜噜噜噜久久久久久91| 国产精品一区www在线观看| 又黄又爽又刺激的免费视频.| 日本成人三级电影网站| 日韩制服骚丝袜av| 女人十人毛片免费观看3o分钟| 午夜福利成人在线免费观看| 国产久久久一区二区三区| 免费观看精品视频网站| 日本-黄色视频高清免费观看| 97超视频在线观看视频| 成人特级av手机在线观看| 最近中文字幕高清免费大全6| 午夜精品在线福利| 欧美丝袜亚洲另类| 久久久久久久久久黄片| 欧美性猛交黑人性爽| 成人性生交大片免费视频hd| 中文字幕av在线有码专区| 国内精品美女久久久久久| 国产亚洲91精品色在线| 国产真实乱freesex| 免费观看精品视频网站| 日韩人妻高清精品专区| 亚洲av成人av| 久久久久九九精品影院| 国产色爽女视频免费观看| 白带黄色成豆腐渣| 五月伊人婷婷丁香| 免费观看人在逋| 免费看日本二区| 亚洲国产精品合色在线| 丰满的人妻完整版| 亚洲婷婷狠狠爱综合网| 国产精品永久免费网站| 小说图片视频综合网站| 深夜精品福利| 午夜老司机福利剧场| 内地一区二区视频在线| 国产精品99久久久久久久久| 在线观看av片永久免费下载| 九色成人免费人妻av| 国产中年淑女户外野战色| 人妻制服诱惑在线中文字幕| 久久99热这里只有精品18| 能在线免费观看的黄片| 日韩欧美在线乱码| 日本色播在线视频| 日韩欧美精品免费久久| 国产激情偷乱视频一区二区| 国内少妇人妻偷人精品xxx网站| 亚洲最大成人av| 免费看日本二区| 大香蕉久久网| 国产欧美日韩精品亚洲av| 欧美中文日本在线观看视频| 九九在线视频观看精品| 国产高清视频在线观看网站| 亚洲自偷自拍三级| 久久久精品欧美日韩精品| 秋霞在线观看毛片| 一个人看视频在线观看www免费| 搡女人真爽免费视频火全软件 | 日本在线视频免费播放| 插逼视频在线观看| 丰满人妻一区二区三区视频av| 最近的中文字幕免费完整| 久久精品国产亚洲av涩爱 | 国产蜜桃级精品一区二区三区| 免费观看的影片在线观看| 国产av麻豆久久久久久久| 国产老妇女一区| 国产蜜桃级精品一区二区三区| 18禁在线播放成人免费| 午夜激情欧美在线| 午夜视频国产福利| 亚洲精品亚洲一区二区| 久久精品综合一区二区三区| 乱系列少妇在线播放| 中文字幕av在线有码专区| 精品人妻偷拍中文字幕| 女生性感内裤真人,穿戴方法视频| 精品一区二区三区视频在线| 欧美潮喷喷水| 两个人的视频大全免费| 国产欧美日韩精品亚洲av| 日韩 亚洲 欧美在线| 看黄色毛片网站| 日本与韩国留学比较| 欧美日韩乱码在线| 亚洲欧美精品综合久久99| 搞女人的毛片| 亚洲乱码一区二区免费版| 日韩欧美免费精品| 一卡2卡三卡四卡精品乱码亚洲| 日本熟妇午夜| 亚洲精品成人久久久久久| 嫩草影院新地址| 人妻久久中文字幕网| 亚洲av五月六月丁香网| 男女下面进入的视频免费午夜| 国产精品日韩av在线免费观看| 日韩 亚洲 欧美在线| 免费av观看视频| 亚洲国产精品sss在线观看| 日韩av在线大香蕉| 一本精品99久久精品77| 国产视频一区二区在线看| 少妇的逼好多水| 国产69精品久久久久777片| 亚洲欧美日韩东京热| 久久精品国产亚洲av香蕉五月| videossex国产| 国产高清三级在线| 成人毛片a级毛片在线播放| 国产成人福利小说| 少妇人妻精品综合一区二区 | 综合色丁香网| 99热这里只有精品一区| 一级黄色大片毛片| 波多野结衣高清作品| 一进一出抽搐动态| 成人永久免费在线观看视频| 大型黄色视频在线免费观看| 中文字幕熟女人妻在线| 免费看日本二区| 国产精品1区2区在线观看.| 久久精品国产清高在天天线| 亚洲熟妇中文字幕五十中出| 国国产精品蜜臀av免费| 三级毛片av免费| 成人高潮视频无遮挡免费网站| 国产精品一及| 亚洲人与动物交配视频| 国产精品国产高清国产av| 国产精品免费一区二区三区在线| 久久久色成人| 午夜免费男女啪啪视频观看 | 给我免费播放毛片高清在线观看| 又黄又爽又刺激的免费视频.| 日韩av不卡免费在线播放| 欧美日韩乱码在线| 国产大屁股一区二区在线视频| 在线观看av片永久免费下载| 黄色视频,在线免费观看| 99久国产av精品| 99视频精品全部免费 在线| 国产一区二区在线av高清观看| 中出人妻视频一区二区| 亚洲成人av在线免费| 精品久久国产蜜桃| 老司机午夜福利在线观看视频| 高清午夜精品一区二区三区 | 日韩精品有码人妻一区| 国产片特级美女逼逼视频| 亚洲无线观看免费| 久久精品国产清高在天天线| 亚洲人与动物交配视频| 成年版毛片免费区| 哪里可以看免费的av片| 极品教师在线视频| 精品一区二区三区视频在线观看免费| 最近视频中文字幕2019在线8| 国产精品久久电影中文字幕| 一区二区三区高清视频在线| 日本欧美国产在线视频| 国语自产精品视频在线第100页| 亚洲熟妇中文字幕五十中出| 成人二区视频| 三级国产精品欧美在线观看| 免费看美女性在线毛片视频| 国产在视频线在精品| 国产精品福利在线免费观看| 能在线免费观看的黄片| h日本视频在线播放| 深爱激情五月婷婷| 老师上课跳d突然被开到最大视频| 精品人妻偷拍中文字幕| 丝袜美腿在线中文| 亚洲中文日韩欧美视频| 看非洲黑人一级黄片| 别揉我奶头~嗯~啊~动态视频| 人妻制服诱惑在线中文字幕| 男女边吃奶边做爰视频| 别揉我奶头 嗯啊视频| 欧美+亚洲+日韩+国产| 最后的刺客免费高清国语| 搡老熟女国产l中国老女人| 亚洲综合色惰| 免费看av在线观看网站| av在线蜜桃| 国产精品久久久久久av不卡| 日韩在线高清观看一区二区三区| 久久久久久久久大av| 免费在线观看成人毛片| 国产亚洲精品久久久久久毛片| 超碰av人人做人人爽久久| 国内精品久久久久精免费| 亚洲精品日韩av片在线观看| 可以在线观看的亚洲视频| 三级男女做爰猛烈吃奶摸视频| 国产精品女同一区二区软件| 天堂网av新在线| 欧美bdsm另类| 色视频www国产| 国产精品综合久久久久久久免费| 日本免费a在线| 国内精品久久久久精免费| 国内揄拍国产精品人妻在线| 国产黄色视频一区二区在线观看 | 狂野欧美激情性xxxx在线观看| 国产在视频线在精品| 国产乱人视频| 久久午夜亚洲精品久久| 看十八女毛片水多多多| 免费看美女性在线毛片视频| 99精品在免费线老司机午夜| 99热网站在线观看| 成熟少妇高潮喷水视频| av天堂在线播放| 简卡轻食公司| 欧美日韩一区二区视频在线观看视频在线 | 精品福利观看| or卡值多少钱| 国产欧美日韩精品一区二区| 国产伦精品一区二区三区视频9| 国产精品一二三区在线看| 成人av一区二区三区在线看| 久久6这里有精品| 亚洲国产日韩欧美精品在线观看| 日本精品一区二区三区蜜桃| 亚洲不卡免费看| 精品人妻视频免费看| 一级黄片播放器| 麻豆久久精品国产亚洲av| 国产成人a∨麻豆精品| 午夜日韩欧美国产| 国产男人的电影天堂91| 日本在线视频免费播放| eeuss影院久久| 麻豆av噜噜一区二区三区| 中文字幕av成人在线电影| 亚洲欧美清纯卡通| 国产精品1区2区在线观看.| 校园人妻丝袜中文字幕| 午夜免费男女啪啪视频观看 | 女的被弄到高潮叫床怎么办| 嫩草影视91久久| 1024手机看黄色片| 欧美三级亚洲精品| 国产成人a∨麻豆精品| 一个人看视频在线观看www免费| 国产精品爽爽va在线观看网站| 中文字幕久久专区| 国产精品一及| 国产av一区在线观看免费| 成人国产麻豆网| 国产精品伦人一区二区| 久久欧美精品欧美久久欧美| 欧美xxxx性猛交bbbb| 欧美一区二区精品小视频在线| 国产精华一区二区三区| 国产熟女欧美一区二区| 成人无遮挡网站| 国产毛片a区久久久久| 最近视频中文字幕2019在线8| 一级毛片aaaaaa免费看小| 亚洲精品色激情综合| 国产大屁股一区二区在线视频| 国产蜜桃级精品一区二区三区| 别揉我奶头 嗯啊视频| 男女下面进入的视频免费午夜| 成年女人毛片免费观看观看9| 国产精品无大码| 欧美丝袜亚洲另类| 国产乱人视频| 久99久视频精品免费| 又粗又爽又猛毛片免费看| 精品人妻熟女av久视频| 国内少妇人妻偷人精品xxx网站| 麻豆成人午夜福利视频| 麻豆久久精品国产亚洲av| 日韩亚洲欧美综合| 99久久久亚洲精品蜜臀av| 精品人妻熟女av久视频| 男人舔奶头视频| 麻豆乱淫一区二区| 欧美日本亚洲视频在线播放| 亚洲成a人片在线一区二区| 亚洲av中文av极速乱| 久久99热6这里只有精品| 欧美日本亚洲视频在线播放| 最新在线观看一区二区三区| 91av网一区二区| 午夜精品国产一区二区电影 | 国产一区二区三区在线臀色熟女| 人人妻人人看人人澡|