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

    Failure mechanism of rock under ultra-high strain rates

    2021-03-07 12:13:08,,,
    土木與環(huán)境工程學報 2021年1期

    , , ,

    State Key Laboratory of Hydraulic Engineering Simulation and Safety; School of Civil Engineering, Tianjin University, Tianjin 300350, P. R. China)

    Abstract: We explored the failure mechanism of rock under ultra-high strain rates using 3D numerical modeling of the light gas gun test. Based on numerical results, it concluded that mesoscopic hydro-compressive failure rather than mesoscopic shear or tensile failure is the main mechanism of rock failure under the condition of shock compressive loading. The shock wave, indicated by the stress signals of two stress gauges in the rock specimen, can be well reproduced by numerical simulation with the quasi-static rather than the dynamic elastic parameters. The simulation results indicate that the compressive shock wave involves a compressive failure loading process similar to that shown in the conventional uniaxial compressive failure test rather than the ultrasonic test. A mesoscopic-rate-dependent failure model was developed to take the dynamic effect into account. Our results revealed that larger rock porosity could result in an decrease in dynamic strength and dynamic effect under shock compressive loading.

    Keywords: dynamic failure; light-gas gun test; shock compressive loading; lattice spring model

    1 Introduction

    The light-gas gun, first developed for hypervelocity research in aeronautics, was later used in the field of material science. The light-gas gun test is mainly used in material science for: 1) obtaining the Hugoniot curve of the material[1-3], 2) measuring dynamic failure strength[4], 3) studying the high-pressure phase transition[5], and 4) investigating impact-induced chemical reactions[6]. The light-gas gun test has also been applied to rock and rock-like materials. Grote et al.[7]conducted a flat plate impact test on cement mortar and concrete. They reported that the average flow stress of cement mortar and concrete increased to 1.3 GPa and 1.7 GPa, respectively, compared with the unconfined quasi-static compressive strength of 46 MPa and 30 MPa. The dynamic effect of reinforced concrete under shock compressive loading was also explored using the light-gas gun test[8]. The bearing capacity of the reinforced concrete was found to be improved with the increase of impact velocity and reinforcement ratio[8]. Recently, similar phenomena were observed by Zhang et al.[1]with marble and gabbro.

    The dynamic elasticity theory for solid materials was adopted to interpret experimental data of the light-gas gun test[1,7-8]. The constitutive model based on the dynamic elasticity theory can be used with the finite element method (FEM) for numerical simulation of the light-gas gun test. Liu et al.[9]simulated the dynamic fracture mode of tungsten alloy by LS-DYNA (an explicit FEM code) and captured the characteristics in their experimental results. Lopatnikov et al.[10]simulated the light-gas gun test of foamed aluminum plates using LS-DYNA and demonstrated that the dynamic deformation-time relationship and kinetic energy changes of the impact plates at different impact velocities predicted by the FEM were consistent with their theoretical model. Nevertheless, as pointed out by Riedel et al.[11], the intrinsic heterogeneity of rock-like material might result in challenges to the determination of the equation of state parameters and the need for a mesomechanical model. Duan et al.[12]incorporated a statistical isotropic elastic micro-crack model into the FEM to simulate the propagation of plane shock waves in soda-lime glass and reproduced their experimental data. Recently, besides FEM with the mesomechanical model[13-14], a number of discontinuum-based numerical methods[15-20]have also been applied in the study of the dynamic failure of rock. The molecular dynamics (MD)[15]and the lattice type model[18-20]were successfully applied to study the dynamic failure of rock. However, most of these studies only reproduced the failure mode of the rock using a mesoscopic constitutive model with tensile failure. The implementation and practicability of the model to simulate the light-gas gun test is still unclear.

    In this work, we used a lattice type model to explore the light-gas gun test conducted on marble by Zhang et al.[1]. Our main purposes were to answer the following questions:

    1) What kind of mesoscopic damage information can be obtained from the experimental data of the light-gas gun test?

    2) What kind of elastic parameters, the dynamic elastic ones or the quasi-static elastic ones, control the shock loading propagation within the rock specimen?

    3) How to describe the rate dependency of rock under shock compressive loading in the light-gas gun test using a mesoscopic constitutive model?

    4) What is the role of the rock’s mesostructure in the dynamic strength and dynamic effects (strain/loading rate dependency) of rock under shock compressive loading?

    In this work, the distinct lattice spring model (DLSM) was adopted as the numerical tool due to its advantages in modeling rock failure based on simple mechanical elements, e.g., Newton’s second law and springs with simple constitutive models. These characteristics of the DLSM make it a suitable numerical tool for failure mechanism study. In addition, the DLSM is computationally appropriate for full 3D simulation. The following sections of this paper include a brief introduction of the physical test and numerical model, a brief description of the development of rate-dependent constitutive models and a comparison of model response with published data of the light-gas gun test conducted by Zhang et al.[1], based on which the failure mechanism of rock under shock compressive loading is explored. Finally, a few conclusions are derived to answer the aforementioned questions.

    2 Methods

    2.1 The light-gas test

    The light-gas test data reported by Zhang et al.[1], obtained using the experimental facilities at the Cavendish Laboratory, Cambridge, UK[21], are adopted in this work. Fig.1 shows the basic working principle of the experimental facilities as well as the components of the marble specimen and the copper flyer. Different from the traditional rock mechanics tests, two stress gauges were placed in the composite specimen made of marble and PMMA (see Fig.1(b)). The two stress gauges were used to obtain the stress history curve of the corresponding location inside the specimen during the test, that is, the waveform of the shock wave. Since there is a certain distance between the two stress gauges, the velocity of the shock wave in the rock specimen can be obtained based on these signals. After the light-gas gun test, both the composite rock specimen and the flyer became powdery and dissipated, as a result, no morphological data was recorded for the failure process, or the failure pattern of the rock specimen. Another important parameter of the light-gas gun test is the impact velocity of the copper flyer, which is the only controllable input of the light-gas gun test. A number of different impact velocities were adopted by Zhang et al.[1]to obtain the Hugoniot parameters of the marble. Basic material properties of the marble, PMMA and copper flyer were also provided in reference [1]. The longitudinal and shear wave velocities were obtained using an ultrasonic transducer, which can be further used to obtain the dynamic elastic parameters. The quasi-static elastic parameters of the marble were obtained by the standard uniaxial compression test, the parameters of which are listed in Table 1. The dynamic elastic parameters of PMMA and copper were close to their quasi-static counterparts, therefore, no quasi-static tests were conducted on the copper and PMMA.

    2.2 Distinct Lattice Spring Model (DLSM)

    The lattice spring model(LSM) was first developed by Hrennikoff in 1941[22]. Its basic principle is to represent the mechanical responses of a solid through a group of spring-like interactions. The LSM is regarded as the ancestor of both the FEM and the discrete element method (DEM). Due to its simplicity, the LSM has been widely used in the study of many fundamental mechanical phenomena of solids[23-25]. The distinct lattice spring model (DLSM) was developed by Zhao et al.[26]to overcome the Poisson’s limitation in the classical LSM. As shown in Fig.2(a), the basic principle of the DLSM is to represent the solid as a group of particles linked through spring bonds, which consist of a normal spring and a shear spring. Defining the normal unit vector:n=(nxnynz)Tis directed from particleito particlej, which are connected by a normal spring, and the normal deformation of the spring is defined as

    (1)

    whereuij=uj-uiis the relative displacement between particlesjandi. In the DLSM[26], the most commonly used constitutive model for the normal spring is

    (2)

    The multi-body shear spring is one of thedistinctive features of the DLSM. Its principle is to calculate the shear deformation through local strain rather than the displacement of the two particles. In the DLSM, the shear deformation is given as

    (3)

    where [ε]bondis the local strain of the spring bond, which can be obtained from an average operation over the local strain defined at the two particles using a least square method. The main feature of the multi-body shear spring is the ability to represent different Poisson’s ratios without violating the rotation invariance[26]. The failure criterion of the shear spring is given as

    (4)

    Eqs.(2) and (4) are the constitutive models adopted in the DLSM for brittle solids. Fig. 3(a) and (b) show the constitutive model for the normal spring and the constitutive model for the shear spring, respectively. Since a solid is represented as a spring network in the DLSM (as shown in Fig. 2(a)), the fracturing process is manifested by the gradual breakage of these springs (mesoscopic failure events). It has the following advantages in simulating the fracturing of brittle materials. First, these constitutive models are easily implemented, and the simulation results are easily explained, as fewer parameters are required compared to other, discontinuous models, e.g., the bonded DEM. Third, the straight forward parallelization of the DLSM[27-28]gives the model relatively high computational performance.

    Fig.2 Basic principle of the DLSM and the computational model for the light-gas gun

    Fig.3 Classical brittle constitutive models used in the DLSM[26] represented in a dimensionless

    Another distinctive feature of the DLSM is that the relationship between the spring parameters(knandks) and the macroscopic material properties (Eandv) is derived from the hyper-elastic theory. No calibration is required for the mesoscopic elastic parameters, as they can be calculated using the following equations.

    (5)

    (6)

    whereα3Dis a lattice coefficient that can be calculated according to

    (7)

    whereliis the initial length of theithspring andVmis the volume representing the computational model. More details of the DLSM and its recent development can be found in reference [25-27].

    2.3 Computational model for the light-gas gun test

    (8)

    whereAgis the area of the stress gauge and sign(·) is the sign operation, which is given as

    (9)

    Eq. (8) gives the time history of the normal stress on the stress gauge.

    The particle size used in this work is 1 mm and there are about 100 000 particles. The input parameters of the light-gas gun test are based on experimental tests. The output of the numerical test is controlled by the selection of the meso-mechanical constitutive model and the corresponding constitutive parameters. A key point of this work is to minimize the difference between the numerical prediction and the physical experimental result. Thus, the failure mechanism of the rock under shock compressive loading may be interpreted from the searching process of the meso-mechanical constitutive model and its parameters.

    2.4 Tri-linear hydro-compressive model

    Many researchers suggest that macroscopic compression and shear failure is essentially tensile failure at the mesoscopic level[28-30]. However, most of those studies are limited to tensile failure of the normal spring[18, 20, 30]. Fig.3(a) shows the simple brittle constitutive model of the normal spring considering only tensile failure. More comprehensive forms, e.g., considering the damage and plastic deformation, can be developed. Jiang and Zhao[31]recently developed a coupled damage plasticity model to describe the dynamic crack propagation of a gypsum-like 3D printing material. When considering the interaction between two particles, shear failure is a natural logical extension. However, there are very few studies on shear failure in classical LSMs due to the absence of shear interaction. In the DLSM, shear failure can be considered due to the introduction of the multi-body shear spring. The corresponding constitutive model is illustrated in Fig.3(b). Shear failure was widely adopted in the bonded DEM, which usually takes into account both the influence of the cohesion of the contact bond between two particles and the frictional angle. The shear constitutive model presented in Fig.3(b) can be viewed as a special case when the frictional angle is zero. Therefore, the prediction of shear fracturing of the DLSM using the constitutive model shown in Fig.3(b) is not conservative.

    (10)

    Fig.4 A tri-linear constitutive model considering the hydro-compressive failure of the normal spring

    To introduce the effects of plasticity and damage, the same approach as that adopted by Jiang and Zhao[31]is used. First, rewrite Eq.(10) as follows.

    F(u)=S(u)knu

    (11)

    whereS(u) is a state parameter of the normal spring, which is defined as follows.

    (12)

    The interpretation of Eq.(10) can be expressed as the following formula if it is based on damage mechanics.

    F(u)=(1-D(u))knu

    (13)

    Then, the damage function can be expressed as follows.

    D(u)=1-S(u)

    (14)

    In the calculation of the spring force, the maximum damage variable corresponding to the loading history is recorded asD*(u), then the damage constitutive equation considering the loading and unloading can be written as

    F(u)=(1-D*(u))knu

    (15)

    The loading and unloading response of this constitutive model is as shown in Fig.4(a). If Eq.(10) is regarded as a plastic response, it can be expressed as

    F(u)=(u-up)kn

    (16)

    Compare Eq.(14), we will have the plastic deformationupas

    up=u-S(u)u

    (17)

    Similarly, to record the maximum plastic deformationup*experienced by the spring at the moment of loading, the pure plastic constitutive model considering loading and unloading can be written as

    F(u)=(u-up*(u))kn

    (18)

    The corresponding pure plastic constitutive model is shown in Fig.4(b). Combined with Eq.(15) and (18), a coupled damage-plastic model can be constructed by introducing the damage-plastic coupling coefficient as follows.

    F(u)=(1-λdp)(u-up*(u))kn+

    λdp(1-D*(u))knu

    (19)

    Whenλdp=1, it represents the pure damage constitutive model, and whenλdp=0, it represents the pure plastic constitutive model. In this work, we consider the copper behavior as pure plastic material, while the rock has a pure damage response.

    Under dynamic loading, the material has obvious dynamic rate effects. To capture this phenomenon, it is convenient to introduce a rate-dependent model. The approach developed by Zhao et al.[18]is adopted in this work to bring rate dependency into the tri-linear constitutive model. The idea is to change the compressive strength parameters of the trilinear constitutive model as a function of the spring deformation rate. It was found that the instantaneous value of the spring deformation rate may not be able to reproduce the correct macroscopic dynamic effect[18]. To solve this problem, Zhao et al.[18]proposed the concept of time non-localization, which uses the average value of the spring deformation rate from the start time to the current time. In this work, following a similar idea, the local average deformation rate is given as

    (20)

    The average deformation rate of the spring is only calculated when the spring is deformed beyond 0.99uc1. Based on this concept, the strength parameters of the corresponding dynamic tri-linear constitutive model are given by the following formula.

    (21)

    (22)

    Fig.5 Dynamic tri-linear hydro-compressive constitutive

    3 Numerical modeling and discussion

    3.1 Elastic parameters selection

    Fig.6 Shock waveform predicted by the DLSM with

    Fig.7 The shock waveform predicted by the DLSM with different mesoscopic failure

    3.2 Failure mechanism of the marble specimen in the light-gas gun test

    Fig.8 Light-gas gun test predicted by the DLSM considering the mesoscopic tensile

    Fig.9 The shock waveform predicted by the DLSM using the tri-linear constitutive

    Fig.10 Influence of the failure constitutive parameters of the tri-linear hydro-compressive constitutive model on the shock waveform predicted by the

    3.3 Dynamic effect of the mesoscopic hydro-compressive failure

    Fig.11 Dynamic strength of the rock specimen in the light-gas gun test under different impact velocities predicted by the DLSM with the rate independent (RI) constitutive model and rate dependent

    3.4 Influence of mesostructure

    In this section, the influence of the mesostructure on the dynamic failure of rock under shock loading is explored. Fig.12 shows the corresponding computational models for the rock specimens with different mesostructures. These computational models are formed by randomly removing a given number of particles from the marble specimen. The porosity, which was used to characterizs different computational models, is defined asn=number of particles being excavated / total number of particles in the original model. As shown in Fig.12, four pore models with a porosity of 5%, 10%, 15%, and 20% are constructed. All the constitutive parameters are the same as those of the previous section. Fig.13 shows the numerical results of the predicted strength of those different porosity models under different impact velocities. Overall, larger porosity results in lower dynamic strength, which seems logically correct. However, there are some variations, for example, when porosity is 5%, a higher strength is obtained. The reason might be that the inter-particle velocities become more violent under shock compressive loading when a small number of particles were removed from the original model. It might trigger the dynamic effect of the constitutive model and consequently, result in a higher strength. Nevertheless, when the porosity continuously increases, the inter-particle velocity increase may be released in the lateral direction due to the existence of a large mesoscopic free surface. This may result in the increase of strength due to the dynamic effect becoming inconspicuous (as shown in Fig.13).

    Fig.12 Computational models of rock specimens

    Fig.13 Numerical prediction of the dynamic strength of rock specimens with different porosities under various

    4 Conclusions

    In this work,light-gas gun tests on marble were numerically investigated using the DLSM to study the failure mechanism of rock under compressive shock loading. Based on a detailed comparison of the numerical response with published experimental data, it was found that the elastic parameters controlling the compressive shock wave propagation in the rock specimen should be determined through quasi-static tests rather than ultrasonic tests. Moreover, mesoscopic tensile failure and shear failure between rock grains are not the mechanisms of rock failure under compressive shock loading. To reproduce the shock stress waveform observed in the light-gas gun test, a mesoscopic hydro-compressive constitutive model is needed. In other words, it is reasonable to conclude that the light-gas gun test may be able to provide information about the hydro-compressive failure of rock under dynamic loading and provide calibration data for the DLSM or other numerical methods to determine the constitutive equations under dynamic hydro-compression loading conditions. It is noted that a rate-dependent hydro-compressive constitutive model is also necessary to reproduce the experimental observation of the strength increase under different impact velocities. Finally, the influence of the mesostructure in terms of porosity was found to decrease both the dynamic strength and the dynamic effect. The findings in this work may provide a better understanding of rock failure under compressive shock loading, which might be useful for a more rational protective design of underground rock engineering structures under blasting loading[33-36].

    Acknowledgements

    This research is financially supported by the National Key R&D Program of China (Grant No. 2018YFC0406804) and the National Natural Science Foundation of China (Grant No. 11772221).

    国产乱人伦免费视频| 在线视频色国产色| 美女免费视频网站| 久久精品国产综合久久久| 99在线人妻在线中文字幕| 嫩草影视91久久| 男女之事视频高清在线观看| 日韩国内少妇激情av| 一个人免费在线观看的高清视频| 国产日韩一区二区三区精品不卡| 亚洲性夜色夜夜综合| 最近最新中文字幕大全电影3 | 黑人巨大精品欧美一区二区蜜桃| 国产成人一区二区三区免费视频网站| 美女午夜性视频免费| 亚洲五月色婷婷综合| 免费搜索国产男女视频| 纯流量卡能插随身wifi吗| 欧美另类亚洲清纯唯美| 国产主播在线观看一区二区| 50天的宝宝边吃奶边哭怎么回事| 免费高清在线观看日韩| 午夜视频精品福利| 正在播放国产对白刺激| 日韩av在线大香蕉| 国产xxxxx性猛交| 99精品在免费线老司机午夜| 亚洲人成电影观看| 国产又色又爽无遮挡免费看| 日韩精品中文字幕看吧| 中文字幕色久视频| 亚洲欧美一区二区三区黑人| 婷婷精品国产亚洲av在线| 成人欧美大片| 国产麻豆69| 19禁男女啪啪无遮挡网站| 亚洲 欧美 日韩 在线 免费| svipshipincom国产片| 免费观看人在逋| 搡老岳熟女国产| 免费久久久久久久精品成人欧美视频| 亚洲成av片中文字幕在线观看| 男男h啪啪无遮挡| 操出白浆在线播放| 在线观看一区二区三区| 亚洲欧美激情在线| 亚洲成av人片免费观看| 一二三四社区在线视频社区8| 日韩中文字幕欧美一区二区| 老熟妇乱子伦视频在线观看| 很黄的视频免费| 国产真人三级小视频在线观看| 怎么达到女性高潮| 中文字幕久久专区| 99精品欧美一区二区三区四区| ponron亚洲| 亚洲精品一区av在线观看| 91成人精品电影| 亚洲少妇的诱惑av| 欧美日本亚洲视频在线播放| 丁香欧美五月| 精品人妻在线不人妻| 久久狼人影院| 97人妻精品一区二区三区麻豆 | 极品人妻少妇av视频| 男人的好看免费观看在线视频 | 一本大道久久a久久精品| 真人做人爱边吃奶动态| 最新美女视频免费是黄的| 黄片播放在线免费| 三级毛片av免费| 中文字幕人妻熟女乱码| 男人操女人黄网站| 禁无遮挡网站| 美女 人体艺术 gogo| 两人在一起打扑克的视频| 自拍欧美九色日韩亚洲蝌蚪91| 18禁黄网站禁片午夜丰满| 操出白浆在线播放| 大香蕉久久成人网| 久久性视频一级片| 久久久久久亚洲精品国产蜜桃av| 久久这里只有精品19| 中文字幕色久视频| 在线播放国产精品三级| 窝窝影院91人妻| 啦啦啦韩国在线观看视频| 国产精品 欧美亚洲| 校园春色视频在线观看| 久久午夜综合久久蜜桃| 两个人看的免费小视频| 亚洲色图综合在线观看| 不卡av一区二区三区| 亚洲va日本ⅴa欧美va伊人久久| 精品国产国语对白av| 中文字幕精品免费在线观看视频| 黑人操中国人逼视频| 中文字幕精品免费在线观看视频| 超碰成人久久| 天天一区二区日本电影三级 | 日本三级黄在线观看| 免费搜索国产男女视频| 男人操女人黄网站| 国产欧美日韩综合在线一区二区| 很黄的视频免费| 成熟少妇高潮喷水视频| 亚洲情色 制服丝袜| 欧美国产日韩亚洲一区| 99精品在免费线老司机午夜| 99国产精品一区二区三区| 国产精品国产高清国产av| 久久精品91蜜桃| 人人妻,人人澡人人爽秒播| 涩涩av久久男人的天堂| 国产高清激情床上av| 日韩欧美三级三区| 19禁男女啪啪无遮挡网站| 日韩欧美一区二区三区在线观看| 男女下面进入的视频免费午夜 | 免费av毛片视频| 日日摸夜夜添夜夜添小说| 国产aⅴ精品一区二区三区波| 夜夜躁狠狠躁天天躁| 国产1区2区3区精品| 一区二区三区国产精品乱码| 777久久人妻少妇嫩草av网站| 亚洲成人免费电影在线观看| 老熟妇仑乱视频hdxx| 人人妻人人澡欧美一区二区 | 精品熟女少妇八av免费久了| 欧美丝袜亚洲另类 | 亚洲精品美女久久av网站| 亚洲成人国产一区在线观看| 免费高清在线观看日韩| 大型av网站在线播放| cao死你这个sao货| 一区在线观看完整版| 午夜老司机福利片| 99久久综合精品五月天人人| 大香蕉久久成人网| 99热只有精品国产| 国产精品二区激情视频| 国产一区在线观看成人免费| 天堂动漫精品| 麻豆av在线久日| 无人区码免费观看不卡| 高清在线国产一区| 99精品在免费线老司机午夜| 老司机在亚洲福利影院| 色老头精品视频在线观看| 一本综合久久免费| 久久精品亚洲熟妇少妇任你| 亚洲精品在线美女| 国产成人系列免费观看| 欧美黄色片欧美黄色片| 成人国产综合亚洲| 亚洲精品国产一区二区精华液| 最新在线观看一区二区三区| 婷婷丁香在线五月| 亚洲成国产人片在线观看| 性欧美人与动物交配| xxx96com| 日日夜夜操网爽| 国产亚洲精品一区二区www| 欧美精品亚洲一区二区| 每晚都被弄得嗷嗷叫到高潮| 久久人人精品亚洲av| 久久久精品欧美日韩精品| a在线观看视频网站| 国产麻豆69| 黄片播放在线免费| 色婷婷久久久亚洲欧美| 女同久久另类99精品国产91| 女人被躁到高潮嗷嗷叫费观| 国产不卡一卡二| 亚洲熟妇熟女久久| 国产伦人伦偷精品视频| 国产精品 国内视频| tocl精华| 十分钟在线观看高清视频www| 国产成人影院久久av| 亚洲人成伊人成综合网2020| 黄色片一级片一级黄色片| 亚洲av电影不卡..在线观看| 国产精品野战在线观看| 亚洲专区国产一区二区| 色在线成人网| 亚洲成人国产一区在线观看| 日韩欧美三级三区| 狠狠狠狠99中文字幕| 男女之事视频高清在线观看| 一夜夜www| 午夜精品在线福利| 亚洲第一电影网av| 亚洲va日本ⅴa欧美va伊人久久| 老司机靠b影院| 亚洲第一av免费看| 变态另类丝袜制服| 精品免费久久久久久久清纯| 一进一出抽搐动态| 亚洲第一欧美日韩一区二区三区| 成人三级做爰电影| 日日夜夜操网爽| 九色亚洲精品在线播放| 亚洲国产精品sss在线观看| 久久精品人人爽人人爽视色| 久久久久国内视频| 日韩一卡2卡3卡4卡2021年| 国内久久婷婷六月综合欲色啪| 久久 成人 亚洲| 在线观看舔阴道视频| 中文亚洲av片在线观看爽| 精品久久久久久,| or卡值多少钱| 久久久国产成人免费| 精品欧美一区二区三区在线| 一级毛片精品| 亚洲专区国产一区二区| 亚洲国产精品999在线| 亚洲精品粉嫩美女一区| 一边摸一边抽搐一进一小说| 亚洲人成电影免费在线| 在线av久久热| 成人欧美大片| 精品不卡国产一区二区三区| 在线观看免费视频网站a站| 久久久久久国产a免费观看| 99riav亚洲国产免费| 日韩欧美在线二视频| 搡老熟女国产l中国老女人| 97人妻天天添夜夜摸| 精品电影一区二区在线| 1024香蕉在线观看| 日韩欧美免费精品| 久久久久九九精品影院| av片东京热男人的天堂| 欧美一级a爱片免费观看看 | 欧美日韩福利视频一区二区| 97超级碰碰碰精品色视频在线观看| 黄色丝袜av网址大全| 亚洲最大成人中文| 午夜精品国产一区二区电影| 老汉色∧v一级毛片| 91字幕亚洲| 老熟妇仑乱视频hdxx| 18禁观看日本| 啦啦啦免费观看视频1| 免费在线观看完整版高清| 黄频高清免费视频| 久久久久久亚洲精品国产蜜桃av| 久久草成人影院| 精品一区二区三区四区五区乱码| 女人精品久久久久毛片| 老司机午夜十八禁免费视频| 两人在一起打扑克的视频| 国产日韩一区二区三区精品不卡| 亚洲一区二区三区不卡视频| 国产99白浆流出| 男人舔女人下体高潮全视频| 亚洲国产精品合色在线| 亚洲国产欧美日韩在线播放| 亚洲专区中文字幕在线| 久久天堂一区二区三区四区| 国产精品自产拍在线观看55亚洲| 久久久久久久久久久久大奶| 国产成人欧美在线观看| 咕卡用的链子| 亚洲精品国产区一区二| 免费在线观看完整版高清| 亚洲精华国产精华精| 成人三级做爰电影| 亚洲一区二区三区不卡视频| 国产精品影院久久| 十八禁网站免费在线| 十八禁人妻一区二区| 狂野欧美激情性xxxx| 香蕉久久夜色| 日韩精品青青久久久久久| 亚洲人成77777在线视频| 亚洲色图 男人天堂 中文字幕| 成熟少妇高潮喷水视频| 熟妇人妻久久中文字幕3abv| 最好的美女福利视频网| 女性生殖器流出的白浆| 一级,二级,三级黄色视频| 午夜成年电影在线免费观看| 欧美日韩瑟瑟在线播放| 午夜福利18| 欧美日韩乱码在线| 在线av久久热| 美国免费a级毛片| 久久人人97超碰香蕉20202| 亚洲熟妇中文字幕五十中出| 国产欧美日韩一区二区精品| 啦啦啦韩国在线观看视频| 波多野结衣高清无吗| а√天堂www在线а√下载| 精品午夜福利视频在线观看一区| 国产免费男女视频| av欧美777| 精品久久久久久久毛片微露脸| 亚洲欧美一区二区三区黑人| 国产一区二区激情短视频| 国产欧美日韩一区二区精品| 午夜久久久久精精品| 又紧又爽又黄一区二区| 人人妻,人人澡人人爽秒播| 视频区欧美日本亚洲| 制服诱惑二区| 久久影院123| 老汉色∧v一级毛片| 在线天堂中文资源库| 人人妻人人澡人人看| 国产精品电影一区二区三区| 国产av在哪里看| 久久久久久久精品吃奶| 久久午夜综合久久蜜桃| 国产亚洲av高清不卡| 岛国视频午夜一区免费看| 国产精华一区二区三区| 看黄色毛片网站| 中文字幕人成人乱码亚洲影| 日韩高清综合在线| 午夜免费鲁丝| 亚洲精品一区av在线观看| 日韩欧美免费精品| 国产亚洲欧美精品永久| 欧美精品啪啪一区二区三区| 美女免费视频网站| 视频在线观看一区二区三区| 久久国产精品人妻蜜桃| 日本黄色视频三级网站网址| 国产成人精品在线电影| 非洲黑人性xxxx精品又粗又长| 欧美日本视频| 一a级毛片在线观看| 在线观看免费日韩欧美大片| 正在播放国产对白刺激| 成人特级黄色片久久久久久久| 99在线人妻在线中文字幕| e午夜精品久久久久久久| 国产亚洲精品一区二区www| 啦啦啦免费观看视频1| 免费观看人在逋| 黄色丝袜av网址大全| 男人舔女人的私密视频| 国产精品自产拍在线观看55亚洲| 桃色一区二区三区在线观看| 久久久久久久久中文| 中文字幕人成人乱码亚洲影| 69精品国产乱码久久久| 久久久久久久久中文| 久久亚洲精品不卡| 久久久久久久久久久久大奶| 亚洲欧洲精品一区二区精品久久久| 777久久人妻少妇嫩草av网站| 国产精品久久久久久精品电影 | 国产av又大| 丁香六月欧美| 精品久久久久久,| 久久性视频一级片| 精品久久久久久久人妻蜜臀av | 国产一区二区三区在线臀色熟女| 久久天躁狠狠躁夜夜2o2o| 精品第一国产精品| 欧美+亚洲+日韩+国产| 大陆偷拍与自拍| 亚洲成a人片在线一区二区| 亚洲av片天天在线观看| 97碰自拍视频| cao死你这个sao货| 啦啦啦观看免费观看视频高清 | 久久久精品国产亚洲av高清涩受| 精品高清国产在线一区| 美女大奶头视频| 校园春色视频在线观看| 久久精品91蜜桃| 午夜亚洲福利在线播放| 国产av一区在线观看免费| 亚洲欧美激情在线| 18禁国产床啪视频网站| 午夜免费激情av| 免费一级毛片在线播放高清视频 | 亚洲五月天丁香| 一级毛片高清免费大全| 久久狼人影院| tocl精华| 香蕉久久夜色| av在线播放免费不卡| 成人手机av| 久久人妻av系列| 午夜影院日韩av| 久久九九热精品免费| 欧美乱色亚洲激情| 91麻豆精品激情在线观看国产| 亚洲人成电影免费在线| 天天一区二区日本电影三级 | 成熟少妇高潮喷水视频| 一进一出抽搐gif免费好疼| 欧美在线一区亚洲| 中文字幕最新亚洲高清| 自线自在国产av| 中文字幕另类日韩欧美亚洲嫩草| 18禁黄网站禁片午夜丰满| 中文字幕精品免费在线观看视频| 涩涩av久久男人的天堂| 日韩大码丰满熟妇| 精品高清国产在线一区| 一a级毛片在线观看| 18禁黄网站禁片午夜丰满| 午夜激情av网站| 国产欧美日韩精品亚洲av| 女人被狂操c到高潮| 久久国产精品影院| 日韩视频一区二区在线观看| 美女扒开内裤让男人捅视频| 免费在线观看影片大全网站| 51午夜福利影视在线观看| 人人澡人人妻人| 国产精品影院久久| 两性夫妻黄色片| 看黄色毛片网站| 国产又色又爽无遮挡免费看| 日韩一卡2卡3卡4卡2021年| 女性被躁到高潮视频| 黄色视频不卡| 首页视频小说图片口味搜索| 不卡av一区二区三区| 搡老岳熟女国产| 国产熟女午夜一区二区三区| netflix在线观看网站| 国产野战对白在线观看| 成年人黄色毛片网站| www日本在线高清视频| 午夜福利,免费看| 性少妇av在线| 久热爱精品视频在线9| 啦啦啦 在线观看视频| 老司机福利观看| 少妇粗大呻吟视频| cao死你这个sao货| 亚洲色图综合在线观看| 精品国产美女av久久久久小说| 桃色一区二区三区在线观看| 女同久久另类99精品国产91| 午夜亚洲福利在线播放| 亚洲情色 制服丝袜| 伊人久久大香线蕉亚洲五| 日日干狠狠操夜夜爽| 少妇的丰满在线观看| 国产精品亚洲美女久久久| 久久久久久国产a免费观看| 桃色一区二区三区在线观看| 国产精品国产高清国产av| av在线天堂中文字幕| 窝窝影院91人妻| 亚洲精品久久成人aⅴ小说| 美女 人体艺术 gogo| 国产一级毛片七仙女欲春2 | 他把我摸到了高潮在线观看| 在线观看66精品国产| 国产精品亚洲一级av第二区| 国产精品av久久久久免费| 日本在线视频免费播放| 婷婷精品国产亚洲av在线| 99在线视频只有这里精品首页| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲精品久久成人aⅴ小说| 国产极品粉嫩免费观看在线| 日韩精品中文字幕看吧| 欧美日韩中文字幕国产精品一区二区三区 | √禁漫天堂资源中文www| 一区二区三区高清视频在线| 一级片免费观看大全| 久久 成人 亚洲| 一进一出好大好爽视频| 国产单亲对白刺激| 国产亚洲精品第一综合不卡| 在线观看免费日韩欧美大片| 一边摸一边抽搐一进一出视频| 亚洲精品美女久久久久99蜜臀| 亚洲欧美激情综合另类| 日韩中文字幕欧美一区二区| 国产一区二区三区在线臀色熟女| 黄色丝袜av网址大全| 巨乳人妻的诱惑在线观看| 午夜精品久久久久久毛片777| 日韩精品中文字幕看吧| 51午夜福利影视在线观看| 男男h啪啪无遮挡| 变态另类丝袜制服| xxx96com| 亚洲激情在线av| 丰满人妻熟妇乱又伦精品不卡| 韩国精品一区二区三区| 一个人观看的视频www高清免费观看 | 亚洲av第一区精品v没综合| 久久精品国产综合久久久| 一二三四社区在线视频社区8| 极品人妻少妇av视频| 丝袜在线中文字幕| 一级片免费观看大全| 国产亚洲精品av在线| 国产一区二区三区综合在线观看| 精品日产1卡2卡| 人人妻人人澡欧美一区二区 | 身体一侧抽搐| 色精品久久人妻99蜜桃| 免费在线观看视频国产中文字幕亚洲| 美女高潮到喷水免费观看| 9色porny在线观看| 午夜a级毛片| 国产激情久久老熟女| 亚洲欧美日韩无卡精品| 女人被狂操c到高潮| 波多野结衣高清无吗| 国内精品久久久久久久电影| 国产欧美日韩一区二区三区在线| cao死你这个sao货| 脱女人内裤的视频| 亚洲,欧美精品.| 色av中文字幕| 国产精品香港三级国产av潘金莲| 三级毛片av免费| 午夜福利高清视频| 一本久久中文字幕| 日本黄色视频三级网站网址| 嫁个100分男人电影在线观看| 亚洲午夜理论影院| 视频区欧美日本亚洲| 9色porny在线观看| 1024香蕉在线观看| 99久久综合精品五月天人人| 亚洲国产精品sss在线观看| 国产精品亚洲av一区麻豆| 国产高清视频在线播放一区| 禁无遮挡网站| 亚洲专区国产一区二区| 亚洲狠狠婷婷综合久久图片| 一级毛片高清免费大全| 久久久国产精品麻豆| 一卡2卡三卡四卡精品乱码亚洲| 日本免费一区二区三区高清不卡 | 少妇 在线观看| 久久天躁狠狠躁夜夜2o2o| 国内精品久久久久精免费| 亚洲男人天堂网一区| 国产片内射在线| 国产蜜桃级精品一区二区三区| 国产亚洲欧美98| 午夜福利欧美成人| 高清在线国产一区| 亚洲欧美日韩无卡精品| 99精品久久久久人妻精品| 欧美乱码精品一区二区三区| videosex国产| 中文字幕av电影在线播放| 国产精品一区二区三区四区久久 | 国产麻豆成人av免费视频| 999久久久国产精品视频| 亚洲色图 男人天堂 中文字幕| 成人国产综合亚洲| 日韩av在线大香蕉| 久久精品91蜜桃| 日本黄色视频三级网站网址| 黄色毛片三级朝国网站| 精品久久久久久成人av| 又紧又爽又黄一区二区| 最新美女视频免费是黄的| 国产精品影院久久| 日韩视频一区二区在线观看| 中文字幕色久视频| 久久精品亚洲熟妇少妇任你| 男女做爰动态图高潮gif福利片 | 色精品久久人妻99蜜桃| 久久青草综合色| 色在线成人网| 成人手机av| 美女国产高潮福利片在线看| videosex国产| 正在播放国产对白刺激| 欧美一区二区精品小视频在线| 亚洲七黄色美女视频| 久久 成人 亚洲| 国产欧美日韩一区二区三| 日韩一卡2卡3卡4卡2021年| 免费无遮挡裸体视频| 亚洲熟妇中文字幕五十中出| 免费久久久久久久精品成人欧美视频| 欧美成人免费av一区二区三区| 亚洲欧美激情在线| 亚洲人成77777在线视频| 亚洲天堂国产精品一区在线| 不卡一级毛片| 日韩av在线大香蕉| 欧美国产精品va在线观看不卡| 国产精品永久免费网站| 中文亚洲av片在线观看爽| 人成视频在线观看免费观看| 啦啦啦韩国在线观看视频| 久久天堂一区二区三区四区| 亚洲免费av在线视频| 免费在线观看影片大全网站| 色播亚洲综合网| 精品久久久精品久久久| 变态另类丝袜制服| 亚洲在线自拍视频| 色综合亚洲欧美另类图片| www.精华液| 日本免费一区二区三区高清不卡 | 色综合欧美亚洲国产小说| 国产欧美日韩一区二区三区在线| 国产精品久久视频播放| 麻豆一二三区av精品| 97人妻精品一区二区三区麻豆 |