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

    應力狀態(tài)和應變率對Q690鋼延性斷裂行為的影響

    2021-03-11 08:49董俊宏沙佳阿里阿里楊波丁浩敏
    土木建筑與環(huán)境工程 2021年1期

    董俊宏 沙佳阿里阿里 楊波 丁浩敏

    摘 要:研究了中國Q690高強度結構鋼在不同低應變率下的真實應力應變關系和延性斷裂特征。加載了3種不同開口尺寸的單軸拉伸平板試件和純剪試件,獲得了廣泛的應力狀態(tài)。用試驗結合數值模擬的方法確定等效塑性斷裂應變、應力三軸度和羅德參數。結果表明,應力三軸度和羅德參數是控制斷裂的主要影響因素。與應力相比,應變對應變率更為敏感,等效塑性斷裂應變隨應變率的增大而增大。此外,由于斷裂前純剪試件標距范圍發(fā)生明顯的變形和應力集中現象,斷裂面以拉應力或拉剪應力為主,純剪試件的設計需要進一步研究。

    關鍵詞:高強度結構鋼;延性斷裂;應力三軸度;羅德參數;應變率

    1 Introduction

    Due to the growing demand for the greater safety and lighter weight of structures, high-strength steel, such as the Chinese Q690 steel, is increasingly utilized in steel structures.

    A fracture in a member of steel buildings could lead to a reduction in the member's loading capacity, and consequently a structural failure could occur. Many studies have investigated the ductile fracture mechanism of steel structures. Researchers have carried out various experimental and numerical studies about the fracture behavior of steel that investigated ductile fractures and damage.

    In many studies, the stress triaxiality and the Lode parameter have been found to be the key parameters governing the ductile fracture of metallic materials. Deole et al.[1] conducted a series of metal fracture tests including tension and shear tests. A fracture locus was built, based on the stress triaxiality and the Lode parameter. They mentioned that the fracture locus depended on the calibration method used. Inspired by the mechanism of nucleation, growth and coalescence of voids, Lou et al.[2-4] developed a ductile fracture model for sheet specimens and extended it to incorporate the stress triaxiality and the Lode parameter. A changeable stress triaxiality cutoff value was also proposed, considering the various microstructures of different metals. Li et al.[5] fully investigated the ductile fracture characteristics of Chinese Q460 high strength structural steel under a quasi-static condition by using mechanical tests of four types of notched specimens. The influence of the stress state on the fracture mechanism of the material was then studied by observing the fracture surfaces of all test specimens using a Scanning Electron Microscope (SEM). Different fracture mechanisms were noticed in different stress triaxiality regions. Bai and Wierzbicki[6] found that ductile fracture was inherently a three-dimensional phenomenon and should be represented in the three-dimensional (3-D) space consisting of the stress triaxiality, the Lode parameter and the equivalent plastic strain to fracture. Sixteen fracture models were evaluated in this 3-D space.

    Other studies have also revealed that the strain rate can affect metal fracture behavior. Song et al.[7] studied the effect of strain rate on the fracture behavior of low-alloy structural carbon steel by investigating monotonic tensile fractures and in-plane shear fractures under different strain rates in a quasi-static loading range. The results showed that the strain rate slightly magnified the equivalent plastic fracture strain. Sjberg et al.[8] presented a methodology for fracture characterization at strain rates up to 1 000/s, at temperatures up to 650 ℃. A high-speed camera combined with digital image correlation was used to evaluate the strain at the fracture. A coupled relationship between the temperature and the stress triaxiality controlling the fracture strain was noticed. Roth and Mohr[9] carried out low, intermediate and high strain rate tensile experiments on flat, smooth, notched and central-hole specimens extracted from high strength steel sheets; an empirical modified Hosford-Coulomb fracture model was proposed and successfully validated. Khan and Liu[10] developed an isotropic ductile fracture criterion based on the magnitude of the stress vector (MSV) to include strain rate and temperature dependences. In this criterion, an exponential term was incorporated to simulate the strain rate's effect on the ductile fracture of an alloy.No.1 DONG Junhong,et al.:Stress state and strain rate dependent ductile fracture behavior of Q690 steelIn this study, uniaxial tension and pure shear tests were performed to investigate the stress state and strain rate dependent ductile behavior of Q690 high strength steel.One unnotched steel specimen and two steel specimens with different geometrical notches were used to conduct the uniaxial tests. Pure shear sheet specimens were adopted for the shear tests. The influence of quasi-static strain rate on the material constitutive model and fracture mechanism was studied. The experimental and numerical results indicated that the shape of pure shear specimens should be enhanced to ensure that fracture happens under the pure shear stress state.

    2 Experimental program

    2.1 Specimen design

    The layout of tensile and pure sheet specimens is illustrated in Fig.1. The tensile specimens include both notched and unnotched tensile sheets; all the sheets were fabricated from a structural steel beam with a nominal yielding strength of 690 MPa. Three samples were tested for each specimen to ensure the reliability of the data. Notch radii of R=2 mm and R=5 mm were used to obtain the results of different stress states; analogous notch radii data can be found in Reference [1] and [4]. The pure shear sheet specimen developed by Bao and Wierzbicki[11] was used in this study. Some modifications (e.g. a smooth arc radius of 5 mm at the gauge area) were made to avoid stress concentration, as depicted in Fig.1(d).

    2.2 Instrumentation

    The tensile and shear tests were performed using a universal testing machine (CMT5105), as shown in Fig.2. The force-displacement curve obtained directly from this machine was unreliable due to the possible slip between the clip and the specimen. Therefore, an extensometer was applied to measure the axial deformation in the gauge area until the specimen fractured. Two different quasi-static strain rates (0.000 1 and 0.01) were chosen. In order to achieve the expected strain rates, the loading rates of the universal testing machine were adjusted to 0.048 mm/min and 4.8 mm/min, respectively. All loading rates were kept constant during the entire test process. Also, a high-speed camera (VIC3D) was used to record the crack evolution of the external surface of the gauge area of the specimens.

    3 Finite element model (FEM)

    3.1 Material properties

    The unnotched tensile specimen (R0) was tested to determine the engineering stress-strain curves shown in Fig.3. The engineering stress is equal to the axial force divided by the initial area of the cross-section at the gauge area, whereas the engineering strain is equal to the axial elongation divided by the original gauge length. The specimens were cut using an accurate cutting machine (Dk7740). The dimension measurement showed that the error was no more than 0.2%. These curves were then converted to true stress-strain curves before necking occurred using Eq. (1) and Eq. (2).

    where σtrue is the true stress,εtrue is the true strain,σ is the engineering stress, and ε is the engineering strain.

    Previous studies have demonstrated that the hardening behavior of steel is significantly influenced by the strain rate. Similar work can be found in Reference[9] and [12]. One of the most famous metal hardening laws is the Johnson-Cook (JC) model[13], as shown in Eq. (3). The strain term of the JC model is always incorporated as a multiplied term in other hardening laws to include the effect of the strain rate[9].

    Another commonly-used hardening law is the Swift law[14], which was used in Reference [1] and [7] for sheet specimens. In this study, the strain rate of the JC model was included in the Swift law. The true stress-strain relationship from yielding to necking was fitted with the modified Swift hardening law as shown in Eq. (4).

    where D is a material parameter, m is the strain hardening exponent,ε0 is the first yield strain,εp is the plastic strain, εp· is the plastic strain rate, and ε*p· is the reference plastic rate. For the steel material used in this investigation, D=1 037.53 MPa, m=0.09 and C=0.002, which were obtained from numerical fitting of the measured data.

    3.2 Solid model

    The finite element software, Abaqus, was utilized to simulate the force-displacement curves. In order to save computational time and simplify the modelling process, only one eighth of the geometry of the tensile specimens (R0, R2 and R5) and half of the pure shear specimens (PS) were modelled due to the symmetrical deformation. The FEMs of the R0 and PS are illustrated in Fig.4 as the representative examples.

    A mesh dependency assessment was performed. Different element sizes were evaluated in terms of computational time and accuracy. The numerical results showed that a fine element size of 0.2 mm at the gauge area was reasonably appropriate for all the specimens. The reduced integration 8-node brick element (C3D8R) was chosen to mesh the specimens, and the implicit dynamic method was used.

    4 Results and discussion

    4.1 Fracture mode

    The experimental results of the fractured specimens showed that the macroscopic fracture morphology of the specimens was almost unaffected by the strain rate. Fig.5 shows the typical fracture modes of the R0, R2, R5 and PS specimens. Figs. 5 (a), (c), (e) and (g) were captured by the high-speed camera. A scanning electron microscope was used to obtain images of the fracture surface for these four types of specimens, as shown in Figs. 5(b), (d), (f) and (h).

    The R0 specimen fractured in the form of an inclined plane, because it was difficult to ensure the verticality of the specimen parallel to the load direction. The cracks initially appeared on the outer surface of the middle part of the specimen, as recorded by the high-speed video camera. In the R5 specimen, a horizontal fracture surface was observed, and cracks also began to appear in the middle part of the specimen. The fracture surface of the R2 specimen was a horizontal plane with the cracks initiating from the root of the notch, which was mainly attributed to the stress concentration by the smaller notch radius. In Figs. 5(b), (d), (f) and (h), we can see that all the fracture surfaces in the tensile specimens were comparatively rough and consistent with the micro-fracture mechanism of nucleation, growth and coalescence of voids. However, the pure shear specimen showed a smoother fracture surface characterized by shear decohesion.

    4.2 Determination of equivalent plastic strain and stress state at fracture

    The numerical results are displayed in Fig.6. It can clearly be seen that there is a difference between the force-displacement curves of the tensile specimens depicted in Figs. 6(a), (b) and (c) and the pure shear specimen's force-displacement curves in Fig.6(d). As for the uniaxial tensile specimens, there was a pronounced material softening behavior before the fracture. Also, different notch radii produced distinct elongations, and a minor notch radius resulted in less elongation. In contrast to the tensile specimens, the pure shear specimen did not experience a material softening process before a fracture, which meant that fracture occurred almost when the critical shear load was reached.

    In steel, the stress triaxiality and the Lode parameter are two key parameters that control the steel's fracture behavior. The stress triaxiality is defined by the ratio of the hydrostatic pressure to the Von Mises equivalent stress, as expressed in Eq. (5).

    The Lode parameter was calculated by using Eq. (6).

    where σ1,σ2,σ3 are the first, the second and the third principal stress, respectively.

    A hybrid experiment-simulation approach was used to determine the equivalent plastic strain, the stress triaxiality and the Lode parameter at each fracture. It is generally assumed that a fracture initiates when the force-displacement curve drops abruptly. It is worth noting that due to the discreteness of the test conditions, three specimens of the same type fractured at different displacements. The points where the experimental force-displacement curves in Fig.6 drop abruptly are known as the average fracture displacement. The maximum sample variance of the fracture displacement was 0.1, which indicates the good repeatability of all test results. Subsequently, the equivalent plastic strain, the stress triaxiality and the Lode parameter were calculated from the first element expected to be deleted in the finite element model at fracture.

    The first element expected to be deleted from the FEM is the hypothetical initial fracture point of the specimen. The determination of the hypothetical initial fracture point is mainly based on experimental crack evolution and the parameters related to the fracture (the equivalent plastic strain, the stress states). The location of this point is clear for the tensile specimens (R0, R2 and R5). Taking R0 as an example, it can be seen from Fig.5(a) that cracks initiated first in the middle part, meaning that the hypothetical initial fracture point was likely to lie at the center of the cross-section. Fig.7 shows the nephograms of the fracture-related parameters, such as the equivalent plastic strain, the maximum principal stress, the second principal stress and the minimal principal stress. For R0, the central element of the cross-section shows the maximum equivalent plastic strain and is most likely to be deleted first. At the same time, this part has the most unfavorable stress state, where the stress triaxiality is highest based on the three principal stresses, shown in Figs. 7(b), (c) and (d). As a result, it can reasonably be assumed that the initial fracture point of R0 is located at the center of the cross-section. The determination of the points for R2 and R5 follow the same principles. The hypothetical initial fracture point is located at the root of notch for R2, whereas it is at the center for R5.

    In terms of the stress state, it is clear from Table 1 that the influence of the strain rate on the stress triaxiality and the Lode parameter was negligible in the R2, R5 and PS specimens, where the stress triaxiality and the Lode parameter remained almost constant. However, there are some variations of the stress state in the R0 specimens under different strain rates. As for the equivalent plastic strain to fracture, it should be noted that the equivalent plastic strain to fracture slightly increases with the rise of strain rate for the tensile specimens, except R5. The small variation may be due to the slight change in the strain rate. Besides, the equivalent plastic strain at the fracture of the PS remained almost unchanged. Therefore, it can be concluded that the equivalent plastic strain to fracture is influenced by strain rate for tensile specimens. The reason why the fracture strain stays constant in the PS specimens could be that there are two different microstructure fracture mechanisms. The fracture behavior for tensile specimens is characterized by the nucleation, growth and coalescence of voids, while it is believed that shear decohesion is responsible for fractures in pure shear tests.

    4.4 Influence of PS specimen design on the stress state

    There are two geometrical defects, that caused the tensile stress on the failure surface of the pure shear specimen. The shear failure surface is expected to be parallel to the external force. However, due to the metal's ductility, a slant failure plane was obtained in most of the pure shear tests[6, 8, 15], where the shear-tensile stress or tension stress dominates rather than the pure shear stress. The deformation history of the PS specimen used in this study is shown in Fig.9; it shows that the PS specimen failed on an inclined plane rather than a horizontal plane.

    Another defect was the stress concentration caused by sudden changes of geometry in that gauge area that impeded the determination of fracture initiation. The pure shear test is used to determine the equivalent plastic fracture strain, stress triaxiality and the Lode parameter in pure shear stress state. However, the cracks started at the point where the shape suddenly changed, as shown in Fig.5(g). This indicates that the fracture was caused by stress concentration instead of the pure shear stress. The assumption that fracture occurs at the center of a cross-section, as previously adopted in many studies, is inappropriate because the fracture practically initiated at the edge points first according to the experimental results. Devising a method for determining the hypothetical initial fracture point of the PS deserves further study.

    In order to avoid the inclined fracture section and stress concentration, which usually lead to a non-pure shear stress state at the fracture point, the gauge area of the PS specimen should be designed carefully. Mohr and Henn[16] recommended using butterfly specimens to avoid the edge effect by using a smooth arc, as shown in Fig.9(a), so stress concentration could be weakened in this way; but this would lead to the slant failure plane as mentioned above. A short and wide gauge area tends to keep a fracture section horizontal. Consequently, the height of the gauge area that is perpendicular to the shear load should be lower. At the same time, a minor smooth arc radius is recommended for the design of gauge area.

    5 Conclusions

    Sheet specimens with different geometrical notches were loaded until fractures occurred to gain data about a wide range of stress states. Two quasi-static loading rates were selected. The effects of stress state and strain rate on the constitutive material and fracture behavior of Chinese Q690 structural steel was investigated. The main conclusions can be summarized as follows:

    1) The effect of the strain rate in the quasi-static range on the stress triaxiality and the Lode parameter can be disregarded.

    2) The equivalent plastic strain to fracture for the tensile specimen is affected by the strain rate. Usually, a high strain rate leads to a high equivalent plastic strain to fracture.

    3) For the pure shear specimen, the equivalent plastic strain at fracture was unaffected by the strain rate.

    4) The pure shear test indicates that the shape of gauge area can significantly influence the stress state and fracture initiation.

    Therefore, the height of the gauge area should be restricted and a small smooth arc should be used as the edge of gauge area to avoid the slant fracture section and stress concentration.

    Acknowledgements

    The authors gratefully acknowledge the financial support provided by the National Natural Science Foundation of China (Grant No. 51778086).

    References:

    [1] DEOLE A D, BARNETT M R, WEISS M. The numerical prediction of ductile fracture of martensitic steel in roll forming [J]. International Journal of Solids and Structures, 2018, 144/145: 20-31.

    [2] LOU Y S, HUH H, LIM S, et al. New ductile fracture criterion for prediction of fracture forming limit diagrams of sheet metals [J]. International Journal of Solids and Structures, 2012, 49(25): 3605-3615.

    [3] LOU Y S, HUH H. Extension of a shear-controlled ductile fracture model considering the stress triaxiality and the Lode parameter [J]. International Journal of Solids and Structures, 2013, 50(2): 447-455.

    [4] LOU Y S, CHEN L, CLAUSMEYER T, et al. Modeling of ductile fracture from shear to balanced biaxial tension for sheet metals [J]. International Journal of Solids and Structures, 2017, 112: 169-184.

    [5] LI W C, LIAO F F, ZHOU T H, et al. Ductile fracture of Q460 steel: Effects of stress triaxiality and Lode angle [J]. Journal of Constructional Steel Research, 2016, 123: 1-17.

    [6] BAI Y L, WIERZBICKI T. A comparative study of three groups of ductile fracture loci in the 3D space [J]. Engineering Fracture Mechanics, 2015, 135: 147-167.

    [7] SONG Q Y, HEIDARPOUR A, ZHAO X L, et al. Experimental and numerical investigation of ductile fracture of carbon steel structural components [J]. Journal of Constructional Steel Research, 2018, 145: 425-437.

    [8] SJBERG T, KAJBERG J, OLDENBURG M. Fracture behaviour of Alloy 718 at high strain rates, elevated temperatures, and various stress triaxialities [J]. Engineering Fracture Mechanics, 2017, 178: 231-242.

    [9] ROTH C C, MOHR D. Effect of strain rate on ductile fracture initiation in advanced high strength steel sheets: Experiments and modeling [J]. International Journal of Plasticity, 2014, 56(4): 19-44.

    [10] KHAN A S, LIU H W. Strain rate and temperature dependent fracture criteria for isotropic and anisotropic metals [J]. International Journal of Plasticity, 2012, 37: 1-15.

    [11] BAO Y B, WIERZBICKI T. On fracture locus in the equivalent strain and stress triaxiality space [J]. International Journal of Mechanical Sciences, 2004, 46(1): 81-98.

    [12] VERLEYSEN P, PEIRS J. Quasi-static and high strain rate fracture behaviour of Ti6Al4V [J]. International Journal of Impact Engineering, 2017, 108: 370-388.

    [13] JOHNSON G R, COOK W H. A constitutive model and data for metals subjected to large strains, high strain rates and high temperatures [C] // 7th International Symposium on Ballistics, 1983: 541-547.

    [14] SWIFT H W. Plastic instability under plane stress [J]. Journal of the Mechanics and Physics of Solids, 1952, 1(1): 1-18.

    [15] SJBERG T, MARTH S, KAJBERG J, et al. Experimental characterisation of the evolution of triaxiality stress state for sheet metal materials [J]. European Journal of Mechanics-A, 2017, 66: 279-286.

    [16] MOHR D, HENN S. Calibration of stress-triaxiality dependent crack formation criteria: A new hybrid experimental-numerical method [J]. Experimental Mechanics, 2007, 47(6): 805-820.

    (編輯 胡英奎)

    国产主播在线观看一区二区| 麻豆成人av在线观看| 97超级碰碰碰精品色视频在线观看| 日韩三级视频一区二区三区| 欧美黑人精品巨大| 日韩成人在线观看一区二区三区| 成年版毛片免费区| 国产精品久久电影中文字幕| 成人亚洲精品av一区二区| 精品久久久久久成人av| svipshipincom国产片| 法律面前人人平等表现在哪些方面| 久久久国产成人免费| 非洲黑人性xxxx精品又粗又长| 首页视频小说图片口味搜索| 欧美成人免费av一区二区三区| 欧美另类亚洲清纯唯美| 一区二区三区精品91| 国产一区二区在线av高清观看| 视频区欧美日本亚洲| 国产av一区二区精品久久| 国产精品,欧美在线| 亚洲真实伦在线观看| 日韩成人在线观看一区二区三区| 久久国产精品男人的天堂亚洲| 韩国精品一区二区三区| 亚洲 国产 在线| av中文乱码字幕在线| 亚洲成人免费电影在线观看| 999精品在线视频| 精品电影一区二区在线| 国产精品影院久久| 亚洲国产毛片av蜜桃av| 草草在线视频免费看| 久久精品aⅴ一区二区三区四区| 欧美黑人精品巨大| 18禁国产床啪视频网站| 国产精品二区激情视频| 99精品欧美一区二区三区四区| 国产三级黄色录像| 听说在线观看完整版免费高清| 国产伦人伦偷精品视频| cao死你这个sao货| 青草久久国产| 欧美日韩一级在线毛片| 97超级碰碰碰精品色视频在线观看| 亚洲国产欧美一区二区综合| 国产成人精品久久二区二区91| 桃红色精品国产亚洲av| 亚洲人成电影免费在线| 亚洲人成伊人成综合网2020| 中出人妻视频一区二区| 国产精品永久免费网站| 最新美女视频免费是黄的| 可以在线观看毛片的网站| 欧美绝顶高潮抽搐喷水| 在线观看66精品国产| 亚洲精品国产区一区二| 亚洲真实伦在线观看| avwww免费| 国产精品,欧美在线| 国产一区二区激情短视频| 午夜亚洲福利在线播放| a在线观看视频网站| 国产精品美女特级片免费视频播放器 | 天天添夜夜摸| 国产亚洲精品久久久久5区| 69av精品久久久久久| 久久九九热精品免费| 在线看三级毛片| 成人18禁在线播放| 人人澡人人妻人| 国产亚洲精品综合一区在线观看 | 成人国产综合亚洲| 国产视频一区二区在线看| 身体一侧抽搐| 欧美日本亚洲视频在线播放| 两个人免费观看高清视频| 国产精品日韩av在线免费观看| 一夜夜www| 一区二区三区精品91| 97超级碰碰碰精品色视频在线观看| 日韩欧美在线二视频| 亚洲精品av麻豆狂野| 91九色精品人成在线观看| 国产黄片美女视频| 亚洲国产欧美日韩在线播放| 日韩欧美国产一区二区入口| 国产精品久久久av美女十八| 久久婷婷人人爽人人干人人爱| 亚洲熟妇中文字幕五十中出| 欧美日韩乱码在线| 久久午夜亚洲精品久久| 99热6这里只有精品| www.精华液| 91在线观看av| 精品不卡国产一区二区三区| 91字幕亚洲| 日韩欧美一区视频在线观看| 天堂√8在线中文| 欧美成人性av电影在线观看| 黄片大片在线免费观看| 成年免费大片在线观看| 日本在线视频免费播放| 看片在线看免费视频| 国产精品98久久久久久宅男小说| 亚洲三区欧美一区| 久久久国产欧美日韩av| 宅男免费午夜| 99国产综合亚洲精品| 午夜福利视频1000在线观看| 亚洲成a人片在线一区二区| 欧美大码av| 91在线观看av| 亚洲熟女毛片儿| 国产av在哪里看| 国产精品一区二区三区四区久久 | 欧美日韩福利视频一区二区| 亚洲男人的天堂狠狠| 99国产精品99久久久久| 欧美日韩精品网址| 国产精品久久久久久亚洲av鲁大| 欧美激情 高清一区二区三区| 少妇裸体淫交视频免费看高清 | 亚洲av第一区精品v没综合| 国产亚洲av嫩草精品影院| 免费看日本二区| 制服人妻中文乱码| 色播在线永久视频| 99精品在免费线老司机午夜| www日本在线高清视频| 亚洲专区国产一区二区| 精品人妻1区二区| 亚洲av成人av| 中文字幕另类日韩欧美亚洲嫩草| 俺也久久电影网| 可以免费在线观看a视频的电影网站| 最近最新中文字幕大全电影3 | 久久久久久人人人人人| 国产又色又爽无遮挡免费看| 天堂√8在线中文| 日韩高清综合在线| 精品一区二区三区四区五区乱码| 成人18禁在线播放| 女人爽到高潮嗷嗷叫在线视频| 亚洲自偷自拍图片 自拍| 欧美亚洲日本最大视频资源| 少妇被粗大的猛进出69影院| 19禁男女啪啪无遮挡网站| 日本 欧美在线| 老司机午夜福利在线观看视频| 麻豆av在线久日| 男人的好看免费观看在线视频 | 久久久国产成人精品二区| 国产精品永久免费网站| 欧美绝顶高潮抽搐喷水| 日韩欧美 国产精品| 男人舔女人下体高潮全视频| 90打野战视频偷拍视频| 久久久久国内视频| 国产区一区二久久| 精品国产乱子伦一区二区三区| 少妇熟女aⅴ在线视频| 精品电影一区二区在线| 久久精品国产清高在天天线| 亚洲av第一区精品v没综合| 18禁黄网站禁片午夜丰满| 观看免费一级毛片| 老司机深夜福利视频在线观看| 国产精品综合久久久久久久免费| 国产蜜桃级精品一区二区三区| 精品国产乱码久久久久久男人| 国产真实乱freesex| 精品日产1卡2卡| 国产亚洲精品av在线| 国产精品久久久久久亚洲av鲁大| 亚洲av成人av| 99riav亚洲国产免费| 久久草成人影院| 一级毛片高清免费大全| 精品国产超薄肉色丝袜足j| 欧美黑人精品巨大| 91av网站免费观看| 黄色片一级片一级黄色片| 男女那种视频在线观看| 欧美 亚洲 国产 日韩一| 国产精品影院久久| 国产三级黄色录像| 真人一进一出gif抽搐免费| 国产又黄又爽又无遮挡在线| 日韩国内少妇激情av| 18禁观看日本| 亚洲国产欧洲综合997久久, | 韩国av一区二区三区四区| 夜夜看夜夜爽夜夜摸| 久热这里只有精品99| 亚洲一码二码三码区别大吗| 国产区一区二久久| 看黄色毛片网站| 在线观看66精品国产| 中文字幕人成人乱码亚洲影| 一个人免费在线观看的高清视频| 欧美色视频一区免费| 国产在线精品亚洲第一网站| 亚洲五月色婷婷综合| 黄色女人牲交| 亚洲一码二码三码区别大吗| 老司机午夜福利在线观看视频| netflix在线观看网站| 男人舔奶头视频| 久久香蕉国产精品| 真人一进一出gif抽搐免费| 久久青草综合色| 中文亚洲av片在线观看爽| 精品国产一区二区三区四区第35| 听说在线观看完整版免费高清| 无人区码免费观看不卡| 国产精品野战在线观看| 亚洲国产欧美日韩在线播放| 欧美午夜高清在线| 亚洲av第一区精品v没综合| 国产97色在线日韩免费| 9191精品国产免费久久| 亚洲全国av大片| 国产午夜精品久久久久久| 最近在线观看免费完整版| 日日夜夜操网爽| 身体一侧抽搐| АⅤ资源中文在线天堂| 一本一本综合久久| 久久午夜综合久久蜜桃| 亚洲五月色婷婷综合| 欧美午夜高清在线| 国产成人影院久久av| 日韩欧美在线二视频| 老司机深夜福利视频在线观看| 91成人精品电影| 男女之事视频高清在线观看| 很黄的视频免费| 精品国产国语对白av| 村上凉子中文字幕在线| 性色av乱码一区二区三区2| 精品久久蜜臀av无| 日本熟妇午夜| 在线视频色国产色| 久久中文字幕一级| 国产日本99.免费观看| 欧美日本视频| 欧美+亚洲+日韩+国产| 久久中文字幕一级| 亚洲人成电影免费在线| 久久午夜综合久久蜜桃| 一级黄色大片毛片| 国产高清激情床上av| 国产精品久久久久久人妻精品电影| 精品久久蜜臀av无| 欧美性长视频在线观看| 美女大奶头视频| 性色av乱码一区二区三区2| 精品欧美一区二区三区在线| 成人永久免费在线观看视频| 色在线成人网| 国产亚洲精品久久久久5区| 91麻豆av在线| 亚洲av五月六月丁香网| 日韩欧美一区视频在线观看| 国产精品久久视频播放| 制服人妻中文乱码| 999久久久精品免费观看国产| 91麻豆精品激情在线观看国产| 高清毛片免费观看视频网站| 国产片内射在线| 久久婷婷成人综合色麻豆| 女性生殖器流出的白浆| 欧美色欧美亚洲另类二区| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美大码av| 制服诱惑二区| 欧美一级a爱片免费观看看 | 伊人久久大香线蕉亚洲五| 琪琪午夜伦伦电影理论片6080| 欧美黄色淫秽网站| 一本大道久久a久久精品| 国产精品香港三级国产av潘金莲| 夜夜爽天天搞| 国产av一区在线观看免费| 日韩中文字幕欧美一区二区| 亚洲成人国产一区在线观看| 99riav亚洲国产免费| 高潮久久久久久久久久久不卡| 久久精品人妻少妇| 国产精品精品国产色婷婷| 成人三级做爰电影| 国产精品一区二区精品视频观看| 国产91精品成人一区二区三区| 欧美日韩一级在线毛片| 一区二区三区国产精品乱码| 亚洲成人精品中文字幕电影| 搡老妇女老女人老熟妇| √禁漫天堂资源中文www| 亚洲中文字幕一区二区三区有码在线看 | 国产日本99.免费观看| 欧美绝顶高潮抽搐喷水| 国产主播在线观看一区二区| 久久性视频一级片| 热re99久久国产66热| 国产又爽黄色视频| 久久久久国内视频| 中文字幕久久专区| 亚洲精品一区av在线观看| 69av精品久久久久久| 精品第一国产精品| 亚洲专区中文字幕在线| av在线天堂中文字幕| 成人三级做爰电影| 亚洲中文字幕日韩| 女性生殖器流出的白浆| 亚洲免费av在线视频| 老熟妇乱子伦视频在线观看| 精品国产美女av久久久久小说| 麻豆国产av国片精品| 国产成+人综合+亚洲专区| 免费av毛片视频| 免费在线观看黄色视频的| 日日夜夜操网爽| 国产一区在线观看成人免费| 精品久久久久久久末码| 欧美最黄视频在线播放免费| 色在线成人网| 国产成人影院久久av| 一区二区三区国产精品乱码| 国产真人三级小视频在线观看| 在线观看免费日韩欧美大片| 热re99久久国产66热| 男女下面进入的视频免费午夜 | 啦啦啦韩国在线观看视频| 国产高清videossex| 熟女电影av网| 精品久久久久久久毛片微露脸| 精品久久久久久成人av| 欧美三级亚洲精品| 亚洲精品粉嫩美女一区| 精品日产1卡2卡| 日本黄色视频三级网站网址| 国产一区二区三区在线臀色熟女| 精品欧美一区二区三区在线| 不卡一级毛片| 性色av乱码一区二区三区2| 亚洲在线自拍视频| av中文乱码字幕在线| 性色av乱码一区二区三区2| 午夜久久久在线观看| 别揉我奶头~嗯~啊~动态视频| 美女 人体艺术 gogo| 久久性视频一级片| 亚洲 欧美一区二区三区| 午夜精品在线福利| 欧美av亚洲av综合av国产av| 国内精品久久久久精免费| 在线观看免费日韩欧美大片| 在线天堂中文资源库| 国产成人精品久久二区二区91| 人人妻人人澡人人看| 国产99久久九九免费精品| 手机成人av网站| 免费在线观看影片大全网站| 精品国产一区二区三区四区第35| av欧美777| 99在线人妻在线中文字幕| 国产91精品成人一区二区三区| 亚洲国产高清在线一区二区三 | 老司机深夜福利视频在线观看| 一级a爱片免费观看的视频| 国产高清有码在线观看视频 | 午夜福利免费观看在线| 亚洲成国产人片在线观看| 国产黄片美女视频| 亚洲国产欧洲综合997久久, | 妹子高潮喷水视频| 日韩三级视频一区二区三区| 欧美色欧美亚洲另类二区| 老熟妇乱子伦视频在线观看| 草草在线视频免费看| 桃红色精品国产亚洲av| 老司机午夜福利在线观看视频| 好看av亚洲va欧美ⅴa在| 国产av在哪里看| 宅男免费午夜| 久久精品国产亚洲av高清一级| 成熟少妇高潮喷水视频| 婷婷亚洲欧美| 午夜福利在线在线| 精品乱码久久久久久99久播| 一夜夜www| 人人妻人人澡人人看| 脱女人内裤的视频| 欧美精品啪啪一区二区三区| 亚洲片人在线观看| 亚洲国产欧美一区二区综合| aaaaa片日本免费| 一本大道久久a久久精品| 香蕉av资源在线| 最近在线观看免费完整版| 黑人欧美特级aaaaaa片| 国产成人欧美| 香蕉av资源在线| 欧美日韩一级在线毛片| 欧美人与性动交α欧美精品济南到| 亚洲,欧美精品.| 黄网站色视频无遮挡免费观看| 欧美绝顶高潮抽搐喷水| 亚洲人成伊人成综合网2020| 在线免费观看的www视频| 久久精品91蜜桃| 国产精品久久久久久亚洲av鲁大| 男人舔女人下体高潮全视频| 亚洲第一青青草原| 搡老妇女老女人老熟妇| 日本一区二区免费在线视频| 国产激情偷乱视频一区二区| av片东京热男人的天堂| 欧美国产精品va在线观看不卡| 两个人免费观看高清视频| 亚洲人成电影免费在线| 日韩国内少妇激情av| 国产精品一区二区三区四区久久 | 99在线视频只有这里精品首页| 国产伦在线观看视频一区| 亚洲片人在线观看| 99久久国产精品久久久| 国产精品亚洲一级av第二区| 国产精品一区二区三区四区久久 | 成人三级做爰电影| 成年版毛片免费区| 波多野结衣巨乳人妻| 曰老女人黄片| 91成人精品电影| 亚洲欧美一区二区三区黑人| 美女扒开内裤让男人捅视频| 国产亚洲av高清不卡| 国产一区二区激情短视频| 欧美成人免费av一区二区三区| 久久精品人妻少妇| 国产高清有码在线观看视频 | 欧美性猛交黑人性爽| 精品无人区乱码1区二区| 啦啦啦韩国在线观看视频| 999久久久精品免费观看国产| 精品国产美女av久久久久小说| 两个人视频免费观看高清| 日本一区二区免费在线视频| 中文字幕高清在线视频| 日本五十路高清| 中亚洲国语对白在线视频| 亚洲av成人一区二区三| 久久久国产欧美日韩av| 午夜精品在线福利| 国产精品久久久久久亚洲av鲁大| 欧美一级毛片孕妇| 国产av一区二区精品久久| 最近最新中文字幕大全电影3 | 级片在线观看| 免费高清视频大片| 一本综合久久免费| 免费看十八禁软件| 12—13女人毛片做爰片一| 午夜老司机福利片| 97人妻精品一区二区三区麻豆 | 精华霜和精华液先用哪个| 国内久久婷婷六月综合欲色啪| 黄片小视频在线播放| 中文字幕久久专区| 欧美激情高清一区二区三区| 国产黄片美女视频| 亚洲一卡2卡3卡4卡5卡精品中文| 日本五十路高清| 亚洲男人的天堂狠狠| 免费电影在线观看免费观看| 丁香六月欧美| 国产伦人伦偷精品视频| 久久国产精品影院| 别揉我奶头~嗯~啊~动态视频| 国产三级在线视频| 法律面前人人平等表现在哪些方面| 美女高潮喷水抽搐中文字幕| 欧美日本视频| 色播亚洲综合网| 不卡一级毛片| 中文字幕精品亚洲无线码一区 | 精品一区二区三区四区五区乱码| 亚洲国产看品久久| 国产精品一区二区免费欧美| 91成年电影在线观看| 欧美国产精品va在线观看不卡| 高清毛片免费观看视频网站| 男人舔奶头视频| 99久久综合精品五月天人人| 国产私拍福利视频在线观看| 亚洲av中文字字幕乱码综合 | 亚洲中文字幕一区二区三区有码在线看 | 不卡av一区二区三区| 亚洲五月天丁香| 久久精品91蜜桃| 日本撒尿小便嘘嘘汇集6| 亚洲国产看品久久| 日韩大尺度精品在线看网址| 午夜亚洲福利在线播放| 欧美三级亚洲精品| 日韩一卡2卡3卡4卡2021年| 国产不卡一卡二| 日本一本二区三区精品| 久久中文字幕人妻熟女| 亚洲精品一卡2卡三卡4卡5卡| 欧美成狂野欧美在线观看| 神马国产精品三级电影在线观看 | 国产乱人伦免费视频| 国产精品精品国产色婷婷| www.熟女人妻精品国产| 波多野结衣巨乳人妻| 日韩大尺度精品在线看网址| 欧美+亚洲+日韩+国产| 精品久久久久久久毛片微露脸| 亚洲天堂国产精品一区在线| 亚洲人成77777在线视频| 国产精品国产高清国产av| 男人的好看免费观看在线视频 | 国产蜜桃级精品一区二区三区| 精品国内亚洲2022精品成人| 欧美日韩亚洲综合一区二区三区_| 欧美日韩精品网址| 别揉我奶头~嗯~啊~动态视频| 18禁黄网站禁片免费观看直播| 中文在线观看免费www的网站 | 性欧美人与动物交配| 亚洲,欧美精品.| 美国免费a级毛片| 中国美女看黄片| 亚洲全国av大片| 午夜久久久在线观看| 亚洲精品国产精品久久久不卡| 欧美性长视频在线观看| 久久九九热精品免费| 午夜影院日韩av| 黑人欧美特级aaaaaa片| 男女之事视频高清在线观看| 亚洲精品在线美女| 国产v大片淫在线免费观看| 亚洲精品色激情综合| 欧美乱码精品一区二区三区| 在线观看免费日韩欧美大片| 亚洲国产毛片av蜜桃av| 一级毛片女人18水好多| 日本a在线网址| 日韩欧美 国产精品| 精品欧美一区二区三区在线| 亚洲久久久国产精品| 成人永久免费在线观看视频| 桃红色精品国产亚洲av| 亚洲国产精品999在线| xxx96com| 欧美成人性av电影在线观看| 久久婷婷成人综合色麻豆| 成年女人毛片免费观看观看9| 亚洲男人天堂网一区| 欧美性长视频在线观看| 欧美日韩一级在线毛片| 两个人免费观看高清视频| 在线国产一区二区在线| 夜夜爽天天搞| av天堂在线播放| 日韩欧美国产在线观看| 中文字幕另类日韩欧美亚洲嫩草| 亚洲电影在线观看av| 最新在线观看一区二区三区| 国产精品免费一区二区三区在线| 搡老熟女国产l中国老女人| 成人手机av| 国产成人系列免费观看| 成人一区二区视频在线观看| 99国产精品一区二区蜜桃av| 欧美+亚洲+日韩+国产| 制服人妻中文乱码| 日韩欧美免费精品| 国产一级毛片七仙女欲春2 | 久久久水蜜桃国产精品网| 老司机在亚洲福利影院| 日本撒尿小便嘘嘘汇集6| 国产成年人精品一区二区| 欧美激情久久久久久爽电影| 亚洲成人久久爱视频| 久久婷婷人人爽人人干人人爱| 欧美日本视频| 制服诱惑二区| 亚洲国产看品久久| 午夜日韩欧美国产| 18禁裸乳无遮挡免费网站照片 | 国产麻豆成人av免费视频| 亚洲激情在线av| 久久精品亚洲精品国产色婷小说| 午夜精品在线福利| 久热这里只有精品99| 好看av亚洲va欧美ⅴa在| 99在线视频只有这里精品首页| 两个人免费观看高清视频| 成人永久免费在线观看视频| 欧美日韩一级在线毛片| 国产91精品成人一区二区三区| 欧美日韩瑟瑟在线播放| 老司机午夜十八禁免费视频| 亚洲av美国av| 日韩精品中文字幕看吧| 99久久国产精品久久久| 久久精品国产清高在天天线| 老司机深夜福利视频在线观看| 国产国语露脸激情在线看|