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

    Effects of thawing-induced softening on fracture behaviors of frozen rock

    2024-03-25 11:06:26TingWngHilingJiQingSunXinjunTnLiyunTng

    Ting Wng,Hiling Ji,Qing Sun,Xinjun Tn,Liyun Tng

    a College of Geology and Environment, Xi’an University of Science and Technology, Xi’an, 710054, China

    b College of Architecture and Civil Engineering, Xi’an University of Science and Technology, Xi’an, 710054, China

    c State Key Laboratory of Geomechanics and Geotechnical Engineering,Institute of Rock and Soil Mechanics,Chinese Academy of Sciences,Wuhan,430064,China

    Keywords: Frozen sandstone Different thawing temperature Fracture toughness Microcrack propagation process Unfrozen water content

    ABSTRACT Due to the presence of ice and unfrozen water in pores of frozen rock,the rock fracture behaviors are susceptible to temperature.In this study,the potential thawing-induced softening effects on the fracture behaviors of frozen rock is evaluated by testing the tension fracture toughness (KIC) of frozen rock at different temperatures(i.e.-20 °C,-15 °C,-12 °C,-10 °C,-8 °C,-6 °C,-4 °C,-2 °C,and 0 °C).Acoustic emission (AE) and digital image correlation (DIC) methods are utilized to analyze the microcrack propagation during fracturing.The melting of pore ice is measured using nuclear magnetic resonance(NMR)method.The results indicate that:(1)The KIC of frozen rock decreases moderately between-20 °C and -4 °C,and rapidly between -4 °C and 0 °C.(2) At-20 °C to-4 °C,the fracturing process,deduced from the DIC results at the notch tip,exhibits three stages: elastic deformation,microcrack propagation and microcrack coalescence.However,at -4 °C-0 °C,only the latter two stages are observed.(3) At-4 °C-0 °C,the AE activities during fracturing are less than that at -20 °C to -4 °C,while more small events are reported.(4) The NMR results demonstrate a reverse variation trend in pore ice content with increasing temperature,that is,a moderate decrease is followed by a sharp decrease and-4 °C is exactly the critical temperature.Next,we interpret the thawing-induced softening effect by linking the evolution in microscopic structure of frozen rock with its macroscopic fracture behaviors as follow:from-20 °C to-4 °C,the thickening of the unfrozen water film diminishes the cementation strength between ice and rock skeleton,leading to the decrease in fracture parameters.From-4 °C to 0 °C,the cementation effect of ice almost vanishes,and the filling effect of pore ice is reduced significantly,which facilitates microcrack propagation and thus the easier fracture of frozen rocks.

    1.Introduction

    Increasing global warming now leads to numerous large-scale,high-potential rock slope failures,rockfalls,and other geological hazards in alpine regions(Gruber and Haeberli,2007;Gobiet et al.,2014),and it brings serious threats to infrastructure,human life,and property safety(Phillips,2006;Huggel et al.,2008).Numerous evidences indicate that melting of the permafrost in high mountains is a major factor for destabilization of rock slope(Huggel et al.,2010;Ravanel et al.,2010;Etzelmüller et al.,2022).Mountain permafrost is generally characterized by thick frozen bedrock,and the mechanical properties of which primarily determine the stability of rock slopes in high mountains (Harris et al.,2009;Phillips et al.,2017).Fracture toughness(KC),a common mechanical parameter of frozen rock,can be considered as a criterion of fracture initiation and propagation to estimate the instability of rock slopes (Fischer et al.,2010;Krautblatter et al.,2013).Due to the coexistence of unfrozen water and ice in pores of frozen rock,its fracture behaviors are susceptible to temperature change.Thus,investigating the fracture behaviors of frozen rock during thawing is crucial for engineering activities in mountain permafrost regions.

    In fracture mechanics,there are basically three modes of rock fracturing (Whittaker et al.,1992): tension (mode I),sliding shear(mode II),and tearing shear mode (mode III).As a result of the comparatively low tensile strength of rocks,tension fracture is the most common in failure analysis of rock slopes (Stead and Wolter,2015;Chen et al.,2021).In addition,tension fracture is encountered in other areas of rock engineering(He et al.,2010;Zhou et al.,2018),including hydraulic fracturing (Pan et al.,2021),shale gas extraction (Li et al.,2019),and deep nuclear waste (Peng et al.,2020).Therefore,numerous researches have been conducted on the tension fracture properties of rocks,including test methods of tension fracture toughnessKIC(Mu?oz-Ibá?ez et al.,2021;Kramarov et al.,2020),various influence factors (Cai et al.,2020;Deng et al.,2021)and fracturing process(Lin et al.,2020;Han et al.,2022).

    Previous researches have provided a basic understanding of the tension fracture process of rock.The International Society of Rock Mechanics and Rock Engineering (ISRM) proposed four standard ways to testKIC,namely,the short rod (SR),chevron bend (CB),cracked chevron notched Brazilian disc(CCNBD)and semi-circular bend(SCB)methods(Kuruppu et al.,2014).Studies on the influence factors ofKICshow that in addition to lithology and water content,temperature also has a significant influence.At high or low temperatures,the fracture toughness of rocks differs greatly from that at room temperature(Dwivedi et al.,2000;Feng et al.,2017;Miao et al.,2020).For the tension fracture,due to the singularity of cracks,a rather extensive plastic zone,termed as the fracture process zone(FPZ),is generated by microcracks around the cracks’tips.Variation of the FPZ directly influences the fracture propagation process(Pan et al.,2021).Methods such as digital image correlation(DIC) and acoustic emission (AE) (Wong and Guo,2019;Lin et al.,2020) have been frequently employed to determine the range of FPZ and investigate its evolution(Wei et al.,2015;Wu et al.,2011;Lin and Labuz,2013).

    While previous researches on frozen rock are mainly focused on the strength,deformability and failure features(Jia et al.,2020a)or the effects of freeze-thaw cycles on its mechanical properties(Jamshidi et al.,2017;Jia et al.,2020b;Zhu et al.,2021;Jamshidi,2021;Wang et al.,2022;Ren et al.,2022),little attention is paid to the fracture behaviors of frozen rock.Heating of frozen rock leads to melting of the pore ice and hence alters the pore structure of the rock,which will have a substantial effect on the fracture mechanical properties of frozen rock (Wang et al.,2022).Investigating the fracture behaviors of frozen rock becomes challenging.However,a proper evaluation of the effects of thawing on the fracture features of frozen rock is rarely reported.For this,a series of semi-circular bending (SCB) tests of frozen rocks at different thawing temperatures(T=-20°C,-15°C,-12°C,-10°C,-8°C,-6°C,-4°C,-2°C and 0°C) were conducted.AE and DIC techniques were then utilized to monitor the AE signals and surface strain of the sample during microcrack propagation,and the impact of thawing was also analyzed.The nuclear magnetic resonance(NMR)method was used to examine the change in unfrozen water content of rock.Combined with the change of unfrozen water content and microcrack propagation process,the mechanism of how thawing modifies the fracture behaviors of frozen rock was discussed.

    2.Methodology

    2.1.Test method of KIC

    The SCB samples are used in this study to evaluate theKICof frozen rock due to the easy preparation and testing procedure(Kuruppu et al.,2014).The size of the SCB sample is determined as shown in Fig.1a,and theKICis calculated according to Eq.(1)-(2)in Kuruppu et al.(2014).

    Fig.1.(a) Size of SCB sample;and (b) SCB rock samples.

    2.2.Digital image correlation (DIC) method

    DIC can be used to accurately measure optical surface deformation.The device consists of a light,a charge coupled device(CCD) camera,and a computing apparatus with relevant software.The correlation calculation approach is used to calculate the deformation information by analyzing the two images before and after deformation.The deformation analysis process of DIC is as follows:

    (1) A pre-deformation image is selected as the reference image and the region of interest is determined (see Fig.2);

    Fig.2.The DIC method to calculate deformation.

    (2) By specifying the subset size and distance,the region of interest is meshed to obtain uniformly spaced virtual grids;

    (3) A group of pixels on each point of the virtual grid is defined as a subset.It is assumed that each subset is in a rigid motion,ignoring the deformation of the subset itself.Then,correlation calculations are performed on each subset using a predefined correlation function and certain search strategies.In the deformed image,the region that has the strongest correlation with the reference subset is determined as the target subset,which marks the position of the reference subset after deformation,and the displacement between them is obtained accordingly(Fig.2) (Pan and Li,2011);and

    (4) The deformation information of the whole field can be obtained by computing all subsets.A zero-mean normalized sum of squared difference is utilized as the similarity measure of the reference and the target subsets in this study due to its robustness to scale and offset variations of illumination as well as its simplicity for optimization (Pan et al.,2010).

    2.3.NMR method

    If a rock sample is in the permanent magnet field of the NMR equipment,the hydrogen (H) protons of pore water generate magnetization vectors along the direction of the magnetic field.While a pulsed electromagnetic field is emitted,these H protons flip to a new equilibrium position,and then back to its original equilibrium position as the alternating electromagnetic field is withdrawn.This process is called relaxation.During the relaxation,the H protons will cut the magnetic flux line and produce an electrical signal,and the electrical signals continues to decay,which is called free induction decay (FID).The NMR signals of FID are positively correlated with the number of water molecules (Zhao et al.,2017).While,the NMR signals of H protons in ice can hardly be detected in low-intensity magnetic field.Therefore,the NMR signals of FID can be used to determine the liquid water (i.e.unfrozen water) content in frozen rock and soils (Tice et al.,1982;Watanabe and Wake,2009):

    The ice content can be calculated by

    whereSrwis the relative unfrozen water content(%),FIDfrozenis the NMR signals of FID at different thawing temperatures,FIDunfrozenis the FID value at 20°C,andSriis the relative ice content (%).

    2.4.Experimental design

    2.4.1.Sample preparation

    (1)Petrographic studies.The sample used in this study is yellow feldspar quartz sandstone bearing a massive structure without obvious bedding layers.In addition,it is composed of 90% quartz,5% feldspar sand and 5% Kaolinite.Fig.3 exhibits the single polarized and orthogonal polarized microscopic thin-section images of the sandstone.

    Fig.3.(a) The single polarized;and (b) orthogonal polarized microscopic thin-section images.

    (2)SCB sample preparation.Twenty two cylindrical samples(diameter×height:100 mm×45 mm)were drilled from the same sandstone block.They were oven-dried at 105°C and placed in a vacuum-saturated cylinder to pump for 4 h under-0.1 MPa.Water was then injected into the saturated cylinder to completely immerse the sample,the water-saturated mass of samples was measured after 72 h,and their porosity was determined.The Pwave test of each sample was also conducted.Among the 22 cylindrical samples,18 samples with similar porosity and P-wave velocity were selected and sliced into 36 semicircle SCB samples(Fig.1b).The dry density of samples was (2.1 ± 0.1) g/cm3,the porosity of samples was(19.18 ±0.2)%,and the P-wave velocity of samples was (2.37 ±0.6) km/s.

    Thirty six saturated SCB samples were wrapped with plastic films and placed in a moisturizing vat(at 20°C,relative humidity of 50%) for a week to ensure uniform water distribution inside the samples.The samples were then frozen for 48 h in an environmental chamber at -20°C.In order to conduct the DIC test,the sample’s surface should be treated before the test to show random speckle patterns.Therefore,the surface of specimens was first painted white;once the white paint dried,black dots were dotted with a marker to create randomly scattered speckle patterns.

    2.4.2.Tests of KIC

    TheKICtests at different temperatures were carried out via a hydraulic servo press,which is equipped with a high-precision environmental chamber (Fig.4a).Loading was controlled by the axial displacement at a rate of 0.12 mm/min,which satisfies the ISRM’s recommendation that the loading rate of the SCB test should not exceed 0.2 mm/min.The chamber has a temperature control accuracy of ±0.1°C.

    Fig.4.(a) The press machine which is fitted out with a high-precision environment chamber,and the AE and DIC strain capture systems;and (b) NMR test system.

    The lower temperatures limit in cold regions (except in polar regions)is usually above-20°C(Ma et al.,2022),thus it is rational to set -20°C as the starting point of thawing for frozen rock.TheKICof 36 SCB saturated samples at nine thawing temperatures(T=-20°C,-15°C,-12°C,-10°C,-8°C,-6°C,-4°C,-2°C,and 0°C)were tested.The testing steps are as follows:a frozen sample is taken out from the-20°C chamber,and four AE probes are fixed on the sample according to the positions given in Fig.1a,and it is then placed back in the chamber equipped to the press,as shown in Fig.4a.Afterwards,the temperature of the chamber was set as a preset test temperature and kept for 4 h,and theKICtest was conducted.

    The DS5-8B AE instrument (Beijing Softland Times Technology Co.,Ltd.)was employed to record AE signals during sample loading.The AE system can collect data at a rate of 3 Mhz (Fig.4a).The preamplifier and trigger threshold of AE acquisition system are both set as 40 dB.

    The variation of surface strain during loading was analyzed using the MatchID-2d/Stereo full field strain measurement system(Fig.4a).A photon FASTCAM Mini UX100 camera was vertically fixed at a tripod 0.5 m away from the press machine.The frame rate is 50 fps(frames per second),and the light source is LED cold light.

    2.4.3.Measurement of the unfrozen water content in frozen rock during thawing

    The MacroMR12-150H NMR instrument (Suzhou Niumai Analyzer Co.,Ltd.)is used in this study(Fig.4b).A low temperature control system is equipped with NMR receiving coil,and it can provide a uniform subzero temperature field.The parameters for NMR test is the same as that of Wang et al.(2021).The following are the steps for measuring the unfrozen water content of frozen rock during thawing:(1)A sample was firstly placed in the receiving coil and the temperature is set as-20°C;(2)After 4 h,the NMR test is conducted;(3)The temperatureTis then set to-15°C,-12°C,-10°C,-8°C,-6°C,-5°C,-4°C,-3°C,-2°C,-1°C,0°C and 20°C,successively,and (4) the NMR signal is measured at each temperature after 4 h.

    3.Results

    3.1.Changes in KIC and fracture energy during thawing

    TheKICrepresents the capacity of a rock to resist unstable propagation of cracks.The fracture energy refers to as the required energy per unit area of crack propagation,which is dissipated during crack propagation.The two parameters are used to describe the initiation and propagation of cracks for assessing the fracture characteristics of rocks.

    TheKICof frozen rock at various thawing temperatures is displayed in Fig.5a.It drops remarkably due to thawing and exhibits a two-stage mode of (a) slow decrease,KICdeclined by 19.5% asTdrops from-20°C to-4°C;and(b)rapid decrease,KICdeclined by 74%as-4°C ≤T<0°C.Therefore,-4°C is a critical temperature that marks the starting of rapid decrease inKICof frozen rock during thawing.

    Fig.5.(a) Changes of KIC and (b) fracture energy during thawing.

    During the SCB sample fracture,the fracture energy can be computed by the work done by the press system on the sample during loading,that is,the area is enclosed by the axial loaddisplacement curve before fracture (Peng et al.,2020).Based on this,the fracture energy of samples at various temperatures is displayed in Fig.5b.It shows that the fracture energy falls dramatically as the temperature increases and also exhibits a twostage mode of being essentially unchanged from-20°C to-12°C and rapidly reduces from-12°C to 0°C.

    3.2.Influence of thawing on fracture process of frozen rock

    Fig.6 illustrates the evolution of horizontal strain (HS) at the notch tip of SCB samples during fracture at various temperatures.It shows that the HS increases with the rise of axial stress,and thawing plays a dominant role.As -20°C ≥T>-4°C (Fig.6a-f),the HS vs.displacement curve exhibits a three-stage mode of (I)changeless stage;(II) moderate increase stage;and (III) rapid increase stage.At-4°C ≤T<0°C(Fig.6g-i),it exhibits a two-stage mode of(II)moderate increase stage,and(III)rapid increase stage.The mode I stage is not observed.

    Fig.6.(a-i)Variations in HS at notch tip during fracture at different temperatures of-20 °C;-15 °C;-12 °C;-10 °C;-8 °C;-6 °C;-4 °C;-2 °C;and 0 °C;and(j)HS of notch tip under peak load at different temperatures.

    (1) At -20°C ≥T>-4°C (Fig.6a-f),in the changeless stage(mode I),the HS is essentially unaffected by increase in the displacement and horizontal strain field (HSF) is basically evenly distributed(HSF at pointAon the HS vs.displacement curve).This suggests that the sample is undergoing a vertically elastic compression stage and the horizontal strain localization does not appear.In the moderate increase stage(mode II),the HS increases linearly as the displacement increases,and the HS around the notch tip is much higher than that in other regions(pointB),and this local region is the FPZ.The presence of FPZ implies the occurrence of strain localization due to expansion of the microcracks.In the rapid increase stage (mode III),the HS grows rapidly and the FPZ expands substantially(pointC).It demonstrates that the FPZ gradually develops,and as the microcrack coalesces,the precut macrocracks grow.Therefore,when the temperature is between-20°C and-4°C,the three stages of HS variation actually correspond to three stages in microcrack evolution at the notch tip of frozen rock: elastic deformation stage,microcrack propagation stage,microcrack coalescence stage.

    (2) At -4°C ≤T<0°C (Fig.6g-i),in mode II stage,the HS increases linearly and the FPZ(pointB)occurs,indicating that microcracks have propagated.At the mode III stage,the HS grows rapidly and the FPZ range expands significantly(pointC),which is the precursor of microcrack propagation.The absence of the changeless stage suggests that microcracks at the notch tip of frozen rock are more likely to propagate at higher temperatures.

    (3) Aiming to understand the effect of thawing on microcrack propagation in frozen rock,the HS at the notch tip of peak load (pointC) is measured at various temperatures,as depicted in Fig.6j.The HS at the notch tip grows approximately linearly with temperature increasing from -20°C to-4°C,while the increase in HS accelerates as temperature surpasses-4°C.This result indicates that thawing induces a significant softening effect on the fracture behaviors of frozen rock,and confirms again that -4°C is a critical temperature for the softening effect.

    According to the above analyses,pointAindicates initiation of the microcrack propagation,pointBdenotes the beginning of microcrack coalescence,and pointCsignifies the occurrence of macrocrack propagation.Therefore,loads at pointA,pointBand pointCcan be marked as the microcrack propagation load,microcrack coalescence load and macrocrack propagation load,respectively.The three loads are extracted and their variation trend with temperature is illustrated in Fig.7.It demonstrates that the thawing leads to a gradual decrease in microcrack propagation load,microcrack coalescence load,and macrocrack propagation load,and they all change in two stages.At-20°C to-4°C,the load at pointAis decreased to 0 kN,and the loads at pointBand pointCare moderately decreased.From-4°C to 0°C,the loads at pointsBandCfall rapidly,while the load at pointAremains constant at 0 kN.

    Fig.7.The variations in microcrack propagation load,microcrack coalescence load,and macrocrack propagation load at different thawing temperatures.

    As shown in Fig.6(pointC),the range of FPZ varies at different temperatures.The FPZ at the peak load depicts the result of microcrack propagation and coalescence prior to macro-fracture.Therefore,many scholars determine the range of FPZ by the variation of the HSF,or displacement field (Wu et al.,2011;Lin et al.,2019).The method of determining the FPZ length by displacement field is shown in Fig.8a and b.

    Fig.8.(a) Incremental horizontal displacement contours at peak load;(b)opening displacements along the notch of the sample at -8 °C;and (c) FPZ length at different thawing temperatures.

    Fig.8a shows the incremental horizontal displacement contours on the sample surface as computed by DIC.To reduce the calculation error,horizontal displacement values on two reference lines along the notch length 2 mm to the side of the notch are extracted(the lines of L1and L2in Fig.8a),and shown in Fig.8b.The distance between the notch tip and the confluence point of two horizontal displacement lines delineates the length of FPZ:Lp(Lin and Labuz,2013).Fig.8c depicts the variation of FPZ length at the peak load with temperatures.Due to thawing,the FPZ length increases and demonstrates a two-stage mode of (a) moderate increase and (b)rapid increase.Unsurprisingly,the increase in FPZ length begins to accelerate at -4°C.

    3.3.The influence of thawing on AE signals during fracturing of frozen rock

    3.3.1.Changes in AE counts and cumulative AE counts

    The changes in AE signals during fracturing of frozen rock under various testing temperatures are shown in Fig.9.At the beginning of loading,the AE counts are weak;with expansion of the microcrack,the AE counts rise.The cumulative AE counts grow in stages during microcrack propagation and exhibit a strong correlation with the change of HS.However,variation of the AE counts is sensitive to the temperature.

    At-20°C ≥T>-4°C,the AE counts gradually intensify during microcracks expansion.The cumulative AE counts rise in three stages (Fig.9a-f).In the elastic deformation stage,microcracks basically do not expand,and the cumulative AE counts rise slowly.In the microcrack propagation stage,increase in the cumulative AE counts is greater;while in the microcrack coalescence stage,the cumulative AE counts rise rapidly.However,as the testing temperature rises over-4°C,the AE counts are much lower than that at-20°C to-4°C(Fig.9g-h).Moreover,the cumulative AE counts grow in two stages: a slow increase during the microcrack propagation stage and a rapid increase during the microcrack coalescence stage.When the ice in the sample is completely melted (0°C),the AE counts during microcrack expansion are nearly zero,and the intensity of the AE signal is extremely low(Fig.9i).It demonstrates that at-4°C-0°C,the expansion of a large number of microcracks releases less elastic energy.

    3.3.2.Change of AE b value

    To quantify the magnitude of seismicity,the G-R formula was used to express the statistical correlation between earthquake magnitude and frequency (Gutenberg and Richter,1944):

    whereMis the magnitude,Nis the number of events for which the magnitude is higher thanM,andaandbare the constants.

    The parametersbandahave specific physical significance.Thebvalue represents the ratio of small-magnitude events to largemagnitude ones.The destruction of a rock sample is comparable to a natural earthquake,with the exception of differences in magnitude and radiation frequency.Consequently,the relationship between the amplitude of an earthquake and its frequency can be utilized to analyze AE events during rock failure.Eq.(3)is expressed as follows (Colombo et al.,2003;Goebel et al.,2013):

    whereAdBis the AE event amplitude,andN′is the hit number of AE events whose amplitude exceedsAdB.Thebvalue is used to analyze the variation in AE events and rock failure.This demonstrates that during rock failure,when microcrack propagation is dominant,higherb-values can be noticed;whereas macrocrack propagation occurs,bvalues fall (Colombo et al.,2003;Nejati and Ghazvinian,2013).

    Fig.10 shows variation of the AEbvalue during thawing.The AEbvalue exhibits two stages with temperatures:basically unchanged and rapidly increasing.At-20°C to-4°C,the AEbvalues remain essentially unchanged,indicating that thawing has slight impact on the amplitude and number of AE events during the loading process.At-4°C-0°C,AEbvalues increase rapidly,indicating that the large amplitude AE events decrease while the small amplitude events rise.In other words,the microcrack propagation is more significant at-4°C-0°C.This is consistent with the results displayed in Fig.9.

    Fig.10.AE b-value at different thawing temperatures.

    3.4.Variation of pore ice during thawing of frozen rock

    The fracture behaviors of frozen rock are susceptible to temperature change,which should be attributed to thawing of the pore ice.Based on the principles in Section 2.3,the variation of pore ice content during thawing can be determined by subtracting the unfrozen water content from overall water content.

    The variation of ice and unfrozen water content is shown in Fig.11.It shows that at-20°C to-4°C,the unfrozen water content gradually rises and the ice content decreases,indicating that a small amount of ice is melted.At-4°C-0°C,the unfrozen water content rapidly grows,and the ice content quickly decreases,showing that a lot of ice is melted.At 0°C,the ice content is zero,which suggests that the ice in the pores of frozen rock is completely melted.Coincidentally,the critical temperature for the variation of ice content is also -4°C,which may hold the interpretation of thawing-induced softening effects on the fracture behaviors of frozen rock.

    Fig.11.The relative unfrozen water content and ice content at varying thawing temperatures.

    4.Discussion

    Thawing can induce a significant softening effect on the fracture behaviors of frozen rock,which is evidenced by the decrease inKIC,fracture energy (Fig.5),and the increases in horizontal strain(Fig.6)and AE signals(Fig.9)with rising temperature from-20°C to 0°C.Moreover,-4°C is recognized as a critical temperature for the softening effect,where the softening effect intensifies (see Figs.5-8,10).Coincidentally,the NMR results demonstrate a reverse variation trend in pore ice content with increase of temperature,that is,a moderate decrease is followed by a sharp decrease and-4°C is exactly the critical point(Fig.11).Accordingly,we interpreted the thawing-induced softening effect by linking evolution in the microscopic structure of frozen rock with its macroscopic fracture behaviors.

    When water-saturated rocks were frozen at -20°C,a vast majority of pore water was froze but the rest remained unfrozen in the form of a thin film between ice and mineral grains(Jia et al.,2019;Wang et al.,2022).Pore ice and unfrozen water film has two strengthening effects on the mineral skeleton (Wang et al.,2021):

    (1) The filling effect of pore ice.A stress concentration effect will occur at a crack tip when the crack is under tension,leading to an increased stress.The stress concentration effect in the crack is positively correlated with its length-diameter ratio(Hertzberg,1996).The presence of ice drastically reduces the effective length of the crack (Fig.12,at -20°C),hence relieves the stress concentration effect at the crack tip.The filling effect depends on the pore ice content.

    Fig.12.Influence mechanisms of thawing on fracture of frozen rock.

    (2) The cementation effect of the unfrozen water film.The ice is cemented to the rock skeleton through the unfrozen water film (Fig.12,at -20°C),hence increasing the resistance to tension failure of the crack.The cementation effect is negatively related to the thickness of unfrozen water film,which depends on temperature (Guerin et al.,2016;Chen et al.,2018;Wang et al.,2021).The melting of ice within the cracks diminishes the two strengthening effects,resulting in a degradation of the fracture mechanical properties of frozen rock.

    The melting of crack ice initiated from its outer boundary,i.e.the unfrozen water film.When the temperature increases from-20°C to-4°C,a small quantity of ice is melted,resulting in thickening of the unfrozen water film (Fig.12,at -4°C).So,the cementation strength between ice and rock skeleton decreases,making microcracks easier to propagate (Figs.6a-g and 7).This process also is responsible for the increase of FPZ length(see Fig.8),and the crack is more likely to fracture.However,as the above two strengthening effects are still significant,high-amplitude AE events dominate during the microcrack propagation process,and the AE signals are relatively strong(Figs.9a-g and 10).The weakening of the fracture mechanical parameters of frozen rocks during thawing is primarily determined by the thickness of the unfrozen water film within this temperature interval.

    As temperature is greater than-4°C,much of the ice is melted(Fig.12,at 0°C),and the unfrozen water content rises rapidly(Fig.11b).As a result,the cementation effect of ice almost vanishes and the filling effect of pore ice is reduced significantly,facilitating the propagation of microcracks(Figs.6g-i and 7).In addition,due to the weakening of the two strengthening effects,small amplitude AE events dominate the microcrack expansion process,and the amplitude of AE signals decreases significantly (Figs.9g-i and 10).After the specimen is completely melted,both strengthening effects disappear,and the AE signals are very small (Fig.9i).

    5.Conclusions

    A series of SCB tests were conducted on frozen rocks at various thawing temperatures,and the microcrack propagation process was monitored by AE and DIC methods.By linking the evolution in microscopic structure of frozen rock with its macroscopic fracture behaviors,the influence mechanisms of thawing were interpreted.

    (1) As the thawing temperature rises,the variation inKICof frozen rock exhibits a two-stage mode of slow decrease stage,KICdeclined by 19.5%from-20°C to-4°C;and rapid decrease,KICdeclined by 74%from-4°C to 0°C.The fracture energy varies in two stages: (a) essentially unchanged(-20°C~-12°C);and(b)rapidly decreased(-12°C-0°C).

    (2) During the loading process,when the thawing temperature is at -20°C to-4°C,the HS at the notch tip varies in three stages:(I)changeless stage;(II)moderate increase stage;and(III) rapid increase stage,which correspond to the three stages: elastic deformation stage,microcrack propagation stage,and the microcrack coalescence stage.However,the change of HS at the notch tip shows a two-stage mode of(II)moderate increase stage and (III) rapid increase stage at -4°C-0°C,which corresponds to the microcrack propagation stage and microcrack coalescence stage.

    (3) The microcrack propagation load,microcrack coalescence load,and macrocrack propagation load all decrease in two stages.At-20°C to-4°C,the microcrack propagation load is reduced to 0 kN,while the microcrack coalescence load and macrocrack propagation load are moderately decreased.At-4°C-0°C,the microcrack coalescence load and macrocrack propagation load fall rapidly,while the microcrack propagation load remains constant at 0 kN.Moreover,when the temperature increases,the FPZ length grows significantly.

    (4) During fracturing of frozen rock,the AE counts change in three stages at -20°C to -4°C and in two stages at 0°C to-4°C,which is consistent with the variation of HS at the notch tip.The AE counts at 0°C to-4°C are much fewer than that at -20°C to -4°C.In addition,at -4°C-0°C,the microcrack propagation at the notch tip is more prominent than that at -20°C to-4°C.

    (5) The weakening mechanism of the fracture mechanical parameters of frozen rocks differs at various thawing stages.At-20°C to-4°C,the thickening of the unfrozen water film diminishes the cementation strength between ice and rock skeleton,so the weakening of fracture mechanical parameters is primarily determined by the thickness of the unfrozen water film.However,at-4°C-0°C,the cementation effect of ice almost vanishes,and the filling effect of pore ice is reduced significantly,which facilitates microcrack propagation.

    Declaration of competing interest

    The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

    Acknowledgments

    We acknowledge the funding support from the National Natural Science Foundation of China (Grant No.42271148).

    一级毛片女人18水好多| 中文字幕人妻丝袜一区二区| 亚洲欧美一区二区三区黑人| 麻豆国产av国片精品| svipshipincom国产片| 无限看片的www在线观看| 少妇粗大呻吟视频| 18禁国产床啪视频网站| 日韩人妻精品一区2区三区| 首页视频小说图片口味搜索| 王馨瑶露胸无遮挡在线观看| 亚洲精品在线观看二区| 久久国产精品大桥未久av| 一级作爱视频免费观看| 亚洲精品av麻豆狂野| 动漫黄色视频在线观看| 亚洲成av片中文字幕在线观看| 久久久久久久精品吃奶| 久久久久国产一级毛片高清牌| 亚洲精品av麻豆狂野| 久久久久久免费高清国产稀缺| 国产av一区二区精品久久| 90打野战视频偷拍视频| 美女 人体艺术 gogo| 国产aⅴ精品一区二区三区波| 成年版毛片免费区| 亚洲视频免费观看视频| 波多野结衣一区麻豆| 国产日韩一区二区三区精品不卡| 变态另类成人亚洲欧美熟女 | 丝袜美足系列| 国产亚洲精品第一综合不卡| а√天堂www在线а√下载 | 91麻豆精品激情在线观看国产 | 午夜久久久在线观看| 老司机亚洲免费影院| av中文乱码字幕在线| 国产精品一区二区在线观看99| 国产欧美日韩综合在线一区二区| 王馨瑶露胸无遮挡在线观看| tube8黄色片| 国产精品自产拍在线观看55亚洲 | av欧美777| 大陆偷拍与自拍| 日韩视频一区二区在线观看| 日韩欧美免费精品| 午夜精品在线福利| 色尼玛亚洲综合影院| 久久 成人 亚洲| 看免费av毛片| 精品久久久精品久久久| 飞空精品影院首页| 久久天堂一区二区三区四区| 人妻 亚洲 视频| 精品视频人人做人人爽| 日韩人妻精品一区2区三区| www.精华液| 国产乱人伦免费视频| 国产精品香港三级国产av潘金莲| 母亲3免费完整高清在线观看| 大片电影免费在线观看免费| 国产视频一区二区在线看| 午夜福利一区二区在线看| 日韩人妻精品一区2区三区| 色尼玛亚洲综合影院| 久热爱精品视频在线9| 两个人免费观看高清视频| 精品福利永久在线观看| 国产亚洲欧美精品永久| 久久香蕉精品热| 国产不卡av网站在线观看| 又大又爽又粗| 久久 成人 亚洲| 亚洲av成人不卡在线观看播放网| 露出奶头的视频| 免费不卡黄色视频| 9热在线视频观看99| 午夜福利欧美成人| 国产成人精品久久二区二区91| 亚洲综合色网址| 色老头精品视频在线观看| av欧美777| 1024视频免费在线观看| 脱女人内裤的视频| tube8黄色片| 亚洲av熟女| 国产精品秋霞免费鲁丝片| 亚洲中文日韩欧美视频| 搡老岳熟女国产| 亚洲黑人精品在线| 亚洲欧洲精品一区二区精品久久久| 欧美成人午夜精品| 18禁国产床啪视频网站| 精品无人区乱码1区二区| 午夜91福利影院| 午夜影院日韩av| 精品免费久久久久久久清纯 | 中文字幕人妻丝袜一区二区| 亚洲人成伊人成综合网2020| 国产亚洲精品一区二区www | 色综合婷婷激情| 免费看a级黄色片| 91大片在线观看| 欧美乱色亚洲激情| 午夜福利一区二区在线看| 久久天躁狠狠躁夜夜2o2o| 黄色丝袜av网址大全| 国产亚洲一区二区精品| 黑丝袜美女国产一区| 亚洲一区高清亚洲精品| 国产精品亚洲一级av第二区| 美女高潮到喷水免费观看| 国产成人一区二区三区免费视频网站| 手机成人av网站| 国产免费av片在线观看野外av| 欧美精品一区二区免费开放| 天天躁日日躁夜夜躁夜夜| 亚洲国产毛片av蜜桃av| 黄片大片在线免费观看| 9191精品国产免费久久| 国产免费av片在线观看野外av| 两性夫妻黄色片| av国产精品久久久久影院| 自线自在国产av| 中文字幕另类日韩欧美亚洲嫩草| 久久人妻av系列| 飞空精品影院首页| 最近最新中文字幕大全免费视频| 两个人看的免费小视频| 亚洲专区中文字幕在线| 久久亚洲真实| 精品少妇久久久久久888优播| 欧美 日韩 精品 国产| 天天添夜夜摸| 美女扒开内裤让男人捅视频| 男人舔女人的私密视频| 大型av网站在线播放| 国产aⅴ精品一区二区三区波| 少妇粗大呻吟视频| 国产精品久久久久久人妻精品电影| 国内久久婷婷六月综合欲色啪| 日韩中文字幕欧美一区二区| 天堂动漫精品| 女人久久www免费人成看片| 搡老熟女国产l中国老女人| 中文字幕精品免费在线观看视频| 亚洲欧美精品综合一区二区三区| 免费看a级黄色片| 国产一区二区三区综合在线观看| 久久久久国内视频| 美女扒开内裤让男人捅视频| 国产在视频线精品| 精品电影一区二区在线| 日韩欧美一区视频在线观看| 国产一区二区三区在线臀色熟女 | 午夜福利在线观看吧| 亚洲精品国产一区二区精华液| 欧美性长视频在线观看| 黑人巨大精品欧美一区二区蜜桃| videosex国产| 夫妻午夜视频| 纯流量卡能插随身wifi吗| 国产有黄有色有爽视频| 人人妻人人添人人爽欧美一区卜| 人人妻,人人澡人人爽秒播| 亚洲av熟女| 淫妇啪啪啪对白视频| 国产精品影院久久| 69av精品久久久久久| 精品第一国产精品| 国产日韩欧美亚洲二区| 黄片大片在线免费观看| 一本综合久久免费| 国产一卡二卡三卡精品| 久久精品国产99精品国产亚洲性色 | 电影成人av| 欧美最黄视频在线播放免费 | 黄片小视频在线播放| 久久久精品国产亚洲av高清涩受| 在线观看66精品国产| 亚洲精品国产精品久久久不卡| 久久草成人影院| 亚洲三区欧美一区| 51午夜福利影视在线观看| 91大片在线观看| 纯流量卡能插随身wifi吗| 精品乱码久久久久久99久播| 精品久久久精品久久久| 水蜜桃什么品种好| 黄色女人牲交| 欧美在线一区亚洲| 欧美成人免费av一区二区三区 | 久热爱精品视频在线9| av中文乱码字幕在线| 日本黄色日本黄色录像| 超色免费av| aaaaa片日本免费| 在线播放国产精品三级| 精品国产国语对白av| 国产xxxxx性猛交| 男人操女人黄网站| 国产精品一区二区精品视频观看| 黄色视频,在线免费观看| 法律面前人人平等表现在哪些方面| 免费在线观看影片大全网站| 麻豆av在线久日| 色综合婷婷激情| 国产蜜桃级精品一区二区三区 | 亚洲美女黄片视频| 一区在线观看完整版| 热99国产精品久久久久久7| 日本精品一区二区三区蜜桃| 午夜成年电影在线免费观看| 国产日韩欧美亚洲二区| 在线观看日韩欧美| 亚洲熟女毛片儿| 丝袜美腿诱惑在线| 免费一级毛片在线播放高清视频 | 极品少妇高潮喷水抽搐| 欧美人与性动交α欧美精品济南到| 欧美国产精品va在线观看不卡| 深夜精品福利| 成人18禁在线播放| 搡老岳熟女国产| 后天国语完整版免费观看| 国产精品久久久久成人av| 欧美黑人欧美精品刺激| 高清在线国产一区| 一级a爱片免费观看的视频| 这个男人来自地球电影免费观看| 成年人午夜在线观看视频| 亚洲精品粉嫩美女一区| 久久热在线av| 国产蜜桃级精品一区二区三区 | 91成人精品电影| 亚洲精品美女久久av网站| 999久久久国产精品视频| 国产乱人伦免费视频| 免费少妇av软件| 欧美激情高清一区二区三区| 亚洲成a人片在线一区二区| 精品国产国语对白av| 国产在视频线精品| 欧美亚洲日本最大视频资源| 啦啦啦 在线观看视频| 精品一区二区三区四区五区乱码| 国产免费男女视频| 国产高清激情床上av| 18在线观看网站| 日日夜夜操网爽| 婷婷精品国产亚洲av在线 | 啦啦啦视频在线资源免费观看| 日韩欧美在线二视频 | 免费看a级黄色片| 99国产极品粉嫩在线观看| 亚洲少妇的诱惑av| 久久精品人人爽人人爽视色| 叶爱在线成人免费视频播放| 国产色视频综合| 99精品在免费线老司机午夜| 1024香蕉在线观看| √禁漫天堂资源中文www| 国产深夜福利视频在线观看| 欧美人与性动交α欧美精品济南到| 国产精品1区2区在线观看. | 色播在线永久视频| 亚洲,欧美精品.| 久久久久国产一级毛片高清牌| 亚洲专区中文字幕在线| 高清毛片免费观看视频网站 | a在线观看视频网站| 国产主播在线观看一区二区| 中文字幕另类日韩欧美亚洲嫩草| 国产成人一区二区三区免费视频网站| 午夜日韩欧美国产| 无遮挡黄片免费观看| 中文字幕人妻丝袜制服| 日韩欧美在线二视频 | 免费观看人在逋| 国产精品永久免费网站| 超碰97精品在线观看| 国产精品.久久久| 首页视频小说图片口味搜索| 久久久久国产精品人妻aⅴ院 | 高清在线国产一区| 老熟妇仑乱视频hdxx| 国产av精品麻豆| 亚洲国产中文字幕在线视频| 日本黄色视频三级网站网址 | 9191精品国产免费久久| 国产淫语在线视频| 香蕉国产在线看| 老司机深夜福利视频在线观看| 午夜福利一区二区在线看| 精品无人区乱码1区二区| 欧美精品亚洲一区二区| 日韩熟女老妇一区二区性免费视频| 免费久久久久久久精品成人欧美视频| 国产欧美日韩一区二区三区在线| 99国产精品一区二区蜜桃av | av有码第一页| 精品午夜福利视频在线观看一区| 成人特级黄色片久久久久久久| 日韩欧美免费精品| 久久久久久久久免费视频了| 999久久久精品免费观看国产| 精品国产一区二区三区久久久樱花| 国产成人欧美在线观看 | 久久久久久久午夜电影 | 一级片'在线观看视频| 在线观看免费视频日本深夜| 国产精品九九99| 中出人妻视频一区二区| 亚洲精品成人av观看孕妇| 亚洲精品美女久久av网站| 国产97色在线日韩免费| 侵犯人妻中文字幕一二三四区| 老汉色av国产亚洲站长工具| 国产成人欧美在线观看 | 欧美亚洲 丝袜 人妻 在线| 国产99久久九九免费精品| 免费看a级黄色片| 亚洲五月婷婷丁香| 日韩免费高清中文字幕av| 日韩欧美三级三区| 免费在线观看亚洲国产| 下体分泌物呈黄色| 亚洲国产精品一区二区三区在线| 色婷婷久久久亚洲欧美| 久久精品国产亚洲av高清一级| 美国免费a级毛片| 欧美激情极品国产一区二区三区| 国精品久久久久久国模美| 亚洲熟妇中文字幕五十中出 | 国产激情欧美一区二区| 亚洲专区国产一区二区| 亚洲中文日韩欧美视频| 欧美av亚洲av综合av国产av| 国产极品粉嫩免费观看在线| 精品少妇一区二区三区视频日本电影| 国产伦人伦偷精品视频| 亚洲在线自拍视频| 波多野结衣av一区二区av| 国产精品自产拍在线观看55亚洲 | 天天躁日日躁夜夜躁夜夜| 波多野结衣av一区二区av| 欧美日韩成人在线一区二区| 久久久久久久精品吃奶| 久久性视频一级片| 国产精品永久免费网站| 亚洲va日本ⅴa欧美va伊人久久| 免费看a级黄色片| 18禁美女被吸乳视频| 男女之事视频高清在线观看| 欧美激情高清一区二区三区| 国产亚洲欧美精品永久| 免费在线观看影片大全网站| 亚洲欧美一区二区三区黑人| 亚洲在线自拍视频| 天天影视国产精品| 波多野结衣一区麻豆| 久久精品亚洲熟妇少妇任你| 久久久久精品人妻al黑| 国产精品亚洲av一区麻豆| av中文乱码字幕在线| 麻豆av在线久日| 精品国产国语对白av| 在线免费观看的www视频| 国产精品亚洲av一区麻豆| 麻豆乱淫一区二区| 91麻豆精品激情在线观看国产 | 国产成人免费无遮挡视频| www.999成人在线观看| 高潮久久久久久久久久久不卡| 午夜福利,免费看| 欧美另类亚洲清纯唯美| 黄色女人牲交| 亚洲情色 制服丝袜| 18禁国产床啪视频网站| 高潮久久久久久久久久久不卡| 欧美黑人欧美精品刺激| 啦啦啦视频在线资源免费观看| 国产精品.久久久| 亚洲中文日韩欧美视频| 欧美激情高清一区二区三区| a级毛片黄视频| 老鸭窝网址在线观看| 精品福利观看| 99riav亚洲国产免费| 一区二区三区激情视频| 久久香蕉激情| 99国产极品粉嫩在线观看| 最新在线观看一区二区三区| 一级a爱片免费观看的视频| 成熟少妇高潮喷水视频| 国产不卡av网站在线观看| 91大片在线观看| 1024香蕉在线观看| 久久青草综合色| 在线天堂中文资源库| 日日摸夜夜添夜夜添小说| 亚洲av日韩精品久久久久久密| 亚洲国产欧美网| 久久人妻福利社区极品人妻图片| 亚洲一区中文字幕在线| 丝袜人妻中文字幕| 每晚都被弄得嗷嗷叫到高潮| 很黄的视频免费| 日韩成人在线观看一区二区三区| 国产亚洲欧美精品永久| 亚洲视频免费观看视频| 免费在线观看日本一区| 色94色欧美一区二区| 亚洲一区二区三区不卡视频| 国产精品久久久久久人妻精品电影| 亚洲专区国产一区二区| 欧美色视频一区免费| 午夜精品国产一区二区电影| 成人国语在线视频| 午夜福利乱码中文字幕| 欧美日韩亚洲国产一区二区在线观看 | 精品午夜福利视频在线观看一区| 中文字幕最新亚洲高清| 日韩中文字幕欧美一区二区| 精品久久久久久久毛片微露脸| 夜夜躁狠狠躁天天躁| 精品国产一区二区三区四区第35| 很黄的视频免费| 黄色 视频免费看| 午夜老司机福利片| tocl精华| 女人爽到高潮嗷嗷叫在线视频| 成年人午夜在线观看视频| 国产精品综合久久久久久久免费 | 性色av乱码一区二区三区2| 丰满迷人的少妇在线观看| 国产淫语在线视频| 夜夜夜夜夜久久久久| 18禁观看日本| 久久性视频一级片| 国产精品98久久久久久宅男小说| 高清黄色对白视频在线免费看| 热re99久久国产66热| 日本黄色日本黄色录像| 最近最新中文字幕大全电影3 | 国产蜜桃级精品一区二区三区 | 美女午夜性视频免费| 99久久精品国产亚洲精品| 每晚都被弄得嗷嗷叫到高潮| 国产精品久久视频播放| 国产亚洲精品久久久久久毛片 | 十八禁高潮呻吟视频| 老司机午夜福利在线观看视频| 久热这里只有精品99| 亚洲在线自拍视频| 黑丝袜美女国产一区| 久久久久国产精品人妻aⅴ院 | 婷婷精品国产亚洲av在线 | 精品人妻熟女毛片av久久网站| 真人做人爱边吃奶动态| 亚洲国产精品一区二区三区在线| 天堂俺去俺来也www色官网| 久99久视频精品免费| 波多野结衣av一区二区av| 国产精品永久免费网站| 大香蕉久久网| 欧美精品啪啪一区二区三区| 日韩大码丰满熟妇| 嫩草影视91久久| 人妻久久中文字幕网| 国产不卡av网站在线观看| 啦啦啦在线免费观看视频4| 少妇裸体淫交视频免费看高清 | 久久精品国产a三级三级三级| 亚洲综合色网址| 黄片小视频在线播放| 一区二区三区国产精品乱码| 69av精品久久久久久| 精品一区二区三区视频在线观看免费 | 老鸭窝网址在线观看| 亚洲精品中文字幕在线视频| 亚洲一区高清亚洲精品| 一进一出好大好爽视频| 99国产精品99久久久久| 人成视频在线观看免费观看| 老司机靠b影院| 女人精品久久久久毛片| 国产不卡av网站在线观看| 久久国产亚洲av麻豆专区| 亚洲七黄色美女视频| 成人特级黄色片久久久久久久| 精品国产一区二区久久| 久热这里只有精品99| 国产激情欧美一区二区| 久久99一区二区三区| 热re99久久精品国产66热6| 精品久久久久久电影网| 少妇被粗大的猛进出69影院| 真人做人爱边吃奶动态| 夜夜爽天天搞| 亚洲专区中文字幕在线| 国产精品久久电影中文字幕 | 中文字幕最新亚洲高清| 免费高清在线观看日韩| 又紧又爽又黄一区二区| 很黄的视频免费| 久久久精品免费免费高清| 免费观看精品视频网站| 满18在线观看网站| 飞空精品影院首页| 成熟少妇高潮喷水视频| 啦啦啦 在线观看视频| 欧美国产精品va在线观看不卡| 久久精品亚洲av国产电影网| 久久久国产欧美日韩av| 久久人人爽av亚洲精品天堂| 中文字幕人妻丝袜制服| 妹子高潮喷水视频| 成年女人毛片免费观看观看9 | 色尼玛亚洲综合影院| 国产精品.久久久| 欧美日韩精品网址| 国产国语露脸激情在线看| 一区福利在线观看| 日日爽夜夜爽网站| 天天操日日干夜夜撸| 欧美日韩亚洲高清精品| 91麻豆av在线| 亚洲熟妇熟女久久| 激情在线观看视频在线高清 | 国产又爽黄色视频| 久久青草综合色| 亚洲国产精品一区二区三区在线| 每晚都被弄得嗷嗷叫到高潮| 欧美日韩亚洲国产一区二区在线观看 | 欧美日韩视频精品一区| 亚洲成人免费av在线播放| 中文字幕最新亚洲高清| 欧美日韩亚洲综合一区二区三区_| 在线观看午夜福利视频| 国产精品影院久久| 中文字幕最新亚洲高清| 叶爱在线成人免费视频播放| 在线观看www视频免费| 国产97色在线日韩免费| 深夜精品福利| 女性被躁到高潮视频| 久久香蕉激情| 国产高清激情床上av| 热re99久久国产66热| 久久久久国内视频| 午夜影院日韩av| 人人澡人人妻人| 日本a在线网址| 亚洲精品中文字幕在线视频| 老司机福利观看| 久久久水蜜桃国产精品网| 欧美不卡视频在线免费观看 | 精品久久久精品久久久| 中文欧美无线码| 欧美av亚洲av综合av国产av| 免费在线观看日本一区| 夜夜爽天天搞| 日本五十路高清| 两个人看的免费小视频| 日本a在线网址| 中文字幕制服av| 黑人猛操日本美女一级片| 天天添夜夜摸| 亚洲精品在线美女| av网站在线播放免费| 国产精品.久久久| 亚洲精品成人av观看孕妇| 欧美丝袜亚洲另类 | 他把我摸到了高潮在线观看| 王馨瑶露胸无遮挡在线观看| 桃红色精品国产亚洲av| 国产又爽黄色视频| 一级毛片高清免费大全| 超色免费av| 国产xxxxx性猛交| 91精品国产国语对白视频| 一边摸一边抽搐一进一小说 | 久久国产精品人妻蜜桃| 亚洲欧美一区二区三区黑人| 18禁观看日本| 久久午夜综合久久蜜桃| 亚洲七黄色美女视频| 91字幕亚洲| 国产高清国产精品国产三级| 日本a在线网址| 伊人久久大香线蕉亚洲五| av天堂久久9| 亚洲七黄色美女视频| www.自偷自拍.com| xxx96com| 免费人成视频x8x8入口观看| 黑人巨大精品欧美一区二区蜜桃| 两个人看的免费小视频| 天天操日日干夜夜撸| 午夜免费成人在线视频| 国产成人av激情在线播放| 母亲3免费完整高清在线观看| 涩涩av久久男人的天堂| 亚洲色图综合在线观看| 久久精品91无色码中文字幕| 在线视频色国产色| 99精品在免费线老司机午夜| 18禁美女被吸乳视频| 午夜两性在线视频| 精品熟女少妇八av免费久了| 飞空精品影院首页| 国产精品国产av在线观看| 法律面前人人平等表现在哪些方面| 欧美老熟妇乱子伦牲交|