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

    Investigation on experimental method of low-impedance materials using modified Hopkinson pressure bar

    2015-04-22 06:17:34MIAOYinggang苗應(yīng)剛LIYulong李玉龍DENGQiong鄧瓊TANGZhongbin湯忠斌HUHaitao胡海濤SUOTao索濤
    關(guān)鍵詞:玉龍海濤

    MIAO Ying-gang(苗應(yīng)剛), LI Yu-long(李玉龍), DENG Qiong(鄧瓊),TANG Zhong-bin(湯忠斌), HU Hai-tao(胡海濤), SUO Tao(索濤)

    (School of Aeronautics, Northwestern Polytechnical University, Xi’an 710072, China)

    ?

    Investigation on experimental method of low-impedance materials using modified Hopkinson pressure bar

    MIAO Ying-gang(苗應(yīng)剛), LI Yu-long(李玉龍), DENG Qiong(鄧瓊),TANG Zhong-bin(湯忠斌), HU Hai-tao(胡海濤), SUO Tao(索濤)

    (School of Aeronautics, Northwestern Polytechnical University, Xi’an 710072, China)

    To increase the detectability of split Hopkinson pressure bar (SHPB) of low-impedance materials, modifications were conducted on traditional SHPB apparatus with a PMMA tube to output transmitted signal, and weak signals were further amplified by semiconductor strain gauges . Experiments on soft rubbers and cushioning foam materials were carried out. In order to analyze the accuracy of the experimental results, the stress equilibrium issues involved in the assumptions of SHPB were investigated. First, by way of re-constructing loading process of incident wave, the stress-strain curve was obtained, along with the stress equilibrium ratio of specimen. Secondly, the influences on the accuracy of stress-strain curves were investigated through the elastic modulus comparisons. And the results illustrate that the bilinear incident wave from experiments can ensure the stress equilibrium deformation of specimen after 2 normalized times, much sooner than ramp incident waves. Moreover, it even facilitates specimen deformation with a constant strain rate. The results confirm that the detectability of the modified SHPB can be down to tens kPa with enough accuracy level.

    SHPB; low-impedance materials; bilinear wave; stress equilibrium; constant strain rate

    Low mechanical impedance materials have gotten more and more popularizations for specific applications. For example, kinds of foam material and sandwiched composites have been widely used in automatic, even aerospace industries for their high specific weight and strength; soft rubbers are employed in portable electronic devices for absorbing impact energy. Thus, for higher level of security and severability of the related devices/structure, it is necessary to model and analyze their mechanical behaviors of materials involved in situations such as impacting. It is indispensable to require a full understanding of their mechanical responses, especially under dynamic loading. The traditional loading technique, such as the split Hopkinson bar, is widely used to measure the dynamic behavior of materials under impact loading. It was originally designed for measuring transient explosive wave by Hopkinson[1], and revolutionarily modified by Davies[2], and Kolsky[3]to test dynamic behavior of materials under high strain rate.

    In the last decades, great progresses had been made in the Hopkinson bar techniques on low impedance materials. Actually, it is very difficult to conduct experiments on these kinds of materials using split Hopkinson pressure bar (SHPB), because their mechanical strength is usually much lower than that of metal specimen, at least by one or two orders, leading to very weak transmitted signals that can hardly be detected. In order to amplify the weak signals, related techniques were utilized to enhance the transmitted signal, which used an aluminum tube to output the transmitted signals as well as higher-sensitivity strain transducers such as the quartz-crystal sensor or semiconductor strain gauges for much low-impedance materials to measure the weak transmitted signals[4-6]. Besides, viscous-elastic Hopkinson bars were also established to measure the low-impedance materials with the aid of analytical solutions to wave dispersion in viscous-elastic bar[7-9]. However, data interpretation involved with complicated solution of wave dispersion in viscous-elastic Hopkinson bar, which confronted its applications.

    An alternative difficulty was to design the incident wave, for ensuring the accuracy of stress equilibrium in the specimen. For low-impendence materials, achieving equilibrium deformation as soon as possible is essential for the accuracy level of experimental results. Ravichandran[10]concluded that the specimen could achieve a stress equilibrium when the ramp loading incident wave reverberated in the specimen for 3 times. But the duration still would be too long, for the reason that the wave speed of soft materials was lower than that of normal engineering materials by one or two orders. After analyzing the effect of the rising edge of incident wave, Yang[11]found that the time duration to achieve specimen uniformity of stress can be shortened if the incident wave had a specific profile. However this kind of wave has not been reported experimentally so far.

    Fig.1 Modified split Hopkinson bar apparatus for low-impedance materials

    Moreover, inertial radial effect under dynamic loading is essential to distinguish from quasi-static loading. In the past decades, Kolsky[3], Davies[12], and Gorham[13]investigated the specimen stress after well considering inertial effects and found that in a constant strain rate it can reduce greatly the inertial effect. Warren[14]commented on the radial inertial effect and found that, for soft materials such as foams, rubber with low flow stress, the radial inertial should be considered. Only in the case of specimen deformation under constant strain rate loading, the inertial term could be reduced, even eliminated[3, 12-13].

    This paper focused mainly on issues involved in measuring soft materials using SHPB techniques.We proposed a modified SHPB apparatus, with the metal-incident bar, and the PMMA tube as transmitted bar. Tests were conducted with the aid of the pulse shaping technique for stress equilibrium and constant strain rate. Meanwhile, stress equilibrium issues were analyzed elaborately and the corresponding numerical calculations were carried out to ascertain the influence of stress equilibrium on the accuracy of experimental results.

    1 Modified SHPB and Experiments

    1.1 Modified SHPB

    In order to amplify the transmitted signals, in the present work we proposed a modified Hopkinson bar configuration as shown in Fig.1. A PMMA transmitted tube as well as high sensitivity-factor semiconductor strain gauge (gauge factor up to 150) was adopted to further improve the detectability of transmitted signals. The semiconductor strain gauges were cemented closely to the interface of specimen/tube, with 50 mm far away. A PMMA pad with size ofφ15 mm×3 mm was fixed on the ends of the transmitted tube, to support the sandwiched specimen deformation. The traditional aluminum incident bar was still used.

    Fig.2 shows the sandwiched section, where the contacting interfaces are labeled as 1 and 2. The stresses and displacements in interfaces 1 and 2 were deduced as follows.

    Fig.2 Sandwiched section of split Hopkinson pressure bar

    The stress, strain and strain rate of the specimen were deduced based on one-dimension stress wave theory[15].

    (1)

    whereεi(t),εr(t) andεt(t) were the incident wave, the reflected wave, and the transmitted wave respectively;CBandCTwere the wave speed of incident bar, and the PMMA tube, respectively; andls. was the specimen length. Thus the strain of specimen could be calculated by

    (2)

    Thentheaveragestressinthespecimen

    couldbeexpressed:

    (3)

    whereAB,AT,Aswere the area of the incident bar, the transmitted tube and the specimen, respectively. AndEB,ETwere the elastic modulus of the incident bar and the transmitted tube, respectively. Then the equations were obtained below:

    (4)

    1.2Experimentdetails

    Tab.1 Dimensions and mechanic parameters of SHPB and specimen

    The specific impedance ratios were the relative impedance values of related bars to the rubber specimen. Specific pulse shapers, superimposed by multilayers tissue and plasticine, were utilized to elongate the rising time of incident wave[16-17]. The incident wave could be trimmed to achieve wave duration about 400 μs with the rising time of about 90 μs, as shown in Fig.3. The corresponding reflected and transmitted waves are also shown in Fig.3.

    Fig.3 Strain signals from modified SHPB

    It can be observed that the transmitted strain signal had been amplified to be in the same order with the incident wave and reflected wave. Similar tests were also conducted using a traditional Hopkinson pressure bar, which used an aluminum bar as the transmitted bar, and the traditional strain gauge to record strain signals. The transmitted strain signal from the traditional aluminum bar was displayed in Fig.3 too. It was hard to distinguish the transmitted strain signal due to its strong weakness. Based on the amplitude of PMMA tube strain from semiconductor strain gauge, it can be calculated that the strain in the aluminum transmitted bar was 7 μεaround. It was in the same order as the noise signal, even covered by the noise.

    Besides, in order to checkthe potential detectability of this SHPB, a lower impedance material, a cushioning foam was selected as specimen, and experiments were performed with the size of 8.6 mm×9.1 mm (as shown in Fig.4), and 2.5 mm in thickness. The loading condition was the same as that of rubber specimen experiments. Namely, the incident wave was also trimmed to be the same as that of Fig.3.

    Fig.4 Cushioning foam specimen

    2 Results and discussions

    2.1 Experimental results

    The true stress-strain curve of rubber specimen and corresponding strain rate were shown in Fig.5. The fitted slope of stress-strain curve was nearly 4 MPa withinεt=0.10,R2=0.962 6 (whereRis correlation coefficient in the least squares method for linear regression). The deformation of specimen had come into nearly a constant strain rate sinceεt=0.10, and the strain rate was around 1 350 s-1.

    Fig.5 True stress-strain curve and accompanying strain rate of soft rubber

    Fig.6 showed the stress-strain curve and accompanying strain rate of the cushioning foam specimen. It came into compacting with the increasing loading. The engineering stress was given for quantifying its mechanical behaviors. It was found that the stress of specimen was lower than 0.06 MPa untilεt=0.45. The fitted slope of stress-strain curve was 0.136 3 MPa withinεt=0.15. The deformation of specimen had come into an approximate constant strain rate sinceεt=0.15, and the strain rate was around 1 880 s-1. Thus it confirmed the detectability of this modified split Hopkinson bar apparatus, which might can test materials with strength down to tens kPa.

    Fig.6 True stress-strain curve and strain rate of the cushioning foam specimen

    In the aforementioned results, although the stress-strain curves and constant strain rates were obtained by utilizing the proposed apparatus, whether the results were the real performances of materials or not are still unknown. Thus the assumptions of the proposed SHPB were needed to be analyzed, especially in the stress equilibrium issue.

    2.2 Formula derivation for analyzing stress equilibrium

    In order to analyze loading procedures of stress wave in specimen, the numerical calculation to re-construct the actual loading process was conducted by computing the stresses and displacements in interfaces 1 and 2, where the reflecting and transmitting of the loading stress wave occurred.

    Fig.7 elaborately shows the actual process of waves’ reflecting and transmitting at both interfaces. Here we defined the normalized timeτ=ls/cs, which was the time needed for the stress wave to propagate between interfaces, wherelsandcsare the length and wave speed of specimen respectively. The reflected coefficientsF, and the transmitted coefficientsTare calculated based on one-dimensional stress wave theory[18]and are listed in the below:

    (5)

    whereTBSmeans the transmitted coefficient of stress wave from bar to specimen, and then the label strategy could be applied to other coefficients. The stress at the interfaces between bars and specimen could be deduced as following:

    (6)

    (7)

    Fig.7 Waves’ reflecting and transmitting at interfaces of specimen and bars

    Here,aparameterR(t) was defined to evaluate the state of stress disequilibrium in specimen as below[10]. It was universally considered to be in a stress equilibrium, whenR(t) were lower than 5%[10]. In the following sections, the stress disequilibrium factor will be used to determine whether the stress equilibrium is achieved in the specimen.

    (8)

    2.3Numericalcalculation

    Thenumericalcalculationswereconductedontheabove-mentionedsoftrubberspecimen.Theincidentwavepresentedabinary-linearconfiguration,whichwasabstractedfromtheincidentwaveofFig.3.TheparametersofHopkinsonbarandspecimenweredeterminedfromTab.1.Thespecimenwasendowedtobelinearelasticperformanceforsimplifyingthecalculation.Thenormalizedtimewhichisabout45 μs was calculated based on its elastic modulus, density and its length. The loading details were plotted in Fig.8 together with related stress equilibrium ratios. This specific loading with the turning point at 2 normalized times facilitated the stress equilibrium ratio coming down to values below 5% quickly. That is, the rubber specimen could come into stress equilibrium deformation quickly just after 2τ.

    Fig.8 Abstracted binary-linear incident wave and its stress equlibrium ratios

    In order to highlight the advantage of bilinear incident wave, additional calculations were also carried out for the PMMA transmitted tubes, PMMA transmitted bar, and aluminum bar, under ramp incident wave loading. The stress equilibrium shared the same tends with the results of the reference[10]. Within normalized time 2τ, the ratios were fully larger than 5%. After that, they were mostly larger than 5% before 5τ. Since then, the specimen could arrive in a stress equilibrium stage. For the involved transmitted bars or tube, there were few differences with the stress equilibrium ratio (as shown in Fig.9), even the specific impedance ratios mismatch dramatically from 515.5 of Al bar to 1 of specimen, as listed in Tab.1. After comparing the ramp wave with bilinear wave, and their influence on the stress equilibrium ratios, it can be concluded that under the bilinear wave loading, the lower slope of the second-linear loading after 2τplays an important role in meeting the requirement of stress equilibrium ratio coming down to values below 5%.

    Fig.9 Ramp loading incident wave’s configuration, and stress equilibrium ratios

    In a word, the rubber specimen could achieve stress equilibrium deformation with this nearly bilinear loading wave. This specific loading wave could facilitate specimen stress equilibrium deformation very sooner.

    2.4 Influence of the stress equilibrium

    The specimen underwent severe deformation of stress disequilibrium state at the very beginning of deformation. The influence of disequilibrium deformation on the stress-strain curves should be worthy of knowing, for it may affect the accuracy of stress strain curve, especially the elastic modulus determination of linear behavior of materials. Thus, based on the reconstructing procedures, the stress-strain curves were re-constructed by Eqs.(4)(6). The parameters and size of bar and specimen were shown in Tab.1. The loading wave was the binary-linear wave from experiments of rubber materials, namely same with the incident strain wave of Fig.8.

    The fitted elastic modulus was obtained by fitting the reconstructed stress-strain curve, which were displayed in Fig.10. For the input elastic modulus 4MPa, the 5% error zone was between 3.8 MPa to 4.2 MPa. Within the loading duration between 2τto 5τ, the modulus were in the 5% error zone. Even with 2τwhen specimens had just undergone severe disequilibrium deformation, the fitted modulus were 4.157 MPa, by only about 4% error more than the input modulus 4 MPa. So with this bilinear incident wave, Eq.(4) could interpret the stress, strain as well as strain rate history with enough accuracy.

    Fig.10 Trends of fitted elastic moduli of rubber and stress equilibrium ratio with normalized time

    Fig.11 Trend of fitted elastic modulus of cushioning foam specimen versus normalized times

    3 Remarking and conclusions

    In this paper, a PMMA tube was selected as a transmitted bar to measure the weak transmitted signal due to low impedance of tested materials. Meanwhile, the semiconductor strain gauge with its gauge factor much higher than the traditional strain gauge by 75 times was adopted to further amplify the transmitted signals. Due to the low-impedance of specimen, the normalized time of specimen would be much longer than that of traditional metal. Thus it would take a long time to achieve the stress equilibrium deformation of specimen. Fortunately, the bilinear incident wave can facilitate the specimen deformation to come into stress equilibrium state after 2τ. For the low-impedance specimen with linear elastic behavior, the stress disequilibrium influence on the stress-strain curve was negligible, with error in 5% from, in the view of elastic modulus determination accuracy. That is, the bilinear loading wave can be adopted in SHPB experiments of the low-impedance materials, with the advantages of both achieving a stress equilibrium sooner with constant strain rate deformation. Therefore, the proposed apparatus could be used to measure the dynamical behavior of some low impedance materials with detectability down to tens kPa.

    [1] Hopkinson B. A method of measuring the pressure in the deformation of high explosives or by the impact of bullets [J]. Philosophical Transactions of the Royal Society of London,1914,A213:437-452.

    [2] Davies R M. A critical study of Hopkinson pressure bar[J]. Philosophical Transactions of the Royal Society of London,1948, A240:375-457.

    [3] Kolsky H. An investigation of the mechanical properties of materials at very high rates of loading[J]. Proceedings of the Royal Society ofLondon, 1949,B62:676-700.

    [4] Chen W, Zhang B, Forrestal M J. A split Hopkinson bar technique for low-impedance materials[J]. Experimental Mechanics, 1999,39(2):81-85.

    [5] Chen W, Lu F, Zhou B. A quartz crystal imbedded split Hopkinson bar for softmaterials[J]. Experimental Mechanics, 2000,40:1-6.

    [6] Trautmann A, Siviour C R, Walley S M, et al. Lubrication of polycarbonate at cryogenic temperatures in the split Hopkinson pressure bar [J]. International Journal Impact Engineering, 2005,31:523-544.

    [7] Zhao H, Gary G, Klepaczko J R. On the use of a viscoelastic split Hopkinson pressure bar[J]. International Journal Impact Engineering, 1997,19:319-330.

    [8] Wang L, Labibes K, Azari Z, et al. Generalization of split Hopkinson bar technique to use viscoelastic bars[J]. International Journal Impact Engineering, 1994,15:669-686.

    [9] Zhao H, Gary G. A three dimensional analytical solution of the longitudinal wave propagation in an infinite linear viscoelastic cylindrical bar: application to experimental techniques[J]. Journal of the Mechanics and Physics of Solids, 1995,43(8):1335-1348.

    [10] Ravichandran G, Subhash G. Critical appraisal of limiting strain rate for compression testing of ceramics in a split Hopkinson pressure bar[J]. Journal of America Ceramic Society, 1994,77(1):263-267.

    [11] Yang L M, Shim V P M. An analysis of stress equilibrium uniformity in split Hopkinson bar test specimens[J]. International Journal Impact Engineering, 2005,28:129-150.

    [12] Davies E D H, Hunter S C. The dynamic compression testing of solids by the method of the split Hopkinson pressure bar[J]. J Mech Phys Solids, 1963,11:155-179.

    [13] Gorham D A. Specimen inertia in high strain-rate compression [J]. J Phys D Appl Phys, 1989,22:1888-1893.

    [14] Warren T L, Forrestal M J. Comments on the effect of radial inertial in the Kolsky bar test for an incompressible material[J]. Experimental Mechanics, 2010,50:1253-1255.

    [15] Lindholm U S. Some experiments with split Hopkinson pressure bar [J]. Journal of the Mechanics and Physics of Solids, 1964,12(3):317-335.

    [16] Frew D J, Forrestal M J, Chen W. Pulse shaping techniques of testing brittle materials with a split Hopkinson pressure bar [J]. Experimental Mechanics, 2002,42(1):93-106.

    [17] Yuan Q L, Li Y L, Li H J, et al. Quasi-static and dynamic compressive fracture behavior of carbon/carbon composites [J]. Carbon, 2008,46:699-703.

    [18] Wang L L. Foundations of stress waves[M]. Amsterdam: Elsevier,2007: 69-70.

    (Edited by Wang Yuxia)

    10.15918/j.jbit1004- 0579.201524.0220

    O 347.3 Document code: A Article ID: 1004- 0579(2015)02- 0269- 08

    Received 2014- 04- 29

    Supported by the National Natural Science Foundation (11272267, 11102168, 10932008);111 Project ( B07050)

    E-mail: suotao@nwpu.edu.cn

    猜你喜歡
    玉龍海濤
    紅山玉龍
    ——中華第一玉龍
    家教世界(2023年28期)2023-11-14 10:08:18
    紅山玉龍
    羅海濤作品
    國畫家(2022年3期)2022-06-16 05:30:06
    圓圓的世界
    感受肌理
    通過反思尋求最優(yōu)解
    玉龍喀什水利樞紐施工
    紅山文化“玉龍”中的“猿”
    大眾考古(2015年6期)2015-06-26 08:27:12
    玉龍山
    海外英語(2013年7期)2013-11-22 08:25:45
    Stochastic simulation of fluid flow in porous media by the complex variable expression method*
    一区二区三区四区激情视频 | 国产日韩欧美在线精品| 只有这里有精品99| 日本撒尿小便嘘嘘汇集6| 一进一出抽搐动态| 国产探花在线观看一区二区| 人妻夜夜爽99麻豆av| 国产精品免费一区二区三区在线| 亚洲aⅴ乱码一区二区在线播放| 成人性生交大片免费视频hd| av专区在线播放| 久久精品91蜜桃| 26uuu在线亚洲综合色| 99国产极品粉嫩在线观看| 国产成人福利小说| 狂野欧美白嫩少妇大欣赏| 日韩视频在线欧美| 国产伦精品一区二区三区四那| 亚洲国产精品合色在线| 18禁黄网站禁片免费观看直播| 床上黄色一级片| 99久久精品国产国产毛片| 欧美成人a在线观看| 麻豆精品久久久久久蜜桃| 精品无人区乱码1区二区| 国产单亲对白刺激| 成人永久免费在线观看视频| 国产一区二区在线观看日韩| av在线观看视频网站免费| 国产爱豆传媒在线观看| 中文字幕熟女人妻在线| 国产单亲对白刺激| 三级经典国产精品| 高清午夜精品一区二区三区 | 国产色爽女视频免费观看| av福利片在线观看| 国产一区亚洲一区在线观看| 国产精品人妻久久久影院| 91麻豆精品激情在线观看国产| 日日摸夜夜添夜夜爱| 国产高清有码在线观看视频| 高清毛片免费看| 成熟少妇高潮喷水视频| 国产高清不卡午夜福利| 亚洲一区二区三区色噜噜| 亚州av有码| 国产精品久久久久久精品电影| 亚洲国产精品成人综合色| 久久99热6这里只有精品| 国产中年淑女户外野战色| 国产伦精品一区二区三区四那| 欧美一区二区精品小视频在线| 国产成人精品婷婷| a级毛片a级免费在线| 可以在线观看的亚洲视频| 精品国内亚洲2022精品成人| 男人舔奶头视频| 丰满乱子伦码专区| 男人狂女人下面高潮的视频| 蜜桃久久精品国产亚洲av| 少妇高潮的动态图| 亚洲av.av天堂| 久久这里只有精品中国| 毛片一级片免费看久久久久| 一级毛片我不卡| 18禁在线无遮挡免费观看视频| 日韩av在线大香蕉| 最近最新中文字幕大全电影3| 在线a可以看的网站| 免费av毛片视频| 欧美成人一区二区免费高清观看| 国产成人精品一,二区 | 亚洲欧美精品综合久久99| 亚洲,欧美,日韩| 亚洲无线观看免费| 一区福利在线观看| 午夜免费激情av| 久久人妻av系列| 免费观看在线日韩| 一区二区三区高清视频在线| a级毛片免费高清观看在线播放| 亚洲av成人av| 一级黄片播放器| 亚洲国产日韩欧美精品在线观看| 久久精品国产亚洲av涩爱 | 一本精品99久久精品77| 久久精品91蜜桃| 久久人人爽人人爽人人片va| 欧美性猛交黑人性爽| 少妇熟女欧美另类| 午夜精品国产一区二区电影 | 久久久久久国产a免费观看| 国产精品美女特级片免费视频播放器| 99久久精品一区二区三区| 男女啪啪激烈高潮av片| 久久99蜜桃精品久久| 国产精品福利在线免费观看| 看非洲黑人一级黄片| 亚洲欧美成人综合另类久久久 | 干丝袜人妻中文字幕| www日本黄色视频网| 色吧在线观看| 99国产极品粉嫩在线观看| 欧美+亚洲+日韩+国产| 欧美又色又爽又黄视频| 欧美bdsm另类| 一区福利在线观看| 久久久久久伊人网av| 又黄又爽又刺激的免费视频.| 一级毛片久久久久久久久女| 精品熟女少妇av免费看| 又黄又爽又刺激的免费视频.| av视频在线观看入口| 1024手机看黄色片| 中文字幕精品亚洲无线码一区| 亚洲在线自拍视频| 亚洲精品日韩在线中文字幕 | 九草在线视频观看| 国产老妇女一区| 又粗又硬又长又爽又黄的视频 | 欧美日韩综合久久久久久| 三级男女做爰猛烈吃奶摸视频| 黑人高潮一二区| 欧美性猛交黑人性爽| 久久精品国产99精品国产亚洲性色| 久久鲁丝午夜福利片| 老熟妇乱子伦视频在线观看| 亚洲欧洲日产国产| 亚洲欧美日韩卡通动漫| 联通29元200g的流量卡| 亚洲最大成人中文| 亚洲熟妇中文字幕五十中出| 女人被狂操c到高潮| 大香蕉久久网| 天堂中文最新版在线下载 | 色综合站精品国产| 免费黄网站久久成人精品| 国产精品伦人一区二区| 一本精品99久久精品77| 免费看a级黄色片| 久久6这里有精品| 天堂影院成人在线观看| 久久99精品国语久久久| 99久久中文字幕三级久久日本| 国产国拍精品亚洲av在线观看| 亚洲无线在线观看| 日本三级黄在线观看| 成年女人看的毛片在线观看| 青青草视频在线视频观看| 91久久精品国产一区二区成人| 亚洲精品久久国产高清桃花| 在线播放无遮挡| 亚洲成人中文字幕在线播放| 免费观看的影片在线观看| 亚洲av中文字字幕乱码综合| 精品国产三级普通话版| 精品欧美国产一区二区三| 97超碰精品成人国产| 亚洲欧美日韩东京热| 丰满的人妻完整版| 久久久久性生活片| 天堂中文最新版在线下载 | 99热这里只有是精品50| 午夜福利在线在线| 一本久久中文字幕| 成人永久免费在线观看视频| 中文精品一卡2卡3卡4更新| av在线天堂中文字幕| 亚洲久久久久久中文字幕| 夜夜看夜夜爽夜夜摸| 久久久久久国产a免费观看| 晚上一个人看的免费电影| 日韩欧美国产在线观看| 亚洲精品日韩在线中文字幕 | 麻豆久久精品国产亚洲av| 国产精品精品国产色婷婷| 欧美高清性xxxxhd video| 亚洲人成网站高清观看| 91av网一区二区| 日韩高清综合在线| 特级一级黄色大片| 一边亲一边摸免费视频| 男女视频在线观看网站免费| 亚洲av中文字字幕乱码综合| av免费在线看不卡| 欧美日韩乱码在线| 国内精品久久久久精免费| 国产三级在线视频| 麻豆国产av国片精品| 亚洲在久久综合| 欧美日韩综合久久久久久| 亚洲激情五月婷婷啪啪| 久久精品国产清高在天天线| 狂野欧美白嫩少妇大欣赏| 国产毛片a区久久久久| 免费不卡的大黄色大毛片视频在线观看 | 久久亚洲国产成人精品v| 欧美最新免费一区二区三区| 亚洲乱码一区二区免费版| 18禁在线无遮挡免费观看视频| 可以在线观看的亚洲视频| 大香蕉久久网| 人人妻人人澡欧美一区二区| 好男人在线观看高清免费视频| 国产午夜精品论理片| 免费av毛片视频| 久久久久性生活片| 国产女主播在线喷水免费视频网站 | 亚洲一区二区三区色噜噜| 亚洲av一区综合| 国产精品永久免费网站| 欧美又色又爽又黄视频| 中国美女看黄片| 少妇的逼水好多| 热99在线观看视频| 亚洲一区高清亚洲精品| 国内精品美女久久久久久| 一级毛片久久久久久久久女| 一进一出抽搐动态| 成人亚洲欧美一区二区av| 亚洲精品亚洲一区二区| 国产精华一区二区三区| 欧美激情国产日韩精品一区| 91午夜精品亚洲一区二区三区| 男人舔女人下体高潮全视频| 51国产日韩欧美| 黄色日韩在线| 岛国在线免费视频观看| kizo精华| 午夜福利高清视频| 久久久久九九精品影院| 三级经典国产精品| 97人妻精品一区二区三区麻豆| 午夜福利视频1000在线观看| 国产黄片美女视频| 亚洲av第一区精品v没综合| 成人综合一区亚洲| 白带黄色成豆腐渣| www.色视频.com| 一级毛片久久久久久久久女| 久久久久久大精品| 麻豆精品久久久久久蜜桃| 久久精品影院6| 中国国产av一级| 日本黄大片高清| 色综合亚洲欧美另类图片| 国产黄片视频在线免费观看| 99国产极品粉嫩在线观看| 有码 亚洲区| 日韩欧美精品免费久久| 深夜精品福利| 欧美成人a在线观看| 成熟少妇高潮喷水视频| 尾随美女入室| 亚洲第一电影网av| 一级黄色大片毛片| 国产成人91sexporn| 九九热线精品视视频播放| 在线播放国产精品三级| 26uuu在线亚洲综合色| 日韩强制内射视频| 国产真实乱freesex| 伦理电影大哥的女人| 又粗又爽又猛毛片免费看| 午夜福利在线观看吧| 久久久久久久久久成人| 观看美女的网站| 日韩大尺度精品在线看网址| 97超视频在线观看视频| 久久久久久久久久久丰满| 级片在线观看| 久久这里只有精品中国| 成人鲁丝片一二三区免费| 国产女主播在线喷水免费视频网站 | 国产极品天堂在线| 久久久久久久久久黄片| 成年av动漫网址| 18+在线观看网站| 国产成年人精品一区二区| 国产精品久久电影中文字幕| 村上凉子中文字幕在线| 国产成人aa在线观看| 变态另类成人亚洲欧美熟女| 日本免费a在线| 乱码一卡2卡4卡精品| www.av在线官网国产| 日韩精品青青久久久久久| 99久久精品国产国产毛片| 九九爱精品视频在线观看| 中文字幕av在线有码专区| 亚洲一级一片aⅴ在线观看| 69av精品久久久久久| 日韩在线高清观看一区二区三区| 在线观看av片永久免费下载| 麻豆国产av国片精品| 国产老妇女一区| 久久6这里有精品| 尾随美女入室| 精品人妻偷拍中文字幕| 麻豆国产av国片精品| 女同久久另类99精品国产91| 亚洲高清免费不卡视频| 亚洲人成网站在线播| 亚洲精品久久久久久婷婷小说 | 99热这里只有是精品在线观看| 99热精品在线国产| 亚洲国产欧洲综合997久久,| 成人美女网站在线观看视频| 美女脱内裤让男人舔精品视频 | 欧美成人免费av一区二区三区| 国产精品1区2区在线观看.| 亚洲精品粉嫩美女一区| 一级毛片aaaaaa免费看小| 亚洲经典国产精华液单| 亚洲欧美日韩无卡精品| 久久99热6这里只有精品| 国产伦精品一区二区三区四那| 丰满人妻一区二区三区视频av| 免费人成视频x8x8入口观看| 成人高潮视频无遮挡免费网站| 成人漫画全彩无遮挡| 国产片特级美女逼逼视频| 国产在视频线在精品| 26uuu在线亚洲综合色| 日韩成人av中文字幕在线观看| 精品人妻熟女av久视频| 可以在线观看毛片的网站| 免费不卡的大黄色大毛片视频在线观看 | 免费看美女性在线毛片视频| 亚洲国产日韩欧美精品在线观看| 日产精品乱码卡一卡2卡三| 亚洲国产欧美人成| 国产伦精品一区二区三区四那| 五月玫瑰六月丁香| 免费观看精品视频网站| 精品久久久久久久久av| 成人av在线播放网站| 精品熟女少妇av免费看| 亚洲熟妇中文字幕五十中出| 联通29元200g的流量卡| 国内少妇人妻偷人精品xxx网站| 午夜精品在线福利| 日本黄大片高清| 午夜福利视频1000在线观看| 国产一级毛片七仙女欲春2| 亚洲欧美清纯卡通| 国产色婷婷99| 又粗又硬又长又爽又黄的视频 | 男女下面进入的视频免费午夜| 嘟嘟电影网在线观看| 国产伦在线观看视频一区| 激情 狠狠 欧美| 欧美xxxx性猛交bbbb| 成人欧美大片| 成年女人永久免费观看视频| 18+在线观看网站| 九九热线精品视视频播放| 麻豆一二三区av精品| 又粗又硬又长又爽又黄的视频 | 菩萨蛮人人尽说江南好唐韦庄 | 啦啦啦韩国在线观看视频| 爱豆传媒免费全集在线观看| 亚洲欧美日韩高清专用| 1000部很黄的大片| 人妻少妇偷人精品九色| 精品久久久久久久人妻蜜臀av| 天堂√8在线中文| 久久久精品大字幕| 亚洲国产欧美人成| 欧美日本视频| 亚洲av一区综合| 亚洲,欧美,日韩| 日韩欧美精品v在线| 亚洲一区高清亚洲精品| 日韩中字成人| 中国美女看黄片| 一级黄色大片毛片| 国内精品宾馆在线| av专区在线播放| 男女边吃奶边做爰视频| 99热6这里只有精品| 国产乱人视频| 精品无人区乱码1区二区| 日本五十路高清| 日韩一区二区三区影片| 国产精品乱码一区二三区的特点| 岛国在线免费视频观看| 欧美最新免费一区二区三区| 国产 一区精品| 亚洲自拍偷在线| 卡戴珊不雅视频在线播放| 在现免费观看毛片| av在线播放精品| 国产成人影院久久av| 看黄色毛片网站| 天堂√8在线中文| av在线老鸭窝| 高清在线视频一区二区三区 | 级片在线观看| 亚洲精品自拍成人| 美女脱内裤让男人舔精品视频 | 亚洲图色成人| av福利片在线观看| 久久久久九九精品影院| 91久久精品电影网| 99热6这里只有精品| 国产三级在线视频| 欧美极品一区二区三区四区| 久久久久久久久大av| 日本一本二区三区精品| 亚洲人成网站在线播放欧美日韩| 丝袜美腿在线中文| 99视频精品全部免费 在线| 高清毛片免费观看视频网站| av在线蜜桃| 国产精品美女特级片免费视频播放器| 最好的美女福利视频网| 中文字幕免费在线视频6| 99久久九九国产精品国产免费| 免费人成视频x8x8入口观看| 国产熟女欧美一区二区| АⅤ资源中文在线天堂| 午夜免费男女啪啪视频观看| 亚洲精品日韩av片在线观看| 精品久久久久久久末码| 不卡一级毛片| 精品无人区乱码1区二区| 国产精品福利在线免费观看| 免费搜索国产男女视频| 全区人妻精品视频| 一进一出抽搐动态| 亚洲国产精品sss在线观看| 精品久久久久久成人av| 亚州av有码| 国产美女午夜福利| 国产真实伦视频高清在线观看| 91在线精品国自产拍蜜月| 欧洲精品卡2卡3卡4卡5卡区| 国产精品,欧美在线| 欧美最黄视频在线播放免费| 嫩草影院精品99| 国产精品乱码一区二三区的特点| 成人无遮挡网站| 免费观看在线日韩| 噜噜噜噜噜久久久久久91| 欧美+日韩+精品| 国产视频内射| 亚洲人成网站在线播放欧美日韩| 超碰av人人做人人爽久久| 国产成人影院久久av| 国语自产精品视频在线第100页| 观看美女的网站| 国产在线精品亚洲第一网站| 国产一区二区在线观看日韩| 精品久久久久久久久亚洲| 免费av毛片视频| 国产精品蜜桃在线观看 | 床上黄色一级片| 国产高清三级在线| 国内揄拍国产精品人妻在线| 国产精品美女特级片免费视频播放器| 99视频精品全部免费 在线| 亚洲aⅴ乱码一区二区在线播放| 美女被艹到高潮喷水动态| 五月玫瑰六月丁香| 亚洲国产精品国产精品| 成人美女网站在线观看视频| av专区在线播放| 久久久久久久久久黄片| 亚洲第一电影网av| 久久精品国产自在天天线| 热99re8久久精品国产| 午夜福利高清视频| 精品欧美国产一区二区三| 午夜福利在线在线| 亚洲图色成人| 午夜福利高清视频| 国产v大片淫在线免费观看| 国产蜜桃级精品一区二区三区| 亚洲欧美日韩高清在线视频| 欧美一区二区亚洲| 三级男女做爰猛烈吃奶摸视频| 亚洲国产色片| 亚洲成人中文字幕在线播放| 老师上课跳d突然被开到最大视频| 国内少妇人妻偷人精品xxx网站| 亚洲一区高清亚洲精品| 国产精品一及| 少妇的逼水好多| 69av精品久久久久久| 亚洲欧美精品自产自拍| 美女xxoo啪啪120秒动态图| 一卡2卡三卡四卡精品乱码亚洲| 日日撸夜夜添| 久久久久性生活片| 九九久久精品国产亚洲av麻豆| 99热精品在线国产| 日日干狠狠操夜夜爽| 亚洲电影在线观看av| 久久精品国产亚洲网站| 婷婷精品国产亚洲av| 国产淫片久久久久久久久| 国产探花极品一区二区| 最后的刺客免费高清国语| 日韩一本色道免费dvd| av在线播放精品| 精品久久久久久久久久免费视频| 欧美性感艳星| 高清午夜精品一区二区三区 | 搡女人真爽免费视频火全软件| 亚洲国产色片| 免费观看a级毛片全部| 男插女下体视频免费在线播放| 亚洲精品自拍成人| 欧美成人a在线观看| 99久久久亚洲精品蜜臀av| 国产一级毛片在线| av在线蜜桃| 免费看日本二区| 日韩三级伦理在线观看| 国产精品伦人一区二区| 亚洲中文字幕一区二区三区有码在线看| 国产视频内射| 亚洲一级一片aⅴ在线观看| 亚洲成a人片在线一区二区| 亚洲欧美成人综合另类久久久 | 在线观看av片永久免费下载| 校园春色视频在线观看| 国产av不卡久久| 看免费成人av毛片| 99在线人妻在线中文字幕| 欧美+日韩+精品| 国产亚洲av嫩草精品影院| 亚洲自拍偷在线| 男女那种视频在线观看| 偷拍熟女少妇极品色| 你懂的网址亚洲精品在线观看 | 特级一级黄色大片| 日本欧美国产在线视频| 免费看光身美女| 亚洲av免费高清在线观看| 国产午夜精品久久久久久一区二区三区| 18禁裸乳无遮挡免费网站照片| 2022亚洲国产成人精品| 亚洲欧美清纯卡通| 搞女人的毛片| 1024手机看黄色片| 色噜噜av男人的天堂激情| 国产欧美日韩精品一区二区| 一个人观看的视频www高清免费观看| 91久久精品电影网| 亚洲精品乱码久久久v下载方式| 22中文网久久字幕| 国产成人a∨麻豆精品| 天堂av国产一区二区熟女人妻| 一个人观看的视频www高清免费观看| 久久热精品热| 18禁在线无遮挡免费观看视频| 在线国产一区二区在线| 日韩欧美精品免费久久| 久99久视频精品免费| 少妇被粗大猛烈的视频| 国产精品一区www在线观看| 国产精品国产高清国产av| 麻豆国产av国片精品| 26uuu在线亚洲综合色| 欧美xxxx性猛交bbbb| 亚洲成人av在线免费| 嫩草影院精品99| 在线天堂最新版资源| 日本黄色片子视频| 国产亚洲91精品色在线| 亚洲精品乱码久久久v下载方式| 亚洲av中文av极速乱| 国产成人影院久久av| 亚洲自拍偷在线| 日韩大尺度精品在线看网址| 免费搜索国产男女视频| 久久欧美精品欧美久久欧美| 乱码一卡2卡4卡精品| av免费观看日本| 国产美女午夜福利| 美女国产视频在线观看| 男人狂女人下面高潮的视频| 一进一出抽搐gif免费好疼| 欧洲精品卡2卡3卡4卡5卡区| 青春草亚洲视频在线观看| 99热这里只有是精品50| 岛国在线免费视频观看| 亚洲精品自拍成人| 亚洲四区av| 大香蕉久久网| 国产一区二区三区在线臀色熟女| 插逼视频在线观看| 岛国在线免费视频观看| 欧美+亚洲+日韩+国产| 真实男女啪啪啪动态图| 国产黄片视频在线免费观看| 中文字幕av在线有码专区| 一本久久精品| 看免费成人av毛片| or卡值多少钱| 国产午夜福利久久久久久| 精品99又大又爽又粗少妇毛片| 免费看a级黄色片| 亚洲欧美日韩高清在线视频| 夜夜爽天天搞| 人人妻人人澡欧美一区二区| 免费观看精品视频网站| 一区二区三区免费毛片| 村上凉子中文字幕在线| 久久这里有精品视频免费| 免费不卡的大黄色大毛片视频在线观看 | 麻豆乱淫一区二区| 日本黄大片高清|