• <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| 成人高潮视频无遮挡免费网站| 国产日本99.免费观看| 最近的中文字幕免费完整| 18禁在线无遮挡免费观看视频 | 搡老妇女老女人老熟妇| 99在线人妻在线中文字幕| 日韩欧美三级三区| 国产精品av视频在线免费观看| 日本色播在线视频| 国产高清视频在线播放一区| 三级国产精品欧美在线观看| 我的老师免费观看完整版| 在线观看免费视频日本深夜| 日韩,欧美,国产一区二区三区 | 色综合站精品国产| 久久中文看片网| 久久久久久久久大av| 久久久久久久久大av| 99热只有精品国产| 久久精品影院6| 精品久久国产蜜桃| 变态另类丝袜制服| 长腿黑丝高跟| 日韩欧美三级三区| 日韩制服骚丝袜av| 黄片wwwwww| 亚洲欧美日韩无卡精品| 亚洲成人中文字幕在线播放| av在线观看视频网站免费| 亚洲精品亚洲一区二区| 国产女主播在线喷水免费视频网站 | 一级黄色大片毛片| 三级经典国产精品| 婷婷精品国产亚洲av| 亚洲最大成人中文| 不卡视频在线观看欧美| 精品福利观看| 亚洲精品一卡2卡三卡4卡5卡| 99久久九九国产精品国产免费| 色综合色国产| 插逼视频在线观看| 成人av一区二区三区在线看| 我的老师免费观看完整版| 简卡轻食公司| 欧美最新免费一区二区三区| 日本一二三区视频观看| 亚洲精品乱码久久久v下载方式| 午夜福利在线观看吧| 97热精品久久久久久| 18禁在线无遮挡免费观看视频 | 国产精品人妻久久久久久| 国产成人一区二区在线| 国产毛片a区久久久久| 少妇的逼水好多| 人妻丰满熟妇av一区二区三区| 精品一区二区三区av网在线观看| 免费观看精品视频网站| 日韩一区二区视频免费看| 99热这里只有是精品50| 欧美国产日韩亚洲一区| av黄色大香蕉| 啦啦啦观看免费观看视频高清| 国产精品野战在线观看| 99久久精品一区二区三区| 精品久久久久久久人妻蜜臀av| 亚洲中文日韩欧美视频| 久久亚洲国产成人精品v| 特级一级黄色大片| 18禁在线播放成人免费| 在线免费观看的www视频| 精品久久国产蜜桃| av在线播放精品| 美女内射精品一级片tv| 国产亚洲av嫩草精品影院| 亚洲欧美清纯卡通| 欧美色视频一区免费| 成年免费大片在线观看| 精品久久久久久成人av| 级片在线观看| 丝袜喷水一区| 看片在线看免费视频| 国产精品亚洲美女久久久| 国产白丝娇喘喷水9色精品| av在线蜜桃| 精品人妻视频免费看| www.色视频.com| 少妇的逼好多水| 中国国产av一级| av视频在线观看入口| 毛片女人毛片| 99热这里只有是精品在线观看| 欧美极品一区二区三区四区| 国产精品久久久久久久久免| 午夜久久久久精精品| 十八禁国产超污无遮挡网站| 午夜视频国产福利| 日韩成人伦理影院| 老熟妇乱子伦视频在线观看| 免费av不卡在线播放| 成年女人看的毛片在线观看| 中文字幕av成人在线电影| 国产一区二区三区av在线 | 日产精品乱码卡一卡2卡三| 深夜a级毛片| 色综合亚洲欧美另类图片| 国产一区二区在线av高清观看| 少妇裸体淫交视频免费看高清| 亚洲五月天丁香| 午夜爱爱视频在线播放| 国产日本99.免费观看| 内地一区二区视频在线| 99热这里只有是精品50| 国产精品亚洲美女久久久| 变态另类成人亚洲欧美熟女| 中文字幕av成人在线电影| 男人舔女人下体高潮全视频| 亚洲一级一片aⅴ在线观看| 欧美一区二区国产精品久久精品| 久久6这里有精品| 久久久久九九精品影院| 亚洲人与动物交配视频| 国产久久久一区二区三区| 日日摸夜夜添夜夜添小说| 日韩在线高清观看一区二区三区| 国产高清不卡午夜福利| 亚洲成人久久性| eeuss影院久久| 亚洲无线在线观看| 国产精品久久视频播放| 亚洲av免费高清在线观看| 一级黄片播放器| 亚洲美女搞黄在线观看 | 国产精品久久久久久精品电影| 国产精品,欧美在线| 色综合站精品国产| 一区二区三区免费毛片| 99久国产av精品| 免费人成在线观看视频色| 亚洲欧美日韩东京热| 国内少妇人妻偷人精品xxx网站| 久久亚洲精品不卡| 欧美性猛交╳xxx乱大交人| 亚洲国产色片| 中文字幕免费在线视频6| 午夜免费男女啪啪视频观看 | 性插视频无遮挡在线免费观看| 桃色一区二区三区在线观看| 亚洲av电影不卡..在线观看| 免费人成在线观看视频色| 婷婷精品国产亚洲av| 嫩草影视91久久| 精华霜和精华液先用哪个| 国产精品日韩av在线免费观看| 国模一区二区三区四区视频| 亚洲成人久久爱视频| 99久国产av精品| 国产探花在线观看一区二区| 特级一级黄色大片| 亚洲精品亚洲一区二区| 久久这里只有精品中国| 日韩欧美免费精品| 日韩欧美精品免费久久| 97超级碰碰碰精品色视频在线观看| 十八禁国产超污无遮挡网站| 久久久久久久亚洲中文字幕| 亚洲,欧美,日韩| 国产精品综合久久久久久久免费| 久久99热这里只有精品18| 国产精品爽爽va在线观看网站| 一级毛片我不卡| 如何舔出高潮| 亚洲欧美成人精品一区二区| 在线观看免费视频日本深夜| 波多野结衣高清无吗| 亚洲一级一片aⅴ在线观看| 少妇猛男粗大的猛烈进出视频 | 1000部很黄的大片| 国产精品无大码| 久久亚洲精品不卡| 校园人妻丝袜中文字幕| 欧美xxxx性猛交bbbb| 久久人妻av系列| 成人高潮视频无遮挡免费网站| ponron亚洲| 亚洲综合色惰| 久久6这里有精品| 成人鲁丝片一二三区免费| 国内精品久久久久精免费| 全区人妻精品视频| 精品一区二区免费观看| 乱人视频在线观看| 国产欧美日韩一区二区精品| 99热这里只有精品一区| 成人永久免费在线观看视频| 永久网站在线| 99久久精品热视频| 欧美性感艳星| 最近的中文字幕免费完整| 精品国内亚洲2022精品成人| 国产毛片a区久久久久| 国产av在哪里看| 国产一区亚洲一区在线观看| 一级毛片久久久久久久久女| 最近在线观看免费完整版| 精品午夜福利视频在线观看一区| 看片在线看免费视频| 中文亚洲av片在线观看爽| 亚洲欧美成人综合另类久久久 | 成人欧美大片| 国产老妇女一区| 少妇人妻一区二区三区视频| 变态另类丝袜制服| 国产精华一区二区三区| 欧美另类亚洲清纯唯美| 亚洲国产日韩欧美精品在线观看| 波多野结衣高清无吗| 久久精品影院6| 99热精品在线国产| 亚洲最大成人手机在线| 久久久a久久爽久久v久久| 特级一级黄色大片| 亚洲精品在线观看二区| 国产不卡一卡二| 搡老妇女老女人老熟妇| 桃色一区二区三区在线观看| 亚洲精华国产精华液的使用体验 | 欧美三级亚洲精品| 久久精品久久久久久噜噜老黄 | 乱码一卡2卡4卡精品| 亚洲丝袜综合中文字幕| 欧美日韩综合久久久久久| 男人舔女人下体高潮全视频| 男女之事视频高清在线观看| 在线免费观看不下载黄p国产| 午夜视频国产福利| 亚洲欧美精品自产自拍| 久久久成人免费电影| 女人十人毛片免费观看3o分钟| 亚洲一区高清亚洲精品| 少妇的逼水好多| 成人av一区二区三区在线看| 99久久中文字幕三级久久日本| 欧美成人一区二区免费高清观看| 午夜精品国产一区二区电影 | a级毛片a级免费在线| 日日干狠狠操夜夜爽| 欧美xxxx性猛交bbbb| 99精品在免费线老司机午夜| 国产极品精品免费视频能看的| 久久亚洲精品不卡| 欧美在线一区亚洲| 免费看日本二区| 人人妻人人看人人澡| 成人毛片a级毛片在线播放| 国产亚洲91精品色在线| 亚洲婷婷狠狠爱综合网| 色视频www国产| 日韩在线高清观看一区二区三区| ponron亚洲| or卡值多少钱| 日本撒尿小便嘘嘘汇集6| 成人精品一区二区免费| 午夜福利18| 在线观看美女被高潮喷水网站| 中国美女看黄片| 内地一区二区视频在线| 日韩国内少妇激情av| 婷婷色综合大香蕉| 一级毛片久久久久久久久女| 美女高潮的动态| 免费看光身美女| 搡老岳熟女国产| 日本a在线网址| 国产91av在线免费观看| 欧美不卡视频在线免费观看| 国产一区二区在线观看日韩| av视频在线观看入口| 一卡2卡三卡四卡精品乱码亚洲| 欧美激情国产日韩精品一区| 97人妻精品一区二区三区麻豆| av在线天堂中文字幕| 97超级碰碰碰精品色视频在线观看| 99精品在免费线老司机午夜| 国产成人精品久久久久久| 黄色欧美视频在线观看| 精品一区二区三区人妻视频| 亚洲熟妇熟女久久| 国产人妻一区二区三区在| 99久久精品一区二区三区| 国产三级中文精品| 99九九线精品视频在线观看视频| 免费看日本二区| 网址你懂的国产日韩在线| 少妇丰满av| 精品久久久久久久末码| 51国产日韩欧美| 欧美另类亚洲清纯唯美| 欧美人与善性xxx| 男人舔女人下体高潮全视频| 成人特级av手机在线观看| 嫩草影院入口| 国产91av在线免费观看| 又黄又爽又刺激的免费视频.| 亚洲欧美精品自产自拍| 日韩三级伦理在线观看| 久久久国产成人免费| 人人妻人人看人人澡| 97人妻精品一区二区三区麻豆| 欧美在线一区亚洲| 日本在线视频免费播放| 久久精品91蜜桃| 深夜精品福利| 欧美成人一区二区免费高清观看| 亚洲人成网站在线播| 一级黄色大片毛片| 欧美精品国产亚洲| 久99久视频精品免费| 欧美日韩国产亚洲二区| 国产高清激情床上av| 国产真实乱freesex| 亚洲国产日韩欧美精品在线观看| 国产成人91sexporn| 亚洲av免费高清在线观看| 久久亚洲国产成人精品v| 精品国产三级普通话版| 一级黄片播放器| 国产亚洲精品久久久久久毛片| 九色成人免费人妻av| 麻豆av噜噜一区二区三区| 日韩在线高清观看一区二区三区| 午夜福利在线在线| 嫩草影院新地址| 中文字幕免费在线视频6| 成年女人看的毛片在线观看| 久久久久久国产a免费观看| 欧美3d第一页| 亚洲专区国产一区二区| 日本精品一区二区三区蜜桃| 简卡轻食公司| 国产精品99久久久久久久久| 久久这里只有精品中国| 毛片女人毛片| 99精品在免费线老司机午夜| 亚洲激情五月婷婷啪啪| 亚洲精品国产成人久久av| 如何舔出高潮| 成年女人看的毛片在线观看| 国产精品日韩av在线免费观看| 亚洲一区高清亚洲精品| 精品不卡国产一区二区三区| 国产成人福利小说| 如何舔出高潮| 精品不卡国产一区二区三区| 日韩三级伦理在线观看| 美女 人体艺术 gogo| 亚洲欧美日韩东京热| 久久久久久久久久黄片| 在线免费观看不下载黄p国产| 欧美+日韩+精品| 97超视频在线观看视频| 成人av在线播放网站| 日韩,欧美,国产一区二区三区 | 老女人水多毛片| 成人毛片a级毛片在线播放| 美女xxoo啪啪120秒动态图| 亚洲第一电影网av| 免费在线观看影片大全网站| 深夜精品福利| 在线免费十八禁| 色5月婷婷丁香| 精品一区二区免费观看| 男人和女人高潮做爰伦理| 人妻丰满熟妇av一区二区三区| 久久久午夜欧美精品| 欧美成人a在线观看| 精品国内亚洲2022精品成人| 国产一区二区三区av在线 | 嫩草影院新地址| 老师上课跳d突然被开到最大视频| 亚洲第一区二区三区不卡| a级毛片a级免费在线| 国产精品一二三区在线看| .国产精品久久| 床上黄色一级片| 色视频www国产| 久久久久国产网址| 免费av毛片视频| 国产麻豆成人av免费视频| 亚洲欧美日韩高清在线视频| 亚洲精品亚洲一区二区| 亚洲av第一区精品v没综合| 乱码一卡2卡4卡精品| 男女啪啪激烈高潮av片| 一夜夜www| 特大巨黑吊av在线直播| 两个人的视频大全免费| 国产精品三级大全| 亚洲欧美精品综合久久99| 国产成人精品久久久久久| 久久久久久大精品| 一级黄片播放器| 国产成人freesex在线 | 欧美极品一区二区三区四区| 精品不卡国产一区二区三区| 一级毛片电影观看 | 大香蕉久久网| 天堂动漫精品| 国产69精品久久久久777片| 亚洲第一区二区三区不卡| 俄罗斯特黄特色一大片| 亚洲精品乱码久久久v下载方式| ponron亚洲| 亚洲国产精品成人综合色| 日韩一区二区视频免费看| avwww免费| 别揉我奶头~嗯~啊~动态视频| 伊人久久精品亚洲午夜| 97人妻精品一区二区三区麻豆| av在线播放精品| 国产激情偷乱视频一区二区| 欧美日本视频| 我的女老师完整版在线观看| 99热全是精品| 亚洲专区国产一区二区| 淫秽高清视频在线观看| 日本三级黄在线观看| 嫩草影视91久久| 亚洲va在线va天堂va国产| 久久久久久久久久成人| 久久久久久伊人网av| 99热这里只有是精品在线观看| 小说图片视频综合网站| 久久综合国产亚洲精品| 淫妇啪啪啪对白视频| 久久国产乱子免费精品| 在线免费十八禁| 国产精品国产高清国产av| 亚洲无线在线观看| 久久久久精品国产欧美久久久| 欧美日本亚洲视频在线播放| 精品一区二区三区视频在线| 久久久久国产精品人妻aⅴ院| 国产男人的电影天堂91| 性欧美人与动物交配| 在线观看66精品国产| 免费看美女性在线毛片视频| 丝袜喷水一区| avwww免费| 又爽又黄无遮挡网站| 成人特级黄色片久久久久久久| а√天堂www在线а√下载| 精品久久久久久久久久免费视频| 在线观看一区二区三区| 国产美女午夜福利| 成人特级黄色片久久久久久久| 亚洲七黄色美女视频| 精品久久久久久成人av| 国产精品福利在线免费观看| 亚洲最大成人手机在线| 一个人看的www免费观看视频| 国产精品三级大全| 又爽又黄a免费视频| 在线观看一区二区三区| 麻豆av噜噜一区二区三区| 国产精品综合久久久久久久免费| 性色avwww在线观看| 伦理电影大哥的女人| 免费看美女性在线毛片视频| 免费一级毛片在线播放高清视频| 国产精品一及| 欧美国产日韩亚洲一区| 国产 一区精品| 国产麻豆成人av免费视频| av.在线天堂| 久久久久国内视频| 色综合色国产| 免费看美女性在线毛片视频| 国产激情偷乱视频一区二区| 日本精品一区二区三区蜜桃| 亚洲不卡免费看| 无遮挡黄片免费观看| 午夜激情福利司机影院| 亚洲成人久久性| 日韩亚洲欧美综合| 国产成人精品久久久久久| 18禁在线无遮挡免费观看视频 | a级毛片免费高清观看在线播放| 日本一本二区三区精品| 国产高清激情床上av| 91精品国产九色| 人人妻人人澡欧美一区二区| 亚洲国产精品成人久久小说 | 国内精品一区二区在线观看| 日韩欧美免费精品| av在线亚洲专区| 国产一区二区三区在线臀色熟女| 国产伦精品一区二区三区四那| 亚洲18禁久久av| 国产aⅴ精品一区二区三区波| 在线观看一区二区三区| 国产单亲对白刺激| 亚洲美女黄片视频| 亚洲人与动物交配视频| 久久久久久九九精品二区国产| 高清毛片免费观看视频网站| 波野结衣二区三区在线| 亚洲av熟女| 成人高潮视频无遮挡免费网站| 少妇裸体淫交视频免费看高清| 国产精品不卡视频一区二区| 国产成人福利小说| 国产精品电影一区二区三区| 亚洲av成人av| 久久人妻av系列| 天堂动漫精品| 黄色欧美视频在线观看| 亚洲av免费在线观看| 成人美女网站在线观看视频| 国产男人的电影天堂91| 精品久久久噜噜| 长腿黑丝高跟| 免费观看精品视频网站| АⅤ资源中文在线天堂| 男女那种视频在线观看| 香蕉av资源在线| 久久精品国产鲁丝片午夜精品| 成人国产麻豆网| 国产av一区在线观看免费| 国产精品久久久久久精品电影| 亚洲不卡免费看| 免费电影在线观看免费观看| 男女那种视频在线观看| 五月玫瑰六月丁香| 成人综合一区亚洲| 欧美精品国产亚洲| 国产高清三级在线| 我的女老师完整版在线观看| 亚洲人成网站高清观看| 精品午夜福利在线看| 美女被艹到高潮喷水动态| 97在线视频观看| 看黄色毛片网站| 嫩草影院精品99| 伦精品一区二区三区| 校园人妻丝袜中文字幕| 婷婷亚洲欧美| 亚洲高清免费不卡视频| 亚洲成人av在线免费| 春色校园在线视频观看| 不卡一级毛片| 51国产日韩欧美| 国产av在哪里看| 欧美色欧美亚洲另类二区| 免费看a级黄色片| 欧美丝袜亚洲另类| 国内精品美女久久久久久| 好男人在线观看高清免费视频| 精品一区二区三区人妻视频| 人妻少妇偷人精品九色| 22中文网久久字幕| 欧美xxxx性猛交bbbb| 自拍偷自拍亚洲精品老妇| 中出人妻视频一区二区| 国产一区二区在线观看日韩| 免费黄网站久久成人精品| 色综合亚洲欧美另类图片| 十八禁国产超污无遮挡网站| 久久热精品热| 国产一区二区激情短视频| 国产精品永久免费网站| 在线免费观看的www视频| 亚洲在线自拍视频| 国产国拍精品亚洲av在线观看| 啦啦啦韩国在线观看视频| 国国产精品蜜臀av免费| 免费一级毛片在线播放高清视频| 九九热线精品视视频播放| 久久久久免费精品人妻一区二区| 亚洲成人精品中文字幕电影| 久久99热6这里只有精品| 日本五十路高清| 久久精品国产99精品国产亚洲性色| 国产不卡一卡二| 老师上课跳d突然被开到最大视频| 久久中文看片网| 欧美潮喷喷水| 中文字幕熟女人妻在线| 国产免费男女视频| 久久午夜福利片| 亚洲第一区二区三区不卡| 熟女人妻精品中文字幕| 亚洲精品在线观看二区| 最新在线观看一区二区三区| 搡老岳熟女国产| 在线播放国产精品三级| 天堂动漫精品| 精品久久久噜噜| 国产国拍精品亚洲av在线观看| 一本精品99久久精品77| 又粗又爽又猛毛片免费看| 欧美成人一区二区免费高清观看| 看黄色毛片网站| 国产一区二区三区av在线 | 国产乱人偷精品视频| 国产熟女欧美一区二区| 国产精品三级大全| 亚洲一区高清亚洲精品| 精品福利观看| 毛片一级片免费看久久久久| 高清午夜精品一区二区三区 |